Bayesian tip dating reveals heterogeneous morphological clocks in Mesozoic birds

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.

Synthetic analysis of trophic diversity and evolution in Enantiornithes with new insights from Bohaiornithidae

Enantiornithines were the dominant birds of the Mesozoic, but understanding of their diet is still tenuous. We introduce new data on the enantiornithine family Bohaiornithidae, famous for their large size and powerfully built teeth and claws. In tandem with previously published data, we comment on the breadth of enantiornithine ecology and potential patterns in which it evolved. Body mass, jaw mechanical advantage, finite element analysis of the jaw, and traditional morphometrics of the claws and skull are compared between bohaiornithids and living birds. We find bohaiornithids to be more ecologically diverse than any other enantiornithine family: Bohaiornis and Parabohaiornis are similar to living plant-eating birds; Longusunguis resembles raptorial carnivores; Zhouornis is similar to both fruit-eating birds and generalist feeders; and Shenqiornis and Sulcavis plausibly ate fish, plants, or a mix of both. We predict the ancestral enantiornithine bird to have been a generalist which ate a wide variety of foods. However, more quantitative data from across the enantiornithine tree is needed to refine this prediction. By the Early Cretaceous, enantiornithine birds had diversified into a variety of ecological niches like crown birds after the K-Pg extinction, adding to the evidence that traits unique to crown birds cannot completely explain their ecological success.

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.

Low morphological disparity and decelerated rate of limb size evolution close to the origin of birds

The origin of birds from theropod dinosaurs involves many changes in musculoskeletal anatomy and epidermal structures, including multiple instances of convergence and homology-related traits that contribute to the refinement of flight capability. Changes in limb sizes and proportions are important for locomotion (for example, the forelimb for bird flight); thus, understanding these patterns is central to investigating the transition from terrestrial to volant theropods. Here we analyse the patterns of morphological disparity and the evolutionary rate of appendicular limbs along avialan stem lineages using phylogenetic comparative approaches. Contrary to the traditional wisdom that an evolutionary innovation like flight would promote and accelerate evolvability, our results show a shift to low disparity and decelerated rate near the origin of avialans that is largely ascribed to the evolutionarily constrained forelimb. These results suggest that natural selection shaped patterns of limb evolution close to the origin of avialans in a way that may reflect the winged forelimb 'blueprint' associated with powered flight.

Complex macroevolution of pterosaurs

Pterosaurs, the earliest flying tetrapods, are the subject of some recent quantitative macroevolutionary analyses from different perspectives.^1-2 Here, we use an integrative approach involving newly assembled phylogenetic and body size datasets, net diversification rates, morphological rates, and morphological disparity to gain a holistic understanding of the pterosaur macroevolution. The first two parameters are important in quantitative analyses of macroevolution, but they have been rarely used in previous pterosaur studies.^{1,3,4,2,5,6,7,8,9,10,11,12} Our study reveals an ∼115-Ma period-from Early Triassic to Early Cretaceous-of multi-wave increasing net diversification rates and disparity, as well as high morphological rates, followed by an ∼65-Ma period-from Early Cretaceous to the end of the Cretaceous-of mostly negative net diversification rates, decreasing disparity, and relatively low morphological rates in pterosaur evolution. Our study demonstrates the following: (1) body size plays an important role in pterosaur lineage diversification during nearly their whole evolutionary history, and the evolution of locomotion, trophic, and ornamental structures also plays a role in different periods; (2) birds, the other major flying tetrapod group at the time, might have affected pterosaur macroevolution for ∼100 Ma; and (3) different mass extinction events might have affected pterosaur evolution differently. Particularly, the revealed decline in pterosaur biodiversity during the Middle and Late Cretaceous periods provides further support for the possible presence of a biodiversity decline of large-sized terrestrial amniotes starting in the mid-Cretaceous,^13,14 which may have been caused by multiple factors including a global land area decrease during these periods.

Evolutionary origins of the prolonged extant squamate radiation

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

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

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

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

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

The Jurassic rise of squamates as supported by lepidosaur disparity and evolutionary rates

The squamates (lizards, snakes, and relatives) today comprise more than 10,000 species, and yet their sister group, the Rhynchocephalia, is represented by a single species today, the tuatara. The explosion in squamate diversity has been tracked back to the Cretaceous Terrestrial Revolution, 100 million years ago (Ma), the time when flowering plants began their takeover of terrestrial ecosystems, associated with diversification of coevolving insects and insect-eating predators such as lizards, birds, and mammals. Squamates arose much earlier, but their long pre-Cretaceous history of some 150 million years (Myr) is documented by sparse fossils. Here, we provide evidence for an initial radiation of squamate morphology in the Middle and Late Jurassic (174-145 Ma), and show that they established their key ecological roles much earlier than had been assumed, and they have not changed them much since.

