Stirred, not shaken: genetic structure of the intermediate snail host Oncomelania hupensis robertsoni in an historically endemic schistosomiasis area

The genetic diversity of Oncomelania hupensis robertsoni, intermediate hosts of Schistosoma japonicum in hilly regions of China, using microsatellite markers

BACKGROUND

The elimination of schistosomiasis remains a challenging task, with current measures primarily focused on the monitoring and control of Oncomelania hupensis (O. hupensis) snail, the sole intermediate host of Schistosome japonicum. Given the emerging, re-emerging, and persistent habitats of snails, understanding their genetic diversity might be essential for their successful monitoring and control. The aims of this study were to analyze the genetic diversity of Oncomelania hupensis robertsoni (O. h. robertsoni) using microsatellite DNA markers; and validate the applicability of previously identified microsatellite loci for O. hupensis in hilly regions.

METHODS

A total of 17 populations of O. h. robertsoni from Yunnan Province in China were selected for analysis of genetic diversity using six microsatellite DNA polymorphic loci (P82, P84, T4-22, T5-11, T5-13, and T6-27).

RESULTS

The number of alleles among populations ranged from 0 to 19, with an average of 5. The average ranges of expected (He) and observed (Ho) heterozygosity within populations were 0.506 to 0.761 and 0.443 to 0.792, respectively. The average fixation index within the population ranged from - 0.801 to 0.211. The average polymorphic information content (PIC) within the population ranged from 0.411 to 0.757, appearing to be polymorphic for all loci (all PIC > 0.5), except for P28 and P48. A total of 68 loci showed significant deviations from Hardy-Weinberg equilibrium (P < 0.05), and pairwise Fst values ranged from 0.051 to 0.379. The analysis of molecular variance indicated that 88% of the variation occurred within snail populations, whereas 12% occurred among snail populations. Phylogenetic trees and principal coordinate analysis revealed two distinct clusters within the snail population, corresponding to "Yunnan North" and "Yunnan South".

CONCLUSIONS

O. h. robertsoni exhibited a relatively high level of genetic differentiation, with variation chiefly existing within snail populations. All snail in this region could be separated into two clusters. The microsatellite loci P82 and P84 might not be suitable for classification studies of O. hupensis in hilly regions. These findings provided important information for the monitoring and control of snail, and for further genetic diversity studies on snail populations.

Population Genetics of Oncomelania hupensis Snails from New-Emerging Snail Habitats in a Currently Schistosoma japonicum Non-Endemic Area

Schistosomiasis is still one of the most significant neglected tropical diseases worldwide, and China is endemic for Schistosoma japonicum. With its great achievement in schistosomiasis control, the government of China has set the goal to eliminate the parasitic disease at the country level by 2030. However, one major challenge is the remaining huge areas of habitats for the intermediate host Oncomelania hupensis. This is further exacerbated by an increasing number of new emerging snail habitats reported each year. Therefore, population genetics on snails in such areas will be useful in evaluation of snail control effect and/or dispersal. We then sampled snails from new emerging habitats in Taicang of Jiangsu, China, a currently S. japonicum non-endemic area from 2014 to 2017, and performed population genetic analyses based on nine microsatellites. Results showed that all snail populations had low genetic diversity, and most genetic variations originated from within snail populations. The estimated effective population size for the 2015 population was infinitive. All snails could be separated into two clusters, and further DIYABC analysis revealed that both the 2016 and the 2017 populations may derive from the 2015, indicating that the 2017 population must have been missed in the field survey performed in 2016. These findings may have implications in development of more practical guidelines for snail monitoring and control.

Genetic diversity and structure of Oncomelania hupensis hupensis in two eco-epidemiological settings as revealed by the mitochondrial COX1 gene sequences

BACKGROUND

Oncomelania hupensis hupensis is the only intermediate host of Schistosoma japonicum, the causative agent of schistosomiasis in China and is therefore of significant medical and veterinary health importance. Although tremendous progress has been achieved, there remains an understudied area of approximately 2.06 billion m² of potential snail habitats. This area could be further increased by annual flooding. Therefore, an understanding of population genetics of snails in these areas may be useful for future monitoring and control activities.

METHODS AND RESULTS

