Direct transmission via households informs models of disease and intervention dynamics in cholera

Abstract

While several basic properties of cholera outbreaks are common to most settings—the pathophysiology of the disease, the waterborne nature of transmission, and others—recent findings suggest that transmission within households may play a larger role in cholera outbreaks than previously appreciated. Important features of cholera outbreaks have long been effectively modeled with mathematical and computational approaches, but little is known about how variation in direct transmission via households may influence epidemic dynamics. In this study, we construct a mathematical model of cholera that incorporates transmission within and between households. We observe that variation in the magnitude of household transmission changes multiple features of disease dynamics, including the severity and duration of outbreaks. Strikingly, we observe that household transmission influences the effectiveness of possible public health interventions (e.g. water treatment, antibiotics, vaccines). We find that vaccine interventions are more effective than water treatment or antibiotic administration when direct household transmission is present. Summarizing, we position these results within the landscape of existing models of cholera, and speculate on its implications for epidemiology and public health.

About the Authors:

Victor A. Meszaros

Contributed equally to this work with: Victor A. Meszaros, Miles D. Miller-Dickson

Roles Formal analysis, Investigation, Methodology, Validation, Visualization, Writing – original draft, Writing – review & editing

Affiliation: Department of Ecology and Evolutionary Biology, Brown University, Providence, RI, United States of America

Miles D. Miller-Dickson

Contributed equally to this work with: Victor A. Meszaros, Miles D. Miller-Dickson

Roles Formal analysis, Investigation, Methodology, Validation, Visualization, Writing – original draft, Writing – review & editing

Affiliation: Department of Ecology and Evolutionary Biology, Brown University, Providence, RI, United States of America

Francis Baffour-Awuah Junior

Roles Formal analysis, Methodology, Writing – original draft, Writing – review & editing

Affiliation: Department of Mathematics, Florida State University, Tallahassee, FL, United States of America

Salvador Almagro-Moreno

Roles Conceptualization, Investigation, Supervision, Writing – original draft, Writing – review & editing

Affiliations Burnett School of Biomedical Sciences, University of Central Florida, Orlando, FL, United States of America, National Center for Integrated Coastal Research, University of Central Florida, Orlando, FL, United States of America

ORCID logo http://orcid.org/0000-0002-1121-1204

C. Brandon Ogbunugafor

Roles Conceptualization, Formal analysis, Funding acquisition, Investigation, Methodology, Project administration, Resources, Supervision, Validation, Writing – original draft, Writing – review & editing

* E-mail: [email protected]

Affiliation: Department of Ecology and Evolutionary Biology, Brown University, Providence, RI, United States of America

ORCID logo http://orcid.org/0000-0002-1581-8345

Introduction

Waterborne infections continue to pose a significant threat to human health. Among the most prominent of the waterborne infections is cholera, caused by the gram negative bacterium Vibrio cholerae, acquired through the consumption of contaminated water or food [1]. New epidemics are often driven by natural disasters or social conflict that facilitate the increased intake of contaminated water [2, 3]. As cholera is so strongly associated with the water supply, many public health intervention efforts have focused on the importance of water quality in preventing outbreaks. Recent studies have demonstrated that in certain settings, direct household transmission can serve an increasingly large role in disease in modern outbreaks [4–7]. We might especially expect direct household transmission to influence disease dynamics in high-density settings, like those associated with displaced populations (e.g. refugee camps) [8, 9].

The notion that cholera outbreaks can be driven or sustained by direct household transmission—which may occur between individuals within and between households—might alter the lens through which we study and potentially intervene in cholera outbreaks. Given how recently epidemiologists have engaged the specifics of direct household transmission, we have yet to fully appreciate the extent to which direct household transmission may alter the dynamics of cholera outbreaks. And while there are many important models of indirect cholera transmission via reservoirs [10–21], and some that examine the theoretical implications of direct household transmission [22, 23], we remain in the dark with regards to how recent insights on direct household transmission may specifically influence our picture of cholera disease dynamics and interventions.

