Atlantic salmon (Salmo salar) is one of the most extensively studied fish species in the world due to its significance in aquaculture, fisheries and ongoing conservation efforts to protect declining populations. Yet, limited genomic resources have hampered our understanding of genetic architecture in the species and the genetic basis of adaptation to the wide range of natural and artificial environments it occupies. In this study, we describe the development of a medium-density Atlantic salmon single nucleotide polymorphism (SNP) array based on expressed sequence tags (ESTs) and genomic sequencing. The array was used in the most extensive assessment of population genetic structure performed to date in this species. A total of 6176 informative SNPs were successfully genotyped in 38 anadromous and freshwater wild populations distributed across the species natural range. Principal component analysis clearly differentiated European and North American populations, and within Europe, three major regional genetic groups were identified for the first time in a single analysis. We assessed the potential for the array to disentangle neutral and putative adaptive divergence of SNP allele frequencies across populations and among regional groups. In Europe, secondary contact zones were identified between major clusters where endogenous and exogenous barriers could be associated, rendering the interpretation of environmental influence on potentially adaptive divergence equivocal. A small number of markers highly divergent in allele frequencies (outliers) were observed between (multiple) freshwater and anadromous populations, between northern and southern latitudes, and when comparing Baltic populations to all others. We also discuss the potential future applications of the SNP array for conservation, management and aquaculture.

Anadromous Atlantic salmon (Salmo salar) is a species of major conservation and management concern in North America, where population abundance has been declining over the past 30 years. Effective conservation actions require the delineation of conservation units to appropriately reflect the spatial scale of intraspecific variation and local adaptation. Towards this goal, we used the most comprehensive genetic and genomic database for Atlantic salmon to date, covering the entire North American range of the species. The database included microsatellite data from 9142 individuals from 149 sampling locations and data from a medium-density SNP array providing genotypes for >3000 SNPs for 50 sampling locations. We used neutral and putatively selected loci to integrate adaptive information in the definition of conservation units. Bayesian clustering with the microsatellite data set and with neutral SNPs identified regional groupings largely consistent with previously published regional assessments. The use of outlier SNPs did not result in major differences in the regional groupings, suggesting that neutral markers can reflect the geographic scale of local adaptation despite not being under selection. We also performed assignment tests to compare power obtained from microsatellites, neutral SNPs and outlier SNPs. Using SNP data substantially improved power compared to microsatellites, and an assignment success of 97% to the population of origin and of 100% to the region of origin was achieved when all SNP loci were used. Using outlier SNPs only resulted in minor improvements to assignment success to the population of origin but improved regional assignment. We discuss the implications of these new genetic resources for the conservation and management of Atlantic salmon in North America.

Identification of discrete and unique assemblages of individuals or populations is central to the management of exploited species. Advances in population genomics provide new opportunities for re-evaluating existing conservation units but comparisons among approaches remain rare. We compare the utility of RAD-seq, a single nucleotide polymorphism (SNP) array and a microsatellite panel to resolve spatial structuring under a scenario of possible trans-Atlantic secondary contact in a threatened Atlantic Salmon, Salmo salar, population in southern Newfoundland. Bayesian clustering indentified two large groups subdividing the existing conservation unit and multivariate analyses indicated significant similarity in spatial structuring among the three data sets. mtDNA alleles diagnostic for European ancestry displayed increased frequency in southeastern Newfoundland and were correlated with spatial structure in all marker types. Evidence consistent with introgression among these two groups was present in both SNP data sets but not the microsatellite data. Asymmetry in the degree of introgression was also apparent in SNP data sets with evidence of gene flow towards the east or European type. This work highlights the utility of RAD-seq based approaches for the resolution of complex spatial patterns, resolves a region of trans-Atlantic secondary contact in Atlantic Salmon in Newfoundland and demonstrates the utility of multiple marker comparisons in identifying dynamics of introgression.

In some wild Atlantic salmon populations, rapid declines in numbers of wild returning adults has been associated with an increase in the prevalence of farmed salmon. Studies of phenotypic variation have shown that interbreeding between farmed and wild salmon may lead to loss of local adaptation. Yet, few studies have attempted to assess the impact of interbreeding at the genome level, especially among North American populations. Here, we document temporal changes in the genetic makeup of the severely threatened Magaguadavic River salmon population (Bay of Fundy, Canada), a population that might have been impacted by interbreeding with farmed salmon for nearly 20 years. Wild and farmed individuals caught entering the river from 1980 to 2005 were genotyped at 112 single-nucleotide polymorphisms (SNPs), and/or eight microsatellite loci, to scan for potential shifts in adaptive genetic variation. No significant temporal change in microsatellite-based estimates of allele richness or gene diversity was detected in the wild population, despite its precipitous decline in numbers over the last two decades. This might reflect the effect of introgression from farmed salmon, which was corroborated by temporal change in linkage-disequilibrium. Moreover, SNP genome scans identified a temporal decrease in candidate loci potentially under directional selection. Of particular interest was a SNP previously shown to be strongly associated with an important quantitative trait locus for parr mark number, which retained its genetic distinctiveness between farmed and wild fish longer than other outliers. Overall, these results indicate that farmed escapees have introgressed with wild Magaguadavic salmon resulting in significant alteration of the genetic integrity of the native population, including possible loss of adaptation to wild conditions.

