Search USGS⌕ Search

Geology topics

Sara J. Oyler-McCance

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

At least 73 records · Page 4Linked to original sources
USGS2016

Latent spatial models and sampling design for landscape genetics

We propose a spatially-explicit approach for modeling genetic variation across space and illustrate how this approach can be used to optimize spatial prediction and sampling design for landscape genetic data. We propose a multinomial data model for categorical microsatellite allele data commonly used in landscape genetic studies, and introduce a latent spatial random effect to allow for spatial correlation between genetic observations. We illustrate how modern dimension reduction approaches to spatial statistics can allow for efficient computation in landscape genetic statistical models covering large spatial domains. We apply our approach to propose a retrospective spatial sampling design for greater sage-grouse ( Centrocercus urophasianus ) population genetics in the western United States.

Annals of Applied Statistics↗
USGS2016

A field ornithologist’s guide to genomics: Practical considerations for ecology and conservation

Vast improvements in sequencing technology have made it practical to simultaneously sequence millions of nucleotides distributed across the genome, opening the door for genomic studies in virtually any species. Ornithological research stands to benefit in three substantial ways. First, genomic methods enhance our ability to parse and simultaneously analyze both neutral and non-neutral genomic regions, thus providing insight into adaptive evolution and divergence. Second, the sheer quantity of sequence data generated by current sequencing platforms allows increased precision and resolution in analyses. Third, high-throughput sequencing can benefit applications that focus on a small number of loci that are otherwise prohibitively expensive, time-consuming, and technically difficult using traditional sequencing methods. These advances have improved our ability to understand evolutionary processes like speciation and local adaptation, but they also offer many practical applications in the fields of population ecology, migration tracking, conservation planning, diet analyses, and disease ecology. This review provides a guide for field ornithologists interested in incorporating genomic approaches into their research program, with an emphasis on techniques related to ecology and conservation. We present a general overview of contemporary genomic approaches and methods, as well as important considerations when selecting a genomic technique. We also discuss research questions that are likely to benefit from utilizing high-throughput sequencing instruments, highlighting select examples from recent avian studies.

The Auk↗
USGS2015

Z chromosome divergence, polymorphism and relative effective population size in a genus of lekking birds

Sex chromosomes contribute disproportionately to species boundaries as they diverge faster than autosomes and often have reduced diversity. Their hemizygous nature contributes to faster divergence and reduced diversity, as do some types of selection. In birds, other factors (mating system and bottlenecks) can further decrease the effective population size of Z-linked loci and accelerate divergence (Fast-Z). We assessed Z-linked divergence and effective population sizes for two polygynous sage-grouse species and compared them to estimates from birds with various mating systems. We found lower diversity and higher F ST for Z-linked loci than for autosomes, as expected. The π Z / π A ratio was 0.38 in Centrocercus minimus, 0.48 in Centrocercus urophasianus and 0.59 in a diverged, parapatric population of C. urophasianus, a broad range given the mating system among these groups is presumably equivalent. The full data set had unequal males and females across groups, so we compared an equally balanced reduced set of C. minimus and individuals pooled from both C. urophasianus subgroups recovering similar estimates: 0.54 for C. urophasianus and 0.38 for C. minimus. We provide further evidence that N eZ / N eA in birds is often lower than expected under random mating or monogamy. The lower ratio in C. minimus could be a consequence of stronger selection or drift acting on Z loci during speciation, as this species differs strongly from C. urophasianus in sexually selected characters with minimal mitochondrial divergence. As C. minimus also exhibited lower genomic diversity, it is possible that a more severe demographic history may contribute to its lower ratio.

Heredity↗
USGS2015

Environmental DNA (eDNA) sampling improves occurrence and detection estimates of invasive Burmese pythons

