Skip to main content
Advertisement
Browse Subject Areas
?

Click through the PLOS taxonomy to find articles in your field.

For more information about PLOS Subject Areas, click here.

  • Loading metrics

Vaccine effectiveness against hospitalization among adolescent and pediatric SARS-CoV-2 cases between May 2021 and January 2022 in Ontario, Canada: A retrospective cohort study

Abstract

Background

Vaccines against SARS-CoV-2 have been shown to reduce risk of infection as well as severe disease among those with breakthrough infection in adults. The latter effect is particularly important as immune evasion by Omicron variants appears to have made vaccines less effective at preventing infection. Therefore, we aimed to quantify the protection conferred by mRNA vaccination against hospitalization due to SARS-CoV-2 in adolescent and pediatric populations.

Methods

We retrospectively created a cohort of reported SARS-CoV-2 case records from Ontario’s Public Health Case and Contact Management Solution among those aged 4 to 17 linked to vaccination records from the COVaxON database on January 19, 2022. We used multivariable logistic regression to estimate the association between vaccination and hospitalization among SARS-CoV-2 cases prior to and during the emergence of Omicron.

Results

We included 62 hospitalized and 27,674 non-hospitalized SARS-CoV-2 cases, with disease onset from May 28, 2021 to December 4, 2021 (Pre-Omicron) and from December 23, 2021 to January 9, 2022 (Omicron). Among adolescents, two mRNA vaccine doses were associated with an 85% (aOR = 0.15; 95% CI: [0.04, 0.53]; p<0.01) lower likelihood of hospitalization among SARS-CoV-2 cases caused by Omicron. Among children, one mRNA vaccine dose was associated with a 79% (aOR = 0.21; 95% CI: [0.03, 0.77]; p<0.05) lower likelihood of hospitalization among SARS-CoV-2 cases caused by Omicron. The calculation of E-values, which quantifies how strong an unmeasured confounder would need to be to nullify our findings, suggest that these effects are unlikely to be explained by unmeasured confounding.

Conclusions

Despite immune evasion by SARS-CoV-2 variants, vaccination continues to be associated with a lower likelihood of hospitalization among adolescent and pediatric Omicron (B.1.1.529) SARS-CoV-2 cases, even when the vaccines do not prevent infection. Continued efforts are needed to increase vaccine uptake among adolescent and pediatric populations.

Introduction

Severe Acute Respiratory Syndrome (SARS-CoV-2) has caused more than 6.5 million deaths globally [1]. Safe and effective vaccines to prevent SARS-CoV-2 infection and severe outcomes have been approved since late 2020 for adults [2]. With the emergence of viral variants B.1.617.2 (Delta) in May 2021 and B.1.1.529 (Omicron) in late November 2021, decreased vaccine effectiveness against infection was observed [36]. Given that the goal of vaccination is to prevent death, severe disease, and overall disease burden, it is important to consider how well vaccines achieve these goals among individuals infected with SARS-CoV-2 [7].

Four recent studies focused on quantifying the real-world effectiveness of SARS-CoV-2 vaccination against hospitalization among adolescent or pediatric populations [811]. All the studies used uninfected controls, and resulting vaccine effectiveness estimates reflect the joint risk of SARS-CoV-2 infection and the risk of hospitalization conditional on infection [12]. Olson et al. [8] showed that BNT162b2 was 94% effective against hospitalization among adolescents infected with the Delta variant (B.1.617.2). After the emergence of the Omicron variant (B.1.1.529), two dose vaccine effectiveness against hospitalization was estimated to be 73% among adolescents ages 12 to 17 and 48% among pediatric patients ages 5 to 11 [9]. In a recent Morbidity and Mortality Weekly Report (MMWR), two dose vaccine effectiveness against hospitalization was between 73% and 94% among adolescent and pediatric populations [10]. In a study among hospitalized patients, two dose vaccine effectiveness against SARS-CoV-2 hospitalization was 93% during the Delta-dominant period and 40% during the Omicron-dominant period among adolescents, and 68% during the Omicron-dominant period among pediatric patients [11]. In contrast to the approach used in these studies, we conditioned upon SARS-CoV-2 infection in our analysis to estimate the direct association between vaccination and hospitalization risk among pediatric and adolescent SARS-CoV-2 cases [13].

