Visual and genetic stock identification of a test fishery to forecast Columbia River spring Chinook salmon stocks 2 weeks into the future

Modern fisheries management strives to balance opposing goals of protection for weak stocks and opportunity for harvesting healthy stocks. Test fisheries can aid management of anadromous fishes if they can forecast the strength and timing of an annual run with adequate time to allow fisheries planning. Integration of genetic stock identification (GSI) can further maximize utility of test fisheries by resolving run forecasts into weak- and healthy-stock subcomponents. Using 5 years (2017-2022) of test fishery data, our study evaluated accuracy, resolution, and lead time of predictions for stock-specific run timing and abundance of Columbia River spring Chinook salmon (Oncorhynchus tshawytscha). We determined if this test fishery (1) could use visual stock identification (VSI) to forecast at the coarse stock resolution (i.e., classification of "lower" vs. "upriver" stocks) upon which current management is based and (2) could be enhanced with GSI to forecast at higher stock resolution. VSI accurately identified coarse stocks (83.3% GSI concordance), and estimated a proxy for abundance (catch per unit effort, CPUE) of the upriver stock in the test fishery that was correlated (R 2 = 0.90) with spring Chinook salmon abundance at Bonneville dam (Rkm 235). Salmon travel rates (~8.6 Rkm/day) provided predictions with 2-week lead time prior to dam passage. Importantly, GSI resolved this predictive ability as finely as the hatchery broodstock level. Lower river stock CPUE in the test fishery was correlated with abundance at Willamette Falls (Rkm 196, R 2 = 0.62), but could not be as finely resolved as achieved for upriver stocks. We described steps to combine VSI and GSI to provide timely in-season information and with prediction accuracy of ~12.4 mean absolute percentage error and high stock resolution to help plan Columbia River mainstem fisheries.

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.

Local environments, not invasive hybridization, influence cardiac performance of native trout under acute thermal stress

Climate-induced expansion of invasive hybridization (breeding between invasive and native species) poses a significant threat to the persistence of many native species worldwide. In the northern U.S. Rocky Mountains, hybridization between native cutthroat trout and non-native rainbow trout has increased in recent decades due, in part, to climate-driven increases in water temperature. It has been postulated that invasive hybridization may enhance physiological tolerance to climate-induced thermal stress because laboratory studies indicate that rainbow trout have a higher thermal tolerance than cutthroat trout. Here, we assessed whether invasive hybridization improves cardiac performance response to acute water temperature stress of native wild trout populations. We collected trout from four streams with a wide range of non-native admixture among individuals and with different temperature and streamflow regimes in the upper Flathead River drainage, USA. We measured individual cardiac performance (maximum heart rate, "MaxHR", and temperature at arrhythmia, "ArrTemp") during laboratory trials with increasing water temperatures (10-28°C). Across the study populations, we observed substantial variation in cardiac performance of individual trout when exposed to thermal stress. Notably, we found significant differences in the cardiac response to thermal regimes among native cutthroat trout populations, suggesting the importance of genotype-by-environment interactions in shaping the physiological performance of native cutthroat trout. However, rainbow trout admixture had no significant effect on cardiac performance (MaxHR and ArrTemp) within any of the three populations. Our results indicate that invasive hybridization with a warmer-adapted species does not enhance the cardiac performance of native trout under warming conditions. Maintaining numerous populations across thermally and hydrologically diverse stream environments will be crucial for native trout to adapt and persist in a warming climate.

Genomics and conservation: Guidance from training to analyses and applications

Environmental change is intensifying the biodiversity crisis and threatening species across the tree of life. Conservation genomics can help inform conservation actions and slow biodiversity loss. However, more training, appropriate use of novel genomic methods and communication with managers are needed. Here, we review practical guidance to improve applied conservation genomics. We share insights aimed at ensuring effectiveness of conservation actions around three themes: (1) improving pedagogy and training in conservation genomics including for online global audiences, (2) conducting rigorous population genomic analyses properly considering theory, marker types and data interpretation and (3) facilitating communication and collaboration between managers and researchers. We aim to update students and professionals and expand their conservation toolkit with genomic principles and recent approaches for conserving and managing biodiversity. The biodiversity crisis is a global problem and, as such, requires international involvement, training, collaboration and frequent reviews of the literature and workshops as we do here.

Reproducibility in ecology and evolution: Minimum standards for data and code

We call for journals to commit to requiring open data be archived in a format that will be simple and clear for readers to understand and use. If applied consistently, these requirements will allow contributors to be acknowledged for their work through citation of open data, and facilitate scientific progress.

An improved germline genome assembly for the sea lamprey Petromyzon marinus illuminates the evolution of germline-specific chromosomes

