Acessibilidade / Reportar erro

Associations of SARS-CoV-2 cycle threshold values with age, gender, sample collection setting, and pandemic period

ABSTRACT

Cycle threshold (Ct) values in COVID-19 reverse-transcription polymerase chain reaction (RT-PCR) tests estimate the viral load in biological samples. Studies have investigated variables associated with SARS-CoV-2 viral load, aiming to identify factors associated with higher transmissibility. Using the results from tests performed between May/2020-July/2022 obtained from the database of a referent hospital in Sao Paulo, Brazil, we investigated associations between Ct values and patient’s age, gender, sample collection setting and pandemic period according to the predominant SARS-CoV-2 variant locally. We also examined variations in Ct values, COVID-19 incidence, mortality, and vaccination coverage over time. The study sample included 42,741 tests. Gender was not significantly associated with Ct values. Age, sample collection setting and the pandemic period were significantly associated with Ct values even after adjustment to the multivariable model. Results showed lower Ct values in older groups, during the Gamma and Delta periods, and in samples collected in emergency units; and higher Ct values in children under 10 years old, home-based tests, during the Omicron period. We found evidence of a linear trend in the association between age and Ct values, with Ct values decreasing as age increases. We found no clear temporal associations between Ct values and local indicators of COVID-19 incidence, mortality, or vaccination between February/2020-November/2022. Our findings suggest that SARS-CoV-2 Ct values, a proxy for viral load and transmissibility, can be influenced by demographic and epidemiological variables.

SARS-CoV-2; COVID-19; Reverse transcriptase polymerase chain reaction; SARS-CoV-2 variants; Viral load; COVID-19 vaccines; Demographic factors

INTRODUCTION

Reverse Transcription Polymerase Chain Reaction (RT-PCR) nucleic acid amplification is the gold standard diagnostic test to document the presence of SARS-CoV-2 in respiratory samples11. Phillips MC, Quintero D, Wald-Dickler N, Holtom P, Butler-Wu SM. SARS-CoV-2 cycle threshold (Ct) values predict future COVID-19 cases. J Clin Virol. 2022;150-151:105153.. RT-PCR has been extensively used from the onset of the COVID-19 pandemic to detect SARS-CoV-2 in both symptomatic and asymptomatic individuals. However, most RT-PCR assays provide not only qualitative but also quantitative results, estimating the viral load in a given sample as a function of the cycle threshold (Ct) value22. Tom MR, Mina MJ. To interpret the SARS-CoV-2 test, consider the cycle threshold value. Clin Infect Dis. 2020;71:2252-4.

3. Rhoads D, Peaper DR, She RC, Nolte FS, Wojewoda CM, Anderson NW, et al. College of American Pathologists (CAP) Microbiology Committee perspective: caution must be used in interpreting the cycle threshold (Ct) value. Clin Infect Dis. 2021;72:e685-6.
-44. Jones TC, Mühlemann B, Veith T, Biele G, Zuchowski M, Hofmann J, et al. Estimating infectiousness throughout SARS-CoV-2 infection course. Science. 2021;373:180.. Ct measures the number of PCR cycles needed for the fluorescence level in a sample, which increases proportionally to the amount of amplified RNA, to surpass the cutoff for positivity; thus, the higher the Ct, the lower the expected RNA amount in a sample55. 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..

Although Ct is not a direct measure of the amount of viable and infecting virus copies in a sample, prior studies indicate a relationship between lower Ct values and higher probability of viral isolation in cell cultures66. La Scola B, Le Bideau M, Andreani J, Hoang VT, Grimaldier C, Colson P, et al. Viral RNA load as determined by cell culture as a management tool for discharge of SARS-CoV-2 patients from infectious disease wards. Eur J Clin Microbiol Infect Dis. 2020;39:1059-61.,77. Bullard J, Dust K, Funk D, Strong JE, Alexander D, Garnett L, et al. Predicting infectious severe acute respiratory syndrome coronavirus 2 from diagnostic samples. Clin Infect Dis. 2020;71:2663-6.. Since identifying and isolating potential SARS-CoV-2 spreaders is an important strategy to control virus transmission, studies have adopted Ct values as proxies for viral loads to identify features that confer higher or lower infectivity to a patient by analyzing associations between Ct values and clinical or demographic characteristics11. Phillips MC, Quintero D, Wald-Dickler N, Holtom P, Butler-Wu SM. SARS-CoV-2 cycle threshold (Ct) values predict future COVID-19 cases. J Clin Virol. 2022;150-151:105153..

To better understand the dynamics and determinants of SARS-CoV-2 Ct values, studies have investigated its association with demographic factors such as age, gender, and race. Thus far, results concerning the effects of age on Ct values have been inconsistent44. Jones TC, Mühlemann B, Veith T, Biele G, Zuchowski M, Hofmann J, et al. Estimating infectiousness throughout SARS-CoV-2 infection course. Science. 2021;373:180.,88. Penney J, Jung A, Koethe B, Doron S. Cycle threshold values and SARS-CoV-2: relationship to demographic characteristics and disease severity. J Med Virol. 2022;94:3978-81.

9. Hurst JH, Heston SM, Chambers HN, Cunningham HM, Price MJ, Suarez L, et al. Severe acute respiratory syndrome coronavirus 2 infections among children in the biospecimens from respiratory virus-exposed kids (BRAVE Kids) study. Clin Infect Dis. 2021;73:e2875-82.

