Large-scale parentage analysis reveals reproductive patterns and heritability of spawn timing in a hatchery population of steelhead (Oncorhynchus mykiss)

Genetic divergence and one-way gene flow influence contemporary evolution and ecology of a partially migratory fish

Recent work has revealed the importance of contemporary evolution in shaping ecological outcomes. In particular, rapid evolutionary divergence between populations has been shown to impact the ecology of populations, communities, and ecosystems. While studies have focused largely on the role of adaptive divergence in generating ecologically important variation among populations, much less is known about the role of gene flow in shaping ecological outcomes. After divergence, populations may continue to interact through gene flow, which may influence evolutionary and ecological processes. Here, we investigate the role of gene flow in shaping the contemporary evolution and ecology of recently diverged populations of anadromous steelhead and resident rainbow trout (Oncorhynchus mykiss). Results show that resident rainbow trout introduced above waterfalls have diverged evolutionarily from downstream anadromous steelhead, which were the source of introductions. However, the movement of fish from above to below the waterfalls has facilitated gene flow, which has reshaped genetic and phenotypic variation in the anadromous source population. In particular, gene flow has led to an increased frequency of residency, which in turn has altered population density, size structure, and sex ratio. This result establishes gene flow as a contemporary evolutionary process that can have important ecological outcomes. From a management perspective, anadromous steelhead are generally regarded as a higher conservation priority than resident rainbow trout, even when found within the same watershed. Our results show that anadromous and resident O. mykiss populations may be connected via gene flow, with important ecological consequences. Such eco-evolutionary processes should be considered when managing recently diverged populations connected by gene flow.

Ancestry and genetic structure of resident and anadromous rainbow trout (Oncorhynchus mykiss) in Argentina

Since the first introduction from North America more than a century ago, rainbow trout (Oncorhynchus mykiss) have rapidly established self-sustaining populations in major river basins of Patagonia. Many generations later, only the freshwater resident life history is expressed in the Chubut and Negro rivers of northern Argentinian Patagonia, whereas both the resident and anadromous life histories are found in the Santa Cruz River of southern Argentina. Despite previous studies that have tried to identify the sources of these introduced populations, uncertainty still exists. Here we combined data from many single-nucleotide polymorphisms and microsatellite loci in O. mykiss populations from Argentina and North America to evaluate putative source populations, gene flow between Argentinian river basins, and genetic diversity differences between Argentinian and North American populations. We found that populations from northern and southern Patagonia are highly differentiated and have limited gene flow between them. Phylogeographic analysis also confirmed that they have separate origins, with the northern populations most closely related to the domesticated rainbow trout strains that are raised worldwide and the Santa Cruz River populations most closely related to North American populations from California and Oregon that have an anadromous component. In addition, fish with different life histories in the Santa Cruz River were found to constitute a single interbreeding population. No evidence was found of reduced genetic variation in introduced rainbow trout, suggesting multiple contributing sources. In spite of these advances in understanding, significant questions remain regarding the origins and evolution of the introduced O. mykiss in Patagonia.

Genetic parentage reveals the (un)natural history of Central Valley hatchery steelhead

Populations composed of individuals descended from multiple distinct genetic lineages often feature significant differences in phenotypic frequencies. We considered hatchery production of steelhead, the migratory anadromous form of the salmonid species Oncorhynchus mykiss, and investigated how differences among genetic lineages and environmental variation impacted life history traits. We genotyped 23,670 steelhead returning to the four California Central Valley hatcheries over 9 years from 2011 to 2019, confidently assigning parentage to 13,576 individuals to determine age and date of spawning and rates of iteroparity and repeat spawning within each year. We found steelhead from different genetic lineages showed significant differences in adult life history traits despite inhabiting similar environments. Differences between coastal and Central Valley steelhead lineages contributed to significant differences in age at return, timing of spawning, and rates of iteroparity among programs. In addition, adaptive genomic variation associated with life history development in this species varied among hatchery programs and was associated with the age of steelhead spawners only in the coastal lineage population. Environmental variation likely contributed to variations in phenotypic patterns observed over time, as our study period spanned both a marine heatwave and a serious drought in California. Our results highlight evidence of a strong genetic component underlying known phenotypic differences in life history traits between two steelhead lineages.

