Elevation of Serum Cytokine Profiles and Liver Metabolomic Normalization in Early Convalescence of COVID-19 Patients

Coronavirus disease 2019 (COVID-19) has become a global public health concern. We aimed to study the cytokine profile during the convalescent phase and its association with liver functions. We performed a retrospective study to investigate the longitudinal dynamic serum cytokine, liver function, and metabolomic profiles, as well as their potential correlations, from the viral replication phase to early convalescence. Our results demonstrated that liver injury was common. Liver injury was significantly associated with higher levels of interleukin (IL)-6 and IL-10 (p < 0.05). However, alanine aminotransferase levels decreased during the first week after hospital discharge (p < 0.01). In parallel, T-cell and B-cell immune response-stimulating cytokine IL-4, but not IL-2, was significantly elevated (p < 0.05). Furthermore, interferon-γ (IFN-γ) and tumor necrosis factor-α (TFN-α) levels increased, in contrast to the decrease in IL-6 and IL-10 levels; liver function returned to normal. The metabolomic analysis supported active recovery during early convalescence of COVID-19 patients that had distinct metabolic profiles associated with the hepatic tricarboxylic acid cycle, amino acid metabolism, and lipid metabolism. In addition, we identified a metabolomic association of IL-4 with liver repair. Our findings suggest that discharged patients continue to recover from the physiological effects of COVID-19, and the association of IL-4, IL-6, and IL-10 levels with metabolic changes and liver function repair may have important implications for clinical manifestations and treatment of COVID-19.

Introduction

Coronavirus disease 2019 (COVID-19) has rapidly spread worldwide and has become a significant public health challenge. As of June 30, 2020, over 10,000,000 cases have been confirmed, with more than 500,000 deaths in 213 countries, and the numbers continue to rise rapidly.

Numerous studies have reported that COVID-19 severity correlates with serum inflammatory cytokine concentrations (1), and mortality often results from cytokine storm (2). Uncontrolled cytokine storm has also been implicated as a central factor contributing to severe acute respiratory syndrome coronavirus (SARS-CoV), and other severe viral infections (3, 4). Several cytokines including interleukin (IL)-6, IL-10, interferon-γ (IFN-γ), and tumor necrosis factor-α (TNF-α) have been observed to increase dramatically during acute infection in COVID-19 patients (1, 5, 6). However, previous serum cytokine profiling has been focused on acute infection of SARS-CoV-2; thus, very limited information has been reported on cytokine levels and kinetics during the convalescent phase of COVID-19.

Liver enzyme abnormality has also been observed in COVID-19 patients, and it appears to be correlated with disease severity (7). Cytokine storm and direct infection of the liver have been suggested to contribute to liver injury, and in some cases, liver failure in patients with COVID-19 (8, 9). To date, no studies have reported on the repair of liver damage in convalescent patients. However, in response to helminthic infections, the type 2 immune response (10), specifically IL-4-dependent macrophage proliferation and activation, is required to promote repair of both the liver and lung.

Furthermore, there is an immense metabolic demand during liver repair and regeneration. Cytokines have been shown to mediate several metabolic changes via a pathway that is commonly initiated through their regulation of the immune system (11). On the other hand, metabolites are required to regulate the homeostasis of cellular activities in hosts (12), such as resolving inflammation from viral infections (13). Metabolomics of H1N1 influenza virus-infected murine lungs identified that metabolic pathways and association clusters were related to inflammatory cytokines (14). Metabolomics has also been successfully used to identify severe drug-induced liver injury, as well as answer important biological questions (15–17).

Therefore, we have performed a retrospective study designed to investigate the longitudinal dynamic serum cytokine profile, liver functions, and metabolomic profiles in infected patients from the viral replication phase through to convalescence.

Materials and Methods

Patients

From January 19 to March 29, 2020, 102 SARS-CoV-2-infection-diagnosed patients were admitted to the First Affiliated Hospital, Zhejiang University School of Medicine. All patients were confirmed to be SARS-CoV-2 nucleic acid positive by real-time fluorescent RT-PCR. Patients were diagnosed in accordance with the World Health Organization's interim guidelines for COVID-19. Data were collected at our hospital. Our study was approved by the Ethics Committee of the First Affiliated Hospital of Zhejiang University School of Medicine (2020 llT-7). The categorization of mild vs. severe COVID-19 patients was conducted according to COVID-19 Diagnosis and Treatment Guideline (Trial 5th version) (Supplementary Material).

Blood Sampling

