Impact of the scale-up of piped water on urogenital schistosomiasis infection in rural South Africa

  1. Frank Tanser  Is a corresponding author
  2. Daniel K Azongo
  3. Alain Vandormael
  4. Till Bärnighausen
  5. Christopher Appleton
  1. Africa Health Research Institute, South Africa
  2. University of KwaZulu-Natal, South Africa
  3. University College London, United Kingdom
  4. Navrongo Health Research Centre, Ghana Health Service, Ghana
  5. Harvard T.H. Chan School of Public Health, United States
  6. University of Heidelberg, Germany

Abstract

Recent work has estimated that sub-Saharan Africa could lose US$3.5 billion of economic productivity every year as a result of schistosomiasis and soil-transmitted helminthiasis. One of the main interventions to control schistosomiasis is the provision of safe water to limit the contact with infected water bodies and break the cycle of transmission. To date, a rigorous quantification of the impact of safe water supplies on schistosomiasis is lacking. Using data from one of Africa’s largest population-based cohorts, we establish the impact of the scale-up of piped water in a typical rural South African population over a seven-year time horizon. High coverage of piped water in the community decreased a child’s risk of urogenital schistosomiasis infection eight-fold (adjusted odds ratio = 0.12, 95% CI 0.06–0.26, p<0.001). The provision of safe water could drive levels of urogenital schistosomiasis infection to low levels of endemicity in rural African settings.

https://doi.org/10.7554/eLife.33065.001

Introduction

Schistosomiasis is a chronic parasitic disease caused by blood flukes of the genus Schistosoma and transmitted by snails found in fresh water bodies that have been contaminated by Schistosoma eggs. Approximately, 240 million people worldwide (WHO, 2010; Stothard et al., 2009) are infected with schistosomiasis, with a disease burden of 3.3 million disability-adjusted life-years (Murray et al., 2012). In sub-Saharan Africa alone, an estimated 163 million people were infected in 2012 (Lai et al., 2015). Following the London declaration (Uniting to Combat Neglected Tropical Diseases, 2012), the World Health Organization (WHO) has set a goal of increasing mass drug administration coverage to 75% of at-risk children in endemic countries (WHO, 2013). Despite progress in many so-called ‘neglected tropical diseases’, the fourth progress report of the London Declaration records that schistosomiasis has the lowest coverage of all helminth diseases treatable by mass drug administration (Uniting to Combat Neglected Tropical Diseases, 2016).

For achieving sustained control, elimination, or eradication of schistosomiasis, improvements of water, sanitation, and hygiene infrastructure and modification of risk behaviour are necessary (Freeman et al., 2013; WHO, 2015; Steinmann et al., 2006; Echazú et al., 2015). One important approach is the provision of safe water supplies to reduce the need for contact with contaminated water bodies and diminish the risk of schistosomiasis transmission. However, quality evidence for the impact of provision of safe water on schistosomiasis infection is currently lacking. Recently, a systematic review concluded that previous studies on this topic were of poor quality and likely confounded by multiple socio-economic and other variables (Grimes et al., 2014). The authors conclude that new research is needed ‘on the relationships between water, sanitation, hygiene, human behaviour, and schistosome transmission'. In this regard, the previously used exposure measures of household or individual access to safe water (Grimes et al., 2014) will underestimate the true health impacts of the provision of safe water at a population level. This is because the introduction of safe water supplies into a community offers protection from schistosomiasis infection in two ways. Firstly, it offers direct protection to individuals with access to safe water by reducing their contact with infected water bodies through household domestic water collection activities. Secondly, the scale-up in access to safe water will confer indirect protection to members of a community in a manner that is analogous to the concept of ‘herd immunity’ (Fine et al., 2011). This form of protection is conferred through a reduction in the number of contacts that infected individuals have with open water bodies leading to a decrease in the overall levels of intensity of infection in the surrounding community. Thus, the greater the proportion of individuals in a community who are not actively contributing to the perpetuation of the schistosomiasis transmission cycle, the smaller the probability that those who interact with open water bodies will come into contact with free swimming cercariae released by infected snails.

Hence, to capture the true overall population-level impact of the introduction of piped water, it is necessary to move beyond a simple household/individual level exposure approach towards a sensitive community-level exposure estimate that quantifies the coverage of piped water in the surrounding local community over an extended period of time preceding the measurement of the disease outcome. To date, no studies have quantified the (likely non-linear) relationship between coverage of piped water (%) at the community level and corresponding reductions in risk of schistosomiasis infection. From an epidemic control and eradication perspective, quantification of such a relationship would provide the best indicator of the true overall impact of the scale-up in piped water and would reveal critical coverage targets for achieving meaningful reductions in risk of schistosomiasis infection.

To accomplish this goal, we nested a Schistosoma haematobium survey within one of Africa’s largest population-based cohorts in which all children and their families (N≈90,000) have been intensively followed up in an endemic area of rural KwaZulu-Natal, South Africa. We exploit the heterogeneity in the space-time scale-up of piped water over seven years and detailed longitudinal risk factor information to provide a rigorous quantification of the impact of the introduction of piped water on risk of urogenital schistosomiasis infection measured in children living in the study community. We then validate the veracity of our findings using an instrumental variable approach.

Results

Of the 2105 participants tested in the parasitological survey, 353 were infected with Schistosoma haematobium. The overall prevalence was 16.8% (95% CI: 15.2–18.4) and the prevalence of heavy infection (≥50 eggs/10 ml of urine) was 9.5% (95% CI: 8.3–10.8). Prevalence was twice as high in boys compared to girls (22.6% versus 11.4%) and peaked at 30% in boys at age 14 years and 15% in girls at age 12 years (Figure 1). Prevalence of heavy infection was 13.5% in boys (and peaked at 21.4% at age 14 years) and 5.7% in girls (peaking at 8.7% at age 12 years). 57% of all infected participants had heavy infections with a geometric mean egg burden of 65 eggs/10 ml urine. Overall, 80% of the egg burden was borne by just 3.4% of all pupils tested (20.4% of all infected participants).

Prevalence (95% CI) of Schistosoma haematobium infection by age and sex among children taking part in the parasitological survey (N = 2,105).

Darker colours represent heavy infections (≥50 eggs per 10 ml urine).

https://doi.org/10.7554/eLife.33065.002

Overall, 67% of participants who were present at school on the day of testing had parental consent and assent and provided a urine specimen. Only 2.2% of linked participants reported receiving praziquantel in the last 12 months (Table 1). Participation in the parasitological survey did not differ significantly by piped water access or any of the environmental and socio-economic control variables used in the analysis. However, there were differences in levels of participation by school (p<0.05) and small differences by age and sex (p<0.1). After standardising for these differences according to age, sex and school, the prevalence remained virtually unchanged at 16.9%.

Table 1
Characteristics of the primary school children within the survey who were linked to the population-based cohort (N=1976).
https://doi.org/10.7554/eLife.33065.003
CovariateTotalInfected (%)(95% CI)
Gender
 Female1016117 (11.5)(7.7–16.9)
 Male960217 (22.6)(17.5–28.7)
Age group
 911814 (11.9)(7.1–19.2)
 1036644 (12.0)(6.8–20.3)
 1159287 (14.7)(11.3–18.9)
 1246992 (19.6)(13.6–27.4)
 1322249 (22.1)(15.8–30.0)
 ≥1420948 (23.0)(15.1–33.3)
Community piped water (quintiles)*
 (Lowest) 1399102 (25.6)(13.0–44.1)
 241587 (21.0)(15.7–27.4)
 339765 (16.4)(10.5–24.6)
 438633 (8.5)(5.4–13.2)
 537947 (12.4)(8.9–16.9)
Household access to water
 No piped water30453 (17.4)(9.9–28.8)
 Piped water1672281 (16.8)(12.9–21.6)
Household assets quintiles
 (Poorest) 139883 (20.9)(13.3–31.1)
 238763 (16.3)(11.7–22.1)
 338963 (16.2)(11.7–22.0)
 436556 (15.3)(11.5–20.1)
 535957 (15.9)(10.4–23.6)
 Missing7812 (15.4)(9.5–23.9)
School grade
 Grade 51039186 (17.9)(12.2–25.5)
 Grade 6937148 (15.8)(10.4–23.2)
Praziquantel in the last 12 months
 No1933321 (16.6)(12.5–21.8)
 Yes4313 (30.2)(17.3–47.3)
Altitude (meters above sea level)
 <507629 (38.2)(19.3–61.4)
 50–100641168 (26.2)(19.2–34.7)
 100–150875108 (12.3)(9.5–16.0)
 150–20029622 (7.4)(4.4–12.4)
 >200887 (8.0)(3.5–17.2)
Distance to water body
 <1 km606112 (18.5)(14.2–23.7)
 1–2 km618122 (19.7)(14.0–27.0)
 2–3 km37668 (18.1)(10.8–28.7)
 >3 km37632 (8.5)(5.4–13.2)
Toilet in household
 No Toilet43864 (14.6)(9.9–21.1)
 Toilet1538270 (17.6)(13.3–22.8)
Land cover classification
 Closed Shrubland787184 (23.3)(17.6–30.3)
 Open Shrubland69683 (11.9)(8.8–15.9)
 Sparse Shrubland43752 (11.9)(8.7–16.0)
 Thickett5615 (26.8)(14.8–43.4)
