Development of eighteen novel microsatellite loci for Masked Bobwhite, Colinus virginianus ridgwayi.

BACKGROUND: Masked Bobwhite (Colinus virginianus ridgwayi) is a critically-endangered New World quail species endemic to Sonoran Desert grasslands of North America. It suffered severe population declines during the nineteenth and twentieth centuries, with its persistence now reliant upon a captive breeding program that requires careful genetic management to maintain extant genetic diversity. Although nuclear microsatellite DNA markers existed for the closely related Northern Bobwhite (C. virginianus), none were available for Masked Bobwhite to inform necessary management decisions. METHODS AND RESULTS: Paired-end Illumina© sequencing was conducted to screen the Masked Bobwhite genome for microsatellite loci. We identified 18 loci exhibiting high polymorphism and limited deviations from genetic equilibrium expectations. These loci were amplified in 78 individuals. Familial relationships were reconstructed via sibship methods and compared to manually-curated pedigree data. Thirteen of fifteen full-sibling groups in the pedigree were exactly reconstructed (86.6%). Three other full-sibling groups partially matched pedigree relationships with high statistical confidence, and likely represented pedigree inaccuracies. Four additional full-sibling pairs were identified with low statistical confidence and likely resulted from analytical artifacts. CONCLUSIONS: The novel microsatellite loci accurately reconstructed parent-offspring and sibling relationships. These loci will be useful for guiding genetic management decisions and identifying pedigree inaccuracies in the captive breeding program.

Genomics and conservation: Guidance from training to analyses and applications.

Environmental change is intensifying the biodiversity crisis and threatening species across the tree of life. Conservation genomics can help inform conservation actions and slow biodiversity loss. However, more training, appropriate use of novel genomic methods and communication with managers are needed. Here, we review practical guidance to improve applied conservation genomics. We share insights aimed at ensuring effectiveness of conservation actions around three themes: (1) improving pedagogy and training in conservation genomics including for online global audiences, (2) conducting rigorous population genomic analyses properly considering theory, marker types and data interpretation and (3) facilitating communication and collaboration between managers and researchers. We aim to update students and professionals and expand their conservation toolkit with genomic principles and recent approaches for conserving and managing biodiversity. The biodiversity crisis is a global problem and, as such, requires international involvement, training, collaboration and frequent reviews of the literature and workshops as we do here.

AdmixPipe v3: facilitating population structure delimitation from SNP data.

Summary: Quantifying genetic clusters (=populations) from genotypic data is a fundamental, but non-trivial task for population geneticists that is compounded by: hierarchical population structure, diverse analytical methods, and complex software dependencies. AdmixPipe v3 ameliorates many of these issues in a single bioinformatic pipeline that facilitates all facets of population structure analysis by integrating outputs generated by several popular packages (i.e. CLUMPAK, EvalAdmix). The pipeline interfaces disparate software packages to parse Admixture outputs and conduct EvalAdmix analyses in the context of multimodal population structure results identified by CLUMPAK. We further streamline these tasks by packaging AdmixPipe v3 within a Docker container to create a standardized analytical environment that allows for complex analyses to be replicated by different researchers. This also grants operating system flexibility and mitigates complex software dependencies. Availability and implementation: Source code, documentation, example files, and usage examples are freely available at https://github.com/stevemussmann/admixturePipeline. Installation is facilitated via Docker container available from https://hub.docker.com/r/mussmann/admixpipe. Usage under Windows operating systems requires the Windows Subsystem for Linux.

Development of twenty-one novel microsatellite loci for Gila topminnow, Poeciliopsis occidentalis occidentalis.

BACKGROUND: Gila topminnow (Poeciliopsis occidentalis occidentalis) was once highly abundant throughout the Lower Colorado River Basin of the southwestern United States. However, this Sonoran Desert endemic suffered extreme population declines over the past century because of habitat degradation and nonnative species introductions. Much of the prior conservation genetic work conducted on the species relied upon a small number of microsatellite loci; many exhibiting low variability in extant populations. Consequently, there was a need for additional microsatellite loci to provide high-resolution delimitation of populations for conservation purposes. METHODS AND RESULTS: Paired-end Illumina sequencing was utilized to screen the Gila topminnow genome for novel microsatellite loci. We identified 21 novel loci that exhibited no deviations from expectations of genetic equilibrium, and cross-amplified in Yaqui topminnow (P. o. sonoriensis). These loci were amplified from 401 samples representing eight populations of Gila topminnow and Yaqui topminnow. Although diversity was low for all populations (observed heterozygosity = 0.12 to 0.45), these novel markers provided ample power to identify population of origin for each individual in Bayesian assignment tests. CONCLUSIONS: This novel set of microsatellite loci provide a useful genetic tool to assess population genetic parameters of the endangered Gila topminnow and delineate populations for identifying conservation priorities. The cross-amplification of these loci in Yaqui topminnow shows promise for application to other Poeciliopsis species of Mexico and Central America.

