Skip to main content

ORIGINAL RESEARCH article

Front. Public Health, 12 August 2024
Sec. Infectious Diseases: Epidemiology and Prevention

Long-term trends in cycle threshold values: a comprehensive analysis of COVID-19 dynamics, viral load, and reproduction number in South Korea

Jungeun ParkJungeun Park1Sung-il ChoSung-il Cho1Sang-Gu KangSang-Gu Kang1Jee-Woun KimJee-Woun Kim1Sunkyung JungSunkyung Jung2Sun-Hwa LeeSun-Hwa Lee2Kyou-Sup HanKyou-Sup Han2Seung-sik Hwang
Seung-sik Hwang1*
  • 1Graduate School of Public Health, Seoul National University, Seoul, Republic of Korea
  • 2Seegene Medical Foundation, Seoul, Republic of Korea

Background: With the emergence of COVID-19 cases, governments quickly responded with aggressive testing, contact tracing, isolation and quarantine measures. South Korea’s testing strategy primarily relied on real-time reverse-transcriptase polymerase chain reaction (real-time RT-PCR), focusing on cycle threshold (Ct) values, indicative of viral load, to determine COVID-19 positivity. This study examined the long-term time series distribution of Ct values measured in the same laboratory using a nationally standardized testing type and sampling method in South Korea. It aimed to link Ct values, new COVID-19 cases, and the reproduction number (Rt), setting the stage for using Ct values effectively.

Methods: This study analyzed nationally collected 296,347 samples Ct values from February 2020 to January 2022 and examined their associations with the number of new cases and Rt trends. The data were categorized into four COVID-19 periods for in-depth analysis. Statistical methods included time series trend analysis, local regression for smoothing, linear regression for association analysis, and calculation of correlation coefficients.

Results: The median Ct values across four COVID-19 periods decreased gradually from 31.71 in the initial period to 21.27 in the fourth period, indicating higher viral load. The comparison of trends between Ct values and the number of new cases revealed that the decline in Ct values preceded the surge in new cases, particularly evident during the initial stages when new cases did not undergo a significant increase. Also, during variant emergence and vaccination rollout, marked shifts in Ct values were observed. Results from linear regression analysis revealed a significant negative relationship between Ct values and new cases (β = −0.33, p < 0.001, R2 = 0.67). This implies that as Ct values decrease, new case numbers increase.

Conclusion: This study demonstrates the potential of Ct values as early indicators for predicting confirmed COVID-19 cases during the initial stages of the epidemic and suggests their relevance in large-scale epidemic monitoring, even when case numbers are similar.

1 Introduction

In South Korea, the first case of coronavirus disease 2019 (COVID-19) was confirmed on January 20, 2020 (1), and the government established its response measures based on aggressive testing, tracing, quarantine, and isolation (2, 3). Of these strategies, aggressive testing is at the core for preventing the spread of infection through the community by enabling early infection detection. To facilitate testing, the government launched multiple screening sites, implemented mobile sample collection teams and mobile health examination services, continually increased diagnostic testing counts, and added COVID-19 to the existing respiratory disease surveillance system (2). As a result, mass COVID-19 tests were done with quick turnaround times, and the government performed tracing, enforced quarantine, and provided treatment based on the test results. As such, South Korea was considered to have managed the COVID-19 pandemic well compared to other countries (4).

In South Korea, COVID-19 tests are primarily performed via the 2019 novel coronavirus real-time reverse-transcriptase polymerase chain reaction (real-time RT-PCR) test to detect viral genes, typically in upper and lower airway samples (5). The assay is designed to detect two or three regions of severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) genes that encode nucleocapsid (N), envelope (E), or spike (S) proteins or the RNA-dependent RNA polymerase (RdRp) gene (5). The cycle threshold (Ct) values, representing the amplification cycles needed to surpass a predetermined threshold, serve as a determinant for COVID-19 positivity (6). A sample is deemed positive when the Ct value of the target gene is equal to or smaller than the cutoff value (6). Because this value can also be understood as the viral load, a lower Ct value is interpreted as a higher viral load (7).

