A community genomics approach to natural hybridization.

Hybridization is a complicated, oft-misunderstood process. Once deemed unnatural and uncommon, hybridization is now recognized as ubiquitous among species. But hybridization rates within and among communities are poorly understood despite the relevance to ecology, evolution and conservation. To clarify, we examined hybridization across 75 freshwater fish communities within the Ozarks of the North American Interior Highlands (USA) by single nucleotide polymorphism (SNP) genotyping 33 species (N = 2865 individuals; double-digest restriction site-associated DNA sequencing (ddRAD)). We found evidence of hybridization (70 putative hybrids; 2.4% of individuals) among 18 species-pairs involving 73% (24/33) of study species, with the majority being concentrated within one family (Leuciscidae/minnows; 15 species; 66 hybrids). Interspecific genetic exchange-or introgression-was evident from 24 backcrossed individuals (10/18 species-pairs). Hybrids occurred within 42 of 75 communities (56%). Four selected environmental variables (species richness, protected area extent, precipitation (May and annually)) exhibited 73-78% accuracy in predicting hybrid occurrence via random forest classification. Our community-level assessment identified hybridization as spatially widespread and environmentally dependent (albeit predominantly within one diverse, omnipresent family). Our approach provides a more holistic survey of natural hybridization by testing a wide range of species-pairs, thus contrasting with more conventional evaluations.

Riverscape community genomics: A comparative analytical approach to identify common drivers of spatial structure.

Genetic differentiation among local groups of individuals, that is, genetic ß-diversity, is a key component of population persistence related to connectivity and isolation. However, most genetic investigations of natural populations focus on a single species, overlooking opportunities for multispecies conservation plans to benefit entire communities in an ecosystem. We present an approach to evaluate genetic ß-diversity within and among many species and demonstrate how this riverscape community genomics approach can be applied to identify common drivers of genetic structure. Our study evaluated genetic ß-diversity in 31 co-distributed native stream fishes sampled from 75 sites across the White River Basin (Ozarks, USA) using SNP genotyping (ddRAD). Despite variance among species in the degree of genetic divergence, general spatial patterns were identified corresponding to river network architecture. Most species (N = 24) were partitioned into discrete subpopulations (K = 2-7). We used partial redundancy analysis to compare species-specific genetic ß-diversity across four models of genetic structure: Isolation by distance (IBD), isolation by barrier (IBB), isolation by stream hierarchy (IBH), and isolation by environment (IBE). A significant proportion of intraspecific genetic variation was explained by IBH (x̄ = 62%), with the remaining models generally redundant. We found evidence for consistent spatial modularity in that gene flow is higher within rather than between hierarchical units (i.e., catchments, watersheds, basins), supporting the generalization of the stream hierarchy model. We discuss our conclusions regarding conservation and management and identify the 8-digit hydrologic unit (HUC) as the most relevant spatial scale for managing genetic diversity across riverine networks.

ClineHelpR: an R package for genomic cline outlier detection and visualization.

BACKGROUND: Patterns of multi-locus differentiation (i.e., genomic clines) often extend broadly across hybrid zones and their quantification can help diagnose how species boundaries are shaped by adaptive processes, both intrinsic and extrinsic. In this sense, the transitioning of loci across admixed individuals can be contrasted as a function of the genome-wide trend, in turn allowing an expansion of clinal theory across a much wider array of biodiversity. However, computational tools that serve to interpret and consequently visualize 'genomic clines' are limited, and users must often write custom, relatively complex code to do so. RESULTS: Here, we introduce the ClineHelpR R-package for visualizing genomic clines and detecting outlier loci using output generated by two popular software packages, bgc and Introgress. ClineHelpR bundles both input generation (i.e., filtering datasets and creating specialized file formats) and output processing (e.g., MCMC thinning and burn-in) with functions that directly facilitate interpretation and hypothesis testing. Tools are also provided for post-hoc analyses that interface with external packages such as ENMeval and RIdeogram. CONCLUSIONS: Our package increases the reproducibility and accessibility of genomic cline methods, thus allowing an expanded user base and promoting these methods as mechanisms to address diverse evolutionary questions in both model and non-model organisms. Furthermore, the ClineHelpR extended functionality can evaluate genomic clines in the context of spatial and environmental features, allowing users to explore underlying processes potentially contributing to the observed patterns and helping facilitate effective conservation management strategies.

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 .

