Oxidative phosphorylation gene expression falls at onset and throughout the development of meningococcal sepsis-induced multi-organ failure in children

A TRANSCRIPTOMIC APPRECIATION OF CHILDHOOD MENINGOCOCCAL AND POLYMICROBIAL SEPSIS FROM A PRO-INFLAMMATORY AND TRAJECTORIAL PERSPECTIVE, A ROLE FOR VASCULAR ENDOTHELIAL GROWTH FACTOR A AND B MODULATION?

ABSTRACT

This study investigated the temporal dynamics of childhood sepsis by analyzing gene expression changes associated with proinflammatory processes. Five datasets, including four meningococcal sepsis shock (MSS) datasets (two temporal and two longitudinal) and one polymicrobial sepsis dataset, were selected to track temporal changes in gene expression. Hierarchical clustering revealed three temporal phases: early, intermediate, and late, providing a framework for understanding sepsis progression. Principal component analysis supported the identification of gene expression trajectories. Differential gene analysis highlighted consistent upregulation of vascular endothelial growth factor A (VEGF-A) and nuclear factor κB1 (NFKB1), genes involved in inflammation, across the sepsis datasets. NFKB1 gene expression also showed temporal changes in the MSS datasets. In the postmortem dataset comparing MSS cases to controls, VEGF-A was upregulated and VEGF-B downregulated. Renal tissue exhibited higher VEGF-A expression compared with other tissues. Similar VEGF-A upregulation and VEGF-B downregulation patterns were observed in the cross-sectional MSS datasets and the polymicrobial sepsis dataset. Hexagonal plots confirmed VEGF-R (VEGF receptor)-VEGF-R2 signaling pathway enrichment in the MSS cross-sectional studies. The polymicrobial sepsis dataset also showed enrichment of the VEGF pathway in septic shock day 3 and sepsis day 3 samples compared with controls. These findings provide unique insights into the dynamic nature of sepsis from a transcriptomic perspective and suggest potential implications for biomarker development. Future research should focus on larger-scale temporal transcriptomic studies with appropriate control groups and validate the identified gene combination as a potential biomarker panel for sepsis.

A community approach to mortality prediction in sepsis via gene expression analysis

Improved risk stratification and prognosis prediction in sepsis is a critical unmet need. Clinical severity scores and available assays such as blood lactate reflect global illness severity with suboptimal performance, and do not specifically reveal the underlying dysregulation of sepsis. Here, we present prognostic models for 30-day mortality generated independently by three scientific groups by using 12 discovery cohorts containing transcriptomic data collected from primarily community-onset sepsis patients. Predictive performance is validated in five cohorts of community-onset sepsis patients in which the models show summary AUROCs ranging from 0.765-0.89. Similar performance is observed in four cohorts of hospital-acquired sepsis. Combining the new gene-expression-based prognostic models with prior clinical severity scores leads to significant improvement in prediction of 30-day mortality as measured via AUROC and net reclassification improvement index These models provide an opportunity to develop molecular bedside tests that may improve risk stratification and mortality prediction in patients with sepsis.

Comprehensive Analysis of Gene Expression Profiles of Sepsis-Induced Multiorgan Failure Identified Its Valuable Biomarkers

Sepsis is an inflammatory-related disease, and severe sepsis would induce multiorgan dysfunction, which is the most common cause of death of patients in noncoronary intensive care units. Progression of novel therapeutic strategies has proven to be of little impact on the mortality of severe sepsis, and unfortunately, its mechanisms still remain poorly understood. In this study, we analyzed gene expression profiles of severe sepsis with failure of lung, kidney, and liver for the identification of potential biomarkers. We first downloaded the gene expression profiles from the Gene Expression Omnibus and performed preprocessing of raw microarray data sets and identification of differential expression genes (DEGs) through the R programming software; then, significantly enriched functions of DEGs in lung, kidney, and liver failure sepsis samples were obtained from the Database for Annotation, Visualization, and Integrated Discovery; finally, protein-protein interaction network was constructed for DEGs based on the STRING database, and network modules were also obtained through the MCODE cluster method. As a result, lung failure sepsis has the highest number of DEGs of 859, whereas the number of DEGs in kidney and liver failure sepsis samples is 178 and 175, respectively. In addition, 17 overlaps were obtained among the three lists of DEGs. Biological processes related to immune and inflammatory response were found to be significantly enriched in DEGs. Network and module analysis identified four gene clusters in which all or most of genes were upregulated. The expression changes of Icam1 and Socs3 were further validated through quantitative PCR analysis. This study should shed light on the development of sepsis and provide potential therapeutic targets for sepsis-induced multiorgan failure.

Transcriptomic data from two primary cell models stimulating human monocytes suggest inhibition of oxidative phosphorylation and mitochondrial function by N. meningitidis which is partially up-regulated by IL-10

Background

Biological interpretation of DNA microarray data may differ depending on underlying assumptions and statistical tests of bioinformatics tools used. We used Gene Set Enrichment Analysis (GSEA) and Ingenuity Pathway Analysis (IPA) to analyze previously generated DNA microarray data from human monocytes stimulated with N. meningitidis and IL-10 (“the model system”), and with meningococcal sepsis plasma before and after immunodepletion of IL-10 (“the patient plasma system”). The objectives were to compare if the two bioinformatics methods resulted in similar biological interpretation of the datasets, and to identify whether GSEA provided additional insight compared with IPA about the monocyte host response to meningococcal activation.

Results

In both experimental models, GSEA and IPA identified genes associated with pro-inflammatory innate immune activation, including TNF-signaling, Toll-like receptor signaling, JAK-STAT-signaling, and type I and type II interferon signaling. GSEA identified genes regulated by the presence of IL-10 with similar gene sets in both the model system and the patient plasma system. In the model system, GSEA and IPA in sum identified 170 genes associated with oxidative phosphorylation/mitochondrial function to be down-regulated in monocytes stimulated with meningococci. In the patient plasma system, GSEA and IPA in sum identified 122 genes associated with oxidative phosphorylation/mitochondrial dysfunction to be down-regulated by meningococcal sepsis plasma depleted for IL-10. Using IPA, we identified IL-10 to up-regulate 18 genes associated with oxidative phosphorylation/mitochondrial function that were down-regulated by N. meningitidis.

Conclusions

Biological processes associated with the gene expression changes in the model system of meningococcal sepsis were comparable with the results found in the patient plasma system. By combining GSEA with IPA, we discovered an inhibitory effect of N. meningitidis on genes associated with mitochondrial function and oxidative phosphorylation, and that IL-10 partially reverses this strong inhibitory effect, thereby identifying, to our knowledge, yet another group of genes where IL-10 regulates the effect of LPS. We suggest that relying on a single bioinformatics tool together with an arbitrarily chosen filtering criteria for data analysis may result in overlooking relevant biological processes and signaling pathways associated with genes differentially expressed between compared experimental conditions.

Electronic supplementary material

The online version of this article (10.1186/s12865-017-0229-5) contains supplementary material, which is available to authorized users.