Many recent studies have utilized Ct values to establish associations between Ct values and demographic factors (8, 9). These studies have explored the correlation between Ct values and various factors such as the occurrence of confirmed cases (10), hospitalization and mortality rates (11, 12), identification of infection stages in patients (13, 14), prediction of clinical characteristics and disease severity (13, 15, 16), differences based on viral mutations and vaccination status (1719), and clinical decisions regarding the need for isolation and quarantine (7). Moreover, a study has utilized Ct values in relation to non-pharmaceutical interventions (20).

Investigating transmission trends and predicting outbreak sizes during an epidemic provide valuable data for formulating effective response measures (21). Currently, the time-varying effective reproductive number (Rt) is employed as an indicator for predicting epidemic trends (22). However, Rt is contingent on the number of confirmed cases, and a significant limitation of this index is its inability to offer real-time predictions of epidemic size due to the delay between the latent period of a particular disease and the reporting of confirmed cases. To address this limitation, some studies have applied the modeling of infection rates and the utilization of surveillance system indicators (2325) to large population distributions of Ct values, demonstrating the ability to predict epidemic trends. Additionally, others have incorporated Ct value distributions into existing Rt calculations to enhance the accuracy of predictions and compare their predictive capabilities (25).

This study examined the long-term time series distribution of Ct values measured in the same laboratory using a nationally standardized testing type and sampling method in South Korea. Furthermore, it sought to analyze the relationship between Ct values, new COVID-19 cases, and the reproduction number (Rt) through comparisons, thus establishing the basis for utilizing Ct values.

2 Materials and methods

2.1 Data

The RT-PCR test results used in this study were obtained from upper respiratory samples collected between February 2020 and January 2022 by the Korea Centers for Disease Control and Prevention. These samples were collected from various sources, including public health centers, temporary screening clinics, medical institutions’ screening clinics, and general medical clinics, spanning 14 out of 17 regions across South Korea.

The individuals tested comprised those exhibiting suspected COVID-19 symptoms, individuals with epidemiological links to confirmed cases, people aged 60 and above, as well as preemptive screening for infection-prone facilities. Positive cases from rapid antigen and emergency screening tests were also included.

Different assays were used in different periods. In the early stages, emergency use assays (February 7 to December 13, 2020: Allplex 2019-nCoV [SARS-CoV-2] Assay to detect the E, RdRP, and N gene targets) were used; the official assays (December 14, 2020, to January 31, 2022: Allplex SARS-CoV-2 Assay to detect the E, RdRP/S, and N gene targets) were used thereafter.

In this study, the total population of the sampling areas was 51,266,914, from which 37,253,185 samples were collected, of which Ct values were determined for 296,347 samples and found to be positive for the RdRp/S target. These positive cases represent 36.3% of all positive cases in South Korea.

The data concerning the number of newly confirmed COVID-19 cases were collected from January 20, 2020 (the date of the first confirmed case) up to January 2022. The number of newly confirmed cases encompasses domestically transmitted and imported cases, as reported in the daily COVID-19 report of the Integrated Management System for Epidemic Diseases provided by the Korea Disease Control and Prevention Agency. The study utilized the daily Rt data for South Korea, extracted from the COVID-19 dataset provided by “Our World In Data”1 a project of the Global Change Data Lab, a registered charity based in England and Wales. The estimation of Rt values was conducted using the Kalman filter, an established technique for Rt estimation. The specific Rt values employed in this study were derived from a research paper that applied the Kalman filter to estimate Rt for 124 countries, including South Korea (26).

2.2 Statistical analysis

In this study, we compared the time-series trends of daily median Ct values with daily the number of newly confirmed COVID-19 cases, and the Rt values. To compare the time-series trends, we employed local regression to smooth the data of Ct values, the number of newly confirmed COVID-19 cases, and Rt values. Local regression, also known as locally weighted scatterplot smoothing (LOWESS), is a non-parametric regression method used to estimate the underlying trend in a dataset by fitting multiple local regressions. Linear regression analysis was used for the correlation analysis. We introduced a time lag between Ct values and new cases and Rt for the association analysis to conduct a detailed analysis of the two indices. The analysis of time-series trends and correlations was conducted by comparing between the entire duration and four specific COVID-19 epidemic periods as delineated by the Korea Disease Control and Prevention Agency (27). The data were analyzed using the SAS 9.4 software (SAS Institute, Inc., Cary, NC, United States).

