Distinct genetic diversity of Oncomelania hupensis, intermediate host of Schistosoma japonicum in mainland China as revealed by ITS sequences

Characterizing genetic variation on the Z chromosome in Schistosoma japonicum reveals host-parasite co-evolution

BACKGROUND

Schistosomiasis is a neglected tropical disease that afflicts millions of people worldwide; it is caused by Schistosoma, the only dioecious flukes with ZW systems. Schistosoma japonicum is endemic to Asia; the Z chromosome of S. japonicum comprises one-quarter of the entire genome. Detection of positive selection using resequencing data to understand adaptive evolution has been applied to a variety of pathogens, including S. japonicum. However, the contribution of the Z chromosome to evolution and adaptation is often neglected.

METHODS

We obtained 1,077,526 high-quality SNPs on the Z chromosome in 72 S. japonicum using re-sequencing data publicly. To examine the faster Z effect, we compared the sequence divergence of S. japonicum with two closely related species, Schistosoma haematobium and S. mansoni. Genetic diversity was compared between the Z chromosome and autosomes in S. japonicum by calculating the nucleotide diversity (π) and Dxy values. Population structure was also assessed based on PCA and structure analysis. Besides, we employed multiple methods including Tajima's D, FST, iHS, XP-EHH, and CMS to detect positive selection signals on the Z chromosome. Further RNAi knockdown experiments were performed to investigate the potential biological functions of the candidate genes.

RESULTS

Our study found that the Z chromosome of S. japonicum showed faster evolution and more pronounced genetic divergence than autosomes, although the effect may be smaller than the variation among genes. Compared with autosomes, the Z chromosome in S. japonicum had a more pronounced genetic divergence of sub-populations. Notably, we identified a set of candidate genes associated with host-parasite co-evolution. In particular, LCAT exhibited significant selection signals within the Taiwan population. Further RNA interference experiments suggested that LCAT is necessary for S. japonicum survival and propagation in the definitive host. In addition, we identified several genes related to the specificity of the intermediate host in the C-M population, including Rab6 and VCP, which are involved in adaptive immune evasion to the host.

CONCLUSIONS

Our study provides valuable insights into the adaptive evolution of the Z chromosome in S. japonicum and further advances our understanding of the co-evolution of this medically important parasite and its hosts.

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.

Comparative proteomics analysis of Schistosoma japonicum developed in different Oncomelania snails as intermediate hosts

Schistosomiasis is a tropical parasitic disease that seriously endangers humans and animals. In this study, two Oncomelania snails, Oncomelania hupensis (O. hupensis) and Oncomelania weishan (O. weishan), were infected with Schistosoma japonicum (S. japonicum) cercariae during the early period, and ICR mice were subsequently infected with two kinds of miracidia that developed in male and female adult worms. In this study, isobaric tags for relative and absolute quantification (iTRAQ) were used to identify four channels: 113, 115, 117, and 119. A total of 2364 adult schistosome proteins were identified, and 1901 proteins were quantitative. Our results revealed 68 differentially expressed proteins (DEPs) in female adult worms, including 24 upregulated proteins and 44 downregulated proteins, and 55 DEPs in male adult worms, including 25 upregulated proteins and 30 downregulated proteins. LC-MS/MS and bioinformatics analysis indicated that these DEPs are mainly concentrated in cellular composition, molecular function, biological function and catabolism pathways. In summary, this proteomics analysis of adult schistosomes that hatched in two intermediate hosts helps to improve our understanding of the growth and developmental mechanisms of S. japonicum.

A chromosome-level genome of the human blood fluke Schistosoma japonicum identifies the genomic basis of host-switching

The evolution and adaptation of S. japonicum, a zoonotic parasite that causes human schistosomiasis, remain unclear because of the lack of whole-genome data. We construct a chromosome-level S. japonicum genome and analyze it together with 72 samples representing six populations of the entire endemic region. We observe a Taiwan zoophilic lineage splitting from zoonotic populations ∼45,000 years ago, consistent with the divergent history of their intermediate hosts. Interestingly, we detect a severe population bottleneck in S. japonicum, largely coinciding with human history in Asia during the last glacial maximum. We identify several genomic regions underlying natural selection, including GATAD2A and Lmln, both showing remarkable differentiation among different areas. RNAi knockdown suggests association of GATAD2A with parasite development and infection in definitive hosts, while Lmln relates to the specificity of the intermediate hosts. Our study provides insights into the evolution of S. japonicum and serves as a resource for further studies.

Spatial variation of life-history traits in Bulinus truncatus, the intermediate host of schistosomes, in the context of field application of niclosamide in Côte d’Ivoire

Background

Control of intermediate host snails using molluscicides for the control and/or elimination of schistosomiasis is a strategy in line with WHO recommendations. Niclosamide is the main chemical molluscicide recommended by WHO. However, except the immediate killing of the snail, the extent of the impact of the molluscicide application on the evolution of snail life-history traits, in relation to recolonization of treated sites is not well known. This study aimed to characterize the spatial variation of life-history traits of populations of the freshwater snail Bulinus truncatus, in relation to niclosamide spraying in the field.

