A new GTSeq resource to facilitate multijurisdictional research and management of walleye Sander vitreus.

Conservation and management professionals often work across jurisdictional boundaries to identify broad ecological patterns. These collaborations help to protect populations whose distributions span political borders. One common limitation to multijurisdictional collaboration is consistency in data recording and reporting. This limitation can impact genetic research, which relies on data about specific markers in an organism's genome. Incomplete overlap of markers between separate studies can prevent direct comparisons of results. Standardized marker panels can reduce the impact of this issue and provide a common starting place for new research. Genotyping-in-thousands (GTSeq) is one approach used to create standardized marker panels for nonmodel organisms. Here, we describe the development, optimization, and early assessments of a new GTSeq panel for use with walleye (Sander vitreus) from the Great Lakes region of North America. High genome-coverage sequencing conducted using RAD capture provided genotypes for thousands of single nucleotide polymorphisms (SNPs). From these markers, SNP and microhaplotype markers were chosen, which were informative for genetic stock identification (GSI) and kinship analysis. The final GTSeq panel contained 500 markers, including 197 microhaplotypes and 303 SNPs. Leave-one-out GSI simulations indicated that GSI accuracy should be greater than 80% in most jurisdictions. The false-positive rates of parent-offspring and full-sibling kinship identification were found to be low. Finally, genotypes could be consistently scored among separate sequencing runs >94% of the time. Results indicate that the GTSeq panel that we developed should perform well for multijurisdictional walleye research throughout the Great Lakes region.

Systematics of the New World bats Eptesicus and Histiotus suggest trans-marine dispersal followed by Neotropical cryptic diversification.

Biodiversity can be boosted by colonization of new habitats such as remote islands and separated continents. Molecular studies have suggested that recently evolved organisms probably colonized already separated continents by dispersal, either via land bridge connections or crossing the ocean. Here we test the on-land and trans-marine dispersal hypotheses by evaluating possibilities of colonization routes over the Bering land bridge and across the Atlantic Ocean in the cosmopolitan bat genus Eptesicus (Chiroptera, Vespertilionidae). Previous molecular studies have found New World Eptesicus more closely related to Histiotus, a Neotropical endemic lineage with enlarged ears, than to Old World Eptesicus. However, phylogenetic relationships within the New World group remained unresolved and their evolutionary history was unclear. Here we studied the systematics of New World Eptesicus and Histiotus using extensive taxonomic and geographic sampling, and genomic data from thousands of ultra-conserved elements (UCEs). We estimated phylogenetic trees using concatenation and multispecies coalescent. All analyses supported four major New World clades and a novel topology where E. fuscus and Histiotus are sister clades that together diverged from two sister clades of Neotropical Eptesicus. Intra-clade divergence suggested cryptic diversity that has been concealed by morphological features, especially in the Neotropics where taxonomic re-evaluations are warranted. Molecular dating estimated that Old World and New World clades diverged around 17 million years ago followed by radiation of major New World clades in the mid-Miocene, when climatic changes might have facilitated global dispersal and radiation events. Biogeographic ancestral reconstruction supported the Neotropical origin of the New World clades, suggesting a trans-Atlantic colonization route from North Africa to the northern Neotropics. We highlight that trans-marine dispersal may be more prevalent than currently acknowledged and may be an important first step to global biodiversification.

Comparing genome-based estimates of relatedness for use in pedigree-based conservation management.

Researchers have long debated which estimator of relatedness best captures the degree of relationship between two individuals. In the genomics era, this debate continues, with relatedness estimates being sensitive to the methods used to generate markers, marker quality, and levels of diversity in sampled individuals. Here, we compare six commonly used genome-based relatedness estimators (kinship genetic distance [KGD], Wang maximum likelihood [TrioML], Queller and Goodnight [Rxy ], Kinship INference for Genome-wide association studies [KING-robust), and pairwise relatedness [RAB ], allele-sharing coancestry [AS]) across five species bred in captivity-including three birds and two mammals-with varying degrees of reliable pedigree data, using reduced-representation and whole genome resequencing data. Genome-based relatedness estimates varied widely across estimators, sequencing methods, and species, yet the most consistent results for known first order relationships were found using Rxy , RAB , and AS. However, AS was found to be less consistently correlated with known pedigree relatedness than either Rxy or RAB . Our combined results indicate there is not a single genome-based estimator that is ideal across different species and data types. To determine the most appropriate genome-based relatedness estimator for each new data set, we recommend assessing the relative: (1) correlation of candidate estimators with known relationships in the pedigree and (2) precision of candidate estimators with known first-order relationships. These recommendations are broadly applicable to conservation breeding programmes, particularly where genome-based estimates of relatedness can complement and complete poorly pedigreed populations. Given a growing interest in the application of wild pedigrees, our results are also applicable to in situ wildlife management.

