Fasciola hepatica demonstrates high levels of genetic diversity,a lack of population structure and high gene flow: possible implications for drug resistance

Fasciola hepatica, the liver fluke, is a trematode parasite of considerable economic importance to the livestock industry and is a re-emerging zoonosis that poses a risk to human health in F. hepatica-endemic areas worldwide. Drug resistance is a substantial threat to the current and future control of F. hepatica, yet little is known about how the biology of the parasite influences the development and spread of resistance. Given that F. hepatica can self-fertilise and therefore inbreed, there is the potential for greater population differentiation and an increased likelihood of recessive alleles, such as drug resistance genes, coming together. This could be compounded by clonal expansion within the snail intermediate host and aggregation of parasites of the same genotype on pasture. Alternatively, widespread movement of animals that typically occurs in the UK could promote high levels of gene flow and prevent population differentiation. We identified clonal parasites with identical multilocus genotypes in 61% of hosts. Despite this, 84% of 1579 adult parasites had unique multilocus genotypes, which supports high levels of genotypic diversity within F. hepatica populations. Our analyses indicate a selfing rate no greater than 2%, suggesting that this diversity is in part due to the propensity for F. hepatica to cross-fertilise. Finally, although we identified high genetic diversity within a given host, there was little evidence for differentiation between populations from different hosts, indicating a single panmictic population. This implies that, once those emerge, anthelmintic resistance genes have the potential to spread rapidly through liver fluke populations.     

Evidence of population structuring following population genetic analyses of Fasciola hepatica from Argentina

Fasciola hepatica, the liver fluke, is a trematode parasite that causes disease of economic importance in livestock. As a zoonosis this parasite also poses a risk to human health in areas where it is endemic. Population genetic studies can reveal the mechanisms responsible for genetic structuring (non-panmixia) within parasite populations and provide valuable insights into population dynamics, which in turn enables theoretical predictions of evolutionary dynamics such as the evolution of drug resistance. Here we genotyped 320 F. hepatica collected from 14 definitive hosts from four provinces in Argentina. STRUCTURE analysis indicated three population clusters, and principal coordinate analysis confirmed this, showing population clustering across provinces. Similarly, pairwise F_ST values amongst all four provinces were significant, with standardised pairwise F_ST (F′_ST) ranging from 0.0754 to 0.6327. Therefore, population genetic structure was evident across these four provinces in Argentina. However, there was no evidence of deviation from Hardy–Weinberg equilibrium, so it appears that within these sub-populations there is largely random mating. We identified 263 unique genotypes, which gave a clonal diversity of 82%. Parasites with identical genotypes, clones, accounted for 26.6% of the parasites studied and were found in 12 of the 14 hosts studied, suggesting some clonemate transmission.     

Identification and differentiation of Fasciola hepatica and Fasciola gigantica using a simple PCR-restriction enzyme method

Accurate morphological differentiation between the liver fluke species Fasciola hepatica and Fasciola gigantica is difficult. We evaluated PCR-restriction enzyme profiles of internal transcribed spacer 1 (ITS1) that could aid in their identification. Fifty F. hepatica and 30 F. gigantica specimens were collected from different hosts in three provinces of Iran. For DNA extraction, we crushed fragments of the worms between two glass slides as a new method to break down the cells. DNA from the crushed materials was then extracted with a conventional phenol-chloroform method and with the newly developed technique, commercial FTA cards. A primer pair was selected to amplify a 463-bp region of the ITS1 sequence. After sequencing 14 samples and in silico analysis, cutting sites of all known enzymes were predicted and TasI was selected as the enzyme that yielded the most informative profile. Crushing produced enough DNA for PCR amplification with both the phenol-chloroform and commercial FTA card method. The DNA extracted from all samples was successfully amplified and yielded a single sharp band of the expected size. Digestion of PCR products with TasI allowed us to distinguish the two species. In all samples, molecular identification was consistent with morphological identification. Our PCR-restriction enzyme profile is a simple, rapid and reliable method for differentiating F. hepatica and F. gigantica, and can be used for diagnostic and epidemiological purposes.     

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