Slope (quintiles)
 (Lowest) 139059 (15.1)(11.0–20.5)
 238672 (18.7)(12.9–26.3)
 340270 (17.4)(11.9–24.8)
 440470 (17.3)(12.4–23.7)
 539463 (16.0)(11.6–21.6)
  1. *Computes the proportion of households having access to piped water in the unique community surrounding each participant in the study (Figure 3). The Quintile (Q) ranges (min–max) are: Q1: 0–36; Q2: 37–59; Q3: 60–75; Q4: 76–92; Q5: 93–100.

The geospatial analysis using the Gaussian kernel approach revealed substantial spatial variation in prevalence of urogenital schistosomiasis across the different communities within the surveillance area (Figure 2). The kernel-derived prevalence of urogenital schistosomiasis varied between 0% and 85% and revealed two clear regions of high transmission where prevalence was >50%. A formal micro-geographical clustering analysis using the Kulldorff spatial scan statistic revealed four significant foci of infection around open water collection sites (relative risk = 2.99–5.63, p<0.05; Figure 2, Table 2). Two clusters of infection were located along the Nsane River in the northeast of the surveillance area (clusters 2 and 3) and the other two located around water bodies towards the south of the study area.

Geographical variations in Schistosoma haematobium prevalence across the surveillance area obtained by a Gaussian kernel applied to participants’ precise household locations.

Approximate locations of participants’ households are shown (incorporating an intentional random error) with white dots representing an infected child. Superimposed on the map are the clusters of infection independently identified by the Kulldorff spatial scan statistic (cluster 5, low relative-risk; cluster 1–4, high relative-risk). The National Road can be seen running along the Eastern boundary of the surveillance area towards Mozambique.

https://doi.org/10.7554/eLife.33065.004
Table 2
Significant spatial clusters (either unexpectedly high or low numbers) of Schistosoma haematobium infections identified by the Kulldorff spatial scan statistic (p<0.05) across the study area (see Figure 2)
https://doi.org/10.7554/eLife.33065.005
Cluster NumberRadius (km)Log-LikelihoodP-valuePrevalence (%)Relative Risk
10.9516.54<0.00184.65.63
22.6865.5<0.001504.66
31.1930.15<0.00149.24.53
41.9614.89 0.00241.92.99
52.9317.2<0.0010.760.05

Between 2001 and 2007, the provision of piped water in the surveillance area nearly doubled from 42% to 80% (Figure 3). The community-based mean coverage of piped water between 2001 and 2007 was 61.7% (Figure 3, bottom right). There is a high degree of temporal and spatial consistency in the patterns of expansion of piped water access across the five respective surveys, providing assurance of the quality of the underlying population-based data.

Time series of maps showing the coverage of piped water (%) between 2001 and 2007 in the study area (as measured using a Gaussian kernel approach) as well as mean piped water coverage over the full study period (bottom right).

Main roads are superimposed.

https://doi.org/10.7554/eLife.33065.006

We identified a clear dose-response between the risk of urogenital schistosomiasis infection and piped water coverage in the surrounding local community (Figure 4a). Participants living in communities with a high coverage of piped water (median coverage = 95.5%) were eight times less likely to be infected (aOR = 0.12, 95% CI 0.06–0.26, p<0.001; Model 1 of Table 3) in comparison to those living in communities with lowest coverage (median coverage = 8.4%). This reduction was greater for girls than boys (aOR = 0.07 versus aOR = 0.17 in boys; Model 1 of Table 4) after controlling for key socio-economic and environmental predictors of infection.

Figure 4 with 1 supplement see all
Adjusted odds ratio of Schistosoma haematobium infection (95% CI) by piped water coverage in the surrounding local community - coverage quintile (Left) and continuous piped water coverage (Right).

The piped water coverage measure (2001–2007) is derived using a Gaussian kernel to calculate the proportion of all households in the unique local community surrounding each participant having access to piped water (Figure 3). Odds ratios are adjusted for age, sex, household assets, toilet in household, landcover class, distance to water body, altitude, slope, treatment in the last 12 months and school grade. Standard errors are adjusted for clustering by school and grade.

https://doi.org/10.7554/eLife.33065.007
Table 3
Logistic regression analysis of the risk factors of Schistosoma haematobium infection.

Model 0 gives the univariate results and Model 1 includes all variables in the model. In Model 2, piped water coverage in the immediate community surrounding each participant has been substituted with household-level piped water covariate.

https://doi.org/10.7554/eLife.33065.009
Model 0: UnivariateModel 1: Community coverageModel 2: Household access
CovariateaOR(95% CI)P-valueaOR(95% CI)P-valueaOR(95% CI)P-value
Community piped water quintiles (vs Lowest)†
 20.77(0.34, 1.75)0.5290.39‡(0.23, 0.66)<0.001
 30.57(0.22, 1.50)0.2500.30(0.15, 0.59)<0.001
 40.27(0.10, 0.71)0.0090.16(0.08, 0.33)<0.001
 50.41(0.17, 0.99)0.0480.12(0.06, 0.26)<0.001
Household access to water (vs No)
 Yes0.96(0.56, 1.64)0.8700.54(0.33, 0.89)0.017
Gender (vs Female)
 Male2.24(1.64, 3.08)<0.0012.62(1.92, 3.59)<0.0012.41(1.77, 3.28)<0.001
Age testing
 Per unit1.19(1.07, 1.31)0.0011.21(1.08, 1.36)0.0011.18(1.06, 1.31)0.002
Grade (vs Grade 5)
 Grade 60.86(0.45, 1.66)0.6480.76(0.51, 1.11)0.1490.77(0.47, 1.28)0.310
Praziquantel in last 12 months (vs No)
 Yes2.18(1.05, 4.52)0.0381.27(0.60, 2.71)0.5291.48(0.69, 3.16)0.307
Altitude Class (vs < 50)
 50–1000.58(0.27, 1.22)0.1470.47(0.23, 0.96)0.0390.50(0.23, 1.09)0.081
 100–1500.23(0.09, 0.60)0.0030.20(0.09, 0.43)<0.0010.20(0.08, 0.51)0.001
 150–2000.13(0.04, 0.39)<0.0010.09(0.03, 0.25)<0.0010.11(0.04, 0.33)<0.001
 >2000.14(0.04, 0.51)0.0040.08(0.03, 0.29)<0.0010.12(0.03, 0.44)0.002
Landcover class (vs Sparse Shrubland)
 Closed Shrubland2.26(1.55, 3.28)<0.0011.56(1.05, 2.31)0.0302.41(1.63, 3.58)<0.001
 Open Shrubland1.00(0.70, 1.44)0.9891.03(0.72, 1.47)0.8631.44(0.96, 2.17)0.079
 Thickett2.71(1.28, 5.75)0.0101.75(0.82, 3.73)0.1452.52(1.23, 5.18)0.012
Slope (square root)
 per unit0.98(0.83, 1.16)0.8181.02(0.89, 1.16)0.7940.91(0.79, 1.04)0.159
Distance to water body (vs < 1 km)
 1–2 km1.08(0.74, 1.58)0.6670.78(0.56, 1.08)0.1310.99(0.69, 1.41)0.946
 2–3 km0.97(0.54, 1.75)0.9280.72(0.47, 1.12)0.1431.04(0.62, 1.77)0.874
 >3 km0.41(0.23, 0.74)0.0030.25(0.12, 0.49)<0.0010.44(0.24, 0.79)0.007
Toilet in household (vs No)
 Yes1.24(0.90, 1.72)0.1861.24(0.87, 1.76)0.2291.20(0.84, 1.72)0.319
Household assets quintile (vs Poorest)
 20.74(0.50, 1.09)0.1230.88(0.60, 1.27)0.4800.87(0.61, 1.25)0.459
 30.73(0.49, 1.09)0.1270.78(0.51, 1.18)0.2350.75(0.48, 1.16)0.186
 40.69(0.42, 1.14)0.1430.81(0.47, 1.40)0.4500.67(0.40, 1.11)0.118
 50.72(0.41, 1.24)0.2280.80(0.48, 1.34)0.3890.62(0.36, 1.05)0.075
 Missing0.69(0.39, 1.21)0.1940.84(0.40, 1.80)0.6580.64(0.29, 1.40)0.258
  1. † Computes the proportion of households having access to piped-water in the unique community surrounding each participant in the study (Figure 3). The Quintile (Q) ranges (min–max) are: Q1: 0–36; Q2: 37–59; Q3: 60–75; Q4: 76–92; Q5: 93–100, ‡ Corresponding values for a model in which community-level piped-water coverage is used as a continuous variable: a 1% increase in the coverage of piped-water in the surrounding community, was independently associated with a 2.5% decrease in the odds of a Schistosoma haematobium infection (aHR=0.975; 95% CI: 0.966, 0.985; p-value<0.001).

Table 4
Multivariable model examining the socio-demographic predictors of Schistosoma haematobium infection stratified by gender.

Model 1 includes all variables in the model. In Model 2, piped water coverage in the immediate community surrounding each participant has been substituted with household-level piped water covariate.

