Genomewide association analyses of fitness traits in captive-reared Chinook salmon: Applications in evaluating conservation strategies

Salmon hatchery strays can demographically boost wild populations at the cost of diversity: quantitative genetic modelling of Alaska pink salmon

Hatcheries are vital to many salmon fisheries, with inherent risks and rewards. While hatcheries can increase the returns of adult fish, the demographic and evolutionary consequences for natural populations interacting with hatchery fish on spawning grounds remain unclear. This study examined the impacts of stray hatchery-origin pink salmon on natural population productivity and resilience. We explored temporal assortative mating dynamics using a quantitative genetic model that assumed the only difference between hatchery- and natural-origin adults was their return timing to natural spawning grounds. This model was parameterized with empirical data from an intensive multi-generational study of hatchery-wild interactions in the world's largest pink salmon fisheries enhancement program located in Prince William Sound, Alaska. Across scenarios of increasing hatchery fish presence on spawning grounds, our findings underscore a trade-off between demographic enhancement and preservation of natural population diversity. While enhancement bolstered natural population sizes towards local carrying capacities, hatchery introgression reduced variation in adult return timing by up to 20%. Results indicated that hatchery-origin alleles can rapidly assimilate into natural populations, despite the reduced fitness of hatchery fish attributable to phenotypic mismatches. These findings elucidate the potential for long-term demographic and evolutionary consequences arising from specific hatchery-wild interactions, emphasizing the need for management strategies that balance demographic enhancement with the conservation of natural diversity.

A single generation in the wild increases fitness for descendants of hatchery-origin Chinook salmon (Oncorhynchus tshawytscha)

Reintroduction is an important tool for the recovery of imperiled species. For threatened Pacific salmonids (Oncorhynchus spp.) species, hatchery-origin (HOR) individuals from a nearby source are often used to reestablish populations in vacant, historically occupied habitat. However, this approach is challenged by the relatively low reproductive success that HOR Pacific salmonids experience when they spawn in the wild, relative to their natural-origin (NOR) counterparts. In this study, we used genetic parentage analysis to compare the reproductive success of three groups of adult Chinook salmon (Oncorhynchus tshawytscha) reintroduced above Cougar Dam on the South Fork McKenzie River, Oregon: HOR Chinook salmon from an integrated stock; first-generation, wild-born descendants (hereafter F 1s) of Chinook salmon produced at the same hatchery; and NOR Chinook salmon that are presumed to have been produced below the dam, on the mainstem McKenzie River, or elsewhere and volitionally entered a trap below Cougar Dam. We found that F 1s produced nearly as many adult offspring as NORs, and 1.8-fold more adult offspring than HORs. This result suggests that, for the South Fork McKenzie reintroduction program, a single generation in the wild increases fitness for the descendants of HOR Chinook salmon. Although these results are encouraging, care must be taken before extrapolating our results to other systems.

Genomic evidence for domestication selection in three hatchery populations of Chinook salmon, Oncorhynchus tshawytscha

Fish hatcheries are widely used to enhance fisheries and supplement declining wild populations. However, substantial evidence suggests that hatchery fish are subject to differential selection pressures compared to their wild counterparts. Domestication selection, or adaptation to the hatchery environment, poses a risk to wild populations if traits specific to success in the hatchery environment have a genetic component and there is subsequent introgression between hatchery and wild fish. Few studies have investigated domestication selection in hatcheries on a genomic level, and even fewer have done so in parallel across multiple hatchery-wild population pairs. In this study, we used low-coverage whole-genome sequencing to investigate signals of domestication selection in three separate hatchery populations of Chinook salmon, Oncorhynchus tshawytscha, after approximately seven generations of divergence from their corresponding wild progenitor populations. We sequenced 192 individuals from populations across Southeast Alaska and estimated genotype likelihoods at over six million loci. We discovered a total of 14 outlier peaks displaying high genetic differentiation (F ST) between hatchery-wild pairs, although no peaks were shared across the three comparisons. Peaks were small (53 kb on average) and often displayed elevated absolute genetic divergence (D xy) and linkage disequilibrium, suggesting some level of domestication selection has occurred. Our study provides evidence that domestication selection can lead to genetic differences between hatchery and wild populations in only a few generations. Additionally, our data suggest that population-specific adaptation to hatchery environments likely occurs through different genetic pathways, even for populations with similar standing genetic variation. These results highlight the need to collect paired genotype-phenotype data to understand how domestication may be affecting fitness and to identify potential management practices that may mitigate genetic risks despite multiple pathways of domestication.

