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.

Single-cell RNA sequencing reveals cell–cell communication and potential biomarker in sepsis and septic shock patients

Co-expression and interaction network analysis reveals dysregulated neutrophil and T-cell activation as the core mechanism associated with septic shock

Bulk and Single-Cell Transcriptome Analysis Reveal Shared Key Genes and Patterns of Immune Dysregulation in Sepsis and Systemic Lupus Erythematosus

Integrative analysis of metabolomics and transcriptomics to uncover biomarkers in sepsis

Bulk RNA Sequencing With Integrated Single-Cell RNA Sequencing Identifies BCL2A1 as a Potential Diagnostic and Prognostic Biomarker for Sepsis

Comprehensive Transcriptome Profiling of Peripheral Blood Mononuclear Cells from Patients with Sepsis

Single-cell transcriptomics reveals the alteration of peripheral blood mononuclear cells driven by sepsis

Dynamic changes in human single-cell transcriptional signatures during fatal sepsis

Identifying Six Chromatin Remodeling-related Genes As Diagnostic Biomarkers in Sepsis Using Bioinformatic Analyses

Single-cell transcriptional gene signature analysis identifies IL-17 signaling pathway as the key pathway in sepsis

Identification of Potential Biomarkers and Immune Features of Sepsis Using Bioinformatics Analysis

Single-cell RNA sequencing and immune repertoire analysis revealed dynamic immune characteristics associated with peripheral blood during sepsis

An immune-cell signature of bacterial sepsis

Deciphering the immune-metabolic nexus in sepsis: a single-cell sequencing analysis of neutrophil heterogeneity and risk stratification

Analysis of bulk RNA-seq data from sepsis patients reveals sepsis-associated lncRNAs and targeted cell death-related genes contributing to immune microenvironment regulation

Integrated analysis reveals NLRC4 as a potential biomarker in sepsis pathogenesis

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

Elucidating common pathogenic transcriptional networks in infective endocarditis and sepsis: integrated insights from biomarker discovery and single-cell RNA sequencing

Screening of prognostic core genes based on cell–cell interaction in the peripheral blood of patients with sepsis

Summary data-based Mendelian randomization and single-cell RNA sequencing analyses identify immune associations with low-level LGALS9 in sepsis

Single-cell RNA sequencing reveals the sustained immune cell dysfunction in the pathogenesis of sepsis secondary to bacterial pneumonia