Discerning structure versus speciation in phylogeographic analysis of Seepage Salamanders (Desmognathus aeneus) using demography, environment, geography, and phenotype.

Numerous mechanisms can drive speciation, including isolation by adaptation, distance, and environment. These forces can promote genetic and phenotypic differentiation of local populations, the formation of phylogeographic lineages, and ultimately, completed speciation. However, conceptually similar mechanisms may also result in stabilizing rather than diversifying selection, leading to lineage integration and the long-term persistence of population structure within genetically cohesive species. Processes that drive the formation and maintenance of geographic genetic diversity while facilitating high rates of migration and limiting phenotypic differentiation may thereby result in population genetic structure that is not accompanied by reproductive isolation. We suggest that this framework can be applied more broadly to address the classic dilemma of "structure" versus "species" when evaluating phylogeographic diversity, unifying population genetics, species delimitation, and the underlying study of speciation. We demonstrate one such instance in the Seepage Salamander (Desmognathus aeneus) from the southeastern United States. Recent studies estimated up to 6.3% mitochondrial divergence and four phylogenomic lineages with broad admixture across geographic hybrid zones, which could potentially represent distinct species supported by our species-delimitation analyses. However, while limited dispersal promotes substantial isolation by distance, microhabitat specificity appears to yield stabilizing selection on a single, uniform, ecologically mediated phenotype. As a result, climatic cycles promote recurrent contact between lineages and repeated instances of high migration through time. Subsequent hybridization is apparently not counteracted by adaptive differentiation limiting introgression, leaving a single unified species with deeply divergent phylogeographic lineages that nonetheless do not appear to represent incipient species.

Speciation Hypotheses from Phylogeographic Delimitation Yield an Integrative Taxonomy for Seal Salamanders (Desmognathus monticola).

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

Congruence and Conflict in the Higher-Level Phylogenetics of Squamate Reptiles: An Expanded Phylogenomic Perspective.

Genome-scale data have the potential to clarify phylogenetic relationships across the tree of life but have also revealed extensive gene tree conflict. This seeming paradox, whereby larger data sets both increase statistical confidence and uncover significant discordance, suggests that understanding sources of conflict is important for accurate reconstruction of evolutionary history. We explore this paradox in squamate reptiles, the vertebrate clade comprising lizards, snakes, and amphisbaenians. We collected an average of 5103 loci for 91 species of squamates that span higher-level diversity within the clade, which we augmented with publicly available sequences for an additional 17 taxa. Using a locus-by-locus approach, we evaluated support for alternative topologies at 17 contentious nodes in the phylogeny. We identified shared properties of conflicting loci, finding that rate and compositional heterogeneity drives discordance between gene trees and species tree and that conflicting loci rarely overlap across contentious nodes. Finally, by comparing our tests of nodal conflict to previous phylogenomic studies, we confidently resolve 9 of the 17 problematic nodes. We suggest this locus-by-locus and node-by-node approach can build consensus on which topological resolutions remain uncertain in phylogenomic studies of other contentious groups. [Anchored hybrid enrichment (AHE); gene tree conflict; molecular evolution; phylogenomic concordance; target capture; ultraconserved elements (UCE).].

Phylogenomics Reveals Ancient Gene Tree Discordance in the Amphibian Tree of Life.

