Skip to main content
Advertisement
  • Loading metrics

Booster vaccination protection against SARS-CoV-2 infections in young adults during an Omicron BA.1-predominant period: A retrospective cohort study

  • Jiayue Wan,

    Roles Formal analysis, Investigation, Methodology, Writing – original draft, Writing – review & editing

    Affiliation School of Operations Research and Information Engineering, Cornell University College of Engineering, Ithaca, New York, United States of America

  • Casey L. Cazer ,

    Contributed equally to this work with: Casey L. Cazer, Marin E. Clarkberg, Shane G. Henderson, Scarlett E. Lee, Genevive R. Meredith, Marwan Osman, David B. Shmoys

    Roles Conceptualization, Formal analysis, Funding acquisition, Investigation, Methodology, Supervision, Validation, Visualization, Writing – review & editing

    Affiliations Department of Public and Ecosystem Health, Cornell University College of Veterinary Medicine, Ithaca, New York, United States of America, Department of Population Medicine and Diagnostic Sciences, Cornell University College of Veterinary Medicine, Ithaca, New York, United States of America

  • Marin E. Clarkberg ,

    Contributed equally to this work with: Casey L. Cazer, Marin E. Clarkberg, Shane G. Henderson, Scarlett E. Lee, Genevive R. Meredith, Marwan Osman, David B. Shmoys

    Roles Conceptualization, Data curation, Investigation, Methodology, Resources, Software, Validation, Writing – review & editing

    Affiliation Institutional Research and Planning, Cornell University, Ithaca, New York, United States of America

  • Shane G. Henderson ,

    Contributed equally to this work with: Casey L. Cazer, Marin E. Clarkberg, Shane G. Henderson, Scarlett E. Lee, Genevive R. Meredith, Marwan Osman, David B. Shmoys

    Roles Formal analysis, Funding acquisition, Methodology, Validation, Writing – review & editing

    Affiliation School of Operations Research and Information Engineering, Cornell University College of Engineering, Ithaca, New York, United States of America

  • Scarlett E. Lee ,

    Contributed equally to this work with: Casey L. Cazer, Marin E. Clarkberg, Shane G. Henderson, Scarlett E. Lee, Genevive R. Meredith, Marwan Osman, David B. Shmoys

    Roles Methodology, Validation, Writing – review & editing

    Affiliation Department of Microbiology and Immunology, Cornell University College of Veterinary Medicine, Ithaca, New York, United States of America

  • Genevive R. Meredith ,

    Contributed equally to this work with: Casey L. Cazer, Marin E. Clarkberg, Shane G. Henderson, Scarlett E. Lee, Genevive R. Meredith, Marwan Osman, David B. Shmoys

    Roles Conceptualization, Investigation, Project administration, Writing – review & editing

    grm79@cornell.edu

    Affiliations Department of Public and Ecosystem Health, Cornell University College of Veterinary Medicine, Ithaca, New York, United States of America, Master of Public Health Program, Cornell University, Ithaca, New York, United States of America

  • Marwan Osman ,

    Contributed equally to this work with: Casey L. Cazer, Marin E. Clarkberg, Shane G. Henderson, Scarlett E. Lee, Genevive R. Meredith, Marwan Osman, David B. Shmoys

    Roles Formal analysis, Visualization, Writing – original draft, Writing – review & editing

    Affiliations Department of Public and Ecosystem Health, Cornell University College of Veterinary Medicine, Ithaca, New York, United States of America, Cornell Atkinson Center for Sustainability, Cornell University, Ithaca, New York, United States of America

  • David B. Shmoys ,

    Contributed equally to this work with: Casey L. Cazer, Marin E. Clarkberg, Shane G. Henderson, Scarlett E. Lee, Genevive R. Meredith, Marwan Osman, David B. Shmoys

    Roles Methodology, Resources, Validation, Writing – review & editing

    Affiliation School of Operations Research and Information Engineering, Cornell University College of Engineering, Ithaca, New York, United States of America

  • Peter I. Frazier

    Roles Funding acquisition, Methodology, Resources, Software, Supervision, Validation, Writing – review & editing

    Affiliation School of Operations Research and Information Engineering, Cornell University College of Engineering, Ithaca, New York, United States of America

Abstract

Background

While booster vaccinations clearly reduce the risk of severe Coronavirus Disease 2019 (COVID-19) and death, the impact of boosters on Severe Acute Respiratory Syndrome Coronavirus 2 (SARS-CoV-2) infections has not been fully characterized: Doing so requires understanding their impact on asymptomatic and mildly symptomatic infections that often go unreported but nevertheless play an important role in spreading SARS-CoV-2. We sought to estimate the impact of COVID-19 booster doses on SARS-CoV-2 infections in a vaccinated population of young adults during an Omicron BA.1-predominant period.

Methods and findings

We implemented a cohort study of young adults in a college environment (Cornell University’s Ithaca campus) from a period when Omicron BA.1 was the predominant SARS-CoV-2 variant on campus (December 5 to December 31, 2021). Participants included 15,800 university students who completed initial vaccination series with vaccines approved by the World Health Organization for emergency use, were enrolled in mandatory at-least-weekly surveillance polymerase chain reaction (PCR) testing, and had no positive SARS-CoV-2 PCR test within 90 days before the start of the study period. Robust multivariable Poisson regression with the main outcome of a positive SARS-CoV-2 PCR test was performed to compare those who completed their initial vaccination series and a booster dose to those without a booster dose.

A total of 1,926 unique SARS-CoV-2 infections were identified in the study population. Controlling for sex, student group membership, date of completion of initial vaccination series, initial vaccine type, and temporal effect during the study period, our analysis estimates that receiving a booster dose further reduces the rate of having a PCR-detected SARS-CoV-2 infection relative to an initial vaccination series by 56% (95% confidence interval [42%, 67%], P < 0.001). While most individuals had recent booster administration before or during the study period (a limitation of our study), this result is robust to the assumed delay over which a booster dose becomes effective (varied from 1 day to 14 days). The mandatory active surveillance approach used in this study, under which 86% of the person-days in the study occurred, reduces the likelihood of outcome misclassification. Key limitations of our methodology are that we did not have an a priori protocol or statistical analysis plan because the analysis was initially done for institutional research purposes, and some analysis choices were made after observing the data.