https://doi.org/10.7554/eLife.33065.010
Model 1: Community-level coverage of piped waterModel 2: Household level access to piped water
FemalesMalesFemalesMales
CovariateaOR§ (95% CI)P-valueaOR§ (95% CI)P-valueaOR§ (95% CI)P-valueaOR§ (95% CI)P-value
Community piped water quintiles (vs Lowest)†
 20.24 (0.10–0.57)‡0.0020.56 (0.35–0.90)‡0.017
 30.21 (0.08–0.55)0.0020.37 (0.17–0.78)0.011
 40.16 (0.06–0.45)0.0010.16 (0.08–0.35)<0.001
 50.07 (0.02–0.20)<0.0010.17 (0.08–0.36)<0.001
Household access to water (vs No)
 Yes0.38 (0.19–0.75)0.0050.76 (0.41–1.42)0.379
Age at Testing
 Per unit9.62 (1.62–57.22)0.0141.24 (1.08–1.42)0.0036.85 (1.56–30.16)0.0111.20 (1.06–1.37)0.005
Age2
 Per unit0.92 (0.85–0.99)0.0240.93 (0.88–0.99)0.017
Toilet in household (vs No)
 Yes1.59 (0.88–2.88)0.1211.12 (0.70–1.86)0.6671.26 (0.75–2.14)0.3711.11 (0.66–1.86)0.687
Household Assets quintile (vs Poorest)
 20.46 (0.24–0.90)0.0261.27 (0.75–2.08)0.3440.55 (0.30–0.99)0.0451.24 (0.76–2.04)0.376
 30.35 (0.13–0.89)0.0321.21 (0.72–1.99)0.4530.35 (0.15–0.84)0.0191.19 (0.71–1.98)0.502
 40.42 (0.18–1.00)0.0611.18 (0.55–2.28)0.6460.40 (0.18–0.91)0.0290.94 (0.46–1.94)0.885
 50.66 (0.32–1.37)0.3070.90 (0.43–1.75)0.7690.54 (0.25–1.17)0.1150.71 (0.34–1.44)0.337
 Missing0.57 (0.20–1.66)0.2881.10 (0.38–2.94)0.8150.40 (0.12–1.35)0.1350.91 (0.34–2.42)0.845
  1. §All estimates simultaneously adjusted for landcover class, distance to water, altitude, slope, treatment in the last 12 months and school grade.

    † Computes the proportion of households having access to piped-water in the unique community surrounding each participant in the study (Figure 3). The Quintile (Q) ranges (min–max) are: Q1: 0–36; Q2: 37–59; Q3: 60–75; Q4: 76–92; Q5: 93–100

Using a continuous measure of community-level piped water coverage, we found that every 1% increase in the coverage of community piped water was independently associated with a 2.5% decrease in the odds of infection (aOR = 0.975, 95% CI 0.966–0.985, p<0.001; see Figure 4b). The relationship is a non-linear one because it is derived from the piped water coverage regression coefficient (–0.025) and then exponentiated to derive an odds ratio. The resulting relationship would imply that relatively modest levels of piped water coverage in a given community could have a large impact on reduction in the odds of infection. For example, a piped water coverage of 30% would correspond to a greater than 50% decrease in the odds of infection for individuals living in that community (Figure 4b). We also found that the result of our primary analysis (ie effect of community-level piped water coverage) was robust to the exact specification of the community (kernel) evaluating the piped water coverage around each member of the cohort. For example, when we increased the size of the kernel’s search radius from 2 km to 3 km then we found that both the continuous (aOR = 0.97, 95% CI: 0.96–0.98, p<0.001) and categorical effect-size estimates (quintile five versus quintile one aOR = 0.12, 95% CI: 0.06–0.24, p<0.001) remained virtually unchanged. Further, we report similar results for a parallel set of analyses (using a Poisson model) to derive risk ratios of infection (Figure 4—figure supplement 1, Supplementary file 2).

For a comparative exercise, we replaced the mean community-level piped water exposure (2001–2007) exposure variable with the piped water coverage only in the year of the parasitological survey (2007). As anticipated, we find in this model that the resulting risk-reduction is attenuated (by 44%) but nevertheless still significant (aOR = 0.986, 95% CI: 0.98–0.99, p=0.004), providing further empirical support for the analytical approach taken in the longitudinal measurement of mean piped water coverage in the years leading up to the parasitological survey.

As a secondary objective, we analysed the impact of household access to piped water on risk of urogenital schistosomiasis infection. We found that this exposure was independently associated with a significant reduction in the odds of infection (aOR = 0.54, 95% CI: 0.33–0.89, p=0.017; Model 2 of Table 3), which was considerably smaller than its community-level counterpart. Once again this reduction was driven mostly by the reduction in risk in girls (aOR = 0.38, 95% CI: 0.19–0.75, p=0.005; Model 2 of Table 4), which was substantially larger than the corresponding impact in boys (aOR = 0.76, 95% CI: 0.41–1.42, p=0.38; Model 2 of Table 4, see also Figure 5). When formally testing this as an interaction effect in the combined-sex model, the reduction in risk was significantly greater for girls (p=0.014).

Comparison of adjusted odds of Schistosoma haematobium infection (95% CI) in participants living in households with access to piped water (relative to participants without household-access to piped water) by sex.

The resulting risk estimates are adjusted for age, household assets, toilet in household, landcover class, distance to water body, altitude, slope, treatment in the last 12 months and school grade. Standard errors are adjusted for clustering by school and grade.

https://doi.org/10.7554/eLife.33065.011

Our risk-reduction estimate for access to piped water at the household level (aOR = 0.54) is close to the Schistosoma haematobium specific aOR of 0.57 generated as part of Grimes et al. systematic review (Grimes et al., 2014) comprising 17 previous studies (Abou-Zeid et al., 2012; Al-Waleedi et al., 2013; Awoke et al., 2013; Dame et al., 2006; Dawet, 2012; Farooq et al., 1966; Howarth et al., 1988; Knopp et al., 2013; Nworie et al., 2012; Reuben et al., 2013; Sady et al., 2013) (Figure 6). The revised overall risk-reduction estimate changes marginally after inclusion of our results and rerunning the meta-analysis, to 0.56 (95% CI: 0.45–0.70). The fact that our effect-size estimate for this secondary analysis is close to that obtained from the systematic review and meta-analysis (based on 17 data-points), lends further credence to our findings - including the risk-reduction result reported from our primary analysis (described above).

Forest plot of Schistosoma haematobium infection according to household/individual level access to a safe water.

The data are taken from Grimes et al. systematic review (Grimes et al., 2014), based on 17 data-points (Abou-Zeid et al., 2012; Al-Waleedi et al., 2013; Awoke et al., 2013; Dame et al., 2006; Dawet, 2012; Farooq et al., 1966; Howarth et al., 1988; Knopp et al., 2013; Nworie et al., 2012; Reuben et al., 2013; Sady et al., 2013) to which we have added our study results. The sizes of the squares represent the weight given to each study, the rhombus is the effect-size with the black lines representing the 95% confidence intervals. The overall rhombus represents the combined effect-size, with the results of this study shown in red.

https://doi.org/10.7554/eLife.33065.012

Finally, our IV analysis to address any potential residual confounding confirms all of our main results. Our instrumental variable met the relevance condition of instrumental variable estimation, as shown by the size of the F-test for instrument strength (F-statistic = 55). In the IV analysis using a linear model, a 1% increase in the coverage of piped water in the surrounding community was independently associated with a 0.48% decrease in the risk of urogenital schistosomiasis infection (Coef= –0.0048, p<0.001; Model 2 of Supplementary file 2). This effect-size estimate is larger than the risk reduction observed without the IV, which was 0.30% (Coef = –0.003, p<0.001; Model 1 of Supplementary file 2). Household access to piped water was associated with a 15.6% reduction in the risk of infection (Coef = –0.156, p=0.035; Model 4 of Supplementary file 2). This effect-size is higher than the estimate in the analysis without the IV, which was 9.2% (Coef = –0.092; p=0.013; Model 3 of Supplementary file 2). These findings demonstrate that the reductions in urogenital schistosomiasis risk are not overestimated as a result of other unmeasured variables or residual confounding.

Discussion

We have linked longitudinal population-based exposures to schistosomiasis outcomes in a large cohort of children to provide a rigorous quantification of the relationship between piped water coverage and the risk of urogenital schistosomiasis infection in a typical rural sub-Saharan African population. This approach provides a strong indicator of the overall impact of the introduction of piped water because it captures both the direct and indirect protection conferred by the scale-up. Our results reveal a large impact of piped water scale-up on urogenital schistosomiasis infection. This impact is substantially greater than previous results had suggested (Grimes et al., 2014), based on measurement of access to safe water supplies at the household or individual level. In our study, children living in communities with a high coverage of piped water were eight times less likely to be infected relative to those living in areas with little or no access to piped water. Every 1% increase in the coverage of piped water in the surrounding local community was associated with a corresponding 2.5% decrease in the odds of infection. As a secondary objective, we also provide a robust household-level piped water risk-reduction estimate (aOR = 0.54) after controlling for key socio-economic and environmental predictors of infection. This result is also important, because the relationship between piped water and schistosomiasis infection at the household or individual level would be confounded in different ways by multiple socio-economic and environmental factors, as the recent systematic review of the existing literature points out (Grimes et al., 2014). Importantly, we validate the veracity of all of the above findings using an instrumental variable approach to control for the influence of residual unobserved confounding that could have biased our effect-size estimates.

