Skip to main content

A prospective pilot study using metabolomics discloses specific fatty acid, catecholamine and tryptophan metabolic pathways as possible predictors for a negative outcome after severe trauma



We wanted to define metabolomic patterns in plasma to predict a negative outcome in severe trauma patients.


A prospective pilot study was designed to evaluate plasma metabolomic patterns, established by liquid chromatography coupled to mass spectrometry, in patients allocated to an intensive care unit (in the University Hospital Arnau de Vilanova, Lleida, Spain) in the first hours after a severe trauma (n = 48). Univariate and multivariate statistics were employed to establish potential predictors of mortality.


Plasma of patients non surviving to trauma (n = 5) exhibited a discriminating metabolomic pattern, involving basically metabolites belonging to fatty acid and catecholamine synthesis as well as tryptophan degradation pathways. Thus, concentration of several metabolites exhibited an area under the receiver operating curve (ROC) higher than 0.84, including 3-indolelactic acid, hydroxyisovaleric acid, phenylethanolamine, cortisol, epinephrine and myristic acid. Multivariate binary regression logistic revealed that patients with higher myristic acid concentrations had a non-survival odds ratio of 2.1 (CI 95% 1.1–3.9).


Specific fatty acids, catecholamine synthesis and tryptophan degradation pathways could be implicated in a negative outcome after trauma. The metabolomic study of severe trauma patients could be helpful for biomarker proposal.


Severe trauma injuries are the first cause of mortality among population younger than 44. With a global mortality higher than 5 million individuals/year its clinical assessment and care is a first order public health question [1]. The use of clinical scores such as APACHE II and SAPS in intensive care units (ICU) have improved the follow-up [2]. However, even considering those scores specifically developed for trauma injuries (e.g. TRISS, TRAM, RISC), these exhibit a large space for improvement in calibration and discriminative power [3].

Accordingly, there is a need for biomarker characterization in order to improve precise measurements of severity and pathophysiological consequences, hereby improving their early, scientifically-guided, therapeutic approaches overall leading to a better prognosis [4, 5]. The search for biomarkers has been usually performed within the sepsis context. Though there have been described more than 170 potential biomarkers, still these have not reached the usefulness for prediction of severity and progression of infective processes [6, 7].

In contrast with the development in sepsis, biomarker searching in trauma patients is less advanced, being mostly restricted to traumatic brain injury (TBI). In this context, some of the reported markers show a good correlation with severity, brain damage and global outcome [8]. However, the use of these markers cannot be extended to other trauma types. In search of novel biomarkers to shed light upon pathophysiology venues, one should consider metabolomics. The use of metabolomics techniques allows the study of the complete set of low-molecular-weight intermediates (metabolites), which vary according to the pathologic state of the cell, tissue, organ, or organism and are context-dependent [9]. This approach offers a high number of potential biomarkers when combining the use of high-resolution chromatographic techniques with mass-spectrometry. The use of metabolomics in critical care management could help to a global understanding on metabolism after injury therefore leading to rationally derived therapeutic measures [6, 10].

Globally, despite there have been previous reports on the use of metabolomics in acute trauma care in the intensive care unit, these have been often restricted to the TBI context and/or non-chromatographic techniques (reviewed in [11]). To overcome these limitations, in the present work we have analyzed the metabolomic profiles of severe trauma patients admitted to an ICU, by comparing the early metabolomicprofiles of those patients who survived the trauma versus the ones belonging to non-survivors. The fact that a single metabolite may predict the fatal outcome after multiple injuries can be questioned. For this reason, we also evaluated if the total profile of metabolites recorded within 48 h after the trauma may define the risk of death within 14 days after admission to the intensive care unit. The results demonstrate the potential existence of early metabolomics markers of trauma-associated mortality.


Study design and participants

A prospective pilot study was designed with patients allocated to an ICU (Arnau de Vilanova University Hospital, Lleida, Spain) after a severe trauma. Protocol was supervised and approved by the Institutional Ethics Committee of the Arnau de Vilanova University Hospital. All participants (or their legal representatives) gave their written consent for the study. Inclusion criteria were all trauma-affected patients with severe status requiring ICU admission caused by a traumatic event in the prior 48 h, age > 16 years, with a follow-up period of 14 days in the ICU, in the period between 2011/10/01 and 2012/11/15. All patients whose samples were not obtained in the first 48 h were excluded.. Similarly, we excluded patients with acute kidney failure, defined as plasma creatinine higher than 1.5 times comparing to initial values or diuresis < 400 ml/24 h. Table 1 shows clinical features of included patients. None of the participants received artificial nutritional support in the first 48 h in ICU. Nutritional specificities were chosen on each patient individual basis, mainly according the possibility of enteral pathway and did not influence survival outcome.

Table 1 Demographic and clinical characteristics according to mortality

Variables registered

Age, gender, length of ICU stay (LOS), Injury Severity Score (ISS) [12], Acute Physiology and Chronic Health Evaluation (APACHE II) [13]classification and mortality at 14 days were recorded. We also registered the severity of loss of respiratory function, measured as the lowest PaO2/FiO2 in the first 24 h of ICU admission, by using arterial gasometry, as well as mechanic ventilation (MV) need. Hemodynamic instability (shock) was diagnosed as the patient having a maintained (> 2 h) low systolic blood pressure level (< 90 mmHg) and/or requiring norepinephrine. Transfusion need was defined as those patients requiring more than one hemoconcentrate, and coagulopathy was defined as Quick score < 70% and/or thrombocyte count < 100,000/μL. All scales were validated by 2 independent evaluators.

