Modeling the early temporal dynamics of viral load in respiratory tract specimens of COVID-19 patients in Incheon, the Republic of Korea.

OBJECTIVE: To investigate the duration and peak of severe acute respiratory syndrome coronavirus 2 shedding as infectivity markers for determining the isolation period.
METHODS: A total of 2,558 upper respiratory tract (URT) and lower respiratory tract (LRT) specimens from 138 patients with laboratory-confirmed coronavirus disease were analyzed. Measurements of sequential viral loads were aggregated using the cubic spline smoothing function of a generalized additive model. The time to negative conversion was compared between symptom groups using survival analysis.
RESULTS: In URT samples, viral RNA levels peaked on day 4 after symptom onset and rapidly decreased until day 10 for both E and RdRp genes, whereas those in LRT samples immediately peaked from symptom onset and decreased until days 15.6 and 15.0 for E and RdRp genes, respectively. Median (interquartile range) time to negative conversion was significantly longer in symptomatic (18.0 [13.0-25.0] days) patients than in asymptomatic (13.0 [9.5-17.5] days) patients. The more types of symptoms a patient had, the longer the time to negative conversion.
CONCLUSIONS: The viral load rapidly changes depending on the time after symptom onset; the viral shedding period may be longer with more clinical symptoms. Different isolation policies should be applied depending on disease severity.

Introduction

From the beginning of the coronavirus disease (COVID-19) outbreak in January 2020 through 31 October 2020, a total of 22,756 patients were infected and 464 deaths were recorded in Korea. Sporadic community outbreaks with unknown sources of infection persisted. The proportion of infected young people aged 20–39 years reached 32.9%, and the vast majority (90.9%) of cases were mild (Korea Centers for Disease Control & Prevention, 2020); this is consistent with the findings from an initial study in China, which reported that 81% of COVID-19 cases were mild (Wu and McGoogan, 2020).

To minimize additional person-to-person transmission of severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2), it is important to identify and immediately isolate people with a confirmed infection until they are no longer infectious. Currently, the diagnostic testing and isolation of COVID-19 patients depend on the result of a real-time polymerase chain reaction (RT-PCR) test for SARS-CoV-2, the result (positive or negative) of which is interpreted using the cycle threshold (Ct) value. Although a positive RT-PCR result may not necessarily mean that a person is infectious, a lower Ct value (inversely related to the viral load) could indicate a higher viral load and, consequently, higher infectivity. Understanding the time series data of Ct values from early infection to disease recovery could guide the prediction of disease progression and clinical decisions.

To date, numerous studies have been published on the viral kinetics of SARS-CoV-2, its clinical features, and the prognosis of COVID-19 patients. Many of the studies have been published in China, mostly focusing on severe to critical cases (Han et al., 2020, He et al., 2020, Qi et al., 2020, To et al., 2020, Xu et al., 2020, Yu et al., 2020, Zheng et al., 2020, Zhu et al., 2020). Although many studies have addressed the temporal aspects of viral kinetics in COVID-19 patients, evidence pertaining to the viral load during the first few days is limited owing to the lack of data on the results of serial RT-PCR testing.

Since the early phase of the epidemic, Korea has performed active surveillance through immediate contact tracing, widespread diagnostic testing, and isolation of confirmed COVID-19 patients in a designated facility. This has enabled the detection of asymptomatic, presymptomatic, and mildly symptomatic patients in the early phase of the disease course. Thus, it has been possible to assess RT-PCR results as a proxy measure of disease severity and clinical progression (COVID-19 National Emergency Response Center et al., 2020).

Incheon Medical Center (IMC) is a public general hospital in Korea located in the Incheon metropolitan area, which has a population of approximately 3 million. Because of its close proximity to Incheon International Airport, IMC has been on the front line of the management of various imported infectious diseases; the first Korean patient with COVID-19 was treated at IMC. As of 31 October 2020, more than 900 COVID-19 cases have been reported in the city, the second largest number among the six metropolitan cities in Korea.

This study aimed to describe the temporal patterns of viral load from symptom onset to negative conversion and to assess the factors associated with the natural history of COVID-19 and delayed negative conversion among 138 mild or asymptomatic patients.

Methods

Study design and participants

