key: cord-0966062-e3ve6w4q
authors: Li, Rui; Li, Yan; Zou, Zhuoru; Liu, Yiming; Li, Xinghui; Zhuang, Guihua; Shen, Mingwang; Zhang, Lei
title: Evaluating the Impact of SARS-CoV-2 Variants on the COVID-19 Epidemic and Social Restoration in the United States: A Mathematical Modelling Study
date: 2022-01-10
journal: Front Public Health
DOI: 10.3389/fpubh.2021.801763
sha: 9d0e718273a17bededa4f570f9a7da13f4dd9d1a
doc_id: 966062
cord_uid: e3ve6w4q

Background: Multiple SARS-CoV-2 variants are still rampant across the United States (US). We aimed to evaluate the impact of vaccination scale-up and potential reduction in the vaccination effectiveness on the COVID-19 epidemic and social restoration in the US. Methods: We extended a published compartmental model and calibrated the model to the latest US COVID-19 data. We estimated the vaccine effectiveness against the variant and evaluated the impact of a potential reduction in vaccine effectiveness on the epidemics. We explored the epidemic trends under different levels of social restoration. Results: We estimated the overall existing vaccine effectiveness against the variant as 88.5% (95% CI: 87.4–89.5%) with the vaccination coverage of 70% by the end of August, 2021. With this vaccine effectiveness and coverage, there would be 498,972 (109,998–885,947) cumulative infections and 15,443 (3,828–27,057) deaths nationwide over the next 12 months, of which 95.0% infections and 93.3% deaths were caused by the variant. Complete social restoration at 60, 65, 70% vaccination coverage would increase cumulative infections to 1.6 (0.2–2.9) million 0.7 (0.1–1.2) million, and 511,159 (110,578–911,740), respectively. At same time it would increase cumulative deaths to 39,040 (5,509–72,570), 19,562 (3,873–35,250), 15,739 (3,841–27,638), respectively. However, if the vaccine effectiveness were reduced to 75%, 50% or 25% due to new SARS-CoV-2 variants, there would be 667,075 (130,682–1,203,468), 1.7 (0.2–3.2) million, 19.0 (5.3–32.7) million new infections and 19,249 (4,281–34,217), 42,265 (5,081–79,448), 426,860 (117,229–736,490) cumulative deaths to occur over the next 12 months. Further, social restoration at a lower vaccination coverage would lead to even greater secondary outbreaks. Conclusion: Current COVID-19 vaccines remain effective against the SARS-CoV-2 variant, and 70% vaccination coverage would be sufficient to restore social activities to a pre-pandemic level. Further reduction in vaccine effectiveness against SARS-CoV-2 variants would result in a potential surge of the epidemic. Multiple measures, including public health interventions, vaccination scale-up and development of a new vaccine booster, should be integrated to counter the new challenges of new SARS-CoV-2 variants.

Conclusion: Current COVID-19 vaccines remain effective against the SARS-CoV-2 variant, and 70% vaccination coverage would be sufficient to restore social activities to a pre-pandemic level. Further reduction in vaccine effectiveness against SARS-CoV-2 variants would result in a potential surge of the epidemic. Multiple measures, including public health interventions, vaccination scale-up and development of a new vaccine booster, should be integrated to counter the new challenges of new SARS-CoV-2 variants.

Keywords: COVID-19 variants, vaccine effectiveness, social restoration, mathematical modelling, COVID-19 vaccination BACKGROUND The COVID-19 pandemic remains a severe public health challenge despite the extensive public health interventions implemented worldwide. By 31st May 2021, accumulatively 170.1 million infected cases and 3.6 million COVID-19-related deaths were reported worldwide (1) . Recent rollouts of effective COVID-19 vaccines have raised hope to control the pandemic. Countries leading the vaccination efforts have seen declining new infections and begun to relax in social distancing and travel bans. Experiences can be drawn from the United Kingdom (UK), Israel, and some states in the US (Vermont, Massachusetts, Hawaii, California), where vaccination coverage has reached 55-70% recently (2) (3) (4) . The vaccination program has successfully vaccinated more than 34 million people in the UK and prevented more than 10,000 deaths, enabling the government to ease previous tough restrictions (5) . Despite the promising progress, the goal of controlling or eliminating the COVID-19 epidemic is still far-reaching, primarily due to the insufficient vaccination coverage worldwide (6) . Besides, the emergence of new variants may hamper the effectiveness of the existing vaccines and shadow the potential population benefits of the vaccination.

