Skip to main content

Molecular diversity of Bulinus species in Madziwa area, Shamva district in Zimbabwe: implications for urogenital schistosomiasis transmission



Bulinus species are freshwater snails that transmit the parasitic trematode Schistosoma haematobium. Despite their importance, the diversity of these intermediate host snails and their evolutionary history is still unclear in Zimbabwe. Bulinus globosus and B. truncatus collected from a urogenital schistosomiasis endemic region in the Madziwa area of Zimbabwe were characterized using molecular methods.


Malacological survey sites were mapped and snails were collected from water contact sites in four communities in the Madziwa area, Shamva district for a period of one year, at three-month intervals. Schistosoma haematobium infections in snails were determined by cercarial shedding and the partial mitochondrial cytochrome c oxidase subunit 1 gene (cox1) was used to investigate the phylogeny and genetic variability of the Bulinus spp. collected.


Among the 1570 Bulinus spp. snails collected, 30 (1.9%) B. globosus were shedding morphologically identified schistosomes. None of the B. truncatus snails were shedding. The mitochondrial cox1 data from 166 and 16 samples for B. globosus and B. truncatus, respectively, showed genetically diverse populations within the two species. Twelve cox1 haplotypes were found from the 166 B. globosus samples and three from the 16 B. truncatus samples with phylogenetic analysis showing that the haplotypes fall into well-supported clusters within their species groups. Both B. truncatus and B. globosus clustered into two distinct lineages. Overall, significant negative values for both Tajima’s D statistic and the Fu’s Fs statistic were observed for B. globosus and B. truncatus.


The study provided new insights into the levels of genetic diversity within B. globosus and additional information on B. truncatus collected from a small geographical area in Zimbabwe. Low prevalence levels of infection observed in the snails may reflect the low transmission level of urogenital schistosomiasis in the area. Our results contribute towards the understanding of the distribution and population genetic structure of Bulinus spp. supporting the mapping of the transmission or risk of transmission of urogenital schistosomiasis, particularly in Zimbabwe.


Freshwater snails of the genus Bulinus act as intermediate hosts for Schistosoma haematobium, the human blood fluke that causes the chronic and debilitating disease, urogenital schistosomiasis [1]. Bulinus species are extensively distributed throughout Africa, countries bordering the Mediterranean and some areas of the Middle East, but there is a considerable variation in compatibility between the different Bulinus species and schistosome parasites in different geographical areas [1, 2]. The host-parasite compatibility is influenced by both biotic and abiotic factors such as temperature, rainfall, water velocity, desiccation, salinity and the genetics of the snail hosts and the parasites [2,3,4,5,6,7].

Among the Bulinus species, B. globosus (Morelet, 1866) and B. truncatus (Audouin, 1827) are recognized as the most important intermediate hosts for S. haematobium and the species are distributed widely in Africa [8]. In Zimbabwe, B. globosus is the only intermediate host for S. haematobium [6]. The intermediate host snail prefers diverse habitat types and is more abundant in the northeast and southeast parts of the country with a patchy distribution in the southwest, correlating with the distribution of urogenital schistosomiasis in Zimbabwe. Bulinus truncatus is more common to the southwest preferring dams and marshy habitats [9].

While both B. truncatus and B. globosus are hermaphroditic, B. globosus preferentially cross-fertilizes if possible whereas B. truncatus preferentially self-fertilizes [10, 11]. Both species inhabit diverse habitats such as rivers, irrigation systems, ponds and lakes and are tolerant to differing water levels and seasonal changes, and have the ability to aestivate during dry seasons [8]. Geographical and temporal variability in ecological factors can cause large fluctuations in snail abundance within and between sites, possibly resulting in population extinction and/or recolonization. Selfing and population bottlenecks reduce genetic diversity within a population while increasing genetic differences between isolated populations [12]. Investigating snail population structure will help understand snail-parasite relationships and parasite transmission dynamics [13,14,15,16,17,18,19].

Due to the high levels of variation in shell morphology within and among Bulinus species, morphological identification can be subjective [8, 20]. Nevertheless, it has been the main method used to characterize Bulinus species in Zimbabwe [6,7,8], except for a few studies that utilized iso and alloenzyme analysis [21, 22] and random amplification of polymorphic nuclear DNA (RAPD) analysis [23]. However, several studies conducted in the region have characterized Bulinus populations using the DNA sequence analysis of the internal transcribed spacer (ITS) rDNA regions [18, 19, 24] and the mitochondrial cytochrome c oxidase subunit 1 (cox1) [13, 14, 17,18,19], with cox1 data being more informative for resolving species level affinities within the genus Bulinus [13, 14, 19]. Abe et al. [19] have reported six unique mitochondrial cox1 haplotypes for B. truncatus in Zimbabwe. However, no information exists regarding cox1 diversity for Zimbabwean B. globosus in the sequence databases. Here, we use mitochondrial cox1 data to identify and characterize Bulinus snails sampled from different water contact sites in a urogenital schistosomiasis endemic area, Madziwa in Zimbabwe and discuss the findings in relation to schistosomiasis transmission.


Study area

