Abstract

Objective. To compare the difference of inflammatory cytokines and lymphocyte subsets between deceased patients and survivors with COVID-19. Methods. This retrospective study included 254 confirmed patients from 10 January to 11 March, 2020, at Tongji Hospital of Wuhan, China. Laboratory and immunologic features were collected and analyzed, and the main outcomes focused on inflammatory cytokines and lymphocyte subsets. Results. A trend of markedly higher levels of inflammatory cytokines as well as lower lymphocyte subset levels in deceased patients was observed compared with survivors. ROC curve analyses indicated that inflammatory cytokines and the decrease levels of T cell, Th (helper T cells) cell, Ts (suppressor T cells) cell, B cell, and NK cell along with Th/Ts ratio increase could be used to predict the death of COVID-19. Multivariate analyses showed that higher levels of IL-6, IL-8, and IL-10 remained significantly correlated with shorter survival time and that the amount of Ts cells was negatively associated with the possibility of death in COVID-19 patients. In conclusion, SARS-CoV-2 would cause lymphopenia and result in decreased lymphocyte subset cells, particularly in Ts cell counts, which further induces hyperinflammatory response and cytokine storm. IL-6, IL-8, IL-10, and Ts cell might be independent predictors for the poor outcome of COVID-19.

1. Introduction

SARS coronavirus 2 (SARS-CoV-2) is a member of β-coronavirus (β-CoV) lineage B belonging to CoV family, which has been confirmed for causing coronavirus disease 2019 (COVID-19) first identified and confirmed in Wuhan, China, in December 2019 [1, 2]. Genomic analysis indicated that SARS-CoV-2 shares about 80% identical genome and almost all encoded proteins with human severe acute respiratory syndrome (SRAS-CoV), while only 50% genomic identity with middle east respiratory syndrome (MERS-CoV) and more similarities to bat CoVs [3, 4]. Considering the highly contagious nature from person-to-person and global epidemic of SARS-CoV-2, WHO firstly declared it as a Public Health Emergency of International Concern at January 30th, 2020 [5], and then characterized it as a pandemic at March 11th, 2020 [6]. Currently, there are more than 50 COVID-19 vaccine candidates in clinical trials globally, and some vaccines have been approved for emergency use in United Kingdom, China, Arab nations (the United Arab Emirates, Bahrain, and Egypt), Republic of Chile, Canada and US, and so on [7].

Previous studies demonstrated that SARS-CoV-2 infection would cause severe respiratory symptoms, lymphopenia, increased levels of peripheral inflammatory cytokines, and cytokine storm which could be related with high severity and mortality of COVID-19 [812]. Besides, lymphocyte subsets play an important functional role in the maintenance of humoral and cytotoxic immunity [13]. Recently, several studies had investigated the association of immunological features and lymphocyte subset alteration with the disease severity in COVID-19, whose results found SARS-CoV-2 might have a main impact on T lymphocytes, especially helper T (Th) and suppressor T (Ts) cells [1315]. Decreased Ts cells might be an independent biomarker for severe COVID-19, and decreased Ts and B cells as well as increased Th/Ts ratio might be the independent predictors of poor treatment efficacy [13]. However, there is still insufficient knowledge about inflammatory cytokines, lymphocyte subsets, and immune response associated with the death of COVID-19.

In our study, we mainly aimed to compare and elucidate the difference and clinical significance of inflammatory cytokines and lymphocyte subsets between deceased COVID-19 patients and survivors, which might help further demonstrate the immunological effect in deceased COVID-19.

2. Material and Methods

2.1. Study Design and Population

We retrospectively recruited 254 laboratory-confirmed COVID-19 patients from 10 January to 11 March, 2020, at Tongji Hospital of Wuhan, China. All patients in the study have been confirmed with no infection of HBV, HCV, and HIV. Ethical Committee of Tongji Hospital of Huazhong University of Science and Technology (No.TJ-IRB20200364) and China-Japan Union Hospital of Jilin University (No.2020032607) approved the study and waived the written informed consent for rapid emerging infections.

2.2. Data Collection

All cases were confirmed on admission by nucleic acid and antibody tests. Information from electronic medical records was retrospectively collected and extracted, including demographic characteristics, medical history, routine blood and urine tests, biochemical detection, and inflammatory cytokines and lymphocyte subset tests. All data were extracted with the unified data collection form by two independent researchers. The severity of patients was classified by the Novel Coronavirus Pneumonia Diagnosis and Treatment Guidance (Seventh revised trial version) [16].

