USGS ScienceSearch

Geology topics

Sara J. Oyler-McCance

Publications and source records attributed to Sara J. Oyler-McCance.

At least 55 records · Page 3Linked to original sources

Evaluation of a Chicken 600K SNP genotyping array in non-model species of grouse

The use of single nucleotide polymorphism (SNP) arrays to generate large SNP datasets for comparison purposes have recently become an attractive alternative to other genotyping methods. Although most SNP arrays were originally developed for domestic organisms, they can be effectively applied to wild relatives to obtain large panels of SNPs. In this study, we tested the cross-species application of the Affymetrix 600K Chicken SNP array in five species of North American prairie grouse ( Centrocercus and Tympanuchus genera). Two individuals were genotyped per species for a total of ten samples. A high proportion (91%) of the total 580 961 SNPs were genotyped in at least one individual (73–76% SNPs genotyped per species). Principal component analysis with autosomal SNPs separated the two genera, but failed to clearly distinguish species within genera. Gene ontology analysis identified a set of genes related to morphogenesis and development (including genes involved in feather development), which may be primarily responsible for large phenotypic differences between Centrocercus and Tympanuchus grouse. Our study provided evidence for successful cross-species application of the chicken SNP array in grouse which diverged ca. 37 mya from the chicken lineage. As far as we are aware, this is the first reported application of a SNP array in non-passerine birds, and it demonstrates the feasibility of using commercial SNP arrays in research on non-model bird species.

Scientific Reports

Evaluation of genetic change from translocation among Gunnison Sage-Grouse (Centrocercus minimus) populations

Maintenance of genetic diversity is important for conserving species, especially those with fragmented habitats or ranges. In the absence of natural dispersal, translocation can be used to achieve this goal, although the success of translocation can be difficult to measure. Here we evaluate genetic change following translocation in Gunnison Sage-Grouse ( Centrocercus minimus ), a species reduced to 7 discrete populations with low levels of gene flow and high levels of genetic differentiation. Between 2000 and 2014, 306 birds from the largest and most genetically diverse population (Gunnison Basin) were translocated to 5 much smaller satellite populations to augment local population size and increase genetic diversity. Although the magnitude of the effect varied by population, we found evidence of increased genetic variation, decreased genetic differentiation from Gunnison Basin, and reproduction between translocated individuals and resident birds. These results suggest that translocations are impacting satellite populations, with current data providing a new baseline for genetic diversity among populations of this imperiled species.

Colorado, Utah

An experimental comparison of composite and grab sampling of stream water for metagenetic analysis of environmental DNA

Use of environmental DNA (eDNA) to assess distributions of aquatic and semi-aquatic macroorganisms is promising, but sampling schemes may need to be tailored to specific objectives. Given the potentially high variance in aquatic eDNA among replicate grab samples, compositing smaller water volumes collected over a period of time may be more effective for some applications. In this study, we compared eDNA profiles from composite water samples aggregated over three hours with grab water samples. Both sampling patterns were performed with identical autosamplers paired at two different sites in a headwater stream environment, augmented with exogenous fish eDNA from an upstream rearing facility. Samples were filtered through 0.8 μm cellulose nitrate filters and DNA was extracted with a cetyl trimethylammonium bromide procedure. Eukaryotic and bacterial community profiles were derived by amplicon sequencing of 12S ribosomal, 16S ribosomal, and cytochrome oxidase I loci. Operational taxa were assigned to genus with a lowest common ancestor approach for eukaryotes and to family with the RDP Classifier software for prokaryotes. Eukaryotic community profiles were more consistent with composite sampling than grab sampling. Downstream, rarefaction curves suggested faster taxon accumulation for composite samples, and estimated richness was higher for composite samples as a set than for grab samples. Upstream, composite sampling produced lower estimated richness than grab samples, but with overlapping standard errors. Furthermore, a bimodal pattern of richness as a function of sequence counts suggested the impact of clumped particles on upstream samples. Bacterial profiles were insensitive to sample method, consistent with the more even dispersion expected for bacteria compared with eukaryotic eDNA. Overall, samples composited over 3 h performed equal to or better than triplicate grab sampling for quantitative community metrics, despite the higher total sequencing effort provided to grab replicates. On the other hand, taxon-specific detection rates did not differ appreciably and the two methods gave similar estimates of the ratio of the common fish genera Salmo and Coregonus at each site. Unexpectedly, Salmo eDNA dropped out substantially faster than Coregonus eDNA between the two sites regardless of sampling method, suggesting that differential settling affects the estimation of relative abundance. We identified bacterial patterns that were associated with eukaryotic diversity, suggesting potential roles as biomarkers of sample representativeness.