Programmed DNA loss is a gene silencing mechanism that is employed by several vertebrate and nonvertebrate lineages, including all living jawless vertebrates and songbirds. Reconstructing the evolution of somatically eliminated (germline-specific) sequences in these species has proven challenging due to a high content of repeats and gene duplications in eliminated sequences and a corresponding lack of highly accurate and contiguous assemblies for these regions. Here, we present an improved assembly of the sea lamprey (Petromyzon marinus) genome that was generated using recently standardized methods that increase the contiguity and accuracy of vertebrate genome assemblies. This assembly resolves highly contiguous, somatically retained chromosomes and at least one germline-specific chromosome, permitting new analyses that reconstruct the timing, mode, and repercussions of recruitment of genes to the germline-specific fraction. These analyses reveal major roles of interchromosomal segmental duplication, intrachromosomal duplication, and positive selection for germline functions in the long-term evolution of germline-specific chromosomes.

Genomic variation across Chinook salmon populations reveals effects of a duplication on migration alleles and supports fine scale structure

The distribution of ecotypic variation in natural populations is influenced by neutral and adaptive evolutionary forces that are challenging to disentangle. This study provides a high-resolution portrait of genomic variation in Chinook salmon (Oncorhynchus tshawytscha) with emphasis on a region of major effect for ecotypic variation in migration timing. With a filtered data set of ~13 million single nucleotide polymorphisms (SNPs) from low-coverage whole genome resequencing of 53 populations (3566 barcoded individuals), we contrasted patterns of genomic structure within and among major lineages and examined the extent of a selective sweep at a major effect region underlying migration timing (GREB1L/ROCK1). Neutral variation provided support for fine-scale structure of populations, while allele frequency variation in GREB1L/ROCK1 was highly correlated with mean return timing for early and late migrating populations within each of the lineages (r²  = .58-.95; p < .001). However, the extent of selection within the genomic region controlling migration timing was much narrower in one lineage (interior stream-type) compared to the other two major lineages, which corresponded to the breadth of phenotypic variation in migration timing observed among lineages. Evidence of a duplicated block within GREB1L/ROCK1 may be responsible for reduced recombination in this portion of the genome and contributes to phenotypic variation within and across lineages. Lastly, SNP positions across GREB1L/ROCK1 were assessed for their utility in discriminating migration timing among lineages, and we recommend multiple markers nearest the duplication to provide highest accuracy in conservation applications such as those that aim to protect early migrating Chinook salmon. These results highlight the need to investigate variation throughout the genome and the effects of structural variants on ecologically relevant phenotypic variation in natural species.

Going the distance to test local adaptation in Coho salmon

The mechanisms underlying local adaptation, where populations evolve traits that confer advantages to the local environment, is a central topic for understanding evolution in natural systems. Conservation goals for species at risk often include defining population boundaries by identifying gene diversity, genetic differentiation, and adaptation to local environments. In this issue of Molecular Ecology, Rougemont et al. (2022) combine genome-wide SNP data with an extensive set of landscape variables to study the genomic mechanisms of local adaptation in the entire North American range of Coho salmon (Oncorhynchus kisutch), representing one of the largest studies of its kind. Migration distance, defined as the distance adult Coho salmon migrate from the ocean to their freshwater spawning ground, was found to be the primary factor driving local adaptation in this species. With climatic changes altering flow regimes and therefore the success of Coho salmon to return to spawning grounds, understanding environmental drivers and the genomic basis for migration is essential in the conservation of anadromous salmonids.

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.

Implications of Large-Effect Loci for Conservation: A Review and Case Study with Pacific Salmon

The increasing feasibility of assembling large genomic datasets for non-model species presents both opportunities and challenges for applied conservation and management. A popular theme in recent studies is the search for large-effect loci that explain substantial portions of phenotypic variance for a key trait(s). If such loci can be linked to adaptations, 2 important questions arise: 1) Should information from these loci be used to reconfigure conservation units (CUs), even if this conflicts with overall patterns of genetic differentiation? 2) How should this information be used in viability assessments of populations and larger CUs? In this review, we address these questions in the context of recent studies of Chinook salmon and steelhead (anadromous form of rainbow trout) that show strong associations between adult migration timing and specific alleles in one small genomic region. Based on the polygenic paradigm (most traits are controlled by many genes of small effect) and genetic data available at the time showing that early-migrating populations are most closely related to nearby late-migrating populations, adult migration differences in Pacific salmon and steelhead were considered to reflect diversity within CUs rather than separate CUs. Recent data, however, suggest that specific alleles are required for early migration, and that these alleles are lost in populations where conditions do not support early-migrating phenotypes. Contrasting determinations under the US Endangered Species Act and the State of California's equivalent legislation illustrate the complexities of incorporating genomics data into CU configuration decisions. Regardless how CUs are defined, viability assessments should consider that 1) early-migrating phenotypes experience disproportionate risks across large geographic areas, so it becomes important to identify early-migrating populations that can serve as reliable sources for these valuable genetic resources; and 2) genetic architecture, especially the existence of large-effect loci, can affect evolutionary potential and adaptability.

Whole-Genome Resequencing to Evaluate Life History Variation in Anadromous Migration of Oncorhynchus mykiss