Contrasting signatures of introgression in North American box turtle (Terrapene spp.) contact zones.

Hybridization occurs differentially across the genome in a balancing act between selection and migration. With the unprecedented resolution of contemporary sequencing technologies, selection and migration can now be effectively quantified such that researchers can identify genetic elements involved in introgression. Furthermore, genomic patterns can now be associated with ecologically relevant phenotypes, given availability of annotated reference genomes. We do so in North American box turtles (Terrapene) by deciphering how selection affects hybrid zones at the interface of species boundaries and identifying genetic regions potentially under selection that may relate to thermal adaptations. Such genes may impact physiological pathways involved in temperature-dependent sex determination, immune system functioning and hypoxia tolerance. We contrasted these patterns across inter- and intraspecific hybrid zones that differ temporally and biogeographically. We demonstrate hybridization is broadly apparent in Terrapene, but with observed genomic cline patterns corresponding to species boundaries at loci potentially associated with thermal adaptation. These loci display signatures of directional introgression within intraspecific boundaries, despite a genome-wide selective trend against intergrades. In contrast, outlier loci for interspecific comparisons exhibited evidence of being under selection against hybrids. Importantly, adaptations coinciding with species boundaries in Terrapene overlap with climatic boundaries and highlight the vulnerability of these terrestrial ectotherms to anthropogenic pressures.

MrBait: universal identification and design of targeted-enrichment capture probes.

