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.

Evidence of an Absence of Inbreeding Depression in a Wild Population of Weddell Seals (Leptonychotes weddellii).

Inbreeding depression can reduce the viability of wild populations. Detecting inbreeding depression in the wild is difficult; developing accurate estimates of inbreeding can be time and labor intensive. In this study, we used a two-step modeling procedure to incorporate uncertainty inherent in estimating individual inbreeding coefficients from multilocus genotypes into estimates of inbreeding depression in a population of Weddell seals (Leptonychotes weddellii). The two-step modeling procedure presented in this paper provides a method for estimating the magnitude of a known source of error, which is assumed absent in classic regression models, and incorporating this error into inferences about inbreeding depression. The method is essentially an errors-in-variables regression with non-normal errors in both the dependent and independent variables. These models, therefore, allow for a better evaluation of the uncertainty surrounding the biological importance of inbreeding depression in non-pedigreed wild populations. For this study we genotyped 154 adult female seals from the population in Erebus Bay, Antarctica, at 29 microsatellite loci, 12 of which are novel. We used a statistical evidence approach to inference rather than hypothesis testing because the discovery of both low and high levels of inbreeding are of scientific interest. We found evidence for an absence of inbreeding depression in lifetime reproductive success, adult survival, age at maturity, and the reproductive interval of female seals in this population.

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.

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.

Candidate adaptive genes associated with lineage divergence: identifying SNPs via next-generation targeted resequencing in mule deer (Odocoileus hemionus).

Mule deer (Odocoileus hemionus) are an excellent nonmodel species for empirically testing hypotheses in landscape and population genomics due to their large population sizes (low genetic drift), relatively continuous distribution, diversity of occupied habitats and phenotypic variation. Because few genomic resources are currently available for this species, we used exon data from a cattle (Bos taurus) reference genome to direct targeted resequencing of 5935 genes in mule deer. We sequenced approximately 3.75 Mbp at minimum 20X coverage in each of the seven mule deer, identifying 23 204 single nucleotide polymorphisms (SNPs) within, or adjacent to, 6886 exons in 3559 genes. We found 91 SNP loci (from 69 genes) with putatively fixed allele frequency differences between the two major lineages of mule deer (mule deer and black-tailed deer), and our estimate of mean genetic divergence (genome-wide FST  = 0.123) between these lineages was consistent with previous findings using microsatellite loci. We detected an over-representation of gamete generation and amino acid transport genes among the genes with SNPs exhibiting potentially fixed allele frequency differences between lineages. This targeted resequencing approach using exon capture techniques has identified a suite of loci that can be used in future research to investigate the genomic basis of adaptation and differentiation between black-tailed deer and mule deer. This study also highlights techniques (and an exon capture array) that will facilitate population genomic research in other cervids and nonmodel organisms.

The population history of endogenous retroviruses in mule deer (Odocoileus hemionus).

Mobile elements are powerful agents of genomic evolution and can be exceptionally informative markers for investigating species and population-level evolutionary history. While several studies have utilized retrotransposon-based insertional polymorphisms to resolve phylogenies, few population studies exist outside of humans. Endogenous retroviruses are LTR-retrotransposons derived from retroviruses that have become stably integrated in the host genome during past infections and transmitted vertically to subsequent generations. They offer valuable insight into host-virus co-evolution and a unique perspective on host evolutionary history because they integrate into the genome at a discrete point in time. We examined the evolutionary history of a cervid endogenous gammaretrovirus (CrERVγ) in mule deer (Odocoileus hemionus). We sequenced 14 CrERV proviruses (CrERV-in1 to -in14), and examined the prevalence and distribution of 13 proviruses in 262 deer among 15 populations from Montana, Wyoming, and Utah. CrERV absence in white-tailed deer (O. virginianus), identical 5' and 3' long terminal repeat (LTR) sequences, insertional polymorphism, and CrERV divergence time estimates indicated that most endogenization events occurred within the last 200000 years. Population structure inferred from CrERVs (F ST = 0.008) and microsatellites (θ = 0.01) was low, but significant, with Utah, northwestern Montana, and a Helena herd being particularly differentiated. Clustering analyses indicated regional structuring, and non-contiguous clustering could often be explained by known translocations. Cluster ensemble results indicated spatial localization of viruses, specifically in deer from northeastern and western Montana. This study demonstrates the utility of endogenous retroviruses to elucidate and provide novel insight into both ERV evolutionary history and the history of contemporary host populations.

Growth and long-term somatic and germline chimerism following fusion of juvenile Botryllus schlosseri.

The colonial ascidian Botryllus schlosseri undergoes a histocompatibility reaction that can result in vascular fusion of distinct genotypes, creating a chimera. Chimerism has both potential benefits, such as an immediate increase in size that may enhance growth rates, and costs. For the latter, the presence of multiple genotypes in a chimera can lead to competition between genetically distinct stem cell lineages, resulting in complete replacement of somatic and germline tissues by a single genotype. Although fusion can occur at any point after metamorphosis, previous studies have focused on chimeras created from sexually mature adults, where no benefit to chimerism has been documented. Here we focus on the costs and benefits of fusion between juveniles, characterizing growth rates and patterns of somatic and germline chimerism after natural and controlled fusion events. We also compared outcomes between low- and high-density growth conditions, the latter more likely representative of what occurs in natural populations. We found that growth rates were density-dependent, and that only chimeras grew under high-density conditions. We also observed a positional component to a post-fusion event called resorption, indicating that extrinsic factors were important in this process. Patterns of germline and somatic chimerism and dominance in chimeras made from fused juveniles were equivalent to those after fusion of sexually mature adults, and there were no age-related differences in these processes. Finally, by using genetic markers that could retrospectively assign genotypes, we also found that the majority of individual testes in a chimera were clonally derived.