From 2016 to 2018, we conducted a trial, using niclosamide to control the intermediate host snails for interrupting the seasonal transmission of urinary schistosomiasis in northern and central Côte d’Ivoire. Five villages (sites) were considered, including three test and two control villages. In the test villages, the molluscicide was sprayed in habitats harboring B. truncatus snails three times a year (November, February–March and June). We sampled six B. truncatus populations: two populations from the control villages without any treatment; one collected before treatment and three sampled 2–3 months after treatment of the site with niclosamide. The snail populations were monitored for several life-history traits, including survival, growth, fecundity and hatchability, under laboratory conditions, over one generation (G1). We tested the population, region (North/Centre) and treatment status (treated/untreated) effects on the variation of the measured life-history traits and correlations between pairs of traits were estimated.

Results

On the whole, the traits varied among populations. The risk of death was lower in northern populations compared to central ones. The age at first reproduction was reached earlier with a smaller size of snails in northern populations. Values of first reproduction features (size and fecundity) were lower in treated snail populations. The overall growth of untreated populations was higher than that of treated ones. The late fecundity and egg hatching were higher in northern than in central snails. At first reproduction, age was negatively correlated with some fecundity parameters. However, growth was positively associated with fecundity.

Conclusions

Our study showed a spatial variation of life-history traits in B. truncatus snails. The mollusciciding seems to have led to the depression of some life-history traits in the snail populations. However, investigations should be carried out over several generations of snails to better clarify the impact of niclosamide on the evolution of the life-history traits.

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.

Oncomelania hupensis quadrasi: Snail intermediate host of Schistosoma japonicum in the Philippines

Oncomelania hupensis quadrasi is the snail intermediate host of Schistosoma japonicum in the Philippines. It was discovered by Dr. Marcos Tubangui in 1932 more than two decades after the discovery of the disease in the country in 1906. This review, the first for O. h. quadrasi, presents past and present works on the taxonomy, biology, ecology, control, possible paleogeographic origin of the snail intermediate host and future in research, control and surveillance of the snail. Extensive references are made of other subspecies of O. hupensis such as the subspecies in China for which majority of the advances has been accomplished. Contrasting views on whether the snail is to be considered an independent species of Oncomelania or as one of several subspecies of Oncomelania hupensis are presented. Snail control methods such as chemical methods using synthetic and botanical molluscicides, environmental manipulation and biological control are reviewed. Use of technologies such as Remote Sensing, Geographical Information System and landscape genetics is stressed for snail surveillance. Control and prevention efforts in the Philippines have consistently focused on mass drug administration which has proved inadequate in elimination of the disease. An integrated approach that includes snail control, environmental sanitation and health education has been proposed. Population movement such as migration for employment and economic opportunities and ecotourism and global climate change resulting in heavy rains and flooding challenge the gains of control and elimination efforts. Concern for possible migration of snails to non-endemic areas is expressed given the various changes both natural and mostly man-made favoring habitat expansion.

Prevalence of Schistosoma mansoni and S. haematobium in Snail Intermediate Hosts in Africa: A Systematic Review and Meta-analysis

Background

Schistosomiasis is caused by Schistosoma mansoni and S. haematobium in Africa. These schistosome parasites use freshwater snail intermediate hosts to complete their lifecycle. Varied prevalence rates of these parasites in the snail intermediate hosts were reported from several African countries, but there were no summarized data for policymakers. Therefore, this study was aimed to systematically summarize the prevalence and geographical distribution of S. mansoni and S. haematobium among freshwater snails in Africa.

Methods

Literature search was carried out from PubMed, Science Direct, and Scopus which reported the prevalence of S. mansoni and S. haematobium among freshwater snails in Africa. The pooled prevalence was determined using a random-effect model, while heterogeneities between studies were evaluated by I² test. The meta-analyses were conducted using Stata software, metan command.

Results

A total of 273,643 snails were examined for the presence of S. mansoni and S. haematobium cercaria in the eligible studies. The pooled prevalence of schistosome cercaria among freshwater snails was 5.5% (95% CI: 4.9–6.1%). The pooled prevalence of S. mansoni and S. haematobium cercaria was 5.6% (95% CI: 4.9–6.3%) and 5.2% (95% CI: 4.6–5.7%), respectively. The highest pooled prevalence was observed from Nigeria (19.0%; 95% CI: 12.7–25.3%), while the lowest prevalence was reported from Chad (0.05%; 95% CI: 0.03–0.13). Higher prevalence of schistosome cercaria was observed from Bulinus globosus (12.3%; 95% CI: 6.2–18.3%) followed by Biomphalaria sudanica (6.7%; 95% CI: 4.5–9.0%) and Biomphalaria pfeifferi (5.1%; 95% CI: 4.1–6.2%). The pooled prevalence of schistosome cercaria obtained using PCR was 26.7% in contrast to 4.5% obtained by shedding cercariae.

Conclusion

This study revealed that nearly 6% of freshwater snails in Africa were infected by either S. haematobium or S. mansoni. The high prevalence of schistosomes among freshwater snails highlights the importance of appropriate snail control strategies in Africa.

Oncomelania hupensis retains its ability to transmit Schistosoma japonicum 13 years after migration from permissive to non-permissive areas

Background

The East Route Project (ERP) of the South-to-North Water Diversion Project (SNWDP) stretches across schistosomiasis endemic and non-endemic areas in China, which may lead to the dispersal of Oncomelania hupensis, the intermediate host of Schistosoma japonicum, from permissive areas along the Yangtze River Basin to non-permissive areas in northern China. A previous survey demonstrated that O. hupensis could survive and breed for 13 years (12 generations) after being transferred to a non-permissive area, and could be infected by S. japonicum. However, it is not clear if the migrated snails will change their ability to transmit S. japonicum.

Methods