Distinct patterns of inheritance shape life-history traits in steelhead trout

Life-history variation is the raw material of adaptation, and understanding its genetic and environmental underpinnings is key to designing effective conservation strategies. We used large-scale genetic pedigree reconstruction of anadromous steelhead trout (Oncorhynchus mykiss) from the Russian River, CA, USA, to elucidate sex-specific patterns of life-history traits and their heritability. SNP data from adults returning from sea over a 14-year period were used to identify 13,474 parent-offspring trios. These pedigrees were used to determine age structure, size distributions and family sizes for these fish, as well as to estimate the heritability of two key life-history traits, spawn date and age at maturity (first reproduction). Spawn date was highly heritable (h2  = 0.73) and had a cross-sex genetic correlation near unity. We provide the first estimate of heritability for age at maturity in ocean-going fish from this species and found it to be highly heritable (h2 from 0.29 to 0.62, depending on sex and method), with a much lower genetic correlation across sexes. We also evaluated genotypes at a migration-associated inversion polymorphism and found sex-specific correlations with age at maturity. The significant heritability of these two key reproductive traits in these imperiled fish, and their patterns of inheritance in the two sexes, is consistent with predictions of both natural and sexually antagonistic selection (sexes experience opposing selection pressures). This emphasizes the importance of anthropogenic factors, including hatchery practices and ecosystem modifications, in shaping the fitness of this species, thus providing important guidance for management and conservation efforts.

Loss of plasticity in maturation timing after ten years of captive spawning in a delta smelt conservation hatchery

Adaptation to captivity in spawning programs can lead to unintentional consequences, such as domestication that results in reduced fitness in the wild. The timing of sexual maturation has been shown to be a trait under domestication selection in fish hatcheries, which affects a fish's access to mating opportunities and aligning their offspring's development with favorable environmental conditions. Earlier maturing fish may be favored in hatchery settings where managers provide artificially optimal growing conditions, but early maturation may reduce fitness in the wild if, for example, there is a mismatch between timing of reproduction and availability of resources that support recruitment. We investigated patterns of maturation timing in a delta smelt (Hypomesus transpacificus) conservation hatchery by quantifying changes to the median age at maturity since the captive spawning program was initiated in 2008. Over the span of a decade, we observed a small, but significant increase in age at maturity among broodstock by 2.2 weeks. This trait had low heritability and was largely controlled by phenotypic plasticity that was dependent on the time of year fish were born. Fish that were born later in the year matured faster, potentially a carryover from selection favoring synchronous spawning in the wild. However, higher DI (domestication index) fish showed a loss of plasticity, we argue, as a result of hatchery practices that breed individuals past peak periods of female ripeness. Our findings suggest that the hatchery setting has relaxed selection pressures for fish to mature quickly at the end of the year and, consequently, has led to a loss of plasticity in age at maturity. Hatchery fish that are re-introduced in the wild may not be able to align maturation with population peaks if their maturation rates are too slow with reduced plasticity, potentially resulting in lower fitness.

Sex-specific heritabilities for length at maturity among Pacific salmonids and their consequences for evolution in response to artificial selection

Artificial selection, whether intentional or coincidental, is a common result of conservation policies and natural resource management. To reduce unintended consequences of artificial selection, conservation practitioners must understand both artificial selection gradients on traits of interest and how those traits are correlated with others that may affect population growth and resilience. We investigate how artificial selection on male body size in Pacific salmon (Oncorhynchus spp.) may influence the evolution of female body size and female fitness. While salmon hatchery managers often assume that selection for large males will also produce large females, this may not be the case-in fact, because the fastest-growing males mature earliest and at the smallest size, and because female age at maturity varies little, small males may produce larger females if the genetic architecture of growth rate is the same in both sexes. We explored this possibility by estimating sex-specific heritability values of and natural and artificial selection gradients on length at maturity in four populations representing three species of Pacific salmon. We then used the multivariate breeder's equation to project how artificial selection against small males may affect the evolution of female length and fecundity. Our results indicate that the heritability of length at maturity is greater within than between the sexes and that sire-daughter heritability values are especially small. Salmon hatchery policies should consider these sex-specific quantitative genetic parameters to avoid potential unintended consequences of artificial selection.

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.

