Using observed incidence to calibrate the transmission level of a mathematical model for Plasmodium vivax dynamics including case management and importation.

In this work, we present a simple and flexible model for Plasmodium vivax dynamics which can be easily combined with routinely collected data on local and imported case counts to quantify transmission intensity and simulate control strategies. This model extends the model from White et al. (2016) by including case management interventions targeting liver-stage or blood-stage parasites, as well as imported infections. The endemic steady state of the model is used to derive a relationship between the observed incidence and the transmission rate in order to calculate reproduction numbers and simulate intervention scenarios. To illustrate its potential applications, the model is used to calculate local reproduction numbers in Panama and identify areas of sustained malaria transmission that should be targeted by control interventions.

Improving access to care and community health in Haiti with optimized community health worker placement.

The national deployment of polyvalent community health workers (CHWs) is a constitutive part of the strategy initiated by the Ministry of Health to accelerate efforts towards universal health coverage in Haiti. Its implementation requires the planning of future recruitment and deployment activities for which mathematical modelling tools can provide useful support by exploring optimised placement scenarios based on access to care and population distribution. We combined existing gridded estimates of population and travel times with optimisation methods to derive theoretical CHW geographical placement scenarios including constraints on walking time and the number of people served per CHW. Four national-scale scenarios that align with total numbers of existing CHWs and that ensure that the walking time for each CHW does not exceed a predefined threshold are compared. The first scenario accounts for population distribution in rural and urban areas only, while the other three also incorporate in different ways the proximity of existing health centres. Comparing these scenarios to the current distribution, insufficient number of CHWs is systematically identified in several departments and gaps in access to health care are identified within all departments. These results highlight current suboptimal distribution of CHWs and emphasize the need to consider an optimal (re-)allocation.

Dynamics of the COVID-19 epidemic in Ireland under mitigation.

BACKGROUND: In Ireland and across the European Union the COVID-19 epidemic waves, driven mainly by the emergence of new variants of the SARS-CoV-2 have continued their course, despite various interventions from governments. Public health interventions continue in their attempts to control the spread as they wait for the planned significant effect of vaccination. METHODS: To tackle this challenge and the observed non-stationary aspect of the epidemic we used a modified SEIR stochastic model with time-varying parameters, following Brownian process. This enabled us to reconstruct the temporal evolution of the transmission rate of COVID-19 with the non-specific hypothesis that it follows a basic stochastic process constrained by the available data. This model is coupled with Bayesian inference (particle Markov Chain Monte Carlo method) for parameter estimation and utilized mainly well-documented Irish hospital data. RESULTS: In Ireland, mitigation measures provided a 78-86% reduction in transmission during the first wave between March and May 2020. For the second wave in October 2020, our reduction estimation was around 20% while it was 70% for the third wave in January 2021. This third wave was partly due to the UK variant appearing in Ireland. In June 2020 we estimated that sero-prevalence was 2.0% (95% CI: 1.2-3.5%) in complete accordance with a sero-prevalence survey. By the end of April 2021, the sero-prevalence was greater than 17% due in part to the vaccination campaign. Finally we demonstrate that the available observed confirmed cases are not reliable for analysis owing to the fact that their reporting rate has as expected greatly evolved. CONCLUSION: We provide the first estimations of the dynamics of the COVID-19 epidemic in Ireland and its key parameters. We also quantify the effects of mitigation measures on the virus transmission during and after mitigation for the three waves. Our results demonstrate that Ireland has significantly reduced transmission by employing mitigation measures, physical distancing and lockdown. This has to date avoided the saturation of healthcare infrastructures, flattened the epidemic curve and likely reduced mortality. However, as we await for a full roll out of a vaccination programme and as new variants potentially more transmissible and/or more infectious could continue to emerge and mitigation measures change silent transmission, challenges remain.

A mechanistic and data-driven reconstruction of the time-varying reproduction number: Application to the COVID-19 epidemic.

