Landscape features along migratory routes influence adaptive genomic variation in anadromous steelhead (Oncorhynchus mykiss)

SNP Polymorphisms Are Associated with Environmental Factors in Sockeye Salmon Populations Across the Northwest Pacific: Insights from Redundancy Analysis

The SNP variation in sockeye salmon across the Asian part of its range was studied in 23 samples from 16 lake-river systems of the West Pacific Coast to improve understanding of genetic adaptation in response to spawning watersheds conditions. Identification of candidate SNPs and environmental factors that can contribute to local adaptations in sockeye salmon populations was carried out using redundancy analysis (RDA), a powerful tool for landscape genetics proven to be effective in genotype-environment association studies. Climatic and hydrographic indices (7 indices in total), reflecting abiotic conditions in freshwater habitats of sockeye salmon and characterizing the temperature regime in the river basin, its variability during the year, the amount of precipitation, as well as the height of the maximum tide in the estuary, were used as predictor factors. Among the 45 analyzed SNPs, several loci (ALDOB-135, HGFA, and RAG3-93) correlated with predictors gradients along the northwest Pacific coast were identified. The putative candidate loci localized in genes involved in the immune and inflammatory responses, as well as genes encoding temperature-sensitive enzymes and some hormones regulating ion homeostasis in fish during the anadromous migration and smoltification, were potentially associated with environmental conditions in natal rivers. The findings could have implications for aquaculture, conservation, and resource management in the context of global climate change.

Historical and contemporary processes driving the origin and structure of an admixed population within a contact zone between subspecies of a north temperate diadromous fish

Hybridization between divergent lineages can result in losses of distinct evolutionary taxa. Alternatively, hybridization can lead to increased genetic variability that may fuel local adaptation and the generation of novel traits and/or taxa. Here, we examined single-nucleotide polymorphisms generated using genotyping-by-sequencing in a population of Dolly Varden char (Pisces: Salmonidae) that is highly admixed within a contact zone between two subspecies (Salvelinus malma malma, Northern Dolly Varden [NDV] and S. m. lordi, Southern Dolly Varden [SDV]) in southwestern Alaska to assess the spatial distribution of hybrids and to test hypotheses on the origin of the admixed population. Ancestry analysis revealed that this admixed population is composed of advanced generation hybrids between NDV and SDV or advanced backcrosses to SDV; no F1 hybrids were detected. Coalescent-based demographic modelling supported the origin of this population about 55,000 years ago by secondary contact between NDV and SDV with low levels of contemporary gene flow. Ancestry in NDV and SDV varies within the watershed and ancestry in NDV was positively associated with distance upstream from the sea, contingent on habitat-type sampled, and negatively associated with the number of migrations that individual fish made to the sea. Our results suggest that divergence between subspecies over hundreds of thousands of years may not be associated with significant reproductive isolation, but that elevated diversity owing to hybridization may have contributed to adaptive divergence in habitat use and life history.

Genomic signatures of climate adaptation in bank voles

Evidence for divergent selection and adaptive variation across the landscape can provide insight into a species' ability to adapt to different environments. However, despite recent advances in genomics, it remains difficult to detect the footprints of climate-mediated selection in natural populations. Here, we analysed ddRAD sequencing data (21,892 SNPs) in conjunction with geographic climate variation to search for signatures of adaptive differentiation in twelve populations of the bank vole (Clethrionomys glareolus) distributed across Europe. To identify the loci subject to selection associated with climate variation, we applied multiple genotype-environment association methods, two univariate and one multivariate, and controlled for the effect of population structure. In total, we identified 213 candidate loci for adaptation, 74 of which were located within genes. In particular, we identified signatures of selection in candidate genes with functions related to lipid metabolism and the immune system. Using the results of redundancy analysis, we demonstrated that population history and climate have joint effects on the genetic variation in the pan-European metapopulation. Furthermore, by examining only candidate loci, we found that annual mean temperature is an important factor shaping adaptive genetic variation in the bank vole. By combining landscape genomic approaches, our study sheds light on genome-wide adaptive differentiation and the spatial distribution of variants underlying adaptive variation influenced by local climate in bank voles.

