On the predictive ability of mechanistic models for the Haitian cholera epidemic

Informing policy via dynamic models: Cholera in Haiti

Public health decisions must be made about when and how to implement interventions to control an infectious disease epidemic. These decisions should be informed by data on the epidemic as well as current understanding about the transmission dynamics. Such decisions can be posed as statistical questions about scientifically motivated dynamic models. Thus, we encounter the methodological task of building credible, data-informed decisions based on stochastic, partially observed, nonlinear dynamic models. This necessitates addressing the tradeoff between biological fidelity and model simplicity, and the reality of misspecification for models at all levels of complexity. We assess current methodological approaches to these issues via a case study of the 2010-2019 cholera epidemic in Haiti. We consider three dynamic models developed by expert teams to advise on vaccination policies. We evaluate previous methods used for fitting these models, and we demonstrate modified data analysis strategies leading to improved statistical fit. Specifically, we present approaches for diagnosing model misspecification and the consequent development of improved models. Additionally, we demonstrate the utility of recent advances in likelihood maximization for high-dimensional nonlinear dynamic models, enabling likelihood-based inference for spatiotemporal incidence data using this class of models. Our workflow is reproducible and extendable, facilitating future investigations of this disease system.

Epidemicity of cholera spread and the fate of infection control measures

The fate of ongoing infectious disease outbreaks is predicted through reproduction numbers, defining the long-term establishment of the infection, and epidemicity indices, tackling the reactivity of the infectious pool to new contagions. Prognostic metrics of unfolding outbreaks are of particular importance when designing adaptive emergency interventions facing real-time assimilation of epidemiological evidence. Our aim here is twofold. First, we propose a novel form of the epidemicity index for the characterization of cholera epidemics in spatial models of disease spread. Second, we examine in hindsight the survey of infections, treatments and containment measures carried out for the now extinct 2010–2019 Haiti cholera outbreak, to suggest that magnitude and timing of non-pharmaceutical and vaccination interventions imply epidemiological responses recapped by the evolution of epidemicity indices. Achieving negative epidemicity greatly accelerates fading of infections and thus proves a worthwhile target of containment measures. We also show that, in our model, effective reproduction numbers and epidemicity indices are explicitly related. Therefore, providing an upper bound to the effective reproduction number (significantly lower than the unit threshold) warrants negative epidemicity and, in turn, a rapidly fading outbreak preventing coalescence of sparse local sub-threshold flare-ups.

Spread and dynamics of the COVID-19 epidemic in Italy: Effects of emergency containment measures

The spread of coronavirus disease 2019 (COVID-19) in Italy prompted drastic measures for transmission containment. We examine the effects of these interventions, based on modeling of the unfolding epidemic. We test modeling options of the spatially explicit type, suggested by the wave of infections spreading from the initial foci to the rest of Italy. We estimate parameters of a metacommunity Susceptible-Exposed-Infected-Recovered (SEIR)-like transmission model that includes a network of 107 provinces connected by mobility at high resolution, and the critical contribution of presymptomatic and asymptomatic transmission. We estimate a generalized reproduction number ([Formula: see text] = 3.60 [3.49 to 3.84]), the spectral radius of a suitable next-generation matrix that measures the potential spread in the absence of containment interventions. The model includes the implementation of progressive restrictions after the first case confirmed in Italy (February 21, 2020) and runs until March 25, 2020. We account for uncertainty in epidemiological reporting, and time dependence of human mobility matrices and awareness-dependent exposure probabilities. We draw scenarios of different containment measures and their impact. Results suggest that the sequence of restrictions posed to mobility and human-to-human interactions have reduced transmission by 45% (42 to 49%). Averted hospitalizations are measured by running scenarios obtained by selectively relaxing the imposed restrictions and total about 200,000 individuals (as of March 25, 2020). Although a number of assumptions need to be reexamined, like age structure in social mixing patterns and in the distribution of mobility, hospitalization, and fatality, we conclude that verifiable evidence exists to support the planning of emergency measures.

Quantitative Microbial Risk Assessment and Infectious Disease Transmission Modeling of Waterborne Enteric Pathogens

PURPOSE OF REVIEW

Waterborne enteric pathogens remain a global health threat. Increasingly, quantitative microbial risk assessment (QMRA) and infectious disease transmission modeling (IDTM) are used to assess waterborne pathogen risks and evaluate mitigation. These modeling efforts, however, have largely been conducted independently for different purposes and in different settings. In this review, we examine the settings where each modeling strategy is employed.

RECENT FINDINGS

