Search

Genomics is narrowing uncertainty in the phylogenetic structure for many amniote groups. For one of the most diverse and species-rich groups, the squamate reptiles (lizards, snakes, and amphisbaenians), an inverse correlation between the number of taxa and loci sampled still persists across all publications using DNA sequence data and reaching a consensus on the relationships among them has been highly problematic. In this study, we use high-throughput sequence data from 289 samples covering 75 families of squamates to address phylogenetic affinities, estimate divergence times, and characterize residual topological uncertainty in the presence of genome-scale data. Importantly, we address genomic support for the traditional taxonomic groupings Scleroglossa and Macrostomata using novel machine-learning techniques. We interrogate genes using various metrics inherent to these loci, including parsimony-informative sites (PIS), phylogenetic informativeness, length, gaps, number of substitutions, and site concordance to understand why certain loci fail to find previously well-supported molecular clades and how they fail to support species-tree estimates. We show that both incomplete lineage sorting and poor gene-tree estimation (due to a few undesirable gene properties, such as an insufficient number of PIS), may account for most gene and species-tree discordance. We find overwhelming signal for Toxicofera, and also show that none of the loci included in this study supports Scleroglossa or Macrostomata. We comment on the origins and diversification of Squamata throughout the Mesozoic and underscore remaining uncertainties that persist in both deeper parts of the tree (e.g., relationships between Dibamia, Gekkota, and remaining squamates; among the three toxicoferan clades Iguania, Serpentes, and Anguiformes) and within specific clades (e.g., affinities among gekkotan, pleurodont iguanians, and colubroid families)., (© The Author(s) 2019. Published by Oxford University Press, on behalf of the Society of Systematic Biologists. All rights reserved. For permissions, please email: journals.permissions@oup.com.)

Processes leading to spectacular diversity of both form and species on islands have been well-documented under island biogeography theory, where distance from source and island size are key factors determining immigration and extinction resistance. But far less understood are the processes governing in situ diversification on the world's mega islands, where large and isolated land masses produced morphologically distinct radiations from related taxa on continental regions. Madagascar has long been recognized as a natural laboratory due to its isolation, lack of influence from adjacent continents, and diversification of spectacular vertebrate radiations. However, only a handful of studies have examined rate shifts of in situ diversification for this island. Here, we examine rates of diversification in the Malagasy snakes of the family Pseudoxyrhophiinae (gemsnakes) to understand if rates of speciation were initially high, enhanced by diversification into distinct biomes, and associated with key dentition traits. Using a genomic sequence-capture data set for 366 samples, we determine that all previously described and newly discovered species are delimitable and therefore useful candidates for understanding diversification trajectories through time. Our analysis detected no shifts in diversification rate between clades or changes in biome or dentition type. Remarkably, we demonstrate that rates of diversification of the gemsnake radiation, which originated in Madagascar during the early Miocene, remained steady throughout the Neogene. However, we do detect a significant slowdown in diversification during the Pleistocene. We also comment on the apparent paradox where most living species originated in the Pleistocene, despite diversification rates being substantially higher during the earlier 15 myr., (© The Author(s) 2019. Published by Oxford University Press, on behalf of the Society of Systematic Biologists. All rights reserved. For permissions, please email: journals.permissions@oup.com.)

Most phylogenies are typically represented as purely bifurcating. However, as genomic data have become more common in phylogenetic studies, it is not unusual to find reticulation among terminal lineages or among internal nodes (deep time reticulation; DTR). In these situations, gene flow must have happened in the same or adjacent geographic areas for these DTRs to have occurred and therefore biogeographic reconstruction should provide similar area estimates for parental nodes, provided extinction or dispersal has not eroded these patterns. We examine the phylogeny of the widely distributed New World kingsnakes (Lampropeltis), determine if DTR is present in this group, and estimate the ancestral area for reticulation. Importantly, we develop a new method that uses coalescent simulations in a machine learning framework to show conclusively that this phylogeny is best represented as reticulating at deeper time. Using joint probabilities of ancestral area reconstructions on the bifurcating parental lineages from the reticulating node, we show that this reticulation likely occurred in northwestern Mexico/southwestern US, and subsequently, led to the diversification of the Mexican kingsnakes. This region has been previously identified as an area important for understanding speciation and secondary contact with gene flow in snakes and other squamates. This research shows that phylogenetic reticulation is common, even in well-studied groups, and that the geographic scope of ancient hybridization is recoverable.

Emerging infectious diseases (EIDs) reduce host population sizes, cause extinction, disassemble communities, and have indirect negative effects on human well-being. Fungal EIDs have reduced population abundances in amphibians and bats across many species over large areas. The recent emergence of snake fungal disease (SFD) may have caused declines in some snake populations in the Eastern United States (EUS), which is home to a phylogenetically and ecologically diverse assembly of 98 taxa. SFD has been documented in only 23 naturally occuring species, although this is likely an underestimate of the number of susceptible taxa. Using several novel methods, including artificial neural networks, we combine phylogenetic and trait-based community estimates from all taxa in this region to show that SFD hosts are both phylogenetically and ecologically randomly dispersed. This might indicate that other species of snakes in the EUS could be currently infected or susceptible to SFD. Our models also indicate that information about key traits that enhance susceptiblity is lacking. Surveillance should consider that all snake species and habitats likely harbor this pathogen.