Motivation: It is a non-trivial task to identify and design capture probes ('baits') for the diverse array of targeted-enrichment methods now available (e.g. ultra-conserved elements, anchored hybrid enrichment, RAD-capture). This often involves parsing large genomic alignments, followed by multiple steps of curating candidate genomic regions to optimize targeted information content (e.g. genetic variation) and to minimize potential probe dimerization and non-target enrichment. Results: In this context, we developed MrBait, a user-friendly, generalized software pipeline for identification, design and optimization of targeted-enrichment probes across a range of target-capture paradigms. MrBait is an open-source codebase that leverages native parallelization capabilities in Python and mitigates memory usage via a relational-database back-end. Numerous filtering methods allow comprehensive optimization of designed probes, including built-in functionality that employs BLAST, similarity-based clustering and a graph-based algorithm that 'rescues' failed probes. Availability and implementation: Complete code for MrBait is available on GitHub (https://github.com/tkchafin/mrbait), and is also available with all dependencies via one-line installation using the conda package manager. Online documentation describing installation and runtime instructions can be found at: https://mrbait.readthedocs.io. Supplementary information: Supplementary data are available at Bioinformatics online.

Hybridization drives genetic erosion in sympatric desert fishes of western North America.

Many species have evolved or currently coexist in sympatry due to differential adaptation in a heterogeneous environment. However, anthropogenic habitat modifications can either disrupt reproductive barriers or obscure environmental conditions which underlie fitness gradients. In this study, we evaluated the potential for an anthropogenically-mediated shift in reproductive boundaries that separate two historically sympatric fish species (Gila cypha and G. robusta) endemic to the Colorado River Basin using ddRAD sequencing of 368 individuals. We first examined the integrity of reproductive isolation while in sympatry and allopatry, then characterized hybrid ancestries using genealogical assignment tests. We tested for localized erosion of reproductive isolation by comparing site-wise genomic clines against global patterns and identified a breakdown in the drainage-wide pattern of selection against interspecific heterozygotes. This, in turn, allowed for the formation of a hybrid swarm in one tributary, and asymmetric introgression where species co-occur. We also detected a weak but significant relationship between genetic purity and degree of consumptive water removal, suggesting a role for anthropogenic habitat modifications in undermining species boundaries or expanding historically limited introgression. In addition, results from basin-wide genomic clines suggested that hybrids and parental forms are adaptively nonequivalent. If so, then a failure to manage for hybridization will exacerbate the long-term extinction risk in parental populations. These results reinforce the role of anthropogenic habitat modification in promoting interspecific introgression in sympatric species by relaxing divergent selection. This, in turn, underscores a broader role for hybridization in decreasing global biodiversity within rapidly deteriorating environments.

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.

Do biofilm communities respond to the chemical signatures of fracking? A test involving streams in North-central Arkansas.

BACKGROUND: Unconventional natural gas (UNG) extraction (fracking) is ongoing in 29 North American shale basins (20 states), with ~6000 wells found within the Fayetteville shale (north-central Arkansas). If the chemical signature of fracking is detectable in streams, it can be employed to bookmark potential impacts. We evaluated benthic biofilm community composition as a proxy for stream chemistry so as to segregate anthropogenic signatures in eight Arkansas River catchments. In doing so, we tested the hypothesis that fracking characteristics in study streams are statistically distinguishable from those produced by agriculture or urbanization. RESULTS: Four tributary catchments had UNG-wells significantly more dense and near to our sampling sites and were grouped as 'potentially-impacted catchment zones' (PICZ). Four others were characterized by significantly larger forested area with greater slope and elevation but reduced pasture, and were classified as 'minimally-impacted' (MICZ). Overall, 46 bacterial phyla/141 classes were identified, with 24 phyla (52%) and 54 classes (38%) across all samples. PICZ-sites were ecologically more variable than MICZ-sites, with significantly greater nutrient levels (total nitrogen, total phosphorous), and elevated Cyanobacteria as bioindicators that tracked these conditions. PICZ-sites also exhibited elevated conductance (a correlate of increased ion concentration) and depressed salt-intolerant Spartobacteria, suggesting the presence of brine as a fracking effect. Biofilm communities at PICZ-sites were significantly less variable than those at MICZ-sites. CONCLUSIONS: Study streams differed by Group according to morphology, land use, and water chemistry but not in biofilm community structure. Those at PICZ-sites covaried according to anthropogenic impact, and were qualitatively similar to communities found at sites disturbed by fracking. The hypothesis that fracking signatures in study streams are distinguishable from those produced by other anthropogenic effects was statistically rejected. Instead, alterations in biofilm community composition, as induced by fracking, may be less specific than initially predicted, and thus more easily confounded by agriculture and urbanization effects (among others). Study streams must be carefully categorized with regard to the magnitude and extent of anthropogenic impacts. They must also be segregated with statistical confidence (as herein) before fracking impacts are monitored.

Stream hierarchy defines riverscape genetics of a North American desert fish.

Global climate change is apparent within the Arctic and the south-western deserts of North America, with record drought in the latter reflected within 640,000 km(2) of the Colorado River Basin. To discern the manner by which natural and anthropogenic drivers have compressed Basin-wide fish biodiversity, and to establish a baseline for future climate effects, the Stream Hierarchy Model (SHM) was employed to juxtapose fluvial topography against molecular diversities of 1092 Bluehead Sucker (Catostomus discobolus). MtDNA revealed three geomorphically defined evolutionarily significant units (ESUs): Bonneville Basin, upper Little Colorado River and the remaining Colorado River Basin. Microsatellite analyses (16 loci) reinforced distinctiveness of the Bonneville Basin and upper Little Colorado River, but subdivided the Colorado River Basin into seven management units (MUs). One represents a cline of three admixed gene pools comprising the mainstem and its lower-gradient tributaries. Six others are not only distinct genetically but also demographically (i.e. migrants/generation <9.7%). Two of these (i.e. Grand Canyon and Canyon de Chelly) are defined by geomorphology, two others (i.e. Fremont-Muddy and San Raphael rivers) are isolated by sharp declivities as they drop precipitously from the west slope into the mainstem Colorado/Green rivers, another represents an isolated impoundment (i.e. Ringdahl Reservoir), while the last corresponds to a recognized subspecies (i.e. Zuni River, NM). Historical legacies of endemic fishes (ESUs) and their evolutionary potential (MUs) are clearly represented in our data, yet their arbiter will be the unrelenting natural and anthropogenic water depletions that will precipitate yet another conservation conflict within this unique but arid region.

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.

The Himalayan uplift and evolution of aquatic biodiversity across Asia: Snowtrout (Cyprininae: Schizothorax) as a test case.

Global biodiversity hotspots are often remote, tectonically active areas undergoing climatic fluctuations, such as the Himalaya Mountains and neighboring Qinghai-Tibetan Plateau (QTP). They provide biogeographic templates upon which endemic biodiversity can be mapped to infer diversification scenarios. Yet, this process can be somewhat opaque for the Himalaya, given substantial data gaps separating eastern and western regions. To help clarify, we evaluated phylogeographic and phylogenetic hypotheses for a widespread fish (Snowtrout: Cyprininae; Schizothorax) by sequencing 1,140 base pair of mtDNA cytochrome-b (cytb) from Central Himalaya samples (Nepal: N = 53; Bhutan: N = 19), augmented with 68 GenBank sequences (N = 60 Schizothorax/N = 8 outgroups). Genealogical relationships (N = 132) were analyzed via maximum likelihood (ML), Bayesian (BA), and haplotype network clustering, with clade divergence estimated via TimeTree. Snowtrout seemingly originated in Central Asia, dispersed across the QTP, then into Bhutan via southward-flowing tributaries of the east-flowing Yarlung-Tsangpo River (YLTR). Headwaters of five large Asian rivers provided dispersal corridors from Central into eastern/southeastern Asia. South of the Himalaya, the YLTR transitions into the Brahmaputra River, facilitating successive westward colonization of Himalayan drainages first in Bhutan, then Nepal, followed by far-western drainages subsequently captured by the (now) westward-flowing Indus River. Two distinct Bhutanese phylogenetic groups were recovered: Bhutan-1 (with three subclades) seemingly represents southward dispersal from the QTP; Bhutan-2 apparently illustrates northward colonization from the Lower Brahmaputra. The close phylogenetic/phylogeographic relationships between the Indus River (Pakistan) and western tributaries of the Upper Ganges (India/Nepal) potentially implicate an historic, now disjunct connection. Greater species-divergences occurred across rather than within-basins, suggesting vicariance as a driver. The Himalaya is a component of the Earth's largest glacial reservoir (i.e., the "third-pole") separate from the Arctic/Antarctic. Its unique aquatic biodiversity must be defined and conserved through broad, trans-national collaborations. Our study provides an initial baseline for this process.

Linking extinction-colonization dynamics to genetic structure in a salamander metapopulation.

Theory predicts that founder effects have a primary role in determining metapopulation genetic structure. However, ecological factors that affect extinction-colonization dynamics may also create spatial variation in the strength of genetic drift and migration. We tested the hypothesis that ecological factors underlying extinction-colonization dynamics influenced the genetic structure of a tiger salamander (Ambystoma tigrinum) metapopulation. We used empirical data on metapopulation dynamics to make a priori predictions about the effects of population age and ecological factors on genetic diversity and divergence among 41 populations. Metapopulation dynamics of A. tigrinum depended on wetland area, connectivity and presence of predatory fish. We found that newly colonized populations were more genetically differentiated than established populations, suggesting that founder effects influenced genetic structure. However, ecological drivers of metapopulation dynamics were more important than age in predicting genetic structure. Consistent with demographic predictions from metapopulation theory, genetic diversity and divergence depended on wetland area and connectivity. Divergence was greatest in small, isolated wetlands where genetic diversity was low. Our results show that ecological factors underlying metapopulation dynamics can be key determinants of spatial genetic structure, and that habitat area and isolation may mediate the contributions of drift and migration to divergence and evolution in local populations.

Limited gene flow and pronounced population genetic structure of Eastern Massasauga (Sistrurus catenatus) in a Midwestern prairie remnant.

As anthropogenic changes continue to ecologically stress wildlife, obtaining measures of gene flow and genetic diversity are crucial for evaluating population trends and considering management and conservation strategies for small, imperiled populations. In our study, we conducted a molecular assessment to expand on previous work to elucidate patterns of diversity and connectivity in the remaining disjunct Eastern Massasauga Rattlesnake (Sistrurus catenatus) hibernacula in Illinois. We assayed genetic data for 327 samples collected during 1999-2015 from the Carlyle Lake study area across 21 microsatellite loci. We found hibernacula formed distinct genetic clusters corresponding to the three main study areas (Dam Recreation Areas, Eldon Hazlet State Park, and South Shore State Park). Genetic structuring and low estimates of dispersal indicated that connectivity among these study areas is limited and each is demographically independent. Hibernacula exhibited moderate levels of heterozygosity (0.60-0.73), but estimates of effective population size (5.2-41.0) were low and track census sizes generated via long-term mark-recapture data. Hibernacula at Carlyle Lake, which represent the only Eastern Massasauga remaining in Illinois, are vulnerable to future loss of genetic diversity through lack of gene flow as well as demographic and environmental stochastic processes. Our work highlights the need to include population-level genetic data in recovery planning and suggests that recovery efforts should focus on managing the three major study areas as separate conservation units in order to preserve and maintain long-term adaptive potential of these populations. Specific management goals should include improving connectivity among hibernacula, maintaining existing wet grassland habitat, and minimizing anthropogenic sources of mortality caused by habitat management (e.g., mowing, prescribed fire) and recreational activities. Our molecular study provides additional details about demographic parameters and connectivity at Carlyle Lake that can be used to guide recovery of Eastern Massasauga in Illinois and throughout its range.

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.

The choices we make and the impacts they have: Machine learning and species delimitation in North American box turtles (Terrapene spp.).

Model-based approaches that attempt to delimit species are hampered by computational limitations as well as the unfortunate tendency by users to disregard algorithmic assumptions. Alternatives are clearly needed, and machine-learning (M-L) is attractive in this regard as it functions without the need to explicitly define a species concept. Unfortunately, its performance will vary according to which (of several) bioinformatic parameters are invoked. Herein, we gauge the effectiveness of M-L-based species-delimitation algorithms by parsing 64 variably-filtered versions of a ddRAD-derived SNP data set collected from North American box turtles (Terrapene spp.). Our filtering strategies included: (i) minor allele frequencies (MAF) of 5%, 3%, 1%, and 0% (= none), and (ii) maximum missing data per-individual/per-population at 25%, 50%, 75%, and 100% (= no filtering). We found that species-delimitation via unsupervised M-L impacted the signal-to-noise ratio in our data, as well as the discordance among resolved clades. The latter may also reflect biogeographic history, gene flow, incomplete lineage sorting, or combinations thereof (as corroborated from previously observed patterns of differential introgression). Our results substantiate M-L as a viable species-delimitation method, but also demonstrate how commonly observed patterns of phylogenetic discordance can seriously impact M-L-classification.

Spatial population genetics in heavily managed species: Separating patterns of historical translocation from contemporary gene flow in white-tailed deer.

Approximately 100 years ago, unregulated harvest nearly eliminated white-tailed deer (Odocoileus virginianus) from eastern North America, which subsequently served to catalyze wildlife management as a national priority. An extensive stock-replenishment effort soon followed, with deer broadly translocated among states as a means of re-establishment. However, an unintended consequence was that natural patterns of gene flow became obscured and pretranslocation signatures of population structure were replaced. We applied cutting-edge molecular and biogeographic tools to disentangle genetic signatures of historical management from those reflecting spatially heterogeneous dispersal by evaluating 35,099 single nucleotide polymorphisms (SNPs) derived via reduced-representation genomic sequencing from 1143 deer sampled statewide in Arkansas. We then employed Simpson's diversity index to summarize ancestry assignments and visualize spatial genetic transitions. Using sub-sampled transects across these transitions, we tested clinal patterns across loci against theoretical expectations of their response under scenarios of re-colonization and restricted dispersal. Two salient results emerged: (A) Genetic signatures from historic translocations are demonstrably apparent; and (B) Geographic filters (major rivers; urban centers; highways) now act as inflection points for the distribution of this contemporary ancestry. These results yielded a statewide assessment of contemporary population structure in deer as driven by historic translocations as well as ongoing processes. In addition, the analytical framework employed herein to effectively decipher extant/historic drivers of deer distribution in Arkansas is also applicable for other biodiversity elements with similarly complex demographic histories.

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.

Parallel introgression, not recurrent emergence, explains apparent elevational ecotypes of polyploid Himalayan snowtrout.

The recurrence of similar evolutionary patterns within different habitats often reflects parallel selective pressures acting upon either standing or independently occurring genetic variation to produce a convergence of phenotypes. This interpretation (i.e. parallel divergences within adjacent streams) has been hypothesized for drainage-specific morphological 'ecotypes' observed in polyploid snowtrout (Cyprinidae: Schizothorax). However, parallel patterns of differential introgression during secondary contact are a viable alternative hypothesis. Here, we used ddRADseq (N = 35 319 de novo and N = 10 884 transcriptome-aligned SNPs), as derived from Nepali/Bhutanese samples (N = 48 each), to test these competing hypotheses. We first employed genome-wide allelic depths to derive appropriate ploidy models, then a Bayesian approach to yield genotypes statistically consistent under the inferred expectations. Elevational 'ecotypes' were consistent in geometric morphometric space, but with phylogenetic relationships at the drainage level, sustaining a hypothesis of independent emergence. However, partitioned analyses of phylogeny and admixture identified subsets of loci under selection that retained genealogical concordance with morphology, suggesting instead that apparent patterns of morphological/phylogenetic discordance are driven by widespread genomic homogenization. Here, admixture occurring in secondary contact effectively 'masks' previous isolation. Our results underscore two salient factors: (i) morphological adaptations are retained despite hybridization and (ii) the degree of admixture varies across tributaries, presumably concomitant with underlying environmental or anthropogenic factors.

Conservation phylogenetics of helodermatid lizards using multiple molecular markers and a supertree approach.

We analyzed both mitochondrial (mt-) and nuclear (n) DNAs in a conservation phylogenetic framework to examine deep and shallow histories of the Beaded Lizard (Heloderma horridum) and Gila Monster (H. suspectum) throughout their geographic ranges in North and Central America. Both mtDNA and intron markers clearly partitioned each species. One intron and mtDNA further subdivided H. horridum into its four recognized subspecies (H. n. alvarezi, charlesbogerti,exasperatum, and horridum). However, the two subspecies of H. suspectum (H. s. suspectum and H. s. cinctum) were undefined. A supertree approach sustained these relationships. Overall, the Helodermatidae is reaffirmed as an ancient and conserved group. Its most recent common ancestor (MRCA) was Lower Eocene [35.4 million years ago (mya)], with a approximately 25 my period of stasis before the MRCA of H. horridum diversified in Lower Miocene. Another approximately 5 my passed before H. h. exasperatum and H. h. horridum diverged, followed by approximately 1.5 my before H. h. alvarezi and H. h. charlesbogerti separated. Heloderma suspectum reflects an even longer period of stasis (approximately 30 my) before diversifying from its MRCA. Both H. suspectum (México) and H. h. alvarezi also revealed evidence of historic range expansion following a recent bottleneck. Our conservation phylogenetic approach emphasizes the origin and diversification of this group, yields information on the manner by which past environmental variance may have impacted its populations and, in turn, allows us to disentangle historic from contemporary impacts that might threaten its long-term persistence. The value of helodermatid conservation resides in natural services and medicinal products, particularly venom constituents, and these are only now being realized.