Genetic variability among Schistosoma japonicum isolates from different endemic regions in China revealed by sequences of three mitochondrial DNA genes

Study of the population genetic structure of Opisthorchis-like eggs in northern Thailand using mitochondrial genes

BACKGROUND

Opisthorchis-like eggs are a public health problem in northern and northeastern Thailand. However, the genetic epidemiology and structure of these parasites in northern Thailand are unknown. Thus, this study investigated their population genetic structure using cytochrome c oxidase subunit 1 (cox1) and NADH dehydrogenase subunit 1 (nad1) nucleotide sequences.

METHODOLOGY/PRINCIPAL FINDINGS

A study was conducted in the hill tribe regions of Chiang Mai Province, northern Thailand. Internal transcribed spacer 2 polymerase chain reaction and restriction fragment length polymorphism were used to distinguish 205 positive feces samples for Opisthorchis-like eggs. The results showed that the prevalence of O. viverrini and Haplorchis taichui was 10.5% and 38.2%, respectively, and the co-infection rate was 37.2%. To determine the genetic structure of O. viverrini and H. taichui using cox1 and nad1 genes, genetic analysis was performed using 30 randomly chosen fecal samples for Opisthorchis-like eggs. Pairwise FST analysis indicated that O. viverrini and H. taichui displayed nonsignificant genetic differentiation within Chiang Mai Province and between interpopulations from different geographic areas. Moreover, within the intrapopulation in Chiang Mai Province, cox1 presented higher gene flow than nad1 in O. viverrini, while nad1 demonstrated higher gene flow than cox1 in H. taichui. The neutrality tests based on Fu's Fs indicated population expansion and selective sweep from bottleneck or hitchhiking in O. viverrini and H. taichui populations, supported by haplotype network patterns. Phylogenetic tree analysis based on cox1 and nad1 revealed the monophyly of O. viverrini and H. taichui and genetic relationships with other isolates collected from Thailand, Lao People's Democratic Republic (PDR), and Vietnam.

CONCLUSIONS/SIGNIFICANCE

This study investigated the molecular discrimination and genetic structure of Opisthorchis-like eggs in northern Thailand. The genetic information derived from this study could be associated with the background, molecular epidemiology, and disease severity of these parasites.

Application of a Genus-Specific LAMP Assay for Schistosome Species to Detect Schistosoma haematobium x Schistosoma bovis Hybrids

Schistosomiasis is a disease of great medical and veterinary importance in tropical and subtropical regions caused by different species of parasitic flatworms of the genus Schistosoma. The emergence of natural hybrids of schistosomes indicate the risk of possible infection to humans and their zoonotic potential, specifically for Schistosoma haematobium and S. bovis. Hybrid schistosomes have the potential to replace existing species, generate new resistances, pathologies and extending host ranges. Hybrids may also confuse the serological, molecular and parasitological diagnosis. Currently, LAMP technology based on detection of nucleic acids is used for detection of many agents, including schistosomes. Here, we evaluate our previously developed species-specific LAMP assays for S. haematobium, S. mansoni, S. bovis and also the genus-specific LAMP for the simultaneous detection of several Schistosoma species against both DNA from pure and, for the first time, S. haematobium x S. bovis hybrids. Proper operation was evaluated with DNA from hybrid schistosomes and with human urine samples artificially contaminated with parasites' DNA. LAMP was performed with and without prior DNA extraction. The genus-specific LAMP properly amplified pure Schistosoma species and different S. haematobium-S. bovis hybrids with different sensitivity. The Schistosoma spp.-LAMP method is potentially adaptable for field diagnosis and disease surveillance in schistosomiasis endemic areas where human infections by schistosome hybrids are increasingly common.

The Complete Mitochondrial Genome of the Caecal Fluke of Poultry, Postharmostomum commutatum, as the First Representative from the Superfamily Brachylaimoidea

Postharmostomum commutatum (Platyhelminthes: Brachylaimoidea), a parasite of the caeca of poultry, has been frequently reported from many countries and regions, including China. However, the molecular epidemiology, population genetics and phylogenetics of this parasite are poorly understood. In the present study, we determined and characterized the complete mitochondrial (mt) genome of P. commutatum, as the first representative from the superfamily Brachylaimoidea. The mt genome of P. commutatum is a circular DNA molecule of 13,799 bp in size and encodes the complete set of 36 genes (12 protein-coding genes, 22 transfer RNA genes, two ribosomal RNA genes) as well as a typical control region. The mt genome of P. commutatum presents a clear bias in nucleotide composition with a negative AT-skew on average (-0.306) and a positive GC-skew on average (0.466). Phylogenetic analyses showed that P. commutatum (superfamily Brachylaimoidea) and other ten members of the order Diplostomida were recovered as sister groups of the order Plagiorchiida, indicating that the order Diplostomida is paraphyletic. This is the first mt genome of any member of the superfamily Brachylaimoidea and should represent a rich source of genetic markers for molecular epidemiological, population genetic and phylogenetic studies of parasitic flukes of socio-economic importance in poultry.

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.

