- Open Access
Geographical genetic structure of Schistosoma japonicum revealed by analysis of mitochondrial DNA and microsatellite markers
Parasites & Vectorsvolume 8, Article number: 150 (2015)
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.
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.
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.
High genetic diversity within and substantial differentiation among populations were demonstrated in S. japonicum.
Schistosomiasis is one of the most serious of human parasitic diseases, infecting an estimated 200 million people in 76 tropical and subtropical countries . Three major species, Schistosoma japonicum, S. mansoni and S. haematobium, can infect humans . Of these, S. mansoni is endemic in Africa, the Middle East and South America and S. haematobium is found in Africa and the Middle East, whereas S. japonicum occurs in China, Indonesia and the Philippines. At present, this parasite is endemic in China in the lake/marshland regions of Anhui, Hubei, Hunan and Jiangxi provinces and in the mountainous areas of Sichuan and Yunnan provinces . Although tremendous efforts over the last 20–30 years have reduced the prevalence of S. japonicum in China by more than 90% , recent indications are that schistosomiasis has re-emerged in certain areas over the last decade . To gain an insight into the transmission and epidemiology of Asian schistosomiasis, a detailed investigation of S. japonicum population genetic structure is required.
Several classes of molecular markers have been applied to investigate the genetic variability of S. japonicum populations. Earlier studies using isoenzymes  and mitochondrial (mt) DNA sequences  reported remarkable genetic similarity between Chinese and Philippine S. japonicum lineages. Later, differences among, but not within, Chinese S. japonicum populations were detected using mt DNA sequences (nad1) . Similarly, Gasser et al.  applied random amplified polymorphism DNA (RAPD) analysis, and found genetic variability among Chinese S. japonicum populations. Subsequently, a set of high-resolution microsatellite markers was used to investigate the genetic structure of S. japonicum populations from eight locations across China, and a high level of polymorphism was detected both between and within populations . Moreover, S. japonicum populations from the lake and mountainous regions were genetically distinct . Similar results from mt DNA data suggested that geographical separation was the major factor accounting for the population divergence of S. japonicum between the lake and mountainous regions in China [3,11]. Previous studies sampled S. japonicum either from mainland China, or a few locations in mainland China and other Asian countries [12,13], but did not explore the genetic population structure of S. japonicum across its entire range.
Most phylogeographical studies have relied on the analysis of mitochondrial DNA variation, because of its rapid evolution and maternal transmission without intermolecular recombination . In contrast, microsatellites can provide information on relatively recent evolutionary processes . Assessing the congruence among independent genetic markers has become an important issue . Indeed, more and more studies have inferred genetic relationships among populations based on the combined analysis of different classes of markers [16-18]. Consistent results provide a robust perspective on species’ evolutionary history , while findings of incongruent patterns from different markers are also valuable to demonstrate important evolutionary processes, such as introgressive hybridization . To date, no study has been undertaken on S. japonicum to investigate the genetic relationships among populations using different classes of genetic markers simultaneously on the same samples.
In this work, we analyzed S. japonicum samples from eight locations in China mainland (covering the lake and mountainous regions) and Taiwan, Indonesia, Japan and the Philippines. We explored the genetic relationships among these populations based on the analysis of three mt genes (nad1, nad4 and 16 s-12 s) and ten microsatellite loci from the same individual worms.
All procedures involving animals were carried out based on the guidelines of the Association for Assessment and Accreditation of Laboratory Animal Care International. The study protocol followed institutional ethical guidelines that were approved by the ethics committee at the National Institute of Parasitic Diseases, Chinese Center for Disease Control and Prevention (NIPD, China CDC; Permit No: IPD2008-4).
Adult individuals of S. japonicum were obtained from twelve different geographical locations (Figure 1), including eight locations across six provinces in mainland China (Anhui, Hunan, Hubei and Jiangxi Provinces in the Yangtze River Basin/lake region, and Sichuan and Yunnan Provinces in the mountainous area of Southwest China) and four other island locations from Asia (Taiwan, Indonesia, Japan and the Philippines). Infected snails (Oncomelania hupensis) from each mainland Chinese location were transported to the laboratory of NIPD, China CDC, Shanghai. Cercariae were released from pooled infected snails from each site and used to infect laboratory-raised rabbits. Finally, 45 days later, the adult schistosomes were perfused from the mesenteric veins of the infected rabbits and washed in saline before being preserved in 95% (v/v) ethanol at 4°C. The S. japonicum isolate from Leyte, the Philippines, was taken into culture originally in 1969 by Dr. Scholice. The material sent to us, as lyophilized adult worms, was provided by Dr. John Bruce, Centre for Tropical Diseases, University of Lowell, USA. The worms from Japan (Kofu) were gifted by Dr. Hiroshi Yamasaki, National Institute of Infectious Diseases, Tokyo, Japan. The other geographical samples (lyophilized adult worms) from Indonesia (Lindu Lake, Sulawesi) and Taiwan (Changhua) were provided by Dr. John Cross, Uniformed Services University of the Health Sciences, Bethesda, USA. These preserved samples were from populations that had been maintained in laboratory animals.
In total, 42 worms were subjected to mitochondrial DNA analysis and 486 individuals were used for the microsatellite analysis (Table 1). Microsatellite data for 401 of the worms came from our previous work .
Genomic DNA was extracted from individual adult schistosomes using the DNeasy Blood & Tissue Kit and Animal Tissues (Spin-Column) protocol from QIAGEN (Hilden, Germany), and stored at −20°C until use.
Mitochondrial gene amplification and sequencing
We amplified three fragments of the mitochondrial genome: 891 bp of the gene encoding NADH dehydrogenase subunit 1 (nad1), 1275 bp of NADH dehydrogenase subunit 4 (nad4) and 1748 bp of 16S-12S mitochondrial ribosomal DNA (16S-12S rRNA). Three pairs of primers were designed based on the complete reference mitochondrial sequence for S. japonicum (ID in GenBank: AF215860) as follows:
16S-12S rRNA-F (5′-ATAATGTTGCGTCTAAGGTC)
and 16S-12S rRNA-R (5′-TAAACACTACCCATCAAATC).
The amplifications were performed in 20 μl reactions containing 0.5 μM of each primer, ~ 50 ng of genomic DNA from an individual worm, Has Taq polymerase (1.25 U, TaKaRa), 1.25 mM MgCl2, 2 μl 10× reaction buffer and 1 μl dNTPs (2.5 mM, TaKaRa). The PCR protocol used a hot-start activation at 94°C for 10 min, followed by 35 cycles of denaturation for 30 s at 94°C, 30 s annealing (55°C for nad1 and nad4, 49°C for 16S-12S rRNA) and extension (90 s for nad1, 150 s for nad4 and 16S-12S rRNA) at 72°C, followed by a final extension at 72°C for 10 min. PCR products were examined using agarose gel electrophoresis (1% w/v) to validate amplification efficiency. Amplification products were then sent to BGI (Beijing Genome Institute; Shanghai, China) for sequencing using an ABI 3730 DNA Analyzer.
The DNA of each individual S. japonicum worm was genotyped at ten microsatellite loci (i.e. Sjp1, 4, 5, 6, 8, 9, 10, 14, 15, and 17), which were characterized previously . The PCRs were performed based on the protocol described in Yin et al. . The PCR products were diluted to an appropriate concentration and analyzed on an ABI 3730 capillary automated sequencer using a LIZ 500 labeled size standard. The allele sizes were read using GeneMapper software (Version 4.0) and checked manually.
Mitochondrial DNA analysis
Unique haplotypes for the mitochondrial DNA sequences (concatenated nad1, nad4 and 16S-12S rRNA sequences; hereafter “combined mtDNAs”) were identified using Arlequin 3.11 . All unique haplotypes were aligned in ClustalW . A phylogenetic tree was constructed using three methods, namely neighbor-joining (NJ), maximum-likelihood (ML) and maximum parsimony (MP), in MEGA 5.10 . A consensus tree was obtained after bootstrap analysis using 104 bootstrap replicates. Schistosoma mansoni (ID in GenBank: NC_002545.1) was used as an outgroup. A median-joining network was generated to infer the relationships among the haplotypes of S. japonicum using Network . Nucleotide diversity (π) per population was calculated in DnaSP Version 5 : comparisons between regions (mainland China vs island region and lake vs mountainous regions in mainland China) were made using unpaired Student’s t test. Hierarchical analysis of molecular variance (AMOVA) was applied to partition the genetic variances into among regions (i.e. mountain and lake regions), among populations within region and within populations, for samples from mainland China, in Arlequin.
Multiple genotypes from a single snail are often very similar to one another and possibly derived from a single miracidium by mutations occurring during the proliferation of cercariae inside the snail . We therefore retained only a single individual from each cluster of near-identical Multi Locus Genotypes (near-identical MLGs) which formed tight single-sex clusters in principal-coordinates analyses (PCoAs). After this, the PCoA was applied to explore the relationships among remaining MLGs. For each population, the genetic diversity was also examined by calculating the expected heterozygosity (H e ) and observed heterozygosity (H o ) in GenAIEx. In this and following analyses, the individuals for which data were missing were excluded, as were populations with N < 6. Third, the inbreeding coefficient (F IS ), which measures the extent of nonrandom mating, was computed in Genepop 3 . Positive values of F IS indicate heterozygote deficiency, whereas, negative F IS values suggest heterozygote excess. Finally, deviations from Hardy-Weinberg equilibrium (HWE) were tested in Genepop, with 104 permutations.
Genetic population differentiation was estimated using Wright’s F-statistics (F ST ) in GenAIEx, and the significance of the F ST value was tested by 104 permutations. AMOVA was applied to quantify intra- and inter-population and region variances in Arlequin. A Mantel test was used to test whether there was a relationship between pairwise geographical distance (km) and pairwise genetic distance (F ST , based on microsatellites) in GenAIEx.
We used STRUCTURE 2.3.4  to cluster the samples on the basis of their microsatellite genotypes. The most likely number of subpopulations was estimated from the same dataset by running the procedure with a value of K from 1 to 13 using the admixture model . In this process, 106 iterations were performed after a burn-in of length 106.
Mitochondrial DNA analysis
We successfully amplified and sequenced 42 S. japonicum individuals for three complete mitochondrial fragments (nad1, nad4 and 16S-12S rRNA). All the sequences were submitted to GenBank under accession numbers KP793743- KP793784 for nad1, KP793785-KP793826 for nad4, KP793827-KP793868 for 16S rRNA and KP793869-KP793910 for 12S rRNA.
Among the 42 samples sequenced, 34 haplotypes were identified (Table 1). No haplotype was found at more than one geographical location. Phylogenetic trees, constructed by different methods (i.e. NJ, ML, and MP), showed identical or similar topologies with only a few differences in bootstrap values (not shown). The Taiwanese sequences formed a group very distant from the others. Lineages from the mountainous region of China (Sichuan and Yunnan) and the other island locations (Indonesia, the Philippines and Japan) were each distinct, forming monophletic clades (almost so in the case of Yunnan) supported by relatively high bootstrap values (Figure 2). In contrast, there was no obvious structure for samples from the lake region in China, with samples from Anhui, Hunan and Jiangxi intermingled across the tree (Figure 2). This pattern was further supported by the median-joining network (Additional file 1: Figure S1).
The nucleotide diversity (π) ranged from 0 to 0.005 (mean = 0.0016). Interestingly, the nucleotide diversities of most (five of eight) populations from mainland China were significantly higher than those from Chinese Taiwan and each of the other island regions (mean for mainland China 0.0022 vs 0.0005; N = 12, t = 2.52, P = 0.03). However, we detected no significant differences in diversity between populations from the mountainous area and lake region of mainland China (N = 8, t = 1.80, P = 0.13). We recognize that sample sizes are small and this might influence the results. The four individuals from Anhui Guichi shared the same haplotype. In a further seven populations (Anhui Tonglin, Hunan Changde, Hunan Yueyang, Sichuan Xichang, Yunnan Eryuan, Indonesia and Japan), each individual had a unique haplotype (Table 1). The AMOVA revealed that most of the genetic variation was at the within-population level (62.53%), however, significant genetic differentiation was found between the mountain and lake regions (6.25%; P < 0.001).
At the 10 microsatellite loci, 317 unique MLGs were detected from 353 individuals (individuals with missing data were excluded: Table 1), and no MLG was shared among populations. According to the PCoA, the MLGs of 262 worms (after exclusion of near-identical MLGs) could be divided into three groups: individuals from the Philippines, Indonesia and Chinese Taiwan were clustered according to their geographical origins and all the worms from mainland China and Japan were grouped together (Figure 3). Identical results were obtained when the PCoA was run using only the 42 individuals for which mitochondrial sequence data were also obtained (not shown). Observed heterozygosity (H o ) was consistently lower than the expected heterozygosity (H e ), and most of the inbreeding coefficient (F IS ) values were positive, with the exception of Indonesia and Japan (Table 1), indicating an excess of homozygotes within subpopulations in comparison to HWE. Indeed, all populations deviated significantly from HWE (Table 1).
Although most of the genetic variation occurred at the within-population level (74.91%) according to hierarchical AMOVA, significant genetic differentiation was detected among the populations. Pairwise F ST values ranged from 0.01 to 0.51 (averaged over 10 loci; Table 2), suggesting low to high population differentiation. Specifically, the F ST between the Taiwanese and Indonesian populations was 0.51 (Table 2), indicating high population differentiation. In contrast, relatively low differentiation was seen among the populations from mainland China (F ST ranged from 0.01 to 0.05; mean = 0.03; Table 2). The Mantel test for isolation by distance revealed a positive, but non-significant, correlation between population genetic differentiation (F ST ) based on the microsatellite data and geographical distance (R 2 = 0.44; P = 0.09; Figure 4).
Based on the highest value of ΔK, the most likely number (K) of clusters was eight (Figure 5A). Using K = 8 in STRUCTURE, we found that samples from Indonesia, the Philippines and Taiwan each formed a unique, distinct cluster. On the other hand, lineages from mainland China and Japan were difficult to separate (Figure 5B). Other values of K were explored, without providing any additional insights (not shown). This clustering was consistent with the PCoA (Figure 3), i.e., both analyses recognized distinct lineages from Indonesia, the Philippines and Taiwan, but not within mainland China and Japan. Populations from mountain regions (Sichuan and Yunnan) were not distinguished from those in the lake regions.
To our knowledge, this report includes the widest geographic coverage to date, surveying 12 S. japonicum populations in four Asian countries to assess the degree of genetic subdivision within the species.
The two classes of markers, mitochondrial DNA and microsatellites, yielded results those were not completely congruent. Based on the phylogenetic tree (constructed using three mitochondrial DNA fragments), the Taiwanese lineage obviously diverged from the others. Moreover, the lineages from the mountain regions of mainland China, Indonesia, Japan and the Philippines were distinct from each other and from populations in the lake region. However, no substructure was detected among the lineages from the Chinese lake region. In agreement with the phylogenetic analysis, the PCoA analysis, based on ten microsatellite loci, showed that individuals from the Philippines, Indonesia and Chinese Taiwan formed distinct clusters according to their geographical origins. This pattern was further supported by a Bayesian approach in STRUCTURE. The S. japonicum lineages from the mountain and lake regions of mainland China and from Japan, that were distinct according to the mitochondrial DNA markers, could not be separated by the microsatellite analysis.
Although some earlier studies failed to detect any evidence of major genetic dissimilarities at the DNA level between Chinese and Philippine S. japonicum , our study, in agreement with the more recent studies using microsatellites , has shown that schistosome populations in the two countries are genetically distinct. Moreover, although S. japonicum was first described from Taiwan in 1914 , there have been few studies investigating genetic variability in this population . Here, we provide the first report that the Taiwanese population is clearly divergent from others.
Genetic divergence observed among regional Chinese and Asian populations in this study may relate to the distribution and morphs of Oncomelania hupensis, the only intermediate host species complex involved in the transmission of S. japonicum . In mainland China, there exist two distinct morphological and allozyme forms of O. hupensis: one with a smooth shell in the mountain region and another with a ribbed shell in the lake region [32,33]. A single phenotype of O. hupensis quadrasi occurs in the Philippines, while O. h. nosophora is present in Japan, O. h. lindoensis in the Indonesian islands, and O.h. formosana and O. h. chiui are present in Chinese Taiwan . Oncomelania hupensis from China is genetically different from O. h. quadrasi in the Philippines , and the O. quadrasi populations in the Philippines have a substructure associated with their geographic origin . These different phenotypic and genotypic morphs of O. hupensis that exist in different regions might contribute to the genetic divergences apparent among populations of its parasite, S. japonicum .
In the present study, F ST values ranged from 0.01 to 0.51, indicating varied levels of pairwise population-genetic differentiation. In fact, strong population genetic differentiation has been detected between the mountain and lake regions of China on a large geographical scale [3,10]. This could be explained by different ecological environments and agricultural practices between the two habitat types. Another explanation may relate to the topographical isolation of the mountainous land from the lower Yangtze River basin, resulting in restrictions to gene flow. In agreement with a previous study which applied another set of microsatellite loci , our Mantel test showed no significant correlation between genetic and geographical distance. However, identical multilocus genotypes were not shared among the populations, which would be direct evidence of gene flow.
Corresponding well with the previous studies indicating high genetic diversity of S. japonicum, based on microsatellites  and mitochondrial DNA , in this study high genetic diversity was detected within most S. japonicum populations, according to both mitochondrial DNA and microsatellite markers. We also observed marked deviations from HWE in most populations, in agreement with findings of a previous study on S. japonicum  and similar to the results of studies on S. mansoni . This significant deviation could have resulted from the effects of non-random mating, inbreeding and population subdivision. Also, the deviation could also have been caused by unintended non-random sampling, given the way in which specimens for analysis were obtained. Furthermore, only low numbers of individuals were available for analysis from some populations.
In conclusion, high genetic diversity within and substantial genetic differentiation among S. japonicum populations were demonstrated using both mitochondrial DNA and microsatellites. Future research should aim to generate a comprehensive link between the population nucleotide diversity/structure differences and corresponding phenotypic variation, such as host preferences and differences in pathology, as well as possible differences in the response to anti-schistosome vaccines and immunodiagnostics. Overall, the present findings provide fundamental biological and evolutionary information on S. japonicum, with significant implications for the control and elimination of this parasite in Asia.
Steinmann P, Keiser J, Bos R, Tanner M, Utzinger J. Schistosomiasis and water resources development: systematic review, meta-analysis, and estimates of people at risk. Lancet Infect Dis. 2006;6(7):411–25.
King CH, Dickman K, Tisch DJ. Reassessment of the cost of chronic helmintic infection: a meta-analysis of disability-related outcomes in endemic schistosomiasis. Lancet. 2005;365(9470):1561–9.
Zhao QP, Sen Jiang M, Dong HF, Nie P. Diversification of Schistosoma japonicum in Mainland China Revealed by Mitochondrial DNA. Plos Neglect Trop D. 2012;6(2):e1503.
Zhou XN, Guo JG, Wu XH, Jiang QW, Zheng J, Dang H, et al. Epidemiology of schistosomiasis in the People’s Republic of China, 2004. Emerg Infect Dis. 2007;13(10):1470–6.
Zhou XN, Bergquist R, Leonardo L, Yang GJ, Yang K, Sudomo M, et al. Schistosomiasis Japonica: Control and Research Needs. In: Zhou XN, Bergquist R, Olveda R, Utzinger J, editors. Advances in Parasitology, Vol 72: Important Helminth Infections in Southeast Asia: Diversity and Potential for Control and Elimination, Pt A, vol. 72. 2010. p. 145–78.
Hope M, Foley DH, McManus DP. Electrophoretically detected allozyme variation reveals only moderate differentiation between Chinese and Philippine Schistosoma japonicum. Acta Trop. 1995;60(2):101–8.
Bowles J, Hope M, Tiu WU, Liu XS, McManus DP. Nuclear and mitochondrial genetic markers highly conserved between Chinese and Philippine Schistosoma japonicum. Acta Trop. 1993;55(4):217–29.
Sorensen E, Drew AC, Brindley PJ, Bogh HO, Gasser RB, Qian BZ, et al. Variation in the sequence of a mitochondrial NADH dehydrogenase I gene fragment among six natural populations of Schistosoma japonicum from China. Int J Parasitol. 1998;28(12):1931–4.
Gasser RB, Qian BZ, Nansen P, Johansen MV, Bogh H. Use of RAPD for the detection of genetic variation in the human blood fluke, Schistosoma japonicum, from mainland China. Mol Cell Probes. 1996;10(5):353–8.
Shrivastava J, Qian BZ, McVean G, Webster JP. An insight into the genetic variation of Schistosoma japonicum in mainland China using DNA microsatellite markers. Mol Ecol. 2005;14(3):839–49.
Zhao GH, Mo XH, Zou FC, Li J, Weng YB, Lin RQ, et al. Genetic variability among Schistosoma japonicum isolates from different endemic regions in China revealed by sequences of three mitochondrial DNA genes. Vet Parasitol. 2009;162(1–2):67–74.
Rudge JW, Carabin H, Balolong Jr E, Tallo V, Shrivastava J, Lu DB, et al. Population Genetics of Schistosoma japonicum within the Philippines Suggest High Levels of Transmission between Humans and Dogs. Plos Neglect Trop D. 2008;2(11):e340.
Woodruff DS, Merenlender AM, Upatham ES, Viyanant V. Genetic variation and differentiation of three Schistosoma species from the Philippines, Laos and Peninsular Malaysia. Am J Trop Med Hyg. 1987;36(2):345–54.
Avise J. Phylogeography: the history and formation of species. Cambridge, MA: Harvard University Press; 2000.
Avise J. Molecular markers, natural history and evolution. 2nd ed. Sunderland, MA: Sinuate; 2004.
Allendorf FW, Seeb LW. Concordance of genetic divergence among sockeye salmon populations at allozyme, nuclear DNA, and mitochondrial DNA markers. Evolution. 2000;54(2):640–51.
Lu G, Basley DJ, Bernatchez L. Contrasting patterns of mitochondrial DNA and microsatellite introgressive hybridization between lineages of lake whitefish (Coregonus clupeaformis); relevance for speciation. Mol Ecol. 2001;10(4):965–85.
Ventura M, Petrusek A, Miro A, Hamrova E, Bunay D, De Meester L, et al. Local and regional founder effects in lake zooplankton persist after thousands of years despite high dispersal potential. Mol Ecol. 2014;23(5):1014–27.
Cummings MP, Otto SP, Wakeley J. Sampling properties of DNA sequence data in phylogenetic analysis. Mol Biol Evol. 1995;12(5):814–22.
Rieseberg LH. Molecular ecology of hybridization. In: Advances in molecular ecology. Amsterdam: Ios Press; 1998.
Yin MB, Hu W, Mo XJ, Wang SY, Brindley PJ, McManus DP, et al. Multiple near-identical genotypes of Schistosoma japonicum can occur in snails and have implications for population-genetic analyses. Int J Parasitol. 2008;38(14):1681–91.
Excoffier L, Laval G, Schneider S. Arlequin (version 3.0): An integrated software package for population genetics data analysis. Evol Bioinform. 2005;1:47–50.
Tamura K, Dudley J, Nei M, Kumar S. MEGA4: Molecular evolutionary genetics analysis (MEGA) software version 4.0. Mol Biol Evol. 2007;24(8):1596–9.
Bandelt HJ, Forster P, Rohl A. Median-joining networks for inferring intraspecific phylogenies. Mol Biol Evol. 1999;16(1):37–48.
Librado P, Rozas J. DnaSP v5: a software for comprehensive analysis of DNA polymorphism data. Bioinformatics. 2009;25(11):1451–2.
Rousset F. GENEPOP ‘ 007: a complete re-implementation of the GENEPOP software for Windows and Linux. Mol Ecol Resour. 2008;8(1):103–6.
Pritchard JK, Stephens M, Donnelly P. Inference of population structure using multilocus genotype data. Genetics. 2000;155(2):945–59.
Evanno G, Regnaut S, Goudet J. Detecting the number of clusters of individuals using the software STRUCTURE: a simulation study. Mol Ecol. 2005;14(8):2611–20.
Fan P-C. The history of schistosomiasis japonica in Taiwan. Kaohsiung J Med Sci. 2006;22(7):309–20.
Anou N, Yanwen X, Youren F. Schistosoma japonicum strains: Differentiation by RAPD and SSR-PCR. Southeast Asian J Trop Med Public Health. 2002;33(4):720–4.
Zhao QP, Sen Jiang M, Littlewood DTJ, Nie P. Distinct Genetic Diversity of Oncomelania hupensis, Intermediate Host of Schistosoma japonicum in Mainland China as Revealed by ITS Sequences. Plos Neglect Trop D. 2010;4(3):e611.
Liu Y, Lou T, Wang Y, Zhang W. Subspecies difference of oncomelaniid snails. Acta Zootaxonm Sin. 1981;6:253–67.
Zhou XN. Population genetics morphology and distribution of Oncomelania hupensis, intermediate host of Schistosoma japonicum in mainland China. Copenhagen: University of Copenhagen; 1994.
Saijuntha W, Jarilla B, Leonardo A, Sunico L, Leonardo L, Andrews R, et al. Genetic Structure Inferred from Mitochondrial 12S Ribosomal RNA Sequence of Oncomelania quadrasi, the Intermediate Snail Host of Schistosoma japonicum in the Philippines. Am J Trop Med Hyg. 2014;90(6):1140–5.
Seto E, Xu B, Liang S, Gong P, Wu WP, Davis G, et al. The use of remote sensing for predictive modeling of schistosomiasis in China. Photogramm Eng Remote Sensing. 2002;68(2):167–74.
Rudge JW, Lu DB, Fang GR, Wang TP, Basanez MG, Webster JP. Parasite genetic differentiation by habitat type and host species: molecular epidemiology of Schistosoma japonicum in hilly and marshland areas of Anhui Province. China Mol Ecol. 2009;18(10):2134–47.
Thiele EA, Sorensen RE, Gazzinelli A, Minchella DJ. Genetic diversity and population structuring of Schistosoma mansoni in a Brazilian village. Int J Parasitol. 2008;38(3–4):389–99.
This research was funded by National Natural Science Foundation of China (No. 91431104) and Foundation of National Science and Technology Major Program (Grant Nos. 2012ZX10004-220 and 2008ZX10004-011).
The authors declare that they have no competing interests.
MY and WH designed the study, HL and BX carried out the molecular work. MY, ZY and HL contributed to data analyses. MY, DM, DB, ZF and WH wrote the manuscript. All authors read and approved the final version.
Mingbo Yin and Hongyan Li contributed equally to this work.
Median-joining Network based on the mitochondrial haplotypes (concatenated nad1 + nad4 + 16-12S) of S. japonicum. Abbreviations: AHGC: Guichi Country, Anhui Province; AHTL: Tongling Country, Anhui Province; HBSH: Shashi City, Hubei Province; HNCD: Changde City, Hunan Province; HNYY: Yueyang City Hunan Province; JXDC: Duchang Country, Jiangxi Province; SCXC: Xichang City, Sichuan Province; YNEY: Dali City, Eryuan Country, Yunnan Province; CTW: Chinese Taiwan; IN: Indonesia JP: Japan; PH: The Philippines. The sizes of nodes are proportional to the frequencies of the haplotypes.