Diversification of Schistosoma japonicum in Mainland China revealed by mitochondrial DNA

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.

Genetic difference between two Schistosoma japonicum isolates with contrasting cercarial shedding patterns revealed by whole genome sequencing

Schistosoma japonicum is one of the major infectious agents of human schistosomiasis, mainly endemic in China and the Philippines. We have previously reported the finding of two schistosome isolates, each with a different cercarial emergence pattern adapted to their different hosts. However, there are currently no whole-genome sequencing studies to investigate the underlining genetics of the adaptive traits. We sampled schistosomes in 2013 and 2020 from a hilly area Shitai (ST) and a marshland area Hexian (HX) of Anhui, China. Ten to 15 male or female adult worms from each site/year were sent for whole genome sequencing. Genetics were analyzed, and selection signals along genomes were detected. Gene enrichment analysis was performed for the genome regions under selection. The results revealed considerable genetic differentiation between the two isolates. The genome "windows" affected by natural selection were fewer in ST (64 windows containing 78 genes) than in HX (318 windows containing 276 genes). Twelve significantly enriched genes were identified in ST, but none in HX. These genes were mainly related to specific DNA binding and intercellular signaling transduction. Some functional region changes identified along the genome of the hilly schistosome may be related to its unique late afternoon cercarial emergence.

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.

Duplex real-time PCR for sexing Schistosoma japonicum cercariae based on W chromosome-specific genes and its applications

As a unique feature among otherwise hermaphroditic trematodes, Schistosoma species are gonochoric parasites whose sex is genetically determined (ZZ for males and ZW for females). However, schistosome larvae are morphologically identical, and sex can only be discriminated by molecular methods. Here, we integrated published Schistosoma. japonicum transcriptome and genome data to identify W chromosome-specific genes as sex biomarkers. Three W chromosome-specific genes of S. japonicum were identified as sex biomarkers from a panel of 12 genes expressed only in females. An efficient duplex real-time PCR (qPCR) method for sexing cercariae was developed which could identify the sex of cercariae within 2 h without DNA extraction. Moreover, this method can be used to identify not only single-sex but also mixed-sex schistosome-infected snails. We observed a nearly equal proportion of single-male, single-female, and mixed-sex schistosome infections in artificially infected snails. Sex-known schistosome-infected snail models can be efficiently constructed with the aid of duplex qPCR. A field study revealed that single-sex schistosome infections were predominant among naturally infected snails. Finally, a schistosomiasis mouse model based on sex-known cercariae infection was shown to be more reliable than a model based on sex-unknown cercariae infection. The developed duplex qPCR method for sexing S. japonicum cercariae can be widely used for schistosomiasis modeling, genetic experiments, and field-based molecular epidemiological studies.

Schistosoma japonicum is a major causative agent of human schistosomiasis. Unlike other parasitic worms, S. japonicum females are determined by the heterogametic sex chromosome (ZW) and males by the homogametic sex chromosome (ZZ). The life cycle of S. japonicum includes the egg, miracidium, mother sporocyst, daughter sporocyst, cercaria, schistosomulum, and adult stages. The sex of adult male and female worms can be morphologically distinguished, whereas the sex of larvae, such as cercariae, can only be discriminated by molecular methods. In this study, we established an efficient duplex real-time PCR method for sexing S. japonicum cercariae based on newly identified W chromosome-specific genes. The established duplex real-time PCR method will facilitate construction of sex-controlled schistosome-infected intermediate host or definitive host models for schistosome-host interplays and schistosomiasis studies. This method is also a powerful tool for investigating the epidemiology of single-sex and mixed-sex schistosome-infected snails in the field.

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.

Comparative analyses of complete ribosomal DNA sequences of Clonorchis sinensis and Metorchis orientalis: IGS sequences may provide a novel genetic marker for intraspecific variation

Both Clonorchis sinensis and Metorchis orientalis are the fish-borne zoonotic trematodes, and have a wide distribution of southeastern Asia, especially in China. Due to the similar morphology, life cycle, and parasitic positions are difficult to differentiate between both metacercariae. In the present study, the complete rDNA sequences of five C. sinensis and five M. orientalis were obtained and compared for the first time. And the IGS rDNA sequences were tested as a genetic marker. The results showed complete rDNA lengths of C. sinensis were range from 8049 bp to 8391 bp, including 1991 bp, 1116 bp, 3854 bp, and 1088-1430 bp belonging to 18S, ITS, 28S and IGS, respectively. And the complete rDNA lengths of M. orientalis were range from 7881 bp to 9355 bp, including 1991 bp, 1077 bp, 3856 bp, and 957-2431 bp belonging to 18S, ITS, 28S and IGS, respectively. Comparative analyses reveal length difference main in IGS, which has higher intraspecific and interspecific variations than other ribosomal regions. Forty four repeat (forward and inverted) sequences were found in the complete rDNAs of C. sinensis and M. orientalis. The phylogenetic analyses showed that the sequences of ITS1, ITS2, 18S and 28S could be used as different level genetic markers. In IGS phylogenetic tree, Opisthorchiidae, Paramphistomidae, Dicrocoeliidae, and Schistosomatidae formed monophyletic groups, and the same length sequences were clustered together in the same species. These findings of the present study provide the new molecular data for studying the complete rDNA of C. sinensis and M. orientalis, and indicate IGS sequences may used as a novel genetic marker for studying intraspecific variation in trematodes.

