Ignoring stratigraphic age uncertainty leads to erroneous estimates of species divergence times under the fossilized birth-death process

Practical guidelines for Bayesian phylogenetic inference using Markov chain Monte Carlo (MCMC)

Phylogenetic estimation is, and has always been, a complex endeavor. Estimating a phylogenetic tree involves evaluating many possible solutions and possible evolutionary histories that could explain a set of observed data, typically by using a model of evolution. Values for all model parameters need to be evaluated as well. Modern statistical methods involve not just the estimation of a tree, but also solutions to more complex models involving fossil record information and other data sources. Markov chain Monte Carlo (MCMC) is a leading method for approximating the posterior distribution of parameters in a mathematical model. It is deployed in all Bayesian phylogenetic tree estimation software. While many researchers use MCMC in phylogenetic analyses, interpreting results and diagnosing problems with MCMC remain vexing issues to many biologists. In this manuscript, we will offer an overview of how MCMC is used in Bayesian phylogenetic inference, with a particular emphasis on complex hierarchical models, such as the fossilized birth-death (FBD) model. We will discuss strategies to diagnose common MCMC problems and troubleshoot difficult analyses, in particular convergence issues. We will show how the study design, the choice of models and priors, but also technical features of the inference tools themselves can all be adjusted to obtain the best results. Finally, we will also discuss the unique challenges created by the incorporation of fossil information in phylogenetic inference, and present tips to address them.

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.

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.

The ClaDS rate-heterogeneous birth-death prior for full phylogenetic inference in BEAST2

Bayesian phylogenetic inference requires a tree prior, which models the underlying diversification process that gives rise to the phylogeny. Existing birth-death diversification models include a wide range of features, for instance, lineage-specific variations in speciation and extinction (SSE) rates. While across-lineage variation in SSE rates is widespread in empirical datasets, few heterogeneous rate models have been implemented as tree priors for Bayesian phylogenetic inference. As a consequence, rate heterogeneity is typically ignored when reconstructing phylogenies, and rate heterogeneity is usually investigated on fixed trees. In this paper, we present a new BEAST2 package implementing the cladogenetic diversification rate shift (ClaDS) model as a tree prior. ClaDS is a birth-death diversification model designed to capture small progressive variations in birth and death rates along a phylogeny. Unlike previous implementations of ClaDS, which were designed to be used with fixed, user-chosen phylogenies, our package is implemented in the BEAST2 framework and thus allows full phylogenetic inference, where the phylogeny and model parameters are co-estimated from a molecular alignment. Our package provides all necessary components of the inference, including a new tree object and operators to propose moves to the Monte-Carlo Markov chain. It also includes a graphical interface through BEAUti. We validate our implementation of the package by comparing the produced distributions to simulated data and show an empirical example of the full inference, using a dataset of cetaceans.

Impacts of Taxon-Sampling Schemes on Bayesian Tip Dating Under the Fossilized Birth-Death Process

Evolutionary timescales can be inferred by molecular-clock analyses of genetic data and fossil evidence. Bayesian phylogenetic methods such as tip dating provide a powerful framework for inferring evolutionary timescales, but the most widely used priors for tree topologies and node times often assume that present-day taxa have been sampled randomly or exhaustively. In practice, taxon sampling is often carried out so as to include representatives of major lineages, such as orders or families. We examined the impacts of different densities of diversified sampling on Bayesian tip dating on unresolved fossilized birth-death (FBD) trees, in which fossil taxa are topologically constrained but their exact placements are averaged out. We used synthetic data generated by simulations of nucleotide sequence evolution, fossil occurrences, and diversified taxon sampling. Our analyses under the diversified-sampling FBD process show that increasing taxon-sampling density does not necessarily improve divergence-time estimates. However, when informative priors were specified for the root age or when tree topologies were fixed to those used for simulation, the performance of tip dating on unresolved FBD trees maintains its accuracy and precision or improves with taxon-sampling density. By exploring three situations in which models are mismatched, we find that including all relevant fossils, without pruning off those that are incompatible with the diversified-sampling FBD process, can lead to underestimation of divergence times. Our reanalysis of a eutherian mammal data set confirms some of the findings from our simulation study, and reveals the complexity of diversified taxon sampling in phylogenomic data sets. In highlighting the interplay of taxon-sampling density and other factors, the results of our study have practical implications for using Bayesian tip dating to infer evolutionary timescales across the Tree of Life. [Bayesian tip dating; eutherian mammals; fossilized birth-death process; phylogenomics; taxon sampling.].