Globally distributed groups may show regionally distinct rates of diversification, where speciation is elevated given timing and sources of ecological opportunity. However, for most organisms, nearly complete sampling at genomic-data scales to reduce topological error in all regions is unattainable, thus hampering conclusions related to biogeographic origins and rates of diversification. We explore processes leading to the diversity of global ratsnakes and test several important hypotheses related to areas of origin and enhanced diversification upon colonizing new continents. We estimate species trees inferred from phylogenomic scale data (304 loci) while exploring several strategies that consider topological error from each individual gene tree. With a dated species tree, we examine taxonomy and test previous hypotheses that suggest the ratsnakes originated in the Old World (OW) and dispersed to New World (NW). Furthermore, we determine if dispersal to the NW represented a source of ecological opportunity, which should show elevated rates of species diversification. We show that ratsnakes originated in the OW during the mid-Oligocene and subsequently dispersed to the NW by the mid-Miocene; diversification was also elevated in a subclade of NW taxa. Finally, the optimal biogeographic region-dependent speciation model shows that the uptick in ratsnake diversification was associated with colonization of the NW. We consider several alternative explanations that account for regionally distinct diversification rates., (Copyright © 2017 Elsevier Inc. All rights reserved.)

Predicting species presence and richness on islands is important for understanding the origins of communities and how likely it is that species will disperse and resist extinction. The equilibrium theory of island biogeography (ETIB) and, as a simple model of sampling abundances, the unified neutral theory of biodiversity (UNTB), predict that in situations where mainland to island migration is high, species-abundance relationships explain the presence of taxa on islands. Thus, more abundant mainland species should have a higher probability of occurring on adjacent islands. In contrast to UNTB, if certain groups have traits that permit them to disperse to islands better than other taxa, then phylogeny may be more predictive of which taxa will occur on islands. Taking surveys of 54 island snake communities in the Eastern Nearctic along with mainland communities that have abundance data for each species, we use phylogenetic assembly methods and UNTB estimates to predict island communities. Species richness is predicted by island area, whereas turnover from the mainland to island communities is random with respect to phylogeny. Community structure appears to be ecologically neutral and abundance on the mainland is the best predictor of presence on islands. With regard to young and proximate islands, where allopatric or cladogenetic speciation is not a factor, we find that simple neutral models following UNTB and ETIB predict the structure of island communities., (© 2015 The Author(s).)

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.)

Next-generation genomic sequencing promises to quickly and cheaply resolve remaining contentious nodes in the Tree of Life, and facilitates species-tree estimation while taking into account stochastic genealogical discordance among loci. Recent methods for estimating species trees bypass full likelihood-based estimates of the multi-species coalescent, and approximate the true species-tree using simpler summary metrics. These methods converge on the true species-tree with sufficient genomic sampling, even in the anomaly zone. However, no studies have yet evaluated their efficacy on a large-scale phylogenomic dataset, and compared them to previous concatenation strategies. Here, we generate such a dataset for Caenophidian snakes, a group with >2500 species that contains several rapid radiations that were poorly resolved with fewer loci. We generate sequence data for 333 single-copy nuclear loci with ∼100% coverage (∼0% missing data) for 31 major lineages. We estimate phylogenies using neighbor joining, maximum parsimony, maximum likelihood, and three summary species-tree approaches (NJst, STAR, and MP-EST). All methods yield similar resolution and support for most nodes. However, not all methods support monophyly of Caenophidia, with Acrochordidae placed as the sister taxon to Pythonidae in some analyses. Thus, phylogenomic species-tree estimation may occasionally disagree with well-supported relationships from concatenated analyses of small numbers of nuclear or mitochondrial genes, a consideration for future studies. In contrast for at least two diverse, rapid radiations (Lamprophiidae and Colubridae), phylogenomic data and species-tree inference do little to improve resolution and support. Thus, certain nodes may lack strong signal, and larger datasets and more sophisticated analyses may still fail to resolve them., (Copyright © 2014 Elsevier Inc. All rights reserved.)

Africa is renowned for its biodiversity and endemicity, yet little is known about the factors shaping them across the continent. African Agama lizards (45 species) have a pan-continental distribution, making them an ideal model for investigating biogeography. Many species have evolved conspicuous sexually dimorphic traits, including extravagant breeding coloration in adult males, large adult male body sizes, and variability in social systems among colorful versus drab species. We present a comprehensive time-calibrated species tree for Agama, and their close relatives, using a hybrid phylogenetic-phylogenomic approach that combines traditional Sanger sequence data from five loci for 57 species (146 samples) with anchored phylogenomic data from 215 nuclear genes for 23 species. The Sanger data are analyzed using coalescent-based species tree inference using (*)BEAST, and the resulting posterior distribution of species trees is attenuated using the phylogenomic tree as a backbone constraint. The result is a time-calibrated species tree for Agama that includes 95% of all species, multiple samples for most species, strong support for the major clades, and strong support for most of the initial divergence events. Diversification within Agama began approximately 23 million years ago (Ma), and separate radiations in Southern, East, West, and Northern Africa have been diversifying for >10Myr. A suite of traits (morphological, coloration, and sociality) are tightly correlated and show a strong signal of high morphological disparity within clades, whereby the subsequent evolution of convergent phenotypes has accompanied diversification into new biogeographic areas., (Copyright © 2014 Elsevier Inc. All rights reserved.)

Both gene-tree discordance and unrecognized diversity are sources of error for accurate estimation of species trees, and can affect downstream diversification analyses by obscuring the correct number of nodes, their density, and the lengths of the branches subtending them. Although the theoretical impact of gene-tree discordance on evolutionary analyses has been examined previously, the effect of unsampled and cryptic diversity has not. Here, we examine how delimitation of previously unrecognized diversity in the milksnake (Lampropeltis triangulum) and use of a species-tree approach affects both estimation of the Lampropeltis phylogeny and comparative analyses with respect to the timing of diversification. Coalescent species delimitation indicates that L. triangulum is not monophyletic and that there are multiple species of milksnake, which increases the known species diversity in the genus Lampropeltis by 40%. Both genealogical and temporal discordance occurs between gene trees and the species tree, with evidence that mitochondrial DNA (mtDNA) introgression is a main factor. This discordance is further manifested in the preferred models of diversification, where the concatenated gene tree strongly supports an early burst of speciation during the Miocene, in contrast to species-tree estimates where diversification follows a birth-death model and speciation occurs mostly in the Pliocene and Pleistocene. This study highlights the crucial interaction among coalescent-based phylogeography and species delimitation, systematics, and species diversification analyses.

