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.

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.

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.

Polygenic Basis and the Role of Genome Duplication in Adaptation to Similar Selective Environments.

Genetic changes underlying adaptation vary greatly in terms of complexity and, within the same species, genetic responses to similar selective pressures may or may not be the same. We examine both complex (supergene) and simple (SNP) genetic variants occurring in populations of rainbow trout (Oncorhynchus mykiss) independently isolated from ocean access and compared them to each other and to an anadromous below-barrier population representing their ancestral source to search for signatures of both parallel and nonparallel adaptation. All landlocked populations displayed an increased frequency of a large inversion on chromosome Omy05, while 3 of the 4 populations exhibited elevated frequencies of another inversion located on chromosome Omy20. In addition, we identified numerous regions outside these 2 inversions that also show significant shifts in allele frequencies consistent with adaptive evolution. However, there was little concordance among above-barrier populations in these specific genomic regions under selection. In part, the lack of concordance appears to arise from ancestral autopolyploidy in rainbow trout that provides duplicate genomic regions of similar functional composition for selection to act upon. Thus, while selection acting on landlocked populations universally favors the resident ecotype, outside of the major chromosomal inversions, the resulting genetic changes are largely distinct among populations. Our results indicate that selection on standing genetic variation is likely the primary mode of rapid adaptation, and that both supergene complexes and individual loci contribute to adaptive evolution, further highlighting the diversity of adaptive genomic variation involved in complex phenotypic evolution.

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.

Rapid parallel evolution of standing variation in a single, complex, genomic region is associated with life history in steelhead/rainbow trout.

Rapid adaptation to novel environments may drive changes in genomic regions through natural selection. Such changes may be population-specific or, alternatively, may involve parallel evolution of the same genomic region in multiple populations, if that region contains genes or co-adapted gene complexes affecting the selected trait(s). Both quantitative and population genetic approaches have identified associations between specific genomic regions and the anadromous (steelhead) and resident (rainbow trout) life-history strategies of Oncorhynchus mykiss. Here, we use genotype data from 95 single nucleotide polymorphisms and show that the distribution of variation in a large region of one chromosome, Omy5, is strongly associated with life-history differentiation in multiple above-barrier populations of rainbow trout and their anadromous steelhead ancestors. The associated loci are in strong linkage disequilibrium, suggesting the presence of a chromosomal inversion or other rearrangement limiting recombination. These results provide the first evidence of a common genomic basis for life-history variation in O. mykiss in a geographically diverse set of populations and extend our knowledge of the heritable basis of rapid adaptation of complex traits in novel habitats.

Characterization and distribution of a 14-Mb chromosomal inversion in native populations of rainbow trout (Oncorhynchus mykiss).

Multiple studies in a range of taxa have found links between structural variants and the development of ecologically important traits. Such variants are becoming easier to find due, in large part, to the increase in the amount of genome-wide sequence data in nonmodel organisms. The salmonids (salmon, trout, and charr) are a taxonomic group with abundant genome-wide datasets due to their importance in aquaculture, fisheries, and variation in multiple ecologically important life-history traits. Previous research on rainbow trout (Oncorhynchus mykiss) has documented a large pericentric (∼55 Mb) chromosomal inversion (CI) on chromosome 5 (Omy05) and a second smaller (∼14 Mb) chromosome inversion on Omy20. While the Omy05 inversion appears to be associated with multiple adaptive traits, the inversion on Omy20 has received far less attention. In this study, we re-analyze RAD-seq and amplicon data from several populations of rainbow trout (O. mykiss) to better document the structure and geographic distribution of variation in the Omy20 CI. Moreover, we utilize phylogenomic techniques to characterize both the age- and the protein-coding gene content of the Omy20 CI. We find that the age of the Omy20 inversion dates to the early stages of O. mykiss speciation and predates the Omy05 inversion by ∼450,000 years. The 2 CIs differ further in terms of the frequency of the homokaryotypes. While both forms of the Omy05 CI are found across the eastern Pacific, the ancestral version of the Omy20 CI is restricted to the southern portion of the species range in California. Furthermore, the Omy20 inverted haplotype is comparable in genetic diversity to the ancestral form, whereas derived CIs typically show substantially reduced genetic diversity. These data contribute to our understanding of the age and distribution of a large CI in rainbow trout and provide a framework for researchers looking to document CIs in other nonmodel species.

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.

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

