No abstract

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.

Documented cases of convergent molecular evolution due to selection are fairly unusual, and examples to date have involved only a few amino acid positions. However, because convergence mimics shared ancestry and is not accommodated by current phylogenetic methods, it can strongly mislead phylogenetic inference when it does occur. Here, we present a case of extensive convergent molecular evolution between snake and agamid lizard mitochondrial genomes that overcomes an otherwise strong phylogenetic signal. Evidence from morphology, nuclear genes, and most sites in the mitochondrial genome support one phylogenetic tree, but a subset of mostly amino acid-altering substitutions (primarily at the first and second codon positions) across multiple mitochondrial genes strongly supports a radically different phylogeny. The relevant sites generally evolved slowly but converged between ancient lineages of snakes and agamids. We estimate that Ϸ44 of 113 predicted convergent changes distributed across all 13 mitochondrial protein-coding genes are expected to have arisen from nonneutral causes-a remarkably large number. Combined with strong previous evidence for adaptive evolution in snake mitochondrial proteins, it is likely that much of this convergent evolution was driven by adaptation. These results indicate that nonneutral convergent molecular evolution in mitochondria can occur at a scale and intensity far beyond what has been documented previously, and they highlight the vulnerability of standard phylogenetic methods to the presence of nonneutral convergent sequence evolution.adaptation ͉ convergence ͉ phylogenetics ͉ reptile

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.

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

No abstract

Demonstrations of interactions between diverse selective forces on bright coloration in defended species are rare. Recent work has suggested that not only do the bright colours of Neotropical poison frogs serve to deter predators, but they also play a role in sexual selection, with females preferring males similar to themselves. These studies report an interaction between the selective forces of mate choice and predation. However, evidence demonstrating phenotypic discrimination by potential predators on these polymorphic species is lacking. The possibility remains that visual (avian) predators possess an inherent avoidance of brightly coloured diurnal anurans and purifying selection against novel phenotypes within populations is due solely to non-random mating. Here, we examine the influence of predation on phenotypic variation in a polymorphic species of poison frog, Dendrobates tinctorius. Using clay models, we demonstrate a purifying role for predator selection, as brightly coloured novel forms are more likely to suffer an attack than both local aposematic and cryptic forms. Additionally, local aposematic forms are attacked, though infrequently, indicating ongoing testing/learning and a lack of innate avoidance. These results demonstrate predator-driven phenotypic purification within populations and suggest colour patterns of poison frogs may truly represent a 'magic trait'.

Aim The species-rich family of true toads (Anura: Bufonidae) has been the focus of several earlier studies investigating the biogeography of geographically widespread taxa. Herein, we employ newly developed Bayesian divergence estimate methods to investigate the biogeographical history of this group. Resulting age estimates are used to test several key temporal hypotheses including that the origin of the bufonid clade pre-dates Gondwanan vicariance (~105 million years ago, Ma). Area cladograms are also invoked to investigate the geographical origin of the family.Location Worldwide, except the Australia-New Guinea plate, Madagascar and the Antarctic.Methods A phylogenetic hypothesis of the relationships among true toads was derived from analysis of 2521 bp of DNA data including fragments from three mitochondrial ( 12S , tRNA val , 16S ) and two nuclear ( RAG-1 , CXCR-4 ) genes. Analysis of multiple, unlinked loci with a Bayesian method for estimating divergence times allowed us to address the timing and biogeographical history of Bufonidae. Resulting divergence estimates permitted the investigation of alternative vicariance/dispersal scenarios that have been proposed for true toads. ResultsOur area cladogram resulting from phylogenetic analysis of DNA data supports a South American origin for Bufonidae. Divergence estimates indicate that the family originated earlier than had been suggested previously (78-99 Ma). The age of the enigmatic Caribbean clade was dated to the late Palaeocene-early Eocene. A return of bufonids to the New World in the Eocene was followed by rapid diversification and secondary expansion into South America by the early Oligocene (Rupelian). Main conclusionsThe South American origin of Bufonidae in the Upper Cretaceous was followed by relatively rapid expansion and radiation around the globe, ending with a return to the Americas via a Eurasian/North American land bridge in the Eocene. Though the exact route of this dispersal (Beringia or North Atlantic) remains unclear, an argument is made for the less frequently invoked North Atlantic connection. The origin of the enigmatic Caribbean lineage was found to be consistent with colonization following the bolide impact at the K/T boundary. These findings provide the first, firm foundation for understanding true toad divergence times and their truly remarkable and global radiation.