Breeding system and demography shape population genetic structure across ecological and climatic zones in the African freshwater snail, Bulinus forskalii (Gastropoda, Pulmonata), intermediate host for schistosomes

Large-scale and small-scale population genetic structure of the medically important gastropod species Bulinus truncatus (Gastropoda, Heterobranchia)

Background

Gastropod snails remain strongly understudied, despite their important role in transmitting parasitic diseases. Knowledge of their distribution and population dynamics increases our understanding of the processes driving disease transmission. We report the first study to use high-throughput sequencing (HTS) to elucidate the population genetic structure of the hermaphroditic snail Bulinus truncatus (Gastropoda, Heterobranchia) on a regional (17–150 km) and inter-regional (1000–5400 km) scale. This snail species acts as an intermediate host of Schistosoma haematobium and Schistosoma bovis, which cause human and animal schistosomiasis respectively.

Methods

Bulinus truncatus snails were collected in Senegal, Cameroon, Egypt and France and identified through DNA barcoding. A single-end genotyping-by-sequencing (GBS) library, comprising 87 snail specimens from the respective countries, was built and sequenced on an Illumina HiSeq 2000 platform. Reads were mapped against S. bovis and S. haematobium reference genomes to identify schistosome infections, and single nucleotide polymorphisms (SNPs) were scored using the Stacks pipeline. These SNPs were used to estimate genetic diversity, assess population structure and construct phylogenetic trees of B. truncatus.

Results

A total of 10,750 SNPs were scored and used in downstream analyses. The phylogenetic analysis identified five clades, each consisting of snails from a single country but with two distinct clades within Senegal. Genetic diversity was low in all populations, reflecting high selfing rates, but varied between locations due to habitat variability. Significant genetic differentiation and isolation by distance patterns were observed at both spatial scales, indicating that gene flow is not strong enough to counteract the effects of population bottlenecks, high selfing rates and genetic drift. Remarkably, the population genetic differentiation on a regional scale (i.e. within Senegal) was as large as that between populations on an inter-regional scale. The blind GBS technique was able to pick up parasite DNA in snail tissue, demonstrating the potential of HTS techniques to further elucidate the role of snail species in parasite transmission.

Conclusions

HTS techniques offer a valuable toolbox to further investigate the population genetic patterns of intermediate schistosome host snails and the role of snail species in parasite transmission.

Graphical Abstract

Supplementary Information

The online version contains supplementary material available at 10.1186/s13071-022-05445-x.

Population genetics of Oncomelania hupensis snails, intermediate hosts of Schistosoma japonium, from emerging, re-emerging or established habitats within China

Schistosomiasis remains one of the world's most significant neglected tropical diseases, second only to malaria in terms of socioeconomic impact. In 2014, China proposed the goal of schistosomiasis japonicum elimination by 2025. However, one major challenge is the widely distributed, and in certain cases potentially increasing, habitats of Oncomelania hupensis, the snail intermediate hosts of S. japonicum. Therefore, an understanding of population genetics of O. hupensis in new or re-emerged habitats, together with that of the established habitats with snail persistence, would be valuable in controlling and predicting the future transmission dynamics of schistosomiasis in China. Using nine microsatellite loci, we conducted population genetic analyses of snails sampled from one habitat where snails were detected for the first time, one (previously eliminated) habitat with re-emerged snails, and one habitat with established snail persistence. Results showed lower diversities, in terms of number of observed alleles per locus (Na), number of effective alleles per locus (NeA), observed (Ho) and expected heterozygosity (He), in snails from new or re-emerged snail habitats than from the habitat with snail persistence. The smallest effective population size was inferred in the re-emerged snail habitat, but the largest was in the new habitat rather than in the habitat with snail persistence. No bottleneck effects were detected in new or re-merged habitats. No or low sub-structure was inferred in new and persistent snail habitats. Snails from the three sites were clearly separated and low gene flow was estimated between sites. We propose that snails at the new habitat may have been introduced through immigration, whereas snails at the re-emerged habitat may be the consequence of those few snails remaining subsequently expanding through reproduction. We discuss our results in terms of their theoretical and applied implications.

Riverscape genetic variation, migration patterns, and morphological variation of the threatened Round Rocksnail, Leptoxis ampla

