Search USGS⌕ Search

SEARCH · Search USGS

Results for “Molecular Ecology Resources”

Search indexed USGS publications on groundwater, aquifers, geologic maps, mineral resources and earthquakes. Explore source records by subject and place.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 19 recordsLinked to original sources
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↗
USGS2021

Time to get real with qPCR controls: The frequency of sample contamination and the informative power of negative controls in environmental DNA studies

Environmental (e)DNA methods have enabled rapid, sensitive and specific inferences of taxa presence throughout diverse fields of ecological study. However, use of eDNA results for decision-making has been impeded by uncertainties associated with false positive tests putatively caused by sporadic or systemic contamination. Sporadic contamination is a process that is inconsistent across samples and systemic contamination occurs consistently over a group of samples. Here, we used empirical data and laboratory experiments to (i) estimate the sporadic contamination rate for each stage of a common, targeted eDNA workflow employing best practice quality control measures under simulated conditions of rare and common target DNA presence, (ii) determine the rate at which negative controls (i.e., “blanks”) detect varying concentrations of systemic contamination, and (iii) estimate the effort that would be required to consistently detect sporadic and systemic contamination. Sporadic contamination rates were very low across all eDNA workflow steps, and, therefore, an intractably high number of negative controls (>100) would be required to determine occurrence of sporadic contamination with any certainty. Contrarily, detection of intentionally introduced systemic contamination was more consistent; therefore, very few negative controls (<5) would be needed to consistently alert to systemic contamination. These results have considerable implications to eDNA study design when resources for sample analyses are constrained.

Molecular Ecology Resources↗
USGS2012

RAD sequencing yields a high success rate for westslope cutthroat and rainbow trout species-diagnostic SNP assays

Hybridization with introduced rainbow trout threatens most native westslope cutthroat trout populations. Understanding the genetic effects of hybridization and introgression requires a large set of high-throughput, diagnostic genetic markers to inform conservation and management. Recently, we identified several thousand candidate single-nucleotide polymorphism (SNP) markers based on RAD sequencing of 11 westslope cutthroat trout and 13 rainbow trout individuals. Here, we used flanking sequence for 56 of these candidate SNP markers to design high-throughput genotyping assays. We validated the assays on a total of 92 individuals from 22 populations and seven hatchery strains. Forty-six assays (82%) amplified consistently and allowed easy identification of westslope cutthroat and rainbow trout alleles as well as heterozygote controls. The 46 SNPs will provide high power for early detection of population admixture and improved identification of hybrid and nonhybridized individuals. This technique shows promise as a very low-cost, reliable and relatively rapid method for developing and testing SNP markers for nonmodel organisms with limited genomic resources.

Molecular Ecology Resources↗
USGS2020

A GT-seq panel for walleye (Sander vitreus) provides important insights for efficient development and implementation of amplicon panels in non-model organisms

Targeted amplicon sequencing methods, such as genotyping-in-thousands by sequencing (GT-seq), facilitate rapid, accurate, and cost-effective analysis of hundreds of genetic loci in thousands of individuals. Development of GT-seq panels is nontrivial, but studies describing trade-offs associated with different steps of GT-seq panel development are rare. Here, we construct a dual-purpose GT-seq panel for walleye ( Sander vitreus ), discuss trade-offs associated with different development and genotyping approaches, and provide suggestions for researchers constructing their own GT-seq panels. Our GT-seq panel was developed using an ascertainment set consisting of restriction site-associated DNA data from 954 individuals sampled from 23 populations in Minnesota and Wisconsin, USA. We conducted simulations to test the utility of all loci for parentage analysis and genetic stock identification and designed 600 primer pairs to maximize joint accuracy for these analyses. We then performed three rounds of primer optimization to remove loci that overamplified and our final panel consisted of 436 loci. We also explored different approaches for DNA extraction, multiplexed polymerase chain reaction (PCR) amplification, and cleanup steps during the GT-seq process and discovered the following: (i) inexpensive Chelex extractions performed well for genotyping; (ii) the exonuclease I and shrimp alkaline phosphatase (ExoSAP) procedure included in some current protocols did not improve results substantially and was probably unnecessary; and (iii) it was possible to PCR amplify panels separately and combine them prior to adapter ligation. Well-optimized GT-seq panels are valuable resources for conservation genetics and our findings and suggestions should aid in their construction in myriad taxa.