Delimiting species and clarifying phylogenetic relationships are the main goals of systematics. For species with questionable taxonomic status, species delimitation approaches using multi-species coalescent models with multiple loci are recommended if morphological data are unavailable or unhelpful. Moreover, these methods will also reduce subjectivity based on genetic distance or requirement of monophyletic genetic lineages. We determine the validity and phylogenetic position of a rare and long controversial species of Chinese reptile, the Szechwan ratsnake (Euprepiophis perlaceus), using multi-locus data from multiple individuals and coalescent-based approaches. Species were first delimited using Bayesian Phylogenetics & Phylogeography (BP&P), Brownie and Bayes Factor model comparison approaches, while relationships among species were estimated using species tree inference in (*)BEAST. Results indicate that Euprepiophis perlaceus is a distinct species sister to Euprepiophis mandarinus. Despite gene tree discrepancy, the coalescent model-based approaches used here demonstrate the taxonomic validity and the phylogenetic position of Euprepiophis perlaceus. These approaches objectively test the validity of questionable species diagnoses based on morphological characters and determine their phylogenetic position., (Copyright © 2013 Elsevier Inc. All rights reserved.)

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.)

Spiny-rayed fishes, or acanthomorphs, comprise nearly one-third of all living vertebrates. Despite their dominant role in aquatic ecosystems, the evolutionary history and tempo of acanthomorph diversification is poorly understood. We investigate the pattern of lineage diversification in acanthomorphs by using a well-resolved time-calibrated phylogeny inferred from a nuclear gene supermatrix that includes 520 acanthomorph species and 37 fossil age constraints. This phylogeny provides resolution for what has been classically referred to as the "bush at the top" of the teleost tree, and indicates acanthomorphs originated in the Early Cretaceous. Paleontological evidence suggests acanthomorphs exhibit a pulse of morphological diversification following the end Cretaceous mass extinction; however, the role of this event on the accumulation of living acanthomorph diversity remains unclear. Lineage diversification rates through time exhibit no shifts associated with the end Cretaceous mass extinction, but there is a global decrease in lineage diversification rates 50 Ma that occurs during a period when morphological disparity among fossil acanthomorphs increases sharply. Analysis of clade-specific shifts in diversification rates reveal that the hyperdiversity of living acanthomorphs is highlighted by several rapidly radiating lineages including tunas, gobies, blennies, snailfishes, and Afro-American cichlids. These lineages with high diversification rates are not associated with a single habitat type, such as coral reefs, indicating there is no single explanation for the success of acanthomorphs, as exceptional bouts of diversification have occurred across a wide array of marine and freshwater habitats.

Background: The extant squamates (>9400 known species of lizards and snakes) are one of the most diverse and conspicuous radiations of terrestrial vertebrates, but no studies have attempted to reconstruct a phylogeny for the group with large-scale taxon sampling. Such an estimate is invaluable for comparative evolutionary studies, and to address their classification. Here, we present the first large-scale phylogenetic estimate for Squamata., Results: The estimated phylogeny contains 4161 species, representing all currently recognized families and subfamilies. The analysis is based on up to 12896 base pairs of sequence data per species (average = 2497 bp) from 12 genes, including seven nuclear loci (BDNF, c-mos, NT3, PDC, R35, RAG-1, and RAG-2), and five mitochondrial genes (12S, 16S, cytochrome b, ND2, and ND4). The tree provides important confirmation for recent estimates of higher-level squamate phylogeny based on molecular data (but with more limited taxon sampling), estimates that are very different from previous morphology-based hypotheses. The tree also includes many relationships that differ from previous molecular estimates and many that differ from traditional taxonomy., Conclusions: We present a new large-scale phylogeny of squamate reptiles that should be a valuable resource for future comparative studies. We also present a revised classification of squamates at the family and subfamily level to bring the taxonomy more in line with the new phylogenetic hypothesis. This classification includes new, resurrected, and modified subfamilies within gymnophthalmid and scincid lizards, and boid, colubrid, and lamprophiid snakes.

Inferring the evolutionary and biogeographic history of taxa occurring in a particular region is one way to determine the processes by which the biodiversity of that region originated. Tree boas of the genus Corallus are an ancient clade and occur throughout Central and South America and the Lesser Antilles, making it an excellent group for investigating Neotropical biogeography. Using sequenced portions of two mitochondrial and three nuclear loci for individuals of all recognized species of Corallus, we infer phylogenetic relationships, present the first molecular analysis of the phylogenetic placement of the enigmatic C. cropanii, develop a time-calibrated phylogeny, and explore the biogeographic history of the genus. We found that Corallus diversified within mainland South America, via over-water dispersals to the Lesser Antilles and Central America, and via the traditionally recognized Panamanian land bridge. Divergence time estimates reject the South American Caribbean-Track as a general biogeographic model for Corallus and implicate a role for events during the Oligocene and Miocene in diversification such as marine incursions and the uplift of the Andes. Our findings also suggest that recognition of the island endemic species, C. grenadensis and C. cookii, is questionable as they are nested within the widely distributed species, C. hortulanus. Our results highlight the importance of using widespread taxa when forming and testing biogeographic hypotheses in complex regions and further illustrate the difficulty of forming broadly applicable hypotheses regarding patterns of diversification in the Neotropical region., (Copyright © 2012 Elsevier Inc. All rights reserved.)

