An insight into the genetic variation of Schistosoma japonicum in mainland China using DNA microsatellite markers

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.

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

BACKGROUND

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

METHODS

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

RESULTS

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

CONCLUSIONS

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

Construction and characterization of microsatellite markers for the Schistosoma japonicum isolate from a hilly area of China based on whole genome sequencing

Schistosoma japonicum had once caused the greatest disease burden in China and has still been transmitted in some hilly areas, for example, in Shitai of Anhui province, where rodents are projected to be the main reservoir. This may lead to a critical need of molecular tools with high efficiency in monitoring the dynamic of the rodent-associated S. japonicum, as an appropriate amount of schistosome input can re-establish its life cycle in a place with snails and then result in the re-emergence of schistosomiasis. Therefore, the goal of this study was to develop high polymorphic microsatellites from the whole genome of rodent-associated S. japonicum strain to monitor its transmission dynamic. We sampled the hilly schistosome isolate from Shitai of Anhui in China and sequenced the parasite with the next-generation sequencing technology. The whole genome was assembled with four different approaches. We then developed 71 microsatellite markers at a genome-wide scale throughout two best assembled genomes. Based on their chromosome mapping and the expected length of targeted sequences, we selected 24 markers for the development of multiplex reactions. Two multiplexes composed of 10 loci were finally developed, and their potential was revealed by their successful application on and capturing the genetic diversity of three schistosome populations. The selected 10 markers, each with clear chromosome location and characteristics, will be greatly useful in tracing the dispersal pathways or/and dynamics of the rodent-associated S. japonicum or others in the hilly area of China or elsewhere.

Evolution of tetraspanin antigens in the zoonotic Asian blood fluke Schistosoma japonicum

BACKGROUND

Despite successful control efforts in China over the past 60 years, zoonotic schistosomiasis caused by Schistosoma japonicum remains a threat with transmission ongoing and the risk of localised resurgences prompting calls for a novel integrated control strategy, with an anti-schistosome vaccine as a core element. Anti-schistosome vaccine development and immunisation attempts in non-human mammalian host species, intended to interrupt transmission, and utilising various antigen targets, have yielded mixed success, with some studies highlighting variation in schistosome antigen coding genes (ACGs) as possible confounders of vaccine efficacy. Thus, robust selection of target ACGs, including assessment of their genetic diversity and antigenic variability, is paramount. Tetraspanins (TSPs), a family of tegument-surface antigens in schistosomes, interact directly with the host's immune system and are promising vaccine candidates. Here, for the first time to our knowledge, diversity in S. japonicum TSPs (SjTSPs) and the impact of diversifying selection and sequence variation on immunogenicity in these protiens were evaluated.

METHODS

SjTSP sequences, representing parasite populations from seven provinces across China, were gathered by baiting published short-read NGS data and were analysed using in silico methods to measure sequence variation and selection pressures and predict the impact of selection on variation in antigen protein structure, function and antigenic propensity.

RESULTS

Here, 27 SjTSPs were identified across three subfamilies, highlighting the diversity of TSPs in S. japonicum. Considerable variation was demonstrated for several SjTSPs between geographical regions/provinces, revealing that episodic, diversifying positive selection pressures promote amino acid variation/variability in the large extracellular loop (LEL) domain of certain SjTSPs. Accumulating polymorphisms in the LEL domain of SjTSP-2, -8 and -23 led to altered structural, functional and antibody binding characteristics, which are predicted to impact antibody recognition and possibly blunt the host's ability to respond to infection. Such changes, therefore, appear to represent a mechanism utilised by S. japonicum to evade the host's immune system.

CONCLUSION

Whilst the genetic and antigenic geographic variability observed amongst certain SjTSPs could present challenges to vaccine development, here we demonstrate conservation amongst SjTSP-1, -13 and -14, revealing their likely improved utility as efficacious vaccine candidates. Importantly, our data highlight that robust evaluation of vaccine target variability in natural parasite populations should be a prerequisite for anti-schistosome vaccine development.

Biophysical description of Bromosulfophthalein interaction with the 28-kDa glutathione transferase from Schistosoma japonicum

Glutathione transferases (GSTs) are major detoxification enzymes vital for the survival and reproduction of schistosomes during infection in humans. Schistosoma encode two GST isoenzymes, the 26- and 28-kDa isoforms, that show different substrate specificities and cellular localisations. Bromosulfophthalein (BSP) has been identified and characterised as a potent 26-kDa Schistosoma japonicum GST (Sj26GST) inhibitor with an anthelmintic potential. This study describes the structure, function, and ligandin properties of the 28-kDa Schistosoma japonicum GST (Sj28GST) towards BSP. Enzyme kinetics show that BSP is a potent enzyme inhibitor, with a specific activity decreases from 60.4 µmol/min/mg to 0.0742 µmol/min/mg and an IC₅₀ in the micromolar range of 0.74 µM. Far-UV circular dichroism confirmed that purified Sj28GST follows a typical GST fold, which is predominantly alpha-helical. Fluorescence spectroscopy suggests that BSP binding occurs at a site distinct from the glutathione-binding site (G-site); however, the binding does not alter the local G-site environment. Isothermal titration calorimetry studies show that the binding of BSP to Sj28GST is exergonic (∆G°= -33 kJ/mol) and enthalpically-driven, with a stoichiometry of one BSP per dimer. The stability of Sj28GST (∆G_(H2O) = 4.7 kcal/mol) is notably lower than Sj26GST, owing to differences in the enzyme's dimeric interfaces. We conclude that Sj28GST shares similar biophysical characteristics with Sj26GST based on its kinetic properties and susceptibility to low concentrations of BSP. The study supports the potential benefits of re-purposing BSP as a potential drug or prodrug to mitigate the scourge of schistosomiasis.

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.

Integrating genomic and epidemiologic data to accelerate progress toward schistosomiasis elimination

The global community has adopted ambitious goals to eliminate schistosomiasis as a public health problem, and new tools are needed to achieve them. Mass drug administration programs, for example, have reduced the burden of schistosomiasis, but the identification of hotspots of persistent and reemergent transmission threaten progress toward elimination and underscore the need to couple treatment with interventions that reduce transmission. Recent advances in DNA sequencing technologies make whole-genome sequencing a valuable and increasingly feasible option for population-based studies of complex parasites such as schistosomes. Here, we focus on leveraging genomic data to tailor interventions to distinct social and ecological circumstances. We consider two priority questions that can be addressed by integrating epidemiological, ecological, and genomic information: (1) how often do non-human host species contribute to human schistosome infection? and (2) what is the importance of locally acquired versus imported infections in driving transmission at different stages of elimination? These questions address processes that can undermine control programs, especially those that rely heavily on treatment with praziquantel. Until recently, these questions were difficult to answer with sufficient precision to inform public health decision-making. We review the literature related to these questions and discuss how whole-genome approaches can identify the geographic and taxonomic sources of infection, and how such information can inform context-specific efforts that advance schistosomiasis control efforts and minimize the risk of reemergence.