Fasciola flukes collected from domestic buffalos and cattle in the Philippines were confirmed as Fasciola gigantica and parthenogenetic Fasciola based on DNA analyses of nuclear pepck and pold genes, and the mitochondrial ND1 gene. This study is the first to elucidate that F. gigantica and parthenogenetic Fasciola coexist in the Philippines with prevalences of 90.6% and 9.4%, respectively. The F. gigantica population showed a high genetic diversity with 25 ND1 haplotypes, suggesting that F. gigantica has existed in the Philippines for a long time. In contrast, parthenogenetic Fasciola flukes showed a single ND1 haplotype (Fsp-ND1-P1), which was identical to the founder haplotype, Fg-C2 of parthenogenetic Fasciola in China. These results indicate that parthenogenetic Fasciola in the Philippines is a recently introduced population from a neighboring continent.     

Molecular characterization of Fasciola flukes using mitochondrial 28S rRNA gene in Naimi Saudi sheep

Fasciolosis is a parasitic disease of medical and economic importance. This retrospective study was conducted on 110 Fasciola flukes collected from livers of 14 infected Naimi sheep slaughtered at Riyadh abattoir in Saudi Arabia during winter season of 2016. Collected specimens were analyzed for their species identification on the basis of partial sequences of mitochondrial 28S rRNA gene. Results have shown the presence of both Fasciola hepatica (F. hepatica) and Fasciola gigantica (F. gigantica) species. Where Fasciola hepatica was predominate (80%). Both intra-species and interspecies genetic distance was studied and results showed that the intraspecific variability among individuals of both species i.e., F. hepatica and F. gigantica, ranging between 0 and 1% while the interspecific diversity between F. hepatica and F. gigantica was only 1%. In conclusion, mitochondrial 28S rRNA gene is a proved as a good marker in identifying Fasciola of different species. Where, the F. hepatica and F. gigantica are present in sheep breed in Riyadh region, Saudi Arabia.     

Substantial genetic divergence between morphologically indistinguishable populations of Fasciola suggests the possibility of cryptic speciation

The liver flukes, Fasciola hepatica and Fasciola gigantica, are considered to be sister species and between them present a major threat worldwide to livestock production. In this study sequence data have been employed from informative regions of the nuclear and mitochondrial genomes of over 200 morphologically F. hepatica-like or F. gigantica-like flukes from Europe, sub-Saharan Africa and South Asia to assess genetic diversity. Evidence is presented for the existence of four well-separated clades: African gigantica-like flukes, Indian gigantica-like flukes, European hepatica-like flukes and African high-altitude hepatica-like flukes. Application of the Biological Species Concept to trematodes is problematic; however, the degree of separation between these groups was sufficient for them to be considered as distinct species using the four times rule for speciation.     

Molecular phylogenetic identification of Fasciola flukes in Nepal

Eighty-one Fasciola flukes collected from 8 districts in Nepal were analyzed for their species identification on the basis of their spermatogenic status and nuclear ribosomal internal transcribed spacer 1 (ITS1) and for their phylogenetic relation with Fasciola flukes from other Asian countries on the basis of the mitochondrial NADH dehydrogenase subunit 1 (nad1) gene. Sixty-one flukes (75.3%) were aspermic Fasciola sp., and 20 flukes (24.7%) were identified as Fasciola gigantica. All of the aspermic flukes displayed the Fh/Fg type in ITS1, which was predominant in aspermic Fasciola sp. from China, and most (60 flukes) displayed the Fsp-ND1-N1 haplotype in the nad1, which had an identical nucleotide sequence to the major haplotype (Fg-C2) of the aspermic flukes from China. These results suggest that aspermic Fasciola sp. was introduced into Nepal from China. Furthermore, the results of the diversity indices, neutrality indices, and median-joining network analysis with reference haplotypes from Asian countries suggest that aspermic Fasciola sp. rapidly expanded its distribution. In contrasts, F. gigantica displayed 10 nad1 haplotypes, which showed higher population diversity indices than the haplotypes of aspermic flukes, indicating that the F. gigantica population was clearly distributed in Nepal earlier than the aspermic Fasciola population. Although the F. gigantica haplotypes from Nepal formed a star-like phylogeny consisting of a main founder haplotype (Fg-ND1-N1), together with some F. gigantica haplotypes from Myanmar and Thailand, the Nepal population differed genetically from F. gigantica populations of neighboring countries as each country had distinct founder haplotype(s).     

