Molecular analysis reveals expansion of Fasciola hepatica distribution from Afghanistan to China

Genetic variation and population structure of Fasciola hepatica: an in silico analysis

Fasciola hepatica is a trematode leading to heavy economic setbacks to the livestock sector globally. The population's genetic information and intimate kinship level are frequently assessed using analysis of mitochondrial DNA. In this analysis, we retrieved cox1 (n = 247) and nad1 (n = 357) sequences of F. hepatica from the NCBI GenBank database and aligned the sequences with the respective reference sequences using MEGA software. The median joining network was drawn using PopArt software while neutrality and diversity indices were estimated with the help of DnaSp software. Neighbor-joining phylogenetic tree was constructed using the MEGA software package. A total of 46 and 98 distinctive haplotypes were observed for cox1 and nad1 genes, respectively. Diversity indices indicated high haplotype and nucleotide diversities in both genes. Positive Tajima's D and Fu's Fs values were found for the entire population of both the genes under study. The cox1 and nad1 gene segments in this study showed high Tajima's D values, suggesting a low likelihood of future population growth. The Tajima's D value of the nad1 gene sequence is lower (2.14910) than that of the cox1 gene sequence (3.40314), which suggests that the former is growing at a slower rate. However, the region-wise analysis revealed that both the cox1 and nad1 genes showed deviation from neutrality suggesting a recent population expansion as a result of an excess of low-frequency polymorphism. Furthermore, the overall host-wise analysis showed positive and significant Tajima's D values for the cox1 and nad1 gene sequences. To the best of our knowledge, this is the first attempt to provide insights into genetic variations and population structure of F. hepatica at a global scale using cox1 and nad1 genes. Our findings suggest the existence of specific variants of F. hepatica in different parts of the world and provide information on the molecular ecology of F. hepatica. The results of this study also mark a critical development in upcoming epidemiological investigations on F. hepatica and will also contribute to understanding the global molecular epidemiology and population structure of F. hepatica.

Human and Animal Fascioliasis: Origins and Worldwide Evolving Scenario

Fascioliasis is a plant- and waterborne zoonotic parasitic disease caused by two trematode species: (i) Fasciola hepatica in Europe, Asia, Africa, the Americas, and Oceania and (ii) F. gigantica, which is restricted to Africa and Asia. Fasciolid liver flukes infect mainly herbivores as ruminants, equids, and camelids but also omnivore mammals as humans and swine and are transmitted by freshwater Lymnaeidae snail vectors. Two phases may be distinguished in fasciolid evolution. The long predomestication period includes the F. gigantica origin in east-southern Africa around the mid-Miocene, the F. hepatica origin in the Near-Middle East of Asia around the latest Miocene to Early Pliocene, and their subsequent local spread. The short postdomestication period includes the worldwide spread by human-guided movements of animals in the last 12,000 years and the more recent transoceanic anthropogenic introductions of F. hepatica into the Americas and Oceania and of F. gigantica into several large islands of the Pacific with ships transporting livestock in the last 500 years. The routes and chronology of the spreading waves followed by both fasciolids into the five continents are redefined on the basis of recently generated knowledge of human-guided movements of domesticated hosts. No local, zonal, or regional situation showing disagreement with historical records was found, although in a few world zones the available knowledge is still insufficient. The anthropogenically accelerated evolution of fasciolids allows us to call them "peridomestic endoparasites." The multidisciplinary implications for crucial aspects of the disease should therefore lead the present baseline update to be taken into account in future research studies.

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.

Multiplex PCR and Sequence Analysis to Investigate Genetic Diversity of Fasciola Isolates from Cattle and Sheep in Turkey

Fasciolosis is a highly prevalent helminthic infection caused by Fasciola hepatica and F. gigantica. With the aim of identifying hybrid Fasciola flukes, multiplex PCR was performed to amplify the pepck gene. Furthermore, to determine Fasciola haplotypes, mitochondrial nad1 gene was amplified and sequenced, and phylogenetic analyses were performed. Adult Fasciola isolates were collected from 51 cattle and 51 sheep, genomic DNA was isolated, and one-step multiplex PCR was subsequently performed to amplify pepck. Isolates that generated a 510 bp band were identified as F. gigantica, those that generated a 241 bp band were identified as F. hepatica, and those that generated both bands were identified as hybrid (aspermic) flukes. Multiplex PCR data identified four isolates as F. gigantica and 84 as F. hepatica. Fourteen hybrid isolates (five cattle and nine sheep) were identified. On unidirectional DNA sequence analysis of nad1 PCR products, three sequences were identified as F. gigantica and 99 as F. hepatica. In addition, only 4 of 87 haplotypes detected for F. hepatica nad1 sequences were found to be previously reported, while the remaining 83 are unique to this study. To summarize, this study is the first to report the existence of hybrid Fasciola flukes and 83 unique haplotypes of F. hepatica in Turkey.