3 Results

Four COVID-19 periods occurred in South Korea between January 20, 2020, and January 31, 2022. We examined the distributions of daily Ct values, new cases, and Rt and analyzed their associations for the entire study period and the four COVID-19 periods.

Ct values declined in general during the entire period, indicating an increased viral load over time. New cases surged after a substantial drop in the Ct values. Upon comparing the periods of reagent changes, vaccination campaigns, and the emergence of variant viruses with potential influences on Ct values, it became evident that Ct values were increasing even before the reagent change. After the vaccination campaigns, a slight decrease in Ct values was observed. In relation to the emergence of variant viruses, it was observed that Ct values exhibited an increasing trend following the emergence of the Alpha, Beta, and Gamma variants. After the emergence of the Delta variant, Ct values remained consistently low with no significant fluctuations, while following the emergence of the Omicron variant, Ct values exhibited an increasing trend (Figure 1).

Figure 1
www.frontiersin.org

Figure 1. Distributions of Ct values and new cases per million population: evolution of the 1-day median Ct values (blue dot) and new cases per million (orange dot), obtained from February 2020, to January 2022 (Korea). Loess curves fitted to Ct median and new cases per million shown as a blue line (alpha = 0.75) and orange line (alpha = 0.75). Skyblue areas surrounding the blue line indicate the 95% confidence limits (CL) associated with median Ct values. The thick black vertical lines indicate the start and end of the four COVID-19 periods. The vertical red dashed lines denote changes in the testing assay and the commencement of vaccination begins. The vertical blue dashed lines indicate the onset of COVID-19 variants.

The study categorized the time-series trends into different periods, each corresponding to different waves of the COVID-19 pandemic. In the first, second, and third periods, notable case spikes were primarily caused by cluster infections within specific regions. During these periods, a clear resemblance was observed between the patterns of Ct value trends and the trends of newly confirmed cases with a particular time lag. However, in the fourth period, characterized by large-scale outbreaks due to community transmission, although there were fluctuations in Ct values, no similarities were evident in the trends of newly confirmed cases (Figure 2).

Figure 2
www.frontiersin.org

Figure 2. Distributions of Ct values and new cases per million during each COVID-19 period: evolution of the 1-day median Ct values (blue dot) and new cases per million (orange dot), obtained from February 2020, to January 2022 (Korea). Loess curves fitted to Ct median and new cases per million shown as a blue line (alpha = 0.75) and orange line (alpha = 0.75). Skyblue areas surrounding the blue line indicate the 95% confidence limits (CL) associated with median Ct values. Vivid orange areas surrounding the orange line indicate the 95% confidence limits (CL) associated with new cases per million.

The study compared the median overall Ct values for each period during the epidemic. The median Ct value was highest (31.71) in the first period and decreased over time, with the lowest value (21.27) in the fourth period (Figure 3).

Figure 3
www.frontiersin.org

Figure 3. Ct values (median) by COVID-19 period.

Next, we compared Ct value trends with Rt. The Ct values trend revealed a pattern similar to the Rt fluctuations trend (Figure 4).

Figure 4
www.frontiersin.org

Figure 4. Distributions of Ct values and Rt: evolution of the 1-day median Ct values (blue dot) and Rt (purple dot), obtained from February 2020, to January 2022 (Korea). Loess curves fitted to Ct median and Rt shown as a blue line (alpha = 0.75) and purple line (alpha = 0.75). Skyblue areas surrounding the blue line indicate the 95% confidence limits (CL) associated with median Ct values. The thick black vertical lines indicate the start and end of the four COVID-19 periods. Vivid orange areas surrounding purple line indicate the 95% confidence limits (CL) associated with Rt.The vertical red dashed lines denote changes in the testing assay and the commencement of vaccination begins. The vertical blue dashed lines indicate the onset of COVID-19 variants.