Selection and demography drive range-wide patterns of MHC-DRB variation in mule deer.

BACKGROUND: Standing genetic variation is important especially in immune response-related genes because of threats to wild populations like the emergence of novel pathogens. Genetic variation at the major histocompatibility complex (MHC), which is crucial in activating the adaptive immune response, is influenced by both natural selection and historical population demography, and their relative roles can be difficult to disentangle. To provide insight into the influences of natural selection and demography on MHC evolution in large populations, we analyzed geographic patterns of variation at the MHC class II DRB exon 2 locus in mule deer (Odocoileus hemionus) using sequence data collected across their entire broad range. RESULTS: We identified 31 new MHC-DRB alleles which were phylogenetically similar to other cervid MHC alleles, and one allele that was shared with white-tailed deer (Odocoileus virginianus). We found evidence for selection on the MHC including high dN/dS ratios, positive neutrality tests, deviations from Hardy-Weinberg Equilibrium (HWE) and a stronger pattern of isolation-by-distance (IBD) than expected under neutrality. Historical demography also shaped variation at the MHC, as indicated by similar spatial patterns of variation between MHC and microsatellite loci and a lack of association between genetic variation at either locus type and environmental variables. CONCLUSIONS: Our results show that both natural selection and historical demography are important drivers in the evolution of the MHC in mule deer and work together to shape functional variation and the evolution of the adaptive immune response in large, well-connected populations.

The relevance of pedigrees in the conservation genomics era.

Over the past 50 years conservation genetics has developed a substantive toolbox to inform species management. One of the most long-standing tools available to manage genetics-the pedigree-has been widely used to characterize diversity and maximize evolutionary potential in threatened populations. Now, with the ability to use high throughput sequencing to estimate relatedness, inbreeding, and genome-wide functional diversity, some have asked whether it is warranted for conservation biologists to continue collecting and collating pedigrees for species management. In this perspective, we argue that pedigrees remain a relevant tool, and when combined with genomic data, create an invaluable resource for conservation genomic management. Genomic data can address pedigree pitfalls (e.g., founder relatedness, missing data, uncertainty), and in return robust pedigrees allow for more nuanced research design, including well-informed sampling strategies and quantitative analyses (e.g., heritability, linkage) to better inform genomic inquiry. We further contend that building and maintaining pedigrees provides an opportunity to strengthen trusted relationships among conservation researchers, practitioners, Indigenous Peoples, and Local Communities.

Nonrandom missing data can bias Principal Component Analysis inference of population genetic structure.

Population genetic studies in non-model systems increasingly use next-generation sequencing to obtain more loci, but such methods also generate more missing data that may affect downstream analyses. Here we focus on the principal component analysis (PCA) which has been widely used to explore and visualize population structure with mean-imputed missing data. We simulated data of different population models with various total missingness (1%, 10%, 20%) introduced either randomly or biased among individuals or populations. We found that individuals biased with missing data would be dragged away from their real population clusters to the origin of PCA plots, making them indistinguishable from true admixed individuals and potentially leading to misinterpreted population structure. We also generated empirical data of the big brown bat (Eptesicus fuscus) using restriction site-associated DNA sequencing (RADseq). We filtered three data sets with 19.12%, 9.87%, and 1.35% total missingness, all showing nonrandom missing data with biased individuals dragged towards the PCA origin, consistent with results from simulations. We highlight the importance of considering missing data effects on PCA in non-model systems where nonrandom missing data are common due to varying sample quality. To help detect missing data effects, we suggest to (1) plot PCA with a colour gradient showing per sample missingness, (2) interpret samples close to the PCA origin with extra caution, (3) explore filtering parameters with and without the missingness-biased samples, and (4) use complementary analyses (e.g., model-based methods) to cross-validate PCA results and help interpret population structure.

Characterizing patterns of genomic variation in the threatened Utah prairie dog: Implications for conservation and management.

