The genetic legacy of 50 years of desert bighorn sheep translocations

Selection leads to remarkable variability in the outcomes of hybridisation across replicate hybrid zones

Hybrid zones have been viewed as an opportunity to see speciation in action. When hybrid zones are replicated, it is assumed that if the same genetic incompatibilities are maintaining reproductive isolation across all instances of secondary contact, those incompatibilities should be identifiable by consistent patterns in the genome. In contrast, changes in allele frequencies due to genetic drift should be idiosyncratic for each hybrid zone. To test this assumption, we simulated 20 replicates of each of 12 hybrid zone scenarios with varied genetic incompatibilities, rates of migration, selection and different starting population size ratios of parental species. We found remarkable variability in the outcomes of hybridisation in replicate hybrid zones, particularly with Bateson-Dobzhansky-Muller incompatibilities and strong selection. We found substantial differences among replicates in the overall genomic composition of individuals, including admixture proportions, inter-specific ancestry complement and number of ancestry junctions. Additionally, we found substantial variation in genomic clines among replicates at focal loci, regardless of locus-specific selection. We conclude that processes other than selection are responsible for some consistent outcomes of hybridisation, whereas selection on incompatibilities can lead to genomically widespread and highly variable outcomes. We highlight the challenge of mapping between pattern and process in hybrid zones and call attention to how selection against incompatibilities will commonly lead to variable outcomes. We hope that this study informs future research on replicate hybrid zones and encourages further development of statistical techniques, theoretical models and exploration of additional axes of variation to understand reproductive isolation.

The potential influence of genome-wide adaptive divergence on conservation translocation outcome in an isolated greater sage-grouse population

Conservation translocations are an important conservation tool commonly employed to augment declining or reestablish extirpated populations. One goal of augmentation is to increase genetic diversity and reduce the risk of inbreeding depression (i.e., genetic rescue). However, introducing individuals from significantly diverged populations risks disrupting coadapted traits and reducing local fitness (i.e., outbreeding depression). Genetic data are increasingly more accessible for wildlife species and can provide unique insight regarding the presence and retention of introduced genetic variation from augmentation as an indicator of effectiveness and adaptive similarity as an indicator of source and recipient population suitability. We used 2 genetic data sets to evaluate augmentation of isolated populations of greater sage-grouse (Centrocercus urophasianus) in the northwestern region of the species range (Washington, USA) and to retrospectively evaluate adaptive divergence among source and recipient populations. We developed 2 statistical models for microsatellite data to evaluate augmentation outcomes. We used one model to predict genetic diversity after augmentation and compared these predictions with observations of genetic change. We used the second model to quantify the amount of observed reproduction attributed to transplants (proof of population integration). We also characterized genome-wide adaptive divergence among source and recipient populations. Observed genetic diversity (HO = 0.65) was higher in the recipient population than predicted had no augmentation occurred (HO = 0.58) but less than what was predicted by our model (HO = 0.75). The amount of shared genetic variation between the 2 geographically isolated resident populations increased, which is evidence of periodic gene flow previously assumed to be rare. Among candidate adaptive genes associated with elevated fixation index (FST) (143 genes) or local environmental variables (97 and 157 genes for each genotype-environment association method, respectively), we found clusters of genes with related functions that may influence the ability of transplants to use local resources and navigate unfamiliar environments and their reproductive potential, all possible reasons for low genetic retention from augmentation.

Restoration of bighorn sheep: History, successes, and remaining conservation issues

Mammals are imperiled worldwide, primarily from habitat loss or modification, and exhibit downward trends in their populations and distributions. Likewise, large-bodied herbivores have undergone a collapse in numbers and are at the highest extinction risk of all mammals. Bighorn sheep (Ovis canadensis) are among those large-bodied herbivores that possess a slow-paced life history, suffer from debilitating diseases, and have experienced range contractions across their historical distribution since the late 1800s. Translocations and reintroductions of these mountain ungulates are key aspects of restoration and often are used to re-establish populations in historical habitat or to supplement declining herds. Millions of US dollars and much effort by state and federal natural resource agencies, as well as public and private organizations, have been expended to restore bighorn sheep. Despite those efforts, translocated populations of bighorn sheep have not always been successful. We assessed restoration of bighorn sheep to provide insights in the context of conservation of populations of bighorn sheep, because this management tool is a frequently used to re-establish populations. We focused briefly on past efforts to restore bighorn sheep populations and followed with updates on the value of habitat enhancements, genetic issues, the importance of ecotypic or phenotypic adaptations when restoring populations, predation, and disease transmission. We also raised issues and posed questions that have potential to affect future decisions regarding the restoration of bighorn sheep. This information will help conservationists improve the success of conserving these iconic large mammals.