Genetic differences among Haplorchis taichui populations in Indochina revealed by mitochondrial COX1 sequences

Haplorchis taichui is an intestinal heterophyid fluke that is pathogenic to humans. It is widely distributed in Asia, with a particularly high prevalence in Indochina. Previous work revealed that the lack of gene flow between three distinct populations of Vietnamese H. taichui can be attributed to their geographic isolation with no interconnected river basins. To test the hypothesis that interconnected river basins allow gene flow between otherwise isolated populations of H. taichui, as previously demonstrated for another trematode, Opisthorchis viverrini, we compared the genetic structures of seven populations of H. taichui from various localities in the lower Mekong Basin, in Thailand and Laos, with those in Vietnam, using the mitochondrial cytochrome c oxidase subunit 1 (COX1) gene. To determine the gene flow between these H. taichui populations, we calculated their phylogenetic relationships, genetic distances and haplotype diversity. Each population showed very low nucleotide diversity at this locus. However, high levels of genetic differentiation between the populations indicated very little gene flow. A phylogenetic analysis divided the populations into four clusters that correlated with the country of origin. The negligible gene flow between the Thai and Laos populations, despite sharing the Mekong Basin, caused us to reject our hypothesis. Our data suggest that the distribution of H. taichui populations was incidentally associated with national borders.

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.

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.

Phylogeography of Tetrancistrum nebulosi (Monogenea, Dactylogyridae) on the host of mottled spinefoot (Siganus fuscescens) in the South China Sea, inferred from mitochondrial COI and ND2 genes

To examine the phylogeographical pattern of Tetrancistrum nebulosi (Monogenea, Dactylogyridae) in the South China Sea, fragments of mitochondrial cytochrome c oxidase subunit I and NADH dehydrogenase subunit 2 genes were obtained for 220 individuals collected from 8 localities along the southeast coast of China and 1 locality in Terengganu, Malaysia. Based on these two genes, two and three distinct clades with geographic signals were revealed on the phylogenetic trees respectively. The divergence between these clades was estimated to occur in the late Pleistocene. Analysis of molecular variance and pairwise F_ST suggested a high rate of gene flow among individuals sampled from the Chinese coast, but with obvious genetic differentiation from the Malaysian population. Mismatch distribution and neutrality tests indicated that the T. nebulosi population experienced expansion in Pleistocene low sea level periods. Vicariance was considered to account for the genetic divergence between Chinese and Malaysian populations, while sea level fluctuations and mainland-island connections during glacial cycles were associated with the slight genetic divergence between the populations along the mainland coast of China and those off Sanya. On the contrary, oceanographic circulations and host migration could lead to genetic homogeneity of populations distributed along the mainland coast of China.

Sequence variation in two mitochondrial DNA regions and internal transcribed spacer among isolates of the nematode Oesophagostomum asperum originating from goats in Hunan Province, China

The present study examined sequence variability in two mitochondrial DNA (mtDNA) regions, namely cytochrome c oxidase subunit 1 (cox1) and NADH dehydrogenase subunit 1 (nad1), and internal transcribed spacer (ITS) of nuclear ribosomal DNA (rDNA) among Oesophagostomum asperum isolates from goats in Hunan Province, China. A portion of the cox1 (pcox1), nad1 (pnad1) genes and the ITS (ITS1+5.8S rDNA+ITS2) rDNA were amplified by polymerase chain reaction (PCR) separately from adult O. asperum individuals and the representative amplicons were subjected to sequencing from both directions. The lengths of pcox1, pnad1 and ITS rDNA were 366 bp, 681 bp and 785 bp, respectively. The A+T contents of gene sequences were 71.5–72% for pcox1, 73.7–74.2% for pnad1 and 58–58.8% for ITS rDNA. Intra-specific sequence variations within O. asperum were 0–1.6% for pcox1, 0–1.9% for pnad1 and 0–1.7% for ITS rDNA, while inter-specific sequence differences among members of the genus Oesophagostomum were significantly higher, being 11.1–12.5%, 13.3–17.7% and 8.5–18.6% for pcox1, pnad1 and ITS rDNA, respectively. Phylogenetic analyses using combined sequences of pcox1 and pnad1, with three different computational algorithms (Bayesian inference, maximum likelihood and maximum parsimony), revealed distinct groups with high statistical support. These findings demonstrated the existence of intra-specific variation in mtDNA and rDNA sequences among O. asperum isolates from goats in Hunan Province, China, and have implications for studying molecular epidemiology and population genetics of O. asperum.