Understanding life history traits is an important first step in formulating effective conservation and management strategies. The use of artificial propagation and supplementation as such a strategy can have numerous effects on the supplemented natural populations and minimizing life history divergence is crucial in minimizing these effects. Here, we use single nucleotide polymorphism (SNP) genotypes for large-scale parentage analysis and pedigree reconstruction in a hatchery population of steelhead, the anadromous form of rainbow trout. Nearly complete sampling of the broodstock for several consecutive years in two hatchery programmes allowed inference about multiple aspects of life history. Reconstruction of cohort age distribution revealed a strong component of fish that spawn at 2 years of age, in contrast to programme goals and distinct from naturally spawning steelhead in the region, which raises a significant conservation concern. The first estimates of variance in family size for steelhead in this region can be used to calculate effective population size and probabilities of inbreeding, and estimation of iteroparity rate indicates that it is reduced by hatchery production. Finally, correlations between family members in the day of spawning revealed for the first time a strongly heritable component to this important life history trait in steelhead and demonstrated the potential for selection to alter life history traits rapidly in response to changes in environmental conditions. Taken together, these results demonstrate the extraordinary promise of SNP-based pedigree reconstruction for providing biological inference in high-fecundity organisms that is not easily achievable with traditional physical tags.

A complex phenotype in salmon controlled by a simple change in migratory timing.

Differentiation between ecotypes is usually presumed to be complex and polygenic. Seasonal patterns of life history in salmon are used to categorize them into ecotypes, which are often considered "distinct" animals. Using whole-genome sequencing and tribal fishery sampling of Chinook salmon, we show that a single, small genomic region is nearly perfectly associated with spawning migration timing but not with adiposity or sexual maturity, traits long perceived as central to salmon ecotypes. Distinct migration timing does not prevent interbreeding between ecotypes, which are the result of a simple, ancient polymorphism segregating within a diverse population. Our finding that a complex migratory phenotype results from a single gene region will facilitate conservation and restoration of this iconic fish.

Publisher Correction: Sex-dependent dominance maintains migration supergene in rainbow trout.

An amendment to this paper has been published and can be accessed via a link at the top of the paper.

Over the falls? Rapid evolution of ecotypic differentiation in steelhead/rainbow trout (Oncorhynchus mykiss).

Adaptation to novel habitats and phenotypic plasticity can be counteracting forces in evolution, but both are key characteristics of the life history of steelhead/rainbow trout (Oncorhynchus mykiss). Anadromous steelhead reproduce in freshwater river systems and small coastal streams but grow and mature in the ocean. Resident rainbow trout, either sympatric with steelhead or isolated above barrier dams or waterfalls, represent an alternative life-history form that lives entirely in freshwater. We analyzed population genetic data from 1486 anadromous and resident O. mykiss from a small stream in coastal California with multiple barrier waterfalls. Based on data from 18 highly variable microsatellite loci (He = 0.68), we conclude that the resident population above one barrier, Big Creek Falls, is the result of a recent anthropogenic introduction from the anadromous population of O. mykiss below the falls. Furthermore, fish from this above-barrier population occasionally descend over the falls and have established a genetically differentiated below-barrier subpopulation at the base of the falls, which appears to remain reproductively isolated from their now-sympatric anadromous ancestors. These results support a hypothesis of rapid evolution of a purely resident life history in the above-barrier population in response to strong selection against downstream movement.

Long-Term Conservation of Ohnologs Through Partial Tetrasomy Following Whole-Genome Duplication in Salmonidae.

Whole-genome duplications (WGDs) have occurred repeatedly and broadly throughout the evolutionary history of eukaryotes. However, the effects of WGD on genome function and evolution remain unclear. The salmonid WGD that occurred approximately 88 million years ago presents an excellent opportunity for studying the effects of WGD as ∼10-15% of each salmonid genome still exhibits tetrasomic inheritance. Herein, we utilized the rainbow trout (Oncorhynchus mykiss) genome assembly and brain transcriptome data to examine the fate of gene pairs (ohnologs) following the salmonid whole-genome duplication. We find higher sequence identity between ohnologs located within known tetrasomic regions than between ohnologs found in disomic regions, and that tetrasomically inherited ohnologs showed greater similarity in patterns of gene expression and per ohnolog were lower expressed, than disomically inherited ohnologs. Enrichment testing for Gene Ontology terms identified 49 over-represented terms in tetrasomically inherited ohnologs compared to disomic ohnologs. However, why these ohnologs are retained as tetrasomic is difficult to answer. It could be that we have identified salmonid specific "dangerous duplicates", that is, genes that cannot take on new roles following WGD. Alternatively, there may be adaptive advantages for retaining genes as functional duplicates in tetrasomic regions, as presumably, movement of these genes into disomic regions would affect both their sequence identity and their gene expression patterns.

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.

