Abstract
Immune waning is key to the timely anticipation of COVID-19 long-term dynamics. We assess the impact of periodic vaccination campaigns using a compartmental epidemiological model with multiple age structures and parameterised using empiric time-dependent vaccine protection data. Despite the uncertainty inherent to such scenarios, we show that vaccination campaigns decreases the yearly number of COVID-19 admissions. However, especially if restricted to individuals over 60 years old, vaccination on its own seems insufficient to prevent thousands of hospital admissions and it suffers the comparison with non-pharmaceutical interventions aimed at decreasing infection transmission. The combination of such interventions and vaccination campaigns appear to provide the greatest reduction in hospital admissions.
1 Introduction
From the beginning, Covid-19 pandemic management had to deal with numerous unknowns and strongly relied on mathematical modelling to guide non-pharmaceutical interventions (NPIs) implementation. The rapid discovery of efficient vaccines led to the hope that public health policy planning could soon return to normal. However, long-term forecasts are difficult because vaccine protection decreases with time [UKHSA, 2022]. Furthermore, the emergence of variants of concerns (VOC) raises major concerns, especially because of lineages from the Omicron variant (Pango lineage BA) that exhibit strong immune escape properties.
In the French context, Bosetti et alii, 2022, studied under what booster administration and NPIs implementation a new epidemic could be contained. However, their work was done before the Omicron VOC, and only on a six-month horizon. Kissler et alii, 2020 investigated the potential long-term effects of the seasonality observed in other epidemics than that of SARS-CoV-2. They highlighted the need for NPIs, but their results were obtained before vaccine implementation and the emergence of VOCs. Saad-Roy et alii, 2021 studied different vaccine administration patterns for several scenarios of immunity duration but, their long-term insights were very uncertain because at the time, there was no data on immune protection waning. Finally, Ghosh et alii [2022] used a non-Markovian setting to capture immunity waning but remain on a few months horizon.
The waning of immune protection is challenging to capture with classical compartmental models because it implies an exponential increase in the number of model compartments, which makes precise parameterization challenging. We extend a non-Markovian approach [Reyné et alii, 2022] that readily accounts for the time spent in each compartment. By using published data, we explore long-term epidemic dynamics in a qualitative way. In our scenarios, we account for immunity waning as well as Omicron-specific phenotypic traits and compare four vaccination campaign strategies. Furthermore, we also investigate the impact of non-pharmaceutical interventions (NPIs), such as mask-wearing or air quality improvement, that can decrease the transmission rate of the infection.
2 Methods
An epidemiological model with time structures
We use an epidemiological compartmental model (Figure 1A) in which susceptible individuals of age a (the density of which is denoted Sa) can either become fully vaccinated (Va), or contract a mild or severe infection . Mildly-infected individuals always recover and move to the compartment Ra. Both recovered and vaccinated individuals can be (re)infected, but at a reduced rate compared to susceptible individuals. If this (re)infection is mild, individuals move to a separate compartment to account for a potential immunity-induced reduction in infectiousness. Vaccinated and recovered individuals may be (re)vaccinated and move to the booster compartment (Ba), where their protection increases. Finally, severely infected and previously immunised mildly infected individuals also end up in the booster compartment upon recovery, assuming high protection against potential new infection. Overall, the boosted compartment consists of individuals with booster vaccination dose(s), two natural infections, one vaccination and one infection, or having recovered from a severe infection.
The model accounts for memory effects, meaning that we record the time spent by the individuals in each compartment. Knowing the time since vaccination (k) for vaccinated individuals, the time since clearance (j) for recovered individuals, and the time since the entry into the booster compartment (l), allows us to readily account for the waning in immune protection.
This model is based on a system of partial differential equations (Appendix S2). Its parameterization reflects the French epidemic and the Omicron VOC (Table 1) [Sofonea et alii, 2022].
Contact rates were allowed to vary with time and fitted using the hospital admissions time series up to May 6, 2022. Following a parsimony principle, we used the last fitted value as the baseline until the rest of the simulations.
External factors, such as the weather, are known to impact transmission dynamics [Ma et alii, 2021]. We included this seasonality by assuming sinusoidal variations such that in summer the contact rate is decreased to –10% and increased by +10% in the winter.
Immune waning
In contrast to earlier models, we could calibrate immunity waning using epidemiological data UKHSA [2022]. More precisely, we used time series of vaccine protection (whether initial doses or boosters) against both symptomatic disease and hospitalisation, for both Delta and Omicron VOC. We fitted linear models on the different time series (Figure 1B).
Vaccination campaign scenarios
Vaccinated individuals were assumed to receive a booster dose 6 months after entering the vaccinated (Va) compartment. Then, we investigate four different scenarios:
In Scenario A, boosted individuals are not vaccinated again.
Scenario B consists in implementing annual vaccination campaigns before winter (in September and October) but only for individuals above 60 years old.
We extend the yearly vaccination campaign to all the population in Scenario C.
Scenario D assumes two national vaccine campaigns per year (in September – October and in March – April) for all the population.
To better assess the impact of vaccination campaigns, we also investigated two additional scenarios in which we implemented a decreased contact rate of -20% compared to the May 2022 value used in the baseline scenario. This value was chosen because it is comparable to others observed in 2021 [Reyné et alii, 2022].
3 Results
Scenario A leads to a high level of daily hospital admissions (Figure 2A) with yearly oscillations attributable to seasonality. Scenario B improves the overall situation but the median number of hospital admissions always remains above 500 per day (Figure 2B). Vaccinating everyone once a year (Scenario C) further lowers the number of daily hospital admissions and it also yields an epidemic wave in the early spring (Figure 2C). Finally, the fourth scenario (Figure 2D).lowers the number of daily hospital admissions even further but it does not prevent two marked epidemic waves per year.
Both Scenarios C and D have the more pronounced epidemic waves. This can be explained by the fact that vaccinating everyone at the same time implies that immunity wane for everyone at the same time.
Figure 2E shows the annual total number of hospital admissions for each scenario. As expected, the more the people are vaccinated, the more the daily hospital admissions lower. However, vaccination alone just seems to contain what was reached with stringent NPIs in 2021 in France (curfews, lockdown, health pass).
Note that, as shown also in Supplementary Figures S1 to S4, simulations yield large 95% confidence intervals (CI) for the total number of hospital admissions.
Finally, we also explored the impact of a 20% reductiong in the contact rate baseline for scenarios A and C (Supplementary Figures S5-S6). On its own the reduction in contact rate can have a stronger impact than yearly vaccination campaigns for the whole population (scenario A’ in Figure 2E). Combining yearly vaccination and a decrease in contact rate provides the strongest decrease in the total number of yearly hospital admissions.
4 Discussion
COVID-19 management now faces the challenge of post-infection and post-vaccination immune waning combined with the emergence risk of new VOCs.
Long term scenarios are obviously extremely uncertain but they can provide valuable insights. By leaning on Omicron-based data on immune protection, we show that immune waning may cause yearly epidemic hospitalization peaks comparable to the largest one seen in 2021. Although vaccination clearly mitigates this impact, it seems insufficient to suppress the epidemic. As shown in the Supplementary Figures S5-S6, variations in these contact rates, i.e. the intensity of NPIs or behavourial changes, has a strong impact on hospitalizations. This is also shown in Figure 2E where NPIs achieving a 20% decrease in contact rate can lead to a decrease in hospitalisations comparable to yearly vaccination of the whole population. Furthermore, we find that the combined use of vaccination along NPIs should be strongly considered to limit the impact of COVID-19 on the hospital burden (Figure 2E).
These results should be taken with caution and regarded as qualitatively prospective due to the numerous sources of uncertainty (Appendix S1). In particular, comparisons should be restricted to our different scenarios that share the same core assumptions.
The sensitivity analyses (Figures S7-S10) highlight the main sources of uncertainty which correspond to factors that are difficult to predict. As discussed by Reyné et alii, 2022, the time-varying contact matrix is impossible to predict — as it depends on government policies, age-specific spontaneous behavioural changes or calendar events such as school holidays— and yet yield an huge uncertainty in the model’s outputs. The seasonality also impacts strongly the results, and also unpredictable. Due to computational constraints, the sensitivity analysis did not contained the baseline contact rate even if it has a strong impact (Figure 2E). Regarding virus-related model parameters, the reduction in contagiousness due to immunity, which is difficult to estimate [Bosetti et alii, 2022; Prunas et alii, 2022], contributes to most of the variance in the model output.
Some factors are not included in the model but could affect the dynamics. For instance, this study neglects virus evolution on a long-term scale although five VOCs have already spread in France in 2021. It also neglects the potential hospitalizations attributable to patients with long-COVID.
Finally, we assumed that the intensity of non-pharmaceutical interventions will remain identical to that enforced in early spring 2022. This neglects any changes in government policies, some of which would probably be necessary to avoid hospitals saturation for some parameter sets.
This work based on empiric estimates of immune protection waning underlines the importance of combining vaccination with other type of interventions, especially NPIs such as improving indoor air quality or mask-wearing. Future work is required to identify the optimal schedule for COVID-19 vaccination campaigns as it will require narrowing many unknowns regarding the biology and spread of the virus.
Data Availability
The scripts are available on demand.
Statements and Declarations
Funding
BR is funded by the Ministère de l’Enseignement Supérieur et de la Recherche.
Competing Interests
The authors have no relevant financial or non-financial interests to disclose.
Authors’ Contributions
All authors contributed to the study conception and design, interpreted the results and wrote the manuscript. Analysis were performed by BR. All authors read and approved the final manuscript.
Data availability
Data and scripts will be made available in a public repository upon acceptance.
S1 Supplementary figures
S2 Model equations
The model partial differential equations system is given by: with for any (t, a, i, j, k) ∈ ℝ+ × [0, amax] × [0, imax] × [0, jmax] × [0, kmax].
The parameter notations are the following:
ρ is the initial vaccination rate,
ρV is the first booster dose vaccination campaign,
ρB is the re-vaccination rate of prospective vaccination campaigns,
γm,s the recovery rates for respectively mildly and severely infected individuals,
µ is the death rate,
εR,V,B the immunity-induced reduction of risk of infection for individuals respectively in the R, V and B compartments,
c the intensity of NPI policy,
K(a, a′) the contact matrix coefficient between age groups a and a′,
βm,s,d,mv the generation time distributions, and
ξ the immunity-induced reduction in transmission (for ‘breakthrough’ infections).
The previous system is coupled with the following boundary conditions: where
pa is the probability of developing a severe form,
ifra is the infection fatality rate, and
νV,R,B is the immunity-induced reduction of virulence for individuals respectively in the R, V, and B compartments.
S3 Omicron related parameters and vaccine efficacy
Regarding the omicron generation interval, we used the data provided by UKSHA [2022]. In particular, we fitted different Gamma distributions on the non-parametric data available. We tested different parameters combinations for the Gamma distribution to explore a range of generation time distributions that reflect the epidemiology of both BA.1 and BA.2 variants. In Figure S11, we show a subset of these Gamma distributions explored within the sensitivity analysis.
Regarding the parameterisation of vaccine effectiveness, we used data provided by the UKHSA [2022] report (https://assets.publishing.service.gov.uk/government/uploads/system/uploads/attachment_data/file/1070356/Vaccine-surveillance-report-week-16.pdf). The raw data was not available so we used the online tool WebPlot Digitalizer (https://apps.automeris.io/wpd/) to retrieve the values. For simplicity, we assumed that individuals received the Pfizer/BioNTech (BNT162b2) vaccine, which was the most widespread in France. We fitted linear models on vaccine effectiveness against symptomatic disease and hospital admissions for the Delta and Omicron VOC. Note that the intrinsic virulence of Omicron is assumed to be divided by 3 (parameter pa in the model) compared to Delta following UK data Nyberg et alii, 2022.
In the sensitivity analysis, we included variations in vaccine protection against symptomatic disease and against hospitalization risk. We used the linear model best fit as our baseline and generated variation by applying a coefficient to the intercept spanning from − 0.05 to 0.05 (Figure 1B).
The reduction of transmission (in so-called ‘breakthrough’ infections) was assumed to be of 50% for vaccinated individuals, as in others modelling works [Bosetti et alii, 2022]. This reduction in transmission was also applied to recovered people that got reinfected.
For the other model parameters, they are identical to Reyné et alii [2022].
Acknowledgements
The authors acknowledge the ISO 9001 certified IRD i-Trop HPC (South Green Platform) at IRD Montpellier for providing HPC resources that have contributed to the research results reported within this paper. We thank all the ETE modelling team for discussions and apologize to Baptiste Elie for that reckless <monospace>rm -r *</monospace>.
Footnotes
Reference added. Equations corrected in Appendix S2. Lines numbered.