This retrospective cohort study included 138 patients with laboratory-confirmed COVID-19 who were admitted to IMC between 19 January and 24 June 2020. IMC is classified as a group B medical institution that treats patients aged 15–65 years with underlying diseases or patients aged > 65 years without underlying diseases. Therefore, most of the study participants were middle-aged patients with mild to moderate COVID-19 symptoms. The patients’ demographic information and clinical data, including data on initially observed signs and symptoms (e.g., fever, cough, sore throat, sputum production, and headache), comorbidities, and date of first symptom onset were collected at hospital admission.

Daily clinical characteristics, such as body temperature, heart rate, respiration rate, blood pressure, type of treatment (e.g., admission to the intensive care unit, noninvasive ventilation, high-flow nasal cannula oxygen therapy, and antiviral treatment), and clinical outcome assessed using the clinical severity score were recorded on a standardized clinical record form during hospitalization. Clinical severity was classified on an eight-category ordinal scale developed by Sung et al. (2020), which is a modified version of the scale developed by Marshall et al. (2020). In this study, clinical severity scores were grouped into three categories: “mild” defined as scores 1–2 (without limitation of daily activities or with limitation of daily activities but without a need for supplemental oxygen therapy); “moderate” defined as scores 3–4 (with a need for supplemental oxygen therapy); and “severe” defined as scores 5–8 (with a need for noninvasive mechanical ventilation or death). By using the daily clinical data recorded for each patient, the median (interquartile range [IQR]) duration of clinical symptoms (i.e., from the initial symptom onset to the last symptom manifestation), real-time RT-PCR positivity (i.e., from the first positive RT-PCR result to the first and second negative results of two consecutive RT-PCR tests), and hospitalization (i.e., from diagnosis or hospital admission to discharge) were calculated.

Cessation of patient isolation was decided on the basis of two consecutive negative results obtained 24 h apart, at 1–2 weeks after the initial diagnosis, following the guidelines of the Korean Center for Disease Control and Prevention (KCDC).

Specimen collection and RT-PCR assays

Patients were confirmed to have the infection if they had a positive RT-PCR result from both upper respiratory tract (URT) and lower respiratory tract (LRT) specimens collected following the guidelines of the KCDC. URT specimens were collected from both nasopharyngeal and oropharyngeal swabs obtained by trained medical staff. Although sputum collection has not been mandatory for mild patients since the publication of the seventh edition of the KCDC guidelines on 02 March 2020, even mild patients admitted to IMC were instructed to collect sputum in a specimen cup, if possible (for patient safety, sputum induction was not performed). Specimen collection was repeatedly performed at a median of 3-day intervals for each patient during hospital admission.

For RT-PCR assays, PowerCheck 2019-nCoV (KogeneBiotech, Seoul, Korea) was used until 14 June 2020, and Real-Q 2019-nCoV (BioSewoom Inc., Seoul, Korea) was used from 15 June 2020. Both kits use E and RNA-dependent RNA polymerase (RdRp) genes as genetic markers. The cutoff values of the cycle threshold (Ct), which is inversely correlated with the amount of viral RNA present, were 35 for the KogeneBiotech kit and 38 for the BioSewoom kit. The test result was interpreted as positive only when the Ct values of both genes were lower than the cutoff value.

Statistical analysis

