Skip to main content
  • Research Article
  • Open access
  • Published:

Comparing alternative cholera vaccination strategies in Maela refugee camp: using a transmission model in public health practice

Abstract

Background

Cholera is a major public health concern in displaced-person camps, which often contend with overcrowding and scarcity of resources. Maela, the largest and longest-standing refugee camp in Thailand, located along the Thai-Burmese border, experienced four cholera outbreaks between 2005 and 2010. In 2013, a cholera vaccine campaign was implemented in the camp. To assist in the evaluation of the campaign and planning for subsequent campaigns, we developed a mathematical model of cholera in Maela.

Methods

We formulated a Susceptible-Infectious-Water-Recovered-based transmission model and estimated parameters using incidence data from 2010. We next evaluated the reduction in cases conferred by several immunization strategies, varying timing, effectiveness, and resources (i.e., vaccine availability). After the vaccine campaign, we generated case forecasts for the next year, to inform on-the-ground decision-making regarding whether a booster campaign was needed.

Results

We found that preexposure vaccination can substantially reduce the risk of cholera even when <50% of the population is given the full two-dose series. Additionally, the preferred number of doses per person should be considered in the context of one vs. two dose effectiveness and vaccine availability. For reactive vaccination, a trade-off between timing and effectiveness was revealed, indicating that it may be beneficial to give one dose to more people rather than two doses to fewer people, given that a two-dose schedule would incur a delay in administration of the second dose. Forecasting using realistic coverage levels predicted that there was no need for a booster campaign in 2014 (consistent with our predictions, there was not a cholera epidemic in 2014).

Conclusions

Our analyses suggest that vaccination in conjunction with ongoing water sanitation and hygiene efforts provides an effective strategy for controlling cholera outbreaks in refugee camps. Effective preexposure vaccination depends on timing and effectiveness. If a camp is facing an outbreak, delayed distribution of vaccines can substantially alter the effectiveness of reactive vaccination, suggesting that quick distribution of vaccines may be more important than ensuring every individual receives both vaccine doses. Overall, this analysis illustrates how mathematical models can be applied in public health practice, to assist in evaluating alternative intervention strategies and inform decision-making.

Peer Review reports

Background

Global conflict, economic plight, and natural disasters interact to displace people on a large scale [1]. Over recent years, an unprecedented increase in refugee and internally displaced populations has led to the largest number of displaced-persons ever on record [2]. Refugee camps are often built in neighboring countries to provide temporary protection and relief for refugees. However, often these crises continue for years or decades and therefore require refugees to live for an indeterminate amount of time in temporary conditions [3]. Because these refugee camps were established to provide temporary rather than permanent shelter, investments in infrastructure have not been prioritized, leading to overcrowding and poor sanitation, in turn increasing the risk of cholera and other infectious diseases [4, 5].

Cholera, a waterborne intestinal infection, causes watery diarrhea and is transmitted through fecal contamination of water and food, as well as person-to-person contact [6]. While the case fatality rate of cholera is low when treated, it can be up to 50% when untreated [7]. Approaches to reduce cholera spread include improvements in water, sanitation and hygiene (WaSH). However, political and economic hurdles can make these longer-term, larger-scale improvements challenging. Vaccination can be used as an complementary approach that results in a substantial reduction in cholera transmission [8, 9]. Further, cholera vaccines can induce herd protection, thereby reducing the risk of disease for both the vaccinated and unvaccinated segments of the population [9, 10]. Although, the effects of vaccination may attenuate over time, recent data has shown that protection remains steady several years after vaccine is administered [11]. Generally, vaccination campaigns can be implemented in conjunction with other health interventions or while more permanent preventative measures (e.g. WaSH) are being put in place.

Maela refugee camp is in northwest Thailand, 8 km east of Burma (Myanmar). In December 2009, its population was 40,009 individuals, who were mostly Burmese refugees. The distribution of ethnic groups across all Thai/Burmese border refugee camps was quite diverse, with ∼61% being of Karen origin [12]. Although cholera has been reported in Thailand in past years [13], a review of the literature between 1982 and 2007 found only 860 cases reported in Thailand (population of over 60 million), with the majority occurring in the northern part of the country [14]. While we do not have a denominator to determine an attack rate for these data, it is clear that compared to the outbreaks reported in the literature, the burden of cholera in Maela had been higher, with more than 1000 cases between 2005 and 2010 among a population of less than 50,000 [15]. In Maela, the majority of individuals have access water and sanitation facilities; however, sociopolitical issues and the mountainous terrain prevent its maintenance or improvement. Therefore, a vaccination campaign was implemented as a critical addition to existing efforts to reduce cholera transmission in Maela. This campaign was the first use of the oral cholera vaccine (OCV) Shanchol in a stable refugee camp [16].

An oral cholera vaccine, Shanchol, was prequalified by the World Health Organization (WHO) in 2011 [17]. The vaccine is administered in two doses, 14 days apart. Efficacy estimates from randomized control trials for the full two-dose series vary by setting and age. A large-scale age-adjusted trial conducted in India found a two-dose efficacy of 65% after follow-up at 5 years [11]. Shorter-term estimates from observational studies have found vaccine effectiveness to be as high as 86.6% after a follow up of 6 months in Guinea [18]. Notably, effectiveness of a reactive vaccination campaign in Haiti was quite close to the efficacy estimate at 63% among individuals self-reporting vaccination [19]. One-dose efficacy was estimated in a randomized trial in Bangladesh and was found to be 52% after 2 years follow up [20]. Observational studies have estimated one-dose vaccine effectiveness to be between 32.5% [21] and 87.3% [9], though the lower bound was not found to be significant. Despite logistical challenges, cholera vaccination in refugee camps has been found to be feasible and acceptable [16, 22, 23] and has been recommended as a potential key intervention to prevent and control cholera transmission in these settings.

In recent years, mathematical modeling has emerged as a useful tool in examining counterfactuals, evaluating intervention strategies, and assisting in policy decision-making [24]. Modeling of cholera transmission, in particular, has been used to guide policy and planning decisions. For example, the US Centers for Disease Control and Prevention (CDC) used real-time modeling to predict the effects of vaccination during the 2010 cholera epidemic in Haiti and to anticipate the total numbers of cases and hospitalizations [25]. More broadly, a wide range of cholera transmission models have been developed [26–31], accounting for different mechanisms, including spatial dynamics [29, 32, 33], age structure [34, 35], environmental drivers [36–39], and disease transmission characteristics such as proportion of asymptomatic individuals [38], hyperinfectiousness [40], dose response effects [27, 29, 41, 42], and multiple transmission pathways [28, 30]. One particularly relevant modeling analysis examined the impact of one compared with two doses of cholera vaccine in Haiti, Zimbabwe, and Guinea when supplies are limited [43]. Another used an agent-based model to examine cholera transmission in a refugee camp setting [44], and one mathematical modeling analysis included displaced-persons to examine cholera dynamics and rainfall in Haiti using a model accounting for environmental transmission [39]. Finally, another used compartmental models to examine cholera transmission and vaccination–induced herd immunity among displaced persons in South Sudan [45]. These models have been useful in explaining different drivers of transmission across settings and evaluating proposed interventions. Modeling of cholera vaccination in displaced person camp settings is growing area of research in need of more attention. A wide range of characteristics such as overcrowding, high levels of population mixing, limited access to clean water and medical care, and political barriers to different interventions make it necessary to develop specific models fit to real-world data to assess the effects of vaccination using population and outbreak data directly from displaced-person settings. These models can provide a useful additional evidence-base to help policymakers make decisions.

In this study, we used the Susceptible-Infectious-Water-Recovered (SIWR) modeling framework [28], expanded to include two-stage vaccination and an adult/child age structure. The SIWR modeling framework accounts for both indirect transmission through environmental water sources as well as a direct pathway representing transmission through food, household water sources, and person-to-person contact. This model is an extension of the classical SIR model (with Susceptible, Infectious, and Recovered compartments to track the total number of individuals at different stages of the disease), with an additional water or environmental compartment representing the concentration of pathogen. The SIWR model has been applied to a range of cholera outbreaks as well as several theoretical studies [28, 29, 32, 42, 46]. It has been integrated with a gravity model to determine how distance and population sizes affect the spread of cholera in Haiti [32] and it was used to estimate the basic reproduction number (\(\mathcal {R}_{0}\)) in a range of settings [29, 30, 32]. Finally, age-structured SIWR models [34, 47, 48] have been used to represent different transmission rates by demographic group.

The primary goal for this study was to develop a transmission model to evaluate a recent vaccine campaign in Maela refugee camp, and use this model to plan for future vaccine campaigns. The Maela refugee camp study setting is particularly relevant given the recent worldwide increase in total numbers of refugees—the approach considered here can be generalized to inform vaccine campaign planning in a wide range of contexts. Because our modeling analysis coincided with preparation and implementation of the cholera vaccination campaign in Maela, we had the opportunity to build and expand our model iteratively based on the vaccination campaign results and implementation. Our model was thus used real-time to help predict the outcome of the campaign, to forecast whether any cholera outbreaks would occur in the near future, and to evaluate the necessity of administering a booster campaign for the years following the campaign. This study illustrates how mathematical modeling can be used iteratively in intervention planning to inform policy and intervention decision-making.

Methods