We infected mice with the cercariae released from the infected transferred snails bred in Jining city of Shandong Province (non-permissive areas) for 13 years. The mice in the control group were infected with cercariae derived from the snails collected in their original habitat (Jiangdu county of Jiangsu Province, permissive areas). Then, we explored the pathogenicity to mice including worm burden, liver egg count and pathology. Additionally, the gene expression profiles of the adult male and female worms recovered from the infected mice were analyzed by RNA sequencing.

Results

The worm burden, liver egg count and pathology of the mice infected with cercariae released from transferred snails bred in non-permissive areas for 13 years showed no significant differences, when compared with the control cercariae. Slight changes occurred at the transcription level between adult male and female worms recovered from mice infected with cercariae derived from snails bred in permissive and non-permissive areas. Only fourteen genes were significantly differentially expressed in the comparison of adult female worms, and no significantly differentially expressed gene was found in the comparison of adult male worms.

Conclusions

Our findings strongly suggest that transferred snails did not change their schistosomiasis transmission ability and the worms derived from them retained the original pathogenicity, even after migrating from permissive to non-permissive areas for 13 years. Therefore, a long-term surveillance system of snails along the SNWDP is urgently needed to prevent the diffusion of O. hupensis and reduce the risk of transmission of schistosomiasis.

Intestinal microbiome profiles in Oncomelania hupensis in mainland China

Oncomelania hupensis plays a significant role in the transmission of schistosomiasis japonica, which remains a major public health concern in China. Understanding the biological characteristics of O. hupensis is a prerequisite for its control; however, there are currently no studies investigating the intestinal microbiota of the O. hupensis snail. This study aimed to profile the intestinal microbiome of O. hupensis across different ecological landscapes in mainland China. DNA was extracted from the intestines of the collected snails and the bacterial communities were detected using 454 pyrosequencing. A total of 3,799 operational taxonomic units (OTUs) were obtained, and Proteobacteria, Firmicutes, and Actinobacteria were identified as the dominant bacterial taxa at the phylum level. Bacillus and Lactococcus were the most common genera in samples obtained from the four ecological landscapes. Snail specimens were clustered into three clades according to microbial community diversity, and thirty-seven genera that contributed to differential microbiota distributions were identified. Co-occurrence network analysis indicated a symbiotic relationship for the intestinal microbiota of O. hupensis, and PICRUSt analysis predicted forty-one metabolic functions in all snail samples, including membrane transport, amino acid metabolism, carbohydrate metabolism, replication and repair, energy metabolism, as well as xenobiotics biodegradation and metabolism. These findings improve our understanding of bacterial ecology in the O. hupensis intestine; further studies will focus on the relationship between O. hupensis intestinal microbiota and the microbiota in their specific ecological environments.

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.

Asian Schistosomiasis: Current Status and Prospects for Control Leading to Elimination

Schistosomiasis is an infectious disease caused by helminth parasites of the genus Schistosoma. Worldwide, an estimated 250 million people are infected with these parasites with the majority of cases occurring in sub-Saharan Africa. Within Asia, three species of Schistosoma cause disease. Schistosoma japonicum is the most prevalent, followed by S. mekongi and S. malayensis. All three species are zoonotic, which causes concern for their control, as successful elimination not only requires management of the human definitive host, but also the animal reservoir hosts. With regard to Asian schistosomiasis, most of the published research has focused on S. japonicum with comparatively little attention paid to S. mekongi and even less focus on S. malayensis. In this review, we examine the three Asian schistosomes and their current status in their endemic countries: Cambodia, Lao People's Democratic Republic, Myanmar, and Thailand (S. mekongi); Malaysia (S. malayensis); and Indonesia, People's Republic of China, and the Philippines (S. japonicum). Prospects for control that could potentially lead to elimination are highlighted as these can inform researchers and disease control managers in other schistosomiasis-endemic areas, particularly in Africa and the Americas.

Molecular Identification and Phylogenetic Analysis of Nuclear rDNA Sequences of Clonorchis sinensis Isolates From Human Fecal Samples in Heilongjiang Province, China

Studying the genetic diversity of parasite is important for understanding their biogeography and molecular epidemiology, as well as for establishing disease prevention and control strategies. Clonorchis sinensis is an important foodborne parasite worldwide. However, despite its epidemiological significance, the genetic diversity of C. sinensis has not been well studied from human in northeastern China. In this study, a total of 342 fecal specimens were collected from residents living in five villages in Heilongjiang Province and analyzed for the presence of C. sinensis by PCR amplification and sequencing of the internal transcribed spacer 1 (ITS1) and ITS2 regions of nuclear ribosomal DNA. 21.64% (74/342) of fecal samples were found to be positive for C. sinensis by PCR. The sequences of the ITS1 region in 34 of the 74 samples (45.95%) matched that of MK179278, Genetic polymorphisms were observed at six nucleotide sites. The ITS2 gene sequence of 37 of the 74 samples (50%) matched that of MK179281. In conclusion, a low degree of genetic diversity between C. sinensis isolates from China and different geographical regions was found at ITS loci. Despite this conservation, sequencing of the rDNA region has provided important data that will be useful for future studies addressing the molecular evolution, biology, medical implications and ecology of C. sinensis.

Transcriptome sequencing and differential gene expression analysis of the schistosome-transmitting snail Oncomelania hupensis inhabiting hilly and marshland regions