Snake diversity in the island of Sri Lanka is extremely high, hosting at least 89 inland (i.e., non-marine) snake species, of which at least 49 are endemic. This includes the endemic genera Aspidura, Balanophis, Cercaspis, Haplocercus, and Pseudotyphlops, which are of uncertain phylogenetic affinity. We present phylogenetic evidence from nuclear and mitochondrial loci showing the relationships of 40 snake species from Sri Lanka (22 endemics) to the remaining global snake fauna. To determine the phylogenetic placement of these species, we create a molecular dataset containing 10 genes for all global snake genera, while also sampling all available species for genera with endemic species occurring in Sri Lanka. Our sampling comprises five mitochondrial genes (12S, 16S, cyt-b, ND2, and ND4) and five nuclear genes (BDNF, c-mos, NT3 RAG-1, and RAG-2), for a total of up to 9582bp per taxon. We find that the five endemic genera represent portions of four independent colonizations of Sri Lanka, with Cercaspis nested within Colubrinae, Balanophis in Natricinae, Pseudotyphlops in Uropeltidae, and that Aspidura+Haplocercus represents a distinct, ancient lineage within Natricinae. We synonymize two endemic genera that render other genera paraphyletic (Haplocercus with Aspidura, and Cercaspis with Lycodon), and discover that further endemic radiations may be present on the island, including a new taxon from the blindsnake family Typhlopidae, suggesting a large endemic radiation. Despite its small size relative to other islands such as New Guinea, Borneo, and Madagascar, Sri Lanka has one of the most phylogenetically diverse island snake faunas in the world, and more research is needed to characterize the island's biodiversity, with numerous undescribed species in multiple lineages., (Copyright © 2012 Elsevier Inc. All rights reserved.)

Numerous taxa show ancient intertropical disjunct distributions. Many can be explained by well-known processes of historical vicariance, such as the breakup of Gondwanaland. Others, such as Asian-Neotropical divergences are not as well understood. To clarify the phylogenetic position and understand biogeographic and temporal origins of the geographically disjunct and morphologically unique genera of hinged-teeth snakes, Scaphiodontophis (n=1) and Sibynophis (n=9; Colubridae), we inferred a time-calibrated phylogeny with additional 107 taxa representing the superfamily Colubroidea using four genes (c-mos, cyt-b, ND2, RAG-1; 3085 bp). We used this tree to estimate ancestral areas for the group. The results show that Scaphiodontophis is sister to Sibynophis, both originated in the late Eocene/Oligocene in Asia and likely dispersed through Beringia to the New World, but unlike other snake groups left no extant species in temperate North America. Current recognition of Scaphiodontophiinae renders Colubrinae paraphyletic, and we resurrect the previously named subfamily Sibynophiinae to encompass both genera and use the tribes Sibynophiini (Sibynophis) and Scaphiodontophiini (Scaphiodontophis) to highlight the geographically distinct areas occupied by these taxa. These results suggest that intercontinental dispersal with extinction in intermediate areas can explain puzzling patterns of ancient intertropical disjunct distributions., (Copyright © 2012 Elsevier Inc. All rights reserved.)

Adaptive radiation (AR) theory predicts that groups sharing the same source of ecological opportunity (EO) will experience deterministic species diversification and morphological evolution. Thus, deterministic ecological and morphological evolution should be correlated with deterministic patterns in the tempo and mode of speciation for groups in similar habitats and time periods. We test this hypothesis using well-sampled phylogenies of four squamate groups that colonized the New World (NW) in the Late Oligocene. We use both standard and coalescent models to assess species diversification, as well as likelihood models to examine morphological evolution. All squamate groups show similar early pulses of speciation, as well as diversity-dependent ecological limits on clade size at a continental scale. In contrast, processes of morphological evolution are not easily predictable and do not show similar pulses of early and rapid change. Patterns of morphological and species diversification thus appear uncoupled across these groups. This indicates that the processes that drive diversification and disparification are not mechanistically linked, even among similar groups of taxa experiencing the same sources of EO. It also suggests that processes of phenotypic diversification cannot be predicted solely from the existence of an AR or knowledge of the process of diversification.

The statistical rigor of species delimitation has increased dramatically over the past decade. Coalescent theory provides powerful models for population genetic inference, and is now increasingly important in phylogenetics and speciation research. By applying probabilistic models, coalescent-based species delimitation provides clear and objective testing of alternative hypotheses of evolutionary independence. As acquisition of multilocus data becomes increasingly automated, coalescent-based species delimitation will improve the discovery, resolution, consistency, and stability of the taxonomy of species. Along with other tools and data types, coalescent-based species delimitation will play an important role in an integrative taxonomy that emphasizes the identification of species limits and the processes that have promoted lineage diversification., (Copyright © 2012 Elsevier Ltd. All rights reserved.)

Spatial and temporal heterogeneity in environmental factors can have profound effects on diversification in species that are tightly linked to their environments. The Caddo Mountain Salamander (Plethodon caddoensis) inhabits a unique physiographic section of the Ouachita Mountains in central North America, a region in which Pleistocene climatic fluctuations have been implicated in driving lineage diversification in two other closely related salamanders. We examined P. caddoensis to determine whether it was similarly impacted by historic climatic changes and test whether physiographic features unique to the area also contributed to its diversification. We found that P. caddoensis is composed of four highly divergent, geographically distinct lineages that abut one another along an east-west axis. Phylogeographic structure was significantly related to both geographic distance and stream drainages, indicating that connectivity of streams and stream-associated habitats (e.g., talus) influence patterns of interpopulation gene flow. Lineages originated during the Middle Miocene and population size decreased in all lineages during the Pleistocene. Surface Geology and precipitation were the most important variables predicting the species distribution. Our results show that the unique physiographic features of the area coupled with species response to climatic factors have driven lineage diversification and phylogeographic structure in P. caddoensis. Variation in responses to historic climatic fluctuations among salamander species in this region underscore the importance of integrating species ecology with other factors such as geology and hydrology in order to better understand the effects of climate change on species with close associations to their environments., (Copyright © 2011 Elsevier Inc. All rights reserved.)

