Study population and data sources
This study was conducted on Qatar’s resident population before and after the emergence of the omicron variant on 19 December 202116 (online supplemental section S1 and online supplemental figure S1). The first analysis estimated the effectiveness of a preomicron infection in preventing reinfection with a preomicron virus from 5 February 2020 (start of the COVID-19 pandemic in Qatar)26 until 18 December 2021 (end of the preomicron period). The second analysis estimated the effectiveness of an omicron infection in preventing reinfection with an omicron virus from 19 December 2021 (onset of the first omicron wave)16 until 12 March 2024 (end of study).
Throughout the study period, Qatar experienced the circulation of SARS-CoV-2 alongside other respiratory viruses, such as influenza, with the patterns of these viruses changing substantially over time due to the non-pharmaceutical interventions implemented during the COVID-19 pandemic.27 28
Data were sourced from the integrated, nationwide, digital health information platforms, which include the national federated databases for COVID-19 laboratory testing, vaccination, hospitalisation and death (online supplemental section S2). These databases contain SARS-CoV-2-related data on testing, with no missing information since the beginning of the pandemic, encompassing all PCR tests regardless of location or facility (online supplemental section S3). Starting 5 January 2022, these databases also incorporated all medically supervised rapid antigen (RA) tests (online supplemental section S3).
Until 31 October 2022, SARS-CoV-2 testing in Qatar was conducted on a mass scale, with around 5% of the population being tested every week, primarily for routine purposes such as screening or meeting travel-related requirements.29 30 Subsequently, testing rates decreased, with less than 1% of the population being tested per week.31 Most infections detected throughout the pandemic were identified through routine testing rather than due to the presence of symptoms.29 30
In December 2020, Qatar initiated its COVID-19 vaccination campaign using mRNA vaccines, prioritising individuals based on coexisting health conditions and age.29 32 Vaccination was tracked nationally and provided to all residents and citizens free of charge.29 32 Demographic information, such as sex, age and nationality, was obtained from the national health registry records. Qatar has a unique demographic profile, with only 9% of its population aged 50 years or older and 89% of its residents being expatriates from more than 150 countries.26 Detailed descriptions of Qatar’s population and national databases have been reported previously.26 29 30 33–36
Study design
A matched, test-negative, case–control study design was employed to assess the effectiveness of SARS-CoV-2 infection against reinfection.3 16 37–39 Cases and controls were identified based on SARS-CoV-2 testing within each analysis period. Individuals with positive tests were designated as cases, while those with negative tests were designated as controls. The odds of having a previous infection were then compared between cases and controls.3 16 37–39
The protection conferred by infection was estimated against five forms of reinfection: asymptomatic, symptomatic, severe24 COVID-19, critical24 COVID-19 and fatal25 COVID-19. SARS-CoV-2 reinfection is conventionally defined as a documented infection occurring at least 90 days after a previous infection to distinguish true reinfections from potential cases of prolonged viral shedding.40–42 Consistent with this definition, individuals with a positive SARS-CoV-2 test within the preceding 90 days were excluded from the analysis to avoid misclassifying cases of prolonged viral positivity as reinfections. The analysis was conducted over the entire study periods, as well as in 3-month intervals since the previous infection.
All PCR tests and a portion of the facility-based RA tests carried out in Qatar, irrespective of setting or location, are categorised based on symptoms and the reason for testing.29 30 These categories include clinical suspicion, contact tracing, surveys or random testing campaigns, individual request, routine healthcare testing, pretravel, post-travel or other. This categorisation allowed differentiation between tests conducted due to asymptomatic or symptomatic infections.
Asymptomatic infection was defined as a positive PCR or RA test result obtained when the reason for testing was a survey, with no symptoms compatible with a respiratory tract infection reported. Conversely, symptomatic infection was defined by a positive PCR or RA test prompted by the presence of symptoms consistent with a respiratory tract infection. Accordingly, only PCR or RA tests from surveys (asymptomatic) or those conducted due to clinical suspicion (symptomatic) were included in the respective analyses.
Classification of severe,24 critical24 and fatal25 COVID-19 followed WHO guidelines (online supplemental section S4). Trained medical personnel, independent of the study investigators, reviewed individual medical records to determine severity.43 As part of the national protocol, all individuals with a positive SARS-CoV-2 test and concurrent hospitalisation underwent infection severity assessments every 3 days until discharge or death, regardless of hospital stay duration or the time between the positive test and the final outcome.43
The severity assessment results were provided to the study investigators in the form of a database containing a categorical variable assessment. However, the database did not include detailed information used to determine severity, such as oxygen use and mechanical ventilation, which are available in the Cerner system that tracks all medical encounters within Qatar’s public healthcare system but to which investigators do not have direct access. The categorical variable assessment is considered complete, as it was conducted in accordance with a national COVID-19 protocol for clinical assessment implemented across the entire public healthcare system.
Individuals experiencing progression to severe, critical or fatal COVID-19 during the follow-up period were classified based on their worst outcome, with death25 being the highest severity, followed by critical24 illness and then severe24 illness.43
All cases and controls that met the inclusion criteria and could be matched were included in the analyses. The test-negative design helps mitigate potential confounding due to differential healthcare-seeking behaviour3 37–39 by restricting the population to individuals who sought testing for documented reasons. This approach ensures comparability between cases and controls regarding their healthcare-seeking patterns. Consequently, only tests with a documented reason for testing were included in the analyses.
The test-negative design applied in this study is an extension of the original test-negative design commonly used in influenza vaccine effectiveness studies.38 39 This extended design leverages the widespread testing conducted during the COVID-19 pandemic, where individuals were tested for various reasons, such as routine requirements, and has been thoroughly investigated and validated through numerous studies, sensitivity analyses, negative control outcomes and mathematical modelling simulations.3 29 32 37 44–46 Importantly, cases and controls are exactly matched by the reason for testing and the testing method in this design. This exact matching minimises bias stemming from differential testing or healthcare-seeking behaviour.
For analyses specific to the omicron era, individuals with a documented preomicron infection were excluded to separately investigate the effects of preomicron and omicron immunities. For the analyses of asymptomatic and symptomatic reinfections, cases and controls were matched on a one-to-one ratio. In contrast, to enhance statistical precision due to the lower number of cases with severe outcomes, a one-to-five matching ratio was employed for analyses of severe,24 critical24 and fatal25 COVID-19 reinfections.
Cases and controls were matched exactly on several factors to minimise the influence of confounding or mediating variables that might affect infection risk26 47–50 and to ensure non-differential healthcare-seeking behaviour.37–39 These factors, informed by our previous research in Qatar,22 29 32 44 45 51 included sex, 10-year age group, nationality, the number of coexisting conditions (0 to ≥6; online supplemental section S5), the number of vaccine doses received by the time of the study (outcome) test (0 to ≥4), calendar week of testing, testing method (PCR or RA) and reason for testing. This exact matching approach ensured that each case–control pair shared identical values for all these variables and controlled for the effect of vaccination on infection risk between cases and controls.
This study aimed to compare the protective effect of infection against various forms of reinfection. However, waning immunity can bias such comparisons if the time between the previous infection and the study test differs across analyses. To address this potential bias, matched case–control pairs identified through surveys were subsequently matched by calendar week of testing to their counterparts identified through clinical suspicion. Only case–control pairs matched on both the initial criteria and this additional calendar week matching were included in the analyses. Consequently, the distributions of time from previous infection to the study test were similar in both asymptomatic and symptomatic reinfection analyses.
Statistical analysis
All SARS-CoV-2 testing records were reviewed to select cases and controls; however, only matched samples were analysed. Cases and controls were described using frequency distributions and measures of central tendency and compared using standardised mean differences (SMDs). An SMD of ≤0.1 indicated adequate matching.52 The median and IQR for the time between the previous infection and the study test were calculated.
Conditional logistic regression was used to estimate ORs and their corresponding 95% CIs, comparing the odds of having a previous infection between cases and controls. The analysis considered the date of the most recent documented previous infection for each participant. Individuals with no documented previous infection served as the reference group for all comparisons. CIs were not adjusted for multiplicity, and interactions were not examined.
Based on the test-negative study design, effectiveness and corresponding 95% CIs were calculated as 1−OR of previous infection among cases versus controls if the OR was ≤137 and as (1/OR)−1 if the OR was >1.35 53 This method ensures a symmetric scale for both negative and positive effectiveness, ranging from −100% to 100%.35 53 In instances where conditional logistic regression failed to converge due to zero events among exposed cases, the 95% CIs were obtained using McNemar’s test. This approach provides only an approximate estimate of the CIs in these specific situations, following a method used in earlier studies.23 54
The effectiveness by time since previous infection was also assessed in 3-month intervals by restricting the study samples to cases and controls with a previous infection falling within the specified 3-month time-interval categories. Subgroup analyses were also performed, restricting the study sample to vaccinated and unvaccinated individuals separately. Statistical analyses were conducted using STATA/SE software V.18.0 (StataCorp).