Magnus W. Jacobsen

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 64oN, we examined the patterns of genome‐wide genetic diversity across locations. We tested for local selection by searching for increased population differentiation using FST‐based outlier tests and by testing for significant associations between allele frequencies and environmental variables. The overall low genetic differentiation found (FST = 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.

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.

Mitogenome sequencing reveals shallow evolutionary histories and recent divergence time between morphologically and ecologically distinct European whitefish (Coregonus spp.)

The advent of second‐generation sequencing has made it possible to quickly and economically generate whole mitochondrial genome (mitogenome) sequences. To date, mitogenome studies of nonmodel organisms have demonstrated increased power for resolving interspecies relationships. We explored an alternate use of such data to recover relationships and population history of closely related lineages with a shallow evolutionary history. Using a GS‐FLX platform, we sequenced 106 mitogenomes from the Coregonus lavaretus (Europe) and Coregonus clupeaformis (North America) species complexes to investigate the evolutionary history of the endangered Danish North Sea houting (NSH) and other closely related Danish and Baltic European lake whitefish (ELW). Two well‐supported clades were found within both ELW and NSH, probably reflecting historical introgression via Baltic migrants. Although ELW and NSH are not reciprocally monophyletic, they share no haplotypes, suggesting recent, but strong, reproductive isolation. The divergence time between NSH and the geographically closest ELW population was estimated using IMa, assuming isolation with migration and a new mutation rate estimate chosen to avoid time‐dependency effects. The estimate of c. 2700 bp was remarkably similar to results obtained using microsatellite markers. Within North American C. clupeaformis, the divergence time between the two lineages (Atlantic and Acadian) was estimated as between 20 000 and 60 000 bp. Under the assumption that NSH and ELW colonized Denmark following the last glacial maximum, Bayesian Serial SimCoal analysis showed consistency with a scenario of long‐term stability, resulting from a rapid initial sixfold population expansion. The findings illustrate the utility of mitogenome data for resolving recent intraspecific divergence events and provide evidence for recent reproductive isolation of the phenotypically divergent NSH.

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.

Genomic footprints of speciation in Atlantic eels (Anguilla anguilla and A. rostrata)

The importance of speciation‐with‐geneflow scenarios is increasingly appreciated. However, the specific processes and the resulting genomic footprints of selection are subject to much discussion. We studied the genomics of speciation between the two panmictic, sympatrically spawning sister species; European (Anguilla anguilla) and American eel (A. rostrata). Divergence is assumed to have initiated more than 3 Ma, and although low gene flow still occurs, strong postzygotic barriers are present. Restriction‐site‐associated DNA (RAD) sequencing identified 328 300 SNPs for subsequent analysis. However, despite the presence of 3757 strongly differentiated SNPs (FST > 0.8), sliding window analyses of FST showed no larger genomic regions (i.e. hundreds of thousands to millions of bases) of elevated differentiation. Overall FST was 0.041, and linkage disequilibrium was virtually absent for SNPs separated by more than 1000 bp. We suggest this to reflect a case of genomic hitchhiking, where multiple regions are under directional selection between the species. However, low but biologically significant gene flow and high effective population sizes leading to very low genetic drift preclude accumulation of strong background differentiation. Genes containing candidate SNPs for positive selection showed significant enrichment for gene ontology (GO) terms relating to developmental processes and phosphorylation, which seems consistent with assumptions that differences in larval phase duration and migratory distances underlie speciation. Most SNPs under putative selection were found outside coding regions, lending support to emerging views that noncoding regions may be more functionally important than previously assumed. In total, the results demonstrate the necessity of interpreting genomic footprints of selection in the context of demographic parameters and life‐history features of the studied species.

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.

Comparative analysis of complete mitochondrial genomes suggests that relaxed purifying selection is driving high nonsynonymous evolutionary rate of the NADH2 gene in whitefish (Coregonus ssp.)

Several studies have recently reported evidence for positive selection acting on the mitochondrial genome (mitogenome), emphasizing its potential role in adaptive divergence and speciation. In this study we searched 107 full mitogenomes of recently diverged species and lineages of whitefish (Coregonus ssp.) for signals of positive selection. These salmonids show several distinct morphological and ecological differences that may be associated with energetics and therefore potentially positive selection at the mitogenome level. We found that purifying selection and genetic drift were the predominant evolutionary forces acting on the analyzed mitogenomes. However, the NADH dehydrogenase 2 gene (ND2) showed a highly elevated dN/dS ratio compared to the other mitochondrial genes, which was significantly higher in whitefish compared to other salmonids. We therefore further examined nonsynonymous evolution in ND2 by (i) mapping amino acid changes to a protein model structure which showed that they were located away from key functional residues of the protein, (ii) locating them in the sequences of other species of fish (Salmonidae, Anguillidae, Scombridae and Percidae) only to find pronounced overlap of nonsynonymous regions. We thus conclude that relaxed purifying selection is driving the evolution of ND2 by affecting mostly regions that have lower functional relevance.

Do North Atlantic eels show parallel patterns of spatially varying selection

BackgroundThe two North Atlantic eel species, the European and the American eel, represent an ideal system in which to study parallel selection patterns due to their sister species status and the presence of ongoing gene flow. A panel of 80 coding-gene SNPs previously analyzed in American eel was used to genotype European eel individuals (glass eels) from 8 sampling locations across the species distribution. We tested for single-generation signatures of spatially varying selection in European eel by searching for elevated genetic differentiation using FST-based outlier tests and by testing for significant associations between allele frequencies and environmental variables.ResultsWe found signatures of possible selection at a total of 11 coding-gene SNPs. Candidate genes for local selection constituted mainly genes with a major role in metabolism as well as defense genes. Contrary to what has been found for American eel, only 2 SNPs in our study correlated with differences in temperature, which suggests that other explanatory variables may play a role. None of the genes found to be associated with explanatory variables in European eel showed any correlations with environmental factors in the previous study in American eel.ConclusionsThe different signatures of selection between species could be due to distinct selective pressures associated with the much longer larval migration for European eel relative to American eel. The lack of parallel selection in North Atlantic eels could also be due to most phenotypic traits being polygenic, thus reducing the likelihood of selection acting on the same genes in both species.

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.

Relationship between amino acid changes in mitochondrial ATP6 and life-history variation in anguillid eels

Mitochondrial genes are part of the oxidative phosphorylation pathway and important for energy production. Although evidence for positive selection at the mitochondrial level exists, few studies have investigated the link between amino acid changes and phenotype. Here we test the hypothesis that differences in two life-history related traits, migratory distance between spawning and foraging areas and larval phase duration, are associated with divergent selection within the mitochondrial ATP6 gene in anguillid eels. We compare amino acid changes among 18 species with the sequence of the putative ancestral species, believed to have shown short migratory distance and larval phase duration. We find positive correlations between both life-history related traits and (i) the number of amino acid changes and (ii) the strength of the combined physico-chemical and structural changes at positions previously identified as candidates for positive selection. This supports a link between genotype and phenotype driven by positive selection at ATP6.