10. Baggio S, L’Huillier AG, Yerly S, Bellon M, Wagner N, Rohr M, et al. Severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) viral load in the upper respiratory tract of children and adults with early acute COVID-19. Clin Infect Dis. 2021;73:148-50.
-1111. 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.. Research exploring associations between Ct values and clinical presentation suggests that higher viral loads may be associated with an increased risk of intubation, overall severity, and mortality; yet, other studies have shown different results88. Penney J, Jung A, Koethe B, Doron S. Cycle threshold values and SARS-CoV-2: relationship to demographic characteristics and disease severity. J Med Virol. 2022;94:3978-81.,1212. Magleby R, Westblade LF, Trzebucki A, Simon MS, Rajan M, Park J, et al. Impact of severe acute respiratory syndrome coronavirus 2 viral load on risk of intubation and mortality among hospitalized patients with coronavirus disease 2019. Clin Infect Dis. 2021;73:e4197-205.

13. Huang JT, Ran RX, Lv ZH, Feng LN, Ran CY, Tong YQ, et al. Chronological changes of viral shedding in adult inpatients with COVID-19 in Wuhan, China. Clin Infect Dis. 2020;71:2158-66.
-1414. Yu X, Sun S, Shi Y, Wang H, Zhao R, Sheng J. SARS-CoV-2 viral load in sputum correlates with risk of COVID-19 progression. Crit Care. 2020;24:170..

While the SARS-CoV-2 Ct dynamics in asymptomatic or pre-symptomatic patients are poorly understood, Ct values typically vary over the course of symptomatic COVID-19 in each individual1515. 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.. Both longitudinal and cross-sectional studies have observed lower values in the first days after symptoms onset, followed by a progressive increase in subsequent days88. Penney J, Jung A, Koethe B, Doron S. Cycle threshold values and SARS-CoV-2: relationship to demographic characteristics and disease severity. J Med Virol. 2022;94:3978-81.,99. Hurst JH, Heston SM, Chambers HN, Cunningham HM, Price MJ, Suarez L, et al. Severe acute respiratory syndrome coronavirus 2 infections among children in the biospecimens from respiratory virus-exposed kids (BRAVE Kids) study. Clin Infect Dis. 2021;73:e2875-82.,1515. 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.

16. Jang S, Rhee JY, Wi YM, Jung BK. Viral kinetics of SARS-CoV-2 over the preclinical, clinical, and postclinical period. Int J Infect Dis. 2021;102:561-5.

17. Kissler SM, Fauver JR, Mack C, Tai CG, Breban MI, Watkins AE, et al. Viral dynamics of SARS-CoV-2 variants in vaccinated and unvaccinated persons. N Engl J Med. 2021;385:2489-91.
-1818. Singanayagam A, Patel M, Charlett A, Lopez Bernal J, Saliba V, Ellis J, et al. Duration of infectiousness and correlation with RT-PCR cycle threshold values in cases of COVID-19, England, January to May 2020. Euro Surveill. 2020;25:2001483.. In studies that do not involve longitudinal follow-up of participants, the sample collection setting can be used as a proxy of clinical presentation, sometimes enabling assumptions about the presence of symptoms and their duration. Samples collected for screening purposes, for example, such as those obtained from athletes and company employees, are presumably more often collected from asymptomatic or pre-symptomatic individuals. Conversely, patients tested in hospitals or emergency services are more likely symptomatic and in the initial phases of disease onset, while those tested in outpatient laboratory facilities presumably had less intense symptoms. On the other hand, patients tested at home are more likely to undergo sample collection several days after symptoms onset due to delays in test scheduling. Hence, the location and sample collection setting could influence Ct values.

Finally, some studies suggest that Ct values can show population-level variations over different periods of the pandemic, with lower mean values coinciding with periods of rising spread, such as those seen during the upsurge of novel SARS-CoV-2 variants11. Phillips MC, Quintero D, Wald-Dickler N, Holtom P, Butler-Wu SM. SARS-CoV-2 cycle threshold (Ct) values predict future COVID-19 cases. J Clin Virol. 2022;150-151:105153.,1515. 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.. Multiple variants have been identified since the beginning of the pandemic, with different pathogenicity and infectivity profiles1919. Volz E, Mishra S, Chand M, Barrett JC, Johnson R, Geidelberg L, et al. Assessing transmissibility of SARS-CoV-2 lineage B.1.1.7 in England. Nature. 2021;593:266-9.,2020. Volz E, Hill V, McCrone JT, Price A, Jorgensen D, O’Toole A, et al. Evaluating the effects of SARS-CoV-2 spike mutation D614G on transmissibility and pathogenicity. Cell. 2021;184:64-75.e11.. In Brazil, predominant lineages included B1.195, B.1.1.128 and B1.1.332121. Voloch CM, Francisco Jr RS, Almeida LG, Cardoso CC, Brustolini OJ, Gerber AL, et al. Genomic characterization of a novel SARS-CoV-2 lineage from Rio de Janeiro, Brazil. J Virol. 2021;95:e00119-21.

22. Resende PC, Delatorre E, Gräf T, Mir D, Motta FC, Appolinario LR, et al. Evolutionary dynamics and dissemination pattern of the SARS-CoV-2 lineage B.1.1.33 during the early pandemic phase in Brazil. Front Microbiol. 2020;11:615280.
-2323. Wolf JM, Kipper D, Borges GR, Streck AF, Lunge VR. Temporal spread and evolution of SARS-CoV-2 in the second pandemic wave in Brazil. J Med Virol. 2022;94:926-36.. P.2 variant may have first appeared in the country in July 2020, and P.1 surge probably occurred in December of the same year. Both P.1 and P.2 seem to be derived from B.1.1.28. P.2 spread throughout the country between September 2020 and January 2021, whereas P.1 saw even faster dissemination, covering the entire Brazilian territory between January and February 20212323. Wolf JM, Kipper D, Borges GR, Streck AF, Lunge VR. Temporal spread and evolution of SARS-CoV-2 in the second pandemic wave in Brazil. J Med Virol. 2022;94:926-36.. Gamma, as P.1 was latter designated, became the most frequently detected variant among all SARS-CoV-2 strains in RT-PCR tests in Sao Paulo State and across Brazil by February 2021; its spread was accompanied by an important surge in the number of COVID-19 cases, corresponding to the “second wave” of the pandemic in the country2323. Wolf JM, Kipper D, Borges GR, Streck AF, Lunge VR. Temporal spread and evolution of SARS-CoV-2 in the second pandemic wave in Brazil. J Med Virol. 2022;94:926-36.

