Kynurenic acid underlies sex-specific immune responses to COVID-19 ================================================================== * Yuping Cai * Daniel J. Kim * Takehiro Takahashi * David I. Broadhurst * Shuangge Ma * Nicholas J.W. Rattray * Arnau Casanovas-Massana * Benjamin Israelow * Jon Klein * Carolina Lucas * Tianyang Mao * Adam J. Moore * M. Catherine Muenker * Jieun Oh * Julio Silva * Patrick Wong * Yale IMPACT Research team * Albert I. Ko * Sajid A. Khan * Akiko Iwasaki * Caroline H. Johnson ## Abstract Coronavirus disease-2019 (COVID-19) has poorer clinical outcomes in males compared to females, and immune responses underlie these sex-related differences in disease trajectory. As immune responses are in part regulated by metabolites, we examined whether the serum metabolome has sex-specificity for immune responses in COVID-19. In males with COVID-19, kynurenic acid (KA) and a high KA to kynurenine (K) ratio was positively correlated with age, inflammatory cytokines, and chemokines and was negatively correlated with T cell responses, revealing that KA production is linked to immune responses in males. Males that clinically deteriorated had a higher KA:K ratio than those that stabilized. In females with COVID-19, this ratio positively correlated with T cell responses and did not correlate with age or clinical severity. KA is known to inhibit glutamate release, and we observed that serum glutamate is lower in patients that deteriorate from COVID-19 compared to those that stabilize, and correlates with immune responses. Analysis of Genotype-Tissue Expression (GTEx) data revealed that expression of kynurenine aminotransferase, which regulates KA production, correlates most strongly with cytokine levels and aryl hydrocarbon receptor activation in older males. This study reveals that KA has a sex-specific link to immune responses and clinical outcomes, in COVID-19 infection. ## Main Sex-related differences in coronavirus disease-2019 (COVID-19) severity and morbidity exist, with the male sex being a risk factor1; male COVID-19 patients have an increased risk of admission (OR 1.68, 95%CI = 1.45-1.90) and in-hospital mortality (OR 1.87, 95%CI = 1.33-2.63)1. It was recently shown that hospitalized patients with moderate SARS-CoV-2 infection have elevated levels of certain inflammatory cytokines and chemokines, and sex-differences exist in these immune responses2. Furthermore, across all ages, female patients at baseline had a more robust T cell activation than males. Loss of T cell activation was correlated with older age in males, and this poorer T cell response was correlated with worse disease outcomes in males only2. Therefore, males and females have clear differences in COVID-19 immune responses that correlate with clinical course. Since immune responses are regulated, in part, by metabolites, it is possible that sex-related differences in metabolism could affect the host immune response to SARS-CoV-2 infection. For instance, specific metabolites are required for macrophage, neutrophil, and T cell functions, enhancing glycolytic and fatty acid synthesis pathways in these cells3. Conversely, immune stimulation can also elicit metabolic reprograming in cells, thereby affecting disease trajectory by altering metabolite abundance4. In addition to the metabolic requirements of the host immune system, viruses also require host-derived metabolites and lipids5. Thus, utilization of metabolic substrates for viral replication could affect metabolite availability required for immune responses. ### Metabolites correlate with COVID-19 To address how metabolites might mediate the sex-related differences in COVID-19 immune response, we first used an untargeted metabolomics workflow with multivariable logistic regression to identify serum metabolites associated with COVID-19. Serum samples were collected from COVID-19 patients (n = 22 females and n = 17 males) on the day of enrollment into the IMPACT study at Yale New Haven Hospital (CT, USA). Samples were taken from patients 3-7 days after hospital admission after confirmation of COVID-19 infection who (1) were not immediately triaged to the intensive care unit, (2) had not received tocilizumab, and (3) had not received high dose corticosteroids (Cohort A described in Takahashi et al.2). Uninfected healthcare worker (HCW) controls (n = 10 females and n = 10 males) were included in the analysis. There was a statistically significant difference in age between the COVID-19 patients and HCWs, which was adjusted for in our models (Extended Data Table 1). We first carried out metabolite identification on detected signals that were present in the serum metabolome of quality control samples pooled from both COVID-19 patients and HCWs. We positively identified 75 metabolites with the highest confidence (Extended Data Table 2). Regression analysis revealed that 17 metabolites were associated with COVID-19 status after adjustment for age, BMI, sex, and multiple comparisons (Extended Data Table 3). Glutamate, cysteine-S-sulfate, palmitoleic acid, arachidonic acid, lysophosphatidylethanolamine (LPE) (22:6), uracil and myristic acid were positively associated with COVID-19. Whereas glutamine, 3-methylxanthine, tryptophan, proline, citrulline, homoserine, 2,3-dihydroxybenzoic acid, lysophosphatidic acid (LPA) (18:2), LPA (20:2), and lysophosphatidylcholine (14:0) were negatively associated with COVID-19. ### Metabolites correlate with immune response in a sex-specific manner Next, we examined how the 75 positively identified serum metabolites from both COVID-19 patients and HCWs might correlate with immune markers (cytokines and chemokines levels in plasma, T cells, B cells, NK-T cells, NK cells, monocytes, macrophages, and dendritic cells in peripheral blood mononuclear cells (PMBCs)) that were previously measured from the same individuals2 (**Fig. 1**). In COVID-19 patients, we observed 36 correlations between immune markers and metabolites with |R|>0.5 and *p* value< 0.05 (**Extended Data Fig. 1,Supplementary Table 1**). However, after stratification by sex, additional correlations were observed between metabolites and immune markers for each sex independently, suggesting that sex-specific immune responses could be linked to metabolism (**Extended Data Fig. 1, Supplementary Table 2**). ![Fig. 1.](http://medrxiv.org/https://www.medrxiv.org/content/medrxiv/early/2020/09/08/2020.09.06.20189159.1/F1.medium.gif) [Fig. 1.](http://medrxiv.org/content/early/2020/09/08/2020.09.06.20189159.1/F1) Fig. 1. Chord diagram of correlations between metabolites and immune markers in COVID-19 patients. Spearman correlations > 0.5 or < −0.5 are displayed and with p< 0.05. Further examination revealed that kynurenic acid (KA), an endogenous ligand of the aryl hydrocarbon receptor (AhR) that regulates immunes responses6, had positive correlations with immune markers (**Fig. 1**). Many of these positive correlations were observed in male patients including IL1RA, IL6, IL10, TNFα, M-CSF, SCF, CX3CL1, CXCL9, CXCL13, CCL1, CCL21, and CCL22. In addition, KA in males was negatively associated with T cell number, naïve CD8 T cells, CD4 effector memory (CD4Tem), and CD8 effector memory (CD8Tem) T cells (**Figs. 2a, 2b**). In female patients, KA was positively associated only with IL12p40, CCL3, CXCL9, and SCF (**Fig. 2b**). In summary, sex-specific differences in correlations between metabolites and immune responses were observed in COVID-19 patients, wherein KA had the most prominent connection to immune response in males. ![Fig. 2.](http://medrxiv.org/https://www.medrxiv.org/content/medrxiv/early/2020/09/08/2020.09.06.20189159.1/F2.medium.gif) [Fig. 2.](http://medrxiv.org/content/early/2020/09/08/2020.09.06.20189159.1/F2) Fig. 2. Tryptophan pathway metabolites and immune responses. **a,** Correlation between kynurenic acid (KA) and immune markers in males with COVID-19 (Pt. M, n = 17) and females with COVID-19 (Pt. F, n = 22). 95% confidence intervals (CIs) for the correlation coefficients are indicated as shaded areas colored according to patient sex. **b**, Tryptophan (T) metabolism pathway schematic. **c**, Heatmap showing correlation between tryptophan metabolites and immune markers in males and females with COVID-19. Spearman correlations > 0.5 or < –0.5 are displayed, p< 0.05. **d**, Correlation between age and KA:kynurenine (K) ratio in patients with COVID-19 and healthcare workers (HCWs). **e**, Correlation between *KYAT3* (expression averaged within each age group) and age in Genotype-Tissue Expression (GTEx) samples (n = 729 males, 1914 females). Metabolites are displayed as ion intensity log10 transformed, cytokines and chemokines are pg/mL log10 transformed, T cell subsets are % in CD3 T cells, T cell number are 10^6 cells/mL, PBMCs are % in live PBMCs. ![Fig. 3.](http://medrxiv.org/https://www.medrxiv.org/content/medrxiv/early/2020/09/08/2020.09.06.20189159.1/F3.medium.gif) [Fig. 3.](http://medrxiv.org/content/early/2020/09/08/2020.09.06.20189159.1/F3) Fig. 3. Tryptophan metabolites, immune markers and disease severity. **a**, Heatmap of correlations between metabolites in the tryptophan pathway and immune markers by disease severity. Spearman correlations > 0.5 or < –0.5 are displayed with a p< 0.05. **b**, Correlation between kynurenic acid (KA) and immune markers by disease severity. 95% confidence intervals (CIs) for the correlation coefficients are indicated in shaded area colored according to disease progression status. **c**, Comparison of the ratio of KA:kynurenine (K) level by disease severity stratified by sex. Stabilized (females n = 16, males = 11), deteriorated (females n = 6, males n = 6). Nonparametric Kruskal–Wallis rank sum test with pairwise Wilcoxon Mann-Whitney U test, *p* values adjusted for false discovery rate (Benjamini-Hochberg). ***p*< 0.01, NS. not significant. **d**, Correlation between the ratio of KA:kynurenine (K) and CXCL9 and CCL1 stratified by disease severity and sex. Metabolites are displayed as ion intensity log10 transformed, cytokines and chemokines are pg/mL log10 transformed, T cell subsets are % in CD3 T cells, T cell number are 10^6 cells/mL, PBMCs are % in live PBMCs. ### Kynurenic acid is associated with a sex-specific immune response To further understand the sex-specific correlation of KA to immune markers, we examined the relationship between KA and kynurenine (**Fig. 2b**). Kynurenine (K) is a product of tryptophan metabolism that is converted to KA by kynurenine aminotransferases (KATs), which are encoded by *KYAT* genes (**Fig. 2c**). Of note, tryptophan was inversely associated with COVID-19 disease, as noted in **Extended Data Table 3**. We examined the ratio of KA:K in patients with COVID-19 as a surrogate for KAT-mediated production of KA from K7. In males, we observed that a high KA:K was positively correlated with IL6, CCL1, CCL21, TNFα, M-CSF, NK cells, and CD8 terminally differentiated effector memory (Temra) T cells (**Fig. 2b**). A high KA:K was negatively correlated with sCD40L, PDGFAA, PDGFAB/BB, monocytes, PD1+TIM3+CD8 T cells, CD38+HLA-DR+CD8+ T cells, naïve T cells, and IL6+CD4 T cells (**Fig. 2c**). Of note, a high KA:K was positively correlated with T cell activation in females, but overall T cell numbers showed a negative correlation with the ratio of KA:K in males with COVID-19 (**Fig. 2c)**. We also observed that KA:K and KA positively correlated with age in males with COVID-19 (**Fig. 2d, Extended Data Fig. 2a**). KA had a low positive correlation to age in females with COVID-19, but the ratio of KA:K was not correlated (**Extended Data Fig. 2a, Fig. 2d**). In HCWs, KA negatively correlated with age only in males (**Extended Data Fig. 2a**), while KA:K negatively correlated with age in both males and females (**Fig. 2d**). Closer examination of other metabolites involved in K and KA metabolism revealed additional correlates of the cellular immune response during COVID-19. The microbial catabolite of tryptophan, indole-3-lactic acid8, was positively associated with IL4+CD4 and CD38+HLADR+CD8 cells in males (Extended Data Fig. 2b). In females, indole-3-lactic acid was negatively associated with plasma levels of G-CSF, M-CSF, and CXCL10; K was positively associated with IL2, CCL15, CXCL13, and SCF; and tryptophan was positively correlated with follicular CD8 T cells (**Fig. 2b**). To evaluate whether the sex-specific association between KA and the immune response, which was observed in males with COVID-19, is a phenomenon also present in healthy individuals, we analyzed tissue-specific expression data from the Genotype-Tissue Expression (GTEx) Project9. While *KYAT* genes generally tended to have more positive correlations with cytokines in males compared to females, *KYAT3* had particularly stronger correlations in a subset of tissues (including the brain and colon), many of which are classically involved in COVID-19 (Extended Data Fig. 3). Within the brain, we found that these positive correlations with cytokines were specific to older males (aged > 60 years old) (Extended Data Fig. 4a). Because KA is a ligand for AhR which regulates immune responses and inflammation6, we analyzed AhR activation using a previously defined score10 and found that AhR activation correlates most positively with *KYAT3* expression in older males in brain and muscle, while closely correlating in colon (**Extended Data Fig. 4b**). Correlations in the brain became even more pronounced when analyzing only the AhR target gene *CYP1B1*, which is classically used to follow AhR activation in the brain (**Extended Data Fig. 4c**)11. We also show that *KYAT3* expression decreases with age in both males and females, which is consistent with the decreasing ratios of KA:K observed in HCW control samples (**Fig. 2e**). In summary, these data suggest that older males are uniquely sensitive to increases in KA levels, reacting disproportionately with increased levels of inflammatory cytokines, likely as part of a broader AhR activation. ![Fig. 4.](http://medrxiv.org/https://www.medrxiv.org/content/medrxiv/early/2020/09/08/2020.09.06.20189159.1/F4.medium.gif) [Fig. 4.](http://medrxiv.org/content/early/2020/09/08/2020.09.06.20189159.1/F4) Fig. 4. Glutamate, immune markers and disease severity. **a**, Comparison of glutamate levels in stablized patients and deteriorated patients (left panel) and stratified by sex (right panel). Stabilized patients (n = 27), deteriorated patients (n = 12), stabilized females (n = 16), deteriorated females (n = 6), stabilized males (n = 11), and deteriorated males (n = 6). Nonparametric Kruskal–Wallis rank sum test with pairwise Wilcoxon Mann-Whitney U test, *p* values adjusted for false discovery rates (FDR) (Benjamini-Hochberg). ***p*< 0.01, NS. not significant. **b**, Correlation between glutamate and Eotaxin2, IL5, IL6, CD4 T cells, CD4rnTreg cells, CD8 T cells and GzB+CD8 cells in stabilized patients and deteriorated patients. **c**, Correlation between glutamate and immune markers eotaxin2, IL5, IL6, CD4 T cells, CD8 T cells, GzB+CD8 cells, and IL6 in stabilized patients and deteriorated patients stratified by sex. 95% confidence intervals (CIs) for the correlation coefficients were indicated as the shadowed area colored according to progression status. Metabolites are displayed as ion intensity log10 transformed, cytokines and chemokines are pg/mL log10 transformed, T cell subsets are % in CD3 T cells, T cell number are 10^6 cells/mL, PBMCs are % in live PBMCs. ### Kynurenic acid correlates with disease severity in a sex-specific manner As KA levels correlate with numerous immune markers, and these immune markers have been previously linked disease progression2, we examined whether KA was directly associated with disease severity. We used previously defined clinical scores to classify disease severity in COVID-19 patients as deteriorated (males n = 6, females n = 6) or stabilized (males n = 11, females n = 16)2. KA levels were not significantly different between deteriorated and stabilized patients, or after additional stratification by sex. However, KA was positively correlated with CXCL9, IL6, IL12p40, CCL1, CCL3, CCL15, CCL21, CCL27, SCF, M-CSF, and G-CSF in patients that deteriorated with COVID-19 (|R|>0.5 and *p* value< 0.05). In stabilized patients, KA was also positively correlated with CXCL9 and CX3CL1 (**Figs. 3a, b**). We further examined whether KA:K was correlated with disease severity by sex. Males who deteriorated had a significantly higher KA:K compared to those that stabilized, whereas there was no difference in KA:K between females based on clinical course (**Fig. 3c**). Furthermore, a high KA:K was positively correlated with CXCL9 and CCL1 in males that deteriorated, but this correlation was not seen in patients that stabilized, or in females (**Fig. 3d**). We also examined whether any of the 17 metabolites associated with COVID-19 status (Extended Data Table 3) were correlated with disease severity. We observed that only glutamate was associated with disease trajectory, and a significantly higher glutamate level was observed in stabilized patients compared to those that deteriorated (**Fig. 4a)**; incidentally, KA is a glutamate receptor antagonist, thus high KA production could inhibit glutamate release12. Correlation analysis revealed that Eotaxin2, IL5, CD4 T cells, CD4rnTreg negatively correlated with glutamate in deteriorated patients, whereas IL6, CD8 T cells and GzB+CD8 T cells positively correlate (**Fig. 4b)**. Further stratification by sex, showed a similar trend where stabilized patients had higher levels of glutamate than deteriorated patients within each sex group (**Fig. 4a)**. Correlation analysis of immune markers with glutamate by sex revealed a negative correlation to Eotaxin2 and positive correlation to IL6, only in females that deteriorate. CD4 T cells negatively correlate and GzB+CD8 T cells positively correlate with glutamate, only in males that deteriorate. IL5 negatively correlates with glutamate in males who deteriorate whereas it negatively correlates with glutamate in females that stabilize. CD8 T cells positively correlate with glutamate in males that deteriorate and negatively correlate in males that stabilize (**Fig. 4c**). These data suggest that low levels of glutamate may contribute to poorer disease progression in patients with COVID-19 by regulating immune responses. In addition, a high KA:K is correlated to poorer prognosis only in male COVID-19 patients. ### Discussion Patients with severe COVID-19 disease experience a “cytokine storm”, characterized by the elevation of pro-inflammatory cytokines and aggressive inflammatory response13, and sex-specificity in immune response has been previously reported that could underlie differences in clinical outcomes2. Our analysis of serum metabolites from COVID-19 patients, reveals that KA and the ratio of KA:K has a strong relationship to sex-specific immune response and clinical disease course. Importantly, a previous study revealed that serum metabolites in tryptophan and K metabolism correlate with IL6 in a sex-aggregated cohort of COVID-19 patients, but sex-specificity was not examined14. Our study shows that, in males, a high ratio of KA:K was positively correlated with increased levels of cytokines/chemokines and portends clinical deterioration. On the other hand, a negative association was observed with Eotaxin, sCD40L, PDGFs, and T cells, indicating that males with a high KA:K may have a poorer response to inflammation associated with COVID-19, including decreased eosinophil recruitment and T cell activation15,16. A previous study showed that PDGFs associate with better prognosis from COVID-19 if the patients have lower cytokine levels, supporting the association between a higher KA:K and poorer outcomes in males17. In females with COVID-19, a high KA:K positively correlated with a small number of cytokines and also T cell activation, but in contradistinction to males, high KA:K was not associated with disease severity. These results therefore support the role of K metabolism in sex-related differences previously reported in immune responses to COVID-192. Analyzing gene expression data from GTEx, we found that older males (but not females or younger males) appear to have exquisite sensitivity to changes in *KYAT* gene expression (which we used as a proxy for KA levels), whereby natural increases in *KYAT* expression are met with concomitant natural increases in tissue cytokine expression. It is worthwhile to note that the tissues exhibiting these sex-specific correlations – including brain, muscle, kidney, and colon – are those that are commonly implicated in symptoms of COVID-19 patients such as anosmia, myalgia, acute kidney injury, and gastrointestinal distress. Given its role in regulating the immune system and inducing pro-inflammatory cytokines like IL66, modulations in the AhR signaling pathway likely account for this differential response among older males. In support of this hypothesis, we show that AhR activation is most strongly associated with *KYAT3* expression in healthy older males. Also, studies have already shown that male rodents have a more toxic response to stereotypical AhR agonists like TCDD18. Furthermore, testosterone-mediated signaling is known to inhibit AhR activity19, and given the decreasing serum levels of testosterone seen in older males20, it seems plausible that healthy older males could be naturally susceptible to greater AhR activation by endogenous ligands. In the context of COVID-19 infection, patients presented with elevated KAT activity (as suggested by high KA:K ratios), especially among deteriorating male patients. A recent study has demonstrated that a similar AhR induction occurs in the context of murine coronavirus infection, inducing IDO-2 expression17. Notably, two major risk factors for COVID-19, type 2 diabetes and obesity, have already been shown to have increased AhR ligand activity21,22. Such an influx of endogenous AhR ligands, combined with an already elevated susceptibility for AhR activation, therefore, would pose a significantly elevated risk of developing a cytokine storm, specifically in older male patients. The analysis in this study also revealed discrete serum metabolites associated with COVID-19 that may account for some of the varying clinical outcomes in these patients. For instance, metabolites that were positively associated with COVID-19 (Extended Data Table 3) have inflammatory (palmitoleic23 and arachidonic acids24) and neurological (glutamate25 and cysteine-S-sulfate26) roles. Metabolites negatively associated with COVID-19, are involved in the urea cycle and nitric oxide (NO) synthesis pathway (proline, citrulline, and glutamine27). The NO synthesis pathway mediate responses to pro-inflammatory cytokines, macrophages, and neutrophils. Low levels of citrulline have been observed in patients with acute respiratory distress syndrome28 and can cause NO synthase uncoupling and decreased NO synthesis, which is important for vascular function and endothelial cell function29. It was recently suggested that therapeutic NO could be used to improve pulmonary vascular function in COVID-1930. Of the 17 metabolites associated with COVID-19 status, only glutamate was associated with disease severity. In addition, sex-specific correlations between immune responses and glutamate were observed by disease severity. Males that deteriorated from COVID-19 had positive correlations between glutamate and CD8 T cells, and a negative correlation with CD4 T cells. A previous study showed that higher innate immune cytokine levels are associated with clinical deterioration in females with COVID-192, here we observed that IL6 was positively correlated with glutamate only in females that deteriorate. Increasing levels of glutamate have also been shown to decrease IL5 secretion31, and we observed a negative correlation between glutamate and IL5 in males that deteriorate and also in females that stabilize, furthermore, Eotaxin2 negatively correlates with glutamate in females that deteriorate. Similar to Eotaxin2, IL5 is also linked to eosinophil activation, therefore in females, glutamate may be important in regulating eosinophilia in COVID-19. Incidentally, KA is a glutamate receptor antagonist. Glutamate receptors are expressed on the surface of T cells and expression of these receptors is triggered by T cell activation32. Glutamate transporters have also been described in various immune cells32, therefore, this correlation could be therefore reflective of the actions of KA on glutamate levels and also immune cell responses to COVID-1933. Because our study did not analyze non-COVID individuals exhibiting similar clinical symptoms to COVID-19 patients, it remains a possibility that elevated KA and KA:K may lead to elevated cytokines and more broadly mediate the inflammatory symptoms of other pathologies. This possibility, however, does not detract from our observations in COVID-19 patients or the potential of AhR as a therapeutic target in COVID-19. In summary, we have identified serum metabolites associated with COVID-19 clinical course, immune response and sex-specific differences. Among these metabolites, perhaps the most salient discovery is the identification of KA as a metabolite associated with sex, age, increased disease severity, and elevated cytokine and chemokine levels. KA is a ligand for AhR, and when activated, AhR is a master regulator of immune responses and inflammation. Sex-specific agonism of AhR has yet to be reported in humans, but appears to be a prominent feature in COVID-19 disease, potentially underlying the cytokine storm and dampening of T cell activation. In addition, KA is known to dampen glutamate release12, and we observed lower levels of glutamate in patients that deteriorate compared to those that stabilize. Further investigation into the relevance of KA, KAT, and AhR activation in COVID-19 and the role of glutamate in clinical outcomes will be of utmost importance, particularly for understanding the sex-specific differences in immune response and patient outcomes. As we learn more about the impact of the metabolome on COVID-19 disease course, clinicians may find that modulating metabolite levels, either through enteral nutrition or targeted metabolic enzymes may alter disease trajectory. ## Data Availability Untargeted metabolomics data, metabolomics protocols, and code is available on the MetaboLights data repository accession number MTBLS1987 ([https://www.ebi.ac.uk/metabolights/](https://www.ebi.ac.uk/metabolights/)). Clinical and immunological data is available from previous publication2. [https://www.ebi.ac.uk/metabolights/](https://www.ebi.ac.uk/metabolights/) ## Methods ### Clinical biospecimens Serum samples were collected from patients enrolled in the IMPACT study from Cohort A as described and stored at −20°C2. Cohort A consisted of 39 patients aged ≥ 18 years old that tested positive for SARS-CoV-2 by RT-PCR from nasopharyngeal and/or oropharyngeal swabs (females n = 22, males n = 17)34. Intersex individuals were not represented in this study. Prior to the serum collection, these patients were not in an intensive care unit, had not received tocilizumab and had not received high dose corticosteroids. Patients on hydroxychloroquine (n = 29) and remdesivir (n = 3) were not excluded. For control groups, we used 20 serum samples collected from COVID-19 uninfected health care workers working at the Yale-New Haven Hospital between April 2nd and April 28th 2020 who enrolled in the IMPACT study (females n = 10, males n = 10). The detailed demographics and clinical characteristics of these study participants and controls are shown in Extended Data Table 1. ### Immune markers and analysis of disease severity An immune panel of markers for each patient was obtained and published in a previous study2. The patients were assessed with a locally developed clinical scoring system for disease severity17; 1: admitted and observed without supplemental oxygen, 2: required ≤ 3L supplemental oxygen via nasal canal to maintain SpO2 > 92%, 3: received tocilizumab, which per hospital treatment protocol required that the patient to require > 3L supplemental oxygen to maintain SpO2 > 92%, or, required > 2L supplemental oxygen to maintain SpO2 > 92% and had a high sensitivity C-reactive protein (CRP) > 70. 4: the patient required intensive care unit (ICU) level care, 5: the patient required intubation and mechanical ventilation. In relation to the WHO scoring, our clinical score 1, 2/3, 4, 5 largely correspond to WHO score 3, 4, 5, 6/7, respectively35. Detailed demographic information is available from2. For the patients who are 90-year-old or older, their ages were protected health information, and 90 was put as the surrogate value for the analyses. Individuals with active chemotherapy against cancers, pregnant patients, patients with background hematological abnormalities, patients with autoimmune diseases and patients with a history of organ transplantation and on immunosuppressive agents, were excluded from this study. ### Serum metabolite extraction Serum samples (50 μL) were thawed and deactivated for COVID-19 in 150 μL acetone:methanol (50:50 v/v) for 60 min at room temperature. Control samples were treated in the exactly the same manner. To precipitate proteins, the samples were incubated for 2 hours at −20 °C, followed by centrifugation at 13,000 rpm (15,000× g) and 4 °C for 15 min. The resulting supernatant was removed and evaporated to dryness for 12 h using a vacuum concentrator (Thermo Fisher Scientific, Waltham, MA, USA). The dry extracts were then reconstituted in 100 μL of ACN:H2O (1:1, v/v), sonicated for 10 min, and centrifuged at 13,000 rpm (15,000× g) and 4 °C for 15 min, to remove insoluble debris. The supernatants were transferred to ultra performance liquid chromatography (UPLC) autosampler vials (Thermo Scientific, Waltham, MA, USA). A pooled quality control (QC) sample was prepared by mixing 5 μL of extracted solution from each sample into a similar UPLC vial. All the vials were then capped and stored at −80 °C prior to UPLC-mass spectrometry (MS) analysis. ### UPLC-MS-based metabolomics analysis To comprehensively analyze the serum metabolome, both hydrophilic interaction chromatography (HILIC)-MS and reverse phase liquid chromatography (RPLC)-MS approaches were used. A UPLC system (H-Class ACQUITY, Waters Corporation, Milford, MA, USA), coupled to a quadrupole time-of flight (QTOF) (Xevo G2-XS QTOF, Waters Corporation, Milford, MA, USA), was used for MS data acquisition. A Waters ACQUITY UPLC BEH Amide column (particle size, 1.7μm; 100 mm (length) × 2.1 mm (i.d.)) and Waters ACQUITY UPLC BEH C18 column (particle size, 1.7 μm; 100 mm (length) × 2.1 mm (i.d.)) were used for the UPLC-based separation of metabolites. The column temperature was kept at 25 °C for HILIC-MS analysis and 30 °C for RPLC-MS analysis. The solvent flow rate was 0.5 mL/min, and the sample injection volume was 4 μL for HILIC-MS and RPLC in positive mode analysis, 2 μL for HILIC-MS in negative mode, and 6 μL for RPLC-MS negative mode. For HILIC-MS analysis, mobile phase A was 25 mM NH4OH and 25 mM NH4OAc in water, while the mobile phase B was acetonitrile, for both electrospray ionization (ESI), positive and negative mode, respectively. The linear gradient was set as follows: 0∼0.5 min: 95% B; 0.5∼7 min: 95% B to 65% B; 7∼8 min: 65% B to 40% B; 8∼9 min: 40% B; 9∼9.1 min: 40% B to 95% B; 9.1∼12 min: 95% B. For RPLC-MS analysis, the mobile phases A was 0.1% formic acid in H2O, while the mobile phases B was 0.1% formic acid in acetonitrile, respectively, for ESI+. Mobile phase A was 5 mM NH4OAc in H2O, while the mobile phases B was acetonitrile for ESI−. The linear gradient was set as follows: 0∼1 min: 1% B; 1∼8 min: 1% B to 100% B; 8∼10 min: 100% B; 10∼10.1 min: 100% B to 1% B; 10.1∼12 min: 1% B. Pooled samples were analyzed every eight injections during the UPLC-MS analysis to monitor the stability of the data acquisition and used for subsequent data normalization. QTOF scan data (300 ms/scan; mass scan range 50–1000 Da) were initially acquired for each biological sample for metabolite quantification. Then, both DDA (data-dependent acquisition) data (QTOF scan time: 100 ms/scan, MSMS scan time 500 ms/scan, collision energy 20 eV, top 5 most intense ions were selected for fragmentation, exclude former target ions (4 s after 2 occurrences)) and MSE data (low energy scan: 300 ms/scan, collision energy 6 eV; high energy scan: 300 ms/scan, collision energy 20 eV, mass scan range 25–1000 Da) were acquired for QC samples to enable metabolite identification. ESI source parameters on the Xevo GS-XS QTOF were set as the following: capillary voltage 1.8 kV, sampling cone 30 V, source temperature 100 °C, desolvation temperature 550 °C, cone gas flow 40 L/h, desolvation gas flow 900 L/h. ### UPLC-MS data processing The raw MS data (.raw) were converted to mzML files using ProteoWizard MSConvert (version 3.0.6150, [http://www.proteowizard.sourceforge.net/](http://www.proteowizard.sourceforge.net/)). The parameters of min SNR and min peak spacing were set as 0.1 for peak picking in ProteoWizard. The files were then processed in R (version 3.4.3), using the XCMS package for feature detection, retention time correction, and alignment36. The XCMS processing parameters were optimized and set as follows: mass accuracy for peak detection = 20 ppm; peak width c = (2, 30); snthresh = 6; bw = 10; mzwid = 0.015; minfrac = 0.5. The CAMERA package was used for subsequent peak annotation. The resulting data were normalized using the support vector regression algorithm in R, to remove an unwanted system error that occurred among intra- and inter-batches37. Initial metabolite identification was performed using the MetDNA algorithm38. Metabolites were further identified by matching retention time with an in-house metabolite standard library. In addition, metabolite identification was carried out by matching accurate mass and experimental MS/MS data against online databases (METLIN and HMDB). ### Multivariable logistic regression Multivariable logistic regression analyses were performed on the R platform (version 3.4.3) using an R function “glm ()” The model for each metabolite was adjusted for age, BMI and sex to discover metabolites associated with COVID-19 disease. Levels of metabolites are log10 transformed ion intensity. *p* values were adjusted for multiple testing with Benjamini-Hochberg-based FDR using an R function “p.adjust ()”. ### Spearman correlation analysis Spearman correlation analyses were performed on R platform (version 4.0.2) using an R package “psych”. Correlation coefficient R and *p* values were calculated using an R function “corr.test ()” Using previously defined interpretations of correlation coefficients, we used an |R| value of 0.5−1.0 to mark moderate-to-very high correlations39. Heatmaps were plotted using an R package “pheatmap”. ### Chord diagram The chord diagrams were plotted on R platform (version 4.0.2) using an R package “circlize”. Correlations between metabolites and immune responses with R > 0.5 or < −0.5, and *p* value < 0.05 were plotted out. ### Gene expression analysis Gene TPMs, subject phenotypes, and sample attributes data were downloaded from GTEX Portal (gtexportal.org, accession phs000424.v8.p2). After TPM values were transformed as log10(TPM+1), composite expression scores were calculated by adding individual expression values together. Patients who were 60 years or older were coded as “Older,” while patients 30 years or younger were coded as “Young.” After loading the expression data into R with the CePa package, Pearson correlation coefficients were calculated for pairs of target genes within each tissue of each sex, and data was visualized as a heatmap displaying the difference between male and female coefficients using the ComplexHeatmap package. Male-specific correlations were validated by scatter plots and linear regressions, which were generated using the ggplot2 R package. ## Data Availability and Code Availability Untargeted metabolomics data, metabolomics protocols, and code is available on the MetaboLights data repository accession number MTBLS1987 ([https://www.ebi.ac.uk/metabolights/](https://www.ebi.ac.uk/metabolights/)). Clinical and immunological data is available from previous publication2. Data processing R code is available in Supplementary Information. ## Ethical statements This study was approved by Yale Human Research Protection Program Institutional Review Boards (FWA00002571, Protocol ID. 2000027690). Informed consent was obtained from all enrolled patients and healthcare workers. ## Author Contributions All authors contributed to the data discussions and writing of this manuscript. Additional contributions: Y.C carried out the metabolomics mass spectrometry and data analysis, and prepared figures and tables, D.J.K carried out design of the gene expression analysis, performed the gene expression analysis and interpreted the results, T.T was involved in acquiring the biospecimens and experimental design, D.B validated statistical analysis design and reproduced the data independently, S.M. assisted with the biostatistics in the manuscript, N.J.W.R was involved in experimental design and figure design, A.C-M, A.J.M, M.C.M acquired the biospecimens and performed viral inactivation, B.I, J.K, C.L, T.M, J.O, J.S, P.W acquired immune cell data, A.I.K, contributed to data and clinical interpretation, S.A.K contributed to data interpretation and clinical discussions on data, A.I was involved in acquiring the biospecimens and experimental design, securing funds and supervising the study, C.H.J was involved in experiment design, data analysis, drafting the initial manuscript, securing funds and supervising the study. Y.C., D.B., C.H.J., and N.J.W.R are members of the COVID-19 Mass Spectrometry Coalition40. ## Yale IMPACT Research Team authors (in alphabetical order) Tara Alpert, Kelly Anastasio, Michael H. Askenase, Maria Batsu, Santos Bermejo, Sean Bickerton, Anderson Brito, Kristina Brower, Molly L. Bucklin, Staci Cahill, Melissa Campbell, Yiyun Cao, Edward Courchaine, Rupak Datta, Giuseppe Deluliis, Charles Dela Cruz, Rebecca Earnest, Shelli Farhadian, Joseph Fauver, Renata Filler, John Fournier, Bertie Geng, Laura Glick, Nathan Grubaugh, Ryan Handoko, Christina Harden, Cole Jensen, Chaney Kalinich, William Khoury-Hanold, Lynda Knaggs, Maxine Kuang, Eriko Kudo, Sarah Lapidus, Joseph Lim, Melissa Linehan, Feimei Liu, Peiwen Lu, Alice Lu-Culligan, Maksym Minasyan, Amyn A. Malik, Anjelica Martin, Irene Matos, David McDonald, Maura Nakahata, Nida Naushad, Allison Nelson, Jessica Nouws, Angela Nuñez, Marcella Nunez-Smith, Abeer Obaid, Camila Odio, Ji Eun Oh, Saad B. Omer, Isabel M. Ott, Annsea Park, Hong-Jai Park, Xiaohua Peng, Mary Petrone, Sarah Prophet, Harold Rahming. Tyler Rice, Aaron Ring, Kadi-Ann Rose, Lorenzo Sewanan, Lokesh Sharma, Albert Shaw, Denise Shepard, Erin Silva, Mikhail Smolgovsky, Eric Song, Nicole Sonnert, Yvette Strong, Codruta Todeasa, Maria Tokuyama, Jordan Valdez, Sofia Velazquez, Arvind Venkataraman, Pavithra Vijayakumar, Chantal B.F. Vogels, Eric Y. Wang, Annie Watkins, Elizabeth B. White, Anne L. Wyllie, Yexin Yang. ## Competing interest declaration The authors declare no competing financial or non-financial interests. ## Additional Information ### Supplementary Information This file contains Supplementary Table 1. Correlations between metabolites and immune markers in patients with COVID-19. Supplementary Table 2. Correlations between metabolites and immune markers in healthcare workers. The two tables include the calculated spearman correlation coefficient and *p* value. Supplementary Information contains data processing code for use in R. ### Extended Data Figures ![Extended Data Figure 1.](http://medrxiv.org/https://www.medrxiv.org/content/medrxiv/early/2020/09/08/2020.09.06.20189159.1/F5.medium.gif) [Extended Data Figure 1.](http://medrxiv.org/content/early/2020/09/08/2020.09.06.20189159.1/F5) Extended Data Figure 1. Scatter plots of correlation coefficient against −log10 (*p* value) between metabolites and immune markers in all patients with COVID-19, males with COVID-19, and females with COVID-19, respectively. ![Extended Data Figure 2.](http://medrxiv.org/https://www.medrxiv.org/content/medrxiv/early/2020/09/08/2020.09.06.20189159.1/F6.medium.gif) [Extended Data Figure 2.](http://medrxiv.org/content/early/2020/09/08/2020.09.06.20189159.1/F6) Extended Data Figure 2. Correlations between metabolites and immune markers in patients with COVID-19 and healthcare workers stratified by sex. **a**, Correlation between age and kynurenic acid levels in patients with COVID-19 (left) and HCWs (right). **b**, Correlation between indole-3-lactic acid and IL4+CD4, CD38+HLA-DR+CD8, G-CSF, M-CSF and CXCL10 in males with COVID-19 and females with COVID-19, respectively. 95% confidence intervals (CIs) for the correlation coefficients were indicated as the shadowed area colored according to sex. ![Extended Data Figure 3.](http://medrxiv.org/https://www.medrxiv.org/content/medrxiv/early/2020/09/08/2020.09.06.20189159.1/F7.medium.gif) [Extended Data Figure 3.](http://medrxiv.org/content/early/2020/09/08/2020.09.06.20189159.1/F7) Extended Data Figure 3. Correlations between KYAT gene expression and cytokines positively associated with either high KA or KA:K. Pearson correlation coefficients were calculated for gene pairs within the indicated tissue for each sex using GTEx data. Differences in the correlations (RMale–RFemale) are presented as heatmaps, with red indicating a more positive correlation in males and blue indicating a more positive correlation in females (n = 729 males, 1914 females). ![Extended Data Figure 4.](http://medrxiv.org/https://www.medrxiv.org/content/medrxiv/early/2020/09/08/2020.09.06.20189159.1/F8.medium.gif) [Extended Data Figure 4.](http://medrxiv.org/content/early/2020/09/08/2020.09.06.20189159.1/F8) Extended Data Figure 4. Correlations between KYAT3, immune markers, and AhR activation in younger and older individuals, stratified by sex. **a,** Correlations between *KYAT3* expression and *IL6, IL10, CXCL9, TNF*, and *M-CSF* in GTEx brain samples. **b**, Correlations between *KYAT3* and AhR activation score in brain, muscle and colon. **c**, Correlation between *KYAT3* and classic AhR target gene *CYP1B1* in brain. ### Extended Data table titles and footnotes View this table: [Extended Data Table 1.](http://medrxiv.org/content/early/2020/09/08/2020.09.06.20189159.1/T1) Extended Data Table 1. Demographic and clinical characteristics of patients with COVID-19 and healthcare workers View this table: [Extended Data Table 2.](http://medrxiv.org/content/early/2020/09/08/2020.09.06.20189159.1/T2) Extended Data Table 2. Positively identified metabolites from both COVID-19 patients and healthcare workers. View this table: [Extended Data Table 3.](http://medrxiv.org/content/early/2020/09/08/2020.09.06.20189159.1/T3) Extended Data Table 3. Multivariable logistic regression analysis of metabolite levels from healthcare workers and COVID-19 patients. ## Acknowledgements We gratefully acknowledge the study participants for their time and commitment to the study. We thank all members of the clinical team at Yale-New Haven Hospital for their dedication and work which made this study possible. This work was in part supported by a gift from the Yale University Rapid Relief Fund, Women’s Health Research at Yale Pilot Project Program, Fast Grant from Emergent Ventures at the Mercatus Center, Mathers Foundation, the Beatrice Kleinberg Neuwirth Fund, and the Ludwig Family Foundation. IMPACT received support from the Yale COVID-19 Research Resource Fund. A.I. is an Investigator of the Howard Hughes Medical Institute. D.J.K. was a Paul and Daisy Soros Fellow and was supported by a grant from the National Cancer Institute (NCI) of the National Institutes of Health (NIH) (F30CA236466) and by a MSTP training grant from the NIH (T32GM007205, T32GM136651). S.A.K was supported by CTSA Grant Number UL1TR001863 from the National Center for Advancing Translational Science (NCATS), components of the NIH, and NIH roadmap for Medical Research, and the Lampman Research Fund in Yale Surgical Oncology. Support was also provided by the Beatrice Kleinberg Neuwirth Fund, Yale Schools of Public Health and Medicine, and NIH U19 AI08992 awarded to A.I.K. ## Footnotes * Labels in Fig.1 revised * Received September 6, 2020. * Revision received September 8, 2020. * Accepted September 8, 2020. * © 2020, Posted by Cold Spring Harbor Laboratory This pre-print is available under a Creative Commons License (Attribution-NonCommercial-NoDerivs 4.0 International), CC BY-NC-ND 4.0, as described at [http://creativecommons.org/licenses/by-nc-nd/4.0/](http://creativecommons.org/licenses/by-nc-nd/4.0/) ## REFERENCES 1. McPadden, J. et al. Clinical Characteristics and Outcomes for 7,995 Patients with SARS-CoV-2 Infection. medRxiv, doi: 10.1101/2020.07.19.20157305 (2020). [Abstract/FREE Full Text](http://medrxiv.org/lookup/ijlink/YTozOntzOjQ6InBhdGgiO3M6MTQ6Ii9sb29rdXAvaWpsaW5rIjtzOjU6InF1ZXJ5IjthOjQ6e3M6ODoibGlua1R5cGUiO3M6NDoiQUJTVCI7czoxMToiam91cm5hbENvZGUiO3M6NzoibWVkcnhpdiI7czo1OiJyZXNpZCI7czoyMToiMjAyMC4wNy4xOS4yMDE1NzMwNXY0IjtzOjQ6ImF0b20iO3M6NTI6Ii9tZWRyeGl2L2Vhcmx5LzIwMjAvMDkvMDgvMjAyMC4wOS4wNi4yMDE4OTE1OS4xLmF0b20iO31zOjg6ImZyYWdtZW50IjtzOjA6IiI7fQ==) 2. Takahashi, T. et al. Sex differences in immune responses that underlie COVID-19 disease outcomes. Nature, doi: 10.1038/s41586-020-2700-3 (2020). [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1038/s41586-020-2700-3&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=32846427&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2020%2F09%2F08%2F2020.09.06.20189159.1.atom) 3. Ganeshan, K. & Chawla, A. Metabolic regulation of immune responses. Annu Rev Immunol 32, 609-634, doi:10.1146/annurev-immunol-032713-120236 (2014). [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1146/annurev-immunol-032713-120236&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=24655299&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2020%2F09%2F08%2F2020.09.06.20189159.1.atom) 4. O’Neill, L. A., Kishton, R. J. & Rathmell, J. A guide to immunometabolism for immunologists. Nature reviews. Immunology 16, 553-565, doi: 10.1038/nri.2016.70 (2016). [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1038/nri.2016.70&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=27396447&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2020%2F09%2F08%2F2020.09.06.20189159.1.atom) 5. Song, J. W. et al. Omics-Driven Systems Interrogation of Metabolic Dysregulation in COVID-19 Pathogenesis. Cell metabolism, doi: 10.1016/j.cmet.2020.06.016 (2020). [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1016/j.cmet.2020.06.016&link_type=DOI) 6. Wirthgen, E., Hoeflich, A., Rebl, A. & Günther, J. Kynurenic Acid: The Janus-Faced Role of an Immunomodulatory Tryptophan Metabolite and Its Link to Pathological Conditions. Frontiers in Immunology 8, doi: 10.3389/fimmu.2017.01957 (2018). [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.3389/fimmu.2017.01957&link_type=DOI) 7. Liu, H. et al. The Metabolic Factor Kynurenic Acid of Kynurenine Pathway Predicts Major Depressive Disorder. Frontiers in psychiatry 9, 552-552, doi: 10.3389/fpsyt.2018.00552 (2018). [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.3389/fpsyt.2018.00552&link_type=DOI) 8. Gao, J. et al. Impact of the Gut Microbiota on Intestinal Immunity Mediated by Tryptophan Metabolism. Front Cell Infect Microbiol 8, 13-13, doi: 10.3389/fcimb.2018.00013 (2018). [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.3389/fcimb.2018.00013&link_type=DOI) 9. Carithers, L. J. et al. A Novel Approach to High-Quality Postmortem Tissue Procurement: The GTEx Project. Biopreservation and biobanking 13, 311-319, doi: 10.1089/bio.2015.0032 (2015). [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1089/bio.2015.0032&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=26484571&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2020%2F09%2F08%2F2020.09.06.20189159.1.atom) 10. Kim, D. J. et al. Vitamin B12 and folic acid alleviate symptoms of nutritional deficiency by antagonizing aryl hydrocarbon receptor. Proceedings of the National Academy of Sciences 117, 15837-15845, doi: 10.1073/pnas.2006949117 (2020). [Abstract/FREE Full Text](http://medrxiv.org/lookup/ijlink/YTozOntzOjQ6InBhdGgiO3M6MTQ6Ii9sb29rdXAvaWpsaW5rIjtzOjU6InF1ZXJ5IjthOjQ6e3M6ODoibGlua1R5cGUiO3M6NDoiQUJTVCI7czoxMToiam91cm5hbENvZGUiO3M6NDoicG5hcyI7czo1OiJyZXNpZCI7czoxMjoiMTE3LzI3LzE1ODM3IjtzOjQ6ImF0b20iO3M6NTI6Ii9tZWRyeGl2L2Vhcmx5LzIwMjAvMDkvMDgvMjAyMC4wOS4wNi4yMDE4OTE1OS4xLmF0b20iO31zOjg6ImZyYWdtZW50IjtzOjA6IiI7fQ==) 11. Jacob, A. et al. Aryl hydrocarbon receptor-dependent upregulation of Cyp1b1 by TCDD and diesel exhaust particles in rat brain microvessels. Fluids Barriers CNS 8, 23-23, doi: 10.1186/2045-8118-8-23 (2011). [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1186/2045-8118-8-23&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=21867498&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2020%2F09%2F08%2F2020.09.06.20189159.1.atom) 12. Réus, G. Z. et al. Kynurenine pathway dysfunction in the pathophysiology and treatment of depression: Evidences from animal and human studies. J Psychiatr Res 68, 316-328, doi: 10.1016/j.jpsychires.2015.05.007 (2015). [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1016/j.jpsychires.2015.05.007&link_type=DOI) 13. Ragab, D., Salah Eldin, H., Taeimah, M., Khattab, R. & Salem, R. The COVID-19 Cytokine Storm; What We Know So Far. Frontiers in Immunology 11, doi: 10.3389/fimmu.2020.01446 (2020). [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.3389/fimmu.2020.01446&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=http://www.n&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2020%2F09%2F08%2F2020.09.06.20189159.1.atom) 14. Thomas, T. et al. COVID-19 infection alters kynurenine and fatty acid metabolism, correlating with IL-6 levels and renal status. JCI Insight 5, doi: 10.1172/jci.insight.140327 (2020). [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1172/jci.insight.140327&link_type=DOI) 15. Paige, L. in Cytokine Effector Functions in Tissues Vol. 1 Ch. 12, (Academic Press, 2017). 16. Chung, H. W. & Lim, J.-B. Clinical significance of elevated serum soluble CD40 ligand levels as a diagnostic and prognostic tumor marker for pancreatic ductal adenocarcinoma. Journal of translational medicine 12, 102, doi:10.1186/1479-5876-12-102 (2014). [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1186/1479-5876-12-102&link_type=DOI) 17. Lucas, C. et al. Longitudinal analyses reveal immunological misfiring in severe COVID-19. Nature, doi: 10.1038/s41586-020-2588-y (2020). [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1038/s41586-020-2588-y&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=32717743&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2020%2F09%2F08%2F2020.09.06.20189159.1.atom) 18. Enan, E., Overstreet, J. W., Matsumura, F., VandeVoort, C. A. & Lasley, B. L. Gender differences in the mechanism of dioxin toxicity in rodents and in nonhuman primates. Reprod Toxicol 10, 401-411, doi:10.1016/0890-6238(96)83995-5 (1996). [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1016/0890-6238(96)83995-5&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=8888412&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2020%2F09%2F08%2F2020.09.06.20189159.1.atom) [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=A1996VJ22400006&link_type=ISI) 19. Jana, N. R. et al. Cross-talk between 2,3,7,8-tetrachlorodibenzo-p-dioxin and testosterone signal transduction pathways in LNCaP prostate cancer cells. Biochemical and biophysical research communications 256, 462-468, doi: 10.1006/bbrc.1999.0367 (1999). [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1006/bbrc.1999.0367&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=10080920&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2020%2F09%2F08%2F2020.09.06.20189159.1.atom) [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000079499000003&link_type=ISI) 20. Harman, S. M., Metter, E. J., Tobin, J. D., Pearson, J. & Blackman, M. R. Longitudinal effects of aging on serum total and free testosterone levels in healthy men. Baltimore Longitudinal Study of Aging. The Journal of clinical endocrinology and metabolism 86, 724-731, doi: 10.1210/jcem.86.2.7219 (2001). [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1210/jc.86.2.724&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=11158037&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2020%2F09%2F08%2F2020.09.06.20189159.1.atom) [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000167004900047&link_type=ISI) 21. Roh, E. et al. Serum aryl hydrocarbon receptor ligand activity is associated with insulin resistance and resulting type 2 diabetes. Acta diabetologica 52, 489-495, doi: 10.1007/s00592-014-0674-z (2015). [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1007/s00592-014-0674-z&link_type=DOI) 22. Kerley-Hamilton, J. S. et al. Obesity is mediated by differential aryl hydrocarbon receptor signaling in mice fed a Western diet. Environmental health perspectives 120, 1252-1259, doi: 10.1289/ehp.1205003 (2012). [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1289/ehp.1205003&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=22609946&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2020%2F09%2F08%2F2020.09.06.20189159.1.atom) [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000308786900021&link_type=ISI) 23. de Souza, C. O. et al. Palmitoleic Acid has Stronger Anti-Inflammatory Potential in Human Endothelial Cells Compared to Oleic and Palmitic Acids. Mol Nutr Food Res 62, e1800322, doi: 10.1002/mnfr.201800322 (2018). [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1002/mnfr.201800322&link_type=DOI) 24. Samuelsson, B. Arachidonic acid metabolism: role in inflammation. Zeitschrift fur Rheumatologie 50 Suppl 1, 3-6 (1991). [PubMed](http://medrxiv.org/lookup/external-ref?access_num=1907059&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2020%2F09%2F08%2F2020.09.06.20189159.1.atom) 25. Zhou, Y. & Danbolt, N. C. Glutamate as a neurotransmitter in the healthy brain. J Neural Transm (Vienna) 121, 799-817, doi:10.1007/s00702-014-1180-8 (2014). [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1007/s00702-014-1180-8&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=24578174&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2020%2F09%2F08%2F2020.09.06.20189159.1.atom) [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000340679500002&link_type=ISI) 26. Abbas, A.-K. et al. S-Sulfo-Cysteine is an Endogenous Amino Acid in Neonatal Rat Brain but an Unlikely Mediator of Cysteine Neurotoxicity. Neurochemical research 33, 301-307, doi:10.1007/s11064-007-9441-7 (2008). [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1007/s11064-007-9441-7&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=17764028&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2020%2F09%2F08%2F2020.09.06.20189159.1.atom) 27. Wu, G. & Morris, S. M., Jr.. Arginine metabolism: nitric oxide and beyond. The Biochemical journal 336 (Pt 1) 1-17, doi: 10.1042/bj3360001 (1998). [Abstract/FREE Full Text](http://medrxiv.org/lookup/ijlink/YTozOntzOjQ6InBhdGgiO3M6MTQ6Ii9sb29rdXAvaWpsaW5rIjtzOjU6InF1ZXJ5IjthOjQ6e3M6ODoibGlua1R5cGUiO3M6NDoiQUJTVCI7czoxMToiam91cm5hbENvZGUiO3M6MTA6InBwYmlvY2hlbWoiO3M6NToicmVzaWQiO3M6NzoiMzM2LzEvMSI7czo0OiJhdG9tIjtzOjUyOiIvbWVkcnhpdi9lYXJseS8yMDIwLzA5LzA4LzIwMjAuMDkuMDYuMjAxODkxNTkuMS5hdG9tIjt9czo4OiJmcmFnbWVudCI7czowOiIiO30=) 28. Ware, L. B. et al. Low plasma citrulline levels are associated with acute respiratory distress syndrome in patients with severe sepsis. Crit Care 17, R10-R10, doi: 10.1186/cc11934 (2013). [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1186/cc11934&link_type=DOI) 29. Wijnands, K. A., Castermans, T. M., Hommen, M. P., Meesters, D. M. & Poeze, M. Arginine and citrulline and the immune response in sepsis. Nutrients 7, 1426-1463, doi: 10.3390/nu7031426 (2015). [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.3390/nu7031426&link_type=DOI) 30. Adusumilli, N. C., Zhang, D., Friedman, J. M. & Friedman, A. J. Harnessing nitric oxide for preventing, limiting and treating the severe pulmonary consequences of COVID-19. Nitric oxide: biology and chemistry 103, 4-8, doi: 10.1016/j.niox.2020.07.003 (2020). [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1016/j.niox.2020.07.003&link_type=DOI) 31. Cai, Y., Guo, T., Wang, Y. & Du, J. Glutamate Metabolism Regulates Immune Escape of Glioma. Madridge Journal of Immunology 2, 53-57 (2018). 32. Pacheco, R., Gallart, T., Lluis, C. & Franco, R. Role of glutamate on T-cell mediated immunity. J Neuroimmunol 185, 9-19, doi: 10.1016/j.jneuroim.2007.01.003 (2007). [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1016/j.jneuroim.2007.01.003&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=17303252&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2020%2F09%2F08%2F2020.09.06.20189159.1.atom) [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000246243400001&link_type=ISI) 33. 1. M. Levite Ganor, Y. & Levite, M. in Nerve-Driven Immunity (ed M. Levite) (Springer, 2012). 34. Moore, H. M. et al. Biospecimen reporting for improved study quality (BRISQ). Cancer Cytopathol 119, 92-101, doi: 10.1002/cncy.20147 (2011). [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1002/cncy.20147&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=21433001&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2020%2F09%2F08%2F2020.09.06.20189159.1.atom) [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000289636800004&link_type=ISI) 35. WHO R&D Blueprint novel Coronavirus COVID-19 Therapeutic Trial Synopsis, <[https://www.who.int/blueprint/priority-diseases/key-action/COVID-19\_Treatment\_Trial\_Design\_Master\_Protocol\_synopsis\_Final\_18022020.pdf](https://www.who.int/blueprint/priority-diseases/key-action/COVID-19\_Treatment\_Trial\_Design\_Master\_Protocol_synopsis_Final_18022020.pdf).> ( 36. Smith, C. A., Want, E. J., O’Maille, G., Abagyan, R. & Siuzdak, G. XCMS: processing mass spectrometry data for metabolite profiling using nonlinear peak alignment, matching, and identification. Analytical chemistry 78, 779-787, doi: 10.1021/ac051437y (2006). [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1021/ac051437y&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=16448051&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2020%2F09%2F08%2F2020.09.06.20189159.1.atom) 37. Shen, X. et al. Normalization and integration of large-scale metabolomics data using support vector regression. Metabolomics: Official journal of the Metabolomic Society 12, 89, doi:10.1007/s11306-016-1026-5 (2016). [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1007/s11306-016-1026-5&link_type=DOI) 38. Shen, X. et al. Metabolic reaction network-based recursive metabolite annotation for untargeted metabolomics. Nature communications 10, 1516, doi:10.1038/s41467-019-09550-x (2019). [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1038/s41467-019-09550-x&link_type=DOI) 39. Mukaka, M. M. Statistics corner: A guide to appropriate use of correlation coefficient in medical research. Malawi medical journal: the journal of Medical Association of Malawi 24, 69-71 (2012). 40. Struwe, W. et al. The COVID-19 MS Coalition-accelerating diagnostics, prognostics, and treatment. The Lancet 395, 1761-1762, doi: 10.1016/S0140-6736(20)31211-3 (2020). [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1016/s0140-6736(20)31211-3&link_type=DOI) 41. Sumner, L. W. et al. Proposed minimum reporting standards for chemical analysis Chemical Analysis Working Group (CAWG) Metabolomics Standards Initiative (MSI). Metabolomics: Official journal of the Metabolomic Society 3, 211-221, doi:10.1007/s11306-007-0082-2 (2007). [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1007/s11306-007-0082-2&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=24039616&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2020%2F09%2F08%2F2020.09.06.20189159.1.atom) [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000249945500006&link_type=ISI)