Introduction

In early December 2019, several pneumonia cases of unknown origin were observed in the city of Wuhan in China. The causative pathogen was identified as a novel enveloped RNA virus belonging to the genus Betacoronavirus that was later named “severe acute respiratory syndrome coronavirus 2” (SARS-CoV-2).The virus showed genetic similarity to both SARS-CoV and Middle East respiratory syndrome coronavirus (MERS-CoV) [1]. Subsequently, it has spread all over the world and was declared a pandemic by the World Health Organization on March 11, 2020 [2]. There are reports of patients with SARS-CoV-2 infections having mild to severe respiratory illness with fever, cough, and shortness of breath as the common symptoms 2–14 days after exposure. However, the majority of patients diagnosed as positive using an RT-PCR test are reported to be either asymptomatic or minimally symptomatic after an incubation period of approximately 5–6 days [3]. Fever, dry cough, shortness of breath, myalgia, fatigue, leucopenia, and radiological signs of progressive pneumonia are the similar clinical and laboratory findings of disease progression seen in SARS-CoV-2, SARS-CoV, and MERS-CoV infections [2].

Severe pneumonia caused by coronaviruses has often been associated with rapid virus replication, infiltration of inflammatory cells, and elevated cytokines, resulting in damage to internal organs and acute respiratory distress syndrome (ARDS) [4]. Virus-induced aberrant and excessive cytokine production has been linked to morbidity and mortality in several viral infections [5]. Cytokine release syndrome (CRS), a systemic inflammatory response, can be caused by an infection, some drugs, or cancers and is characterized by a major increase in the level of a large number of pro-inflammatory cytokines. CRS was found to be the major cause of morbidity in patients infected with SARS or MERS-CoV [6]. In SARS-CoV infection, ARDS is the ultimate result of a cytokine storm in which large amounts of pro-inflammatory cytokines and chemokines released by immune effector cells sustain an aberrant systemic inflammatory response [1]. The observed peculiar ability of SARS-CoV to counteract the antiviral IFN response suggests that the severity of disease might be due to immune dysregulation, which could be characterized by an insufficient type I interferon response, paralleled by aberrant secretion of pro-inflammatory chemokines by alveolar macrophages, dendritic cells, and pneumocytes [1]. Data suggest that MERS-CoV infection shares certain immunological aspects with SARS-CoV infection in terms of involvement of the chemokine system [1]. Estimation of innate immune cell counts has indicated that activated dendritic cells, activated mast cells, and neutrophils are more abundant in patients with viral pneumonia than in healthy controls. A recent study has suggested the neutrophil-to-lymphocyte (NLR) ratio to be associated with the disease severity of COVID-19. In patients with ARDS, the alveolar spaces are reported to be occupied by infiltrating neutrophils and monocytes, suggesting a pathogenic role of those cells. SARS-CoV-2 triggers a robust IFN response, hallmarked by the expression of numerous ??interferon-stimulated genes?? (ISGs), thus indicating a role of the innate immune response in COVID-19 infection. The protective potential of ISGs may account for the lower proportion of severe cases and the lower case-fatality rate in COVID-19 compared to SARS [7]. There has been documentation of IFN treatment as an antiviral therapy against SARS-CoV-2 infection (ChiCTR: ChiCTR2000029600). However, the timing and the dose required need careful consideration.