Parentage-based tagging combined with genetic stock identification is a cost-effective and viable replacement for coded-wire tagging in large-scale assessments of marine Chinook salmon fisheries in British Columbia, Canada

Wild Pacific salmon, including Chinook salmon Oncorhynchus tshawytscha, have been supplemented with hatchery propagation for over 50 years in support of increased ocean harvest, mitigation for hydroelectric development, and conservation of threatened populations. In Canada, the Wild Salmon Policy for Pacific salmon was established with the goal of maintaining and restoring healthy and diverse Pacific salmon populations, making conservation of wild salmon and their habitats the highest priority for resource management decision-making. For policy implementation, a new approach to the assessment and management of Chinook salmon and the associated hatchery production and fisheries management are needed. Implementation of genetic stock identification (GSI) and parentage-based tagging (PBT) for marine fisheries assessment may overcome problems associated with coded-wire tag-based (CWT) assessment and management of Chinook salmon fisheries, providing at a minimum information equivalent to that derived from the CWT program. GSI and PBT were used to identify Chinook salmon sampled in 2018 and 2019 marine fisheries (18,819 individuals genotyped) in British Columbia to specific conservation units (CU), populations, and broodyears. Individuals were genotyped at 391 single nucleotide polymorphisms via direct sequencing of amplicons. Very high accuracy of assignment to population and age (>99.5%) via PBT was observed for 1994 Chinook salmon of ages 2-4 years, with a 105,722-individual, 380-population baseline available for assignment. Application of a GSI-PBT system of identification to individuals in 2019 fisheries provided high-resolution estimates of stock composition, catch, and exploitation rate by CU or population, with fishery exploitation rates directly comparable to those provided by CWTs for 13 populations. GSI and PBT provide an alternate, cheaper, and more effective method in the assessment and management of Canadian-origin Chinook salmon relative to CWTs, and an opportunity for a genetics-based system to replace the current CWT system for salmon assessment.

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.

Differential DNA methylation in somatic and sperm cells of hatchery vs wild (natural-origin) steelhead trout populations

Environmental factors such as nutrition, stress, and toxicants can influence epigenetic programming and phenotypes of a wide variety of species from plants to humans. The current study was designed to investigate the impacts of hatchery spawning and rearing on steelhead trout (Oncorhynchus mykiss) vs the wild fish on a molecular level. Additionally, epigenetic differences between feeding practices that allow slow growth (2 years) and fast growth (1 year) hatchery trout were investigated. The sperm and red blood cells (RBC) from adult male slow growth/maturation hatchery steelhead, fast growth/maturation hatchery steelhead, and wild (natural-origin) steelhead were collected for DNA preparation to investigate potential alterations in differential DNA methylation regions (DMRs) and genetic mutations, involving copy number variations (CNVs). The sperm and RBC DNA both had a large number of DMRs when comparing the hatchery vs wild steelhead trout populations. The DMRs were cell type specific with negligible overlap. Slow growth/maturation compared to fast growth/maturation steelhead also had a larger number of DMRs in the RBC samples. A number of the DMRs had associated genes that were correlated to various biological processes and pathologies. Observations demonstrate a major epigenetic programming difference between the hatchery and wild natural-origin fish populations, but negligible genetic differences. Therefore, hatchery conditions and growth/maturation rate can alter the epigenetic developmental programming of the steelhead trout. Interestingly, epigenetic alterations in the sperm allow for potential epigenetic transgenerational inheritance of phenotypic variation to future generations. The impacts of hatchery exposures are not only important to consider on the fish exposed, but also on future generations and evolutionary trajectory of fish in the river populations.