The study analyzed the relationship between Ct values, the number of newly confirmed cases, and Rt. While comparing Ct value trends with the occurrence of newly confirmed cases and Rt trends, a time lag was observed. Based on the correlation analysis between the daily median Ct values and the occurrence of newly confirmed cases, the highest correlation coefficient between the two variables was observed with a time lag of 7 days. Notably, Ct values and the number of newly confirmed cases occurring 7 days later exhibited a negative correlation (Figure 5).

Figure 5
www.frontiersin.org

Figure 5. Association between Ct values and number of new cases.

Next, we analyzed the association between Ct values and new cases per COVID-19 period. During all four periods, new cases declined with increasing Ct values (Figure 6).

Figure 6
www.frontiersin.org

Figure 6. Association between Ct values and the number of new cases during each COVID-19 period.

Based on the correlation analysis between the daily median Ct values and Rt, the highest correlation coefficient between the two variables was observed with a time lag of 1 day. Notably, the analysis of the relationship between Ct values and Rt 1 day later exhibited a negative correlation (Figure 7).

Figure 7
www.frontiersin.org

Figure 7. Association between Ct values and Rt.

In conclusion, the analysis of the relationship between Ct values and Rt, divided by epidemic periods, revealed a consistent pattern. In all periods, it was observed that as Ct values increased, Rt decreased, indicating a negative correlation between the two variables (Figure 8).

Figure 8
www.frontiersin.org

Figure 8. Association between Ct values and Rt during each COVID-19 period.

4 Discussion

Since the first case of COVID-19 in the country on January 20, 2020, South Korea established an aggressive screening testing strategy to ensure early detection of individuals who developed the infection. Patients were traced, quarantined, and treated as necessary. Based on these implemented strategies, South Korea is considered to have done well in COVID-19 management compared to other countries. COVID-19 testing was performed via real-time RT-PCR, with Ct values, and the number of amplification cycles needed to cross a certain threshold, used to determine positivity. This value can also be understood as corresponding to the patient’s viral load.

Recently, there has been a shift away from the binary use of Ct values for determining COVID-19 positivity toward research that explores their utilization in understanding demographic variations, correlations with hospitalization and mortality rates, identification of infection stages in patients, and their relevance to non-pharmaceutical interventions (820). Moreover, some studies have delved into applying Ct values to predict future epidemic trends and monitor infectious diseases (2325).

In this current study, we sought to examine the long-term time series distribution of Ct values measured in the same laboratory using nationally standardized testing types and sampling methods in South Korea. Additionally, we aimed to contribute to the establishment of the basis for utilizing Ct values by analyzing their relationship through comparisons with new COVID-19 cases and the Rt.

From January 20, 2020, to January 31, 2022, there were four COVID-19 periods in Korea. Here, we analyzed the associations between daily Ct values, the number of new cases, and the Rt for the entire COVID-19 period and each COVID-19 period. Ct values generally decreased during the entire period analyzed, indicating an increased viral load over time. We also observed that new cases surged after a substantial drop in Ct values. Regarding the four different COVID-19 periods, the average daily number of new cases was markedly higher during the second period (142.8) than during the first period (71.5), even though the cumulative number of cases was similar between the two periods, showing that the virus spreads more rapidly when Ct values are low.

The trajectory of Ct values displayed a resemblance to the trend observed in Rt, a predictive indicator of incidence. Analysis of their association with new cases, incorporating time lag, revealed that Ct values anticipate the onset of cases marginally sooner than Rt.

Exploring the factors influencing Ct values, notably variant emergence and vaccination endeavors, delineated distinct shifts pre- and post-events. The ascendancy of the Alpha variant in December 2020 marked the onset of variant dominance, persisting until mid-July. The subsequent emergence of the Beta and Gamma variants in January 2021 coincided with an observable uptrend in Ct values. The Delta variant’s dominance from April 2021 to December 2021 ushered in a phase of consistently low Ct values devoid of significant fluctuations. However, the emergence of the Omicron variant in late November 2021 heralded a noticeable uptick in Ct values. Further analysis underscored a decline in Ct values following the inception of vaccination efforts.

