Incorporating the mutational landscape of SARS-COV-2 variants and case-dependent vaccination rates into epidemic models ======================================================================================================================= * Mohammad Mihrab Chowdhury * Md Rafiul Islam * Md Sakhawat Hossain * Nusrat Tabassum * Angela Peace ## Abstract Coronavirus Disease (COVID-19), which began as a small outbreak in Wuhan, China in December 2019, became a global pandemic within months due to its high transmissibility. In the absence of pharmaceutical treatment, various non-pharmaceutical interventions (NPIs) to contain the spread of COVID-19 brought the entire world to a halt. After almost a year of seemingly returning to normalcy with the world’s quickest vaccine development, the emergence of more infectious and vaccine resistant coronavirus variants is bringing the situation back to where it was a year ago. In the light of this new situation, we conducted a study to portray the possible scenarios based on the three key factors : impact of interventions (pharmaceutical and NPIs), vaccination rate, and vaccine efficacy. In our study, we assessed two of the most crucial factors, transmissibility and vaccination rate, in order to reduce the spreading of COVID-19 in a simple but effective manner. In order to incorporate the time-varying mutational landscape of COVID-19 variants, we estimated a weighted transmissibility composed of the proportion of existing strains that naturally vary over time. Additionally, we consider time varying vaccination rates based on the number of daily new cases. Our method for calculating the vaccination rate from past active cases is an effective approach in forecasting probable future scenarios as it actively tracks people’s attitudes toward immunization as active cases change. Our simulations show that if a large number of individuals cannot be vaccinated in a short period of time, adopting NPIs is the best approach to manage disease transmission with the emergence of new vaccine breakthrough variants and more infectious variants. Keywords * Break-Through Variants * Waning Immunity * Reinfection * Non-Pharmaceutical Interventions (NPIs) * Vaccination Rate ## 1. Introduction The catastrophic effect of the H1N1 virus (Influenza Flu) in the twentieth century has returned with the worldwide spread of coronavirus. Coronavirus illness (COVID-19) is an excellent illustration of what may happen if the proper preventive measures are not taken in a timely way, even in this advanced age of technology and medical understanding (Rossen et al., 2020). Although NPIs have played a critical role in reducing the death toll and number of infected people in the absence of medical treatment, delayed interventions and underestimating the severity of the situation have resulted in over 205 million infected cases, 4.3 million deaths, and billions of dollars in economic losses worldwide (Cutler and Summers, 2020). Due to delayed measures during the initial wave, coronavirus is responsible for between 593,126 and 752,284 excess deaths in the USA (Center for Disease Control, 2021b). To counteract the economic losses and death toll, the world’s quickest vaccine was invented (Andreadakis et al., 2020). With the proper allocation of the vaccine, the number of infected persons and disease-related mortality has reduced over time, allowing people to return to their normal lifestyles without using NPIs after more than a year (Islam et al., 2021). However, the advent of new variants (delta, lambda, delta plus, omicron and so on) is prompting a new pandemic wave to sweep the world right now. With the mass vaccination program, there was an almost exponential drop in the number of infected persons between January 8, 2021 to June 7, 2021 but with the emergence of new variants, the number of infected people is exponentially growing again. If proper steps are not implemented quickly, the new variants could produce a considerably worse situation than previous coronavirus waves. The false sense of security can reenact the historic super-spreader incident of the Philadelphia Parade-1918, which resulted in 13000 more deaths in Philadelphia in three months owing to the Spanish flu pandemic (Stetler, 2017). In light of this situation, in our study, we sought to model scenarios by varying the overall population adoption of NPIs in conjunction with vaccine efficacy, antibody waning, and change in prevalence of the co-circulating variants of the coronavirus over time. One of the most pressing questions during the COVID-19 pandemic was how to correctly parameterize the value of transmission rate (*β*). It has been quite challenging to keep up with the evolution of coronavirus in estimating *β* value. Due to the continual appearance of newer strains of coronavirus with distinct characteristics, estimating a transmissibility value without taking the strains into consideration has now become redundant. In our model, we have considered a time-varying variant dependent transmission parameter (*β*) by adapting the approach from Islam et al. (2021). To assess the transmissibility, our approach takes into consideration a crucial factor: the changing nature of coronavirus predominant strains over time, making it useful in the present. On the other hand, boosting the number of vaccinated people is one of the most important elements in reducing COVID-19’s devastating effects (Mahmud et al., 2021). As people’s perceptions fluctuate over time, especially with the number of active cases (Center for Disease Control, 2021a), fixed vaccination rate will give a false hope to curb the disease. This method gives a fair picture of how the vaccination rate changes when the number of current cases changes, which is one of the key reasons why individuals become vaccinated (Hamel et al., 2021). Furthermore, with the emergence of new variants, this strategy can be very valuable for planning the third dose and booster dose vaccine based just on historical data of daily cases. As vaccine protection lasts at least six months (Pfizer Inc, 2021), we observed from our simulations that the best way to reduce disease-induced death is to increase vaccination rates, presuming the vaccine is extremely effective against all variants. However, owing to people’s concern over vaccination safety, NPIs are the only thing that can prevent unvaccinated individuals from falling ill. Interestingly, we found that NPIs had an effect comparable to raising vaccination rates to control the epidemic. As a result, the best method for reducing the spread and death toll is to boost immunization and implement NPIs until herd immunity is achieved. With more data available, this model can provide more accurate projections. ## 2. Methodology ### 2.1. Model development We developed an extension of the basic SEIR type compartmental deterministic model to project future disease dynamics in the USA. Our model consists of two distinct tracks for vaccinated and unvaccinated persons, which take into consideration vaccination rates, waning immunity, viral variantions, and reinfection. The compartmental model diagram is shown in figure 1. The model divides the population into seven compartments: unvaccinated susceptible individuals (*S*), vaccinated susceptible individuals (*V*), unvaccinated exposed individuals (*E*), vaccinated exposed individuals (*E**v*), unvaccinated infected individuals (*I*), vaccinated infected individuals (*I**v*), and recovered individuals (*R*). ![Figure 1:](http://medrxiv.org/https://www.medrxiv.org/content/medrxiv/early/2021/12/19/2021.11.28.21266882/F1.medium.gif) [Figure 1:](http://medrxiv.org/content/early/2021/12/19/2021.11.28.21266882/F1) Figure 1: SEIR (Susceptible, Exposed, Infected, Recovered) type Compartmental Model. Here people moves from susceptible (*S*) and vaccinated (*V*) compartment to exposed compartment (*E* and *E**v*) at *βσ* and *βρσ* respectively. From *E* and *E**v* people become infected (*I* and *I**v*) after incubation period *α*. Then from *I* and *I**v*, infected people either go to the recovery compartment *R* or die. With the waning of antibodies, individuals move back to the *S* compartment from the *R* and *V* compartments. Here, we assume that successful transmission has occurred for both exposed vaccinated and unvaccinated individuals (*E* and *E**v*). Unlike existing models Byambasuren et al. (2020); Buitrago-Garcia et al. (2020); Aguilar et al. (2020); Kronbichler et al. (2020), this model does not explicitly track asymptomatic infections separately from symptomatic infections. This will not impact the result of our model as the asymptomatic cases go to the recovery without getting identified. In fact, there isn’t enough information on asymptomatic cases of the variants. Further to that, due to the increased virulence of the new variants, we focused our research on how the changing proportions of variants affects symptomatic people, for whom the cases may be more deadly than the original strain. ### 2.2. Disease Dynamics The susceptible population is exposed to the illness at a rate of *β* and vaccinated at a rate of *η*. Although, the vaccinated susceptible population is exposed to the disease at the same rate as the unvaccinated people, but only a fraction of them will become infected due to the vaccine’s protection. After the incubation period (*α*), the exposed population goes to *I* and *I**v* compartment. Then infected persons transfer from *I* and *I**v* compartments by recovering at *γ* and *γ**v* or by passing away at *μ* and *μ**v* rates respectively. We are assuming both vaccinated and non-vaccinated people will take same time to recover. The vaccinated and recovered population again become susceptible owing to the decrease in immunity over time at a rates *ω**v* and *ω* respectively. The following system of differential equations is our model equations which keep tracks how the individual moves from one compartment to another compartment. ![Formula][1] ### 2.3. Co-circulating Variants After nearly a year of pandemic, the most crucial thing is to track the strains precisely as the fate of the pandemic is driven by these emerging and circulating strains. The real-time monitoring of the prevalent strain among co-circulating strains is fundamental in implementing optimal pharmacological and non-pharmaceutical treatments as the characteristics of the pandemic will change with the change in prevalent strain. Moreover, this will help to determine more infectious strains, viral load, and fatal variants. The US Department of Health and Human Services (HHS) established a SARS-CoV-2 Interagency Group (SIG) that included the CDC, National Institutes of Health (NIH), Food and Drug Administration (FDA), Biomedical Advanced Research and Development Authority (BARDA), and the Department of Defense (DoD) in order to keep track of co-circulating variants. The SIG classifies the circulating variants into four categories: VBM (Variants Being Monitored), VOC (Variants of Concern), VOI (Variant of Interest), and VOHC (Variant of High Concern). This classification of variants is based on the proportion of co-circulating variants at the national and regional levels, as well as the impact of mutations on intervention strategies, the severity of illness, and disease transmissibility. Ten variants are currently classified as VBM, two as VOI, and none as VOC or VOHC (Centers for Disease Control and Prevention, 2021b). The table 1 contains the co-circulating variants classified as VBM and VOC in the United States, along with their relative transmissibility. View this table: [Table 1:](http://medrxiv.org/content/early/2021/12/19/2021.11.28.21266882/T1) Table 1: The SARS-CoV-2 co-circulating variants in the United States that were classified as VBM and VOC by SIG and their relative transmissibility (as of Decemeber 16, 2021) View this table: [Table 2:](http://medrxiv.org/content/early/2021/12/19/2021.11.28.21266882/T2) Table 2: Model Parameters and their description ### 2.4. Model Parameterization Our model consists of ten parameters. Depending on the availability of data, we estimated some and took some of the well established parameters value based on literature. Determining the value of transmissibility, *β*, is the most challenging parameter of all due to the coronavirus’s ever-changing nature and the abundance of circulating variants. To take into account circulating variants and their relative transmissibility, we used the data from Center for Disease Control (2021a) to generate a weighted average of transmissibility based on relative infectivity and circulating variants prevalence over time by incorporating the approach from Islam et al. (2021). Then *β* becomes time dependent that takes into account the emergence of newer virulent variants. By plotting the daily infection and vaccination data over time from the Center for Disease Control (2021a), we see that the number of people who have been vaccinated per day fluctuates according to the number of infected people per day(figure 3a). We assume that the time-varying vaccination rate (*η*) is proportional to the number of active cases based on this phenomenon. If the number of cases rises, more individuals will be interested to take the vaccine. Using the data from Center for Disease Control (2021a), we have calculated the vaccination rate (*η*) by using the formula, ![Graphic][2], where for total population we used the most recent value from the U.S. Census Bureau (2020). ![Figure 2:](http://medrxiv.org/https://www.medrxiv.org/content/medrxiv/early/2021/12/19/2021.11.28.21266882/F2.medium.gif) [Figure 2:](http://medrxiv.org/content/early/2021/12/19/2021.11.28.21266882/F2) Figure 2: (a) Based on CDC’s Prevalence statistics, the change in prevalence of variants of coronavirus over time from May 29 to Aug 28 (Center for Disease Control, 2021a). (b) The relative infectivity of circulating viral variants is shown based on 50% increase in infectivity for alpha (B.1.1.7)with respect to gamma and other co-circulating variants and a 50% increase in infectivity for delta (B.1.617.2) in comparison to alpha (B.1.1.7) (blue circles). To estimate future relative infectivity, a logistic equation is fitted with the model (pink curve). (See the Appendix C for fitted logistic equation) ![Figure 3:](http://medrxiv.org/https://www.medrxiv.org/content/medrxiv/early/2021/12/19/2021.11.28.21266882/F3.medium.gif) [Figure 3:](http://medrxiv.org/content/early/2021/12/19/2021.11.28.21266882/F3) Figure 3: (a) Daily data on infections and vaccination rate from the Center for Disease Control (2021a). (b) The red solid line here represents the linear function which estimates the vaccination rate. Using the vaccination data from the Center for Disease Control (2021a), we utilized statistical techniques to fit the vaccination rate to the number of daily cases shown in figure 3b. To estimate the vaccination rate using the linear regression approach, we used the number of daily cases as a predictor to estimate the vaccination rate using the linear regression approach. This fitted function will predict future vaccination rates based on the values of past cases. We have assumed the incubation period (*α*) is the same for both vaccinated and unvaccinated people to become infected (Center for Disease Control, 2021c). The recovery rate (*γ* and *γ**v*) and disease induced death rate (*μ*) for the non vaccinated persons are adopted from Usherwood et al. (2021). We have used disease induced death rate (*μ**v*) adopted from Vahidy et al. (2021) for vaccinated people. Based on the recent study, we assumed that the vaccination would be effective for at least 6 months (Pfizer Inc, 2021) and the recovered patients would have the same period of protection as vaccinated persons. ## 3. Results ### 3.1. Analytic Results The basic reproduction number (*R*) is an epidemiologic statistic used to describe the transmissibility of infectious agents. Although the vaccination of susceptible individuals of the community will limit the number of effective contacts between infectious and susceptible people, this activity will not reduce the *R* value as *R* assumes a completely susceptible population (Delamater et al., 2019). As a result, in addition to calculating the basic reproduction number (*R*) without the vaccination compartment, we have also estimated the effective reproduction number (*R**e*). The analytic expression of the effective reproduction number (*R**e*) is, ![Formula][3] The detailed calculation using the next generation matrix approach of the basic reproduction number and effective reproduction number is shown in the Appendix A and B. The *R**e* value, rather than the *R* value, would be reduced if the number of susceptible people in a population were reduced through vaccination. If *R**e* can be decreased to a value *<* 1, immunization might potentially halt an epidemic (Anderson and May, 1992; Anderson, 1992; Rubió, 2012). ### 3.2. Co-Circulating Variants, Vaccine Efficacy and NPIs Keeping a steady rate of viral transmission, we simulated the scenarios based on the adoption of NPIs and vaccine effectiveness (figure: 4). The reduction in vaccine efficacy occurs with time as well as with the advent of new vaccine breakthrough variants. People, on the other hand, are adopting NPIs as they see fit because there is no set of mandates for individuals to follow when it comes to NPI strategy, making it impossible to get information on which NPI approach is employed by how many people. Instead of concentrating on which NPIs are used, it is now more appropriate to take into account the proportion of the population that uses NPIs as a measure of prevention. ![Figure 4:](http://medrxiv.org/https://www.medrxiv.org/content/medrxiv/early/2021/12/19/2021.11.28.21266882/F4.medium.gif) [Figure 4:](http://medrxiv.org/content/early/2021/12/19/2021.11.28.21266882/F4) Figure 4: The effect of NPIs adoption and vaccine efficacy on the daily number of infected cases is depicted in this graph. Figures a, b, and c depict the change in daily case count with fixed vaccination rate (70%,80%,95%) and changing adoption of NPIs and vaccine effectiveness, while figures d, e, and f depict the corresponding cumulative case count.From a,b,c,d,e,f, it is clear that if the majority of the population uses NPIs, NPIs will have the greatest influence on lowering the daily number of instances in the short and long term. If NPI adoption is not achievable, the only way to minimize the number of cases is to have a high vaccination rate with high effectiveness. (Check out Appendix D for a comparison of Peak values.) To simulate the scenarios of waning efficacy of vaccine, we varied the effectiveness of the vaccine from 70% to 95% (figure 4a, 4b, 4c). In the simulation, as the efficacy rises, the number of persons afflicted decreases significantly, implying that a highly effective vaccine is needed to reduce the number of illnesses in the short and long term. On the other hand, if the number of infected people decreases, people will be less likely to use NPIs which is why we have varied the total populations adoption of NPIs to simulate the situations with different adoption rate. Figure 4 indicates that even if 80% of the population uses NPIs, the number of sick people will remain relatively low, regardless of vaccine effectiveness. It implies that using NPIs, until herd immunity is obtained or high vaccine efficacy is maintained indefinitely, is the best strategy. ### 3.3. Effective Reproduction Number The effective reproduction number *R**e* is affected by changes in immunization rate, vaccine efficacy, and NPI adoption. While vaccine efficacy cannot be modified quickly, NPI adoption and vaccination rate can change drastically. Also, against the new variants, vaccine loses its efficacy (Center for Disease Control, 2021d) for whcih we calculated the *R**e* value by setting the vaccine effectiveness to 80% instead of 90% while varying the vaccination rate, *η* from .0001 to .003, as well as the NPIs adoption rate from 0% to 80%. The effective reproduction number fluctuates between 3.1953 (best scenario) and 0.4553 (worst scenario). ![Figure 5:](http://medrxiv.org/https://www.medrxiv.org/content/medrxiv/early/2021/12/19/2021.11.28.21266882/F5.medium.gif) [Figure 5:](http://medrxiv.org/content/early/2021/12/19/2021.11.28.21266882/F5) Figure 5: Heat-map here depict how with fixed vaccine efficiency the effective reproduction number (*R**e*), varies with vaccine rates and adoption of NPIs. With vaccine efficacy 80%, the most efficient method to halt the spread of the disease is to increase immunization rates paired with good implementation of NPIs until herd immunity is reached. ## 4. Discussion The adoption of NPIs, effective allocation, and deployment by assuring the high efficacy of the vaccine against the virus throughout time are the most significant components in controlling the spread of the COVID-19. With the introduction of novel coronavirus variants and declining immunity over time, along with an increasing perception of safety, has produced an ideal environment for the superspreading of highly infectious and lethal coronavirus variants. In this situation, we developed an extension of the standard compartmental SEIR type model that dynamically tracks changes in co-circulating variants and dynamic vaccination rates. As time passes, the emergence of new coronavirus variants renders a fixed transmission rate obsolete, because each variant differs in terms of transmission potential. The overall coronavirus transmission rate changes over time depending on the transmissibility of the co-circulating variants. As a result, tracking the change in *β* value with the change in co-circulating variants over time is a useful technique to track the transmission rate of coronavirus at any given time. Furthermore, the vaccination data was not uniform at first due to people’s anxiety and misinformation about the COVID-19 vaccinated. Based on our data fitting, it is clear that the vaccination rate currently follows the trend in daily number of cases with a lag, allowing us to create a time-varying case-dependent vaccination rate. This strategy particularly accounts for people’s shifting attitudes toward vaccination as the number of cases grows, making it a practical way to account for variants and vaccination rates in the current situation. To reduce the complexity caused by the continual change in variants, we did not explicitly include asymptomatic infected people in our model, instead focusing on the shifting landscape of variant strains and parameterizing relative transmissibility using empirical time series weighted prevalence data for symptomatic people (figure 2). Furthermore, rather than indicating which NPIs are most beneficial, we now state what percentage of the population has implemented NPIs. This is the logical thing to do right now, because individuals are no longer adopting a predetermined set of NPIs plan, but rather choosing to protect themselves against the COVID-19 in their own way, as states no longer have a mandate. Moreover, to avoid the ambiguity about the efficacy of different NPIs measures, we concentrated on the overall adoption of NPIs by the proportion of the population. Based on the findings of our model, we can definitely claim that vaccination rate, vaccine effectiveness, and NPIs are significant factors in reducing the number of infected patients and controlling the disease. According to our simulation, if vaccination efficacy declines slightly, the number of sick people increases considerably, demonstrating that vaccine efficacy is crucial for viral containment. As a result, the race to create and deliver booster doses and third dose is critical right now in order to maintain the vaccine’s efficacy against new strains. On the other hand, if mass vaccination is not possible in a timely way or if vaccine breakthrough variants become prevalent, NPIs are the only option for controlling the spread of COVID-19. The project’s next phase will be to add asymptomatic individuals with strain-based data. Furthermore, vaccination data will grow more consistent over time, improving the accuracy of this model. Therefore, this modeling technique highlights the importance of the individuals to become vaccinated, as well as the need of planning and distributing booster dosages in a timely manner in order to mitigate the spread of coronavirus. ## Data Availability All data produced in the present study are available upon reasonable request to the authors. ## Supplementary Materials ### A. Base Reproduction Number The Jacobian matrix J is obtained from linearizing the system 1. The disease-free equilibria of the model at *X* = (*S*,‥…0, 0); where *S* is the initial population of susceptible individuals. Evaluating the Jacobian matrix at *X* yields ![Formula][4] The linearized infected subsystem of 1 evolves according to the following set of equations for minor perturbations *z* = (*S, E, I, V, E**v*, *I**v*, *R*) near the *X* ![Formula][5] Here, ![Formula][6] We decompose the matrix M into transmission (T) and transition (Σ) matrices respectively, obtaining ![Formula][7] where ![Formula][8] and ![Formula][9] K is a large-domain next-generation matrix. Because T is ranked first, the NGM K is likewise ranked first. As a result, non-zero entries are only found in the first row of K. Therefore, K’s spectral radius is the first item on the diagonal, i.e., *K*11, which is identical to *R*. Therefore, ![Formula][10] ### B. Effective Reproduction Number The Jacobian matrix J is obtained from linearizing the system 1. The disease-free equilibria of the model at ![Graphic][11]; where *S* is the initial population of susceptible individuals. Now, reorganizing the DFE using, N = S + V as at DFE the population are either in susceptible compartment or in vaccinated compartment, we get, ![Graphic][12] Evaluating the Jacobian matrix at *X* yields ![Formula][13] The linearized infected subsystem of 1 evolves according to the following set of equations for minor perturbations *z* = (*S, E, I, V, E**v*, *I**v*, *R*) near the *X* ![Formula][14] Here, ![Formula][15] We decompose the matrix M into transmission (T) and transition (Σ) matrices respectively, obtaining ![Formula][16] where ![Formula][17] and ![Formula][18] K is a large-domain next-generation matrix. Because T is ranked first, the NGM K is likewise ranked first. As a result, non-zero entries are only found in the first row of K. Therefore, K’s spectral radius is the first item on the diagonal, i.e., *K*11, which is identical to *R**e*. Therefore, ![Formula][19] ### C. Data Fitting #### C.1. Relative Infectivity Relative infectivity is fitted against the number of days as time from May 29,2021 to August 28,2021 using the following logistic model. ![Formula][20] where as the p-values are within range, all of the parameters *a*1, *b*1,and *c*1 are statistically significant (table 3). View this table: [Table 3:](http://medrxiv.org/content/early/2021/12/19/2021.11.28.21266882/T3) Table 3: Model Parameters and their description. The fitted model: ![Formula][21] #### C.2. Vaccination Rate Based on the daily number of cases, we built a simple linear regression model to estimate the vaccination rate. The fitted model is ![Formula][22] where all the estimate coefficient of the daily instances is statistically significant as all the *p − values* are small. The fitted model: ![Formula][23] ### D. Peak Infections The number of infections at their peak fluctuates dramatically depending on the circumstance. This statistic is critical for the health sector to take adequate precautions in order to prepare for the worst-case scenario in terms of affected persons. ![Figure 6:](http://medrxiv.org/https://www.medrxiv.org/content/medrxiv/early/2021/12/19/2021.11.28.21266882/F6.medium.gif) [Figure 6:](http://medrxiv.org/content/early/2021/12/19/2021.11.28.21266882/F6) Figure 6: The variation in the peak value of daily infected cases as a function of NPI uptake and vaccination effectiveness. With an increase in NPI uptake or vaccine efficacy, the peak value is much lower. ## 5. Acknowledgement I would like to thank Dr. Linda J.S. Allen, Professor, Department of Mathematics and Statistics, Texas Tech University, for helping me with the project. * Received November 28, 2021. * Revision received December 19, 2021. * Accepted December 19, 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. J. B. Aguilar, J. S. Faust, L. M. Westafer, and J. B. Gutierrez. Investigating the impact of asymptomatic carriers on COVID-19 transmission. MedRxiv, 2020. 2. M. Alkhatib, V. Svicher, R. Salpini, F. A. Ambrosio, M. C. Bellocchi, L. Carioti, L. Piermatteo, R. Scutari, G. Costa, A. Artese, et al. Sars-cov-2 variants and their relevant mutational profiles: Update summer 2021. Microbiology spectrum, 9(3):e01096–21, 2021. 3. R. M. Anderson. The concept of herd immunity and the design of community-based immunization programmes. Vaccine, 10(13):928–935, 1992. [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1016/0264-410X(92)90327-G&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=1471414&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2021%2F12%2F19%2F2021.11.28.21266882.atom) [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=A1992JV24600010&link_type=ISI) 4. R. M. Anderson and R. M. May. Infectious diseases of humans: dynamics and control. Oxford University Press, 1992. 5. Z. Andreadakis, A. Kumar, R. G. Román, S. Tollefsen, M. Saville, and S. Mayhew. The COVID-19 vaccine development landscape. Nature reviews. Drug discovery, 19(5):305–306, 2020. [PubMed](http://medrxiv.org/lookup/external-ref?access_num=http://www.n&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2021%2F12%2F19%2F2021.11.28.21266882.atom) 6. R. K. Borchering, C. Viboud, E. Howerton, C. P. Smith, S. Truelove, M. C. Runge, N. G. Reich, L. Contamin, J. Levander, J. Salerno, et al. Modeling of future COVID-19 cases, hospitalizations, and deaths, by vaccination rates and nonpharmaceutical intervention scenarios—United States, April–September 2021. Morbidity and Mortality Weekly Report, 70(19):719, 2021. 7. D. Buitrago-Garcia, D. Egli-Gany, M. J. Counotte, S. Hossmann, H. Imeri, A. M. Ipekci, G. Salanti, and N. Low. Occurrence and transmission potential of asymptomatic and presymptomatic SARS-CoV-2 infections: A living systematic review and meta-analysis. PLoS medicine, 17(9):e1003346. 2020. 8. O. Byambasuren, M. Cardona, K. Bell, J. Clark, M.-L. McLaws, and P. Glasziou. Estimating the extent of asymptomatic COVID-19 and its potential for community transmission: systematic review and meta-analysis. Official Journal of the Association of Medical Microbiology and Infectious Disease Canada, 5(4):223–234, 2020. 9. F. Campbell, B. Archer, H. Laurenson-Schafer, Y. Jinnai, F. Konings, N. Batra, B. Pavlin, K. Vandemaele, M. D. Van Kerkhove, T. Jombart, et al. Increased transmissibility and global spread of sars-cov-2 variants of concern as at june 2021. Eurosurveillance, 26(24):2100509, 2021. 10. Center for Disease Control. Covid-data-tracker. [https://covid.cdc.gov/covid-data-tracker/#datatracker-home](https://covid.cdc.gov/covid-data-tracker/#datatracker-home) (Last accessed Sep 12, 2021), 2021a. 11. Center for Disease Control. Excess deaths associated with COVID-19. [https://www.cdc.gov/nchs/nvss/vsrr/covid19/excess\_deaths.htm](https://www.cdc.gov/nchs/nvss/vsrr/covid19/excess_deaths.htm) (Sep 12, 2021), 2021b. 12. Center for Disease Control. Interim Clinical Guidance for Management of Patients with Confirmed Coronavirus Disease COVID-19. [https://www.cdc.gov/coronavirus/2019-ncov/hcp/clinical-guidance-management-patients.html](https://www.cdc.gov/coronavirus/2019-ncov/hcp/clinical-guidance-management-patients.html) x(Feb. 16, 2021), 2021c. 13. Center for Disease Control. Science brief: COVID-19 vaccines and vaccination. [https://www.cdc.gov/coronavirus/2019-ncov/science/science-briefs/fully-vaccinated-people.html#ref15](https://www.cdc.gov/coronavirus/2019-ncov/science/science-briefs/fully-vaccinated-people.html#ref15) (July 27, 2021), 2021d. 14. Centers for Disease Control and Prevention. Prioritizing COVID-19 contact tracing mathematical modeling methods and findings. [https://www.cdc.gov/coronavirus/2019-ncov/php/contact-tracing/contact-tracing-plan/prioritization/mathematicalmodeling.html](https://www.cdc.gov/coronavirus/2019-ncov/php/contact-tracing/contact-tracing-plan/prioritization/mathematicalmodeling.html)(Nov 23, 2020), 2020. 15. Centers for Disease Control and Prevention. Delta variant: What we know about the science. [https://www.cdc.gov/coronavirus/2019-ncov/variants/delta-variant.html](https://www.cdc.gov/coronavirus/2019-ncov/variants/delta-variant.html), 2021a. 16. Centers for Disease Control and Prevention. Sars-cov-2 variant classifications and definitions. [https://www.cdc.gov/coronavirus/2019-ncov/variants/variant-classifications.html](https://www.cdc.gov/coronavirus/2019-ncov/variants/variant-classifications.html), 2021b. 17. D. M. Cutler and L. H. Summers. The COVID-19 Pandemic and the $16 Trillion Virus. JAMA, 324 (15):1495–1496, 10 2020. ISSN 0098-7484. doi: 10.1001/jama.2020.19759. URL [https://doi.org/10.1001/jama.2020.19759](https://doi.org/10.1001/jama.2020.19759). [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1001/jama.2020.19759&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=33044484&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2021%2F12%2F19%2F2021.11.28.21266882.atom) 18. N. G. Davies, S. Abbott, R. C. Barnard, C. I. Jarvis, A. J. Kucharski, J. D. Munday, C. A. Pearson, T. W. Russell, D. C. Tully, A. D. Washburne, et al. Estimated transmissibility and impact of sars-cov-2 lineage b. 1.1. 7 in england. Science, 372(6538), 2021. 19. P. L. Delamater, E. J. Street, T. F. Leslie, Y. T. Yang, and K. H. Jacobsen. Complexity of the basic reproduction number (R0). Emerging infectious diseases, 25(1):1, 2019. [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.3201/eid2501.171901&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=30560777&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2021%2F12%2F19%2F2021.11.28.21266882.atom) 20. L. Hamel, L. Lopes, G. Sparks, A. Kirzinger, A. Kearney, M. Stokes, and M. Brodie. Kff COVID-19 vaccine monitor: September 2021, 2021. 21. M. R. Islam, T. Oraby, A. McCombs, M. M. Chowdhury, M. Al-Mamun, M. G. Tyshenko, and C. Kadelka. Evaluation of the United States COVID-19 vaccine allocation strategy. PLoS ONE, 16(11):e0259700, 2021. 22. K. Katella. Comparing the COVID-19 vaccines: how are they different? Yale Medicine. Family Health. [https://www.yalemedicine.org/news/COVID-19-vaccine-comparison](https://www.yalemedicine.org/news/COVID-19-vaccine-comparison), 2021. 23. A. Kronbichler, D. Kresse, S. Yoon, K. H. Lee, M. Effenberger, and J. I. Shin. Asymptomatic patients as a source of COVID-19 infections: A systematic review and meta-analysis. International journal of infectious diseases, 98:180–186, 2020. [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=1001016/j.ijid.2020.06.052&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=http://www.n&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2021%2F12%2F19%2F2021.11.28.21266882.atom) 24. M. S. Mahmud, M. Kamrujjaman, M. M. I. Y. Adan, M. A. Hossain, M. M. Rahman, M. S. Islam, M. Mohebujjaman, and M. M. Molla. Vaccine efficacy and sars-cov-2 control in california and u.s. during the session 2020–2026: A modeling study. Infectious Disease Modelling, 2021. ISSN 2468-0427. doi: [https://doi.org/10.1016/j.idm.2021.11.002](https://doi.org/10.1016/j.idm.2021.11.002). URL [https://www.sciencedirect.com/science/article/pii/S2468042721000701](https://www.sciencedirect.com/science/article/pii/S2468042721000701). 25. Pfizer Inc. Pfizer and BioNTech confirm high efficacy and no serious safety concerns through up to six months following second dose in updated top line analysis of landmark COVID-19 vaccine study. Pfizer, New York, April 1, 2021. ([https://www.pfizer.com/news/press-release/press-release-detail/pfizer-and-biontech-confirm-high-efficacy-and-no-serious.opensinnewtab](https://www.pfizer.com/news/press-release/press-release-detail/pfizer-and-biontech-confirm-high-efficacy-and-no-serious.opensinnewtab)) (press release)., 2021. 26. L. M. Rossen, A. M. Branum, F. B. Ahmad, P. Sutton, and R. N. Anderson. Excess deaths associated with COVID-19, by age and race and ethnicity—United States, January 26–October 3, 2020. Morbidity and Mortality Weekly Report, 69(42):1522, 2020. 27. P.P. Rubió. Is the basic reproductive number (R0) for measles viruses observed in recent outbreaks lower than in the pre-vaccination era? Eurosurveillance, 17(31):20233, 2012. 28. C. M. Stetler. The 1918 Spanish influenza: three months of horror in Philadelphia. Pennsylvania History, 84(4):462–487, 2017. 29. U.S. Census Bureau. QuickFacts United States. [https://www.census.gov/quickfacts/fact/table/US/POP010220](https://www.census.gov/quickfacts/fact/table/US/POP010220), 2020. 30. T. Usherwood, Z. LaJoie, and V. Srivastava. A model and predictions for COVID-19 considering population behavior and vaccination. Scientific Reports, 11(1):1–11, 2021. 31. F. S. Vahidy, L. Pischel, M. E. Tano, A. P. Pan, M. L. Boom, H. D. Sostman, K. Nasir, and S. B. Omer. Real World Effectiveness of COVID-19 mRNA Vaccines against Hospitalizations and Deaths in the United States. medRxiv, 2021. doi: 10.1101/2021.04.21.21255873. [Abstract/FREE Full Text](http://medrxiv.org/lookup/ijlink/YTozOntzOjQ6InBhdGgiO3M6MTQ6Ii9sb29rdXAvaWpsaW5rIjtzOjU6InF1ZXJ5IjthOjQ6e3M6ODoibGlua1R5cGUiO3M6NDoiQUJTVCI7czoxMToiam91cm5hbENvZGUiO3M6NzoibWVkcnhpdiI7czo1OiJyZXNpZCI7czoyMToiMjAyMS4wNC4yMS4yMTI1NTg3M3YyIjtzOjQ6ImF0b20iO3M6NTA6Ii9tZWRyeGl2L2Vhcmx5LzIwMjEvMTIvMTkvMjAyMS4xMS4yOC4yMTI2Njg4Mi5hdG9tIjt9czo4OiJmcmFnbWVudCI7czowOiIiO30=) 32. W. Yang and J. Shaman. Development of a model-inference system for estimating epidemiological characteristics of sars-cov-2 variants of concern. Nature communications, 12(1):1–9, 2021. 33. W. Yang, S. K. Greene, E. R. Peterson, W. Li, R. Mathes, L. Graf, R. Lall, S. Hughes, J. Wang, and A. Fine. Epidemiological characteristics of the b. 1.526 sars-cov-2 variant. medRxiv, 2021. [1]: /embed/graphic-2.gif [2]: /embed/inline-graphic-1.gif [3]: /embed/graphic-7.gif [4]: /embed/graphic-10.gif [5]: /embed/graphic-11.gif [6]: /embed/graphic-12.gif [7]: /embed/graphic-13.gif [8]: /embed/graphic-14.gif [9]: /embed/graphic-15.gif [10]: /embed/graphic-16.gif [11]: /embed/inline-graphic-2.gif [12]: /embed/inline-graphic-3.gif [13]: /embed/graphic-17.gif [14]: /embed/graphic-18.gif [15]: /embed/graphic-19.gif [16]: /embed/graphic-20.gif [17]: /embed/graphic-21.gif [18]: /embed/graphic-22.gif [19]: /embed/graphic-23.gif [20]: /embed/graphic-24.gif [21]: /embed/graphic-26.gif [22]: /embed/graphic-27.gif [23]: /embed/graphic-28.gif