A total-evidence approach to dating with fossils, applied to the early radiation of the hymenoptera

Middle Jurassic insect mines on gymnosperms provide missing links to early mining evolution

We investigated the mining mode of insect feeding, involving larval consumption of a plant's internal tissues, from the Middle Jurassic (165 million years ago) Daohugou locality of Northeastern China. Documentation of mining from the Jurassic Period is virtually unknown, and results from this time interval would address mining evolution during the temporal gap of mine-seed plant diversifications from the previous Late Triassic to the subsequent Early Cretaceous. Plant fossils were examined with standard microscopic procedures for herbivory and used the standard functional feeding group-damage-type system of categorizing damage. All fossil mines were photographed and databased. We examined 2014 plant specimens, of which 27 occurrences on 14 specimens resulted in eight, new, mine damage types (DTs) present on six genera of bennettitalean, ginkgoalean, and pinalean gymnosperms. Three conclusions emerge from this study. First, these mid-Mesozoic mines are morphologically conservative and track plant host anatomical structure rather than plant phylogeny. Second, likely insect fabricators of these mines were three basal lineages of polyphagan beetles, four basal lineages of monotrysian moths, and a basal lineage tenthredinoid sawflies. Third, the nutrition hypothesis, indicating that miners had greater access to nutritious, inner tissues of new plant lineages, best explains mine evolution during the mid-Mesozoic.

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.

Fast Bayesian Inference of Phylogenies from Multiple Continuous Characters

Time-scaled phylogenetic trees are an ultimate goal of evolutionary biology and a necessary ingredient in comparative studies. The accumulation of genomic data has resolved the tree of life to a great extent, yet timing evolutionary events remain challenging if not impossible without external information such as fossil ages and morphological characters. Methods for incorporating morphology in tree estimation have lagged behind their molecular counterparts, especially in the case of continuous characters. Despite recent advances, such tools are still direly needed as we approach the limits of what molecules can teach us. Here, we implement a suite of state-of-the-art methods for leveraging continuous morphology in phylogenetics, and by conducting extensive simulation studies we thoroughly validate and explore our methods' properties. While retaining model generality and scalability, we make it possible to estimate absolute and relative divergence times from multiple continuous characters while accounting for uncertainty. We compile and analyze one of the most data-type diverse data sets to date, comprised of contemporaneous and ancient molecular sequences, and discrete and continuous morphological characters from living and extinct Carnivora taxa. We conclude by synthesizing lessons about our method's behavior, and suggest future research venues.

Jurassic shuotheriids show earliest dental diversification of mammaliaforms

Shuotheriids are Jurassic mammaliaforms that possess pseudotribosphenic teeth in which a pseudotalonid is anterior to the trigonid in the lower molar, contrasting with the tribosphenic pattern of therian mammals (placentals, marsupials and kin) in which the talonid is posterior to the trigonid1-4. The origin of the pseudotribosphenic teeth remains unclear, obscuring our perception of shuotheriid affinities and the early evolution of mammaliaforms1,5-9. Here we report a new Jurassic shuotheriid represented by two skeletal specimens. Their complete pseudotribosphenic dentitions allow reidentification of dental structures using serial homology and the tooth occlusal relationship. Contrary to the conventional view1,2,6,10,11, our findings show that dental structures of shuotheriids can be homologized to those of docodontans and partly support homologous statements for some dental structures between docodontans and other mammaliaforms6,12. The phylogenetic analysis based on new evidence removes shuotheriids from the tribosphenic ausktribosphenids (including monotremes) and clusters them with docodontans to form a new clade, Docodontiformes, that is characterized by pseudotribosphenic features. In the phylogeny, docodontiforms and 'holotherians' (Kuehneotherium, monotremes and therians)13 evolve independently from a Morganucodon-like ancestor with triconodont molars by labio-lingual widening their posterior teeth for more efficient food processing. The pseudotribosphenic pattern passed a cusp semitriangulation stage9, whereas the tribosphenic pattern and its precursor went through a stage of cusp triangulation. The two different processes resulted in complex tooth structures and occlusal patterns that elucidate the earliest diversification of mammaliaforms.

Origin and Early Evolution of Hydrocharitaceae and the Ancestral Role of Stratiotes

The combined morphological features of Stratiotes (Hydrocharitaceae) pollen, observed with light and electron microscopy, make it unique among all angiosperm pollen types and easy to identify. Unfortunately, the plant is (and most likely was) insect-pollinated and produces relatively few pollen grains per flower, contributing to its apparent absence in the paleopalynological record. Here, we present fossil Stratiotes pollen from the Eocene of Germany (Europe) and Kenya (Africa), representing the first reliable pre-Pleistocene pollen records of this genus worldwide and the only fossils of this family discovered so far in Africa. The fossil Stratiotes pollen grains are described and compared to pollen from a single modern species, Stratiotes aloides L. The paleophytogeographic significance and paleoecological aspects of these findings are discussed in relation to the Hydrocharitaceae fossil records and molecular phylogeny, as well as the present-day distribution patterns of its modern genera.

Climatic and biogeographic processes underlying the diversification of the pantropical flowering plant family Annonaceae

Tropical forests harbor the richest biodiversity among terrestrial ecosystems, but few studies have addressed the underlying processes of species diversification in these ecosystems. We use the pantropical flowering plant family Annonaceae as a study system to investigate how climate and biogeographic events contribute to diversification. A super-matrix phylogeny comprising 835 taxa (34% of Annonaceae species) based on eight chloroplast regions was used in this study. We show that global temperature may better explain the recent rapid diversification in Annonaceae than time and constant models. Accelerated accumulation of niche divergence (around 15 Ma) lags behind the increase of diversification rate (around 25 Ma), reflecting a heterogeneous transition to recent diversity increases. Biogeographic events are related to only two of the five diversification rate shifts detected. Shifts in niche evolution nevertheless appear to be associated with increasingly seasonal environments. Our results do not support the direct correlation of any particular climatic niche shifts or historical biogeographical event with shifts in diversification rate. Instead, we suggest that Annonaceae diversification can lead to later niche divergence as a result of increasing interspecific competition arising from species accumulation. Shifts in niche evolution appear to be associated with increasingly seasonal environments. Our results highlight the complexity of diversification in taxa with long evolutionary histories.