Conclusions

We observed that boosters are effective, relative to completion of initial vaccination series, in further reducing the rate of SARS-CoV-2 infections in a college student population during a period when Omicron BA.1 was predominant; booster vaccinations for this age group may play an important role in reducing incidence of COVID-19.

Author summary

Why was this study done?

  • In late 2021, government and health organizations were developing recommendations for Severe Acute Respiratory Syndrome Coronavirus 2 (SARS-CoV-2) booster vaccination for at-risk groups and the general public.
  • University students are a generally healthy population, but high rates of SARS-CoV-2 infections occur on college campuses.
  • Understanding whether booster vaccination reduces SARS-CoV-2 infections, specifically the Omicron variant, among university students will help universities and others decide whether to recommend or require booster vaccination in this population to prevent spread of disease to the local community and disruptions to the school year.

What did the researchers do and find?

  • We collected data on SARS-CoV-2 infections and vaccination status of over 15,000 students enrolled in mandatory testing at one university.
  • We created a statistical model to compare SARS-CoV-2 infection rates among university students who had received a booster vaccination prior to an Omicron variant outbreak and students who had completed an initial vaccination series but did not receive a booster vaccination.
  • The statistical model accounted for differences in infection risk due to participation in social organizations and athletic teams, sex, time during the study period, and initial vaccination date.
  • Students who received a booster were approximately half as likely to get infected as students who had completed their initial vaccination series but not received a booster.

What do these findings mean?

  • Booster vaccination may be a useful tool to reduce SARS-CoV-2 infections on university campuses.
  • Our analysis was initially performed as part of an institutional response and thus did not have an a priori protocol or statistical analysis plan, which should be considered in interpreting our results. Additional studies are needed to confirm these findings in the context of new SARS-CoV-2 variants and for new booster formulations.

Introduction

Coronavirus Disease 2019 (COVID-19) vaccines reduce Severe Acute Respiratory Syndrome Coronavirus 2 (SARS-CoV-2) infections and symptom severity [1], yet breakthrough infections occur, especially with the Omicron (B.1.1.529) variant [2,3]. The effectiveness of Food and Drug Administration (FDA)-authorized or approved vaccines BNT162b2, mRNA-1273, and Ad26.COV2.S in preventing SARS-CoV-2 infections has dropped dramatically due to immune evasion and waning of vaccine-induced immunity over time [47]. The Omicron variant exhibits immune system escape as the result of several mutations [8,9]; this, and the high transmissibility of the Omicron variant, are leading to higher infection rates, strain on healthcare systems, and increased mortality [10,11]. With the emergence of new variants combined with waning immunity, the CDC recommends a booster 6 months after an initial mRNA vaccine series or 2 months after Ad26.COV2.S vaccination to prevent symptomatic and severe outcomes of COVID-19 [7,12]. The booster dose elicits an increase in antibody neutralization titers against the Omicron variant and causes affinity maturation leading to a better antibody response, maintaining long-term protection against severe COVID-19 outcomes [1215].

While boosters are understood to be effective against severe disease, hospitalization [16], and symptomatic infections [12,17] resulting from the Omicron variant, limited information is available about their effectiveness against asymptomatic and mild symptomatic infections that may go unreported in the absence of asymptomatic surveillance testing. Asymptomatic and mild symptomatic infections play an important role in spreading SARS-CoV-2 [18]. Moreover, the clinical course of COVID-19 is understood to vary by age [19]; current estimates of booster effectiveness based on the general population may not apply to cohorts whose age distribution differs substantially from that of the general population.

We leveraged polymerase chain reaction (PCR) testing data derived from a SARS-CoV-2 surveillance program that required mandatory routine SARS-CoV-2 testing of students (through December 14, 2021) and departure tests before students left the Ithaca campus. During December 5 to December 31, 2021, there was an explosive increase in Omicron cases to the point where the Omicron variant (BA.1) became predominant at Cornell University’s Ithaca campus [20], and 1,926 SARS-CoV-2 infections were identified in a vaccinated population of students. We sought to estimate the effectiveness of COVID-19 boosters in reducing SARS-CoV-2 infections in the vaccinated population of students during this Omicron-predominant period.

Methods

We used a retrospective cohort study to estimate the effectiveness of COVID-19 boosters in reducing SARS-CoV-2 infections. Specific to our context was a highly vaccinated study population, and a period where Omicron BA.1 was the dominant viral strain. We did not have an a priori protocol or statistical analysis plan. As detailed in Section A in S1 Appendix, our initial analysis [21] came as part of Cornell University’s institutional response to COVID-19 during an Omicron-driven outbreak among a highly vaccinated population. This plan was made after observing aggregate infection rates that compared those who had uploaded early proof of booster vaccination (less than 17% of those in the study population who would eventually provide documentation) against those who had not. We subsequently developed a study plan that modified two of the covariates in this analysis (adding sex and changing how primary vaccine series completion date was coded) because we hypothesized that they are confounders. We then made additional modifications in response to reviewer requests. As described in Section A in S1 Appendix, estimates of vaccine booster effectiveness were not sensitive to these changes in the analysis.

We utilized deidentified student data from Cornell University’s COVID-19 surveillance database. We included PCR-positive cases identified at Cornell surveillance testing sites, the campus student health service, and a sampling site operated by the local hospital system during the period of December 5, 2021 to December 31, 2021 (referred to as the “Omicron outbreak” or the study period), in which the Omicron variant (BA.1) was the predominant strain. Samples collected at Cornell surveillance sites and the local hospital’s sampling site were tested using EZ-SARS-CoV-2 Real-Time RT-PCR [22] within a two-stage Dorfman procedure [23] using pools of size 5. The assay targets two highly conserved regions of the SARS-CoV-2 N gene, and the Omicron variant is not expected to significantly affect test performance. Samples collected at the campus student health service were tested using the Cepheid Xpert Xpress SARS-CoV-2 Assay [24]. The FDA reports that the introduction of the Omicron variant does not appear to have resulted in a significant change in the test performance of this assay [25]. This work was completed as a part of Cornell’s institutional planning and preparedness, designated as exempt from IRB review by the Cornell IRB.

Study context

