Sphenodontian phylogeny and the impact of model choice in Bayesian morphological clock estimates of divergence times and evolutionary rates

A juvenile pleurosaurid (Lepidosauria: Rhynchocephalia) from the Tithonian of the Mörnsheim Formation, Germany

Late Jurassic rhynchocephalians from the Solnhofen Archipelago have been known for almost two centuries. The number of specimens and taxa is constantly increasing, but little is known about the ontogeny of these animals. The well-documented marine taxon Pleurosaurus is one of such cases. With over 15 described (and many more undescribed) specimens, there were no unambiguous juveniles so far. Some authors have argued that Acrosaurus, another common component of the Solnhofen Archipelago herpetofauna, might represent an early ontogenetic stage of Pleurosaurus, but the lack of proper descriptions for this taxon makes this assignment tentative, at best. Here, we describe the first unambiguous post-hatchling juvenile of Pleurosaurus and tentatively attribute it to Pleurosaurus cf. P. ginsburgi. The new specimen comes from the Lower Tithonian of the Mörnsheim Formation, Germany. This specimen is small, disarticulated, and incomplete, but preserves several of its craniomandibular bones and presacral vertebrae. It shares with Pleurosaurus a set of diagnostic features, such as an elongated and triangular skull, a low anterior flange in its dentition, and an elongated axial skeleton. It can be identified as a juvenile due to the presence of an unworn dentition, well-spaced posteriormost dentary teeth, a large gap between the last teeth and the coronoid process of the dentary, and poorly ossified vertebrae with unfused neural arches. Acrosaurus shares many anatomical features with both this specimen and Pleurosaurus, which could indicate that the two genera are indeed synonyms. The early ontogenetic stage inferred for the new Pleurosaurus specimen argues for an even earlier ontogenetic placement for specimens referred to Acrosaurus, the latter possibly pertaining to hatchlings.

Assessing the Adequacy of Morphological Models Using Posterior Predictive Simulations

Reconstructing the evolutionary history of different groups of organisms provides insight into how life originated and diversified on Earth. Phylogenetic trees are commonly used to estimate this evolutionary history. Within Bayesian phylogenetics a major step in estimating a tree is in choosing an appropriate model of character evolution. While the most common character data used is molecular sequence data, morphological data remains a vital source of information. The use of morphological characters allows for the incorporation fossil taxa, and despite advances in molecular sequencing, continues to play a significant role in neontology. Moreover, it is the main data source that allows us to unite extinct and extant taxa directly under the same generating process. We therefore require suitable models of morphological character evolution, the most common being the Mk Lewis model. While it is frequently used in both palaeobiology and neontology, it is not known whether the simple Mk substitution model, or any extensions to it, provide a sufficiently good description of the process of morphological evolution. In this study we investigate the impact of different morphological models on empirical tetrapod datasets. Specifically, we compare unpartitioned Mk models with those where characters are partitioned by the number of observed states, both with and without allowing for rate variation across sites and accounting for ascertainment bias. We show that the choice of substitution model has an impact on both topology and branch lengths, highlighting the importance of model choice. Through simulations, we validate the use of the model adequacy approach, posterior predictive simulations, for choosing an appropriate model. Additionally, we compare the performance of model adequacy with Bayesian model selection. We demonstrate how model selection approaches based on marginal likelihoods are not appropriate for choosing between models with partition schemes that vary in character state space (i.e., that vary in Q-matrix state size). Using posterior predictive simulations, we found that current variations of the Mk model are often performing adequately in capturing the evolutionary dynamics that generated our data. We do not find any preference for a particular model extension across multiple datasets, indicating that there is no "one size fits all" when it comes to morphological data and that careful consideration should be given to choosing models of discrete character evolution. By using suitable models of character evolution, we can increase our confidence in our phylogenetic estimates, which should in turn allow us to gain more accurate insights into the evolutionary history of both extinct and extant taxa.

Cladistic estimates of evolutionary rates focused on palaeontological datasets using TNT

We describe a protocol for estimating evolutionary rates from phylogenetic trees based on parsimony character optimization. The rate estimation is conducted through a TNT script and the results are analysed in a script for the software environment R. The TNT script allows analysing multiple optimal topologies, considering optimization ambiguity, and alternative time-calibrations or pre-calibrated trees. The R script summarizes estimated rates on a consensus tree and plots the variation of evolutionary rates through time, jointly with the phylogenetic diversity and a new metric (clade completeness index) that measures the distribution of missing data along the tree. We present results for simulated and empirical analyses, and evaluate the impact of missing data and alternative calibration methods in rate estimates. We found that while missing data can lower the nominal values of evolutionary rates, the overall pattern of rate variation through time remained robust. Empirical cases highlight different scenarios, such as datasets in which peaks of evolutionary rates can be coupled or decoupled from diversification dynamics (phylogenetic diversity) and cases in which missing data may influence the variation of estimated evolutionary rates. We conclude with recommendations for using this protocol and interpreting the results of parsimony-based rate estimates.

Early-middle Permian Mediterranean gorgonopsian suggests an equatorial origin of therapsids

