Bayesian Total-Evidence Dating Revisits Sloth Phylogeny and Biogeography: A Cautionary Tale on Morphological Clock Analyses

Combining morphological and molecular characters through Bayesian total-evidence dating allows inferring the phylogenetic and timescale framework of both extant and fossil taxa, while accounting for the stochasticity and incompleteness of the fossil record. Such an integrative approach is particularly needed when dealing with clades such as sloths (Mammalia: Folivora), for which developmental and biomechanical studies have shown high levels of morphological convergence whereas molecular data can only account for a limited percentage of their total species richness. Here, we propose an alternative hypothesis of sloth evolution that emphasizes the pervasiveness of morphological convergence and the importance of considering the fossil record and an adequate taxon sampling in both phylogenetic and biogeographic inferences. Regardless of different clock models and morphological datasets, the extant sloth Bradypus is consistently recovered as a megatherioid, and Choloepus as a mylodontoid, in agreement with molecular-only analyses. The recently extinct Caribbean sloths (Megalocnoidea) are found to be a monophyletic sister-clade of Megatherioidea, in contrast to previous phylogenetic hypotheses. Our results contradict previous morphological analyses and further support the polyphyly of "Megalonychidae," whose members were found in five different clades. Regardless of taxon sampling and clock models, the Caribbean colonization of sloths is compatible with the exhumation of islands along Aves Ridge and its geological time frame. Overall, our total-evidence analysis illustrates the difficulty of positioning highly incomplete fossils, although a robust phylogenetic framework was recovered by an a posteriori removal of taxa with high percentages of missing characters. Elimination of these taxa improved topological resolution by reducing polytomies and increasing node support. However, it introduced a systematic and geographic bias because most of these incomplete specimens are from northern South America. This is evident in biogeographic reconstructions, which suggest Patagonia as the area of origin of many clades when taxa are underrepresented, but Amazonia and/or Central and Southern Andes when all taxa are included. More generally, our analyses demonstrate the instability of topology and divergence time estimates when using different morphological datasets and clock models and thus caution against making macroevolutionary inferences when node support is weak or when uncertainties in the fossil record are not considered.

A 50-million-year-old, three-dimensionally preserved bat skull supports an early origin for modern echolocation

Bats are among the most recognizable, numerous, and widespread of all mammals. But much of their fossil record is missing, and bat origins remain poorly understood, as do the relationships of early to modern bats. Here, we describe a new early Eocene bat that helps bridge the gap between archaic stem bats and the hyperdiverse modern bat radiation of more than 1,460 living species. Recovered from ∼50 million-year-old cave sediments in the Quercy Phosphorites of southwestern France, Vielasia sigei's remains include a near-complete, three-dimensionally preserved skull-the oldest uncrushed bat cranium yet found. Phylogenetic analyses of a 2,665 craniodental character matrix, with and without 36.8 kb of DNA sequence data, place Vielasia outside modern bats, with total evidence tip-dating placing it sister to the crown clade. Vielasia retains the archaic dentition and skeletal features typical of early Eocene bats, but its inner ear shows specializations found in modern echolocating bats. These features, which include a petrosal only loosely attached to the basicranium, an expanded cochlea representing ∼25% basicranial width, and a long basilar membrane, collectively suggest that the kind of laryngeal echolocation used by most modern bats predates the crown radiation. At least 23 individuals of V. sigei are preserved together in a limestone cave deposit, indicating that cave roosting behavior had evolved in bats by the end of the early Eocene; this period saw the beginning of significant global climate cooling that may have been an evolutionary driver for bats to first congregate in caves.

Anatomy and phylogeny of a new small macraucheniid (Mammalia: Litopterna) from the Bahía Inglesa Formation (late Miocene), Atacama Region, Northern Chile

We describe a new macraucheniine macraucheniid, Micrauchenia saladensis gen. et sp. nov., from the late Miocene (Huayquerian SALMA). This is the first litoptern from Bahía Inglesa Formation, Chile. The specimen includes a partial mandible, cervical and thoracic vertebrae fragments, and portions of the forelimbs (a scapula fragment, an ulna-radius fragment, seven carpals, three metapodials, two proximal phalanges and four intermediate phalanges). The postcranial anatomy of Micrauchenia saladensis is consistent with terrestrial and cursorial locomotion, which suggests an allochthonous position of this specimen within the marine Bahía Inglesa Formation. The fusion of the ulna and radius and the presence of a radial aliform expansion align Micrauchenia with other macraucheniines, with which it shares these features. We interpret the fusion of the ulna and radius as a cursorial specialization and the aliform expansion as an adaptation for strong flexion movements and to resist higher transverse stresses during locomotion. In addition, Micrauchenia saladensis is the smallest member of the subfamily Macraucheniinae. To test the systematics and phylogenetics of this specimen, we expanded previous morphological matrices of macraucheniids by adding one dental and eight postcranial characters and scoring Micrauchenia saladensis. We performed maximum parsimony and Bayesian phylogenetic analyses, the latter applied for the first time to macraucheniid phylogeny. Our analyses confirm Micrauchenia saladensis as a member of the subfamily Macraucheniinae, although with uncertain affinities within this subfamily.

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.

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.