Molecular characterization revealed Fasciola specimens in Ecuador are all Fasciola hepatica,none at all of Fasciola gigantica or parthenogenic Fasciola species

All 225 Fasciola flukes obtained from domestic animals (73 cattle, 7 sheep and 1 pig) of 18 distinct geographic areas in Ecuador-South America, were identified as Fasciola hepatica, based on molecular analyses of nuclear pepck and pold genes, and mitochondrial nad1gene as well as the morphological observation of sperm within the seminal vesicles. Fasciola gigantica and parthenogenic Fasciola forms endemic to Asian countries were not found in this study, although zebu cattle and water buffalos have introduced into South America from Asia; this could be due to the absence of suitable intermediate host snails. The results of pepck analysis using multiplex PCR developed previously showed that 32 of the flukes could not be confirmed as F. hepatica, suggesting that the method is unreliable for the accurate discrimination of F. hepatica, and that pepck gene of the species consists of multiple loci, not a single locus. The results of genetic diversity, phylogenetic, and network analyses based on mitochondrial nad1 sequences suggest that F. hepatica populations in South America, including Ecuador, formed from the ancestral F. hepatica individuals introduced into the continent along with anthropogenic movement of livestock infected with the species.     

A Case of Fasciola hepatica Infection Mimicking Cholangiocarcinoma and ITS-1 Sequencing of the Worm

Fascioliasis is a zoonotic infection caused by Fasciola hepatica or Fasciola gigantica. We report an 87-year-old Korean male patient with postprandial abdominal pain and discomfort due to F. hepatica infection who was diagnosed and managed by endoscopic retrograde cholangiopancreatography (ERCP) with extraction of 2 worms. At his first visit to the hospital, a gallbladder stone was suspected. CT and magnetic retrograde cholangiopancreatography (MRCP) showed an intraductal mass in the common bile duct (CBD) without proximal duct dilatation. Based on radiological findings, the presumed diagnosis was intraductal cholangiocarcinoma. However, in ERCP which was performed for biliary decompression and tissue diagnosis, movable materials were detected in the CBD. Using a basket, 2 living leaf-like parasites were removed. The worms were morphologically compatible with F. hepatica. To rule out the possibility of the worms to be another morphologically close species, in particular F. gigantica, 1 specimen was processed for genetic analysis of its ITS-1 region. The results showed that the present worms were genetically identical (100%) with F. hepatica but different from F. gigantica.     

Identification of new polymorphic positions in rDNA sequences of the “intermediate” Fasciola forms

Fascioliasis is a foodborne zoonotic disease generally caused by the parasitic flukes Fasciola gigantica and Fasciola hepatica in class Trematoda. An “intermediate” Fasciola forms between F. gigantica and F. hepatica has been shown to exist. However, the relationships among F. gigantica, F. hepatica, and “intermediate” Fasciola forms remain unclear. In this study, we found five new polymorphic positions in 18S and 28S rDNAs sequences of “intermediate” Fasciola forms. According to the high-throughput sequencing results, all known 16 polymorphic positions of “intermediate” Fasciola forms show a clear and consistent tendency for F. gigantica or F. hepatica, and the percentages of the most frequently occurring bases were different in specimens. In the three ITS sequence fragments, hybrid-type base combinations of the polymorphic positions were detected, and the percentages of the most frequent base combinations were different in specimens too. In addition, interestingly, the newly detected ITS-802 position was not a traditional polymorphic position in “intermediate” Fasciola forms, and the bases in ITS-802 position are not same as the allele bases of F. gigantica or F. hepatica. Our results will be helpful to investigations into the molecular taxonomy, population genetics, and ecology of F. gigantica, F. hepatica, and “intermediate” Fasciola forms.     

Characterisation of Pseudosuccinea columella and Radix natalensis (Gastropoda: Lymnaeidae) in Egypt using shell and molecular data

ABSTRACT