Molecular phylogenies have yielded strong support for many parts of the amphibian Tree of Life, but poor support for the resolution of deeper nodes, including relationships among families and orders. To clarify these relationships, we provide a phylogenomic perspective on amphibian relationships by developing a taxon-specific Anchored Hybrid Enrichment protocol targeting hundreds of conserved exons which are effective across the class. After obtaining data from 220 loci for 286 species (representing 94% of the families and 44% of the genera), we estimate a phylogeny for extant amphibians and identify gene tree-species tree conflict across the deepest branches of the amphibian phylogeny. We perform locus-by-locus genealogical interrogation of alternative topological hypotheses for amphibian monophyly, focusing on interordinal relationships. We find that phylogenetic signal deep in the amphibian phylogeny varies greatly across loci in a manner that is consistent with incomplete lineage sorting in the ancestral lineage of extant amphibians. Our results overwhelmingly support amphibian monophyly and a sister relationship between frogs and salamanders, consistent with the Batrachia hypothesis. Species tree analyses converge on a small set of topological hypotheses for the relationships among extant amphibian families. These results clarify several contentious portions of the amphibian Tree of Life, which in conjunction with a set of vetted fossil calibrations, support a surprisingly younger timescale for crown and ordinal amphibian diversification than previously reported. More broadly, our study provides insight into the sources, magnitudes, and heterogeneity of support across loci in phylogenomic data sets.[AIC; Amphibia; Batrachia; Phylogeny; gene tree-species tree discordance; genomics; information theory.].

Evolutionary legacies in contemporary tetrapod imperilment.

The Tree of Life will be irrevocably reshaped as anthropogenic extinctions continue to unfold. Theory suggests that lineage evolutionary dynamics, such as age since origination, historical extinction filters and speciation rates, have influenced ancient extinction patterns - but whether these factors also contribute to modern extinction risk is largely unknown. We examine evolutionary legacies in contemporary extinction risk for over 4000 genera, representing ~30,000 species, from the major tetrapod groups: amphibians, birds, turtles and crocodiles, squamate reptiles and mammals. We find consistent support for the hypothesis that extinction risk is elevated in lineages with higher recent speciation rates. We subsequently test, and find modest support for, a primary mechanism driving this pattern: that rapidly diversifying clades predominantly comprise range-restricted, and extinction-prone, species. These evolutionary patterns in current imperilment may have important consequences for how we manage the erosion of biological diversity across the Tree of Life.

Amphibian Speciation Rates Support a General Role of Mountains as Biodiversity Pumps.

AbstractContinental mountain areas cover <15% of global land surface, yet these regions concentrate >80% of global terrestrial diversity. One prominent hypothesis to explain this pattern proposes that high mountain diversities could be explained by higher diversification rates in regions of high topographic complexity (HTC). While high speciation in mountains has been detected for particular clades and regions, the global extent to which lineages experience faster speciation in mountains remains unknown. Here we addressed this issue using amphibians as a model system (>7,000 species), and we found that families showing high speciation rates contain a high proportion of species distributed in mountains. Moreover, we found that lineages inhabiting areas of HTC speciate faster than lineages occupying areas that are topographically less complex. When comparing across regions, we identified the same pattern in five biogeographical realms where higher speciation rates are associated with higher levels of complex topography. Low-magnitude differences in speciation rates between some low and high complex topographies suggest that high mountain diversity is also affected by low extinction and/or high colonization rates. Nevertheless, our results bolster the importance of mountains as engines of speciation at different geographical scales and highlight their importance for the conservation of global biodiversity.

Interrogating Genomic-Scale Data for Squamata (Lizards, Snakes, and Amphisbaenians) Shows no Support for Key Traditional Morphological Relationships.

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

Phylogenetic and spatial distribution of evolutionary diversification, isolation, and threat in turtles and crocodilians (non-avian archosauromorphs).