We sampled snails from Hexian (HX), Zongyang (ZY) and Shitai (ST) in Anhui (schistosomiasis transmission control), and from Hengtang (HT), Taicang (TC), Dongsan (DS) and Xisan (XS) in Jiangsu (schistosomiasis transmission interrupted), downstream of Anhui. ST, DS and XS are classified as hilly and mountainous areas, and HX, ZY, TC and HT as lake and marshland areas. The mitochondrial cytochrome c oxidase subunit I gene were sequenced. Out of 115 snails analyzed, 29 haplotypes were identified. We observed 56 (8.72%) polymorphic sites consisting of 51 transitions, four transversions and one multiple mutational change. The overall haplotype and nucleotide diversity were 0.899 and 0.01569, respectively. Snail populations in Anhui had higher genetic diversity than in Jiangsu. 73.32% of total variation was distributed among sites and 26.68% within sites. Snails were significantly separated according to eco-epidemiological settings in both network and phylogenetic analyses.

CONCLUSION

Our results could provide important guidance towards assessing coevolutionary interactions of snails with S. japonicum, as well as for future molluscan host monitoring and control activities.

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.

Genetic Evidence of Contemporary Dispersal of the Intermediate Snail Host of Schistosoma japonicum: Movement of an NTD Host Is Facilitated by Land Use and Landscape Connectivity

BACKGROUND

While the dispersal of hosts and vectors-through active or passive movement-is known to facilitate the spread and re-emergence of certain infectious diseases, little is known about the movement ecology of Oncomelania spp., intermediate snail host of the parasite Schistosoma japonicum, and its consequences for the spread of schistosomiasis in East and Southeast Asia. In China, despite intense control programs aimed at preventing schistosomiasis transmission, there is evidence in recent years of re-emergence and persistence of infection in some areas, as well as an increase in the spatial extent of the snail host. A quantitative understanding of the dispersal characteristics of the intermediate host can provide new insights into the spatial dynamics of transmission, and can assist public health officials in limiting the geographic spread of infection.

METHODOLOGY/PRINCIPAL FINDINGS

Oncomelania hupensis robertsoni snails (n = 833) were sampled from 29 sites in Sichuan, China, genotyped, and analyzed using Bayesian assignment to estimate the rate of recent snail migration across sites. Landscape connectivity between each site pair was estimated using the geographic distance distributions derived from nine environmental models: Euclidean, topography, incline, wetness, land use, watershed, stream use, streams and channels, and stream velocity. Among sites, 14.4% to 32.8% of sampled snails were identified as recent migrants, with 20 sites comprising >20% migrants. Migration rates were generally low between sites, but at 8 sites, over 10% of the overall host population originated from one proximal site. Greater landscape connectivity was significantly associated with increased odds of migration, with the minimum path distance (as opposed to median or first quartile) emerging as the strongest predictor across all environmental models. Models accounting for land use explained the largest proportion of the variance in migration rates between sites. A greater number of irrigation channels leading into a site was associated with an increase in the site's propensity to both attract and retain snails.

CONCLUSIONS/SIGNIFICANCE

Our findings have important implications for controlling the geographic spread of schistosomiasis in China, through improved understanding of the dispersal capacity of the parasite's intermediate host.

A Complex System of Glacial Sub-Refugia Drives Endemic Freshwater Biodiversity on the Tibetan Plateau