The superfamily Colubroidea (> 2500 species) includes the majority of snake species and is one of the most conspicuous and well-known radiations of terrestrial vertebrates. However, many aspects of the phylogeny of the group remain contentious, and dozens of genera have yet to be included in molecular phylogenetic analyses. We present a new, large-scale, likelihood-based phylogeny for the colubroids, including 761 species sampled for up to five genes: cytochrome b (93% of 761 species sampled), ND4 (69%), ND2 (28%), c-mos (54%), and RAG-1 (13%), totaling up to 5814bp per species. We also compare likelihood bootstrapping and a recently proposed ultra-fast measure of branch support (Shimodaira-Hasegawa-like [SHL] approximate likelihood ratio), and find that the SHL test shows strong support for several clades that were weakly-supported by bootstrapping in this or previous analyses (e.g., Dipsadinae, Lamprophiidae). We find that SHL values are positively related to branch lengths, but show stronger support for shorter branches than bootstrapping. Despite extensive missing data for many taxa (mean=67% per species), neither bootstrap nor SHL support values for terminal species are related to their incompleteness, and that most highly incomplete taxa are placed in the expected families from previous taxonomy, typically with very strong support. The phylogeny indicates that the Neotropical colubrine genus Scaphiodontophis represents an unexpectedly ancient lineage within Colubridae. We present a revised higher-level classification of Colubroidea, which includes a new subfamily for Scaphiodontophis (Scaphiodontophiinae). Our study provides the most comprehensive phylogeny of Colubroidea to date, and suggests that SHL values may provide a useful complement to bootstrapping for estimating support on likelihood-based trees., (Copyright © 2010 Elsevier Inc. All rights reserved.)

Climatic changes associated with Pleistocene glacial cycles profoundly affected species distributions, patterns of interpopulation gene flow, and demography. In species restricted to montane habitats, ranges may expand and contract along an elevational gradients in response to environmental fluctuations and create high levels of genetic variation among populations on different mountains. The salamander Plethodon fourchensis is restricted to high-elevation, mesic forest on five montane isolates in the Ouachita Mountains. We used DNA sequence data along with ecological niche modelling and coalescent simulations to test several hypotheses related to the effects of Pleistocene climatic fluctuations on species in montane habitats. Our results revealed that P. fourchensis is composed of four well-supported, geographically structured lineages. Geographic breaks between lineages occurred in the vicinity of major valleys and a narrow high-elevation pass. Ecological niche modelling predicted that environmental conditions in valleys separating most mountains are suitable; however, interglacial periods like the present are predicted to be times of range expansion in P. fourchensis. Divergence dating and coalescent simulations indicated that lineage diversification occurred during the Middle Pleistocene via the fragmentation of a wide-ranging ancestor. Bayesian skyline plots showed gradual decreases in population size in three of four lineages over the most recent glacial period and a slight to moderate amount of population growth during the Holocene. Our results not only demonstrate that climatic changes during the Pleistocene had profound effects on species restricted to montane habitats, but comparison of our results for P. fourchensis with its parapatric, sister taxon, P. ouachitae, also emphasizes how responses can vary substantially even among closely related, similarly distributed taxa.

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.

Many studies have revealed that lineages currently inhabiting formerly glaciated areas were pushed into southern glacial refugia and have expanded into their modern range since the last glacial maximum. There have been few studies that compare the effects of glacial cycles on lineage diversification and historical demography in closely related species with overlapping ranges. In this study we compare phylogeographic structure, historical demography, and approximate lineage age in two closely related and broadly co-occurring venomous snakes in eastern North America, the cottonmouth (Agkistrodon piscivorus) and copperhead (A. contortrix) using sequences from the mtDNA gene cytochrome b. We inferred three geographic lineages of A. contortrix and two of A. piscivorus with no common geographic or temporal pattern of lineage diversification identified for these species. Lineage diversification occurred in the Late Pliocene for A. piscivorus (approximately 2.5mya) and in the Early Pleistocene for A. contortrix ( approximately 1.5mya). Demographic estimates revealed population expansion following the last glacial maximum (approximately 20,000 years ago) in two lineages of A. contortrix (the Central clade and Eastern clade) and one lineage of A. piscivorus (the Continental clade). The Florida clade of A. piscivorus is the only lineage for which constant population size through time was inferred, possibly due to stable populations persisting in areas unaffected by glacial advances. Our data suggest that unique habitat preferences may have shaped both the phylogeographic and demographic histories of each species.

Most phylogeographic studies examine organisms that do not have transcontinental distributions and therefore the genetic and temporal effects of barriers across an entire continent cannot be assessed with respect to a single species. We examined the phylogeographic structure, lineage age, and historical demography using sequences from the mtDNA cytochrome b gene of the widespread North American racer (Coluber constrictor), one of the few abundant transcontinental snakes that occurs throughout many diverse biomes. Our results indicate that this complex is comprised of six lineages differing greatly in geographic extent, with the largest (a central US clade) being approximately 26 times greater than the smallest (a lineage restricted to the Florida Panhandle and nearby portions of adjacent States). Most of the six lineages appear to be separated at previously identified genetic barriers for several vertebrates with similar ranges. Lineage diversification in this species began in the late Miocene, separating populations in the Florida Peninsula from the remainder of the US. Diversification of lineages continued throughout the Pliocene and early Pleistocene. Four of the six lineages occur east of the Mississippi River, with only two distinctly young ( approximately 1.5 mya) lineages found west of the Mississippi River (one occurs west of Continental Divide). All methods of demographic inference, including the mismatch distribution, Fu and Li's D and Tajima's D , and Bayesian skyline plots revealed population expansion occurring in the mid-to-late Pleistocene for every lineage, regardless of size or proximity to formerly glaciated areas. Population expansion for lineages found east of the Mississippi River occurred earlier and was much greater than those found west of the River.

