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

Organisms typically show evidence of adaptation to features within their local environment. However, many species undergo long-distance dispersal or migration across larger geographic regions that consist of highly heterogeneous habitats. Therefore, selection may influence adaptive genetic variation associated with landscape features at residing sites and along migration routes in migratory species. We tested for genomic adaptation to landscape features at natal spawning sites and along migration paths to the ocean of anadromous steelhead trout (Oncorhynchus mykiss) in the Columbia River Basin. Results from multivariate ordination, gene-environment association and outlier analyses using 24,526 single nucleotide polymorphisms (SNPs) provided evidence that adaptive allele frequencies were more commonly associated with landscape features along migration paths than features at natal sites (91.8% vs. 8.2% of adaptive loci, respectively). Among the 45 landscape variables tested, migration distance to the ocean and mean annual precipitation along migration paths were significantly associated with adaptive genetic variation in three distinct genetic groups. Additionally, variables such as minimum migration water temperature and mean migration slope were significant only in inland stocks of steelhead that migrate up to 1,200 km farther than those near the coast, indicating regional differences in migratory selective pressures. This study provides novel approaches for investigating migratory corridors and some of the first evidence that environment along migration paths can lead to substantial divergent selection. Consequently, our approach to understand genetic adaptation to migration conditions can be applied to other migratory species when migration or dispersal paths are generally known.

Restricted gene flow between resident Oncorhynchus mykiss and an admixed population of anadromous steelhead.

The species Oncorhynchus mykiss is characterized by a complex life history that presents a significant challenge for population monitoring and conservation management. Many factors contribute to genetic variation in O. mykiss populations, including sympatry among migratory phenotypes, habitat heterogeneity, hatchery introgression, and immigration (stray) rates. The relative influences of these and other factors are contingent on characteristics of the local environment. The Rock Creek subbasin in the middle Columbia River has no history of hatchery supplementation and no dams or artificial barriers. Limited intervention and minimal management have led to a dearth of information regarding the genetic distinctiveness of the extant O. mykiss population in Rock Creek and its tributaries. We used 192 SNP markers and collections sampled over a 5-year period to evaluate the temporal and spatial genetic structures of O. mykiss between upper and lower watersheds of the Rock Creek subbasin. We investigated potential limits to gene flow within the lower watershed where the stream is fragmented by seasonally dry stretches of streambed, and between upper and lower watershed regions. We found minor genetic differentiation within the lower watershed occupied by anadromous steelhead (FST = 0.004), and evidence that immigrant influences were prevalent and ubiquitous. Populations in the upper watershed above partial natural barriers were highly distinct (FST = 0.093) and minimally impacted by apparent introgression. Genetic structure between watersheds paralleled differences in local demographics (e.g., variation in size), migratory restrictions, and habitat discontinuity. The evidence of restricted gene flow between putative remnant resident populations in the upper watershed and the admixed anadromous population in the lower watershed has implications for local steelhead productivity and regional conservation.

Climate variables explain neutral and adaptive variation within salmonid metapopulations: the importance of replication in landscape genetics.

Understanding how environmental variation influences population genetic structure is important for conservation management because it can reveal how human stressors influence population connectivity, genetic diversity and persistence. We used riverscape genetics modelling to assess whether climatic and habitat variables were related to neutral and adaptive patterns of genetic differentiation (population-specific and pairwise FST ) within five metapopulations (79 populations, 4583 individuals) of steelhead trout (Oncorhynchus mykiss) in the Columbia River Basin, USA. Using 151 putatively neutral and 29 candidate adaptive SNP loci, we found that climate-related variables (winter precipitation, summer maximum temperature, winter highest 5% flow events and summer mean flow) best explained neutral and adaptive patterns of genetic differentiation within metapopulations, suggesting that climatic variation likely influences both demography (neutral variation) and local adaptation (adaptive variation). However, we did not observe consistent relationships between climate variables and FST across all metapopulations, underscoring the need for replication when extrapolating results from one scale to another (e.g. basin-wide to the metapopulation scale). Sensitivity analysis (leave-one-population-out) revealed consistent relationships between climate variables and FST within three metapopulations; however, these patterns were not consistent in two metapopulations likely due to small sample sizes (N = 10). These results provide correlative evidence that climatic variation has shaped the genetic structure of steelhead populations and highlight the need for replication and sensitivity analyses in land and riverscape genetics.

Environmental adaptation in Chinook salmon (Oncorhynchus tshawytscha) throughout their North American range.

Landscape genomics is a rapidly growing field with recent advances in both genotyping efficiency and statistical analyses that provide insight towards local adaptation of populations under varying environmental and selective pressure. Chinook salmon (Oncorhynchus tshawytscha) are a broadly distributed Pacific salmon species, occupying a diversity of habitats throughout the northeastern Pacific with pronounced variation in environmental and climate features but little is understood regarding local adaptation in this species. We used a multivariate method, redundancy analysis (RDA), to identify polygenic correlations between 19,703 SNP loci and a suite of environmental variables in 46 collections of Chinook salmon (1956 total individuals) distributed throughout much of its North American range. Models in RDA were conducted on both rangewide and regional scales by hierarchical partitioning of the populations into three distinct genetic lineages. Our results indicate that between 5.8 and 21.8% of genomic variation can be accounted for by environmental features, and 566 putatively adaptive loci were identified as targets of environmental adaptation. The most influential drivers of adaptive divergence included precipitation in the driest quarter of the year (Rangewide and North Coastal Lineage, anova P = 0.002 and 0.01, respectively), precipitation in the wettest quarter of the year (Interior Columbia River Stream-Type Lineage, anova P = 0.03), variation in mean diurnal range in temperature (South Coastal Lineage, ANOVA P = 0.005), and migration distance (Rangewide, anova P = 0.001). Our results indicate that environmental features are strong drivers of adaptive genomic divergence in this species, and provide a foundation to investigate how Chinook salmon might respond to global environmental change.

Relative contributions of neutral and non-neutral genetic differentiation to inform conservation of steelhead trout across highly variable landscapes.

Mounting evidence of climatic effects on riverine environments and adaptive responses of fishes have elicited growing conservation concerns. Measures to rectify population declines include assessment of local extinction risk, population ecology, viability, and genetic differentiation. While conservation planning has been largely informed by neutral genetic structure, there has been a dearth of critical information regarding the role of non-neutral or functional genetic variation. We evaluated genetic variation among steelhead trout of the Columbia River Basin, which supports diverse populations distributed among dynamic landscapes. We categorized 188 SNP loci as either putatively neutral or candidates for divergent selection (non-neutral) using a multitest association approach. Neutral variation distinguished lineages and defined broad-scale population structure consistent with previous studies, but fine-scale resolution was also detected at levels not previously observed. Within distinct coastal and inland lineages, we identified nine and 22 candidate loci commonly associated with precipitation or temperature variables and putatively under divergent selection. Observed patterns of non-neutral variation suggest overall climate is likely to shape local adaptation (e.g., potential rapid evolution) of steelhead trout in the Columbia River region. Broad geographic patterns of neutral and non-neutral variation demonstrated here can be used to accommodate priorities for regional management and inform long-term conservation of this species.