Skip to main content

ORIGINAL RESEARCH article

Front. Immunol., 07 June 2021
Sec. Inflammation

Cytokine/Chemokine Expression Is Closely Associated Disease Severity of Human Adenovirus Infections in Immunocompetent Adults and Predicts Disease Progression

Jin LiJin Li1Jinli WeiJinli Wei2Zhixiang XuZhixiang Xu2Chunmei JiangChunmei Jiang3Mianhuan LiMianhuan Li2Jie ChenJie Chen4Yanjie LiYanjie Li2Minghui YangMinghui Yang1Yuchen GuYuchen Gu2Fuxiang WangFuxiang Wang2Yuelong Shu&#x;Yuelong Shu1†Yang Yang&#x;Yang Yang2†Litao Sun&#x;Litao Sun1†Yingxia Liu*&#x;Yingxia Liu2*†
  • 1School of Public Health (Shenzhen), Sun Yat-sen University, Shenzhen, China
  • 2Shenzhen Key Laboratory of Pathogen and Immunity, National Clinical Research Center for Infectious Disease, State Key Discipline of Infectious Disease, Shenzhen Third People’s Hospital, Second Hospital Affiliated to Southern University of Science and Technology, Shenzhen, China
  • 3Department of Infectious Disease, The People’s Hospital of Longhua, Shenzhen, China
  • 4Research and Development Department, Guangzhou Sagene Biotech Co., Ltd., Guangzhou, China

Increasing human Adenovirus (HAdV) infections complicated with acute respiratory distress syndrome (ARDS) even fatal outcome were reported in immunocompetent adolescent and adult patients. Here, we characterized the cytokine/chemokine expression profiles of immunocompetent patients complicated with ARDS during HAdV infection and identified biomarkers for disease severity/progression. Forty-eight cytokines/chemokines in the plasma samples from 19 HAdV-infected immunocompetent adolescent and adult patients (ten complicated with ARDS) were measured and analyzed in combination with clinical indices. Immunocompetent patients with ARDS caused by severe acute respiratory disease coronavirus (SARS-CoV)-2, 2009 pandemic H1N1 (panH1N1) or bacteria were included for comparative analyses. Similar indices of disease course/progression were found in immunocompetent patients with ARDS caused by HAdV, SARS-CoV-2 or panH1N infections, whereas the HAdV-infected group showed a higher prevalence of viremia, as well as increased levels of aspartate aminotransferase (AST), alanine aminotransferase (ALT) and creatine kinase (CK). Expression levels of 33 cytokines/chemokines were increased significantly in HAdV-infected patients with ARDS compared with that in healthy controls, and many of them were also significantly higher than those in SARS-CoV-2-infected and panH1N1-infected patients. Expression of interferon (IFN)-γ, interleukin (IL)-1β, hepatocyte growth factor (HGF), monokine induced by IFN-γ (MIG), IL-6, macrophage-colony stimulating factor (M-CSF), IL-10, IL-1α and IL-2Ra was significantly higher in HAdV-infected patients with ARDS than that in those without ARDS, and negatively associated with the ratio of the partial pressure of oxygen in arterial blood/fraction of inspired oxygen (PaO2/FiO2). Analyses of the receiver operating characteristic curve (ROC) showed that expression of IL-10, M-CSF, MIG, HGF, IL-1β, IFN-γ and IL-2Ra could predict the progression of HAdV infection, with the highest area under the curve (AUC) of 0.944 obtained for IL-10. Of note, the AUC value for the combination of IL-10, IFN-γ, and M-CSF reached 1. In conclusion, the “cytokine storm” occurred during HAdV infection in immunocompetent patients, and expression of IL-10, M-CSF, MIG, HGF, IL-1β, IFN-γ and IL-2Ra was closely associated with disease severity and could predict disease progression.

Introduction