2.3. Statistical Analysis

Categorical variables were described as percentages and compared by chi-squared or Fisher’s exact test; continuous variables were described as means (SD) or medians (interquartile range (IQR)) and compared by nonparametric comparative tests or -tests. Spearman rank correlation coefficient was used to analyze the correlations between paired data. Area under the receiver operating characteristic curve (AUROC) was performed to evaluate the predictive values of peripheral inflammatory cytokines and lymphocyte subsets in predicting the mortality of COVID-19. Kaplan-Meier analyses with Log-rank test were performed to assess the clinical significance for survival time between high and low levels of inflammatory cytokines. Multivariate Cox regression model was performed to determine the predictive factors for deceased COVID-19. All statistical analyses were accomplished by the SPSS 22.0 and R 3.6.1 software, and results with were considered statistical significant.

3. Results

3.1. Demographics and Baseline Characteristic of Deceased Cases of COVID-19

254 patients were included with 51 (20.1%) deceased. The median age of all cases was 67 years (IQR, 58-73) and 70 years (IQR, 64-76) for deceased patients. Male patients were 130 (51.2%) including 32 (62.7%) deceased cases. The occurrence of any comorbidities in all patients reached 57.4% and 58.8% in deceased patients. Hypertension (39.6%) and diabetes (19.8%) were the most common comorbidities. Fever (85.1%), cough (85.1%), dyspnea (66.0%), expectoration (64.4%), myalgia or fatigue (63.4%), and chest tightness (53.9%) were most common symptoms; less common manifestations included nausea or vomiting (29.7%), sore throat (29.7%), abdominal pain (27.7%), and hemoptysis (27.5%). The median of hospitalization time was 22 days (IQR 12-29) in all patients, 9 days (IQR 5-16) in deceased patients, and 25 days (IQR 16-30) in survivors.

3.2. Laboratory Findings of Deceased COVID-19 Patients on Admission

Compared with normal range, levels of leucocytes and neutrophils in deceased patients were significantly elevated above upper limit of normal range (Table 1). The counts of lymphocytes were all significantly decreased below lower limit of normal range in deceased patients. Leukocytes, neutrophils, lymphocytes, and platelets were decreased in 14 (5.5%), 10 (3.9%), 151 (59.4%), and 38 (15.0%) patients. The number of patients who developed lymphopenia in deceased cases was 49 (96.1%). Counts of leucocytes, neutrophils, monocytes, and neutrophil-to-lymphocyte ratio (NLR) in deceased patients were markedly higher than survivors (all ), while levels of lymphocytes and platelets were all significantly lower in deceased patients (Table 1).

Levels of Glu, ALT, AST, ALP, and GGT were all significantly elevated in deceased cases (all ) (Table 1). TP, albumin, and TC levels were all markedly lower in deceased cases (all ). LDH, CRP, D-dimer, NT-proBNP, PCT, and cTnI concentrations were all markedly increased in deceased cases (all ) (Table 1).

In routine urine tests, the positive rates of urine protein, urine ketone body (KET), and urobilinogen (URO) were 54.6%, 16.4%, and 5.3%, respectively (Table 1). Compared with survivors, those positive rates of deceased patients were significantly higher (all ). Urine specific gravity (USG) and pH levels of all patients were within the normal range. However, the level of USG in deceased cases was significantly elevated (all ) while pH levels decreased ().

Blood gas analyses showed that the median arterial partial pressure of oxygen level in deceased patients was 68 mmHg, and the median of partial pressure of carbon dioxide was 31.9 mmHg. Total partial pressure of carbon dioxide, arterial partial pressure of oxygen, standard bicarbonate, base excess of extracellular fluid, total carbon dioxide, and oxygen saturation were significantly lower in deceased patients than in survivors (all ) (Table 1).

3.3. Peripheral Inflammatory Cytokines and Lymphocyte Subset Levels with Deceased COVID-19

The level of interleukin-1β (IL-1β) was undetectable (<5 pg/mL) in almost all patients so that no statistical significance was observed between groups (Table 2). Levels of interleukin-2R (IL-2R), interleukin-6 (IL-6), interleukin-8 (IL-8), interleukin-10 (IL-10), and tumor necrosis factor α (TNF-α) on admission were significantly elevated in 120 (47.2%), 136 (53.5%), 20 (7.9%), 49 (19.3%), and 93 (36.6%) patients (Table 2). Compared with survivors, these levels were significantly higher in deceased patients (all ) (Table 2, Figure 1).