Study Cohort and Data Collection

Maela refugee camp is situated in the northwest of Thailand, 8 km east of Burma (Fig. 1). In December 2010, the population was 43,645 individuals, who were mostly Burmese refugees. In total, approximately one-third of the residents were children and two-thirds were adults. As per local nongovernmental organization data collected by Première Urgence Aide Medicale Internationale (PU-AMI) in the camp, residents who were 15 years of age or older were considered adults because that was the working age in the camp [15]. Four cholera outbreaks occurred in Maela between 2005 and 2010. The second-largest and most recent outbreak occurred in 2010, when 362 cholera cases were confirmed. Maela has a long-established medical and public health infrastructure (Maela camp was established in 1984 and cholera surveillance has been ongoing since 2005 [15]). Under normal operations (i.e., non-outbreak conditions), laboratory testing for cholera was performed for three groups: inpatients (at the camp hospital) meeting the criteria for suspected cholera (acute watery diarrhea and moderate or severe dehydration); any other patient for whom a clinician strongly suspected cholera based on his or her clinical judgment (such as death associated with acute watery diarrhea especially in a patient ≥5 years of age); and, per week, 5–10 randomly selected outpatients (at either of the two camp outpatient clinics) ≥5 years of age with suspected cholera whenever the number of acute watery diarrhea cases in the preceding week surpassed the 5-year median for that same week. During the 2010 outbreak (after the first test from one of the above came back positive), all patients (inpatient or outpatient) seeking care for suspected cholera received confirmatory testing. Confirmatory testing was performed by isolation of toxigenic Vibrio cholerae O1, and laboratory-confirmed cholera positive cases were defined as ’incident cholera cases’ (also termed confirmed cases). Within 2 days of initial case identification, mobile teams visited the households and neighbors of the case to conduct active case finding and early referral of suspected cases. Further details on the surveillance and laboratory procedures can be found in [15].

Fig. 1
figure 1

Maela Map and Model Schematic. Left: Map showing the location of Maela refugee camp. Map figure was generated in R version 3.2.4 [49] using the ‘raster’ package [50] which uses data from [51] photograph taken by JH. Right: Flow diagram of Age-structured SIWR-based model of cholera transmission in Maela. Vital dynamics and migration rates are represented by dotted lines. All infected classes (second row) shed pathogen into the common water source (W), which can subsequently infect susceptible individuals depending on their vaccination status

Following WHO prequalification of Shanchol cholera vaccine in September 2011, the Thailand Ministry of Public Health sponsored a campaign for the population of Maela. The campaign was implemented in 2013 by PU-AMI with technical support from CDC. Pregnant women and children <1 year old were not given vaccine (approximately 2000 refugees were excluded in total) per the manufacturer’s recommendations [52]. Overall approximately 81% of refugees were given at least one dose of vaccine, while 64% of refugees were given two doses [15]. A baseline census was conducted and then cholera vaccine was subsequently administered. Residents of Maela have since been and will continue to be followed prospectively for new suspected cholera cases, as well as regular laboratory testing of a subset of watery diarrhea cases in the camp [15].

Our analysis considers the total population at the time of the campaign: 45,233 individuals, of whom 27,901 are adults and 17,332 are children. We fit our model to incidence data from the 2010 outbreak, the most recent outbreak at the time of the campaign, as shown in Fig. 2.

Fig. 2
figure 2

Model fit to data. Model fit to Maela cholera data from the 2010 epidemic for adults (left) and children (right). Data points for confirmed cases are indicated as circles, and simulation output is indicated as lines

Model Structure and Parameter Estimation

Age-structured SIWR Model

The model structure is shown in Fig. 1. The population of Maela was separated into six model classes, non-vaccinated-adults (variables marked with a), once-vaccinated-adults (variables marked with Va), twice-vaccinated-adults (variables marked with VVa), non-vaccinated-children (variables marked with c), once-vaccinated children (variables marked with Vc) and twice-vaccinated-children (variables marked with VVc). Each class is further broken into Susceptible-Infectious-Recovered compartments. Any infectious individual can infect any susceptible individual regardless of model class. The force of infection, λ, is the rate at which susceptible individuals become infectious (incidence rate). Children were defined as individuals <15 years old. Environmental water source contamination is tracked in a separate compartment (W), into which any infected individual can shed. Additionally, the pathogen concentration in the water, (W), contributes to the calculation of λ.

The equations for the full model are given in the Supporting Information. Below are the force of infection equations for each demographic class and the environmental pathogen equations which are not explicitly represented in model schematic are shown in Fig. 1.

Force of Infection Equations

$$\begin{array}{@{}rcl@{}} \begin{aligned} \lambda_{a} &=\beta_{aa}I_{a} + \beta_{ca}I_{c} + \beta_{wa}W + \beta_{aa}IV_{a} + \beta_{ca}IV_{c}\\&\quad + \beta_{aa}IVV_{a} + \beta_{ca}IVV_{c}\\ \lambda_{va} &= (1-V\!E_{1}) (\beta_{aa} I_{a} + \beta_{ca} I_{c} + \beta_{wa} W + \beta_{aa} IV_{a}\\&\quad + \beta_{ca} IV_{c} + \beta_{aa} IVV_{a} + \beta_{ca} IVV_{c})\\ \lambda_{vva} &= (1-V\!E_{2}) (\beta_{aa} I_{a} + \beta_{ca} I_{c} + \beta_{wa} W + \beta_{aa} IV_{a}\\&\quad + \beta_{ca} IV_{c} + \beta_{aa} IVV_{a} + \beta_{ca} IVV_{c})\\ \lambda_{c} &= \beta_{ac}I_{a} + \beta_{cc}I_{c} + \beta_{wc}W + \beta_{ac}IV_{a} + \beta_{cc}IV_{c}\\&\quad + \beta_{ac}IVV_{a} + \beta_{cc}IVV_{c}\\ \lambda_{vc} &= (1-V\!E_{1}) (\beta_{ac} I_{a} + \beta_{cc} I_{c} + \beta_{wc} W + \beta_{ac} IV_{a}\\&\quad + \beta_{cc} IV_{c} + \beta_{ac} IVV_{a} + \beta_{cc} IVV_{c})\\ \lambda_{vvc} &= (1-V\!E_{2})(\beta_{ac} I_{a} + \beta_{cc} I_{c} + \beta_{wc} W + \beta_{ac} IV_{a}\\&\quad + \beta_{cc} IV_{c} + \beta_{ac} IVV_{a} + \beta_{cc} IVV_{c}) \\ \end{aligned} \end{array} $$
(1)

Environmental Pathogen

$$\begin{array}{@{}rcl@{}} \begin{aligned} \lambda_{w} &= I_{a} + \sigma I_{c} + IV_{a} + \sigma IV_{c} + IVV_{a} + \sigma IVV_{c}\\ \dot{W} &= \alpha \lambda_{w} - \xi W \end{aligned} \end{array} $$
(2)

Within the force of infection equations, the βij’s represent transmission from i to j. Infectious individuals shed into W based on σ, the child to adult shedding ratio of pathogen, the number of infected individuals in each class, and α, the rate at individuals shed into the environment. Furthermore, ξ is the rate at which bacteria decay in the environment. Within the remaining equations (shown in the Supporting Information and Fig. 1), the ν’s represent the rates of vaccination for susceptible and recovered individuals. Therefore, individuals can move from nonvaccinated to once-vaccinated and then onto twice-vaccinated depending on which vaccine scenario is being simulated (see below for details on vaccination scenarios). The μ’s represent the combined rates for out-migration and death, and the M’s are the in-migration rate of susceptible or recovered individuals in Maela. B is the birth rate of individuals into the population. Individuals are born as susceptible children. The γ’s are the recovery rates of infectious individuals. All parameter definitions, values and sources are given in Table 1.

Table 1 Model parameter values and estimates. Estimated parameters are indicated in bold

Measurement Model

The disease surveillance data from the 2010 epidemic measure weekly cholera incidence among adults (≥15 years old) and children (<15 years old). To accurately represent this in the model, simulated weekly case counts are scaled by estimated reporting fraction parameters, ka and kc. These scaling parameters represent a combination of factors including the cholera reporting rate (to account for missed cases, e.g. cases with milder symptoms where the individual did not seek medical attention), the fraction of asymptomatic cases (as these would not be reported), and a correction for any errors in the population size (i.e., if the population at risk is larger or smaller than the total recorded population of Maela).

Parameter Estimation from the 2010 Epidemic

Conducting an identifiability analysis is a necessary step before estimating parameter values from a fitted model. Identifiability is typically broken into two broad categories. First, structural identifiability is used to examine how the structure of the model and measured variables can affect what parameters are capable of being estimated assuming perfect, noise-free data [30, 55–59]. Once this is established, practical identifiability is used to examine the identifiability of parameter values given the actual data set being used. This step accounts for real-world data issues such as noise and sampling frequency [60, 61].We examined both structural and practical identifiability before parameter estimation. For more information on the techniques and subsequent assumptions made to reconcile identifiability issues, see Supporting Information.

We initially fit the model to estimate unknown parameters using 2010 cholera incidence data from Maela. Because the vaccine campaign was not implemented until 2013, the vaccinated compartments and relevant parameters in the full model were not included in the parameter estimation. Excluding vaccination reduced the model substantially to just three compartments for adults, three compartments for children, and the environmental pathogen compartment (W).