Pseudosuccinea columella and Radix natalensis live in the same habitat in Egypt and are important intermediate hosts of Fasciola hepatica and F. gigantica. Our study aimed to characterise both snail species using molecular analysis and shell measurements. The ranges of morphometric parameters overlapped in the two lymnaeids, indicating that they do not clearly differentiate the two species. PCR-sequence analysis of the nuclear ribosomal small subunit rRNA and the polymorphic mitochondrial cytochrome oxidase subunit 1 (CO1) genes were used to determine the genetic identity and the potential diversity of the snails. Little intrasequence variations were detected in the sequences of both gene loci, indicating the potential homogeneity of lymnaeid populations in Egypt. Generated sequences of the mitochondrial CO1 gene locus for R. natalensis showed obvious heterogeneity compared to other sequences in GenBank. Molecular characterisation of these lymnaeids might help to understand the snails’ biodiversity in a bid to control these populations and their related diseases.     

Molecular characterization of Fasciola spp. from the endemic area of northern Iran based on nuclear ribosomal DNA sequences

Fasciolosis caused by Fasciola spp. (Platyhelminthes: Trematoda: Digenea) is considered as the most important helminth infection of ruminants in tropical countries, causing considerable socioeconomic problems. In the endemic regions of the North of Iran, Fasciola hepatica and Fasciola gigantica have been previously characterized on the basis of morphometric differences, but the use of molecular markers is necessary to distinguish exactly between species and intermediate forms. Samples from buffaloes and goats from different localities of northern Iran were identified morphologically and then genetically characterized by sequences of the first (ITS-1) and second (ITS-2) Internal Transcribed Spacers (ITS) of nuclear ribosomal DNA (rDNA). Comparison of the ITS of the northern Iranian samples with sequences of Fasciola spp. from GenBank showed that the examined specimens had sequences identical to those of the most frequent haplotypes of F. hepatica (n = 25, 48.1%) and F. gigantica (n = 20, 38.45%), which differed from each other in different variable nucleotide positions of ITS region sequences, and their intermediate forms (n = 7, 13.45%), which had nucleotides overlapped between the two Fasciola species in all the positions. The ITS sequences from populations of Fasciola isolates in buffaloes and goats had experienced introgression/hybridization as previously reported in isolates from other ruminants and humans. Based on ITS-1 and ITS-2 sequences, flukes are scattered in pure F. hepatica, F. gigantica and intermediate Fasciola clades, revealing that multiple genotypes of Fasciola are able to infect goats and buffaloes in North of Iran. Furthermore, the phylogenetic trees based upon the ITS-1 and ITS-2 sequences showed a close relationship of the Iranian samples with isolates of F. hepatica and F. gigantica from different localities of Africa and Asia. In the present study, the intergenic transcribed spacers ITS-1 and ITS-2 showed to be reliable approaches for the genetic differentiation of Fasciola spp., providing bases for further studies on F. hepatica, F. gigantica and their intermediate forms in the endemic areas in Asia.     

Sequence heterogeneity in a repetitive dna element of Fasciola

A repetitive DNA sequence used as a specific probe for Fasciola hepatica infections in snails was examined in F. hepatica and Fasciola gigantica (36 individuals) specimens from five continents. The degree of intraspecific identity ranged from 79 to 99% in F. hepatica and from 93 to 99% in F. gigantica. The interspecific identity ranged from 81 to 100%, confirming the suggestion that the DNA probe sequence could be used worldwide as an epidemiological tool for the examination of the intermediate host snail in fasciolosis. Differentiation between F. hepatica and F. gigantica infections in snails was not possible using the probe.     

Molecular characterization of Fasciola gigantica from Mauritania based on mitochondrial and nuclear ribosomal DNA sequences