Haplocauda, a New Genus of Fireflies Endemic to the Amazon Rainforest (Coleoptera: Lampyridae)

Simple Summary

The Amazon Rainforest is still a frontier in the study of insect biodiversity, housing several species yet to be described and studied. Yet, its continuous deforestation urges scientists to study insect diversity there before it is too late. Here, we identified and described new species and a new genus of fireflies from the Amazon Rainforest. These new species are unique among other fireflies due to their unique abdominal morphology, supposed here to involve a clamping mechanism used during mating.

Abstract

Most firefly genera have poorly defined taxonomic boundaries, especially in the Neotropics, where they are more diverse and more difficult to identify. Recent advances that shed light on the diversity of fireflies in South America have focused mainly on Atlantic Rainforest taxa, whereas lampyrids in other biomes remained largely unstudied. We found three new firefly species endemic to the Amazon basin that share unique traits of the male abdomen where sternum VIII and the pygidium are modified and likely work as a copulation clamp. Here we test and confirm the hypothesis that these three species form a monophyletic lineage and propose Haplocaudagen. nov. to accommodate the three new species. Both maximum parsimony and probabilistic (Bayesian and maximum likelihood) phylogenetic analyses confirmed Haplocaudagen. nov. monophyly, and consistently recovered it as the sister group to Scissicauda, fireflies endemic to the Atlantic Rainforest that also feature a copulation clamp on abdominal segment VIII, although with a different configuration. We provide illustrations, diagnostic descriptions, and keys to species based on males and females. The three new species were sampled from different regions, and are likely allopatric, a common pattern among Amazonian taxa.

Morphology and Bayesian tip-dating recover deep Cretaceous-age divergences among major chrysidid lineages (Hymenoptera: Chrysididae)

We integrated phylogenetic, biogeographic and palaeontological data to reconstruct the evolutionary history of the cuckoo wasps. We propose a phylogenetic hypothesis based on a comprehensive morphological study resulting in 300 characters coded for both living and extinct species. Phylogenetic relationships and divergence time estimation were simultaneously inferred in a Bayesian tip-dating framework, applying a relaxed morphological clock. Results unequivocally indicate Chrysididae to be monophyletic, as well as all traditionally recognized subfamilies and tribes. Within the Chrysidinae, Elampini was placed as the sister-group of the other three chrysidine tribes, with Parnopini as sister to the clade including Allocoeliini and Chrysidini. Dating analysis indicates that the major lineages started to differentiate around 130 Mya during the Early Cretaceous. The clades recognized as subfamilies started differentiating during the Palaeogene and the Neogene. Our results reveal an intricate process on the geographic evolution of chrysidid wasps and dispute previous ideas that Cretaceous-old splits in their early history could be associated with vicariant events related to the breakup between Africa and South America. The present-day southern disjunctions of some groups are interpreted as the outcome of more recent dispersals and extinctions of representatives from Nearctic and Palaearctic faunas during the Neogene, when northern continents became significantly colder.

Investigating Morphological Complexes Using Informational Dissonance and Bayes Factors: A Case Study in Corbiculate Bees

It is widely recognized that different regions of a genome often have different evolutionary histories and that ignoring this variation when estimating phylogenies can be misleading. However, the extent to which this is also true for morphological data is still largely unknown. Discordance among morphological traits might plausibly arise due to either variable convergent selection pressures or else phenomena such as hemiplasy. Here, we investigate patterns of discordance among 282 morphological characters, which we scored for 50 bee species particularly targeting corbiculate bees, a group that includes the well-known eusocial honeybees and bumblebees. As a starting point for selecting the most meaningful partitions in the data, we grouped characters as morphological modules, highly integrated trait complexes that as a result of developmental constraints or coordinated selection we expect to share an evolutionary history and trajectory. In order to assess conflict and coherence across and within these morphological modules, we used recently developed approaches for computing Bayesian phylogenetic information allied with model comparisons using Bayes factors. We found that despite considerable conflict among morphological complexes, accounting for among-character and among-partition rate variation with individual gamma distributions, rate multipliers, and linked branch lengths can lead to coherent phylogenetic inference using morphological data. We suggest that evaluating information content and dissonance among partitions is a useful step in estimating phylogenies from morphological data, just as it is with molecular data. Furthermore, we argue that adopting emerging approaches for investigating dissonance in genomic datasets may provide new insights into the integration and evolution of anatomical complexes. [Apidae; entropy; morphological modules; phenotypic integration; phylogenetic information.].