The genetic variation of different developmental stages of Schistosoma japonicum: do the distribution in snails and pairing preference benefit the transmission?

BACKGROUND

Schistosoma japonicum is a waterborne parasite that causes schistosomiasis in humans and in more than 40 animal species. Schistosoma japonicum shows distinct genetic differentiation among geographical populations and multiple hosts, but the genetic diversity of different developmental stages of S. japonicum from is less studied. Such studies could elucidate ecological mechanisms in disease transmission by analysing feedbacks in individual physiology and population state.

METHODS

After infection using cercariae from a pool of snails shedding together (Method I) and infection using mixed equal numbers of cercariae from individually shed snails (Method II), different developmental stages of S. japonicum were genotyped with microsatellite loci, including 346 cercariae, 701 adult worms and 393 miracidia. Genetic diversity and molecular variation were calculated at different population levels. Kinships (I') among cercariae at intra-snail and inter-snail levels were evaluated. Genetic distance (Dsw) was compared between paired and unpaired worms, and partner changing was investigated through paternity identification for miracidia.

RESULTS

The cercaria clones in individual snails varied from 1 to 8 and the kinship of cercariae within individual snails was significant higher (P < 0.001) than that among different snails after deleting near-identical multi-locus genotypes (niMLGs). The allelic diversity of worms in Method I was lower (P < 0.001) than that in Method II, and allele frequency among mice in Method I was also less consistent. The parents of some miracidia were worms that were not paired when collected. The Dsw between each female of paired and unpaired males was much larger (P < 0.001) than that between the female and male in each pair.

CONCLUSIONS

Most of the infected snails contained multiple miracidia clones. The aggregation of genetically similar S. japonicum miracidia in individual snails and the unbalanced distribution of miracidia among snails suggests a non-uniform genetic distribution of cercariae among snails in the field. This further influenced the genetic structure of adult worms from infections with different cercariae sampling methods. Schistosoma japonicum in mice can change paired partner, preferring to mate with genetically similar worms. These characteristics provide implications for understanding the balance in genetic diversity of S. japonicum related to the transmission of schistosomiasis.

Genetic variability of Taenia solium cysticerci recovered from experimentally infected pigs and from naturally infected pigs using microsatellite markers

The adult Taenia solium, the pork tapeworm, usually lives as a single worm in the small intestine of humans, its only known definitive host. Mechanisms of genetic variation in T. solium are poorly understood. Using three microsatellite markers previously reported [1], this study explored the genetic variability of T. solium from cysts recovered from experimentally infected pigs. It then explored the genetic epidemiology and transmission in naturally infected pigs and adult tapeworms recovered from human carriers from an endemic rural community in Peru. In an initial study on experimental infection, two groups of three piglets were each infected with proglottids from one of two genetically different tapeworms for each of the microsatellites. After 7 weeks, pigs were slaughtered and necropsy performed. Thirty-six (92.3%) out of 39 cysts originated from one tapeworm, and 27 (100%) out of 27 cysts from the other had exactly the same genotype as the parental tapeworm. This suggests that the microsatellite markers may be a useful tool for studying the transmission of T. solium. In the second study, we analyzed the genetic variation of T. solium in cysts recovered from eight naturally infected pigs, and from adult tapeworms recovered from four human carriers; they showed genetic variability. Four pigs had cysts with only one genotype, and four pigs had cysts with two different genotypes, suggesting that multiple infections of genetically distinct parental tapeworms are possible. Six pigs harbored cysts with a genotype corresponding to one of the identified tapeworms from the human carriers. In the dendrogram, cysts appeared to cluster within the corresponding pigs as well as with the geographical origin, but this association was not statistically significant. We conclude that genotyping of microsatellite size polymorphisms is a potentially important tool to trace the spread of infection and pinpoint sources of infection as pigs spread cysts with a shared parental genotype.

Taenia solium, the pork tapeworm, is a major cause of epilepsy in developing countries. Although it has been deemed a potentially eradicable pathogen, it remains prevalent in rural communities. This two-part study aims to evaluate the utility of three microsatellite markers previously reported, to identify parasites and to establish relationships among them. In the first study, we evaluated the genetic variability of the progeny of two individual tapeworms by infecting two groups of three pigs each. We found variation of 8% and 0% in the two groups with respect to the parental tapeworm, indicating that the cysts source may be identifiable. Next, in the second study we described the genetic relationships among tapeworms obtained from four carriers and cysts obtained from eight naturally infected pigs in a rural community. We demonstrated that pigs can have two types of cysts, suggesting multiple infections. In addition, we found relatedness between 6 pigs and one tapeworm identified in the community. Our results indicate the potential for microsatellite markers to identify genetic relationships between parasites and thereby establish routes of transmission. It is likely that the limited number of microsatellites prevented us from establishing relatedness with more precision. Therefore, further evaluation of additional microsatellites is recommended.

Geographic strain differentiation of Schistosoma japonicum in the Philippines using microsatellite markers

BACKGROUND

Microsatellites have been found to be useful in determining genetic diversities of various medically-important parasites which can be used as basis for an effective disease management and control program. In Asia and Africa, the identification of different geographical strains of Schistosoma japonicum, S. haematobium and S. mansoni as determined through microsatellites could pave the way for a better understanding of the transmission epidemiology of the parasite. Thus, the present study aims to apply microsatellite markers in analyzing the populations of S. japonicum from different endemic areas in the Philippines for possible strain differentiation.

METHODOLOGY/ PRINCIPAL FINDINGS

Experimental mice were infected using the cercariae of S. japonicum collected from infected Oncomelania hupensis quadrasi snails in seven endemic municipalities. Adult worms were harvested from infected mice after 45 days of infection and their DNA analyzed against ten previously characterized microsatellite loci. High genetic diversity was observed in areas with high endemicity. The degree of genetic differentiation of the parasite population between endemic areas varies. Geographical separation was considered as one of the factors accounting for the observed difference between populations. Two subgroups have been observed in one of the study sites, suggesting that co-infection with several genotypes of the parasite might be present in the population. Clustering analysis showed no particular spatial structuring between parasite populations from different endemic areas. This result could possibly suggest varying degrees of effects of the ongoing control programs and the existing gene flow in the populations, which might be attributed to migration and active movement of infected hosts from one endemic area to another.

CONCLUSIONS/ SIGNIFICANCE

Based on the results of the study, it is reasonable to conclude that genetic diversity could be one possible criterion to assess the infection status in highly endemic areas. Genetic surveillance using microsatellites is therefore important to predict the ongoing gene flow and degree of genetic diversity, which indirectly reflects the success of the control program in schistosomiasis-endemic areas.

Genetic diversity and selection of three nuclear genes in Schistosoma japonicum populations

BACKGROUND

The blood fluke, Schistosoma japonicum still causes severe disease in China, the Philippines and Indonesia. Although there have been some studies the molecular epidemiology of this persistent and harmful parasite, few have explored the possibility and implications of selection in S. japonicum populations.

METHODS

We analyzed diversity and looked for evidence of selection at three nuclear genes (SjIpp2, SjFabp and SjT22.6) in 13 S. japonicum populations.

RESULTS

SjT22.6 was found to exhibit high nucleotide diversity and was under positive selection in the mountainous region of mainland China. As a tegumental protein, its secondary and tertiary structure differed between S. japonicum strains from the mountainous and lakes regions. In contrast, SjIpp2 and SjFabp had relatively low levels of nucleotide diversity and did not show significant departure from neutrality.

CONCLUSIONS

As a tegument-associated antigen-encoding gene of S. japonicum, SjT22.6 has high nucleotide diversity and appears to be under positive selection in the mountainous region of mainland China.

The basic reproductive ratio of Barbour's two-host schistosomiasis model with seasonal fluctuations

Background

Motivated by the first mathematical model for schistosomiasis proposed by Macdonald and Barbour’s classical schistosomiasis model tracking the dynamics of infected human population and infected snail hosts in a community, in our previous study, we incorporated seasonal fluctuations into Barbour’s model, but ignored the effect of bovine reservoir host in the transmission of schistosomiasis. Inspired by the findings from our previous work, the model was further improved by integrating two definitive hosts (human and bovine) and seasonal fluctuations, so as to understand the transmission dynamics of schistosomiasis japonica and evaluate the ongoing control measures in Liaonan village, Xingzi County, Jiangxi Province.

Methods

The basic reproductive ratio R₀ and its computation formulae were derived by using the operator theory in functional analysis and the monodromy matrix theory. The mathematical methods for global dynamics of periodic systems were used in order to show that R₀ serves as a threshold value that determines whether there was disease outbreak or not. The parameter fitting and the ratio calculation were performed with surveillance data obtained from the village of Liaonan using numerical simulation. Sensitivity analysis was carried out in order to understand the impact of R₀ on seasonal fluctuations and snail host control. The modified basic reproductive ratios were compared with known results to illustrate the infection risk.

Results

The Barbour’s two-host model with seasonal fluctuations was proposed. The implicit expression of R₀ for the model was given by the spectral radius of next infection operator. The R₀s for the model ranged between 1.030 and 1.097 from 2003 to 2010 in the village of Liaonan, Xingzi County, China, with 1.097 recorded as the maximum value in 2005 but declined dramatically afterwards. In addition, we proved that the disease goes into extinction when R₀ is less than one and persists when R₀ is greater than one. Comparisons of the different improved models were also made.

Conclusions

Based on the mechanism and characteristics of schistosomiasis transmission, Barbour’s model was improved by considering seasonality. The implicit formula of R₀ for the model and its calculation were given. Theoretical results showed that R₀ gave a sharp threshold that determines whether the disease dies out or not. Simulations concluded that: (i) ignoring seasonality would overestimate the transmission risk of schistosomiasis, and (ii) mollusiciding is an effective control measure to curtail schistosomiasis transmission in Xingzi County when the removal rate of infected snails is small.

Electronic supplementary material

The online version of this article (doi:10.1186/s13071-017-1983-1) 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.

An investigation into the potential effects of infrapopulation structure and other sources of sampling error, on population genetic studies of the transmission of Schistosoma japonicum (Trematoda: Digenea)

Background

Schistosoma japonicum remains a major challenge to human and animal health. Earlier microsatellite-based studies reported possible definitive-host-specific private alleles within S. japonicum, opening the possibility that different definitive hosts might harbour different parasite strains. Previous investigations have also detected near-identical multilocus genotypes in populations of adult worms - possibly the result of mutations occurring during the asexual (intramolluscan) phase of clonal expansion. Research has also revealed extensive deviations from Hardy-Weinberg Proportions (HWP) and conflicting results among studies. The present study was performed to examine some of the potential effects of infrapopulation structure on microsatellite-based studies of the transmission ecology of S. japonicum. Potential sources of bias considered included organotropic distribution of worms, non-random mating and corrections for clonal expansion.

Results

Stool samples from naturally infected hosts were used to infect snails in the laboratory and thereby expose mice. 274 individual worms were typed at seven microsatellite loci. Removal of individuals bearing duplicate MLGs (as a correction for presumed clonal expansion) had an impact on both HWP and organotropic genetic differentiation. The study found no evidence that heterozygote deficiencies were caused by a Wahlund effect. Female-male pairings appeared to be random and there was no evidence for mate choice by heterozygosity. There was some indication that excess heterozygosity, induced by clonal expansion, can offset heterozygote deficiencies caused by small population size or populations fragmented by parasite control efforts.

Conclusions

The view is supported that miracidia are preferable to adult worms in investigations into host-specific parasite lineages. Where adults must be used, extreme care should be taken with regard to sampling if infrapopulations of small animals are compared with those of larger animals; this is because of organotropic patterns in genetic variation and the tendency to sample from different organs in differently sized hosts. As corrections for clones may accentuate signals of population subdivision, corrections should only be made if tests for clonal expansion prove positive. Finally, evidence for heterozygote deficiency caused by small sample size, calls for carefully designed random and comprehensive sampling strategies for S. japonicum in China, where control efforts have greatly fragmented parasite populations.

Electronic supplementary material

The online version of this article (doi:10.1186/s13071-016-1454-0) contains supplementary material, which is available to authorized users.

Temporal genetic diversity of Schistosoma japonicum in two endemic sites in China revealed by microsatellite markers

BACKGROUND

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.

METHODS

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.

CONCLUSIONS

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.

Identification and Characterization of Microsatellite Markers Derived from the Whole Genome Analysis of Taenia solium

Background

Infections with Taenia solium are the most common cause of adult acquired seizures worldwide, and are the leading cause of epilepsy in developing countries. A better understanding of the genetic diversity of T. solium will improve parasite diagnostics and transmission pathways in endemic areas thereby facilitating the design of future control measures and interventions. Microsatellite markers are useful genome features, which enable strain typing and identification in complex pathogen genomes. Here we describe microsatellite identification and characterization in T. solium, providing information that will assist in global efforts to control this important pathogen.

Methods

For genome sequencing, T. solium cysts and proglottids were collected from Huancayo and Puno in Peru, respectively. Using next generation sequencing (NGS) and de novo assembly, we assembled two draft genomes and one hybrid genome. Microsatellite sequences were identified and 36 of them were selected for further analysis. Twenty T. solium isolates were collected from Tumbes in the northern region, and twenty from Puno in the southern region of Peru. The size-polymorphism of the selected microsatellites was determined with multi-capillary electrophoresis. We analyzed the association between microsatellite polymorphism and the geographic origin of the samples.

Results

The predicted size of the hybrid (proglottid genome combined with cyst genome) T. solium genome was 111 MB with a GC content of 42.54%. A total of 7,979 contigs (>1,000 nt) were obtained. We identified 9,129 microsatellites in the Puno-proglottid genome and 9,936 in the Huancayo-cyst genome, with 5 or more repeats, ranging from mono- to hexa-nucleotide. Seven microsatellites were polymorphic and 29 were monomorphic within the analyzed isolates. T. solium tapeworms were classified into two genetic groups that correlated with the North/South geographic origin of the parasites.

Conclusions/Significance

The availability of draft genomes for T. solium represents a significant step towards the understanding the biology of the parasite. We report here a set of T. solium polymorphic microsatellite markers that appear promising for genetic epidemiology studies.

Taenia solium, the pork tapeworm, is an important pathogen as it is a major cause of acquired epilepsy in developing countries. The parasite was eliminated from most developed countries decades ago due to improvement in sanitary conditions but it remains a common infection across Asia, Africa and Latin America. Identification of genetic variants within T. solium will enable to study the genetic epidemiology, distribution and movement of this parasite within endemic communities, which will ultimately facilitate the design of control strategies to reduce the health and economic burden of disease.

Microsatellites have been used in other parasites to identify genetic variants. In this study, we partially sequenced the genome of T. solium and identified microsatellites widely distributed in the genome using bioinformatics tools. We evaluated the distribution of these microsatellites collected from 20 tapeworms from the north and 20 tapeworms from the south of Peru. We identified seven polymorphic microsatellites, and evaluated their capacity to differentiate genetic variants of T. solium. Interestingly, tapeworms from the North and South of Peru showed different genotypes, suggesting its use as a potential marker to differentiate geographic origin.

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.

Exploring molecular variation in Schistosoma japonicum in China

Schistosomiasis is a neglected tropical disease that affects more than 200 million people worldwide. The main disease-causing agents, Schistosoma japonicum, S. mansoni and S. haematobium, are blood flukes that have complex life cycles involving a snail intermediate host. In Asia, S. japonicum causes hepatointestinal disease (schistosomiasis japonica) and is challenging to control due to a broad distribution of its snail hosts and range of animal reservoir hosts. In China, extensive efforts have been underway to control this parasite, but genetic variability in S. japonicum populations could represent an obstacle to eliminating schistosomiasis japonica. Although a draft genome sequence is available for S. japonicum, there has been no previous study of molecular variation in this parasite on a genome-wide scale. In this study, we conducted the first deep genomic exploration of seven S. japonicum populations from mainland China, constructed phylogenies using mitochondrial and nuclear genomic data sets, and established considerable variation between some of the populations in genes inferred to be linked to key cellular processes and/or pathogen-host interactions. Based on the findings from this study, we propose that verifying intraspecific conservation in vaccine or drug target candidates is an important first step toward developing effective vaccines and chemotherapies against schistosomiasis.

Comparative Phylogenetic Studies on Schistosoma japonicum and Its Snail Intermediate Host Oncomelania hupensis: Origins, Dispersal and Coevolution

BACKGROUND

Schistosoma japonicum causes major public health problems in China and the Philippines; this parasite, which is transmitted by freshwater snails of the species Oncomelania hupensis, causes the disease intestinal schistosomiasis in humans and cattle. Researchers working on Schistosoma in Africa have described the relationship between the parasites and their snail intermediate hosts as coevolved or even as an evolutionary arms race. In the present study this hypothesis of coevolution is evaluated for S. japonicum and O. hupensis. The origins and radiation of the snails and the parasite across China, and the taxonomic validity of the sub-species of O. hupensis, are also assessed.

METHODOLOGY/PRINCIPAL FINDINGS

The findings provide no evidence for coevolution between S. japonicum and O. hupensis, and the phylogeographical analysis suggests a heterochronous radiation of the parasites and snails in response to different palaeogeographical and climatic triggers. The results are consistent with a hypothesis of East to West colonisation of China by Oncomelania with a re-invasion of Japan by O. hupensis from China. The Taiwan population of S. japonicum appears to be recently established in comparison with mainland Chinese populations.

CONCLUSIONS/SIGNIFICANCE

The snail and parasite populations of the western mountain region of China (Yunnan and Sichuan) appear to have been isolated from Southeast Asian populations since the Pleistocene; this has implications for road and rail links being constructed in the region, which will breach biogeographical barriers between China and Southeast Asia. The results also have implications for the spread of S. japonicum. In the absence of coevolution, the parasite may more readily colonise new snail populations to which it is not locally adapted, or even new intermediate host species; this can facilitate its dispersal into new areas. Additional work is required to assess further the risk of spread of S. japonicum.

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.

Wetlands explain most in the genetic divergence pattern of Oncomelania hupensis

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

Reductions in genetic diversity of Schistosoma mansoni populations under chemotherapeutic pressure: the effect of sampling approach and parasite population definition

Detecting potential changes in genetic diversity in schistosome populations following chemotherapy with praziquantel (PZQ) is crucial if we are to fully understand the impact of such chemotherapy with respect to the potential emergence of resistance and/or other evolutionary outcomes of interventions. Doing so by implementing effective, and cost-efficient sampling protocols will help to optimise time and financial resources, particularly relevant to a disease such as schistosomiasis currently reliant on a single available drug. Here we explore the effect on measures of parasite genetic diversity of applying various field sampling approaches, both in terms of the number of (human) hosts sampled and the number of transmission stages (miracidia) sampled per host for a Schistosoma mansoni population in Tanzania pre- and post-treatment with PZQ. In addition, we explore population structuring within and between hosts by comparing the estimates of genetic diversity obtained assuming a 'component population' approach with those using an 'infrapopulation' approach. We found that increasing the number of hosts sampled, rather than the number of miracidia per host, gives more robust estimates of genetic diversity. We also found statistically significant population structuring (using Wright's F-statistics) and significant differences in the measures of genetic diversity depending on the parasite population definition. The relative advantages, disadvantages and, hence, subsequent reliability of these metrics for parasites with complex life-cycles are discussed, both for the specific epidemiological and ecological scenario under study here and for their future application to other areas and schistosome species.

A circular analysis of chronobiology of Schistosoma japonicum cercarial emergence from hilly areas of Anhui, China

About 46 mammal species have been suspected as reservoir hosts for Schistosoma japonicum and therefore the track of the target parasites, in relation to definitive host species, may be of great importance in terms of theoretical and practical implications. The circadian rhythm of cercariae emergence, a genetically controlled behavior for parasites to adapt to their definitive hosts, may seem to be a perfect biological marker for S. japonicum. In this study, a late (or nocturnal) cercarial emergence pattern was observed on the parasites from one hilly region in Anhui of China, where rodents serve as reservoirs, and on the first generation of the parasites. Moreover, by using the circular statistics, the homogeneity of parasites in such trait was also demonstrated. All these provide evidence for the genetically controlled biological trait, which seems essential in the investigation of macro- or micro-dynamics of parasite transmission of interest. This is particularly true in the case of S. japonicum when multiple parasite isolates or strains are more likely to exist.

Low genetic diversity in wide-spread Eurasian liver fluke Opisthorchis felineus suggests special demographic history of this trematode species

Opisthorchis felineus or Siberian liver fluke is a trematode parasite (Opisthorchiidae) that infects the hepato-biliary system of humans and other mammals. Despite its public health significance, this wide-spread Eurasian species is one of the most poorly studied human liver flukes and nothing is known about its population genetic structure and demographic history. In this paper, we attempt to fill this gap for the first time and to explore the genetic diversity in O. felineus populations from Eastern Europe (Ukraine, European part of Russia), Northern Asia (Siberia) and Central Asia (Northern Kazakhstan). Analysis of marker DNA fragments from O. felineus mitochondrial cytochrome c oxidase subunit 1 and 3 (cox1, cox3) and nuclear rDNA internal transcribed spacer 1 (ITS1) sequences revealed that genetic diversity is very low across the large geographic range of this species. Microevolutionary processes in populations of trematodes may well be influenced by their peculiar biology. Nevertheless, we suggest that lack of population genetics structure observed in O. felineus can be primarily explained by the Pleistocene glacial events and subsequent sudden population growth from a very limited group of founders. Rapid range expansion of O. felineus through Asian and European territories after severe bottleneck points to a high dispersal potential of this trematode species.

Introgressive hybridization of Schistosoma haematobium group species in Senegal: species barrier break down between ruminant and human schistosomes

Background

Schistosomes are dioecious parasitic flatworms, which live in the vasculature of their mammalian definitive hosts. They are the causative agent of schistosomiasis, a disease of considerable medical and veterinary importance in tropical and subtropical regions. Schistosomes undergo a sexual reproductive stage within their mammalian host enabling interactions between different species, which may result in hybridization if the species involved are phylogenetically close. In Senegal, three closely related species in the Schistosoma haematobium group are endemic: S. haematobium, which causes urogenital schistosomiasis in humans, and S. bovis and S. curassoni, which cause intestinal schistosomiasis in cows, sheep and goats.

Methodology/Principal Findings

Large-scale multi-loci molecular analysis of parasite samples collected from children and domestic livestock across Senegal revealed that interactions and hybridization were taking place between all three species. Evidence of hybridization between S. haematobium/S. curassoni and S. haematobium/S. bovis was commonly found in children from across Senegal, with 88% of the children surveyed in areas of suspected species overlap excreting hybrid miracidia. No S. haematobium worms or hybrids thereof were found in ruminants, although S. bovis and S. curassoni hybrid worms were found in cows. Complementary experimental mixed species infections in laboratory rodents confirmed that males and females of each species readily pair and produce viable hybrid offspring.

Conclusions/Significance

These data provide indisputable evidence for: the high occurrence of bidirectional hybridization between these Schistosoma species; the first conclusive evidence for the natural hybridisation between S. haematobium and S. curassoni; and demonstrate that the transmission of the different species and their hybrids appears focal. Hybridization between schistosomes has been known to influence the disease epidemiology and enhance phenotypic characteristics affecting transmission, morbidity and drug sensitivity. Therefore, understanding and monitoring such inter-species interactions will be essential for optimizing and evaluating control strategies across such potential hybrid zones.

Schistosome blood flukes are transmitted through water contact and cause a severe debilitating disease in humans and their livestock. Understanding the biology and epidemiology of these parasites is essential to enable the development of better control strategies in endemic areas. Several species of schistosome exist and species barriers are normally maintained by differences in ecology, host specificity, and evolutionary history. However, hybridization between closely related species can occur if parasites infect the same definitive host. Here we report on the introgression between a human (S. haematobium) and two ruminant schistosomes (S. bovis and S. curassoni), the prevalence and distribution of their hybrids and the novel evidence for the presence of S. haematobium/S. curassoni hybrids in Senegalese children. Modern sampling and genotyping techniques have brought to light the extent of these hybridization events which could have been facilitated by the natural progression in farming, population (both human and livestock) movements and expansion, as well as changes in snail distribution, creating areas of sympatric transmission. Hybridization can lead to phenotypic characteristics that can influence disease epidemiology and control success, highlighting the importance of monitoring these evolving hybrid zones.

Population genetic structuring in Opisthorchis viverrini over various spatial scales in Thailand and Lao PDR

Khon Kaen Province in northeast Thailand is known as a hot spot for opisthorchiasis in Southeast Asia. Preliminary allozyme and mitochondrial DNA haplotype data from within one endemic district in this Province (Ban Phai), indicated substantial genetic variability within Opisthorchis viverrini. Here, we used microsatellite DNA analyses to examine the genetic diversity and population structure of O. viverrini from four geographically close localities in Khon Kaen Province. Genotyping based on 12 microsatellite loci yielded a mean number of alleles per locus that ranged from 2.83 to 3.7 with an expected heterozygosity in Hardy–Weinberg equilibrium of 0.44–0.56. Assessment of population structure by pairwise F_ST analysis showed inter-population differentiation (P<0.05) which indicates population substructuring between these localities. Unique alleles were found in three of four localities with the highest number observed per locality being three. Our results highlight the existence of genetic diversity and population substructuring in O. viverrini over a small spatial scale which is similar to that found at a larger scale. This provides the basis for the investigation of the role of parasite genetic diversity and differentiation in transmission dynamics and control of O. viverrini.

Infection with the liver fluke (Opisthorchis viverrini) is a risk factor for cholangiocarcinoma (CCA), which is highly prevalent in Khon Kaen Province, Thailand. Within this Province, there is considerable variation in the epidemiology of opisthorchiasis among districts. Preliminary allozyme and mitochondrial DNA data indicate that genetic variation in O. viverrini occurs even over a small endemic area within the province. Here, we used microsatellite DNA analyses to examine the population genetic structure of O. viverrini from four geographically close localities. Analyses of adult worms based on 12 microsatellite loci revealed varying levels of genetic diversity and population substructuring between the four localities. Worms originating from one locality (Phu Wiang) had significantly higher genetic diversity than the other three localities. Data on the population genetic structure observed in these localities are similar to those found at a larger geographic scale. This provides background data to further investigate the biological and epidemiological significance of these genetic variants.

Retrotransposon-microsatellite amplified polymorphism, an electrophoretic approach for studying genetic variability among Schistosoma japonicum geographical isolates

In the present study, retrotransposon-microsatellite amplified polymorphism (REMAP) was used to examine genetic variability among Schistosoma japonicum isolates from different endemic provinces in mainland China, using S. japonicum from Japan and the Philippines for comparison. Of the 50 primer combinations screened, eight produced highly reproducible REMAP fragments. Using these primers, 190 distinct DNA fragments were generated in total, of which 147 (77.37%) were polymorphic, indicating considerable genetic variation among the 43 S. japonicum isolates examined. The percentage of polymorphic bands (PPB) among S. japonicum isolates from mainland China, Japan, and the Philippines was 77.37%; PPB values of 18.42% and 53.68% were found among isolates from southwestern (SW) China and the lower Yangtze/Zhejiang province in eastern (E) China, respectively. Based on REMAP profiles, unweighted pair-group method with arithmetic averages (UPGMA) dendrogram analysis revealed that all of the S. japonicum samples grouped into three distinct clusters: parasites from mainland China, Japan, and the Philippines were clustered in each individual clade. Within the mainland China cluster, SW China isolates (from Sichuan and Yunnan provinces) grouped together, whereas worms from E China (Zhejiang, Anhui, Jiangxi, Jiangsu, Hunan, and Hubei provinces) grouped together. These results demonstrated that the REMAP marker system provides a reliable electrophoretic technique for studying genetic diversity and population structures of S. japonicum isolates from mainland China, and could be applied to other pathogens of human and animal health significance.

Differences in genomic architecture between two distinct geographical strains of the blood fluke Schistosoma japonicum reveal potential phenotype basis

The Chinese (SjC) and Philippine (SjP) strains of the blood fluke Schistosoma japonicum have been shown to present clearly different phenotypes in fecundity, pathology, drug sensitivity and immunology. We used microarray based comparative genomic hybridisation (aCGH) to investigate structural differences in the genomes of the two strains and identified seven distinct regions of the S. japonicum genome that present differential aCGH representing either deletion or duplication regions in SjP. Within these regions, genes predicted to be associated with the recognised phenotypic differences were identified and that may provide new insights into the biology and evolution of the two strains, with implications for the epidemiology and control of schistosomiasis japonica in China and the Philippines.

A research agenda for helminth diseases of humans: basic research and enabling technologies to support control and elimination of helminthiases

Successful and sustainable intervention against human helminthiases depends on optimal utilisation of available control measures and development of new tools and strategies, as well as an understanding of the evolutionary implications of prolonged intervention on parasite populations and those of their hosts and vectors. This will depend largely on updated knowledge of relevant and fundamental parasite biology. There is a need, therefore, to exploit and apply new knowledge and techniques in order to make significant and novel gains in combating helminthiases and supporting the sustainability of current and successful mass drug administration (MDA) programmes. Among the fields of basic research that are likely to yield improved control tools, the Disease Reference Group on Helminth Infections (DRG4) has identified four broad areas that stand out as central to the development of the next generation of helminth control measures: 1) parasite genetics, genomics, and functional genomics; 2) parasite immunology; 3) (vertebrate) host–parasite interactions and immunopathology; and 4) (invertebrate) host–parasite interactions and transmission biology. The DRG4 was established in 2009 by the Special Programme for Research and Training in Tropical Diseases (TDR). The Group was given the mandate to undertake a comprehensive review of recent advances in helminthiases research in order to identify notable gaps and highlight priority areas. This paper summarises recent advances and discusses challenges in the investigation of the fundamental biology of those helminth parasites under the DRG4 Group's remit according to the identified priorities, and presents a research and development agenda for basic parasite research and enabling technologies that will help support control and elimination efforts against human helminthiases.

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.