Mutations in the receptor-binding domain (RBD) of SARS-CoV-2 enabled the virus to bind more effectively with the host receptor-Angiotensin converting enzyme 2 (ACE2) and enhanced better integration of the virus into the host (7) . The B.1.1.7 strain, with an N501Y substitution in RBD, has shown a 59% higher transmissibility and a 45% higher mortality rate compared to the wild type (8, 9) . Laboratory and clinical studies demonstrated that most COVID-19 vaccines remained effective against B.1.1.7 (10) (11) (12) (13) . However, the additional E484K substitution of RBD in B.1.1.7 and the emergence of more transmissible Delta variant B.1.167.2 cast doubts on the effectiveness of the existing vaccines (14) (15) (16) .

The US has experienced multiple waves of severe COVID-19 epidemics in the past (1). In response, the US has developed effective vaccines against COVID-19 (17) (18) (19) . The US has seen a rapid rollout of a comprehensive vaccination program since 13th December 2020 (20) . By 31st May 2021, the vaccination coverage in the US population has reached 51% (21) . The US has established a comprehensive surveillance system that closely tracked the spread of the COVID-19 variants (22 (24) (25) (26) . Studies also indicated that a population vaccination coverage of at least 50-80% is required to reduce the effective reproductive number to below one and enables restoration of social activities to a pre-pandemic level (27, 28) . Some studies have evaluated the transmission of hypothetical SARS-CoV-2 variants before the actual report of B.1.1.7 in the US (29, 30) . They have substantially underestimated either the transmissibility of the variant or the update of vaccination in the US in these models, leading to a biassed projection of the epidemic trends. A more precise evaluation of the impact of vaccination and new variants is warranted to inform the future epidemic and relevant public health interventions.

The emergence of new COVID-19 variants has complex implications on the epidemic and social restoration due to its potential risks of reducing vaccine effectiveness and protection. Our study aims to evaluate the epidemic trend of the COVID-19 epidemic amid increasing vaccination coverage in the US. We further explore the potential population impact on the epidemic trends if further new variants emerge with higher transmissibility, mortality, and lower vaccine effectiveness. Our findings will help inform public health measures for epidemic management of new variants in the US.

We collected publicly available reportable epidemiological data in the US from Johns Hopkins University Coronavirus Resource Center (1) and the Centers for Disease Control and Prevention (CDC) (21, 23) . These websites provided daily confirmed COVID-19 infection cases and death cases from 1st March 2020 to 31st May 2021, and the number of COVID-19 vaccine uptakes and variant proportions from 13th December 2020 to 31st May 2021. All four types of data were used to calibrate the model (details in Supplementary Materials).

We extended a previously published dynamic compartmental model (27) to describe the circulation of two SARS-CoV-2 strains (the wild type and the variant) in the US. Our model consisted of 18 compartments (Supplementary Figure 1) . A susceptible or vaccinated individual (S, V) may be infected by a SARS-CoV-2 strain (either the wild type or the variant) and entered a latent infection stage (E, E m , m denotes the variant). After a mean incubation period of 5.2 (4.1-7.0) days (31), a proportion of infected individuals developed symptoms (I 1 , I m,1 ) before being diagnosed and documented (T 1 , T m,1 ). The remaining asymptomatic infections (A, A m ) would spontaneously recover (R, R m ). Undiagnosed and diagnosed infected individuals may progress to severe/critical stage (I 2 , I m,2 , T 2 , T m,2 ) and die (D, D m ) or recover (R, R m ). The description of progression rates, derivation of the force of infection and mathematical expressions were documented in the Supplementary Materials.

We assumed that the dominant variant was 59% (56-63%) more transmissible and 45% (18-78%) more deadly than the wild type SARS-CoV-2 (8, 9) . We also assumed the recovered individuals could not be reinfected by any strains of SARS-CoV-2, suggesting a complete cross-protection and no co-infection of SARS-CoV-2 strains. In this study, we calculated the weightedaverage effectiveness (91.4%) of the three available vaccines, Pfizer (92.6%), Moderna (92.1%), and J&J Jensen (66.9%), based on their population coverage in the US (17) (18) (19) 32) .