Population connectivity in voles (Microtus sp.) as a gauge for tall grass prairie restoration in midwestern North America.

Ecological restoration can promote biodiversity conservation in anthropogenically fragmented habitats, but effectiveness of these management efforts need to be statistically validated to determine 'success.' One such approach is to gauge the extent of recolonization as a measure of landscape permeability and, in turn, population connectivity. In this context, we estimated dispersal and population connectivity in prairie vole (Microtus ochrogaster; N = 231) and meadow vole (M. pennsylvanicus; N = 83) within five tall-grass prairie restoration sites embedded within the agricultural matrix of midwestern North America. We predicted that vole dispersal would be constrained by the extent of agricultural land surrounding restored habitat patches, spatially isolating vole populations and resulting in significant genetic structure. We first employed genetic assignment tests based on 15 microsatellite DNA loci to validate field-derived species-designations, then tested reclassified samples with multivariate and Bayesian clustering to assay for spatial and temporal genetic structure. Population connectivity was further evaluated by calculating pairwise FST, then potential demographic effects explored by computing migration rates, effective population size (Ne), and average relatedness (r). Genetic species assignments reclassified 25% of initial field identifications (N = 11 M. ochrogaster; N = 67 M. pennsylvanicus). In M. ochrogaster population connectivity was high across the study area, reflected in little to no spatial or temporal genetic structure. In M. pennsylvanicus genetic structure was detected, but relatedness estimates identified it as kin-clustering instead, underscoring social behavior among populations rather than spatial isolation as the cause. Estimates of Ne and r were stable across years, reflecting high dispersal and demographic resilience. Combined, these metrics suggest the agricultural matrix is highly permeable for voles and does not impede dispersal. High connectivity observed confirms that the restored landscape is productive and permeable for specific management targets such as voles and also demonstrates population genetic assays as a tool to statistically evaluate effectiveness of conservation initiatives.

Taxonomic Uncertainty and the Anomaly Zone: Phylogenomics Disentangle a Rapid Radiation to Resolve Contentious Species (Gila robusta Complex) in the Colorado River.

Species are indisputable units for biodiversity conservation, yet their delimitation is fraught with both conceptual and methodological difficulties. A classic example is the taxonomic controversy surrounding the Gila robusta complex in the lower Colorado River of southwestern North America. Nominal species designations were originally defined according to weakly diagnostic morphological differences, but these conflicted with subsequent genetic analyses. Given this ambiguity, the complex was re-defined as a single polytypic unit, with the proposed "threatened" status under the U.S. Endangered Species Act of two elements being withdrawn. Here we re-evaluated the status of the complex by utilizing dense spatial and genomic sampling (n = 387 and >22 k loci), coupled with SNP-based coalescent and polymorphism-aware phylogenetic models. In doing so, we found that all three species were indeed supported as evolutionarily independent lineages, despite widespread phylogenetic discordance. To juxtapose this discrepancy with previous studies, we first categorized those evolutionary mechanisms driving discordance, then tested (and subsequently rejected) prior hypotheses which argued phylogenetic discord in the complex was driven by the hybrid origin of Gila nigra. The inconsistent patterns of diversity we found within G. robusta were instead associated with rapid Plio-Pleistocene drainage evolution, with subsequent divergence within the "anomaly zone" of tree space producing ambiguities that served to confound prior studies. Our results not only support the resurrection of the three species as distinct entities but also offer an empirical example of how phylogenetic discordance can be categorized within other recalcitrant taxa, particularly when variation is primarily partitioned at the species level.

Defining relictual biodiversity: Conservation units in speckled dace (Leuciscidae: Rhinichthys osculus) of the Greater Death Valley ecosystem.

The tips in the tree of life serve as foci for conservation and management, yet clear delimitations are masked by inherent variance at the species-population interface. Analyses using thousands of nuclear loci can potentially sort inconsistencies, yet standard categories applied to this parsing are themselves potentially conflicting and/or subjective [e.g., DPS (distinct population segments); DUs (Diagnosable Units-Canada); MUs (management units); SSP (subspecies); ESUs (Evolutionarily Significant Units); and UIEUs (uniquely identified evolutionary units)]. One potential solution for consistent categorization is to create a comparative framework by accumulating statistical results from independent studies and evaluating congruence among data sets. Our study illustrates this approach in speckled dace (Leuciscidae: Rhinichthys osculus) endemic to two basins (Owens and Amargosa) in the Death Valley ecosystem. These fish persist in the Mojave Desert as isolated Plio-Pleistocene relicts and are of conservation concern, but lack formal taxonomic descriptions/designations. Double digest RAD (ddRAD) methods identified 14,355 SNP loci across 10 populations (N = 140). Species delimitation analyses [multispecies coalescent (MSC) and unsupervised machine learning (UML)] delineated four putative ESUs. F ST outlier loci (N = 106) were juxtaposed to uncover the potential for localized adaptations. We detected one hybrid population that resulted from upstream reconnection of habitat following contemporary pluvial periods, whereas remaining populations represent relics of ancient tectonism within geographically isolated springs and groundwater-fed streams. Our study offers three salient conclusions: a blueprint for a multifaceted delimitation of conservation units; a proposed mechanism by which criteria for intraspecific biodiversity can be potentially standardized; and a strong argument for the proactive management of critically endangered Death Valley ecosystem fishes.