Rapture facilitates inexpensive and high-throughput parent-based tagging in salmonids

Accurate methods for tracking individuals are crucial to the success of fisheries and aquaculture management. Management of migratory salmonid populations, which are important for the health of many economies, ecosystems, and indigenous cultures, is particularly dependent on data gathered from tagged fish. However, the physical tagging methods currently used have many challenges including cost, variable marker retention, and information limited to tagged individuals. Genetic tracking methods combat many of the problems associated with physical tags, but have their own challenges including high cost, potentially difficult marker design, and incompatibility of markers across species. Here we show the feasibility of a new genotyping method for parent-based tagging (PBT), where individuals are tracked through the inherent genetic relationships with their parents. We found that Rapture sequencing, a combination of restriction-site associated DNA and capture sequencing, provides sufficient data for parentage assignment. Additionally, the same capture bait set, which targets specific restriction-site associated DNA loci, can be used for both Rainbow Trout Oncorhynchus mykiss and Chinook Salmon Oncorhynchus tshawytscha. We input 248 single nucleotide polymorphisms from 1,121 samples to parentage assignment software and compared parent-offspring relationships of the spawning pairs recorded in a hatchery. Interestingly, our results suggest sperm contamination during hatchery spawning occurred in the production of 14% of offspring, further confirming the need for genetic tagging in accurately tracking individuals. PBT with Rapture successfully assigned progeny to parents with a 98.86% accuracy with sufficient genetic data. Cost for this pilot study was approximately $3 USD per sample. As costs vary based on the number of markers used and individuals sequenced, we expect that when implemented at a large-scale, per sample costs could be further decreased. We conclude that Rapture PBT provides a cost-effective and accurate alternative to the physical coded wire tags, and other genetic-based methods.

RAPTURE (RAD capture) panel facilitates analyses characterizing sea lamprey reproductive ecology and movement dynamics

Genomic tools are lacking for invasive and native populations of sea lamprey (Petromyzon marinus). Our objective was to discover single nucleotide polymorphism (SNP) loci to conduct pedigree analyses to quantify reproductive contributions of adult sea lampreys and dispersion of sibling larval sea lampreys of different ages in Great Lakes tributaries. Additional applications of data were explored using additional geographically expansive samples. We used restriction site-associated DNA sequencing (RAD-Seq) to discover genetic variation in Duffins Creek (DC), Ontario, Canada, and the St. Clair River (SCR), Michigan, USA. We subsequently developed RAD capture baits to genotype 3,446 RAD loci that contained 11,970 SNPs. Based on RAD capture assays, estimates of variance in SNP allele frequency among five Great Lakes tributary populations (mean F ST 0.008; range 0.00-0.018) were concordant with previous microsatellite-based studies; however, outlier loci were identified that contributed substantially to spatial population genetic structure. At finer scales within streams, simulations indicated that accuracy in genetic pedigree reconstruction was high when 200 or 500 independent loci were used, even in situations of high spawner abundance (e.g., 1,000 adults). Based on empirical collections of larval sea lamprey genotypes, we found that age-1 and age-2 families of full and half-siblings were widely but nonrandomly distributed within stream reaches sampled. Using the genomic scale set of SNP loci developed in this study, biologists can rapidly genotype sea lamprey in non-native and native ranges to investigate questions pertaining to population structuring and reproductive ecology at previously unattainable scales.

Sex-dependent dominance maintains migration supergene in rainbow trout

