Skip to main content

Metabolomic diferences between COVID-19 and H1N1 influenza induced ARDS



Acute respiratory distress syndrome (ARDS) is a type of respiratory failure characterized by lung inflammation and pulmonary edema. Coronavirus disease 2019 (COVID-19) is associated with ARDS in the more severe cases. This study aimed to compare the specificity of the metabolic alterations induced by COVID-19 or Influenza A pneumonia (IAP) in ARDS.


Eighteen patients with ARDS due to COVID-19 and twenty patients with ARDS due to IAP, admitted to the intensive care unit. ARDS was defined as in the American-European Consensus Conference. As compared with patients with COVID-19, patients with IAP were younger and received more often noradrenaline to maintain a mean arterial pressure > 65 mm Hg. Serum samples were analyzed by Nuclear Magnetic Resonance Spectroscopy. Multivariate Statistical Analyses were used to identify metabolic differences between groups. Metabolic pathway analysis was performed to identify the most relevant pathways involved in ARDS development.


ARDS due to COVID-19 or to IAP induces a different regulation of amino acids metabolism, lipid metabolism, glycolysis, and anaplerotic metabolism. COVID‐19 causes a significant energy supply deficit that induces supplementary energy-generating pathways. In contrast, IAP patients suffer more marked inflammatory and oxidative stress responses. The classificatory model discriminated against the cause of pneumonia with a success rate of 100%.


Our findings support the concept that ARDS is associated with a characteristic metabolomic profile that may discriminate patients with ARDS of different etiologies, being a potential biomarker for the diagnosis, prognosis, and management of this condition.

Graphical Abstract


The current Coronavirus disease 2019 (COVID-19) pandemic, caused by the Severe Acute Respiratory Syndrome Coronavirus 2 (SARS-CoV-2) [1], strained critical care resources in many countries, and the management of lung injury in these patients posed a tremendous challenge for clinicians [2]. Several patients with COVID-19 developed severe acute respiratory distress syndrome (ARDS) with mortality rates around 30% [3] in the first pandemic wave. ARDS is characterized by lung inflammation and hyperpermeability pulmonary edema. Currently, the diagnosis of ARDS is based on the presence of clinical, physiological, and radiological criteria [4,5,6]. Unlike other clinical conditions, to date, there are no specific molecular markers that help in the prognosis of this condition. Advancement in the understanding of the pathogenesis of ARDS is necessary for designing innovative and effective therapeutic approaches.

Molecular approaches are needed to understand the mechanisms of ARDS induced by COVID-19. Of particular interest is the use of metabolomics for the characterization of this condition. The metabolome reflects early and specific alterations in the pathophysiological state of biological systems. In this context, Magnetic Resonance Spectroscopy (MRS) emerges as a highly potential tool for studying metabolic disorders in respiratory diseases [7]. Several studies have proved the potential of MRS-based metabolomics to monitor patients with ARDS induced by respiratory infections [8,9,10]. A previous study compared the specific metabolic fingerprint of ARDS patients with either influenza A pneumonia (IAP) caused by the H1N1-2009 influenza virus or pneumonia caused by Streptococcus pneumoniae [11]. Here, we used a similar approach for the characterization of the metabolic fingerprint of COVID-19-induced ARDS. We compared COVID-19 and IAP patients to identify the metabolic reprogramming involved in these two conditions. The identification of metabolic pathways involved in ARDS caused by the H1N1-2009 influenza virus or by SARS-CoV-2 will improve our understanding of the pathogenesis of COVID-19. Finally, as a proof of concept of the diagnostic potential of these metabolic biomarkers, we developed a predictive model to identify the etiological pathogens responsible for ARDS.


Experimental design

Patients with Coronavirus disease 2019 (COVID-19, n = 18) were recruited from March 1, 2020, to June 31, 2020, in Hospital Universitario de Getafe, Madrid. Metabolic profile from ARDS patients with H1N1-2009 influenza pneumonia was acquired and analyzed in our previous study [12]. Serum samples from ARDS patients with H1N1-2009 influenza virus pneumonia were obtained in Hospital Universitario de Getafe, Madrid, Spain (n = 10) and Hospital del Mar, Barcelona, Spain (n = 10) during the 2009 pandemic and stored at − 80 °C until NMR analysis in September of 2017. H1N1-2009 infection was confirmed by RT-PCR or either nasopharyngeal swab samples or tracheal secretions. In all cases, serum samples were obtained within 24 h of presentation to the emergency department. Blood samples were collected in BD Vacutainer tubes after each participant signed informed consent. After collection, the sample was left at room temperature for 30 min to clot. The clot was removed by centrifugation at 1500 × g for 10 min at 4 °C. The resulting supernatant was immediately transferred into 2 ml Eppendorf tubes and stored at − 80 °C.

Inclusion criteria for both studies were: age ≥ 18 years, diagnosis of ARDS, confirmed infection by SARS-CoV-2 or H1N1-2009 influenza virus by real-time reverse transcription-polymerase chain reaction (RT-PCR) of nasopharyngeal swab samples, and admission to the Intensive Care Unit (ICU). ARDS was defined as in the American-European Consensus Conference (AECC) [13]. All patients were mechanically ventilated.

Clinical information was obtained by retrospective chart review, and data of the Sequential Organ Failure Assessment (SOFA) and the Simplified Acute Physiologic Score-II (SAPS II) scores on admission, the presence of renal or cardiovascular failure (SOFA score of the respective component > 2) [14] and status at hospital discharge (hospital mortality) were collected.