Total counts of T cells and suppressor T (Ts) cells were both below lower limit of normal ranges, while counts of B cells, helper T (Th) cells, Ts cells, NK cells, and Th/Ts ratio were within the normal limits (Table 2). Counts of lymphocytes and T cell subsets were significantly decreased in deceased cases (all ). Th/Ts ratios were significantly elevated beyond upper limit of normal range (3.32 vs. 2.00/μL, ). Counts of T, B, Th, Ts, and NK cells were decrease in 71/137 (51.8%), 39/137 (28.5%), 65/137 (47.4%), 85/137 (62.0%), and 64/137 (46.7%) cases. Th/Ts ratio was increased in 29 (21.2%) patients. Meanwhile, counts of T, B, Th, Ts, and NK cells were markedly lower in deceased cases, but Th/Ts ratio was higher (all ) (Table 2, Figure 2).

3.4. Correlation Analysis between Cytokine Profiles and Lymphocyte Subsets in COVID-19

Total counts of T, Ts, B, and NK cells were negatively related with levels of LDH, CRP, IL-2R, IL-6, IL-8, IL-10, and TNF-α (all ). Except for TNF-α (), counts of Th cells was negatively correlated with numbers of neutrophils () and platelets () and concentrations of CRP (), LDH (), CRP (), IL-2R (), IL-6 (), IL-8 (), and IL-10 (). Total numbers of T, Ts, and NK cells were also negatively correlated with neutrophil counts (all ), and Ts, B, and NK cells with leucocytes (all ). Numbers of B and NK cells were negatively correlated with glucose level (). NK cell count was also negatively correlated with numbers of monocytes () and platelets () (Figure 3). Details of the correlation between all the collected data are shown in Supplementary file see (available here).

3.5. Predictive Values of Potential Markers for Distinguishing Deceased COVID-19

ROC curve analyses indicated that the plasma levels of inflammatory cytokines could predict the mortality of COVID-19 as all AUCs were above 0.7 (Figure 4(a)). ROC curve analyses in lymphocyte subsets indicated that the decrease of T, Th, Ts, B, and NK cell counts along with Th/Ts ratio increase could be used to predict the occurrence of death in COVID-19 patients (all ) (Figure 4(b)).

Kaplan-Meier survival curves showed that higher IL-2R, IL-6, IL-8, IL-10, and TNF-α levels were significantly correlated with reduced survival time in COVID-19 (all ) (Figure 5).

Multivariate Cox model analyses showed that higher levels of IL-6, IL-8, and IL-10 remained significantly correlated with shorter survival time (IL-6: HR 26.22, ; IL-8: HR 3.11, ; IL-10: HR 2.94, ), and Ts cell count (/μL) was negatively associated with COVID-19-related death (HR 0.98, ) (Figure 6(b)), indicating that IL-6, IL-8, and IL-10 levels and Ts cell count could be the independent risk factors to predict the possibility of death of COVID-19.

4. Discussion

Previous studies suggested that several coronavirus infections could cause sustained responses of cytokine storm and dysregulation in lymphocyte subset levels, resulting in high incidence of immune disorders and mortality [13, 17]. Although a number of researches on clinical and epidemiological characteristics have been published [9, 11, 18], there is still insufficient knowledge of the potential mechanism in COVID-19. Thus, elucidating the features of peripheral inflammatory cytokines and lymphocyte subsets would be important for exploring the immune response mechanism of COVID-19.

Our laboratory findings showed that most patients presented lymphopenia, increased neutrophil count, and elevated levels of glucose, NT-proBNP, procalcitonin, cTnI, ALT, AST, LDH, CRP, and D-dimer. In deceased cases, there were significantly higher numbers of leucocytes and neutrophils and markedly higher levels of NLR, glucose, NT-proBNP, procalcitonin, cTnI, ALT, AST, LDH, CRP, and D-dimer than survivors but significantly lower numbers of lymphocytes, monocyte, and platelets. Leukocytosis, lymphopenia, and increased concentration of acute phase reactants were more commonly seen in deceased cases, which could promote death from COVID-19. Patients with uncontrolled blood glucose or diabetes had longer median length of stay and significantly higher mortality than those with well-controlled blood glucose or nondiabetes [1921]. The increased neutrophil count, NLR, and D-dimer might indicate higher incidence of poor clinical outcomes [10, 14, 22] and myocardial injury which might be induced by inflammation was related with cardiac insufficiency, arrhythmias, and adverse endpoints [23]. Recently, one study using machine learning tools has identified 3 indicators (LDH, lymphocyte, and hs-CRP) which could predict the death of COVID-19 more than 10 days in advance with the accuracy of over 90% and found that high levels of LDH require immediate medical attention in patients [24]. Previous studies had demonstrated that acute phase reactants play key roles in the diagnosis and prediction of systemic inflammation and infection as well as the prognosis in severe patients of SARS and MERS [2527], which were consistent with our findings. Further research on acute phase reactants would help with early diagnosis, prediction of disease process, and clinical outcomes.