Sample processing

Blood was obtained in all cases by standard venipuncture or central venous catheterisation, between 24 and 48 h after the traumatic event in sodium citrate tubes. Immediately after sampling, diethylpentaacetic acid (1 mM) and butyl-hydroxy-toluene (10 μM) were added to the sample, in order to avoid artefactual oxidation. Plasma was obtained after centrifugation at 3000 x g, 4 °C for 10 min, and aliquoted into cryovials for immediate storage at − 80 °C. Samples were processed in a double-blind fashion, being aleatorized for extraction and injection. Metabolites from plasma were extracted as previously described [14]: samples were thawed at 4 °C, and 300 μl of cold methanol (containing 1 μg/ml of 13C-phenylalanine as internal standard) were added to 100 μl of plasma for deproteinization, followed by incubation at − 20 °C for 1 h and then, centrifuged at 12000 g for 3 min. The supernatants were recovered, evaporated using a Speed Vac (Thermo Fisher Scientific, Barcelona, Spain) and re-suspended in water 0.4% acetic acid/methanol (50/50).

Metabolomic analyses

For the metabolomic study, an Agilent 1290 liquid chromatography system coupled to an ESI-Q-TOF MS/MS 6520 instrument (Agilent Technologies, Santa Clara, CA, US) was used. In all cases, 2 μL of extracted sample was applied onto a reversed-phase column (Zorbax SB-Aq 1.8 μm 2.1 × 50 mm; Agilent Technologies) equipped with a precolumn (Zorbax-SB-C8 Rapid Resolution Cartridge 2.1 × 30 mm 3.5 μm; Agilent Technologies) with a column temperature of 60 °C. The flow rate was 0.6 mL/min. Solvent A was composed of water containing 0.2% acetic acid and solvent B was composed of methanol 0.2% acetic acid. The gradient started at 2% B and increased to 98% B in 13 min and held at 98% B for 6 min. Post-time was established in 5 min.

Data were collected in positive electrospray mode time of flight operated in full-scan mode at 100–3000 m/z in an extended dynamic range (2 GHz), using N2 as the nebulizer gas (5 L/min, 350 °C). The capillary voltage was 3500 V with a scan rate of 1 scan/s. The ESI source used a separate nebulizer for the continuous, low-level (10 L/min) introduction of reference mass compounds: 121.050873, 922.009798 (positive ion mode) and 119.036320, 966.000725 (negative ion mode), which were used for continuous, online mass calibration. MassHunter Data Analysis Software (Agilent Technologies) was used to collect the results, and MassHunter Qualitative Analysis Software (Agilent Technologies) to obtain the molecular features of the samples, representing different, co-migrating ionic species of a given molecular entity using the Molecular Feature Extractor algorithm (Agilent Technologies), as described [14, 15]. Finally, MassHunter Mass Profiler Professional Software (Agilent Technologies) and Metaboanalyst platform [16] were used to perform a non-targeted metabolomic analysis of the extracted features. We selected samples with a minimum of 2 ions. Multiple charge states were not considered. Compounds from different samples were aligned using a retention time window of 0.1% ± 0.25 min and a mass window of 10.0 ppm ±2.0 mDa. We selected only those features in the 80-100th quartile of relative abundance and corrected for individual bias.

Statistical analyses

The data are presented either as the mean ± standard deviation, the median (interquartile range) or as a percentage. Differences between survivor and non-survivors (NS) were evaluated with Chi-squared (for discrete variables) or with Mann-Whitney test by using the SPSS software ver 24 (IBM Corp, Armonk, NY, USA).

For metabolomic-derived variables, after log transformation and auto scaling of their abundances, multivariate statistics including Hierarchical Clustering, principal component analyses (PCA) and partial least square discriminant analyses (PLS-DA) were done using the Metaboanalyst platform [16]. The same platform was employed to evaluate differences between survivors and non surviving patients were evaluated by Student’s t test (p < 0.05), to generate receiver operating characteristic (ROC) curves, and to explore pathway impact (using hypergeometric test of over representation analyses). In the case of ROC curves, data were not log transformed nor autoscaled, being used the raw values of MS counts.

The identity of differential metabolites were annotated by searching their characteristics in the PCDL database from Agilent (Agilent Technologies, Barcelona, Spain), which uses retention times in a standardized chromatographic system as an orthogonal parameter to complement accurate mass data (accurate mass retention time approach) according to previously published works [16].

Mortality risk was calculated with binary regression logistic (both univariate and multivariate) by calculating odds ratio (OR) and confidence intervals (95%). Those markers with better discriminatory power according ROC curves were used. For the multivariate model we included those variables showing significance in the univariate approach, by using an stepwise selection type.