Tip dating and Bayes factors provide insight into the divergences of crown bird clades across the end-Cretaceous mass extinction

The origin of crown birds (Neornithes) remains contentious owing to conflicting divergence time hypotheses obtained from alternative sources of data. The fossil record suggests limited diversification of Neornithes in the Late Mesozoic and a substantial radiation in the aftermath of the Cretaceous-Palaeogene (K-Pg) mass extinction, approximately 66 Ma. Molecular clock studies, however, have yielded estimates for neornithine origins ranging from the Early Cretaceous (130 Ma) to less than 10 Myr before the K-Pg. We use Bayes factors to compare the fit of node ages from different molecular clock studies to an independent morphological dataset. Our results allow us to reject scenarios of crown bird origins deep in the Early Cretaceous, as well as an origin of crown birds within the last 10 Myr of the Cretaceous. The scenario best supported by our analyses is one where Neornithes originated between the Early and Late Cretaceous (ca 100 Ma), while numerous divergences within major neoavian clades either span or postdate the K-Pg. This study affirms the importance of the K-Pg on the diversification of modern birds, and the potential of combined-evidence tip-dating analyses to illuminate recalcitrant 'rocks versus clocks' debates.

Embracing uncertainty: The way forward in plant fossil phylogenetics

Although molecular phylogenetics remains the most widely used method of inferring the evolutionary history of living groups, the last decade has seen a renewed interest in morphological phylogenetics, mostly driven by the promises that integrating the fossil record in phylogenetic trees offers to our understanding of macroevolutionary processes and dynamics and the possibility that the inclusion of fossil taxa could lead to more accurate phylogenetic hypotheses. The plant fossil record presents some challenges to its integration in a phylogenetic framework. Phylogenies including plant fossils often retrieve uncertain relationships with low support, or lack of resolution. This low support is due to the pervasiveness of morphological convergence among plant organs and the fragmentary nature of many plant fossils, and it is often perceived as a fundamental weakness reducing the utility of plant fossils in phylogenetics. Here I discuss the importance of uncertainty in morphological phylogenetics and how we can identify important information from different patterns and types of uncertainty. I also review a set of methodologies that can allow us to understand the causes underpinning uncertainty and how these practices can help us to further our knowledge of plant fossils. I also propose that a new visual language, including the use of networks instead of trees, represents an improvement on the old visualization based on consensus trees and more adequately serves phylogeneticists working with plant fossils. This set of methods and visualization tools represents an important way forward in a fundamental field for our understanding of the evolutionary history of plants.

A Palaeogene stem crotaphytid (Aciprion formosum) and the phylogenetic affinities of early fossil pleurodontan iguanians

Pleurodonta is an ancient, diverse clade of iguanian lizard distributed primarily in the Western Hemisphere. Although the clade is a frequent subject of systematic research, phylogenetic resolution among the major pleurodontan clades is elusive. That uncertainty has complicated the interpretations of many fossil pleurodontans. I describe a fossil skull of a pleurodontan lizard from the Palaeogene of Wyoming that was previously allocated to the puzzling taxon Aciprion formosum, and provide an updated morphological matrix for iguanian lizards. Phylogenetic analyses using Bayesian inference demonstrate that the fossil skull is the oldest and first definitive stem member of Crotaphytidae (collared and leopard lizards), establishing the presence of that clade in North America during the Palaeogene. I also discuss new or revised hypotheses for the relationships of several early pleurodontans. In particular, I examine potential evidence for crown-Pleurodonta in the Cretaceous of Mongolia (Polrussia), stem Pleurodonta in the Cretaceous of North America (Magnuviator) and a stem anole in the Eocene of North America (Afairiguana). I suggest that the placement of the fossil crotaphytid is stable to the uncertain phylogeny of Pleurodonta, but recognize the dynamic nature of fossil diagnosis and the potential for updated systematic hypotheses for the other fossils analysed here.

Skyline Fossilized Birth-Death Model is Robust to Violations of Sampling Assumptions in Total-Evidence Dating

Several total-evidence dating studies under the fossilized birth-death (FBD) model have produced very old age estimates, which are not supported by the fossil record. This phenomenon has been termed "deep root attraction (DRA)." For two specific data sets, involving divergence time estimation for the early radiations of ants, bees, and wasps (Hymenoptera) and of placental mammals (Eutheria), it has been shown that the DRA effect can be greatly reduced by accommodating the fact that extant species in these trees have been sampled to maximize diversity, so-called diversified sampling. Unfortunately, current methods to accommodate diversified sampling only consider the extreme case where it is possible to identify a cut-off time such that all splits occurring before this time are represented in the sampled tree but none of the younger splits. In reality, the sampling bias is rarely this extreme and may be difficult to model properly. Similar modeling challenges apply to the sampling of the fossil record. This raises the question of whether it is possible to find dating methods that are more robust to sampling biases. Here, we show that the skyline FBD (SFBD) process, where the diversification and fossil-sampling rates can vary over time in a piecewise fashion, provides age estimates that are more robust to inadequacies in the modeling of the sampling process and less sensitive to DRA effects. In the SFBD model we consider, rates in different time intervals are either considered to be independent and identically distributed or assumed to be autocorrelated following an Ornstein-Uhlenbeck (OU) process. Through simulations and reanalyses of Hymenoptera and Eutheria data, we show that both variants of the SFBD model unify age estimates under random and diversified sampling assumptions. The SFBD model can resolve DRA by absorbing the deviations from the sampling assumptions into the inferred dynamics of the diversification process over time. Although this means that the inferred diversification dynamics must be interpreted with caution, taking sampling biases into account, we conclude that the SFBD model represents the most robust approach currently available for addressing DRA in total-evidence dating.

