Estimating effectiveness of case-area targeted response interventions against cholera in Haiti
Abstract
Case-area targeted interventions (CATIs) against cholera are conducted by rapid response teams, and may include various activities like water, sanitation, hygiene measures. However, their real-world effectiveness has never been established. We conducted a retrospective observational study in 2015–2017 in the Centre department of Haiti. Using cholera cases, stool cultures and CATI records, we identified 238 outbreaks that were responded to. After adjusting for potential confounders, we found that a prompt response could reduce the number of accumulated cases by 76% (95% confidence interval, 59 to 86) and the outbreak duration by 61% (41 to 75) when compared to a delayed response. An intense response could reduce the number of accumulated cases by 59% (11 to 81) and the outbreak duration by 73% (49 to 86) when compared to a weaker response. These results suggest that prompt and repeated CATIs were significantly effective at mitigating and shortening cholera outbreaks in Haiti.
Introduction
On October 2017, Global Task Force on Cholera Control (GTFCC) partners committed to reduce cholera deaths by 90% and to eliminate the disease in 20 countries by 2030, through a multi-sectoral approach (Anon, 2017). This new global strategy planned to combine long-term disease prevention in cholera hotspots with sustainable WaSH (water sanitation and hygiene) solutions and large-scale use of oral cholera vaccine (OCV), with the short-term strengthening of early detection of outbreaks and immediate and effective response through reactive OCV campaigns and rapid response teams (RRTs) (Anon, 2017). RRTs, also referred as mobile teams, have been successfully implemented against polio or Ebola outbreaks (Global Polio Eradication Initiative, 2017; World Health Organization (WHO), 2014). However, response interventions targeted to neighbours of cholera cases (case-area targeted interventions [CATIs]) using combinations of water, sanitation, and hygiene measures, and/or prophylactic antibiotics have rarely been documented, evaluated or promoted against cholera in the published literature (Voelckel, 1971; Piarroux and Bompangue, 2011; Deepthi et al., 2013; Taylor et al., 2015; Mwambi et al., 2016; Finger et al., 2018).
In practice, CATIs are supported by the frequent household transmission of Vibrio cholerae O1 (Weil et al., 2009; Weil et al., 2014; Taylor et al., 2015; Domman et al., 2018), the increased cholera risk among neighbours living within a few dozen meters of cases during the few days following disease onset (Debes et al., 2016; Azman et al., 2018), and the significant protection of household contacts of cholera patients by promoting hand washing with soap and treatment of water (George et al., 2016). A micro-simulation modelling study suggests that early CATIs can be more resource-efficient than mass interventions against cholera (Finger et al., 2018). However, CATI effectiveness has never been evaluated in a real-world setting.
Haiti has implemented CATIs as a national coordinated strategy against cholera since July 2013 (Rebaudet et al., 2019a). After the disease was accidentally imported in October 2010 (Piarroux et al., 2011), the country experienced a massive epidemic, with a total of 820,085 suspected cases and 9792 cholera-related deaths recorded by April 20, 2019 by the Haitian Ministry of Public Health and Population (MSPP) (http://mspp.gouv.ht/newsite/documentation.php, accessed Jul 1, 2019). In 2013, only 68% of Haitian households drank from improved water sources, 26% had access to improved sanitation facilities and 34% had water and soap available for hand washing (République d’Haïti, Ministère de la Santé Publique (MSPP), 2013). But little of the $1.5 billion USD designated by the Plan for the Elimination of Cholera in Haiti 2013–2022 to develop water and sanitation infrastructures has been expended or pledged so far (Ministry of Public Health and Population, National Directorate for Water Supply and Sanitation, 2013). Two pilot OCV reactive campaigns vaccinated approximately 100,000 people in 2012 and to date, additional campaigns have targeted about 10% of the Haitian population (Ivers, 2017; Poncelet, 2015). UNICEF thus backed the MSPP and the Haitian National Directorate for Water and Sanitation (DINEPA) to launch a complementary nationwide coordinated cholera alert-response strategy aiming to interrupt local cholera outbreaks at an early stage (Rebaudet et al., 2019a). This program planned to rapidly send multisectoral rapid response teams to every patient household and neighbourhood in order to identify additional cases, to decontaminate patient premises, to educate on risk factors and methods of prevention and management, to distribute soap and oral rehydration salts (ORS), to chlorinate water at the household level or directly at collection points, and to propose prophylactic antibiotics to close contacts of cholera cases.
This response CATI strategy was implemented gradually from mid-2013 and became an essential pillar of the fight when the national cholera elimination plan was updated in mid-2016 (République d’Haïti, Ministère de la Santé Publique et de la Population, Direction Nationale de l’Eau Potable et de l’Assainissement, 2016). Implementation of this strategy offers a unique opportunity to evaluate the effectiveness of CATIs against cholera outbreaks. Based on available data, we conducted a retrospective observational study over 3 years in the Centre department of Haiti addressing the outcome of local cholera outbreaks according to the response promptness and intensity. We present here the first effectiveness estimates for rapid and targeted response interventions against cholera.
Results
Outbreak and response characteristics
From January 1, 2015, to December 31, 2017, the line-listing of the Centre department reported a total of 10,931 patients with suspected cholera, including 10,428 with a comprehensive location. Details on cholera cases are summarized in Appendix 1—table 1. Intravenous (IV) rehydration was mentioned for 2144 of them. These patients originated from 1497 localities and their time distribution exhibited a marked seasonality (Figure 1A and B). Concomitantly, 1070 stools sampled in Centre department were cultured for V. cholerae O1, of which 509 (48%) were positive (Figure 1A), including 360 with a comprehensive location. Additional details on cholera cultures are summarized in Appendix 1—table 1.
Defining outbreaks by the occurrence of at least two suspected cholera cases with at least one severely dehydrated case or positive culture, within the same locality, during a three-day time window, and after a refractory period of at least 21 days with no case, we identified 452 cholera outbreaks (Figure 2), which mainly occurred during case incidence peaks (Figure 1C) and were distributed across 290 localities (Figure 3). The median cumulative number of cases per outbreak was 3 (Interquartile range [IQR], 4), and the median duration of outbreaks was 5 days (IQR, 18).
Over the same period, 3,887 CATIs were notified in the Centre department by non-governmental organization (NGO) rapid response teams, including 3,533 CATIs (91%) with a comprehensive location, and 2719 (70%) conducted in tandem with staff of EMIRA (Equipe mobile d’intervention rapide, Rapid intervention mobile team, i.e. cholera rapid response team of the MSPP) (Figure 1D). Based on CATI activities summarized in Appendix 1—table 1, a total of 3,596 CATIs (93%) were categorized as complete (at least decontamination, education and distribution of chlorine tablets), and 1922 (49%) also included a reported antibiotic prophylaxis. Overall, 633 complete CATIs (18%) were conducted in localities experiencing an identified cholera outbreak (Figure 1D).
Analysis of confounders
Baseline characteristics of outbreaks and comparisons between the four classes of response promptness are presented in Table 1. The time to the first complete CATI (response promptness) significantly improved during the six semesters of the study, was significantly higher in more densely populated localities, and was lower in localities targeted by a previous OCV campaign. Outbreaks with prompter responses exhibited significantly more positive cultures during the first three days than outbreaks with delayed responses. None of the other covariates were significantly associated with response intensity (Table 1).
Baseline characteristics of outbreaks and comparisons between the four classes of response intensity are presented in Table 2. The numbers of complete CATIs per week and per case (response intensity) significantly improved during the six semesters of the study. Outbreaks receiving more CATIs per case exhibited significantly fewer cases during the first three days than outbreaks receiving less intense responses (Table 2).
CATI effectiveness according to the response promptness
There was a positive association between the time to the first complete CATI after outbreak onset, and the number of cases recorded from the fourth day of the outbreak (Figure 4A, Table 3). Consequently, the prompter the response, the higher the CATI effectiveness on the reduction of outbreak size (Table 3). Compared to a first complete CATI > 7 days after outbreak onset, the crude effectiveness of a first complete CATI ≤ 1 day (cCE1) was 83% (95% CI, 71 to 90), and after adjusting for potential confounders (aCE1), 76% (59 to 86).
Similarly, there was a positive association between the time to the first complete CATI after outbreak onset, and the duration of outbreaks (Figure 4B, Table 4). Consequently, the prompter the response, the higher the CATI effectiveness on the reduction of outbreak duration (Table 4). Compared to a first complete CATI > 7 days after outbreak onset, the crude effectiveness of a first complete CATI ≤ 1 day (cCE2) was 59% (36 to 74), and after adjusting for potential confounders, (aCE2) 61% (41 to 75).
CATI effectiveness according to the response intensity
In addition, there was a negative association between the number of complete CATIs per week of outbreak, and the number of cases recorded from the fourth day of outbreak (Figure 4C, Table 5). Consequently, the more intense the response, the significantly higher the CATI effectiveness was estimated to be on the reduction of outbreak size (Table 5). Compared to a number of complete CATIs < 0.25 per week, the crude effectiveness of a number of complete CATIs ≥ 1 per week (cCE3) was 74% (95% CI, 44 to 88), and after adjusting for potential confounders (aCE3), 59% (11 to 81).
Similarly, there was a negative association between the number of complete CATIs per case, and the duration of outbreaks (Figure 4D, Table 6). Consequently, the more intense the response, the significantly higher the CATI effectiveness on the reduction of outbreak duration (Table 6). Compared to a number of complete CATIs < 0.25 per case, the crude effectiveness of a number of complete CATIs ≥ 1 per case (cCE4) was 76% (95% CI, 54 to 88), and after adjusting for potential confounders (aCE4), 73% (49 to 86).
Several sensitivity analyses using alternative definitions of cholera outbreak, alternative definitions of CATIs and alternative methods of covariate selection for adjustment yielded consistent estimates of CATI effectiveness according to response promptness and response intensity (Appendix 3).
Effectiveness of antibiotic prophylaxis
Finally, stratified analyses showed that three estimates of CATI effectiveness out of four appeared higher in the subgroup of outbreaks that were only responded to by complete CATIs with antibiotic prophylaxis (ATB) than in the subgroup of outbreaks only responded to by complete CATIs that never included ATB (Table 7). More precisely, the adjusted effectiveness of a prompt response on outbreak size (aCE1) was 63% (24 to 82) when all CATIs included antibiotic prophylaxis, and 39% (-38 to 73) when no CATI did. The adjusted effectiveness of a prompt response on outbreak duration (aCE2) was 74% (43 to 88) when all CATIs included antibiotic prophylaxis, and 58% (11 to 80) when no CATI did. Similarly, the adjusted effectiveness of an intense response on outbreak duration (aCE4) was 90% (72 to 96) when all CATIs included antibiotic prophylaxis, and 79% (46 to 92) when no CATI did. Conversely, the adjusted effectiveness of an intense response on outbreak size (aCE3) was 62% (3 to 85) when all CATIs included antibiotic prophylaxis, and 76% (12 to 94) when no CATI did (Table 7).
Discussion
Our quasi-experimental study, based on epidemiological and intervention records over three years in one administrative department of Haiti, showed that prompt and repeated response CATIs conducted by rapid response teams were significantly associated with shortening of cholera outbreaks and mitigating of outbreak case load. Of note, numerous suspected cholera outbreaks spontaneously ended before any response could be conducted. But when taking into account this significant confounding by indication (Remschmidt et al., 2015), the prompter the first complete CATI was implemented and the more complete CATIs were conducted, the fewer cases were recorded and the shorter the outbreak lasted.
While many mild suspected outbreaks may spontaneously end without any response intervention, prompt and repeated CATIs appear difficult to sustain during the largest outbreaks. As suggested by the slow increase in the number of CATIs observed during the study period, CATIs may be logistically complex to implement, and response teams can become overwhelmed when they try to simultaneously address a large number of cases (Finger et al., 2018; Rebaudet et al., 2019a). Such strategy certainly is most relevant at the beginning of epidemics, or during trough periods or tails of epidemics (Finger et al., 2018; Rebaudet et al., 2013).
Our study comes with a number of limitations. Because CATIs were not randomized, response effectiveness may have been biased by unmeasured confounders. As we observed a significant confounding by indication on the probability for an outbreak to receive a CATI response, we limited our analysis on outbreaks that were responded to. We subsequently did not observe any consistent residual difference of initial severity between classes of response promptness and response intensity. However, our models were adjusted for potential confounders and took into account the heterogeneity between localities. This quasi-experimental study was also stratified on response promptness and on response intensity, which yielded consistent response effectiveness estimates (Shadish et al., 2002).
Analyses may also have been biased by missing epidemiological data. Indeed, some patients do not seek care, even when they experience severe dehydration. Besides, stool sampling for confirmation culture was not systematic, which certainly led us to overlook several authentic outbreaks and mis-select clusters of non-cholera diarrhoeas. It may have led us to misdate several outbreak onset and outbreak end. Depending on the differential distribution of these potential biases among classes of response promptness and intensity, these limits could have led to over- or under-estimation of the effectiveness of prompt and intense CATIs. Nevertheless, our outbreak definition aimed to deal with those missing data and be specific in order to analyse CATI effectiveness on definite outbreaks. Like for many diseases, no standardized cholera outbreak criteria exists, and several definitions may be more or less suitable depending on interventions and analyses objectives (Brady et al., 2015). Our retrospectively defined outbreaks may be an approximate unit of analysis in terms of space, time and population, which may also have biased effectiveness results. We therefore conducted a sensitivity analysis using alternative definitions, including systematically lab-confirmed cholera outbreaks, which showed consistent and robust estimates (Appendix 3.1). We also used mixed models in order to take into account heterogeneity between localities in the random effect (Berridge and Crouchley, 2011). Additional CATI effectiveness studies at the household and at the administrative commune levels are underway in Haiti.
Our study analysed 3887 CATIs prospectively notified by rapid response teams to UNICEF. But some additional CATIs may have been omitted, while other CATIs remained unrecorded because they were implemented by the EMIRA alone. Nevertheless, many of their respective CATIs actually overlapped, and we thus believe our response database to be reasonably exhaustive. Conversely, only 16% of complete CATIs were conducted in a locality experiencing a current outbreak. The remaining CATIs were implemented in response to sporadic cases that did not meet outbreak definition criteria, as illustrated by much higher rates with less stringent outbreak definitions (Appendix 3.1). Sporadic CATIs may have prevented, delayed or attenuated the emergence of outbreaks. They may also be associated with the propensity of future outbreak response. We thus included the frequency of previous complete CATIs in our analysis but found no significant association with response promptness or intensity.
Our study aimed to assess the overall effectiveness of a CATI strategy. It neither aimed to estimate the respective effectiveness of each response components, nor the optimal radius of intervention, which would warrant dedicated field studies comparing different types of interventions. We thus chose a conservative definition of complete CATIs and performed a sensitivity analysis with alternative CATI definitions that exhibited consistent results (Appendix 3.2). Because nearly all CATIs included house decontamination, education and chlorine distribution, stratified analyses on these activities were not possible. However, three effectiveness estimates out of four appeared higher when all CATIs included antibiotic prophylaxis than when no CATI did. Several trials have also suggested that chemoprophylaxis has a protective effect among household contacts of people with cholera (Reveiz et al., 2011), and a micro-simulation model suggested that administration of antibiotics in CATIs could effectively avert secondary cases (Finger et al., 2018). But considering the risk of resistance selection (Mhalu et al., 1979; Dromigny et al., 2002), the selected distribution of antibiotic prophylaxis to close contacts is usually not recommended (Global Task Force on Cholera Control (GTFCC), 2018) and must, at the minimum, be used with caution and close monitoring of antibiotic susceptibility. In Haiti, all clinical V. cholerae O1 isolates have remained susceptible to doxycycline between 2013 and 2019 (Haitian Ministry of Public Health and Population, MSPP). As suggested by previous field or modelling studies (Ali et al., 2016; Parker et al., 2017b; Parker et al., 2017a; Finger et al., 2018), adding the administration of a single-dose OCV during CATIs could be an effective, but likely logistically complex, strategy.
Overall, our results suggest that case-area targeted interventions are significantly effective to mitigate and shorten local cholera outbreaks. Household water treatment, sanitation and hygiene promotion, as well as antibiotic prophylaxis theoretically prevent both human-to-human and environment-to-human cholera transmission pathways. Regardless of their respective role, which has been much debated (Morris, 2011; Kupferschmidt, 2017; Rebaudet et al., 2019b), our results thus confirm the relevance of promoting rapid response teams as a key component of the new global strategy for cholera control (Global Task Force on Cholera Control (GTFCC), 2017; Anon, 2017). Such findings need to be replicated in other settings and at other spatial and time scales. It will be critical to understand where CATIs should be prioritized, which radius is optimal, and which intervention components are most effective.
Materials and methods
Study design, setting and cholera surveillance
Request a detailed protocolTo assess CATI effectiveness, we conducted a retrospective observational study, which compared the outcome of cholera outbreaks according to the promptness or intensity of response CATIs. This corresponded to a quasi-experimental study using a post test-only design with stratified groups (Shadish et al., 2002). The study was conducted from January 1, 2015, to December 31, 2017 in the Centre department, one of the 10 administrative districts of Haiti. Centre department covers an area of 3487 km2, with an altitude ranging from 69 m to 1959 m, and is administratively subdivided in 12 communes. In 2015, the Centre population was estimated to be 746,236 inhabitants, including 20% living in urban neighbourhoods, and 80% in numerous rural settlements (Institut Haitien de Statistique et d’Informatique (IHSI), 2015). For the purpose of this study, we designate urban neighbourhoods and rural settlements as ‘localities’.
In 2015–2017, 17 cholera treatment centres, cholera treatment units and acute diarrhoea treatment centres officially treated and recorded suspected cholera cases and associated deaths to the MSPP. A probable suspected cholera case was defined as a patient who develops acute watery diarrhoea with or without vomiting. Daily cases and deaths tolls aged <or ≥ five years old were separately notified to the department health directorates. From 2014, the health directorate of the Centre department established a line-listing of all suspected cholera cases, mentioning sex, age, date of admission, address and use of IV rehydration (a surrogate for severe dehydration). Finally, routine bacteriological confirmation of a subset of suspected cholera cases was performed at the National Laboratory of Public Health (LNSP) in Port-au-Prince Metropolitan Area, using stool sampling with Carry-Blair transport medium and standard culture and phenotyping methods (Centers for Disease Control and Prevention (CDC), 1999).
Procedures: rapid case-area targeted interventions (CATIs)
Request a detailed protocolFrom July-2013, the nationwide case-area targeted rapid response strategy to eliminate cholera in Haiti was laboriously but increasingly implemented throughout the country (Rebaudet et al., 2019a). In the Centre department between 2015 and 2017, UNICEF established a partnership with Zanmi Lasante, Oxfam, ACTED and IFRC (International Federation of Red Cross and Red Crescent), four NGOs that hired WaSH rapid response teams composed of local Haitian staff. MSPP also established its own teams called EMIRAs, which included healthcare workers (nurses, auxiliary nurses). Staff of the NGO rapid response teams and EMIRA worked together and deployed mixed teams, which were requested to respond to every suspected cholera case or death within 48 hr after admission at healthcare facility. For this purpose, rapid response teams were encouraged to get epidemiological cholera data on a daily basis from departmental health directorates and treatment centres (Rebaudet et al., 2019a). The core methodology of response CATIs had been established with the MSPP and its partners and included: (i) door-to-door visits to affected families and their neighbours (minimum five households depending on the local geography), who were proposed house decontamination by chlorine spraying of latrines and other potentially contaminated surfaces; (ii) on-site organization of education sessions about cholera and hygiene promotion; (iii) and distribution of one cholera kit per household (composed of five soaps, five sachets of ORS, and chlorine tablets [80 Aquatabs33 mg in urban settings or 150 Aquatabs in rural areas]). EMIRA staff also provided (iv) prophylactic antibiotics to contacts living in the same house as cholera cases with one dose of doxycycline 300 mg for non-pregnant adults only. When appropriate, rapid response teams also: (v) established manual bucket chlorination at drinking water collection points during one or more weeks, by hiring and instructing local volunteers; (vi) chlorinated water supply systems and reported potential malfunctions to DINEPA; (vii) supervised safe funeral practices for cholera casualties; and (viii) provided primary care to cholera cases found in the community. CATIs were prospectively documented and transmitted by WaSH rapid response teams to UNICEF with date, location (i.e., commune, communal section, locality) and implemented activities, including specific activities of embedded EMIRA staff.
Response CATIs were defined as complete if rapid response teams reported at least education, decontamination and distribution of chlorine tablets. A sensitivity analysis of CATI effectiveness estimates using alternative CATI definitions is provided in Appendix 3.2.
Outbreaks identification and characterization
Request a detailed protocolIn order to identify cholera outbreaks, we first cleaned the anonymised case line-listing provided by the health directorate of the Centre department, the anonymised stool culture database provided by the LNSP and the response database provided by UNICEF. We manually corrected date errors and duplicates. Using repeated field investigations, GPS coordinates provided by rapid response teams, and several geographic repositories (http://ihsi.ht/publication_cd_atlas.htm, https://www.indexmundi.com/zp/ha/, https://www.openstreetmap.org/, https://www.google.fr/maps, accessed Jul 1, 2019), we corrected case, culture and response addresses with unified and geolocated locality names. We included every suspected case, every stool culture positive for V. cholerae O1 and every complete CATI of a WaSH rapid response team reported in the Centre department between January 2015 and December 2017.
To assess response effectiveness, we needed to escape the double pitfall of an overly restrictive definition of outbreaks, for example by requiring a bacteriological documentation for each suspected case and, on the contrary, of an unspecific definition, in which a large number of non-cholera diarrhoea cases would have been included. In addition, we had to deal with the fact that some patients with a positive culture were missing in the line-listing. Considering the median and the maximum O1-serogroup cholera incubation period are about 1.5 and 7 days, respectively (Azman et al., 2013), we thus defined outbreaks by the occurrence of at least two suspected cholera cases with at least one severely dehydrated case or positive culture, within the same locality, during a three-day time window, and after a refractory period of at least 21 days with no case. Outbreak onset was defined as the date of the first suspected case or positive culture, and outbreak end as the date of the last case or positive culture before a refractory period of at least 21 days. We conducted a sensitivity analysis using alternative outbreak definitions (Appendix 3.1).
For each identified outbreak, we then counted the numbers of cases and positive culture during the first three days as surrogates of initial severity. With a median incubation period of 1.5 days (Azman et al., 2013), we considered that a response – even a prompt one – would have little impact on the occurrence of additional cases during the two days following detection of the first case. Using a geographic information system (GIS), we extracted locality characteristics such as median population density (Sorichetta et al., 2015) and travel time to the nearest town (Weiss et al., 2018), using 1000 m radius buffer zones. Because cholera transmission and CATI response against cholera in Haiti were found to be influenced by rainfall (Eisenberg et al., 2013; Rebaudet et al., 2019a), we obtained NASA satellite estimates of daily-accumulated rainfall (TRMM_3B42_daily v7, area-averaged with 0.25° x 0.25° accuracy) (https://giovanni.gsfc.nasa.gov/giovanni/, accessed Jul 1, 2019). We gathered vaccine coverage of OCV campaigns conducted between 2012 and 2014, and accumulated incidence rates of suspected cholera cases between 2010 and 2014 (Haitian Ministry of Public Health and Population, MSPP), as surrogates of the population immunity against cholera. In order to better take into account the propensity of localities to experience outbreaks and receive response CATIs, we also counted the number of previous cases per year and the number of previous complete CATIs per year in the same locality since the beginning of the study. To take into account the possible variation of CATI implementation and effectiveness over time, we divided the three-year study period into six semesters (first and last six months of every year).
We then considered that outbreaks were responded to if at least one complete CATI was implemented within seven days after the last recorded case of the outbreak. In order to characterize the response promptness in this subgroup of outbreaks, we first counted the number of days between outbreak onset and the first complete CATI, and split outbreaks that were responded to between four classes of response promptness: >7 days, 3 to 7 days, 2 days and ≤1 day. In order to characterize the response intensity in the subgroup of outbreaks that were responded to, we also counted the number of complete CATIs per outbreak, divided this number by the outbreak duration (in week), and split outbreaks that were responded to between four classes of response intensity: <0.25, 0.25 to 0.5, 0.5 to 1 and ≥1 CATIs per week. We also divided the number of complete CATIs per outbreak by the number of accumulated cases per outbreak, and split outbreaks that were responded to between four classes of response intensity: <0.25, 0.25 to 0.5, 0.5 to 1 and ≥1 CATIs per case. Finally, we calculated two surrogates of outbreak outcome: the number of accumulated suspected cases from the fourth day of outbreak (outbreak size), and the number of days between the first and the last reported case or culture (outbreak duration). A sensitivity analysis of CATI effectiveness using alternative response time windows and categories is provided in Appendix 3: Sensitivity analyses of CATI effectiveness.
Statistical analysis
Analysis of confounders
Request a detailed protocolThe assessment of a possible confounding by indication is detailed in Appendix 2 (Remschmidt et al., 2015). We found that CATI response was more likely in more severe outbreaks. To handle this major bias, we therefore assessed CATI effectiveness (CE) by analysing the outcome of outbreaks that were responded to, according to the response promptness and according to response intensity. In two separate analyses, we compared two endpoints between the four classes of response promptness and between the four classes of response intensity (exposure): the number of cases from the fourth day of outbreak (CE represented the proportion of averted cases); and the outbreak duration (CE represented the proportion of averted days).
As an initial univariate step, we looked for possible confounders among baseline outbreak characteristics of response groups. First, each possible confounder was modelled as an independent variable, and time to the first complete CATI as a recurrent time-to-event outcome, using Cox survival models for Andersen-Gill counting process (AG-CP). This AG-CP survival model was chosen to take into account the correlated repetitions of outbreaks within localities (Andersen and Gill, 1982). Each possible confounder was also modelled as a fixed effect variable, the number of CATIs per week or the number of CATIs per case as dependent variables, and localities as a random effect, using generalized linear mixed models (GLMMs) with a negative-binomial distribution. The mixed model approach aimed to take into account the homogeneous pattern within localities, and the negative-binomial distribution to take into account overdispersion (Berridge and Crouchley, 2011).
CATI effectiveness according to response promptness
Request a detailed protocolThe first evaluation of CATI effectiveness (CE1) was then performed by comparing the outbreak size (number of cases from the fourth day of outbreak) between the four classes of response promptness (time to the first complete CATI). For this, we used GLMMs with cases from the fourth day of outbreak as a dependent variable, localities as a random effect, and a negative-binomial distribution (Berridge and Crouchley, 2011). For each class of response promptness, we estimated the crude CATI effectiveness (cCE1) as: 1 – Incidence ratio. We then obtained adjusted estimates of CATI effectiveness (aCE1) by adjusting for confounders for which p-values were less than 0.25 at the initial univariate step (Mickey and Greenland, 1989). A sensitivity analysis of CATI effectiveness using alternative methods of covariate selection is provided in Appendix 3.3.
A second evaluation of CATI effectiveness (CE2) was performed by comparing the outbreak duration between the four classes of response promptness, using survival analyses censoring outbreak extinction. We assessed time-to-event by Kaplan-Meier analysis to illustrate the cumulative probability of outbreak end between the different response promptness classes. In order to estimate CATI effectiveness according to response promptness and take into account the correlated repetitions of outbreaks within localities, we then fitted Andersen-Gill (AG-CP) survival models (Andersen and Gill, 1982). For each class of response promptness, we estimated the crude CATI effectiveness (cCE2) as: 1 – (1/Hazard ratio). We then obtained adjusted estimates of CATI effectiveness (aCE2) using the same methodology.
CATI effectiveness according to response intensity
Request a detailed protocolWe then estimated the effectiveness of the response intensity, by comparing the outbreak size or duration between different classes of response intensity, using the same methodology as for the effectiveness according to response promptness. In order to avoid that cases or duration be included both within outcome and exposure variables, we approximated response intensity by the number of complete CATIs per week ratio when comparing the number of cases accumulated from the fourth day of outbreak (CE3). Conversely, we used the number of complete CATIs per case ratio when comparing the duration of outbreak (CE4).
For all effectiveness analyses, a p-value of less than 0.05 (two-sided) was considered to indicate statistical significance.
Effectiveness of antibiotic prophylaxis
Request a detailed protocolIn order to assess the effectiveness of antibiotic prophylaxis, we conducted similar comparisons of outbreak size using GLMMs or outbreak duration using Andersen-Gill (AG-CP) survival models according to the response promptness or to the response intensity, stratified by whether all complete CATIs or none of the complete CATIs included antibiotic prophylaxis. We adjusted estimates of CATI effectiveness for the same confounders as in previous analyses. Using a Bonferroni correction for multiple comparisons, a p-value of less than 0.025 (two-sided) was considered to indicate statistical significance.
Software
The GIS and the map were done using QGIS software v3.03 and layers obtained from Haiti Centre National de l’Information Géospatiale (CNIGS) (http://cnigs.ht/, accessed Jul 1, 2019). Analyses and graphs were done using RStudio version 1.0.136 for Mac with R version 3.4.2 and the {ggplot2}, {lme4}, {survival} and {survminer} packages.
Appendix 1
Data characteristics
The study was conducted from January 1, 2015, to December 31, 2017 in the Centre department, where the cholera started in October 2010 (Piarroux et al., 2011). Like the rest of the country, this administrative district experienced a massive epidemic in 2010–2011 (Gaudart et al., 2013; Figure 1—figure supplement 1). Incidence then gradually decreased from 2012 to 2014, but in the following years, the Centre department remained one of the most affected area of the country (Rebaudet et al., 2019a). In order to better guide the response CATI strategy, the health directorate of the Centre department established case line-listing from 2014. This gave us the opportunity to retrospectively assess the impact of these interventions at the level of localities between 2015 and 2017 (Figure 1—figure supplement 1).
We thus analysed suspected cholera cases, cholera stool cultures and case-area targeted interventions (CATIs) recorded between January 1, 2015, and December 31, 2017. Baseline characteristics of these data are summarized in Appendix 1—table 1.
Appendix 2
Assessment of confounding by indication
In the event that case-area targeted interventions (CATIs) were significantly more likely implemented in more severe outbreaks, estimates of CATI effectiveness could be underestimated by a confounding by indication (Remschmidt et al., 2015).
Therefore, we initially compared baseline characteristics of outbreaks that were responded to (≥1 complete CATI implemented within 7 days after the last recorded case of the outbreak) and outbreaks that were not. This outcome variable following a binomial distribution, univariate logistic mixed models were used to estimate odds ratios associated with each covariate. Localities were modelled as a random effect variable, in order to take into account the heterogeneity between localities.
We then evaluated CATI effectiveness (CE5) by comparing the number of cases from the 4th day of outbreak between outbreaks that were and were not responded to. For this, we used logistic mixed models with cases from the 4th day of outbreak as a dependent variable (binomial distribution) (Berridge and Crouchley, 2011). We estimated the crude CATI effectiveness (cCE5) as: 1 – Odds ratio. We then obtained adjusted estimates of CATI effectiveness (aCE5) by adjusting for the following covariates: the number of cases and the number of positive cultures during the first 3 days of outbreak, rainfall, population density, travel time to the nearest town, accumulated case incidence between 2010 and 2014, coverage of OCV campaigns between 2012 and 2014 OCV campaigns, and the number of semesters since the beginning of the study.
A second evaluation of CATI effectiveness (CE6) was performed by comparing the duration of outbreaks between outbreaks that were and were not responded to, using survival analyses censoring outbreak extinction. We fitted Cox models for Andersen-Gill counting process (AG-CP) (Andersen and Gill, 1982). We estimated the crude CATI effectiveness (cCE6) as: 1 – (1/Hazard ratio). We then obtained adjusted estimates of CATI effectiveness (aCE6) by adjusting for confounders for which p-values were less than 0.25 at the initial univariate step (Mickey and Greenland, 1989).
For all effectiveness analyses, a p-value of less than 0.05 (two-sided) was considered to indicate statistical significance.
Overall, 238 identified outbreaks (53%) received a complete CATI within 7 days after the last recorded case, while 214 did not (Figure 2, Appendix 2—table 1). The proportion of outbreaks that were responded to progressively increased along the study semesters. These outbreaks occurred in localities which were significantly more densely populated, were significantly closer to a town, had significantly been less targeted by a previous mass OCV campaign than localities of outbreaks were non-responded to (Appendix 2—table 1). Outbreaks that were responded to exhibited a more severe onset (numbers of suspected cholera cases and positive cultures during the first 3 days were significantly higher), and they tended to be preceded by more cases than outbreaks that received no response (Appendix 2—table 1). This indicated a significant confounding by indication. Outbreaks that were responded to appeared to be more frequently preceded by complete CATIs than outbreaks that received no response (Appendix 2—table 1), which suggested a higher propensity of response.
Outbreaks that were responded to exhibited a paradoxically worse outcome than outbreaks that received no response. The median number of cases from the 4th day was 1 (interquartile range [IQR], 5) and 0 (IQR, 2) in outbreaks that were and were not responded to, respectively (Appendix 2—figure 1A, Appendix 2—table 2). Whereas the median duration of outbreaks that received at least one complete CATI was 11 days (IQR, 26.8), it was 3 days (IQR, 8) for outbreaks that were not responded to (Appendix 2—figure 1B, Appendix 2—table 3). The distribution of the number of cases from the 4th day of outbreaks that received no response (Appendix 2—figure 1A) looked like the distribution of the number of cases from the 4th day of outbreaks that were responded to within ≤ 1 days (Figure 4A), and looked like the distribution of the number of cases from the 4th day of outbreaks that received ≥ 1 CATI per week (Figure 4C). The Kaplan-Meier curve of outbreaks that received no response (Appendix 2—figure 1B) also looked like the Kaplan-Meier curve of outbreaks that were responded to within ≤ 1 days (Figure 4B), and looked like the Kaplan-Meier curve of outbreaks that received ≥ 1 CATI per case (Figure 4D). This illustrates the effect of the confounding by indication: because outbreaks that were not responded to were initially less severe, their outcome appeared better than the outcome of outbreaks that were responded to; the outcome of outbreaks that were not responded to also appeared close to the outcome of outbreaks that received a prompt and intense response.
Consequently, the crude CATI effectiveness in reducing the number of cases from the 4th day of outbreak (cCE5) was estimated to be −228% (95% CI, −353 to −138), and after adjusting for potential confounders (aCE5), −411% (−638 to −254) (Appendix 2-table 2). The crude CATI effectiveness on the duration of outbreak (cCE6) was −319% (−457 to −216), and after adjusting for potential confounders (aCE6), −300% (−441 to −196) (Appendix 2—table 3). This confirmed the significant confounding by indication, which explained why outbreaks that were responded to paradoxically exhibited a worse outcome than outbreaks that received no response (Remschmidt et al., 2015). This may be explained by the fact that numerous little outbreaks ended automatically, often before mobile teams arrived for the response. In absence of randomization, response teams likely tended to give priority to initially more severe outbreaks.
Appendix 3
Sensitivity analyses of CATI effectiveness
To assess the potential impact of our choices of definitions for outbreaks and for complete CATIs, we conducted several sensitivity analyses of CATI effectiveness according to different outbreak definitions, response definitions, and methods of model adjustment.
Appendix 3.1
Alternative outbreak definitions
In the main analyses, a cholera outbreak was defined by the occurrence, within the same locality, of at least two suspected cholera cases with at least one severely dehydrated case or positive culture, during a three-day time window, and after a refractory period of at least 21 days with no case. Outbreak onset was defined as the date of the first suspected case or positive culture, and outbreak end as the date of the last case or positive culture before a refractory period of at least 21 days. The three-day time window was initially chosen as it roughly corresponds to twice the median time incubation period of cholera (Azman et al., 2013). The 21 day refractory period was initially chosen as it roughly corresponds to twice the maximum incubation period (Azman et al., 2013) after the end of symptoms in the last case.
Nevertheless, using various thresholds of suspected cases, severely dehydrated cases, positive cultures, time window, and refractory period, we thus tested several alternative outbreak definitions summarized in Appendix 3—table 1.
In order to estimate the confounding by indication, we then compared outbreak outcome between outbreaks that were and were not responded to, as described in Appendix 2. In order to estimate CATI effectiveness, we compared outbreak outcome between classes of response promptness and between classes of response intensity, as described in the main manuscript.
Overall, all outbreak definitions led to a significant confounding by indication (Appendix 3—table 2). When considering only outbreaks that were responded to, CATI effectiveness according to response promptness and response intensity on the reduction of accumulated cases and on the reduction of outbreak duration remained consistent, irrespective of the adopted outbreak definition (Appendix 3—table 2). Some alternative outbreak definitions even brought higher effectiveness estimates than the definition used in the main manuscript.
Appendix 3.2
Alternative response definitions
In the main analyses, cholera outbreaks were considered as responded to if they received at least one complete CATI (i.e. a case-area targeted intervention with at least education, house decontamination by spraying and distribution of chlorine tablets) within 7 days after the last recorded case of the outbreak.
We thus tested several alternative response definitions summarized in Appendix 3—table 3.
In order to estimate the confounding by indication, we then compared outbreak outcome between outbreaks that were and were not responded to, as described in Appendix 2. In order to estimate CATI effectiveness, we compared outbreak outcome between classes of response promptness and between classes of response intensity, as described in the main manuscript.
Overall, all response definitions led to a significant confounding by indication (Appendix 3—table 4). When considering only outbreaks that were responded to, CATI effectiveness according to response promptness and response intensity on the reduction of accumulated cases and on the reduction of outbreak duration remained consistent, irrespective of the adopted response definition (Appendix 3—table 4). Some alternative response definitions even brought higher effectiveness estimates than the definition used in the main manuscript.
Appendix 3.3
Alternative adjustment methods for effectiveness estimates
In the main analyses, we adjusted CATI effectiveness estimates for covariates for which p-values were less than 0.25 at the initial univariate step (Mickey and Greenland, 1989).
We thus tested two alternative methods of confounder selection: adjusting on all eight covariates (number of cases and number of positive cultures during the first 3 days of outbreak, rainfall, population density, travel time to the nearest town, accumulated case incidence between 2010 and 2014, coverage of OCV campaigns between 2012 and 2014 OCV campaigns, and number of semesters since the beginning of the study); and minimizing the Akaike information criterion (AIC) of the models in order to avoid overfitting (Appendix 3—table 5).
Overall, all adjustment methods of models led to consistent CATI effectiveness estimates (Appendix 3—table 5). Inclusion of all covariates did not bring important overfitting.
Data availability
Data generated or analysed during this study are included in the manuscript and supporting files. Source data files have been provided for Figures 1 and 4.
References
-
Cox's Regression Model for Counting Processes: A Large Sample StudyThe Annals of Statistics 10:1100–1120.https://doi.org/10.1214/aos/1176345976
-
The incubation period of cholera: a systematic reviewJournal of Infection 66:432–438.https://doi.org/10.1016/j.jinf.2012.11.013
-
Micro-Hotspots of risk in urban cholera epidemicsThe Journal of Infectious Diseases 218:1164–1168.https://doi.org/10.1093/infdis/jiy283
-
ReportLaboratory Methods for the Diagnosis of Epidemic Dysentery and CholeraNational Center for Infectious Diseases.
-
Cholera cases cluster in time and space in matlab, Bangladesh: implications for targeted preventive interventionsInternational Journal of Epidemiology 38:2134–2139.https://doi.org/10.1093/ije/dyw267
-
Cholera outbreak in a village in south India - Timely action saved livesJournal of Infection and Public Health 6:35–40.https://doi.org/10.1016/j.jiph.2012.05.003
-
Emergence and rapid spread of tetracycline-resistant Vibrio cholerae strains, MadagascarEmerging Infectious Diseases 8:336–338.https://doi.org/10.3201/eid0803.010258
-
Spatio-temporal dynamics of cholera during the first year of the epidemic in HaitiPLOS Neglected Tropical Diseases 7:e2145.https://doi.org/10.1371/journal.pntd.0002145
-
Randomized controlled trial of Hospital-Based hygiene and water treatment intervention (CHoBI7) to reduce choleraEmerging Infectious Diseases 22:233–241.https://doi.org/10.3201/eid2202.151175
-
ReportStantard Operating Procedures. Responding to a Poliovirus Event or Outbreak. Part 1: General SOPsGeneva, Switzerland: World Health Organization (WHO).
-
ReportPopulation Totale, Population De 18 Ans Et Plus, Ménages Et Densités Estimés En 2015Direction Des Statistiques Démographiques Et Sociales, Ministère De l’Economie Et Des Finances, République d’Haïti.
-
Eliminating cholera transmission in HaitiNew England Journal of Medicine 376:101–103.https://doi.org/10.1056/NEJMp1614104
-
The impact of confounder selection criteria on effect estimationAmerican Journal of Epidemiology 129:125–137.https://doi.org/10.1093/oxfordjournals.aje.a115101
-
ReportNational Plan for the Elimination of Cholera in Haiti 2013-2022Port-au-Prince, Haiti: Ministry of Public Health and Population.
-
Cholera—Modern Pandemic Disease of Ancient LineageEmerging Infectious Diseases 17:2099–2104.https://doi.org/10.3201/eid1711.111109
-
Timely response and containment of 2016 cholera outbreak in northern ZambiaMedical Journal of Zambia 43:64–69.
-
Understanding the cholera epidemic, HaitiEmerging Infectious Diseases 17:1161–1168.https://doi.org/10.3201/eid1707.110059
-
Encyclopedia of Infectious Diseases : Modern Methodologies639–653, Needs for an integrative approach of epidemics: The Example of Cholera, Encyclopedia of Infectious Diseases : Modern Methodologies, Wiley & Sons, Inc.
-
Cholera in Haiti: successes and challenges in a low income countryWkly Epidemiol Rec 90:542–544.
-
The case-area targeted rapid response strategy to control cholera in Haiti: a four-year implementation studyPLOS Neglected Tropical Diseases 13:e0007263.https://doi.org/10.1371/journal.pntd.0007263
-
ReportHaïti. Enquête Mortalité, Morbidité Et Utilisation Des Services (EMMUS-V)Pétionville, Haïti: Institut Haïtien de l’Enfance.
-
ReportPlan National D’élimination Du Choléra. Développement Du Moyen Terme. Juillet 2016 - Décembre 2018République d’Haïti.
-
Experimental and Quasi-Experimental Designsdesigns for Generalizedgeneralized Causalcausal Inferenceinference103–134, Quasi-Experimental designs that either lack a control group or lack pretest observations on the outcome, Experimental and Quasi-Experimental Designsdesigns for Generalizedgeneralized Causalcausal Inferenceinference, Houghton Mifflin Company.
-
Clinical outcomes in household contacts of patients with cholera in BangladeshClinical Infectious Diseases 49:1473–1479.https://doi.org/10.1086/644779
-
Bacterial shedding in household contacts of cholera patients in Dhaka, BangladeshThe American Journal of Tropical Medicine and Hygiene 91:738–742.https://doi.org/10.4269/ajtmh.14-0095
-
ReportEbola and Marburg Virus Disease Epidemics: Preparedness, Alert, Control, and EvaluationGeneva, Switzerland: World Health Organization.
Article and author information
Author details
Funding
UNICEF
- Samuel Beaulieu
- Gregory Bulit
- Stanislas Rebaudet
Assistance Publique - Hôpitaux de Marseille
- Stanislas Rebaudet
Assistance Publique - Hôpitaux de Paris
- Renaud Piarroux
The funders of this study (UNICEF, APHM, APHP) had staff (co-authors of this manuscript) who had a role in data collection, analyses and writing of the report. However, the funders had no role in study design, data collection and analysis, decision to publish, or preparation of the manuscript.
Acknowledgements
We are grateful to the staff of MSPP, UNICEF, DINEPA and NGOs, who cared for patients, implemented and coordinated field responses, analysed stool cultures, gathered epidemiological and intervention data. We thank Bevan Hurley for editing the manuscript.
Ethics
Human subjects: All analyses retrospectively included routinely collected cholera surveillance and control data. The study protocol received authorization #1718-24 from the National Bioethics Committee of Haiti Ministry of Health (MSPP). The study only analysed anonymised data. Informed consent from patients and from people who benefited from a response intervention was therefore not required for this study.
Copyright
© 2019, Michel et al.
This article is distributed under the terms of the Creative Commons Attribution License, which permits unrestricted use and redistribution provided that the original author and source are credited.
Metrics
-
- 1,704
- views
-
- 209
- downloads
-
- 24
- citations
Views, downloads and citations are aggregated across all versions of this paper published by eLife.
Download links
Downloads (link to download the article as PDF)
Open citations (links to open the citations from this article in various online reference manager services)
Cite this article (links to download the citations from this article in formats compatible with various reference manager tools)
Further reading
-
- Epidemiology and Global Health
Artificially sweetened beverages containing noncaloric monosaccharides were suggested as healthier alternatives to sugar-sweetened beverages. Nevertheless, the potential detrimental effects of these noncaloric monosaccharides on blood vessel function remain inadequately understood. We have established a zebrafish model that exhibits significant excessive angiogenesis induced by high glucose, resembling the hyperangiogenic characteristics observed in proliferative diabetic retinopathy (PDR). Utilizing this model, we observed that glucose and noncaloric monosaccharides could induce excessive formation of blood vessels, especially intersegmental vessels (ISVs). The excessively branched vessels were observed to be formed by ectopic activation of quiescent endothelial cells (ECs) into tip cells. Single-cell transcriptomic sequencing analysis of the ECs in the embryos exposed to high glucose revealed an augmented ratio of capillary ECs, proliferating ECs, and a series of upregulated proangiogenic genes. Further analysis and experiments validated that reduced foxo1a mediated the excessive angiogenesis induced by monosaccharides via upregulating the expression of marcksl1a. This study has provided new evidence showing the negative effects of noncaloric monosaccharides on the vascular system and the underlying mechanisms.
-
- Epidemiology and Global Health
- Microbiology and Infectious Disease
Influenza viruses continually evolve new antigenic variants, through mutations in epitopes of their major surface proteins, hemagglutinin (HA) and neuraminidase (NA). Antigenic drift potentiates the reinfection of previously infected individuals, but the contribution of this process to variability in annual epidemics is not well understood. Here, we link influenza A(H3N2) virus evolution to regional epidemic dynamics in the United States during 1997—2019. We integrate phenotypic measures of HA antigenic drift and sequence-based measures of HA and NA fitness to infer antigenic and genetic distances between viruses circulating in successive seasons. We estimate the magnitude, severity, timing, transmission rate, age-specific patterns, and subtype dominance of each regional outbreak and find that genetic distance based on broad sets of epitope sites is the strongest evolutionary predictor of A(H3N2) virus epidemiology. Increased HA and NA epitope distance between seasons correlates with larger, more intense epidemics, higher transmission, greater A(H3N2) subtype dominance, and a greater proportion of cases in adults relative to children, consistent with increased population susceptibility. Based on random forest models, A(H1N1) incidence impacts A(H3N2) epidemics to a greater extent than viral evolution, suggesting that subtype interference is a major driver of influenza A virus infection ynamics, presumably via heterosubtypic cross-immunity.