Although only relatively few freshwater invertebrate families are reported from the Tibetan Plateau, the degree of endemism may be high. Many endemic lineages occur within permafrost areas, raising questions about the existence of isolated intra-plateau glacial refugia. Moreover, if such refugia existed, it might be instructive to learn whether they were associated with lakes or with more dynamic ecosystems such as ponds, wetlands, or springs. To study these hypotheses, we used pulmonate snails of the plateau-wide distributed genus Radix as model group and the Lake Donggi Cona drainage system, located in the north-eastern part of the plateau, as model site. First, we performed plateau-wide phylogenetic analyses using mtDNA data to assess the overall relationships of Radix populations inhabiting the Lake Donggi Cona system for revealing refugial lineages. We then conducted regional phylogeographical analyses applying a combination of mtDNA and nuclear AFLP markers to infer the local structure and demographic history of the most abundant endemic Radix clade for identifying location and type of (sub-)refugia within the drainage system. Our phylogenetic analysis showed a high diversity of Radix lineages in the Lake Donggi Cona system. Subsequent phylogeographical analyses of the most abundant endemic clade indicated a habitat-related clustering of genotypes and several Late Pleistocene spatial/demographic expansion events. The most parsimonious explanation for these patterns would be a scenario of an intra-plateau glacial refugium in the Lake Donggi Cona drainage system, which might have consisted of isolated sub-refugia. Though the underlying processes remain unknown, an initial separation of lake and watershed populations could have been triggered by lake-level fluctuations before and during the Last Glacial Maximum. This study inferred the first intra-plateau refugium for freshwater animals on the Tibetan Plateau. It thus sheds new light on the evolutionary history of its endemic taxa and provides important insights into the complex refugial history of a high-altitude ecosystem.

Wetlands explain most in the genetic divergence pattern of Oncomelania hupensis

Understanding the divergence patterns of hosts could shed lights on the prediction of their parasite transmission. No effort has been devoted to understand the drivers of genetic divergence pattern of Oncomelania hupensis, the only intermediate host of Schistosoma japonicum. Based on a compilation of two O. hupensis gene datasets covering a wide geographic range in China and an array of geographical distance and environmental dissimilarity metrics built from earth observation data and ecological niche modeling, we conducted causal modeling analysis via simple, partial Mantel test and local polynomial fitting to understand the interactions among isolation-by-distance, isolation-by-environment, and genetic divergence. We found that geography contributes more to genetic divergence than environmental isolation, and among all variables involved, wetland showed the strongest correlation with the genetic pairwise distances. These results suggested that in China, O. hupensis dispersal is strongly linked to the distribution of wetlands, and the current divergence pattern of both O. hupensis and schistosomiasis might be altered due to the changed wetland pattern with the accomplishment of the Three Gorges Dam and the South-to-North water transfer project.

Molecular characterization of sympatrically distributed Neotricula aperta-like snails in the Mekong River, Kratie, Cambodia

Fifty-six samples of Neotricula aperta-like snails were collected from six locations in Cambodia. Their mitochondrial cytochrome c oxidase subunit 1 (cox1) sequences were examined using haplotype network and neighbor-joining (NJ) tree analysis. Twenty-seven haplotypes (H1-H27) were observed and were divided into two different groups/lineages. Of 27, 17 haplotypes (H11-H27) were clustered with the reference samples of the γ-race N. aperta. The remaining 10 haplotypes (H1-H10) were clustered in a separate group/lineage, differing from the reference samples of the α-, β-, and γ-race N. aperta, suggesting a new lineage belonging the genus Neotricula. Our results show that both the γ-race and a new lineage were sympatrically present approximately 60 km upstream of the Mekong River near the Kratie port, Cambodia. Further morphological and molecular studies are required to confirm the taxonomic status of this new, unidentified lineage.

Genetic structure inferred from mitochondrial 12S ribosomal RNA sequence of Oncomelania quadrasi, the intermediate snail host of Schistosoma japonicum in the Philippines

Species and subspecies of the Oncomelania hupensis species complex are recognized as intermediate hosts of Schistosoma japonicum. Of these species and subspecies, O. quadrasi is distributed throughout the Philippines. This study used 12S ribosomal RNA sequences to explore the genetic structure of O. quadrasi populations in the Philippines. Three subspecies, O. h. hupensis, O. h. formosana, and O. h. chiui of this group were also examined. The phylogenetic tree and haplotypes network showed that O. quadrasi separated from the subspecies. Ten O. quadrasi haplotypes (Oq1-Oq10) clustered in relation to their geographic origin. Genetic differentiation (FST) and estimated gene flow (Nm) among populations showed significant differences, ranging from 0.556-1.000 to 0.00-0.74, respectively. Genetic differences among groups (FCT = 0.466), populations within a group (FSC = 0.727), and populations (FST = 0.854) were observed. These results indicate that the O. quadrasi populations in the Philippines have a substructure associated with their geographic origin.

Spatially explicit modeling of schistosomiasis risk in eastern China based on a synthesis of epidemiological, environmental and intermediate host genetic data

