Efficient population representation with more genetic markers increases performance of a steelhead (Oncorhynchus mykiss) genetic stock identification baseline.

Genetic stock identification (GSI) is an important fisheries management tool to identify the origin of fish harvested in mixed stock fisheries. Periodic updates of genetic baselines can improve performance via the addition of unsampled or under-sampled populations and the inclusion of more informative markers. We used a combination of baselines to evaluate how population representation, marker number, and marker type affected the performance and accuracy of genetic stock assignments (self-assignment, bias, and holdout group tests) for steelhead (Oncorhynchus mykiss) in the Snake River basin. First, we compared the performance of an existing genetic baseline with a newly developed one which had a reduced number of individuals from more populations using the same set of markers. Self-assignment rates were significantly higher (p < 0.001; +5.4%) for the older, larger baseline, bias did not differ significantly between the two, but there was a significant improvement in performance for the new baseline in holdout results (p < 0.001; mean increase of 25.0%). Second, we compared the performance of the new baseline with increased numbers of genetic markers (~2x increase of single-nucleotide polymorphisms; SNPs) for the same set of baseline individuals. In this comparison, results produced significantly higher rates of self-assignment (p < 0.001; +9.7%) but neither bias nor leave-one-out were significantly affected. Third, we compared 334 SNPs versus opportunistically discovered microhaplotypes from the same amplicons for the new baseline, and showed the latter produced significantly higher rates of self-assignment (p < 0.01; +2.6%), similar bias, but slightly lower holdout performance (-0.1%). Combined, we show the performance of genetic baselines can be improved via representative and efficient sampling, that increased marker number consistently improved performance over the original baseline, and that opportunistic discovery of microhaplotypes can lead to small improvements in GSI performance.

Contrasting patterns of sequence variation in steelhead populations reflect distinct evolutionary processes.

Multiple evolutionary processes influence genome-wide allele frequencies and quantifying effects of genetic drift, and multiple forms of selection remain challenging in natural populations. Here, we investigate variation at major effect loci in contrast to patterns of neutral drift across a wide collection of steelhead (Oncorhynchus mykiss) populations that have declined in abundance due to anthropogenic impacts. Whole-genome resequencing of 74 populations of steelhead revealed genome-wide patterns (~8 million SNPs) consistent with expected neutral population structure. However, allelic variation at major effect loci associated with adult migration timing (chromosome 28: GREB1L/ROCK1) and age at maturity (chromosome 25: SIX6) reflected how selection has acted on phenotypic variation in contrast with neutral structure. Variation at major effect loci was influenced by evolutionary processes with differing signals between the strongly divergent Coastal and Inland lineages, while allele frequencies within and among populations within the Inland lineage have been driven by local natural selection as well as recent anthropogenic influences. Recent anthropogenic effects appeared to have influenced the frequency of major effect alleles including artificial selection for specific traits in hatchery stocks with subsequent gene flow into natural populations. Selection from environmental factors at various scales has also likely influenced variation for major effect alleles. These results reveal evolutionary mechanisms that influence allele frequencies at major effect loci that are critical for conservation of phenotypic traits and life history variation of this protected species.

Single-parentage assignments reveal negative-assortative mating in an endangered salmonid.

Understanding reproductive patterns in endangered species is critical for supporting their recovery efforts. In this study we use a combination of paired-parent and single-parent assignments to examine the reproductive patterns in an endangered population of sockeye salmon (Oncorhynchus nerka) that uses Redfish Lake in central Idaho as a spawning and nursery lake. Recovery efforts include the release of maturing adults into the lake for volitional spawning. The lake is also inhabited by a population of resident O. nerka that is genetically indistinguishable, but phenotypically smaller, to the maturing adults released into the lake. The resident population is difficult to sample and the reproductive patterns between the two groups are unknown. We used results of paired- and single-parentage assignments to specifically examine the reproductive patterns of male fish released into the lake under an equal sex ratio and a male-biased sex ratio. Assignment results of offspring leaving the lake indicated a reproductive shift by males under the two scenarios. Males displayed an assortative mating pattern under an equal sex ratio and spawned almost exclusively with the released females. Under a male-biased sex ratio most males shifted to a negative-assortative mating pattern and spawned with smaller females from the resident population. These males were younger and smaller than males that spawned with released females suggesting they were unable to compete with larger males for spawning opportunities with the larger, released females. The results provided insights into the reproductive behavior of this endangered population and has implications for recovery efforts.