Development of novel DNA marker for species discrimination of Fasciola flukes based on the fatty acid binding protein type I gene

Background

Multiplex polymerase chain reaction (PCR) and PCR-restriction fragment length polymorphism (RFLP) for nuclear phosphoenolpyruvate carboxykinase (pepck) and polymerase delta (pold), respectively, have been used to differentiate Fasciola hepatica, F. gigantica, and hybrid Fasciola flukes. However, discrimination errors have been reported in both methods. This study aimed to develop a multiplex PCR based on a novel nuclear marker, the fatty acid binding protein type I (FABP) type I gene.

Methods

Nucleotide sequence variations of FABP type I were analyzed using DNA samples of F. hepatica, F. gigantica, and hybrid Fasciola flukes obtained from 11 countries in Europe, Latin America, Africa, and Asia. A common forward primer for F. hepatica and F. gigantica and two specific reverse primers for F. hepatica and F. gigantica were designed for multiplex PCR.

Results

Specific fragments of F. hepatica (290 bp) and F. gigantica (190 bp) were successfully amplified using multiplex PCR. However, the hybrid flukes contained fragments of both species. The multiplex PCR for FABP type I could precisely discriminate the 1312 Fasciola samples used in this study. Notably, no discrimination errors were observed with this novel method.

Conclusions

Multiplex PCR for FABP type I can be used as a species discrimination marker in place of pepck and pold. The robustness of the species-specific primer should be continuously examined using a larger number of Fasciola flukes worldwide in the future since nucleotide substitutions in the primer regions may cause amplification errors.

Graphical abstract

Supplementary Information

The online version contains supplementary material available at 10.1186/s13071-022-05538-7.

Population structure, molecular characterization, and phylogenetic analysis of Fasciola gigantica from two locations in Uganda

Fasciola gigantica is a major pathogen that causes fasciolosis in Africa. A recent study in Uganda demonstrated that Fasciola flukes were present in 65.7% of slaughtered cattle. However, molecular identification of Fasciola species has not yet been performed in the country. In the present study, 292 Fasciola flukes were collected from Kampala and Gulu, Uganda. The samples were identified as F. gigantica using a multiplex polymerase chain reaction (PCR) assay for phosphoenolpyruvate carboxykinase (pepck) and a PCR-restriction fragment length polymorphism (RFLP) assay for DNA polymerase delta (pold). A significant genetic difference between F. gigantica obtained from cattle slaughtered at Kampala and Gulu was observed by analyzing the mitochondrial markers NADH dehydrogenase subunit 1 (nad1) and cytochrome C oxidase subunit 1 (cox1). Fasciola collected from Gulu had a more diversified population than that collected from Kampala, probably because of differences in livestock management systems. One of the possible reasons for this observation is that cattle slaughtered in Gulu were reared under an extensive communal grazing system, which is suitable for maintaining parasite diversity, whereas cattle slaughtered in Kampala mainly originated from fenced/closed farms, which limits parasite diversity. However, the cause of the difference between these two locations was not clearly defined by the results of this study. The F. gigantica population from Uganda was related to that obtained from Zambia. A star-like phylogeny was detected in a median-joining network analysis, which indicated rapid population expansion and suggested that the F. gigantica populations from both countries are maintained by domestic ruminants in eastern Africa. Interestingly, the F. gigantica population from Uganda was not related to those from Egypt and Nigeria. The results of the present study suggest that F. gigantica populations in African countries are indigenous to each country or region.

Proteomic Analysis of Extracellular Vesicles From Fasciola hepatica Hatching Eggs and Juveniles in Culture