Fasciolosis caused by Fasciola hepatica and Fasciola gigantica (Platyhelminthes: Trematoda: Digenea) is considered the most important helminth infection of ruminants in tropical countries, causing considerable socioeconomic problems. From Africa, F. gigantica has been previously characterized from Burkina Faso, Senegal, Kenya, Zambia and Mali, while F. hepatica has been reported from Morocco and Tunisia, and both species have been observed from Ethiopia and Egypt on the basis of morphometric differences, while the use of molecular markers is necessary to distinguish exactly between species. Samples identified morphologically as F. gigantica (n = 60) from sheep and cattle from different geographical localities of Mauritania were genetically characterized by sequences of the first (ITS-1), the 5.8S, and second (ITS-2) Internal Transcribed Spacers (ITS) of nuclear ribosomal DNA (rDNA) genes and the mitochondrial Cytochrome c Oxidase I (COI) gene.Comparison of the sequences of the Mauritanian samples with sequences of Fasciola spp. from GenBank confirmed that all samples belong to the species F. gigantica. The nucleotide sequencing of ITS rDNA of F. gigantica showed no nucleotide variation in the ITS-1, 5.8S, and ITS-2 rDNA sequences among all samples examined and those from Burkina Faso, Kenya, Egypt and Iran. The phylogenetic trees based on the ITS-1 and ITS-2 sequences showed a close relationship of the Mauritanian samples with isolates of F. gigantica from different localities of Africa and Asia. The COI genotypes of the Mauritanian specimens of F. gigantica had a high level of diversity, and they belonged to the F. gigantica phylogenically distinguishable clade. The present study is the first molecular characterization of F. gigantica in sheep and cattle from Mauritania, allowing a reliable approach for the genetic differentiation of Fasciola spp. and providing basis for further studies on liver flukes in the African countries.     

Similar yet different: co-analysis of the genetic diversity and structure of an invasive nematode parasite and its invasive mammalian host

Animal parasitic nematodes can cause serious diseases and their emergence in new areas can be an issue of major concern for biodiversity conservation and human health. Their ability to adapt to new environments and hosts is likely to be affected by their degree of genetic diversity, with gene flow between distinct populations counteracting genetic drift and increasing effective population size. The raccoon roundworm (Baylisascaris procyonis), a gastrointestinal parasite of the raccoon (Procyon lotor), has increased its global geographic range after being translocated with its host. The raccoon has been introduced multiple times to Germany, but not all its populations are infected with the parasite. While fewer introduced individuals may have led to reduced diversity in the parasite, admixture between different founder populations may have counteracted genetic drift and bottlenecks. Here, we analyse the population genetic structure of the roundworm and its raccoon host at the intersection of distinct raccoon populations infected with B. procyonis. We found evidence for two parasite clusters resulting from independent introductions. Both clusters exhibited an extremely low genetic diversity, suggesting small founding populations subjected to inbreeding and genetic drift with no, or very limited, genetic influx from population admixture. Comparison of the population genetic structures of both host and parasite suggested that the parasite spread to an uninfected raccoon founder population. On the other hand, an almost perfect match between cluster boundaries also suggested that the population genetic structure of B. procyonis has remained stable since its introduction, mirroring that of its raccoon host.     

Genotypic characterization and species identification of Fasciola spp. with implications regarding the isolates infecting goats in Vietnam

Ribosomal RNA sequences (361 or 362 bp) of the second internal transcribed spacer 2 (ITS-2) and a portion of mitochondrial cox1 (423 bp) for Fasciola spp. obtained from specimens collected in indigenous and hybrid goats and sheep in Vietnam were characterized for genotypic status and hybridization/introgression. Alignment of 48 ITS-2 sequences (also those from goats and sheep in this study) indicates that F. gigantica and F. hepatica differ typically from each other at seven sites whereas one of these is a distinguishing deletion (T) at the 327th position in F. gigantica relative to F. hepatica. The isolates from the mountainous goats in the North of Vietnam (Yen Bai province) showed the ITS-2 composition relatively identical to that of F. hepatica. The ITS-2 sequences from populations of Fasciola isolates in goats had probably experienced introgression/hybridization as reported previously in other ruminants and humans. All Vietnamese goat-of-origin specimens had high pairwise percentage of mitochondrial cox1 sequences to F. gigantica (97-100%), and very low identity to F. hepatica (91-93%), suggesting their maternal linkage to be traced to F. gigantica. The presence of hybrid and/or introgressed populations of liver flukes bearing genetic material from both F. hepatica and F. gigantica in the goats/sheep in Vietnam, regardless of indigenous or imported hosts, appears to be the first demonstration from a tropical country.     

Immunodiagnosis of Fasciola gigantica Infection Using Monoclonal Antibody-Based Sandwich ELISA and Immunochromatographic Assay for Detection of Circulating Cathepsin L1 Protease

Background

