Sex-dependent dominance at a single locus maintains variation in age at maturity in salmon.

Males and females share many traits that have a common genetic basis; however, selection on these traits often differs between the sexes, leading to sexual conflict. Under such sexual antagonism, theory predicts the evolution of genetic architectures that resolve this sexual conflict. Yet, despite intense theoretical and empirical interest, the specific loci underlying sexually antagonistic phenotypes have rarely been identified, limiting our understanding of how sexual conflict impacts genome evolution and the maintenance of genetic diversity. Here we identify a large effect locus controlling age at maturity in Atlantic salmon (Salmo salar), an important fitness trait in which selection favours earlier maturation in males than females, and show it is a clear example of sex-dependent dominance that reduces intralocus sexual conflict and maintains adaptive variation in wild populations. Using high-density single nucleotide polymorphism data across 57 wild populations and whole genome re-sequencing, we find that the vestigial-like family member 3 gene (VGLL3) exhibits sex-dependent dominance in salmon, promoting earlier and later maturation in males and females, respectively. VGLL3, an adiposity regulator associated with size and age at maturity in humans, explained 39% of phenotypic variation, an unexpectedly large proportion for what is usually considered a highly polygenic trait. Such large effects are predicted under balancing selection from either sexually antagonistic or spatially varying selection. Our results provide the first empirical example of dominance reversal allowing greater optimization of phenotypes within each sex, contributing to the resolution of sexual conflict in a major and widespread evolutionary trade-off between age and size at maturity. They also provide key empirical evidence for how variation in reproductive strategies can be maintained over large geographical scales. We anticipate these findings will have a substantial impact on population management in a range of harvested species where trends towards earlier maturation have been observed.

Life-history genomic regions explain differences in Atlantic salmon marine diet specialization.

Animals employ various foraging strategies along their ontogeny to acquire energy, and with varying degree of efficiencies, to support growth, maturation and subsequent reproduction events. Individuals that can efficiently acquire energy early are more likely to mature at an earlier age, as a result of faster energy gain which can fuel maturation and reproduction. We aimed to test the hypothesis that heritable resource acquisition variation that covaries with efficiency along the ontogeny would influence maturation timing of individuals. To test this hypothesis, we utilized Atlantic salmon as a model which exhibits a simple, hence trackable, genetic control of maturation age. We then monitored the variation in diet acquisition (quantified as stomach fullness and composition) of individuals with different ages, and linked it with genomic regions (haploblocks) that were previously identified to be associated with age-at-maturity. Consistent with the hypothesis, we demonstrated that one of the life-history genomic regions tested (six6) was indeed associated with age-dependent differences in stomach fullness. Prey composition was marginally linked to six6, and suggestively (but non-significantly) to vgll3 genomic regions. We further showed Atlantic salmon switched to the so-called 'feast and famine' strategy along the ontogeny, where older age groups exhibited heavier stomach content, but that came at the expense of running on empty more often. These results suggest genetic variation underlying resource utilization may explain the genetic basis of age structure in Atlantic salmon. Given that ontogenetic diet has a genetic component and the strong spatial diversity associated with these genomic regions, we predict populations with diverse maturation age will have diverse evolutionary responses to future changes in marine food web structures.

Co-inheritance of sea age at maturity and iteroparity in the Atlantic salmon vgll3 genomic region.