Males and females often differ in their fitness optima for shared traits that have a shared genetic basis, leading to sexual conflict. Morphologically differentiated sex chromosomes can resolve this conflict and protect sexually antagonistic variation, but they accumulate deleterious mutations. However, how sexual conflict is resolved in species that lack differentiated sex chromosomes is largely unknown. Here we present a chromosome-anchored genome assembly for rainbow trout (Oncorhynchus mykiss) and characterize a 55-Mb double-inversion supergene that mediates sex-specific migratory tendency through sex-dependent dominance reversal, an alternative mechanism for resolving sexual conflict. The double inversion contains key photosensory, circadian rhythm, adiposity and sex-related genes and displays a latitudinal frequency cline, indicating environmentally dependent selection. Our results show sex-dependent dominance reversal across a large autosomal supergene, a mechanism for sexual conflict resolution capable of protecting sexually antagonistic variation while avoiding the homozygous lethality and deleterious mutations associated with typical heteromorphic sex chromosomes.

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

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

Artificial selection on reproductive timing in hatchery salmon drives a phenological shift and potential maladaptation to climate change

The timing of breeding migration and reproduction links generations and substantially influences individual fitness. In salmonid fishes, such phenological events (seasonal return to freshwater and spawning) vary among populations but are consistent among years, indicating local adaptation in these traits to prevailing environmental conditions. Changing reproductive phenology has been observed in many populations of Atlantic and Pacific salmon and is sometimes attributed to adaptive responses to climate change. The sockeye salmon spawning in the Cedar River near Seattle, Washington, USA, have displayed dramatic changes in spawning timing over the past 50 years, trending later through the early 1990s, and becoming earlier since then. We explored the patterns and drivers of these changes using generalized linear models and mathematical simulations to identify possible environmental correlates of the changes, and test the alternative hypothesis that hatchery propagation caused inadvertent selection on timing. The trend toward later spawning prior to 1993 was partially explained by environmental changes, but the rapid advance in spawning since was not. Instead, since its initiation in 1991, the hatchery has, on average, selected for earlier spawning, and, depending on trait heritability, could have advanced spawning by 1-3 weeks over this period. We estimated heritability of spawning date to be high (h 2 ~0.8; 95% CI: 0.5-1.1), so the upper end of this range is not improbable, though at lower heritabilities a smaller effect would be expected. The lower reproductive success of early spawners and relatively low survival of early emerging juveniles observed in recent years suggest that artificial and natural selection are acting in opposite directions. The fitness costs of early spawning may be exacerbated by future warming; thus, the artificially advanced phenology could reduce the population's productivity. Such artificial selection is known in many salmon hatcheries, so there are broad consequences for the productivity of wild populations comingled with hatchery-produced fish.

The future of parentage analysis: From microsatellites to SNPs and beyond

Parentage analysis is a cornerstone of molecular ecology that has delivered fundamental insights into behaviour, ecology and evolution. Microsatellite markers have long been the king of parentage, their hypervariable nature conferring sufficient power to correctly assign offspring to parents. However, microsatellite markers have seen a sharp decline in use with the rise of next-generation sequencing technologies, especially in the study of population genetics and local adaptation. The time is ripe to review the current state of parentage analysis and see how it stands to be affected by the emergence of next-generation sequencing approaches. We find that single nucleotide polymorphisms (SNPs), the typical next-generation sequencing marker, remain underutilized in parentage analysis but are gaining momentum, with 58 SNP-based parentage analyses published thus far. Many of these papers, particularly the earlier ones, compare the power of SNPs and microsatellites in a parentage context. In virtually every case, SNPs are at least as powerful as microsatellite markers. As few as 100-500 SNPs are sufficient to resolve parentage completely in most situations. We also provide an overview of the analytical programs that are commonly used and compatible with SNP data. As the next-generation parentage enterprise grows, a reliance on likelihood and Bayesian approaches, as opposed to strict exclusion, will become increasingly important. We discuss some of the caveats surrounding the use of next-generation sequencing data for parentage analysis and conclude that the future is bright for this important realm of molecular ecology.

Comparison of coded-wire tagging with parentage-based tagging and genetic stock identification in a large-scale coho salmon fisheries application in British Columbia, Canada