Similar to the original SIWR model [30], the structural identifiability analysis of our simplified model indicated that waterborne transmission parameters and α were not separately identifiable for our model, and instead formed an identifiable combination. We therefore re-scaled W and βW to enable all unknown model parameters (βI,βW,σ,ξ, and the k’s) to be structurally identifiable (see Supporting Information for more details). From this point forward, we have only used the re-scaled versions of βW and W. We next estimated the structurally identifiable set of unknown parameters from the data, using Poisson maximum likelihood with Nelder-Mead optimization (using fminsearch in Matlab).

Prior to fitting, we also fixed the demographic parameters μa,μc,Ma,Mc and B based on PU-AMI data collected in the camp. The infectious period, γ, was fixed to 4 days [53, 54]. Further, because of apparent practical identifiability limitations discovered during the initial parameter estimation (see Supporting Information for more details), we set all human-human transmission parameters denoted βI, equal to each other, and separately we set all human-water transmission parameters, denoted βW, equal to each other. We subsequently set σ, the child-to-adult shedding ratio of pathogen into the environment, to 1 since the βW’s were equal (similar to [42]). Thus the remaining parameters (ka,kc,βI,βW,ξ) were assumed unknown.

Basic Reproduction Number (\(\mathcal {R}_{0}\))

The basic reproduction number, \(\mathcal {R}_{0}\), is the number of secondary infections caused by a introduction of infectious material (individuals or pathogens) in a completely susceptible population [62] and is a commonly used measure of disease transmissibility [63]. \(\mathcal {R}_{0}\) for the model was determined using the next-generation approach [64]. For the simplified and scaled model without vaccination, \(\mathcal {R}_{0}\) is given by

$$\begin{array}{@{}rcl@{}} \begin{aligned} \mathcal{R}_0 &= \frac{N_{A} (\beta_{I} + \beta_{W})}{\gamma+\mu_{A}} + \frac{N_{C}(\beta_{I} + \sigma \beta_{W})}{\gamma + \mu_{C}}\\ \end{aligned} \end{array} $$
(3)

where NA and NC are the total adult and child populations in the camp. As we are assuming that the simulation starts near the disease free equilibrium, the total population is equivalent to the completely susceptible population (apart from the initial infected cases which begin the epidemic). Equation (3) shows that the overall \(\mathcal {R}_{0}\) is given by a sum of the contributions by adults and children, with waterborne transmission among children weighted by the ratio of child-to-adult shedding rates (σ). Note that the individual terms for adult and child contributions to \(\mathcal {R}_{0}\) follow the same general form as for the original SIWR model [28].

Vaccination Strategies: Exploration of Dynamics

To explore the effects of Shanchol on cholera transmission dynamics in this setting, we conducted hypothetical vaccination scenarios varying timing and dosage while assuming a limited amount of vaccine was administered. We incorporated different vaccination effectiveness estimates (VEs) into the transmission parameters of the vaccinated groups, and used the best-fit parameter estimates as our baseline parameters for all scenarios; see Table 1 for all parameter values. We assume that vaccination reduces the susceptibility of individuals, but that if infected, individuals are equally infectious regardless of past vaccination or cholera exposure history. We set the vaccine effectiveness in our model to be the lower bound of all estimates from recent Shanchol studies to ensure that our results do not over-estimate the impact of Shanchol. Specifically, we assumed that two doses confer an effectiveness of 63% (corresponding to the estimate from Haiti [19]) and one dose confers an effectiveness of 32.5% (corresponding to the estimate from India [21]). In all scenarios, we start with one initial observed infected individual in both Ia and Ic (i.e. \(\frac {I_{a}}{K_{a}}\) and \(\frac {I_{c}}{K_{a}}\), respectively) to establish consistency between scenarios and because reactive vaccination scenarios could occur only after cholera has been observed. We examine alternative pre-vaccination seeding scenarios in Supporting Information Additional file 1: Table A2 and in the forecasting section below. All scenarios are simulated for 365 days and do not include waning immunity. Since there was only 1 death from cholera during the 2010 outbreak in Maela, we set infectious and non-infectious mortality rates equal to each other [15]. In the baseline scenario, we simulate a cholera outbreak with no vaccination campaign.

We next examined all possible one and two-dose pre-vaccination combinations (wherein all doses of vaccine are administered before the outbreak begins) by varying the proportion of the population that was given one or two-doses. As a test case for subsequent analyses, we assumed that only 20,000 doses of Shanchol were administered during the outbreak. We chose 20,000 because it was a rounded approximation of 50% of the population of Maela, which allows us to evaluate the disease dynamics in the presence of potential logistical or resource limitations. We considered two approaches: pre-vaccination and reactive vaccination (in which vaccine is administered after cholera is detected). For both approaches, we considered four general strategies:

  • Two-dose: two doses of vaccine administered to 10,000 individuals. For reactive vaccination, this is implemented as two 4-day campaigns, the first of which begins one week after the first detection of cholera and the second which begins two weeks later, in accordance with Shanchol administration guidance [6].

  • One-dose: single dose of vaccine administered to 20,000 individuals. For reactive vaccination, this is implemented as an 8-day campaign which begins one week after the first detection of cholera.

  • Mixed: single dose of vaccine is administered to 10,000 individuals, and two doses of vaccine are administered to another 5000 individuals.

  • First come, first served: This strategy is structured similarly to the two-dose strategy, except that we do not track or control the proportion who receive one or two-doses—the 20,000 doses are administered on a first come, first served basis in two 4-day campaigns (with no prespecified number of single doses or two doses). This results in ∼14,000 individuals receiving one dose and ∼3000 individuals receiving two doses.

Unless otherwise indicated, we distributed vaccine doses proportionally among adults and children. The durations of the reactive vaccine campaigns were determined based on estimates of the number of doses that could be administered per day in Maela informed by the 2013 campaign. In all reactive vaccination strategies, we assume that it takes 1 week for antibodies to confer protection. One week was chosen because a significant increase in titer levels occurred 7 days post immunization, and no significant increase in titers levels occurred beyond that time point in Shanchol seroconversion studies [65, 66]. Therefore, vaccine administered on day 7 of a cholera outbreak will not confer protection until day 14.

Given the uncertainty and variability in effectiveness estimates (they depend greatly on study design features like length of follow up), we conducted two additional analyses. First, we simultaneously varied one-dose effectiveness from 0% to 63% and number of total doses administered from 0 to twice the total camp population for all pre-vaccination scenarios. The maximum possible one-dose effectiveness was assumed to correspond to the two-dose effectiveness estimate. Second, we examined different reactive vaccination delays from 0 to 8 weeks and simulataneously varied the one-dose effectiveness from 0% to 63%.

Finally, because VE estimates in children are lower [67] but our model uses pooled VE estimates across all ages, as a sensitivity analysis we also tested a lower VE among children, by dividing the one- and two-dose VE among children by 2 to examine whether a lower VE among children would affect our overall results [67].

Model simulation and parameter estimation used the ode23tb solver and the fminsearch optimization function in MATLAB [68].

Maela Campaign Forecasting

We next used the actual campaign coverage information for Maela to forecast the potential effects of the vaccine campaign should cholera be introduced. We forecasted the total case counts that would occur in the event of an introduction of cholera in two scenarios: (1) early summer 2013 immediately after the OCV campaign had ended, and (2) 2014 assuming no cholera introduction in 2013. We simulated the 2013 scenario using OCV campaign coverage and for the counterfactual scenario of no 2013 OCV campaign. The 2014 scenario was used to assess the need for a booster campaign.

As we were simulating over a longer time than in the vaccination scenario runs, we updated our simplified model to include waning immunity. For this, we supposed that the highest level of protection (Ra and Rc) comes from having recovered from a natural infection. We then supposed this immunity can wane into two compartments: V2a and V2c, equivalent to receiving two doses of vaccine, which then wane into V1a and V1c, equivalent to receiving only one dose of vaccine. The V1a and V1c compartments then wane into the susceptible classes, Sa and Sc. In this model, vaccinated are considered in the appropriate V2 or V1 compartment, rather than the SVV and SV compartments, and we assume that if a vaccinated individual is infected, they simply join the overall Ia and Ic classes, so that once recovered they are fully immune (Ra and Rc). For the resulting model equations, see the Supporting Information.

For initial conditions, we used the vaccination coverage data from the actual campaign in the total camp population (i.e., included and excluded individuals) to determine initial vaccination status, wherein 22.3% of adults were vaccinated with only one dose, 49.3% of adults were vaccinated with two doses, 18.7% of children were vaccinated with only one dose, and 63.6% of children were vaccinated with two doses. We note that we adjusted the percent coverage of individuals included in the campaign from that reported in [15] to account for the total camp coverage of both included and excluded individuals. See the Supporting Information for details on these calculations. In addition to the vaccination campaign, we also considered the possibility of pre-existing immunity (because there had been cholera epidemics previously in the camp). We therefore ran all scenarios with two options: assuming a fully susceptible population, and alternatively a partially immune population (up to 50%). We denote the fraction of the initial population which is immune to be Fimm (split evenly across the immune compartments). We also allowed the fraction of immune individuals (assumed to all be at V1-level immunity) in incoming migrations, denoted Mimm, to range from 0 to 1. Finally, to model the introduction of cholera, we added a single adult case and a single child case into the model. For a cholera introduction in early summer 2013, each of the following scenarios was run:

  • No OCV campaign in a fully susceptible population

  • OCV campaign in a fully susceptible population

  • No OCV campaign in a partially immune population

  • OCV campaign in a partially immune population