Anadromous fish experience physiological modifications necessary to migrate between vastly different freshwater and marine environments, but some species such as Oncorhynchus mykiss demonstrate variation in life history strategies with some individuals remaining exclusively resident in freshwater, whereas others undergo anadromous migration. Because there is limited understanding of genes involved in this life history variation across populations of this species, we evaluated the genomic difference between known anadromous (n = 39) and resident (n = 78) Oncorhynchus mykiss collected from the Klickitat River, WA, USA, with whole-genome resequencing methods. Sequencing of these collections yielded 5.64 million single-nucleotide polymorphisms that were tested for significant differences between resident and anadromous groups along with previously identified candidate gene regions. Although a few regions of the genome were marginally significant, there was one region on chromosome Omy12 that provided the most consistent signal of association with anadromy near two annotated genes in the reference assembly: COP9 signalosome complex subunit 6 (CSN6) and NACHT, LRR, and PYD domain–containing protein 3 (NLRP3). Previously identified candidate genes for anadromy within the inversion region of chromosome Omy05 in coastal steelhead and rainbow trout were not informative for this population as shown in previous studies. Results indicate that the significant region on chromosome Omy12 may represent a minor effect gene for male anadromy and suggests that this life history variation in Oncorhynchus mykiss is more strongly driven by other mechanisms related to environmental rearing such as epigenetic modification, gene expression, and phenotypic plasticity. Further studies into regulatory mechanisms of this trait are needed to understand drivers of anadromy in populations of this protected species.

Applying genomics in assisted migration under climate change: Framework, empirical applications, and case studies

The rate of global climate change is projected to outpace the ability of many natural populations and species to adapt. Assisted migration (AM), which is defined as the managed movement of climate-adapted individuals within or outside the species ranges, is a conservation option to improve species' adaptive capacity and facilitate persistence. Although conservation biologists have long been using genetic tools to increase or maintain diversity of natural populations, genomic techniques could add extra benefit in AM that include selectively neutral and adaptive regions of the genome. In this review, we first propose a framework along with detailed procedures to aid collaboration among scientists, agencies, and local and regional managers during the decision-making process of genomics-guided AM. We then summarize the genomic approaches for applying AM, followed by a literature search of existing incorporation of genomics in AM across taxa. Our literature search initially identified 729 publications, but after filtering returned only 50 empirical studies that were either directly applied or considered genomics in AM related to climate change across taxa of plants, terrestrial animals, and aquatic animals; 42 studies were in plants. This demonstrated limited application of genomic methods in AM in organisms other than plants, so we provide further case studies as two examples to demonstrate the negative impact of climate change on non-model species and how genomics could be applied in AM. With the rapidly developing sequencing technology and accumulating genomic data, we expect to see more successful applications of genomics in AM, and more broadly, in the conservation of biodiversity.

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 region associated with run timing has similar haplotypes and phenotypic effects across three lineages of Chinook salmon

Conserving life-history variation is a stated goal of many management programs, but the most effective means by which to accomplish this are often far from clear. Early- and late-migrating forms of Chinook salmon (Oncorhynchus tshawytscha) face unequal pressure from natural and anthropogenic forces that may alter the impacts of genetic variation underlying heritable migration timing. Genomic regions of chromosome 28 are known to be strongly associated with migration variation in adult Chinook salmon, but it remains unclear whether there is consistent association among diverse lineages and populations in large basins such as the Columbia River. With high-throughput genotyping (GT-seq) and phenotyping methods, we examined the association of genetic variation in 28 markers (spanning GREB1L to ROCK1 of chromosome 28) with individual adult migration timing characteristics gleaned from passive integrated transponder recordings of over 5000 Chinook salmon from the three major phylogeographic lineages that inhabit the Columbia River Basin. Despite the strong genetic differences among them in putatively neutral genomic regions, each of the three lineages exhibited very similar genetic variants in the chromosome 28 region that were significantly associated with adult migration timing phenotypes. This is particularly notable for the interior stream-type lineage, which exhibits an earlier and more constrained freshwater entry than the other lineages. In both interior stream-type and interior ocean-type lineages of Chinook salmon, heterozygotes of the most strongly associated linkage groups had largely intermediate migration timing relative to homozygotes, and results indicate codominance or possibly marginal partial dominance of the allele associated with early migration. Our results lend support to utilization of chromosome 28 variation in tracking and predicting run timing in these lineages of Chinook salmon in the Columbia River.

An evaluation of the potential factors affecting lifetime reproductive success in salmonids

Lifetime reproductive success (LRS), the number of offspring produced over an organism's lifetime, is a fundamental component of Darwinian fitness. For taxa such as salmonids with multiple species of conservation concern, understanding the factors affecting LRS is critical for the development and implementation of successful conservation management practices. Here, we reviewed the published literature to synthesize factors affecting LRS in salmonids including significant effects of hatchery rearing, life history, and phenotypic variation, and behavioral and spawning interactions. Additionally, we found that LRS is affected by competitive behavior on the spawning grounds, genetic compatibility, local adaptation, and hybridization. Our review of existing literature revealed limitations of LRS studies, and we emphasize the following areas that warrant further attention in future research: (1) expanding the range of studies assessing LRS across different life-history strategies, specifically accounting for distinct reproductive and migratory phenotypes; (2) broadening the variety of species represented in salmonid fitness studies; (3) constructing multigenerational pedigrees to track long-term fitness effects; (4) conducting LRS studies that investigate the effects of aquatic stressors, such as anthropogenic effects, pathogens, environmental factors in both freshwater and marine environments, and assessing overall body condition, and (5) utilizing appropriate statistical approaches to determine the factors that explain the greatest variation in fitness and providing information regarding biological significance, power limitations, and potential sources of error in salmonid parentage studies. Overall, this review emphasizes that studies of LRS have profoundly advanced scientific understanding of salmonid fitness, but substantial challenges need to be overcome to assist with long-term recovery of these keystone species in aquatic ecosystems.