Dynamic climatic oscillations during the Pleistocene had profound effects on the distributions of species across North America. Although the role of historical climate change on speciation remains controversial, the impact on genetic variation within species has been well documented. We examined mtDNA sequences from the cytochrome b gene (1117 bp) and a portion of the NADH-4 gene (659 bp) for 286 individuals of Diadophis punctatus to infer phylogeographic patterns and population structure and to examine historical demographic patterns in both glaciated and unglaciated regions of North America. We inferred 14 lineages that replace each other geographically across the United States. Several of these lineages appear to be confined to specific habitats (floodplains, grasslands, montane environments) and traverse previously identified genetic barriers for terrestrial vertebrates including the Mississippi and Apalachicola Rivers, the Appalachian Mountains, and the western continental divide. We also observed overlapping ranges between some haplotype groups and several instances of secondary contact associated with ecological transition zones in eastern South Carolina, southern Oklahoma and central California. Within the US, diversification began during the late Miocene and continued into the mid-Pleistocene, suggesting these lineages pre-dated the last glacial maximum. Coalescent and non-coalescent demographic analyses indicate that independent lineages currently occupying previously glaciated or unsuitable areas in eastern, central and western US underwent post-glacial population expansion likely from southern refugia during the late Pleistocene/early Holocene. Conversely, southern lineages display patterns consistent with long-term population stability. Such long-term persistence of genetic structure may be due to the competitive effects between lineages or ecosystem stability in more southern latitudes.

The systematic relationships of snakes in the Lampropeltis mexicana complex (L. mexicana, L. alterna, and L. ruthveni) are poorly known despite several taxonomic studies over the last 80 years. Mitochondrial DNA sequences were used to infer the phylogeny of the L. mexicana complex. At least one representative sample from the nine currently recognized species of Lampropeltis was sequenced. Our results suggest that a deep basal split resulted in the divergence of two groups of Lampropeltis, with one group occupying the upland areas of western United States and most of western and central Mexico, and the other northeastern Mexico and the lowland areas of the southern United States. Results also revealed that the L. mexicana complex and Lampropeltis triangulum are polyphyletic, with taxa from both groups nested together in deeply divergent northern and southern clades. These results are incongruent with previous hypotheses of phylogenetic relationships based on morphology, and suggest that morphological characters shared among the various tri-colored Lampropeltis (e.g., hemipenal structure and tri-colored pattern) may be difficult to interpret phylogenetically.

To examine Holarctic snake dispersal, we inferred a phylogenetic tree from four mtDNA genes and one scnDNA gene for most species of the Old World (OW) and New World (NW) colubrid group known as ratsnakes. Ancestral area distributions are estimated for various clades using divergence-vicariance analysis and maximum likelihood on trees produced using Bayesian inference. Dates of divergence for the same clades are estimated using penalized likelihood with statistically crosschecked calibration references obtained from the Miocene fossil record. With ancestral areas and associated dates estimated, various hypotheses concerning the age and environment associated with the origin of ratsnakes and the dispersal of NW taxa from OW ancestors were tested. Results suggest that the ratsnakes originated in tropical Asia in the late Eocene and subsequently dispersed to the Western and Eastern Palearctic by the early Oligocene. These analyses also suggest that the monophyletic NW ratsnakes (the Lampropeltini) diverged from OW ratsnakes and dispersed through Beringia in the late Oligocene/early Miocene when this land bridge was mostly composed of deciduous and coniferous forests.

Snakes of the subfamily Boinae are found in Madagascar, the Papuan-Pacific Islands, and the Neotropics. It has been suggested that genera within each of these particular areas do not form monophyletic groups. Further, it was proposed that the New World Boa constrictor is more closely related to boine genera in Madagascar than to boines in the Neotropics. Along with inferring the relationship of all boine genera using data from the cytochrome b gene and morphology, the placement of Boa was also examined. Phylogenetic inferences using maximum likelihood and Bayesian (BI) methods for combined data analyses and separate analyses of DNA sequence and morphological data were conducted. Priors, parametric bootstraps, and the Shimodaira-Hasegawa test were used to examine the previously proposed placement of Boa with Madagascan taxa using these DNA data. DNA data and combined data analyses strongly reject the hypothesis that Boa is more closely related to Old World genera than to other New World genera. Additionally, strong tree support suggests that all species within Madagascar, the Papuan-Pacific Islands, and the Neotropics each form a monophyletic group with respect to their geographic region.

Albinism is an uncommon phenomenon and inherited condition in animals characterized by a partial or complete lack of melanin. The family Xenodermidae Gray, 1849, is a group of caenophidian snakes widely distributed in South, East, and Southeast Asia, including five recognized genera and 36 species. However, there are currently no reports of albinism in any species in Xenodermidae. Achalinus sheni Ma, Xu, Qi, Wang, Tang, Huang & Jiang, 2023 was first described based on five male specimens from Loudi City and Nanyue District, Hunan Province, China. At the time, there were no descriptions on female individuals. In this study, we report in detail a collected albinistic specimen of A. sheni, which is the first discovery of wild albinism in the family Xenodermidae. We also provide photographs and descriptions of the first three female specimens of A. sheni and extend the diagnosis of this species. [ABSTRACT FROM AUTHOR]

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. [ABSTRACT FROM AUTHOR]

