Introgressive hybridization and latitudinal admixture clines in North Atlantic eels
© Wielgoss et al.; licensee BioMed Central Ltd. 2014
Received: 19 July 2013
Accepted: 13 March 2014
Published: 28 March 2014
Hybridization, the interbreeding of diagnosably divergent species, is a major focus in evolutionary studies. Eels, both from North America and Europe migrate through the Atlantic to mate in a vast, overlapping area in the Sargasso Sea. Due to the lack of direct observation, it is unknown how these species remain reproductively isolated. The detection of inter-species hybrids in Iceland suggests on-going gene flow, but few studies to date have addressed the influence of introgression on genetic differentiation in North Atlantic eels.
Here, we show that while mitochondrial lineages remain completely distinct on both sides of the Atlantic, limited hybridization is detectable with nuclear DNA markers. The nuclear hybridization signal peaks in the northern areas and decreases towards the southern range limits on both continents according to Bayesian assignment analyses. By simulating increasing proportions of both F1 hybrids and admixed individuals from the southern to the northern-most locations, we were able to generate highly significant isolation-by-distance patterns in both cases, reminiscent of previously published data for the European eel. Finally, fitting an isolation-with-migration model to our data supports the hypothesis of recent asymmetric introgression and refutes the alternative hypothesis of ancient polymorphism.
Fluctuating degrees of introgressive hybridization between Atlantic eel species are sufficient to explain temporally varying correlations of geographic and genetic distances reported for populations of the European eel.
KeywordsAnguilla anguilla Anguilla rostrata Gene flow Isolation-by-distance Simulation Migration barriers
Hybridization, the interbreeding of diagnosably divergent species, is a major focus in evolutionary studies [1–4] as it is a key concept for understanding the demographic and evolutionary cohesiveness of natural populations [5, 6]. Hybridization has great potential to rapidly introduce variability into a recipient population, if barriers to recombination can be overcome. Cooper  put forth strong experimental evidence that sexual recombination speeds up the origin and spread of adaptations in an artificial environment compared to purely clonal strains of Escherichia coli, thus supporting the theoretical predictions from the Fisher-Muller model [8, 9]. Based on the evidence gathered from empirical data, including invasive species  and signatures of massive horizontal gene transfer in a paradigmatic long-term asexual species , it might be generalized that in order to quickly adapt, any mechanism of lateral gene transfer or recombination is highly favored in unstable or novel environments. Importantly, the same might apply for populations overlapping after a phase of allopatry, as is assumed to be the case for North Atlantic eels .
Summary of the molecular and population genetics literature for European and American eels
Sample Sizes (n)
Fine et al.
Candidate markers for eel species differentiation: transferrins
Fine et al.
Heterogeneity among North Atlantic eels (not significant**)
Sick et al.
Polymorphism in American eels only, monomorphy in European eels
Pantelouris et al.
Differentiation among European continental and Icelandic eels (significant**)
Pantelouris et al.
Differentiation among American continental and Icelandic eels (significant**)
de Ligny & Pantelouris 
First available diagnostic marker: MDH;
Differentiation among American and European continental eels (significant**);
No differentiation among eels from Azores and Europe.
Williams et al.
ADH, PHI, SDH, MDH, EST
Latitudinal clines at three allozyme markers (MDH invariable)
ADH & SDH clines establish at larval American eel stages
PHI cline establishes during freshwater residency of American eels
Koehn & Williams 
ADH, PHI, SDH
Latitudinal clines at SDH & PHI loci temporally stable
ADH cline unstable and allele frequencies vary among years
Comparini & Rodinò 
Evidence for two eel species at spawning grounds in the Sargasso Sea
Williams et al.
First indication of genetic hybrids in Iceland
Avise et al.
Evidence for an eel hybrid zone: cyto-nuclear disequilibrium in Iceland
Maes & Volckaert 
Evidence against panmixia in European eels: IBD (r =0.78; P = 0.030)
Maes et al.
No interannual differentiation in European eels: no IBT (r = 0.0050, P > 0.05)
Avise et al.
Strong evidence for two eel species in the North Atlantic
Avise et al.
Evidence for an eel hybrid zone: cyto-nuclear disequilibrium in Iceland
Lintas et al.
Extensive variability in European eels
Daemen et al.
Latitudinal haplotype diversity cline in European eels
Albert et al.
Quantification of total fraction of hybrid eels in Iceland (15.5%);
Latitudinal gradient of hybrid portions in Iceland;
Evidence for high portion of later generation hybrids (30%);
Indication of higher survival rates of hybrid eels in Iceland.
Gagnaire et al.
Evidence for selection and non-neutral introgression
Daemen et al.
Low, significant genetic differentiation in European eels (F ST = 0.040; P < 0.050)
Wirth & Bernatchez 
Evidence against panmixia in European eels: IBD (r = 0.46; P < 0.0070);
Low, significant genetic differentiation in European eels (F ST = 0.0017; P = 0.0014);
Genetic intermediacy of Icelandic eels among North Atlantic locations.
Wirth & Bernatchez 
Evidence for long-term population decline in North Atlantic eels;
Differentiation among North Atlantic eels (F ST = 0.018; P < 0.0010);
No evidence against panmixia in American eels: no IBD (r = 0.0030; P > 0.40).
Mank & Avise 
Mild genetic differentiation among North Atlantic eels (G ST = 0.055; SE = 0.0049);
Genetic intermediacy of Icelandic eels among North Atlantic locations.
Dannewitz et al.
Temporal instability of IBD pattern in European eels;
Low genetic differentiation in European eels (F ST = 0.0014; P < 0.010);
Temporal genetic variation exceeds geographic variation.
Maes et al.
Evidence for interannual differentiation in Europe: IBT (r = 0.18; P = 0.043);
Contradicts allozyme pattern: no IBT, but IBD.
Palm et al.
No genetic differentiation among same-aged silvering eels between a northern and a southern European locality (F ST = −0.00003; P =0.61);
Als et al.
No evidence against panmixia in American or European leptocephali
(F ST = 0.00019, P = 0.4755): no IBD, nor IBT;
Evidence for inter-species hybridization in the Sargasso Sea
Côté et al.
No evidence against panmixia in the American eel (F ST = 0.00009; P = 0.998);
Effective population size for American eels: N E =10 532 (CI95%: 9 312–11 752)
Atlantic eel stocks have rapidly declined by over 95% compared to the levels prior to 1980 . As a consequence the European eel has been listed as critically endangered on the IUCN red list . This dramatic decrease is attributable to a combination of factors, including habitat destruction, pollution, and over-fishing , climate induced changes in the Gulf Stream circulation , as well as the recent introduction and spread of Anguillicola crassus, an exotic swimbladder nematode that was initially introduced along with infected Japanese eels [50–55]. Disentangling the population structure of both North Atlantic eel species thus has implications on international eel stock management both from conservation and fisheries perspectives.
Advanced molecular tools permitted the collection of larger and more sensitive data sets at the population level (for a schematic, more exhaustive review refer to Table 1). However, despite huge sampling efforts these studies have come to different and even contradictory conclusions. Three independent studies reported isolation-by-distance (IBD) in A. anguilla which provided evidence against the long-held paradigm of panmixia for this species based on the measurement of very low, but nevertheless highly significant genetic structure within European eels [33, 36, 39]. Wirth and Bernatchez  suggested that IBD patterns are due to a stable temporal delay of spawning migration in eels from Northern habitats, as distances are markedly extended compared to Western and Southern European populations. Alternatively, more than one spawning area might be used by different populations and, also, different currents might carry the leptocephali larvae back from the Sargasso Sea to their parent’s original freshwater habitats in Europe. Thirdly, albeit seemingly less likely, assortative mating among regional groups might be responsible for the observed IBD patterns. Finally, Maes et al.  suggested that, given the random factors affecting spawning success in the open ocean, a sweepstake strategy  might explain genetically patchy recruitment in sampling locations across Europe , leading to a weak but significant isolation-by-time (IBT) signal.
However, most recent studies did not find statistical support for population structure and rejected both IBD and IBT altogether, thus favoring panmixia in both Atlantic eel species [42–45]. In addition, temporal variation among recruits between different arrival waves within the same years significantly exceeded both the geographic and inter-annual genetic differentiation reported above [57–59]. Several points have been raised to explain the discrepancies between early evidence against panmixia and later studies. Palm et al. argued that genetic differentiation in the European eel might be largely explained by uncontrolled temporal variation between juvenile glass eel samples. Dannewitz et al.  speculated that the use of different cohorts or life stages in the initial studies with strongest support against panmixia could have lead to artifactual isolation-by-distance patterns. Lastly, after evaluating the weak genetic differentiation among North Atlantic eels using basic summary statistics, Mank and Avise  concluded that the large overlap in their allelic frequencies is generated by extensive homoplasy associated with a mutation-driven saturation effect. This argument casted doubts on the usefulness of rapidly evolving microsatellite loci for short term evolutionary and hybridization studies in eels.
Here, we investigate the extent of population genetic differentiation, and by performing simulations, ask whether IBD patterns might be produced by increasing inter-species hybridization and admixture clines. Such a scenario would also be able to explain the emergence of an IBT signal and fluctuating genetic structure signatures over time. By combining and extending two available data sets [39, 40] genotyped at nine microsatellite markers, we investigated eels sampled from the whole North Atlantic distribution area for signals of inter-species hybridization. In a first step, we inferred the admixture proportions for each individual and tested if the most extreme IBD signal detected in European eels  can be explained by latitudinal admixture and introgression clines using a simulation approach. Second, we quantified the amount of gene flow that would be necessary to generate the observed correlations, and, finally, we tested if incomplete lineage sorting or recent introgression are most likely to explain the observed lack of monophyly in most nuclear markers observed so far.
Nine microsatellite loci were used to infer intra-specific relationships among North-Atlantic eels. Inter-specific comparison of genetic differentiation for both F ST and R ST estimates were highly congruent with one another at 0.0146 (P < 0.001) and 0.0147 (P < 0.001), respectively. Thus, microsatellite markers used appear to essentially evolve by addition or removal of simple repeats in a step-wise manner. Overall, levels of observed polymorphism were high in North Atlantic eels, ranging from H o = 0.38 at locus Aro121 to H o = 0.90 at locus Ang101, with a mean of H o = 0.78. The same was true for the average numbers of alleles when correcting for sample sizes, ranging from A R = 10.9 in Aro054, to A R = 17.9 at locus Ang114. None of the 36 pairwise tests for linkage disequilibrium among loci was significant after Bonferroni correction . However, Micro-Checker revealed the presence of null alleles for the two newly added loci Ang075 and Aro146. This observation might explain the moderately high F IS values for most sampling localities (see Additional file 1) and a systematic deviation from HWE for all populations. However, after correction for null alleles using the algorithm FreeNA, all but three out of 21 sampling locations were compatible with HWE after Bonferroni correction  (see Additional file 2) and in each case where a significant deviation from HWE was detected, only a single locus contributed to the effect.
Detection of hybrid eels in Iceland
Geographic admixture and introgression clines and their impact on IBD patterns
However, these numbers exceed by far the literature values reported on clearly evidenced first generation hybrids both at continental ranges  as well as in the vicinity of the Sargasso Sea . Thus, our simulation results are not fully consistent with observation of hybrids in nature, and it seems that our assumptions of increasing clines of F1 hybrids alone represents a clear oversimplification.
Population size and migration rate parameters
Effective population size (4N i μ)
Migration rate (M/μ)
Population migration rate (2NM)
Northern A. rostrata/Icelandic A. anguilla
Northern A. rostrata/Northern A. anguilla
Intermediate A. rostrata/Western A. anguilla
Southern A. rostrata/Southern A. anguilla
The fascinating life history of Atlantic eels that involves a rare catadromous spawning migration, their economic impact and the concomitant international trade has prompted numerous studies in the past 30 years. Therefore it is astonishing that, despite the large efforts made, so many questions still remain unanswered . Our empirical and simulated data suggest that the paradigm of panmixia that was recently wavering due to major population genetics advances might still hold [42–44, 61]. Here we propose a scenario consisting of two randomly mating populations with fluctuating and clinal introgression rates; this result alone will affect the large-scale management of both endangered species since first- and later generation hybrids apparently accumulate at Northern habitats.
Hybridization pattern among North-Atlantic eels
In our attempt to detect hybrids relying on nine microsatellite markers, we came to two major conclusions. First, 16 Icelandic eels sharing an American mtDNA are genetically intermediate (i.e. hybrids) and significantly depart from pure European or American samples (P < 0.05). Second, the admixture clines observed on both sides of the Atlantic (Figure 4A) are very unlikely to be obtained by chance alone. Combinational statistics tells us that obtaining a hierarchic rank of 1, 2, 3 from three samples (American eel) has a probability P = (1/3!); obtaining a rank of 1, 2, 3, 4 from 4 samples (European eel) has a probability P = (1/4!). Therefore by combining these occurrences, we have only one chance in 144 (P = 0.007) of obtaining the result shown in Figure 4A. Though, one cannot fully exclude the presence of homoplasy, noise did not erase genetic information as previously assumed .
Here we argue that latitudinal clines of admixed individuals between European and American eels might be sufficient to explain IBD patterns observed in the European eel using neutral microsatellite markers . The latitudinal hybrid clines are most likely due to a very recent onset of gene flow after a secondary overlap and could reflect superior hybrid fitness in the northern parts of the Atlantic. The following line of argument supports this inference.
First, we confirm that Icelandic samples display high admixture proportions when compared to the continental samples. Hybrid frequencies in Iceland are high and the proportion of inferred American haplotypes is in good agreement with reported cyto-nuclear disequilibria [12, 37, 75, 88], which is a hallmark of hybrid zones [1, 2, 74]. Second, while mitochondrial lineages in eels remain 100% distinct on both sides of the Atlantic, the hybridization signal expands further to continental stocks in the nuclear genes, with decreasing latitudinal allelic richness and admixture proportions. Thus, this diffusion most likely depicts recent on-going gene flow introducing new alleles into each continental nuclear gene pool by back-crossing in the absence of maternal lineage mixing. Third, when simulating declining proportions of both F1 hybrids or more realistic fractions of admixed individuals from North to South, a stepwise decrease of approximately 5% of foreign alleles per 1,000 km of coastline would explain the IBD signal reported in Wirth and Bernatchez , whereas absent or low gene flow will fail to do so. Fourth, an isolation-with-migration model fits the data better then the alternative hypothesis of clinal ancestral polymorphism due to incomplete lineage sorting. The latter analysis also provides strong evidence for directional gene flow from American toward European eels, which is in good agreement with the absence of isolation-by-distance in the American eel [40, 45].
Importantly, previous studies suggest that actual F1 and later generation hybrids are both rare on the continents  as well as in the vicinity of the spawning grounds in the Sargasso Sea . Thus, the portions of F1 hybrids needed to explain our finding are not fully consistent with observations in nature. However, assuming a more realistic scenario, we could show that increasing latitudinal introgression clines impacted IBD in an identical way as using actual F1 hybrid clines. This can be explained by the fact that the distribution of Q-values of admixed individuals averages a value very close to the one apparent in F1 hybrids (Q = 0.55 versus 0.50, respectively).
Given semi-permeable barriers to gene flow (such as is apparent from our results) neutral markers that show intermediate linkage disequilibria with loci involved in reproductive isolation exhibit reduced gene flow compared to other neutral but unlinked markers. Such markers with reduced effective migration can display clining gene flow, as evidenced in tropical eel species which showed both admixture and introgression clines to produce an IBD pattern . While we did not explicitly distinguish between these two different clines in our simulations, we put forth strong evidence that gene flow among Atlantic eel species influence measures of population structure such as IBD.
Reunifying IBD and IBT signals
In light of our results, we are now in the position to consolidate some of the seemingly contradictory results of recent population genetics studies. Catchment areas for the earliest larval stages (group 0) are clearly overlapping among Atlantic eels ([31, 44, 90]; but see ). In contrast, the number of hybrids sampled in the Sargasso Sea was very low, and centred in the spawning areas of the European eel . Thus, there are reasons to believe that actual spawning areas are partially “allopatric”, possibly due to strong oceanic fronts that separate spawning eels within distinct areas . Interestingly, Albert et al. showed that the proportion of hybrids reaching Icelandic waters peaked in 1999 and 2000, while it decreased from 2000 to 2003. This trend could explain the strong patterns of isolation-by-distance in two independent studies that relied on samples from 1997 to 1999 [33, 39]; while later studies did not find support for population structure at all [42–44, 57–59]. The former could then be expected to lead to the clear pattern of isolation-by-time in genotype data that spans the temporal time-frame in question . Thus, as a working hypothesis, sudden bursts of hybridization and hybrid arrivals could explain temporally unstable patterns of IBD, with a clear signal of IBT. Accordingly, the fluctuations in hybrid recruit proportions would reflect changes in the degree and timing of overlap of spawning grounds in an unstable oceanic environment . Since only a very small percentage (0.8%) of hatched eel larvae is expected to reach suitable habitats within two years of dispersal , the lower the overlap of each species spawning grounds, the lower and less significant the correlation coefficients of the IBD signal on either continent would be, such as found in , or . A recent modelling approach has highlighted that even minimal mixing among any existing sub-populations in the spawning grounds would entirely erase detectable signatures of genetic divergence in arriving eel recruits at the continents . Thus, among the alternative explanations they offer, we assume that admixture and introgression clines can completely explain our findings.
Evidence for cyto-nuclear incompatabilites and asymmetric gene flow
Even though Iceland is at the intersect of the North Atlantic eel distributions, Icelandic eels show a clear affinity towards the European eel gene pool: the suspect individuals carrying American haplotypes in Iceland are hybrids, and no pure American migrants were recovered there. This observation is in accordance with the results collected by Albert et al. who could not identify pure American eels in Iceland based on a much larger sample size and denser geographic coverage. Recently, Frankowski et al. reported the regional sampling of eels with American mitochondria in Europe. However, these eels were sampled from both aquacultures and frequently stocked German rivers. In the absence of nuclear genetic evidence, these eels were most likely introduced American eels, which escaped into natural habitats, reminiscent of the situation in Japan, where imported foreign eels are frequently detected in rivers . If anything, we would assume that stocking activities would rather erase our fine signal of clinal gene flow.
While American and European eels are clearly separated according to mitochondrial data, interbreeding signatures gradually increased toward Northern habitats in our data. This cyto-nuclear disequilibrium had already been noticed earlier by Avise et al.. Importantly, using an RNA-seq approach, Gagnaire and colleagues  delivered a very intriguing molecular mechanism that could explain this finding, and has implications for asymmetric inter-species hybridization. Two components of the ATP synthase complex, one encoded in the mitochondrium (mt-atp6), the other in the nucleus (atp5c1), show significant signs of diversifying selection between the two Atlantic eel species possibly due to co-adaptation in at least one of the species. Thus, conflicts between certain cyto-nuclear combinations could render energy metabolism dysfunctional, diminishing the number of possible hybrid combinations as a result. Interestingly, the authors deliver strong evidence for asymmetric gene flow in several nuclear genes, most of which are directed into the European gene pool. This observation is matching our data based on neutral microsatellite markers.
Dispersal time and development
The most obvious difference between the two North Atlantic eel species is their divergent larval dispersal strategy, coupled with differences in the onset of metamorphosis [46, 96, 97]. Data for the timing of metamorphosis and recruitment are based on otolith daily increments . However, these calculations provided contrasting results between research groups, and, even more worrying, “back-calculated” larval migration times did not match the field observations . These inconsistencies suggest that daily increments are not suitable for calculating the timing of dispersal, but rather that they provide a descent proxy for the timing of metamorphosis of leptocephali into glass eels, which takes place on the verge of open ocean and continental shelves . Kettle and Haines  predicted a minimum of two years for larvae to successfully cross the Atlantic solely by passive drift, which was concordant with the early estimates based on larval growth rates by Schmidt . This would suggest a three- to six-fold longer migration time for the European eel compared to its American congener (A. rostrata), which arrives within some 7–12 months ([76, 96, 98, 100–102]; but see ). That fact might explain why hybrid eels, are also most abundant at intermediate ranges (with on average intermediate arrival times at the continental shelf; ).
Hybrid fitness and ecological peculiarities in Iceland
In light of the frequent detection of hybrids in Iceland, we speculate on the possibility that Nordic habitats, especially Iceland and Greenland  might depict ecotone habitats compared to the much warmer continents. Freshwater habitats in these areas were definitely uninhabitable during extensive glaciation events in the Pleistocene and must have been colonized afterward (i.e., not earlier than 10,000 years ago). It is suggested that environmental peculiarities characterize eel habitats in Iceland, as freshwater temperatures are typically much lower there as compared to most potential continental habitats [37, 76]. This anomaly is directly reflected in Icelandic glass eel’s otolith microstructure and microchemistry, which lack both the usual sharp decrease in Sr:Ca ratios and elevated increment accumulations . In addition, it is worth noting that the diffuse otolith increment zone after metamorphosis has never been observed in any other eel species outside Iceland. Thus, environmental opportunity might favour F1 hybrids in Iceland and other Nordic habitats.
In support of this argument, Albert et al. quantified the hybrid proportions in both recruiting and resident eel stages over several years. A total of 70% of putative hybrids fell into the first generation category, whereas 30% belonged to later generation hybrids. They observed an approximately two-fold increase in hybrid proportions from the recruiting glass eel to the resident yellow eel stages. These results suggest a higher hybrid survival upon residency. The presence of second and later generation hybrids indeed demonstrates that hybrids transmitted their genes to the next generations and this would explain why the admixture extends further South on both continents.
The calculations presented here are based on the strongest IBD detected and therefore represent the highest gene flow estimates that could be extrapolated from the genetic data. The age-dependent fitness differences  were not taken into account here. Therefore, the hybridization rate is also somewhat biased, given that the northern A. anguilla samples in the analysed data set are collected from yellow eels fin clips and not from glass eels. However, as has been mentioned elsewhere , temporal inter-cohort variation is not expected to lead to a clinally structured partition of genetic variation, but rather to unordered differentiation patterns. The rather weak traces of admixture detected by Albert et al. in the continental samples only moderately fit with our results. However, it is worth noting that the few American eels with a later generation hybrid signature belonged to the most northern samples (Medomak River and Boston Harbour) and as such follow the cline presented here. Moreover, Als et al. only detected one hybrid leptocephalus out of 388 eel larvae in the Sargasso Sea. This very low hybrid prevalence is strikingly different from what is observed at later developmental stages in Iceland . Therefore, selection gradients and landscape genetics must be the main evolutionary forces shaping the latitudinal cline.
Leaving aside the statistical issues (running MCMCs with or without prior), the nature of the genetic markers (co-dominant versus dominant), their relative merits and the quantitative aspects, we now have a couple of clear scenarios and hypotheses to test. The fast developing next-generation sequencing field has already provided us with primer molecular data, that open up the opportunity to enlighten the ominous marine part of the eel’s life cycle (SNPs ). Our understanding of the evolutionary dynamics of eel stocks and the selective factors that shape introgressive hybridization in North-Atlantic eels is thus currently going through a quantum leap. Especially if the Northern distribution range can be more systematically and extensively sampled from now on, we will be able to appreciate the actual rate of hybridization among the two Atlantic eel species and follow its fluctuation in space and time. This would also have important practical implications in the monitoring of conservation strategies that have been enacted in Regulation (EC) No 1100/2007 by the European Union  in response to the dramatic decline in eel recruitment in the past decades . As a prerequisite for appropriate actions, the fact that admixed individuals and hybrids preferentially tend to accumulate in Northern habitats must be accounted for. For example, bursts of increased inter-species hybridization in certain years might lead to increased numbers of northward-bound recruits, and reduced catchment in more Southern ranges.
Here, we explicitly focused on the influence of hybridization on genetic differentiation signatures in North Atlantic eels. We relied on the available nuclear microsatellite genotypes and mitochondrial sequence data of 1,263 samples from both continental American and European eels as well as from intermediate Iceland. When simulating continuously increasing proportions of F1 hybrid individuals from the southern to the northern-most locations in Europe, highly significant isolation-by-distance patterns arose, that are reminiscent of previously published data . Therefore, introgressive hybridization alone is sufficient to explain the correlation of geographic and genetic distances reported for the European freshwater eel. Moreover, contrasting signals among nuclear and mitochondrial lineages suggest a recent onset of gene flow, most likely after glacial retreat following the last Ice Age (vicariant scenario ). Importantly, our findings are in agreement with previous results on genetic isolation patterns in European eels, either based on geography (IBD [33, 39]), or inter-annual genetic composition (IBT ). If we can assume that the known overlap of the two species’ spawning grounds  is annually changing, the IBD signal should decline in some years (with low overlap), and increase in others. This will then automatically lead to an even higher IBT signal among annual recruitment waves within species. The clear separation of mitochondrial lineages might be explained by cyto-nuclear disequilibria recently detected . Thus, the lack of pure American expatriates in Iceland (see also ), could be due to the American eel’s much faster ontogenetic development and metamorphosis which might prevent its settlement in this northern region (“early ripe, early rotten”). All in all, evidence for hybrid survival in Iceland  favours the introgression hypothesis followed by subsequent backcrossing. This might not only hold for Iceland but for other Nordic regions as well, and might generate a North-to-South-hybrid gradient in both Atlantic eel species.
Samples and data sets combination
A total of 1263 North Atlantic eels were collected in 1999. The study includes 12 European (n = 561), one Icelandic (n = 300) and eight North American (n = 402) samples (see [39, 40]); see Additional file 1 for an overview of all eels samples.
Microsatellite genotyping and summary statistics
Nearly all North Atlantic eel specimens (n = 1042) have been successfully genotyped using nine microsatellite markers. Original genotypes for seven microsatellite loci [39, 40] were supplemented with two additional loci, Ang075 (Genbank AF237903; Primer sequences, Ang075-F, TATCAGGAACTCGATACGCC, and, Ang075-R, ACGCATCACCAGCCCTTGC), and Aro146 (AF237904; Aro146-F, CAGTTATCCATCTACAGGTG, and, Aro146-R, GAAATAAGAGAATGAGACTCTG). The same genotyping procedure was applied to the other eel species. The fragment sizes were determined by reference to a size standard using the software Genescan version 2.1 and Genotyper version 2.0 (Applied Biosystems Inc., Foster City, CA). Allelic diversity, genetic variation and deviation from Hardy-Weinberg Equilibrium (HWE) were calculated with Genepop on the web  and Genetix version 4.05 . All microsatellites were tested for null alleles using Micro-Checker. Allelic diversity and private allelic richness were also inferred after correcting for unequal sampling sizes using HP-rare. Pairwise genetic differentiation was calculated with Arlequin version 3.1  and statistical significance was inferred after 10,000 permutations.
Mitochondrial DNA sequencing and species identification
All North Atlantic eel samples were screened by PCR-RFLP analysis of a 362 bp segment of the Cytochrome b gene (cytb) . This test is based on a diagnostic HinfI restriction site, specific to American eels. In order to confirm this quick screening approach, the cytb amplicon was directly sequenced in all Iceland individuals with American haplotypes using an ABI 377 automated sequencer (GenBank accession numbers: KJ546041, KJ546042, KJ546043, KJ546044, KJ546045, KJ546046, KJ546047, KJ546048, KJ546049, KJ546050, KJ546051, KJ5x46052, KJ546053, KJ546054, KJ546055, KJ546056). Incorporating known sequence data from Genbank, a haplotype network based on maximum parsimony was constructed in Tcs version 1.20  relying on an alignment of 278 bp (GenBank accession numbers: AB021767, AB021776, AF006714, AF006715, AF006716, AF006717, AF165069, AF172394, AF368238, AF368239, AF368240, AF368241, AF368242, AF368243, AF368244, AF368245, AF368246, AF368247, AF368248, AF368249, AF368250, AF368251, AF368252, AF368253, AF368254, AF485271, AF485272, AF485273, AF485274, AF485275, AF485276, D28775, D84302, EF427617, EF427618, EU223996, EU223997, EU315235, EU315236, EU315237, EU315238, EU315239, EU315240, EU315241, EU315242, EU315243, EU492326, EU492327, M85080). This network was compared to a maximum likelihood tree to correct for ambiguous, multiple connections. In brief, all redundant sequences were removed from the data set, and the best-fit model of sequence evolution (HKY + I + G) was chosen based on the agreement of all Information Criteria (cAIC, AIC2, BIC) used in Modelgenerator version 0.85 . Based on the estimated shape parameter of the γ-distribution, α = 0.02, the proportion of invariable sites, Pinvar = 0.8849, and an expected transition-transversion ratio of 6.28, a phylogenetic maximum-likelihood tree was inferred in PhyML Online version 3.0  and compared to the haplotype network derived by the parsimony approach in Tcs.
Assessing individual ancestry proportions
Individual ancestry proportions in North Atlantic eel species were estimated using Structure version 2.3.2[80–83] performing 100,000 burn-in steps followed by 1,000,000 MCMC repeats and three iterations to check for Markov chain convergence. The most likely number of populations was assessed using Evanno’s ad-hoc statistic ∆K. First, we ran an admixture model without priors to infer the individual ancestry proportions and to detect putative clines of admixture over the whole sampling area. Second, an admixture model including prior information on sampling localities was used (except for Icelandic eels which were included without population information), to estimate the admixture proportions more precisely. A threshold level for the posterior probability of Q = 0.9 was used, in order to reach maximum assignment efficiency. Due to the low degree of genetic differentiation between the two North Atlantic eel species (Table 1), the efficiency and performance of hybrid identification may have been hampered by the limited number of markers available . Thus, to test if those Icelandic eels with American mitochondrial haplotypes in our sample (n = 16) are actually intermediate rather than pure expatriates, another test, analogous to an urn model, was performed. Groups of 16 individuals were drawn 1000 times at random from either Anguilla gene pools, excluding Iceland. Frequencies of average ancestry proportions were plotted, and the average value for the 16 suspected hybrids was compared to each gene pool, and to the average of the remaining Icelandic samples.
Inference of geographic admixture clines
Support for up to three geographical groupings was apparent from distance-based phenograms [39–41], the European eel samples were clustered roughly according to the 7-year-average of sea surface temperatures (SST) in the North Atlantic Ocean , into Northern (Baltic Sea, Elbe, Imsa), Western (Grand-Lieu, Couesnon, Severn) and Southern (Mediterranean Sea, Minho, Adour) groups. Likewise, the American eels were split into three latitudinal groups, Southern (St. Johns River, South Edisto); intermediate (Wye River, Hudson River, Boston Harbour,); and Northern (Prince Edwards Island, Trinité, Medomak River) samples. How biologically meaningful this categorization is was investigated statistically in R 2.10.1, on the basis of genotypes, by extracting each individual’s ancestry proportions from the most highly supported scenario in Structure and assigning them to their respective geographic group. A Boxcox transformation of the data was applied to find the transformation term that maximizes normal distribution of the data, and thus to increase the statistical power of the test. Normality was tested using the Shapiro-Wilk test, and significant differentiation between pairs of continental groups within species was checked for by applying pairwise non-parametric Wilcoxon tests.
Simulation of hybridization clines and gene flow quantification
We applied a simulation approach to test our hypothesis, that the IBD signal detected in Wirth and Bernatchez  for the European eel could have been generated by the steady increase either inter-species F1 hybrids or admixture clines from Southern to Northern Europe. Instead of the real dataset, which comprised eels from 12 sampling locations along the coastlines of Europe, we used a simple mating scheme to generate 12 virtual hybrid populations differing by a constantly increased level of inter-species hybridization. In essence, for assessing the impact of F1 hybrids on IBD signals, we generated the data in five steps: A) We identified the purest sampling localities of each eel species by considering average admixture coefficients as derived from Structure, thereby reducing any confounding effect due to actual admixture. B) We generated 12 synthetic hybrid populations. For the first hybrid population, we did not allow for admixture, and thus generated it by randomly mating 50 “parental” pairs drawn at random from the purest European eel sample (sampling with replacement). C) For the second and each subsequent population, we increased the proportion of American eel parents sampled at random by 1%, that is, we added one additional American eel parent per population for the 50 crossings. D) We followed the same approach for a total of four different gene flow scenarios: linear increase of interspecies hybridization from one population to the next by 1%, 2% 3% and 4%. For each of the four datasets, pairwise genetic D CE distances  were calculated among populations using Phylip version 3.68 . E) Given both significant IBD and that populations follow a cline of admixture from South to North in the real data, we assigned the hypothetical pairwise genetic distance data to the real geographic distances among the genotyped populations. Finally, in order to estimate the gene flow necessary to establish the degree of correlation (and significance) observed in the original dataset , both correlation coefficients (r) and the slope of the trend line (b) were matched against the simulated gene flow levels.
In a second more realistic simulation approach, we created different degrees of admixed genotypes by sampling alleles randomly from the two species gene pools and this in direct proportion given by the actual admixture values apparent from the STRUCTURE output. Following this approach, we obtained virtual admixed genotypes nearly satisfying an average admixture proportion of 50% (expected for F1 hybrids), at Q = 0.55. We performed this simulation using an R script relying on packages adegenet [119, 120], ade4 , ape [122, 123], pegas  and plyr . The simulation consisted of these principle steps: i) We first sorted the genotype dataset according to Q-values inferred using the aforementioned best run in STRUCTURE and set cut-offs for “pure” American (Q < 0.1) and “pure” European (Q > 0.9) eels. From this we calculated the respective allele frequency spectrum for both species. ii) We used the remaining individuals to estimate the distribution of all admixture proportions, i.e., we assessed the probability density of Q, as discretized in bins of equal sizes of 0.05 ranging [0.10; 0.90]. In order to create admixed individuals, we randomly sampled the Q-proportions using multinomial sampling from the observed proportions of Q-values in each bin. More specifically, for any given admixed genotype with a certain Q-value the appropriate proportion of alleles inherited from either A. rostrata, or A. anguilla were determined for all nine loci using binomial sampling, and then, using multinomial sampling, alleles were randomly drawn from the allele frequency spectrum of each locus in each “pure species” pool. iii) Finally, we assigned linearly increasing fractions of admixed genotypes of 1%, 2%, 3% or 4%, analogous to the simplified simulation using F1 hybrids, and this again for 100 virtual individuals per population with European eels being the focal species receiving introgressed American alleles.
Lastly, we used the Isolation-with-Migration model implemented in the software IMa2  to test the hypothesis that the pattern observed in Wirth and Bernatchez  results from varying gene flow rates between European and American populations as opposed to ancestral polymorphism due to incomplete lineage sorting. By means of Markov Chain Monte Carlo simulations, migration rates (m = M/μ with M being the mutation rate per generation per gene copy and μ the mutation rate), effective population sizes (q = 4N i μ with N i , the effective size of population i) and divergence time (t = Tμ where T is the time since the common ancestry) were estimated, which allowed us to distinguish gene flow from the other two estimated demographic parameters that are dependent on ancestral polymorphism. IMa2 can deal with more than two populations, but it requires a very large number of loci . Thus, we favoured a two-population model and analysed four pairs of population to estimate migration rates between: a) northern A. rostrata and Icelandic populations; b) northern populations; c) intermediate populations; and d) southern populations of A. rostrata and A. anguilla. Data sets included all genetic markers, namely, the nine microsatellite loci, and the cytb fragment reduced to the part containing the recognition site, because this 4-base pair stretch was available for the total dataset. Mutation models used for analyses were a stepwise mutation model for microsatellite loci and the HKY model for the cytb fragment. We used an exponential prior on migration to explore both low and high migration rates, even with the modest number of microsatellites available. Moreover, this model makes tests of migration more conservative than a uniform prior [85, 126]. Preliminary runs were performed to set the upper bounds (or the mean for migration rate) of the prior distribution of each parameter. The upper limit for the prior for the population size and divergence time parameters were set to 500 and 100, respectively (except for the population pair northern A. rostrata-Icelandic population for which the prior for t was 30) and the mean of the exponential distribution for the migration rate parameter was set to 2. Since none of the tested values gave satisfying swapping rates between successive chains, we did not run multiple Metropolis-coupled chains and instead, conducted runs with burn-in periods of 30 million followed by a chain of more than 200 million steps, and this for three independent runs with different random number seeds. To assure consistency of results across runs and test for convergence of chains, we compared autocorrelations of values over time, as well as parameter values generated for both the first and second halves of each run .
Availability of supporting data
The data sets supporting the results of this article are available in TreeBASE, http://purl.org/phylo/treebase/phylows/study/TB2:S15621, and LabArchives, http://dx.doi.org/10.6070/H4FX77C8.
Cavalli-Sforza & Edwards’s genetic distance
Hasegawa Kishino and Yano’s molecular genetic substitution model
Isolation-with-migration program version 2
Number of populations
Evanno’s ad-hoc statistic
Markov chain monte Carlo
Slatkin’s fixation index
Sea surface temperature.
We owe thanks to the editor and an anonymous referee for their constructive comments, as well as our collaborators who contributed specimens or otherwise helped with sampling. We thank Annabel Whibley, Pascale Chesselet and Louis Bernatchez for stimulating comments on the manuscript. This work was supported by grants (WI-2710) from the Deutsche Forschungsgemeinschaft (DFG) to T.W. and A.M. as well as grants from the Landesgraduiertenförderung Baden-Württemberg to S.W., the University Konstanz to S.W. and A.M, and an EU Marie Curie PEOPLE postdoctoral fellowship for Career Development to S.W. (grant number FP7-PEOPLE-2012-IEF-331824). The project was also supported by the network “Bibliothèque du vivant” funded by the CNRS, the Muséum National d’Histoire Naturelle and the CEA (Centre National de Séquençage).
- Barton NH, Hewitt GM: Analysis of hybrid zones. Annu Rev Ecol Syst. 1985, 16: 113-148. 10.1146/annurev.es.16.110185.000553.Google Scholar
- Barton NH, Hewitt GM: Adaptation, speciation and hybrid zones. Nature. 1989, 341 (6242): 497-503. 10.1038/341497a0.PubMedGoogle Scholar
- Harrison RG: Pattern and process in a narrow hybrid zone. Heredity. 1986, 56: 337-349. 10.1038/hdy.1986.55.Google Scholar
- Harrison RG: Hybrid Zones: windows on Evolutionary Processes. Oxford surveys in Evolutionary Biology. Edited by: Futuyma D, Antonovics J. 1990, 69-128. vol. 7Google Scholar
- Waples RS, Gaggiotti O: What is a population? An empirical evaluation of some genetic methods for identifying the number of gene pools and their degree of connectivity. Mol Ecol. 2006, 15 (6): 1419-1439. 10.1111/j.1365-294X.2006.02890.x.PubMedGoogle Scholar
- Wirth T, Le Guellec R, Vancassel M, Veuille M: Molecular and reproductive characterization of sibling species in the European earwig (Forficula auricularia). Evolution. 1998, 52 (1): 260-265. 10.2307/2410942.Google Scholar
- Cooper TF: Recombination speeds adaptation by reducing competition between beneficial mutations in populations of Escherichia coli. PloS Biol. 2007, 5 (9): 1899-1905.Google Scholar
- Fisher RA: The Genetical Theory of Natural Selection. 1930, Oxford, U.K.: Oxford Univ. PressGoogle Scholar
- Muller HJ: Some genetic aspects of sex. Am Nat. 1932, 66: 118-138. 10.1086/280418.Google Scholar
- Cox GW: Alien Species and Evolution. 2003, Washington, DC, U.S.A.: Island PressGoogle Scholar
- Gladyshev EA, Meselson M, Arkhipova IR: Massive horizontal gene transfer in bdelloid rotifers. Science. 2008, 320 (5880): 1210-1213. 10.1126/science.1156407.PubMedGoogle Scholar
- Avise JC, Nelson WS, Arnold J, Koehn RK, Williams GC, Thorsteinsson V: The evolutionary genetic status of Icelandic eels. Evolution. 1990, 44 (5): 1254-1262. 10.2307/2409286.Google Scholar
- Tesch FW: The eel. 2003, Oxford, U.K.: Blackwell PublishingGoogle Scholar
- Knights B: A review of the possible impacts of long-term oceanic and climate changes and fishing mortality on recruitment of anguillid eels of the Northern Hemisphere. Sci Total Environ. 2003, 310 (1–3): 237-244.PubMedGoogle Scholar
- Kleckner RC, McCleave JD: The Northern limit of spawning by Atlantic eels (Anguilla spp.) in the Sargasso Sea in relation to the thermal fronts and surface water masses. J Mar Res. 1988, 46 (3): 647-667. 10.1357/002224088785113469.Google Scholar
- Schmidt J: The breeding places of the eel. Annu Rep Smithsonian Institute. 1924, 1925: 279-316.Google Scholar
- Schoth M, Tesch FW: Spatial distribution of small 0-group Anguilla larvae in the Sargasso Sea. Helgolander Meeresuntersuchungen. 1982, 35 (3): 309-320. 10.1007/BF02006139.Google Scholar
- Tesch FW, Kracht R, Schoth M, Smith DG, Wegner G: Report on the eel expedition of FRV (Anton Dohrn) and R. K. “Friedrich Heincke” to the Sargasso Sea 1979. ICES CM. 1979, 1979/ M6: 1-29.Google Scholar
- Tucker DW: A new solution to the Atlantic eel problem. Nature. 1959, 183 (4660): 495-501. 10.1038/183495a0.Google Scholar
- de Ligny WD, Pantelouris EM: Origin of the European eel. Nature. 1973, 246 (5434): 518-519. 10.1038/246518a0.Google Scholar
- Avise JC, Helfman GS, Saunders NC, Hales LS: Mitochondrial DNA differentiation in North Atlantic eels: population genetic consequences of an unusual life history pattern. Proc Natl Acad Sci U S A. 1986, 83 (12): 4350-4354. 10.1073/pnas.83.12.4350.PubMedPubMed CentralGoogle Scholar
- Minegishi Y, Aoyama J, Inoue JG, Miya M, Nishida M, Tsukamoto K: Molecular phylogeny and evolution of the freshwater eels genus Anguilla based on the whole mitochondrial genome sequences. Mol Phylogenet Evol. 2005, 34 (1): 134-146. 10.1016/j.ympev.2004.09.003.PubMedGoogle Scholar
- Aarestrup K, Okland F, Hansen MM, Righton D, Gargan P, Castonguay M, Bernatchez L, Howey P, Sparholt H, Pedersen MI, McKinley RS: Oceanic Spawning Migration of the European Eel (Anguilla anguilla). Science. 2009, 325 (5948): 1660-1660. 10.1126/science.1178120.PubMedGoogle Scholar
- Fine JM, Drilhon A, Amouch P, Boffa G: Existence of serum groups in Anguilla anguilla L. Demonstration by electrophoresis and autoradiography of several transferrin types. C R Acad Sci Hebd Seances Acad Sci D. 1964, 258: 753-756.PubMedGoogle Scholar
- Fine JM, Drilhon A, Ridgway CJ, Amouch P, Boffa G: Groups of transferrins in the Anguilla species. Differences in the phenotypic frequencies of transferrins in Anguilla anguilla and Anguilla rostrata. C R Acad Sci Hebd Seances Acad Sci D. 1967, 265: 58-60.PubMedGoogle Scholar
- Sick K, Bahn E, Frydenberg O, Nielsen JT, Von Wettstein D: Haemoglobin polymorphism of the American fresh water eel Anguilla. Nature. 1967, 214: 1141-1142. 10.1038/2141141a0.PubMedGoogle Scholar
- Pantelouris EM, Árnason A, Tesch F-W: Genetic Variation in the Eel. II. Transferrins, haemoglobins and esterases in the eastern North Atlantic. Possible interpretations of phenotypic frequency differences. Genet Res. 1970, 16: 277-284. 10.1017/S0016672300002536.PubMedGoogle Scholar
- Pantelouris EM, Árnason A, Tesch F-W: Genetic variation in the eel. III. Comparisons of Rhode Island and icelandic populations. Implications for the Atlantic eel problem. Mar Biol. 1971, 9: 242-249. 10.1007/BF00351386.Google Scholar
- Williams GC, Koehn RK, Mitton JB: Genetic differentiation without isolation in the American eel, Anguilla rostrata. Evolution. 1973, 27 (2): 192-204. 10.2307/2406960.Google Scholar
- Koehn RK, Williams GC: Genetic differentiation without isolation in the american eel, anguilla rostrata. II. Temporal stability of geographic patterns. Evolution. 1978, 32 (3): 624-637. 10.2307/2407728.Google Scholar
- Comparini A, Rodino E: Electrophoretic evidence for two species of Anguilla leptocephali in the Sargasso Sea. Nature. 1980, 287 (5781): 435-437. 10.1038/287435a0.Google Scholar
- Williams GC, Koehn RK: Icelandic eels: evidence for a single species of Anguilla in the North Atlantic. Copeia. 1984, 1: 221-223.Google Scholar
- Maes GE, Volckaert FAM: Clinal genetic variation and isolation by distance in the European eel Anguilla anguilla (L.). Biol J Linn Soc. 2002, 77 (4): 509-521. 10.1046/j.1095-8312.2002.00124.x.Google Scholar
- Maes GE, Pujolar JM, Hellemans B, Volckaert FAM: Evidence for isolation by time in the European eel (Anguilla anguilla L.). Mol Ecol. 2006, 15 (8): 2095-2107. 10.1111/j.1365-294X.2006.02925.x.PubMedGoogle Scholar
- Lintas C, Hirano J, Archer S: Genetic variation of the European eel (Anguilla anguilla). Mol Mar Biol Biotechnol. 1998, 7: 263-269.PubMedGoogle Scholar
- Daemen E, Cross T, Ollevier F, Volckaert F: Analysis of the genetic structure of European eel (Anguilla anguilla) using microsatellite DNA and mtDNA markers. Mar Biol. 2001, 139: 755-764. 10.1007/s002270100616.Google Scholar
- Albert V, Jonsson B, Bernatchez L: Natural hybrids in Atlantic eels (Anguilla anguilla, A. rostrata): evidence for successful reproduction and fluctuating abundance in space and time. Mol Ecol. 2006, 15 (7): 1903-1916. 10.1111/j.1365-294X.2006.02917.x.PubMedGoogle Scholar
- Gagnaire PA, Albert V, JÓNsson B, Bernatchez L: Natural selection influences AFLP intraspecific genetic variability and introgression patterns in Atlantic eels. Mol Ecol. 2009, 18 (8): 1678-1691. 10.1111/j.1365-294X.2009.04142.x.PubMedGoogle Scholar
- Wirth T, Bernatchez L: Genetic evidence against panmixia in the European eel. Nature. 2001, 409 (6823): 1037-1040. 10.1038/35059079.PubMedGoogle Scholar
- Wirth T, Bernatchez L: Decline of North Atlantic eels: a fatal synergy?. Proc R Soc B Biol Sci. 2003, 270 (1516): 681-688. 10.1098/rspb.2002.2301.Google Scholar
- Mank JE, Avise JC: Microsatellite variation and differentiation in North Atlantic eels. J Hered. 2003, 94 (4): 310-314. 10.1093/jhered/esg062.PubMedGoogle Scholar
- Dannewitz J, Maes GE, Johansson L, Wickström H, Volckaert FAM, Järvi T: Panmixia in the European eel: a matter of time. Proc R Soc B Biol Sci. 2005, 272 (1568): 1129-1137. 10.1098/rspb.2005.3064.Google Scholar
- Palm S, Dannewitz J, Prestegaard T, Wickstrom H: Panmixia in European eel revisited: no genetic difference between maturing adults from southern and northern Europe. Heredity. 2009, 103 (1): 82-89. 10.1038/hdy.2009.51.PubMedGoogle Scholar
- Als TD, Hansen MM, Maes GE, Castonguay M, Riemann L, Aarestrup KIM, Munk P, Sparholt H, Hanel R, Bernatchez L: All roads lead to home: panmixia of European eel in the Sargasso Sea. Mol Ecol. 2011, 20 (7): 1333-1346. 10.1111/j.1365-294X.2011.05011.x.PubMedGoogle Scholar
- Côté CL, Gagnaire PA, Bourret V, Verreault G, Castonguay M, Bernatchez L: Population genetics of the American eel (Anguilla rostrata): FST = 0 and North Atlantic Oscillation effects on demographic fluctuations of a panmictic species. Mol Ecol. 2013, 22 (7): 1763-1776. 10.1111/mec.12142.PubMedGoogle Scholar
- van Ginneken VJT, Maes GE: The european eel (Anguilla anguilla, Linnaeus), its lifecycle, evolution and reproduction: a literature review. Rev Fish Biol Fisher. 2005, 15 (4): 367-398. 10.1007/s11160-006-0005-8.Google Scholar
- Freyhof J, Kottelat M: Anguilla anguilla. IUCN 2013 IUCN Red List of Threatened Species. 2013, http://www.iucnredlist.org,Google Scholar
- Feunteun E: Management and restoration of european eel population (anguilla anguilla): an impossible bargain. Ecol Eng. 2002, 18 (5): 575-591. 10.1016/S0925-8574(02)00021-6.Google Scholar
- Durif CM, Gjosaeter J, Vollestad LA: Influence of oceanic factors on Anguilla anguilla (L.) over the twentieth century in coastal habitats of the Skagerrak, southern Norway. Proc Biol Sci. 2011, 278 (1704): 464-473. 10.1098/rspb.2010.1547.PubMedPubMed CentralGoogle Scholar
- Lefebvre F, Wielgoss S, Nagasawa K, Moravec F: On the origin of Anguillicoloides crassus, the invasive nematode of anguillid eels. Aquat Invasions. 2012, 7 (4): 443-453. 10.3391/ai.2012.7.4.001.Google Scholar
- Sasal P, Taraschewski H, Valade P, Grondin H, Wielgoss S, Moravec F: Parasite communities in eels of the Island of reunion (indian ocean): a lesson in parasite introduction. Parasitol Res. 2008, 102 (6): 1343-1350. 10.1007/s00436-008-0916-5.PubMedGoogle Scholar
- Wielgoss S, Hollandt F, Wirth T, Meyer A: Genetic signatures in an invasive parasite of Anguilla anguilla correlate with differential stock management. J Fish Biol. 2010, 77 (1): 191-210. 10.1111/j.1095-8649.2010.02670.x.PubMedGoogle Scholar
- Wielgoss S, Sanetra M, Meyer A, Wirth T: Isolation and characterization of short tandem repeats in an invasive swimbladder nematode, parasitic in Atlantic freshwater eels, Anguillicola crassus. Mol Ecol Notes. 2007, 7 (6): 1051-1053. 10.1111/j.1471-8286.2007.01773.x.Google Scholar
- Wielgoss S, Taraschewski H, Meyer A, Wirth T: Population structure of the parasitic nematode Anguillicola crassus, an invader of declining North Atlantic eel stocks. Mol Ecol. 2008, 17 (15): 3478-3495. 10.1111/j.1365-294X.2008.03855.x.PubMedGoogle Scholar
- Wielgoss S: Population genetic aspects of a newly established parasite-host system between the nematode invader Anguillicola crassus, and the North Atlantic freshwater eels, Anguilla sp. Fachbereich Biologie. 2009, Universität Konstanz, Online resource?(KOPS):?http://kops.ub.uni-konstanz.de/volltexte/2009/8134Google Scholar
- Hedgecock D: Temporal and spatial genetic-structure of marine animal populations in the California current. Calif Cooperative Oceanic Fisheries Investig Reports. 1994, 35: 73-81.Google Scholar
- Pujolar JM, Maes GE, Volckaert FAM: Genetic patchiness among recruits in the European eel Anguilla anguilla. Mar Ecol Prog Ser. 2006, 307: 209-217.Google Scholar
- Pujolar JM, De Leo GA, Ciccotti E, Zane L: Genetic composition of Atlantic and Mediterranean recruits of European eel Anguilla anguilla based on EST-linked microsatellite loci. J Fish Biol. 2009, 74 (9): 2034-2046. 10.1111/j.1095-8649.2009.02267.x.PubMedGoogle Scholar
- Pujolar JM, Maes GE, Volckaert FAM: Genetic and morphometric heterogeneity among recruits of the European eel, Anguilla anguilla. Bull Mar Sci. 2007, 81 (2): 297-308.Google Scholar
- Andrello M, Bevacqua D, Maes GE, De Leo GA: An integrated genetic-demographic model to unravel the origin of genetic structure in European eel (Anguilla anguilla L.). Evol Appl. 2011, 4 (4): 517-533. 10.1111/j.1752-4571.2010.00167.x.PubMedPubMed CentralGoogle Scholar
- Gagnaire P-A, Normandeau E, Bernatchez L: Comparative genomics reveals adaptive protein evolution and a possible cytonuclear incompatibility between european and american eels. Mol Biol Evol. 2012, 29 (10): 2909-2919. 10.1093/molbev/mss076.PubMedGoogle Scholar
- Arnold ML: Natural Hybridization and Evolution. 1997, New-York: Oxford University PressGoogle Scholar
- Darwin C: The Origin of Species by Means of Natural Selection. 1859, London, UK: J. MurrayGoogle Scholar
- Dobzhansky T: Genetics and the Origin of Species. 1937, New York, U.S.A.: Columbia Univ. PressGoogle Scholar
- Mayr E: Systematics and the Origin of Species. 1942, New York, U.S.A.: Columbia Univ. PressGoogle Scholar
- Wilson EO: The challenge from related species. The Genetics of Colonizing Species. Edited by: Baker HG, Stebbins GL. 1965, Orlando, USA: Academic Press, 7-24.Google Scholar
- Coyne JA, Orr HA: Speciation. 2004, Sunderland, MA, U.S.A: Sinauer AssociatesGoogle Scholar
- Arnold ML, Meyer A: Natural hybridization in primates: one evolutionary mechanism. Zoology. 2006, 109 (4): 261-276. 10.1016/j.zool.2006.03.006.PubMedGoogle Scholar
- Arnold ML, Bennett BD, Zimmer EA: Natural hybridization between Iris fulva and Iris hexagona - pattern of ribosomal DNA variation. Evolution. 1990, 44 (6): 1512-1521. 10.2307/2409333.Google Scholar
- Arnold ML, Buckner CM, Robinson JJ: Pollen-mediated introgression and hybrid speciation in Louisiana Irises. Proc Natl Acad Sci U S A. 1991, 88 (4): 1398-1402. 10.1073/pnas.88.4.1398.PubMedPubMed CentralGoogle Scholar
- Rieseberg LH, Carter R, Zona S: Molecular tests of the hypothesized hybrid origin of two diploid Helianthus species (Asteraceae). Evolution. 1990, 44 (6): 1498-1511. 10.2307/2409332.Google Scholar
- Rieseberg LH, Vanfossen C, Desrochers AM: Hybrid speciation accompanied by genomic reorganization in wild sunflowers. Nature. 1995, 375 (6529): 313-316. 10.1038/375313a0.Google Scholar
- Xie J, Fu Y, Jiang D, Li G, Huang J, Li B, Hsiang T, Peng Y: Intergeneric transfer of ribosomal genes between two fungi. BMC Evol Biol. 2008, 8: 87-10.1186/1471-2148-8-87.PubMedPubMed CentralGoogle Scholar
- Arnold ML: Natural hybridization as an evolutionary process. Annu Rev Ecol Syst. 1992, 23: 237-261. 10.1146/annurev.es.23.110192.001321.Google Scholar
- Asmussen MA, Arnold J: The effects of admixture and population subdivision on cytonuclear disequilibria. Theor Popul Biol. 1991, 39 (3): 273-300. 10.1016/0040-5809(91)90024-A.PubMedGoogle Scholar
- Kuroki M, Kawai M, Jonsson B, Aoyama J, Miller MJ, Noakes DLG, Tsukamoto K: Inshore migration and otolith microstructure/microchemistry of anguillid glass eels recruited to Iceland. Environ Biol Fishes. 2008, 83 (3): 309-325. 10.1007/s10641-008-9341-y.Google Scholar
- Sokal RR, Rohlf FJ: Biometry. 1995, New York, U.S.A: W. H. Freeman and CompanyGoogle Scholar
- Van Oosterhout C, Hutchinson WF, Wills DPM, Shipley P: MICRO-CHECKER: software for identifying and correcting genotyping errors in microsatellite data. Mol Ecol Notes. 2004, 4 (3): 535-538. 10.1111/j.1471-8286.2004.00684.x.Google Scholar
- Chapuis M-P, Estoup A: Microsatellite null alleles and estimation of population differentiation. Mol Biol Evol. 2007, 24 (3): 621-631.PubMedGoogle Scholar
- Falush D, Stephens M, Pritchard JK: Inference of population structure using multilocus genotype data: dominant markers and null alleles. Mol Ecol Notes. 2007, 7 (4): 574-578. 10.1111/j.1471-8286.2007.01758.x.PubMedPubMed CentralGoogle Scholar
- Hubisz MJ, Falush D, Stephens M, Pritchard JK: Inferring weak population structure with the assistance of sample group information. Mol Ecol Resour. 2009, 9 (5): 1322-1332. 10.1111/j.1755-0998.2009.02591.x.PubMedPubMed CentralGoogle Scholar
- Pritchard JK, Stephens M, Donnelly P: Inference of population structure using multilocus genotype data. Genetics. 2000, 155 (2): 945-959.PubMedPubMed CentralGoogle Scholar
- Falush D, Stephens M, Pritchard JK: Inference of population structure using multilocus genotype data: linked loci and correlated allele frequencies. Genetics. 2003, 164 (4): 1567-1587.PubMedPubMed CentralGoogle Scholar
- Mantel N: The detection of disease clustering and a generalized regression approach. Cancer Res. 1967, 27 (2): 209-220.PubMedGoogle Scholar
- Hey J: Isolation with migration models for more than two populations. Mol Biol Evol. 2010, 27 (4): 905-920. 10.1093/molbev/msp296.PubMedPubMed CentralGoogle Scholar
- Hey J, Nielsen R: Integration within the Felsenstein equation for improved Markov chain Monte Carlo methods in population genetics. Proc Natl Acad Sci U S A. 2007, 104 (8): 2785-2790. 10.1073/pnas.0611164104.PubMedPubMed CentralGoogle Scholar
- Avise JC: Catadromous eels continue to be slippery research subjects. Mol Ecol. 2011, 20 (7): 1317-1319. 10.1111/j.1365-294X.2011.05012.x.PubMedGoogle Scholar
- Arnold J: Cytonuclear disequilibria in hybrid zones. Annu Rev Ecol Syst. 1993, 24: 521-554. 10.1146/annurev.es.24.110193.002513.Google Scholar
- Gagnaire PA, Minegishi Y, Zenboudji S, Valade P, Aoyama J, Berrebi P: Within-population structure highlighted by differential introgression across semipermeable barriers to gene flow in Anguilla marmorata. Evolution. 2011, 65 (12): 3413-3427. 10.1111/j.1558-5646.2011.01404.x.PubMedGoogle Scholar
- Tsukamoto K: Spawning of eels near a seamount. Nature. 2006, 439 (7079): 929-929. 10.1038/439929a.PubMedGoogle Scholar
- McCleave JD: Contrasts between spawning times of Anguilla species estimated from larval sampling at sea and from otolith analysis of recruiting glass eels. Mar Biol. 2008, 155 (3): 249-262. 10.1007/s00227-008-1026-8.Google Scholar
- Munk P, Hansen MM, Maes GE, Nielsen TG, Castonguay M, Riemann L, Sparholt H, Als TD, Aarestrup K, Andersen NG, Bachler M: Oceanic fronts in the Sargasso Sea control the early life and drift of Atlantic eels. Proc R Soc B-Biol Sci. 2010, 277 (1700): 3593-3599. 10.1098/rspb.2010.0900.Google Scholar
- Kettle AJ, Haines K: How does the European eel (Anguilla anguilla) retain its population structure during its larval migration across the North Atlantic Ocean?. Can J Fish Aquat Sci. 2006, 63 (1): 90-106. 10.1139/f05-198.Google Scholar
- Frankowski J, Jennerich S, Schaarschmidt T, Ubl C, Juerss K, Bastrop R: Validation of the occurrence of the American eel Anguilla rostrata (Le Sueur, 1817) in free-draining European inland waters. Biol Invasions. 2009, 11 (6): 1301-1309. 10.1007/s10530-008-9337-8.Google Scholar
- Okamura A, Zhang H, Mikawa N, Kotake A, Yamada Y, Utoh T, Horie N, Tanaka S, Oka HP, Tsukamoto K: Decline in non-native freshwater eels in Japan: ecology and future perspectives. Environ Biol Fishes. 2008, 81 (3): 347-358. 10.1007/s10641-007-9205-x.Google Scholar
- Arai T, Otake T, Tsukamoto K: iming of metamorphosis and larval segregation of the Atlantic eels Anguilla rostrata and A. anguilla, as revealed by otolith microstructure and microchemistry. Mar Biol. 2000, 137: 39-45. 10.1007/s002270000326.Google Scholar
- Palumbi SR: Genetic divergence, reproductive isolation, and marine speciation. Annu Rev Ecol Syst. 1994, 25: 547-572. 10.1146/annurev.es.25.110194.002555.Google Scholar
- Lecomte-Finiger R: The early life of the European eel. Nature. 1994, 370 (6489): 424-424. 10.1038/370424a0.Google Scholar
- Schmidt J: The breeding places of the eel. Philos Trans R Soc B Biol Sci. 1923, 211: 179-U138. 10.1098/rstb.1923.0004.Google Scholar
- Boetius J, Harding EF: A re-examination of Johannes Schmidt ’s Atlantic eel investigations. Dana. 1985, 4: 129-162.Google Scholar
- Wang CH, Tzeng WN: Interpretation of geographic variation in size of American eel Anguilla rostrata elvers on the Atlantic coast of North America using their life history and otolith ageing. Mar Ecol Prog Ser. 1998, 168: 35-43.Google Scholar
- Wang CH, Tzeng WN: The timing of metamorphosis and growth rates of American and European eel leptocephali: a mechanism of larval segregative migration. Fish Res. 2000, 46 (1–3): 191-205.Google Scholar
- Bonhommeau S, Castonguay M, Rivot E, Sabatie R, Le Pape O: The duration of migration of Atlantic Anguilla larvae. Fish Fish. 2010, 11 (3): 289-306. 10.1111/j.1467-2979.2010.00362.x.Google Scholar
- Boetius J: Greenland eels, Anguilla rostrata LeSueur. Dana. 1985, 4: 41-48.Google Scholar
- COM: Establishing measures for the recovery of the stock of European eel. Council Regulation (EC) No 1100/2007. Official J European Union L. 2007, 248: 17-23.Google Scholar
- Dekker W: A Conceptual Management Framework for the Restoration of the Declining European eel Stock. Eels at the Edge: Science, Status, and Conservation Concerns. Edited by: Casselman JM, Cairns DK. 2009, Bethesda, Maryland: American Fisheries Society Symposium 58, 3-19.Google Scholar
- Raymond M, Rousset F: GENEPOP (version 1.2): population genetics software for exact tests and ecumenicism. J Hered. 1995, 86 (3): 248-249.Google Scholar
- Dawson KJ, Belkhir K: A Bayesian approach to the identification of panmictic populations and the assignment of individuals. Genet Res. 2001, 78 (1): 59-77.PubMedGoogle Scholar
- Kalinowski ST: HP-RARE 1.0: a computer program for performing rarefaction on measures of allelic richness. Mol Ecol Notes. 2005, 5 (1): 187-189. 10.1111/j.1471-8286.2004.00845.x.Google Scholar
- Excoffier L, Laval G, Schneider S: Arlequin (version 3.0): an integrated software package for population genetics data analysis. Evol Bioinforma. 2005, 1: 47-50.Google Scholar
- Tagliavini J, Harrison IJ, Gandolfi G: Discrimination between Anguilla anguilla and A. rostrata by polymerase chain reaction-restriction fragment length polymorphism analysis. J Fish Biol. 1995, 47 (4): 741-743.Google Scholar
- Clement M, Posada D, Crandall KA: TCS: a computer program to estimate gene genealogies. Mol Ecol. 2000, 9 (10): 1657-1659. 10.1046/j.1365-294x.2000.01020.x.PubMedGoogle Scholar
- Keane TM, Creevey CJ, Pentony MM, Naughton TJ, McInerney JO: Assessment of methods for amino acid matrix selection and their use on empirical data shows that ad hoc assumptions for choice of matrix are not justified. BMC Evol Biol. 2006, 6: 29-10.1186/1471-2148-6-29.PubMedPubMed CentralGoogle Scholar
- Ocean color browser. http://oceancolor.gsfc.nasa.gov/cgi/browse.pl,
- Evanno G, Regnaut S, Goudet J: Detecting the number of clusters of individuals using the software STRUCTURE: a simulation study. Mol Ecol. 2005, 14 (8): 2611-2620. 10.1111/j.1365-294X.2005.02553.x.PubMedGoogle Scholar
- Vaha JP, Primmer CR: Efficiency of model-based Bayesian methods for detecting hybrid individuals under different hybridization scenarios and with different numbers of loci. Mol Ecol. 2006, 15 (1): 63-72.PubMedGoogle Scholar
- Cavallis L, Edwards AWF: Phylogenetic analysis: models and estimation procedures. Am J Hum Genet. 1967, 19 (3P1): 233-Google Scholar
- Felsenstein J: PHYLIP: phylogenetic inference package version 3.68. Distributed over the Internet. 2008, http://evolution.genetics.washington.edu/phylip.html,Google Scholar
- Jombart T: adegenet: a R package for the multivariate analysis of genetic markers. Bioinformatics. 2008, 24 (11): 1403-1405. 10.1093/bioinformatics/btn129.PubMedGoogle Scholar
- Jombart T, Ahmed I: adegenet 1.3-1: new tools for the analysis of genome-wide SNP data. Bioinformatics. 2011, 27 (21): 3070-3071. 10.1093/bioinformatics/btr521.PubMedPubMed CentralGoogle Scholar
- Dray S, Dufour AB: The ade4 package: implementing the duality diagram for ecologists. J Stat Softw. 2007, 22 (4): 1-20.Google Scholar
- Paradis E, Claude J, Strimmer K: APE: analyses of phylogenetics and evolution in R language. Bioinformatics. 2004, 20 (2): 289-290. 10.1093/bioinformatics/btg412.PubMedGoogle Scholar
- Popescu AA, Huber KT, Paradis E: ape 3.0: New tools for distance-based phylogenetics and evolutionary analysis in R. Bioinformatics. 2012, 28 (11): 1536-1537. 10.1093/bioinformatics/bts184.PubMedGoogle Scholar
- Paradis E: pegas: an R package for population genetics with an integrated-modular approach. Bioinformatics. 2010, 26 (3): 419-420. 10.1093/bioinformatics/btp696.PubMedGoogle Scholar
- Wickham H: The split-apply-combine strategy for data analysis. J Stat Softw. 2011, 40 (1): 1-29.Google Scholar
- Hey J: The divergence of chimpanzee species and subspecies as revealed in multipopulation isolation-with-migration analyses. Mol Biol Evol. 2010, 27 (4): 921-933. 10.1093/molbev/msp298.PubMedPubMed CentralGoogle Scholar
This article is published under license to BioMed Central Ltd. This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/2.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly credited.