Therapsids were a dominant component of middle-late Permian terrestrial ecosystems worldwide, eventually giving rise to mammals during the early Mesozoic. However, little is currently known about the time and place of origin of Therapsida. Here we describe a definitive therapsid from the lower-?middle Permian palaeotropics, a partial skeleton of a gorgonopsian from the island of Mallorca, western Mediterranean. This specimen represents, to our knowledge, the oldest gorgonopsian record worldwide, and possibly the oldest known therapsid. Using emerging relaxed clock models, we provide a quantitative timeline for the origin and early diversification of therapsids, indicating a long ghost lineage leading to the evolutionary radiation of all major therapsid clades within less than 10 Myr, in the aftermath of Olson's Extinction. Our findings place this unambiguous early therapsid in an ancient summer wet biome of equatorial Pangaea, thus suggesting that the group originated in tropical rather than temperate regions.

Bayesian tip-dated timeline for diversification and major biogeographic events in Muroidea (Rodentia), the largest mammalian radiation

BACKGROUND

Extinct organisms provide vital information about the time of origination and biogeography of extant groups. The development of phylogenetic methods to study evolutionary processes through time has revolutionized the field of evolutionary biology and led to an unprecedented expansion of our knowledge of the tree of life. Recent developments applying Bayesian approaches, using fossil taxa as tips to be included alongside their living relatives, have revitalized the use of morphological data in evolutionary tree inferences. Eumuroida rodents represent the largest group of mammals including more than a quarter of all extant mammals and have a rich fossil record spanning the last ~ 45 million years. Despite this wealth of data, our current understanding of the classification, major biogeographic patterns, and divergence times for this group comes from molecular phylogenies that use fossils only as a source of node calibrations. However, node calibrations impose several constraints on tree topology and must necessarily make a priori assumptions about the placement of fossil taxa without testing their placement in the tree.

RESULTS

We present the first morphological dataset with extensive fossil sampling for Muroidea. By applying Bayesian morphological clocks with tip dating and process-based biogeographic models, we provide a novel hypothesis for muroid relationships and revised divergence times for the clade that incorporates uncertainty in the placement of all fossil species. Even under strong violation of the clock model, we found strong congruence between results for divergence times, providing a robust timeline for muroid diversification. This new timeline was used for biogeographic analyses, which revealed a dynamic scenario mostly explained by dispersal events between and within the Palearctic and North African regions.

CONCLUSIONS

Our results provide important insights into the evolution of Muroidea rodents and clarify the evolutionary pathways of their main lineages. We exploited the advantage of tip dating Bayesian approaches in morphology-based datasets and provided a classification of the largest superfamily of mammals resulting from robust phylogenetic inference, inferring the biogeographical history, diversification, and divergence times of its major lineages.

Phylogenetic relationships of Neogene hamsters (Mammalia, Rodentia, Cricetinae) revealed under Bayesian inference and maximum parsimony

There is an ongoing debate about the internal systematics of today's group of hamsters (Cricetinae), following new insights that are gained based on molecular data. Regarding the closely related fossil cricetids, however, most studies deal with only a limited number of genera and statements about their possible relationships are rare. In this study, 41 fossil species from the Late Miocene to the Pliocene, belonging to seven extinct cricetine genera, Collimys, Rotundomys, Neocricetodon, Pseudocricetus, Cricetulodon, Apocricetus and Hattomys are analysed in a phylogenetic framework using traditional maximum parsimony and Bayesian inference approaches. Following thorough model testing, a relaxed-clock Bayesian inference analysis is performed under tip-dating to estimate divergence times simultaneously. Furthermore, so-called 'rogue' taxa are identified and excluded from the final trees to improve the informative value of the shown relationships. Based on these resulting trees, the fit of the topologies to the stratigraphy is assessed and the ancestral states of the characters are reconstructed under a parsimonious approach and stochastic character mapping. The overall topologies resulting from Bayesian and parsimonious approaches are largely congruent to each other and confirm the monophyly of most of the genera. Additionally, synapomorphies can be identified for each of these genera based on the ancestral state reconstructions. Only Cricetulodon turns out to be paraphyletic, while 'Cricetulodon' complicidens is a member of Neocricetodon. Lastly, this work makes a contribution to a debate that went on for decades, as the genus Kowalskia can be confirmed as junior synonym of Neocricetodon.

A new abelisaurid dinosaur from the end Cretaceous of Patagonia and evolutionary rates among the Ceratosauria

