Explore the vast range of titles available.

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.

Single nucleotide polymorphisms (SNPs) are useful markers for phylogenetic studies owing in part to their ubiquity throughout the genome and ease of collection. Restriction site associated DNA sequencing (RADseq) methods are becoming increasingly popular for SNP data collection, but an assessment of the best practises for using these data in phylogenetics is lacking. We use computer simulations, and new double digest RADseq (ddRADseq) data for the lizard family Phrynosomatidae, to investigate the accuracy of RAD loci for phylogenetic inference. We compare the two primary ways RAD loci are used during phylogenetic analysis, including the analysis of full sequences (i.e., SNPs together with invariant sites), or the analysis of SNPs on their own after excluding invariant sites. We find that using full sequences rather than just SNPs is preferable from the perspectives of branch length and topological accuracy, but not of computational time. We introduce two new acquisition bias corrections for dealing with alignments composed exclusively of SNPs, a conditional likelihood method and a reconstituted DNA approach. The conditional likelihood method conditions on the presence of variable characters only (the number of invariant sites that are unsampled but known to exist is not considered), while the reconstituted DNA approach requires the user to specify the exact number of unsampled invariant sites prior to the analysis. Under simulation, branch length biases increase with the amount of missing data for both acquisition bias correction methods, but branch length accuracy is much improved in the reconstituted DNA approach compared to the conditional likelihood approach. Phylogenetic analyses of the empirical data using concatenation or a coalescent-based species tree approach provide strong support for many of the accepted relationships among phrynosomatid lizards, suggesting that RAD loci contain useful phylogenetic signal across a range of divergence times despite the presence of missing data. Phylogenetic analysis of RAD loci requires careful attention to model assumptions, especially if downstream analyses depend on branch lengths.

Anolis lizards (anoles) are textbook study organisms in evolution and ecology. Although several topics in evolutionary biology have been elucidated by the study of anoles, progress in some areas has been hampered by limited phylogenetic information on this group. Here, we present a phylogenetic analysis of all 379 extant species of Anolis, with new phylogenetic data for 139 species including new DNA data for 101 species. We use the resulting estimates as a basis for defining anole clade names under the principles of phylogenetic nomenclature and to examine the biogeographic history of anoles. Our new taxonomic treatment achieves the supposed advantages of recent subdivisions of anoles that employed ranked Linnaean-based nomenclature while avoiding the pitfalls of those approaches regarding artificial constraints imposed by ranks. Our biogeographic analyses demonstrate complexity in the dispersal history of anoles, including multiple crossings of the Isthmus of Panama, two invasions of the Caribbean, single invasions to Jamaica and Cuba, and a single evolutionary dispersal from the Caribbean to the mainland that resulted in substantial anole diversity. Our comprehensive phylogenetic estimate of anoles should prove useful for rigorous testing of many comparative evolutionary hypotheses.

The correlation between species diversification and morphological evolution has long been of interest in evolutionary biology. We investigated the relationship between these processes during the radiation of 250+ scincid lizards that constitute Australia's most species-rich clade of terrestrial vertebrates. We generated a time-calibrated phylogenetic tree for the group that was more than 85% complete at the species level and collected multivariate morphometric data for 183 species. We reconstructed the dynamics of species diversification and trait evolution using a Bayesian statistical framework (BAMM) that simultaneously accounts for variation in evolutionary rates through time and among lineages. We extended the BAMM model to accommodate time-dependent phenotypic evolution, and we describe several new methods for summarizing and visualizing macroevolutionary rate heterogeneity on phylogenetic trees. Two major clades (Lerista, Ctenotus; > 90 spp. each) are associated with high rates of species diversification relative to the background rate across Australian sphenomorphine skinks. The Lerista clade is characterized by relatively high lability of body form and has undergone repeated instances of limb reduction, but Ctenotus is characterized by an extreme deceleration in the rate of body shape evolution. We estimate that rates of phenotypic evolution decreased by more than an order of magnitude in the common ancestor of the Ctenotus clade. These results provide evidence for a modal shift in phenotypic evolutionary dynamics and demonstrate that major axes of morphological variation can be decoupled from species diversification. More generally, the Bayesian framework described here can be used to identify and characterize complex mixtures of dynamic processes on phylogenetic trees.