QMRA research has focused on food contamination and recreational water in high-income countries (HICs) and drinking water and wastewater in low- and middle-income countries (LMICs). IDTM research has focused on large outbreaks (predominately LMICs) and vaccine-preventable diseases (LMICs and HICs). Human ecology determines the niches that pathogens exploit, leading researchers to focus on different risk assessment research strategies in different settings. To enhance risk modeling, QMRA and IDTM approaches should be integrated to include dynamics of pathogens in the environment and pathogen transmission through populations.

Conditions for transient epidemics of waterborne disease in spatially explicit systems

Waterborne diseases are a diverse family of infections transmitted through ingestion of-or contact with-water infested with pathogens. Outbreaks of waterborne infections often show well-defined spatial signatures that are typically linked to local eco-epidemiological conditions, water-mediated pathogen transport and human mobility. In this work, we apply a spatially explicit network model describing the transmission cycle of waterborne pathogens to determine invasion conditions in metacommunities endowed with a realistic spatial structure. Specifically, we aim to define conditions under which pathogens can temporarily colonize a set of human communities, thus triggering a transient epidemic outbreak. To that end, we apply generalized reactivity analysis, a recently developed methodological framework for the study of transient dynamics in ecological systems subject to external perturbations. The study of pathogen invasion is complemented by the detection of the spatial signatures associated with the perturbations to a disease-free system that are expected to be amplified the most over different time scales. Understanding the drivers of waterborne disease dynamics over time scales that are relevant to epidemic and/or endemic transmission is a crucial, cross-disciplinary challenge, as large portions of the developing world still struggle to cope with the burden of these infections.

Rainfall as a driver of epidemic cholera: Comparative model assessments of the effect of intra-seasonal precipitation events

The correlation between cholera epidemics and climatic drivers, in particular seasonal tropical rainfall, has been studied in a variety of contexts owing to its documented relevance. Several mechanistic models of cholera transmission have included rainfall as a driver by focusing on two possible transmission pathways: either by increasing exposure to contaminated water (e.g. due to worsening sanitary conditions during water excess), or water contamination by freshly excreted bacteria (e.g. due to washout of open-air defecation sites or overflows). Our study assesses the explanatory power of these different modeling structures by formal model comparison using deterministic and stochastic models of the type susceptible-infected-recovered-bacteria (SIRB). The incorporation of rainfall effects is generalized using a nonlinear function that can increase or decrease the relative importance of the large precipitation events. Our modelling framework is tested against the daily epidemiological data collected during the 2015 cholera outbreak within the urban context of Juba, South Sudan. This epidemic is characterized by a particular intra-seasonal double peak on the incidence in apparent relation with particularly strong rainfall events. Our results show that rainfall-based models in both their deterministic and stochastic formulations outperform models that do not account for rainfall. In fact, classical SIRB models are not able to reproduce the second epidemiological peak, thus suggesting that it was rainfall-driven. Moreover we found stronger support across model types for rainfall acting on increased exposure rather than on exacerbated water contamination. Although these results are context-specific, they stress the importance of a systematic and comprehensive appraisal of transmission pathways and their environmental forcings when embarking in the modelling of epidemic cholera.

Near real-time forecasting for cholera decision making in Haiti after Hurricane Matthew

Computational models of cholera transmission can provide objective insights into the course of an ongoing epidemic and aid decision making on allocation of health care resources. However, models are typically designed, calibrated and interpreted post-hoc. Here, we report the efforts of a team from academia, field research and humanitarian organizations to model in near real-time the Haitian cholera outbreak after Hurricane Matthew in October 2016, to assess risk and to quantitatively estimate the efficacy of a then ongoing vaccination campaign. A rainfall-driven, spatially-explicit meta-community model of cholera transmission was coupled to a data assimilation scheme for computing short-term projections of the epidemic in near real-time. The model was used to forecast cholera incidence for the months after the passage of the hurricane (October-December 2016) and to predict the impact of a planned oral cholera vaccination campaign. Our first projection, from October 29 to December 31, predicted the highest incidence in the departments of Grande Anse and Sud, accounting for about 45% of the total cases in Haiti. The projection included a second peak in cholera incidence in early December largely driven by heavy rainfall forecasts, confirming the urgency for rapid intervention. A second projection (from November 12 to December 31) used updated rainfall forecasts to estimate that 835 cases would be averted by vaccinations in Grande Anse (90% Prediction Interval [PI] 476-1284) and 995 in Sud (90% PI 508-2043). The experience gained by this modeling effort shows that state-of-the-art computational modeling and data-assimilation methods can produce informative near real-time projections of cholera incidence. Collaboration among modelers and field epidemiologists is indispensable to gain fast access to field data and to translate model results into operational recommendations for emergency management during an outbreak. Future efforts should thus draw together multi-disciplinary teams to ensure model outputs are appropriately based, interpreted and communicated.