WE calibrated the model using a non-linear least-squares method that minimised the Root Mean Squared Error (RMSE) between model-simulated and reported data. Based on the "calibrated" scenario, we perturbed model parameters around the "calibrated" parameter set to generate a band of curves that best describe the data variations and retain a minimal level of RMSE. We randomly generated 200 small "perturbating factors." For each of the perturbing factors, we randomly sampled 100 parameter sets based on Latin Hypercube Sampling (LHS) between the parameter range generated by a random walk (adding the perturbing factor in both positive and negative directions). Hence, we obtained 200 groups of various perturbance, and each group has 100 randomised parameter sets. For each of the 200 groups, we calculated the number of data points covered by the band of curves simulated by the 100 parameter sets and their RMSE. We hence selected the one with the minimal RMSE across 200 bands as the set of simulations that best explained the observed data. The 100 curves in the selected band were used to calculate the 95% CI of the model outcomes. The vaccine effectiveness against the variant was estimated spontaneously during model calibration. We validated that the estimated effectiveness of 88.5% produced the lowest RMSE in Supplementary Figure 2 .

The emergence of new variants (e.g., E484K substitution in the B.1.1.7, B.1.167.2 variant) may potentially reduce the effectiveness of the existing vaccines. We evaluate its impact on the COVID-19 epidemic when the vaccine effectiveness (1) adopts the model-estimated value for the variant (baseline scenario); reduces to (2) 75%; (3) 50% and (4) 25%. For each scenario, we calculated the cumulative infections and deaths due to COVID-19 over the next 12 months (1st June 2021 to 31st May 2022). Since we did not know the viral properties of potential new variants, we also simulated nine scenarios with varying viral transmissibility and mortality as a sensitivity analysis (Figure 2 ).

Social restoration (public person-to-person contact rates recovered to 100% of the pre-pandemic level and Mask usage returned to 0) would significantly enlarge the force of infection for both wild type and the variant. To assess the timing of social restoration and herd immunity, we conducted three scenarios of social restoration at 60, 65, and 70% vaccination coverage levels and calculated the daily and cumulative COVID-19 infections and deaths over the next 12 months (Figure 3 ).

Based on the selected 100 parameter sets in model calibration, we produced the sensitivity analysis to accommodate the uncertainty of model parameters and determine the 95% CI of the cumulative COVID-19 infections and deaths. In addition, we also explore the impact on the epidemic trends of COVID-19 in several scenarios (new variants emerge with higher transmissibility, mortality, and lower vaccine effectiveness; social restoration; Figures 2, 3) . All analyses and simulations were performed in MATLAB R 2019b.

Our calibrated model estimates the overall existing vaccine effectiveness against the variant to be 88.5% (95% CI: 87.4-89.5%) (Figure 1 , additional validation in Supplementary Figure 2) , slightly lower than the weighted average effectiveness (91.4%) among the three vaccines available in the US. Our model also demonstrates that the estimated variant trend closely resembled the reported data ( Figure 1C ).

Despite the rise of the variant, our model predicts that the increasing vaccination and persistent public health interventions would control the COVID-19 epidemic over the next 12 months (1st June 2021 to 31st May 2022). If the current trend of vaccination is sustained, vaccination coverage (having received ≥1 doses) in the US will exceed 70% by the end of August 2021 and a plateau around 75% at the end of November Figure 1D) 