Our results show that nearly 1 in 5 primary school children (median age = 11 years) were infected with urogenital schistosomiasis and nearly 1 in 10 were heavily infected. These levels of urogenital schistosomiasis infection are markedly lower than the overall prevalence of 60–70% recorded in the general area thirty years earlier (Gear et al., 1980) and 69% prevalence recorded in a more recent survey undertaken in a neighbouring district at similar altitude (<300 metres above sea level) but with limited access to piped water (Saathoff et al., 2004). Importantly, our geographical analyses revealed remarkable spatial heterogeneity in urogenital schistosomiasis prevalence across the relatively small surveillance area. For example, the prevalence of infection was >80% in some communities where the roll-out of piped water had not taken place. At a micro-geographical level, we identified clear spatial clustering of infections around four specific open water collection sites. This heterogeneity in schistosomiasis transmission is consistent with the results of a 2004 study conducted 200 km north of the surveillance area (Saathoff et al., 2004).

Another important finding to emerge from the work is that the greater burden of urogenital schistosomiasis infection is borne by boys who had more than twice the prevalence of infection in comparison to girls. In addition, there were clear sex differentials in the pathways of infection and the impact of piped water at a household and community level. Whilst environmental variables (distance to water-body, altitude and landcover class) were highly predictive of infection in both sexes, higher socio-economic status was independently associated with a decrease in the odds of infection for girls but not for boys. At the community level, piped water coverage was associated with significantly greater reductions in the infection risk for girls. Similarly, household access to piped water was associated with a significant reduction in risk in girls but not in boys. These sex differences at both the community and household level could reflect the cultural context in which girls are responsible for the majority of the domestic water collection and/or a situation where boys swim more often in infected water bodies (Ngorima et al., 2008; Thomassen Morgas et al., 2010).

The differential in the community/household risk-reduction result for boys suggests that the scale-up of piped water has conferred substantial indirect protection to boys through the lowering of the overall parasite reservoir in a particular community. Had we only analysed access to piped water as a single household-level exposure, we would have failed to demonstrate that the scale-up this intervention offered any significant protection to boys. Boys are exposed in the community to free swimming cercariae released by snails infected not just by themselves but also by their family members and neighbours. Therefore, in order to fully capture the impact of piped water in reducing schistosomiasis transmission risk, it is crucial to look beyond the household and consider the piped water coverage level in the surrounding local community. This view is also expressed by other authors (Esrey and Habicht, 1986). The greater than four-fold ratio in the community-level (aOR = 0.12) versus household-level (aOR = 0.54) risk-reduction estimate provides empirical support for this viewpoint. Whilst no such community-level analyses have previously been conducted for Schistosoma haematobium, a study in China found Schistosoma japonicum infection rates to be significantly lower in communities where more than 50% of people used ‘hygienic lavatories’ (Yang et al., 2009).

Our study design has several strengths that overcome many of the limitations typically associated with school-based surveys in which the exposures and disease outcomes are measured at the same point in time. By nesting the parasitological survey within a population cohort followed over seven years preceding the measurement of the disease, we utilise longitudinal data collected both on the participants of the parasitological survey, as well as the entire population to provide a rigorous quantification of the impact of the introduction of piped water on risk of urogenital schistosomiasis infection. A key strength of the work is that rather than treating access to piped water as a household/individual level exposure as in previous work (Abou-Zeid et al., 2012; Al-Waleedi et al., 2013; Awoke et al., 2013; Dame et al., 2006; Dawet, 2012; Farooq et al., 1966; Howarth et al., 1988; Knopp et al., 2013; Nworie et al., 2012; Reuben et al., 2013; Sady et al., 2013), we were able to quantify the clear dose-response in the relationship between community coverage of piped water and risk of infection. In addition, we used a novel spatial methodology to perform a micro-geographical clustering analysis, based on the exact location of the participants’ households, to pinpoint clear transmission foci around four open water collection points. Crucially, through the use of high resolution satellite imagery, ancillary environmental data and variables collated as part of the ongoing population cohort, we controlled for the key individual, social, environmental and community determinants of schistosomiasis infection which are vital to quantifying the independent impact of piped water. Notwithstanding the above, we wanted to completely exclude the possibility of residual confounding (on the basis of unmeasured social and environmental factors) leading to an over-estimate of the impact of the introduction of piped water on risk of infection. To address this possibility, we performed an instrumental variables (IV) analysis using the year that piped water was first introduced into the community as the instrument. Our results remained robust in this analysis: all effect-size estimates remained significant and increased to some extent in comparison to the main results. These findings demonstrate that the reductions in risk reported in this paper are robust and may be regarded as being marginally conservative.

However, the study has some limitations. The results are based on a single urine sample and as a result we would likely have missed some light infections (Kosinski et al., 2011). The true prevalence of infection is therefore probably higher than 16.8% reported in this study. A further limitation is that, while South Africa is a country with near-universal primary school attendance (Lund, 2011), we could be missing some children who did not attend school on the day of testing because they have schistosomiasis-related morbidity. This would lead to an under-estimate of the prevalence of urogenital schistosomiasis. In addition, the high mobility of this population (Dobra et al., 2017) means that children would be moving between communities, resulting in an ascertainment error leading ultimately to a reduction in the true effect-size, which would bias the result towards the null hypothesis of no piped water impact. In addition, infrequent, ad-hoc schistosomiasis ‘test and treat’ activities conducted in schools in the highest-risk areas (where piped water coverage was lowest), could theoretically bias the effect-size towards the null hypothesis of no piped water impact. However, introduction of such a bias by these activities is unlikely given the large piped water effect-sizes that we report and the fact that only 2% of participants in our survey reported receiving praziquantel in the previous 12 months (which we use as a control variable in the statistical model). In any event, any related ‘residual confounding’ effect if it did exist, would be controlled for in the IV approach which confirmed all of our main findings.

Schistosomiasis and soil-transmitted helminthiasis affect more than 1.5 billion of the world's poorest people (Lai et al., 2015). Without a change in strategy to control these diseases, the population of sub-Saharan Africa will lose an estimated US$3.5 billion of economic productivity every year, which is comparable to recent acute epidemics, including the 2014 Ebola and 2015 Zika epidemics (Lo et al., 2017). Our analytical approach (quantifying the combined direct and indirect impact of the scale-up in piped water) has captured a large measurable association between the introduction of piped water and the reduction in risk of urogenital schistosomiasis in a rural African community. This reduction is substantially greater than had been suggested by previous work (Grimes et al., 2014), which focused on household or individual level access to safe water. The continued expansion of safe water supplies, decreasing prevalence and highly focal nature of disease transmission, suggest that the disease is retreating and could be vulnerable to eradication in similar African contexts. Recent modelling work indicates that the current WHO treatment goals (which aim to mass-treat 75% of at-risk children) should be successful in eradicating schistosomiasis transmission if continued for a decade in areas of low to moderate transmission (Anderson et al., 2015). However, once mass treatment is suspended, the prevalence of infection can rebound to pre-control levels over a period of 25–30 years (Gurarie et al., 2015), emphasising the need for a fully integrated approach that includes the provision of safe water supplies, community health education as well as the possibility of focal snail control (Knopp et al., 2013; King et al., 2015).