24. Sarsomics. SARS-OMICS: atlas multi-ômico da COVID-19 brasileira. [cited 2023 Jul 17]. Available from: https://sarsomics.com
https://sarsomics.com...
-2525. Fundação Oswaldo Cruz. Dashboard genomic network. [cited 2023 Jul 17]. Available from: https://www.genomahcov.fiocruz.br/dashboard-en/
https://www.genomahcov.fiocruz.br/dashbo...
. By August 2021, when reported deaths and cases were declining, Delta had replaced Gamma as the most prominent variant in Brazil2626. Giovanetti M, Fonseca V, Wilkinson E, Tegally H, San EJ, Althaus CL, et al. Replacement of the Gamma by the Delta variant in Brazil: impact of lineage displacement on the ongoing pandemic. Virus Evol. 2022;8:veac024.. Finally, the Omicron variant, first reported in South Africa, surpassed the detection rate of the Delta variant in Sao Paulo around November 2021, and was associated with a fast increase in the number of COVID-19 cases2424. Sarsomics. SARS-OMICS: atlas multi-ômico da COVID-19 brasileira. [cited 2023 Jul 17]. Available from: https://sarsomics.com
https://sarsomics.com...
,2727. World Health Organization. Classification of Omicron (B.1.1.529): SARS-CoV-2 variant of concern. [cited 2023 Jul 17]. Available from: https://www.who.int/news/item/26-11-2021-classification-of-omicron- (b.1.1.529)-sars-cov-2-variant-of-concern
https://www.who.int/news/item/26-11-2021...
,2828. World Health Organization. WHO Coronavirus (COVID-19) dashboard. [cited 2023 Jul 17]. Available from: https://covid19.who.int
https://covid19.who.int...
.

In this study, we used information collected from the database of a reference laboratory in Sao Paulo, Brazil (23° 32’ 56”S, 46° 38’ 20”W) to investigate associations between Ct values and patient’s age, gender, sample collection setting and pandemic period, categorized according to the predominant SARS-CoV-2 variant locally. We also explored variations in Ct values, incidence of COVID-19 cases, mortality, and vaccination coverage over time.

MATERIALS AND METHODS

This cross-sectional study included data from patients with SARS-CoV-2 infection confirmed by a positive RT-PCR test performed in nasal and/or oral cavity swabs at the Hospital Israelita Albert Einstein laboratory. This tertiary hospital has a reference laboratory that supports 17 inpatient and outpatient facilities located in 11 different municipalities in the Sao Paulo State with diagnostic resources. Participants included in the study underwent RT-PCR testing as requested by their medical providers, without additional tests for the purpose of this research. Qualified personnel processed the samples according to institution’s standards and the test manufacturer’s instructions, storing the test results, including Ct values, in a centralized database along with demographic information. We initially selected all patients who underwent SARS-CoV-2 testing from May 14, 2020 to July 26, 2022. From this initial sample, we used Python (version 3.9.7, Python Software Foundation, Delaware, USA) and libraries pandas 1.3.3 and NumPy 1.20.3 to identify all patients whose test had been done using the Cobas-6800 System and Cobas® SARS-CoV-2 test (Roche Diagnostics, Basel, Switzerland). For each unique patient, we maintained only the first test in the analysis to eliminate test results that could correspond to recurrent testing, and to focus on early disease onset.

Ct measurements are influenced by factors such as the source and type of specimen, collection technique, and technical characteristics of the assay33. Rhoads D, Peaper DR, She RC, Nolte FS, Wojewoda CM, Anderson NW, et al. College of American Pathologists (CAP) Microbiology Committee perspective: caution must be used in interpreting the cycle threshold (Ct) value. Clin Infect Dis. 2021;72:e685-6.. Thus, we opted for restricting our analysis to Ct values generated with a single RT-PCR equipment to avoid potential biases generated by different manufacturers and models. We selected the Cobas® SARS-CoV-2 assay because it had the highest and most stable volume of results throughout the study period. The assay generates two Ct values for each sample: one for gene E and another for ORF1a/b. Only the results for gene ORF1a/b were included in this study.

We compared Ct values across different sample collection settings categorized as those made primarily with screening purposes (athletes and company employees); and those collected in laboratory units; emergency units; hospitals; outpatient clinics; and home-based tests. Our analysis included these categories to explore differences in Ct values associated with the collection setting, which could indirectly indicate the presence/severity of symptoms and time since their onset if applicable.

Participant characteristics are presented as counts and percentages for categorical variables, and medians and interquartile ranges (IQR) for numerical variables. Age is presented both as a numeric variable (years) and as a categorical variable (decade). We used data from the SARS-OMICs database, which includes test results from Hospital Israelita Albert Einstein, to define four pandemic periods according to the predominant SARS-CoV-2 variant: October 2020 to December 2020 (pre-Gamma); January to July 2021 (Gamma); August 2021 to October 2021 (Delta); and November 2021 to September 2022 (Omicron)2424. Sarsomics. SARS-OMICS: atlas multi-ômico da COVID-19 brasileira. [cited 2023 Jul 17]. Available from: https://sarsomics.com
https://sarsomics.com...
.