In this study, we construct a mathematical model of cholera in a high-density setting, with random-mixing populations. Unlike prior studies, our model combines features of models that include waterborne transmission with parameters informed by recent empirical findings on direct household transmission. We observe that an outbreak in a high-density setting typified by a large amount of direct household transmission can resemble an outbreak driven exclusively by waterborne transmission, but with several critical differences, especially with regards to the effectiveness of interventions. Consequently, high-density cholera epidemics with direct household transmission suggest different intervention strategies in order to meaningfully suppress transmission. Our findings emphasize the need for more data on direct household transmission in high-density settings, as it can have important consequences on disease dynamics and optimal intervention measures.

Model of disease dynamics with household transmission

Description

We present a dynamical model of cholera based on a series of ordinary differential equations (ODEs). The model incorporates dynamics of both direct household transmission and disease flow through water reservoirs, as with other mathematical models of waterborne infections [22, 24]. We’ve codified our iteration as the Waterborne, Abiotic and other Indirectly Transmitted or (W.A.I.T.) framework, and will refer to it as such throughout this study.

Our model structure is adapted from a mathematical model of the 2010 outbreak in Haiti [10]. The model describes a cholera epidemic outbreak and is set with a population of 500000 (approximating the number of individuals in a mid-sized city in many countries affected by cholera).

In order to interrogate the under-examined role of direct household transmission within a cholera epidemic, the model allows us to investigate the impact of direct transmission of the susceptible population as a proxy for transmission within and between households. Throughout the text, we refer to all non-waterborne, environmental transmission in this model as “direct,” “household,” or “direct household” transmission. We use the terms interchangeably. By constructing a framework which allows us to combine different routes of transmission, we are afforded the opportunity to observe how disease burden and ideal intervention strategies are influenced by different transmission strategies.(1)(2)(3)(4)(5)(6)(7)

The compartmental diagram

The model can be visualized using a modified compartmental diagram that contains information on the spread of disease agents and hosts through the system, as seen in Fig 1. The model features seven ordinary differential equations—one for each compartment—and twenty-two parameters. The compartments S, I, A, R and V represent populations of susceptible, symptomatically infected, asymptomatically infected, recovered and vaccinated individuals respectively. The WH and WL compartments represent total V. cholerae bacteria, of each type (high or low infectious) within a single water reservoir accessible to the sum total of individuals within the population.

[Figure omitted. See PDF.]

Fig 1. A Waterborne, Abiotic and other Indirectly Transmitted (W.A.I.T.) dynamical compartmental model of cholera.

Top: This diagram is a variation of a standard compartmental diagram. Green arrows highlight the flow of the population of hosts through the system. Bottom: here the red arrows highlight flow of disease through the system.

https://doi.org/10.1371/journal.pone.0229837.g001

Some model features, which have also reflected findings in prior studies [10, 25, 26], include the following:

* Disease is acquired by susceptible individuals S, via infected drinking water from either of the two reservoirs, or through direct household transmission. The fraction of S-individuals that move either to A or I is then given by p and 1 − p respectively

* The reservoir compartments WH and WL differ in how much V. cholerae is required from each to infect an individual, with the high-infectious reservoir (WH) having a lower minimum infectious dose. With time, populations of V. cholerae within the high-infectious reservoir naturally decay to a population with lower infectiousness.

* There are two groups of infected individuals: I-infected individuals have significantly higher excretion rates and can be treated with antibiotics, decreasing their excretion rate. A-infected individuals have lower natural excretion rates.

* Individuals enter and exit the system via a constant natural birth rate as well as infected and uninfected death rates.

* Infection is removed from the system either via bacterial decay in the low-infectious reservoir, or via the death of infected individuals.

* The number of S-individuals that are vaccinated is modeled as a fraction of the total population size, and individuals return to the susceptible compartment as their vaccine induced immunity wanes.

* I and A-individuals recover from the disease at a natural rate, I-individuals treated with antibiotics recover at a faster rate, and the recovered state provides a temporary immunity which, after waning, will feed individuals back into the susceptible population.

* Importantly, the direct transmission parameters are based on an empirical study of household transmission in Bangladesh [7].

Analytic equations and parameters

The set of ordinary differential equations (Eqs 1–7) define the dynamics of the system. Direct household transmission is implemented through the term η. This term translates to the total daily fraction of all possible interactions between susceptible individuals and infected individuals (symptomatic and asymptomatic) that lead to an infection. As stated previously, household transmission is implemented using an empirical study that estimated rates of transmission of cholera via households in Bangladesh [7].

