Search

Significant advances have been made in species delimitation and numerous methods can test precisely defined models of speciation, though the synthesis of phylogeography and taxonomy is still sometimes incomplete. Emerging consensus treats distinct genealogical clusters in genome-scale data as strong initial evidence of speciation in most cases, a hypothesis that must therefore be falsified under an explicit evolutionary model. We can now test speciation hypotheses linking trait differentiation to specific mechanisms of divergence with increasingly large data sets. Integrative taxonomy can, therefore, reflect an understanding of how each axis of variation relates to underlying speciation processes, with nomenclature for distinct evolutionary lineages. We illustrate this approach here with Seal Salamanders (Desmognathus monticola) and introduce a new unsupervised machine-learning approach for species delimitation. Plethodontid salamanders are renowned for their morphological conservatism despite extensive phylogeographic divergence. We discover 2 geographic genetic clusters, for which demographic and spatial models of ecology and gene flow provide robust support for ecogeographic speciation despite limited phenotypic divergence. These data are integrated under evolutionary mechanisms (e.g., spatially localized gene flow with reduced migration) and reflected in emergent properties expected under models of reinforcement (e.g., ethological isolation and selection against hybrids). Their genetic divergence is prima facie evidence for species-level distinctiveness, supported by speciation models and divergence along axes such as behavior, geography, and climate that suggest an ecological basis with subsequent reinforcement through prezygotic isolation. As data sets grow more comprehensive, species-delimitation models can be tested, rejected, or corroborated as explicit speciation hypotheses, providing for reciprocal illumination of evolutionary processes and integrative taxonomies. [Desmognathus; integrative taxonomy; machine learning; species delimitation.]., (© The Author(s) 2022. Published by Oxford University Press, on behalf of the Society of Systematic Biologists. All rights reserved. For permissions, please email: journals.permissions@oup.com.)

A fundamental assumption of evolutionary biology is that phylogeny follows a bifurcating process. However, hybrid speciation and introgression are becoming more widely documented in many groups. Hybrid inference studies have been historically limited to small sets of taxa, while exploration of the prevalence and trends of reticulation at deep time scales remains unexplored. We study the evolutionary history of an adaptive radiation of 109 gemsnakes in Madagascar (Pseudoxyrhophiinae) to identify potential instances of introgression. Using several network inference methods, we find 12 reticulation events within the 22-million-year evolutionary history of gemsnakes, producing 28% of the diversity for the group, including one reticulation that resulted in the diversification of an 18 species radiation. These reticulations are found at nodes with high gene tree discordance and occurred among parental lineages distributed along a north-south axis that share similar ecologies. Younger hybrids occupy intermediate contact zones between the parent lineages showing that post-speciation dispersal in this group has not eroded the spatial signatures of introgression. Reticulations accumulated consistently over time, despite drops in overall speciation rates during the Pleistocene. This suggests that while bifurcating speciation rates may decline as the result of species accumulation and environmental change, speciation by hybridization may be more robust to these processes., (© The Author(s) 2023. Published by Oxford University Press on behalf of The Society for the Study of Evolution (SSE). All rights reserved. For permissions, please e-mail: journals.permissions@oup.com.)

Inferring the history of divergence between species in a framework that permits the presence of gene flow has been crucial for characterizing the "gray zone" of speciation, which is the period of time where lineages have diverged but have not yet achieved strict reproductive isolation. However, estimates of both divergence times and rates of gene flow often ignore spatial information, for example when considering the location and width of hybrid zones with respect to changes in the environment between lineages. Using population genomic data from the North American ratsnake complex (Pantherophis obsoletus), we connected phylogeographic estimates of lineage structure, migration, historical demography, and timing of divergence with hybrid zone dynamics. We examined the spatial context of diversification by linking migration and timing of divergence to the location and widths of hybrid zones. Artificial neural network approaches were applied to understand how landscape features and past climate have influenced population genetic structure among these lineages. We found that rates of migration between lineages were associated with the overall width of hybrid zones. Timing of divergence was not related to migration rate or hybrid zone width across species pairs but may be related to the number of alleles weakly introgressing through hybrid zones. This research underscores how incomplete reproductive isolation can be better understood by considering differential allelic introgression and the effects of historical and contemporary landscape features on the formation of lineages as well as overall genomic estimates of migration rates through time., (© 2020 The Authors. Evolution © 2020 The Society for the Study of Evolution.)

