Modelling long-term COVID-19 hospital admission dynamics using empirical immune protection waning data ====================================================================================================== * Bastien Reyné * Mircea T. Sofonea * Samuel Alizon ## 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 embedded multiple age structures and empiric time-dependent vaccine protection kinetics. 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 [Kucharski et alii, 2020]. The first steps of the COVID-19 management due to the SARS-CoV-2 emergence led to exploring some possibilities. For instance, some considered the possibility to reach herd immunity — thus terminating the epidemic provided a sufficient proportion of the population is immune to the disease — as a perennial replacement of the emergency use of lock-downs and other strong non-pharmaceutical interventions (NPIs) to alleviate the hospital burden. These considerations happened at a time when vaccines were not available. Given the estimates of the proportion of the population needed to be immune to reach the said threshold, it would have required a massive infection of the population [Kwok et alii, 2020; Randolph and Barreiro, 2020; Fontanet and Cauchemez, 2020], and hence a non-negligible death toll in different populations. In late 2020, the availability of safe and efficient vaccines provided a powerful tool for pandemic management. Furthermore, it was realised that in addition to protecting against severe symptoms, vaccination also protected against infection itself and even reduced transmission rate leading to the hope that we could eradicate, or at least easily control, the virus. However, the hope was a little dampened by the emergence of the Alpha variant of concern (VOC) that was both more transmissible and more virulent than the ancestral lineages [Davies, Abbott, et alii, 2021; Davies, Jarvis, et alii, 2021]. And it also reminded virus evolution cannot be taken out of the equation [Alizon and Sofonea, 2021]. The emergence of the Delta VOC a few months later — once again both more transmissible and more virulent [Fisman and Tuite, 2021; Alizon, Haim-Boukobza, et alii, 2021] — combined with the vaccine efficacy *ab initio* cast some doubt on the strategy to be adopted. Later on, the overall SARS-CoV-2 reported time-induced loss of immunity entombed the possibility for herd immunity. In addition to this emergence of variant issue, field data rapidly showed that immunity against SARS-CoV-2 wanes over time. For instance, protection against Omicron VOC severe forms decreases from ∼90% to ∼73% in six months for individuals with a full vaccination schedule (2 doses with Pfizer for 65+ y.o. individuals) [UKHSA, 2022]. For the protection against infection, the decrease is even more pronounced [UKHSA, 2022]. This is consistent with numerous studies on antibody titres that show a decline with time [Peng et alii, 2022; Pérez-Alós et alii, 2022; Levin et alii, 2021]. This problem is exacerbated by the fact that some VOCs are particularly prone to evading pre-existing immune responses [Dejnirattisai et alii, 2022]. Note that individual histories become particularly important to anticipate the efficiency of the immune response [Reynolds et alii, 2022]. In France, the current management strategy consists in vaccinating the population as much as possible and implementing NPIs when epidemic waves when hospital capacity approaches its limits. Mathematical modelling has provided in many countries and all along the pandemic useful insights and good approximations to assess what could happen in the near future if the current conditions remain the same. For instance, the first waves and the emergence of the Alpha, Delta and Omicron VOCs were quite rapidly followed by modelling studies providing useful insights to guide public policies [Salje et alii, 2020; Domenico et alii, 2021; Alizon, Haim-Boukobza, et alii, 2021; Sofonea, Roquebert, et alii, 2022]. However, past one month, the projections made become less precise and past few months some mathematical and biological assumptions might rapidly alter the general trends [Sofonea and Alizon, 2021]. That being said, we believe such long-term trends could still be of interest to the authorities that would consider the possibility of long-term management. Another proof-of-concept point would be to correct some limitations of many mathematical models based on the usual set of ordinary differential equations (ODEs) [Sofonea, Reyné, et alii, 2021]. Thanks to the work and dedication of some people, e.g. at UKHSA, this pandemic generates data that is unique in the history of epidemics, such as the vaccine efficacy over time since vaccination for the different variants circulating in the population. Such data can help us anticipate long-term trends. However, usual ODEs-based models are inadequate to integrate empirical data that depart too much from the exponential distribution, which is the case for immune protection waning. Because of these technical difficulties, long-term projections are scarce in the literature. 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. However, 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 their scenarios also remained on a few months’ horizon. In this work, by extending a non-Markovian approach [Reyné et alii, 2022] that readily accounts for the time spent in each compartment, we leverage the field data on immune protection from UKHSA and qualitatively explore long-term epidemic dynamics. We focus on the case of France but it reflects that of many high-income countries in terms of vaccination coverage and age distributions. In our scenarios, we account for immunity waning as well as Omicron-specific phenotypic traits and compare three vaccination campaign strategies. Furthermore, we also investigate the impact of non-pharmaceutical interventions (NPIs) at the whole population level, 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 build on an existing non-Markovian epidemiological compartmental model that accounts for time since entry in the different compartments [Reyné et alii, 2022]. The main extension is that we implement vaccine booster doses and loss of immunity over time for vaccinated, boosted, or recovered individuals (Figure 1A). Susceptible individuals of age *a* (the density of which is denoted *S**a*) can either become fully vaccinated (*V**a*), or contract a mild ![Graphic][1] or severe infection ![Graphic][2] Mildly-infected individuals always recover and move to the compartment *R**a*. Both recovered and vaccinated individuals can be (re)infected, but at a reduced rate compared to susceptible individuals (see also the section below about parameterisation for further details about the biological assumptions made). If this (re)infection is mild, individuals move to a separate compartment ![Graphic][3] to account for a potential immunity-induced reduction in infectiousness. Vaccinated and recovered individuals may be (re)vaccinated and move to the booster compartment (*B**a*), where their protection increases. Finally, severely infected ![Graphic][4] and previously immunised mildly infected individuals ![Graphic][5] 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. ![Figure 1:](http://medrxiv.org/https://www.medrxiv.org/content/medrxiv/early/2022/08/11/2022.05.16.22275130/F1.medium.gif) [Figure 1:](http://medrxiv.org/content/early/2022/08/11/2022.05.16.22275130/F1) Figure 1: Description of the memory-based model with immune waning. **A)** Flowchart of the compartmental model. Arrows show transitions between compartments. *S**a* stands for Susceptible individuals of age ![Graphic][6] and ![Graphic][7] stand respectively for mildly/severely/severely-that-will-die/mildly-partly-immune infected individuals of age *a* infected since *i* days, *V**ak* stands for Vaccinated individuals of age *a* vaccinated *k* days ago, *R**aj* stands for Recovered individuals of age *a* that cleared the disease *j* days ago and finally *B**aℓ* stands for individuals of age *a* that received a booster vaccine dose *ℓ* days ago. Orange arrows show some of the transitions that depend on the time spent in the compartment (here *ℓ* days) that are parametrized through real immune waning data shown on panel B. **B)** Waning of vaccine efficacy against infection. Dots corresponds to real data from UKHSA [2022] for the Pfizer/BioNTech vaccine (BNT162b2) after a booster dose. The lines correspond to the baseline of the immunity decrease model implementation and the shaded areas to the uncertainty used within the sensitivity analysis. Here, we show the protection against an Omicron/BA.1 VOC infection for individuals that received a booster dose. Large grey arrows indicate where the waning data shown acts in the compartmental model. 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 (*ℓ*), allows us to readily account for the waning in immune protection. This particularity is modelled using partial differential equations, that provide for each compartment the rate at which individuals enter the compartment with a boundary equation and the time-dependent depart of the compartment through the differential equation. For instance, for the vaccinated compartment we have the boundary equation ![Formula][8] coupled with the differential equation ![Formula][9] where *ρ*(*t, a*) is the initial vaccination rate (1st and 2nd doses), *ε**V* (*k*) is the vaccine efficacy *k* days after the first dose, *λ*(*t, a*) is the force of infection and *ρ**V* (*k*) corresponds to the third dose vaccination rate. The complete set of equations is available in Appendix S2. ### Fitting procedure and transmission rates To simulate long-term scenarios, we need to generate a population with a heterogeneous exposition background to SARS-CoV-2 infection and vaccination. This is done by fitting the dynamics of the daily COVID-19 hospital admissions, which correspond to the severe cases in our model. The fit procedure is based on ordinary least squares and is done piecewise (the length of pieces varies, depending on important changes in transmission rates). Transmission rates were fitted from the beginning of the simulations (Jan 1, 2021) up to May 6, 2022. The last fitted value corresponds to the baseline transmission rate at the end of the simulations. This is not completely neutral since the fitting account for NPIs or spontaneous behavioural changes in the reflection of the epidemic. In France, at that time the restrictions were quite low, only mask-wearing in public transportations was required so it reflects a situation with very scarce epidemic control amongst the population. External factors, such as the weather, are known to impact transmission dynamics up to -17% [Ma et alii, 2021]. We decided to include a seasonality by assuming sinusoidal variations such that in summer the transmission rate is decreased to –10% and increased by +10% in the winter. Since this seasonality is largely driven by behavioural changes, we assume that it applies similarly to all variants. Finally, in some scenarios, we further decrease the baseline transmission rate by 20%. This is done to assess both the potential impact of long-term NPIs (such as improving indoor air quality) and to test the sensitivity of the model regarding this parameter. ### Virus related parameters The presented model is not multi-strain, each time a new VOC becomes more prevalent it provokes a switch in the model parameters. We also neglect virus evolution meaning the current resident strain (Omicron BA.2 VOC) will remain till the end of simulations. The virus-related parameters are available in Table 1. View this table: [Table 1:](http://medrxiv.org/content/early/2022/08/11/2022.05.16.22275130/T1) Table 1: Description of the Omicron-related model parameters and their default value. For each parameter, we indicate the default value used, the range in the sensitivity analyses, and the references. ### Vaccinal efficacy and immunity waning In contrast to earlier models, we could calibrate immunity waning using epidemiological data (Figure 1B). The UKHSA reports provide time series for vaccine efficacy over time for the different variants UKHSA [2022]. They also provide estimates of vaccine effectiveness against symptomatic disease, which we used as protection against infection in our model (we assumed no difference between age groups for this parameter). We also used their vaccine effectiveness time series against hospitalisation to parameterize variations in protection against severe cases in our model. This waning was assumed to differ between individuals younger or older than 60 y.o. The implementation was done by assuming a linear model for the decrease in immune protection, as detailed in Appendix S4. Note also that although we had access to vaccine efficacy for both two-doses vaccinated and boosted individuals, we could not find similar data for recovered individuals. For simplicity purposes, we assumed that recovered individuals had the same level of immunity as two-doses vaccinated individuals. This hypothesis may overestimate the level of protection of natural infections against infection by another variant. We more-over assumed that previously hospitalised individuals had higher protection than individuals that developed mild symptoms, as it is reported they have higher antibody titers [Servellita et alii, 2022]. ### Past vaccine rollout and prospective future vaccination campaigns As in many high-income countries, in France, the first vaccination campaign started in early January 2021. In our model, the vaccine rollout was fitted by age groups accordingly to real data (see Reyné et alii [2022] for details). The youngest individuals (0-10 years old) were not eligible for vaccination in this model. For this initial vaccination campaign, we assumed moreover an automatic third (booster) dose six months after the second one, meaning that every individual leave the *V**a* compartment to enter the boosted compartment (*B**a*). The future (and prospective) vaccine campaigns — presented hereafter — concern only individuals in the booster compartment (*B**a*). We assumed perfect population compliance to future vaccination from the individuals already vaccinated, meaning that everyone eligible receives his booster vaccine dose. This assumption is likely to be optimistic but it is also counterbalanced by the fact that we also assumed unvaccinated individuals will remain as such. We investigated three 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**. ## 3 Results Scenario A leads to a high level of daily hospital admissions (Figure 2A) with yearly oscillations attributable to seasonality. ![Figure 2:](http://medrxiv.org/https://www.medrxiv.org/content/medrxiv/early/2022/08/11/2022.05.16.22275130/F2.medium.gif) [Figure 2:](http://medrxiv.org/content/early/2022/08/11/2022.05.16.22275130/F2) Figure 2: Median number of COVID-19 hospital admissions in four vaccination campaign scenarios. Blue dots corresponds to real data and green shaded areas to vaccination campaigns periods. The four panels correspond to scenarios without additional booster **(A)**, with an annual vaccination campaign in September-October for every individuals above 60 y.o. **(B)**, an annual vaccination campaign in September-October for the whole population **(C)**. See Supplementary Figures S1 to S3 for the 95% confidence intervals and the interquartile range. 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). Scenario C exhibits the most pronounced epidemic waves, which appear to occur twice a year. This can be explained by the fact that vaccinating everyone at the same time implies that immunity wane for everyone at the same time. This generates two periodicities, one related to seasonality, and the second related to vaccination campaigns. Figure 3 shows the annual total number of hospital admissions for each scenario. As expected, the more 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). ![Figure 3:](http://medrxiv.org/https://www.medrxiv.org/content/medrxiv/early/2022/08/11/2022.05.16.22275130/F3.medium.gif) [Figure 3:](http://medrxiv.org/content/early/2022/08/11/2022.05.16.22275130/F3) Figure 3: Comparison of the total number of hospital admissions per year for four vaccination scenarios. Blue line corresponds to real data for 2021. The different type of boxplots correspond to the different scenarios. Scenarios **A’, B’** and **C’**, in the gray shaded area, corresponds to Scenarios A, B and C where the baseline transmission rate was decreased by 20%. Complete dynamics for those alternate scenarios are shown in Supplementary Figures S4-S5. Note that, as shown also in Supplementary Figures S1 to S3, simulations yield large 95% confidence intervals (CI) for the total number of hospital admissions. Finally, we explored the impact of a 20% reduction in the baseline transmission rate for the three scenarios (Supplementary Figures S4-S5). On its own, such a reduction in transmission rates 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 transmission rate provides the strongest decrease in the total number of yearly hospital admissions. ## 4 Discussion COVID-19 management now faces two major challenges: first the post-infection and post-vaccination immune waning and, second, the emergence risk of new VOCs. These two factors amplify the level of uncertainty and unpredictability inherent to long-term predictions. However, more than two years after the first outbreak, five, six, or seven epidemic waves later depending on the countries and six million deaths world-wide, we believe it necessary to better apprehend long-term dynamics. It would provide valuable insights for health policy decision-making in comparing different strategies. Although the nature of future VOCs is largely unpredictable, mathematical modelling can incorporate the heterogeneity in immune protection observed in vaccinated populations by using alternate formalism that accounts for the time spent in each compartment and thus implements a realistic loss of immunity over time. We apply this approach to the case of France, which resembles that of many high-income countries, by leveraging empirical Omicron-based data on immune protection. Despite a huge uncertainty inherent to this approach, we obtain long-term qualitative COVID-19 hospital admissions trajectories with noticeable differences between scenarios. In most of our scenarios, we find that immune waning may cause yearly epidemic hospitalisation peaks comparable to the largest one seen in 2021. Although, as expected, vaccination decreases the total number of hospital admissions, our simulations suggest they are insufficient to suppress the epidemic and the public health burden by itself. By comparing scenarios with different levels of baseline NPIs, we show that all other things being equal, policies decreasing the transmission rate by 20% can lead to a decrease in hospitalisations comparable to yearly vaccination of the whole population (Figure 3). In comparison, mask-wearing is reported to have a similar order of magnitude, with a reported reproduction number reduction of -19% [Leech et alii, 2022]. We also show an effect in that the combined implementation of vaccination and NPIs have the greatest impact on reducing the hospital burden (Figure 3). These results should be taken with caution and regarded as qualitatively prospective due to the numerous sources of uncertainty (Appendix S1). Moreover, this study presents limitations induced by the strong biological assumptions summed up in Box 1. For these reasons, comparisons should be restricted to our different scenarios that share the same core assumptions. The sensitivity analyses (Figures S6-S8) 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 yields a huge uncertainty in the model’s outputs but is one of the most difficult model components to parameterize —as it depends on government policies, age-specific spontaneous behavioural changes or calendar events such as school holidays. Seasonality also impacts strongly the results and its precise effect is still under investigation. 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], does contribute to a large proportion of the variance in the model output (which is the daily number of newly hospitalised individuals). Some factors are not included in the model but could affect the dynamics. For instance, this study does not address virus evolution on a long-term scale although five VOCs have already spread in France in 2021. It also neglects the potential hospitalisations 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 hospital saturation for some parameter sets. Overall, in spite of its limitations, this work is, to our knowledge, one of the first to build on empiric estimates of immune protection waning to provide us with long-term perspectives. Its results underline the potential interest of combining vaccination with other types of interventions, especially NPIs such as improving indoor air quality or mask-wearing, to minimise the COVID-19 burden on hospitals in high-income countries. Future work could help to identify the optimal schedule for COVID-19 vaccination campaigns, which will require narrowing many unknowns regarding the biology and spread of the virus. Box 1: **Main biological assumptions, simplifications, and limitations of the model** As with any mathematical model, the one developed in this study makes several assumptions regarding the biological processes at work to be able to interpret the results in a meaningful manner. The main limitations are hereafter listed: * The **life histories of immunised hosts are limited to few patterns**. Indeed, we assume that all individuals in a given compartment have the same protection against the infection and the disease but, in reality, there is heterogeneity driven by the number of vaccine doses, infections, and their order [Reynolds et alii, 2022]. However, thanks to the non-Markovian structure of the model, we do capture temporal heterogeneity (i.e. protection waning). * We assume **a perfect compliance of population to vaccination**, thus ignoring a potential ‘fatigue’ in the population [Bodas et alii, 2022; Di Domenico et alii, 2021]. We also assume an **efficient vaccine rollout**, even though logistical issues might emerge. * There is **no virus evolution** in the model. However, on such a long-term scale, new phenotypic variants will likely emerge. * There is **no spatial structure**, which could artificially increase the magnitude of the yearly epidemic waves [Thomine et alii, 2021]. ## Data Availability The scripts will be available upon publication. ## 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 available upon publication in the following git repository [https://gitlab.in2p3.fr/ete/long\_term\_covid](https://gitlab.in2p3.fr/ete/long_term_covid). ## S1 Supplementary figures ![Figure S1:](http://medrxiv.org/https://www.medrxiv.org/content/medrxiv/early/2022/08/11/2022.05.16.22275130/F4.medium.gif) [Figure S1:](http://medrxiv.org/content/early/2022/08/11/2022.05.16.22275130/F4) Figure S1: Detailed output of Scenario A (no additional vaccination). Blue dots corresponds to real data. The black line correspond to the median trajectory. Lighter shaded area correspond to the 95% confidence interval while the darker area correspond to the interquartile range. The horizontal blue line indicates the highest national incidence in hospital admissions reached during the first COVID-19 wave in France. ![Figure S2:](http://medrxiv.org/https://www.medrxiv.org/content/medrxiv/early/2022/08/11/2022.05.16.22275130/F5.medium.gif) [Figure S2:](http://medrxiv.org/content/early/2022/08/11/2022.05.16.22275130/F5) Figure S2: Detailed output of Scenario B (yearly vaccination of individuals of more than 60 years old). Green shaded areas correspond to vaccination campaigns periods for individuals above 60 y.o. See Figure S1 for additional details. ![Figure S3:](http://medrxiv.org/https://www.medrxiv.org/content/medrxiv/early/2022/08/11/2022.05.16.22275130/F6.medium.gif) [Figure S3:](http://medrxiv.org/content/early/2022/08/11/2022.05.16.22275130/F6) Figure S3: Detailed output of Scenario C (yearly vaccination of all the population). See Figure S2 for additional details. ![Figure S4:](http://medrxiv.org/https://www.medrxiv.org/content/medrxiv/early/2022/08/11/2022.05.16.22275130/F7.medium.gif) [Figure S4:](http://medrxiv.org/content/early/2022/08/11/2022.05.16.22275130/F7) Figure S4: Detailed output of Scenario A (yearly vaccination of all the population) with a 20% reduction of the transmission rate. See Figure S2 for additional details. ![Figure S5:](http://medrxiv.org/https://www.medrxiv.org/content/medrxiv/early/2022/08/11/2022.05.16.22275130/F8.medium.gif) [Figure S5:](http://medrxiv.org/content/early/2022/08/11/2022.05.16.22275130/F8) Figure S5: Detailed output of Scenario C (yearly vaccination of all the population) with a 20% reduction of the transmission rate. See Figure S2 for additional details. ![Figure S6:](http://medrxiv.org/https://www.medrxiv.org/content/medrxiv/early/2022/08/11/2022.05.16.22275130/F9.medium.gif) [Figure S6:](http://medrxiv.org/content/early/2022/08/11/2022.05.16.22275130/F9) Figure S6: Global sensitivity analysis for Scenario A. The graphs shows the origin of the variance captured by each model parameter using Sobol indices (see Reyné et alii, 2022 for methodological details). A large part of the variations originates from the unknowns in the contact matrix between ages, as in Reyné et alii, 2022. The magnitude of seasonality also matters for the long-term trends, as well as the transmission rate (i.e. the intensity of the NPIs). ![Figure S7:](http://medrxiv.org/https://www.medrxiv.org/content/medrxiv/early/2022/08/11/2022.05.16.22275130/F10.medium.gif) [Figure S7:](http://medrxiv.org/content/early/2022/08/11/2022.05.16.22275130/F10) Figure S7: Global sensitivity analysis for Scenario B. See S6 for details. ![Figure S8:](http://medrxiv.org/https://www.medrxiv.org/content/medrxiv/early/2022/08/11/2022.05.16.22275130/F11.medium.gif) [Figure S8:](http://medrxiv.org/content/early/2022/08/11/2022.05.16.22275130/F11) Figure S8: Global sensitivity analysis for Scenario C. See S6 for details. ## S2 Model equations The model partial differential equations system is given by: ![Formula][10] ![Formula][11] ![Formula][12] ![Formula][13] ![Formula][14] ![Formula][15] ![Formula][16] ![Formula][17] with ![Formula][18] for any (*t, a, i, j, k*) ∈ ℝ+ × [0, *a*max] × [0, *i*max] × [0, *j*max] × [0, *k*max]. 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, * *β*(*t*) transmission rate, accounting for NPIs policies, * *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: ![Formula][19] ![Formula][20] ![Formula][21] ![Formula][22] ![Formula][23] ![Formula][24] ![Formula][25] with * *p**a* is the probability of developing a severe form, * ifr*a* 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 Model parameters For the other model parameters, they are identical to Reyné et alii [2022]. We sum up here the different parameters used in the model. However, we encourage readers to check Reyné et alii [2022] for the details regarding the contact matrixes, that will not be described in detail here. Long story short, we use 39 different France-specific contact matrixes corresponding to diverse NPIs implementation thath happened in France. For the others parameters, you can refer to the Table S1. View this table: [Table S1:](http://medrxiv.org/content/early/2022/08/11/2022.05.16.22275130/T2) Table S1: Model parameters. For each parameter, we indicate the default value used and the references used. ## S4 Omicron related parameters 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 parameter 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 S9, we show a subset of these Gamma distributions explored within the sensitivity analysis. The selected baseline is Gamma distribution with a 1.84 shape parameter and a 0.53 rate parameter. Note also that the intrinsic virulence of Omicron is assumed to be divided by 3 (parameter *p**a* in the model) compared to Delta following UK data. This is mostly qualitative and was not age-differentiated [Nyberg et alii, 2022]. Finally, note that we did not set up an intrinsic basic reproduction number (ℛ) for the Omicron VOC. While it was still possible for the Alpha VOC and arguably for the Delta VOC to have some estimates, the ℛ only refers to a completely naive population. Hence, when the transmission advantage of a VOC is partly due to an escape of the already built population immunity, it becomes hazardous to retrieve ‘the true ℛ’ for such a VOC. This is the case for the Omicron VOC. Therefore, we only fit the transmission rate on hospital admissions time series, rendering it difficult to dissociate the temporal reproduction number and the NPIs/changes in the behaviour of the population. ![Figure S9:](http://medrxiv.org/https://www.medrxiv.org/content/medrxiv/early/2022/08/11/2022.05.16.22275130/F12.medium.gif) [Figure S9:](http://medrxiv.org/content/early/2022/08/11/2022.05.16.22275130/F12) Figure S9: Omicron generation time. Shaded areas data originate from the UKSHA [2022] report. The lines correspond to a subset of the Gamma distributions used in the sensitivity analysis. The thicker line corresponds to the baseline Gamma distribution. ## S5 Vaccine efficacy The vaccine acts in three different ways in this model. First, the vaccine protects against infection. (This corresponds to the *ε**V/R/B* parameters in the model, and accounts for the time since the entry in the compartments.) 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](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/](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. On each time series (protection against Delta or Omicron, whether the individual received two-dose or a supplementary booster dose), we fitted a linear model (cf. Figure S10). Second, the vaccine protects also against severe forms if individuals are nonetheless infected. (In the model, the probability of developing a severe case *p**a* is decreased by (1−*ν**V/R/B*), which also accounts for the entry in the *V/R/B* compartments.) The same methodology is done as previsously, with the difference that UKHSA data show the total protection against virulence, *i*.*e*. [1−*ε**V/R/B*(·)] [1−*ν**V/R/B*(·)]. Another difference is that we distinguish between individuals above and below 60 y.o. as shown in Figure S11. One last thing to note is that we assumed that the last data point for individuals below 60 y.o. that received a booster was the upper limit provided and not the baseline value since the baseline for boosted individuals was below tho one with two-doses only. In the sensitivity analysis, we included variations in the total protection against virulence (cumulative effect of *ε* and *ν*). The variations correspond to add to the intercept of the linear model a coefficient between -0.05 and +0.05. Finally, the vaccine reduces the infectiousness of infected individuals. The reduction of transmission (in so-called ‘breakthrough’ infections) was assumed to be 50% for vaccinated individuals, as in other modelling works [Bosetti et alii, 2022]. This reduction in transmission was also applied to recovered people that got reinfected. Note that other studies tend to find different estimates [Prunas et alii, 2022]. ![Figure S10:](http://medrxiv.org/https://www.medrxiv.org/content/medrxiv/early/2022/08/11/2022.05.16.22275130/F13.medium.gif) [Figure S10:](http://medrxiv.org/content/early/2022/08/11/2022.05.16.22275130/F13) Figure S10: Model implementation of the decrease in immunity following booster vaccination. This is the vaccine efficacy against an Omicron infection after receiving a booster dose compared to fully susceptible individuals. Dots corresponds to real data from UKHSA [2022] for the Pfizer/BioNTech vaccine (BNT162b2) after a booster dose. The line correspond to the baseline of the immunity decrease model implementation. ![Figure S11:](http://medrxiv.org/https://www.medrxiv.org/content/medrxiv/early/2022/08/11/2022.05.16.22275130/F14.medium.gif) [Figure S11:](http://medrxiv.org/content/early/2022/08/11/2022.05.16.22275130/F14) Figure S11: Model implementation of the decrease in immunity following booster vaccination. This is the vaccine efficacy against an Omicron severe case after receiving a booster dose compared to fully susceptible individuals. Dots corresponds to real data from UKHSA [2022] for the Pfizer/BioNTech vaccine (BNT162b2) after a booster dose. Blue color is for individuals below 60 y.o. and red for older people. The lines correspond to the baseline of the immunity decrease model implementation. ## 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 rm -r *. ## Footnotes * The text has been expanded. The results didn't change that much. * Received May 16, 2022. * Revision received August 11, 2022. * Accepted August 11, 2022. * © 2022, Posted by Cold Spring Harbor Laboratory This pre-print is available under a Creative Commons License (Attribution-NonCommercial 4.0 International), CC BY-NC 4.0, as described at [http://creativecommons.org/licenses/by-nc/4.0/](http://creativecommons.org/licenses/by-nc/4.0/) ## References 1. Alizon, Samuel, Stéphanie Haim-Boukobza, et al.ii (July 2021). “Rapid spread of the SARS-CoV-2 Delta variant in some French regions, June 2021”. en. In: Eurosurveillance 26.28. Publisher: European Centre for Disease Prevention and Control, page 2100573. ISSN: 1560-7917. doi: 10.2807/1560-7917.ES.2021.26.28.2100573. URL: [https://www.eurosurveillance.org/content/10.2807/1560-7917.ES.2021.26.28.2100573](https://www.eurosurveillance.org/content/10.2807/1560-7917.ES.2021.26.28.2100573) (visited on 07/29/2021). [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.2807/1560-7917.ES.2021.26.28.2100573&link_type=DOI) 2. Alizon, Samuel and Mircea T. Sofonea (2021). “SARS-CoV-2 virulence evolution: Avirulence theory, immunity and trade-offs”. en. In: Journal of Evolutionary Biology n/a.n/a. ISSN: 1420-9101. doi: 10.1111/jeb.13896. URL: [https://onlinelibrary.wiley.com/doi/abs/10.1111/jeb.13896](https://onlinelibrary.wiley.com/doi/abs/10.1111/jeb.13896) (visited on 10/05/2021). [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1111/jeb.13896&link_type=DOI) 3. Bodas, Moran et al.ii (2022). “Overcoming the effect of pandemic fatigue on vaccine hesitancy—Will belief in science triumph?” en. In: Journal of Nursing Scholarship n/a.n/a (). ISSN: 1547-5069. doi: 10.1111/jnu.12778. URL: [https://onlinelibrary.wiley.com/doi/abs/10.1111/jnu.12778](https://onlinelibrary.wiley.com/doi/abs/10.1111/jnu.12778) (visited on 07/12/2022). [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1111/jnu.12778&link_type=DOI) 4. Bosetti, Paolo et al.ii (Jan. 2022). “Impact of booster vaccination on the control of COVID-19 Delta wave in the context of waning immunity: application to France in the winter 2021/22”. en. In: Eurosurveillance 27.1. Publisher: European Centre for Disease Prevention and Control, page 2101125. ISSN: 1560-7917. doi: 10.2807/1560-7917.ES.2022.27.1.2101125. URL: [https://www.eurosurveillance.org/content/10.2807/1560-7917.ES.2022.27.1.2101125](https://www.eurosurveillance.org/content/10.2807/1560-7917.ES.2022.27.1.2101125) (visited on 03/17/2022). [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.2807/1560-7917.ES.2022.27.1.2101125&link_type=DOI) 5. Davies, Nicholas G., Sam Abbott, et al.ii (Apr. 2021). “Estimated transmissibility and impact of SARS-CoV-2 lineage B.1.1.7 in England”. en. In: Science 372.6538. Publisher: American Association for the Advancement of Science Section: Research Article. ISSN: 0036-8075, 1095-9203. doi: 10.1126/science.abg3055. URL: [https://science.sciencemag.org/content/372/6538/eabg3055](https://science.sciencemag.org/content/372/6538/eabg3055) (visited on 05/31/2021). [Abstract/FREE Full Text](http://medrxiv.org/lookup/ijlink/YTozOntzOjQ6InBhdGgiO3M6MTQ6Ii9sb29rdXAvaWpsaW5rIjtzOjU6InF1ZXJ5IjthOjQ6e3M6ODoibGlua1R5cGUiO3M6NDoiQUJTVCI7czoxMToiam91cm5hbENvZGUiO3M6Mzoic2NpIjtzOjU6InJlc2lkIjtzOjE3OiIzNzIvNjUzOC9lYWJnMzA1NSI7czo0OiJhdG9tIjtzOjUwOiIvbWVkcnhpdi9lYXJseS8yMDIyLzA4LzExLzIwMjIuMDUuMTYuMjIyNzUxMzAuYXRvbSI7fXM6ODoiZnJhZ21lbnQiO3M6MDoiIjt9) 6. Davies, Nicholas G., Christopher I. Jarvis, et al.ii (May 2021). “Increased mortality in community-tested cases of SARS-CoV-2 lineage B.1.1.7”. en. In: Nature 593.7858. Number: 7858 Publisher: Nature Publishing Group, pages 270–274. ISSN: 1476-4687. doi: 10.1038/s41586-021-03426-1. URL: [https://www.nature.com/articles/s41586-021-03426-1](https://www.nature.com/articles/s41586-021-03426-1) (visited on 06/02/2021). [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1038/s41586-021-03426-1&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=33723411&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2022%2F08%2F11%2F2022.05.16.22275130.atom) 7. Dejnirattisai, Wanwisa et al.ii (Jan. 2022). “Reduced neutralisation of SARS-CoV-2 omicron B.1.1.529 variant by post-immunisation serum”. English. In: The Lancet 399.10321. Publisher: Elsevier, pages 234–236. ISSN: 0140-6736, 1474-547X. doi: 10.1016/S0140-6736(21)02844-0. URL: [https://www.thelancet.com/journals/lancet/article/PIIS0140-6736(21)02844-0/fulltext](https://www.thelancet.com/journals/lancet/article/PIIS0140-6736(21)02844-0/fulltext) (visited on 07/08/2022). [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1016/S0140-6736(21)02844-0&link_type=DOI) 8. Di Domenico, Laura et al.ii (Dec. 2021). “Adherence and sustainability of interventions informing optimal control against the COVID-19 pandemic”. en. In: Communications Medicine 1.1. Number: 1 Publisher: Nature Publishing Group, pages 1–13. ISSN: 2730-664X. doi: 10.1038/s43856-021-00057-5. URL: [https://www.nature.com/articles/s43856-021-00057-5](https://www.nature.com/articles/s43856-021-00057-5) (visited on 07/11/2022). [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1038/s43856-021-00057-5&link_type=DOI) 9. Domenico, Laura Di et alii (Apr. 2021). “Impact of January 2021 curfew measures on SARS-CoV-2 B.1.1.7 circulation in France”. en. In: Eurosurveillance 26.15. Publisher: European Centre for Disease Prevention and Control, page 2100272. ISSN: 1560-7917. doi: 10.2807/1560-7917.ES.2021.26.15.2100272. URL: [https://www.eurosurveillance.org/content/10.2807/1560-7917.ES.2021.26.15.2100272](https://www.eurosurveillance.org/content/10.2807/1560-7917.ES.2021.26.15.2100272) (visited on 07/08/2022). [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.2807/1560-7917.ES.2021.26.15.2100272&link_type=DOI) 10. Fisman, David N. and Ashleigh R. Tuite (Oct. 2021). “Evaluation of the relative virulence of novel SARS-CoV-2 variants: a retrospective cohort study in Ontario, Canada”. en. In: CMAJ 193.42. Publisher: CMAJ Section: Research, E1619–E1625. ISSN: 0820-3946, 1488-2329. doi: 10.1503/cmaj.211248. URL: [https://www.cmaj.ca/content/193/42/E1619](https://www.cmaj.ca/content/193/42/E1619) (visited on 11/03/2021). [Abstract/FREE Full Text](http://medrxiv.org/lookup/ijlink/YTozOntzOjQ6InBhdGgiO3M6MTQ6Ii9sb29rdXAvaWpsaW5rIjtzOjU6InF1ZXJ5IjthOjQ6e3M6ODoibGlua1R5cGUiO3M6NDoiQUJTVCI7czoxMToiam91cm5hbENvZGUiO3M6NDoiY21haiI7czo1OiJyZXNpZCI7czoxMjoiMTkzLzQyL0UxNjE5IjtzOjQ6ImF0b20iO3M6NTA6Ii9tZWRyeGl2L2Vhcmx5LzIwMjIvMDgvMTEvMjAyMi4wNS4xNi4yMjI3NTEzMC5hdG9tIjt9czo4OiJmcmFnbWVudCI7czowOiIiO30=) 11. Fontanet, Arnaud and Simon Cauchemez (Oct. 2020). “COVID-19 herd immunity: where are we?” en. In: Nature Reviews Immunology 20.10, pages 583–584. ISSN: 1474-1741. doi: 10.1038/s41577-020-00451-5. URL: [https://www.nature.com/articles/s41577-020-00451-5](https://www.nature.com/articles/s41577-020-00451-5) (visited on 02/08/2022). [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1038/s41577-020-00451-5&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=http://www.n&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2022%2F08%2F11%2F2022.05.16.22275130.atom) 12. Ghosh, Samiran, Malay Banerjee, and Vitaly Volpert (2022). “Immuno-Epidemiological Model-Based Prediction of Further Covid-19 Epidemic Outbreaks Due to Immunity Waning”. en. In: Mathematical Modelling of Natural Phenomena 17. Publisher: EDP Sciences, page 9. ISSN: 0973-5348, 1760-6101. doi: 10.1051/mmnp/2022017. URL: [https://www.mmnp-journal.org/articles/mmnp/abs/2022/01/mmnp220088/mmnp220088.html](https://www.mmnp-journal.org/articles/mmnp/abs/2022/01/mmnp220088/mmnp220088.html) (visited on 05/31/2022). [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1051/mmnp/2022017&link_type=DOI) 13. Kissler, Stephen M. et al.ii (May 2020). “Projecting the transmission dynamics of SARS-CoV-2 through the postpandemic period”. In: Science 368.6493. Publisher: American Association for the Advancement of Science, pages 860–868. doi: 10.1126/science.abb5793. URL: [https://www.science.org/doi/10.1126/science.abb5793](https://www.science.org/doi/10.1126/science.abb5793) (visited on 03/21/2022). [Abstract/FREE Full Text](http://medrxiv.org/lookup/ijlink/YTozOntzOjQ6InBhdGgiO3M6MTQ6Ii9sb29rdXAvaWpsaW5rIjtzOjU6InF1ZXJ5IjthOjQ6e3M6ODoibGlua1R5cGUiO3M6NDoiQUJTVCI7czoxMToiam91cm5hbENvZGUiO3M6Mzoic2NpIjtzOjU6InJlc2lkIjtzOjEyOiIzNjgvNjQ5My84NjAiO3M6NDoiYXRvbSI7czo1MDoiL21lZHJ4aXYvZWFybHkvMjAyMi8wOC8xMS8yMDIyLjA1LjE2LjIyMjc1MTMwLmF0b20iO31zOjg6ImZyYWdtZW50IjtzOjA6IiI7fQ==) 14. Kucharski, Adam J. et al.ii (May 2020). “Early dynamics of transmission and control of COVID-19: a mathematical modelling study”. English. In: The Lancet Infectious Diseases 20.5. Publisher: Elsevier, pages 553–558. ISSN: 1473-3099, 1474-4457. doi: 10.1016/S1473-3099(20)30144-4. URL: [https://www.thelancet.com/article/S1473-3099(20)30144-4/fulltext](https://www.thelancet.com/article/S1473-3099(20)30144-4/fulltext) (visited on 06/16/2022). [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1016/S1473-3099(20)30144-4&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=32171059&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2022%2F08%2F11%2F2022.05.16.22275130.atom) 15. Kwok, Kin On et al.ii (June 2020). “Herd immunity – estimating the level required to halt the COVID-19 epidemics in affected countries”. English. In: Journal of Infection 80.6. Publisher: Elsevier, e32–e33. ISSN: 0163-4453, 1532-2742. doi: 10.1016/j.jinf.2020.03.027. URL: [https://www.journalofinfection.com/article/S0163-4453(20)30154-7/fulltext](https://www.journalofinfection.com/article/S0163-4453(20)30154-7/fulltext) (visited on 02/08/2022). [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1016/j.jinf.2020.03.027&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=http://www.n&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2022%2F08%2F11%2F2022.05.16.22275130.atom) 16. Leech, Gavin et al.ii (2022). “Mask wearing in community settings reduces SARS-CoV-2 transmission”. In: Proceedingsofthe National Academyof Sciences 119.23, e2119266119. doi: 10.1073/pnas.2119266119. URL: [https://www.pnas.org/doi/abs/10.1073/pnas.2119266119](https://www.pnas.org/doi/abs/10.1073/pnas.2119266119). [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1073/pnas.2119266119&link_type=DOI) 17. Levin, Einav G. et al.ii (Dec. 2021). “Waning Immune Humoral Response to BNT162b2 Covid-19 Vaccine over 6 Months”. In: New England Journal of Medicine 385.24. Publisher: Massachusetts Medical Society _eprint: [https://doi.org/10.1056/NEJMoa2114583](https://doi.org/10.1056/NEJMoa2114583), e84. ISSN: 0028-4793. doi: 10.1056/NEJMoa2114583. URL: [https://doi.org/10.1056/NEJMoa2114583](https://doi.org/10.1056/NEJMoa2114583) (visited on 07/08/2022). [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1056/NEJMoa2114583&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=34614326&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2022%2F08%2F11%2F2022.05.16.22275130.atom) 18. Ma, Yiqun et al.ii (June 2021). “Role of meteorological factors in the transmission of SARS-CoV-2 in the United States”. en. In: Nature Communications 12.1, page 3602. ISSN: 2041-1723. doi: 10.1038/s41467-021-23866-7. URL: [https://www.nature.com/articles/s41467-021-23866-7](https://www.nature.com/articles/s41467-021-23866-7) (visited on 09/25/2021). [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1038/s41467-021-23866-7&link_type=DOI) 19. Nyberg, Tommy et al.ii (2022). “Comparative Analysis of the Risks of Hospitalisation and Death Associated with SARS-CoV-2 Omicron (B.1.1.529) and Delta (B.1.617.2) Variants in England”. en. In: SSRN Electronic Journal. ISSN: 1556-5068. doi: 10.2139/ssrn.4025932. URL: [https://www.ssrn.com/abstract=4025932](https://www.ssrn.com/abstract=4025932) (visited on 03/15/2022). [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.2139/ssrn.4025932&link_type=DOI) 20. Peng, Qiaoli et al.ii (Mar. 2022). “Waning immune responses against SARS-CoV-2 variants of concern among vaccinees in Hong Kong”. English. In: eBioMedicine 77. Publisher: Elsevier. ISSN: 2352-3964. doi: 10.1016/j.ebiom.2022.103904. URL: [https://www.thelancet.com/journals/ebiom/article/PIIS2352-3964(22)00088-3/fulltext](https://www.thelancet.com/journals/ebiom/article/PIIS2352-3964(22)00088-3/fulltext) (visited on 07/08/2022). [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1016/j.ebiom.2022.103904&link_type=DOI) 21. Pérez-Alós, Laura et al.ii (Mar. 2022). “Modeling of waning immunity after SARS-CoV-2 vaccination and influencing factors”. en. In: Nature Communications 13.1. Number: 1 Publisher: Nature Publishing Group, page 1614. ISSN: 2041-1723. doi: 10.1038/s41467-022-29225-4. URL: [https://www.nature.com/articles/s41467-022-29225-4](https://www.nature.com/articles/s41467-022-29225-4) (visited on 07/08/2022). [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1038/s41467-022-29225-4&link_type=DOI) 22. Prunas, Ottavia et al.ii (Mar. 2022). “Vaccination with BNT162b2 reduces transmission of SARS-CoV-2 to household contacts in Israel”. In: Science 375.6585. Publisher: American Association for the Advancement of Science, pages 1151–1154. doi: 10.1126/science.abl4292. URL: [https://www.science.org/doi/10.1126/science.abl4292](https://www.science.org/doi/10.1126/science.abl4292) (visited on 03/24/2022). [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1126/science.abl4292&link_type=DOI) 23. Randolph, Haley E. and Luis B. Barreiro (May 2020). “Herd Immunity: Understanding COVID-19”. en. In: Immunity 52.5, pages 737–741. ISSN: 1074-7613. doi: 10.1016/j.immuni.2020.04.012. URL: [https://www.sciencedirect.com/science/article/pii/S1074761320301709](https://www.sciencedirect.com/science/article/pii/S1074761320301709) (visited on 02/08/2022). [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1016/j.immuni.2020.04.012&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=32433946&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2022%2F08%2F11%2F2022.05.16.22275130.atom) 24. Reyné, Bastien et al.ii (Mar. 2022). “Non-Markovian modelling highlights the importance of age structure on Covid-19 epidemiological dynamics”. en. In: Mathematical Modelling of Natural Phenomena. Publisher: EDP Sciences. ISSN: 0973-5348, 1760-6101. doi: 10.1051/mmnp/2022008. URL: [https://www.mmnp-journal.org/articles/mmnp/abs/forth/mmnp220048/mmnp220048.html](https://www.mmnp-journal.org/articles/mmnp/abs/forth/mmnp220048/mmnp220048.html) (visited on 03/15/2022). [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1051/mmnp/2022008&link_type=DOI) 25. Reynolds, Catherine J. et al.ii (June 2022). “Immune boosting by B.1.1.529 (Omicron) depends on previous SARS-CoV-2 exposure”. In: Science .. Publisher: American Association for the Advancement of Science, eabq1841. doi: 10.1126/science.abq1841. URL: [https://www.science.org/doi/10.1126/science.abq1841](https://www.science.org/doi/10.1126/science.abq1841) (visited on 07/08/2022). [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1126/science.abq1841&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=35699621&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2022%2F08%2F11%2F2022.05.16.22275130.atom) 26. Saad-Roy, Chadi M. et al.ii (Mar. 2021). “Epidemiological and evolutionary considerations of SARS-CoV-2 vaccine dosing regimes”. en. In: Science, eabg8663. ISSN: 0036-8075, 1095-9203. doi: 10.1126/science.abg8663. URL: [https://www.sciencemag.org/lookup/doi/10.1126/science.abg8663](https://www.sciencemag.org/lookup/doi/10.1126/science.abg8663) (visited on 03/23/2021). [Abstract/FREE Full Text](http://medrxiv.org/lookup/ijlink/YTozOntzOjQ6InBhdGgiO3M6MTQ6Ii9sb29rdXAvaWpsaW5rIjtzOjU6InF1ZXJ5IjthOjQ6e3M6ODoibGlua1R5cGUiO3M6NDoiQUJTVCI7czoxMToiam91cm5hbENvZGUiO3M6Mzoic2NpIjtzOjU6InJlc2lkIjtzOjEyOiIzNzIvNjU0MC8zNjMiO3M6NDoiYXRvbSI7czo1MDoiL21lZHJ4aXYvZWFybHkvMjAyMi8wOC8xMS8yMDIyLjA1LjE2LjIyMjc1MTMwLmF0b20iO31zOjg6ImZyYWdtZW50IjtzOjA6IiI7fQ==) 27. Salje, Henrik et al.ii (July 2020). “Estimating the burden of SARS-CoV-2 in France”. In: Science 369.6500. Publisher: American Association for the Advancement of Science, pages 208–211. doi: 10.1126/science.abc3517. URL: [https://www.science.org/doi/10.1126/science.abc3517](https://www.science.org/doi/10.1126/science.abc3517) (visited on 09/25/2021). [Abstract/FREE Full Text](http://medrxiv.org/lookup/ijlink/YTozOntzOjQ6InBhdGgiO3M6MTQ6Ii9sb29rdXAvaWpsaW5rIjtzOjU6InF1ZXJ5IjthOjQ6e3M6ODoibGlua1R5cGUiO3M6NDoiQUJTVCI7czoxMToiam91cm5hbENvZGUiO3M6Mzoic2NpIjtzOjU6InJlc2lkIjtzOjEyOiIzNjkvNjUwMC8yMDgiO3M6NDoiYXRvbSI7czo1MDoiL21lZHJ4aXYvZWFybHkvMjAyMi8wOC8xMS8yMDIyLjA1LjE2LjIyMjc1MTMwLmF0b20iO31zOjg6ImZyYWdtZW50IjtzOjA6IiI7fQ==) 28. Servellita, Venice et al.ii (Apr. 2022). “Neutralizing immunity in vaccine breakthrough infections from the SARS-CoV-2 Omicron and Delta variants”. English. In: Cell 185.9. Publisher: Elsevier, 1539–1548.e5. ISSN: 0092-8674, 1097-4172. doi: 10.1016/j.cell.2022.03.019. URL: [https://www.cell.com/cell/abstract/S0092-8674(22)00329-4](https://www.cell.com/cell/abstract/S0092-8674(22)00329-4) (visited on 07/12/2022). [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1016/j.cell.2022.03.019&link_type=DOI) 29. Sofonea, Mircea T. and Samuel Alizon (Aug. 2021). “Anticipating COVID-19 intensive care unit capacity strain: A look back at epidemiological projections in France”. en. In: Anaesthesia Critical Care & Pain Medicine 40.4, page 100943. ISSN: 2352-5568. doi: 10.1016/j.accpm.2021.100943. URL: [https://www.sciencedirect.com/science/article/pii/S2352556821001478](https://www.sciencedirect.com/science/article/pii/S2352556821001478) (visited on 07/06/2022). [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1016/j.accpm.2021.100943&link_type=DOI) 30. Sofonea, Mircea T., Bastien Reyné, et al.ii (June 2021). “Memory is key in capturing COVID-19 epidemiological dynamics”. en. In: Epidemics 35, page 100459. ISSN: 1755-4365. doi: 10.1016/j.epidem.2021.100459. URL: [https://www.sciencedirect.com/science/article/pii/S1755436521000189](https://www.sciencedirect.com/science/article/pii/S1755436521000189) (visited on 11/17/2021). [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1016/j.epidem.2021.100459&link_type=DOI) 31. Sofonea, Mircea T., Bénédicte Roquebert, et al.ii (2022). “Analyzing and Modeling the Spread of SARS-CoV-2 Omicron Lineages BA.1 and BA.2, France, September 2021–February 2022 - Volume 28, Number 7—July 2022 - Emerging Infectious Diseases journal - CDC”. en-us. In: doi: 10.3201/eid2807.220033. URL: [https://www.nc.cdc.gov/eid/article/28/7/22-0033_article](https://www.nc.cdc.gov/eid/article/28/7/22-0033_article) (visited on 07/08/2022). [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.3201/eid2807.220033&link_type=DOI) 32. Thomine, Olivier et al.ii (Oct. 2021). “Emerging dynamics from high-resolution spatial numerical epidemics”. In: eLife 10.Edited by Talía Malagón and Eduardo Franco. Publisher: eLife Sciences Publications, Ltd, e71417. ISSN: 2050-084X. doi: 10.7554/eLife.71417. URL: [https://doi.org/10.7554/eLife.71417](https://doi.org/10.7554/eLife.71417) (visited on 12/10/2021). [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.7554/eLife.71417&link_type=DOI) 33. UKHSA (2022). “COVID-19 vaccine surveillance report - week 16”. en. In: Report, page 22. URL: [https://www.gov.uk/government/publications/covid-19-vaccine-weekly-surveillance-reports](https://www.gov.uk/government/publications/covid-19-vaccine-weekly-surveillance-reports). 34. UKSHA (2022). “SARS-CoV-2 variants of concern and variants under investigation”. en. In: page 33. URL: [https://www.gov.uk/government/publications/investigation-of-sars-cov-2-variants-technical-briefings](https://www.gov.uk/government/publications/investigation-of-sars-cov-2-variants-technical-briefings). ## References 1. Béraud, Guillaume et al.ii (July 2015). “The French Connection: The First Large Population-Based Contact Survey in France Relevant for the Spread of Infectious Diseases”. en. In: PLOS ONE 10.7. Publisher: Public Library of Science, e0133203. ISSN: 1932-6203. doi: 10.1371/journal.pone.0133203. URL: [https://journals.plos.org/plosone/article?id=10.1371/journal.pone.0133203](https://journals.plos.org/plosone/article?id=10.1371/journal.pone.0133203) (visited on 06/03/2021). [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1371/journal.pone.0133203&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=26176549&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2022%2F08%2F11%2F2022.05.16.22275130.atom) 2. Bosetti, Paolo et al.ii (Jan. 2022). “Impact of booster vaccination on the control of COVID-19 Delta wave in the context of waning immunity: application to France in the winter 2021/22”. en. In: Eurosurveillance 27.1. Publisher: European Centre for Disease Prevention and Control, page 2101125. ISSN: 1560-7917. doi: 10.2807/1560-7917.ES.2022.27.1.2101125. URL: [https://www.eurosurveillance.org/content/10.2807/1560-7917.ES.2022.27.1.2101125](https://www.eurosurveillance.org/content/10.2807/1560-7917.ES.2022.27.1.2101125) (visited on 03/17/2022). [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.2807/1560-7917.ES.2022.27.1.2101125&link_type=DOI) 3. Challen, Robert et al.ii (Mar. 2021). “Risk of mortality in patients infected with SARS-CoV-2 variant of concern 202012/1: matched cohort study”. en. In: BMJ 372. Publisher: British Medical Journal Publishing Group Section: Research, n579. ISSN: 1756-1833. doi: 10.1136/bmj.n579. URL: [https://www.bmj.com/content/372/bmj.n579](https://www.bmj.com/content/372/bmj.n579) (visited on 06/02/2021). [Abstract/FREE Full Text](http://medrxiv.org/lookup/ijlink/YTozOntzOjQ6InBhdGgiO3M6MTQ6Ii9sb29rdXAvaWpsaW5rIjtzOjU6InF1ZXJ5IjthOjQ6e3M6ODoibGlua1R5cGUiO3M6NDoiQUJTVCI7czoxMToiam91cm5hbENvZGUiO3M6MzoiYm1qIjtzOjU6InJlc2lkIjtzOjE3OiIzNzIvbWFyMDlfMTAvbjU3OSI7czo0OiJhdG9tIjtzOjUwOiIvbWVkcnhpdi9lYXJseS8yMDIyLzA4LzExLzIwMjIuMDUuMTYuMjIyNzUxMzAuYXRvbSI7fXM6ODoiZnJhZ21lbnQiO3M6MDoiIjt9) 4. Davies, Nicholas G. et al.ii (May 2021). “Increased mortality in community-tested cases of SARS-CoV-2 lineage B.1.1.7”. en. In: Nature 593.7858. Number: 7858 Publisher: Nature Publishing Group, pages 270–274. ISSN: 1476-4687. doi: 10.1038/s41586-021-03426-1. URL: [https://www.nature.com/articles/s41586-021-03426-1](https://www.nature.com/articles/s41586-021-03426-1) (visited on 06/02/2021). [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1038/s41586-021-03426-1&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=33723411&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2022%2F08%2F11%2F2022.05.16.22275130.atom) 5. Ferretti, Luca et al.ii (May 2020). “Quantifying SARS-CoV-2 transmission suggests epidemic control with digital contact tracing”. en. In: Science 368.6491. Publisher: American Association for the Advancement of Science Section: Research Article. ISSN: 0036-8075, 1095-9203. doi: 10.1126/science.abb6936. URL: [https://science.sciencemag.org/content/368/6491/eabb6936](https://science.sciencemag.org/content/368/6491/eabb6936) (visited on 06/01/2021). [Abstract/FREE Full Text](http://medrxiv.org/lookup/ijlink/YTozOntzOjQ6InBhdGgiO3M6MTQ6Ii9sb29rdXAvaWpsaW5rIjtzOjU6InF1ZXJ5IjthOjQ6e3M6ODoibGlua1R5cGUiO3M6NDoiQUJTVCI7czoxMToiam91cm5hbENvZGUiO3M6Mzoic2NpIjtzOjU6InJlc2lkIjtzOjE3OiIzNjgvNjQ5MS9lYWJiNjkzNiI7czo0OiJhdG9tIjtzOjUwOiIvbWVkcnhpdi9lYXJseS8yMDIyLzA4LzExLzIwMjIuMDUuMTYuMjIyNzUxMzAuYXRvbSI7fXM6ODoiZnJhZ21lbnQiO3M6MDoiIjt9) 6. Fisman, David N. and Ashleigh R. Tuite (Oct. 2021). “Evaluation of the relative virulence of novel SARS-CoV-2 variants: a retrospective cohort study in Ontario, Canada”. en. In: CMAJ 193.42. Publisher: CMAJ Section: Research, E1619–E1625. ISSN: 0820-3946, 1488-2329. doi: 10.1503/cmaj.211248. URL: [https://www.cmaj.ca/content/193/42/E1619](https://www.cmaj.ca/content/193/42/E1619) (visited on 11/03/2021). [Abstract/FREE Full Text](http://medrxiv.org/lookup/ijlink/YTozOntzOjQ6InBhdGgiO3M6MTQ6Ii9sb29rdXAvaWpsaW5rIjtzOjU6InF1ZXJ5IjthOjQ6e3M6ODoibGlua1R5cGUiO3M6NDoiQUJTVCI7czoxMToiam91cm5hbENvZGUiO3M6NDoiY21haiI7czo1OiJyZXNpZCI7czoxMjoiMTkzLzQyL0UxNjE5IjtzOjQ6ImF0b20iO3M6NTA6Ii9tZWRyeGl2L2Vhcmx5LzIwMjIvMDgvMTEvMjAyMi4wNS4xNi4yMjI3NTEzMC5hdG9tIjt9czo4OiJmcmFnbWVudCI7czowOiIiO30=) 7. Hozé, Nathanaël et al.ii (June 2021). “Monitoring the proportion of the population infected by SARS-CoV-2 using age-stratified hospitalisation and serological data: a modelling study”. en. In: The Lancet Public Health 6.6, e408–e415. ISSN: 2468-2667. doi: 10.1016/S2468-2667(21)00064-5. URL: [https://www.sciencedirect.com/science/article/pii/S2468266721000645](https://www.sciencedirect.com/science/article/pii/S2468266721000645) (visited on 06/07/2021). [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1016/S2468-2667(21)00064-5&link_type=DOI) 8. Lefrancq, Noémie et al.ii (2021). “Evolution of outcomes for patients hospitalised during the first 9 months of the SARS-CoV-2 pandemic in France: A retrospective national surveillance data analysis”. en. In: The Lancet Regional Health - Europe 5, page 100087. ISSN: 26667762. doi: 10.1016/j.lanepe.2021.100087. URL: [https://linkinghub.elsevier.com/retrieve/pii/S2666776221000648](https://linkinghub.elsevier.com/retrieve/pii/S2666776221000648) (visited on 04/30/2021). [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1016/j.lanepe.2021.100087&link_type=DOI) 9. Nyberg, Tommy et al.ii (2022). “Comparative Analysis of the Risks of Hospitalisation and Death Associated with SARS-CoV-2 Omicron (B.1.1.529) and Delta (B.1.617.2) Variants in England”. en. In: SSRN Electronic Journal. ISSN: 1556-5068. doi: 10.2139/ssrn.4025932. URL: [https://www.ssrn.com/abstract=4025932](https://www.ssrn.com/abstract=4025932) (visited on 03/15/2022). [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.2139/ssrn.4025932&link_type=DOI) 10. Prunas, Ottavia et al.ii (Mar. 2022). “Vaccination with BNT162b2 reduces transmission of SARS-CoV-2 to household contacts in Israel”. In: Science 375.6585. Publisher: American Association for the Advancement of Science, pages 1151–1154. doi: 10.1126/science.abl4292. URL: [https://www.science.org/doi/10.1126/science.abl4292](https://www.science.org/doi/10.1126/science.abl4292) (visited on 03/24/2022). [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1126/science.abl4292&link_type=DOI) 11. Reyné, Bastien et al.ii (Mar. 2022). “Non-Markovian modelling highlights the importance of age structure on Covid-19 epidemiological dynamics”. en. In: Mathematical Modelling of Natural Phenomena. Publisher: EDP Sciences. ISSN: 0973-5348, 1760-6101. doi: 10.1051/mmnp/2022008. URL: [https://www.mmnp-journal.org/articles/mmnp/abs/forth/mmnp220048/mmnp220048.html](https://www.mmnp-journal.org/articles/mmnp/abs/forth/mmnp220048/mmnp220048.html) (visited on 03/15/2022). [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1051/mmnp/2022008&link_type=DOI) 12. Salje, Henrik et al.ii (July 2020). “Estimating the burden of SARS-CoV-2 in France”. In: Science 369.6500. Publisher: American Association for the Advancement of Science, pages 208–211. doi: 10.1126/science.abc3517. URL: [https://www.science.org/doi/10.1126/science.abc3517](https://www.science.org/doi/10.1126/science.abc3517) (visited on 09/25/2021). [Abstract/FREE Full Text](http://medrxiv.org/lookup/ijlink/YTozOntzOjQ6InBhdGgiO3M6MTQ6Ii9sb29rdXAvaWpsaW5rIjtzOjU6InF1ZXJ5IjthOjQ6e3M6ODoibGlua1R5cGUiO3M6NDoiQUJTVCI7czoxMToiam91cm5hbENvZGUiO3M6Mzoic2NpIjtzOjU6InJlc2lkIjtzOjEyOiIzNjkvNjUwMC8yMDgiO3M6NDoiYXRvbSI7czo1MDoiL21lZHJ4aXYvZWFybHkvMjAyMi8wOC8xMS8yMDIyLjA1LjE2LjIyMjc1MTMwLmF0b20iO31zOjg6ImZyYWdtZW50IjtzOjA6IiI7fQ==) 13. UKHSA (2022). “COVID-19 vaccine surveillance report - week 16”. en. In: Report, page 22. URL: [https://www.gov.uk/government/publications/covid-19-vaccine-weekly-surveillance-reports](https://www.gov.uk/government/publications/covid-19-vaccine-weekly-surveillance-reports). 14. UKSHA (2022). “SARS-CoV-2 variants of concern and variants under investigation”. en. In: page 33. URL: [https://www.gov.uk/government/publications/investigation-of-sars-cov-2-variants-technical-briefings](https://www.gov.uk/government/publications/investigation-of-sars-cov-2-variants-technical-briefings). 15. Verity, Robert et al.ii (June 2020). “Estimates of the severity of coronavirus disease 2019: a model-based analysis”. English. In: The Lancet Infectious Diseases 20.6. Publisher: Elsevier, pages 669–677. ISSN: 1473-3099, 1474-4457. doi: 10.1016/S1473-3099(20)30243-7. URL: [https://www.thelancet.com/journals/laninf/article/PIIS1473-3099(20)30243-7/abstract](https://www.thelancet.com/journals/laninf/article/PIIS1473-3099(20)30243-7/abstract) (visited on 06/02/2021). [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1016/S1473-3099(20)30243-7&link_type=DOI) [1]: /embed/inline-graphic-1.gif [2]: /embed/inline-graphic-2.gif [3]: /embed/inline-graphic-3.gif [4]: /embed/inline-graphic-4.gif [5]: /embed/inline-graphic-5.gif [6]: F1/embed/inline-graphic-6.gif [7]: F1/embed/inline-graphic-7.gif [8]: /embed/graphic-2.gif [9]: /embed/graphic-3.gif [10]: /embed/graphic-15.gif [11]: /embed/graphic-16.gif [12]: /embed/graphic-17.gif [13]: /embed/graphic-18.gif [14]: /embed/graphic-19.gif [15]: /embed/graphic-20.gif [16]: /embed/graphic-21.gif [17]: /embed/graphic-22.gif [18]: /embed/graphic-23.gif [19]: /embed/graphic-24.gif [20]: /embed/graphic-25.gif [21]: /embed/graphic-26.gif [22]: /embed/graphic-27.gif [23]: /embed/graphic-28.gif [24]: /embed/graphic-29.gif [25]: /embed/graphic-30.gif