Understanding Local Adaptation to Prepare Populations for Climate Change

Adaptation within species to local environments is widespread in nature. Better understanding this local adaptation is critical to conserving biodiversity. However, conservation practices can rely on species’ trait averages or can broadly assume homogeneity across the range to inform management. Recent methodological advances for studying local adaptation provide the opportunity to fine-tune efforts for managing and conserving species. The implementation of these advances will allow us to better identify populations at greatest risk of decline because of climate change, as well as highlighting possible strategies for improving the likelihood of population persistence amid climate change. In the present article, we review recent advances in the study of local adaptation and highlight ways these tools can be applied in conservation efforts. Cutting-edge tools are available to help better identify and characterize local adaptation. Indeed, increased incorporation of local adaptation in management decisions may help meet the imminent demands of managing species amid a rapidly changing world.

Genomic differentiation in Pacific cod using Pool-Seq

Patterns of genetic differentiation across the genome can provide insight into selective forces driving adaptation. We used pooled whole genome sequencing, gene annotation, and environmental covariates to evaluate patterns of genomic differentiation and to investigate mechanisms responsible for divergence among proximate Pacific cod (Gadus macrocephalus) populations from the Bering Sea and Aleutian Islands and more distant Washington Coast cod. Samples were taken from eight spawning locations, three of which were replicated to estimate consistency in allele frequency estimation. A kernel smoothing moving weighted average of relative divergence (F ST) identified 11 genomic islands of differentiation between the Aleutian Islands and Bering Sea samples. In some islands of differentiation, there was also elevated absolute divergence (d XY) and evidence for selection, despite proximity and potential for gene flow. Similar levels of absolute divergence (d XY) but roughly double the relative divergence (F ST) were observed between the distant Bering Sea and Washington Coast samples. Islands of differentiation were much smaller than the four large inversions among Atlantic cod ecotypes. Islands of differentiation between the Bering Sea and Aleutian Island were associated with SNPs from five vision system genes, which can be associated with feeding, predator avoidance, orientation, and socialization. We hypothesize that islands of differentiation between Pacific cod from the Bering Sea and Aleutian Islands provide evidence for adaptive differentiation despite gene flow in this commercially important marine species.

Gene Pool Boundaries for the Yosemite Toad (Anaxyrus canorus) Reveal Asymmetrical Migration Within Meadow Neighborhoods

The Yosemite toad (Anaxyrus [Bufo] canorus) is a federally threatened species of meadow-specializing amphibian endemic to the high-elevation Sierra Nevada Mountains of California. The species is one of the first amphibians to undergo a large demographic collapse that was well-documented, and is reputed to remain in low abundance throughout its range. Recent phylogeographic work has demonstrated that Pleistocene toad lineages diverged and then admixed to differing extents across an elevational gradient. Although lineage divisions may have significant effects on evolutionary trajectories over large spatial and temporal scales, present-day population dynamics must be delineated in order to manage and conserve the species effectively. In this study, we used a double-digest RADseq dataset to address three primary questions: (1) Are single meadows or neighborhoods of nearby meadows most correlated with population boundaries? (2) Does asymmetrical migration occur among neighborhoods of nearby meadows? (3) What topographic or hydrological variables predict such asymmetrical migration in these meadow neighborhoods? Hierarchical STRUCTURE and AMOVA analyses suggested that populations are typically circumscribed by a single meadow, although 84% of meadows exist in neighborhoods of at least two meadows connected by low levels of migration, and over half (53%) of neighborhoods examined display strong asymmetrical migration. Meadow neighborhoods often contain one or more large and flat “hub” meadows that experience net immigration, surrounded by smaller and topographically rugged “satellite” meadows with net emigration. Hubs tend to contain more genetic diversity and could be prioritized for conservation and habitat management and as potential sources for reestablishment efforts.

IIb-RAD-sequencing coupled with random forest classification indicates regional population structuring and sex-specific differentiation in salmon lice (Lepeophtheirus salmonis)