Wild Pacific salmon, including Coho salmon Onchorynchus kisutch, have been supplemented with hatchery propagation for over 50 years in support of increased ocean harvest and conservation of threatened populations. In Canada, the Wild Salmon Policy for Pacific salmon was established with the goal of maintaining and restoring healthy and diverse Pacific salmon populations, making conservation of wild salmon and their habitats the highest priority for resource management decision-making. A new approach to the assessment and management of wild coho salmon, and the associated hatchery production and fishery management is needed. Implementation of parentage-based tagging (PBT) may overcome problems associated with coded-wire tag-based (CWT) assessment and management of coho salmon fisheries, providing at a minimum information equivalent to that derived from the CWT program. PBT and genetic stock identification (GSI) were used to identify coho salmon sampled in fisheries (8,006 individuals) and escapements (1,692 individuals) in British Columbia to specific conservation units (CU), populations, and broodyears. Individuals were genotyped at 304 single nucleotide polymorphisms (SNPs) via direct sequencing of amplicons. Very high accuracy of assignment to population (100%) via PBT for 543 jack (age 2) assigned to correct age and collection location and 265 coded-wire tag (CWT, age 3) coho salmon assigned to correct age and release location was observed, with a 40,774-individual, 267-population baseline available for assignment. Coho salmon from un-CWTed enhanced populations contributed 65% of the catch in southern recreational fisheries in 2017. Application of a PBT-GSI system of identification to individuals in 2017 fisheries and escapements provided high-resolution estimates of stock composition, catch, and exploitation rate by CU or population, providing an alternate and more effective method in the assessment and management of Canadian-origin coho salmon relative to CWTs, and an opportunity for a genetic-based system to replace the current CWT system for coho salmon assessment.

Anadromy Redux? Genetic Analysis to Inform Development of an Indigenous American River Steelhead Broodstock

The construction of dams and water diversions has severely limited access to spawning habitat for anadromous fishes. To mitigate for these impacts, hatchery programs rear and release millions of juvenile salmonids, including steelhead, the anadromous ecotype of the species Oncorhynchus mykiss. These programs sometimes use nonindigenous broodstock sources that may have negative effects on wild populations. In California, however, only one anadromous fish hatchery program currently uses nonnative broodstock: the steelhead program at Nimbus Fish Hatchery on the American River, a tributary of the Sacramento River in the California Central Valley. The goal of this study was to determine if potentially appropriate sources to replace the broodstock for the Nimbus Hatchery steelhead program exist in the Upper American River, above Nimbus and Folsom dams. We show that all Upper American River O. mykiss sampled share ancestry with other populations in the Central Valley steelhead distinct population segment, with limited introgression from out-of-basin sources in some areas. Furthermore, some Upper American River populations retain adaptive genomic variation associated with a migratory life history, supporting the hypothesis that these populations display adfluvial migratory behavior. Together, these results provide insights into the evolution of trout populations above barrier dams. We conclude that some Upper American River O. mykiss populations represent genetically appropriate sources from which fisheries managers could potentially develop a new broodstock for the Nimbus Hatchery steelhead program to reestablish a native anadromous population in the Lower American River and contribute to recovery of the threatened Central Valley steelhead distinct population segment.

Saving the spandrels? Adaptive genomic variation in conservation and fisheries management

As highlighted by many of the papers in this issue, research on the genomic basis of adaptive phenotypic variation in natural populations has made spectacular progress in the past few years, largely due to the advances in sequencing technology and analysis. Without question, the resulting genomic data will improve the understanding of regions of the genome under selection and extend knowledge of the genetic basis of adaptive evolution. What is far less clear, but has been the focus of active discussion, is how such information can or should transfer into conservation practice to complement more typical conservation applications of genetic data. Before such applications can be realized, the evolutionary importance of specific targets of selection relative to the genome-wide diversity of the species as a whole must be evaluated. The key issues for the incorporation of adaptive genomic variation in conservation and management are discussed here, using published examples of adaptive genomic variation associated with specific phenotypes in salmonids and other taxa to highlight practical considerations for incorporating such information into conservation programmes. Scenarios are described in which adaptive genomic data could be used in conservation or restoration, constraints on its utility and the importance of validating inferences drawn from new genomic data before applying them in conservation practice. Finally, it is argued that an excessive focus on preserving the adaptive variation that can be measured, while ignoring the vast unknown majority that cannot, is a modern twist on the adaptationist programme that Gould and Lewontin critiqued almost 40 years ago.