ADMIXPIPE: population analyses in ADMIXTURE for non-model organisms.

BACKGROUND: Research on the molecular ecology of non-model organisms, while previously constrained, has now been greatly facilitated by the advent of reduced-representation sequencing protocols. However, tools that allow these large datasets to be efficiently parsed are often lacking, or if indeed available, then limited by the necessity of a comparable reference genome as an adjunct. This, of course, can be difficult when working with non-model organisms. Fortunately, pipelines are currently available that avoid this prerequisite, thus allowing data to be a priori parsed. An oft-used molecular ecology program (i.e., STRUCTURE), for example, is facilitated by such pipelines, yet they are surprisingly absent for a second program that is similarly popular and computationally more efficient (i.e., ADMIXTURE). The two programs differ in that ADMIXTURE employs a maximum-likelihood framework whereas STRUCTURE uses a Bayesian approach, yet both produce similar results. Given these issues, there is an overriding (and recognized) need among researchers in molecular ecology for bioinformatic software that will not only condense output from replicated ADMIXTURE runs, but also infer from these data the optimal number of population clusters (K). RESULTS: Here we provide such a program (i.e., ADMIXPIPE) that (a) filters SNPs to allow the delineation of population structure in ADMIXTURE, then (b) parses the output for summarization and graphical representation via CLUMPAK. Our benchmarks effectively demonstrate how efficient the pipeline is for processing large, non-model datasets generated via double digest restriction-site associated DNA sequencing (ddRAD). Outputs not only parallel those from STRUCTURE, but also visualize the variation among individual ADMIXTURE runs, so as to facilitate selection of the most appropriate K-value. CONCLUSIONS: ADMIXPIPE successfully integrates ADMIXTURE analysis with popular variant call format (VCF) filtering software to yield file types readily analyzed by CLUMPAK. Large population genomic datasets derived from non-model organisms are efficiently analyzed via the parallel-processing capabilities of ADMIXTURE. ADMIXPIPE is distributed under the GNU Public License and freely available for Mac OSX and Linux platforms at: https://github.com/stevemussmann/admixturePipeline .

Multi-targeted management of upland game birds at the agroecosystem interface in midwestern North America.

Despite its imperative, biodiversity conservation is chronically underfunded, a deficiency that often forces management agencies to prioritize. Single-species recovery thus becomes a focus (often with socio-political implications), whereas a more economical approach would be the transition to multi-targeted management (= MTM). This challenge is best represented in Midwestern North America where biodiversity has been impacted by 300+ years of chronic anthropogenic disturbance such that native tall-grass prairie is now supplanted by an agroecosystem. Here, we develop an MTM with a population genetic metric to collaboratively manage three Illinois upland gamebirds: common pheasant (Phasianus colchicus; pheasant), northern bobwhite quail (Colinus virginianus; quail), and threatened-endangered (T&E) greater prairie chicken (Tympanuchus cupido pinnatus; prairie chicken). We first genotyped our study pheasant at 19 microsatellite DNA loci and identified three captive breeding stocks (N = 143; IL Department of Natural Resources) as being significantly bottlenecked, with relatedness >1st-cousin (µR = 0.158). 'Wild' (non-stocked) pheasant [N = 543; 14 Pheasant-Habitat-Areas (PHAs)] were also bottlenecked, significantly interrelated (µR = 0.150) and differentiated (µFST = 0.047), yet distinct from propagation stock. PHAs that encompassed significantly with larger areas also reflected greater effective population sizes (µNE = 43; P<0.007). We juxtaposed these data against previously published results for prairie chicken and quail, and found population genetic structure driven by drift, habitat/climate impacts, and gender-biased selection via hunter-harvest. Each species (hunter-harvested or T&E) is independently managed, yet their composite population genetic baseline provides the quantitative criteria needed for an upland game bird MTM. Its implementation would require agricultural plots to be rehabilitated/reclaimed using a land-sharing/sparing portfolio that differs markedly from the Conservation Reserve Program (CRP), where sequestered land decreases as agricultural prices escalate. Cost-savings for an MTM would accrue by synchronizing single-species management with a dwindling hunter-harvest program, and by eliminating propagation/stocking programs. This would sustain not only native grasslands and their resident species, but also accelerate conservation at the wildlife-agroecosystem interface.