Table 1 shows clinical features of enclosed patients. APACHE II scores were significantly different between both groups (trauma surviving and non-surviving patients) analyzed. This agrees with previous data, showing that scores higher than 15 points have an overall mortality near 25%, and those scores higher than 25 a mortality of 50% [17] . When evaluating differences between trauma surviving and non-surviving patients there were almost significant differences in age, VM and LOS. All 5 non-surviving patients were victims of a high energy impact (traffic crash casualties or high level fall). In 4 cases, there was generalized hypoperfusion (haemorragic shock and hypoxia) and 1 case it was almost exclusively due to TBI. Of note, one of non-surviving patients had a very differential metabolomic signature (Addittional file 1: Figure S1). As it is shown, this patient presents very different relative abundance of all the metabolites and arouse as a potential outlier. When clinical parameters were check out we discovered that this patient suffered from amyotrophic lateral sclerosis so we discarded the corresponding metabolomics profile for further multivariate analyses, with TBI as predominant mechanism. No patients with sepsis were present in our cohort.

By employing metabolomics we were able to detect 2084 molecular features in plasma (see Declaration section for data availability). As shown in Fig. 1, plasma from non-surviving patients exhibited a specific metabolomic signature in plasma (Fig. 1a) by both PLS-DA and hierarchical clustering analyses. The most important metabolites defining this signature are shown in Fig. 1c. Among them, we were able to identify L-Tyrosine, Uric acid, 2-Amino-3-methyl-1-butanol, hypoxanthine and L-Isoleucine. As shown in Fig. 1b the discrimination between groups using hierarchical clustering algorithm with 25 molecules with the lowest p value (T-Student Test) had a non-optimal performance, with some surviving patients clustered together with non-surviving patients.

Fig. 1
figure 1

Patients non surviving severe trauma exhibit a metabolomic signature. a PLS-DA analyses shows that metabolome is able to discriminate between those patients who survive and those patients who do not survive, with those metabolites contributing to separation shown in the VIP score graph (c). b Heatmap hierarchical clustering analyses using the 25 metabolites with the lowest p value (T-Student test) indicates that there is not a perfect clusterization between groups. PLS-DA cross-validation details (4 components): Accuracy: 0.92, R2: 0.97, Q2: -0.189. The negative value of Q2 means that the model is not all predictive or is overfitted, probably because the low number of not surviving patients

When applying univariate analyses (Student T-Test) 84 significantly different molecules (p < 0.05, 10 upregulated in surviving patients, 74 in non-surviving patients) appeared (Table 2 for annotated molecules and Additional file 1: Table S1 for full details). In order to define the capacity of specific metabolites as potential biomarkers, ROC curves were also performed using those metabolites in 80-100th percentile of abundance. Among the ROC curves with higher area under the curve (AUC), we selected those annotated with a potential identity basing on exact mass and retention time (Fig. 2). All these selected metabolites (3-Indolelactic acid (AUC: 0.938), hydroxyisovaleric acid (AUC: 0.904), phenylethanolamine (AUC: 0.875), cortisol (AUC: 0.865), epinephrine (AUC: 0.865), and myristic acid (AUC: 0.846)) were increased in plasma of non-surviving patients. To further evaluate the capacity of these metabolites we combined their values. Figure 2g shows that combining all these selected metabolites we achieved an AUC = 0.9 and, interestingly, using only the values of cortisol and myristic acid the AUC values reach 0.965. The discrimination capacity was lower than those of APACHE II or ISS (Table 1). Further, no significant correlation was found between these clinical scores and chosen metabolites for further analyses (cortisol, myristic acid), as shown in Table 3. Nonetheless, as shown in Table 4, these metabolomic biomarkers, as well as APACHE II score, were considered risk factors in univariate mortality risk. Noteworthy, myristic acid level was the only independent risk factor in multivariate models of mortality.

Fig. 2
figure 2

ROC curves of those metabolites with a potential identity and higher AUC values (a-f). g ROC curve using the values of all the metabolites from (a) to (f). h ROC curve using Cortisol and Myristic acid

Table 2 Annotated differential metabolites according to mortality
Table 3 Correlation between potential biomarkers of mortality (n = 48)
Table 4 Mortality Risk Model (binary logistic regression)

Pathway analyses reveal the 11 annotated metabolites clustered (p = 0.11) along several amino-acid related pathways (Fig. 3), including fatty acid biosynthesis, and tryptophan and tyrosine metabolism.

Fig. 3
figure 3

Pathway impact of differential metabolites according survival. The 11 differential metabolites annotated (Indolelactic acid, Epinephrine, Phenylethanolamine, 2-Hydroxyisovaleric acid,Cortisol, L-Tryptophan, Myristic acid, Pyridoxal, Bilirubin, Erythrono-1,4-lactone, Elaidic Acid) were searched against KEGG pathway library of H Sapiens using an hypergeometric test for over representation analyses. P value is encoded in intensity of red color of circles (representing pathways) and size of the circles represent the diversity of the pathway, based on the number of metabolites contained in the pathway