Co-inheritance in life-history traits may result in unpredictable evolutionary trajectories if not accounted for in life-history models. Iteroparity (the reproductive strategy of reproducing more than once) in Atlantic salmon (Salmo salar) is a fitness trait with substantial variation within and among populations. In the Teno River in northern Europe, iteroparous individuals constitute an important component of many populations and have experienced a sharp increase in abundance in the last 20 years, partly overlapping with a general decrease in age structure. The physiological basis of iteroparity bears similarities to that of age at first maturity, another life-history trait with substantial fitness effects in salmon. Sea age at maturity in Atlantic salmon is controlled by a major locus around the vgll3 gene, and we used this opportunity demonstrate that these two traits are co-inherited around this genome region. The odds ratio of survival until second reproduction was up to 2.4 (1.8-3.5 90% CI) times higher for fish with the early-maturing vgll3 genotype (EE) compared to fish with the late-maturing genotype (LL). The L allele was dominant in individuals remaining only one year at sea before maturation, but the dominance was reversed, with the E allele being dominant in individuals maturing after two or more years at sea. Post hoc analysis indicated that iteroparous fish with the EE genotype had accelerated growth prior to first reproduction compared to first-time spawners, across all age groups, whereas this effect was not detected in fish with the LL genotype. These results broaden the functional link around the vgll3 genome region and help us understand constraints in the evolution of life-history variation in salmon. Our results further highlight the need to account for genetic correlations between fitness traits when predicting demographic changes in changing environments.

Atlantic salmon populations reveal adaptive divergence of immune related genes - a duplicated genome under selection.

BACKGROUND: Populations of Atlantic salmon display highly significant genetic differences with unresolved molecular basis. These differences may result from separate postglacial colonization patterns, diversifying natural selection and adaptation, or a combination. Adaptation could be influenced or even facilitated by the recent whole genome duplication in the salmonid lineage which resulted in a partly tetraploid species with duplicated genes and regions. RESULTS: In order to elucidate the genes and genomic regions underlying the genetic differences, we conducted a genome wide association study using whole genome resequencing data from eight populations from Northern and Southern Norway. From a total of ~4.5 million sequencing-derived SNPs, more than 10 % showed significant differentiation between populations from these two regions and ten selective sweeps on chromosomes 5, 10, 11, 13-15, 21, 24 and 25 were identified. These comprised 59 genes, of which 15 had one or more differentiated missense mutation. Our analysis showed that most sweeps have paralogous regions in the partially tetraploid genome, each lacking the high number of significant SNPs found in the sweeps. The most significant sweep was found on Chr 25 and carried several missense mutations in the antiviral mx genes, suggesting that these populations have experienced differing viral pressures. Interestingly the second most significant sweep, found on Chr 5, contains two genes involved in the NF-KB pathway (nkap and nkrf), which is also a known pathogen target that controls a large number of processes in animals. CONCLUSION: Our results show that natural selection acting on immune related genes has contributed to genetic divergence between salmon populations in Norway. The differences between populations may have been facilitated by the plasticity of the salmon genome. The observed signatures of selection in duplicated genomic regions suggest that the recently duplicated genome has provided raw material for evolutionary adaptation.

Low but significant genetic differentiation underlies biologically meaningful phenotypic divergence in a large Atlantic salmon population.

Despite decades of research assessing the genetic structure of natural populations, the biological meaning of low yet significant genetic divergence often remains unclear due to a lack of associated phenotypic and ecological information. At the same time, structured populations with low genetic divergence and overlapping boundaries can potentially provide excellent models to study adaptation and reproductive isolation in cases where high-resolution genetic markers and relevant phenotypic and life history information are available. Here, we combined single nucleotide polymorphism (SNP)-based population inference with extensive phenotypic and life history data to identify potential biological mechanisms driving fine-scale subpopulation differentiation in Atlantic salmon (Salmo salar) from the Teno River, a major salmon river in Europe. Two sympatrically occurring subpopulations had low but significant genetic differentiation (FST  = 0.018) and displayed marked differences in the distribution of life history strategies, including variation in juvenile growth rate, age at maturity and size within age classes. Large, late-maturing individuals were virtually absent from one of the two subpopulations, and there were significant differences in juvenile growth rates and size at age after oceanic migration between individuals in the respective subpopulations. Our findings suggest that different evolutionary processes affect each subpopulation and that hybridization and subsequent selection may maintain low genetic differentiation without hindering adaptive divergence.