The freshwater snail Oncomelania hupensis is the unique intermediate host of the blood fluke Schistosoma japonicum, which is the major cause of schistosomiasis. The snail inhabits two contrasting environments: the hilly and marshland regions. The hilly snails are smaller in size and have the typical smooth shell, whereas the marshland snails are larger and possess the ribbed shell. To reveal the differences in gene expression between the hilly and marshland snails, a total of six snails, three per environment, were individually examined by RNA sequencing technology. All paired-end reads were assembled into contigs from which 34,760 unigenes were predicted. Based on single nucleotide polymorphisms, principal component analysis and neighbor-joining clustering revealed two distinct clusters of hilly and marshland snails. Analysis of expression changes between environments showed that upregulated genes relating to immunity and development were enriched in hilly snails, while those associated with reproduction were over-represented in marshland snails. Eight differentially expressed genes between the two types of snails were validated by qRT-PCR. Our study identified candidate genes that could be targets for future functional studies, and provided a link between expression profiling and ecological adaptation of the snail that may have implications for schistosomiasis control.

Schistosoma japonicum transmission risk maps at present and under climate change in mainland China

BACKGROUND

The South-to-North Water Diversion (SNWD) project is designed to channel fresh water from the Yangtze River north to more industrialized parts of China. An important question is whether future climate change and dispersal via the SNWD may synergistically favor a northward expansion of species involved in hosting and transmitting schistosomiasis in China, specifically the intermediate host, Oncomelania hupensis.

METHODOLOGY/ PRINCIPAL FINDINGS

In this study, climate spaces occupied by the four subspecies of O. hupensis (O. h. hupensis, O. h. robertsoni, O. h. guangxiensis and O. h. tangi) were estimated, and niche conservatism tested among each pair of subspecies. Fine-tuned Maxent (fMaxent) and ensemble models were used to anticipate potential distributions of O. hupensis under future climate change scenarios. We were largely unable to reject the null hypothesis that climatic niches are conserved among the four subspecies, so factors other than climate appear to account for the divergence of O. hupensis populations across mainland China. Both model approaches indicated increased suitability and range expansion in O. h. hupensis in the future; an eastward and northward shift in O. h. robertsioni and O. h. guangxiensis, respectively; and relative distributional stability in O. h. gangi.

CONCLUSIONS/SIGNIFICANCE

The southern parts of the Central Route of SNWD will coincide with suitable areas for O. h. hupensis in 2050-2060; its suitable areas will also expand northward along the southern parts of the Eastern Route by 2080-2090. Our results call for rigorous monitoring and surveillance of schistosomiasis along the southern Central Route and Eastern Route of the SNWD in a future, warmer China.

Assessing environmental factors associated with regional schistosomiasis prevalence in Anhui Province, Peoples' Republic of China using a geographical detector method

Background

Schistosomiasis is a water-borne disease caused by trematode worms belonging to genus Schistosoma, which is prevalent most of the developing world. Transmission of the disease is usually associated with multiple biological characteristics and social factors but also factors can play a role. Few studies have assessed the exact and interactive influence of each factor promoting schistosomiasis transmission.

Methods

We used a series of different detectors (i.e., specific detector, risk detector, ecological detector and interaction detector) to evaluate separate and interactive effects of the environmental factors on schistosomiasis prevalence. Specifically, (i) specific detector quantifies the impact of a risk factor on an observed spatial disease pattern, which were ranked statistically by a value of Power of Determinate (PD) calculation; (ii) risk detector detects high risk areas of a disease on the condition that the study area is stratified by a potential risk factor; (iii) ecological detector explores whether a risk factor is more significant than another in controlling the spatial pattern of a disease; (iv) interaction detector probes whether two risk factors when taken together weaken or enhance one another, or whether they are independent in developing a disease. Infection data of schistosomiasis based on conventional surveys were obtained at the county level from the health authorities in Anhui Province, China and used in combination with information from Chinese weather stations and internationally available environmental data.

Results

The specific detector identified various factors of potential importance as follows: Proximity to Yangtze River (0.322) > Land cover (0.285) > sunshine hours (0.256) > population density (0.109) > altitude (0.090) > the normalized different vegetation index (NDVI) (0.077) > land surface temperature at daytime (LST_day) (0.007). The risk detector indicated that areas of schistosomiasis high risk were located within a buffer distance of 50 km from Yangtze River. The ecological detector disclosed that the factors investigated have significantly different effects. The interaction detector revealed that interaction between the factors enhanced their main effects in most cases.

Conclusion

Proximity to Yangtze River had the strongest effect on schistosomiasis prevalence followed by land cover and sunshine hours, while the remaining factors had only weak influence. Interaction between factors played an even more important role in influencing schistosomiasis prevalence than each factor on its own. High risk regions influenced by strong interactions need to be targeted for disease control intervention.

Electronic supplementary material

The online version of this article (doi:10.1186/s40249-017-0299-x) contains supplementary material, which is available to authorized users.

Genetic variation between Schistosoma japonicum lineages from lake and mountainous regions in China revealed by resequencing whole genomes