NMR data acquisition

Serum samples were collected after each participant signed informed consent within 24 h of ICU admission and examined (40 µl of serum) by high-resolution magic angle spinning (HR-MAS) NMR operating at 4 °C to reduce metabolic degradation. HR-MAS NMR was performed at 500.13 MHz using a Bruker AMX500 spectrometer 11.7 T. HR-MAS NMR has several strengths for clinical studies [15]; (1) signals in NMR spectrum have the same sensitivity independently of the properties of the metabolite; (2) a combination of NMR techniques enables the unambiguous identification of metabolic signals; (3) HR-MAS NMR enables analysis of intact samples, which is essential considering that factors associated with sample preparation contribute to analytical variability; (4) metabolite profiles obtained by NMR are virtually independent of the operator and instrument used, which provides a high degree of reliability to the derived results.

Samples were placed into a 50 μl zirconium oxide rotor using a rinsed with a cylindrical insert, together with 15 µl of 0.1 mM solution of TSP in deuterium water (D2O), and spun at 4000 Hz spinning rate to remove the effects of spinning sidebands from the acquired spectra. Several bidimensional homonuclear and heteronuclear experiments such as standard gradient-enhanced correlation spectroscopy (COSY), 1H–1H total correlated spectroscopy (TOCSY), and gradient-selected heteronuclear single quantum correlation (HSQC) protocols were performed to carry out component assignments. Between consecutive two-dimensional (2D) spectra, a control 1H NMR spectrum was continuously measured to detect metabolic degradation or microbiologically contamination. No metabolic differences were noted in the signals of multiple spectra acquired under the same conditions. Standard solvent-suppressed spectra were grouped into 32,000 data points, averaged over 256 acquisitions [16]. The data acquisition lasted 13 min using a sequence based on the first increment of the nuclear Overhauser effect spectroscopy (NOESY) pulse sequence to effect suppression of the water. Sample acquisitions were performed using a spectral width of 8333.33 Hz before Fourier transformation, and the free induction decay (FID) signals were multiplied by an exponential weight function corresponding to a line broadening of 0.3 Hz. Spectra were referenced to the TSP singlet at 0 ppm chemical shift.

NMR data were processed for statistical analysis and metabolic quantification. The chemical shifts region from 5.00 to 5.20 ppm was excluded from the analysis to remove the random and known effects of variation in the water resonance suppression. Similarly, the chemical shifts region from 0 to 0.04 ppm containing the internal reference (TSP) was excluded from the statistical analysis. Phase and baseline correction (Whittaker method) were performed automatically using MestRenova v. 8.1 software (Mestrelab Research S.L., Santiago de Compostela, Spain). For statistical analysis, 1H NMR spectra were automatically data-reduced to integral segments or buckets of equal length (δ = 0.01 ppm) to compensate for variations in resonance positions [17], and they were normalized to the total sum of the spectral regions. For metabolic quantification, full resolution spectra were normalized to the total sum of the spectral regions too. Relative intensity was calculated as the initial intensity normalized to the total sum of the spectral regions.

Statistical analysis

Quantitative and qualitative variables were compared by the Student´s t-test or the Chi-square test, respectively. A p value less than 0.05 was considered statistically significant. The statistical package SPSS IBM Statistics 19.0 was used for the analysis. Descriptive data are presented as mean (SD) for continuous variables and percentages for discrete variables.

Principal Components Analysis (PCA) [18] was performed over binned NMR spectral data using the Metabonomic package (rel.3.3.1) [19] to analyze in ARDS patients the differences between SARS-CoV-2 and H1N1-2009 influenza infection. In PCA, the data collected on a set of samples are resolved into principal components. The first principal component is defined by the spectral profile (loadings) in the data that describes most of the variation; the second principal component, orthogonal to the first, is the second-best profile describing the variation, and so on. The principal components are composed of so-called scores and loadings. Loadings contain information about the variables (chemical shifts) in the dataset, and scores hold information on samples (intensities) in the dataset. Before PCA, NMR processed data were centered, and Pareto scaled [20]. Hotteling’s T2 test [21] identified the NMR areas (NMR bins) from the PCA loading matrix responsible for group clustering. The NMR signals in the identified NMR areas were individually integrated into full resolution spectra for metabolic quantification using the Global Spectral Deconvolution algorithm of MestRenova v. 8.1 (Mestrelab Research S.L., Santiago de Compostela, Spain). Metabolites identification was performed manually using Chenomx Profiler tool [22]. Metabolites assignments were confirmed by analyzing 2D-NMR spectra using MestRenova software and the Human Metabolome Database [23].

For the metabolic quantification, statistical significance was determined using a Bonferroni corrected Student’s t-test [24], assuming unequal variance with p < 0.05 considered significant.

Partial Least Square Discriminant Analysis (PLS-DA) [25] was developed as a classificatory model using MetaboAnalyst v.5.0 [26] to differentiate in ARDS patients those due to IAP from those due to COVID-19. PLS-DA models have commonly used classification methods for analyzing high-dimensional data. The number of latent variables used to develop the PLS-DA model were evaluated by R2 and Q2 robustness parameters. R2 can be considered a metric of how the algorithm fits the training data, and Q2 is a metric of algorithm performance on test data [27]. Q2 parameter, which evaluates the classification functions derived from the probability of belonging to each group, was computed by leave-one-out cross-validation (LOOCV) to minimize the variance in training [28]. Three PLS components were selected to develop a classificatory model based on the best robustness results (R2 = 0.94; Q2 = 0.89). Model performance was evaluated by Prediction accuracy during cross-validation and the Area Under de Curve Receiver Operating Characteristic (AUC-ROC). AUC-ROC curve is a performance measurement for classification problems at various threshold settings. The ROC curve is plotted with True Positive Rate (sensitivity) against the False Positive Rate (1-specificity).