In our modeling scheme (W.A.I.T.), environmental dynamics are also realized within their own set of differential equations. The terms quantify the likelihood that individuals become infected with V. cholerae given a daily water consumption rate α and a minimum infectious dose κ for each reservoir.

For either route of transmission (household or waterborne) the parameter p scales the rate of transmission by serving as the fraction of individuals that move to the A-compartment (asymptomatic), whereas 1 − p serves as the fraction that move to the I-compartment (symptomatic). p is, therefore, a population-level expression of host factors.

Other notable terms include τ, the vaccination rate, which represents the daily fraction of S-individuals that receive vaccinations; ξs, the symptomatic individuals excretion rate, which is three orders of magnitude larger than that of the asymptomatic case and can be reduced via antibiotic administration; λ, the relative recovery rate of those receiving antibiotics; and θ, the proportion of symptomatic individuals receiving antibiotics. Table 1 lists the parameters and values used in this simulation.

[Figure omitted. See PDF.]

Table 1. Cholera household (direct) transmission model parameters.

https://doi.org/10.1371/journal.pone.0229837.t001

The initial conditions used across all simulations included a population of 500 000 susceptible individuals with a small initial amount of Vibrio cholerae present in each reservoir to prime the outbreak (5 cells at t = 0). Additionally, the values of the symptomatic, asymptomatic, vaccinated, and recovered compartments were set to zero at the start of each run.

Parameter sensitivity analysis and basic reproductive ratio

The basic reproductive number, R0 is a contentious measure in epidemiology, but can be interpreted as a signature for the average infectiousness of a pathogen in a given setting [27, 28]. In this model, given the nominal parameters, we find that R0 = 8.37, this value is in agreement with reproductive ratio ranges (R0 = 3–19) determined from past epidemics and computational models of cholera [19, 25]. Estimates for the reproductive ratio of cholera vary within the literature, which correspond to different epidemic settings and contexts. In the S1 File, we outline how the numerical value for R0 was computed in the model presented in this study.

The relative sensitivity of R0 to changes in the model parameters was examined. Both an independent parameter sensitivity analysis as well as a partial rank correlation coefficient (PRCC) analysis [45] were conducted. The PRCC was performed at 500 iterations with 50 samples per iteration. Fig 2 shows the results of these investigations. Here we find that p, κL, W and α have the greatest independent (uncorrelated) sensitivity, while p, π, W, and α hold the greatest interrelated sensitivity by mean value. These parameters include two parameters strongly related to environmental aspects of the model: the rate of water consumption α, and the total water reservoir size W. Both the independent and the correlation sensitivity analyses are in agreement with regards to the sensitivity of the model parameters on the R0.

[Figure omitted. See PDF.]

Fig 2. R0 sensitivity analysis.

Left: A tornado plot showing the independent sensitivity of R0 to changes in the model parameters. Black bars indicate the value of R0 when the associated parameter is increased by 25% from its nominal value. White bars indicate the value of R0 when the associated parameter is decreased by 25%. Right: A Partial Rank Correlation Coefficient (PRCC) analysis was preformed with respect to R0 highlighting the inter-correlated sensitivities of each of the model parameters. The blue bars show the mean value of each PRCC, with error bars at one standard deviation. This analysis was preformed by sampling over uniform distributions of 15% around the nominal model parameter values in line with the methods provided by Blower and Dowlatabadi (1994) [45].

https://doi.org/10.1371/journal.pone.0229837.g002

Of note, p (proportion of asymptomatic cases) is the most sensitive parameter within the model by a large margin. The fraction of infected individuals that become asymptomatic can be considered as an average characteristic of the model population. By our assessment, the high sensitivity in R0 due to p is caused by the relatively large distinction between the excretion rates of symptomatic (ξS) and asymptomatic (ξA) individuals (ξS/ξA ∼ 103). Given this, if the fraction of individuals that move into the asymptomatic compartment, p, is increased, the basic reproductive ratio is noticeably reduced.

At the disease-free equilibrium—the regime where R0 is calculated—the direct transmission term η shows low sensitivity to the basic reproductive number, indicating that variations in the household transmission rate may not lead to significant variations in the R0 value. Thus, whereas the rate η may have a noticeable impact on the dynamical peak of the infection, or on the time scales associated with the outbreak of the infection, it may not largely impact how explosive (in terms of R0) the contagion is initially. More generally, this can be interpreted as a limitation of the R0 metric with regards to the epidemic features that it captures: it is not all-encompassing. In addition to R0, a full and proper understanding of the epidemic dynamics requires a thorough analysis of model-specific features (e.g. the behavior of the peak of the infection).