Schistosoma infection is a major cause of morbidity and mortality worldwide. Schistosomiasis japonica is endemic in mainland China along the Yangtze River, typically distributed in two geographical categories of lake and mountainous regions. Study on schistosome genetic diversity is of interest in respect of understanding parasite biology and transmission, and formulating control strategy. Certain genetic variations may be associated with adaptations to different ecological habitats. The aim of this study is to gain insight into Schistosoma japonicum genetic variation, evolutionary origin and associated causes of different geographic lineages through examining homozygous Single Nucleotide Polymorphisms (SNPs) based on resequenced genome data. We collected S. japonicum samples from four sites, three in the lake regions (LR) of mid-east (Guichi and Tonglin in Anhui province, Laogang in Hunan province) and one in mountainous region (MR) (Xichang in Sichuan province) of south-west of China, resequenced their genomes using Next Generation Sequencing (NGS) technology, and made use of the available database of S. japonicum draft genomic sequence as a reference in genome mapping. A total of 14,575 SNPs from 2059 genes were identified in the four lineages. Phylogenetic analysis confirmed significant genetic variation exhibited between the different geographical lineages, and further revealed that the MR Xichang lineage is phylogenetically closer to LR Guich lineage than to other two LR lineages, and the MR lineage might be evolved from LR lineages. More than two thirds of detected SNPs were nonsynonymous; functional annotation of the SNP-containing genes showed that they are involved mainly in biological processes such as signaling and response to stimuli. Notably, unique nonsynonymous SNP variations were detected in 66 genes of MR lineage, inferring possible genetic adaption to mountainous ecological condition.

The genetic diversity and geographical separation study of Oncomelania hupensis populations in mainland China using microsatellite loci

Background

Oncomelania hupensis is the unique intermediate host of Schistosoma japonicum, which plays a crucial role in the transmission of schistosomiasis. The endemic area of S. japonicum is strictly consistent with the geographical distribution of O. hupensis.

Methods

A total of 24 populations of O. hupensis from four ecological landscapes were selected for analysis of genetic diversity by screening eight microsatellite DNA polymorphic loci.

Results

The number of alleles per locus ranged from 29 to 70 with an average of 45.625 and that of effective alleles were 18.5 to 45.8 with an average of 27.4. The observed (Ho) and expected (He) heterozygosities varied from 0.331 to 0.57 and from 0.888 to 0.974, respectively. The mean of polymorphism information content (PIC) for all populations was 0.940, appearing polymorphic for all loci. For the fixation index of F-Statistics, Fit and Fst were 54.95 and 37.62 %, respectively. Variation of O. hupensis chiefly exists among individuals, accounting for 60.58 % of the total variation determined by Analysis of Molecular Variation (AMOVA). Variation among individuals within populations, among populations within groups and among groups only accounted for 26.60, 8.04 and 4.78 %, respectively. This distribution of variation suggests that genetic differences principally originate from within-populations rather than among-populations. Moreover, UPGMA cluster analysis showed that the populations spreading within middle and lower reaches of the Yangtze River (HBWH, JSYZ, JXNC, HNHS, JXJJ, AHWW, HBJL, JXDC, HNNX, JSYZJZ, ZJJH, AHNG and AHWJ) clustered together first, then gathered with the populations in the high mountains (SCMS, SCYA, SCPJ, YNEY, SCLS, YNWS and SCXC), coastal hills (FJFQ and FJFZ) and Karst landform (GXBS and GXYZ) successively.

Conclusion

This study provides novel insight into the theoretical source of genetic differentiation of Oncomelania hupensis in mainland China, which is critical for the epidemiological investigation and surveillance of S. japonicum.

Co-dispersal of the blood fluke Schistosoma japonicum and Homo sapiens in the Neolithic Age

The global spread of human infectious diseases is of considerable public health and biomedical interest. Little is known about the relationship between the distribution of ancient parasites and that of their human hosts. Schistosoma japonicum is one of the three major species of schistosome blood flukes causing the disease of schistosomiasis in humans. The parasite is prevalent in East and Southeast Asia, including the People's Republic of China, the Philippines and Indonesia. We studied the co-expansion of S. japonicum and its human definitive host. Phylogenetic reconstruction based on complete mitochondrial genome sequences showed that S. japonicum radiated from the middle and lower reaches of the Yangtze River to the mountainous areas of China, Japan and Southeast Asia. In addition, the parasite experienced two population expansions during the Neolithic agriculture era, coinciding with human migration and population growth. The data indicate that the advent of rice planting likely played a key role in the spread of schistosomiasis in Asia. Moreover, the presence of different subspecies of Oncomelania hupensis intermediate host snails in different localities in Asia allowed S. japonicum to survive in new rice-planting areas, and concurrently drove the intraspecies divergence of the parasite.

Genetic Diversity of Schistosoma haematobium Eggs Isolated from Human Urine in Sudan

The genetic diversity of Schistosoma haematobium remains largely unstudied in comparison to that of Schistosoma mansoni. To characterize the extent of genetic diversity in S. haematobium among its definitive host (humans), we collected S. haematobium eggs from the urine of 73 infected schoolchildren at 5 primary schools in White Nile State, Sudan, and then performed a randomly amplified polymorphic DNA marker ITS2 by PCR-RFLP analysis. Among 73 S. haematobium egg-positive cases, 13 were selected based on the presence of the S. haematobium satellite markers A4 and B2 in their genomic DNA, and used for RFLP analysis. The 13 samples were subjected to an RFLP analysis of the S. haematobium ITS2 region; however, there was no variation in size among the fragments. Compared to the ITS2 sequences obtained for S. haematobium from Kenya, the nucleotide sequences of the ITS2 regions of S. haematobium from 4 areas in Sudan were consistent with those from Kenya (> 99%). In this study, we demonstrate for the first time that most of the S. haematobium population in Sudan consists of a pan-African S. haematobium genotype; however, we also report the discovery of Kenyan strain inflow into White Nile, Sudan.