Identification of differentially expressed genes in Oncomelania hupensis chronically infected with Schistosoma japonicum

Oncomelania hupensis is the unique intermediate host of Schistosoma japonicum. The schistosome-snail interaction is biomedically important. To identify differentially expressed transcripts in O. hupensis chronically infected with S. japonicum, suppression subtractive hybridization (SSH) was used to construct a cDNA library in each direction for transcripts that are more abundantly enriched in head-foot part of the infected O. hupensis and for those that are more abundantly enriched in the uninfected, as head-foot part contains hemocytes and hemolymph which are associated with the snail internal defense system. After differential screening, 39 transcripts were identified, including nine and 30 transcripts enriched in infected and uninfected snails, respectively. Some of the transcripts have similar homology to available sequences in current databases, including transposase, caveolin-like protein, pancreatic trypsin inhibitor-like protein, prosaposin, glutathione s-transferase (GST), and several hypothetical proteins, while most of the transcripts do not match with any sequences in available databases. The identified transcripts were involved functionally in cell growth, metabolism, signal transduction, and immune responses. Two forward library transcripts and 11 reverse library transcripts were selected for real-time PCR, and 10 of them were confirmed to be consistent with the SSH results. It is intriguing to continue functional studies for some genes such as pancreatic trypsin inhibitor; a hypothetical protein (HS576367) related to calcium ion binding; GST; and several unknown proteins (HS576353 and HS576355). These identified differentially expressed genes may be key targets for understanding the molecular mechanism of co-existence during which the snail is unable to rid itself of the schistosome in chronic infection stage.