The southeastern Nearctic is a biodiversity hotspot that is also rich in cryptic species. Numerous hypotheses (e.g., vicariance, local adaptation, and Pleistocene speciation in glacial refugia) have been tested in an attempt to explain diversification and the observed pattern of extant biodiversity. However, previous phylogeographic studies have both supported and refuted these hypotheses. Therefore, while data support one or more of these diversification hypotheses, it is likely that taxa are forming within this region in species-specific ways. Here, we generate a genomic data set for the cornsnakes (Pantherophis guttatus complex), which are widespread across this region, spanning both biogeographic barriers and climatic gradients. We use phylogeographic model selection combined with hindcast ecological niche models to determine regions of habitat stability through time. This combined approach suggests that numerous drivers of population differentiation explain the current diversity of this group of snakes. The Mississippi River caused initial speciation in this species complex, with more recent divergence events linked to adaptations to ecological heterogeneity and allopatric Pleistocene refugia. Lastly, we discuss the taxonomy of this group and suggest there may be additional cryptic species in need of formal recognition., (© 2020 John Wiley & Sons Ltd.)

Genetic structure can be influenced by local adaptation to environmental heterogeneity and biogeographic barriers, resulting in discrete population clusters. Geographic distance among populations, however, can result in continuous clines of genetic divergence that appear as structured populations. Here, we evaluate the relevant importance of these three factors over a landscape characterized by environmental heterogeneity and the presence of a hypothesized biogeographic barrier in producing population genetic structure within 13 codistributed snake species using a genomic data set. We demonstrate that geographic distance and environmental heterogeneity across western North America contribute to population genomic divergence. Surprisingly, landscape features long thought to contribute to biogeographic barriers play little role in divergence community wide. Our results suggest that isolation by environment is the most important contributor to genomic divergence. Furthermore, we show that models of population clustering that incorporate spatial information consistently outperform nonspatial models, demonstrating the importance of considering geographic distances in population clustering. We argue that environmental and geographic distances as drivers of community-wide divergence should be explored before assuming the role of biogeographic barriers., (© 2019 John Wiley & Sons Ltd.)

Two main hypotheses have been proposed to explain the diversification of the Caatinga biota. The riverine barrier hypothesis (RBH) claims that the São Francisco River (SFR) is a major biogeographic barrier to gene flow. The Pleistocene climatic fluctuation hypothesis (PCH) states that gene flow, geographic genetic structure and demographic signatures on endemic Caatinga taxa were influenced by Quaternary climate fluctuation cycles. Herein, we analyse genetic diversity and structure, phylogeographic history, and diversification of a widespread Caatinga lizard (Cnemidophorus ocellifer) based on large geographical sampling for multiple loci to test the predictions derived from the RBH and PCH. We inferred two well-delimited lineages (Northeast and Southwest) that have diverged along the Cerrado-Caatinga border during the Mid-Late Miocene (6-14 Ma) despite the presence of gene flow. We reject both major hypotheses proposed to explain diversification in the Caatinga. Surprisingly, our results revealed a striking complex diversification pattern where the Northeast lineage originated as a founder effect from a few individuals located along the edge of the Southwest lineage that eventually expanded throughout the Caatinga. The Southwest lineage is more diverse, older and associated with the Cerrado-Caatinga boundaries. Finally, we suggest that C. ocellifer from the Caatinga is composed of two distinct species. Our data support speciation in the presence of gene flow and highlight the role of environmental gradients in the diversification process., (© 2015 John Wiley & Sons Ltd.)

Phylogenetic niche conservatism (PNC) typically refers to the tendency of closely related species to be more similar to each other in terms of niche than they are to more distant relatives. This has been implicated as a potential driving force in speciation and other species-richness patterns, such as latitudinal gradients. However, PNC has not been very well defined in most previous studies. Is it a pattern or a process? What are the underlying endogenous (e.g. genetic) and exogenous (e.g. ecological) factors that cause niches to be conserved? What degree of similarity is necessary to qualify as PNC? Is it possible for the evolutionary processes causing niches to be conserved to also result in niche divergence in different habitats? Here, we revisit these questions, codifying a theoretical and operational definition of PNC as a mechanistic evolutionary process resulting from several factors. We frame this both from a macroevolutionary and population-genetic perspective. We discuss how different axes of physical (e.g. geographic) and environmental (e.g. climatic) heterogeneity interact with the fundamental process of PNC to produce different outcomes of ecological speciation. We also review tests for PNC, and suggest ways that these could be improved or better utilized in future studies. Ultimately, PNC as a process has a well-defined mechanistic basis in organisms, and future studies investigating ecological speciation would be well served to consider this, and frame hypothesis testing in terms of the processes and expected patterns described herein. The process of PNC may lead to patterns where niches are conserved (more similar than expected), constrained (divergent within a limited subset of available niches), or divergent (less similar than expected), based on degree of phylogenetic relatedness between species., (© 2014 Cambridge Philosophical Society.)