The emergence of a more transmissible COVID-19 variant (e.g., B.1.167.2, up by 90% of the wild type) than the current dominant variant (up by 59% of the wild type) but share a similar mortality rate and vaccine effectiveness will result in 694,193 (143,823-1,244,563) cumulative infections and 19,861 (4,579-35,143) cumulative deaths over the next 12 months (Figure 2) . In contrast, if its transmissibility is only half of the variant (up by 30% of the wild type), the estimated number of cumulative infections and deaths will only be 381,684 (89,438-673,930) and 12,824 (3,371-22,276) over the next 12 months. However, we acknowledge the latter scenario may not happen due to the less competitive nature of the variant than the current one.

The emergence of a more deadly COVD-19 variant (up by 75% of the wild type) than the current variant (up by 45% of the wild type) but share a similar transmissibility rate and vaccine effectiveness will result in more deaths but fewer infections, estimated 495,830 (109,604-882,056) cumulative infections and 18,197 (4,393-32,002) cumulative deaths over the next 12 months (Figure 2) . In contrast, a reduced mortality (up by 15% of the wild type) will only cause a similar number of new infections and less deaths (12,615 [3,250-21,980] ).

Retaining high vaccine effectiveness enables a sooner social restoration (Figure 3) . We project that, with the estimated effectiveness (88.5% against the variant), restoring social activity to the pre-pandemic level at the 60%, 65% and 70% vaccination coverage will result in 1.6 (0.2-2.9) million, 0.7 (0. 

Our study assesses the potential impact of vaccination on the COVID-19 epidemic over the next 12 months, with the consideration of an increasing spread of the potential variants. We identified several key findings. First, our model estimated the current vaccine effectiveness against the variant is about 88.5%. With this effectiveness, we project that the number of new infections would be depleted by 98% in August 2021, as the population vaccination coverage increases to 70% during the same period. This enables potential social restoration to the pre-pandemic level. However, if the vaccine effectiveness against the variant is reduced to 25% due to the E484K mutation or B.1.167.2 variant, our model predicted a further wave of the epidemic may be inevitable. Further, if a 30% more transmissible variant emerges, then vaccine effectiveness needs to be above 50% to prevent a major resurge of the epidemic. Our study narrows the gap of reality and simulation after considering mass vaccination and emerging COVID-19 variants based on previous publications (33) (34) (35) (36) (37) (38) .

Our finding demonstrates the overall vaccine effectiveness for the variant is slightly lower than that for the wild type. This is consistent with previous findings. Laboratory studies based on serum and plasma neutralisation assay with authentic and pseudoviruses of the variant show no significant change in neutralising activity of current vaccines against the variant, despite the presence of RBD mutations (11, 12, 14, (39) (40) (41) (42) . However, Emary et al. (10) find in vitro that neutralising antibody titres generated by AstraZeneca are lower for the variant, but the clinical vaccine effectiveness remains at 70.4%, indicating a lower activity may still be sufficient to confer protection or elicit host immunity. Further studies provide evidence of vaccine protective effects for the variant in a population (13, (43) (44) (45) (46) . Abu-Raddad et al. (43) reported an estimate of 87.0% effectiveness of Pfizer against the variant. Consistently, Haas et al. (44) reported an overall effectiveness of Pfizer vaccine to be 95.3% in Israel, where the epidemic is dominated by B.1.1.7 (94.5%). Lopez Bernal et al. (46) estimated a combined 60% vaccine effectiveness for both Pfizer and AstraZeneca in England, where the variant is first reported. Nevertheless, with this high effectiveness for the variant, our study supports that current vaccination effort will control the COVID-19 epidemic in the US.

Our finding also simulates that a new variant (e.g., Delta) or an E484K-mutated variant may reduce the overall vaccine effectiveness and become dominant in the US. In this case, we project further waves of the epidemic with its size proportional to the level of vaccine effectiveness reduced. To cater for these scenarios, we recommend several potential strategies. First, maintaining the current vaccination effort is of vital importance for COVID-19 control even with reduced vaccine effectiveness. This is because a weakened vaccine still reduces the overall number of new infections and the risk of new variants. Besides, even though a vaccine may not completely protect against the infection, it still significantly reduces the clinical severities of the infected individuals (18, 32) . Second, the pharmaceutical industry should prioritise its efforts in developing new vaccine boosters for the emerging new variants of concerns (47) (48) (49) . At a population level, epidemiological studies need to focus on whether a booster shot, or third dose of the existing vaccines, or a combination of both should be implemented to maximise vaccination protection. Third, orchestrated vaccination efforts worldwide are necessary to ensure less-developed regions can receive sufficient vaccines for epidemic control and hence reduce the risk of new COVID-19 variants.

Our finding indicated that social restoration is strongly dependent not only on the vaccine effectiveness but also on its coverage. We discovered that 70% vaccination coverage with the current vaccine effectiveness might allow safe restoration of social activities to the pre-pandemic level. However, the decline of vaccine effectiveness would lead to further waves of epidemic and reduce the likelihood of social restoration. In fact, if the vaccine effectiveness reduces to 25%, current public health interventions in the US would need to be further strengthened to control the epidemic. As the widespread SARS-CoV-2 would result in a higher likelihood of variants and reduce vaccine effectiveness, which in turn causes further spread of the virus (7, 50, 51), a high vaccination coverage will limit the viral spread and break the vicious cycle. Therefore, encouraging the community to vaccinate is essential to reduce variants and retain vaccine effectiveness for social restoration. Authorities should enhance health promotion to reduce the public's misunderstanding about vaccination and provide free and accessible vaccination (52) (53) (54) . Besides, relaxing restrictions for fully vaccinated individuals in public places and travelling is also an alternative to encourage the community to vaccinate.

This study has several limitations. First, our model did not consider age structure and variations in the risk of infection and mortality across age groups. Second, our model structure only investigated one SARS-CoV-2 variant. We have mostly focused on the B.1.1.7 variant since B.1.1.7 is currently the most reported variant in the US but only included other variants as a part of sensitivity analysis. The latest report has indicated that Beta (B.1.351), Gamma (P.1) variants and Delta (B.1.617.2) have accounted for 11.9% of new infections, and the proportion is growing (23) . Modelling the competition of variants would be our future investigation. Third, we assumed that the natural immunity and immunity elicited by vaccines was not wane during our simulations. If the neutralising antibody reduction results in reduced vaccine protection, more strict public health interventions or a new vaccine booster for enhancing the immune response may be necessary. Finally, the model did not consider the impact of the individual willingness of vaccination on herd immunity.

In conclusion, our modelling exercise indicates that the current vaccines in the US remain effective for the variant, and 70% of vaccination coverage would be sufficient to restore social activities to a pre-pandemic level if no new, more transmissible variants emerge. The emergence of new, more transmissible variants accompanied by the uncertain impact on vaccine effectiveness would potentially result in further waves of the epidemic. Our findings confirm that multiple measures, including public health interventions, vaccination scale-up and development of a new vaccine booster, should be integrated to counter the new challenges of new SARS-CoV-2 variants.

The datasets presented in this study can be found in online repositories. The names of the repository/repositories and accession number(s) can be found in the article/Supplementary Material.

LZ conceived the study. LZ, MS, YL, and RL designed and constructed the model. RL performed the modelled analyses and graphed and interpreted the results. RL, ZZ, YML, and XL contributed to the collection of data and model parameters. RL drafted the manuscript. LZ, MS, and GZ critically revised the manuscript. All authors reviewed the manuscript and approved the final version. 

The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fpubh. 2021.801763/full#supplementary-material

JHU. Coronavirus Resource Center: COVID-19 Tracking

Restrictions Are Easing Across the UK. (2021)

Available online at

Israel to Lift Most COVID-19 Restrictions in June-Health Minister

COVID-19 Vaccination and Blood Clotting

Available online at: https:// ourworldindata.org/covid-vaccinations

The variant gambit: COVID-19's next move

Estimated transmissibility and impact of SARS-CoV-2 lineage B.1.1.7 in England

Mortality risk in patients infected with SARS-CoV-2 of the lineage B.1.1.7 in the UK

Efficacy of ChAdOx1 nCoV-19 (AZD1222) vaccine against SARS-CoV-2 variant of concern 202012/01 (B.1.1.7): an exploratory analysis of a randomised controlled trial

Neutralizing activity of BNT162b2-elicited serum

Antibody resistance of SARS-CoV-2 variants B.1.351 and B.1.1.7. Nature

BNT162b2 Vaccination effectively prevents the rapid rise of SARS-CoV-2 variant B.1.1.7 in high risk populations in Israel

Sensitivity of SARS-CoV-2 B117 to mRNA vaccine-elicited antibodies

Serum neutralizing activity elicited by mRNA-1273 vaccine

Variant That First Appeared in India Arrives in the US

Safety and efficacy of the BNT162b2 mRNA Covid-19 vaccine

Efficacy and Safety of the mRNA-1273 SARS-CoV-2 vaccine

Safety and efficacy of single-dose Ad26.COV2.S vaccine against Covid-19

COVID-19 Vaccination Program Operational Guidance

COVID-19 Vaccinations in the United States

Tracking of Variants

Variant Proportions in the US

The impact of vaccination on COVID-19 outbreaks in the United States

The Impact of Vaccination to Control COVID-19 Burden in the United States: a Simulation Modeling Approach

Epidemiological differences in the impact of COVID-19 vaccination in the United States and China. Vaccines

Projected COVID-19 epidemic in the United States in the context of the effectiveness of a potential vaccine and implications for social distancing and face mask use

Stability analysis of SEIR model related to efficiency of vaccines for COVID-19 situation

Accelerated vaccine rollout is imperative to mitigate highly transmissible COVID-19 variants

Modeling of future COVID-19 cases, hospitalizations, and deaths, by vaccination rates and nonpharmaceutical intervention scenarios -United States

Early transmission dynamics in Wuhan, China, of novel coronavirus-infected pneumonia

Safety and efficacy of the BNT162b2 mRNA Covid-19 vaccine

Forecasting the daily and cumulative number of cases for the COVID-19 pandemic in India. Chaos: an interdisciplinary

Mathematical modeling of the COVID-19 pandemic with intervention strategies

Mathematical modeling and optimal intervention strategies of the COVID-19 outbreak

A mathematical model for COVID-19 transmission dynamics with a case study of India

Effects of New York's executive order on face mask use on COVID-19 infections and mortality: a modeling study

What is required to prevent a second major outbreak of sars-cov-2 upon lifting quarantine in wuhan city

Multiple SARS-CoV-2 variants escape neutralization by vaccine-induced humoral immunity

Neutralization of SARS-CoV-2 lineage B117 pseudovirus by BNT162b2 vaccine-elicited human sera

Inactivated COVID-19 vaccine BBV152/COVAXIN effectively neutralizes recently emerged B 1.1.7 variant of SARS-CoV-2

SARS-CoV-2 variant B.1.1.7 is susceptible to neutralizing antibodies elicited by ancestral spike vaccines

Effectiveness of the BNT162b2 Covid-19 vaccine against the B.1.1.7 and B.1.351 variants

Impact and effectiveness of mRNA BNT162b2 vaccine against SARS-CoV-2 infections and COVID-19 cases, hospitalisations, and deaths following a nationwide vaccination campaign in Israel: an observational study using national surveillance data

Effectiveness of BNT162b2 vaccine against the B.1.1.7 variant of SARS-CoV-2 among healthcare workers in Brescia

Effectiveness of the Pfizer-BioNTech and Oxford-AstraZeneca vaccines on covid-19 related symptoms, hospital admissions, and mortality in older adults in England: test negative case-control study

Booster Shot Generates Promising Immune Response Against Variants Found in South Africa, Brazil

People With Underlying Conditions to Be First to Get Covid Vaccine Booster Shots

Cost-effectiveness analysis of BNT162b2 COVID-19 booster vaccination in the United States

Novel SARS-CoV-2 variants: the pandemics within the pandemic

Evaluation, and Treatment of Coronavirus (COVID-19)

Researchers Propose 5 Strategies to Encourage COVID-19 Immunization

Achieving positive mental health and wellbeing on the COVID-19 frontline. Innovation (2020) 1:100024

Systemic approach in burnout prevention: a qualitative study on wellbeing of china deployed healthcare workers during COVID-19 pandemic

The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.