Genetic diversity of Schistosoma japonicum miracidia from individual rodent hosts

Schistosoma japonicum is an important parasite in terms of clinical, veterinary and socio-economic impacts, and rodents, a long neglected reservoir for the parasite, have recently been found to act as reservoir hosts in some endemic areas of China. Any difference in the host's biological characteristics and/or associated living habitats among rodents may result in different environments for parasites, possibly resulting in a specific population structure of parasites within hosts. Therefore knowledge of the genetic structure of parasites within individual rodents could improve our understanding of transmission dynamics and hence our ability to develop effective control strategies. In this study, we aimed to describe a host-specific structure for S. japonicum and its potential influencing factors. The results showed a significant genetic differentiation among hosts. Two factors, including sampling seasons and the number of miracidia genotyped per host, showed an effect on the genetic diversity of an infrapopulation through a univariable analysis but not a multivariable analysis. A possible scenario of clustered infection foci and the fact of multiple definitive host species, the latter of which is unique to S. japonicum compared with other schistosomes, were proposed to explain the observed results and practical implications for control strategies are recommended.

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.

Neuroschistosomiasis

Schistosomiasis (bilharzia) is a neglected tropical disease caused by digenetic trematode platyhelminths of the genus Schistosoma. Neuroschistosomiasis is one of the most severe clinical outcomes associated with schistosome infection. Neurological complications early during the course of infection are thought to occur through in situ egg deposition following aberrant migration of adult worms to the brain or spinal cord. The presence of eggs in the CNS induces a cell-mediated Th2-driven periovular granulomatous reaction. The mass effect of thousands of eggs and the large granulomas concentrated within the brain or spinal cord explain the signs and symptoms of increased intracranial pressure, myelopathy, radiculopathy and subsequent clinical sequelae. Myelopathy (acute transverse myelitis and subacute myeloradiculopathy) of the lumbosacral region is the most common neurological manifestation of S. mansoni or S. haematobium infection, whereas acute encephalitis of the cortex, subcortical white matter, basal ganglia or internal capsule is typical of S. japonicum infection. Cerebral complications include encephalopathy with headache, visual impairment, delirium, seizures, motor deficits and ataxia, whereas spinal symptoms include lumbar pain, lower limb radicular pain, muscle weakness, sensory loss and bladder dysfunction. The finding of eggs in the stool or a positive serology, provides supportive but not direct evidence of neuroschistosomiasis. A definitive diagnosis can only be made with histopathological study showing Schistosoma eggs and granulomas. Schistosomicidal drugs (notably praziquantel), steroids and surgery are currently used for the treatment of neuroschistosomiasis. During the 'acute phase' of the disease, neuroschistosomiasis is treated with corticosteroids which are augmented with a course of praziquantel once female worm ovipositioning commences. Surgery should be reserved for special cases such as in those with evidence of medullary compression and in those who deteriorate despite clinical management.