Phylogenies are used to estimate rates of speciation and extinction, reconstruct historical diversification scenarios, and link these to ecological and evolutionary factors, such as climate or organismal traits. Recent models can now estimate the effects of binary, multistate, continuous, and biogeographic characters on diversification rates. Others test for diversity dependence (DD) in speciation and extinction, which has become recognized as an important process in numerous clades. A third class incorporates flexible time-dependent functions, enabling reconstruction of major periods of both expanding and contracting diversity. Although there are some potential problems (particularly for estimating extinction), these methods hold promise for answering many classic questions in ecology and evolution, such as the origin of adaptive radiations, and the latitudinal gradient in species richness., (Copyright © 2013 Elsevier Ltd. All rights reserved.)

Phylogeographic inference can determine the timing of population divergence, historical demographic processes, patterns of migration, and when extended to multiple species, the history of communities. Single-locus analyses can mislead interpretations of the evolutionary history of taxa and comparative analyses. It is therefore important to revisit previous single-locus phylogeographic studies, particularly those that have been used to propose general patterns for regional biotas and the processes responsible for generating inferred patterns. Here, we employ a multilocus statistical approach to re-examine the phylogeography of Lampropeltis zonata. Using nonparametic and Bayesian species delimitation, we determined that there are two well-supported species within L. zonata. Ecological niche modelling supports the delimitation of these taxa, suggesting that the two species inhabit distinct climatic environments. Gene flow between the two taxa is low and appears to occur unidirectionally. Further, our data suggest that gene flow was mediated by females, a rare pattern in snakes. In contrast to previous analyses, we determined that the divergence between the two lineages occurred in the late Pliocene (c. 2.07 Ma). Spatially and temporally, the divergence of these lineages is associated with the inundation of central California by the Monterey Bay. The effective population sizes of the two species appear to have been unaffected by Pleistocene glaciation. Our increased sampling of loci for L. zonata, combined with previously published multilocus analyses of other sympatric species, suggests that previous conclusions reached by comparative phylogeographic studies conducted within the California Floristic Province should be reassessed., (© 2013 John Wiley & Sons Ltd.)

For studies investigating trait evolution, there are at least two important questions. First, have traits under consideration influenced cladogenesis and extinction in the group? Second, how do fossil data alter inferences about trait evolution or diversification-rate dynamics? However, relatively few studies have assessed these questions. Here, we use recently developed methods to test for trait-dependent diversification in the New World colubrid snake tribe Lampropeltini. We also integrate data from fossil taxa into phylogenetic estimation of evolutionary parameters using a simple Monte Carlo randomization test. These analyses suggest that ecological conditions in temperate regions are tied to higher rates of cladogenesis, but that body size is not related to diversification in the group. We also find that the inclusion of fossil taxa alters absolute estimates of size and the rate of size evolution, but not the overall pattern of ecomorphological diversification, as well as estimates of evolutionary rates, particularly extinction., (© 2012 The Authors. Journal of Evolutionary Biology © 2012 European Society For Evolutionary Biology.)

Studies of speciation and taxon delimitation are usually decoupled. Combining these methods provides a stronger theoretical ground for recognizing new taxa and understanding processes of speciation. Using coalescent methods, we examine speciation, post-speciation population demographics, and taxon delimitation in the Arizona Mountain Kingsnake (Lampropeltis pyromelana), a species restricted to high elevations in southwestern United States and northern Mexico (SW). These methods provide a solid foundation for understanding how biogeographic barriers operate at the regional scale in the SW. Bayesian species delimitation methods, using three loci from samples of L. pyromelana taken throughout their range, show strong support for the existence of two species that are separated by low elevation habitats found between the Colorado Plateau/ Mogollon Rim and the Sierra Madre Occidental. Our results suggest an allopatric mode of speciation given the near absence of gene flow over time, which resulted in two lineages of unequal population sizes. Speciation likely occurred prior to the Pleistocene, during the aridification of the SW and/or the uplift of the Colorado Plateau, and while these species occupy similar high-elevation niches, they are isolated by xeric conditions found in the intervening low deserts. Furthermore, post-speciation demographics suggest that populations of both lineages were not negatively impacted by climate change throughout the Pleistocene. Finally, our results suggest that at least for this group, where divergence is old and gene flow is low, Bayesian species delimitation performs well., (Copyright © 2011 Elsevier Inc. All rights reserved.)