It is believed that a dysregulated host immune response and production of inflammatory cytokines, known as a “cytokine storm”, correlate with disease severity and poor prognosis during both SARS-CoV and MERS-CoV infections [8]. Angiotensin-converting enzyme 2 (ACE2), an essential regulatory enzyme of cardiovascular function, has become of extreme interest during the SARS-CoV-2 pandemic, as reduction in ACE2 expression leads to the release of pro-inflammatory cytokines, ARDS, myocarditis, and hypercoagulability, with the possibility of exacerbation of acute coronary syndrome and induction of pulmonary embolism [9]. During SARS-CoV-2 infection, many patients develop ARDS, leading to pulmonary edema and respiratory failure, and these symptoms are associated with a cytokine storm [10]. The cytokine and chemokine profiles in SARS-CoV-2-infected patients with different clinical manifestations remain unclear. In the present study, we have investigated the plasma cytokine and chemokine levels in hospitalized asymptomatic (n = 39) and mildly symptomatic (n = 35) SARS-CoV-2-infected patients, recovered individuals 45-60 days after infection with SARS-CoV-2 (n = 40), and healthy control individuals (n = 36) from Pune, Maharashtra, India.

Materials and methods

Study population

The current study was carried out between March 2020 and June 2020. The prevention and control plan for the management of close contacts of individuals with COVID-19 have evolved, keeping the focus on emphasizing identification and quarantine of asymptomatic patients. To identify asymptomatic patients, reverse transcription polymerase chain reaction (RT-PCR) screening of throat and nasal swabs for SARS-CoV-2 for close contacts of index cases was done, and positive individuals were transferred to designated hospitals for centralized isolation [11]. Blood samples from asymptomatic (n = 39) and mildly symptomatic (n = 35) patients with SARS-CoV-2 infection were collected in the acute phase during hospitalization (≤6 days). The current study population consisted of a subset of the patients whose nasal/throat swabs tested positive for SARS-CoV-2 by RT-PCR to ICMR-NIV, Pune. They were admitted to the state-run Naidu Hospital and the YCM Hospital, which are tertiary-care hospitals in Pune, Maharashtra. Asymptomatic patients who developed symptoms within 6 days after admission were excluded from the study. Finally, 39 asymptomatic cases, defined by a SARS-CoV-2-positive nucleic acid test result but without any relevant clinical symptoms during hospitalization, were included in the current study. The symptoms of mildly symptomatic COVID-19 patients were cough, fever, fatigue, headache, myalgias, and diarrhoea. A total of 35 mildly symptomatic patients were selected for comparison with the asymptomatic individuals. All patients from both categories recovered. The recovered individuals (n = 40) included in the study had a previous history of symptomatic/asymptomatic SARS-CoV-2 infection and had become negative for SARS-CoV-2 by RT-PCR. Blood samples from recovered individuals were collected between 45 and 60 days post-diagnosis. The recovered individuals and the patient groups (both asymptomatic and mild symptomatic) were all distinct.

Thirty-six apparently healthy individuals from blood donation camps organized by Sassoon General Hospital, Pune, in December 2019, were recruited as healthy controls for cytokine comparison (Table 1). The healthy controls were screened for SARS-CoV-2 by ELISA. Plasma samples were used to detect the presence of IgG antibodies using a commercial ELISA (COVID Kavach-Anti-SARS-CoV-2 IgG Antibody Detection ELISA, M/s Cadila Healthcare Limited, Ahmedabad). Only those who were negative for anti-SARS-CoV-2 IgG antibody were included in the study.

Table 1 Characteristics of the study population

The study was approved by the Institutional Ethical Committee for Research on Humans, based on the guidelines set by the Indian Council of Medical Research, New Delhi. Informed written consent was obtained from all study participants.

Plasma separation and storage

About 3 to 4 ml of blood from each study participant was collected in a K3 EDTA tube, and the plasma was separated within 4 hours of collection by the Ficoll-Hypaque (Sigma, USA) density gradient centrifugation method and stored at −80 °C. Plasma samples were thawed at the time of estimation of cytokine, chemokine, and growth factors levels.

Estimation of cytokine, chemokine, and growth factor levels

