Epidemiological application of the cycle threshold value of RT-PCR for estimating infection period in cases of SARS-CoV-2

Article information

J Med Life Sci. 2023;20(3):107-114
Publication date (electronic) : 2023 September 30
doi : https://doi.org/10.22730/jmls.2023.20.3.107
1WHO Health Emergencies Programme, WHO Western Pacific Regional Office, Manila, Philippines
2Jeju Center for Infectious Diseases Control and Prevention, Jeju, Republic of Korea
Correspondence to Jong-Myon Bae Jeju Center for Infectious Diseases Control and Prevention, 102 Jejudaehak-ro, Jeju 63243, Republic of Korea Tel: 82-64-755-5567 Fax: 82-64-725-2593 E-mail: jmbae@jejunu.ac.kr
Received 2023 May 18; Revised 2023 July 30; Accepted 2023 September 4.

Abstract

Epidemiological control of coronavirus disease 2019 (COVID-19) is needed to estimate the infection period of confirmed cases and identify potential cases. The present study, targeting confirmed cases for which the time of COVID-19 symptom onset was disclosed, aimed to investigate the relationship between intervals (day) from symptom onset to testing the cycle threshold (CT) values of real-time reverse transcription-polymerase chain reaction. Of the COVID-19 confirmed cases, those for which the date of suspected symptom onset in the epidemiological investigation was specifically disclosed were included in this study. Interval was defined as the number of days from symptom onset (as disclosed by the patient) to specimen collection for testing. A locally weighted regression smoothing (LOWESS) curve was applied, with intervals as explanatory variables and CT values (CTR for RdRp gene and CTE for E gene) as outcome variables. After finding its non-linear relationship, a polynomial regression model was applied to estimate the 95% confidence interval values of CTR and CTE by interval. The application of LOWESS in 331 patients identified a U-shaped curve relationship between the CTR and CTE values according to the number of interval days, and both CTR and CTE satisfied the quadratic model for interval days. Active application of these results to epidemiological investigations would minimize the chance of failing to identify individuals who are in contact with COVID-19 confirmed cases, thereby reducing the potential transmission of the virus to local communities.

INTRODUCTION

Coronavirus disease 2019 (COVID-19), officially named by the World Health Organization (WHO) on February 11, 2020, is an infectious respiratory disease caused by severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2), a single-stranded ribonucleic acid genome virus [1,2].

Regarding the main epidemiological characteristics of COVID-19, its infection period begins even before symptom onset, and the highest viral load is found at an early stage [3,4]. Accordingly, the Korea Disease Control and Prevention Agency (KDCA) has implemented testing and self-isolation for persons who were in close contact 2 days prior to the clinical symptom onset of newly confirmed cases during the epidemiological investigation [1,5]. In addition, the suspected symptoms of COVID-19 are non-specific, and some cases show mild symptoms [6,7]. Moreover, in some cases, patients are asymptomatic, and the disease is only confirmed after testing for other reasons. During the epidemiological investigation, some patients started to recollect that they manifested COVID-19-related symptoms, while others were asymptomatic [8]. If confirmed COVID-19 patients fail to recognize the suspected symptoms, it becomes a substantial hindrance in the efforts to control local community transmission. Thus, a rapid and accurate test is of the utmost importance [9,10].

To confirm COVID-19, the WHO and the KDCA recommend the use of real-time reverse transcription-polymerase chain reaction (RT-PCR) to amplify RNA-dependent RNA polymerase (RdRp) and envelope (E) genes of SARS-CoV-2 [1,4,7,9,11,12], in which cycle threshold (CT) values for the two genes are obtained [7,13]. The RT-PCR CT values determine the viral load of the corresponding specimen [7,14,15]. Thus, the CT values can be used to estimate the infection period of the confirmed cases, which can then be used for the epidemiological investigation [13,16,17]. Therefore, the symptom onset time of confirmed cases who identified themselves as asymptomatic using the CT values can be assumed. The present study aimed to determine the relationship between the intervals (days) from symptom onset to the testing of CT values using RT-PCR, with confirmed cases that specifically disclosed the time of symptom onset. Once it becomes possible to predict the time of symptom onset more accurately in individual confirmed cases using the derived relationship, highly contagious contacts of patients can be completely identified, and viral transmission to local communities can be effectively prevented.

