Molecular characterization and phylogenetic analysis of Fasciola gigantica from western Java, Indonesia

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.

Prevalence, risk factors, and infection intensity of fasciolosis in dairy cattle in Boyolali, Indonesia

Background and Aim: Fasciolosis is a significant problem in veterinary and public health, causing huge economic losses. Epidemiological studies of fasciolosis in dairy cattle in Indonesia are few and existing reports primarily focus on prevalence. This study aimed to determine the prevalence, risk factors, and infection intensity of fasciolosis in dairy cattle in Boyolali, Indonesia. Materials and Methods: This cross-sectional study included 400 dairy cattle from 72 household farms in eight subdistricts. Fecal samples (n=400) were examined using the Flukefinder® kit and the simple sedimentation technique was the gold standard for fasciolosis. In-person interviews using questionnaires collected data on farmers, farms, and animal characteristics. Chi-square and logistic regression analyses were performed to evaluate the associated risk factors for fasciolosis, and p < 0.05 was considered statistically significant. Results: The overall prevalence of fasciolosis in dairy cattle in Boyolali, Indonesia, was 16.50% (95% confidence interval [CI] 12.85-20.15) at the animal level (n = 400), whereas 40.28% at household farms (n = 72) level (95% CI 18.67-51.88). The relative sensitivity and specificity of the Flukefinder® kit compared with those of the gold standard were 79.49% and 92.52%, respectively, with a moderate agreement (kappa=0.59; p < 0.001). Fasciolosis was more likely in cattle originating from the Mojosongo subdistrict than from other subdistricts (odds ratio (OR)=5.28, 95% CI 1.22-22.94); from farms that did not process manure versus from those that did (OR = 3.03, 95% CI 1.43-4.71); and with farmers that had never attended extension programs compared with those who had (OR = 4.72, 95% CI 1.99-11.19). Studied cattle were mostly affected by light Fasciola spp. infections (92.4%, 95% CI 77.8-100%) followed by moderate (6.1%, 95% CI 0-22.2%) and heavy (1.5%, 95% CI 0-5.6%) infections. Conclusion: Fasciolosis is prevalent in dairy cattle in Boyolali, Indonesia. Control efforts should target the high-risk Mojosongo subdistrict, emphasize the importance of processing manure, and encourage farmers to attend extension programs. Flukefinder® is a practical on-site diagnostic kit for fasciolosis in Indonesian dairy farms. Parasite species identification and a malacological survey of intermediate hosts of Fasciola spp. in the farming environment are required for further research.

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.

Towards the comprehension of fasciolosis (re-)emergence: an integrative overview

The increasing distribution and prevalence of fasciolosis in both human and livestock are concerning. Here, we examine the various types of factors influencing fasciolosis transmission and burden and the interrelations that may exist between them. We present the arsenal of molecules, 'adjusting' capabilities and parasitic strategies of Fasciola to infect. Such features define the high adaptability of Fasciola species for parasitism that facilitate their transmission. We discuss current environmental perturbations (increase of livestock and land use, climate change, introduction of alien species and biodiversity loss) in relation to fasciolosis dynamics. As Fasciola infection is directly and ultimately linked to livestock management, living conditions and cultural habits, which are also changing under the pressure of globalization and climate change, the social component of transmission is also discussed. Lastly, we examine the implication of increasing scientific and political awareness in highlighting the current circulation of fasciolosis and boosting epidemiological surveys and novel diagnostic techniques. From a joint perspective, it becomes clear that factors weight differently at each place and moment, depending on the biological, environmental, social and political interrelating contexts. Therefore, the analyses of a disease as complex as fasciolosis should be as integrative as possible to dissect the realities featuring each epidemiological scenario. Such a comprehensive appraisal is presented in this review and constitutes its main asset to serve as a fresh integrative understanding 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.

Prevalence and risk factors of trematode infection in swamp buffaloes reared under different agro-climatic conditions in Java Island of Indonesia

Background and Aims:

This study was conducted to estimate the prevalence and intensity and to identify the associated risk factors and impact of trematode infection in swamp buffaloes reared under different agro-climatic conditions in Java, Indonesia.

Materials and Methods:

A total of 580 fecal samples were collected from swamp buffaloes in five different agro-climatic areas in Banten and Central Java Provinces, Indonesia. The fecal samples were examined using the Danish Bilharziasis Laboratory Technique to determine the prevalence and intensity of trematode infection. The risk factors for infection were determined from an in-depth interview of owners/keepers, and the results were analyzed using Chi-square tests and multiple logistic regression. The infection was also correlated with swamp buffalo production parameters (body weight and body condition score [BCS]).

Results:

From all fecal samples, the overall prevalence of trematode infection was 64.83%, which comprised Fasciola spp. (16.03%; mean eggs per gram [EPG]±SD: 1.02±0.43) and Paramphistomatidae (62.93%; mean EPG±SD: 1.01±0.66). The main risk factor for trematode infection was feeding animals with rice straw (odds ratio [OR]: 40.124); the risk of trematode infection was 40.142 times higher in buffaloes that consumed rice straw. Other risk factors included the frequency of anthelmintic treatment (OR: 4.666), age (OR: 0.449), and drinking water source (OR: 0.358). Trematode infection did not significantly affect the body weight or BCS of swamp buffaloes.

Conclusion:

Although the prevalence of trematode infection was high in swamp buffaloes, the intensity of infection was low, and the infection did not affect the animals’ physical parameters.

Genetic diversity and multiplicity of infection in Fasciola gigantica isolates of Pakistani livestock

Fasciola spp. are responsible for over 3 billion US dollars of production loss annually in livestock and cause widespread zoonotic disease. Nevertheless, understating of the emergence and spread of the trematode species is poor. The multiplicity of F. gigantica infection and its spread is potentially influenced by multiple factors, including the abundance of suitable intermediate hosts, climatic conditions favouring the completion of the parasite's lifecycle, and translocation of infected animals, or free-living parasite stages between regions. Here we describe the development of a 'tremabiome' metabarcoding sequencing method to explore the numbers of F. gigantica genotypes per infection and patterns of parasite spread, based on genetic characteristics of the mitochondrial NADH dehydrogenase 1 (mt-ND-1) locus. We collected F. gigantica from three abattoirs in the Punjab and Balochistan provinces of Pakistan, and our results show a high level of genetic diversity in 20 F. gigantica populations derived from small and large ruminants consigned to slaughter in both provinces. This implies that F. gigantica can reproduce in its definitive hosts through meiosis involving cross- and self-breeding, as described in the closely related species, Fasciola hepatica. The genetic diversity between the 20 populations derived from different locations also illustrates the impact of animal movements on gene flow. Our results demonstrate the predominance of single haplotypes, consistent with a single introduction of F. gigantica infection in 85% of the hosts from which the parasite populations were derived. This is consistent with clonal reproduction in the intermediate snail hosts.

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 Characterization of Cambodian Fasciola gigantica and Dispersal Direction of the Species in Asia

Fasciola gigantica and hybrid Fasciola are distributed throughout Asia. Herein, we investigated the species of the Fasciola fluke distributed in three hotspots of fascioliasis in Cambodia. A total of 92 flukes collected from 21 slaughtered cattle from Kandal (44), Battambang (41), and Kratie (7) Provinces were identified as F. gigantica using multiplex PCR for a nuclear phosphoenolpyruvate carboxykinase (PEPCK) gene. The overall prevalence of F. gigantica infestation was 7.14% (21/294). Phylogenetic as well as population genetics analyses were performed using the mitochondrial NADH dehydrogenase subunit 1 (ND1). The 19 ND1 haplotypes were identified from Cambodian F. gigantica (haplotype diversity, 0.83). All of the haplotypes were classified into F. gigantica haplogroup C, which includes ND1 haplotypes detected from Thailand, Vietnam, Indonesia, Myanmar, and China. Among haplogroup C, novel and unique haplotypes of Cambodia were found in the Battambang and Kandal Provinces, and the nucleotide diversity of the Cambodian population (0.00532) was the highest. Pairwise fixation indices among the F. gigantica populations from these countries indicated that the Cambodian and Thailand populations were related to each other. The highest genetic diversity in the Cambodian population suggests that F. gigantica in Cambodia may be the ancestor of the populations in Southeast Asian countries. Most likely, livestock movement, including Zebu cattle, played an important role in the transmission of F. gigantica. In this study, the hybrid Fasciola flukes that are commonly found in neighboring countries, were not found in Cambodia. Further comprehensive investigations of Fasciola prevalence should be conducted by analyzing a wider range of hosts throughout Cambodia to reach a more solid conclusion about the absence of hybrid flukes.

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

Recently, phosphoenol pyruvate carboxykinase (pepck) and DNA polymerase delta (pold) were established as reliable nuclear markers for species identification of Fasciola spp. in multiplex polymerase chain reaction (PCR) and PCR-restriction fragment length polymorphism-based assays, respectively. Currently, little is known about Fasciola species distribution in Central Asia. Therefore, the objective of the present study was to perform precise molecular species identification of liver flukes from Afghanistan and to reveal their dispersal route(s) via phylogenetic analysis based on mitochondrial nad1 haplotypes. Ninety-two Fasciolaflukes collected from sheep in Kabul, Afghanistan, were identified as F. hepatica based on pepck and pold screening. Although the pepck fragment pattern obtained via multiplex PCR analysis could not distinguish the species of the seven Fasciola flukes, the pepck nucleotide sequence data confirmed that they were F. hepatica.The 20 nad1 haplotypes detected among the Afghani liver flukes were closely related to those from China and Egypt, with the F_STvalue (-0.003, P = .41) between the F. hepatica populations from Afghanistan and China confirming a very close relationship. Nucleotide diversity was greater in the population from Afghanistan compared with that from China, indicating that the Afghani population was older, and that the dispersal direction of F. hepatica was from Afghanistan to China. The results of the present study contribute to our understanding of the dispersal of F. hepatica from its predicted origin, the Fertile Crescent.