As we collect range-wide genetic data for morphologically-defined species, we increasingly unearth evidence for cryptic diversity. Delimiting this cryptic diversity is challenging, both because the divergences span a continuum and because the lack of overt morphological differentiation suggests divergence has proceeded heterogeneously. Herein, we address these challenges as we diagnose and describe species in three co-occurring species groups of Australian lizards. By integrating genomic and morphological data with data on hybridization and introgression from contact zones, we explore several approaches—and their relative benefits and weaknesses—for testing the validity of cryptic lineages. More generally, we advocate that genetic delimitations of cryptic diversity must consider whether these lineages are likely to be durable and persistent through evolutionary time.

Broad paradigms of vertebrate genomic repeat element evolution have been largely shaped by analyses of mammalian and avian genomes. Here, based on analyses of genomes sequenced from over 60 squamate reptiles (lizards and snakes), we show that patterns of genomic repeat landscape evolution in squamates challenge such paradigms. Despite low variance in genome size, squamate genomes exhibit surprisingly high variation among species in abundance (ca. 25–73% of the genome) and composition of identifiable repeat elements. We also demonstrate that snake genomes have experienced microsatellite seeding by transposable elements at a scale unparalleled among eukaryotes, leading to some snake genomes containing the highest microsatellite content of any known eukaryote. Our analyses of transposable element evolution across squamates also suggest that lineage-specific variation in mechanisms of transposable element activity and silencing, rather than variation in species-specific demography, may play a dominant role in driving variation in repeat element landscapes across squamate phylogeny. Large-scale patterns of genomic repeat element evolution have been studied mainly in birds and mammals. Here, the authors analyze the genomes of over 60 squamate reptiles and show high variation in repeat elements compared to mammals and birds, and particularly high microsatellite seeding in snakes.

Targeted sequence capture is becoming a widespread tool for generating large phylogenomic data sets to address difficult phylogenetic problems. However, this methodology often generates data sets in which increasing the number of taxa and loci increases amounts of missing data. Thus, a fundamental (but still unresolved) question is whether sampling should be designed to maximize sampling of taxa or genes, or to minimize the inclusion of missing data cells. Here, we explore this question for an ancient, rapid radiation of lizards, the pleurodont iguanians. Pleurodonts include many well-known clades (e.g., anoles, basilisks, iguanas, and spiny lizards) but relationships among families have proven difficult to resolve strongly and consistently using traditional sequencing approaches. We generated up to 4921 ultraconserved elements with sampling strategies including 16, 29, and 44 taxa, from 1179 to approximately 2.4 million characters per matrix and approximately 30% to 60% total missing data. We then compared mean branch support for interfamilial relationships under these 15 different sampling strategies for both concatenated (maximum likelihood) and species tree (NJst) approaches (after showing that mean branch support appears to be related to accuracy). We found that both approaches had the highest support when including loci with up to 50% missing taxa (matrices with ~40-55% missing data overall). Thus, our results show that simply excluding all missing data may be highly problematic as the primary guiding principle for the inclusion or exclusion of taxa and genes. The optimal strategy was somewhat different for each approach, a pattern that has not been shown previously. For concatenated analyses, branch support was maximized when including many taxa (44) but fewer characters (1.1 million). For species-tree analyses, branch support was maximized with minimal taxon sampling (16) but many loci (4789 of 4921). We also show that the choice of these sampling strategies can be critically important for phylogenomic analyses, since some strategies lead to demonstrably incorrect inferences (using the same method) that have strong statistical support. Our preferred estimate provides strong support for most interfamilial relationships in this important but phylogenetically challenging group.