Utah prairie dogs (Cynomys parvidens) are federally threatened due to eradication campaigns, habitat destruction, and outbreaks of plague. Today, Utah prairie dogs exist in small, isolated populations, making them less demographically stable and more susceptible to erosion of genetic variation by genetic drift. We characterized patterns of genetic structure at neutral and putatively adaptive loci in order to evaluate the relative effects of genetic drift and local adaptation on population divergence. We sampled individuals across the Utah prairie dog species range and generated 2955 single nucleotide polymorphisms using double digest restriction site-associated DNA sequencing. Genetic diversity was lower in low-elevation sites compared to high-elevation sites. Population divergence was high among sites and followed an isolation-by-distance model. Our results indicate that genetic drift plays a substantial role in the population divergence of the Utah prairie dog, and colonies would likely benefit from translocation of individuals between recovery units, which are characterized by distinct elevations, despite the detection of environmental associations with outlier loci. By understanding the processes that shape genetic structure, better informed decisions can be made with respect to the management of threatened species to ensure that adaptation is not stymied.

Phylogeography of Cedros and Tiburón Island Mule Deer in North America's Desert Southwest.

Though mule deer (Odocoileus hemionus) persist in robust populations throughout most of their North American distribution, habitat loss, unregulated hunting, and other factors have reduced their historical range in México. Two of the 6 putative subspecies inhabiting México's deserts and Baja California peninsula are of conservation concern, occupying islands in the Pacific Ocean (Odocoileus hemionus cerrosensis on Cedros Island: endangered) and Sea of Cortés (Odocoileus hemionus sheldoni on Tiburón Island: threatened). Focusing on the desert southwest (n = 448), we sampled Tiburón (n = 22) and Cedros (n = 15) Island mule deer using contemporary samples and natural history museum specimens to complete a phylogeographic evaluation of the species complex, and assess the phylogeography of these insular subspecies. Both insular subspecies formed endemic haplotype lineages, consistent with island biogeographic theory. Bayesian skyline plots were consistent with Holocene demographic expansion. Cedros Island deer were genetically most similar to adjacent mainland Baja California deer, but exhibited a suite of unique haplotypes and reduced genetic variation. Tiburón Island deer haplotypes unexpectedly nested within a mainland lineage found in distant New Mexico, rather than the adjacent mainland Sonoran lineage. Such findings suggest the importance of postglacial climate fluctuations and biotic community turnover in the phylogeographic history of mule deer in the desert southwest. Our genetic data corroborates cultural, archaeological, and phenotypic evidence supporting Cedros and Tiburón deer endemicity and subspecies status. Reduced genetic variation, divergence from mainland populations, and demographic trends on both islands indicate that conservation, monitoring, and management are critical to ensure persistence of these endemic insular subspecies.

Major histocompatibility complex variation is similar in little brown bats before and after white-nose syndrome outbreak.

White-nose syndrome (WNS), caused by the fungal pathogen Pseudogymnoascus destructans (Pd), has driven alarming declines in North American hibernating bats, such as little brown bat (Myotis lucifugus). During hibernation, infected little brown bats are able to initiate anti-Pd immune responses, indicating pathogen-mediated selection on the major histocompatibility complex (MHC) genes. However, such immune responses may not be protective as they interrupt torpor, elevate energy costs, and potentially lead to higher mortality rates. To assess whether WNS drives selection on MHC genes, we compared the MHC DRB gene in little brown bats pre- (Wisconsin) and post- (Michigan, New York, Vermont, and Pennsylvania) WNS (detection spanning 2014-2015). We genotyped 131 individuals and found 45 nucleotide alleles (27 amino acid alleles) indicating a maximum of 3 loci (1-5 alleles per individual). We observed high allelic admixture and a lack of genetic differentiation both among sampling sites and between pre- and post-WNS populations, indicating no signal of selection on MHC genes. However, post-WNS populations exhibited decreased allelic richness, reflecting effects from bottleneck and drift following rapid population declines. We propose that mechanisms other than adaptive immunity are more likely driving current persistence of little brown bats in affected regions.

Integrating genomics into conservation management.

Captive populations play a vital role in conservation, supporting both in situ and ex situ conservation projects. A healthy, self-sustaining captive population serves as valuable insurance against extinctions and a source of individuals for reintroduction into the wild. Genetic and genomic data can advance conservation and management across a range of scenarios. Although expanding, genomic data do not exist for most captive populations and there remain few examples of how these data can be integrated directly into conservation efforts. In this issue of Molecular Ecology Resources, Humble et al. (2020) present a chromosomal-level genome assembly for the extinct-in-the-wild scimitar-horned oryx (Oryx dammah), and resequencing data to survey genetic variation across captive populations. These data are exciting for reintroduction planning by offering guidance for selecting source populations to maximize global genetic variation and valuable resources for effective post-release monitoring. The genomic resources generated also are valuable tools with the power to yield novel and important insights into the evolution of scimitar-horned oryx and other antelope species.