Human adenoviruses (HAdVs) are double-stranded, non-enveloped icosahedral viruses that typically cause mild infections involving the upper or lower respiratory tract, gastrointestinal tract, or conjunctiva (13). Until recently, a total of 104 subtypes within A–G species were identified, categorized either by virus neutralization assay or genomic/bioinformatics analyses (http://hadvwg.gmu.edu/). HAdV infections occurred mainly in children and immunocompromised patients, owing to a lack of humoral immunity or impaired immunity (1). In closed or crowded settings, HAdVs are highly contagious and epidemics may occur even in healthy children and adults (e.g., military recruits) (47). HAdV infections in the respiratory tract in immunocompetent adults can also occur, typically with mild and self-limiting symptoms in most cases (1, 2). As a result, HAdVs are rarely detected by clinicians and public-health surveillance in such settings. However, increasingly severe HAdV infections complicated by acute respiratory distress syndrome (ARDS) in immunocompetent adolescents and adults have been reported with a mortality prevalence of ≤80% (816), which poses a serious threat to public health. The ongoing global pandemic due to severe acute respiratory coronavirus-2 (SARS-CoV-2) infection and seasonally prevalent 2009 pandemic H1N1 (panH1N1) infection can cause severe infections in immunocompetent adults (1720). Due to distinct prognoses and treatments among these respiratory viruses, differentiating HAdV infections accurately and in a timely manner is crucial for clinicians and epidemiologists.

The “cytokine storm” (CS) is the release of circulating cytokines/chemokines during an acute infection. The CS has been shown to be correlated directly with acute lung injury and the development of ARDS during viral infections (21, 22), including influenza viruses and coronaviruses (2330). Moreover, some cytokines/chemokines have been found to predict the disease progression accurately (24, 26, 3133). HAdV has also been shown to be a proinflammatory virus that can trigger the release of high levels of inflammatory cytokines and chemokines in children patients, and the expression levels differed among different disease severity (3436). However, the specific expression profiles of cytokines/chemokines in immunocompetent adolescent and adult patients, especially those complicated with ARDS caused by HAdV infection, have rarely been reported.

In this study, we characterized the cytokine/chemokine expression profiles of immunocompetent patients complicated with ARDS during HAdV infection, and compared them against other common causative agents of ARDS, including SARS-CoV-2, panH1N1 and bacteria. Furthermore, we identified biomarkers for disease severity and progression in HAdV infections.

Materials and Methods

Ethical Approval of the Study Protocol

The study protocol was approved (SZTHEC2016001) by the Ethics Committees of Shenzhen Third People’s Hospital (Shenzen, China). Written informed consent was obtained from patients infected with HAdV, panH1N1, or suffering from bacterial pneumonia. Verbal informed consent was obtained from coronavirus 2019 (COVID-19) patients because pen and paper were not allowed in containment facilities. The study was conducted in accordance with the International Conference on Harmonization Guidelines for Good Clinical Practice, the Declaration of Helsinki 1965 and later amendments, and the ethical guidelines of Shenzhen Third People’s Hospital.

Patient Information and Data Collection

This was a retrospective study. Study participants were immunocompetent patients hospitalized in Shenzhen Third People’s Hospital from 2018 to 2020. Nineteen laboratory analyses-confirmed HAdV-infected patients (10 complicated with ARDS and nine without ARDS) were enrolled. Patients complicated with ARDS caused by infection with SARS-CoV-2 (n = 28), panH1N1 (n = 21), suffering from bacterial pneumonia (n = 10) and healthy controls (HCs) (n = 8) were also included for analyses. All patients complicated with ARDS met the criteria of the Berlin definition for the diagnosis of ARDS (37). Plasma and respiratory specimens were collected from patients and HCs. Clinical information, including complete blood counts and blood biochemistry, were collected at the earliest time point after hospitalization. Immunocompromised patients infected with the human immunodeficiency virus, suffering from neutropenia, receiving immunosuppressive chemotherapy, with malignant neoplasms, or pregnant or breastfeeding women, were excluded, as reported previously (38).

Real-Time Reverse Transcription-Quantitative Polymerase Chain Reaction (RT-qPCR) and Next-Generation Sequencing (NGS)

Sputum or nasopharyngeal swabs were collected from patients at various time-points after hospitalization. RT-qPCR was undertaken as described previously (39), and using commercial qRT-PCR kits for the detection of HAdV (Mabsky Biotech Co., Ltd., Shenzhen, China). Samples positive for HAdVs were subjected to virus isolation using HEP-2 cells and NGS (Guangzhou Sagene Biotech, Guangzhou, China).

Measurement of Expression of Cytokines and Chemokines

The plasma samples of laboratory analyses-confirmed patients were collected at the earliest possible time-point after hospitalization and thereafter. The plasma of HCs was included as the negative control. The concentrations of 48 cytokines including adaptive immunity cytokines, pro-inflammatory cytokines, and anti-inflammatory cytokines associated with infectious diseases and inflammation were measured using the Bio-Plex Pro Human Cytokine Screening Panel (Bio-Rad Laboratories, Hercules, CA, USA) on a Luminex™ 200 instrument (Merck Millipore, Burlington, MA, USA) following manufacturer instructions, as reported previously (31).

Quantification of Hypoxia and Lung Injury

The partial pressure of oxygen in arterial blood (PaO2) was measured by a blood-gas analyzer (ABL90; Radiometer, Copenhagen, Denmark) at various time-points after hospitalization, as reported previously (27, 31). The fraction of inspired oxygen (FiO2) was calculated using the following formula: FiO2 = (21 + oxygen flow [in units of L/minute] × 4)/100. The PaO2/FiO2 ratio (in mmHg) was obtained by dividing the PaO2 value with the FiO2 value. A PaO2/FiO2 ratio ≤100 mmHg is considered one of the criteria for severe ARDS, as defined previously (37).

Statistical Analysis

The Mann–Whitney U-test was employed to determine the differences between groups of continuous variables. The Fisher exact test was used for categorical variables. The Spearman correlation coefficient was employed to analyze the linear correlation. The area under the receiver operating characteristic (ROC) curve (AUC) of plasma cytokine levels was estimated for patients developing ARDS or not developing ARDS. Moreover, the combined values for the prediction of developing ARDS was calculated using binary logistic regression. Statistical values were calculated using SPSS 20.0 (IBM, Armonk, NY, USA). A P-value of 0.01–0.05, 0.001–0.01 and <0.001 was considered significant, very significant, and extremely significant, respectively.

Results

Distinct Clinical Characteristics of Immunocompetent Patients Complicated With ARDS Infected by HAdV, SARS-CoV-2, panH1N1 or Bacteria

Nineteen patients with laboratory analyses-confirmed HAdVs with ten patients complicated by ARDS (case No. 01-10) were included in this study, The patients were predominantly male (78.9%) and <60 years of age (94.7%) (Table S1). Patients initially showed influenza-like symptoms that developed rapidly into pneumonia and ARDS within 6 days after illness onset (d.a.o) in most cases of the ARDS group (Table S1). Moreover, only HAdV-7 (8/10) and HAdV-55 (2/10) were found in the ARDS group, and seven patients received antiviral treatment with cidofovir (Table S1). To compare the distinct characteristics of HAdV infections with those of other important respiratory infections, clinical information was first collected and analyzed from immunocompetent patients with ARDS caused by SARS-CoV-2 (n = 28, collected in 2020), panH1N1 (n = 21, collected during 2018–2019) and bacteria (n = 10, collected in 2019) (Table S2). The median age of the HAdV-infected group was significantly younger than that of the other groups, and also a lower instance of underlying diseases. The median number of days from disease onset to hospital admission and hospitalization duration were similar among the four groups. There was no significant difference in the proportion of patients developing severe ARDS, needing mechanical ventilation or transfer to the intensive unit care unit (ICU), or in-hospital prevalence of death (Table S2). However, compared with patients infected with SARS-CoV-2 or panH1N1, a higher prevalence of hepatic insufficiency and cardiac failure was found in HAdV-infected patients, though the difference was not significant in some cases. Notably, a significantly higher prevalence of viremia in HAdV-infected patients (90%) was found (Table S2). The prevalence of increased levels of aspartate aminotransferase (AST), alanine aminotransferase (ALT) and creatine kinase (CK) in HAdV-infected patients upon hospital admission was significantly higher than that in the other groups, especially compared with SARS-CoV-2-infected and panH1N1-infected patients (Table 1). Moreover, the prevalence of lymphopenia in HAdV-infected patients was significantly higher than that in SARS-CoV-2-infected cases and bacterial-pneumonia patients (Table 1). Also, the prevalence of leukopenia and thrombocytopenia in HAdV-infected patients was significantly higher than that in panH1N1-infected patients (Table 1).

TABLE 1
www.frontiersin.org

Table 1 Laboratory results of hospitalized patients complicated with ARDS caused by HAdV, SARS-CoV-2, panH1N1 and bacteria.

Comparison of Cytokines/Chemokines Expression Profiles Amongst Immunocompetent Patients Complicated With ARDS Infected by HAdV, SARS-CoV-2, panH1N1 or Bacteria

Expression levels of cytokines/chemokines indicate disease severity in some other respiratory infections (24, 26, 27, 31). Hence, the concentrations of 48 cytokines/chemokines including adaptive immunity cytokines, pro-inflammatory cytokines, and anti-inflammatory cytokines associated with infectious diseases and inflammation were measured and compared in patients infected with HAdV, SARS-CoV-2, panH1N1, or suffering from bacterial pneumonia complicated with ARDS (Figures 1 and S1). Expression of 33 cytokines/chemokines was significantly higher in HAdV-infected patients complicated with ARDS than that in HCs (Figures 1 and S1), and expression levels of most of these cytokines were significantly positively correlated with the HAdV load (Figure S2). Among these 33 cytokines/chemokines, significantly increased expression of interleukin (IL)-1α, IL-1β, IL-2Ra, IL-10, IL-18, granulocyte-colony stimulating factor (G-CSF), interferon (IFN)-α2, IFN-γ, interferon gamma-induced protein (IP)-10, macrophage inflammatory protein (MIP)-1α, cutaneous T cell-attracting chemokine (CTACK), macrophage-colony stimulating factor (M-CSF), and monokine induced by IFN-γ (MIG) were also found in the remaining three groups, whereas concentrations in the HAdV-infected group were significantly higher than those in the SARS-CoV-2-infected group and panH1N1-infected group (Figures 1 and S1). Interestingly, significantly increased expression of monocyte chemoattractant protein (MCP)-1, tumor necrosis factor (TNF)-α and stromal cell-derived factor (SDF)-1α was found only in HAdV-infected patients among the four groups, and significantly increased levels of IL-17, Skp, Cullin, F-box containing complex (SCF)-β and tumor necrosis factor-related apoptosis inducing ligand (TRAIL) were found specifically in the HAdV-infected group among the three viral infection groups (Figures 1 and S1). These results indicated that the CS occurred in HAdV infection, and that it was possibly even more severe than that in SARS-CoV-2 and panH1N1 infections in immunocompetent patients.

FIGURE 1
www.frontiersin.org

Figure 1 Comparison of cytokines/chemokines expression profiles among patients infected with HAdV, SARS-CoV-2, panH1N1 or bacteria complicated with ARDS upon admission. Samples from patients infected with HAdV (n = 9), SARS-CoV-2 (n = 28) and panH1N1 (n = 21) were collected at the earliest possible time-point after hospitalization and subjected to cytokines/chemokines measurement. Healthy controls (n = 8) and patients with bacterial pneumonia (n = 10) were included as controls. Results of 33 cytokines/chemokines showing significantly increased expression in the HAdV group are shown. Statistical analyses with the reference group of HAdV and HCs are shown in black and red, respectively. P = 0.01–0.05, 0.001–0.01 and <0.001 was considered significant (*), very significant (**) and extremely significant (***), respectively, whereas ns represents not significant.

Differential Expression of Cytokines/Chemokines in HAdV Patients Complicated With or Not Complicated With ARDS

We wished to investigate the cytokines/chemokines that might associated with the disease severity of HAdV infection. Hence, the expression profile of these 33 cytokines/chemokines in patients with or without ARDS at different phases of disease were analyzed further (Figure 2). The date of sample collection was stratified into three groups according to disease progression, as reported previously (32), the first 7 days after illness onset (0~7 d.a.o), between 8 and 14 days following illness onset (8~14 d.a.o), and during the recovery phase from 15 days after disease onset (≥ 15 d.a.o). During the first 7 d.a.o, expression of all 33 cytokines/chemokines was increased in HAdV-infected patients with ARDS, whereas expression of only IFN-γ, IL-12, IL-18, fibroblast growth factor (FGF) basic, G-CSF, IP-10 and CTACK was increased in patients not suffering from ARDS. During 8~14 d.a.o, expression of TNF-α, SDF-1α, TRAIL, IL-5, IL-7, IL-13 and IL-17 in patients with ARDS returned to normal, and expression of 10 cytokines/chemokines (IFN-γ, IL-1α, IL-1β, IL-18, IL-4, MIP-1α, G-CSF, IP-10, CTACK and SCF-β) was increased in patients without ARDS. Furthermore, expression of nine cytokines/chemokines (hepatocyte growth factor (HGF), MIG, IFN-γ, IL-1β, IL-6, M-CSF, IL-10, IL-1α and IL-2Ra) was significantly higher in HAdV-infected patients with ARDS than in those without ARDS during disease progression. Moreover, higher expression of MIG, M-CSF, IL-6, IL-1α and IL-2Ra was mainly found within 14 d.a.o, and higher expression of IFN-γ, IL-1β and IL-10 was found during the entire disease progression. Notably, most of these cytokines/chemokines showing significantly increased expression in patients with ARDS showed an obvious trend of decrease after 14 d.a.o. These results suggested a possible role of cytokines/chemokines in the pathogenesis of HAdV infection, especially the differentially expressed cytokines/chemokines between patients with or without ARDS.

FIGURE 2
www.frontiersin.org

Figure 2 Comparison of cytokines/chemokines showing significant increases in HAdV patients with or without ARDS during disease progression. Expression of 33 cytokines/chemokines measured at different days after illness onset were compared between HAdV patients with or without ARDS. P = 0.01–0.05, 0.001–0.01 and 0.0001–0.001 was considered significant (*), very significant (**) and extremely significant (***), respectively, whereas ns represents not significant.

Expression of IL-10, M-CSF, MIG, HGF, IL-1β, IFN-γ and IL-2Ra Was Closely Associated With Disease Severity and Could Predict Disease Progression in HAdV Infection

It has been shown that concentrations of specific cytokines/chemokines are associated with disease severity in H7N9, H5N1, H5N6 and SARS-CoV-2 infections (24, 26, 27, 31, 40, 41). Spearman coefficient correlation analysis was done to analyze the association between expression of the nine differentially expressed cytokines/chemokines and the PaO2/FiO2 ratio in HAdV-infected patients (Figure 3A). The concentration of the nine cytokines/chemokines was negatively correlated with the PaO2/FiO2 ratio, with the r values ranging from −0.5519 (IL-2Ra) to −0.6915 (MIG) and P < 0.05, which indicated that these cytokines/chemokines were associated with disease severity in HAdV infection. The dynamic change in expression of these nine cytokines/chemokines in HAdV-infected patients complicated with ARDS was also measured and analyzed (Figure S3). The concentrations of these cytokines/chemokines in most cases showed an obvious decreasing trend during disease progression. In addition, we calculated the AUC to test whether the nine cytokines/chemokines could be biomarkers for the prediction of developing ARDS in HAdV-infected patients. These 17 HAdV-infected patients were divided into an ARDS group and non-ARDS group. Expression of the nine cytokines/chemokines at the earliest time-point after hospital admission was used for the calculation (Figures 3B and S4A). The AUC for IL-10, M-CSF, MIG, HGF, IL-1β, IFN-γ and IL-2Ra was >0.8 and P < 0.05, and IL-10 possessed the highest AUC (0.944). Then we further analyzed the AUC values of different combination of IL-10 with the other ones (Figure 3C, Figure S4B, C), and found that combination of IL-10 and M-CSF increased the value of AUC to 0.958 (Figure 3C). More surprisingly, the AUC value for the combination of IL-10, IFN-γ, and M-CSF reached 1 (Figure 3C).

FIGURE 3
www.frontiersin.org

Figure 3 Expression of IL-10, M-CSF, MIG, HGF, IL-1β, IFN-γ and IL-2Ra was closely associated with disease severity and could predict disease progression. (A) Spearman rank correlation analysis between expression of HGF, MIG, IFN-γ, IL-1β, IL-6, M-CSF, IL-10, IL-1α and IL-2Ra measured from plasma samples and the corresponding PaO2/FiO2 ratio on the same day. (B) ROC curves of expression levels of HGF, MIG, IFN-γ, IL-1β, M-CSF, IL-10, and IL-2Ra upon hospital admission for HAdV-infected patients with and without ARDS during hospitalization. (C) ROC curves of different combination of IL-10, IFN-γ, and M-CSF expression levels upon hospital admission for HAdV-infected patients with and without ARDS during hospitalization. All P-values for the ROC curves were <0.05.

Discussion

Initially, most respiratory infections show influenza-like symptoms. ARDS (which carries high mortality and morbidity) serves as a common complication during severe respiratory infections (37). Compared with severe infection with SARS-CoV-2 or panH1N1, the overall disease course of severe HAdV infection (including the interval from disease onset to hospital admission and hospitalization duration) and disease progression (including developing severe ARDS, utilization of mechanical ventilation, requiring ICU transfer, and the in-hospital mortality prevalence) were similar (Table S2). These data indicated that the severity of ARDS caused by these three viruses was similar. However, unlike influenza viruses and SARS-CoV-2, laboratory testing of HAdVs is not a routine examination for patients with pneumonia in most hospitals. Thus, it is easy to ignore HAdV-induced severe infections, especially for immunocompetent adolescents and adults (42). Apart from the laboratory-based diagnosis, distinct clinical characteristics could also provide some hints for clinicians. For example, severe infection with SARS-CoV-2 or panH1N1 was found mainly in elderly patients with underlying diseases, whereas most HAdV-infected patients with ARDS reported in this study and previous studies were under 60 years of age (12, 43, 44), much younger than those infected with SARS-CoV-2 or panH1N1 (Table S2). Meanwhile, with regard to complete blood counts and blood biochemistry, severe HAdV infections were characterized by significantly higher levels of AST, ALT, and CK than those in SARS-CoV-2-infected and panH1N1-infected patients (Table 1). Therefore, immunocompetent adolescents and young adults with respiratory symptoms and significantly elevated AST, ALT, and CK were proposed to be at higher risk of HAdV infection when compared with other respiratory viruses. Currently, options for efficacious antiviral treatment of HAdV infections are limited. Studies have shown that cidofovir and its lipid ester analog, brincidofovir, may be the most promising drugs (2, 45). Fortunately, HAdV infections were confirmed by laboratory analyses within 5 days after hospital admission for all patients with ARDS, and six patients received cidofovir treatment and survived in our study cohort (Table S1). These results highlight the importance of a timely laboratory-based diagnosis of HAdV infection and the benefit of cidofovir treatment for severe HAdV infection. Similar to data from other studies (2, 46), HAdV-7 and HAdV-55 in species B were the main serotypes that caused severe pneumonia in immunocompetent adolescent and adults, and the underlying mechanism merits further investigation.

The CS with uncontrolled proinflammatory responses has important roles in the immunopathogenesis of infection and is associated with the severity and outcome of disease (21, 22, 47, 48). Increased expression of cytokines/chemokines, and several clinical abnormalities reflecting hyperinflammation and tissue damage in CS disorders (22), (e.g., increased levels of C-reactive protein and lactate dehydrogenase, hypoalbuminemia, thrombocytopenia, hepatic and renal insufficiency) were also observed in HAdV-infected patients (Table 1), indicating that the CS occurred following HAdV infection. Totally, 33 tested cytokines/chemokines were found to be significantly elevated in HAdV patients with ARDS, among which elevation of IP-10, TNF-α, IL-1α, IL-1β, IL-6, MIP-1α, MIG, and IFN-γ have also been found from children patients in previous studies (34, 49). When comparing with SARS-CoV-2 or panH1N1 infections, some cytokines were commonly elevated, while elevated cytokines including MCP-1, TNF-α, SDF-1α, IL-17, SCF-β, and TRAIL were found specifically in the HAdV-infected group (Figure 1 and Figure S1). Moreover, many of the cytokines/chemokines showing increased expression in HAdV-infected patients complicated with ARDS also showed significantly higher expression than those with SARS-CoV-2 or panH1N1 infections (Figures 1 and S1). These results suggested the common and specific roles of some cytokines in different respiratory infections and possibly a more severe CS in HAdV-infected patients complicated with ARDS than SARS-CoV-2 or panH1N1, which merits further investigation. Furthermore, the concentrations of nine cytokines including IL-6, M-CSF, IL-1β, MIG which are monocyte/macrophage activation associated biomarkers (41), IFN-γ which is the sole type II IFN and important in early host defense against infection (50), IL-1α which activates the inflammatory process (51), IL-10 which regulates and suppresses the expression of proinflammatory cytokines (52), HGF, and IL-2Ra were significantly higher in HAdV-infected patients complicated with ARDS, and negatively correlated with the PaO2/FiO2 ratio (Figures 2, 3). Hence, these cytokines might have crucial roles in the pathophysiology of HAdV infection. Moreover, studies have aimed to establish prediction models for disease progression based on expression of cytokines/chemokines in plasma (24, 26, 3133). Specific cytokines/chemokines could accurately predict disease progression, such as IP-10 for COVID-19, and Ang II for H7N9 infection (24, 26, 3133). Unlike the other viruses, IL-10 in HAdV infection serves as the best predictor for disease progression with the highest area under the curve (AUC) of 0.944 (Figure 3B), indicating that predictors for different viral infections were virus specific. Notably, we also found the AUC value for the combination of IL-10, IFN-γ, and M-CSF reached 1 (Figure 3C), therefore, expression levels the three cytokines might fully predict the disease severity of HAdV infection. Given the crucial role of the CS in severe inflammation and damage to vital organs, short courses of treatment with corticosteroids at low-to-moderate doses might be beneficial against severe infections by HAdVs, panH1N1, H7N9 and SARS-CoV-2 (53, 54).

In summary, the CS occurred during HAdV infection. Expression of IL-10, M-CSF, MIG, HGF, IL-1β, IFN-γ and IL-2Ra was closely associated with disease severity and could predict disease progression. These results aid understanding of the clinical characteristics and immunopathologic mechanisms of severe HAdV infection in immunocompetent adolescents and adults. Our data also highlight the increasing public-health threat of emerging (or re-emerging) HAdV subtypes, especially HAdV-7 and HAdV-55.

Age Specific Reference Ranges Used to Define Abnormalities in Blood Results

Leukopenia (× 109/L): 2 months-2 years: <5, >2 years: <4; Lymphopenia (× 109/L): 2-11 months: <4.0, 1-11 years: <1.5, 12+ years: <1; Neutropenia (× 109/L): All ages: <1.5; Neutrophilia (× 109/L): All ages: >8.5; Thrombocytopenia (× 109/L): All ages: <150; Elevated AST(U/L): All ages: >50; Elevated ALT (U/L): All ages: >50; Elevated CRE (µmol/L): All ages: >120; Elevated CK (U/L): All ages: >200; Elevated CRP (nmol/L): All ages: >10; Elevated LDH (U/L): <7 years: >400, 7-15 years: >300, 16+ years: >250.

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.

Ethics Statement

The studies involving human participants were reviewed and approved by Ethics Committees of Shenzhen Third People’s Hospital. Written informed consent to participate in this study was provided by the participants’ legal guardian/next of kin.

Author Contributions

JL, YY, and JC contributed to the analysis and interpretation of data. JW, ZX, ML, MY, CJ, YJL, YG, and FW enrolled the patients and collected the data. YS, LS, and YXL were responsible for the concept and design of the study; JL drafted the article. YXL, LS, YS, and YY were responsible for the critical revision for important intellectual content and for the final approval of the article. All authors agree to be accountable for the content of the work. All authors contributed to the article and approved the submitted version.

Funding

This work was supported by grants from the Shenzhen Science and Technology Research and Development Project [grant no. 202002073000001 and JSGG20200225152008136], National Natural Science Foundation of China [grant no. 31971147] and Shenzhen Science and Technology Innovation Commission [grant no. JCYJ20190807155011406 and KQTD20180411143323605] and National Science and Technology Major Project [grant no. 2018ZX10713001], the Ministry of Science and Technology (2020YFC0846300).

Conflict of Interest

JC was employed by the company Guangzhou Sagene Biotech.

The remaining 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 thank the participants of this study and colleges in Shenzhen Third People’s hospital for the data and sample collection.

Supplementary Material

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

Supplementary Figure 1 | HeatMap of differentially expressed cytokines/chemokines. The HeatMap of the differentially expressed cytokines/chemokines between HAdV-infected patients against (A) healthy controls, (B) SARS-CoV-2-infected patients, (C) panH1N1-infected patients (C, D) bacterial-pneumonia patients.

Supplementary Figure 2 | Correlation between cytokines/chemokines expression and HAdV viral load. Expression of the 26 cytokines/chemokines showing significant increases measured from plasma samples collected upon hospital admission and thereafter and the corresponding viral load (indicated as log (copies/mL)) on the same day were evaluated using Spearman rank correlation analysis.

Supplementary Figure 3 | Dynamic changes in expression of nine differentially expressed cytokines/chemokines measured in HAdV patients with or without ARDS at the indicated time-points from HAdV patients complicated with ARDS. The fatal case (case 06) is marked in red.

Supplementary Figure 4 | (A) ROC curves of expression levels of IL-6 and IL-1α. (B) ROC curves of the combination of IL-10 with HGF, MIG, IL-1β, M-CSF, IL-10, and IL-2Ra expression levels upon hospital admission for HAdV-infected patients with and without ARDS during hospitalization. (C) ROC curves of the combination of IL-10, IFN-γ with HGF, MIG, IL-1β, and IL-2Ra expression levels upon hospital admission for HAdV-infected patients with and without ARDS during hospitalization. All P-values for the ROC curves were <0.05 except IL-6 and IL-1α.

Abbreviations

WBC, white blood cell; LYM, lymphocyte; NEU, neutrophil; PLT, platelet; AST, aspartate aminotransferase; ALT, alanine aminotransferase; CRE, serum creatinine; BUN, blood urea nitrogen; CK, creatine kinase; CRP, C-reactive protein; LDH, lactate dehydrogenase.

References

1. Lynch JP 3rd, Fishbein M, Echavarria M. Adenovirus. Semin Respir Crit Care Med (2011) 32(4):494–511. doi: 10.1055/s-0031-1283287

PubMed Abstract | CrossRef Full Text | Google Scholar

2. Ison MG, Hayden RT. Adenovirus. Microbiol Spectr (2016) 4(4). doi: 10.1128/microbiolspec.DMIH2-0020-2015

PubMed Abstract | CrossRef Full Text | Google Scholar

3. Green M, Pina M, Kimes R, Wensink PC, MacHattie LA, Thomas CA Jr. Adenovirus DNA. I. Molecular Weight and Conformation. Proc Natl Acad Sci U S A (1967) 57(5):1302–9. doi: 10.1073/pnas.57.5.1302

PubMed Abstract | CrossRef Full Text | Google Scholar

4. Sanchez JL, Binn LN, Innis BL, Reynolds RD, Lee T, Mitchell-Raymundo F, et al. Epidemic of Adenovirus-Induced Respiratory Illness Among US Military Recruits: Epidemiologic and Immunologic Risk Factors in Healthy, Young Adults. J Med Virol (2001) 65(4):710–8. doi: 10.1002/jmv.2095

PubMed Abstract | CrossRef Full Text | Google Scholar

5. Kajon AE, Moseley JM, Metzgar D, Huong HS, Wadleigh A, Ryan MA, et al. Molecular Epidemiology of Adenovirus Type 4 Infections in US Military Recruits in the Postvaccination Era (1997-2003). J Infect Dis (2007) 196(1):67–75. doi: 10.1086/518442

PubMed Abstract | CrossRef Full Text | Google Scholar

6. Ryan MA, Gray GC, Smith B, McKeehan JA, Hawksworth AW, Malasig MD. Large Epidemic of Respiratory Illness Due to Adenovirus Types 7 and 3 in Healthy Young Adults. Clin Infect Dis (2002) 34(5):577–82. doi: 10.1086/338471

PubMed Abstract | CrossRef Full Text | Google Scholar

7. Centers for Disease Control and Prevention (CDC). Civilian Outbreak of Adenovirus Acute Respiratory Disease–South Dakota, 1997. MMWR Morb Mortal Wkly Rep (1998) 47(27):567–70.

PubMed Abstract | Google Scholar

8. Yu Z, Zeng Z, Zhang J, Pan Y, Chen M, Guo Y, et al. Fatal Community-Acquired Pneumonia in Children Caused by Re-Emergent Human Adenovirus 7d Associated With Higher Severity of Illness and Fatality Rate. Sci Rep (2016) 6:37216. doi: 10.1038/srep37216

PubMed Abstract | CrossRef Full Text | Google Scholar

9. Kajon AE, Ison MG. Severe Infections With Human Adenovirus 7d in 2 Adults in Family, Illinois, USA, 2014. Emerg Infect Dis (2016) 22(4):730–3. doi: 10.3201/eid2204.151403

PubMed Abstract | CrossRef Full Text | Google Scholar

10. Scott MK, Chommanard C, Lu X, Appelgate D, Grenz L, Schneider E, et al. Human Adenovirus Associated With Severe Respiratory Infection, Oregon, USA, 2013-2014. Emerg Infect Dis (2016) 22(6):1044–51. doi: 10.3201/eid2206.151898

PubMed Abstract | CrossRef Full Text | Google Scholar

11. Lafolie J, Mirand A, Salmona M, Lautrette A, Archimbaud C, Brebion A, et al. Severe Pneumonia Associated With Adenovirus Type 55 Infection, France, 2014. Emerg Infect Dis (2016) 22(11):2012–4. doi: 10.3201/eid2211.160728

PubMed Abstract | CrossRef Full Text | Google Scholar

12. Zhang SY, Luo YP, Huang DD, Fan H, Lu QB, Wo Y, et al. Fatal Pneumonia Cases Caused by Human Adenovirus 55 in Immunocompetent Adults. Infect Dis (Lond) (2016) 48(1):40–7. doi: 10.3109/23744235.2015.1055585

PubMed Abstract | CrossRef Full Text | Google Scholar

13. Yamamoto D, Okamoto M, Lupisan S, Suzuki A, Saito M, Tamaki R, et al. Impact of Human Adenovirus Serotype 7 in Hospitalized Children With Severe Fatal Pneumonia in the Philippines. Jpn J Infect Dis (2014) 67(2):105–10. doi: 10.7883/yoken.67.105

PubMed Abstract | CrossRef Full Text | Google Scholar

14. Pfortmueller CA, Barbani MT, Schefold JC, Hage E, Heim A, Zimmerli S. Severe Acute Respiratory Distress Syndrome (ARDS) Induced by Human Adenovirus B21: Report on 2 Cases and Literature Review. J Crit Care (2019) 51:99–104. doi: 10.1016/j.jcrc.2019.02.019

PubMed Abstract | CrossRef Full Text | Google Scholar

15. Gu L, Liu Z, Li X, Qu J, Guan W, Liu Y, et al. Severe Community-Acquired Pneumonia Caused by Adenovirus Type 11 in Immunocompetent Adults in Beijing. J Clin Virol (2012) 54(4):295–301. doi: 10.1016/j.jcv.2012.04.018

PubMed Abstract | CrossRef Full Text | Google Scholar

16. Tan D, Zhu H, Fu Y, Tong F, Yao D, Walline J, et al. Severe Community-Acquired Pneumonia Caused by Human Adenovirus in Immunocompetent Adults: A Multicenter Case Series. PloS One (2016) 11(3):e0151199. doi: 10.1371/journal.pone.0151199

PubMed Abstract | CrossRef Full Text | Google Scholar

17. Skarbinski J, Jain S, Bramley A, Lee EJ, Huang J, Kirschke D, et al. Hospitalized Patients With 2009 Pandemic Influenza a (H1N1) Virus Infection in the United States–September-October 2009. Clin Infect Dis (2011) 52 Suppl 1:S50–9. doi: 10.1093/cid/ciq021

PubMed Abstract | CrossRef Full Text | Google Scholar

18. Jain S, Kamimoto L, Bramley AM, Schmitz AM, Benoit SR, Louie J, et al. Hospitalized Patients With 2009 H1N1 Influenza in the United States, April-June 2009. N. Engl J Med (2009) 361(20):1935–44. doi: 10.1056/NEJMoa0906695

CrossRef Full Text | Google Scholar

19. Wu J, Li J, Zhu G, Zhang Y, Bi Z, Yu Y, et al. Clinical Features of Maintenance Hemodialysis Patients With 2019 Novel Coronavirus-Infected Pneumonia in Wuhan, China. Clin J Am Soc Nephrol (2020) 15(8):1139–45. doi: 10.2215/CJN.04160320

PubMed Abstract | CrossRef Full Text | Google Scholar

20. Wang C, Horby PW, Hayden FG, Gao GF. A Novel Coronavirus Outbreak of Global Health Concern. Lancet (2020) 395(10223):470–3. doi: 10.1016/S0140-6736(20)30185-9

PubMed Abstract | CrossRef Full Text | Google Scholar

21. D’Elia RV, Harrison K, Oyston PC, Lukaszewski RA, Clark GC. Targeting the “Cytokine Storm” for Therapeutic Benefit. Clin Vaccine Immunol (2013) 20(3):319–27. doi: 10.1128/CVI.00636-12

PubMed Abstract | CrossRef Full Text | Google Scholar

22. Fajgenbaum DC, June CH. Cytokine Storm. N Engl J Med (2020) 383(23):2255–73. doi: 10.1056/NEJMra2026131

PubMed Abstract | CrossRef Full Text | Google Scholar

23. Zeng H, Belser JA, Goldsmith CS, Gustin KM, Veguilla V, Katz JM, et al. A(H7N9) Virus Results in Early Induction of Proinflammatory Cytokine Responses in Both Human Lung Epithelial and Endothelial Cells and Shows Increased Human Adaptation Compared With Avian H5N1 Virus. J Virol (2015) 89(8):4655–67. doi: 10.1128/JVI.03095-14

PubMed Abstract | CrossRef Full Text | Google Scholar

24. Guo J, Huang F, Liu J, Chen Y, Wang W, Cao B, et al. The Serum Profile of Hypercytokinemia Factors Identified in H7N9-Infected Patients can Predict Fatal Outcomes. Sci Rep (2015) 5:10942. doi: 10.1038/srep10942

PubMed Abstract | CrossRef Full Text | Google Scholar

25. Yang Y, Wong G, Yang L, Tan S, Li J, Bai B, et al. Comparison Between Human Infections Caused by Highly and Low Pathogenic H7N9 Avian Influenza Viruses in Wave Five: Clinical and Virological Findings. J Infect (2019) 78(3):241–8. doi: 10.1016/j.jinf.2019.01.005

PubMed Abstract | CrossRef Full Text | Google Scholar

26. de Jong MD, Simmons CP, Thanh TT, Hien VM, Smith GJ, Chau TN, et al. Fatal Outcome of Human Influenza a (H5N1) is Associated With High Viral Load and Hypercytokinemia. Nat Med (2006) 12(10):1203–7. doi: 10.1038/nm1477

PubMed Abstract | CrossRef Full Text | Google Scholar

27. Bi Y, Tan S, Yang Y, Wong G, Zhao M, Zhang Q, et al. Clinical and Immunological Characteristics of Human Infections With H5N6 Avian Influenza Virus. Clin Infect Dis (2019) 68(7):1100–9. doi: 10.1093/cid/ciy681

PubMed Abstract | CrossRef Full Text | Google Scholar

28. Jiang Y, Xu J, Zhou C, Wu Z, Zhong S, Liu J, et al. Characterization of Cytokine/Chemokine Profiles of Severe Acute Respiratory Syndrome. Am J Respir Crit Care Med (2005) 171(8):850–7. doi: 10.1164/rccm.200407-857OC

PubMed Abstract | CrossRef Full Text | Google Scholar

29. Mubarak A, Alturaiki W, Hemida MG. Middle East Respiratory Syndrome Coronavirus (MERS-Cov): Infection, Immunological Response, and Vaccine Development. J Immunol Res (2019) 2019:6491738. doi: 10.1155/2019/6491738

PubMed Abstract | CrossRef Full Text | Google Scholar

30. Mahallawi WH, Khabour OF, Zhang Q, Makhdoum HM, Suliman BA. MERS-Cov Infection in Humans is Associated With a Pro-Inflammatory Th1 and Th17 Cytokine Profile. Cytokine (2018) 104:8–13. doi: 10.1016/j.cyto.2018.01.025

PubMed Abstract | CrossRef Full Text | Google Scholar

31. Yang Y, Shen C, Li J, Yuan J, Wei J, Huang F, et al. Plasma IP-10 and MCP-3 Levels are Highly Associated With Disease Severity and Predict the Progression of COVID-19. J Allergy Clin Immunol (2020) 146(1):119–27 e4. doi: 10.1016/j.jaci.2020.04.027

PubMed Abstract | CrossRef Full Text | Google Scholar

32. Huang F, Guo J, Zou Z, Liu J, Cao B, Zhang S, et al. Angiotensin II Plasma Levels are Linked to Disease Severity and Predict Fatal Outcomes in H7N9-Infected Patients. Nat Commun (2014) 5:3595. doi: 10.1038/ncomms4595

PubMed Abstract | CrossRef Full Text | Google Scholar

33. Tang NL, Chan PK, Wong CK, To KF, Wu AK, Sung YM, et al. Early Enhanced Expression of Interferon-Inducible Protein-10 (CXCL-10) and Other Chemokines Predicts Adverse Outcome in Severe Acute Respiratory Syndrome. Clin Chem (2005) 51(12):2333–40. doi: 10.1373/clinchem.2005.054460

PubMed Abstract | CrossRef Full Text | Google Scholar

34. Moro MR, Bonville CA, Suryadevara M, Cummings E, Faddoul D, Kobayaa H, et al. Clinical Features, Adenovirus Types, and Local Production of Inflammatory Mediators in Adenovirus Infections. Pediatr Infect Dis J (2009) 28(5):376–80. doi: 10.1097/INF.0b013e31819075a6

PubMed Abstract | CrossRef Full Text | Google Scholar

35. Chen RF, Lee CY. Adenoviruses Types, Cell Receptors and Local Innate Cytokines in Adenovirus Infection. Int Rev Immunol (2014) 33(1):45–53. doi: 10.3109/08830185.2013.823420

PubMed Abstract | CrossRef Full Text | Google Scholar

36. Qi L, Wang Y, Wang H, Deng J. Adenovirus 7 Induces Interlukin-6 Expression in Human Airway Epithelial Cells Via P38/NF-Kappab Signaling Pathway. Front Immunol (2020) 11:551413. doi: 10.3389/fimmu.2020.551413

PubMed Abstract | CrossRef Full Text | Google Scholar

37. Thompson BT, Chambers RC, Liu KD. Acute Respiratory Distress Syndrome. N Engl J Med (2017) 377(6):562–72. doi: 10.1056/NEJMra1608077

PubMed Abstract | CrossRef Full Text | Google Scholar

38. Gu L, Qu J, Sun B, Yu X, Li H, Cao B. Sustained Viremia and High Viral Load in Respiratory Tract Secretions are Predictors for Death in Immunocompetent Adults With Adenovirus Pneumonia. PloS One (2016) 11(8):e0160777. doi: 10.1371/journal.pone.0160777

PubMed Abstract | CrossRef Full Text | Google Scholar

39. Bil-Lula I, De Franceschi N, Pawlik K, Wozniak M. Improved Real-Time PCR Assay for Detection and Quantification of All 54 Known Types of Human Adenoviruses in Clinical Samples. Med Sci Monit (2012) 18(6):BR221–8. doi: 10.12659/msm.882898

PubMed Abstract | CrossRef Full Text | Google Scholar

40. Laing AG, Lorenc A, Del Molino Del Barrio I, Das A, Fish M, Monin L, et al. A Dynamic COVID-19 Immune Signature Includes Associations With Poor Prognosis. Nat Med (2020) 26(10):1623–35. doi: 10.1038/s41591-020-1038-6

PubMed Abstract | CrossRef Full Text | Google Scholar

41. Abers MS, Delmonte OM, Ricotta EE, Fintzi J, Fink DL, de Jesus AAA, et al. An Immune-Based Biomarker Signature is Associated With Mortality in COVID-19 Patients. JCI Insight (2021) 6(1):e144455. doi: 10.1172/jci.insight.144455

CrossRef Full Text | Google Scholar

42. Cederwall S, Pahlman LI. Respiratory Adenovirus Infections in Immunocompetent and Immunocompromised Adult Patients. Epidemiol Infect (2020) 147:e328. doi: 10.1017/S0950268819002176

PubMed Abstract | CrossRef Full Text | Google Scholar

43. Yoon H, Jhun BW, Kim SJ, Kim K. Clinical Characteristics and Factors Predicting Respiratory Failure in Adenovirus Pneumonia. Respirology (2016) 21(7):1243–50. doi: 10.1111/resp.12828

PubMed Abstract | CrossRef Full Text | Google Scholar

44. Sun B, He H, Wang Z, Qu J, Li X, Ban C, et al. Emergent Severe Acute Respiratory Distress Syndrome Caused by Adenovirus Type 55 in Immunocompetent Adults in 2013: A Prospective Observational Study. Crit Care (2014) 18(4):456. doi: 10.1186/s13054-014-0456-6

PubMed Abstract | CrossRef Full Text | Google Scholar

45. Waye MMY, Sing CW. Anti-Viral Drugs for Human Adenoviruses. Pharmaceuticals (2010) 3(10):3343–54. doi: 10.3390/ph3103343

CrossRef Full Text | Google Scholar

46. Lion T. Adenovirus Infections in Immunocompetent and Immunocompromised Patients. Clin Microbiol Rev (2014) 27(3):441–62. doi: 10.1128/CMR.00116-13

PubMed Abstract | CrossRef Full Text | Google Scholar

47. Jose RJ, Manuel A. COVID-19 Cytokine Storm: The Interplay Between Inflammation and Coagulation. Lancet Respir Med (2020) 8(6):e46–7. doi: 10.1016/S2213-2600(20)30216-2

PubMed Abstract | CrossRef Full Text | Google Scholar

48. Guo XJ, Thomas PG. New Fronts Emerge in the Influenza Cytokine Storm. Semin Immunopathol (2017) 39(5):541–50. doi: 10.1007/s00281-017-0636-y

PubMed Abstract | CrossRef Full Text | Google Scholar

49. Biserni GB, Dondi A, Masetti R, Bandini J, Dormi A, Conti F, et al. Immune Response Against Adenovirus in Acute Upper Respiratory Tract Infections in Immunocompetent Children. Vaccines (Basel) (2020) 8(4):602. doi: 10.3390/vaccines8040602

CrossRef Full Text | Google Scholar

50. Schroder K, Hertzog PJ, Ravasi T, Hume DA. Interferon-Gamma: An Overview of Signals, Mechanisms and Functions. J Leukoc Biol (2004) 75(2):163–89. doi: 10.1189/jlb.0603252

PubMed Abstract | CrossRef Full Text | Google Scholar

51. Di Paolo NC, Shayakhmetov DM. Interleukin 1α and the Inflammatory Process. Nat Immunol (2016) 17(8):906–13. doi: 10.1038/ni.3503

PubMed Abstract | CrossRef Full Text | Google Scholar

52. Rojas JM, Avia M, Martín V, Sevilla N. IL-10: A Multifunctional Cytokine in Viral Infections. J Immunol Res (2017) 2017:6104054. doi: 10.1155/2017/6104054

PubMed Abstract | CrossRef Full Text | Google Scholar

53. Cao B, Gao H, Zhou B, Deng X, Hu C, Deng C, et al. Adjuvant Corticosteroid Treatment in Adults With Influenza a (H7N9) Viral Pneumonia. Crit Care Med (2016) 44(6):e318–28. doi: 10.1097/CCM.0000000000001616

PubMed Abstract | CrossRef Full Text | Google Scholar

54. Shang L, Zhao J, Hu Y, Du R, Cao B. On the Use of Corticosteroids for 2019-Ncov Pneumonia. Lancet (2020) 395(10225):683–4. doi: 10.1016/S0140-6736(20)30361-5

PubMed Abstract | CrossRef Full Text | Google Scholar

Keywords: human adenoviruses, cytokine storm, prediction, immunocompetent adults, disease progression

Citation: Li J, Wei J, Xu Z, Jiang C, Li M, Chen J, Li Y, Yang M, Gu Y, Wang F, Shu Y, Yang Y, Sun L and Liu Y (2021) Cytokine/Chemokine Expression Is Closely Associated Disease Severity of Human Adenovirus Infections in Immunocompetent Adults and Predicts Disease Progression. Front. Immunol. 12:691879. doi: 10.3389/fimmu.2021.691879

Received: 07 April 2021; Accepted: 20 May 2021;
Published: 07 June 2021.

Edited by:

Daniela Novick, Weizmann Institute of Science, Israel

Reviewed by:

Soraya Mezouar, Aix-Marseille Université, France
Laura Croce, University of Pavia, Italy

Copyright © 2021 Li, Wei, Xu, Jiang, Li, Chen, Li, Yang, Gu, Wang, Shu, Yang, Sun and Liu. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.

*Correspondence: Yingxia Liu, yingxialiu@hotmail.com

These authors have contributed equally to this work and share senior authorship

Disclaimer: All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article or claim that may be made by its manufacturer is not guaranteed or endorsed by the publisher.