Identification of key genes as potential diagnostic biomarkers in sepsis by bioinformatics analysis

Abstract

Background Sepsis, an infection-triggered inflammatory syndrome, poses a global clinical challenge with limited therapeutic options. Our study is designed to identify potential diagnostic biomarkers of sepsis onset in critically ill patients by bioinformatics analysis. Methods Gene expression profiles of GSE28750 and GSE74224 were obtained from the Gene Expression Omnibus (GEO) database. These datasets were merged, normalized and de-batched. Weighted gene co-expression network analysis (WGCNA) was performed and the gene modules most associated with sepsis were identified as key modules. Functional enrichment analysis of the key module genes was then conducted. Moreover, differentially expressed gene (DEG) analysis was conducted by the “limma” R package. Protein-protein interaction (PPI) network was created using STRING and Cytoscape, and PPI hub genes were identified with the cytoHubba plugin. The PPI hub genes overlapping with the genes in key modules of WGCNA were determined to be the sepsis-related key genes. Subsequently, the key overlapping genes were validated in an external independent dataset and sepsis patients recruited in our hospital. In addition, CIBERSORT analysis evaluated immune cell infiltration and its correlation with key genes. Results By WGCNA, the greenyellow module showed the highest positive correlation with sepsis (0.7, p = 2e − 19). 293 DEGs were identified in the merged datasets. The PPI network was created, and the CytoHubba was used to calculate the top 20 genes based on four algorithms (Degree, EPC, MCC, and MNC). Ultimately, LTF, LCN2, ELANE, MPO and CEACAM8 were identified as key overlapping genes as they appeared in the PPI hub genes and the key module genes of WGCNA. These sepsis-related key genes were validated in an independent external dataset (GSE131761) and sepsis patients recruited in our hospital. Additionally, the immune infiltration profiles differed significantly between sepsis and non-sepsis critical illness groups. Correlations between immune cells and these five key genes were assessed, revealing that plasma cells, macrophages M0, monocytes, T cells regulatory, eosinophils and NK cells resting were simultaneously and significantly associated with more than two key genes. Conclusion This study suggests a critical role of LTF, LCN2, ELANE, MPO and CEACAM8 in sepsis and may provide potential diagnostic biomarkers and therapeutic targets for the treatment of sepsis.