Genome-wide SNP analysis reveals a genetic basis for sea-age variation in a wild population of Atlantic salmon (Salmo salar).

Delaying sexual maturation can lead to larger body size and higher reproductive success, but carries an increased risk of death before reproducing. Classical life history theory predicts that trade-offs between reproductive success and survival should lead to the evolution of an optimal strategy in a given population. However, variation in mating strategies generally persists, and in general, there remains a poor understanding of genetic and physiological mechanisms underlying this variation. One extreme case of this is in the Atlantic salmon (Salmo salar), which can show variation in the age at which they return from their marine migration to spawn (i.e. their 'sea age'). This results in large size differences between strategies, with direct implications for individual fitness. Here, we used an Illumina Infinium SNP array to identify regions of the genome associated with variation in sea age in a large population of Atlantic salmon in Northern Europe, implementing individual-based genome-wide association studies (GWAS) and population-based FST outlier analyses. We identified several regions of the genome which vary in association with phenotype and/or selection between sea ages, with nearby genes having functions related to muscle development, metabolism, immune response and mate choice. In addition, we found that individuals of different sea ages belong to different, yet sympatric populations in this system, indicating that reproductive isolation may be driven by divergence between stable strategies. Overall, this study demonstrates how genome-wide methodologies can be integrated with samples collected from wild, structured populations to understand their ecology and evolution in a natural context.

Stable isotope analyses on archived fish scales reveal the long-term effect of nitrogen loads on carbon cycling in rivers.

Stable isotope analysis of organic matter in sediment records has long been used to track historical changes in productivity and carbon cycling in marine and lacustrine ecosystems. While flow dynamics preclude stratigraphic measurements of riverine sediments, such retrospective analysis is important for understanding biogeochemical cycling in running waters. Unique collections of riverine fish scales were used to analyse δ(15) N and δ(13) C variations in the food web of two European rivers that experience different degrees of anthropogenic pressure. Over the past four decades, dissolved inorganic N loading remained low and constant in the Teno River (70°N, Finland); in contrast, N loading increased fourfold in the Scorff River (47°N, France) over the same period. Archived scales of Atlantic salmon parr, a riverine life-stage that feeds on aquatic invertebrates, revealed high δ(15) N values in the Scorff River reflecting anthropogenic N inputs to that riverine environment. A strong correlation between dissolved inorganic N loads and δ(13) C values in fish scales was observed in the Scorff River, whereas no trend was found in the Teno River. This result suggests that anthropogenic N-nutrients enhanced atmospheric C uptake by primary producers and its transfer to fish. Our results illustrate for the first time that, as for lakes and marine ecosystems, historical changes in anthropogenic N loading can affect C cycling in riverine food webs, and confirm the long-term interactions between N and C biogeochemical cycles in running waters.

Basin-scale phenology and effects of climate variability on global timing of initial seaward migration of Atlantic salmon (Salmo salar).

Migrations between different habitats are key events in the lives of many organisms. Such movements involve annually recurring travel over long distances usually triggered by seasonal changes in the environment. Often, the migration is associated with travel to or from reproduction areas to regions of growth. Young anadromous Atlantic salmon (Salmo salar) emigrate from freshwater nursery areas during spring and early summer to feed and grow in the North Atlantic Ocean. The transition from the freshwater ('parr') stage to the migratory stage where they descend streams and enter salt water ('smolt') is characterized by morphological, physiological and behavioural changes where the timing of this parr-smolt transition is cued by photoperiod and water temperature. Environmental conditions in the freshwater habitat control the downstream migration and contribute to within- and among-river variation in migratory timing. Moreover, the timing of the freshwater emigration has likely evolved to meet environmental conditions in the ocean as these affect growth and survival of the post-smolts. Using generalized additive mixed-effects modelling, we analysed spatio-temporal variations in the dates of downstream smolt migration in 67 rivers throughout the North Atlantic during the last five decades and found that migrations were earlier in populations in the east than the west. After accounting for this spatial effect, the initiation of the downstream migration among rivers was positively associated with freshwater temperatures, up to about 10 °C and levelling off at higher values, and with sea-surface temperatures. Earlier migration occurred when river discharge levels were low but increasing. On average, the initiation of the smolt seaward migration has occurred 2.5 days earlier per decade throughout the basin of the North Atlantic. This shift in phenology matches changes in air, river, and ocean temperatures, suggesting that Atlantic salmon emigration is responding to the current global climate changes.