A monotreme-like auditory apparatus in a Middle Jurassic haramiyidan

Among extant vertebrates, mammals are distinguished by having a chain of three auditory ossicles (the malleus, incus and stapes) that transduce sound waves and promote an increased range of audible-especially high-frequencies¹. By contrast, the homologous bones in early fossil mammals and relatives also functioned in chewing through their bony attachments to the lower jaw². Recent discoveries of well-preserved Mesozoic mammals have provided glimpses into the transition from the dual (masticatory and auditory) to the single auditory function for the ossicles, which is now widely accepted to have occurred at least three times in mammal evolution^3-6. Here we report a skull and postcranium that we refer to the haramiyidan Vilevolodon diplomylos (dating to the Middle Jurassic epoch (160 million years ago)) and that shows excellent preservation of the malleus, incus and ectotympanic (which supports the tympanic membrane). After comparing this fossil with other Mesozoic and extant mammals, we propose that the overlapping incudomallear articulation found in this and other Mesozoic fossils, in extant monotremes and in early ontogeny in extant marsupials and placentals is a morphology that evolved in several groups of mammals in the transition from the dual to the single function for the ossicles.

A new family of diprotodontian marsupials from the latest Oligocene of Australia and the evolution of wombats, koalas, and their relatives (Vombatiformes)

We describe the partial cranium and skeleton of a new diprotodontian marsupial from the late Oligocene (~26–25 Ma) Namba Formation of South Australia. This is one of the oldest Australian marsupial fossils known from an associated skeleton and it reveals previously unsuspected morphological diversity within Vombatiformes, the clade that includes wombats (Vombatidae), koalas (Phascolarctidae) and several extinct families. Several aspects of the skull and teeth of the new taxon, which we refer to a new family, are intermediate between members of the fossil family Wynyardiidae and wombats. Its postcranial skeleton exhibits features associated with scratch-digging, but it is unlikely to have been a true burrower. Body mass estimates based on postcranial dimensions range between 143 and 171 kg, suggesting that it was ~5 times larger than living wombats. Phylogenetic analysis based on 79 craniodental and 20 postcranial characters places the new taxon as sister to vombatids, with which it forms the superfamily Vombatoidea as defined here. It suggests that the highly derived vombatids evolved from wynyardiid-like ancestors, and that scratch-digging adaptations evolved in vombatoids prior to the appearance of the ever-growing (hypselodont) molars that are a characteristic feature of all post-Miocene vombatids. Ancestral state reconstructions on our preferred phylogeny suggest that bunolophodont molars are plesiomorphic for vombatiforms, with full lophodonty (characteristic of diprotodontoids) evolving from a selenodont morphology that was retained by phascolarctids and ilariids, and wynyardiids and vombatoids retaining an intermediate selenolophodont condition. There appear to have been at least six independent acquisitions of very large (>100 kg) body size within Vombatiformes, several having already occurred by the late Oligocene.

An extended model for phylogenetic maximum likelihood based on discrete morphological characters

Maximum likelihood is a common method of estimating a phylogenetic tree based on a set of genetic data. However, models of evolution for certain types of genetic data are highly flawed in their specification, and this misspecification can have an adverse impact on phylogenetic inference. Our attention here is focused on extending an existing class of models for estimating phylogenetic trees from discrete morphological characters. The main advance of this work is a model that allows unequal equilibrium frequencies in the estimation of phylogenetic trees from discrete morphological character data using likelihood methods. Possible extensions of the proposed model will also be discussed.

Homoplasy-Based Partitioning Outperforms Alternatives in Bayesian Analysis of Discrete Morphological Data