For a cholera introduction in 2014, we ran the model without cholera introduction for 2013 (to allow the waning immunity, migration, birth, and death dynamics to continue), and then for 2014 each of the following scenarios was run twice (once seeding with observed cases and once seeding with actual cases):

  • 2013 OCV campaign in a fully susceptible population

  • 2013 OCV campaign in a partially immune population

To incorporate uncertainty in our parameter values into our projections, we used Latin hyper-cube sampling (LHS). We sampled 1000 parameter sets for each forecasting scenario using the parameter ranges shown in Table 2. For most parameters, we based our ranges on the maximum and minimum values observed from the Maela demographic data from 2009 to 2013, as well as parameter ranges evaluated from the literature (see Table 2). Although the ranges of one compared with two-dose effectiveness overlap in studies, we assume that the maximum one-dose effectiveness is equal to the minimum two-dose effectiveness. For the remaining transmission and reporting rate parameters, we used broad ±50% ranges. The ranges used for the estimated parameters were generally similar to or wider than their estimated confidence bounds but may better reflect the additional uncertainties coming from the fact that each introduction of cholera may be different (in time of year, contact patterns among the population, reporting, ongoing interventions, etc.).

Table 2 LHS ranges among fitted parameters for forecasting scenarios

Results

Parameter Estimation and Uncertainty

Figure 2 shows the model fitted to weekly incidence data for adults and children, with parameter estimates given in Table 1. To formally examine parameter uncertainty and practical identifiability, we plotted profile likelihoods for fitted parameters ka,kc,βW,βI,ξ, shown in Supporting Information Additional file 1: A1. All parameters were shown to be identifiable (finite confidence bounds), with clear minima in each case, although the uncertainty was comparatively high for βW and ξ, consistent with previous studies showing that these two parameters are often practically unidentifiable [30, 42].

Vaccination Strategies: Exploration of Dynamics

We assessed different roll-out vaccination campaign strategies to determine the most effective method of preventing a cholera outbreak in Maela.

Baseline Scenario – No Vaccination

The baseline scenario in which no vaccination is given yields a total of 395.4 cases with an attack rate of 8.7 per 1000 people. See Table 3 for case counts and attack rates from the theoretical exploration of dynamics analysis.

Table 3 Attack rates and numbers of cases for each theoretical scenario (assuming 20,000 doses) using best fit parameter values

Pre-Vaccination Scenarios

We first examined how variation in pre-vaccination coverage of one compared with two doses affects cumulative case counts. As one-dose or two-dose coverages increase, case counts decrease. The cutoff for one-dose coverage alone to result in less than 50 total cases is 96% and the corresponding cutoff for two-dose coverage alone is 49%. All the pre-vaccination scenarios we considered administering 20,000 doses result in similar cumulative case counts, while the real-world Maela scenario (with higher coverage) was more effective. See Fig. 3 for details.

Fig. 3
figure 3

Cumulative cases varying pre-vaccination coverage. Cumulative cases resulting from methodically varying one-dose vs. two-dose coverage for pre-vaccination. With labels indicating pre-vaccination scenarios with 20,000 doses. (FCFS = first-come, first-served)

We then simulated pre-vaccination scenarios administering 20,000 doses to Maela and using the one-dose, two-dose, mixed, and first come, first served strategies. All four pre-vaccination campaigns yielded similar results. An outbreak occurs, but it is substantially smaller than the baseline scenario. Total case counts range from 247.6 with an attack rate of 5.6 per 1000 people in the two-dose scenario to 255.3 with an attack rate of 5.6 per 1000 people in the one-dose scenario. See Fig. 4 and Table 3 for details.

Fig. 4
figure 4

Cumulative cases for alternative pre-vaccination scenarios after administering 20,000 doses. Cumulative cholera cases in adults and children, for alternative pre-vaccination scenarios. Where ’nv’ is the baseline, no vaccination scenario; ’p1d’ the one-dose scenario; ’p2d’ is the two-dose scenario; ’pmixed’ is the mixed scenario, and pfcfs is the first-come first-served scenario

Finally, we varied one-dose effectiveness (from 0 to the full two-dose effectiveness) and total number of doses (from 0 to twice the total camp population (full coverage with two doses)) in all pre-vaccination scenarios. Intuitively, when one-dose effectiveness is low, the two-dose scenario is most effective. For instance, administering the full two-dose series to 22,000 people is sufficient to achieve <50 cases. However, the one-dose strategy alone can achieve large reductions in case counts if the combination of effectiveness and number of doses is sufficient. For instance, either an effectiveness of 31% combined with 45,000 doses or an effectiveness of 63% combined with 22,000 doses are minimally sufficient to achieve <50 cumulative cases. If one-dose effectiveness is very high and approximately equal to two-dose effectiveness, the two scenarios behave similarly. See Fig. 5 for details.

Fig. 5
figure 5

Cumulative cases varying one-dose vaccine effectiveness and total doses for different pre-vaccination scenarios. Cumulative cases resulting from methodically varying one-dose vaccine effectiveness and total doses administered in campaign implementation for the all pre-vaccination scenarios. (FCFS = first-come, first-served)

Reactive Vaccination Scenarios

Next, we examined reactive vaccination scenarios, again administering 20,000 doses given in separate campaigns starting 7 days after cholera is initially detected. The same one-dose, two-dose, mixed, and first come, first served scenarios were simulated, shown in Fig. 6 and Table 3.

Fig. 6
figure 6

Cumulative cases for alternative reactive vaccination scenarios after administering 20,000 doses. Cumulative cholera cases in adults and children, for alternative reactive vaccination scenarios. Where ’nv’ is the baseline, no vaccination scenario; ’rv1d’ is the one-dose scenario; ’rv2d’ is the two-dose scenario; ’rvMixed’ is the mixed scenario; and ’fcfs’ is the first come, first served scenario

In all vaccination test-case scenarios (shown in Fig. 6), an outbreak occurs but it has substantially lower numbers of cases than the baseline scenario. The one-dose, two-dose, mixed, and first come, first served scenarios all yield similar results, with total case counts ranging from 267.6 with an attack rate of 5.9 per 1000 people in the one-dose scenario to 273.7 with an attack rate of 6.1 per 1000 people in the first come, first served scenario. Given the fact that the one-dose effectiveness is almost exactly half of the two-dose effectiveness, the one-dose scenario provides approximately the same amount of coverage as the other scenarios.

The case counts are higher and the differences between scenarios are lower in the reactive vaccination scenarios than in the pre-vaccination scenarios, indicating that pre-vaccination is more effective at reducing the spread of cholera.

To evaluate the effects of vaccine campaign timing, we methodically varied one-dose vaccine effectiveness from 0% to 63% and initial vaccination campaign onset from 0 days to 55 days (or 8 weeks) after cholera was first identified for all reactive vaccination scenarios. Results are shown in Fig. 7. As with the pre-vaccination scenarios, when one-dose effectiveness is low or uncertain, the two-dose strategy is preferable. For instance, if one-dose effectiveness is 20%, the minimum number of cases (i.e. assuming no delay) for the one-dose strategy is 320 while the minimum for the two-dose strategy is 291.8. Once one-dose effectiveness is sufficiently high, all campaigns are comparable given relatively short delays. This is illustrated by the 20,000 dose reactive vaccination scenarios (see Fig. 6). A one-dose effectiveness of 27% is minimally sufficient for it to be the preferable strategy across all campaigns and given any delay. When reaction time is important, one-dose is often the preferred strategy because all doses of vaccine are administered in a single 8 day campaign instead of across two 4 day campaigns separated by 14 days. This results in more population level protection earlier in the epidemic. Lastly, beyond a certain point in the epidemic, one-dose effectiveness has little impact on total case counts across all campaigns. This interplay between vaccination effectiveness and timing highlights the importance of minimizing delays in reactive vaccination campaigns with pre-vaccination being the most effective way to reduce the spread of cholera (see Fig. 7).

Fig. 7
figure 7

Cumulative cases varying one-dose vaccine effectiveness and vaccine campaign delay for reactive vaccination scenarios. Cumulative cases resulting from methodically varying one-dose vaccine effectiveness and delay in campaign implementation for all reactive vaccination scenarios

Finally, to account for lower VE estimates in children [67], we we examined the relationship between one-dose effectiveness and total number of doses for pre-vaccination and the relationship between one-dose effectiveness and delay in campaign implementation for reactive vaccination, with a reduced VE among children (divided by 2). Overall, the results did not change with reduced VE estimates among children (results shown in Supporting Information Additional file 1)

Maela Campaign Forecasting

Table 4 shows the forecasted numbers of cases and attack rate for each of the 2013 and 2014 scenarios using the parameter estimates in Table 1.

Table 4 Forecasts for the 2013 and 2014 cholera seasons

Forecasts for the 2013 OCV Campaign