Evidence for adaptive introgression of exons across a hybrid swarm in deer.

BACKGROUND: Secondary contact between closely related lineages can result in a variety of outcomes, including hybridization, depending upon the strength of reproductive barriers. By examining the extent to which different parts of the genome introgress, it is possible to infer the strength of selection and gain insight into the evolutionary trajectory of lineages. Following secondary contact approximately 8000 years ago in the Pacific Northwest, mule deer (Odocoileus hemionus hemionus) and black-tailed deer (O. h. columbianus) formed a hybrid swarm along the Cascade mountain range despite substantial differences in body size (up to two times) and habitat preference. In this study, we examined genetic population structure, extent of introgression, and selection pressures in freely interbreeding populations of mule deer and black-tailed deer using mitochondrial DNA sequences, 9 microsatellite loci, and 95 SNPs from protein-coding genes. RESULTS: We observed bi-directional hybridization and classified approximately one third of the 172 individuals as hybrids, almost all of which were beyond the F1 generation. High genetic differentiation between black-tailed deer and mule deer at protein-coding genes suggests that there is positive divergent selection, though selection on these loci is relatively weak. Contrary to predictions, there was not greater selection on protein-coding genes thought to be associated with immune function and mate choice. Geographic cline analyses were consistent across genetic markers, suggesting long-term stability (over hundreds of generations), and indicated that the center of the hybrid swarm is 20-30 km to the east of the Cascades ridgeline, where there is a steep ecological transition from wet, forested habitat to dry, scrub habitat. CONCLUSIONS: Our data are consistent with a genetic boundary between mule deer and black-tailed deer that is porous but maintained by many loci under weak selection having a substantial cumulative effect. The absence of clear reproductive barriers and the consistent centering of geographic clines at a sharp ecotone suggests that ecology is a driver of hybrid swarm dynamics. Adaptive introgression in this study (and others) promotes gene flow and provides valuable insight into selection strength on specific genes and the evolutionary trajectory of hybridizing taxa.

Guidelines for planning genomic assessment and monitoring of locally adaptive variation to inform species conservation.

Identifying and monitoring locally adaptive genetic variation can have direct utility for conserving species at risk, especially when management may include actions such as translocations for restoration, genetic rescue, or assisted gene flow. However, genomic studies of local adaptation require careful planning to be successful, and in some cases may not be a worthwhile use of resources. Here, we offer an adaptive management framework to help conservation biologists and managers decide when genomics is likely to be effective in detecting local adaptation, and how to plan assessment and monitoring of adaptive variation to address conservation objectives. Studies of adaptive variation using genomic tools will inform conservation actions in many cases, including applications such as assisted gene flow and identifying conservation units. In others, assessing genetic diversity, inbreeding, and demographics using selectively neutral genetic markers may be most useful. And in some cases, local adaptation may be assessed more efficiently using alternative approaches such as common garden experiments. Here, we identify key considerations of genomics studies of locally adaptive variation, provide a road map for successful collaborations with genomics experts including key issues for study design and data analysis, and offer guidelines for interpreting and using results from genomic assessments to inform monitoring programs and conservation actions.

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.

Hidden in plain sight: Cryptic and endemic malaria parasites in North American white-tailed deer (Odocoileus virginianus).

Malaria parasites of the genus Plasmodium are diverse in mammal hosts, infecting five mammalian orders in the Old World, but were long considered absent from the diverse deer family (Cervidae) and from New World mammals. There was a description of a Plasmodium parasite infecting a single splenectomized white-tailed deer (WTD; Odocoileus virginianus) in 1967 but none have been reported since, which has proven a challenge to our understanding of malaria parasite biogeography. Using both microscopy and polymerase chain reaction, we screened a large sample of native and captive ungulate species from across the United States for malaria parasites. We found a surprisingly high prevalence (up to 25%) and extremely low parasitemia of Plasmodium parasites in WTD throughout the eastern United States. We did not detect infections in the other ungulate species nor in western WTD. We also isolated the parasites from the mosquito Anopheles punctipennis. Morphologically, the parasites resemble the parasite described in 1967, Plasmodium odocoilei. Our analysis of the cytochrome b gene revealed two divergent Plasmodium clades in WTD representative of species that likely diverged 2.3 to 6 million years ago, concurrent with the arrival of the WTD ancestor into North America across Beringia. Multigene phylogenetic analysis placed these clades within the larger malaria parasite clade. We document Plasmodium parasites to be common in WTD, endemic to the New World, and as the only known malaria parasites from deer (Cervidae). These findings reshape our knowledge of the phylogeography of the malaria parasites and suggest that other mammal taxa may harbor infection by endemic and occult malaria parasites.