Population connectivity patterns of genetic diversity, immune responses and exposure to infectious pneumonia in a metapopulation of desert bighorn sheep

Habitat fragmentation is an important driver of biodiversity loss and can be remediated through management actions aimed at maintenance of natural connectivity in metapopulations. Connectivity may protect populations from infectious diseases by preserving immunogenetic diversity and disease resistance. However, connectivity could exacerbate the risk of infectious disease spread across vulnerable populations. We tracked the spread of a novel strain of Mycoplasma ovipneumoniae in a metapopulation of desert bighorn sheep Ovis canadensis nelsoni in the Mojave Desert to investigate how variation in connectivity among populations influenced disease outcomes. M. ovipneumoniae was detected throughout the metapopulation, indicating that the relative isolation of many of these populations did not protect them from pathogen invasion. However, we show that connectivity among bighorn sheep populations was correlated with higher immunogenetic diversity, a protective immune response and lower disease prevalence. Variation in protective immunity predicted infection risk in individual bighorn sheep and was associated with heterozygosity at genetic loci linked to adaptive and innate immune signalling. Together, these findings may indicate that population connectivity maintains immunogenetic diversity in bighorn sheep populations in this system and has direct effects on immune responses in individual bighorn sheep and their susceptibility to infection by a deadly pathogen. Our study suggests that the genetic benefits of population connectivity could outweigh the risk of infectious disease spread and supports conservation management that maintains natural connectivity in metapopulations.

Genomic and common garden approaches yield complementary results for quantifying environmental drivers of local adaptation in rubber rabbitbrush, a foundational Great Basin shrub

The spatial structure of genomic and phenotypic variation across populations reflects historical and demographic processes as well as evolution via natural selection. Characterizing such variation can provide an important perspective for understanding the evolutionary consequences of changing climate and for guiding ecological restoration. While evidence for local adaptation has been traditionally evaluated using phenotypic data, modern methods for generating and analyzing landscape genomic data can directly quantify local adaptation by associating allelic variation with environmental variation. Here, we analyze both genomic and phenotypic variation of rubber rabbitbrush (Ericameria nauseosa), a foundational shrub species of western North America. To quantify landscape genomic structure and provide perspective on patterns of local adaptation, we generated reduced representation sequencing data for 17 wild populations (222 individuals; 38,615 loci) spanning a range of environmental conditions. Population genetic analyses illustrated pronounced landscape genomic structure jointly shaped by geography and environment. Genetic-environment association (GEA) analyses using both redundancy analysis (RDA) and a machine-learning approach (Gradient Forest) indicated environmental variables (precipitation seasonality, slope, aspect, elevation, and annual precipitation) influenced spatial genomic structure and were correlated with allele frequency shifts indicative of local adaptation at a consistent set of genomic regions. We compared our GEA-based inference of local adaptation with phenotypic data collected by growing seeds from each population in a greenhouse common garden. Population differentiation in seed weight, emergence, and seedling traits was associated with environmental variables (e.g., precipitation seasonality) that were also implicated in GEA analyses, suggesting complementary conclusions about the drivers of local adaptation across different methods and data sources. Our results provide a baseline understanding of spatial genomic structure for E. nauseosa across the western Great Basin and illustrate the utility of GEA analyses for detecting the environmental causes and genetic signatures of local adaptation in a widely distributed plant species of restoration significance.

Spatial population genetics in heavily managed species: Separating patterns of historical translocation from contemporary gene flow in white-tailed deer

Approximately 100 years ago, unregulated harvest nearly eliminated white-tailed deer (Odocoileus virginianus) from eastern North America, which subsequently served to catalyze wildlife management as a national priority. An extensive stock-replenishment effort soon followed, with deer broadly translocated among states as a means of re-establishment. However, an unintended consequence was that natural patterns of gene flow became obscured and pretranslocation signatures of population structure were replaced. We applied cutting-edge molecular and biogeographic tools to disentangle genetic signatures of historical management from those reflecting spatially heterogeneous dispersal by evaluating 35,099 single nucleotide polymorphisms (SNPs) derived via reduced-representation genomic sequencing from 1143 deer sampled statewide in Arkansas. We then employed Simpson's diversity index to summarize ancestry assignments and visualize spatial genetic transitions. Using sub-sampled transects across these transitions, we tested clinal patterns across loci against theoretical expectations of their response under scenarios of re-colonization and restricted dispersal. Two salient results emerged: (A) Genetic signatures from historic translocations are demonstrably apparent; and (B) Geographic filters (major rivers; urban centers; highways) now act as inflection points for the distribution of this contemporary ancestry. These results yielded a statewide assessment of contemporary population structure in deer as driven by historic translocations as well as ongoing processes. In addition, the analytical framework employed herein to effectively decipher extant/historic drivers of deer distribution in Arkansas is also applicable for other biodiversity elements with similarly complex demographic histories.

