Abstract
BACKGROUND The SARS-CoV-2 Omicron (B.1.1.529) variant has two main sub-lineages, BA.1 and BA.2 with significant genetic distance between them. This study investigated protection of infection with one sub-lineage against reinfection with the other sub-lineage in Qatar during a large BA.1 and BA.2 Omicron wave, from December 19, 2021 to February 21, 2022.
METHODS Two national matched, retrospective cohort studies were conducted to estimate effectiveness of BA.1 infection against reinfection with BA.2 (N=20,197; BA.1-against-BA.2 study), and effectiveness of BA.2 infection against reinfection with BA.1 (N=100,925; BA.2-against-BA.1 study). Associations were estimated using Cox proportional-hazards regression models.
RESULTS In the BA.1-against-BA.2 study, cumulative incidence of infection was estimated at 0.03% (95% CI: 0.01-0.07%) for the BA.1-infected cohort and at 0.62% (95% CI: 0.51-0.75%) for the uninfected-control cohort, 15 days after the start of follow-up. Effectiveness of BA.1 infection against reinfection with BA.2 was estimated at 94.9% (95% CI: 88.4-97.8%). In the BA.2-against-BA.1 study, cumulative incidence of infection was estimated at 0.03% (95% CI: 0.02-0.04%) for the BA.2-infected cohort and at 0.17% (95% CI: 0.15-0.21%) for the uninfected-control cohort, 15 days after the start of follow-up. Effectiveness of BA.2 infection against reinfection with BA.1 was estimated at 85.6% (95% CI: 77.4-90.9%).
CONCLUSIONS Infection with an Omicron sub-lineage appears to induce strong, but not full protection against reinfection with the other sub-lineage, for at least several weeks after the initial infection.
Introduction
Reinfections with the severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) variants that can evade immune response are a concern, potentially challenging the global response to the pandemic.1 This is especially true of the Omicron2 (B.1.1.529) variant and its two main sub-lineages, BA.1 and BA.2, which harbor multiple mutations that can mediate immune evasion.2-4 While SARS-CoV-2 infection with earlier variants elicits >85% protection against reinfection with the Alpha2 (B.1.1.7),5-8 Beta2 (B.1.351),5,7,8 and Delta2 (B.1.617.2)7,9 variants, protection against reinfection with the Omicron BA.1 sub-lineage is inferior at <60%.7
Qatar has been experiencing a large Omicron wave that started on December 19, 2021 and peaked in mid-January, 2022.7,10-12 Initially, the BA.1 sub-lineage was predominant, but within days, the BA.2 sub-lineage predominated (Figure 1). Considering the significant genetic distance between BA.1 and BA.2, we aimed to investigate and estimate protection of prior infection with each sub-lineage against the other.
Methods
Data sources and study design
This study analyzed the national, federated databases for coronavirus disease 2019 (COVID-19), retrieved from the integrated nationwide digital-health information platform. Databases include all SARS-CoV-2-related data and associated demographic information, with no missing information, since pandemic onset. These include all polymerase chain reaction (PCR) testing and more recently, rapid antigen testing (RAT) conducted at healthcare facilities (from January 5, 2022 onwards). These also include all COVID-19 vaccination records, COVID-19 hospitalizations, infection severity and mortality classifications per World Health Organization (WHO) guidelines,13,14 in addition to sex, age, and nationality information retrieved from the national registry. Further description of these national databases can be found in previous publications.5,15-18
During this study, from December 19, 2021 to February 21, 2022, nearly all infection incidence was due to the Omicron variant. A total of 315 random SARS-CoV-2-positive specimens collected between December 19, 2021 and January 22, 2022 were viral whole-genome sequenced. Of these, 300 (95.2%) were confirmed as Omicron infections and 15 (4.8%) as Delta infections.7,10-12 Of 286 Omicron infections with confirmed sub-lineage 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. Further details about viral genome sequencing and real-time, reverse-transcription PCR (RT-qPCR) genotyping of random specimens during the study are found in Supplementary Appendix Section S1.
Informed by viral genome sequencing and RT-qPCR genotyping, a SARS-CoV-2 infection with the BA.1 sub-lineage was proxied as an S-gene “target failure” (SGTF) case using the TaqPath COVID-19 Combo Kit (Thermo Fisher Scientific, USA).19 Conversely, an infection with the BA.2 sub-lineage was proxied as a non-SGTF case using the same assay.
We assessed effectiveness of BA.1 infection against reinfection with BA.2 (denoted as ; BA.1-against-BA.2 study), and effectiveness of BA.2 infection against reinfection with BA.1 (denoted as ; BA.2-against-BA.1 study), using two matched, retrospective cohort studies. PES was defined as the proportional reduction in susceptibility to documented infection, regardless of symptoms, among those with the prior sub-lineage infection versus those without.7,8 The BA.1-against-BA.2 study followed a cohort of individuals with documented BA.1 infections and compared incidence of BA.2 infection in this cohort with that in a control cohort of individuals with no record of prior SARS-CoV-2 infection. The BA.2-against-BA.1 study followed a cohort of individuals with documented BA.2 infections and compared incidence of BA.1 infection in this cohort with that in a control cohort of individuals with no record of prior SARS-CoV-2 infection.
To optimize specificity in defining the cohorts, the BA.1-infected and BA.2-infected cohorts were defined based on existence of an infection documented only using PCR and with a PCR cycle threshold value <30, between December 19, 2021 and February 21, 2022. In all cohorts of the two studies, the two case and two control cohorts, persons with a record of a prior infection before December 19, 2021 were excluded. This is to ensure that estimated and are not affected by immunity induced by prior infections with earlier variants. Record of COVID-19 vaccination was not an exclusion criterion, but the regression analyses adjusted for vaccination status (unvaccinated, one dose, two doses, or three doses at the start of the follow-up). The control cohorts in the two studies were defined on the basis of PCR-negative tests between November 1, 2021 and December 18, 2021 (Figure 2), to ensure that all persons in these cohorts have a record of a recent active residence in Qatar.
Ideally, SARS-CoV-2 reinfection is defined as a documented infection ≥90 days after an earlier infection, to avoid misclassification of prolonged infections as reinfections, if a shorter time interval is used.20-22 Since the Omicron wave started only few weeks ago, this definition could not be used. Analysis of durations of follow-up was conducted to identify the longest time interval possible while maintaining adequate durations of follow-up and precision of estimates. Informed by the analysis, reinfection was defined as a documentation of infection ≥35 days after the prior infection. At this interval, only a small number of documented reinfections could have been prolonged prior infections rather than true reinfections.20-24 Cohorts were thus followed after completion of 35 days since documentation of the BA.1 (or BA.2) infection.
Individuals in each of the BA.1-infected and BA.2-infected cohorts were exact-matched in a 1:1 ratio by sex, 10-year age group, and nationality to uninfected individuals in control cohorts (Figure 2), to control for known differences in the risk of exposure to SARS-CoV-2 infection in Qatar.15,25-28 Matching was performed through an iterative process that ensured that each control was alive and infection-free at the start of follow-up. Follow-up was defined, for each matched pair, at ≥35 days after documentation of the BA.1 infection in the BA.1-infected cohort and BA.2 infection in the BA.2-infected cohort. Cohorts were followed up until the first of the following events: a PCR-documented BA.1 infection, a PCR-documented BA.2 infection, other PCR-documented infection (documented with an assay other than TaqPath), RAT-documented infection, death, and end of study censoring (February 21, 2022).
Laboratory methods
Laboratory methods for the RT-qPCR testing, rapid antigen testing, and viral genome sequencing are found in Section S1.
Statistical analysis
Frequency distributions and measures of central tendency were used to describe full and matched cohorts. Group comparisons were performed using standardized mean differences (SMDs), with an SMD <0.1 indicating adequate matching.29 Cumulative incidence of infection was defined as the proportion of individuals at risk whose primary endpoint was an incident infection during follow-up, and was estimated in each cohort using the Kaplan–Meier estimator method.30 Incidence rate of infection in each cohort, which was defined as the number of identified infections divided by the number of person-weeks contributed by all individuals in the cohort, was estimated, along with its 95% confidence interval (CI), using a Poisson log-likelihood regression model with the STATA 17.0 stptime command.
The hazard ratio comparing incidence of infection in case versus control cohorts and corresponding 95% CI were calculated using Cox regression adjusted for matching factors and COVID-19 vaccination status (unvaccinated, one dose, two doses, or three doses at the start of the follow-up) with the STATA 17.0 stcox command. Shoenfeld residuals and log-log plots for survival curves were used to test the proportional-hazards assumption and to investigate its adequacy. 95% CIs were not adjusted for multiplicity and should not be used to infer definitive differences between cohorts. Interactions were not considered. and were estimated using the equation: Effectiveness = 1−adjusted hazrad ratio. Statistical analyses were conducted in STATA/SE version 17.0.
Due to the large Omicron wave, use of rapid antigen testing was expanded rapidly to supplement PCR testing starting from January 5, 2022, precluding ascertainment of the Omicron sub-lineage in these tests. While 70.2% of all PCR tests (positive or negative) during the study were conducted using an assay that targets the S-gene, a minority of infections were documented with other commercial PCR kits/platforms that are not affected by the del69/70 mutation in the S-gene (Section S1), also precluding ascertainment of the Omicron sub-lineage in these tests.
The main estimate for each of and was generated by randomly assigning a sub-lineage status (BA.1 or BA.2) for each RAT-documented infection and non-TaqPath PCR-documented infection, diagnosed in a specific calendar day, on the basis of the probability of the infection being BA.1 or BA.2 in that specific day. This probability was determined by the observed distribution of identified BA.1 and BA.2 infections in each calendar day (Figure 1). was then estimated by setting the study outcome in the analysis of the BA.1-against-BA.2 study as a BA.2 infection. was estimated by setting the study outcome in the analysis of the BA.2-against-BA.1 study as a BA.1 infection.
In a sensitivity analysis, a second estimate was provided for each of and , by randomly assigning a sub-lineage status (BA.1 or BA.2) for each RAT-documented infection and non-TaqPath PCR-documented infection, in each cohort of the two studies, on the basis of the observed distribution of identified BA.1 and BA.2 infections among the infections diagnosed in the specific considered cohort.
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 STROBE guidelines. The STROBE checklist is found in Table S1.
Results
BA.1-against-BA.2 study
Figure 2 shows the population selection process for the BA.1-against-BA.2 study. Table 1 shows baseline characteristics of full and matched cohorts. The study was based on the total population of Qatar and is broadly representative of the diverse (international), but young and predominantly male, total population of Qatar (Table S2).
The median time of follow-up was 14 days (interquartile range (IQR), 12-17 days) for the BA.1-infected and 14 days (interquartile range (IQR), 11-17 days) for the uninfected-control cohorts (Figure 3). The proportion of individuals who had a PCR or RAT test during follow-up was 9.5% for the BA.1-infected cohort and 12.2% for the uninfected-control cohort.
One PCR-documented BA.2 infection, 9 PCR-documented BA.1 infections, 3 other PCR-documented infections, and 4 RAT-documented infections were recorded in the BA.1-infected cohort ≥35 days after the BA.1 infection (Figure 2). Thirty-three PCR-documented BA.2 infections, 11 PCR-documented BA.1 infections, 15 other PCR-documented infections, and 101 RAT-documented infections were recorded during the corresponding time of follow-up for the uninfected-control cohort.
In the main analysis for , cumulative incidence of infection was estimated at 0.03% (95% CI: 0.01-0.07%) for the BA.1-infected cohort and at 0.62% (95% CI: 0.51-0.75%) for the uninfected-control cohort, 15 days after the start of follow-up (Figure 3). The hazard ratio for infection, adjusted for sex, 10-year age group, nationality group, and vaccination status, was estimated at 0.05 (95% CI: 0.02-0.12) (Table 2). The effectiveness of BA.1 infection against reinfection with BA.2 was estimated at 94.9% (95% CI: 88.4-97.8%). In the sensitivity analysis (Table 2 and Figure S1), the effectiveness of BA.1 infection against reinfection with BA.2 was estimated at 98.4% (95% CI: 93.3-99.6%).
BA.2-against-BA.1 study
Figure 2 shows the population selection process for the BA.2-against-BA.1 study. Table 1 shows baseline characteristics of full and matched cohorts. The study population was representative of the population of Qatar (Table S2).
The median time of follow-up was 12 days (IQR, 8-15 days) for both the BA.2-infected cohort and the uninfected-control cohort (Figure 3). The proportion of individuals who had a PCR or RAT test during follow-up was 6.4% for the BA.2-infected cohort and 8.8% for the uninfected-control cohort.
Nineteen PCR-documented BA.1 infections, 41 PCR-documented BA.2 infections, 17 other PCR-documented infections, and 16 RAT-documented infections were recorded in the BA.2-infected cohort ≥35 days after the BA.2 infection (Figure 2). Thirty-nine PCR-documented BA.1 infections, 94 PCR-documented BA.2 infections, 80 other PCR-documented infections, and 338 RAT-documented infections were recorded during the corresponding time of follow-up for the uninfected-control cohort.
In the main analysis for , cumulative incidence of infection was estimated at 0.03% (95% CI: 0.02-0.04%) for the BA.2-infected cohort and at 0.17% (95% CI: 0.15-0.21%) for the uninfected-control cohort, 15 days after the start of follow-up (Figure 3). The adjusted hazard ratio for infection was estimated at 0.14 (95% CI: 0.09-0.23) (Table 2). The effectiveness of BA.2 infection against reinfection with BA.1 was estimated at 85.6% (95% CI: 77.4-90.9%). In the sensitivity analysis (Table 2 and Figure S1), the effectiveness of BA.2 infection against reinfection with BA.1 was estimated at 82.2% (95% CI: 73.7-88.0%).
Discussion
Reinfections with BA.2 (or BA.1) shortly after infection with BA.1 (or BA.2) have been observed in Qatar during a large Omicron wave in which both sub-lineages were intensely circulating. Indeed, 0.9% (1,062 cases) of all individuals who had a PCR-positive test with a known sub-lineage status during the Omicron wave, between December 19, 2021 and February 21, 2022, had also a subsequent PCR-positive test with the other counter sub-lineage within the same duration.
However, it is remarkable that incidence of reinfection, regardless of sub-lineage, was much lower in the BA.1-infected and BA.2-infected cohorts than incidence of infection in the corresponding uninfected-control cohorts (Figure 2), consistent with strong protection against reinfection regardless of sub-lineage. Our findings indicate that infection with an Omicron sub-lineage appears to elicit strong protection against reinfection with the other sub-lineage at an effectiveness that exceeds 85%, similar to the protection observed for infection with original virus or earlier variants (Alpha, Beta, or Delta) against reinfection with original virus or earlier variants.2,5-9,20,24
These findings, in the context of broader evidence for natural immunity,2,5-9,20-22 suggest that natural immunity of SARS-CoV-2 variants cluster into two groups: early non-Omicron variants, and Omicron BA.1 and BA.2 sub-lineages. Within each group, there appears to be strong protection against reinfection with an effectiveness that exceeds 85%. However, across groups, the protection may not exceed 60%, as was observed recently.7 This conclusion is also supported by evidence on sensitivity of variants to SARS-CoV-2 antibodies.2-4,31,32
This study has limitations. Since the Omicron wave was initially dominated with BA.1 (Figure 1), the follow-up in the BA.2-against-BA.1 study was shifted in calendar time to after the follow-up in the BA.1-against-BA.2 study. With the high intensity of infection transmission, followed by rapid decline of the Omicron wave, more of the uninfected-controls in the BA.2-against-BA.1 study may have experienced an undocumented Omicron infection compared to the uninfected-controls in the BA.1-against-BA.2 study. This would bias to a lower value and may explain why was lower than that of .
Effectiveness against reinfection was estimated for only few weeks after the primary infection. A longer duration of follow-up may identify differences not yet seen given the recency of the Omicron wave. However, evidence has been consistent that natural immunity, unlike vaccine immunity, wanes slowly with minimal waning for at least several months after primary infection.2,5-9,20-22
BA.1 and BA.2 ascertainment was based on proxy criteria, presence or absence of SGTF using the TaqPath PCR assay, but this method of ascertainment is well established not only for Omicron sub-lineages, but also for other variants such as Alpha.6,19,33 BA.1 and BA.2 ascertainment was not possible for infections diagnosed using RAT or other PCR testing. This limitation was mitigated by basing the main analysis estimate on the distribution of known BA.1 and BA.2 cases for each calendar day, and by providing a sensitivity analysis where the estimate was made based on the distribution of known BA.1 and BA.2 cases in each cohort of the two studies.
Some Omicron infections may have been misclassified Delta infections, but this is not likely, as Delta incidence was limited during the time of follow-up (Section S1). With the recency of the Omicron wave, we had to use a short interval of 35 days to define reinfection, perhaps introducing bias due to misclassification of prolonged infections as reinfections. However, such potential bias is less likely to affect and , but may affect estimates of effectiveness of BA.1 (or BA.2) infection against reinfection with BA.1 (or BA.2); that is, when both the primary infection and the reinfection are both due to the same sub-lineage. Such effectiveness estimates are not reported in this study (but found in a separate analysis to be comparable to and ). Regardless, such bias leads to underestimation of effectiveness, as it would inflate incident cases only in the BA.1-infected or BA.2 infected cohorts, thereby further supporting our finding of strong protection against reinfection.
As an observational study, the investigated cohorts were neither blinded nor randomized, so unmeasured or uncontrolled confounding cannot be excluded. While matching was done for sex, age, and nationality, this was not possible for other factors, such as comorbidities, occupation, or geography, as such data were not available. However, matching was done to control for factors that affect infection exposure in Qatar.15,25-28 Matching by age may have reduced potential bias due to comorbidities. The number of individuals with severe chronic conditions is also small in Qatar’s young population.15,34 Matching by nationality may have partially controlled for differences in occupational risk or socio-economic status, given the association between nationality and occupation in Qatar.15,25-28 Qatar is essentially a city state and infection incidence was broadly distributed across neighborhoods/areas; that is, geography is not likely to have been a confounding factor. Lastly, matching by the considered factors has been shown to provide adequate control of bias in studies that used control groups in Qatar.16,35-38 These control groups included unvaccinated cohorts versus vaccinated cohorts within two weeks of the first dose,16,35-37 when vaccine protection is negligible,39,40 and mRNA-1273-versus BNT162b2-vaccinated cohorts, also in the first two weeks after the first dose.38 A strength of this study is exclusion of those with a documented prior infection before the Omicron wave, to minimize potential confounding introduced by natural immunity due to earlier variants.
In conclusion, infection with an Omicron sub-lineage appears to induce strong, but not full protection against reinfection with the other sub-lineage, for at least several weeks after the initial infection.
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.
Sources of support and 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.
Author contributions
HC co-designed the study, performed 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. HY, HAK, and MS conducted viral genome sequencing. PT and MRH conducted the multiplex, RT-qPCR variant screening and 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 1. Laboratory methods and variant ascertainment
Real-time reverse-transcription polymerase chain reaction testing
Nasopharyngeal and/or oropharyngeal swabs were collected for 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.
Rapid antigen testing
SARS-CoV-2 antigen tests were performed on nasopharyngeal swabs using one of the following lateral flow antigen tests: Panbio COVID-19 Ag Rapid Test Device (Abbott, USA); SARS-CoV-2 Rapid Antigen Test (Roche, Switzerland); Standard Q COVID-19 Antigen Test (SD Biosensor, Korea); or CareStart COVID-19 Antigen Test (Access Bio, USA). All antigen tests were performed point-of-care according to each manufacturer’s instructions at public or private hospitals and clinics throughout Qatar with prior authorization and training by the Ministry of Public Health (MOPH). Antigen test results were electronically reported to the MOPH in real time using the Antigen Test Management System which is integrated with the national COVID-19 database.
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 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 infections and 15 (4.8%) as Delta (B.1.617.2)9 infections.4,10 Of 286 Omicron infections with confirmed sub-lineage 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,11 The study duration coincided with the intense Omicron wave where Delta incidence was limited. Accordingly, any PCR or rapid antigen positive test during the study duration, between December 19, 2021 and February 21, 2022, was assumed to be an Omicron infection.
Informed by the viral genome sequencing and the RT-qPCR genotyping, a SARS-CoV-2 infection with the BA.1 sub-lineage was proxied as an S-gene “target failure” (SGTF) case using the TaqPath COVID-19 Combo Kits (Thermo Fisher Scientific, USA12).13 A SARS-CoV-2 infection with the BA.2 sub-lineage was proxied as a non-SGTF case using the TaqPath assay.