Epidemiological impact of SARS-CoV-2 vaccination: mathematical modeling analyses ================================================================================ * Monia Makhoul * Houssein H. Ayoub * Hiam Chemaitelly * Shaheen Seedat * Ghina R Mumtaz * Sarah Al-Omari * Laith J. Abu-Raddad ## Abstract **Background** Several SARS-CoV-2 vaccine candidates are currently in the pipeline. This study aims to inform SARS-CoV-2 vaccine development, licensure, decision-making, and implementation by determining key preferred vaccine product characteristics and associated population-level impact. **Methods** Vaccination impact was assessed at various efficacies using an age-structured mathematical model describing SARS-CoV-2 transmission and disease progression, with application for China. **Results** A prophylactic vaccine with efficacy against acquisition (*VE**S*) of ≥70% is needed to eliminate this infection. A vaccine with *VE**S* <70% will still have a major impact, and may control the infection if it reduces infectiousness or infection duration among those vaccinated who acquire the infection, or alternatively if supplemented with a moderate social-distancing intervention (<20% reduction in contact rate), or complemented with herd immunity. Vaccination is cost-effective. For a vaccine with *VE**S* of 50%, number of vaccinations needed to avert one infection is only 2.4, one severe disease case is 25.5, one critical disease case is 33.2, and one death is 65.1. Gains in effectiveness are achieved by initially prioritizing those ≥60 years. Probability of a major outbreak is virtually zero with a vaccine with *VE**S* ≥70%, regardless of number of virus introductions. Yet, an increase in social contact rate among those vaccinated (behavior compensation) can undermine vaccine impact. **Conclusions** Even a partially-efficacious vaccine can offer a fundamental solution to control SARS-CoV-2 infection and at high cost-effectiveness. In addition to the primary endpoint on infection acquisition, developers should assess natural history and disease progression outcomes and/or proxy biomarkers, since such secondary endpoints may prove critical in licensure, decision-making, and vaccine impact. Keywords * SARS-CoV-2 * COVID-19 * coronavirus * epidemiology * vaccine * mathematical model ## Introduction Following the Severe Acute Respiratory Syndrome (SARS) epidemic in 2002 and the Middle East Respiratory Syndrome (MERS) epidemic in 2012 [1], a novel coronavirus, severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2), emerged in late December 2019 in Wuhan, Hubei province, China [2, 3]. While the earlier coronavirus epidemics were rather limited in scope and scale [1], SARS-CoV-2 rapidly spread [4] and evolved into a pandemic [5]. In absence of a vaccine, even if partially efficacious [6], containment of the epidemic in China necessitated large-scale contact tracing and testing through deployment of thousands of healthcare fieldworkers along with severe quarantine measures [4]. The strain put on healthcare systems [7], and the global human [8, 9] and economic [10] losses caused by the virus and resulting disease, designated as Coronavirus Disease 2019 (COVID-2019) [11], accelerated efforts towards vaccine development [6, 12]. While multiple vaccine candidates are currently in the pipeline, they are still in early stage of development [6, 12, 13]. Assessment of the population-level impact of vaccine candidates through mathematical modeling is a critical component of the process of vaccine development, value proposition, licensure, decision-making, and pathways and costs of vaccine administration, and has been utilized for a wide range of infectious diseases [14-28]. In early stages of development, modeling is used to define the vaccine’s key preferred product characteristics, by estimating levels of efficacy necessary to observe significant population-level impact, determining necessary duration of protection/immunity incurred by the vaccine, and identifying priority populations for optimal effectiveness [21, 29, 30]. These parameters provide early guidance to developers, manufacturers, regulators, and decision makers about candidates that are likely to be optimal through specifying vaccine characteristics that will maximize public health impact and cost-effectiveness [21, 28, 29, 31, 32]. Once key attributes are established, modeling plays an integral role in building the case for investment in vaccine development, and in ensuring rapid roll-out post-licensing, through assessment of risks, costs, and predicted returns associated with different immunization strategies [29, 33]. Post-vaccination, modeling is used to inform design and interpretation of surveillance studies [25-27]. We aimed in this study to provide the scientific evidence necessary to inform and accelerate SARS-CoV-2 vaccine development, licensure, decision-making, and implementation by determining key preferred vaccine product characteristics and associated population-level impact, at a critical time for such development [6, 12, 13]. ## Methods ### Mathematical Model A deterministic model was constructed to describe SARS-CoV-2 transmission dynamics in a given population, here China as an illustrative example, in presence of vaccination (Figure S1 of Supplementary Material (SM)). The model extended a recently-developed age-structured model focused on analyzing SARS-CoV-2 epidemiology in China [34], the country with the most complete epidemic cycle to date. Model structure was informed by current understanding of SARS-CoV-2 natural history and epidemiology, and consisted of a set of coupled nonlinear differential equations that stratify the population into compartments based on vaccination status, age group, infection status, infection stage, and disease stage (Section A of SM). For each of vaccinated and unvaccinated populations, nine age groups were considered, each representing a 10-year age band except for the last category (0-9, 10-19, …, ≥80 years). Susceptible individuals were at risk of being exposed to the infection at varying hazard rates depending on their age group and vaccination status. Following a latency period, infected individuals develop asymptomatic or mild infection followed by recovery, or severe infection followed by severe disease and then recovery, or critical infection followed by critical disease and either recovery or disease mortality. Mixing between individuals of different age groups was determined by an age mixing matrix that allows a range of mixing. Details on model structure are in SM. The model was coded, fitted, and analysed using MATLAB R2019a [35]. ### Model parametrization and fitting The model was parameterized and calibrated using empirical data on SARS-CoV-2 natural history and epidemiology. Age-specific distributions of infected individuals across the mild, severe, or critical infection stages were based on case-severity levels observed in China [4, 36, 37]. Critical disease cases were at risk of disease mortality, with the relative mortality rate in each age group informed by the age-specific crude case fatality rate observed in China [3, 38]. Population size, demographic structure (age distribution), and life expectancy, as of 2020, were obtained from the United Nations World Population Prospects database [39]. The model was fitted to empirical time series data for the daily and cumulative numbers of diagnosed SARS-CoV-2 cases and deaths [40], number of recovered individuals [40], as well as the age-specific attack rate [38, 41]. Details of model parameters, values, and justifications are in Tables S1 and S2 and Section B of SM. ### Product characteristics of candidate vaccines We assessed the impact of a prophylactic vaccine that reduces susceptibility to infection. However, since the first available vaccine may only be partially efficacious against infection acquisition, we also assessed the impact of the vaccine assuming additional “breakthrough” effects, that is effects that modulate natural history of infection for those vaccinated but still acquire the infection. Specifically, we assumed that vaccination may reduce infectiousness per one contact (by reducing viral load), infection duration (by faster clearance with vaccine-induced immunity), and likelihood of developing severe or critical disease (by rapid immune response that prevents disease progression). Definitions of these vaccine efficacies are summarized in Table 1. View this table: [Table 1.](http://medrxiv.org/content/early/2020/04/23/2020.04.19.20070805/T1) Table 1. Key vaccine product characteristics used to assess impact of a vaccine against SARS-CoV-2 Other relevant characteristics include duration of protection elicited by the vaccine and vaccination effect on adherence to social distancing—we investigated impact of increasing social contact rate following vaccination with the perception of protection. ### Measures of vaccine impact Population-level impact of SARS-CoV-2 vaccination was assessed by quantifying incidence, cumulative incidence, and reduction in incidence of each of infections, severe disease cases, critical disease cases, and deaths arising in presence of vaccination compared to the counter-factual scenario of no-vaccination. Vaccination impact was further assessed by quantifying *effectiveness*, that is number of vaccinations needed to avert one infection or one adverse disease outcome (ratio of number of vaccinations relative to that of averted outcomes). The latter measure is essentially cost-effectiveness with no costs included, as they are not yet available. Vaccination impact was assessed at: 1) *VES* = 50% but ![Graphic][1] but ![Graphic][2] but ![Graphic][3] but ![Graphic][4] and 5) ![Graphic][5]but. Vaccine was assumed to elicit protection over 10 years. ### Vaccination program scenarios Two vaccination program scenarios were considered. In both programs, it was assumed that vaccination is introduced in absence of social-distancing interventions, as the purpose of vaccination is to replace such interventions. The first program scenario assumes vaccine introduction and scale-up to 80% coverage before epidemic onset. This scenario is relevant for assessing impact of vaccination on future SARS-CoV-2 introductions in countries where the epidemic has been contained or at low level, such as in China. The scenario is also relevant to assess *maximum* potential impact of vaccination regardless of current epidemic status. The second program scenario assumes vaccine introduction during the epidemic’s exponential growth phase, with scale-up to 80% coverage within one month. ### Additional analyses Incidence of new infections was assessed at various levels of *VE**S* to determine the *minimum* efficacy needed to fully control the infection. Incidence was also assessed in a scenario where vaccination was introduced with a social-distancing intervention to estimate level of social distancing needed to complement vaccination to control the infection. Incidence was assessed in another scenario where those vaccinated increased their social contacts (behavior compensation), to assess consequences on vaccination impact. Lastly, we derived and estimated likelihood of occurrence of a major outbreak following infection introduction in a vaccinated but infection-free population (Section D of SM). ### Uncertainty analysis A multivariable uncertainty analysis was conducted to determine range of uncertainty around model predictions using five-hundred model runs. At each run, Latin Hypercube sampling [42, 43] was applied in selecting the natural history and disease progression parameter values from ranges specified by assuming ±30% uncertainty around parameters’ point estimates. The model was then refitted to input data and vaccine impact assessed in the new fitted model. Resulting distribution for vaccine impact across all 500 runs was used to calculate predicted means and 95% uncertainty intervals (UIs). ## Results Figures 1 and 2 illustrate impact of vaccination assuming different vaccine product characteristics (efficacies) for each vaccination program roll-out scenario. In the first scenario (Figure 1 and Figures S2 and S3 of SM), where vaccination was scaled up to 80% coverage before epidemic onset, the epidemic in absence of vaccination peaked at 158 days after virus introduction, but at 286 days when ![Graphic][6], 452 days when *VE**I* = 50%, and 462 days when *VE**S* = 50%. There was no epidemic when ![Graphic][7]. A vaccine with *VE**S* = 50% reduced peak infection incidence by 84.4% and cumulative/total infections by 52.8%, peak severe disease incidence by 83.9% and cumulative severe disease cases by 53.4%, peak critical disease incidence by 82.1% and cumulative critical disease cases by 46.7%, peak death incidence by 79.0% and cumulative deaths by 44.4%. A vaccine with *VE**I* = 50% yielded slightly lower reductions in incidence of infection and adverse outcomes, while a vaccine with ![Graphic][8] was less impactful though still achieved considerable reductions. A vaccine with ![Graphic][9] had no impact on infection incidence, but reduced peak incidence of each of severe and critical disease by 38.5%, and deaths by 40.0% (Figure 1), and the cumulative incidence of the latter three outcomes by ∼39% (Figures S2 and S3 of SM). ![Figure 1.](http://medrxiv.org/https://www.medrxiv.org/content/medrxiv/early/2020/04/23/2020.04.19.20070805/F1.medium.gif) [Figure 1.](http://medrxiv.org/content/early/2020/04/23/2020.04.19.20070805/F1) Figure 1. Impact of SARS-CoV-2 vaccination on number of A) new infections, B) new severe disease cases, C) new critical disease cases, and D) new deaths in the scenario assuming vaccine scale-up to 80% coverage before epidemic onset. Duration of vaccine protection is 10 years. Impact was assessed at ![Graphic][10]. ![Figure 2.](http://medrxiv.org/https://www.medrxiv.org/content/medrxiv/early/2020/04/23/2020.04.19.20070805/F2.medium.gif) [Figure 2.](http://medrxiv.org/content/early/2020/04/23/2020.04.19.20070805/F2) Figure 2. Impact of SARS-CoV-2 vaccination on number of A) new infections, B) new severe disease cases, C) new critical disease cases, and D) new deaths in the scenario assuming vaccine introduction during the exponential growth phase of the epidemic, with scale-up to 80% coverage within one month. Duration of vaccine protection is 10 years. Impact was assessed at ![Graphic][11]. In the second scenario (Figure 2 and Figures S4 and S5 of SM), where vaccination was rapidly scaled up to 80% coverage during the exponential growth phase, the epidemic peaked earlier and at lower values for incidence of infection and adverse outcomes. Impact of a vaccine with *VE**S* = 50% was initially similar to that of a vaccine with ![Graphic][12], however, over time, the latter was more impactful in reducing infection and adverse outcomes. Reduction in cumulative number of new infections (at end of epidemic cycle) was highest for ![Graphic][13] at 53.4%, followed by *VE**S* = 50% at 41.2%, *VE**I* = 50% at 28.2%, and ![Graphic][14] at 23.1%, with no reduction for ![Graphic][15]. Reduction in cumulative number of new deaths for these efficacies was, respectively, 47.2%, 34.8%, 22.5%, 18.0%, and 30.0%. Figure 3 illustrates vaccine effectiveness in averting infection and adverse outcomes by end of epidemic cycle for first program scenario. For *VE**S* = 50%, 2.4 vaccinations were needed to avert one infection, 25.4 to avert one severe disease case, 33.2 to avert one critical disease case, and 65.1 to avert one death. Effectiveness was nearly comparable for *VE**I* = 50%, whereas more vaccinations were needed to avert one infection or one adverse outcome for ![Graphic][16] and ![Graphic][17]. Best effectiveness was for ![Graphic][18] where only 1.3 vaccinations were needed to avert one infection, 13.6 to avert one severe disease case, 15.5 to avert one critical disease case, and 28.9 to avert one death. Graphs illustrating temporal evolution of vaccine effectiveness for vaccination program scenarios 1 and 2 are in Figures S6 and S7 of SM, respectively. ![Figure 3.](http://medrxiv.org/https://www.medrxiv.org/content/medrxiv/early/2020/04/23/2020.04.19.20070805/F3.medium.gif) [Figure 3.](http://medrxiv.org/content/early/2020/04/23/2020.04.19.20070805/F3) Figure 3. SARS-CoV-2 vaccine effectiveness. Number of vaccinations needed to avert A) one infection, B) one severe disease case, C) one critical disease case, and D) one death, by end of epidemic cycle. Scenario assumes vaccine scale-up to 80% coverage before epidemic onset. Duration of vaccine protection is 10 years. Impact was assessed at![Graphic][19]. Panel A does not include the result for ![Graphic][20], as this efficacy has no impact on number of infections—it affects only severe and critical disease and death. Figure 4 shows effectiveness of age-group prioritization by end of epidemic cycle using a vaccine with *VE**S* = 50%, introduced before epidemic onset. Prioritizing adults ≥20 years of age was most effective in reducing infection incidence, with ≤3 vaccinations needed to avert one infection. Meanwhile, prioritizing adults ≥60 years of age was most effective in reducing new deaths, with ≤36 vaccinations needed to avert one death. Prioritizing children was least effective with large number of vaccinations needed to avert one infection or one adverse outcome. Of note that there are minor differences in effectiveness over time. For instance, in initial phases of the epidemic, prioritizing those 60-69 years of age was slightly more effective than prioritizing those 40-49 years of age (Figure S8 of SM). Meanwhile, towards end of epidemic cycle, the inverse was true. ![Figure 4.](http://medrxiv.org/https://www.medrxiv.org/content/medrxiv/early/2020/04/23/2020.04.19.20070805/F4.medium.gif) [Figure 4.](http://medrxiv.org/content/early/2020/04/23/2020.04.19.20070805/F4) Figure 4. Effectiveness of age-group prioritization using a SARS-CoV-2 vaccine with *VE**s* of 50%. Number of vaccinations needed to avert A) one infection, B) one severe disease case, C) one critical disease case, and D) one death, by prioritizing different age groups for vaccination. Scenario assumes vaccine scale-up to 80% coverage before epidemic onset, and duration of vaccine protection of 10 years. Effectiveness is assessed at end of epidemic cycle. Figure 5A shows cumulative number of infections (at end of epidemic cycle) at various *VE**S* levels for a vaccine introduced before epidemic onset. A gradual decrease is observed as *VE**S* increases, with an accelerated reduction as *VE**S* approaches 60%—the level beyond which number of infections approaches zero. Epidemic onset is prevented at *VE**S* = 69%. Figure 5B illustrates the gains in effectiveness as *VE**S* increases, with 5.3 vaccinations needed to avert one infection at *VE**S* = 30%, but only 1.3 at *VE**S* = 60%. ![Figure 5.](http://medrxiv.org/https://www.medrxiv.org/content/medrxiv/early/2020/04/23/2020.04.19.20070805/F5.medium.gif) [Figure 5.](http://medrxiv.org/content/early/2020/04/23/2020.04.19.20070805/F5) Figure 5. Impact of varying levels of vaccine efficacy in reducing susceptibility (*VE**s*) on A) cumulative number of new SARS-CoV-2 infections and B) number of vaccinations needed to avert one SARS-CoV-2 infection. Scenario assumes vaccine scale-up to 80% coverage before epidemic onset. Duration of vaccine protection is 10 years. Measures are assessed at end of epidemic cycle. While a vaccine with *VE**S* = 50% cannot fully control the epidemic, Figure S9 of SM shows the impact when vaccination is supplemented with a social-distancing intervention that reduces the contact rate. Less than 20% reduction in contact rate would be sufficient to fully control the epidemic. Vaccinated individuals may increase their contact rate with the perception of protection. Figure 6 shows the consequences of behavior compensation. A 20% increase in contact rate among those vaccinated lowered the reduction in cumulative incidence from 52.8% to only 21.0%. A 41.8% increase in contact rate nullified impact of vaccination in reducing incidence. ![Figure 6.](http://medrxiv.org/https://www.medrxiv.org/content/medrxiv/early/2020/04/23/2020.04.19.20070805/F6.medium.gif) [Figure 6.](http://medrxiv.org/content/early/2020/04/23/2020.04.19.20070805/F6) Figure 6. Impact of vaccination with reduced adherence to social distancing for those vaccinated. Figure shows impact of varying levels of behavior compensation post-vaccination on the vaccine-induced reduction in cumulative number of new SARS-CoV-2 infections by end of epidemic cycle. Scenario assumes vaccine scale-up to 80% coverage before epidemic onset, *VE**s* is 50%, and duration of vaccine protection is 10 years. Figure 7 reports probability of occurrence of a major outbreak at varying levels of vaccine efficacy. A gradual decrease is noted with the increase in efficacy that accelerates close to *VE**S* or *VE**I* of ![Graphic][21] 85%, and ![Graphic][22] of ∼33% beyond which no major outbreak is expected to occur. Of note that this figure shows (conservatively) the upper bound of the probability of a major outbreak. Results for the lower bound are in Figure S10 of SM and derivations in Section D of SM. ![Figure 7.](http://medrxiv.org/https://www.medrxiv.org/content/medrxiv/early/2020/04/23/2020.04.19.20070805/F7.medium.gif) [Figure 7.](http://medrxiv.org/content/early/2020/04/23/2020.04.19.20070805/F7) Figure 7. Probability of occurrence of a major outbreak following vaccination. Probability of occurrence of a major outbreak upon virus introduction at varying levels of A) *VES*, B) *VEI*, C)![Graphic][23], and D)![Graphic][24]. Scenario assumes vaccine scale-up to 80% coverage before epidemic onset. Duration of vaccine protection is 10 years. Figure does not include the result for ![Graphic][25], as this efficacy has no impact on probability of occurrence of a major outbreak. Uncertainty analysis demonstrated robustness of model predictions to a wide range of uncertainty in input parameters (Figure S11 of SM). ## Discussion Above results indicate that even a partially-efficacious vaccine can offer a fundamental solution to the SARS-CoV-2 pandemic—the vaccine does not need to have sterilizing immunity to fully control the infection. Indeed, a vaccine with *VE**S* ≥70% could be sufficient to control the pandemic at ≥80% coverage (Figure 5A). Even a vaccine with *VE**S* <70% may still control the infection if it additionally (and plausibly) induces “breakthrough” effects such as reduction in viral load (reduction in infectiousness; *VE**I*) or faster infection clearance (reduction in infection duration; ![Graphic][26]) among those vaccinated who still acquire the infection. The latter effects have *individually* (that is in absence of protection against acquisition) a comparable impact on transmission to that of a prophylactic vaccine that reduces infection acquisition (Figures 1 and 2). Even in absence of such effects, infection control can still be achieved if vaccination is supplemented with only a moderate social-distancing intervention (Figure S9 of SM), or complemented with partial herd immunity—a considerable fraction of the population could have acquired the infection and developed protective antibodies by the time of vaccine roll-out. Even a vaccine that does not prevent infection yet only mitigates disease progression (reduction in severe or critical disease and death; ![Graphic][27]), could still yield significant gains by curbing disease burden (Figures 1 and 2). Results also indicated that vaccine impact depends on time of vaccine introduction whether before (Figure 1 and Figures S2 and S3 of SM) or after (Figure 2 and Figures S4 and S5 of SM) epidemic onset and/or growth—maximal gains are achieved with earlier introduction. Early introduction optimally defers epidemic growth, flattens the incidence curve, and reduces number of infections and disease outcomes (Figures 1 and 2 and Figures S2-S5 of SM). The vaccine will likely be cost-effective over a broad range of efficacy levels (Figure 5B). For a vaccine with *VE**S* of 50%, number of vaccinations needed to avert one infection is only 2.4, one severe disease case is 25.5, one critical disease case is 33.2, and one death is 65.1 (Figure 3). Return on effectiveness is also rapid for such respiratory infection with fast-growing epidemic scale (Figures S6 and S7 of SM). Effectiveness can be further enhanced by prioritizing vaccination for those ≥60 years of age, for an optimal reduction in disease cases and deaths (Figure 4). Conversely, prioritizing children is least effective, with their lower risk for developing adverse outcomes (Figure 4) [4, 44, 45]. This being said, prioritizing vaccination for any single age group, regardless of that age group, has overall lower effectiveness than extending vaccination to all age groups—vaccinating only one age group reduces the reproduction number (*R*) only marginally, whereas vaccinating all age groups reduces *R* to an epidemic domain where small reductions in *R* can have more substantial impact on epidemic size (Figure S12 of SM, also Figure 3 versus Figure 4). Consequently, roll-out strategies should initially prioritize individuals ≥60 years of age, but then incrementally cover younger age cohorts, and eventually the entire population. Vaccination will also reduce the likelihood of a major outbreak following virus introduction/reintroduction into a population (Figure 7). With a vaccine of *VE**S* ≥70%, infection transmission chains may not be sustainable, regardless of number of virus introductions. Of concern, however, is the potential increase in social contact rate among those vaccinated (behavior compensation)—a 42% increase in contact rate can virtually nullify the gains of a vaccine with *VE**S* of 50% at coverage of 80% (Figure 6). Roll-out of a vaccine with intermediate efficacy should be coupled with public health communication that stresses caution in social mixing following vaccination. This study has limitations. Model estimations are contingent on validity and generalizability of input data. While we used available evidence for SARS-CoV-2 natural history and epidemiology, our understanding of the epidemiology is still evolving. We assessed vaccine impact using China as an illustrative example, given the advanced epidemic cycle, yet evidence suggests that many infections may have been undocumented in this country, particularly in the early epidemic phase [46]. This may affect some estimates, such as mortality, probably towards overestimation [34, 47]. While the *absolute* impact of the vaccine on disease severity and mortality may have been overestimated, the *relative* impact (reduction rate) is less likely to have been affected. Our baseline *R* for China was 2.1 [34], but *R* may vary across settings, thus affecting estimates for the minimum efficacy needed for infection elimination. For instance, for an *R* of 3, the minimum *VE**S* needed for elimination is about 90% (Figure S12 of SM). We assessed vaccine impact for one epidemic cycle, with no assessment of seasonality or future cycles. We assumed a long duration of vaccine protection (10 years), but this has limited impact on the predictions for one epidemic cycle, provided the duration of vaccine protection is greater than one year. Despite these limitations, our model was complex enough to factor the different key vaccine product characteristics, but also parsimonious enough to be tailored to the nature of available data. The model also generated results that are valid to a wide range of model assumptions. ## Conclusions With most of the world’s population remaining susceptible to SARS-CoV-2 and the need to impose disruptive social-distancing interventions, vaccination is a reliable intervention in the long-term. Findings show that even a partially-efficacious vaccine provides a fundamental solution to the SARS-CoV-2 pandemic and at high cost-effectiveness. Vaccine impact and cost-effectiveness will not only depend on its efficacy in preventing infection, but can be enhanced if those vaccinated who still acquire the infection have reduced infectiousness, duration of infection, and disease severity. Vaccine developers should thus not only assess the primary endpoint of reduction in acquisition, but also other outcomes and/or proxy biomarkers including reductions in viral load and disease outcomes and speed of infection clearance for those vaccinated and unvaccinated. Totality of these primary and secondary endpoints may prove critical in the licensure process, decision-making, and vaccine impact once introduced into a population. ## Data Availability All data are available within the manuscript and its supplementary material. ## Contributors MM and HHA constructed, coded, and parameterized the mathematical model. MM conducted the analyses. HC supported model parametrization and analyses and wrote the first draft of the paper. LJA conceived and led the design of the study, construct and parameterization of the mathematical model, and drafting of the article. All authors contributed to discussion and interpretation of the results and to the writing of the manuscript. All authors have read and approved the final manuscript. ## Declaration of interests We declare no competing interests. ## Supplementary Material ### I. SARS-CoV-2 mathematical vaccine model #### A. Model structure We extended a recently-developed age-structured deterministic compartmental model to describe the impact of vaccination on severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) transmission dynamics and progression of the resulting disease, Coronavirus Disease 2019 (COVID-2019), in a given population. The model stratifies the unvaccinated and vaccinated populations into compartments according to age group (0-9, 10-19, 20-29,…, ≥80 years), infection status (uninfected, infected), infection stage (mild, severe, critical), and disease stage (severe, critical). Transmission and disease progression dynamics in vaccinated and unvaccinated cohorts are described using age-specific sets of nonlinear ordinary differential equations, where each age group *a* (*a* = 1, 2,…9) refers to a 10-year age band (0-9,10-19,..70-79) apart from the last group including all those aged ≥ 80 years. The model is illustrated in Figure S1. ![Figure S1.](http://medrxiv.org/https://www.medrxiv.org/content/medrxiv/early/2020/04/23/2020.04.19.20070805/F8.medium.gif) [Figure S1.](http://medrxiv.org/content/early/2020/04/23/2020.04.19.20070805/F8) Figure S1. Schematic diagram describing the basic structure of the SARS-CoV-2 vaccine model. The following set of equations describes the transmission dynamics among unvaccinated and vaccinated populations in the first age group: #### Unvaccinated population aged 0-9 years ![Formula][28] #### Vaccinated population aged 0-9 years ![Formula][29] For subsequent age groups, the following set of equations was used: #### Unvaccinated populations aged 10+ years ![Formula][30] #### Vaccinated populations aged 10+ years ![Formula][31] The definitions of population variables and symbols used in the equations are in Table S1. View this table: [Table S1.](http://medrxiv.org/content/early/2020/04/23/2020.04.19.20070805/T2) Table S1. Definitions of population variables and symbols used in the model. The force of infection (hazard rate of infection) experienced by the unvaccinated susceptible populations *S* (*a*) is given by ![Formula][32] while that of vaccinated susceptible populations *V* (*a*) is given by ![Formula][33] where *β* is the overall infectious contact rate. The mixing among the different age groups is dictated by the mixing matrix ℋ*a,a*”. This matrix provides the probability that an individual in the *a* age group will mix with an individual in the *a*′ age group (regardless of vaccination status). The mixing matrix is given by ![Formula][34] Here, *δ**a,a*′ is the identity matrix. *e**Age* *∈* [0,1] measures the degree of assortativeness in the mixing. At the extreme *e**Age* = 0, the mixing is fully proportional. Meanwhile, at the other extreme, *e**Age* = 1, the mixing is fully assortative, that is individuals mix only with members in their own age group. #### B. Parameter values The input parameters of the model were chosen based on current empirical data for SARS-CoV-2 natural history and epidemiology. The parameter values are listed in Table S2. View this table: [Table S2.](http://medrxiv.org/content/early/2020/04/23/2020.04.19.20070805/T3) Table S2. Model assumptions in terms of parameter values. #### C. The basic reproduction number *R* Using the second generation matrix method described by Heffernan *et al*. [12], the basic reproduction number in absence of vaccination is given by ![Formula][35] Where *w*(*a*) is the proportion of the population in each age group. Meanwhile, the basic reproduction number in presence of vaccination is given by ![Formula][36] Where *z**unvac* (*a*) and *z**vac* (*a*) are, respectively, the proportions of unvaccinated and vaccinated populations in each age group. #### D. Probability of a major outbreak Based on Whittle’s method [13], and by constructing Bailey’s ratios [13], the probability of a major outbreak was derived, that is the probability that the fraction of susceptible individuals that become infected is ≥ *γ*, where *γ* is a specific chosen level of the final attack rate. We defined *ρ**U**M, N* (*a*), *ρ**U* *S, N* (*a*), and *ρ**U* *C, N* (*a*) to be, respectively, Bailey’s ratios corresponding to the unvaccinated infectious classes *I**M* (*a*), *I**S* (*a*), and *I**C* (*a*) for each age group. Similarly, we defined *ρ**V* *M, N* (*a*), *ρ**V* *S, N* (*a*), and *ρ**V* *C, N* (*a*) to be, respectively, Bailey’s ratios corresponding to the vaccinated infectious classes ![Graphic][37], and ![Graphic][38] for each age group. The number of initial infections is denoted by *b*. We define ![Formula][39] In the case of no vaccination, the probability of a major outbreak is given by the following three cases: #### Case 1 For ![Graphic][40], and ![Graphic][41], the probability of a major outbreak lies between ![Graphic][42] and ![Graphic][43] #### Case 2 For ![Graphic][44], and ![Graphic][45], the probability of a major outbreak lies between 0 and ![Graphic][46]. #### Case 3 For ![Graphic][47], the probability of a major outbreak is 0. Meanwhile, in the case of vaccinating a fraction of the population, the probability of a major outbreak is given by the following three cases: #### Case 1 For ![Graphic][48] and ![Graphic][49], the probability of a major outbreak lies between ![Graphic][50] and ![Graphic][51]. #### Case 2 For ![Graphic][52] and ![Graphic][53], the probability of a major outbreak lies between 0 and ![Graphic][54]. #### Case 3 For ![Graphic][55], the probability of a major outbreak is 0. ![Figure S2.](http://medrxiv.org/https://www.medrxiv.org/content/medrxiv/early/2020/04/23/2020.04.19.20070805/F9.medium.gif) [Figure S2.](http://medrxiv.org/content/early/2020/04/23/2020.04.19.20070805/F9) Figure S2. Impact of SARS-CoV-2 vaccination on the cumulative number of A) new infections, B) new severe disease cases, C) new critical disease cases, and D) new deaths in the scenario assuming vaccine scale-up to 80% coverage before epidemic onset. Duration of vaccine protection is 10 years. Impact was assessed at ![Graphic][56]. ![Figure S3.](http://medrxiv.org/https://www.medrxiv.org/content/medrxiv/early/2020/04/23/2020.04.19.20070805/F10.medium.gif) [Figure S3.](http://medrxiv.org/content/early/2020/04/23/2020.04.19.20070805/F10) Figure S3. Role of SARS-CoV-2 vaccination in reducing the cumulative number of A) new infections, B) new severe disease cases, C) new critical disease cases, and D) new deaths in the scenario assuming vaccine scale-up to 80% coverage before epidemic onset. Duration of vaccine protection is 10 years. Impact was assessed at ![Graphic][57]. ![Figure S4.](http://medrxiv.org/https://www.medrxiv.org/content/medrxiv/early/2020/04/23/2020.04.19.20070805/F11.medium.gif) [Figure S4.](http://medrxiv.org/content/early/2020/04/23/2020.04.19.20070805/F11) Figure S4. Impact of SARS-CoV-2 vaccination on the cumulative number of A) new infections, B) new severe disease cases, C) new critical disease cases, and D) new deaths in the scenario assuming vaccine introduction during the exponential growth phase of the epidemic, with scale-up to 80% coverage within one month. Duration of vaccine protection is 10 years. Impact was assessed at ![Graphic][58]. ![Figure S5.](http://medrxiv.org/https://www.medrxiv.org/content/medrxiv/early/2020/04/23/2020.04.19.20070805/F12.medium.gif) [Figure S5.](http://medrxiv.org/content/early/2020/04/23/2020.04.19.20070805/F12) Figure S5. Role of SARS-CoV-2 vaccination in reducing the cumulative number of A) new infections, B) new severe disease cases, C) new critical disease cases, and D) new deaths in the scenario assuming vaccine introduction during the exponential growth phase of the epidemic, with scale-up to 80% coverage within one month. Duration of vaccine protection is 10 years. Impact was assessed at ![Graphic][59]. ![Figure S6.](http://medrxiv.org/https://www.medrxiv.org/content/medrxiv/early/2020/04/23/2020.04.19.20070805/F13.medium.gif) [Figure S6.](http://medrxiv.org/content/early/2020/04/23/2020.04.19.20070805/F13) Figure S6. Temporal evolution of SARS-CoV-2 vaccine effectiveness in the scenario assuming vaccine scale-up to 80% coverage before epidemic onset. Number of vaccinations needed to avert A) one new infection, B) one new severe disease case, C) one new critical disease case, and D) one new death, depending on time into the epidemic. Duration of vaccine protection is 10 years. Impact was assessed ![Graphic][60]. Panel A does not include the result for ![Graphic][61], as this efficacy has no impact on number of infections—it affects only severe and critical disease and death. ![Figure S7.](http://medrxiv.org/https://www.medrxiv.org/content/medrxiv/early/2020/04/23/2020.04.19.20070805/F14.medium.gif) [Figure S7.](http://medrxiv.org/content/early/2020/04/23/2020.04.19.20070805/F14) Figure S7. Temporal evolution of SARS-CoV-2 vaccine effectiveness in the scenario assuming vaccine introduction during the exponential growth phase of the epidemic, with scale-up to 80% coverage within one month. Number of vaccinations needed to avert A) one new infection, B) one new severe disease case, C) one new critical disease case, and D) one new death, depending on time into the epidemic. Duration of vaccine protection is 10 years. Impact was assessed at ![Graphic][62]. Panel A does not include the result for ![Graphic][63], as this efficacy has no impact on number of infections—it affects only severe and critical disease and death. ![Figure S8.](http://medrxiv.org/https://www.medrxiv.org/content/medrxiv/early/2020/04/23/2020.04.19.20070805/F15.medium.gif) [Figure S8.](http://medrxiv.org/content/early/2020/04/23/2020.04.19.20070805/F15) Figure S8. Temporal evolution of effectiveness of age-group prioritization using a SARS-CoV-2 vaccine with *VE**s* of 50%. Number of vaccinations needed to avert A) one infection and B) one death, by prioritizing individuals aged 40-49 or 60-69 years. Scenario assumes vaccine scale-up to 80% coverage before epidemic onset. Duration of vaccine protection is 10 years. ![Figure S9.](http://medrxiv.org/https://www.medrxiv.org/content/medrxiv/early/2020/04/23/2020.04.19.20070805/F16.medium.gif) [Figure S9.](http://medrxiv.org/content/early/2020/04/23/2020.04.19.20070805/F16) Figure S9. Impact of a social-distancing intervention reducing the contact rate in the population on the cumulative number of new SARS-CoV-2 infections, when introduced to supplement the impact of a vaccine that has 50% efficacy in reducing susceptibility, *VE**S*. Scenario assumes vaccine scale-up to 80% coverage before epidemic onset. Duration of vaccine protection is 10 years. ![Figure S10.](http://medrxiv.org/https://www.medrxiv.org/content/medrxiv/early/2020/04/23/2020.04.19.20070805/F17.medium.gif) [Figure S10.](http://medrxiv.org/content/early/2020/04/23/2020.04.19.20070805/F17) Figure S10. Probability of occurrence of a major outbreak following vaccination. Upper (solid line) and lower (dashed line) bounds of the probability of occurrence of a major outbreak upon virus introduction at varying levels of A) *VES*, B) *VEI*, C)![Graphic][64], and D) ![Graphic][65]. Scenario assumes vaccine scale-up to 80% coverage before epidemic onset. Duration of vaccine protection is 10 years. Figure does not include the result for ![Graphic][66], as this efficacy has no impact on probability of occurrence of a major outbreak. ![Figure S11.](http://medrxiv.org/https://www.medrxiv.org/content/medrxiv/early/2020/04/23/2020.04.19.20070805/F18.medium.gif) [Figure S11.](http://medrxiv.org/content/early/2020/04/23/2020.04.19.20070805/F18) Figure S11. Uncertainty analysis. Model predictions for the mean cumulative number of new infections and associated 95% uncertainty interval (UI) at various levels of vaccine efficacy in reducing susceptibility (*VE**S*) generated through 500 simulation runs. Scenario assumes vaccine scale-up to 80% coverage before epidemic onset. Duration of vaccine protection is 10 years. The solid black line, dashed lines, and shades show respectively, the mean, 95% uncertainty interval, and individual estimates for the cumulative number of new infections across all 500 uncertainty runs. ![Figure S12.](http://medrxiv.org/https://www.medrxiv.org/content/medrxiv/early/2020/04/23/2020.04.19.20070805/F19.medium.gif) [Figure S12.](http://medrxiv.org/content/early/2020/04/23/2020.04.19.20070805/F19) Figure S12. Vaccine effectiveness of age-group prioritization and the reproduction number *R*. Model predictions for SARS CoV-2 attack rate at various levels of *R*, indicating also the impact on *R* of prioritizing those 60-69 years of age for vaccination or vaccinating all age groups. The blue dashed-dotted line shows the model-predicted attack rate at various levels of *R*. The blue, red, and green stars show, respectively, the model-predicted attack rates in absence of vaccination, by prioritizing vaccination at 80% coverage for those 60-69 years of age, and by extending vaccination at 80% coverage to all age groups. The figure highlights how effectiveness of the vaccine (number of vaccinations needed to avert one infection) depends on the position on the *R* curve—prioritizing vaccination for any single age group, regardless of that age group, has overall lower effectiveness than extending vaccination to all age groups. The reason is that vaccinating one age group reduces *R* only marginally, whereas vaccinating all age groups reduces *R* to an epidemic domain where small reductions in *R* can have more substantial impact on epidemic size. ![Figure S13.](http://medrxiv.org/https://www.medrxiv.org/content/medrxiv/early/2020/04/23/2020.04.19.20070805/F20.medium.gif) [Figure S13.](http://medrxiv.org/content/early/2020/04/23/2020.04.19.20070805/F20) Figure S13. Impact of varying levels of vaccine efficacy in reducing susceptibility, *VE**S*, on the cumulative number of new SARS-CoV-2 infections when the reproduction number *R* is 3. Scenario assumes vaccine scale-up to 80% coverage before epidemic onset. Duration of vaccine protection is 10 years. ## Acknowledgements This publication was made possible by NPRP grant number 9-040-3-008 and NPRP grant number 12S-0216-190094 from the Qatar National Research Fund (a member of Qatar Foundation). The authors are also grateful for support provided by the Biomedical Research Program and the Biostatistics, Epidemiology, and Biomathematics Research Core, both at Weill Cornell Medicine-Qatar. GM acknowledges support by UK Research and Innovation as part of the Global Challenges Research Fund, grant number ES/P010873/1. The statements made herein are solely the responsibility of the authors. ## Footnotes * Funding: The Qatar National Research Fund (NPRP 9-040-3-008 and NPRP 12S-0216-190094) and UK Research and Innovation as part of the Global Challenges Research Fund, grant number ES/P010873/1 * Disclose funding received for this work: others * Financial disclosure: This publication was made possible by NPRP grant number 9-040-3-008 and NPRP grant number 12S-0216-190094 from the Qatar National Research Fund (a member of Qatar Foundation). The authors are also grateful for support provided by the Biomedical Research Program and the Biostatistics, Epidemiology, and Biomathematics Research Core, both at Weill Cornell Medicine-Qatar. GM acknowledges support by UK Research and Innovation as part of the Global Challenges Research Fund, grant number ES/P010873/1. The statements made herein are solely the responsibility of the authors. * Received April 19, 2020. * Revision received April 19, 2020. * Accepted April 23, 2020. * © 2020, Posted by Cold Spring Harbor Laboratory The copyright holder for this pre-print is the author. All rights reserved. The material may not be redistributed, re-used or adapted without the author's permission. ## References 1. 1.Yang, Y., et al., The deadly coronaviruses: The 2003 SARS pandemic and the 2020 novel coronavirus epidemic in China. J Autoimmun, 2020: p. 102434. 2. 2.Lauer, S.A., et al., The Incubation Period of Coronavirus Disease 2019 (COVID-19) From Publicly Reported Confirmed Cases: Estimation and Application. Ann Intern Med, 2020. 3. 3.Wu, Z. and J.M. McGoogan, Characteristics of and Important Lessons From the Coronavirus Disease 2019 (COVID-19) Outbreak in China: Summary of a Report of 72314 Cases From the Chinese Center for Disease Control and Prevention. JAMA, 2020. 4. 4.World Health Organization, Report of the WHO-China Joint Mission on Coronavirus Disease 2019 (COVID-19). Available from :[https://www.who.int/docs/default-source/coronaviruse/who-china-joint-mission-on-covid-19-final-report.pdf](https://www.who.int/docs/default-source/coronaviruse/who-china-joint-mission-on-covid-19-final-report.pdf). Accessed on March 10, 2020. 2020. 5. 5.World Health Organization (WHO), WHO Director-General’s opening remarks at the media briefing on COVID-19 - 11 March 2020. Available from: [https://www.who.int/dg/speeches/detail/who-director-general-s-opening-remarks-at-the-media-briefing-on-covid-1911-march-2020](https://www.who.int/dg/speeches/detail/who-director-general-s-opening-remarks-at-the-media-briefing-on-covid-1911-march-2020). Accessed on March 14, 2020. 2020. 6. 6.Lu, S., Timely development of vaccines against SARS-CoV-2. Emerg Microbes Infect, 2020. 9(1): p. 542–544. [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1080/22221751.2020.1737580&link_type=DOI) 7. 7.Legido-Quigley, H., et al., Are high-performing health systems resilient against the COVID-19 epidemic? Lancet, 2020. 395(10227): p. 848–850. [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1016/S0140-6736(20)30551-1&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=http://www.n&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2020%2F04%2F23%2F2020.04.19.20070805.atom) 8. 8.Baud, D., et al., Real estimates of mortality following COVID-19 infection. Lancet Infect Dis, 2020. 9. 9.Remuzzi, A. and G. Remuzzi, COVID-19 and Italy: what next? Lancet, 2020. 10. 10.McKibbin W.J. F.R., The global macroeconomic impacts of COVID-19: seven scenarios. Available from: [http://dx.doi.org/10.2139/ssrn.3547729](http://dx.doi.org/10.2139/ssrn.3547729) 2020. 11. 11. World Health Organization (WHO), Naming the coronavirus disease (COVID-19) and the virus that causes it. Available from: [https://www.who.int/emergencies/diseases/novel-coronavirus-2019/technical-guidance/naming-the-coronavirus-disease-(covid-2019)-and-the-virus-that-causes-it](https://www.who.int/emergencies/diseases/novel-coronavirus-2019/technical-guidance/naming-the-coronavirus-disease-(covid-2019)-and-the-virus-that-causes-it). Accessed on March 11, 2020. 2020. 12. 12.National Institute of Allergy and Infectious Diseases (NIH), NIH clinical trial of investigational vaccine for COVID-19 begins. Available from: [https://www.nih.gov/news-events/news-releases/nih-clinical-trial-investigational-vaccine-covid-19-begins](https://www.nih.gov/news-events/news-releases/nih-clinical-trial-investigational-vaccine-covid-19-begins). 2020. 13. 13.Ernst D. Pipeline: Investigational Therapies for COVID-19. Available from: [https://www.infectiousdiseaseadvisor.com/home/topics/respiratory/pipeline-investigational-therapies-for-covid-19/](https://www.infectiousdiseaseadvisor.com/home/topics/respiratory/pipeline-investigational-therapies-for-covid-19/). Published on March 17, 2020. Accessed on March 18,2020. 14. 14.McLean, A.R. and S.M. Blower, Modelling HIV vaccination. Trends Microbiol, 1995. 3(12): p. 458–62. 15. 15.Blower, S.M. and A.R. McLean, Prophylactic vaccines, risk behavior change, and the probability of eradicating HIV in San Francisco. Science, 1994. 265(5177): p. 1451–4. [Abstract/FREE Full Text](http://medrxiv.org/lookup/ijlink/YTozOntzOjQ6InBhdGgiO3M6MTQ6Ii9sb29rdXAvaWpsaW5rIjtzOjU6InF1ZXJ5IjthOjQ6e3M6ODoibGlua1R5cGUiO3M6NDoiQUJTVCI7czoxMToiam91cm5hbENvZGUiO3M6Mzoic2NpIjtzOjU6InJlc2lkIjtzOjEzOiIyNjUvNTE3Ny8xNDUxIjtzOjQ6ImF0b20iO3M6NTA6Ii9tZWRyeGl2L2Vhcmx5LzIwMjAvMDQvMjMvMjAyMC4wNC4xOS4yMDA3MDgwNS5hdG9tIjt9czo4OiJmcmFnbWVudCI7czowOiIiO30=) 16. 16.McLean, A.R. and S.M. Blower, Imperfect vaccines and herd immunity to HIV. Proc Biol Sci, 1993. 253(1336): p. 9–13. [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1098/rspb.1993.0075&link_type=DOI) 17. 17.Andersson, K.M., A.D. Paltiel, and D.K. Owens, The potential impact of an HIV vaccine with rapidly waning protection on the epidemic in Southern Africa: examining the RV144 trial results. Vaccine, 2011. 29(36): p. 6107–12. 18. 18.Andersson, K.M., et al., Predicting the impact of a partially effective HIV vaccine and subsequent risk behavior change on the heterosexual HIV epidemic in low- and middle-income countries: A South African example. J Acquir Immune Defic Syndr, 2007. 46(1): p. 78–90. [PubMed](http://medrxiv.org/lookup/external-ref?access_num=17589368&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2020%2F04%2F23%2F2020.04.19.20070805.atom) 19. 19.Davenport, M.P., et al., Predicting the impact of a nonsterilizing vaccine against human immunodeficiency virus. J Virol, 2004. 78(20): p. 11340–51. [Abstract/FREE Full Text](http://medrxiv.org/lookup/ijlink/YTozOntzOjQ6InBhdGgiO3M6MTQ6Ii9sb29rdXAvaWpsaW5rIjtzOjU6InF1ZXJ5IjthOjQ6e3M6ODoibGlua1R5cGUiO3M6NDoiQUJTVCI7czoxMToiam91cm5hbENvZGUiO3M6MzoianZpIjtzOjU6InJlc2lkIjtzOjExOiI3OC8yMC8xMTM0MCI7czo0OiJhdG9tIjtzOjUwOiIvbWVkcnhpdi9lYXJseS8yMDIwLzA0LzIzLzIwMjAuMDQuMTkuMjAwNzA4MDUuYXRvbSI7fXM6ODoiZnJhZ21lbnQiO3M6MDoiIjt9) 20. 20.Wells, C., et al., Harnessing case isolation and ring vaccination to control Ebola. PLoS Negl Trop Dis, 2015. 9(5): p. e0003794. [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1371/journal.pntd.0003794&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=26024528&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2020%2F04%2F23%2F2020.04.19.20070805.atom) 21. 21.1. W. Morrow Boily, M.-C., et al., The role of mathematical models in vaccine development and public health decision making, W. Morrow, et al., Editors. 2012, Wiley-Blackwell. p. 480–508. 22. 22.Abu-Raddad, L.J., et al., Analytic insights into the population level impact of imperfect prophylactic HIV vaccines. J Acquir Immune Defic Syndr, 2007. 45(4): p. 454–67. [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1097/QAI.0b013e3180959a94&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=17554215&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2020%2F04%2F23%2F2020.04.19.20070805.atom) [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000248028100012&link_type=ISI) 23. 23.Boily, M.C., et al., Measuring the public-health impact of candidate HIV vaccines as part of the licensing process. Lancet Infect Dis, 2008. 8(3): p. 200–7. [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1016/S1473-3099(07)70292-X&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=18291341&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2020%2F04%2F23%2F2020.04.19.20070805.atom) [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000253544300022&link_type=ISI) 24. 24.Alsallaq, R.A., et al., Population level impact of an imperfect prophylactic vaccine for herpes simplex virus-2. Sex Transm Dis, 2010. 37(5): p. 290–7. [PubMed](http://medrxiv.org/lookup/external-ref?access_num=20351622&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2020%2F04%2F23%2F2020.04.19.20070805.atom) 25. 25.Gay, N.J., et al., Interpretation of serological surveillance data for measles using mathematical models: implications for vaccine strategy. Epidemiol Infect, 1995. 115(1): p. 139–56. [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1017/S0950268800058209&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=7641827&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2020%2F04%2F23%2F2020.04.19.20070805.atom) 26. 26.Michael, E., M.N. Malecela-Lazaro, and J.W. Kazura, Epidemiological modelling for monitoring and evaluation of lymphatic filariasis control. Adv Parasitol, 2007. 65: p. 191–237. [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1016/S0065-308X(07)65003-9&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=18063097&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2020%2F04%2F23%2F2020.04.19.20070805.atom) 27. 27.Basanez, M.G., et al., A research agenda for helminth diseases of humans: modelling for control and elimination. PLoS Negl Trop Dis, 2012. 6(4): p. e1548. [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1371/journal.pntd.0001548&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=22545162&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2020%2F04%2F23%2F2020.04.19.20070805.atom) 28. 28.Hill, E.M., et al., Seasonal influenza: Modelling approaches to capture immunity propagation. PLoS Comput Biol, 2019. 15(10): p. e1007096. [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1371/journal.pcbi.1007096&link_type=DOI) 29. 29.Sah, P., et al., Future epidemiological and economic impacts of universal influenza vaccines. Proc Natl Acad Sci U S A, 2019. 116(41): p. 20786–20792. [Abstract/FREE Full Text](http://medrxiv.org/lookup/ijlink/YTozOntzOjQ6InBhdGgiO3M6MTQ6Ii9sb29rdXAvaWpsaW5rIjtzOjU6InF1ZXJ5IjthOjQ6e3M6ODoibGlua1R5cGUiO3M6NDoiQUJTVCI7czoxMToiam91cm5hbENvZGUiO3M6NDoicG5hcyI7czo1OiJyZXNpZCI7czoxMjoiMTE2LzQxLzIwNzg2IjtzOjQ6ImF0b20iO3M6NTA6Ii9tZWRyeGl2L2Vhcmx5LzIwMjAvMDQvMjMvMjAyMC4wNC4xOS4yMDA3MDgwNS5hdG9tIjt9czo4OiJmcmFnbWVudCI7czowOiIiO30=) 30. 30.Gottlieb, S.L., et al., Advancing vaccine development for gonorrhoea and the Global STI Vaccine Roadmap. Sex Health, 2019. 16(5): p. 426–432. 31. 31.Gottlieb, S.L., et al., Modelling efforts needed to advance herpes simplex virus (HSV) vaccine development: Key findings from the World Health Organization Consultation on HSV Vaccine Impact Modelling. Vaccine, 2019. 37(50): p. 7336–7345. 32. 32.Peasah, S.K., et al., Cost-effectiveness of increased influenza vaccination uptake against readmissions of major adverse cardiac events in the US. PLoS One, 2019. 14(4): p. e0213499. 33. 33.Spicknall, I.H., et al., Review of mathematical models of HSV-2 vaccination: Implications for vaccine development. Vaccine, 2019. 37(50): p. 7396–7407. 34. 34.Ayoub, H.H., et al., Characterizing key attributes of the epidemiology of COVID-19 in China: Model-based estimations. medRxiv, 2020: p. 2020.04.08.20058214. 35. 35.MATLAB®, The Language of Technical Computing. The MathWorks, Inc. 2019. 36. 36.Guan, W.J., et al., Clinical Characteristics of Coronavirus Disease 2019 in China. N Engl J Med, 2020. 37. 37.Huang, C., et al., Clinical features of patients infected with 2019 novel coronavirus in Wuhan, China. Lancet, 2020. 395(10223): p. 497–506. [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1016/S0140-6736(20)30183-5&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=http://www.n&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2020%2F04%2F23%2F2020.04.19.20070805.atom) 38. 38.Novel Coronavirus Pneumonia Emergency Response Epidemiology, T., [The epidemiological characteristics of an outbreak of 2019 novel coronavirus diseases (COVID-19) in China]. Zhonghua Liu Xing Bing Xue Za Zhi, 2020. 41(2): p. 145–151. [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.3760/cma.j.issn.0254-6450.2020.02.003&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=32064853&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2020%2F04%2F23%2F2020.04.19.20070805.atom) 39. 39.United Nations Department of Economic and Social Affairs Population Dynamics, The 2019 Revision of World Population Prospects. Available from [https://population.un.org/wpp/](https://population.un.org/wpp/). Accessed on March 1st, 2020. 2020. 40. 40.COVID-19 Outbreak Live Update. Available from: [https://www.worldometers.info/coronavirus/](https://www.worldometers.info/coronavirus/). Accessed on March 14, 2020. 2020. 41. 41.Makhoul M., A.H.H., Chemaitelly H., Seedat S., Mumtaz G.R., and Abu-Raddad L.J.,, Epidemiological impact of SARS-CoV-2 vaccination: mathematical modeling analyses. submitted for publication. 42. 42.Mckay, M.D., R.J. Beckman, and W.J. Conover, A Comparison of Three Methods for Selecting Values of Input Variables in the Analysis of Output from a Computer Code. Technometrics, 1979. 21(2): p. 239–245. [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.2307/1268522&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=WOS:A1979GW1&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2020%2F04%2F23%2F2020.04.19.20070805.atom) [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=A1979GW16700012&link_type=ISI) 43. 43.Sanchez, M.A. and S.M. Blower, Uncertainty and sensitivity analysis of the basic reproductive rate - Tuberculosis as an example. American Journal of Epidemiology, 1997. 145(12): p. 1127–1137. [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1093/oxfordjournals.aje.a009076&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=9199543&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2020%2F04%2F23%2F2020.04.19.20070805.atom) [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=A1997XG00800008&link_type=ISI) 44. 44.Davies, N.G., et al., Age-dependent effects in the transmission and control of COVID-19 epidemics. medRxiv, 2020: p. 2020.03.24.20043018. 45. 45.Zhu, Y., et al., Children are unlikely to have been the primary source of household SARS-CoV-2 infections. medRxiv, 2020: p. 2020.03.26.20044826. 46. 46.Li, R., et al., Substantial undocumented infection facilitates the rapid dissemination of novel coronavirus (SARS-CoV2). Science, 2020. 47. 47.Verity R., et al., Estimates of the severity of coronavirus disease 2019: a model-based analysis. The Lancet, 2020. 20: p. 30243–7. ## References 1. 1.Li, R., et al., Substantial undocumented infection facilitates the rapid dissemination of novel coronavirus (SARS-CoV2). Science, 2020. 2. 2.Lauer, S.A., et al., The Incubation Period of Coronavirus Disease 2019 (COVID-19) From Publicly Reported Confirmed Cases: Estimation and Application. Ann Intern Med, 2020. 3. 3.Zou, L., et al., SARS-CoV-2 Viral Load in Upper Respiratory Specimens of Infected Patients. N Engl J Med, 2020. 4. 4.Rothe, C., et al., Transmission of 2019-nCoV Infection from an Asymptomatic Contact in Germany. N Engl J Med, 2020. 382(10): p. 970–971. [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1056/NEJMc2001468&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=http://www.n&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2020%2F04%2F23%2F2020.04.19.20070805.atom) 5. 5.World Health Organization, Report of the WHO-China Joint Mission on Coronavirus Disease 2019 (COVID-19). Available from :[https://www.who.int/docs/default-source/coronaviruse/who-china-joint-mission-on-covid-19-final-report.pdf](https://www.who.int/docs/default-source/coronaviruse/who-china-joint-mission-on-covid-19-final-report.pdf). Accessed on March 10, 2020. 2020. 6. 6.United Nations Department of Economic and Social Affairs Population Dynamics, The 2019 Revision of World Population Prospects. Available from [https://population.un.org/wpp/](https://population.un.org/wpp/). Accessed on March 1st, 2020. 2020. 7. 7.Novel Coronavirus Pneumonia Emergency Response Epidemiology, T., [The epidemiological characteristics of an outbreak of 2019 novel coronavirus diseases (COVID-19) in China]. Zhonghua Liu Xing Bing Xue Za Zhi, 2020. 41(2): p. 145–151. [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.3760/cma.j.issn.0254-6450.2020.02.003&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=32064853&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2020%2F04%2F23%2F2020.04.19.20070805.atom) 8. 8.Wu, Z. and J.M. McGoogan, Characteristics of and Important Lessons From the Coronavirus Disease 2019 (COVID-19) Outbreak in China: Summary of a Report of 72314 Cases From the Chinese Center for Disease Control and Prevention. JAMA, 2020. 9. 9.Guan, W.J., et al., Clinical Characteristics of Coronavirus Disease 2019 in China. N Engl J Med, 2020. 10. 10.Huang, C., et al., Clinical features of patients infected with 2019 novel coronavirus in Wuhan, China. Lancet, 2020. 395(10223): p. 497–506. [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1016/S0140-6736(20)30183-5&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=http://www.n&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2020%2F04%2F23%2F2020.04.19.20070805.atom) 11. 11.Ayoub, H.H., et al., Characterizing key attributes of the epidemiology of COVID-19 in China: Model-based estimations. medRxiv, 2020: p. 2020.04.08.20058214. 12. 12.Heffernan, J.M., R.J. Smith, and L.M. Wahl, Perspectives on the basic reproductive ratio. Journal of the Royal Society Interface, 2005. 2(4): p. 281–293. 13. 13.Whittle, P., The outcome of a stochastic epidemic—a note on Bailey’s paper. Biometrika, 1955. 42(1-2): p. 116–122. [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1093/biomet/42.1-2.116&link_type=DOI) [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=A1955XF69000009&link_type=ISI) [1]: /embed/inline-graphic-1.gif [2]: /embed/inline-graphic-2.gif [3]: /embed/inline-graphic-3.gif [4]: /embed/inline-graphic-4.gif [5]: /embed/inline-graphic-5.gif [6]: /embed/inline-graphic-6.gif [7]: /embed/inline-graphic-7.gif [8]: /embed/inline-graphic-8.gif [9]: /embed/inline-graphic-9.gif [10]: F1/embed/inline-graphic-10.gif [11]: F2/embed/inline-graphic-11.gif [12]: /embed/inline-graphic-12.gif [13]: /embed/inline-graphic-13.gif [14]: /embed/inline-graphic-14.gif [15]: /embed/inline-graphic-15.gif [16]: /embed/inline-graphic-16.gif [17]: /embed/inline-graphic-17.gif [18]: /embed/inline-graphic-18.gif [19]: F3/embed/inline-graphic-19.gif [20]: F3/embed/inline-graphic-20.gif [21]: /embed/inline-graphic-21.gif [22]: /embed/inline-graphic-22.gif [23]: F7/embed/inline-graphic-23.gif [24]: F7/embed/inline-graphic-24.gif [25]: F7/embed/inline-graphic-25.gif [26]: /embed/inline-graphic-26.gif [27]: /embed/inline-graphic-27.gif [28]: /embed/graphic-10.gif [29]: /embed/graphic-11.gif [30]: /embed/graphic-12.gif [31]: /embed/graphic-13.gif [32]: /embed/graphic-16.gif [33]: /embed/graphic-17.gif [34]: /embed/graphic-18.gif [35]: /embed/graphic-21.gif [36]: /embed/graphic-22.gif [37]: /embed/inline-graphic-28.gif [38]: /embed/inline-graphic-29.gif [39]: /embed/graphic-23.gif [40]: /embed/inline-graphic-30.gif [41]: /embed/inline-graphic-31.gif [42]: /embed/inline-graphic-32.gif [43]: /embed/inline-graphic-33.gif [44]: /embed/inline-graphic-34.gif [45]: /embed/inline-graphic-35.gif [46]: /embed/inline-graphic-36.gif [47]: /embed/inline-graphic-37.gif [48]: /embed/inline-graphic-38.gif [49]: /embed/inline-graphic-39.gif [50]: /embed/inline-graphic-40.gif [51]: /embed/inline-graphic-41.gif [52]: /embed/inline-graphic-42.gif [53]: /embed/inline-graphic-43.gif [54]: /embed/inline-graphic-44.gif [55]: /embed/inline-graphic-45.gif [56]: F9/embed/inline-graphic-46.gif [57]: F10/embed/inline-graphic-47.gif [58]: F11/embed/inline-graphic-48.gif [59]: F12/embed/inline-graphic-49.gif [60]: F13/embed/inline-graphic-50.gif [61]: F13/embed/inline-graphic-51.gif [62]: F14/embed/inline-graphic-52.gif [63]: F14/embed/inline-graphic-53.gif [64]: F17/embed/inline-graphic-54.gif [65]: F17/embed/inline-graphic-55.gif [66]: F17/embed/inline-graphic-56.gif