Introduction

Since the outbreak in Wuhan, Hubei Province, China in December 2019, Coronavirus disease 2019 (COVID-19) caused by severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) has spread rapidly around the world, resulting in a very high number of infections and deaths, threatening public security, and placing a heavy burden on global healthcare systems and economy (Guan et al. 2020; Organization 2020b; Zheng 2020). Although it is less virulent than Middle East respiratory syndrome coronavirus (MERS-CoV) and severe acute respiratory syndrome coronavirus (SARS-CoV), a greater number of infections and deaths have been caused by SARS-CoV-2 because of its higher transmissibility (Cheng et al. 2007; Chan et al. 2015; Li Q et al. 2020). Unfortunately, to date, there are no specific antiviral drugs or licensed vaccines for COVID-19, and isolation is the most effective way to control the spread of infection (Skurnik et al. 2020). Therefore, serial data on viral RNA dynamics and antibody responses in COVID-19 patients are urgently needed to guide the management of the pandemic and the development of antiviral therapies and vaccines.

Although there has been research on the kinetics of viral shedding and antibody responses during SARS-CoV-2 infection, the data are limited by the short monitoring period (around 40 days) (Long et al. 2020a, b; Okba et al. 2020; To et al. 2020; Zhao et al. 2020). As such, long-term changes in antibody profiles in cases of SARS-CoV-2 infection and how long these antibody responses will maintained remain unknown. In order to gain a greater understanding of viral shedding and the immune response in humans following SARS-CoV-2 infection, here we described the long-term profile of viral nucleic acid shedding in respiratory tract and the longitudinal antibody response against the nucleocapsid protein (NP) and the surface spike protein receptor-binding domain (RBD) of SARS-CoV-2 in 112 patients admitted to the Honghu People’s Hospital, Hubei Province, China. Besides, 18 asymptomatic individuals together with 17 symptomatic COVID-19 patients were also included in our study, and the asymptomatic individuals were the close contacts of these 17 symptomatic cases.

Materials and Methods

Patients

A total of 112 patients who were diagnosed by real-time quantitative PCR (RT-qPCR) testing between January and March 2020 and admitted to Honghu People’s Hospital, Hubei Province were included in this retrospective cohort study. These patients were followed up biweekly until 80-days post onset. Throat swabs and serum samples were collected for viral load monitoring and measurement of antibody levels, respectively. Demographic information and clinical examination results were obtained from the electronic medical records. The diagnosis of COVID-19 and evaluation of disease severity is based on the 6th edition of the New Coronavirus Pneumonia Prevention and Control Program released by the National Health Commission of the People’s Republic of China (2020). As described, SARS-CoV-2 infection was diagnosed based primarily on clinical manifestations, chest computerized tomography, and detection of viral nucleic acid in respiratory tract swab samples by RT-q PCR, which is currently the gold standard diagnostic test. Patients presenting with fever, mild clinical symptoms, or mild pneumonia by chest imaging were defined as mild cases, while severe cases were those with respiration rate ≥ 30 times/min, blood oxygen saturation ≤ 93% in the resting state, arterial partial pressure of oxygen/fraction of inspiration oxygen ≤ 300 mmHg, or chest imaging showing > 50% lesion progression within 24–48 h. The date of symptoms onset was set as the day that fever, cough, chills, dizziness, headache, or fatigue appeared.

In late July, we conducted another round of follow-up in which 54 convalescent patients participated. The follow-up time ranged from 158 to 194 days post onset of symptoms. Serum samples were collected for antibodies measurement.

Additionally, 18 individuals who were confirmed to be infected with SARS-CoV-2 by RT-qPCR but did not show any signs or symptoms and had normal chest imaging were included in the analysis together with 17 symptomatic COVID-19 patients, and the asymptomatic individuals were the close contacts of these 17 symptomatic cases.

SARS-CoV-2 RNA Measurement

SARS-CoV-2 RNA load in throat swabs was measured by RT-qPCR using primers and probes targeting the NP and the open reading frame1b genes of SARS-CoV-2 according to the World Health Organization (WHO) guidelines (Organization 2020a). Specimens with cycle threshold values (Ct) < 40 were considered positive; Ct > 40 was assigned a value of 40. The diagnostic criteria were according to the recommendations of the National Institute for Viral Disease Control and Prevention in China (China. 2020).