Several limitations need to be acknowledged in our study. Firstly, while our study data were collected at the national level, they only represented 40% of the total data, with test results missing for two out of 17 cities and provinces in Korea. This incomplete coverage may introduce potential biases.

Furthermore, there are inherent limitations associated with utilizing Ct values and confirmed case occurrences as supplementary variables. The lack of precise timing for virus infection and sampling during specimen collection for Ct value determination introduces uncertainty. This uncertainty regarding the exact timing of virus infection and sampling points may hinder our ability to fully and accurately account for differences in Ct values and the occurrence of confirmed cases.

In addition, our use of daily median values for Ct values has its own limitations. This approach does not consider the sampling variability inherent in the median itself nor explores various trends and analyses that could account for factors like data skewness.

Despite these limitations, our study’s data, collected at the national level of viral load, holds potential applicability in diverse settings, especially within the context of surveillance. This dataset could prove valuable in nationwide applications, including cluster epidemiological investigations and implementing COVID-19 control measures such as testing, contact tracing, and isolation/quarantine at the national level.

Moreover, there is a pressing need for comparative studies examining the impact of non-pharmaceutical interventions, such as mask-wearing and social distancing, on Ct values, the incidence of newly confirmed cases, and the effective Rt, depending on the timing of their implementation.

Finally, the significant association between Ct values and COVID-19 case occurrence 1 week later, as revealed by our study, underscores the importance of adjusting for potential confounding effects of confirmed cases from 1 week prior when analyzing these variables. Such adjustments will facilitate a clearer understanding of the relationship between Ct values and confirmed cases and highlight their utility as predictive factors.

As observed in this study, the negative correlation between Ct values and the number of new cases in the overall population is consistent with previous reports (11, 24, 28). Similar correlations were observed in studies focusing on hospital settings and local populations (29, 30). Furthermore, Avadhanula et al. (29) found that the correlation between Ct values and the number of new cases after 1 week was strongest among the factors studied.

In this study, by analyzing the trends of Ct values measured in the same laboratory using nationally standardized testing types and sampling methods, we were able to confirm the potential utility of Ct values for predicting the occurrence of confirmed cases during the early stages of the epidemic when the number of cases is low. Furthermore, even when the scale of confirmed cases is similar, we discovered variations in the extent of virus spread based on Ct values. As suggested by Yin et al. (31) and Musalkova et al. (32), this study also proposes the utilization of large-scale epidemic and infectious disease monitoring incorporating Ct value information.

Data availability statement

The datasets presented in this article are not readily available because the access rights to some of the datasets used in this study are owned by “Seegene Medical Foundation.” Requests to access the datasets should be directed to SL, lshkim@mf.seegene.com.

Ethics statement

The studies involving humans were approved by the Institutional Review Board (IRB) of Seoul National University (IRB No. E2207/001-002). The studies were conducted in accordance with the local legislation and institutional requirements. The ethics committee/institutional review board waived the requirement of written informed consent for participation from the participants or the participants’ legal guardians/next of kin because The institution providing the data has already submitted approved IRB documents.

Author contributions

JP: Writing – original draft, Writing – review & editing. SC: Writing – original draft, Writing – review & editing. S-GK: Writing – review & editing. J-WK: Writing – review & editing. SJ: Writing – review & editing. SL: Writing – review & editing. K-SH: Writing – review & editing. SH: Writing – original draft, Writing – review & editing.

Funding

The author(s) declare that financial support was received for the research, authorship, and/or publication of this article. This study was supported by the Seegene Medical Foundation (grant 200-20210107).

Conflict of interest

The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.

Publisher’s note

All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.

Abbreviations

Ct, Cycle threshold; E, Envelope protein; LOWESS, Locally weighted scatterplot smoothing; N, Nucleocapsid; Rt, Reproduction number; RdRp, RNA-dependent RNA polymerase; Real-time RT-PCR, Real-time reverse-transcriptase polymerase chain reaction; S, Spike protein

Footnotes

References

1. Korea Disease Control and Prevention Agency. Korea disease control and prevention agency. (2022). Available at: https://www.kdca.go.kr (Accessed May 9, 2023).