Estimating the Age of Poorly Dated Fossil Specimens and Deposits Using a Total-Evidence Approach and the Fossilized Birth-Death Process

Bayesian total-evidence approaches under the fossilized birth-death model enable biologists to combine fossil and extant data while accounting for uncertainty in the ages of fossil specimens, in an integrative phylogenetic analysis. Fossil age uncertainty is a key feature of the fossil record as many empirical data sets may contain a mix of precisely dated and poorly dated fossil specimens or deposits. In this study, we explore whether reliable age estimates for fossil specimens can be obtained from Bayesian total-evidence phylogenetic analyses under the fossilized birth-death model. Through simulations based on the example of the Baltic amber deposit, we show that estimates of fossil ages obtained through such an analysis are accurate, particularly when the proportion of poorly dated specimens remains low and the majority of fossil specimens have precise dates. We confirm our results using an empirical data set of living and fossil penguins by artificially increasing the age uncertainty around some fossil specimens and showing that the resulting age estimates overlap with the recorded age ranges. Our results are applicable to many empirical data sets where classical methods of establishing fossil ages have failed, such as the Baltic amber and the Gobi Desert deposits. [Bayesian phylogenetic inference; fossil age estimates; fossilized birth-death; Lagerstätte; total-evidence.].

The importance of the Andes in the evolutionary radiation of Sigmodontinae (Rodentia, Cricetidae), the most diverse group of mammals in the Neotropics

The Andean mountains stand out for their striking species richness and endemicity that characterize many emblematic Neotropical clades distributed in or around these mountains. The radiation of the Sigmodontinae subfamily, the most diversified mammalian group in the Neotropics, has been historically related to Andean orogenesis. We aim to evaluate this interplay between geological processes and biological responses through the diversification dynamics, the biogeographical history, and the range evolution of the subfamily. For these, we built the most comprehensive phylogeny and gathered 14,836 occurrences for the subfamily. We identified one shift in the speciation rate in the genus Akodon, which suffered their Andean radiation after the arrival of non-Andean ancestors. Our biogeographic analyses show multiple dispersal paths throughout the evolution that allowed this subfamily to colonize all Neotropics. The Northern Andes and Central-Southern Andes were the most important sources of diversity. In addition, the Central-Southern Andes were the most relevant sink, receiving the highest number of lineages. The Andean region exhibited higher speciation and turnover rates than non-Andean regions. Thus, our results support the crucial role of the Andean Mountains in the Sigmodontinae radiation, acting as a "macroevolutionary cradle" and "species attractor" for several sigmodontine lineages at different times, and as a "species pump" becoming the biogeographic source of multiple widely distributed neotropical lineages. Then, complex macroevolutionary dynamics would explain these rodents' high extant Andean diversity and their wide distribution in the Neotropics.

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.

The Occurrence Birth-Death Process for combined-evidence analysis in macroevolution and epidemiology

