Patterns of genomic variation in Coho salmon following reintroduction to the interior Columbia River

Coho salmon were extirpated in the mid‐20th century from the interior reaches of the Columbia River but were reintroduced with relatively abundant source stocks from the lower Columbia River near the Pacific coast. Reintroduction of Coho salmon to the interior Columbia River (Wenatchee River) using lower river stocks placed selective pressures on the new colonizers due to substantial differences with their original habitat such as migration distance and navigation of six additional hydropower dams. We used restriction site‐associated DNA sequencing (RAD‐seq) to genotype 5,392 SNPs in reintroduced Coho salmon in the Wenatchee River over four generations to test for signals of temporal structure and adaptive variation. Temporal genetic structure among the three broodlines of reintroduced fish was evident among the initial return years (2000, 2001, and 2002) and their descendants, which indicated levels of reproductive isolation among broodlines. Signals of adaptive variation were detected from multiple outlier tests and identified candidate genes for further study. This study illustrated that genetic variation and structure of reintroduced populations are likely to reflect source stocks for multiple generations but may shift over time once established in nature.     

Adaptive divergence in embryonic thermal plasticity among Atlantic salmon populations

In the context of global changes, the long‐term viability of populations of endangered ectotherms may depend on their adaptive potential and ability to cope with temperature variations. We measured responses of Atlantic salmon embryos from four populations to temperature variations and used a Q_ST–F_ST approach to study the adaptive divergence among these populations. Embryos were reared under two experimental conditions: a low temperature regime at 4 °C until eyed‐stage and 10 °C until the end of embryonic development and a high temperature regime with a constant temperature of 10 °C throughout embryonic development. Significant variations among populations and population × temperature interactions were observed for embryo survival, incubation time and length. Q_ST was higher than F_ST in all but one comparison suggesting an important effect of divergent selection. Q_ST was also higher under the high‐temperature treatment than at low temperature for length and survival due to a higher variance among populations under the stressful warmer treatment. Interestingly, heritability was lower for survival under high temperature in relation to a lower additive genetic variance under that treatment. Overall, these results reveal an adaptive divergence in thermal plasticity in embryonic life stages of Atlantic salmon suggesting that salmon populations may differentially respond to temperature variations induced by climate change. These results also suggest that changes in temperature may alter not only the adaptive potential of natural populations but also the selection regimes among them.     

Inter- and intrapopulation variation in temperature sum requirements at hatching in Norwegian Atlantic salmon

The duration of the ontogenetic development of the eggs up to hatching time was studied under stable water temperatures for a total of 47 families (the offspring of a single pair of spawners) from eight wild populations and one domesticated stock of Norwegian Atlantic salmon. Eggs from each family were incubated at stable water temperatures of 7·3 and 9·4° C. The number of days required to 50% hatching differed significantly between populations for both experimental temperatures together and for each temperature separately. The interpopulation and intrapopulation variation each accounted for about 50% of the overall variation. Significant differences were found for the days required to 50% hatching between eggs of the domestic stock from The Norwegian Salmon Breeding Company and two natural populations and between some of the wild populations. It is concluded that the time required from fertilisation to hatching in Atlantic salmon is a population-specific adaptive trait.     

A critical review of adaptive genetic variation in Atlantic salmon: implications for conservation

Here we critically review the scale and extent of adaptive genetic variation in Atlantic salmon (Salmo salar L.), an important model system in evolutionary and conservation biology that provides fundamental insights into population persistence, adaptive response and the effects of anthropogenic change. We consider the process of adaptation as the end product of natural selection, one that can best be viewed as the degree of matching between phenotype and environment. We recognise three potential sources of adaptive variation: heritable variation in phenotypic traits related to fitness, variation at the molecular level in genes influenced by selection, and variation in the way genes interact with the environment to produce phenotypes of varying plasticity. Of all phenotypic traits examined, variation in body size (or in correlated characters such as growth rates, age of seaward migration or age at sexual maturity) generally shows the highest heritability, as well as a strong effect on fitness. Thus, body size in Atlantic salmon tends to be positively correlated with freshwater and marine survival, as well as with fecundity, egg size, reproductive success, and offspring survival. By contrast, the fitness implications of variation in behavioural traits such as aggression, sheltering behaviour, or timing of migration are largely unknown. The adaptive significance of molecular variation in salmonids is also scant and largely circumstantial, despite extensive molecular screening on these species. Adaptive variation can result in local adaptations (LA) when, among other necessary conditions, populations live in patchy environments, exchange few or no migrants, and are subjected to differential selective pressures. Evidence for LA in Atlantic salmon is indirect and comes mostly from ecological correlates in fitness-related traits, the failure of many translocations, the poor performance of domesticated stocks, results of a few common-garden experiments (where different populations were raised in a common environment in an attempt to dissociate heritable from environmentally induced phenotypic variation), and the pattern of inherited resistance to some parasites and diseases. Genotype x environment interactions occurr for many fitness traits, suggesting that LA might be important. However, the scale and extent of adaptive variation remains poorly understood and probably varies, depending on habitat heterogeneity, environmental stability and the relative roles of selection and drift. As maladaptation often results from phenotype-environment mismatch, we argue that acting as if populations are not locally adapted carries a much greater risk of mismanagement than acting under the assumption for local adaptations when there are none. As such, an evolutionary approach to salmon conservation is required, aimed at maintaining the conditions necessary for natural selection to operate most efficiently and unhindered. This may require minimising alterations to native genotypes and habitats to which populations have likely become adapted, but also allowing for population size to reach or extend beyond carrying capacity to encourage competition and other sources of natural mortality.     