Assessing the permeability of landscape features to animal movement: using genetic structure to infer functional connectivity.

Human-altered environments often challenge native species with a complex spatial distribution of resources. Hostile landscape features can inhibit animal movement (i.e., genetic exchange), while other landscape attributes facilitate gene flow. The genetic attributes of organisms inhabiting such complex environments can reveal the legacy of their movements through the landscape. Thus, by evaluating landscape attributes within the context of genetic connectivity of organisms within the landscape, we can elucidate how a species has coped with the enhanced complexity of human altered environments. In this research, we utilized genetic data from eastern chipmunks (Tamias striatus) in conjunction with spatially explicit habitat attribute data to evaluate the realized permeability of various landscape elements in a fragmented agricultural ecosystem. To accomplish this we 1) used logistic regression to evaluate whether land cover attributes were most often associated with the matrix between or habitat within genetically identified populations across the landscape, and 2) utilized spatially explicit habitat attribute data to predict genetically-derived Bayesian probabilities of population membership of individual chipmunks in an agricultural ecosystem. Consistency between the results of the two approaches with regard to facilitators and inhibitors of gene flow in the landscape indicate that this is a promising new way to utilize both landscape and genetic data to gain a deeper understanding of human-altered ecosystems.

Range-wide analysis of genetic structure in a widespread, highly mobile species (Odocoileus hemionus) reveals the importance of historical biogeography.

Highly mobile species that thrive in a wide range of habitats are expected to show little genetic differentiation across their range. A limited but growing number of studies have revealed that patterns of broad-scale genetic differentiation can and do emerge in vagile, continuously distributed species. However, these patterns are complex and often shaped by both historical and ecological factors. Comprehensive surveys of genetic variation at a broad scale and at high resolution are useful for detecting cryptic spatial genetic structure and for investigating the relative roles of historical and ecological processes in structuring widespread, highly mobile species. In this study, we analysed 10 microsatellite loci from over 1900 samples collected across the full range of mule deer (Odocoileus hemionus), one of the most widely distributed and abundant of all large mammal species in North America. Through both individual- and population-based analyses, we found evidence for three main genetic lineages, one corresponding to the 'mule deer' morphological type and two to the 'black-tailed deer' type. Historical biogeographic events likely are the primary drivers of genetic divergence in this species; boundaries of the three lineages correspond well with predictions based on Pleistocene glacial cycles, and substructure within each lineage demonstrates island vicariance. However, across large geographic areas, including the entire mule deer lineage, we found that genetic variation fit an isolation-by-distance pattern rather than discrete clusters. A lack of genetic structure across wide geographic areas of the continental west indicates that ecological processes have not resulted in restrictions to gene flow sufficient for spatial genetic structure to emerge. Our results have important implications for our understanding of evolutionary mechanisms of divergence, as well as for taxonomy, conservation and management.

Identification of novel single nucleotide polymorphisms (SNPs) in deer (Odocoileus spp.) using the BovineSNP50 BeadChip.

Single nucleotide polymorphisms (SNPs) are growing in popularity as a genetic marker for investigating evolutionary processes. A panel of SNPs is often developed by comparing large quantities of DNA sequence data across multiple individuals to identify polymorphic sites. For non-model species, this is particularly difficult, as performing the necessary large-scale genomic sequencing often exceeds the resources available for the project. In this study, we trial the Bovine SNP50 BeadChip developed in cattle (Bos taurus) for identifying polymorphic SNPs in cervids Odocoileus hemionus (mule deer and black-tailed deer) and O. virginianus (white-tailed deer) in the Pacific Northwest. We found that 38.7% of loci could be genotyped, of which 5% (n = 1068) were polymorphic. Of these 1068 polymorphic SNPs, a mixture of putatively neutral loci (n = 878) and loci under selection (n = 190) were identified with the F(ST)-outlier method. A range of population genetic analyses were implemented using these SNPs and a panel of 10 microsatellite loci. The three types of deer could readily be distinguished with both the SNP and microsatellite datasets. This study demonstrates that commercially developed SNP chips are a viable means of SNP discovery for non-model organisms, even when used between very distantly related species (the Bovidae and Cervidae families diverged some 25.1-30.1 million years before present).