In the 2013 forecasting results, we see a larger spread of total case numbers for runs in the scenario without the OCV campaign compared to the scenario with the OCV campaign. Case counts range from 0 to approximately 1000 in the fully susceptible population. The partially immune population runs generally have lower case counts. Furthermore, for the scenarios that consider the OCV campaign, we see the vast majority of runs having case counts close to 0. For forecasts among a fully susceptible population, see Fig. 8, and for forecasts among a partially immune population, see Supporting Information Additional file 1: Figure A7.

Fig. 8
figure 8

2013 forecasting results among a fully susceptible population. 2013 forecast with a single actual case in adults and children as seeding. Simulations with no OCV campaign on left and with OCV campaign on right

Forecasts for the 2014 Cholera Season

The 2014 forecasting results are quite similar to the 2013 runs. The partially immune population results in more simulations with 0 total cases, compared with the fully susceptible population. Because population immunity wanes between 2013 and 2014, there is higher proportion of larger outbreaks for the 2014 forecasting scenarios, but the vast majority of runs remain close to 0 for both populations. For details, see Supporting Information Additional file 1: Figure A8.

Discussion

Using a mathematical model of cholera transmission in a refugee camp, we have shown the dramatic effect that vaccination can have on reducing the risk of cholera outbreaks in refugee settings. Our analyses suggest that pre-vaccination campaigns, even if not fully covering the whole camp population, can lead to substantial reductions in the number of observed cases in the event of an outbreak. Of course, the best pre-vaccination strategy depends on both one-dose effectiveness and the amount of doses available. If resources are limited and one-dose offers acceptable protection, single dose-based strategies would be preferred. But if there are enough vaccine doses available or the effectiveness of a single dose is questionable, then full vaccination strategies are better. We also found that reactive vaccination strategies can result in a moderate reduction in the number of infections. Further, the small differences in attack rates in our model results suggest that, at the best-fit parameters, vaccinating more individuals with at least one dose (even if the full vaccination sequence is not completed) may be nearly or more effective than vaccinating fewer individuals with a complete two-dose program assuming conservative vaccine effectiveness estimates. This consideration is particularly important if vaccination occurs later in the outbreak, where time is of the essence as the outbreak is already ongoing—the full two-dose sequence requires a delay between doses, hindering its effectiveness in reactive vaccination campaigns that begin later in an epidemic. Additionally, if logistical constraints limit the ability to follow up and provide patients with a second dose, a one-dose strategy may be preferable. However, the relative one vs. two-dose effectiveness, timing and number of doses available should all be considered. Finally, our projections of future outbreaks after the 2013 vaccination campaign in the Maela refugee camp suggest that vaccination may have prevented outbreaks in 2013 and 2014, as no cases were observed in either year, both in the model and during follow-up in the camp. This ability to consider counterfactual scenarios and generate projections highlights the potential of modeling to help guide real-time public health decision-making.

Our mathematical model of cholera transmission reproduces the dynamics observed in Maela’s 2010 cholera outbreak (Fig. 2). We conducted a theoretical exploration of the dynamics of our model to examine different vaccination scenarios with only 20,000 doses of Shanchol distributed in the camp to reflect potential public health intervention strategies. Doses were distributed proportionally among children and adults. These results provide insight into the most effective strategies for vaccination when logistics might impede complete coverage of a population with one or two doses. We found that the two-dose pre-vaccination strategy was marginally the most effective with an attack rate of 5.5 cases per 1000 people, but others provided comparable protection. On the other hand, the most effective reactive vaccination scenarios are the one- and two-dose strategies with attack rates of 5.9 cases per 1000 people (Table 3), while again others provided similar protection. All similarities between scenarios are the result of the conservative effectiveness estimates we used such that one-dose effectiveness is nearly exactly half of two-dose effectiveness.

Cholera vaccines provide protection for individuals who do not directly receive vaccination [9, 10]. As with many transmission models, our analysis implicitly accounted for indirect protection (herd immunity). Specifically, the direct vaccination of individuals in the model reduces transmission in the population which in turn provides indirect effects protecting nonvaccinated individuals from becoming infected. To illustrate this, we calculated the number of cases that occur among nonvaccinated individuals in the two-dose pre-vaccination strategy with the VE used in our analysis and compared this to the number of cases that occur among nonvaccinated individuals in the same scenario when the VE was set to 0%. We found that when we set the VE to 0%, 308.4 cases occurred among nonvaccinated individuals, but when the VE was set to our values from the literature (see Table 1), only 220.9 cases occurred among nonvaccinated individuals. Thus, our model shows a herd protection effect for individuals who did not receive cholera vaccine.

For pre-vaccination, all scenarios can achieve large reductions in case counts. If one-dose effectiveness is low or uncertain, the two-dose scenario is preferable since it can guarantee larger or more certain reductions in case counts i.e., 44,000 doses (administered to 22,000 people) are sufficient to achieve <50 cumulative cases. On the other hand, if one-dose effectiveness is only ∼50% that of two-dose effectiveness, 45,000 doses (administered to 45,000 people) are sufficient (see Fig. 5 for details). In the reactive vaccination scenarios, two doses may be preferable if one-dose effectiveness is low or uncertain. Further, we see a crucial interplay between timing and vaccine effectiveness on transmission. Past a certain point in the outbreak, one-dose effectiveness does not substantially change the total case counts. For instance at day 50 in the one-dose scenario, the largest potential difference in cases is 65. (total cases: 395.3 and 333.2 for a one-dose effectiveness of 0% and 63%, respectively) while at day 20 the largest potential difference in cases is 202.2 (total cases: 395.5 and 202.3 for a one-dose effectiveness of 0% and 63%, respectively). Additionally, when one-dose effectiveness is sufficiently high, the one-dose strategy is preferable, with one-dose able to achieve the lowest case counts during an idealized situation in which there is no delay in reactive vaccination administration and the one-dose effectiveness is nearly equal to the two-dose effectiveness. Although this scenario may not be possible in real-world settings, it underscores the necessity of considering delays and relative one vs. two-dose effectiveness in reactive vaccination campaigns. For the one-dose scenario considered here, more individuals are given vaccine over an 8 day campaign resulting in more people having some protection early in the outbreak. On the other hand, in the two-dose scenario individuals who have been vaccinated have a higher level of protection, but achieving the same amount of population level protection will take a longer period of time i.e., two 4 day campaigns occurring 14 days apart (see Fig. 7 for details). The results of the forecasting analysis for 2013 show that regardless of population level immunity, the OCV campaign using a partially mixed strategy (coverage levels shown in Additional file 1: Table A4) prevents a majority of outbreaks that might otherwise have occurred, shown in Fig. 8 and Additional file 1: Figure A8. Additionally, even if we assume a fully susceptible population prior to the 2013 OCV campaign, the vast majority of post-OCV campaign runs in 2014 still result in no outbreak. Indeed, the median of total cases across all runs is <1 (see Table 4 and Additional file 1: Figure A8). This suggests that an introduction of cholera would likely not have resulted in a significant outbreak even given a conservative assumption about the population level of immunity. Thus, it was determined that there was no need for a booster campaign in Maela (and indeed there was no cholera outbreak that year).

A key limitation of our model is the uncertainty in parameter values. Because we fit our model to Maela data from 2010, certain parameter values (e.g., birth and death rates) might not be completely representative of Maela in 2013 and 2014. However, to assess the uncertainty in our parameter values, we conducted global sensitivity analyses using LHS and assessed both practical and structural identifiability of the model. The quantitative results are heavily dependant on the vaccine effectiveness estimates however, we chose the lower bounds of values from recent studies to ensure that the model generated conservative results. Furthermore, the qualitative conclusions i.e. the interplay between vaccine timing and vaccine effectiveness will occur for most realistic effectiveness estimates. We also examined the effects of reducing child VE estimates and found that this did not substantially change any of our conclusions. Other weaknesses of this analysis are inherent in the model assumptions. For example, one key assumption for the vaccination scenario simulations is that there is no waning immunity because we are simulating over a short time course. We do, however, incorporate waning immunity into the forecasting analysis. Another assumption is that each infectious individual is equally infectious regardless of previous vaccination, cholera exposure, or time since infection. In any case, this assumption would result in an overestimate of the total number of cases because we are ignoring the fact that individuals who have been vaccinated might be less infectious. We are also ignoring hyperinfectiousness, as incorporating this into our model as including it would require tracking of pathogen through the human host and a shorter time scale because it decays after 18 hours [78]. Although we are not explicitly modeling asymptomatic infections or errors in disease reporting, these are accounted for in the scaling factors, ka for adults and kc for children. Additionally, in our analyses we used a deterministic model which will not capture stochastic fluctuations. Stochasticity may play a role, particularly in the early phase of an outbreak following a new introduction of cholera—our model may therefore miss some of the stochastic die-out of epidemics in our forecasts, which again, would lead to an overestimate of the number of cases. We assume that the mortality rate for individuals with cholera is the same as that for individuals without cholera. Although this is likely not true, the data we fit our model to did not have a sufficient number of deaths to calculate case fatality rates for diseased compared with non-diseased population groups (there was only one death among individuals with cholera). As individuals in Maela have easy access to medical care, they are likely to be treated and have a high rate of survival. For the sake of parsimony [79], we are also assuming that human-human transmission parameters are equal across demographic groups and separately human-water transmission parameters are equal since we obtained similar fits for a range of β values in the practical identifiability analysis (see Supporting Information). We are using data from a specific setting, which may limit the external validity of our results; however, the qualitative dynamic results obtained should be consistent for other refugee camp settings and have been seen in similar mathematical modeling analyses [43, 80].