The changing face of genome assemblies: Guidance on achieving high-quality reference genomes

The quality of genome assemblies has improved rapidly in recent years due to continual advances in sequencing technology, assembly approaches, and quality control. In the field of molecular ecology, this has led to the development of exceptional quality genome assemblies that will be important long-term resources for broader studies into ecological, conservation, evolutionary, and population genomics of naturally occurring species. Moreover, the extent to which a single reference genome represents the diversity within a species varies: pan-genomes will become increasingly important ecological genomics resources, particularly in systems found to have considerable presence-absence variation in their functional content. Here, we highlight advances in technology that have raised the bar for genome assembly and provide guidance on standards to achieve exceptional quality reference genomes. Key recommendations include the following: (a) Genome assemblies should include long-read sequencing except in rare cases where it is effectively impossible to acquire adequately preserved samples needed for high molecular weight DNA standards. (b) At least one scaffolding approach should be included with genome assembly such as Hi-C or optical mapping. (c) Genome assemblies should be carefully evaluated, this may involve utilising short read data for genome polishing, error correction, k-mer analyses, and estimating the percent of reads that map back to an assembly. Finally, a genome assembly is most valuable if all data and methods are made publicly available and the utility of a genome for further studies is verified through examples. While these recommendations are based on current technology, we anticipate that future advances will push the field further and the molecular ecology community should continue to adopt new approaches that attain the highest quality genome assemblies.

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.

Whole genome resequencing reveals genomic regions associated with thermal adaptation in redband trout

Adaptation to local environments involves evolution of ecologically important traits and underlying physiological processes. Here, we used low coverage whole-genome resequencing (lcWGR) on individuals to identify genome regions involved in thermal adaptation in wild redband trout Oncorhynchus mykiss gairdneri, a subspecies of rainbow trout that inhabits ecosystems ranging from cold montane forests to high elevation deserts. This study includes allele frequency-based analyses for selective sweeps among populations, followed by multiple association tests for specific sets of phenotypes measured under thermal stress (acute and chronic survival/mortality; high or low cardiac performance groups). Depending on the groups in each set of analyses, sequencing reads covered 43%-75% of the genome at ≥15× and each analysis included millions of SNPs across the genome. In tests for selective sweeps among populations, a total of six chromosomal regions were significant. The further association tests for specific phenotypes revealed that the region on chromosome 4 was consistently the most significant and contains the cerk gene (ceramide kinase). This study provides insight into a potential genetic mechanism of local thermal adaptation and suggests cerk may be an important candidate gene. However, further validation of this cerk gene is necessary to determine if the association with cardiac performance results in a functional role to influence thermal performance when exposed to high water temperatures and hypoxic conditions.

Validation and association of candidate markers for adult migration timing and fitness in Chinook Salmon

Recent studies have begun to elucidate the genetic basis for phenotypic traits in salmonid species, but many questions remain before these candidate genes can be directly incorporated into conservation management. In Chinook Salmon (Oncorhynchus tshawytscha), a region of major effect for migration timing has been discovered that harbors two adjacent candidate genes (greb1L, rock1), but there has been limited work to examine the association between these genes and migratory phenotypes at the individual, compared to the population, level. To provide a more thorough test of individual phenotypic association within lineages of Chinook Salmon, 33 candidate markers were developed across a 220 Kb region on chromosome 28 previously associated with migration timing. Candidate and neutral markers were genotyped in individuals from representative collections that exhibit phenotypic variation in timing of arrival to spawning grounds from each of three lineages of Chinook Salmon. Association tests confirmed the majority of markers on chromosome 28 were significantly associated with arrival timing and the strongest association was consistently observed for markers within the rock1 gene and the intergenic region between greb1L and rock1. Candidate markers alone explained a wide range of phenotypic variation for Lower Columbia and Interior ocean-type lineages (29% and 78%, respectively), but less for the Interior stream-type lineage (5%). Individuals that were heterozygous at markers within or upstream of rock1 had phenotypes that suggested a pattern of dominant inheritance for early arrival across populations. Finally, previously published fitness estimates from the Interior stream-type lineage enabled tests of association with arrival timing and two candidate markers, which revealed that fish with homozygous mature genotypes had slightly higher fitness than fish with premature genotypes, while heterozygous fish were intermediate. Overall, these results provide additional information for individual-level genetic variation associated with arrival timing that may assist with conservation management of this species.

