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.

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.

A semi-automated spectral approach to analyzing cyclical growth patterns using fish scales.

We introduce a new semi-automated approach to analyzing growth patterns recorded on fish scales. After manually specifying the center of the scale, the algorithm radially unwraps the scale patterns along a series of transects from the center to the edge of the scale. A sliding window Fourier transform is used to produce a spectrogram for each sampled transect of the scale image. The maximum frequency over all sampled transects of the average spectrogram yields a well-discriminated peak frequency trace that can then serve as a growth template for that fish. The spectrogram patterns of individual fish scales can be adjusted to a common period accounting for differences in date of return or size of fish at return without biasing the growth profile of the scale. We apply the method to 147 Atlantic salmon scale images sampled from 3 years and contrast the information derived with this automated approach to what is obtained using classical human operator measurements. The spectrogram analysis quantifies growth patterns using the entire scale image rather than just a single transect and provides the possibility of more robustly analyzing individual scale growth patterns. This semi-automated approach that removes essentially all the human operator interventions provides an opportunity to process large datasets of fish scale images and combined with advanced analyses such as deep learning methods could lead to a greater understanding of salmon marine migration patterns and responses to variations in ecosystem conditions.

Estimating multiple years, tributary-specific, and overall Atlantic salmon smolt abundance in a large Canadian catchment using capture-mark-recapture experiments.

Population monitoring of Atlantic salmon (Salmo salar L.) abundance is an essential element to understand annual stock variability and inform fisheries management processes. Smolts are the life stage marking the transition from the freshwater to the marine phase of anadromous Atlantic salmon. Estimating smolt abundance allows for subsequent inferences on freshwater and marine survival rates. Annual abundances of out-migrating Atlantic salmon smolts were estimated using Bayesian models and an 18-year capture-mark-recapture time series from two to five trapping locations within the Restigouche River (Canada) catchment. Some of the trapping locations were at the outlet of large upstream tributaries, and these sampled a portion of the total out-migrating population of smolts for the watershed, whereas others were located just above the head of tide of the Restigouche River and sampled the entire run of salmon smolts. Due to logistic and environmental conditions, not all trapping locations were operational each year. Additionally, recapture rates were relatively low (<5%), and the absolute number of recaptures was relatively few (most often a few dozen), leading to incoherent and highly uncertain estimates of tributary-specific and whole catchment abundance estimates when the data were modeled independently among trapping locations and years. Several models of increasing complexity were tested using simulated data, and the best-performing model in terms of bias and precision incorporated a hierarchical structure among years on the catchability parameters and included an explicit spatial structure to account for the annual variations in the number of sampled locations within the watershed. When the best model was applied to the Restigouche River catchment dataset, the annual smolt abundance estimates varied from 250,000 to 1 million smolts, and the subbasin estimates of abundance were consistent with the spatial structure of the monitoring programme. Ultimately, increasing the probabilities of capture and the absolute number of recaptures at the different traps will be required to improve the precision and reduce the bias of the estimates of smolt abundance for the entire basin and within subbasins of the watershed. The model and approach provide a significant improvement in the models used to date based on independent estimates of abundance by trapping location and year. Total abundance and relative production in discrete spawning, nesting, or rearing areas provide critical information to appropriately understand and manage the threats to species that can occur at subpopulation spatial scales.

Spatial synchrony in the response of a long range migratory species (Salmo salar) to climate change in the North Atlantic Ocean.

A major challenge in understanding the response of populations to climate change is to separate the effects of local drivers acting independently on specific populations, from the effects of global drivers that impact multiple populations simultaneously and thereby synchronize their dynamics. We investigated the environmental drivers and the demographic mechanisms of the widespread decline in marine survival rates of Atlantic salmon (Salmo salar) over the last four decades. We developed a hierarchical Bayesian life cycle model to quantify the spatial synchrony in the marine survival of 13 large groups of populations (called stock units, SU) from two continental stock groups (CSG) in North America (NA) and Southern Europe (SE) over the period 1971-2014. We found strong coherence in the temporal variation in postsmolt marine survival among the 13 SU of NA and SE. A common North Atlantic trend explains 37% of the temporal variability of the survivals for the 13 SU and declines by a factor of 1.8 over the 1971-2014 time series. Synchrony in survival trends is stronger between SU within each CSG. The common trends at the scale of NA and SE capture 60% and 42% of the total variance of temporal variations, respectively. Temporal variations of the postsmolt survival are best explained by the temporal variations of sea surface temperature (SST, negative correlation) and net primary production indices (PP, positive correlation) encountered by salmon in common domains during their marine migration. Specifically, in the Labrador Sea/Grand Banks for populations from NA, 26% and 24% of variance is captured by SST and PP, respectively and in the Norwegian Sea for populations from SE, 21% and 12% of variance is captured by SST and PP, respectively. The findings support the hypothesis of a response of salmon populations to large climate-induced changes in the North Atlantic simultaneously impacting populations from distant continental habitats.

Spawning-strategy-dependent diets in two North American populations of Atlantic salmon Salmo salar.

The diet of repeat-spawner Atlantic salmon Salmo salar was investigated using carbon and nitrogen stable-isotope values from the outer growth band of scales, which reflect the fish's consumption and growth during their most recent marine phase. Isotope values for S. salar displaying different spawning strategies were compared between and within the Miramichi and Nashwaak Rivers, New Brunswick, Canada and a Bayesian mixing model was used to infer dietary contributions from potential prey items. Significant differences in the stable-isotope values were found among spawning strategies and between rivers, indicating differences in diet and feeding area, consistent with hypotheses. Bayesian mixing model results inferred the main prey items consumed during marine feeding by S. salar to consist of hyperiid amphipods and capelin Mallotus villosus for repeat alternate spawners from both rivers, sandlance Ammodytes sp. for repeat consecutive spawners from the Miramichi River and amphipods for repeat consecutive spawners from the Nashwaak River. These results demonstrate the diversity of feeding tactics among S. salar spawning strategies from the same river and between populations from different rivers. Accounting for differences in prey availability and the subsequent impact on S. salar diet and spawner return rates (i.e., marine survival) will facilitate the application of ecosystem-based management practices, such as ensuring that fisheries for forage species do not indirectly adversely affect S. salar return rates.

Conservation genomics of anadromous Atlantic salmon across its North American range: outlier loci identify the same patterns of population structure as neutral loci.

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.

Development of a short-term in situ caging methodology to assess long-term effects of industrial and municipal discharges on salmon smolts.

Trial experiments to develop an in situ method for determining effects of short-term exposure to contaminants on salmon during the sensitive smolt stage were carried out. Wild Atlantic salmon (Salmo salar) smolts were caged at various estuarine locations in northeast New Brunswick, Canada at different levels of anthropogenic impact in 2000 and 2001. Survival, growth (weight and K), and feeding (delta(13)C and delta(15)N) parameters were measured following summer grow-out at a hatchery. A summary of distributions of the four parameters indicated that smolts caged at locations receiving lower levels of industrial and municipal discharges feed and grow better than those caged at sites receiving higher levels. The findings, however, were not repeatable between the 2 years studied. The observed inconsistency between years may be a consequence of the relatively low concentrations of alkylphenolic contaminants (putative causative agents) and overall steroidogenic potency in river water at the time of caging. Differences in temperature and salinity from 2000 to 2001 may have further confounded comparisons across treatments and between years. In future studies, caging in closer proximity to industrial and municipal discharges or in systems with higher concentrations of waterborne contaminants or impoundments would help further the assessment and applicability of this methodology and allow a more robust comparison of smolt feeding and growth among reference and exposure sites.