In July 2020, as a part of its reopening strategy, Cornell implemented a robust COVID-19 surveillance testing program for students, faculty, and staff (1,637,394 PCR tests performed as of February 2, 2022). At the start of December 2021, 23,389 students were enrolled in full-time academic programs at Cornell’s Ithaca campus. Masks were required in all on-campus buildings. COVID-19 vaccination was required for all students (97% vaccinated), though medical or religious exemptions were granted. Boosters were encouraged (October 21, 2021) before being required (January 31, 2022). Through December 14, 2021, all undergraduates and professional program (veterinary, business, and law) students (n = 17,017) were enrolled in at-least-weekly mandatory PCR surveillance (compliance: 99.8% for surveillance tests scheduled between November 29, 2021 and December 14, 2021) using an anterior nasal swab sample. Approximately 86% of person-days included in the study (see the Study population section below) fell within this period. Students leaving the Ithaca campus before the end of the study period (82% of students in the study population indicated that they left before the end of the study period through a Cornell-managed checklist) were required to seek an additional test shortly before departure. Additionally, students were prompted to seek a PCR test if they were identified as a close contact of a case or if they developed COVID-like symptoms. Free PCR testing was available to any constituent of the university 7 days per week. Whole genome sequencing confirmed the presence of the Omicron variant (BA.1) in Cornell community COVID-19 surveillance samples collected on December 1, 2021; by December 11, 2021, Omicron was the predominant variant on the Ithaca campus [20].

Study population

Our study included 15,800 Cornell students. Study inclusion criteria were as follows: being an active student at Cornell’s Ithaca campus enrolled in mandatory surveillance testing, being vaccinated with a vaccine approved by the World Health Organization (WHO) for emergency use at the time of the study (2 doses of BNT162b2, 2 doses of mRNA-1273, 1 dose of Ad26.COV2.S, or 2 doses of other WHO-approved vaccines), and having no positive SARS-CoV-2 PCR test within 90 days before the start of the study period (Fig 1). Students were required to upload a photo of a vaccine card documenting their booster status via an online form or request an exemption by February 1, 2022 [26]. According to vaccination records that were uploaded and validated by February 8, 2022, 11.9% of the study population (n = 1,876) received a booster before the study period, with the earliest booster dose administered on August 3, 2021. Our analysis does not distinguish the booster vaccine manufacturer, although this information is available in Table B in S1 Appendix. Students were excluded from the study population if they were not subject to mandatory surveillance testing, i.e., graduate students (n = 6,372), tested PCR-positive within 90 days before the start of the study period (n = 214), did not have at least one PCR surveillance test between November 5, 2021 and December 4, 2021 (n = 130), had not completed a WHO-approved initial vaccine series (n = 701), had an unspecified sex (n = 9), or had invalid vaccination records (n = 163). Students with no PCR test records during the study period (n = 285) were considered lost to follow-up.

thumbnail
Fig 1. Enrollment of the study population.

Students were included if they were enrolled in a full-time academic program based at Cornell’s Ithaca campus in Fall 2021, were subject to mandatory surveillance testing, had no positive SARS-CoV-2 PCR test within 90 days before the start of the study period, had at least one PCR test between November 5, 2021 and December 4, 2021, had a valid vaccination record indicating completion of WHO-approved initial vaccination series at least 14 days before the study period, and specified their sex in university records. PCR, polymerase chain reaction; SARS-CoV-2, Severe Acute Respiratory Syndrome Coronavirus 2; WHO, World Health Organization.

https://doi.org/10.1371/journal.pmed.1004153.g001

Statistical analysis

The primary outcome of our study was PCR-confirmed SARS-CoV-2 infection. To estimate the effectiveness of COVID-19 boosters against SARS-CoV-2 infections during this Omicron-predominant period, we calculated the person-days that each student contributed to the boosted and nonboosted population during the study period. The total number of person-days contributed by a student is the number of days between December 5, 2021 and either their final test date in the study period or their first PCR-positive test date, whichever comes first. Of the study population, 38,023 PCR tests were performed during the study period. Approximately 98.2% of the students (n = 15,515) had at least one PCR test (mean = 2.4 tests, median = 2.0 tests, SD = 1.2 tests), hence contributing at least one person-day (Fig A in S1 Appendix). We assume that booster vaccinations become effective 7 days after administration, based on results from a study of the Israeli general population [27]. As a result, each person-day has its booster status labeled as boosted (1) or control (0) based on whether the associated individual received their booster dose at least 7 days before that day (Fig B in S1 Appendix).

We performed a multivariable Poisson regression to estimate the effect of receiving a booster dose on having a positive PCR-based diagnosis during the study period. The unit of observation is an exposure period of an individual (measured in person days). An offset equal to log(person-days), i.e., the natural logarithm of the number of days in the exposure period, is added in the regression model for each exposure period. Each exposure period lasts until the end of the week (Week 0 is December 5 to December 11, etc.), a change in booster status (in which case a week is split into two exposure periods), or the person completes their last test in the study, whichever comes first. This model can be seen as a Cox proportional hazards analysis [28], in which it is assumed that exposure risk accrues at a constant rate within each exposure period. To account for possible within-individual correlation among observations over time (i.e., the risk of some individual contracting the disease is likely persistent over time due to individual behavior), we used generalized estimating equations (GEEs) [29,30] with an exchangeable correlation structure over weeks in the study period to obtain consistent coefficient estimates and robust variance estimates for the regression model.