The presented results highlight the usefulness of metabolomic approaches in the characterization of pathophysiological events after severe trauma, as suggested by previous reviews [11]. To the best of our knowledge, this is one the first reports using untargeted liquid chromatography coupled to mass spectrometry (LC-MS) in plasma for a non-favorable outcome in a severe trauma population. Previous results, using nuclear magnetic resonance (NMR) have shown that there are clear age and gender-dependent differences in energy metabolism and oxidative stress response in case-control approaches [18]. In a similar case-control study, Brodie et al. showed metabolomics changes in oxidative stress, nucleotide synthesis and muscle catabolism [10], but in this case, they employed targeted LC-MS. However, the study was not addressed for mortality. In a similar study, comprising 95 severely injured patients, targeted metabolomics for lactate and succinate revealed its potential as markers for negative outcome [19]. Similarly, an NMR-based study in the context of combat injuries revealed also succinate as marker of mortality, in addition to malonate [20]. However, the number of metabolites potentially detectable are lower than in our study, for technical limitations. Other studies have focused on the use of plasma metabolomics in TBI management [21, 22], even using cerebrospinal fluid [23] or urine [24] as the sample of reference. Similarly, in a case-control approach, other researchers indicate that a low number of metabolites, one of them stearic acid, allowed for differentiation of TBI from controls, and even helped to classify TBI patients according severity [25].

In our cohort, six annotated metabolites defined survival chances with a very high accuracy. The finding of cortisol agrees with previous reported data [26]. However, these contrast with other reports [27], where vasopressor-dependent critical illness individuals do not show the same behavior. Differences in gender, organ dysfunction and other clinical features could explain this difference. Pathophysiologically, this could be linked to poor glycemic control attributable to high cortisol levels [28]. Interestingly, cortisol is one of the physiological mediators of hyperaminoacidemia [29] and increased branched chain aminoacid levels [30] found in the postinjury phase by previous metabolomic approaches [10]. Of note, epinephrine also predicted low chances of survival. This result agree with previous series, where high catecholamine values were associated to adverse cardiac events [31], or to higher mortality, via endothelial damage and glycocalix degradation, as well as hyperfibrinolysis [32]. The increase of the trace amine phenylethanolamine can be ascribed to the same phenomenon, as it is considered a by-product of catecholamine synthesis [33, 34]. Nonetheless, no previous report on its usefulness as mortality biomarker was present.

3-indolelactic acid is a metabolite of tryptophan degradation [35]. No previous results related this metabolite to non-favorable outcome in trauma. Recent metabolomics reports indicate that this marker is decreased in plasma of patients with cachexia [36]. Nonetheless, previous results with experimental models of cardiac arrest, reveal that tryptophan catabolism was directly related to non-survival outcome [37]. Similarly, high catabolism of tryptophan has been also associated to mortality in bacteremia [38]. Levels of quinolinic acid in cerebrospinal fluid, other of the tryptophan metabolites, are strongly related to low survival chances in TBI. The fact that tryptophan metabolism is increased in trauma patients, even after intravenous aminoacid infusion, is recognized since long time ago [39].

Concerning myristic acid, this long chain fatty acid has marked effect in hepatic cells [40], leading to cellular stress and steatosis. Independent results have revealed that it also confers cardiovascular risk [41]. Myristic acid in blood may be the result of increased lypolisis in adipose tissue, one of the physiological effects of cortisol and epinephrine. However, why only this fatty acid, and not other shorter or longer chain fatty acid correlate with mortality is not known. In cellular terms, myristolation is a key event [42] in immune response, especially in innate immune response, lymphopoiesis for T cells, and the formation of the immunological synapse. Of note, recent evidences link trauma to immune disarrangements, through the release of the so-called damage-associated molecular patterns (DAMPs). DAMPs are considered danger signals which, paradoxically, lead to propagation of injuries to remote tissues, contributing to multiple organ failure and even death (reviewed in [43]). Whether the changes in myristic acid are associated to DAMPs buildup, release or signaling is yet unknown. Other metabolomic works discovered octanoic and decanoic fatty acids, structurally related to myristic, as potential biomarkers of severity in TBI patients [44]. Previous works, employing NMR also uncovered lipids (in their case monounsaturated fatty acids, triacylglycerol and phospholipids) as predictors for survival in the context of trauma [45].

Myristic acid belongs to long chain free fatty acids (FFAs). The relationship of FFAs with a negative outcome are not new, as their concentrations in plasma are directly associated to trauma severity scores [46]. Further, they are increased after bone fracture [47]. Further, non-trauma related intensive-care conditions such as burn or acute pancreatitis injuries also lead to increased concentrations of FFAs [48,49,50]. FFAs could trigger myocardial injury through TLR4 activation [51], and they have been directly related with cardiovascular mortality [52]. In the context of trauma, they inhibit the production of an anti-inflammatory interleukin, IL-10 [53] In addition, they directly induce necrosis [54] through membrane modifications, which could also include increased susceptibility to hemolysis, a well-known condition with negative consequences [55]. Noteworthy, recent data show that differences in albumin FFA-binding properties could explain different capacities of albumin infusions to modulate cell damage in ICU environment [56]. Of note, previous works from our work revealed an impact of FFAs as determinants of lifespan in comparative physiology approaches [57], highlighting the potential influence of circulating FFAs in insulin signaling and other relevant physiological processes [58].