High throughput data analyses of the immune characteristics of Microtus fortis infected with Schistosoma japonicum

Microtus fortis exhibits natural resistance against Schistosoma japonicum, and the parasite cannot grow and develop in M. fortis. Extensive research has been carried out, however, the associated mechanism remains unclear. In the present study, we analysed the combined data obtained from a cytokine chip assay, transcriptome, and metabolome. The cytokine profile from C57BL/6 and M. fortis mice was assessed before and after infection. Several cytokines increased during the second and third week post-infection. Some transcripts related to cytokine genes and associated proteins were also highly expressed (i.e., Hgf, C3, and Lbp). The liver metabolism of M. fortis following infection with S. japonicum was assessed. We identified 25 different metabolites between the uninfected and infected M. fortis, and 22 different metabolites between infected M. fortis and C57BL/6 mice. The metabolomic pathways of these differential metabolites were then analysed with MetPA, revealing that they were involved in histidine metabolism, valine, leucine, and isoleucine biosyntheses, and lysine degradation. Thus, the elevated expression of these metabolites and pathways may promote the phagocytic function of the neutrophils and natural killer cell activity following TLR activation. These results provide novel insight into the resistance mechanism of M. fortis against S. japonicum.

Signatures of mito-nuclear discordance in Schistosoma turkestanicum indicate a complex evolutionary history of emergence in Europe

High levels of molecular diversity were identified in mitochondrial cytochrome c oxidase (cox1) gene sequences of Schistosoma turkestanicum from Hungary. These cox1 sequences were all specific to Hungary which contrasted with the low levels of diversity seen in the nuclear internal transcribed spacer region (ITS) sequences, the majority of which were shared between China and Iran isolates. Measures of within and between host molecular variation within S. turkestanicum showed there to be substantial differences in molecular diversity, with cox1 being significantly more diverse than the ITS. Measures of haplotype frequencies revealed that each host contained its own subpopulation of genetically unique parasites with significant levels of differentiation. Pairwise mismatch analysis of cox1 sequences indicated S. turkestanicum populations to have a bimodal pairwise difference distribution and to be stable unlike the ITS sequences, which appeared to have undergone a recent population expansion event. Positive selection was also detected in the cox1 sequences, and biochemical modelling of the resulting protein illustrated significant mutational events causing an alteration to the isoelectric point of the cox1 protein, potentially altering metabolism. The evolutionary signature from the cox1 indicates local adaptation and long establishment of S. turkestanicum in Hungary with continual introgression of nuclear genes from Asian isolates. These processes have led to the occurrence of mito-nuclear discordance in a schistosome population.

Characterization of the complete nuclear ribosomal DNA sequences of Eurytrema pancreaticum

Eurytrema pancreaticum is one of the most common trematodes of cattle and sheep, and also infects humans occasionally, causing great economic losses and medical costs. In this study, the sequences of the complete nuclear ribosomal DNA (rDNA) repeat units of five E. pancreaticum individuals were determined for the first time. They were 8306–8310 bp in length, including the small subunit (18S) rDNA, internal transcribed spacer 1 (ITS1), 5.8S rDNA, internal transcribed spacer 2 (ITS2), large subunit (28S) rDNA and intergenic spacer (IGS). There were no length variations in any of the investigated 18S (1996 bp), ITS1 (1103 bp), 5.8S (160 bp), ITS2 (231 bp) or 28S (3669 bp) rDNA sequences, whereas the IGS rDNA sequences of E. pancreaticum had a 4-bp length variation, ranging from 1147 to 1151 bp. The intraspecific variations within E. pancreaticum were 0–0.2% for 18S rDNA, 0–0.5% for ITS1, 0% for 5.8S rDNA and ITS2, 0–0.2% for 28S rDNA and 2.9–20.2% for IGS. There were nine types of repeat sequences in ITS1, two types in 28S rDNA, but none in IGS. A phylogenetic analysis based on the 18S rDNA sequences classified E. pancreaticum in the family Dicrocoeliidae of Plagiorchiata, closely related to the suborder Opisthorchiata. These results provide useful information for the further study of Dicrocoeliidae trematodes.

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 dynamic changes of CD3e-CD11c+ dendritic cells in spleens and bone marrow of mice infected with Schistosoma japonicum