Phylodynamic models generally aim at jointly inferring phylogenetic relationships, model parameters, and more recently, the number of lineages through time, based on molecular sequence data. In the fields of epidemiology and macroevolution, these models can be used to estimate, respectively, the past number of infected individuals (prevalence) or the past number of species (paleodiversity) through time. Recent years have seen the development of “total-evidence” analyses, which combine molecular and morphological data from extant and past sampled individuals in a unified Bayesian inference framework. Even sampled individuals characterized only by their sampling time, that is, lacking morphological and molecular data, which we call occurrences, provide invaluable information to estimate the past number of lineages. Here, we present new methodological developments around the fossilized birth–death process enabling us to (i) incorporate occurrence data in the likelihood function; (ii) consider piecewise-constant birth, death, and sampling rates; and (iii) estimate the past number of lineages, with or without knowledge of the underlying tree. We implement our method in the RevBayes software environment, enabling its use along with a large set of models of molecular and morphological evolution, and validate the inference workflow using simulations under a wide range of conditions. We finally illustrate our new implementation using two empirical data sets stemming from the fields of epidemiology and macroevolution. In epidemiology, we infer the prevalence of the coronavirus disease 2019 outbreak on the Diamond Princess ship, by taking into account jointly the case count record (occurrences) along with viral sequences for a fraction of infected individuals. In macroevolution, we infer the diversity trajectory of cetaceans using molecular and morphological data from extant taxa, morphological data from fossils, as well as numerous fossil occurrences. The joint modeling of occurrences and trees holds the promise to further bridge the gap between traditional epidemiology and pathogen genomics, as well as paleontology and molecular phylogenetics. [Birth–death model; epidemiology; fossils; macroevolution; occurrences; phylogenetics; skyline.]

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.

Inferring the Total-Evidence Timescale of Marattialean Fern Evolution in the Face of Model Sensitivity

Phylogenetic divergence-time estimation has been revolutionized by two recent developments: 1) total-evidence dating (or "tip-dating") approaches that allow for the incorporation of fossils as tips in the analysis, with their phylogenetic and temporal relationships to the extant taxa inferred from the data and 2) the fossilized birth-death (FBD) class of tree models that capture the processes that produce the tree (speciation, extinction, and fossilization) and thus provide a coherent and biologically interpretable tree prior. To explore the behavior of these methods, we apply them to marattialean ferns, a group that was dominant in Carboniferous landscapes prior to declining to its modest extant diversity of slightly over 100 species. We show that tree models have a dramatic influence on estimates of both divergence times and topological relationships. This influence is driven by the strong, counter-intuitive informativeness of the uniform tree prior, and the inherent nonidentifiability of divergence-time models. In contrast to the strong influence of the tree models, we find minor effects of differing the morphological transition model or the morphological clock model. We compare the performance of a large pool of candidate models using a combination of posterior-predictive simulation and Bayes factors. Notably, an FBD model with epoch-specific speciation and extinction rates was strongly favored by Bayes factors. Our best-fitting model infers stem and crown divergences for the Marattiales in the mid-Devonian and Late Cretaceous, respectively, with elevated speciation rates in the Mississippian and elevated extinction rates in the Cisuralian leading to a peak diversity of ${\sim}$2800 species at the end of the Carboniferous, representing the heyday of the Psaroniaceae. This peak is followed by the rapid decline and ultimate extinction of the Psaroniaceae, with their descendants, the Marattiaceae, persisting at approximately stable levels of diversity until the present. This general diversification pattern appears to be insensitive to potential biases in the fossil record; despite the preponderance of available fossils being from Pennsylvanian coal balls, incorporating fossilization-rate variation does not improve model fit. In addition, by incorporating temporal data directly within the model and allowing for the inference of the phylogenetic position of the fossils, our study makes the surprising inference that the clade of extant Marattiales is relatively young, younger than any of the fossils historically thought to be congeneric with extant species. This result is a dramatic demonstration of the dangers of node-based approaches to divergence-time estimation, where the assignment of fossils to particular clades is made a priori (earlier node-based studies that constrained the minimum ages of extant genera based on these fossils resulted in much older age estimates than in our study) and of the utility of explicit models of morphological evolution and lineage diversification. [Bayesian model comparison; Carboniferous; divergence-time estimation; fossil record; fossilized birth-death; lineage diversification; Marattiales; models of morphological evolution; Psaronius; RevBayes.].

An early burst in brachiopod evolution corresponding with significant climatic shifts during the Great Ordovician Biodiversification Event

