Prediction of sepsis mortality using metabolite biomarkers in the blood: a meta-analysis of death-related pathways and prospective validation
BMC Medicine volume 18, Article number: 83 (2020)
Sepsis is a leading cause of death in intensive care units (ICUs), but outcomes of individual patients are difficult to predict. The recently developed clinical metabolomics has been recognized as a promising tool in the clinical practice of critical illness. The objective of this study was to identify the unique metabolic biomarkers and their pathways in the blood of sepsis nonsurvivors and to assess the prognostic value of these pathways.
We searched PubMed, EMBASE, Cochrane, Web of Science, CNKI, Wangfang Data, and CQVIP from inception until July 2019. Eligible studies included the metabolomic analysis of blood samples from sepsis patients with the outcome. The metabolic pathway was assigned to each metabolite biomarker. The meta-analysis was performed using the pooled fold changes, area under the receiver operating characteristic curve (AUROC), and vote-counting of metabolic pathways. We also conducted a prospective cohort metabolomic study to validate the findings of our meta-analysis.
The meta-analysis included 21 cohorts reported in 16 studies with 2509 metabolite comparisons in the blood of 1287 individuals. We found highly limited overlap of the reported metabolite biomarkers across studies. However, these metabolites were enriched in several death-related metabolic pathways (DRMPs) including amino acids, mitochondrial metabolism, eicosanoids, and lysophospholipids. Prediction of sepsis death using DRMPs yielded a pooled AUROC of 0.81 (95% CI 0.76–0.87), which was similar to the combined metabolite biomarkers with a merged AUROC of 0.82 (95% CI 0.78–0.86) (P > 0.05). A prospective metabolomic analysis of 188 sepsis patients (134 survivors and 54 nonsurvivors) using the metabolites from DRMPs produced an AUROC of 0.88 (95% CI 0.78–0.97). The sensitivity and specificity for the prediction of sepsis death were 80.4% (95% CI 66.9–89.4%) and 78.8% (95% CI 62.3–89.3%), respectively.
DRMP analysis minimizes the discrepancies of results obtained from different metabolomic methods and is more practical than blood metabolite biomarkers for sepsis mortality prediction.
Sepsis is defined as the potentially life-threatening condition caused by the body’s extreme response to infection. Sepsis is one of the leading causes of death worldwide, and nearly 6 million people die of sepsis all over the world annually [1, 2]. Despite the increasing use of advanced technology for its treatment, such as bundled early goal-directed therapy (EGDT) , the prognosis of sepsis remains poor. The high mortality rates of sepsis are partially due to the lack of an effective approach to predict sepsis outcomes. In addition, only a small number of studies had investigated the molecular mechanisms of sepsis-induced organ failure and death [4, 5].
Several traditional sepsis outcome prediction approaches are currently used in clinical practice, such as the Sequential Organ Failure Assessment (SOFA) , the Acute Physiology and Chronic Health Evaluation II (APACHE II) , and the Simplified Acute Physiology Score II (SAPS II) . However, their performance (sensitivity, specificity) has not been shown to be adequate for all cases . There is a critical need to understand the mechanisms of sepsis-induced death and to identify better prognostic models that would facilitate the development of adapted strategies for different cases of sepsis.
The recently developed omics techniques facilitate high-throughput screening of disease-specific biomarkers in biological fluids, of which metabolomics is one of the most promising approaches . Metabolomics aims to measure the small molecules (metabolites) within cells, biofluids, tissues, or organisms using various analytical techniques. Unlike genomics, transcriptomics, and proteomics, metabolomics represents the molecular phenotype of an organism because metabolites and their concentrations are the direct functional “readout” of cellular activity and the state of cells and tissues . Therefore, clinical metabolomics offers a strategic advantage for the elucidation of the new roles of metabolism in disease, the identification of biomarkers, and the development of new therapeutics [11,12,13,14,15,16].
Recent studies have highlighted the potential prognostic role of metabolomics in sepsis patients [17,18,19]. Although promising, the existing literature of metabolomic studies on sepsis mortality prediction is limited by the use of individual cohorts with small sample sizes and low statistical power, as well as the use of varying analytical pipelines that can make it challenging to synthesize findings. Furthermore, the lack of validation in independent cohorts limits the clinical utility of metabolomic profiling in sepsis outcome prediction [20,21,22]. Meta-analysis can overcome these limitations by increasing the sample size and thus statistical power to generate the best estimation. Unfortunately, no studies have performed either descriptive or quantitative meta-analysis of metabolomics for predicting sepsis outcomes.
In this study, we systematically reviewed the literature to identify all eligible clinical metabolomic studies containing the prognosis of sepsis published before July 2019. We retrospectively generated a comprehensive dataset and performed both descriptive and quantitative meta-analysis using the curated dataset. In addition, we conducted a prospective metabolomic cohort study to validate the findings in the meta-analysis. The primary goals of this study were to identify the metabolic biomarkers and their pathways in the blood of sepsis nonsurvivors and to assess the prognostic significance of metabolomic profiling in sepsis patients.
We performed a comprehensive literature search of articles through the following databases without date limitation: PubMed, EMBASE, the Cochrane Library, Web of Science, China National Knowledge Infrastructure (CNKI), Wanfangdata, and CQVIP. The search was updated to July 1, 2019, and not restricted by language. The main search terms included the following: “Sepsis” (e.g., “Severe Sepsis” and “Sepsis, Severe” and “Pyemia” and “Pyemias” and “Pyohemia” and “Pyohemias” and “Pyaemia” and “Pyaemias” and “Septicemia” and “Septicemias” and “Poisoning, Blood” and “Blood Poisoning” and “Poisonings, Blood” and “Bacteremia” and “Endotoxemia” and “Fungemia” and “Candidemia” and “Parasitemia” and “Viremia”) and “Metabolomics” (e.g., “Metabolomic” and “Metabonomics” and “Metabonomic”). The reference list was also checked for relevant articles. The detailed search strategy is listed in Additional file 1: Supplemental methods.
Inclusion and exclusion criteria
The inclusion criteria for selecting the studies for this meta-analysis were as follows: (1) metabolomic profiling performed in patients with sepsis or septic shock according to the published diagnosis criteria , (2) metabolites were measured in plasma or serum, (3) sepsis survivors versus sepsis nonsurvivors, (4) mortality data and the correlation of metabolites with sepsis survivors and/or nonsurvivors were reported, and (5) the type of study design was not restricted. The exclusion criteria were as follows: (1) patients < 18 years old and studies involving animals, (2) studies without the sepsis outcome, and (3) meeting abstracts, letters to the editor, case reports, and reviews.
Data extraction and quality assessment
All potentially eligible articles were independently evaluated, and the information was extracted by two authors (JW and ST). Disagreements were resolved by discussion with a third person (YS). For each study, the following items were extracted: first author, year of publication, country, type of patients, total number of cases and gender, follow-up time, mortality, clinical scores (i.e., SOFA and APACHE II), the analytical platforms, metabolites with significant changes, metabolite fold change (FC, nonsurvivors/survivor), adjusted P value, area under receiver operating characteristic (ROC) curve (AUROC), and validation. Because of the inconsistency of the names for the same reported metabolites between publications, especially for the nomenclature of lipids, we then used the software OpenRefine (https://openrefine.org/) and the ID conversion tool in MetaboAnalyst 4.0 (www.metaboanalyst.ca) to match the names in the publications to the names in HMDB or PubChem. We removed the ambiguity by adding the identifiers to each metabolite if available (CAS, HMDB, or KEGG). We performed cross-data quality checks between reviewers at each step and reviewed all the included references after dataset construction.
Newcastle-Ottawa Scale (NOS) was used to assess the risk of bias for each of the included studies by two authors independently. The NOS has three parts: selection (0–4 points), comparability (0–2 points), and outcome assessment (0–3 points). We made slight modifications in the parts of selection and outcome assessment in original NOS based on the guidelines for a reliable metabolomic study [24, 25]. Both original and modified NOS scores are reported (Additional file 1: Table S1 and Table S2). Details for NOS and modified NOS are available in Additional file 1: Supplemental methods.
Data synthesis and meta-analysis
For pathway analysis, we assigned both the chemical class and the biochemical pathway to each metabolite according to HMDB and our in-house database. Cochran’s Q test and Higgins I-squared statistic were used to assess the heterogeneity across the included studies. The pooled AUROC and FC were calculated using a random model when the heterogeneity was statistically significant across the studies (P < 0.10 or I2 > 50%). Otherwise, a fixed-effect model was used. Sensitivity analyses were conducted by serially excluding each study to determine the influence of individual studies on the pooled AUROC. Publication bias was evaluated using funnel plots. The meta-analysis was performed using MedCalc 19.0.7. The frequency of a particular pathway (chemical class) reported to be significant across the studies was analyzed by the vote-counting method.
Validation of meta-analysis results using a prospective cohort metabolomic study
This prospective cohort study was approved by the Institutional Review Board (IRB) of Yantai Yuhuangding Hospital (11) and registered on the Chinese Clinical Trial Registry (Registry ID: ChiCTR1800015321). All the protocols conformed to the World Medical Association Declaration of Helsinki-Ethical Principles for Medical Research Involving Human Subjects. Briefly, the diagnostic evaluation was performed upon admission to the intensive care unit (ICU) of Yantai Yuhuangding Hospital, and patients who met the Third International Consensus Definitions for Sepsis and Septic Shock were eligible for selection . The exclusion criteria included the following: (1) younger than 18 years or older than 85 years of age, (2) diabetes and other metabolic-related diseases, (3) AIDS, (4) pregnant women, and (5) incomplete clinical data.
A total of 188 patients were enrolled between June 2017 and May 2018, and written permission was obtained from all the patients or their guardians. SOFA and APACHE II scores were assessed during the first 24 h of ICU admission. Heparinized plasma samples were collected at the time of ICU admission and stored at − 80 °C until analysis. The patients were followed up for the survival status weekly for 28 days. On day 28, out of 188 patients enrolled, 134 survived and 54 died. We then divided the plasma samples collected during admission into two groups: survivors (N = 134) and nonsurvivors (N = 54).
Metabolomic analysis was performed on extracted metabolites as described previously [27,28,29]. The chromatographic peaks were identified using MultiQuant (v3.0, AB SCIEX), confirmed by manual inspection, and the peak areas were integrated. The data were log2 transformed before statistical analysis. Partial least squares discriminant analysis (PLS-DA) was conducted in MetaboAnalyst 4.0. Seven metabolites including isoleucine (amino acids), alanine (amino acids), acetylcarnitine (mitochondrial metabolism), lactic acid (mitochondrial metabolism), pyruvic acid (mitochondrial metabolism), LysoPG (22:0) (lysophospholipid metabolism), and LysoPC (24:0) (lysophospholipid metabolism) were selected based on the results of the meta-analysis. Multivariate ROC analysis was conducted using MetaboAnalyst 4.0 (https://www.metaboanalyst.ca). The ROC curve was generated based on Monte Carlo cross-validation of random forest models . Repeated random cross-validation (rdCV) and permutation test were used for internal validation of the established classification model.
Other statistical analyses were conducted in GraphPad Prism 8.0. All statistical tests were 2-tailed, and the significance threshold (P) was set at 0.05.
Study inclusion and characteristics
The initial search strategies retrieved a total of 1814 articles. After meticulous inspection of the articles, 16 clinical metabolomic studies published between 2003 and 2019 were finally enrolled in our meta-analysis, in which sepsis survivors were compared to the sepsis nonsurvival group (Additional File 1: Table S3) [31,32,33,34,35,36,37,38,39,40,41,42,43,44,45,46]. The processes of study selection are summarized in a flow diagram (Fig. 1).
Studies related to sepsis mortality prediction using metabolomic approaches were mostly published after 2013, underscoring the emerging nature of this field of research. Among them, 7 studies were from the USA; 3 studies were performed in China; 2 studies were conducted in France and Canada, respectively; and 1 in Germany and Italy, respectively (Additional File 1: Table S3). APACHE and SOFA scores for mortality prediction were both reported in 7 studies. One study had both SOFA and SAPS scores, 2 studies reported SOFA scores, and 4 studies reported APACHE scores. Traditional prognostic scores were missing in 2 studies.
Sepsis is a severe critical illness syndrome with multi-organ dysfunction, and tissue biopsy is generally not performed for patients with sepsis. In addition, the urine metabolome shows more interindividual and technical variability than that of the blood . Therefore, in our meta-analysis, we chose only blood-based (plasma and serum) metabolomic studies (10 studies with plasma and 6 studies with serum) (Table 1). Overall, the resulting collection contains curated quality-checked data of 21 cohorts reported in 16 studies and over 2509 metabolite comparisons in blood from 1287 individuals (Additional file 1: Table S4).
Quality assessment of clinical metabolomic studies
Good quality control and proper validation are the prerequisites for the success of a clinical metabolomic study . The pooled QC samples were used in all MS-driven metabolomic studies of the prognosis of sepsis. Only 6 out of 16 studies reported validations either by animal studies, independent cohorts, or other omics techniques (Additional file 1: Table S5). Therefore, this meta-analysis is useful for validating the findings of individual studies. The details of quality assessment are described in Additional file 1: Supplemental results.
Meta-analysis of chemical classes and biochemical pathways using the vote-counting method
We first checked the overlap in metabolite biomarkers used for prediction of sepsis mortality across the studies but found few in common (Additional file 1: Figure S1). We then explored the unique metabolic features in sepsis nonsurvivors compared to the survivors by analyzing the chemical classes and metabolic pathways of the biomarkers.
We assigned the chemical class to each of the biomarkers based on the classification in HMDB (Additional file 1: Table S6). The main chemical classes for metabolites with prognostic values for sepsis are amino acid and derivatives, lipids and lipid-like molecules, and organic acids and derivatives (Fig. 2a).
We next analyzed the metabolic pathways for the differential metabolites between sepsis survivors and nonsurvivors (Fig. 2b). Interestingly, we found that these metabolite biomarkers are enriched in several death-related metabolic pathways (DRMPs). In detail, in sepsis nonsurvivors, lysophospholipid metabolism and phospholipid metabolism were significantly downregulated in nonsurvivors compared to the survivors (Fig. 2b and Additional file 1: Table S7). Acylcarnitines were observed to be increased dramatically in the plasma and serum of the sepsis nonsurvivors, indicating the downregulation of fatty acid oxidation. In addition, metabolites in the TCA cycle were increased in sepsis nonsurvivors. Together, these findings indicated the mitochondrial dysfunction in sepsis nonsurvivors. Furthermore, branched-chain amino acid metabolism, amino-sugar metabolism, and eicosanoids were dramatically upregulated in sepsis nonsurvivors (Fig. 2b).
Meta-analysis of AUROCs
ROC curve analysis is an effective method for evaluating the accuracy of diagnostic tests in modern medicine. Out of 16 studies enrolled, 12 conducted ROC curve analysis using the biomarker metabolites and six performed the direct comparisons between metabolite biomarkers and traditional scores (Additional file 1: Table S8). We next pooled the AUROCs from each study and performed a meta-analysis to analyze the prediction accuracy for sepsis prognosis by metabolomics. The sensitivity analysis revealed that the overall AUROC estimation and conclusions were not affected by each study.
Prediction performance of combined metabolite biomarkers without summarizing pathways
We first performed the meta-analysis of combined biomarker metabolites without summarizing pathways. Visual inspection of the funnel plot of the included studies revealed significant asymmetry (Additional file 1: Figure S2). The random-effect model was selected for this meta-analysis due to the heterogeneity (Q = 620.5, I2 = 95.2%, and P < 0.0001). The combined metabolite biomarkers yielded a pooled AUROC of 0.82 (95% CI 0.78–0.86, P < 0.001) (Fig. 3).
Prediction performance of each DRMP
We then investigated the prediction accuracy of each DRMP on sepsis outcome, including lysophospholipid metabolism, amino acid, and mitochondrial metabolism.
Lysophospholipid metabolism. There was no significant publication bias for the studies reporting lysophospholipids based on the visualization of the funnel plot (Additional file 1: Figure S3). Cochran’s Q was 3.45, and Higgins’ I2 was 13.5% (P = 0.33). The fixed-effect model was then used for this analysis. Meta-analysis showed a pooled AUROC of 0.77 (95% CI 0.72–0.82) for lysophospholipid metabolism (Fig. 4a).
Amino acid metabolism. Visual inspection of the funnel plot revealed significant asymmetry, suggesting the presence of potential publication bias in the studies using amino acids as biomarkers for sepsis mortality prediction (Additional file 1: Figure S4). Strong evidence of heterogeneity between the included studies was observed (Q = 177.04, I2 = 95.5%, and P < 0.01). The random-effect model was then chosen for the meta-analysis. Metabolites from amino acid metabolism produced a combined AUROC of 0.84 (95% CI 0.75–0.93) (Fig. 4b).
Mitochondrial metabolism. The funnel plot analysis indicated that there was no significant publication bias for the studies using metabolite biomarkers from mitochondrial metabolism (Additional file 1: Figure S5). Cochran’s Q was 27.75, and Higgins’ I2 was 82.0% (P < 0.001). We then used the random-effect model for the analysis. Meta-analysis showed a pooled AUROC of 0.78 (95% CI 0.71–0.86) for the prediction of sepsis death (Fig. 4c).
Prediction performance of the combined DRMPs
We next investigated the performance of metabolite biomarkers from combined death-related metabolic pathways. Funnel plot analysis revealed no obvious publication bias (Additional file 1: Figure S6). We observed significant evidence of heterogeneity between the included studies (Q = 212.2, I2 = 91.5%, and P < 0.0001), and the random-effect model was then chosen for the meta-analysis. A merged AUROC of 0.81 (95% CI 0.76–0.87) was obtained using the combined DRMPs (Fig. 5).
We compared the prediction accuracy for sepsis outcomes using different biomarkers. As shown in Fig. 6, the combined DRMPs have the similar prediction accuracy as merged metabolite biomarkers, and there was no significant difference in AUROCs between two approaches (P > 0.05).
Validation of the meta-analysis using a prospective metabolomic study
Finally, we conducted a prospective clinical study to confirm our findings in the meta-analysis. Metabolomic analysis was performed using the plasma samples collected from 188 sepsis patients, including 134 survivors and 54 nonsurvivors. The patients’ characteristics are listed in Table 2. The 3D PLS-DA plot revealed a distinct separation of the plasma metabolomic profiles between sepsis survivors and nonsurvivors (Additional file 1: Figure S7). We selected seven metabolites from the identified DRMPs including isoleucine (amino acid), alanine (amino acid), acetylcarnitine (mitochondrial metabolism), lactic acid (mitochondrial metabolism), pyruvic acid (mitochondrial metabolism), LysoPG (22:0) (lysophospholipid metabolism), and LysoPC (24:0) (lysophospholipid metabolism) as the potential biomarkers for sepsis outcome prediction. Interestingly, these metabolites were all statistically different between sepsis survivors and nonsurvivors, according to the Mann-Whitney U test (P < 0.05, 0.01, or 0.001, Additional file 1: Figure S8). ROC analysis of selected metabolites based on Monte Carlo cross-validation of random forest models yielded an AUROC of 0.88 (95% CI 0.78–0.97), which was significantly higher than the SOFA score (AUROC of 0.56, 95% CI 0.45–0.66) and APACHE II score (AUROC of 0.66, 95% CI 0.58–0.83) (Fig. 7). The accuracy, sensitivity, and specificity for the prediction of sepsis death were 80.1% (95% CI 69.2–88.0%), 80.4% (95% CI 66.9–89.4%), and 78.8% (95% CI 62.3–89.3%), respectively (Additional file 1: Table S9). The negative likelihood ratio (NLR), positive likelihood ratio (PLR), and diagnostic odds ratio (DOR) were 0.25 (95% CI 0.13–0.46), 3.79 (95% CI 1.93–7.43), and 15.2 (95% CI 5.04–46.2), respectively (Additional file 1: Table S9). Cross-validation revealed the rdCV accuracy of 77.3% (Additional file 1: Table S9) and permutation P value (500 times) of 0.0002 (Additional file 1: Figure S9).
Sepsis is one of the major causes of death in US hospitals. The identification of biomarkers that distinguish patients at high risk for poor outcomes will likely help the development of new treatment strategies. This meta-analysis of 21 cohorts with 1287 septic patients revealed that despite the heterogeneity of patients, there are unique metabolic signatures in sepsis nonsurvivors. Our study also showed that the analysis of DRMPs minimizes the discrepancies between metabolomic methods and is more practical than metabolite biomarkers themselves for prognostic prediction in patients with sepsis. To our knowledge, this is the first meta-analysis of the clinical utility of metabolomics for sepsis mortality prediction.
In this study, we made slight modifications for NOS based on the widely accepted guidelines for a reliable metabolomic study [24, 25]. For example, the sample size is a critical factor for metabolomic studies, and a good clinical metabolomic study generally requires at least 20 subjects per group to receive enough statistical power . We thus replaced question 4 in the section of “selection” in original NOS with sample size check in the modified NOS (Additional file 1: Supplemental methods and Table S1 and S2).
Our study identified death-related metabolic patterns in sepsis (Fig. 8). Briefly, infection-induced sepsis causes acute kidney injury (AKI) in the host , followed by ischemia and hypoxia in other organs such as the liver and lung. Acute respiratory distress syndrome (ARDS) in the lung contributes to the systemic metabolic responses [51, 52]. The aberrant metabolic responses of greatest importance to death in sepsis are mitochondrial dysfunction, breakdown of proteins and DNA, and uncontrolled inflammatory and immune responses, which produce DRMPs. The aggregation of metabolites in DRMPs leads to organ failure and eventually death.
The dramatic changes of metabolites in the serum and plasma of sepsis nonsurvivors include the increase of eicosanoids (eicosanoid metabolism), TCA cycle metabolites (lactate, pyruvate, and citric acid) (mitochondrial metabolism), acylcarnitines (mitochondrial metabolism), and amino acids (amino acid metabolism) and a sharp decrease of lysophospholipids (lysophospholipid metabolism).
Eicosanoid “storms” and uncontrolled inflammation
Eicosanoids are essential lipid mediators involved in the onset of inflammation and the innate immune responses . Our analysis revealed higher levels of circulating eicosanoids in sepsis nonsurvivors compared to the survivors. This observation was in line with the gene expression-based prognostic model developed by Sweeney and coauthors, in which genes from inflammation-related pathways were largely activated in sepsis nonsurvivors . Both hyperinflammation and immune hypoactivation lead to patient morbidity and death [55, 56]. The combination of multi-omic results might provide a comprehensive view of death-related changes at the molecular level in sepsis nonsurvivors.
Mitochondrial dysfunction in energy production
The accumulation of lactic acid, pyruvate, citric acid, and acylcarnitines in the serum and plasma of sepsis nonsurvivors indicates that death in sepsis might be associated with profound mitochondrial dysfunction in energy production. In contrast, sepsis survivors have better preservation of ATP, mitochondrial function, and biogenesis markers. Citric acid, lactic acid, and pyruvic acid are critical energetic substrates used by mitochondria for aerobic catabolism. Sepsis patients who ultimately die have poor aerobic catabolism due to mitochondrial dysfunction, displaying elevated concentrations of TCA cycle metabolites in plasma and serum. Acylcarnitines are essential for beta-oxidation of fatty acids and play essential roles in maintaining energy homeostasis in the human body. Acylcarnitines that are not utilized for energy production in fatty acid beta-oxidation due to the mitochondrial dysfunction will be reversely transported from mitochondria to the cytoplasm and then into the plasma and serum .
Systemic uncontrolled proteolysis
In our analysis, we found that sepsis nonsurvivors usually have significantly higher levels of circulating branched-chain amino acids in their plasma and serum. This difference is likely caused by systemic uncontrolled proteolysis in sepsis. Recent studies had suggested uncontrolled proteolysis as the fundamental pathological mechanism in septic shock that contributes to cell injury and organ dysfunction [58, 59]. Additionally, the transcriptomic analysis showed the upregulation of genes associated with proteasome degradation in sepsis nonsurvivors .
Profound defects of organ healing capability
In the nonsurvivors of sepsis patients, the levels of lysophospholipids were significantly lower than those in the survivors. Lysophosphatidic acid (LPA) had been reported to be overproduced in response to tissue injury, and it can promote healing in multiple organs such as the lung, skin, gastrointestinal tract, and cornea [60,61,62]. This suggested that the organ healing capability in nonsurvivors is lower than survivors.
Prognostic model for mortality prediction based on DRMPs
Our analysis revealed limited overlap in the reported metabolite biomarkers for sepsis mortality prediction across studies. This is likely caused by the discrepancies between analytical methods used for metabolomic analysis in these studies, such as different analytical instruments, instrument sensitivity, and metabolome coverage, and different sample preparation approaches. Certain metabolite biomarkers might not be well measured in one metabolomic method compared to another. However, these biomarkers were enriched in several DRMPs, and the changes in DRMPs are highly consistent no matter what metabolites are measured in these pathways. The DRMPs and metabolite biomarkers had a similar accuracy for predicting sepsis outcomes. Therefore, the analysis of the combined changes in DRMPs might be more practical than specific metabolite biomarkers themselves for prognosis prediction in patients with sepsis.
Our study has some limitations. First, currently, there are no established guidelines for the experimental design, analytical procedures, and data analysis of clinical metabolomic studies. The data were obtained from studies that had different experimental designs, analytical platforms, and patients’ characteristics. The substantial heterogeneity among the studies may affect the interpretation of the results. Second, due to the lack of full datasets for metabolomic analysis of sepsis , we performed the meta-analysis using the reported metabolites. This approach is then susceptible to publication bias.
Our meta-analysis of clinical metabolomic studies of sepsis prognosis reveals new roles of metabolism in sepsis-induced death and highlights the potential value of death-related metabolic pathways as the biomarkers in the prediction of sepsis mortality. These results will serve as a benchmark for future prognostic model development using metabolomics.
Availability of data and materials
All data generated or analyzed during this study were included in this published article and its additional file.
Acquired immune deficiency syndrome
Acute kidney injury
Analysis of variance
Acute Physiology and Chronic Health Evaluation
Acute respiratory distress syndrome
Area under receiver operating characteristic
Chemical Abstracts Service
China National Knowledge Infrastructure
Diagnostic odds ratio
Death-related metabolic pathways
Early goal-directed therapy
Human Metabolome Database
Intensive care unit
Institutional Review Board
Kyoto Encyclopedia of Genes and Genomes
Negative likelihood ratio
Positive likelihood ratio
Partial least squares discriminant analysis
Repeated random cross-validation
Receiver operating characteristic
Simplified Acute Physiology Score
Sequential Organ Failure Assessment
Tricarboxylic acid cycle
Fleischmann C, Scherag A, Adhikari NK, Hartog CS, Tsaganos T, Schlattmann P, Angus DC, Reinhart K, International Forum of Acute Care T. Assessment of global incidence and mortality of hospital-treated sepsis. Current estimates and limitations. Am J Respir Crit Care Med. 2016;193(3):259–72.
Kyriacou DN. Government regulation of sepsis care. JAMA. 2019;322(3):250–1.
Rivers EP, Coba V, Visbal A, Whitmill M, Amponsah D. Management of sepsis: early resuscitation. Clin Chest Med. 2008;29(4):689–704 ix-x.
Pool R, Gomez H, Kellum JA. Mechanisms of organ dysfunction in sepsis. Crit Care Clin. 2018;34(1):63–80.
Chalkias A, Xanthos T. Letter to the editor: Sepsis-associated in-hospital cardiac arrest: epidemiology, pathophysiology, and potential therapies. J Critical Care. 2017;40:314.
Innocenti F, Tozzi C, Donnini C, De Villa E, Conti A, Zanobetti M, Pini R. SOFA score in septic patients: incremental prognostic value over age, comorbidities, and parameters of sepsis severity. Intern Emerg Med. 2018;13(3):405–12.
Ho KM, Dobb GJ, Knuiman M, Finn J, Lee KY, Webb SA. A comparison of admission and worst 24-hour Acute Physiology and Chronic Health Evaluation II scores in predicting hospital mortality: a retrospective cohort study. Crit Care. 2006;10(1):R4.
Moreno RP, Metnitz PG, Almeida E, Jordan B, Bauer P, Campos RA, Iapichino G, Edbrooke D, Capuzzo M, Le Gall JR, et al. SAPS 3--from evaluation of the patient to evaluation of the intensive care unit. Part 2: development of a prognostic model for hospital mortality at ICU admission. Intensive Care Med. 2005;31(10):1345–55.
Walley KR. Biomarkers in sepsis. Curr Infect Dis Rep. 2013;15(5):413–20.
Liu X, Ren H, Peng D. Sepsis biomarkers: an omics perspective. Front Med. 2014;8(1):58–67.
Han J, Xia Y, Lin L, Zhang Z, Tian H, Li K. Next-generation metabolomics in the development of new antidepressants: using albiflorin as an example. Curr Pharm Des. 2018;24(22):2530–40.
Yu L, Li K, Zhang X. Next-generation metabolomics in lung cancer diagnosis, treatment and precision medicine: mini review. Oncotarget. 2017;8(70):115774–86.
Cui S, Li K, Ang L, Liu J, Cui L, Song X, Lv S, Mahmud E. Plasma phospholipids and sphingolipids identify stent restenosis after percutaneous coronary intervention. JACC Cardiovasc Interv. 2017;10(13):1307–16.
Koen N, Du Preez I, du Loots T. Metabolomics and personalized medicine. Adv Protein Chem Struct Biol. 2016;102:53–78.
Li S, Todor A, Luo R. Blood transcriptomics and metabolomics for personalized medicine. Comput Struct Biotechnol J. 2016;14:1–7.
Patti GJ, Tautenhahn R, Rinehart D, Cho K, Shriver LP, Manchester M, Nikolskiy I, Johnson CH, Mahieu NG, Siuzdak G. A view from above: cloud plots to visualize global metabolomic data. Anal Chem. 2013;85(2):798–804.
Kiehntopf M, Nin N, Bauer M. Metabolism, metabolome, and metabolomics in intensive care: is it time to move beyond monitoring of glucose and lactate? Am J Respir Crit Care Med. 2013;187(9):906–7.
Zurfluh S, Baumgartner T, Meier MA, Ottiger M, Voegeli A, Bernasconi L, Neyer P, Mueller B, Schuetz P. The role of metabolomic markers for patients with infectious diseases: implications for risk stratification and therapeutic modulation. Expert Rev Anti-Infect Ther. 2018;16(2):133–42.
Dos Santos CC. Shedding metabo‘light’ on the search for sepsis biomarkers. Crit Care. 2015;19:277.
Basoglu A, Sen I, Meoni G, Tenori L, Naseri A. NMR-based plasma metabolomics at set intervals in newborn dairy calves with severe sepsis. Mediat Inflamm. 2018;2018:8016510.
Whelan SP, Carchman EH, Kautza B, Nassour I, Mollen K, Escobar D, Gomez H, Rosengart MA, Shiva S, Zuckerbraun BS. Polymicrobial sepsis is associated with decreased hepatic oxidative phosphorylation and an altered metabolic profile. J Surg Res. 2014;186(1):297–303.
Neugebauer S, Giamarellos-Bourboulis EJ, Pelekanou A, Marioli A, Baziaka F, Tsangaris I, Bauer M, Kiehntopf M. Metabolite profiles in sepsis: developing prognostic tools based on the type of infection. Crit Care Med. 2016;44(9):1649–62.
Dellinger RP, Levy MM, Carlet JM, Bion J, Parker MM, Jaeschke R, Reinhart K, Angus DC, Brun-Buisson C, Beale R, et al. Surviving Sepsis Campaign: international guidelines for management of severe sepsis and septic shock: 2008. Crit Care Med. 2008;36(1):296–327.
Playdon MC, Joshi AD, Tabung FK, Cheng S, Henglin M, Kim A, Lin T, van Roekel EH, Huang J, Krumsiek J, et al. Metabolomics analytics workflow for epidemiological research: perspectives from the consortium of metabolomics studies (COMETS). Metabolites. 2019;9(7):E145.
Zeng X, Zhang Y, Kwong JS, Zhang C, Li S, Sun F, Niu Y, Du L. The methodological quality assessment tools for preclinical and clinical studies, systematic review and meta-analysis, and clinical practice guideline: a systematic review. J Evid Based Med. 2015;8(1):2–10.
Singer M, Deutschman CS, Seymour CW, Shankar-Hari M, Annane D, Bauer M, Bellomo R, Bernard GR, Chiche JD, Coopersmith CM, et al. The third international consensus definitions for sepsis and septic shock (Sepsis-3). JAMA. 2016;315(8):801–10.
Li K, Wang X, Pidatala VR, Chang CP, Cao X. Novel quantitative metabolomic approach for the study of stress responses of plant root metabolism. J Proteome Res. 2014;13(12):5879–87.
Yuan M, Breitkopf SB, Yang X, Asara JM. A positive/negative ion-switching, targeted mass spectrometry-based metabolomics platform for bodily fluids, cells, and fresh and fixed tissue. Nat Protoc. 2012;7(5):872–81.
Shi Y, Lin P, Wang X, Zou G, Li K. Sphingomyelin phosphodiesterase 1 (SMPD1) mediates the attenuation of myocardial infarction-induced cardiac fibrosis by astaxanthin. Biochem Biophys Res Commun. 2018;503(2):637–43.
Chong J, Wishart DS, Xia J. Using MetaboAnalyst 4.0 for comprehensive and integrative metabolomics data analysis. Curr Protoc Bioinformatics. 2019;68(1):e86.
Drobnik W, Liebisch G, Audebert FX, Frohlich D, Gluck T, Vogel P, Rothe G, Schmitz G. Plasma ceramide and lysophosphatidylcholine inversely correlate with mortality in sepsis patients. J Lipid Res. 2003;44(4):754–61.
Seymour CW, Yende S, Scott MJ, Pribis J, Mohney RP, Bell LN, Chen YF, Zuckerbraun BS, Bigbee WL, Yealy DM, et al. Metabolomics in pneumonia and sepsis: an analysis of the GenIMS cohort study. Intensive Care Med. 2013;39(8):1423–34.
Langley RJ, Tsalik EL, van Velkinburgh JC, Glickman SW, Rice BJ, Wang C, Chen B, Carin L, Suarez A, Mohney RP, et al. An integrated clinico-metabolomic model improves prediction of death in sepsis. Sci Transl Med. 2013;5(195):195ra195.
Su L, Huang Y, Zhu Y, Xia L, Wang R, Xiao K, Wang H, Yan P, Wen B, Cao L, et al. Discrimination of sepsis stage metabolic profiles with an LC/MS-MS-based metabolomics approach. BMJ Open Respir Res. 2014;1(1):e000056.
Rogers AJ, McGeachie M, Baron RM, Gazourian L, Haspel JA, Nakahira K, Fredenburgh LE, Hunninghake GM, Raby BA, Matthay MA, et al. Metabolomic derangements are associated with mortality in critically ill adult patients. PLoS One. 2014;9(1):e87538.
Mickiewicz B, Duggan GE, Winston BW, Doig C, Kubes P, Vogel HJ, Alberta SN. Metabolic profiling of serum samples by 1H nuclear magnetic resonance spectroscopy as a potential diagnostic approach for septic shock. Crit Care Med. 2014;42(5):1140–9.
Mickiewicz B, Tam P, Jenne CN, Leger C, Wong J, Winston BW, Doig C, Kubes P, Vogel HJ, Alberta SN. Integration of metabolic and inflammatory mediator profiles as a potential prognostic approach for septic shock in the intensive care unit. Crit Care. 2015;19:11.
Kamisoglu K, Haimovich B, Calvano SE, Coyle SM, Corbett SA, Langley RJ, Kingsmore SF, Androulakis IP. Human metabolic response to systemic inflammation: assessment of the concordance between experimental endotoxemia and clinical cases of sepsis/SIRS. Crit Care. 2015;19:71.
Su L, Li H, Xie A, Liu D, Rao W, Lan L, Li X, Li F, Xiao K, Wang H, et al. Dynamic changes in amino acid concentration profiles in patients with sepsis. PLoS One. 2015;10(4):e0121933.
Liu Z, Yin P, Amathieu R, Savarin P, Xu G. Application of LC-MS-based metabolomics method in differentiating septic survivors from non-survivors. Anal Bioanal Chem. 2016;408(27):7641–9.
Ferrario M, Cambiaghi A, Brunelli L, Giordano S, Caironi P, Guatteri L, Raimondi F, Gattinoni L, Latini R, Masson S, et al. Mortality prediction in patients with severe septic shock: a pilot study using a target metabolomics approach. Sci Rep. 2016;6:20391.
Mogensen KM, Lasky-Su J, Rogers AJ, Baron RM, Fredenburgh LE, Rawn J, Robinson MK, Massarro A, Choi AM, Christopher KB. Metabolites associated with malnutrition in the intensive care unit are also associated with 28-day mortality. JPEN J Parenter Enteral Nutr. 2017;41(2):188–97.
Dalli J, Colas RA, Quintana C, Barragan-Bradford D, Hurwitz S, Levy BD, Choi AM, Serhan CN, Baron RM. Human sepsis eicosanoid and proresolving lipid mediator temporal profiles: correlations with survival and clinical outcomes. Crit Care Med. 2017;45(1):58–68.
Wang L, Ko ER, Gilchrist JJ, Pittman KJ, Rautanen A, Pirinen M, Thompson JW, Dubois LG, Langley RJ, Jaslow SL, et al. Human genetic and metabolite variation reveals that methylthioadenosine is a prognostic biomarker and an inflammatory regulator in sepsis. Sci Adv. 2017;3(3):e1602096.
Chung KP, Chen GY, Chuang TY, Huang YT, Chang HT, Chen YF, Liu WL, Chen YJ, Hsu CL, Huang MT, et al. Increased plasma acetylcarnitine in sepsis is associated with multiple organ dysfunction and mortality: a multicenter cohort study. Crit Care Med. 2019;47(2):210–8.
Liu Z, Triba MN, Amathieu R, Lin X, Bouchemal N, Hantz E, Le Moyec L, Savarin P. Nuclear magnetic resonance-based serum metabolomic analysis reveals different disease evolution profiles between septic shock survivors and non-survivors. Crit Care. 2019;23(1):169.
Kim K, Mall C, Taylor SL, Hitchcock S, Zhang C, Wettersten HI, Jones AD, Chapman A, Weiss RH. Mealtime, temporal, and daily variability of the human urinary and plasma metabolomes in a tightly controlled environment. PLoS One. 2014;9(1):e86223.
Broadhurst D, Goodacre R, Reinke SN, Kuligowski J, Wilson ID, Lewis MR, Dunn WB. Guidelines and considerations for the use of system suitability and quality control samples in mass spectrometry assays applied in untargeted clinical metabolomic studies. Metabolomics. 2018;14(6):72.
Blaise BJ, Correia G, Tin A, Young JH, Vergnaud AC, Lewis M, Pearce JT, Elliott P, Nicholson JK, Holmes E, et al. Power analysis and sample size determination in metabolic phenotyping. Anal Chem. 2016;88(10):5179–88.
Rudiger A, Singer M. Acute kidney injury. Lancet. 2012;380(9857):1904.
Yadav H, Thompson BT, Gajic O. Fifty years of research in ARDS. Is acute respiratory distress syndrome a preventable disease? Am J Respir Crit Care Med. 2017;195(6):725–36.
Griffiths M, Proudfoot A. ARDS, up close and personal. Thorax. 2016;71(12):1073–5.
Dennis EA, Norris PC. Eicosanoid storm in infection and inflammation. Nat Rev Immunol. 2015;15(8):511–23.
Sweeney TE, Perumal TM, Henao R, Nichols M, Howrylak JA, Choi AM, Bermejo-Martin JF, Almansa R, Tamayo E, Davenport EE, et al. A community approach to mortality prediction in sepsis via gene expression analysis. Nat Commun. 2018;9(1):694.
Hotchkiss RS, Opal S. Immunotherapy for sepsis - a new approach against an ancient foe. New Engl J Med. 2010;363(1):87–9.
Said EA, Dupuy FP, Trautmann L, Zhang YW, Shi Y, El-Far M, Hill BJ, Noto A, Ancuta P, Peretz Y, et al. Programmed death-1-induced interleukin-10 production by monocytes impairs CD4(+) T cell activation during HIV infection. Nat Med. 2010;16(4):452–459.
Houten SM, Wanders RJ. A general introduction to the biochemistry of mitochondrial fatty acid beta-oxidation. J Inherit Metab Dis. 2010;33(5):469–77.
DeLano FA, Hoyt DB, Schmid-Schonbein GW. Pancreatic digestive enzyme blockade in the intestine increases survival after experimental shock. Sci Transl Med. 2013;5(169):169ra11.
Karnad DR, Bhadade R, Verma PK, Moulick ND, Daga MK, Chafekar ND, Iyer S. Intravenous administration of ulinastatin (human urinary trypsin inhibitor) in severe sepsis: a multicenter randomized controlled study. Intens Care Med. 2014;40(6):830–8.
Tager AM, LaCamera P, Shea BS, Campanella GS, Selman M, Zhao Z, Polosukhin V, Wain J, Karimi-Shah BA, Kim ND, et al. The lysophosphatidic acid receptor LPA1 links pulmonary fibrosis to lung injury by mediating fibroblast recruitment and vascular leak. Nat Med. 2008;14(1):45–54.
Liliom K, Guan Z, Tseng JL, Desiderio DM, Tigyi G, Watsky MA. Growth factor-like phospholipids generated after corneal injury. Am J Phys. 1998;274(4):C1065–74.
Demoyer JS, Skalak TC, Durieux ME. Lysophosphatidic acid enhances healing of acute cutaneous wounds in the mouse. Wound Repair Regen. 2000;8(6):530–7.
Sweeney TE, Khatri P. Generalizable biomarkers in critical care: toward precision medicine. Crit Care Med. 2017;45(6):934–9.
We would like to thank Dr. Jingbo Zhai from the Research Center of Tianjin University of Traditional Chinese Medicine (TCM) for his assistance on meta-analysis methodology.
This work was supported by the Shandong Provincial Natural Science Foundation (no. ZR2017MH075).
Ethics approval and consent to participate
This prospective cohort study was approved by the Institutional Review Board (IRB) of Yantai Yuhuangding Hospital (11) and registered on the Chinese Clinical Trial Registry (registry ID: ChiCTR1800015321). Written consent was obtained from all participants.
Consent for publication
The authors declare that they have no competing interests.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Supplemental methods. Supplemental results. Table S1. The risk of bias assessment of included studies using NOS. Table S2. The risk of bias assessment of included studies using modified NOS. Table S3. The clinical metabolomic studies selected in the meta-analysis. Table S4. Dataset summary for the included studies. Table S5. Quality assessment of clinical metabolomic studies for the mortality prediction of sepsis. Table S6. The identified metabolites biomarkers for the prediction of sepsis death and their chemical classes, metabolic pathways, Log2 fold change and P-value. Table S7. The vote count, pooled P-value, and Log2 fold change of the significantly altered metabolic pathways. Table S8. Studies containing the direct comparisons of prediction accuracy for sepsis outcomes between metabolomics and traditional scores. Table S9. Prediction accuracy of selected 7-analyte from DRMPs for sepsis death. Figure S1. Venn diagram showing the poor overlap of reported metabolite biomarkers across the studies for sepsis mortality prediction. Figure S2. The visualization of publication bias by the funnel plot for studies using metabolite biomarkers. Figure S3. Assessment of publication bias by funnel plot for studies using lysophospholipids as biomarkers for sepsis outcome prediction. Figure S4. Evaluation of publication bias by funnel plot for studies using amino acids as biomarkers for sepsis outcome prediction. Figure S5. Assessment of publication bias by funnel plot for studies using metabolites from mitochondrial metabolism as biomarkers for sepsis outcome prediction. Figure S6. Assessment of publication bias by funnel plot for studies using metabolites from DRMPs as biomarkers for sepsis outcome prediction. Figure S7. 3-D PLS-DA revealed the distinct separation of the plasma metabolome of sepsis nonsurvivors from sepsis survivors. Figure S8. The significant differences of seven selected biomarkers from DRMPs between sepsis survivors and nonsurvivors. Figure S9. The permutation test for the validation of the classification model robustness.
About this article
Cite this article
Wang, J., Sun, Y., Teng, S. et al. Prediction of sepsis mortality using metabolite biomarkers in the blood: a meta-analysis of death-related pathways and prospective validation. BMC Med 18, 83 (2020). https://doi.org/10.1186/s12916-020-01546-5
- Death-related metabolic pathways