Environmental DNA (eDNA) methods are used to detect DNA that is shed into the aquatic environment by cryptic or low density species. Applied in eDNA studies, occupancy models can be used to estimate occurrence and detection probabilities and thereby account for imperfect detection. However, occupancy terminology has been applied inconsistently in eDNA studies, and many have calculated occurrence probabilities while not considering the effects of imperfect detection. Low detection of invasive giant constrictors using visual surveys and traps has hampered the estimation of occupancy and detection estimates needed for population management in southern Florida, USA. Giant constrictor snakes pose a threat to native species and the ecological restoration of the Florida Everglades. To assist with detection, we developed species-specific eDNA assays using quantitative PCR (qPCR) for the Burmese python ( Python molurus bivittatus ), Northern African python ( P . sebae ), boa constrictor ( Boa constrictor ), and the green ( Eunectes murinus ) and yellow anaconda ( E . notaeus ). Burmese pythons, Northern African pythons, and boa constrictors are established and reproducing, while the green and yellow anaconda have the potential to become established. We validated the python and boa constrictor assays using laboratory trials and tested all species in 21 field locations distributed in eight southern Florida regions. Burmese python eDNA was detected in 37 of 63 field sampling events; however, the other species were not detected. Although eDNA was heterogeneously distributed in the environment, occupancy models were able to provide the first estimates of detection probabilities, which were greater than 91%. Burmese python eDNA was detected along the leading northern edge of the known population boundary. The development of informative detection tools and eDNA occupancy models can improve conservation efforts in southern Florida and support more extensive studies of invasive constrictors. Generic sampling design and terminology are proposed to standardize and clarify interpretations of eDNA-based occupancy models.

Florida↗
USGS2015

Landscape characteristics influencing the genetic structure of greater sage-grouse within the stronghold of their range: a holistic modeling approach

Given the significance of animal dispersal to population dynamics and geographic variability, understanding how dispersal is impacted by landscape patterns has major ecological and conservation importance. Speaking to the importance of dispersal, the use of linear mixed models to compare genetic differentiation with pairwise resistance derived from landscape resistance surfaces has presented new opportunities to disentangle the menagerie of factors behind effective dispersal across a given landscape. Here, we combine these approaches with novel resistance surface parameterization to determine how the distribution of high- and low-quality seasonal habitat and individual landscape components shape patterns of gene flow for the greater sage-grouse ( Centrocercus urophasianus ) across Wyoming. We found that pairwise resistance derived from the distribution of low-quality nesting and winter, but not summer, seasonal habitat had the strongest correlation with genetic differentiation. Although the patterns were not as strong as with habitat distribution, multivariate models with sagebrush cover and landscape ruggedness or forest cover and ruggedness similarly had a much stronger fit with genetic differentiation than an undifferentiated landscape. In most cases, landscape resistance surfaces transformed with 17.33-km-diameter moving windows were preferred, suggesting small-scale differences in habitat were unimportant at this large spatial extent. Despite the emergence of these overall patterns, there were differences in the selection of top models depending on the model selection criteria, suggesting research into the most appropriate criteria for landscape genetics is required. Overall, our results highlight the importance of differences in seasonal habitat preferences to patterns of gene flow and suggest the combination of habitat suitability modeling and linear mixed models with our resistance parameterization is a powerful approach to discerning the effects of landscape on gene flow.

Ecology and Evolution↗
USGS2015

Genomic single-nucleotide polymorphisms confirm that Gunnison and Greater sage-grouse are genetically well differentiated and that the Bi-State population is distinct