Natural populations of Bulinus species were collected from four rural communities in Mashonaland province, in a district highly endemic for urogenital schistosomiasis, as described in a previous study [25]. The communities are in semiarid areas and they rely on local rivers and streams for most of the household activities including bathing, fishing, swimming, washing, gardening and subsistence farming. There are only a few boreholes that are used by the residents for provision of drinking water but these are located ~ 5 km away. Within the four communities, there are distinct periods of high and low rainfall patterns, with the main rainfall period being from late October to April.

Sample collection, morphological analysis and patent testing

This study was part of a larger project investigating the burden of urogenital schistosomiasis in a highly endemic area of Zimbabwe, during which snail survey sites were selected by asking the residents to identify the sites which they frequently used for human water contact activity. The water contact sites are described in Fig. 1. Spatial variation in snail populations was investigated by collecting snails at each water contact point in the communities and all the surveyed sites were mapped using a global positioning system (GPS) (Trimble Navigation Ltd, California, USA). At each water contact site, sampling was performed at four different time points representing the rainy season, post-rainy season, winter and hot-dry season to investigate the temporal effects of season on the snail populations. Snail sampling at all sites was performed using a metal scoop or by handpicking the snails for 30 min at both the main parts of the water body and at the water edges. The snails were transported to the field laboratory where they were identified using shell morphology as described by Brown [8]. Snails with globose, ovate shells of small to medium size, and sinistral with a pseudobranch were grouped as belonging to the genus Bulinus. From the snails identified as Bulinus spp., B. globosus was identified by their truncate columella and microsculpture of nodules with short ridges or corrugations. Bulinus snails with a straight or evenly concave or twisted columella with no truncation were identified as B. truncatus. Both species were tested for patent parasitic infections by cercarial shedding. The snails were individually placed in flat-bottomed glass vials containing dechlorinated water and exposed to artificial light for a maximum of 4 h [26]. During the 4-hour period, the emergence of cercariae was checked at regular short intervals. Cercarial identification was based on morphology using a binocular microscope as described by Frandsen & Christensen [27]. Bifurcate cercariae were considered to be of mammalian origin. After shedding, the snails were counted and preserved in absolute ethanol for molecular analysis.

Fig. 1

Location of the water contact sites where Bulinus snails were collected in Madziwa area, Shamva district, Zimbabwe

DNA extraction and amplification

Among the collected and morphologically identified snails, a subset of 166 B. globosus and all the collected 16 B. truncatus were molecularly analysed. Up to eight adult snails were randomly selected per site for molecular analysis [23]. The preserved snails were transported to the Key Schistosomiasis Laboratory, National Institute of Parasitic Diseases, Shanghai, China, for molecular analysis. Snails were removed from ethanol and their soft parts were removed from the snail shell using forceps. For each snail, the soft part was then transferred to a clean Petri dish and immersed in TE buffer (10 mM Tris, 0.1 mM EDTA) pH 7.4 for 1 h. DNA extraction was performed from the head and foot region of the snail tissue using the DNeasy Blood and Tissue Kit (Qiagen, Crawley, UK) following the manufacturer’s instructions. DNA quantification and purity were measured using a Nanodrop ND-1000 spectrophotometer (Nanodrop Technologies Inc., Willington, USA).

Amplification and sequencing of mitochondrial DNA

The cox1 fragment was amplified by PCR, using the BulCox 5 forward (5′-CCT TTA AGA GGN CCT ATT GC-3′) and BulCox 14 reverse (5′-GGA AAT CAG TAM AYA AAA CCA GC-3′) primers previously described by Kane et al. [18] in a C1000TM thermal cycler (Bio-Rad, California, USA). Amplification was carried out in 25 μl PCR reactions containing a premixed master mix (2.5 μl of 20 mM MgCl2, 2.5 μl, 5× buffer, 2.5 μl of 20 mM dNTPs, one unit of Taq DNA polymerase), 1 μl (10 pmol) of each forward and reverse primer and 1 μl (10–100 ng) of the DNA template. The PCR conditions were as follows: denaturation at 95 °C for 5 min followed by 45 cycles of 95 °C for 5 min, 45 °C for 30 s and 72 °C for 1 min and a final extension step at 72 °C for 10 min. Negative controls (no template DNA) were included with each set of reactions. PCR products were viewed by 1% gel electrophoresis (Electrophoresis power supplies E455, CONSORT, Turnhout, Belgium) and read with a Molecular Imager, Gel DocTM XRt Imaging system (Bio-Rad). Amplicons were purified and Sanger sequenced in both forward and reverse directions using the BigDye v3.1 Terminator Cycle Sequencing Kit and read with the automated DNA fragment analyzer ABI-377 (Applied Biosystems, Carlsbad, USA) at Sangon Biotech, Shanghai, China.

Data analysis