Gondwanan dinosaur faunae during the 20 Myr preceding the Cretaceous-Palaeogene (K/Pg) extinction included several lineages that were absent or poorly represented in Laurasian landmasses. Among these, the South American fossil record contains diverse abelisaurids, arguably the most successful groups of carnivorous dinosaurs from Gondwana in the Cretaceous, reaching their highest diversity towards the end of this period. Here we describe Koleken inakayali gen. et sp. n., a new abelisaurid from the La Colonia Formation (Maastrichtian, Upper Cretaceous) of Patagonia. Koleken inakayali is known from several skull bones, an almost complete dorsal series, complete sacrum, several caudal vertebrae, pelvic girdle and almost complete hind limbs. The new abelisaurid shows a unique set of features in the skull and several anatomical differences from Carnotaurus sastrei (the only other abelisaurid known from the La Colonia Formation). Koleken inakayali is retrieved as a brachyrostran abelisaurid, clustered with other South American abelisaurids from the latest Cretaceous (Campanian-Maastrichtian), such as Aucasaurus, Niebla and Carnotaurus. Leveraging our phylogeny estimates, we explore rates of morphological evolution across ceratosaurian lineages, finding them to be particularly high for elaphrosaurine noasaurids and around the base of Abelisauridae, before the Early Cretaceous radiation of the latter clade. The Noasauridae and their sister clade show contrasting patterns of morphological evolution, with noasaurids undergoing an early phase of accelerated evolution of the axial and hind limb skeleton in the Jurassic, and the abelisaurids exhibiting sustained high rates of cranial evolution during the Early Cretaceous. These results provide much needed context for the evolutionary dynamics of ceratosaurian theropods, contributing to broader understanding of macroevolutionary patterns across dinosaurs.

Toward transparent taxonomy: an interactive web-tool for evaluating competing taxonomic arrangements

Informative and consistent taxonomy above the species level is essential to communication about evolution, biodiversity and conservation, and yet the practice of taxonomy is considered opaque and subjective by non-taxonomist scientists and the public alike. While various proposals have tried to make the basis for the ranking and inclusiveness of taxa more transparent and objective, widespread adoption of these ideas has lagged. Here, we present TaxonomR, an interactive online decision-support tool to evaluate alternative taxonomic classifications. This tool implements an approach that quantifies the criteria commonly used in taxonomic treatments and allows the user to interactively manipulate weightings for different criteria to compare scores for taxonomic groupings under those weights. We use the butterfly taxon Argynnis to demonstrate how different weightings applied to common taxonomic criteria result in fundamentally different genus-level classifications that are predominantly used in different continents and geographic regions. These differences are objectively compared and quantified using TaxonomR to evaluate the kinds of criteria that have been emphasized in earlier classifications, and the nature of the support for current alternative taxonomic arrangements. The main role of TaxonomR is to make taxonomic decisions transparent via an explicit prioritization scheme. TaxonomR is not a prescriptive application. Rather, it aims to be a tool for facilitating our understanding of alternative taxonomic classifications that can, in turn, potentially support global harmony in biodiversity assessments through evidence-based discussion and community-wide resolution of historically entrenched taxonomic tensions.

Virtual endocasts of Clevosaurus brasiliensis and the tuatara: Rhynchocephalian neuroanatomy and the oldest endocranial record for Lepidosauria

Understanding the origins of the vertebrate brain is fundamental for uncovering evolutionary patterns in neuroanatomy. Regarding extinct species, the anatomy of the brain and other soft tissues housed in endocranial spaces can be approximated by casts of these cavities (endocasts). The neuroanatomical knowledge of Rhynchocephalia, a reptilian clade exceptionally diverse in the early Mesozoic, is restricted to the brain of its only living relative, Sphenodon punctatus, and unknown for fossil species. Here, we describe the endocast and the reptilian encephalization quotient (REQ) of the Triassic rhynchocephalian Clevosaurus brasiliensis and compare it with an ontogenetic series of S. punctatus. To better understand the informative potential of endocasts in Rhynchocephalia, we also examine the brain-endocast relationship in S. punctatus. We found that the brain occupies 30% of its cavity, but the latter recovers the general shape and length of the brain. The REQ of C. brasiliensis (0.27) is much lower than S. punctatus (0.84-1.16), with the tuatara being close to the mean for non-avian reptiles. The endocast of S. punctatus is dorsoventrally flexed and becomes more elongated throughout ontogeny. The endocast of C. brasiliensis is mostly unflexed and tubular, possibly representing a more plesiomorphic anatomy in relation to S. punctatus. Given the small size of C. brasiliensis, the main differences may result from allometric and heterochronic phenomena, consistent with suggestions that S. punctatus shows peramorphic anatomy compared to Mesozoic rhynchocephalians. Our results highlight a previously undocumented anatomical diversity among rhynchocephalians and provide a framework for future neuroanatomical comparisons among lepidosaurs.

A new sphenodontian (Diapsida: Lepidosauria) from the Upper Triassic (Norian) of Germany and its implications for the mode of sphenodontian evolution

The Arnstadt Formation of Saxony-Anhalt, Germany has yielded some of Germany's most substantial finds of Late Triassic tetrapods, including the sauropodomorph Plateosaurus and the stem-turtle Proganochelys quenstedti. Here, we describe an almost complete skull of a new sphenodontian taxon from this formation (Norian, 227-208 Ma), making it the oldest known articulated sphenodontian skull from Europe and one of the oldest in the world. The material is represented by the dermal skull roof and by the complete maxilla and temporal region, as well as parts of the palate, braincase, and lower jaw. A phylogenetic assessment recovers it as a basal sphenodontian closely related to Planocephalosaurus robinsonae and to Eusphenodontia, making it the earliest-diverging sphenodontian known with an articulated skull. Its cranial anatomy is generally similar to the well-known Diphydontosaurus avonis from the Rhaetian of England, showing that this successful phenotype was already established in the clade around 10 myr earlier than assumed. An analysis of evolutionary change rates recovers high rates of evolution in basal sphenodontians, with decreasing rates throughout the evolution of the group. However, contrary to previous studies, reversals in this trend were identified, indicating additional peaks of evolutionary change. These results improve our understanding of the early sphenodontian diversity in Europe, providing critical information on evolutionary trends throughout the history of the clade and sparking renewed interest in its evolution.