Sage-grouse are iconic, declining inhabitants of sagebrush habitats in western North America, and their management depends on an understanding of genetic variation across the landscape. Two distinct species of sage-grouse have been recognized, Greater ( Centrocercus urophasianus ) and Gunnison sage-grouse ( C. minimus ), based on morphology, behavior, and variation at neutral genetic markers. A parapatric group of Greater Sage-Grouse along the border of California and Nevada ("Bi-State") is also genetically distinct at the same neutral genetic markers, yet not different in behavior or morphology. Because delineating taxonomic boundaries and defining conservation units is often difficult in recently diverged taxa and can be further complicated by highly skewed mating systems, we took advantage of new genomic methods that improve our ability to characterize genetic variation at a much finer resolution. We identified thousands of single-nucleotide polymorphisms (SNPs) among Gunnison, Greater, and Bi-State sage-grouse and used them to comprehensively examine levels of genetic diversity and differentiation among these groups. The pairwise multilocus fixation index ( F ST ) was high (0.49) between Gunnison and Greater sage-grouse, and both principal coordinates analysis and model-based clustering grouped samples unequivocally by species. Standing genetic variation was lower within the Gunnison Sage-Grouse. The Bi-State population was also significantly differentiated from Greater Sage-Grouse, albeit more weakly ( F ST = 0.09), and genetic clustering results were consistent with reduced gene flow with Greater Sage-Grouse. No comparable genetic divisions were found within the Greater Sage-Grouse sample, which spanned the southern half of the range. Thus, we provide much stronger genetic evidence supporting the recognition of Gunnison Sage-Grouse as a distinct species with low genetic diversity. Further, our work confirms that the Bi-State population is differentiated from other Greater Sage-Grouse. The level of differentiation is much less than the divergence between Greater and Gunnison sage-grouse, supporting the idea that the Bi-State represents a unique population within the Greater Sage-Grouse. New genomic methods like the restriction-site-associated DNA (RAD-tag) method used here illustrate how increasing the number of markers and coverage of the genome can better characterize patterns of genetic variation, particularly among recently diverged taxa, providing vital information for conservation and management.

The Condor↗
USGS2015

Potential demographic and genetic effects of a sterilant applied to wild horse mares

Wild horse populations on western ranges can increase rapidly, resulting in the need for the Bureau of Land Management (BLM) to remove animals in order to protect the habitat that horses share with numerous other species. As an alternative to removals, BLM has sought to develop a long-term, perhaps even permanent, contraceptive to aid in reducing population growth rates. With long-term (perhaps even permanent) efficacy of contraception, however, comes increased concern about the genetic health of populations and about the potential for local extirpation. We used simulation modeling to examine the potential demographic and genetic consequences of applying a mare sterilant to wild horse populations. Using the VORTEX software package, we modeled the potential effects of a sterilant on 70 simulated populations having different initial sizes (7 values), growth rates (5 values), and genetic diversity (2 values). For each population, we varied the treatment rate of mares from 0 to 100 percent in increments of 10 percent. For each combination of these treatment levels, we ran 100 stochastic simulations, and we present the results in the form of tables and graphs showing mean population size after 20 years, mean number of removals after 20 years, mean probability of extirpation after 50 years, and mean heterozygosity after 50 years. By choosing one or two combinations of initial population size, population growth rate, and genetic diversity that best represent a herd of interest, a manager can assess the likely effects of a contraceptive program by examining the output tables and graphs representing the selected conditions.

Open-File Report↗
USGS2014

Development of 13 microsatellites for Gunnison Sage-grouse ( Centrocercus minimus ) using next-generation shotgun sequencing and their utility in Greater Sage-grouse ( Centrocercus urophasianus )

Gunnison Sage-grouse are an obligate sagebrush species that has experienced significant population declines and has been proposed for listing under the U.S. Endangered Species Act. In order to examine levels of connectivity among Gunnison Sage-grouse leks, we identified 13 novel microsatellite loci though next-generation shotgun sequencing, and tested them on the closely related Greater Sage-grouse. The number of alleles per locus ranged from 2 to 12. No loci were found to be linked, although 2 loci revealed significant departures from Hardy–Weinberg equilibrium or evidence of null alleles. While these microsatellites were designed for Gunnison Sage-grouse, they also work well for Greater Sage-grouse and could be used for numerous genetic questions including landscape and population genetics.

Conservation Genetics Resources↗
USGS2014

Two low coverage bird genomes and a comparison of reference-guided versus de novo genome assemblies