The infected populations in our model (both asymptomatic and symptomatic) initially increase—moving away from the disease-free equilibrium—and eventually reach stable non-zero (although small) equilibrium values (see Fig 3(a)). This is consistent with R0 > 1, as well as the eigenvalue analysis of the Jacobian matrix at the disease-free equilibrium (shown in the Supplemental Appendix). The eigenvalues of the Jacobian matrix for the ODE system at the disease-free equilibrium are all real numbers, with two values greater than zero, indicating that the disease-free equilibrium is unstable. We present the Jacobian matrix of the system and its corresponding eigenvalues in the Supplemental Appendix.

[Figure omitted. See PDF.]

Fig 3. Basic model dynamics as a function of household transmission and water treatment: Of the (a) infected and (b) recovered populations, as well as the (c) combined bacterial load of the water reservoirs, normalized to maximum bacterial load value across all 3 cases.

Lines correspond to different epidemic conditions, including a no household transmission scenario, one with household transmission and one where household transmission is coupled with a water treatment intervention (which decreases the amount of transmission through consumption of contaminated water).

https://doi.org/10.1371/journal.pone.0229837.g003

Model dynamics: The added effects of direct household transmission

Introduction of infection between individuals within and between households into models of cholera has a substantial impact on the dynamics of outbreaks. In order to showcase this impact, three scenarios were investigated with simulations: (i) No household transmission (η = 0), (ii) household transmission with standard water consumption (η & α = nominal values), and (iii) household transmission with water treatment (decreased consumption of contaminated water) (η = nominal value and 0.5α). Simulations were run at a time resolution of 10,000 timesteps over a period of 60 days for all three cases. Figs 3–5 depict various aspects of the dynamics of disease as well as the effects of potential interventions in each scenario. Fig 3(a) demonstrates a 5.63x increase in the peak number of total infected individuals when realistic household transmission rates are introduced [7]. Fig 3(b) shows the increased peak values of the total number of recovering individuals, as more agents are drawn into the recovered compartments due to the increase in the population of infected individuals. This correlation is an example of how household transmission increases the intensity of the epidemic within the model. The total flow of individuals moving through the system increases in correspondence with an increase in household transmission. In the Supplemental Appendix, we provide information on how the different transmission scenarios influence the total number of infections across the duration of the model.

[Figure omitted. See PDF.]

Fig 4. Peak infected population behavior: A heat map of peak total infected population as a function of household (direct) transmission terms (1 − p)η and pη.

Along the X and Y axis, the household transmission terms are notated as multiples of their nominal model values as seen in Table 1. The top left corner denotes the region where pη = (1 − p)η = 0, and no household transmission occurs. While the bottom right corner denotes the region where η = nominal model value. The gray-scale intensity within the heat map represents the peak total infected population. Comparing the intensity of the top left corner to the bottom right we see an increase in the peak total infected population as household transmission is introduced and its effects increased.

https://doi.org/10.1371/journal.pone.0229837.g004

[Figure omitted. See PDF.]

Fig 5. Comparing the effectiveness of interventions in the setting of household transmission.

There are three hypothetical intervention scenarios, and their effect on the dynamics of the total population of infected individuals in the cases when household (direct) transmission is (a) absent and (b) present. The primary Y-axes of fig(a) and (b) are identical. The inset in fig(a) highlights the regions, x = 0 to 40 (days) and y = 0 to 6 ⋅ 103 in order to better distinguish the differences in the dynamics of the curves.

https://doi.org/10.1371/journal.pone.0229837.g005

Fig 4 highlights properties of the peak total infected population, and how it is modified by household transmission. While there are several features of an epidemic that one could use to capture this influence, we consider the peak number of infected individuals because it signifies the state where the infection is exacting its greatest toll. Each of the 250 000 grey-scale pixels that compose the heat-map represents a single run of the cholera household transmission model. For each run, we plot the peak of the total (combined I + A) infected population, modulating the rate of transmission into the A and I compartments via the direct transmission terms, pη and (1 − p)η respectively. All other model parameters are kept constant, held at their nominal values as shown in Table 1. The x and y axes show fractions of the nominal values of pη and (1 − p)η, respectively, from Table 1. At the top left corner where pη = (1 − p)η = 0, there is no household transmission present. The peak total infected population here is at its lowest value, 5647. This is equal to the same value we see at the peak of the solid curve in Fig 3(a). As we increase η and thus the amount of household transmission present, we see a dramatic increase in the peak number of total infected individuals (lower right corner).