Phylogenetic analysis of Pomacea canaliculata isolates collected from rice fields in different origins of China by combined mitochondrial 12S and 16S genes

To study the genetic relationships of Pomacea canaliculata collected from rice fields in China, the mitochondrial (mt) 12S and 16S of 9 P. canaliculata isolates from 5 southern provinces in China were sequenced and analyzed. The intra-specific sequence variations of P. canaliculata were 0-1.1% for 12S and 0--0.6% for 16S, while the inter-specific variations among common Pomacea species in mt 12S and 16S were 3.0-11.7% and 2.3-10.1%, respectively. Phylogenetic analysis based on combined sequences of mt 12S and 16S revealed complex genetic structure of P. canaliculata in China. Two phylogenetic groups of P. canaliculata were indicated in China with one group sistered to P. canaliculata isolates from USA, and two groups were even found in the same province. The phylogenetic relationships of Pomacea spp. also could be effectively inferred by combined sequences of mt 12S and 16S. These findings provided basic information for further study of population genetics and diffusion pattern of P. canaliculata in China as well as in the world.

Characterization of Haemaphysalis flava (Acari: Ixodidae) from Qingling subspecies of giant panda (Ailuropoda melanoleuca qinlingensis) in Qinling Mountains (Central China) by morphology and molecular markers

Tick is one of important ectoparasites capable of causing direct damage to their hosts and also acts as vectors of relevant infectious agents. In the present study, the taxa of 10 ticks, collected from Qinling giant pandas (Ailuropoda melanoleuca qinlingensis) in Qinling Mountains of China in April 2010, were determined using morphology and molecular markers (nucleotide ITS2 rDNA and mitochondrial 16S). Microscopic observation demonstrated that the morphological features of these ticks were similar to Haemaphysalis flava. Compared with other Haemaphysalis species, genetic variations between Haemaphysalis collected from A. m. qinlingensis and H. flava were the lowest in ITS2 rDNA and mitochondrial 16S, with sequence differences of 2.06%–2.40% and 1.30%–4.70%, respectively. Phylogenetic relationships showed that all the Haemaphysalis collected from A. m. qinlingensis were grouped with H. flava, further confirmed that the Haemaphysalis sp. is H. flava. This is the first report of ticks in giant panda by combining with morphology and molecular markers. This study also provided evidence that combining morphology and molecular tools provide a valuable and efficient tool for tick identification.

Genetic variability in the three mitochondrial genes among Oesophagostomum asperum isolates from different regions in Shaanxi and Hunan Provinces, China

The present study examined sequence variations in three mitochondrial DNA (mtDNA) regions, namely, NADH dehydrogenase subunit 5 (nad5), adenosine triphosphate subunit 6 (atp6) and cytochrome c oxidase subunit 3 (cox3), among Oesophagostomum asperum isolates from different regions in Shaanxi and Hunan provinces, China. The lengths for partial sequences of nad5 (pnad5), atp6 (patp6) and cox3 (pcox3) were 427 bp, 381 bp and 337 bp, respectively. The intra-specific sequence variations among all O. asperum samples were 0-2.11%, 0-1.84% and 0-1.48% for pnad5, patp6 and pcox3, respectively, while the inter-specific sequence differences among Oesophagostomum species in pig and small ruminants were 18-21.3% for pnad5, 18.3-24.5% for patp6 and 10.6-13.7% for pcox3. A phylogenetic analysis based on combined sequences of three mtDNA fragments indicated that all O. asperum isolates were grouped in one solid clade, and the Oesophagostomum spp. from pig were located in another clade. However, these mtDNA fragments could not reveal genetic relationships of geographical isolates of O. asperum in China. These results provided valuable information for studying population genetics of Oesophagostomum spp., and for controlling Oesophagostomum infection in animals as well as humans.

Genetic variability of Baylisascaris schroederi from the Qinling subspecies of the giant panda in China revealed by sequences of three mitochondrial genes