The Canadian province of Ontario is a large (population 14.6 million) and diverse jurisdiction with high levels of SARS-CoV-2 vaccine coverage [14]. Approximately 87% of Ontario residents ages 12 to 17 and 53% of Ontario residents ages 5 to 11 received at least one SARS-CoV-2 dose as of January 2022 [15]. Most individuals ages 12–17 were eligible to receive the BNT162b2 (Pfizer-BioNTech, Comirnaty) vaccine beginning on May 28, 2021 [16], followed by ages 5 to 11 on November 28, 2021 [17]. Health Canada authorized Moderna Spikevax on August 27, 2021 for ages 12 to 17, but it was not authorized for ages 6 to 11 until after our study period [18]. We calculated one dose and two dose mRNA vaccine effectiveness against hospitalization among adolescent and pediatric SARS-CoV-2 cases to examine if vaccines prevent hospitalization when they fail to prevent SARS-CoV-2 infection.

Methods

Data

Robust public health surveillance systems in Ontario enable individual-level linkage of the SARS-CoV-2 vaccination database and the reportable disease database. Confirmed SARS-CoV-2 cases were identified in Ontario’s Public Health Case and Contact Management Solution (CCM). The CCM includes patient demographics (i.e., sex, age, comorbidities), geographic location (i.e., health region [5 units at the sub-provincial level]), and case characteristics (i.e., test date, symptom onset date, hospital admission and discharge dates) for all laboratory-confirmed SARS-CoV-2 cases in Ontario. SARS-CoV-2 vaccination information was identified from the provincial COVaxON database. COVaxON data includes vaccine administration information (i.e., dose dates, dose locations, dose indication, vaccine product) for Ontario residents with a provincial health insurance number. We linked data from CCM and COVaxON using a unique pseudo identifier present in both datasets. Median household income and percent visible minority at the census subdivision level were from the 2016 Census of Population [19].

Study design

We conducted a retrospective cohort study of individuals aged 4–17 years who had a confirmed SARS-CoV-2 infection in Ontario, Canada. These individuals had a positive reverse transcription real-time polymerase chain reaction (PCR) test between May 28, 2021 and December 4, 2021 (pre-Omicron), or between December 23, 2021 and January 9, 2022 (Omicron). Our study period ended because Ontario restricted publicly funded PCR testing to only select members of the population [20].

Our analyses only included incident SARS-CoV-2 infections (re-infections were excluded). We excluded SARS-CoV-2 cases where individuals had received a vaccine not approved for use in this population in Canada (e.g., Johnson & Johnson/Janssen COVID-19 vaccine, Oxford-AstraZeneca COVID-19 vaccine), or had received three vaccine doses (ages 12–17) or two vaccine doses (ages 4–11). During the study period, adolescents were not eligible for third doses (i.e., booster dose), and the recommended interval between first and second doses for children was eight weeks [21]. We extracted the data on January 19, 2022, but only included test dates up to January 9, 2022, to account for delays between testing, hospitalization, and reporting. SARS-CoV-2 cases aged 4 years were included in the analysis because cases among pediatric patients aged 4 and aged 5 were in the same 2-year age group in these data, and pediatric patients who turned 5 in the 2021 calendar year were eligible for vaccination.

Measures

The outcome was hospitalization due to SARS-CoV-2. Hospitalizations were identified by a reported hospital admission date, or a reported hospitalization or ICU admission due to SARS-CoV-2 in the CCM. The exposure was SARS-CoV-2 vaccination. We considered individuals one-dose vaccinated 14 or more days after the date the first vaccine dose was administered; individuals were considered two-dose vaccinated 14 or more days after the date the second vaccine dose was administered. Individuals infected with SARS-CoV-2 within 13 days of being vaccinated were excluded from our analyses.

Case onset date, age, sex, asthma, immunocompromising condition, and health region were included as confounders [11, 22]. A Directed Acyclic Graph (DAG) outlining the main hypothesized relationships between the variables is presented in S1 Fig [23]. Case onset date was the date of symptom onset for symptomatic cases and the specimen collection date for asymptomatic cases. In the CCM, age was reported in 2-year age groups, and asthma and immunocompromising condition were reported by the patient or provider. Health region included northern, eastern, central, or western Ontario, or Toronto [24]. Median household income and percent visible minority in 2016 at the census subdivision level were included as confounders in a sensitivity analysis.

Statistical analysis