Genomic islands of divergence infer a phenotypic landscape in Pacific lamprey

High rates of dispersal can breakdown coadapted gene complexes. However, concentrated genomic architecture (i.e., genomic islands of divergence) can suppress recombination to allow evolution of local adaptations despite high gene flow. Pacific lamprey (Entosphenus tridentatus) is a highly dispersive anadromous fish. Observed trait diversity and evidence for genetic basis of traits suggests it may be locally adapted. We addressed whether concentrated genomic architecture could influence local adaptation for Pacific lamprey. Using two new whole genome assemblies and genotypes from 7,716 single nucleotide polymorphism (SNP) loci in 518 individuals from across the species range, we identified four genomic islands of divergence (on chromosomes 01, 02, 04, and 22). We determined robust phenotype-by-genotype relationships by testing multiple traits across geographic sites. These trait associations probably explain genomic divergence across the species' range. We genotyped a subset of 302 broadly distributed SNPs in 2,145 individuals for association testing for adult body size, sexual maturity, migration distance and timing, adult swimming ability, and larval growth. Body size traits were strongly associated with SNPs on chromosomes 02 and 04. Moderate associations also implicated SNPs on chromosome 01 as being associated with variation in female maturity. Finally, we used candidate SNPs to extrapolate a heterogeneous spatiotemporal distribution of these predicted phenotypes based on independent data sets of larval and adult collections. These maturity and body size results guide future elucidation of factors driving regional optimization of these traits for fitness. Pacific lamprey is culturally important and imperiled. This research addresses biological uncertainties that challenge restoration efforts.

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.

Long-term evaluation of fitness and demographic effects of a Chinook Salmon supplementation program

While the goal of supplementation programs is to provide positive, population-level effects for species of conservation concern, these programs can also present an inherent fitness risk when captive-born individuals are fully integrated into the natural population. In order to evaluate the long-term effects of a supplementation program and estimate the demographic and phenotypic factors influencing the fitness of a threatened population of Chinook Salmon (Oncorhynchus tshawytscha), we genotyped tissue samples spanning a 19-year period (1998-2016) to generate pedigrees from adult fish returning to Johnson Creek, Idaho, USA. We expanded upon previous estimates of relative reproductive success (RRS) to include grandparentage analyses and used generalized linear models to determine whether origin (hatchery or natural) or phenotypic traits (timing of arrival to spawning grounds, body length, and age) significantly predicted reproductive success (RS) across multiple years. Our results provide evidence that this supplementation program with 100% natural-origin broodstock provided a long-term demographic boost to the population (mean of 4.56 times in the first generation and mean of 2.52 times in the second generation). Overall, when spawning in nature, hatchery-origin fish demonstrated a trend toward lower RS compared to natural-origin fish (p < 0.05). However, when hatchery-origin fish successfully spawned with natural-origin fish, they had similar RS compared to natural by natural crosses (first-generation mean hatchery by natural cross RRS = 1.11 females, 1.13 males; second-generation mean hatchery by natural cross RRS = 1.03 females, 1.08 males). While origin, return year, and body length were significant predictors of fitness for both males and females (p < 0.05), return day was significant for males but not females (p > 0.05). These results indicate that supplementation programs that reduce the potential for genetic adaptation to captivity can be effective at increasing population abundance while limiting long-term fitness effects on wild populations.

Physiological and genomic signatures of evolutionary thermal adaptation in redband trout from extreme climates

Temperature is a master environmental factor that limits the geographical distribution of species, especially in ectotherms. To address challenges in biodiversity conservation under ongoing climate change, it is essential to characterize relevant functional limitations and adaptive genomic content at population and species levels. Here, we present evidence for adaptive divergence in cardiac function and genomic regions in redband trout (Oncorhynchus mykiss gairdneri) populations from desert and montane streams. Cardiac phenotypes of individual fish were measured in the field with a custom-built electrocardiogram apparatus. Maximum heart rate and its rate limiting temperature during acute warming were significantly higher in fish that have evolved in the extreme of a desert climate compared to a montane climate. Association mapping with 526,301 single nucleotide polymorphisms (SNPs) across the genome revealed signatures of thermal selection both within and among ecotypes. Among desert and montane populations, 435 SNPs were identified as putative outliers under natural selection and 20 of these loci showed significant association with average summer water temperatures among populations. Phenotypes for cardiac performance were variable within each ecotype, and 207 genomic regions were strongly associated with either maximum heart rate or rate limiting temperatures among individuals. Annotation of significant loci provided candidate genes that underlie thermal adaptation, including pathways associated with cardiac function (IRX5, CASQ1, CAC1D, and TITIN), neuroendocrine system (GPR17 and NOS), and stress response (SERPH). By integrating comparative physiology and population genomics, results here advance our knowledge on evolutionary processes of thermal adaptation in aquatic ectotherms.

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.

Genomic variation underlying complex life-history traits revealed by genome sequencing in Chinook salmon