Cost-effective genome-wide estimation of allele frequencies from pooled DNA in Atlantic salmon (Salmo salar L.).

BACKGROUND: New sequencing technologies have tremendously increased the number of known molecular markers (single nucleotide polymorphisms; SNPs) in a variety of species. Concurrently, improvements to genotyping technology have now made it possible to efficiently genotype large numbers of genome-wide distributed SNPs enabling genome wide association studies (GWAS). However, genotyping significant numbers of individuals with large number of SNPs remains prohibitively expensive for many research groups. A possible solution to this problem is to determine allele frequencies from pooled DNA samples, such 'allelotyping' has been presented as a cost-effective alternative to individual genotyping and has become popular in human GWAS. In this article we have tested the effectiveness of DNA pooling to obtain accurate allele frequency estimates for Atlantic salmon (Salmo salar L.) populations using an Illumina SNP-chip. RESULTS: In total, 56 Atlantic salmon DNA pools from 14 populations were analyzed on an Atlantic salmon SNP-chip containing probes for 5568 SNP markers, 3928 of which were bi-allelic. We developed an efficient quality control filter which enables exclusion of loci showing high error rate and minor allele frequency (MAF) close to zero. After applying multiple quality control filters we obtained allele frequency estimates for 3631 bi-allelic loci. We observed high concordance (r > 0.99) between allele frequency estimates derived from individual genotyping and DNA pools. Our results also indicate that even relatively small DNA pools (35 individuals) can provide accurate allele frequency estimates for a given sample. CONCLUSIONS: Despite of higher level of variation associated with array replicates compared to pool construction, we suggest that both sources of variation should be taken into account. This study demonstrates that DNA pooling allows fast and high-throughput determination of allele frequencies in Atlantic salmon enabling cost-efficient identification of informative markers for discrimination of populations at various geographical scales, as well as identification of loci controlling ecologically and economically important traits.

Fish scales and SNP chips: SNP genotyping and allele frequency estimation in individual and pooled DNA from historical samples of Atlantic salmon (Salmo salar).

BACKGROUND: DNA extracted from historical samples is an important resource for understanding genetic consequences of anthropogenic influences and long-term environmental change. However, such samples generally yield DNA of a lower amount and quality, and the extent to which DNA degradation affects SNP genotyping success and allele frequency estimation is not well understood. We conducted high density SNP genotyping and allele frequency estimation in both individual DNA samples and pooled DNA samples extracted from dried Atlantic salmon (Salmo salar) scales stored at room temperature for up to 35 years, and assessed genotyping success, repeatability and accuracy of allele frequency estimation using a high density SNP genotyping array. RESULTS: In individual DNA samples, genotyping success and repeatability was very high (> 0.973 and > 0.998, respectively) in samples stored for up to 35 years; both increased with the proportion of DNA of fragment size > 1000 bp. In pooled DNA samples, allele frequency estimation was highly repeatable (Repeatability = 0.986) and highly correlated with empirical allele frequency measures (Mean Adjusted R2 = 0.991); allele frequency could be accurately estimated in > 95% of pooled DNA samples with a reference group of at least 30 individuals. SNPs located in polyploid regions of the genome were more sensitive to DNA degradation: older samples had lower genotyping success at these loci, and a larger reference panel of individuals was required to accurately estimate allele frequencies. CONCLUSIONS: SNP genotyping was highly successful in degraded DNA samples, paving the way for the use of degraded samples in SNP genotyping projects. DNA pooling provides the potential for large scale population genetic studies with fewer assays, provided enough reference individuals are also genotyped and DNA quality is properly assessed beforehand. We provide recommendations for future studies intending to conduct high-throughput SNP genotyping and allele frequency estimation in historical samples.