The largest peak value that we observe, 31787, appears in the lower right corner of the plot where the household transmission terms are their largest. This value corresponds to the peak of the dashed-dotted line in Fig 3(a). Implementing direct household transmission in this model thus increases the peak total infected population by approximately 5.63x. Of note: p is expected to be an approximately fixed property of hosts, whereas direct household transmission is modulated solely by the η-term, representing varying levels of inter- and intra-household transmission. Thus, we only consider the peak infection values along the diagonal elements of the heat map as representing relevant (or realistic) scenarios. The rest of the space on the heat map showcases the extent to which host-level affinities for symptomatic or asymptomatic infection-types can affect the peak levels of infection.

Like several other models, the W.A.I.T. modeling scheme utilized in this study allows us to independently model the dynamics of the environmental reservoirs via a set of separate differential equations. In Fig 3(c), we observe that the normalized total reservoir bacterial load peak is 2.95x larger for the case where household transmission is present. Additionally, the peak bacterial load is shifted 0.63x of a day later in the simulation when household transmission is introduced. This “late arrival” phenomenon is also observed, within the infected and recovered populations. This indicates, within the scope of the model, that the overall duration of the epidemic is also slightly increased when household transmission is introduced.

Intervention schemes

Having explored the impact of household transmission on disease dynamics, we then analyzed a range of potential sample interventions: vaccination, antibiotic administration, and water treatment (decreasing contaminated water consumption). Each were realized by modifying model parameters relevant to that intervention. We compare the impact of these interventions in cases with and without household transmission. Note that these intervention schemes assume implementation occurs at the onset of the infection.

Fig 5 depicts the total number of both symptomatic and asymptomatic infected individuals in settings utilizing various intervention schemes. In the case where no direct household transmission is present (Fig 5a), decreasing the consumption of contaminated water (α) has a strong effect on disease dynamics: a 0.5x decrease in contaminated water consumption provides an about 0.5x decrease in the peak number of infected individuals. In a similar setting with no direct transmission, a 1.5x increase in vaccination and antibiotic interventions provide a smaller relative decrease in the peak number of infected individuals. This demonstrates that for cholera epidemics driven largely by the waterborne transmission, and not through household transmission, one would expect (unsurprisingly) that targeting the water consumption would be the best way to attenuate the epidemic.

Alternatively, our explorations of settings featuring high levels of household transmission (Fig 5b) show that vaccine interventions emerge as the most effective strategy. In settings featuring high levels of household transmission, we predict that emergency vaccination programs would be more effective than water treatment-based interventions. The reasoning behind this: vaccination directly modulates the fraction of person-to-person contacts that involve an infected individual by removing infected individuals from the population, and preventing an explosive outbreak from getting going. In settings with little household transmission, vaccination indirectly mitigates the spread of infection by removing infected individuals that could contribute to the transmission cycle by infecting the water reservoir.

Discussion

In this study, we propose a modeling framework for cholera in high-density, well-mixed population settings that fully encompasses two different routes of transmission: standard transmission via an environmental reservoir intermediate (waterborne), and directly between hosts in and between households. The question of how these two routes interact is important on the backdrop of relatively new findings that highlight the importance of household transmission [4–7]. Importantly, our model structure allows one to examine how different magnitudes of household transmission influence disease dynamics. A recent model of high-density settings simulated these dynamics effectively, but did so for a single refugee camp with a specific set of parameter values [9]. Our model, alternatively, explores the space of possible direct transmission parameters that may pertain to different high-density cholera settings.