Measurement theory and paleobiology

Measurement theory, a branch of applied mathematics, offers guiding principles for extracting meaning from empirical observations and is applicable to any science involving measurements. Measurement theory is highly relevant in paleobiology because statistical approaches assuming ratio-scaled variables are commonly used on data belonging to nominal and ordinal scale types. We provide an informal introduction to representational measurement theory and argue for its importance in robust scientific inquiry. Although measurement theory is widely applicable in paleobiology research, we use the study of disparity to illustrate measurement theoretical challenges in the quantitative study of the fossil record. Respecting the inherent properties of different measurements enables meaningful inferences about evolutionary and ecological processes from paleontological data.

Phylogeny and paleobiogeography of the enigmatic North American primate Ekgmowechashala illuminated by new fossils from Nebraska (USA) and Guangxi Zhuang Autonomous Region (China)

Ekgmowechashala is a poorly documented but very distinctive primate known only from the late early Oligocene (early Arikareean) of western North America. Because of its highly autapomorphous dentition and spatiotemporal isolation, the phylogenetic and biogeographic affinities of Ekgmowechashala have long been debated. Here, we describe the oldest known fossils of Ekgmowechashala from the Brown Siltstone Beds of the Brule Formation, White River Group of western Nebraska. We also describe a new ekgmowechashaline taxon from the Nadu Formation (late Eocene) in the Baise Basin of Guangxi Zhuang Autonomous Region in southern China. Phylogenetic analysis suggests that North American Ekgmowechashala and the new Chinese taxon are sister taxa that are nested within a radiation of southern Asian adapiforms that also includes Gatanthropus, Muangthanhinius, and Bugtilemur. The new Chinese ekgmowechashaline helps fill the considerable disparity in dental morphology between Ekgmowechashala and more primitive ekgmowechashalids known from southern Asia. Our study underscores the fundamental role of southern Asia as a refugium for multiple primate clades during the cooler and drier climatic regime that prevailed after the Eocene-Oligocene transition. The colonization of North America by Ekgmowechashala helps define the beginning of the Arikareean Land Mammal Age and corresponds to an example of the Lazarus effect, whereby a taxon (in this case, the order Primates) reappears suddenly in the fossil record after a lengthy hiatus.

Geologically calibrated mammalian tree and its correlation with global events, including the emergence of humans

A robust timetree for Mammalia was constructed using the time calibration function of BEAST v1.10.4 and MEGA 11. The analysis involved the application of times of the most recent common ancestors, including a total of 19 mammalian fossil calibration ages following Benton et al. (Palaeontologia Electronica, 2015, 1-106) for their minimum ages. Additionally, fossil calibration ages for Gorilla, Pan, and a geologic event calibration age for otters were incorporated. Using these calibration ages, I constructed a geologically calibrated tree that estimates the age of the Homo and Pan splitting to be 5.69 Ma. The tree carries several significant implications. First, after the initial rifting at 120 Ma, the Atlantic Ocean expanded by over 500 km around Chron 34 (84 Ma), and vicariant speciation between Afrotheria (Africa) and Xenarthra (South America) appears to have commenced around 70 Ma. Additionally, ordinal level differentiations began immediately following the K-Pg boundary (66.0 Ma), supporting previous hypothesis that mammalian radiation rapidly filled ecological niches left vacant by non-avian dinosaurs. I constructed a diagram depicting the relationship between base substitution rate and age using an additional function in BEAST v1.10.4. The diagram reveals an exponential increase in the base substitution rate approaching recent times. This increased base substitution rate during the Neogene period may have contributed to the expansion of biodiversity, including the extensive adaptive radiation that led to the evolution of Homo sapiens. One significant driving factor behind this radiation could be attributed to the emergence and proliferation of C4 grasses since 20 Ma. These grasses have played a role in increasing carbon fixation, reducing atmospheric CO2 concentration, inducing global cooling, and initiating Quaternary glacial-interglacial cycles, thereby causing significant climatic changes.

Dating the evolution of the complex thalloid liverworts (Marchantiopsida): total-evidence dating analysis supports a Late Silurian-Early Devonian origin and post-Mesozoic morphological stasis

Divergence times based on molecular clock analyses often differ from those derived from total-evidence dating (TED) approaches. For bryophytes, fossils have been excluded from previous assessments of divergence times, and thus, their utility in dating analyses remains unexplored. Here, we conduct the first TED analyses of the complex thalloid liverworts (Marchantiopsida) that include fossils and evaluate macroevolutionary trends in morphological 'diversity' (disparity) and rates. Phylogenetic analyses were performed on a combined dataset of 130 discrete characters and 11 molecular markers (sampled from nuclear, plastid and mitochondrial genomes). Taxon sampling spanned 56 extant species - representing all the orders within Marchantiophyta and extant genera within Marchantiales - and eight fossil taxa. Total-evidence dating analyses support the radiation of Marchantiopsida during Late Silurian-Early Devonian (or Middle Ordovician when the outgroup is excluded) and that of Ricciaceae in the Middle Jurassic. Morphological change rate was high early in the history of the group, but it barely increased after Late Cretaceous. Disparity-through-time analyses support a fast increase in diversity until the Middle Triassic (c. 250 Ma), after which phenotypic evolution slows down considerably. Incorporating fossils in analyses challenges previous assumptions on the affinities of extinct taxa and indicates that complex thalloid liverworts radiated c. 125 Ma earlier than previously inferred.