BACKGROUND: The origin of turtles and crocodiles and their easily recognized body forms dates to the Triassic and Jurassic. Despite their long-term success, extant species diversity is low, and endangerment is extremely high compared to other terrestrial vertebrate groups, with ~ 65% of ~ 25 crocodilian and ~ 360 turtle species now threatened by exploitation and habitat loss. Here, we combine available molecular and morphological evidence with statistical and machine learning algorithms to present a phylogenetically informed, comprehensive assessment of diversification, threat status, and evolutionary distinctiveness of all extant species. RESULTS: In contrast to other terrestrial vertebrates and their own diversity in the fossil record, the recent extant lineages of turtles and crocodilians have not experienced any global mass extinctions or lineage-wide shifts in diversification rate or body-size evolution over time. We predict threat statuses for 114 as-yet unassessed or data-deficient species and identify a concentration of threatened turtles and crocodilians in South and Southeast Asia, western Africa, and the eastern Amazon. We find that unlike other terrestrial vertebrate groups, extinction risk increases with evolutionary distinctiveness: a disproportionate amount of phylogenetic diversity is concentrated in evolutionarily isolated, at-risk taxa, particularly those with small geographic ranges. Our findings highlight the important role of geographic determinants of extinction risk, particularly those resulting from anthropogenic habitat-disturbance, which affect species across body sizes and ecologies. CONCLUSIONS: Extant turtles and crocodilians maintain unique, conserved morphologies which make them globally recognizable. Many species are threatened due to exploitation and global change. We use taxonomically complete, dated molecular phylogenies and various approaches to produce a comprehensive assessment of threat status and evolutionary distinctiveness of both groups. Neither group exhibits significant overall shifts in diversification rate or body-size evolution, or any signature of global mass extinctions in recent, extant lineages. However, the most evolutionarily distinct species tend to be the most threatened, and species richness and extinction risk are centered in areas of high anthropogenic disturbance, particularly South and Southeast Asia. Range size is the strongest predictor of threat, and a disproportionate amount of evolutionary diversity is at risk of imminent extinction.

Phylogenomic data reveal reticulation and incongruence among mitochondrial candidate species in Dusky Salamanders (Desmognathus).

Gene flow between evolutionarily distinct lineages is increasingly recognized as a common occurrence. Such processes distort our ability to diagnose and delimit species, as well as confound attempts to estimate phylogenetic relationships. A conspicuous example is Dusky Salamanders (Desmognathus), a common model-system for ecology, evolution, and behavior. Only 22 species are described, 7 in the last 40 years. However, mitochondrial datasets indicate the presence of up to 45 "candidate species" and multiple paraphyletic taxa presenting a complex history of reticulation. Some authors have even suggested that the search for species boundaries in the group may be in vain. Here, we analyze nuclear and mitochondrial data containing 161 individuals from at least 49 distinct evolutionary lineages that we treat as candidate species. Concatenated and species-tree methods do not estimate fully resolved relationships among these taxa. Comparing topologies and applying methods for estimating phylogenetic networks, we find strong support for numerous instances of hybridization throughout the history of the group. We suggest that these processes may be more common than previously thought across the phylogeography-phylogenetics continuum, and that while the search for species boundaries in Desmognathus may not be in vain, it will be complicated by factors such as crypsis, parallelism, and gene-flow.

Novel Approaches for Phylogenetic Inference from Morphological Data and Total-Evidence Dating in Squamate Reptiles (Lizards, Snakes, and Amphisbaenians).

Here, I combine previously underutilized models and priors to perform more biologically realistic phylogenetic inference from morphological data, with an example from squamate reptiles. When coding morphological characters, it is often possible to denote ordered states with explicit reference to observed or hypothetical ancestral conditions. Using this logic, we can integrate across character-state labels and estimate meaningful rates of forward and backward transitions from plesiomorphy to apomorphy. I refer to this approach as MkA, for "asymmetric." The MkA model incorporates the biological reality of limited reversal for many phylogenetically informative characters, and significantly increases likelihoods in the empirical data sets. Despite this, the phylogeny of Squamata remains contentious. Total-evidence analyses using combined morphological and molecular data and the MkA approach tend toward recent consensus estimates supporting a nested Iguania. However, support for this topology is not unambiguous across data sets or analyses, and no mechanism has been proposed to explain the widespread incongruence between partitions, or the hidden support for various topologies in those partitions. Furthermore, different morphological data sets produced by different authors contain both different characters and different states for the same or similar characters, resulting in drastically different placements for many important fossil lineages. Effort is needed to standardize ontology for morphology, resolve incongruence, and estimate a robust phylogeny. The MkA approach provides a preliminary avenue for investigating morphological evolution while accounting for temporal evidence and asymmetry in character-state changes.