As a greater number and diversity of high-quality vertebrate reference genomes become available, it is increasingly feasible to use these references to guide new draft assemblies for related species. Reference-guided assembly approaches may substantially increase the contiguity and completeness of a new genome using only low levels of genome coverage that might otherwise be insufficient for de novo genome assembly. We used low-coverage (~3.5–5.5x) Illumina paired-end sequencing to assemble draft genomes of two bird species (the Gunnison Sage-Grouse, Centrocercus minimus , and the Clark's Nutcracker, Nucifraga columbiana ). We used these data to estimate de novo genome assemblies and reference-guided assemblies, and compared the information content and completeness of these assemblies by comparing CEGMA gene set representation, repeat element content, simple sequence repeat content, and GC isochore structure among assemblies. Our results demonstrate that even lower-coverage genome sequencing projects are capable of producing informative and useful genomic resources, particularly through the use of reference-guided assemblies.

PLoS ONE↗
USGS2014

Hierarchical spatial genetic structure in a distinct population segment of greater sage-grouse

Greater sage-grouse ( Centrocercus urophasianus ) within the Bi-State Management Zone (area along the border between Nevada and California) are geographically isolated on the southwestern edge of the species’ range. Previous research demonstrated that this population is genetically unique, with a high proportion of unique mitochondrial DNA (mtDNA) haplotypes and with significant differences in microsatellite allele frequencies compared to populations across the species’ range. As a result, this population was considered a distinct population segment (DPS) and was recently proposed for listing as threatened under the U.S. Endangered Species Act. A more comprehensive understanding of the boundaries of this genetically unique population (where the Bi-State population begins) and an examination of genetic structure within the Bi-State is needed to help guide effective management decisions. We collected DNA from eight sampling locales within the Bi-State (N = 181) and compared those samples to previously collected DNA from the two most proximal populations outside of the Bi-State DPS, generating mtDNA sequence data and amplifying 15 nuclear microsatellites. Both mtDNA and microsatellite analyses support the idea that the Bi-State DPS represents a genetically unique population, which has likely been separated for thousands of years. Seven mtDNA haplotypes were found exclusively in the Bi-State population and represented 73 % of individuals, while three haplotypes were shared with neighboring populations. In the microsatellite analyses both STRUCTURE and FCA separate the Bi-State from the neighboring populations. We also found genetic structure within the Bi-State as both types of data revealed differences between the northern and southern part of the Bi-State and there was evidence of isolation-by-distance. STRUCTURE revealed three subpopulations within the Bi-State consisting of the northern Pine Nut Mountains (PNa), mid Bi-State, and White Mountains (WM) following a north–south gradient. This genetic subdivision within the Bi-State is likely the result of habitat loss and fragmentation that has been exacerbated by recent human activities and the encroachment of singleleaf pinyon ( Pinus monophylla ) and juniper ( Juniperus spp.) trees. While genetic concerns may be only one of many priorities for the conservation and management of the Bi-State greater sage-grouse, we believe that they warrant attention along with other issues (e.g., quality of sagebrush habitat, preventing future loss of habitat). Management actions that promote genetic connectivity, especially with respect to WM and PNa, may be critical to the long-term viability of the Bi-State DPS.

California, Nevada↗
USGS2014

The historical distribution of Gunnison Sage-Grouse in Colorado

The historical distribution of Gunnison Sage-Grouse ( Centrocercus minimus ) in Colorado is described based on published literature, observations, museum specimens, and the known distribution of sagebrush ( Artemisia spp.). Historically, Gunnison Sage-Grouse were widely but patchily distributed in up to 22 counties in south-central and southwestern Colorado. The historical distribution of this species was south of the Colorado-Eagle river drainages primarily west of the Continental Divide. Potential contact areas with Greater Sage-Grouse ( C. urophasianus ) were along the Colorado-Eagle river system in Mesa, Garfield, and Eagle counties, west of the Continental Divide. Gunnison Sage-Grouse historically occupied habitats that were naturally highly fragmented by forested mountains and plateaus/mesas, intermountain basins without robust species of sagebrush, and river systems. This species adapted to use areas with more deciduous shrubs (i.e., Quercus spp., Amelanchier spp., Prunus spp.) in conjunction with sagebrush. Most areas historically occupied were small, linear, and patchily distributed within the overall landscape matrix. The exception was the large intermountain basin in Gunnison, Hinsdale, and Saguache counties. The documented distribution east of the Continental Divide within the large expanse of the San Luis Valley (Alamosa, Conejos, Costilla, and Rio Grande counties) was minimal and mostly on the eastern, northern, and southern fringes. Many formerly occupied habitat patches were vacant by the mid 1940s with extirpations continuing to the late 1990s. Counties from which populations were recently extirpated include Archuleta and Pitkin (1960s), and Eagle, Garfield, Montezuma, and Ouray (1990s).

Colorado↗
USGS2014

Genetic tools for wildlife management

Granted interim status in November, 2013, The Wildlife Society’s (TWS) Molecular Ecology Working Group aims to promote scientific advancement by applying molecular techniques to wildlife ecology, management, and conservation. The working group—composed of sci - entists from diverse backgrounds—met for the first time in Pittsburgh at the TWS Annual Conference held in October. Our overarching goal is to enhance awareness of molecular ecology and genetic applica - tions to wildlife biology and act as an informational and networking resource. During the group’s interim status, which runs for three years, we intend to focus on a broad scope of molecular ecology that is applicable to wildlife including genetic and ge - nomic methods, conservation genetics, non-invasive genetic population monitoring, landscape genetics, evolutionary genetics, and molecular forensics

The Wildlife Professional↗
USGS2013

Genetic characterization of the Pacific sheath-tailed bat (Emballonura semicaudata rotensis) using mitochondrial DNA sequence data

Emballonura semicaudata occurs in the southwestern Pacific and populations on many islands have declined or disappeared. One subspecies (E. semicaudata rotensis) occurs in the Northern Mariana Islands, where it has been extirpated from all but 1 island (Aguiguan). We assessed genetic similarity between the last population of E. s. rotensis and 2 other subspecies, and examined genetic diversity on Aguiguan. We sampled 12 E. s. rotensis, sequenced them at 3 mitochondrial loci, and compared them with published sequences from 2 other subspecies. All 12 E. s. rotensis had identical sequences in each of the 3 regions. Using cytochrome-b (Cytb) data E. s. rotensis was sister to E. s. palauensis in a clade separate from E. s. semicaudata. 12S ribosomal RNA (12S) sequences grouped all E. s. semicaudata in 1 clade with E. s. rotensis in a clade by itself. Genetic distances among the 3 subspecies at Cytb were smallest between E. s. palauensis and E. s. rotensis. Distance between E. s. semicaudata and the other 2 subspecies was not different from the distance between E. s. semicaudata and the full species E. raffrayana. A similar relationship was found using the 12S data. These distances are larger than those typically reported for mammalian subspecies using Cytb sequence and within the range of sister species.

Journal of Mammalogy↗
USGS2012

Sample design effects in landscape genetics

An important research gap in landscape genetics is the impact of different field sampling designs on the ability to detect the effects of landscape pattern on gene flow. We evaluated how five different sampling regimes (random, linear, systematic, cluster, and single study site) affected the probability of correctly identifying the generating landscape process of population structure. Sampling regimes were chosen to represent a suite of designs common in field studies. We used genetic data generated from a spatially-explicit, individual-based program and simulated gene flow in a continuous population across a landscape with gradual spatial changes in resistance to movement. Additionally, we evaluated the sampling regimes using realistic and obtainable number of loci (10 and 20), number of alleles per locus (5 and 10), number of individuals sampled (10-300), and generational time after the landscape was introduced (20 and 400). For a simulated continuously distributed species, we found that random, linear, and systematic sampling regimes performed well with high sample sizes (>200), levels of polymorphism (10 alleles per locus), and number of molecular markers (20). The cluster and single study site sampling regimes were not able to correctly identify the generating process under any conditions and thus, are not advisable strategies for scenarios similar to our simulations. Our research emphasizes the importance of sampling data at ecologically appropriate spatial and temporal scales and suggests careful consideration for sampling near landscape components that are likely to most influence the genetic structure of the species. In addition, simulating sampling designs a priori could help guide filed data collection efforts.

Conservation Genetics↗
USGS2012

Development and characterization of thirteen microsatellite loci in Clark's nutcracker (Nucifraga columbiana)

Clark’s nutcrackers are important seed dispersers for two widely-distributed western North American conifers, whitebark pine and limber pine, which are declining due to outbreaks of mountain pine beetle and white pine blister rust. Because nutcracker seed dispersal services are key to maintaining viable populations of these imperiled pines, knowledge of movement patterns of Clark’s nutcrackers helps managers understand local extinction risks for these trees. To investigate population structure within Clark’s nutcracker, we developed primers for and characterized 13 polymorphic microsatellite loci. In a screen of 22 individuals from one population, levels of variability ranged from 6 to 15 alleles. No loci were found to be linked, although 4 loci revealed significant departures from Hardy–Weinberg equilibrium and evidence of null alleles. These microsatellite loci will enable population genetic analyses of Clark’s nutcrackers, which could provide insights into the spatial relationships between nutcrackers and the trees they help disperse.

Conservation Genetics Resources↗
USGS2012

Rapid microsatellite identification from Illumina paired-end genomic sequencing in two birds and a snake

Identification of microsatellites, or simple sequence repeats (SSRs), can be a time-consuming and costly investment requiring enrichment, cloning, and sequencing of candidate loci. Recently, however, high throughput sequencing (with or without prior enrichment for specific SSR loci) has been utilized to identify SSR loci. The direct "Seq-to-SSR" approach has an advantage over enrichment-based strategies in that it does not require a priori selection of particular motifs, or prior knowledge of genomic SSR content. It has been more expensive per SSR locus recovered, however, particularly for genomes with few SSR loci, such as bird genomes. The longer but relatively more expensive 454 reads have been preferred over less expensive Illumina reads. Here, we use Illumina paired-end sequence data to identify potentially amplifiable SSR loci (PALs) from a snake (the Burmese python, Python molurus bivittatus ), and directly compare these results to those from 454 data. We also compare the python results to results from Illumina sequencing of two bird genomes (Gunnison Sage-grouse, Centrocercus minimus , and Clark's Nutcracker, Nucifraga columbiana ), which have considerably fewer SSRs than the python. We show that direct Illumina Seq-to-SSR can identify and characterize thousands of potentially amplifiable SSR loci for as little as $10 per sample – a fraction of the cost of 454 sequencing. Given that Illumina Seq-to-SSR is effective, inexpensive, and reliable even for species such as birds that have few SSR loci, it seems that there are now few situations for which prior hybridization is justifiable.

PLoS ONE↗
USGS2011

Effects of sample size, number of markers, and allelic richness on the detection of spatial genetic pattern

The influence of study design on the ability to detect the effects of landscape pattern on gene flow is one of the most pressing methodological gaps in landscape genetic research. To investigate the effect of study design on landscape genetics inference, we used a spatially-explicit, individual-based program to simulate gene flow in a spatially continuous population inhabiting a landscape with gradual spatial changes in resistance to movement. We simulated a wide range of combinations of number of loci, number of alleles per locus and number of individuals sampled from the population. We assessed how these three aspects of study design influenced the statistical power to successfully identify the generating process among competing hypotheses of isolation-by-distance, isolation-by-barrier, and isolation-by-landscape resistance using a causal modelling approach with partial Mantel tests. We modelled the statistical power to identify the generating process as a response surface for equilibrium and non-equilibrium conditions after introduction of isolation-by-landscape resistance. All three variables (loci, alleles and sampled individuals) affect the power of causal modelling, but to different degrees. Stronger partial Mantel r correlations between landscape distances and genetic distances were found when more loci were used and when loci were more variable, which makes comparisons of effect size between studies difficult. Number of individuals did not affect the accuracy through mean equilibrium partial Mantel r , but larger samples decreased the uncertainty (increasing the precision) of equilibrium partial Mantel r estimates. We conclude that amplifying more (and more variable) loci is likely to increase the power of landscape genetic inferences more than increasing number of individuals.

Molecular Ecology Resources↗