Efficient species identification for Pacific salmon genetic monitoring programs

Genetic monitoring of Pacific salmon in the Columbia River basin provides crucial information to fisheries managers that is otherwise challenging to obtain using traditional methods. Monitoring programs such as genetic stock identification (GSI) and parentage-based tagging (PBT) involve genotyping tens of thousands of individuals annually. Although rare, these large sample collections inevitably include misidentified species, which exhibit low genotyping success on species-specific Genotyping-in-Thousands by sequencing (GT-seq) panels. For laboratories involved in large-scale genotyping efforts, diagnosing non-target species and reassigning them to the appropriate monitoring program can be costly and time-consuming. To address this problem, we identified 19 primer pairs that exhibit consistent cross-species amplification among salmonids and contain 51 species informative variants. These genetic markers reliably discriminate among 11 salmonid species and two subspecies of Cutthroat Trout and have been included in species-specific GT-seq panels for Chinook Salmon, Coho Salmon, Sockeye Salmon, and Rainbow Trout commonly used for Pacific salmon genetic monitoring. The majority of species-informative amplicons (16) were newly identified from the four existing GT-seq panels, thus demonstrating a low-cost approach to species identification when using targeted sequencing methods. A species-calling script was developed that is tailored for routine GT-seq genotyping pipelines and automates the identification of non-target species. Following extensive testing with empirical and simulated data, we demonstrated that the genetic markers and accompanying script accurately identified species and are robust to missing genotypic data and low-frequency, shared polymorphisms among species. Finally, we used these tools to identify Coho Salmon incidentally caught in the Columbia River Chinook Salmon sport fishery and used PBT to determine their hatchery of origin. These molecular and computing resources provide a valuable tool for Pacific salmon conservation in the Columbia River basin and demonstrate a cost-effective approach to species identification for genetic monitoring programs.

Genetic variation associated with adult migration timing in lineages of Steelhead and Chinook Salmon in the Columbia River

With the discovery of a major effect region (GREB1L, ROCK1) for adult migration timing in genomes of both Chinook Salmon and Steelhead, several subsequent studies have investigated the effect size and distribution of early and late migration alleles among populations in the Columbia River. Here, we synthesize the results of these studies for the major lineages of Chinook Salmon and Steelhead that include highly distinct groups in the interior Columbia River that exhibit atypical life histories from most coastal lineage populations of these two species. Whole-genome studies with high marker density have provided extensive insight into SNPs most associated with adult migration timing, and suites of markers for each species have been genotyped in large numbers of individuals to further validate phenotypic effects. For Steelhead, the largest phenotypic effect sizes have been observed in the coastal lineage (36% of variation for passage timing at Bonneville Dam; 43% of variation for tributary arrival timing) compared to the inland lineage (7.5% of variation for passage timing at Bonneville Dam; 8.4% of variation for tributary arrival timing) that overwinter in freshwater prior to spawning. For Chinook Salmon, large effect sizes have been observed in all three lineages for multiple adult migration phenotypes (Coastal lineage: percentage of variation of 27.9% for passage timing at Bonneville Dam, 28.7% for arrival timing for spawning; Interior ocean type: percentage of variation of 47.6% for passage timing at Bonneville Dam, 39.6% for tributary arrival timing, 77.9% for arrival timing for spawning; Interior stream type: percentage of variation of 35.3% for passage at Bonneville Dam, 9.8% for tributary arrival timing, 4.7% for arrival timing for spawning). Together, these results have extended our understanding of genetic variation associated with life history diversity in distinct populations of the Columbia River, however, much research remains necessary to determine the causal mechanism for this major effect region on migration timing in these species.

Contrasting patterns of sequence variation in steelhead populations reflect distinct evolutionary processes

