Evaluation of COVID-19 vaccination strategies with a delayed second dose

Model structure

We extended our previous model of agent-based COVID-19 transmission to include vaccination [13]. The model encapsulates the natural history of COVID-19 with classes of individuals including: susceptible; vaccinated; latently infected (not yet infectious); asymptomatic (and infectious); pre-symptomatic (and infectious); symptomatic with either mild or severe illness; recovered; and dead. Model population was stratified into six age groups of 0-4, 5-19, 20-49, 50-64, 65-79, and 80+ years based on US census data [14]. We sampled daily contacts within and between age groups from a negative-binomial distribution parameterized using an empirically-determined contact network (Appendix, Table A1) [15].

Disease dynamics

In our agent-based model, the risk of infection for people susceptible to COVID-19 depended on contact with infectious individuals that could be in asymptomatic, pre-symptomatic, or symptomatic stages of the disease. Using recent estimates, we parameterized the infectivity of asymptomatic, mild symptomatic, and severe symptomatic individuals to be 26%, 44%, and 89% relative to the pre-symptomatic stage [16–18]. For each infected individual, the incubation period was sampled from a Gamma distribution with a mean of 5.2 days [19]. A proportion of infected individuals developed symptomatic disease following a pre-symptomatic stage [20]. The duration of the pre-symptomatic stage and infectious period following symptom onset was sampled from a Gamma distribution with a mean of 2.3 days and 3.2 days, respectively [17,20,21]. Those who did not develop symptomatic disease remained asymptomatic until recovery, with an infectious period that was sampled from a Gamma distribution with a mean of 5 days [21,22]. We assumed that recovered individuals are immune against reinfection for the remainder of simulation timelines. Model parameters are summarized in Appendix, Table A2.

Infection outcomes

A proportion of severe symptomatic cases were hospitalized within 2-5 days of symptom onset [13,23], and thereafter did not contribute to the spread of infection. We assumed that all mild symptomatic cases and severely ill individuals self-isolated within 24 hours of symptom onset. The daily number of contacts during self-isolation was reduced by an average of 74%, based on a matrix derived from a representative sample population during COVID-19 lockdown [24]. Non-ICU and ICU admissions were parameterized based on age-stratified data for COVID-19 hospitalizations, and the presence of comorbidities [25,26]. The lengths of non-ICU and ICU stays were sampled from Gamma distributions with means of 12.4 and 14.4 days, respectively [27,28].

Vaccination

We implemented a two-dose vaccination campaign achieving 40% coverage of the entire population within one year. Prioritization was sequentially set to: (i) healthcare workers (5% of the total population) [29], adults with comorbidities, and those aged 65 and older; and (ii) other individuals aged 18-64 [30]. We assumed that 70% was the maximum achievable coverage in any age group, with an age-dependent distribution (Appendix, Table A3) mirroring seasonal influenza vaccination in the US [31]. We simulated a roll-out strategy offering 30 vaccine doses daily per 10,000 population, corresponding to 6.93 million vaccine doses per week for the entire US population. This rate corresponds to the goal of administering ∼100 million vaccine doses in the first 100 days, as outlined by the Biden administration [32].

For Pfizer-BioNTech vaccines administered with no delay, the interval between the first and second doses was 21 days [2]. This interval was 28 days for Moderna vaccines [3]. Vaccine efficacy against symptomatic and severe disease was set to be 52% and 80% for vaccines from Pfizer-BioNTech and Moderna, respectively, 14 days after the first dose [2,3]. The corresponding efficacies following the second dose were 95% and 94%, respectively, 1 week after the second dose [2,3]. In the base case, we assumed that vaccine protection against infection was the same as efficacy against symptomatic disease after each dose of the vaccine. In a scenario analysis, we reduced the protection against infection to 50% of the efficacy against disease.

In order to evaluate the impact of vaccination with DSD relative to the tested schedules in clinical trials, we assumed that the estimated efficacies of the first and second doses of both Pfizer-BioNTech and Moderna vaccines would be maintained if the maximum time interval between doses did not exceed 6 weeks [33]. Starting from week 7, we applied a waning rate of 5% per week (against both infection and symptomatic disease). As sensitivity analysis, we also simulated the model without reduction of vaccine efficacy up to 18 weeks after the first dose. Further, we investigated a scenario in which the same reduction factor of 5% per week was also applied to the efficacy of the second dose, if the second dose was received later than 6 weeks after the first dose (Appendix, Figure A3). In this scenario, we assumed that with a sufficiently long delay between doses, administration of the second dose only restored the vaccine efficacy equivalent to that conferred by one dose (Appendix, Figures A4-A5).

Model scenarios