Sex-specific markers undetected in green sunfish Lepomis cyanellus using restriction-site associated DNA sequencing.

We used restriction-site associated DNA sequencing for SNP discovery and genotyping of known-sex green sunfish Lepomis cyanellus DNA samples to search for sex-diagnostic single nucleotide polymorphisms (SNPs) and restriction-site associated sequences present in one sex and absent in the other. The bioinformatic analyses discovered candidate SNPs and sex-specific restriction-site associated sequences that fit patterns of male or female heterogametic sex determination systems. However, when primers were developed and tested, no candidates reliably identified phenotypic sex. The top performing SNP candidate (ZW_218) correlated with phenotypic sex 63.0% of the time and the presence-absence loci universally amplified in both sexes. We recommend further investigations that interrogate a larger fraction of the L. cyanellus genome. Additionally, studies on the effect of temperature and rearing density on sex determination, as well as breeding of sex-reversed individuals, could provide more insights into the sex determination system of L. cyanellus.

Mating systems and predictors of relative reproductive success in a Cutthroat Trout subspecies of conservation concern.

Mating systems and patterns of reproductive success in fishes play an important role in ecology and evolution. While information on the reproductive ecology of many anadromous salmonids (Oncorhynchus spp.) is well detailed, there is less information for nonanadromous species including the Yellowstone Cutthroat Trout (O. clarkii bouvieri), a subspecies of recreational angling importance and conservation concern. Using data from a parentage-based tagging study, we described the genetic mating system of a migratory population of Yellowstone Cutthroat Trout, tested for evidence of sexual selection, and identified predictors of mating and reproductive success. The standardized variance in mating success (i.e., opportunity for sexual selection) was significantly greater for males relative to females, and while the relationship between mating success and reproductive success (i.e., Bateman gradient) was significantly positive for both sexes, a greater proportion of reproductive success was explained by mating success for males (r 2 = 0.80) than females (r 2 = 0.59). Overall, the population displayed a polygynandrous mating system, whereby both sexes experienced variation in mating success due to multiple mating, and sexual selection was variable across sexes. Tests for evidence of sexual selection indicated the interaction between mating success and total length best-predicted relative reproductive success. We failed to detect a signal of inbreeding avoidance among breeding adults, but the group of parents that produced progeny were on average slightly less related than adults that did not produce progeny. Lastly, we estimated the effective number of breeders (N b) and effective population size (N e) and identified while N b was lower than N e, both are sufficiently high to suggest Yellowstone Cutthroat Trout in Burns Creek represent a genetically stable and diverse population.

Contemporary genetic structure affects genetic stock identification of steelhead trout in the Snake River basin.

Genetic stock identification is a widely applied tool for the mixed-stock management of salmonid species throughout the North Pacific Rim. The effectiveness of genetic stock identification is dependent on the level of differentiation among stocks which is often high due to the life history of these species that involves high homing fidelity to their natal streams. However, the utility of this tool can be reduced when natural genetic structuring has been altered by hatchery translocation and/or supplementation. We examined the genetic population structure of ESA-listed steelhead in the Snake River basin of the United States. We analyzed 9,613 natural-origin adult steelhead returning to Passive Integrated Transponder detection sites throughout the basin from 2010 through 2017. Individuals were genotyped at 180 single nucleotide polymorphic genetic markers and grouped into 20 populations based on their return location. While we expected to observe a common pattern of hierarchical genetic structuring due to isolation by distance, we observed low genetic differentiation between populations in the upper Salmon River basin compared to geographically distant populations in the lower Snake River basin. These results were consistent with lower genetic stock assignment probabilities observed for populations in this upper basin. We attribute these patterns of reduced genetic structure to the translocation of lower basin steelhead stocks and ongoing hatchery programs in the upper Salmon River basin. We discuss the implications of these findings on the utility of genetic stock identification in the basin and discuss opportunities for increasing assignment probabilities in the face of low genetic structure.