Plasma concentrations of cytokines, chemokines, and growth factors were determined using a Bio-plex Multiplex Immunoassay System (Bio-Rad, Hercules, CA, USA) using a Bio-plex ProTM Human Cytokine 27-plex assay kit as reported previously [12, 13] as per the manufacturer’s instructions. Levels of 15 cytokines, including the pro-inflammatory (IL-1β, IL-5, IL-6, IL-7, IL-9, IL-15, IL-17, TNF-α), anti-inflammatory (IL1-RA, IL-4, IL-10, IL-13), and Th1 (IL-2, IFN-γ, IL-12 p70) cytokines along with seven chemokines (eotaxin, CCL-2,CCL-3,CCL-4,CCL-5, IL-8, CXCL-10) and five growth factors (basic fibroblast growth factor [FGF], G-CSF, GM-CSF, vascular endothelial growth factor [VEGF], platelet-derived growth factor-bb [PDGF-bb]), were estimated.

Statistical analysis

Levels of all analytes were analysed after log10 transforming the observed concentrations of individual cytokines/chemokines/growth factors. A value of 0.1 pg/ml was used in cases where the concentrations of the cytokines, chemokines, and growth factors in the tested samples were too low to be detected. The Mann–Whitney U test was used for numerical data for comparisons among the study groups. A p-value of <0.05 was considered significant. Statistical analysis was performed using IBM SPSS Statistics 25 software (SPSS Inc., Chicago, IL, USA). A global-view heat map was generated using the heatmap.2 package in R. Receiver operating characteristic analysis (ROC) was performed using GraphPad Prism 8 software (GraphPad, San Diego, CA, USA).

Results

Characteristics of the study population

The characteristics of the study groups are presented in Table 1. The patients investigated in this study included (a) 74 patients with SARS-CoV-2 infection (45 males, 29 females) with age ranging from 15 to 79 years, which included 39 asymptomatic and 35 mildly symptomatic patients, (b) 40 recovered individuals, and (c) 36 healthy controls. The patients did not have any coinfections with other respiratory infections and had an uneventful recovery. The patients in the acute phase of infection were positive for SARS-CoV-2 by qRT-PCR [11, 14], while the recovered individuals had a reported history of SARS-CoV-2 infection. The clinical features of the mildly symptomatic patients are described in Table 2.

Table 2 Clinical features of mildly symptomatic patients

Peripheral cytokines and chemokines in asymptomatic and mildly symptomatic patients with SARS-CoV-2 infection

IL-1β, IL-6, IL-7, IL-9, TNF-α, IL1-RA, IL-4, IL-13, CCL-4 (MIP-1α), and CXCL-10 (IP-10) were significantly elevated in both asymptomatic and mildly symptomatic patients compared to controls (p < 0.05). Importantly, the levels of IL-17, eotaxin, and PDGF-BB were significantly elevated in asymptomatic patients only, when compared to controls. IL-5, IFN-γ, IL-12 p70, CCL-2, CCL-3, IL-8, and VEGF were significantly lower in both asymptomatic and mildly symptomatic patients than in controls (p < 0.05) (Table 3). Importantly, the level of G-CSF was significantly lower only in mildly symptomatic patients when compared to controls. Levels of IL-15, IL-10, IL-2, CCL-5 (RANTES), FGF, and GM-CSF were comparable among all groups (Table 3 and Figs. 1 and 2).

Table 3 Levels of cytokines/chemokines in patients with recent SARS-CoV-2 infection and in healthy controls
Fig. 1
figure 1

Global heat map of levels of cytokines/chemokines in SARS-CoV-2-infected patients, recovered individuals, and healthy controls

Fig. 2
figure 2

Comparison of concentrations of cytokines and chemokines expressed in Log10(pg/mL) in different study groups. a IL-6, b TNF-α, c CXCL-10, d IL-4, e IL-17. *p <0.05; **p < 0.005; ***p < 0.0001

Comparison among the disease categories revealed no difference in the levels of the chemokines/cytokines except for (a) CXCL-10 being significantly higher and (b) IL-17, IL-4, and VEGF being lower in mildly symptomatic patients (p < 0.01) (Table 3 and Figs. 1 and 2).