Categorical variables are presented as frequencies and proportions. Continuous variables are summarized as medians and IQRs. Differences in means were compared using Student’s t-test, whereas differences in medians were compared using the Kruskal–Wallis rank-sum test. The sequential trends of Ct values were based on the patients’ demographic and clinical characteristics, including (a) age, (b) changes in clinical severity status between hospital admission and discharge, (c) number and types of clinical symptoms that appeared after hospital admission, and (d) whether negative conversion was delayed or not. Sequential data of viral loads obtained from both URT and LRT specimens were aggregated using the cubic spline smoothing function of a generalized additive model: Log(μwhere t refers to the day of observation, Y denotes the observed Ct value at time t, s denotes a smoothing function, and time denotes the number of days from symptom onset. The time at which the peak viral load was observed and the related confidence interval (CI) were estimated through breakpoint analysis, which detects a change point in spline curves by fitting piecewise linear regressions (Muggeo, 2003, Muggeo, 2017). This analysis indicated the time when the viral load reached a peak and the time when it no longer decreased by a statistically significant amount. Survival analysis (Kaplan–Meier method) and log-rank tests were used to compare the time from hospital admission to negative conversion between the different symptom groups. The symptom groups were defined according to (a) the presence of symptoms (i.e., asymptomatic, presymptomatic, and symptomatic) and (b) the number and types of symptoms.

A P-value of < 0.05 was considered to indicate statistical significance when interpreting the results. All statistical analyses were performed using R software (version 3.6.3).

Results

Clinical characteristics and prognosis

The median age of the patients was 38.0 (12.0–68.0) years. Nearly half (46.4%) of the patients were aged 20–39 years, followed by those aged 40–59 years (34.8%); 67 (48.6%) patients were female. Of the 138 patients, seven were referred to a tertiary hospital because of symptom exacerbation, 131 were discharged, and none died (Table 1 ). Six patients had at least one comorbidity. The patients’ comorbidities included obesity (N = 4, 2.9%), diabetes without complications (N = 3, 2.2%), and hypertension (N = 3, 2.2%); one patient had both obesity and hypertension. Seven patients were current smokers (5.1%). Fever (N = 57, 41.3%) and cough (N = 39, 28.3%) were the most common symptoms at hospital admission, followed by sore throat (N = 22, 15.9%) and sputum production (N = 17, 12.3%) (Table 1). The median (IQR) duration between symptom onset and admission was 2.0 (0.0–7.0) days, and that between symptom onset and the last clinical symptom was 15.0 (11.0–25.0) days. The median (IQR) durations from hospital admission to the first and last negative RT-PCR results were 18.0 (13.0–24.0) and 19.0 (14.0–26.0) days, respectively. In general, the median (IQR) duration from hospital admission to discharge was 19.5 (15.0–26.0) days, and the median (IQR) number of RT-PCR tests per patient was 8.0 (6.0–12.0), with a 3-day median test interval (Table 1).

Table 1

Baseline characteristics of the patients.

	Number of patients(N = 138)
Age, years
Median (minimum–maximum)	36.0 (12.0–68.0)
0–9	0 (0.0%)
10–19	13 (9.4%)
20–39	64 (46.4%)
40–59	48 (34.8%)
≥60	13 (9.4%)
Sex
Female	67 (48.6%)
Male	71 (51.4%)
Symptoms at hospital admission
History of fever	57 (41.3%)
Cough	39 (28.3%)
Sore throat	22 (15.9%)
Sputum production	17 (12.3%)
Headache	14 (10.1%)
Runny nose	12 (8.7%)
Myalgia	11 (8.0%)
Diarrhea	4 (2.9%)
Chest pain	2 (1.4%)
Dyspnea	2 (1.4%)
Abdominal pain	1 (0.7%)
Duration (median, IQR)
Days from exposure to symptom onset	4.0 (3.0–7.0)
Days from symptom onset to admission	2.0 (0.0–7.0)
Days from symptom onset to the last symptom	15.0 (11.0–25.0)
Days from admission to the first negative result of two consecutive RT-PCR tests	18.0 (13.0–24.0)
Days from admission to the second negative of two consecutive RT-PCR tests	19.0 (14.0–26.0)
Days from admission to discharge	19.5 (15.0–26.0)
SARS-CoV-2 RT-PCR assay
Total tests (oronasopharyngeal sample)	1,286
Total tests (sputum sample)	1,272
Number of samples per patient (median, IQR)	8.0 (6.0–12.0)
Clinical severity score (maximum score during hospitalization)
1 No limitation in activity	37 (26.8%)
2 Limitation in activity but no O2 requirement (temperature ≥ 37.5 °C)	76 (55.1%)
3 O2 administration with a nasal prong	15 (10.9%)
4 O2 administration with a facial mask	3 (2.2%)
5 Noninvasive ventilation (high-flow O2, FiO2 > 0.4)	7 (5.1%)
6 Invasive ventilation	0 (0.0%)
7 Multi-organ failure or ECMO	0 (0.0%)
8 Death	0 (0.0%)

IQR, interquartile range; ECMO, extracorporeal membrane oxygenation. The clinical severity scores were based on the eight-category ordinal scale developed by Sung et al. (2020), which was modified from the scale developed by Marshall et al. (2020).

Of the 138 patients, 23 (16.7%) were asymptomatic at hospital admission, 16 (11.6%) developed symptoms later (presymptomatic), and three showed progression from a moderate to severe state (Figure 1 ). Of the 111 patients who were mildly symptomatic during the first 24 h of admission, 18 showed progression to a moderate or severe state (13.1%). In summary, most (N = 106, 76.8%) of the patients were classified as having mild disease, defined as having no activity limitations or having activity limitations but without oxygen requirements (Table 1). Fifteen patients (10.9%) were classified as having moderate disease, whereas 10 (7.2%) had severe disease. None of the patients included in this study died or required invasive ventilation or extracorporeal membrane oxygenation (ECMO).

Figure 1

Disease progression in the study patients.

Clinical severity was classified according to the eight-category ordinal scale developed by Sung et al. (2020). “Mild” means 1) without limitations in daily activities or 2) with limitations in daily activities but without a need for supplemental oxygen therapy. “Moderate to severe” refers to the need for 3) supplemental oxygen therapy through a nasal cannula, 4) supplemental oxygen therapy through a facial mask, or 5) high-flow supplemental oxygen therapy or noninvasive mechanical ventilation. None of the patients in this study 6) needed invasive ventilation, 7) multi-organ failure or extracorporeal membrane oxygenation (ECMO) or 8) died.

