Abstract
In 2021, Qatar experienced considerable incidence of SARS-CoV-2 infection that was dominated sequentially by the Alpha, Beta, and Delta variants. Using the cycle threshold (Ct) value of an RT-qPCR-positive test to proxy the inverse of infectiousness, we investigated infectiousness of SARS-CoV-2 infections by variant, age, sex, vaccination status, prior infection status, and reason for testing in a random sample of 18,355 RT-qPCR-genotyped infections. Regression analyses were conducted to estimate associations with the Ct value of RT-qPCR-positive tests. Compared to Beta infections, Alpha and Delta infections demonstrated 2.56 higher Ct cycles (95% CI: 2.35-2.78), and 4.92 fewer cycles (95% CI: 4.67-5.16), respectively. The Ct value declined gradually with age and was especially high for children <10 years of age, signifying lower infectiousness of small children. Children <10 years of age had 2.18 higher Ct cycles (95% CI: 1.88-2.48) than those 10-19 years of age. Compared to unvaccinated individuals, the Ct value was higher among individuals who had received one or two vaccine doses, but the Ct value decreased gradually with time since the second-dose vaccination. Ct value was 2.07 cycles higher (95% CI: 1.42-2.72) for those with a prior infection than those without prior infection. The Ct value was lowest among individuals tested because of symptoms and was highest among individuals tested as a travel requirement. Delta was substantially more infectious than Beta. Prior immunity, whether due to vaccination or prior infection, is associated with lower infectiousness of breakthrough infections, but infectiousness increases gradually with time since the second-dose vaccination.
Introduction
The severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) pandemic continues with progressive viral evolution more than two years after it first emerged.(1) Between January 18, 2021 and May 31, 2021, Qatar experienced a SARS-CoV-2 Alpha(2) (B.1.1.7) variant wave(3) that was immediately followed by a Beta(2) (B.1.351) variant wave.(4) Starting in June 2021, the Delta(2) (B.1.617.2) variant dominated a prolonged low-incidence phase that persisted until November of 2021.(5-7) We investigated the effects of SARS-CoV-2 variant, age, vaccination, and prior infection on infectiousness of SARS-CoV-2 infections.
The real-time (quantitative) reverse transcription-polymerase chain reaction (RT-qPCR) cycle threshold (Ct) value of an RT-qPCR-positive SARS-CoV-2 test is a measure of the inverse of viral load and correlates strongly with culturable virus.(8) Therefore, it can be used to proxy inverse of SARS-CoV-2 infectiousness.(8-13) Higher Ct values signify lower infectiousness.(8-13) We assessed differences in Ct values in a random sample of 18,355 RT-qPCR-genotyped SARS-CoV-2 infections in relation to variant status, age, sex, vaccination status, prior infection status, and reason for testing.
Methods
Study population, data sources, and study design
This cross sectional study was conducted in the resident population of Qatar, applying a methodology used recently to investigate effects of the BA.1/BA.2 subvariant, vaccination, and prior infection on infectiousness of SARS-CoV-2 Omicron(2) (B.1.1.529) infections.(14) Similarly, several effects on the infectiousness of SARS-CoV-2 infections were investigated including pre-Omicron variants (Alpha, Beta, and Delta), mRNA COVID-19 vaccination status (BNT162b2(15) and mRNA-1273(16)), prior infection status, reason for RT-qPCR testing, study-period month of the RT-qPCR test (to account for the evolving phase of SARS-CoV-2 incidence), and demographic factors, including sex, age, and nationality.
The present study was conducted on a sample of 18,355 SARS-CoV-2 RT-qPCR-positive swabs that were collected randomly on a weekly basis from among all RT-qPCR-confirmed infections in Qatar between March 23, 2021 and November 6, 2021. These documented infections were RT-qPCR genotyped as part of a national project for surveillance of SARS-CoV-2 variants in Qatar.(7, 17-19) Details of laboratory methods for RT-qPCR testing and variant ascertainment are provided in Section 1 of the Supplementary Material. Coronavirus disease 2019 (COVID-19) laboratory testing, vaccination, clinical infection, and demographic data for this population were extracted from the national, federated SARS-CoV-2 databases, which include all RT-qPCR testing, reason for RT-qPCR testing, COVID-19 vaccinations, and related demographic details since the start of the pandemic. Further description of Qatar’s national COVID-19 databases has been reported previously.(4, 6, 7, 13, 20)
Every SARS-CoV-2 RT-qPCR test conducted in Qatar is classified based on the reason for testing (clinical symptoms, contact tracing, surveys or random testing campaigns, individual requests, routine healthcare testing, pre-travel, at port of entry, or other). RT-qPCR testing is performed at a mass scale and most infections are diagnosed not for appearance of symptoms, but because of routine testing.(6) Qatar has unusually young, diverse demographics, in that only 9% of its residents are ≥50 years of age, and 89% are expatriates from over 150 countries.(20, 21) Nearly all individuals were vaccinated in Qatar; however, vaccinations performed elsewhere were recorded in the health system at the port of entry upon arrival to Qatar, per national requirements.
For standardization of RT-qPCR Ct values, we analyzed only RT-qPCR-confirmed infections diagnosed with the TaqPath COVID-19 Combo Kit (Thermo Fisher Scientific, USA(22)). For each individual, all RT-qPCR-positive swabs during the study period were included, provided that at least 90 days had elapsed between two consecutive positive swabs to avoid inclusion of positive results from the same infectious episode.(23, 24) A summary measure was derived for the primary outcome, the RT-qPCR Ct value,(25) by averaging Ct values of the N, ORF1ab, and S gene targets. This average Ct value was used as the dependent variable in all analyses.
Both vaccination status and prior infection status were ascertained at the time of the RT-qPCR test. Vaccination status was defined by the number of administered vaccine doses and months elapsed since the last vaccine dose, with one month defined as 30 days. Only individuals vaccinated with BNT162b2(15) or mRNA-1273(16) vaccines were included in the analyses, as these have been the vaccines of choice in the COVID-19 immunization program in Qatar.(26-28) Rare occurrences of mixed vaccination regimens were excluded. Nearly all vaccinated persons received their second vaccine dose per protocol(15, 16) within 30 days of first dose. History of prior infection was defined as an RT-qPCR-positive test that occurred ≥90 days before the study RT-qPCR-positive test.(4, 5, 29, 30) An RT-qPCR-positive test that occurred <90 days prior to the study RT-qPCR-positive test was not considered a prior infection, but was considered a category of its own. This is because the prior RT-qPCR-positive test and the study RT-qPCR-positive test may both reflect the same infection.(23, 24) A small number of RT-qPCR tests that had no recorded Ct value were excluded from the analysis, but these constituted only 0. 5% of all RT-qPCR-genotyped infections. Otherwise, data for the remaining study variables were complete.
Oversight
Hamad Medical Corporation and Weill Cornell Medicine-Qatar Institutional Review Boards approved this retrospective study with a waiver of informed consent. The study was reported following the Strengthening the Reporting of Observational Studies in Epidemiology (STROBE) guidelines. The STROBE checklist is found in Supplementary Table 1.
Statistical analysis
Frequency distributions and measures of central tendency were used to describe the study population with respect to a priori determined factors. These included SARS-CoV-2 variant, vaccination status (factoring dose number and months since vaccination), prior infection status, reason for RT-qPCR testing, study-period month of the RT-qPCR test, and demographic factors, age, sex, and nationality.
Association of each of these factors with Ct value was assessed using univariable linear regression analyses. Unadjusted β coefficients, 95% confidence intervals (CIs), and the F-test of overall covariate significance were reported. Adjusted β coefficients and associated 95% CIs and p-values were estimated using multivariable linear regression analyses that included all covariates in the model.
Two-sided p-value <0.05 indicated statistical significance. Interactions were not considered. Statistical analyses were conducted in STATA/SE version 16.(31)
Results
Figure 1 shows the process of selecting the study population and Table 1 describes study population characteristics. This was a national study involving a random sample of 18,355 RT-qPCR-confirmed SARS-CoV-2 infections. Therefore, the study population is broadly representative of the population of Qatar. The sample included 3,347 (18.2%) Alpha infections, 5,576 (30.4%) Beta infections, and 9,432 (51.4%) Delta infections (Table 1).
Characteristics of the 18,355 RT-qPCR-genotyped SARS-CoV-2 infections between March 23 and November 6, 2021.
Abbreviations: COVID-19, coronavirus disease 2019; RT-qPCR, real-time reverse-transcription polymerase chain reaction; SARS-CoV-2, severe acute respiratory syndrome coronavirus 2.
Compared to Beta infections, Alpha and Delta infections were associated with 2.56 higher Ct cycles (95% CI: 2.35-2.78), and 4.92 fewer cycles (95% CI: 4.67-5.16) (Table 2), respectively, indicating the highest infectiousness for the Delta variant.
Associations with RT-qPCR Ct value among the 18,355 RT-qPCR-genotyped SARS-CoV-2 infections between March 23 and November 6, 2021.
Ct value declined gradually with age and was especially high for children <10 years of age, signifying lower infectiousness of small children. Children <10 years of age had 2.18 higher Ct cycles (95% CI: 1.88-2.48) than those 10-19 years of age (Table 2). The 10-19 age group was chosen as a reference, and not the <10 age group, because of the different manifestations of this infection in small children. Males had higher Ct values than females and there were some differences in Ct value by nationality.
Compared to unvaccinated individuals, Ct value was higher among individuals who received one or two vaccine doses (Table 2). However, Ct value decreased gradually with time since second-dose vaccination. Very few individuals received a booster dose during the study period (Table 1) to allow for estimation of effect of booster vaccination on Ct value. Ct value was 2.07 cycles higher (95% CI: 1.42-2.72) for those with a prior infection compared to those without prior infection.
The Ct value was lowest when testing was performed due to suspicion of infection exposure (Table 2), such as appearance of symptoms or recent exposure to an infected person (contact tracing). The Ct value was highest for infections diagnosed because of routine testing for reasons unrelated to infection exposure, such as in a random survey or because of travel re quirements. Stratified analyses for Alpha (Table 3), Beta (Table 4), and Delta (Table 5) infections suggested similar findings.
Associations with RT-qPCR Ct value among 3,347 Alpha RT-qPCR-genotyped SARS-CoV-2 infections between March 23 and November 6, 2021.
Associations with RT-qPCR Ct value among 5,576 Beta RT-qPCR-genotyped SARS-CoV-2 infections between March 23 and November 6, 2021.
Associations with RT-qPCR Ct value among 9,432 Delta RT-qPCR-genotyped SARS-CoV-2 infections between March 23 and November 6, 2021.
Discussion
Delta infections were associated with considerably lower Ct values than Beta infections, indicating higher infectiousness of this variant, perhaps because of higher viral load and/or longer duration of infection. This appears to be the first direct comparison of the infectiousness of Delta versus Beta infections, and supports the high infectiousness of the Delta variant compared to pre-Omicron variants such as Alpha, as reported previously.(32, 33) This difference in viral load between Delta and Beta appears also to extend to severity of infection, as Delta infections were found associated with higher severity than Beta infections.(34) Of note that Beta infections were also found earlier to be associated with higher severity than Alpha infections.(35)
Prior immunity against SARS-CoV-2 infection, whether due to vaccination or prior infection, was associated with higher Ct value at infection, and thus lower infectiousness of breakthrough infections. This confirms earlier findings,(13, 14) and suggests that strength of immunity is manifest not only in protection against infection, but also against the infectiousness, if a breakthrough infection occurs.(13) However, this effect appeared to depend on the time since the prior immunological event (Table 2). Ct values decreased gradually with time since second-dose vaccination, paralleling the established pattern of waning of vaccine effectiveness after the second dose.(6, 36, 37)
Ct values decreased with age, perhaps reflecting slower virus clearance with aging (38) and confirming our earlier findings.(14) Ct values were particularly high for infections among small children <10 years of age. This finding supports a lesser role for small children than adults in the transmission of infection, as suggested in studies of secondary transmission within households.(39-41) There were differences in Ct value by sex and nationality, but these may be a consequence of different test-seeking behaviors for different socio-economic groups in Qatar’s diverse population. Ct values also varied by reason for testing, with lower Ct values of infections diagnosed because of suspicion of infection, and higher Ct values of infections diagnosed through routine testing unrelated to infection exposure. This finding also confirms our earlier finding for Omicron infections.(14)
The study has limitations. While both Delta and Beta infections had lower Ct values than Alpha infections, this finding should be interpreted with caution. RT-qPCR genotyping in Qatar started only after the Alpha wave peaked in the first week of March 2021. Many Alpha infections may have been older or prolonged infections, rather than recent infections, explaining the higher Ct values of these infections.
Eligible individuals were selected from the pool of those who had documented RT-qPCR-confirmed infections, but other infections may have occurred that were never documented. It is possible that infections in those with prior infections or those vaccinated are less likely to be diagnosed, perhaps because of minimal or no symptoms.(13) Nevertheless, RT-qPCR testing in Qatar is done at a mass scale, where a significant proportion of the population is being tested every week.(6) The majority of infections are identified not because of symptoms, but because of routine testing for other reasons (Table 1).(6) The date of symptom onset was not available for symptomatic cases. Therefore, an analysis factoring the duration between symptom onset and RT-qPCR test was not possible.
A small number of RT-qPCR tests had missing Ct values and were excluded from the analysis, but these constituted only 0.5% of all RT-qPCR-genotyped infections. The study population consisted mostly of working-age adults; thus, the results may not be generalizable to the elderly. Too few individuals received the booster dose during the study period to allow for estimation of effect of booster vaccination on Ct values.
In conclusion, the Delta variant appears substantially more infectious than the Beta variant, explaining its global reach in the pre-Omicron era. Infectiousness of SARS-CoV-2 infections increases with age, apparently reflecting slower virus clearance with aging. Prior immunity against SARS-CoV-2 infection, whether due to vaccination or prior infection, is associated with lower infectiousness of breakthrough infections. However, infectiousness of breakthrough infections increases gradually with time since second-dose vaccination, paralleling the waning of vaccine effectiveness after the second dose.
Data Availability
The dataset of this study is a property of the Qatar Ministry of Public Health that was provided to the researchers through a restricted-access agreement that prevents sharing the dataset with a third party or publicly. Future access to this dataset can be considered through a direct application for data access to Her Excellency the Minister of Public Health (https://www.moph.gov.qa/english/Pages/default.aspx). Aggregate data are available within the manuscript and its Supplementary information.
Author contributions
SHQ co-designed the study, performed the statistical analyses, and co-wrote the first draft of the article. PT and MRH conducted the multiplex, RT-qPCR variant screening and viral genome sequencing. HC co-designed the study and supported the statistical analyses. LJA conceived and co-designed the study, led the statistical analyses, and co-wrote the first draft of the article. All authors contributed to data collection and acquisition, database development, discussion and interpretation of the results, and to the writing of the manuscript. All authors have read and approved the final manuscript.
Competing interests
Dr. Butt has received institutional grant funding from Gilead Sciences unrelated to the work presented in this paper. Otherwise, we declare no competing interests.
Supplementary Material
1. Laboratory methods and variant ascertainment
1.1 Real-time reverse-transcription polymerase chain reaction testing
Nasopharyngeal and/or oropharyngeal swabs were collected for polymerase chain reaction (PCR) testing and placed in Universal Transport Medium (UTM). Aliquots of UTM were: 1) extracted on KingFisher Flex (Thermo Fisher Scientific, USA), MGISP-960 (MGI, China), or ExiPrep 96 Lite (Bioneer, South Korea) followed by testing with real-time reverse-transcription PCR (RT-qPCR) using TaqPath COVID-19 Combo Kits (Thermo Fisher Scientific, USA) on an ABI 7500 FAST (Thermo Fisher Scientific, USA); 2) tested directly on the Cepheid GeneXpert system using the Xpert Xpress SARS-CoV-2 (Cepheid, USA); or 3) loaded directly into a Roche cobas 6800 system and assayed with the cobas SARS-CoV-2 Test (Roche, Switzerland). The first assay targets the viral S, N, and ORF1ab gene regions. The second targets the viral N and E-gene regions, and the third targets the ORF1ab and E-gene regions.
All PCR testing was conducted at the Hamad Medical Corporation Central Laboratory or Sidra Medicine Laboratory, following standardized protocols. All PCR testing was performed with extensively used, investigated, and validated commercial platforms having essentially 100% sensitivity and specificity.
1.2 Classification of infections by variant type
Surveillance for SARS-CoV-2 variants in Qatar is mainly based on viral genome sequencing and multiplex RT-qPCR variant screening(1) of random positive clinical samples,(2-7) complemented by deep sequencing of wastewater samples.(4, 8)
The accuracy of the RT-qPCR genotyping was verified against either Sanger sequencing of the receptor-binding domain (RBD) of SARS-CoV-2 surface glycoprotein (S) gene, or by viral whole-genome sequencing on a Nanopore GridION sequencing device. From 236 random samples (27 Alpha-like, 186 Beta-like, and 23 “other” variants), PCR genotyping results for Alpha-like, Beta-like, and ‘other’ variants were in 88.8% (23 out of 27), 99.5% (185 out of 186), and 100% (23 out of 23) agreement with the SARS-CoV-2 lineages assigned by sequencing.
Within the “other” variant category, Sanger sequencing and/or Illumina sequencing of the RBD of SARS-CoV-2 spike gene on 728 random samples confirmed that 701 (96.3%) were Delta cases and 17 (2.3%) were other variant cases, with 10 (1.4%) samples failing lineage assignment.6,8 Accordingly, a Delta case was proxied as any “other” case identified through the RT-qPCR based variant screening.
All the variant RT-qPCR screening was conducted at the Sidra Medicine Laboratory following standardized protocols.
2. Supplementary tables
Acknowledgements
We acknowledge the many dedicated individuals at Hamad Medical Corporation, the Ministry of Public Health, the Primary Health Care Corporation, Qatar Biobank, Sidra Medicine, and Weill Cornell Medicine-Qatar for their diligent efforts and contributions to make this study possible. The authors are grateful for institutional salary support from the Biomedical Research Program and the Biostatistics, Epidemiology, and Biomathematics Research Core, both at Weill Cornell Medicine-Qatar, as well as for institutional salary support provided by the Ministry of Public Health, Hamad Medical Corporation, and Sidra Medicine. The authors are also grateful for the Qatar Genome Programme and Qatar University Biomedical Research Center for institutional support for the reagents needed for the viral genome sequencing. The funders of the study had no role in study design, data collection, data analysis, data interpretation, or writing of the article. Statements made herein are solely the responsibility of the authors.
References