IRAP: An efficient retrotransposon-based electrophoretic technique for studying genetic variability among geographical isolates of Schistosoma japonicum

In the present study, a inter-retrotransposon-amplified polymorphism (IRAP) technique, based on retrotransposons, was used to examine genetic variability among Schistosoma japonicum isolates from different provinces in mainland China. Of the 15 primers screened, 5 produced highly reproducible IRAP patterns. Using these primers, 54 discernible DNA fragments were generated with 40 (74.07%) being polymorphic, indicating considerable genetic variation among the examined S. japonicum isolates. The primer LTR-11 was found to be able to differentiate male and female parasites, producing one constant specific band for female S. japonicum isolates. The percentages of polymorphic bands (PPB) among all parasites, among isolates from mountainous provinces and among those from the lake/marshland areas were 74.07, 48.15, and 66.67%, respectively. UPGMA analysis revealed that the IRAP profiles could group S. japonicum isolates in mainland China into two clades (mountainous and lake/marshland types), and samples from the same geographical origins clustered together. These results demonstrated that the IRAP technique is suitable for studying genetic diversity and population structures, and also provides an effective technique for studying sex differentiation of S. japonicum.

Electrophoretic detection of genetic variability among Schistosoma japonicum isolates by sequence-related amplified polymorphism

