PREDICTIVE MODELING OF CHOLERA OUTBREAKS IN BANGLADESH

Predictive modeling for infectious diarrheal disease in pediatric populations: A systematic review

Introduction

Diarrhea is still a significant global public health problem. There are currently no systematic evaluation of the modeling areas and approaches to predict diarrheal illness outcomes. This paper reviews existing research efforts in predictive modeling of infectious diarrheal illness in pediatric populations.

Methods

We conducted a systematic review via a PubMed search for the period 1990-2021. A comprehensive search query was developed through an iterative process and literature on predictive modeling of diarrhea was retrieved. The following filters were applied to the search results: human subjects, English language, and children (birth to 18 years). We carried out a narrative synthesis of the included publications.

Results

Our literature search returned 2671 articles. After manual evaluation, 38 of these articles were included in this review. The most common research topic among the studies were disease forecasts 14 (36.8%), vaccine-related predictions 9 (23.7%), and disease/pathogen detection 5 (13.2%). Majority of these studies were published between 2011 and 2020, 28 (73.7%). The most common technique used in the modeling was machine learning 12 (31.6%) with various algorithms used for the prediction tasks. With change in the landscape of diarrheal etiology after rotavirus vaccine introduction, many open areas (disease forecasts, disease detection, and strain dynamics) remain for pathogen-specific predictive models among etiological agents that have emerged as important. Additionally, the outcomes of diarrheal illness remain under researched. We also observed lack of consistency in the reporting of results of prediction models despite the available guidelines highlighting the need for common data standards and adherence to guidelines on reporting of predictive models for biomedical research.

Conclusions

Our review identified knowledge gaps and opportunities in predictive modeling for diarrheal illness, and limitations in existing attempts whilst advancing some precursory thoughts on how to address them, aiming to invigorate future research efforts in this sphere.

Effect of cross-immunity in a two-strain cholera model with aquatic component

The bacteria Vibrio cholerae relies heavily upon an aquatic reservoir as a transmission route with two distinct serotypes observed in many recent outbreaks. In this paper, we extend previously studied ordinary differential equation epidemiological models to create a two-strain SIRP (susceptible-infectious-recovered-pathogen) system which incorporates both partial cross-immunity between disease strains and environmental pathogen transmission. Of particular interest are undamped anti-phase periodic solutions, as these display a type of coexistence where strains routinely switch dominance, and understanding what drives this switch can optimize the efficiency of the host population's control measures against the disease. We derive the basic reproduction number R0 and use stability analysis to examine the disease free and single-strain equilibria. We formulate a unique coexistence equilibrium and prove uniform persistence of both strains when R0>1. In addition, we simulate solutions to this system, along with seasonally forced versions of the model with and without host coinfection. Cross-immunity and transmission pathways influence damped or sustained oscillatory dynamics, where the presence of seasonality can modify, amplify or synchronize the period and phase of serotypes, driving epidemic waves. Cycling of serotypes over large time intervals, similar to observed data, is found for a range of cross-immunity levels, and the inclusion of coinfection in the model contributes to sustained anti-phase periodic solutions.

A linear noise approximation for stochastic epidemic models fit to partially observed incidence counts

Stochastic epidemic models (SEMs) fit to incidence data are critical to elucidating outbreak dynamics, shaping response strategies, and preparing for future epidemics. SEMs typically represent counts of individuals in discrete infection states using Markov jump processes (MJPs), but are computationally challenging as imperfect surveillance, lack of subject-level information, and temporal coarseness of the data obscure the true epidemic. Analytic integration over the latent epidemic process is impossible, and integration via Markov chain Monte Carlo (MCMC) is cumbersome due to the dimensionality and discreteness of the latent state space. Simulation-based computational approaches can address the intractability of the MJP likelihood, but are numerically fragile and prohibitively expensive for complex models. A linear noise approximation (LNA) that approximates the MJP transition density with a Gaussian density has been explored for analyzing prevalence data in large-population settings, but requires modification for analyzing incidence counts without assuming that the data are normally distributed. We demonstrate how to reparameterize SEMs to appropriately analyze incidence data, and fold the LNA into a data augmentation MCMC framework that outperforms deterministic methods, statistically, and simulation-based methods, computationally. Our framework is computationally robust when the model dynamics are complex and applies to a broad class of SEMs. We evaluate our method in simulations that reflect Ebola, influenza, and SARS-CoV-2 dynamics, and apply our method to national surveillance counts from the 2013-2015 West Africa Ebola outbreak.