Distribution Status of Hybrid Types in Large Liver Flukes, Fasciola Species (Digenea: Fasciolidae), from Ruminants and Humans in Vietnam

The aim of this study is to delineate ‘admixed hybrid’ and ‘introgressive’ Fasciola genotypes present in the Fasciola population in Vietnam. Adult liver flukes collected from ruminants in 18 Provinces were morphologically sorted out by naked eyes for small (S), medium (M) and large (L) body shapes; and human samples (n=14) from patients. Nuclear ribosomal (rDNA) ITS1 and ITS2, and mitochondrial (mtDNA) nad1 markers were used for determination of their genetic status. Total 4,725 worm samples of ruminants were tentatively classified by their size: 6% (n=284) small (S)-, 13% (n=614) medium (M)-, and 81% (n=3,827) large (L)-forms. All the representative (n=120, as 40 each group) and 14 human specimens, possessed maternal mtDNA of only F. gigantica and none of F. hepatica. Paternally, all (100%) of the L-(n=40) and 77.5% (n=31) of the M-flukes had single F. gigantica rDNA indicating ‘pure’ F. gigantica. A majority (90%, n=36) of the S- and 15% (n=6) of the M-worms had single F. hepatica rDNA, indicating their introgressive; the rest (10%, n=4) of the S- and 7.5% (n=3) of the M-flukes had mixture of both F. gigantica and F. hepatica rDNAs, confirming their admixed hybrid genetic status. Fourteen human samples revealed 9 (64%) of pure F. gigantica, 3 (22%) of introgressive and 2 (14%) of admixed hybrid Fasciola spp. By the present study, it was confirmed that the small worms, which are morphologically identical with F. hepatica, are admixed and/or introgressive hybrids of Fasciola spp., and able to be the pathogens of human fascioliasis.

First genetic characterization of Fasciola hepatica in Argentina by nuclear and mitochondrial gene markers

Fasciola hepatica is a trematode showing genetic variation among isolates from different regions of the world. The objective of this work was to characterize for the first time F. hepatica isolates circulating in different regions of Argentina. Twenty-two adult flukes were collected from naturally infected bovine livers in different areas from Argentina and used for DNA extraction. We carried out PCR amplification and sequence analysis of the ribosomal internal transcribed spacer 1 (ITS1), mitochondrial nicotinamide adenine dinucleotide dehydrogenase subunits 4 and 5 (nad4 and nad5) and mitochondrial cytochrome c oxidase subunit I (cox1) genes as genetic markers. Phylogenies were reconstructed using maximum parsimony algorithm. A total of 6 haplotypes were found for cox1, 4 haplotypes for nad4 and 3 haplotypes for nad5. The sequenced ITS1 fragment was identical in all samples. The analyzed cox1 gene fragment is the most variable marker and is recommended for future analyses. No geographic association was found in the Argentinean samples.

Identity of Fasciola spp. in sheep in Egypt

Background

In Egypt, liver flukes, Fasciola spp. (Digenea: Fasciolidae), have a serious impact on the farming industry and public health. Both Fasciola hepatica and Fasciola gigantica are known to occur in cattle, providing the opportunity for genetic recombination. Little is known on the identity and genetic variability of Fasciola populations in sheep.

Methods

This study was performed to determine the prevalence of liver flukes in sheep in Menofia Province as a representative area of the delta region in Egypt, as measured by postmortem examination of slaughtered animals at three abattoirs. The identity and genetic variability of Fasciola spp. in slaughtered animals were determined by PCR-sequence analysis of the nuclear ribosomal internal transcribed spacer 1 (ITS1) and the mitochondrial NADH dehydrogenase subunit 1 (nad1) genes.

Results

Physical inspection of the liver indicated that 302 of 2058 (14.7%) slaughtered sheep were infected with Fasciola spp. Sequence analysis of the ITS1 and nad1 genes of liver flukes from 17 animals revealed that 11 animals were infected with F. hepatica, four with F. gigantica, and two with both species. Seventy eight of 103 flukes genetically characterized from these animals were F. hepatica, 23 were F. gigantica, and two had ITS1 sequences identical to F. hepatica but nad1 sequences identical to F. gigantica. nad1 sequences of Egyptian isolates of F. gigantica showed pronounced differences from those in the GenBank database. Egyptian F. gigantica haplotypes formed haplogroup D, which clustered in a sister clade with haplogroups A, B and C circulating in Asia, indicating the existence of geographic isolation in the species.

Conclusions

Both F. hepatica and F. gigantica are prevalent in sheep in Egypt and an introgressed form of the two occurs as the result of genetic recombination. In addition, a geographically isolated F. gigantica population is present in the country. The importance of these observations in epidemiology of fascioliasis needs to be examined in future studies.