In the present study, sequence-related amplification polymorphism (SRAP) was utilized to study the genetic variability among Schistosoma japonicum isolates from different provinces in China, using Schistosoma mansoni from Puerto Rico for comparison. Five out of ten tested SRAP primer combinations displayed significant polymorphisms among S. japonicum isolates from China, namely ME2/EM1, ME4/EM1, ME4/EM6, ME5/EM4 and ME5/EM5. Analysis of the 61 S. japonicum samples from China with five SRAP primer combinations identified a total of 83 reproducible polymorphic fragments. The number of fragments using each primer combination ranged from 14 to 19, with an average of 16 polymorphic bands per primer pair, and the size of fragment ranged approximately from 100 to 1000 bp. Representative-specific SRAP fragments were excised from the gels, and confirmed by PCR amplification of genomic DNA using primers designed and based on the sequences of these SRAP fragments. Based on SRAP profiles, unweighted pair-group method with arithmetic averages (UPGMA) dendrogram was constructed. UPGMA clustering algorithm categorized S. japonicum isolates from China into nine clades and two lineages (representing the mountainous and lake/marshland regions). These results indicate the usefulness of the SRAP technique for revealing genetic variability among S. japonicum isolates from China, and the SRAP technique should be applicable to other living organisms.

Transmission dynamics of two strains of Schistosoma mansoni utilizing novel intermediate and definitive hosts