Venous whole blood samples were collected from patients during early morning, prior to breakfast. All samples were immediately shipped to the biosafety level 3 laboratory in our hospital. The convalescent blood samples were collected from discharged patients who had met the official hospital discharge criteria, specifically a negative result from two consecutive COVID-19 nucleic acid tests and the disappearance of major clinical signs.

Serum Cytokine Measurement

Serum was separated and stored at −80°C for cytokine detection. The IL-2/IL-4/IL-6/IL-10/TNF-α/IFN-γ cytokine assay kits (Cat No#: 8930960) were provided by Agilent Biosciences (Agilent Technologies, California, USA). The cytokines were detected using flow cytometry (ACEA NovoCyte, Agilent Technologies, California, USA) and analyzed using Novocyte kit software, according to the manufacturer's instructions.

Metabolomic Analysis

A total of 123 blood samples were included in our final analysis, covering both the viral replication (72 samples) and convalescent (41 samples) phases, and 10 plasma samples were collected from healthy participants to serve as controls. The collected samples were centrifuged at 3,000 rpm for 5 min and the supernatant was stored at −80°C before analysis. Metabolites were extracted from plasma, and UPLC-MS/MS analysis was performed using a Waters Acquity UPLC coupled with a Xevo G2-Q-Tof (Waters, Milford, MA, USA) in both positive and negative modes. The obtained raw data were pre-processed using Progenesis QI ver. 2.2 (Nonlinear Dynamic). Metabolites were identified by searching the HMDB library (https://hmdb.ca/spectra/ms/search). Pathway analysis was performed using the MetaboAnalyst 4.0 online tool (http://www.metaboanalyst.ca/). Detailed information is outlined in supporting documents.

Statistical Analyses

Statistical analyses were performed using the Statistical Package for the Social Sciences (SPSS, v. 18.0; SPSS, Inc., Chicago, IL, USA). All tests were two-tailed, and p-values < 0.05 were considered indicative of statistical significance. Student's t-test or the nonparametric Mann–Whitney U-test, as appropriate, was used for comparisons of continuous data. Categorical data were compared by Fisher's exact test and correlation analysis was calculated by Pearson correlation or Spearman correlation test.

Results

Demographic and Clinical Characteristics of COVID-19 Patients

We conducted a retrospective study of all patients (102) with confirmed SARS-CoV-2 infection in our hospital. Two patients were subsequently excluded from the final analysis because they were discharged from the hospital the next day. The remaining 100 patients were included in the analysis. In total, 123 blood samples was involved in metabolomics analysis including both discovery group (N = 110) and validation group (N = 13) (Figure 1). The median age of the patients was 54.28 years. Most of the infected patients were men (61; 61%); less than half had underlying diseases (48; 48%), including diabetes (15; 15%), hypertension (35; 35%), cardiovascular disease (8; 8%), pulmonary disease (4; 4%), fatty liver (4; 4%), chronic kidney disease (5; 5%), and HBV (4; 4%). On admission, most patients had fever (83; 83%), cough (69; 69%), and other common symptoms including phlegm (40; 40%), chest distress (25; 25%), myalgia (18; 18%), and fatigue (14; 14%). There were only a small number of patients presenting with headache (9; 9%), dizziness (7; 7%), diarrhea (8; 8%), and nausea or vomiting (7; 7%) (Table 1).

FIGURE 1

Figure 1. Design of the metabolomics study. This study, involving 123 plasma samples, included both discovery group and validation group. Through the analysis of the data from the discovery set, we screened and identified metabolites that were closely related to COVID-19. To validate the reliability of results from discovery set, an independent set of 13 plasma samples was collected and all samples were correctly identified by a random forest machine learning model based on the identified metabolites (Supplementary Table 9).

TABLE 1

Table 1. Demographics and baseline characteristics of patients infected with COVID-19.

Liver Injury and Recovery in COVID-19 Patients

To investigate whether liver enzyme abnormality occurred in this COVID-19 patient cohort, multiple liver function indicators, including serum alanine aminotransferase (ALT), aspartate aminotransferase (AST), alkaline phosphatase (ALP), and total bilirubin (TBiL), were collected from the onset of the disease. Among the 100 COVID-19 patients in our hospital, 62 (62%) had increased levels of ALT, AST, ALP, and TBiL (Table 1 and Supplementary Table 3). These results suggest that liver injury is a common clinical feature of SARS-CoV-2 infection. To determine the course of liver recovery, liver function indicators were measured weekly after patients were discharged from hospital. Serum ALT levels decreased within the first week after discharge (p < 0.01) (Figure 2). Levels of ALT continued to decline significantly through the second (p < 0.05) and fourth weeks (p < 0.001) (Figure 2). Significant changes to AST, ALP, and TBiL levels also occurred toward the normal range (Figure 2), indicating the active recovery of liver functions during early convalescence of COVID-19.

FIGURE 2

Figure 2. The dynamic changes of ALT, AST, ALP, and TBiL in virus replication phase and early convalescence. Levels of significance: p-value: ***, <0.001; **, <0.01; *, <0.05.

Association of Elevated IL-6 and IL-10 With Liver Injury

Cytokine analysis revealed that cytokines IL-6 and IL-10 were elevated in COVID-19 patients on admission and continued to increase significantly as the disease progressed (Figure 3), consistent with previous reports (1). IL-6 and IL-10 levels gradually declined into normal ranges as the viral RNA became undetectable (p < 0.05) (Figure 3). Furthermore, as shown in Supplementary Table 5, IL-6 showed sustained higher levels in the liver impaired group, compared with the normal liver function group. Elevated levels of IL-6 became significant 7–9 days after disease onset (p = 0.034) and became even more pronounced as the disease progressed (p = 0.001 at days 13–15 and ≥16) (Supplementary Table 5). These results suggest a possible correlation between liver damage and the inflammatory responses induced by SARS-CoV-2 infection. In addition, increased levels of IL-6 and IL-10 were also significantly correlated with disease severity and patient age (Supplementary Table 5).

FIGURE 3

Figure 3. The dynamic changes of IL-6, IL-10, IL-4, TNF-α, IFN-γ, and IL-2 in virus replication phase and early convalescence. Levels of significance: p-value: ***, <0.001; **, <0.01; *, <0.05.

Elevation of Serum IL-4 Level in Early Convalescence

Cytokine profiles in convalescent patients remain less characterized; thus, we performed serum cytokine analysis weekly after patients were discharged from hospital. Significantly, IL-4 was elevated in convalescent patients within the first week after discharge (p = 0.036) and continued to stay high at week 2 (p = 0.022) (Figure 3), even though IL-4 had no significant fluctuation through the viral replication phase as previously reported (1). IL-4 plays a key role in promoting naïve T cells to develop into Th2-like cells (18) and immunoglobulin class-switching from IgM to IgE and IgG (19). The observed elevation of IL-4 might suggest that there is a differential set of regulatory T cell and B cell immune responses occurring in early convalescence of COVID-19 patients.

Antiviral cytokines IFN-γ and TNF-α continued to stay at high levels as patients transitioned from the viral replication phase to convalescence, whereas IL-6 and IL-10 gradually declined and normalized (Figure 3). In particular, the serum levels of both IFN-γ and TNF-α showed meaningful increases above normal ranges from weeks 1 to 2 after hospital discharge, suggesting an additional boost of T cell immune responses (IFN-γ: 10.98 to 14.33 pg/ml; TNF-α: 23.27 to 37.96 pg/ml; Supplementary Table 4). In contrast, IL-2 levels remained within normal ranges throughout the course of the study (Figure 3).

Plasma Metabolomic Alternations Associated With Viral Replication and Convalescent Phases

Metabolites were extracted from plasma, and UPLC-QTOF MS analysis was performed for metabolomic analysis. From a total of 2,077 metabolic features detected in metabolomic analysis, 243 were selected based on the following selection criteria: absolute log2 FC > 1; p-value < 0.05; projection variable VIP > 1; the area of ROC curve > 0.8. Among them, in comparison with replication phase, 50 metabolites showed differential profiles in convalescent patients, mostly with upregulation back toward the normal levels (Figures 4A–C), demonstrating that there was a significant change in metabolites during recovery from COVID-19. Using the results obtained in the viral replication phase/convalescence comparison as a reference, we ranked metabolites from the highest fold change value to the lowest value and drew a Cleveland plot showing fold changes of each metabolite across the three pairwise comparisons (Figure 4D). Among the three groups, the metabolites in the viral replication phase notably decreased.

FIGURE 4

Figure 4. Metabolomic profiling of COVID-19 patient plasma. (A) Principal components analysis of the metabolites across the three study groups based on the data of ESI+ and ESI– modes. Colors display the subjects from different groups (blue, normal controls; yellow, virus replication group; green, convalescent group; and red, QC samples). (B) OPLS-DA plot of COVID-19 patients in virus replication phase and early convalescence. (C) Heatmap of 50 selected metabolites, the colors changing from red (upregulated metabolites) to blue (downregulated metabolites). (D) Cleveland plots. Fifty metabolites are ranked according to their fold changes in the virus replication group vs. convalescent group (the highest fold changes on the top; the lowest on the bottom). Fold changes in the virus replication group vs. normal controls and the convalescent group vs. normal controls are also shown. (E) MetPA. Pathway analysis was performed using MetaboAnalyst 4.0 for the 50 differential metabolites. The color and size of each circle are based on its p-value and pathway impact value, respectively.

As shown in Figure 4E, pathway analysis identified changes in 18 pathways during convalescence (Supplementary Table 8). Three pathways, namely, linoleic acid metabolism, biosynthesis of unsaturated fatty acids, and sphingolipid metabolism, were significantly altered. They included 32 downregulated metabolites, such as triglycerides, decanoylcarnitine, and diglycerides, and 43 upregulated metabolites, such as sphingosine 1-phosphate and cholic acid (Supplementary Figure 7). Notably, many metabolites returned to the levels observed in healthy controls (Supplementary Figure 6). Together, the metabolomic analysis indicated that discharged patients continued to recover from the physiological impacts of COVID-19.

Metabolomic Association of Inflammatory Cytokines With Liver Functions

Since liver injury from SARS-CoV-2 was associated with the extent of cytokine expression (Supplementary Table 6), we examined the interaction between metabolic changes of metabolites and cytokines on liver functions. As shown in Figure 5 and Supplementary Figures 4, 5, most of the upregulated metabolites in the convalescent stage were negatively correlated with IL-6, IL-10, ALT, and TBiL, while downregulated metabolites were positively correlated with IL-4. There were 28 metabolites, including amino acids, glycerophospholipids, ceramides, and unsaturated fatty acids, which were significantly correlated with cytokines (IL-4, IL-6, and IL-10) (p < 0.05). In addition, 10 metabolites, including amino acids, glycerophospholipids, steroids, and steroid derivatives, were significantly correlated with liver function indicators (ALT and TBiL) (p < 0.05) (Supplementary Figures 4, 5).

FIGURE 5

Figure 5. Change of expression level metabolites related to the cytokines.

Discussion

In this study, we aimed to analyze the liver function repair and cytokine profiles during convalescence of COVID-19, which has barely been characterized. Cytokine storm and direct infection of the liver has been suggested to contribute to liver injury, and in some cases, liver failure in patients with COVID-19 (8, 9). Our results supported the notion that liver injury is a common clinical feature of SARS-CoV-2 infection, and we observed the active recovery of liver functions during the early convalescent stage of COVID-19. Importantly, serum ALT level decreased within the first week after discharge (p < 0.01) and continued to significantly decline through the second (p < 0.05) and fourth weeks (p < 0.001) (Figure 2). Significant shifts of AST, ALP, and TBiL levels toward the normal range also occurred (Figure 2), indicating the active recovery of liver functions during early convalescence.

Pro-inflammatory cytokines, such as TNF-α and IL-6, are primarily involved in the promotion of inflammatory processes and have an important role in liver injury (20). Cytokine profile analysis confirmed that both IL-6 and IL-10 became significantly elevated (p < 0.05) during the early viral replication phase and then gradually declined into the normal range as viral RNA became undetectable (Figure 3). IL-6 and IL-10 levels were higher in patients with abnormal liver function than in those with normal liver function (Supplementary Table 5). Furthermore, after discharge from hospital, liver function normalized during the first 2 weeks of convalescence, and both inflammation-regulating cytokines continued to decline significantly (p < 0.05), indicating clinical improvement.

We found that immune response-stimulating cytokine IL-4, but not IL-2, was significantly elevated during the first 2 weeks of convalescence. This was a surprising observation, as no significant fluctuation of IL-4 was seen during the viral replication phase (1), which was also confirmed in our study. IL-4 is known to regulate a variety of immune responses, including differentiation of naïve T cell into Th2 cells, and immunoglobulin class switching to IgG1 and IgE in B cells. For example, IL-4 could signal through IL-4Rα to trigger specialized macrophage activation, promoting the mitigation of helminthic infection and tissue repair in the liver and lung (21), or reduce the production of C-reactive protein (CRP) by human primary hepatocytes (22). Furthermore, IL-4 polymorphism has also been associated with an increased risk of liver disease (23) and severe respiratory syncytial virus (RSV) infection (24). All things considered, it is reasonable to postulate that an increase of IL-4 levels may play an important role in the convalescence of COVID-19 through either T cells, B cells, or other type 2 immunity-associated cells, such as macrophages.

Furthermore, abnormal liver function became normalized during the first 2 weeks of convalescence (Figure 2). Additionally, the antiviral cytokines IFN-γ and TNF-α, in contrast to the decline of IL-6 and IL-10, were observed to significantly increase during convalescence of those patients (Figure 3 and Supplementary Table 4). Both IL-4 and IFN-γ were shown to downregulate the expression of SARS coronavirus receptor angiotensin-converting enzyme 2 (ACE2), in order to inhibit viral infection (25). Together, these clinical features further support the notion that type 2 immunity may contribute to liver and lung repair following injury by SARS-CoV-2. Finally, IL-13 has already been shown to share many biological functions with IL-4 and can inhibit IL-6 production through peripheral blood mononuclear cells (26). Therefore, it is recommended to include IL-13 in future cytokine profile analysis (27).

To further corroborate the association of liver injury and repair with cytokine profile changes, metabolic changes were determined in comparison with healthy controls. Metabolomic analysis shows significantly distinct profiles of metabolites and cytokines between the viral replication and the convalescent phases, including the liver-associated amino acid, TCA cycle, steroid hormone biosynthesis, and lipid metabolism (Figure 6). Remarkably, all of these metabolic changes appear to support liver repair. For example, during convalescence, the saturated fatty acid palmitic acid decreased, potentially mitigating the apoptosis of hepatocytes (28). In contrast, the unsaturated fatty acids docosapentaenoic acid (DPA) and docosahexaenoic acid (DHA) increased, possibly promoting liver repair by inhibiting ALT (29). An increase of tryptophan may reverse liver injury by maintaining protein synthesis activity (30). Importantly, during the preparation of our manuscript, part of the metabolomic profile we identified in the viral replication phase was also observed in other studies on SARS-CoV-2 (31, 32). These results are consistent with the general notion that cellular metabolites play key roles in resolving inflammation resulting from various viral infections, including H1N1 influenza (13, 14). Together, these metabolic changes strongly indicate that the liver is undergoing active recovery.

FIGURE 6

Figure 6. Diagram to illustrate potential relationships among liver injury repair, inflammatory cytokine reactions, and metabolic changes in early convalescence. The red color represents upregulated metabolites and blue represents downregulated metabolites. The black arrow represents the changes of cytokines in early convalescence compared to virus replication phase. Lyso PC, lysophosphatidylcholines; PC, phosphatidylcholine; DHEAS, dehydroepiandrosterone sulfate; DHA, docosahexaenoic acid; DPA, docosapentaenoic acid; PA, phosphatidic acid; PI: phosphatidylinositol.

The immune system is particularly sensitive to metabolite availability (21). On the other hand, cytokines have been shown to mediate several metabolic changes via a pathway that is commonly initiated through their regulation of the immune system (11). Correlation analysis identified that many metabolites were associated with proinflammatory IL-6 and anti-inflammatory IL-10 cytokine levels, suggesting a potential cytokine-mediated metabolic dysfunction from COVID-19 pathogenesis (Figure 6). IL-6 can signal through three main pathways: JAK-STAT3, SHP-2-MAPK, and PI3K-AKT. At least 14 upregulated metabolites, such as oleic acid and DHA, and 3 downregulated metabolites, such as palmitic acid and lactosylceramide, were revealed, each of which are known to inhibit these three pathways (Figure 6 and Supplementary Material). IL-10 is a major anti-inflammatory cytokine secreted by macrophages, and it exerts its effects via the JAK-STAT3 pathway (33). Palmitic acid and DHA levels were downregulated and upregulated, respectively, both of which may inhibit STAT3 signaling to promote liver repair. IL-4 is one of the best-known anti-inflammatory cytokines, mediating its biological roles predominantly via the JAK-STAT6 pathway (34). Our metabolic analysis identified four upregulated metabolites, including sphingosine 1-phosphate and DHA, and two downregulated metabolites, palmitic acid and lactosylceramide, which could stimulate T cells to produce IL-4 in early convalescence (Figure 6). The diagram in Figure 6 summarizes the potential interactions among liver injury repair, inflammatory cytokine reactions, and metabolic changes.

Conclusion

In conclusion, by analyzing liver function repair and cytokine profiles in early convalescence of COVID-19, in tandem with the viral replication phase, we have identified that liver injury is a common clinical feature in COVID-19 patients, and it is associated with the increase of cytokine IL-6 and IL-10 levels. Importantly, serum levels of IL-4 were significantly elevated during early convalescence, suggesting a potentially important role of Th2 immune response in liver injury repair. The correlation of liver injury and repair with cytokines was further corroborated by metabolomic analysis, which identified a series of related biomarkers for the recovery of COVID-19 patients. Collectively, our new findings may have important implications in analysis of clinical manifestations and the potential therapeutic treatment of COVID-19.

Data Availability Statement

The original contributions presented in the study are included in the article/Supplementary Material, further inquiries can be directed to the corresponding author/s.

Ethics Statement

The studies involving human participants were reviewed and approved by the Ethics Committee of the First Affiliated Hospital of Zhejiang University School of Medicine (2020 llT-7). The patients/participants provided their written informed consent to participate in this study.

Author Contributions

YL wrote the manuscript and designed the overall scheme. XHe analyzed the medical record information of these patients. MD and LL were responsible for the collection and interpretation of data. XHu, XY, and LH were mainly responsible for the metabolomic analysis by UPLC-MS/MS. YH and QZ participated in cytokine measurement and analysis. JW and LZ were responsible for the collection of the blood samples and shipped the samples to the biosafety level 3 laboratory. XL designed the retrospective study. YQ was the general manager of the project and designed the research. All authors have read and approved the final version of this manuscript.

Funding

This work was supported by the pharmaceutical innovation research funding project of China Pharmaceutical Association—Servier Youth Hospital, Zhejiang Provincial Science and Technology Department Key R & D Plan Emergency Project (Nos. 2020c03123-8, LGF18H310002, and 2020C03046) and the National Natural Science Foundation of China (Nos. 81872935, 81971982, and 81703612). The funding body was not involved in study design, collection/analysis/interpretation of data, and manuscript writing.

Conflict of Interest

The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.

Acknowledgments

We acknowledge Dr. Guofeng Cheng for his technical support to this study. Dr. Cheng was given the opportunity to review the manuscript for medical and scientific accuracy.

Supplementary Material

The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fmed.2021.626633/full#supplementary-material

References

1. Liu J, Li S, Liu J, Liang B, Wang X, Wang H, et al. Longitudinal characteristics of lymphocyte responses and cytokine profiles in the peripheral blood of SARS-CoV-2 infected patients. EBioMedicine. (2020) 55:102763. doi: 10.1016/j.ebiom.2020.102763

PubMed Abstract | CrossRef Full Text | Google Scholar

2. Mehta P, McAuley DF, Brown M, Sanchez E, Tattersall RS, Manson JJ, et al. COVID-19: consider cytokine storm syndromes and immunosuppression. Lancet. (2020) 395:1033–4. doi: 10.1016/S0140-6736(20)30628-0

PubMed Abstract | CrossRef Full Text | Google Scholar

3. Chousterman BG, Swirski FK, Weber GF. Cytokine storm and sepsis disease pathogenesis. Semin Immunopathol. (2017) 39:517–28. doi: 10.1007/s00281-017-0639-8

PubMed Abstract | CrossRef Full Text | Google Scholar

4. Wauquier N, Becquart P, Padilla C, Baize S, Leroy EM. Human fatal zaire ebola virus infection is associated with an aberrant innate immunity and with massive lymphocyte apoptosis. PLoS Negl Trop Dis. (2010) 4:e837. doi: 10.1371/journal.pntd.0000837

PubMed Abstract | CrossRef Full Text | Google Scholar

5. Chen G, Wu D, Guo W, Cao Y, Huang D, Wang H, et al. Clinical and immunological features of severe and moderate coronavirus disease 2019. J Clin Investig. (2020) 130:2620–9. doi: 10.1172/JCI137244

PubMed Abstract | CrossRef Full Text | Google Scholar

6. Huang C, Wang Y, Li X, Ren L, Zhao J, Hu Y, et al. Clinical features of patients infected with 2019 novel coronavirus in Wuhan, China. Lancet. (2020) 395:497–506. doi: 10.1016/S0140-6736(20)30183-5

PubMed Abstract | CrossRef Full Text | Google Scholar

7. Sun J, Aghemo A, Forner A, Valenti L. COVID-19 and liver disease. Liver Int. (2020) 40:1278–81. doi: 10.1111/liv.14470

CrossRef Full Text | Google Scholar

8. Li J, Fan JG. Characteristics and mechanism of liver injury in 2019 coronavirus disease. J Clin Transl Hepatol. (2020) 8:13–7. doi: 10.14218/JCTH.2020.00019

PubMed Abstract | CrossRef Full Text | Google Scholar

9. Wang Y, Liu S, Liu H, Li W, Lin F, Jiang L, et al. SARS-CoV-2 infection of the liver directly contributes to hepatic impairment in patients with COVID-19. J Hepatol. (2020) 73:807–16. doi: 10.1016/j.jhep.2020.05.002

PubMed Abstract | CrossRef Full Text | Google Scholar

10. Minutti CM, Jackson-Jones LH, Garcia-Fojeda B, Knipper JA, Sutherland TE, Logan N, et al. Local amplifiers of IL-4Ralpha-mediated macrophage activation promote repair in lung and liver. Science. (2017) 356:1076–80. doi: 10.1126/science.aaj2067

PubMed Abstract | CrossRef Full Text | Google Scholar

11. Newell MK, Villalobos-Menuey E, Schweitzer SC, Harper ME, Camley RE. Cellular metabolism as a basis for immune privilege. J Immune Based Ther Vaccines. (2006) 4:1. doi: 10.1186/1476-8518-4-1

PubMed Abstract | CrossRef Full Text | Google Scholar

12. Wang YP, Lei QY. Metabolite sensing and signaling in cell metabolism. Signal Transduct Target Ther. (2018) 3:30. doi: 10.1038/s41392-018-0024-7

PubMed Abstract | CrossRef Full Text | Google Scholar

13. Sanchez EL, Lagunoff M. Viral activation of cellular metabolism. Virology. (2015) 479–480:609–18. doi: 10.1016/j.virol.2015.02.038

CrossRef Full Text | Google Scholar

14. Chandler JD, Hu X, Ko EJ, Park S, Lee YT, Orr M, et al. Metabolic pathways of lung inflammation revealed by high-resolution metabolomics (HRM) of H1N1 influenza virus infection in mice. Am J Physiol Regul, Integr Comp Physiol. (2016) 311:R906–R16. doi: 10.1152/ajpregu.00298.2016

PubMed Abstract | CrossRef Full Text | Google Scholar

15. Li S, Sullivan NL, Rouphael N, Yu T, Banton S, Maddur MS, et al. Metabolic phenotypes of response to vaccination in humans. Cell. (2017) 169:862–77. e17. doi: 10.1016/j.cell.2017.04.026

PubMed Abstract | CrossRef Full Text | Google Scholar

16. Patin F, Baranek T, Vourc'h P, Nadal-Desbarats L, Goossens JF, Marouillat S, et al. Combined Metabolomics and Transcriptomics Approaches to Assess the IL-6 Blockade as a Therapeutic of ALS: deleterious alteration of lipid metabolism. Neurotherapeutics. (2016) 13:905–17. doi: 10.1007/s13311-016-0461-3

PubMed Abstract | CrossRef Full Text | Google Scholar

17. Xie Z, Chen E, Ouyang X, Xu X, Ma S, Ji F, et al. Metabolomics and Cytokine Analysis for Identification of Severe Drug-Induced Liver Injury. J Proteome Res. (2019) 18:2514–24. doi: 10.1021/acs.jproteome.9b00047

PubMed Abstract | CrossRef Full Text | Google Scholar

18. Zhao J, Jiang L, Deng L, Xu W, Cao Y, Chen C, et al. Important roles of CD32 in promoting suppression of IL-4 induced immune responses by a novel anti-IL-4Ralpha therapeutic antibody. MAbs. (2019) 11:837–47. doi: 10.1080/19420862.2019.1601985

PubMed Abstract | CrossRef Full Text | Google Scholar

19. Guo B, Rothstein TL. IL-4 upregulates Igalpha and Igbeta protein, resulting in augmented IgM maturation and B cell receptor-triggered B cell activation. J Immunol. (2013) 191:670–7. doi: 10.4049/jimmunol.1203211

PubMed Abstract | CrossRef Full Text | Google Scholar

20. Ge X, Feng Z, Xu T, Wu B, Chen H, Xu F, et al. A novel imidazopyridine derivative, X22, attenuates sepsis-induced lung and liver injury by inhibiting the inflammatory response in vitro and in vivo. Drug Des Dev Ther. (2016) 10:1947–59. doi: 10.2147/DDDT.S101449

PubMed Abstract | CrossRef Full Text | Google Scholar

21. Gieseck RL 3rd, Wilson MS, Wynn TA. Type 2 immunity in tissue repair and fibrosis. Nat Rev Immunol. (2018) 18:62–76. doi: 10.1038/nri.2017.90

PubMed Abstract | CrossRef Full Text | Google Scholar

22. Gabay C, Porter B, Guenette D, Billir B, Arend WP. Interleukin-4 (IL-4) and IL-13 enhance the effect of IL-1beta on production of IL-1 receptor antagonist by human primary hepatocytes and hepatoma HepG2 cells: differential effect on C-reactive protein production. Blood. (1999) 93:1299–307. doi: 10.1182/blood.V93.4.1299.404k26_1299_1307

PubMed Abstract | CrossRef Full Text | Google Scholar

23. Wu Z, Qin W, Zeng J, Huang C, Lu Y, Li S. Association between IL-4 polymorphisms and risk of liver disease: an updated meta-analysis. Medicine. (2015) 94:e1435. doi: 10.1097/MD.0000000000001435

PubMed Abstract | CrossRef Full Text | Google Scholar

24. Choi EH, Lee HJ, Yoo T, Chanock SJ. A common haplotype of interleukin-4 gene IL4 is associated with severe respiratory syncytial virus disease in Korean children. J Infect Dis. (2002) 186:1207–11. doi: 10.1086/344310

PubMed Abstract | CrossRef Full Text | Google Scholar

25. de Lang A, Osterhaus AD, Haagmans BL. Interferon-gamma and interleukin-4 downregulate expression of the SARS coronavirus receptor ACE2 in Vero E6 cells. Virology. (2006) 353:474–81. doi: 10.1016/j.virol.2006.06.011

PubMed Abstract | CrossRef Full Text | Google Scholar

26. Minty A, Chalon P, Derocq JM, Dumont X, Guillemot JC, Kaghad M, et al. Interleukin-13 is a new human lymphokine regulating inflammatory and immune responses. Nature. (1993) 362:248–50. doi: 10.1038/362248a0

PubMed Abstract | CrossRef Full Text | Google Scholar

27. Aversa G, Punnonen J, Cocks BG, de Waal Malefyt R, Vega F Jr, Zurawski SM, et al. An interleukin 4 (IL-4) mutant protein inhibits both IL-4 or IL-13-induced human immunoglobulin G4 (IgG4) and IgE synthesis and B cell proliferation: support for a common component shared by IL-4 and IL-13 receptors. J Exp Med. (1993) 178:2213–8. doi: 10.1084/jem.178.6.2213

PubMed Abstract | CrossRef Full Text | Google Scholar

28. Martinez L, Torres S, Baulies A, Alarcon-Vila C, Elena M, Fabrias G, et al. Myristic acid potentiates palmitic acid-induced lipotoxicity and steatohepatitis associated with lipodystrophy by sustaning de novo ceramide synthesis. Oncotarget. (2015) 6:41479–96. doi: 10.18632/oncotarget.6286

PubMed Abstract | CrossRef Full Text | Google Scholar

29. Guo XF, Sinclair AJ, Kaur G, Li D. Differential effects of EPA, DPA and DHA on cardio-metabolic risk factors in high-fat diet fed mice. Prostaglandins Leukot Essent Fatty Acids. (2018) 136:47–55. doi: 10.1016/j.plefa.2017.09.011

PubMed Abstract | CrossRef Full Text | Google Scholar

30. Ohta Y, Sahashi D, Sasaki E, Ishiguro I. Alleviation of carbon tetrachloride-induced chronic liver injury and related dysfunction by L-tryptophan in rats. Ann Clin Biochem. (1999) 36 (Pt 4):504–10. doi: 10.1177/000456329903600415

PubMed Abstract | CrossRef Full Text | Google Scholar

31. Shen B, Yi X, Sun Y, Bi X, Du J, Zhang C, et al. Proteomic and metabolomic characterization of COVID-19 patient sera. Cell. (2020) 182:59–72. e15. doi: 10.1016/j.cell.2020.05.032

PubMed Abstract | CrossRef Full Text | Google Scholar

32. Wu D ST, Yang X, Song JX, Zhang M, Yao C, Zhou X, et al. Plasma metabolomic and lipidomic alterations associated with COVID-19. Natl Sci Rev. (2020) 7:1157–168. doi: 10.1101/2020.04.05.20053819

CrossRef Full Text | Google Scholar

33. Dumoutier L, Renauld JC. Viral and cellular interleukin-10 (IL-10)-related cytokines: from structures to functions. Eur Cytokine Netw. (2002) 13:5–15. doi: 10.1016/B978-044450973-4/50135-9

PubMed Abstract | CrossRef Full Text | Google Scholar

34. Busch-Dienstfertig M, Gonzalez-Rodriguez S. IL-4, JAK-STAT signaling, and pain. Jak-stat. (2013) 2:e27638. doi: 10.4161/jkst.27638

CrossRef Full Text | Google Scholar