The present study examined the variations in three mitochondrial (mt) DNA sequences, namely cytochrome b (cytb), cytochrome c oxidase subunit 3 (cox3) and NADH dehydrogenase subunit 5 (nad5), among Baylisascaris schroederi isolates from the Qinling subspecies of the giant panda in Shaanxi province, northwestern China. No differences in length were detected in the three mt fragments from different isolates. The intra-specific sequence variations within all B. schroederi samples were 0-2.6% for pcytb, 0-1.8% for pcox3 and 0-2.1% for pnad5, while the inter-specific sequence differences among members of the genus Baylisascaris were 8.2-15.2%, 6.2-15.9% and 8.4-16.0% for pcytb, pcox3, pnad5, respectively. A phylogenetic analysis of the combined sequences of pcytb, pcox3 and pnad 5 showed that all B. schroederi samples in the present study were located in two large clusters, with one cluster containing samples from giant pandas in Sichuan province. These findings provide basic information for further study of molecular epidemiology and control of B. schroederi infection in the Qinling subspecies of the giant panda and throughout China.

Genetic variability among Dicrocoelium dendriticum isolates from different regions in Shaanxi Province, China revealed by sequences of three mitochondrial genes

The genetic variations in three mitochondrial DNA (mtDNA) regions, namely portion of cytochrome c oxidase subunit 1 (pcox1), NADH dehydrogenase subunit 1 (pnad1) and cytochrome b (pcytb), were examined for Dicrocoelium dendriticum samples isolated from different origins in Shaanxi Province, northwestern China. The intra-specific sequence differences within D. dendriticum samples were 0-0.52% for pcox1, 0-0.73% for pnad1 and 0-0.58% for pcytb. Phylogenetic analyses based on combined sequences of three mtDNA showed that all D. dendriticum samples were clustered together in same clade of Paragonimus westermani. But the phylogenetic trees could not reveal geographically genetic relationships of D. dendriticum isolates in this province. These findings will provide basic information for further study of molecular epidemiology and control of D. dendriticum infection in this province as well as in China.

Molecular characterization of Cyclospora-like organisms from golden snub-nosed monkeys in Qinling Mountain in Shaanxi province, northwestern China

Cyclospora spp. have been identified as one of the most important intestinal pathogens causing protracted diarrhea in animals and human beings. To determine the Cyclospora species in the non-human primate Rhinopithecus roxellanae, a total of 71 fecal samples from 19 endangered snub-nosed monkeys in Shaanxi province were collected and examined using Sheater's sugar flotation technique and by sequencing the fragments of 18S rDNA. Only two Cyclospora isolates from 2 golden snub-nosed monkeys (R. roxellanae) were obtained and identified between July 2011 and August of 2012. The sequences of the 18S rDNA for the two Cyclospora isolates were 477 bp, with no nucleotide variation between them. Phylogenetic analysis based on the 18S rDNA sequences revealed that the two Cyclospora isolates were posited into the clade Cyclospora spp. and sistered to C. colobi. These results first showed that Cyclospora infection occurred in R. roxellanae in hot and rainy weather, which would provide useful information for further understanding the molecular epidemiology of Cyclospora spp. and the control of Cyclospora infection in non-human primates as well as in human beings.

A foreign invader or a reclusive native? DNA bar coding reveals a distinct European lineage of the zoonotic parasite Schistosoma turkestanicum (syn. Orientobilharzia turkestanicum ())

Natural foci of Schistosoma turkestanicum (syn. Orientobilharzia turkestanicum) has been identified in the Gemenc Forest regions of Hungary utilising red deer as the definitive host. In order to identify the origins of this parasite in Europe standard DNA bar coding techniques were employed to sequence fragments of the cytochrome oxidase 1 (cox1) and the nuclear ribosomal internal transcribed region (ITS) from 10 individual adult male worms. Phylogenetic reconstruction using maximum likelihood phylogenetic reconstruction and haplotype networks of the cox1 showed all the worms to be of a distinct unique Hungarian lineage although some ITS haplotypes were shared with worms from populations in China and Iran. Molecular clock analysis suggests an early divergence event around 270,000years before present (YBP) between all S. turkestanicum populations giving rise to the Chinese, Iranian and Hungarian lineages. However, divergence of the sequences within the Hungarian population appears to have occurred approximately 63,000 YBP suggesting a long established population of S. turkestanicum in Europe. This suggests that the Hungarian population of S. turkestanicum has been native since the Ice Age and probably established itself during the last interglacial period as red deer moved into Europe from North Africa and the Middle East. This may also indicate that the parasite may have unknown populations established in several other countries in Eastern, Central and Southern Europe.