METHODS

Of the confirmed cases in Jeju Province for which the nasal or oral specimens for COVID-19 testing were collected before or on March 9, 2021, individuals who specifically disclosed the day of the onset of suspected symptoms during the epidemiological investigation were included in the present study. From basic and in-depth epidemiological reports, data for sex, date of birth, symptom onset day, specimen collection day, and CT values of RdRp and E genes (CTR and CTE, respectively) were obtained. The RT-PCR was conducted at the Jeju Special Self-Governing Province Institute of Environment Research, a test institution officially designated by the KDCA.

Following the calculation of age based on birth date and collection day, the cases were categorized into the following four groups based on patient age: ≤19, 20-39, 40-59, and ≥60 years. The interval was defined as the number of days between the onset of COVID-19 suspected symptoms based on the information provided by the confirmed cases and specimen collection. To apply the regression analysis to the estimation, interval days with only one case were excluded from the analysis.

The Pearson correlation coefficient was calculated for the correlation analysis between CTR and CTE. A paired t-test was used to examine any statistical differences between the CTR and CTE that were paired by case. To test CTR and CTE differences by sex, t-tests and Mann-Whitney tests were used. In addition, analysis of variance and the Kruskal-Wallis test were used to analyze CTR and CTE differences by age group and number of interval days. Through such univariate analyses, variables with a 10% or lower P-value were defined as potential explanatory variables, followed by a general regression analysis using CTR and CTE values as outcome variables. A non-linear relationship between the interval days and CTR and CTE was tested using a locally weighted regression smoothing (LOWESS) curve [18]. In the case of a U-shaped parabolic curve, the polynomial regression model was applied with a quadratic term [19], and then the 95% confidence interval (CI) values of the CTR and CTE by interval day were estimated using the derived polynomial regression equation [20]. These analyses were performed using STATA 15 (StataCorp, College Station, TX, USA) statistical program, in which the statistical significance level was set as 0.05. This study was approved by the Institutional Review Board of Jeju National University (approval No. 2021-034).

RESULT

Of the 335 cases that satisfied the inclusion criteria, there was only one case each for the 11, 12, 13, and 14 interval days from clinical symptom onset until the test date. Thus, these four cases were excluded, resulting in a total of 331 subjects. The Pearson correlation coefficient between CTR and CTE in each case was 0.985, and the paired t-test revealed a significant difference (P<0.001). Thus, a statistical analysis was performed for each CTR and CTE.

Distributions of CTR and CTE according to sex, age, and interval are shown in Table 1. There was no difference in CTR and CTE distribution according to sex (150 men and 181 women subjects) (P=0.306 for CTR; P=0.177 for CTE). The distributions of CTR and CTE in the four age groups exhibited borderline statistical significance (P=0.072 in CTR and P=0.066 in CTE). In addition, the CTR and CTE distributions by interval showed statistically significant differences (P<0.001). Thus, a general regression analysis was applied to the CTR and CTE outcome values, with each interval day as an independent variable, with adjustment for age group.

Distribution of the cycle threshold in the RdRp gene (CTR) and E gene (CTE)

The LOWESS curves between the interval days and CTR and CTE showed U-shaped relationships, and both had inflection points on the second day of the interval (Fig. 1). Accordingly, a polynomial regression analysis was performed by applying a quadratic term, as shown in Table 2, in which a statistically significant contribution from the nonlinear terms was found. The predicted values of the 95% CI for the CTR and CTE values by interval based on the regression models are shown in Table 3 and Figure 2.

Figure 1.

The locally weighted regression smoothing (LOWESS) curves of critical values (yhatlowess) of the RdRp gene (CRT) (A) and critical values (yhat) of the E gene (CRE) (B) by interval showing a U-shape pattern.

Outputs of polynomial regression analysis applying the quadratic term of intervals