Subspecies hybridization as a potential conservation tool in species reintroductions

Reintroductions are a powerful tool for the recovery of endangered species. However, their long-term success is strongly influenced by the genetic diversity of the reintroduced population. The chances of population persistence can be improved by enhancing the population's adaptive ability through the mixing of individuals from different sources. However, where source populations are too diverse the reintroduced population could also suffer from outbreeding depression or unsuccessful admixture due to behavioural or genetic barriers. For the reintroduction of Asiatic wild ass Equus hemionus ssp. in Israel, a breeding core was created from individuals of two different subspecies (E. h. onager & E. h. kulan). Today the population comprises approximately 300 individuals and displays no signs of outbreeding depression. The aim of this study was a population genomic evaluation of this conservation reintroduction protocol. We used maximum likelihood methods and genetic clustering analyses to investigate subspecies admixture and test for spatial autocorrelation based on subspecies ancestry. Further, we analysed heterozygosity and effective population sizes in the breeding core prior to release and the current wild population. We discovered high levels of subspecies admixture in the breeding core and wild population, consistent with a significant heterozygote excess in the breeding core. Furthermore, we found no signs of spatial autocorrelation associated with subspecies ancestry in the wild population. Inbreeding and variance effective population size estimates were low. Our results indicate no genetic or behavioural barriers to admixture between the subspecies and suggest that their hybridization has led to greater genetic diversity in the reintroduced population. The study provides rare empirical evidence of the successful application of subspecies hybridization in a reintroduction. It supports use of intraspecific hybridization as a tool to increase genetic diversity in conservation translocations.

Evaluating wildlife translocations using genomics: A bighorn sheep case study

Wildlife restoration often involves translocation efforts to reintroduce species and supplement small, fragmented populations. We examined the genomic consequences of bighorn sheep (Ovis canadensis) translocations and population isolation to enhance understanding of evolutionary processes that affect population genetics and inform future restoration strategies. We conducted a population genomic analysis of 511 bighorn sheep from 17 areas, including native and reintroduced populations that received 0-10 translocations. Using the Illumina High Density Ovine array, we generated datasets of 6,155 to 33,289 single nucleotide polymorphisms and completed clustering, population tree, and kinship analyses. Our analyses determined that natural gene flow did not occur between most populations, including two pairs of native herds that had past connectivity. We synthesized genomic evidence across analyses to evaluate 24 different translocation events and detected eight successful reintroductions (i.e., lack of signal for recolonization from nearby populations) and five successful augmentations (i.e., reproductive success of translocated individuals) based on genetic similarity with the source populations. A single native population founded six of the reintroduced herds, suggesting that environmental conditions did not need to match for populations to persist following reintroduction. Augmentations consisting of 18-57 animals including males and females succeeded, whereas augmentations of two males did not result in a detectable genetic signature. Our results provide insight on genomic distinctiveness of native and reintroduced herds, information on the relative success of reintroduction and augmentation efforts and their associated attributes, and guidance to enhance genetic contribution of augmentations and reintroductions to aid in bighorn sheep restoration.

Genetic assignment of illegally trafficked neotropical primates and implications for reintroduction programs

The black and gold howler monkey (Alouatta caraya) is a neotropical primate threatened by habitat loss and capture for illegal trade in Argentina. Using multilocus microsatellite genotypes from 178 A. caraya individuals sampled from 15 localities in Argentina, we built a genotype reference database (GRDB). Bayesian assignment methods applied to the GRDB allowed us to correctly re-assign 73% of individuals to their true location of origin and 93.3% to their cluster of origin. We used the GRDB to assign 22 confiscated individuals (17 of which were reintroduced), and 3 corpses to both localities and clusters of origin. We assigned with a probability >70% the locality of origin of 14 individuals and the cluster of origin of 21. We found that most of the confiscated individuals were assigned to one cluster (F-Ch-C) and two localities included in the GRDB, suggesting that trafficked A. caraya primarily originated in this area. Our results reveal that only 4 of 17 reintroduced individuals were released in sites corresponding to their cluster of origin. Our findings illustrate the applicability of genotype databases for inferring hotspots of illegal capture and for guiding future reintroduction efforts, both of which are essential elements of species protection and recovery programs.