Our results showed that all inflammation cytokine levels were substantially elevated, which is characterized as cytokine storm, to be significantly higher in most of the deceased cases similar with SARS and MERS results [28, 29]. Previous studies had suggested that inflammation cytokines and chemokines hold vital positions in cytokine storm [30]. Elevated concentrations of IL-6, IL-8, and IL-10 remained significantly correlated with shorter survival time on multivariate analyses, indicating that the association between baselines of IL-6, IL-8, and IL-10 and shorter survival time in COVID-19 might be the independent predictor of poor clinical outcomes and an addition to prognostic factors. Elevated concentration of IL-6 had been indicated as a stable indicator for the progression and adverse endpoint of COVID-19, and IL-6 and IL-10 were significantly related with the positivity duration of SARS-CoV-2 [31, 32], which were consistent with our findings. Previous studies showed that IL-6R antagonist tocilizumab could reverse the cytokine storm [33], and specifically blockading the signaling pathways regulated by IL-6 might be a promising method to reduce inflammation-related injuries and the mortality [34, 35]. IL-10 has a major anti-inflammatory effect regulated by JAK-STAT pathway [36], produced by antigen presentation cells (APC) including dendritic cells, macrophages, and Th cells [36, 37]. Higher IL-8 levels were found in deceased cases that might induce the hyperinflammatory response, while only S protein in all SARS proteins could induce the activation of IL-8 promoter [38, 39]. Anti-Spike-IgG (Anti-S-IgG) in SARS-CoV-2 infection would promote proinflammatory monocyte or macrophage accumulation and IL-8 production in lungs [40]. Additionally, a single center, randomised, open-label, phase 2 trial of anti-IL-8 therapy (BMS-986253) compared with Standard of Care in treating Hospitalized severe patients with COVID-19 (NCT04347226) is ongoing for improving the health condition of patients.

Hyperinflammatory cytokines may be associated with tissue injury, respiratory failure, multiple organ failure, or shock [31]. Besides, cytokines released under the innate immune responses against virus infection could also induce the release of glucocorticoids and other peptides by the neuroendocrine system, resulting in impaired immune response [31]. SARS-CoV-2 viral proteins might target several innate immune signaling proteins from pathways including the interferon (IFN) pathway, NF-κB pathway, and JAK/STAT signaling pathway [15, 41] that play crucial roles in excessive inflammation and immunopathology against virus infection, causing increased inflammatory cytokines/chemokines levels. Additionally, previous studies indicated that therapies using cytokines or cytokine inhibitors have been successively successful [42], and interferon α of antipathogenic cytokines is an officially recommended drug for treating COVID-19 in China [43]. The dysregulation of immune response would induce the hyperinflammatory even cause mortality. However, the underling role of cytokine storms in COVID-19 pathogenesis is still unclear and needs to further research.

