Ecology and Chronic Wasting Disease Epidemiology Shape Prion Protein Gene Variation in Rocky Mountain Elk (Cervus elaphus nelsoni).

As chronic wasting disease (CWD) continues to spread across North America, the relationship between CWD and host genetics has become of interest. In Rocky Mountain elk (Cervus elaphus nelsoni), one or two copies of a leucine allele at codon 132 of the prion protein gene (132L*) has been shown to prolong the incubation period of CWD. Our study examined the relationship between CWD epidemiology and codon 132 evolution in elk from Wyoming, USA, from 2011 to 2018. Using PCR and Sanger sequencing, we genotyped 997 elk and assessed the relationship between genotype and CWD prevalence estimated from surveillance data. Using logistic regression, we showed that each 1% increase in CWD prevalence is associated with a 9.6% increase in the odds that an elk would have at least one copy of leucine at codon 132. In some regions, however, 132L* variants were found in the absence of CWD, indicating that evolutionary and epidemiologic patterns can be heterogeneous across space and time. We also provide evidence that naturally occurring CWD is not rare in 132L* elk, which merits the study of shedding kinetics in 132L* elk and the influence of genotype on CWD strain diversity. The management implications of cervid adaptations to CWD are difficult to predict. Studies that investigate the degree to which evolutionary outcomes are shaped by host spatial structure can provide useful epidemiologic insight, which can in turn aid management by informing scale and extent of mitigation actions.

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.

Multi-population puma connectivity could restore genomic diversity to at-risk coastal populations in California.

Urbanization is decreasing wildlife habitat and connectivity worldwide, including for apex predators, such as the puma (Puma concolor). Puma populations along California's central and southern coastal habitats have experienced rapid fragmentation from development, leading to calls for demographic and genetic management. To address urgent conservation genomic concerns, we used double-digest restriction-site associated DNA (ddRAD) sequencing to analyze 16,285 genome-wide single-nucleotide polymorphisms (SNPs) from 401 pumas sampled broadly across the state. Our analyses indicated support for 4-10 geographically nested, broad- to fine-scale genetic clusters. At the broadest scale, the four genetic clusters had high genetic diversity and exhibited low linkage disequilibrium, indicating that pumas have retained genomic diversity statewide. However, multiple lines of evidence indicated substructure, including 10 finer-scale genetic clusters, some of which exhibited fixed alleles and linkage disequilibrium. Fragmented populations along the Southern Coast and Central Coast had particularly low genetic diversity and strong linkage disequilibrium, indicating genetic drift and close inbreeding. Our results demonstrate that genetically at risk populations are typically nested within a broader-scale group of interconnected populations that collectively retain high genetic diversity and heterogenous fixations. Thus, extant variation at the broader scale has potential to restore diversity to local populations if management actions can enhance vital gene flow and recombine locally sequestered genetic diversity. These state- and genome-wide results are critically important for science-based conservation and management practices. Our nested population genomic analysis highlights the information that can be gained from population genomic studies aiming to provide guidance for the conservation of fragmented populations.

Spatio-temporal analyses reveal infectious disease-driven selection in a free-ranging ungulate.

Infectious diseases play an important role in wildlife population dynamics by altering individual fitness, but detecting disease-driven natural selection in free-ranging populations is difficult due to complex disease-host relationships. Chronic wasting disease (CWD) is a fatal infectious prion disease in cervids for which mutations in a single gene have been mechanistically linked to disease outcomes, providing a rare opportunity to study disease-driven selection in wildlife. In Wyoming, USA, CWD has gradually spread across mule deer (Odocoileus hemionus) populations, producing natural variation in disease history to evaluate selection pressure. We used spatial variation and a novel temporal comparison to investigate the relationship between CWD and a mutation at codon 225 of the mule deer prion protein gene that slows disease progression. We found that individuals with the 'slow' 225F allele were less likely to test positive for CWD, and the 225F allele was more common in herds exposed to CWD longer. We also found that in the past 2 decades, the 225F allele frequency increased more in herds with higher CWD prevalence. This study expanded on previous research by analysing spatio-temporal patterns of individual and herd-based disease data to present multiple lines of evidence for disease-driven selection in free-ranging wildlife.

Pronghorn population genomics show connectivity in the core of their range.

Preserving connectivity in the core of a species' range is crucial for long-term persistence. However, a combination of ecological characteristics, social behavior, and landscape features can reduce connectivity among wildlife populations and lead to genetic structure. Pronghorn (Antilocapra americana), for example, exhibit fluctuating herd dynamics and variable seasonal migration strategies, but GPS tracking studies show that landscape features such as highways impede their movements, leading to conflicting hypotheses about expected levels of genetic structure. Given that pronghorn populations declined significantly in the early 1900s, have only partially recovered, and are experiencing modern threats from landscape modification, conserving connectivity among populations is important for their long-term persistence in North America. To assess the genetic structure and diversity of pronghorn in the core of their range, we genotyped 4,949 genome-wide single-nucleotide polymorphisms and 11 microsatellites from 398 individuals throughout the state of Wyoming. We found no evidence of genetic subdivision and minimal evidence of isolation by distance despite a range that spans hundreds of kilometers, multiple mountain ranges, and three interstate highways. In addition, a rare variant analysis using putatively recent mutations found no genetic division between pronghorn on either side of a major highway corridor. Although we found no evidence that barriers to daily and seasonal movements of pronghorn impede gene flow, we suggest periodic monitoring of genetic structure and diversity as a part of management strategies to identify changes in connectivity.

Accuracy of de novo assembly of DNA sequences from double-digest libraries varies substantially among software.

Advances in DNA sequencing have made it feasible to gather genomic data for non-model organisms and large sets of individuals, often using methods for sequencing subsets of the genome. Several of these methods sequence DNA associated with endonuclease restriction sites (various RAD and GBS methods). For use in taxa without a reference genome, these methods rely on de novo assembly of fragments in the sequencing library. Many of the software options available for this application were originally developed for other assembly types and we do not know their accuracy for reduced representation libraries. To address this important knowledge gap, we simulated data from the Arabidopsis thaliana and Homo sapiens genomes and compared de novo assemblies by six software programs that are commonly used or promising for this purpose (ABySS, CD-HIT, Stacks, Stacks2, Velvet and VSEARCH). We simulated different mutation rates and types of mutations, and then applied the six assemblers to the simulated data sets, varying assembly parameters. We found substantial variation in software performance across simulations and parameter settings. ABySS failed to recover any true genome fragments, and Velvet and VSEARCH performed poorly for most simulations. Stacks and Stacks2 produced accurate assemblies of simulations containing SNPs, but the addition of insertion and deletion mutations decreased their performance. CD-HIT was the only assembler that consistently recovered a high proportion of true genome fragments. Here, we demonstrate the substantial difference in the accuracy of assemblies from different software programs and the importance of comparing assemblies that result from different parameter settings.