We considered a range of 5%-30% pre-existing immunity in the population, with 10% for the base-case scenario [34,35]. To parameterize the model at a given level of pre-existing immunity, we ran simulations in the absence of vaccination, and derived the infection rates in different age groups once the overall attack rate reached that pre-specified level. The corresponding age distributions of recovered individuals were used for the initial population at the start of vaccination. We simulated the model with a 10,000 population for a time horizon of one year to evaluate the impact of DSD vaccination relative to the recommended schedule on reducing attack rate, hospitalizations, and deaths for both Pfizer-BioNTech and Moderna vaccines. For results presented here, outcomes were averaged over 1000 independent replications of each scenario.

DSD vaccination with waning efficacy of the first dose

We found that delaying the second dose of the Moderna vaccines would lead to a larger reduction of infections, hospitalizations, and deaths, as compared to the recommended 28-day interval between doses (Figure 1). Hospitalizations and deaths can be minimized by delaying the second dose by 9-12 weeks. For example, at 10% pre-existing immunity and equal vaccine efficacy against infection and symptomatic disease, we projected that delaying the second dose of Moderna vaccines by 9-12 weeks would avert an additional 0.75 − 0.87 hospitalizations and 0.08 − 0.10 deaths per 10,000 population compared with the recommended vaccination schedule. This DSD schedule remains the most effective strategy across a range of pre-existing levels from 5% to 30% and even when vaccine efficacy against infection is 50% of the efficacy against symptomatic disease (Appendix, Figure A1).

• Download figure
• Open in new tab

Figure 1.

Projected number of infections, hospitalizations, and deaths averted per 10,000 population for a DSD vaccination program compared to the recommended schedule of two-doses of Pfizer-BioNTech (with a 21-day interval) and Moderna (with a 28-day interval) vaccines. The efficacy of vaccines in blocking infection transmission was assumed to be the same as the efficacy against symptomatic and severe disease. Waning of vaccine efficacy was 5% per week, starting from week 7 after the first dose. Efficacy after the second dose was 94% for Moderna and 95% for Pfizer-BioNTech vaccines.

For the Pfizer-BioNTech vaccines, we project that the DSD strategy would lead to an increase in the number of infections compared to the 21-day recommended schedule between doses in most scenarios of delay interval, pre-existing immunity levels and relative vaccine efficacy against infection and disease (Figure 1, Appendix: Figure A1). However, depending on the level of pre-existing immunity, DSD may avert more hospitalizations and deaths. These additional benefits are due to prioritization of elderly and individuals with comorbidities receiving the first dose; thus, reducing severe outcomes among these high-risk individuals.

DSD vaccination without waning efficacy of the first dose

When the efficacy of the first dose did not wane for up to 18 weeks after being administered, we found that the DSD strategy averted more infections, hospitalizations, and deaths, compared to the recommended schedule for both Moderna and Pfizer-BioNTech vaccines (Figure 2). Similar to our results of waning vaccine efficacy, we found that delaying the second dose of Moderna vaccine by 9-12 weeks would lead to largest reductions in hospitalizations and deaths. At 10% pre-existing immunity, for example, such a strategy would avert 0.99 − 1.15 hospitalizations and 0.10 − 0.12 deaths per 10,000 population. The performance of DSD vaccination with Pfizer-BioNTech vaccines improved with an increasing delay of the second dose within 15 weeks from the recommended schedule. Although the scale of reduction varied, we observed a similar increasing trend for different levels of pre-existing immunity. Notably, DSD also outperformed the recommended vaccination schedule, even when there was a 50% reduction in the vaccine efficacy against infection compared to symptomatic disease (Appendix, Figure A2).

• Download figure
• Open in new tab

Figure 2.

Projected number of infections, hospitalizations, and deaths averted per 10,000 population for a DSD vaccination program compared to the recommended schedule of two-doses of Pfizer-BioNTech (with a 21-day interval) and Moderna (with a 28-day interval) vaccines. The efficacy of vaccines in blocking infection transmission was assumed to be the same as the efficacy against symptomatic and severe disease. There was no waning of vaccine efficacy within 15 weeks from the recommended schedule for the second dose. Efficacy after the second dose was 94% for Moderna and 95% for Pfizer-BioNTech vaccines.

DSD vaccination with waning efficacy of the first and second doses

Even with waning efficacy for both doses, DSD with Moderna vaccines would still lead to a larger reduction of infections, hospitalizations, and deaths, as compared to the recommended schedule, irrespective of whether vaccine efficacy against infection was the same or 50% lower than efficacy against disease (Appendix, Figures A4-A5). The lowest number of hospitalizations and deaths were achieved with a 9-12 week delay of the second dose. For Pfizer-BioNTech vaccines, we observed no clear advantage of DSD over the recommended schedule (Appendix, Figures A4-A5). Across the range of pre-existing immunity levels considered, DSD led to a higher number of infections. However, similar to the scenario of waning efficacy for only the first dose, DSD may avert more hospitalizations and deaths depending on the level of pre-existing immunity.