Introduction

A novel β-coronavirus, formally named as the severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) by the International Committee on Taxonomy of Viruses, was identified as the cause of an outbreak that occurred in Wuhan in early December 2019. The epidemic of SARS-CoV-2 has now spread to over 210 countries, territories or areas with more than 7,000,000 confirmed cases worldwide [1]. The disease caused by SARS-CoV-2 has been recently named as coronavirus disease-2019 (COVID-19) by the World Health Organization. Patients with COVID-19 may develop mild or severe symptoms of acute respiratory infection. The mild and moderate cases show symptoms of fever, dry cough, fatigue, with or without abnormal chest CT findings, but with a good prognosis [2,3,4]. In contrast, some patients develop severe pneumonia, acute respiratory distress syndrome (ARDS) or multiple organ failure, with death rates ranging from 4.3 to 15% according to different study reports [2, 5,6,7]. Factors leading to severe COVID-19 have not been determined yet. Several viral and immunological parameters may provide clues for the underlying mechanisms. Among them, the viral loads may reflect the magnitude of SARS-CoV-2 infection in patients. Lymphopenia and inflammatory cytokine storm have been found to be associated with the disease severities in highly pathogenic severe acute respiratory syndrome coronavirus (SARS-CoV) and Middle East respiratory syndrome coronavirus (MERS-CoV) infections [8, 9]. In addition, decreased peripheral blood lymphocyte counts and elevated serum inflammatory cytokine levels have been reported in COVID-19 patients [5, 10]. However, a thorough investigation into the dynamic changes in different lymphocyte subsets, the kinetics of inflammatory cytokines changes in the peripheral blood of COVID-19 patients, and the relationships between these parameters and disease severity is still lacking.

In this study, a cohort of 154 clinically confirmed patients with positive SARS-CoV-2 detection by RT-PCR were enrolled and clinically treated and monitored for up to 4 weeks after admission. The clinical and laboratory results were collected and analyzed. Particularly, the viral loads in patients were assessed by using the Ct values of specific RT-PCR tests. The changes in lymphocyte subsets and cytokines profiles in the peripheral blood of COVID-19 patients were analyzed. The results of our analysis demonstrated that the initial SARS-CoV-2 RNA loads varied in patients, but were comparable in different patient groups stratified by age, gender, comorbidities and disease severity. In contrast, reduced counts of T cells and subsets were found in patients with old age and comorbidities at admission and were characteristic of severe COVID-19. Thus, severe COVID-19 developed preferentially in patients with compromised immunity, but was not associated with initial viral loads.

Materials and methods

Data collection

Written informed consents were obtained from all 154 patients upon admission into the Second People's Hospital, Fuyang, China. This retrospective single-center study was approved by the Ethics Committee of the Second People's Hospital, Fuyang (Number 2020020004). Medical record information including epidemiological, demographic, clinical manifestation, laboratory data, and outcome data was obtained for all patients.

The diagnosis of COVID-19 for patients was performed according to the Guidelines for the Diagnosis and Treatment of New Coronavirus Pneumonia (version 7) published by the National Health Commission of China. Mild and moderate cases met the following criteria: (1) epidemiological history, (2) fever or other respiratory symptoms, (3) frequently with typical CT image abnormities of viral pneumonia, and (4) positive RT-PCR result for SARS-CoV-2 RNA. Severe patients additionally met at least one of the following conditions: (1) shortness of breath, RR ≥ 30 times/min, (2) oxygen saturation (resting state) ≤ 93%, or (3) PaO2/FiO2 ≤ 300 mmHg.

Laboratory examination

