Assessing pre- and post-zygotic barriers between North Atlantic eels (Anguilla anguilla and A. rostrata).

Elucidating barriers to gene flow is important for understanding the dynamics of speciation. Here we investigate pre- and post-zygotic mechanisms acting between the two hybridizing species of Atlantic eels: Anguilla anguilla and A. rostrata. Temporally varying hybridization was examined by analyzing 85 species-diagnostic single-nucleotide polymorphisms (SNPs; FST ⩾0.95) in eel larvae sampled in the spawning region in the Sargasso Sea in 2007 (N=92) and 2014 (N=460). We further investigated whether genotypes at these SNPs were nonrandomly distributed in post-F1 hybrids, indicating selection. Finally, we sequenced the mitochondrial ATP6 and nuclear ATP5c1 genes in 19 hybrids, identified using SNP and restriction site associated DNA (RAD) sequencing data, to test a previously proposed hypothesis of cytonuclear incompatibility leading to adenosine triphosphate (ATP) synthase dysfunction and selection against hybrids. No F1 hybrids but only later backcrosses were observed in the Sargasso Sea in 2007 and 2014. This suggests that interbreeding between the two species only occurs in some years, possibly controlled by environmental conditions at the spawning grounds, or that interbreeding has diminished through time as a result of a declining number of spawners. Moreover, potential selection was found at the nuclear and the cytonuclear levels. Nonetheless, one glass eel individual showed a mismatch, involving an American ATP6 haplotype and European ATP5c1 alleles. This contradicted the presence of cytonuclear incompatibility but may be explained by that (1) cytonuclear incompatibility is incomplete, (2) selection acts at a later life stage or (3) other genes are important for protein function. In total, the study demonstrates the utility of genomic data when examining pre- and post-zyotic barriers in natural hybrids.

Signatures of natural selection between life cycle stages separated by metamorphosis in European eel.

BACKGROUND: Species showing complex life cycles provide excellent opportunities to study the genetic associations between life cycle stages, as selective pressures may differ before and after metamorphosis. The European eel presents a complex life cycle with two metamorphoses, a first metamorphosis from larvae into glass eels (juvenile stage) and a second metamorphosis into silver eels (adult stage). We tested the hypothesis that different genes and gene pathways will be under selection at different life stages when comparing the genetic associations between glass eels and silver eels. RESULTS: We used two sets of markers to test for selection: first, we genotyped individuals using a panel of 80 coding-gene single nucleotide polymorphisms (SNPs) developed in American eel; second, we investigated selection at the genome level using a total of 153,423 RAD-sequencing generated SNPs widely distributed across the genome. Using the RAD approach, outlier tests identified a total of 2413 (1.57%) potentially selected SNPs. Functional annotation analysis identified signal transduction pathways as the most over-represented group of genes, including MAPK/Erk signalling, calcium signalling and GnRH (gonadotropin-releasing hormone) signalling. Many of the over-represented pathways were related to growth, while others could result from the different conditions that eels inhabit during their life cycle. CONCLUSIONS: The observation of different genes and gene pathways under selection when comparing glass eels vs. silver eels supports the adaptive decoupling hypothesis for the benefits of metamorphosis. Partitioning the life cycle into discrete morphological phases may be overall beneficial since it allows the different life stages to respond independently to their unique selection pressures. This might translate into a more effective use of food and niche resources and/or performance of phase-specific tasks (e.g. feeding in the case of glass eels, migrating and reproducing in the case of silver eels).

Speciation and demographic history of Atlantic eels (Anguilla anguilla and A. rostrata) revealed by mitogenome sequencing.

Processes leading to speciation in oceanic environments without obvious physical barriers remain poorly known. European and American eel (Anguilla anguilla and A. rostrata) spawn in partial sympatry in the Sargasso Sea. Larvae are advected by the Gulf Stream and other currents towards the European/North African and North American coasts, respectively. We analyzed 104 mitogenomes from the two species along with mitogenomes of other Anguilla and outgroup species. We estimated divergence time between the two species to identify major events involved in speciation. We also considered two previously stated hypotheses: one where the ancestral species was present in only one continent but was advected across the Atlantic by ocean current changes and another where population declines during Pleistocene glaciations led to increasing vicariance, facilitating speciation. Divergence time was estimated to ∼3.38 Mya, coinciding with the closure of the Panama Gateway that led to reinforcement of the Gulf Stream. This could have advected larvae towards European/North African coasts, in which case American eel would be expected to be the ancestral species. This scenario could, however, not be unequivocally confirmed by analyses of dN/dS, nucleotide diversity and effective population size estimates. Extended bayesian skyline plots showed fluctuations of effective population sizes and declines during glaciations, and thus also lending support to the importance of vicariance during speciation. There was evidence for positive selection at the ATP6 and possibly ND5 genes, indicating a role in speciation. The findings suggest an important role of ocean current changes in speciation of marine organisms.

