Investigating Vaccine-induced Immunity And Its Effect in Mitigating SARS-CoV-2 Epidemics in China
Nov 10, 2023
Abstract
Background: To allow a return to a pre-COVID-19 lifestyle, virtually every country has initiated a vaccination program to mitigate severe disease burden and control transmission. However, it remains to be seen whether herd immunity will be within reach of these programs.
Methods: We developed a compartmental model of SARS-CoV-2 transmission for China, a population with low prior immunity from natural infection. Two vaccination programs were tested and model-based estimates of the immunity level in the population were provided.
cistanche tubulosa-improve immune system
Results: We found that it is unlikely to reach herd immunity for the Delta variant given the relatively low efficacy of the vaccines used in China throughout 2021 and the lack of prior natural immunity. We estimated that assuming a vaccine efficacy of 90% against the infection, vaccine-induced herd immunity would require coverage of 93% or higher of the Chinese population. However, even when vaccine-induced herd immunity is not reached, we estimated that vaccination programs can reduce SARS-CoV-2 infections by 50–62% in case of an all-or-nothing vaccine model and an epidemic starts to unfold on December 1, 2021.
Conclusions: Efforts should be taken to increase the population’s confidence and willingness to be vaccinated and to develop highly efficacious vaccines for a wide age range.
cistanche tubulosa-improve immune system
Keywords: COVID-19, Herd immunity, Vaccination program, Delta variant, SLIR model
Background
The first wave of novel coronavirus disease 2019 (COVID-19) in China subsided quickly after the implementation of strict containment measures and travel restrictions starting in March 2020 [1–4]. As of November 12, 2021, the COVID-19 pandemic has caused over 251 million reported cases and 5 million deaths globally [5]. The pandemic is far from over, as severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) has undergone some significant mutations and a number of variants have become widespread due to increased transmissibility and/or immune escape characteristics—e.g., variants Alpha [6–12], Beta [13, 14], Gamma [13, 15], and Delta [16–18]. Throughout the globe, a rapid surge of Delta variant cases suggests a clear competitive advantage compared with Alpha, Beta, and Gamma [16]; more than 90% of daily sequences from global initiatives on sharing all influenza data (GISAID) are ascribable to the Delta variant since July 2021 [19]. Despite no major epidemics, China has been experiencing several minor local outbreaks caused by imported cases of the Delta variant, including the outbreaks in Guangzhou, Nanjing, and Zhengzhou city [20–22]. To suppress transmission, a large share of the world needs to have immunity to SARS-CoV-2, especially to the Delta variant. Effective vaccines against COVID-19 represent the most viable option to suppress SARS-CoV-2 transmission globally. The effectiveness of vaccination programs depends on several key factors, including vaccine supply, willingness to receive the vaccine, vaccine efficacy, and the age groups targeted by the vaccination effort. However, current vaccination programs are all based on vaccines developed against the original SARS-CoV-2 lineage, and the efficacy seems to be reduced against the Delta variant [23]. In China, home to about 1.4 billion people (~18% of the world population), 2.37 billion doses have been administered as of November 12, 2021 [24]; 76.5% of the whole population has been vaccinated with two doses, corresponding to 82.4% of the target population (i.e., individuals aged 3 years and older). However, it remains to be seen if the vaccine coverage may reach a level sufficient to achieve herd immunity. Countries around the globe are facing the same question. The classical herd immunity level is defined as 1-1/R0, where R0 is the basic reproduction number—the average number of infections generated by a typical infectious individual in a fully susceptible population [25]. For a vaccine with efficacy VE that gives life-long protection, the level of herd immunity required to stop transmission is (1-1/R0)/VE. However, this estimate is an oversimplification of a complex phenomenon as it ignores the heterogeneities of the actual human population (e.g., social mixing patterns, age-specific susceptibility to infection) [25, 26] as well as of vaccination (e.g., lifelong immunity, sterilizing vaccine). To overcome this limitation, here we integrate a contact survey specific to the Chinese population [27] as well as official demographic statistics to develop an age-structured stochastic model to simulate SARS-CoV-2 transmission (Additional file 1: Fig. S1). We then use this model to evaluate whether herd immunity is achievable against the Delta variant or not via mass vaccination.
cistanche tubulosa-improve immune system
Click here to view Cistanche Enhance Immunity products
【Ask for more】 Email:cindy.xue@wecistanche.com / Whats App: 0086 18599088692 / Wechat: 18599088692
Methods
SARS-CoV-2 transmission and vaccination model
We built a compartmental model of SARS-CoV-2 transmission and vaccination, based on an age-structured stochastic susceptible-latent-infectious-removed (SLIR) scheme, accounting for heterogeneous contact patterns by age [27] and heterogeneous susceptibility to infection by age as estimated using contact tracing data in Hunan province of China [28]. In the model, the population is divided into four epidemiological categories: susceptible, latent, infectious, and removed, stratified by 16 age groups. Susceptible individuals can become infected after contact with an infectious individual according to the age-specific force of infection. The rate at which contacts occur is determined by the mixing patterns of each age group. The latent period and average generation time were set to be 4.4 [28–30] and 7 [31] days, respectively. We consider a basic reproductive number (R0) of 6.0 according to estimates for the SARS-CoV-2 Delta variant [1–4, 6–12, 16–18]. Simulations are initiated with 40 infectious individuals [32], corresponding to the number of cases first detected in a local outbreak in Beijing on June 11, 2020. We consider a 2-dose vaccine that only susceptible individuals are eligible for vaccination (we recall that natural immunity is close to 0 in China as of November 2021 [33]) and that the duration of vaccine-induced immunity lasts longer than the time horizon considered in this study (i.e., 1 year). Details about the model and parameters are reported in Additional file 1: Sec. 1 and Tab. S1.
Baseline scenario
As the baseline scenario, we considered the following assumptions:
i) Epidemic seeding: An epidemic is assumed to be triggered by 40 SARS-CoV-2 infectious individuals on December 1, 2021 [32].
ii) Vaccination strategy: Vaccines have been rolling out in China since November 30, 2020 [24], which is the earliest date reported by the government, and have been extended to individuals aged 3+ years since early November, and we tested two different vaccination strategies:
a) Strategy 1—random distribution of vaccines to individuals aged 12+ years, then extended to individuals aged 3+ years starting from November 1, 2021;
b) Strategy 2—random distribution of vaccines to individuals aged 3+ years since the start of the vaccination program, namely November 30, 2020.
We considered that a fraction of the population (about 2%—Additional file 1: Tab. S2) is not eligible to receive the vaccine as pregnant women and individuals with allergies or other conditions are excluded from the vaccination campaign in China as of November 2021 [34–40] (see Additional file 1: Sec. 2 for detail).
iii) Vaccine capacity: We used the historical data on daily administrated doses in China until November 2, 2021 (Additional file 1: Fig. S2), then projected the future daily capacity based on the average daily doses administrated over the period October 2– November 2, 2021 [24]. Thus, from November 3, 2021, and beyond, we estimated a daily vaccine administration capacity of 2.30 million doses for the China population (details are reported in Additional file 1: Sec. 3).
iv) Vaccine efficacy: The vaccine schedule requires two doses with a 21-day interval. VE against infection for individuals aged 18–59 years old reaches its maximum value 14 days after vaccinating 2 doses and is estimated at 54.3% for the Delta variant [41–43]. This estimate is based on the efficacy measured against the original lineages and the reduction of neutralizing antibodies estimated for the Delta variant in clinical studies (see Additional file 1: Tab. S1 for details). The relative VE against infection within 0– 13 days after the second dose compared with maximum protection is 83.8% for the Delta variant [44]. VE against death for individuals aged 18–59 years old is 93% for Delta variant [45–47]. We explored higher VE values against infection [48] and tested a two-dose schedule with a 14-day interval as sensitivity analyses (Additional file 1: Tab. S1). In addition, COVID-19 vaccines may not be equally effective across age groups in preventing infection. To understand the impact of this assumption, we also tested a relative VE of 50% and 75% for individuals aged 3–17 and 60+ years as compared to VE for individuals aged 18–59 years.
v) Vaccine action: We considered two mechanisms to model vaccine efficacy: an “all-or-nothing” vaccine (baseline analysis), in which the vaccine provides full protection to a fraction VE of individuals who are vaccinated and no protection to the remaining 1-VE vaccinated individuals. The second option we considered is a “leaky” vaccine in which all vaccinated individuals have a certain level of protection from the infection corresponding to VE [49].
vi) Initial immunity: As of November 2021, there is essentially no population immunity from natural infection in China [33]. For the sake of generalizability of the results to other countries that had widespread transmission, we explored a scenario where 30% of the population is initially immune to the infection, and the fraction of immune individuals by age group is proportional to the population size.
vii) Susceptibility to infection by age: Children under 15 years of age are estimated to have a lower susceptibility to SARS-CoV-2 infection as compared to adults (i.e., individuals aged 15 to 64 years), while individuals aged 65+ years have the highest susceptibility to infection [28].
viii)Immunity duration: We let the transmission model run for 1 year, assuming life-long protection from natural infection or vaccination.
ix) Disease burden: The infection fatality ratio for original lineages manifests in 0.0923% for individuals aged 0–19, rising to 6.7959% for individuals aged over 80 years [50, 51]. The risk of death associated with the Delta variant compared to the original lineages is 2.37 [52].
Comprehensive sensitivity analyses to evaluate the impact of the baseline assumptions on our results are carried out as well (Additional file 1: Tab. S1).
cistanche tubulosa-improve immune system
Alternative vaccination scenarios
We tested three alternative scenarios to explore the potential for vaccination-induced herd immunity, where (i) the start of the epidemic is delayed from December 1, 2021, to January 1, 2022, and February 1, 2022; (ii) the value of the reproduction number in a fully susceptible population and under a certain level of non-pharmaceutical interventions (NPIs), denoted as RNPIs 0, varies between 1.1 and 6; (iii) combinations of scenarios I and ii. For scenario (ii), we did not explicitly model single NPIs such as case isolation, contact tracing, wearing masks, social distancing, and improved hygiene. Instead, the synergetic effect of these measures was considered as a reduction of the reproduction number.
Data analysis
For each scenario, 100 stochastic simulations were performed, and the mean and 95% confidence interval (95% CI) were then estimated.
We used the next-generation matrix (NGM) [53] approach to estimate the effective reproduction number, Re. Herd immunity is considered achievable when Re <1. Details are reported in Additional file 1: Sec. 4 and 5.
Results
Baseline scenario
By forward simulating 1 year of epidemic and assuming no vaccine hesitancy, continued vaccination efforts would lead to a final coverage of 90.7% of the target population, which corresponds to 88.6% of the total population for strategy 1 (Fig. 1a). For strategy 2, the estimated coverage of the total population is 88.7% (Fig. 1b). Under any scenario, the mean daily incidence never reaches 250 over 10,000 residents (Fig. 1c, d). We estimated that the effective reproduction number at the time the infection is seeded (Re) is still well above the epidemic threshold, namely 4.03 (95% CI 3.19–4.70) and 3.18 (95% CI 3.15–3.24) for strategy 1 and 2, respectively (Fig. 1e, f). These estimates suggest that the vaccine coverage on December 1, 2021, is not enough to prevent onward transmission, regardless of the vaccination strategy. He is estimated to cross the epidemic threshold (i.e., 1) on January 31 and February 5, 2022, for strategies 1 and 2, respectively, due to the accumulation of immune individuals both through continued vaccination efforts and natural infections (Fig. 1g, h). The estimated infection attack rates (which include all SARS-CoV-2 infections, independent of whether an individual develops symptoms or not) are 45.2% (95% CI 37.2–48.5%) and 45.7% (95% CI 40.6–48.8%) for strategies 1 and 2, respectively (Fig. 1g, h). Note that the proportion of vaccine-immune individuals stops increasing while the proportion of naturally immune individuals is still increasing. In fact, when all individuals are either vaccinated or infected, the proportion of vaccine-immune individuals will stop increasing. However, in this situation, unprotected/partially protected vaccinated individuals can still be infected, which leads to an increase in the proportion of naturally immune individuals.
Fig. 1 Time series of vaccine coverage, daily incidence, effective reproductive number, and proportion of immune individuals. an Age-specific vaccine coverage over time for strategy 1. The dotted lines correspond to the start of the epidemic. The inserted table shows the age-specific coverage for the two key time points (the start of the epidemic (i.e., December 1, 2021) and the time that the coverage remains constant (i.e., March 11)). The line corresponds to the mean value, while the shaded area represents 95% CI. b As a, but for strategy 2. c Daily incidence per 10,000 for strategy 1 (mean and 95% CI). d As c, but for strategy 2. e Effective reproduction number Re over time (mean and 95% CI) for strategy 1. The shaded area in gray indicates the epidemic threshold Re =1. The numbers around the shaded area indicate when Re crosses this threshold (i.e., January 31) for strategy 1. f As e, but for strategy 2. g Proportion of immune individuals due to either natural infection or vaccination over time for strategy 1. h As g, but for strategy 2
Although vaccine-induced immunity is not enough to prevent viral circulation, all the scenarios considered are associated with substantial mitigation of the COVID-19 burden. We estimate the infection attack rates for the two vaccination strategies to decrease by more than 50% with respect to a reference scenario with no interventions (Fig. 2a, b). Both strategies lead to more than 90% reduction in the number of deaths (Fig. 2c, d). These results were based on the assumption of an “all-or-nothing” vaccine. To test the robustness of our findings to this assumption, we tested a “leaky” vaccine. In this case, we estimated a lower reduction of the infection attack rate (12% as compared to about 50%); however, we estimated a similar reduction in the number of deaths (about 85% as compared to about 90%), Fig. 2e–h. The obtained results show that herd immunity cannot be reached by December 1, 2021, regardless of the adopted vaccination strategy when the R0 is set at 5 or 7 (Additional file 1: Fig. S3), when the initial number of seeds is varied in the range from 10 to 100 (Additional file 1: Fig. S4), and when equal susceptibility to infection by age is assumed (Additional file 1: Fig. S5). The same conclusion is obtained when we considered a more parsimonious model with 3 age groups (Additional file 1: Fig. S6). Finally, we also conducted a counterfactual analysis where we assume that a part of the population was already immune before the start of the vaccination campaign (similar to the situation in Western countries). Under this assumption, we found that a 30% initial immunity proportion would not lead to Re below the epidemic threshold for two strategies before December 1, 2021, (Additional file 1: Fig. S7). As regards the parameters regulating the vaccination process, we found that the vaccine efficacy 14 days after the second dose has the largest impact, followed by the vaccine efficacy of individuals aged 3–17 and 60+ relative to individuals aged 18–59 years (Additional file 1: Fig. S8 and S9). On the other hand, the relative vaccine efficacy within 0–13 days after the second dose and the time interval between the first and second dose have a more moderate impact on the overall effectiveness of the analyzed vaccination strategies (Additional file 1: Fig. S10 and S11).
Fig. 2 Disease burdens of COVID-19 in the baseline scenario. a Cumulative number of infections per 10,000 individuals after 1 simulated year for reference scenario and two vaccination strategies using “all-or-nothing” vaccine model (mean and 95% CI). b Reduction in infections (mean and 95% CI) with respect to the reference scenario in different age groups and the total population. The 95% CI of the reduction may cross 0 as the burden between the reference scenario and the vaccination scenario is approximately the same in some simulations. We thus trimmed the lower limit of 95% CI at 0 throughout the manuscript. c, d as for a, b, but for death. e–h as for a–d, but for the “leaky” vaccine model
Scenario 1: Delaying the start of the epidemic
The findings presented thus far suggest that herd immunity against the Delta variant cannot be built through vaccination by December 1, 2021. Next, we tested to what extent the start of a new epidemic wave needs to be delayed (e.g., by keeping strict restrictions for international travel) to allow immunity to build up in the population, potentially reaching herd immunity levels. According to the daily vaccine capacity used in the baseline scenario (based on the history of daily vaccination capacity data up to November 2, 2021), we estimated that Re remains above the epidemic threshold for both two strategies even if the seeding of an epidemic is delayed to February 1, 2022 (Fig. 3a), while the reduction in infections increases to 56.8% and 57.4% for strategies 1–2, respectively. It is important to stress that the source of uncertainty in our estimates of Re is the boot-strapped contact matrix by age and the posterior distribution of the susceptibility to infection by age. This explains why the estimated confidence interval of Re for strategy 1 is wider than for strategy 2 (which, in Fig. 3a, is smaller than the size of the dot). In fact, for strategy 2, the vaccination is essentially uniform by age and thus the uncertainty on age-dependent parameters is negligible. On the contrary, for strategy 1, the young population is vaccinated at a later stage, which implies that the uncertainty on age-dependent parameters reflects a larger uncertainty on Re. We also reported the impact of delaying the start of the epidemic on vaccine coverage and daily incidence in Additional file 1: Fig. S12.
Fig. 3 Impact of delaying the start of the epidemic and adopting NPIs. an Effective reproduction number Re (mean and 95% CI) as a function of vaccine coverage at the time when infection is seeded. Colors refer to the scenario of delaying the start of the epidemic to a different date. The shaded area in gray indicates Re ≤1. b Cumulative number of infections per 10,000 individuals after 1 simulated year for reference scenario and two vaccination strategies (mean and 95% CI). c Reduction in infections (mean and 95% CI) with respect to the reference scenario. d As a, but for net reproduction number Rt (mean and 95% CI) adopting different intensities of NPIs. e As b, but for the scenario of adopting different intensities of NPIs. f As c, but for the scenario of adopting different intensities of NPIs
Scenario 2: Adopting NPIs in case of a new outbreak
The results presented so far suggest that herd immunity against the Delta variant is not achievable at any time point. Adopting NPIs as a response to an epidemic outbreak can lower the transmission potential of the virus. It is thus worth investigating the synergetic effect of vaccination programs combined with NPIs of different intensities. It is important to note that we do not explicitly model every single measure to limit transmission (e.g., case isolation, contact tracing, wearing masks, social distancing, improved hygiene). These measures are implicit as concerted strategies that result in a decreased reproduction number. We explored RNPIs 0 in the range 1.1–6.0 corresponding to different intensities of interventions. Values between 1 and 2 are shown in the main text, while larger values are shown in Additional file 1: Fig. S13. We also reported the impact of adopting NPIs in case of a new outbreak on daily incidence in Additional file 1: Fig. S14. The mean net reproduction number (defined as the reproduction number accounting both for immunity and interventions) on December 1, 2021, for strategy 1 can be reduced to below 1 only when RNPIs 0 ≤1.5, while for strategy 2, RNPIs 0 can be up to 1.8 (Fig. 3d). By forward vaccinating and simulating 1 year of epidemic, substantial infections could be reduced (close to 100%) thanks to the synergetic effect of vaccination and NPIs (Fig. 3e, f). Note that the reductions in infections are obviously smaller than 100% for RNPIs 0 ≤1.3, as the number of cumulative infections is extremely low in the reference scenario.
Scenario 3: Delaying the start of the epidemic and adopting NPIs
To further improve the potential for vaccination-induced herd immunity and reduce the COVID-19 burden, here we tested the combination of the two scenarios mentioned above: delaying the start of the epidemic and adopting NPIs of different levels of intensity in response to a new outbreak. Should an epidemic start in December 2021– February 2022, strategies 1 and 2 can succeed in blocking transmission only if moderate NPIs (RNPIs 0 in the range 1.5–2.0) are adopted (Fig. 4). The results of reduction in infections compared with a reference scenario for two strategies are showed in Additional file 1: Fig. S15. The effectiveness of age-targeted vaccination strategies depends on the age-mixing patterns of the population [54]. To test the robustness of our findings, we tested an alternative contact matrix for China [55] and found consistent results (Additional file 1: Fig. S16 and S17).
Herd immunity threshold
Till now, herd immunity is unattainable for any vaccination strategy considering the relatively low efficacy (54.3%) of the analyzed vaccine in preventing the infection from the Delta. We thus explored the potential of herd immunity for the two vaccination strategies given a higher efficacy (95%) (Additional file 1: Fig. S18). We estimated that Re can decrease below 1.0 for two strategies (Additional file 1: Fig. S18a). The estimated herd immunity thresholds under these two strategies are 91.3% and 84.5% respectively, which suggests that the level of immunity needed to lead the effective reproduction number below the epidemic threshold is lower if vaccination is extended to individuals aged 3+ years early on. We also estimated the infection attack rate under different vaccination coverages under the assumption that vaccination stops at the time the epidemic is seeded. This purely hypothetical scenario shows that when individuals aged 12+ years are prioritized (strategy 1), despite a fairly high estimated reproduction number when vaccine coverage equals 80% (3.2), the estimated infection attack rate is relatively low (10.0%) (Additional file 1: Fig. S18b). In fact, given the age-targeted vaccination program and the lack of natural immunity, the susceptible population is mostly concentrated in the young population. The high number of contacts in younger age groups, combined with the high vaccination coverage in the rest of the population, lead to a fairly high reproduction number but, at the same time, the infections are focused on a small segment of the population only (young individuals) and thus the overall infection attack rate remains fairly low. We also explored whether herd immunity is achievable or not and what is the herd immunity threshold by estimating Re under the assumption that all individuals are eligible to be vaccinated and have vaccinated 2 doses before the epidemic started. We considered vaccine efficacy in the range of 60–100% and explored different scenarios on vaccination coverage. Our results show that, for a vaccine with an efficacy lower than 85%, herd immunity is unattainable, even in the extreme case where the vaccine coverage is 100% (Fig. 5). Vaccine-induced herd immunity may only be achievable with higher VE and coverage. For example, for a vaccine with 90% efficacy against infection from the Delta variant, more than 93% of the population would need to be vaccinated to reach herd immunity (Fig. 5). In the presence of NPIs, the net reproduction number can be reduced below the unit for lower vaccine efficacy and coverage values (Additional file 1: Fig. S19).
Fig. 4 Impact of delaying the start of the epidemic start and adopting NPIs on net reproduction number. a Net reproduction number Rt as a function of RNPIs 0 and epidemic start date for strategy 1. The bold line in black indicates Rt =1. b As a, but for strategy 2
Discussion
Our study evaluated the feasibility of reaching herd immunity against the SARS-CoV-2 Delta variant through vaccination, considering heterogeneity in population age structure, age-specific contact patterns, vaccine efficacy, and biological characteristics of SARS-CoV-2, including the basic reproduction number, susceptibility to infection by age, and key time-to-event periods (e.g., latent period, generation time). Our findings show that herd immunity is unlikely to be reached against the Delta variant given the relatively low efficacy of the current vaccines (developed against the original SARS-CoV-2 lineage), also in the presence of prior natural immunity up to 30%. Even considering vaccines with higher efficacy, our results show that extending the vaccination program to young children as soon plays a key role to increase the potential of reaching herd immunity and reduce the infection attack rate. If we consider a protection against the Delta variant of 90% (which goes beyond current vaccines), herd immunity would require the vaccination of 93% of the whole population. The adoption of NPIs could prevent the spread of a major epidemic wave even when the herd immunity level is not reached, but such an option obviously entails social and economic costs. Further, both strategies considered in this study would mitigate the overwhelming majority of infections. Our study explored if and when vaccination-induced herd immunity can be reached in China. Under the hypothesis that the circulating strain has the same transmissibility as the Delta variant and that the vaccination campaign will not slow down due to vaccine hesitancy, herd immunity seems to remain unreachable even in the extreme case where the vaccine coverage is 100%. Nonetheless, it is important to remark that the effectiveness of the vaccination program is impacted both by the natural immunity accumulated in the population (which is close to 0 in China as of November 2021) and the age structure of the population. In fact, in populations with a higher natural immunity level and a lower proportion of children, herd immunity may be achievable.
Our findings pointed to the importance of adopting NPIs and/or self-precautionary measures until herd immunity is reached or the burden of the epidemic becomes manageable. These measures can either help delay the seeding of the infection (e.g., strict border control measures) or should an epidemic start to unfold, mitigate its burden (e.g., social distancing, contact tracing, testing, wearing masks, hygiene practices, limiting contacts). However, questions remain about which NPIs need to be implemented, their intensity, and timing. Future studies are needed to address these questions.
cistanche tubulosa-improve immune system
A key role in determining the success of a vaccination campaign is played by the willingness to vaccine the population. According to previous surveys on COVID-19 vaccine hesitancy, vaccine acceptance in China was estimated to vary between 60.4 and 91.3% for the general population aged 18 years and above [56–59] and maybe even lower for older adults [59]. Similar estimates were obtained for several other countries including the UK (71.5%) [60] and the USA (75.4%) [60]. Given these levels of vaccine hesitancy, achieving high levels of coverage may remain an elusive target. Efforts to increase the population’s confidence and willingness to be vaccinated will thus be of paramount importance to allow a return to a pre-COVID-19 lifestyle. Our study shows that the spread of the more transmissible Delta variant has substantially increased the herd immunity threshold to a level that may not be feasible in any population so mitigation strategies become even more relevant. Previous studies have estimated the herd immunity threshold either through natural infection or vaccination under the assumption of a homogenously mixed population [61–63], but heterogeneity in contact structure, age structure of the population, susceptibility to infection by age, and order in which individuals are vaccinated are all key factors shaping the herd immunity level [25]. To explore the impact of the heterogeneities included in the model on the obtained results, we tested an alternative model based on a fully homogeneous population, thus neglecting the contact structure, age structure of the population, susceptibility to infection by age, and order in which individuals are vaccinated that are accounted for in the main analysis (Additional file 1: Fig. S20). When considering R0=6.0 and vaccine efficacy against the infection = 95%, we estimated the theoretical herd immunity threshold (i.e., not accounting for the waning of immunity and emergence of the new variant with immune escape) to be 87.7% for the homogeneous model as compared to 91.3% and 84.5% of vaccination strategies 1 and 2 for the heterogeneous model. Our developed model is based on social mixing patterns estimated for the Shanghai population [27] and on China-specific data on COVID-19 epidemiology (population immunity, etc.). Nevertheless, the introduced modeling framework is flexible and can be tailored to other countries. We tested a scenario somehow resembling the situation in the USA, where we considered natural immunity [64] and the adoption of the BNT162b2/ Pfizer vaccine, whose efficacy against the Delta variant was estimated at 79% [48]. Also, in this scenario, we estimated that herd immunity may not be reached (Additional file 1: Fig. S21). Moreover, vaccination hesitancy may jeopardize the vaccination effort in the USA and other Western countries as well.
This study is prone to the limitations pertaining to modeling exercises. First, VE against infections from the Delta variant was inferred instead of directly measured from epidemiological observations. Moreover, VE for children has not been estimated for the vaccines in use in China; therefore, we have assumed the same VE as in adults based on immunogenicity studies [65]. Given such a lack of field evidence, we have conducted a sensitivity analysis where a lower vaccine efficacy is assumed for children. The overall conclusions of the study do not change. Still, further data on age-specific vaccine efficacy could help refine the obtained estimates. Second, we assumed that immunity induced either from infection or vaccination lasts more than the time horizon considered in the simulations (i.e., 1 year). There is evidence from laboratory studies and the field suggesting that the protection lasts several months [66]. Despite these preliminary pieces of evidence, the duration of the immunity remains a research area of paramount importance and intrinsically linked to viral evolution. It is also possible that waning immunity will continue to provide protection against severe disease but only partially against infection or transmission, which affects the herd immunity threshold. Overall, the duration and quality of immunity will determine the periodicity of COVID-19 outbreaks globally [67, 68]. Moreover, booster vaccination may be an efficient way to improve the vaccine effectiveness [69–71]. For example, in Chile’s report about the effectiveness of booster dose [69–71], the vaccine effectiveness of CoronaVac against infection increases from 50.18 to 70.89% after booster shot. The increased effectiveness of vaccination associated with the booster shot may contribute to increased immunity in the population and deserves further investigation.
Third, in the baseline scenario, we referred to an inactivated SARS-CoV-2 vaccine (BBIBP-CorV) taken to be 54.3% efficacious against the Delta variant infection. However, several other vaccines (including CoronaVac, WBIP-CorV, Ad5-nCoV, and ZF2001) are licensed and have been used in China. We varied vaccine efficacy up to 79% in sensitivity analyses. The main conclusion about the potential of herd immunity and the need to extend the vaccination campaign to children early as well as to use more efficacious vaccines is unaltered.
Conclusion
In conclusion, based on the current evidence, reaching vaccine-induced herd immunity in a population with little/no natural immunity is challenging. A key step has been made in early November 2021 with the authorization of a vaccine for 3+ years old individuals. Minimizing vaccine hesitancy in all age groups will be another key step to increase the immunity level of the population. These, together with highly efficacious vaccines or booster vaccinations, will be even more crucial given the possible emergence of new SARS-CoV-2 variants that are more transmissible or with immune escape. Importantly, even if herd immunity is unlikely to be reached due to the waning of immunity and the emergence of new variants, vaccination will continue to dramatically reduce the COVID-19 burden.
References
1. Wu JT, Leung K, Leung GM. Nowcasting and forecasting the potential domestic and international spread of the 2019-nCoV outbreak originating in Wuhan, China: a modeling study. Lancet. 2020;395(10225):689–97. https:// doi.org/10.1016/S0140-6736(20)30260-9.
2. Pan A, Liu L, Wang C, Guo H, Hao X, Wang Q, et al. Association of public health interventions with the epidemiology of the COVID-19 outbreak in Wuhan. China. JAMA. 2020;323(19):1915–23. https://doi.org/10.1001/jama.2 020.6130.
3. Li Q, Guan X, Wu P, Wang X, Zhou L, Tong Y, et al. Early transmission dynamics in Wuhan, China, of novel coronavirus-infected pneumonia. N Engl J Med. 2020;382(13):1199–207. https://doi.org/10.1056/NEJMoa2001316.
4. Chinazzi M, Davis JT, Ajelli M, Gioannini C, Litvinova M, Merler S, et al. The effect of travel restrictions on the spread of the 2019 novel coronavirus (COVID-19) outbreak. Science. 2020;368(6489):395–400. https://doi.org/1 0.1126/science.aba9757.
5. World Health Organization. WHO coronavirus (COVID-19) dashboard. 2021. https://covid19.who.int/. Accessed 12 November 2021.
6. Washington NL, Gangavarapu K, Zeller M, Bolze A, Cirulli ET, Schiabor Barrett KM, et al. Emergence and rapid transmission of SARS-CoV-2 B.1.1.7 in the United States. Cell. 2021;184(10):2587–94.e2587. https://doi.org/10.1016/j. cell.2021.03.052.
7. Gaymard A, Bosetti P, Feri A, Destras G, Enouf V, Andronico A, et al. Early assessment of diffusion and possible expansion of SARS-CoV-2 Lineage 20I/ 501Y.V1 (B.1.1.7, variant of concern 202012/01) in France, January to March 2021. Euro Surveill. 2021;26(9):2100133. https://doi.org/10.2807/1560-7917. ES.2021.26.9.2100133.
8. Davies NG, Abbott S, Barnard RC, Jarvis CI, Kucharski AJ, Munday JD, et al. Estimated transmissibility and impact of SARS-CoV-2 lineage B.1.1.7 in England. Science. 2021;372(6538):eabg3055. https://doi.org/10.1126/ science.abg3055.
9. Leung K, Shum MH, Leung GM, Lam TT, Wu JT. Early transmissibility assessment of the N501Y mutant strains of SARS-CoV-2 in the United Kingdom, October to November 2020. Euro Surveill. 2021;26(1). https://doi. org/10.2807/1560-7917.ES.2020.26.1.2002106.
10. Zhao S, Lou J, Cao L, Zheng H, Chong MKC, Chen Z, et al. Quantifying the transmission advantage associated with N501Y substitution of SARS-CoV-2 in the UK: an early data-driven analysis. J Travel Med. 2021;28(2). https://doi. org/10.1093/jtm/taab011.
11. Graham MS, Sudre CH, May A, Antonelli M, Murray B, Varsavsky T, et al. Changes in symptomatology, reinfection, and transmissibility associated with the SARS-CoV-2 variant B.1.1.7: an ecological study. Lancet Public Health. 2021;6(5):e335–45. https://doi.org/10.1016/S2468-2667(21)00055-4.
12. Volz E, Mishra S, Chand M, Barrett JC, Johnson R, Geidelberg L, et al. Assessing transmissibility of SARS-CoV-2 lineage B.1.1.7 in England. Nature. 2021;593(7858):266–9. https://doi.org/10.1038/s41586-021-03470-x.
13. Hoffmann M, Arora P, Groß R, Seidel A, Hörnich BF, Hahn AS, et al. SARS CoV-2 variants B.1.351 and P.1 escape from neutralizing antibodies. Cell. 2021;184(9):2384–93.e2312. https://doi.org/10.1016/j.cell.2021.03.036.
14. Wibmer CK, Ayres F, Hermanus T, Madzivhandila M, Kgagudi P, Oosthuysen B, et al. SARS-CoV-2 501Y.V2 escapes neutralization by South African COVID-19 donor plasma. Nat Med. 2021;27(4):622–5. https://doi.org/10.1038/s41 591-021-01285-x.
15. Faria NR, Mellan TA, Whittaker C, Claro IM, DdS C, Mishra S, et al. Genomics and epidemiology of the P.1 SARS-CoV-2 lineage in Manaus, Brazil. Science. 2021:eabh2644. https://doi.org/10.1126/science.abh2644.
16. Campbell F, Archer B, Laurenson-Schafer H, Jinnai Y, Konings F, Batra N, et al. Increased transmissibility and global spread of SARS-CoV-2 variants of concern as of June 2021. Euro Surveill. 2021;26(24). https://doi.org/10.2807/1 560-7917.Es.2021.26.24.2100509.
17. Challen R, Dyson L, Overton CE, Guzman-Rincon LM, Hill EM, Stage HB, et al. Early epidemiological signatures of novel SARS-CoV-2 variants: establishment of B.1.617.2 in England. medRxiv. 2021;2021(21258365):2006– 5. https://doi.org/10.1101/2021.06.05.21258365.
18. Dagpunar J. Interim estimates of increased transmissibility, growth rate, and reproduction number of the COVID-19 B.1.617.2 variant of concern in the United Kingdom. medRxiv. 2021;2021(21258293):2006–3. https://doi.org/1 0.1101/2021.06.03.21258293.
19. Outbreak.info. Variant of Concern Reports. https://outbreak.info/situation reports. Accessed August 10, 2021.
20. Chinese Center for Disease Control and Prevention. Distribution of novel coronavirus disease 2019. http://2019ncov.chinacdc.cn/2019-nCoV/. Accessed November 12, 2021.
21. Zhang M, Xiao J, Deng A, Zhang Y, Zhuang Y, Hu T, et al. Transmission dynamics of an outbreak of the COVID-19 delta variant B.1.617.2 – Guangdong Province, China, May – June 2021. China CDC Weekly. 2021; 3(27):584–6. https://doi.org/10.46234/ccdcw2021.148.
22. Wang YP, Chen RC, Hu FY, Lan Y, Yang ZW, Zhan C, et al. Transmission, viral kinetics and clinical characteristics of the emergent SARS-CoV-2 Delta VOC in Guangzhou, China. Clinical medicine. 2021;40. https://doi.org/10.1016/j. eclinm.2021.101129.
23. Lopez Bernal J, Andrews N, Gower C, Gallagher E, Simmons R, Thelwall S, et al. Effectiveness of Covid-19 vaccines against the B.1.617.2 (Delta) Variant. N Engl J Med. 2021. https://doi.org/10.1056/NEJMoa2108891.
24. National Health Commission of the People’s Republic of China. Update on the doses of COVID-19 vaccine administration. http://www.nhc.gov.cn/xcs/ yqjzqk/list_gzbd.shtml. Accessed November 12, 2021.
25. Britton T, Ball F, Trapman P. A mathematical model reveals the influence of population heterogeneity on herd immunity to SARS-CoV-2. Science. 2020; 369(6505):846–9. https://doi.org/10.1126/science.abc6810.
26. May RM, Anderson RM. Spatial heterogeneity and the design of immunization programs. Math Biosci. 1984;72(1):83–111. https://doi.org/10.1 016/0025-5564(84)90063-4.
27. Zhang J, Klepac P, Read JM, Rosello A, Wang X, Lai S, et al. Patterns of human social contact and contact with animals in Shanghai, China. Sci Rep. 2019;9(1):15141. https://doi.org/10.1038/s41598-019-51609-8.
28. Hu S, Wang W, Wang Y, Litvinova M, Luo K, Ren L, et al. Infectivity, susceptibility, and risk factors associated with SARS-CoV-2 transmission under intensive contact tracing in Hunan, China. Nat Commun. 2021;12(1): 1533. https://doi.org/10.1038/s41467-021-21710-6.
29. Xin H, Li Y, Wu P, Li Z, Lau EHY, Qin Y, et al. Estimating the latent period of coronavirus disease 2019 (COVID-19). Clin Infect Dis. 2021. https://doi.org/1 0.1093/cid/ciab746.
30. Zhao S, Tang B, Musa SS, Ma S, Zhang J, Zeng M, et al. Estimating the generation interval and inferring the latent period of COVID-19 from the contact tracing data. Epidemics. 2021;36:100482. https://doi.org/10.1016/j. epidem.2021.100482.
31. Sun K, Wang W, Gao L, Wang Y, Luo K, Ren L, et al. Transmission heterogeneities, kinetics, and controllability of SARS-CoV-2. Science. 2021; 371(6526). https://doi.org/10.1126/science.abe2424.
32. Pang XH, Ren LL, Wu SS, Ma WT, Yang J, Di L, et al. Cold-chain food contamination is the possible origin of the COVID-19 resurgence in Beijing. Natl Sci Rev. 2020;7(12):1861–4. https://doi.org/10.1093/nsr/nwaa264.
33. Chen X, Chen Z, Azman AS, Deng X, Sun R, Zhao Z, et al. Serological evidence of human infection with SARS-CoV-2: a systematic review and meta-analysis. Lancet Glob Health. 2021;9(5):e598–609. https://doi.org/10.101 6/S2214-109X(21)00026-7.
34. National Bureau of Statistics. China Population & Employment Statistics Yearbook 2019. https://navi.cnki.net/knavi/yearbooks/YZGRL/detail. Accessed March 8, 2021.
35. United Nations. World Population Prospects 2019. https://population.un.org/ wpp/Download/Standard/Population/. Accessed March 8, 2021.
36. Wang HH, Wang JJ, Wong SY, Wong MC, Li FJ, Wang PX, et al. Epidemiology of multimorbidity in China and implications for the healthcare system: a cross-sectional survey among 162,464 community household residents in southern China. BMC Med. 2014;12(1):188. https:// doi.org/10.1186/s12916-014-0188-0.
37. World Health Organization. WHO SAGE values framework for the allocation and prioritization of COVID-19 vaccination. https://www.who.int/publications/i/item/who-sage-values-framework-for-the-allocation-and-prioritization of-covid-19-vaccination. Accessed March 8, 2021.
38. World Health Organization. Interim recommendations for use of the Pfizer– BioNTech COVID-19 vaccine, BNT162b2, under Emergency Use Listing. https://www.who.int/publications/i/item/WHO-2019-nCoV-vaccines-SAGE_ recommendation-BNT162b2-2021.1. Accessed March 8, 2021.
39. World Health Organization. Pfizer BioNTech COVID-19 vaccine: What you need to know. https://www.who.int/news-room/feature-stories/detail/who-can-take the-pfizer-biontech-covid-19%2D%2Dvaccine. Accessed March 8, 2021.
40. World Health Organization. Interim recommendations for use of the ChAdOx1-S [recombinant] vaccine against COVID-19 (AstraZeneca COVID-19 vaccine AZD1222, SII Covishield, SK Bioscience). https://www.who.int/ publications/i/item/WHO-2019-nCoV-vaccines-SAGE_recommendation-A ZD1222-2021.1. Accessed March 8, 2021.
41. Chen X, Chen Z, Azman AS, Sun R, Lu W, Zheng N, et al. Neutralizing antibodies against SARS-CoV-2 variants induced by natural infection or vaccination: a systematic review and pooled meta-analysis. Clin Infect Dis. 2021. https://doi.org/10.1093/cid/ciab646.
42. Al Kaabi N, Zhang Y, Xia S, Yang Y, Al Qahtani MM, Abdulrazzaq N, et al. Effect of 2 inactivated SARS-CoV-2 vaccines on symptomatic COVID-19 infection in adults: a randomized clinical trial. JAMA. 2021;326(1):35–45. https://doi.org/10.1001/jama.2021.8565.
43. Khoury DS, Cromer D, Reynaldi A, Schlub TE, Wheatley AK, Juno JA, et al. Neutralizing antibody levels are highly predictive of immune protection from symptomatic SARS-CoV-2 infection. Nat Med. 2021;27(7):1205–11. https://doi.org/10.1038/s41591-021-01377-8.
44. Palacios R, Batista AP, Albuquerque CSN, Patiño EG, JdP S, MTRP C, et al. Efficacy and safety of a COVID-19 inactivated vaccine in healthcare professionals in Brazil: The PROFISCOV Study. SSRN. 2021. https://doi.org/1 0.2139/ssrn.3822780.
45. Li XN, Huang Y, Wang W, Jing QL, Zhang CH, Qin PZ, et al. Effectiveness of inactivated SARS-CoV-2 vaccines against the Delta variant infection in Guangzhou: a test-negative case-control real-world study. Emerg Microbes Infect. 2021;10(1):1751–9. https://doi.org/10.1080/22221751.2021.1969291.
46. Ma K, Yi Y, Li Y, Sun L, Deng A, Hu T, et al. Effectiveness of inactivated COVID-19 vaccines against COVID-19 pneumonia and severe illness caused by the B.1.617.2 (Delta) variant: evidence from an outbreak in Guangdong, China. SSRN. 2021. https://doi.org/10.2139/ssrn.3895639.
47. Hu Z, Tao B, Li Z, Song Y, Yi C, Li J, et al. Effectiveness of inactive COVID-19 vaccines against severe illness in B.1.617.2 (Delta) variant-infected patients in Jiangsu, China. medRxiv. 2021;2021(21263010):2009–2. https://doi.org/10.11 01/2021.09.02.21263010.
48. Sheikh A, McMenamin J, Taylor B, Robertson C. SARS-CoV-2 Delta VOC in Scotland: demographics, risk of hospital admission, and vaccine effectiveness. Lancet. 2021;397(10293):2461–2. https://doi.org/10.1016/S014 0-6736(21)01358-1.
49. Halloran ME, Longini IM, Struchiner CJ. Modes of action and time-varying VES. In: Design and Analysis of Vaccine Studies. New York, NY: Springer New York; 2010. p. 131–51.
50. Yang J, Chen X, Deng X, Chen Z, Gong H, Yan H, et al. Disease burden and clinical severity of the first pandemic wave of COVID-19 in Wuhan, China. Nat Commun. 2020;11(1):5411. https://doi.org/10.1038/s41467-020-19238-2.
51. Poletti P, Tirani M, Cereda D, Trentini F, Guzzetta G, Sabatino G, et al. Association of age with the likelihood of developing symptoms and critical disease among close contacts exposed to patients with confirmed SARS CoV-2 infection in Italy. JAMA Netw Open. 2021;4(3):e211085. https://doi. org/10.1001/jamanetworkopen.2021.1085.
52. Fisman DN, Tuite AR. Progressive increase in virulence of novel SARS-CoV-2 variants in Ontario, Canada. medRxiv. 2021;2021(21260050):2007–5. https:// doi.org/10.1101/2021.07.05.21260050.
53. Diekmann O, Heesterbeek JA, Metz JA. On the definition and the computation of the basic reproduction ratio R0 in models for infectious diseases in heterogeneous populations. J Math Biol. 1990;28(4):365–82. https://doi.org/10.1007/BF00178324.
54. Yang J, Marziano V, Deng X, Guzzetta G, Zhang J, Trentini F, et al. Despite vaccination, China needs non-pharmaceutical interventions to prevent widespread outbreaks of COVID-19 in 2021. Nat Hum Behav. 2021;5(8):1009– 20. https://doi.org/10.1038/s41562-021-01155-z.
55. Mistry D, Litvinova M, Pastore YPA, Chinazzi M, Fumanelli L, Gomes MFC, et al. Inferring high-resolution human mixing patterns for disease modeling. Nat Commun. 2021;12(1):323. https://doi.org/10.1038/s41467-020-20544-y.
56. Wang J, Jing R, Lai X, Zhang H, Lyu Y, Knoll MD, et al. Acceptance of COVID-19 Vaccination during the COVID-19 pandemic in China. Vaccines (Basel). 2020;8(3). https://doi.org/10.3390/vaccines8030482.
57. Chen M, Li Y, Chen J, Wen Z, Feng F, Zou H, et al. An online survey of the attitude and willingness of Chinese adults to receive COVID-19 vaccination. Hum Vaccin Immunother. 2021;17(7):1–10. https://doi.org/10.1080/2164551 5.2020.1853449.
58. Wang C, Han B, Zhao T, Liu H, Liu B, Chen L, et al. Vaccination willingness, vaccine hesitancy, and estimated coverage at the first round of COVID-19 vaccination in China: a national cross-sectional study. Vaccine. 2021;39(21): 2833–42. https://doi.org/10.1016/j.vaccine.2021.04.020.
59. Gan L, Chen Y, Hu P, Wu D, Zhu Y, Tan J, et al. Willingness to receive SARS CoV-2 vaccination and associated factors among Chinese adults: a cross-sectional survey. Int J Environ Res Public Health. 2021;18(4). https://doi.org/1 0.3390/ijerph18041993.
60. Lazarus JV, Ratzan SC, Palayew A, Gostin LO, Larson HJ, Rabin K, et al. A global survey of potential acceptance of a COVID-19 vaccine. Nat Med. 2021;27(2):225–8. https://doi.org/10.1038/s41591-020-1124-9.
61. Hodgson D, Flasche S, Jit M, Kucharski AJ, Group CC-W. Centre for Mathematical Modelling of Infectious Disease C-WG. The potential for vaccination-induced herd immunity against the SARS-CoV-2 B.1.1.7 variant. Euro Surveill. 2021;26(20). https://doi.org/10.2807/1560-7917.ES.2021.26.20.2100428.
62. Kwok KO, Lai F, Wei WI, Wong SYS, Tang JWT. Herd immunity - estimating the level required to halt the COVID-19 epidemics in affected countries. J Inf. 2020;80(6):e32–3. https://doi.org/10.1016/j.jinf.2020.03.027.
63. Omer SB, Yildirim I, Forman HP. Herd immunity and implications for SARS CoV-2 control. JAMA. 2020;324(20):2095–6. https://doi.org/10.1001/jama.202 0.20892.
64. Centers for Disease Control and Prevention. Nationwide Commercial Laboratory Seroprevalence Survey. https://covid.cdc.gov/covid-data-tracker/ #national-lab. Accessed August 8, 2021.
65. Zhang Y, Zeng G, Pan H, Li C, Hu Y, Chu K, et al. Safety, tolerability, and immunogenicity of an inactivated SARS-CoV-2 vaccine in healthy adults aged 18-59 years: a randomized, double-blind, placebo-controlled, phase 1/ 2 clinical trial. Lancet Infect Dis. 2021;21(2):181–92. https://doi.org/10.1016/ S1473-3099(20)30843-4.
66. Wang Z, Muecksch F, Schaefer-Babajew D, Finkin S, Viant C, Gaebler C, et al. Naturally enhanced neutralizing breadth against SARS-CoV-2 one year after infection. Nature. 2021;595(7867):426–31. https://doi.org/10.1038/s41586- 021-03696-9.
67. Lavine JS, Bjornstad ON, Antia R. Immunological characteristics govern the transition of COVID-19 to endemicity. Science. 2021;371(6530):741. https:// doi.org/10.1126/science.abe6522.
68. Saad-Roy CM, Wagner CE, Baker RE, Morris SE, Farrar J, Graham AL, et al. Immune life history, vaccination, and the dynamics of SARS-CoV-2 over the next 5 years. Science. 2020;370(6518):811–8. https://doi.org/10.1126/ science.abd7343.
69. Barda N, Dagan N, Cohen C, Hernan MA, Lipsitch M, Kohane IS, et al. Effectiveness of a third dose of the BNT162b2 mRNA COVID-19 vaccine for preventing severe outcomes in Israel: an observational study. Lancet. 2021; 398(10316):2093–100. https://doi.org/10.1016/S0140-6736(21)02249-2.
70. Saciuk Y, Kertes J, Shamir Stein N, Ekka ZA. Effectiveness of a third dose of BNT162b2 mRNA vaccine. J Infect Dis. 2021;225(1):30–3. https://doi.org/10.1 093/infdis/jiab556.
71. Rafael Araos AJ, vCovid-Ministry of Health. COVID-19 vaccine effectiveness assessment in Chile. https://cdn.who.int/media/docs/default-source/blue print/chile_rafael-araos_who-vr-call_25oct2021.pdf?sfvrsn=7a7ca72a_7. Accessed October 10, 2021.
