Timing of Breakthrough Infection Risk After Vaccination Against SARS-CoV-2
==========================================================================

* David N. Fisman
* Nelson Lee
* Ashleigh R. Tuite

## Abstract

**Background** Provision of safe and effective vaccines has been a remarkable public health achievement during the SARS-CoV-2 pandemic. The effectiveness and durability of protection of the first two doses of SARS-CoV-2 vaccines is an important area for study, as are questions related to optimal dose combinations and dosing intervals.

**Methods** We performed a case-cohort study to generate real-world evidence on efficacy of first and second dose of SARS-CoV-2 vaccines, using a population-based case line list and vaccination database for the province of Ontario, Canada between December 2020 and October 2021. Risk of infection after vaccination was evaluated in all laboratory-confirmed vaccinated SARS-CoV-2 cases, and a 2% sample of vaccinated controls, evaluated using survival analytic methods, including construction of Cox proportional hazards models. Vaccination status was treated as a time-varying covariate.

**Results** First and second doses of SARS-CoV-2 vaccine markedly reduced risk of infection (first dose efficacy 68%, 95% CI 67% to 69%; second dose efficacy 88%, 95% CI 87 to 88%). In multivariable models, extended dosing intervals were associated with lowest risk of breakthrough infection (HR for redosing 0.64 (95% CI 0.61 to 0.67) at 6-8 weeks). Heterologous vaccine schedules that mixed viral vector vaccine first doses with mRNA second doses were significantly more effective than mRNA only vaccines. Risk of infection largely vanished during the time period 4-6 months after the second vaccine dose, but rose markedly thereafter.

**Interpretation** A case-cohort design provided an efficient means to identify strong protective effects associated with SARS-CoV-2 vaccination, particularly after the second dose of vaccine. However, this effect appeared to wane once more than 6 months had elapsed since vaccination. Heterologous vaccination and extended dosing intervals improved the durability of immune response.

## Introduction

The rapid introduction of highly effective vaccines against SARS-CoV-2 represents a remarkable scientific achievement amidst the current global pandemic (1). It has been estimated that these vaccines have averted hundreds of thousands of deaths (2, 3). While mRNA vaccines represent the first application of a novel and highly immunogenic vaccine construct, viral vector vaccines have also been extremely effective in preventing severe illness and death (4). The introduction of novel vaccines against a virus that has been recognized for only two years (5), means that there is of necessity a great deal of uncertainty related to their application. The relative protection afforded by vaccines against emerging variants of concern (6), vaccine safety (7, 8), the immunogenicity of mixed vaccine schedules (9), optimal spacing between administration of doses (10), the number of doses required for a full primary series (11), and the necessity of future boosting (12), are all the subject of active study.

The Canadian province of Ontario represents a large (population 14.6 million) and diverse jurisdiction, with high levels of 2-dose SARS-CoV-2 vaccine coverage (approximately 77% as of December 2021) (13, 14). The province has made use of several different vaccine formulations, predominantly BNT162b2 (BNT/Pfizer, “Cominrty”), mRNA-1273 (Moderna, “SpikeVax”), and ChAdOX1-S (AstraZeneca, “Covishield”). Ontario’s vaccination program has prioritized risk groups over time. For example, the program initially focussed on vaccinating residents of long-term care facilities and health-care workers (15). Once higher risk groups were vaccinated, eligibility was expanded using a descending age-based approach. The province also consciously chose to offer extended dosing intervals to lower risk populations, due to supply constraints and the potential advantages of rapid, maximal coverage with first doses of vaccine (16). Pediatric vaccination programs for those aged 5 to 11 years were introduced in in late November 2021 (17).

Beginning in autumn 2021, third doses have been offered to a growing list of individuals identified as being at elevated risk due to age, medical comorbidity or congregate living residence. Although there was some evidence suggesting that the use of prolonged dosing intervals in the primary series had resulted in durable immunity that made third doses unnecessary (18), the province expanded access to third vaccine to all aged 50 and over December 2021 (19). In late December, the rapid spread of the Omicron variant of concern in the province resulted in an expansion of third dose eligibility to all aged 18 and over who received their second dose at least 3 months previously.