The effective reproduction number Reff is a critical epidemiological parameter that characterizes the transmissibility of a pathogen. However, this parameter is difficult to estimate in the presence of silent transmission and/or significant temporal variation in case reporting. This variation can occur due to the lack of timely or appropriate testing, public health interventions and/or changes in human behavior during an epidemic. This is exactly the situation we are confronted with during this COVID-19 pandemic. In this work, we propose to estimate Reff for the SARS-CoV-2 (the etiological agent of the COVID-19), based on a model of its propagation considering a time-varying transmission rate. This rate is modeled by a Brownian diffusion process embedded in a stochastic model. The model is then fitted by Bayesian inference (particle Markov Chain Monte Carlo method) using multiple well-documented hospital datasets from several regions in France and in Ireland. This mechanistic modeling framework enables us to reconstruct the temporal evolution of the transmission rate of the COVID-19 based only on the available data. Except for the specific model structure, it is non-specifically assumed that the transmission rate follows a basic stochastic process constrained by the observations. This approach allows us to follow both the course of the COVID-19 epidemic and the temporal evolution of its Reff(t). Besides, it allows to assess and to interpret the evolution of transmission with respect to the mitigation strategies implemented to control the epidemic waves in France and in Ireland. We can thus estimate a reduction of more than 80% for the first wave in all the studied regions but a smaller reduction for the second wave when the epidemic was less active, around 45% in France but just 20% in Ireland. For the third wave in Ireland the reduction was again significant (>70%).

Parallel trends in the transmission of SARS-CoV-2 and retail/recreation and public transport mobility during non-lockdown periods.

Recent literature strongly supports the hypothesis that mobility restriction and social distancing play a crucial role in limiting the transmission of severe acute respiratory syndrome coronavirus-2 (SARS-CoV-2). During the first wave of the coronavirus disease 2019 (COVID-19) pandemic, it was shown that mobility restriction reduced transmission significantly. This study found that, in the period between the first two waves of the COVID-19 pandemic, there was high positive correlation between trends in the transmission of SARS-CoV-2 and mobility. These two trends oscillated simultaneously, and increased mobility following the relaxation of lockdown rules was significantly associated with increased transmission. From a public health perspective, these results highlight the importance of tracking changes in mobility when relaxing mitigation measures in order to anticipate future changes in the spread of SARS-CoV-2.

Predicting the impact of outdoor vector control interventions on malaria transmission intensity from semi-field studies.

BACKGROUND: Semi-field experiments with human landing catch (HLC) measure as the outcome are an important step in the development of novel vector control interventions against outdoor transmission of malaria since they provide good estimates of personal protection. However, it is often infeasible to determine whether the reduction in HLC counts is due to mosquito mortality or repellency, especially considering that spatial repellents based on volatile pyrethroids might induce both. Due to the vastly different impact of repellency and mortality on transmission, the community-level impact of spatial repellents can not be estimated from such semi-field experiments. METHODS: We present a new stochastic model that is able to estimate for any product inhibiting outdoor biting, its repelling effect versus its killing and disarming (preventing host-seeking until the next night) effects, based only on time-stratified HLC data from controlled semi-field experiments. For parameter inference, a Bayesian hierarchical model is used to account for nightly variation of semi-field experimental conditions. We estimate the impact of the products on the vectorial capacity of the given Anopheles species using an existing mathematical model. With this methodology, we analysed data from recent semi-field studies in Kenya and Tanzania on the impact of transfluthrin-treated eave ribbons, the odour-baited Suna trap and their combination (push-pull system) on HLC of Anopheles arabiensis in the peridomestic area. RESULTS: Complementing previous analyses of personal protection, we found that the transfluthrin-treated eave ribbons act mainly by killing or disarming mosquitoes. Depending on the actual ratio of disarming versus killing, the vectorial capacity of An. arabiensis is reduced by 41 to 96% at 70% coverage with the transfluthrin-treated eave ribbons and by 38 to 82% at the same coverage with the push-pull system, under the assumption of a similar impact on biting indoors compared to outdoors. CONCLUSIONS: The results of this analysis of semi-field data suggest that transfluthrin-treated eave ribbons are a promising tool against malaria transmission by An. arabiensis in the peridomestic area, since they provide both personal and community protection. Our modelling framework can estimate the community-level impact of any tool intervening during the mosquito host-seeking state using data from only semi-field experiments with time-stratified HLC.

Aedes mosquitoes acquire and transmit Zika virus by breeding in contaminated aquatic environments.

Zika virus (ZIKV) is a mosquito-borne flavivirus that predominantly circulates between humans and Aedes mosquitoes. Clinical studies have shown that Zika viruria in patients persists for an extended period, and results in infectious virions being excreted. Here, we demonstrate that Aedes mosquitoes are permissive to ZIKV infection when breeding in urine or sewage containing low concentrations of ZIKV. Mosquito larvae and pupae, including from field Aedes aegypti can acquire ZIKV from contaminated aquatic systems, resulting in ZIKV infection of adult females. Adult mosquitoes can transmit infectious virions to susceptible type I/II interferon receptor-deficient (ifnagr-/-) C57BL/6 (AG6) mice. Furthermore, ZIKV viruria from infected AG6 mice can causes mosquito infection during the aquatic life stages. Our studies suggest that infectious urine could be a natural ZIKV source, which is potentially transmissible to mosquitoes when breeding in an aquatic environment.

