Advertisement
Research Article| Volume 66, P90-95, January 2018

Download started.

Ok

Use of oral cholera vaccine as a vaccine probe to define the geographical dimensions of person-to-person transmission of cholera

Open AccessPublished:November 21, 2017DOI:https://doi.org/10.1016/j.ijid.2017.11.020

      Highlights

      • Cholera, a waterborne disease, is known to be transmitted from person to person.
      • However, the geographical dimensions of this transmission are uncertain.
      • Oral cholera vaccines were used as probes to define these dimensions.
      • Dimensions were 500 m in rural Bangladesh and 150 m in Kolkata, India.
      • Person-to-person transmission of cholera can be sustained over long distances.

      Abstract

      Background

      Cholera is known to be transmitted from person to person, and inactivated oral cholera vaccines (OCVs) have been shown to confer herd protection via interruption of this transmission. However, the geographic dimensions of chains of person-to-person transmission of cholera are uncertain. The ability of OCVs to confer herd protection was used to define these dimensions in two cholera-endemic settings, one in rural Bangladesh and the other in urban India.

      Methods

      Two large randomized, placebo-controlled trials of inactivated OCVs, one in rural Matlab, Bangladesh and the other in urban Kolkata, India, were reanalyzed. Vaccine herd protection was evaluated by relating the risk of cholera in placebo recipients to vaccine coverage of surrounding residents residing within concentric rings. In Matlab, concentric rings in 100-m increments up to 700 m were evaluated; in Kolkata, 50-m increments up to 350 m were evaluated.

      Results

      One hundred and eight cholera cases among 24 667 placebo recipients were detected during 1 year of post-vaccination follow-up at Matlab; 128 cholera cases among 34 968 placebo recipients were detected during 3 years of follow-up in Kolkata. Consistent inverse relationships were observed between vaccine coverage of the ring and the risk of cholera in the central placebo recipient for rings with radii up to 500 m in Matlab and up to 150 m in Kolkata.

      Conclusions

      These results suggest that the dimensions of chains of person-to-person transmission in endemic settings can be quite large and may differ substantially from setting to setting. Using OCVs as ‘probes’ to define these dimensions can inform geographical targeting strategies for the deployment of these vaccines in endemic settings.

      Keywords

      Introduction

      Killed oral cholera vaccines (OCVs) are now stockpiled by the World Health Organization (WHO) and are recommended public health tools for the control of cholera outbreaks (
      • Martin S.
      • Costa A.
      • Perea W.
      Stockpiling oral cholera vaccine.
      ). It is generally accepted that the rational use of such vaccines, either for epidemic or endemic cholera, will require that vaccination be targeted to geographically circumscribed populations at greatest risk of cholera. Knowledge of the geographic dimensions of chains of person-to-person cholera transmission will be important for effective geographic targeting, as OCVs have been demonstrated to confer both direct protection to vaccinees and herd protection to populations, the latter operating via interruption of person-to-person transmission.
      The geographic dimensions of chains of person-to-person cholera transmission were estimated in this study. A geographic information system (GIS)-based method that has been employed previously to evaluate whether OCVs confer vaccine herd protection was used (
      • Ali M.
      • Emch M.
      • von Seidlein L.
      • Yunus M.
      • Sack D.A.
      • Rao M.
      • et al.
      Herd immunity conferred by killed oral cholera vaccines in Bangladesh: a reanalysis.
      ). It was reasoned that because vaccine herd protection results from interruption of person-to-person cholera transmission, delineation of the dimensions of OCV herd protection should demarcate the dimensions of chains of person-to-person transmission.

      Methods

      Overview

      The vaccine herd protective effects of inactivated OCVs were analyzed in two randomized, placebo-controlled trials, one in Matlab, Bangladesh (
      • Clemens J.
      • Sack D.A.
      • Harris J.
      • van Loon F.
      • Chakraborty J.
      • Ahmed F.
      • et al.
      Field trial of oral cholera vaccines in Bangladesh: results from long-term follow-up.
      ) and the other in Kolkata, India (
      • Sur D.
      • Lopez A.L.
      • Kanungo S.
      • Paisley A.
      • Manna B.
      • Niyogi S.K.
      • et al.
      Efficacy and safety of a modified killed whole-cell oral cholera vaccine in India: an interim analysis of a cluster-randomized, double-blind, placebo-controlled trial.
      ), using a GIS-based method (
      • Ali M.
      • Emch M.
      • von Seidlein L.
      • Yunus M.
      • Sack D.A.
      • Rao M.
      • et al.
      Herd immunity conferred by killed oral cholera vaccines in Bangladesh: a reanalysis.
      ,
      • Ali M.
      • Sur D.
      • You Y.A.
      • 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.
      ). In this approach, a ‘virtual cluster’ is defined as persons whose residences are within a specified radius of the residence of each person under analysis (termed the ‘focal person’), and the risk of cholera for each focal person under analysis is related to the vaccine coverage in surrounding virtual clusters. An inverse relationship suggests vaccine herd protection, and when the focal person under analysis has received placebo, indirect vaccine protection is measured (
      • Ali M.
      • Emch M.
      • von Seidlein L.
      • Yunus M.
      • Sack D.A.
      • Rao M.
      • et al.
      Herd immunity conferred by killed oral cholera vaccines in Bangladesh: a reanalysis.
      ,
      • Ali M.
      • Sur D.
      • You Y.A.
      • 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.
      ,
      • Clemens J.
      • Shin S.
      • Ali M.
      New approaches to the assessment of vaccine herd protection in clinical trials.
      ). Conceptually, the maximum radius in which indirect vaccine protection is demonstrated should demarcate the geographical size of a surrounding population that puts an unimmunized individual at risk of becoming infected via person-to-person transmission. In this study, the indirect OCV protection of each focal person (placebo recipient) by vaccination in successive rings of persons in surrounding residences was analyzed, thus identifying the dimensions of chains of person-to-person cholera transmission in the two study settings.

      OCV trials under analysis

      The two randomized, placebo-controlled trials analyzed in this study were conducted in rural Matlab, Bangladesh, an area bisected by the Dhonagoda River, and in urban Kolkata, India, comprising wards 29, 30, and 33 (Figure 1); these trials have been described in detail elsewhere (
      • Clemens J.
      • Sack D.A.
      • Harris J.
      • van Loon F.
      • Chakraborty J.
      • Ahmed F.
      • et al.
      Field trial of oral cholera vaccines in Bangladesh: results from long-term follow-up.
      ,
      • Sur D.
      • Lopez A.L.
      • Kanungo S.
      • Paisley A.
      • Manna B.
      • Niyogi S.K.
      • et al.
      Efficacy and safety of a modified killed whole-cell oral cholera vaccine in India: an interim analysis of a cluster-randomized, double-blind, placebo-controlled trial.
      ,
      • Clemens J.
      • Sack D.
      • Harris J.
      • Chakraborty J.
      • Khan M.R.
      • Stanton B.
      • et al.
      Field trial of oral cholera vaccines in Bangladesh.
      ). Dosing with inactivated OCV or placebo was conducted in 1985 in Matlab and in 2006 in Kolkata. In Matlab, children aged 2–14 years and non-pregnant female adults (≥15 years) were eligible to participate in the trial. In Kolkata, non-pregnant persons aged ≥1 year were eligible. In Matlab, eligible persons were individually randomized to a three-dose regimen of an oral cholera toxin B subunit-killed whole cell (BS-WC) vaccine, oral killed whole cell (WC)-only vaccine, or oral placebo. In Kolkata, eligible persons were randomized by residential dwelling to a two-dose regimen of an oral killed WC-only vaccine or oral placebo.

      Diarrhea surveillance

      In Matlab, surveillance was conducted for all diarrheal patients from the study area who attended either the International Centre for Diarrhoeal Disease Research, Bangladesh (icddr,b) Matlab hospital or two community-operated treatment centers. Diarrhea was defined as the presence of at least three loose or liquid motions in the 24 h before presentation, or one to two or an indeterminate number of loose or liquid stools in the 24 h before presentation with at least two signs of dehydration (poor skin turgor, sunken eyes, dry mucous membranes, weakened radial pulse) on presentation. Stools or rectal swabs were collected from these patients and were tested for Vibrio cholerae O1 using conventional microbiological methods (
      • Bopp C.
      • Ries A.
      • Wells J.
      Laboratory methods for the diagnosis of epidemic dysentery and cholera.
      ).
      In Kolkata, diarrhea surveillance was conducted in nine project health clinics and two governmental hospitals. Diarrhea was defined as having three or more loose stools in the 24 h before presentation, or one to two or an indeterminate number of loose or liquid stools in the 24 h before presentation together with moderate or severe dehydration, according to WHO criteria, on presentation (
      • WHO
      The treatment of diarrhea: a manual for physicians and other senior health workers.
      ). Rectal swabs were collected from all diarrheal patients and tested in the project laboratory for V. cholerae O1 and V. cholerae O139 using conventional methods (
      • Bopp C.
      • Ries A.
      • Wells J.
      Laboratory methods for the diagnosis of epidemic dysentery and cholera.
      ).
      In both studies, cholera was defined as non-bloody diarrhea in which V. cholerae O1 was isolated. V. cholerae O139 was not isolated in the Kolkata trial. In the present analyses, cholera cases occurring during 1 year after the time of vaccination in Matlab and during 3 years after the time of vaccination in Kolkata were considered; these were both intervals in which indirect herd OCV protection was demonstrated in the earlier analyses (
      • Ali M.
      • Emch M.
      • von Seidlein L.
      • Yunus M.
      • Sack D.A.
      • Rao M.
      • et al.
      Herd immunity conferred by killed oral cholera vaccines in Bangladesh: a reanalysis.
      ,
      • Ali M.
      • Sur D.
      • You Y.A.
      • 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.
      ).

      Defining virtual clusters

      In the earlier published analyses of these trials (
      • Ali M.
      • Emch M.
      • von Seidlein L.
      • Yunus M.
      • Sack D.A.
      • Rao M.
      • et al.
      Herd immunity conferred by killed oral cholera vaccines in Bangladesh: a reanalysis.
      ,
      • Ali M.
      • Sur D.
      • You Y.A.
      • 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.
      ), a statistical criterion was used to define the radius around each placebo recipient to demarcate the virtual cluster of people for determination of vaccine coverage. This criterion was related to the variability of the variances of vaccine coverage across virtual clusters, and was not based on whether vaccine coverage was inversely related to the risk of cholera. This criterion yielded a radius of 500 m for Matlab and 250 m for Kolkata.
      In the present analyses, virtual clusters were defined a priori as concentric rings of residential populations in 100-m increments up to 700 m (seven rings centering each placebo recipient) in Matlab, and in 50-m increments up to 350 m (also seven rings centering each placebo recipient) in Kolkata. The difference in the size of increments was related to the dispersed population in Matlab versus the congested urban setting in Kolkata. Vaccine coverage in each ring was computed as the number of recipients of at least two doses of vaccine divided by the number of all eligible residents at the time of the first dose.

      Statistical analysis

      The earlier analyses found no design effect due to the residence-based allocation in Kolkata; therefore this trial was analyzed as if it had been individually randomized (
      • Sur D.
      • Lopez A.L.
      • Kanungo S.
      • Paisley A.
      • Manna B.
      • Niyogi S.K.
      • et al.
      Efficacy and safety of a modified killed whole-cell oral cholera vaccine in India: an interim analysis of a cluster-randomized, double-blind, placebo-controlled trial.
      ). Indirect vaccine protection was assessed by relating the risk of cholera in recipients of at least two doses of placebo (focal persons) to vaccine coverage of the defined rings around these recipients. In crude analyses the relationship for each ring between the quintile of vaccine coverage around each focal person (quintiles were ascertained from the distribution of coverages for all rings of the same radius) and the risk of cholera in the focal person, was evaluated, and this relationship was assessed with the Cochran–Armitage test for trend. In multivariable logistic regression models, the relationship between the levels of vaccine coverage of the defined ring (expressed dimensionally and fitted as an independent variable) and the occurrence of cholera for each focal person (the dependent variable) was assessed, controlling for potentially confounding variables. Based on the findings of the earlier work, candidate variables for inclusion as covariates in the models for Matlab were age at the date of first dose, distance to the icddr,b hospital in Matlab, and distance to the Dhonagoda River. Age at date of first dose, sex, owning at least one luxury item, per-capita household expenditure, washing hands with soap and water after defecation, and living in an owned house were considered as candidate covariates in the Kolkata models. To select the covariates for the final models, those candidate variables found to be associated with the risk of cholera at p< 0.10 in the bivariate models were included, and these were retained in the final model if they remained significant at p< 0.20 in a backward elimination process. Multivariable odds ratios (OR) relating each percentage increase in vaccine coverage to the occurrence of cholera were estimated by exponentiation of the coefficient for the vaccine coverage variable in the models; p-values and 95% confidence intervals (CI) for these OR were estimated using the standard errors of these coefficients. All p-values and 95% CI were two-sided.

      Ethics

      The project in Matlab was approved by the Ethical Review Committees of the icddr,b, Dhaka, Bangladesh and WHO (Geneva). The project in Kolkata was approved by the ethics committees of the National Institute of Cholera and Enteric Diseases (NICED) and the Health Ministry Screening Committee of India, and the International Vaccine Institute Institutional Review Board. Informed consent was obtained from all participants in these studies.

      Results

      There were 124 035 individuals eligible for vaccination in Matlab and 107 774 individuals eligible for vaccination in Kolkata. Among these populations, 108 cholera cases occurred among 24 667 recipients of at least two doses of placebo during the 1 year of follow-up in Matlab; 128 cholera cases occurred among 34 968 two-dose placebo recipients during the 3 years of follow-up in Kolkata.
      The previously published reports that evaluated vaccine coverage from 0 m to 500 m around placebo recipients in Matlab (
      • Ali M.
      • Emch M.
      • von Seidlein L.
      • Yunus M.
      • Sack D.A.
      • Rao M.
      • et al.
      Herd immunity conferred by killed oral cholera vaccines in Bangladesh: a reanalysis.
      ) and from 0 m to 250 m around placebo recipients in Kolkata (
      • Ali M.
      • Sur D.
      • You Y.A.
      • 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.
      ) each found significant inverse relationships between vaccine coverage and the risk of cholera in placebo recipients (Table 1).
      Table 1Relationships between inactivated oral cholera vaccine coverage and the risk of cholera in placebo recipients in earlier published analyses (
      • Ali M.
      • Emch M.
      • von Seidlein L.
      • Yunus M.
      • Sack D.A.
      • Rao M.
      • et al.
      Herd immunity conferred by killed oral cholera vaccines in Bangladesh: a reanalysis.
      ,
      • Ali M.
      • Sur D.
      • You Y.A.
      • 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.
      ).
      KolkataMatlab
      HR
      Hazard ratio for the relative reduction in the risk of cholera for each percentage increase in vaccine coverage. Adjusted for age, sex, religion, distance from household to nearest river, distance from household to nearest treatment center (kilometers), and experienced dysentery during follow-up.
      p-Value95% CIOR
      Odds ratio for the relative reduction in the rate of cholera for each percentage increase in vaccine coverage. Adjusted for age, individuals living in a larger cluster specified in the stratification, wards, monthly per-capita expenditure of the household, individuals living in a household always washing hands with soap and water after defecation, individuals living in their own house, individuals living in a household owning at least one luxury item, and distance from the household to the nearest health clinic.
      p-Value95% CI
      Cluster level cholera vaccine coverage0.97<0.00010.96–0.980.96<0.00010.94–0.98
      HR, hazard ratio; CI, confidence interval; OR, odds ratio.
      a Hazard ratio for the relative reduction in the risk of cholera for each percentage increase in vaccine coverage. Adjusted for age, sex, religion, distance from household to nearest river, distance from household to nearest treatment center (kilometers), and experienced dysentery during follow-up.
      b Odds ratio for the relative reduction in the rate of cholera for each percentage increase in vaccine coverage. Adjusted for age, individuals living in a larger cluster specified in the stratification, wards, monthly per-capita expenditure of the household, individuals living in a household always washing hands with soap and water after defecation, individuals living in their own house, individuals living in a household owning at least one luxury item, and distance from the household to the nearest health clinic.
      Cholera vaccine coverage for each concentric ring in the present analysis is presented in Table 2. Average vaccine coverage in each ring ranged from 41% to 45% in Matlab and from 26% to 30% in Kolkata. In the crude analyses in Matlab, a consistent inverse relationship between the quintile of vaccine coverage within the concentric ring and the risk of cholera in the focal person at the center of the ring was observed for all rings with radii up to 500 m. However, such a relationship was observed only for rings with radii up to 200 m for focal persons in Kolkata (Table 3).
      Table 2Descriptive statistics for inactivated oral cholera vaccine coverage of populations residing in the different concentric rings surrounding placebo recipients in the Matlab (Bangladesh) and Kolkata (India) trials.
      Ring size (meters)Mean %Median %Minimum %Maximum %
      Matlab, Bangladesh
      0–10045.3546.60085.71
      101–20043.0145.080100.00
      201–30042.2644.850100.00
      301–40043.0444.870100.00
      401–50042.3244.740100.00
      501–60041.9844.32080.00
      601–70042.2544.42082.35
      Kolkata, India
      0–5026.5526.69071.47
      51–10029.7229.19062.62
      101–15030.0129.27068.92
      151–20030.7029.76060.67
      201–25031.0530.03074.42
      250–30030.8329.48080.00
      301–35030.0229.5512.1768.18
      Table 3Risk of cholera among placebo recipients by quintile of vaccine coverage for each concentric ring around the placebo recipients in the Matlab (Bangladesh) and Kolkata (India) trials.
      MatlabKolkata
      % coverage by distancePlacebo recipientsCasesRisk/1000% coverage by distancePlacebo recipientsCasesRisk/1000
      0–100 m0–50 m
      ≤36.365389305.57≤21.3810 633444.14
      36.37–44.235089254.9121.39–27.337641364.71
      44.24–50.595057244.7527.34–32.556613274.08
      50.60–57.224841193.9232.56–40.495711152.63
      ≥57.234291102.33≥40.50437061.37
      Z-statistic2.432.81
      (p-Value)
      The Z-statistics and p-values were derived using the Cochran–Armitage trend test, as described in the Methods section.
      (0.0158)(0.0049)
      101–200 m51–100 m
      ≤32.354911306.11≤23.137059385.38
      32.36–41.674968295.8423.14–27.187410354.72
      41.68–48.514919224.4727.19–32.036995223.15
      48.52–55.425001153.0032.04–37.896755213.11
      ≥55.434868122.47≥37.906749121.78
      Z-statistic3.403.85
      (p-Value)
      The Z-statistics and p-values were derived using the Cochran–Armitage trend test, as described in the Methods section.
      (0.0007)(0.0001)
      201–300 m101–150 m
      ≤30.864891346.95≤24.657410374.99
      30.87–41.034978244.8224.66–28.047430344.58
      41.04–47.984962153.0228.05–31.577116253.51
      47.99–55.004962255.0431.58–36.056453192.94
      ≥55.014874102.05≥36.066559131.98
      Z-statistic3.203.34
      (p-Value)
      The Z-statistics and p-values were derived using the Cochran–Armitage trend test, as described in the Methods section.
      (0.0014)(0.0008)
      301–400 m151–200 m
      ≤32.794882275.53≤24.876900405.80
      32.80–41.514922265.2824.88–28.327332243.27
      41.52–47.974975204.0228.33–31.717066192.69
      47.98–54.484945265.2631.72–36.346519223.37
      ≥54.49494391.82≥36.357151233.22
      Z-statistic2.502.18
      (p-Value)
      The Z-statistics and p-values were derived using the Cochran–Armitage trend test, as described in the Methods section.
      (0.0125)(0.0292)
      401–500 m201–250 m
      ≤32.774966295.84≤25.266284264.14
      32.78–41.714968275.4325.27–28.107039283.98
      41.72–47.344940255.0628.11–31.216225243.86
      47.35–53.374952173.4331.22–35.357493253.34
      ≥53.384841102.07≥35.367927253.15
      Z-statistic3.201.16
      (p-Value)
      The Z-statistics and p-values were derived using the Cochran–Armitage trend test, as described in the Methods section.
      (0.0014)(0.24)
      501–600 m251–300 m
      ≤32.394862285.76≤25.706547294.43
      32.40–41.304923142.8425.71–28.237293253.43
      41.31–46.694940183.6428.24–30.866899304.35
      46.70–52.434927295.8930.87–35.126833213.07
      ≥52.445015193.79≥35.137396233.11
      Z-statistic0.301.28
      (p-Value)
      The Z-statistics and p-values were derived using the Cochran–Armitage trend test, as described in the Methods section.
      (0.78)(0.20)
      601–700 m301–350 m
      ≤32.344827214.35≤25.437134233.22
      32.35–41.275019254.9825.44–28.296827223.22
      41.28–47.014985275.4228.30–30.886887294.21
      47.02–52.994886204.0930.89–34.196881304.36
      ≥53.004950153.03≥34.207239243.32
      Z-statistic1.20−0.56
      (p-Value)
      The Z-statistics and p-values were derived using the Cochran–Armitage trend test, as described in the Methods section.
      (0.24)(0.57)
      a The Z-statistics and p-values were derived using the Cochran–Armitage trend test, as described in the Methods section.
      The backward elimination process for selecting covariates for the final model resulted in no covariates retained for Matlab and only age retained for Kolkata. Multivariable logistic regression models showed a significant association between OCV coverage in the ring and the risk of cholera in focal persons for rings with radii up to 500 m in Matlab and up to 150 m for placebo recipients in Kolkata (Table 4). To assess whether the relationships between OCV coverage and cholera risk in focal persons, by distance from the focal person, could be explained simply by variations in population density, population density was added as a covariate in these models. Models that included population density as a covariate again showed significant inverse associations between OCV coverage and the risk of cholera in focal persons for distances up to 500 m in Matlab and 150 m in Kolkata.
      Table 4Odds ratios relating the risk of cholera in placebo recipients to inactivated oral cholera vaccine coverage in successive concentric rings of populations in the Matlab (Bangladesh) and Kolkata (India) trials.
      MatlabKolkata
      Ring size (meters)OR
      Adjusted for age at the date of first dose.
      95% CIp-ValueRing size (meters)OR
      Adjusted for age at the date of first dose.
      95% CIp-Value
      0–1000.9850.971–0.9990.03300–500.9840.970–0.9990.0315
      101–2000.9840.973–0.9960.007651–1000.9740.957–0.9930.0062
      201–3000.9840.973–0.9960.0070101–1500.9720.949–0.9960.0223
      301–4000.9850.973–0.9970.0182151–2000.9890.965–1.0130.3673
      401–5000.9860.974–0.9990.0283201–2501.0020.978–1.0280.8496
      501–6000.9990.985–1.0130.8695251–3000.9790.951–1.0070.1383
      601–7000.9950.981–1.0080.4343301–3501.0040.973–1.0350.8184
      OR, odds ratio; CI, confidence interval.
      a Adjusted for age at the date of first dose.

      Discussion

      Using analyses of the indirect herd protection conferred by OCVs as a probe to define the geographical dimensions of chains of person-to-person transmission of cholera in Matlab, a rural site, and Kolkata, an urban site, an inverse relationship was found between vaccine coverage of the surrounding population and the risk of cholera for placebo recipients for virtual clusters with radii up to 500 m in Matlab and up to 150 m in Kolkata. These results suggest that the dimensions of chains of person-to-person transmission in endemic settings can be quite large and may differ substantially from setting to setting.
      Cholera has long been known to occur in geographic clusters. Early studies of space–time clustering of cholera noted clustering of cases within the same household or within small groups of nearby households (
      • Craig M.
      Time-space clustering of Vibrio cholerae 01 in Matlab, Bangladesh, 1970–1982.
      ,
      • Glass R.I.
      • Lee J.V.
      • Huq M.I.
      • Hossain K.M.
      • Khan M.R.
      Phage types of Vibrio cholerae O1 biotype El Tor isolated from patients and family which vaccination contacts in Bangladesh: epidemiologic implications.
      ,
      • Sugimoto J.D.
      • Koepke A.A.
      • Kenah E.E.
      • Halloran M.E.
      • Chowdhury F.
      • Khan A.I.
      • et al.
      Household transmission of Vibrio cholerae in Bangladesh.
      ,
      • Giebultowicz S.
      • Ali M.
      • Yunus M.
      • Emch M.
      A comparison of spatial and social clustering of cholera in Matlab, Bangladesh.
      ,
      • Mosley W.H.
      • Alvero M.G.
      • Joseph P.R.
      • Tamayo J.F.
      • Gomez C.Z.
      • et al.
      Studies of cholera El Tor in the Philippines. 4. Transmission of infection among neighbourhood and community contacts of cholera patients.
      ,
      • Philippines Cholera Committee
      Study on the transmission of E1 Tor cholera during an outbreak in Can-Itom community in the Philippines.
      ,
      • Tamayo J.F.
      • Mosley W.H.
      • Alvero M.G.
      • Joseph P.R.
      • Gomez C.Z.
      • et al.
      Studies of cholera El Tor in the Philippines. 3. Transmission of infection among household contacts of cholera patients.
      ), leading to the conventional wisdom that the person-to-person spread of cholera is usually limited to household transmission. More recently, two analyses performed in Matlab and Kolkata found an increased risk of clinically significant cholera shortly after the detection of index cases for populations within 450 m of the residence of index cases in Matlab, and within 50 m of residences of index cases in Kolkata (
      • Ali M.
      • Debes A.K.
      • Luquero F.J.
      • Kim D.R.
      • Park J.Y.
      • Digilio L.
      • et al.
      Potential for controlling cholera using a ring vaccination strategy: re-analysis of data from a cluster-randomized clinical trial.
      ,
      • Debes A.K.
      • Ali M.
      • Azman A.S.
      • Yunus M.
      • Sack D.A.
      Hoclrea cases clusterin time and space in Matlab, Bngaladesh: implications for targeted preventive interventions.
      ).
      It is interesting to compare the present study estimates with those of these more recent estimates, particularly since the present study and the earlier studies were done at the same study sites (
      • Ali M.
      • Debes A.K.
      • Luquero F.J.
      • Kim D.R.
      • Park J.Y.
      • Digilio L.
      • et al.
      Potential for controlling cholera using a ring vaccination strategy: re-analysis of data from a cluster-randomized clinical trial.
      ,
      • Debes A.K.
      • Ali M.
      • Azman A.S.
      • Yunus M.
      • Sack D.A.
      Hoclrea cases clusterin time and space in Matlab, Bngaladesh: implications for targeted preventive interventions.
      ). All three analyses identified surrounding populations much larger than the immediate household or clusters of households as being relevant to cholera transmission, and all identified substantially larger zones in the rural Bangladesh site than in the urban India site. However, in contrast to the two earlier analyses in Matlab and Kolkata, which evaluated the population at risk for clinically significant cholera around identified cholera cases, regardless of the type of transmission, the present study analyses specifically addressed the size of surrounding zones in which the use of OCVs had a demonstrable effect on person-to-person transmission.
      Whether a simple proxy measure – population density – might serve in lieu of analyses of OCV herd protection to define the geographic dimensions of person-to-person transmission was examined. However, inclusion of population density as an additional covariate in the multivariable models failed to erase the relationships between vaccine coverage and cholera risk in focal persons, by distance from the focal person, indicating that population density cannot be used as a simple proxy measure for estimating the geographical dimensions of chains of person-to-person transmission.
      Several potential limitations of this analysis require discussion. First, the impact of vaccine coverage calculated at the time that vaccine doses were given was analyzed, rather than vaccine coverage during post-vaccination follow-up, which would have changed over time due to population births, deaths, and migrations. Such misclassification of coverage over time would most likely have been non-differential with respect to the risk of cholera and would have tended to make the estimates of the radii of persons contributing to person-to-person transmission of cholera conservative. Second, in the analyses the ORs were in the range of 0.974–0.999 for Matlab and 0.949–0.996 for Kolkata, which may seem modest in magnitude. However, it should be appreciated that these ORs reflect relative reductions of risk of cholera per one percentage increase in vaccine coverage. Third, for simplicity, the analyses considered concentric circular geographic rings of populations. It may be that zones of person-to-person transmission in the two settings did not conform geometrically to circles.
      Fourth, although randomized clinical trials were analyzed, the present analyses of herd protection were observational in design and were subject to the limitations of such a non-randomized design. Fifth, the maximum radius for herd protection may vary not only from site to site, but also from vaccine to vaccine; had the studies been performed with other cholera vaccines, the results may have been different, depending on the capacity of the vaccines studied to interrupt transmission. Nonetheless, the present analyses are of relevance to the deployment of inactivated OCVs, which are currently stockpiled by the WHO for deployment in cholera-affected settings. Sixth, the analyses addressed cholera occurring in endemic settings, in which the population has some level of natural immunity; whether the findings can be generalized to epidemic settings in immunologically naïve populations, where the global OCV stockpile is occasionally deployed, is uncertain and constitutes an important subject for future research. Finally, caution should be exercised in generalizing the different geographical dimensions of person-to-person transmission in Kolkata versus Dhaka to other urban versus rural areas. The Matlab and Kolkata sites differed in multiple ways, not all of which were related to their urban versus rural contexts.
      This study has both theoretical and practical implications. The maximum size of virtual clusters for which indirect vaccine protection was detected in the two trials differed from the sizes of virtual clusters that were analyzed in the earlier published analyses of the trials, in which a radius for virtual clusters was determined with an arbitrary statistical criterion (
      • Ali M.
      • Emch M.
      • von Seidlein L.
      • Yunus M.
      • Sack D.A.
      • Rao M.
      • et al.
      Herd immunity conferred by killed oral cholera vaccines in Bangladesh: a reanalysis.
      ,
      • Ali M.
      • Sur D.
      • You Y.A.
      • 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.
      ). This discrepancy suggests that work is required to improve the earlier method for sizing virtual clusters in analyzing vaccine herd protection. On the other hand, it was reassuring that the present study analyses showed a pattern of declining vaccine herd protection with increasing distance of the vaccinated population from the central placebo recipient under analysis, providing evidence of the construct validity of the method. From a practical perspective, the approach used herein may provide guidance to geographically targeted OCV strategies for the control of outbreaks of cholera occurring in endemic settings.
      In summary, the results of this study suggest that the dimensions of chains of person-to-person transmission in endemic settings can be quite large and may differ substantially from setting to setting. Using OCVs as ‘probes’ to define these dimensions can inform geographical targeting strategies for the deployment of these vaccines in endemic settings.

      Contributors

      MA and JC contributed to the study design. MA, DS, SK, SD, FM, SKB, and JC contributed to the implementation and supervision of the study. MA, BM, and DRK analyzed the data and took responsibility for the accuracy of the data analysis. All authors participated in the writing of the manuscript and approved the final version of the manuscript.

      Funding

      Bill & Melinda Gates Foundation. The icddr,b receives core funding from the governments of Bangladesh, Canada, Sweden, and the UK. The International Vaccine Institute receives core funding from the governments of South Korea and Sweden.

      Conflict of interest

      None declared.

      Acknowledgements

      We are grateful to the people of the study areas and the field staff for whom the project was successfully conducted. This study was funded by the Bill & Melinda Gates Foundation in grants to the icddr,b, the International Vaccine Institute (IVI), and the DOVE project of Johns Hopkins Bloomberg School of Public Health (OPP1148763). The study was also supported by core grants to the icddr,b from the governments of Bangladesh, Canada, Sweden, and the UK and to the IVI from the governments of Sweden and South Korea. The sponsors had no role in the study design, in the collection, analysis, and interpretation of the data, in the writing of the manuscript, or in the decision to submit the manuscript for publication.

      References

        • Ali M.
        • Emch M.
        • von Seidlein L.
        • Yunus M.
        • Sack D.A.
        • Rao M.
        • et al.
        Herd immunity conferred by killed oral cholera vaccines in Bangladesh: a reanalysis.
        Lancet. 2005; 366: 44-49
        • Ali M.
        • Sur D.
        • You Y.A.
        • 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: 1123-1131
        • Ali M.
        • Debes A.K.
        • Luquero F.J.
        • Kim D.R.
        • Park J.Y.
        • Digilio L.
        • et al.
        Potential for controlling cholera using a ring vaccination strategy: re-analysis of data from a cluster-randomized clinical trial.
        PLoS Med. 2016; 13e1002120
        • Bopp C.
        • Ries A.
        • Wells J.
        Laboratory methods for the diagnosis of epidemic dysentery and cholera.
        Centers for Disease Control and Prevention, Atlanta1999
        • Clemens J.
        • Sack D.
        • Harris J.
        • Chakraborty J.
        • Khan M.R.
        • Stanton B.
        • et al.
        Field trial of oral cholera vaccines in Bangladesh.
        Lancet. 1986; 2: 124-127
        • Clemens J.
        • Sack D.A.
        • Harris J.
        • van Loon F.
        • Chakraborty J.
        • Ahmed F.
        • et al.
        Field trial of oral cholera vaccines in Bangladesh: results from long-term follow-up.
        Lancet. 1990; 335: 270-273
        • Clemens J.
        • Shin S.
        • Ali M.
        New approaches to the assessment of vaccine herd protection in clinical trials.
        Lancet Infect Dis. 2011; 11: 482-487https://doi.org/10.1016/S1473-3099(10)70318-2
        • Craig M.
        Time-space clustering of Vibrio cholerae 01 in Matlab, Bangladesh, 1970–1982.
        Soc Sci Med. 1988; 26: 5-13
        • Debes A.K.
        • Ali M.
        • Azman A.S.
        • Yunus M.
        • Sack D.A.
        Hoclrea cases clusterin time and space in Matlab, Bngaladesh: implications for targeted preventive interventions.
        Int J Epidemiol. 2016; 45: 2134-2139
        • Giebultowicz S.
        • Ali M.
        • Yunus M.
        • Emch M.
        A comparison of spatial and social clustering of cholera in Matlab, Bangladesh.
        Health Place. 2011; 17: 490-497
        • Glass R.I.
        • Lee J.V.
        • Huq M.I.
        • Hossain K.M.
        • Khan M.R.
        Phage types of Vibrio cholerae O1 biotype El Tor isolated from patients and family which vaccination contacts in Bangladesh: epidemiologic implications.
        J Infect Dis. 1983; 148: 998-1004
        • Martin S.
        • Costa A.
        • Perea W.
        Stockpiling oral cholera vaccine.
        Bull World Health Organ. 2012; 90: 714
        • Mosley W.H.
        • Alvero M.G.
        • Joseph P.R.
        • Tamayo J.F.
        • Gomez C.Z.
        • et al.
        Studies of cholera El Tor in the Philippines. 4. Transmission of infection among neighbourhood and community contacts of cholera patients.
        Bull World Health Organ. 1965; 33: 651-660
        • Philippines Cholera Committee
        Study on the transmission of E1 Tor cholera during an outbreak in Can-Itom community in the Philippines.
        Bull World Health Organ. 1970; 43 ([No authors listed]): 413-419
        • Sugimoto J.D.
        • Koepke A.A.
        • Kenah E.E.
        • Halloran M.E.
        • Chowdhury F.
        • Khan A.I.
        • et al.
        Household transmission of Vibrio cholerae in Bangladesh.
        PLoS Negl Trop Dis. 2014; 8: e3314
        • Sur D.
        • Lopez A.L.
        • Kanungo S.
        • Paisley A.
        • Manna B.
        • Niyogi S.K.
        • et al.
        Efficacy and safety of a modified killed whole-cell oral cholera vaccine in India: an interim analysis of a cluster-randomized, double-blind, placebo-controlled trial.
        Lancet. 2009; 374: 1694-1702
        • Tamayo J.F.
        • Mosley W.H.
        • Alvero M.G.
        • Joseph P.R.
        • Gomez C.Z.
        • et al.
        Studies of cholera El Tor in the Philippines. 3. Transmission of infection among household contacts of cholera patients.
        Bull World Health Organ. 1965; 33: 645-649
        • WHO
        The treatment of diarrhea: a manual for physicians and other senior health workers.
        4th revision. World Health Organization, Geneva2005