Multiple evolutionary processes influence genome-wide allele frequencies and quantifying effects of genetic drift, and multiple forms of selection remain challenging in natural populations. Here, we investigate variation at major effect loci in contrast to patterns of neutral drift across a wide collection of steelhead (Oncorhynchus mykiss) populations that have declined in abundance due to anthropogenic impacts. Whole-genome resequencing of 74 populations of steelhead revealed genome-wide patterns (~8 million SNPs) consistent with expected neutral population structure. However, allelic variation at major effect loci associated with adult migration timing (chromosome 28: GREB1L/ROCK1) and age at maturity (chromosome 25: SIX6) reflected how selection has acted on phenotypic variation in contrast with neutral structure. Variation at major effect loci was influenced by evolutionary processes with differing signals between the strongly divergent Coastal and Inland lineages, while allele frequencies within and among populations within the Inland lineage have been driven by local natural selection as well as recent anthropogenic influences. Recent anthropogenic effects appeared to have influenced the frequency of major effect alleles including artificial selection for specific traits in hatchery stocks with subsequent gene flow into natural populations. Selection from environmental factors at various scales has also likely influenced variation for major effect alleles. These results reveal evolutionary mechanisms that influence allele frequencies at major effect loci that are critical for conservation of phenotypic traits and life history variation of this protected species.

Intraspecific variation among Chinook Salmon populations indicates physiological adaptation to local environmental conditions

Understanding interpopulation variation is important to predicting species responses to climate change. Recent research has revealed interpopulation variation among several species of Pacific salmonids; however, the environmental drivers of population differences remain elusive. We tested for local adaptation and countergradient variation by assessing interpopulation variation among six populations of fall-run Chinook Salmon from the western United States. Juvenile fish were reared at three temperatures (11, 16 and 20°C), and five physiological metrics were measured (routine and maximum metabolic rate, aerobic scope, growth rate and critical thermal maximum). We then tested associations between these physiological metrics and 15 environmental characteristics (e.g. rearing temperature, latitude, migration distance, etc.). Statistical associations between the five physiological metrics and 15 environmental characteristics supported our hypotheses of local adaptation. Notably, latitude was a poor predictor of population physiology. Instead, our results demonstrate that populations from warmer habitats exhibit higher thermal tolerance (i.e. critical thermal maxima), faster growth when warm acclimated and greater aerobic capacity at high temperatures. Additionally, populations with longer migrations exhibit higher metabolic capacity. However, overall metabolic capacity declined with warm acclimation, indicating that future climate change may reduce metabolic capacity, negatively affecting long-migrating populations. Linking physiological traits to environmental characteristics enables flexible, population-specific management of disparate populations in response to local conditions.

Long-distance migration is a major factor driving local adaptation at continental scale in Coho salmon

Inferring the genomic basis of local adaptation is a long-standing goal of evolutionary biology. Beyond its fundamental evolutionary implications, such knowledge can guide conservation decisions for populations of conservation and management concern. Here, we investigated the genomic basis of local adaptation in the Coho salmon (Oncorhynchus kisutch) across its entire North American range. We hypothesized that extensive spatial variation in environmental conditions and the species' homing behaviour may promote the establishment of local adaptation. We genotyped 7829 individuals representing 217 sampling locations at more than 100,000 high-quality RADseq loci to investigate how recombination might affect the detection of loci putatively under selection and took advantage of the precise description of the demographic history of the species from our previous work to draw accurate population genomic inferences about local adaptation. The results indicated that genetic differentiation scans and genetic-environment association analyses were both significantly affected by variation in recombination rate as low recombination regions displayed an increased number of outliers. By taking these confounding factors into consideration, we revealed that migration distance was the primary selective factor driving local adaptation and partial parallel divergence among distant populations. Moreover, we identified several candidate single nucleotide polymorphisms associated with long-distance migration and altitude including a gene known to be involved in adaptation to altitude in other species. The evolutionary implications of our findings are discussed along with conservation applications.

Whole genome resequencing identifies local adaptation associated with environmental variation for redband trout