A broad portfolio of phenotypic diversity in natural organisms can buffer against exploitation and increase species persistence in disturbed ecosystems. The study of genomic variation that accounts for ecological and evolutionary adaptation can represent a powerful approach to extend understanding of phenotypic variation in nature. Here we present a chromosome-level reference genome assembly for Chinook salmon (Oncorhynchus tshawytscha; 2.36 Gb) that enabled association mapping of life-history variation and phenotypic traits for this species. Whole-genome re-sequencing of populations with distinct life-history traits provided evidence that divergent selection was extensive throughout the genome within and among phylogenetic lineages, indicating that a broad portfolio of phenotypic diversity exists in this species that is related to local adaptation and life-history variation. Association mapping with millions of genome-wide SNPs revealed that a genomic region of major effect on chromosome 28 was associated with phenotypes for premature and mature arrival to spawning grounds and was consistent across three distinct phylogenetic lineages. Our results demonstrate how genomic resources can enlighten the genetic basis of known phenotypes in exploited species and assist in clarifying phenotypic variation that may be difficult to observe in naturally occurring organisms.

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.

Mechanisms of thermal adaptation and evolutionary potential of conspecific populations to changing environments

Heterogeneous and ever-changing thermal environments drive the evolution of populations and species, especially when extreme conditions increase selection pressure for traits influencing fitness. However, projections of biological diversity under scenarios of climate change rarely consider evolutionary adaptive potential of natural species. In this study, we tested for mechanistic evidence of evolutionary thermal adaptation among ecologically divergent redband trout populations (Oncorhynchus mykiss gairdneri) in cardiorespiratory function, cellular response and genomic variation. In a common garden environment, fish from an extreme desert climate had significantly higher critical thermal maximum (p < .05) and broader optimum thermal window for aerobic scope (>3°C) than fish from cooler montane climate. In addition, the desert population had the highest maximum heart rate during warming (20% greater than montane populations), indicating improved capacity to deliver oxygen to internal tissues. In response to acute heat stress, distinct sets of cardiac genes were induced among ecotypes, which helps to explain the differences in cardiorespiratory function. Candidate genomic markers and genes underlying these physiological adaptations were also pinpointed, such as genes involved in stress response and metabolic activity (hsp40, ldh-b and camkk2). These markers were developed into a multivariate model that not only accurately predicted critical thermal maxima, but also evolutionary limit of thermal adaptation in these specific redband trout populations relative to the expected limit for the species. This study demonstrates mechanisms and limitations of an aquatic species to evolve under changing environments that can be incorporated into advanced models to predict ecological consequences of climate change for natural organisms.

Vive la résistance: genome-wide selection against introduced alleles in invasive hybrid zones