Our study shows some limitations. Regarding the selection of possible biomarkers based on the Student t-test can be questioned, because its use to study differences between groups different in size and variance is prone to errors. This is evident in the box and whisker plots of each of the proposed biomarkers, where some of the highest readings in the surviving group are very similar to the results of the non-survivors. Nonetheless, the ROC curves of the selected biomarkers showed high diagnostic accuracy. PLS-DA model statistics for non-survival could show some overfitting, due to relative low number of samples. In our cohort, mortality was low (ca 10%). In addition, almost no mortality was evident by TBI (which shows a high mortality in the first 96 h after trauma), and many non-surviving patients overcoming the first 48 h could show multiorgan failure, with probably a marked influence in metabolome. Further, those patients with higher severity could require fluid therapy, thereby also conditioning metabolic profile. Similarly, differences in metabolomics profiles due to fasting-fed state before trauma could be also a limitation, as suggested by data in preclinical models [59]. Additionally, age differences, mechanisms of trauma, sex, and body mass index could also impinge changes in metabolome. Of note, this study does not offer a confirmatory cohort, which may be useful for enhancing robustness of findings. Further meta-analyses resulting from exploration with larger datasets/higher mortality rates/highly stratified injury severities are also warranted. However, even with these limitations, similarly sized studies have been reported on the discovery of metabolome-based biomarkers of mortality [20].

We recognize that finding similar candidates within other studies which may show some underpowered sampling does not completely validate our findings. However, these reported studies have been obtained with different analytical techniques, across different geographic environments, thereby suggesting a common pathophysiological foundations. In this respect, the routine use of LC-MS based metabolomics is still not achievable in most hospitals, due to cost-benefit constraints. Nonetheless, LC-MS based metabolomics are helpful for proposal of novel potential biomarkers, enhancing patient’s stratification. Further, its application could pave the way for discovery of pathophysiologically relevant pathways, hereby leading to rationally designed and personalized therapeutics which should enhance the prognosis of these patients.


Overall, revealed metabolic pathways disclose a catecholamine and cortisol stress response with a potential role of fatty acid metabolism and degradation of specific amino acids. Further, the fact that multivariate and ROC analyses show a better behavior than consolidated clinical markers could be helpful in the search for measures enhancing prognosis in these patients.



Acute Physiology and Chronic Health Evaluation


Area Under the Curve


Damage-Associated Molecular Patterns


Electrospray Ionization


Free Fatty Acids


Fraction of Inspired Oxygen


Intensive Care Unit


Injury Severity Score


Kyoto Encyclopedia of Genes and Genomes


Liquid Chromatography Coupled to Mass Spectrometry


Length of ICU Stay




Mechanic Ventilation


Nuclear Magnetic Resonance


Odds Ratio


Arterial Partial Pressure of Oxygen


Principal Component Analyses


Partial Least Square Discriminant Analyses (PLS-DA)


Quadrupole Time-Of-Flight


Revised Injury Severity Classification


Receiver Operating Curve


Simplified Acute Physiologic Score


Traumatic Brain Injury


Trauma Risk Adjustment Model