Serology Testing

IgM and IgG antibodies against SARS-CoV-2 NP and RBD were detected using enzyme-linked immunosorbent assay (ELISA) kits (Guangzhou Darui Biotechnology Co., Guangzhou, China) according to the manufacturer’s instructions. Briefly, these assays employed the quantitative sandwich enzyme immunoassay technique. Serum samples were diluted at 1:100 and added to 96-well polystyrene microplates coated with recombinant SARS-CoV-2 protein antigens, followed by secondary antibody labelled with the conjugate. The optical density (OD) at 450 nm of each well was measured to quantify the amount of antibody in the serum sample. Positivity was determined based on the manufacturer’s recommended cutoff values of 0.4 for IgM and 0.3 for IgG. The total binding IgG and IgM antibodies against SARS-CoV-2 NP in serum were measured by a chemiluminescence microparticle immunoassay (CIMA). Antibody levels were expressed as the ratio of the chemiluminescence signal over the cut-off value (S/CO). A S/CO value of > 1.0 was considered positive.

Statistical Analysis

Prism v8.0 (GraphPad, La Jolla, CA, USA) and SPSS v26.0 (SPPS Inc, Chicago, IL, USA) were used for statistical analyses. Statistical comparisons between groups were performed using the Mann–Whitney U test for continuous variables, and the χ2 test for categorical variables. Correlations were assessed using a Spearman’s Rank correlation coefficient (r). All tests were 2-sided, and P values < 0.05 were considered statistically significant.

Results

Demographic and Clinical Characteristics of the Study Population

Of the 112 patients confirmed as positive for SARS-CoV-2 infection between January and March 2020, 87 (77.7%) were mild cases and 25 (22.3%) were severe. The demographic and clinical characteristics of each group are summarized in Table 1. The median age of the severe group was older than the mild group (median age, 65 vs 50 years, P < 0.001), and showed a tendency for higher male-to-female ratio (68.0% vs 51.7%, P = 0.176). Fever (76.8%) was the most common symptom among the 112 patients, followed by dry cough (54.5%), fatigue (33.0%), and shortness of breath (19.6%). Significantly higher proportions of patients in the severe disease group experienced fatigue (52.0% vs 27.6%, P = 0.030), chill (24.0% vs 5.8%, P = 0.015), and palpitation (20.0% vs 3.5%, P = 0.013), and showed greater abnormalities in biochemical indices at the time of admission than the milder patients group, including higher levels of C-reactive protein (51.57 vs 6.27 mg/L, P < 0.001), d-dimer (1.09 vs 0.7 mg/L, P = 0.033), and lactate dehydrogenase (306 vs 226.5 U/L, P = 0.003). Leukocyte and neutrophil counts were higher in the severe group than in the mild group (7.26 × 109 vs 5.48 × 109/L, P = 0.016 and 5.72 × 109 vs 3.33 × 109/L, P < 0.001, respectively), while the lymphocyte count was lower (0.93 × 109 vs 1.35 × 109/L, P = 0.003) in patients with more severe disease.

Table 1 Clinical characteristics of COVID-19 patients.

SARS-CoV-2 RNA Dynamics

We collected 1065 throat swabs from 112 COVID-19 patients for viral load monitoring (mean of 9.5 throat specimens per patient). The median of viral RNA following-up time was 57 days (interquartile range [IQR], 49–63.75 days), and the longest was 80 days. The temporal profile of the viral load for all patients is shown (Fig. 1A, 1B). The peak viral load in throat swabs was significantly higher in the severe group than in the mild group but did not differ by sex and showed no correlation with age (Fig. 1C–1E). The median duration of viral shedding was 15 days (IQR, 9–27 days). In five patients (4.5%; four mild and one severe case), viral RNA had turn negative for two consecutive tests and then recurrence occurred 45 days after the initial symptom onset. In another patient with mild disease, viral RNA was detected by RT-qPCR up to 45 days after symptom onset.

Fig. 1
figure 1