De Novo transcriptome analysis of Oncomelania hupensis after molluscicide treatment by next-generation sequencing: implications for biology and future snail interventions

The freshwater snail Oncomelania hupensis is the only intermediate host of Schistosoma japonicum, which causes schistosomiasis. This disease is endemic in the Far East, especially in mainland China. Because niclosamide is the only molluscicide recommended by the World Health Organization, 50% wettable powder of niclosamide ethanolamine salt (WPN), the only chemical molluscicide available in China, has been widely used as the main snail control method for over two decades. Recently, a novel molluscicide derived from niclosamide, the salt of quinoid-2',5-dichloro-4'-nitro-salicylanilide (Liu Dai Shui Yang An, LDS), has been developed and proven to have the same molluscicidal effect as WPN, with lower cost and significantly lower toxicity to fish than WPN. The mechanism by which these molluscicides cause snail death is not known. Here, we report the next-generation transcriptome sequencing of O. hupensis; 145,008,667 clean reads were generated and assembled into 254,286 unigenes. Using GO and KEGG databases, 14,860 unigenes were assigned GO annotations and 4,686 unigenes were mapped to 250 KEGG pathways. Many sequences involved in key processes associated with biological regulation and innate immunity have been identified. After the snails were exposed to LDS and WPN, 254 unigenes showed significant differential expression. These genes were shown to be involved in cell structure defects and the inhibition of neurohumoral transmission and energy metabolism, which may cause snail death. Gene expression patterns differed after exposure to LDS and WPN, and these differences must be elucidated by the identification and annotation of these unknown unigenes. We believe that this first large-scale transcriptome dataset for O. hupensis will provide an opportunity for the in-depth analysis of this biomedically important freshwater snail at the molecular level and accelerate studies of the O. hupensis genome. The data elucidating the molluscicidal mechanism will be of great benefit in future snail control efforts.

Geographical genetic structure of Schistosoma japonicum revealed by analysis of mitochondrial DNA and microsatellite markers

BACKGROUND

Schistosoma japonicum is a significant public health risk in parts of China and elsewhere in Southeast Asia. To gain an insight into the epidemiology of schistosomiasis japonica, a detailed investigation of S. japonicum genetic population structure is needed.

METHODS

Using three mitochondrial DNA fragments and ten microsatellite loci, we investigated the genetic diversity within and structure among twelve populations of S. japonicum sampled on a geographical scale covering most major endemic areas.

RESULTS

Schistosoma japonicum lineages from Indonesia, the Philippines and Chinese Taiwan were clearly distinct from each other and from those in mainland China. Within mainland China, there was some evidence for genetic divergence between populations from the mountain and lake regions. However, the analysis inferred no clear sub-population structure in the lake region of mainland China. High genetic diversity was found among S. japonicum populations of mainland China and this was significantly higher than those from island regions.

CONCLUSIONS

High genetic diversity within and substantial differentiation among populations were demonstrated in S. japonicum.

Resistance to niclosamide in Oncomelania hupensis, the intermediate host of Schistosoma japonicum: should we be worried?

As the currently only available molluscicide, niclosamide has been widely used for snail control for over 2 decades in China. There is therefore a concern about the emergence of niclosamide-resistant snail populations following repeated, extensive use of the chemical. The purpose of this study was to investigate the likelihood of niclosamide resistance in Oncomelania hupensis in China. Active adult O. hupensis snails derived from 20 counties of 10 schistosomiasis-endemic provinces of China, of 10 snails in each drug concentration, were immersed in solutions of 1, 0.5, 0.25, 0.125, 0.063, 0.032, 0.016 and 0.008 mg L-1 of a 50% wettable powder of niclosamide ethanolamine salt (WPN) for 24 and 48 h at 25 °C, and the median lethal concentration (LC50) was estimated. Then, the 24- and 48-h WPN LC50 values were compared with those determined in the same sampling sites in 2002. The results indicated that the 24- and 48-h WPN LC50 values for O. hupensis were not significantly different from those determined in 2002 (P = 0.202 and 0.796, respectively). It is concluded that the current sensitivity of O. hupensis to niclosamide has not changed after more than 2 decades of repeated, extensive application in the main endemic foci of China, and there is no evidence of resistance to niclosamide detected in O. hupensis.

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.

Spatial pattern of schistosomiasis in Xingzi, Jiangxi Province, China: the effects of environmental factors

Background

The recent rebounds of schistosomiasis in the middle and lower reaches of the Yangtze River pose a challenge to the current control strategies. In this study, identification of potential high risk snail habitats was proposed, as an alternative sustainable control strategy, in Xingzi County, China. Parasitological data from standardized surveys were available for 36,208 locals (aged between 6–65 years) from 42 sample villages across the county and used in combination with environmental data to investigate the spatial pattern of schistosomiasis risks.

Methods

Environmental factors measured at village level were examined as possible risk factors by fitting a logistic regression model to schsitosomiasis risk. The approach of ordinary kriging was then used to predict the prevalence of schistosomiasis over the whole county.

Results

Risk analysis indicated that distance to snail habitat and wetland, rainfall, land surface temperature, hours of daylight, and vegetation are significantly associated with infection and the residual spatial pattern of infection showed no spatial correlation. The predictive map illustrated that high risk regions were located close to Beng Lake, Liaohuachi Lake, and Shixia Lake.

Conclusions