Bayesian analysis of morphological data is becoming increasingly popular mainly (but not only) because it allows for time-calibrated phylogenetic inference using relaxed morphological clocks and tip dating whenever fossils are available. As with molecular data, recent studies have shown that modeling among-character rate variation (ACRV) in morphological matrices greatly improves phylogenetic inference. In a likelihood framework this may be accomplished, for instance, by employing a hidden Markov model to assign characters to rate categories drawn from a (discretized) $\Gamma$ distribution and/or by partitioning data sets according to rate heterogeneity and estimating per-partition branch lengths, conditioned on a single topology. While the first approach is available in many phylogenetic analysis software, there is still no clear consensus on how to partition data, except perhaps in the simplest cases (e.g., "by codon" partitioning of coding sequences). Additionally, there is a trade-off between improvement in likelihood scores and the number of free parameters in the analysis, which rises quickly with the number of partitions. This trade-off may be dealt with by employing statistics that penalize overfitting of complex models, such as Akaike or Bayesian information criteria, or the more recently introduced stepping-stone method for marginal likelihood approximation. We applied the latter to three distinct matrices of discrete morphological data and demonstrated that sorting characters by homoplasy scores (obtained from implied weighting parsimony analysis) outperformed other partitioning strategies (anatomically-based and PartitionFinder2). The method was in fact so efficient in segregating characters by rates of evolution that no within-partition ACRV modeling was necessary, while among-partition rate variation was adequately accommodated by rate multipliers. We conclude that partitioning by homoplasy is a powerful and easy-to-implement strategy to address ACRV in complex data sets. We provide some guidelines focusing on morphological matrices, although this approach may be also applicable to molecular data sets.

A Systematist's Guide to Estimating Bayesian Phylogenies From Morphological Data

Phylogenetic trees are crucial to many aspects of taxonomic and comparative biology. Many researchers have adopted Bayesian methods to estimate their phylogenetic trees. In this family of methods, a model of morphological evolution is assumed to have generated the data observed by the researcher. These models make a variety of assumptions about the evolution of morphological characters, and these assumptions are translated into mathematics as parameters. The incorporation of prior distributions further allows researchers to quantify their prior beliefs about the value any one parameter can take. How to translate biological knowledge into mathematical language is difficult, and can be confusing to many biologists. This review aims to help systematics researchers understand the biological meaning of common models and assumptions. Using examples from the insect fossil record, I will demonstrate empirically what assumptions mean in concrete terms, and discuss how researchers can use and understand Bayesian methods for phylogenetic estimation.

Morphological Data Sets Fit a Common Mechanism Much More Poorly than DNA Sequences and Call Into Question the Mkv Model

The Mkv evolutionary model, based on minor modifications to models of molecular evolution, is being increasingly used to infer phylogenies from discrete morphological data, often producing different results from parsimony. The critical difference between Mkv and parsimony is the assumption of a "common mechanism" in the Mkv model, with branch lengths determining that probability of change for all characters increases or decreases at the same tree branches by the same exponential factor. We evaluate whether the assumption of a common mechanism applies to morphology, by testing the implicit prediction that branch lengths calculated from different subsets of characters will be significantly correlated. Our analysis shows that DNA (38 data sets tested) is often compatible with a common mechanism, but morphology (86 data sets tested) generally is not, showing very disparate branch lengths for different character partitions. The low levels of branch length correlation demonstrated for morphology (fitting models without a common mechanism) suggest that the Mkv model is too unrealistic and inadequate for the analysis of most morphological data sets. [Bayesian analysis; Mkv model; morphological data; phylogenetics.].

Exceptional preservation of mid-Cretaceous marine arthropods and the evolution of novel forms via heterochrony

Evolutionary origins of novel forms are often obscure because early and transitional fossils tend to be rare, poorly preserved, or lack proper phylogenetic contexts. We describe a new, exceptionally preserved enigmatic crab from the mid-Cretaceous of Colombia and the United States, whose completeness illuminates the early disparity of the group and the origins of novel forms. Its large and unprotected compound eyes, small fusiform body, and leg-like mouthparts suggest larval trait retention into adulthood via heterochronic development (pedomorphosis), while its large oar-like legs represent the earliest known adaptations in crabs for active swimming. Our phylogenetic analyses, including representatives of all major lineages of fossil and extant crabs, challenge conventional views of their evolution by revealing multiple convergent losses of a typical "crab-like" body plan since the Early Cretaceous. These parallel morphological transformations may be associated with repeated invasions of novel environments, including the pelagic/necto-benthic zone in this pedomorphic chimera crab.

Early bursts of disparity and the reorganization of character integration