Comparison of stochastic and deterministic frameworks in dengue modelling.

We perform estimations of compartment models for dengue transmission in rural Cambodia with increasing complexity regarding both model structure and the account for stochasticity. On the one hand, we successively account for three embedded sources of stochasticity: observation noise, demographic variability and environmental hazard. On the other hand, complexity in the model structure is increased by introducing vector-borne transmission, explicit asymptomatic infections and interacting virus serotypes. Using two sources of case data from dengue epidemics in Kampong Cham (Cambodia), models are estimated in the bayesian framework, with Markov Chain Monte Carlo and Particle Markov Chain Monte Carlo. We highlight the advantages and drawbacks of the different formulations in a practical setting. Although in this case the deterministic models provide a good approximation of the mean trajectory for a low computational cost, the stochastic frameworks better reflect and account for parameter and simulation uncertainty.

Dengue modeling in rural Cambodia: Statistical performance versus epidemiological relevance.

Dengue dynamics are shaped by the complex interplay between several factors, including vector seasonality, interaction between four virus serotypes, and inapparent infections. However, paucity or quality of data do not allow for all of these to be taken into account in mathematical models. In order to explore separately the importance of these factors in models, we combined surveillance data with a local-scale cluster study in the rural province of Kampong Cham (Cambodia), in which serotypes and asymptomatic infections were documented. We formulate several mechanistic models, each one relying on a different set of hypotheses, such as explicit vector dynamics, transmission via asymptomatic infections and coexistence of several virus serotypes. Models are confronted with the observed time series using Bayesian inference, through Markov chain Monte Carlo. Model selection is then performed using statistical information criteria, and the coherence of epidemiological characteristics (reproduction numbers, incidence proportion, dynamics of the susceptible classes) is assessed in each model. Our analyses on transmission dynamics in a rural endemic setting highlight that two-strain models with interacting effects better reproduce the long term data, but they are difficult to parameterize when relying on incidence cases only. On the other hand, considering the available data, incorporating vector and asymptomatic components seems of limited added-value when seasonality and underreporting are already accounted for.

Accounting for non-stationarity in epidemiology by embedding time-varying parameters in stochastic models.

The spread of disease through human populations is complex. The characteristics of disease propagation evolve with time, as a result of a multitude of environmental and anthropic factors, this non-stationarity is a key factor in this huge complexity. In the absence of appropriate external data sources, to correctly describe the disease propagation, we explore a flexible approach, based on stochastic models for the disease dynamics, and on diffusion processes for the parameter dynamics. Using such a diffusion process has the advantage of not requiring a specific mathematical function for the parameter dynamics. Coupled with particle MCMC, this approach allows us to reconstruct the time evolution of some key parameters (average transmission rate for instance). Thus, by capturing the time-varying nature of the different mechanisms involved in disease propagation, the epidemic can be described. Firstly we demonstrate the efficiency of this methodology on a toy model, where the parameters and the observation process are known. Applied then to real datasets, our methodology is able, based solely on simple stochastic models, to reconstruct complex epidemics, such as flu or dengue, over long time periods. Hence we demonstrate that time-varying parameters can improve the accuracy of model performances, and we suggest that our methodology can be used as a first step towards a better understanding of a complex epidemic, in situation where data is limited and/or uncertain.

Structure in the variability of the basic reproductive number (R0) for Zika epidemics in the Pacific islands.

Before the outbreak that reached the Americas in 2015, Zika virus (ZIKV) circulated in Asia and the Pacific: these past epidemics can be highly informative on the key parameters driving virus transmission, such as the basic reproduction number (R0). We compare two compartmental models with different mosquito representations, using surveillance and seroprevalence data for several ZIKV outbreaks in Pacific islands (Yap, Micronesia 2007, Tahiti and Moorea, French Polynesia 2013-2014, New Caledonia 2014). Models are estimated in a stochastic framework with recent Bayesian techniques. R0 for the Pacific ZIKV epidemics is estimated between 1.5 and 4.1, the smallest islands displaying higher and more variable values. This relatively low range of R0 suggests that intervention strategies developed for other flaviviruses should enable as, if not more effective control of ZIKV. Our study also highlights the importance of seroprevalence data for precise quantitative analysis of pathogen propagation, to design prevention and control strategies.