Mathematical models tend to oversimplify the biological details of vaccine or infection-derived immunity effectiveness. Yet, epidemiological outcomes may diverge when assuming polarised immunity-individuals are either fully susceptible or completely immune-compared to leaky immunity-where all individuals are partially protected. We explore the differences between the two by taking advantage of a non-Markovian framework, which allows us to explicitly record the 'age' of the immunity and vary its effectiveness accordingly. A basic scenario reveals that leaky immunity leads to a shorter time between reinfections. A more data-driven scenario based on SARS-CoV-2 data finds that leaky immunity yields substantially more reinfections than polarised immunity and a higher number of infected individuals, yet with a lower probability of hospitalisation. Our findings emphasize the critical role of immune memory modelling assumptions, especially for long-term epidemiological dynamics and public health policies.

Projects such as the European Covid-19 Forecast Hub publish forecasts on the national level for new deaths, new cases, and hospital admissions, but not direct measurements of hospital strain like critical care bed occupancy at the sub-national level, which is of particular interest to health professionals for planning purposes. We present a sub-national French framework for forecasting hospital strain based on a non-Markovian compartmental model, its associated online visualisation tool and a retrospective evaluation of the real-time forecasts it provided from January to December 2021 by comparing to three baselines derived from standard statistical forecasting methods (a naive model, auto-regression, and an ensemble of exponential smoothing and ARIMA). In terms of median absolute error for forecasting critical care unit occupancy at the two-week horizon, our model only outperformed the naive baseline for 4 out of 14 geographical units and underperformed compared to the ensemble baseline for 5 of them at the 90% confidence level (n = 38). However, for the same level at the 4 week horizon, our model was never statistically outperformed for any unit despite outperforming the baselines 10 times spanning 7 out of 14 geographical units. This implies modest forecasting utility for longer horizons which may justify the application of non-Markovian compartmental models in the context of hospital-strain surveillance for future pandemics.

Human Papillomaviruses (HPV) are one of the most prevalent sexually transmitted infections (STI) and the most oncogenic viruses known to humans. The vast majority of HPV infections clear in less than 3 years, but the underlying mechanisms, especially the involvement of the immune response, are still poorly known. Building on earlier work stressing the potential importance of stochasticity in HPV infection clearance, we develop a stochastic mathematical model of HPV dynamics with an explicit description of the intracellular level. We show that the random partitioning of virus episomes upon stem cell division and the occurrence of symmetric divisions greatly affect viral persistence. These results call for more detailed within-host studies to better understand the relative importance of stochasticity and immunity in HPV infection clearance.

In an epidemic, individuals can widely differ in the way they spread the infection depending on their age or on the number of days they have been infected for. In the absence of pharmaceutical interventions such as a vaccine or treatment, non-pharmaceutical interventions (e.g. physical or social distancing) are essential to mitigate the pandemic. We develop an original approach to identify the optimal age-stratified control strategy to implement as a function of the time since the onset of the epidemic. This is based on a model with a double continuous structure in terms of host age and time since infection. By applying optimal control theory to this model, we identify a solution that minimizes deaths and costs associated with the implementation of the control strategy itself. We also implement this strategy for three countries with contrasted age distributions (Burkina-Faso, France, and Vietnam). Overall, the optimal strategy varies throughout the epidemic, with a more intense control early on, and depending on host age, with a stronger control for the older population, except in the scenario where the cost associated with the control is low. In the latter scenario, we find strong differences across countries because the control extends to the younger population for France and Vietnam 2 to 3 months after the onset of the epidemic, but not for Burkina Faso. Finally, we show that the optimal control strategy strongly outperforms a constant uniform control exerted over the whole population or over its younger fraction. This improved understanding of the effect of age-based control interventions opens new perspectives for the field, especially for age-based contact tracing.

Interactions within a population shape the spread of infectious diseases but contact patterns between individuals are difficult to access. We hypothesised that key properties of these patterns can be inferred from multiple infection data in longitudinal follow-ups. We developed a simulator for epidemics with multiple infections on networks and analysed the resulting individual infection time series by introducing similarity metrics between hosts based on their multiple infection histories. We find that, depending on infection multiplicity and network sampling, multiple infection summary statistics can recover network properties such as degree distribution. Furthermore, we show that by mining simulation outputs for multiple infection patterns, one can detect immunological interference between pathogens (i.e. the fact that past infections in a host condition future probability of infection). The combination of individual-based simulations and analysis of multiple infection histories opens promising perspectives to infer and validate transmission networks and immunological interference for infectious diseases from longitudinal cohort data.

Infections of stratified epithelia contribute to a large group of common diseases, such as dermatological conditions and sexually transmitted diseases. To investigate how epithelial structure affects infection dynamics, we develop a general ecology-inspired model for stratified epithelia. Our model allows us to simulate infections, explore new hypotheses and estimate parameters that are difficult to measure with tissue cell cultures. We focus on two contrasting pathogens: Chlamydia trachomatis and Human papillomaviruses (HPV). Using cervicovaginal parameter estimates, we find that key infection symptoms can be explained by differential interactions with the layers, while clearance and pathogen burden appear to be bottom-up processes. Cell protective responses to infections (e.g. mucus trapping) generally lowered pathogen load but there were specific effects based on infection strategies. Our modeling approach opens new perspectives for 3D tissue culture experimental systems of infections and, more generally, for developing and testing hypotheses related to infections of stratified epithelia.