Schistosomiasis japonica is a major parasitic disease threatening millions of people in China. Though overall prevalence was greatly reduced during the second half of the past century, continued persistence in some areas and cases of re-emergence in others remain major concerns. As many regions in China are approaching disease elimination, obtaining quantitative data on Schistosoma japonicum parasites is increasingly difficult. This study examines the distribution of schistosomiasis in eastern China, taking advantage of the fact that the single intermediate host serves as a major transmission bottleneck. Epidemiological, population-genetic and high-resolution ecological data are combined to construct a predictive model capable of estimating the probability that schistosomiasis occurs in a target area (“spatially explicit schistosomiasis risk”). Results show that intermediate host genetic parameters are correlated with the distribution of endemic disease areas, and that five explanatory variables—altitude, minimum temperature, annual precipitation, genetic distance, and haplotype diversity—discriminate between endemic and non-endemic zones. Model predictions are correlated with human infection rates observed at the county level. Visualization of the model indicates that the highest risks of disease occur in the Dongting and Poyang lake regions, as expected, as well as in some floodplain areas of the Yangtze River. High risk areas are interconnected, suggesting the complex hydrological interplay of Dongting and Poyang lakes with the Yangtze River may be important for maintaining schistosomiasis in eastern China. Results demonstrate the value of genetic parameters for risk modeling, and particularly for reducing model prediction error. The findings have important consequences both for understanding the determinants of the current distribution of S. japonicum infections, and for designing future schistosomiasis surveillance and control strategies. The results also highlight how genetic information on taxa that constitute bottlenecks to disease transmission can be of value for risk modeling.

Schistosomiasis is considered the second most devastating parasitic disease after malaria. In China, it is transmitted to humans, cattle and other vertebrate hosts by a single intermediate snail host. It has long been suggested that the close co-evolutionary relationship between parasite and intermediate host makes the snail a major transmission bottleneck in the disease life cycle. Here, we use a novel approach to model the disease distribution in eastern China based on a combination of epidemiological, ecological, and genetic information. We found four major high risk areas for schistosomiasis occurrence in the large lakes and flood plain regions of the Yangtze River. These regions are interconnected, suggesting that the disease may be maintained in eastern China in part through the annual flooding of the Yangtze River, which drives snail transport and admixture of genotypes. The novel approach undertaken yielded improved prediction of schistosomiasis disease distribution in eastern China. Thus, it may also be of value for the predictive modeling of other host- or vector-borne diseases.

Impact of gastrointestinal parasitic nematodes of sheep, and the role of advanced molecular tools for exploring epidemiology and drug resistance - an Australian perspective

Parasitic nematodes (roundworms) of small ruminants and other livestock have major economic impacts worldwide. Despite the impact of the diseases caused by these nematodes and the discovery of new therapeutic agents (anthelmintics), there has been relatively limited progress in the development of practical molecular tools to study the epidemiology of these nematodes. Specific diagnosis underpins parasite control, and the detection and monitoring of anthelmintic resistance in livestock parasites, presently a major concern around the world. The purpose of the present article is to provide a concise account of the biology and knowledge of the epidemiology of the gastrointestinal nematodes (order Strongylida), from an Australian perspective, and to emphasize the importance of utilizing advanced molecular tools for the specific diagnosis of nematode infections for refined investigations of parasite epidemiology and drug resistance detection in combination with conventional methods. It also gives a perspective on the possibility of harnessing genetic, genomic and bioinformatic technologies to better understand parasites and control parasitic diseases.

Assessment of the genetic relationship between Dictyocaulus species from Bos taurus and Cervus elaphus using complete mitochondrial genomic datasets

BACKGROUND

Dictyocaulus species are strongylid nematodes of major veterinary significance in ruminants, such as cattle and cervids, and cause serious bronchitis or pneumonia (dictyocaulosis or "husk"). There has been ongoing controversy surrounding the validity of some Dictyocaulus species and their host specificity. Here, we sequenced and characterized the mitochondrial (mt) genomes of Dictyocaulus viviparus (from Bos taurus) with Dictyocaulus sp. cf. eckerti from red deer (Cervus elaphus), used mt datasets to assess the genetic relationship between these and related parasites, and predicted markers for future population genetic or molecular epidemiological studies.