Sequence variability in four mitochondrial genes among Bunostomum trigonocephalum isolates from four provinces in China

The present study examined sequence variability in four mitochondrial genes, namely cytochrome c oxidase subunit (cox1), cytochrome b (cytb) and NADH dehydrogenase subunits 1 and 5 (nad1 and nad5), among Bunostomum trigonocephalum isolates from four different geographic regions in China. Ten B. trigonocephalum samples were collected from each of the four provinces (Heilongjiang, Jilin, Shaanxi and Yunnan), China. A part of the cox1 (pcox1), cytb (pcytb), nad1 and nad5 genes (pnad1 and pnad5) were amplified separately from individual hookworms by polymerase chain reaction (PCR) and were subjected to direct sequencing in order to define sequence variations and their phylogenetic relationships. The intra-specific sequence variations within B. trigonocephalum were 0–1.9% for pcox1, 0–2.0% for pcytb, 0–1.6% for pnad1 and 0–1.7% for pnad5. The A+T contents of the sequences were 69.6–70.4% (pcox1), 71.9–72.7 (pcytb), 70.4–71.1% (pnad1) and 72.0–72.6% (pnad5). However, the inter-specific sequence differences among members of the family Ancylostomatidae were significantly higher, being 12.1–14.2% for pcox1, 13.7–16.0 for cytb, 17.6–19.4 for nad1 and 16.0–21.6 for nad5. Phylogenetic analyses based on the combined partial sequences of cox1, cytb, nad1 and nad5 using three inference methods, namely Bayesian inference (Bayes), maximum likelihood (ML) and maximum parsimony (MP), revealed that all the B. trigonocephalum samples form monophyletic groups, but samples from the same geographical origin did not always cluster together, suggesting that there was no obvious geographical distinction within B. trigonocephalum based on sequences of the four mtDNA genes. These results demonstrated the existence of low-level intra-specific variation in mitochondrial DNA (mtDNA) sequences among B. trigonocephalum isolates from different geographic regions.

Sequence variation in three mitochondrial DNA genes among isolates of Ascaridia galli originating from Guangdong, Hunan and Yunnan provinces, China

The present study examined sequence variation in three mitochondrial DNA (mtDNA) genes, namely cytochrome c oxidase subunit 3 (cox3) and NADH dehydrogenase subunits 1 and 4 (nad1 and nad4), among Ascaridia galli isolates from different geographical localities in China. A portion of cox3 (pcox3), nad1 (pnad1) and nad4 (pnad4) genes were amplified by polymerase chain reaction (PCR) separately from adult A. galli individuals and the amplicons were subjected to sequencing from both directions. The length of the sequences of pcox3, pnad1 and pnad4 were 408 bp, 471 bp and 333 bp, respectively. The intraspecific sequence variations within A. galli were 0-1.7% for pcox3, 0-2.8% for pnad1 and 0-3.4% for pnad4. The A+T contents of the sequences were 67.16-67.65% (pcox3), 67.09-67.94% (pnad1) and 69.91-71.77% (pnad4). The interspecific sequence differences among members of the Ascaridida were significantly higher, being 13.2-30.9%, 12.8-29.0% and 15.1-34.1% for pcox3, pnad1 and pnad4, respectively. Phylogenetic analyses using combined sequences of pcox3, pnad1 and pnad4, with three different computational algorithms (Bayesian analysis, maximum likelihood and maximum parsimony), all revealed distinct groups with high statistical support. These findings demonstrated the existence of intraspecific variation in mitochondrial DNA (mtDNA) sequences among A. galli isolates from different geographical regions in China, and have implications for studying molecular epidemiology and population genetics of A. galli.

Sequence variability in three mitochondrial genes between the two pig nodule worms Oesophagostomum dentatum and O. quadrispinulatum