'Early bursts' of morphological disparity (i.e. diversity of anatomical types) are common in the fossil record. We typically model such bursts as elevated early rates of independent character change. Developmental theory predicts that modules of linked characters can change together, which would mimic the effects of elevated independent rates on disparity. However, correlated change introducing suboptimal states should encourage breakup (parcellation) of character suites allowing new (or primitive) states to evolve until new suites arise (relinkage). Thus, correlated change-breakup-relinkage presents mechanisms for early bursts followed by constrained evolution. Here, I analyse disparity in 257 published character matrices of fossil taxa. For each clade, I use inverse-modelling to infer most probably rates of independent change given both time-homogeneous and separate 'early versus late' rates. These rates are used to estimate expected disparity given both independent change models. The correlated change-breakup-relinkage model also predicts elevated frequencies of compatible character state-pairs appearing out of order in the fossil record (e.g. 01 appearing after 00 and 11; = low stratigraphic compatibility), as one solution to suboptimal states induced by correlated change is a return to states held before that change. As predicted by the correlated change-breakup-relinkage model, early disparity in the majority of clades both exceeds the expectations of either independent change model and excess early disparity correlates with low stratigraphic compatibility among character-pairs. Although it is possible that other mechanisms for linking characters contribute to these patterns, these results corroborate the idea that reorganization of developmental linkages is often associated with the origin of groups that biologists recognize as new higher taxa and that such reorganization offers a source of new disparity throughout the Phanerozoic.

The origin of squamates revealed by a Middle Triassic lizard from the Italian Alps

Modern squamates (lizards, snakes and amphisbaenians) are the world's most diverse group of tetrapods along with birds ¹ and have a long evolutionary history, with the oldest known fossils dating from the Middle Jurassic period-168 million years ago^2-4. The evolutionary origin of squamates is contentious because of several issues: (1) a fossil gap of approximately 70 million years exists between the oldest known fossils and their estimated origin^5-7; (2) limited sampling of squamates in reptile phylogenies; and (3) conflicts between morphological and molecular hypotheses regarding the origin of crown squamates^6,8,9. Here we shed light on these problems by using high-resolution microfocus X-ray computed tomography data from the articulated fossil reptile Megachirella wachtleri (Middle Triassic period, Italian Alps ¹⁰ ). We also present a phylogenetic dataset, combining fossils and extant taxa, and morphological and molecular data. We analysed this dataset under different optimality criteria to assess diapsid reptile relationships and the origins of squamates. Our results re-shape the diapsid phylogeny and present evidence that M. wachtleri is the oldest known stem squamate. Megachirella is 75 million years older than the previously known oldest squamate fossils, partially filling the fossil gap in the origin of lizards, and indicates a more gradual acquisition of squamatan features in diapsid evolution than previously thought. For the first time, to our knowledge, morphological and molecular data are in agreement regarding early squamate evolution, with geckoes-and not iguanians-as the earliest crown clade squamates. Divergence time estimates using relaxed combined morphological and molecular clocks show that lepidosaurs and most other diapsids originated before the Permian/Triassic extinction event, indicating that the Triassic was a period of radiation, not origin, for several diapsid lineages.

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

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

Mosasauroid phylogeny under multiple phylogenetic methods provides new insights on the evolution of aquatic adaptations in the group

Mosasauroids were a successful lineage of squamate reptiles (lizards and snakes) that radiated during the Late Cretaceous (95-66 million years ago). They can be considered one of the few lineages in the evolutionary history of tetrapods to have acquired a fully aquatic lifestyle, similarly to whales, ichthyosaurs and plesiosaurs. Despite a long history of research on this group, their phylogenetic relationships have only been tested so far using traditional (unweighted) maximum parsimony. However, hypotheses of mosasauroid relationships and the recently proposed multiple origins of aquatically adapted pelvic and pedal features in this group can be more thoroughly tested by methods that take into account variation in branch lengths and evolutionary rates. In this study, we present the first mosasauroid phylogenetic analysis performed under different analytical methods, including maximum likelihood, Bayesian inference, and implied weighting maximum parsimony. The results indicate a lack of congruence in the topological position of halisaurines and Dallasaurus. Additionally, the genus Prognathodon is paraphyletic under all hypotheses. Interestingly, a number of traditional mosasauroid clades become weakly supported, or unresolved, under Bayesian analyses. The reduced resolutions in some consensus trees create ambiguities concerning the evolution of fully aquatic pelvic/pedal conditions under many analyses. However, when enough resolution was obtained, reversals of the pelvic/pedal conditions were favoured by parsimony and likelihood ancestral state reconstructions instead of independent origins of aquatic features in mosasauroids. It is concluded that most of the observed discrepancies among the results can be associated with different analytical procedures, but also due to limited postcranial data on halisaurines, yaguarasaurines and Dallasaurus.

