Summary
Background Waning of natural infection protection and vaccine protection highlight the need to evaluate changes in population immunity over time. Population immunity of previous SARS-CoV-2 infection or of COVID-19 vaccination are defined, respectively, as the overall protection against reinfection or against breakthrough infection at a given point in time in a given population.
Methods We estimated these population immunities in Qatar’s population between July 1, 2020 and November 30, 2022, to discern generic features of the epidemiology of SARS-CoV-2. Effectiveness of previous infection, mRNA primary-series vaccination, and mRNA booster (third-dose) vaccination in preventing infection were estimated, month by month, using matched, test-negative, case-control studies.
Findings Previous-infection effectiveness against reinfection was strong before emergence of Omicron, but declined with time after a wave and rebounded after a new wave. Effectiveness dropped immediately after Omicron emergence from 88.3% (95% CI: 84.8-91.0%) in November 2021 to 51.0% (95% CI: 48.3-53.6%) in December 2021. Primary-series effectiveness against infection was 84.0% (95% CI: 83.0-85.0%) in April 2021, soon after introduction of vaccination, before waning gradually to 52.7% (95% CI: 46.5-58.2%) by November of 2021. Effectiveness declined linearly by ∼1 percentage point every 5 days. After Omicron emergence, effectiveness dropped suddenly from 52.7% (95% CI: 46.5-58.2%) in November 2021 to negligible levels in December 2021. Booster effectiveness dropped immediately after Omicron emergence from 83.0% (95% CI: 65.6 -91.6%) in November 2021 to 32.9% (95% CI: 26.7-38.5%) in December 2021, and continued to decline thereafter. Effectiveness of previous infection and vaccination against severe, critical, or fatal COVID-19 were generally >80% throughout the study duration.
Interpretation High population immunity may not be sustained beyond a year. This creates fertile grounds for repeated waves of infection to occur, but these waves may increasingly exhibit a benign pattern of infection.
Funding The Biomedical Research Program and the Biostatistics, Epidemiology, and the Biomathematics Research Core, both at Weill Cornell Medicine-Qatar, Ministry of Public Health, Hamad Medical Corporation, Sidra Medicine, Qatar Genome Programme, Qatar University Biomedical Research Center, and Qatar University Internal Grant ID QUCG-CAS-23/24-114.
Evidence before this study SARS-CoV-2 infection induces protection against reinfection, but this protection wanes with time since last infection. Similarly, COVID-19 primary-series and booster vaccination induce protection against SARS-CoV-2 infection, but this protection also wanes with time since last dose. These immunity patterns demonstrate the need for the concept of population immunity to track evolution of overall immune protection over time in a given population. Previous-infection and vaccine population immunities in a specific country can be defined as the overall protection against infection at a given point in time in the full national population. A search of PubMed, Google Scholar, and the International Vaccine Access Center’s VIEW-hub databases up to April 21, 2023 using the keywords “vaccination”, “infection”, “immunity”, “protection”, “SARS-CoV-2”, and “COVID-19” did not identify studies that investigated this epidemiological concept for a national population throughout the COVID-19 pandemic.
Added value of this study This study analyzed the national federated databases for SARS-CoV-2 infection and COVID-19 vaccination in Qatar, a country that experienced SARS-CoV-2 waves dominated by different pre-Omicron variants and Omicron subvariants. Using a matched, test-negative study design, population immunity against infection of each of previous infection, primary-series vaccination, and booster vaccination were characterized at the national level month by month for two calendar years to discern generic features of the epidemiology of SARS-CoV-2. The three forms of population immunity showed rapid variation over time driven by waning of protection. Vaccine-derived population immunity declined by 1 absolute percentage point every 5 days. Omicron introduction immensely reduced the three forms of population immunity within one month by about 50 absolute percentage points. Meanwhile, previous-infection and vaccine population immunities against severe COVID-19 were durable with slow waning even after Omicron emergence.
Implications of all available evidence Both previous-infection and vaccine population immunities vary rapidly at a national level creating fertile grounds for repeated waves of infection to occur even within months of each other. High levels of population immunity may not be sustained for more than a year or so. Preventing infection/reinfection, transmission, or future waves of infection cannot sustainably be done with current vaccines nor by the entire population being infected. Timely administration of boosters for those vulnerable to severe COVID-19 may remain essential for years to come. Repeated waves of infection may also facilitate further evolution of the virus and continual immune evasion. Emergence of a new variant that is substantially different from circulating variants can suddenly and immensely reduce population immunity leading to large epidemic waves. However, the durability of population immunity against severe COVID-19 will likely curtail the severity of future waves.
Introduction
Although immune protection of the primary series of coronavirus disease 2019 (COVID-19) vaccines is high against severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) infection immediately after the second dose,1, 2 protection wanes with time and may not last for more than 1 year after the second dose.3–5 The emergence of the immune-evasive Omicron subvariants reduced vaccine effectiveness immediately after the second dose to only about 50% and accelerated the waning of protection.6–8 Vaccine protection against Omicron subvariants may not last for more than 6 months after the second dose.6–8 Booster vaccination restores vaccine protection to the level observed immediately after the second dose,8, 9 but this protection also wanes with time and in a similar pattern to that of the primary series.6–10
Though protection of natural infection against reinfection is high immediately after infection,11, 12 the protection wanes with time and is not expected to last for more than 3 years.13 Introduction of Omicron subvariants reduced pre-Omicron infection protection against Omicron reinfection to only about 50%,14, 15 and accelerated the waning of this protection.13 Protection of a pre-Omicron infection against Omicron reinfection may not last for more than 1 year.13, 16 While protection of Omicron infection against Omicron reinfection is high,16–18 evidence suggests that this protection is declining fairly rapidly due to combined effect of waning of protection over time and progressive immune evasiveness of Omicron subvariants.16
These immune protection patterns demonstrate the need for the concept of population immunity to track evolution of population-level immune protection over time in a given population. Population immunity of primary-series vaccination is defined as the overall protection against infection at a given point in time among persons with only primary-series vaccination relative to those unvaccinated. Booster and previous-infection population immunities are defined analogously. These three forms of population immunity are expected to vary fairly rapidly even within a time horizon of only 1 year due to the combined effects of waning of vaccine and natural immunity and viral immune evasion. Such rapid variation may create fertile grounds for repeated waves of infection to occur even within few months of each other.
We aimed to characterize variation of these three forms of population immunity over time, month by month, in Qatar, a country that experienced SARS-CoV-2 waves dominated sequentially by the index virus,19 Alpha,20 Beta,21 Omicron BA.1 and BA.2,14 Omicron BA.4 and BA.5,18 Omicron BA.2.75*,16 and currently Omicron XBB*, in addition to a prolonged low-incidence phase dominated by Delta22 (Figure 1).
Methods
Study population and data sources
This study was conducted on the population of Qatar including data between July 1, 2020 and November 30, 2022. It analyzed the national, federated databases for COVID-19 laboratory testing, vaccination, hospitalization, and death, retrieved from the integrated, nationwide, digital-health information platform (Section S1 of the Supplementary Appendix). Databases include all SARS-CoV-2-related data with no missing information since the onset of the pandemic, including all polymerase chain reaction (PCR) tests, and from January 5, 2022 onward, all medically supervised rapid antigen tests. SARS-CoV-2 testing in Qatar was done at large scale, and up to October 31, 2022, was mostly done for routine reasons such as for screening or travel-related purposes, with infections primarily diagnosed not because of appearance of symptoms, but because of routine testing.3, 14 Qatar launched its COVID-19 vaccination program in December of 2020 using almost exclusively the BNT162b2 and mRNA-1273 mRNAvaccines.23 Detailed descriptions of Qatar’s population and of the national databases have been reported previously.3, 9, 10, 14, 19, 24
Study design
Effectiveness of previous infection, two-dose (primary-series) vaccination, and third-dose (booster) vaccination in preventing infection, regardless of symptoms, were estimated for Qatar’s population by calendar month. The calendar time over which these effectiveness measures were estimated depended on availability of primary-series and booster vaccination in the population and on having sufficient case numbers to generate estimates.
Each of the effectiveness measures was estimated using the test-negative, case-control, study design, a standard design for assessing immune protection of natural infection14, 15, 25 and vaccination.3, 26–28 In this design, effectiveness estimates are derived by comparing odds of previous infection or vaccination among positive tests (cases) compared to negative tests (controls).25–27
Effectiveness of previous infection in preventing reinfection was defined as the proportional reduction in susceptibility to infection among those with a previous infection versus those without.12, 14, 25 Effectiveness of primary-series (or booster) vaccination in preventing infection was defined as the proportional reduction in susceptibility to infection among those vaccinated versus those unvaccinated.3, 26–28
All cases (SARS-CoV-2-positive tests) and controls (SARS-CoV-2-negative tests) identified in Qatar during each calendar month of the study duration were matched exactly one-to-one by sex, 10-year age group, nationality, number of coexisting conditions (none, 1, 2, or ≥3), method of testing (PCR or rapid antigen), reason for SARS-CoV-2 testing, and (by design) calendar month of testing. Additionally, cases and controls were matched by COVID-19 vaccine type and number of vaccine doses at time of SARS-CoV-2 test in all analyses of effectiveness of previous infection, and by status of most recent prior infection (no documented prior infection, documented pre-Omicron prior infection, or documented Omicron prior infection) in all vaccine effectiveness analyses.
Matching was done to balance observed confounders between exposure groups that are related to risk of infection in Qatar.19, 29–32 Matching by the considered factors was informed by results of prior studies that used matching to control for differences in infection risk in Qatar, including test-negative, case-control studies.2–4, 23, 33 For each calendar month, only the first SARS-CoV-2-positive test for each case and the first SARS-CoV-2-negative test for each control were included. Persons qualified as controls if they had at least one SARS-CoV-2-negative test and no record of a SARS-CoV-2-positive test in the considered calendar month.
All persons who received a vaccine other than BNT162b2 or mRNA-1273, or who received a different mix of vaccines, were excluded. Individuals who received a fourth vaccine dose (second booster dose) prior to the study SARS-CoV-2 test were excluded. Bivalent boosters have only been recently introduced in Qatar and remain at low coverage, thus they were not included in this study. Tests occurring within 14 days of a second dose or 7 days of a third dose were excluded. SARS-CoV-2 tests conducted because of travel-related requirements were excluded.3, 4 These inclusion and exclusion criteria were implemented to allow for build-up of immunity after vaccination,1, 9 and to minimize different types of potential bias, as informed by earlier analyses on the same population.3, 4 Every case (or control) that met the inclusion criteria and that could be matched to a control (case) was included in the analysis. Previous infection status and COVID-19 vaccination status were ascertained at the time of the SARS-CoV-2 test.
SARS-CoV-2 reinfection is conventionally defined as a documented infection ≥90 days after an earlier infection, to avoid misclassification of prolonged SARS-CoV-2 positivity as reinfection, if a shorter time interval is used.34, 35 Previous infection was thus defined as a SARS-CoV-2-positive test ≥90 days before the study’s SARS-CoV-2 test.14, 15, 25 Cases or controls with SARS-CoV-2-positive tests <90 days before the study’s SARS-CoV-2 test were excluded.
Effectiveness was also estimated against severe, critical, or fatal COVID-19 due to SARS-CoV-2 infection using the same methodology, but for 6-calendar-month durations owing to the small number of cases with severe forms of COVID-19. Cases and controls were matched one-to-five to increase precision of estimates. Classification of COVID-19 case severity (acute-care hospitalizations),36 criticality (intensive-care-unit hospitalizations),36 and fatality37 followed World Health Organization (WHO) guidelines, and assessments were made by trained medical personnel using individual chart reviews (Section S3).
Each person who had a SARS-CoV-2-positive test and COVID-19 hospital admission was subject to an infection severity assessment every three days until discharge or death, regardless of the length of hospital stay or the time between the SARS-CoV-2-positive test and the final disease outcome. Individuals who progressed to severe,36 critical,36 or fatal37 COVID-19 between the SARS-CoV-2-positive test and the end of study were classified based on their worst outcome, starting with death, followed by critical disease, and then severe disease.
Oversight
The institutional review boards at Hamad Medical Corporation and Weill Cornell Medicine– Qatar approved this retrospective study with a waiver of informed consent. The study was reported according to the Strengthening the Reporting of Observational Studies in Epidemiology (STROBE) guidelines (Table S1).
Statistical analysis
While all records of SARS-CoV-2 testing were examined for selection of cases and controls, only matched samples were analyzed. Cases and controls were described using frequency distributions and measures of central tendency and compared using standardized mean differences (SMDs). An SMD of ≤0.1 indicated adequate matching.38 The median and interquartile range (IQR) of the duration between the immunological event (previous infection, primary-series vaccination, or booster vaccination) and SARS-CoV-2 test were calculated for cases and controls in each calendar-month analysis.
The odds ratio (and its associated 95% confidence interval (CI)), comparing odds of the immunological event among cases to that among controls, was estimated using conditional logistic regression, that is factoring the matching in the study design. This analytical approach was implemented to reduce potential bias due to variation in epidemic phase,26, 39 gradual vaccination roll-out,26, 39 and other confounders.19, 40, 41 CIs did not factor multiplicity and interactions were not examined.
Vaccine effectiveness and associated 95% CI were estimated as 1-odds ratio (OR) of the immunological event among cases versus controls if the OR was <1,25–27 and as (1/OR)-1 if the OR was ≥1.10, 42 The latter was done to ensure symmetric scale for both negative and positive effectiveness, ranging from -100%-100%.10, 42 Statistical analyses were conducted in STATA/SE version 17.0 (Stata Corporation, College Station, TX, USA).
Role of the funding source
The funders of the study had no role in study design, data collection, data analysis, data interpretation, or writing of the report. The corresponding author had full access to all the data in the study and had final responsibility for the decision to submit for publication.
Results
Effectiveness of previous infection against reinfection
Characteristics of matched cases and controls of the study population are shown in Table S2. Effectiveness of previous infection against reinfection was high from July of 2020 to November of 2021, a duration coinciding with incidence of pre-Omicron variants (Figure 1), at a level that exceeded 70% (Table S3 and Figure 2A). Effectiveness declined slowly over time after a wave, but rebounded to a higher level after a new wave. Note that since cases and controls with SARS-CoV-2-positive tests <90 days before the study’s SARS-CoV-2 test were excluded to avoid misclassification of prolonged SARS-CoV-2 positivity as reinfection,34, 35 estimates of the effectiveness of previous infections represent infections that occurred ≥90 days prior and not that of recent infections.
Effectiveness dropped immediately and immensely with the emergence of Omicron in December 2021 (Table S3 and Figure 2A). While effectiveness was 88.3% (95% CI: 84.8-91.0%) in November 2021, it dropped to 51.0% (95% CI: 48.3-53.6%) in December 2021; a drop of ∼40 absolute percentage points within one month associated with the change in the circulating variant. Effectiveness continued at this level during the large BA.1/BA.2 Omicron wave, up to March 2022.
Effectiveness rebounded again after end of the BA.1/BA.2 wave and reached 80.8% (95% CI: 77.6-83.5%) in May of 2022 before gradually declining in subsequent months (Table S3 and Figure 2A). Effectiveness was low at ∼30% by end of the study in November of 2022, a time coinciding with incidence of BA.2.75* and XBB* subvariants (Figure 1).
Effectiveness of primary-series vaccination against infection
Characteristics of matched cases and controls of the study population are shown in Table S2. Effectiveness of primary-series vaccination against infection was very high at >80% right after introduction of vaccination (Table S4 and Figure 2B). Effectiveness was 84.5% (95% CI: 77.8-89.2%) in February of 2021 and remained high up to June of 2021, coinciding with the rapid scale-up of vaccination in Qatar3 and incidence of pre-Omicron variants (Figure 1).
However, starting from July of 2021, when the mass vaccination campaigns slowed down along with a decline in the recentness of the second dose for much of the population, effectiveness started to decline rapidly (Table S4 and Figure 2B). Effectiveness was 80.7% (95% CI: 77.2-83.7%) in June of 2021 before dropping to 76.2% (95% CI: 71.8-79.9%) in July of 2021, 63.0% (95% CI: 57.8-67.5%) in August of 2021, 61.9% (95% CI: 54.7-68.0%) in September of 2021, 44.6% (95% CI: 32.8-54.4%) in October of 2021, and 52.7% (95% CI: 46.5-58.2%) in November of 2021, right before introduction of Omicron (Figure 1).
Effectiveness declined linearly with calendar time by ∼1 absolute percentage point every 5 days (Figure 3A). Effectiveness declined linearly with median time since the second dose also by ∼1 absolute percentage point every 5 days (Figure 3B). Extrapolating this linear trend indicated that effectiveness would reach 0% in 13.4 months after the second dose.
Effectiveness dropped suddenly and massively with introduction of Omicron in December of 2021 (Table S4 and Figure 2B). While effectiveness was 52.7% (95% CI: 46.5-58.2%) in November of 2021, it became negligible during the BA.1/BA.2 wave starting from December of 2021, that is a drop of ∼50 absolute percentage points within one month associated with the change in circulating variant. Effectiveness remained negligible in subsequent months and was notably negative during the BA.4/BA.5 wave in June and July of 2022.
Effectiveness of booster vaccination against infection
Characteristics of matched cases and controls of the study population are shown in Table S2. Effectiveness of booster (third dose) vaccination against infection was very high at 83.0% (95% CI: 65.6-91.6%) in November of 2021 (Table S5 and Figure 2C), when booster vaccination was being scaled up9 and Delta was circulating (Figure 1). But effectiveness dropped suddenly and immensely with the emergence of Omicron in December of 2021. Effectiveness dropped from 83.0% (95% CI: 65.6-91.6%) in November of 2021 to 32.9% (95% CI: 26.7-38.5%) in December of 2021, that is a drop of ∼50 absolute percentage points within one month associated with the change in the circulating variant. Effectiveness continued to decline thereafter and became negative between April and November of 2022, particularly during the BA.4/BA.5 wave.
Effectiveness against severe, critical, or fatal COVID-19
Unlike effectiveness against infection, effectiveness against severe, critical, or fatal COVID-19 of previous infection (Table S6 and Figure 4A), primary-series vaccination (Table S6 and Figure 4B), and booster vaccination (Table S6 and Figure 4C) were high throughout the study duration at generally >80%. There appeared to be some waning in effectiveness of primary-series vaccination and booster vaccination over time, but the number of severe COVID-19 cases was too small to allow precise estimation of potential waning in effectiveness.
Discussion
This study investigated three population immunity patterns of relevance to the future of this pandemic. First, both natural and vaccine population immunity wane with time, with the waning being particularly rapid for vaccine population immunity. Even before Omicron introduction, vaccine protection was waning by ∼1 absolute percentage point every 5 days. The protection induced by vaccination appeared to fully wane within 14 months. These findings suggest that repeated waves of infection will become a generic feature of SARS-CoV-2 epidemiology. High levels of population immunity may not be sustained for more than a year or so.
Second, sudden emergence of a new variant that is substantially different from previously circulating variants can immensely and immediately affect population immunity. Protection of previous infection dropped by ∼40 absolute percentage points within only one month after Omicron introduction. Similarly, each of primary-series and booster protections dropped by ∼50 absolute percentage points. As a result, Qatar experienced its largest SARS-CoV-2 epidemic wave to date (Figure 1). This demonstrates how the long-term epidemiology of this infection can be strongly sensitive to major shifts in virus evolution. Active surveillance of emerging variants and development of an early warning system are critical to mitigate the consequences of immune-evasive variants.43
Third, although population immunity against infection waned rapidly, population immunity against severe COVID-19 was durable and showed slow waning even with introduction of Omicron. Infection with common-cold coronaviruses,44 and perhaps influenza,45 induces only a year-long protection against infection, but life-long immunity against severe reinfection.46 While it is premature to make long-term predictions, this finding suggests that SARS-CoV-2 epidemiology may exhibit a similar pattern to that of common-cold coronaviruses. Long-term immune protection against severe COVID-19 could contribute to a benign pattern of infection that is perhaps not dissimilar to that of common-cold coronaviruses.
Some of the vaccine effectiveness measures post Omicron introduction, particularly for booster vaccination, were negative in value, perhaps suggesting negative immune imprinting. This effect was pronounced during the BA.4/BA.5 wave. This finding supports similar recent findings in this same population.10, 47 Such imprinting effects have been observed for other infections such as influenza.48, 49 It remains to be seen whether these effects will be of consequence in the future epidemiology of SARS-CoV-2 infection.
This study has limitations. Since SARS-CoV-2 reinfection is defined as an infection ≥90 days after an earlier infection, estimates of previous-infection protection lagged actual level of protection by 3 months. This has underestimated previous-infection population immunity during or right after large waves, when population immunity was building up rapidly. With the relatively young population of Qatar, our findings may not be generalizable to other countries where elderly citizens constitute a larger proportion of the total population.
At this stage of the pandemic, the reference groups of persons with no previous infection and unvaccinated persons may not be representative due to depletion of susceptibles and undocumented infections and vaccinations. Vaccinations outside of Qatar may not have been registered in the national database after easing vaccine requirements in 2021. Bias due to depletion of susceptibles may lead to underestimation of previous-infection or vaccine protections,50 even in the test-negative study design which is less prone to effect of this bias.3, 25 Misclassification of previous-infection status can lead to underestimation of previous-infection effectiveness particularly if the majority of the population had already been infected.25 These sources of bias may explain some of the negative vaccine effectiveness values observed in this study post Omicron introduction and the lower than expected previous-infection protection in the last few months of the study. Estimates after the first large Omicron wave may provide qualitative findings rather than precise quantitative estimates.
Bias can arise in real-world observational data in unexpected ways, or from unknown sources, such as subtle differences in test-seeking behavior or changes in the pattern of testing due to policy changes, tests’ accessibility, or behavioral differences (Section S1). The reason for testing varies with time even if there are no policy changes. During a wave, a large proportion of testing occurs because of clinical suspicion of infection. At times of low incidence, a large proportion of testing occurs because of routine reasons.
The study did not investigate population immunity of a second booster dose or of bivalent boosters due to the low coverage of these vaccinations in Qatar. In the first few weeks of the pandemic, PCR workflows were still under development and there was a chance of false-positive tests. While this may have affected only a very small number of PCR tests, reinfections were so rare at this time that few such false-positive cases could have resulted in underestimation of previous-infection protection in the first few months of the pandemic.
While matching was done for several factors, it was not possible for other factors such as geography or occupation, as such data were unavailable. However, Qatar is essentially a city state and infection incidence was broadly distributed across neighborhoods. Nationality, age, and sex provide a powerful proxy for socio-economic status in this country,19, 29–32 and thus matching by these factors may have partially controlled for other factors such as occupation. The matching prescription had already been investigated in previous studies of different epidemiologic designs, and using control groups to test for null effects.2–4, 23, 33 These studies have supported that this prescription provides adequate control of the differences in infection exposure.2–4, 23, 33 The study was implemented on Qatar’s total population, perhaps minimizing the likelihood of bias.
In conclusion, both previous-infection and vaccine population immunities wane with time and the waning is rapid for vaccine protection. High levels of population immunity may not be sustained for more than a year. The waning of population immunity facilitates fertile grounds for repeated waves of infection to occur even within few months of each other. Timely administration of boosters for those vulnerable to COVID-19 may remain essential for years to come. Emergence of a new variant that is substantially different from previously circulating variants can suddenly and immensely reduce population immunity leading to large epidemic waves. Repeated waves of infection may also facilitate further evolution of the virus and continual immune evasion. However, with the durability of population immunity against severe COVID-19, the repeated epidemic waves may increasingly exhibit an overall benign pattern of 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. The data are available under restricted access for preservation of confidentiality of patient data. Access can be obtained through a direct application for data access to Her Excellency the Minister of Public Health (https://www.moph.gov.qa/english/OurServices/eservices/Pages/Governmental-HealthCommunication-Center.aspx). The raw data are protected and are not available due to data privacy laws. Aggregate data are available within the paper and its supplementary information.
Contributors
SHQ, HC, and LJA co-designed the study. SHQ performed the statistical analyses and co-wrote the first draft of the article. LJA conceived the study, led the statistical analyses, and co-wrote the first draft of the article. PVC designed mass PCR testing to allow routine capture of SGTF variants and conducted viral genome sequencing. PT and MRH designed and conducted multiplex, RT-qPCR variant screening and viral genome sequencing. HY, AAA, and HAK 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.
Data sharing
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. The data are available under restricted access for preservation of confidentiality of patient data. Access can be obtained through a direct application for data access to Her Excellency the Minister of Public Health (https://www.moph.gov.qa/english/OurServices/eservices/Pages/Governmental-HealthCommunication-Center.aspx). The raw data are protected and are not available due to data privacy laws. Aggregate data are available within the paper and its supplementary information.
Declaration of 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: Further details on methods
Data sources and testing
Qatar’s national and universal public healthcare system uses the Cerner-system advanced digital health platform to track all electronic health record encounters of each individual in the country, including all citizens and residents registered in the national and universal public healthcare system. Registration in the public healthcare system is mandatory for citizens and residents.
The databases analyzed in this study are data-extract downloads from the Cerner-system that have been implemented on a regular (twice weekly) schedule since the onset of the pandemic by the Business Intelligence Unit at Hamad Medical Corporation. Hamad Medical Corporation (HMC) is the national public healthcare provider in Qatar. At every download all tests, coronavirus disease 2019 (COVID-19) vaccinations, hospitalizations related to COVID-19, and all death records regardless of cause are provided to the authors through .csv files. These databases have been analyzed throughout the pandemic not only for study-related purposes, but also to provide policymakers with summary data and analytics to inform the national response.
Every health encounter in the Cerner-system is linked to a unique individual through the HMC Number that links all records for this individual at the national level. Databases were merged and analyzed using the HMC Number to link all records whether for testing, vaccinations, hospitalizations, and deaths. All deaths in Qatar are tracked by the public healthcare system. All COVID-19-related healthcare was provided only in the public healthcare system. No private entity was permitted to provide COVID-19-related healthcare. COVID-19 vaccination was also provided only through the public healthcare system. These health records were tracked throughout the COVID-19 pandemic using the Cerner system. This system has been implemented in 2013, before the onset of the pandemic.
Demographic details for every HMC Number (individual) such as sex, age, and nationality are collected upon issuing of the universal health card, based on the Qatar Identity Card, which is a mandatory requirement by the Ministry of Interior to every citizen and resident in the country.
Severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) testing in the healthcare system in Qatar is done at a mass scale, and up to October 31, 2022, was mostly done for routine reasons, where about 5% of the population were tested every week.1, 2 All SARS-CoV-2 testing in any facility in this country is tracked nationally in one database, the national testing database. This database covers all testing in all locations and facilities throughout the country, whether public or private. Every polymerase chain reaction (PCR) test and an increasing proportion of the medically supervised rapid antigen tests conducted in Qatar, regardless of location or setting, are classified on the basis of symptoms and 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). Based on the distribution of the reason for testing up to October 31, 2022, most of the tests that have been conducted in Qatar were conducted for routine reasons, such as being travel-related. About 75% of those diagnosed are also diagnosed not because of appearance of symptoms, but because of routine testing.1, 2
The first large Omicron wave that peaked in January of 2022 was massive and strained the testing capacity in the country.1, 3 Accordingly, rapid antigen testing was introduced to relieve the pressure on PCR testing. Implementation of this change in testing occurred quickly precluding incorporation of reason for testing in a large proportion of the rapid antigen tests for several months. While the reason for testing is available for all PCR tests, it is not available for all rapid antigen tests. Availability of reason for testing for the rapid antigen tests also varied with time. Rapid antigen test kits are available for purchase in pharmacies in Qatar, but outcome of home-based testing is not reported nor documented in the national databases. Since SARS-CoV-2-test outcomes are linked to specific public health measures, restrictions, and privileges, testing policy and guidelines stress facility-based testing as the core testing mechanism in the population. While facility-based testing is provided free of charge or at low subsidized costs, depending on the reason for testing, home-based rapid antigen testing is de-emphasized and not supported as part of national policy.
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.4, 5 Further descriptions of the study population and these national databases were reported previously.1, 2, 5–7
Comorbidity classification
Comorbidities were ascertained and classified based on the ICD-10 codes as recorded in the electronic health record encounters of each individual in the Cerner-system national database that includes all citizens and residents registered in the national and universal public healthcare system. The public healthcare system provides healthcare to the entire resident population of Qatar free of charge or at heavily subsidized costs, including prescription drugs.
All encounters for each individual were analyzed to determine the comorbidity classification for that individual, as part of a recent national analysis to assess healthcare needs and resource allocation. The Cerner-system national database includes encounters starting from 2013, after this system was launched in Qatar. As long as each individual had at least one encounter with a specific comorbidity diagnosis since 2013, this person was classified with this comorbidity. Individuals who have comorbidities but never sought care in the public healthcare system, or seek care exclusively in private healthcare facilities, were classified as individuals with no comorbidity due to absence of recorded encounters for them.
Section S2: 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.
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 Coronavirus Disease 2019 (COVID-19) database.
Classification of infections by variant type
Surveillance for SARS-CoV-2 variants in Qatar is based on viral genome sequencing and multiplex RT-qPCR variant screening8 of random positive clinical samples,2, 9–13 complemented by deep sequencing of wastewater samples.11, 14, 15 Further details on the viral genome sequencing and multiplex RT-qPCR variant screening throughout the SARS-CoV-2 waves in Qatar can be found in previous publications.1–3, 7, 9–13, 16–20
Section S3: COVID-19 severity, criticality, and fatality classification
Classification of COVID-19 case severity (acute-care hospitalizations),21 criticality (intensive-care-unit hospitalizations),21 and fatality22 followed World Health Organization (WHO) guidelines. Assessments were made by trained medical personnel independent of study investigators and using individual chart reviews, as part of a national protocol applied to every hospitalized COVID-19 patient. Each hospitalized COVID-19 patient underwent an infection severity assessment every three days until discharge or death. We classified individuals who progressed to severe, critical, or fatal COVID-19 between the time of the documented infection and the end of the study based on their worst outcome, starting with death,22 followed by critical disease,21 and then severe disease.21
Severe COVID-19 disease was defined per WHO classification as a SARS-CoV-2 infected person with “oxygen saturation of <90% on room air, and/or respiratory rate of >30 breaths/minute in adults and children >5 years old (or ≥60 breaths/minute in children <2 months old or ≥50 breaths/minute in children 2-11 months old or ≥40 breaths/minute in children 1–5 years old), and/or signs of severe respiratory distress (accessory muscle use and inability to complete full sentences, and, in children, very severe chest wall indrawing, grunting, central cyanosis, or presence of any other general danger signs)”.21 Detailed WHO criteria for classifying SARS-CoV-2 infection severity can be found in the WHO technical report.21
Critical COVID-19 disease was defined per WHO classification as a SARS-CoV-2 infected person with “acute respiratory distress syndrome, sepsis, septic shock, or other conditions that would normally require the provision of life sustaining therapies such as mechanical ventilation (invasive or non-invasive) or vasopressor therapy”.21 Detailed WHO criteria for classifying SARS-CoV-2 infection criticality can be found in the WHO technical report.21
COVID-19 death was defined per WHO classification as “a death resulting from a clinically compatible illness, in a probable or confirmed COVID-19 case, unless there is a clear alternative cause of death that cannot be related to COVID-19 disease (e.g. trauma). There should be no period of complete recovery from COVID-19 between illness and death. A death due to COVID-19 may not be attributed to another disease (e.g. cancer) and should be counted independently of preexisting conditions that are suspected of triggering a severe course of COVID-19”. Detailed WHO criteria for classifying COVID-19 death can be found in the WHO technical report.22
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. HHA acknowledges the support of Qatar University Internal Grant ID QUCG-CAS-23/24-114. 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.