The aquaculture industry has been dealing with salmon lice problems forming serious threats to salmonid farming. Several treatment approaches have been used to control the parasite. Treatment effectiveness must be optimized, and the systematic genetic differences between subpopulations must be studied to monitor louse species and enhance targeted control measures. We have used IIb-RAD sequencing in tandem with a random forest classification algorithm to detect the regional genetic structure of the Norwegian salmon lice and identify important markers for sex differentiation of this species. We identified 19,428 single nucleotide polymorphisms (SNPs) from 95 individuals of salmon lice. These SNPs, however, were not able to distinguish the differential structure of lice populations. Using the random forest algorithm, we selected 91 SNPs important for geographical classification and 14 SNPs important for sex classification. The geographically important SNP data substantially improved the genetic understanding of the population structure and classified regional demographic clusters along the Norwegian coast. We also uncovered SNP markers that could help determine the sex of the salmon louse. A large portion of the SNPs identified to be under directional selection was also ranked highly important by random forest. According to our findings, there is a regional population structure of salmon lice associated with the geographical location along the Norwegian coastline.

Heterogeneous genetic basis of age at maturity in salmonid fishes

Understanding the genetic basis of repeated evolution of the same phenotype across taxa is a fundamental aim in evolutionary biology and has applications in conservation and management. However, the extent to which interspecific life-history trait polymorphisms share evolutionary pathways remains underexplored. Here, we address this gap by studying the genetic basis of a key life-history trait, age at maturity, in four species of Pacific salmonids (genus Oncorhynchus) that exhibit intra- and interspecific variation in this trait-Chinook Salmon, Coho Salmon, Sockeye Salmon, and Steelhead Trout. We tested for associations in all four species between age at maturity and two genome regions, six6 and vgll3, that are strongly associated with the same trait in Atlantic Salmon (Salmo salar). We also conducted a genome-wide association analysis in Steelhead to assess whether additional regions were associated with this trait. We found the genetic basis of age at maturity to be heterogeneous across salmonid species. Significant associations between six6 and age at maturity were observed in two of the four species, Sockeye and Steelhead, with the association in Steelhead being particularly strong in both sexes (p = 4.46 × 10^-9 after adjusting for genomic inflation). However, no significant associations were detected between age at maturity and the vgll3 genome region in any of the species, despite its strong association with the same trait in Atlantic Salmon. We discuss possible explanations for the heterogeneous nature of the genetic architecture of this key life-history trait, as well as the implications of our findings for conservation and management.

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.

Y-chromosome haplotypes are associated with variation in size and age at maturity in male Chinook salmon

Variation in size and age at maturity is an important component of life history that is influenced by both environmental and genetic factors. In salmonids, large size confers a direct reproductive advantage through increased fecundity and egg quality in females, while larger males gain a reproductive advantage by monopolizing access to females. In addition, variation in size and age at maturity in males can be associated with different reproductive strategies; younger smaller males may gain reproductive success by sneaking among mating pairs. In both sexes, there is a trade-off between older age and increased reproductive success and increased risk of mortality by delaying reproduction. We identified four Y-chromosome haplogroups that showed regional- and population-specific variation in frequency using RADseq data for 21 populations of Alaska Chinook salmon. We then characterized the range-wide distribution of these haplogroups using GT-seq assays. These haplogroups exhibited associations with size at maturity in multiple populations, suggesting that lack of recombination between X and Y-chromosomes has allowed Y-chromosome haplogroups to capture different alleles that influence size at maturity. Ultimately, conservation of life history diversity in Chinook salmon may require conservation of Y-chromosome haplotype diversity.

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

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

Power of a dual-use SNP panel for pedigree reconstruction and population assignment

The use of high-throughput, low-density sequencing approaches has dramatically increased in recent years in studies of eco-evolutionary processes in wild populations and domestication in commercial aquaculture. Most of these studies focus on identifying panels of SNP loci for a single downstream application, whereas there have been few studies examining the trade-offs for selecting panels of markers for use in multiple applications. Here, we detail the use of a bioinformatic workflow for the development of a dual-purpose SNP panel for parentage and population assignment, which included identifying putative SNP loci, filtering for the most informative loci for the two tasks, designing effective multiplex PCR primers, optimizing the SNP panel for performance, and performing quality control steps for downstream applications. We applied this workflow to two adjacent Alaskan Sockeye Salmon populations and identified a GTseq panel of 142 SNP loci for parentage and 35 SNP loci for population assignment. Only 50-75 panel loci were necessary for >95% accurate parentage, whereas population assignment success, with all 172 panel loci, ranged from 93.9% to 96.2%. Finally, we discuss the trade-offs and complexities of the decision-making process that drives SNP panel development, optimization, and testing.