We controlled for sex, student group (undergraduate or professional), fraternity/sorority participation (yes or no), varsity athletic team participation (yes or no), initial vaccine type (BNT162b2, mRNA-1273, Ad26.COV2.S, or other WHO-approved vaccines), initial vaccination series completion date (date of receiving 1 dose of Ad26.COV2.S or 2 doses of other WHO-approved vaccines, grouped by months from January 2021 to November 2021), and temporal effect during the study period (grouped by week). We also created univariable models for each of these covariates to obtain unadjusted incidence rate ratios. We included student group, fraternity/sorority participation, and athletic team participation as covariates because case investigation of data from before the study period suggested that these covariates explained much of the heterogeneity in the risk of infection across students [31]. Undergraduate students were divided into 3 subgroups based on fraternity/sorority and athletic team participation. Students who were in both fraternities/sororities and athletic teams were classified in the fraternity/sorority group. Thus, we had six student groups: (1) undergraduate fraternity/sorority participants; (2) undergraduate athletes not in Group 1; (3) undergraduates not in Groups 1 or 2; (4) law students; (5) postbaccalaureate business students; (6) veterinary students. Initial vaccination series completion date was included as a covariate to adjust for heterogeneous social behavior, inclination to receive vaccination, and waning vaccine immune response. The week during the study period was included as a categorical covariate to control for time-varying prevalence and the resulting risk of exposure. This mitigates bias that would have otherwise been created by temporal variation in the fraction of on-campus students boosted. We used a categorical covariate rather than a continuous one because prevalence did not change linearly with time. We did not include age as a covariate because there is little age variation in the study population (Fig C in S1 Appendix). The dependent variable was whether a student tests PCR-positive for COVID-19 in a particular exposure period. The regression model was where λ is the incidence rate, and greek letters on the right-hand side of the specification are the regression coefficients. The p-values for the estimated coefficients and the 95% confidence intervals for the adjusted incidence rate ratios were adjusted using the Bonferroni correction [32].

To estimate the effectiveness of a booster dose, we used

We performed several robustness checks. To assess the effect of students leaving campus at the end of the semester and the accordingly paused mandatory surveillance, we performed a secondary GEE Poisson regression with the same specification as above, but for a shorter study period (from December 5, 2021 to December 18, 2021, the last day of exams).

To assess the possibility of bias remaining in our analysis despite the included controls, e.g., bias due to days since booster dose administration, we performed a GEE Poisson regression with multiple classes for the booster status (unboosted, 0 to 6 days after booster administration, ≥7 days after booster dose administration). We hypothesized that booster vaccination has negligible effect 0 to 6 days after administration.

We further performed a sensitivity analysis on the delay to boosted status, allowing person-days to count as boosted only after this delay after administration had elapsed, varying this delay from the day after booster administration (day 1) to 14 days (Fig B in S1 Appendix). We emphasize that varying this parameter only affects the booster status of person days from people that were boosted within or shortly before the start of the study period. We hypothesized that the shorter delays would result in smaller but still statistically significant booster effectiveness estimates as person-days from 0 to 6 days after booster administration would dilute the effect observed in ≥7 days after booster administration group. We expected to see similar estimates of booster effectiveness for all delays ≥7 days. We emphasize that this analysis is not designed to comment on how quickly boosters become effective, but instead only on the robustness of our main findings to assumptions about this delay parameter.

We also performed a logistic regression with GEE to assess the robustness of our conclusions to model specification. The unit of observation is a person-day. We assumed that booster vaccinations become effective 7 days after administration. The dependent variable is whether a student tests PCR-positive for SARS-CoV-2 on a particular day. The regression model was where p is the probability of infection, and greek letters on the right-hand side of the specification are the regression coefficients. This logistic model can be seen as a form of discrete time interval–censored survival model, in which a person’s endpoint is the first day on which they tested positive, and the right-censoring date is the date of their last negative test. As with the Poisson regression model, information bias is possible because testing did not occur every day. We approximated the incidence rate ratio using the adjusted odds ratio (aOR) when estimating the booster effectiveness.

A subset analysis of undergraduate students living on campus was performed to assess the robustness of our results to additional kinds of student clustering, i.e., residence, through which virus could spread. The description of this analysis is included in Section B in S1 Appendix.