Trauma and Injury Severity Score


  1. Norton R, Kobusingye O. Injuries. N Engl J Med. 2013;368:1723–30.

    Article  CAS  PubMed  Google Scholar 

  2. Vincent J-L, Moreno R. Clinical °review: scoring systems in the critically ill. Crit Care. 2010;14:207.

    Article  PubMed  PubMed Central  Google Scholar 

  3. Serviá L, Badia M, Montserrat N, Trujillano J. Severity scores in trauma patients admitted to ICU. Physiological and anatomic models. Med Intensiva. 2019;43(1):26–34.

    Article  PubMed  Google Scholar 

  4. Legrand M, Januzzi JL, Mebazaa A. Critical research on biomarkers: what’s new? Intensive Care Med. 2013;39:1824–8.

    Article  PubMed  Google Scholar 

  5. Garcia-Simon M. Prognosis biomarkers of severe sepsis and septic shock by 1H NMR urine metabolomics in the intensive care unit. PLoS One. 2015;10(11):e0140993. eCollection 2015.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  6. Serkova NJ, Standiford TJ, Stringer KA. The emerging field of quantitative blood metabolomics for biomarker discovery in critical illnesses. Am J Respir Crit Care Med. 2011;184:647–55.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  7. Pierrakos C, Vincent J-L. Sepsis biomarkers: a review. Crit Care. 2010;14:R15.

    Article  PubMed  PubMed Central  Google Scholar 

  8. Gordillo-Escobar E, Egea-Guerrero JJ, Rodríguez-Rodríguez A, Murillo-Cabezas F. Usefulness of biomarkers in the prognosis of severe head injuries. Med Int. 2016;40:105–12.

    Article  CAS  Google Scholar 

  9. Jové M, Portero-Otín M, Naudí A, Ferrer I, Pamplona R. Metabolomics of human brain aging and age-related neurodegenerative diseases. J Neuropathol Exp Neurol. 2014;73:640–57.

    Article  CAS  PubMed  Google Scholar 

  10. Parent BA, Seaton M, Sood RF, Gu H, Djukovic D, Raftery D, et al. Use of metabolomics to trend recovery and therapy after injury in critically ill trauma patients. JAMA Surg. 2016;151:e160853.

    Article  PubMed  Google Scholar 

  11. Jayaraman SP, Anand RJ, DeAntonio JH, Mangino M, Aboutanos MB, Kasirajan V, et al. Metabolomics and precision medicine in trauma: the state of the field. Shock. 2017.

    Article  CAS  Google Scholar 

  12. Baker SP, O’Neill B, Haddon W, Long WB. The injury severity score: a method for describing patients with multiple injuries and evaluating emergency care. J Trauma. 1974;14:187–96.

    Article  CAS  Google Scholar 

  13. Knaus WA, Wagner DP, Draper EA, Zimmerman JE, Bergner M, Bastos PG, et al. The APACHE III prognostic system. Risk prediction of hospital mortality for critically ill hospitalized adults. Chest. 1991;100:1619–36.

    Article  CAS  Google Scholar 

  14. Jové M, Mauri-Capdevila G, Suárez I, Cambray S, Sanahuja J, Quílez A, et al. Metabolomics predicts stroke recurrence after transient ischemic attack. Neurology. 2015;84:36–45.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  15. Sana TR, Roark JC, Li X, Waddell K, Fischer SM. Molecular formula and METLIN personal metabolite database matching applied to the identification of compounds generated by LC/TOF-MS. J Biomol Tech. 2008;19:258–66.

    PubMed  PubMed Central  Google Scholar 

  16. Xia J, Sinelnikov IV, Han B, Wishart DS. MetaboAnalyst 3.0--making metabolomics more meaningful. Nucleic Acids Res. 2015;43:W251–7.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  17. Knaus WA, Draper EA, Wagner DP, Zimmerman JE. APACHE II: a severity of disease classification system. Crit Care Med. 1985;13:818–29.

    Article  CAS  Google Scholar 

  18. Lusczek ER, Myers C, Popovsky K, Mulier K, Beilman G, Sawyer R. Plasma metabolomics pilot study suggests age and sex-based differences in the metabolic response to traumatic injury. Injury. 2018.

    Article  Google Scholar 

  19. D’Alessandro A, Moore HB, Moore EE, Reisz JA, Wither MJ, Ghasabyan A, et al. Plasma succinate is a predictor of mortality in critically injured patients. J Trauma Acute Care Surg. 2017.

    Article  Google Scholar 

  20. Lusczek ER, Muratore SL, Dubick MA, Beilman GJ. Assessment of key plasma metabolites in combat casualties. J Trauma Acute Care Surg. 2017;82:309–16.

    Article  CAS  PubMed  Google Scholar 

  21. Dickens AM, Posti JP, Takala RS, Ala-Seppälä HM, Mattila I, Coles JC, et al. Serum metabolites associate with CT findings following TBI. J Neurotrauma. 2018;35:2673–83.

    Article  PubMed  Google Scholar 

  22. Wolahan SM, Lebby E, Mao HC, McArthur D, Real C, Vespa PM, et al. Novel metabolomic comparison of arterial and jugular venous blood in severe adult TBI patients and the impact of pentobarbital infusion. J Neurotrauma. 2018.

    Article  Google Scholar 

  23. Hagos FT, Empey PE, Wang P, Ma X, Poloyac SM, Bayır H, et al. Exploratory application of neuropharmacometabolomics in severe childhood traumatic brain injury. Crit Care Med. 2018.

    Article  CAS  Google Scholar 

  24. Ottens AK, Stafflinger JE, Griffin HE, Kunz RD, Cifu DX, Niemeier JP. Post-acute brain injury urinary signature: a new resource for molecular diagnostics. J Neurotrauma. 2014;31:782–8.

    Article  PubMed  PubMed Central  Google Scholar 

  25. Fiandaca MS, Mapstone M, Mahmoodi A, Gross T, Macciardi F, Cheema AK, et al. Plasma metabolomic biomarkers accurately classify acute mild traumatic brain injury from controls. PLoS One. 2018;13:e0195318.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  26. Pandya U, Polite N, Wood T, Lieber M. Increased total serum random cortisol levels predict mortality in critically ill trauma patients. Am Surg. 2014;80:1112–8.

    PubMed  Google Scholar 

  27. Rady MY, Johnson DJ, Patel B, Larson J, Helmers R. Cortisol levels and corticosteroid administration fail to predict mortality in critical illness: the confounding effects of organ dysfunction and sex. Arch Surg. 2005;140:661–8; discussion 669.

    Article  PubMed  Google Scholar 

  28. Gale SC, Sicoutris C, Reilly PM, Schwab CW, Gracias VH. Poor glycemic control is associated with increased mortality in critically ill trauma patients. Am Surg. 2007;73:454–60.

    PubMed  Google Scholar 

  29. Gelfand RA, Matthews DE, Bier DM, Sherwin RS. Role of counterregulatory hormones in the catabolic response to stress. J Clin Invest. 1984;74:2238–48.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  30. Mao H, Wang H, Wang B, Liu X, Gao H, Xu M, et al. Systemic metabolic changes of traumatic critically ill patients revealed by an NMR-based metabonomic approach. J Proteome Res. 2009;8:5423–30.

    Article  CAS  PubMed  Google Scholar 

  31. Naganathar S, De’Ath HD, Wall J, Brohi K. Admission biomarkers of trauma-induced secondary cardiac injury predict adverse cardiac events and are associated with plasma catecholamine levels. J Trauma Acute Care Surg. 2015;79:71–7.

    Article  CAS  PubMed  Google Scholar 

  32. Johansson PI, Stensballe J, Rasmussen LS, Ostrowski SR. High circulating adrenaline levels at admission predict increased mortality after trauma. J Trauma Acute Care Surg. 2012;72:428–36.

    Article  CAS  Google Scholar 

  33. Ammar NM, Farag MA, Kholeif TE, Metwally NS, El-Sheikh NM, El Gendy AN, et al. Serum metabolomics reveals the mechanistic role of functional foods and exercise for obesity management in rats. J Pharm Biomed Anal. 2017;142:91–101.

    Article  CAS  PubMed  Google Scholar 

  34. Shannon HE, Cone EJ, Yousefnejad D. Physiologic effects and plasma kinetics of phenylethanolamine and its N-methyl homolog in the dog. J Pharmacol Exp Ther. 1981;217:379–85.

    CAS  PubMed  Google Scholar 

  35. Morita I, Kawamoto M, Hattori M, Eguchi K, Sekiba K, Yoshida H. Determination of tryptophan and its metabolites in human plasma and serum by high-performance liquid chromatography with automated sample clean-up system. J Chromatogr. 1990;526:367–74.

    Article  CAS  Google Scholar 

  36. Cala MP, Agulló-Ortuño MT, Prieto-García E, González-Riano C, Parrilla-Rubio L, Barbas C, et al. Multiplatform plasma fingerprinting in cancer cachexia: a pilot observational and translational study. J Cachexia Sarcopenia Muscle. 2018;9:348–57.

    Article  PubMed  PubMed Central  Google Scholar 

  37. Ristagno G, Fries M, Brunelli L, Fumagalli F, Bagnati R, Russo I, et al. Early kynurenine pathway activation following cardiac arrest in rats, pigs, and humans. Resuscitation. 2013;84:1604–10.

    Article  CAS  PubMed  Google Scholar 

  38. Huttunen R, Syrjänen J, Aittoniemi J, Oja SS, Raitala A, Laine J, et al. High activity of indoleamine 2,3 dioxygenase enzyme predicts disease severity and case fatality in bacteremic patients. Shock. 2010;33:149–54.

    Article  CAS  PubMed  Google Scholar 

  39. Moyer ED, McMenamy RH, Cerra FB, Reed RA, Yu L, Chenier R, et al. Multiple systems organ failure: III contrasts in plasma amino acid profiles in septic trauma patients who subsequently survive and do not survive-effects of intravenous amino acids. J Trauma. 1981;21:263–74.

    Article  CAS  Google Scholar 

  40. Speziali G, Liesinger L, Gindlhuber J, Leopold C, Pucher B, Brandi J, et al. Myristic acid induces proteomic and secretomic changes associated with steatosis, cytoskeleton remodeling, endoplasmic reticulum stress, protein turnover and exosome release in HepG2 cells. J Proteome. 2018;181:118–30.

    Article  CAS  Google Scholar 

  41. Ebbesson SOE, Voruganti VS, Higgins PB, Fabsitz RR, Ebbesson LO, Laston S, et al. Fatty acids linked to cardiovascular mortality are associated with risk factors. Int J Circumpolar Health. 2015;74:28055.

    Article  PubMed  Google Scholar 

  42. Udenwobele DI, Su R-C, Good SV, Ball TB, Varma Shrivastav S, Shrivastav A. Myristoylation: an important protein modification in the immune response. Front Immunol. 2017;8:751.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  43. Bortolotti P, Faure E, Kipnis E. Inflammasomes in tissue damages and immune disorders after trauma. Front Immunol. 2018;9:1900.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  44. Posti JP, Dickens AM, Orešič M, Hyötyläinen T, Tenovuo O. Metabolomics profiling as a diagnostic tool in severe traumatic brain injury. Front Neurol. 2017;8:398.

    Article  PubMed  PubMed Central  Google Scholar 

  45. Cohen MJ, Serkova NJ, Wiener-Kronish J, Pittet J-F, Niemann CU. 1H-NMR-based metabolic signatures of clinical outcomes in trauma patients--beyond lactate and base deficit. J Trauma. 2010;69:31–40.

    Article  PubMed  Google Scholar 

  46. Stoner HB, Frayn KN, Barton RN, Threlfall CJ, Little RA. The relationships between plasma substrates and hormones and the severity of injury in 277 recently injured patients. Clin Sci. 1979;56:563–73.

    Article  CAS  Google Scholar 

  47. Nixon JR, Brock-Utne JG. Free fatty acid and arterial oxygen changes following major injury: a correlation between hypoxemia and increased free fatty acid levels. J Trauma. 1978;18:23–6.

    Article  CAS  Google Scholar 

  48. Kamolz LP, Andel H, Mittlböck M, Winter W, Haslik W, Meissl G, et al. Serum cholesterol and triglycerides: potential role in mortality prediction. Burns. 2003;29:810–5.

    Article  Google Scholar 

  49. Sztefko K, Panek J. Serum free fatty acid concentration in patients with acute pancreatitis. Pancreatology. 2001;1:230–6.

    Article  CAS  PubMed  Google Scholar 

  50. Qi P, Abdullahi A, Stanojcic M, Patsouris D, Jeschke MG. Lipidomic analysis enables prediction of clinical outcomes in burn patients. Sci Rep. 2016;6:38707.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  51. Wang Y, Qian Y, Fang Q, Zhong P, Li W, Wang L, et al. Saturated palmitic acid induces myocardial inflammatory injuries through direct binding to TLR4 accessory protein MD2. Nat Commun. 2017;8:13997.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  52. Miedema MD, Maziarz M, Biggs ML, Zieman SJ, Kizer JR, Ix JH, et al. Plasma-free fatty acids, fatty acid-binding protein 4, and mortality in older adults (from the cardiovascular health study). Am J Cardiol. 2014;114:843–8.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  53. Ikeoka DT, Pachler C, Mader JK, Bock G, Neves AL, Svehlikova E, et al. Lipid-heparin infusion suppresses the IL-10 response to trauma in subcutaneous adipose tissue in humans. Obesity (Silver Spring). 2011;19:715–21.

    Article  CAS  Google Scholar 

  54. Qin X, Dong W, Sharpe SM, Sheth SU, Palange DC, Rider T, et al. Role of lipase-generated free fatty acids in converting mesenteric lymph from a noncytotoxic to a cytotoxic fluid. Am J Physiol Gastrointest Liver Physiol. 2012;303:G969–78.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  55. Hardaway RM, Vasquez Y. A shock toxin that produces disseminated intravascular coagulation and multiple organ failure. Am J Med Sci. 2001;322:222–8.

    Article  CAS  Google Scholar 

  56. Penn AH, Dubick MA, Torres Filho IP. Fatty acid saturation of albumin used in resuscitation fluids modulates cell damage in shock: in vitro results using a novel technique to measure fatty acid binding capacity. Shock. 2017.

    Article  CAS  Google Scholar 

  57. Jové M, Naudí A, Aledo JC, Cabré R, Ayala V, Portero-Otin M, et al. Plasma long-chain free fatty acids predict mammalian longevity. Sci Rep. 2013;3:3346.

    Article  PubMed  PubMed Central  Google Scholar 

  58. Karpe F, Dickmann JR, Frayn KN. Fatty acids, obesity, and insulin resistance: time for a reevaluation. Diabetes. 2011;60:2441–9.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  59. Witowski N, Lusczek E, Determan C, Lexcen D, Mulier K, Ostrowski B, et al. A four-compartment metabolomics analysis of the liver, muscle, serum, and urine response to polytrauma with hemorrhagic shock following carbohydrate prefeed. PLoS One. 2015;10:e0124467.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