Aquatic ectotherms are predicted to harbour genomic signals of local adaptation resulting from selective pressures driven by the strong influence of climate conditions on body temperature. We investigated local adaptation in redband trout (Oncorhynchus mykiss gairdneri) using genome scans for 547 samples from 11 populations across a wide range of habitats and thermal gradients in the interior Columbia River. We estimated allele frequencies for millions of single nucleotide polymorphism loci (SNPs) across populations using low-coverage whole genome resequencing, and used population structure outlier analyses to identify genomic regions under divergent selection between populations. Twelve genomic regions showed signatures of local adaptation, including two regions associated with genes known to influence migration and developmental timing in salmonids (GREB1L, ROCK1, SIX6). Genotype-environment association analyses indicated that diurnal temperature variation was a strong driver of local adaptation, with signatures of selection driven primarily by divergence of two populations in the northern extreme of the subspecies range. We also found evidence for adaptive differences between high-elevation desert vs. montane habitats at a smaller geographical scale. Finally, we estimated vulnerability of redband trout to future climate change using ecological niche modelling and genetic offset analyses under two climate change scenarios. These analyses predicted substantial habitat loss and strong genetic shifts necessary for adaptation to future habitats, with the greatest vulnerability predicted for high-elevation desert populations. Our results provide new insight into the complexity of local adaptation in salmonids, and important predictions regarding future responses of redband trout to climate change.

Population genomic monitoring provides insight into conservation status but no correlation with demographic estimates of extinction risk in a threatened trout

The current extinction crisis requires effective assessment and monitoring tools. Genetic approaches are appealing given the relative ease of field sampling required to estimate genetic diversity characteristics assumed related to population size, evolutionary potential, and extinction risk, and to evaluate hybridization with non-native species simultaneously. However, linkages between population genetic metrics of diversity from survey-style field collections and demographic estimates of population size and extinction risk are still in need of empirical examples, especially for remotely distributed species of conservation concern where the approach might be most beneficial. We capitalized on an exceptional opportunity to evaluate congruence between genetic diversity metrics and demographic-based estimates of abundance and extinction risk from a comprehensive Multiple Population Viability Analysis (MPVA) in a threatened fish, the Lahontan cutthroat trout (LCT). We sequenced non-native trout reference samples and recently collected and archived tissue samples of most remaining populations of LCT (N = 60) and estimated common genetic assessment metrics, predicting minimal hybridization with non-native trout, low diversity, and declining diversity over time. We further hypothesized genetic metrics would correlate positively with MPVA-estimated abundance and negatively with extinction probability. We uncovered several instances of hybridization that pointed to immediate management needs. After removing hybridized individuals, cautious interpretation of low effective population sizes (2-63) suggested reduced evolutionary potential for many LCT populations. Other genetic metrics did not decline over time nor correlate with MPVA-based estimates of harmonic mean abundance or 30-year extinction probability. Our results demonstrate benefits of genetic monitoring for efficiently detecting hybridization and, though genetic results were disconnected from demographic assessment of conservation status, they suggest reduced evolutionary potential and likely a higher conservation risk than currently recognized for this threatened fish. We emphasize that genetic information provides essential complementary insight, in addition to demographic information, for evaluating species status.

Influence of environmental conditions at spawning sites and migration routes on adaptive variation and population connectivity in Chinook salmon