Associations between Ct values and patient’s age, gender, sample collection setting and pandemic period were calculated using descriptive analysis, as well as univariable and multivariable linear regression models with robust standard errors and a 0.05 significance level. In the univariable models, we present coefficients (effect estimates) for each category of the independent variables compared to a referent category, along with its 95% confidence intervals. The multivariable model included all variables to explore the association of each independent variable adjusted for the remaining covariates studied. We used a linear trend test to assess if Ct values varied linearly across age groups.

Finally, we evaluated the relation between Ct values, COVID-19 incidence, mortality, and vaccination coverage over time between February 2020 and November 2022, using graphical methods and health indicators made available online by the state health department2929. São Paulo. Governo do Estado. Fundação Sistema Estadual de Análise de Dados. SP contra o novo coronavírus: boletim completo. [cited 2023 Jul 17]. Available from: https://www.seade.gov.br/coronavirus/
https://www.seade.gov.br/coronavirus/...
,3030. São Paulo. Governo do Estado. #Vacina Já. [cited 2023 Jul 17]. Available from: https://www.vacinaja.sp.gov.br/vacinometro/
https://www.vacinaja.sp.gov.br/vacinomet...
. COVID-19 incidence and mortality in the Sao Paulo State (total population of 46,649,132 residents) were presented as frequencies/counts per month3131. Pujadas E, Chaudhry F, McBride R, Richter F, Zhao S, Wajnberg A, et al. SARS-CoV-2 viral load predicts COVID-19 mortality. Lancet Respir Med. 2020;8:e70.. Vaccination coverage in the Sao Paulo extended metropolitan area was described as cumulative percentage per month. The state’s extended metropolitan region includes ten municipalities: Sao Paulo, Cotia, Franco da Rocha, Guarulhos, Mogi das Cruzes, Osasco, Santana do Parnaiba, Sao Bernardo do Campo, Suzano, and Taboao da Serra, with an estimated population of 16,970,273 inhabitants3131. Pujadas E, Chaudhry F, McBride R, Richter F, Zhao S, Wajnberg A, et al. SARS-CoV-2 viral load predicts COVID-19 mortality. Lancet Respir Med. 2020;8:e70..

Statistical package Stata (version 15.1, StataCorp LLC, Texas, USA) and Microsoft® Excel (version 16.68, Microsoft Corporation, Washington, USA) were used for all analyses.

The institution’s ethics committee reviewed and approved the study with exemption of informed consent (CAAE 34882620.7.0000.0071, SGPP 4280-20). We kept all identifiable patient information confidential throughout the study.

RESULTS

A total of 1,088,682 RT-PCR tests were collected between May 2020 and July 2022 using the Cobas-6800 System and Cobas® SARS-CoV-2 test and registered in the institution’s database. Of these, we excluded 933,671 tests with a negative result and/or with a Ct value higher than 38 (n=155,011); 62,764 positive results without available corresponding Ct values (n=92,247); 163 tests collected in settings outside the study categories (n=92,084); 45,941 results for gene E (n=46,143); 2,202 test repetitions for the same patient (n=43,941); one test with inconsistent Ct value (n=43,940); four tests from patients younger than one-year-old inconsistently categorized as athletes or employees (n=43,946); and 1,194 results from small municipalities where vaccination and testing dynamics could differ from those in the Sao Paulo extended metropolitan region (n=42,741).

Hence, the final sample included 42,741 tests. Table 1 summarizes the demographic characteristics of study participants, sample collection settings, and pandemic periods. Most participants were female (54.7%), aged 20 to 59 years old (66.5%), and tested during the Omicron (49.3%) or pre-Gamma (39.3%) variant predominance. Laboratories were the most frequent sample collection setting (24.4%), followed by athletes and employees (20.7%), outpatient clinics (18.5%) and emergency units (17.8%).

The 50-59 and 60-69 age groups (21.5), emergency units (20.8), athletes and employees (21.4) and Delta variant period (20.9) showed the lowest mean Ct values. In turn, age group 0-9 years (23.6), home-based tests (24.0) and Omicron period (22.5) presented the highest mean Ct values (Table 1). Figures 1A, 1B and 1C show the distributions of Ct values by age group, pandemic period and sample collection setting, respectively.

Table 1
Characteristics of study participants.

Figure 1
(A) Cycle threshold values distributed by age groups. Squares indicate mean values; (B) Cycle threshold values distributed by pandemic period. Squares indicate mean values; (C) Cycle threshold values distributed by sample collection setting. Squares indicate mean values.

Table 2 summarizes the results from the univariable and multivariable linear regression models. We found no statistically significant differences between Ct values of male and female participants. Both univariable and multivariable analyses showed a consistent, statistically significant association between age and Ct values, with decreasing Ct values associated with older age. Accordingly, the multivariable model indicated a statistically significant linear trend in the association between age categories and Ct values (p<0.0001). Sample collection setting was also significantly associated with Ct values, with the lowest Ct values observed in samples collected in emergency units. Samples from home-based tests had, on average, a 3.10-unit higher Ct value (95% confidence interval 2.90 to 3.30, p<0.001) adjusted for age, gender, and pandemic period.

Table 2
Univariable and multivariable linear regression models for the associations between age group, gender, sample collection setting, and pandemic period with SARS-CoV-2 RT-PCR cycle thresholds.

Finally, the pandemic period showed statistically significant association with Ct values, with higher Ct values found in samples from the Omicron predominance period (Table 2).

Figure 2 presents a graphical evaluation of Ct results over time along with indicators of COVID-19 incidence (absolute number of cases per month), mortality (absolute number of deaths per month), and vaccination coverage (cumulative percentage per month). We observed no clear changes in Ct values according to these factors, neither a consistent trend over time; however, we observed a noticeable decrease in the mortality-to-cases ratio following the increase in vaccination coverage.

Figure 2
Cycle threshold values, COVID-19 incidence and mortality, and vaccination coverage over time.

DISCUSSION