A Review of the Environmental Trigger and Transmission Components for Prediction of Cholera

Climate variables influence the occurrence, growth, and distribution of Vibrio cholerae in the aquatic environment. Together with socio-economic factors, these variables affect the incidence and intensity of cholera outbreaks. The current pandemic of cholera began in the 1960s, and millions of cholera cases are reported each year globally. Hence, cholera remains a significant health challenge, notably where human vulnerability intersects with changes in hydrological and environmental processes. Cholera outbreaks may be epidemic or endemic, the mode of which is governed by trigger and transmission components that control the outbreak and spread of the disease, respectively. Traditional cholera risk assessment models, namely compartmental susceptible-exposed-infected-recovered (SEIR) type models, have been used to determine the predictive spread of cholera through the fecal–oral route in human populations. However, these models often fail to capture modes of infection via indirect routes, such as pathogen movement in the environment and heterogeneities relevant to disease transmission. Conversely, other models that rely solely on variability of selected environmental factors (i.e., examine only triggers) have accomplished real-time outbreak prediction but fail to capture the transmission of cholera within impacted populations. Since the mode of cholera outbreaks can transition from epidemic to endemic, a comprehensive transmission model is needed to achieve timely and reliable prediction with respect to quantitative environmental risk. Here, we discuss progression of the trigger module associated with both epidemic and endemic cholera, in the context of the autochthonous aquatic nature of the causative agent of cholera, V. cholerae, as well as disease prediction.

Statistical modeling of computer malware propagation dynamics in cyberspace

Modeling cyber threats, such as the computer malicious software (malware) propagation dynamics in cyberspace, is an important research problem because models can deepen our understanding of dynamical cyber threats. In this paper, we study the statistical modeling of the macro-level evolution of dynamical cyber attacks. Specifically, we propose a Bayesian structural time series approach for modeling the computer malware propagation dynamics in cyberspace. Our model not only possesses the parsimony property (i.e. using few model parameters) but also can provide the predictive distribution of the dynamics by accommodating uncertainty. Our simulation study shows that the proposed model can fit and predict the computer malware propagation dynamics accurately, without requiring to know the information about the underlying attack-defense interaction mechanism and the underlying network topology. We use the model to study the propagation of two particular kinds of computer malware, namely the Conficker and Code Red worms, and show that our model has very satisfactory fitting and prediction accuracies.

Methods for Model Calibration under High Uncertainty: Modeling Cholera in Bangladesh

Background. Published data on a disease do not always correspond directly to the parameters needed to simulate natural history. Several calibration methods have been applied to computer-based disease models to extract needed parameters that make a model's output consistent with available data. Objective. To assess 3 calibration methods and evaluate their performance in a real-world application. Methods. We calibrated a model of cholera natural history in Bangladesh, where a lack of active surveillance biases available data. We built a cohort state-transition cholera natural history model that includes case hospitalization to reflect the passive surveillance data-generating process. We applied 3 calibration techniques: incremental mixture importance sampling, sampling importance resampling, and random search with rejection sampling. We adapted these techniques to the context of wide prior uncertainty and many degrees of freedom. We evaluated the resulting posterior parameter distributions using a range of metrics and compared predicted cholera burden estimates. Results. All 3 calibration techniques produced posterior distributions with a higher likelihood and better fit to calibration targets as compared with prior distributions. Incremental mixture importance sampling resulted in the highest likelihood and largest number of unique parameter sets to better inform joint parameter uncertainty. Compared with naïve uncalibrated parameter sets, calibrated models of cholera in Bangladesh project substantially more cases, many of which are not detected by passive surveillance, and fewer deaths. Limitations. Calibration cannot completely overcome poor data quality, which can leave some parameters less well informed than others. Calibration techniques may perform differently under different circumstances. Conclusions. Incremental mixture importance sampling, when adapted to the context of high uncertainty, performs well. By accounting for biases in data, calibration can improve model projections of disease burden.

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.

PHENOMENOLOGICAL FORECASTING OF DISEASE INCIDENCE USING HETEROSKEDASTIC GAUSSIAN PROCESSES: A DENGUE CASE STUDY