The strengths of this model include the novelty of considering vaccination scenarios while explicitly accounting for environmental transmission of cholera in a refugee camp. Additionally, our identifiability and sensitivity analyses methodically considered parameter uncertainty. Another strength is the fact that we used real-world data to inform our model which in turn, provided insight for the public health response. Specifically, we liaised with the CDC as well as local nongovernmental organizations in real-time and used data from the OCV campaign to directly inform the model. Further, the model was then used to explore counterfactual forecasting scenarios to help answer outstanding questions among trial investigators about whether or not a booster campaign was necessary. Finally, the consistency of our results with other analyses indicates the robustness of our findings.

In general, our results indicate that vaccination should be considered in conjunction with WaSH with the caveat that immunity may wane over time. Furthermore, the trade-off between vaccine effectiveness (i.e., one-dose compared with two-dose) and timing of reactive vaccination should be carefully considered. The WHO currently holds a stockpile of over 3 million doses of oral cholera vaccine to allow countries or institutions to request doses of vaccine during cholera outbreaks. The average time from when requests were approved to receipt in country has been 14.4 days with an additional 9.5 days until vaccination actually started [81, 82]. As seen in our vaccination scenario results, timing is crucial to the impact of vaccination. Recent mathematical modeling work has examined how best to allocate global stockpile reserves [83].

Ideally, pre-vaccination should be considered as a short-term transmission reduction strategy, compared to the potentially longer-lasting effects of improved WaSH. A recent study used a static model fit to data from Malawi to estimate cases averted in Haiti by implementation of oral cholera vaccine and/or WaSH and found that a combined implementation of WaSH and vaccine resulted in the greatest reductions in cases [84]. Furthermore, WaSH reduces transmission for a wide range of infectious diseases and if maintained is more permanent, while Shanchol targets cholera, and its effects do not last as long—although vaccine campaigns may be easier to implement, as they do not require sustained maintenance. Thus, both WaSH and vaccination may have their roles to play in an effective intervention strategy.

Our analyses suggest that vaccination provides an effective strategy for preventing cholera outbreaks in refugee camps and that cholera vaccination should be considered, even in the absence of an ongoing outbreaks. Given the dramatic increases in displaced populations and refugee settlements across the world it is critical that vaccination be considered with water sanitation and hygiene improvements. If a camp is facing an outbreak, delayed distribution of vaccines can substantially alter the effectiveness of a reactive vaccine campaign, suggesting that quick distribution of vaccines (e.g., using a first come, first served approach) may be more important than ensuring that every individual gets both vaccine doses.

Conclusions

We developed an age-structured SIWR-based transmission model to consider different cholera vaccination strategies in Maela, the largest and most long-standing refugee camp in Thailand. Our model was fit to cholera incidence data from 2010 and was parameterized using demographic data collected from the camp. Because there was no outbreak in Maela following the campaign, it would be difficult to evaluate the effectiveness of the vaccine without a mathematical model. Indeed, our analysis therefore provides new evidence showing that the campaign in Maela was effective. We additionally considered multiple scenarios, including both a theoretical exploration of the effects of variation in timing, effectiveness and supply, as well as the real-world coverage of vaccine in Maela. Our modeling analysis allowed us to examine a wider range of intervention strategies and approaches than would be practically possible in real-world trials, uncovering which strategies are most advantageous under resource limited circumstances. We found that the preferred number of doses per person and timing of vaccination campaigns should be considered in the context of one vs. two dose effectiveness and logistical constraints. Importantly, our analysis coincided with an actual cholera vaccination campaign in the camp and was used to evaluate the campaign and to help determine that there was no need for a follow-up booster campaign. The setting of our analysis is particularly relevant given the recent worldwide increase in total numbers of refugees. Results from our model highlight the utility of vaccination to prevent cholera. Vaccination campaigns can be combined with more permanent water, sanitation, and hygiene infrastructure improvements to reduce the risk of cholera and other enteric disease epidemics. Overall, this study demonstrates that mathematical modeling can generate useful insights into real-time intervention decisions.

Availability of data and materials

The datasets used and/or analysed during the current study are available from the corresponding author on reasonable request.

Abbreviations

CDC:

Centers for Disease Control and Prevention

OCV:

Oral Cholera Vaccine

PU-AMI:

Première Urgence Aide Medicale Internationale

SIWR:

Susceptible-Infectious-Water-Recovered model

VE:

Vaccination Effectiveness

WaSH:

Water, Sanitation and Hygiene

WHO:

World Health Organization