Ancestry and adaptive evolution of anadromous, resident, and adfluvial rainbow trout (Oncorhynchus mykiss) in the San Francisco bay area: application of adaptive genomic variation to conservation in a highly impacted landscape

The streams draining of into San Francisco Bay, California, have been impacted by habitat alteration for over 150 years, and roads, dams, water diversions, and other impediments now block the paths of many aquatic migratory species. These changes can affect the genetic structure of fish populations, as well as driving adaptive evolution to novel environmental conditions. Here, we determine the evolutionary relationships of San Francisco Bay Area steelhead/rainbow trout (Oncorhynchus mykiss) populations and show that (i) they are more closely related to native coastal steelhead than to the California Central Valley lineage, with no evidence of introgression by domesticated hatchery rainbow trout, (ii) populations above and below barriers within watersheds are each other's closest relatives, and (iii) adaptive genomic variation associated with migratory life-history traits in O. mykiss shows substantial evolutionary differences between fish above and below dams. These findings support continued habitat restoration and protection of San Francisco Bay Area O. mykiss populations and demonstrate that ecological conditions in novel habitats above barriers to anadromy influence life-history evolution. We highlight the importance of considering the adaptive landscape in conservation and restoration programs for species living in highly modified habitats, particularly with respect to key life-history traits.

Effective number of breeders from sibship reconstruction: empirical evaluations using hatchery steelhead

Effective population size (N_e) is among the most important metrics in evolutionary biology. In natural populations, it is often difficult to collect adequate demographic data to calculate N_e directly. Consequently, genetic methods to estimate N_e have been developed. Two N_e estimators based on sibship reconstruction using multilocus genotype data have been developed in recent years: sibship assignment and parentage analysis without parents. In this study, we evaluated the accuracy of sibship reconstruction using a large empirical dataset from five hatchery steelhead populations with known pedigrees and using 95 single nucleotide polymorphism (SNP) markers. We challenged the software COLONY with 2,599,961 known relationships and demonstrated that reconstruction of full‐sib and unrelated pairs was greater than 95% and 99% accurate, respectively. However, reconstruction of half‐sib pairs was poor (<5% accurate). Despite poor half‐sib reconstruction, both estimators provided accurate estimates of the effective number of breeders (N_b) when sample sizes were near or greater than the true N_b and when assuming a monogamous mating system. We further demonstrated that both methods provide roughly equivalent estimates of N_b. Our results indicate that sibship reconstruction and current SNP panels provide promise for estimating N_b in steelhead populations in the region.

Shifting Thresholds: Rapid Evolution of Migratory Life Histories in Steelhead/Rainbow Trout, Oncorhynchus mykiss

Expression of phenotypic plasticity depends on reaction norms adapted to historic selective regimes; anthropogenic changes in these selection regimes necessitate contemporary evolution or declines in productivity and possibly extinction. Adaptation of conditional strategies following a change in the selection regime requires evolution of either the environmentally influenced cue (e.g., size-at-age) or the state (e.g., size threshold) at which an individual switches between alternative tactics. Using a population of steelhead (Oncorhynchus mykiss) introduced above a barrier waterfall in 1910, we evaluate how the conditional strategy to migrate evolves in response to selection against migration. We created 9 families and 917 offspring from 14 parents collected from the above- and below-barrier populations. After 1 year of common garden-rearing above-barrier offspring were 11% smaller and 32% lighter than below-barrier offspring. Using a novel analytical approach, we estimate that the mean size at which above-barrier fish switch between the resident and migrant tactic is 43% larger than below-barrier fish. As a result, above-barrier fish were 26% less likely to express the migratory tactic. Our results demonstrate how rapid and opposing changes in size-at-age and threshold size contribute to the contemporary evolution of a conditional strategy and indicate that migratory barriers may elicit rapid evolution toward the resident life history on timescales relevant for conservation and management of conditionally migratory species.