PeerJ

Characterizing range-wide population divergence in an alpine-endemic bird: A comparison of genetic and genomic approaches

The delineation of intraspecific units that are evolutionarily and demographically distinct is an important step in the development of species-specific management plans. Neutral genetic variation has served as the primary data source for delineating “evolutionarily significant units,” but with recent advances in genomic technology, we now have an unprecedented ability to utilize information about neutral and adaptive variation across the entire genome. Here, we use traditional genetic markers (microsatellites) and a newer reduced-representation genomic approach (single nucleotide polymorphisms) to delineate distinct groups of white-tailed ptarmigan (Lagopus leucura), an alpine-obligate species that is distributed in naturally fragmented habitats from Alaska to New Mexico. Five subspecies of white-tailed ptarmigan are currently recognized but their distinctiveness has not been verified with molecular data. Based on analyses of 436 samples at 12 microsatellite loci and 95 samples at 14,866 single nucleotide polymorphism loci, we provide strong support for treating two subspecies as distinct intraspecific units—L. l. altipetens, found in Colorado and neighboring states; and L. l. saxatilis, found on British Columbia’s Vancouver Island—but our findings reveal more moderate patterns of divergence within the remainder of the species’ range. Results based on genetic and genomic datasets generally agreed with one another, indicating that in many cases microsatellite loci may be sufficient for describing major patterns of genetic structure across species’ ranges. This work will inform future conservation and management decisions for the white-tailed ptarmigan, a species that may be vulnerable to future changes in climate.

Conservation Genetics

Quantifying functional connectivity: The role of breeding habitat, abundance, and landscape features on range‐wide gene flow in sage‐grouse

Functional connectivity, quantified using landscape genetics, can inform conservation through the identification of factors linking genetic structure to landscape mechanisms. We used breeding habitat metrics, landscape attributes, and indices of grouse abundance, to compare fit between structural connectivity and genetic differentiation within five long‐established Sage‐Grouse Management Zones (MZ) I‐V using microsatellite genotypes from 6,844 greater sage‐grouse ( Centrocercus urophasianus ) collected across their 10.7 million‐km 2 range. We estimated structural connectivity using a circuit theory‐based approach where we built resistance surfaces using thresholds dividing the landscape into “habitat” and “nonhabitat” and nodes were clusters of sage‐grouse leks (where feather samples were collected using noninvasive techniques). As hypothesized, MZ‐specific habitat metrics were the best predictors of differentiation. To our surprise, inclusion of grouse abundance‐corrected indices did not greatly improve model fit in most MZs. Functional connectivity of breeding habitat was reduced when probability of lek occurrence dropped below 0.25 (MZs I, IV) and 0.5 (II), thresholds lower than those previously identified as required for the formation of breeding leks, which suggests that individuals are willing to travel through undesirable habitat. The individual MZ landscape results suggested terrain roughness and steepness shaped functional connectivity across all MZs. Across respective MZs, sagebrush availability (<10%–30%; II, IV, V), tree canopy cover (>10%; I, II, IV), and cultivation (>25%; I, II, IV, V) each reduced movement beyond their respective thresholds. Model validations confirmed variation in predictive ability across MZs with top resistance surfaces better predicting gene flow than geographic distance alone, especially in cases of low and high differentiation among lek groups. The resultant resistance maps we produced spatially depict the strength and redundancy of range‐wide gene flow and can help direct conservation actions to maintain and restore functional connectivity for sage‐grouse.

Evolutionary Applications

The genetic network of greater sage-grouse: Range-wide identification of keystone hubs of connectivity

Genetic networks can characterize complex genetic relationships among groups of individuals, which can be used to rank nodes most important to the overall connectivity of the system. Ranking allows scarce resources to be guided toward nodes integral to connectivity. The greater sage-grouse (Centrocercus urophasianus) is a species of conservation concern that breeds on spatially discrete leks that must remain connected by genetic exchange for population persistence. We genotyped 5,950 individuals from 1,200 greater sage-grouse leks distributed across the entire species’ geographic range. We found a small-world network composed of 458 nodes connected by 14,481 edges. This network was composed of hubs—that is, nodes facilitating gene flow across the network—and spokes—that is, nodes where connectivity is served by hubs. It is within these hubs that the greatest genetic diversity was housed. Using indices of network centrality, we identified hub nodes of greatest conservation importance. We also identified keystone nodes with elevated centrality despite low local population size. Hub and keystone nodes were found across the entire species’ contiguous range, although nodes with elevated importance to network-wide connectivity were found more central: especially in northeastern, central, and southwestern Wyoming and eastern Idaho. Nodes among which genes are most readily exchanged were mostly located in Montana and northern Wyoming, as well as Utah and eastern Nevada. The loss of hub or keystone nodes could lead to the disintegration of the network into smaller, isolated subnetworks. Protecting both hub nodes and keystone nodes will conserve genetic diversity and should maintain network connections to ensure a resilient and viable population over time. Our analysis shows that network models can be used to model gene flow, offering insights into its pattern and process, with application to prioritizing landscapes for conservation.