Viral load kinetics of SARS-CoV-2 infection

A total of 2,558 samples from 138 patients were collected, comprising 1,286 URT samples and 1,272 sputum samples. The mean Ct value was lower in the LRT samples than in the URT samples, indicating a higher viral load in the LRT samples; however, the difference was not statistically significant for both E and RdRp genes (Supplementary Table 1). The differences in mean Ct values of were not significant between the symptomatic (including presymptomatic patients) and asymptomatic groups (Supplementary Table 2).

Spline curves suggested that in the URT samples, the E gene RNA levels peaked on day 4 (95% CI 2.6–5.4 days) after symptom onset and rapidly decreased until day 10 (95% CI 9.2–11.9 days) (Figure 2 A). Similar patterns were observed for the RdRp gene in the URT samples (Figure 2B). In contrast, the viral RNA levels in sputum samples immediately peaked from disease onset and rapidly decreased until 15.6 days (95% CI 13.2–18.0 days) after symptom onset for the E gene and 15.0 days (95% CI 12.6–17.3 days) after symptom onset for the RdRp gene (Figure 2C and D).

Figure 2

Time trends of viral loads for the E and RdRp genes in upper and lower respiratory tract specimens from all patients.

Red numbers indicate the time point in days (95% confidence interval) of the spline curve peak to the lowest cycle threshold value (inversely related to the highest viral load). Blue numbers indicate the time point when the viral load shows no further significant decrease. For data consistency, only real-time polymerase chain reaction assays performed up to 14 June 2020 using the same type of test kit were included.

The temporal trends and peak viral loads differed according to age group (Figure 3 A). For both URT and LRT specimens, the viral load observed in the early phase of the disease course was higher in patients aged ≥ 60 years than in those < 60 years (Figure 3A). The peak viral load occurred 2–3 days after symptom onset in mild and stable patients, whereas patients with aggravated symptoms during hospitalization (mild to severe) showed an immediate peak in viral load on the day of symptom onset (Figure 3B). In addition, patients with fever, respiratory symptoms (e.g., cough, sputum production, sore throat, runny nose, and dyspnea), and other symptoms (e.g., diarrhea, headache, and vomiting) during hospitalization showed higher viral loads in the early disease course than those with respiratory symptoms alone (Figure 3C). Higher levels of SARS-CoV-2 RNA in the early course of the disease were associated with a longer time (> 17 days) from hospital admission to the first negative RT-PCR result (Figure 3D).

Figure 3

Comparisons of temporal patterns of viral load according to age group, clinical severity, clinical symptoms, and days from admission to negative conversion.