The intimate host-parasite relationship mandates adaptation to the genetic and phenotypic variability of their counterparts. Here, inbred and outcrossed strains of Schistosoma mansoni were challenged with "local" and "novel" intermediate and definitive hosts to examine effects of genetic variability and novelty on infection success and dynamics. Genetically distinct lines of Biomphalaria glabrata intermediate hosts exposed to inbred and outcrossed S. mansoni larvae were assessed for differences in both snail and parasite life-history parameters. Cercariae from each parasite-snail treatment were used to infect "local" and "novel" Mus musculus definitive hosts to assess parasite infectivity and fitness. Outcrossed parasites significantly reduced snail growth, were more productive, and induced greater host mortality than inbred parasites. Mouse strain did not influence parasite infectivity or reproduction, but parasite and snail host genetic background did, affecting both sex-specific infectivity and parasite productivity. Overall, genetic background of S. mansoni and its intermediate snail host altered life history traits and transmission dynamics of the parasite throughout its life cycle.

Polymorphic microsatellites in the human bloodfluke, Schistosoma japonicum, identified using a genomic resource

Re-emergence of schistosomiasis in regions of China where control programs have ceased requires development of molecular-genetic tools to track gene flow and assess genetic diversity of Schistosoma populations. We identified many microsatellite loci in the draft genome of Schistosoma japonicum using defined search criteria and selected a subset for further analysis. From an initial panel of 50 loci, 20 new microsatellites were selected for eventual optimization and application to a panel of worms from endemic areas. All but one of the selected microsatellites contain simple tri-nucleotide repeats. Moderate to high levels of polymorphism were detected. Numbers of alleles ranged from 6 to 14 and observed heterozygosity was always >0.6. The loci reported here will facilitate high resolution population-genetic studies on schistosomes in re-emergent foci.

Schistosomiasis in the People's Republic of China: the era of the Three Gorges Dam

The potential impact of the Three Gorges Dam (TGD) on schistosomiasis transmission in China has invoked considerable global concern. The TGD will result in changes in the water level and silt deposition downstream, favoring the reproduction of Oncomelania snails. Combined with blockages of the Yangtze River's tributaries, these changes will increase the schistosomiasis transmission season within the marshlands along the middle and lower reaches of the Yangtze River. The changing schistosome transmission dynamics necessitate a comprehensive strategy to control schistosomiasis. This review discusses aspects of the epidemiology and transmission of Schistosoma japonicum in China and considers the pathology, clinical outcomes, diagnosis, treatment, immunobiology, and genetics of schistosomiasis japonica together with an overview of current progress in vaccine development, all of which will have an impact on future control efforts. The use of synchronous praziquantel (PZQ) chemotherapy for humans and domestic animals is only temporarily effective, as schistosome reinfection occurs rapidly. Drug delivery requires a substantial infrastructure to regularly cover all parts of an area of endemicity. This makes chemotherapy expensive and, as compliance is often low, a less than satisfactory control option. There is increasing disquiet about the possibility that PZQ-resistant schistosomes will develop. Consequently, as mathematical modeling predicts, vaccine strategies represent an essential component in the future control of schistosomiasis in China. With the inclusion of focal mollusciciding, improvements in sanitation, and health education into the control scenario, China's target of reducing the level of schistosome infection to less than 1% by 2015 may be achievable.