In 2015 the US federal government sponsored a dengue forecasting competition using historical case data from Iquitos, Peru and San Juan, Puerto Rico. Competitors were evaluated on several aspects of out-of-sample forecasts including the targets of peak week, peak incidence during that week, and total season incidence across each of several seasons. our team was one of the winners of that competition, outperforming other teams in multiple targets/locales. In this paper we report on our methodology, a large component of which, surprisingly, ignores the known biology of epidemics at large-for example, relationships between dengue transmission and environmental factors-and instead relies on flexible nonparametric nonlinear Gaussian process (GP) regression fits that "memorize" the trajectories of past seasons, and then "match" the dynamics of the unfolding season to past ones in real-time. Our phenomenological approach has advantages in situations where disease dynamics are less well understood, or where measurements and forecasts of ancillary covariates like precipitation are unavailable, and/or where the strength of association with cases are as yet unknown. In particular, we show that the GP approach generally outperforms a more classical generalized linear (autoregressive) model (GLM) that we developed to utilize abundant covariate information. We illustrate variations of our method(s) on the two benchmark locales alongside a full summary of results submitted by other contest competitors.

Efficient Data Augmentation for Fitting Stochastic Epidemic Models to Prevalence Data

Stochastic epidemic models describe the dynamics of an epidemic as a disease spreads through a population. Typically, only a fraction of cases are observed at a set of discrete times. The absence of complete information about the time evolution of an epidemic gives rise to a complicated latent variable problem in which the state space size of the epidemic grows large as the population size increases. This makes analytically integrating over the missing data infeasible for populations of even moderate size. We present a data augmentation Markov chain Monte Carlo (MCMC) framework for Bayesian estimation of stochastic epidemic model parameters, in which measurements are augmented with subject-level disease histories. In our MCMC algorithm, we propose each new subject-level path, conditional on the data, using a time-inhomogeneous continuous-time Markov process with rates determined by the infection histories of other individuals. The method is general, and may be applied to a broad class of epidemic models with only minimal modifications to the model dynamics and/or emission distribution. We present our algorithm in the context of multiple stochastic epidemic models in which the data are binomially sampled prevalence counts, and apply our method to data from an outbreak of influenza in a British boarding school.

Improving immunization approaches to cholera

INTRODUCTION

Cholera's impact is greatest in resource-limited countries. In the last decade several large epidemics have led to a global push to improve and implement the tools for cholera prevention and control. Areas covered: PubMed, Google Scholar and the WHO website were searched to review the literature and summarize the current status of cholera vaccines to make recommendations on improving immunization approaches to cholera. Oral cholera vaccines (OCVs) have demonstrated their effectiveness in endemic, outbreak response and emergency settings, highlighting their potential for wider adoption. While two doses of the currently available OCVs are recommended by manufacturers, a single dose would be easier to implement. Encouragingly, recent studies have shown that cold chain requirements may no longer be essential. The establishment of the global OCV stockpile in 2013 has been a major advance in cholera preparedness. New killed and live-attenuated vaccines are being actively explored as candidate vaccines for endemic settings and/or as a traveller's vaccine. The recent advances in cholera vaccination approaches should be considered in the global cholera control strategy. Expert commentary: The development of affordable cholera vaccines is a major success to improve cholera control. New vaccines and country specific interventions will further reduce the burden of this disease globally.

PREDICTIVE MODELING OF CHOLERA OUTBREAKS IN BANGLADESH

Despite seasonal cholera outbreaks in Bangladesh, little is known about the relationship between environmental conditions and cholera cases. We seek to develop a predictive model for cholera outbreaks in Bangladesh based on environmental predictors. To do this, we estimate the contribution of environmental variables, such as water depth and water temperature, to cholera outbreaks in the context of a disease transmission model. We implement a method which simultaneously accounts for disease dynamics and environmental variables in a Susceptible-Infected-Recovered-Susceptible (SIRS) model. The entire system is treated as a continuous-time hidden Markov model, where the hidden Markov states are the numbers of people who are susceptible, infected, or recovered at each time point, and the observed states are the numbers of cholera cases reported. We use a Bayesian framework to fit this hidden SIRS model, implementing particle Markov chain Monte Carlo methods to sample from the posterior distribution of the environmental and transmission parameters given the observed data. We test this method using both simulation and data from Mathbaria, Bangladesh. Parameter estimates are used to make short-term predictions that capture the formation and decline of epidemic peaks. We demonstrate that our model can successfully predict an increase in the number of infected individuals in the population weeks before the observed number of cholera cases increases, which could allow for early notification of an epidemic and timely allocation of resources.