(A) Age group; (B) Changes in clinical severity score from the first 24 h of hospital admission and the maximum score during hospitalization; (C) Type of symptoms observed from hospital admission to discharge (respiratory symptoms included cough, sputum production, sore throat, runny nose, and dyspnea; other symptoms included those listed in Table 1 other than fever and respiratory symptoms); (D) Delayed time to negative conversion (17 days based on the median duration from admission to the first negative real-time polymerase chain reaction [RT-PCR] result). The dashed line refers to the cutoff value for negative RT-PCR results. Only results for the E gene from upper respiratory tract (URT) specimens are presented here, as URT specimens were collected from all patients, and E gene is used as a screening test in the RT-PCR protocols.

Comparison of time to negative conversion between the groups

The median (IQR) duration from hospital admission to negative results of two consecutive RT-PCR tests was 18.0 (13.0–24.0) days, and no statistically significant difference in duration was found according to sex (P = 0.63) and age (P = 0.58) (Supplementary Table 3).

According to the Kaplan–Meier analysis, the proportions of patients with negative conversion at days 14 and 21 from diagnosis were 57.1% and 85.2% in the asymptomatic group and 32.8% and 59.5% in the symptomatic (including presymptomatic patients) group, respectively (P = 0.10) (Figure 4 A). The median (IQR) time from hospital admission (diagnosis) to the first negative conversion was 13.0 (9.5–17.5) days in asymptomatic patients and 18.0 (13.0–25.0) days in symptomatic (including presymptomatic) patients; however, the difference was not significant (P = 0.10).

Figure 4

Comparison of proportions of negative conversion between symptom groups.

Comparison of the days from hospital admission to the first negative result of two consecutive RT-PCR results, depending on the (A) presence of symptoms and (B) type of symptoms.

When stratified by types of clinical symptoms, more than half of the patients (N = 81, 58.7%) showed at least both fever and respiratory symptoms. The median (IQR) duration from admission to the first negative RT-PCR result was longer in patients with fever, respiratory symptoms, and other symptoms (20.0 [15.0–29.0] days) than in those who had fever alone (11 [6.5–18.5] days) or respiratory symptoms alone (14.5 [12.8–18.0] days), whereas asymptomatic patients had an even shorter time to negative conversion (13 [9.5–17.5] days) (P = 0.01) (Supplementary Table 3). A Kaplan–Meier plot showed that the proportions of patients with negative conversion at days 14 and 21 from hospital admission were 20.8% and 54.2% among those with all types of symptoms and 50.0% and 87.5% in those with respiratory symptoms alone, respectively (P = 0.04) (Figure 4B).

Discussion

Although there have been numerous studies on the temporal dynamics of SARS-CoV-2 transmission, evidence is lacking on the viral load in the early days before or after the onset of symptoms and its relationship to clinical progression, owing to the time delays from symptom onset to diagnosis and hospital admission. In Korea, aggressive contact tracing, widespread diagnostic testing, and immediate quarantine since the beginning of the COVID-19 pandemic have enabled the sequential observation of asymptomatic or mildly symptomatic patients in a designated facility. This has enabled researchers to investigate the association between viral load and disease severity in mildly symptomatic patients from the early disease course.

This study revealed that the SARS-CoV-2 viral load peaked on day 4 after symptom onset for both the E and RdRp genes. Several other studies have reported that viral loads from respiratory specimens peaked a few days after symptom onset or hospital admission (He et al., 2020, Huang et al., 2020, Tan et al., 2020, Widders et al., 2020, Zhurakivska et al., 2020); however, these patterns were reported based on the aggregation of relatively small datasets. Therefore, the trends in viral load during the first days of the COVID-19 course remain uncertain. Comparisons of viral RNA levels between URT and LRT samples showed longer detection of viral RNA in LRT samples (decreased until days 15.6 and 15 for the E and RdRp genes, respectively) than in URT samples (decreased until days 10.6 and 10.5 for the E and RdRp genes, respectively), which is consistent with previous findings (Walsh et al., 2020). Additionally, whereas several previous studies have reported that sputum samples had higher viral loads and showed viral load peaks later than URT samples (Kim et al., 2020, Liu et al., 2020a, Liu et al., 2020b, Wölfel et al., 2020, Zheng et al., 2020), this study observed that the viral loads in sputum samples peaked at the time of symptom onset, and there was no difference in the mean viral load between the two sample types.