Google Scholar

2. Park, Y, Huh, IS, Lee, J, Kang, CR, Cho, S, Ham, HJ, et al. Application of testing-tracing-treatment strategy in response to the COVID-19 outbreak in Seoul, Korea. J Korean Med Sci. (2020) 35:e396. doi: 10.3346/jkms.2020.35.e396

PubMed Abstract | Crossref Full Text | Google Scholar

3. The Goverment of the Republic of Korea.Tackling COVID-19: Health, quarantine and economic measures of South Korea. (2020).

Google Scholar

4. Lee, JY. (2020).Testing-tracing-treatment (3T) strategy for COVID-19 response: Field experience and research. Seoul National University. Available at: https://www.snu.ac.kr/coronavirus/research?md=v&bbsidx=128906 (Accessed November 30, 2022).

Google Scholar

5. Lim, J, and Lee, J. Current laboratory diagnosis of coronavirus disease 2019. Korean J Intern Med. (2020) 35:741–8. doi: 10.3904/kjim.2020.257

PubMed Abstract | Crossref Full Text | Google Scholar

6. al Bayat, S, Mundodan, J, Hasnain, S, Sallam, M, Khogali, H, Ali, D, et al. Can the cycle threshold (Ct) value of RT-PCR test for SARS CoV2 predict infectivity among close contacts? J Infect Public Health. (2021) 14:1201–5. doi: 10.1016/j.jiph.2021.08.013

PubMed Abstract | Crossref Full Text | Google Scholar

7. Tom, MR, and Mina, MJ. To interpret the SARS-CoV-2 test, consider the cycle threshold value. Clin Infect Dis. (2020) 71:2252–4. doi: 10.1093/cid/ciaa619

Crossref Full Text | Google Scholar

8. Ade, C, Pum, J, Abele, I, Raggub, L, Bockmühl, D, and Zöllner, B. Analysis of cycle threshold values in SARS-CoV-2-PCR in a long-term study. J Clin Virol. (2021) 138:104791. doi: 10.1016/j.jcv.2021.104791

PubMed Abstract | Crossref Full Text | Google Scholar

9. Mahallawi, WH, Alsamiri, AD, Dabbour, AF, Alsaeedi, H, and Al-Zalabani, AH. Association of Viral Load in SARS-CoV-2 patients with age and gender. Front Med (Lausanne). (2021) 8:608215. doi: 10.3389/fmed.2021.608215

PubMed Abstract | Crossref Full Text | Google Scholar

10. Valent, F, Mazzilis, G, Doimo, A, Marzinotto, S, Bulfoni, M, Moretti, V, et al. High viral load in upper respiratory specimens of persons with SARS-CoV-2 infection in a northern Italian area. Infect Dis Clin Pract (Baltim Md). (2021) 29:e356–60. doi: 10.1097/IPC.0000000000001068

PubMed Abstract | Crossref Full Text | Google Scholar

11. Tso, CF, Garikipati, A, Green-Saxena, A, Mao, Q, and Das, R. Correlation of population SARS-CoV-2 cycle threshold values to local disease dynamics: exploratory observational study. JMIR Public Health Surveill. (2021) 7:e28265. doi: 10.2196/28265

PubMed Abstract | Crossref Full Text | Google Scholar

12. Dehesh, P, Baradaran, HR, Eshrati, B, Motevalian, SA, Salehi, M, and Donyavi, T. The relationship between population-level SARS-CoV-2 cycle threshold values and trend of COVID-19 infection: longitudinal study. JMIR Public Health Surveill. (2022) 8:e36424. doi: 10.2196/36424

PubMed Abstract | Crossref Full Text | Google Scholar

13. Jacot, D, Greub, G, Jaton, K, and Opota, O. Viral load of SARS-CoV-2 across patients and compared to other respiratory viruses. Microbes Infect. (2020) 22:617–21. doi: 10.1016/j.micinf.2020.08.004

PubMed Abstract | Crossref Full Text | Google Scholar