We find that variation in household transmission changes how the disease is spread. Interestingly, the presence of direct transmission not only influences R0, but rather, it changes specific features of disease that cannot be captured by the R0. For example, variation in household transmission influences both the peak number of cases in an outbreak and the duration of a given outbreak. Whereas, the R0 speaks more directly to the “explosiveness” of an outbreak, how rapidly it progresses initially. This result has relevance for understanding specifics of cholera outbreaks and for general canon in mathematical epidemiology. With regards to cholera outbreaks, household transmission can perniciously influence disease outcomes by extending the severity and duration of a given outbreak. With regards to broader conversations in mathematical epidemiology, the household transmission case highlights the limits of R0 as a tool for assessing the severity of epidemics. While it remains a powerful proxy for certain characteristics of an outbreak, undue focus on it as the defining characteristic of an epidemic can be misleading.

As cholera remains a major global health problem in high-density, temporary settlements, our results suggest that the public health sector should be mindful of how sensitive outbreak dynamics are to different transmission characteristics. Specifically, the nature and amount of household transmission can have a meaningful impact, not only on the trajectory of disease, but on which intervention strategies might be most effective in preventing or attenuating a given epidemic. Our results support recent findings suggesting that vaccines can be especially effective in high-density settings with large amounts of household transmission [9]. That is, water treatment and antibiotics—both of which may be relevant in other settings—are likely to be less effective in subverting epidemics in refugee camps, and related settings where household transmission is high. One can explain this by highlighting the nature of high-density settings: transmission can occur too quickly (between individuals) for more structural-based (e.g. water treatment) or post-infection (e.g. antibiotics) interventions to have as meaningful an effect. Vaccines, however, can lower the peak of infection by decreasing the number of individuals entering a high-density setting, preventing both the number of individuals capable of direct transmission and the contribution of disease into the common water supply.

In closing, we urge epidemiologists to appreciate that epidemics are diverse ecological phenomena, where different settings may require peculiar and specific interventions. Our findings implore the research community to study epidemics with an increased level of granularity, as essential properties can be defined by particulars.

Supporting information

[Figure omitted. See PDF.]

S1 File. Supplemental appendix.

Here we describe a set of additional analyses and calculations relevant to content in the main text. These include:

* Calculation of the sum total of cases as a function of direct infection

* An outline of the method for determining the reproductive ratio (R0) for the Cholera household (direct) infection model.

* Discussion of the stability of the system at the disease free equilibrium.

https://doi.org/10.1371/journal.pone.0229837.s001

(PDF)

Acknowledgments

The authors would like to thank T.Togun, N.Harris, S.Robinson, A.Murillo, M.Rodriguez, C.Castillo-Chavez, L.Gomez, and S.Scarpino for various contributions on the topic. The authors would also like to thank 2 anonymous reviewers for their helpful feedback on a draft of the manuscript.

Citation: Meszaros VA, Miller-Dickson MD, Baffour-Awuah F Junior, Almagro-Moreno S, Ogbunugafor CB (2020) Direct transmission via households informs models of disease and intervention dynamics in cholera. PLoS ONE 15(3): e0229837. https://doi.org/10.1371/journal.pone.0229837

References

1. Clemens JD, Nair GB, Ahmed T, Qadri F, Holmgren J. Cholera. The Lancet. 2017 Sep;390(10101):1539–1549.

2. Watson JT, Gayer M, Connolly MA. Epidemics after natural disasters. Emerging infectious diseases. 2007 Jan;13(1):1. pmid:17370508

3. Dureab FA, Shibib K, Al-Yousufi R, Jahn A. Yemen: cholera outbreak and the ongoing armed conflict. The Journal of Infection in Developing Countries. 2018 May;12(05):397–403. pmid:31865306

4. Rahman Z, Rahman MD, Rashid M, Monira S, Johura FT, Mustafiz M, et al. Vibrio cholerae Transmits through Water among the Household Contacts of Cholera Patients in Cholera Endemic Coastal Villages of Bangladesh, 2015-2016 (CHoBI7 Trial). Frontiers in Public Health. 2018 Oct;6(1):1635.

5. Rafique R, Rashid M, Monira A, Rahman Z, Mahmud MD, Mustafiz M, et al. Transmission of infectious Vibrio cholerae through drinking water among the household contacts of cholera patients (CHoBI7 trial). Frontiers in Microbiology. 2016 Oct;7(1):1635. pmid:27803695

6. George CM, Hasan K, Monira S, Rahman Z, Saif-Ur-Rahman KM, Rashid M, et al. A prospective cohort study comparing household contact and water Vibrio cholerae isolates in households of cholera patients in rural Bangladesh. PLoS Neglected Tropical Diseases. 2018 Jul;12(7):e0006641. pmid:30052631