Identification of quantitative genetic components of fitness variation in farmed,hybrid and native salmon in the wild

Feral animals represent an important problem in many ecosystems due to interbreeding with wild conspecifics. Hybrid offspring from wild and domestic parents are often less adapted to local environment and ultimately, can reduce the fitness of the native population. This problem is an important concern in Norway, where each year, hundreds of thousands of farm Atlantic salmon escape from fish farms. Feral fish outnumber wild populations, leading to a possible loss of local adaptive genetic variation and erosion of genetic structure in wild populations. Studying the genetic factors underlying relative performance between wild and domesticated conspecific can help to better understand how domestication modifies the genetic background of populations, and how it may alter their ability to adapt to the natural environment. Here, based upon a large-scale release of wild, farm and wild x farm salmon crosses into a natural river system, a genome-wide quantitative trait locus (QTL) scan was performed on the offspring of 50 full-sib families, for traits related to fitness (length, weight, condition factor and survival). Six QTLs were detected as significant contributors to the phenotypic variation of the first three traits, explaining collectively between 9.8 and 14.8% of the phenotypic variation. The seventh QTL had a significant contribution to the variation in survival, and is regarded as a key factor to understand the fitness variability observed among salmon in the river. Interestingly, strong allelic correlation within one of the QTL regions in farmed salmon might reflect a recent selective sweep due to artificial selection.     

Family‐assisted inference of the genetic architecture of major histocompatibility complex variation

With their direct link to individual fitness, genes of the major histocompatibility complex (MHC) are a popular system to study the evolution of adaptive genetic diversity. However, owing to the highly dynamic evolution of the MHC region, the isolation, characterization and genotyping of MHC genes remain a major challenge. While high‐throughput sequencing technologies now provide unprecedented resolution of the high allelic diversity observed at the MHC, in many species, it remains unclear (i) how alleles are distributed among MHC loci, (ii) whether MHC loci are linked or segregate independently and (iii) how much copy number variation (CNV) can be observed for MHC genes in natural populations. Here, we show that the study of allele segregation patterns within families can provide significant insights in this context. We sequenced two MHC class I (MHC‐I) loci in 1267 European barn owls (Tyto alba), including 590 offspring from 130 families using Illumina MiSeq technology. Coupled with a high per‐individual sequencing coverage (~3000×), the study of allele segregation patterns within families provided information on three aspects of the architecture of MHC‐I variation in barn owls: (i) extensive sharing of alleles among loci, (ii) strong linkage of MHC‐I loci indicating tandem architecture and (iii) the presence of CNV in the barn owl MHC‐I. We conclude that the additional information that can be gained from high‐coverage amplicon sequencing by investigating allele segregation patterns in families not only helps improving the accuracy of MHC genotyping, but also contributes towards enhanced analyses in the context of MHC evolutionary ecology.     

A mobile sex‐determining region,male‐specific haplotypes and rearing environment influence age at maturity in Chinook salmon

Variation in age at maturity is an important contributor to life history and demographic variation within and among species. The optimal age at maturity can vary by sex, and the ability of each sex to evolve towards its fitness optimum depends on the genetic architecture of maturation. Using GWAS of RAD sequencing data, we show that age at maturity in Chinook salmon exhibits sex‐specific genetic architecture, with age at maturity in males influenced by large (up to 20 Mb) male‐specific haplotypes. These regions showed no such effect in females. We also provide evidence for translocation of the sex‐determining gene between two different chromosomes. This has important implications for sexually antagonistic selection, particularly that sex linkage of adaptive genes may differ within and among populations based on chromosomal location of the sex‐determining gene. Our findings will facilitate research into the genetic causes of shifting demography in Chinook salmon as well as a better understanding of sex determination in this species and Pacific salmon in general.     