14. Weinbergerova, B, Mayer, J, Hrabovsky, S, Novakova, Z, Pospisil, Z, Martykanova, L, et al. COVID-19’s natural course among ambulatory monitored outpatients. Sci Rep. (2021) 11:10124. doi: 10.1038/s41598-021-89545-1

PubMed Abstract | Crossref Full Text | Google Scholar

15. Suh, HJ, Kim, DH, Heo, EY, Lee, HW, Lee, JK, Lee, CS, et al. Clinical characteristics of COVID-19: clinical dynamics of mild severe acute respiratory syndrome coronavirus 2 infection detected by early active surveillance. J Korean Med Sci [Internet]. (2020) 35:e297. doi: 10.3346/jkms.2020.35.e297

Crossref Full Text | Google Scholar

16. Marks, M, Millat-Martinez, P, Ouchi, D, Roberts, C, Alemany, A, Corbacho-Monné, M, et al. Transmission of COVID-19 in 282 clusters in Catalonia, Spain: a cohort study. Lancet Infect Dis. (2021) 21:629–36. doi: 10.1016/S1473-3099(20)30985-3

PubMed Abstract | Crossref Full Text | Google Scholar

17. Singanayagam, A, Hakki, S, Dunning, J, Madon, KJ, Crone, MA, Koycheva, A, et al. Community transmission and viral load kinetics of the SARS-CoV-2 delta (B.1.617.2) variant in vaccinated and unvaccinated individuals in the UK: a prospective, longitudinal, cohort study. Lancet Infect Dis. (2022) 22:183–95. doi: 10.1016/S1473-3099(21)00648-4

PubMed Abstract | Crossref Full Text | Google Scholar

18. Earnest, R, Uddin, R, Matluk, N, Renzette, N, Siddle, KJ, Loreth, C, et al. Comparative transmissibility of SARS-CoV-2 variants Delta and alpha in New England, USA. medRxiv. (2021). doi: 10.1101/2021.10.06.21264641

Crossref Full Text | Google Scholar

19. Wollschläger, P, Todt, D, Gerlitz, N, Pfaender, S, Bollinger, T, Sing, A, et al. SARS-CoV-2 N gene dropout and N gene Ct value shift as indicator for the presence of B.1.1.7 lineage in a commercial multiplex PCR assay. Clin Microbiol Infect. (2021) 27:1353.e1–5. doi: 10.1016/j.cmi.2021.05.025

PubMed Abstract | Crossref Full Text | Google Scholar

20. Satlin, MJ, Zucker, J, Baer, BR, Rajan, M, Hupert, N, Schang, LM, et al. Changes in SARS-CoV-2 viral load and mortality during the initial wave of the pandemic in new York City. PLoS One. (2021) 16:e0257979. doi: 10.1371/journal.pone.0257979

PubMed Abstract | Crossref Full Text | Google Scholar

21. Das, A. An approximation-based approach for periodic estimation of effective reproduction number: a tool for decision-making in the context of coronavirus disease 2019 (COVID-19) outbreak. Public Health. (2020) 185:199–201. doi: 10.1016/j.puhe.2020.06.047

PubMed Abstract | Crossref Full Text | Google Scholar

22. Lin, Y, Yang, B, Cobey, S, Lau, EHY, Adam, DC, Wong, JY, et al. Incorporating temporal distribution of population-level viral load enables real-time estimation of COVID-19 transmission. Nat Commun. (2022) 13:1155. doi: 10.1038/s41467-022-28812-9

PubMed Abstract | Crossref Full Text | Google Scholar

23. Andriamandimby, SF, Brook, CE, Razanajatovo, N, Randriambolamanantsoa, TH, Rakotondramanga, JM, Rasambainarivo, F, et al. Cross-sectional cycle threshold values reflect epidemic dynamics of COVID-19 in Madagascar. Epidemics. (2022) 38:100533. doi: 10.1016/j.epidem.2021.100533

PubMed Abstract | Crossref Full Text | Google Scholar

24. Walker, AS, Pritchard, E, House, T, Robotham, JV, Birrell, PJ, Bell, I, et al. Ct threshold values, a proxy for viral load in community SARS-CoV-2 cases, demonstrate wide variation across populations and over time. eLife. (2021) 10:e64683. doi: 10.7554/eLife.64683