Most existing methods for modeling trait evolution are univariate, although researchers are often interested in investigating evolutionary patterns and processes across multiple traits. Principal components analysis (PCA) is commonly used to reduce the dimensionality of multivariate data so that univariate trait models can be fit to individual principal components. The problem with using standard PCA on phylogenetically structured data has been previously pointed out yet it continues to be widely used in the literature. Here we demonstrate precisely how using standard PCA can mislead inferences: The first few principal components of traits evolved under constant-rate multivariate Brownian motion will appear to have evolved via an \"early burst\" process. A phylogenetic PCA (pPCA) has been proprosed to alleviate these issues. However, when the true model of trait evolution deviates from the model assumed in the calculation of the pPCA axes, we find that the use of pPCA suffers from similar artifacts as standard PCA. We show that data sets with high effective dimensionality are particularly likely to lead to erroneous inferences. Ultimately, all of the problems we report stem from the same underlying issue—by considering only the first few principal components as univariate traits, we are effectively examining a biased sample of a multivariate pattern. These results highlight the need for truly multivariate phylogenetic comparative methods. As these methods are still being developed, we discuss potential alternative strategies for using and interpreting models fit to univariate axes of multivariate data.

Biological invasions are both a pressing environmental challenge and an opportunity to investigate fundamental ecological processes, such as the role of top predators in regulating biodiversity and food-web structure. In whole-ecosystem manipulations of small Caribbean islands on which brown anole lizards ( Anolis sagrei ) were the native top predator, we experimentally staged invasions by competitors (green anoles, Anolis smaragdinus ) and/or new top predators (curly-tailed lizards, Leiocephalus carinatus ). We show that curly-tailed lizards destabilized the coexistence of competing prey species, contrary to the classic idea of keystone predation. Fear-driven avoidance of predators collapsed the spatial and dietary niche structure that otherwise stabilized coexistence, which intensified interspecific competition within predator-free refuges and contributed to the extinction of green-anole populations on two islands. Moreover, whereas adding either green anoles or curly-tailed lizards lengthened food chains on the islands, adding both species reversed this effect—in part because the apex predators were trophic omnivores. Our results underscore the importance of top-down control in ecological communities, but show that its outcomes depend on prey behaviour, spatial structure, and omnivory. Diversity-enhancing effects of top predators cannot be assumed, and non-consumptive effects of predation risk may be a widespread constraint on species coexistence. Whole-ecosystem manipulations of Caribbean islands occupied by brown anoles, involving the addition of competitors (green anoles) and/or top predators (curly-tailed lizards), demonstrate that predator introductions can alter the ecological niches and destabilize the coexistence of competing prey species.

Squamate reptiles (lizards and snakes) are a pivotal group whose relationships have become increasingly controversial. Squamates include >9000 species, making them the second largest group of terrestrial vertebrates. They are important medicinally and as model systems for ecological and evolutionary research. However, studies of squamate biology are hindered by uncertainty over their relationships, and some consider squamate phylogeny unresolved, given recent conflicts between molecular and morphological results. To resolve these conflicts, we expand existing morphological and molecular datasets for squamates (691 morphological characters and 46 genes, for 161 living and 49 fossil taxa, including a new set of 81 morphological characters and adding two genes from published studies) and perform integrated analyses. Our results resolve higher-level relationships as indicated by molecular analyses, and reveal hidden morphological support for the molecular hypothesis (but not vice-versa). Furthermore, we find that integrating molecular, morphological, and paleontological data leads to surprising placements for two major fossil clades (Mosasauria and Polyglyphanodontia). These results further demonstrate the importance of combining fossil and molecular information, and the potential problems of estimating the placement of fossil taxa from morphological data alone. Thus, our results caution against estimating fossil relationships without considering relevant molecular data, and against placing fossils into molecular trees (e.g. for dating analyses) without considering the possible impact of molecular data on their placement.