Development of a genotype‐by‐sequencing immunogenetic assay as exemplified by screening for variation in red fox with and without endemic rabies exposure

Pathogens are recognized as major drivers of local adaptation in wildlife systems. By determining which gene variants are favored in local interactions among populations with and without disease, spatially explicit adaptive responses to pathogens can be elucidated. Much of our current understanding of host responses to disease comes from a small number of genes associated with an immune response. High‐throughput sequencing (HTS) technologies, such as genotype‐by‐sequencing (GBS), facilitate expanded explorations of genomic variation among populations. Hybridization‐based GBS techniques can be leveraged in systems not well characterized for specific variants associated with disease outcome to “capture” specific genes and regulatory regions known to influence expression and disease outcome. We developed a multiplexed, sequence capture assay for red foxes to simultaneously assess ~300‐kbp of genomic sequence from 116 adaptive, intrinsic, and innate immunity genes of predicted adaptive significance and their putative upstream regulatory regions along with 23 neutral microsatellite regions to control for demographic effects. The assay was applied to 45 fox DNA samples from Alaska, where three arctic rabies strains are geographically restricted and endemic to coastal tundra regions, yet absent from the boreal interior. The assay provided 61.5% on‐target enrichment with relatively even sequence coverage across all targeted loci and samples (mean = 50×), which allowed us to elucidate genetic variation across introns, exons, and potential regulatory regions (4,819 SNPs). Challenges remained in accurately describing microsatellite variation using this technique; however, longer‐read HTS technologies should overcome these issues. We used these data to conduct preliminary analyses and detected genetic structure in a subset of red fox immune‐related genes between regions with and without endemic arctic rabies. This assay provides a template to assess immunogenetic variation in wildlife disease systems.     

Allozyme variation in Atlantic salmon from the British Isles: associations with geography and the environment

Data on geographical variation in allele frequencies at enzyme coding loci in Atlantic salmon from the British Isles were collated from published and unpublished sources. Statistically significant differences in allele frequencies were found among samples both within and among river systems, suggesting that the Atlantic salmon in the British Isles is not a panmictic population and that even within major river systems it cannot be treated as a single genetic stock for fisheries management purposes. Although there was some evidence of regional differences in the frequency of some rare alleles, most single‐locus variation did not show strong geographic patterns, with the exception of the AAT‐4 * locus at which allele frequencies had a significant latitudinal cline. There was some evidence for the existence of genetically‐distinct celtic and boreal races of Atlantic salmon in the British Isles as previously has been suggested. Multiple regression analyses revealed associations between genetic variation and local environmental conditions ( i.e . between variation at MEP‐2 * and both temperature and local river gradient), providing additional evidence for adaptive population divergence in the species.     

Thermal adaptation and phenotypic plasticity in a warming world: Insights from common garden experiments on Alaskan sockeye salmon

An important unresolved question is how populations of coldwater‐dependent fishes will respond to rapidly warming water temperatures. For example, the culturally and economically important group, Pacific salmon (Oncorhynchus spp.), experience site‐specific thermal regimes during early development that could be disrupted by warming. To test for thermal local adaptation and heritable phenotypic plasticity in Pacific salmon embryos, we measured the developmental rate, survival, and body size at hatching in two populations of sockeye salmon (Oncorhynchus nerka) that overlap in timing of spawning but incubate in contrasting natural thermal regimes. Using a split half‐sibling design, we exposed embryos of 10 families from each of two populations to variable and constant thermal regimes. These represented both experienced temperatures by each population, and predicted temperatures under plausible future conditions based on a warming scenario from the downscaled global climate model (MIROC A1B scenario). We did not find evidence of thermal local adaptation during the embryonic stage for developmental rate or survival. Within treatments, populations hatched within 1 day of each other, on average, and among treatments, did not differ in survival in response to temperature. We did detect plasticity to temperature; embryos developed 2.5 times longer (189 days) in the coolest regime compared to the warmest regime (74 days). We also detected variation in developmental rates among families within and among temperature regimes, indicating heritable plasticity. Families exhibited a strong positive relationship between thermal variability and phenotypic variability in developmental rate but body length and mass at hatching were largely insensitive to temperature. Overall, our results indicated a lack of thermal local adaptation, but a presence of plasticity in populations experiencing contrasting conditions, as well as family‐specific heritable plasticity that could facilitate adaptive change.     

