skip to main content
US FlagAn official website of the United States government
dot gov icon
Official websites use .gov
A .gov website belongs to an official government organization in the United States.
https lock icon
Secure .gov websites use HTTPS
A lock ( lock ) or https:// means you've safely connected to the .gov website. Share sensitive information only on official, secure websites.


Title: Contrasting the depths of divergence between gene-tree and coalescent estimates in the North American racers (Colubridae: Coluber constrictor )
Abstract The North American racers (Coluber constrictor) are widely distributed across the Nearctic and numerous studies have demonstrated extensive variation in morphology, ecology, and population genetic structure. Here we take an integrative approach to understand lineage diversification within this taxon by combining genomic sequence capture data, mtDNA sequence data, morphometrics, and ecological niche models. Both the genomic data and mtDNA phylogeographic analyses support five lineages distributed across the range of this species. However, demographic model selection based on these two datasets strongly conflict in both the model of divergence and estimates of timing of lineage divergence. While mtDNA and concatenated genomic data suggest a Miocene origin of these distinct groups, coalescent-based demographic models with the sequence capture data suggest lineage diversification occurred at ~33 kya in allopatry without gene flow. Using linear morphological measurements of head shape we demonstrate that lineages distributed largely east and west of the Mississippi River are distinguishable. Furthermore, ecological niche models demonstrate that lineages distributed in subtropical habitats have environmental niche space that is significantly differentiated from lineages distributed across the continent. Taken together, these results suggest that ecology is an important axis of lineage divergence within this group and that more fine-scale analyses may find even greater differentiation between the populations identified here. This abstract translated to Spanish is avaliable in the Supporting Infromation section (Este resumen traducido al español está disponible en la sección, Supporting Infromation).  more » « less
Award ID(s):
2323125 1831241
PAR ID:
10568602
Author(s) / Creator(s):
; ; ; ; ;
Publisher / Repository:
Oxford University Press
Date Published:
Journal Name:
Zoological Journal of the Linnean Society
Volume:
203
Issue:
1
ISSN:
0024-4082
Format(s):
Medium: X
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 that 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. 
    more » « less
  2. 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 (but not preventing) gene flow and through parapatric differentiation along an ecological gradient. 
    more » « less
  3. Summary In the hyperdiverse fungi, the process of speciation is virtually unknown, including for the > 20 000 species of ectomycorrhizal mutualists. To understand this process, we investigated patterns of genome‐wide differentiation in the ectomycorrhizal porcini mushroom, Boletus edulis , a globally distributed species complex with broad ecological amplitude. By whole‐genome sequencing 160 individuals from across the Northern Hemisphere, we genotyped 792 923 single nucleotide polymorphisms to characterize patterns of genome‐wide differentiation and to identify the adaptive processes shaping global population structure. We show that B. edulis exhibits contrasting patterns of genomic divergence between continents, with multiple lineages present across North America, while a single lineage dominates Europe. These geographical lineages are inferred to have diverged 1.62–2.66 million years ago, during a period of climatic upheaval and the onset of glaciation in the Pliocene–Pleistocene boundary. High levels of genomic differentiation were observed among lineages despite evidence of substantial and ongoing introgression. Genome scans, demographic inference, and ecological niche models suggest that genomic differentiation is maintained by environmental adaptation, not physical isolation. Our study uncovers striking patterns of genome‐wide differentiation on a global scale and emphasizes the importance of local adaptation and ecologically mediated divergence, rather than prezygotic barriers such as allopatry or genomic incompatibility, in fungal population differentiation. 
    more » « less
  4. Abstract Divergence dating analyses in systematics provide a framework to develop and test biogeographic hypotheses regarding speciation. However, as molecular datasets grow from multilocus to genomic, sample sizes decrease due to computational burdens, and the testing of fine-scale biogeographic hypotheses becomes difficult. In this study, we use coalescent demographic models to investigate the diversification of poorly known rice paddy snakes from Southeast Asia (Homalopsidae:Hypsiscopus), which have conflicting dates of origin based on previous studies. We use coalescent modeling to test the hypothesis thatHypsiscopusdiversified 2.5 mya during the Khorat Plateau uplift in Thailand. Additionally, we use ecological niche analyses to identify potential differences in the niche space of the two most widely distributed species in the past and present. Our results suggestHypsiscopusdiversified ~ 2.4 mya, supporting that the Khorat Plateau may have initiated the diversification of rice paddy snakes. We also find significant niche differentiation and shifts between species ofHypsiscopus, indicating that environmental differences may have sustained differentiation of this genus after the Khorat Plateau uplift. Our study expands on the diversification history of snakes in Southeast Asia, and highlights how results from smaller multilocus datasets can be useful in developing and testing biogeographic hypotheses alongside genomic datasets. 
    more » « less
  5. Abstract Hybrid zones can be studied by modeling clines of trait variation (e.g., morphology, genetics) over a linear transect. Yet, hybrid zones can also be spatially complex, can shift over time, and can even lead to the formation of hybrid lineages with the right combination of dispersal and vicariance. We reassessed Sibley’s (1950) gradient between Collared Towhee (Pipilo ocai) and Spotted Towhee (Pipilo maculatus) in Central Mexico to test whether it conformed to a typical tension-zone cline model. By comparing historical and modern data, we found that cline centers for genetic and phenotypic traits have not shifted over the course of 70 years. This equilibrium suggests that secondary contact between these species, which originally diverged over 2 million years ago, likely dates to the Pleistocene. Given the amount of mtDNA divergence, parental ends of the cline have very low autosomal nuclear differentiation (FST = 0.12). Dramatic and coincident cline shifts in mtDNA and throat color suggest the possibility of sexual selection as a factor in differential introgression, while a contrasting cline shift in green back color hints at a role for natural selection. Supporting the idea of a continuum between clinal variation and hybrid lineage formation, the towhee gradient can be analyzed as one population under isolation-by-distance, as a two-population cline, and as three lineages experiencing divergence with gene flow. In the middle of the gradient, a hybrid lineage has become partly isolated, likely due to forested habitat shrinking and fragmenting as it moved upslope after the last glacial maximum and a stark environmental transition. This towhee system offers a window into the potential outcomes of hybridization across a dynamic landscape including the creation of novel genomic and phenotypic combinations and incipient hybrid lineages. 
    more » « less