Background: Using molecular data generated by high throughput next generation sequencing (NGS) platforms to infer phylogeny is becoming common as costs go down and the ability to capture loci from across the genome goes up. While there is a general consensus that greater numbers of independent loci should result in more robust phylogenetic estimates, few studies have compared phylogenies resulting from smaller datasets for commonly used genetic markers with the large datasets captured using NGS. Here, we determine how a 5-locus Sanger dataset compares with a 377-locus anchored genomics dataset for understanding the evolutionary history of the pseudoxyrhophiine snake radiation centered in Madagascar. The Pseudoxyrhophiinae comprise ∼86 % of Madagascar's serpent diversity, yet they are poorly known with respect to ecology, behavior, and systematics. Using the 377-locus NGS dataset and the summary statistics species-tree methods STAR and MP-EST, we estimated a well-supported species tree that provides new insights concerning intergeneric relationships for the pseudoxyrhophiines. We also compared how these and other methods performed with respect to estimating tree topology using datasets with varying numbers of loci. Methods: Using Sanger sequencing and an anchored phylogenomics approach, we sequenced datasets comprised of 5 and 377 loci, respectively, for 23 pseudoxyrhophiine taxa. For each dataset, we estimated phylogenies using both gene-tree (concatenation) and species-tree (STAR, MP-EST) approaches. We determined the similarity of resulting tree topologies from the different datasets using Robinson-Foulds distances. In addition, we examined how subsets of these data performed compared to the complete Sanger and anchored datasets for phylogenetic accuracy using the same tree inference methodologies, as well as the program *BEAST to determine if a full coalescent model for species tree estimation could generate robust results with fewer loci compared to the summary statistics species tree approaches. We also examined the individual gene trees in comparison to the 377-locus species tree using the program MetaTree. (Continued on next page) Results: Using the full anchored dataset under a variety of methods gave us the same, well-supported phylogeny for pseudoxyrhophiines. The African pseudoxyrhophiine Duberria is the sister taxon to the Malagasy pseudoxyrhophiines genera, providing evidence for a monophyletic radiation in Madagascar. In addition, within Madagascar, the two major clades inferred correspond largely to the aglyphous and opisthoglyphous genera, suggesting that feeding specializations associated with tooth venom delivery may have played a major role in the early diversification of this radiation. The comparison of tree topologies from the concatenated and species-tree methods using different datasets indicated the 5-locus dataset cannot beused to infer a correct phylogeny for the pseudoxyrhophiines under any method tested here and that summary statistics methods require 50 or more loci to consistently recover the species-tree inferred using the complete anchored dataset. However, as few as 15 loci may infer the correct topology when using the full coalescent species tree method BEAST. MetaTree analyses of each gene tree from the Sanger and anchored datasets found that none of the individual gene trees matched the 377-locus species tree, and that no gene trees were identical with respect to topology. Conclusions: Our results suggest that ≥50 loci may be necessary to confidently infer phylogenies when using summaryspecies-tree methods, but that the coalescent-based method BEAST consistently recovers the same topology using only 15 loci. These results reinforce that datasets with small numbers of markers may result in misleading topologies, and further, that the method of inference used to generate a phylogeny also has a major influence on the number of loci necessary to infer robust species trees. [ABSTRACT FROM AUTHOR]

The genus Pseustes Fitzinger, 1843 is composed of three recognized species, Pseustes poecilonotus, P. shropshirei and P. sulphureus, which may be the largest sized colubrid snake in the New World. The group has a complex systematic history that has yet to be untangled using modern molecular phylogenetic approaches. The systematic position, within-group diversity and distribution are therefore uncertain. We obtained samples of four species from multiple specimens across their distribution and analysed one nuclear and two mitochondrial genes to determine the phylogenetic placement of the genus and infer relationships among Pseustes lineages. We find strong support for the paraphyly of Pseustes with respect to the monotypic genus Spilotes, both of which are nested within a clade of at least 23 other New World Colubrinae genera. Based on our results, we formally revise the taxonomy of P. poecilonotus and P. sulphureus, resurrecting the taxon P. polylepis for populations of P. poecilonotus from South America and allocating P. sulphureus to the genus Spilotes which renders both genera monophyletic. Additionally, we identify two lineages that are putatively new and currently unrecognized species. Finally, the placement of P. sulphureus, the type species of Pseustes, in the genus Spilotes, requires the allocation of the senior synonym Phrynonax be considered for the remaining Pseustes taxa. [ABSTRACT FROM AUTHOR]

Aim Island biogeography predicts positive relationships between area and species richness resulting from the ecological processes of colonization in smaller areas and the evolutionary process of in situ diversification in larger areas. However, the relative contribution of these ecological and evolutionary processes to diversity on islands over time is unclear without associated measures of phylogenetic diversity ( PD), which account for the evolutionary distinctiveness of component taxa. Location Global. Methods We assemble a global dataset of island richness for snakes, along with a phylogeny for the group, to examine the origins of species richness and PD on islands. Results Ecological models driven by colonization explain a majority of species richness in island snake faunas, while in situ diversification is rare, restricted to a few large islands and rarely contributing significantly to island richness. Similarly, PD on islands is tied to variables that affect colonization probability such as isolation and climate, but not area, indicating that in situ diversification does not exert a significant influence on island richness. Main conclusions Even large islands with high in situ diversification do not produce phylogenetically distinct faunas, but are instead dependent on colonization of new lineages from the continent, as predicted by ecological theory. Evolutionary processes such as in situ diversification are relatively insignificant drivers of both species richness and PD compared with colonization from regional species pools, even over evolutionary time-scales on large islands. [ABSTRACT FROM AUTHOR]