Dasycerine rove beetles: Cretaceous diversification, phylogeny and historical biogeography (Coleoptera: Staphylinidae: Dasycerinae)

Within the hyperdiverse beetle family Staphylinidae, Dasycerinae is one of the smallest and most cryptic subfamilies, comprising a sole extant genus characterized by a latridiid beetle-like body form. Little has been known about their early diversification, character evolution, phylogeny and historical biogeography because of limited fossil material and lack of a phylogeny integrating extant and extinct representatives. Here we report an unexpectedly diverse dasycerine fauna from the mid-Cretaceous of northern Myanmar, including a new genus and four new species. To reconstruct the early evolutionary history of Dasycerinae, we present a phylogenetic framework of the subfamily based on a dataset integrating all extant and extinct taxa using parsimony, maximum-likelihood and Bayesian methods. †Cedasyrus gen. n., characterized by distinct sexual dimorphism in antennal and elytral lengths, is recovered as the basal-most lineage, sister to the remaining two extinct genera and all living Dasycerus species. †Vetudasycerus is recovered as sister to †Protodasycerus + Dasycerus. Among all extinct taxa, †Protodasycerus bears distinctly longer elytra, and appears to represent a transitional form from †Vetudasycerus to Dasycerus. Phylogenetic inferences and ancestral distribution reconstruction support an "Out-of-Orient" model for Dasycerinae. Either the Bering- or North Atlantic Land Bridge may have served as a passageway for dasycerine dispersal between Eurasian and North American continents. An elevation-reconstruction analysis indicated that the ancestor of the extant Dasycerus probably lived at a high altitude and stayed at this elevation through the end of the Miocene. We propose that the extinction of dasycerine ancestors living on the Tethyan islands at low altitude was likely caused by sea-level rise and climatic warming during the Late Cretaceous. The high-altitude areas might have played the role of refugia that harboured subalpine derivatives which eventually gave rise to the extant Dasycerus.

Massive cranium from Harbin in northeastern China establishes a new Middle Pleistocene human lineage

It has recently become clear that several human lineages coexisted with Homo sapiens during the late Middle and Late Pleistocene. Here, we report an archaic human fossil that throws new light on debates concerning the diversification of the Homo genus and the origin of H. sapiens. The fossil was recovered in Harbin city in northeastern China, with a minimum uranium-series age of 146 ka. This cranium is one of the best preserved Middle Pleistocene human fossils. Its massive size, with a large cranial capacity (∼1,420 mL) falling in the range of modern humans, is combined with a mosaic of primitive and derived characters. It differs from all the other named Homo species by presenting a combination of features, such as long and low cranial vault, a wide and low face, large and almost square orbits, gently curved but massively developed supraorbital torus, flat and low cheekbones with a shallow canine fossa, and a shallow palate with thick alveolar bone supporting very large molars. The excellent preservation of the Harbin cranium advances our understanding of several less-complete late Middle Pleistocene fossils from China, which have been interpreted as local evolutionary intermediates between the earlier species Homo erectus and later H. sapiens. Phylogenetic analyses based on parsimony criteria and Bayesian tip-dating suggest that the Harbin cranium and some other Middle Pleistocene human fossils from China, such as those from Dali and Xiahe, form a third East Asian lineage, which is a part of the sister group of the H. sapiens lineage. Our analyses of such morphologically distinctive archaic human lineages from Asia, Europe, and Africa suggest that the diversification of the Homo genus may have had a much deeper timescale than previously presumed. Sympatric isolation of small populations combined with stochastic long-distance dispersals is the best fitting biogeographical model for interpreting the evolution of the Homo genus.

•

More than 100,000 years ago, several human species coexisted in Asia, Europe, and Africa

•

A completely preserved fossil human cranium discovered in the Harbin area provides critical evidence for understanding the evolution of humans and the origin of our species

•

The Harbin cranium has a large cranial capacity (∼1,420 mL) falling in the range of modern humans, but is combined with a mosaic of primitive and derived characters

•

Our comprehensive phylogenetic analyses suggest that the Harbin cranium represents a new sister lineage for Homo sapiens

•

A multi-directional “shuttle dispersal model” is more likely to explain the complex phylogenetic connections among African and Eurasian Homo species/populations

Sustained high rates of morphological evolution during the rise of tetrapods

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

Deep time diversity and the early radiations of birds

