Abstract
The results of a simulation-based evaluation of several policies for vaccine rollout are reported. In the presence of limited vaccine supply, this policy choice is a pressing issue for several countries worldwide, and the adopted course of action will affect the extension or easing of non-pharmaceutical interventions in the next months. We employ a suitably generalised, age-structure, stochastic SEIR (Susceptible → Exposed → Infectious → Removed) epidemic model that can accommodate quantitative descriptions of the major effects resulting from distinct vaccination strategies. The different rates of social contacts among distinct age-groups are informed by a recent survey conducted in Greece. The results are summarised and evaluated in terms of the total number of deaths and infections as well as life years lost. The optimal strategy is found to be one based on fully vaccinating the elderly/at risk as quickly as possible, while extending the time-interval between the two vaccine doses to 12 weeks for all individuals below 75 years old, in agreement with epidemic theory which suggests targeting a combination of susceptibility and infectivity. This policy, which is similar to the approaches adopted in the UK and in Canada, is found to be effective in reducing deaths and life years lost in the period while vaccination is still being carried out.
1 Introduction
Since December 2019, COVID-19 has presented a global threat to public health and to the worldwide economy, and it will likely continue to disrupt livelihoods until a high percentage of the population is vaccinated. High vaccination rates will be necessary to reach herd immunity in a short period of time. Standard theory (e.g. Andersson and Britton, 2000) suggests that a proportion approximately equal to 1-1/R0 of the population will have to become immune (either through vaccination or previous infection) in order to effectively suppress disease transmission, where R0 is the virus’ basic reproduction number. The actual vaccination coverage required is likely to vary due to population heterogeneity, previous levels of spread of infection, and other local factors. In addition, the exact value of R0 for SARS-CoV-2 under “ normal” conditions remains quite uncertain since there has been very little disease spread without some mitigation effort due to non-pharmaceutical interventions. Therefore, constrained scenarios are likely to give a realistic estimate of the effect of distinct vaccination policies and this approach is adopted in the present paper.
Assuming a vaccination coverage between 60%-80% of the population, 3.1-4.1 billion people worldwide will need to be vaccinated (Wang et al., 2020). With several seemingly highly efficacious vaccines available (efficacy estimated at 94.1%, 95% and 62% for Moderna, Pfizer-BioNTech and Oxford-AstraZeneca respectively) against COVID-19 disease (FDA, 2020; Polack et al., 2020; Voysey et al., 2021) it appears that a return to near-normality for society and for the economy may soon be possible. Unfortunately, limited supply is currently an impediment to achieving high vaccination coverage rapidly (Bollyky, 2021).
In addition to social distancing (Lewnard et al., 2020) and mass testing (Taipale et al., 2021), the fair allocation of scarce medical interventions such as vaccines presents ethical challenges as there are different allocation principles – treating people equally, favouring the worst-off, maximising total benefits, and promoting and rewarding social usefulness – and no single principle can address all morally relevant considerations (Emanuel et al., 2020; Persad et al., 2009). Modelling studies broadly agree that, when vaccine supply is limited, prioritising the elderly is a necessary strategy to reduce COVID-19 mortality, whereas the prioritisation of younger individuals would have an impact on reducing transmission (Bubar et al., 2021; Matrajt et al., 2020). This is in agreement with epidemic theory (e.g. Andersson and Britton, 2000) which suggests that the focus for disease control should be based on a combination of targeting susceptibility and infectivity. Therefore, assuming very scarce resources, it makes sense to focus upon the most vulnerable individuals in the population. On the other extreme is the presence of nearly unlimited vaccine supply, whence aiming for achieving herd immunity is straightforward. In this work we focus on the intermediate problem which many European countries are currently facing, and prioritisation of vaccines is of the essence.
Due to supply constraints, it was decided in the UK to delay the administration of the second dose based on the rationale that SARS-CoV-2 vaccines offer considerable protection after the first dose and that more people could benefit. Although this approach seems appealing, the impact of delaying the second dose is not straightforward as it depends on several parameters such as the efficacy of the first dose in time, the levels of transmission in the population, vaccination rollout, and the vaccine profile (reduction in symptoms or in symptoms and infection) (Matrajt et al., 2021; Paltiel et al., 2021a, 2021b). Country-specific information on the age-distribution of the population and social mixing patterns are also necessary to obtain realistic estimates.
The main contribution of this work is the evaluation of different vaccination strategies and their potential benefits in Greece, a typical country of the EU area in terms of vaccine availability and administration. The current strategy (strategy I) is to give the second vaccine dose three weeks after the first for the Pfizer vaccine, which currently consists of the largest portion of the available vaccines in the EU. We consider an alternative policy (strategy II) where, after the vaccination of medical personnel and those over 75, a portion of the available vaccines is distributed with a three-month time interval between the two doses. This aims for faster partial coverage of economically active individuals, therefore offering indirect protection to a larger proportion of the population and ultimately for potentially reducing the pandemic cost to public health and the society. This is implicitly informed by aiming for a combined effect of reducing susceptibility and infectivity in the population. Different scenarios of vaccine availability and transmission rates are considered, as well as different scenarios for the acquired immunity after the first dose for strategy II. We assess our results through simulation of an age-structured stochastic SEIR (Susceptible → Exposed → Infectious → Removed) epidemic model, suitably modified to account for the number of vaccinated individuals with different protocols.
2 The multitype S(V)EIR model and simulation description
The model used for the simulations for the different vaccination strategies is an age-structured stochastic SEIR model that accounts for different vaccinated populations, termed S(V)EIR henceforth. A schematic representation of the model is given in Figure 1. In order to evaluate the effects of different vaccination strategies to the Greek population, this model also accounts for the age composition of the population, the social mixing rates of different age groups, the intention to get vaccinated, as well as the different risk of death of each age group. The code for simulating the model is made freely available at: https://github.com/pbarmpounakis/Evaluating-the-effects-of-vaccine-rollout-policies-in-European-countries-A-simulation-study.
A detailed description of the model follows, while a table summarising the quantitative assumptions made is given in Appendix B.
States and vaccination effect assumptions
Two groups are considered for the vaccinated people representing the two distinct vaccination categories. In vaccination group 1, individuals receive the 2nd dose of the vaccine after 3 months while in vaccination group 2 it is given after 3 weeks. In both vaccination groups, individuals who received the 1st dose of the vaccine move to states V11 and V21 respectively and remain fully susceptible. Two weeks after the 1st dose individuals from both vaccination groups move to the second stage (V12, V22 respectively), whence immunity jumps to 68%. Individuals from vaccination group 2 remain at V22 for another 7 days when they take their 2nd dose and move to state 3 (V23) with their immunity jumping at 95%. Individuals from vaccination group 1 take their second dose 78 days after entering V12 and then move to state V15 with their immunity dropping linearly (from 68% to 50%). They move to V16 14 days later when their immunity jumps to 95%.
Transmission model assumptions
New infections from each state s and age group i follow a Binomial distribution with size given by the number of people in state s and age group i, and infection probability , where immunitys is the level of immunity at stage s, Ij is the number of infectious individuals at age group j, Nj is the total number of individuals at age group j, λij is the i, j entry in the transmission matrix λ, and n_ groups is the total number of different age groups. Following infection, individuals of age group i follow the Exposed (Ei) → Infectious (Ii) → Removed (Ri) path with a constant exposure time of 2 days and constant infectious period set at 6 days for non-vaccinated individuals and 2 or 3 days for vaccinated ones, depending upon the scenario of immunity waning. The total number of deaths is computed by multiplying the number in Ri with the infection fatality ratio (IFR) of each age group for the unvaccinated individuals and with [IFR x 5%] for those vaccinated (Haas et al., 2021).
Different scenarios for R0 and immunity
Transmission levels corresponding to R0=1.2 and R0=1.4 are considered along with various levels of immunity at each stage of vaccination for group 1. We ran 1000 simulations for each scenario and computed the median as well as 90% equal-tailed uncertainty intervals.
Worst Case Scenario: It is assumed that, during the time interval of 3 months between the first and second dose (strategy II), the acquired immunity drops linearly to 34% (Supplementary Figure 1, Appendix A). The effective infectious period of those vaccinated is reduced by 50% to 3 days.
Baseline Scenario: Here it is assumed that during the 3 months between the two doses (strategy II) the acquired immunity drops linearly to 50% and the infectious period of those vaccinated is set at 2 days (Figure 1).
Optimistic Scenario: In this case a constant immunity of 68% is assumed for the interval between the first and the second dose and the infectious period lasts 2 days (Suppl. Figure 2, Appendix A).
Fraction of vaccines given to general population
Different percentages are considered for the proportion of the available vaccines distributed under strategy II. These are set to 0% (strategy I), 20%, 50% and 100% (strategy II), and the resulting number of deaths and years lost are calculated for each of these cases.
Vaccine Availability
Two levels of vaccine availability are used, with a baseline level and a limited level with reduced number of vaccines; see Appendix B.
3 Results
Our main finding is that the optimal strategy in terms of the reduction in deaths and number of years lost, is the one that where all available vaccine doses are given under strategy II, using a time interval of three months between the two doses.
The results vary, depending upon the adopted scenarios regarding immunity waning and the value of R0, but they are robust in the sense that, in all cases, the optimal strategy, is found to be the one that allocates 100% of the available doses under strategy II.
Next, we present plots of the results for the baseline immunity waning scenario with R0=1.2, and with standard vaccine availability (Figures 3-5); detailed tables are given in Appendix C and several additional results from the optimistic and worst-case scenarios are summarised in a web supplement at: https://github.com/pbarmpounakis/Evaluating-the-effects-of-vaccine-rollout-policies-in-European-countries-A-simulation-study
4 Discussion
After vaccination of medical personnel, individuals with high risk of fatal outcome from COVID-19, and people aged over 75 years old with a time interval of 3-4 weeks between the two vaccine doses, the strategy of vaccinating the rest of the population with an interval of three months between the two doses can result in a significantly reduced number of deaths and years of life lost in the Greek population. When only 20% or 50% of the available vaccines are distributed using strategy II, the results are not significantly different to strategy I in terms of deaths, although they do provide an improvement in the number of life years saved. In conclusion, rolling out 100% of the available vaccines using the delayed second dose strategy appears to be the most effective option.
In the absence of detailed social contact data between different groups, we accounted for age groups as a surrogate for population composition, and we used the contact rate data between different age-groups from the recent survey (Sypsa et al., 2021). Therefore, the results reported here offer a conservative assessment since no attempt is made to prioritize individuals with many contacts such as mass transit employees, those working in the hospitality industry, super-markets and so on. Consequently, in practice, the benefits are expected to be even greater if a more targeted approach is adopted.
We used a multitype, age-structured, stochastic epidemic model with constant transmission rate and constant exposed and infectious periods. This approach has some limitations which are not expected to materially affect the results. First, in our model, we assumed that vaccine efficacy was mediated by a reduction in infections and not just in clinical disease. Recent modelling studies suggests that, if vaccines reduce symptomatic infection only, then the optimal protection for minimising deaths is prioritise older individuals (e.g., Matrajt et al., 2021). Our assumption seems realistic as recent data suggest that COVID-19 vaccines are effective in the prevention of infection (Amit et al., 2021; Thompson et al., 2021). Second, we assessed two scenarios for viral transmission rates (R0=1.2 and R0=1.4). For higher transmission levels, a recent study similarly found that vaccinating high-risk groups first constituted the optimal use of available vaccines (Matrajt et al., 2021). On the other hand, moderate transmission levels are a more realistic scenario as most counties continue to implement social distancing measures during vaccination. Alternative scenarios may be considered for the transmission rate, but the overall outcomes are not expected to be substantially influenced as the current assumptions regarding R0 may be thought of as an “ average” version of a time-varying rate. In addition, it is known (e.g. Andersson and Britton, 2000) that the final size of a stochastic epidemic is invariant to the presence of an exposed period and to different distributional assumptions on the infectious period duration. Hence, these assumptions will not alter the conclusions of this work. Other recent relevant results supporting our assumptions include (Tuite et al., 2021; Vasiliou et al., 2021).
Although we have chosen to primarily emphasize the results of the proposed approach in terms of quantities of interest in public health, additional gains are to be expected in terms of social and economic aspects of public life by offering faster vaccine coverage to the economically active population. An empirical application of the proposed approach is effectively followed in the United Kingdom, and the outcome seems to be a significantly faster overall reduction in SARS-CoV-2 circulation.
Overall, our results clearly indicate that, in the presence of a limited vaccine supply, distributing all available doses with a 3-month intermediate time interval could offer important advantages in terms of public health as well as to the wider society and the economy.
Data Availability
All data and code for this project is available at https://github.com/pbarmpounakis/Evaluating-the-effects-of-vaccine-rollout-policies-in-European-countries-A-simulation-study