Unraveling historical introgression and resolving phylogenetic discord within Catostomus (Osteichthys: Catostomidae).

BACKGROUND: Porous species boundaries can be a source of conflicting hypotheses, particularly when coupled with variable data and/or methodological approaches. Their impacts can often be magnified when non-model organisms with complex histories of reticulation are investigated. One such example is the genus Catostomus (Osteichthys, Catostomidae), a freshwater fish clade with conflicting morphological and mitochondrial phylogenies. The former is hypothesized as reflecting the presence of admixed genotypes within morphologically distinct lineages, whereas the latter is interpreted as the presence of distinct morphologies that emerged multiple times through convergent evolution. We tested these hypotheses using multiple methods, to including multispecies coalescent and concatenated approaches. Patterson's D-statistic was applied to resolve potential discord, examine introgression, and test the putative hybrid origin of two species. We also applied naïve binning to explore potential effects of concatenation. RESULTS: We employed 14,007 loci generated from ddRAD sequencing of 184 individuals to derive the first highly supported nuclear phylogeny for Catostomus. Our phylogenomic analyses largely agreed with a morphological interpretation,with the exception of the placement of Xyrauchen texanus, which differs from both morphological and mitochondrial phylogenies. Additionally, our evaluation of the putative hybrid species C. columbianus revealed a lack introgression and instead matched the mitochondrial phylogeny. Furthermore, D-statistic tests clarified all discrepancies based solely on mitochondrial data, with agreement among topologies derived from concatenation and multispecies coalescent approaches. Extensive historic introgression was detected across six species-pairs. Potential endemism in the Virgin and Little Colorado Rivers was also apparent, and the former genus Pantosteus was derived as monophyletic, save for C. columbianus. CONCLUSIONS: Complex reticulated histories detected herein support the hypothesis that introgression was responsible for conflicts that occurred within the mitochondrial phylogeny, and explains discrepancies found between it and previous morphological phylogenies. Additionally, the hybrid origin of C. columbianus was refuted, but with the caveat that more fine-grain sampling is still needed. Our diverse phylogenomic approaches provided largely concordant results, with naïve binning useful in exploring the single conflict. Considerable diversity was found within Catostomus across southwestern North America, with two drainages [Virgin River (UT) and Little Colorado River (AZ)] reflecting unique composition.

Sequence-based molecular phylogenetics and phylogeography of the American box turtles (Terrapene spp.) with support from DNA barcoding.

The classification of the American box turtles (Terrapene spp.) has remained enigmatic to systematists. Previous comprehensive phylogenetic studies focused primarily on morphology. The goal of this study was to re-assess the classification of Terrapene spp. by obtaining DNA sequence data from a broad geographic range and from all four recognized species and 11 subspecies within the genus. Tissue samples were obtained for all taxa except for Terrapene nelsoni klauberi. DNA was extracted, and the mitochondrial DNA (mtDNA) cytochrome b (Cytb) and nuclear DNA (nucDNA) glyceraldehyde-3-phosphate-dehydrogenase (GAPD) genes were amplified via polymerase chain reaction and sequenced. In addition, the mtDNA gene commonly used for DNA barcoding (cytochrome oxidase c subunit I; COI) was amplified and sequenced to calculate pairwise percent DNA sequence divergence comparisons for each Terrapene taxon. The sequence data were analyzed using maximum likelihood and Bayesian phylogenetic inference, a molecular clock, AMOVAs, SAMOVAs, haplotype networks, and pairwise percent sequence divergence comparisons. Terrapene carolina mexicana and T. c. yucatana formed a monophyletic clade with T. c. triunguis, and this clade was paraphyletic to the rest of T. carolina. Terrapene ornata ornata and T. o. luteola lacked distinction phylogenetically, and Terrapene nelsoni was confirmed to be the sister taxon of T. ornata. Terrapene c. major, T. c. bauri, and Terrapene coahuila were not well resolved for some of the analyses. The DNA barcoding results indicated that all taxa were different species (>2% sequence divergence) except for T. c. triunguis - T. c. mexicana and T. o. ornata - T. o. luteola. The results suggest that T. c. triunguis should be elevated to species status (Terrapene mexicana), and mexicana and yucatana should be included in this group as subspecies. In addition, T. o. ornata and T. o. luteola should not be considered separate subspecies. The DNA barcoding data support these recommended taxonomic revisions. Because conservation efforts are typically species-based, these results will be important for facilitating successful conservation management strategies.