Whole Transcriptomic Analysis of Key Genes and Signaling Pathways in Endogenous ARDS

OBJECTIVE: To analyze the differentially expressed genes (DEGs) in rats with endogenous acute respiratory distress syndrome (ARDS) lung injury and explore the pathogenesis and early diagnostic molecular markers using whole transcriptomic data. METHODS: Twelve 8-week-old male Sprague Dawley rats were selected and randomly and equally divided into ARDS lung injury group and normal control group. RNA was extracted from the left lung tissues of both the groups and sequenced using the paired-end sequencing mode of the Illumina Hiseq sequencing platform. The DEGs of miRNA, cirRNA, lncRNA, and mRNA were screened using DESeq2 software, and the ceRNA regulatory network was constructed using Cytoscape. Gene ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) enrichment analysis were performed using the mRNA DEGs. STRING and Cytoscape software were used to construct the protein interaction network and identify the 15 key genes, which were verified using quantitative real-time polymerase chain reaction (qRT-PCR). RESULTS: Based on different screening conditions, and compared with the control group, the ARDS lung injury group showed 836 mRNA DEGs (386 upregulated and 450 downregulated), 110 lncRNA DEGs (53 upregulated and 57 downregulated), 19 circRNA DEGs (3 upregulated and 16 downregulated), and 6 miRNA DEGs (5 upregulated and 1 downregulated gene). GO showed that the DEGs of mRNA were mainly involved in biological processes, such as defense response to lipopolysaccharide and other organisms, leukocyte chemotaxis, neutrophil chemotaxis, and cytokine-mediated signaling. KEGG enrichment analysis showed that the DEGs played their biological roles mainly by participating in IL-17, TNF, and chemokine signaling pathways. The PPI analysis showed a total of 281 node proteins and 634 interaction edges. The top 15 key genes, which were screened, included Cxcl10, Mx1, Irf7, Isg15, Ifit3, Ifit2, Rsad2, Ifi47, Oasl, Dhx58, Usp18, Cmpk2, Herc6, Ifit1, and Gbp4. The ceRNA network analysis showed 69 nodes and 73 correlation pairs, where the key gene nodes were miR-21-3p, Camk2g, and Stx2. CONCLUSIONS: The chemotaxis, migration, and degranulation of inflammatory cells, cytokine immune response, autophagy, and apoptosis have significant biological functions in the occurrence and development of endogenous acute lung injury during ARDS. Thus, the camk2g/miR-21-3p/lncRNA/circRNA network, CXCL10/CXCR3, and IL-17 signaling pathways might provide novel insights and targets for further studying the lung injury mechanism and clinical treatment.