Environmental and biotic factors affecting freshwater snail intermediate hosts in the Ethiopian Rift Valley region
Parasites & Vectors volume 13, Article number: 292 (2020)
Knowledge of the distribution and habitat preference of freshwater snail intermediate hosts can provide information to initiate and set-up effective snail control programmes. However, there is limited research conducted on the factors driving the occurrence and abundance of freshwater snail intermediate hosts in the Ethiopian Rift Valley. Hence, in this study, we investigated how environmental and biotic factors influence the occurrence and abundance of the snail intermediate hosts in Ethiopian Rift Valley region.
Data on freshwater snails, physico-chemical water quality parameters, physical characteristics of habitat, predators and competitors, and anthropogenic activity variables were collected from 174 sampling sites during the wet season of 2017 and 2018. Generalized linear models were used to identify the main environmental and biotic factors affecting the occurrence and abundance of the snail species.
It was found that Bulinus globosus (31.7%) was the most abundant snail species followed by Lymnaea natalensis (21.6%), Lymnaea truncatula (15.1%) and Biomphalaria pfeifferi (14.6%). Generalized linear models indicated that physico-chemical parameters (water temperature, turbidity, chlorophyll-a, dissolved oxygen, chemical oxygen demand, alkalinity, calcium, magnesium, nitrate and ammonia), physical habitat characteristics (water depth, canopy cover, macrophyte cover and substrate type) and biotic factors (abundance of predators and competitors) were found to be the main variables determining the occurrence and abundance of snail species in the Ethiopian Rift Valley region. In terms of anthropogenic activities, human settlement, farming, bathing and swimming, clothes washing, grazing, drainage of land, car washing, boating, fishing and silviculture were also important variables determining the occurrence and abundance of snail species in the region.
The findings reported herein suggest that integrated snail control strategies should be considered to control snails via protection of water bodies from disturbance by anthropogenic activities. In this way, it is possible to reduce the concentration of organic matter and dissolved ions in aquatic ecosystems which are conducive for the presence of snails.
Snails are invertebrate animals of the class Gastropoda found in freshwater and other aquatic habitats around the world . Approximately 5000 species of snails are found to inhabit different habitats worldwide . Freshwater snails receive considerable attention as they are intermediate hosts of several trematodes that can cause diseases in humans and animals . Among snail-borne diseases, schistosomiasis causes a major public health concern. It is estimated that 779 million people are at risk of schistosomiasis , with 250 million people at risk of infection, mainly in sub-Saharan Africa . Schistosomiasis has a widespread distribution in several parts of Ethiopia usually at an altitude between 1200–2000 m above sea level , where an estimated 38.3 million people are living in schistosomiasis endemic areas . In endemic areas, children have a greater risk of becoming infected as well as re-infected [8,9,10]. Fascioliasis is another important parasitic snail-borne disease in tropical and subtropical countries . It has been estimated that 2.4 million people are infected worldwide with fascioliasis while 180 million people are at risk of infection . The diseases are acquired through repeated contact with freshwater during washing, swimming or fishing [4, 5].
Transmission of snail-borne diseases depends on the distribution of specific freshwater snails that act as intermediate hosts and are the first step for a trematode parasite to reach the development stage to infect humans or animals [13, 14]. When parasite eggs are released into freshwater bodies with urine or feces, miracidia hatch and infect the intermediate host snails [15, 16]. In the snails, miracidium develops into a mother sporocyst. In schistosomes, the sporocyst develops into the second generation sporocysts, in which the infective larvae cercariae are formed. In some hermaphroditic trematodes (e.g. liver flukes), the mother sporocyst develops into rediae which produce cercariae [16, 17]. Once the cercariae are released into the water, they either penetrate the skin of the definitive host (e.g. schistosomes) or are ingested after encysting as metacercariae in or on edible plants or animals. After entering the definitive host, the schistosome larvae mature into adult worms in the blood vessels of the liver, intestine and bladder. The worms lay thousands of eggs that cause damage as they grow through tissues and consequently, infection occurs accordingly [17, 18]. Most malacologists studying freshwater snails focus mainly on snail species belonging to the genera Biomphalaria, Bulinus and Lymnaea which are known to serve as intermediate hosts of schistosomes and play a crucial role in the transmission of the disease in tropical and subtropical regions of the world . Other freshwater snails that have no prominent medical or veterinary importance receive less or no attention. Generally, these are herbivores removing vegetation biomass, which may affect the standing crop and distribution of primary producers in an aquatic ecosystem . In Africa, there are several snail species belonging to the genera Biomphalaria, Bulinus and Lymnaea that act as intermediate hosts of trematode parasites and that are of medical and veterinary importance .
A countrywide snail survey made in Ethiopia showed that planorbid (Biomphalaria spp.), bulinid and lymnaeid snails of medical and veterinary importance have a wide geographical distribution [21, 22]. It is recognized that two species of the genus Biomphalaria (B. pfeifferi and B. sudanica) are the sole intermediate hosts transmitting Schistosoma mansoni , whereas two species of the genus Bulinus (Bu. abyssinicus and Bu. africanus) are the intermediate hosts transmitting Schistosoma heamatobium . The distribution of B. sudanica is restricted to the three Rift Valley areas of Ethiopia, along the shores of lakes Ziway and Abaya and the interface between Lake Hawassa and its tributary Tikur Wuha River , whereas B. pfeifferi has a ubiquitous distribution . The distribution of Bu. abyssinicus and Bu. africanus is limited to the lowland areas, including the Awash and Wabe Shebele Valleys and along the Ethiopian-Sudan border . Lymnaea natalensis and L. truncatula are the two snail species belonging to the genus Lymnaea found to transmit Fasciola parasites causing fascioliasis in humans [25, 26] and ruminant animals in Ethiopia [27, 28] and other African countries [29, 30]. Identifying factors that influence the distribution and habitat preference of the snail intermediate hosts is critical to snail-borne disease prevention and control efforts. In Africa, several ecological studies have shown that biotic and abiotic factors affect the distribution and habitat preference of freshwater snails [31,32,33,34,35]. Similarly, in Ethiopia, a number of studies have been carried out on the ecology of freshwater snails. These studies suggested that local environmental and biotic factors determine the occurrence and abundance of freshwater snails in each aquatic habitat due to environmental heterogeneity among ecological zones [22, 23, 35]. As in many other geographical regions in the country, the endemicity of snail-borne diseases has long been established in the Ethiopian Rift Valley [6, 36, 37]. However, there has been little research conducted on the factors driving the occurrence and abundance of freshwater snail intermediate hosts in the Ethiopian Rift Valley. Therefore, this study aimed to (i) determine the local distribution and diversity of freshwater snail intermediate hosts and (ii) examine the influence of environmental and biotic factors on their occurrence and abundance in the Ethiopian Rift Valley region. The findings of this study could be useful for the development of appropriate preventive and control measures against snail intermediate hosts in the Ethiopian Rift Valley.
Study area and location of sampling sites
Data was collected from two spatially distinct agro-ecological zones, warm temperate rainy and tropical rainy climate zones, which are located in the Ethiopian Rift Valley (Fig.1). These two studied zones are 125 km apart. In the warm temperate rainy climate zone, sampling was carried out along the shores of Lake Ziway which is situated between latitudes 07°51″N and 08°07″N and longitudes 38°43″E and 38°50″E at an elevation of 1636 m above sea level . The mean annual rainfall and temperature of the area are 837 mm and 19.3 °C, respectively . In the tropical rainy climate zone, sampling was carried out along the shores of Lake Hawassa and in Shallo wetland where the Tikur Wuha River drains before it joins Lake Hawassa. Lake Hawassa is situated between latitudes 06°58″N and 07°14″N and longitudes 38°22″E and 38°28″E at an elevation of 1685 m above sea level . The area receives a mean annual rainfall of 950 mm and has a mean annual temperature of 19.8 °C . The aquatic habitats surveyed from both agro-ecological zones are surrounded by local communities, and they serve as an important source of water for surrounding rural communities for various purposes such as domestic use, irrigation, livestock watering, fishing, recreation and alike.
Field sampling was carried out during two consecutive years, both during the wet season (November, 2017 and 2018). Sampling sites were selected from a range of levels of accessibility to anthropogenic activities and existing variations in ecological characteristics of the habitats (Fig. 2).
Physico-chemical water quality measurements were performed both onsite during sampling and in the laboratory. Water temperature, pH, day-time dissolved oxygen and electrical conductivity were measured onsite using a multi-probe meter (HQ40d Single-Input Multi-parameter Digital Meter; Hach Company, Loveland, USA). Chlorophyll-a was measured onsite using a hand-held fluorometer (Aqua Fluor; Turner Designs, San Jose, USA). Turbidity was measured using a turbidity meter (Wag-WT3020; Halma PLC Company, Amersham, UK). In addition, water samples were collected from each sampling site in polyethylene bottles for analysis of other physico-chemical parameters. Afterwards, the samples were transported (in the dark) to the laboratory in an ice cooled box. Concentrations of total suspended solids, chemical oxygen demand, nitrate, ammonia, total phosphorus, orthophosphate, chloride, alkalinity, total hardness, calcium and magnesium were determined according to the standard method  at the laboratory of the Environmental Health Science and Technology of Jimma University, Ethiopia.
Physical characteristics of snail habitats, such as water depth, transparency and ambient temperature were measured at each sampling site using the United States Environmental Protection Agency habitat assessment protocol . Water depth was measured using a graduated stick calibrated in centimeters. Transparency of water was determined with a Secchi disk 30 cm in diameter attached to a calibrated cord. Ambient temperature was measured using a mercury-in-glass thermometer. Percentage of aquatic macrophyte (emergent, submerged and floating) cover was visually estimated by a simple estimation of the proportion of the site covered by aquatic plants within a 500 m stretch considering the sampling site as the center. The percentage of the aquatic macrophyte cover was classified into five groups: very low (< 10%); low (10–35%); moderate (35–65%); high (65–90%); and very high (> 90%) . Canopy cover was estimated visually based on the percentage of shade . At each sampling site, substrate was carefully assessed by observation and classified into grass, silt, detritus, sand, gravel, boulders and bedrock . Presence or absence (1/0) of anthropogenic activities in and around each sampling site was assessed by observation following the method used by Mereta et al. .
A hand-held global positioning system (GPS) unit was used to record geographical coordinate readings at each sampling site. The map of the study area showing the locations of sampling sites was created using the Geographic Information System (GIS) software ArcGIS10.3.
Freshwater snails survey
Snails were collected using a standard scoop net with a mesh size of 300 µm supported by a metal frame . Water samples were collected for the determination of physico-chemical parameters at the sampling sites where snails were collected. At each sampling site, all substrates were thoroughly searched to collect snails. Sampling time was fixed at 30 min per site. Sampling area per site was approximately 5 m2, whereas lengths of 10 m along the Lake shoreline and wetland were sampled. Afterwards, snails collected from each sampling site were transferred to plastic vials containing 80% ethanol to transport to the laboratory of the Department of Environmental Health Science and Technology, Jimma University, for identification. In the laboratory, each snail (except for specimens of the genera Physa, Ferrisia and Burnupia) was identified to species level according to the shell morphology using the field guide to African freshwater snails  and individuals were counted to determine the number of each snail species collected per sampling site.
Predators and competitors survey
To collect data on invertebrate predators and competitors of the freshwater snail species, macroinvertebrates were sampled from the same sampling site where snails were collected. Macroinvertebrate sampling was carried out using a standard hand-held rectangular drag (20 × 30 cm) with a cutting metal frame covered with a 300 µm mesh net. Macroinvertebrate sampling was carried out for 10 min over a distance of 10 m by disturbing the bottom sediment by foot and by sweeping the kick net between macrophytes to dislodge the macroinvertebrates . Samples of macroinvertebrates collected from each sampling site were sorted in the field, placed in labeled plastic vials containing 80% ethanol and transported to the laboratory for identification. In the laboratory, macroinvertebrates were identified to the family level based on morphology using a stereomicroscope (10×) and identification key  and then categorized according to functional feeding groups: predators, scrapers, gatherer-collectors, filterer-collectors, and shredders . Afterwards, scrapers and macroinvertebrates belonging to the family Physidae  were considered as competitors of snails, whereas invertebrates such as Psychodidae , Dytiscidae , Hydrophilidae , Belostomatidae , Odonata  and Glossiphoniidae  were considered as predators of snails. Presence or absence (1/0) and abundance of predators and competitors were considered as biotic factors to assess the occurrence and abundance of freshwater snail species.
Data exploration and regression modeling were performed using R software (Version 3.5.1) . A Shapiro–Wilk normality test showed that data on abundance of each snail species did not show a normal distribution (P < 0.05). Therefore, a non-parametric Kruskal–Wallis ANOVA was used to test whether significant differences in the abundance of snail species existed between the agro-ecological zones, and in relation to presence/absence of predators and competitors. All statistical tests were performed at the significance level of 0.05.
Generalized linear models were used to identify the variables affecting the occurrence and abundance of the snail species. Logistic regression models were used to predict the occurrence of the snail species and zero-inflated Poisson regression models were used to predict the abundance. Because a large number of zero values reduce the power of the test statistic , only snail species with frequency of occurrence greater than 20% were used for regression analysis. Therefore, five snail species (i.e. B. pfeifferi, B. sudanica, L. natalensis, L. truncatula and Bu. globosus) were selected for regression analysis in this study.
During the preliminary analysis, the dataset was explored to detect outliers and collinearity among predictor variables to decrease uncertainty of the model . To find a set of predictor variables that does not contain collinearity, Spearman’s rank order correlation coefficient was determined to generate a matrix of pairwise correlations between all predictor variables. As correlation coefficients only show pairwise correlations, we calculated variance inflation factors (VIFs) to assess which predictor variables are collinear and should be dropped before starting the analysis. This procedure continued until no further collinearity existed .
A stepwise backward selection procedure was followed to build the model starting from the full model. The model with the lowest Akaike information criterion (AIC) value was retained as the optimal model . The goodness-of-fit of the models was assessed using the relation between the residuals and predictor variables and the normality of the residuals was tested using a QQ-plot (probability plot). Retained models were only considered reliable if no relations between the residuals and the predictor variables were visually observed and residuals were normally distributed .
Occurrence and abundance of freshwater snails
A total of 2700 individual snails were collected from the 174 sampling sites. All collected snails were of the subclass Pulmonata which were grouped into 4 families: Planorbidae, Lymnaeidae, Physidae and Ancylidae. The snails were encountered in 137 sampling sites (79%).
The most abundant family was Planorbidae (n = 1494) which accounted for 55.3% of the total number of snails followed by the family Lymnaeidae, Physidae and Ancyclidae accounting for 36.6%, 8% and 0.48% of the total number of snails, respectively. The spatial distribution of the snail species shows that the tropical rainy climate zone had a higher relative abundance (57%) compared to the warm temperate rainy climate zone (43%). The frequency of occurrence of snails was higher in the warm temperate rainy climate zone accounting for 80% of the total sampling sites compared to the tropical rainy climate zone (77.5%). In addition, the warm temperate rainy climate zone supported higher snail richness (7 species and 3 genera) compared to the tropical rainy climate zone (7 species). Planorbidae and Lymnaeidae families were represented in both agro-ecological zones, whereas Physidae and Ancyclidae were collected only from the warm temperate rainy climate zone. The frequency of occurrence and abundance of the freshwater snail species collected during the study period are presented in Table 1.
Among all genera of snail species collected, Bu. globosus was the most abundant snail species accounting for 32% of the total number of snails and it was encountered in 45% of the sampling sites. The species was mainly found in the tropical rainy climate zone. Lymnaea natalensis was the second most abundant snail species accounting for 22% of the total number of snails. It was recorded at 44% of the sampling sites and most of which were encountered in the warm temperate rainy climate zone. Biomphalaria pfeifferi and B. sudanica were collected from both the tropical rainy climate and warm temperate rainy climate zones and were recorded at 36% and 24% of the sampling sites, respectively. Burnupia was the least abundant snail genus collected from 1% of the sampling sites. A Kruskal–Wallis test indicated that the tropical rainy climate zone supported a significantly higher abundance of B. pfeifferi (χ2 = 9.71, df = 1, P = 0.0018), L. truncatula (χ2 = 36.47, df = 1, P < 0.0001), Bu. globosus (χ2 = 7.89, df = 1, P = 0.0049) and a higher abundance of the total number of snails (χ2 = 5.37, df = 1, P = 0.0204). A significantly higher abundance of predators (χ2 = 11.19, df = 1, P = 0008) was encountered in the tropical rainy climate zone, whereas a significantly higher abundance of competitors (χ2 = 13.30, df = 1, P = 0.0003) was observed in the warm temperate rainy climate zone.
Factors affecting the occurrence of freshwater snail species
In the present study, the logistic regression model included water depth, turbidity, emergent macrophyte cover and human settlement as important factors determining the occurrence of B. pfeifferi and the model explained 23.4% of the variation (Additional file 1: Table S1). The occurrence of B. pfeifferi increased with increasing water depth, but decreased with increasing turbidity. Emergent macrophyte cover showed a positive association with regard to the occurrence of B. pfeifferi in the presence of human settlements (Additional file 1: Table S1). The logistic regression model determining the occurrence of B. sudanica selected water temperature, concentration of dissolved oxygen, submerged macrophyte cover and human settlements as determining factors (13.7% of variance explained) (Additional file 1: Table S1). The occurrence of B. sudanica increased with increasing water temperature, but decreased with increasing concentration of dissolved oxygen. Based on model output, submerged macrophyte cover promotes the occurrence of B. sudanica at sites situated near human settlements.
Water temperature, calcium concentration, magnesium concentration and agro-ecological zone were selected by the logistic regression model as important factors determining the occurrence of L. truncatula explaining 82.8% of the variation (Additional file 1: Table S1). The model demonstrated that increasing calcium concentration promotes the occurrence of L. truncatula, while increasing magnesium concentration and water temperature restrict its occurrence. Lymnaea truncatula is more likely to occur in the warm temperate rainy climate zone (Additional file 1: Table S1). There was no variable selected by the logistic regression model to explain the occurrence of L. natalensis.
According to the logistic regression model, water depth, alkalinity, chloride concentration, substrate type, agro-ecological zone and fishing explained the occurrence of Bu. globosus (36.7% of variance explained) (Additional file 1: Table S1). Based on the output of the model we found that increasing water depth and alkalinity positively influences the occurrence of Bu. globosus, while increasing chloride concentration leads to a decrease in the occurrence of Bu. globosus. Bulinus globosus frequently occurred in the warm temperate rainy climate zone, but was less likely to occur at sites with silt and in the presence of fishing activity (Additional file 1: Table S1).
Factors affecting the abundance of freshwater snail species
According to the zero-inflated Poisson regression model, the abundance of B. pfeifferi was positively correlated with physico-chemical water quality parameters such as water temperature and chemical oxygen demand, but negatively correlated with turbidity and dissolved oxygen (Additional file 2: Table S2). A higher abundance of B. pfeifferi was found at sites with submerged macrophyte cover and emergent macrophyte cover and at sites characterized by the presence of anthropogenic activities (i.e. human settlements, farming and clothes washing). The species was found at lower abundances at sites where bathing and swimming, drainage of land and car washing took place.
The output of the model showed that the abundance of B. sudanica increased with increasing water temperature, alkalinity and canopy cover, but decreased with increasing dissolved oxygen concentration, chlorophyll-a and chloride concentration in the presence of human settlements. Biomphalaria sudanica was abundant at sites characterized by submerged macrophyte cover (Additional file 2: Table S2). Abundance of L. natalensis increased with increasing alkalinity and competitor abundance, but decreased with increasing ammonia concentration and dissolved oxygen (Additional file 2: Table S2). A zero-inflated Poisson regression model revealed that L. natalensis was found at higher abundances at sites where grazing by cattle was observed. The abundance of L. truncatula was positively related to nitrate concentration and magnesium concentration, but negatively related to water temperature, chlorophyll-a, ammonia concentration and calcium concentration (Additional file 2: Table S2). Anthropogenic activities (i.e. clothes washing and boating) reduced the abundance of L. truncatula, but car washing and the presence of a silviculture promoted the abundance of L. truncatula (Additional file 2: Table S2).
Water temperature, water depth, alkalinity, nitrate concentration, chloride concentration, ammonia concentration, predator abundance, canopy cover, floating macrophyte cover, emergent macrophyte cover, substrate type, agro-ecological zone and fishing were selected by the model as the main variables determining the abundance of Bu. globosus (Additional file 2: Table S2). Abundance of Bu. globosus increased with increasing water temperature, water depth, nitrate concentration, alkalinity concentration, predator abundance, but decreased with increasing chloride concentration, ammonia concentration and canopy cover (Additional file 2: Table S2). Bulinus globosus was found at higher abundances at sites with floating macrophyte cover in the presence of silt, but was less abundant at sites with emergent macrophyte cover and where detritus was present (Additional file 2: Table S2).
This ecological investigation of gastropods focused on identifying environmental and biotic factors significantly affecting the occurrence and abundance of freshwater snail intermediate hosts in the Ethiopian Rift Valley. This study revealed that the habitat preference of freshwater snail species depends on physico-chemical water quality parameters, physical characteristics of habitat, biological factors and anthropogenic activities.
Regression analysis between physico-chemical water quality parameters and the occurrence and abundance of snail species showed that B. pfeifferi, B. sudanica, L. truncatula and Bu. globosus were abundant at sites with a relatively higher water temperature. This might be due to the fact that a higher temperature increases the food availability . Higher temperature also crucially increases the snail metabolic rate and therefore increases the size of the snail population by reducing the duration of the development periods [35, 63]. On the other hand, L. truncatula showed an opposite trend and less likely occurred at higher temperatures in this study. This could be due to the species being more sensitive to thermal stress compared to the other species . Previous studies pointed out that water temperature is one of the important factors determining the occurrence and abundance of snail species [26, 62, 64, 65]. This study revealed that B. pfeifferi, B. sudanica and L. natalensis were more abundant at sites with a low dissolved oxygen. Biomphalaria pfeifferi was more abundant at sites with a high chemical oxygen demand, whereas it was less likely to occur and less abundant at sites with a high turbidity. A positive correlation of snail species to low values of dissolved oxygen and chemical oxygen demand could be explained by the ability of snail species to occupy water bodies with a high content in organic matter . Outputs of the models indicated that B. sudanica, L. natalensis and Bu. globosus were abundant at a high water alkalinity, whereas B. sudanica and Bu. globosus were less abundant at sites with a high chloride concentration. Our findings match with the results of the study by Oloyede et al. , who recorded high abundance of snails in alkaline water. Calcium and magnesium concentrations were also identified as important factors significantly affecting the occurrence and abundance of L. truncatula. The ratio of calcium to magnesium is important for calcium uptake for the development of snail shells [67, 68]. A higher nitrate concentration is an indication of eutrophication and thus favours the occurrence of the snails .
Previous studies reported that water depth was an important ecological factor determining the distribution of snail species [33, 35, 67]. Similarly, in our study, we found a positive relationship between water depth and the occurrence of both B. pfeifferi and Bu. globosus and a negative relationship between water depth and the occurrence of L. truncatula. On the other hand, the abundance of B. pfeifferi was positively associated to water depth. Moreover, we found that B. sudanica was abundant at sites characterized by a high canopy cover, which supplies the snails with shade and shelter . However, Bu. globosus was less abundant at sites characterized by a high canopy cover which is probably related to indirect effects of canopy cover on snails through a negative impact on the ability of sunlight to reach the bottom and decrease primary production, which represents the food supply for snails . A negative correlation between chlorophyll-a content and the abundance of both B. sudanica and L. truncatula observed in our study is in line with previous studies [69, 72, 73]. Chlorophyll-a is necessary for the productivity of plants and algae in a water body which is used as a food source for snails. However, the level of chlorophyll-a measured in this study, 14.03 μg/l, is above 10 μg/l that support high concentrations of harmful algae to aquatic organisms . The increase in presence of B. pfeifferi and the increase in abundance of B. sudanica and Bu. globosus at sites with a high macrophyte cover observed in our study, indicate that food availability and the accessibility of aquatic weeds as a suitable surface on which egg masses of snails are deposited play an important role [21, 62, 75, 76]. On the contrary, the decrease in abundance of Bu. globosus at sites with macrophyte cover could be due to the fact that a dense emergent macrophyte cover prevents oxygen which promotes growth of snails .
A positive relationship between Bu. globosus abundance and predator abundance and between L. natalensis abundance and competitor abundance observed in the present study might be due to the presence of macrophyte cover that may suppress the predation and competition activities . This finding is consistent with the earlier reports by Dejenie et al. , who reported that predator abundance was positively associated with the increasing number of aquatic invertebrates. This finding, however, contradicts that of Yigezu et al. , who found that the abundance of snails decreases as the abundance of predators increases (Tables 2, 3).
According to the survey, it was observed that the occurrence and abundance of snail species is associated with anthropogenic activities such as human settlement, farming, bathing and swimming, clothes washing, grazing, drainage of land, car washing, boating, fishing and silviculture in the study area. The increased presence and high abundance of snail species in habitats disturbed by anthropogenic activities could be due to the high concentration of organic matter and dissolved ions which are conducive for snail species [30, 34, 35, 78,79,80]. The results of the present study show that Bu. globosus was likely to occur, and more abundant, in the warm temperate rainy climate zone. However, L. truncatula frequently occurred, and was less abundant, in the warm temperate rainy climate zone compared to the tropical rainy climate zone. Differences in ecological and climatic conditions between the agro-ecological zones may be responsible for this result.
Implication for snail control
Results of this study revealed that snail species occurred frequently and were abundant in habitats disturbed by anthropogenic activities such as settlement, farming, clothes washing, car washing, grazing, silviculture, and bathing and swimming. Hence, integrated snail control strategies should be considered as a priority to reduce the presence and density of freshwater snail intermediate hosts and thus to control the spread of snail-borne diseases at a local scale. Protection of water bodies from disturbance by anthropogenic activities (i.e. through regulation of human settlement in the areas) may be a good strategy to reduce the concentration of organic matter and dissolved ions in aquatic ecosystem, which are conducive for presence of snails. However, contact with cercariae-infested water is a decisive factor related to the risk of acquiring snail-borne diseases. Hence, a safe water supply may be an alternative approach to reduce the likelihood of cercariae infecting a final host. On the other hand, even if safe water supplies reduce such water contact, they may not completely prevent it. Therefore, a reduction in infection could be attempted through behavioral changes, such as boiling water before drinking it, preventing animals from getting to the water bodies, and the use of commercial soap and a locally available soapberry endod (Phytolacca dodecandra) as a detergent during water contact, which appears to provide some protection from infection.
Generalized linear models indicated that physico-chemical water quality parameters, physical habitat characteristics and biotic factors were found to be the main variables determining the occurrence and abundance of snail species in the Ethiopian Rift Valley region. This study also revealed that snail species frequently occurred and were abundant in habitats disturbed by anthropogenic activities. Therefore, the findings reported herein suggest that integrated snail control strategies should be considered to control snails via protection of water bodies from disturbance by anthropogenic activities. In this way, it is possible to reduce the concentration of organic matter and dissolved ions in aquatic ecosystem which are conducive to the presence of snails. This study may be used as baseline for further studies on trematode infections in snail intermediate hosts from these water bodies and epidemiology of schistosomiasis among lakeshore communities to ascertain whether active transmission is occurring within these areas.
Availability of data and materials
Data supporting the conclusions of this article are included within the article and its additional files. The dataset generated and/or analyzed during the present study is available from the corresponding author.
nephelometric turbidity unit
Pyron M, Brown KM. Introduction to mollusca and the class Gastropoda. In: Thorph JH, Rogers DC, editors. Ecology and general biology: Thorp and Covich’s freshwater invertebrates. Cambridge: Academic Press; 2015. p. 383–421.
Soldánová M, Selbach C, Kalbe M, Kostadinova A, Sures B. Swimmer’s itch: etiology, impact, and risk factors in Europe. Trends Parasitol. 2013;29:65–74.
WHO. Working to overcome the global impact of neglected tropical diseases. The first WHO report on neglected tropical diseases. Geneva: World Health Organization; 2010.
Steinmann P, Keiser J, Bos R, Tanner M, Utzinger J. Schistosomiasis and water resources development: systematic review, meta-analysis, and estimates of people at risk. Lancet Infect Dis. 2006;6:411–25.
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:1665–75.
Kloos H, Lo CT, Birrie H, Ayele T, Tedla S, Tsegay F. Schistosomiasis in Ethiopia. Soc Sci Med. 1988;26:803–27.
Arancha A, Anegagrie M, Zewdie D, Benito A. Schistosoma mansoni in a rural community of Ethiopia. The 27th European congress of clinical microbiology and infectious diseases, 22–25 April 2017, Veinna: Austria; 2017.
Alebie G, Erko B, Aemero M, Petros B. Epidemiological study on Schistosoma mansoni infection in Sanja area, Amhara region, Ethiopia. Parasit Vectors. 2014;7:15.
Amsalu G, Mekonnen Z, Erko E. A new focus of schistosomiasis mansoni in Hayk town, northeastern Ethiopia. BMC Res Notes. 2015;8:22.
Bajiro M, Dana D, Levecke B. Prevalence and intensity of Schistosoma mansoni infections among schoolchildren attending primary schools in an urban setting in southwest, Ethiopia. BMC Res Notes. 2017;10:677.
Singh KL, Singh D, Singh VK. Toxicity of Bauhinia variegata and Mimusops elengi with plant mollucicides agaainst Lymnaea acuminata. J Biol Earth Sci. 2012;2:B76–82.
Degheidy NS, Al-Malki JS. Epidemiological studies of fasciolosis in human and animals at Taif, Saudi Arabia. World Appl Sci J. 2012;19:1099–104.
Akinwale O, Kane R, Rollinson D, Stothard J, Ajayi M, Akande D, et al. Molecular approaches to the identification of Bulinus species in southwest Nigeria and observations on natural snail infections with schistosomes. J Helminthol. 2011;85:283–93.
Hanington PC, Forys MA, Loker ES. A somatically diversified defense factor, FREP3, is a determinant of snail resistance to schistosome infection. PLoS Neglect Trop Dis. 2012;6:e1591.
Paniker CJ, Ghost S. Paniker’s textbook of medical parasitology, 7th ed. In: Arora S, editor. New Delhi: JP Medical Ltd; 2017.
Brown SP, Renaud F, Guegan JF, Thomas F. Evolution of trophic transmission in parasites: the need to reach a mating place? J Evol Biol. 2001;14:815–20.
Doughty BL. Schistosomes and other trematodes. Medical microbiology. 4th ed. Galveston: University of Texas Medical Branch at Galveston; 1996.
Abhay SR, Simon LG, Hotez JP, Tsuji M. Medical parasitology. Vademecum Parasitology Landes Bioscience, Austin: Texas USA; 2009:195–221.
Tolley-Jordan L, Huryn AD, Bogan AE. Effects of land-use change on a diverse pleurocerid snail assemblage. Aquat Conserv. 2015;25:235–49.
Bogan AE, Dillo RT. The ecology of freshwater mollusks. J N Am Benthol Soc. 2001;20:325–8.
Kloos H, Souza CD, Gazzinelli A, Soares Filho BS, Temba PD, 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.
Yusuf Z, Dagne K, Erko B, Siemuri O. Polyploidy in Bulinid Snails, with emphasis on Bulinus truncatus/tropicus complex (Planorbidae: pulmonate mollusks) from various localities in Ethiopia. World J Cell Biol Gen. 2017;3:11–20.
Erko B, Balcha F, Kifle D. The ecology of Biomphalaria sudanica in Lake Ziway, Ethiopia. Afr J Ecol. 2006;44:347–52.
Erko B, Medhin G, Berhe N, Abebe F, Gebre-Michael T, Gundersen S. Epidemiological studies on intestinal schistosomiasis in Wondo Genet, southern Ethiopia. Ethiop Med J. 2002;40:29–39.
Dida GO, Gelder FB, Anyona DN, Matano AS, Abuom PO, Adoka SO, et al. Distribution and abundance of schistosomiasis and fascioliasis host snails along the Mara river in Kenya and Tanzania. Infect Ecol Epidemiol. 2014;4:24281.
Utzinger J, Tanner M. Microhabitat preferences of Biomphalaria pfeifferi and Lymnaea natalensis in a natural and a man-made habitat in southeastern Tanzania. Mem Inst Oswaldo Cruz. 2000;95:287–94.
Abdisa T, Jilo K. Review on ovine fasciolosis in Ethiopian. J Vet Sci Res. 2017;2:1–10.
Ayalneh B, Bogale B, Dagnachew S. Review on Ovine fasciolosis in Ethiopia. Acta Parasitol Glob. 2018;9:7–14.
Aliyu AA, Ajogi I, Ajanusi O, Reuben R. Epidemiological studies of Fasciola gigantica in cattle in Zaria, Nigeria using coprology and serology. J Public Health Epidemiol. 2014;6:85–91.
Usman A, Babeker E, Malann Y. Effects of some physico-chemical parameters on prevalence of intermediate host of animal trematodes in Bauchi State, Nigeria. Sci World J. 2017;12:94–7.
Barkia H, Barkia A, Yacoubi R, Guemri Y, El Madhi Y, Belghyti D. Trematode infection among freshwater gastropods in the Gharb area, Morocco. Environ Earth Sci. 2015;5:174–81.
Ofulla AV, Adoka SO, Anyona DN, Abuom PO, Karanja D, Vulule JM, et al. Spatial distribution and habitat characterization of schistosomiasis host snails in lake and land habitats of western Kenya. Lakes Reserv Res Manag. 2013;18:197–215.
Oloyede OO, Otarigho B, Morenikeji O. Diversity, distribution and abundance of freshwater snails in Eleyele dam, Ibadan, southwest Nigeria. Zool Ecol. 2017;27:35–43.
Omonijo AO, Asaolu SO, Ofoezie IE. Ecology of Schistosoma snail vectors in Ado-Ekiti local government area, Ekiti State, Nigeria. Int J Pure Appl Zool. 2016;4:77–84.
Yigezu G, Mandefro B, Mengesha Y, Yewhalaw D, Beyene A, Ahmednur M, et al. Habitat suitability modelling for predicting potential habitats of freshwater snail intermediate hosts in Omo-Gibe river basin, southwest Ethiopia. Ecol Inform. 2018;45:70–80.
Gundersen SG, Birrie H, Torvik HP, Scherbaum H. Control of Schistosoma mansoni in the Blue Nile Valley of western Ethiopia by mass chemotherapy and focal snail control: a primary health care experience. Trans R Soc Trop Med Hyg. 1990;84:819–25.
Mitiku H, Legesse M, Teklemariam Z, Erko B. Transmission of Schistosoma mansoni in Tikur Wuha area, southern Ethiopia. Ethiop J Health Dev. 2010;24:180–4.
Mengesha G, Mamo Y, Sahle K, Elphick C, Bekele A. Effects of land-use on birds diversity in and around Lake Zeway, Ethiopia. J Sci Dev. 2014;2:5–22.
Deganovsky AM, Getahun BA. Water balance and level regime of Ethiopian lakes as integral indicators of climate change. The 12th World Lake Conference, Jaipir, Rajasthan, India; 2008.
Worako AW. Evaluation of the water quality status of lake Hawassa by using water quality index, southern Ethiopia. J Water Res Environ Eng. 2015;7:58–65.
APHA, AWWA, WPCF. Standard methods for the examination of water and wastewater. 19th ed. Washington, DC: American Public Health Association; 1995.
Barbour MT, Gerritsen J, Snyder BD, Stribling JB. Rapid bioassessment protocols for use in streams and wadeable rivers: periphyton, benthic macroinvertebrates and fish. 2nd ed. Washington, DC: US Environmental Protection Agency, Office of Water; 1999.
Parsons M, Ransom G, Thoms M, Norris R. Australian river assessment system: AusRivAS physical and chemica assessment module. Canberra: Environmental Australia; 2002. p. 47.
Posa MRC, Sodhi NS. Effects of anthropogenic land use on forest birds and butterflies in Subic Bay, Philippines. Biol Conserv. 2006;129:256–70.
McCreadie JW, Adler PH. The roles of abiotic factors, dispersal, and species interactions in structuring stream assemblages of black flies (Diptera: Simuliidae). Aquat Biosyst. 2012;8:14.
Mereta ST, Boets P, Bayih AA, Malu A, Ephrem Z, Sisay A, et al. Analysis of environmental factors determining the abundance and diversity of macroinvertebrate taxa in natural wetlands of southwest Ethiopia. Ecol Inform. 2012;7:52–61.
Opisa S, Odiere MR, Jura WG, Karanja DM, Mwinzi PN. Malacological survey and geographical distribution of vector snails for schistosomiasis within informal settlements of Kisumu city, western Kenya. Parasit Vectors. 2011;4:226.
Mandahl-Barth G. Key to the identification of east and central African freshwater snails of medical and veterinary importance. Bull World Health Organ. 1962;27:135–50.
Gabriels W, Lock K, De Pauw N, Goethals PL. Multimetric macroinvertebrate index Flanders (MMIF) for biological assessment of rivers and lakes in Flanders (Belgium). Limnologica. 2010;40:199–207.
Bouchard RW Jr. Guide to aquatic invertebrates of the Upper Midwest. In: Ferrington LC, Karius ML, editors. Identification manual for students, citizen monitors, and aquatic resource professionals. St. Paul: Water Resources Center, University of Minnesota; 2004. p. 208.
Tomanova S, Goitia E, Helešic J. Trophic levels and functional feeding groups of macroinvertebrates in Neotropical streams. Hydrobiologia. 2006;556:251–64.
Zukowski S, Walker KF. Freshwater snails in competition: alien Physa acuta (Physidae) and native Glyptophysa gibbosa (Planorbidae) in the river Murray, south Australia. Mar Freshw Res. 2009;60:999–1005.
El Bardicy S, Tadros M, Yousif F, Hafez S. Predatory activity of Psychoda alternata say (Diptera: Psychodidae) larvae on Biomphalaria glabrata and Lymnaea natalensis snails and the free-living larval stages of Schistosoma mansoni. Aust Bas Appl Sci J. 2009;3:4503–9.
Dillon RT. The ecology of freshwater molluscs. 1st ed. Cambridge: Cambridge University Press; 2000.
Inoda T, Inoda Y, Rullan JK. Larvae of the water scavenger beetle, Hydrophilus acuminatus (Coleoptera: Hydrophilidae) are specialist predators of snails. Eur J Entomol. 2015;112:145–50.
Younes A, El-Sherief H, Gawish F, Mahmoud M. Sphaerodema urinator Duforas (Hemiptera: Belostomatidae) as a predator of Fasciola intermediate host, Lymnaea natalensis Krauss. Egypt J Biol Pest Control. 2016;26:191–6.
Younes A, El-Sherif H, Gawish F, Mahmoud M. Potential of Hemianax ephippiger (Odonata-Aeshnidae) nymph as predator of Fasciola intermediate host, Lymnaea natalensis. Asian Pac J Trop Biomed. 2015;5:671–5.
Brönmark C, Malmqvist B. Interactions between the leech Glossiphonia complanata and its gastropod prey. Oecologia. 1986;69:268–76.
Venables WN, Smith DM. An introduction to R: notes on R, a programming environment for data analysis and graphics. Vienna: R Development Core Team; 2008.
Pramual P, Wongpakam K. Seasonal variation of black fly (Diptera: Simuliidae) species diversity and community structure in tropical streams of Thailand. Entomol Sci. 2010;13:17–28.
Zuur AF, Leona EN, Walker NJ, Saveliev AA, Smith GM. Mixed effects models and extensions in ecology with R. New York: Springer Science and Business Media; 2009.
Malek EA. Factors conditioning the habitat of bilharziasis intermediate hosts of the family Planorbidae. Bull World Health Organ. 1958;18:785.
De La Rocque S, Rioux JA, Slingenbergh J. Climate change: effects on animal disease systems and implications for surveillance and control. Rev Sci Tech. 2008;27:339–54.
Hofkin BV, Mkoji GM, Koech DK, Loker ES. Control of schistosome-transmitting snails in Kenya by the north American crayfish Procambarus clarkii. J Trop Med Hyg. 1991;45:339–44.
Kazibwe F, Makanga B, Rubaire-Akiiki C, Ouma J, Kariuki C, Kabatereine N, 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:433–44.
Gallardo-Mayenco A, Toja J. Spatio-temporal distribution of Simuliids (Diptera) and associated environmental factors in two mediterranean basins of southern Spain. Limnetica. 2002;21:47–57.
Abdulkadir F, Maikaje D, Umar Y. The influence of physico-chemical and ecological factors on the distribution of freshwater snails in Manchok water intake in Kaduna state, Nigeria. Niger J Chem Res. 2013;18:58–66.
Ofoezie IE. Distribution of freshwater snails in the man-made Oyan reservoir, Ogun state, Nigeria. Hydrobiologia. 1999;416:181–91.
Daldorph P, Thomas J. The effect of nutrient enrichment on a freshwater community dominated by macrophytes and molluscs and its relevance to snail control. J Appl Biol. 1991;28:685–702.
Walker PD, Wijnhoven S, Van der Velde G. Macrophyte presence and growth form influence macroinvertebrate community structure. Aquat Bot. 2013;104:80–7.
Aguiar F, Ferreira M, Pinto P. Relative influence of environmental variables on macroinvertebrate assemblages from an Iberian basin. J N Am Benthol Soc. 2002;21:43–53.
Hussein MA, Obuid-Allah AH, Mahmoud AA, Fangary HM. Population dynamics of freshwater snails (Mollusca: Gastropoda) at Qena governorate, upper Egypt. Egypt Acad J Biol Sci. 2011;3:11–22.
Tchakonté S, Ajeagah GA, Diomandé D, Camara AI, Ngassam P. Diversity, dynamic and ecology of freshwater snails related to environmental factors in urban and suburban streams in Douala-Cameroon (central Africa). Aquat Ecol. 2014;48:379–95.
Kansas Departement of Health and Environment (KDHE). Water quality standards white paper: chlorophyll-a criteria for public water supply lakes or reservoirs. KDHE bureau of water. Topeka: KDHE; 2011.
El Deeb FA, El-Shenawy NS, Soliman MF, Mansour SA. Freshwater snail distribution related to physico-chemical parameters and aquatic macrophytes in Giza and Kafr El-Shiekh governorates, Egypt. J Vet Sci Res. 2017;3:8–13.
Sturrock RF. The schistosomes and their intermediate hosts. In: Mahmoud AAF, editor. Schistosomiasis. London: Imperial College Press; 2001. p. 7–83.
Dejenie T, Yohannes M, Assmelash T. Characterization of mosquito breeding sites in and in the vicinity of Tigray microdams. Ethiop J Health Sci. 2011;21:57–66.
Camara I, Bony Y, Diomandé D, Edia O, Konan F, Kouassi C, et al. Freshwater snail distribution related to environmental factors in Banco national park, an urban reserve in the Ivory Coast (west Africa). Afr Zool. 2012;47:160–8.
Ernould JC, Kaman A, Labbo R, Couret D, Chippaux JP. Recent urban growth and urinary schistosomiasis in Niamey, Niger. Trop Med Int Health. 2000;5:431–7.
Madsen H, Coulibaly G, Furu P. Distribution of freshwater snails in the river Niger basin in Mali with special reference to the intermediate hosts of schistosomes. Hydrobiologia. 1987;146:77–88.
We would like to acknowledge the VLIR-UOS Network Ethiopia Programme for financial and logistic support. Our gratitude goes to laboratory assistants at the Department of Environmental Health Science and Technology, Jimma University, for support during field and laboratory activities. Finally, the surveys would not have been completed without the support of the local community who allowed us access to sampling sites.
This work was financially supported by Flemish Interuniversity Council (VLIR-UOS Network Ethiopia Programme).
Ethics approval and consent to participate
Consent for publication
The authors declare that they have no competing interests.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
About this article
Cite this article
Olkeba, B.K., Boets, P., Mereta, S.T. et al. Environmental and biotic factors affecting freshwater snail intermediate hosts in the Ethiopian Rift Valley region. Parasites Vectors 13, 292 (2020). https://doi.org/10.1186/s13071-020-04163-6