SARS-CoV-2 RNA profiles in throat swabs of COVID-19 patients. A Temporal profile of viral load (inversely related to cycle threshold [Ct] value) in all patients (n = 112). Each line represents an individual patient. The thick line shows the trend in viral load using smoothing splines. B SARS-CoV-2 RNA positive rate at different time intervals. C Comparison of peak viral load according to sex. D Relationship between age and peak viral load. E Comparison of peak viral load between mild and severe patients. Results are shown as median and interquartile range; P values were derived from Mann–Whitney U test, or Spearman's test; r values were derived from Spearman's test.

Antibody Response Against SARS-CoV-2 NP and RBD

Longitudinal antibody responses were evaluated in our cohort. A total of 474 serum samples were obtained from the 112 patients (mean of 4.23 serum specimens per patient). These patients were followed up biweekly up to 80-days post disease onset. IgM and IgG antibodies against SARS-CoV-2 NP were detected by ELISA. To further understand the correlation between antibody responses to SARS-CoV-2 and clinical parameters, the antibody responses were stratified according to sex, age, and disease severity (Fig. 2A). Overall, the levels of IgM antibodies varied in the two groups, while there was no significant difference in IgG levels when patients were stratified according to sex, disease severity, and age. Anti-NP IgM levels peaked at the first 2 weeks after symptom onset and decreased steadily thereafter, and were undetectable in most patients after 4–6 weeks. Compared to females, male patients had significantly higher serum anti-NP IgM level; and the anti-NP IgM level was higher in older patients and the severe patients group compared to the younger group and the mild disease group, respectively, but the differences were not statistically significant. Anti-RBD IgM level reached its peak at 3–4 weeks post disease onset, which occurred later than anti-NP IgM level.

Fig. 2
figure 2

Profiles of antibodies against SARS-CoV-2 NP and RBD in COVID-19 patients (n = 112). A Antibodies against SARS-CoV-2 NP and RBD in patients according to sex, age, and disease severity. Results are shown as median and interquartile range; P values were derived from the Mann–Whitney U test. B Seropositive rates of IgM and IgG against SARS-CoV-2 NP at different intervals. C, D Relationship between levels of antibodies against SARS-CoV-2 NP (C) and RBD (D) measured by ELISA. The P and r values were derived from Spearman’s test. *P < 0.05; **P < 0.01.

We then evaluated the seropositive rates of anti-NP IgM and IgG antibodies based on the number of days after symptom onset (Fig. 2B; Supplementary Table S1). The positivity rate for IgM peaked in the second week with 86.49% and declined steadily afterward. The seroconversion rate of anti-NP IgG was 69.2% in the first week after symptom onset, reached 100% in the fifth week, and then remained 100% up to 80 days post disease onset.

Next, we estimated Spearman’s correlation among the antibody titres against NP and RBD (Fig. 2C, 2D). ELISA OD values for anti-RBD IgM showed a weak but significant positive correlation with anti-NP IgM levels (r = 0.285; P < 0.001; Fig. 2C), and anti-RBD IgG levels showed a significant positive association with anti-NP IgG levels (r = 0.623; P < 0.001; Fig. 2D). While there was no correlation between the persistence of SARS-CoV-2 RNA and antibodies levels (Supplement Fig. S1).

Follow-Up Study on Convalescent Patients

In late July, we conducted another round of follow-up in which included 54 recovered patients. The follow-up time ranged from 158 to 194 post onset of symptoms. All the samples collected in July tested positive (S/CO > 1) for antibodies against NP; the temporal profile is shown in Fig. 3A. At 158–194 days post onset of symptoms, the antibody levels decreased approximately 46.17% overall (median, 151.61 S/CO, IQR, 287.11–453.01 S/CO; P = 0.128; Fig. 3B) compared with those at the acute phase (the period when the viral RNA can be found in a respiratory specimen). Moreover, the median percentage of decrease was 46.31% (IQR 5.37%–68.96%) for antibodies against NP in the mild group, whereas it was 46.02% (IQR 11.01%–65.75%) in the severe group.

Fig. 3
figure 3