Tropical fasciolosis caused by Fasciola gigantica infection is one of the major diseases infecting ruminants in the tropical regions of Africa and Asia including Thailand. Parasitological diagnosis of fasciolosis is often unreliable and possesses low sensitivity. Therefore, the detection of circulating parasite antigens is thought to be a better alternative for diagnosis of fasciolosis, as it reflects the real parasite burden.

Methods

In this study, we have produced a monoclonal antibody (MoAb) against recombinant F. gigantica cathepsin L1 (rFgCatL1), and developed both sandwich enzyme-linked immunosorbent assay (sandwich ELISA) and immunochromatographic (IC) test for rapid detection of circulating cathepsin L1 protease (CatL1) in the sera from mice experimentally and cattle naturally infected with Fasciola gigantica. MoAb 4E3 and biotinylated rabbit anti-recombinant CatL1 antibody were selected due to their high reactivities and specificities.

Results

The lower detection limits of sandwich ELISA and IC test were 3 pg/ml and 0.256 ng/ml, respectively. Sandwich ELISA and IC test could detect F. gigantica infection from day 1 to 35 post infection. In experimental mice, the sensitivity, specificity and accuracy were 95%, 100% and 98.6% (for sandwich ELISA), and 93%, 100% and 98.2% (for IC test), while in natural cattle they were 98.3%, 100% and 99.5% (for sandwich ELISA), and 96.7%, 100% and 99.1% (for IC test).

Conclusions

These two assay methods showed high efficiencies and precisions for diagnosis of fasciolosis by F. gigantica.     

Occurrence of spermic diploid and aspermic triploid forms of Fasciola in Vietnam and their molecular characterization based on nuclear and mitochondrial DNA

Fasciola spp. found in Asian countries are diversified in nature, and they should therefore be characterized by spermatogenesis, ploidy and genetic differentiation as well as morphology. The present study showed that spermic diploid and aspermic triploid forms of Fasciola occurred in Vietnam. The spermic diploid specimens were accurately identified as F. gigantica, while the aspermic triploids could not be identified on the basis of their morphology by the ratio of body length and width and DNA sequences of nuclear ribosomal ITS1 and mitochondrial NDI and COI genes. The molecular data also indicated that Vietnamese aspermic triploids might be hybrids and/or their offspring between Fasciola hepatica and F. gigantica, because they showed the ITS1-Fh/Fg haplotype, which had chimeric sequences of the two species. Furthermore, the aspermic triploids seem to have originated in countries other than Vietnam and to have rapidly spread to that country with infected animals.     

Towards delineating functions within the fasciola secreted cathepsin l protease family by integrating in vivo based sub-proteomics and phylogenetics

Background

Fasciola hepatica, along with Fasciola gigantica, is the causative agent of fasciolosis, a foodborne zoonotic disease affecting grazing animals and humans worldwide. Pathology is directly related to the release of parasite proteins that facilitate establishment within the host. The dominant components of these excretory-secretory (ES) products are also the most promising vaccine candidates, the cathepsin L (Cat L) protease family.

Methodology/Principal Findings

The sub-proteome of Cat L proteases from adult F. hepatica ES products derived from in vitro culture and in vivo from ovine host bile were compared by 2-DE. The individual Cat L proteases were identified by tandem mass spectrometry with the support of an in-house translated liver fluke EST database. The study reveals plasticity within the CL1 clade of Cat L proteases; highlighted by the identification of a novel isoform and CL1 sub-clade, resulting in a new Cat L phylogenetic analysis including representatives from other adult Cat L phylogenetic clades. Additionally, for the first time, mass spectrometry was shown to be sufficiently sensitive to reveal single amino acid polymorphisms in a resolved 2-DE protein spot derived from pooled population samples.

Conclusions/Significance

We have investigated the sub-proteome at the population level of a vaccine target family using the Cat L proteases from F. hepatica as a case study. We have confirmed that F. hepatica exhibits more plasticity in the expression of the secreted CL1 clade of Cat L proteases at the protein level than previously realised. We recommend that superfamily based vaccine discovery programmes should screen parasite populations from different host populations and, if required, different host species via sub-proteomic assay in order to confirm the relative expression at the protein level prior to the vaccine development phase.