Predicted confidence intervals of cycle threshold of the RdRp gene (CTR) and the E gene (CTE) by intervals based on the data found in Table 2

Figure 2.

Predicted confidence intervals of critical values of the RdRp gene (CRT) (A) and E gene (CRE) (B) by interval are displayed.

DISCUSSION

In cases where the interval between the onset of COVID-19 suspected symptoms and specimen collection was identified during the epidemiological investigation, a polynomial regression equation was generated to determine the relationship between the CTR and CTE values by interval after adjusting for age groups. As shown in Table 2, the quadratic polynomial regression equations were generated for each of CTR and CTE, based on which 95% CI values of estimated CTR and CTE values by interval were obtained as in Table 3 and Figure 2.

In an analysis of CT values for spike genes in 1,077 confirmed cases of COVID-19 cases in Germany, Ade et al. [21] reported that older age was associated with a higher mean CT value, which was particularly prominent in patients aged ≥80 years. In contrast, the present study found that the CT values of the RdRp and E genes were lower when the age of patients increased (Table 1). Further studies are needed to determine whether these differences could be attributed to differences in the genes used for testing.

To control COVID-19, a definitive test for SARS-CoV-2 is important for the following reasons. First, early detection and isolation of an infected individual is a core measure for preventing the spread of SARS-CoV-2 because the virus has rapid viral shedding during the early stages of infection [12,22]. Second, the introduction of an accurate test is essential to prevent the global spread of COVID-19, as some patients clinically present with no specific symptoms or are asymptomatic [7,23]. Lastly, COVID-19 symptoms are similar to those of seasonal influenza; hence, screening during the early stages of COVID-19 is necessary, especially for individuals without specific symptoms [24].

One of the advantages of the present study is its attempt to establish a relationship between epidemiological information and quantitative RT-PCR values, leading to the identification of estimated CTR and CTE CIs by interval, as shown in Table 3, which can be used for epidemiological investigations for the control of COVID-19 as follows.

First, it can be used in cases in which one of the CTR and CTE values does not fall within the CIs of the CTR and CTE of the interval that the patients have disclosed. For example, a patient wakes up in the morning and manifests and recognizes COVID-19 symptoms when he/ she decides to be tested and has his/her specimen collected in the afternoon. In this hypothetical case, the results were 22.0 each for the CTR and CTE. As shown in Table 3, the CTE was within the CIs of interval day 0, whereas the CTR was out of the range. Thereafter, the epidemiological investigation team (EIT) determined that the symptoms were not recognized by the patient until the next day due to sleep. A re-evaluation was performed on interval day 1 when both CTE and CTR values were between the CIs. Based on these outcomes, the EIT can then conclude the necessity to extend the range of the contact tracing from -2 days to -3 days based on the day of symptom onset.

Second, it can be used for cases where the CIs of both the CTR and CTE by interval are out of range. For example, a confirmed case disclosed the interval day as one day during the first epidemiological investigation, and both the CTR and CTE values were 23.0. These cases do not fall within the CIs of the 1-day interval in Table 3, but they satisfy the CIs of 3-5 days. In this case, the EIT can predict the presence of mild COVID-19 symptoms, which will allow the patient to recognize the symptoms after two or more days, depending on the difference in the individual level of recognition. Thus, the EIT should reassess the day of symptom onset with the confirmed case through a second epidemiological investigation and then adjust the scope of the epidemiological investigation according to the information disclosed by the case. Because it allows for the additional identification of contacts, it should prevent transmission in the local community.

Finally, another example could be when a person was identified as a contact and underwent an RT-PCR test, obtaining a 30.0. The patient disclosed no COVID-19 symptoms during the epidemiological investigation. In this case, the EIT can consider the following two possibilities for deciding this confirmed case as an asymptomatic case: (1) a case where 8-10 or longer days had passed since the testing day according to the CIs in Table 3, and (2) a case immediately after the infection period has just started. To determine which of these possibilities should be applied, contacts within -2 days of the test day of the corresponding asymptomatic confirmed case should be tested. If these tests confirm new cases, it should be assumed that the infection period began 8-10 days prior, and the case was demonstrated to be a preceding infection; therefore, epidemiological investigations should consider expanding its range.