Laboratory confirmation of SARS-CoV-2 infection was performed with real-time RT-PCR (Bojie Ltd, Shanghai, China) per Chinese CDC recommendation. The RT-PCR tests use two primer pairs and corresponding probes targeting the ORF1ab and N regions of the SARS-CoV-2 genome [11]. Sputum and/or throat-swab specimens were collected from all patients and the samples were maintained in viral-transport medium prior to laboratory testing. A test result was considered positive if the Ct values of both RT-PCR reactions were below 38 according to the guidelines published by the National Health Commission of China. Influenza A and B infections were excluded by real-time RT-PCR. Specific IgG and IgM for coxsackievirus, respiratory syncytial virus, parainfluenza virus, enterovirus, Chlamydia and Mycoplasma pneumonia were tested by chemiluminescence immunoassay. Specimens, including sputum or alveolar lavatory fluid, blood, urine, and feces, were cultured to identify pathogenic bacteria or fungi that may be associated with the SARS-CoV-2 infection. A lymphocyte test kit (Mindray Bio-medical Inc., Shenzhen, China) was used for lymphocyte subset analysis. Plasma cytokines (IL-6) were measured with Human Th1/2 Cytokine Kit II (Jinzhun Ltd., Shenzhen, China). All tests were performed according to product manuals.

Statistical analyses

Categorical variables were expressed in frequency or percentage, and significance was detected by Chi square or Fisher’s exact test. Continuous variables were expressed as mean ± standard deviation or mean ± standard error, and the significance was tested by t test. Nonparametric variables were expressed in median and quartile intervals or mean ± standard error, and significance was tested by Mann–Whitney U, ANOVA or Kruskal–Wallis test. P < 0.05 was considered statistically significant in all statistical analyses. Principal component analysis (PCA) was performed to identify the major contributing factors among clinical parameters to distinguish mild and severe cases of COVID-19 patients. The diagnostic values of selected parameters for differentiating mild and severe cases of COVID-19 patients were assessed by receiver operating characteristic (ROC) and the area under the ROC curve (AUC). SPSS statistical software (Macintosh version 26.0, IBM, Armonk, NY, USA) and R package were used for statistical analysis.

Results

Demographic and clinical characteristics of COVID-19 patients

A total of 154 patients were enrolled in this study and were divided into two groups according to the published guidelines, including 122 mild/moderate cases (79.2%, referred to as mild cases thereafter) and 32 severe/critical cases (20.7%, referred to as severe cases). Among 122 patients with mild symptoms, eight patients had no chest CT imaging abnormality. Among the 32 patients with severe symptoms, 4 patients required ICU care and ventilator-aided breathing.

