skip to main content

This content will become publicly available on January 1, 2023

Title: Genome-Scale Data Reveal Deep Lineage Divergence and a Complex Demographic History in the Texas Horned Lizard ( Phrynosoma cornutum ) throughout the Southwestern and Central United States
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 that more » 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
; ; ; ; ; ; ; ; ; ;
Baldauf, Sandra
Award ID(s):
1929679 2023723
Publication Date:
Journal Name:
Genome Biology and Evolution
Sponsoring Org:
National Science Foundation
More Like this
  1. 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 treemore »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 (but not preventing) gene flow and through parapatric differentiation along an ecological gradient.« 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 providemore »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 investigate 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 New study systems and tools are needed to understand how divergence and speciation occur between lineages with gene flow. Migratory birds often exhibit divergence despite seasonal migration, which brings populations into contact with one another. We studied divergence between 2 subspecies of Northern Saw-whet Owl (Aegolius acadicus), in which a sedentary population on the islands of Haida Gwaii, British Columbia (A. a. brooksi), exists in the presence of the other form (A. a. acadicus) during migration but not during the breeding season. Prior research showed fixed mtDNA divergence but left open the question of nuclear gene flow. We usedmore »2,517 ultraconserved element loci to examine the demographic history of this young taxon pair. Although we did not observe fixed single nucleotide polymorphism differences between populations among our genotyped individuals, 100% of the birds were diagnosable and δaδI analyses suggested the demographic model best fitting the data was one of split-bidirectional-migration (i.e. speciation with gene flow). We dated the split between brooksi and acadicus to ~278 Kya, and our analyses suggested gene flow between groups was skewed, with ~0.7 individuals per generation coming from acadicus into brooksi and ~4.4 going the opposite direction. Coupled with an absence of evidence of phenotypic hybrids and the birds’ natural history, these data suggest brooksi may be a young biological species arising despite historic gene flow.« less
  4. Kubatko, Laura (Ed.)
    Abstract Evidence from natural systems suggests that hybridization between animal species is more common than traditionally thought, but the overall contribution of introgression to standing genetic variation within species remains unclear for most animal systems. Here, we use targeted exon capture to sequence thousands of nuclear loci and complete mitochondrial genomes from closely related chipmunk species in the Tamias quadrivittatus group that are distributed across the Great Basin and the central and southern Rocky Mountains of North America. This recent radiation includes six overlapping, ecologically distinct species (Tamias canipes, Tamias cinereicollis, Tamias dorsalis, T. quadrivittatus, Tamias rufus, and Tamias umbrinus)more »that show evidence for widespread introgression across species boundaries. Such evidence has historically been derived from a handful of markers, typically focused on mitochondrial loci, to describe patterns of introgression; consequently, the extent of introgression of nuclear genes is less well characterized. We conducted a series of phylogenomic and species-tree analyses to resolve the phylogeny of six species in this group. In addition, we performed several population-genomic analyses to characterize nuclear genomes and infer coancestry among individuals. Furthermore, we used emerging quartets-based approaches to simultaneously infer the species tree (SVDquartets) and identify introgression (HyDe). We found that, in spite of rampant introgression of mitochondrial genomes between some species pairs (and sometimes involving up to three species), there appears to be little to no evidence for nuclear introgression. These findings mirror other genomic results where complete mitochondrial capture has occurred between chipmunk species in the absence of appreciable nuclear gene flow. The underlying causes of recurrent massive cytonuclear discordance remain unresolved in this group but mitochondrial DNA appears highly misleading of population histories as a whole. Collectively, it appears that chipmunk species boundaries are largely impermeable to nuclear gene flow and that hybridization, while pervasive with respect to mtDNA, has likely played a relatively minor role in the evolutionary history of this group. [Cytonuclear discordance; hyridization; introgression, phylogenomics; SVDquartets; Tamias.]« less
  5. Aim: To review the histories of the Colorado River and North American monsoon system to ascertain their effects on the genetic divergence of desert‐adapted animals. Location: Lower Colorado River region, including Mojave and Sonoran deserts, United States. Methods: We synthesized recent geological literature to summarize initiation phases of lower Colorado River evolution, their discrepancies, and potential for post‐vicariance dispersal of animals across the river. We simulated data under geological models and performed a meta‐analysis of published and unpublished genetic data including population diversity metrics, relatedness and historical migration rates to assess alternative divergence hypotheses. Results: The two models for arrivalmore »of the Colorado River into the Gulf of California impose east‐west divergence ages of 5.3 and 4.8 Ma, respectively. We found quantifiable river‐associated differentiation in the lower Colorado River region in reptiles, arachnids and mammals relative to flying insects. However, topological statistics, historical migration rates and cross‐river extralimital populations suggest that the river should be considered a leaky barrier that filters, rather than prevents, gene flow. Most markers violated neutrality tests. Differential adaptation to monsoon‐based precipitation differences may contribute to divergence between Mojave and Sonoran populations and should be tested. Main Conclusions: Rivers are dynamic features that can both limit and facilitate gene flow through time, the impacts of which are mitigated by species‐specific life history and dispersal traits. The Southwest is a geo‐climatically complex region with the potential to produce pseudocongruent patterns of genetic divergence, offering a good setting to evaluate intermediate levels of geological‐biological (geobiological) complexity.« less