COVID-19 vaccine impact on rates of SARS-CoV-2 cases and post vaccination strain sequences among healthcare workers at an urban academic medical center: a prospective cohort study =================================================================================================================================================================================== * Tara C. Bouton * Sara Lodi * Jacquelyn Turcinovic * Sarah E. Weber * Emily Quinn * Cathy Korn * Jacqueline Steiner * Elissa M. Schechter-Perkins * Elizabeth Duffy * Elizabeth J. Ragan * Bradford P. Taylor * Beau Schaeffer * Nancy Miller * Ravin Davidoff * William P. Hanage * John Connor * Cassandra Pierre * Karen R. Jacobson ## Abstract **Background** COVID-19 vaccine trials and post-implementation data suggest vaccination decreases SARS-CoV-2 infections. **Objective** Estimate COVID-19 vaccination’s impact on SARS-CoV-2 case rates and viral diversity among healthcare workers (HCW) during a high community prevalence period. **Design, Setting, Participants** A prospective cohort study from Boston Medical Center (BMC)’s HCW vaccination program, where staff received two doses of BNT162b2 or mRNA-1273. **Measurements** PCR-confirmed SARS-CoV-2 cases among HCWs from December 09, 2020 to February 23, 2021. Weekly SARS-CoV-2 rates per 100,000 person-day overall and by time from first injection (1-14 and >14 days) were compared with surrounding community rates. Viral genome sequences from SARS CoV-2 positive samples. **Results** SARS-CoV-2 cases occurred in 1.4% (96/7109) of HCWs given at least a first dose and 0.3% (17/5913) of HCWs given both vaccine doses. Adjusted SARS-CoV-2 infection rate ratios were 0.73 (95% CI 0.53-1.00) 1-14 days and 0.18 (0.10-0.32) >14 days from first dose. HCW SARS-CoV-2 cases >14 days from initial dose compared to within 14 days were more often older (46 versus 38 years, p=0.007), Latinx (10% versus 8%, p=0.03), and asymptomatic (48% versus 11%, p=0.0002). SARS-CoV-2 rates among HCWs fell below those of the surrounding community, with a 18% versus 11% weekly decrease respectively (p=0.14). Comparison of 48 SARS-CoV-2 genomes sequenced from post-first dose cases did not indicate selection pressure towards known spike-antibody escape mutations. **Limitations** Unable to adjust for infection risk outside of the workplace. Lack follow up on symptoms post SARS-CoV-2 diagnosis. Small number of vaccinated HCW cases. **Conclusion** Our results indicate a positive impact of COVID-19 vaccines on SARS-CoV-2 case rates. Post-vaccination isolates did not show unusual genetic diversity or selection for mutations of concern. ## Introduction On December 15, 2020, Boston Medical Center (BMC), an urban, safety-net hospital, started offering its 10,590 healthcare workers (HCWs) severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) vaccines BNT162b2 (Pfizer/BioNTech) and then mRNA-1273 (Moderna). The vaccination campaign coincided with Massachusetts’ second surge, which peaked at 6,000 new daily COVID-19 cases [1]. Vaccine was not available to anyone in the general population until February 1, 2021. This allowed examination of vaccine effectiveness during a period of higher prevalence and increased SARS-CoV-2 viral diversity than in the initial clinical trials [2, 3]. Concerns exist that newer SARS-CoV-2 variants have increased infectivity, modest decrease in neutralizing activity, and may impact vaccine effectiveness via escape from vaccine-induced immunity, specifically by mutations in the spike protein [4]. Assessment of viral genome sequencing of time-matched cases from vaccinated and unvaccinated individuals is needed to see whether such selection is evident. Our aims were threefold: 1) to compare infection rates among HCWs who did and did not receive a SARS-CoV-2 vaccine, 2) to compare infection rates over time between HCWs and the surrounding community in the months following the BMC vaccination initiative and, 3) to compare genomic and spike protein mutations between cases detected post-vaccination and among unvaccinated cases. ## Methods Boston University Medical Campus’ institutional review board approved the study. At BMC HCWs are screened daily for COVID-19 symptoms and tested if symptomatic. Asymptomatic testing is available to HCWs for workplace exposures, following out-of-state travel, and per request. All HCWs diagnosed with SARS-CoV-2 complete contact tracing and clinical questionnaires. As BMC was a BNT162b2 trial site, 66 HCWs had been vaccinated prior to the vaccine initiative and were included in analyses. We identified all HCWs with SARS-CoV-2 by reverse transcriptase-polymerase chain reaction (RT-PCR) between December 9, 2020 and February 23, 2021. We compared demographics and characteristics of cases detected post-vaccination by time since first dose to RT-PCR positive (1-14 days and >14 days from first dose to RT-PCR) using Pearson chi-squared, Fisher’s exact test, and Student t-tests. RT-PCR cycle threshold was dichotomized as above (negative) or below (positive) 24 across 5 instruments (Appendix Table 1) [5, 6, 7]. We computed the crude weekly case rates by vaccination status at the time of the tests (unvaccinated, 1-14 days, and >14 days from first dose to RT-PCR) as the number of weekly cases divided by person-days of follow-up during that week (Appendix text). To control for confounding from community trends in infection, we adjusted the rates using direct standardization to the weekly rates in Middlesex, Norfolk, and Suffolk counties in Massachusetts [1]. We compared the weekly decline in rates after December 30, 2020 (14 days after the vaccine initiative started) for BMC HCWs and the community using a negative binomial regression model including an interaction term between week and group (BMC HCW versus community) and an offset for person-day at risk. Residual isolates available from HCW SARS-CoV-2 cases tested at BMC were amplified using a modified ARTIC-primer based protocol and sequenced on an Illumina platform. Nucleotide substitutions, insertions, and deletions were identified with LoFreq [8] following alignment to the Wuhan Hu-1 reference sequence (NC_045512.2) [9] with Bowtie2 [10]. We then restricted to cases from between January 1 and February 23, 2021 to time-match vaccinated and unvaccinated cases, which controlled for the expected accumulation of mutations due to the virus’s high mutation rate. Significance calculations for synonymous and nonsynonymous mutations between vaccination groups used a Wilcoxon non-parametric test. ## Results During the study period, 67% (7,109/10,590) of eligible HCWs were vaccinated with at least one dose. Post vaccination SARS-CoV-2 cases occurred in 96/7,109 (1.3%) HCWs who received at least one dose, 17/5,913 (0.3%) HCWs given both doses, and 329/3,481 (9.5%) unvaccinated HCWs. Seventy percent (67/96) of post vaccination SARS-CoV-2 cases occurred within 14 days of the initial dose (Figure 1). HCWs diagnosed with SARS-CoV-2 >14 days after first dose were more often older (46 vs 38 years, p=0.007), Latinx (10% vs 8%, p=0.03), asymptomatic (48% vs 11%, p=0.0002), afebrile (100% vs 84%, p=0.02), and trended toward receipt of BNT162b2 (69% vs 48%, p=0.06), compared to HCWs diagnosed within 14 days of first dose (Table 1). There was no difference in RT-PCR cycle threshold. The adjusted rate ratios of SARS-CoV-2 infection were 0.73 (95% CI 0.53,1.00) and 0.18 (0.10,0.32) for 1-14 days and >14 days respectively from first dose compared to unvaccinated follow-up (Table 2). View this table: [Table 1:](http://medrxiv.org/content/early/2021/03/31/2021.03.30.21254655/T1) Table 1: Employee Characteristics Stratified by Vaccination-to-Diagnosis Timing View this table: [Table 2:](http://medrxiv.org/content/early/2021/03/31/2021.03.30.21254655/T2) Table 2: SARS-CoV-2 rate reductions among BMC HCW by vaccination status ![Figure 1:](http://medrxiv.org/https://www.medrxiv.org/content/medrxiv/early/2021/03/31/2021.03.30.21254655/F1.medium.gif) [Figure 1:](http://medrxiv.org/content/early/2021/03/31/2021.03.30.21254655/F1) Figure 1: Time elapsed from first dose of SARS-CoV-2 vaccination to positive RT-PCR* *Abbreviations: SARS-CoV-2, Severe acute respiratory syndrome coronavirus 2; RT-PCR, reverse transcriptase-polymerase chain reaction Following COVID-19 vaccine rollout to HCWs, case rates fell rapidly to below those of the community (Appendix Figure 1). Though not statistically significant, the observed weekly decrease in case rates was faster among BMC HCW than the community (HCW RR: 0.82; Community RR: 0.89; p=0.14). We sequenced isolates from 50% (48/96) of SARS-CoV-2 cases diagnosed after at least the first vaccine dose. Sequences from 36 post-first vaccine SARS-CoV-2 cases and 51 unvaccinated cases were included in the time-matched analysis. Whole genome analysis identified 97 total nonsynonymous single nucleotide variants (SNVs) in the spike protein, with 64 and 68 SNVs among unvaccinated and post at least initial vaccine isolates (Figure 2A). Of these SNVs, 35 were seen in both the unvaccinated and vaccinated cases, consistent with these infections occurring in an environment where similar variants were circulating. To determine whether there was greater accumulation of SNVs in SARS-CoV-2 genomes from either unvaccinated or vaccinated groups, we plotted the distribution of nonsynonymous SNVs between genomes in each population (Figure 2B). This analysis showed that most SNVs were found in one single genome from an individual case, with little accumulation of multiple mutations in either population. The exception to this trend were three mutations associated with the D614G variant which appeared in all genomes analyzed. ![Figure 2:](http://medrxiv.org/https://www.medrxiv.org/content/medrxiv/early/2021/03/31/2021.03.30.21254655/F2.medium.gif) [Figure 2:](http://medrxiv.org/content/early/2021/03/31/2021.03.30.21254655/F2) Figure 2: SARS-CoV-2 genome sequence variation among time-matched vaccinated and unvaccinated positive SARS-CoV-2 RT-PCR healthcare worker cases. A) Venn diagram illustrating the nonsynonymous SNVs in spike are distributed. A total of 97 SNVs were distributed over 87 time-matched sequenced genomes. 68 different SNVs were distributed over 36 genomes from individuals with a positive SARS-CoV-2 PCR following at least their first vaccination. 64 SNVs were distributed among 51 genomes from unvaccinated individuals. 35 SNVs were identified in genomes from both vaccinated and unvaccinated individuals. B) Distribution of spike SNVs between genomes in vaccinated and unvaccinated populations. C) Box-and-whiskers plots illustrating the number of nonsynonymous mutations found in infected HCW that were unvaccinated, infected less than 14 days after vaccination (≤ 14) or infected 15 or more days after vaccination as compared to the Wuhan-Hu-1 sequence. Significance calculations were calculated using a Wilcox non-parametric test. Analysis of the amino acid changes seen in vaccinated cases identified none with the antibody-evading 484K mutation. There were two amino acid changes (T1117A and N121D) seen in two separate infections (Table 3). Phylogenetic analysis showed that in both cases these were two separate infections with the same closely related lineage, and not multiple separate acquisitions of the same mutation, which would suggest convergent evolution. While genome sequences from cases among vaccines >14 days post first dose appeared more divergent from the original Wuhan-Hu-1 reference when compared with unvaccinated, this was no longer the case after controlling for the date of sampling; more recent samples are more divergent because more time has elapsed in which to accumulate mutations (Figure 2C). View this table: [Table 3.](http://medrxiv.org/content/early/2021/03/31/2021.03.30.21254655/T3) Table 3. Nonsynonymous mutations found in the Spike protein of vaccinated individuals who became infected (N=48). ## Conclusions Our findings show SARS-CoV-2 vaccines decreased case rates at a time of high community prevalence and domestic circulation of concerning variants [11]. While post vaccination SARS-CoV-2 cases do occur, the majority are in the two weeks following the first dose. We found adjusted rate reduction of new cases in vaccinated HCWs in days 1-14 and >14 days post first dose versus unvaccinated of 27% and 82% respectively, similar to findings in the original clinical trials and more recently reported in Israel and other parts of the United States [12, 13, 14, 15, 16]. Vaccine protection starts to have a greater impact in two weeks after the first vaccination dose. Individuals infected more than two weeks after first vaccination dose compared to before two weeks from first dose were older and reported fewer symptoms. While we observed a strong decrease in rates among vaccinated versus unvaccinated individual HCWs, we did not observe a significant decrease in rates among HCWs versus the community in the months following the HCW vaccination initiative. This suggests that HCW case rates fell steeply from their peak during the surge, but then flattened to falling community rates, even at a time when few members of the community were eligible for vaccine. HCWs at greatest risk for contracting SARS-CoV-2 infection outside of the workplace may have chosen not to receive vaccine, something we could not adjust for with our data, but could lead to the plateau in population-level effectiveness. We were unable to adjust for markers such as a social vulnerability index; however, disparities in vaccination uptake exist [17] and may lead to reduced vaccine effectiveness when compared with expectations from trial results. Even when access is not an obstacle, as in the case of HCWs, vaccine hesitancy and other barriers will need to be better understood and addressed to reach the levels of population immunity necessary to slow the spread of SARS-CoV-2 infection. In the meantime, mitigation measures such as masking, physical distancing, ventilation, and regular testing are still needed in hospital settings to protect both staff and patients. This is supported by recent modeling data suggesting that to avoid future outbreaks some mitigation measures will need to be maintained in the community at large even when high vaccine coverage is attained, as a result of the circulation of more transmissible variants [18]. Our data lacked follow-up from initial SARS-CoV-2 diagnosis, so initially asymptomatic cases may have progressed to symptomatic disease in all groups. Additionally, though HCWs diagnosed with SARS-CoV-2 more than 2 weeks from their initial dose were more likely to have received the BNT162b2 vaccine, this result likely reflects that BNT162b2 was given early when the force of infection was highest. There is intense interest in, and anxiety over, the potential of vaccination to select for mutants against which vaccines are less effective [19]. We found no early evidence of specific spike mutations or mutations associated with neutralizing vaccine escape, such as the E484K mutation. We observed overall genome divergence relative to the Wuhan-Hu-1 sequence as expected, with no evident increase of greater divergence in the Spike protein. Our study compared cases among vaccinees with those of contemporaneous unvaccinated individuals in the same population, which reflect the pool of virus to which vaccinated individuals were exposed, in order to avoid spurious associations between SNVs that happen to be locally common in different populations. It is a strength of this study that such time matched samples were available. We were also able to stratify participant samples by weeks from initial dose, not simply dichotomous vaccinated or unvaccinated, allowing observation of whether building immunity impacts the virus. Although we did not see clear trends, variation in the Spike sequence among vaccinated cases should be closely monitored by studies with larger sample size with careful consideration given to vaccine timing and comparison to concurrent unvaccinated controls. Widespread vaccination does appear to impact case rates in real world campaigns. We must continue to improve vaccine coverage for all globally, educate about the need for vigilance in social distancing in the weeks following first dose, and carefully monitor for any ongoing viral evolution that impacts upon vaccine efficacy. ## Data Availability A limited dataset is available from the authors by request. ## Acknowledgements The authors would like to thank John Goldie for helping with the BMC vaccination data, Boston Medical Center for providing financial support for this study, and all the patients and staff involved in the BMC COVID-19 efforts this year. The biologic samples used for the analyses presented here were obtained from the BMC/BU COVID-19 Biorepository. TCB is supported by the National Institutes of Health [grant number T32, DA013911] and the Burroughs Wellcome Fund/American Society for Tropical Medicine and Hygiene Postdoctoral Fellowship in Tropical Infectious Diseases. KJ was funded by the National Institutes of Health [1UL1TR001430]. ## Appendix ### SARS-CoV-2 Rate Reductions Among BMC HCW by Vaccination Status: Full Methodology Number of cases while unvaccinated, 1-14 days post vaccination and >14 days post vaccination were specified by week between December 9th 2020 and February 24th 2021. For each BMC employee, person-day of follow-up at risk while unvaccinated started on December 9th, 2020 and ended at the earliest of infection date, date of first vaccination shot or February 24th, 2021. Person-day follow-up at risk 1-14 days post vaccination started on the day after the first vaccination shot and ended at the earliest of infection date, the date corresponding to 14 after the first vaccination shot or end of the study. Person-day follow-up at risk >14 days post vaccination started 15 days after the first vaccination shot and ended at the earlier of infection date, or end of the study. Infection rates by vaccination status (unvaccinated, 1-14 post vaccination and >14 days post-vaccination) were expressed as number of cases per 100,000 person-day follow-up at risk overall and by week. Adjustment for temporal trends of infection in the community was performed by direct standardization to the weekly rates in 3 counties including and surrounding Boston where most BMC employees live. Directly standardized rates of the study population (unvaccinated BMC employees) are weighted averages of the weekly rates, where the weights are corresponding strata population size in the community. Weekly cases and estimated population in the 3 counties were obtained from the Massachusetts Department of Public Health publicly available dashboard and US Census Bureau, respectively and were used to compute the weekly infection rates in the community as number of cases per 100,000 person-days at risk. ![Appendix Figure 1:](http://medrxiv.org/https://www.medrxiv.org/content/medrxiv/early/2021/03/31/2021.03.30.21254655/F3.medium.gif) [Appendix Figure 1:](http://medrxiv.org/content/early/2021/03/31/2021.03.30.21254655/F3) Appendix Figure 1: Weekly Rates of SARS-CoV-2 Cases per 100,000 Individuals View this table: [Appendix Table 1:](http://medrxiv.org/content/early/2021/03/31/2021.03.30.21254655/T4) Appendix Table 1: SARS-CoV-2 Platforms Used To Diagnose via RT-PCR, Cycle Threshold Values* * Received March 30, 2021. * Revision received March 30, 2021. * Accepted March 31, 2021. * © 2021, Posted by Cold Spring Harbor Laboratory The copyright holder for this pre-print is the author. All rights reserved. The material may not be redistributed, re-used or adapted without the author's permission. ## References 1. 1.Massachusetts Department of Public Health. COVID-19 Response Reporting. Mass.gov. 2021. [https://www.mass.gov/info-details/covid-19-response-reporting%0D](https://www.mass.gov/info-details/covid-19-response-reporting). Accessed 17 Feb 2021. 2. 2.Wang Z, Schmidt F, Weisblum Y, Muecksch F, Barnes CO, Finkin S, et al. mRNA vaccine-elicited antibodies to SARS-CoV-2 and circulating variants. Nature. 2021. doi:10.1038/s41586-021-03324-6. [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1038/s41586-021-03324-6&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=33567448&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2021%2F03%2F31%2F2021.03.30.21254655.atom) 3. 3.Greaney AJ, Loes AN, Crawford KHD, Starr TN, Malone KD, Chu HY, et al. Comprehensive mapping of mutations in the SARS-CoV-2 receptor-binding domain that affect recognition by polyclonal human plasma antibodies. Cell Host Microbe. 2021;S1931-3128:00082–2. doi:10.1016/j.chom.2021.02.003. [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1016/j.chom.2021.02.003&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=33592168&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2021%2F03%2F31%2F2021.03.30.21254655.atom) 4. 4.Abdool Karim SS, de Oliveira T. New SARS-CoV-2 Variants - Clinical, Public Health, and Vaccine Implications. N Engl J Med. 2021 Mar 24. doi: 10.1056/NEJMc2100362. Epub ahead of print. [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1056/NEJMc2100362&link_type=DOI) 5. 5.Bullard J, Dust K, Funk D, Strong JE, Alexander D, Garnett L, et al. Predicting infectious SARS-CoV-2 from diagnostic samples. Clin Infect Dis. 2020 May 22:ciaa638. doi: 10.1093/cid/ciaa638. [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1093/cid/ciaa638&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=32442256&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2021%2F03%2F31%2F2021.03.30.21254655.atom) 6. 6.Basile K, McPhie K, Carter I, Alderson S, Rahman H, Donovan L, et al. Cell-based culture of SARS-CoV-2 informs infectivity and safe de-isolation assessments during COVID-19. Clin Infect Dis. 2020 Oct 24:ciaa1579. doi: 10.1093/cid/ciaa1579. [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1093/cid/ciaa1579&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=33098412&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2021%2F03%2F31%2F2021.03.30.21254655.atom) 7. 7.Kim MC, Cui C, Shin KR, Bae JY, Kweon OJ, Lee MK, et al. Duration of Culturable SARS-CoV-2 in Hospitalized Patients with Covid-19. N Engl J Med. 2021 Feb 18;384(7):671–673. doi: 10.1056/NEJMc2027040. [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1056/NEJMc2027040&link_type=DOI) 8. 8.Wilm A, Aw PPK, Bertrand D, Yeo GHT, Ong SH, Wong CH, et al. LoFreq: a sequence-quality aware, ultra-sensitive variant caller for uncovering cell-population heterogeneity from high-throughput sequencing datasets. Nucleic Acids Res. 2012;40:11189–201. doi:10.1093/nar/gks918. [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1093/nar/gks918&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=23066108&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2021%2F03%2F31%2F2021.03.30.21254655.atom) [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000313414800010&link_type=ISI) 9. 9.Severe acute respiratory syndrome coronavirus 2 isolate Wuhan-Hu-1, complete genome. 2020. [http://www.ncbi.nlm.nih.gov/nuccore/NC\_045512.2](http://www.ncbi.nlm.nih.gov/nuccore/NC_045512.2). Accessed 28 Oct 2020. 10. 10.Langmead B, Salzberg SL. Fast gapped-read alignment with Bowtie 2. Nat Methods. 2012;9:357–9. doi:10.1038/nmeth.1923. [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1038/nmeth.1923&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=22388286&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2021%2F03%2F31%2F2021.03.30.21254655.atom) [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000302218500017&link_type=ISI) 11. 11.United States Centers for Disease Control and Prevention. US COVID-19 Cases Caused by Variants. cdc.gov. 2021. [https://www.cdc.gov/coronavirus/2019-ncov/transmission/variant-cases.html](https://www.cdc.gov/coronavirus/2019-ncov/transmission/variant-cases.html). Accessed 9 Mar 2021. 12. 12.Polack FP, Thomas SJ, Kitchin N, Absalon J, Gurtman A, Lockhart S, et al. Safety and Efficacy of the BNT162b2 mRNA Covid-19 Vaccine. N Engl J Med. 2020;383:2603–15. doi:10.1056/NEJMoa2034577. [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1056/NEJMoa2034577&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=33301246&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2021%2F03%2F31%2F2021.03.30.21254655.atom) 13. 13.Baden LR, El Sahly HM, Essink B, Kotloff K, Frey S, Novak R, et al. Efficacy and Safety of the mRNA-1273 SARS-CoV-2 Vaccine. N Engl J Med. 2021;384:403–16. doi:10.1056/NEJMoa2035389. [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1056/NEJMoa2035389&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=33378609&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2021%2F03%2F31%2F2021.03.30.21254655.atom) 14. 14.Amit S, Regev-Yochay G, Afek A, Kreiss Y, Leshem E. Early rate reductions of SARS-CoV-2 infection and COVID-19 in BNT162b2 vaccine recipients. Lancet. 2021 Mar 6;397(10277):875–877. doi: 10.1016/S0140-6736(21)00448-7. [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1016/S0140-6736(21)00448-7&link_type=DOI) 15. 15.Keehner J, Horton LE, Pfeffer MA, Longhurst CA, Schooley RT, Currier JS, Abeles SR, Torriani FJ. SARS-CoV-2 Infection after Vaccination in Health Care Workers in California. N Engl J Med. 2021 Mar 23. doi: 10.1056/NEJMc2101927. Epub ahead of print. [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1056/NEJMc2101927&link_type=DOI) 16. 16.Daniel W, Nivet M, Warner J, Podolsky DK. Early Evidence of the Effect of SARS-CoV-2 Vaccine at One Medical Center. N Engl J Med. 2021 Mar 23. doi: 10.1056/NEJMc2102153. Epub ahead of print. [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1056/NEJMc2102153&link_type=DOI) 17. 17.Burger AE, Reither EN, Mamelund S-E, Lim S. Black-white disparities in 2009 H1N1 vaccination among adults in the United States: A cautionary tale for the COVID-19 pandemic. Vaccine. 2021;39:943– 51. doi:10.1016/j.vaccine.2020.12.069. [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1016/j.vaccine.2020.12.069&link_type=DOI) 18. 18.Moore S, Hill EM, Tildesley MJ, et al. Vaccination and non-pharmaceutical interventions for COVID-19: a mathematical modelling study. Lancet Infect Dis. 2021 Mar 18:S1473-3099(21)00143-2. doi: 10.1016/S1473-3099(21)00143-2. Epub ahead of print. [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1016/S1473-3099(21)00143-2&link_type=DOI) 19. 19.Saad-Roy CM, Morris SE, Metcalf CJE, et al. Epidemiological and evolutionary considerations of SARS-CoV-2 vaccine dosing regimes. Science. 2021 Mar 9:eabg8663. doi: 10.1126/science.abg8663. Epub ahead of print. [Abstract](http://medrxiv.org/lookup/ijlink/YTozOntzOjQ6InBhdGgiO3M6MTQ6Ii9sb29rdXAvaWpsaW5rIjtzOjU6InF1ZXJ5IjthOjQ6e3M6ODoibGlua1R5cGUiO3M6NDoiQUJTVCI7czoxMToiam91cm5hbENvZGUiO3M6Mzoic2NpIjtzOjU6InJlc2lkIjtzOjE3OiJzY2llbmNlLmFiZzg2NjN2MSI7czo0OiJhdG9tIjtzOjUwOiIvbWVkcnhpdi9lYXJseS8yMDIxLzAzLzMxLzIwMjEuMDMuMzAuMjEyNTQ2NTUuYXRvbSI7fXM6ODoiZnJhZ21lbnQiO3M6MDoiIjt9)