Within riverine systems, headwater populations are hypothesized to harbour higher amounts of genetic distinctiveness than populations in the main stem of a river and display increased genetic diversity in large, downstream habitats. However, these hypotheses were mostly developed with insects and fish, and they have not been tested on many invertebrate lineages. Pleuroceridae gastropods are of particular ecological importance to rivers of eastern North America, sometimes comprising over 90% of macroinvertebrate biomass. Yet, virtually nothing is known of pleurocerid landscape genetics, including whether genetic diversity follows predictions made by hypotheses developed on more mobile species. Moreover, the commonly repeated hypothesis that intraspecific morphological variation in gastropods results from ecophenotypic plasticity has not been well tested on pleurocerids. Using 2bRAD-seq to discover single nucleotide polymorphisms, we show that the threatened, Cahaba River endemic pleurocerid, Leptoxis ampla, has limited gene flow among populations and that migration is downstream-biased, conflicting with previous hypotheses. Both tributary and main stem populations harbour unique genomic profiles, and genetic diversity was highest in downstream populations. Furthermore, L. ampla shell morphology was more correlated with genetic differences among individuals and populations than habitat characteristics. We anticipate similar genetic and demographic patterns to be seen in other pleurocerids, and hypotheses about gene flow and population demographics that were based on more mobile taxa often, but not always, apply to freshwater gastropods. From a conservation standpoint, genetic structure of L. ampla populations suggests distinctive genetic diversity is lost with localized extirpation, a phenomenon common across the range of Pleuroceridae.

Population genetics of the Schistosoma snail host Bulinus truncatus in Egypt

The tropical freshwater snail Bulinus truncatus serves as an important intermediate host of several human and cattle Schistosoma species in many African regions. Despite some ecological and malacological studies, there is no information on the genetic diversity of B. truncatus in Egypt. Here, we sampled 70-100 snails in ten localities in Upper Egypt and the Nile Delta. Per locality, we sequenced 10 snails at a partial fragment of the cytochrome c oxidase subunit 1 gene (cox1) and we genotyped 25-30 snails at six microsatellite markers. A total of nine mitochondrial haplotypes were detected, of which five were unique to the Nile Delta and three were unique to Upper Egypt, indicating that snail populations may have evolved independently in both regions. Bayesian clustering and hierarchical F-statistics using microsatellite markers further revealed strong population genetic structure at the level of locality. Observed heterozygosity was much lower compared to what is expected under random mating, which could be explained by high selfing rates, population size reductions and to a lesser extent by the Wahlund effect. Despite these observations, we found signatures of gene flow and cross-fertilization, even between snails from the Nile Delta and Upper Egypt, indicating that B. truncatus can travel across large distances in Egypt. These observations could have serious consequences for disease epidemiology, as it means that infected snails from one region could rapidly and unexpectedly spark a new epidemic in another distant region. This could be one of the factors explaining the rebound of human Schistosoma infections in the Nile Delta, despite decades of sustained schistosomiasis control.

Urogenital schistosomiasis and soil-transmitted helminthiasis (STH) in Cameroon: An epidemiological update at Barombi Mbo and Barombi Kotto crater lakes assessing prospects for intensified control interventions

BACKGROUND

The crater lakes of Barombi Mbo and Barombi Kotto are well-known transmission foci of schistosomiasis and soil-transmitted helminthiasis having had several important control initiatives previously. To collect contemporary epidemiological information, a cross-sectional survey was undertaken inclusive of: signs and symptoms of disease, individual treatment histories, local water, sanitation and hygiene (WASH)-related factors and malacological surveillance, with molecular characterisation of specimens.

METHODS

At each lake, a community cross-sectional survey was undertaken using a combination of stool and urine parasitological sampling, and interview with pro-forma questionnaires. A total of 338 children and adults participated. Material from snail and parasite species were characterised by DNA methods.

RESULTS

Egg-patent prevalence of urogenital schistosomiasis was 8.7% at Barombi Mbo (all light-intensity infections) and 40.1% at Barombi Kotto (21.2% heavy-intensity infections). Intestinal schistosomiasis was absent. At Barombi Kotto, significantly more women reported signs and symptoms associated with female genital schistosomiasis. While there had been extensive recent improvement in WASH-related infrastructure at Barombi Mbo, water contact risk scores were higher among schistosomiasis-infected participants (P < 0.001) and at Barombi Kotto in general (P < 0.001). Across both lakes, mean prevalence of STH was very low (6.3%) evidencing an impressive decrease of 79.0% over the last decade; neither Strongyloides stercoralis nor Ascaris lumbricoides were found. A total of 29 freshwater sampling sites were inspected for snails, 13 in Barombi Mbo and 16 in Barombi Kotto; water chemistry differed significantly (P < 0.0001) between lakes for both mean pH (7.9 v. 9.6) and mean conductivity (64.3 μS v. 202.1 μS) respectively. Only two Bulinus camerunensis found on the central island of Barombi Kotto were observed to shed schistosome cercariae, but schistosome DNA was later detected in Bulinus sampled from both lakes as well as in Indoplanorbis exustus, an invasive species from Asia.

