COVID-19 infection alters kynurenine and fatty acid metabolism, correlating with IL-6 levels and renal status

Reprogramming of host metabolism supports viral pathogenesis by fueling viral proliferation, by providing, for example, free amino acids and fatty acids as building blocks. To investigate metabolic effects of SARS-COV-2 infection, we evaluated serum metabolites of COVID-19 patients (n = 33; diagnosed by nucleic acid testing), as compared to COVID-19-negative controls (n = 16). Targeted and untargeted metabolomics analyses identified altered tryptophan metabolism into the kynurenine pathway, which regulates inflammation and immunity. Indeed, these changes in tryptophan metabolism correlated with interleukin-6 (IL-6) levels. Widespread dysregulation of nitrogen metabolism was also seen in infected patients, with altered levels of most amino acids, along with increased markers of oxidant stress (e.g., methionine sulfoxide, cystine), proteolysis, and renal dysfunction (e.g., creatine, creatinine, polyamines). Increased circulating levels of glucose and free fatty acids were also observed, consistent with altered carbon homeostasis. Interestingly, metabolite levels in these pathways correlated with clinical laboratory markers of inflammation (i.e., IL-6 and C-reactive protein) and renal function (i.e., blood urea nitrogen). In conclusion, this initial observational study identified amino acid and fatty acid metabolism as correlates of COVID-19, providing mechanistic insights, potential markers of clinical severity, and potential therapeutic targets.


C L I N I C A L M E D I C I N E
stratified by circulating levels of IL-6, and correlated to inflammatory markers and renal function. Though observational, this study identifies substantial alterations in fatty acid metabolism in infected patients and unanticipated alterations of glucose homeostasis and amino acid metabolism; specifically, dysregulated tryptophan, arginine, and sulfur metabolism are potential therapeutic targets.