METHODS

The mt genomes were amplified from single adult males of D. viviparus and Dictyocaulus sp. cf. eckerti (from red deer) by long-PCR, sequenced using 454-technology and annotated using bioinformatic tools. Amino acid sequences inferred from individual genes of each of the two mt genomes were compared, concatenated and subjected to phylogenetic analysis using Bayesian inference (BI), also employing data for other strongylids for comparative purposes.

RESULTS

The circular mt genomes were 13,310 bp (D. viviparus) and 13,296 bp (Dictyocaulus sp. cf. eckerti) in size, and each contained 12 protein-encoding, 22 transfer RNA and 2 ribosomal RNA genes, consistent with other strongylid nematodes sequenced to date. Sliding window analysis identified genes with high or low levels of nucleotide diversity between the mt genomes. At the predicted mt proteomic level, there was an overall sequence difference of 34.5% between D. viviparus and Dictyocaulus sp. cf. eckerti, and amino acid sequence variation within each species was usually much lower than differences between species. Phylogenetic analysis of the concatenated amino acid sequence data for all 12 mt proteins showed that both D. viviparus and Dictyocaulus sp. cf. eckerti were closely related, and grouped to the exclusion of selected members of the superfamilies Metastrongyloidea, Trichostrongyloidea, Ancylostomatoidea and Strongyloidea.

CONCLUSIONS

Consistent with previous findings for nuclear ribosomal DNA sequence data, the present analyses indicate that Dictyocaulus sp. cf. eckerti (red deer) and D. viviparus are separate species. Barcodes in the two mt genomes and proteomes should serve as markers for future studies of the population genetics and/or epidemiology of these and related species of Dictyocaulus.

Molecular characterisation of Galba truncatula, Lymnaea neotropica and L. schirazensis from Cajamarca, Peru and their potential role in transmission of human and animal fascioliasis

Background

Human and animal fascioliasis is emerging in many world regions, among which Andean countries constitute the largest regional hot spot and Peru the country presenting more human endemic areas. A survey was undertaken on the lymnaeid snails inhabiting the hyperendemic area of Cajamarca, where human prevalences are the highest known among the areas presenting a "valley transmission pattern", to establish which species are present, genetically characterise their populations by comparison with other human endemic areas, and discuss which ones have transmission capacity and their potential implications with human and animal infection.

Methods

Therefore, ribosomal DNA ITS-2 and ITS-1, and mitochondrial DNA 16S and cox1 were sequenced by the dideoxy chain-termination method.

Results

Results indicate the presence of three, morphologically similar, small lymnaeid species belonging to the Galba/Fossaria group: Galba truncatula, Lymnaea neotropica and L. schirazensis. Only one combined haplotype for each species was found. The ITS-1, 16S and cox1 haplotypes of G. truncatula are new. No new haplotypes were found in the other two species. This scenario changes previous knowledge, in which only L. viator (= L. viatrix) was mentioned. Galba truncatula appears to be the most abundant, with high population densities and evident anthropophyly including usual presence in human neighbourhood. Infection by Fasciola hepatica larval stages were molecularly confirmed in two populations of this species. The nearness between G. truncatula populations presenting liver fluke infection and both human settings and schools for children, together with the absence of populations of other lymnaeid species in the locality, suggest a direct relationship with human infection.

Conclusions

The geographical overlap of three lymnaeid species poses problems for epidemiological studies and control action. First, a problem in classifying lymnaeid specimens in both field and laboratory activities, given their transmission capacity differences: G. truncatula mainly involved in transmission to humans, L neotropica typically responsible for livestock infection, and L. schirazensis unable for transmission. Although several phenotypic characteristics may be helpful for a preliminary specimen classification, a definitive classification can only be obtained by marker sequencing. Aditionally, L. schirazensis increases the confusion, owing to its ability to mix with other Galba/Fossaria species and distort fascioliasis data such as transmission capacity and infection susceptibility. Second, a problem for epidemiological analysis, surveillance and control by methods as mathematical modelling and Remote Sensing - Geographical Information Systems. In Cajamarca, low resolution mapping may be insufficient, as already verified in Andean areas where different lymnaeid species overlap.