ABSTRACT
Recent data have shown elevated infection rates in several subpopulations at risk of SARS-CoV-2 infection and COVID-19, including immunocompromised (IC) individuals. Previous research suggests that IC persons have reduced risks of hospitalization and medically-attended COVID-19 with 2 doses of mRNA-1273 (SpikeVax; Moderna) compared to two doses of BNT162b2 (Comirnaty; Pfizer/BioNTech). The main objective of this retrospective cohort study was to compare real-world effectiveness of third doses of mRNA-1273 versus BNT162b2 at multiple time points on occurrence of COVID-19 hospitalization and medically-attended COVID-19 among IC adults in the US. The HealthVerity (HV) medical and pharmacy claims database, which contains data from >330 million patients, was the data source. Both subgroup and sensitivity analyses were conducted in addition to the core comparisons noted. In propensity score-adjusted analyses, receiving mRNA-1273 vs BNT162b2 as third dose was associated with 32% (relative risk [RR] 0.68; 95% confidence interval [CI] 0.51-0.89), 29% (0.71; 0.57-0.86), and 23% (0.77; 0.62-0.93) lower risk of COVID-19 hospitalization after 90, 180, and 270 days, respectively. Corresponding reductions in medically-attended COVID-19 were 8% (0.92; 0.86-0.98), 6% (0.94; 0.90-0.98), and 2% (0.98; 0.94-1.02), respectively. Our findings suggest a third dose of mRNA-1273 is more effective than a third dose of BNT162b2 in preventing COVID-19 hospitalization and breakthrough medically-attended COVID-19 among IC adults in the US.
INTRODUCTION
By September 2023, the United States (US) had registered >6 million cumulative hospitalizations due to cases of severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) infection and novel coronavirus disease 2019 (COVID-19), and >1 million associated deaths [1]. Although the availability of two mRNA vaccines, mRNA-1273 (SpikeVax; Moderna) and BNT162b2 (Comirnaty; Pfizer/BioNTech), has proved to be a critical tool against SARS-CoV-2 infection and COVID-19, there is a paucity of real-world data on the long-term effectiveness of vaccination. Recent data have shown increasing rates of infection in certain subpopulations at elevated risk of COVID-19 related morbidity and mortality, including immunocompromised individuals [2].
The safety and efficacy of two doses of both mRNA-1273 and BNT162b2 in reducing the risk of SARS-CoV-2 infection and of severe outcomes from COVID-19 in the general population and in adolescents has been demonstrated in phase 3 clinical trials [3–5] and in real-world studies [6, 7]. These results subsequently have been replicated in further observational studies, providing continued real-world evidence of the protective effects of two doses of both mRNA vaccines against SARS-CoV-2 infection, COVID-19 related hospitalization, and in-hospital mortality, despite the emergence of variants [6, 8, 9].
However, although high vaccine efficacy (VE) has been reported for both mRNA vaccines in the general population, these trials excluded immunocompromised (IC) individuals, such as those with underlying immunocompromising conditions and those prescribed immune-modifying therapies [10, 11]. Subsequent real-world data in the vaccinated IC population receiving two doses showed attenuated VE, with a higher risk of infection, hospitalization, death, persistent infection and shedding, viral evolution, reduced antibody and neutralization titers, and infection of household contacts [12]. In addition, observational data have suggested that there is a progressive dose- and time-dependent reduction in protection against SARS-CoV-2 infection among persons receiving a primary series of two doses [13–15].
In view of these increased risks, third and fourth doses of mRNA-based vaccines have been investigated in this vulnerable group [16] and, accordingly, both mRNA-1273 (50 or 100 μg; individuals aged ≥18 years) and BNT162b2 (30 μg; individuals aged ≥12 years) have been authorized and recommended for administration of a third primary dose ≥1 month after completion of the primary series in moderately to severely IC individuals [17, 18].
Our previous study showed that IC persons had reduced risks of hospitalization and medically-attended COVID-19 with two doses of mRNA-1273 compared to two doses of BNT162b2 [19], and the present study is an update assessing the comparative effectiveness of a third dose of both mRNA vaccines, regardless of number of doses included in the primary series. A matched cohort study showed that three doses of mRNA-1273 were associated with a significantly higher relative vaccine effectiveness (rVE) against SARS-CoV-2 infection and severe COVID-19 related outcomes compared to two doses, not only across subgroups of demographic and clinical characteristics but also in immunocompromised individuals [20, 21], and other studies have suggested that a third dose of BNT162b2 is safe and effective in preventing severe COVID-19-related outcomes [22, 23].
Previous influenza vaccine trials have demonstrated that high dose vaccines lead to improved immune responses in IC individuals compared with standard dose vaccines, which suggests that a high dose vaccine offers greater effectiveness for IC populations [24]. Although both mRNA-1273 and BNT162b2 employ the same mode of action, the composition of each vaccine is different, with the mRNA dosage and type of lipid nanoparticles used in the delivery system differing between vaccines. Also, the mRNA-1273 primary series contained 100 μg of mRNA vaccine and 50 μg for the booster [17], whereas BNT162b2 contained 30 μg of mRNA vaccine for each primary and booster dose [18]. Observational studies have shown differences between the two mRNA vaccines both in terms of immune response [25] and clinical effectiveness [19] in IC populations [24].
Study objective
The primary objective of this study was to compare real-world effectiveness of a third dose of mRNA-1273 versus a third dose of BNT162b2 at multiple time points (90, 180, and 270 days after start of follow-up) on occurrence of COVID-19 hospitalization and of medically-attended breakthrough COVID-19 among adults in the US who are immunocompromised and have received two doses of mRNA-based vaccine.
METHODS
Study design
This was a retrospective cohort study designed to compare real-world effectiveness of a third dose of mRNA-1273 versus a third dose of BNT162b2 on COVID-19 hospitalization and medically-attended breakthrough COVID-19 among US adults who are immunocompromised and have already received two doses of mRNA-based COVID-19 vaccines. The data source for this study was the HealthVerity (HV) aggregated and anonymized medical and pharmacy claims database, which holds data from >150 US healthcare payers for >330 million patients. Claims were drawn from the time period December 11, 2020 through August 31, 2022, a period which saw the emergence of dominant SARS-CoV-2 variants [26]. All available data of eligible patients were used in the derivation of baseline variables which were not otherwise defined or specified. At the time of data abstraction, the HV database contained approximately 3 million individuals with a third dose of mRNA-1273 and 4 million individuals with a third dose of BNT162b2, identified in HV closed-claim data feeds.
The sample of patients identified from closed claims (specifically private source 20) was enriched with an additional claims data feed (private source 17) linked through HealthVerity. This additional data feed includes adjudicated pharmacy claims sourced from a pharmacy benefit manager (PBM) and an associated pharmacy enrollment file. For a proportion (approximately 23%) of patients in private source 17, adjudicated medical claims are linked to the pharmacy enrollment file. For patients with linked medical claims in private source 17, 100% of adjudicated medical claims (i.e., all claims for those patients) are observable over the period in which a patient is enrolled in the PBM plan.
The index date (day 0) was defined as the day when study participants received a third dose of COVID-19 vaccine, either BNT162b2 or mRNA-1273. The baseline period was defined as the one-year period prior to index date. The start of follow-up was defined as 14 days after the index date. Patients were followed until the earliest of (a) receiving any additional COVID-19 vaccine, (b) observed occurrence of study outcome(s), (c) the date of a break in continuous enrollment in a study participant’s health plan, or (d) August 31, 2022 (Figure 1).
Study Population
The study population included IC adults in the US who had received three doses of mRNA-based COVID-19 vaccines (mRNA-1273 and BNT162b2) irrespective of the combination of the vaccines (e.g., three homologous mRNA-1273 or BNT162b2, or a heterologous combination of these), with the third dose considered the vaccine dose of interest, and the date of that dose establishing the index date. Immunocompromised individuals also were required to have been continuously enrolled in both medical and pharmacy insurance for at least one year before receiving the third dose and identified as immunocompromised on or before receiving the third dose.
Inclusion criteria
Individuals were eligible for inclusion in the analysis if they were aged ≥ 18 years on the index date, had continuous enrollment in medical or pharmacy plans during a 365-day baseline period prior to and including index date (i.e., no gaps were allowed), and had evidence of having received three doses of mRNA-based COVID-19, the first dose being on or after December 11, 2020. The required window between the first and second dose was determined by the combination of the first two doses: 28 ± 5-days for homologous 2 doses of mRNA-1273, 21 ± 5-days for homologous two doses of BNT162b2, and 28 ± 5-days for heterologous two doses. The third dose must have been received at least 42 days after the second dose. In addition, individuals had to be identified as immunocompromised based on evidence of one or more of the following criteria [19]: (a) blood or stem cell transplant in the two years prior to the index date; (b) history of organ transplant and taking immunosuppressive therapy within the 60 days prior to index date; (c) active cancer treatment in the 180 days prior to index date with an active cancer diagnosis in the 365 days prior to treatment; (d) prior history of a primary immunodeficiency disorder (e.g., for conditions such as DiGeorge syndrome and Wiskott-Aldrich syndrome); (e) history of an HIV diagnosis code prior to index date; or (f) immunosuppressive therapy in the 60 days prior to the index date.
Exclusion criteria
Individuals were excluded from the analysis if there was any evidence of the following: (a) receiving a non-mRNA COVID-19 vaccine up to and including the index date; (b) receipt of >1 mRNA vaccine dose on day 0; (c) receipt of any COVID-19 vaccine during day 1 through day 13, inclusive; (d) COVID-19 or related infectious disease during day 0 through day 13, inclusive; (e) missing or unknown sex or age on the index date; or (f) missing or unknown geographic region.
Exposure
The comparison cohort was defined as those individuals who received the mRNA-1273 vaccine as the third dose after completion of two primary doses of mRNA-based COVID-19 vaccines within the time window specified above. The third dose could be either 50 µg or 100 µg of vaccine, depending on whether it was administered as part of a three-dose primary series or as a two-dose primary series with 100 µg dosing followed by a 50 µg booster dose. The reference cohort was defined as individuals who received the BNT162b vaccine as their third dose after completing two primary mRNA-based COVID-19 vaccines. Vaccine administrations were identified using Current Procedural Terminology (CPT) and/or National Drug Code (NDC) codes, with time windows between doses based on health authority recommendation [27]. The list of relevant codes is provided in Supplemental Table 1. For the first two doses, only the primary series codes were used. As noted previously, the third mRNA based COVID-19 vaccine dose was allowed to be either part of a three-dose primary series (as recommended for the IC population) or a booster dose administered following a two-dose primary series.
Outcomes
Breakthrough COVID-19 hospitalization was identified based on presence of an International Classification of Diseases, 10 edition (ICD-10) diagnosis code for COVID-19 (U07.1) in any diagnosis field coincident with a hospitalization episode. HV does not include an inpatient indicator in medical claims for de-identification purposes. Inpatient claims were distinguished from outpatient claims by place of service, revenue code and bill type codes. The hospitalization episode was created by combining concurrent/adjacent inpatient claims with a 4-day allowable gap. To increase specificity, we further required evidence of respiratory distress during the same episode, such as bronchitis, cough, or use of supplemental oxygen [28,29]. The full list of diseases and codes are provided in Supplemental Table 2. Medically-attended breakthrough COVID-19 was identified using the same U07.1 ICD-10 diagnosis code in any place of service setting (i.e., inpatient, outpatient, emergency room, or urgent care).
Baseline Covariates
In addition to age, sex, geographic region, insurance type, and calendar month of index date, the following healthcare resource utilization measures and comorbid conditions were identified as potential confounders through literature review: (a) number of prior hospitalizations; (b) number of outpatient visits during baseline period; (c) number of unique immunosuppressive therapies used during baseline period; (d) the Charlson Quan Comorbidity Score during baseline period [30]; (e) the Kim claims-based frailty score during baseline period [31]; (f) recent medically-attended breakthrough COVID-19 infection (within six months prior to the index date); and (g) prior history of comorbid conditions. A full list of conditions is presented in Table 1. The algorithm used to derive these variables from claims data was documented in supplemental material published by Mues et al. [19].
Statistical analysis
Distributions of baseline variables within each vaccine cohort were described as numbers and percentages for categorical variables, and as mean (standard deviation), and median (interquartile range) for continuous variables. Absolute standardized differences (ASDs) were calculated to determine the balance of baseline variables (ASD ≤0.10) between the two cohorts prior and post weighting [32, 33]. Crude incidence rates for the two endpoints were calculated for each vaccine cohort and reported on a per thousand person-year (TPY) basis with associated 95% confidence intervals (CI). Crude Kaplan-Meier curves with 95% CIs were plotted.
Propensity score (PS) based methods were employed to account for measured baseline confounders [32]. Logistic regression was used to calculate the PS. Baseline covariates that were hypothesized to be confounding variables were included. Since participants had a negligible percentage of missing baseline demographic data (age, sex, and geographic region), persons with missing data for one or more of these data elements were excluded without imputation (801 out of 6.42 million; Figure 2). The inverse probability of treatment weighting (IPTW) method was applied in the main analysis. The weights were calculated as 1/PS for the mRNA-1273 group and 1/(1-PS) for the BNT162b2 group. A Cox proportional hazard model was initially used to estimate hazard ratios comparing two cohorts on both endpoints. However, the proportional hazard assumption was not met. A non-parametric method was substituted, and crude and weighted Kaplan-Meier curves were generated to estimate event-free survival probabilities at multiple time points (90, 180 and, 270 days). The comparison then was conducted on both relative and absolute scales: relative risk (RR) and risk difference (RD) to compare event-free survival probabilities between two cohorts at multiple time points. The 95% CI was obtained by bootstrapping with 1000 replicates.
* The entire HV dataset includes multiple sources of both closed and open claims during 2015 - 2022. Individuals from different sources and years can be counted duplicate and since many did not have enrollment file (such as open claims).
Subgroup analysis
Descriptive analyses were conducted for the following pre-specified subgroups of interest: participants aged 65+; participants without medically-attended COVID-19 in the 6 months prior to index day; participants with primary immunodeficiency; participants with active cancer; and participants diagnosed with HIV.
Sensitivity Analysis
A series of sensitivity analyses were performed: 1) a more sensitive definition of COVID-19 hospitalization was tested, specifically requiring only a ICD-10 diagnosis code (U07.1) in an inpatient setting; 2) a second alternative definition of COVID-19 hospitalization in which a COVID-19 hospitalization was captured using more restrictive definition of respiratory distress (cf. Supplemental Table 2); 3) using both open and closed claims from HV for outcome identification purpose; and 4) repeating the analyses using PS matching with a nearest neighbor approach and caliper = 0.01.
Software
Data cleaning and cohort selection was performed using the Aetion Evidence Platform: Substantiate software application for real-world data analysis, which has been scientifically validated for observational cohort studies using large healthcare databases [34]. All statistical analyses were conducted on SAS 9.4 (SAS Institute Inc.); plots were generated using R (4.3.1).
Ethical considerations
The study was designed to facilitate an understanding of real-world clinical practice, and no tests, treatments, or investigations were conducted. Since the study utilized anonymized retrospective data, no identifiable protected health information was abstracted.
RESULTS
A total of 113,027 immunocompromised adults were identified for inclusion in the primary analysis; of these, 52,943 and 60,084 received a third dose vaccination with mRNA-1273 and BNT162b2, respectively (Figure 2; Table 1). Individuals in the mRNA-1273 cohort were slightly older than in the BNT162b2 cohort (mean ages of 53.6 and 52.3 years, respectively) and more than one-half were aged 50-64 years. There were slightly fewer women in the mRNA-1273 cohort (54.4% vs. 53.5%, respectively; Table 1). Most patients were commercially insured (more than 60% among both cohorts). This immunocompromised population sample had a variety of comorbid conditions: for instance, in each cohort more than 60% had cardiovascular disease and more than one-half had hypertension. They also were actively seeking healthcare attention: on average, they had 38 outpatient claims during a year. In general, most demographics and comorbidities were comparable between these two cohorts except for calendar month of index date and geographic region. Since the BNT162b2 booster was approved earlier than mRNA-1273, the calendar month/year of index date tended to be earlier among individuals receiving BNT162b2 vs. mRNA-1273: 44.7% patients in the BNT162b2 cohort received their third dose during October-December 2021, while 45.8% patients in the mRNA-1273 cohort received their third dose during November-December 2021.
The median follow-up time of the two cohorts was comparable for both study endpoints (Table 2). The crude rate of COVID-19 hospitalization was lower among individuals who received a third dose of mRNA-1273 (7.89 per thousand person-years [TPY]; 95% CI, 6.88, 8.91) compared to those who received a third dose of BNT162b2 (9.99 per TPY; 95% CI, 8.93, 11.04). The crude rate of medically-attended breakthrough COVID-19 infection was also lower among the mRNA-1273 cohort (166.46 per TPY; 95% CI, 161.68, 171.25) compared to the BNT162b2 cohort (178.54 per TPY; 95% CI, 173.94, 183.15).
In the main comparison analyses, after IPTW, the measured baseline variables were considered balanced between the two cohorts (ASD ≤0.10) (Table 1). The weighted RR showed that, compared to BNT162b2, receiving mRNA-1273 as the third dose was associated with 32% (RR = 0.68, 95% CI: 0.51, 0.89), 29% (RR= 0.71, 95% CI: 0.57, 0.86), and 23% (RR= 0.77, 95% CI: 0.63, 0.93) lower risk of COVID-19 hospitalization after 90, 180, and 270 days, respectively (Table 3). The associations were statistically significant at all three time points. For medically-attended COVID-19, the magnitude of association was lower but consistent: receiving mRNA-1273 as the third dose was associated with 8% (RR= 0.92, 95% CI: 0.86, 0.98), 6% (RR= 0.94, 95% CI: 0.90, 0.98), and 2% (RR= 0.98, 95% CI: 0.94, 1.02) lower risk of medically-attended breakthrough COVID-19 after 90, 180, and 270 days. The association was statistically significant at 90 and 180 days, but not at 270 days. The absolute scale showed similar results: compared to BNT162b, receiving mRNA-1273 as the third dose was associated with 0.8, 1.5, and 1.7 fewer COVID-19 hospitalizations per thousand patients vaccinated after 90, 180, and 270 days, and 3.1, 5.0, and 3.0 fewer medically-attended breakthrough COVID-19 cases per thousand patients vaccinated after 90, 180, and 270 days (Table 3).
Results from sensitivity analyses provided further confidence to the main results. Findings were robust in response to differing definitions of COVID-19 hospitalization: adjusted RRs and RDs were comparable to the main analyses at all three different timepoints using either more sensitive or more specific definitions of COVID-19 hospitalization. Including open claims for outcome identification also provided similar results as the main analyses for both endpoints (Supplemental Table 3). We also tested the sensitivity of results to use of different analytical methods. All analyses, including the alternative endpoint definitions, were repeated using PS matching instead of IPTW. All measured baseline variables were balanced between cohorts after matching (Supplemental Table 4). The RRs and RDs were very similar between these methods on different endpoints at the three time points (Supplemental Tables 3 & 5).
Results also were similar among all five pre-specified subgroups: elderly adults (age 65+), three groups with important immunocompromising conditions (primary immunodeficiency, active cancer, and HIV), and a COVID-19 naïve group (i.e., persons with no evidence of medically-attended COVID-19 within 6 months prior to the index date), with the rate of COVID-19 hospitalization and medically-attended COVID-19 numerically lower in the mRNA-1273 cohort than in the BNT162b2 cohort (Supplemental Table 6). The statistical comparison was not conducted due to limited sample size and statistical power.
DISCUSSION
Observational studies have shown that a third dose of the mRNA-1273 or BNT162b2 vaccines is associated with a reduced risk of SARS-CoV-2 infection and associated severe COVID-19 related outcomes in immunocompetent adults [9, 20, 21]. In this large retrospective cohort study, we examined the real-world effectiveness of a third dose of mRNA-1273 versus a third dose of BNT162b2 in IC adults. It extended the work of a previous study that compared the effectiveness of two doses of mRNA-1273 and BNT162b2 vaccines among IC adults in the US [19]. This previous study showed that receiving two doses of mRNA-1273 was associated with lower medically-attended breakthrough COVID-19, and breakthrough COVID-19 hospitalization compared to receiving two doses of BNT162b2 among a population of IC adults.
Our findings are consistent with the previously established evidence of real-world effectiveness of both the mRNA-1273 and BNT162b2 vaccines among immunocompromised individuals compared to unvaccinated or partially vaccinated individuals [35]. The results suggested that rates of medically-attended COVID-19 and of COVID-19 hospitalization were lower among individuals who received a third dose of mRNA-1273 than those who received a third dose of BNT162b2. The evidence that a third dose of the mRNA-1273 vaccine is effective in preventing medically-attended COVID-19 outcomes and subsequent hospitalizations in IC individuals also supports recommendations from other studies [12, 15].
The results reported here have strength and uniqueness compared with prior studies that can provide valuable research insights in this fast-evolving COVID-19 field. To reflect the practice of mixed vaccination, this study did not require a homologous vaccine series as inclusion criteria. Furthermore, the third vaccine dose was inclusive of two common regimens: a three-dose primary series as recommended for immunocompromised patients, and a two-dose primary series plus booster. This generates a study sample reflective of a broader exposed population and may mimic more closely the real-world scenario that has evolved during the pandemic. Also, the methodology we used allows exploration of the effect of rVE on endpoints at different time points during follow-up. We investigated rVE at 90, 180, and 270 days after start of follow-up and results showed that mRNA1273 was more effective than BNT162b2 across different durations of follow-up.
From a statistical methodology standpoint, the results would confer causal interpretation if all three causal inference assumptions are met: consistency, exchangeability, and positivity [36]. We identified important potential confounders from published literature, derived them from the claims data based on established algorithms, and implemented two PS-based methods to account for measured ones. Even before adjustment, most measured baseline variables were comparable between cohorts and all of them were considered balanced after weighting/matching (Table 1). We empirically assessed loss to follow-up between two comparison cohorts to identify potential selection bias. The censoring reasons were distributed similarly between cohorts (Supplemental Table 7), and we found no evidence of differential loss to follow-up which potentially might lead to selection bias. Comparative effectiveness was estimated by a non-parametric method. Although less commonly used than Cox proportional hazard regression, it has fewer underlying assumptions and can estimate RDs and RRs which have more straightforward interpretation then the hazard ratio [37]. PS-based methods could be easily incorporated and has been applied in similar studies [22, 23, 38].
There are several limitations of this study. First, due to the lack of randomization inherent in claims data, there might be residual confounding due to unmeasured confounders, such as lifestyle or occupation among others. Second, all variables, including exposure, outcomes, and baseline covariates, were derived from medical and/or pharmacy claims. Random errors in these claims could lead to measurement error. However, it is believed to be non-differential between cohorts and would only bias results towards the null [39]. To mitigate this risk, this study employed a COVID-19 ICD-10 code to identify claims associated with COVID-19; this code has been shown to have high sensitivity, specificity, positive predictive value, and negative predictive value in various data sources [40, 41, 42]. We further conducted a series of sensitivity analyses using more sensitive or more specific outcome definitions which yielded consistent results. Lastly, there were a few limitations of the data sources we employed. To protect the privacy of the participants, during the anonymization process, people who were aged 85 years or older were coded as age 94 years. We believe the impact of this systematic misclassification was negligible since the prevalence of this age group in the overall sample was minimal. The HV database captured mainly commercially insured patients with some Medicaid and Medicare Advantage enrollees. Extrapolation of the conclusion to other populations, such as the uninsured, general Medicaid or Medicare fee-for-service enrollees, should be undertaken with caution.
Despite such limitations, real-world studies such as that reported herein play an important part in highlighting areas of concern not addressed in clinical trials and provide evidence of the comparative effectiveness of interventions in individuals seen in routine clinical practice.
CONCLUSION
The results from this observational comparative VE database study provide evidence that a third dose of mRNA-1273 is more effective than a third dose of BNT162b2 in preventing breakthrough medically-attended COVID-19 and subsequent hospitalizations among IC adults in the US.
Conflict of Interest
TS, LL, MG, JM, EB and NVV are employees of Moderna, Inc. and may hold stock/stock options in the company.
Other Disclosures
None
Funding
This study was funded by Moderna, Inc. Employees of Moderna participated in the design and conduct of the study; collection, management, analysis, and interpretation of the data; preparation, review, or approval of the manuscript; or the decision to submit the manuscript for publication.
Data sharing statement
Individual-level data reported in this study are not publicly shared. Upon request, and subject to review, Aetion, Inc. may provide the de-identified aggregate-level data that support the findings of this study. De-identified data (including participant data as applicable) may be shared upon approval of an analysis proposal and a signed data access agreement.
AUTHOR CONTRIBUTIONS
ACKNOWLEDGEMENTS AND DISCLOSURES
Medical writing and editorial assistance were provided by K Ian Johnson BSc, MBPS, SRPharmS of MEDiSTRAVA, in accordance with Good Publication Practice (GPP) guidelines, funded by Moderna, Inc., and under the direction of the authors.