Identification of four mitochondria-related genes in sepsis based on RNA sequencing technology

OBJECTIVES

The purpose of this study was to identify and analyze the mitochondrial genes associated with sepsis patients in order to elucidate the underlying mechanism of sepsis immunity and provide new ideas for the clinical treatment of sepsis.

METHODS

The hospitalized cases of sepsis (n = 20) and systemic inflammatory response syndrome (SIRS) (n = 12) admitted to the Emergency Intensive Care Unit (EICU) of the Affiliated Hospital of Southwest Medical University from January 2019 to December 2019 were collected consecutively. RNA-seq was used to sequence the RNA (mRNA) of peripheral blood cells. Bioinformatics techniques were used to screen and identify differentially expressed RNAs, with an absolute value of fold change (FC) greater than or equal to 1.2 and a false discovery rate (FDR) less than 0.05. At the same time, mitochondrial genes were obtained from the MitoCarta 3.0 database. Differential genes were then intersected with mitochondrial genes. The resulting crossover genes were subjected to GO, KEGG, and PPI analysis. Subsequently, the GSE65682 dataset was downloaded from the GEO database for survival analysis to assess the prognostic value of core genes, and GSE67652 was downloaded for ROC curve analysis to validate the diagnostic value of core genes. Finally, the localization of core genes was clarified through 10X single-cell sequencing.

RESULTS

The crossing of 314 sepsis differential genes and 1136 mitochondrial genes yielded 28 genes. GO and KEGG analysis showed that the crossover genes were mainly involved in the mitochondrion, mitochondrial matrix, and mitochondrial inner membrane. Survival analysis screened four genes that were significantly negatively associated with the prognosis of sepsis, namely FIS1, FKBP8, GLRX5, and GUK1. A comparison of peripheral blood RNA-seq results between the sepsis group and the SIRS group showed that the expression levels of these four genes were significantly decreased in the sepsis group compared to the SIRS group. ROC curve analysis based on GSE67652 indicates these four genes' high sensitivity and specificity for sepsis detection. Additionally, single-cell RNA sequencing found that the core genes were mainly expressed in macrophages, T cells, and B cells.

CONCLUSIONS

Mitochondria-related genes (FIS1, FKBP8, GLRX5, GUK1) were underexpressed in the sepsis group, negatively correlated with survival, and mainly distributed in immune cells. This finding may guide studying the immune-related mechanisms of sepsis. This study protocol was reviewed by the Ethics Committee of the Affiliated Hospital of Southwest Medical University (ethics number: KY2018029), the clinical trial registration number is ChiCTR1900021261, and the registration date is February 4, 2019.

Identification of sepsis-associated mitochondrial genes through RNA and single-cell sequencing approaches

BACKGROUND

Sepsis ranks among the most formidable clinical challenges, characterized by exorbitant treatment costs and substantial demands on healthcare resources. Mitochondrial dysfunction emerges as a pivotal risk factor in the pathogenesis of sepsis, underscoring the imperative to identify mitochondrial-related biomarkers. Such biomarkers are crucial for enhancing the accuracy of sepsis diagnostics and prognostication.

METHODS

In this study, adhering to the SEPSIS 3.0 criteria, we collected peripheral blood within 24 h of admission from 20 sepsis patients at the ICU of the Southwest Medical University Affiliated Hospital and 10 healthy volunteers as a control group for RNA-seq. The RNA-seq data were utilized to identify differentially expressed RNAs. Concurrently, mitochondrial-associated genes (MiAGs) were retrieved from the MitoCarta3.0 database. The differentially expressed genes were intersected with MiAGs. The intersected genes were then subjected to GO (Gene Ontology), and KEGG (Kyoto Encyclopedia of Genes and Genomes) analyses and core genes were filtered using the PPI (Protein-Protein Interaction) network. Subsequently, relevant sepsis datasets (GSE65682, GSE28750, GSE54514, GSE67652, GSE69528, GSE95233) were downloaded from the GEO (Gene Expression Omnibus) database to perform bioinformatic validation of these core genes. Survival analysis was conducted to assess the prognostic value of the core genes, while ROC (Receiver Operating Characteristic) curves determined their diagnostic value, and a meta-analysis confirmed the accuracy of the RNA-seq data. Finally, we collected 5 blood samples (2 normal controls (NC); 2 sepsis; 1 SIRS (Systemic Inflammatory Response Syndrome), and used single-cell sequencing to assess the expression levels of the core genes in the different blood cell types.

RESULTS

Integrating high-throughput sequencing with bioinformatics, this study identified two mitochondrial genes (COX7B, NDUFA4) closely linked with sepsis prognosis. Survival analysis demonstrated that patients with lower expression levels of COX7B and NDUFA4 exhibited a higher day survival rate over 28 days, inversely correlating with sepsis mortality. ROC curves highlighted the significant sensitivity and specificity of both genes, with AUC values of 0.985 for COX7B and 0.988 for NDUFA4, respectively. Meta-analysis indicated significant overexpression of COX7B and NDUFA4 in the sepsis group in contrast to the normal group (P < 0.01). Additionally, single-cell RNA sequencing revealed predominant expression of these core genes in monocytes-macrophages, T cells, and B cells.

CONCLUSION

The mitochondrial-associated genes (MiAGs) COX7B and NDUFA4 are intimately linked with the prognosis of sepsis, offering potential guidance for research into the mechanisms underlying sepsis.