7. Sugimoto JD, Koepke AA, Kenah EE, Halloran ME, Chowdhury F, Khan AI, et al. Household Transmission of Vibrio cholerae in Bangladesh. PLOS Neglected Tropical Diseases. 2014 Nov;8(11):e3314. pmid:25411971

8. Kuhl ME, Crooks A, Hailegiorgis A. Disease modeling within refugee camps: A multi-agent systems approach. Proceedings of the 2013 Winter Simulation Conference. Simulation: Making Decisions in a Complex World. 2013 Dec:1697–1706.

9. Havumaki J, Meza R, Phares CR, Date K, Eisenberg MC. Comparing alternative cholera vaccination strategies in Maela refugee camp: using a transmission model in public health practice. BMC Infectious Diseases. 2019 Dec;19(1075):1–17.

10. Andrews JR, Basu S. Transmission dynamics and control of cholera in Haiti: an epidemic model. The Lancet. 2011 Apr;377(9773):1248–1255.

11. Codeço CT. Endemic and epidemic dynamics of cholera: the role of the aquatic reservoir. BMC Infectious Diseases. 2001 Feb;1(1):1–14. pmid:11208258

12. Bertuzzo E, Casagrandi R, Gatto M, Rodriguez-Iturbe I, Rinaldo A. On spatially explicit models of cholera epidemics. Journal of the Royal Society Interface. 2009 Jul;7(43):321–333.

13. Cholera epidemic in Haiti, 2010: using a transmission model to explain spatial spread of disease and identify optimal control interventions. Annals of internal medicine. 2011 May;154(9):593–601. pmid:21383314

14. Mari L, Bertuzzo E, Righetto L, Casagrandi R, Gatto M, Rodriguez-Iturbe I, et al. Modelling cholera epidemics: the role of waterways, human mobility and sanitation. Journal of the Royal Society Interface. 2012 Feb;9(67):376–388.

15. Eisenberg MC, Shuai Z, Tien JH, Van den Driessche P. A cholera model in a patchy environment with water and human movement. Mathematical Biosciences. 2013 Nov;246(1):105–112. pmid:23958383

16. Sun GQ, Xie JH, Huang SH, Jin Z, Li MT, Liu L. Transmission dynamics of cholera: Mathematical modeling and control strategies. Communications in Nonlinear Science and Numerical Simulation. 2017 Apr;45(1):235–244.

17. Akman O, Schaefer E. An evolutionary computing approach for parameter estimation investigation of a model for cholera. Journal of Biological Dynamics. 2016 May;9(1):147–158.

18. Akman O, Romadan MC, Schaefer E. Examination of models for cholera: insights into model comparison methods. Letters in Biomathematics. 2016 Jul;3(1):93–118.

19. Neilan RM, Gaff H, Schaefer E, Fister KR. Modeling Optimal Intervention Strategies for Cholera. Bulletin of Mathematical Biology. 2010 Mar;72(8):2004–2018.

20. Mwasa A, Tchuenche JM. Mathematical analysis of a cholera model with public health interventions. Journal of Biosystems. 2011 Apr;105(3):190–200. pmid:21540075

21. Al-Arydah M, Mwasa A, Jean TM, Smith RJ. Modeling Cholera Disease With Education And Chlorination. Journal of Biological Systems. 2013 Dec;21(4):20.

22. Tien JH, Earn D. Multiple transmission pathways and disease dynamics in a waterborne pathogen model. Bulletin of Mathematical Biology. 2010 Aug;72(6):1506–1533. pmid:20143271

23. Shuai Z, Van den Driessche P. Global dynamics of cholera models with differential infectivity. Mathematical Biosciences. 2011 Dec;234(2):118–126. pmid:22001141

24. Li M, Ma J, van den Driessche P. Model for disease dynamics of a waterborne pathogen on a random network. Journal of Mathematical Biology. 2015 Oct;71(4):961–977. pmid:25326654

25. Hartley DM, Morris JG Jr, Smith DL. Hyperinfectivity: A Critical Element in the Ability of V. cholerae to Cause Epidemics? PLOS Medicine. 2005 Dec;3(1):e7. pmid:16318414