The enrolled COVID-19 patients consisted of 86 males (55.8%) and 68 females (44.1%) (Table 1). The mean age of the patients was 42.4 years. The medium age of the severe symptom group (52.2 ± 14.0 years) was older than that of the mild group (39.9 ± 15.0 years). A total of 52 (33.8%) patients in both groups had baseline diseases, including hypertension (23 [14.9%]), diabetes (11 [7.1%]), chronic liver diseases (12 [7.8%]), coronary heart diseases (3 [1.9%]), malignant tumors (3 [1.3%]), cerebrovascular diseases (2 [1.3%], COPD (2 [1.3%]) and AIDS (1 [0.65%] (Table 1). Two severe group patients had mixed fungal and six bacterial infections, respectively. At admission, 31 (96.9%) patients of the severe group and 86.4% of mild group patients had fever, with a significant difference in body temperature (Table 1). The levels of the inflammatory markers ESR (P = 0.008), C-reactive protein (CRP, P < 0.001), serum amyloid A protein (SAA, P = 0.039), interleukin-6 (IL-6, P < 0.001), and other markers including albumin (P < 0.001), aspartate transaminase (P = 0.003), alanine transaminase (P = 0.043), lactate dehydrogenase (P < 0.001), fibrinogen (P = 0.004) and d-dimer (P < 0.001) in the peripheral blood of the severe group patients were significantly higher at admission compared to the mild group patients (Table 2). No significant differences in the numbers of red blood cells and platelets, or in the serum levels of procalcitonin, total bilirubin, creatinine, creatine kinase and PT were observed between the two groups (Table 2).

Table 1 Baseline characteristics of patients infected with SARS-CoV-2
Table 2 Laboratory findings of patients infected with SARS-CoV-2 on admission to hospital

The association of the severity of COVID-19 with comorbidities in our patient cohort was also analyzed (Fig. 1). Overall, patients with underlying diseases were more prone to developing severe COVID-19. However, patients with comorbidities such as chronic hepatitis B virus infection, liver cirrhosis, chronic respiratory diseases, HIV infection, and cancers (comorbidity group 1, CG-1) which are known to impair immunological functions were significantly more frequently affected compared to those with comorbidities unrelated to immune functions (comorbidity group 2, CG-2) (Fig. 1a), in the age-matched groups (Fig. 1b).

Fig. 1
figure 1

The association of the severity of COVID-19 with underlying diseases in the patient cohort. a The association of the severity of COVID-19 with comorbidities. The comorbidity group 1 (CG 1) includes chronic hepatitis B virus infection, liver cirrhosis, chronic respiratory diseases, HIV infection, and cancers which are known to impair immunological functions. The comorbidity group 2 (CG2) includes psychological diseases and fatty liver diseases which are primarily not related to host immune functions. The risk to develop severe COVID-19 was shown for patients without known comorbidities (No CG), CG1, and CG2. b The ranges of the ages of patients with and without comorbidities. The means ± SD are given in the graph. Results were tested for significance by Chi square test (a) or t test (b). A P value < 0.05 is considered as significant

Viral levels and the correlation with disease severities

The Ct values of RT-PCR reactions performed were used to assess viral RNA levels in patient samples. An absolute quantification of viral RNA loads could not be obtained due to a lack of proved standard material. The amplification reaction for the N region regularly resulted in lower Ct values compared with that for the ORF1ab region, indicating a better sensitivity with N region amplification. Therefore, the Ct values of N-specific RT-PCR tests were analyzed to assess the dynamic changes in viral RNA levels (Fig. 2). For the detection of SARS-CoV-2 RNA, sputum or throat swab specimens were used. Overall, sputum specimens gave higher viral levels indicated by lower Ct values (Ct values up to 17) than throat swab specimens (Ct value up to 21), which was consistent with an early report [12]. Thus, the viral load data generated by using sputum or throat swab specimens were analyzed separately. We compared the measured SARS-CoV-2 RNA levels in sputum specimens from COVID-19 patients at admission among groups divided according to age, sex, underlying diseases and disease severity (Fig. 2a). The Ct values varied from > 38 up to 17 and no difference between patient groups was found. Consistently, the results obtained with throat swab samples showed no difference between patient groups (Fig. 2b). However, viral RNA levels in severe cases were found to be significantly higher in sputum specimens, but not in throat swab specimens between days 3 to 5 after disease onset compared to those in mild cases (Fig. 2c). Nevertheless, a number of mild cases had relatively high viral loads with a Ct value of up to 20. The viral RNA levels in the patients decreased and became undetectable (Ct > 38) over time (Fig. 2d). The mild cases became negative for viral RNAs with a median of 21 days after the disease onset, which was significantly earlier than severe cases with a mean of 35 days based on the test results using sputum specimens (Fig. 2e). No significant difference was found between the mild and severe cases at 28 and 25 days after disease onset if throat swab specimens were used (Fig. 2e).

Fig. 2
figure 2

Detection of SARS-CoV-2 RNA in patient samples by using a specific RT-PCR test. The Ct values of RT-PCR test for the SARS-CoV-2N region were used to assess viral RNA levels in patient samples. Sputum or throat swab specimens were used. a, b The measured SARS-CoV-2 RNAs levels in sputum (a) and throat swab (b) specimens from COVID-19 patients at admission were compared according to the age, sex, comorbidity, and the disease severity. c SARS-CoV-2 RNA levels in sputum and throat swab specimens from the mild and severe cases between the days 3 and 5 after disease onset. d Kinetics of SARS-CoV-2 RNA levels (Ct values) measured by RT PCR with sputum and throat-swap (swap) specimens. ***P < 0.001. e The Meier–Kaplan survival analysis for the detection of SARS-CoV RNAs based on the test results using sputum specimens and throat swab specimens. The median and quartile intervals (ac) or means ± SD (d) are given in the graph. Results were tested for significance by t test, Mann–Whitney U or Kruskal–Wallis test (ad), and log-rank (Mantel–Cox) test (e). A P value < 0.05 is considered as significant

Kinetic analysis of lymphocyte subsets in the peripheral blood of COVID-19 patients

Lymphopenia was observed in 47.4% (73/154) of patients in our cohort, 40.2% (49/122) and 75% (24/32) in mild and severe cases at the onset of disease, respectively (Table 2). As shown in Table 2, the absolute counts of lymphocytes in the peripheral blood of patients subsequently developing severe COVID-19 was significantly lower at hospital admission and disease onset, respectively (Table 2, Fig. 3). Sustained lower lymphocyte counts in severe cases were observed, compared to those in mild cases (Fig. 3c). The absolute counts of total white blood cells (WBC) and neutrophils (NE) were initially comparable in both patient groups, but significantly increased only in severe cases at day 9 up to day 17 after disease onset (Fig. 3a, b). A significant difference in monocyte counts was observed between the two groups only at days 3–5, but not seen afterward (Table 2, Fig. 3d). The kinetic changes in different lymphocyte subsets in the peripheral blood of COVID-19 patients were monitored from disease onset up to 18 days. Flow cytometric analysis for CD3+ T cells, CD4+ and CD8+ T cell subsets, B cells and NK cells was performed (Fig. 4). Firstly, we analyzed initial lymphocyte counts and CD3+, CD4+ and CD8+ T cell counts at admission in patients by the age and presence of comorbidities (Fig. 4a, b). Clearly, patients of older age and with comorbidities of group 1 (CG1) had significantly lower counts for lymphocytes and T cell subsets. Consistent with the findings about lymphocytes, sustained significantly lower CD3+, CD8+ and CD4+ T cell counts were found in patients with severe disease compared to those of the mild cases during clinical observation (Fig. 4c). The lowest CD3+, CD4+ and CD8+ T cell counts in patients with severe COVID-19 were generally found 3–5 days after disease onset. However, they recovered and reached comparable levels to those in the mild cases after day 18 (Fig. 4d, Figure Suppl. 1). As severe COVID-19 cases had high SARS-CoV-2 RNA levels and low T, CD4+, and CD8+ cell counts, these parameters were tested for possible correlations. However, our analyses did not show any significant correlation between SARS-CoV-2 RNA level and T, CD4+or CD8+ T cell counts (Figure Suppl 2).

Fig. 3
figure 3

Kinetic analysis of cell counts of different populations of WBCs in COVID-19 patients. The absolute numbers of white blood cells (a), neutrophils (b), lymphocytes (c) and monocytes (d) in the peripheral blood of mild (blue line) and severe (red line) COVID-19 patients were analyzed at different time points after the disease onset. *P < 0.05; **P < 0.01; ***P < 0.001. The upper (purple line) and lower (green line) limits of normal values and the means ± SD are given in the graph. Results were tested for significance by t test or Mann–Whitney U test. A P value < 0.05 is considered as significant

Fig. 4
figure 4

Kinetic analysis of cell counts of different lymphocyte subsets in COVID-19 patients. The kinetic changes of different lymphocyte subsets in the peripheral blood of COVID-19 patients were measured from the disease onset up to 18 days later. Flow cytometric analysis to for CD3+ T cells, CD4+ and CD8+ T cell subsets, B cells and NK cells was performed. The initial lymphocyte counts and CD3+, CD4+ and CD8+ T cell counts at admission in patients according to age (a) and the presence of comorbidities (b). c The CD3+, CD8+ and CD4+ T cell counts in patients with mild and severe disease at 3–5 days after the disease onset. d The absolute numbers of CD3+ T cells, CD8+ T cells, CD4+ T cells, B cells and NK cells in the peripheral blood of mild (blue line) and severe (red line) COVID-19 patients were analyzed at different time points after the disease onset. *P < 0.05; **P < 0.01; ***P < 0.001. The median and quartile intervals (ac) or means ± SD (d) are given in the graph. Results were tested for significance by two-way ANOVA (a), t test (b), Mann–Whitney U test (c) and t test or Mann–Whitney U test (d). A P value < 0.05 is considered as significant

No significant differences in B cell and NK cell counts were observed between the two groups during the whole disease course (Fig. 4d).

Kinetic analysis of serum inflammatory markers in COVID-19 patients

A previous study demonstrated a marked elevation of inflammatory cytokine levels such as IL-6 in the sera of COVID-19 patients [2]. Therefore, we further characterized the kinetic changes of inflammatory markers, including CRP, IL-6, SAA, and PCT in the sera of the patients. The serum levels of these markers in patients with mild disease showed minor fluctuations, while significant elevations of these markers were found severe cases (Fig. 5). CRP reached its peak in serum at 6–8 days after disease onset (Fig. 5a). IL-6 levels showed sustained increases in the severe group compared to the mild group over 18 days, even when some other markers returned to normal levels (Fig. 5b). The IL-6 levels in some patients with severe COVID-19 remained high as late as day 12–14 after disease onset, even though these patients had clinically recovered (Fig. 5e).

Fig. 5
figure 5

Kinetic analysis of serum levels of inflammatory markers in COVID-19 patients. The concentrations of C-reactive protein (a), IL-6 (b), serum amyloid A (c), procalcitonin (d), serum IL-6 concentrations at days 3–5 and 12–14 (e) in mild (blue line) and severe (red line) COVID-19 patients were analyzed at different time points after the disease onset. *P < 0.05; **P < 0.01; ***P < 0.001. The means ± SD (ad) or median and quartile intervals (e) are given in the graph. Results were tested for significance by the Mann–Whitney U test. A P value < 0.05 is considered as significant

Significant increases in serum SAA and PCT levels were observed in the severe group after disease onset, but became only significantly different to those in the mild group at certain time points (Fig. 5c, d). The three examined markers CRP, SAA and PCT reached similar levels between the severe and mild patient groups at day 18 after disease onset (Fig. 5). A statistical analysis showed that SARS-CoV-2 RNA levels in patients did not correlate with the levels of serum inflammatory markers. No correlation was found between the initial levels of lymphocyte subsets in peripheral blood and serum inflammatory markers in the late period between days 12 and 14 after disease onset (not shown).

Factors associated with severe COVID-19

Next, we examined the possibility of using the above-mentioned parameters as prognostic factors for identifying severe cases in COVID-19 patients. PCA was firstly performed with R package “factoextra” to identify correlating variables for distinguishing severe patients from mild patients (Figure Suppl 3A). The lymphocyte and CD3+ T cell counts were found to be the two most important contributing variables among immunological markers. Previously, other factors such as neutrophil-to-CD8 + T cell ratio (N8R), neutrophil-to-lymphocyte ratio (NLR), neutrophil counts (NEC) and white blood cells counts (WBCC) had been proposed as potential prognostic factors for severe COVID-19 [13]. To assess the diagnostic value of these parameters, the receiver operating characteristic (ROC) curve and the area under ROC curve (AUC) were calculated by R package “pROC” (Figure Suppl 3B), and the AUC values for lymphocyte and CD3+ T cell counts were found to be 0.739 and 0.793, respectively. At the same time, the cutoff values were calculated from the ROC curves, with a value of 0.85 for lymphocyte counts (L) (specificity: 60.6%, sensitivity: 77.9%) and 665 for T cell counts (T) (665, 83.3%, 67.9%) (Figure Suppl 3B). Further, the ROC curves were calculated for other parameters: CD4+ T cell counts (CD4), CD8+ cell counts (CD8), and combined CD4+ and CD8+ T cell counts (CD4xCD8) with the AUC values 0.767 (cutoff: 269, specificity: 60.0%, sensitivity: 87.2%), 0.772 (236, 73.3%, 73.4%), and 0.800 (115,520, 90.0%, 62.4%), respectively (Figure Suppl 3B). Overall, the immune related parameters were found to have associations with the chance to develop severe COVID-19, but their diagnostic values need to be evaluated in larger cohorts.

Discussion

In this study, we analyzed the clinical features including SARS-CoV-2 RNA load and immunological characteristics of peripheral blood in a patient cohort with COVID-19 from Anhui Province, China. The major clinical features of this patient cohort are comparable with published data. The patients with severe COVID-19 were characteristic of older age and underlying diseases.

Our data revealed that the initial SARS-CoV-2 RNA levels in patients at admission did not differ with age, sex, comorbidities and disease severity at later time points. The viral loads in some mild cases were also high and comparable with those in severe cases. Thus, high viral load may be not the primary cause of severe COVID-19.

Recent studies reported that the lymphocyte counts were normal in COVID-19 patients with mild disease, and in contrast, the majority of patients with severe disease had lymphopenia [5, 13, 14]. Our study suggested that a number of patients, particularly those of old age and with underlying diseases might have initially low lymphocyte and T cell counts (Fig. 4). There was no correlation between initial SARS-CoV-2 RNA load and lymphocyte and T cell counts in patients at admission (Figure Suppl 2). A plausible explanation for these observations is that SARS-CoV-2 may replicate more efficiently in patients with compromised immunity and thus causes severe disease in these patients. Consistently, the viral loads in both mild and severe cases increased rapidly and reached peak levels within 3–5 days after disease onset. However, in severe cases, the viral loads were significantly higher than those in mild cases, especially at that time point. Strikingly, all patients requiring ICU treatment had high SARS-CoV-2 loads at days 3–5 after the onset of disease (Figure Suppl 4). Thus, accurate quantitative detection of SARS-CoV-2 loads may be an important parameter to guide the clinical handling of COVID patients.

Previous research on SARS-CoV and MERS-CoV infections has demonstrated the correlation between marked lymphopenia and the severity of the disease [15,16,17,18,19]. Altered antigen presenting cell (APC) function and impaired dendritic cell migration during SARS-CoV infection may result in insufficient T cell priming and thus contribute to decreased numbers of virus-specific T cells in the lungs [20, 21]. Moreover, delayed type I interferon response or inflammatory monocyte–macrophages derived pro-inflammatory cytokines could also sensitize T cells to undergo apoptosis [22]. In contrast, MERS-CoV was found to be able to infect large numbers of human immune cells, including dendritic cells, macrophages, and T cells [23,24,25]. MERS-CoV infection of T cells results in apoptosis mediated by a combination of extrinsic and intrinsic apoptosis pathways, which is believed to contribute to viral spread and severe immunopathology [25]. It remains unclear whether SARS-CoV-2 may induce T cell contraction by specific mechanisms. Yet, a correlation of initial SARS-Cov-2 RNA loads in patients with degrees of lymphopenia was not seen in our study.

Severe COVID-19 patients have significant and sustained higher levels of inflammatory markers including serum IL-6 concentrations than in mild cases, again suggesting that the magnitude of the cytokine storm is associated with disease severity. T cells are important for dampening overactive innate immune responses during viral infection [26, 27]. Thus, low T cell counts in patients with SARS-CoV-2 infection may result in aggravated inflammatory responses, while normal T cell levels and functions may maintain the immunological regulation and prevent overshooting inflammation. Early identification of risk factors for severe COVID-19 patients may facilitate appropriate supportive care and prompt access to intensive care unit if necessary. Some studies, such as one in a 61-patient cohort [13], reported that NLR was the most useful prognostic factor affecting the prognosis for severe COVID-19. However, the diagnostic value of such parameters to predict disease severity was not proven in our cohort involving a larger number of patients.

Finally, the differences in T cells, CRP, and IL6 in the patient groups classified by underlying disease were analyzed. The results showed that the differences between the No CG and CG1 groups were only compared in this study und significant at some points, as shown in Figure Suppl 5. A comparison of the No CG with CG2 did not show any significant difference. This is rather inconclusive due to the small sample size of CG2.

In conclusion, our analyses demonstrated that initial SARS-CoV-2 RNA load varied in patients, but were comparable in different patient groups classified according to age, gender, comorbidities and disease severity. In contrast, reduced counts of T cells and subsets were measured in patients with older age and comorbidities at admission and were characteristic of severe COVID-19. Thus, severe COVID-19 developed preferentially in patients with compromised immunity and was not associated with initial viral load. Our results highlighted the importance of close surveillance of SARS-CoV-2 infection for specific risk groups including elderly people and those with underlying diseases with impaired immunological functions.