Many species that undergo long breeding migrations, such as anadromous fishes, face highly heterogeneous environments along their migration corridors and at their spawning sites. These environmental challenges encountered at different life stages may act as strong selective pressures and drive local adaptation. However, the relative influence of environmental conditions along the migration corridor compared with the conditions at spawning sites on driving selection is still unknown. In this study, we performed genome-environment associations (GEA) to understand the relationship between landscape and environmental conditions driving selection in seven populations of the anadromous Chinook salmon (Oncorhynchus tshawytscha)-a species of important economic, social, cultural, and ecological value-in the Columbia River basin. We extracted environmental variables for the shared migration corridors and at distinct spawning sites for each population, and used a Pool-seq approach to perform whole genome resequencing. Bayesian and univariate GEA tests with migration-specific and spawning site-specific environmental variables indicated many more candidate SNPs associated with environmental conditions at the migration corridor compared with spawning sites. Specifically, temperature, precipitation, terrain roughness, and elevation variables of the migration corridor were the most significant drivers of environmental selection. Additional analyses of neutral loci revealed two distinct clusters representing populations from different geographic regions of the drainage that also exhibit differences in adult migration timing (summer vs. fall). Tests for genomic regions under selection revealed a strong peak on chromosome 28, corresponding to the GREB1L/ROCK1 region that has been identified previously in salmonids as a region associated with adult migration timing. Our results show that environmental variation experienced throughout migration corridors imposed a greater selective pressure on Chinook salmon than environmental conditions at spawning sites.

Genomic data support management of anadromous Arctic Char fisheries in Nunavik by highlighting neutral and putatively adaptive genetic variation

Distinguishing neutral and adaptive genetic variation is one of the main challenges in investigating processes shaping population structure in the wild, and landscape genomics can help identify signatures of adaptation to contrasting environments. Arctic Char (Salvelinus alpinus) is an anadromous salmonid and the most harvested fish species by Inuit people, including in Nunavik (Québec, Canada), one of the most recently deglaciated regions in the world. Unlike many other anadromous salmonids, Arctic Char occupy coastal habitats near their natal rivers during their short marine phase restricted to the summer ice-free period. Our main objective was to document putatively neutral and adaptive genomic variation in anadromous Arctic Char populations from Nunavik and bordering regions to inform local fisheries management. We used genotyping by sequencing (GBS) to genotype 18,112 filtered single nucleotide polymorphisms (SNP) in 650 individuals from 23 sampling locations along >2000 km of coastline. Our results reveal a hierarchical genetic structure, whereby neighboring hydrographic systems harbor distinct populations grouped by major oceanographic basins: Hudson Bay, Hudson Strait, Ungava Bay, and Labrador Sea. We found genetic diversity and differentiation to be consistent both with the expected postglacial recolonization history and with patterns of isolation-by-distance reflecting contemporary gene flow. Results from three gene-environment association methods supported the hypothesis of local adaptation to both freshwater and marine environments (strongest associations with sea surface and air temperatures during summer and salinity). Our results support a fisheries management strategy at a regional scale, and other implications for hatchery projects and adaptation to climate change are discussed.

The Impacts of Dam Construction and Removal on the Genetics of Recovering Steelhead (Oncorhynchus mykiss) Populations across the Elwha River Watershed

Dam construction and longitudinal river habitat fragmentation disrupt important life histories and movement of aquatic species. This is especially true for Oncorhynchus mykiss that exhibits both migratory (steelhead) and non-migratory (resident rainbow) forms. While the negative effects of dams on salmonids have been extensively documented, few studies have had the opportunity to compare population genetic diversity and structure prior to and following dam removal. Here we examine the impacts of the removal of two dams on the Elwha River on the population genetics of O. mykiss. Genetic data were produced from >1200 samples collected prior to dam removal from both life history forms, and post-dam removal from steelhead. We identified three genetic clusters prior to dam removal primarily explained by isolation due to dams and natural barriers. Following dam removal, genetic structure decreased and admixture increased. Despite large O. mykiss population declines after dam construction, we did not detect shifts in population genetic diversity or allele frequencies of loci putatively involved in migratory phenotypic variation. Steelhead descendants from formerly below and above dammed populations recolonized the river rapidly after dam removal, suggesting that dam construction did not significantly reduce genetic diversity underlying O. mykiss life history strategies. These results have significant evolutionary implications for the conservation of migratory adaptive potential in O. mykiss populations above current anthropogenic barriers.

Investigating Diadromy in Fishes and Its Loss in an -Omics Era

