- Open Access
Circulating N-formylmethionine and metabolic shift in critical illness: a multicohort metabolomics study
Critical Care volume 26, Article number: 321 (2022)
Cell stress promotes degradation of mitochondria which release danger-associated molecular patterns that are catabolized to N-formylmethionine. We hypothesized that in critically ill adults, the response to N-formylmethionine is associated with increases in metabolomic shift-related metabolites and increases in 28-day mortality.
We performed metabolomics analyses on plasma from the 428-subject Correction of Vitamin D Deficiency in Critically Ill Patients trial (VITdAL-ICU) cohort and the 90-subject Brigham and Women’s Hospital Registry of Critical Illness (RoCI) cohort. In the VITdAL-ICU cohort, we analyzed 983 metabolites at Intensive Care Unit (ICU) admission, day 3, and 7. In the RoCI cohort, we analyzed 411 metabolites at ICU admission. The association between N-formylmethionine and mortality was determined by adjusted logistic regression. The relationship between individual metabolites and N-formylmethionine abundance was assessed with false discovery rate correction via linear regression, linear mixed-effects, and Gaussian graphical models.
Patients with the top quartile of N-formylmethionine abundance at ICU admission had a significantly higher adjusted odds of 28-day mortality in the VITdAL-ICU (OR, 2.4; 95%CI 1.5–4.0; P = 0.001) and RoCI cohorts (OR, 5.1; 95%CI 1.4–18.7; P = 0.015). Adjusted linear regression shows that with increases in N-formylmethionine abundance at ICU admission, 55 metabolites have significant differences common to both the VITdAL-ICU and RoCI cohorts. With increased N-formylmethionine abundance, both cohorts had elevations in individual short-chain acylcarnitine, branched chain amino acid, kynurenine pathway, and pentose phosphate pathway metabolites.
The results indicate that circulating N-formylmethionine promotes a metabolic shift with heightened mortality that involves incomplete mitochondrial fatty acid oxidation, increased branched chain amino acid metabolism, and activation of the pentose phosphate pathway.
The study of the underlying biological processes of critical illness is hindered by patient heterogeneity . Defining critical illness based on common molecular mechanisms is useful for risk prediction and intervention response and can illuminate shared etiological and pathogenic pathways [2,3,4]. Common molecular mechanisms are described in specific populations of the critically ill using transcriptomics of COVID-19 , muscle transcriptomics and metabolomics in multiorgan failure , and genomics of sepsis due to pneumonia . But these observations are not generalizable beyond the specific population studied. We argue that the molecular response to evolutionarily conserved danger stimuli may reflect common pathological mechanisms in heterogeneous critically ill patients .
Danger-associated molecular patterns (DAMPs) are defense signals that communicate to the whole organism that individual cells have failed to adapt to stressors . In humans, degrading mitochondria yield DAMPs in the form of mitochondrial DNA and N-formyl-methionyl peptides. We and others have found that mitochondrial DAMPs are associated with significant dysregulation of metabolic homeostasis, potent activation of innate immunity, and adverse outcomes in the critically ill [10,11,12]. The further breakdown of the DAMP N-formyl-methionyl peptide (fMet-Leu-Phy) produces circulating N-formylmethionine that signals the substantial disruption of homeostasis .
Plasma metabolic profiles in the critically ill provide a way to measure point-in-time changes that reflect patients’ genetics as well as the cellular response to physiological stress, inflammation, and substrate abundance . But these metabolic profiles are not routinely used to provide common molecular insights into the response to critical illness. Therefore, we studied circulating N-formylmethionine abundance during early critical illness to identify patients with a metabolic shift who are at the greatest risk for poor outcomes. To test the hypothesis that common differences related to plasma metabolomic shift exist in the response to N-formylmethionine, we performed metabolomics cohort study on plasma samples from two distinct critically ill biorepositories with clinical heterogeneity.
Study design and patients
The Graz, Austria-based VITdAL-ICU trial randomized 492 critically ill adults with 25-hydroxyvitamin D (25(OH)D) ≤ 20 ng/ml to high-dose oral vitamin D3 or placebo. Patients were enrolled from medical, neurological, cardiac surgical, and anesthesia ICUs between 2010 and 2012 and followed for 180 days . Further details of the trial, consent, inclusion, and exclusion criteria are published and described in Additional file 1: Methods . Plasma was obtained at days 0 (prior to randomization), 3, and 7. The 428 VITdAL-ICU trial subjects with plasma samples available comprised the VITdAL-ICU cohort [3, 16].
The Boston, USA-based Registry of Critical Illness (RoCI) is a registry and biorepository of adult medical ICU patients with SIRS, sepsis, and sepsis ARDS at the Brigham and Women’s Hospital. Details of the cohort, patient recruitment, and consent are published [17, 18] and described in e-Appendix 1. Between 2008 and 2010, 90 medical ICU patients had plasma collected within 72 h of ICU admission (defined as day 0) and these patients comprise the RoCI cohort .
Metabolomics and study design
Metabolomics methods and analyses have been previously published [3, 16, 18] and are described in Additional file 1: Methods. In the VITdAL-ICU cohort, profiling for 983 metabolites in 1215 plasma samples from 428 subjects collected during the VITdAL-ICU trial was performed in 2017 [3, 15]. In the RoCI cohort, 411 metabolites were determined in 2011 on 90 patients . Metabolomics data from each cohort separately underwent cube root transformation and Pareto scaling to generate data on the same scale and to approximate normal distribution.
The exposure of interest was the relative abundance of N-formylmethionine determined at day 0 by liquid chromatography coupled mass spectrometry by Metabolon, Inc., and analyzed as a continuous variable. We chose to study N-formylmethionine as it is a known catabolic product of the danger-associated molecular pattern N-formyl-methionyl peptide and was measured in both the VITdAL-ICU cohort and RoCI cohort analyzed by Metabolon, Inc. Abundance is defined as the total ion count for the given mass-to-charge ratio assigned to N-formylmethionine. N-formylmethionine was also analyzed as quartiles of abundance with elevated N-formylmethionine defined as abundance in the top quartile. The primary outcome was all-cause 28-day mortality determined by hospital records and patient follow-up in the VITdAL-ICU cohort and hospital records and the US Social Security Administration Death Master File in the RoCI cohort. Distributions of and differences in crude survival between quartiles of N-formylmethionine were determined via the Kaplan–Meier survival curve and the log-rank test, respectively. The association between quartiles of N-formylmethionine and 28-day mortality was determined by logistic regression adjusted for age, sex, SAPS II, admission diagnosis, and baseline 25(OH)D in the VITdAL-ICU cohort and logistic regression adjusted for age, sex, race, and APACHE II score in the RoCI cohort.
For day 0 data in the VITdAL-ICU and RoCI cohorts, Student’s t test was performed to determine if significant differences exist with elevated N-formylmethionine (highest quartile) using MetaboAnalyst . We corrected for multiple testing via the Benjamini–Hochberg procedure to adjust the false discovery rate (FDR) to 0.05, producing a q value which was used to identify all significant differences . Day 0 data were also analyzed using orthogonal partial least square-discriminant analysis (OPLS-DA) to assess the significance of classification discrimination (SIMCA 15.0 Umetrics, Sweden) . OPLS-DA is a powerful statistical modeling tool that provides insights into metabolite differences between patients with and without elevated N-formylmethionine .
For single time point data in the VITdAL-ICU and RoCI cohorts, correlations between individual metabolites and N-formylmethionine abundance (continuous) at day 0 were determined utilizing linear regression models utilizing robust standard errors. For the VITdAL-ICU cohort, the linear regression model was corrected for age, sex, baseline 25(OH)D, SAPS II, and admission diagnosis. For the RoCI cohort, the linear regression model was adjusted for age, sex, race, and APACHE II. We performed linear regression diagnostics to evaluate model assumptions including variance inflation factor, Breusch–Pagan/Cook–Weisberg test, normality of residuals, and Cook's distance.
In the VITdAL-ICU cohort, correlations between N-formylmethionine and individual metabolites over days 0, 3, and 7 were also determined utilizing linear mixed-effects models correcting for age, sex, baseline 25(OH)D, absolute increase in 25(OH)D at day 3, SAPS II, plasma day, admission diagnosis, and an individual subject-specific random intercept. All linear regression and mixed-effects models were analyzed using STATA 16.1MP . A q value threshold of 0.05 was used to identify all significant associations . Rain plots were produced in R-3.6.2 .
Networks are useful in illustrating relationships in biological and pathophysiological pathways . To identify N-formylmethionine-specific metabolite connections, we estimated Gaussian graphical models (GGMs) for day 0 in the VITdAL-ICU cohort . Gaussian graphical models were used as they can demonstrate metabolic reactions and functional dependence between metabolites following adjustment for confounding factors . We inferred a N-formylmethionine-specific network for relative metabolite abundance including age, sex, SAPS II, admission diagnosis, and baseline 25(OH)D as covariates into the model . Edges between metabolites were allotted if both their Pearson correlations and partial correlations remained statistically significant at a q value threshold of 0.05 . Gaussian graphical models were produced using the GeneNet R package, version 1.2.13 in R-3.6.2 .
The baseline characteristics of the 428 subject VITdAL-ICU cohort are published and summarized in Table 1 [3, 16]. Baseline characteristics differed between subjects stratified by N-formylmethionine quartiles for age, sepsis, creatinine at day 0, and ICU type (Table 2; Additional file 2). The overall 28-day mortality of the 428 subject VITdAL-ICU cohort was 22%. Significant differences exist in survival between quartiles of N-formylmethionine abundance at day 0 (log-rank P < 0.001, Fig. 1A). Patients with the top quartile of plasma N-formylmethionine abundance at day 0 had a significantly higher odds of 28-day mortality (OR = 2.4 (95%CI 1.5–4.0; P = 0.001)) compared with patients in the other three quartiles, after adjustment for age, sex, 25(OH)D level, SAPS II, and admission diagnosis. Further adjustment for bilirubin, creatinine, or propofol did not materially alter the association between plasma N-formylmethionine and 28-day mortality (See Additional file 3).
In day 0 VITdAL-ICU cohort plasma samples, significant differences by t test exist for 517 individual metabolites (q value threshold of 0.05) in subjects with elevated N-formylmethionine relative to those without. t test results show that 430 metabolites were increased and 87 metabolites decreased in subjects with elevated N-formylmethionine relative to those without (See Additional file 4: Day 0 t test VITdAL-ICU tab). Selected plots of normalized unadjusted metabolite data highlight the day 0 increases in short-chain acylcarnitine (C2–C7), branched chain amino acids (BCAA), purine and pentose phosphate pathway metabolites all relative to N-formylmethionine quartiles (See Additional file 5). Regarding subject metabolomic profiles, though the multivariable OPLS-DA model has marginal predictability (Q2 = 0.35), the model is valid (y axis Q2 intercept of − 0.102) and is significantly different in patients with elevated N-formylmethionine compared to those without (CV-ANOVA p value < 0.001) (See Additional file 6) [28, 29].
In the VITdAL-ICU cohort linear regression analysis, there were significant associations at day 0 between abundance of 399 individual metabolites (q value threshold of 0.05) and increased N-formylmethionine abundance (analyzed as continuous) following adjustment for age, sex, SAPS II, admission diagnosis, and baseline 25(OH)D (See Additional file 7: Highlights and Additional file 4: Day 0 VITdAL-ICU tab for full data). Metabolites associated with increased N-formylmethionine at day 0 included elevated lactate, short-chain acylcarnitines, as well as BCAA, kynurenine pathway, and pentose phosphate pathway metabolites and decreased lysophosphatidylcholines. The same pattern was observed with further adjustment for bilirubin, creatinine, propofol exposure, or sepsis as well as removal of the SAPS II term (See Additional file 8).
Mixed-effects modeling of 1215 total days 0, 3, and 7 plasma samples from the VITdAL-ICU cohort shows 432 metabolites were significantly positively associated with increased N-formylmethionine abundance highlighted by increases in lactate, short-chain acylcarnitines, BCAA, purine and pentose phosphate pathway metabolites (See Additional file 9: Highlights and Additional file 4: Days 0, 3, and 7 VITdAL-ICU tab for full data). One hundred and three metabolites had significant negative associations with increased N-formylmethionine abundance, primarily by decreased lysophosphatidylcholines (See Additional file 9: Highlights and Additional file 4: Days 0, 3, and 7 VITdAL-ICU tab for full data). Similar patterns of association are found when the VITdAL-ICU cohort was limited to days 0, 3, and 7 data from the 216 VITdAL-ICU cohort subjects who received placebo (See Additional file 10: Highlights and Additional file 4: Days 0, 3, and 7 VITdAL-ICU Placebo tab for full data).
In the VITdAL-ICU cohort, we investigated N-formylmethionine-specific metabolic networks via Gaussian graphical models. The Gaussian graphical model analysis revealed eight N-formylmethionine-specific functional modules at day 0 (See Additional file 11). Similar to the day 0 linear regression analyses, short-chain acylcarnitine, BCAA, purine, kynurenine pathway, and pentose phosphate pathway metabolites were prominently featured in the N-formylmethionine-specific Gaussian graphical model analysis. Metabolites present in each functional module were all elevated with increased N-formylmethionine abundance and had biological or functional similarity.
Baseline characteristics of the 90-subject RoCI cohort were previously published and are again presented in Table 1 . Baseline characteristics differed between subjects stratified by N-formylmethionine quartile for APACHE II score, sepsis, and creatinine at day 0 (Table 2). The 28-day mortality of the RoCI cohort was 33%. There was a significantly worse survival with the top quartile of N-formylmethionine abundance (log-rank P < 0.001, Fig. 1B). Patients with the top quartile of plasma N-formylmethionine abundance had a fivefold higher odds of 28-day mortality (OR = 5.1 95%CI 1.4–18.7; P = 0.015) following adjustment for age, sex, race, and APACHE II score, compared with the other three quartiles. Further adjustment for bilirubin or creatinine did not materially alter the plasma N-formylmethionine mortality associations (See Additional file 3). In the 73 patients who had an existing measurement of plasma mitochondrial DNA levels [11, 30], there was a significant positive correlation between plasma mitochondrial DNA levels and N-formylmethionine abundance (Pearson correlation coefficient r2 = 0.47; P value < 0.001).
In day 0 RoCI cohort plasma samples, significant differences by t test were found for 127 individual metabolites (q value threshold of 0.05) in subjects with elevated N-formylmethionine relative to those without (See Additional file 4: Day 0 t test RoCI tab). Day 0 t test differences were present with increased individual short-chain acylcarnitine, BCAA, and pentose phosphate pathway metabolites and decreased lysophosphatidylcholine in subjects with the highest quartile of plasma N-formylmethionine abundance relative to the lower quartiles. Similar to the VITdAL-ICU cohort, the multivariable OPLS-DA RoCI cohort model is valid (Q2 intercept of -0.271) and has marginal predictability (Q2 = 0.35), but shows significant differences (CV-ANOVA P value < 0.001) in metabolomic profiles at day 0 in patients with highest quartile of plasma N-formylmethionine abundance relative to the lower quartiles (See Additional file 6) [28, 29].
In adjusted linear regression analysis of day 0 data from the RoCI cohort, significant differences exist in 82 individual metabolites (q value threshold 0.05) with increased N-formylmethionine abundance following adjustment for age, sex, race, and APACHE II. Day 0 differences were present with increased individual short-chain acylcarnitine, BCAA, and pentose phosphate pathway species in subjects with increased N-formylmethionine abundance (See Additional file 12: Highlights and Additional file 4: Day 0 RoCI for full data). Similar patterns were observed with further adjustment for bilirubin, creatinine, or sepsis (See Additional file 8).
Linear regression diagnostics
As determined by the variance inflation factor, in both cohorts the independent variables do not share a perfect, linear relationship. The Breusch–Pagan/Cook–Weisberg tests indicate that heteroskedasticity was present in both cohorts, and thus, robust standard errors were utilized in multiple linear regression. Individual metabolites at day 0 were shown to have approximately normal residuals when subjected to multiple linear regression with N-formylmethionine as the exposure in both the VITdAL-ICU and RoCI cohorts. Slight deviations from normal were generally noted at the upper tail. The residuals of xenobiotics (i.e., medications) flagrantly violated normality and thus linear regression estimates are not reported. In both cohorts, Cook's distance values were all less than 0.5, indicating that the multiple linear regression associations were not driven by outliers.
Metabolites common to both VITdAL-ICU and RoCI cohorts
In the 337 metabolites analyzed in both the VITdAL-ICU and RoCI cohorts, adjusted linear regression of day 0 plasma shows that with increases in N-formylmethionine abundance, 80 metabolites common to both the VITdAL-ICU and RoCI cohorts were changed (q value threshold of 0.10) with increases in N-formylmethionine abundance. Fifty-five metabolites common to both the VITdAL-ICU and RoCI cohorts had significant changes (q value threshold of 0.05) with increases in N-formylmethionine abundance. For each of the 55 individual metabolites, the direction of the relationship (increased or decreased) was identical in both cohorts (Table 3 for highlights; Additional file 13: Full data). With increased N-formylmethionine abundance, both cohorts had elevated individual short-chain acylcarnitine, BCAA, kynurenine pathway, and pentose phosphate pathway metabolites. The rain plot shows highlighted day 0 metabolites that increased in the VITdAL-ICU and validation cohorts relative to increases in N-formylmethionine abundance (Fig. 2). While differences in rain plot circle size reflective of study size exist for the VITdAL-ICU and RoCI cohorts, the effect size and direction of effect followed similar patterns for individual metabolites (Table 3).
With elevated physiologic stress, a metabolic shift in energy production occurs away from mitochondrial beta fatty oxidation to favor glycolysis for energy production, redox balance, and substrate need for biosynthesis [31, 32]. Although previous work suggests that endogenous danger signals are important in the critical illness response, our robust data argue that distinct nuanced energy metabolism patterns, and elevated mortality risks are present with increased N-formylmethionine abundance [11, 33]. We identified N-formylmethionine-specific metabolite networks using Gaussian graphical models which highlight the same groups identified in our regression analyses . The analyses in the VITdAL-ICU cohort and confirmed in the RoCI cohort highlight the potential of N-formylmethionine as a metabolic shift biomarker reflective of cellular performance and mortality outcome.
Though the response to critical illness is complex and requires highly demanding cellular functions, the initial signals of physiological stress are crude and widely conserved . Due to evolutionary resemblance, mitochondrial DNA and N-formyl-methionyl peptides act similarly to bacterial pathogen-associated molecular patterns [36, 37]. Circulating N-formylmethionine in critical illness is produced by the degradation of mitochondrial N-formyl-methionyl peptides and is a marker for mitochondrial-based protein turnover [13, 38]. Though not required for activation of N-formylpeptide receptor (FPR) subtypes, we argue that the abundance of N-formylmethionine is a biomarker of cell stress and signals ongoing disruption of homeostasis .
In both cohorts, we found evidence of a metabolic shift (Fig. 3). First, the observed rise in lactate with increased N-formylmethionine abundance indicates an increase in aerobic glycolysis or tissue hypoxia . Second, we show that with increased N-formylmethionine abundance, both cohorts have significant increases in modified purine nucleosides known to be liberated during oxidative stress . Third, our observation of increased circulating plasma short-chain acylcarnitines (C2–C7) with increased N-formylmethionine abundance is evidence of incomplete mitochondrial fatty acid β-oxidation and reflects a decrease in mitochondrial energy production . Fourth, our data show that circulating BCAA metabolites are elevated with increased N-formylmethionine abundance. BCAA metabolites are known to be utilized for energy production when mitochondrial fatty acid β-oxidation is incomplete . Further, the observed elevated C3- and C5-acylcarnitines are known to be derived from BCAA metabolites . In our study cohorts, the observed increases in plasma lactate, purine nucleotides, short-chain acylcarnitines, and BCAA catabolic metabolites with increased N-formylmethionine abundance are thus indicative of a metabolic shift.
Our results show novel increases in pentose phosphate pathway metabolites with increased N-formylmethionine abundance. The pentose phosphate pathway is important for redox balance via NADPH production through a oxidative branch as well as nucleic acid synthesis through a non-oxidative branch (See Additional file 14) . Though quiescent at homeostasis and rarely reported in critical illness, evidence suggests that the pentose phosphate pathway acts as a “metabolic redox sensor” that is upregulated during oxidative stress [46, 47]. Transcriptomic profiling of whole blood in critically ill patients has shown increases in pentose phosphate pathway oxidative branch activity, suggesting a metabolic shift away from mitochondrial beta oxidation . We show that with increased N-formylmethionine abundance, both cohorts have evidence of significant increases in the pentose phosphate pathway metabolites of the non-oxidative branch, which is further indicative of a metabolic shift (Table 3).
Our novel study approach has several strengths. Both of our cohorts were analyzed using the same Metabolon platform. Further, our use of clinical data allows for modeling of metabolite abundance with adjustment for subject characteristics. Our findings are consistent with known metabolic shift patterns increasing the biological plausibility of our work. Finally, we validated our findings in an external critical illness cohort.
There are potential limitations to our approach. Despite our RoCI cohort being diverse, our VITdAL-ICU cohort subjects are from a white population. Our ability to determine causality is limited as we performed two metabolomics cohort studies. Our approach is subject to bias and confounding despite our covariate and multiplicity adjustment. With only 337 metabolites common to both cohorts, some important individual metabolites noted in the VITdAL-ICU cohort cannot be analyzed in the RoCI cohort. As we did not measure the N-formylmethionine by concentration, our abundance quartiles are unlikely to have the same N-formylmethionine concentration in both the VITdAL-ICU and RoCI cohorts. Thus, our comparison of top quartiles of abundance between the two cohorts cannot provide an accurate cut point for the metabolic shift responses to cell stress. Further, we do not have extensive data on nutritional support which may confound the metabolite differences observed. Finally, our similar metabolic shift findings in our two cohorts should be considered as a hypothesis that requires further confirmation and careful interpretation.
Metabolomics can bridge the knowledge gap between clinical heterogeneity and severity of critical illness . Finding consistent metabolic shift patterns in clinically heterogeneous cohorts has the potential to identify biologically meaningful phenotypes . We find that with severe cell stress indicated by elevated circulating N-formylmethionine, a shared metabolic shift exists in both the VITdAL-ICU and RoCI cohorts despite clinical heterogeneity. As different clinical presentations can produce different metabolomic patterns at different times, the finding of a shared metabolic shift associated with cell stress may be helpful in detecting a severity of illness phenotype . Next research steps include confirming our metabolic profiling and mortality results with determination of N-formylmethionine concentrations via ELISA . To be generalizable and robust, such further work will require multiple cohorts with clinical heterogeneity and adequate sample sizes. If confirmed, such pathobiology-driven study of the clinical heterogeneity of severe critical illness has future clinical implications in providing insights beyond the limitations of illness severity scores and outcome prediction models .
Our findings offer a deep and nuanced dive into the early critical illness response. In addition to elevated mortality, we found that patients with elevated N-formylmethionine have increases in metabolites related to short-chain acylcarnitines, BCAA catabolism, and the pentose phosphate pathway. The metabolomic patterns related to energy utilization with increases in N-formylmethionine allow for a common metabolomic shift to be identified in the heterogeneous critically ill population. Defining a common metabolic shift is a crucial first step toward a structural framework to increase mechanistic understanding of early critical illness.
Availability of data and materials
The datasets supporting the conclusions of this article are included within the article and its additional files.
Danger-associated molecular patterns
Intensive Care Unit
- VITdAL-ICU trial:
Correction of Vitamin D Deficiency in Critically Ill Patients trial
The Registry of Critical Illness
Orthogonal partial least square-discriminant analysis
Branched chain amino acid
Gaussian graphical models
Acute Physiology and Chronic Health Evaluation
Simplified Acute Physiology Score
Analysis of variance of the cross-validated residuals
Scicluna BP. Resolving patient heterogeneity in critical illness requires multi-scale approaches. EBioMedicine. 2022;77:103918.
Braga D, Barcella M, Herpain A, Aletti F, Kistler EB, Bollen Pinto B, Bendjelid K, Barlassina C. A longitudinal study highlights shared aspects of the transcriptomic response to cardiogenic and septic shock. Crit Care. 2019;23(1):414.
Amrein K, Lasky-Su JA, Dobnig H, Christopher KB. Metabolomic basis for response to high dose vitamin D in critical illness. Clin Nutr. 2021;40(4):2053–60.
Initiative C-H. Mapping the human genetic architecture of COVID-19. Nature. 2021;600(7889):472–7.
Frishberg A, Kooistra E, Nuesch-Germano M, Pecht T, Milman N, Reusch N, Warnat-Herresthal S, Bruse N, Handler K, Theis H, et al. Mature neutrophils and a NF-kappaB-to-IFN transition determine the unifying disease recovery dynamics in COVID-19. Cell Rep Med. 2022;2022:100652.
Carre JE, Orban JC, Re L, Felsmann K, Iffert W, Bauer M, Suliman HB, Piantadosi CA, Mayhew TM, Breen P, et al. Survival in critical illness is associated with early activation of mitochondrial biogenesis. Am J Respir Crit Care Med. 2010;182(6):745–51.
Rautanen A, Mills TC, Gordon AC, Hutton P, Steffens M, Nuamah R, Chiche JD, Parks T, Chapman SJ, Davenport EE, et al. Genome-wide association study of survival from sepsis due to pneumonia: an observational cohort study. Lancet Respir Med. 2015;3(1):53–60.
Russell CD, Baillie JK. Treatable traits and therapeutic targets: Goals for systems biology in infectious disease. Curr Opin Syst Biol. 2017;2:140–6.
Chen GY, Nunez G. Sterile inflammation: sensing and reacting to damage. Nat Rev Immunol. 2010;10(12):826–37.
Nakahira K, Hisata S, Choi AM. The roles of mitochondrial damage-associated molecular patterns in diseases. Antioxid Redox Signal. 2015;23(17):1329–50.
Johansson PI, Nakahira K, Rogers AJ, McGeachie MJ, Baron RM, Fredenburgh LE, Harrington J, Choi AMK, Christopher KB. Plasma mitochondrial DNA and metabolomic alterations in severe critical illness. Crit Care. 2018;22(1):360.
Faust HE, Reilly JP, Anderson BJ, Ittner CAG, Forker CM, Zhang P, Weaver BA, Holena DN, Lanken PN, Christie JD, et al. Plasma mitochondrial DNA levels are associated with ARDS in trauma and sepsis patients. Chest. 2020;157(1):67–76.
Raabe CA, Groper J, Rescher U. Biased perspectives on formyl peptide receptors. Biochim Biophys Acta Mol Cell Res. 2019;1866(2):305–16.
Johnson CH, Ivanisevic J, Siuzdak G. Metabolomics: beyond biomarkers and towards mechanisms. Nat Rev Mol Cell Biol. 2016;17(7):451–9.
Amrein K, Schnedl C, Holl A, Riedl R, Christopher KB, Pachler C, Urbanic Purkart T, Waltensdorfer A, Munch A, Warnkross H, et al. Effect of high-dose vitamin D3 on hospital length of stay in critically ill patients with vitamin D deficiency: the VITdAL-ICU randomized clinical trial. JAMA. 2014;312(15):1520–30.
Chary S, Amrein K, Lasky-Su J, Dobnig H, Christopher KB. The sex-specific metabolic response to critical illness: a post-hoc metabolomics study of the VITdAL-ICU trial. Sci Rep. 2021;11:3951.
Dolinay T, Kim YS, Howrylak J, Hunninghake GM, An CH, Fredenburgh L, Massaro AF, Rogers A, Gazourian L, Nakahira K, et al. Inflammasome-regulated cytokines are critical mediators of acute lung injury. Am J Respir Crit Care Med. 2012;185(11):1225–34.
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.
Chong J, Xia J. Using MetaboAnalyst 4.0 for metabolomics data analysis, interpretation, and integration with other omics data. Methods Mol Biol. 2020;2104:337–60.
Benjamini Y, Hochberg Y. Controlling for false discovery rate: a practical and powerful approach to multiple testing. J R Stat Soc Ser B (Methodol). 1995;57:289–300.
Matthew B, William R. Partial least squares for discrimination. J Chemom. 2003;17(3):166–73.
Worley B, Powers R. PCA as a practical indicator of OPLS-DA model reliability. Curr Metab. 2016;4(2):97–103.
StataCorp: Stata Statistical Software: Release 16. In. College Station, TX: StataCorp LP (2019)
Henglin M, Niiranen T, Watrous JD, Lagerborg KA, Antonelli J, Claggett BL, Demosthenes EJ, von Jeinsen B, Demler O, Vasan RS, et al. A single visualization technique for displaying multiple metabolite-phenotype associations. Metabolites. 2019;9(7):128.
Zhang B, Tian Y, Zhang Z. Network biology in medicine and beyond. Circ Cardiovasc Genet. 2014;7(4):536–47.
Krumsiek J, Suhre K, Illig T, Adamski J, Theis FJ. Gaussian graphical modeling reconstructs pathway reactions from high-throughput metabolomics data. BMC Syst Biol. 2011;5:21.
Do KT, Pietzner M, Rasp DJ, Friedrich N, Nauck M, Kocher T, Suhre K, Mook-Kanamori DO, Kastenmuller G, Krumsiek J. Phenotype-driven identification of modules in a hierarchical map of multifluid metabolic correlations. NPJ Syst Biol Appl. 2017;3:28.
Westerhuis JA, Hoefsloot HCJ, Smit S, Vis DJ, Smilde AK, van Velzen EJJ, van Duijnhoven JPM, van Dorsten FA. Assessment of PLSDA cross validation. Metabolomics. 2008;4(1):81–9.
Eriksson L, Trygg J, Wold S. CV-ANOVA for significance testing of PLS and OPLS models. J Chemom. 2008;22(11–12):594–600.
Nakahira K, Kyung SY, Rogers AJ, Gazourian L, Youn S, Massaro AF, Quintana C, Osorio JC, Wang Z, Zhao Y, et al. Circulating mitochondrial DNA in patients in the ICU as a marker of mortality: derivation and validation. PLoS Med. 2013;10(12):e1001577.
Preiser JC, Ichai C, Orban JC, Groeneveld AB. Metabolic response to the stress of critical illness. Br J Anaesth. 2014;113(6):945–54.
Puthucheary ZA, Astin R, McPhail MJW, Saeed S, Pasha Y, Bear DE, Constantin D, Velloso C, Manning S, Calvert L, et al. Metabolic phenotype of skeletal muscle in early critical illness. Thorax. 2018;73(10):926–35.
Venereau E, Ceriotti C, Bianchi ME. DAMPs from cell death to new life. Front Immunol. 2015;6:422.
Mittelstrass K, Ried JS, Yu Z, Krumsiek J, Gieger C, Prehn C, Roemisch-Margl W, Polonikov A, Peters A, Theis FJ, et al. Discovery of sexual dimorphisms in metabolic and genetic biomarkers. PLoS Genet. 2011;7(8):e1002215.
Jacob F (1970) La Logique du vivant. Une histoire de l’hérédité. Paris: Éditions Gallimard
Hemmi H, Takeuchi O, Kawai T, Kaisho T, Sato S, Sanjo H, Matsumoto M, Hoshino K, Wagner H, Takeda K, et al. A Toll-like receptor recognizes bacterial DNA. Nature. 2000;408(6813):740–5.
Waller JP. The Nh2-terminal residues of the proteins from cell-free extracts of E. coli. J Mol Biol. 1963;7:483–96.
Garcia-Contreras M, Tamayo-Garcia A, Pappan KL, Michelotti GA, Stabler CL, Ricordi C, Buchwald P. Metabolomics study of the effects of inflammation, hypoxia, and high glucose on isolated human pancreatic islets. J Proteome Res. 2017;16(6):2294–306.
Gao JL, Becker EL, Freer RJ, Muthukumaraswamy N, Murphy PM. A high potency nonformylated peptide agonist for the phagocyte N-formylpeptide chemotactic receptor. J Exp Med. 1994;180(6):2191–7.
Levy B, Desebbe O, Montemont C, Gibot S. Increased aerobic glycolysis through beta2 stimulation is a common mechanism involved in lactate formation during shock states. Shock. 2008;30(4):417–21.
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.
Koves TR, Ussher JR, Noland RC, Slentz D, Mosedale M, Ilkayeva O, Bain J, Stevens R, Dyck JR, Newgard CB, et al. Mitochondrial overload and incomplete fatty acid oxidation contribute to skeletal muscle insulin resistance. Cell Metab. 2008;7(1):45–56.
Neinast M, Murashige D, Arany Z. Branched chain amino acids. Annu Rev Physiol. 2019;81:139–64.
Newgard CB. Interplay between lipids and branched-chain amino acids in development of insulin resistance. Cell Metab. 2012;15(5):606–14.
Kletzien RF, Harris PK, Foellmi LA. Glucose-6-phosphate dehydrogenase: a “housekeeping” enzyme subject to tissue-specific regulation by hormones, nutrients, and oxidant stress. FASEB J. 1994;8(2):174–81.
Kruger A, Gruning NM, Wamelink MM, Kerick M, Kirpy A, Parkhomchuk D, Bluemlein K, Schweiger MR, Soldatov A, Lehrach H, et al. The pentose phosphate pathway is a metabolic redox sensor and regulates transcription during the antioxidant response. Antioxid Redox Signal. 2011;15(2):311–24.
Haji-Michael PG, Ladriere L, Sener A, Vincent JL, Malaisse WJ. Leukocyte glycolysis and lactate output in animal sepsis and ex vivo human blood. Metabolism. 1999;48(6):779–85.
Nalos M, Parnell G, Robergs R, Booth D, McLean AS, Tang BM. Transcriptional reprogramming of metabolic pathways in critically ill patients. Intensive Care Med Exp. 2016;4(1):21.
Thistlethwaite LR, Li X, Burrage LC, Riehle K, Hacia JG, Braverman N, Wangler MF, Miller MJ, Elsea SH, Milosavljevic A. Clinical diagnosis of metabolic disorders using untargeted metabolomic profiling and disease-specific networks learned from profiling data. Sci Rep. 2022;12(1):6556.
Patti GJ, Yanes O, Siuzdak G. Innovation: Metabolomics: the apogee of the omics trilogy. Nat Rev Mol Cell Biol. 2012;13(4):263–9.
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.
Duvvuri B, Baddour AA, Deane KD, Feser ML, Nelson JL, Demoruelle MK, Lood C. Mitochondrial N-formyl methionine peptides associate with disease activity as well as contribute to neutrophil activation in patients with rheumatoid arthritis. J Autoimmun. 2021;119:102630.
Desai N, Gross J. Scoring systems in the critically ill: uses, cautions, and future directions. BJA Educ. 2019;19(7):212–8.
This work is dedicated to the memory of our dear friend and colleague Nathan Edward Hellman, MD, PhD. The authors thank Shawn Murphy and Henry Chueh and the Mass General Brigham Research Patient Data Registry group for facilitating use of their database. Graphical abstract and Figure 3 are created with the BioRender.com application.
KN is supported by Foundation for the National Institutes of Health (NIH)/National Center for Advancing Translational Sciences grant KL2-TR-002385, R01 HL123915. AJR is supported by NIH grant R01 HL152083. LEF is supported by NIH grant R01 HL114839. RMB is supported by NIH grants R01 HL142093 and R01 GM115605. KBC is supported by NIH grant R01 GM115774. The VITdAL-ICU trial was supported by the European Society for Clinical Nutrition and Metabolism (ESPEN), a research grant including provision of study medication from Fresenius Kabi (Germany), and the Austrian National Bank (Jubiläumsfonds, Project Nr. 14143). Landspitali University Hospital Science Fund: A2021-03
Ethics approval and consent to participate
The study was conducted according to the guidelines of the Declaration of Helsinki and approved by Mass General Brigham Human Research Committee (Protocol numbers 2008P000495 and 2015P002766) and the Ethikkommission der Medizinischen Universität Graz, Austria (Ethikkommission number 21-214). Prior to the Registry of Critical Illness cohort enrollment, written informed consent was obtained which included permission for blood samples to be utilized and saved for research purposes in the future. At VITdAL-ICU trial enrollment, written informed consent was obtained, if possible, directly from the patient or from a legal surrogate. VITdAL-ICU trial consent included permission for blood specimens to be saved for future research studies.
Consent for publication
Dr. Amrein reports receiving lecture fees from Fresenius Kabi. Dr. Baron has served on Advisory Boards for Merck and Genentech. Drs. Sigurdsson, Kobayashi, Lasky-Su, Nakahira, Rogers, Pinilla-Vera, Baron, Fredenburgh, and Christopher report no financial or other relationships that might lead to a conflict of interest.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Additional file 1.
Additional file 2.
Admission Diagnosis category of VITdAL-ICU cohort patients by day 0 N-Formylmethionine abundance quartiles.
Additional file 3.
Associations between plasma N-formylmethionine abundance quartiles and 28-day mortality. Logistic regression models shown for crude and adjusted VITdAL-ICU and RoCI cohort survival data. Additional models shown with adjustment for total bilirubin, creatinine, and/or propofol exposure.
Additional file 4.
Full metabolomics results from VITdAL-ICU and RoCI Cohorts. Six Individual Worksheets include Day 0 t test VITdAL-ICU: t test results of Day 0 data from the VITdAL-ICU cohort; Day 0 VITdAL-ICU: Linear regression results of Day 0 data from the VITdAL-ICU cohort adjusted for age, sex, baseline 25(OH)D, SAPS II, and admission diagnosis; Day 0,3,7 VITdAL-ICU: Mixed effects regression results of Day 0,3,7 data from the VITdAL-ICU cohort adjusted for age, sex, baseline 25(OH)D, absolute increase in 25(OH)D at day 3, SAPS II, plasma day, admission diagnosis, and an individual subject-specific random-intercept; Day 0,3,7 VITdAL-ICU Placebo: Mixed effects regression results of Day 0,3,7 data from the VITdAL-ICU cohort in patients who received Placebo adjusted for age, sex, baseline 25(OH)D, absolute increase in 25(OH)D at day 3, SAPS II, plasma day, admission diagnosis, and an individual subject-specific random-intercept; Day 0 t test RoCI: t test results of Day 0 data from the RoCI cohort; Day 0 RoCI: Linear regression results of Day 0 data from the RoCI cohort adjusted for age, sex, race, and APACHE II score. T-statistics are presented for t tests. Beta coefficients are presented for linear regression and mixed effects analysis. In all worksheets, raw p values are presented. A false discovery rate adjusted p value (q value) threshold of 0.05 was used to identify all significant differences.
Additional file 5.
Dot Boxplots of metabolite Sub Pathways by N-formylmethionine quartiles at Day 0. Unadjusted normalized abundance of metabolites from 428 VITdAL-ICU cohort subjects at day 0 by quartiles of N-formylmethionine abundance. Sub Pathways shown include short-chain acylcarnitine, branched chain amino acid, pentose phosphate pathway and purine metabolites. Box plots show the data between the first and third quartile (box), the median and whiskers at 1.5 times the interquartile range. All individual data points are visualized using a bee swarm plot.
Additional file 6.
Day 0 OPLS-DA model goodness of fit, predictive ability and model significance in both VITdAL-ICU and RoCI Cohorts. R2 and Q2 quality metrics, the permutation diagnostics and overall model significance are presented.
Additional file 7.
Highlighted Significantly Different Metabolites with increased N-formylmethionine abundance in the VITdAL-ICU Cohort at day 0.
Additional file 8.
Additional metabolomics results from VITdAL-ICU and RoCI Cohorts. Nine individual worksheets include Linear regression results of Day 0 data from the VITdAL-ICU cohort adjusted for age, sex, baseline 25(OH)D, SAPS II, admission diagnosis, with additional adjustment for total bilirubin, creatinine, sepsis or propofol exposure. Day 0 data from the VITdAL-ICU cohort is also shown adjusted for age, sex, baseline 25(OH)D, and admission diagnosis but without SAPS II. Linear regression results of Day 0 data from the RoCI cohort are shown adjusted for age, sex, baseline 25(OH)D, SAPS II, admission diagnosis, with additional adjustment for total bilirubin, creatinine, or sepsis. In all worksheets, raw p values are presented. A false discovery rate adjusted p value (q value) threshold of 0.05 was used to identify all significant differences.
Additional file 9.
Metabolites significantly changed with increased N-formylmethionine abundance in the VITdAL-ICU Cohort over days 0, 3 and 7.
Additional file 10.
Metabolites significantly changed with increased N-formylmethionine abundance in VITdAL-ICU Cohort placebo patients over days 0, 3 and 7.
Additional file 11.
Day 0 N-formylmethionine-specific Metabolic Networks with similar effects via Gaussian graphical models.
Additional file 12.
Significantly Different Metabolites with increased N-formylmethionine abundance in the RoCI Cohort at day 0.
Additional file 13.
Significantly Different Metabolites with increased N-formylmethionine at day 0 in both VITdAL-ICU and RoCI Cohorts.
Additional file 14.
Increased Pentose Phosphate Pathway Metabolites with N-formylmethionine. Metabolites increased in the VITdAL-ICU and RoCI cohorts with increased N-formylmethionine abundance are highlighted in red, metabolites increased in the VITdAL-ICU cohort with increased N-formylmethionine abundance are highlighted in green. The Pentose Phosphate Pathway produces NADPH and ribose 5-phosphate for redox regulation and biosynthesis, respectively.
Rights and permissions
Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecommons.org/licenses/by/4.0/. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated in a credit line to the data.
About this article
Cite this article
Sigurdsson, M.I., Kobayashi, H., Amrein, K. et al. Circulating N-formylmethionine and metabolic shift in critical illness: a multicohort metabolomics study. Crit Care 26, 321 (2022). https://doi.org/10.1186/s13054-022-04174-y
- Critical illness
- Metabolic shift
- Pentose phosphate pathway
- Branched chain amino acids