Environmental extremes drive population structure at the northern range limit of Atlantic salmon in North America.

Conservation of exploited species requires an understanding of both genetic diversity and the dominant structuring forces, particularly near range limits, where climatic variation can drive rapid expansions or contractions of geographic range. Here, we examine population structure and landscape associations in Atlantic salmon (Salmo salar) across a heterogeneous landscape near the northern range limit in Labrador, Canada. Analysis of two amplicon-based data sets containing 101 microsatellites and 376 single nucleotide polymorphisms (SNPs) from 35 locations revealed clear differentiation between populations spawning in rivers flowing into a large marine embayment (Lake Melville) compared to coastal populations. The mechanisms influencing the differentiation of embayment populations were investigated using both multivariate and machine-learning landscape genetic approaches. We identified temperature as the strongest correlate with genetic structure, particularly warm temperature extremes and wider annual temperature ranges. The genomic basis of this divergence was further explored using a subset of locations (n = 17) and a 220K SNP array. SNPs associated with spatial structuring and temperature mapped to a diverse set of genes and molecular pathways, including regulation of gene expression, immune response, and cell development and differentiation. The results spanning molecular marker types and both novel and established methods clearly show climate-associated, fine-scale population structure across an environmental gradient in Atlantic salmon near its range limit in North America, highlighting valuable approaches for predicting population responses to climate change and managing species sustainability.

Transatlantic secondary contact in Atlantic Salmon, comparing microsatellites, a single nucleotide polymorphism array and restriction-site associated DNA sequencing for the resolution of complex spatial structure.

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.

Variable parallelism in the genomic basis of age at maturity across spatial scales in Atlantic Salmon.

Complex traits often exhibit complex underlying genetic architectures resulting from a combination of evolution from standing variation, hard and soft sweeps, and alleles of varying effect size. Increasingly, studies implicate both large-effect loci and polygenic patterns underpinning adaptation, but the extent that common genetic architectures are utilized during repeated adaptation is not well understood. Sea age or age at maturation represents a significant life history trait in Atlantic Salmon (Salmo salar), the genetic basis of which has been studied extensively in European Atlantic populations, with repeated identification of large-effect loci. However, the genetic basis of sea age within North American Atlantic Salmon populations remains unclear, as does the potential for a parallel trans-Atlantic genomic basis to sea age. Here, we used a large single-nucleotide polymorphism (SNP) array and low-coverage whole-genome resequencing to explore the genomic basis of sea age variation in North American Atlantic Salmon. We found significant associations at the gene and SNP level with a large-effect locus (vgll3) previously identified in European populations, indicating genetic parallelism, but found that this pattern varied based on both sex and geographic region. We also identified nonrepeated sets of highly predictive loci associated with sea age among populations and sexes within North America, indicating polygenicity and low rates of genomic parallelism. Despite low genome-wide parallelism, we uncovered a set of conserved molecular pathways associated with sea age that were consistently enriched among comparisons, including calcium signaling, MapK signaling, focal adhesion, and phosphatidylinositol signaling. Together, our results indicate parallelism of the molecular basis of sea age in North American Atlantic Salmon across large-effect genes and molecular pathways despite population-specific patterns of polygenicity. These findings reveal roles for both contingency and repeated adaptation at the molecular level in the evolution of life history variation.

Body length changes for Atlantic salmon (Salmo salar) over five decades exhibit weak spatial synchrony over a broad latitudinal gradient.

Understanding the factors that drive spatial synchrony among populations or species is important for management and recovery of populations. The range-wide declines in Atlantic salmon (Salmo salar) populations may be the result of broad-scale changes in the marine environment. Salmon undergo rapid growth in the ocean; therefore changing marine conditions may affect body size and fecundity estimates used to evaluate whether stock reference points are met. Using a dataset that spanned five decades, 172,268 individuals, and 19 rivers throughout Eastern Canada, we investigated the occurrence of spatial synchrony in changes in the body size of returning wild adult Atlantic salmon. Body size was then related to conditions in the marine environment (i.e., climate indices, thermal habitat availability, food availability, density-dependence, and fisheries exploitation rates) that may act on all populations during the ocean feeding phase of their life cycle. Body size increased during the 1980s and 1990s for salmon that returned to rivers after one (1SW) or two winters at sea (2SW); however, significant changes were only observed for 1SW and/or 2SW in some mid-latitude and northern rivers (10/13 rivers with 10 of more years of data during these decades) and not in southern rivers (0/2), suggesting weak spatial synchrony across Eastern Canada. For 1SW salmon in nine rivers, body size was longer when fisheries exploitation rates were lower. For 2SW salmon, body size was longer when suitable thermal habitat was more abundant (significant for 3/8 rivers) and the Atlantic Multidecadal Oscillation was higher (i.e., warmer sea surface temperatures; significant for 4/8 rivers). Overall, the weak spatial synchrony and variable effects of covariates on body size across rivers suggest that changes in Atlantic salmon body size may not be solely driven by shared conditions in the marine environment. Regardless, body size changes may have consequences for population management and recovery through the relationship between size and fecundity.

Evaluating the effect of dorsal muscle biopsies on adult Atlantic salmon growth and marine return rates.

Increasing conservation and animal-welfare concerns have driven the development of non-lethal sampling of fish populations, with the use of muscle tissue biopsies now being routinely applied as a sampling method in the wild. Crucial to the success of non-lethal sampling, however, is an evaluation of the short- and long-term consequences of the treatment and ultimately the determination of how these may affect organism mortality and other fitness-related traits. The current study evaluated the use of a dorsal muscle biopsies on post-spawned Atlantic salmon emigrating to sea and undertaking a 2-month long-feeding migration before returning to spawn. Using mark-recapture, return rates and growth were compared between fish that were biopsied and externally tagged, and a control group tagged only with external tags. The biopsy treatment showed no lasting effects on fish as estimated from the two key fitness-related parameters. Results, therefore, suggest the technique can be more widely applied to gather information on marine migrating Atlantic salmon and other anadromous fishes that can be intercepted as they descend and ascend rivers during seasonal migrations. Coupled with modern tagging technologies, the use of biopsies may facilitate an improved understanding of movement and its consequences in terms of feeding patterns and growth.