26. Department of Economic and Social Affairs, Population Division. World Population Prospects The 2008 Revision, Highlights, Working Paper. The United Nations. 2008.

27. the United Nations Children’s Fund (UNICEF). Ministère De La Santédu Cameroun, Données Choléra UNICEF. 2014.

28. Rahaman MM, Majid MA, Alam AK, Islam MR. Effects of doxycycline in actively purging cholera patients: A double-blind clinical trial. American Society for Microbiology. 1976 Oct;10(4):610–612.

29. Kolaye G, Bowong S, Houe R, Aziz-Alaoui M, Cadivel M. Mathematical Assessment of the Role of Environmental Factors on the Dynamical transmission of Cholera. Communications in Nonlinear Science and Numerical Simulation. 2019 Feb;67(1):203–222.

30. Benenson AS, Mosley WH, Fahimuddin M, Oseasohn RO. Cholera Vaccine Field Trials in East Pakistan. 2. Effectiveness in the Field. Bull World Health Organization. 1968;38(3):359–372.

31. Kaper JB, Morris JG Jr, Levine MM. Cholera. Clin Microbiol Rev. 1995 Jan;8(1):48–86. pmid:7704895

32. Levin MM, Kaper JB, Herrington D, Losonsky G, Morris JG, Clements ML, et al. Volunteer studies of deletion mutants of Vibrio cholerae O1 prepared by recombinant techniques Infect Immun. 1988 Jan; 56(1): 161–167.

33. Woodward WE, Mosley WH, McCormack WM. The spectrum of cholera in rural East Pakistan. I. Correlation of bacteriologic and serologic results. J Infect Dis. 1970 May;121:Suppl 121:10+.

34. Thiem VD, Deen JL, Von L, Canh DG, Anh DD, Park JK, et al. Long-term effectiveness against cholera of oral killed whole-cell vaccine produced in Vietnam. Vaccine. 2006 May 15;24(20):4297–303. pmid:16580760

35. Sur D, Lopez AL, Kanungo S, Paisley A, Manna B, Ali M, et al. Efficacy and safety of a modified killed-whole-cell oral cholera vaccine in India: An interim analysis of a cluster-randomized, double-blind, placebo-controlled trial. PLoS Neglected Tropical Diseases. 2011 Oct;5(10):e1289.

36. Hartly DM, Morris JG, Smith DL. Cholera and Post Earthquake Response in Haiti. World Health Organization: Pan American Health Organization. 2011;1(18):1–9.

37. Single-dose azithromycin for childhood cholera. Indian Pediatrics. 2010 Apr; 47(4):305–306. pmid:20431157

38. US Army Corps of Engineers Mobile District and Topographic Engineering Center Water Resources Assessment of Haiti. US Army. 1999 Aug.

39. Almagro-Moreno S, Taylor R. Cholera: environmental reservoirs and impact on disease transmission. Microbiol Spectr. 2013 Dec;1(2):0003–2012.

40. Blower S, Dowlatabadi H. Sensitivity and Uncertainty Analysis of Complex Models of Disease Transmission: An HIV Model, as an Example. International Statistical Review. 1994 Aug;62(2):229–243.

41. Holland J. Notes on R0. Unplublished. Stanford University: Department of Anthropological Sciences. May 2007.

42. Ayeni AO. Domestic Water Source, Sanitation and High Risk of Bacteriological Diseases in the Urban Slum: Case of Cholera in Makoko, Lagos, Nigeria. Journal of Environment Pollution and Human Health. 2014 Jan;2(1):12–15.

43. Diekmann O, Heesterbeek J, Roberts M. The construction of next-generation matrices for compartmental epidemic models. Journal of Royal Society Interface. 2010 Jun 6;7(47):873–85.

44. Alam A, LaRocque R, Harris J, Vanderspurt C, Ryan E, Qadri F, et al. Hyperinfectivity of Human-Passaged Vibrio cholerae Can Be Modeled by Growth in the Infant Mouse. Infection and Immunity. 2005 Oct;73(10):6674–6679. pmid:16177344

45. Merrell DS, Butler SM, Qadri F, Dolganov NA, Alam A, Cohen MB, et al. Host-induced epidemic spread of the cholera bacterium. Nature. 2002 Jun;417(1):642–645. pmid:12050664