- Open Access
Genetic interrelationships of North American populations of giant liver fluke Fascioloides magna
Parasites & Vectorsvolume 8, Article number: 288 (2015)
Population structure and genetic interrelationships of giant liver fluke Fascioloides magna from all enzootic North American regions were revealed in close relation with geographical distribution of its obligate definitive cervid hosts for the first time.
Variable fragments of the mitochondrial cytochrome c oxidase subunit I (cox1; 384 bp) and nicotinamide dehydrogenase subunit I (nad1; 405 bp) were applied as a tool. The concatenated data set of both cox1 and nad1 sequences (789 bp) contained 222 sequences that resulted in 50 haplotypes. Genetic data were analysed using Bayesian Inference (BI), Maximum Likelihood (ML) and Analysis of Molecular Variance (AMOVA).
Phylogenetic analysis revealed two major clades of F. magna, which separated the parasite into western and eastern populations. Western populations included samples from Rocky Mountain trench (Alberta) and northern Pacific coast (British Columbia and Oregon), whereas, the eastern populations were represented by individuals from the Great Lakes region (Minnesota), Gulf coast, lower Mississippi, and southern Atlantic seaboard region (Mississippi, Louisiana, South Carolina, Georgia, Florida) and northern Quebec and Labrador. Haplotype network and results of AMOVA analysis confirmed explicit genetic separation of western and eastern populations of the parasite that suggests long term historical isolation of F. magna populations.
The genetic makeup of the parasite’s populations correlates with data on historical distribution of its hosts. Based on the mitochondrial data there are no signs of host specificity of F. magna adults towards any definitive host species; the detected haplotypes of giant liver fluke are shared amongst several host species in adjacent populations.
Spatial distribution of parasites is closely coupled with distribution of their hosts; parasites require suitable hosts for nutrients and other resources, and thus they cannot occur in territories where susceptible hosts are absent. Therefore, the maximum geographical range of a parasite must equal the combined ranges of its hosts. The hosts also represent the only mean of dispersal for the majority of parasitic taxa and in the case of parasites with multi-host lifecycles, the parasite’s dispersal is driven by its most mobile host. Furthermore, geographical distribution and relationships of the parasites of terrestrial hosts may be affected by historical contractions and expansions of their hosts’ distribution, including recent bottlenecks caused by humans [1, 2]. On the larger scale, a positive relationship between host and parasite species richness is inevitable and generally supported, since most host species harbour at least one host specific parasite .
Influence of hosts and parasites on their respective biodiversity is of particular interest, since each partner in a host-parasite association potentially exerts a selective pressure on the other . The relationship between host and parasite represents an intimate interaction between at least two genetic systems . The host-parasite interaction is based on subtle interplay between parasite survival strategies and host defence mechanisms . Revealing population genetic structure and host specificity of the parasitic species in question is the first natural step towards understanding the underlying processes of natural selection.
Giant liver fluke, Fascioloides magna, represents a very interesting parasitic model characterized by a wide spectrum of intermediate and definitive hosts, large spatial distribution, potential to colonize new territories and adapt to new host species. It parasitizes a wide range of definitive ruminant hosts, especially cervids. The fluke has a strong predilection to liver parenchyma where it is localized in thin-walled fibrous pseudocysts. It is generally accepted that F. magna is of North American origin and may have co-evolved with ancestral cervids, Odocoileus spp. . Currently, F. magna is enzootic in five major areas across the United States and Canada: (1) the northern Pacific coast, (2) the Rocky Mountain trench, (3) northern Quebec and Labrador, (4) the Great Lakes region, and (5) Gulf coast, lower Mississippi, and southern Atlantic seaboard.
Concerning the host spectrum of giant liver fluke in North America, a variety of free-living (e.g., wapiti Cervus elaphus canadensis, white-tailed deer Odocoileus virginianus, caribou Rangifer tarandus, black-tailed deer Odocoileus hemionus columbianus and mule deer Odocoileus hemionus hemionus) and domestic ruminants (cattle Bos taurus, goat Capra hircus, and sheep Ovis aries) are known to be susceptible to F. magna infection but only white-tailed deer, caribou, and wapiti contribute significantly to maintaining its population . It is assumed that the parasite was widespread in white-tailed deer in major wetland habitats throughout North America; however, the interrelationships of this host with F. magna were finely tuned due to low number of flukes within individual deer. On the other hand, potential for translocation of liver flukes in wapiti is higher due to increased F. magna egg production and subsequent release to the environment .
In Europe in the 19th and 20th centuries, game keeping and hunting had a tendency to increase the biodiversity of game species in local hunting grounds by extensive introduction and restocking of “exotic” game. As a consequence of the introduction of wapiti from North America, the giant liver fluke was unintentionally introduced to Europe, where it established three natural foci of infection; (1) northern Italy, (2) Czech Republic and south-western Poland, and (3) Danube floodplain forests (for review see ). In Europe, F. magna is considered to be an introduced invasive species with high potential to colonize new geographic territories and establish local populations. In particular, the Danube floodplain forests represent an expanding natural focus of fascioloidosis with inevitable spread of the parasite down the Danube River . In Europe, giant liver fluke has shown high capability to adapt to different intermediate aquatic snail hosts, such as species of Galba, Lymnaea, and Radix [9–11] and definitive ruminant species, in particular red deer Cervus elaphus elaphus, roe deer Capreolus capreolus, and fallow deer Dama dama [10, 12].
The origin of European populations of F. magna and the subsequent course of colonisation and migratory routes of this alien parasite in Europe recently were unravelled . Phylogenetic analyses based on concatenated cox1 + nad1 datasets revealed two genetically separated clades of European F. magna. The Italian population represented one phylogenetic clade while the second one included populations from the Czech focus and the Danube floodplain forests. Results clearly indicated that F. magna was introduced from North America to Europe at least twice; first to Italy and independently to the Czech Republic. As comparative samples, several F. magna specimens from North America were included in the above-mentioned work. The overall analysis of North American and European giant liver fluke populations clearly confirmed a western North American origin of the Italian population; these flukes clustered with specimens from Alberta (Canada) and Oregon (USA). On the other hand, representatives of the Czech Republic and Danube floodplain forests displayed close genetic relationships with parasites from the south-eastern USA. These data showed for the first time relatively high genetic molecular diversity of North American F. magna individuals . However, the limited number of the flukes from North America did not cover all five enzootic regions and were insufficient for detailed phylogeographic analysis. Therefore, the aim of the current work was to determine population structure of F. magna from a more extensive dataset including all five enzootic North American regions using previously applied mitochondrial cox1 and nad1 molecular markers [13, 14] in order to reveal genetic interrelationships of giant liver fluke on its original continent. Results were assessed in relation to the geographical distribution of its obligate definitive cervid hosts.
A total of 248 F. magna from 37 livers originating from all five North American enzootic regions were included in the analysis (Table 1). From the 248 parasites, 140 samples were newly analysed, while the remaining ones were already included in our previous analysis  (see Table 1, superscript a). Flukes were isolated directly from parenchymatous cysts and immediately rinsed in PBS buffer and fixed in 96 % ethanol.
Flukes from the enzootic region of northern Quebec and Labrador (NQL) were collected from muskox (Ovibos moschatus) and caribou (Rangifer tarandus), while specimens from the Rocky Mountain trench (RMT) originated from wapiti. In the northern Pacific coast (NPC), F. magna samples were isolated from Roosevelt elk (Cervus elaphus roosevelti) and black-tailed deer. In Minnesota, belonging to the Great Lakes region (GLR), all flukes came from white-tailed deer. Finally, the Gulf coast, lower Mississippi, and southern Atlantic seaboard (SAS) enzootic region was represented by F. magna from wapiti and white-tailed deer; in particular flukes were obtained from US states of Mississippi (MS), Louisiana (LA), Georgia (GA), Florida (FL), and South Carolina (SC).
DNA isolation, PCR amplification and sequencing
Genomic DNA was isolated from 20 mg of adult flukes using phenol:chlorophorm:isoamyl alcohol extraction and ethanol precipitation . In order to completely remove remaining PCR inhibitors, such as divalent cations and proteins, two additional wash steps using the QIAamp® DNA Kit (QIAGEN, Hilden, Germany) were employed in the DNA purification procedure. Finally, genomic DNA was diluted in deionised water and stored at −20 °C.
For amplification and sequencing of partial mitochondrial cytochrome c oxidase subunit I (cox1; 384 bp) and nicotinamide dehydrogenase subunit I (nad1; 405 bp), the previously designed primers and PCR conditions were applied [13, 14]. The PCR products were loaded on the 1 % agarose gel and purified either using the Wizard PCR purification Kit (Promega, Madison, Wisconsin) or with exonuclease I and shrimp alkaline phosphatase . Sequencing was performed using an automatic genetic analyzer Applied Biosystems 3130xl (Applied Biosystems, Foster City, California) and BigDye Terminator v3.1 Cycle sequencing kit (Applied Biosystems). Contiguous sequences were assembled and inspected for errors using Geneious version 7.1.7 (Biomatters, Auckland, New Zealand). The translation into the amino-acid sequence was performed according to the trematode mitochondrial code [16, 17]. Concatenation of cox1 and nad1 datasets was performed in SeaView 4.2 . Distribution maps showing sampling locations were prepared using the Inkscape version 0.485.0 (free download from https://inkscape.org).
Computational, statistical and phylogenetic analyses of genetic data
Phylogenetic reconstruction of the relationships between mtDNA haplotypes was performed using concatenated data with Maximum Likelihood algorithms (ML) in PhyML 3.0  and Bayesian Inference (BI) in MrBayes 3.2.3 . We used a concatenated dataset based on the experience from a previous analysis , where individual gene networks did not provide conflicting patterns. Statistical support of the resulting topologies was obtained by 1000 bootstrap replications in PhyML. MrBayes was run in two parallel runs with four chains each and two million MCMC replications sampled every 1000 generations. Twenty percent of the run was discarded as burn-in. Convergence of the parameters obtained in the two runs was inspected in Tracer 1.6 (http://tree.bio.ed.ac.uk/software/tracer). Model of molecular evolution for the BI analysis was selected in PartitionFinder  allowing separate parameter estimation for each gene and each codon position. HKY + I model was selected for the first two codon positions, whereas GTR + G was selected for the third position. TN93 model was selected in PartitionFinder for the PhyML analysis, where the usage of separate models for codon positions is not supported. Phylogenetic trees were rooted using cox1 and nad1 sequences of the closest available relative, Fasciola hepatica (GenBank Accession No. NC002546). The topology used to position a root for the trees from the ML and BI analyses was obtained in PhyML by a ML analysis of translated amino acid sequences under the LG model. Amino acid sequence analysis was used due to a relatively deep level of genetic variation between the two genera, which created a very long branch when analysed in the DNA mode.
Genealogical information contained in the two mtDNA genes was visualized using haplotype network in the software TCS 1.21 . To explore the effect of geographical distribution on the structuring of the populations, analysis of molecular variance (AMOVA) was performed using the concatenated dataset in Arlequin 188.8.131.52 . The hierarchical grouping of samples for the analysis was designed in two alternative schemes: 1) populations were grouped into five enzootic foci  of F. magna (see Table 1); 2) population groups were reorganized into six sets to better reflect the distribution and sharing of haplotypes between populations seen in TCS network. Regions in the north-west (NPC and RMT) were pooled, whereas the SAS region in the south-east was divided into three separate groups (Table 2). Significance of the obtained results was tested with 10,000 permutations of the data. Analyses of haplotype diversity (Hd), nucleotide diversity (Pi) and neutrality tests (Tajima’s D, Fu and Li’s D, Fu and Li’s F) for population sets used in the second AMOVA analysis were calculated for concatenated data in DNASP 5.10.1 . Significance of the neutrality tests was obtained with 10,000 coalescent simulations.
The analysis of 384 bp cox1 mtDNA fragment (128 amino acids, aa) and 405 bp nad1 fragment (134 aa + stop codon) revealed 32 cox1 (CO1-Ha) and 28 nad1 (ND1-Ha) mitochondrial haplotypes (Table 3). The numbering of CO1-Ha and ND1-Ha haplotypes, as presented in Table 3, adopted the strategy of numbering applied in our recently published study on genetic interrelationship of European populations of F. magna . The newly determined haplotypes of North American F. magna respected the numbering of the mentioned work and continuously proceeded in numbering new haplotypes. Since three cox1 (CO1-Ha2, 4, and 5) and four nad1 (ND1-Ha1, 2, 5, and 7) haplotypes were detected exclusively in the European populations of the parasite  they are not presented in Table 3.
Of 43 polymorphic sites detected in cox1 (transitions (ts)/transversions (tv) ratio; 40/3), 36 substitutions were not responsible for change in amino acid sequence while seven mutations underwent the non-synonymous substitutions. As for nad1, 16 substitutions out of 23 (ts/tv ratio; 19/4) were silent whereas seven mutations changed the amino acid sequence. The concatenated dataset of both cox1 and nad1 sequences (789 bp) contained 222 sequences that resulted in 50 haplotypes (Table 4; Fig. 1). The concatenated dataset contained 63 variable characters, of which 42 characters were parsimony-informative.
Phylogenetic analysis revealed two major clades (Fig. 2), one comprising samples from western enzootic regions, Rocky Mountain trench (RMT) and northern Pacific coast (NPC) (Canadian provinces BC, AB, and US state OR) (Fig. 2, clade B) whereas the second clade comprised samples from eastern enzootic regions, specifically the Great Lakes region (GLR) (US state MN), Gulf coast, lower Mississippi, and southern Atlantic seaboard enzootic region (SAS) (US states MS, LA, SC, GA, FL) and northern Quebec and Labrador (NQL) (Canadian provinces QC and NL) (Fig. 2, clade A).
The spectrum of definitive hosts sampled in the western enzootic regions was limited to three cervid hosts: wapiti, Roosevelt elk, and black-tailed deer. The respective samples of the western clade (US state OR, Canadian provinces AB and BC) were characterized by single cervid hosts; AB/wapiti, BC/Roosevelt elk, and OR/black-tailed deer (Table 1). However, as evident from the internal structure of the phylogenetic clade B, the interrelationships of western F. magna population were not related to cervid host species.
The second dominant phylogenetic clade (clade A) of eastern enzootic regions was polyphyletic. Some of the respective sublineages within clade A mirror the geography/host pattern. The most homogenous internal lineages were no. 2 (GLR-Minnesota; white-tailed deer), nos. 4 and 6 (SAS-Louisiana, Mississippi; wapiti), and no. 7 (SAS-South Carolina; white-tailed deer) (Fig. 2). One specimen from Georgia (SAS) (Ha50, lineage 3), clustered separately. Lineage no. 9 was created by haplotypes specific to SAS (US states SC, FL, GA) and two hosts; wapiti and white-tailed deer. On the other hand, lineage nos. 1 and 8 were heterologous; haplotypes within these lineages corresponded to F. magna from geographically distant regions – SAS and NQL. Even more diverse was lineage no. 5, which included representatives of all eastern populations: SAS, NQL and GLR. Additionally, besides heterogeneity in geography, the three heterologous lineages (1, 5 and 8) were characterized by the most diverse spectrum of definitive hosts – wapiti, caribou, muskox and white-tailed deer (Fig. 2).
The mitochondrial network also revealed frequent sharing of haplotypes between western enzootic regions NPC and RMT, and relative isolation of populations from Minnesota and Florida (Fig. 3). Although haplotypes specific for SAS and NQL did not create separate clusters in the network, haplotypes of these two eastern enzootic regions were not shared among different geographic areas. The only exception was Ha24, which was shared between MN (GLR) and QC (NQL) (Table 4). Regardless, samples from different locations within SAS had very little overlap. Haplotypes were shared within but not among three groups (SC + FL, GA, MS + LA) with one exception of shared haplotype between Florida and Georgia. On the contrary, populations in the north-eastern enzootic focus (NQL) shared most of their haplotypes between Quebec and Labrador despite different host origin (caribou and muskox). Surprising overlap in haplotypes was detected for individuals from Minnesota (GLR) and Quebec (NQL). Other haplotypes did not cluster according to host origin either. For example haplotypes were shared among wapiti and black-tailed deer in Alberta and Oregon.
The AMOVA analysis reflects the pattern obtained from phylogenetic analysis and haplotype network. Populations grouped according to their enzootic membership explained very little variance in the data, whereas variance in data within the areas was much higher (Table 2). On the contrary, populations organized into six groups according to the distribution of the haplotypes explained almost 73 % of the variance at the highest level (among areas), whereas less than 3 % was left among populations within the areas. Characterization of the same six groups performed in DNASP revealed some differences. Most of the population groups had similar levels of Hd and Pi diversities and non-significant results of neutrality tests (Table 5). Only the western population (AB, BC, OR) and the SAS2 group (FL, GA) had lower values of Hd and Pi. The SAS2 group also had significantly negative values in all three tests of neutrality. Negative values of neutrality tests usually are interpreted as an indication of population expansion after a bottleneck or a selective sweep (e.g., ).
Current data provide a complete picture of the genetic structure of F. magna populations from all enzootic regions in North America (NA); the Rocky Mountain trench (RMT), the northern Pacific coast (NPC), northern Quebec and Labrador (NQL), the Great Lakes region (GLR), and the Gulf coast, lower Mississippi, and southern Atlantic seaboard (SAS). The most straightforward finding was explicit genetic separation of western (NPC and RMT) and eastern (NQL, GLR and SAS) populations of the parasite.
Parasites in general rely almost entirely on their hosts for dispersal ranging from small to large scales . The geographical distributions of most parasite species are limited by the distributions of potential host species or by environmental constraints on the parasite’s rates of development . The population genetic structure of a parasite, and consequently its ability to adapt to a given host, is strongly linked to its own life history as well as the life history of its host . From this point of view, relationships of genetically and demographically variable hosts and parasites should be studied in direct collaboration .
The genetic structure of NA giant liver fluke populations needs to be related to the historical and current distribution of obligate definitive NA cervid hosts of F. magna, in particular wapiti, white-tailed deer, and caribou. According to Pybus , current NA populations of F. magna are separated in detached pockets across North America and the parasite may have co-evolved with the ancestral Odocoileus spp. Giant liver fluke has originally been widespread in white-tailed deer in major wetland habitats throughout NA where wapiti and caribou sympatric with white-tailed deer encountered F. magna in overlapping contaminated regions.
The distinction of two mitochondrial lineages of F. magna corresponds very probably to historical distribution, extirpation, re-introduction, and current occurrence of two obligate definitive hosts of giant liver fluke; white-tailed deer in eastern and wapiti in the western part of NA. Odocoileus spp. is native to North America and its populations were well established in the south-eastern USA during the Pleistocene epoch . In the early part of the 19th century unrestricted hunting, deforestation and extensive agricultural development led to dramatic declines in the white-tailed deer populations in the south-eastern USA . Implementation of an extensive restocking programme during the late 19th and in 20th century increased the white-tailed deer populations, especially in the south-east where re-establishing deer populations was considered essential .
On the other hand, wapiti is of Eurasian origin and represents phylogenetically the most derived “old world deer”, which displayed different migratory routes to the NA continent with the postglacial reopening of the passage southward from Beringia into the mid-continent . Generally, wapiti were extirpated from the Rocky Mountains of Alberta and British Columbia by the early 1900s and reduced to a few remnant populations in isolated areas. In the first half of the 20th century, wapiti were restocked from the Yellowstone National Park and expanded in number and geographic distribution within the Banff National Park (BNP). A significant exchange of wapiti between BNP and Kootenay National Park through Vermilion Pass provided a natural dispersal mechanism for F. magna from British Columbia into Alberta (see  and references therein).
Western populations of wapiti and populations of white-tailed deer in the eastern part of NA are separated by an expansive area of dry grasslands (the Great Plains region) throughout the core of Canada and USA, with very limited ecological conditions suitable for maintaining the giant liver fluke. Consequently, western and eastern populations of F. magna might be separated on a historical timeframe and evolved distinct genetic structure. Alternatively, the apparent western and eastern F. magna population structure may be a result of the lack of suitable cervid hosts following widespread extirpation of ungulate populations in eastern and central NA following European colonization that may also reduce the opportunity for genetic admixing among local F. magna populations . The effect on population segregation has been documented by studies of mtDNA phylogeography in white-tailed deer [30, 33] as well as two other hosts of F. magna, the black-tailed deer  and caribou . However, we find the historical separation more probable. The effect of such a recent bottleneck (<200 years ago) would have to be extremely strong to produce reciprocal monophyly between the two population clades A and B. Furthermore, we would expect to see strongly reduced numbers of haplotypes in local populations of F. magna. Neither the diversity of obtained haplotypes (Fig. 3) nor the results of neutrality tests (Table 5) point to such a scenario.
Phylogenetic reconstruction of the relationships between concatenated mtDNA haplotypes using Maximum Likelihood and Bayesian Inference, as well as haplotype network reflecting genealogical information congruently revealed genetic admixture of F. magna individuals from geographically distant eastern populations of the parasite (NQL, GLR, SAS), where overlapping distribution of white-tailed deer and caribou may have played an important role. Caribou were present in North America as early as the glacial periods of the middle Pleistocene in Beringia . In the 19th century, peripheral populations of caribou within the United States were extirpated and populations that occurred from Minnesota to Maine, in New York, Wisconsin and Michigan disappeared [37, 38]. While white-tailed deer was widespread in major wetland habitats throughout NA, eastern populations of caribou ranged as far south as Alabama in eastern NA  and overlapped with white-tailed deer in the Great Lakes region [40, 41]. Consequently, caribou sympatric with white-tailed deer could have encountered F. magna in overlapping contaminated wetland habitats, thus facilitating its movement into new areas and establishing new parasite populations. The present fluke population in caribou in the northeast (NQL) is most probably a residual population that survived the caribou extirpations in more southern regions. In addition, giant liver fluke further spilled over into local muskox populations that were overlapping with infected caribou in NQL region.
In contrast to the general pattern of mixing in the eastern populations, some populations seem to retain genetic distinctiveness. For example, all but one haplotype from Florida and Georgia created a separate cluster in the haplotype network (Fig. 3). Despite the samples originated from wapiti, which was probably introduced to the area recently, white-tailed deer was the naturally occurring host in the past. In the mitochondrial analysis by Ellsworth et al.  the populations of O. virginianus from south Florida showed genetic differentiation attributable to pleistocene climatic oscillations . Thus, current genetic patterns of the parasite in the same region may be remnants of past population structure the host. In relation to that it is noteworthy that the FL + GA fluke population was the only one where significantly negative values of neutrality tests were seen (Table 5).
Pertinent to the genetic relatedness of western populations of giant liver fluke (NPC and RMT) confirmed by statistical testing (AMOVA analysis) and phylogenetic analyses (Maximum Likelihood and Bayesian Inference), the distribution of six populations (subspecies) of Cervus elaphus in North America  needs to be considered. The Roosevelt elk population (C. e. roosevelti) along the Pacific coast of British Columbia and US states Washington/Oregon (NPC enzootic region) is immediately adjacent to the Rocky Mountain elk (C. e. canadensis) population in the RMT enzootic region. The Coastal Mountains to the west (where the Roosevelt elk originated) are separated from the Rocky Mountains to the east (origin of Rocky Mountain elk) only by a broad lowland area (the Interior Plateau) , which provides an extensive network of contiguous forests, wetlands, lakes and rivers. This biotope offers shared habitats for different populations of wapiti in western enzootic regions NPC and RMT and also suitable environmental conditions for perpetuation of F. magna populations. The current shared genetic makeup of the western populations of F. magna could be explained by admixing of the fluke population in association with distribution of different wapiti subspecies in this region.
Based on the mitochondrial data, it can be concluded that there are no signs of host specificity of F. magna adults towards any definitive host species; the detected haplotypes of giant liver fluke are shared amongst several host species in adjacent populations. Similar genetic patterns of geographic isolation with shared hosts broadly displayed e.g., winter ticks Dermacentor albipictus on cervids in North America , suggests widespread patterns of historic factors directly affecting multiple host parasite relationships.
Genetically diverse NA populations of F. magna reflect historical distribution, past extirpation, subsequent re-introduction, and current occurrence of the obligate definitive cervid hosts of giant liver fluke. The present study provides missing pieces of the puzzle and completes the comprehensive dataset on mitochondrial structure and population diversity of all (NA and European) F. magna populations. Comparison of original results with published data  revealed the following. Since North America is the original continent of F. magna, a high level of molecular diversity in mitochondrial haplotypes evident between and within respective enzootic regions was anticipated. While the total number of concatenated haplotypes in NA populations ranged from three (RMT) and eight (each in NQL and NPC) to 26 (SAS) (present study) only two haplotypes were determined in Danube floodplain forests and four in each of Italy and the Czech Republic . As expected, lower genetic heterogeneity was determined in newly established natural foci after introduction of non-indigenous species and the bottleneck effect is evident (data not shown).
The only haplotype common for parasites from both continents was haplotype 1, assessed in F. magna from Italy, Alberta, and Oregon, thus confirming the western NA origin of fascioloidosis in the first European focus – Italy . After comparison of present data with data achieved for European F. magna samples  two additional haplotypes (Ha40, Ha44; Table 4) detected in specimens from South Carolina (SAS) were identical with European samples from the Czech Republic and Danube floodplain forests (data not shown). This provides the first indication of the likely origin of the Czech focus of fascioloidosis. However, detailed study of genetic interrelationships of global F. magna populations should be assessed by multilocus population genetic markers, such as microsatellites, polymorphic and codominant markers which were recently designed specifically for giant liver fluke  and can provide more detailed population structuring.
Prugnolle F, Liua H, de Meeu T, Ballouxa F. Population genetics of complex life-cycle parasites: an illustration with trematodes. Int J Parasitol. 2005;35:255–63.
Štefka J, Hypša V, Scholz T. Interplay of host specificity and biogeography in the population structure of a cosmopolitan endoparasite: microsatellite study of Ligula intestinalis (Cestoda). Mol Ecol. 2009;18:1187–206.
Poulin R, Krasnov BR, Mouillot D, Thieltges DW. The comparative ecology and biogeography of parasites. Philos Trans R Soc Lond B Biol Sci. 2011;366:2379–90.
Renaud F, Clayton D, de Meeüs T. Biodiversity and evolution in host-parasite associations. Biodivers Conserv. 1996;5:963–74.
Salzet M, Capron A, Stefano GB. Molecular crosstalk in host-parasite relationships: schistosome and leech host interactions. Parasitol Today. 2000;16:536–40.
Pybus MJ. Liver flukes. In: Samuel WM, Pybus MJ, Kocan AA, editors. Parasitic Diseases of Wild Mammals. 2nd ed. Iowa: Iowa State University Press, Ames; 2001.
Králová-Hromadová I, Bazsalovicsová E, Demiaszkiewicz A. Molecular characterization of Fascioloides magna (Trematoda: Fasciolidae) from south-western Poland based on mitochondrial markers. Acta Parasitol. 2015;60:544–7.
Marinković D, Kukolj V, Aleksić-Kovačević S, Jovanović M, Knežević M. The role of hepatic myofibroblasts in liver cirrhosis in fallow deer (Dama dama) naturally infected with giant liver fluke (Fascioloides magna). BMC Vet Res. 2013;9:45.
Erhardová-Kotrlá B. The occurrence of Fascioloides magna (Bassi, 1875) in Czechoslovakia. Prague: Academia Publishing House of the Czechoslovak Academy of Sciences; 1971.
Rajský D, Patus A, Bukovjan K. Prvý nález Fascioloides magna (Bassi, 1875) na Slovensku. Slovenský veterinárny časopis. 1994;19:29–30. In Slovak.
Faltýnková A, Horáčková E, Hirtová L, Novobilský A, Modrý D, Scholz T. Is Radix peregra a new intermediate host of Fascioloides magna (Trematoda) in Europe? Field and experimental evidence. Acta Parasitol. 2006;51:87–90.
Rajský D, Čorba J, Várady M, Špakulová M, Cabadaj R. Control of fascioloidosis (Fascioloides magna Bassi, 1875) in red deer and roe deer. Helminthologia. 2002;39:67–70.
Králová-Hromadová I, Bazsalovicsová E, Štefka J, Špakulová M, Vávrová S, Szemes T, et al. Multiple origins of European populations of the giant liver fluke Fascioloides magna (Trematoda: Fasciolidae), a liver parasite of ruminants. Int J Parasitol. 2011;41:373–83.
Králová-Hromadová I, Špakulová M, Horáčková E, Turčeková Ľ, Novobilský A, Beck R, et al. Sequence analysis of ribosomal and mitochondrial genes of the giant liver fluke Fascioloides magna (Trematoda: Fasciolidae): intraspecific variation and differentiation from Fasciola hepatica. J Parasitol. 2008;94:58–67.
Werle E, Schneider C, Renner M, Volker M, Fiehn W. Convenient singlestep, one tube purification of PCR products for direct sequencing. Nucleic Acids Res. 1994;22:4354–5.
Garey JR, Wolstenholme DR. Platyhelminth mitochondrial DNA: Evidence for early evolutionary origin of a tRNA (serAGN) that contains a dihydrouridine arm replacement loop, and of serine-specifying AGA and AGG codons. J Mol Evol. 1989;28:374–87.
Ohama T, Osawa S, Watanabe K, Jukes TH. Evolution of the mitochondrial genetic code. IV. AAA as an asparagine codon in some animal mitochondria. J Mol Evol. 1990;30:329–32.
Gouy M, Guindon S, Gascuel O. SeaView Version 4: a multiplatform graphical user interface for sequence alignment and phylogenetic tree building. Mol Biol Evol. 2010;27:221–4.
Guindon S, Dufayard JF, Lefort V, Anisimova M, Hordijk W, Gascuel O. New Algorithms and Methods to Estimate Maximum-Likelihood Phylogenies: Assessing the Performance of PhyML 3.0. Syst Biol. 2010;59:307–21.
Ronquist F, Teslenko M, Van Der Mark P, Ayres DL, Darling A, Höhna S, et al. MrBayes 3.2: Efficient Bayesian phylogenetic inference and model choice across a large model space. Syst Biol. 2012;61:539–42.
Lanfear R, Calcott B, Ho SYW, Guindon S. PartitionFinder: combined selection of partitioning schemes and substitution models for phylogenetic analyses. Mol Biol Evol. 2012;29:1695–701.
Clement M, Posada D, Crandall KA. TCS: a computer program to estimate gene genealogies. Mol Ecol. 2000;9:1657–9.
Excoffier L, Lischer HEL. Arlequin suite ver 3.5: a new series of programs to perform population genetics analyses under Linux and Windows. Mol Ecol Resour. 2010;10:564–7.
Librado P, Rozas J. DnaSP v5: a software for comprehensive analysis of DNA polymorphism data. Bioinformatics. 2009;25:1451–2.
Fu YX. Statistical tests of neutrality of mutations against population growth, hitchhiking, and background selection. Genetics. 1997;147:915–25.
Dobson A, Carper R. Global Warming and Potential changes in Host-Parasite and Disease-Vector Relationships. In: Peters RL, Lovejoy TE, editors. Global warming and biodiversity. New Haven, CT: Yale University Press; 1992. p. 201–17.
van Schaik J, Kerth G, Bruyndonckx N, Christe P. The effect of host social system on parasite population genetic structure: comparative population genetics of two ectoparasitic mites and their bat hosts. BMC Evol Biol. 2014;14:18.
Roth JA, Laerm J. A late Pleistocene vertebrate assemblage from Edisto Island. SC Brimleyana. 1980;3:1–29.
Noble RE. Mississippi deer herd: past and present. Mississippi Game and Fish. 1966;29:14–5.
Ellsworth DL, Honeycutt RL, Silvy NJ, Bickham JW, Klimstra WD. Historical biogeography and contemporary patterns of mitochondrial DNA variation in white-tailed deer from the Southeastern United States. Evolution. 1994;48:122–36.
Burns JA. Mammalian faunal dynamics in Late Pleistocene Alberta. Canada Quatern Int. 2010;217:37–42.
Pybus MJ, Butterworth EW, Woodes JG. An expanding population of giant liver fluke (Fascioloides magna) in elk (Cervus canadensis) and other ungulates in Canada. J Wildl Dis. 2015;51:431–5.
Cullingham CI, Merrill EH, Pybus MJ, Bollinger TK, Wilson GA, Coltman DW. Broad and fine-scale analysis of white-tailed deer populations: estimating the relative risk of chronic wasting disease spread. Evol Appl. 2011;4:116–31.
Latch EK, Heffelfinger JR, Fike JA, Rhodes Jr OE. Species-wide phylogeography of North American mule deer (Odocoileus hemionus): cryptic glacial refugia and postglacial recolonization. Mol Ecol. 2009;18:1730–45.
McDevitt AD, Mariani S, Hebblewhite M, Decesare NJ, Morgantini L, Seip D, et al. Survival in the Rockies of an endangered hybrid swarm from diverged caribou (Rangifer tarandus) lineages. Mol Ecol. 2009;18:665–79.
Kurtén B, Anderson E. Pleistocene mammals of North America. New York: Columbia University Press; 1980.
Cringan AT. History, food habits and range requirements of the woodland caribou of continental North America. T N Am Wildl Nat Res. 1957;22:485–501.
Heinselman ML. Fire in the Virgin Forests of the Boundary Waters Canoe Area. Minnesota Quaternary Res. 1973;3:329–82.
Churcher CS, Parmalee PW, Bell GL, Lamb JP. Caribou from the late Pleistocene of northwestern Alabama. Can J Zool. 1989;67:1210–6.
Banfield AWF. The mammals of Canada. Toronto: University of Toronto Press; 1974.
Bryant LD, Maser C. Classification and distribution. In: Thomas JW, Toweill DE, editors. Elk of North America. Harrisburg: PA: Stackpole Books; 1982. p. 1–59.
Leo S, Samuel WM, Pybus MJ, Sperling FD. Origin of Dermacentor albipictus (Acari: Ixodidae) on elk in the Yukon. Canada J Wildl Dis. 2014;50:544–51.
Minárik G, Bazsalovicsová E, Zvijáková Ľ, Štefka J, Pálková L, Králová-Hromadová I. Development and characterization of multiplex panels of polymorphic microsatellite loci in giant liver fluke Fascioloides magna (Trematoda: Fasciolidae), using next-generation sequencing approach. Mol Biochem Parasitol. 2014;195:30–3.
We wish to convey our acknowledgement and indebtedness to the following colleagues for providing F. magna samples: Frank Phillips (Government Newfoundland, Quebec, Canada); Joseph Townley (Nunatsiavut Government, Newfoundland and Labrador, Canada); Don Robertson, Lorrie Dyer, Lyle Wilcox and Reg Dyer (Vancouver Island University, Canada); Vasyl Tkach (University of North Dakota, Grand Forks, North Dakota, USA); Kevin M. Keel (College of Veterinary Medicine, University of Georgia, Athens, USA); Dale Pulis (St.-Cloud, Minnesota, USA). We also express our gratitude to the Warden Service in Banff National Park, Alberta, Canada. This work was supported by the Slovak Grand Agency VEGA (project VEGA no. 2/0133/13), EU 7th Framework programme MODBIOLIN project No. 316304, and as implementation of the project Centre of Excellence for Parasitology (ITMS26220120022) supported by the Research and Development Operational Programme funded by the European Regional Development Fund (rate 0.1).
Note Nucleotide sequence data reported in this paper are available in the GenBank™, EMBL and DDBJ databases under the accession numbers KP635011 – KP635063.
The authors declare that they have no competing interests.
EB and IKH conceived the research outline and aims of the work. EB performed experimental work (DNA extraction, PCR amplification, PCR product purification, analyses of sequence data) and wrote the draft of the manuscript. IKH designed and managed the whole study, participated in analyses of sequence data, and critically reviewed the whole manuscript. JS performed statistical and phylogenetic analyses, constructed the phylogenetic tree and network phylogram, and wrote the corresponding parts of the manuscript. GM and SB sequenced the PCR products and revised the manuscript. MP provided the parasitic material, contributed to the outline of the manuscript in relation to distribution of cervid hosts and provided English language corrections. All authors read and approved the final manuscript.