This study investigated associations between age, gender, sample collection setting, and pandemic period defined by the predominant SARS-CoV-2 variant locally and Ct values, used as a proxy of viral loads in respiratory secretions. Age, sample collection setting and pandemic period were significantly associated with Ct values even after adjustment in a multivariable model. We observed lower Ct values in samples from participants tested for screening purposes (athletes and company employees) and in emergency units, possibly because these individuals are likely pre-symptomatic or seeking care early in the disease onset, as seen in other studies1515. 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.

16. Jang S, Rhee JY, Wi YM, Jung BK. Viral kinetics of SARS-CoV-2 over the preclinical, clinical, and postclinical period. Int J Infect Dis. 2021;102:561-5.

17. Kissler SM, Fauver JR, Mack C, Tai CG, Breban MI, Watkins AE, et al. Viral dynamics of SARS-CoV-2 variants in vaccinated and unvaccinated persons. N Engl J Med. 2021;385:2489-91.
-1818. Singanayagam A, Patel M, Charlett A, Lopez Bernal J, Saliba V, Ellis J, et al. Duration of infectiousness and correlation with RT-PCR cycle threshold values in cases of COVID-19, England, January to May 2020. Euro Surveill. 2020;25:2001483.. Higher Ct values among participants tested at home could also be explained by this phenomenon, given that these tests are usually scheduled after symptoms onset, leading to delays compared with other testing modalities. Different disease presentation patterns may play a role in these results, as individuals undergoing testing in emergency units are expected to be more symptomatic than patients in other settings.

Previous studies have shown conflicting results concerning the correlation between SARS-CoV-2 viral load and disease severity. While a lower Ct value could reflect higher viral replication and direct tissue damage, tests done at later timepoints in the disease course—in patients with escalating symptoms and exacerbated inflammatory response—could also show higher Ct values88. Penney J, Jung A, Koethe B, Doron S. Cycle threshold values and SARS-CoV-2: relationship to demographic characteristics and disease severity. J Med Virol. 2022;94:3978-81.,1212. Magleby R, Westblade LF, Trzebucki A, Simon MS, Rajan M, Park J, et al. Impact of severe acute respiratory syndrome coronavirus 2 viral load on risk of intubation and mortality among hospitalized patients with coronavirus disease 2019. Clin Infect Dis. 2021;73:e4197-205.,3232. Salvatore PP, Dawson P, Wadhwa A, Rabold EM, Buono S, Dietrich EA, et al. Epidemiological correlates of polymerase chain reaction cycle threshold values in the detection of severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2). Clin Infect Dis. 2021;72:e761-7.. Shah et al.’s systematic review and meta-analysis showed that hospitalized patients with lower Ct values had a higher mortality risk and disease severity3333. Shah VP, Farah WH, Hill JC, Hassett LC, Binnicker MJ, Yao JD, et al. Association between SARS-CoV-2 cycle threshold values and clinical outcomes in patients with COVID-19: a systematic review and meta-analysis. Open Forum Infect Dis. 2021;8:ofab453.. Another study found an association between higher Ct values and increased odds of severe presentation88. Penney J, Jung A, Koethe B, Doron S. Cycle threshold values and SARS-CoV-2: relationship to demographic characteristics and disease severity. J Med Virol. 2022;94:3978-81.. Finally, a third study failed to find statistically significant associations between Ct values and poor outcomes3434. Martínez MJ, Basile L, Sisó-Almirall A, Cristino V, Cuesta G, Hurtado JC, et al. Lack of prognostic value of SARS-CoV2 RT-PCR cycle threshold in the community. Infect Dis Ther. 2022;11:587-93..

As observed in other publications, we found no consistent association between gender and Ct values88. Penney J, Jung A, Koethe B, Doron S. Cycle threshold values and SARS-CoV-2: relationship to demographic characteristics and disease severity. J Med Virol. 2022;94:3978-81.,3535. Heald-Sargent T, Muller WJ, Zheng X, Rippe J, Patel AB, Kociolek LK. Age-related differences in nasopharyngeal severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) levels in patients with mild to moderate coronavirus disease 2019 (COVID-19). JAMA Pediatr. 2020;174:902-3.; however, age and Ct values showed a statistically significant association, suggesting lower SARS-CoV-2 viral loads in younger age groups. Observed in both the univariable and multivariable analyses, this finding was reinforced by a significant linear trend test. The effect of age on Ct values is not well established in the literature88. Penney J, Jung A, Koethe B, Doron S. Cycle threshold values and SARS-CoV-2: relationship to demographic characteristics and disease severity. J Med Virol. 2022;94:3978-81.,3636. Brasil. Ministério da Saúde. Agência Nacional de Vigilância Sanitária. Resolução - RDC Nº 377, de 28 de abril de 2020. Autoriza, em caráter temporário e excepcional, a utilização de “testes rápidos” (ensaios imunocromatográficos) para a COVID-19 em farmácias, suspende os efeitos do § 2º do art. 69 e do art. 70 da Resolução de Diretoria Colegiada - RDC nº 44, de 17 de agosto de 2009. Diário Oficial da União, Brasília, 29 abr. 2020, Seção 1:56. [cited 2023 Jul 17]. Available from: https://www.in.gov.br/en/web/dou/-/resolucao-rdc-n-377-de-28-de-abril-de-2020-254429215
https://www.in.gov.br/en/web/dou/-/resol...
. Viral loads analysis of 260 patients younger than 21 years old found similar Ct values across age groups, regardless of symptoms onset99. Hurst JH, Heston SM, Chambers HN, Cunningham HM, Price MJ, Suarez L, et al. Severe acute respiratory syndrome coronavirus 2 infections among children in the biospecimens from respiratory virus-exposed kids (BRAVE Kids) study. Clin Infect Dis. 2021;73:e2875-82.. Another study with 405 participants, all within the first five days after symptoms onset, found no associations between age and Ct values1010. Baggio S, L’Huillier AG, Yerly S, Bellon M, Wagner N, Rohr M, et al. Severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) viral load in the upper respiratory tract of children and adults with early acute COVID-19. Clin Infect Dis. 2021;73:148-50.. Conversely, a German research including 3,303 COVID-19 cases found statistically significant differences between some age groups, albeit limited by the use of different assays, small number of pediatric patients, and absence of clinical data44. Jones TC, Mühlemann B, Veith T, Biele G, Zuchowski M, Hofmann J, et al. Estimating infectiousness throughout SARS-CoV-2 infection course. Science. 2021;373:180.. A larger study including 21,831 participants found lower Ct values in young adults compared with other age groups, but these effects lost significance after adjustment for symptoms and number of positive genes1111. 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.. Associations between age and Ct values could be related to behavioral factors, such as health care-seeking patterns, or age-associated differences in immunological and physiological characteristics, including clinical presentation. Understanding the mechanisms underlying such effects would require further studies based on more detailed clinical data.