Lipid bilayer properties potentially contributed to the evolutionary disappearance of betaine lipids in seed plants

BACKGROUND

Many organisms rely on mineral nutrients taken directly from the soil or aquatic environment, and therefore, developed mechanisms to cope with the limitation of a given essential nutrient. For example, photosynthetic cells have well-defined responses to phosphate limitation, including the replacement of cellular membrane phospholipids with non-phosphorous lipids. Under phosphate starvation, phospholipids in extraplastidial membranes are replaced by betaine lipids in microalgae. In higher plants, the synthesis of betaine lipid is lost, driving plants to other strategies to cope with phosphate starvation where they replace their phospholipids by glycolipids.

RESULTS

The aim of this work was to evaluate to what extent betaine lipids and PC lipids share physicochemical properties and could substitute for each other. By neutron diffraction experiments and dynamic molecular simulation of two synthetic lipids, the dipalmitoylphosphatidylcholine (DPPC) and the dipalmitoyl-diacylglyceryl-N,N,N-trimethylhomoserine (DP-DGTS), we found that DP-DGTS bilayers are thicker than DPPC bilayers and therefore are more rigid. Furthermore, DP-DGTS bilayers are more repulsive, especially at long range, maybe due to unexpected unscreened electrostatic contribution. Finally, DP-DGTS bilayers could coexist in the gel and fluid phases.

CONCLUSION

The different properties and hydration responses of PC and DGTS provide an explanation for the diversity of betaine lipids observed in marine organisms and for their disappearance in seed plants.

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

Reconciling fossils with phylogenies reveals the origin and macroevolutionary processes explaining the global cycad biodiversity

The determinants of biodiversity patterns can be understood using macroevolutionary analyses. The integration of fossils into phylogenies offers a deeper understanding of processes underlying biodiversity patterns in deep time. Cycadales are considered a relict of a once more diverse and globally distributed group but are restricted to low latitudes today. We still know little about their origin and geographic range evolution. Combining molecular data for extant species and leaf morphological data for extant and fossil species, we study the origin of cycad global biodiversity patterns through Bayesian total-evidence dating analyses. We assess the ancestral geographic origin and trace the historical biogeography of cycads with a time-stratified process-based model. Cycads originated in the Carboniferous on the Laurasian landmass and expanded in Gondwana in the Jurassic. Through now-vanished continental connections, Antarctica and Greenland were crucial biogeographic crossroads for cycad biogeography. Vicariance is an essential speciation mode in the deep and recent past. Their latitudinal span increased in the Jurassic and restrained toward subtropical latitudes in the Neogene in line with biogeographic inferences of high-latitude extirpations. We show the benefits of integrating fossils into phylogenies to estimate ancestral areas of origin and to study evolutionary processes explaining the global distribution of present-day relict groups.

The evolutionary history of bees in time and space

Bees are the most significant pollinators of flowering plants. This partnership began ca. 120 million years ago, but the uncertainty of how and when bees spread across the planet has greatly obscured investigations of this key mutualism. We present a novel analysis of bee biogeography using extensive new genomic and fossil data to demonstrate that bees originated in Western Gondwana (Africa and South America). Bees likely originated in the Early Cretaceous, shortly before the breakup of Western Gondwana, and the early evolution of any major bee lineage is associated with either the South American or African land masses. Subsequently, bees colonized northern continents via a complex history of vicariance and dispersal. The notable early absences from large landmasses, particularly in Australia and India, have important implications for understanding the assembly of local floras and diverse modes of pollination. How bees spread around the world from their hypothesized Southern Hemisphere origin parallels the histories of numerous flowering plant clades, providing an essential step to studying the evolution of angiosperm pollination syndromes in space and time.

Integrating Fossil Flowers into the Angiosperm Phylogeny Using Molecular and Morphological Evidence

Fossils are essential to infer past evolutionary processes. The assignment of fossils to extant clades has traditionally relied on morphological similarity and on apomorphies shared with extant taxa. The use of explicit phylogenetic analyses to establish fossil affinities has so far remained limited. In this study, we built a comprehensive framework to investigate the phylogenetic placement of 24 exceptionally preserved fossil flowers. For this, we assembled a new species-level data set of 30 floral traits for 1201 extant species that were sampled to capture the stem and crown nodes of all angiosperm families. We explored multiple analytical approaches to integrate the fossils into the phylogeny, including different phylogenetic estimation methods, topological-constrained analyses, and combining molecular and morphological data of extant and fossil species. Our results were widely consistent across approaches and showed minor differences in the support of fossils at different phylogenetic positions. The placement of some fossils agrees with previously suggested relationships, but for others, a new placement is inferred. We also identified fossils that are well supported within particular extant families, whereas others showed high phylogenetic uncertainty. Finally, we present recommendations for future analyses combining molecular and morphological evidence, regarding the selection of fossils and appropriate methodologies, and provide some perspectives on how to integrate fossils into the investigation of divergence times and the temporal evolution of morphological traits. [Angiosperms; fossil flowers; phylogenetic uncertainty; RoguePlots.].

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

Molecular timetrees using relaxed clocks and uncertain phylogenies

A common practice in molecular systematics is to infer phylogeny and then scale it to time by using a relaxed clock method and calibrations. This sequential analysis practice ignores the effect of phylogenetic uncertainty on divergence time estimates and their confidence/credibility intervals. An alternative is to infer phylogeny and times jointly to incorporate phylogenetic errors into molecular dating. We compared the performance of these two alternatives in reconstructing evolutionary timetrees using computer-simulated and empirical datasets. We found sequential and joint analyses to produce similar divergence times and phylogenetic relationships, except for some nodes in particular cases. The joint inference performed better when the phylogeny was not well resolved, situations in which the joint inference should be preferred. However, joint inference can be infeasible for large datasets because available Bayesian methods are computationally burdensome. We present an alternative approach for joint inference that combines the bag of little bootstraps, maximum likelihood, and RelTime approaches for simultaneously inferring evolutionary relationships, divergence times, and confidence intervals, incorporating phylogeny uncertainty. The new method alleviates the high computational burden imposed by Bayesian methods while achieving a similar result.

