Identification of ferroptosis-related genes in male mice with sepsis-induced acute lung injury based on transcriptome sequencing

Abstract

Abstract Background Sepsis can result in acute lung injury (ALI). Studies have shown that pharmacological inhibition of ferroptosis can treat ALI. However, the regulatory mechanisms of ferroptosis in sepsis-induced ALI remain unclear. Methods Transcriptome sequencing was performed on lung tissue samples from 10 sepsis-induced mouse models of ALI and 10 control mice. After quality control measures, clean data were used to screen for differentially expressed genes (DEGs) between the groups. The DEGs were then overlapped with ferroptosis-related genes (FRGs) to obtain ferroptosis-related DEGs (FR-DEGs). Subsequently, least absolute shrinkage and selection operator (Lasso) and Support Vector Machine-Recursive Feature Elimination (SVM-RFE) were used to obtain key genes. In addition, Ingenuity Pathway Analysis (IPA) was employed to explore the disease, function, and canonical pathways related to the key genes. Gene set enrichment analysis (GSEA) was used to investigate the functions of the key genes, and regulatory miRNAs of key genes were predicted using the NetworkAnalyst and StarBase databases. Finally, the expression of key genes was validated with the GSE165226 and GSE168796 datasets sourced from the Gene Expression Omnibus (GEO) database and using quantitative real-time polymerase chain reaction (qRT-PCR). Results Thirty-three FR-DEGs were identified between 1843 DEGs and 259 FRGs. Three key genes, Ncf2, Steap3, and Gclc, were identified based on diagnostic models established by the two machine learning methods. They are mainly involved in infection, immunity, and apoptosis, including lymphatic system cell migration and lymphocyte and T cell responses. Additionally, the GSEA suggested that Ncf2 and Steap3 were similarly enriched in mRNA processing, response to peptides, and leukocyte differentiation. Furthermore, a key gene-miRNA network including 2 key genes (Steap3 and Gclc) and 122 miRNAs, and a gene-miRNA network with 1 key gene (Steap3) and 3 miRNAs were constructed using NetworkAnalyst and StarBase, respectively. Both databases predicted that mmu-miR-15a-5p was the target miRNA of Steap3. Finally, Ncf2 expression was validated using both datasets and qRT-PCR, and Steap3 was validated using GSE165226 and qRT-PCR. Conclusions This study identified two FR-DEGs (Ncf2 and Steap3) associated with sepsis-induced ALI via transcriptome analyses, as well as their functional and metabolic pathways.