Idaho, Wyoming, Montana, Utah, Nevada

Long-term persistence of horse fecal DNA in the environment makes equids particularly good candidates for non-invasive sampling

Fecal DNA collected noninvasively can provide valuable information about genetic and ecological characteristics. This approach has rarely been used for equids, despite the need for conservation of endangered species and management of abundant feral populations. We examined factors affecting the efficacy of using equid fecal samples for conservation genetics. First, we evaluated two fecal collection methods (paper bag vs. ethanol). Then, we investigated how time since deposition and month of collection impacted microsatellite amplification success and genotyping errors. Between May and November 2014, we collected feral horse fecal samples of known age each month in a feral horse Herd Management Area in western Colorado and documented deterioration in the field with photographs. Samples collected and dried in paper bags had significantly higher amplification rates than those collected and stored in ethanol. There was little difference in the number of loci that amplified per sample between fresh fecal piles and those that had been exposed to the environment for up to 2 months (in samples collected in paper bags). After 2 months of exposure, amplification success declined. When comparing fresh (0–2 months) and old (3–6 months) fecal piles, samples from fresh piles had more matching genotypes across samples, better amplification success and less allelic dropout. Samples defecated during the summer and collected within 2 months of deposition had highest number of genotypes matching among samples, and lowest rates of amplification failure and allelic dropout. Due to the digestive system and amount of fecal material produced by equids, as well as their occurrence in arid ecosystems, we suggest that they are particularly good candidates for noninvasive sampling using fecal DNA.

Colorado

Post-release breeding of translocated sharp-tailed grouse and an absence of artificial insemination effects

Context: Translocation has become a widely used method to restore wildlife populations following extirpation. For some species, such as lekking grouse, which breed at traditional mating grounds, reproduction is linked to culturally established geographic locations. Cultural centres are lost upon extirpation, making restoration into otherwise rehabilitated habitats especially challenging. The process by which species with culturally dependent reproduction sometimes become re-established is poorly understood and merits investigation to improve conservation strategies. Aims: We reintroduced CSTG to vacant habitat in north-central Nevada, USA, from 2013 to 2017, with concordant goals of promoting females to nest and males to lek. We tested the utility of performing artificial insemination (AI) on females before translocation and we conducted paternity analyses to understand male reproduction. Methods: We monitored females for the effects of AI on nest initiation, nest survival and egg fertility. We used post-hatch extra-embryonic membranes and other tissues to evaluate paternity of chicks produced at the restoration site. Key results: Artificial insemination had no effect on female survival or nest initiation, and did not fertilise any eggs within nine sampled clutches ( n = 102 eggs). Most paternity was attributable to male residents that had survived for ≥1 year at the restoration site before the arrival of translocated females. Conclusions: Artificial insemination neither aided nor harmed female reproduction. A small number of translocated, resident adult males sired reproduction following female release. Implications: The presence of resident males at restoration sites may be more likely to result in post-translocation reproduction than is pre-translocation AI. Restoring CSTG to vacant habitat should focus on translocating females into suitable nesting habitat while simultaneously ensuring that reproductively capable males are available within or adjacent to the nesting habitat.

Wildlife Research

Genetic mark–recapture improves estimates of maternity colony size for Indiana bats