A large, highly-vaccinated population like Ontario’s, which has received a several different vaccination regimens, with variable dosing intervals, provides an opportunity to generate real-world evidence on efficacy of vaccination against breakthrough infections. The province’s vaccination database and COVID-19 case line provide detailed information on vaccine administration as well as individual-level covariates that can be used to adjust for confounding by factors associated with both infection risk and vaccination priority. Our objective was to use these data to generate estimates of 1- and 2-dose vaccination protection over time, against breakthrough SARS-CoV-2 infection; we also performed exploratory analyses to generate hypotheses related to optimal dosing intervals and the relative efficacy of different vaccine combinations.

## Methods

### Data Sources

We created a cohort of individuals who received vaccination against SARS-CoV-2 in Ontario, without record of infection prior to the date of receipt of first vaccine dose. We included all reported cases of infection meeting these criteria, as well as a 2% random sample of vaccinated individuals with no record of infection. Cases were identified in the Province’s Case and Contact Management (CCM) database as described elsewhere (20, 21); we included only cases with a unique “pseudo-health card number”, which permitted linkage with the provincial vaccination database. Vaccination information on cases and controls was extracted from the Province’s COVaxON dataset (21), which includes dosage dates, vaccines used, basic demographic information, and fields on indication for vaccination, including residence in a congregate setting (including residence in a nursing home, retirement community, or other congregate setting), priority health condition, and healthcare worker status. As with cases, the pool of potential controls was limited to individuals with a pseudo-health card number. A flow diagram outlining creation of the cohort is presented in **Figure 1**.