References

  1. Burrows K, Kinney PL. Exploring the climate change, migration and conflict nexus. Int J Environ Res Public Health. 2016; 13(4):443.

    Article  PubMed  PubMed Central  Google Scholar 

  2. UNHCR Figures at a Glance. Available from: http://www.unhcr.org/en-us/figures-at-a-glance.html. Accessed 13 Dec 2019.

  3. Elias CJ, Alexander BH, Sokly T. Infectious disease control in a long-term refugee camp: the role of epidemiologic surveillance and investigation. Am J Public Health. 1990; 80(7):824–8.

    Article  PubMed  PubMed Central  CAS  Google Scholar 

  4. Taylor DL, Kahawita TM, Cairncross S, Ensink JH. The impact of water, sanitation and hygiene interventions to control cholera: a systematic review. PLoS One. 2015; 10(8):e0135676.

    Article  PubMed  PubMed Central  CAS  Google Scholar 

  5. Shultz A, Omollo JO, Burke H, Qassim M, Ochieng JB, Weinberg M, et al.Cholera outbreak in Kenyan refugee camp: risk factors for illness and importance of sanitation. Am J Trop Med Hyg. 2009; 80(4):640–5.

    Article  PubMed  Google Scholar 

  6. World Health Organization, et al.Cholera vaccines: WHO position paper. Wkly Epidemiol Rec. 2010; 85(13):117–28.

    Google Scholar 

  7. Sack RB, Siddique AK. Corpses and the spread of cholera. Lancet. 1998; 352(9140):1570.

    Article  PubMed  CAS  Google Scholar 

  8. World Health Organization, et al.Cholera vaccine: WHO position paper, August 2017–Recommendations. Vaccine. 2018; 36(24):3418–20.

    Article  PubMed  CAS  Google Scholar 

  9. Azman AS, Parker LA, Rumunu J, Tadesse F, Grandesso F, Deng LL, et al.Effectiveness of one dose of oral cholera vaccine in response to an outbreak: a case-cohort study. Lancet Glob Health. 2016; 4(11):e856–63.

    Article  PubMed  Google Scholar 

  10. Ali M, Sur D, You YA, Kanungo S, Sah B, Manna B, et al.Herd protection by a bivalent killed whole–cell oral cholera vaccine in the slums of Kolkata, India. Clin Infect Dis. 2013; 56(8):1123–31.

    Article  PubMed  CAS  Google Scholar 

  11. Bhattacharya SK, Sur D, Ali M, Kanungo S, You YA, Manna B, et al.5 year efficacy of a bivalent killed whole-cell oral cholera vaccine in Kolkata, India: a cluster-randomised, double-blind, placebo-controlled trial. Lancet Inf Dis. 2013; 13(12):1050–6.

    Article  CAS  Google Scholar 

  12. Thailand Burma Border Consortium; Burmese border refugee sites with population figures. 2009. Available from: http://www.theborderconsortium.org/media/11846/2009-12-dec-map-tbbc-unhcr-1-.pdf. Accessed 13 Dec 2019.

  13. World Health Organization; Cholera Map. 2016. Available from: http://gamapserver.who.int/mapLibrary/Files/Maps/Global_Cholera_2010_2015.png. Accessed 13 Dec 2019.

  14. Wiwanitkit V. Cholera outbreak in Thailand during the past 25-year period, a summary on epidemiology. Rev Esp Enferm Dig organo oficial de la Sociedad Espanola de Patologia Digestiva. 2008; 100(4):244–5.

    CAS  Google Scholar 

  15. Phares CR, Date K, Travers P, Déglise C, Wongjindanon N, Ortega L, et al.Mass vaccination with a two-dose oral cholera vaccine in a long-standing refugee camp, Thailand. Vaccine. 2016; 34:128–33.

    Article  PubMed  Google Scholar 

  16. Martin S, Lopez AL, Bellos A, Deen J, Ali M, Alberti K, et al.Post-licensure deployment of oral cholera vaccines: a systematic review. Bull World Health Organ. 2014; 92(12):881–93.

    Article  PubMed  PubMed Central  Google Scholar 

  17. World Health Organization, et al.WHO Consultation on oral cholera vaccine (OCV) stockpile strategic framework: potential objectives and possible policy options: 18-20 September 2011, Geneva, Switzerland. Geneva: World Health Organization; 2012.

    Google Scholar 

  18. Luquero FJ, Grout L, Ciglenecki I, Sakoba K, Traore B, Heile M, et al.Use of Vibrio cholerae vaccine in an outbreak in Guinea. N Engl J Med. 2014; 370(22):2111–20.

    Article  PubMed  CAS  Google Scholar 

  19. Ivers LC, Hilaire IJ, Teng JE, Almazor CP, Jerome JG, Ternier R, et al.Effectiveness of reactive oral cholera vaccination in rural Haiti: a case-control study and bias-indicator analysis. Lancet Glob Health. 2015; 3(3):e162–8.

    Article  PubMed  PubMed Central  Google Scholar 

  20. Qadri F, Ali M, Lynch J, Chowdhury F, Khan AI, Wierzba TF, et al.Efficacy of a single–dose regimen of inactivated whole-cell oral cholera vaccine: results from 2 years of follow–up of a randomised trial. Lancet Inf Dis. 2018; 18(6):666–74.

    Article  CAS  Google Scholar 

  21. Wierzba TF, Kar SK, Mogasale VV, Kerketta AS, You YA, Baral P, et al.Effectiveness of an oral cholera vaccine campaign to prevent clinically–significant cholera in Odisha State, India. Vaccine. 2015; 33:2463–9.

    Article  PubMed  Google Scholar 

  22. Legros D, Paquet C, Perea W, Marty I, Mugisha NK, Royer H, et al.Mass vaccination with a two-dose oral cholera vaccine in a refugee camp. Bull World Health Organ. 1999; 77(10):837.

    PubMed  PubMed Central  CAS  Google Scholar 

  23. Msyamboza KP, Hausi H, Chijuwa A, Nkukumila V, Kubwalo HW, Im J, et al.Feasibility and acceptability of oral cholera vaccine mass vaccination campaign in response to an outbreak and floods in Malawi. Pan Afr Med J. 2016; 23. https://doi.org/10.11604/pamj.2016.23.203.8346.

  24. Heesterbeek H, Anderson RM, Andreasen V, Bansal S, De Angelis D, Dye C, et al.Modeling infectious disease dynamics in the complex landscape of global health. Science. 2015; 347(6227):aaa4339.

    Article  PubMed  PubMed Central  CAS  Google Scholar 

  25. Date KA, Vicari A, Hyde TB, Mintz E, Danovaro-Holliday MC, Henry A, et al.Considerations for oral cholera vaccine use during outbreak after earthquake in Haiti, 2010–2011. Emerg Infect Dis. 2011; 17(11):2105.

    Article  PubMed  PubMed Central  Google Scholar 

  26. Chao DL, Halloran ME, Longini IM. Vaccination strategies for epidemic cholera in Haiti with implications for the developing world. Proc Natl Acad Sci. 2011; 108(17):7081–5.

    Article  PubMed  PubMed Central  Google Scholar 

  27. Codeço CT. Endemic and epidemic dynamics of cholera: the role of the aquatic reservoir. BMC Infect Dis. 2001; 1(1):1.

    Article  PubMed  PubMed Central  Google Scholar 

  28. Tien JH, Earn DJ. Multiple transmission pathways and disease dynamics in a waterborne pathogen model. Bull Math Biol. 2010; 72(6):1506–33.

    Article  PubMed  Google Scholar 

  29. Mukandavire Z, Liao S, Wang J, Gaff H, Smith DL, Morris JG. Estimating the reproductive numbers for the 2008–2009 cholera outbreaks in Zimbabwe. Proc Natl Acad Sci. 2011; 108(21):8767–72.

    Article  PubMed  PubMed Central  Google Scholar 

  30. Eisenberg MC, Robertson SL, Tien JH. Identifiability and estimation of multiple transmission pathways in cholera and waterborne disease. J Theor Biol. 2013; 324:84–102.

    Article  PubMed  Google Scholar 

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

    Article  PubMed  PubMed Central  Google Scholar 

  32. Tuite AR, Tien J, Eisenberg M, Earn DJ, Ma J, Fisman DN. Cholera epidemic in Haiti, 2010: using a transmission model to explain spatial spread of disease and identify optimal control interventions. Ann Intern Med. 2011; 154(9):593–601.

    Article  PubMed  Google Scholar 

  33. Bertuzzo E, Finger F, Mari L, Gatto M, Rinaldo A. On the probability of extinction of the Haiti cholera epidemic. Stoch Env Res Risk A. 2016; 30(8):2043–55.

    Article  Google Scholar 

  34. Brauer F, Shuai Z, Van Den Driessche P. Dynamics of an age-of-infection cholera model. Math Biosci Eng. 2013; 10(5-6):1335–49.

    PubMed  Google Scholar 

  35. Fister KR, Gaff H, Lenhart S, Numfor E, Schaefer E, Wang J. Optimal control of vaccination in an age-structured cholera model. Springer; 2016. pp. 221–48. https://doi.org/10.1007/978-3-319-40413-4_14.

    Chapter  Google Scholar 

  36. Pascual M, Rodó X, Ellner SP, Colwell R, Bouma MJ. Cholera dynamics and El Nino-southern oscillation. Science. 2000; 289(5485):1766–9.

    Article  PubMed  CAS  Google Scholar 

  37. Rinaldo A, Bertuzzo E, Mari L, Righetto L, Blokesch M, Gatto M, et al.Reassessment of the 2010–2011 Haiti cholera outbreak and rainfall-driven multiseason projections. Proc Natl Acad Sci. 2012; 109(17):6602–7.

    Article  PubMed  PubMed Central  Google Scholar 

  38. King AA, Ionides EL, Pascual M, Bouma MJ. Inapparent infections and cholera dynamics. Nature. 2008; 454(7206):877.

    Article  PubMed  CAS  Google Scholar 

  39. Eisenberg MC, Kujbida G, Tuite AR, Fisman DN, Tien JH. Examining rainfall and cholera dynamics in Haiti using statistical and dynamic modeling approaches. Epidemics. 2013; 5(4):197–207.

    Article  PubMed  Google Scholar 

  40. Hartley DM, Morris Jr JG, Smith DL. Hyperinfectivity: a critical element in the ability of V. cholerae to cause epidemics?. PLoS Med. 2005; 3(1):e7.

    Article  PubMed  PubMed Central  Google Scholar 

  41. Dunworth JB. Nonlinear incidence of waterborne diseases. Columbus: The Ohio State University; 2011.

    Google Scholar 

  42. Lee EC, Kelly MR, Ochocki BM, Akinwumi SM, Hamre KE, Tien JH, et al.Model distinguishability and inference robustness in mechanisms of cholera transmission and loss of immunity. J Theor Biol. 2017; 420:68–81.

    Article  PubMed  PubMed Central  Google Scholar 

  43. Azman AS, Luquero FJ, Ciglenecki I, Grais RF, Sack DA, Lessler J. The impact of a one-dose versus two-dose oral cholera vaccine regimen in outbreak settings: a modeling study. PLoS Med. 2015; 12(8):e1001867.

    Article  PubMed  PubMed Central  CAS  Google Scholar 

  44. Crooks AT, Hailegiorgis AB. An agent-based modeling approach applied to the spread of cholera. Environ Mod Softw. 2014; 62:164–77.

    Article  Google Scholar 

  45. Peak CM, Reilly AL, Azman AS, Buckee CO. Prolonging herd immunity to cholera via vaccination: Accounting for human mobility and waning vaccine effects. PLoS Negl Trop Dis. 2018; 12(2):e0006257.

    Article  PubMed  PubMed Central  Google Scholar 

  46. Eisenberg MC, Shuai Z, Tien JH, Van den Driessche P. A cholera model in a patchy environment with water and human movement. Math Biosci. 2013; 246(1):105–12.

    Article  PubMed  Google Scholar 

  47. Alexanderian A, Gobbert MK, Fister KR, Gaff H, Lenhart S, Schaefer E. An age-structured model for the spread of epidemic cholera: analysis and simulation. Nonlinear Anal Real World Appl. 2011; 12(6):3483–98.

    Article  Google Scholar 

  48. Dimitrov DT, Troeger C, Halloran ME, Longini IM, Chao DL. Comparative effectiveness of different strategies of oral cholera vaccination in Bangladesh: a modeling study. PLoS Negl Trop Dis. 2014; 8(12):e3343.

    Article  PubMed  PubMed Central  Google Scholar 

  49. R Core Team. R: A Language and Environment for Statistical Computing. Vienna, Austria. 2017. Available from: https://www.R-project.org/.

  50. Hijmans RJ, van Etten J, Cheng J, Mattiuzzi M, Sumner M, Greenberg JA, Lamigueiro OP, et al.Package ‘raster’. R package. 2015.

  51. Global Administrative Areas; GADM database of Global Administrative Areas, version 2.0. 2012. Available from: https://www.gadm.org.

  52. Shanchol Package insert. Available from: http://www.who.int/immunization_standards/vaccine_quality/pq_250_cholera_1dose_shantha_insert.pdf?ua=1. Accessed 13 Dec 2019.

  53. Weil AA, Khan AI, Chowdhury F, LaRocque RC, Faruque A, Ryan ET, et al.Clinical outcomes in household contacts of patients with cholera in Bangladesh. Clin Infect Dis. 2009; 49(10):1473–9.

    Article  PubMed  Google Scholar 

  54. World Health Organization; Cholera. 2017. Available from: http://www.who.int/mediacentre/factsheets/fs107/en/. Accessed 13 Dec 2019.

  55. Audoly S, Bellu G, D’Angio L, Saccomani MP, Cobelli C. Global identifiability of nonlinear models of biological systems. IEEE Trans Biomed Eng. 2001; 48(1):55–65.

    Article  PubMed  CAS  Google Scholar 

  56. Evans ND, White LJ, Chapman MJ, Godfrey KR, Chappell MJ. The structural identifiability of the susceptible infected recovered model with seasonal forcing. Math Biosci. 2005; 194(2):175–97.

    Article  PubMed  Google Scholar 

  57. Chapman JD, Evans ND. The structural identifiability of susceptible–infective–recovered type epidemic models with incomplete immunity and birth targeted vaccination. Biomed Sig Process Control. 2009; 4(4):278–284.

    Article  Google Scholar 

  58. Bellman R, Åström KJ. On structural identifiability. Math Biosci. 1970; 7(3-4):329–339.

    Article  Google Scholar 

  59. Chis OT, Banga JR, Balsa-Canto E. Structural identifiability of systems biology models: a critical comparison of methods. PloS One. 2011; 6(11):e27755.

    Article  PubMed  PubMed Central  CAS  Google Scholar 

  60. Raue A, Kreutz C, Maiwald T, Bachmann J, Schilling M, Klingmüller U, et al.Structural and practical identifiability analysis of partially observed dynamical models by exploiting the profile likelihood. Bioinformatics. 2009; 25(15):1923–9.

    Article  PubMed  CAS  Google Scholar 

  61. Eisenberg MC, Jain HV. A confidence building exercise in data and identifiability: Modeling cancer chemotherapy as a case study. J Theor Biol. 2017; 431:63–78.

    Article  PubMed  PubMed Central  Google Scholar 

  62. Macdonald G. The analysis of equilibrium in malaria. Trop Dis Bull. 1952; 49(9):813.

    PubMed  CAS  Google Scholar 

  63. Heesterbeek J, Dietz K. The concept of Ro in epidemic theory. Statistica Neerlandica. 1996; 50(1):89–110.

    Article  Google Scholar 

  64. Van den Driessche P, Watmough J. Reproduction numbers and sub-threshold endemic equilibria for compartmental models of disease transmission. Math Biosci. 2002; 180(1):29–48.

    Article  PubMed  CAS  Google Scholar 

  65. Charles RC, Hilaire IJ, Mayo-Smith LM, Teng JE, Jerome JG, Franke MF, et al.Immunogenicity of a killed bivalent (O1 and O139) whole cell oral cholera vaccine, Shanchol, in Haiti. PLoS Negl Trop Dis. 2014; 8(5):e2828.

    Article  PubMed  PubMed Central  Google Scholar 

  66. Saha A, Chowdhury MI, Khanam F, Bhuiyan MS, Chowdhury F, Khan AI, et al.Safety and immunogenicity study of a killed bivalent (O1 and O139) whole-cell oral cholera vaccine Shanchol, in Bangladeshi adults and children as young as 1 year of age. Vaccine. 2011; 29(46):8285–92.

    Article  PubMed  Google Scholar 

  67. Bi Q, Ferreras E, Pezzoli L, Legros D, Ivers LC, Date K, et al.Protection against cholera from killed whole-cell oral cholera vaccines: a systematic review and meta-analysis. Lancet Infect Dis. 2017; 17(10):1080–8.

    Article  PubMed  PubMed Central  CAS  Google Scholar 

  68. MATLAB and Statistics Toolbox Release. Natick: The MathWorks, Inc.; 2016.

  69. World Health Organization; Cholera Fact Sheet. Available from: http://www.who.int/mediacentre/factsheets/fs107/en/. Accessed 13 Dec 2019.

  70. Levine M, Black R, Clements M, Cisneros L, Nalin D, Young C. Duration of infection-derived immunity to cholera. J Infect Dis. 1981; 143(6):818–20.

    Article  PubMed  CAS  Google Scholar 

  71. Koelle K, Pascual M. Disentangling extrinsic from intrinsic factors in disease dynamics: a nonlinear time series approach with an application to cholera. Am Nat. 2004; 163(6):901–13.

    Article  PubMed  Google Scholar 

  72. Atlas RM. Sanitation and Disease-Health Aspects of Excreta and Waste-water Management. JAWRA J Am Water Resour Assoc. 1984; 20(5):803.

    Article  Google Scholar 

  73. Feachem R, Bradley D, Garelick H, Mara D. Sanitation and Disease-Health Aspects of Excreta and Waste-water Management. Vol. 501. Hoboken: Chichester: John Whiley & Sons; 1983.

    Google Scholar 

  74. Xu HS, Roberts N, Singleton F, Attwell R, Grimes D, Colwell R. Survival and viability of nonculturable Escherichia coli and Vibrio cholerae in the estuarine and marine environment. Microb Ecol. 1982; 8(4):313–23.

    Article  PubMed  CAS  Google Scholar 

  75. Nelson EJ, Harris JB, Morris Jr JG, Calderwood SB, Camilli A. Cholera transmission: the host, pathogen and bacteriophage dynamic. Nat Rev Microbiol. 2009; 7(10):693.

    Article  PubMed  CAS  Google Scholar 

  76. Qadri F, Wierzba TF, Ali M, Chowdhury F, Khan AI, Saha A, et al.Efficacy of a single-dose, inactivated oral cholera vaccine in Bangladesh. N Engl J Med. 2016; 374(18):1723–32.

    Article  PubMed  CAS  Google Scholar 

  77. 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-randomised, double-blind, placebo-controlled trial. The Lancet. 2009; 374:1694–702.

    Article  CAS  Google Scholar 

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

    Article  PubMed  PubMed Central  CAS  Google Scholar 

  79. Stoica P, SÖDERSTRÖM T. On the parsimony principle. Int J Control. 1982; 36(3):409–18.

    Article  Google Scholar 

  80. Conway JM, Tuite AR, Fisman DN, Hupert N, Meza R, Davoudi B, et al.Vaccination against 2009 pandemic H1N1 in a population dynamical model of Vancouver, Canada: timing is everything. BMC Public Health. 2011; 11(1):932.

    Article  PubMed  PubMed Central  Google Scholar 

  81. WHO, et al.Guidance on how to access the oral cholera vaccine (OCV) from the ICG emergency stockpile. Geneva: Switzerland: World Health Organization.

  82. World Health Organization et al.International coordinating group on vaccine provision for cholera control: annual meeting, 12–13 July 2016: World Health Organization; 2017.

  83. Moore SM, Lessler J. Optimal allocation of the limited oral cholera vaccine supply between endemic and epidemic settings. J R Soc Interface 12.111. 2015:20150703.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  84. Fung ICH, Fitter DL, Borse RH, Meltzer MI, Tappero JW. Modeling the effect of water, sanitation, and hygiene and oral cholera vaccine implementation in Haiti. Am J Trop Med Hyg. 2013; 89(4):633–40.

    Article  PubMed  Google Scholar 