Epidemicity thresholds for water-borne and water-related diseases

Determining the conditions that favor pathogen establishment in a host community is key to disease control and eradication. However, focusing on long-term dynamics alone may lead to an underestimation of the threats imposed by outbreaks triggered by short-term transient phenomena. Achieving an effective epidemiological response thus requires to look at different timescales, each of which may be endowed with specific management objectives. In this work we aim to determine epidemicity thresholds for some prototypical examples of water-borne and water-related diseases, a diverse family of infections transmitted either directly through water infested with pathogens or by vectors whose lifecycles are closely associated with water. From a technical perspective, while conditions for endemicity are determined via stability analysis, epidemicity thresholds are defined through generalized reactivity analysis, a recently proposed method that allows the study of the short-term instability properties of ecological systems. Understanding the drivers of water-borne and water-related disease dynamics over timescales that may be relevant to epidemic and/or endemic transmission is a challenge of the utmost importance, as large portions of the developing world are still struggling with the burden imposed by these infections.

River networks as ecological corridors: A coherent ecohydrological perspective

This paper draws together several lines of argument to suggest that an ecohydrological framework, i.e. laboratory, field and theoretical approaches focused on hydrologic controls on biota, has contributed substantially to our understanding of the function of river networks as ecological corridors. Such function proves relevant to: the spatial ecology of species; population dynamics and biological invasions; the spread of waterborne disease. As examples, we describe metacommunity predictions of fish diversity patterns in the Mississippi-Missouri basin, geomorphic controls imposed by the fluvial landscape on elevational gradients of species' richness, the zebra mussel invasion of the same Mississippi-Missouri river system, and the spread of proliferative kidney disease in salmonid fish. We conclude that spatial descriptions of ecological processes in the fluvial landscape, constrained by their specific hydrologic and ecological dynamics and by the ecosystem matrix for interactions, i.e. the directional dispersal embedded in fluvial and host/pathogen mobility networks, have already produced a remarkably broad range of significant results. Notable scientific and practical perspectives are thus open, in the authors' view, to future developments in ecohydrologic research.

Transmission dynamics of cholera in Yemen, 2017: a real time forecasting

BACKGROUND

A large epidemic of cholera, caused by Vibrio cholerae, serotype Ogawa, has been ongoing in Yemen, 2017. To improve the situation awareness, the present study aimed to forecast the cholera epidemic, explicitly addressing the reporting delay and ascertainment bias.

METHODS

Using weekly incidence of suspected cases, updated as a revised epidemic curve every week, the reporting delay was explicitly incorporated into the estimation model. Using the weekly case fatality risk as calculated by the World Health Organization, ascertainment bias was adjusted, enabling us to parameterize the family of logistic curves (i.e., logistic and generalized logistic models) for describing the unbiased incidence in 2017.

RESULTS

The cumulative incidence at the end of the epidemic, was estimated at 790,778 (95% CI: 700,495, 914,442) cases and 767,029 (95% CI: 690,877, 871,671) cases, respectively, by using logistic and generalized logistic models. It was also estimated that we have just passed through the epidemic peak by week 26, 2017. From week 27 onwards, the weekly incidence was predicted to decrease.

CONCLUSIONS

Cholera epidemic in Yemen, 2017 was predicted to soon start to decrease. If the weekly incidence is reported in the up-to-the-minute manner and updated in later weeks, not a single data point but the entire epidemic curve must be precisely updated.

Modeling Key Drivers of Cholera Transmission Dynamics Provides New Perspectives for Parasitology

Hydroclimatological and anthropogenic factors are key drivers of waterborne disease transmission. Information on human settlements and host mobility on waterways along which pathogens and hosts disperse, and relevant hydroclimatological processes, can be acquired remotely and included in spatially explicit mathematical models of disease transmission. In the case of epidemic cholera, such models allowed the description of complex disease patterns and provided insight into the course of ongoing epidemics. The inclusion of spatial information in models of disease transmission can aid in emergency management and the assessment of alternative interventions. Here, we review the study of drivers of transmission via spatially explicit approaches and argue that, because many parasitic waterborne diseases share the same drivers as cholera, similar principles may apply.

Big-data-driven modeling unveils country-wide drivers of endemic schistosomiasis