Peripheral cytokines and chemokines, post-recovery

Among the patient groups

IL-1β, IL-6, IL-7, IL-9, IL-17, TNF-α, IL1-RA, IL-4, IL-10, IL-13, IFN-γ, eotaxin, CCL-2, CCL-3, CCL-4, CCL-5, IL-8, CXCL-10, G-CSF, and PDGF-BB were significantly lower in recovered individuals than in both asymptomatic and mildly symptomatic patients. However, Basic FGF was significantly lower only in the recovered group than in the asymptomatic group, while VEGF was significantly lower only in the recovered group than in the mildly symptomatic group. Levels of IL-5, IL-15, IL-2, IL-12 p70, and GM-CSF were comparable among the recovered, asymptomatic, and mildly symptomatic groups (Table 3).

Recovered vs. control groups

IL-5, IL-9, IL-17, IL1-RA, IL-10, IFN-γ, IL-12p70, eotaxin, CCL-2, CCL-3, CCL-4, CCL-5, IL-8, basic FGF, G-CSF, VEGF, and PDGF-BB were significantly lower in the recovered group than in the control group. Levels of IL-1β, IL-6, IL-7, IL-15, TNF-α, IL-4, IL-2, and GM-CSF were comparable among the recovered and control groups. Importantly, the level of CXCL-10 was significantly higher in the recovered group than in the control group (Table 3 and Fig. 2).

ROC analysis of cytokines

To generate a cutoff for the analytes among the patient categories, ROC curve analysis was performed (Fig. 3 and Table 4). The ROC analysis generated cutoffs for individual cytokines/chemokines that could be used to distinguish healthy controls from patients with SARS-CoV-2 infection and might have future diagnostic importance.

Fig. 3
figure 3

Receiver operating characteristic (ROC) curves for IL-1β, TNF-α, CXCL-10, and IL-4, validating their applicability as biomarkers of recent infection. a ROC curve generated by comparing concentrations of TNF-α (pg/mL) in SARS-CoV-2-infected patients with their concentrations in healthy controls reveals the inefficiency of TNF-α as a biomarker (cutoff, 29.19; sensitivity, 98.65%; specificity, 83.33%), b ROC curve generated by comparing concentrations of CXCL-10 (pg/mL) in SARS-CoV-2-infected patients with their concentrations in healthy controls reveals the efficiency of CXCL-10 as a biomarker (cutoff, 124.9; sensitivity, 97.3%; specificity, 91.67%), c ROC curve generated by comparing concentrations of IL-1β (pg/mL) in SARS-CoV-2-infected patients with their concentrations in healthy controls reveals the efficiency of IL-1β as a biomarker of infection (cutoff, 0.545; sensitivity, 91.89%; specificity, 86.11%). d ROC curve generated by comparing concentrations of IL-4 (pg/mL) in SARS-CoV-2-infected patients with their concentrations in healthy controls reveals the efficiency of IL-4 as a biomarker of infection (cutoff, 1.395; sensitivity, 97.30%; specificity, 66.67%)

Table 4 ROC characteristics of the cytokines TNF-α, CXCL-10, IL-1β, and IL-4 in plasma of SARS-CoV-2-infected patients and healthy controls

The ROC characteristics of CXCL-10 included a cutoff value of 124.9 pg/mL (sensitivity, 97.30%; specificity, 91.67%). Similarly, the ROC characteristics of TNF-α, IL-1β, and IL-4 showed cutoff values of 29.19 pg/mL, 0.545 pg/mL, and 1.395 pg/mL, respectively (sensitivity, 98.65%, 91.89%, and 97.30%, respectively; specificity, 83.33%, 86.11%, and 66.67%, respectively). The AUC values of the above analytes were greater than 0.9 (p <0.0001), which is indicative of higher diagnostic value. The cytokines IL-1β, IL-4, and TNF-α and the chemokine CXCL-10 were found to be predictive of early SARS-CoV-2 infection, which was confirmed by ROC analysis (AUC = 0.9917, 0.9202, 0.9336, and 0.9508, respectively). The high sensitivity and specificity of the analytes’ cutoffs suggest that they can potentially act as a biomarker of SARS-CoV-2 infection.