In this study, sequence variation in three mitochondrial DNA regions, namely cytochrome c oxidase subunit (cox1) and NADH dehydrogenase subunits 1 and 4 (nad1 and nad4), between Oesophagostomum dentatum and O. quadrispinulatum isolated from pigs in different geographical origins in Mainland China was examined, and their phylogenetic relationships were reconstructed. A partial of the cox1 (pcox1), nad1, and nad4 genes (pnad1 and pnad4) were amplified separately from individual nodule worms by PCR and were subjected to direct sequencing in order to define sequence variations. While the intraspecific sequence variations within each of the two species were 0.3-5.2% for pcox1, 0-4.9% for pnad1, and 0-7.1% for pnad4, the interspecific sequence differences were significantly higher, being 10.7-13.4% for pcox1, 11-14.6% for pnad1, and 14.9-18% for pnad4, respectively. There were a number of nucleotide positions in the pcox1, pnad1, and pnad4 sequences with no apparent intraspecific variation but distinct interspecific differences among those samples of Oesophagostomum spp. examined, which may be used as genetic makers for the identification and differentiation of the Oesophagostomum spp. Phylogenetic analyses using three inference methods, namely Bayesian inference, maximum likelihood, and maximum parsimony based on the combined sequences of pcox1, pnad1, and pnad4 revealed that the O. dentatum and O. quadrispinulatum form monophyletic groups, respectively. These findings demonstrated clearly the usefulness of the three mitochondrial sequences for studying systematics, population genetic structures, and the molecular ecology of Oesophagostomum spp.

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.

Genetic variability among Trichuris ovis isolates from different hosts in Guangdong Province, China revealed by sequences of three mitochondrial genes

This study examined sequence variation in three mitochondrial DNA (mtDNA) regions, namely cytochrome c oxidase subunit 1 (cox1), NADH dehydrogenase subunit 5 (nad5) and cytochrome b (cytb), among Trichuris ovis isolates from different hosts in Guangdong Province, China. A portion of the cox1 (pcox1), nad5 (pnad5) and cytb (pcytb) genes was amplified separately from individual whipworms by PCR, and was subjected to sequencing from both directions. The size of the sequences of pcox1, pnad5 and pcytb was 618, 240 and 464 bp, respectively. Although the intra-specific sequence variations within T. ovis were 0-0.8% for pcox1, 0-0.8% for pnad5 and 0-1.9% for pcytb, the inter-specific sequence differences among members of the genus Trichuris were significantly higher, being 24.3-26.5% for pcox1, 33.7-56.4% for pnad5 and 24.8-26.1% for pcytb, respectively. Phylogenetic analyses using combined sequences of pcox1, pnad5 and pcytb, with three different computational algorithms (maximum likelihood, maximum parsimony and Bayesian inference), indicated that all of the T. ovis isolates grouped together with high statistical support. These findings demonstrated the existence of intra-specific variation in mtDNA sequences among T. ovis isolates from different hosts, and have implications for studying molecular epidemiology and population genetics of T. ovis.

The second transcribed spacer rDNA sequence: an effective genetic marker for inter-species phylogenetic analysis of trematodes in the order Strigeata

In the present study, the second nuclear internal transcribed spacer (ITS-2) rDNA of Schistosoma japonicum isolates in mainland China was amplified, sequenced, and assessed for inferring the intra- and inter-species phylogenetic relationships of trematodes in the order Strigeata. The fragment containing ITS-2 rDNA was obtained from 24 S. japonicum isolates from eight epidemic provinces in mainland China. The length polymorphisms were observed among these ITS-2 rDNA sequences, ranging from 343 to 346 bp, and the intra- and inter-population variations in ITS-2 sequence were 0.0-2.1% among S. japonicum isolates in China. Phylogenetic analyses using the maximum parsimony and maximum likelihood methods revealed that the ITS-2 rDNA sequence is not a suitable marker for studying inter- and intra-population variation in S. japonicum. However, phylogenetic analysis of trematodes in the order Strigeata indicated that the ITS-2 rDNA sequence provides an effective molecular marker for studying inter-species phylogenetic relationships among trematodes in this order.

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.

Genetic variation among Clonorchis sinensis isolates from different geographic regions in China revealed by sequence analyses of four mitochondrial genes

The present study examined sequence variation in four mitochondrial (mt) genes, namely cytochrome c oxidase subunits 1 (cox1) and 2 (cox2), and NADH dehydrogenase subunits 1 and 2 (nad1 and nad2) among Clonorchis sinensis isolates from different endemic regions in China, and their phylogenetic relationships with other zoonotic trematodes were reconstructed. A portion of the cox1 and cox2 genes (pcox1 and pcox2), and nad1 and nad2 genes (pnad1 and pnad2) were amplified separately from individual liver flukes by polymerase chain reaction (PCR) and the amplicons were subjected to sequencing from both directions. The intra-specific sequence variations within C. sinensis were 0–1.6% for pcox1, 0–1.4% for pcox2, 0–0.9% for pnad1 and 0–1.0% for pnad2. Phylogenetic analyses based on the combined sequences of pcox1, pcox2, pnad1 and pnad2 revealed that all the C. sinensis isolates grouped together and were closely related to Opisthorchis felineus. These findings revealed the existence of intra-specific variation in mitochondrial DNA (mtDNA) sequences among C. sinensis isolates from different geographic regions, and demonstrated that mtDNA sequences provide reliable genetic markers for phylogenetic studies of zoonotic trematodes.