The major limitations of the present study are as follows. First, the symptom onset days identified by the epidemiological investigation were based only on the description of the confirmed cases. Therefore, the onset day of COVID-19 suspected symptoms can differ depending on the recognition level of the confirmed case; hence, it is subjectively determined [4]. Moreover, the interval days obtained from the cases in the present study might be different by a few days, as described in the first two hypothetical cases, which was also supported by the statistically significant differences in the paired t-test for CTR and CTE values for each case. To overcome this, additional studies should be considered including a large number of test cases who manifested symptoms and were confirmed as having COVID-19 on the second RT-PCR test conducted within 24 hours because of COVID-19-related symptoms during the 14-day self-isolation period after obtaining negative RT-PCR results from contact tracing should be considered [25]. Furthermore, a more valid CTR and CTE range for day 0 can be obtained with more sensitive recognition of the onset of even mild symptoms. However, to overcome the limitation of subjectivity for other interval days, specific interval days obtained from the first epidemiological investigation should be subjected to a test if they fall within the estimated CIs of CTR and CTE, as in the first two hypothetical cases, and reflected in the second epidemiological investigation. Second, the interval in the present study was limited to 0-10 days. Thus, further studies should be conducted at intervals of 11 days or longer. Third, the subjects were limited to residents of a specific region; therefore, the results may differ from those of residents of other provinces or countries. However, it may not be expected to show a marked difference because all subjects were exposed to COVID-19 for the first time, and the CT value is a measure of the SARS-CoV-2 gene. The population vaccination program was not initiated until the study subjects were selected. Finally, the period covered in this study was before the COVID-19 variant became the dominant species. In other words, as mutations such as delta and omicron become the dominant species, the direct application of the results of this study is limited. However, if the same analysis methodology used in this study is applied to future strains, the infection period can be estimated and reflected in epidemiological investigations. Future studies should attempt to estimate the CTR and CTE values more accurately through further analysis of more confirmed cases.

Acknowledgements

We deeply appreciate the staff of the infectious disease team in six community health centers in Jeju Province who performed epidemiological investigations on confirmed cases of COVID-19; the RT-PCR test staff of the Jeju Institute of Health and Environment who conducted RTPCR tests to confirm COVID-19; and the staff of Team for Infection Control of Jeju Province Office who analyzed epidemiological investigation reports and decided on the countermeasures.

References