Those significant environmental factors can perfectly explain spatial variation in infection and the high risk snail habitats delineated by the predicted map of schistosomiasis risks will help local decision-makers to develop a more sustainable control strategy.

Water-related parasitic diseases in China

Water-related parasitic diseases are directly dependent on water bodies for their spread or as a habitat for indispensable intermediate or final hosts. Along with socioeconomic development and improvement of sanitation, overall prevalence is declining in the China. However, the heterogeneity in economic development and the inequity of access to public services result in considerable burden due to parasitic diseases in certain areas and populations across the country. In this review, we demonstrated three aspects of ten major water-related parasitic diseases, i.e., the biology and pathogenicity, epidemiology and recent advances in research in China. General measures for diseases control and special control strategies are summarized.

A specific PCR assay for the identification and differentiation of Schistosoma japonicum geographical isolates in mainland China based on analysis of mitochondrial genome sequences

In the present study, near-complete mt genome sequences for eight representative Schistosoma japonicum samples from seven endemic provinces in mainland China were analyzed. Sequence differences among the eight mt genomes of S. japonicum samples were 0.20-2.51%. Variation in protein-coding genes was greater than that in rRNA genes. The mt DNA sequences of S. japonicum samples from south-western (SW) China were 2 bp [position 11727-11728 within tRNA-Cys, microsatellite (AG) indel] longer than those of the parasites from the lower Yangtze/Zhejiang areas. Representative DNA sequencing confirmed that such (AG) indel could be exploited for identification and differentiation of S. japonicum populations in SW China's Yunnan and Sichuan province which have two (AG) repeats from those in all remaining endemic provinces along the Yangtze River below the Three Gorges regions or close to the east coast of China (e.g., Zhejiang) which have only one (AG) repeat. Phylogenetic analyses based on the concatenated amino acids of 12 protein-coding genes also showed that samples from SW China (Sichuan and Yunnan provinces), above the Three Gorges Dam, formed a distinct cluster. Based on this indel polymorphism, a pair of specific primers was designed and used to develop a specific-PCR polyacrylamide gel detection assay. There was an obvious length difference in the amplified PCR products between S. japonicum samples from the two endemic types. The specific-PCR assay allowed the specific identification of S. japonicum, with no amplicons being amplified from other closely related trematodes, and the minimum amount of DNA detectable was 0.05 ng. This approach is inexpensive, easy to perform and the whole detection process can be completed within 4h. Examination of 81 S. japonicum samples from SW China's Yunnan and Sichuan provinces, and 264 samples from the lower Yangtze provinces (Hubei, Jiangsu, Jiangxi, Anhui and Hunan) and from Zhejiang validated the value of the specific PCR assay and proved its reliability. These findings indicate that the specific PCR assay would provide a useful tool for the epidemiological surveillance and for tracing the source of S. japonicum infection in humans and animals in China.

Diversification of Schistosoma japonicum in Mainland China revealed by mitochondrial DNA

Background

Schistosoma japonicum still causes severe parasitic disease in mainland China, but mainly in areas along the Yangtze River. However, the genetic diversity in populations of S. japonicum has not been well understood across its geographical distribution, and such data may provide insights into the epidemiology and possible control strategies for schistosomiasis.

Methodology/Principal Findings

In this study infected Oncomelania snails were collected from areas in the middle and lower (ML) reaches of the Yangtze River, including Hubei, Hunan, Anhui, Jiangxi and Jiangsu provinces, and in the upper reaches of the river, including Sichuan and Yunnan provinces in southwest (SW) China. The adult parasites obtained from experimentally infected mice using isolated cercariae were sequenced individually for several fragments of mitochondrial regions, including Cytb-ND4L-ND4, 16S-12S and ND1. Populations in the ML reaches exhibited a relatively high level of diversity in nucleotides and haplotypes, whereas a low level was observed for populations in the SW, using either each single fragment or the combined sequence of the three fragments. Pairwise analyses of F-statistics (Fst) revealed a significant genetic difference between populations in the ML reaches and those in the SW, with limited gene flow and no shared haplotypes in between. It is rather obvious that genetic diversity in the populations of S. japonicum was significantly correlated with the geographical distance, and the geographical separation/isolation was considered to be the major factor accounting for the observed difference between populations in the ML reaches and those in the SW in China.

Conclusions

S. japonicum in mainland China exhibits a high degree of genetic diversity, with a similar pattern of genetic diversity as observed in the intermediate host snails in the same region in China.

Despite the existing threat of schistosomiasis in some rural areas along the Yangtze River, the genetic diversity of Schistosoma japonicum has not been investigated across its wide geographical distribution in China, and such information may provide insight into the disease epidemiology and the development of its control measures. In this study, the adult parasites, obtained through infecting mice with cercariae from snails of the genus Oncomelania collected from a wide range of localities in currently endemic areas of schistosomiasis in the middle and lower (ML) reaches of the Yangtze River, and in Sichuan and Yunnan provinces in the upper reaches of the river in southwest (SW) China, were sequenced individually for mitochondrial genes. In general, a relatively high degree of genetic variation was observed in populations in the ML reaches in terms of nucleotide and haplotype diversity, but a low level was observed in populations in the SW. The significant difference in genetic diversity as revealed by F-statistics, and the existence of no shared haplotypes, were observed between populations in the ML reaches and those in the SW, indicating the effect of geographical separation/isolation upon the schistosomes and probably the parasite-snail system in China.

Controlled Synthesis of Water Soluble Silver Nanoparticles and the Molluscicidal Effect