Profile of a flower: How rates of morphological evolution drive floral diversification in Ericales and angiosperms

PREMISE

Recent studies of floral disparity in the asterid order Ericales have shown that flowers vary strongly among families and that disparity is unequally distributed between the three flower modules (perianth, androecium, gynoecium). However, it remains unknown whether these patterns are driven by heterogeneous rates of morphological evolution or other factors.

METHODS

Here, we compiled a data set of 33 floral characters scored for 414 species of Ericales sampled from 346 genera and all 22 families. We conducted ancestral state reconstructions using an equal-rates Markov model for each character. We estimated rates of morphological evolution for Ericales and for a separate angiosperm-wide data set of 19 characters and 792 species, creating "rate profiles" for Ericales, angiosperms, and major angiosperm subclades. We compared morphological rates among flower modules within each data set separately and between data sets, and we compared rates among angiosperm subclades using the angiosperm data set.

RESULTS

The androecium exhibits the highest evolutionary rates across most characters, whereas most perianth and gynoecium characters evolve more slowly in both Ericales and angiosperms. Both high and low rates of morphological evolution can result in high floral disparity in Ericales. Analyses of an angiosperm-wide floral data set reveal that this pattern appears to be conserved across most major angiosperm clades.

CONCLUSIONS

Elevated rates of morphological evolution in the androecium of Ericales may explain the higher disparity reported for this floral module. Comparing rates of morphological evolution through rate profiles proves to be a powerful tool in understanding floral evolution.

An eosimiid primate of South Asian affinities in the Paleogene of Western Amazonia and the origin of New World monkeys

Recent fossil discoveries in Western Amazonia revealed that two distinct anthropoid primate clades of African origin colonized South America near the Eocene/Oligocene transition (ca. 34 Ma). Here, we describe a diminutive fossil primate from Brazilian Amazonia and suggest that, surprisingly, a third clade of anthropoids was involved in the Paleogene colonization of South America by primates. This new taxon, Ashaninkacebus simpsoni gen. et sp. nov., has strong dental affinities with Asian African stem anthropoids: the Eosimiiformes. Morphology-based phylogenetic analyses of early Old World anthropoids and extinct and extant New World monkeys (platyrrhines) support relationships of both Ashaninkacebus and Amamria (late middle Eocene, North Africa) to the South Asian Eosimiidae. Afro-Arabia, then a mega island, played the role of a biogeographic stopover between South Asia and South America for anthropoid primates and hystricognathous rodents. The earliest primates from South America bear little adaptive resemblance to later Oligocene-early Miocene platyrrhine monkeys, and the scarcity of available paleontological data precludes elucidating firmly their affinities with or within Platyrrhini. Nonetheless, these data shed light on some of their life history traits, revealing a particularly small body size and a diet consisting primarily of insects and possibly fruit, which would have increased their chances of survival on a natural floating island during this extraordinary over-water trip to South America from Africa. Divergence-time estimates between Old and New World taxa indicate that the transatlantic dispersal(s) could source in the intense flooding events associated with the late middle Eocene climatic optimum (ca. 40.5 Ma) in Western Africa.

Phylogenetic congruence, conflict and consilience between molecular and morphological data

Morphology and molecules are important data sources for estimating evolutionary relationships. Modern studies often utilise morphological and molecular partitions alongside each other in combined analyses. However, the effect of combining phenomic and genomic partitions is unclear. This is exacerbated by their size imbalance, and conflict over the efficacy of different inference methods when using morphological characters. To systematically address the effect of topological incongruence, size imbalance, and tree inference methods, we conduct a meta-analysis of 32 combined (molecular + morphology) datasets across metazoa. Our results reveal that morphological-molecular topological incongruence is pervasive: these data partitions yield very different trees, irrespective of which method is used for morphology inference. Analysis of the combined data often yields unique trees that are not sampled by either partition individually, even with the inclusion of relatively small quantities of morphological characters. Differences between morphology inference methods in terms of resolution and congruence largely relate to consensus methods. Furthermore, stepping stone Bayes factor analyses reveal that morphological and molecular partitions are not consistently combinable, i.e. data partitions are not always best explained under a single evolutionary process. In light of these results, we advise that the congruence between morphological and molecular data partitions needs to be considered in combined analyses. Nonetheless, our results reveal that, for most datasets, morphology and molecules can, and should, be combined in order to best estimate evolutionary history and reveal hidden support for novel relationships. Studies that analyse only phenomic or genomic data in isolation are unlikely to provide the full evolutionary picture.

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

Origin of adaptations to open environments and social behaviour in sabretoothed cats from the northeastern border of the Tibetan Plateau

The iconic sabretooth Homotherium is thought to have hunted cooperatively, but the origin of this behaviour and correlated morphological adaptations are largely unexplored. Here we report the most primitive species of Amphimachairodus (Amphimachairodus hezhengensis sp. nov.), a member of Machairodontini basal to Homotherium, from the Linxia Basin, northeastern border of the Tibetan Plateau (9.8-8.7 Ma). The long snout, laterally oriented and posteriorly located orbit of Amphimachairodus suggest a better ability to observe the surrounding environment, rather than targeting single prey, pointing to an adaptation to the open environment or social behaviour. A pathological forepaw of Amphimachairodus provides direct evidence of partner care. Our analyses of trait evolutionary rates support that traits correlated with killing behaviour and open environment adaptation evolved prior to other traits, suggesting that changes in hunting behaviour may be the major evolutionary driver in the early evolution of the lineage. A. hezhengensis represents one of the most important transitions in the evolution of Machairodontini, leading to adaptation in open environments and contributing to their further dispersal and radiation worldwide. This rapid morphological change is likely to be correlated with increasingly arid environments caused by the rise of the Tibetan Plateau, and competition from abundant large carnivores in this area.