The overall moderate levels of infection (significantly lower than historical surveys (Gear et al., 1980) combined with the extreme inter-host and geographical heterogeneity as well as ever increasing access to safe water suggest that sustained mass treatment campaigns in such populations could drive schistosomiasis towards low levels of endemicity (Crompton and WHO, 2006). Our results emphasise the need to rapidly increase the provision of safe water supplies to schistosomiasis endemic areas, where the resulting impact will extend well beyond prevention of transmission of schistosomiasis to diarrhoeal diseases, for example (Esrey and Habicht, 1986); as well as having substantial related economic benefits (Lo et al., 2017). Scale-up in the provision of safe water and sanitation, along with the small cost of drugs combined with the renewed commitment to increase funding for neglected tropical diseases, including schistosomiasis, means that there is a real opportunity to push back the disease from areas where it has long been endemic.

Materials and methods

Schistosomiasis in South Africa

Request a detailed protocol

In South Africa, Schistosoma haematobium is endemic over an area of 320 000 km2 including parts of six of the country’s nine provinces to the east of latitude 26°E (Doumenge et al., 1987). Highest prevalence and intensities occur in the lower lying areas of Limpopo, Mpumalanga and KwaZulu-Natal provinces. Estimates of the number of children infected have been based largely on the 38 year old Atlas of Bilharzia (Gear et al., 1980), and range from 1.5 million to 4 million (Moodley et al., 2003). The most frequently infected age group is 5 to 15 years after which prevalence declines to a generally low level in adulthood. Morbidity of the urinary tract, which is greatest in children with high worm and egg burdens, follows a similar pattern. Considerable gynaecological morbidity, also due to Schistosoma haematobium and known as female genital schistosomiasis, has been demonstrated in KwaZulu-Natal (Hegertun et al., 2013) and is probably present over much of the endemic area. Schistosoma haematobium has a significantly wider distribution in South Africa than its congener Schistosoma mansoni which causes intestinal or rectal schistosomiasis (Appleton and Madsen, 2012). KwaZulu-Natal is the best researched of the endemic provinces in South Africa and historically Schistosoma haematobium prevalence and intensities have reached 70% or more and 200 eggs/10 ml (geometric mean) respectively between the coast and an altitude of 300 meters above sea level. Both prevalence and intensity have been observed to decrease with increasing altitude to around 20% and <10 eggs/10 ml above 800 metres above sea level. Schistosoma mansoni is less common and limited to patches of transmission, mostly in the province’s lowlands (Appleton and Kvalsvig, 2006). In the coastal plain of north-eastern KwaZulu-Natal (where the survey was conducted), Schistosoma mansoni is largely absent (Gear et al., 1980), for reasons discussed in detail elsewhere (Appleton and Madsen, 2012).

Study area

Request a detailed protocol

The Africa Health Research Institute (AHRI) runs one of Africa’s largest ongoing population-based cohorts in the Umkhanyakude district of the northern KwaZulu-Natal province (Figure 7). Approximately 250 km north of the city of Durban, the surveillance area is 438 km2 in size with a population of approximately 90,000 isiZulu-speaking people and 10,000 households. AHRI has collected detailed demographic, socio-economic and behavioural data every four to six months from residents in the surveillance area since 2000 (Tanser et al., 2008). All households have been positioned to an accuracy of <2 m using differential global positioning systems (Tanser et al., 2001).

Location of the study area in South Africa.
https://doi.org/10.7554/eLife.33065.013

The surveillance area is designated as a coastal lowland (<300 metres above sea level) and falls within the endemic schistosomiasis belt of South Africa (Gear et al., 1980; Moodley et al., 2003). The medication praziquantel is available in some of the local clinics to treat schistosomiasis symptoms. In the period leading up to the parasitological survey, there were no mass treatment campaigns, but the local Department of Health did undertake occasional, ad-hoc ‘test and treat’ activities in a small number of schools in high-risk areas at different points in time. The population has been exposed to a devastating HIV epidemic. In 2010, the HIV prevalence was 24% among all adults and approaching 60% in some age-groups (Tanser et al., 2013; Vandormael et al., 2018).

Sample size

Request a detailed protocol

The sample size calculation was based on the number of participants required to test the impact of access to piped water on infection with Schistosoma haematobium stratified by sex, as well as provide good geographical coverage of the surveillance area. Given two grades sampled at each of the 33 schools (50 pupils per grade) in the study area and an assumption of 70% consent and 90% linkage to the population cohort, we would attain an estimated sample size of 2079. We then assumed that the piped water exposure occurred in 50% of participants and that 30% of pupils tested for Schistosoma haematobium were infected. This prevalence estimate was based on historical parasitological surveys in surrounding areas (Gear et al., 1980), more recent surveys in the province of KwaZulu-Natal (Jinabhai et al., 2001; Taylor et al., 2001) and expected reductions in prevalence as a result of the piped water roll-out in the study area during the preceding decade. Under these assumptions, we would attain 99% power to detect a relative risk of ≥1.5 (95% CI) in boys and girls respectively as well as provide good geographical coverage of the study area to facilitate identification of micro-geographical clustering of infections.

Ethics approval

Request a detailed protocol

Ethical clearance for the study was obtained from the Biomedical Research Ethics Committee of the University of KwaZulu-Natal (reference number E165/05). Prior to testing, we were required to obtain both written informed consent from parents and written assent from children participating in the survey (described further below).

The Schistosoma haematobium survey

Request a detailed protocol

We conducted a parasitological survey among 2105 consenting children (median age = 11) in their fifth and sixth year of primary school at all of the 33 primary schools in the surveillance area between March and December 2007. We chose to sample these two school grades based on likely highest levels of Schistosoma haematobium infection as well as to facilitate a follow-up survey one year later in the same cohort of children (before completion of their primary school education at the end of year 7). Prior to the study, we initiated an extensive community preparation process that included consultations with headmasters and the governing bodies of participating schools, and delivered presentations to parents and teachers at school meetings. We distributed consent forms to the parents through their children the day before testing. Assent forms were signed by the children on the day of testing. We excluded children who did not have both parental informed consent and child assent from the parasitological survey.

After receiving the written consent and assent forms, we gave participants a cup of juice to drink and a 500 ml honey jar to provide a full urine sample between the hours of 10h00 and 12h00. We labelled each urine sample with a unique bar code and performed a rapid specimen test for the presence of micro haematuria using urinalysis reagent strips (Bayer Uristix). Following national guidelines, we treated those with detectable micro haematuria with a single oral dose of praziquantel (Biltricide, Bayer) on the basis of body weight (40 mg per kg) (Saathoff et al., 2004; WHO Expert Committee, 2002). We sent the urine samples for microscopic analysis (described in detail below) and participants found to be falsely negative at baseline were treated at a follow-up visit. Treatment was administered and strictly supervised by qualified nurses. After at least 30 days, we revisited the schools and gave participants (who were treated at the first visit) a second dose of praziquantel if micro haematuria was again detected. We administered a brief questionnaire after testing and asked participants if they had been treated with praziquantel in the last 12 months.

Laboratory procedures

Request a detailed protocol

We aliquoted the urine sample into two 10 ml sub-samples and diluted with a preservation solution (Formalin Methiolate concentration), and uniquely bar-coded each aliquot. The microscopic analysis was done using the filtration technique as described by the World Health Organization (WHO, 1991). All eggs on the slides were counted and expressed as eggs per centilitre (10 ml) of urine. Specimens of less than 10 ml were measured before filtration and the number of eggs per 10 ml calculated. Samples with <50 eggs/10 ml of urine specimen were categorised as light infection while those with ≥50 eggs/10 ml of urine specimen were categorised as heavy infection (WHO Expert Committee, 2002).

Derivation of community-level piped water coverage

Request a detailed protocol

We used all five separate population-based surveys (2001, 2003/2004, 2005, 2006, 2007) conducted up to the year of parasitological survey (2007) in the surveillance area to measure access to piped water (Tanser et al., 2008). There were on average 10,267 households per survey. All households in the surveillance area were asked; ‘which is the most often used source of drinking-water in this household?’ The household key informant was asked to select the ‘most common’ source from a list of 11 different sources of water listed. Households were deemed to have access to piped water if they reported their most common source of drinking water to be piped water delivered either directly to the household or from a public tap.

We used this population-based data to derive a continuous measure of the mean piped water coverage in the unique community surrounding each participant’s household (between 2001 and 2007). This community-level piped water measure was first constructed by means of a moving two-dimensional Gaussian kernel of 2 km search radius (Tanser, 2006; Tanser et al., 2009). For comparative purposes, we also produced an equivalent analysis using a kernel search radius of 3 km. First, the piped water data were superimposed on a geographic representation of the surveillance area consisting of a grid of 30m × 30 m pixels. Next, the kernel moves systematically across the grid and calculates a Gaussian weighted estimate of the piped water coverage for the unique neighbourhood around each and every pixel on the grid. The methodology therefore provides a sensitive and realistic estimate of piped water coverage in each household’s unique surrounding ‘virtual community’. We located each participant to an exact homestead of residence and extracted the corresponding mean piped water coverage (2001–2007) exposure estimate in the surrounding local community.

The community-level piped water coverage exposure measure (2001–2007) fully captures the dynamic nature of the roll-out and availability of piped water in the seven years leading up to the parasitological survey. This is important for two reasons. Firstly, it ensures that the exposure measurement (piped water coverage) precedes the measurement of the disease outcome (urogenital schistosomiasis). This is vital for a chronic infection like schistosomiasis in which children may have been infected for several years before the disease outcome was measured. Secondly, there are occasional instances of water schemes being implemented in the surveillance area but not working consistently for periods of time. In this regard, a community having access to piped water in one year does not necessarily guarantee access to piped water in the following year. Thus, our community-level piped water coverage measure fully captures the vagaries of access to piped water in the years leading up to the parasitological survey.

Statistical analysis

Request a detailed protocol

Of the 2105 children who participated in our parasitological survey, 94% (N = 1,976) could be linked to their corresponding records in the population-based cohort and were included in the statistical analysis. From this we obtained the precise geographical coordinates of each participant’s place of residence as well as detailed individual-level and household-level predictors of infection. Children were considered to be infected with Schistosoma haematobium on the basis of the microscopy result. To facilitate direct comparison with previous work (Grimes et al., 2014), we used a logistic regression as our main statistical model to quantify the impact of access to piped water both at the community and household level respectively, whilst simultaneously controlling for the key environmental and socio-economic predictors of urogenital schistosomiasis infection (described below). However, in light of the expected relatively high prevalence of infection, we also report results for a Poisson regression model to obtain risk ratios of Schistosoma haematobium infection. We undertook the logistic and Poisson regression analyses in Stata 14.1, and utilised the svy suite of commands to account for the clustering of participants in the 33 schools (two grades) from 66 primary sampling units (Kreuter and Valliant, 2007). The failure to correct for clustering violates the independence assumption and can lead to distorted parameter estimates and standard errors (Guo and Zhao, 2000). We stratified our analyses by sex to test for any differential impacts in pathways of infection between males and females.

Finally, we formally compared the resulting risk-reduction estimate from our secondary analysis (ie exposure to piped water at the household level) to the effect-size estimate generated as part of Grimes et al. systematic review and meta-analysis (Grimes et al., 2014) comprising previous studies (Abou-Zeid et al., 2012; Al-Waleedi et al., 2013; Awoke et al., 2013; Dame et al., 2006; Dawet, 2012; Farooq et al., 1966; Howarth et al., 1988; Knopp et al., 2013; Nworie et al., 2012; Reuben et al., 2013; Sady et al., 2013). To do this, we captured the results (i.e., sample sizes, odds ratios, and 95% confidence intervals) from the 17 data-points along with the results from our analysis described above. We then ran a random-effects meta-analysis on this data using the metan package in Stata (Harris et al., 2008).

Spatial analysis

Request a detailed protocol

We used two different and complementary approaches to assess the spatial distribution of urogenital schistosomiasis infection in the study population. Firstly, we used the Gaussian kernel methodology described above (using the precise location of every participant’s household) to produce a prevalence map that varies in continuous geographical space. Using this approach, a median of 166 (IQR = 114–248) participants were evaluated for each cell’s unique neighbourhood (search radius = 3 km).

Secondly, we used a Kulldorff spatial scan statistic (Kulldorff et al., 2005; Kulldorff, 1997) to formally identify and map micro-geographical clusters of unusually high numbers of urogenital schistosomiasis infections that were unlikely to have arisen by chance. The Kulldorff spatial scan statistic is implemented within the SaTScan spatial cluster detection programme (Kulldorff, 2006). We have used this technique used in our previous work in this population to identify spatial clusters of HIV infections (Tanser et al., 2009). A spatial scan statistic is a cluster detection test that is able to both detect the location of clusters and evaluate their statistical significance without problems associated with multiple testing. This is done by gradually scanning a window across space. The statistical theory behind spatial scan statistics is described elsewhere (Kulldorff et al., 2005). Briefly, the spatial scan statistic imposes a circular window on a map and it allows the centre of the circle to move systematically across the study region. For any given position of the centre, the size of the circle changes continuously so that it can take any value between zero and a specified maximum value. For every potential cluster, a likelihood ratio test statistic is used to determine if the number of urogenital schistosomiasis cases within the potential cluster was higher than expected. Expected numbers of cases were calculated on the basis of the null hypothesis of spatial randomness by assuming that the number of cases in each circle is an independent Bernoulli random variable with constant prevalence. The circle with the maximum likelihood is defined as the most likely cluster (implying that it is least likely to have occurred by chance). The maximum observed value of the test statistic for each possible cluster is then compared to the overall distribution of maximum values (Kulldorff, 2006). The p-value of the statistic is obtained through Monte Carlo hypothesis testing (9999 iterations), where the null hypothesis of no clustering is rejected if the simulated p-value is <0.05. We allowed the clusters to overlap by <50% and set the maximum search radius of the circle to be 3 km to facilitate comparison with the Gaussian kernel smoothing approach described above.

Adjustment for observed social and environmental predictors of infection

Request a detailed protocol

Schistosomiasis is a complex parasitic disease characterised by key environmental and social predictors which vary in time and space (Shiff, 2001). It is therefore not possible to accurately quantify the independent impact of piped water without taking into account the extent to which the local environmental and social conditions in a particular area are conducive to transmission. To control for key environmental differences, we extracted altitude, slope and distance to the nearest water body for every household in the surveillance area (Figure 8). To extract altitude and slope, we used the 30 m raster digital elevation model obtained from the South African Surveyor General. Distance to the nearest water body was derived using a combination of 1:50 000 digital data and additional water bodies digitised from 1:10 000 orthorectified aerial photographs. We also used LandSat TM bands 1–7 to map five different ecological units occurring in the surveillance area (Figure 8). The derivation of this landcover classification is described in further detail below. For socio-demographic control variables, we used age, sex, praziquantel in the last 12 months, school grade, presence of the toilet in the household and a household socio-economic assets index (generated in a principal component analysis based on household ownership of 24 assets).

Environmental control variables used in the statistical analysis with locations of 33 schools superimposed (yellow diamonds).

(Top left) Altitude in metres above sea level (MASL) (Top right) Distance to nearest water body (km) (Bottom left) Slope (degrees) (Bottom right) Satellite-derived landcover classification.

https://doi.org/10.7554/eLife.33065.014

To map the ecological units across the study area, we obtained 2007 LandSat multi-spectral, high-resolution (30 m × 30 m) satellite imagery which captures reflectance values at seven regions of electromagnetic spectrum. The data were corrected for radiometric and geometric distortions (Richards and Jia, 2006). We used an iterative self-organising unsupervised classifier ISOCLUST implemented within IDRISI Taiga (Clark University, MA, USA) to classify the image into five landcover classes. The clustering procedure implements a histogram peak cluster analysis technique that searches for frequency peaks higher than those of its immediate neighbours. We used the ‘broad level of generalisation’ option which allocates all pixels to major frequency peaks in the multidimensional histogram (minor frequency peaks are ignored) (Eastman, 2009). The module uses a true maximum likelihood clustering procedure that leads to a more efficient and accurate placement of clusters than either random or systematic placement. These landcover classes were subsequently described using a combination of orthorectified 1:10 000 aerial photographs and physical on-the-ground inspection of 40 sites.

A common problem in the use of high resolution satellite imagery is the high degree of spectral variation and spectral scatter within small areas leading to inconsistencies in the classification process. In addition, from a disease exposure perspective, it makes sense to consider an area greater than immediate the 30 m × 30 m area around the place of residence. We therefore employed an image segmentation procedure in IDRISI Taiga that uses a watershed delineation approach to partition the image into a series of ecological units (segments) on the basis of spectral similarity (using a 30 pixel similarity tolerance) (Eastman, 2009). Each segment in the resulting image was then classified as one of the five broad landcover classes identified in the previous classification result (using a majority rule). This process improves the accuracy of the initial classification and producing a smoother classification result while preserving the boundaries between segments (Richards and Jia, 2006). Two of the resulting landcover classes were then combined into a single class due to their similarity and small numbers of individuals living within them.

Adjustment for unobserved social and environmental predictors of infection

Request a detailed protocol

To test the veracity of our findings, we extended the statistical analysis to address the possibility of residual confounding (via unobserved factors) by carrying out an instrumental variable (IV) analysis. IV analyses allow a quasi-random exposure assignment by isolating the random component of an exposure that is not completely randomly assigned (Baiocchi et al., 2014; Bärnighausen et al., 2017; Angrist et al., 1996). IV analyses require a variable, the instrument, that is significantly associated with the exposure of interest (here: piped water exposure) but is not independently associated with the outcome of interest (here: urogenital schistosomiasis infection). We used the year when piped water was first provided in a traditional Zulu ward, the isigodi, as an instrumental variable to isolate the exogenous or as-good-as random variation in the exposure of interest (piped water). A confluence of factors, each of which is fast-changing in this community, led to the delivery of piped water to a particular local area at a particular point in time, including complex political constellations at the local and higher levels, availability of construction capacity, and presence of particular household members. The joint presence of this constantly-changing multitude of factors leading to piped water delivery at a particular time-point was likely as-good-as random regarding other potential determinants of infection, such as wealth or coverage with health services (UWP Engineers, 2002). It is thus unlikely to be independently correlated with the outcome (urogenital schistosomiasis infection) and likely meets the exclusion restriction (Angrist et al., 1996). For the IV analysis, we used a one-stage maximum likelihood estimation with a linear probability model as implemented in the ivregress command in Stata 14.1 (Baum, 2006). Specifically, we added to the model our instrumental variable—the year that piped water was introduced into the Zulu ward—along with the same environmental and socio-economic control variables used in our primary analysis. To facilitate direct comparison with the IV analysis, we repeated our main analysis using a linear probability model.

Role of the funding source

Request a detailed protocol

The funders (The International Collaboration in Infectious Disease Research and the Wellcome Trust) had no role in the design and conduct of the study, interpretation of the data, or approval of the report. The corresponding author had access to all study data and made the final decision to submit for publication.

References

  1. Book
    1. Baum CF
    (2006)
    An Introduction to Modern Econometrics Using Stata
    Texas: Stata press.
  2. Book
    1. Crompton DWT
    2. WHO
    (2006)
    Preventive Chemotherapy in Human Helminthiasis: Coordinated Use of Anthelminthic Drugs in Control Interventions: A Manual for Health Professionals and Programme Managers
    Geneva: World Health Organization.
    1. Dame J
    2. Banwat E
    3. Egah D
    4. Shabi M
    (2006)
    Schistosomiasis among students in a local government area of Kaduna State in Northern Nigeria
    Highland Medical Research Journal 4:60–69.
  3. Book
    1. Doumenge J-P
    2. Mott KE
    3. Cheung C
    (1987)
    Atlas of the Global Distribution of Schistosomiasis Continued
    Presses universitaires de Bordeaux.
  4. Book
    1. Eastman J
    (2009)
    IDRISI Taiga Guide to GIS and Image Processing
    Worcester: Clark University.
    1. Farooq M
    2. Nielsen J
    3. Samaan SA
    4. Mallah MB
    5. Allam AA
    (1966)
    The epidemiology of Schistosoma haematobium and S. mansoni infections in the Egypt-49 project area. 2. Prevalence of bilharziasis in relation to personal attributes and habits
    Bulletin of the World Health Organization 35:293.
  5. Book
    1. Gear J
    2. Pitchford R
    3. Van Eeden J
    (1980)
    Atlas of Bilharzia in South Africa
    Johannesburg: South African Institute for Medical Research.
    1. Harris R
    2. Bradburn M
    3. Deeks J
    (2008)
    Metan: fixed-and random-effects meta-analysis
    The Stata Journal 8:3–28.
    1. Howarth SE
    2. Wilson JM
    3. Ranaivoson E
    4. Crook SE
    5. Denning AM
    6. Hutchings MS
    (1988)
    Worms, wells and water in western Madagascar
    The Journal of Tropical Medicine and Hygiene 91:255–264.
    1. Kreuter F
    2. Valliant R
    (2007)
    A survey on survey statistics: What is done and can be done in Stata
    The Stata Journal 1:1–21.
    1. Kulldorff M
    (1997) A spatial scan statistic
    Communications in Statistics - Theory and Methods 26:1481–1496.
    https://doi.org/10.1080/03610929708831995
    1. Murray CJ
    2. Vos T
    3. Lozano R
    4. Naghavi M
    5. Flaxman AD
    6. Michaud C
    7. Ezzati M
    8. Shibuya K
    9. Salomon JA
    10. Abdalla S
    11. Aboyans V
    12. Abraham J
    13. Ackerman I
    14. Aggarwal R
    15. Ahn SY
    16. Ali MK
    17. Alvarado M
    18. Anderson HR
    19. Anderson LM
    20. Andrews KG
    21. Atkinson C
    22. Baddour LM
    23. Bahalim AN
    24. Barker-Collo S
    25. Barrero LH
    26. Bartels DH
    27. Basáñez MG
    28. Baxter A
    29. Bell ML
    30. Benjamin EJ
    31. Bennett D
    32. Bernabé E
    33. Bhalla K
    34. Bhandari B
    35. Bikbov B
    36. Bin Abdulhak A
    37. Birbeck G
    38. Black JA
    39. Blencowe H
    40. Blore JD
    41. Blyth F
    42. Bolliger I
    43. Bonaventure A
    44. Boufous S
    45. Bourne R
    46. Boussinesq M
    47. Braithwaite T
    48. Brayne C
    49. Bridgett L
    50. Brooker S
    51. Brooks P
    52. Brugha TS
    53. Bryan-Hancock C
    54. Bucello C
    55. Buchbinder R
    56. Buckle G
    57. Budke CM
    58. Burch M
    59. Burney P
    60. Burstein R
    61. Calabria B
    62. Campbell B
    63. Canter CE
    64. Carabin H
    65. Carapetis J
    66. Carmona L
    67. Cella C
    68. Charlson F
    69. Chen H
    70. Cheng AT
    71. Chou D
    72. Chugh SS
    73. Coffeng LE
    74. Colan SD
    75. Colquhoun S
    76. Colson KE
    77. Condon J
    78. Connor MD
    79. Cooper LT
    80. Corriere M
    81. Cortinovis M
    82. de Vaccaro KC
    83. Couser W
    84. Cowie BC
    85. Criqui MH
    86. Cross M
    87. Dabhadkar KC
    88. Dahiya M
    89. Dahodwala N
    90. Damsere-Derry J
    91. Danaei G
    92. Davis A
    93. De Leo D
    94. Degenhardt L
    95. Dellavalle R
    96. Delossantos A
    97. Denenberg J
    98. Derrett S
    99. Des Jarlais DC
    100. Dharmaratne SD
    101. Dherani M
    102. Diaz-Torne C
    103. Dolk H
    104. Dorsey ER
    105. Driscoll T
    106. Duber H
    107. Ebel B
    108. Edmond K
    109. Elbaz A
    110. Ali SE
    111. Erskine H
    112. Erwin PJ
    113. Espindola P
    114. Ewoigbokhan SE
    115. Farzadfar F
    116. Feigin V
    117. Felson DT
    118. Ferrari A
    119. Ferri CP
    120. Fèvre EM
    121. Finucane MM
    122. Flaxman S
    123. Flood L
    124. Foreman K
    125. Forouzanfar MH
    126. Fowkes FG
    127. Fransen M
    128. Freeman MK
    129. Gabbe BJ
    130. Gabriel SE
    131. Gakidou E
    132. Ganatra HA
    133. Garcia B
    134. Gaspari F
    135. Gillum RF
    136. Gmel G
    137. Gonzalez-Medina D
    138. Gosselin R
    139. Grainger R
    140. Grant B
    141. Groeger J
    142. Guillemin F
    143. Gunnell D
    144. Gupta R
    145. Haagsma J
    146. Hagan H
    147. Halasa YA
    148. Hall W
    149. Haring D
    150. Haro JM
    151. Harrison JE
    152. Havmoeller R
    153. Hay RJ
    154. Higashi H
    155. Hill C
    156. Hoen B
    157. Hoffman H
    158. Hotez PJ
    159. Hoy D
    160. Huang JJ
    161. Ibeanusi SE
    162. Jacobsen KH
    163. James SL
    164. Jarvis D
    165. Jasrasaria R
    166. Jayaraman S
    167. Johns N
    168. Jonas JB
    169. Karthikeyan G
    170. Kassebaum N
    171. Kawakami N
    172. Keren A
    173. Khoo JP
    174. King CH
    175. Knowlton LM
    176. Kobusingye O
    177. Koranteng A
    178. Krishnamurthi R
    179. Laden F
    180. Lalloo R
    181. Laslett LL
    182. Lathlean T
    183. Leasher JL
    184. Lee YY
    185. Leigh J
    186. Levinson D
    187. Lim SS
    188. Limb E
    189. Lin JK
    190. Lipnick M
    191. Lipshultz SE
    192. Liu W
    193. Loane M
    194. Ohno SL
    195. Lyons R
    196. Mabweijano J
    197. MacIntyre MF
    198. Malekzadeh R
    199. Mallinger L
    200. Manivannan S
    201. Marcenes W
    202. March L
    203. Margolis DJ
    204. Marks GB
    205. Marks R
    206. Matsumori A
    207. Matzopoulos R
    208. Mayosi BM
    209. McAnulty JH
    210. McDermott MM
    211. McGill N
    212. McGrath J
    213. Medina-Mora ME
    214. Meltzer M
    215. Mensah GA
    216. Merriman TR
    217. Meyer AC
    218. Miglioli V
    219. Miller M
    220. Miller TR
    221. Mitchell PB
    222. Mock C
    223. Mocumbi AO
    224. Moffitt TE
    225. Mokdad AA
    226. Monasta L
    227. Montico M
    228. Moradi-Lakeh M
    229. Moran A
    230. Morawska L
    231. Mori R
    232. Murdoch ME
    233. Mwaniki MK
    234. Naidoo K
    235. Nair MN
    236. Naldi L
    237. Narayan KM
    238. Nelson PK
    239. Nelson RG
    240. Nevitt MC
    241. Newton CR
    242. Nolte S
    243. Norman P
    244. Norman R
    245. O'Donnell M
    246. O'Hanlon S
    247. Olives C
    248. Omer SB
    249. Ortblad K
    250. Osborne R
    251. Ozgediz D
    252. Page A
    253. Pahari B
    254. Pandian JD
    255. Rivero AP
    256. Patten SB
    257. Pearce N
    258. Padilla RP
    259. Perez-Ruiz F
    260. Perico N
    261. Pesudovs K
    262. Phillips D
    263. Phillips MR
    264. Pierce K
    265. Pion S
    266. Polanczyk GV
    267. Polinder S
    268. Pope CA
    269. Popova S
    270. Porrini E
    271. Pourmalek F
    272. Prince M
    273. Pullan RL
    274. Ramaiah KD
    275. Ranganathan D
    276. Razavi H
    277. Regan M
    278. Rehm JT
    279. Rein DB
    280. Remuzzi G
    281. Richardson K
    282. Rivara FP
    283. Roberts T
    284. Robinson C
    285. De Leòn FR
    286. Ronfani L
    287. Room R
    288. Rosenfeld LC
    289. Rushton L
    290. Sacco RL
    291. Saha S
    292. Sampson U
    293. Sanchez-Riera L
    294. Sanman E
    295. Schwebel DC
    296. Scott JG
    297. Segui-Gomez M
    298. Shahraz S
    299. Shepard DS
    300. Shin H
    301. Shivakoti R
    302. Singh D
    303. Singh GM
    304. Singh JA
    305. Singleton J
    306. Sleet DA
    307. Sliwa K
    308. Smith E
    309. Smith JL
    310. Stapelberg NJ
    311. Steer A
    312. Steiner T
    313. Stolk WA
    314. Stovner LJ
    315. Sudfeld C
    316. Syed S
    317. Tamburlini G
    318. Tavakkoli M
    319. Taylor HR
    320. Taylor JA
    321. Taylor WJ
    322. Thomas B
    323. Thomson WM
    324. Thurston GD
    325. Tleyjeh IM
    326. Tonelli M
    327. Towbin JA
    328. Truelsen T
    329. Tsilimbaris MK
    330. Ubeda C
    331. Undurraga EA
    332. van der Werf MJ
    333. van Os J
    334. Vavilala MS
    335. Venketasubramanian N
    336. Wang M
    337. Wang W
    338. Watt K
    339. Weatherall DJ
    340. Weinstock MA
    341. Weintraub R
    342. Weisskopf MG
    343. Weissman MM
    344. White RA
    345. Whiteford H
    346. Wiebe N
    347. Wiersma ST
    348. Wilkinson JD
    349. Williams HC
    350. Williams SR
    351. Witt E
    352. Wolfe F
    353. Woolf AD
    354. Wulf S
    355. Yeh PH
    356. Zaidi AK
    357. Zheng ZJ
    358. Zonies D
    359. Lopez AD
    360. AlMazroa MA
    361. Memish ZA
    (2012) Disability-adjusted life years (DALYs) for 291 diseases and injuries in 21 regions, 1990-2010: a systematic analysis for the Global Burden of Disease Study 2010
    The Lancet 380:2197–2223.
    https://doi.org/10.1016/S0140-6736(12)61689-4
    1. Nworie O
    2. Nya O
    3. Anyim C
    4. Okoli C
    5. Okonkwo E
    (2012)
    Prevalence of urinary schistosomiasis among primary school children in Afikpo North Local government area of Ebonyi State
    Annals of Biological Research 3:3894–3897.
    1. Reuben RC
    2. Tanimu H
    3. Musa JA
    (2013)
    Epidemiology of urinary schistosomiasis among secondary school students in Lafia, Nasarawa State, Nigeria
    Journal of Biology, Agriculture and Healthcare  3:73–82.
  6. Book
    1. Richards JA
    2. Jia X
    (2006)
    Remote Sensing Digital Image Analysis: An Introduction
    Berlin: Springer-Verlag.
  7. Book
    1. Shiff C
    (2001)
    Epidemiology of helminth infections
    In: Nelson K, Masters-Williams C, Graham N, editors. Infectious Disease Epidemiology. Gaithersburg: Aspen Publishers. pp. 711–726.
  8. Report
    1. Uniting to Combat Neglected Tropical Diseases
    (2016)
    Reaching the Unreached: Fourth Progress Report of the London Declaration
    Uniting to Combat Neglected Tropical Diseases.
  9. Report
    1. UWP Engineers
    (2002)
    Umkhanyakude District Municipality Water Services Development Plan
    Umkhanyakude District Municipality .
    1. WHO Expert Committee
    (2002)
    Prevention and control of schistosomiasis and soil-transmitted helminthiasis
    World Health Organization Technical Report Series 912:.
  10. Report
    1. WHO
    (1991)
    Basic Laboratory Methods in Medical Parasitology
    Geneva: World Health Organization.
  11. Report
    1. WHO
    (2010)
    Informal Consultation on Expanding Schistosomiasis Control in Africa
    Geneva: World Health Organisation.
  12. Report
    1. WHO
    (2013)
    Sustaining the Drive to Overcome the Global Impact of Neglected Tropical Diseases: Second WHO Report on Neglected Tropical Diseases
    World Health Organisation.
  13. Report
    1. WHO
    (2015)
    Water Sanitation and Hygiene for Accelerating and Sustaining Progress on Neglected Tropical Diseases a Global Strategy 2015–2020
    Geneva: World Health Organisation.

Article and author information

Author details

  1. Frank Tanser

    1. Africa Health Research Institute, KwaZulu-Natal, South Africa
    2. School of Nursing and Public Health, University of KwaZulu-Natal, Durban, South Africa
    3. Centre for the AIDS Programme of Research in South Africa (CAPRISA), University of KwaZulu-Natal, Durban, South Africa
    4. Institute of Epidemiology and Health Care, University College London, London, United Kingdom
    Contribution
    Conceptualization, Data curation, Formal analysis, Supervision, Funding acquisition, Validation, Investigation, Methodology, Writing—original draft, Project administration, Writing—review and editing
    For correspondence
    ftanser@gmail.com
    Competing interests
    No competing interests declared
    ORCID icon "This ORCID iD identifies the author of this article:" 0000-0001-9797-0000
  2. Daniel K Azongo

    Navrongo Health Research Centre, Ghana Health Service, Navrongo, Upper East Region, Ghana
    Contribution
    Investigation, Methodology, Writing—original draft, Project administration, Writing—review and editing
    Competing interests
    No competing interests declared
  3. Alain Vandormael

    1. Africa Health Research Institute, KwaZulu-Natal, South Africa
    2. School of Nursing and Public Health, University of KwaZulu-Natal, Durban, South Africa
    3. KwaZulu-Natal Research Innovation and Sequencing Platform (KRISP), College of Health Sciences, University of KwaZulu-Natal, Durban, South Africa
    Contribution
    Formal analysis, Writing—original draft, Writing—review and editing
    Competing interests
    No competing interests declared
    ORCID icon "This ORCID iD identifies the author of this article:" 0000-0002-5742-0511
  4. Till Bärnighausen

    1. Africa Health Research Institute, KwaZulu-Natal, South Africa
    2. Institute of Epidemiology and Health Care, University College London, London, United Kingdom
    3. Department of Global Health and Population, Harvard T.H. Chan School of Public Health, Boston, United States
    4. Institute for Public Health, University of Heidelberg, Heidelberg, Germany
    Contribution
    Methodology, Writing—original draft, Writing—review and editing
    Competing interests
    No competing interests declared
  5. Christopher Appleton

    School of Life Sciences, University of KwaZulu-Natal, Durban, South Africa
    Contribution
    Conceptualization, Funding acquisition, Investigation, Methodology, Writing—original draft, Writing—review and editing
    Competing interests
    No competing interests declared

Funding

National Institutes of Health

  • Frank Tanser
  • Christopher Appleton

Wellcome Trust

  • Frank Tanser

The funders had no role in study design, data collection and interpretation, or the decision to submit the work for publication.

Acknowledgements

The study was funded through the National Institutes of Health via the International Collaboration in Infectious Disease Research (ICIDR). Funding for the Africa Health Research Institute’s population-based cohort was received from the Wellcome Trust, UK. The authors are indebted to the School-Health team at Hlabisa Hospital for their invaluable assistance in conducting the parasitological survey. The authors wish to express their grateful thanks to Colleen Archer (University of KwaZulu-Natal) for conducting the microscopic analysis and Colin Newell for database support.

Ethics

Human subjects: Ethical clearance for the study was obtained from the Biomedical Research Ethics Committee of the University of KwaZulu-Natal (reference number E165/05). Prior to testing, we were required to obtain both written informed consent from parents and written assent from children participating in the survey

Version history

  1. Received: October 24, 2017
  2. Accepted: January 31, 2018
  3. Version of Record published: February 20, 2018 (version 1)

Copyright

© 2018, Tanser 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,869
    Page views
  • 181
    Downloads
  • 15
    Citations

Article citation count generated by polling the highest count across the following sources: Crossref, PubMed Central, Scopus.

Download links

A two-part list of links to download the article, or parts of the article, in various formats.

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)

  1. Frank Tanser
  2. Daniel K Azongo
  3. Alain Vandormael
  4. Till Bärnighausen
  5. Christopher Appleton