Gene flow from domesticated escapes alters the life history of wild Atlantic salmon.

Interbreeding between domesticated and wild animals occurs in several species. This gene flow has long been anticipated to induce genetic changes in life-history traits of wild populations, thereby influencing population dynamics and viability. Here, we show that individuals with high levels of introgression (domesticated ancestry) have altered age and size at maturation in 62 wild Atlantic salmon Salmo salar populations, including seven ancestral populations to breeding lines of the domesticated salmon. This study documents widespread changes to life-history traits in wild animal populations following gene flow from selectively bred, domesticated conspecifics. The continued high abundance of escaped, domesticated Atlantic salmon thus threatens wild Atlantic salmon populations by inducing genetic changes in fitness-related traits. Our results represent key evidence and a timely warning concerning the potential ecological impacts of the globally increasing use of domesticated animals.

Single nucleotide polymorphisms to discriminate different classes of hybrid between wild Atlantic salmon and aquaculture escapees.

Many wild Atlantic salmon (Salmo salar) populations are threatened by introgressive hybridization from domesticated fish that have escaped from aquaculture facilities. A detailed understanding of the hybridization dynamics between wild salmon and aquaculture escapees requires discrimination of different hybrid classes; however, markers currently available to discriminate the two types of parental genome have limited power to do this. Using a high-density Atlantic salmon single nucleotide polymorphism (SNP) array, in combination with pooled-sample allelotyping and an Fst outlier approach, we identified 200 SNPs that differentiated an important Atlantic salmon stock from the escapees potentially hybridizing with it. By simulating multiple generations of wild-escapee hybridization, involving wild populations in two major phylogeographic lineages and a genetically diverse set of escapees, we showed that both the complete set of SNPs and smaller subsets could reliably assign individuals to different hybrid classes up to the third hybrid (F3) generation. This set of markers will be a useful tool for investigating the genetic interactions between native wild fish and aquaculture escapees in many Atlantic salmon populations.

Finding markers that make a difference: DNA pooling and SNP-arrays identify population informative markers for genetic stock identification.

Genetic stock identification (GSI) using molecular markers is an important tool for management of migratory species. Here, we tested a cost-effective alternative to individual genotyping, known as allelotyping, for identification of highly informative SNPs for accurate genetic stock identification. We estimated allele frequencies of 2880 SNPs from DNA pools of 23 Atlantic salmon populations using Illumina SNP-chip. We evaluated the performance of four common strategies (global F ST, pairwise F ST, Delta and outlier approach) for selection of the most informative set of SNPs and tested their effectiveness for GSI compared to random sets of SNP and microsatellite markers. For the majority of cases, SNPs selected using the outlier approach performed best followed by pairwise F ST and Delta methods. Overall, the selection procedure reduced the number of SNPs required for accurate GSI by up to 53% compared with randomly chosen SNPs. However, GSI accuracy was more affected by populations in the ascertainment group rather than the ranking method itself. We demonstrated for the first time the compatibility of different large-scale SNP datasets by compiling the largest population genetic dataset for Atlantic salmon to date. Finally, we showed an excellent performance of our top SNPs on an independent set of populations covering the main European distribution range of Atlantic salmon. Taken together, we demonstrate how combination of DNA pooling and SNP arrays can be applied for conservation and management of salmonids as well as other species.

Temporally stable population-specific differences in run timing of one-sea-winter Atlantic salmon returning to a large river system.