Schistosoma japonicum as a pathogeny requires dendritic cells to activate immune response. So, the research is to study the dynamic changes of CD3e^-CD11c⁺ dendritic cells in mice infected with S. japonicum. Zero, 7, 28, 35, and 63 days were selected to study the variation of dendritic cells, and the proportions of CD3e^-CD11c⁺ dendritic cells and CD86⁺ mature dendritic cells in spleens and bone marrow were tested by flow cytometry. As a result, the variation trends of dendritic cells in spleen and bone marrow are similar as follows: the proportions of CD3e^-CD11c⁺ dendritic cells increased first and then decreased from day 35, but the percentages of CD86⁺ mature dendritic cells decreased from day 28 and increased in day 63. In vitro, cultured dendritic cells treated with SEA and SAWA were tested by flow cytometry, the variation trends of CD86 on dendritic cells are consistent with the results in days 28 and 63. Besides CD86, the expression of MHC-II also hints immune regulation. In conclusion, it is speculated that dendritic cells play a role of immune regulation through MHC-II and CD86 in S. japonicum infection. Immune regulation of dendritic cells is not only in favor of the survival of host and parasite but also can be used in the therapy for immune diseases.

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.

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.

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.

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.

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

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

Molecular epidemiology and genetic diversity of Entamoeba species in a chelonian collection

Veterinary medicine has focused recently on reptiles, due to the existence of captive collections in zoos and an increase in the acquisition of reptiles as pets. The protozoan parasite, Entamoeba can cause amoebiasis in various animal species and humans. Although amoebiasis disease is remarkably rare in most species of chelonians and crocodiles, these species may serve as Entamoeba species carriers that transmit parasites to susceptible reptile species, such as snakes and lizards, which can become sick and die. In this study, we identified the Entamoeba species in a population of healthy (disease-free) chelonians, and evaluated their diversity through the amplification and sequencing of a small subunit rDNA region. Using this procedure, three Entamoeba species were identified: Entamoeba invadens in 4.76 % of chelonians, Entamoeba moshkovskii in 3.96 % and Entamoeba terrapinae in 50 %. We did not detect mixed Entamoeba infections. Comparative analysis of the amplified region allowed us to determine the intra-species variations. The E. invadens and E. moshkovskii strains isolated in this study did not exhibit marked differences with respect to the sequences reported in GenBank. The analysis of the E. terrapinae isolates revealed three different subgroups (A, B and C). Although subgroups A and C were very similar, subgroup B showed a relatively marked difference with respect to subgroups A and C (Fst = 0.984 and Fst = 1.000, respectively; 10-14 % nucleotide variation, as determined by blast) and with respect to the sequences reported in GenBank. These results suggested that E. terrapinae subgroup B may be either in a process of speciation or belong to a different lineage. However, additional research is necessary to support this statement conclusively.

A deep analysis of the small non-coding RNA population in Schistosoma japonicum eggs

Background

Schistosoma japonicum is a parasitic flatworm that causes zoonotic schistosomiasis. The typical outcome of schistosomiasis is hepatic granuloma and fibrosis, which is primarily induced by soluble egg-derived antigens. Although schistosomal eggs represent an important pathogenic stage to the host, the biology of this critical stage is largely unknown. We previously investigated the expression profiles of sncRNAs during different developmental stages of this parasite. However, using small RNA extracted from egg-deposited liver tissues generated limited information about sncRNAs in eggs. Here, we characterized the complete small RNAome in this stage of the parasite after optimization of RNA purification.

Methodology and Principal Findings

A library, SjE, was constructed with the small RNA extracted from S. japonicum eggs and subjected to high-throughput sequencing. The data were depicted by comprehensive bioinformatic analysis to explore the expression features of sncRNAs in the egg stage. MicroRNAs accounted for about one quarter of the total small RNA population in this stage, with a strongly biased expression pattern of certain miRNA family members. Sja-miR-71, sja-miR-71-5p, and sja-miR-36-3p were suggested to play important roles in embryo development. A panel of transfer RNA fragments (tRFs) precisely processed from the 5′ end of mature tRNAs was identified for the first time, which represented a strong egg stage-biased expression. The tRNA-Ala derived small RNAs were the most highly expressed Sj-tRFs in eggs. Further, the expression of siRNAs from 29 types of well-defined transposable elements (TEs) was observed to be relatively stable among different developmental stages.

Conclusions and Significance

In this study, we characterized the sncRNA profile in the egg stage of S. japonicum. Featured expression of sncRNAs, especially the tRNA-derived small RNAs, was identified, which was further compared with that of other developmental stages. These novel findings would facilitate a deeper understanding of the biology of schistosomal parasites.

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.