Diadromy, the predictable movements of individuals between marine and freshwater environments, is biogeographically and phylogenetically widespread across fishes. Thus, despite the high energetic and potential fitness costs involved in moving between distinct environments, diadromy appears to be an effective life history strategy. Yet, the origin and molecular mechanisms that underpin this migratory behavior are not fully understood. In this review, we aim first to summarize what is known about diadromy in fishes; this includes the phylogenetic relationship among diadromous species, a description of the main hypotheses regarding its origin, and a discussion of the presence of non-migratory populations within diadromous species. Second, we discuss how recent research based on -omics approaches (chiefly genomics, transcriptomics, and epigenomics) is beginning to provide answers to questions on the genetic bases and origin(s) of diadromy. Finally, we suggest future directions for -omics research that can help tackle questions on the evolution of diadromy.

Steelhead (Oncorhynchus mykiss) lineages and sexes show variable patterns of association of adult migration timing and age-at-maturity traits with two genomic regions

As life history diversity plays a critical role in supporting the resilience of exploited populations, understanding the genetic basis of those life history variations is important for conservation management. However, effective application requires a robust understanding of the strength and universality of genetic associations. Here, we examine genetic variation of single nucleotide polymorphisms in genomic regions previously associated with migration phenology and age-at-maturity in steelhead (Oncorhynchus mykiss) from the Columbia River. We found chromosome 28 markers (GREB1L, ROCK1 genes) explained significant variance in migration timing in both coastal and inland steelhead. However, strength of association was much greater in coastal than inland steelhead (R 2 0.51 vs. 0.08), suggesting that genomic background and challenging inland migration pathways may act to moderate effects of this region. Further, we found that chromosome 25 candidate markers (SIX6 gene) were significantly associated with age and size at first return migration for inland steelhead, and this pattern was mediated by sex in a predictable pattern (males R 2 = 0.139-0.170; females R 2 = 0.096-0.111). While this encourages using these candidate regions in predicting life history characteristics, we suggest that stock-specific associations and haplotype frequencies will be useful in guiding implementation of genetic assays to inform management.

Distribution of genetic variation underlying adult migration timing in steelhead of the Columbia River basin

Fish migrations are energetically costly, especially when moving between freshwater and saltwater, but are a viable strategy for Pacific salmon and trout (Oncorhynchus spp.) due to the advantageous resources available at various life stages. Anadromous steelhead (O. mykiss) migrate vast distances and exhibit variation for adult migration phenotypes that have a genetic basis at candidate genes known as greb1L and rock1. We examined the distribution of genetic variation at 13 candidate markers spanning greb1L, intergenic, and rock1 regions versus 226 neutral markers for 113 populations (n = 9,471) of steelhead from inland and coastal lineages in the Columbia River. Patterns of population structure with neutral markers reflected genetic similarity by geographic region as demonstrated in previous studies, but candidate markers clustered populations by genetic variation associated with adult migration timing. Mature alleles for late migration had the highest frequency overall in steelhead populations throughout the Columbia River, with only 9 of 113 populations that had a higher frequency of premature alleles for early migration. While a single haplotype block was evident for the coastal lineage, we identified multiple haplotype blocks for the inland lineage. The inland lineage had one haplotype block that corresponded to candidate markers within the greb1L gene and immediately upstream in the intergenic region, and the second block only contained candidate markers from the intergenic region. Haplotype frequencies had similar patterns of geographic distribution as single markers, but there were distinct differences in frequency between the two haplotype blocks for the inland lineage. This may represent multiple recombination events that differed between lineages where phenotypic differences exist between freshwater entry versus arrival timing as indicated by Micheletti et al. (2018a). Redundancy analyses were used to model environmental effects on allelic frequencies of candidate markers, and significant variables were migration distance, temperature, isothermality, and annual precipitation. This study improves our understanding of the spatial distribution of genetic variation underlying adult migration timing in steelhead as well as associated environmental factors and has direct conservation and management implications.