CONCLUSIONS

STH is currently at very low levels while urogenital schistosomiasis is of greatest concern at Barombi Kotto. This assessment highlights a unique opportunity for further study of the epidemiological dynamics at these crater lakes, to evaluate future intensified interventions both in terms of gaining and sustaining control at Barombi Kotto or in moving towards local interruption of transmission of both diseases at Barombi Mbo.

Population genetic structure of the freshwater snail, Bulinus globosus, (Gastropoda: Planorbidae) from selected habitats of KwaZulu-Natal, South Africa

The freshwater snail Bulinus globosus is an important intermediate host of Schistosoma haematobium, the causative agent of urinary schistosomiasis. This disease is of major health concern, especially in Africa where the majority of cases have been reported. In this study the inter- and intra-genetic diversity and population genetic structure of B. globosus from nine locations in the KwaZulu-Natal province of South Africa was studied using four polymorphic microsatellite loci (BgZ1-BgZ4). Moderate genetic diversity was detected within populations with a mean diversity (HE) of 0.49±0.09. The majority of populations significantly deviated from Hardy-Weinberg equilibrium (p<0.05), due to a deficit of heterozygotes. Such deviations may be due to founder events that were caused by bottlenecks that occurred as a result of frequent droughts and flooding that these snails' habitats are exposed to. Overall, the populations studied seem to be partially inbreeders/selfers with mean estimates of 0.24/0.38. A discernable genetic structure was elucidated among populations as evident by the mean pairwise FST of 0.58±0.13. There was no significant association between genetic and geographical distance among populations, an indication of limited gene flow. This increases the chances of populations losing alleles due to genetic drift. Populations in close proximity demonstrated high genetic differentiation (58.77% total variation) due to allelic differences between them. The sample populations fell into 12 clusters, however, the populations from uMkhanyakude and uThungulu exhibited no discernable genetic structure. Genetically, the Bhobhoyi site found within the uGu district was equidistant to the two main sampling regions.

Analysis of genetic diversity and population structure of Bellamya quadrata from lakes of middle and lower Yangtze River

As an endemic species of freshwater gastropods in China, Bellamya quadrata plays an important role in ecosystem service provision and commercial importance. However, the species is overharvested and its natural habitats are under severe threat due to fragmentation and loss. To estimate the genetic diversity and population structure of B. quadrata, 285 individuals from eight lake populations across middle and lower Yangtze River were sampled. Seven microsatellite loci were genotyped. Our results showed that (i) the genetic diversity of B. quadrata was high in most of the studied populations, yet effective population sizes appear to be rather small in some populations; (ii) low levels of genetic differentiation exists among populations but gene flow was generally high; (iii) no clear geographic or genetic structure was observed in the studied region, implying mechanisms (zoochoric dispersal and anthropogenic translocations) that enhance dispersal and gene flow have promoted population connectivity. However, the comparatively high genetic diversity of B. quadrata could be attributed to a lag phase, suggesting that the genetic diversity of this species may be lost in the future and the priorities for conservation of B. quadrata are necessary.

Phylogeny of seven Bulinus species originating from endemic areas in three African countries, in relation to the human blood fluke Schistosoma haematobium

Background

Snails species belonging to the genus Bulinus (Planorbidae) serve as intermediate host for flukes belonging to the genus Schistosoma (Digenea, Platyhelminthes). Despite its importance in the transmission of these parasites, the evolutionary history of this genus is still obscure. In the present study, we used the partial mitochondrial cytochrome oxidase subunit I (cox1) gene, and the nuclear ribosomal ITS, 18S and 28S genes to investigate the haplotype diversity and phylogeny of seven Bulinus species originating from three endemic countries in Africa (Cameroon, Senegal and Egypt).

Results