Ecological segregation moderates a climactic conclusion to trout hybridization.

For decades, it has been assumed that introgressive hybridization between introduced rainbow trout and native cutthroat trout in western North America will lead to genomic extinction of the latter. A broad-scale re-examination of their interaction indicates that ecological differences between these species and demographic processes are dictating the location and extent of their hybrid zones, and that runaway introgression between these taxa is unlikely.

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

Effective population size (Ne ) is among the most important metrics in evolutionary biology. In natural populations, it is often difficult to collect adequate demographic data to calculate Ne directly. Consequently, genetic methods to estimate Ne have been developed. Two Ne 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 (Nb ) when sample sizes were near or greater than the true Nb and when assuming a monogamous mating system. We further demonstrated that both methods provide roughly equivalent estimates of Nb . Our results indicate that sibship reconstruction and current SNP panels provide promise for estimating Nb in steelhead populations in the region.

Climate, Demography, and Zoogeography Predict Introgression Thresholds in Salmonid Hybrid Zones in Rocky Mountain Streams.

Among the many threats posed by invasions of nonnative species is introgressive hybridization, which can lead to the genomic extinction of native taxa. This phenomenon is regarded as common and perhaps inevitable among native cutthroat trout and introduced rainbow trout in western North America, despite that these taxa naturally co-occur in some locations. We conducted a synthetic analysis of 13,315 genotyped fish from 558 sites by building logistic regression models using data from geospatial stream databases and from 12 published studies of hybridization to assess whether environmental covariates could explain levels of introgression between westslope cutthroat trout and rainbow trout in the U.S. northern Rocky Mountains. A consensus model performed well (AUC, 0.78-0.86; classification success, 72-82%; 10-fold cross validation, 70-82%) and predicted that rainbow trout introgression was significantly associated with warmer water temperatures, larger streams, proximity to warmer habitats and to recent sources of rainbow trout propagules, presence within the historical range of rainbow trout, and locations further east. Assuming that water temperatures will continue to rise in response to climate change and that levels of introgression outside the historical range of rainbow trout will equilibrate with those inside that range, we applied six scenarios across a 55,234-km stream network that forecast 9.5-74.7% declines in the amount of habitat occupied by westslope cutthroat trout populations of conservation value, but not the wholesale loss of such populations. We conclude that introgression between these taxa is predictably related to environmental conditions, many of which can be manipulated to foster largely genetically intact populations of westslope cutthroat trout and help managers prioritize conservation activities.

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

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

A new species of myxozoan (Myxosporea) from the brain and spinal cord of rainbow trout (Oncorhynchus mykiss) from Idaho.

A new species of Myxosporea, Myxobolus neurotropus n. sp., is described from the brain and spinal cord of rainbow trout (Oncorhynchus mykiss) from Duncan Creek, Owyhee County, Idaho. Spores are oval, have 2 pyriform polar capsules, and possess a thick spore wall (sutural rim) with a short intracapsular offshoot. The mean spore dimensions are length 11.8 microm, width 10.8 microm, and thickness 8.8 microm. This myxozoan is compared to other described Myxobolus species found in cranial tissues of salmonids in terms of spore morphology and phylogenetic analysis. Because it is found in brain and spinal cord, it is encountered while performing screening tests for Myxobolus cerebralis, the causative agent of salmonid whirling disease. Where chronic inflammation and granulomatous lesions are associated with M. cerebralis, histological examination shows no host response to M. neurotropus n. sp. A diagnostic polymerase chain reaction (PCR) test is included as an aid in properly identifying the species.