- Open Access
Modelling the spatial and seasonal distribution of suitable habitats of schistosomiasis intermediate host snails using Maxent in Ndumo area, KwaZulu-Natal Province, South Africa
Parasites & Vectors volume 9, Article number: 572 (2016)
Schistosomiasis is a snail-borne disease endemic in sub-Saharan Africa transmitted by freshwater snails. The distribution of schistosomiasis coincides with that of the intermediate hosts as determined by climatic and environmental factors. The aim of this paper was to model the spatial and seasonal distribution of suitable habitats for Bulinus globosus and Biomphalaria pfeifferi snail species (intermediate hosts for Schistosoma haematobium and Schistosoma mansoni, respectively) in the Ndumo area of uMkhanyakude district, South Africa.
Maximum Entropy (Maxent) modelling technique was used to predict the distribution of suitable habitats for B. globosus and B. pfeifferi using presence-only datasets with ≥ 5 and ≤ 12 sampling points in different seasons. Precipitation, maximum and minimum temperatures, Normalised Difference Vegetation Index (NDVI), Normalised Difference Water Index (NDWI), pH, slope and Enhanced Vegetation Index (EVI) were the background variables in the Maxent models. The models were validated using the area under the curve (AUC) and omission rate.
The predicted suitable habitats for intermediate snail hosts varied with seasons. The AUC for models in all seasons ranged from 0.71 to 1 and the prediction rates were between 0.8 and 0.9. Although B. globosus was found at more localities in the Ndumo area, there was also evidence of cohabiting with B. pfiefferi at some of the locations. NDWI had significant contribution to the models in all seasons.
The Maxent model is robust in snail habitat suitability modelling even with small dataset of presence-only sampling sites. Application of the methods and design used in this study may be useful in developing a control and management programme for schistosomiasis in the Ndumo area.
Schistosomiasis is a snail-borne disease prevalent in humans . The disease ranks second to malaria in terms of the negative socio-economic effects it has in endemic communities [2–4] mostly in the rural parts of sub-Saharan Africa . The spatial and temporal distribution of intermediate snail hosts that transmit the disease determines the distribution of the disease in endemic areas. People living in rural or semi-rural communities are in constant contact with schistosome-infested water . In South Africa, particularly in KwaZulu-Natal Province both urinary (Schistosoma haematobium) and intestinal (S. mansoni) schistosomiasis are endemic indicating the presence of the intermediate snail hosts [7–9]. The prevalence of urinary schistosomiasis has been reported to be high (68–80 %) in uMkhanyakude district [10, 11] including Ndumo area. However, spatial and temporal modelling of the snail habitats mainly at a micro-geographical scale to get insights of schistosomiasis transmission dynamics has not been done.
It is essential that the distribution of schistosome intermediate hosts be known for effective design, implementation and evaluation of schistosomiasis control programs . The population dynamics of the intermediate snail hosts and parasite transmission patterns in situ are still to be adequately studied . Even with inexpensive and effective anthelminthic medication (Praziquantel) the knowledge of the local dynamics of snail populations is still important for timing of mass drug administration, the global strategy endorsed by the World Health Organization . This will help to follow the periods when re-infection is very low [14–17] so as to reduce the chances of reinfection. The sustainability of this control strategy has been challenged, as there is rapid re-infection after deworming [3, 18]. There has been a shift from morbidity control to transmission control and local elimination . Hence, there is a stronger focus on intermediate snail hosts and transmission sites, along with primary prevention tailored to specific socio-ecological systems [16, 17, 20]. Understanding the dynamics of transmission of schistosomiasis could help to identify hot spots where transmission may be intense and build towards effective local intervention programmes.
Specific habitat requirements of intermediate snail hosts are governed by environmental factors . The intermediate snail hosts need an aquatic environment and thrive even in small water bodies (SWBs), such as ponds, ditches and other humid areas consisting of open water, aquatic vegetation and/or inundated grass . Although snails may reproduce through selfing and aestivate during the dry season  triggered by the drying of water pools, live snails are limited to locations with standing water or with enough moisture for survival .
Habitats of intermediate snail hosts can be mapped by extensive ground surveys requiring considerable amounts of time, manpower and money. Thus the use of remotely sensed imagery is a useful alternative for habitat detection as it reduces costs, time and manpower . The use of satellite remote sensing data and techniques for risk profiling of environment-related diseases, including schistosomiasis, has increased considerably over the past 30 years [25, 26]. Remote sensing data have been mainly used to relate schistosomiasis prevalence at the school level to remote sensing measurements such as Normalised Difference Vegetation Index (NDVI) and Normalised Difference Water Index (NDWI), to model and spatially predict the risk of infection [15, 25–29]. However, the aim of the application of remote sensing data is to characterize the environmental conditions of potential disease transmission sites, which are in many cases spatially disjunct from the school location where epidemiological surveys are usually being conducted .
Modelling habitats of intermediate snail hosts using the Maximum Entropy (Maxent) modelling technique could substantially improve our understanding of the temporal and spatial distribution of current risk of schistosomiasis and create novel possibilities for improved schistosomiasis control and management [22, 25, 30–32]. Maxent is a recently developed ecological modelling method capable of achieving high predictive performance  using the presence-only data [33–35] in contrast to the background environmental conditions. For planning of successful interventions against schistosomiasis and to target populations living in high risk areas, it is of great importance to determine the current spatial distribution of infection at a reasonably fine scale, including the distribution of parasites and host species . Detailed maps of possible distribution of habitats of intermediate snail hosts provide valuable information for the prediction of infection risk zones but are currently lacking for most parts of the world  There is need for micro-geographical studies on the spatial and temporal distribution of these species to guide the control and management of schistosomiasis especially at the community level. Therefore, the purpose of this paper was to model the spatial and seasonal distribution of suitable habitats of intermediate snail hosts of Schistosoma spp. based on climatic and non-climatic factors using Maxent in the Ndumo area, uMkhanyakude district in the KwaZulu-Natal province of South Africa. This method has only been used by Stensgaard  in modelling schistosomiasis/snails in Africa at the continental level and Pedersen  at the national level.
Ndumo area is located in the uMkhanyakude Health District in the KwaZulu-Natal (KZN) province, South Africa (Fig. 1). uMkhanyakude is located in the northenmost eastern part of the KwaZulu-Natal province bordering Mozambique and Swaziland to the north and north-west, respectively. The area is approximately 40 × 30 km mostly characterised by seasonal streams flowing towards the Pongola flood plain. There are two main dams in this area: Nsunduza and Namaneni. The climate of the area ranges from tropical to subtropical  experiencing low precipitation averaging at 690 mm per year. The year was divided into four seasons according to temperature and rainfall; rainy (December to February), post-rainy (March to May), cold/dry (June to August) and hot/dry (September to November) based on previous studies in the same region [39–41]. The presence of schistosome intermediate snail hosts (B. globosus and B. pfeifferi) in uMkhanyakude has been reported by Appleton , De Kock et al.  and De Kock & Wolmarans .
Snail survey data
The snail sampling points were fairly distributed in the study area (Fig. 1) and were located along rivers, streams and dams targeting suspected transmission sites as advised by Appleton & Miranda . A total of 16 sites were sampled for a period of 1 year (May 2014 to April 2015). The number of sites sampled was considered adequate as the main species modelling technique (Maxent) used in this study has been proven to perform well with a minimum of 10 sampling points  and Pearson et al.  supported the use of Maxent when sample sizes are very low (<10 but ≥ 5). Maxent also requires presence-only data [33, 34] hence only sites where snails were present in a particular season were used in Maxent seasonal models (Table 1). During the period of sampling B. globosus was found in 12 sites and at 7 of those sites both B. globosus and B. pfeifferi were found. Four sites including S3, S5, S11 and S15 (Fig. 1) had neither B. globosus nor B. pfeifferi.
Climatic and environmental factors
When considering relevant candidate variables to be used in a model, the ecology of the species in question needs to be taken into consideration, hence the climatic and non-climatic predictor variables were selected based on their perceived biological relevance for host snail distributions [36, 37, 45–47]. Specifics and sources of the climatic and environmental variables used in this study are listed in Table 2. We focussed on 2 overall classes of environmental variables that have been shown to influence host snail distribution patterns [48, 49] namely (i) climatic variables (temperature and precipitation) and (ii) natural habitat conditions (water bodies and soil conditions) . The survival and reproduction rates of snails in relation to temperature have been described in a number of studies [47, 48, 50–53]. The diurnal temperature range was chosen to account for the demonstrated importance of fluctuating temperatures, as previously shown for B. pfeifferi . Temperature of the warmest and coldest quarter was considered to account for the sensitivity of snails to temperature extremes [53, 55]. Seasonal precipitation was used as a measure of the availability of suitable temporary water bodies that snails are known to inhabit . Soil pH was considered as it has been shown to influence pH in water bodies [36, 56, 57]. A 3 month average (March-May) for maximum temperature (Tmax), minimum temperature (Tmin), Normalised Difference Vegetation Index (NDVI), Normalised Difference Water Index (NDWI), Enhanced Vegetation Index (EVI) were used as they have been used by Pedersen et al. . In this study we considered these factors at 2 levels: annual averages and seasonal averages to model the variation of suitable habitats of the host snails at micro-geographical scale.
The Normalised Difference Water Index (NDWI) was calculated based on Moderate Resolution Imaging Spectroradiometer (MODIS) reflectance and temperature amplitude was calculated from MODIS minimum and maximum temperature. Climate Hazards Group InfraRed Precipitation with Station data (CHIRPS) was recorded for 30 years (January 1981 to August 2014) and other datasets from MODIS started between 2000 and 2003 − 2014. We used these datasets as they have higher spatial resolution; most of the datasets were accessed through the International Research Institute for Climate and Society (IRI) data library portal (http://iridl.ldeo.columbia.edu/SOURCES/).
We performed a descriptive statistical analysis (means and standard deviations, SD) of climatic and environmental variables to determine their spatial variation annually and in different seasons. We also removed highly correlated variables (>0.9)  as the multicollinearity may violate statistical assumptions and may alter model predictions . NDVI was highly correlated with other variables hence was removed from the models in other seasons.
Predictive model implementation
Bulinus globosus and Biomphalaria pfeifferi species distribution models were developed using the Maximum Entropy (Maxent) [33–35] approach using the species distribution models (SDM) toolbox developed by Brown et al.  and implemented in ArcGIS 10.2.
Model performance was expressed as the area under (the receiver operator characteristic) curve (AUC)  supported by sensitivity and specificity . An AUC value of 0.5 indicates that the model predicts no better than a random model, while AUC values of > 0.75 are considered in the “best” model category . However, comparing models across species using AUC scores is problematic, as AUC is influenced by species’ prevalence . This issue was alleviated by only comparing AUC values among models within species . The spatial jackknifing or geographically structured k-fold validation which test evaluation performance of spatially segregated localities was also used in this study [60, 63]. Hence we considered the omission rate and model feature class complexity instead of the specificity and sensitivity since we used small sample sizes (<25) which could have inflated the performance of the model. The spatial jackknifing script in the SDM toolbox chooses the best model by evaluating each model's omission rates (OR), AUC and model feature class complexity. A jackknife procedure, implemented in Maxent, was used to quantify the explanatory power of each environmental variable. “Maximum training sensitivity plus specificity statistics” output from Maxent as a threshold criterion was used in partitioning of the observations into suitable and unsuitable habitats following the recommendation by Hu & Jiang . We then calculated the area for the suitable habitats and also expressed it as a percentage of the total habitat for each model to quantify the differences in habitats in different seasons and between B. globosus and B. pfeifferi.
Variation of climatic and environmental factors
The distribution of habitats of intermediate snail hosts is influenced by the variation or changes in environmental and climatic factors. Tables 3 and 4 indicate the spatial variability of these factors annually and seasonally, respectively in the Ndumo area of uMkhanyakude district in South Africa.
Spatial distribution of schistosome host snail habitats
The predicted distribution of B. globosus and B. pfeifferi habitats suitability based on annual averages of climatic variables and non-climatic variables (shown in Table 3) is depicted in Fig. 2. The model predictions indicate that Ndumo area may have more locations with relatively higher probabilities of finding suitable habitats for B. globosus compared to B. pfeifferi. Figure 3 shows the predicted snail habitats suitable for the two snail species (B. globosus and B. pfeifferi) in different seasons in the Ndumo area. The possibility of cohabiting of the two species and focality in suitable habitats is also depicted in the two model outputs (Figs. 2 and 3).
The habitat threshold values used to determine the suitable and unsuitable areas were different for each model (Table 5). Figure 4 shows the variation in suitable areas for the host snails in different seasons classified based on maximum training sensitivity plus specificity from Maxent model.
The AUC values ranged from 0.71 to 1 for B. globosus models and 0.80 to 0.87 for B. pfeifferi respectively, indicating good model performance. The weighted prediction for the best performing models was high, ranging from 0.80 to 0.89 for B. globosus and 0.89 to 0.90 for B. pfeifferi models. The performance of Maxent models in modelling the spatial distribution of B. globosus and B. pfeifferi and the related factors as well as the variation in the suitable habitats is shown in the Table 5.
NDWI had a consistent high contribution to all of the models in Table 5 regardless of the season. Slope only showed significant contribution on B. globosus annual and hot/dry season models. Maximum temperature had a higher contribution to both B. globosus and B. pfeifferi cold/dry season models. Temperature amplitude had higher significance in B. globosus and annual and hot/dry season models. Bulinus globosus had more suitable areas in the cold/dry season compared to B. pfeifferi (Table 5). However, the two species shared most of the localities. Using annual averages we found that there are more areas suitable for B. globosus compared to B. pfeifferi (Table 5). Bulinus globosus showed the highest percentage of suitable area in the cold/dry season compared to post-rainy and hot/dry seasons, with lowest values in the hot/dry season. This is explained by the tolerance of these snail species (B. globosus and B. pfeifferi) to the variations in climatic and environmental factors (Tables 3 and 4).
The objective of this study was to model and predict the distribution of suitable habitats for B. globosus and B. pfeifferi at the micro-geographical scale. Maxent models [33–35] do not predict the actual limits of a species’ range but can identify localities with similar conditions for occurrence . Models presented in our study have indicated a good estimation of the distribution of suitable habitats of schistosome intermediate snail hosts at a micro-geographical scale based on climatic and environmental factors. The results of this study have supported the well-known establishment that Maxent is an efficient tool in modelling species distribution when the datasets are small [43, 44, 65]. In our case we had other seasons with very low numbers of sampling points (<5) and were considered not suitable for obtaining meaningful results [44, 60]. According to the model performance evaluation criteria by Phillips & Dudík,  the Maxent model used in this study performed satisfactorily as indicated by the high AUC values which exceeded 0.75 except in 1 case where it was 0.71. However the spatial jackknifing method  used in the model evaluation in this study considered omission rate first as the key measure of model performance hence the same model had higher prediction rate (0.8). The Maxent model performance was also appraised by Pedersen et al.  and Stensgaard et al.  on intermediate snail host habitat suitability modelling. The two species have different foci but overlap at most of the locations in our study area. Similar observations were made by Pedersen et al.  based on a national scale study in Zimbabwe. The suitable habitats contract and expand as determined by the environmental and climatic factors as the cold/dry season models showed a wider niche for B. globosus compared to post-rainy and hot/dry seasons. However, there are still important conceptual uncertainties in these models which need to be investigated, especially identification of causal relationships between species distribution and predictors . The small sample size and the jacknifing validation method used in this study may inflate the accuracy of the models, hence there is need to assess the performance of this model in a different setting.
We also observed that snail presence probability varied by locations indicating differences in location suitability for the two species, B. globous and B. pfeifferi. The two snail species showed different levels of sensitivity to different climatic and environmental factors in terms of their suitable habitats. However, NDWI which is an indication of surface water, was the most consistent and significant variable in both species models in all seasons. There are large areas in the eastern parts of our study area where there are water bodies but do not seem to be places where lots of snails are expected. This might be due to factors other than surface water (as detected with NDWI) such as the slope and temperatures which do not favour the presence of host snails. The use of “presence only” data may exclude from further consideration certain habitat types that are deemed falsely unsuitable thereby limiting the creation of models that accurately discriminate between suitable and unsuitable habitats . Our results are slightly different from those of Pedersen et al.  who conducted a similar study at the national scale in Zimbabwe which is experiencing almost similar climatic conditions as our study area. Their study showed that B. globosus distribution was more influenced by the maximum temperature in the post-rain season (March to May) while B. pfeifferi was more influenced by the minimum temperature in that same season. The difference is mainly because Pedersen et al.  only considered these variables from 1 season unlike in this current study where we considered different seasons. pH did not show any contribution to the spatial distribution of the host snails. This is contrary to observations by Pedersen et al.  where pH showed a significant contribution to the presence of these species. Our study is mainly local and may have lacked significant spatial variation in terms of soil and vegetation types which determine the spatial variation in pH.
Precipitation indirectly affects snails as it indicates the probability of an aquatic environment. However, this variable did not show higher significance levels in the models in our study. Precipitation is related to NDVI, EVI and NDWI but they did not show high correlations. After the rainy season, there are higher chances of finding snails in their preferred habitats since the water bodies are reduced in volume and snail densities are at their peak with many snails having grown to full size . Since suitable habitats mainly for B. globosus have shown a seasonal peak (cold/dry season in our study) and is strongly influenced by NDWI which is a measure of surface water, the inclusion of multi-temporal classification of remote sensing images for surface water detection as noted by De Roeck et al.  could improve spatial distribution model outputs. The determined relationships may be used to predict possible spatial and temporal changes or variation in snail habitats and snail densities over the past and future projections complementary to Maxent models. It is unfortunate that we did not have enough sampling sites to run Maxent models for all the seasons for both species as some sites became dry during the dry season. In this study the rainy season (December to February) had too few positive data points for further Maxent analysis and it was anticipated that there could be more suitable areas based on increased rainfall which may increase surface water which is critical for snail habitats. However, Anderson , noted that it is difficult to measure the spatial relationship between rainfall and snail population dynamics and infection transmission since the effect of rainfall varies depending on the species of snail and the geographical location. Although snails do not thrive without water as they tend to aestivate and may not be detected easily for such studies, too much water may also reduce snail populations  especially in fast flowing rivers or streams. However these results give insight into the spatial and temporal dynamics of the suitable habitats for intermediate snail hosts which is critical for control, monitoring and management of the disease.
Although the presence of snails is necessary to determine an area of schistosomiasis transmission the mere presence of snails is not sufficient to define a transmission site. If the snails are not infected and there is no human-water contact, transmission does not occur. Thus, while useful, this model alone may not be adequate for developing a schistosomiasis control, monitoring and management scheme. Notwithstanding this, our model was able to confirm observations made by Moodley et al.  and Pitchford  based on minimum temperature suitable for schistosomiasis transmission; and by Manyangadze et al.  and Saathoff et al.  through parasitological surveys, mainly for S. haematobium. Thus our model complements the efforts made using different methodologies to understand the dynamics of schistosomiasis in this area. Combining the spatial distribution models of schistosomiasis based on environmental and socio-economic factors  and the current model (distribution of potential habitats for schistosome intermediate snail hosts) may help to develop effective schistosomiasis control strategies. This could also be supported by the data on the distribution of infected snails and human-water contact sites.
This study intended to model the spatial distribution of the suitable habitats for B. globosus and B. pfeifferi at a micro-geographical scale. The method (Maxent) used is robust in modelling suitable habitats for the host snails even with a small dataset of presence-only sampling sites. The results showed that suitable habitats of the schistosome intermediate snail hosts B. globosus and B. pfiefferi may vary at the micro-geographical scale. Although B. globosus was found at more localities in the Ndumo area, there was also evidence of cohabiting with B. pfiefferi at some of the locations. NDWI, which is a proxy for surface water was more significant and consistent in its contribution to the models in all seasons. The methods and design used in this study give informative results which may help in control, monitoring and management of schistosomiasis in the area.
Area under (the receiver operator characteristic) curve
Climate Hazards Group InfraRed Precipitation with Station data
Enhanced Vegetation Index
International Research Institute for Climate and Society
International Soil Reference and Information Centre - World Inventory of Soil Emission Potentials
Moderate resolution imaging spectroradiometer
National aeronautics and space administration
Normalised difference vegetation index
Normalised difference water index
Species distribution model
University of California Santa Barbara
United States geological survey
Moser W, Greter H, Schindler C, Allan F, Ngandolo BNR, Moto DD, et al. The spatial and seasonal distribution of Bulinus truncatus, Bulinus forskalii and Biomphalaria pfeifferi, the intermediate host snails of schistosomiasis, in N’Djamena, Chad. Geospat Health. 2014;9(1):109–18.
Sleigh A, Li X, Jackson S, Huang K. Eradication of schistosomiasis in Guangxi, China. Part 3. Community diagnosis of the worst-affected areas and maintenance strategies for the future. Bull World Health Org. 1998;76(6):581.
Utzinger J, Bergquist R, Shu-Hua X, Singer BH, Tanner M. Sustainable schistosomiasis control - the way forward. Lancet. 2003;362(9399):1932–4.
Kazibwe F, Makanga B, Rubaire-Akiiki C, Ouma J, Kariuki C, Kabatereine NB, et al. Ecology of Biomphalaria (Gastropoda: Planorbidae) in Lake Albert, Western Uganda: snail distributions, infection with schistosomes and temporal associations with environmental dynamics. Hydrobiologia. 2006;568(1):433–44.
World Health Organisation. WHO schistosomiasis fact sheet;2014; 2015.
Stothard JR, Kabatereine NB, Tukahebwa EM, Kazibwe F, Mathieson W, Webster JP, et al. Field evaluation of the Meade Readiview handheld microscope for diagnosis of intestinal schistosomiasis in Ugandan school children. T Roy Soc Trop Med H. 2005;73(5):949–55.
De Kock K, Wolmarans C. Distribution and habitats of the Bulinus africanus species group, snail intermediate hosts of Schistosoma haematobium and S. mattheei in South Africa. Water SA. 2005;31(1):117–25.
De Kock K, Wolmarans C, Bornman M. Distribution and habitats of Biomphalaria pfeifferi, snail intermediate host of Schistosoma mansoni, in South Africa. Water SA. 2004;30(1):29–36.
Appleton CC. Freshwater molluscs of Southern Africa with a chapter on Bilharzia and its snail hosts. Durban: University of KwaZulu Natal Press; 1996.
Saathoff E, Olsen A, Magnussen P, Kvalsvig JD, Becker W, Appleton CC. Patterns of Schistosoma haematobium infection, impact of praziquantel treatment and re-infection after treatment in a cohort of schoolchildren from rural KwaZulu-Natal/South Africa. BMC Infect Dis. 2004;4:40.
Lankford B, Pringle C, Dickens C, Lewis F, Chhotray V, Mander M, et al. The impacts of ecosystem services and environmental governance on human well-being in the Pongola region, South Africa. Report to the Natural Environment Research Council. London, Norwich, and Pietermaritzburg: University of East Anglia and Institute of Natural Resources; 2010.
Kloos H, Souza C, Gazzinelli A, Soares Filho BS, Temba PC, Bethony J, et al. The distribution of Biomphalaria spp. in different habitats in relation to physical, biological, water contact and cognitive factors in a rural area in Minas Gerais, Brazil. Mem Inst Oswaldo Cruz. 2001;96:57–66.
WHO. Prevention and Control of Schistosomiasis and Soil-Transmitted Helminthiasis: Report of a WHO Expert Committee. Geneva: World Health Organization; 2002.
Stothard JR, Mgeni AF, Khamis S, Seto E, Ramsan M, Rollinson D. Urinary schistosomiasis in schoolchildren on Zanzibar Island (Unguja), Tanzania: a parasitological survey supplemented with questionnaires. T Roy Soc Trop Med H. 2002;96(5):507–14.
Brooker S. Spatial epidemiology of human schistosomiasis in Africa: risk models, transmission dynamics and control. T Roy Soc Trop Med H. 2007;101(1):1–8.
Stothard JR, Chitsulo L, Kristensen TK, Utzinger J. Control of schistosomiasis in sub-Saharan Africa: progress made, new opportunities and remaining challenges. Parasitology. 2009;136(13):1665–75.
Rollinson D, Knopp S, Levitz S, Stothard JR, Tchuem Tchuenté LA, Garba A, et al. Time to set the agenda for schistosomiasis elimination. Acta Trop. 2013;128(2):423–40.
Clements A, Bosqué-Oliva E, Sacko M, Landouré A, Dembélé R, Traoré M, et al. A comparative study of the spatial distribution of schistosomiasis in Mali in 1984–1989 and 2004–2006. PLoS Negl Trop Dis. 2009;3(5):e431.
Walz Y, Wegmann M, Dech S, Vounatsou P, Poda JN, N'Goran EK, et al. Modeling and validation of environmental suitability for schistosomiasis transmission using remote sensing. PLoS Negl Trop Dis. 2015;9(11):e0004217. doi:10.1371/journal.pntd.0004217.
Utzinger J, N’Goran EK, Caffrey CR, Keiser J. From innovation to application: social-ecological context, diagnostics, drugs and integrated control of schistosomiasis. Acta Trop. 2011;120:S121–37.
Brown D. Freshwater snails of Africa and their medical importance. London: Taylor and Francis; 1994.
De Roeck E, Van Coillie F, De Wulf R, Soenen K, Charlier J, Vercruysse J, et al. Fine-scale mapping of vector habitats using very high resolution satellite imagery: a liver fluke case-study. Geospat Health. 2014;8(3):671–83.
McCullough F, Eyakuze VM, Msinde J, Nditi H. Water resources and bilharziasis transmission in the Misungwi area, Mwanza District, north-west Tanzania. E Afr Med J. 1968;45(5):295–308.
Clennon JA, King CH, Muchiri EM, Kitron U. Hydrological modelling of snail dispersal patterns in Msambweni, Kenya and potential resurgence of Schistosoma haematobium transmission. Parasitology. 2007;134(5):683–93.
Simoonga C, Kazembe LN, Kristensen TK, Olsen A, Appleton CC, Mubita P, et al. Remote sensing, geographical information system and spatial analysis for schistosomiasis epidemiology and ecology in Africa. Parasitology. 2009;136(13):1683–93.
Manyangadze T, Chimbari MJ, Gebreslasie M, Mukaratirwa S. Application of geo-spatial technology in schistosomiasis modelling in Africa: a review. Geospat Health. 2015;10(326):99–110.
Seto E, Xu B, Liang S, Gong P, Wu W, Davis G, et al. The use of remote sensing for predictive modeling of schistosomiasis in China. Photogramm Eng Rem S. 2002;68(2):167–74.
Stensgaard A, Jorgensen A, Kabatareine NB, Malone JB, Kristensen TK. Modeling the distribution of Schistosoma mansoni and host snails in Uganda using satellite sensor data and Geographical Information Systems. Parasitologia. 2005;47(1):115.
Walz Y, Wegmann M, Dech S, Raso G, Utzinger J. Risk profiling of schistosomiasis using remote sensing: approaches, challenges and outlook. Parasit Vectors. 2015;8(1):163.
Lacaux J, Tourre YM, Vignolles C, Ndione JA, Lafaye M. Classification of ponds from high-spatial resolution remote sensing: Application to Rift Valley Fever epidemics in Senegal. Remote Sens Environ. 2007;106(1):66–74.
Charlier J, Bennema SC, Caron Y, Counotte M, Ducheyne E, Hendrickx G, et al. Towards assessing fine-scale indicators for the spatial transmission risk of Fasciola hepatica in cattle. Geospat Health. 2011;5(2):239–45.
Estallo EL, Lamfri MA, Scavuzzo CM, Almeida FFL, Introini MV, Zaidenberg M, et al. Models for predicting Aedes aegypti larval indices based on satellite images and climatic variables. J Am Mosquito Contr. 2008;24(3):368–76.
Phillips SJ, Dudík M. Modeling of species distributions with Maxent: new extensions and a comprehensive evaluation. Ecography. 2008;31(2):161–75.
Phillips SJ, Anderson RP, Schapire RE. Maximum entropy modeling of species geographic distributions. Ecol Model. 2006;190(3):231–59.
Elith J, Phillips SJ, Hastie T, Dudík M, Chee YE, Yates CJ. A statistical explanation of MaxEnt for ecologists. Divers Distrib. 2011;17(1):43–57.
Pedersen UB, Midzi N, Mduluza T, Soko W, Stensgaard AS, Vennervald BJ, et al. Modelling spatial distribution of snails transmitting parasitic worms with importance to human and animal health and analysis of distributional changes in relation to climate. Geospat Health. 2014;8(2):335–51.
Stensgaard AS, Utzinger J, Vounatsou P, Hurlimann E, Schur N, Saarnak CF, et al. Large-scale determinants of intestinal schistosomiasis and intermediate host snail distribution across Africa: does climate matter? Acta Trop. 2013;128(2):378–90.
Schulze B. Climate of South Africa: General Survey by BR Schulze. Government Printer: South Africa; 1965.
Pfukenyi D, Mukaratirwa S, Willingham AL, Monrad J. Epidemiological studies of Fasciola gigantica infections in cattle in the highveld and lowveld communal grazing areas of Zimbabwe. Onderstepoort J Vet. 2006;73(1):37–51.
Chingwena G, Mukaratirwa S, Chimbari M, Kristensen TK, Madsen H. Population dynamics and ecology of freshwater gastropods in the highveld and lowveld regions of Zimbabwe, with emphasis on schistosome and amphistome intermediate hosts. Afr Zool. 2004;39(1):55–62.
Chandiwana S, Christensen N, Frandsen F. Seasonal patterns in the transmission of Schistosoma haematobium, S. mattheei and S. mansoni in the highveld region of Zimbabwe. Acta Trop. 1987;44(4):433–44.
Appleton CC, Miranda N. Locating bilharzia transmission sites in South Africa: guidelines for public health personnel. S Afr J Infect Dis. 2015;30(3):95–102.
Wisz MS, Hijmans RJ, Li J, Peterson AT, Graham CH, Guisan A. Effects of sample size on the performance of species distribution models. Divers Distrib. 2008;14(5):763–73.
Pearson RG, Raxworthy CJ, Nakamura M, Townsend PA. Predicting species distributions from small numbers of occurrence records: a test case using cryptic geckos in Madagascar. J Biogeogr. 2007;34(1):102–17.
Malone J, Yilma JM, McCarroll JC, Erko B, Mukaratirwa S, Zhou X. Satellite climatology and the environmental risk of Schistosoma mansoni in Ethiopia and east Africa. Acta Trop. 2001;79(1):59–72.
Stensgaard A-S, Jørgensen A, Kabatereine NB, Rahbek C, Kristensen TK. Modeling freshwater snail habitat suitability and areas of potential snail-borne disease transmission in Uganda. Geospat Health. 2006;1(1):93–104.
Zhou XN, Yang GJ, Yang K, Wang XH, Hong QB, Sun LP, et al. Potential impact of climate change on schistosomiasis transmission in China. T Roy Soc Trop Med H. 2008;78(2):188–94.
Appleton CC. Review of literature on abiotic factors influencing the distribution and life cycles of bilharziasis intermediate host snails. Malacol Rev. 1978;11:1–25.
Brooker S, Michael E. The potential of geographical information systems and remote sensing in the epidemiology and control of human helminth infections. Adv Parasitol. 2000;47:245–88.
Woolhouse M, Chandiwana S. Population dynamics model for Bulinus globosus, intermediate host for Schistosoma haematobium, in river habitats. Acta Trop. 1990;47(3):151–60.
Woolhouse M, Chandiwana S. Population biology of the freshwater snail Bulinus globosus in the Zimbabwe highveld. J Appl Ecol. 1990;27:41–59.
Mas-Coma S, Valero MA, Bargues MD. Climate change effects on trematodiases, with emphasis on zoonotic fascioliasis and schistosomiasis. Vet Parasitol. 2009;163(4):264–80.
Pflüger W. Experimental epidemiology of schistosomiasis. Z Parasitenkd. 1981;66(2):221–9.
Appleton CC, Eriksson IM. The influence of fluctuating above-optimal temperature regimes on the fecundity of Biomphalaria pfeifferi (Mollusca: Planorbidae). T Roy Soc Trop Med H. 1984;78(1):49–54.
Appleton CC. The influence of temperature on the life-cycle and distribution of Biomphalaria pfeifferi (Krauss, 1948) in South-Eastern Africa. Int J Parasitol. 1977;7(5):335–45.
Bavia ME, Hale LF, Malone JB, Braud DH, Shane SM. Geographic information systems and the environmental risk of schistosomiasis in Bahia, Brazil. T Roy Soc Trop Med H. 1999;60(4):566–72.
Kloos H, Passos LKJ, LoVerde P, Oliveira RC, Gazzinelli A. Distribution and Schistosoma mansoni infection of Biomphalaria glabrata in different habitats in a rural area in the Jequitinhonha Valley, Minas Gerais, Brazil: environmental and epidemiological aspects. Mem I Oswaldo Cruz. 2004;99(7):673–81.
Fourcade Y, Engler JO, Rödder D, Secondi J. Mapping species distributions with MAXENT using a geographically biased sample of presence data: a performance assessment of methods for correcting sampling bias. PLoS ONE. 2014;9(5):e97122.
Heikkinen RK, Luoto M, Araújo MB, Virkkala R, Thuiller W, Sykes MT. Methods and uncertainties in bioclimatic envelope modelling under climate change. Prog Phys Geog. 2006;30(6):751–77.
Brown JL. SDMtoolbox: a python‐based GIS toolkit for landscape genetic, biogeographic and species distribution model analyses. Methods Ecol Evol. 2014;5(7):694–700.
Lobo JM, Jiménez‐Valverde A, Real R. AUC: a misleading measure of the performance of predictive distribution models. Global Ecol Biogeography. 2008;17(2):145–51.
Liu C, Berry PM, Dawson TP, Pearson RG. Selecting thresholds of occurrence in the prediction of species distributions. Ecography. 2005;28(3):385–93.
Shcheglovitova M, Anderson RP. Estimating optimal complexity for ecological niche models: A jackknife approach for species with small sample sizes. Ecol Model. 2013;269:9–17.
Hu J, Jiang Z. Climate change hastens the conservation urgency of an endangered ungulate. PLoS ONE. 2011;6(8):e22873.
Suárez-Seoane S, de la Morena ELG, Prieto MBM, Osborne PE, de Juana E. Maximum entropy niche-based modelling of seasonal changes in little bustard (Tetrax tetrax) distribution. Ecol Model. 2008;219(1):17–29.
Barry S, Elith J. Error and uncertainty in habitat models. J Appl Ecol. 2006;43(3):413–23.
Anderson R. Determinants of infection in human schistosomiasis. Baillieres Clin Trop Med Comm Dis. 1987;2(2):279–300.
Moodley I, Kleinschmidt I, Sharp B, Craig M, Appleton CC. Temperature-suitability maps for schistosomiasis in South Africa. Ann Trop Med Parasit. 2003;97(6):617–27.
Pitchford R. Temperature and schistosome distribution in South-Africa. S Afr J Sci. 1981;77(6):252–61.
Manyangadze T, Chimbari MJ, Gebreslasie M, Mukaratirwa S. Risk factors and micro-geographical heterogeneity of Schistosoma haematobium in Ndumo area, uMkhanyakude district, KwaZulu-Natal, South Africa. Acta Trop. 2016;159:176–84.
The authors would like to express their appreciation to Malaria and Bilharzia in Southern Africa (MABISA) team members mainly in South Africa for field assistance.
The research has been supported by College of Health Sciences scholarship program at University of KwaZulu-Natal and Malaria and Bilharzia in Southern Africa (MABISA) project funded by WHO Special Programme for Research and Training in Tropical Diseases (TDR) and the Canadian International Development Research Centre (IDRC) which facilitated the field data collection and analysis and the contribution from the NASA SERVIR project: NNX12AQ70G for facilitating access to satellite data.
Availability of data and material
The data supporting the conclusions of this article are included within the article.
TM, MJC, SM and MG conceptualised the study. TM collected the data and worked on data access and analysis together with PC and MG. MJC and SM contributed on snail ecology. TM drafted the manuscript and all authors read and approved the final manuscript.
The authors declare that they have no competing interests.
Consent for publication
Ethics approval and consent to participate