Standardized differences (SD) [25, 26] and statistical tests (Fisher’s exact tests and unpaired t-tests) were used to compare differences in baseline characteristics between hospitalized and non-hospitalized cases. Standardized mean differences are not impacted by sample size, in contrast to the common statistical tests used in descriptive tables. A standardized difference greater than 0.10 represents a lack of covariate balance [27].

We used multivariable logistic regression to calculate the association between vaccination and hospitalization among adolescent (ages 12–17) and pediatric (ages 4–11) SARS-CoV-2 cases, while adjusting for case onset date, age, sex, asthma, immunocompromising condition, and health region. We examined one-dose and two-dose vaccine effectiveness among the adolescent population prior to the emergence (May 28, 2021 to December 4, 2021) and after the emergence of Omicron (December 23, 2021 to January 9, 2022), and one-dose vaccine effectiveness among the pediatric population after the emergence of Omicron. Ontario residents under age 12 were ineligible for SARS-CoV-2 vaccination during most of the pre-Omicron time period. Dates were chosen to align with the rise in the prevalence of the Omicron variant in Ontario. In the pre-Omicron period, less than 10% of SARS-CoV-2 samples had S-gene target failure and in the latter period, more than 90% of SARS-CoV-2 samples had S-gene target failure [28]. Vaccine effectiveness (VE) was calculated using the formula VE = (1-aOR)*100%.

We performed additional analyses to quantify the impact of unmeasured confounding and to explore the robustness of our results to our assumptions. First, we calculated E-values based on the results of our multivariable logistic regression models [29]. An E-value quantifies how strong an unmeasured confounder would need to be to explain away the association between SARS-CoV-2 vaccination and hospitalization, conditional on other measured covariates [30]. We calculated an E-values using the following formula [29]: E-value = . Although there is no accepted threshold in the literature for an E-value that demonstrates robustness to unmeasured confounding [31], an E-value is a useful measure to quantify the strength of unmeasured confounding needed to alter our findings [32]. We conducted four model-based sensitivity analyses. First, we allowed for a 14-day delay between case onset date and hospitalization, where SARS-CoV-2 cases with onset dates after January 5, 2022 were excluded. In our main analysis, we allowed for a 10-day delay between case onset date and hospitalization. Second, we repeated our analyses with adjustment for specimen collection date as opposed to case onset date to lessen the impact of differential case onset dates by symptom status. Third, we only included symptomatic individuals as opposed to all SARS-CoV-2 cases. Fourth, we examined the impact of including median household income and percent visible minority at the census subdivision level as confounders in addition to sex, asthma, immunocompromising condition, age, and case onset date in a multivariable logistic regression model. Individual-level measures of socioeconomic status and race were not available in our data and are therefore unmeasured confounders in our study [11]. This study is reported as per the Strengthening the Reporting of Observational Studies in Epidemiology (STROBE) guideline.

Ethics