Evaluating the Accuracy of Methods for Detecting Correlated Rates of Molecular and Morphological Evolution

Determining the link between genomic and phenotypic change is a fundamental goal in evolutionary biology. Insights into this link can be gained by using a phylogenetic approach to test for correlations between rates of molecular and morphological evolution. However, there has been persistent uncertainty about the relationship between these rates, partly because conflicting results have been obtained using various methods that have not been examined in detail. We carried out a simulation study to evaluate the performance of 5 statistical methods for detecting correlated rates of evolution. Our simulations explored the evolution of molecular sequences and morphological characters under a range of conditions. Of the methods tested, Bayesian relaxed-clock estimation of branch rates was able to detect correlated rates of evolution correctly in the largest number of cases. This was followed by correlations of root-to-tip distances, Bayesian model selection, independent sister-pairs contrasts, and likelihood-based model selection. As expected, the power to detect correlated rates increased with the amount of data, both in terms of tree size and number of morphological characters. Likewise, greater among-lineage rate variation in the data led to improved performance of all 5 methods, particularly for Bayesian relaxed-clock analysis when the rate model was mismatched. We then applied these methods to a data set from flowering plants and did not find evidence of a correlation in evolutionary rates between genomic data and morphological characters. The results of our study have practical implications for phylogenetic analyses of combined molecular and morphological data sets, and highlight the conditions under which the links between genomic and phenotypic rates of evolution can be evaluated quantitatively.

A Morrison stem gekkotan reveals gecko evolution and Jurassic biogeography

Geckos are a speciose and globally distributed clade of Squamata (lizards, including snakes and amphisbaenians) that are characterized by a host of modifications for nocturnal, scansorial and insectivorous ecologies. They are among the oldest divergences in the lizard crown, so understanding the origin of geckoes (Gekkota) is essential to understanding the origin of Squamata, the most species-rich extant tetrapod clade. However, the poor fossil record of gekkotans has obscured the sequence and timing of the assembly of their distinctive morphology. Here, we describe the first North American stem gekkotan based on a three-dimensionally preserved skull from the Morrison Formation of western North America. Despite its Late Jurassic age, the new species already possesses several key characteristics of the gekkotan skull along with retained ancestral features. We show that this new stem gekkotan, and several previously named species of uncertain phylogenetic relationships, comprise a widespread clade of early crown lizards, substantiating faunal homogeneity in Laurasia during the Late Jurassic that extended across disparate ecological, body-size and physiological classes.

The affinities of the Late Triassic Cryptovaranoides and the age of crown squamates

Most living reptile diversity is concentrated in Squamata (lizards, including snakes), which have poorly known origins in space and time. Recently, †Cryptovaranoides microlanius from the Late Triassic of the United Kingdom was described as the oldest crown squamate. If true, this result would push back the origin of all major lizard clades by 30-65 Myr and suggest that divergence times for reptile clades estimated using genomic and morphological data are grossly inaccurate. Here, we use computed tomography scans and expanded phylogenetic datasets to re-evaluate the phylogenetic affinities of †Cryptovaranoides and other putative early squamates. We robustly reject the crown squamate affinities of †Cryptovaranoides, and instead resolve †Cryptovaranoides as a potential member of the bird and crocodylian total clade, Archosauromorpha. Bayesian total evidence dating supports a Jurassic origin of crown squamates, not Triassic as recently suggested. We highlight how features traditionally linked to lepidosaurs are in fact widespread across Triassic reptiles. Our study reaffirms the importance of critically choosing and constructing morphological datasets and appropriate taxon sampling to test the phylogenetic affinities of problematic fossils and calibrate the Tree of Life.

An early-diverging iguanodontian (Dinosauria: Rhabdodontomorpha) from the Late Cretaceous of North America