We employ modified tip-dating methods to date divergence times within the Strophomenoidea, one of the most abundant and species-rich brachiopod clades to radiate during the Great Ordovician Biodiversification Event (GOBE), to determine if significant environmental changes at this time correlate with the diversification of the clade. Models using origination, extinction and sampling rates to estimate prior probabilities of divergence times strongly support both high rates of anatomical change per million years and rapid divergences shortly before the clade first appears in the fossil record. These divergence times indicate much higher rates of cladogenesis than are typical of brachiopods during this interval. The correspondence of high speciation rates and high anatomical disparity suggests punctuated (speciational) change drove the high frequencies of early anatomical change, which in turn suggests increased ecological opportunities rather than shifting developmental constraints account for high rates of anatomical change. The pulse of rapid evolution began coincident with cooling temperatures, the start of major oscillations in sea level and increased levels of atmospheric oxygen. Our results suggest that these factors permitted major geographical and ecological expansion of strophomenoids with intervals of geographical isolation, resulting in elevated speciation rates and corresponding elevated frequencies of punctuated change.

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.

Including fossils in phylogeny: a glimpse into the evolution of the superfamily Evanioidea (Hymenoptera: Apocrita) under tip-dating and the fossilized birth–death process

Using a fossilized birth–death model, a new phylogeny of the superfamily Evanioidea (including ensign wasps, nightshade wasps and hatchet wasps) is proposed, with estimates of divergence times for its constitutive families and for corroborating the monophyly of Evanioidea. Additionally, our Bayesian analyses demonstrate the monophyly of †Anomopterellidae, †Othniodellithidae, †Andreneliidae, Aulacidae, Gasteruptiida and Evaniidae, whereas †Praeaulacidae and †Baissidae appear to be paraphyletic. Vectevania vetula and Hyptiogastrites electrinus are transferred to Aulacidae. We estimate the divergence time of Evanioidea to be in the Late Triassic (~203 Mya). Additionally, three new othniodellithid wasps are described and figured from mid-Cretaceous Burmese amber as the new genus Keratodellitha, with three new species: Keratodellitha anubis sp. nov., Keratodellitha basilisci sp. nov. and Keratodellitha kirin sp. nov. We also document a temporal shift in relative species richness between Ichneumonoidea and Evanioidea.

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.

The Making of Calibration Sausage Exemplified by Recalibrating the Transcriptomic Timetree of Jawed Vertebrates

Molecular divergence dating has the potential to overcome the incompleteness of the fossil record in inferring when cladogenetic events (splits, divergences) happened, but needs to be calibrated by the fossil record. Ideally but unrealistically, this would require practitioners to be specialists in molecular evolution, in the phylogeny and the fossil record of all sampled taxa, and in the chronostratigraphy of the sites the fossils were found in. Paleontologists have therefore tried to help by publishing compendia of recommended calibrations, and molecular biologists unfamiliar with the fossil record have made heavy use of such works (in addition to using scattered primary sources and copying from each other). Using a recent example of a large node-dated timetree inferred from molecular data, I reevaluate all 30 calibrations in detail, present the current state of knowledge on them with its various uncertainties, rerun the dating analysis, and conclude that calibration dates cannot be taken from published compendia or other secondary or tertiary sources without risking strong distortions to the results, because all such sources become outdated faster than they are published: 50 of the (primary) sources I cite to constrain calibrations were published in 2019, half of the total of 280 after mid-2016, and 90% after mid-2005. It follows that the present work cannot serve as such a compendium either; in the slightly longer term, it can only highlight known and overlooked problems. Future authors will need to solve each of these problems anew through a thorough search of the primary paleobiological and chronostratigraphic literature on each calibration date every time they infer a new timetree, and that literature is not optimized for that task, but largely has other objectives.

Fossil palm reading: using fruits to reveal the deep roots of palm diversity

PREMISE

Fossils are essential for understanding evolutionary history because they provide direct evidence of past diversity and geographic distributions. However, resolving systematic relationships between fossils and extant taxa, an essential step for many macroevolutionary studies, requires extensive comparative work on morphology and anatomy. While palms (Arecaceae) have an excellent fossil record that includes numerous fossil fruits, many are difficult to identify due in part to limited comparative data on modern fruit structure.

METHODS