We received ethics approval for this study from the Research Ethics Board at the University of Toronto (#00031358). This study includes secondary analysis of public health surveillance data collected under Ontario’s Health Protection and Promotion Act. Therefore, the need for informed consent was waived by the ethics committee.

Results

In total, 62 hospitalized SARS-CoV-2 cases and 27,674 non-hospitalized SARS-CoV-2 cases were included (Fig 1; Table 1). In the full cohort across both age groups and the two time periods, individuals hospitalized with SARS-CoV-2 were more likely to be unvaccinated, immunocompromised, have an asthma diagnosis, live in the West or North health regions, and have an earlier case onset date (Table 1). We were unable to present frequencies, percentages, means, and standard deviations for the cohort stratified by age group and time period due to small cell sizes, however the relationships were quantified using standardized differences and statistical tests (S1 Table).

thumbnail
Table 1. Description of adolescent and pediatric SARS-CoV-2 cases by hospitalization status (n = 27,736).

https://doi.org/10.1371/journal.pone.0283715.t001

In the pre-Omicron period after adjustment for sex, asthma, immunocompromising condition, age, health region, and case onset date, we observed that one and two SARS-CoV-2 vaccine doses were not significantly associated (adjusted OR (aOR) one dose = 3.09 [95% CI: 0.31, 16.84], p = 0.26; aOR two doses = 0.90 [95% CI: 0.22, 2.86), p = 0.15) with additional protection against hospitalization among adolescent SARS-CoV-2 cases (Table 2). In the Omicron period after adjustment for sex, asthma, immunocompromising condition, age, health region, and case onset date, one SARS-CoV-2 dose was not significantly associated (aOR = 0.90 [95% CI: 0.05, 5.83], p = 0.93) with additional protection against hospitalization, but two SARS-CoV-2 doses were associated with an 85% (aOR = 0.15, [95% CI: 0.04, 0.53], p<0.01) lower likelihood of hospitalization among adolescent SARS-CoV-2 cases. After adjustment for sex, asthma, immunocompromising condition, age, health region, and case onset date, one SARS-CoV-2 vaccine dose was associated with a 79% (aOR = 0.21, [95% CI: 0.03, 0.77], p<0.05) lower likelihood of hospitalization among pediatric SARS-CoV-2 cases (ages 4–11) in the Omicron period (Table 2).

thumbnail
Table 2. Adjusted odds ratios of the association between vaccination status and hospitalization among adolescent and pediatric SARS-CoV-2 cases (n = 27,736).

https://doi.org/10.1371/journal.pone.0283715.t002

The E-value for the adjusted association between vaccination with two doses and the risk of hospitalization in our primary analysis among the adolescent population in the Omicron period (i.e., aOR = 0.15) was 12.8 [29]. An unmeasured confounder would need to be independently associated with vaccination and hospitalization by a 12.8-fold risk ratio to result in a null association, after adjusting for sex, asthma, immunocompromising condition, age, health region, and case onset date. Similarly, the E-value for the upper 95% CI bound (i.e., 95% CI upper bound = 0.53) is 3.2. A set of unmeasured confounders would need to be independently associated with vaccination and hospitalization by a 3.2-fold risk ratio for the 95% CI to encompass the null value after controlling sex, asthma, immunocompromising condition, age, health region, and case onset date. The E-value for the adjusted association between vaccination with one dose and the risk of hospitalization among the pediatric population (i.e., aOR = 0.21) was 9.0 and the E-value for the upper 95% CI bound (i.e., 95% CI upper bound = 0.77) was 1.9.

Our vaccine effectiveness estimates did not substantially change when we matched on specimen collection date, when we only included symptomatic cases, and when we included income and visible minority status at the census subdivision level (S2 Table). In our sensitivity analysis where only cases with onset dates prior to January 5, 2022 were included, the association between one vaccine dose and hospitalization among pediatric Omicron SARS-CoV-2 cases was invalid (due to too few hospitalizations in the time period), but the association among adolescent Omicron SARS-CoV-2 cases remained unchanged compared to our main analysis (S2 Table).

Discussion

In this population-based analysis from the Canadian province of Ontario, we find evidence that mRNA vaccines against SARS-CoV-2 prevent hospitalization when adolescents receive a full two dose series and children received one dose in the context of an epidemic dominated by the immune-evasive Omicron variants. Among adolescent SARS-CoV-2 cases one vaccine dose was not significantly associated with a lower likelihood of protection against hospitalization in the pre-Omicron and Omicron periods, and two doses were not significantly protective in the pre-Omicron period. Our vaccine effectiveness estimates demonstrate robustness to unmeasured confounding, as evidenced by the high E-values for the estimates. As all subjects in this study had SARS-CoV-2 infection, the protection afforded by vaccination against hospitalization relates to attenuation of disease and is independent of protection provided against infection. This may in part explain our finding that vaccination was ineffective at preventing hospitalization in adolescents conditional upon infection with the Delta variant. SARS-CoV-2 infection with the Delta variant was less common among vaccinated individuals, and our study cohort during the Delta time period likely had a higher proportion of high risk population members compared to during the Omicron period [33].

The studies that have focused on SARS-CoV-2 vaccine effectiveness against hospitalization in adolescent and pediatric populations are not directly comparable to ours due to differences in control selection. Our estimated vaccine effectiveness estimates against hospitalization are lower than most estimates from studies using negative controls in the United States [810], likely because we isolated the impact of vaccination of hospitalization risk, independent of susceptibility to SARS-CoV-2 infection. In a case-control study conducted between May and October 2021 (i.e., prior to widespread infection with Omicron) with test-negative and syndrome-negative controls, vaccine effectiveness against hospitalization was 94% among those with two doses and 97% among those with one dose [8]. A study with 164 hospitals in the United States from April 2021 to January 2022 found that two dose vaccine effectiveness against hospitalization was 74% among children ages 5–11 years and between 73% and 94% among adolescents [10].

In contrast, two dose vaccine effectiveness was calculated to be 40% against hospitalization among adolescents during the circulation of Omicron in a study including 31 hospitals across the United States [11]. This lower estimate but may be due to earlier adolescent vaccination dates in the United States compared to in Ontario, as well as the impact of vaccine waning [11]. Unlike prior studies, our control group consisted of non-hospitalized SARS-CoV-2 cases, not uninfected individuals. This allowed for the isolation of the effectiveness of vaccination against hospitalization independent of infection risk, which explains our lower vaccine effectiveness estimates against hospitalization compared to the majority of prior studies [12].

Our study has several notable strengths. Due to Ontario’s high quality public health surveillance data, we isolated the impact of vaccination on hospitalization risk with control for individual level demographic and health related factors. Our study includes the population of a diverse region with publicly funded healthcare [34, 35]. We used a quantitative bias analysis to quantify the susceptibility of our results to unmeasured confounding [29]. Our study also had a few limitations. First, we did not consider time since vaccination in our analysis. In a recent study in the United States, time since vaccination was not shown to significantly impact SARS-CoV-2 vaccine effectiveness estimates against hospitalization in younger individuals [10]. We were unable to assess two dose effectiveness against hospitalization specifically among cases ages 4–11 years due to few hospitalizations among two-dose vaccinated individuals. Additionally, with only 62 total hospitalizations included in our study, we were likely underpowered to detect all true vaccine effects. Vaccination may be protective where our study showed a lack of statistical significance (i.e., Type II error), and a larger study over a longer period is needed. Some individuals aged four in our study were ineligible for vaccination. If they had a higher incidence of hospitalization, our vaccine effectiveness estimates would be biased away from the null. The availability of testing, and the propensity to get tested may have differed between vaccinated and unvaccinated individuals, and the comorbidities may have been more completely ascertained among hospitalized SARS-CoV-2 cases. Finally, comorbidities (i.e., asthma and immunocompromising conditions) in the CCM did not have a standardized definition which leads to residual confounding.

With the continued emergence of variants that may further decrease SARS-CoV-2 vaccine effectiveness against infection, it is vital to consider how effective these vaccines are at preventing severe outcomes [36]. Future studies should examine the impact of subsequent vaccine doses in adolescent and pediatric populations on SARS-CoV-2 hospitalization risk with emerging variants.

Conclusions

In this evaluation of the effectiveness of SARS-CoV-2 vaccination in adolescent and pediatric Omicron cases ages 4 to 17 in Ontario, Canada, we found that vaccination with two doses among adolescents and one dose among children is associated with a decreased likelihood of hospitalization, even when the vaccines do not prevent infection. SARS-CoV-2 vaccines remain an effective intervention to prevent severe outcomes in adolescent and pediatric populations, and continued efforts are needed to increase vaccine uptake in these populations.

Supporting information

S1 Fig. Directed Acyclic Graph (DAG).

DAG of the relationship between vaccination status and hospitalization among adolescent and pediatric SARS-CoV-2 cases.

https://doi.org/10.1371/journal.pone.0283715.s001

(TIF)

S1 Table. Description of adolescent and pediatric SARS-CoV-2 cases by hospitalization status, stratified by age and time period (n = 27,736).

https://doi.org/10.1371/journal.pone.0283715.s002

(DOCX)

S2 Table. Adjusted odds ratios of the association between vaccination status and hospitalization among adolescent and pediatric SARS-CoV-2 cases with varying assumptions.

https://doi.org/10.1371/journal.pone.0283715.s003

(DOCX)

Acknowledgments

We wish 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. We are also appreciative to the reviewers for providing comprehensive and constructive reviews.

References

  1. 1. World Health Organization. WHO Coronavirus (COVID-19) Dashboard. 2021 Mar 19 [cited 2022 March 20]. Available from: https://covid19.who.int.
  2. 2. 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. Open Forum Infect Dis. 2022. pmid:35611346
  3. 3. Buchan SA, Chung H, Brown KA, Austin PC, Fell DB, Gubbay JB, et al. Estimated Effectiveness of COVID-19 Vaccines Against Omicron or Delta Symptomatic Infection and Severe Outcomes. JAMA Network Open. 2022;5(9):e2232760–e. pmid:36136332
  4. 4. Ontario Agency for Health Protection and Promotion (Public Health Ontario). COVID-19 Public Health Measures Related to the Delta Variant. 2021 Jul 12 [cited 2022 March 20]. Available from: https://www.publichealthontario.ca/-/media/documents/ncov/voc/2021/07/covid-19-public-health-measures-delta-variant.pdf?sc_lang=en.
  5. 5. Ontario Agency for Health Protection and Promotion (Public Health Ontario). COVID-19 Variant of Concern Omicron (B.1.1.529): Risk Assessment, February 23, 2022. 2022 Feb 25 [cited 2022 March 15]. Available from: https://www.publichealthontario.ca/-/media/Documents/nCoV/voc/2022/02/covid-19-omicron-b11529-risk-assessment-feb-23.pdf?sc_lang=en.
  6. 6. Fowlkes AL, Yoon SK, Lutrick K, Gwynn L, Burns J, Grant L, et al. Effectiveness of 2-Dose BNT162b2 (Pfizer BioNTech) mRNA Vaccine in Preventing SARS-CoV-2 Infection Among Children Aged 5–11 Years and Adolescents Aged 12–15 Years—PROTECT Cohort, July 2021–February 2022. MMWR Morb Mortal Wkly Rep. 2022;ePub. pmid:35298453
  7. 7. World Health Organization. WHO, UN set out steps to meet world COVID vaccination targets. 2021 Oct 7 [cited 2022 March 15]. Available from: https://www.who.int/news/item/07-10-2021-who-un-set-out-steps-to-meet-world-covid-vaccination-targets.
  8. 8. Olson SM, Newhams MM, Halasa NB, Price AM, Boom JA, Sahni LC, et al. Effectiveness of BNT162b2 Vaccine against Critical Covid-19 in Adolescents. N Engl J Med. 2022;386(8):713–23. pmid:35021004
  9. 9. Dorabawila V, Hoefer D, Bauer UE, Bassett MT, Lutterloh E, Rosenberg ES. Effectiveness of the BNT162b2 vaccine among children 5–11 and 12–17 years in New York after the Emergence of the Omicron Variant. medRxiv. 2022:2022.02.25.22271454.
  10. 10. Klein NP, Stockwell MS, Demarco M, Gaglani M, Kharbanda AB, Irving SA, et al. Effectiveness of COVID-19 Pfizer-BioNTech BNT162b2 mRNA Vaccination in Preventing COVID-19–Associated Emergency Department and Urgent Care Encounters and Hospitalizations Among Nonimmunocompromised Children and Adolescents Aged 5–17 Years—VISION Network, 10 States, April 2021–January 2022. MMWR Morb Mortal Wkly Rep. 2022;71:352–8. pmid:35239634
  11. 11. Price AM, Olson SM, Newhams MM, Halasa NB, Boom JA, Sahni LC, et al. BNT162b2 Protection against the Omicron Variant in Children and Adolescents. N Engl J Med. 2022. pmid:35353976
  12. 12. Chua H, Feng S, Lewnard JA, Sullivan SG, Blyth CC, Lipsitch M, et al. The Use of Test-negative Controls to Monitor Vaccine Effectiveness: A Systematic Review of Methodology. Epidemiology. 2020;31(1):43–64. pmid:31609860
  13. 13. Hudgens MG, Halloran ME. Causal Vaccine Effects on Binary Postinfection Outcomes. J Am Stat Assoc. 2006;101(473):51–64. pmid:19096723
  14. 14. Statistics Canada. Table 17-10-0134-01 Estimates of population (2016 Census and administrative data), by age group and sex for July 1st, Canada, provinces, territories, health regions (2018 boundaries) and peer groups. 2022.
  15. 15. Public Health Ontario. COVID-19 Vaccine Uptake in Ontario: December 14, 2020 to January 30, 2022. 2022 Feb 4 [cited 2022 March 15]. Available from: https://www.publichealthontario.ca/-/media/documents/ncov/epi/covid-19-vaccine-uptake-ontario-epi-summary.pdf?la=en.
  16. 16. Ontario Ministry of Health. COVID-19 Vaccine Booking Expanding to Youth 12+ Ahead of Schedule. 2021 May 21 [cited 2022 Mar 15]. Available from: https://news.ontario.ca/en/release/1000185/covid-19-vaccine-booking-expanding-to-youth-12-ahead-of-schedule.
  17. 17. Ontario Ministry of Health. COVID-19 Vaccine Bookings to Open For All Children Aged Five to 11. 2021 Nov 22 [cited 2022 March 15]. Available from: https://news.ontario.ca/en/release/1001195/covid-19-vaccine-bookings-to-open-for-all-children-aged-five-to-11.
  18. 18. Health Canada. Health Canada authorizes use of Moderna COVID-19 vaccine in children 6 months to 5 years of age 2022 [cited 2022 Dec 8]. Available from: https://www.canada.ca/en/health-canada/news/2022/07/health-canada-authorizes-use-of-moderna-covid-19-vaccine-in-children-6-months-to-5-years-of-age.html.
  19. 19. Statistics Canada. Data products, 2016 Census. 2021 [cited 2022 Mar 15]. Available from: https://www12.statcan.gc.ca/census-recensement/2016/dp-pd/index-eng.cfm.
  20. 20. Ontario Ministry of Health. Ontario Updating Public Health Measures and Guidance in Response to Omicron. 2021 Dec 30 [cited 2022 March 15]. Available from: https://news.ontario.ca/en/release/1001386/ontario-updating-public-health-measures-and-guidance-in-response-to-omicron.
  21. 21. Public Health Agency of Canada. Summary of the National Advisory Committee on Immunization (NACI) Statement of November 19, 2021: Recommendation on the use of the Pfizer- BioNTech COVID-19 vaccine (10 mcg) in children 5–11 years of age. 2021.
  22. 22. Graff K, Smith C, Silveira L, Jung S, Curran-Hays S, Jarjour J, et al. Risk Factors for Severe COVID-19 in Children. Pediatr Infect Dis J. 2021;40(4). pmid:33538539
  23. 23. Greenland S, Pearl J, Robins JM. Causal diagrams for epidemiologic research. Epidemiology. 1999;10(1):37–48. pmid:9888278
  24. 24. Ontario Health. Ontario Health Regions. 2022 Sep 27 [cited 2022 October 3]. Available from: https://www.ontariohealth.ca/about-us/our-programs/ontario-health-regions.
  25. 25. Cohen J. Statistical Power Analysis for the Behavioral Sciences. 2nd ed. Hillsdale, NJ: Erlbaum; 1988.
  26. 26. Yang D, Dalton JE. A unified approach to measuring the effect size between two groups using SAS. SAS Global Forum, 2012.
  27. 27. Austin PC. Using the Standardized Difference to Compare the Prevalence of a Binary Variable Between Two Groups in Observational Research. Commun Stat Simul Comput. 2009;38(6):1228–34.
  28. 28. Public Health Ontario. Early Dynamics of Omicron in Ontario, November 1 to December 23, 2021. 2022 Jan [cited 2022 Oct 5]. Available from: https://www.publichealthontario.ca/-/media/Documents/nCoV/epi/covid-19-early-dynamics-omicron-ontario-epi-summary.pdf?sc_lang=en.
  29. 29. VanderWeele TJ, Ding P. Sensitivity analysis in observational research: introducing the E-value. Ann Intern Med. 2017;167(4):268–74. pmid:28693043
  30. 30. Linden A, Mathur MB, VanderWeele TJ. Conducting sensitivity analysis for unmeasured confounding in observational studies using E-values: The evalue package. Stata J. 2020;20(1):162–75.
  31. 31. Ioannidis JPA, Tan YJ, Blum MR. Limitations and Misinterpretations of E-Values for Sensitivity Analyses of Observational Studies. Ann Intern Med. 2019;170(2):108–11. pmid:30597486
  32. 32. VanderWeele TJ, Mathur MB, Ding P. Correcting Misinterpretations of the E-Value. Ann Intern Med. 2019;170(2):131–2. pmid:30597489
  33. 33. Juthani PV, Gupta A, Borges KA, Price CC, Lee AI, Won CH, et al. Hospitalisation among vaccine breakthrough COVID-19 infections. The Lancet Infectious Diseases. 2021;21(11):1485–6. pmid:34506735
  34. 34. Ontario Ministry of Health and Ontario Ministry of Long-Term Care. Ontario Health Insurance Plan. 2022 Jul 06 [cited 2022 Oct 03]. Available from: https://www.health.gov.on.ca/en/pro/programs/ohip/sob/.
  35. 35. Ontario Ministry of Finance. Fact Sheet 9: Ethnic origin and visible minorities. 2022 Apr 08 [cited 2022 October 05]. Available from: https://www.ontario.ca/document/2016-census-highlights/fact-sheet-9-ethnic-origin-and-visible-minorities.
  36. 36. Callaway E. Why does the Omicron sub-variant spread faster than the original? Nature. 2022;602(7898):556–7. Epub 2022/02/19. pmid:35177843