Characterization of Baylisascaris schroederi from Qinling subspecies of giant panda in China by the first internal transcribed spacer (ITS-1) of nuclear ribosomal DNA

In the present study, a total of 20 nematode isolates, (including 10 male and 10 female worms) representing Baylisascaris schroederi from 5 Qinling subspecies of giant pandas (Ailuropoda melanoleuca) in Shaanxi Province of China, were characterized and grouped genetically by the first internal transcribed spacer (ITS-1) of nuclear ribosomal DNA (rDNA). The rDNA fragment spanning 3' end of 18S rDNA, complete ITS-1 rDNA, and 5' end of 5.8S rDNA were amplified and sequenced. The sequence variability in ITS-1 rDNA was examined within B. schroederi and among parasites in order Ascaridata available in GenBank™, and their phylogenetic relationships were also reconstructed. The sequences of ITS-1 rDNA for all the B. schroederi isolates were 427 bp in length, with no genetic variation detected among these isolates. Phylogenetic analyses based on the ITS-1 rDNA sequences revealed that all the male and female B. schroederi isolates sequenced in the present study were posited into the clade of genus Baylisascaris, sistered to zoonotic nematodes in genus Ascaris, and the ITS-1 rDNA sequence could distinguish different species in order Ascaridata. These results showed that the ITS-1 rDNA provides a suitable molecular marker for the inter-species phylogenetic analysis and differential identification of nematodes in order Ascaridata.

Sequence variability in three mitochondrial DNA regions of Spirometra erinaceieuropaei spargana of human and animal health significance

Sequence variability in three mitochondrial DNA (mtDNA) regions, namely cytochrome c oxidase subunit 3 (cox3), NADH dehydrogenase subunits 1 and 4 (nad1 and nad4) in Spirometra erinaceieuropaei spargana from different geographical regions in China was examined. A portion of each of the cox3 (pcox3), nad1 (pnad1) and nad4 genes (pnad4) were amplified separately from individual S. erinaceieuropaei spargana by polymerase chain reaction (PCR). Representative amplicons were subjected to sequencing in order to estimate sequence variability. The sequences of pcox3, pnad1 and pnad4 were 541, 607 and 847 bp in length, respectively. The A+T contents of the sequences were 68.39–68.76% (pcox3), 63.76–64.91% (pnad1) and 67.18–67.77% (pnad4), respectively, while the intra-specific sequence variations within each of the S. erinaceieuropaei spargana were 0–1.5% for pcox3, 0–2.8% for pnad1 and 0–2.7% for pnad4. Phylogenetic analysis using neighbour joining (NJ), maximum likelihood (ML) and maximum parsimony (MP) methods, indicated that all the spargana isolates in Hunan Province represented S. erinaceieuropaei. These findings demonstrated clearly the usefulness of the three mtDNA sequences for population genetics studies of S. erinaceieuropaei spargana of human and animal health significance.

Sequence variability in two mitochondrial DNA regions and internal transcribed spacer among three cestodes infecting animals and humans from China

Sequence variability in two mitochondrial DNA (mtDNA) regions, namely cytochromecoxidase subunit 1 (cox1) and NADH dehydrogenase subunit 4 (nad4), and internal transcribed spacer (ITS) of rDNA among and within three cestodes,Spirometra erinaceieuropaei,Taenia multicepsandTaenia hydatigena, from different geographical origins in China was examined. A portion of thecox1 (pcox1),nad4 genes (pnad4) and the ITS (ITS1+5.8S rDNA+ITS2) were amplified separately from individual cestodes by polymerase chain reaction (PCR). Representative amplicons were subjected to sequencing in order to estimate sequence variability. While the intra-specific sequence variations within each of the tapeworm species were 0–0.7% for pcox1, 0–1.7% for pnad4 and 0.1–3.6% for ITS, the inter-specific sequence differences were significantly higher, being 12.1–17.6%, 18.7–26.2% and 31–75.5% for pcox1, pnad4 and ITS, respectively. Phylogenetic analyses based on the pcox1 sequence data revealed thatT. multicepsandT. hydatigenawere more closely related to the other members of theTaeniagenus, andS. erinaceieuropaeiwas more closely related to the other members of theSpirometragenus. These findings demonstrated clearly the usefulness of mtDNA and rDNA sequences for population genetic studies of these cestodes of socio-economic importance.

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.