The identification of extracellular vesicles (EVs) in Fasciola hepatica has provided a new way to understand parasite-host communication. Most of the studies on EVs have focused on the adult stage of F. hepatica, but recently, the presence of EVs from different developmental stages has been reported. To better understand the potential role of EVs in the biology of the parasite and in the infection process, the protein cargo of EVs from embryonated eggs and newly-excysted juvenile (NEJs) flukes cultured up to 28 days, has been analyzed. EVs were isolated by size exclusion chromatography and evaluated by nanoparticle tracking analysis and transmission electron microscopy. LC-MS/MS proteomic analysis of EVs revealed the presence of 23 different proteins from embryonated egg-derived EVs and 29 different proteins from NEJ-derived EVs. Most of the identified proteins had been previously described in EVs from F. hepatica adults, including cytoskeletal proteins, glycolytic enzymes, stress-related proteins and tetraspanins. Nevertheless, EVs from hatching eggs and NEJs exhibited qualitative differences in composition, when compared to EVs form adults, including the absence of cathepsin cysteine peptidases. The differential content of the EVs released by the different developmental stages of the parasite reflect the intense activity of NEJs at this early stage, with several proteins involved in membrane traffic and cell physiology. This new set of identified proteins could help to understand key metabolic, biochemical and molecular mechanisms mediated by EVs that take place upon egg hatching and after parasite excystment.

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.

On the arrival of fasciolosis in the Americas

Fasciola hepatica is a worldwide emerging and re-emerging parasite heavily affecting several regions in South America. Some lymnaeid snail species of American origin are among the major hosts of F. hepatica worldwide. Recent paleoparasitological findings detected its DNA in a 2300-year-old sample in Patagonia, countering the common hypothesis of the recent arrival of F. hepatica in the Americas during European colonization. Thus, the theory of an initial introduction in the 1500s can no longer be sustained. This article discusses how it was possible for F. hepatica to reach and spread in the Americas in relation to the availability and compatibility of hosts through natural and incidental introductions. Our study will serve to better understand the ongoing Neotropical scenario of fasciolosis.

Molecular characterization and phylogenetic analyses of Fasciola gigantica of buffaloes and goats in Punjab, Pakistan

Fasciola gigantica is considered to be a major pathogen causing fasciolosis in the Indian subcontinent, resulting in production losses of millions of dollars in the livestock industry. Understading the dispersal origin and the patterns of spread of F. gigantica is important. A total of 53 Fasciola flukes collected from buffaloes and goats in Punjab, Pakistan between 2017 and 2018 were identified as F. gigantica based on the multiplex PCR for the phosphoenolpyruvate carboxykinase (pepck) and the PCR-restriction fragment length polymorphism (RFLP) for DNA polymerase delta (pold). A significant genetic difference between F. gigantica from buffaloes and goats was indicated by the genetic analyses of mitochondrial markers, NADH dehydrogenase subunit 1 (nad1) and cytochrome C oxidase subunit 1 (cox1). Phylogenetic analysis of the seventeen nad1 haplotypes of F. gigantica from Pakistan with those in neighbouring countries of the Indian subcontinent revealed that all the haplotypes identified in Pakistan were clustered in haplogroup A. fasciola gigantica with the eight haplotypes might be expanded in Pakistan from Indian origin, along with the migration of the domestic animals, since they were related to Indian haplotypes. In contrast, the remaining nine haplotypes were not shared with any neighbouring countries, suggesting independent origin, probably from neighbouring Middle East countries. However, cautious interpretation is required due to the very limited samples size of this study. Our study provides a proof of concept for a method that could be used to investigate the epidemiology of F. gigantica.

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.

Molecular characterization and phylogenetic analysis of Fasciola hepatica from high- plateau and steppe areas in Algeria

A previous study based on mitochondrial DNA markers reported the presence of Fasciola hepatica in Algeria. However, a precise species identification is still required. In this report, a total of 68 Fasciola isolates, collected from high-plateau (Bordj-Bou-Arreridj) and steppe (Djelfa) areas of Algeria, were identified at the species level by multiplex PCR and PCR-restriction fragment length polymorphism (RFLP) for nuclear phosphoenolpyruvate carboxykinase (pepck) and DNA polymerase delta (pold), respectively. The result of the multiplex PCR conflicted with that of the PCR-RFLP; however, subsequent nucleotide sequencing of pepck clearly showed that all isolates should be classified as F. hepatica. The two mitochondrial markers, NADH dehydrogenase subunit I (nad1) and cytochrome c oxidase subunit 1 (cox1), revealed a close relationship between the parasite populations from the plateau and those from the steppe. A dispersal direction from the high plateau to the steppe was indicated because the former population was more diversified than the latter. Moreover, these populations were more closely related to populations from Spain than those from Egypt or Afghanistan. Given the population characteristic of F. hepatica in Spain and the history of cattle trade, it seems likely that the parasite was introduced to Algeria from Europe through a route across the Mediterranean Sea.