The viral RNA levels in URT and LRT samples decreased until days 10 and 15 after symptom onset, respectively, which has important implications for the release of patients from isolation. As prolonged isolation of patients and repeated diagnostic tests impose social and economic burdens, the World Health Organization and the governments of each participating country have modified the strategies for discontinuing the isolation of confirmed COVID-19 patients from a laboratory-based to a symptom-based approach. For example, France has reduced the quarantine period from 10 to 7 days (Atlani-Duault et al., 2020). Despite growing evidence showing that the period of detection of viral RNA by RT-PCR may be longer than the period of infectivity and that the viral load considerably decreases after the first week of illness, there is still a lack of quality data on the actual time when the viral load peaks and when it no longer decreases. The current findings confirmed that viral loads from respiratory specimens that peaked at 4 days after symptom onset showed a rapid decrease and stopped significantly decreasing beyond day 10 for URT specimens and day 15 for LRT specimens.

It was also discovered that a higher viral load at the beginning of the disease course could indicate aggravated clinical disease, clinical manifestation of more types of symptoms, and longer (delayed) time to negative conversion (> 17 days). Recent studies have revealed a longer time to negative conversion in patients with severe disease (Lee et al., 2020a, Xu et al., 2020); however, the current data showed no significant difference in the time to negative conversion between the mild and severe groups. This is possibly because severe or critical patients, defined as those who needed invasive mechanical ventilation or ECMO (Sung et al., 2020), were not included in this study. Patients showing progression to a severe or critical status were immediately referred to a tertiary hospital.

The duration of viral shedding was longer in symptomatic patients than in asymptomatic patients, and the time to negative conversion was significantly different between the two groups, which is consistent with the results of previous studies. (Lee et al., 2020b) reported negative conversion rates of 33.7% and 29.6% at 2 weeks after diagnosis among asymptomatic and symptomatic patients, respectively. During this period, the viral load in the symptomatic group tended to decrease more slowly than that in the asymptomatic group. Additionally, the more types of symptoms a patient had, the longer the duration to negative conversion.

This study had some limitations. First, the patients included in this study were mostly limited to those who satisfied the criteria for admission to IMC, which is a Group B medical institution where only patients with mild to moderate symptoms are admitted. Because severe cases were immediately referred to a tertiary hospital, it was difficult to investigate the association between viral load and clinical severity among patients with severe or critical disease. Second, as asymptomatic cases may be underrepresented, these findings cannot be generalized to the community at large, given that other studies have reported that 30–40% of patients in Korea were asymptomatic (Lee et al., 2020b, Workman et al., 2020).

Despite these limitations, a total of 2,558 specimens obtained from both the URT and LRT were included in this study, nearly one-fifth (N = 492, 19.2%) of which were collected within the first 5 days after symptom onset. This represents a reasonably large amount of data contributing to evidence on the early period of the disease course. This was possible because the median duration from symptom onset to hospital admission was 2 days, which seems to be the result of aggressive testing and quarantine, as mentioned above.

In summary, this study quantitatively assessed the time trends of viral loads in COVID-19 patients from symptom onset to disease recovery, by aggregating RT-PCR results obtained from 2,558 respiratory specimens from 138 patients. The viral load of COVID-19 patients peaked in the first few days (2–4 days) after symptom onset and continued to decrease until days 10–15 of the illness. Viral loads represented by Ct values could be used for monitoring clinical progression and for making decisions regarding intervention. Although the symptom-based strategy, rather than the previous laboratory-based strategy, is recommended for deciding whether to discontinue the isolation of COVID-19 patients, assessing the viral load at any point of the disease course could provide necessary evidence when determining the need for intervention or prolonged hospitalization.

Conflicts of interest

The authors have no conflicts of interest to declare.

Funding source

This project was supported by the Government-wide R&D Fund Project for Infectious Disease Research (GFID), Republic of Korea (grant No. HG18C0088).

Ethical approval

This study was approved by the institutional review board of IMC (approval no. 115288-202010-HR-048-01).

Author contributions

Conceptualization: HKC, JYK, JHK; data curation: AYL, YJO, JKL, JBS, HJK, BH, SWP, YJ, CYY, YOP, JYK; formal analysis: AYL, JHK; methodology: AYL, HKC, JYK, JHK; visualization: AYL, JHK; writing – original draft and editing: AYL, HKC, JYK, JHK.