Genetic mark–recapture methods are increasingly being used to estimate demographic parameters in species where traditional techniques are problematic or imprecise. The federally endangered Indiana bat Myotis sodalis has declined dramatically and threats such as white-nose syndrome continue to afflict this species. To date, important demographic information for Indiana bats has been difficult to estimate precisely using traditional techniques such as emergence counts. Successful management and protection of Indiana bats requires better methods to estimate population sizes and survival rates throughout the year, particularly during summer when these bats reproduce and are widely dispersed away from their winter hibernacula. In addition, the familial makeup of maternity colonies is unknown, yet important for understanding local and regional population dynamics. We had four objectives in this study. For the first two objectives we investigated the potential use of DNA from fecal samples (fecal DNA) collected at roosts to obtain genetically based mark–recapture estimates of 1) colony size and 2) survival rates, for an Indiana bat maternity colony in Indianapolis, Indiana. The third objective was to compare our genetically based colony-size estimates with emergence counts conducted at the same roost tree to evaluate the genetic mark–recapture method. Our fourth objective was to use fecal DNA to estimate levels of relatedness among individuals sampled at the roost. In the summer of 2008, we collected fecal pellets and conducted emergence counts at a prominent roost tree during three time periods each lasting 7 or 8 d. We genotyped fecal DNA using five highly polymorphic microsatellite loci to identify individuals and used a robust-design mark–recapture approach to estimate survival rates as well as colony size at the roost tree. Emergence count estimates at the roost tree ranged from 100 to 215, whereas genetic mark–recapture estimates were higher, ranging from 122 to 266 and more precise. Apparent survival was 0.994 (SE = 0.04) between sampling periods suggesting that few bats died or permanently emigrated during the course of the study. Relatedness estimates, r, between all pairs of individuals averaged 0.055 ranging from 0 to 0.779, indicating that most individuals were not closely related. We demonstrate here the promise of using fecal DNA to estimate demographic information for Indiana bats and potentially other bat species.

Journal of Fish and Wildlife Management

Extensive shared polymorphism at non-MHC immune genes in recently diverged North American prairie grouse

Gene polymorphisms shared between recently diverged species are thought to be widespread and most commonly reflect introgression from hybridization or retention of ancestral polymorphism through incomplete lineage sorting. Shared genetic diversity resulting from incomplete lineage sorting is usually maintained for a relatively short period of time, but under strong balancing selection it may persist for millions of years beyond species divergence (balanced trans-species polymorphism), as in the case of the major histocompatibility complex (MHC) genes. However, balancing selection is much less likely to act on non-MHC immune genes. The aim of this study was to investigate the patterns of shared polymorphism and selection at non-MHC immune genes in five grouse species from Centrocercus and Tympanuchus genera. For this purpose, we genotyped five non-MHC immune genes that do not interact directly with pathogens, but are involved in signaling and regulate immune cell growth. In contrast to previous studies with MHC, we found no evidence for balancing selection or balanced trans-species polymorphism among the non-MHC immune genes. No haplotypes were shared between genera and in most cases more similar allelic variants sorted by genus. Between species within genera, however, we found extensive shared polymorphism, which was most likely attributable to introgression or incomplete lineage sorting following recent divergence and large ancestral effective population size (i.e., weak genetic drift). Our study suggests that North American prairie grouse may have attained relatively low degree of reciprocal monophyly at nuclear loci and reinforces the rarity of balancing selection in non-MHC immune genes.

Immunogenetics

Developing approaches for linear mixed modeling in landscape genetics through landscape-directed dispersal simulations

Dispersal can impact population dynamics and geographic variation, and thus, genetic approaches that can establish which landscape factors influence population connectivity have ecological and evolutionary importance. Mixed models that account for the error structure of pairwise datasets are increasingly used to compare models relating genetic differentiation to pairwise measures of landscape resistance. A model selection framework based on information criteria metrics or explained variance may help disentangle the ecological and landscape factors influencing genetic structure, yet there are currently no consensus for the best protocols. Here, we develop landscape-directed simulations and test a series of replicates that emulate independent empirical datasets of two species with different life history characteristics (greater sage-grouse; eastern foxsnake). We determined that in our simulated scenarios, AIC and BIC were the best model selection indices and that marginal R 2 values were biased toward more complex models. The model coefficients for landscape variables generally reflected the underlying dispersal model with confidence intervals that did not overlap with zero across the entire model set. When we controlled for geographic distance, variables not in the underlying dispersal models (i.e., nontrue) typically overlapped zero. Our study helps establish methods for using linear mixed models to identify the features underlying patterns of dispersal across a variety of landscapes.

Ecology and Evolution

Polygamy slows down population divergence in shorebirds

