Relative contributions of vaccination and previous infection to population-level SARS-CoV-2 immunity over time: a simulation modelling study ============================================================================================================================================ * Joshua Szanyi * Tim Wilson * Hassan Andrabi * Tony Blakely ## Abstract Population-level immunity to SARS-CoV-2 directly impacts the incidence of COVID-19 morbidity and mortality. Understanding how this immunity is likely to change over time in the context of future vaccination schedules and emerging SARS-CoV-2 variants is critical to inform pandemic policy. This study simulates population-level COVID-19 immunity (including relative contributions of vaccination and previous infection) in Victoria, Australia over 18 months using an agent-based model and logistic regression equations that predict immunity and waning following vaccination and/or infection. Previous infection was found to drive most immunity against infection even with ongoing regular vaccination, however a greater proportion of overall immunity against mortality was accounted for by vaccination. Although previous infection appears to be driving a substantial component of population-level COVID-19 immunity currently, improved vaccines providing longer lasting (and better sterilizing) immunity are likely to be a critical component of the future pandemic response given the risks associated with SARS-CoV-2 infection. Keywords * COVID-19 * vaccine * immunity * infection * model ## Introduction The emergence of immune-evasive variants, waning immunity, and relatively poor sterilizing immunity afforded by vaccination mean that currently available vaccines have had less of an impact on the long-term course of the COVID-19 pandemic than hoped. Protection against symptomatic Omicron BA.1 infection following a third dose of an mRNA vaccine, for example, peaks at approximately 60% – 70% and wanes rapidly thereafter; 1 protection against more severe clinical outcomes is higher but still decreases over time. 2 The Omicron sub-variants BA.4/BA.5 possess even greater immune evasion capacity than earlier sub-variants, 3 highlighting the challenge of responding to an ever-evolving virus under strong selection pressure to evade pre-existing immune responses in the population. 4 Where the primary infection and subsequent re-exposure are due to the same SARS-CoV-2 variant, the level of immunity that develops following infection appears to be high, particularly among younger people. 5 It also appears to wane at a slower rate than vaccine-derived immunity. 5 Hybrid immunity (vaccination augmented by natural infection) provides additional protection above either sensitizing event alone.6 COVID-19-related morbidity and mortality is influenced by the level of immunity in a population at any given point in time. Accordingly, it is important to understand total population-level immunity, the relative contributions of vaccination and previous infection to this immunity, and how this is likely to change to inform ongoing pandemic policy such as vaccine scheduling. This is a challenge, however, given the current dynamic and mixed picture of vaccination and infection-derived immunity in many jurisdictions in addition to imperfect knowledge regarding expected protection against future variants. We previously developed a logistic regression model of vaccine effectiveness (VE) that provides estimates of protection and waning against symptomatic infection and hospitalisation with the Omicron BA.1 variant. 7 We subsequently extended this equation, incorporating it into an agent-based model (ABM) to assess optimal pandemic policy in the context of ongoing SARS-CoV-2 evolution over an 18-month period in Victoria, Australia. 8 The current study aims to use this ABM to simulate total population-level immune protection against SARS-CoV-2 in Victoria over the same period, and the relative contributions of vaccination and previous infection to this protection given four illustrative COVID-19 vaccination schedules and eight potential SARS-CoV-2 variant scenarios. ## Methods VE estimates were based on a previously published logistic regression model that estimated waning protection against symptomatic infection and hospitalisation from the Omicron BA.1 variant. 7 We extended this initial equation to account for age using data from a study reporting age-stratified VE against the Delta variant. 9 We then added the outcomes of death and any infection (using UK Health Security Agency data). 2 Finally, we developed a similar model to represent protection following previous infection (that is, protection against reinfection with the same SARS-CoV-2 variant upon re-exposure) using the same relative differences in protection across clinical outcomes as for VE and waning set halfway between nil waning and that for VE against hospitalisation. 7 Hybrid immunity was calculated as 1 – (1 – VE) × (1 – protection from natural infection) based on evidence suggesting that vaccine- and infection-derived immunity act independently to provide an overall level of protection. 6 For further details, please see the Supplementary Materials of 8. These immunity estimates were then assigned to agents in an ABM that was developed to assess optimal pandemic policy in the state of Victoria, Australia (6.6 million people), in the 18 months from 1 April 2022.8 For the current study, four vaccine schedules were included – nil further vaccination beyond 1 April 2022; three 6-monthly boosters of current-generation mRNA vaccines (with rollout from July 2022, January 2023 and July 2023); two doses of an illustrative next-generation Omicron-targeted vaccine (with rollout from October 2022 and April 2023); and two doses of an Omicron-targeted vaccine with its waning rate set at 25% of that for current vaccines (Table 1). Relative VE of the Omicron-targeted vaccine against the SARS-CoV-2 variants modelled in the current study was set using an odds ratio of 2 applied to current-generation VE estimates (Supplementary Table 3; by way of illustration, a doubling of VE on the OR scale means a VE of 90% becomes 95%, and 40% become 57%). View this table: [Table 1:](http://medrxiv.org/content/early/2022/08/19/2022.08.18.22278963/T1) Table 1: Modelled vaccination schedules At each round of vaccination, 80% of agents vaccinated in the previous round were vaccinated with the new vaccine. Vaccine schedules were modelled in the context of eight potential future SARS-CoV-2 variants, all with incursion into Victoria 91 days into the model run, with two potential levels of innate transmissibility (R0 11 or 14), immune escape capacity and virulence (low, approximating the Omicron variant, or high) 8. Immune escape possibilities for each SARS-CoV-2 variant are shown in Supplementary Table 3 and Supplementary Table 4. The virus circulating at model outset, prior to the emergence of a new variant, approximated Omicron BA.2. In the ABM, public health and social measures (including physical distancing and mask wearing) (de)escalated through five stages based on expected health service pressure (Supplementary Table 1 and Supplementary Table 2). Over the 18-month period modelled, weekly SARS-CoV-2 infections and deaths due to COVID-19 were averaged across 100 separate draws of model input parameters, each run four times to dampen stochastic variation. To illustrate variability across the 100 model runs, 5th and 95th percentiles of weekly infections and deaths are also presented in the results; the purpose of this study was not to precisely predict epidemic dynamics in the short-term but explore the cumulative impact of various pandemic policies given a highly uncertain future. Population-level immunity was quantified as the average percentage reduction (compared to nil vaccination and nil previous infection) of the probability of the event of interest (being infected, death if infected) across all agents. This average immunity was calculated from vaccines alone, previous infection alone, and the total from both vaccination and previous infection (including hybrid immunity). ## Results Figure 1 shows total population immunity over 18 months for all modelled viral variants, in the context of 6-monthly boosters with a current generation mRNA vaccine. Total immunity varied depending on the variant, particularly for protection against infection. Lower levels of immunity against infection were seen with more highly virulent variants; higher levels of total immunity were seen with more transmissible variants possessing lower virulence. ![Figure 1:](http://medrxiv.org/https://www.medrxiv.org/content/medrxiv/early/2022/08/19/2022.08.18.22278963/F1.medium.gif) [Figure 1:](http://medrxiv.org/content/early/2022/08/19/2022.08.18.22278963/F1) Figure 1: Simulated total immunity against (a) any infection and (b) death for 6-monthly boosters with a current generation vaccine across eight modelled SARS-CoV-2 variants*. * Variant incursion at 13 weeks. See Table 1 for vaccination timing. Irrespective of the vaccination scenario, population-level immunity reached a relatively steady state, but with more peaks and troughs in scenarios with ongoing vaccination. This is evident in Figure 2 and Figure 3, which show simulated total immunity, vaccine-derived immunity, and immunity from previous infection against any infection (symptomatic and asymptomatic infection) and death for one future ![](http://medrxiv.org/https://www.medrxiv.org/content/medrxiv/early/2022/08/19/2022.08.18.22278963/F2/graphic-4.medium.gif) [](http://medrxiv.org/content/early/2022/08/19/2022.08.18.22278963/F2/graphic-4) ![](http://medrxiv.org/https://www.medrxiv.org/content/medrxiv/early/2022/08/19/2022.08.18.22278963/F2/graphic-5.medium.gif) [](http://medrxiv.org/content/early/2022/08/19/2022.08.18.22278963/F2/graphic-5) Figure 2: Simulated total immunity, vaccine-derived immunity, and immunity from previous infection against any infection over an 18-month period in Victoria, Australia for one future SARS-CoV-2 variant* and (a) nil further vaccination, (b) 6-monthly boosters with a current generation vaccine, (c) two doses of an Omicron-targeted next-generation vaccine with the same waning as current vaccines, and (d) two doses of an Omicron-targeted vaccine with 75% less waning than current vaccines. * variant with R 11, antigenically Omicron-like, virulence approximating the Omicron variant, emerging in July 2022 (13 weeks), with high immune escape capacity (halving vaccine effectiveness and/or natural protection on the odds ratio scale). See Table 1 for vaccination timing. CG: current generation; OT: Omicron-targeted. ![](http://medrxiv.org/https://www.medrxiv.org/content/medrxiv/early/2022/08/19/2022.08.18.22278963/F3/graphic-6.medium.gif) [](http://medrxiv.org/content/early/2022/08/19/2022.08.18.22278963/F3/graphic-6) ![](http://medrxiv.org/https://www.medrxiv.org/content/medrxiv/early/2022/08/19/2022.08.18.22278963/F3/graphic-7.medium.gif) [](http://medrxiv.org/content/early/2022/08/19/2022.08.18.22278963/F3/graphic-7) Figure 3: Simulated total immunity, vaccine-derived immunity, and immunity from previous infection against death over an 18-month period in Victoria, Australia for one future SARS-CoV-2 variant* and (a) nil further vaccination, (b) 6-monthly boosters with a current generation vaccine, (c) two doses of an Omicron-targeted next-generation vaccine with the same waning as current vaccines, and (d) two doses of an Omicron-targeted vaccine with 75% less waning than current vaccines. * variant with R 11, antigenically Omicron-like, virulence approximating the Omicron variant, emerging in July 2022 (13 weeks), with high immune escape capacity (halving vaccine effectiveness and/or natural protection on the odds ratio scale). See Table 1 for vaccination timing. CG: current generation; OT: Omicron-targeted. SARS-CoV-2 variant (R0 11, virulence approximating that of the Omicron variant, antigenically similar to the Omicron variant, with high immune escape capacity that halves VE or natural protection on the OR scale) across the four vaccination scenarios. With respect to total infections (Figure 2), previous infection drove most total immunity even with ongoing regular vaccination. Vaccination events provided boosts to total population immunity against infection which waned rapidly. Modelled Omicron-targeted vaccines resulted in higher peaks of total population-level immunity compared to current-generation vaccines. Waning of vaccine-derived immunity was accompanied by waves of infections followed by boosts in natural-infection derived immunity, keeping total immunity to infection relatively steady. When no further vaccines were provided, vaccine-derived immunity waned to negligible levels within months and was compensated for by natural immunity to stabilize total population-level immunity against infection. Mean cumulative infections over 18 months were 6.5 million, 6.3 million, 6 million, and 5.6 million for the nil, current generation, Omicron-targeted, and Omicron-targeted (25% waning) vaccine scenarios respectively. For deaths (Figure 3), a greater proportion of overall immunity at the population level was accounted for by vaccination compared to that for infection (noting that waning of VE against deaths is less than that against infection). Like that seen for protection against infection, as vaccine-derived immunity to death fell, natural immunity compensated to keep overall population immunity levels relatively stable. In the scenario where no further vaccines were given beyond 1 April 2022, total immunity against death fell to 67% by 18 months; this was lower than in all scenarios of ongoing vaccination (81% at 18 months for the scenario of an Omicron-targeted vaccine with reduced waning, for example). Mean cumulative deaths over 18 months were 5,670, 3,655, 3,497, and 3,205 for the nil further, current generation, Omicron-targeted, and Omicron-targeted (25% waning) vaccine scenarios respectively. ## Discussion Forecasting future changes in population-level immunity to SARS-CoV-2 is a complex and inherently uncertain task, made challenging by factors such as the development of new vaccines, the impact of ongoing infection, the effects of waning immunity, and the emergence of novel variants. Such analyses are critical, however, to inform decision-making. This study illustrates how population immunity is likely to vary over an 18-month period in the setting of a new SARS-CoV-2 variant and several illustrative vaccination schedules. Our model predicts that population-level immunity will vary with SARS-CoV-2 variant characteristics. In terms of infection risk, less protection was observed at the aggregate level when the newly emergent viral variant was more virulent than current Omicron sub-variants. This reflects that for highly virulent variants hospital capacity rapidly exceeds a tolerable limit, resulting in the imposition of more stringent public health and social measures to reduce transmission. While crucial for avoiding acute morbidity and mortality and pressure on health care services from a highly virulent variant, this affords the population less protection against transmission in the future. Issues such as these highlight the complex trade-offs that characterise pandemic policy making. For the illustrative SARS-CoV-2 variant shown in Figure 2 and Figure 3, a steady state of immunity is reached no matter the vaccination policy modelled. For protection against infection, the contribution of vaccination to overall protection in the longer-term is minimal. This differs for protection against mortality, where vaccination is responsible for the bulk of immunity early in the modelled period and provides ongoing high levels of protection in scenarios where regular vaccination continues. The marginal benefit of vaccination was observed to decrease over time, as only 80% of those who were vaccinated in the previous round receive the new vaccine (which we believe is a reasonable, if not even conservative, assumption). In Australia, 96% of eligible people received two COVID-19 vaccine doses while only 72% of those eligible have received a third dose10. In situations of no further vaccination, the fall in vaccination-derived immunity is partially compensated for by natural infection immunity. This, however, is accompanied by a steadier and more sustained pattern of infections and deaths compared to the wave-like patterns for infection and mortality seen with ongoing regular vaccination. Accordingly, our modelling supports the assertion that vaccination is the preferred method of achieving maximal population immunity to SARS-CoV-2 given the risks associated with SARS-CoV-2 infection, including acute morbidity, mortality, long COVID and additional sequelae including increased cardiovascular disease risk. 11 Vaccination, by contrast, is associated with a very low risk of adverse effects and has the additional benefit of reducing the risk of developing long COVID. 12 Predictably, modelled vaccines with improved VE and less rapid waning were able to achieve similar levels of immunity to current generation vaccines but with less frequent dosing. ## Conclusions This simulation modelling study suggests that in the absence of a new COVID-19 vaccine which provides substantially better sterilizing immunity than currently available products, the bulk of protection at the population level against SARS-CoV-2 infection is likely to continue to come from natural infection. However, vaccination appears to be responsible for a substantial amount of population-level protection against serious outcomes such as mortality. Given the risks associated with SARS-CoV-2 infection, the development of improved vaccines with enhanced protection against infection and longer lasting immunity is likely to be a critical component of the future response to the pandemic. Increasing the proportion of individuals vaccinated with each rollout would also likely assist in maintaining vaccine-mediated protection at maximum achievable levels. Finally, there is an urgent need for better estimates of protection following primary SARS-CoV-2 infection, and how these change over time and between variants. In combination with VE and waning estimates of next-generation vaccines, this will assist in informing ongoing policy decisions regarding optimal vaccination scheduling in Australia and internationally. ## Supporting information Supplementary material [[supplements/278963_file03.pdf]](pending:yes) ## Data Availability Requests for additional model output data may be granted upon reasonable request to the researchers. Access to model code is generally not available but negotiated access may be possible. Please contact the researchers for further information. ## Funding An anonymous philanthropic donation and strategic University of Melbourne funding. ## Acknowledgements The authors would like to acknowledge Shania Rossiter, Samantha Howe and Jessie Zeng for their contributions during the development of the model upon which this study is based and during the conceptualization of this manuscript. * Received August 18, 2022. * Revision received August 18, 2022. * Accepted August 19, 2022. * © 2022, Posted by Cold Spring Harbor Laboratory This pre-print is available under a Creative Commons License (Attribution-NonCommercial-NoDerivs 4.0 International), CC BY-NC-ND 4.0, as described at [http://creativecommons.org/licenses/by-nc-nd/4.0/](http://creativecommons.org/licenses/by-nc-nd/4.0/) ## References 1. 1.Andrews N, Stowe J, Kirsebom F, et al. Covid-19 vaccine effectiveness against the omicron (B. 1.1. 529) variant. New England Journal of Medicine 2022. 2. 2.UK Health Security Agency. COVID-19 vaccine surveillance report - week 15, 2022. 3. 3.Hachmann NP, Miller J, Collier A-rY, et al. Neutralization Escape by SARS-CoV-2 Omicron Subvariants BA.2.12.1, BA.4, and BA.5. New England Journal of Medicine 2022; 387(1): 86–8. 4. 4.Markov PV, Katzourakis A, Stilianakis NI. Antigenic evolution will lead to new SARS-CoV-2 variants with unpredictable severity. Nature Reviews Microbiology 2022; 20(5): 251–2. [PubMed](http://medrxiv.org/lookup/external-ref?access_num=http://www.n&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2022%2F08%2F19%2F2022.08.18.22278963.atom) 5. 5.Hansen CH, Michlmayr D, Gubbels SM, Mølbak K, Ethelberg S. Assessment of protection against reinfection with SARS-CoV-2 among 4 million PCR-tested individuals in Denmark in 2020: a population-level observational study. Lancet 2021; 397(10280): 1204–12. [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1016/S0140-6736(21)00575-4&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=33743221&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2022%2F08%2F19%2F2022.08.18.22278963.atom) 6. 6.Altarawneh HN, Chemaitelly H, Ayoub HH, et al. Effects of Previous Infection and Vaccination on Symptomatic Omicron Infections. New England Journal of Medicine 2022. 7. 7.Szanyi J, Wilson T, Scott N, Blakely T. A log-odds system for waning and boosting of COVID-19 vaccine effectiveness. Vaccine 2022; 40(28): 3821–4. 8. 8.Szanyi J, Wilson T, Howe S, Zeng J, Andrabi H, Blakely T. An integrated epidemiologic and economic model to assess optimal COVID-19 pandemic policy. medRxiv 2022: 2022.08.01.22278262. 9. 9.Andrews N, Tessier E, Stowe J, et al. Duration of Protection against Mild and Severe Disease by Covid-19 Vaccines. New England Journal of Medicine 2022. 10. 10.Australian Government Department of Health and Aged Care. Vaccination numbers and statistics. 2022. [https://www.health.gov.au/initiatives-and-programs/covid-19-vaccines/numbers-statistics](https://www.health.gov.au/initiatives-and-programs/covid-19-vaccines/numbers-statistics) (accessed August 17 2022). 11. 11.Xie Y, Xu E, Bowe B, Al-Aly Z. Long-term cardiovascular outcomes of COVID-19. Nature Medicine 2022; 28(3): 583–90. 12. 12.Antonelli M, Penfold RS, Merino J, et al. Risk factors and disease profile of post-vaccination SARS-CoV-2 infection in UK users of the COVID Symptom Study app: a prospective, community-based, nested, case-control study. The Lancet Infectious Diseases 2022; 22(1): 43–55. [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1016/S1473-3099(21)00460-6&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=34480857&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2022%2F08%2F19%2F2022.08.18.22278963.atom)