Prior research suggests that a decrease in average Ct values in a population may occur during or preceding an increase in the incidence of COVID-19 cases, a pattern that could be used to predict changes in disease incidence over time11. Phillips MC, Quintero D, Wald-Dickler N, Holtom P, Butler-Wu SM. SARS-CoV-2 cycle threshold (Ct) values predict future COVID-19 cases. J Clin Virol. 2022;150-151:105153.,1111. 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.,1515. 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.. Such variation could result from the introduction and dissemination of a new variant for which most individuals have limited immunity, thus presenting with higher viral load following infection. Moreover, lower Ct values in a population are presumably associated with higher transmissibility and a rise in disease incidence. Evidence for this hypothesis comes from a study in which Ct values were inversely correlated with outbreak indicators including daily changes in COVID-19 percent positivity, transmission rates, and hospitalizations3737. Tso CF, Garikipati A, Green-Saxena A, Mao Q, 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.. Our analysis found no clear association between average Ct values and the emergence of new variants. By selecting only the first RT-PCR test for patients tested more than once, we excluded individuals at later phases of the disease course whose augmented prevalence in populations evaluated by cross-sectional sampling during a waning epidemic could potentially increase average Ct values1515. 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..

The changes in average Ct values over time observed could be partially related to SARS-CoV-2 variant local predominance, vaccination coverage, as well as behavioral and testing dynamics. For instance, point-of-care rapid tests became increasingly popular in Brazil over time and their widespread use may have resulted in those with early and mild symptoms being tested at home and seeking medical care only when presenting severe symptoms3838. Instituto Brasileiro de Geografia e Estatística. Brasil/São Paulo: panorama. [cited 2023 Jul 17]. Available from: https://cidades.ibge.gov.br/brasil/sp/panorama
https://cidades.ibge.gov.br/brasil/sp/pa...
,3939. Brasil. Ministério da Saúde. Agência Nacional de Vigilância Sanitária. Resolução - RDC Nº 595, de 28 de janeiro de 2022. Dispõe sobre os requisitos e procedimentos para a solicitação de registro, distribuição, comercialização e utilização de dispositivos médicos para diagnóstico in vitro como autoteste para detecção de antígeno do SARS-CoV-2, em consonância ao Plano Nacional de Expansão da Testagem para Covid-19 (PNE-Teste), e dá outras providências. Diário Oficial da União, Brasília, 28 jan. 2022, Seção 1:Extra A. [cited 2023 Jul 17]. Available from: https://www.in.gov.br/en/web/dou/-/resolucao-rdc-n-595-de-28-de-janeiro-de-2022-376825970
https://www.in.gov.br/en/web/dou/-/resol...
. We hypothesized that a rise in vaccination coverage could lead individuals to being less fearful of severe disease complications and therefore being less prone to seek medical care early in the disease course, with both behaviors leading to higher average Ct values over time. Nevertheless, we found no clear evidence of such trends.

Our study had a few limitations. The cross-sectional analysis of retrospective data limited our ability to draw causal inferences. We lacked access to clinical information including presence of symptoms, time since symptoms onset, vaccination status, or comorbidities. As we restricted our analysis to tests done in the Cobas-6800 platform and Cobas® SARS-CoV-2 test, the Delta variant period had a smaller number of samples. Our dataset included tests performed in a single institution, and even though they comprised samples collected in ten different cities in the Sao Paulo State, generalizability to other populations or settings may be limited. Finally, the absence of individual-level information on the infecting variant lead us to estimate the most likely predominance using external sources, which could result in residual confounding.

Despite these limitations, the strengths of this study include its large volume of samples, obtained from a reference laboratory located in one of the most affected countries worldwide, and analyzed by only one type of RT-PCR assay. The lower Ct values observed in samples from older adults during the Gamma and Delta variants predominance periods, and in samples collected at emergency units is relevant for developing control strategies against disease spread aimed at individuals with characteristics associated with higher viral loads and, thus, higher risk of virus transmission.

CONCLUSION

Our findings point to a significant association between SARS-CoV-2 Ct values and age, pandemic period and sample collection setting, with lower Ct values associated with older age groups, the Gamma and Delta periods, and samples collected in emergency units. We failed to find statistically significant associations between Ct values and gender. Moreover, we found no clear association between average SARS-CoV-2 Ct values and COVID-19 incidence, mortality, or vaccination coverage. Our results suggest that SARS-CoV-2 Ct values, a proxy for viral load and transmissibility, can be influenced by some demographic and epidemiological variables. Identifying factors associated with higher COVID-19 transmissibility as well as other infections could help to develop targeted prevention strategies.

ACKNOWLEDGMENTS

We would like thank Denise dos Anjos Laurentis S. Campos for contributing with data extraction.