We studied fruits of 207 palm species, representing nearly every modern genus, using X-ray microcomputed tomography. We then developed a morphological data set to test whether the fossil record of fruits can improve our understanding of palm diversity in the deep past. To evaluate the accuracy with which this data set recovers systematic relationships, we performed phylogenetic pseudofossilization analyses. We then used the data set to investigate the phylogenetic relationships of five previously published fossil palm fruits.

RESULTS

Phylogenetic analyses of fossils and pseudofossilization of extant taxa show that fossils can be placed accurately to the tribe and subtribe level with this data set, but node support must be considered. The phylogenetic relationships of the fossils suggest origins of many modern lineages in the Cretaceous and early Paleogene. Three of these fossils are suitable as new node calibrations for palms.

CONCLUSIONS

This work improves our knowledge of fruit structure in palms, lays a foundation for applying fossil fruits to macroevolutionary studies, and provides new insights into the evolutionary history and early diversification of Arecaceae.

Mind the Outgroup and Bare Branches in Total-Evidence Dating: a Case Study of Pimpliform Darwin Wasps (Hymenoptera, Ichneumonidae)

Taxon sampling is a central aspect of phylogenetic study design, but it has received limited attention in the context of total-evidence dating, a widely used dating approach that directly integrates molecular and morphological information from extant and fossil taxa. We here assess the impact of commonly employed outgroup sampling schemes and missing morphological data in extant taxa on age estimates in a total-evidence dating analysis under the uniform tree prior. Our study group is Pimpliformes, a highly diverse, rapidly radiating group of parasitoid wasps of the family Ichneumonidae. We analyze a data set comprising 201 extant and 79 fossil taxa, including the oldest fossils of the family from the Early Cretaceous and the first unequivocal representatives of extant subfamilies from the mid-Paleogene. Based on newly compiled molecular data from ten nuclear genes and a morphological matrix that includes 222 characters, we show that age estimates become both older and less precise with the inclusion of more distant and more poorly sampled outgroups. These outgroups not only lack morphological and temporal information but also sit on long terminal branches and considerably increase the evolutionary rate heterogeneity. In addition, we discover an artifact that might be detrimental for total-evidence dating: "bare-branch attraction," namely high attachment probabilities of certain fossils to terminal branches for which morphological data are missing. Using computer simulations, we confirm the generality of this phenomenon and show that a large phylogenetic distance to any of the extant taxa, rather than just older age, increases the risk of a fossil being misplaced due to bare-branch attraction. After restricting outgroup sampling and adding morphological data for the previously attracting, bare branches, we recover a Jurassic origin for Pimpliformes and Ichneumonidae. This first age estimate for the group not only suggests an older origin than previously thought but also that diversification of the crown group happened well before the Cretaceous-Paleogene boundary. Our case study demonstrates that in order to obtain robust age estimates, total-evidence dating studies need to be based on a thorough and balanced sampling of both extant and fossil taxa, with the aim of minimizing evolutionary rate heterogeneity and missing morphological information. [Bare-branch attraction; ichneumonids; fossils; morphological matrix; phylogeny; RoguePlots.].

The probability distribution of the ancestral population size conditioned on the reconstructed phylogenetic tree with occurrence data

We consider a homogeneous birth-death process with three different sampling schemes. First, individuals can be sampled through time and included in a reconstructed phylogenetic tree. Second, they can be sampled through time and only recorded as a point 'occurrence' along a timeline. Third, extant individuals can be sampled and included in the reconstructed phylogenetic tree with a fixed probability. We further consider that sampled individuals can be removed or not from the process, upon sampling, with fixed probability. We derive the probability distribution of the population size at any time in the past conditional on the joint observation of a reconstructed phylogenetic tree and a record of occurrences not included in the tree. We also provide an algorithm to simulate ancestral population size trajectories given the observation of a reconstructed phylogenetic tree and occurrences. This distribution can be readily used to draw inferences about the ancestral population size in the field of epidemiology and macroevolution. In epidemiology, these results will allow data from epidemiological case count studies to be used in conjunction with molecular sequencing data (yielding reconstructed phylogenetic trees) to coherently estimate prevalence through time. In macroevolution, it will foster the joint examination of the fossil record and extant taxa to reconstruct past biodiversity.

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

