Molecular analyses confirm the coexistence of Fasciola gigantica and parthenogenetic Fasciola in the Philippines

Assessment of genetic markers for multilocus sequence typing (MLST) of Fasciola isolates from Iran

BACKGROUND

Several markers have been described to characterise the population structure and genetic diversity of Fasciola species (Fasciola hepatica (F. hepatica) and Fasciola gigantica (F. gigantica). However, sequence analysis of a single genomic locus cannot provide sufficient resolution for the genetic diversity of the Fasciola parasite whose genomes are ∼1.3 GB in size.

OBJECTIVES

To gain a better understanding of the gene diversity of Fasciola isolates from western Iran and to identify the most informative markers as candidates for epidemiological studies, five housekeeping genes were evaluated using a multilocus sequence typing (MLST) approach.

METHODS

MLST analysis was developed based on five genes (ND1, Pepck, Pold, Cyt b and HSP70) after genomic DNA extraction, amplification and sequencing. Nucleotide diversity and phylogeny analysis were conducted on both concatenated MLST loci and each individual locus. A median joining haplotype network was created to examine the haplotypes relationship among Fasciola isolates.

RESULTS

Thirty-three Fasciola isolates (19 F. hepatica and 14 F. gigantica) were included in the study. A total of 2971 bp was analysed for each isolate and 31 sequence types (STs) were identified among the 33 isolates (19 for F. hepatica and 14 for F. gigantica isolates). The STs produced 44 and 42 polymorphic sites and 17 and 14 haplotypes for F. hepatica and F. gigantica, respectively. Haplotype diversity was 0.982 ± 0.026 and 1.000 ± 0.027 and nucleotide diversity was 0.00200 and 0.00353 ± 0.00088 for F. hepatica and F. gigantica, respectively. There was a high degree of genetic diversity with a Simpson's index of diversity of 0.98 and 1 for F. hepatica and F. gigantica, respectively. While HSP70 and Pold haplotypes from Fasciola species were separated by one to three mutational steps, the haplotype networks of ND1 and Cyt b were more complex and numerous mutational steps were found, likely due to recombination.

CONCLUSIONS

Although HSP70 and Pold genes from F. gigantica were invariant over the entire region of sequence coverage, MLST was useful for investigating the phylogenetic relationship of Fasciola species. The present study also provided insight into markers more suitable for phylogenetic studies and the genetic structure of Fasciola parasites.

Dispersal direction of Malaysian Fasciola gigantica from neighboring southeast Asian countries inferred using mitochondrial DNA analysis

Fasciola gigantica and hybrid Fasciola flukes, responsible for the disease fasciolosis, are found in Southeast Asian countries. In the present study, we performed molecular species identification of Fasciola flukes distributed in Terengganu, Malaysia using multiplex PCR for phosphoenolpyruvate carboxykinase (pepck) and PCR-restriction fragment length polymorphism (RFLP) for DNA polymerase delta (pold). Simultaneously, phylogenetic analysis based on mitochondrial NADH dehydrogenase subunit 1 (nad1) was performed for the first time on Malaysian Fasciola flukes to infer the dispersal direction among neighboring countries. A total of 40 flukes used in this study were identified as F. gigantica. Eight nad1 haplotypes were identified in the F. gigantica population of Terengganu. Median-joining network analysis revealed that the Malaysian population was related to those obtained from bordering countries such as Thailand and Indonesia. However, genetic differentiation was detected using population genetics analyses. Nevertheless, the nucleotide diversity (π) value suggested that F. gigantica with the predominant haplotypes was introduced into Malaysia from Thailand and Indonesia. The dispersal direction suggested by population genetics in the present study may not be fully reliable since Fasciola flukes were collected from a single location in one state of Malaysia. Further studies analyzing more samples from many locations are required to validate the dispersal direction proposed herein.

The causative agents of fascioliasis in animals and humans: Parthenogenetic Fasciola in Asia and other regions

Parthenogenetic Fasciola is the causative agent of fascioliasis in animals and humans and is widely distributed in Asian countries, such as Japan, South Korea, China, Vietnam, Thailand, the Philippines, Myanmar, Bangladesh, Nepal, and India. Parthenogenetic Fasciola geographically originated from central and eastern China, where it exists between the habitats of Fasciola hepatica and Fasciola gigantica; it likely appeared thousands of years ago following hybridization between F. hepatica and F. gigantica. Parthenogenetic Fasciola consists of diploids and triploids that possess nuclear genome of both F. hepatica and F. gigantica and mitochondrial genome of either F. hepatica or F. gigantica. Maternal parents of parthenogenetic Fasciola are either F. hepatica having Fh-C4 haplotype or F. gigantica having Fg-C2 haplotype in mitochondrial NADH dehydrogenase subunit 1 (ND1) nucleotide sequences. Parthenogenetic Fasciola flukes with the Fh-C4 haplotype have spread from China to South Korea and Japan, whereas the flukes with the Fg-C2 haplotype have not only spread to Korea and Japan but also southward to Vietnam, Thailand, the Philippines, Myanmar, Bangladesh, Nepal, and India. Parthenogenetic Fasciola can be distinguished from F. hepatica and F. gigantica using combinational DNA sequence analysis of nuclear phosphoenolpyruvate carboxykinase (pepck) and DNA polymerase delta (pold) along with mitochondrial ND1 markers. The establishment of parthenogenetic Fasciola is expected as follows: parthenogenetic diploids with the Fh-C4 and Fg-C2 haplotypes first appeared based on single or multiple interspecific hybridization events; subsequently, parthenogenetic triploids emerged via backcross events between the maternal parthenogenetic diploid and either paternal bisexual F. hepatica or F. gigantica. Parthenogenetic Fasciola diploids and triploids then survived for thousands of years by clonal parthenogenetic reproduction, and generated descendants with ND1 haplotypes, which were derived from the Fh-C4 and Fg-C2 due to nucleotide substitution. Thus, the emergence of parthenogenetic Fasciola may be due to extremely uncommon and accidental events. Parthenogenetic Fasciola should be treated as a new asexual hybrid species.