Key innovations and the diversification of Hymenoptera

The order Hymenoptera (wasps, ants, sawflies, and bees) represents one of the most diverse animal lineages, but whether specific key innovations have contributed to its diversification is still unknown. We assembled the largest time-calibrated phylogeny of Hymenoptera to date and investigated the origin and possible correlation of particular morphological and behavioral innovations with diversification in the order: the wasp waist of Apocrita; the stinger of Aculeata; parasitoidism, a specialized form of carnivory; and secondary phytophagy, a reversal to plant-feeding. Here, we show that parasitoidism has been the dominant strategy since the Late Triassic in Hymenoptera, but was not an immediate driver of diversification. Instead, transitions to secondary phytophagy (from parasitoidism) had a major influence on diversification rate in Hymenoptera. Support for the stinger and the wasp waist as key innovations remains equivocal, but these traits may have laid the anatomical and behavioral foundations for adaptations more directly associated with diversification.

Diversity, Form, and Postembryonic Development of Paleozoic Insects

While Mesozoic, Paleogene, and Neogene insect faunas greatly resemble the modern one, the Paleozoic fauna provides unique insights into key innovations in insect evolution, such as the origin of wings and modifications of postembryonic development including holometaboly. Deep-divergence estimates suggest that the majority of contemporary insect orders originated in the Late Paleozoic, but these estimates reflect divergences between stem groups of each lineage rather than the later appearance of the crown groups. The fossil record shows the initial radiations of the extant hyperdiverse clades during the Early Permian, as well as the specialized fauna present before the End Permian mass extinction. This review summarizes the recent discoveries related to the documented diversity of Paleozoic hexapods, as well as current knowledge about what has actually been verified from fossil evidence as it relates to postembryonic development and the morphology of different body parts.

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.

Combining palaeontological and neontological data shows a delayed diversification burst of carcharhiniform sharks likely mediated by environmental change

Estimating deep-time species-level diversification processes remains challenging. Both the fossil record and molecular phylogenies allow the estimation of speciation and extinction rates, but each type of data may still provide an incomplete picture of diversification dynamics. Here, we combine species-level palaeontological (fossil occurrences) and neontological (molecular phylogenies) data to estimate deep-time diversity dynamics through process-based birth-death models for Carcharhiniformes, the most speciose shark order today. Despite their abundant fossil record dating back to the Middle Jurassic, only a small fraction of extant carcharhiniform species is recorded as fossils, which impedes relying only on the fossil record to study their recent diversification. Combining fossil and phylogenetic data, we recover a complex evolutionary history for carcharhiniforms, exemplified by several variations in diversification rates with an early low diversity period followed by a Cenozoic radiation. We further reveal a burst of diversification in the last 30 million years, which is partially recorded with fossil data only. We also find that reef expansion and temperature change can explain variations in speciation and extinction through time. These results pinpoint the primordial importance of these environmental variables in the evolution of marine clades. Our study also highlights the benefit of combining the fossil record with phylogenetic data to address macroevolutionary questions.

The origin of modern patterns of continental diversity in Mauritiinae palms: the Neotropical museum and the Afrotropical graveyard

While the latitudinal diversity gradient has received much attention, biodiversity and species richness also vary between continents across similar latitudes. Fossil information can be used to understand the evolutionary mechanisms that generated such variation between continents of similar latitudes. We integrated fossil data into a phylogenetic analysis of the Mauritiinae palms, whose extant diversity is restricted to the Neotropics, but extended across Africa and India during most of the Cenozoic. Mauritiinae diverged from its sister lineage Raphiinae ca 106 Ma. Using ancestral state estimation and a lineage through time analysis, we found that diversity arose globally during the late Cretaceous and Palaeocene across South America, Africa and India. The Palaeocene-Eocene transition (ca 56 Ma) marked the end of global Mauritiinae expansion, and the beginning of their decline in both Africa and India. Mauritiinae disappeared from the Indian subcontinent and Africa at the end of the Eocene and the Miocene, respectively. By contrast, Neotropical diversity steadily increased over the last 80 Myr. Taken together, our results suggest that the Neotropics functioned as a continental-scale refuge for Mauritiinae palms, where lineages survived and diversified while global climatic changes that drastically reduced rainforests led to their demise on other continents.

The complete mitochondrial genome of the woodwasp Euxiphydria potanini (Hymenoptera, Xiphydrioidea) and phylogenetic implications for symphytans

The long-necked woodwasp superfamily Xiphydrioidea belongs to the suborder Symphyta (Hymenoptera). Here we newly characterize the complete mitochondrial genome of the South Korean Euxiphydria potanini (Xiphydriidae) using next-generation sequencing: 16,500 bp long with 84.27% A + T content and 37 typical mitochondrial genes including those encoding 13 PCGs, 2 rRNAs, 22 tRNAs, and one A + T rich region. We compare the patterns of symphytan mitochondrial gene arrangement with those of an ancestral insect form and found some synapomorphic rearrangements in phylogenetic context. We use a variety of nucleotide and amino acid sequence alignments (thirteen mtPCGs and/or eight nDNAs) alongside step-by-step exclusions of long-branched taxa to elucidate the phylogenetic position of Xiphydrioidea and phylogenetic relationships among the seven symphytan superfamilies, except for Anaxyeloidea of which no mtgenome was available. The monophyly of symphytan superfamilies (with weak support for Pamphilioidea), sister-group relationship of Xiphydrioidea and Cephoidea, and Symphyta being paraphyletic to Apocrita, etc. are consistently supported by maximum likelihood and Bayesian inference trees. We also discuss the problematic phylogenetic positions of Orussoidea and Siricoidea and propose a hypothetical scenario of morphological character transition during hymenopteran evolution based on morphological key characteristics, such as the cenchrus and the wasp-waist.