Reconstructing the history of biodiversity has been hindered by often-separate analyses of stem and crown groups of the clades in question that are not easily understood within the same unified evolutionary framework. Here, we investigate the evolutionary history of birds by analyzing three supertrees that combine published phylogenies of both stem and crown birds. Our analyses reveal three distinct large-scale increases in the diversification rate across bird evolutionary history. The first increase, which began between 160 and 170 Ma and reached its peak between 130 and 135 Ma, corresponds to an accelerated morphological evolutionary rate associated with the locomotory systems among early stem birds. This radiation resulted in morphospace occupation that is larger and different from their close dinosaurian relatives, demonstrating the occurrence of a radiation among early stem birds. The second increase, which started ∼90 Ma and reached its peak between 65 and 55 Ma, is associated with rapid evolution of the cranial skeleton among early crown birds, driven differently from the first radiation. The third increase, which occurred after ∼40 to 45 Ma, has yet to be supported by quantitative morphological data but gains some support from the fossil record. Our analyses indicate that the bird biodiversity evolution was influenced mainly by long-term climatic changes and also by major paleobiological events such as the Cretaceous-Paleogene (K-Pg) extinction.

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

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

The patterns and modes of the evolution of disparity in Mesozoic birds

The origin of birds from non-avian theropod dinosaurs is one of the greatest transitions in evolution. Shortly after diverging from other theropods in the Late Jurassic, Mesozoic birds diversified into two major clades-the Enantiornithes and Ornithuromorpha-acquiring many features previously considered unique to the crown group along the way. Here, we present a comparative phylogenetic study of the patterns and modes of Mesozoic bird skeletal morphology and limb proportions. Our results show that the major Mesozoic avian groups are distinctive in discrete character space, but constrained in a morphospace defined by limb proportions. The Enantiornithines, despite being the most speciose group of Mesozoic birds, are much less morphologically disparate than their sister clade, the Ornithuromorpha-the clade that gave rise to living birds, showing disparity and diversity were decoupled in avian history. This relatively low disparity suggests that diversification of enantiornithines was characterized in exhausting fine morphologies, whereas ornithuromorphs continuously explored a broader array of morphologies and ecological opportunities. We suggest this clade-specific evolutionary versatility contributed to their sole survival of the end-Cretaceous mass extinction.

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.

Decelerated dinosaur skull evolution with the origin of birds

The evolutionary radiation of birds has produced incredible morphological variation, including a huge range of skull form and function. Investigating how this variation arose with respect to non-avian dinosaurs is key to understanding how birds achieved their remarkable success after the Cretaceous-Paleogene extinction event. Using a high-dimensional geometric morphometric approach, we quantified the shape of the skull in unprecedented detail across 354 extant and 37 extinct avian and non-avian dinosaurs. Comparative analyses reveal fundamental differences in how skull shape evolved in birds and non-avian dinosaurs. We find that the overall skull shape evolved faster in non-avian dinosaurs than in birds across all regions of the cranium. In birds, the anterior rostrum is the most rapidly evolving skull region, whereas more posterior regions-such as the parietal, squamosal, and quadrate-exhibited high rates in non-avian dinosaurs. These fast-evolving elements in dinosaurs are strongly associated with feeding biomechanics, forming the jaw joint and supporting the jaw adductor muscles. Rapid pulses of skull evolution coincide with changes to food acquisition strategies and diets, as well as the proliferation of bony skull ornaments. In contrast to the appendicular skeleton, which has been shown to evolve more rapidly in birds, avian cranial morphology is characterised by a striking deceleration in morphological evolution relative to non-avian dinosaurs. These results may be due to the reorganisation of skull structure in birds-including loss of a separate postorbital bone in adults and the emergence of new trade-offs with development and neurosensory demands. Taken together, the remarkable cranial shape diversity in birds was not a product of accelerated evolution from their non-avian relatives, despite their frequent portrayal as an icon of adaptive radiations.

Exact Distribution of Divergence Times from Fossil Ages and Tree Topologies

Being given a phylogenetic tree of both extant and extinct taxa in which the fossil ages are the only temporal information (namely, in which divergence times are considered unknown), we provide a method to compute the exact probability distribution of any divergence time of the tree with regard to any speciation (cladogenesis), extinction, and fossilization rates under the Fossilized Birth–Death model. We use this new method to obtain a probability distribution for the age of Amniota (the synapsid/sauropsid or bird/mammal divergence), one of the most-frequently used dating constraints. Our results suggest an older age (between about 322 and 340 Ma) than has been assumed by most studies that have used this constraint (which typically assumed a best estimate around 310–315 Ma) and provide, for the first time, a method to compute the shape of the probability density for this divergence time. [Divergence times; fossil ages; fossilized birth–death model; probability distribution.]