Download references

Acknowledgements

We would like to thank Dr. Nuttapong Wongjindanon Thailand Ministry of Public Health – U.S. Centers for Disease Control and Prevention Collaboration, Nonthaburi, Thailand for his helpful comments, advice, and support. Furthermore, we would like to thank PU-AMI and Niamh de Loughry in Maesot, Thailand, for providing us with key data and information for our analyses and for helping us gain access to and navigate Maela.

Funding

This work was supported by the National Science Foundation grant OCE-1115881 (to JH, and MCE), the National Institute of General Medical Sciences of the National Institutes of Health under Award Number U01GM110712 (JH, RM, and MCE), and the International Health Travel Award Grant from the Department of Epidemiology, University of Michigan, Ann Arbor (JH). The funding bodies had no role in the design of the study and collection, analysis, and interpretation of data or writing the manuscript.

Author information

Authors and Affiliations

Authors

Contributions

JH, RM and MCE came up with the research plan, analyzed the data, and wrote the paper. CRP provided guidance on the dataset and KD and CRP provided guidance on the analysis plan and paper. All authors provided edits and approved the final manuscript.

Corresponding author

Correspondence to Marisa C. Eisenberg.

Ethics declarations

Ethics approval and consent to participate

This study was determined to be nonregulated by the University of Michigan institutional review board (HUM00061865) because only deidentified data were used. Therefore, the need for consent was not required for this study, however in the original data collection (see [15] for more details), entry screeners obtained consent verbally (as illiteracy is high). Deidentified data was shared with the authors as part of an ongoing collaboration between JH, RM, and MCE with the CDC and its partners. However, CDC and its partners retain ownership of the data and approved its use in publication.

Consent for publication

Not Applicable.

Competing interests

The authors declare that they have no competing interests.

Additional information

Publisher’s Note

Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Supplementary information

Rights and permissions

Open Access This article is distributed under the terms of the Creative Commons Attribution 4.0 International License (http://creativecommons.org/licenses/by/4.0/), which permits unrestricted use, distribution, and reproduction in any medium, provided you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated.

Reprints and permissions

About this article

Check for updates. Verify currency and authenticity via CrossMark

Cite this article

Havumaki, J., Meza, R., Phares, C.R. et al. Comparing alternative cholera vaccination strategies in Maela refugee camp: using a transmission model in public health practice. BMC Infect Dis 19, 1075 (2019). https://doi.org/10.1186/s12879-019-4688-6

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: https://doi.org/10.1186/s12879-019-4688-6

Keywords