Chromosomal fusion and life history‐associated genomic variation contribute to within‐river local adaptation of Atlantic salmon

Chromosomal inversions have been implicated in facilitating adaptation in the face of high levels of gene flow, but whether chromosomal fusions also have similar potential remains poorly understood. Atlantic salmon are usually characterized by population structure at multiple spatial scales; however, this is not the case for tributaries of the Miramichi River in North America. To resolve genetic relationships between populations in this system and the potential for known chromosomal fusions to contribute to adaptation, we genotyped 728 juvenile salmon using a 50 K SNP array. Consistent with previous work, we report extremely weak overall population structuring (Global F_ST = 0.004) and failed to support hierarchical structure between the river's two main branches. We provide the first genomic characterization of a previously described polymorphic fusion between chromosomes 8 and 29. Fusion genomic characteristics included high LD, reduced heterozygosity in the fused homokaryotes, and strong divergence between the fused and the unfused rearrangement. Population structure based on fusion karyotype was five times stronger than neutral variation (F_ST = 0.019), and the frequency of the fusion was associated with summer precipitation supporting a hypothesis that this rearrangement may contribute local adaptation despite weak neutral differentiation. Additionally, both outlier variation among populations and a polygenic framework for characterizing adaptive variation in relation to climate identified a 250‐Kb region of chromosome 9, including the gene six6 that has previously been linked to age‐at‐maturity and run‐timing for this species. Overall, our results indicate that adaptive processes, independent of major river branching, are more important than neutral processes for structuring these populations.     

Temporal and spatial variation in growth of juvenile Atlantic salmon

Spatial and temporal variation in length‐at‐age and environmental factors affecting variation in growth rate of juvenile Atlantic salmon Salmo salar were studied using data from a long‐term study in the River Stjørdalselva, central Norway. Mean annual instantaneous growth rate among 1+ and 2+year juvenile Atlantic salmon varied between 0·59 and 1·50 g g⁻¹ year⁻¹ and mean instantaneous daily growth rate of young‐of‐the‐year (YOY) varied between 0·013 and 0·033 g g⁻¹ day⁻¹. Between year variation in growth was larger than the within year intra‐watercourse spatial variation. For YOY and 1+year Atlantic salmon, a major part of the observed between year variation in growth rates was explained by variation in mean daily water discharge and spring temperature. For 2+year juvenile Atlantic salmon, mean daily water discharge and cohort density were the only variables to significantly explain variation in growth rates. A large part of the within water‐course spatial variation could not be explained by temperature variations and juvenile Atlantic salmon in the uppermost areas of the river, experiencing the lowest ambient temperatures during the growth period, displayed the highest growth rates. Within the baselines set by temperature, biotic and abiotic factors connected to water flow regime and variation in food availability are suggested to be a major determinants of the temporal and spatial variation in juvenile Atlantic salmon growth rates.     

Chromosome polymorphisms track trans‐Atlantic divergence and secondary contact in Atlantic salmon

Pleistocene glaciations drove repeated range contractions and expansions shaping contemporary intraspecific diversity. Atlantic salmon (Salmo salar) in the western and eastern Atlantic diverged >600,000 years before present, with the two lineages isolated in different southern refugia during glacial maxima, driving trans‐Atlantic genomic and karyotypic divergence. Here, we investigate the genomic consequences of glacial isolation and trans‐Atlantic secondary contact using 108,870 single nucleotide polymorphisms genotyped in 80 North American and European populations. Throughout North America, we identified extensive interindividual variation and discrete linkage blocks within and between chromosomes with known trans‐Atlantic differences in rearrangements: Ssa01/Ssa23 translocation and Ssa08/Ssa29 fusion. Spatial genetic analyses suggest independence of rearrangements, with Ssa01/Ssa23 showing high European introgression (>50%) in northern populations indicative of post‐glacial trans‐Atlantic secondary contact, contrasting with low European ancestry genome‐wide (3%). Ssa08/Ssa29 showed greater intrapopulation diversity, suggesting a derived chromosome fusion polymorphism that evolved within North America. Evidence of potential selection on both genomic regions suggests that the adaptive role of rearrangements warrants further investigation in Atlantic salmon. Our study highlights how Pleistocene glaciations can influence large‐scale intraspecific variation in genomic architecture of northern species.