Minnesota, Wisconsin↗
USGS2023

A new genomic resource to enable standardized surveys of SNPs across the native range of brook trout (Salvelinus fontinalis)

Understanding how genetic diversity is distributed across spatiotemporal scales in species of conservation or management concern is critical for identifying large-scale mechanisms affecting local conservation status and implementing large-scale biodiversity monitoring programmes. However, cross-scale surveys of genetic diversity are often impractical within single studies, and combining datasets to increase spatiotemporal coverage is frequently impeded by using different sets of molecular markers. Recently developed molecular tools make surveys based on standardized single-nucleotide polymorphism (SNP) panels more feasible than ever, but require existing genomic information. Here, we conduct the first survey of genome-wide SNPs across the native range of brook trout (Salvelinus fontinalis), a cold-adapted species that has been the focus of considerable conservation and management effort across eastern North America. Our dataset can be leveraged to easily design SNP panels that allow datasets to be combined for large-scale analyses. We performed restriction site-associated DNA sequencing for wild brook trout from 82 locations spanning much of the native range and domestic brook trout from 24 hatchery strains used in stocking efforts. We identified over 24,000 SNPs distributed throughout the brook trout genome. We explored the ability of these SNPs to resolve relationships across spatial scales, including population structure and hatchery admixture. Our dataset captures a wide spectrum of genetic diversity in native brook trout, offering a valuable resource for developing SNP panels. We highlight potential applications of this resource with the goal of increasing the integration of genomic information into decision-making for brook trout and other species of conservation or management concern.

Molecular Ecology Resources↗
USGS2020

Characterization of a Y-specific duplication/insertion of the anti-Mullerian hormone type II receptor gene based on a chromosome-scale genome assembly of yellow perch, Perca flavescens

Yellow perch, Perca flavescens , is an ecologically and economically important species native to a large portion of the northern United States and southern Canada and is also a promising candidate species for aquaculture. However, no yellow perch reference genome has been available to facilitate improvements in both fisheries and aquaculture management practices. By combining Oxford Nanopore Technologies long-reads, 10X Genomics Illumina short linked reads and a chromosome contact map produced with Hi-C, we generated a high-continuity chromosome-scale yellow perch genome assembly of 877.4 Mb. It contains, in agreement with the known diploid chromosome yellow perch count, 24 chromosome-size scaffolds covering 98.8% of the complete assembly (N50 = 37.4 Mb, L50 = 11). We also provide a first characterization of the yellow perch sex determination locus that contains a male-specific duplicate of the anti-Mullerian hormone type II receptor gene ( amhr2by ) inserted at the proximal end of the Y chromosome (chromosome 9). Using this sex-specific information, we developed a simple PCR genotyping assay which accurately differentiates XY genetic males ( amhr2by + ) from XX genetic females ( amhr2by − ). Our high-quality genome assembly is an important genomic resource for future studies on yellow perch ecology, toxicology, fisheries and aquaculture research. In addition, characterization of the amhr2by gene as a candidate sex-determining gene in yellow perch provides a new example of the recurrent implication of the transforming growth factor beta pathway in fish sex determination, and highlights gene duplication as an important genomic mechanism for the emergence of new master sex determination genes.

Molecular Ecology Resources↗
USGS2023

A novel assembly pipeline and functional annotations for targeted sequencing: A case study on the globally threatened Margaritiferidae (Bivalvia: Unionida)