Integrating Paleontological and Phylogenetic Approaches to Macroevolution

With proliferation of molecular phylogenies and advances in statistical modeling, phylogeneticists can now address macroevolutionary questions that had traditionally been the purview of paleontology. Interest has focused on three areas at the intersection of phylogenetic and paleontological research: time-scaling phylogenies, understanding trait evolution, and modeling species diversification. Fossil calibrations have long been crucial for scaling phylogenies to absolute time, but recent advances allow more equal integration of extinct taxa. Simulation and empirical studies have shown that fossil data can markedly improve inferences about trait evolution, especially for models with heterogeneous temporal dynamics and in clades for which the living forms are unrepresentative remnants of their larger clade. Recent years have also seen a productive cross-disciplinary conversation about the nature and uncertainties of inferring diversification dynamics. Challenges remain, but the present time represents a flowering of interest in integrating these two views on the history of life.

The first iguanian lizard from the Mesozoic of Africa

The fossil record shows that iguanian lizards were widely distributed during the Late Cretaceous. However, the biogeographic history and early evolution of one of its most diverse and peculiar clades (acrodontans) remain poorly known. Here, we present the first Mesozoic acrodontan from Africa, which also represents the oldest iguanian lizard from that continent. The new taxon comes from the Kem Kem Beds in Morocco (Cenomanian, Late Cretaceous) and is based on a partial lower jaw. The new taxon presents a number of features that are found only among acrodontan lizards and shares greatest similarities with uromastycines, specifically. In a combined evidence phylogenetic dataset comprehensive of all major acrodontan lineages using multiple tree inference methods (traditional and implied weighting maximum-parsimony, and Bayesian inference), we found support for the placement of the new species within uromastycines, along with Gueragama sulamericana (Late Cretaceous of Brazil). The new fossil supports the previously hypothesized widespread geographical distribution of acrodontans in Gondwana during the Mesozoic. Additionally, it provides the first fossil evidence of uromastycines in the Cretaceous, and the ancestry of acrodontan iguanians in Africa.

Morphological Phylogenetics in the Genomic Age

Evolutionary trees underpin virtually all of biology, and the wealth of new genomic data has enabled us to reconstruct them with increasing detail and confidence. While phenotypic (typically morphological) traits are becoming less important in reconstructing evolutionary trees, they still serve vital and unique roles in phylogenetics, even for living taxa for which vast amounts of genetic information are available. Morphology remains a powerful independent source of evidence for testing molecular clades, and - through fossil phenotypes - the primary means for time-scaling phylogenies. Morphological phylogenetics is therefore vital for transforming undated molecular topologies into dated evolutionary trees. However, if morphology is to be employed to its full potential, biologists need to start scrutinising phenotypes in a more objective fashion, models of phenotypic evolution need to be improved, and approaches for analysing phenotypic traits and fossils together with genomic data need to be refined.

RevBayes: Bayesian Phylogenetic Inference Using Graphical Models and an Interactive Model-Specification Language

Programs for Bayesian inference of phylogeny currently implement a unique and fixed suite of models. Consequently, users of these software packages are simultaneously forced to use a number of programs for a given study, while also lacking the freedom to explore models that have not been implemented by the developers of those programs. We developed a new open-source software package, RevBayes, to address these problems. RevBayes is entirely based on probabilistic graphical models, a powerful generic framework for specifying and analyzing statistical models. Phylogenetic-graphical models can be specified interactively in RevBayes, piece by piece, using a new succinct and intuitive language called Rev. Rev is similar to the R language and the BUGS model-specification language, and should be easy to learn for most users. The strength of RevBayes is the simplicity with which one can design, specify, and implement new and complex models. Fortunately, this tremendous flexibility does not come at the cost of slower computation; as we demonstrate, RevBayes outperforms competing software for several standard analyses. Compared with other programs, RevBayes has fewer black-box elements. Users need to explicitly specify each part of the model and analysis. Although this explicitness may initially be unfamiliar, we are convinced that this transparency will improve understanding of phylogenetic models in our field. Moreover, it will motivate the search for improvements to existing methods by brazenly exposing the model choices that we make to critical scrutiny. RevBayes is freely available at http://www.RevBayes.com [Bayesian inference; Graphical models; MCMC; statistical phylogenetics.].