Wild populations of Atlantic salmon have declined worldwide. While the causes for this decline may be complex and numerous, increased mortality at sea is predicted to be one of the major contributing factors. Examining the potential changes occurring in the genome-wide composition of populations during this migration has the potential to tease apart some of the factors influencing marine mortality. Here, we genotyped 5568 SNPs in Atlantic salmon populations representing two distinct regional genetic groups and across two cohorts to test for differential allelic and genotypic frequencies between juveniles (smolts) migrating to sea and adults (grilses) returning to freshwater after 1 year at sea. Given the complexity of the traits potentially associated with sea mortality, we contrasted the outcomes of a single-locus F(ST) based genome scan method with a new multilocus framework to test for genetically based differential mortality at sea. While numerous outliers were identified by the single-locus analysis, no evidence for parallel, temporally repeated selection was found. In contrast, the multilocus approach detected repeated patterns of selection for a multilocus group of 34 covarying SNPs in one of the two populations. No significant pattern of selective mortality was detected in the other population, suggesting different causes of mortality among populations. These results first support the hypothesis that selection mainly causes small changes in allele frequencies among many covarying loci rather than a small number of changes in loci with large effects. They also point out that moving away from the a strict 'selective sweep paradigm' towards a multilocus genetics framework may be a more useful approach for studying the genomic signatures of natural selection on complex traits in wild populations.

Little is known about the genetic basis differentiating resident and anadromous forms found in many salmonid species. Using a medium-density SNP array, we documented genomic diversity and divergence at 2336 genetically mapped loci among three pairs of North American anadromous and freshwater Atlantic salmon populations. Our results show that across the genome, freshwater populations have lower diversity and a smaller proportion of private polymorphism relative to anadromous populations. Moreover, differentiation was more pronounced among freshwater than among anadromous populations at multiple spatial scales, suggesting a large effect of genetic drift in these isolated freshwater populations. Using nonhierarchical and hierarchical genome scans, we identified hundreds of markers spread across the genome that are potentially under divergent selection between anadromous and freshwater populations, but few outlier loci were repeatedly found in all three freshwater–anadromous comparisons. Similarly, a sliding window analysis revealed numerous regions of high divergence that were nonparallel among the three comparisons. These last results show little evidence for the parallel evolution of alleles selected for in freshwater populations, but suggest nonparallel adaptive divergence at many loci of small effects distributed through the genome. Overall, this study emphasizes the important role of genetic drift in driving genome-wide reduction in diversity and divergence in freshwater Atlantic salmon populations and suggests a complex multigenic basis of adaptation to resident and anadromous strategies with little parallelism.

We performed population genetic analyses on the American eel (Anguilla rostrata) with three main objectives. First, we conducted the most comprehensive analysis of neutral genetic population structure to date to revisit the null hypothesis of panmixia in this species. Second, we used this data to provide the first estimates of contemporary effective population size (Ne ) and to document temporal variation in effective number of breeders (Nb ) in American eel. Third, we tested for statistical associations between temporal variation in the North Atlantic Oscillation (NAO), the effective number of breeders and two indices of recruit abundance. A total of 2142 eels from 32 sampling locations were genotyped with 18 microsatellite loci. All measures of differentiation were essentially zero, and no evidence for significant spatial or temporal genetic differentiation was found. The panmixia hypothesis should thus be accepted for this species. Nb estimates varied by a factor of 23 among 12 cohorts, from 473 to 10,999. The effective population size Ne was estimated at 10,532 (95% CI, 9312-11,752). This study also showed that genetically based demographic indices, namely Nb and allelic richness (Ar), can be used as surrogates for the abundance of breeders and recruits, which were both shown to be positively influenced by variation during high (positive) NAO phases. Thus, long-term genetic monitoring of American glass eels at several sites along the North American Atlantic coast would represent a powerful and efficient complement to census monitoring to track demographic fluctuations and better understand their causes.

Fisheries targeting mixtures of populations risk the overutilization of minor stock constituents unless harvests are monitored and managed. We evaluated stock composition and exploitation of Atlantic salmon (Salmo salar) in a subsistence fishery in coastal Labrador, Canada, using genetic mixture analysis and individual assignment with a microsatellite baseline (15 loci, 11 829 individuals, 12 regional groups) encompassing the species’ western Atlantic range. Bayesian and maximum likelihood mixture analyses of fishery samples over 6 years (2006–2011; 1772 individuals) indicate contributions of adjacent stocks of 96%–97%. Estimates of fishery-associated exploitation were highest for Labrador salmon (4.2%–10.6% per year) and generally <1% for other regions. Individual assignment of fishery samples indicated nonlocal contributions to the fishery (e.g., Quebec, Newfoundland) were rare and primarily in southern Labrador, consistent with migration pathways utilizing the Strait of Belle Isle. This work illustrates how genetic analysis of mixed stock Atlantic salmon fisheries in the Northwest Atlantic using this new baseline can disentangle exploitation and reveal complex migratory behaviours.