Population genomic diversity and structure at the discontinuous southern range of the Great Gray Owl in North America

Species' distributions are often discontinuous near the edge of the range where the environment may be more variable than the core of the range. Range discontinuity can reduce or cut off gene flow to small peripheral populations and lead to genetic drift and subsequent loss of genetic diversity. The southern extent of the Great Gray Owl (Strix nebulosa) range in North America is discontinuous, unlike their northern core range across the boreal forests. We sampled owls from five different locations on the periphery of the range across the western US (Wyoming, Idaho, California, northern Oregon, and southern Oregon) to investigate genetic population structure and genetic diversity. Using a reduced-representation genomic sequencing approach to genotype 123 individuals at 4817 single nucleotide polymorphic loci, we identified four genetically differentiated populations: California, southern Oregon, northern Oregon, and Wyoming and Idaho grouped together as a single Rocky Mountain population. The four genetically differentiated populations of Great Gray Owls identified in this study display high differentiation and low genetic variation, which is suggestive of long-term isolation and lack of connectivity, potentially caused by range discontinuity. The populations that lack habitat connectivity to the rest of the breeding range (i.e. those in California and Oregon) had lower genetic diversity than the Rocky Mountain population that is connected to the core of the range. These factors and other risks (such as disease and human-caused mortality) heighten susceptibility of these range-edge populations to future habitat and climate changes, genetic diversity erosion, and potential extinction vortex. For these reasons, protecting and monitoring this species on the southern edge of their range is vital

Assessing captive spawning strategies for supplementation production of Delta Smelt

Objective: To support the declining wild population of Delta Smelt, a conservation hatchery has expanded its mission from maintaining a backup population as insurance against extinction to also producing fish for release into the wild. The substantially higher production demands require a balance between producing large numbers of fish while adhering to conservation genetic principles that maximize retention of effective population size (Ne) and thus overall diversity. Methods: We performed spawning experiments at the hatchery to evaluate the genetic consequences of two spawning strategies: (1) a pooled strategy where we fertilized premixed eggs from three dams with premixed milt from three sires and (2) a partial-factorial strategy where eggs from three dams were mixed and then apportioned among three containers, each container then receiving milt from one sire. We used genetic parentage analysis of larval offspring to determine the reproductive success of spawners in 10 replicate crosses of each strategy. Result: The contributions of parents to offspring were more even in partial-factorial crosses and consequently resulted in higher Ne (average Ne = 5.50 ± 0.38; expected Ne = 6.0), suggesting its potential for maintaining genetic diversity over time. In contrast, our pooled spawning experiment produced lower and more variable Ne values (average Ne = 3.86 ± 1.30), demonstrating that this more efficient method of production entails high costs in terms of long-term genetic management. Treating our experiments as hypothetical pools of fish for release, we combined the Ne values for pooled or partial-factorial crosses to calculate the effective size of a release population (NeR). Unequal family sizes reduced NeR for our pooled experiment to half of the expected value, whereas the partial-factorial experiment NeR was 88% of the expected value. Conclusion: We discuss the benefits and risks of each method and how these can be considered when designing a spawning strategy for Delta Smelt supplementation

Novel hybrid finds a peri-urban niche: Allen’s Hummingbirds in southern California

Species range expansions and contractions can have ecological and genetic consequences, and thus are important areas of study for conservation. Hybridization and introgression are not uncommon in closely related populations that experience secondary contact during a range expansion. Allen’s Hummingbird (Selasphorus sasin) in California comprises two subspecies: the migratory S. s. sasin, which winters in central Mexico and breeds in central and northern California, and the resident S. s. sedentarius, which lives and breeds year-round on several of the Channel Islands off the California coast. Within recent decades, Allen’s Hummingbirds have been found living and breeding year-round in the southern California peri-urban mainland near Los Angeles. Ornithologists assumed that the L.A. birds were an expansion of the island subspecies, S. s. sedentarius due to similar but very subtle morphological characteristics. However, the genetic relationships among the three putative populations of Allen's hummingbird—migratory, southern California mainland, and island—are unknown. We investigated these relationships by analyzing variation of single nucleotide polymorphisms from the three geographic regions where S. sasin are present. Our population genomic analyses indicate that S. sasin hummingbirds inhabiting mainland southern California are a hybrid population resulting from admixture between S. s. sasin and S. s. sedentarius. From one perspective, these results may be interpreted as a positive development for S. s. sasin as the growing population represent an overall increase in the S. sasin population, and the expanding population contains a significant representation of S. s. sasin alleles

Novel hybrid finds a peri-urban niche: Allen’s Hummingbirds in southern California

Species range expansions and contractions can have ecological and genetic consequences, and thus are important areas of study for conservation. Hybridization and introgression are not uncommon in closely related populations that experience secondary contact during a range expansion. Allen’s Hummingbird (Selasphorus sasin) in California comprises two subspecies: the migratory S. s. sasin, which winters in central Mexico and breeds in central and northern California, and the resident S. s. sedentarius, which lives and breeds year-round on several of the Channel Islands off the California coast. Within recent decades, Allen’s Hummingbirds have been found living and breeding year-round in the southern California peri-urban mainland near Los Angeles. Ornithologists assumed that the L.A. birds were an expansion of the island subspecies, S. s. sedentarius due to similar but very subtle morphological characteristics. However, the genetic relationships among the three putative populations of Allen's hummingbird—migratory, southern California mainland, and island—are unknown. We investigated these relationships by analyzing variation of single nucleotide polymorphisms from the three geographic regions where S. sasin are present. Our population genomic analyses indicate that S. sasin hummingbirds inhabiting mainland southern California are a hybrid population resulting from admixture between S. s. sasin and S. s. sedentarius. From one perspective, these results may be interpreted as a positive development for S. s. sasin as the growing population represent an overall increase in the S. sasin population, and the expanding population contains a significant representation of S. s. sasin alleles

Functional connectivity in a continuously distributed, migratory species as revealed by landscape genomics

Maintaining functional connectivity is critical for the long-term conservation of wildlife populations. Landscape genomics provides an opportunity to assess long-term functional connectivity by relating environmental variables to spatial patterns of genomic variation resulting from generations of movement, dispersal and mating behaviors. Identifying landscape features associated with gene flow at large geographic scales for highly mobile species is becoming increasingly possible due to more accessible genomic approaches, improved analytical methods and enhanced computational power. We characterized the genetic structure and diversity of migratory mule deer Odocoileus hemionus using 4051 single nucleotide polymorphisms in 406 individuals sampled across multiple habitats throughout Wyoming, USA. We then identified environmental variables associated with genomic variation within genetic groups and statewide using a stepwise approach to first evaluate nonlinear relationships of landscape resistance with genetic distances and then use mixed-effects modeling to choose top landscape genomic models. We identified three admixed genetic groups of mule deer and found that environmental variables associated with gene flow varied among genetic groups, revealing scale-dependent and regional variation in functional connectivity. At the statewide scale, more gene flow occurred in areas with low elevation and mixed habitat. In the southern genetic group, more gene flow occurred in areas with low elevation. In the northern genetic group, more gene flow occurred in grassland and forest habitats, while highways and energy infrastructure reduced gene flow. In the western genetic group, the null model of isolation by distance best represented genetic patterns. Overall, our findings highlight the role of different seasonal ranges on mule deer genetic connectivity, and show that anthropogenic features hinder connectivity. This study demonstrates the value of combining a large, genome-wide marker set with recent advances in landscape genomics to evaluate functional connectivity in a wide-ranging migratory species

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

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