Results
Forty-nine subjects were enrolled in this observational study: 33 COVID-19-positive subjects and 16 COVID-19-negative subjects; the latter included "never positive" subjects and convalescent plasma donors ( Figure 1A). Owing to the skewed demographics of the COVID-19-positive patients (prevalence and mortality being ~2-fold higher in men) (28), 76% of the subjects in this group were male, 56.5 ± 18.1 years old (mean ± standard deviation) versus 38% male, 37.8 ± 11.6 years old for the control group. Of the COVID-19-positive patients, 5 subjects had IL-6 levels of less than 10 pg/mL; 10 subjects had IL-6 levels of 10-65 pg/mL, and 18 subjects had IL-6 levels of more than 90 pg/mL. Subjects in these 3 subgroups had no significant differences in age or sex. Reference ranges for IL-6 levels in healthy controls are less than 5 pg/mL ( Figure 1B). Of all the COVID-19-positive subjects, only 6 had renal function parameters within the reference range (blood urea nitrogen [BUN] of 7-26 mg/dL and creatinine of 0.5-0.95 mg/dL), of which 3 were in the low and 2 in the medium IL-6 group (Supplemental Table 1 and Supplemental Figure 1A; supplemental material available online with this article; https://doi.org/10.1172/jci.insight.140327DS1). Targeted metabolomics analyses were performed on sera using ultra-high-pressure liquid chromatography-mass spectrometry (UHPLC-MS), and the results are reported extensively in Supplemental Table 1. The metabolic phenotypes of sera from COVID-19-positive patients differed substantially from controls, when examined by partial least squares discriminant analysis (PLS-DA; Figure 1C). Hierarchical clustering analysis ( Figure  1D) highlighted significant associations of COVID-19 and IL-6 levels with amino acid metabolism, purines, acylcarnitines, and fatty acids; a vectorial version of this figure is provided in Supplemental Figure 2. Volcano plots ( Figure 1E) from targeted metabolomics analyses highlighted the top metabolites that increased (blue) or decreased (red) in the sera of COVID-19-positive patients, as compared with controls. To expand on these observations, untargeted metabolomics analyses were performed based on high-resolution, accurate intact mass, isotopic patterns, and MS2 fragmentation patterns, as described previously (29); the "top" correlates with clinical laboratory parameters are reported in Supplemental Table 2. In Figure 2A, volcano plot analyses highlighted 3034 and 2484 differential metabolites for negative and positive ion modes between COVID-19-positive and -negative subjects, respectively. PLS-DA analysis based on untargeted metabolomics data supported the targeted results and further separated COVID-19-positive and COVID-19-negative subjects ( Figure 2B), the former group separating from the latter across principal component 1 (18.7% of the total variance) as a function of IL-6 levels. Interestingly, volcano plot analyses comparing COVID-19-positive subjects with low, medium, and high serum IL-6 levels, as compared with controls, showed that the low IL-6 group was most comparable to controls (199 and 218 metabolites significantly decreasing and increasing), whereas the medium (218 and 765) and high IL-6 (206 and 789) groups differed the most ( Figure 2C). Metabolite Set Enrichment Analysis (30) of merged targeted and untargeted metabolomics data revealed a significant impact of COVID-19 on amino acid metabolism, especially the pathways involving tryptophan, aspartate, arginine, tyrosine, and lysine ( Figure 2D). Top hits from these pathways were mapped against Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway map hsa01100 in Supplemental Figure 3, and tryptophan metabolism is highlighted in Figure 2E. Targeted absolute quantitative measurements, as determined by stable isotope-labeled internal standards, are provided in Supplemental Table 1 for a subset of metabolites from the most significantly affected pathways.
Alteration of tryptophan metabolism in sera of COVID-19-positive patients. Tryptophan metabolism was the top pathway affected by COVID-19 in the analysis of targeted and untargeted metabolomics data ( Figure  2E). As such, focused analysis of this pathway highlighted significant decreases (inversely proportional to IL-6 concentration) in tryptophan, serotonin, and indolepyruvate levels ( Figure 3A). In contrast, increases in kynurenine, kynurenic acid, picolinic acid, and nicotinic acid, but not anthranilate, suggested hyperactivation of the kynurenine pathway ( Figure 3A). Although it is beyond the scope of this work to identify diagnostic metabolic markers of COVID-19 positivity, receiver operating characteristic (ROC) curves were calculated for the absolute quantitative measurements of serum tryptophan and kynurenine in these sera ( Figure 3B); interestingly, tryptophan levels lower than 105 μM and kynurenine levels higher than 5.3 μM had areas under the curves more than 95% in distinguishing between the 2 groups (positive vs. negative). insight.jci.org https://doi.org/10.1172/jci.insight.140327

C L I N I C A L M E D I C I N E
The potential relevance of this observation will require additional study in the future.
Dysregulation of amino acid metabolism in patients with COVID-19 as a function of IL-6 levels. Beyond tryptophan, metabolomics analyses revealed a significant impact of SARS-CoV-2 infection on serum amino acid levels ( Figure 4). Specifically, gluconeogenic (e.g., alanine, glycine, serine, glutamine, and histidine) and (A) Forty-nine subjects were studied, of which 16 were COVID-19-negative and 33 were COVID-19-positive patients, as determined by nucleic acid testing of nasopharyngeal swabs. IL-6 levels were determined during routine clinical care using a clinically validated ELISA (B), and the results were used to divide COVID-19-positive patients into groups with low (≤10 pg/mL), medium (10-65 pg/mL) and high (>90 pg/mL) IL-6 levels (although IL-6 levels were treated as a continuous variable, no patients had a result of >65 but ≤90 pg/mL). Sera were obtained from these subjects for metabolomics analyses. Asterisks indicate significance by ANOVA (1-way ANOVA with Tukey's multiple comparisons, *P < 0.05; **P < 0.01; ***P < 0.001; ****P < 0.0001). sulfur-containing amino acids (e.g., cysteine, taurine) tended to decrease, especially in the moderate-high IL-6 group ( Figure 4A). In contrast, oxidized forms of sulfur-containing amino acids (e.g., methionine sulfoxide, cystine), as well as of arginine, increased ( Figure 4A). Although no significant changes were noted in methionine levels, increases in acetyl-methionine and hydroxyproline, by-products of proteolysis and collagen catabolism, respectively, were observed, particularly in patients with COVID-19 with the highest IL-6 levels ( Figure 4B). Despite increases in arginine, the urea cycle metabolic intermediates ornithine and citrulline decreased in patients with COVID-19 ( Figure 4C). In addition, creatine and creatinine, as well as the polyamines spermidine and acetyl-spermidine, increased significantly, particularly in patients with medium and high IL-6 levels ( Figure 4, C and D); this is suggestive of renal dysfunction in these groups, which was confirmed by clinical laboratory measurements of creatinine and BUN (Supplemental Table 1). Evidence of hyperglycemia, hemolysis, and lipid abnormalities in sera of patients with COVID-19. Because altered amino acid metabolism suggested dysregulated nitrogen metabolism, we determined whether this was also accompanied by significant alterations in carbon metabolism. Notably, all patients with COVID-19 in this cohort, independent of IL-6 levels, exhibited hyperglycemia ( Figure 5 -also confirmed by clinical laboratory measurements -Supplemental Table 1). Some metabolic intermediates of the glycolysis and pentose phosphate pathways were increased in subjects with the highest IL-6 levels; this was suggestive of potential hemolysis (e.g., ribose phosphate; Figure 5), although no significant increases in lactate dehydrogenase (LDH) were noted across the 3 groups of patients with COVID-19 (Supplemental Figure 1A). However, unexpected decreases in lactate levels were noted in the same subjects, in the absence of differences in the levels of all carboxylic acids except α-ketoglutarate ( Figure 5); this suggested altered transamination homeostasis consistent with an altered nitrogen balance.
In addition, sera from patients with COVID-19 demonstrated significant changes in levels of acylcarnitines and free fatty acids ( Figure 6). Specifically, all short-and medium-chain acylcarnitines, but not acyl-C18:3, were significantly decreased in all patients with COVID-19, independent of IL-6 levels. Finally, all fatty acids, except for nonanoic acid, were increased in all patients with COVID-19, independent of IL-6 levels ( Figure 6).

Metabolic correlates of laboratory markers of inflammation and renal function.
To determine the potential clinical relevance of these metabolic findings, metabolite levels were correlated to clinical laboratory markers of inflammation (i.e., IL-6 and C-reactive protein [CRP], Figure 7, A and B) and renal function (i.e., BUN and creatinine, Figure 7, C and D). Notably, several acylcarnitines, kynurenine, and methionine sulfoxide were among the top correlates to IL-6 ( Figure 7A). Free fatty acid and tryptophan levels were among the top negative correlates to CRP levels, which were also directly proportional to levels of picolinic acid, inosine, and short-chain acylcarnitines ( Figure 7B). Interestingly, acylcarnitines, but not free carnitine, ranked among the top positive correlates to BUN levels, along with creatinine (as expected), carboxylic acids (e.g., fumarate, itaconate), and oxidized purines (e.g., inosine; Figure 7C). As a supportive internal validation, creatinine, as measured by mass spectrometry, had a significant positive correlation (Spearman's r = 0.93) with clinical laboratory measurements of creatinine ( Figure 7D). In addition, purine deamination products (inosine and hypoxanthine) and several acylcarnitines ranked among the top correlates to creatinine ( Figure 7D).

Discussion
The present study provides what is likely the first observational metabolic characterization of sera from COVID-19-positive patients. The results show marked alterations in nitrogen metabolism (especially amino acid homeostasis, catabolism, and transamination) and carbon metabolism (especially affecting levels of glucose and free fatty acids). Metabolites in these pathways significantly correlated with circulating levels of inflammatory markers (i.e., IL-6 and CRP) and renal function (i.e., BUN and creatinine). Interestingly, the major effects of COVID-19 were on tryptophan metabolism and the kynurenine pathway. The rate-limiting step in this pathway is indole 2,3-dioxygenase (IDO1), a pleiotropic enzyme with complex effects, including immunoregulation, typically as a negative regulator of inflammation and immunization (31,32). The kynurenine/tryptophan ratio is a general measure of IDO activity, although its meaning alters depending upon metabolic activity and its medium, and high serum IL-6 levels, as compared with controls. (D) Pathway analysis of untargeted metabolomics data identified significant effects of COVID-19 on amino acids, especially regarding tryptophan, aspartate, arginine, tyrosine, and lysine metabolism. (E) Tryptophan metabolism, one of the top hits from pathway analysis, is mapped against KEGG pathway map hsa01100. insight.jci. effects on kynurenine and tryptophan half-lives. Of particular relevance to our study, IDO helps limit lung inflammation in the context of IL-6, and IDO deletion severely exacerbates inflammatory lung pathology in mice (33). IDO regulates immune responses by depleting tryptophan and by generating tryptophan metabolites that activate the aryl hydrocarbon receptor; these metabolic alterations were observed in the current report. Moreover, IDO activity (inferred from metabolite levels) inversely correlated with IL-6 levels. Although causes and effects cannot be inferred by our results, it is reasonable to hypothesize that IDO activity is mechanistically involved in COVID-19 inflammation. Consistent with this, IDO is expressed by the lung mucosa, including antigen-presenting and epithelial cells (34). IDO1 is downstream of interferon signaling (involving STAT6 and NF-κB) (31), which is activated by viral infection (35). Intriguingly, however, type I interferon production is suppressed in SARS-CoV-2 infection; in addition, interferon-γ production correlates with disease severity (36)(37)(38)(39). In one predicted mechanism, interferon induces ACE2 receptor expression (40), which also binds SARS-CoV-2 spike protein S to mediate infection of host cells (11)(12)(13). Thus, interferon and its signaling partners can have complex and potentially opposing effects in SARS-CoV-2 pathology and the ensuing inflammation. Importantly, IDO activity is also regulated by oxidative stress, being enhanced by superoxide and NO. In summary, multiple factors may affect how IDO limits inflammation and affects adaptive immunity. Therefore, IDO activity may increase survival from COVID-19 by limiting pulmonary inflammation while simultaneously decreasing development of protective adaptive immunity to SARS-CoV-2.
Recently, we identified IDO pathway activation in individuals with trisomy 21 (Down syndrome) (31). Indeed, 4 of the 6 receptors for interferon-α and -γ are encoded by genes on chromosome 21. In addition, individuals with Down syndrome have respiratory comorbidities (e.g., pulmonary hypertension, chronic Similarly, significant increases were seen in the levels of some intermediates of the glycolytic and pentose phosphate pathways in sera of COVID-19-positive patients. No significant changes were noted in serum levels of carboxylic acids across the groups. Asterisks indicate significance by ANOVA (1-way ANOVA with Tukey's multiple comparisons, *P < 0.05; **P < 0.01; ***P < 0.001; ****P < 0.0001). The y axis in the dot plots indicates AU. PPP, pentose phosphate pathway. obstructive pulmonary disease), gastrointestinal disorders, diabetes, and accelerated aging (31), all risk factors for, or manifestations of, severe COVID-19 illness. In the Down syndrome population, therapies targeting IDO1 or its upstream pathways (i.e., Janus-activated kinase [JAK]), were proposed to ameliorate some of the accompanying comorbidities (41). Thus, the current findings suggest the possible relevance of JAK inhibitors in COVID-19 (42). Epidemiological studies of COVID-19 in the Down syndrome population could also provide clues regarding interferon signaling in SARS-CoV-2 infection.
Nonetheless, blunted interferon type I and III responses were reported in COVID-19 (43), which could, in part, be explained by SARS-CoV-2 proteins ORF3B, ORF6, and nsp nuclease; these were implicated in antagonizing interferon and cleaving host mRNA to prevent ribosomal loading and causing host shutoff, respectively (44,45). Similarly, impaired type I interferon activity and exacerbated inflammatory responses were also described in patients with severe COVID-19 (46), although progressive increases in disease severity correlated with transcript levels of JAK1, STAT1 and STAT2, interferon-α2, interferon-α receptors 1 and 2, and interferon regulatory factors 1, 4, 5, and 7. Thus, strategies were designed to supplement interferon-α/β to treat COVID-19 (47), which significantly reduce persistence of virus in the upper respiratory tract and reduce the duration of elevated blood levels of IL-6 and CRP. Although speculative, metabolites in the kynurenine pathway are critical immunomodulators, which contribute, along with arginine (one of the most affected pathways identified herein), to immunosuppressive activity of dendritic cells (48) and to CD8 + T cell suppression in breast cancer (49,50). Therefore, activating this pathway may allow SARS-CoV-2 to evade immunity. Moreover, dysregulated tryptophan metabolism is seen in inflammation (51) and aging, resulting in decreased NAD synthesis in older animals (52) and humans, an observation suggesting nicotinamide riboside as a dietary supplement (53); this provides a potential metabolic explanation for the severity of COVID-19 in older subjects. It is worth noting that proinflammatory signaling favors proteolysis and amino acid catabolism (represented herein by increased acetyl-methionine and hydroxyproline and altered levels of free amino acids), which can be antagonized, in part, by antiinflammatory cytokines (e.g., IL-37) (54) or inflammasome inhibitors (55). These antiinflammatory cytokines show promise in regulating glycemia in aging mice (56), which is relevant to finding increased glycemia in our patients with COVID-19; this may also help explain why diabetes is a risk factor for severe COVID-19 illness (57) and outcomes (58). Interestingly, similar increases in circulating glucose levels are found in critically ill patients with diabetes (59)(60)(61). Similarly, in trauma, extreme lipolysis is seen minutes to hours after trauma/shock, producing increased circulating free fatty acids (59), which we also identified in patients with COVID-19 (e.g., long-chain free fatty acids), especially those with high IL-6 levels. Phospholipase A2 (PLA2), which cleaves fatty acids in the sn-2 position of phospholipids, is required for replication of other coronaviruses, such as Middle Eastern Respiratory Syndrome-CoV (MERS-CoV) and human coronavirus 229E (26). Elevated levels of long-chain polyunsaturated fatty acids (PUFAs), which are highly enriched in the sn-2 position of phospholipids, in the absence of elevated longchain acylcarnitines, which could reflect impaired mitochondrial fatty acid oxidation, suggest that this pathway may be upregulated in COVID-19. Increased PLA2 activity is associated with increased production of bioactive lipids by metabolism of omega-6 PUFAs and lysophosphatidylcholine. Some resulting metabolites (e.g., eicosanoids, oxylipins) may help propagate infection and thromboinflammatory complications (62), by directly supporting viral capsid membrane synthesis (25)(26)(27), by targeting proinflammatory immune cells, or by activating platelets (63). Thus, inhibiting PLA2 activity or increasing omega-3 PUFA levels, which are metabolized to bioactive lipids with antiinflammatory activity, might reduce disease severity.
Inflammatory stimuli promote platelet hyperreactivity in older individuals (23). Notably, polyamine by-products of nitrogen metabolism (e.g., spermidine, acetyl-spermidine) were greatly increased in COVID-19 patient sera. Primary amines, such as tranexamic acid, prevent hyperfibrinolysis in actively bleeding patients by interacting with tissue plasminogen activator (64). One can speculate that inflammation-induced alterations in nitrogen metabolism may participate in dysregulated coagulation leading to thrombus formation in severe COVID-19. In contrast, other metabolites associated with hypercoagulability, such as homocysteine (65) (which increases in Down syndrome, ref. 66), were decreased in patients with COVID-19. The latter may be explained by increased oxidant stress in patients with COVID-19, suggested by increased serum levels of methionine sulfoxide and cystine, along with decreased antioxidants, such as cysteine and taurine (67) (but, interestingly, not reduced glutathione). Conversely, decreased serum acylcarnitine levels, which may have an anticoagulant function (68), are consistent with a potential hypercoagulable phenotype.
Interestingly, increased serum sphingosine 1-phosphate levels were observed in patients with COVID-19. Circulating levels of this metabolite are significantly influenced by red blood cell (RBC) sphingosine kinase 1 activity in response to hypoxia (69), whereas alterations to this pathway mediate responses to angiotensin II-induced stimulation in the hypoxic kidney in chronic kidney disease (70); the latter may be relevant in light of renal dysfunction in some patients with COVID-19 in our study. In contrast, one would have predicted significantly increased serum levels of lactate and carboxylic acid markers of hypoxia (e.g., succinate) in patients with COVID-19, similar to that seen in patients following ischemic (71) or hemorrhagic shock (72). Surprisingly, no significant increases in serum lactate or succinate levels were observed, which may be because our patients were either normoxic or receiving intense respiratory therapy at the time of blood draw. Although unrelated to IL-6, fumarate and itaconate positively correlated with BUN and inosine, suggesting potential crosstalk between deaminated purine salvage reactions and carboxylic acid metabolism in renal dysfunction in COVID-19, similar to acute kidney ischemia (73). In addition, hemolysis markers (e.g., metabolic intermediates of the glycolytic and pentose phosphate pathways) correlated with IL-6 levels, suggesting interactions between inflammation, disease severity (e.g., correlation between LDH and BUN; Supplemental Figure 1B), and, ultimately, the effect of therapeutic interventions and mechanical/oxidant stress on circulating RBCs. As such, future studies will assess COVID-19 effects on RBCs.
The present study has several limitations. First, analyses were performed on a small set of sera obtained as a by-product of routine clinical laboratory testing. Although immediately refrigerated and stored for less than 24 hours, technical bias may have been introduced owing to collection procedures. Protocols are being implemented to prospectively collect and bank fresh samples of serum and plasma, RBCs, and buffy coats. In addition, males were disproportionally represented (75% of patients with COVID- 19). As such, no sex-specific analysis was performed, and future investigations on larger, prospectively enrolled cohorts are underway to address this issue. Similarly, the mean age of COVID-19-positive patients was 56, whereas the control population was generally younger (mean: 33 years old). Nonetheless, although aging is associated with increased systemic inflammation, comparisons to prior studies on the effect of aging on the plasma/ whole-blood metabolome in humans (74,75) and animal models (52,54,56) does not suggest that the observations reported here are due to aging alone. In addition, the patients with COVID-19 in this study had significant disease, by definition, because they were all inpatients; future studies will investigate patients with milder disease and asymptomatic infected patients. Conversely, the control group did not include patients requiring hospitalization for respiratory infections distinct from SARS-CoV-2. As such, some of the biosignatures reported herein (e.g., fatty acids, refs. 25, 26) may be explained by the effects of responses to viral infection, in general, rather than being specifically associated with SARS-CoV-2. In addition, all patients with COVID-19 in this study exhibited some degree of renal dysfunction, which may affect serum metabolism by compromising metabolite filtration capacity. Nonetheless, it is worth noting that the present study did not identify an effect of the inflammatory state (e.g., IL-6 levels) on the characteristic increases in purine breakdown and deamination/oxidation products seen in the context of chronic kidney disease (70) or acute kidney ischemia (73,76). However, accumulation of metabolites in this pathway did correlate with BUN and creatinine, consistent with a role of renal dysfunction in regulating circulating levels of oxidized purines. In addition, the patients in this study with more severe disease were most likely simultaneously receiving multiple therapeutic interventions, which may have confounded some of the findings in this group. Finally, it will be important to study serial samples from patients throughout their clinical course to evaluate the serum metabolome as a function of clinical status (e.g., stage of infection, on or off a ventilator, receiving hemodialysis or not) and in response to various therapies (e.g., treatment with IL-6 receptor inhibitors).
In conclusion, we provide the first observational report of the comprehensive serum metabolome in patients with COVID-19 as a function of IL-6 level, the latter as a proxy for disease severity and as a promising therapeutic target. These results suggest potential therapeutic targets in systems metabolism, including (a) IDO1 or its upstream modulators (e.g., JAK, STAT, and/or interferon signaling); (b) amino acid catabolism as a function of oxidant stress (e.g., taurine) or inflammation-induced proteolysis (through dietary supplementation with taurine/N-acetylcysteine, using antibodies against proinflammatory cytokines or antagonists to their receptors, such as, anakinra, or using antiinflammatory cytokines or inflammasome inhibitors) (54)(55)(56); and (c) free fatty acid levels and unsaturation, which are amenable to dietary or pharmacological intervention (e.g., fish oil in the diet or fatty acid desaturase inhibitors).

Methods
Blood collection and processing. Subjects seen at NewYork-Presbyterian/Columbia University Irving Medical Center included 33 COVID-19-positive patients, as determined by SARS-CoV-2 nucleic acid testing of nasopharyngeal swabs; in this group, the severity of the disease was inferred from serum IL-6 levels, which were determined by a CLIA-certified ELISA-based assay. The control group included 16 subjects, all of whom were SARS-CoV-2 negative by nasopharyngeal swab at the time of the blood draw. Some patients in this group were "never positive" subjects, and some were COVID-19 convalescent patients who were previously positive but currently negative, as determined by testing nasopharyngeal swabs, and at least 14 days postresolution of symptoms. Serum was obtained from freshly drawn blood after an overnight hold at 4°C and all samples were then deidentified. Sera were extracted via a modified Folch method (chloroform/ methanol/water 8:4:3), which completely inactivates other coronaviruses, such as MERS-CoV (77). Briefly, 20 μL of serum was diluted in 130 μL of LC-MS-grade water, 600 μL of ice-cold chloroform/methanol (2:1) was added, and the samples were vortexed for 10 seconds. Samples were then incubated at 4°C for 5 minutes, quickly vortexed (5 seconds), and centrifuged at 14,000 g for 10 minutes at 4°C. The top (i.e., aqueous) phase was transferred to a new tube for metabolomics analysis.
UHPLC-MS metabolomics. Metabolomics analyses were performed using a Vanquish UHPLC coupled online to a Q Exactive mass spectrometer (Thermo Fisher Scientific, Bremen, Germany). Samples were analyzed using 5-, 15-, and 17-minute gradients, as described (78,79). For targeted quantitative experiments, extraction solutions were supplemented with stable isotope-labeled standards, and endogenous metabolite concentrations were quantified against the areas calculated for heavy isotopologs for each internal standard (78,79). Data were analyzed using Maven (Princeton University, Princeton, New Jersey, USA) and Compound Discoverer 2.1 (Thermo Fisher Scientific).
Statistics. Graphs and statistical analyses were prepared with GraphPad Prism 8.0 (GraphPad Software, Inc, La Jolla, California, USA), GENE E (Broad Institute, Cambridge, Massachusetts, USA), and Metabo-Analyst 4.0 (80). Spearman's correlations and related P values were calculated with RStudio (Boston, Massachusetts, USA). One-way ANOVA with Tukey's multiple comparisons follow-up test was used for comparing metabolite concentrations across different groups. P < 0.05 was considered statistically significant. Data are shown as mean ± SD unless indicated.
Study approval. This observational study was conducted according to the Declaration of Helsinki, in accordance with good clinical practice guidelines, and approved by the Columbia University Institutional Review Board. Residual samples were obtained following the completion of routine clinical laboratory testing for each subject and deidentified before evaluation for research purposes.
Additional approval was obtained from the University of Colorado -Anschutz Medical Campus Institutional Review Board as part of the COVID Omics project (protocol 1300v3).

Author contributions
TT, ROF, SLS, and EAH designed the study. TT, ROF, and EAH collected and processed the samples. DS, LB, JAR, TN, and ADA performed metabolomics analyses and prepared the figures. JAR, TN, and ADA performed data analysis and prepared figures and tables. ADA wrote the first draft of the manuscript, which was significantly revised by SLS, TT, JCZ, and KEH and finally approved by DS, LB, JAR, TN, TT, JCZ, KEH, ROF, EAH, KCH, ADA, and SLS.