BACKGROUND

The vast majority of all life that ever existed on earth is now extinct and several aspects of their evolutionary history can only be assessed by using morphological data from the fossil record. Sphenodontian reptiles are a classic example, having an evolutionary history of at least 230 million years, but currently represented by a single living species (Sphenodon punctatus). Hence, it is imperative to improve the development and implementation of probabilistic models to estimate evolutionary trees from morphological data (e.g., morphological clocks), which has direct benefits to understanding relationships and evolutionary patterns for both fossil and living species. However, the impact of model choice on morphology-only datasets has been poorly explored.

RESULTS

Here, we investigate the impact of a wide array of model choices on the inference of evolutionary trees and macroevolutionary parameters (divergence times and evolutionary rates) using a new data matrix on sphenodontian reptiles. Specifically, we tested different clock models, clock partitioning, taxon sampling strategies, sampling for ancestors, and variations on the fossilized birth-death (FBD) tree model parameters through time. We find a strong impact on divergence times and background evolutionary rates when applying widely utilized approaches, such as allowing for ancestors in the tree and the inappropriate assumption of diversification parameters being constant through time. We compare those results with previous studies on the impact of model choice to molecular data analysis and provide suggestions for improving the implementation of morphological clocks. Optimal model combinations find the radiation of most major lineages of sphenodontians to be in the Triassic and a gradual but continuous drop in morphological rates of evolution across distinct regions of the phenotype throughout the history of the group.

CONCLUSIONS

We provide a new hypothesis of sphenodontian classification, along with detailed macroevolutionary patterns in the evolutionary history of the group. Importantly, we provide suggestions to avoid overestimated divergence times and biased parameter estimates using morphological clocks. Partitioning relaxed clocks offers methodological limitations, but those can be at least partially circumvented to reveal a detailed assessment of rates of evolution across the phenotype and tests of evolutionary mosaicism.

First monk seal from the Southern Hemisphere rewrites the evolutionary history of true seals

Living true seals (phocids) are the most widely dispersed semi-aquatic marine mammals, and comprise geographically separate northern (phocine) and southern (monachine) groups. Both are thought to have evolved in the North Atlantic, with only two monachine lineages-elephant seals and lobodontins-subsequently crossing the equator. The third and most basal monachine tribe, the monk seals, have hitherto been interpreted as exclusively northern and (sub)tropical throughout their entire history. Here, we describe a new species of extinct monk seal from the Pliocene of New Zealand, the first of its kind from the Southern Hemisphere, based on one of the best-preserved and richest samples of seal fossils worldwide. This unanticipated discovery reveals that all three monachine tribes once coexisted south of the equator, and forces a profound revision of their evolutionary history: rather than primarily diversifying in the North Atlantic, monachines largely evolved in the Southern Hemisphere, and from this southern cradle later reinvaded the north. Our results suggest that true seals crossed the equator over eight times in their history. Overall, they more than double the age of the north-south dichotomy characterizing living true seals and confirms a surprisingly recent major change in southern phocid diversity.

Tip dating with fossil sites and stratigraphic sequences

Tip dating, a method of phylogenetic analysis in which fossils are included as terminals and assigned an age, is becoming increasingly widely used in evolutionary studies. Current implementations of tip dating allow fossil ages to be assigned as a point estimate, or incorporate uncertainty through the use of uniform tip age priors. However, the use of tip age priors has the unwanted effect of decoupling the ages of fossils from the same fossil site. Here we introduce a new Markov Chain Monte Carlo (MCMC) proposal, which allows fossils from the same site to have linked ages, while still incorporating uncertainty in the age of the fossil site itself. We also include an extension, allowing fossil sites to be ordered in a stratigraphic column with age bounds applied only to the top and bottom of the sequence. These MCMC proposals are implemented in a new open-source BEAST2 package, palaeo. We test these new proposals on a dataset of early vertebrate fossils, concentrating on the effects on two sites with multiple acanthodian fossil taxa but wide age uncertainty, the Man On The Hill (MOTH) site from northern Canada, and the Turin Hill site from Scotland, both of Lochkovian (Early Devonian) age. The results show an increased precision of age estimates when fossils have linked tip ages compared to when ages are unlinked, and in this example leads to support for a younger age for the MOTH site compared with the Turin Hill site. There is also a minor effect on the tree topology of acanthodians. These new MCMC proposals should be widely applicable to studies that employ tip dating, particularly when the terminals are coded as individual specimens.