Sexual selection may act as a promotor of speciation since divergent mate choice and competition for mates can rapidly lead to reproductive isolation. Alternatively, sexual selection may also retard speciation since polygamous individuals can access additional mates by increased breeding dispersal. High breeding dispersal should hence increase gene flow and reduce diversification in polygamous species. Here, we test how polygamy predicts diversification in shorebirds using genetic differentiation and subspecies richness as proxies for population divergence. Examining microsatellite data from 79 populations in 10 plover species (Genus: Charadrius ) we found that polygamous species display significantly less genetic structure and weaker isolation-by-distance effects than monogamous species. Consistent with this result, a comparative analysis including 136 shorebird species showed significantly fewer subspecies for polygamous than for monogamous species. By contrast, migratory behavior neither predicted genetic differentiation nor subspecies richness. Taken together, our results suggest that dispersal associated with polygamy may facilitate gene flow and limit population divergence. Therefore, intense sexual selection, as occurs in polygamous species, may act as a brake rather than an engine of speciation in shorebirds. We discuss alternative explanations for these results and call for further studies to understand the relationships between sexual selection, dispersal, and diversification.

Evolution

Sage grouse

Sage grouse are a group of chicken-sized birds with a unique breeding behavior and dependence on sagebrush shrubs (genus Artemisia) for food and shelter throughout their life cycle. In the last century, human population expansion throughout western North America has reduced the amount of sagebrush and degraded and fragmented the remaining areas. Vanishing sagebrush has resulted in sage grouse (genus Centrocercus) population declines and elevated conservation concern. Western Colorado is home to both species of sage grouse: greater sage grouse (Centrocercus urophasianus) and Gunnison sage grouse (Centrocercus minimus). Populations in the state, and throughout their range, have declined sufficiently to warrant consideration for federal protection for both species under the Endangered Species Act.

Book chapter

Integration of genetic and demographic data to assess population risk in a continuously distributed species

The identification and demographic assessment of biologically meaningful populations is fundamental to species’ ecology and management. Although genetic tools are used frequently to identify populations, studies often do not incorporate demographic data to understand their respective population trends. We used genetic data to define subpopulations in a continuously distributed species. We assessed demographic independence and variation in population trends across the distribution. Additionally, we identified potential barriers to gene flow among subpopulations. We sampled greater sage-grouse ( Centrocercus urophasianus ) leks from across their range (≈175,000 Km 2 ) in Wyoming and amplified DNA at 14 microsatellite loci for 1761 samples. Subsequently, we assessed population structure in unrelated individuals ( n = 872) by integrating results from multiple Bayesian clustering approaches and used the boundaries to inform our assessment of long-term population trends and lek activity over the period of 1995–2013. We identified four genetic clusters of which two northern ones showed demographic independence from the others. Trends in population size for the northwest subpopulation were statistically different from the other three genetic clusters and the northeast and southwest subpopulations demonstrated a general trend of increasing proportion of inactive leks over time. Population change from 1996 to 2012 suggested population growth in the southern subpopulations and decline, or neutral, change in the northern subpopulations. We suggest that sage-grouse subpopulations in northern Wyoming are at greater risk of extirpation than the southern subpopulations due to smaller census and effective population sizes and higher variability within subpopulations. Our research is an example of incorporating genetic and demographic data and provides guidance on the identification of subpopulations of conservation concern.

Conservation Genetics

Differential influences of local subpopulations on regional diversity and differentiation for greater sage-grouse ( Centrocercus urophasianus )

The distribution of spatial genetic variation across a region can shape evolutionary dynamics and impact population persistence. Local population dynamics and among-population dispersal rates are strong drivers of this spatial genetic variation, yet for many species we lack a clear understanding of how these population processes interact in space to shape within-species genetic variation. Here, we used extensive genetic and demographic data from 10 subpopulations of greater sage-grouse to parameterize a simulated approximate Bayesian computation (ABC) model and (i) test for regional differences in population density and dispersal rates for greater sage-grouse subpopulations in Wyoming, and (ii) quantify how these differences impact subpopulation regional influence on genetic variation. We found a close match between observed and simulated data under our parameterized model and strong variation in density and dispersal rates across Wyoming. Sensitivity analyses suggested that changes in dispersal (via landscape resistance) had a greater influence on regional differentiation, whereas changes in density had a greater influence on mean diversity across all subpopulations. Local subpopulations, however, varied in their regional influence on genetic variation. Decreases in the size and dispersal rates of central populations with low overall and net immigration (i.e. population sources) had the greatest negative impact on genetic variation. Overall, our results provide insight into the interactions among demography, dispersal and genetic variation and highlight the potential of ABC to disentangle the complexity of regional population dynamics and project the genetic impact of changing conditions.

Molecular Ecology

Critical considerations for the application of environmental DNA methods to detect aquatic species