Profile of antibodies against SARS-CoV-2 NP measured by CIMA. A Temporal pattern of antibodies to NP in all patients (n = 112). Each line represents an individual patient. The thick line shows the trend using smoothing splines. B Dynamic changes in antibodies against NP in the acute phase (the period when the viral RNA can be found in a respiratory specimen) and 158–194 days post symptoms onset (n = 54). The P and r values were derived from Spearman’s test.

Comparison Between Asymptomatic Individuals and the related Symptomatic Cases

We compared SARS-CoV-2 RNA dynamics and antibody responses between 18 asymptomatic individuals and 17 symptomatic COVID-19 patients, and the asymptomatic individuals were the close contacts of these symptomatic cases. The clinical characteristics of these individuals are summarized in Table 2. The symptomatic group tended to have a higher median age and included more males than the asymptomatic group.

Table 2 Characteristics of asymptomatic individuals and the related symptomatic cases.

In total, 276 throat swabs (127 from 18 asymptomatic individuals and 149 from 17 symptomatic cases) were collected for viral RNA load monitoring. There were no statistically significant differences in the peak Ct value and the persistence of viral RNA between the asymptomatic individuals and the symptomatic cases (Fig. 4A–4C). A further 116 serum samples (2–5 per patient) were obtained at different time points and tested for antibodies against NP. The date of diagnosis was set as the starting point to compare the antibody dynamics between the asymptomatic individuals and symptomatic cases. The temporal profile of antibodies is shown in Supplementary Figure S2. Antibody levels against NP were significantly lower in the asymptomatic group (Fig. 4D). The serologic conversion of IgG occurred earlier in the symptomatic cases than in the asymptomatic carriers, and the IgM seroconversion rate in the asymptomatic group was much lower (16.7% vs 70.6%, P = 0.002; Fig. 4E).

Fig. 4
figure 4

SARS-CoV-2 RNA and antibody dynamics in asymptomatic patients (n = 18) and symptomatic cases (n = 17). A Temporal profile of viral RNA load (inversely related to cycle threshold [Ct] value). Each line represents an individual patient. The thick line shows the trend in viral load using smoothing splines. B, C Peak viral RNA load (B) and persistence of viral RNA (C) in asymptomatic patients and symptomatic cases. D Antibodies against NP measured by CIMA. E Seropositivity of antibodies against NP at different intervals measured by ELISA. Results are shown as median and interquartile range; P values were derived from the Mann–Whitney U test.

Discussion

This study illustrates the kinetics of SARS-CoV-2 RNA levels and the antibody responses to SARS-CoV-2 NP and RBD proteins in COVID-19 patients with different illness severities, including mild, severe, and asymptomatic. Compared to the peak viral load of SARS-CoV occurring at approximately 10 days and that of MERS-CoV occurring in the second week after symptom onset (Peiris et al. 2003a; Oh et al. 2015; Oh et al. 2016), in the present study, peak viral loads in throat swabs were in most cases already detectable at the initial presentation of symptoms and decreased during hospitalization. This finding is consistent with previous reports (Huang et al. 2020; To et al. 2020) and supports the view that peak transmissibility occurs at or before symptom onset (He et al. 2020; Rothe et al. 2020). Thus, strategies used to control the spread of SARS-CoV may not be applicable to SARS-CoV-2, which require stricter measures based on the viral dynamics. An association between higher viral load and worse outcome has been reported for SARS-CoV and MERS-CoV (Peiris et al. 2003b; Hung et al. 2004; Oh et al. 2016). Similarly, we observed much higher peak viral RNA loads and stronger antibody responses in patients with severe COVID-19 symptoms compared to patients with mild disease. The correlation of antibody responses with disease severity is still unclear. Higher viral load may lead to more severe disease severity and may stimulate the production of antibodies. In turn, a robust antibody response may result in tissue injury and drive the symptoms of COVID-19.