Olson's Gap or Olson's Extinction? A Bayesian tip-dating approach to resolving stratigraphic uncertainty

Adaptive radiations and mass extinctions are of critical importance in structuring terrestrial ecosystems. However, the causes and progress of these transitions often remain controversial, in part because of debates surrounding the completeness of the fossil record and biostratigraphy of the relevant fossil-bearing formations. The early-middle Permian, when a substantial faunal turnover in tetrapods coincided with a restructuring of the trophic structure of ecosystems, is such a time. Some have suggested the transition is obscured by a gap in the tetrapod fossil record (Olson's Gap), while others suggest a correlation between North American and Russian tetrapod-bearing formations allows the interval to be documented in detail. The latter biostratigraphic scheme has been used to support a mass extinction at this time (Olson's Extinction). Bayesian tip-dating methods used frequently in phylogenetics are employed to resolve this debate. Bayes factors are used to compare the results of analyses incorporating tip age priors based on different stratigraphic hypotheses, to show which stratigraphic scheme best fits the morphological data and phylogeny. Olson's Gap is rejected, and the veracity of Olson's Extinction is given further support. Tip-dating approaches have great potential to resolve debates surrounding the stratigraphic ages of critical formations where appropriate morphological data is available.

Estimating the evolutionary rates in mosasauroids and plesiosaurs: discussion of niche occupation in Late Cretaceous seas

Observations of temporal overlap of niche occupation among Late Cretaceous marine amniotes suggest that the rise and diversification of mosasauroid squamates might have been influenced by competition with or disappearance of some plesiosaur taxa. We discuss that hypothesis through comparisons of the rates of morphological evolution of mosasauroids throughout their evolutionary history with those inferred for contemporary plesiosaur clades. We used expanded versions of two species-level phylogenetic datasets of both these groups, updated them with stratigraphic information, and analyzed using the Bayesian inference to estimate the rates of divergence for each clade. The oscillations in evolutionary rates of the mosasauroid and plesiosaur lineages that overlapped in time and space were then used as a baseline for discussion and comparisons of traits that can affect the shape of the niche structures of aquatic amniotes, such as tooth morphologies, body size, swimming abilities, metabolism, and reproduction. Only two groups of plesiosaurs are considered to be possible niche competitors of mosasauroids: the brachauchenine pliosaurids and the polycotylid leptocleidians. However, direct evidence for interactions between mosasauroids and plesiosaurs is scarce and limited only to large mosasauroids as the predators/scavengers and polycotylids as their prey. The first mosasauroids differed from contemporary plesiosaurs in certain aspects of all discussed traits and no evidence suggests that early representatives of Mosasauroidea diversified after competitions with plesiosaurs. Nevertheless, some mosasauroids, such as tylosaurines, might have seized the opportunity and occupied the niche previously inhabited by brachauchenines, around or immediately after they became extinct, and by polycotylids that decreased their phylogenetic diversity and disparity around the time the large-sized tylosaurines started to flourish.

A Simulation-Based Evaluation of Tip-Dating Under the Fossilized Birth-Death Process