Temporal dynamics of migration-linked genetic variation are driven by streamflows and riverscape permeability

Landscape permeability is often explored spatially, but may also vary temporally. Landscape permeability, including partial barriers, influences migratory animals that move across the landscape. Partial barriers are common in rivers where barrier passage varies with streamflow. We explore the influence of partial barriers on the spatial and temporal distribution of migration‐linked genotypes of Oncorhynchus mykiss, a salmonid fish with co‐occurring resident and migratory forms, in tributaries to the South Fork Eel River, California, USA, Elder and Fox Creeks. We genotyped >4,000 individuals using RAD‐capture and classified individuals as resident, heterozygous or migratory genotypes using life history‐associated loci. Across four years of study (2014–2017), the permeability of partial barriers varied across dry and wet years. In Elder Creek, the largest waterfall was passable for adults migrating up‐river 4–39 days each year. In this stream, the overall spatial pattern, with fewer migratory genotypes above the waterfall, remained true across dry and wet years (67%–76% of migratory alleles were downstream of the waterfall). We also observed a strong relationship between distance upstream and proportion of migratory alleles. In Fox Creek, the primary barrier is at the mouth, and we found that the migratory allele frequency varied with the annual timing of high flow events. In years when rain events occurred during the peak breeding season, migratory allele frequency was high (60%–68%), but otherwise it was low (30% in two years). We highlight that partial barriers and landscape permeability can be temporally dynamic, and this effect can be observed through changing genotype frequencies in migratory animals.

Anadromy, potamodromy and residency in brown trout Salmo trutta: the role of genes and the environment

Brown trout Salmo trutta is endemic to Europe, western Asia and north-western Africa; it is a prominent member of freshwater and coastal marine fish faunas. The species shows two resident (river-resident, lake-resident) and three main facultative migratory life histories (downstream-upstream within a river system, fluvial-adfluvial potamodromous; to and from a lake, lacustrine-adfluvial (inlet) or allacustrine (outlet) potamodromous; to and from the sea, anadromous). River-residency v. migration is a balance between enhanced feeding and thus growth advantages of migration to a particular habitat v. the costs of potentially greater mortality and energy expenditure. Fluvial-adfluvial migration usually has less feeding improvement, but less mortality risk, than lacustrine-adfluvial or allacustrine and anadromous, but the latter vary among catchments as to which is favoured. Indirect evidence suggests that around 50% of the variability in S. trutta migration v. residency, among individuals within a population, is due to genetic variance. This dichotomous decision can best be explained by the threshold-trait model of quantitative genetics. Thus, an individual's physiological condition (e.g., energy status) as regulated by environmental factors, genes and non-genetic parental effects, acts as the cue. The magnitude of this cue relative to a genetically predetermined individual threshold, governs whether it will migrate or sexually mature as a river-resident. This decision threshold occurs early in life and, if the choice is to migrate, a second threshold probably follows determining the age and timing of migration. Migration destination (mainstem river, lake, or sea) also appears to be genetically programmed. Decisions to migrate and ultimate destination result in a number of subsequent consequential changes such as parr-smolt transformation, sexual maturity and return migration. Strong associations with one or a few genes have been found for most aspects of the migratory syndrome and indirect evidence supports genetic involvement in all parts. Thus, migratory and resident life histories potentially evolve as a result of natural and anthropogenic environmental changes, which alter relative survival and reproduction. Knowledge of genetic determinants of the various components of migration in S. trutta lags substantially behind that of Oncorhynchus mykiss and other salmonines. Identification of genetic markers linked to migration components and especially to the migration-residency decision, is a prerequisite for facilitating detailed empirical studies. In order to predict effectively, through modelling, the effects of environmental changes, quantification of the relative fitness of different migratory traits and of their heritabilities, across a range of environmental conditions, is also urgently required in the face of the increasing pace of such changes.

Genomics and conservation units: The genetic basis of adult migration timing in Pacific salmonids