A suite of microsatellite markers optimized for amplification of DNA from Addax (Addax nasomaculatus) blood preserved on FTA cards.

The addax (Addax nasomaculatus) is a critically endangered antelope that is currently maintained in zoos through regional, conservation breeding programs. As for many captive species, incomplete pedigree data currently impedes the ability of addax breeding programs to confidently manage the genetics of captive populations and to select appropriate animals for reintroduction. Molecular markers are often used to improve pedigree resolution, thereby improving the long-term effectiveness of genetic management. When developing a suite of molecular markers, it is important to consider the source of DNA, as the utility of markers may vary across DNA sources. In this study, we optimized a suite of microsatellite markers for use in genotyping captive addax blood samples collected on FTA cards. We amplified 66 microsatellite loci previously described in other Artiodactyls. Sixteen markers amplified a single product in addax, but only 5 of these were found to be polymorphic in a sample of 37 addax sampled from a captive herd at Fossil Rim Wildlife Center in the US. The suite of microsatellite markers developed in this study provides a new tool for the genetic management of captive addax, and demonstrates that FTA cards can be a useful means of sample storage, provided appropriate loci are used in downstream analyses.

Fine-scale analysis reveals cryptic landscape genetic structure in desert tortoises.

Characterizing the effects of landscape features on genetic variation is essential for understanding how landscapes shape patterns of gene flow and spatial genetic structure of populations. Most landscape genetics studies have focused on patterns of gene flow at a regional scale. However, the genetic structure of populations at a local scale may be influenced by a unique suite of landscape variables that have little bearing on connectivity patterns observed at broader spatial scales. We investigated fine-scale spatial patterns of genetic variation and gene flow in relation to features of the landscape in desert tortoise (Gopherus agassizii), using 859 tortoises genotyped at 16 microsatellite loci with associated data on geographic location, sex, elevation, slope, and soil type, and spatial relationship to putative barriers (power lines, roads). We used spatially explicit and non-explicit Bayesian clustering algorithms to partition the sample into discrete clusters, and characterize the relationships between genetic distance and ecological variables to identify factors with the greatest influence on gene flow at a local scale. Desert tortoises exhibit weak genetic structure at a local scale, and we identified two subpopulations across the study area. Although genetic differentiation between the subpopulations was low, our landscape genetic analysis identified both natural (slope) and anthropogenic (roads) landscape variables that have significantly influenced gene flow within this local population. We show that desert tortoise movements at a local scale are influenced by features of the landscape, and that these features are different than those that influence gene flow at larger scales. Our findings are important for desert tortoise conservation and management, particularly in light of recent translocation efforts in the region. More generally, our results indicate that recent landscape changes can affect gene flow at a local scale and that their effects can be detected almost immediately.

Hybrid swarm between divergent lineages of mule deer (Odocoileus hemionus).

Studies of hybrid zones have revealed an array of evolutionary outcomes, yet the underlying structure is typically characterized as one of three types: a hybrid zone, a hybrid swarm or a hybrid taxon. Our primary objective was to determine which of these three structures best characterizes a zone of hybridization between two divergent lineages of mule deer (Odocoileus hemionus), mule deer and black-tailed deer. These lineages are morphologically, ecologically and genetically distinct, yet hybridize readily along a zone of secondary contact between the east and west slopes of the Cascade Mountains (Washington and Oregon, USA). Using microsatellite and mitochondrial DNA, we found clear evidence for extensive hybridization and introgression between lineages, with varying degrees of admixture across the zone of contact. The pattern of hybridization in this region closely resembles a hybrid swarm; based on data from 10 microsatellite loci, we detected hybrids that extend well beyond the F1 generation, did not detect linkage disequilibrium at the centre of the zone and found that genotypes were associated randomly within the zone of contact. Introgression was characterized as bidirectional and symmetric, which is surprising given that the zone of contact occurs along a sharp ecotone and that lineages are characterized by large differences in body size (a key component of mating success). Regardless of the underlying mechanisms promoting hybrid swarm maintenance, it is clear that the persistence of a hybrid swarm presents unique challenges for management in this region.