The silver nanoparticles with 42 nm in diameter were synthesized under control of cetyltrimethyl ammonium bromide (CTAB). CTAB acted not only as reactant but also as dispersing agent. The efficient molluscicidal effect of the silver nanoparticles showed that a new way of the schistosomiasis prevention has been explored.

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

BACKGROUND

Oncomelania hupensis robertsoni is the sole intermediate host for Schistosoma japonicum in western China. Given the close co-evolutionary relationships between snail host and parasite, there is interest in understanding the distribution of distinct snail phylogroups as well as regional population structures. Therefore, this study focuses on these aspects in a re-emergent schistosomiasis area known to harbour representatives of two phylogroups - the Deyang-Mianyang area in Sichuan Province, China. Based on a combination of mitochondrial and nuclear DNA, the following questions were addressed: 1) the phylogeography of the two O. h. robertsoni phylogroups, 2) regional and local population structure in space and time, and 3) patterns of local dispersal under different isolation-by-distance scenarios.

RESULTS

The phylogenetic analyses confirmed the existence of two distinct phylogroups within O. h. robertsoni. In the study area, phylogroups appear to be separated by a mountain range. Local specimens belonging to the respective phylogroups form monophyletic clades, indicating a high degree of lineage endemicity. Molecular clock estimations reveal that local lineages are at least 0.69-1.58 million years (My) old and phylogeographical analyses demonstrate that local, watershed and regional effects contribute to population structure. For example, Analyses of Molecular Variances (AMOVAs) show that medium-scale watersheds are well reflected in population structures and Mantel tests indicate isolation-by-distance effects along waterways.

CONCLUSIONS

The analyses revealed a deep, complex and hierarchical structure in O. h. robertsoni, likely reflecting a long and diverse evolutionary history. The findings have implications for understanding disease transmission. From a co-evolutionary standpoint, the divergence of the two phylogroups raises species level questions in O. h. robertsoni and also argues for future studies relative to the distinctness of the respective parasites. The endemicity of snail lineages at the regional level supports the concept of endemic schistosomiasis areas and calls for future geospatial analyses for a better understanding of respective boundaries. Finally, local snail dispersal mainly occurs along waterways and can be best described by using cost distance, thus potentially enabling a more precise modelling of snail, and therefore, parasite dispersal.

Compensatory density feedback of Oncomelania hupensis populations in two different environmental settings in China

Background

The most recent strategy for schistosomiasis control in the People's Republic of China aims to reduce the likelihood of environmental contamination of schistosome eggs. Despite considerable progress, it is believed that achievements would be further consolidated with additional intermediate host snail control measures. We provide an empirical framework for discerning the relative contribution of intrinsic effects (density feedback) from other extrinsic drivers of snail population dynamics.

Methods

We set up experiments in two study locations to collect reproduction data of Oncomelania hupensis, the intermediate host snail of Schistosoma japonicum. We applied a set of four population dynamic models that have been widely used to study phenomenological time-series data to examine the properties of demographic density feedback patterns from abundance data. We also contrasted the obtained results with the component feedback of density on survival rate to determine whether adult survival was the principal driver of the demographic feedback observed.

Results

Demographic density feedback models (Ricker- and Gompertz-logistic) accounted for > 99% of Akaike's information criterion model weight, with the Gompertz ranking highest in all O. hupensis population groups. We found some evidence for stronger compensatory feedback in the O. hupensis population from Sichuan compared to a Jiangsu population. Survival rates revealed strong component feedback, but the log-linear relationships (i.e. Gompertz) had less support in the demographic feedback analysis.

Conclusions

Our findings indicate that integrated schistosomiasis control measures must continue to reduce parasite abundance further because intermediate host snail populations tend to grow exponentially at low densities, especially O. hupensis populations in mountainous regions. We conclude that density feedback in adult survival is the principal component contribution to the demographic phenomenon observed in the population fitness (r)-abundance relationship.

Conservation and variation in mitochondrial genomes of gastropods Oncomelania hupensis and Tricula hortensis, intermediate host snails of Schistosoma in China

The complete mitochondrial genomes of intermediate host snails for Schistosoma in China were sequenced, including the sub-species Oncomelania hupensis hupensis in two types, and O. hupensis robertsoni, intermediate hosts for S. japonicum, and Tricula hortensis, the intermediate host of S. sinensium. Four genomes have completely the same gene order as in other caenogastropods, containing 13 protein-coding genes and 22 transfer RNAs. The gene size, start codon and termination codon are mostly the same for all protein-coding genes. However, pairwise sequence alignments revealed quite different degrees of variation. The ribbed-shelled O. hupensis hupensis and the smooth-shelled but with varix O. hupensis hupensis had a lower level of genetic distance (3.1% for protein-coding genes), but the coden usages differed obviously in the mitochondrial genomes of these two types of snails, implying that their genetic difference may be larger than previously recognized. The mean genetic distance between O. hupensis hupensis and O. hupensis robertsoni was 12% for protein-coding genes, indicating a higher degree of genetic difference. In consideration of the difference in morphology and distribution, we considered that O. hupensis hupensis and O. hupensis robertsoni can be considered as separate species. The ribbed-shelled O. hupensishupensis and smooth-shelled O. hupensis robertsoni were phylogenetically clustered together within a same clade, which was then clustered with T. hortensis, confirming their close relationship. However, species or sub-species in the Oncomelania from southeastern Asian countries should be included in future study in order to resolve the phylogenetic relationship and origination of all snails in the genus.