Specific detection of Angiostrongylus cantonensis in the snail Achatina fulica using a loop-mediated isothermal amplification (LAMP) assay

Angiostrongylus cantonensis, a rat lungworm, can cause eosinophilic meningitis and angiostrongyliasis in humans following ingestion of contaminated foods or intermediate/paratenic hosts with infective larvae. The snail Achatina fulica is one of the important intermediate hosts of A. cantonensis and is commonly eaten by humans in some countries. In the present study, we developed a loop-mediated isothermal amplification (LAMP) method for the specific detection of A. cantonensis in Ac. fulica. Primers for LAMP were designed based on the first internal transcribed spacer (ITS-1) of nuclear ribosomal DNA (rDNA) of A. cantonensis. Specificity tests showed that only the products of A. cantonensis were detected when DNA samples of A. cantonensis and the heterologous control samples Anisakis simplex s.s, Trichuris trichiura, Toxocara canis, Trichinella spiralis and Ascaris lumbricoides were amplified by LAMP. Sensitivity evaluation indicated that the LAMP assay is 10 times more sensitive than the conventional polymerase chain reaction (PCR) assay. The established LAMP assay is rapid, inexpensive and easy to be performed. It can be used in clinical applications for rapid and sensitive detection of A. cantonensis in snails, which has implications for the effective control of angiostrongyliasis.

Genetic diversity and relatedness of Fasciola spp. isolates from different hosts and geographic regions revealed by analysis of mitochondrial DNA sequences

The present study examined sequence variability in a portion of the mitochondrial cytochrome c oxidase subunit 1 (pcox1) and NADH dehydrogenase subunits 4 and 5 (pnad4 and pnad5) among 39 isolates of Fasciola spp., from different hosts from China, Niger, France, the United States of America, and Spain; and their phylogenetic relationships were re-constructed. Intra-species sequence variations were 0.0-1.1% for pcox1, 0.0-2.7% for pnad4, and 0.0-3.3% for pnad5 for Fasciola hepatica; 0.0-1.8% for pcox1, 0.0-2.5% for pnad4, and 0.0-4.2% for pnad5 for Fasciola gigantica, and 0.0-0.9% for pcox1, 0.0-0.2% for pnad4, and 0.0-1.1% for pnad5 for the intermediate Fasciola form. Whereas, nucleotide differences were 2.1-2.7% for pcox1, 3.1-3.3% for pnad4, and 4.2-4.8% for pnad5 between F. hepatica and F. gigantica; were 1.3-1.5% for pcox1, 2.1-2.9% for pnad4, 3.1-3.4% for pnad5 between F. hepatica and the intermediate form; and were 0.9-1.1% for pcox1, 1.4-1.8% for pnad4, 2.2-2.4% for pnad5 between F. gigantica and the intermediate form. Phylogenetic analysis based on the combined sequences of pcox1, pnad4 and pnad5 revealed distinct groupings of isolates of F. hepatica, F. gigantica, or the intermediate Fasciola form irrespective of their origin, demonstrating the usefulness of the mtDNA sequences for the delineation of Fasciola species, and reinforcing the genetic evidence for the existence of the intermediate Fasciola form.

The ribosomal intergenic spacer (IGS) region in Schistosoma japonicum: structure and comparisons with related species

The intergenic spacer (IGS) between the 28S and 18S ribosomal RNA genes was PCR-amplified, sequenced and characterized for Schistosoma japonicum from mainland China, and compared with those of other Schistosoma species. Excluding flanking portions of the 28S and 18S genes, the IGS in the longest sequenced amplicon from S. japonicum IGS was 1457bp in length. However, intra-specific and intra-individual variation was noted. The IGS region of S. japonicum is strikingly different in structure from those of African Schistosoma species for which data are available. S. japonicum has a shorter IGS and largely lacks a long region of complex repeats seen in the African species. However, careful comparisons with African species highlighted the presence of a few shared repeat motifs that were not apparent from study of African species only. Such motifs presumably have functional significance. Discovery of such motifs may in general be aided by comparisons of relatively distant taxa rather than of sibling taxa.