Inferring epidemiological parameters such as the R0 from time-scaled phylogenies is a timely challenge. Most current approaches rely on likelihood functions, which raise specific issues that range from computing these functions to finding their maxima numerically. Here, we present a new regression-based Approximate Bayesian Computation (ABC) approach, which we base on a large variety of summary statistics intended to capture the information contained in the phylogeny and its corresponding lineage-through-time plot. The regression step involves the Least Absolute Shrinkage and Selection Operator (LASSO) method, which is a robust machine learning technique. It allows us to readily deal with the large number of summary statistics, while avoiding resorting to Markov Chain Monte Carlo (MCMC) techniques. To compare our approach to existing ones, we simulated target trees under a variety of epidemiological models and settings, and inferred parameters of interest using the same priors. We found that, for large phylogenies, the accuracy of our regression-ABC is comparable to that of likelihood-based approaches involving birth-death processes implemented in BEAST2. Our approach even outperformed these when inferring the host population size with a Susceptible-Infected-Removed epidemiological model. It also clearly outperformed a recent kernel-ABC approach when assuming a Susceptible-Infected epidemiological model with two host types. Lastly, by re-analyzing data from the early stages of the recent Ebola epidemic in Sierra Leone, we showed that regression-ABC provides more realistic estimates for the duration parameters (latency and infectiousness) than the likelihood-based method. Overall, ABC based on a large variety of summary statistics and a regression method able to perform variable selection and avoid overfitting is a promising approach to analyze large phylogenies.

Predicting the emergence of new pathogenic strains is a key goal of evolutionary epidemiology. However, the majority of existing studies have focussed on emergence at the population level, and not within a host. In particular, the coexistence of pre-existing and mutated strains triggers a heightened immune response due to the larger total pathogen population; this feedback can smother mutated strains before they reach an ample size and establish. Here, we extend previous work for measuring emergence probabilities in non-equilibrium populations, to within-host models of acute infections. We create a mathematical model to investigate the emergence probability of a fitter strain if it mutates from a self-limiting strain that is guaranteed to go extinct in the long-term. We show that ongoing immune cell proliferation during the initial stages of infection causes a drastic reduction in the probability of emergence of mutated strains; we further outline how this effect can be accurately measured. Further analysis of the model shows that, in the short-term, mutant strains that enlarge their replication rate due to evolving an increased growth rate are more favoured than strains that suffer a lower immune-mediated death rate ('immune tolerance'), as the latter does not completely evade ongoing immune proliferation due to inter-parasitic competition. We end by discussing the model in relation to within-host evolution of human pathogens (including HIV, hepatitis C virus, and cancer), and how ongoing immune growth can affect their evolutionary dynamics.

The contact structure between hosts shapes disease spread. Most network-based models used in epidemiology tend to ignore heterogeneity in the weighting of contacts between two individuals. However, this assumption is known to be at odds with the data for many networks (e.g. sexual contact networks) and to have a critical influence on epidemics' behavior. One of the reasons why models usually ignore heterogeneity in transmission is that we currently lack tools to analyze weighted networks, such that most studies rely on numerical simulations. Here, we present a novel framework to estimate key epidemiological variables, such as the rate of early epidemic expansion (r0) and the basic reproductive ratio (R0), from joint probability distributions of number of partners (contacts) and number of interaction events through which contacts are weighted. These distributions are much easier to infer than the exact shape of the network, which makes the approach widely applicable. The framework also allows for a derivation of the full time course of epidemic prevalence and contact behaviour, which we validate with numerical simulations on networks. Overall, incorporating more realistic contact networks into epidemiological models can improve our understanding of the emergence and spread of infectious diseases.

Many pathogens associated with chronic infections evolve so rapidly that strains found late in an infection have little in common with the initial strain. This raises questions at different levels of analysis because rapid within-host evolution affects the course of an infection, but it can also affect the possibility for natural selection to act at the between-host level. We present a nested approach that incorporates within-host evolutionary dynamics of a rapidly mutating virus (hepatitis C virus) targeted by a cellular cross-reactive immune response, into an epidemiological perspective. The viral trait we follow is the replication rate of the strain initiating the infection. We find that, even for rapidly evolving viruses, the replication rate of the initial strain has a strong effect on the fitness of an infection. Moreover, infections caused by slowly replicating viruses have the highest infection fitness (i.e., lead to more secondary infections), but strains with higher replication rates tend to dominate within a host in the long-term. We also study the effect of cross-reactive immunity and viral mutation rate on infection life history traits. For instance, because of the stochastic nature of our approach, we can identify factors affecting the outcome of the infection (acute or chronic infections). Finally, we show that anti-viral treatments modify the value of the optimal initial replication rate and that the timing of the treatment administration can have public health consequences due to within-host evolution. Our results support the idea that natural selection can act on the replication rate of rapidly evolving viruses at the between-host level. It also provides a mechanistic description of within-host constraints, such as cross-reactive immunity, and shows how these constraints affect the infection fitness. This model raises questions that can be tested experimentally and underlines the necessity to consider the evolution of quantitative traits to understand the outcome and the fitness of an infection.