Evaluating the effects of second-dose vaccine-delay policies in European countries: A simulation study based on data from Greece ================================================================================================================================ * P. Barmpounakis * N. Demiris * I. Kontoyiannis * G. Pavlakis * V. Sypsa ## Abstract The results of a simulation-based evaluation of several policies for vaccine rollout are reported, particularly focusing on the effects of delaying the second dose of two-dose vaccines. In the presence of limited vaccine supply, the specific 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 (as well as some other model parameters) are informed by a recent survey conducted in Greece, but the conclusions are much more widely applicable. 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. ## 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 [1] suggests that a proportion approximately equal to 1-1/*R* of the population will have to become immune (either through vaccination or previous infection) in order to effectively suppress disease transmission, where *R* 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 *R* 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, and also due to the appearance of new variants. 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 [2]. 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 [3-5] 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 [6]. In addition to social distancing [7] and mass testing [8], 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 [9,10]. 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 [11,12]. This agrees with epidemic theory [1] 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 on 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 and Canada to delay the administration of the second dose of all vaccines, based on the rationale that SARS-CoV-2 vaccination offers considerable protection after the first dose and that more people could benefit [13,14]. 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) [15-17]. 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, primarily based on data from Greece, a typical country of the EU area in terms of vaccine availability and administration, with a population of around 10.8 million people. 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. The prioritization of the medical personnel and those over 75 years old is kept constant for every strategy. Our methodology examines scenarios where the two different vaccination schedules are combined in different proportions, allowing us to explore the optimal portion of the population that should be vaccinated using the extended three-month time interval between the two doses. This is something that has not been extensively explored in the literature, since earlier studies primarily focus on finding the optimal timing of the second dose, considering that the entire population will follow the same schedule [18-20]. Extending the time interval between doses to three months 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. We opted for a stochastic model because the most effective way to describe the spread of a disease is stochastic, based on the specification of the probability of disease transmission between two individuals. One may incorporate additional sources of stochasticity in the length of latent and infectious periods, but it is well known that such uncertainty is immaterial in terms of its effect on the outcome of an epidemic, particularly in a large population setting such as country-level studies; see for example [21]. Since the writing of this paper, related modelling techniques using deterministic dynamics have been proposed, that use optimization techniques to infer the model parameters and find the optimal dosing schedule [18,19]. An extension of a discrete-time, deterministic susceptible-infected-recovered model was used in [22] to plan the scheduling of first and second vaccine doses, with the underlying objective of the optimization problem being the concurrent minimization of both the healthcare impact of the epidemic and of the socio-economic impact due to the implementation of non-pharmaceutical interventions. A similarly extended SEIR model was also used to measure (through simulation) whether the effect of a standard vaccination schedule for different uptake scenarios is enough to stop the epidemic without the need of non-pharmaceutical interventions [20]. ## Materials and methods ### The multitype S(V)EIR model and simulation description The model used for the simulation of different vaccination strategies is an age-structured stochastic SEIR model that accounts for different vaccinated populations and vaccinated states, termed S(V)EIR henceforth. A schematic representation of the model is given in Fig 1. ![Fig 1.](http://medrxiv.org/https://www.medrxiv.org/content/medrxiv/early/2021/11/20/2021.05.19.21257486/F1.medium.gif) [Fig 1.](http://medrxiv.org/content/early/2021/11/20/2021.05.19.21257486/F1) Fig 1. Schematic representation of the S(V)EIR epidemic model for the baseline scenario of immunity waning. In order to evaluate the effects of different vaccination strategies, 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 publicly available at: [https://github.com/pbarmpounakis/Evaluating-the-effects-of-vaccine-rollout-policies-in-European-countries-A-simulation-study](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 summary of the quantitative assumptions made is given in S1 Table in Supporting information 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% [5]. 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% after two weeks. Individuals from vaccination group 1 take their second dose 78 days after entering V12 and then move to state V15 with their immunity changing based on different waning immunity scenarios described below. The stages V13, V14 and V15 account for the drop of the immunity due to waning vaccine efficacy 26, 52 and 78 days after the first dose, respectively. They move to V16 14 days later when their immunity jumps to 95%, pertaining to the reported efficacy of the mRNA-based vaccines which are mostly used in European countries [3-5] ### 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 ![Graphic][1], where *immunity**s* is the level of immunity at stage *s, I**j* is the number of infectious individuals at age group *j, N**j* 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 (E*i*) → Infectious (I*i*) → Removed (R*i*) path with a constant exposure time of 2 days, based on an average incubation time of approximately 5 days [23-25] and assuming that infectiousness starts approximately 2 days prior to the occurrence of symptoms [26-28]. The infectious period is also assumed constant and set at 6 days for non-vaccinated individuals [29-31] and 2 or 3 days for vaccinated ones, depending upon the scenario of immunity waning and vaccine efficacy in reducing the infectious period [32,33]. The choices of these values are conservative assuming that in reality people experiencing influenza-like symptoms will get tested and self-isolate, resulting in lower effective infectious period. The total number of deaths is computed by multiplying the number in R*i* with the infection fatality ratio (IFR) of each age group for the unvaccinated individuals as it was reported by [34] and with [IFR x 5%] for those vaccinated [35]. ### Different scenarios for *R**t* and immunity Transmission levels corresponding to a constant effective reproduction number *R**t*=1.2 and *R**t*=1.4 are considered along with various levels of immunity at each stage of vaccination for group 1. These choices of *R**t* suggest moderate transmission levels without the presence of a ‘hard lockdown’ for extended periods of time, closely resembling the policy that most European counties implement regarding social distancing measures during the vaccination period. *R**t* is calculated as the largest eigenvalue of the next generation matrix, using an appropriate contact matrix. The contact matrices used for the calculations of the values of *R**t* were based on a social contacts survey assessing contacts in Greece; we have used the data collected in the second half of September 2020. This contact matrix informs the relative infectivity between age-groups but, importantly, the scale is set by the value of *R**t*. We ran 1000 simulations for each scenario and computed the median as well as 90% equal-tailed uncertainty intervals. As precise data are not available for the precise course of infectivity and acquired immunity, three scenarios are considered. #### Worst Case Scenario It is assumed that, during the three months between the first and second dose (strategy II), the acquired immunity drops linearly to 34% (S1 Fig, Supporting information A). The effective infectious period of those vaccinated is reduced by 50% to 3 days. #### Baseline Scenario Here it is assumed that during the three 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 (Fig 1). #### Optimistic Scenario In this case a constant immunity of 68% is assumed for the entire time between the first and second dose, and the infectious period lasts 2 days (S2 Fig, Supporting information 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); the resulting number of deaths and life years lost are computed in each of these cases. Moreover, in the model we assume that individuals due for the second dose have priority over those waiting to have their first dose, keeping the time interval between the two doses intact. The remaining doses available each day are given to unvaccinated individuals. ### Vaccine availability Two levels of vaccine availability are considered; a baseline level and a limited level with reduced number of vaccines, see S1 Table in Supporting information B. ### Intention to get vaccinated The populations’ intention to get vaccinated is informed by a telephone survey contacted by [36]; see S1 Table in Supporting Information B. We assume that after the vaccination coverage of an age group reaches the percentage of people answering “*Probably/Definitely Yes”* to whether they intend to get vaccinated, the vaccination rollout continues to the next (younger) age group. The percentage of people answering “*Probably/Definitely No”* remains unvaccinated. The individuals answering “*Don’t know/Don’t answer”* are distributed among the two other groups according to their respected adjusted percentages. ### Ethics Statement This is a simulation study so no human participants were involved. For this reason, review and approval of this study from an ethical committee were not necessary. Since no human participants were recruited there are no available consent data or details. Regarding the data survey based on the [36,37], since this survey was conducted by one of the authors on the current manuscript, we also provide an ethics statement for this. The protocol of the survey was approved by the Institutional Review Board of the Hellenic Scientific Committee for the Study of AIDS and STDs. Participants provided oral informed consent. The survey was conducted over the phone and participants were asked for their consent at the beginning. ## 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 (Table 1). View this table: [Table 1.](http://medrxiv.org/content/early/2021/11/20/2021.05.19.21257486/T1) Table 1. Comparison of strategy I (0% of population with 3 month interval) and strategy II (100% of population with 3 month interval) Comparison of strategy I and strategy II (at 100% doses given) for 2021 with *R**t*=1.2, under the baseline immunity scenario and standard vaccine availability. “Gain” refers to the number of *fewer* deaths and life years lost under strategy II (extended interval between doses). The results vary between different immunity waning scenarios and different values of *R**t*, but they are robust in that the optimal strategy is always found to be the one that allocates 100% of the available doses under strategy II. The intermediate allocations of 20% or 50% of the available doses for strategy II show similar mortality for the whole population with Strategy I. But because fewer younger people get infected and die, there is a reduction in total life years lost (S2-S4 Tables). Next, we examine the resulting figures for daily deaths, life years lost, and daily infections, under the baseline immunity waning scenario, with *R**t*=1.2, and with standard vaccine availability (Figs 3-5); detailed summaries of simulation results are presented in S2-S9 Tables given in Supporting information C; several additional results based on both 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](https://github.com/pbarmpounakis/Evaluating-the-effects-of-vaccine-rollout-policies-in-European-countries-A-simulation-study) ![Fig 2.](http://medrxiv.org/https://www.medrxiv.org/content/medrxiv/early/2021/11/20/2021.05.19.21257486/F2.medium.gif) [Fig 2.](http://medrxiv.org/content/early/2021/11/20/2021.05.19.21257486/F2) Fig 2. Cumulative number of deaths. Cumulative number of deaths over time when different percentages of doses are allocated under strategy II, with *Rt*=1.2, immunity drop between the two vaccine doses is at the baseline scenario, and with standard vaccine availability. ![Fig 3.](http://medrxiv.org/https://www.medrxiv.org/content/medrxiv/early/2021/11/20/2021.05.19.21257486/F3.medium.gif) [Fig 3.](http://medrxiv.org/content/early/2021/11/20/2021.05.19.21257486/F3) Fig 3. Number of new daily deaths. Number of new daily deaths, when different percentages of doses are allocated under strategy II, *Rt*=1.2, immunity drop is at the baseline scenario, and with standard vaccine availability. [A zoomed-in version of this figure is available in the supplement, offering higher resolution locally.] ![Fig 4.](http://medrxiv.org/https://www.medrxiv.org/content/medrxiv/early/2021/11/20/2021.05.19.21257486/F4.medium.gif) [Fig 4.](http://medrxiv.org/content/early/2021/11/20/2021.05.19.21257486/F4) Fig 4. Total number of years of life lost Total number of years of life lost when different percentages of doses are allocated under strategy II, *Rt*=1.2, immunity drop is at the baseline scenario, and with standard vaccine availability. ![Fig 5.](http://medrxiv.org/https://www.medrxiv.org/content/medrxiv/early/2021/11/20/2021.05.19.21257486/F5.medium.gif) [Fig 5.](http://medrxiv.org/content/early/2021/11/20/2021.05.19.21257486/F5) Fig 5. New daily infections. New daily infections when different percentages of doses are allocated under strategy II, *Rt*=1.2, immunity drop is at the baseline scenario, and with standard vaccine availability. ## Discussion After the vaccination of medical personnel, high-risk individuals, and people aged over 75 years old with a time interval of 3-4 weeks between doses, the strategy of vaccinating the rest of the population with an interval of three months between the two doses (strategy II) can result in a significantly reduced number of deaths and years of life lost. 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 [37]. 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. Although this approach of course has some limitations, they 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 suggest that, if vaccines reduce symptomatic infection only, then the optimal protection for minimising deaths is obtained by prioritising older individuals [15]. This assumption is realistic especially in view of recent data suggesting that COVID-19 vaccines are indeed effective in the prevention of infection at least before the occurrence of the Delta strain [38,39]. Second, we assessed two scenarios for viral transmission rates (*R**t*=1.2 and *R**t*=1.4). For higher transmission levels, a recent study similarly found that vaccinating high-risk groups first constituted the optimal use of available vaccines [15]. On the other hand, moderate transmission levels are a more realistic scenario as most counties continue to implement moderate 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 *R**t* may be thought of as an “average” version of a time-varying rate. In addition, it is known [1] 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 [40,41]. The main conclusions of the present study and all relevant assumptions made about vaccine efficacy are in broad agreement with the results obtained using optimization techniques for a model calibrated using data from Italy [18], despite the fact that the authors only allowed for the vaccine to protect against transmission and not disease. Recent results from Israel suggest that there is also protection against hospitalization and death [35] and therefore these results may be conservative. Similar conclusions are reported when varying vaccine availability and using alternative efficacy assumptions [19]. Assuming that efficacy remains constant after the first dose, [42] used simulation and showed that the effectiveness of vaccination programmes in reducing infections, hospitalizations and deaths is maximized with a delay of 12 to 15 weeks for both the Pfizer and Moderna vaccines. Similar recommendations on delaying the second dose for individuals below the age of 65 are made in [43]. 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 Canada, and the outcome seems to be a significantly faster 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](https://github.com/pbarmpounakis/Evaluating-the-effects-of-vaccine-rollout-policies-in-European-countries-A-simulation-study) [https://github.com/pbarmpounakis/Evaluating-the-effects-of-vaccine-rollout-policies-in-European-countries-A-simulation-study](https://github.com/pbarmpounakis/Evaluating-the-effects-of-vaccine-rollout-policies-in-European-countries-A-simulation-study) ## Supporting information ### Supporting information A Graphical presentation of different scenarios ![S1 Fig.](http://medrxiv.org/https://www.medrxiv.org/content/medrxiv/early/2021/11/20/2021.05.19.21257486/F6.medium.gif) [S1 Fig.](http://medrxiv.org/content/early/2021/11/20/2021.05.19.21257486/F6) S1 Fig. Schematic representation of the S(V)EIR epidemic model for the *worst-case* scenario regarding immunity waning. ![S2 Fig.](http://medrxiv.org/https://www.medrxiv.org/content/medrxiv/early/2021/11/20/2021.05.19.21257486/F7.medium.gif) [S2 Fig.](http://medrxiv.org/content/early/2021/11/20/2021.05.19.21257486/F7) S2 Fig. Schematic representation of the S(V)EIR epidemic model for the *optimistic* scenario regarding immunity waning. ### Supporting information B View this table: [S1 Table.](http://medrxiv.org/content/early/2021/11/20/2021.05.19.21257486/T2) S1 Table. Model Assumptions ### Supporting information C Baseline Scenario - Vaccine Availability - *R**t*=1.2 View this table: [S2 table.](http://medrxiv.org/content/early/2021/11/20/2021.05.19.21257486/T3) S2 table. Cumulative number of deaths, when 0% of vaccines allocated to ages 18-74 View this table: [S3 Table.](http://medrxiv.org/content/early/2021/11/20/2021.05.19.21257486/T4) S3 Table. Cumulative number of deaths, when 20% of vaccines allocated to ages 18-74 View this table: [S4 Table.](http://medrxiv.org/content/early/2021/11/20/2021.05.19.21257486/T5) S4 Table. Cumulative number of deaths, when 50% of vaccines allocated to ages 18-74 View this table: [S5 Table.](http://medrxiv.org/content/early/2021/11/20/2021.05.19.21257486/T6) S5 Table. Cumulative number of deaths, when 100% of vaccines allocated to ages 18-74 View this table: [S6 Table.](http://medrxiv.org/content/early/2021/11/20/2021.05.19.21257486/T7) S6 Table. Cumulative number of infections, when 0% of vaccines allocated to ages 18-74 View this table: [S7 Table.](http://medrxiv.org/content/early/2021/11/20/2021.05.19.21257486/T8) S7 Table. Cumulative number of infections, when 20% of vaccines allocated to ages 18-74 View this table: [S8 Table.](http://medrxiv.org/content/early/2021/11/20/2021.05.19.21257486/T9) S8 Table. Cumulative number of infections, when 50% of vaccines allocated to ages 18-74 View this table: [S9 Table.](http://medrxiv.org/content/early/2021/11/20/2021.05.19.21257486/T10) S9 Table. Cumulative number of infections, when 100% of vaccines allocated to ages 18-74 * Received May 19, 2021. * Revision received November 20, 2021. * Accepted November 20, 2021. * © 2021, Posted by Cold Spring Harbor Laboratory This pre-print is available under a Creative Commons License (Attribution-NoDerivs 4.0 International), CC BY-ND 4.0, as described at [http://creativecommons.org/licenses/by-nd/4.0/](http://creativecommons.org/licenses/by-nd/4.0/) ## References 1. 1.Andersson, H., & Britton, T. (2000). Stochastic Epidemic Models and Their Statistical Analysis. Lecture Notes in Statistics, vol 151, Springer, New York, NY. 2. 2.Wang, W., et al. (2020). Global, regional, and national estimates of target population sizes for covid-19 vaccination: descriptive study. BMJ 2020;371:m4704. 3. 3.FDA (2020). FDA Briefing Document, Moderna COVID-19 Vaccine, [https://www.fda.gov/media/144434/download](https://www.fda.gov/media/144434/download). Last accessed December 27, 2020. 4. 4.Polack, F.P., et al. (2020). Safety and Efficacy of the BNT162b2 mRNA Covid-19 Vaccine. New England Journal of Medicine, 383(27), 2603–2615. [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1056/NEJMoa2034577&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=http://www.n&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2021%2F11%2F20%2F2021.05.19.21257486.atom) 5. 5.Voysey, M. et al. (2021). Safety and efficacy of the ChAdOx1 nCoV-19 vaccine (AZD1222) against SARS-CoV-2: an interim analysis of four randomised controlled trials in Brazil, South Africa, and the UK, The Lancet. 397:10269, 99–111. [PubMed](http://medrxiv.org/lookup/external-ref?access_num=http://www.n&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2021%2F11%2F20%2F2021.05.19.21257486.atom) 6. 6.Bollyky, T.J. (2021). U.S. COVID-19 Vaccination Challenges Go Beyond Supply. Annals of Internal Medicine. 174:558–559. 7. 7.Lewnard, J.A. et al., (2020). Scientific and ethical basis for social-distancing interventions against COVID-19. The Lancet Infectious Diseases, volume 20, issue 6, 631–633. [PubMed](http://medrxiv.org/lookup/external-ref?access_num=http://www.n&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2021%2F11%2F20%2F2021.05.19.21257486.atom) 8. 8.Taipale, J., Kontoyiannis, I., Linnarsson, S. (2021). Population-scale testing can suppress the spread of infectious disease. arxiv:2104.06857 [q-bio.PE]. 9. 9.Emanuel, E.J., Persad, G., Upshur, R., Thome, B., Parker, M., Glickman, A., et al., (2020). Fair Allocation of Scarce Medical Resources in the Time of Covid-19. New England Journal of Medicine 382(21), 2049–2055. [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1056/NEJMsb2005114&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=32202722&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2021%2F11%2F20%2F2021.05.19.21257486.atom) 10. 10.Persad, G., Wertheimer, A., Emanuel, E.J., (2009). Principles for allocation of scarce medical interventions. The Lancet, 373(9661), 423–431. [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1016/S0140-6736(09)60137-9&link_type=DOI) 11. 11.Bubar KM, Reinholt K, Kissler SM, Lipsitch M, Cobey S, Grad YH, Larremore DB. Model-informed COVID-19 vaccine prioritization strategies by age and serostatus. Science. 2021 Feb 26;371(6532):916–921. doi: 10.1126/science.abe6959. Epub 2021 Jan 21. PMID: 33479118; PMCID: PMC7963218. [Abstract/FREE Full Text](http://medrxiv.org/lookup/ijlink/YTozOntzOjQ6InBhdGgiO3M6MTQ6Ii9sb29rdXAvaWpsaW5rIjtzOjU6InF1ZXJ5IjthOjQ6e3M6ODoibGlua1R5cGUiO3M6NDoiQUJTVCI7czoxMToiam91cm5hbENvZGUiO3M6Mzoic2NpIjtzOjU6InJlc2lkIjtzOjEyOiIzNzEvNjUzMi85MTYiO3M6NDoiYXRvbSI7czo1MDoiL21lZHJ4aXYvZWFybHkvMjAyMS8xMS8yMC8yMDIxLjA1LjE5LjIxMjU3NDg2LmF0b20iO31zOjg6ImZyYWdtZW50IjtzOjA6IiI7fQ==) 12. 12.Matrajt, L., Eaton, J., Leung, T., Brown, E.R. (2020). Vaccine optimization for COVID-19: Who to vaccinate first? Science Advances, 7(6). 13. 13.Official statement for the delay of 2nd dose vaccines United Kingdom: [https://www.gov.uk/government/news/statement-from-the-uk-chief-medical-officers-on-the-prioritisation-of-first-doses-of-covid-19-vaccines](https://www.gov.uk/government/news/statement-from-the-uk-chief-medical-officers-on-the-prioritisation-of-first-doses-of-covid-19-vaccines) 14. 14.Official statement for the delay of 2nd dose vaccines Canada: [https://www.canada.ca/en/public-health/services/immunization/national-advisory-committee-on-immunization-naci/recommendations-use-covid-19-vaccines.html](https://www.canada.ca/en/public-health/services/immunization/national-advisory-committee-on-immunization-naci/recommendations-use-covid-19-vaccines.html) 15. 15.Matrajt, L., Eaton, J., Leung, T., Dimitrov, D., Schiffer, J.T., Swan, D.A., et al. (2021). Optimizing vaccine allocation for COVID-19 vaccines: potential role of single-dose vaccination. medRxiv, 2020.2012.2031. 20249099. 16. 16.Paltiel, A. D., Schwartz, J. L., Zheng, A., Walensky, R. P. (2021a). Clinical Outcomes of a COVID-19 Vaccine: Implementation Over Efficacy. Health Affairs, 10–1377. 17. 17.Paltiel, A.D., Zheng, A., Schwartz, J.L. (2021b). Speed Versus Efficacy: Quantifying Potential Tradeoffs in COVID-19 Vaccine Deployment. Annals of Internal Medicine, Volume 174, Issue 4, pp. 568–570. 18. 18.Silva, P. J. S., Sagastizábal, C, Nonato, L. G., Struchiner, C. J., Pereira, T., (2021) Optimized delay of the second COVID-19 vaccine dose reduces ICU admissions. Proceedings of the National Academy of Sciences. 19. 19.Ferreira et al. (2021). Assessing optimal time between doses in two-dose vaccination regimen in an ongoing epidemic of SARS-CoV-2. medRxiv 2021.07.28.21261200. 20. 20.Moore, S. & Hill, E., Tildesley, M., Dyson, L., Keeling, M., (2021). Vaccination and non-pharmaceutical interventions for COVID-19: a mathematical modelling study. The Lancet Infectious Diseases. 21. doi:10.1016/S1473-3099(21)00143-2. [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1016/S1473-3099(21)00143-2&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=33743847&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2021%2F11%2F20%2F2021.05.19.21257486.atom) 21. 21.Diekmann, O., Heesterbeek, H., Britton, T., Mathematical Tools for Understanding Infectious Disease Dynamics. Princeton University Press, pp. 516, ISBN 978-0-691-15539-5. 22. 22.Parino, F, Zino, L, Calafiore, GC, Rizzo, A. A model predictive control approach to optimally devise a two-dose vaccination rollout: A case study on COVID-19 in Italy. Int J Robust Nonlinear Control. 2021; 1–16. [https://doi.org/10.1002/rnc.5728](https://doi.org/10.1002/rnc.5728). 23. 23.Li, R., et al. (2020). Substantial undocumented infection facilitates the rapid dissemination of novel coronavirus (SARS-CoV2). Science, Vol. 368, Issue 6490, pp. 489–493. [Abstract/FREE Full Text](http://medrxiv.org/lookup/ijlink/YTozOntzOjQ6InBhdGgiO3M6MTQ6Ii9sb29rdXAvaWpsaW5rIjtzOjU6InF1ZXJ5IjthOjQ6e3M6ODoibGlua1R5cGUiO3M6NDoiQUJTVCI7czoxMToiam91cm5hbENvZGUiO3M6Mzoic2NpIjtzOjU6InJlc2lkIjtzOjEyOiIzNjgvNjQ5MC80ODkiO3M6NDoiYXRvbSI7czo1MDoiL21lZHJ4aXYvZWFybHkvMjAyMS8xMS8yMC8yMDIxLjA1LjE5LjIxMjU3NDg2LmF0b20iO31zOjg6ImZyYWdtZW50IjtzOjA6IiI7fQ==) 24. 24.Lauer, S.A., Grantz, K.H., Bi, Q. et al. (2020). The Incubation Period of Coronavirus Disease 2019 (COVID-19) From Publicly Reported Confirmed Cases: Estimation and Application. Annals of Internal Medicine, 172(9): 577-582. 25. 25.Liu, T., et al. (2020) Transmission Dynamics of 2019 Novel Coronavirus (2019-nCoV) (1/25/2020). Available at SSRN: [https://ssrn.com/abstract=3526307](https://ssrn.com/abstract=3526307) or [http://dx.doi.org/10.2139/ssrn.3526307](http://dx.doi.org/10.2139/ssrn.3526307). 26. 26.He, X., Lau, E.H.Y., Wu, P. et al. (2020). Temporal dynamics in viral shedding and transmissibility of COVID-19. Nature Medicine 26(5):672–675. [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.7326/M20-3012&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=http://www.n&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2021%2F11%2F20%2F2021.05.19.21257486.atom) 27. 27.Ganyani, T., Kremer, C., Chen, D., et al. (2020). Estimating the generation interval for coronavirus disease (COVID-19) based on symptom onset data. Eurosurveillance, 25(17):2000257. 28. 28.Li Q., Guan X., Wu P., Wang X., Zhou L., Tong Y., et al. (2020). Early Transmission Dynamics in Wuhan, China, of Novel Coronavirus-Infected Pneumonia. New England Journal of Medicine, 382(13):1199–207. [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1056/NEJMoa2001316&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=http://www.n&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2021%2F11%2F20%2F2021.05.19.21257486.atom) 29. 29.Cereda, D. et al. (2020). The early phase of the COVID-19 outbreak in Lombardy, Italy. arXiv e-prints; arxiv:2003.09320. 30. 30.Bi, Q. et al. (2020). Epidemiology and transmission of COVID-19 in 391 cases and 1286 of their close contacts in Shenzhen, China: a retrospective cohort study. The Lancet Infectious Diseases, 20(8): 911–919. [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1016/s1473-3099(20)30287-5&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=http://www.n&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2021%2F11%2F20%2F2021.05.19.21257486.atom) 31. 31.Lavvezzo, E., Franchin, E., Ciavarella, C., Cuomo-Dannenburg, G., Barzon, L., Del Vecchio, C., et al. (2020). Suppression of a SARS-CoV-2 outbreak in the Italian municipality of Vo’. Nature. 584, 425–429. [PubMed](http://medrxiv.org/lookup/external-ref?access_num=http://www.n&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2021%2F11%2F20%2F2021.05.19.21257486.atom) 32. 32.de Gier, B., et al. (2021). Vaccine effectiveness against SARS-CoV-2 transmission to household contacts during dominance of Delta variant (B.1.617.2), August-September 2021, the Netherlands, medRxiv: doi:10.1101/2021.10.14.21264959. [Abstract/FREE Full Text](http://medrxiv.org/lookup/ijlink/YTozOntzOjQ6InBhdGgiO3M6MTQ6Ii9sb29rdXAvaWpsaW5rIjtzOjU6InF1ZXJ5IjthOjQ6e3M6ODoibGlua1R5cGUiO3M6NDoiQUJTVCI7czoxMToiam91cm5hbENvZGUiO3M6NzoibWVkcnhpdiI7czo1OiJyZXNpZCI7czoyMToiMjAyMS4xMC4xNC4yMTI2NDk1OXYxIjtzOjQ6ImF0b20iO3M6NTA6Ii9tZWRyeGl2L2Vhcmx5LzIwMjEvMTEvMjAvMjAyMS4wNS4xOS4yMTI1NzQ4Ni5hdG9tIjt9czo4OiJmcmFnbWVudCI7czowOiIiO30=) 33. 33.Singanayagam, A., et al. (2021). Community transmission and viral load kinetics of the SARS-CoV-2 delta (B.1.617.2) variant in vaccinated and unvaccinated individuals in the UK: a prospective, longitudinal, cohort study, The Lancet Infectious Diseases, Elsevier: 1473–3099, [https://doi.org/10.1016/S1473-3099(21)00648-4](https://doi.org/10.1016/S1473-3099(21)00648-4). 34. 34.CDC (2021). Covid-19 Pandemic and Planning Scenarios, US Centers for Disease Control and Prevention [https://www.cdc.gov/coronavirus/2019-ncov/hcp/planning-scenarios.html](https://www.cdc.gov/coronavirus/2019-ncov/hcp/planning-scenarios.html), March 19, 2021 update. 35. 35.Haas, J.E., Angulo, F.J., McLaughlin, J.M. et al. (2021). Impact and effectiveness of mRNA BNT162b2 vaccine against SARS-CoV-2 infections and COVID-19 cases, hospitalisations, and deaths following a nationwide vaccination campaign in Israel: an observational study using national surveillance. The Lancet, 397(10287): 1819–1829. 36. 36.Sypsa, V., Roussos, S., Engeli, V., Paraskevis, D., Tsiodras, S., Hatzakis, A. (2021) Trends in COVID-19 vaccination intent, determinants and reasons for vaccine hesitancy: results from repeated cross-sectional surveys in the adult general population of Greece during November 2020-June 2021, medRxiv 2021.06.23.21259376. 37. 37.Sypsa, V., et al. (2021). Effects of Social Distancing Measures during the First Epidemic Wave of Severe Acute Respiratory Syndrome Infection, Greece. Emerging Infectious Diseases. 27(2): 452–462. [PubMed](http://medrxiv.org/lookup/external-ref?access_num=http://www.n&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2021%2F11%2F20%2F2021.05.19.21257486.atom) 38. 38.Amit, S., Regev-Yochay, G., Afek, A., Kreiss, Y., Leshem, E. (2021). Early rate reductions of SARS-CoV-2 infection and COVID-19 in BNT162b2 vaccine recipients. The Lancet, 397 (10277), 875–877. 39. 39.Thompson, M.G. et al. (2021). Interim Estimates of Vaccine Effectiveness of BNT162b2 and mRNA-1273 COVID-19 Vaccines in Preventing SARS-CoV-2 Infection Among Health Care Personnel, First Responders, and Other Essential and Frontline Workers - Eight U.S. Locations, December 2020-March 2021. MMWR. Morbidity and Mortality Weekly Report 70(13), 495–500. [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.15585/MMWR.MM7013E3&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=http://www.n&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2021%2F11%2F20%2F2021.05.19.21257486.atom) 40. 40.Tuite, A.R, Zhu, L., Fisman, D.N., Salomon, J.A. (2021). Alternative dose allocation strategies to increase benefits from constrained COVID-19 vaccine supply. Annals of internal medicine, 174(4), 570–572. 41. 41.Vasileiou, E., et al. (2021). Effectiveness of first dose of COVID-19 vaccines against hospital admissions in Scotland: national prospective cohort study of 5.4 million people. The Lancet preprints, February 2021, available at SSRN: [https://ssrn.com/abstract=3789264](https://ssrn.com/abstract=3789264). 42. 42.Moghadas SM, Vilches TN, Zhang K, Nourbakhsh S, Sah P, et al. (2021) Evaluation of COVID-19 vaccination strategies with a delayed second dose. PLOS Biology 19(4): e3001211. [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1371/journal.pbio.3001211&link_type=DOI) 43. 43.Romero-Brufau, S., (2021) Public health impact of delaying second dose of BNT162b2 or mRNA-1273 covid-19 vaccine: simulation agent based modeling study. 44. 44.Koh, W.C. et al. (2020). What do we know about SARS-CoV-2 transmission? A systematic review and meta-analysis of the secondary attack rate and associated risk factors. PLoS ONE 15(10): e0240205. [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1371/journal.pone.0240205&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=33031427&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2021%2F11%2F20%2F2021.05.19.21257486.atom) 45. 45.ELSTAT (2021), Greek Statistics Authority, [https://www.statistics.gr/](https://www.statistics.gr/). 46. 46.Skowronski D.M., De Serres G. (2020). Safety and efficacy of the BNT162b2 mRNA Covid-19 vaccine. New England Journal of Medicine, 2021. 47. 47.Baden, L. R. et al. (2021). Efficacy and safety of the mRNA-1273 SARS-CoV-2 vaccine. New England. Journal of Medicine, 384:403–416. [1]: /embed/inline-graphic-1.gif