Abstract
BACKGROUND Qatar experienced a large SARS-CoV-2 Omicron (B.1.1.529) wave that started on December 19, 2021 and peaked in mid-January, 2022. We investigated effects of Omicron subvariant (BA.1 and BA.2), previous vaccination, and prior infection on infectiousness of Omicron infections, between December 23, 2021 and February 20, 2022.
METHODS Univariable and multivariable regression analyses were conducted to estimate the association between the RT-qPCR cycle threshold (Ct) value of PCR tests (a proxy for SARS-CoV-2 infectiousness) and each of the Omicron subvariants, mRNA vaccination, prior infection, reason for RT-qPCR testing, calendar week of RT-qPCR testing (to account for phases of the rapidly evolving Omicron wave), and demographic factors.
RESULTS Compared to BA.1, BA.2 was associated with 3.53 fewer cycles (95% CI: 3.46-3.60), signifying higher infectiousness. Ct value decreased with time since second and third vaccinations. Ct values were highest for those who received their boosters in the month preceding the RT-qPCR test—0.86 cycles (95% CI: 0.72-1.00) higher than for unvaccinated persons. Ct value was 1.30 (95% CI: 1.20-1.39) cycles higher for those with a prior infection compared to those without prior infection, signifying lower infectiousness. Ct value declined gradually with age. Ct value was lowest for those who were tested because of symptoms and was highest for those who were tested for travel-related purposes. Ct value was lowest during the exponential-growth phase of the Omicron wave and was highest after the wave peaked and was declining.
CONCLUSIONS The BA.2 subvariant appears substantially more infectious than the BA.1 subvariant. This may reflect higher viral load and/or longer duration of infection, thereby explaining the rapid expansion of this subvariant in Qatar.
Introduction
Qatar experienced a large severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) Omicron (B.1.1.529)1 wave that started on December 19, 2021 and peaked in mid-January, 2022.2-5 We investigated effects of Omicron subvariant (BA.1 and BA.2), previous vaccination, and prior infection on infectiousness of Omicron infections, between December 23, 2021 and February 20, 2022. Incidence was initially dominated by BA.1, but within a few days, BA.2 predominated (Figure 1). Incidence of the Delta variant was minimal and no other variants were detected in viral genome sequencing and real-time reverse-transcription polymerase chain reaction (RT-qPCR) genotyping of randomly collected samples (Section S1 of Supplementary Appendix).
Methods
Study population, data sources, and study design
The RT-qPCR cycle threshold (Ct) value is a measure of the inverse of viral load and correlates strongly with culturable virus;6 thus, it can be used as a proxy for SARS-CoV-2 infectiousness.6-11 We investigated several effects on the RT-qPCR Ct value of Omicron variant infections in the resident population of Qatar using a cross-sectional study design. These included: subvariant (BA.1 versus BA.2), mRNA (BNT162b212 and mRNA-127313) vaccination status, prior infection status, reason for RT-qPCR testing, study-period week of the RT-qPCR test (to account for the evolving phase of the rapid Omicron wave), and demographic factors including sex, age, and nationality.
The study population included all individuals with an RT-qPCR-confirmed SARS-CoV-2 infection in Qatar between December 23, 2021 and February 20, 2022. 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 can be found in previous publications.11,14-17
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.15 About 75% of those diagnosed over recent months were diagnosed not because of appearance of symptoms, but because of routine testing.15 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.14,18 Nearly all individuals were vaccinated in Qatar, however, vaccinations performed elsewhere were still recorded in the health system at the port of entry upon arrival to Qatar per country requirements.
Informed by the viral genome sequencing and the RT-qPCR genotyping (Section S1), a SARS-CoV-2 infection with the BA.1 subvariant was proxied as an S-gene “target failure” (SGTF) case using the TaqPath COVID-19 Combo Kit (Thermo Fisher Scientific, USA19) that tests for the S-gene and is affected by the del69/70 mutation in the S-gene.20 A SARS-CoV-2 infection with the BA.2 subvariant was proxied as a non-SGTF case using this TaqPath Kit. For ascertainment of subvariant status and standardization of RT-qPCR Ct values, we analyzed only the RT-qPCR-confirmed infections diagnosed with this TaqPath Kit.
For each individual, we selected only the first positive RT-qPCR-confirmed swab during the study period. We subsequently derived a summary measure for our primary outcome, the RT-qPCR Ct value, by averaging the Ct values of the N, ORF1ab, and S (if not an S-gene “target failure” case) genes. 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 factoring the number of administered vaccine doses and months elapsed since the last vaccine dose, with one month defined as 30 days. Only vaccination with BNT162b212 and mRNA-127313 vaccines were considered in the analyses, as these have been the vaccines of choice in the COVID-19 immunization program in Qatar.21-23 Rare occurrences of mixed vaccination regimens were excluded. Prior infection was defined as an RT-qPCR-positive test that occurred ≥90 days before the study RT-qPCR-positive test.3,17,24-31 An RT-qPCR-positive test that occurred <90 days prior to the study RT-qPCR-positive test was still factored in the analysis, but was not considered a prior infection. This is because this RT-qPCR-positive test and the study RT-qPCR-positive test may both reflect the same prolonged infection.32-34 A small number of RT-qPCR tests had no recorded Ct value and were thus excluded from the analysis, but these constituted only 0.1% of all RT-qPCR tests. Otherwise, data on the remaining study variables were complete.
Oversight
Hamad Medical Corporation and Weill Cornell Medicine-Qatar Institutional Review Boards approved this retrospective study with waiver of informed consent. The study was reported following STROBE guidelines. The STROBE checklist is found in Table S1.
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 Omicron infection subvariant, vaccination status (factoring dose number and months since vaccination), prior infection status, reason for RT-qPCR testing, study-period week of the RT-qPCR test, and demographic factors, namely sex, age, and nationality. Differences between BA.1 and BA.2 infections were estimated using standardized mean differences (SMDs).
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.
The 95% CIs were not adjusted for multiplicity. Two-sided p-value <0.05 indicated statistical significance. Interactions were not considered. Statistical analyses were conducted in STATA/SE version 16.35
Results
Figure 2 shows the process of selecting the study population and Table 1 describes the study population characteristics. This was a national study involving 156,202 individuals infected with Omicron who are broadly representative of the population of Qatar.
Compared to BA.1, BA.2 was associated with 3.53 fewer cycles (95% CI: 3.46-3.60), signifying higher infectiousness (Table 2). Ct value decreased with time since second and third vaccinations, mirroring the established pattern of waning vaccine effectiveness.4 Ct values were highest for those who received their boosters in the month preceding the RT-qPCR test—0.86 cycles (95% CI: 0.72-1.00) higher than for unvaccinated persons. Ct value was 1.30 (95% CI: 1.20-1.39) cycles higher for those with a prior infection compared to those without prior infection, signifying lower infectiousness.
Ct value declined gradually with age (Table 2), perhaps reflecting slower virus clearance with aging. There were differences in Ct value by sex and nationality, but these may reflect different test-seeking behaviors for different socio-economic groups in Qatar’s diverse population, or differences in the rates of prior undocumented infection by nationality.14,36-39 Ct value was lowest for those who were tested because of symptoms and was highest for those who were tested for travel-related purposes. Ct value was lowest during the exponential-growth phase of the Omicron wave, as a large proportion of infections were recent, and was highest after the wave peaked and was declining, as a small proportion of infections were recent. Stratified analyses for each of BA.1 (Table 3) and BA.2 (Table 4) showed similar findings.
Discussion
The BA.2 subvariant appears substantially more infectious than the BA.1 subvariant, consistent with recent findings of a household study from Denmark.40 This may reflect higher viral load and/or longer duration of infection, thereby explaining the rapid expansion of this subvariant in Qatar (Figure 1). Natural immunity from previous infection and strength of vaccine immunity correlate with less infectious breakthrough infections, as observed for earlier SARS-CoV-2 variants.11 Symptomatic infection and older age are associated with higher infectiousness.
Limitations and caveats
A small number of RT-qPCR tests had no available Ct value and were thus excluded from the analysis, but these constituted only 0.1% of all RT-qPCR tests. The study was implemented on documented RT-qPCR-confirmed infections, but other infections may have occurred but were never documented. It is possible that infections in those with prior infection or those vaccinated are less likely to be diagnosed, perhaps because of minimal or no symptoms. Nevertheless, RT-qPCR testing in Qatar is done at a mass scale, where a significant proportion of the population is being tested every week.15 The majority of infections are identified not because of symptoms, but because of routine testing for other reasons (Table 1).15 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.
The study population consists of mostly working-age adults and thus the results may not necessarily be generalizable to other population groups, such as the elderly. The analyses controlled for sex, age, and nationality but it was not possible to control for other factors, such as comorbidities or socio-economic factors, as data on these factors were not available. Of note that the number of individuals with severe chronic conditions is small in Qatar’s young population.14,41 The national list of vaccine prioritization included only 19,800 individuals of all age groups with serious co-morbid conditions to be prioritized in the first phase of vaccine roll-out.15 Factoring nationality in the analyses may have (partially) controlled socio-economic differences/occupational risk, in consideration of the association between nationality and occupation in Qatar.14,36-39
BA.1 and BA.2 ascertainment was based on proxy criteria; presence or absence of SGTF using the TaqPath Kit, but this method of ascertainment is well established not only for Omicron subvariants, but also for other variants such as Alpha.20,30,42 Some Omicron infections may have been misclassified Delta infections, but this is not likely, as Delta incidence was limited during the study (Section S1).
Time since vaccination was associated with lower Ct value, mirroring the established pattern of waning vaccine effectiveness.4 Unexpectedly, however, the Ct value for those who had their second dose >6 months earlier was lower than that among unvaccinated persons. With the high vaccine coverage in Qatar (exceeding 85%), the unvaccinated group may be different in other uncontrolled attributes from the vaccinated group. For example, those unvaccinated may have chosen not to receive the vaccine because of undocumented prior infection and thus are not truly immune naïve.
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. HC co-designed the study, supported the statistical analyses, and co-wrote the first draft of the article. LJA conceived and co-designed the study, led 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. HY, HAK, and MS conducted viral genome sequencing. 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 Appendix
Section S1. Laboratory methods and variant ascertainment
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.
Classification of infections by variant type
Surveillance for severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) variants in Qatar is mainly based on viral genome sequencing and multiplex RT-qPCR variant screening1 of random positive clinical samples,2-7 complemented by deep sequencing of wastewater samples.4,8
A total of 315 random SARS-CoV-2-positive specimens collected between December 19, 2021 and January 22, 2022 were viral whole-genome sequenced on a Nanopore GridION sequencing device. Of these, 300 (95.2%) were confirmed as Omicron (B.1.1.529)9 infections and 15 (4.8%) as Delta (B.1.617.2)9 infections.4,10,11 Of 286 Omicron infections with confirmed subvariant status, 68 (23.8%) were BA.1 cases and 218 (76.2%) were BA.2 cases. No Delta case was detected in sequencing after January 8, 2022, nor were other variants.
Additionally, a total of 1,315 random SARS-CoV-2-positive specimens collected between December 22, 2021 and January 1, 2022 were RT-qPCR genotyped. The RT-qPCR genotyping identified 1 B.1.617.2-like Delta case, 366 BA.1-like Omicron cases, 898 BA.2-like Omicron cases, and 50 were undetermined cases where the genotype could not be assigned.
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 147 random SARS-CoV-2-positive specimens all collected in December of 2021, RT-qPCR genotyping was able to assign a genotype in 129 samples. The agreement between RT-qPCR genotyping and sequencing was 100% for Delta (n=82), 100% for Omicron BA.1 (n=18), and 93% for Omicron BA.2 (27 of 29 were correctly assigned to BA.2 and remaining 2 specimens genotyped as BA.2 were B.1.617.2 by sequencing). Of the remaining 18 specimens: 10 failed PCR amplification and sequencing, 8 could not be assigned a genotype by RT-qPCR (4 of 8 were B.1.617.2 by sequencing, and the remaining 4 failed sequencing). All the variant RT-qPCR genotyping was conducted at the Sidra Medicine Laboratory following standardized protocols.
The large Omicron-wave exponential-growth phase in Qatar started on December 19, 2021 and peaked in mid-January, 2022.4,10-13 The study duration coincided with the intense Omicron wave where Delta incidence was limited. Accordingly, any PCR-positive test during the study duration, between December 23, 2021 and February 20, 2022, was assumed to be an Omicron infection. Of note that the study duration started on December 23, 2021, and not on December 19, 2021, to minimize the occurrence of residual Delta incidence during the first few days of the Omicron wave.
Informed by the viral genome sequencing and the RT-qPCR genotyping, a SARS-CoV-2 infection with the BA.1 subvariant was proxied as an S-gene “target failure” (SGTF) case using the TaqPath COVID-19 Combo Kit (Thermo Fisher Scientific, USA14) that tests for the S-gene and is affected by the del69/70 mutation in the S-gene.15 A SARS-CoV-2 infection with the BA.2 subvariant was proxied as a non-SGTF case using this TaqPath Kit. For ascertainment of subvariant status and standardization of RT-qPCR cycle threshold values, we analyzed only the RT-qPCR-confirmed infections diagnosed with this TaqPath Kit.
Acknowledgements
We acknowledge the many dedicated individuals at Hamad Medical Corporation, the Ministry of Public Health, the Primary Health Care Corporation, the 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 support from the Biomedical Research Program and the Biostatistics, Epidemiology, and Biomathematics Research Core, both at Weill Cornell Medicine-Qatar, as well as for 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. Statements made herein are solely the responsibility of the authors. The funders of the study had no role in study design, data collection, data analysis, data interpretation, or writing of the article.
References