New Pamphiliids with Varying Venations from Lower Cretaceous Yixian Formation of Northeast China (Hymenoptera, Pamphiliidae)

One new genus and three new species of Pamphiliidae, Dolicholyda obtusata gen. et sp. nov., Dolicholyda confluens sp. nov., and Dolicholyda angusta sp. nov. are described and illustrated. All of them were collected from the Lower Cretaceous Yixian Formation in Duolun County, Inner Mongolia, China. The new genus is established based on the following characters: body surface without punctations; forewing with pterostigma lanceolate and sclerotized around the margins; angle between 1-M and 1-Cu nearly 90°; cell 1mcu long and obviously longer than length of pterostigma. In most cases, the ap-Cu is present, and its length varied. Additionally, we revise the diagnostic characteristics of Juralydinae based on the new specimens. New findings enhance our understanding of the wing venation characteristics of fossil pamphiliids and expand the diversity of Pamphiliidae in the Mesozoic.

Was There a Cambrian Explosion on Land? The Case of Arthropod Terrestrialization

Arthropods, the most diverse form of macroscopic life in the history of the Earth, originated in the sea. Since the early Cambrian, at least ~518 million years ago, these animals have dominated the oceans of the world. By the Silurian-Devonian, the fossil record attests to arthropods becoming the first animals to colonize land, However, a growing body of molecular dating and palaeontological evidence suggests that the three major terrestrial arthropod groups (myriapods, hexapods, and arachnids), as well as vascular plants, may have invaded land as early as the Cambrian-Ordovician. These dates precede the oldest fossil evidence of those groups and suggest an unrecorded continental "Cambrian explosion" a hundred million years prior to the formation of early complex terrestrial ecosystems in the Silurian-Devonian. We review the palaeontological, phylogenomic, and molecular clock evidence pertaining to the proposed Cambrian terrestrialization of the arthropods. We argue that despite the challenges posed by incomplete preservation and the scarcity of early Palaeozoic terrestrial deposits, the discrepancy between molecular clock estimates and the fossil record is narrower than is often claimed. We discuss strategies for closing the gap between molecular clock estimates and fossil data in the evolution of early ecosystems on land.

Phylogenomic Analyses of the Tenthredinoidea Support the Familial Rank of Athaliidae (Insecta, Tenthredinoidea)

The systematic status of the genus Athalia and related genera is a perennial controversy in sawfly taxonomy. Several authors have hypothesized that the placement of Athalia within the Tenthredinidae is artificial, but no studies have focused on this topic. If the hypothesis that Athalia does not belong to Tenthredinidae can be supported, the taxonomic framework of Tenthredinoidea needs revision. We present a comprehensive phylogenomic study of Tenthredinoidae, focusing on the positions of Athalia and related genera by sampling 80 representatives mainly of the Tenthredinoidea, including Heptamelinae and Blasticotomidae. Our phylogenetic reconstructions based on nuclear genes and mitochondrial (mt) sequences support Athalia and related genera as a distinct clade sister to Tenthredinidae + (Cimbicidae + Diprionidae). A comparison of symphytan mitochondrial genomes reveals an innovative gene rearrangement pattern in Athaliidae, in which Dentathalia demonstrates a more ancestral pattern than Athalia and Hypsathalia. The lineage specificity of mt rRNA secondary structures also provides sufficient support to consider Athaliidae as a separate family. In summary, the phylogeny and genomic structural changes unanimously support the taxonomic treatment of Athaliidae as a family and the re-establishment of Dentathalia as a valid genus.

Early diversifications of angiosperms and their insect pollinators: were they unlinked?

The present-day ubiquity of angiosperm-insect pollination has led to the hypothesis that these two groups coevolved early in their evolutionary history. However, recent fossil discoveries and fossil-calibrated molecular dating analyses challenge the notion that early diversifications of angiosperms and insects were inextricably linked. In this article, we examine (i) the discrepancies between dates of emergence for angiosperms and major clades of insects; (ii) the long history of gymnosperm-insect pollination modes, which likely shaped early angiosperm-insect pollination mutualisms; and (iii) how the K-Pg (Cretaceous-Paleogene) mass extinction event was vital in propelling modern angiosperm-insect mutualisms. We posit that the early diversifications of angiosperms and their insect pollinators were largely decoupled until the end of the Cretaceous.

Total evidence time-scaled phylogenetic and biogeographic models for the evolution of sea cows (Sirenia, Afrotheria)