Multiple decades of stocking has resulted in limited hatchery introgression in wild brook trout (Salvelinus fontinalis) populations of Nova Scotia

Many populations of freshwater fishes are threatened with losses, and increasingly, the release of hatchery individuals is one strategy being implemented to support wild populations. However, stocking of hatchery individuals may pose long-term threats to wild populations, particularly if genetic interactions occur between wild and hatchery individuals. One highly prized sport fish that has been heavily stocked throughout its range is the brook trout (Salvelinus fontinalis). In Nova Scotia, Canada, hatchery brook trout have been stocked since the early 1900s, and despite continued stocking efforts, populations have suffered declines in recent decades. Before this study, the genetic structure of brook trout populations in the province was unknown; however, given the potential negative consequences associated with hatchery stocking, it is possible that hatchery programs have adversely affected the genetic integrity of wild populations. To assess the influence of hatchery supplementation on wild populations, we genotyped wild brook trout from 12 river systems and hatchery brook trout from two major hatcheries using 100 microsatellite loci. Genetic analyses of wild trout revealed extensive population genetic structure among and within river systems and significant isolation-by-distance. Hatchery stocks were genetically distinct from wild populations, and most populations showed limited to no evidence of hatchery introgression (<5% hatchery ancestry). Only a single location had a substantial number of hatchery-derived trout and was located in the only river where a local strain is used for supplementation. The amount of hatchery stocking within a watershed did not influence the level of hatchery introgression. Neutral genetic structure of wild populations was influenced by geography with some influence of climate and stocking indices. Overall, our study suggests that long-term stocking has not significantly affected the genetic integrity of wild trout populations, highlighting the variable outcomes of stocking and the need to evaluate the consequences on a case-by-case basis.

Genomic and phenotypic effects of inbreeding across two different hatchery management regimes in Chinook salmon

Genomic approaches permit direct estimation of inbreeding and its effect on fitness. We used genomic-based estimates of inbreeding to investigate their relationship with eight adult traits in a captive-reared Pacific salmonid that is released into the wild. Estimates were also used to determine whether alternative broodstock management approaches reduced risks of inbreeding. Specifically, 1,100 unlinked restriction-site associated (RAD) loci were used to compare pairwise relatedness, derived from a relationship matrix, and individual inbreeding, estimated by comparing observed and expected homozygosity, across four generations in two hatchery lines of Chinook salmon that were derived from the same source. The lines are managed as "integrated" with the founding wild stock, with ongoing gene flow, and as "segregated" with no gene flow. While relatedness and inbreeding increased in the first generation of both lines, possibly due to population subdivision caused by hatchery initiation, the integrated line had significantly lower levels in some subsequent generations (relatedness: F₂ -F₄ ; inbreeding F₂ ). Generally, inbreeding was similar between the lines despite large differences in effective numbers of breeders. Inbreeding did not affect fecundity, reproductive effort, return timing, fork length, weight, condition factor, and daily growth coefficient. However, it delayed spawn timing by 1.75 days per one standard deviation increase in F (~0.16). The results indicate that integrated management may reduce inbreeding but also suggest that it is relatively low in a small, segregated hatchery population that maximized number of breeders. Our findings demonstrate the utility of genomics to monitor inbreeding under alternative management strategies in captive breeding programs.

Conservation and Management of Salmon in the Age of Genomics

Salmon were among the first nonmodel species for which systematic population genetic studies of natural populations were conducted, often to support management and conservation. The genomics revolution has improved our understanding of the evolutionary ecology of salmon in two major ways: (a) Large increases in the numbers of genetic markers (from dozens to 10⁴-10⁶) provide greater power for traditional analyses, such as the delineation of population structure, hybridization, and population assignment, and (b) qualitatively new insights that were not possible with traditional genetic methods can be achieved by leveraging detailed information about the structure and function of the genome. Studies of the first type have been more common to date, largely because it has taken time for the necessary tools to be developed to fully understand the complex salmon genome. We expect that the next decade will witness many new studies that take full advantage of salmonid genomic resources.

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.

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.