The understanding of migration patterns can significantly contribute to conservation and management. The spawning migrations of Atlantic salmon (Salmo salar) cover thousands of kilometers from the feeding areas at sea to their natal rivers to reproduce. Migrating salmon are exposed to intensive harvest, but little is known of the population-specific differences in migration behavior. In this study, timing of return migration was investigated among one-sea-winter Atlantic salmon within a river system. By utilizing knowledge of the genetic population structure, population of origin was reliably identified for c. 1500 fish caught in mixed stock fisheries after adopting an approach to minimize the complications arising from potential nonsampled populations. Results demonstrated significant and temporally stable differences among populations as well as between sexes. Generally, female salmon from tributary populations entered fresh water first. Run timing was not however related to in-river migration distance. Rather, one-sea-winter salmon from larger populations and with a higher proportion of multi-sea-winter females arrived later in the season. These findings are a significant step toward a more thorough understanding of the salmon migration behavior and behavioral ecology, providing concrete tools for the management and conservation of the remaining indigenous Atlantic salmon stocks.

Temporally stable genetic structure and low migration in an Atlantic salmon population complex: implications for conservation and management.

The evolutionary potential of a population is closely related to two key population genetic parameters, namely the effective population size (N e) and migration rate (m). Furthermore, knowledge of these parameters is required in order to assess potential constraints on local adaptation and for the development of biologically sound management strategies. We addressed these key issues by investigating the temporal and spatial genetic structure of over 2000 adult Atlantic salmon (Salmo salar) collected from 17 sites in the Teno and Näätämö rivers in northernmost Europe with up to five time points spanning temporal intervals up to 24 years (∼4 generations). In all cases except one, local populations were found to be temporally stable within the river system. Estimates of N e were generally a magnitude larger for the mainstem and headwater populations (MS+HW, N e∼340-1200) than for the tributary populations (N e∼35-160), thus explaining the higher genetic diversity and lower divergence of the MS+HW populations compared to tributaries. The overall migration rates to tributaries were low, and in some cases, low enough for local adaptations to potentially evolve, despite their lower N e. Signs of a population bottleneck and natural recruitment from nearby populations were detected in one local population. This highlights a fact which is relevant for the conservation and management of highly substructured population systems in general: that even when the overall census size is large, local populations can be vulnerable to perturbations. To preserve the current and to regain the historical distribution of salmon within the river system, we propose that the status of the total population complex should be evaluated at the local population level rather than from descriptive statistics at the system level.

Life-history and habitat features influence the within-river genetic structure of Atlantic salmon.

Defining populations and identifying ecological and life-history characteristics affecting genetic structure is important for understanding species biology and hence, for managing threatened or endangered species or populations. In this study, populations of the world's largest indigenous Atlantic salmon (Salmo salar) stock were first inferred using model-based clustering methods, following which life-history and habitat variables best predicting the genetic diversity of populations were identified. This study revealed that natal homing of Atlantic salmon within the Teno River system is accurate at least to the tributary level. Generally, defining populations by main tributaries was observed to be a reasonable approach in this large river system, whereas in the mainstem of the river, the number of inferred populations was fewer than the number of distinct sampling sites. Mainstem and headwater populations were genetically more diverse and less diverged, while each tributary fostered a distinct population with high genetic differentiation and lower genetic diversity. Population structure and variation in genetic diversity among populations were poorly explained by geographical distance. In contrast, age-structure, as estimated by the proportion of multisea-winter spawners, was the most predictive variable in explaining the variation in the genetic diversity of the populations. This observation, being in agreement with theoretical predictions, emphasizes the essence of large multisea-winter females in maintaining the genetic diversity of populations. In addition, the unique genetic diversity of populations, as estimated by private allele richness, was affected by the ease of accessibility of a site, with more difficult to access sites having lower unique genetic diversity. Our results show that despite this species' high capacity for migration, tributaries foster relatively closed populations with little gene flow which will be important to consider when developing management strategies for the system.