Intensifying macrovertebrate reconnaissance together with refined age-dating of mid-Cretaceous assemblages in recent decades is producing a more nuanced understanding of the impact of the Cretaceous Thermal Maximum on terrestrial ecosystems. Here we report discovery of a new early-diverging ornithopod, Iani smithi gen. et sp. nov., from the Cenomanian-age lower Mussentuchit Member, Cedar Mountain Formation of Utah, USA. The single known specimen of this species (NCSM 29373) includes a well-preserved, disarticulated skull, partial axial column, and portions of the appendicular skeleton. Apomorphic traits are concentrated on the frontal, squamosal, braincase, and premaxilla, including the presence of three premaxillary teeth. Phylogenetic analyses using parsimony and Bayesian inference posit Iani as a North American rhabdodontomorph based on the presence of enlarged, spatulate teeth bearing up to 12 secondary ridges, maxillary teeth lacking a primary ridge, a laterally depressed maxillary process of the jugal, and a posttemporal foramen restricted to the squamosal, among other features. Prior to this discovery, neornithischian paleobiodiversity in the Mussentuchit Member was based primarily on isolated teeth, with only the hadrosauroid Eolambia caroljonesa named from macrovertebrate remains. Documentation of a possible rhabdodontomorph in this assemblage, along with published reports of an as-of-yet undescribed thescelosaurid, and fragmentary remains of ankylosaurians and ceratopsians confirms a minimum of five, cohabiting neornithischian clades in earliest Late Cretaceous terrestrial ecosystems of North America. Due to poor preservation and exploration of Turonian-Santonian assemblages, the timing of rhabdodontomorph extirpation in the Western Interior Basin is, as of yet, unclear. However, Iani documents survival of all three major clades of Early Cretaceous neornithischians (Thescelosauridae, Rhabdodontomorpha, and Ankylopollexia) into the dawn of the Late Cretaceous of North America.

Syngnathoid Evolutionary History and the Conundrum of Fossil Misplacement

Seahorses, pipefishes, trumpetfishes, shrimpfishes, and allies are a speciose, globally distributed clade of fishes that have evolved a large number of unusual body plans. The clade that includes all these forms, Syngnathoidei, has become a model for the study of life history evolution, population biology, and biogeography. Yet, the timeline of syngnathoid evolution has remained highly contentious. This debate is largely attributable to the nature of the syngnathoid fossil record, which is both poorly described and patchy for several major lineages. Although fossil syngnathoids have been used to calibrate molecular phylogenies, the interrelationships of extinct species and their affinities to major living syngnathoid clades have scarcely been quantitatively tested. Here, I use an expanded morphological dataset to reconstruct the evolutionary relationships and clade ages of fossil and extant syngnathoids. Phylogenies generated using different analytical methodologies are largely congruent with molecular phylogenetic trees of Syngnathoidei but consistently find novel placements for several key taxa used as fossil calibrators in phylogenomic studies. Tip-dating of the syngnathoid phylogeny finds a timeline for their evolution that differs slightly from the one inferred using molecular trees but is generally congruent with a post-Cretaceous diversification event. These results emphasize the importance of quantitatively testing the relationships of fossil species, particularly when they are critical to assessing divergence times.

Total evidence phylogeny of platyrrhine primates and a comparison of undated and tip-dating approaches

There have been multiple published phylogenetic analyses of platyrrhine primates (New World monkeys) using both morphological and molecular data, but relatively few that have integrated both types of data into a total evidence approach. Here, we present phylogenetic analyses of recent and fossil platyrrhines, based on a total evidence data set of 418 morphological characters and 10.2 kilobases of DNA sequence data from 17 nuclear genes taken from previous studies, using undated and tip-dating approaches in a Bayesian framework. We compare the results of these analyses with molecular scaffold analyses using maximum parsimony and Bayesian approaches, and we use a formal information theoretic approach to identify unstable taxa. After a posteriori pruning of unstable taxa, the undated and tip-dating topologies appear congruent with recent molecular analyses and support largely similar relationships, with strong support for Stirtonia as a stem alouattine, Neosaimiri as a stem saimirine, Cebupithecia as a stem pitheciine, and Lagonimico as a stem callitrichid. Both analyses find three Greater Antillean subfossil platyrrhines (Xenothrix, Antillothrix, and Paralouatta) to form a clade that is related to Callicebus, congruent with a single dispersal event by the ancestor of this clade to the Greater Antilles. They also suggest that the fossil Proteropithecia may not be closely related to pitheciines, and that all known platyrrhines older than the Middle Miocene are stem taxa. Notably, the undated analysis found the Early Miocene Panamacebus (currently recognized as the oldest known cebid) to be unstable, and the tip-dating analysis placed it outside crown Platyrrhini. Our tip-dating analysis supports a late Oligocene or earliest Miocene (20.8-27.0 Ma) age for crown Platyrrhini, congruent with recent molecular clock analyses.

A Triassic crown squamate

Mammals, birds, and squamates (lizards, snakes, and relatives) are key living vertebrates, and thus understanding their evolution underpins important questions in biodiversity science. Whereas the origins of mammals and birds are relatively well understood, the roots of squamates have been obscure. Here, we report a modern-type lizard from the Late Triassic of England [202 million years (Ma)], comprising a partial skeleton, skull, and mandibles. It displays at least 15 unique squamate traits and further shares unidentatan and anguimorph apomorphies. The new discovery fixes the origin of crown Squamata as much older than had been thought, and the revised dating shows substantial diversification of modern-type squamates following the Carnian Pluvial Episode, 232 Ma ago.

Evolutionary origins of the prolonged extant squamate radiation