All statistical analyses were performed in Python (V3.7.11), using the statsmodels package (V0.13.2). All code is publicly available (https://github.com/jiayuewan/booster_effectiveness). The data are protected by an institutional data requirement at Cornell that covers the use of data collected during COVID-19 surveillance and are maintained by Cornell University’s Office of the Vice President for Research and Innovation. These can be requested from that entity.

Results

Boosted students were more likely to be female, professional students, and have an earlier initial vaccine series than unboosted individuals (Table 1). During the study period, a total of 1,926 PCR-positive cases were identified and reported out of 15,800 students included in this analysis (overall infection risk of 12.2%). None of the cases required hospitalization. The overall infection risk among students that received a booster dose before December 5, 2021 is 6.2% (117 PCR-positives out of 1,876 boosted students), whereas the infection risk among students that did not receive a booster before December 5, 2021 is 13.0% (1,809 PCR-positives out of 13,924 students). The cumulative incidence among students who received a booster before December 5, 2021 was approximately half the cumulative incidence among students who were unboosted at that time (Fig 2). Similarly, the incidence rate was 0.6 per 100 person-days (127 PCR-positive cases/19,842 person-days) among students who had received the booster dose ≥7 days earlier, compared to 1.3 per 100 person-days (1,799 PCR-positive cases/135,214 person-days) among those who had not received the booster dose or had received the booster dose <7 days earlier.

thumbnail
Fig 2. SARS-CoV-2 infection cumulative incidence rate (number of infections per person) and its 95% confidence interval during the study period, broken out by booster dose status.

https://doi.org/10.1371/journal.pmed.1004153.g002

thumbnail
Table 1. Distribution of person-day data by sex, student group, initial vaccination series completion date, initial vaccine type, and week in the study period.

https://doi.org/10.1371/journal.pmed.1004153.t001

The booster effectiveness estimate from the main GEE Poisson model is 56% (95% CI: [42%, 67%]). In addition to boosters, the incidence rate in our fitted GEE Poisson model depends on the student group and the date of completing the initial vaccine series (Table 2). The incidence rate was significantly lower among students who completed the initial vaccine series after May 1, 2021 compared with those completing the initial vaccine series earlier. During the study period, the incidence rate was significantly higher among undergraduate students participating in fraternity and sorority activities (aIRR 2.16 [1.84, 2.54]) or belonging to athletic teams (aIRR 2.02 [1.57, 2.60]), and among professional students at the business school (aIRR 1.64 [1.15, 2.32]), when compared with other undergraduates. Students vaccinated with Ad26.COV2.S had a higher risk of infection relative to other vaccines, but the difference was not statistically significant in our model.

thumbnail
Table 2. Summary of the full GEE Poisson regression model with covariates for sex, student group, initial vaccination series completion date, initial vaccine type, and week in the study period.

https://doi.org/10.1371/journal.pmed.1004153.t002

To study the robustness to population size variation as students left campus at the end of the semester, we performed a GEE Poisson regression (Table C in S1 Appendix) for a shorter period (December 5, 2021 to December 18, 2021). This yields a similar booster effectiveness of 57% (95% CI: [42%, 68%]), suggesting that our estimate of booster effectiveness from the full regression model is robust to the effect of students leaving campus at the end of the semester and the accordingly paused mandatory surveillance. As an additional robustness check, we performed a GEE Poisson regression (Table D in S1 Appendix) with multiple classes for the booster status of a person-day (unboosted, 0 to 6 days after booster administration, ≥7 days after booster dose administration). As expected, a booster status corresponding to 0 to 6 days after booster dose administration was not found to have a statistically significant effect on PCR-diagnosed infection (P > 0.99) relative to being unboosted.

We further performed a sensitivity analysis to assess the impact of the assumed delay between booster administration and booster effectiveness. Although 5,993 students received their booster dose during the study period (Fig D in S1 Appendix), 5,210 of them received it after leaving the active surveillance program at the end of the semester. These students did not contribute any person-days to the study after their last surveillance test, so they contributed only unboosted person-days. For the 783 students who received a booster dose while in active surveillance, we varied this delay from 1 day to 14 days and calculated the booster effectiveness against infection using the full model (Figs B and E in S1 Appendix). The estimated booster effectiveness is consistently above 48% and is relatively flat after 7 days, suggesting that our estimate of booster effectiveness is robust to variation in the time required to mount an effective immune response after booster vaccination. Similarly, the incidence in the boosted population does not vary substantially as we change the assumed delay between booster administration and effectiveness (Table E in S1 Appendix).

In a separate sensitivity analysis using logistic regression with GEE (Table F in S1 Appendix), the estimated booster effectiveness is 57% (95% CI: [43%, 68%]). The estimate from the logistic regression with GEE is close to the estimate from the full GEE Poisson regression model.

In a subset analysis (Section B in S1 Appendix) on the undergraduate population living in on-campus housing, we did not observe a significant change in the estimated booster effectiveness on the subpopulation whether or not an additional covariate for housing building was included in the regression. This subset analysis justifies the use of our regression model without controlling for housing type on the full study population, given that housing information is only available for the on-campus undergraduates.

Discussion

This study is one of the first community studies to quantify booster vaccine effectiveness from an actively surveilled population of young adults. Our study provides evidence that booster vaccinations significantly reduced infections in university settings during periods when the Omicron BA.1 variant was predominant, supporting the implementation of booster vaccination requirements to minimize community transmission.

In this retrospective analysis of SARS-CoV-2 tests performed at Cornell University’s Ithaca campus during a 27-day period when Omicron BA.1 was the predominant variant, the incidence of COVID-19 infections was approximately halved among participants vaccinated with a booster dose of a COVID-19 vaccine approved by WHO for emergency use when compared with vaccinated participants without a booster dose. The calculated booster effectiveness was 56% (95% CI: [42%, 67%]), which is slightly lower than a previously reported effectiveness against symptomatic COVID-like illness in adults (66%; 95% CI: [64%, 68%]) in the same time period [17]. Our estimate is similar to the reduction in cumulative incidence of reported infections associated with booster vaccination during an Omicron wave in Los Angeles County [33] and consistent with a pooled estimate of booster effectiveness (47%; 95% CI: [19%, 65%]) from a meta-analysis of studies completed during periods of Omicron predominance across the globe [34]. Importantly, our analysis includes both symptomatic and asymptomatic infections because the student population was under active surveillance, whereas none of the studies included in the meta-analysis examined a population under active surveillance [7,3540]. Interestingly, between two studies on the effectiveness of a second booster dose, i.e., fourth vaccine dose [41,42], one that used data collected during active surveillance [42] also found lower effectiveness than one that did not [41]. Students vaccinated with Ad26.COV2.S had a higher risk of infection relative to other vaccines, similar to other studies [43], though the difference was not statistically significant, likely because a small number of students had Ad26.COV2.S initial doses (Table 1, Table B in S1 Appendix). We did not find a booster status corresponding to 0 to 6 days after booster dose administration to have a statistically significant effect on PCR-diagnosed infection relative to being unboosted. This could be because of an insufficient immune response to the booster within 6 days of administration, consistent with previous findings [22], or because of a lack of statistical power resulting from a small number of person days in the 0-to-6-day window after booster administration.

We observed lower incidence rates among students with more recent initial vaccinations, consistent with waning protection against SARS-CoV-2 infections observed within a few months of completing the initial vaccination series [44]. Months in which fewer individuals completed their initial vaccine series did not have a statistically significant effect, consistent with fewer data points resulting in less statistical power. Varsity athletes and students in fraternities/sororities at institutions of higher education may have more social contact than other undergraduates; thus, consistent with previous studies, we found that this population may be at higher risk for SARS-CoV-2 infection [31,45,46]. As in other higher-education institutions [46], contact tracing at Cornell identified fraternity/sorority gatherings as significant spreading events. Local case investigation efforts also pointed to additional events, including post-Thanksgiving break travel and a series of end-of-semester gatherings, as contributing to Omicron spread. With the emergence of highly transmissible variants, travel and social gatherings may put students at increased risk of SARS-CoV-2 infection.

Our Poisson regression analysis is similar to Bar-On and colleagues’ [47], a previous study of vaccine booster effectiveness against COVID-19. See also Bar-on and colleagues’ study [41] for another similar approach to estimating vaccine booster effectiveness against COVID-19. There are two common alternative statistical approaches to estimating vaccine effectiveness. The first approach is the emulated trial design (e.g., [48]), which matches subjects with the same covariates in an observational study to emulate a randomized trial. This design offers more robustness against misspecification of the functional form by which booster status and infection risk depend on covariates, but requires a larger dataset to enable matching simultaneously on all covariates. We did not have a sufficiently large dataset to match on all the covariates included in our model. A second alternative approach is the test-negative case–control design [17], which controls bias by including those individuals who sought testing, and improvements that model the propensity to test [49]. Since testing was mandatory and compliance was high, it is not necessary to model the propensity to test in this study population.

Limitations and conclusion

Several limitations of our approach should be considered in interpreting our results. First, we did not have an a priori protocol or statistical analysis plan because the analysis was initially done for institutional research purposes, and some analysis choices were made after observing the data (see Section A in S1 Appendix). Thus, confidence intervals and hypothesis tests would be invalidated by correlation that could have been unintentionally introduced between the data and the choice of analysis plan. Second, a limitation of the Poisson regression is information bias due to the fact that individuals were not tested every day. Delays in reporting a positive test would cause the period of the delay to have a misclassified infection status, whether this period is within one exposure period or spans multiple exposure periods. Third, there may be additional confounding bias because, even after controlling for covariates, people who avoid risk could also be more likely to get a booster vaccination, which could result in the boosted group having decreased exposure to SARS-CoV-2. This would cause our regression analyses to overestimate vaccine effectiveness. Conversely, if people had engaged in more risky behavior in response to perceived protection offered by booster vaccination, our analyses would underestimate vaccine effectiveness. Fourth, our data do not allow distinguishing between booster doses and additional doses for immunocompromised individuals, or distinguishing between symptomatic and asymptomatic infections for confirmed PCR-positive cases. As mentioned earlier, not all samples were tested for S-gene dropout, and sequencing of every PCR-positive case was not performed. If boosters are more effective against the Delta variant than the Omicron variant, then Delta infections during the study period would lead our analysis to overestimate booster effectiveness against Omicron. Fifth, enough COVID-19 cases were observed to provide sufficient statistical power for moderately precise booster effectiveness estimates, but not to estimate how booster effectiveness varies with the manufacturer of the booster dose or the original vaccine. Most students were boosted with BNT162b2 in this study (Table B in S1 Appendix). Sixth, we did not include information on previous COVID-19 infections and had no information on the students’ underlying medical conditions or use of nonpharmaceutical interventions, although students were required to wear a mask (cloth, medical, or respirator) while inside campus buildings. This may lead to confounding bias if medical history differentially affected booster choice or if booster choice affected the use of nonpharmaceutical interventions. Seventh, gaps in time between a student’s last test and their departure from Ithaca (or the end of the study period, whichever came first) could lead to attrition bias. Eighth, although nearly 90% of positive samples tested for S-gene dropout during the study period exhibited this marker of the Omicron BA.1 variant [20], and sequencing confirmed the presence of the Omicron variant in all samples sequenced during this period, not all samples were tested for S-gene dropout, and sequencing of each PCR-positive case was not performed. Finally, not every individual who was boosted before or during the study period may have uploaded their vaccination record, which may have led to information bias towards the null.

In conclusion, our findings are consistent with the notion that booster vaccine doses, relative to being vaccinated with the initial series, are effective in reducing infections in young adults during Omicron-predominant periods, thereby helping universities and other institutions to remain open safely. In the 2 months subsequent to the period covered in this study, almost the entire Cornell student body, if eligible, received at least 1 booster dose.

Supporting information

S1 Appendix. Supplementary Appendix.

Section A. Analysis Timeline. Section B. Additional Subset Analysis. Table A. Summary of the GEE Poisson regression model for the on-campus undergraduate population with and without housing building as an additional covariate, and other covariates for sex, student group, initial vaccine series completion date, initial vaccine type, week during the study period. Table B. Distribution of booster dose type among students in the study population that received their booster dose before or during the study period, broken out by initial vaccination and booster dose type. Table C. Summary of the GEE Poisson regression model in the shortened study period (December 5, 2021 to December 18, 2021) with covariates for sex, student group, initial vaccine series completion date, initial vaccine type, and week during the study period. Table D. Summary of the GEE Poisson regression model with multiclass booster status (divided into unboosted, 0–6 days after booster dose administration, and ≥7 days after booster dose administration) and covariates for sex, student group, initial vaccine series completion date, initial vaccine type, and week during the study period. Table E. Number of unboosted and boosted person-days, PCR-positive cases, incidence rate with respect to different assumed delays for the booster to become effective after administration. Table F. Summary of the GEE logistic regression model with covariates for sex, student group, initial vaccine series completion date, initial vaccine type, and week during the study period. Fig A. Distribution of the number of person-days contributed by each student in the study population. Fig B. Allocation of person-days to the control and booster cohorts. Fig C. Age distribution of students in the study population. Fig D. Cumulative number of students receiving COVID-19 booster dose, over time. Fig E. Mean and 95% confidence interval for the booster effectiveness against infection during the Omicron predominance period, as we vary the assumed delay for the booster dose to become effective after booster administration

https://doi.org/10.1371/journal.pmed.1004153.s002

(DOCX)

Acknowledgments

We would like to thank the large number of people at Cornell University who supported the university’s response to COVID-19.

References

  1. 1. Moghadas SM, Vilches TN, Zhang K, Wells CR, Shoukat A, Singer BH, et al. The impact of vaccination on coronavirus disease 2019 (COVID-19) outbreaks in the United States. Clin Infect Dis. 2021;73(12):2257–2264. pmid:33515252
  2. 2. Brown CM, Vostok J, Johnson H, Burns M, Gharpure R, Sami S, et al. Outbreak of SARS-CoV-2 infections, including COVID-19 vaccine breakthrough infections, associated with large public gatherings—Barnstable County, Massachusetts, July 2021. MMWR Morb Mortal Wkly Rep. 2021;70(31):1059–1062. pmid:34351882
  3. 3. Gupta RK, Topol EJ. COVID-19 vaccine breakthrough infections. Science. 2021;374(6575):1561–1562. pmid:34941414
  4. 4. Cele S, Jackson L, Khoury DS, Khan K, Moyo-Gwete T, Tegally H, et al. Omicron extensively but incompletely escapes Pfizer BNT162b2 neutralization. Nature. 2022;602(7898):654–656. pmid:35016196
  5. 5. Bruxvoort KJ, Sy LS, Qian L, Ackerson BK, Luo Y, Lee GS, et al. Effectiveness of mRNA-1273 against delta, mu, and other emerging variants of SARS-CoV-2: test negative case-control study. BMJ. 2021;375:e068848. pmid:34911691
  6. 6. Tartof SY, Slezak JM, Fischer H, Hong V, Ackerson BK, Ranasinghe ON, et al. Effectiveness of mRNA BNT162b2 COVID-19 vaccine up to 6 months in a large integrated health system in the USA: a retrospective cohort study. Lancet. 2021;398(10309):1407–1416. pmid:34619098
  7. 7. Ferdinands JM, Rao S, Dixon BE, Mitchell PK, DeSilva MB, Irving SA, et al. Waning 2-Dose and 3-Dose Effectiveness of mRNA Vaccines Against COVID-19-Associated Emergency Department and Urgent Care Encounters and Hospitalizations Among Adults During Periods of Delta and Omicron Variant Predominance—VISION Network, 10 States, August 2021-January 2022. MMWR Morb Mortal Wkly Rep. 2022;71(7):255–263. pmid:35176007
  8. 8. Lu L, Mok BW-Y, Chen L-L, Chan JM-C, Tsang OT-Y, Lam BH-S, et al. Neutralization of SARS-CoV-2 Omicron variant by sera from BNT162b2 or Coronavac vaccine recipients. Clin Infect Dis. 2021. pmid:34915551
  9. 9. Wang L, Cheng G. Sequence analysis of the emerging SARS-CoV-2 variant Omicron in South Africa. J Med Virol. 2022;94(4):1728–1733. pmid:34897752
  10. 10. Han P, Su C, Zhang Y, Bai C, Zheng A, Qiao C, et al. Molecular insights into receptor binding of recent emerging SARS-CoV-2 variants. Nat Commun. 2021;12(1):6103. pmid:34671049
  11. 11. Karim SSA, Karim QA. Omicron SARS-CoV-2 variant: a new chapter in the COVID-19 pandemic. Lancet. 2021;398(10317):2126–2128. pmid:34871545
  12. 12. Thompson MG, Natarajan K, Irving SA, Rowley EA, Griggs EP, Gaglani M, et al. Effectiveness of a third dose of mRNA vaccines against COVID-19-associated emergency department and urgent care encounters and hospitalizations among adults during periods of Delta and Omicron variant predominance—VISION Network, 10 States, August 2021-January 2022. MMWR Morb Mortal Wkly Rep. 2022;71(4):139–145. pmid:35085224
  13. 13. Choi A, Koch M, Wu K, Chu L, Ma L, Hill A, et al. Safety and immunogenicity of SARS-CoV-2 variant mRNA vaccine boosters in healthy adults: an interim analysis. Nat Med. 2021;27(11):2025–2031. pmid:34526698
  14. 14. Bajema KL, Dahl RM, Evener SL, Prill MM, Rodriguez-Barradas MC, Marconi VC, et al. Comparative effectiveness and antibody responses to Moderna and Pfizer-BioNTech COVID-19 vaccines among hospitalized veterans—five veterans affairs medical centers, United States, February 1-September 30, 2021. MMWR Morb Mortal Wkly Rep. 2021;70(49):1700–1705. pmid:34882654
  15. 15. Barda N, Dagan N, Cohen C, Hernán 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–2100. pmid:34756184
  16. 16. Collie S, Champion J, Moultrie H, Bekker L-G, Gray G. Effectiveness of BNT162b2 Vaccine against Omicron Variant in South Africa. N Engl J Med. 2021;386(5):494–496. 429 pmid:34965358
  17. 17. Accorsi EK, Britton A, Fleming-Dutra KE, Smith ZR, Shang N, Derado G, et al. Association between 3 doses of mRNA COVID-19 vaccine and symptomatic infection caused by the SARS-CoV-2 Omicron and Delta variants. JAMA. 2022;327(7):639–651. pmid:35060999
  18. 18. Johansson MA, Quandelacy TM, Kada S, Prasad PV, Steele M, Brooks JT, et al. SARS-CoV-2 transmission from people without COVID-19 Symptoms. JAMA Netw Open. 2021;4(1):e2035057. pmid:33410879
  19. 19. Bhattacharyya RP, Hanage WP. Challenges in Inferring Intrinsic Severity of the SARS-CoV-2 Omicron Variant. N Engl J Med. 2022;386(7):e14. pmid:35108465
  20. 20. Meredith GR, Diel DG, Frazier PI, Henderson SG, Koretzky GA, Wan J, et al. Routine Surveillance and Vaccination on a University Campus During the Spread of the SARS-CoV-2 Omicron Variant. JAMA Netw Open. 2022;5(5):e2212906-e. pmid:35583871
  21. 21. Cornell COVID-19 Modeling Team. Booster Effectiveness Against Omicron. 2022. Available from: https://covid.cornell.edu/_assets/files/booster-effectiveness-against-omicron-jan-27.pdf.
  22. 22. Laverack M, Tallmadge RL, Venugopalan R, Cronk B, Zhang X, Rauh R, et al. Clinical evaluation of a multiplex real-time RT-PCR assay for detection of SARS-CoV-2 in individual and pooled upper respiratory tract samples. Arch Virol. 2021;166(9):2551–2561. pmid:34259914
  23. 23. Cleary B, Hay JA, Blumenstiel B, Harden M, Cipicchio M, Bezney J, et al. Using viral load and epidemic dynamics to optimize pooled testing in resource-constrained settings. Sci Transl Med. 2021;13(589):eabf1568. pmid:33619080
  24. 24. Loeffelholz MJ, Alland D, Butler-Wu SM, Pandey U, Perno CF, Nava A, et al. Multicenter Evaluation of the Cepheid Xpert Xpress SARS-CoV-2 Test. J Clin Microbiol. 2020;58(8):e00926–e00920. pmid:32366669
  25. 25. Food and Drug Administration. SARS-CoV-2 Viral Mutations: Impact on COVID-19 Tests 2021. Available from: https://www.fda.gov/medical-devices/coronavirus-covid-19-and-medical-devices/sars-cov-2-viral-mutations-impact-covid-19-tests.
  26. 26. Cornell University. Booster requirement for spring semester 2021. Available from: https://covid.cornell.edu/updates/20212112-booster-requirement.cfm.
  27. 27. Patalon T, Gazit S, Pitzer VE, Prunas O, Warren JL, Weinberger DM. Odds of Testing Positive for SARS-CoV-2 Following Receipt of 3 vs 2 Doses of the BNT162b2 mRNA Vaccine. JAMA Intern Med. 2022;182(2):179–184. pmid:34846533
  28. 28. Crowther MJ, Riley RD, Staessen JA, Wang J, Gueyffier F, Lambert PC. Individual patient data meta-analysis of survival data using Poisson regression models. BMC Med Res Methodol. 2012;12(1):34. pmid:22443286
  29. 29. Hardin J, Hilbe J. Generalized Estimating Equations. 1st ed. New York. 2022. doi: 10.1201/9781420035285.
  30. 30. Liang KY, Zeger SL. Longitudinal data analysis using generalized linear models. Biometrika. 1986;73(1):13–22.
  31. 31. Frazier PI, Cashore JM, Duan N, Henderson SG, Janmohamed A, Liu B, et al. Modeling for COVID-19 college reopening decisions: Cornell, a case study. Proc Natl Acad Sci U S A. 2022;119(2):e2112532119. pmid:34969678
  32. 32. Dunn OJ. Multiple Comparisons among Means. J Am Stat Assoc. 1961;56(293):52–64.
  33. 33. Danza P, Koo TH, Haddix M, Fisher R, Traub E, OY K, et al. SARS-CoV-2 Infection and Hospitalization Among Adults Aged >/ = 18 Years, by Vaccination Status, Before and During SARS-CoV-2 B.1.1.529 (Omicron) Variant Predominance—Los Angeles County, California, November 7, 2021-January 8, 2022. MMWR Morb Mortal Wkly Rep 2022;71 (5):177–181. pmid:35113851
  34. 34. Zhu Y, Liu S, Zhang D. Effectiveness of COVID-19 Vaccine Booster Shot Compared with Non-Booster: A Meta-Analysis. Vaccines (Basel). 2022;10(9):1396. pmid:36146474
  35. 35. Drawz PE, DeSilva M, Bodurtha P, Vazquez Benitez G, Murray A, Chamberlain AM, et al. Effectiveness of BNT162b2 and mRNA-1273 Second Doses and Boosters for Severe Acute Respiratory Syndrome Coronavirus 2 (SARS-CoV-2) Infection and SARS-CoV-2-Related Hospitalizations: A Statewide Report from the Minnesota Electronic Health Record Consortium. Clin Infect Dis. 2022;75(5):890–892. pmid:35137021
  36. 36. Sharma A, Oda G, Holodniy M. Effectiveness of mRNA-based vaccines during the emergence of SARS-CoV-2 Omicron variant. Clin Infect Dis. 2022:ciac325. pmid:35475889
  37. 37. Butt AA, Talisa VB, Shaikh OS, Omer SB, Mayr FB. Relative Vaccine Effectiveness of a SARS-CoV-2 mRNA Vaccine Booster Dose Against the Omicron Variant. Clin Infect Dis. 2022:ciac328. pmid:35511586
  38. 38. Suah JL, Tng BH, Tok PSK, Husin M, Thevananthan T, Peariasamy KM, et al. Real-world effectiveness of homologous and heterologous BNT162b2, CoronaVac, and AZD1222 booster vaccination against Delta and Omicron SARS-CoV-2 infection. Emerg Microbes Infect. 2022;11(1):1343–1345. pmid:35499301
  39. 39. Andrews N, Stowe J, Kirsebom F, Toffa S, Sachdeva R, Gower C, et al. Effectiveness of COVID-19 booster vaccines against COVID-19-related symptoms, hospitalization and death in England. Nat Med. 2022;28(4):831–837. pmid:35045566
  40. 40. Sritipsukho P, Khawcharoenporn T, Siribumrungwong B, Damronglerd P, Suwantarat N, Satdhabudha A, et al. Comparing real-life effectiveness of various COVID-19 vaccine regimens during the delta variant-dominant pandemic: a test-negative case-control study. Emerg Microbes Infect. 2022;11(1):585–592. pmid:35114893
  41. 41. Bar-On YM, Goldberg Y, Mandel M, Bodenheimer O, Freedman L, Kalkstein N, et al. Protection of BNT162b2 vaccine booster against Covid-19 in Israel. N Engl J Med. 2021;385(15):1393–1400. pmid:34525275
  42. 42. Regev-Yochay G, Gonen T, Gilboa M, Mandelboim M, Indenbaum V, Amit S, et al. Efficacy of a Fourth Dose of Covid-19 mRNA Vaccine against Omicron. N Engl J Med. 2022;386(14):1377–1380. pmid:35297591
  43. 43. Garcia-Beltran WF, et al. mRNA-based COVID-19 vaccine boosters induce neutralizing immunity against SARS-CoV-2 Omicron variant. Cell. 2022;185(3):457–466.e4. pmid:34995482
  44. 44. Goldberg Y, Mandel M, Bar-On YM, Bodenheimer O, Freedman L, Haas EJ, et al. Waning immunity after the BNT162b2 vaccine in Israel. N Engl J Med. 2021;385(24):e85. pmid:34706170
  45. 45. Sasser P, McGuine TA, Haraldsdottir K, Biese KM, Goodavish L, Stevens B, et al. Reported COVID-19 incidence in Wisconsin high school athletes in Fall 2020. J Athl Train. 2021;57(1):59–64. pmid:34129671
  46. 46. Vang KE, Krow-Lucal ER, James AE, Cima MJ, Kothari A, Zohoori N, et al. Participation in fraternity and sorority activities and the spread of COVID-19 among residential university communities—Arkansas, August 21-September 5, 2020. MMWR Morb Mortal Wkly Rep. 2021;70(1):20–23. pmid:33411698
  47. 47. Bar-On YM, Goldberg Y, Mandel M, Bodenheimer O, Amir O, Freedman L, et al. Protection by a Fourth Dose of BNT162b2 against Omicron in Israel. N Engl J Med. 2022;386(18):1712–1720. pmid:35381126
  48. 48. 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–1423. pmid:33626250
  49. 49. Jackson ML, Phillips CH, Benoit J, Kiniry E, Madziwa L, Nelson JC, et al. The impact of selection bias on vaccine effectiveness estimates from test-negative studies. Vaccine. 2018;36(5):751–757. pmid:29254838