Temporal genetic diversity of Schistosoma japonicum in two endemic sites in China revealed by microsatellite markers
Parasites & Vectors volume 9, Article number: 36 (2016)
Schistosomiasis is one of the neglected tropical diseases. The causative agent of schistosomiasis in China, Schistosoma japonicum, has long been a major public health problem. An understanding of fundamental evolutionary and genetic processes in this species has major implications for its control and elimination. Intensive control efforts have greatly reduced the incidence of schistosomiasis in China, but little is known about the genetic consequences of these efforts.
To investigate this, we sampled twice (years 2003 and 2011) from two endemic regions where populations of S. japonicum had persisted despite control efforts and genotyped these samples using ten microsatellite markers. Our main hypothesis was that parasite genetic diversity would be greatly reduced across this time period.
There was no apparent reduction in allelic diversity, and a non-significant reduction in clonal diversity in these parasite populations between 2003 and 2011. We did, however, detect temporal genetic differentiation among the samples. Such a significant temporal genetic variation of S. japonicum populations has not been reported before.
Schistosomiasis is a serious parasitic disease, infecting over 200 million people and threatening the health of about 779 million people around the world . Schistosomiasis due to infection with Schistosoma japonicum has long been a major public health problem in China . In the last decade, the disease has been intensively controlled in most endemic regions in China . Specifically, in 2004, a national program for schistosomiasis control was launched, aiming to reduce the transmission of S. japonicum from cattle and humans to snails. Effective interventions were implemented from 2005 to 2007, including keeping cattle away from snail-infested grasslands, providing farmers with mechanised farming tools and improving public environmental sanitation by supplying tap water and building lavatories. Annual synchronous rounds of chemotherapy are also routinely used . It has been proposed that evolutionary theory should have an important role in the design, application and interpretation of such programs . Control interventions are expected to reduce the genetic diversity within the parasite population. Associated with this, intensive selection is likely, driving rapid evolutionary changes in natural parasite populations . For example, under strong drug selective pressure, resistance could emerge and become fixed, since the drug-resistant phenotypes and alleles have an advantage over the wild type in these adverse circumstances . Although we cannot test for such selection here, we can look for evidence of reduced genetic diversity in parasite populations following intervention.
Genetic diversity is of great importance for any organism to adapt in a changing environment and, in the case of parasites, to respond to the pressure of intervention programs . Different types of molecular markers, for instance, mitochondrial DNAs  and microsatellites , have been applied to evaluate genetic variation in S. japonicum populations among different geographical regions , as well as the genetic differentiation among host species  or host individuals . Some recent studies have reported spatial genetic variation [10, 12, 13] among natural S. japonicum populations. There have also been reports on spatio-temporal modeling for the prevalence of S. japonicum infection [14, 15]. However, to date, no study has investigated temporal changes in genetic structure of S. japonicum populations.
In this study, we sampled again in 2011 from seven endemic regions where S. japonicum populations were present in 2003. In 2011, the parasite was only detected in two of the seven locations, suggesting a high degree of intervention success at several locations. For the two locations at which parasites had persisted, we hypothesised that a considerable reduction in genetic diversity (allelic richness and clonal diversity) would have occurred between these time points as a consequence of the interventions. In addition, we applied a range of population-genetic analyses under the hypothesis that there should be no temporal changes between the time points.
All procedures involving animals were carried out according to the guidelines of the Association for Assessment and Accreditation of Laboratory Animal Care International. Our 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).
In 2011, we re-sampled seven endemic sites in China that had been investigated in 2003 for our previous project . In the interim, the National Program for Schistosomiasis Control had reduced the transmission of S. japonicum  to such an extent that infected snails were detected in only two out of seven locations: Yueyang city in Hunan Province and Shashi city in Hubei Province (Table 1). These two sites are both close to the Yangtze River and about 129 km apart.
Ten laboratory-raised Kunming mice were each infected with about 1,000 S. japonicum cercariae isolated from infected snails (Oncomelania hupensis) collected from two locations (i.e. Yueyang and Shashi) in 2011. For each location, we screened 5,000 snails in 2011 and detected three infected snails from Shashi and 25 infected snails from Yueyang. After 45 days, adult worms were obtained by perfusing the hepatic portal system and mesenteric veins of each infected mouse. All worms were washed at least three times with normal saline (0.9 % NaCl) to remove the host tissues before being stored in 95 % ethanol at 4 °C. The gender of each worm was noted.
DNA extraction, PCR amplification and microsatellite genotyping
In this study, we randomly selected 44 worms per location from the year 2011 samples for the genetic analyses. Genomic DNA was extracted individually from adult schistosomes using a DNeasy Blood & Tissue Kit (QIAGEN, Germany), following the manufacturer’s instructions, and then stored at -20 °C until use.
The DNA of each individual S. japonicum was genotyped at ten microsatellite loci (i.e. Sjp1, 4, 5, 6, 8, 9, 10, 14, 15 and 17) from our previous work . The PCRs were performed in a total volume of 12 μl following the protocol in . Then the PCR products were diluted to appropriate concentrations and analysed on an ABI 3730 capillary automated sequencer using a LIZ 500-labeled size standard. Allele sizes were read using GeneMapper software (version 4.0) and checked carefully.
To compare the genetic variation of S. japonicum populations between 2003 and 2011, we included microsatellite data from 51 worms from Yueyang and 23 worms from Shashi, which were collected in 2003 and published previously . This yielded four population categories for analysis (two per time point and per location). For some analyses, these categories were further subdivided by sex (termed subpopulations here). We excluded all individuals for which data were missing.
Clonal diversity (R) was calculated manually as R = (G-1)/(N-1) , where G is the number of genotypes and N represents sample size. Clonal diversity can vary from 0 (all individuals belong to one clone) to 1 (each individual is unique). All 162 individuals were included in this analysis. Principal coordinates analyses (PCoA) were then performed to display the relatedness of all unique multi-locus genotypes (MLGs) in GenAlEx 6.5 . Somatic mutations occurring during the proliferation of cercariae from a single miracidium may produce multiple near-identical MLGs . Inclusion of clusters of near-identical siblings may bias results. We therefore removed near-identical MLGs in the subsequent analyses, retaining only a single representative of each cluster . PCoAs were then done to show the relatedness of the 82 retained MLGs. All subsequent analyses were based on these 82 MLGs. Genetic diversity indices (allelic richness A r ; sample size was corrected) were estimated using FSTAT V126.96.36.199 . As the sample size was low, males and females were not treated separately in the calculations of allelic richness.
Genetic differentiation between pairs of samples was calculated as Wright’s F ST in Arlequin 3.11 , and the significance of the F ST value was tested using 104 permutations. Hierarchical analysis of molecular variance (AMOVA) was applied to quantify genetic variances into spatial, temporal and within-sample components in Arlequin 3.11.
In total, 138 unique MLGs were detected among 162 S. japonicum individuals based on ten microsatellite loci. No shared MLG was found among the eight subpopulations (location, year and gender). The PCoA revealed that the 138 unique MLGs fell into a main cluster (Cluster 1) with small outlying clusters (Fig. 1a). Most (6 of 8) MLGs in the male subpopulation from Yueyang in 2011 (YY11M) were assigned to a small cluster (Cluster 2), whereas the other two were in Cluster 1 (Fig. 1a). Additionally, some (8 of 23) unique MLGs in the female subpopulation from Yueyang in 2011 (YY11F) formed another small cluster (Cluster 3) while the remainder fell into Cluster 1 (Fig. 1a). Overall, 56 individuals exhibited near-identical MLGs (Table 1), indicating that they might be siblings. When only a single representative of each near-identical MLG group was retained, no visually obvious temporal or spatial structuring was apparent among the remaining 82 MLGs (Fig. 1b).
Importantly, allelic richness among the 82 retained MLGs and corrected for sample size did not change dramatically at either site between 2003 and 2011 (Table 1). Indeed, allelic richness at Shashi was greater in 2011 than in 2003. The clonal diversity (R) of subpopulations ranged from 0.35 to 1.00 (average R = 0.85). The value of 1.00 (detected in three subpopulations; Table 1) indicated that all the individuals possess different MLGs. The clonal diversity of the 2011 male subpopulation in Yueyang was particularly low (0.35), compared with 0.95 in 2003. By 2011, the overall clonal diversity had decreased (average R = 0.76), but not significantly (P = 0.06).
Genetic differentiation was significant for all tested pairs of four S. japonicum samples and the F ST values ranged from 0.009 to 0.046 (averaged over all loci, Table 2). AMOVA analyses indicated that within-population variance explained most of the observed genetic variation (97.45 %). The remaining genetic variation was due to additional variance across space (0.02 %) and across time within space (2.53 %; Table 3). The temporal component was significant but spatial was not (Table 3).
Using microsatellite markers, we noted a non-significant decrease in clonal diversity between the two time periods, but no obvious change in allelic richness. Measures of gene flow and population structure indicated substantial temporal genetic change (between 2003 and 2011) according to population differentiation tests (F ST) and the AMOVA. Spatial structure (between Yueyang and Shashi) was also indicated according to F ST values, but not according to the AMOVA results.
Clonal diversities were in agreement with previous studies using microsatellites [8, 20]. A high clonal diversity (average R = 0.94) was detected in S. japonicum populations from the year 2003. However, by 2011, this had decreased but not significantly. The clonal diversity in subpopulation YY11M was very low (0.35), perhaps indicating that the sample contained many sibling cercariae from a single male miracidium. Such a temporal change in clonal diversity of S. japonicum populations has not been observed previously.
The maintenance of high allelic richness across time was not expected. Infected snails were more difficult to find in 2011 than in 2003 (and indeed were not found at all in five of seven locations in 2011). This suggests a reduction in overall population size of S. japonicum in the region, but apparently not to an extent sufficient to impact markedly on allelic richness. Similarly, the extent of temporal population-genetic structure was unexpected. Such temporal changes are not previously known for S. japonicum, although structure at small spatial scales has been reported [8, 11, 16]. Little is known about the metapopulation dynamics of S. japonicum over the vast region of the mid-Yangtze basin, from which our samples came. Local extinction followed by repopulation from elsewhere in the region might partly explain the extent of genetic differentiation observed. Given the known fluxes of water in the region , especially in connection with annual floods, there may be a large regional population of S. japonicum subject to substantial dispersal and gene flow.
At our study locations, intensive control measures (annual chemotherapy) had been implemented during the investigation period (spanning 9 years), which might have acted as a strong selective pressure and been responsible for the temporal changes observed. There were also marked ecological changes during this period, in particular the completion of the Three Gorges Dam (TGD), which significantly reduced the amount of water reaching our sampling locations (Dongting and Poyang lakes) and has moderated the annual floods. Given the lack of baseline data from prior to the construction of the TGD, and the lack of understanding of metapopulation dynamics in S. japonicum, few strong conclusions can be drawn.
There are several limitations of our study. One is the use of different mammal hosts between 2003 (rabbits) and 2011 (mice). LoVerde et al.  noted that use of rodents (as opposed to primates) as experimental hosts led to increased reduction in genetic diversity over several generations in S. mansoni. However, we used mammal hosts only for a single passage (field-collected snails were used to infect laboratory hosts from which adult worms were harvested) and not serial passages such as used by LoVerde et al. . In effect, our worms are analogous to those subjected to electrophoresis from the M0 generation of laboratory hosts, in which very little deviation from the parental population was observed . Another limitation relates to the small number of infected snails (three) found in Shashi in 2011. Despite this, the clonal and allelic diversities within the parasite population there were high. The inference is that at least one of the infected snails had been invaded by multiple unrelated miracidia. A further point to consider is that adult S. japonicum can be long-lived in the mammal host, 47 years being reported in one human case . Long-lived individuals can carry substantial genetic diversity through a bottleneck event. The interventions have taken place over a relative handful of years, perhaps insufficient time for any genetic bottleneck to have caused genetic erosion. Similarly, the effects of regional ecological change, such as caused by the TGD, may take more time to become apparent in the genetic structure of schistosome populations.
In conclusion, using ten microsatellite loci, we rejected the hypothesis that genetic diversity in two S. japonicum populations would be reduced following control programs between the years of 2003 and 2011. We did, however, find significant genetic differentiation across time in these populations. Overall, our findings provide a deeper understanding of molecular epidemiology and population genetics of S. japonicum, which may be of value for effective control and elimination of schistosomiasis.
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.
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.
Li SZ, Luz A, Wang XH, Xu LL, Wang Q, Qian YJ, et al. Schistosomiasis in China: acute infections during 2005-2008. Chin Med J. 2009;122(9):1009–14.
Wang LD, Guo JG, Wu XH, Chen HG, Wang TP, Zhu SP, et al. China’s new strategy to block Schistosoma japonicum transmission: experiences and impact beyond schistosomiasis. Trop Med Int Health. 2009;14(12):1475–83.
Webster JP, Gower CM, Norton AJ. Evolutionary concepts in predicting and evaluating the impact of mass chemotherapy schistosomiasis control programmes on parasites and their hosts. Evol Appl. 2008;1(1):66–83.
Rollinson D, Kaukas A, Johnston DA, Simpson AJG, Tanaka M. Some molecular insights into Schistosome evolution. Int J Parasitol. 1997;27(1):11–28.
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.
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.
Chen F, Li J, Sugiyama H, Zhou DH, Song HQ, Zhao GH, et al. Genetic variability among Schistosoma japonicum isolates from the Philippines, Japan and China revealed by sequence analysis of three mitochondrial genes. Mitochondrial DNA. 2015;26(1):35–40.
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.
Yin M, Zheng HX, Su J, Feng Z, McManus DP, Zhou XN, et al. Co-dispersal of the blood fluke Schistosoma japonicum and Homo sapiens in the Neolithic age. Sci Rep. 2015;5:18058.
Zhou YB, Liang S, Chen GX, Rea C, Han SM, He ZG, et al. Spatial-temporal variations of Schistosoma japonicum distribution after an integrated national control strategy: a cohort in a marshland area of China. BMC Public Health. 2013;13:297.
Yang GJ, Vounatsou P, Zhou XN, Tanner M, Utzinger J. A Bayesian-based approach for spatio-temporal modeling of county level prevalence of Schistosoma japonicum infection in Jiangsu province, China. Int J Parasitol. 2005;35(2):155–62.
Yang K, Zhou XN, Wang XH, Yang GJ, Jia TW, He WL. Establishment and application of spatio-temporal model of Schistosomiasis japonica in a county in marshland region. Chin J Schistosom Control. 2012;24(1):14–20.
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.
Collins C, Xu J, Tang S. Schistosomiasis control and the health system in PR China. Infect Dis Poverty. 2012;1:8.
Dorken ME, Eckert CG. Severely reduced sexual reproduction in northern populations of a clonal plant, Decodon verticillatus (Lythraceae). J Ecol. 2001;89(3):12.
Peakall R, Smouse PE. GENALEX 6: genetic analysis in Excel. Population genetic software for teaching and research. Mol Ecol Notes. 2006;6(1):288–95.
Yin M, Li H, McManus DP, Blair D, Su J, Yang Z, et al. Geographical genetic structure of Schistosoma japonicum revealed by analysis of mitochondrial DNA and microsatellite markers. Parasite Vector 2015;8:150.
Goudet J. FSTAT (Version 1.2): a computer program to calculate F-statistics. J Hered. 1995;86(6):485–6.
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.
Yu F, Chen Z, Ren X, Yang G. Analysis of historical floods on the Yangtze River, China: characteristics and explanations. Geomorphology. 2009;113(3-4):210–6.
LoVerde PT, Dewald J, Minchella DJ, Bosshardt SC, Damian RT. Evidence for host-induced selection in Schistosoma mansoni. J Parasitol. 1985;71(3):297–301.
Hall SC, Kehoe EL. Prolonged survival of Schistosoma japonicum. Calif Med. 1970;113(2):75–7.
This research was funded by the National Natural Science Foundation of China (No. 91431104) and Foundation of National Science, International S&T Cooperation Program of China (2014DFA31130) 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, DB and HL contributed to data analyses. MY, DB and WH wrote the manuscript. All authors read and approved the final version.
About this article
Cite this article
Yin, M., Li, H., Blair, D. et al. Temporal genetic diversity of Schistosoma japonicum in two endemic sites in China revealed by microsatellite markers. Parasites Vectors 9, 36 (2016). https://doi.org/10.1186/s13071-016-1326-7
- Schistosoma japonicum
- Genetic diversity
- Population structure
- Temporal variation