Direct cytopathic effects and immune evasion have been proved to play important roles in disease severity of coronavirus infection [17]. Here, we demonstrated that lymphocytes and T cell subset counts showed significant reduction in most cases of COVID-19, but B cell counts decreased only in 28.5% of them. T cell and Ts cell counts decreased beyond the lower limit of normal range, while numbers of Th, B, and NK cells remained within normal range. In deceased cases, all lymphocytes and T cell subset counts reduced markedly below the lower limit of normal range, which were significantly lower than survivors. One previous study suggested that lymphocyte subset levels in severe cases had significantly reduced, excluding only B cells [15] which was different with our results. This might due to the fact that there were only 6 out of 8 severe cases in their study and no available immunological data in moderate cases. In our study, Th/Ts ratio elevated in 21.2% patients and was markedly higher in deceased cases than survivors. The reduction in counts of all lymphocyte subsets was correlated with both severity and mortality of COVID-19 [44], but Th/Ts ratio only correlated with mortality. However, only Ts cell was found significantly associated with deceased COVID-19 on multivariate analyses, indicating that Ts cell might be the potential predictor for deceased COVID-19, which was consistent with previous studies [13, 44]. We also found that counts of lymphocyte subsets were all negatively correlated with IL-6, IL-8, and IL-10 levels. Previous studies indicated that IL-6 could suppress normal T cell activation while IL-8 would attract neutrophil granulocytes and T lymphocytes, which might partly explain the lymphopenia [38, 39]. Besides, IL-6 and IL-8 induced by SARS-CoV would impair T cell functions on prime dendritic cells as well as viral clearance of macrophages and dendritic cells, thus causing impairment of adaptive immune response [33]. As for IL-10, it could indirectly inhibit T cell activation and production of T cell-derived cytokines by the inhibitory effects on APCs as well as directly suppress production of some cytokines [45]. IL-10 could also mediate direct effects on T cells, B cells, NK cells, mast cells, and eosinophils [45]. It was reported that IL-10 signaling blockade at the priming stage could induce stronger Ts cell responses, which is beneficial for treatment of chronic viruses and cancer [37]. Our study demonstrated that the depletion of lymphocyte subsets might cause hyperinflammatory response in deceased cases, resulting in overproduction of proinflammatory cytokines and cytokine storm in SARS-CoV-2 infection, which might be part of the immunogenicity of COVID-19.

Mechanism of lymphocyte reduction in deceased COVID-19 is still unclear. Previous studies revealed that lymphopenia might be correlated with apoptosis activation and P53 signaling pathway in lymphocytes [46]. Similar with SARS-CoV and MERS-CoV [47, 48], these infections could induce apoptosis in T cells and cause significantly decrease in numbers of peripheral lymphocyte subsets, which would further induce apoptosis in multiple organs and T lymphocytes. Moreover, different exposure doses of SARS-CoV-2 could also affect lymphocyte response. High dose exposure would cause adverse endpoints and viral clearance delay, which might be related to the inefficient T and B cell responses caused by lymphopenia, subsequently inducing cytokine storm and destructive tissue inflammation [41]. It is believed that lymphocytes, especially Th and Ts cells, serve major roles in attenuating overactivated innate immune response in SARS-CoV-2 infection [15, 49]. In this process, Th cell helps promoting virus-specific antibodies production, regulating immune responses, and mediating the activation, proliferation, and deletion of immune cells (especially Ts cells), which may result in the finding that high-level Th cell count had a strong correlation with the severity of SARS-CoV infection [47]. Besides, Ts cells are critical for regulating viral clearance and immune-mediated injury, which can be strongly affected by lymphopenia; Ts and NK cells are essential for facilitating appropriate antiviral responses, which might be associated with the elevated level of NK group 2 member A (NKG2A) [49, 50]. Glucocorticoid treatment in COVID-19 might induce the decrease in lymphocytes [13, 15], and drugs on lymphocyte proliferation or preventing apoptosis might be beneficial for preventing lymphopenia or restoring lymphocyte numbers in severe cases [51].

Several limitations may exist in the study. Firstly, this is a retrospective study in which we only estimated the number of lymphocyte and T cell subsets on admission with severity and mortality in COVID-19. The percentage, temporal change, and function of lymphocyte subsets together with the dynamics of cellular immune response after SARS-CoV-2 infection might need to be elucidated by the large-sample study. Secondly, there is an impairment of immune system after SARS-CoV-2 infection that might cause coinfection with various pathogens or secondary infection, possibly influencing the effect of immune response in COVID-19.

In conclusion, SARS-CoV-2 would induce lymphopenia with decreased levels of lymphocyte subset cells, particularly Ts cells, and subsequently induce hyperinflammatory response and cytokine storm. Levels of IL-6, IL-8, IL-10, and Ts cell count might be the independent predictors of disease severity and adverse endpoints and additions to prognostic factors. Further researches on the immune pathogenesis and effect are of importance for the clinical management and development of therapeutics and vaccines of COVID-19.

Data Availability

The datasets analyzed in the study are available from our corresponding author by request.

Conflicts of Interest

The authors declare no conflicts of interest.

Authors’ Contributions

NJ, ZL, TL, and QY designed the article. ZL, BY, MJ, YS, YH, YL, and YW collected and analyzed the data. NJ, ZL, and BY drafted of the manuscript. DS, PM, JZ, TL, and QY discussed and revised the manuscript. Nan Jiang and Zhijun Li contributed equally to the work.

Acknowledgments

The authors are grateful to all health workers.

Supplementary Materials

Correlation coefficient and value between cytokine profiles and lymphocyte subsets in COVID-19. (Supplementary Materials)