Evolutionary and ecological consequences of hybridization between native and invasive species are notoriously complicated because patterns of selection acting on non-native alleles can vary throughout the genome and across environments. Rapid advances in genomics now make it feasible to assess locus-specific and genome-wide patterns of natural selection acting on invasive introgression within and among natural populations occupying diverse environments. We quantified genome-wide patterns of admixture across multiple independent hybrid zones of native westslope cutthroat trout and invasive rainbow trout, the world's most widely introduced fish, by genotyping 339 individuals from 21 populations using 9380 species-diagnostic loci. A significantly greater proportion of the genome appeared to be under selection favouring native cutthroat trout (rather than rainbow trout), and this pattern was pervasive across the genome (detected on most chromosomes). Furthermore, selection against invasive alleles was consistent across populations and environments, even in those where rainbow trout were predicted to have a selective advantage (warm environments). These data corroborate field studies showing that hybrids between these species have lower fitness than the native taxa, and show that these fitness differences are due to selection favouring many native genes distributed widely throughout the genome.

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.

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.

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 (F_ST = 0.004), and evidence that immigrant influences were prevalent and ubiquitous. Populations in the upper watershed above partial natural barriers were highly distinct (F_ST = 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.

Genetic basis of adult migration timing in anadromous steelhead discovered through multivariate association testing

Migration traits are presumed to be complex and to involve interaction among multiple genes. We used both univariate analyses and a multivariate random forest (RF) machine learning algorithm to conduct association mapping of 15 239 single nucleotide polymorphisms (SNPs) for adult migration-timing phenotype in steelhead (Oncorhynchus mykiss). Our study focused on a model natural population of steelhead that exhibits two distinct migration-timing life histories with high levels of admixture in nature. Neutral divergence was limited between fish exhibiting summer- and winter-run migration owing to high levels of interbreeding, but a univariate mixed linear model found three SNPs from a major effect gene to be significantly associated with migration timing (p < 0.000005) that explained 46% of trait variation. Alignment to the annotated Salmo salar genome provided evidence that all three SNPs localize within a 46 kb region overlapping GREB1-like (an oestrogen target gene) on chromosome Ssa03. Additionally, multivariate analyses with RF identified that these three SNPs plus 15 additional SNPs explained up to 60% of trait variation. These candidate SNPs may provide the ability to predict adult migration timing of steelhead to facilitate conservation management of this species, and this study demonstrates the benefit of multivariate analyses for association studies.

Genomic patterns of diversity and divergence of two introduced salmonid species in Patagonia, South America

Invasive species have become widespread in aquatic environments throughout the world, yet there are few studies that have examined genomic variation of multiple introduced species in newly colonized environments. In this study, we contrast genomic variation in two salmonid species (anadromous Chinook Salmon, Oncorhynchus tshawytscha, 11,579 SNPs and resident Brook Charr Salvelinus fontinalis, 13,522 SNPs) with differing invasion success after introduction to new environments in South America relative to populations from their native range in North America. Estimates of genetic diversity were not significantly different between introduced and source populations for either species, indicative of propagule pressure that has been shown to maintain diversity in founding populations relative to their native range. Introduced populations also demonstrated higher connectivity and gene flow than those in their native range. Evidence for candidate loci under divergent selection was observed, but was limited to specific introduced populations and was not widely evident. Patterns of genomic variation were consistent with general dispersal potential of each species and therefore also the notion that life history variation may contribute to both invasion success and subsequent genetic structure of these two salmonids in Patagonia.

Genomic signatures among Oncorhynchus nerka ecotypes to inform conservation and management of endangered Sockeye Salmon

Conservation of life history variation is an important consideration for many species with trade-offs in migratory characteristics. Many salmonid species exhibit both resident and migratory strategies that capitalize on benefits in freshwater and marine environments. In this study, we investigated genomic signatures for migratory life history in collections of resident and anadromous Oncorhynchus nerka (Kokanee and Sockeye Salmon, respectively) from two lake systems, using ~2,600 SNPs from restriction-site-associated DNA sequencing (RAD-seq). Differing demographic histories were evident in the two systems where one pair was significantly differentiated (Redfish Lake, F_ST = 0.091 [95% confidence interval: 0.087 to 0.095]) but the other pair was not (Alturas Lake, F_ST = -0.007 [-0.008 to -0.006]). Outlier and association analyses identified several candidate markers in each population pair, but there was limited evidence for parallel signatures of genomic variation associated with migration. Despite lack of evidence for consistent markers associated with migratory life history in this species, candidate markers were mapped to functional genes and provide evidence for adaptive genetic variation within each lake system. Life history variation has been maintained in these nearly extirpated populations of O. nerka, and conservation efforts to preserve this diversity are important for long-term resiliency of this species.

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.

Differential Expression of Genes that Control Respiration Contribute to Thermal Adaptation in Redband Trout (Oncorhynchus mykiss gairdneri)

Organisms can adapt to local environmental conditions as a plastic response or become adapted through natural selection on genetic variation. The ability to adapt to increased water temperatures will be of paramount importance for many fish species as the climate continues to warm and water resources become limited. Because increased water temperatures will reduce the dissolved oxygen available for fish, we hypothesized that adaptation to low oxygen environments would involve improved respiration through oxidative phosphorylation (OXPHOS). To test this hypothesis, we subjected individuals from two ecologically divergent populations of inland (redband) rainbow trout (Oncorhynchus mykiss gairdneri) with historically different temperature regimes (desert and montane) and their F1 progeny to diel cycles of temperature stress and then examined gene expression data for 80 nuclear- and mitochondrial-encoded OXPHOS subunits that participate in respiration. Of the 80 transcripts, 7 showed ≥ 2-fold difference in expression levels in gill tissue from desert fish under heat stress whereas the montane fish had none and the F1 only had one differentially expressed gene. A structural analysis of the proteins encoded by those genes suggests that the response could coordinate the formation of supercomplexes and oligomers. Supercomplexes may increase the efficiency of respiration because complexes I, III, and IV are brought into close proximity and oligomerization of complex V alters the macrostructure of mitochondria to improve respiration. Significant differences in gene expression patterns in response to heat stress in a common environment indicate that the response was not due to plasticity but had a genetic basis.

Transcriptomic response to heat stress among ecologically divergent populations of redband trout

Background

As ectothermic organisms have evolved to differing aquatic climates, the molecular basis of thermal adaptation is a key area of research. In this study, we tested for differential transcriptional response of ecologically divergent populations of redband trout (Oncorhynchus mykiss gairdneri) that have evolved in desert and montane climates. Each pure strain and their F1 cross were reared in a common garden environment and exposed over four weeks to diel water temperatures that were similar to those experienced in desert climates within the species’ range. Gill tissues were collected from the three strains of fish (desert, montane, F1 crosses) at the peak of heat stress and tested for mRNA expression differences across the transcriptome with RNA-seq.

Results

Strong differences in transcriptomic response to heat stress were observed across strains confirming that fish from desert environments have evolved diverse mechanisms to cope with stressful environments. As expected, a large number of total transcripts (12,814) were differentially expressed in the study (FDR ≤ 0.05) with 2310 transcripts in common for all three strains, but the desert strain had a larger number of unique differentially expressed transcripts (2875) than the montane (1982) or the F1 (2355) strain. Strongly differentiated genes (>4 fold change and FDR ≤ 0.05) were particularly abundant in the desert strain (824 unique contigs) relative to the other two strains (montane = 58; F1 = 192).

Conclusions

This study demonstrated patterns of acclimation (i.e., phenotypic plasticity) within strains and evolutionary adaptation among strains in numerous genes throughout the transcriptome. Key stress response genes such as molecular chaperones (i.e., heat shock proteins) had adaptive patterns of gene expression among strains, but also a much higher number of metabolic and cellular process genes were differentially expressed in the desert strain demonstrating these biological pathways are critical for thermal adaptation to warm aquatic climates. The results of this study further elucidate the molecular basis for thermal adaptation in aquatic ecosystems and extend the potential for identifying genes that may be critical for adaptation to changing climates.

Electronic supplementary material

The online version of this article (doi:10.1186/s12864-015-1246-5) contains supplementary material, which is available to authorized users.

Selection for upper thermal tolerance in rainbow trout (Oncorhynchus mykiss Walbaum)

Rainbow trout (Oncorhynchus mykiss Walbaum) in southern Western Australia have undergone passive selection for over 19 generations to survive high water temperatures. Based on the conceptual model of 'oxygen- and capacity-limited thermal tolerance', we measured critical thermal maximum (CTmax), maximum heart rate (fH,max) and aerobic scope to test the hypothesis that these rainbow trout can maintain aerobic scope at high temperatures through a robust cardiac performance supporting oxygen delivery. Across five family groups CTmax averaged 29.0±0.02°C. Aerobic scope was maximized at 15.8±0.3°C (Topt), while the upper pejus temperature (Tpej, set at 90% of maximum aerobic scope) was 19.9±0.3°C. Although aerobic scope decreased at temperatures above Topt, the value at 25°C remained well over 40% of the maximum. Furthermore, pharmacologically stimulated fH,max increased with temperature, reaching a peak value between 23.5±0.4 and 24.0±0.4°C (Tmax) for three family groups. The Arrhenius breakpoint temperature (TAB) for fH,max was 20.3±0.3 to 20.7±0.4°C, while the average Q10 breakpoint temperature (TQB, when the incremental Q10<1.6) for fH,max was 21.6±0.2 to 22.0±0.4°C. Collectively, fH,max progressively became less temperature dependent beyond 20°C (TAB and TQB), which coincides with the upper Tpej for aerobic scope. Although upper thermal performance indices for both aerobic scope and fH,max were compared among family groups in this population, appreciable differences were not evident. Compared with other populations of rainbow trout, the present assessment is consistent with the prediction that this strain has undergone selection and shows the ability to tolerate higher water temperatures.

Genes predict long distance migration and large body size in a migratory fish, Pacific lamprey

Elucidation of genetic mechanisms underpinning migratory behavior could help predict how changes in genetic diversity may affect future spatiotemporal distribution of a migratory species. This ability would benefit conservation of one such declining species, anadromous Pacific lamprey (Entosphenus tridentatus). Nonphilopatric migration of adult Pacific lamprey has homogenized population-level neutral variation but has maintained adaptive variation that differentiates groups based on geography, run-timing and adult body form. To investigate causes for this adaptive divergence, we examined 647 adult lamprey sampled at a fixed location on the Columbia River and radiotracked during their subsequent upstream migration. We tested whether genetic variation [94 neutral and adaptive single nucleotide polymorphisms (SNPs) previously identified from a genomewide association study] was associated with phenotypes of migration distance, migration timing, or morphology. Three adaptive markers were strongly associated with morphology, and one marker also correlated with upstream migration distance and timing. Genes physically linked with these markers plausibly influence differences in body size, which is also consistently associated with migration distance in Pacific lamprey. Pacific lamprey conservation implications include the potential to predict an individual's upstream destination based on its genotype. More broadly, the results suggest a genetic basis for intrapopulation variation in migration distance in migratory species.

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.

Use of genotyping by sequencing data to develop a high-throughput and multifunctional SNP panel for conservation applications in Pacific lamprey

Next-generation sequencing data can be mined for highly informative single nucleotide polymorphisms (SNPs) to develop high-throughput genomic assays for nonmodel organisms. However, choosing a set of SNPs to address a variety of objectives can be difficult because SNPs are often not equally informative. We developed an optimal combination of 96 high-throughput SNP assays from a total of 4439 SNPs identified in a previous study of Pacific lamprey (Entosphenus tridentatus) and used them to address four disparate objectives: parentage analysis, species identification and characterization of neutral and adaptive variation. Nine of these SNPs are FST outliers, and five of these outliers are localized within genes and significantly associated with geography, run-timing and dwarf life history. Two of the 96 SNPs were diagnostic for two other lamprey species that were morphologically indistinguishable at early larval stages and were sympatric in the Pacific Northwest. The majority (85) of SNPs in the panel were highly informative for parentage analysis, that is, putatively neutral with high minor allele frequency across the species' range. Results from three case studies are presented to demonstrate the broad utility of this panel of SNP markers in this species. As Pacific lamprey populations are undergoing rapid decline, these SNPs provide an important resource to address critical uncertainties associated with the conservation and recovery of this imperiled species.