Six potential biomarkers in septic shock: a deep bioinformatics and prospective observational study

Abstract

Background: Septic shock occurs when sepsis is related to severe hypotension and leads to a remarkable high number of deaths. The early diagnosis of septic shock is essential to reduce mortality. High-quality biomarkers can be objectively measured and evaluated as indicators to accurately predict disease diagnosis. However, single-gene prediction efficiency is inadequate; therefore, we identified a risk-score model based on gene signature to elevate predictive efficiency.

Methods: The gene expression profiles of GSE33118 and GSE26440 were downloaded from the Gene Expression Omnibus (GEO) database. These two datasets were merged, and the differentially expressed genes (DEGs) were identified using the limma package in R software. Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway enrichments of DEGs were performed. Subsequently, Lasso regression and Boruta feature selection algorithm were combined to identify the hub genes of septic shock. GSE9692 was then subjected to weighted gene co-expression network analysis (WGCNA) to identify the septic shock-related gene modules. Subsequently, the genes within such modules that matched with septic shock-related DEGs were identified as the hub genes of septic shock. To further understand the function and signaling pathways of hub genes, we performed gene set variation analysis (GSVA) and then used the CIBERSORT tool to analyze the immune cell infiltration pattern of diseases. The diagnostic value of hub genes in septic shock was determined using receiver operating characteristic (ROC) analysis and verified using quantitative PCR (qPCR) and Western blotting in our hospital patients with septic shock.

Results: A total of 975 DEGs in the GSE33118 and GSE26440 databases were obtained, of which 30 DEGs were remarkably upregulated. With the use of Lasso regression and Boruta feature selection algorithm, six hub genes (CD177, CLEC5A, CYSTM1, MCEMP1, MMP8, and RGL4) with expression differences in septic shock were screened as potential diagnostic markers for septic shock among the significant DEGs and were further validated in the GSE9692 dataset. WGCNA was used to identify the co-expression modules and module-trait correlation. Enrichment analysis showed significant enrichment in the reactive oxygen species pathway, hypoxia, phosphatidylinositol 3-kinases (PI3K)/Protein Kinase B (AKT)/mammalian target of rapamycin (mTOR) signaling, nuclear factor-κβ/tumor necrosis factor alpha (NF-κβ/TNF-α), and interleukin-6 (IL-6)/Janus Kinase (JAK)/Signal Transducers and Activators of Transcription 3 (STAT3) signaling pathways. The receiver operating characteristic curve (ROC) of these signature genes was 0.938, 0.914, 0.939, 0.956, 0.932, and 0.914, respectively. In the immune cell infiltration analysis, the infiltration of M0 macrophages, activated mast cells, neutrophils, CD8 T cells, and naive B cells was more significant in the septic shock group. In addition, higher expression levels of CD177, CLEC5A, CYSTM1, MCEMP1, MMP8, and RGL4 messenger RNA (mRNA) were observed in peripheral blood mononuclear cells (PBMCs) isolated from septic shock patients than from healthy donors. Higher expression levels of CD177 and MMP8 proteins were also observed in the PBMCs isolated from septic shock patients than from control participants.

Conclusions: CD177, CLEC5A, CYSTM1, MCEMP1, MMP8, and RGL4 were identified as hub genes, which were of considerable value in the early diagnosis of septic shock patients. These preliminary findings are of great significance for studying immune cell infiltration in the pathogenesis of septic shock, which should be further validated in clinical studies and basic studies.

Keywords: WGCNA; bioinformatic analysis; biomarker; prognosis; septic shock.