The proliferation of genomic sequencing approaches has significantly impacted the field of phylogenetics. Target capture approaches provide a cost-effective, fast and easily applied strategy for phylogenetic inference of non-model organisms. However, several existing target capture processing pipelines are incapable of incorporating whole genome sequencing (WGS). Here, we develop a new pipeline for capture and de novo assembly of the targeted regions using whole genome re-sequencing reads. This new pipeline captured targeted loci accurately, and given its unbiased nature, can be used with any target capture probe set. Moreover, due to its low computational demand, this new pipeline may be ideal for users with limited resources and when high-coverage sequencing outputs are required. We demonstrate the utility of our approach by incorporating WGS data into the first comprehensive phylogenomic reconstruction of the freshwater mussel family Margaritiferidae. We also provide a catalogue of well-curated functional annotations of these previously uncharacterized freshwater mussel-specific target regions, representing a complementary tool for scrutinizing phylogenetic inferences while expanding future applications of the probe set.

Molecular Ecology Resources↗
USGS2013

Transcriptome resources for the frogs Lithobates clamitans and Pseudacris regilla , emphasizing antimicrobial peptides and conserved loci for phylogenetics

We developed genetic resources for two North American frogs, Lithobates clamitans and Pseudacris regilla , widespread native amphibians that are potential indicator species of environmental health. For both species, mRNA from multiple tissues was sequenced using 454 technology. De novo assemblies with Mira3 resulted in 50 238 contigs (N50 = 687 bp) and 48 213 contigs (N50 = 686 bp) for L. clamitans and P. regilla , respectively, after clustering with CD-Hit-EST and purging contigs below 200 bp. We performed BLASTX similarity searches against the Xenopus tropicalis proteome and, for predicted ORFs, HMMER similarity searches against the Pfam-A database. Because there is broad interest in amphibian immune factors, we manually annotated putative antimicrobial peptides. To identify conserved regions suitable for amplicon resequencing across a broad taxonomic range, we performed an additional assembly of public short-read transcriptome data derived from two species of the genus Rana and identified reciprocal best TBLASTX matches among all assemblies. Although P. regilla , a hylid frog, is substantially more diverged from the ranid species, we identified 56 genes that were sufficiently conserved to allow nondegenerate primer design with Primer3. In addition to providing a foundation for comparative genomics and quantitative gene expression analysis, our results enable quick development of nuclear sequence-based markers for phylogenetics or population genetics.

Molecular Ecology Resources↗
USGS2025

Development of high-throughput genomic resources to inform white-tailed deer population and disease management

White-tailed deer ( Odocoileus virginianus ) are the most abundant and widespread cervid in North America. Genetic data are used as a tool to monitor populations and make management decisions for this game species. However, the development and use of genomic tools that can generate a set of markers suitable for longitudinal genomic data collection, whether for management purposes or to study the demographic and evolutionary processes of widely distributed species, have been challenging. This is mainly due to the cost required to fully implement and interpret the data produced. Here, we generated whole genome resequencing data for 44 free-ranging deer from three regions in their central and eastern North American range and identified over 89 million single nucleotide polymorphisms (SNPs). We used a subset of these SNPs to develop two nested SNP tools, a high-density array (702,183 SNPs) and a medium-density array (72,723 SNPs) to support deer and chronic wasting disease (CWD) management and research. SNPs were selected to ensure an even distribution across scaffolds of the reference genome and include SNPs associated with CWD susceptibility. Using genotyping results for 469 deer from 15 states in the US and Mexico generated by the high-density array and 1335 deer from 18 states generated by the medium-density array, we assessed genotyping success across different populations and explored some insights into population structure. These genomic tools offer a standard set of markers that will enable researchers and managers to address important questions related to white-tailed deer and CWD management. Our SNP arrays also offer the opportunity to examine aspects of white-tailed deer ecology and evolutionary history that were previously difficult to address.

Molecular Ecology Resources↗
USGS2022

Haploid gynogens facilitate disomic marker development in paleotetraploid sturgeons