It is now routinely possible to generate genomics-scale datasets for nonmodel species; however, many questions remain about how best to use these data for conservation and management. Some recent genomics studies of anadromous Pacific salmonids have reported a strong association between alleles at one or a very few genes and a key life history trait (adult migration timing) that has played an important role in defining conservation units. Publication of these results has already spurred a legal challenge to the existing framework for managing these species, which was developed under the paradigm that most phenotypic traits are controlled by many genes of small effect, and that parallel evolution of life history traits is common. But what if a key life history trait can only be expressed if a specific allele is present? Does the current framework need to be modified to account for the new genomics results, as some now propose? Although this real-world example focuses on Pacific salmonids, the issues regarding how genomics can inform us about the genetic basis of phenotypic traits, and what that means for applied conservation, are much more general. In this perspective, we consider these issues and outline a general process that can be used to help generate the types of additional information that would be needed to make informed decisions about the adequacy of existing conservation and management frameworks.

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

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

Selection at a genomic region of major effect is responsible for evolution of complex life histories in anadromous steelhead

BACKGROUND

Disparity in the timing of biological events occurs across a variety of systems, yet the understanding of genetic basis underlying diverse phenologies remains limited. Variation in maturation timing occurs in steelhead trout, which has been associated with greb1L, an oestrogen target gene. Previous techniques that identified this gene only accounted for about 0.5-2.0% of the genome and solely investigated coastal populations, leaving uncertainty on the genetic basis of this trait and its prevalence across a larger geographic scale.

RESULTS

We used a three-tiered approach to interrogate the genomic basis of complex phenology in anadromous steelhead. First, fine scale mapping with 5.3 million SNPs from resequencing data covering 68% of the genome confirmed a 309-kb region consisting of four genes on chromosome 28, including greb1L, to be the genomic region of major effect for maturation timing. Second, broad-scale characterization of candidate greb1L genotypes across 59 populations revealed unexpected patterns in maturation phenology for inland fish migrating long distances relative to those in coastal streams. Finally, genotypes from 890 PIT-tag tracked steelhead determined associations with early versus late arrival to spawning grounds that were previously unknown.

CONCLUSIONS

This study clarifies the genetic bases for disparity in phenology observed in steelhead, determining an unanticipated trait association with premature versus mature arrival to spawning grounds and identifying multiple candidate genes potentially contributing to this variation from a single genomic region of major effect. This illustrates how dense genome mapping and detailed phenotypic characterization can clarify genotype to phenotype associations across geographic ranges of species.

Utility of pooled sequencing for association mapping in nonmodel organisms

High-density genome-wide sequencing increases the likelihood of discovering genes of major effect and genomic structural variation in organisms. While there is an increasing availability of reference genomes across broad taxa, the greatest limitation to whole-genome sequencing of multiple individuals continues to be the costs associated with sequencing. To alleviate excessive costs, pooling multiple individuals with similar phenotypes and sequencing the homogenized DNA (Pool-Seq) can achieve high genome coverage, but at the loss of individual genotypes. Although Pool-Seq has been an effective method for association mapping in model organisms, it has not been frequently utilized in natural populations. To extend bioinformatic tools for rapid implementation of Pool-Seq data in nonmodel organisms, we developed a pipeline called PoolParty and illustrate its effectiveness in genetic association mapping. Alignment expectations based on five pooled Chinook salmon (Oncorhynchus tshawytscha) libraries showed that approximately 48% genome coverage per library could be achieved with reasonable sequencing effort. We additionally examined male and female O. tshawytscha libraries to illustrate how Pool-Seq techniques can successfully map known genes associated with functional differences among sexes such as growth hormone 2. Finally, we compared pools of individuals of different spawning ages for each sex to discover novel genes involved with age at maturity in O. tshawytscha such as opsin4 and transmembrane protein19. While not appropriate for every system, Pool-Seq data processed by the PoolParty pipeline is a practical method for identifying genes of major effect in nonmodel organisms when high genome coverage is necessary and cost is a limiting factor.