Squamata is the most diverse clade of terrestrial vertebrates. Although the origin of pan-squamates lies in the Triassic, the oldest undisputed members of extant clades known from nearly complete, uncrushed material come from the Cretaceous. Here, we describe three-dimensionally preserved partial skulls of two new crown lizards from the Late Jurassic of North America. Both species are placed at the base of the skink, girdled, and night lizard clade Pan-Scincoidea, which consistently occupies a position deep inside the squamate crown in both morphological and molecular phylogenies. The new lizards show that several features uniting pan-scincoids with another major lizard clade, the pan-lacertoids, in trees using morphology were convergently acquired as predicted by molecular analyses. Further, the palate of one new lizard bears a handful of ancestral saurian characteristics lost in nearly all extant squamates, revealing an underappreciated degree of complex morphological evolution in the early squamate crown. We find strong evidence for close relationships between the two new species and Cretaceous taxa from Eurasia. Together, these results suggest that early crown squamates had a wide geographic distribution and experienced complicated morphological evolution even while the Rhynchocephalia, now solely represented by the tuatara, was the dominant clade of lepidosaurs.

Successive climate crises in the deep past drove the early evolution and radiation of reptiles

Climate change-induced mass extinctions provide unique opportunities to explore the impacts of global environmental disturbances on organismal evolution. However, their influence on terrestrial ecosystems remains poorly understood. Here, we provide a new time tree for the early evolution of reptiles and their closest relatives to reconstruct how the Permian-Triassic climatic crises shaped their long-term evolutionary trajectory. By combining rates of phenotypic evolution, mode of selection, body size, and global temperature data, we reveal an intimate association between reptile evolutionary dynamics and climate change in the deep past. We show that the origin and phenotypic radiation of reptiles was not solely driven by ecological opportunity following the end-Permian extinction as previously thought but also the result of multiple adaptive responses to climatic shifts spanning 57 million years.

Integrative Phylogenetics: Tools for Palaeontologists to Explore the Tree of Life

The modern era of analytical and quantitative palaeobiology has only just begun, integrating methods such as morphological and molecular phylogenetics and divergence time estimation, as well as phenotypic and molecular rates of evolution. Calibrating the tree of life to geological time is at the nexus of many disparate disciplines, from palaeontology to molecular systematics and from geochronology to comparative genomics. Creating an evolutionary time scale of the major events that shaped biodiversity is key to all of these fields and draws from each of them. Different methodological approaches and data employed in various disciplines have traditionally made collaborative research efforts difficult among these disciplines. However, the development of new methods is bridging the historical gap between fields, providing a holistic perspective on organismal evolutionary history, integrating all of the available evidence from living and fossil species. Because phylogenies with only extant taxa do not contain enough information to either calibrate the tree of life or fully infer macroevolutionary dynamics, phylogenies should preferably include both extant and extinct taxa, which can only be achieved through the inclusion of phenotypic data. This integrative phylogenetic approach provides ample and novel opportunities for evolutionary biologists to benefit from palaeontological data to help establish an evolutionary time scale and to test core macroevolutionary hypotheses about the drivers of biological diversification across various dimensions of organisms.

Species delimitation and coexistence in an ancient, depauperate vertebrate clade

Background

A major challenge to understanding how biodiversity has changed over time comes from depauperons, which are long-lived lineages with presently low species diversity. The most famous of these are the coelacanths. This clade of lobe-finned fishes occupies a pivotal position on the vertebrate tree between other fishes and tetrapods. Yet only two extant species and fewer than 100 extinct forms are known from the coelacanth fossil record, which spans over 400 million years of time. Although there is evidence for the existence of additional genetically isolated extant populations, a poor understanding of morphological disparity in this clade has made quantifying coelacanth species richness difficult.

Results

Here, we quantify variation in a sample of skulls and skeletons of the Triassic eastern North American coelacanth †Diplurus that represents the largest assemblage of coelacanth individuals known. Based on the results of these quantitative comparisons, we identify a diminutive new species and show that multiple lacustrine ecosystems in the Triassic rift lakes of the Atlantic coastline harbored at least three species of coelacanths spanning two orders of magnitude in size.

Conclusions

Conceptions about the distribution of species diversity on the tree of life may be fundamentally misguided when extant diversity is used to gauge signals of extinct diversity.

Our results demonstrate how specimen-based assessments can be used to illuminate hidden biodiversity and show the utility of the fossil record for answering questions about the hidden richness of currently species-poor lineages.

Early cephalopod evolution clarified through Bayesian phylogenetic inference

BACKGROUND

Despite the excellent fossil record of cephalopods, their early evolution is poorly understood. Different, partly incompatible phylogenetic hypotheses have been proposed in the past, which reflected individual author's opinions on the importance of certain characters but were not based on thorough cladistic analyses. At the same time, methods of phylogenetic inference have undergone substantial improvements. For fossil datasets, which typically only include morphological data, Bayesian inference and in particular the introduction of the fossilized birth-death model have opened new possibilities. Nevertheless, many tree topologies recovered from these new methods reflect large uncertainties, which have led to discussions on how to best summarize the information contained in the posterior set of trees.

RESULTS