Molecular phylogenies are often used to test hypotheses about the tempo and mode of speciation and extinction. One commonly used statistic is Pybus and Harvey's γ, which measures the density of ordered internode distances on an ultrametric tree to infer earlier (negative γ) or later (positive γ) bursts of diversification. However, coalescent theory predicts that γ might be biased toward negative values (inferring early bursts of diversification) when using gene trees rather than species trees. Gene divergences predate species divergences, increasingly so at higher effective population sizes (N(e)), and proportionally more so toward the tips of the tree. Thus, gene trees will have a higher density of older nodes in many cases (particularly at higher N(e)), due to the disproportionate lengthening of terminal branches. This will yield an artifactual signature of early bursts of diversification when estimating γ from gene trees. We simulate gene trees within species trees under both Yule (pure-birth) and birth-death processes, and demonstrate support for these predictions. However, for most realistic estimates of θ in natural populations, gene trees provide relatively good estimates of γ, despite the disproportionate overestimation of younger node ages. This is corroborated with an empirical dataset of North American fence lizards (Sceloporus)., (© 2011 The Author(s). Evolution© 2011 The Society for the Study of Evolution.)

Ecological adaptive radiation theory predicts an increase in both morphological and specific diversification when organisms colonize new environments. Accordingly, bursts of morphological diversification, characterized by low within-subclade morphological disparity, may be associated with these increases in speciation rates. Conversely, increasing species density, reduction in available habitat, or increasing extinction rates are expected to cause rates of diversification to decline. We test these hypotheses by examining the tempo and mode of speciation in the lampropeltinine snakes, a morphologically variable group that colonized the New World approximately 24 million years ago and radiated throughout the Miocene. We show that specific diversification increased early in the history of the group, and that most morphological variation is partitioned among, rather than within subclades. These patterns provide further evidence for the hypothesis that morphological variation tends to be strongly partitioned among lineages when clades undergo early bursts of species diversification. A reduction in speciation rates may be indicative of density dependent effects due to a saturation of available ecological opportunity, rather than increases in extinction rates at the onset of the Pleistocene/Pliocene glacial cycles. This evidence runs counter to the general Pleistocene species pump model.

Niche conservatism and niche divergence are both important ecological mechanisms associated with promoting allopatric speciation across geographical barriers. However, the potential for variable responses in widely distributed organisms has not been fully investigated. For allopatric sister lineages, three patterns for the interaction of ecological niche preference and geographical barriers are possible: (i) niche conservatism at a physical barrier; (ii) niche divergence at a physical barrier; and (iii) niche divergence in the absence of a physical barrier. We test for the presence of these patterns in a transcontinentally distributed snake species, the common kingsnake (Lampropeltis getula), to determine the relative frequency of niche conservatism or divergence in a single species complex inhabiting multiple distinct ecoregions. We infer the phylogeographic structure of the kingsnake using a range-wide data set sampled for the mitochondrial gene cytochrome b. We use coalescent simulation methods to test for the presence of structured lineage formation vs. fragmentation of a widespread ancestor. Finally, we use statistical techniques for creating and evaluating ecological niche models to test for conservatism of ecological niche preferences. Significant geographical structure is present in the kingsnake, for which coalescent tests indicate structured population division. Surprisingly, we find evidence for all three patterns of conservatism and divergence. This suggests that ecological niche preferences may be labile on recent phylogenetic timescales, and that lineage formation in widespread species can result from an interaction between inertial tendencies of niche conservatism and natural selection on populations in ecologically divergent habitats.

Sky islands provide ideal opportunities for understanding how climatic changes associated with Pleistocene glacial cycles influenced species distributions, genetic diversification, and demography. The salamander Plethodon ouachitae is largely restricted to high-elevation, mesic forest on six major mountains in the Ouachita Mountains. Because these mountains are separated by more xeric, low-elevation valleys, the salamanders appear to be isolated on sky islands where gene flow among populations on different mountains may be restricted. We used DNA sequence data along with ecological niche modelling and coalescent simulations to test several hypotheses related to diversifications in sky island habitats. Our results revealed that P. ouachitae is composed of seven well-supported lineages structured across six major mountains. The species originated during the Late Pliocene, and lineage diversification occurred during the Middle Pleistocene in a stepping stone fashion with a cyclical pattern of dispersal to a new mountain followed by isolation and divergence. Diversification occurred primarily on an east-west axis, which is likely related to the east-west orientation of the Ouachita Mountains and the more favourable cooler and wetter environmental conditions on north slopes compared to south-facing slopes and valleys. All non-genealogical coalescent methods failed to detect significant population expansion in any lineages. Bayesian skyline plots showed relatively stable population sizes over time, but indicated a slight to moderate amount of population growth in all lineages starting approximately 10 000-12 000 years ago. Our results provide new insight into sky island diversifications from a previously unstudied region, and further demonstrate that climatic changes during the Pleistocene had profound effects on lineage diversification and demography, especially in species from environmentally sensitive habitats in montane regions.