REFERENCES

  • 1
    Phillips MC, Quintero D, Wald-Dickler N, Holtom P, Butler-Wu SM. SARS-CoV-2 cycle threshold (Ct) values predict future COVID-19 cases. J Clin Virol. 2022;150-151:105153.
  • 2
    Tom MR, Mina MJ. To interpret the SARS-CoV-2 test, consider the cycle threshold value. Clin Infect Dis. 2020;71:2252-4.
  • 3
    Rhoads D, Peaper DR, She RC, Nolte FS, Wojewoda CM, Anderson NW, et al. College of American Pathologists (CAP) Microbiology Committee perspective: caution must be used in interpreting the cycle threshold (Ct) value. Clin Infect Dis. 2021;72:e685-6.
  • 4
    Jones TC, Mühlemann B, Veith T, Biele G, Zuchowski M, Hofmann J, et al. Estimating infectiousness throughout SARS-CoV-2 infection course. Science. 2021;373:180.
  • 5
    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.
  • 6
    La Scola B, Le Bideau M, Andreani J, Hoang VT, Grimaldier C, Colson P, et al. Viral RNA load as determined by cell culture as a management tool for discharge of SARS-CoV-2 patients from infectious disease wards. Eur J Clin Microbiol Infect Dis. 2020;39:1059-61.
  • 7
    Bullard J, Dust K, Funk D, Strong JE, Alexander D, Garnett L, et al. Predicting infectious severe acute respiratory syndrome coronavirus 2 from diagnostic samples. Clin Infect Dis. 2020;71:2663-6.
  • 8
    Penney J, Jung A, Koethe B, Doron S. Cycle threshold values and SARS-CoV-2: relationship to demographic characteristics and disease severity. J Med Virol. 2022;94:3978-81.
  • 9
    Hurst JH, Heston SM, Chambers HN, Cunningham HM, Price MJ, Suarez L, et al. Severe acute respiratory syndrome coronavirus 2 infections among children in the biospecimens from respiratory virus-exposed kids (BRAVE Kids) study. Clin Infect Dis. 2021;73:e2875-82.
  • 10
    Baggio S, L’Huillier AG, Yerly S, Bellon M, Wagner N, Rohr M, et al. Severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) viral load in the upper respiratory tract of children and adults with early acute COVID-19. Clin Infect Dis. 2021;73:148-50.
  • 11
    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.
  • 12
    Magleby R, Westblade LF, Trzebucki A, Simon MS, Rajan M, Park J, et al. Impact of severe acute respiratory syndrome coronavirus 2 viral load on risk of intubation and mortality among hospitalized patients with coronavirus disease 2019. Clin Infect Dis. 2021;73:e4197-205.
  • 13
    Huang JT, Ran RX, Lv ZH, Feng LN, Ran CY, Tong YQ, et al. Chronological changes of viral shedding in adult inpatients with COVID-19 in Wuhan, China. Clin Infect Dis. 2020;71:2158-66.
  • 14
    Yu X, Sun S, Shi Y, Wang H, Zhao R, Sheng J. SARS-CoV-2 viral load in sputum correlates with risk of COVID-19 progression. Crit Care. 2020;24:170.
  • 15
    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.
  • 16
    Jang S, Rhee JY, Wi YM, Jung BK. Viral kinetics of SARS-CoV-2 over the preclinical, clinical, and postclinical period. Int J Infect Dis. 2021;102:561-5.
  • 17
    Kissler SM, Fauver JR, Mack C, Tai CG, Breban MI, Watkins AE, et al. Viral dynamics of SARS-CoV-2 variants in vaccinated and unvaccinated persons. N Engl J Med. 2021;385:2489-91.
  • 18
    Singanayagam A, Patel M, Charlett A, Lopez Bernal J, Saliba V, Ellis J, et al. Duration of infectiousness and correlation with RT-PCR cycle threshold values in cases of COVID-19, England, January to May 2020. Euro Surveill. 2020;25:2001483.
  • 19
    Volz E, Mishra S, Chand M, Barrett JC, Johnson R, Geidelberg L, et al. Assessing transmissibility of SARS-CoV-2 lineage B.1.1.7 in England. Nature. 2021;593:266-9.
  • 20
    Volz E, Hill V, McCrone JT, Price A, Jorgensen D, O’Toole A, et al. Evaluating the effects of SARS-CoV-2 spike mutation D614G on transmissibility and pathogenicity. Cell. 2021;184:64-75.e11.
  • 21
    Voloch CM, Francisco Jr RS, Almeida LG, Cardoso CC, Brustolini OJ, Gerber AL, et al. Genomic characterization of a novel SARS-CoV-2 lineage from Rio de Janeiro, Brazil. J Virol. 2021;95:e00119-21.
  • 22
    Resende PC, Delatorre E, Gräf T, Mir D, Motta FC, Appolinario LR, et al. Evolutionary dynamics and dissemination pattern of the SARS-CoV-2 lineage B.1.1.33 during the early pandemic phase in Brazil. Front Microbiol. 2020;11:615280.
  • 23
    Wolf JM, Kipper D, Borges GR, Streck AF, Lunge VR. Temporal spread and evolution of SARS-CoV-2 in the second pandemic wave in Brazil. J Med Virol. 2022;94:926-36.
  • 24
    Sarsomics. SARS-OMICS: atlas multi-ômico da COVID-19 brasileira. [cited 2023 Jul 17]. Available from: https://sarsomics.com
    » https://sarsomics.com
  • 25
    Fundação Oswaldo Cruz. Dashboard genomic network. [cited 2023 Jul 17]. Available from: https://www.genomahcov.fiocruz.br/dashboard-en/
    » https://www.genomahcov.fiocruz.br/dashboard-en/
  • 26
    Giovanetti M, Fonseca V, Wilkinson E, Tegally H, San EJ, Althaus CL, et al. Replacement of the Gamma by the Delta variant in Brazil: impact of lineage displacement on the ongoing pandemic. Virus Evol. 2022;8:veac024.
  • 27
    World Health Organization. Classification of Omicron (B.1.1.529): SARS-CoV-2 variant of concern. [cited 2023 Jul 17]. Available from: https://www.who.int/news/item/26-11-2021-classification-of-omicron- (b.1.1.529)-sars-cov-2-variant-of-concern
    » https://www.who.int/news/item/26-11-2021-classification-of-omicron- (b.1.1.529)-sars-cov-2-variant-of-concern
  • 28
    World Health Organization. WHO Coronavirus (COVID-19) dashboard. [cited 2023 Jul 17]. Available from: https://covid19.who.int
    » https://covid19.who.int
  • 29
    São Paulo. Governo do Estado. Fundação Sistema Estadual de Análise de Dados. SP contra o novo coronavírus: boletim completo. [cited 2023 Jul 17]. Available from: https://www.seade.gov.br/coronavirus/
    » https://www.seade.gov.br/coronavirus/
  • 30
    São Paulo. Governo do Estado. #Vacina Já. [cited 2023 Jul 17]. Available from: https://www.vacinaja.sp.gov.br/vacinometro/
    » https://www.vacinaja.sp.gov.br/vacinometro/
  • 31
    Pujadas E, Chaudhry F, McBride R, Richter F, Zhao S, Wajnberg A, et al. SARS-CoV-2 viral load predicts COVID-19 mortality. Lancet Respir Med. 2020;8:e70.
  • 32
    Salvatore PP, Dawson P, Wadhwa A, Rabold EM, Buono S, Dietrich EA, et al. Epidemiological correlates of polymerase chain reaction cycle threshold values in the detection of severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2). Clin Infect Dis. 2021;72:e761-7.
  • 33
    Shah VP, Farah WH, Hill JC, Hassett LC, Binnicker MJ, Yao JD, et al. Association between SARS-CoV-2 cycle threshold values and clinical outcomes in patients with COVID-19: a systematic review and meta-analysis. Open Forum Infect Dis. 2021;8:ofab453.
  • 34
    Martínez MJ, Basile L, Sisó-Almirall A, Cristino V, Cuesta G, Hurtado JC, et al. Lack of prognostic value of SARS-CoV2 RT-PCR cycle threshold in the community. Infect Dis Ther. 2022;11:587-93.
  • 35
    Heald-Sargent T, Muller WJ, Zheng X, Rippe J, Patel AB, Kociolek LK. Age-related differences in nasopharyngeal severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) levels in patients with mild to moderate coronavirus disease 2019 (COVID-19). JAMA Pediatr. 2020;174:902-3.
  • 36
    Brasil. Ministério da Saúde. Agência Nacional de Vigilância Sanitária. Resolução - RDC Nº 377, de 28 de abril de 2020. Autoriza, em caráter temporário e excepcional, a utilização de “testes rápidos” (ensaios imunocromatográficos) para a COVID-19 em farmácias, suspende os efeitos do § 2º do art. 69 e do art. 70 da Resolução de Diretoria Colegiada - RDC nº 44, de 17 de agosto de 2009. Diário Oficial da União, Brasília, 29 abr. 2020, Seção 1:56. [cited 2023 Jul 17]. Available from: https://www.in.gov.br/en/web/dou/-/resolucao-rdc-n-377-de-28-de-abril-de-2020-254429215
    » https://www.in.gov.br/en/web/dou/-/resolucao-rdc-n-377-de-28-de-abril-de-2020-254429215
  • 37
    Tso CF, Garikipati A, Green-Saxena A, Mao Q, 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.
  • 38
    Instituto Brasileiro de Geografia e Estatística. Brasil/São Paulo: panorama. [cited 2023 Jul 17]. Available from: https://cidades.ibge.gov.br/brasil/sp/panorama
    » https://cidades.ibge.gov.br/brasil/sp/panorama
  • 39
    Brasil. Ministério da Saúde. Agência Nacional de Vigilância Sanitária. Resolução - RDC Nº 595, de 28 de janeiro de 2022. Dispõe sobre os requisitos e procedimentos para a solicitação de registro, distribuição, comercialização e utilização de dispositivos médicos para diagnóstico in vitro como autoteste para detecção de antígeno do SARS-CoV-2, em consonância ao Plano Nacional de Expansão da Testagem para Covid-19 (PNE-Teste), e dá outras providências. Diário Oficial da União, Brasília, 28 jan. 2022, Seção 1:Extra A. [cited 2023 Jul 17]. Available from: https://www.in.gov.br/en/web/dou/-/resolucao-rdc-n-595-de-28-de-janeiro-de-2022-376825970
    » https://www.in.gov.br/en/web/dou/-/resolucao-rdc-n-595-de-28-de-janeiro-de-2022-376825970
  • FUNDING: BMF received funding from the Institutional Scientific Initiation Scholarship Program (IC-PIBIC), grant no. 143205/2022-6, funded by the Brazilian National Council for Scientific and Technological Development (CNPq), from September 2022 to August 2023. This research received grant from ITpS, Termo de Cooperacao - Parceiros, Hospital Albert Einstein, Nº C2022-0751.

Publication Dates

  • Publication in this collection
    20 Oct 2023
  • Date of issue
    2023

History

  • Received
    14 Apr 2023
  • Accepted
    17 July 2023
Instituto de Medicina Tropical de São Paulo Av. Dr. Enéas de Carvalho Aguiar, 470, 05403-000 - São Paulo - SP - Brazil, Tel. +55 11 3061-7005 - São Paulo - SP - Brazil
E-mail: revimtsp@usp.br