Abstract
Vaccination against SARS-CoV-2 has been effective in reducing the burden of severe disease and death from COVID-19. Third doses of mRNA-based vaccines have provided a way to address waning immunity and broaden protection against emerging SARS-CoV-2 variants. However, their comparative effectiveness for a range of COVID-19 outcomes across diverse populations is unknown. We emulated a target trial using electronic health records of US veterans who received a third dose of either BNT162b2 or mRNA-1273 vaccines between 20 October 2021 and 8 February 2022, during a period that included Delta- and Omicron-variant waves. Eligible veterans had previously completed an mRNA vaccine primary series. We matched recipients of each vaccine in a 1:1 ratio according to recorded risk factors. Each vaccine group included 65,196 persons. The excess number of events over 16 weeks per 10,000 persons for BNT162b2 compared with mRNA-1273 was 45.4 (95% CI: 19.4, 84.7) for documented infection, 3.7 (2.2, 14.1) for symptomatic COVID-19, 10.6 (5.1, 19.7) for COVID-19 hospitalization, 2.0 (−3.1, 6.3) for COVID-19 intensive care unit admission and 0.2 (−2.2, 4.0) for COVID-19 death. After emulating a second target trial of veterans who received a third dose between 1 January and 1 March 2022, during a period restricted to Omicron-variant predominance, the excess number of events over 9 weeks per 10,000 persons for BNT162b2 compared with mRNA-1273 was 63.2 (95% CI: 15.2, 100.7) for documented infection. The 16-week risks of COVID-19 outcomes were low after a third dose of mRNA-1273 or BNT162b2, although risks were lower with mRNA-1273 than with BNT162b2, particularly for documented infection.
Similar content being viewed by others
Main
Additional doses of messenger RNA (mRNA)-based vaccines are being widely deployed to counter waning immunity and broaden protection against emerging variants of severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2)1,2. However, head-to-head comparisons of the effectiveness of a third dose of different mRNA-based vaccines have been lacking.
In a previous head-to-head comparison of breakthrough COVID-19 outcomes after the first dose of the BNT162b2 (Pfizer-BioNTech) and mRNA-1273 (Moderna) vaccines among 439,684 US veterans aged ≥18 yr, we found a low risk of documented infection and severe COVID-19 outcomes (for example, hospitalization and death) in a period marked by SARS-CoV-2 Alpha-variant predominance in both vaccine groups, but lower risk for the mRNA-1273 vaccine (with a similar pattern for documented infection in a period marked by Delta-variant predominance)3. It is unclear whether the same findings apply to third doses and a range of COVID-19 outcomes in periods marked by other SARS-CoV-2 variants, including Omicron.
A report from Spain estimated a third dose of mRNA-1273 to be 13% more effective than a third dose of BNT162b2 for the prevention of documented SARS-CoV-2 infection during a period of Omicron-variant predominance4. In other reports, effectiveness of mRNA-based vaccines against symptomatic COVID-19 was analysed separately for each third dose compared with no vaccination5 (or with the primary vaccination series6). These comparisons provide indirect evidence for comparative effectiveness but only if the comparators from each analysis had similar characteristics. Additional studies are needed to compare the effectiveness of a third dose of these vaccines head-to-head for severe COVID-19 outcomes to inform the choice of vaccine in coordinated public health responses. Ideally, any comparative effectiveness study should include racially diverse groups, evaluate potential differences in effectiveness according to time since the completion of the primary vaccination series and separately address time frames that include predominance of different SARS-CoV-2 variants.
Here we analysed data from the national healthcare databases of the Department of Veterans Affairs (VA), the largest integrated healthcare system in the US, to compare the effectiveness of a third dose of either the BNT162b2 or the mRNA-1273 vaccine among US veterans who had completed an mRNA vaccine primary series and received a third dose between (1) 20 October 2021 and 8 February 2022 (a period spanning Delta- and Omicron-variant predominance) or (2) 1 January and 1 March 2022 (a period restricted to Omicron-variant predominance). Recipients of each vaccine were matched in a 1:1 ratio according to their risk factors. Comparative effectiveness was estimated over 16 weeks in the combined Delta-Omicron period and over 9 weeks in the Omicron period for five COVID-19 outcomes: documented SARS-CoV-2 infection, symptomatic COVID-19 and COVID-19-related hospitalization, admission to an intensive care unit (ICU) and death.
Results
Primary analyses during a period spanning SARS-CoV-2 Delta- and Omicron-variant predominance
Study population and follow-up
Among 147,553 eligible recipients of a BNT162b2 third dose and 214,728 eligible recipients of an mRNA-1273 third dose, 65,196 BNT162b2 recipients were matched to 65,196 mRNA-1273 recipients (Fig. 1). The baseline characteristics of this matched population (Table 1) were similar to those of the eligible population (Supplementary Table 1). The median age was 70 yr (interquartile range, 62–74 yr), 96% of individuals were male and 24% were Black. The two vaccine groups had similar distributions of demographics, coexisting conditions and markers of healthcare utilization (see also Extended Data Fig. 1).
The median follow-up was 77 d (interquartile range, 61 to 97 d). Over a 16-week-follow-up, 2,994 SARS-CoV-2 infections were documented, of which 200 were detected as symptomatic COVID-19 within the VA healthcare system, 194 required hospitalization, 52 required ICU admission and 22 resulted in death.
Comparative effectiveness
Over a 16-week period spanning Delta- and Omicron-variant predominance, the estimated risk of documented infection was 353.9 (95% CI: 326.7, 373.2) events per 10,000 persons for the BNT162b2 third dose and 308.5 (95% CI: 276.8, 320.9) events per 10,000 persons for the mRNA-1273 third dose (Fig. 2). As expected, we found a nearly identical risk pattern in the two vaccine groups in the evaluations of two negative outcome controls7: symptomatic COVID-19 during the first 7 d after the third vaccine dose (Extended Data Fig. 2) and death from causes other than COVID-19 during follow-up (Extended Data Fig. 3).
The estimated 16-week risk ratios (95% CI) for recipients of a third dose of BNT162b2 as compared with mRNA-1273 were 1.15 (1.06, 1.30) for documented SARS-CoV-2 infection, 1.21 (1.12, 2.14) for symptomatic COVID-19, 1.64 (1.27, 2.79) for COVID-19 hospitalization, 1.37 (0.67, 3.14) for COVID-19 ICU admission and 1.08 (0.46, 6.39) for COVID-19 death (Table 2). The estimated risk differences (BNT162b2 minus mRNA-1273), expressed as events over 16 weeks per 10,000 persons (95% CI), were 45.4 (19.4, 84.7) for documented SARS-CoV-2 infection, 3.7 (2.2, 14.1) for symptomatic COVID-19, 10.6 (5.1, 19.7) for COVID-19 hospitalization, 2.0 (−3.1, 6.3) for COVID-19 ICU admission and 0.2 (−2.2, 4.0) for COVID-19 death.
Estimates were similar across subgroups defined at baseline according to age, race, months since completion of COVID-19 vaccine primary series (Table 3) and by vaccine type of primary series (risk ratio for documented infection comparing a third dose of BNT162b2 vs mRNA-1273, 1.19 (0.98, 2.65) among recipients of a BNT162b2 primary series, 1.04 (0.84, 2.32) among recipients of an mRNA-1273 primary series). Estimates were also similar in sensitivity analyses that excluded eligible individuals with previously documented SARS-CoV-2 infection (6.3%), those who did not receive a homologous COVID-19 vaccine primary series (3.3%) and those whose third dose could not be identified as a booster dose (17.2%) (data not shown).
Secondary analyses during a period of SARS-CoV-2 Omicron-variant predominance
Study population and follow-up
Among 25,557 eligible recipients of a BNT162b2 third dose and 36,809 eligible recipients of an mRNA-1273 third dose, 7,894 BNT162b2 recipients were matched to 7,894 mRNA-1273 recipients. Compared with the eligible population, the matched population was generally similar with respect to baseline demographic and clinical characteristics but included a higher percentage of men, White individuals, individuals with an urban residence and individuals who had received no SARS-CoV-2 tests in the past year (Supplementary Table 1). Baseline characteristics of the matched persons are shown in Supplementary Table 2. All variables were well-balanced between the two vaccine groups (Extended Data Fig. 1).
As compared with the matched population during a period spanning Delta- and Omicron-variant predominance, the matched population during a period of Omicron-variant predominance was, on average, younger and included a higher percentage of individuals who had received no SARS-CoV-2 tests in the previous year and no influenza vaccinations in the previous 5 yr at a VA facility (Supplementary Table 1).
The median follow-up was 54 d (interquartile range, 40–60 d). Over a 9-week follow-up, 214 SARS-CoV-2 infections were documented.
Comparative effectiveness
Over a 9-week period of Omicron-variant predominance, the estimated risk of documented SARS-CoV-2 infection was also higher with a third dose of the BNT162b2 vaccine than with a third dose of the mRNA-1273 vaccine; the estimated risk ratio was 1.57 (1.12, 2.10) and the estimated risk difference, expressed as events over 9 weeks per 10,000 persons, was 63.2 (15.2, 100.7) (see also Fig. 3).
Discussion
We quantified the comparative effectiveness of a third dose of the BNT162b2 and mRNA-1273 vaccines for the prevention of COVID-19 outcomes in the largest integrated healthcare system in the United States. Although the risks of all five measured COVID-19 outcomes (documented infection, symptomatic COVID-19 and COVID-19-related hospitalization, ICU admission and death) over 16 weeks were low for both vaccines during a period spanning Delta- and Omicron-variant predominance (risks <4% for documented infection and <0.03% for death, within each vaccine group), recipients of a third dose of BNT162b2 had an excess per 10,000 persons of 45 documented infections and 11 hospitalizations compared with recipients of a third dose of mRNA-1273. We also found a higher risk of documented infection among recipients of a third dose of BNT162b2 compared with mRNA-1273 over 9 weeks of follow-up during a period of Omicron-variant predominance, although this estimate was less precise due to a smaller number of eligible persons.
Few head-to-head comparisons of a third dose of mRNA vaccines are available. A report from Spain estimated mRNA-1273 boosters to be 13% more effective than BNT162b2 boosters for the prevention of documented SARS-CoV-2 infection during a period of Omicron-variant predominance, but more severe outcomes could not be evaluated4. In two other reports, effectiveness against symptomatic COVID-19 was evaluated separately for each mRNA booster either compared with no vaccination (in two test-negative case-control studies in England5) or compared with the primary series (in two matched cohort studies in Qatar6). In each report, an indirect comparison of these results requires that the controls from each study had similar characteristics, and more severe outcomes could not be evaluated for both mRNA booster groups.
The differential effectiveness we report in our study might be explained by the higher mRNA content of mRNA-1273 (50 μg for booster doses, 100 μg for third primary doses) compared with of BNT162b2 (30 μg for booster and third primary doses). Our primary analysis considered any third dose of these vaccines, of which 83% were distinguished as booster doses, and results were similar in sensitivity analyses that were restricted to these.
The strengths of our study are that first, the VA healthcare databases contain rich data on demographics and medical history, which allowed us to carefully match recipients of each vaccine type according to key confounders. Second, the databases contain detailed information on laboratory test results and healthcare encounters, which allowed us to capture outcomes related to COVID-19 in both outpatient and inpatient settings. Third, the large size of the study population allowed us to evaluate less common COVID-19 outcomes (hospitalizations, ICU admissions, deaths). Fourth, the demographic composition of the US veteran population allowed us to provide evidence for a diverse cohort (24% Black, 8% Hispanic) and to conduct subgroup analyses among older persons (≥70 yr of age) and Black persons.
Our study also has several potential limitations. First, as in any observational analysis of comparative effectiveness, the vaccine groups could differ with respect to risk factors for the outcomes. However, after rigorously matching recipients of each vaccine type, the two groups had similar demographics, comorbidities, time since completion of the primary vaccination series, number of SARS-CoV-2 tests in the previous year and markers of healthcare utilization (for example, number of primary care visits and of influenza vaccinations in the previous 5 yr). Further, much less confounding is expected when comparing recipients of different third doses of vaccines than when comparing recipients of a third dose with, for example, unvaccinated persons. In addition, our two analyses involving negative-outcome controls suggested little confounding.
Second, the possibility of outcome misclassification cannot be ruled out if veterans sought testing or care outside the VA healthcare system or if testing was differential between groups. However, our use of the VA COVID-19 National Surveillance Tool allowed us to integrate laboratory data with clinical notes to capture infections documented inside and outside the system, and the SARS-CoV-2 testing frequencies over follow-up were similar in both groups (Table 2). Further, we selected regular VA users who had a known residential address and who had routinely collected information (body-mass index, smoking status) to improve outcome ascertainment. We expect any residual incomplete ascertainment to be non-differential between the vaccination groups under comparison, although absolute risks may be slightly biased downwards. Finally, our study population was mostly made up of men (96%) and older persons (95% were >50 yr old), which may limit the generalizability of our findings.
Given the high effectiveness of a third dose of both the BNT162b2 and mRNA-1273 vaccines, either vaccine is strongly recommended to any individual. This study provides evidence of clear and comparable benefits of these vaccines for the most severe outcomes (the difference in estimated 16-week risk of death between the two groups was two-thousandths of one percent). While the differences in estimated risk for less severe outcomes between the two groups were small on the absolute scale, they may be meaningful when considering the population scale at which these vaccines are deployed. However, decisions about vaccination campaigns are complex and must incorporate considerations that extend beyond the scope of the present study. Therefore, our findings need to be combined with those from previous studies to inform future decisions about choice of vaccine in coordinated public health responses.
In summary, although the absolute risks of breakthrough SARS-CoV-2 infection and more severe COVID-19 outcomes were low regardless of the third mRNA-based COVID-19 dose received, this study involving a nationwide cohort of US veterans provides evidence of a lower 16-week risk of COVID-19-related outcomes among recipients of a third dose of the mRNA-1273 vaccine compared with recipients of a third dose of the BNT162b2 vaccine, particularly for documented infection. This pattern was consistent across a period spanning Delta- and Omicron-variant predominance and a period restricted to Omicron-variant predominance. Further evaluation of the comparative effectiveness and safety of additional doses of these vaccines is needed.
Methods
Specification of the target trials
We designed an observational analysis to emulate a hypothetical pragmatic trial of a third dose of BNT162b2 compared with mRNA-1273 and risks of COVID-19 outcomes in the VA healthcare system. Supplementary Table 3 summarizes the key protocol components.
Eligibility criteria included veteran status; age ≥65 yr or 18–64 yr with high risk of severe COVID-19 between 20 October and 18 November 2021 (based on the presence of at least one coexisting condition listed in Table 1)8 or ≥18 yr between 19 November 2021 and 8 February 2022 (based on national guidelines for third dose deployment)9,10; receipt of the second dose of an mRNA vaccine primary series at least 6 months earlier (based on the same guidelines); known residential address outside of a long-term care facility; and known smoking status and known body-mass index within the previous year. Individuals had to have used the VA healthcare system during the previous year (defined as receiving care at a station eligible to administer the vaccines under study and having at least one primary care visit) but not within the previous 3 d (which may indicate the start of symptomatic disease).
The interventions of interest were a third dose of either the BNT162b2 or the mRNA-1273 vaccine. To ensure balance of important characteristics across groups, eligible veterans in the target trial would be randomly assigned to one of these two vaccines within strata defined by calendar date of the third dose (5 d bins), calendar month of the second mRNA vaccine dose, age (5 yr bins), sex (male, female), race (White, Black, other, unknown), urbanicity of residence (urban, not urban), geographic location (coded as 19 categories of Veterans Integrated Services Network) and number of SARS-CoV-2 tests performed in the past 12 months (0, 1, ≥2).
The five outcomes of interest were documented SARS-CoV-2 infection, documented symptomatic COVID-19 and COVID-19-related hospital admission, ICU admission and death. For each eligible individual, follow-up started on the day the third dose of vaccine was received (baseline) and ended on the day of the outcome of interest, death, 112 d (16 weeks) after baseline, or the end of the study period (15 February 2022), whichever happened first.
Our target trial evaluates comparative effectiveness of a third dose of the vaccines during a period spanning times when SARS-CoV-2 Delta- and Omicron-variants were circulating. The Delta variant had decreased to a share of 26% of circulating variants in the United States as of 25 December 2021, as it was rapidly displaced by the Omicron variant, which rose to a 100% share as of 12 February 202211. To evaluate the comparative effectiveness during a period of Omicron-variant predominance, we considered a second target trial that was identical to the first trial except that the recruitment period was 1 January to 1 March 2022, and the only outcome of interest was documented SARS-CoV-2 infection because the period was too short to accrue a sufficient number of rarer, more severe outcomes.
Emulation of the target trials
We emulated both target trials using the VA healthcare databases3. Vaccination was identified using records from the Corporate Data Warehouse and the VA COVID-19 Shared Data Resource. SARS-CoV-2 infections were identified using the VA COVID-19 National Surveillance Tool3,12, which integrates data on laboratory tests with natural language processing of clinical notes to capture diagnoses documented inside and outside the VA healthcare system. Symptomatic COVID-19 was defined as ≥1 of the following symptoms documented within 4 d of SARS-CoV-2 infection: fever, chills, cough, shortness of breath or difficulty breathing, sore throat, loss of taste or smell, headache, myalgia/muscle pain, diarrhoea and vomiting. COVID-19 hospitalization was defined as a hospitalization within 21 d after documented SARS-CoV-2 infection, COVID-19 ICU admission was defined as an ICU admission during a COVID-19 hospitalization, and COVID-19 death was defined as a death within 30 d after documented SARS-CoV-2 infection. Supplementary Table 4 provides detailed information on all study variables and their ascertainment.
We attempted to mimic the stratified randomization of the target trial by matching persons who received a third dose of BNT162b2 and of mRNA-1273 on the basis of the calendar date of the third dose, calendar month of the second dose, age, sex, race, urbanicity of residence, geographic location, and number of SARS-CoV-2 tests performed in the past 12 months, using the same matching algorithm described in our previous study3.
To explore the possibility of residual confounding (for example, by underlying health status or healthcare-seeking behaviour), we incorporated two negative outcome controls7. First, we evaluated the risk of symptomatic COVID-19 in the first 7 d after the third vaccine dose. Second, we evaluated the risk of death from causes other than COVID-19 during the follow-up period.
Statistical analyses
Covariate balance after matching was evaluated by plotting the mean differences between variable values (standardized for continuous variables) for the vaccination groups, with a difference of 0.1 or less considered to be acceptable13.
Cumulative incidence (risk) curves for the vaccine groups were estimated using the Kaplan-Meier estimator14. We then calculated 16-week risks of each outcome and compared them between the vaccine groups via differences and ratios. We conducted subgroup analyses by age (<70 or ≥70 yr), race (Black or White), time since completion of the COVID-19 vaccine primary series (6–7, 8, or ≥9 months) and vaccine type of the primary series (BNT162b2 or mRNA-1273). We conducted sensitivity analyses that excluded eligible individuals who (1) had previously documented SARS-CoV-2 infection, (2) did not receive a homologous third dose compared with the COVID-19 vaccine primary series and (3) had a third dose that could not be identified as a booster dose on the basis of procedural codes available in the ‘Inpatient’ and ‘Outpatient’ domains, to compare recipients of vaccines at known doses, as the dose of mRNA-1273 differs for booster doses (50 μg) vs third primary doses (100 μg).
We used a nonparametric bootstrapping procedure (including both matching and subsequent analyses) with 500 iterations to calculate percentile-based 95% confidence intervals for all estimates.
Analyses were performed using R software version 3.6.0 (R Foundation for Statistical Computing) and SAS Enterprise Guide version 8.2 (SAS Institute).
Reporting summary
Further information on research design is available in the Nature Portfolio Reporting Summary linked to this article.
Data availability
The data that support the findings of this study are available from the VA. VA data are made freely available to researchers (behind the VA firewall) with an approved VA study protocol. More information is available at https://www.virec.research.va.gov or by contacting the VA Information Resource Center (VIReC) at VIReC@va.gov. Source data are provided with this paper.
Code availability
Access to the computer code used in this research is available by request to the corresponding author.
References
Garcia-Beltran, W. F. et al. mRNA-based COVID-19 vaccine boosters induce neutralizing immunity against SARS-CoV-2 Omicron variant. Cell 185, 457–466.e4 (2022).
Gruell, H. et al. mRNA booster immunization elicits potent neutralizing serum activity against the SARS-CoV-2 Omicron variant. Nat. Med. 28, 477–480 (2022).
Dickerman, B. A. et al. Comparative effectiveness of BNT162b2 and mRNA-1273 vaccines in U.S. veterans. N. Engl. J. Med. 386, 105–115 (2022).
Monge, S. et al. Effectiveness of mRNA vaccine boosters against infection with the SARS-CoV-2 omicron (B.1.1.529) variant in Spain: a nationwide cohort study. Lancet Infect. Dis. 22, 1313–1320 (2022).
Andrews, N. et al. Covid-19 vaccine effectiveness against the Omicron (B.1.1.529) variant. N. Engl. J. Med. 386, 1532–1546 (2022).
Abu-Raddad, L. J. et al. Effect of mRNA vaccine boosters against SARS-CoV-2 Omicron infection in Qatar. N. Engl. J. Med. 386, 1804–1816 (2022).
Lipsitch, M., Tchetgen Tchetgen, E. & Cohen, T. Negative controls: a tool for detecting confounding and bias in observational studies. Epidemiology 21, 383–388 (2010).
Underlying Medical Conditions Associated with Higher Risk for Severe COVID-19: Information for Healthcare Professionals (Centers for Disease Control and Prevention, accessed 3 March 2022).
FDA Authorizes Booster Dose of Pfizer-BioNTech COVID-19 Vaccine for Certain Populations [Press release] (US FDA, 22 September 2021).
Coronavirus (COVID-19) Update: FDA Takes Additional Actions on the Use of a Booster Dose for COVID-19 Vaccines [Press release] (US FDA, 20 October 2021).
COVID Data Tracker. Variant Proportions https://covid.cdc.gov/covid-data-tracker/#variant-proportions (Centers for Disease Control and Prevention, accessed 3 March 2022).
Chapman, A. et al. A natural language processing system for national COVID-19 surveillance in the US Department of Veterans Affairs. In Proc. 1st Workshop on NLP for COVID-19 at ACL 2020 (Association for Computational Linguistics, 2020; https://www.aclweb.org/anthology/2020.nlpcovid19-acl.10
Harder, V. S., Stuart, E. A. & Anthony, J. C. Propensity score techniques and the assessment of measured covariate balance to test causal associations in psychological research. Psychol. Methods 15, 234–249 (2010).
Kaplan, E. L. & Meier, P. Nonparametric estimation from incomplete observations. J. Am. Stat. Assoc. 53, 457–481 (1958).
Acknowledgements
This research was supported by the US Department of Veterans Affairs (VA) Office of Research and Development (ORD) Cooperative Studies Program (CSP) Epidemiology Center at the VA Boston Healthcare System through CSP #2032, by resources and the use of facilities at the VA Boston Healthcare System and VA Informatics and Computing Infrastructure (VINCI) (VA HSR RES 13-457) and by the use of data from the VA COVID-19 Shared Data Resource. B.A.D., A.L.M. and M.A.H. were supported by CSP #2032. B.A.D. was supported by a grant (R00 CA248335) from the National Institutes of Health. H.G. and B.R.F. were supported by a grant (MVP000) from the VA Million Veteran Program. M.J.F.M. was supported by a grant (T32 GM140972) from the National Institute of General Medical Sciences Interdisciplinary Training Program for Biostatisticians. We thank D. C. Posner, Y.-L. (A.) Ho and K. E. Lynch for insights on COVID-19 data extraction and phenotype definitions; C. A. Hoag for management of the administrative and regulatory aspects of the project; J. M. Gaziano for advice and support; the VA COVID-19 Shared Data Resource team for their contributions and support; and the VA healthcare providers, employees and volunteers for their dedication to caring for our veterans through this pandemic.
Author information
Authors and Affiliations
Contributions
B.A.D., H.G. and M.A.H. designed the study. H.G. and K.E.K. designed the data extraction and assembly workflow. H.G., B.R.F., M.J.F.M., N.A. and J.K.W. gathered the data. B.A.D., H.G. and A.L.M. designed the analytic pipeline. H.G. analysed the data. B.A.D., H.G. and M.A.H. vouch for the data and analysis. B.A.D. wrote the first draft of the manuscript. All authors critically reviewed the manuscript and decided to proceed with publication.
Corresponding author
Ethics declarations
Competing interests
B.R.F. reported completing a 6-month paid internship with Moderna before joining the Department of Veterans Affairs and owning five shares of Moderna, Inc. stock purchased in February 2020 outside the submitted work. M.A.H. reported receiving personal fees from Cytel and ProPublica outside the submitted work. The other authors declare no competing interests.
Ethical compliance
This work was approved by the VA Boston Healthcare System Research & Development Committee and Human Subjects Subcommittee, and an exemption of informed consent and HIPAA authorization was granted.
Peer review
Peer review information
Nature Microbiology thanks the anonymous reviewers for their contribution to the peer review of this work.
Additional information
Publisher’s note Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Disclaimer: The contents of this article do not represent the views of the US Department of Veterans Affairs or the US Government. The views and opinions expressed in this article are those of the authors and do not necessarily reflect the official policy or position of the Department of Health and Human Services and its agencies, including Biomedical Advanced Research and Development Authority and the Food and Drug Administration, as well as any other agency of the US Government. Assumptions made within and interpretations from the analysis do not necessarily reflect the position of any US Government entity.
Extended data
Extended Data Fig. 1 Covariate Balance (Love) Plot.
Shows the difference in means for conditions identified by the Centers for Disease Control and Prevention (CDC) as risk factors for severe COVID-19 in the matched population from the (a) primary analysis during a period spanning Delta- and Omicron-variant predominance and (b) secondary analysis during a period of Omicron-variant predominance. A strict balance cut-off was set at 0.1. An asterisk denotes the standardization of continuous variables.
Extended Data Fig. 2 Negative control 1: Cumulative Incidence of Symptomatic COVID-19 in the First 7 Days After the Third Vaccine Dose (October 20, 2021–February 15, 2022).
Solid blue line represents the risk curve for BNT162b2, dashed orange line represents the risk curve for mRNA-1273 and shaded areas represent pointwise 95% confidence intervals.
Extended Data Fig. 3 Negative control 2: Cumulative Incidence of Non-COVID-19 Death Over the Follow-up (October 20, 2021–February 15, 2022).
Solid blue line represents the risk curve for BNT162b2, dashed orange line represents the risk curve for mRNA-1273 and shaded areas represent pointwise 95% confidence intervals.
Supplementary information
Supplementary Information
Supplementary Tables 1–4.
Source data
Source Data Fig. 2
Statistical source data.
Source Data Fig. 3
Statistical source data.
Source Data Extended Data Fig. 1
Statistical source data.
Source Data Extended Data Fig. 2
Statistical source data.
Source Data Extended Data Fig. 3
Statistical source data.
Rights and permissions
About this article
Cite this article
Dickerman, B.A., Gerlovin, H., Madenci, A.L. et al. Comparative effectiveness of third doses of mRNA-based COVID-19 vaccines in US veterans. Nat Microbiol 8, 55–63 (2023). https://doi.org/10.1038/s41564-022-01272-z
Received:
Accepted:
Published:
Issue Date:
DOI: https://doi.org/10.1038/s41564-022-01272-z
This article is cited by
-
Comparative effectiveness of alternative intervals between first and second doses of the mRNA COVID-19 vaccines
Nature Communications (2024)