All cox1 sequences were manually checked and edited using Sequencher v5.1 (http://www.Genecodes) to remove any ambiguities between forward and reverse strands. The consensus sequences of each sample were aligned in Sequencher v5.1 and any ambiguities between sequences were checked by visualisation of the original sequence chromatograms. Sequences were identified by BLAST searching [28] against the GenBank database.

Haplotype analysis

The consensus sequences from all samples were grouped and aligned in MacClade 4.05 and then collapsed together using Collapse v1.2 ( to identify samples with identical sequences/haplotypes. Within each site, and also overall, unique sequences and any group of identical sequences represented unique individual haplotypes. A group was defined as a collection of sequences regardless of the site where they were collected. Haplotypes were given a site-haplotype identifier code, consisting of a site number, letter H (haplotype for B. globosus) or T (haplotype for B. truncatus) and a number representing the different haplotypes in the area. Overall and in each site, the number of individual Bulinus snails presenting the same haplotype was recorded. Haplotype data were submitted to the EMBL/GenBank database under the accession numbers MN397785–MN397822.

Phylogenetic analysis

The haplotype data were exported to MEGA 6.0 [29] and aligned using the Clustal W algorithm [30]. Phylogenetic relationships between the haplotypes were inferred using the Neighbor-Joining (NJ), Minimum Evolution (ME) and Maximum Likelihood (ML) methods and the tree topologies were tested using 1000 bootstrap replicates in MEGA 6.0 [29]. The analysis was run using the Tamura-Nei with Gamma distribution nucleotide substitution model, which was the best-fit model for the data, inferred using the model test function in MEGA 6.0 [29, 31, 32]. A sequence for Biomphalaria glabrata (GenBank: AY380531) was used as the outgroup. Additional sequences for Bulinus spp. (B. globosus, B. africanus and B. truncatus) from Uganda, South Africa, Egypt, Tanzania, Kenya, Senegal, Portugal and Zimbabwe available on GenBank were also included in the analysis (Additional file 1: Table S1). The net nucleotide divergence between the main haplotype groups found within the study was calculated in MEGA 6.0 [29] with the Juke-Cantor model [33].

Population genetic analysis

The data were exported into DnaSP v6 [34] to determine the interpopulation and intrapopulation diversity of the Bulinus populations analysed. A population was defined as a collection of sequences from the same site or community or overall. Only 16 B. truncatus were collected, with not enough representatives per site and so were analysed as a single population. For B. globosus the overall diversity, diversity within the sites and also within the communities was measured. Haplotype diversity (h) and nucleotide diversity (π), was calculated using the Juke-Cantor model [33]. The diversification index (Fst) between each pair of sites and between B. globosus and B. truncatus was estimated in DnaSP v6 [34].

Test for selection

To investigate if there was significant selection occurring, Tajima’s D [35] and Fu’s Fs statistic neutrality tests for selection were conducted in DnaSP v6. Nucleotide divergence within and between the communities was calculated in DnaSP v6 [34].


Morphological identification and infection status of Bulinus species

From the collected snails, 1558 were putatively identified as Bulinus spp. based on their sinistral, ovate, globose shells. Among these, 1542 were morphologically identified as B. globosus and 16 as B. truncatus, based on the shape of the columella. Of these, a subset of 166 B. globosus and all 16 B. truncatus were molecularly characterized with a 100% match between the morphological and molecular species identification.

Among the B. globosus collected and tested for patent infections, 30/1542 (1.9%) were shedding schistosome cercariae (Table 1). The breakdown of infected snails by time point was as follows: 18 snails at baseline (end of February 2018; rainy season) and 12 snails at 6 months follow-up (September 2018; summer). No snails were shedding at 3 (winter season), 9 (early rainy season), or 12 (rainy season) months follow-up. None of the B. truncatus were found to be infected.

Table 1 Demographic information of the Bulinus samples analysed in Madziwa, Zimbabwe

Analysis of sequence data

After alignment and editing, the query length of the cox1 sequences was 607 bp. The query sequences for B. globosus were highly similar (> 98.0%) to B. globosus from Pietermaritzburg, South Africa (GenBank: AM286290 and AM286289). The B. truncatus query sequences had a high similarity (99.34%) to B. truncatus (GenBank: MG759464) from Zimbabwe.

Haplotype and population genetic diversity, and genetic divergence of B. globosus

Twelve unique B. globosus cox1 haplotypes were found from the 166 samples sequenced. The overall haplotype diversity exceeded 0.5 with a nucleotide diversity estimate of 0.002. The haplotypes were not highly divergent (pairwise distance of 0.002–0.007). Within the sites and within communities as expected, haplotype diversity was generally low but ranging between 0–0.786 (Table 2). The percentage occurrence of each haplotype differed among sites and communities. One common haplotype (H1) was found in all but one locality representing 114 (68.7%) of the 166 samples analysed. Some haplotypes were locally restricted being observed in one site only but for those represented by a considerable number of snails, there was a clear pattern of haplotype sharing within the river systems with high haplotype sharing being observed in sites that are close to each other or within the same river (Table 1). In Mupfure River, which had the highest number of haplotypes among the rivers, sites that were close together shared the same haplotypes (Tables 1, 2). Some sites had complete B. globosus homogeneity and this was observed mainly in Kamoyo stream where all sites had only one haplotype (H1). The homogeneity in Nduna/Kamoyo stream was also supported by lack of haplotype diversity. However, among the sites that had more than one haplotype, overall intra-site haplotype diversity estimates ranged between 0.250–0.786 and a considerable number of sites recorded a haplotype diversity of more than 0.5. The maximum number of haplotypes that could be found per site did not exceed four. The mean overall diversity for B. globosus was 0.002.

Table 2 Mitochondrial cytochrome c oxidase 1 diversity of B. globosus and B. truncatus collected from the Madziwa area, Zimbabwe

Nucleotide divergence between and within communities was low ranging between 0–0.003 (Table 3). Expectedly, the genetic divergence between B. globosus and B. truncatus was high (Fst = 0.920) supporting the species delimitation. When the pairwise fixation index (Fst) among the sites was estimated between the populations, a significant number of the populations showed significant levels of differentiation from each other whilst others showed no genetic differentiation. The pairwise Fst values ranged between 0–0.429 for B. globosus. A high genetic differentiation of > 30% and above was mainly found among different river systems demonstrating geographically defined gene population structuring (see pairwise Fst values for B. globosus in Additional file 1: Table S2). This suggests that the river system maybe acting as a barrier to outbreeding. However, genetic differentiation was very low, ranging between 0.004–0.061. Within the communities and river systems, evolutionary divergence ranged between 0–0.003 (Table 4).

Table 3 Nucleotide divergences among B. globosus populations in different communities
Table 4 Estimation of evolutionary divergence within localities

Haplotype, population genetic diversity and genetic divergence of B. truncatus

Bulinus truncatus populations were lower in number compared to B. globosus and they were only found in the Nyamaruru, Mupfure rivers and Kamoyo stream. Out of the 16 samples of B. truncatus analysed there were three unique cox1 haplotypes. The overall haplotype diversity was 0.575 with an estimated nucleotide diversity of 0.017. The average haplotype diversity for B. truncatus was also 0.017.

Within the river systems, haplotype diversity ranged between 0–0.286 (Table 2). All the three haplotypes were recorded in Mupfure River while the other two rivers exhibited a complete homogeneity (Table 2). The mean pairwise divergence among the haplotypes ranged from 0.017 to 0.034. The pairwise Fst value of B. truncatus was generated only for Mupfure River and Kamoyo stream, which had more than one sample (0.922). Nucleotide divergence between the Mupfure River and Kamoyo stream was 0.031 (Table 3).

Neutrality tests

Tajima’s D statistic and Fu’s Fs statistic showed no selection among the B. globosus populations among the sites. However, overall there were significant negative values for both Tajima’s D statistic and Fu’s Fs statistic. In the case of B. truncatus, overall, both neutrality tests were not significant, but Tajima’s D was significant for Mupfure River (Table 2).

Phylogenetic structuring

Bulinus globosus

Eleven of the 12 B. globosus haplotypes clustered closely with B. globosus from South Africa. H8 also had strong support for divergence away from the main haplotype group. One haplotype (H6) was an outlier and clustered more closely with a sample from Tanzania (Figs. 2, 3, 4). The topology of the three phylogenetic trees (Figs 2, 3, 4; ML, NJ and ME, respectively) was the same. H6 was found in two rivers (Nyamaruru and Mupfure), while H8 occurred in Mupfure River only. H4 was present in Mupfure River and Zvisokwe, a small stream branching into Mupfure River (Fig. 5). These water systems are also mapped in Fig. 1. The accession numbers corresponding to each sequence for all the haplotypes are found in Additional file 1: Table S3.

Fig. 2

Maximum likelihood phylogenetic analysis of the B. globosus and B. truncatus haplotypes collected in Madziwa, Zimbabwe, in relation to other Bulinus species. Each terminal bar is marked using the species name and the code given to the haplotype, e.g. H1 or T1. Bulinus globosus and B. truncatus haplotypes are represented by the letter H and T, respectively. The bootstrap values below 80% are not shown. There were a total of 607 positions in the final dataset. The accession numbers of the individual sequences forming the haplotypes are shown in Additional file 1: Table S3. The distribution of sequences among B. globosus haplotypes was as follows: haplotype 1 (n = 114 samples); haplotype 2 (n = 14); haplotype 3 (n = 11); haplotype 4 (n = 9); haplotype 5 (n = 8); haplotype 6 (n = 3); haplotype 8 (n = 2); haplotypes 7, 9, 10, 11 and 12 (n = 1 sample each). For B. truncatus, there were only three haplotypes. Among these three, the most common haplotype (T1) represented 10 samples, T2 representing 5 samples and T3 representing 1 sample

Fig. 3

Neighbor-joining phylogenetic tree of the B. globosus and B. truncatus haplotypes collected in Madziwa, Zimbabwe, in relation to other Bulinus species. Each terminal bar is marked using the species name and the code given to the haplotype, e.g. H1 or T1. Bulinus globosus and B. truncatus haplotypes are represented by the letter H and T, respectively. The bootstrap values below 80% are not shown. The accession numbers of the individual sequences forming the haplotypes are shown in Additional file 1: Table S3

Fig. 4

Minimum evolution phylogenetic tree of the B. globosus and B. truncatus haplotypes collected in Madziwa, Zimbabwe, in relation to other Bulinus species. Each terminal bar is marked using the species name and the code given to the haplotype, e.g. H1 or T1. Bulinus globosus and B. truncatus haplotypes are represented by the letter H and T, respectively. The significant bootstrap values for 1000 replicates are shown next to the branches. The accession numbers of the individual sequences forming the haplotypes are shown in Additional file 1: Table S3

Fig. 5

Detailed Maximum Likelihood phylogenetic tree of B. globosus and B. truncatus showing distribution of haplotypes by site, community and water system. Each terminal bar is marked using the species name, site number and the code given to the haplotype, e.g. H1 or T1. Bulinus globosus and B. truncatus haplotypes are represented by the letter H and T, respectively. The significant bootstrap values for 1000 replicates are shown next to the branches. The accession numbers of the individual sequences forming the haplotypes are shown in Additional file 1: Table S3

Bulinus truncatus

The three B. truncatus haplotypes clustered with data from previously characterized samples from Zimbabwe, but there was a considerable phylogenetic separation between the groups (Figs. 2, 3, 4). The topology of the three phylogenetic trees (Figs. 2, 3, 4; ML, NJ and ME, respectively) was the same. There was a strong geographical clustering as the clusters were made up of haplotypes from localities that were close together or from a single locality (Fig. 5). While T2 and T3 were found in Mupfure River, T1 was present in Kamoyo stream and Nyamaruru River. The accession numbers corresponding to each sequence for all the haplotypes are found in Additional file 1: Table S3.


The present results show that Bulinus population diversity exists in Madziwa, Zimbabwe. While the shell microsculpture clearly differed between the two species, as previously observed by Brown [8], identification can be confusing due to the twisting of the columella in B. truncatus and B. africanus-group. Within the B. africanus-group, it is difficult to clearly distinguish between B. africanus and B. globosus because the main morphologically variable feature, the penis, varies in size and shape due to parasitic infection, thus was not used in the present study to differentiate the two species. Likewise, as corroborated by Brown [8], it was difficult to construct a dichotomous key within the B. tropicus/truncatus group because of the continuous variation in shell shape between the species in this group.

Amplification of the cox1 shows that, within the species morphologically identified as B. globosus and B. truncatus, complete homogeneity for both was only recorded in Kamoyo stream, Nduna. In other sites, on average within the same site, approximately 25% of the samples had a different haplotype. Although the diversity recorded is low, it is expected within species. However, the diversity recorded in this study is greater than that recorded in Cameroon by Nalugwa et al. [36] in B. truncatus and B. forskalii populations where on average 16% of the snails sampled had a different haplotype. The ability of Bulinus spp. to self-fertilize and cross-breed induces some genetic consequences [37]. Bulinus globosus preferentially outbreeds increasing its diversity, therefore, higher genetic variability in the species is expected compared to that of B. truncatus, which preferentially self-fertilizes [38]. However, a previous study has shown that diversification is high in polyploid species such as B. truncatus [39].

Except for the dominant B. globosus cox1 haplotype, which spanned across almost all the sites where B. globosus was found, the diversity in haplotypes among the sites showed that genetic drift most likely due to bottleneck effects is a probable factor in determining diversity within and between species in this area. Our results showed homogeneity within the B. globosus originating from Nduna, Kamoyo stream, compared to the other communities. The complete homogeneity could be a result of natural selection, drift or perhaps a recent introduction of the intermediate host snail into these sites. Nevertheless, as presented by Schmid-hempel & Stauffer [39], host-parasite susceptibility can increase with loss of genetic variation; the snails in this river had the highest infection rates compared to the other rivers regardless of the fact that very few sites were sampled. On the contrary, a schistosomiasis prevalence study in the area has shown that Nduna community has a significantly low prevalence of the disease compared to Chihuri and Mupfure communities [40]. Both Mupfure and Chihuri had a higher genetic diversity of the intermediate host as compared to Nduna. Furthermore, the two areas had some infected snails at some point during the year. Thus, the contradicting scenario might be explained by other factors related to water contact behaviour of the human communities, given that schistosomiasis infection is a function of duration of exposure [41]. We may also hypothesize that the transmission of the disease in Nduna started relatively recently compared to other communities due to the absence of the intermediate host snail. Alternatively, it is highly probable that there is a risk of occurrence of schistosomiasis-transmitting genotypes becoming widely successful with selection of traits independent of parasitic infection [42]. The expanding population theory is corroborated by the significant test of neutrality for both Tajima’s D statistic and Fu’s Fs statistic for the whole B. globosus population or the B. truncatus population in Mupfure River. Similar results have also been reported by Zein-Eddine et al. [13].

Samples that were collected and analysed in this study, were from sites within close proximity to one another and they had an overall low genetic variation contrary to expectations due to limited gene flow between the populations. However, reproductive isolation, due to the different river systems, may have allowed independent evolution and divergence of these populations leading to formation of different haplotypes. The reproductive isolation might result in complete homogeneity in Kamoyo stream and the presence of ‘unique’ haplotypes in some of the rivers. Like in the previous study by Standley et al. [43] for Biomphalaria choanomphala in Lake Victoria where the ‘unique’ haplotypes were locally restricted to the central and western regions of the lakeshore in Uganda, in the present study, the ‘unique’ haplotypes are also restricted to Nyamaruru and Mupfure rivers.

To our knowledge, the present study provides the first genetic assessment of combined data for B. globosus and B. truncatus from Zimbabwe, where urogenital schistosomiasis is endemic. However, more sophisticated statistical analyses of these haplotypes should be performed to examine the levels of divergence which may lead to species delimitation.


A detailed understanding of genetic variability within intermediate host snail species is important for understanding disease epidemiology. A novel insight from our study is the mitochondrial cox1 variability within B. globosus and additional B. truncatus haplotypes in Madziwa area, Zimbabwe. Further investigations using new tools for detection of infection will help investigate the snail populations involved in transmission. The identification of Bulinus in these areas will support targeted control efforts for schistosomiasis in the study area. Further studies over a large spatial area with higher sample size are also required to understand the genetic demographic patterns of these intermediate host snail species.

Availability of data and materials

Data supporting the conclusions of this article are included within the article and its additional files. The newly generated sequences were submitted to the GenBank database under the accession numbers MN397785–MN397822. The datasets analysed during the present study are available from the corresponding author upon reasonable request.



cytochrome oxidase subunit 1


  1. 1.

    Gryseels B, Polman K, Clerinx J, Kestens L. Human schistosomiasis. Lancet. 2006;368:1106–18.

    PubMed  Article  Google Scholar 

  2. 2.

    Rollinson D, Stothard JR, Southgate VR. Interactions between intermediate snail hosts of the genus Bulinus and schistosomes of the Schistosoma haematobium group. Parasitology. 2001;123(Suppl):245–60.

    Article  Google Scholar 

  3. 3.

    Kariuki HC, Clennon JA, Brady MS, Kitron U, Sturrock RF, Ouma JH, et al. Distribution patterns and cercarial shedding of Bulinus nasutus and other snails in the Msambweni area, Coast Province, Kenya. Am J Trop Med Hyg. 2004;70:449–56.

    PubMed  Article  Google Scholar 

  4. 4.

    Allan F, Dunn AM, Emery AM, Stothard JR, Johnston DA, Kane RA, et al. Use of sentinel snails for the detection of Schistosoma haematobium transmission on Zanzibar and observations on transmission patterns. Acta Trop. 2013;128:234–40.

    PubMed  Article  Google Scholar 

  5. 5.

    Opisa S, Odiere MR, Jura WGZO, 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.

    PubMed  PubMed Central  Article  Google Scholar 

  6. 6.

    Woolhouse ME, Chandiwana SK. Temporal patterns in the epidemiology of schistosome infections of snails: a model for field data. Parasitology. 1990;100:247–53.

    PubMed  Article  Google Scholar 

  7. 7.

    Woolhouse ME, Chandiwana SK. Spatial and temporal heterogeneity in the population dynamics of Bulinus globosus and Biomphalaria pfeifferi and in the epidemiology of their infection with schistosomes. Parasitology. 1989;98:21–34.

    PubMed  Article  Google Scholar 

  8. 8.

    Brown D. Freshwater snails of Africa and their medical importance. 2nd ed. London: Taylor and Francis; 1994.

    Google Scholar 

  9. 9.

    Makura O, Kristensen TK. National freshwater snail survey of Zimbabwe. In: The tenth international malacological congress. Tübingen: University of Tübingen Press; 1991. p. 227–232.

  10. 10.

    Jarne P, Finot L, Delay B, Thaler L. Self-fertilization versus cross-fertilization in the hermaphroditic freshwater snail Bulinus globosus. Evolution. 1991;45:1136–46.

    PubMed  Article  Google Scholar 

  11. 11.

    Jarne P, Charlesworth D. The evolution of the selfing rate in functionally hermaphrodite plants and animals. Annu Rev Ecol Syst. 1993;24:441–66.

    Article  Google Scholar 

  12. 12.

    Viard F, Bremond P, Labbo R, Justy F, Delay B, Jarne P. Microsatellites and the genetics of highly selfing populations in the freshwater snail Bulinus truncatus. Genetics. 1996;142:1237–47.

    CAS  PubMed  PubMed Central  Google Scholar 

  13. 13.

    Zein-Eddine R, Djuikwo-Teukeng FF, Al-Jawhari M, Senghor B, Huyse T, Dreyfuss G. Phylogeny of seven Bulinus species originating from endemic areas in three African countries, in relation to the human blood fluke Schistosoma haematobium. BMC Evol Biol. 2014;14:271.

    PubMed  PubMed Central  Article  Google Scholar 

  14. 14.

    Zein-Eddine R, Djuikwo-Teukeng FF, Darc Y, Dreyfussa G, Van den Broeckd F. Population genetics of the Schistosoma snail host Bulinus truncatus in Egypt. Acta Trop. 2017;172:36–43.

    PubMed  Article  Google Scholar 

  15. 15.

    Charbonnel N, Angers B, Rasatavonjizay R, Remond PB, Debain C, Jarne P. Evolutionary aspects of the metapopulation dynamics of Biomphalaria pfeifferi, the intermediate host of Schistosoma mansoni. J Evol Biol. 2002;15:248–61.

    Article  Google Scholar 

  16. 16.

    Thiele AE, Corrêa-Oliveira G, Gazzinelli A, Minchella JD. Elucidating the temporal and spatial dynamics of Biomphalaria glabrata genetic diversity in three Brazilian villages. Trop Med Int Health. 2013;18:1164–73.

    PubMed  PubMed Central  Article  Google Scholar 

  17. 17.

    Stothard JR, Rollinson D. Partial DNA sequences from the mitochondrial cytochrome oxidase subunit I (COI) gene can differentiate the intermediate snail hosts Bulinus globosus and B. nasutus (Gastropoda: Planorbidae). J Nat Hist. 1997;31:727–37.

    Article  Google Scholar 

  18. 18.

    Kane RA, Stothard JR, Emery AM, Rollinson D. Molecular characterization of freshwater snails in the genus Bulinus: a role for barcodes? Parasit Vectors. 2008;1:15.

    PubMed  PubMed Central  Article  CAS  Google Scholar 

  19. 19.

    Abe EM, Guo YH, Shen H, Mutsaka-Makuvaza MJ, Habib MR, Xue JB, et al. Phylogeography of Bulinus truncatus (Audouin, 1827) (Gastropoda: Planorbidae) in selected African countries. Trop Med Infect Dis. 2018;3:E127.

    PubMed  Article  Google Scholar 

  20. 20.

    Mandahl-Barth G. The species of the genus Bulinus, intermediate hosts of Schistosoma. Bull World Health Organ. 1965;33:33–44.

    CAS  PubMed  PubMed Central  Google Scholar 

  21. 21.

    Mukaratirwa S, Sieglsmund HR, Kristensen TK, Chandiwana SK. Population genetics and genetic variability of Bulinus globosus (Gastropoda: Planorbidae) from the two main river systems in Zimbabwe. J Hered. 1996;87:288–94.

    Article  Google Scholar 

  22. 22.

    Mukaratirwa S, Siegismund HR, Kristensen TK, Chandiwana SK. Genetic structure and parasite compatibility of Bulinus globosus (Gastropoda: Planorbidae) from two areas of different endemicity of Schistosoma haematobium in Zimbabwe. Int J Parasitol. 1996;26:269–80.

    CAS  PubMed  Article  Google Scholar 

  23. 23.

    Davies CM, Webster JP, Krüger O, Munatsi A, Ndamba J, Woolhouse ME. Host-parasite population genetics: a cross-sectional comparison of Bulinus globosus and Schistosoma haematobium. Parasitology. 1999;119:295–302.

    PubMed  Article  Google Scholar 

  24. 24.

    Stothard JR, Hughes S, Rollinson D. Variation within the internal transcribed spacer (ITS) of ribosomal DNA genes of intermediate snail hosts within the genus Bulinus (Gastropoda: Planorbidae). Acta Trop. 1996;61:19–29.

    CAS  PubMed  Article  Google Scholar 

  25. 25.

    Mutsaka-Makuvaza MJ, Matsena-Zingoni Z, Tshuma C, Ray S, Zhou XN, Webster B, et al. Reinfection of urogenital schistosomiasis in pre-school children in a highly endemic district in northern Zimbabwe: compliance assessment in a 12 months longitudinal study. Infect Dis Poverty. 2018;7:102.

    PubMed  PubMed Central  Article  Google Scholar 

  26. 26.

    WHO. Field use of molluscicides in schistosomiasis control programmes: an operational manual for programme managers. Geneva: World Health Organization; 2017.

    Google Scholar 

  27. 27.

    Frandsen F, Christensen NØ. An introductory guide to the identification of cercariae from African freshwater snails with special reference to cercariae of trematode species of medical and veterinary importance. Acta Trop. 1984;41:181–202.

    CAS  PubMed  Google Scholar 

  28. 28.

    Altschul SF, Gish W, Miller W, Myers EW, Lipman DJ. Basic local alignment search tool. J Mol Biol. 1990;215:403–10.

    CAS  PubMed  Article  Google Scholar 

  29. 29.

    Tamura K, Stecher G, Peterson D, Filipski A, Kumar S. MEGA6: molecular evolutionary genetics analysis version 6.0. Mol Biol Evol. 2013;30:2725–9.

    CAS  PubMed  PubMed Central  Article  Google Scholar 

  30. 30.

    Thompson JD, Higgins DG, Gibson T. CLUSTAL W: improving the sensitivity of progressive multiple sequence alignment through sequence weighting, position-specific gap penalties and weight matrix choice. Nucleic Acids Res. 1994;22:4673–80.

    CAS  PubMed  PubMed Central  Article  Google Scholar 

  31. 31.

    Tamura K, Nei M. Estimation of the number of nucleotide substitutions in the control region of mitochondrial DNA in humans and chimpanzees. Mol Biol Evol. 1993;10:512–26.

    CAS  PubMed  Google Scholar 

  32. 32.

    Felsenstein J. Confidence limits on phylogenies: an approach using the bootstrap. Evolution. 1985;39:783–91.

    PubMed  Article  Google Scholar 

  33. 33.

    Jukes TH, Cantor CR. Evolution of protein molecules. In: Munro HN, editor. Mammalian protein metabolism. New York: Academic Press; 1969. p. 21–132.

    Google Scholar 

  34. 34.

    Rozas J, Ferrer-Mata A, Sánchez-DelBarrio JC, Guirao-Rico S, Librado P, Ramos-Onsins SE, et al. DnaSP 6: DNA sequence polymorphism analysis of large data sets. Mol Biol Evol. 2017;34:3299–302.

    CAS  PubMed  Article  Google Scholar 

  35. 35.

    Tajima F. Statistical methods for testing the neutral mutation hypothesis by DNA polymorphism. Genetics. 1989;123:585–95.

    CAS  PubMed  PubMed Central  Google Scholar 

  36. 36.

    Nalugwa A, Jørgensen A, Nyakaana S, Kristensen TK. Molecular phylogeny of Bulinus (Gastropoda: Planorbidae) reveals the presence of three species complexes in the Albertine Rift freshwater bodies. Int J Genet Mol Biol. 2010;2:130–9.

    CAS  Google Scholar 

  37. 37.

    Jarne P. Resistance genes at the population level. Parasitol Today. 1993;9:216–7.

    CAS  PubMed  Article  Google Scholar 

  38. 38.

    Brown DS, Shaw KM. Freshwater snails of the Bulinus truncatus/tropicus complex in Kenya: tetraploid species. J Mollus Stud. 1989;55:509–32.

    Article  Google Scholar 

  39. 39.

    Schmid-hempel S, Stauffer HP. Parasites and flower choice of bumblebees. Anim Behav. 1998;55:819–25.

    CAS  PubMed  Article  PubMed Central  Google Scholar 

  40. 40.

    Mutsaka-Makuvaza MJ, Matsena-Zingoni Z, Katsidzira A, Tshuma C, Chin’ombe N, Zhou XN, et al. Urogenital schistosomiasis and risk factors of infection in mothers and preschool children in an endemic district in Zimbabwe. Parasit Vectors. 2019;12:427.

    PubMed  PubMed Central  Article  Google Scholar 

  41. 41.

    Dalton PR, Pole D. Water-contact patterns in relation to Schistosoma haematobium infection. Bull World Health Organ. 1978;56:417–26.

    CAS  PubMed  PubMed Central  Google Scholar 

  42. 42.

    Campbell G, Noble LR, Rollinson D, Southgate VR, Webster JP, Jones CS. Low genetic diversity in a snail intermediate host (Biomphalaria pfeifferi Krass, 1848) and schistosomiasis transmission in the Senegal River Basin. Mol Ecol. 2010;19:241–56.

    CAS  PubMed  Article  Google Scholar 

  43. 43.

    Standley CJ, Goodacre SL, Wade CM, Stothard JR. The population genetic structure of Biomphalaria choanomphala in Lake Victoria, East Africa: implications for schistosomiasis transmission. Parasit Vectors. 2014;7:524.

    PubMed  PubMed Central  Article  Google Scholar 

Download references


Our appreciation goes to the Madziwa residents and their leaders who allowed us to work in the area and also assisted in identification of water contact sites.


Sample collection was funded by individuals while laboratory analysis was funded by National Institute of Parasitic Diseases, Shanghai.

Author information




MJMM, NM, CT and XNZ designed the study protocol. MJMM and NM collected specimens and GPS data. MJMM, EA and NM performed laboratory analysis of specimens. MJMM, BW, JM and TM performed the analysis. MJMM wrote the initial draft of the manuscript. NM, BW, NC, EA and FA reviewed the manuscript. All authors read and approved the final manuscript.

Corresponding author

Correspondence to Nicholas Midzi.

Ethics declarations

Ethics approval and consent to participate

Prior to the commencement of the study, approval was obtained from the community leaders in the area. Ethical clearance of the study was sought from the Joint Research Ethics Committee (JREC) for the University of Zimbabwe, College of Health Sciences and the Parirenyatwa Group of Hospitals (JREC/251/16). Further approval was sought from the national ethical review board, the Medical Research Council of Zimbabwe (MRCZ) (MRCZ/A/2149). The communities were sensitized about the nature of the study and study objectives in local language (Shona).

Consent for publication

Not applicable.

Competing interests

The authors declare that they have no competing interests.

Additional information

Publisher's Note

Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Supplementary information

Additional file 1: Table S1.

List of sequences drawn from the Genbank incorporated for phylogenetic analysis. Table S2. Fst values for B. globosus populations. Table S3. Accession numbers of individual sequences for the haplotypes recorded in Madziwa, Zimbabwe.

Rights and permissions

Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit The Creative Commons Public Domain Dedication waiver ( applies to the data made available in this article, unless otherwise stated in a credit line to the data.

Reprints and Permissions

About this article

Verify currency and authenticity via CrossMark

Cite this article

Mutsaka-Makuvaza, M.J., Zhou, X., Tshuma, C. et al. Molecular diversity of Bulinus species in Madziwa area, Shamva district in Zimbabwe: implications for urogenital schistosomiasis transmission. Parasites Vectors 13, 14 (2020).

Download citation


  • Bulinus globosus
  • Bulinus truncatus
  • Cox1
  • Diversity
  • Phylogenetics
  • Zimbabwe