Discussion

A good proportion of patients with SARS-CoV-2 infection are reported to have mild to moderate respiratory illness. However, in parallel, other patients (identified mostly through contact tracing) remain asymptomatic. Information such as epidemiological characteristics, virus levels, and immune responses in this subset of patients is scarce. To gain a better understanding of the host factors modulating the course of the infection in individuals with SARS-CoV-2 infection, we have compared the levels of plasma cytokines, chemokines, and growth factors in a group of hospitalized asymptomatic and mildly symptomatic patients with SARS-CoV-2 infection and in recovered individuals.

Comparative analysis of cytokine levels in MERS-CoV-infected patients and normal healthy controls had shown a Th2-skewed immune response in infected patients, which was somewhat unique for MERS-CoV infection when compared to SARS-CoV, which has been shown to induce a Th1-type immune response [15]. The asymptomatic and mildly symptomatic patient groups from the current study did not exhibit any such specific pattern. The observations of the present study and published data on the cytokine profiles in acute MERS-CoV and SARS-CoV infections appear to be conflicting. It is possible that SARS-CoV-2 may have evolved a different mechanism to mitigate host defences.

Infection with coronaviruses results in the activation of monocytes, macrophages, and dendritic cells, leading to the secretion of prominent pro-inflammatory cytokines such as IL-6, along with other inflammatory cytokines [16]. Elevated IL-6 can activate either classical cis-signaling or trans-signaling, leading to cytokine release syndrome/cytokine storm [16]. It is understood that the timing of IL-6 activation is critical, as it is necessary for the early immune response for fighting SARS-CoV-2 and other viral infections. Low IL-6 levels in the early stages of infection could result in uncontrolled viral replication. Thus, rather than blocking IL-6 release, regulating it would have therapeutic value [17]. Several studies have found elevated serum concentrations of inflammatory cytokines, including interleukin IL-6, in cases of severe COVID-19. However, these elevated levels of inflammatory cytokines, including IL-6, were much lower than those reported in patients with ARDS unrelated to COVID-19, such as sepsis and chimeric antigen receptor (CAR) T-cell-induced cytokine release syndrome [18] The dynamic role of IL-6 in the cytokine storm and CRS in the context of SARS-CoV-2 infection has been elucidated by Zhang et al. [19]. It is believed that IL-6 blockade is a promising strategy for treatment of COVID-induced CRS, as elevated IL-6 levels were consistently reported in several studies of COVID-19 and might serve as a predictive biomarker for disease severity [20]. Treatment of the cytokine storm has become an imperative part of rescuing severe COVID-19 patients, as mild or severe cytokine storms are reported in severe cases. For this purpose, the IL-6R antagonist tocilizumab has been suggested as an important drug to save patients' lives [19]. On similar lines, regulated higher levels of IL-6 in both asymptomatic and mildly symptomatic patients in the current study may be correlated with disease outcome. The lower IL-6 levels in the recovered individuals compared to both patient groups could be used as a marker of recovery, which goes hand in hand with the reported beneficial effects of IL-6 blockade therapy [8].