We present a large, newly compiled morphological character matrix of Cambrian and Ordovician cephalopods to conduct a comprehensive phylogenetic analysis and resolve existing controversies. Our results recover three major monophyletic groups, which correspond to the previously recognized Endoceratoidea, Multiceratoidea, and Orthoceratoidea, though comprising slightly different taxa. In addition, many Cambrian and Early Ordovician representatives of the Ellesmerocerida and Plectronocerida were recovered near the root. The Ellesmerocerida is para- and polyphyletic, with some of its members recovered among the Multiceratoidea and early Endoceratoidea. These relationships are robust against modifications of the dataset. While our trees initially seem to reflect large uncertainties, these are mainly a consequence of the way clade support is measured. We show that clade posterior probabilities and tree similarity metrics often underestimate congruence between trees, especially if wildcard taxa are involved.

CONCLUSIONS

Our results provide important insights into the earliest evolution of cephalopods and clarify evolutionary pathways. We provide a classification scheme that is based on a robust phylogenetic analysis. Moreover, we provide some general insights on the application of Bayesian phylogenetic inference on morphological datasets. We support earlier findings that quartet similarity metrics should be preferred over the Robinson-Foulds distance when higher-level phylogenetic relationships are of interest and propose that using a posteriori pruned maximum clade credibility trees help in assessing support for phylogenetic relationships among a set of relevant taxa, because they provide clade support values that better reflect the phylogenetic signal.

An exceptionally preserved Sphenodon-like sphenodontian reveals deep time conservation of the tuatara skeleton and ontogeny

Sphenodontian reptiles are an extremely old evolutionary lineage forming the closest relatives to squamates (lizards and snakes) and were globally distributed and more diverse than squamates during the first half of their evolutionary history. However, the majority of their fossils are highly fragmentary, especially within sphenodontines-the group including its single surviving species, Sphenodon punctatus (the tuatara of New Zealand)-thus severely hampering our understanding on the origins of the tuatara. Here, we present a new sphenodontian species from the Early Jurassic of North America (Arizona, USA) represented by a nearly complete articulated skeleton and dozens of upper and lower jaws forming the most complete ontogenetic series in the sphenodontian fossil record. CT-scanning provides plentitude of data that unambiguously place this new taxon as one of the earliest evolving and oldest known sphenodontines. Comparisons with Sphenodon reveal that fundamental patterns of mandibular ontogeny and skeletal architecture in Sphenodon may have originated at least ~190Mya. In combination with recent findings, our results suggest strong morphological stability and an ancient origin of the modern tuatara morphotype.

Reappraisal and new material of the holotype of Draconyx loureiroi (Ornithischia: Iguanodontia) provide insights on the tempo and modo of evolution of thumb-spiked dinosaurs

The Upper Jurassic Lourinhã Formation is well known for its rich assemblage of fossil vertebrates. In this formation, ornithopod dinosaurs are represented by two iguanodontian species, Eousdryosaurus nanohallucis and Draconyx loureiroi. We recently became aware of unreported material belonging to the holotype of Draconyx loureiroi, consisting of partially articulated manual elements. We here re-describe the holotype specimen ML 357, including the newly discovered material. The specimen was subjected to CT-scanning and its surface data used to assess anatomical characters. Linear measurements of metatarsal III were used to estimate the body length of the specimen. The Draconyx loureiroi holotype was included in two datasets and analysed with maximum parsimony and Bayesian inference approaches to estimate evolutionary rates among Iguanodontia. We present evidence that Draconyx loureiroi is a valid taxon nested in Styracosterna and is clearly diagnosable by a unique combination of characters. Both maximum parsimony and Bayesian inference indicate high evolutionary rates across the Jurassic/Cretaceous transition for the base of Iguanodontia. Length estimation suggests that Draconyx loureiroi was a relatively small, bipedal and possibly cursorial animal. Given its basal phyletic position, we interpret this bauplan was the ancestral condition for Styracosterna, that only later in the Cretaceous evolved into giant quadrupedal forms.

A Triassic stem lepidosaur illuminates the origin of lizard-like reptiles

The early evolution of diapsid reptiles is marked by a deep contrast between our knowledge of the origin and early evolution of archosauromorphs (crocodiles, avian and non-avian dinosaurs) to that of lepidosauromorphs (squamates (lizards, snakes) and sphenodontians (tuataras)). Whereas the former include hundreds of fossil species across various lineages during the Triassic period¹, the latter are represented by an extremely patchy early fossil record comprising only a handful of fragmentary fossils, most of which have uncertain phylogenetic affinities and are confined to Europe^1-3. Here we report the discovery of a three-dimensionally preserved reptile skull, assigned as Taytalura alcoberi gen. et sp. nov., from the Late Triassic epoch of Argentina that is robustly inferred phylogenetically as the earliest evolving lepidosauromorph, using various data types and optimality criteria. Micro-computed tomography scans of this skull reveal details about the origin of the lepidosaurian skull from early diapsids, suggesting that several traits traditionally associated with sphenodontians in fact originated much earlier in lepidosauromorph evolution. Taytalura suggests that the strongly evolutionarily conserved skull architecture of sphenodontians represents the plesiomorphic condition for all lepidosaurs, that stem and crown lepidosaurs were contemporaries for at least ten million years during the Triassic, and that early lepidosauromorphs had a much broader geographical distribution than has previously been thought.