Patterns, biases and prospects in the distribution and diversity of Neotropical snakes.

Motivation: We generated a novel database of Neotropical snakes (one of the world's richest herpetofauna) combining the most comprehensive, manually compiled distribution dataset with publicly available data. We assess, for the first time, the diversity patterns for all Neotropical snakes as well as sampling density and sampling biases. Main types of variables contained: We compiled three databases of species occurrences: a dataset downloaded from the Global Biodiversity Information Facility (GBIF), a verified dataset built through taxonomic work and specialized literature, and a combined dataset comprising a cleaned version of the GBIF dataset merged with the verified dataset. Spatial location and grain: Neotropics, Behrmann projection equivalent to 1° × 1°. Time period: Specimens housed in museums during the last 150 years. Major taxa studied: Squamata: Serpentes. Software format: Geographical information system (GIS). Results: The combined dataset provides the most comprehensive distribution database for Neotropical snakes to date. It contains 147,515 records for 886 species across 12 families, representing 74% of all species of snakes, spanning 27 countries in the Americas. Species richness and phylogenetic diversity show overall similar patterns. Amazonia is the least sampled Neotropical region, whereas most well-sampled sites are located near large universities and scientific collections. We provide a list and updated maps of geographical distribution of all snake species surveyed. Main conclusions: The biodiversity metrics of Neotropical snakes reflect patterns previously documented for other vertebrates, suggesting that similar factors may determine the diversity of both ectothermic and endothermic animals. We suggest conservation strategies for high-diversity areas and sampling efforts be directed towards Amazonia and poorly known species.

Using phylogenomics to understand the link between biogeographic origins and regional diversification in ratsnakes.

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.

Methodological congruence in phylogenomic analyses with morphological support for teiid lizards (Sauria: Teiidae).

A well-known issue in phylogenetics is discordance among gene trees, species trees, morphology, and other data types. Gene-tree discordance is often caused by incomplete lineage sorting, lateral gene transfer, and gene duplication. Multispecies-coalescent methods can account for incomplete lineage sorting and are believed by many to be more accurate than concatenation. However, simulation studies and empirical data have demonstrated that concatenation and species tree methods often recover similar topologies. We use three popular methods of phylogenetic reconstruction (one concatenation, two species tree) to evaluate relationships within Teiidae. These lizards are distributed across the United States to Argentina and the West Indies, and their classification has been controversial due to incomplete sampling and the discordance among various character types (chromosomes, DNA, musculature, osteology, etc.) used to reconstruct phylogenetic relationships. Recent morphological and molecular analyses of the group resurrected three genera and created five new genera to resolve non-monophyly in three historically ill-defined genera: Ameiva, Cnemidophorus, and Tupinambis. Here, we assess the phylogenetic relationships of the Teiidae using "next-generation" anchored-phylogenomics sequencing. Our final alignment includes 316 loci (488,656bp DNA) for 244 individuals (56 species of teiids, representing all currently recognized genera) and all three methods (ExaML, MP-EST, and ASTRAL-II) recovered essentially identical topologies. Our results are basically in agreement with recent results from morphology and smaller molecular datasets, showing support for monophyly of the eight new genera. Interestingly, even with hundreds of loci, the relationships among some genera in Tupinambinae remain ambiguous (i.e. low nodal support for the position of Salvator and Dracaena).

Predicting community structure in snakes on Eastern Nearctic islands using ecological neutral theory and phylogenetic methods.

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.

Biogeography of worm lizards (Amphisbaenia) driven by end-Cretaceous mass extinction.

Worm lizards (Amphisbaenia) are burrowing squamates that live as subterranean predators. Their underground existence should limit dispersal, yet they are widespread throughout the Americas, Europe and Africa. This pattern was traditionally explained by continental drift, but molecular clocks suggest a Cenozoic diversification, long after the break-up of Pangaea, implying dispersal. Here, we describe primitive amphisbaenians from the North American Palaeocene, including the oldest known amphisbaenian, and provide new and older molecular divergence estimates for the clade, showing that worm lizards originated in North America, then radiated and dispersed in the Palaeogene following the Cretaceous-Palaeogene (K-Pg) extinction. This scenario implies at least three trans-oceanic dispersals: from North America to Europe, from North America to Africa and from Africa to South America. Amphisbaenians provide a striking case study in biogeography, suggesting that the role of continental drift in biogeography may be overstated. Instead, these patterns support Darwin and Wallace's hypothesis that the geographical ranges of modern clades result from dispersal, including oceanic rafting. Mass extinctions may facilitate dispersal events by eliminating competitors and predators that would otherwise hinder establishment of dispersing populations, removing biotic barriers to dispersal.

Contrasting models of parity-mode evolution in squamate reptiles.

Recent analyses using large-scale phylogenies suggest a radically different history for the evolution of live birth and egg laying in squamate reptiles (lizards and snakes) than traditionally understood. What is the ancestral condition for lizards and snakes? How frequently does live bearing evolve in egg-laying lineages? Can the eggshell ever re-evolve in live-bearing lineages? Answering these fundamental questions about the evolution of key physiological processes will require additional data from genomic, developmental, and fossil data.

Biogeographic analysis reveals ancient continental vicariance and recent oceanic dispersal in amphibians.

Amphibia comprises over 7000 extant species distributed in almost every ecosystem on every continent except Antarctica. Most species also show high specificity for particular habitats, biomes, or climatic niches, seemingly rendering long-distance dispersal unlikely. Indeed, many lineages still seem to show the signature of their Pangaean origin, approximately 300 Ma later. To date, no study has attempted a large-scale historical-biogeographic analysis of the group to understand the distribution of extant lineages. Here, I use an updated chronogram containing 3309 species (∼ 45% of extant diversity) to reconstruct their movement between 12 global ecoregions. I find that Pangaean origin and subsequent Laurasian and Gondwanan fragmentation explain a large proportion of patterns in the distribution of extant species. However, dispersal during the Cenozoic, likely across land bridges or short distances across oceans, has also exerted a strong influence. Finally, there are at least three strongly supported instances of long-distance oceanic dispersal between former Gondwanan landmasses during the Cenozoic. Extinction from intervening areas seems to be a strong factor in shaping present-day distributions. Dispersal and extinction from and between ecoregions are apparently tied to the evolution of extraordinarily adaptive expansion-oriented phenotypes that allow lineages to easily colonize new areas and diversify, or conversely, to extremely specialized phenotypes or heavily relictual climatic niches that result in strong geographic localization and limited diversification.

Coalescent species delimitation in milksnakes (genus Lampropeltis) and impacts on phylogenetic comparative analyses.

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.

Early origin of viviparity and multiple reversions to oviparity in squamate reptiles.

Viviparity has putatively evolved 115 times in squamates (lizards and snakes), out of only ~ 140 origins in vertebrates, and is apparently related to colder climates and other factors such as body size. Viviparity apparently evolves from oviparity via egg-retention, and such taxa may thus still have the machinery to produce thick-shelled eggs. Parity mode is also associated with variable diversification rates in some groups. We reconstruct ancestral parity modes accounting for state-dependent diversification in a large-scale phylogenetic analysis, and find strong support for an early origin of viviparity at the base of Squamata, and a complex pattern of subsequent transitions. Viviparous lineages have higher rates of speciation and extinction, and greater species turnover through time. Viviparity is associated with lower environmental and body temperatures in lizards and amphisbaenians, but not female mass. These results suggest that parity mode is a labile trait that shifts frequently in response to ecological conditions.