Genome-wide single-generation signatures of local selection in the panmictic European eel.

Next-generation sequencing and the collection of genome-wide data allow identifying adaptive variation and footprints of directional selection. Using a large SNP data set from 259 RAD-sequenced European eel individuals (glass eels) from eight locations between 34 and 64(o) N, we examined the patterns of genome-wide genetic diversity across locations. We tested for local selection by searching for increased population differentiation using F(ST) -based outlier tests and by testing for significant associations between allele frequencies and environmental variables. The overall low genetic differentiation found (F(ST) = 0.0007) indicates that most of the genome is homogenized by gene flow, providing further evidence for genomic panmixia in the European eel. The lack of genetic substructuring was consistent at both nuclear and mitochondrial SNPs. Using an extensive number of diagnostic SNPs, results showed a low occurrence of hybrids between European and American eel, mainly limited to Iceland (5.9%), although individuals with signatures of introgression several generations back in time were found in mainland Europe. Despite panmixia, a small set of SNPs showed high genetic differentiation consistent with single-generation signatures of spatially varying selection acting on glass eels. After screening 50 354 SNPs, a total of 754 potentially locally selected SNPs were identified. Candidate genes for local selection constituted a wide array of functions, including calcium signalling, neuroactive ligand-receptor interaction and circadian rhythm. Remarkably, one of the candidate genes identified is PERIOD, possibly related to differences in local photoperiod associated with the >30° difference in latitude between locations. Genes under selection were spread across the genome, and there were no large regions of increased differentiation as expected when selection occurs within just a single generation due to panmixia. This supports the conclusion that most of the genome is homogenized by gene flow that removes any effects of diversifying selection from each new generation.

Assessing patterns of hybridization between North Atlantic eels using diagnostic single-nucleotide polymorphisms.

The two North Atlantic eel species, the European eel (Anguilla anguilla) and the American eel (Anguilla rostrata), spawn in partial sympatry in the Sargasso Sea, providing ample opportunity to interbreed. In this study, we used a RAD (Restriction site Associated DNA) sequencing approach to identify species-specific diagnostic single-nucleotide polymorphisms (SNPs) and design a low-density array that combined with screening of a diagnostic mitochondrial DNA marker. Eels from Iceland (N=159) and from the neighboring Faroe Islands (N=29) were genotyped, along with 94 larvae (49 European and 45 American eel) collected in the Sargasso Sea. Our SNP survey showed that the majority of Icelandic eels are pure European eels but there is also an important contribution of individuals of admixed ancestry (10.7%). Although most of the hybrids were identified as F1 hybrids from European eel female × American eel male crosses, backcrosses were also detected, including a first-generation backcross (F1 hybrid × pure European eel) and three individuals identified as second-generation backcrosses originating from American eel × F1 hybrid backcrosses interbreeding with pure European eels. In comparison, no hybrids were observed in the Faroe Islands, the closest bodies of land to Iceland. It is possible that hybrids show an intermediate migratory behaviour between the two parental species that ultimately brings hybrid larvae to the shores of Iceland, situated roughly halfway between the Sargasso Sea and Europe. Only two hybrids were observed among Sargasso Sea larvae, both backcrosses, but no F1 hybrids, that points to temporal variation in the occurrence of hybridization.

A resource of genome-wide single-nucleotide polymorphisms generated by RAD tag sequencing in the critically endangered European eel.

Reduced representation genome sequencing such as restriction-site-associated DNA (RAD) sequencing is finding increased use to identify and genotype large numbers of single-nucleotide polymorphisms (SNPs) in model and nonmodel species. We generated a unique resource of novel SNP markers for the European eel using the RAD sequencing approach that was simultaneously identified and scored in a genome-wide scan of 30 individuals. Whereas genomic resources are increasingly becoming available for this species, including the recent release of a draft genome, no genome-wide set of SNP markers was available until now. The generated SNPs were widely distributed across the eel genome, aligning to 4779 different contigs and 19,703 different scaffolds. Significant variation was identified, with an average nucleotide diversity of 0.00529 across individuals. Results varied widely across the genome, ranging from 0.00048 to 0.00737 per locus. Based on the average nucleotide diversity across all loci, long-term effective population size was estimated to range between 132,000 and 1,320,000, which is much higher than previous estimates based on microsatellite loci. The generated SNP resource consisting of 82,425 loci and 376,918 associated SNPs provides a valuable tool for future population genetics and genomics studies and allows for targeting specific genes and particularly interesting regions of the eel genome.