1. Hong KH, Lee SW, Kim TS, Huh HJ, Lee J, Kim SY, et al. Guidelines for laboratory diagnosis of coronavirus disease 2019 (COVID-19) in Korea. Ann Lab Med 2020;40:351–60.
2. Lim J, Lee J. Current laboratory diagnosis of coronavirus disease 2019. Korean J Intern Med 2020;35:741–8.
3. Sethuraman N, Jeremiah SS, Ryo A. Interpreting diagnostic tests for SARS-CoV-2. JAMA 2020;323:2249–51.
4. Fang FC, Naccache SN, Greninger AL. The laboratory diagnosis of coronavirus disease 2019- frequently asked questions. Clin Infect Dis 2020;71:2996–3001.
5. Korea Diseases Control and Prevention Agency (KDCA). The Guideline of COVID-19 control and management [Internet]. Cheongju: KDCA; c2021. [cited 2021 May 1]. Available from: http://kdca.go.kr/board/board.es?mid=a20507020000&bid=0019.
6. Kim SW, Kim SM, Kim YK, Kim JY, Lee YM, Kim BO, et al. Clinical characteristics and outcomes of COVID-19 cohort patients in Daegu metropolitan city outbreak in 2020. J Korean Med Sci 2021;36e12.
7. Song GS, Lee YR, Kim S, Kim W, Choi J, Yoo D, et al. Laboratory diagnosis of coronavirus disease 19 (COVID-19) in Korea: current status, limitation, and challenges. Korean J Clin Lab Sci 2020;52:284–95.
8. Bai Y, Yao L, Wei T, Tian F, Jin DY, Chen L, et al. Presumed asymptomatic carrier transmission of COVID-19. JAMA 2020;323:1406–7.
9. Huh HJ, Hong KH, Kim TS, Song SH, Roh KH, Lee H, et al. Surveillance of coronavirus disease 2019 (COVID-19) testing in clinical laboratories in Korea. Ann Lab Med 2021;41:225–9.
10. Ji T, Liu Z, Wang G, Guo X, Khan SA, Lai C, et al. Detection of COVID-19: a review of the current literature and future perspectives. Biosens Bioelectron 2020;166:112455.
11. Corman VM, Landt O, Kaiser M, Molenkamp R, Meijer A, Chu DK, et al. Detection of 2019 novel coronavirus (2019-nCoV) by real-time RT-PCR. Euro Surveill 2020;25:2000045.
12. Sung H, Roh KH, Hong KH, Seong MW, Ryoo N, Kim HS, et al. COVID-19 molecular testing in Korea: practical essentials and answers from experts based on experiences of emergency use authorization assays. Ann Lab Med 2020;40:439–47.
13. Drew RJ, O’Donnell S, LeBlanc D, McMahon M, Natin D. The importance of cycle threshold values in interpreting molecular tests for SARS-CoV-2. Diagn Microbiol Infect Dis 2020;98:115130.
14. Tom MR, Mina MJ. To interpret the SARS-CoV-2 test, consider the cycle threshold value. Clin Infect Dis 2020;71:2252–4.
15. Schnuriger A, Perrier M, Marinho V, Michel Y, Saloum K, Boukli N, et al. Caution in interpretation of SARS-CoV-2 quantification based on RT-PCR cycle threshold value. Diagn Microbiol Infect Dis 2021;100:115366.
16. Cevik M, Tate M, Lloyd O, Maraolo AE, Schafers J, Ho A. SARSCoV-2, SARS-CoV, and MERS-CoV viral load dynamics, duration of viral shedding, and infectiousness: a systematic review and meta-analysis. Lancet Microbe 2021;2:e13–22.
17. Engelmann I, Alidjinou EK, Ogiez J, Pagneux Q, Miloudi S, Benhalima I, et al. Preanalytical issues and cycle threshold values in SARS-CoV-2 real-time RT-PCR testing: should test results include these? ACS Omega 2021;6:6528–36.
18. Cleveland WS, Devlin SJ. Locally weighted regression: an approach to regression analysis by local fitting. J Am Stat Assoc 1988;83:596–610.
19. Ostertagová E. Modelling using polynomial regression. Procedia Engineering 2012;48:500–6.
20. Foster G. Interpreting and visualizing regression models using Stata, by Michael N. Mitchell (Stata Press, College Station, Texas, 2012), pp. xxix + 558. Econ Rec 2013;89:132–4.
21. Ade C, Pum J, Abele I, Raggub L, Bockmühl D, Zöllner B. Analysis of cycle threshold values in SARS-CoV-2-PCR in a long-term study. J Clin Virol 2021;138:104791.
22. Kang J, Jang YY, Kim J, Han SH, Lee KR, Kim M, et al. South Korea’s responses to stop the COVID-19 pandemic. Am J Infect Control 2020;48:1080–6.
23. Cheng MP, Papenburg J, Desjardins M, Kanjilal S, Quach C, Libman M, et al. Diagnostic testing for severe acute respiratory syndrome-related coronavirus 2: a narrative review. Ann Intern Med 2020;172:726–34.
24. Binnicker MJ. Emergence of a novel coronavirus disease (COVID-19) and the importance of diagnostic testing: why partnership between clinical laboratories, public health agencies, and industry is essential to control the outbreak. Clin Chem 2020;66:664–6.
25. Jung J, Jang H, Kim HK, Kim J, Kim A, Ko KP. The importance of mandatory COVID-19 diagnostic testing prior to release from quarantine. J Korean Med Sci 2020;35e314.

Article information Continued

Figure 1.

