|Home | About | Journals | Submit | Contact Us | Français|
Conceived and designed the experiments: AA FJL RFG BTG JL. Performed the experiments: FJL AR PPP RFG CNB. Analyzed the data: AA FJL BTG JL. Wrote the paper: AA FJL AR PPP RFG CNB BTG JL.
Use of cholera vaccines in response to epidemics (reactive vaccination) may provide an effective supplement to traditional control measures. In Haiti, reactive vaccination was considered but, until recently, rejected in part due to limited global supply of vaccine. Using Bissau City, Guinea-Bissau as a case study, we explore neighborhood-level transmission dynamics to understand if, with limited vaccine and likely delays, reactive vaccination can significantly change the course of a cholera epidemic.
We fit a spatially explicit meta-population model of cholera transmission within Bissau City to data from 7,551 suspected cholera cases from a 2008 epidemic. We estimated the effect reactive vaccination campaigns would have had on the epidemic under different levels of vaccine coverage and campaign start dates. We compared highly focused and diffuse strategies for distributing vaccine throughout the city. We found wide variation in the efficiency of cholera transmission both within and between areas of the city. “Hotspots”, where transmission was most efficient, appear to drive the epidemic. In particular one area, Bandim, was a necessary driver of the 2008 epidemic in Bissau City. If vaccine supply were limited but could have been distributed within the first 80 days of the epidemic, targeting vaccination at Bandim would have averted the most cases both within this area and throughout the city. Regardless of the distribution strategy used, timely distribution of vaccine in response to an ongoing cholera epidemic can prevent cases and save lives.
Reactive vaccination can be a useful tool for controlling cholera epidemics, especially in urban areas like Bissau City. Particular neighborhoods may be responsible for driving a city's cholera epidemic; timely and targeted reactive vaccination at such neighborhoods may be the most effective way to prevent cholera cases both within that neighborhood and throughout the city.
Cholera remains a major public health threat, causing 3–5 million cases and 100,000–120,000 deaths each year. In 2010, data on vaccine performance and the changing epidemiology of cholera prompted the WHO's Strategic Advisory Group to recommend that reactive vaccination be considered in specific areas. We built a spatially explicit stochastic model of cholera transmission and fit it to data from a 2008 epidemic in Bissau City, Guinea Bissau. Using this model we examined the potential effectiveness of reactive vaccination for controlling cholera transmission in Bissau City, comparing strategies for distributing limited vaccine. In simulations, early targeting of a single transmission “hotspot”, Bandim, was the most effective strategy, and led to the greatest reduction in cases both within Bandim and in areas where no vaccine was distributed. This finding has implications for cholera control in urban settings in general: public health officials will often know which areas of a city were hotspots of cholera transmission in the past or where conditions promote efficient transmission. When there is limited vaccine, our work suggests that targeting reactive vaccination at these areas will lead to the greatest reduction in cases both in these areas and elsewhere in the city.
With the introduction of inexpensive, easy to administer, and effective oral vaccines against cholera, vaccination in response to an epidemic (reactive vaccination) may be an effective supplement to conventional control measures. Two safe and internationally licensed oral cholera vaccines are currently available, Dukoral and Shanchol. Both protect against clinical cholera two or more years after vaccination, but neither confers long lasting immunity –. On an epidemic timescale, these vaccines have efficacies ranging from 66 to 86% , .
Vaccination against cholera has been used preventatively , –, but before 2012, we know of only two instances, in The Federated States of Micronesia in 2000 and Vietnam in 2008, where vaccination commenced during an epidemic , . Vaccine efficacy estimates ranged from 76 to 80%, however, no analysis on how vaccination affected the course of the epidemic was reported for either case , .
New data on vaccine performance and the changing epidemiology of cholera prompted the WHO's Strategic Advisory Group to recommend in 2010 that reactive vaccination be considered in specific areas . In order to facilitate rapid procurement and deployment of an oral cholera vaccine, some have proposed the creation of a revolving global stockpile , . While discussions of the global stockpile proceed, countries that use reactive vaccination must contend with a limited supply that may arrive after a significant delay.
Spatial heterogeneities may influence how cholera vaccine can best be distributed in a reactive campaign. The effectiveness of a campaign and optimal allocation strategy will depend upon local cholera transmission dynamics, vaccine supply, and logistical delays , . Human movement, water and sewerage infrastructure, and natural waterways facilitate cholera transmission across a city. Within neighborhoods, there can be marked variation in the efficiency of transmission.
One country that may benefit from reactive vaccination is Guinea-Bissau, where outbreaks have occurred every three to four years since 1994. Sector Autónomo de Bissau (SAB), or Bissau City, the capital, consistently reports the most cholera cases within the country (unpublished data, Guinea-Bissau Ministry of Health). In 2008, 67% of reported cases occurred in SAB while only 25% of the national population live within its boundaries . Reactive vaccination in SAB may be possible in future epidemics given the concentration of cases within the city and the Ministry of Health's experience with vaccination campaigns.
Here, we explore the possible effectiveness of different reactive vaccination strategies using SAB as a case study. We fit a neighborhood-based meta-population model to the 2008 cholera epidemic. Using this model, we characterize the spatio-temporal dynamics of cholera transmission within the city and estimate the impact that different reactive vaccination strategies could have had on the course of the epidemic.
During the 2008 epidemic, the Guinea-Bissau Ministry of Health, the WHO, and Mèdecins Sans Frontières implemented a clinic-based cholera surveillance system, which has been described previously . In brief, upon arrival at either the cholera treatment center in the Hospital National Simao Mendes or one of five cholera treatment units (Figure 1C and 1D), health care providers entered patients into a surveillance registry. A patient's age, sex, area of residence, treatment facility, date of presentation, and clinical diagnosis were recorded.
Modified WHO cholera case definitions were used . A suspected case was any person suffering from acute watery diarrhea, and a confirmed case was a suspected case with a positive stool sample containing Vibrio cholerae O1 or O139. We included all suspected and confirmed cases with complete information on their presentation date and home sanitary area in this analysis. The population for each sanitary area within the city was extrapolated from 1991 census data using a constant linear growth rate estimated by the Direcção-Geral Saúde. To estimate the population density in each sanitary area we traced the residential areas using Google Earth (v220.127.116.117), then divided each sanitary area's population by its estimated residential area.
We fit a discrete-time Susceptible-Infectious-Recovered meta-population model to the confirmed and suspected cases reported during the 2008 epidemic with each of 14 sanitary areas in SAB treated as a distinct population. We assume the epidemic follows a first-order Markov process with a fixed generation time of five days. At each time step, the incidence in each area follows a Poisson distribution with a mean determined by the number infected in the last time step in all areas and the proportion of the area's population remaining susceptible. After infection, individuals were assumed to remain immune for the duration of the epidemic (See Text S1 for model details).
We considered models of cholera transmission with and without seasonality assuming (A) equal transmission coefficients between and within all areas of SAB; (B) different transmission coefficients within each area and equal transmission coefficients between all areas; (C) different transmission coefficients within each area and unique symmetric transmission coefficients between each pair of areas; and, (D) different transmission coefficients within each area and unique asymmetric transmission coefficients between each pair of areas in the city. We chose the best model based on Deviance Information Criteria (Text S1). To assess fit we simulated 300,000 epidemics predicting five, fifteen, and fifty days ahead drawing new parameters from the posterior distribution every 1000 simulations.
Posterior distributions were approximated using Markov Chain Monte Carlo methods using JAGS 3.1.0 and R 2.14.0 with non-informative priors , . We ran 3 chains of 400,000 iterations with a burn in of 50,000, and assessed convergence using the potential scale reduction factor and through visual inspection .
We assume every vaccinated individual receives two doses in a vaccine campaign over a 20 day period and that 75% are fully protected ( ) , , , . In our model vaccinees get no protection until 10 days after the second dose , . Hence, 75% of the susceptible vaccinees are considered immune starting 30 days after their first dose, with no protection before (Table 1).
We considered campaigns with 50,000, 75,000, or 100,000 doses (i.e. 25,000, 37,500, and 50,000 individuals vaccinated) and targeted vaccination at one, two, three, or all (14) areas (Table 3). When the proposed number of vaccinees in a specific area exceeded the population size, we distributed vaccine to the other vaccination areas or, in the campaigns with one vaccination area, we dispersed the vaccine throughout the city with each person having equal probability of getting vaccinated. We varied the starting time of the vaccination campaign between 20 and 120 days after the first case was detected.
We considered targeted and diffuse (city-wide) campaigns. In diffuse campaigns, vaccine was distributed throughout all areas of SAB. In targeted campaigns, we considered three different strategies to select vaccination areas. In the population-based strategy, we selected the areas with the largest population. In the connectivity-based strategy, we vaccinated in areas estimated to be most “connected” to other areas. In the attack rate-based strategy, we chose the areas with the highest attack rate in the 2008 epidemic. We allocated vaccine proportional to population size in all simulations.
For each vaccination scenario we ran 5,000 simulations calculating the difference between the final epidemic size with and without vaccination. Epidemics were assumed to follow the observed 2008 epidemic course until 30 days after the first dose. In each simulation we drew new parameters from the joint posterior distribution. As a sensitivity analysis, we ran simulations with different generation times (3–10 days) and vaccine efficacies (65%–85%). Additional simulation study details are available in Text S1.
Original data collection was approved by the Mèdecins Sans Frontières ERB and the National Ethical Review Board of Guinea-Bissau . The analyses presented in this article were conducted on de-identified data and deemed to be non-human subject research by the Johns Hopkins Bloomberg School of Public Health IRB.
The first case in SAB was reported on June 5, 2008 in Bairro-Militar, the most populated area of the city (Figures 1A, 1B), one month after the first reported case in Guinea-Bissau. Within three weeks, all 14 areas had reported cases (Figure 1C). The Ministry of Health officially declared an epidemic one month after the first case report from SAB. The National Laboratory of Microbiology and the Pasteur Laboratory in Dakar, Senegal identified all positive specimens analyzed as Vibrio cholerae O1 El Tor Ogawa.
Nationally, 14,226 suspected cases and 228 deaths were reported with 67% (9,393) of cases and 32% (73) of deaths reported in SAB. The last case in the country was reported in SAB on January 11, 2009. Individual-level data in SAB was collected between June 5, 2008 and October 28, 2008, over which time 8,024 (85%) suspected and confirmed cases were reported. These analyses focus on 7,551 suspected and confirmed cases with complete information on date of presentation, home area, and clinical diagnosis (Figure S1).
In SAB, weekly incidence ranged from 14 to 755. Within-area attack rates ranged from 9.1 to 40.6 per 1,000 (Table 2, Figure 1D), with Bandim having both the most cases (1,816) and the highest attack rate.
The final model fit both the overall and area-specific epidemic curves well, even when predicting as far as 50 days (i.e. 10 time steps) ahead (Figures 2A,2B). To understand how transmission varied through time, we calculated the odds that an incident case was caused locally (i.e. attributable to transmission between people in the same area) for each area throughout the course of the epidemic (Figure 3). Only Bandim, Plaque, and Santa-Luzia have an odds consistently greater than 1, suggesting internally driven epidemics in these areas.
We define the effective internal basic reproductive number () as the expected number of cases caused within a given area by one infected individual, within the same area, at the beginning of the epidemic. Only areas with can sustain an epidemic absent infections introduced from other areas. The strength of internal epidemics varied with estimates of ranging from 0.01 (95% Credible Interval (CI) 0.00–0.07) in Ajuda to 1.17 (95% CI 0.99–1.33) in Bandim (Figure 4). We found no significant correlation between and either estimated population size or population density.
Bandim is the only area where we estimate , and it appears to have played a necessary role in driving the epidemic. With Bandim removed, simulated introductions of cases fail to cause epidemics. In contrast, city-wide epidemics occur with removal of any other single area.
In simulated epidemics based upon our best-fit model, we find that, on average, at least 10% of cases in each area are caused by cases in other areas (Figure 2C, Text S1). External transmission coefficients represent epidemic connectivity between areas, and our estimates suggest heterogeneity in inter-area transmission (Text S1). Based on simulations, we estimate that Bandim contributed over 10% of the cases to over half (7/13) of the other areas (Figure 2D), highlighting the crucial role it played in the epidemic.
The sum of the external transmission coefficients for any area provides an estimate of the effective external basic reproductive number (). This number is the estimated number of cases a single infectious case in that area would cause in all other areas of SAB given the pre-epidemic level of population immunity. Estimates of ranged from 0.37 (95% CI 0.16–0.71) in Belem to 7.32 (95% CI 6.29–8.37) in Missira (Figure 4).
Vaccination in the area(s) with the highest attack rate leads to larger reduction in cases than all other targeted and city-wide campaigns at all starting times. Targeting vaccination at Bandim only, the area with the highest attack rate, within the first 80 days of the epidemic averts more cases than other strategies regardless of vaccine quantity (Figure 5). Targeted vaccination in Bandim starting on day 20 is expected to reduce the final size of the epidemic by 41% (95% Predictive Interval (PI) 0.21–0.69), 56% (95% PI 0.30–0.85), and 67% (95% PI 0.40–0.89) with 25,000, 37,500, and 50,000 vaccinees, respectively. In comparison, a city-wide campaign starting on the same day is expected to reduce the epidemic size by 21% (95% PI 0.07–0.34), 30% (95% PI 0.17–0.44), and 40% (95% PI 0.27–0.55) for 25,000, 37,500, and 50,000 vaccinees (Tables 4,S1,S2).
We found wide variability in the outcomes using different targeting strategies, with the differences diminishing as vaccination is delayed (Figure 5). Under the population-based strategy, only a targeted campaign in the three most populated areas averts more cases than a city-wide campaign (Figure 5, Table 4). Targeting the areas estimated to be most “connected” to others averts fewer cases than city-wide campaigns regardless of vaccination starting time and doses.
Starting day has a profound impact on the effect of all vaccination campaigns: the sooner vaccination begins, the more cases are averted. With 37,500 vaccinees, each day delay in vaccination results in an average of 39.5 (95% CI 37.7–44.2) fewer cases averted when targeting based on attack rate. Increasing the size of a vaccination campaign early on in the epidemic can significantly improve case prevention, however, the marginal benefit of additional vaccine diminishes as vaccination is delayed. On average, each additional person vaccinated as part of a targeted campaign in Bandim starting on day 20 averts 7.5 cases compared to 1.7 cases averted per vaccinee in campaigns starting two months later.
In simulations, early targeted vaccination leads to fewer cases both within the targeted area and throughout the city when compared to diffuse campaigns. When starting vaccination on day 20 (Figure 6A), targeting Bandim averts more cases both in Bandim (1,173) and in all the other areas combined (2,265) when compared to a city-wide campaign (341 averted in Bandim and 1,741 in all other areas). As the vaccination campaign is delayed, these differences shrink (Figure 6).
Using a simple spatially explicit model of cholera transmission, we captured the essential dynamics of the 2008 cholera epidemic in SAB, Guinea-Bissau. This model suggests that there was significant transmission between areas in SAB and that one area, Bandim, drove the epidemic. Our simulations show that early distribution of vaccine is the most important determinant of the number of cases prevented. For example, vaccinating 25,000 individuals in Bandim on epidemic day 20 would have averted more cases (3,109, 95% PI 1,475–5,198) than vaccinating 50,000 in the same area just 40 days later (2,732, 95% PI 1,630–3,738). Our simulations suggest that an early vaccination campaign targeted at Bandim alone would have outperformed distributing the same vaccine quantity throughout the city. Not only are more cases prevented overall, but more are prevented in areas outside of Bandim.
Our results suggest that rapid small-scale vaccination may be more effective than a delayed larger-scale vaccination campaign. For example, on average, each day delay results in an additional 39.5 cases when targeting 37,500 people in the areas with the highest attack rate. Applying the average case fatality ratio from the 2008 epidemic (1.58 per 100 cases ) we estimate that each week delay in vaccination would have resulted in an average of 4.4 cholera-related deaths.
Transmission hotspots for other infectious diseases have been exploited to devise novel prevention and control approaches , . For example, targeted interventions in hotspots may be key to effective malaria control and elimination . Similarly, cholera hotspots can serve as targets for both reactive and preventative interventions. Identification of hotspots during an epidemic may be challenging. In the case of SAB, Bandim is an area which has had high attack rates in previous epidemics and few improvements in water and sanitation infrastructure. Such historical information may be useful in targeting vaccination; however, more research on combining historical and real-time surveillance data is needed.
In our model, vaccination campaigns lasted 20 days, but in reality the duration will vary by the number of vaccinees targeted and the vaccine used. If Shanchol were used with the recommended inter-dose period of 14 days, the campaign would likely exceed 20 days. While this suggests that our results underestimate the speed by which Shanchol vaccination would occur, these differences would be offset by partial immunity conferred before a second dose .
As the time to distribute vaccine doses increases, we expect to avert fewer cases. However, there is some evidence that a single dose of oral cholera vaccine may be sufficient for reactive vaccination , . If one dose is sufficient to elicit a strong protective response for the time-scale of an epidemic, more people could be vaccinated quickly.
Cholera's generation time is not well characterized and varies widely with the concentration of bacteria in the environment, its survival rate, and the route of transmission –. We ran analyses with alternate generation times of 3, 7, and 10 days and got the same qualitative results (Figures S3, S4, S5, S6, S7, S8 and Tables S3, S4, S5). We also found that varying the vaccine efficacy to 65% and 85% changed the number of cases averted, but preserved the relative performance of each strategy over time (Figure S2 and Tables S7,S6).
There are a number of limitations to this work. We focus on a single epidemic in Guinea-Bissau. A longer time series would provide insight into variability in transmission across epidemics. The data came from an intensified surveillance effort from both Mèdecins Sans Frontières and the Guinea-Bissau Ministry of Health, however suspected cases that presented after October 28, 2008 were only captured by the national surveillance system without details on timing and home sanitary area.
There are several possible alternative explanations for the elevated attack rate in Bandim. The cholera case definition used is not 100% specific, and some cholera cases may be false positives. People may be more likely to seek care if their neighbors do, hence clinic visits may cluster even if cholera does not. In addition, Bandim has been the location to several surveillance programs and public health interventions through the Bandim Health Project , perhaps leading to increased awareness. However, if these phenomena were consistent throughout the epidemic they would not lead to elevated estimates of the local transmission rate under our algorithm.
We found that how rapidly vaccine can be distributed during a cholera epidemic is the most important determinant of the effectiveness of a reactive vaccination program; and that a single area of SAB was an essential driver of the epidemic. Hence, early targeting of this area would have been the most effective way to reactively distribute vaccine. These results may apply to urban cholera epidemics more generally. It seems reasonable that cholera epidemics in other urban settings, particularly in Africa, may be disproportionally driven by specific parts of the city. If these hotspots can be identified, targeted reactive vaccination may be an effective way to prevent cases both within that area and throughout the city, especially when vaccine supply is limited. Regardless of the distribution strategy used, timely distribution of vaccine in response to an ongoing cholera epidemic can prevent cases and save lives.
5-day aggregated case counts for all sanitary areas during the 2008 epidemic. Data collected from cholera treatment center and cholera treatment units throughout the city from June 5, 2008 to October 28, 2008.
Vaccine efficacy sensitivity analysis.Comparison of proportion of epidemic averted with different 65%, 75% (as in main analysis), and 85% vaccine efficacy over different vaccination starting times. All scenarios shown use attack rate based targeting.
Comparison of transmission parameters with different generation times. Posterior means and standard deviation for transmission coefficients, ('s on diagonals and 's on off-diagonals) with 3, 5, and 7 day generation times.
Comparison of internal and external effective reproductive numbers for different generation time aggregations.
Proportion of cases caused in each area by others from 3, 5, 7, and 10-day generation time models. The sum of each row is equal to one, representing 100% of the area's epidemic.
Vaccination simulation results with 3-day generation time, 75% vaccine efficacy, and 75,000 doses.
Vaccination simulation results with 7-day generation time, 75% vaccine efficacy, and 75,000 doses.
Vaccination simulation results with 10-day generation time, 75% vaccine efficacy, and 75,000 doses.
Vaccination simulation results with 50,000 doses and 75% vaccine efficacy. Proportion and number of cases averted in 5,000 simulations under different vaccination strategies (Median and 95% Predictive Interval).
Vaccination simulation results with 100,000 doses and 75% vaccine efficacy. Proportion and number of cases averted in 5,000 simulations under different vaccination strategies (Median and 95% Predictive Interval).
Vaccination simulation results from 3-day generation time model, 75,000 doses. Proportion and number of cases averted in 5,000 simulations under different vaccination strategies (Median and 95% Predictive Interval).
Vaccination simulation results from 7-day generation time model, 75,000 doses. Proportion and number of cases averted in 5,000 simulations under different vaccination strategies (Median and 95% Predictive Interval).
Vaccination simulation results from 10-day generation time model. Proportion and number of cases averted in 5,000 simulations under different vaccination strategies (Median and 95% Predictive Interval).
Vaccination simulation results with 75,000 doses and 65% vaccine efficacy. Proportion and number of cases averted in 5,000 simulations under different vaccination strategies (Median and 95% Predictive Interval).
Vaccination simulation results with 75,000 doses and 85% vaccine efficacy. Proportion and number of cases averted in 5,000 simulations under different vaccination strategies (Median and 95% Predictive Interval).
Details on final model, model selection, and simulations.
ASA and JL would like to thank the JHSPH Infectious Disease Dynamics Group for helpful feedback. The authors would like to thank Laurence Carbonetti and Kenneth C. Cline for the useful comments.