Acipenseriformes (sturgeons and paddlefishes) are of substantial conservation concern, and development of genomic resources for these species is difficult due to past whole genome duplication. Development of disomic markers for polyploid organisms can be challenging due to difficulty in resolving alleles at a single locus from those among duplicated loci. In this study, we detail the development of disomic markers for the endangered pallid sturgeon ( Scaphirhynchus albus ) found in North America. One of the strategies for pallid sturgeon conservation is to stock U.S. rivers with offspring of pure pallid sturgeon, but introgression with the sympatric shovelnose sturgeon ( S. platorynchus ) threatens pallid sturgeon genetic integrity. Currently, 19 microsatellite loci are used to differentiate between both species and their hybrids, but the markers are insufficient to robustly identify backcrosses. We performed double digest restriction site-associated DNA sequencing (ddRADseq) on shovelnose sturgeon haploid gynogens to produce a reduced-representation genomic reference. Contiguous sequences that were heterozygous within a haploid individual were flagged as potentially encompassing multiple loci. Approximately 60 individuals of each species from two management units were sequenced, and reads were mapped to the haploid reference to identify single nucleotide polymorphisms (SNPs) at individual loci. The final data set contained 11,082 microhaplotyped loci which offer at least an order of magnitude greater resolution for species discrimination than the current panel of 19 microsatellites. These markers will be used to examine a larger sample of Scaphirhynchus individuals throughout their ranges to determine the extent and trajectory of hybridization.

Molecular Ecology Resources↗
USGS2013

Factors influencing detection of eDNA from a stream-dwelling amphibian

Environmental DNA (eDNA) methods for detecting and estimating abundance of aquatic species are emerging rapidly, but little is known about how processes such as secretion rate, environmental degradation, and time since colonization or extirpation from a given site affect eDNA measurements. Using stream-dwelling salamanders and quantitative PCR (qPCR) analysis, we conducted three experiments to assess eDNA: (i) production rate; (ii) persistence time under different temperature and light conditions; and (iii) detectability and concentration through time following experimental introduction and removal of salamanders into previously unoccupied streams. We found that 44–50 g individuals held in aquaria produced 77 ng eDNA/h for 2 h, after which production either slowed considerably or began to equilibrate with degradation. eDNA in both full-sun and shaded treatments degraded exponentially to <1% of the original concentration after 3 days. eDNA was no longer detectable in full-sun samples after 8 days, whereas eDNA was detected in 20% of shaded samples after 11 days and 100% of refrigerated control samples after 18 days. When translocated into unoccupied streams, salamanders were detectable after 6 h, but only when densities were relatively high (0.2481 individuals/m 2 ) and when samples were collected within 5 m of the animals. Concentrations of eDNA detected were very low and increased steadily from 6–24 h after introduction, reaching 0.0022 ng/L. Within 1 h of removing salamanders from the stream, eDNA was no longer detectable. These results suggest that eDNA detectability and concentration depend on production rates of individuals, environmental conditions, density of animals, and their residence time.

Molecular Ecology Resources↗
USGS2014

A likelihood-based approach for assessment of extra-pair paternity and conspecific brood parasitism in natural populations

Genotypes are frequently used to assess alternative reproductive strategies such as extra-pair paternity and conspecific brood parasitism in wild populations. However, such analyses are vulnerable to genotyping error or molecular artifacts that can bias results. For example, when using multilocus microsatellite data, a mismatch at a single locus, suggesting the offspring was not directly related to its putative parents, can occur quite commonly even when the offspring is truly related. Some recent studies have advocated an ad-hoc rule that offspring must differ at more than one locus in order to conclude that they are not directly related. While this reduces the frequency with which true offspring are identified as not directly related young, it also introduces bias in the opposite direction, wherein not directly related young are categorized as true offspring. More importantly, it ignores the additional information on allele frequencies which would reduce overall bias. In this study, we present a novel technique for assessing extra-pair paternity and conspecific brood parasitism using a likelihood-based approach in a new version of program cervus. We test the suitability of the technique by applying it to a simulated data set and then present an example to demonstrate its influence on the estimation of alternative reproductive strategies.