The cox1 region showed much more variation than the ribosomal markers within Bulinus sequences. High levels of genetic diversity were detected at all loci in the seven studied species, with clear segregation between individuals and appearance of different haplotypes, even within same species from the same locality. Sequences clustered into two lineages; (A) groups Bulinus truncatus, B. tropicus, B. globosus and B. umbilicatus; while (B) groups B. forskalii, B. senegalensis and B. camerunensis. Interesting patterns emerge regarding schistosome susceptibility: Bulinus species with lower genetic diversity are predicted to have higher infection prevalence than those with greater diversity in host susceptibility.

Conclusion

The results reported in this study are very important since a detailed understanding of the population genetic structure of Bulinus is essential to understand the epidemiology of many schistosome parasites.

Eutrophication and Dreissena invasion as drivers of biodiversity: a century of change in the mollusc community of Oneida Lake

Changes in nutrient loading and invasive species are among the strongest human-driven disturbances in freshwater ecosystems, but our knowledge on how they affect the biodiversity of lakes is still limited. We conducted a detailed historical analysis of the mollusc community of Oneida Lake based on our comprehensive lakewide study in 2012 and previous surveys dating back to 1915. In the early 20th century, the lake had a high water clarity, with abundant macrophytes and benthic algae, and hosted the most diverse molluscan community in New York State, including 32 gastropod and 9 unionid species. By the 1960s, lake turbidity increased during a period of anthropogenic eutrophication, resulting in a 38% decline in species richness and a 95% reduction in abundance of native gastropods grazing on benthic algae. Following the invasion of Dreissena spp. in 1991 and subsequent increases in water clarity, native gastropod species richness expanded by 37% and abundance increased 20-fold by 2012. In contrast, filter-feeding unionids were unaffected by increased turbidity during the period of eutrophication but were extirpated by dreissenids. Through contrasting effects on turbidity, eutrophication and Dreissena spp. have likely driven the observed changes in native grazing gastropods by affecting the abundance of light-limited benthic algae. Given the high species richness and ecological importance of benthic grazers, monitoring and managing turbidity is important in preserving molluscan diversity.

Sympatric and allopatric experimental infections of the planorbid snail Gyraulus chinensis with miracidia of Euparyphium albuferensis (Trematoda: Echinostomatidae)

An experimental infection with echinostomatid miracidia in sympatric or 'local' vs. allopatric or 'away' snail combinations, as a model to examine parasite compatibility, was carried out. We employed Euparyphium albuferensis miracidia to infect Gyraulus chinensis snails, from three different natural parks: Albufera (Valencia, Spain); the Ebro Delta (Tarragona, Spain) and Coto de Doñana (Huelva, Spain). Insignificant differences between the three snail strains were noted for the infection rate and the rhythm of daily cercarial production. However, a significantly higher total cercarial production per snail, patent period and life span were observed in local snails. The different infection characteristics in the three G. chinensis strains considered reveal that E. albuferensis miracidia demonstrate local adaptation.

Low genetic diversity in a snail intermediate host (Biomphalaria pfeifferi Krass, 1848) and schistosomiasis transmission in the Senegal River Basin

Population genetic perturbations of intermediate hosts, often a consequence of human pressure on environmental resources, can precipitate unexpectedly severe disease outbreaks. Such disturbances are set to become increasingly common following range changes concomitant with climate shifts, dwindling natural resources and major infrastructure changes such as hydroprojects. Construction of the Diama dam in the Senegal River Basin (SRB) reduced river salinity, enabling the freshwater snail intermediate host Biomphalaria pfeifferi to rapidly expand its distribution. A serious public health problem ensued, with an epidemic of intestinal schistosomiasis occurring in the previously schistosome-free Richard-Toll region within 2 years. The current study aimed to assess the population variability of B. pfeifferi in the SRB, and speculate upon its subsequent impact on host-parasite interactions following such engineered ecological change. Genetic variation at nine polymorphic microsatellite loci revealed little population differentiation in SRB snails compared with those from natural habitats in Zimbabwe, where Schistosoma mansoni transmission is much lower. 'Open' SRB habitats are associated with greater water contact, smaller population sizes and less genetic diversity, with sites downstream of Richard-Toll showing greater inter- and intrapopulation variation, concomitant with less frequent human contact. These observations may be explained by rapid expansion into pristine habitat selecting for high fecundity genotypes at the expense of schistosome resistance, presenting S. mansoni with genetically homogenous highly fecund susceptible populations around the focal point, promoting development of a highly compatible host-parasite relationship. Longitudinal study of such systems may prove important in predicting public health risks engendered by future environmental engineering projects.