A few studies have reported prolonged shedding of SARS-CoV-2 RNA (Li N et al. 2020; Liu et al. 2020; Xiao et al. 2020; Yang et al. 2020). Notably, in our study, 5.36% (6/112) of patients showed prolonged viral shedding (> 45 days in throat swabs) and five patients among them had two successive swabs that were negative but experienced recurrence of COVID-19, although this may have been due to false-negative RT-qPCR results rather than actual recurrence or reinfection (Xiao et al. 2020). Another patient tested positive for SARS-CoV-2 RNA > 45 days post-symptom onset. However, antibody responses in these six patients did not differ from those in other patients. These data indicate that self-quarantine and regular monitoring of viral RNA levels are necessary in patients who have recovered from COVID-19, even after they have had 2 consecutive negative nucleic acid tests and met the discharge criteria.

All patients in this study achieved antibody seroconversion within a month post symptom onset. IgM is usually produced before IgG during viral infections (Aljumaili et al. 2014; Singh et al. 2015; Samanta and Willis 2016). This was also true in the case of SARS-CoV-2 infection: the median time for seroconversion was 17 days post symptom onset for anti-NP IgM and 20 days after the onset for anti-NP IgG. However, this finding contradicts earlier observations that the IgG response occurred simultaneously with or slightly earlier than the IgM response (To et al. 2020; Zhao et al. 2020). Some factors might have contributed to the variation in antibody response in different studies, including the frequency of specimen collection and different test kits with diverse specificity and sensitivity. Previous studies investigating the kinetics of antibodies to SARS-CoV-2 are limited to around 40 days after the onset of symptoms (Long et al. 2020a, 2020b; Okba et al. 2020; To et al. 2020; Zhao et al. 2020). The antibody response in patients following MERS-CoV infection could last more than 34 months after onset (Payne et al. 2016), and antibodies could still be detected up to 36 months post infection in 89% of individuals with SARS-CoV infection (Cao et al. 2007).

In this study, we described for the first time the 6-month longitudinal dynamics of the antibody response to SARS-CoV-2 in a cohort of COVID-19 patients. Our data showed that antibodies against SARS-CoV-2 NP can persist more than 194 days post illness onset, but most patients showed a variable degree of reduction in antibody levels. Although the antibody level was higher in patients with more severe symptoms in the acute phase, the degree of reduction was around 46% in the late convalescent stage (158–194 days after the onset) in both mild and severe patients. Whether the magnitude of antibody levels in the convalescent sera could prevent re-infection of SARS-CoV-2 is unknown and needs further research. Moreover, due to the decreasing trend observed in present study, antibodies to SARS-CoV-2 may disappear over a period of 18 months, suggesting the necessity of vaccine development.

The number of cases of asymptomatic SARS-CoV-2 infection is increasing rapidly (Chang et al. 2020; Li R et al. 2020; Team 2020). Asymptomatic individuals may be unaware of their infection status and therefore unknowingly transmit the virus to others if they are not isolated for medical observation. We monitored changes in viral RNA loads in throat swabs and serum antibody levels in several asymptomatic cases and symptomatic cases. In the former group, SARS-CoV-2 RNA persisted for a shorter time compared to the latter group and the anti-NP antibody levels were lower. Delayed IgG seroconversion and a lower IgM seropositive rates were observed in asymptomatic individuals, suggesting that asymptomatic individuals had a weaker immune response to SARS-CoV-2 infection and the antibody levels had a great potential to predict the progression of COVID-19.

Our study had some limitations. Firstly, we reported the viral RNA cycle threshold value, which does not distinguish between active and dead virus. Secondly, the time points of the sample collection were irregular and mainly relied on physician’s judgement, with some specimens collected after long intervals. Finally, given the small sample size, our results must be interpreted with caution and validated in a larger cohort.

In summary, we have determined the dynamics of the antibody response to SARS-CoV-2 up to 6 months after the onset of COVID-19 at a population level. IgM was detected earlier than IgG, and disappeared in most cases at around 4–6 weeks after symptoms onset, while IgG could last for more than 194 days after the onset, but patients showed approximately 46% reduction in antibody levels compared with the acute phase. These data imply that the antibody may persist for quite a long time and the vaccine could provide protection against SARS-CoV-2 infection. Delayed IgG seroconversion and a lower IgM seropositive rate were observed in asymptomatic individuals, indicating that asymptomatic individuals had a weak immune response to SARS-COV-2 infection.