Identification of potential biomarkers and pathways for sepsis using RNA sequencing technology and bioinformatic analysis

Abstract

Long non-coding RNAs (lncRNAs) has been proven by many to play a crucial part in the process of sepsis. To obtain a better understanding of sepsis, the molecular biomarkers associated with it, and its possible pathogenesis, we obtained data from RNA-sequencing analysis using serum from three sepsis patients and three healthy controls (HCs). Using edgeR (one of the Bioconductor software package), we identified 1118 differentially expressed mRNAs (DEmRNAs) and 1394 differentially expressed long noncoding RNAs (DElncRNAs) between sepsis patients and HCs. We identified the biological functions of these disordered genes using Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) signaling pathway analyses. The GO analysis showed that the homophilic cell adhesion via plasma membrane adhesion molecules was the most significantly enriched category. The KEGG signaling pathway analysis indicated that the differentially expressed genes (DEGs) were most significantly enriched in retrograde endocannabinoid signaling. Using STRING, a protein-protein interaction network was also created, and Cytohubba was used to determine the top 10 hub genes. To examine the relationship between the hub genes and sepsis, we examined three datasets relevant to sepsis that were found in the gene expression omnibus (GEO) database. PTEN and HIST2H2BE were recognized as hub gene in both GSE4607, GSE26378, and GSE9692 datasets. The receiver operating characteristic (ROC) curves indicate that PTEN and HIST2H2BE have good diagnostic value for sepsis. In conclusion, this two hub genes may be biomarkers for the early diagnosis of sepsis, our findings should deepen our understanding of the pathogenesis of sepsis.