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.
The relative influence of Neogene geomorphological events and Quaternary climatic changes as causal mechanisms on Neotropical diversification remains largely speculative, as most divergence timing ...inferences are based on a single locus and have limited taxonomic or geographic sampling. To investigate these influences, we use a multilocus (two mitochondrial and 11 nuclear genes) range-wide sampling of Phyllopezus pollicaris, a gecko complex widely distributed across the poorly studied South American 'dry diagonal' biomes. Our approach couples traditional and model-based phylogeography with geospatial methods, and demonstrates Miocene diversification and limited influence of Pleistocene climatic fluctuations on P. pollicaris. Phylogeographic structure and distribution models highlight that persistence across multiple isolated regions shaped the diversification of this species complex. Approximate Bayesian computation supports hypotheses of allopatric and ecological/sympatric speciation between lineages that largely coincide with genetic clusters associated with Chaco, Cerrado, and Caatinga, standing for complex diversification between the 'dry diagonal' biomes. We recover extremely high genetic diversity and suggest that eight well-supported clades may be valid species, with direct implications for taxonomy and conservation assessments. These patterns exemplify how low-vagility species complexes, characterized by strong genetic structure and pre-Pleistocene divergence histories, represent ideal radiations to investigate broad biogeographic histories of associated biomes.
Species delimitation is a major research focus in evolutionary biology because accurate species boundaries are a prerequisite for the study of speciation. New species delimitation methods (SDMs) can ...accommodate nonmonophyletic species and gene tree discordance as a result of incomplete lineage sorting via the coalescent model, but do not explicitly accommodate gene flow after divergence. Approximate Bayesian computation (ABC) can incorporate gene flow and estimate other relevant parameters of the speciation process while testing alternative species delimitation hypotheses. We evaluated the accuracy of BPP, SpeDeSTEM, and ABC for delimiting species using simulated data and applied these methods to empirical data from lizards of the Liolaemus darwinii complex. Overall, BPP was the most accurate, ABC showed an intermediate accuracy, and SpeDeSTEM was the least accurate under most simulated conditions. All three SDMs showed lower accuracy when speciation occurred despite gene flow, as found in previous studies, but ABC was the method with the smallest decrease in accuracy. All three SDMs consistently supported the distinctness of southern and northern lineages within L. darwinii. These SDMs based on genetic data should be complemented with novel SDMs based on morphological and ecological data to achieve truly integrative and statistically robust approaches to species discovery.
Display omitted
► We perform a phylogenetic analysis of iguanian lizards using 29 nuclear loci. ► Most acrodont iguanian relationships are resolved with high support. ► Pleurodont relationships ...recovered with high Bayesian and variable likelihood support. ► Pleurodonts are inferred to have originated in North America and dispersed southward. ►
Anolis removed from Polychrotidae and placed in Dactyloidae.
Iguanian lizards form a diverse clade whose members have been the focus of many comparative studies of ecology, behavior, and evolution. Despite the importance of phylogeny to such studies, interrelationships among many iguanian clades remain uncertain. Within the Old World clade Acrodonta, Agamidae is sometimes found to be paraphyletic with respect to Chamaeleonidae, and recent molecular studies have produced conflicting results for many major clades. Within the largely New World clade Pleurodonta, relationships among the 12 currently recognized major subclades (mostly ranked as families) have been largely unresolved or poorly supported in previous studies. To clarify iguanian evolutionary history, we first infer phylogenies using concatenated maximum-likelihood (ML) and Bayesian analyses of DNA sequence data from 29 nuclear protein-coding genes for 47 iguanian and 29 outgroup taxa. We then estimate a relaxed-clock Bayesian chronogram for iguanians using BEAST. All three methods produce identical topologies. Within Acrodonta, we find strong support for monophyly of Agamidae with respect to Chamaeleonidae, and for almost all relationships within agamids. Within Pleurodonta, we find strong Bayesian support for almost all relationships, and strong ML support for some interfamilial relationships and for monophyly of almost all families (excepting Polychrotidae). Our phylogenetic results suggest a non-traditional biogeographic scenario in which pleurodonts originated in the Northern Hemisphere and subsequently spread southward into South America. The pleurodont portion of the tree is characterized by several very short, deep branches, raising the possibility of deep coalescences that may confound concatenated analyses. We therefore also use 27 of these genes to implement a coalescent-based species-tree approach for pleurodonts. Although this analysis strongly supports monophyly of the pleurodont families, interfamilial relationships are generally different from those in the concatenated tree, and support is uniformly poor. However, a species-tree analysis using only the seven most variable loci yields higher support and more congruence with the concatenated tree. This suggests that low support in the 27-gene species-tree analysis may be an artifact of the many loci that are uninformative for very short branches. This may be a general problem for the application of species-tree methods to rapid radiations, even with phylogenomic data sets. Finally, we correct the non-monophyly of Polychrotidae by recognizing the pleurodont genus
Anolis (
sensu lato) as a separate family (Dactyloidae), and we correct the non-monophyly of the agamid genus
Physignathus by resurrection of the genus
Istiurus for the former
Physignathus lesueurii.
Lizards have been model organisms for ecological and evolutionary studies from individual to community levels at multiple spatial and temporal scales. Here we highlight lizards as models for ...phylogeographic studies, review the published population genetics/phylogeography literature to summarize general patterns and trends and describe some studies that have contributed to conceptual advances. Our review includes 426 references and 452 case studies: this literature reflects a general trend of exponential growth associated with the theoretical and empirical expansions of the discipline. We describe recent lizard studies that have contributed to advances in understanding of several aspects of phylogeography, emphasize some linkages between phylogeography and speciation and suggest ways to expand phylogeographic studies to test alternative pattern-based modes of speciation. Allopatric speciation patterns can be tested by phylogeographic approaches if these are designed to discriminate among four alternatives based on the role of selection in driving divergence between populations, including: (i) passive divergence by genetic drift; (ii) adaptive divergence by natural selection (niche conservatism or ecological speciation); and (iii) socially-mediated speciation. Here we propose an expanded approach to compare patterns of variation in phylogeographic data sets that, when coupled with morphological and environmental data, can be used to to discriminate among these alternative speciation patterns. Correction made after online publication (28/07/2010): (minor deletion in the last line of the abstract).
Molecular data offer great potential to resolve the phylogeny of living taxa but can molecular data improve our understanding of relationships of fossil taxa? Simulations suggest that this is ...possible, but few empirical examples have demonstrated the ability of molecular data to change the placement of fossil taxa. We offer such an example here. We analyze the placement of snakes among squamate reptiles, combining published morphological data (363 characters) and new DNA sequence data (15,794 characters, 22 nuclear loci) for 45 living and 19 fossil taxa. We find several intriguing results. First, some fossil taxa undergo major changes in their phylogenetic position when molecular data are added. Second, most fossil taxa are placed with strong support in the expected clades by the combined data Bayesian analyses, despite each having >98% missing cells and despite recent suggestions that extensive missing data are problematic for Bayesian phylogenetics. Third, morphological data can change the placement of living taxa in combined analyses, even when there is an overwhelming majority of molecular characters. Finally, we find strong but apparently misleading signal in the morphological data, seemingly associated with a burrowing lifestyle in snakes, amphisbaenians, and dibamids. Overall, our results suggest promise for an integrated and comprehensive Tree of Life by combining molecular and morphological data for living and fossil taxa.
Molecular phylogenetics has entered a new era in which species trees are estimated from a collection of gene trees using methods that accommodate their heterogeneity and discordance with the species ...tree. Empirical evaluation of species trees is necessary to assess the performance (i.e., accuracy and precision) of these methods with real data, which consists of gene genealogies likely shaped by different historical and demographic processes. We analyzed 20 loci for 16 species of the South American lizards of the Liolaemus darwinii species group and reconstructed a species tree with *BEAST, then compared the performance of this method under different sampling strategies of loci, individuals, and sequence lengths. We found an increase in the accuracy and precision of species trees with the number of loci, but for any number of loci, accuracy substantially decreased only when using only one individual per species or 25% of the full sequence length (∼ 147 bp). In addition, locus "informativeness" was an important factor in the accuracy/precision of species trees when using a few loci, but it became increasingly irrelevant with additional loci. Our empirical results combined with the previous simulation studies suggest that there is an optimal range of sampling effort of loci, individuals, and sequence lengths for a given speciation history and information content of the data. Future studies should be directed toward further assessment of other factors that can impact performance of species trees, including gene flow, locus "informativeness," tree shape, missing data, and errors in species delimitation.
Display omitted
•We use an integrative approach to test species limits in Patagonian lizards.•We found support for a candidate species.•We detected genome-wide signs of hybridization.•More data is ...needed to assess the extent of hybridization and taxonomic implications.
Recent conceptual and methodological advances have enabled an increasing number of studies to address the problem of species delimitation in a comprehensive manner. This is of particular interest in cases of species whose divergence times are recent and/or effective population sizes are large, where the conclusions obtained from a single source of evidence may lead to erroneous estimations of true species numbers or incorrect assignment of individuals to species. Iguanian lizards of the Liolaemus kingii group (13 species) comprise an important component of the endemic fauna of Patagonia. The southernmost species of this group (namely L. baguali, L. escarchadosi, L. sarmientoi, and L. tari) show widely overlapping distributions across southern Patagonia, also, their phylogenetic relationships are ambiguous and species boundaries have not been explicitly tested. Here we use a comprehensive approach to assess species limits through the use of molecular and morphological information (mitochondrial cytb, nuclear sequences collected by ddRADseq, and linear, meristic and landmark-based morphometrics). We found support for the current taxonomy given that the different analyses recognized the nominal species (4 entities), also a candidate species was supported by mitochondrial and morphological data. In addition, we detected signs of admixture between some of the species. Our results indicate that the L. kingii group can serve as a model system in studies of diversification accompanied by hybridization in nature, which in turn might have been promoted by past climatic oscillations and generalist morphologies. We emphasize the importance of using multiple lines of evidence in order to solve evolutionary stories, and minimizing potential erroneous results that may arise when relying on a single source of information.
Aim: To investigate the historical distribution of the Cerrado across Quaternary climatic fluctuations and to generate historical stability maps to test: (1) whether the 'historical climate' ...stability hypothesis explains squamate reptile richness in the Cerrado; and (2) the hypothesis of Pleistocene connections between savannas located north and south of Amazonia. Location: The Cerrado, a savanna biome and a global biodiversity hotspot distributed mainly in central Brazil. Methods: We generated occurrence datasets from 1000 presence points randomly selected from the entire distribution of the Cerrado, as determined by two spatial definitions. We modelled the potential Cerrado distribution by implementing a maximum-entropy machine-learning algorithm across four time projections: current, mid-Holocene (6 ka), Last Glacial Maximum (LGM, 21 ka) and Last Interglacial (LIG, 120 ka). We generated historical stability maps (refugiai areas) by overlapping presence/absence projections of all scenarios, and checked consistencies with qualitative comparisons with available fossil pollen records. We built a spatially explicit simultaneous autoregressive model to explore the relationship between current climate, climatic stability, and squamate species richness. Results: Models predicted the LGM and LIG as the periods of narrowest and widest Cerrado distributions, respectively, and were largely corroborated by palynological evidence. We found evidence for two savanna corridors (eastern coastal during the LIG, and Andean during the LGM) and predicted a large refugiai area in the north-eastern Cerrado (Serra Geral de Goiás refugium). Variables related to climatic stability predicted squamate richness better than present climatic variables did. Main conclusions: Our results indicate that Bolivian savannas should be included within the Cerrado range and that the Cerrado's biogeographical counterparts are not Chaco and Caatinga but rather the disjunct savannas of the Guyana shield plateaus. Climatic stability is a good predictor of Cerrado squamate richness, and our stability maps could be used in future studies to test diversity patterns and genetic signatures of different taxonomic groups and as a higher-order landscape biodiversity surrogate for conservation planning.
Squamate reptiles (lizards and snakes) are one of the most diverse groups of terrestrial vertebrates. Recent molecular analyses have suggested a very different squamate phylogeny relative to ...morphological hypotheses, but many aspects remain uncertain from molecular data. Here, we analyse higher-level squamate phylogeny with a molecular dataset of unprecedented size, including 161 squamate species for up to 44 nuclear genes each (33 717 base pairs), using both concatenated and species-tree methods for the first time. Our results strongly resolve most squamate relationships and reveal some surprising results. In contrast to most other recent studies, we find that dibamids and gekkotans are together the sister group to all other squamates. Remarkably, we find that the distinctive scolecophidians (blind snakes) are paraphyletic with respect to other snakes, suggesting that snakes were primitively burrowers and subsequently re-invaded surface habitats. Finally, we find that some clades remain poorly supported, despite our extensive data. Our analyses show that weakly supported clades are associated with relatively short branches for which individual genes often show conflicting relationships. These latter results have important implications for all studies that attempt to resolve phylogenies with large-scale phylogenomic datasets.