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.

Home ground advantage: Local Atlantic salmon have higher reproductive fitness than dispersers in the wild.

A long-held, but poorly tested, assumption in natural populations is that individuals that disperse into new areas for reproduction are at a disadvantage compared to individuals that reproduce in their natal habitat, underpinning the eco-evolutionary processes of local adaptation and ecological speciation. Here, we capitalize on fine-scale population structure and natural dispersal events to compare the reproductive success of local and dispersing individuals captured on the same spawning ground in four consecutive parent-offspring cohorts of wild Atlantic salmon (Salmo salar). Parentage analysis conducted on adults and juvenile fish showed that local females and males had 9.6 and 2.9 times higher reproductive success than dispersers, respectively. Our results reveal how higher reproductive success in local spawners compared to dispersers may act in natural populations to drive population divergence and promote local adaptation over microgeographic spatial scales without clear morphological differences between populations.

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.

Genomic signatures of fine-scale local selection in Atlantic salmon suggest involvement of sexual maturation, energy homeostasis and immune defence-related genes.

Elucidating the genetic basis of adaptation to the local environment can improve our understanding of how the diversity of life has evolved. In this study, we used a dense SNP array to identify candidate loci potentially underlying fine-scale local adaptation within a large Atlantic salmon (Salmo salar) population. By combining outlier, gene-environment association and haplotype homozygosity analyses, we identified multiple regions of the genome with strong evidence for diversifying selection. Several of these candidate regions had previously been identified in other studies, demonstrating that the same loci could be adaptively important in Atlantic salmon at subdrainage, regional and continental scales. Notably, we identified signals consistent with local selection around genes associated with variation in sexual maturation, energy homeostasis and immune defence. These included the large-effect age-at-maturity gene vgll3, the known obesity gene mc4r, and major histocompatibility complex II. Most strikingly, we confirmed a genomic region on Ssa09 that was extremely differentiated among subpopulations and that is also a candidate for local selection over the global range of Atlantic salmon. This region colocalized with a haplotype strongly associated with spawning ecotype in sockeye salmon (Oncorhynchus nerka), with circumstantial evidence that the same gene (six6) may be the selective target in both cases. The phenotypic effect of this region in Atlantic salmon remains cryptic, although allelic variation is related to upstream catchment area and covaries with timing of the return spawning migration. Our results further inform management of Atlantic salmon and open multiple avenues for future research.

First record of proliferative kidney disease agent Tetracapsuloides bryosalmonae in wild brown trout and European grayling in Finland.

The myxozoan endoparasite Tetracapsuloides bryosalmonae causes temperature-driven proliferative kidney disease (PKD) in salmonid fishes. Despite the economic and ecological importance of PKD, information about the distribution of the parasite is still scarce. Here, we report for the first time the occurrence of T. bryosalmonae in wild brown trout Salmo trutta and European grayling Thymallus thymallus populations in Finland. We detected T. bryosalmonae at high prevalence in both brown trout and European grayling from the transboundary Finnish-Russian River Koutajoki system (Rivers Oulankajoki, Kuusinkijoki, Kitkajoki, Maaninkajoki, and Juumajoki) in north-eastern Finland. In southern Finland, T. bryosalmonae was detected in River Siuntionjoki young-of-the-year brown trout collected both in 2015 and 2016 (100% prevalence), while the parasite was not observed in fish from 3 other rivers (Ingarskila, Mustajoki, and Vantaanjoki) flowing to the Gulf of Finland. Our results, together with those from recent studies of Atlantic salmon, indicate that T. bryosalmonae is distributed over much higher latitudes in northern Europe than previously appreciated. We expect that increasing water temperatures will likely cause new PKD outbreaks in these more northerly regions in the future.

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.

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.

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.

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.

A proteomics approach reveals divergent molecular responses to salinity in populations of European whitefish (Coregonus lavaretus).

Osmoregulation is a vital physiological function for fish, as it helps maintain a stable intracellular concentration of ions in environments of variable salinities. We focused on a primarily freshwater species, the European whitefish (Coregonus lavaretus), to investigate the molecular mechanisms underlying salinity tolerance and examine whether these mechanisms differ between genetically similar populations that spawn in freshwater vs. brackishwater environments. A common garden experiment involving 27 families in two populations and five salinity treatments together with a large-scale, high-resolution mass spectrometry experiment that quantified 1500 proteins was conducted to assess phenotypic and proteomic responses during early development, from fertilization until hatching, in the studied populations. The populations displayed drastically different phenotypic and proteomic responses to salinity. Freshwater-spawning whitefish showed a significantly higher mortality rate in higher salinity treatments. Calcium, an ion involved in osmotic stress sensing, had a central role in the observed proteomic responses. Brackishwater-spawning fish were capable of viable osmoregulation, which was modulated by cortisol, an important seawater-adaptation hormone in teleost fish. Several proteins were identified to play key roles in osmoregulation, most importantly a highly conserved cytokine, tumour necrosis factor, whereas calcium receptor activities were associated with salinity adaptation. These results imply that individuals from these populations are most likely adapted to their local environments, even though the baseline level of genetic divergence between them is low (F(ST)=0.049). They also provide clues for choosing candidate loci for studying the molecular basis of salinity adaptation in other species. Further, our approach provides an example of how proteomic methods can be successfully used to obtain novel insights into the molecular mechanisms behind adaptation in non-model organism.

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.

Efficiency of model-based Bayesian methods for detecting hybrid individuals under different hybridization scenarios and with different numbers of loci.

Accurate detection of offspring resulting from hybridization between individuals of distinct populations has a range of applications in conservation and population genetics. We assessed the hybrid identification efficiency of two methods (implemented in the STRUCTURE and NEWHYBRIDS programs) which are tailored to identifying hybrid individuals but use different approaches. Simulated first- and second-generation hybrids were used to assess the performance of these two methods in detecting recent hybridization under scenarios with different levels of genetic divergence and varying numbers of loci. Despite the different approaches of the methods, the hybrid detection efficiency was generally similar and neither of the two methods outperformed the other in all scenarios assessed. Interestingly, hybrid detection efficiency was only minimally affected by whether reference population allele frequency information was included or not. In terms of genotyping effort, efficient detection of F1 hybrid individuals requires the use of 12 or 24 loci with pairwise F(ST) between hybridizing parental populations of 0.21 or 0.12, respectively. While achievable, these locus numbers are nevertheless higher than the number of loci currently commonly applied in population genetic studies. The method of STRUCTURE seemed to be less sensitive to the proportion of hybrids included in the sample, while NEWHYBRIDS seemed to perform slightly better when individuals from both backcross and F1 hybrid classes were present in the sample. However, separating backcrosses from purebred parental individuals requires a considerable genotyping effort (at least 48 loci), even when divergence between parental populations is high.