Aim We determine whether trait-driven diversification yields similar body size distributions for snakes in local, regional and phylogenetic assemblages. Location United States, North America. Methods Using total length and mass, we examine body size frequency distributions ( BSFD) across 79 sites and respective biomes to determine if these areas represent random subsamples from the source pools of taxon body sizes. Using QuaSSE, we determine if the most probable model of trait-driven diversification in the three most common groups of snakes in North America, the ratsnakes, pitvipers and watersnakes, is similar to the predicted regional BSFD. Results BSFD of snakes at the community, biome, regional and clade scales show symmetric distributions of body size. These patterns may simply be generated from random statistical subsampling. Speciation rates are not highest at or near the modal body size and simulations show that linear trait-driven models can still yield highly symmetric distributions of body size. Main conclusions In this study region, processes such as competition due to size do not alter BSFD from one scale to the other. It is likely that rates of speciation are not highest at the mode in snakes and trait-driven diversification is not likely to account for a regional pool of body sizes from which local communities are drawn, although persistence of modal body sizes through time could yield regional BSFD. [ABSTRACT FROM AUTHOR]

Most phylogeographic studies have used maximum likelihood or maximum parsimony to infer phylogeny and bootstrap analysis to evaluate support for trees. Recently, Bayesian methods using Marlov chain Monte Carlo to search tree space and simultaneously estimate tree support have become popular due to its fast search speed and ability to create a posterior distribution of parameters of interest. Here, I present a study that utilizes Bayesian methods to infer phylogenetic relationships of the cornsnake (Elaphe guttata) complex using cytochrome b sequences. Examination of the posterior probability distributions confirms the existence of three geographic lineages. Additionally, there is no support for the monophyly of the subspecies of E. guttata. Results suggest the three geographic lineages partially conform to the ranges of previously defined subspecies, although Shimodaira–Hasegawa tests suggest that subspecies-constrained trees produce significantly poorer likelihood estimates than the most likely trees reflecting the evolution of three geographic assemblages. Based on molecular support, these three geographic assemblages are recognized as species using evolutionary species criteria: E. guttata, Elaphe slowinskii, and Elaphe emoryi [phylogeographic, maximum likelihood, maximum parsimony, bootstrap, Bayesian, Markov chain Monte Carlo, cornsnake, Cytochrome b, geographic lineages, E. guttta, E. slowinskii, and E. emoryi]. [Copyright &y& Elsevier]

[Display omitted] • UCE data resolves previously problematic relationships in Thamnophiini. • Support for reversing taxonomic changes of Haldea striatula and Liodytes pygaea. • Thamnophis is the only major genus to have crossed westward multiple times. • High level of gene tree discordance with species tree despite high nodal support. North American Thamnophiini (gartersnakes, watersnakes, brownsnakes, and swampsnakes) are an ecologically and phenotypically diverse temperate clade of snakes representing 61 species across 10 genera. In this study, we estimate phylogenetic trees using ∼3,700 ultraconserved elements (UCEs) for 76 specimens representing 75% of all Thamnophiini species. We infer phylogenies using multispecies coalescent methods and time calibrate them using the fossil record. We also conducted ancestral area estimation to identify how major biogeographic boundaries in North America affect broadscale diversification in the group. While most nodes exhibited strong statistical support, analysis of concordant data across gene trees reveals substantial heterogeneity. Ancestral area estimation demonstrated that the genus Thamnophis was the only taxon in this subfamily to cross the Western Continental Divide, even as other taxa dispersed southward toward the tropics. Additionally, levels of gene tree discordance are overall higher in transition zones between bioregions, including the Rocky Mountains. Therefore, the Western Continental Divide may be a significant transition zone structuring the diversification of Thamnophiini during the Neogene and Pleistocene. Here we show that despite high levels of discordance across gene trees, we were able to infer a highly resolved and well-supported phylogeny for Thamnophiini, which allows us to understand broadscale patterns of diversity and biogeography. [ABSTRACT FROM AUTHOR]

[Display omitted] • >4500 ultraconserved elements loci resolved the phylogeny of the snake superfamily Elapoidea. • Anatomical data and phylogeny supported establishment of a new Afro-Asian family of snakes, Micrelapidae. • Different phylogeny estimation methods inferred largely congruent topologies. • Time calibrated phylogeny supported the hypothesis of a rapid elapoid radiation in the Eocene. The highly diverse snake superfamily Elapoidea is considered to be a classic example of ancient, rapid radiation. Such radiations are challenging to fully resolve phylogenetically, with the highly diverse Elapoidea a case in point. Previous attempts at inferring a phylogeny of elapoids produced highly incongruent estimates of their evolutionary relationships, often with very low statistical support. We sought to resolve this situation by sequencing over 4,500 ultraconserved element loci from multiple representatives of every elapoid family/subfamily level taxon and inferring their phylogenetic relationships with multiple methods. Concatenation and multispecies coalescent based species trees yielded largely congruent and well-supported topologies. Hypotheses of a hard polytomy were not retained for any deep branches. Our phylogenies recovered Cyclocoridae and Elapidae as diverging early within Elapoidea. The Afro-Malagasy radiation of elapoid snakes, classified as multiple subfamilies of an inclusive Lamprophiidae by some earlier authors, was found to be monophyletic in all analyses. The genus Micrelaps was consistently recovered as sister to Lamprophiidae. We establish a new family, Micrelapidae fam. nov. , for Micrelaps and assign Brachyophis to this family based on cranial osteological synapomorphy. We estimate that Elapoidea originated in the early Eocene and rapidly diversified into all the major lineages during this epoch. Ecological opportunities presented by the post-Cretaceous-Paleogene mass extinction event may have promoted the explosive radiation of elapoid snakes. [ABSTRACT FROM AUTHOR]