PubMed Abstract | Crossref Full Text | Google Scholar

25. Hay, JA, Kennedy-Shaffer, L, Kanjilal, S, Lennon, NJ, Gabriel, SB, Lipsitch, M, et al. Estimating epidemiologic dynamics from cross-sectional viral load distributions. Science. (2021) 373:eabh0635. doi: 10.1126/science.abh0635

Crossref Full Text | Google Scholar

26. Arroyo-Marioli, F, Bullano, F, Kucinskas, S, and Rondón-Moreno, C. Tracking R of COVID-19: a new real-time estimation using the Kalman filter. PLoS One. (2021) 16:e0244474. doi: 10.1371/journal.pone.0244474

Crossref Full Text | Google Scholar

27. Yang, S, Jang, J, and Park, SY. Two-year report of COVID-19 outbreak from January 20, 2020 to January 19, 2022 in the Republic of Korea. Korea Disease Control and Prevention Agency (KDCA); (2022).

Google Scholar

28. Mishra, B, Ranjan, J, Purushotham, P, Saha, S, Payal, P, Kar, P, et al. High proportion of low cycle threshold value as an early indicator of COVID-19 surge. J Med Virol. (2022) 94:240–5. doi: 10.1002/jmv.27307

PubMed Abstract | Crossref Full Text | Google Scholar

29. Avadhanula, V, Nicholson, EG, Ferlic-Stark, L, Piedra, FA, Blunck, BN, Fragoso, S, et al. Viral load of severe acute respiratory syndrome coronavirus 2 in adults during the first and second wave of coronavirus disease 2019 pandemic in Houston, Texas: the potential of the Superspreader. J Infect Dis. (2021) 223:1528–37. doi: 10.1093/infdis/jiab097

Crossref Full Text | Google Scholar

30. Calistri, P, Danzetta, ML, Amato, L, Cito, F, Di Giuseppe, A, Zenobio, V, et al. Epidemiological significance of SARS-CoV-2 RNA dynamic in Naso-pharyngeal swabs. Microorganisms. (2021) 9:1264. doi: 10.3390/microorganisms9061264

Crossref Full Text | Google Scholar

31. Yin, N, Dellicour, S, Daubie, V, Franco, N, Wautier, M, Faes, C, et al. Leveraging of SARS-CoV-2 PCR Cycle cycle thresholds values to forecast COVID-19 trends. Front Med. (2021) 8:743988. doi: 10.3389/fmed.2021.743988

Crossref Full Text | Google Scholar

32. Musalkova, D, Piherova, L, Kwasny, O, Dindova, Z, Stancik, L, Hartmannova, H, et al. Trends in SARS-CoV-2 cycle threshold values in the Czech Republic from April 2020 to April 2022. Sci Rep. (2023) 13:6156. doi: 10.1038/s41598-023-32953-2

PubMed Abstract | Crossref Full Text | Google Scholar

Keywords: cycle threshold value, COVID-19, reverse transcription polymerase chain reaction, testing, viral load

Citation: Park J, Cho S, Kang S-G, Kim J-W, Jung S, Lee SH, Han K-S and Hwang S (2024) Long-term trends in cycle threshold values: a comprehensive analysis of COVID-19 dynamics, viral load, and reproduction number in South Korea. Front. Public Health. 12:1394565. doi: 10.3389/fpubh.2024.1394565

Received: 01 March 2024; Accepted: 25 July 2024;
Published: 12 August 2024.

Edited by:

Reza Lashgari, Shahid Beheshti University, Iran

Reviewed by:

Muhammad Fawad Rasool, Bahauddin Zakariya University, Pakistan
Nicolas Yin, Laboratoire Hospitalier Universitaire de Bruxelles (LHUB-ULB), Belgium

Copyright © 2024 Park, Cho, Kang, Kim, Jung, Lee, Han and Hwang. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.

*Correspondence: Seung-sik Hwang, cyberdoc@snu.ac.kr

Disclaimer: All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article or claim that may be made by its manufacturer is not guaranteed or endorsed by the publisher.