Schistosomiasis is a parasitic infection that is widespread in sub-Saharan Africa, where it represents a major health problem. We study the drivers of its geographical distribution in Senegal via a spatially explicit network model accounting for epidemiological dynamics driven by local socioeconomic and environmental conditions, and human mobility. The model is parameterized by tapping several available geodatabases and a large dataset of mobile phone traces. It reliably reproduces the observed spatial patterns of regional schistosomiasis prevalence throughout the country, provided that spatial heterogeneity and human mobility are suitably accounted for. Specifically, a fine-grained description of the socioeconomic and environmental heterogeneities involved in local disease transmission is crucial to capturing the spatial variability of disease prevalence, while the inclusion of human mobility significantly improves the explanatory power of the model. Concerning human movement, we find that moderate mobility may reduce disease prevalence, whereas either high or low mobility may result in increased prevalence of infection. The effects of control strategies based on exposure and contamination reduction via improved access to safe water or educational campaigns are also analyzed. To our knowledge, this represents the first application of an integrative schistosomiasis transmission model at a whole-country scale.

Modeling Biphasic Environmental Decay of Pathogens and Implications for Risk Analysis

As the appreciation for the importance of the environment in infectious disease transmission has grown, so too has interest in pathogen fate and transport. Fate has been traditionally described by simple exponential decay, but there is increasing recognition that some pathogens demonstrate a biphasic pattern of decay-fast followed by slow. While many have attributed this behavior to population heterogeneity, we demonstrate that biphasic dynamics can arise through a number of plausible mechanisms. We examine the identifiability of a general model encompassing three such mechanisms: population heterogeneity, hardening off, and the existence of viable-but-not-culturable states. Although the models are not fully identifiable from longitudinal sampling studies of pathogen concentrations, we use a differential algebra approach to determine identifiable parameter combinations. Through case studies using Cryptosporidium and Escherichia coli, we show that failure to consider biphasic pathogen dynamics can lead to substantial under- or overestimation of disease risks and pathogen concentrations, depending on the context. More reliable models for environmental hazards and human health risks are possible with an improved understanding of the conditions in which biphasic die-off is expected. Understanding the mechanisms of pathogen decay will ultimately enhance our control efforts to mitigate exposure to environmental contamination.

Mobile phone data highlights the role of mass gatherings in the spreading of cholera outbreaks

The spatiotemporal evolution of human mobility and the related fluctuations of population density are known to be key drivers of the dynamics of infectious disease outbreaks. These factors are particularly relevant in the case of mass gatherings, which may act as hotspots of disease transmission and spread. Understanding these dynamics, however, is usually limited by the lack of accurate data, especially in developing countries. Mobile phone call data provide a new, first-order source of information that allows the tracking of the evolution of mobility fluxes with high resolution in space and time. Here, we analyze a dataset of mobile phone records of ∼150,000 users in Senegal to extract human mobility fluxes and directly incorporate them into a spatially explicit, dynamic epidemiological framework. Our model, which also takes into account other drivers of disease transmission such as rainfall, is applied to the 2005 cholera outbreak in Senegal, which totaled more than 30,000 reported cases. Our findings highlight the major influence that a mass gathering, which took place during the initial phase of the outbreak, had on the course of the epidemic. Such an effect could not be explained by classic, static approaches describing human mobility. Model results also show how concentrated efforts toward disease control in a transmission hotspot could have an important effect on the large-scale progression of an outbreak.

Disease and disaster: Optimal deployment of epidemic control facilities in a spatially heterogeneous population with changing behaviour

Epidemics of water-borne infections often follow natural disasters and extreme weather events that disrupt water management processes. The impact of such epidemics may be reduced by deployment of transmission control facilities such as clinics or decontamination plants. Here we use a relatively simple mathematical model to examine how demographic and environmental heterogeneities, population behaviour, and behavioural change in response to the provision of facilities, combine to determine the optimal configurations of limited numbers of facilities to reduce epidemic size, and endemic prevalence. We show that, if the presence of control facilities does not affect behaviour, a good general rule for responsive deployment to minimise epidemic size is to place them in exactly the locations where they will directly benefit the most people. However, if infected people change their behaviour to seek out treatment then the deployment of facilities offering treatment can lead to complex effects that are difficult to foresee. So careful mathematical analysis is the only way to get a handle on the optimal deployment. Behavioural changes in response to control facilities can also lead to critical facility numbers at which there is a radical change in the optimal configuration. So sequential improvement of a control strategy by adding facilities to an existing optimal configuration does not always produce another optimal configuration. We also show that the pre-emptive deployment of control facilities has conflicting effects. The configurations that minimise endemic prevalence are very different to those that minimise epidemic size. So cost-benefit analysis of strategies to manage endemic prevalence must factor in the frequency of extreme weather events and natural disasters.