(2018)
Impact of the scale-up of piped water on urogenital schistosomiasis infection in rural South Africa
eLife 7:e33065.
https://doi.org/10.7554/eLife.33065

Share this article

https://doi.org/10.7554/eLife.33065

Further reading

    1. Medicine
    2. Microbiology and Infectious Disease
    3. Epidemiology and Global Health
    4. Immunology and Inflammation
    Edited by Jos WM van der Meer et al.
    Collection

    eLife has published articles on a wide range of infectious diseases, including COVID-19, influenza, tuberculosis, HIV/AIDS, malaria and typhoid fever.

    1. Ecology
    2. Epidemiology and Global Health
    Aleksandra Kovacevic, David RM Smith ... Lulla Opatowski
    Research Article

    Non-pharmaceutical interventions implemented to block SARS-CoV-2 transmission in early 2020 led to global reductions in the incidence of invasive pneumococcal disease (IPD). By contrast, most European countries reported an increase in antibiotic resistance among invasive Streptococcus pneumoniae isolates from 2019 to 2020, while an increasing number of studies reported stable pneumococcal carriage prevalence over the same period. To disentangle the impacts of the COVID-19 pandemic on pneumococcal epidemiology in the community setting, we propose a mathematical model formalizing simultaneous transmission of SARS-CoV-2 and antibiotic-sensitive and -resistant strains of S. pneumoniae. To test hypotheses underlying these trends five mechanisms were built into the model and examined: (1) a population-wide reduction of antibiotic prescriptions in the community, (2) lockdown effect on pneumococcal transmission, (3) a reduced risk of developing an IPD due to the absence of common respiratory viruses, (4) community azithromycin use in COVID-19 infected individuals, (5) and a longer carriage duration of antibiotic-resistant pneumococcal strains. Among 31 possible pandemic scenarios involving mechanisms individually or in combination, model simulations surprisingly identified only two scenarios that reproduced the reported trends in the general population. They included factors (1), (3), and (4). These scenarios replicated a nearly 50% reduction in annual IPD, and an increase in antibiotic resistance from 20% to 22%, all while maintaining a relatively stable pneumococcal carriage. Exploring further, higher SARS-CoV-2 R0 values and synergistic within-host virus-bacteria interaction mechanisms could have additionally contributed to the observed antibiotic resistance increase. Our work demonstrates the utility of the mathematical modeling approach in unraveling the complex effects of the COVID-19 pandemic responses on AMR dynamics.