Chemokines are known to play a role in the cytokine storm related to coronavirus infection [1]. In an analysis of 23 patients with SARS infection, Jiang et al. showed that serum CXCL-10 was markedly elevated during the early stage of SARS and remained at a high level until convalescence [21]. Analysis of serum cytokines and chemokines in patients with MERS-CoV infection revealed that levels of IL-6 and CXCL-10 were significantly elevated in MERS-CoV patients who developed severe disease [22]. A recent report indicated that increased production of pro-inflammatory chemokines, mainly CXCL-10 and CXCL-8, is a characteristic of COVID-19. Interestingly, CXCL-8 transcription is upregulated by SARS-CoV infection but not by SARS-CoV-2 infection, while the opposite has been observed for CXCL-10 [1]. This coincides with the current study’s finding of elevated CXCL-10 in both patient categories. The levels of CXCL-10 have been reported to be high in SARS-CoV-2-infected patients with severe disease compared to healthy controls, indicating that the levels of CXCL-10 increase as the disease progresses in severity [5]. CXCL10/IP10 and MCP-1 have also been reported as biomarkers associated with the severity of COVID-19 and can be related to the risk of death in COVID-19 patients. Studies have indicated CXCL10 to be associated with thrombosis and responsible for inhibition of endothelial healing. Since the serum CXCL10 level correlates directly with the severity of COVID-19, anti-CXCL10 antibody treatment might represent a new regimen for COVID-19 patients, especially those with thrombotic events. Anti-CXCL10 antibody is reported to be in a clinical trial to evaluate its ability to prevent cardiovascular events [23]. In the same line, the current study suggests that CXCL-10 could be associated with pathogenesis of SARS-CoV-2 infection. At the same time, recovered individuals at 45-60 days post-onset had higher levels of CXCL-10 than the controls, suggesting that the immune system of SARS-CoV-2 patients may take a longer time to normalize.

SARS CoV-2 infection is associated with the depletion of antiviral defences as well as enhanced production of pro-inflammatory cytokines [24]. Our patient population exhibited a similar scenario, with downregulation of antiviral cytokines and higher levels of pro-inflammatory cytokines. Lower levels of IFN-γ, the key antiviral cytokine, could be attributed to the observed lower levels of IL-12 p70, as reported elsewhere [25]. A study carried out in a Chinese population showed that IFN-γ, IL-12, and the inflammatory cytokines IL-1β and IL-6 can induce hyper-innate inflammatory response due to invasion of the respiratory tract by SARS-CoV, leading to the activation of Th1-cell-mediated immunity by the stimulation of NK and cytotoxic T lymphocytes (CTLs) [26]. Our study patients exhibited increased levels of the inflammatory cytokines IL-1β, IL-6, and TNF-α, which got normalized during the recovery phase with the activation of the CTL response (our unpublished data).

Comparison of the immunological profiles of children with symptomatic and asymptomatic SARS-CoV-2 infection indicated that those with symptomatic infection had higher IL-6, IL-10, TNF-α, and IFN-γ levels [27]. Comparison of serum cytokine and chemokine levels in asymptomatic, mild, moderate, severe, and convalescent SARS-CoV-2 infection cases in China showed that serum concentrations of IL-1RA, IL-1β, IL-6, and IP-10 were lower in asymptomatic cases than in symptomatic cases [5]. However, in the current study, comparison of asymptomatic patients and those with mild symptomatic disease revealed no difference in the levels of chemokines/cytokines, except for CXCL-10 being significantly higher and IL-17, IL-4, and VEGF being lower in mildly symptomatic patients, suggesting that these molecules could be used as biomarkers of the clinical severity of SARS-CoV-2 infections in Indian patients, which may help to evaluate the disease severity and outcome of COVID-19 patients. The overall variations in the cytokines TNF-α, IL-1β, IL-4, and CXCL-10 among the patients and healthy controls in the current study suggest their applicability as biomarkers of recent SARS-CoV-2 infection (Fig. 3 and Tables 3 and 4).

The following limitations of the current study should be taken into consideration: A major limitation was that we could not analyse immunological responses in patients with severe COVID-19. It is possible that the differences between the asymptomatic and mildly symptomatic patients would have been more pronounced, as reported elsewhere [3, 5].

In conclusion, the immunological reaction triggered by infection with SARS-CoV-2 in Indian patients is mainly of pro-inflammatory character and could be due to a reduced antiviral cytokine response against a background of CXCL-10 chemokine upregulation.