Quantitative genetics of migration-related traits in rainbow and steelhead trout

Rainbow trout (Oncorhynchus mykiss) exhibit remarkable life history diversity throughout their native range, and among the most evident is variation in migratory propensity. Although some populations and ecotypes will remain resident in freshwater habitats throughout their life history, others have the ability to undertake tremendous marine migrations. Those that migrate undergo a suite of behavioral, morphological, and physiological adaptations in a process called smoltification. We describe a quantitative genetic analysis of 22 growth, size, and morphological traits in addition to overall life history classification (resident or migrant) over the temporal process of smoltification in a large multi-generation experimental pedigree (n = 16,139) of migratory and resident rainbow trout derived from a wild population, which naturally segregates for migratory propensity. We identify significant additive genetic variance and covariance among the suite of traits that make up a component of the migratory syndrome in this species. Additionally, we identify high heritability estimates for the life history classifications and observe a strong negative genetic correlation between the migratory and resident life history trajectories. Given the large heritability estimates of all of the traits that segregate between migratory and resident rainbow trout, we conclude that these traits can respond to selection. However, given the high degree of genetic correlation between these traits, they do not evolve in isolation, but rather as a suite of coordinated characters in a predictable manner.

Time scale matters: genetic analysis does not support adaptation-by-time as the mechanism for adaptive seasonal declines in kokanee reproductive life span

Seasonal declines of fitness-related traits are often attributed to environmental effects or individual-level decisions about reproductive timing and effort, but genetic variation may also play a role. In populations of Pacific salmon (Oncorhynchus spp.), seasonal declines in reproductive life span have been attributed to adaptation-by-time, in which divergent selection for different traits occurs among reproductively isolated temporal components of a population. We evaluated this hypothesis in kokanee (freshwater obligate Oncorhynchus nerka) by testing for temporal genetic structure in neutral and circadian-linked loci. We detected no genetic differences in presumably neutral loci among kokanee with different arrival and maturation dates within a spawning season. Similarly, we detected no temporal genetic structure in OtsClock1b, Omy1009uw, or OmyFbxw11, candidate loci associated with circadian function. The genetic evidence from this study and others indicates a lack of support for adaptation-by-time as an important evolutionary mechanism underlying seasonal declines in reproductive life span and a need for greater consideration of other mechanisms such as time-dependent, adaptive adjustment of reproductive effort.

Genetic assessment of a summer chum salmon metapopulation in recovery

Programs to rebuild imperiled wild fish populations often include hatchery-born fish derived from wild populations to supplement natural spawner abundance. These programs require monitoring to determine their demographic, biological, and genetic effects. In 1990s in Washington State, the Summer Chum Salmon Conservation Initiative developed a recovery program for the threatened Hood Canal summer chum salmon Evolutionarily Significant Unit (ESU) (the metapopulation) that used in-river spawners (wild fish) for each respective supplementation broodstock in six tributaries. Returning spawners (wild-born and hatchery-born) composed subsequent broodstocks, and tributary-specific supplementation was limited to three generations. We assessed impacts of the programs on neutral genetic diversity in this metapopulation using 16 microsatellite loci and a thirty-year dataset spanning before and after supplementation, roughly eight generations. Following supplementation, differentiation among subpopulations decreased (but not significantly) and isolation by distance patterns remained unchanged. There was no decline in genetic diversity in wild-born fish, but hatchery-born fish sampled in the same spawning areas had significantly lower genetic diversity and unequal family representation. Despite potential for negative effects from supplementation programs, few were detected in wild-born fish. We hypothesize that chum salmon natural history makes them less vulnerable to negative impacts from hatchery supplementation.