Molecular Ecology Resources↗
USGS2016

Detecting the movement and spawning activity of bigheaded carps with environmental DNA

Bigheaded carps are invasive fishes threatening to invade the Great Lakes basin and establish spawning populations, and have been monitored using environmental DNA (eDNA). Not only does eDNA hold potential for detecting the presence of species, but may also allow for quantitative comparisons like relative abundance of species across time or space. We examined the relationships among bigheaded carp movement, hydrography, spawning and eDNA on the Wabash River, IN, USA. We found positive relationships between eDNA and movement and eDNA and hydrography. We did not find a relationship between eDNA and spawning activity in the form of drifting eggs. Our first finding demonstrates how eDNA may be used to monitor species abundance, whereas our second finding illustrates the need for additional research into eDNA methodologies. Current applications of eDNA are widespread, but the relatively new technology requires further refinement.

Molecular Ecology Resources↗
USGS2016

SNP discovery in candidate adaptive genes using exon capture in a free-ranging alpine ungulate

Identification of genes underlying genomic signatures of natural selection is key to understanding adaptation to local conditions. We used targeted resequencing to identify SNP markers in 5321 candidate adaptive genes associated with known immunological, metabolic and growth functions in ovids and other ungulates. We selectively targeted 8161 exons in protein-coding and nearby 5′ and 3′ untranslated regions of chosen candidate genes. Targeted sequences were taken from bighorn sheep ( Ovis canadensis ) exon capture data and directly from the domestic sheep genome ( Ovis aries v. 3; oviAri3). The bighorn sheep sequences used in the Dall's sheep ( Ovis dalli dalli ) exon capture aligned to 2350 genes on the oviAri3 genome with an average of 2 exons each. We developed a microfluidic qPCR-based SNP chip to genotype 476 Dall's sheep from locations across their range and test for patterns of selection. Using multiple corroborating approaches ( lositan and bayescan ), we detected 28 SNP loci potentially under selection. We additionally identified candidate loci significantly associated with latitude, longitude, precipitation and temperature, suggesting local environmental adaptation. The three methods demonstrated consistent support for natural selection on nine genes with immune and disease-regulating functions (e.g. Ovar-DRA, APC, BATF2, MAGEB18), cell regulation signalling pathways (e.g. KRIT1, PI3K, ORRC3), and respiratory health (CYSLTR1). Characterizing adaptive allele distributions from novel genetic techniques will facilitate investigation of the influence of environmental variation on local adaptation of a northern alpine ungulate throughout its range. This research demonstrated the utility of exon capture for gene-targeted SNP discovery and subsequent SNP chip genotyping using low-quality samples in a nonmodel species.

Molecular Ecology Resources↗
USGS2016

Detection limits of quantitative and digital PCR assays and their influence in presence-absence surveys of environmental DNA

A set of universal guidelines is needed to determine the limit of detection (LOD) in PCR-based analyses of low concentration DNA. In particular, environmental DNA (eDNA) studies require sensitive and reliable methods to detect rare and cryptic species through shed genetic material in environmental samples. Current strategies for assessing detection limits of eDNA are either too stringent or subjective, possibly resulting in biased estimates of species’ presence. Here, a conservative LOD analysis grounded in analytical chemistry is proposed to correct for overestimated DNA concentrations predominantly caused by the concentration plateau, a nonlinear relationship between expected and measured DNA concentrations. We have used statistical criteria to establish formal mathematical models for both quantitative and droplet digital PCR. To assess the method, a new Grass Carp ( Ctenopharyngodon idella ) TaqMan assay was developed and tested on both PCR platforms using eDNA in water samples. The LOD adjustment reduced Grass Carp occupancy and detection estimates while increasing uncertainty – indicating that caution needs to be applied to eDNA data without LOD correction. Compared to quantitative PCR, digital PCR had higher occurrence estimates due to increased sensitivity and dilution of inhibitors at low concentrations. Without accurate LOD correction, species occurrence and detection probabilities based on eDNA estimates are prone to a source of bias that cannot be reduced by an increase in sample size or PCR replicates. Other applications also could benefit from a standardized LOD such as GMO food analysis, and forensic and clinical diagnostics.