Sustained high rates of morphological evolution during the rise of tetrapods

The fish-to-tetrapod transition is one of the most iconic events in vertebrate evolution, yet fundamental questions regarding the dynamics of this transition remain unresolved. Here, we use advances in Bayesian morphological clock modelling to reveal the evolutionary dynamics of early tetrapodomorphs (tetrapods and their closest fish relatives). We show that combining osteological and ichnological calibration data results in major shifts on the time of origin of all major groups of tetrapodomorphs (up to 25 million years) and that low rates of net diversification, not fossilization, explain long ghost lineages in the early tetrapodomorph fossil record. Further, our findings reveal extremely low rates of morphological change for most early tetrapodomorphs, indicating widespread stabilizing selection upon their 'fish' morphotype. This pattern was broken only by elpistostegalians (including early tetrapods), which underwent sustained high rates of morphological evolution for ~30 Myr during the deployment of the tetrapod body plan.

Bayesian Morphological Clock versus Parsimony: An insight into the relationships and dispersal events of postvacuum Cricetidae (Rodentia, Mammalia)

Establishing an evolutionary timeline is fundamental for tackling a great variety of topics in evolutionary biology, including the reconstruction of patterns of historical biogeography, coevolution and diversification. However, the tree of life is pruned by extinction and molecular data cannot be gathered for extinct lineages. Until recently methodological challenges have prevented the application of tip-dating Bayesian approaches in morphology-based fossil-only datasets. Herein, we present a morphological dataset for a group of cricetid rodents to which we apply an array of methods fairly new in palaeontology that can be used by palaeontologists for the analysis of entirely extinct clades. We compare the tree topologies obtained by traditional parsimony, time-calibrated and non-calibrated Bayesian inference phylogenetic approaches and calculate stratigraphic congruence indices for each. Bayesian tip-dated clock methods outperform parsimony in the case of our dataset, which includes highly homoplastic morphological characters. Regardless, all three topologies support the monophyly of Megacricetodontinae, Democricetodontinae and Cricetodontinae. Dispersal and speciation events inferred through Bayesian Binary Markov chain Monte Carlo and biodiversity analyses provide evidence for a correlation between biogeographic events, climatic changes and diversification in cricetids.

Empirical and Methodological Challenges to the Model-Based Inference of Diversification Rates in Extinct Clades

Changes in speciation and extinction rates are key to the dynamics of clade diversification, but attempts to infer them from phylogenies of extant species face challenges. Methods capable of synthesizing information from extant and fossil species have yielded novel insights into diversification rate variation through time, but little is known about their behavior when analyzing entirely extinct clades. Here, we use empirical and simulated data to assess how two popular methods, PyRate and Fossil BAMM, perform in this setting. We inferred the first tip-dated trees for ornithischian dinosaurs, and combined them with fossil occurrence data to test whether the clade underwent an end-Cretaceous decline. We then simulated phylogenies and fossil records under empirical constraints to determine whether macroevolutionary and preservation rates can be teased apart under paleobiologically realistic conditions. We obtained discordant inferences about ornithischian macroevolution including a long-term speciation rate decline (BAMM), mostly flat rates with a steep diversification drop (PyRate) or without one (BAMM), and episodes of implausibly accelerated speciation and extinction (PyRate). Simulations revealed little to no conflation between speciation and preservation, but yielded spuriously correlated speciation and extinction estimates while time-smearing tree-wide shifts (BAMM) or overestimating their number (PyRate). Our results indicate that the small phylogenetic datasets available to vertebrate paleontologists and the assumptions made by current model-based methods combine to yield potentially unreliable inferences about the diversification of extinct clades. We provide guidelines for interpreting the results of the existing approaches in light of their limitations, and suggest how the latter may be mitigated.

Sphenofontis velserae gen. et sp. nov., a new rhynchocephalian from the Late Jurassic of Brunn (Solnhofen Archipelago, southern Germany)

The Solnhofen Archipelago is well known for its fossil vertebrates of Late Jurassic age, among which figure numerous rhynchocephalian specimens, representing at least six and up to nine genera. A new taxon, named Sphenofontis velserae gen. et sp. nov., increases rhynchocephalian diversity in the Solnhofen Archipelago and is herein described based on a single, well-preserved specimen originating from the Late Kimmeridgian of the Brunn quarry, near Regensburg. The exquisite preservation of the holotype allowed a detailed description of the animal, revealing a skeletal morphology that includes both plesiomorphic and derived features within rhynchocephalians. Sphenofontis is herein referred to Neosphenodontia and tentatively to sphenodontine sphenodontids. It notably differs from all other rhynchocephalians known from the Jurassic of Europe, showing instead closer resemblance with the Middle Jurassic Cynosphenodon from Mexico and especially the extant Sphenodon. This is evidence for a wide distribution of taxa related to the extant tuatara early in the Mesozoic, and also for the presence of less-specialized rhynchocephalians coexisting with more derived forms during the earliest time in the history of the Solnhofen Archipelago.