Species detection using environmental DNA (eDNA) has tremendous potential for contributing to the understanding of the ecology and conservation of aquatic species. Detecting species using eDNA methods, rather than directly sampling the organisms, can reduce impacts on sensitive species and increase the power of field surveys for rare and elusive species. The sensitivity of eDNA methods, however, requires a heightened awareness and attention to quality assurance and quality control protocols. Additionally, the interpretation of eDNA data demands careful consideration of multiple factors. As eDNA methods have grown in application, diverse approaches have been implemented to address these issues. With interest in eDNA continuing to expand, supportive guidelines for undertaking eDNA studies are greatly needed. Environmental DNA researchers from around the world have collaborated to produce this set of guidelines and considerations for implementing eDNA methods to detect aquatic macroorganisms. Critical considerations for study design include preventing contamination in the field and the laboratory, choosing appropriate sample analysis methods, validating assays, testing for sample inhibition and following minimum reporting guidelines. Critical considerations for inference include temporal and spatial processes, limits of correlation of eDNA with abundance, uncertainty of positive and negative results, and potential sources of allochthonous DNA. We present a synthesis of knowledge at this stage for application of this new and powerful detection method.

Methods in Ecology and Evolution

Contrasting evolutionary histories of MHC class I and class II loci in grouse—Effects of selection and gene conversion

Genes of the major histocompatibility complex (MHC) encode receptor molecules that are responsible for recognition of intracellular and extracellular pathogens (class I and class II genes, respectively) in vertebrates. Given the different roles of class I and II MHC genes, one might expect the strength of selection to differ between these two classes. Different selective pressures may also promote different rates of gene conversion at each class. Despite these predictions, surprisingly few studies have looked at differences between class I and II genes in terms of both selection and gene conversion. Here, we investigated the molecular evolution of MHC class I and II genes in five closely related species of prairie grouse ( Centrocercus and Tympanuchus ) that possess one class I and two class II loci. We found striking differences in the strength of balancing selection acting on MHC class I versus class II genes. More than half of the putative antigen-binding sites (ABS) of class II were under positive or episodic diversifying selection, compared with only 10 % at class I. We also found that gene conversion had a stronger role in shaping the evolution of MHC class II than class I. Overall, the combination of strong positive (balancing) selection and frequent gene conversion has maintained higher diversity of MHC class II than class I in prairie grouse. This is one of the first studies clearly demonstrating that macroevolutionary mechanisms can act differently on genes involved in the immune response against intracellular and extracellular pathogens.

Heredity

Rangewide genetic analysis of Lesser Prairie-Chicken reveals population structure, range expansion, and possible introgression

The distribution of the Lesser Prairie-Chicken ( Tympanuchus pallidicinctus ) has been markedly reduced due to loss and fragmentation of habitat. Portions of the historical range, however, have been recolonized and even expanded due to planting of conservation reserve program (CRP) fields that provide favorable vegetation structure for Lesser Prairie-Chickens. The source population(s) feeding the range expansion is unknown, yet has resulted in overlap between Lesser and Greater Prairie-Chickens ( T. cupido ) increasing the potential for hybridization. Our objectives were to characterize connectivity and genetic diversity among populations, identify source population(s) of recent range expansion, and examine hybridization with the Greater Prairie-Chicken. We analyzed 640 samples from across the range using 13 microsatellites. We identified three to four populations corresponding largely to ecoregions. The Shinnery Oak Prairie and Sand Sagebrush Prairie represented genetically distinct populations ( F ST > 0.034 and F ST > 0.023 respectively). The Shortgrass/CRP Mosaic and Mixed Grass ecoregions appeared admixed ( F ST = 0.009). Genetic diversity was similar among ecoregions and N e ranged from 142 (95 % CI 99&ndash;236) for the Shortgrass/CRP Mosaic to 296 (95 % CI 233&ndash;396) in the Mixed Grass Prairie. No recent migration was detected among ecoregions, except asymmetric dispersal from both the Mixed Grass Prairie and to a lesser extent the Sand Sagebrush Prairie north into adjacent Shortgrass/CRP Mosaic ( m = 0.207, 95 % CI 0.116&ndash;0.298, m = 0.097, 95 % CI 0.010&ndash;0.183, respectively). Indices investigating potential hybridization in the Shortgrass/CRP Mosaic revealed that six of the 13 individuals with hybrid phenotypes were significantly admixed suggesting hybridization. Continued monitoring of diversity within and among ecoregions is warranted as are actions promoting genetic connectivity and range expansion.

Colorado, Kansas, New Mexico, Texas