skip to main content

Title: Phylogenomic Assessment of Biodiversity Using a Reference-Based Taxonomy: An Example With Horned Lizards (Phrynosoma)
Phylogenomic investigations of biodiversity facilitate the detection of fine-scale population genetic structure and the demographic histories of species and populations. However, determining whether or not the genetic divergence measured among populations reflects species-level differentiation remains a central challenge in species delimitation. One potential solution is to compare genetic divergence between putative new species with other closely related species, sometimes referred to as a reference-based taxonomy. To be described as a new species, a population should be at least as divergent as other species. Here, we develop a reference-based taxonomy for Horned Lizards ( Phrynosoma ; 17 species) using phylogenomic data (ddRADseq data) to provide a framework for delimiting species in the Greater Short-horned Lizard species complex ( P. hernandesi ). Previous species delimitation studies of this species complex have produced conflicting results, with morphological data suggesting that P. hernandesi consists of five species, whereas mitochondrial DNA support anywhere from 1 to 10 + species. To help address this conflict, we first estimated a time-calibrated species tree for P. hernandesi and close relatives using SNP data. These results support the paraphyly of P. hernandesi; we recommend the recognition of two species to promote a taxonomy that is consistent with species monophyly. more » There is strong evidence for three populations within P. hernandesi , and demographic modeling and admixture analyses suggest that these populations are not reproductively isolated, which is consistent with previous morphological analyses that suggest hybridization could be common. Finally, we characterize the population-species boundary by quantifying levels of genetic divergence for all 18 Phrynosoma species. Genetic divergence measures for western and southern populations of P. hernandesi failed to exceed those of other Phrynosoma species, but the relatively small population size estimated for the northern population causes it to appear as a relatively divergent species. These comparisons underscore the difficulties associated with putting a reference-based approach to species delimitation into practice. Nevertheless, the reference-based approach offers a promising framework for the consistent assessment of biodiversity within clades of organisms with similar life histories and ecological traits. « less
; ; ; ; ;
Award ID(s):
Publication Date:
Journal Name:
Frontiers in Ecology and Evolution
Sponsoring Org:
National Science Foundation
More Like this
  1. Baldauf, Sandra (Ed.)
    Abstract The southwestern and central United States serve as an ideal region to test alternative hypotheses regarding biotic diversification. Genomic data can now be combined with sophisticated computational models to quantify the impacts of paleoclimate change, geographic features, and habitat heterogeneity on spatial patterns of genetic diversity. In this study, we combine thousands of genotyping-by-sequencing (GBS) loci with mtDNA sequences (ND1) from the Texas horned lizard (Phrynosoma cornutum) to quantify relative support for different catalysts of diversification. Phylogenetic and clustering analyses of the GBS data indicate support for at least three primary populations. The spatial distribution of populations appears concordant with habitat type, with desert populations in AZ and NM showing the largest genetic divergence from the remaining populations. The mtDNA data also support a divergent desert population, but other relationships differ and suggest mtDNA introgression. Genotype–environment association with bioclimatic variables supports divergence along precipitation gradients more than along temperature gradients. Demographic analyses support a complex history, with introgression and gene flow playing an important role during diversification. Bayesian multispecies coalescent analyses with introgression (MSci) analyses also suggest that gene flow occurred between populations. Paleo-species distribution models support two southern refugia that geographically correspond to contemporary lineages. We find thatmore »divergence times are underestimated and population sizes are overestimated when introgression occurred and is ignored in coalescent analyses, and furthermore, inference of ancient introgression events and demographic history is sensitive to inclusion of a single recently admixed sample. Our analyses cannot refute the riverine barrier or glacial refugia hypotheses. Results also suggest that populations are continuing to diverge along habitat gradients. Finally, the strong evidence of admixture, gene flow, and mtDNA introgression among populations suggests that P. cornutum should be considered a single widespread species under the General Lineage Species Concept.« less
  2. Smith, Stephen (Ed.)
    Abstract Understanding how gene flow affects population divergence and speciation remains challenging. Differentiating one evolutionary process from another can be difficult because multiple processes can produce similar patterns, and more than one process can occur simultaneously. Although simple population models produce predictable results, how these processes balance in taxa with patchy distributions and complicated natural histories is less certain. These types of populations might be highly connected through migration (gene flow), but can experience stronger effects of genetic drift and inbreeding, or localized selection. Although different signals can be difficult to separate, the application of high-throughput sequence data can provide the resolution necessary to distinguish many of these processes. We present whole-genome sequence data for an avian species group with an alpine and arctic tundra distribution to examine the role that different population genetic processes have played in their evolutionary history. Rosy-finches inhabit high elevation mountaintop sky islands and high-latitude island and continental tundra. They exhibit extensive plumage variation coupled with low levels of genetic variation. Additionally, the number of species within the complex is debated, making them excellent for studying the forces involved in the process of diversification, as well as an important species group in which to investigatemore »species boundaries. Total genomic variation suggests a broadly continuous pattern of allele frequency changes across the mainland taxa of this group in North America. However, phylogenomic analyses recover multiple distinct, well supported, groups that coincide with previously described morphological variation and current species-level taxonomy. Tests of introgression using D-statistics and approximate Bayesian computation reveal significant levels of introgression between multiple North American taxa. These results provide insight into the balance between divergent and homogenizing population genetic processes and highlight remaining challenges in interpreting conflict between different types of analytical approaches with whole-genome sequence data. [ABBA-BABA; approximate Bayesian computation; gene flow; phylogenomics; speciation; whole-genome sequencing.]« less
  3. Abstract

    Conservation efforts rely on robust taxonomic assessments that should be based on critical assessment of interspecific boundaries, infraspecific variation, and potentially distinctive peripheral populations. The meadow vole (Microtus pennsylvanicus) is widely distributed across North America, including 28 morphologically defined subspecies and numerous isolated populations. Because some subspecies are of high conservation concern, we examined geographic variation across the range of the species to test existing infraspecific taxonomy in terms of local and regional diversification. We sequenced mitochondrial DNA (mtDNA) from 20 subspecies of M. pennsylvanicus and contextualized infraspecific variation through comparison of pairwise genetic distances derived from an extended data set of 63 species of Microtus. We found strong support for at least three divergent clades within M. pennsylvanicus, with observed intraspecific clade divergence exceeding that between several pairwise comparisons of sister species within Microtus. Six nuclear genes were then sequenced to test the validity of mtDNA structure and to further evaluate the possibility of cryptic, species-level diversity using Bayes factor species delimitation (BFD) analyses. BFD consistently and decisively supported multiple species based on the multilocus approach. We propose that taxonomic revision of the meadow vole is required, with the eastern clade now identified as M. pennsylvanicus (Ord 1815),more »the western clade as M. drummondii (Audubon and Bachman 1853), and the coastal Florida clade as M. dukecampbelli (Woods, Post, and Kilpatrick 1982). We suggest that such an arrangement would more closely reflect evolutionary history and provide critical context for further examination of distinctive southern peripheral populations that harbor novel evolutionary legacies and adaptive potential.

    « less
  4. Sosa-Calvo, Jeffrey (Ed.)
    Abstract Using genetic, morphological, and geographical evidence, we investigate the species-level taxonomy and evolutionary history of the Pseudomyrmex elongatulus group, a clade of ants distributed from southwestern United States to Costa Rica. Through targeted enrichment of 2,524 UCE (ultraconserved element) loci we generate a phylogenomic data set and clarify the phylogenetic relationships and biogeographic history of these ants. The crown group is estimated to have originated ~8 Ma, in Mexico and/or northern Central America, and subsequently expanded into southern Central America and the southwestern Nearctic. The P. elongatulus group contains a mix of low- and high-elevation species, and there were apparently multiple transitions between these habitat types. We uncover three examples of one species—of restricted or marginal geographical distribution—being embedded phylogenetically in another species, rendering the latter paraphyletic. One of these cases involves an apparent workerless social parasite that occurs sympatrically with its parent species, with the latter serving as host. This suggests a sympatric origin of the parasite species within the distribution range of its host. Species boundaries are tested using three molecular delimitation approaches (SODA, bPTP, BPP) but these methods generate inflated species estimates (26–46 species), evidently because of a failure to distinguish population structure from species differences.more »In a formal taxonomic revision of the P. elongatulus group, based on almost 3,000 specimens from ~625 localities, we allow for geographic variation within species and we employ distinctness-in-sympatry criteria for testing hypotheses about species limits. Under these guidelines we recognize 13 species, of which nine are new: P. arcanus, sp. nov. (western Mexico); P. capillatus, sp. nov. (western Mexico); P. cognatus, sp. nov. (Chiapas, Mexico to Nicaragua); P. comitator, sp. nov. (Chiapas, Mexico); P. ereptor, sp. nov. (Veracruz, Mexico); P. exoratus, sp. nov. (southeastern Mexico, Honduras); P. fasciatus, sp. nov. (Chiapas, Mexico to Costa Rica); P. nimbus, sp. nov. (Costa Rica); and P. veracruzensis, sp. nov. (Veracruz, Mexico). Our study highlights the value of combining phylogenomic, phenotypic, and geographical data to resolve taxonomic and evolutionary questions.« less
  5. Identifying the evolutionary and ecological mechanisms that drive lineage diversification in the species-rich tropics is of broad interest to evolutionary biologists. Here, we use phylogeographic and demographic analyses of genomic scale RADseq data to assess the impact of a large geographic feature, the Amazon River, on lineage formation in a venomous pitviper, Bothrops atrox. We compared genetic differentiation in samples from four sites near Santarem, Brazil that spanned the Amazon and represented major habitat types. A species delimitation analysis identified each population as a distinct evolutionary lineage while a species tree analysis with populations as taxa revealed a phylogenetic tree consistent with dispersal across the Amazon from north to south. Phylogenetic analyses of mtDNA variation confirmed this pattern and suggest that all lineages originated during the mid- to late-Pleistocene. Historical demographic analyses support a population model of lineage formation through isolation between lineages with low ongoing migration between large populations and reject a model of differentiation through isolation by distance alone. Our results provide a rare example of a phylogeographic pattern demonstrating dispersal over evolutionary time scales across a large tropical river and suggest a role for the Amazon River as a driver of in-situ divergence by both impeding (butmore »not preventing) gene flow and through parapatric differentiation along an ecological gradient.« less