Molecular phylogenetic studies that have included sirenians from the genera Trichechus, Dugong, and Hydrodamalis have resolved their interrelationships but have yielded divergence age estimates that are problematically discordant. The ages of these lineage splits have profound implications for how to interpret the sirenian fossil record-including clade membership, biogeographic patterns, and correlations with Earth history events. In an effort to address these issues, here we present a total evidence phylogenetic analysis of Sirenia that includes living and fossil species and applies Bayesian tip-dating methods to estimate their interrelationships and divergence times. In addition to extant sirenians, our dataset includes 56 fossil species from 106 dated localities and numerous afrotherian outgroup taxa. Genetic, morphological, temporal, and biogeographic data are assessed simultaneously to bring all available evidence to bear on sirenian phylogeny. The resulting time-tree is then used for Bayesian geocoordinates reconstruction analysis, which models ancestral geographic areas at splits throughout the phylogeny, thereby allowing us to infer the direction and timing of dispersals. Our results suggest that Pan-Sirenia arose in North Africa during the latest Paleocene and that the Eocene evolution of stem sirenians was primarily situated in the Tethyan realm. In the late Eocene, some lineages moved into more northern European latitudes, an area that became the source region for a key trans-Atlantic dispersal towards the Caribbean and northern-adjacent west Atlantic. This event led to the phylogenetic and biogeographic founding of crown Sirenia with the Dugongidae-Trichechidae split occurring at the Eocene-Oligocene boundary (~33.9 Ma), temporally coincident with the onset of dropping global sea levels and temperatures. This region became the nexus of sirenian diversification and supported taxonomically-rich dugongid communities until the earliest Pliocene. The Dugonginae-Hydrodamalinae split occurred near Florida during the early Miocene (~21.2 Ma) and was followed by a west-bound dispersal that gave rise to the Pacific hydrodamalines. The late middle Miocene (~12.2 Ma) split of Dugong from all other dugongines also occurred near Florida and our analyses suggest that the Indo-Pacific distribution of modern dugongs is the result of a trans-Pacific dispersal. From at least the early Miocene, trichechid evolution was based entirely in South America, presumably within the Pebas Wetlands System. We infer that the eventual establishment of Amazon drainage into the South Atlantic allowed the dispersal of Trichechus out of South America no earlier than the mid-Pliocene. Our analyses provide a new temporal and biogeographic framework for understanding major events in sirenian evolution and their possible relationships to oceanographic and climatic changes. These hypotheses can be further tested with the recovery and integration of new fossil evidence.

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.

Bayesian phylogenetic inference using relaxed-clocks and the multispecies coalescent

The multispecies coalescent (MSC) model accommodates both species divergences and within-species coalescent and provides a natural framework for phylogenetic analysis of genomic data when the gene trees vary across the genome. The MSC model implemented in the program bpp assumes a molecular clock and the Jukes–Cantor model, and is suitable for analyzing genomic data from closely related species. Here we extend our implementation to more general substitution models and relaxed clocks to allow the rate to vary among species. The MSC-with-relaxed-clock model allows the estimation of species divergence times and ancestral population sizes using genomic sequences sampled from contemporary species when the strict clock assumption is violated, and provides a simulation framework for evaluating species tree estimation methods. We conducted simulations and analyzed two real datasets to evaluate the utility of the new models. We confirm that the clock-JC model is adequate for inference of shallow trees with closely related species, but it is important to account for clock violation for distant species. Our simulation suggests that there is valuable phylogenetic information in the gene-tree branch lengths even if the molecular clock assumption is seriously violated, and the relaxed-clock models implemented in bpp are able to extract such information. Our Markov chain Monte Carlo algorithms suffer from mixing problems when used for species tree estimation under the relaxed clock and we discuss possible improvements. We conclude that the new models are currently most effective for estimating population parameters such as species divergence times when the species tree is fixed.

What is the age of flowering plants?

The origin of flowering plants (angiosperms) was one of the most transformative events in the history of our planet. Despite considerable interest from multiple research fields, numerous questions remain, including the age of the group as a whole. Recent studies have reported a perplexing range of estimates for the crown-group age of angiosperms, from ~140 million years (Ma; Early Cretaceous) to 270 Ma (Permian). Both ends of the spectrum are now supported by both macroevolutionary analyses of the fossil record and fossil-calibrated molecular dating analyses. Here, we first clarify and distinguish among the three ages of angiosperms: the age of their divergence with acrogymnosperms (stem age); the age(s) of emergence of their unique, distinctive features including flowers (morphological age); and the age of the most recent common ancestor of all their living species (crown age). We then demonstrate, based on recent studies, that fossil-calibrated molecular dating estimates of the crown-group age of angiosperms have little to do with either the amount of molecular data or the number of internal fossil calibrations included. Instead, we argue that this age is almost entirely conditioned by its own prior distribution (typically a calibration density set by the user in Bayesian analyses). Lastly, we discuss which future discoveries or novel types of analyses are most likely to bring more definitive answers. In the meantime, we propose that the age of angiosperms is best described as largely unknown (140-270 Ma) and that contrasting age estimates in the literature mostly reflect conflicting prior distributions. We also suggest that future work that depends on the time scale of flowering plant diversification be designed to integrate over this vexing uncertainty.

Sexual selection promotes giraffoid head-neck evolution and ecological adaptation

The long neck of the giraffe has been held as a classic example of adaptive evolution since Darwin's time. Here we report on an unusual fossil giraffoid, Discokeryx xiezhi, from the early Miocene, which has an unusual disk-shaped headgear and the most complicated head-neck joints in known mammals. The distinctive morphology and our finite element analyses indicate an adaptation for fierce head-butting behavior. Tooth enamel isotope data suggest that D. xiezhi occupied a niche different from that of other herbivores, comparable to the characteristic high-level browsing niche of modern giraffes. The study shows that giraffoids exhibit a higher headgear diversity than other ruminants and that living in specific ecological niches may have fostered various intraspecific combat behaviors that resulted in extreme head-neck morphologies in different giraffoid lineages.

Ecological radiations of insects in the Mesozoic

The Mesozoic is a key era for the rise of the modern insect fauna. Among the most important evolutionary events in Mesozoic insects are the radiation of holometabolous insects, the origin of eusocial and parasitoid insects, diversification of pollinating insects, and development of advanced mimicry and camouflage. These events are closely associated with the diversification of insect ecological behaviors and colonization of new ecospaces. At the same time, insects had evolved more complex and closer ecological associations with various plants and animals. Mesozoic insects played a key and underappreciated ecological role in reconstructing and maintaining terrestrial ecosystems. A greater understanding of the history of insects may help to mitigate future changes in insect diversity and abundance.