Abstract
Effect of prior SARS-CoV-2 infection on vaccine protection remains poorly understood. Here, we investigated whether persons vaccinated after a prior infection have better protection against future infection than those vaccinated without prior infection. Effect of prior infection was assessed in Qatar’s population, where the Alpha (B.1.1.7) and Beta (B.1.351) variants dominate incidence, using two national retrospective, matched-cohort studies, one for the BNT162b2 (Pfizer-BioNTech) vaccine, and one for the mRNA-1273 (Moderna) vaccine. Incidence rates of infection among BNT162b2-vaccinated persons, with and without prior infection, were estimated, respectively, at 1.66 (95% CI: 1.26-2.18) and 11.02 (95% CI: 9.90-12.26) per 10,000 person-weeks. The incidence rate ratio was 0.15 (95% CI: 0.11-0.20). Analogous incidence rates among mRNA-1273-vaccinated persons were estimated at 1.55 (95% CI: 0.86-2.80) and 1.83 (95% CI: 1.07-3.16) per 10,000 person-weeks. The incidence rate ratio was 0.85 (95% CI: 0.34-2.05). Prior infection enhanced protection of those BNT162b2-vaccinated, but not those mRNA-1273-vaccinated. These findings may have implications for dosing, interval between doses, and potential need for booster vaccination.
Main text
Effect of prior acute respiratory syndrome coronavirus 2 (SARS-CoV-2) infection on vaccine protection against acquisition of infection remains poorly understood1-3. Qatar launched Coronavirus Disease 2019 (COVID-19) immunization in December 21, 2020, first using the BNT162b24 (Pfizer-BioNTech) vaccine and subsequently adding the mRNA-12735 (Moderna) vaccine6,7. As vaccination was scaled up following the FDA-approved protocol, the country experienced two back-to-back SARS-CoV-2 waves from January-June, 2021, which were dominated by the Alpha8 (B.1.1.7) and Beta8 (B.1.351) variants6,7,9-11 (Methods). This provided an opportunity to assess whether persons vaccinated after a prior SARS-CoV-2 infection have better protection against future infection than those vaccinated without prior infection.
Leveraging the national, federated databases that have captured all SARS-CoV-2 vaccinations and PCR testing since the epidemic onset (Methods), we investigated this question using two retrospective, matched-cohort studies. We compared incidence of documented SARS-CoV-2 infection in the national cohort of individuals who completed ≥14 days after the second BNT162b2 vaccine dose, but who had experienced a prior PCR-confirmed infection, with incidence among individuals who completed ≥14 days after the second BNT162b2 dose, but who had not experienced a prior infection, between December 21, 2020-June 6, 2021 (Figure 1). The same comparison was made for the mRNA-1273 vaccine (Figure 2). Cohorts were matched in a 1:1 ratio by sex, 5-year age group, nationality, and calendar week of the first vaccine dose, to control for differences in exposure risk12,13 and variant exposure6,7,9-11. Reporting of the study followed the STROBE guidelines (Supplementary Table 1).
The process for identifying SARS-CoV-2 infections in the national cohort of individuals who completed ≥14 days after the second BNT162b2 vaccine dose and who had experienced a PCR-confirmed infection before the first dose, compared with the process for identifying SARS-CoV-2 infections in the national cohort of individuals who completed ≥14 days after the second BNT162b2 vaccine dose, but who had experienced no PCR-confirmed infection before the first dose. Cohorts were matched in a 1:1 ratio by sex, 5-year age group, nationality, and calendar week of the first vaccine dose. Total follow-up time among BNT162b2-vaccinated persons, with and without prior infection, was 308,086.0 and 305,891.9 person-weeks, respectively.
The process for identifying SARS-CoV-2 infections in the national cohort of individuals who completed ≥14 days after the second mRNA-1273 vaccine dose and who had experienced a PCR-confirmed infection before the first dose, compared with the process for identifying SARS-CoV-2 infections in the national cohort of individuals who completed ≥14 days after the second mRNA-1273 vaccine dose, but who had experienced no PCR-confirmed infection before the first dose. Cohorts were matched in a 1:1 ratio by sex, 5-year age group, nationality, and calendar week of the first vaccine dose. Total follow-up time among mRNA-1273-vaccinated persons, with and without prior infection, was 70,729.9 and 70,872 person-weeks, respectively.
Kaplan-Meier curves showing the cumulative incidence of documented SARS-CoV-2 infection in the national cohort of individuals who completed ≥14 days after the second vaccine dose and who had a prior PCR-confirmed infection, compared to the cumulative incidence of documented SARS-CoV-2 infection in the matched national cohort of individuals who completed ≥14 days after the second vaccine dose, but without prior PCR-confirmed infection. The curves compare vaccination with A) the BNT162b2 (Pfizer-BioNTech) vaccine and B) the mRNA-1273 vaccine. Cohorts were matched in a 1:1 ratio by sex, 5-year age group, nationality, and calendar week of the first vaccine dose. The curves for a longer time of follow up for only the BNT162b2 vaccine are in Supplementary Figure 1. Vaccination with BNT162b2 started few weeks before vaccination with mRNA-1273.
Figures 1-2 show the process for identifying infections in these cohorts, and Table 1 presents their demographic characteristics. Using the Kaplan–Meier estimator14, cumulative infection incidence among BNT162b2-vaccinated persons, with and without prior infection, was estimated at 0.14% (95% CI: 0.11-0.19%) and 0.93% (95% CI: 0.83-1.04%), respectively, after 63 days of follow-up (Figure 1). Incidence rates of infection were estimated, respectively, at 1.66 (95% CI: 1.26-2.18) and 11.02 (95% CI: 9.90-12.26) per 10,000 person-weeks. The incidence rate ratio was estimated at 0.15 (95% CI: 0.11-0.20).
Cumulative infection incidence among mRNA-1273-vaccinated persons, with and without prior infection, was estimated at 0.06% (95% CI: 0.03-0.12%) and 0.08% (95% CI: 0.04-0.15%), respectively, after 63 days of follow-up (Figure 1). Incidence rates were estimated, respectively, at 1.55 (95% CI: 0.86-2.80) and 1.83 (95% CI: 1.07-3.16) per 10,000 person-weeks. The incidence rate ratio was estimated at 0.85 (95% CI: 0.34-2.05).
Infection incidence was low in these cohorts during a time of intense incidence in Qatar6,7,15, indicating that both vaccines were highly effective against the Alpha and Beta variants6,7, which dominated incidence9 (Methods). Still, prior infection of those BNT162b2-vaccinated further enhanced protection and reduced the incidence rate by 85% (6.6-fold) compared to those without prior infection. No evidence for such an effect was found for those mRNA-1273-vaccinated.
These findings are perhaps explained by the observed differences in effectiveness of these two vaccines against the Alpha and Beta variants, estimated in Qatar at 89.5% (95% CI: 85.9-92.3%) and 75.0% (95% CI: 70.5-78.9%) for BNT162b2, respectively6, and at 100% (95% CI: 91.8-100.0%) and 96.4% (95% CI: 91.9-98.7%) for mRNA-1273, respectively7.
The differences in effectiveness could have risen for a variety of reasons, such as differences in dosing, interval between doses, or the biology of both vaccines and their mechanisms of action. The dose of each of these two vaccines differed—it was 30-μg per dose for BNT162b24 and 100 μg per dose for mRNA-12735. This may have resulted in a more activated immune response for the mRNA-1273 vaccine than the BNT162b2 vaccine, and made the existence of prior immunity due to natural infection of no additional benefit for the mRNA-1273 vaccine. The interval between doses also differed and was one week longer for mRNA-12735. Evidence suggests that a longer dose interval could be associated with improved protection after receiving the second dose16.
Limitations include identifying prior infection based on a record of a PCR-positive result, thereby missing those who may have been infected, but were unaware of their infection, or who did not seek testing by PCR to document the infection. Misclassification of prior infection status could lead to underestimation of the effect size of prior infection on vaccine protection. Depletion of the cohorts with prior infection due to COVID-19 mortality at time of the prior infection may have biased these cohorts toward healthier individuals with stronger immune responses. However, COVID-19 mortality has been low in Qatar’s predominantly young and working-age population12,17, and no evidence for such bias was found in the mRNA-1273 vaccine results, where the incidence rate was similar for those with and without prior infection.
We assessed risk of only documented infections, but other infections may have occurred and gone undocumented, perhaps because of minimal/mild or no symptoms. Our cohorts predominantly included working-age adults; therefore, results may not necessarily be generalizable to other population groups, such as children or the elderly. Matching was done for age, sex, nationality, and calendar week of the first vaccine dose, and could not be done for other factors, such as comorbidities or additional socio-demographic factors, as these were not available to study investigators. However, matching by age and sex may have served as a proxy given that co-morbidities are associated with older age and may be different between women and men. Matching by nationality may have also captured some of the occupational risk given the distribution of the labor force in Qatar18-20.
Imperfect assay sensitivity and specificity of PCR or antibody testing could have affected current or prior infection ascertainment. However, all PCR and serological testing was performed with extensively used, investigated, and validated commercial platforms with essentially 100% sensitivity and specificity (Methods). Unlike blinded, randomized clinical trials, the investigated observational cohorts were neither blinded nor randomized.
Our results demonstrate low infection incidence among those vaccinated with BNT162b2 or mRNA-1273, but among those vaccinated with BNT162b2, protection against infection was further enhanced and infection incidence was further reduced by prior infection. In contrast, those vaccinated with mRNA-1273 were as well protected as those who received the vaccine after a prior infection. These findings may have implications for the potential need of a booster vaccination.
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
LJA conceived and co-designed the study, led the statistical analyses, and co-wrote the first draft of the article. HC co-designed the study, performed the statistical analyses, and co-wrote the first draft of the article. All authors contributed to data collection and acquisition, database development, discussion and interpretation of the results, and to the writing of the manuscript.
All authors have read and approved the final manuscript.
Competing interests
Dr. Butt has received institutional grant funding from Gilead Sciences unrelated to the work presented in this paper. Otherwise, we declare no competing interests.
Methods
Data sources and study design
Analyses were conducted using the centralized, integrated, and standardized national severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) databases compiled at Hamad Medical Corporation (HMC), the main public healthcare provider and the nationally designated provider for all Coronavirus Disease 2019 (COVID-19) healthcare needs. Through a nation-wide digital health information platform, these databases have captured all SARS-CoV-2-related data along with related-demographic details with no missing information since the start of the epidemic, including all records of polymerase chain reaction (PCR) testing, antibody testing, COVID-19 hospitalizations, vaccinations, infection severity classification per World Health Organization (WHO) guidelines21 (performed by trained medical personnel through individual chart reviews), and COVID-19 deaths, also assessed per WHO guidelines22. Every PCR test conducted in Qatar, regardless of location (outpatient clinic, drive-thru, or hospital, etc.), is classified on the basis of symptoms and the reason for testing (clinical symptoms, contact tracing, random testing campaigns (surveys), individual requests, routine healthcare testing, pre-travel, and port of entry). Qatar has unique demographics by sex and nationality, since expatriates from over 150 countries comprise 89% of the population12,23.
The nature of circulating SARS-CoV-2 virus was informed by weekly rounds of viral genome sequencing and multiplex, quantitative, reverse-transcription PCR (RT-qPCR) variant screening24 of randomly collected clinical samples6,7,9-11, as well as by the results of deep sequencing of wastewater samples9. The weekly rounds of viral genome sequencing from January 1-May 19, 2021 identified Beta (n=623; 50.9%), Alpha (n=193; 15.8%), Delta (n=43; 3.5%), and wild-type/undetermined variants (n=366; 29.9%) in 1,225 randomly collected, PCR-positive specimens9,10. Meanwhile, the weekly rounds of multiplex RT-qPCR variant screening from March 23-May 10, 2021 identified Beta-like (n=2,605; 66.4%), Alpha-like (n=970; 24.7%), and “other” variants (n=349; 8.9%) in 3,924 randomly collected PCR-positive specimens9,11. Sanger sequencing of the receptor binding domain of SARS-CoV-2 spike protein on 109 “other” specimens confirmed that 103 were Delta-like, 3 were B.1-like, and 3 were undetermined9,11.
All records of PCR testing in Qatar were examined in this study. Every individual that met the inclusion criteria in the national database, that is being vaccinated with BNT162b2 or mRNA-1273 and completing ≥14 days after the second vaccine dose, for each of these cohort studies, was classified based on infection status (with or without PCR-positive swab before the start of the study). Individuals were matched based on infection status on a 1:1 ratio by sex, 5-year age group, nationality (>75 nationality groups), and calendar week of first vaccine dose to control for differences in exposure risk12,13 and variant exposure6,7,9-11. Only matched samples were included in the analysis.
Further background on Qatar’s epidemic, such as on reinfections25,26, national seroprevalence surveys12,18-20, PCR surveys12, and other epidemiological studies can be found in previous publications on this epidemic6,7,12,13,27-34.
Laboratory methods
Nasopharyngeal and/or oropharyngeal swabs (Huachenyang Technology, China) were collected for PCR testing and placed in Universal Transport Medium (UTM). Aliquots of UTM were: extracted on a QIAsymphony platform (QIAGEN, USA) and tested with real-time reverse-transcription PCR (RT-qPCR) using TaqPath™ COVID-19 Combo Kits (100% sensitivity and specificity35; Thermo Fisher Scientific, USA) on an ABI 7500 FAST (ThermoFisher, USA); extracted using a custom protocol36 on a Hamilton Microlab STAR (Hamilton, USA) and tested using AccuPower SARS-CoV-2 Real-Time RT-PCR Kits (100% sensitivity and specificity37; Bioneer, Korea) on an ABI 7500 FAST; or loaded directly into a Roche cobas® 6800 system and assayed with a cobas® SARS-CoV-2 Test (95% sensitivity, 100% specificity38; Roche, Switzerland). The first assay targets the viral S, N, and ORF1ab regions. The second targets the viral RdRp and E-gene regions, and the third targets the ORF1ab and E-gene regions. Antibodies against SARS-CoV-2 in serological samples were detected using a Roche Elecsys® Anti-SARS-CoV-2 assay (99.5% sensitivity39, 99.8% specificity39,40; Roche, Switzerland), an electrochemiluminescence immunoassay that uses a recombinant protein representing the nucleocapsid (N) antigen for antibody binding. Results were interpreted according to the manufacturer’s instructions (reactive: optical density (proxy for antibody titer41) cutoff index ≥1.0 vs. non-reactive: optical density cutoff index <1.0).
All PCR tests were conducted at the Hamad Medical Corporation Central Laboratory or Sidra Medicine Laboratory, following standardized protocols.
Statistical analysis
Descriptive statistics (frequency distributions and measures of central tendency) were used to characterize study samples. Significant associations were determined using two-sided p-values. The Kaplan–Meier estimator method14 was used to estimate the cumulative risk of documented infection. Cumulative risk was defined as the proportion of individuals identified with an infection during the study period among all eligible individuals in each cohort.
Incidence rates of documented infection in each cohort were calculated by dividing the number of infection cases identified during the study by the number of person-weeks contributed by all eligible individuals in the cohort. Incidence rates and corresponding 95% CIs were estimated using a Poisson log-likelihood regression model with the STATA 17.042 stptime command. Follow-up person-time was calculated from the day each person completed 14 days after the second vaccine dose up to the infection swab, all-cause death, or end-of-study censoring (June 6, 2021). The incidence rate ratio and corresponding 95% CI were calculated using the exact method.
Statistical analyses were conducted in STATA/SE version 17.042.
Ethical approvals
The study was approved by the Hamad Medical Corporation and Weill Cornell Medicine-Qatar Institutional Review Boards with waiver of informed consent.
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.
Supplementary Material
Kaplan-Meier curves showing the cumulative incidence of documented SARS-CoV-2 infection in the national cohort of individuals who completed ≥14 days after the second vaccine dose and who had a prior PCR-confirmed infection, compared to the cumulative incidence of documented SARS-CoV-2 infection in the matched national cohort of individuals who completed ≥14 days after the second vaccine dose, but without prior PCR-confirmed infection. The curves compare vaccination with A) the BNT162b2 (Pfizer-BioNTech) vaccine and B) the mRNA-1273 vaccine. Cohorts were matched in a 1:1 ratio by sex, 5-year age group, nationality, and calendar week of the first vaccine dose. The cumulative infection incidence among the BNT162b2-vaccinated persons, with and without prior infection, was estimated at 0.16% (95% CI: 0.11-0.23%) and 1.45% (95% CI: 1.20-1.76%), respectively, after 132 days of follow-up.
Acknowledgements
We acknowledge the many dedicated individuals at Hamad Medical Corporation, the Ministry of Public Health, the Primary Health Care Corporation, and the Qatar Biobank for their diligent efforts and contributions to make this study possible. The authors are grateful for support from the Biomedical Research Program, the Biostatistics, Epidemiology, and Biomathematics Research Core, and the Genomics Core, all at Weill Cornell Medicine-Qatar, as well as for support provided by the Ministry of Public Health and Hamad Medical Corporation. The authors are also grateful for the Qatar Genome Programme for supporting 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.