Genetic diversity of Fasciola hepatica in Spain and Peru

In the present study, molecular characterization of Fasciola flukes from Spain was performed to reveal the relation with the previously reported Peruvian F. hepatica population. The nuclear DNA markers, phosphoenolpyruvate carboxykinase (pepck) and DNA polymerase delta (pold), were used for species identification of Fasciola flukes. A total of 196 Fasciola flukes were identified as F. hepatica by pepck and pold, and 26 haplotypes were detected in mitochondrial NADH dehydrogenase subunit 1 (nad1). Only one of them was previously found in Spanish samples; which indicates the existence of high genetic diversity and population structure in F. hepatica from Spain. Three haplotypes were identical to those from Peruvian F. hepatica. The pairwise fixation index value confirmed a relatively close relationship between the Spanish and Peruvian F. hepatica samples. The Spanish samples showed clearly higher genetic variability than the Peruvian population. These results are discussed in relation with the hypothesis of the introduction of the parasite in America from Europe and recent evidence of pre-Hispanic F. hepatica from Argentina revealed by ancient DNA.

Multilocus approach reveals discordant molecular markers and corridors for gene flow between North African populations of Fasciola hepatica

Fasciolosis is a foodborne trematodosis characterised by a worldwide distribution. Various approaches have been developed for the study of the causative agents of this parasitic infection: Fasciola hepatica, Fasciola gigantica and the aspermic intermediated forms (hybrid and introgressed). In the present study, novel and common molecular markers (pepck and pold, ITS, CO1, ND1 and CO1-trnT-rrnL) were used to characterise Fasciola flukes from the Tunisian-Algerian border, to estimate the gene flow between these populations and to evaluate the reliability of different molecular markers. All nuclear and mitochondrial markers, apart from pepck, supported the monophyly of the studied flukes identified as F. hepatica. Multiplex PCR for pepck revealed three different genotypes corresponding to F. hepatica (pepck-Fh), F. gigantica (pepck-Fg) and the aspermic Fasciola flukes (pepck-Fh/Fg). Sequence analysis of pepck revealed high polymorphism, length variation, within this intronic marker. The observed inconsistencies were due to the position of the forward primer within the intronic region. Pepck sequences showed different level of heterozygosity and homozygosity with length polymorphisms in the introns. Pepck multiplex PCR patterns could not differentiate between Fasciola species. All studies based on only pepck multiplex PCR with mitochondrial markers should be revised. Nuclear and mitochondrial markers revealed an important gene flow between Tunisian and Algerian populations of F. hepatica. The combination of nuclear and mitochondrial sequence analysis is still the best method to distinguish these taxa. Effective measures are needed in order to better control cross-country illegal trade of vector.

Genetic diversity and distribution of Fasciola hepatica haplotypes in Iran: Molecular and phylogenetic studies

Fasciolosis is a zoonotic disease caused by Fasciola hepatica and Fasciola gigantica. Over the last decade, diagnostic tools to detect and differentiate Fasciola species have improved, but our knowledge of the distribution of haplotypes and gene flow of this parasite is not comprehensive yet. The purpose of this study was to investigate this gap in the epidemiology of F. hepatica in different provinces of Iran between 2015 and 2017. Isolated Fasciola were collected from abattoirs in 9 provinces. The partial sequence of mitochondrial NADH dehydrogenase subunit 1 (ND1) gene was used for the identification and molecular analysis of F. hepatica isolates. The amplified PCR products were purified and subjected to direct sequencing for subsequent construction of phylogenetic tree and network analysis. In the 130 subjects analyzed, 37 ND1 haplotypes were detected. This is the first study in Iran which investigates F. hepatica population and its genetic structure, based on mitochondrial ND1 marker in different geographical regions of Iran.