Molecular epidemiology and phylogeography of Schistosoma mansoni around Lake Victoria

Intestinal schistosomiasis continues to be a major public health problem in sub-Saharan Africa, and is endemic in communities around Lake Victoria. Interest is growing in the molecular evolution and population genetic structure of Schistosoma mansoni and we describe a detailed analysis of the molecular epidemiology and phylogeography of S. mansoni from Lake Victoria. In total, 388 cytochrome oxidase 1 (COI) sequences were obtained from 25 sites along the Ugandan, Tanzanian and Kenyan shorelines of Lake Victoria, and 122 unique barcodes were identified; 9 corresponded to previously discovered barcodes from Lakes Victoria and Albert. A subset of the data, composed of COI sequences from miracidia from 10 individual children, was used for population genetics analyses; these results were corroborated by microsatellite analysis of 4 isolates of lab-passaged adult worms. Overall, 12 barcodes were found to be shared across all 3 countries, whereas the majority occurred singly and were locally restricted. The population genetics analyses were in agreement in revealing high diversity at the level of the human host and negligible population structuring by location. The lack of correlation between genetic distance and geographical distance in these data may be attributed to the confounding influence of high intra-individual diversity as well as human migration between communities.

Combined mitochondrial 16S and 12S rDNA sequences: an effective genetic marker for inter-species phylogenetic analysis of zoonotic trematodes

The present study studied the genetic variation among Schistosoma japonicum isolates from different endemic regions in mainland China and examined the phylogenetic relationships of zoonotic trematodes using the combined mitochondrial 16S and 12S ribosomal DNA sequences. The fragments of 16S and 12S rDNA were amplified from 22 S. japonicum isolates, and sequenced, and the relevant sequences of other nine trematode species belonging to six genera in four families were downloaded from GenBank, and their phylogenetic relationships were re-constructed by unweighted pair-group method with arithmetic averages analyses using the combined 16S and 12S rDNA sequences, with Trichinella spiralis as outgroup. The results showed that the partial sequences of mitochondrial 16S and 12S rDNA of S. japonicum were 757 and 797 bp, respectively, and they were quite conserved among the S. japonicum isolates. Phylogenetic analysis revealed that the combined 16S and 12S rDNA sequences were not able to distinguish S. japonicum isolates in mountainous areas from those in lake/marshland areas in mainland China. However, the combined sequences could distinguish different species of zoonotic trematodes. Therefore, the combined mitochondrial 16S and 12S rDNA sequences provide an effective molecular marker for the inter-species phylogenetic analysis and differential identification of zoonotic trematodes.

Applying evolutionary genetics to schistosome epidemiology

We review how molecular markers and evolutionary analysis have been applied to the study of schistosome parasites, important pathogens that infect over 200 million people worldwide. Topics reviewed include phylogenetics and biogeography, hybridization, infection within snails, mating systems, and genetic structure. Some interesting generalizations include that schistosome species hybridize frequently and have switched definitive hosts repeatedly in evolutionary time. We show that molecular markers can be used to infer epidemiologically relevant processes such as spatial variation in transmission, or to reveal complex patterns of mate choice. Analysis of genetic structure data shows that transmission foci can be structured by watershed boundaries, habitat types, and host species. We also discuss sampling and analytical problems that arise when using larvae to estimate genetic parameters of adult schistosome populations. Finally, we review pitfalls in methodologies such as genotyping very small individuals, statistical methods for identifying clonemates or for identifying sibling groups, and estimating allele frequencies from pooled egg samples.

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

Background

Oncomelania hupensis is the unique intermediate host of Schistosoma japonicum, which causes schistosomiasis endemic in the Far East, and especially in mainland China. O. hupensis largely determines the parasite's geographical range. How O. hupensis's genetic diversity is distributed geographically in mainland China has never been well examined with DNA sequence data.

Methodology/Principal Findings

In this study we investigate the genetic variation among O. hupensis from different geographical origins using the combined complete internal transcribed spacer 1 (ITS1) and ITS2 regions of nuclear ribosomal DNA. 165 O. hupensis isolates were obtained in 29 localities from 7 provinces across mainland China: lake/marshland and hill regions in Anhui, Hubei, Hunan, Jiangxi and Jiangsu provinces, located along the middle and lower reaches of Yangtze River, and mountainous regions in Sichuan and Yunnan provinces. Phylogenetic and haplotype network analyses showed distinct genetic diversity and no shared haplotypes between populations from lake/marshland regions of the middle and lower reaches of the Yangtze River and populations from mountainous regions of Sichuan and Yunnan provinces. The genetic distance between these two groups is up to 0.81 based on Fst, and branch time was estimated as 2–6 Ma. As revealed in the phylogenetic tree, snails from Sichuan and Yunnan provinces were also clustered separately. Geographical separation appears to be an important factor accounting for the diversification of the two groups of O. hupensis in mainland China, and probably for the separate clades between snails from Sichuan and Yunnan provinces. In lake/marshland and hill regions along the middle and lower reaches of the Yangtze River, three clades were identified in the phylogenetic tree, but without any obvious clustering of snails from different provinces.

Conclusions

O. hupensis in mainland China may have considerable genetic diversity, and a more complex population structure than expected. It will be of significant importance to consider the genetic diversity of O. hupensis when assessing co-evolutionary interactions with S. japonicum.

The intermediate host of Schistosoma japonicum in Asia is the snail Oncomelania hupensis, which can be separated phenotypically into ribbed- and smooth-shelled morphotypes. In China, the typical morphotype is ribbed-shelled, with its distribution restricted to mainland China. Smooth-shelled snails with varix are also distributed in China, which are considered to belong to the same subspecies as the ribbed-shelled snails. In this study we investigate the genetic variation among O. hupensis from different geographical origins using combined complete ITS1 and ITS2 regions of nuclear ribosomal DNA. Snails including ribbed-shelled and smooth-shelled (but with varix on the shell) from the lake/marshland region of the middle and lower reaches of the Yangtze River, and smooth-shelled snails from mountainous regions of Sichuan and Yunnan provinces, were genetically distinct with no shared haplotypes detected. Furtheremore, the snails from Sichuan and Yunnan provinces were clustered in separate clades in the phylogenetic tree, and three clades were observed for snails from the middle and lower reaches of the Yangtze River. The population diversity of O. hupensis in China is thus considered large, and evolutionary relationships in the host-parasite system of O. hupensis-S. japonicum may be of interest for further research.