The locally weighted regression smoothing (LOWESS) curves of critical values (yhatlowess) of the RdRp gene (CRT) (A) and critical values (yhat) of the E gene (CRE) (B) by interval showing a U-shape pattern.

Figure 2.

Predicted confidence intervals of critical values of the RdRp gene (CRT) (A) and E gene (CRE) (B) by interval are displayed.

Table 1.

Distribution of the cycle threshold in the RdRp gene (CTR) and E gene (CTE)

Interval (days) Number CTR
CTE
Mean SD Median Min Max P-value Mean SD Median Min Max P-value
All 331 23.53 5.59 22.4 11.01 37.37 23.01 5.46 21.98 10.65 37.84
Sex
Men 150 23.87 5.47 23.05 15.36 36.79 0.306* 23.45 5.34 22.69 15.11 37.84 0.177*
Women 181 23.24 5.69 21.96 11.01 37.37 0.246 22.64 5.55 21.42 10.65 36.19 0.148
Age
≤19 years 43 25.40 5.96 24.91 15.61 37.37 0.072 24.82 5.60 23.69 14.49 35.70 0.066
20-39 years 96 23.60 5.18 23.35 15.36 36.63 0.087§ 22.97 5.11 22.49 14.26 35.86 0.063§
40-59 years 129 23.34 5.80 21.17 11.01 36.79 22.95 5.73 21.33 10.65 37.84
≥60 years 63 22.53 5.30 21.94 14.63 37.10 21.94 5.11 21.43 14.46 36.19
Interval
0 55 24.50 6.12 22.94 15.61 37.37 <0.001 23.83 5.90 21.79 14.49 36.14 <0.001
1 76 21.73 4.82 20.04 14.63 36.79 <0.001§ 21.06 4.67 20.11 14.46 36.28 <0.001§
2 56 21.86 5.19 20.44 14.66 36.35 21.41 5.17 20.22 14.59 36.19
3 47 23.37 4.99 23.24 16.19 34.44 23.02 4.79 22.50 15.27 34.45
4 28 22.08 3.58 21.53 16.89 34.40 21.84 3.31 22.10 17.08 33.96
5 15 24.29 5.04 23.80 16.28 35.60 23.83 4.89 23.30 15.94 34.10
6 16 26.36 6.00 26.59 15.80 35.48 25.72 5.99 25.66 14.26 35.67
7 21 26.08 6.77 26.84 11.01 37.10 25.51 6.51 25.85 10.65 35.86
8 5 31.18 4.27 31.04 26.2 35.73 30.49 4.57 30.51 25.72 35.94
9 7 29.05 3.25 29.48 23.17 33.60 28.79 2.69 29.09 23.9 32.46
10 5 30.96 3.98 31.22 27.20 36.71 30.58 4.75 30.18 25.66 37.84

SD: standard deviation.

*

t-test value,

Mann-Whitney test value,

analysis of variance value,

§

Kruskal-Wallis test value.

Table 2.

Outputs of polynomial regression analysis applying the quadratic term of intervals

Variable Cycle threshold of RdRp gene Cycle threshold of E gene
R-squared value 0.119 (0.000) 0.126 (0.000)
Coefficient of interval -0.814 (0.023) -0.701 (0.044)
Coefficient of interval2 0.176 (0.000) 0.166 (0.000)
Constant 23.35 (0.000) 22.67 (0.000)

Values are presented as number (P-value).

Table 3.

Predicted confidence intervals of cycle threshold of the RdRp gene (CTR) and the E gene (CTE) by intervals based on the data found in Table 2

Intervals (day) 95% CI of CTR
95% CI of CTE
Min Max Min Max
0 23.11 25.89 22.48 25.18
1 20.55 22.91 19.91 22.21
2 20.48 23.23 20.08 22.75
3 21.87 24.87 21.56 24.48
4 20.13 24.02 19.95 23.72
5 21.4 26.95 21.25 26.42
6 23.78 28.92 23.21 28.22
7 23.83 28.32 23.33 27.70
8 26.57 35.78 26.02 34.96
9 25.16 32.94 25.01 32.57
10 26.36 35.56 26.11 35.06

CI: confidence interval.