![Figure 1.](http://medrxiv.org/https://www.medrxiv.org/content/medrxiv/early/2022/01/05/2022.01.04.22268773/F1.medium.gif)

[Figure 1.](http://medrxiv.org/content/early/2022/01/05/2022.01.04.22268773/F1)

Figure 1. Flow Diagram for Creation of Case-Cohort Sample

We considered individuals to have been vaccinated with a first dose of vaccine during time at risk 14 or more days after the date of the first vaccine dose; individuals were considered vaccinated with two doses of vaccine during time at risk 14 or more days after second vaccine dose. While some researchers, including ourselves (21-23), have previously defined two dose-vaccinated time as commencing 7 or more days after second dose, we use the longer window for consistency with Canadian public health authorities (24).

### Analysis

We performed a case-cohort analysis using survival analytic methods. As person time data were available, we were able to estimate hazard ratios directly using Cox proportional hazards models, rather than approximate such ratios via risk set sampling. Controls were reweighted by a factor of 50 to account for 2% selection probability. We stratified within-person time at risk according to an individual’s vaccination status, which was treated as a time-varying covariate.

We created Cox proportional hazards models including only vaccination as an exposure (unvaccinated, single-dose vaccinated, or fully vaccinated), as well as adjusted models that adjusted for variables that we expected to confound observed vaccine effect. We evaluated the proportional hazards assumption both graphically using log-log plots, and statistically based on Schoenfeld residuals (25). As we found hazards of infection following vaccination to be non-proportional through inspection of log-log plots and testing of Schoenberg residuals (**Supplementary Figure**), we re-ran our Cox proportional hazards models with time at risk further stratified by 2-month blocks, within a given vaccination status.

Lastly, we performed restriction analyses in which we restricted cases individuals with infection due to a non-variant of concern viral strain, to infection with N501Y+ variants of concern (alpha, beta and gamma variants), and to infection with likely delta infection. Classification of likely variant status of the infecting strain was as described previously (21). We conducted all analyses in Stata version 15.1 using sts, stox and related commands. The study was conducted in accordance with the STROBE guidelines for observational research (26), and received ethics approval from the Research Ethics Board at the University of Toronto.

## Results

Our final cohort consisted of 52,948 cases and 211,373 controls. Cases and controls differed significantly according to age distribution, sex, vaccine combination received, documented priority health condition, residence in a congregate setting, and healthcare worker status, as well as interval between first and second doses (**Table 1**). Crude time to breakthrough infection by vaccination status is presented as a Kaplan-Meier curve (**Figure 2**); vaccination was associated with significant reduction in infection risk, with 2-dose vaccination reducing risk relative to single dose vaccination (P < 0.001 by log-rank test for all comparisons).

View this table:
[Table 1:](http://medrxiv.org/content/early/2022/01/05/2022.01.04.22268773/T1)

Table 1: Characteristics of SARS-CoV-2-Vaccinated Cases and Controls, Ontario, Canada

![Figure 2.](http://medrxiv.org/https://www.medrxiv.org/content/medrxiv/early/2022/01/05/2022.01.04.22268773/F2.medium.gif)

[Figure 2.](http://medrxiv.org/content/early/2022/01/05/2022.01.04.22268773/F2)

Figure 2. Time to Infection Without SARS-CoV-2 Vaccination, or After 1 or 2 Vaccine Doses
Vaccination status is treated as a time-varying covariate and defined by time since vaccination and cumulative doses, as described in the text. X-axis, person-time; Y-axis, cumulative incidence of infection. Shaded areas represent 95% confidence intervals.

We constructed a Cox proportional hazards model that included only vaccination status; the hazard ratio for breakthrough infection after a single dose of vaccine was 0.32 (95% CI 0.31 to 0.33); after a second dose of vaccine the hazard ratio declined to 0.12 (95% CI 0.12 to 0.13), corresponding to effectiveness of 68% and 88% against infection. We found little change in these after adjusting for potential confounders in multivariable models (**Table 2**). In multivariable analyses we identified a strong U-shaped relationship between dosing interval and risk of breakthrough infection, with risk lowest when vaccine doses were administered 6-8 weeks apart. Among individuals who received two doses of vaccine, increased risk of breakthrough infection was seen in those who received two doses of ChAdOx1-S, relative to two doses of mRNA vaccine. However, for individuals who received an initial dose of ChAdOx1-S, there was a significant reduction in breakthrough risk following a second dose of mRNA vaccine, in contrast to receipt of two doses of mRNA vaccine. Risk of breakthrough was highest in younger adults (aged 20-49), and was elevated in healthcare workers, those residing in congregate living settings, and in those identified as having priority medical conditions.

View this table:
[Table 2.](http://medrxiv.org/content/early/2022/01/05/2022.01.04.22268773/T2)

Table 2. Results of Multivariable Cox Proportional Hazards Model

Models that stratified time at risk both by vaccination status and 60-day time blocks since vaccination generated similar results to initial multivariable models, demonstrating high efficacy of first and second vaccine doses (**Supplementary Table**). The lowest hazard ratio for breakthrough infection was seen between four and 6 months after second vaccine dose but increased after 6 months of observation (**Figure 3**). In exploratory restriction analyses in which we included only cases with non-variant of concern infection, cases with N501Y+ infection (alpha, beta or gamma variants), or cases with delta variant infection, we found significant protection against breakthrough infection with both 1 and 2 doses of vaccine, against all variants. However, there was significant heterogeneity across variants (P for heterogeneity < 0.001 for both 1 and 2 doses), with highest efficacy after two doses seen against N501Y+ infection (96.9%, 95% CI 96.7% to 97.1%), and lower efficacy seen against delta variant (66.1%, 95% CI 62.5% to 69.4%) (**Table 3**).

View this table:
[Table 3.](http://medrxiv.org/content/early/2022/01/05/2022.01.04.22268773/T3)

Table 3. Strain- and Dose-Specific Efficacy of SARS-CoV-2 Vaccines in Restriction Analyses

![Figure 3.](http://medrxiv.org/https://www.medrxiv.org/content/medrxiv/early/2022/01/05/2022.01.04.22268773/F3.medium.gif)

[Figure 3.](http://medrxiv.org/content/early/2022/01/05/2022.01.04.22268773/F3)

Figure 3. Relative Hazard of Infection by Time Since First or Second SARS-CoV-2 Vaccine Dose
Hazard ratios for infection by vaccination status, with time stratified by 60-day blocks since vaccination. Circles, hazard ratios; vertical lines 95% CI.

## Discussion

We find that in a large cohort of vaccinated individuals in Ontario, Canada vaccination provided remarkable protection against breakthrough infection, with increased protection afforded by a second vaccine dose relative to a single dose. We also find that maximal protection against infection is provided 4-6 months after the second dose of vaccine, at which time protection against infection was nearly complete. However, we found that infection risk increased again after 6 months, suggesting the need for third vaccine doses at this interval (11). Consistent with effects demonstrated using viral neutralization assays (10), we find that the extended interval between first and second vaccine doses has resulted in enhanced immune protection, with a relative reduction in risk of approximately 40% when second vaccine doses were provided 6-8 weeks after first doses. We also identified significant protection against infection associated with heterologous vaccine schedules, in particular the use of initial dose ChAdOx1-S, followed by a second dose of mRNA vaccine, though viral vector vaccines are no longer in wide use in Canada due to concerns about risk of vaccine induced thrombotic thrombocytopenia (VITT)(7). However, the observed effects have substantial face validity and are consistent with the previously described immune enhancement seen with heterologous vaccination schedules as well as the high efficacy of mRNA booster vaccination following non-mRNA SARS-CoV-2 vaccination (27-31).

Our findings on vaccine efficacy are consistent with those reported by others using more typical cohort designs for study of vaccine efficacy. Our application of a case-cohort design, with both cases and controls contributing person-time to unvaccinated, and single- and double-dose vaccinated exposures, allowed us to generate these estimates without access to an unvaccinated control group. The fact that all individuals in our study chose to be vaccinated potentially eliminates confounding by factors associated with vaccination choices.

Like any observational study, ours is potentially subject to confounding by measured and unmeasured factors, though inclusion of available covariates in multivariable models did not significantly change our estimates of vaccine effectiveness. We were also surprised at the sparse nature of data on healthcare worker status, congregate living status and priority health conditions in the CoVaxON dataset. The small number of individuals with these classifications suggests that the dataset is insensitive for these elements. However, inasmuch as CoVaxON data would have been recorded without knowledge of a case’s subsequent infection status these missing data should not result in bias.

In summary, we use a case-cohort study with stratified person-time analysis of vaccination status to demonstrate the strong protection afforded by SARS-CoV-2 vaccines in Ontario, Canada. Despite this strong protection, we do find that an increase in hazard greater than 6 months after second dose of vaccine; our analysis also suggests that increased intervals between vaccine doses, as widely implemented in Canada, provides more durable protection against infection than shorter dosing intervals.

## Data Availability

Data are property of the Government of Ontario. Questions about data and analysis should be directed to the study authors.

View this table:
[Supplementary Table.](http://medrxiv.org/content/early/2022/01/05/2022.01.04.22268773/T4)

Supplementary Table. Multivariable Cox Proportional Hazards Model with Stratified Time Since Vaccination

## Acknowledgements

The authors with to thank the staff at Public Health Ontario and Ontario’s public health units for collecting, sequencing, analyzing, and providing access to the data used for this analysis.

## Footnotes

*   Funding statement: The research was supported by a grant to DNF from the Canadians Institutes for Health Research (2019 COVID-19 rapid researching funding OV4-170360).

*   Declaration of competing interests: DNF has served on advisory boards related to influenza and SARS-CoV-2 vaccines for Seqirus, Pfizer, Astrazeneca and Sanofi-Pasteur Vaccines, and has served as a legal expert on issues related to COVID-19 epidemiology for the Elementary Teachers Federation of Ontario and the Registered Nurses Association of Ontario. ART was employed by the Public Health Agency of Canada when the research was conducted. The work does not represent the views of the Public Health Agency of Canada.

*   Received January 4, 2022.
*   Revision received January 4, 2022.
*   Accepted January 5, 2022.


*   © 2022, Posted by Cold Spring Harbor Laboratory

This pre-print is available under a Creative Commons License (Attribution-NonCommercial-NoDerivs 4.0 International), CC BY-NC-ND 4.0, as described at [http://creativecommons.org/licenses/by-nc-nd/4.0/](http://creativecommons.org/licenses/by-nc-nd/4.0/)

## References

1.  1.Kreier F. ‘Unprecedented achievement’: who received the first billion COVID vaccinations? Nature. 2021.
    
    

2.  2.Mesle MM, Brown J, Mook P, Hagan J, Pastore R, Bundle N, et al. Estimated number of deaths directly averted in people 60 years and older as a result of COVID-19 vaccination in the WHO European Region, December 2020 to November 2021. Euro Surveill. 2021;26(47).
    
    

3.  3.Gupta S, Cantor J, Simon KI, Bento AI, Wing C, Whaley CM. Vaccinations Against COVID-19 May Have Averted Up To 140,000 Deaths In The United States. Health Aff (Millwood). 2021;40(9):1465–72.
    
    

4.  4.Higdon MM, Wahl B, Jones CB, Rosen JG, Truelove SA, Baidya A, et al. A systematic review of COVID-19 vaccine efficacy and effectiveness against SARS-CoV-2 infection and disease. medRxiv. 2021:2021.09.17.21263549.
    
    

5.  5.Wu P, Hao X, Lau EHY, Wong JY, Leung KSM, Wu JT, et al. Real-time tentative assessment of the epidemiological characteristics of novel coronavirus infections in Wuhan, China, as at 22 January 2020. Euro Surveill. 2020;25(3).
    
    

6.  6.Jablonska K, Aballea S, Toumi M. The real-life impact of vaccination on COVID-19 mortality in Europe and Israel. Public Health. 2021;198:230–7.
    
    

7.  7.Greinacher A, Selleng K, Palankar R, Wesche J, Handtke S, Wolff M, et al. Insights in ChAdOx1 nCoV-19 vaccine-induced immune thrombotic thrombocytopenia. Blood. 2021;138(22):2256–68.
    
    

8.  8.Witberg G, Barda N, Hoss S, Richter I, Wiessman M, Aviv Y, et al. Myocarditis after Covid-19 Vaccination in a Large Health Care Organization. N Engl J Med. 2021;385(23):2132–9.
    
    

9.  9.Schmidt T, Klemis V, Schub D, Mihm J, Hielscher F, Marx S, et al. Immunogenicity and reactogenicity of heterologous ChAdOx1 nCoV-19/mRNA vaccination. Nat Med. 2021;27(9):1530–5.
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=http://www.n&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2022%2F01%2F05%2F2022.01.04.22268773.atom) 

10. 10.Payne RP, Longet S, Austin JA, Skelly DT, Dejnirattisai W, Adele S, et al. Immunogenicity of standard and extended dosing intervals of BNT162b2 mRNA vaccine. Cell. 2021;184(23):5699–714 e11.
    
    

11. 11.Barda N, Dagan N, Cohen C, Hernan MA, Lipsitch M, Kohane IS, et al. Effectiveness of a third dose of the BNT162b2 mRNA COVID-19 vaccine for preventing severe outcomes in Israel: an observational study. Lancet. 2021;398(10316):2093–100.
    
    

12. 12.Burki T. Booster shots for COVID-19-the debate continues. Lancet Infect Dis. 2021;21(10):1359–60.
    
    

13. 13.Statistics Canada. Population estimates, quarterly. Table: 17-10-0009-01 (formerly CANSIM 051-0005). Available via the Internet at [https://www150.statcan.gc.ca/t1/tbl1/en/tv.action?pid=1710000901](https://www150.statcan.gc.ca/t1/tbl1/en/tv.action?pid=1710000901). Last accessed May 29, 2020. 2020.
    
    

14. 14.Ontario Agency for Health Protection and Promotion (Public Health Ontario). COVID-19 vaccine uptake in Ontario: December 14, 2020 to November 28, 2021. Available via the Internet at publichealthontario.ca/-/media/documents/ncov/epi/covid-19-vaccine-uptake-ontario-epi-summary.pdf?la=en. Last accessed December 7, 2021. 2021.
    
    

15. 15.Ontario Ministry of Health. Ontario’s COVID-19 vaccination plan. Available via the Internet at [https://covid-19.ontario.ca/ontarios-covid-19-vaccination-plan#our-three-phased-vaccination-plan](https://covid-19.ontario.ca/ontarios-covid-19-vaccination-plan#our-three-phased-vaccination-plan). Last accessed December 8, 2021. 2021.
    
    

16. 16.Tuite AR, Zhu L, Fisman DN, Salomon JA. Alternative Dose Allocation Strategies to Increase Benefits From Constrained COVID-19 Vaccine Supply. Ann Intern Med. 2021;174(4):570–2.
    
    

17. 17.Ontario Ministry of Health. COVID-19 Vaccine Bookings to Open For All Children Aged Five to 11. Families can book appointments through a variety of channels starting November 23rd. Available via the Internet at [https://news.ontario.ca/en/release/1001195/covid-19-vaccine-bookings-to-open-for-all-children-aged-five-to-11](https://news.ontario.ca/en/release/1001195/covid-19-vaccine-bookings-to-open-for-all-children-aged-five-to-11). Last accessed December 8, 2021. 2021.
    
    

18. 18.Skowronski DM, Setayeshgar S, Febriani Y, Ouakki M, Zou M, Talbot D, et al. Two-dose SARS-CoV-2 vaccine effectiveness with mixed schedules and extended dosing intervals: test-negative design studies from British Columbia and Quebec, Canada. medRxiv. 2021:2021.10.26.21265397.
    
    

19. 19.Health OMo. COVID-19 vaccine third dose recommendations. Available via the Internet at [https://www.health.gov.on.ca/en/pro/programs/publichealth/coronavirus/docs/vaccine/COVID-19\_vaccine\_third\_dose\_recommendations.pdf](https://www.health.gov.on.ca/en/pro/programs/publichealth/coronavirus/docs/vaccine/COVID-19\_vaccine_third_dose_recommendations.pdf). Last accessed December 8, 2021. 2021.
    
    

20. 20.Fisman DN, Greer AL, Hillmer M, O’Brien SF, Drews SJ, Tuite AR. COVID-19 case age distribution: correction for differential testing by age. medRxiv. 2020:2020.09.15.20193862.
    
    

21. 21.Fisman DN, Tuite AR. Evaluation of the relative virulence of novel SARS-CoV-2 variants: a retrospective cohort study in Ontario, Canada. CMAJ. 2021;193(42):E1619–E25.
    
    [Abstract/FREE Full Text](http://medrxiv.org/lookup/ijlink/YTozOntzOjQ6InBhdGgiO3M6MTQ6Ii9sb29rdXAvaWpsaW5rIjtzOjU6InF1ZXJ5IjthOjQ6e3M6ODoibGlua1R5cGUiO3M6NDoiQUJTVCI7czoxMToiam91cm5hbENvZGUiO3M6NDoiY21haiI7czo1OiJyZXNpZCI7czoxMjoiMTkzLzQyL0UxNjE5IjtzOjQ6ImF0b20iO3M6NTA6Ii9tZWRyeGl2L2Vhcmx5LzIwMjIvMDEvMDUvMjAyMi4wMS4wNC4yMjI2ODc3My5hdG9tIjt9czo4OiJmcmFnbWVudCI7czowOiIiO30=) 

22. 22.Dagan N, Barda N, Kepten E, Miron O, Perchik S, Katz MA, et al. BNT162b2 mRNA Covid-19 Vaccine in a Nationwide Mass Vaccination Setting. N Engl J Med. 2021;384(15):1412–23.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1056/NEJMOA2101765&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=http://www.n&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2022%2F01%2F05%2F2022.01.04.22268773.atom) 

23. 23.Chung H, He S, Nasreen S, Sundaram ME, Buchan SA, Wilson SE, et al. Effectiveness of BNT162b2 and mRNA-1273 COVID-19 vaccines against symptomatic SARS-CoV-2 infection and severe COVID-19 outcomes in Ontario, Canada. medRxiv. 2021:2021.05.24.21257744.
    
    

24. 24.Ontario Ministry of Health. COVID-19 Fully Vaccinated Status in Ontario. Available via the Internet at [https://www.health.gov.on.ca/en/pro/programs/publichealth/coronavirus/docs/vaccine/COVID-19\_fully\_vaccinated\_status\_ontario.pdf](https://www.health.gov.on.ca/en/pro/programs/publichealth/coronavirus/docs/vaccine/COVID-19\_fully_vaccinated_status_ontario.pdf). Last accessed December 2, 2021. 2021.
    
    

25. 25.Grambsch PM, Therneau TM. Proportional hazards tests and diagnostics based on weighted residuals. Biometrika. 1994; 81:515–26.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1093/biomet/81.3.515&link_type=DOI) 
    
    [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=A1994PP36700006&link_type=ISI) 

26. 26.von Elm E, Altman DG, Egger M, Pocock SJ, Gotzsche PC, Vandenbroucke JP, et al. The Strengthening the Reporting of Observational Studies in Epidemiology (STROBE) statement: guidelines for reporting observational studies. Bull World Health Organ. 2007;85(11):867–72.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.2471/BLT.07.045120&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=18038077&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2022%2F01%2F05%2F2022.01.04.22268773.atom) 
    
    [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000250858300014&link_type=ISI) 

27. 27.Ostadgavahi AT, Booth R, Sisson G, McMullen N, Warhuus M, Robertson P, et al. Heterologous immunization with Covishield and Pfizer vaccines against SARS-CoV-2 elicits a robust humoral immune response. J Infect Dev Ctries. 2021;15(5):653–6.
    
    

28. 28.Munro APS, Janani L, Cornelius V, Aley PK, Babbage G, Baxter D, et al. Safety and immunogenicity of seven COVID-19 vaccines as a third dose (booster) following two doses of ChAdOx1 nCov-19 or BNT162b2 in the UK (COV-BOOST): a blinded, multicentre, randomised, controlled, phase 2 trial. Lancet. 2021.
    
    

29. 29.European Medicines Agency. Heterologous primary and booster COVID-19 vaccination. Evidence based regulatory considerations. Available via the Internet at [https://www.ema.europa.eu/en/documents/report/heterologous-primary-booster-covid-19-vaccination-evidence-based-regulatory-considerations\_en.pdf](https://www.ema.europa.eu/en/documents/report/heterologous-primary-booster-covid-19-vaccination-evidence-based-regulatory-considerations_en.pdf). Last accessed January 2, 2022.; 2021.
    
    

30. 30.Nordstrom P, Ballin M, Nordstrom A. Effectiveness of heterologous ChAdOx1 nCoV-19 and mRNA prime-boost vaccination against symptomatic Covid-19 infection in Sweden: A nationwide cohort study. Lancet Reg Health Eur. 2021;11:100249.
    
    

31. 31.World Health Organization. Interim recommendations for heterologous COVID-19 vaccine schedules. Available via the Internet at [https://www.who.int/publications/i/item/WHO-2019-nCoV-vaccines-SAGE-recommendation-heterologous-schedules](https://www.who.int/publications/i/item/WHO-2019-nCoV-vaccines-SAGE-recommendation-heterologous-schedules). Last accessed January 2, 2022. 2021.