Download references


The authors are in debt with patients and families for allowing this research. We are also indebted to ICU nurses for its help in sample obtention.


This work has been partially supported by the IRBLleida biobank and RETICS BIOBANCOS RD09/0076/00059, the Spanish Ministry of Economy and Competitiveness, Institute of Health Carlos III (PI 17–00134) and the Generalitat of Catalonia (2017SGR969). FEDER Funds are also acknowledged: “A way to make Europe”. MJ is a professor under the Serra Hunter program (Generalitat de Catalunya).

Availability of data and materials

The datasets generated and/or analysed during the current study are available in the repository,

Author information

Authors and Affiliations



LS and JT designed the experimental work, drafted the manuscript and analyzed data. MJ, JS and RP performed and analysed metabolomics analyses. NM, MB and LS obtainded samples. JT and MPO performed global analyses, planned general structure, prepared figures and wrote and revised the whole manuscript. All authors read and approved the final manuscript.

Corresponding authors

Correspondence to Manuel Portero-Otin or Javier Trujillano.

Ethics declarations

Ethics approval and consent to participate

Protocol was supervised and approved by the Institutional Ethics Committee of the Arnau de Vilanova University Hospital. All participants (or their legal representatives) gave their written consent for the study.

Consent for publication

Not applicable.

Competing interests