Molecular Ecology Resources↗
USGS2017

eDNAoccupancy: An R package for multi-scale occupancy modeling of environmental DNA data

In this article we describe eDNAoccupancy, an R package for fitting Bayesian, multi-scale occupancy models. These models are appropriate for occupancy surveys that include three, nested levels of sampling: primary sample units within a study area, secondary sample units collected from each primary unit, and replicates of each secondary sample unit. This design is commonly used in occupancy surveys of environmental DNA (eDNA). eDNAoccupancy allows users to specify and fit multi-scale occupancy models with or without covariates, to estimate posterior summaries of occurrence and detection probabilities, and to compare different models using Bayesian model-selection criteria. We illustrate these features by analyzing two published data sets: eDNA surveys of a fungal pathogen of amphibians and eDNA surveys of an endangered fish species.

Molecular Ecology Resources↗
USGS2019

The future is now: Amplicon sequencing and sequence capture usher in the conservation genomics era

The genomics revolution has initiated a new era of population genetics where genome-wide data are frequently used to understand complex patterns of population structure and selection. However, the application of genomic tools to inform management and conservation has been somewhat rare outside a few well studied species. Fortunately, two recently developed approaches, amplicon sequencing and sequence capture, have the potential to significantly advance the field of conservation genomics. Here, amplicon sequencing refers to highly multiplexed PCR followed by high-throughput sequencing (e.g., GTseq), and sequence capture refers to using capture probes to isolate loci from reduced-representation libraries (e.g., Rapture). Both approaches allow sequencing of thousands of individuals at relatively low costs, do not require any specialized equipment for library preparation, and generate data that can be analyzed without sophisticated computational infrastructure. Here, we discuss the advantages and disadvantages of each method and provide a decision framework for geneticists who are looking to integrate these methods into their research programme. While it will always be important to consider the specifics of the biological question and system, we believe that amplicon sequencing is best suited for projects aiming to genotype <500 loci on many individuals (>1,500) or for species where continued monitoring is anticipated (e.g., long-term pedigrees). Sequence capture, on the other hand, is best applied to projects including fewer individuals or where >500 loci are required. Both of these techniques should smooth the transition from traditional genetic techniques to genomics, helping to usher in the conservation genomics era.

Molecular Ecology Resources↗
USGS2020

Attack of the PCR clones: Rates of clonality have little effect on RAD-seq genotype calls

Interpretation of high-throughput sequence data requires an understanding of how decisions made during bioinformatic data processing can influence results. One source of bias that is often cited is PCR clones (or PCR duplicates). PCR clones are common in restriction site-associated sequencing (RAD-seq) data sets, which are increasingly being used for molecular ecology. To determine the influence PCR clones and the bioinformatic handling of clones have on genotyping, we evaluate four RAD-seq data sets. Data sets were compared before and after clones were removed to estimate the number of clones present in RAD-seq data, quantify how often the presence of clones in a data set causes genotype calls to change compared to when clones were removed, investigate the mechanisms that lead to genotype call changes and test whether clones bias heterozygosity estimates. Our RAD-seq data sets contained 30%–60% PCR clones, but 95% of RAD-tags had five or fewer clones. Relatively few genotypes changed once clones were removed (5%–10%), and the vast majority of these changes (98%) were associated with genotypes switching from a called to no-call state or vice versa. PCR clones had a larger influence on genotype calls in individuals with low read depth but appeared to influence genotype calls at all loci similarly. Removal of PCR clones reduced the number of called genotypes by 2% but had almost no influence on estimates of heterozygosity. As such, while steps should be taken to limit PCR clones during library preparation, PCR clones are likely not a substantial source of bias for most RAD-seq studies.

Molecular Ecology Resources↗