A New Single Nucleotide Polymorphism Database for Rainbow Trout Generated Through Whole Genome Resequencing.

Single-nucleotide polymorphisms (SNPs) are highly abundant markers, which are broadly distributed in animal genomes. For rainbow trout (Oncorhynchus mykiss), SNP discovery has been previously done through sequencing of restriction-site associated DNA (RAD) libraries, reduced representation libraries (RRL) and RNA sequencing. Recently we have performed high coverage whole genome resequencing with 61 unrelated samples, representing a wide range of rainbow trout and steelhead populations, with 49 new samples added to 12 aquaculture samples from AquaGen (Norway) that we previously used for SNP discovery. Of the 49 new samples, 11 were double-haploid lines from Washington State University (WSU) and 38 represented wild and hatchery populations from a wide range of geographic distribution and with divergent migratory phenotypes. We then mapped the sequences to the new rainbow trout reference genome assembly (GCA_002163495.1) which is based on the Swanson YY doubled haploid line. Variant calling was conducted with FreeBayes and SAMtools mpileup, followed by filtering of SNPs based on quality score, sequence complexity, read depth on the locus, and number of genotyped samples. Results from the two variant calling programs were compared and genotypes of the double haploid samples were used for detecting and filtering putative paralogous sequence variants (PSVs) and multi-sequence variants (MSVs). Overall, 30,302,087 SNPs were identified on the rainbow trout genome 29 chromosomes and 1,139,018 on unplaced scaffolds, with 4,042,723 SNPs having high minor allele frequency (MAF > 0.25). The average SNP density on the chromosomes was one SNP per 64 bp, or 15.6 SNPs per 1 kb. Results from the phylogenetic analysis that we conducted indicate that the SNP markers contain enough population-specific polymorphisms for recovering population relationships despite the small sample size used. Intra-Population polymorphism assessment revealed high level of polymorphism and heterozygosity within each population. We also provide functional annotation based on the genome position of each SNP and evaluate the use of clonal lines for filtering of PSVs and MSVs. These SNPs form a new database, which provides an important resource for a new high density SNP array design and for other SNP genotyping platforms used for genetic and genomics studies of this iconic salmonid fish species.

Evolutionary restoration potential evaluated through the use of a trait-linked genetic marker.

Human-driven evolution can impact the ecological role and conservation value of impacted populations. Most evolutionary restoration approaches focus on manipulating gene flow, but an alternative approach is to manipulate the selection regime to restore historical or desired trait values. Here we examined the potential utility of this approach to restore anadromous migratory behavior in coastal California steelhead trout (Oncorhynchus mykiss) populations. We evaluated the effects of natural and anthropogenic environmental variables on the observed frequency of alleles at a genomic marker tightly associated with migratory behavior across 39 steelhead populations from across California, USA. We then modeled the potential for evolutionary restoration at sites that have been impacted by anthropogenic barriers. We found that complete barriers such as dams are associated with major reductions in the frequency of anadromy-associated alleles. The removal of dams is therefore expected to restore anadromy significantly. Interestingly, accumulations of large numbers of partial barriers (passable under at least some flow conditions) were also associated with significant reductions in migratory allele frequencies. Restoration involving the removal of partial barriers could be evaluated alongside dam removal and fishway construction as a cost-effective tool to restore anadromous fish migrations. Results encourage broader consideration of in situ evolution during the development of habitat restoration projects.

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.

Estimating population structure under nonequilibrium conditions in a conservation context: continent-wide population genetics of the giant Amazon river turtle, Podocnemis expansa (Chelonia; Podocnemididae).

Giant Amazon river turtles, Podocnemis expansa, are indigenous to the Amazon, Orinoco, and Essequibo River basins, and are distributed across nearly the entire width of the South American continent. Although once common, their large size, high fecundity, and gregarious nesting, made P. expansa especially vulnerable to over-harvesting for eggs and meat. Populations have been severely reduced or extirpated in many areas throughout its range, and the species is now regulated under Appendix II of the Convention on International Trade in Endangered Species. Here, we analyse data from mitochondrial DNA sequence and multiple nuclear microsatellite markers with an array of complementary analytical methods. Results show that concordance from multiple data sets and analyses can provide a strong signal of population genetic structure that can be used to guide management. The general lack of phylogeographic structure but large differences in allele and haplotype frequencies among river basins is consistent with fragmented populations and female natal-river homing. Overall, the DNA data show that P. expansa populations lack a long history of genetic differentiation, but that each major tributary currently forms a semi-isolated reproductive population and should be managed accordingly.