Bayesian molecular dating is widely used to study evolutionary timescales. This procedure usually involves phylogenetic analysis of nucleotide sequence data, with fossil-based calibrations applied as age constraints on internal nodes of the tree. An alternative approach is tip-dating, which explicitly includes fossil data in the analysis. This can be done, for example, through the joint analysis of molecular data from present-day taxa and morphological data from both extant and fossil taxa. In the context of tip-dating, an important development has been the fossilized birth-death process, which allows non-contemporaneous tips and sampled ancestors while providing a model of lineage diversification for the prior on the tree topology and internal node times. However, tip-dating with fossils faces a number of considerable challenges, especially, those associated with fossil sampling and evolutionary models for morphological characters. We conducted a simulation study to evaluate the performance of tip-dating using the fossilized birth-death model. We simulated fossil occurrences and the evolution of nucleotide sequences and morphological characters under a wide range of conditions. Our analyses of these data show that the number and the maximum age of fossil occurrences have a greater influence than the degree of among-lineage rate variation or the number of morphological characters on estimates of node times and the tree topology. Tip-dating with the fossilized birth-death model generally performs well in recovering the relationships among extant taxa but has difficulties in correctly placing fossil taxa in the tree and identifying the number of sampled ancestors. The method yields accurate estimates of the ages of the root and crown group, although the precision of these estimates varies with the probability of fossil occurrence. The exclusion of morphological characters results in a slight overestimation of node times, whereas the exclusion of nucleotide sequences has a negative impact on inference of the tree topology. Our results provide an overview of the performance of tip-dating using the fossilized birth-death model, which will inform further development of the method and its application to key questions in evolutionary biology.

Evolutionary Models for the Diversification of Placental Mammals Across the KPg Boundary

Deciphering the timing of the placental mammal radiation is a longstanding problem in evolutionary biology, but consensus on the tempo and mode of placental diversification remains elusive. Nevertheless, an accurate timetree is essential for understanding the role of important events in Earth history (e.g., Cretaceous Terrestrial Revolution, KPg mass extinction) in promoting the taxonomic and ecomorphological diversification of Placentalia. Archibald and Deutschman described three competing models for the diversification of placental mammals, which are the Explosive, Long Fuse, and Short Fuse Models. More recently, the Soft Explosive Model and Trans-KPg Model have emerged as additional hypotheses for the placental radiation. Here, we review molecular and paleontological evidence for each of these five models including the identification of general problems that can negatively impact divergence time estimates. The Long Fuse Model has received more support from relaxed clock studies than any of the other models, but this model is not supported by morphological cladistic studies that position Cretaceous eutherians outside of crown Placentalia. At the same time, morphological cladistics has a poor track record of reconstructing higher-level relationships among the orders of placental mammals including the results of new pseudoextinction analyses that we performed on the largest available morphological data set for mammals (4,541 characters). We also examine the strengths and weaknesses of different timetree methods (node dating, tip dating, and fossilized birth-death dating) that may now be applied to estimate the timing of the placental radiation. While new methods such as tip dating are promising, they also have problems that must be addressed if these methods are to effectively discriminate among competing hypotheses for placental diversification. Finally, we discuss the complexities of timetree estimation when the signal of speciation times is impacted by incomplete lineage sorting (ILS) and hybridization. Not accounting for ILS results in dates that are older than speciation events. Hybridization, in turn, can result in dates than are younger or older than speciation dates. Disregarding this potential variation in "gene" history across the genome can distort phylogenetic branch lengths and divergence estimates when multiple unlinked genomic loci are combined together in a timetree analysis.

Bayesian tip dating reveals heterogeneous morphological clocks in Mesozoic birds

Recently, comprehensive morphological datasets including nearly all the well-recognized Mesozoic birds became available, making it feasible for statistically rigorous methods to unveil finer evolutionary patterns during early avian evolution. Here, we exploited the advantage of Bayesian tip dating under relaxed morphological clocks to estimate both the divergence times and evolutionary rates while accounting for their uncertainties. We further subdivided the characters into six body regions (i.e. skull, axial skeleton, pectoral girdle and sternum, forelimb, pelvic girdle and hindlimb) to assess evolutionary rate heterogeneity both along the lineages and across partitions. We observed extremely high rates of morphological character changes during early avian evolution, and the clock rates are quite heterogeneous among the six regions. The branch subtending Pygostylia shows an extremely high rate in the axial skeleton, while the branches subtending Ornithothoraces and Enantiornithes show notably high rates in the pectoral girdle and sternum and moderately high rates in the forelimb. The extensive modifications in these body regions largely correspond to refinement of the flight capability. This study reveals the power and flexibility of Bayesian tip dating implemented in MrBayes to investigate evolutionary dynamics in deep time.