The aim of the present study was to identify differentially expressed genes (DEGs) and their related functions and pathways of major burn injuries, and to prevent the occurrence of complications. The expression profiling of E-GEOD-37069 was downloaded from ArrayExpress Archive. The DEGs of major burn injuries were identified. Gene Ontology (GO) and Kyoto Encyclopedia of Gene and Genomes (KEGG) functional enrichment analysis were then performed for the DEGs. Based on the KEGG database, a pathway relationship network was constructed, and DEGs in significant GO terms and pathways were investigated. Gene signal network and gene co-expression network of these inserted DEGs were constructed. A total of 3,328 DEGs of major burn injuries were identified, including 1,337 up- and 1,991 downregulated DEGs. These DEGs were mainly enriched into various GO terms, including transcription, DNA-dependent, signal transduction and blood coagulation. Moreover, they were also enriched into different pathways, such as hematopoietic cell lineage, metabolic pathway and chemokine signaling pathway. The pathway relationship network was constructed with 72 nodes. The MAPK signaling pathway was the hub node. Based on the same gene symbol, 702 DEGs were obtained, identified in both GO terms and pathways. Finally, the gene signaling network and gene co-expression network were constructed with 391 and 128 nodes, respectively. These identified DEGs, including GNB2, LILRA2, ARRB2 and ARHGEF2, may be potential key genes involved in the treatment of major burn injuries and prevention of complications.

The aim of the study was to identify key long non-coding RNAs (lncRNA) and related subpathways following severe burn injuries and research their functions. The miRNA-mRNA and lncRNA-miRNA interactions were downloaded from starBase v2.0 database. In addition, mRNA-miRNA interactions were obtained from TarBase, mirTarBase, mir2Disease, miRecords (V4.0) databases. The relationships of lncRNA-miRNA-mRNA were constructed. Genes of expression profiling were intersected with mRNA and lncRNA in lncRNA-mRNA interaction. Screened mRNAs were enriched into various pathways and screened lncRNAs were embedded into candidate pathways. Wallenius approximation methods were used to calculate the false discovery rate value of each sub-pathway. Based on the results of significant sub-pathways, the related lncRNA-mRNA network was constructed. A total of 18,081 genes were obtained. The lncRNA-mRNA intersections including 835 lncRNAs, 1,749 mRNAs and 7,693 interacting pairs were constructed. The enriched mRNAs were further enriched into various candidate pathways such as ribosome biogenesis in eukaryotes. Several sub-pathways were screened, including ribosome biogenesis in eukaryotes and MAPK signaling pathway. The network of pathway-lncRNA-mRNA was constructed. Hub-genes were identified, including C14orf169 and YLPM1. Several hub-lncRNAs were obtained, including PRKAG2 antisense RNA 1 and LEF1 antisense RNA 1. Several hub-lncRNAs including C14orf169, YLPM1, TTTY15, and PCBP1-AS1 were screened. The sub-pathways regulated by these lncRNAs were identified, and functions were predicted.