(C) PLOS One
This story was originally published by PLOS One and is unaltered.
. . . . . . . . . .
COVID-19: A complex disease with a unique metabolic signature [1]
['Veronica Ghini', 'Department Of Chemistry', 'Ugo Schiff', 'University Of Florence', 'Sesto Fiorentino Florence', 'Magnetic Resonance Center', 'Cerm', 'Sesto Fiorentino', 'Florence', 'Walter Vieri']
Date: 2023-11
Plasma of COVID-19 patients contains a strong metabolomic/lipoproteomic signature, revealed by the NMR analysis of a cohort of >500 patients sampled during various waves of COVID-19 infection, corresponding to the spread of different variants, and having different vaccination status. This composite signature highlights common traits of the SARS-CoV-2 infection. The most dysregulated molecules display concentration trends that scale with disease severity and might serve as prognostic markers for fatal events. Metabolomics evidence is then used as input data for a sex-specific multi-organ metabolic model. This reconstruction provides a comprehensive view of the impact of COVID-19 on the entire human metabolism. The human (male and female) metabolic network is strongly impacted by the disease to an extent dictated by its severity. A marked metabolic reprogramming at the level of many organs indicates an increase in the generic energetic demand of the organism following infection. Sex-specific modulation of immune response is also suggested.
Metabolites and lipoproteins are the main components of human plasma and their concentration can be determined by nuclear magnetic resonance. In COVID-19 patients there are significant alterations in the concentration of several of these molecules. Using the plasma of more than 600 subjects (510 patients sampled in the acute phase of the infection plus 95 independent recovered subjects), we demonstrate that the dysregulation of these molecules is a function of the disease severity but it is not affected by either the SARS-CoV-2 variants or the vaccination status. The disease signature is particularly evident in those cases that subsequently evolve towards a fatal outcome, and could have prognostic value. Building on this large amount of data, we propose a metabolic reconstruction of the disease using a sex-specific multi-organ metabolic model.
Funding: This work was funded by Regione Toscana, project COMETA, Bando COVID-19 (
https://www.regione.toscana.it/-/bando-ricerca-covid-19-toscana ) to GL and PT. The project also covers the salaries of VG, WV, VP, TC and NB. The funders had no role in study design, data collection and analysis, decision to publish, or preparation of the manuscript.
Data Availability: This study is available at the NIH Common Fund's National Metabolomics Data Repository (NMDR) website, the Metabolomics Workbench,
https://www.metabolomicsworkbench.org where it has been assigned Study ID ST002404. The data can be accessed directly via its Project DOI:
http://dx.doi.org/10.21228/M89T2Q .
Copyright: © 2023 Ghini et al. This is an open access article distributed under the terms of the Creative Commons Attribution License , which permits unrestricted use, distribution, and reproduction in any medium, provided the original author and source are credited.
Here, we report a detailed and comprehensive characterization of the metabolomic and lipoproteomic fingerprint of plasma samples of > 500 hospitalized COVID-19 patients, with different disease severities, infected with different viral variants and with different vaccination status. Sex-specific differences as well as the contribution of several comorbidities are also analysed. Our data deeply extend a first metabolomic/lipoproteomic characterization of the disease published at the beginning of 2022 [ 18 ], performed on a smaller number of subjects of the same cohort, infected before a significant spread of the δ variant and before widespread COVID-19 vaccination. While confirming previously identified severity markers [ 18 – 28 ], we establish for the first time a correlation between the levels of a few metabolites and lipoproteins and the fatal outcome of the disease. These molecules can therefore be proposed as predictive and prognostic biomarkers. The observed changes are interpreted through simulations of the overall metabolic state of the human body with a recently developed sex-specific multi-organ metabolic model [ 30 ], which until now has been used to predict known biomarkers of inherited metabolic diseases in different biofluids. Based on this reconstruction and using as input the changes in the metabolome observed in our cohort, we obtained a comprehensive view of the impact of COVID-19 on the entire human metabolism, which represents a step forward with respect to previous models based on not cell type-specific and not sex-specific flux balance analyses [ 31 ].
Massive worldwide efforts by research groups using omics sciences have been made to unravel the disease mechanisms and to identify biomarkers of the disease severity [ 9 – 13 ]. In this framework, 1 H NMR [ 14 – 17 ] plays a role for its ability to reveal a complex blood plasma signature exhibiting the presence of a strong fingerprint of the COVID-19 disease. Highly reproducible alterations of a large number of blood metabolites and lipoprotein parameters were identified as markers of the disease, suggesting the reprogramming of important metabolic pathways aimed at the energy supply for viral replication and for host immunological response [ 18 – 29 ]. However, how the administration of vaccines and the development of the different viral variants impact on the disease fingerprint, remain poorly characterized aspects.
From the year 2020, the resilience of worldwide national health systems was profoundly challenged by the Coronavirus disease 2019 (COVID-19) pandemic caused by the severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2). The viral infection is characterized by a broad spectrum of clinical manifestations from an asymptomatic or pauci-symptomatic disease (in more than 80% of subjects), to interstitial pneumonia and acute respiratory distress syndrome requiring hospitalization and even ventilation of the patients [ 1 – 3 ]. Moreover, it has been demonstrated that a significant subset of patients developed concomitant multi-organ dysfunctions, including acute kidney and liver injuries, thromboembolism and sepsis that contribute to a fatal outcome [ 4 , 5 ]. Additionally, some comorbidities have been proposed as risk factors of disease severity and fatal outcome [ 6 , 7 ]. Although vaccines are now available and have demonstrated high efficacy in decreasing the severity of SARS-CoV-2, vaccination does not prevent SARS-CoV-2 transmission. The diffusion of the different viral variants as well as the phenomenon of the so called “Long-COVID”, i.e. long-term effects associated to COVID-19 infection [ 8 ], further complicate this scenario.
Results
In this study, EDTA-plasma samples from 510 COVID-19 positive patients (COVID-19 group) were collected during the period 20/06/2020-17/06/2022. The demographic and clinical data of the cohort are reported in S1 Table. This population well represents the incidence of the disease in Tuscan hospitals along the course of the pandemic, in terms of sex, age, severity, and main risk factors/comorbidities. One hundred and fifty-four of these individuals were also evaluated at a follow-up visit; this group is hereafter named as the “follow-up group”, S2 Table. The reference group COVID-19-R (S1 Table) is constituted by 95 recovered subjects who had contracted the infection during the first wave and did not show any symptom of long COVID at the follow-up visit, when they were sampled. Consistently their metabolite and lipoprotein levels fall within the normality ranges [19].
Fingerprint of the main COVID-19 variants To evaluate whether the metabolic and lipoprotein fingerprint of the disease was significantly changed as a function of the different variants of the virus, we selected 3 groups of samples. Based on the data of the statistical incidence of the various SARs-Cov-2 variants in the Tuscany Region (TR) (TR does not systematically determine the variants of each patient through genomic sequencing), only patients infected in periods characterized by defined dominant variants were included. In particular, we considered: i) 87 samples collected when the Wild-Type (WT), α, and β variants were prevalent (01/2020-06/2021), the “wt-α-β group”–these samples are new with respect to those included in our previous study [18]; ii) 91 samples collected when the δ variant was dominant, i.e. the “δ group” (07/2021-12/2021); iii) 44 samples collected when the ο variant was dominant (01/2022-02/2022), the “ο group”. Together, the three groups account for a total of 222 samples. The detailed demographic and clinical characteristics of these subjects are reported in Fig 1A. PPT PowerPoint slide
PNG larger image
TIFF original image Download: Fig 1. Metabolomic profiling of the main COVID-19 variants. (A) Main demographic and clinical characteristics of the enrolled subjects (right panel); for each group of subjects, the sample collection period was reported (left panel). (B) List of the metabolites quantified in plasma samples. The p-values and Cliff’s Delta effects size are reported for the comparison between each of the three groups of variants with respect to the COVID-19-R group; p-values <0.05 are highlighted. (C) Values of Log 2 fold change (FC) of quantified metabolites. Positive/negative values have higher/lower concentration in plasma samples from each group of variants with respect to the COVID-19-R group. Colour coding: wt-α-β group (cyan); δ group (red); o group (orange).
https://doi.org/10.1371/journal.ppat.1011787.g001 From the multivariate statistics (see Methods) no clear differences emerge among the metabolomic fingerprints of the three groups of variants; instead, all of them can be almost perfectly discriminated from COVID-19 recovered subjects (S1 Fig). The three groups show very similar trends for both metabolite and lipoprotein levels, Figs 1B–1C and 2. Conversely, eleven out of 25 metabolites and 16 out of 30 lipoprotein main parameters and main fractions result to be significantly different from the recovered subjects, regardless of the variant, Figs 1B–1C and 2. PPT PowerPoint slide
PNG larger image
TIFF original image Download: Fig 2. Lipoproteomic profiling of the main COVID-19 variants. A) List of lipoprotein parameters (main parameters, calculated features and main fractions) quantified in plasma samples. The p-values and Cliff’s Delta effects size are reported for the comparison between each of the three groups of variants with respect to the COVID-19-R group; p-values <0.05 are highlighted. B) Values of Log 2 fold change (FC) of quantified lipoprotein parameters (main parameters, calculated features and main fractions). Positive/negative values have higher/lower concentration in plasma samples from each group of variants with respect to COVID-19-R group. Colour coding: wt-α-β group (cyan); δ group (red); o group (orange).
https://doi.org/10.1371/journal.ppat.1011787.g002 The three ketone bodies (3-hydroxybutyrate, acetoacetic acid and acetone) along with the amino acids Phe, Met and Ile, the sugars glucose and mannose, and the glycoproteins GlycA and GlycB, are significantly increased in all COVID-19 groups. Citric acid and acetic acid, instead, show decreased levels in all COVID-19 groups. Histidine is the only molecule that shows a different direction of changes among the three groups with respect to the COVID-19-R group; it decreased in the wt-α-β group while it increased in the δ and ο groups. For a few other metabolites (Val, Tyr, lactic acid, creatine, Leu, Gln, formic acid, Ala) we observe the same trends of changes among the three groups, although their variations are not consistently significant. The lack of significance could result from different number of samples in the various groups combined with relatively small fold changes (|Log 2 (FC)| < 1). As for lipoproteins, all the parameters, except free cholesterol (Chol) and apolipoprotein B100 (Apo B100) associated to IDL, show the same trends along all the variants, Fig 2. Significantly decreased level of LDL- and HDL-Chol, total Apo A1 and Apo A2, Apo A1-HDL, Apo A2-HDL, triglycerides (TG)-HDL and phospholipids (PhL)-LDL, -IDL and -HDL were observed in all the COVID-19 groups with respect to the recovered subjects. Increased values of the ratio Apo B100/Apo A1 and of TG-LDL are also coherently observed. Moreover, significantly increased levels of free Chol-LDL are observed only in wt-α-β and δ groups, whereas increased levels of the Apo B100-VLDL are observed only in the wt-α-β and o variants. A significant decrement of ApoB-LDL is instead only reported for the δ group. Lipoprotein sub-fractions and particle numbers are reported in S2 Fig.
Fingerprint of the disease as a function of the vaccination status Among the 510 COVID-19 positive hospitalized subjects, 71 were vaccinated (with two or more doses of DNA or RNA vaccine, administered at least 15 days before the first positive swab). With the aim of characterizing if the vaccine introduces some important alterations in the metabolomic/lipoproteomic fingerprint of the disease, we compared these samples (hereafter named “VAX group”) with the samples of 80 COVID-19 positive patients unvaccinated collected in the same period (“NO-VAX group”). The demographic and clinical data of these groups are reported in Fig 3A. PPT PowerPoint slide
PNG larger image
TIFF original image Download: Fig 3. Metabolomic and lipoproteomic profiling of the VAX and NO-VAX groups. (A) Main demographic and clinical characteristics of the subjects. (B-C) Values of Log 2 fold change (FC) of quantified metabolites (B) and lipoprotein parameters (main parameters, calculated features and main fractions) (C). Positive/negative values have higher/lower concentration in plasma samples from the VAX or NO-VAX groups with respect to COVID-19-R group; p-values <0.05 are highlighted with coloured triangles. Colour coding: VAX group (green); NO-VAX (grey).
https://doi.org/10.1371/journal.ppat.1011787.g003 Neither a significant clustering between VAX and NO-VAX groups is observed by multivariate statistics (PCA and RF, S3 Fig) nor significantly different levels of metabolites or lipoproteins are found between the two groups. The two groups show the same metabolic and lipoproteomic alterations when compared to the COVID-19-R group (Fig 3B–3C). Only succinate shows an opposite trend, being significantly increased when considering the VAX group vs. the recovered subjects and decreased (but not significantly, p-value > 0.05) for the comparison between the NO-VAX group vs. the recovered subjects. Lipoprotein sub-fractions and particle numbers are reported in S4 Fig.
Sex-specific differences in COVID-19 disease Sex differences in the disease fingerprint are hereafter evaluated. Since our data did not show significant differences in the metabolomic and lipoproteomic fingerprints associated with the main variants of the virus or with the vaccination status, all the 510 samples (274 from the previous study plus 236 newly enrolled) collected from COVID-19 patients are analyzed together, S1 Table. The presence of sex-specific differences in the plasma metabolomic and lipoproteomic profiles of healthy subjects has been characterized in detail [32]. Thus, to extract only COVID-19 related differences, 287 male and 223 female COVID-19 patients are compared to 52 male and 43 female recovered subjects, respectively, Fig 4A. PPT PowerPoint slide
PNG larger image
TIFF original image Download: Fig 4. Sex-related metabolomic and lipoproteomic profiling of COVID-19 subjects. (A) Number of samples used for the analysis. (B-C) Values of Log 2 fold change (FC) of quantified metabolites (B) and lipoprotein parameters (main parameters, calculated features and main fractions) (C). Positive/negative values have higher/lower concentration in plasma samples from the male (M) or female (F) groups with respect to the M or F COVID-19-R group, respectively; p-values <0.05 are highlighted with coloured triangles. Colour coding: male group (blue); female group (pink).
https://doi.org/10.1371/journal.ppat.1011787.g004 Again, the strong and solid metabolic fingerprint induced by the infection can be detected in both male and female COVID-19 subjects, Fig 4B–4C. No significant differences in the dysregulation trends of either metabolites or lipoprotein parameters emerge between the two groups, with the only exceptions of formic acid and Gln that are slightly but significantly decreased only for male patients. Lipoprotein sub-fractions and particle numbers are reported in S5 Fig.
Comorbidities-dependent variations Finally, the contribution of the patient’s comorbidities on the dysregulated levels of metabolites and lipoproteins was also evaluated. S1 Table lists the number of subjects affected by the comorbidities that can be considered as the main risk-factors of severe prognosis [6,7], i.e. asthma or chronic obstructive pulmonary disease (COPD), cardiovascular diseases (such as coronary artery disease (CAD), congestive heart failure (CHF), hypertension, type 2 diabetes (T2DM), dyslipidaemia, chronic kidney disease (CKD) and immune deficiency). Not unexpectedly, subjects affected by T2DM are characterized by higher glucose and mannose concentrations and those affected by CKD have high levels of creatinine (S8 Fig). Whereas creatinine was not found significantly different between the COVID-19 and COVID-19-R groups, high levels of mannose and glucose are very important markers of the COVID-19 signature and even for the discrimination between severe and fatal patients. Importantly, the observed differences in the levels of these two sugars are preserved when diabetic subjects are discarded from the analysis (S9 Fig). It is also important that the trend of increasing concentration for mannose and even more for glucose when going from mild to fatal is more evident for the T2DM subjects that in the no-T2DM patients (S9 Fig).
The follow-up cohort For 154 out of 510 COVID-19 patients, a plasma sample was also collected during the follow-up visit, 2–6 months after the first negative swab (follow-up group S2 Table). Thus, for these patients two plasma samples were available, one at the onset of the acute infection (T1) and the second months after the negative swab (T2). They include 111 patients from the wt-α-β group (including 72 subjects with T1 analyzed already in the previous study), 33 from the δ group, 7 from the ο group, and 3 not attributable to a specific variant. None of these subjects was diagnosed with long-COVID. The availability of two samples from the same individual allowed us to use a paired test. The same results emerge as from the above-described comparisons between the two independent groups of COVID-19 positive subjects and the COVID-19-R subjects. Consistently, no significant differences in metabolite and lipoprotein concentrations are detected when comparing the T2 samples with those from the COVID-19-R subjects. Indeed, the dysregulated molecules during the infection (i.e., at T1), are essentially reverted to within the normality range in samples collected at T2, S10 Fig. This behaviour is independent of the SARS-CoV-2 variant.
[END]
---
[1] Url:
https://journals.plos.org/plospathogens/article?id=10.1371/journal.ppat.1011787
Published and (C) by PLOS One
Content appears here under this condition or license: Creative Commons - Attribution BY 4.0.
via Magical.Fish Gopher News Feeds:
gopher://magical.fish/1/feeds/news/plosone/