Metabolic pathways analyses

Metabolic pathway analysis was performed using The Pathway Analysis module [29] of Metaboanalyst v.5.0 [26] that combines results from robust pathway enrichment analysis [30] with pathway topology analysis [31] to help researchers identify the most relevant pathways involved in the conditions under study.

Briefly, pathway enrichment analysis examines whether metabolites in predefined pathways are at the top or bottom of a ranked list. In contrast, pathway topology analysis applies graph theory to measure the importance of an experimentally identified metabolite in a predefined metabolic pathway. KEGG metabolic pathways were used as the backend knowledgebase, the selected pathway enrichment analysis method was GlobalAncova [32], node importance measure for topological analysis was out-degree centrality. Centrality is a standard metric used in graph theory to estimate the relative importance of individual nodes to the overall network [33]. Out-degree is the number of outgoing links or the number of successor nodes.


Characteristics of study patients and laboratory findings on admission

We compared ARDS patients with COVID-19 (n = 18) with ARDS patients with IAP (n = 20). Compared to patients with COVID-19, patients with IAP were younger and received more often noradrenaline to maintain a mean arterial pressure > 65 mm Hg (Table 1). Disease severity as measured by SAPS II score, SOFA score ad mortality was similar in both groups.

Table 1 Clinical characteristics of the study groups

Metabolomic analysis

An unsupervised classification study with PCA was carried out to analyze in ARDS patients the differences between SARS-CoV-2 and H1N1-2009 influenza infection. PCA (Fig. 1A) provided nearly perfect discrimination between the two groups of subjects. The resonances identified as significantly different by PCA loadings analysis (Fig. 1B) were individually integrated for metabolic quantification. The resonances were identified according to Chenomx Profiler database (Fig. 2) and characteristic cross-peaks from 2D spectra to help in the unequivocal assignation of these metabolites. ARDS patients induced by COVID-19 showed higher metabolic concentrations of free fatty acids, acetone, creatinine, and lactate, and lower metabolic concentrations of valine, 2-hydroxybutyrate, proline, methyl-guanidine, glucose, and tyrosine (Table 2; Additional file 1: Figure S1).

Fig. 1

Principal component analysis (PCA) performed on the NMR data of serum samples from ARDS patients diagnosed with Coronavirus disease 2019 (COVID-19) and H1N1-2009 influenza pneumonia (IAP). A Score plot discriminated between groups along PC2. B PCA 2 loading plot identified the resonances that induce the clustering between IAP and COVID-19 groups. The figure highlights the significant metabolites

Fig. 2

Representative 1H HR-MAS NMR spectrum serum sample from COVID-19 ARDS patient. The figure highlights metabolic assignments

Table 2 Relative change in the concentration of the identified metabolites

The specific metabolic fingerprint (unassigned 1H NMR spectra) was then used to develop a partial least squares discriminant analysis (PLS-DA) to identify SARS-CoV-2 infection. PLS-DA successfully discriminated COVID-19 from IAP (Fig. 3) (Prediction accuracy during cross validation = 100%; AUCROC = 1).

Fig. 3

A Predicted class probabilities plot (average of the cross-validation) for each sample using the best classifier (based on R2 = 0.94 and Q2 = 0.89 PLS robustness parameters) of the PLS model (3 PLS components). As a balanced subsampling approach is used for model training, the classification boundary is always at the center (x = 0.5, the dotted line). B Area under de curve receiver operating characteristic (AUCROC) of PLS-DA. COVID, COVID-19; H1N1, influenza A pneumonia

Metabolic pathway analysis was performed to identify the most relevant pathways involved in ARDS development (Fig. 4; Additional file 1: Table S1). This pathway analysis identified alterations in amino acids biosynthesis and metabolism, glycerolipid metabolism and fatty acid degradation, glycolysis, and anaplerotic metabolism.

Fig. 4

Summary of metabolic pathway analysis. Y axis represents the statistical p values from enrichment analysis, and the X axis represents the pathway impact value calculated from pathway topology analysis


The description of the metabolic alterations induced by SARS-CoV-2 infection in ICU patients is fundamental for a better understanding of the pathobiology of the disease. In the present study, we compared the metabolomic profile of ARDS due either to IAP or to COVID-19 by MRS using untargeted multivariate statistical analysis and metabolic pathway analysis. We found that the activation of many metabolic pathways was different between ARDS patients with COVID-19 or IAP. Furthermore, the serum metabolite profile of patients with ARDS discriminates the specific virus infection (H1N1-2009 influenza pneumonia versus SARS-CoV-2 pneumonia). PLS-DA model provided a classification accuracy of 100%. These findings are helpful for the understanding of the pathogenesis of severe COVID-19. Specifically, the metabolomic profile of ARDS in these patients suggests alterations in energy pathways, inflammatory response, and oxidative stress.

Previous studies that have analyzed the metabolism of patients with COVID-19 [34,35,36,37,38,39] were designed to compare the metabolic profile of COVID-19 patients with healthy controls or to evaluate the metabolic differences between patients with a positive or negative outcome. Thus they cannot discern between the metabolic dysregulation due to SARS-CoV-2 infection or due to ARDS development. To the best of our knowledge, this is the first study designed to compare the metabolic profile of ICU patients with similar severity of ARDS due to COVID-19 or to other viral respiratory infections, e.g. IAP.

We have found specific metabolic differences between ARDS patients induced by COVID-19 or IAP. Most of these metabolic alterations have been previously reported as biomarkers of ARDS or ARDS severity [8,9,10]. For example, a similar serum metabolic profile including proline, glutamate, phenylalanine, and valine was reported as a sensitive biomarker of ARDS severity (mild, moderate, and severe) [40]. However, we have to consider that, unlike previous studies, in which metabolic changes induced by the viral infection itself or by the occurrence of ARDS cannot be distinguished, IAP and COVID-19 patients in the present study all met the criteria for the diagnosis of ARDS, and disease severity was similar. Thus differences in the metabolic profile herein reported are better explained by virus-specific pathogenetic mechanisms rather than by the occurrence of ARDS or by disease severity.

The ability of patients to normalize energy metabolism has been reported as one of the critical factors determining the COVID-19 progression [39]. Compared with IAP patients, COVID-19 patients showed up-regulation of energy-generating pathways, i.e. glycolysis, fatty acid degradation, CoA biosynthesis, glycerolipids, and glycerophospholipids metabolism. The increase of lactate-to-glucose ratio found in COVID-19 patients is a biomarker of the up-regulation in the glycolysis pathway [41]. In the same context, dysregulation of the choline metabolism and elevated levels of free poly-unsaturated fatty acids are biomarkers of the energy deficiency reported in COVID-19 patients [42]. However, other authors [34, 43] have associated the dysregulation of lipid metabolism with a higher atherosclerotic risk in COVID-19 patients. In the same line, COVID-19 patients showed a higher phenylalanine-tyrosine ratio that has been associated with an adverse outcome [39] and may indicate a higher cardiovascular risk [44] in COVID-19 patients [45]. Other supplementary energy-generating pathways were also up-regulated. The excess of ketone bodies such as acetone suggests that they are used as an alternative energy source. Ketosis could be explained in the context of acute illness and lack of adequate caloric intake.

Alteration in amino acidic metabolism has been reported as one of the key features of ARDS development [8, 9], and it has also been found significantly up-regulated in COVID-19 patients [34, 39, 46, 47]. However, when we compared the metabolic profile of patients with ARDS due to COVID-19 or to IAP, we found that the amino acidic metabolism was decreased in COVID-19 patients. The serum concentrations of branched-chain amino acids (BCAAs), including isoleucine and valine, were decreased in COVID-19 compared with IAP patients. As elevated circulating BCAAs may promote oxidative stress [48], the lower levels of BCAAs in patients with COVID-19 may result in less intense inflammatory response as compared to patients with influenza A [49]. Downregulated BCAAs may also be considered a potential marker of the infection and its further involvement in the dysregulation of pantothenate and CoA biosynthesis [47], as confirmed by the enrichment analysis. Pantothenic acid (vitamin B5) is required for coenzyme A formation and is also essential for α-ketoglutarate and pyruvate dehydrogenase complexes as well as fatty acid oxidation, compromising the mitochondrial energy metabolism [50]. The increase in 2-hydroxybutyric acid, a readout of hepatic glutathione synthesis and marker of oxidative stress [50], and essential amino acids such as proline [51] confirmed more marked inflammatory and oxidative stress responses [52] in IAP than in COVID-19 patients. A previous study also identified dysregulation of propanoate metabolism as a novel pathway in the progression of COVID-19 [46], suggesting potential roles played by gut microbiota in the immune response [53]. Finally, the methyl-guanidine-to-creatinine ratio is an index of hydroxyl radical formation in the lung, and it was identified previously as a specific metabolic pattern of IAP [11].

Several limitations of the study should be acknowledged. First, the metabolic concentrations reported are relative to the total metabolic concentration, and baseline clinical differences among groups should be taken into consideration when interpreting the results. We did not perform absolute quantification because of limitations in sample manipulation for HR-MAS NMR analysis. External validation is required before the application of the specific metabolic fingerprint in clinical practice. Second, despite the overall similarity in disease severity (as measured by the SAPS II score, the SOFA score, and the mortality), patients differed in some characteristics, such as age and the requirement of noradrenaline, that could have an impact on the metabolic profile. Also, oxygenation impairment differed in the two groups, although the difference did not reach statistical significance. Third, some aspects of patients management could differ in the cohorts, as they span some years apart. Specifically, the way patients were mechanically ventilated could have had an impact on the metabolic profile. After the publication of the ARMA trial in 2000 [54], different studies reported changes in the way mechanical ventilation was used, i.e., a lower tidal volume and slightly higher PEEP levels [55, 56]. However, other studies have failed to show significant changes after 2010 [57, 58]. Thus, it is unlikely that the metabolic changes reported in the present study are due to different mechanical ventilation strategies in the two groups. Fourth, ARDS was diagnosed according to the AECC, followed when the first cohort was recruited.


In summary, we have characterized a specific metabolomic fingerprint that allows the discrimination between ARDS due to SARS-CoV-2 or to H1N1-2009 influenza virus in ARDS patients. The description of the metabolic alterations herein reported will help better understand the pathobiology of ARDS and its different causes and may have translational implications for biomarker discovery and the design of novel therapeutic targets.

Availability of data and materials

The datasets generated and/or analysed during the current study will be available in a online repository after acceptance.


  1. 1.

    Albarello F, Pianura E, Di Stefano F, Cristofaro M, Petrone A, Marchioni L, et al. 2019-novel coronavirus severe adult respiratory distress syndrome in two cases in Italy: an uncommon radiological presentation. Int J Infect Dis. 2020;93:192–7.

    CAS  PubMed  PubMed Central  Google Scholar 

  2. 2.

    White DB, Lo B. A framework for rationing ventilators and critical care beds during the COVID-19 pandemic. JAMA. 2020;323(18):1773–4.

    CAS  PubMed  Google Scholar 

  3. 3.

    Penuelas O, Del Campo-Albendea L, de Aledo ALG, Anon JM, Rodriguez-Solis C, Mancebo J, et al. Long-term survival of mechanically ventilated patients with severe COVID-19: an observational cohort study. Ann Intensive Care. 2021;11(1):143.

    CAS  PubMed  PubMed Central  Google Scholar 

  4. 4.

    Force ADT, Ranieri VM, Rubenfeld GD, Thompson BT, Ferguson ND, Caldwell E, et al. Acute respiratory distress syndrome: the Berlin Definition. JAMA. 2012;307(23):2526–33.

    Google Scholar 

  5. 5.

    Lorente JA, Cardinal-Fernandez P, Munoz D, Frutos-Vivar F, Thille AW, Jaramillo C, et al. Acute respiratory distress syndrome in patients with and without diffuse alveolar damage: an autopsy study. Intensive Care Med. 2015;41(11):1921–30.

    PubMed  Google Scholar 

  6. 6.

    Thille AW, Esteban A, Fernandez-Segoviano P, Rodriguez JM, Aramburu JA, Penuelas O, et al. Comparison of the Berlin definition for acute respiratory distress syndrome with autopsy. Am J Respir Crit Care Med. 2013;187(7):761–7.

    PubMed  Google Scholar 

  7. 7.

    Snowden S, Dahlén S-E, Wheelock CE. Application of metabolomics approaches to the study of respiratory diseases. Bioanalysis. 2012;4(18):2265–90.

    CAS  PubMed  Google Scholar 

  8. 8.

    Metwaly SM, Winston BW. Systems biology ARDS research with a focus on metabolomics. Metabolites. 2020;10(5):207.

    CAS  PubMed Central  Google Scholar 

  9. 9.

    Viswan A, Singh C, Kayastha AM, Azim A, Sinha N. An NMR based panorama of the heterogeneous biology of acute respiratory distress syndrome (ARDS) from the standpoint of metabolic biomarkers. NMR Biomed. 2020;33(2):e4192.

    PubMed  Google Scholar 

  10. 10.

    Stringer KA, McKay RT, Karnovsky A, Quémerais B, Lacy P. Metabolomics and its application to acute lung diseases. Front Immunol. 2016;7:44.

    PubMed  PubMed Central  Google Scholar 

  11. 11.

    Izquierdo-García JL, Nin N, Cardinal-Fernandez P, Ruiz-Cabello J. Lorente JÁJIcm: metabolomic profile of acute respiratory distress syndrome of different etiologies. Intensive Care Med. 2019;45(9):1318–20.

    PubMed  PubMed Central  Google Scholar 

  12. 12.

    Izquierdo-Garcia JL, Nin N, Jimenez-Clemente J, Horcajada JP, Arenas-Miras MDM, Gea J, et al. Metabolomic profile of ards by nuclear magnetic resonance spectroscopy in patients with H1N1 influenza virus pneumonia. Shock. 2017;50:504.

    Google Scholar 

  13. 13.

    Bernard GR, Artigas A, Brigham KL, Carlet J, Falke K, Hudson L et al. The American-European Consensus Conference on ARDS: Definitions, mechanisms, relevant outcomes, and clinical trial coordination. Am J Respir Critical Care Med. 1994; 149:818–824.

  14. 14.

    Vincent JL, Moreno R, Takala J, Willatts S, De Mendonca A, Bruining H, et al. The SOFA (Sepsis-related Organ Failure Assessment) score to describe organ dysfunction/failure: On behalf of the Working Group on Sepsis-Related Problems of the European Society of Intensive Care Medicine. Intensive Care Med. 1996;22(7):707–10.

    CAS  PubMed  Google Scholar 

  15. 15.

    Gowda GN, Raftery D. Can NMR solve some significant challenges in metabolomics? J Magn Reson. 2015;260:144–60.

    PubMed Central  Google Scholar 

  16. 16.

    Yoon D, Kim YJ, Lee WK, Choi BR, Oh SM, Lee YS, et al. Metabolic changes in serum metabolome of beagle dogs fed black ginseng. Metabolites. 2020;10(12):517.

    CAS  PubMed Central  Google Scholar 

  17. 17.

    Holmes E, Foxall P, Nicholson J, Neild G, Brown S, Beddell C, et al. Automatic data reduction and pattern recognition methods for analysis of 1H nuclear magnetic resonance spectra of human urine from normal and pathological states. Anal Biochem. 1994;220(2):284–96.

    CAS  PubMed  Google Scholar 

  18. 18.

    Hotelling H. Analysis of a complex of statistical variables into principal components. J Educ Psychol. 1933;24(6):417.

    Google Scholar 

  19. 19.

    Izquierdo-Garcia JL, Rodriguez I, Kyriazis A, Villa P, Barreiro P, Desco M, et al. A novel R-package graphic user interface for the analysis of metabonomic profiles. BMC Bioinform. 2009;10:363.

    Google Scholar 

  20. 20.

    van den Berg RA, Hoefsloot HC, Westerhuis JA, Smilde AK, van der Werf MJ. Centering, scaling, and transformations: improving the biological information content of metabolomics data. BMC Genomics. 2006;7(1):1–15.

    Google Scholar 

  21. 21.

    Hotelling H. The generalization of student’s ratio. Ann Math Stat. 1931;2(3):360–78.

    Google Scholar 

  22. 22.

    Weljie AM, Newton J, Mercier P, Carlson E, Slupsky CM. Targeted profiling: quantitative analysis of 1H NMR metabolomics data. Anal Chem. 2006;78(13):4430–42.

    CAS  PubMed  Google Scholar 

  23. 23.

    Wishart DS, Jewison T, Guo AC, Wilson M, Knox C, Liu Y, et al. HMDB 3.0—the human metabolome database in 2013. Nucleic Acids Res. 2013;41:D801-807.

    CAS  PubMed  Google Scholar 

  24. 24.

    Vinaixa M, Samino S, Saez I, Duran J, Guinovart JJ, Yanes O. A guideline to univariate statistical analysis for LC/MS-based untargeted metabolomics-derived data. Metabolites. 2012;2(4):775–95.

    CAS  PubMed  PubMed Central  Google Scholar 

  25. 25.

    Lindon JC, Holmes E, Nicholson JK. Pattern recognition methods and applications in biomedical magnetic resonance. Prog Nucl Magn Reson Spectrosc. 2001;39(1):1–40.

    CAS  Google Scholar 

  26. 26.

    Chong J, Wishart DS, Xia J. Using MetaboAnalyst 4.0 for comprehensive and integrative metabolomics data analysis. Curr Protocols Bioinform. 2019;68(1):e86.

    Google Scholar 

  27. 27.

    Eriksson L, Byrne T, Johansson E, Trygg J, Vikström C: Multi-and megavariate data analysis basic principles and applications, vol. 1. Umetrics Academy; 2013.

  28. 28.

    Lachenbruch PA, Mickey MR. Estimation of error rates in discriminant analysis. Technometrics. 1968;10(1):1–11.

    Google Scholar 

  29. 29.

    Xia J, Wishart DS. MetPA: a web-based metabolomics tool for pathway analysis and visualization. Bioinformatics. 2010;26(18):2342–4.

    CAS  PubMed  Google Scholar 

  30. 30.

    Kankainen M, Gopalacharyulu P, Holm L, Orešič M. MPEA—metabolite pathway enrichment analysis. Bioinformatics. 2011;27(13):1878–9.

    CAS  PubMed  Google Scholar 

  31. 31.

    Aittokallio T, Schwikowski B. Graph-based methods for analysing networks in cell biology. Brief Bioinform. 2006;7(3):243–55.

    CAS  PubMed  Google Scholar 

  32. 32.

    Hummel M, Meister R, Mansmann U. GlobalANCOVA: exploration and assessment of gene group effects. Bioinformatics. 2008;24(1):78–85.

    CAS  PubMed  Google Scholar 

  33. 33.

    Liu G, Lee DP, Schmidt E, Prasad GL. Pathway analysis of global metabolomic profiles identified enrichment of caffeine, energy, and arginine metabolism in smokers but not moist snuff consumers. Bioinform Biol Insights. 2019;13:1177932219882961–1177932219882961.

    PubMed  PubMed Central  Google Scholar 

  34. 34.

    Bruzzone C, Bizkarguenaga M, Gil-Redondo R, Diercks T, Arana E, GarcíadeVicuña A, et al. SARS-CoV-2 infection dysregulates the metabolomic and lipidomic profiles of serum. Science. 2020;23(10):101645.

    CAS  Google Scholar 

  35. 35.

    Shen B, Yi X, Sun Y, Bi X, Du J, Zhang C, et al. Proteomic and metabolomic characterization of COVID-19 patient sera. Cell. 2020;182(1):59-72.e15.

    CAS  PubMed  PubMed Central  Google Scholar 

  36. 36.

    Wu D, Shu T, Yang X, Song J-X, Zhang M, Yao C, et al. Plasma metabolomic and lipidomic alterations associated with COVID-19. Natl Sci Rev. 2020;7(7):1157–68.

    CAS  PubMed  PubMed Central  Google Scholar 

  37. 37.

    Lodge S, Nitschke P, Kimhofer T, Coudert JD, Begum S, Bong S-H, et al. NMR spectroscopic windows on the systemic effects of SARS-CoV-2 infection on plasma lipoproteins and metabolites in relation to circulating cytokines. J Proteome Res. 2021;20(2):1382–96.

    CAS  PubMed  Google Scholar 

  38. 38.

    Lodge S, Nitschke P, Kimhofer T, Wist J, Bong S-H, Loo RL, et al. Diffusion and relaxation edited proton nmr spectroscopy of plasma reveals a high-fidelity supramolecular biomarker signature of SARS-CoV-2 infection. Anal Chem. 2021;93(8):3976–86.

    CAS  PubMed  PubMed Central  Google Scholar 

  39. 39.

    Baranovicova E, Bobcakova A, Vysehradsky R, Dankova Z, Halasova E, Nosal V, et al. The ability to normalise energy metabolism in advanced COVID-19 disease seems to be one of the key factors determining the disease progression—a metabolomic NMR study on blood plasma. Appl Sci. 2021;11(9):4231.

    CAS  Google Scholar 

  40. 40.

    Viswan A, Ghosh P, Gupta D, Azim A, Sinha N. Distinct metabolic endotype mirroring acute respiratory distress syndrome (ARDS) subphenotype and its heterogeneous biology. Sci Rep. 2019;9(1):1–13.

    CAS  Google Scholar 

  41. 41.

    Serkova NJ, Van Rheen Z, Tobias M, Pitzer JE, Wilkinson JE, Stringer KA. Utility of magnetic resonance imaging and nuclear magnetic resonance-based metabolomics for quantification of inflammatory lung injury. Am J Physiol Lung Cel Mol Physiol. 2008;295(1):L152–61.

    CAS  Google Scholar 

  42. 42.

    Schwarz B, Sharma L, Roberts L, Peng X, Bermejo S, Leighton I, et al. Cutting edge: Severe SARS-CoV-2 infection in humans is defined by a shift in the serum lipidome, resulting in dysregulation of eicosanoid immune mediators. J Immunol. 2021;206(2):329–34.

    CAS  PubMed  Google Scholar 

  43. 43.

    Darwesh AM, Bassiouni W, Sosnowski DK, Seubert JM. Can N-3 polyunsaturated fatty acids be considered a potential adjuvant therapy for COVID-19-associated cardiovascular complications? Pharmacol Ther. 2021;219:107703.

    CAS  PubMed  Google Scholar 

  44. 44.

    Murr C, Grammer TB, Meinitzer A, Kleber ME, März W, Fuchs D. Immune activation and inflammation in patients with cardiovascular disease are associated with higher phenylalanine to tyrosine ratios: the ludwigshafen risk and cardiovascular health study. J Amino Acids. 2014;2014:783730.

    PubMed  PubMed Central  Google Scholar 

  45. 45.

    Meoni G, Ghini V, Maggi L, Vignoli A, Mazzoni A, Salvati L, et al. Metabolomic/lipidomic profiling of COVID-19 and individual response to tocilizumab. PLoS Pathog. 2021;17(2):e1009243.

    CAS  PubMed  PubMed Central  Google Scholar 

  46. 46.

    Pang Z, Zhou G, Chong J, Xia J. Comprehensive meta-analysis of COVID-19 global metabolomics datasets. Metabolites. 2021;11(1):44.

    CAS  PubMed  PubMed Central  Google Scholar 

  47. 47.

    Barberis E, Timo S, Amede E, Vanella VV, Puricelli C, Cappellano G, et al. Large-scale plasma analysis revealed new mechanisms and molecules associated with the host response to SARS-CoV-2. Int J Mol Sci. 2020;21(22):8623.

    CAS  PubMed Central  Google Scholar 

  48. 48.

    Zhenyukh O, Civantos E, Ruiz-Ortega M, Sánchez MS, Vázquez C, Peiró C, et al. High concentration of branched-chain amino acids promotes oxidative stress, inflammation and migration of human peripheral blood mononuclear cells via mTORC1 activation. Free Radical Biol Med. 2017;104:165–77.

    CAS  Google Scholar 

  49. 49.

    Mudd PA, Crawford JC, Turner JS, Souquette A, Reynolds D, Bender D, et al. Distinct inflammatory profiles distinguish COVID-19 from influenza with limited contributions from cytokine storm. Sci Adv. 2020;6(50):eabe3024.

    CAS  PubMed  PubMed Central  Google Scholar 

  50. 50.

    Depeint F, Bruce WR, Shangari N, Mehta R, O’Brien PJ. Mitochondrial function and toxicity: role of the B vitamin family on mitochondrial energy metabolism. Chem Biol Interact. 2006;163(1):94–112.

    CAS  PubMed  Google Scholar 

  51. 51.

    Becker DF, Zhu W, Natarajan S. Proline metabolism and protection against oxidative stress. FASEB J. 2009;23:678.677.

    Google Scholar 

  52. 52.

    Phang JM, Pandhare J, Liu Y. The metabolism of proline as microenvironmental stress substrate. J Nutr. 2008;138(10):2008S-2015S.

    CAS  PubMed  Google Scholar 

  53. 53.

    Dhar D, Mohanty A. Gut microbiota and Covid-19- possible link and implications. Virus Res. 2020;285:198018.

    CAS  PubMed  PubMed Central  Google Scholar 

  54. 54.

    Network ARDS. Ventilation with lower tidal volumes as compared with traditional tidal volumes for acute lung injury and the acute respiratory distress syndrome. N Engl J Med. 2000;342(18):1301–8.

    Google Scholar 

  55. 55.

    Esteban A, Ferguson ND, Meade MO, Frutos-Vivar F, Apezteguia C, Brochard L, et al. Evolution of mechanical ventilation in response to clinical research. Am J Respir Crit Care Med. 2008;177(2):170–7.

    PubMed  Google Scholar 

  56. 56.

    Esteban A, Frutos-Vivar F, Muriel A, Ferguson ND, Peñuelas O, Abraira V, et al. Evolution of mortality over time in patients receiving mechanical ventilation. Am J Respir Crit Care Med. 2013;188(2):220–30.

    PubMed  Google Scholar 

  57. 57.

    Peñuelas O, Muriel A, Abraira V, Frutos-Vivar F, Mancebo J, Raymondos K, et al. Inter-country variability over time in the mortality of mechanically ventilated patients. Intensive Care Med. 2020;46(3):444–53.

    PubMed  PubMed Central  Google Scholar 

  58. 58.

    Peñuelas O, Frutos-Vivar F, Muriel A, Mancebo J, García-Jiménez A, de Pablo R, et al. Mechanical ventilation in Spain, 1998–2016: epidemiology and outcomes. Medicina Intensiva (English Edition). 2021;45(1):3–13.

    Google Scholar 

Download references


This research was supported by grants: (i) PID2019-10656RJ-I00 from the Spanish Ministry of Science and Innovation; (ii) S2017/BMD-3727-EXOHEP-CM from Comunidad de Madrid and Fondos FEDER, Madrid, Spain; (iii) B2017/BMD3875 from the Comunidad de Madrid, Madrid, Spain; (iv) SAF2017-84494-C2-1-R from the Spanish Ministry of Economy, Industry, and Competitiveness (MEIC-AEI); (v) Fondo Sectorial de Salud (155219), Agencia Nacional de Investigación e Innovación (vi) KK-2019/bmG19 from the Gobierno Vasco, Dpto. Industria, Innovación, Comercio y Turismo, under the ELKARTEK program and (vii) the European Union’s Horizon 2020 Research and Innovation Program, under the Marie Skłodowska-Curie grant agreement no. 823854 (INNOVA4TB). JRC received grants from the BBVA Foundation (Ayudas a Equipos de Investigación Científica de Biomedicina 2018) and from La Caixa Foundation (Health Research Call 2020 / HR20-00075). CIC biomaGUNE is supported by the Maria de Maeztu Units of Excellence Programme from the Spanish State Research Agency (Grant No. MDM-2017-0720). The funders were not involved at any stage, from study design to submission of the manuscript for publication. NMR studies were performed at the CAI BIOIMAC (Centro de Asistencia a la Investigación Bioimagen Complutense) node of the ICTS (Infraestructura Científicas y Técnicas Singulares) REDiB. The authors would also like to thank the nursing staff of the ICU and the medical team for their collaboration in managing these patients under unusually stressful conditions.


This research was supported by grants: (i) PID2019-10656RJ-I00 from the Spanish Ministry of Science and Innovation; (ii) S2017/BMD-3727-EXOHEP-CM from Comunidad de Madrid and Fondos FEDER, Madrid, Spain; (iii) B2017/BMD3875 from the Comunidad de Madrid, Madrid, Spain; (iv) SAF2017-84494-C2-1-R from the Spanish Ministry of Economy, Industry, and Competitiveness (MEIC-AEI); (v) Fondo Sectorial de Salud (155219), Agencia Nacional de Investigación e Innovación (vi) KK-2019/bmG19 from the Gobierno Vasco, Dpto. Industria, Innovación, Comercio y Turismo, under the ELKARTEK program and (vii) the European Union’s Horizon 2020 Research and Innovation Program, under the Marie Skłodowska-Curie grant agreement no. 823854 (INNOVA4TB). JRC received grants from the BBVA Foundation (Ayudas a Equipos de Investigación Científica de Biomedicina 2018) and from La Caixa Foundation (Health Research Call 2020/HR20-00075). CIC biomaGUNE is supported by the Maria de Maeztu Units of Excellence Programme from the Spanish State Research Agency (Grant No. MDM-2017-0720). The funders were not involved at any stage, from study design to submission of the manuscript for publication.

Author information




Conceptualization: JAL, NN, JRC, and JLIG; methodology: JAL, PV, and JLIG; software: JAL, NN and JLIG; formal analysis: JAL, PV, ABMC and JLIG; investigation: DV, RH, and OP; resources: JAL, NN, JRC and JLIG; data curation: JAL, and JLIG; writing—original draft preparation: JLIG; writing—review and editing: JAL, NN, JRC and JLIG; supervision: JAL, NN, JRC and JLIG; funding acquisition: JAL, NN, JRC, and JLIG. All authors have read and agreed to the published version of the manuscript. All authors read and approved the final manuscript.

Corresponding author

Correspondence to Jose L. Izquierdo-Garcia.

Ethics declarations

Ethics approval and consent to participate

The study was conducted according to the guidelines of the Declaration of Helsinki, and approved by the Ethics Committee (CEIm) of Hospital Universitario de Getafe (A03/20. April 28, 2020).

Consent for publication

Not applicable.

Competing interests

The authors declare no competing interests.

Additional information

Publisher's Note

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

Supplementary Information

Additional file 1.

Table S1: Metabolic Pathway analysis results. Figure S1: Boxplot figures of the metabolites found significant different (p-value < 0.05) between H1N1 influenza and COVID-19 patients.

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 The Creative Commons Public Domain Dedication waiver ( applies to the data made available in this article, unless otherwise stated in a credit line to the data.

Reprints and Permissions

About this article

Verify currency and authenticity via CrossMark

Cite this article

Lorente, J.A., Nin, N., Villa, P. et al. Metabolomic diferences between COVID-19 and H1N1 influenza induced ARDS. Crit Care 25, 390 (2021).

Download citation


  • ARDS
  • COVID-19
  • H1N1 influenza
  • Metabolomics
  • NMR
  • Diagnosis