The authors declare that they have no competing interests.

Publisher’s Note

Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Additional file

Additional file 1:

Table S1. Differential metabolites according mortality (full list). Figure S1. Clinical specificity for metabolomics signature. A. PLS-DA analyses shows that metabolome is able to discriminate between those patients who survive and those patients who do not survive. Although the metabolomic profile of survived patients is homogenous, one of the non-surviving patients had a specific metabolomic profile (black arrow). B. Heatmap hierarchical clustering analyses using the 25 metabolites with the lowest p value (T-Student test) confirms that one patient has a specific metabolomic profile (black arrow). PLS-DA cross-validation details (4 components): Accuracy: 0.91, R2: 0.96, Q2: -0.141. The negative value of Q2 means that the model is not all predictive or is overfitted, probably because the low number of not surviving patients. (DOCX 239 kb)

Rights and permissions

Open Access This article is distributed under the terms of the Creative Commons Attribution 4.0 International License (, which permits unrestricted use, distribution, and reproduction in any medium, provided you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. The Creative Commons Public Domain Dedication waiver ( applies to the data made available in this article, unless otherwise stated.

Reprints and permissions

About this article

Check for updates. Verify currency and authenticity via CrossMark

Cite this article

Servià, L., Jové, M., Sol, J. et al. A prospective pilot study using metabolomics discloses specific fatty acid, catecholamine and tryptophan metabolic pathways as possible predictors for a negative outcome after severe trauma. Scand J Trauma Resusc Emerg Med 27, 56 (2019).

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: