Bayesian estimation of speciation and extinction from incomplete fossil occurrence data

PhyloJunction: A Computational Framework for Simulating, Developing, and Teaching Evolutionary Models

We introduce PhyloJunction, a computational framework designed to facilitate the prototyping, testing, and characterization of evolutionary models. PhyloJunction is distributed as an open-source Python library that can be used to implement a variety of models, thanks to its flexible graphical modeling architecture and dedicated model specification language. Model design and use are exposed to users via command-line and graphical interfaces, which integrate the steps of simulating, summarizing, and visualizing data. This article describes the features of PhyloJunction-which include, but are not limited to, a general implementation of a popular family of phylogenetic diversification models-and, moving forward, how it may be expanded to not only include new models, but to also become a platform for conducting and teaching statistical learning.

Trait-mediated speciation and human-driven extinctions in proboscideans revealed by unsupervised Bayesian neural networks

Species life-history traits, paleoenvironment, and biotic interactions likely influence speciation and extinction rates, affecting species richness over time. Birth-death models inferring the impact of these factors typically assume monotonic relationships between single predictors and rates, limiting our ability to assess more complex effects and their relative importance and interaction. We introduce a Bayesian birth-death model using unsupervised neural networks to explore multifactorial and nonlinear effects on speciation and extinction rates using fossil data. It infers lineage- and time-specific rates and disentangles predictor effects and importance through explainable artificial intelligence techniques. Analysis of the proboscidean fossil record revealed speciation rates shaped by dietary flexibility and biogeographic events. The emergence of modern humans escalated extinction rates, causing recent diversity decline, while regional climate had a lesser impact. Our model paves the way for an improved understanding of the intricate dynamics shaping clade diversification.

Late Cretaceous ammonoids show that drivers of diversification are regionally heterogeneous

Palaeontologists have long sought to explain the diversification of individual clades to whole biotas at global scales. Advances in our understanding of the spatial distribution of the fossil record through geological time, however, has demonstrated that global trends in biodiversity were a mosaic of regionally heterogeneous diversification processes. Drivers of diversification must presumably have also displayed regional variation to produce the spatial disparities observed in past taxonomic richness. Here, we analyse the fossil record of ammonoids, pelagic shelled cephalopods, through the Late Cretaceous, characterised by some palaeontologists as an interval of biotic decline prior to their total extinction at the Cretaceous-Paleogene boundary. We regionally subdivide this record to eliminate the impacts of spatial sampling biases and infer regional origination and extinction rates corrected for temporal sampling biases using Bayesian methods. We then model these rates using biotic and abiotic drivers commonly inferred to influence diversification. Ammonoid diversification dynamics and responses to this common set of diversity drivers were regionally heterogeneous, do not support ecological decline, and demonstrate that their global diversification signal is influenced by spatial disparities in sampling effort. These results call into question the feasibility of seeking drivers of diversity at global scales in the fossil record.

A quantitative test of the "Ecomorphotype Hypothesis" for fossil true seals (Family Phocidae)

The fossil record of true seals (Family Phocidae) is mostly made up of isolated bones, some of which are type specimens. Previous studies have sought to increase referral of non-overlapping and unrelated fossils to these taxa using the 'Ecomorphotype Hypothesis', which stipulates that certain differences in morphology between taxa represent adaptations to differing ecology. On this basis, bulk fossil material could be lumped to a specific ecomorphotype, and then referred to species in that ecomorphotype, even if they are different bones. This qualitative and subjective method has been used often to expand the taxonomy of fossil phocids, but has never been quantitatively tested. We test the proposed ecomorphotypes using morphometric analysis of fossil and extant northern true seal limb bones, specifically principal components analysis and discriminant function analysis. A large amount of morphological overlap between ecomorphotypes, and poor discrimination between them, suggests that the 'Ecomorphotype Hypothesis' is not a valid approach. Further, the analysis failed to assign fossils to ecomorphotypes designated in previous studies, with some fossils from the same taxa being designated as different ecomorphotypes. The failure of this approach suggests that all fossils referred using this method should be considered to have unknown taxonomic status. In light of this, and previous findings that phocid limb bones have limited utility as type specimens, we revise the status of named fossil phocid species. We conclude that the majority of named fossil phocid taxa should be considered nomina dubia.

Testing extinction events and temporal shifts in diversification and fossilization rates through the skyline Fossilized Birth-Death (FBD) model: The example of some mid-Permian synapsid extinctions

In the last decade, the Fossilized Birth-Death (FBD) process has yielded interesting clues about the evolution of biodiversity through time. To facilitate such studies, we extend our method to compute the probability density of phylogenetic trees of extant and extinct taxa in which the only temporal information is provided by the fossil ages (i.e. without the divergence times) in order to deal with the piecewise constant FBD process, known as the "skyline FBD", which allows rates to change between pre-defined time intervals, as well as modelling extinction events at the bounds of these intervals. We develop approaches based on this method to assess hypotheses about the diversification process and to answer questions such as "Does a mass extinction occur at this time?" or "Is there a change in the fossilization rate between two given periods?". Our software can also yield Bayesian and maximum-likelihood estimates of the parameters of the skyline FBD model under various constraints. These approaches are applied to a simulated dataset in order to test their ability to answer the questions above. Finally, we study an updated dataset of Permo-Carboniferous synapsids to get additional insights into the dynamics of biodiversity change in three clades (Ophiacodontidae, Edaphosauridae and Sphenacodontidae) in the Pennsylvanian (Late Carboniferous) and Cisuralian (Early Permian), and to assess support for end-Sakmarian (or Artinskian) and end-Cisuralian mass extinction events discussed in previous studies.

Imbalanced speciation pulses sustain the radiation of mammals

The evolutionary histories of major clades, including mammals, often comprise changes in their diversification dynamics, but how these changes occur remains debated. We combined comprehensive phylogenetic and fossil information in a new "birth-death diffusion" model that provides a detailed characterization of variation in diversification rates in mammals. We found an early rising and sustained diversification scenario, wherein speciation rates increased before and during the Cretaceous-Paleogene (K-Pg) boundary. The K-Pg mass extinction event filtered out more slowly speciating lineages and was followed by a subsequent slowing in speciation rates rather than rebounds. These dynamics arose from an imbalanced speciation process, with separate lineages giving rise to many, less speciation-prone descendants. Diversity seems to have been brought about by these isolated, fast-speciating lineages, rather than by a few punctuated innovations.

The Angiosperm Terrestrial Revolution buffered ants against extinction

With ~14,000 extant species, ants are ubiquitous and of tremendous ecological importance. They have undergone remarkable diversification throughout their evolutionary history. However, the drivers of their diversity dynamics are not well quantified or understood. Previous phylogenetic analyses have suggested patterns of diversity dynamics associated with the Angiosperm Terrestrial Revolution (ATR), but these studies have overlooked valuable information from the fossil record. To address this gap, we conducted a comprehensive analysis using a large dataset that includes both the ant fossil record (~24,000 individual occurrences) and neontological data (~14,000 occurrences), and tested four hypotheses proposed for ant diversification: co-diversification, competitive extinction, hyper-specialization, and buffered extinction. Taking into account biases in the fossil record, we found three distinct diversification periods (the latest Cretaceous, Eocene, and Oligo-Miocene) and one extinction period (Late Cretaceous). The competitive extinction hypothesis between stem and crown ants is not supported. Instead, we found support for the co-diversification, buffered extinction, and hyper-specialization hypotheses. The environmental changes of the ATR, mediated by the angiosperm radiation, likely played a critical role in buffering ants against extinction and favoring their diversification by providing new ecological niches, such as forest litter and arboreal nesting sites, and additional resources. We also hypothesize that the decline and extinction of stem ants during the Late Cretaceous was due to their hyper-specialized morphology, which limited their ability to expand their dietary niche in changing environments. This study highlights the importance of a holistic approach when studying the interplay between past environments and the evolutionary trajectories of organisms.

The angiosperm radiation played a dual role in the diversification of insects and insect pollinators

Interactions with angiosperms have been hypothesised to play a crucial role in driving diversification among insects, with a particular emphasis on pollinator insects. However, support for coevolutionary diversification in insect-plant interactions is weak. Macroevolutionary studies of insect and plant diversities support the hypothesis that angiosperms diversified after a peak in insect diversity in the Early Cretaceous. Here, we used the family-level fossil record of insects as a whole, and insect pollinator families in particular, to estimate diversification rates and the role of angiosperms on insect macroevolutionary history using a Bayesian process-based approach. We found that angiosperms played a dual role that changed through time, mitigating insect extinction in the Cretaceous and promoting insect origination in the Cenozoic, which is also recovered for insect pollinator families only. Although insects pollinated gymnosperms before the angiosperm radiation, a radiation of new pollinator lineages began as angiosperm lineages increased, particularly significant after 50 Ma. We also found that global temperature, increases in insect diversity, and spore plants were strongly correlated with origination and extinction rates, suggesting that multiple drivers influenced insect diversification and arguing for the investigation of different explanatory variables in further studies.

Drivers of the artiodactyl turnover in insular western Europe at the Eocene-Oligocene Transition

Simultaneously investigating the effects of abiotic and biotic factors on diversity dynamics is essential to understand the evolutionary history of clades. The Grande Coupure corresponds to a major faunal turnover at the Eocene-Oligocene transition (EOT) (~34.1 to 33.55 Mya) and is defined in western Europe as an extinction of insular European mammals coupled with the arrival of crown clades from Asia. Here, we focused on the species-rich group of endemic European artiodactyls to determine the drivers of the Grande Coupure during the major environmental disruptions at the EOT. Using Bayesian birth-death models, we analyzed an original high-resolution fossil dataset (90 species, >2,100 occurrences) from southwestern France (Quercy area) and estimated the regional diversification and diversity dynamics of endemic and immigrant artiodactyls. We show that the endemic artiodactyl radiation was mainly related to the Eocene tropical conditions, combined with biotic controls on speciation and clade-related diversity dependence. We further highlight that the major environmental changes at the transition (77% of species became extinct) and the concurrent increase in seasonality in Europe during the Oligocene were likely the main drivers of their decline. Surprisingly, our results do not support the widely-held hypothesis of active competition between endemic and immigrant artiodactyls but rather suggest a passive or opportunistic replacement by immigrants, which is further supported by morphological clustering of specific ecological traits across the Eocene-Oligocene transition. Our analyses provide insights into the evolutionary and ecological processes driving the diversification and decline of mammalian clades during a major biological and climatic crisis.

PhyloJunction: a computational framework for simulating, developing, and teaching evolutionary models

We introduce PhyloJunction, a computational framework designed to facilitate the prototyping, testing, and characterization of evolutionary models. PhyloJunction is distributed as an open-source Python library that can be used to implement a variety of models, through its flexible graphical modeling architecture and dedicated model specification language. Model design and use are exposed to users via command-line and graphical interfaces, which integrate the steps of simulating, summarizing, and visualizing data. This paper describes the features of PhyloJunction - which include, but are not limited to, a general implementation of a popular family of phylogenetic diversification models - and, moving forward, how it may be expanded to not only include new models, but to also become a platform for conducting and teaching statistical learning.

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.

Fluctuation in the diversity of mayflies (Insecta, Ephemerida) as documented in the fossil record

Due to their aquatic larvae, the evolution of mayflies is intricately tied to environmental changes affecting lakes and rivers. Despite a rich fossil record, little is known about the factors shaping the pattern of diversification of mayflies in deep time. We assemble an unprecedented dataset encompassing all fossil occurrences of mayflies and perform a Bayesian analysis to identify periods of increased origination or extinction. We provide strong evidence for a major extinction of mayflies in the mid-Cretaceous. This extinction and subsequent faunal turnover were probably connected with the rise of angiosperms. Their dominance caused increased nutrient input and changed the chemistry of the freshwater environments, a trend detrimental mainly to lacustrine insects. Mayflies underwent a habitat shift from hypotrophic lakes to running waters, where most of their diversity has been concentrated from the Late Cretaceous to the present.

Bayesian analyses indicate bivalves did not drive the downfall of brachiopods following the Permian-Triassic mass extinction

Certain times of major biotic replacement have often been interpreted as broadly competitive, mediated by innovation in the succeeding clades. A classic example was the switch from brachiopods to bivalves as major seabed organisms following the Permian-Triassic mass extinction (PTME), ~252 million years ago. This was attributed to competitive exclusion of brachiopods by the better adapted bivalves or simply to the fact that brachiopods had been hit especially hard by the PTME. The brachiopod-bivalve switch is emblematic of the global turnover of marine faunas from Palaeozoic-type to Modern-type triggered by the PTME. Here, using Bayesian analyses, we find that unexpectedly the two clades displayed similar large-scale trends of diversification before the Jurassic. Insight from a multivariate birth-death model shows that the extinction of major brachiopod clades during the PTME set the stage for the brachiopod-bivalve switch, with differential responses to high ocean temperatures post-extinction further facilitating their displacement by bivalves. Our study strengthens evidence that brachiopods and bivalves were not competitors over macroevolutionary time scales, with extinction events and environmental stresses shaping their divergent fates.

Diversification Models Conflate Likelihood and Prior, and Cannot be Compared Using Conventional Model-Comparison Tools

Time-calibrated phylogenetic trees are a tremendously powerful tool for studying evolutionary, ecological, and epidemiological phenomena. Such trees are predominantly inferred in a Bayesian framework, with the phylogeny itself treated as a parameter with a prior distribution (a "tree prior"). However, we show that the tree "parameter" consists, in part, of data, in the form of taxon samples. Treating the tree as a parameter fails to account for these data and compromises our ability to compare among models using standard techniques (e.g., marginal likelihoods estimated using path-sampling and stepping-stone sampling algorithms). Since accuracy of the inferred phylogeny strongly depends on how well the tree prior approximates the true diversification process that gave rise to the tree, the inability to accurately compare competing tree priors has broad implications for applications based on time-calibrated trees. We outline potential remedies to this problem, and provide guidance for researchers interested in assessing the fit of tree models. [Bayes factors; Bayesian model comparison; birth-death models; divergence-time estimation; lineage diversification].

Global impact and selectivity of the Cretaceous-Paleogene mass extinction among sharks, skates, and rays

The Cretaceous-Paleogene event was the last mass extinction event, yet its impact and long-term effects on species-level marine vertebrate diversity remain largely uncharacterized. We quantified elasmobranch (sharks, skates, and rays) speciation, extinction, and ecological change resulting from the end-Cretaceous event using >3200 fossil occurrences and 675 species spanning the Late Cretaceous-Paleocene interval at global scale. Elasmobranchs declined by >62% at the Cretaceous-Paleogene boundary and did not fully recover in the Paleocene. The end-Cretaceous event triggered a heterogeneous pattern of extinction, with rays and durophagous species reaching the highest levels of extinction (>72%) and sharks and nondurophagous species being less affected. Taxa with large geographic ranges and/or those restricted to high-latitude settings show higher survival. The Cretaceous-Paleogene event drastically altered the evolutionary history of marine ecosystems.

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.

Multiple drivers and lineage-specific insect extinctions during the Permo-Triassic

The Permo-Triassic interval encompasses three extinction events including the most dramatic biological crisis of the Phanerozoic, the latest Permian mass extinction. However, their drivers and outcomes are poorly quantified and understood for terrestrial invertebrates, which we assess here for insects. We find a pattern with three extinctions: the Roadian/Wordian (≈266.9 Ma; extinction of 64.5% insect genera), the Permian/Triassic (≈252 Ma; extinction of 82.6% insect genera), and the Ladinian/Carnian boundaries (≈237 Ma; extinction of 74.8% insect genera). We also unveil a heterogeneous effect of these extinction events across the major insect clades. Because extinction events have impacted Permo-Triassic ecosystems, we investigate the influence of abiotic and biotic factors on insect diversification dynamics and find that changes in floral assemblages are likely the strongest drivers of insects' responses throughout the Permo-Triassic. We also assess the effect of diversity dependence between three insect guilds; an effect ubiquitously found in current ecosystems. We find that herbivores held a central position in the Permo-Triassic interaction network. Our study reveals high levels of insect extinction that profoundly shaped the evolutionary history of the most diverse non-microbial lineage.

Mapping species diversification metrics in macroecology: Prospects and challenges

The intersection of macroecology and macroevolution is one of today’s most active research in biology. In the last decade, we have witnessed a steady increment of macroecological studies that use metrics attempting to capture macroevolutionary processes to explain present-day biodiversity patterns. Evolutionary explanations of current species richness gradients are fundamental for understanding how diversity accumulates in a region. Although multiple hypotheses have been proposed to explain the patterns we observe in nature, it is well-known that the present-day diversity patterns result from speciation, extinction, colonization from nearby areas, or a combination of these macroevolutionary processes. Whether these metrics capture macroevolutionary processes across space is unknown. Some tip-rate metrics calculated directly from a phylogenetic tree (e.g., mean root distance -MRD-; mean diversification rate -mDR-) seem to return very similar geographical patterns regardless of how they are estimated (e.g., using branch lengths explicitly or not). Model-based tip-rate metrics —those estimated using macroevolutionary mixtures, e.g., the BAMM approach— seem to provide better net diversification estimates than only speciation rates. We argue that the lack of appropriate estimates of extinction and dispersal rates in phylogenetic trees may strongly limit our inferences about how species richness gradients have emerged at spatial and temporal scales. Here, we present a literature review about this topic and empirical comparisons between select taxa with several of these metrics. We implemented a simple null model approach to evaluate whether mapping of these metrics deviates from a random sampling process. We show that phylogenetic metrics by themselves are relatively poor at capturing speciation, extinction, and dispersal processes across geographical gradients. Furthermore, we provide evidence of how parametric biogeographic methods can improve our inference of past events and, therefore, our conclusions about the evolutionary processes driving biodiversity patterns. We recommend that further studies include several approaches simultaneously (e.g., spatial diversification modeling, parametric biogeographic methods, simulations) to disentangle the relative role of speciation, extinction, and dispersal in the generation and maintenance of species richness gradients at regional and global scales.

The role of paleontological data in bryophyte systematics

Systematics reconstructs tempo and mode in biological evolution by resolving the phylogenetic fabric of biodiversity. The staggering duration and complexity of evolution, coupled with loss of information (extinction), render exhaustive reconstruction of the evolutionary history of life unattainable. Instead, we sample its products-phenotypes and genotypes-to generate phylogenetic hypotheses, which we sequentially reassess and update against new data. Current consensus in evolutionary biology emphasizes fossil integration in total-evidence analyses, requiring in-depth understanding of fossils-age, phenotypes, and systematic affinities-and a detailed morphological framework uniting fossil and extant taxa. Bryophytes present a special case: deep evolutionary history but sparse fossil record and phenotypic diversity encompassing small dimensional scales. We review how these peculiarities shape fossil inclusion in bryophyte systematics. Paucity of the bryophyte fossil record, driven primarily by phenotypic (small plant size) and ecological constraints (patchy substrate-hugging populations), and incomplete exploration, results in many morphologically isolated, taxonomically ambiguous fossil taxa. Nevertheless, instances of exquisite preservation and pioneering studies demonstrate the feasibility of including bryophyte fossils in evolutionary inference. Further progress will arise from developing extensive morphological matrices for bryophytes, continued exploration of the fossil record, re-evaluation of previously described fossils, and training specialists in identification and characterization of bryophyte fossils, and in bryophyte morphology.

Diversity dynamics of microfossils from the Cretaceous to the Neogene show mixed responses to events

Microfossils have a ubiquitous and well-studied fossil record with temporally and spatially fluctuating diversity, but how this arises and how major events affect speciation and extinction is uncertain. We present one of the first applications of PyRate to a micropalaeontological global occurrence dataset, reconstructing diversification rates within a Bayesian framework from the Mesozoic to the Neogene in four microfossil groups: planktic foraminiferans, calcareous nannofossils, radiolarians and diatoms. Calcareous and siliceous groups demonstrate opposed but inconsistent responses in diversification. Radiolarian origination increases from c. 104 Ma, maintaining high rates into the Cenozoic. Calcareous microfossil diversification rates significantly declines across the Cretaceous-Palaeogene boundary, while rates in siliceous microfossil groups remain stable until the Paleocene-Eocene transition. Diversification rates in the Cenozoic are largely stable in calcareous groups, whereas the Palaeogene is a turbulent time for diatoms. Diversification fluctuations are driven by climate change and fluctuations in sea surface temperatures, leading to different responses in the groups generating calcareous or siliceous microfossils. Extinctions are apparently induced by changes in anoxia, acidification and stratification; speciation tends to be associated with upwelling, productivity and ocean circulation. These results invite further micropalaeontological quantitative analysis and study of the effects of major transitions in the fossil record. Despite extensive occurrence data, regional diversification events were not recovered; neither were some global events. These unexpected results show the need to consider multiple spatiotemporal levels of diversity and diversification analyses and imply that occurrence datasets of different clades may be more appropriate for testing some hypotheses than others.

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.

Studying speciation and extinction dynamics from phylogenies: addressing identifiability issues

A lot of what we know about past speciation and extinction dynamics is based on statistically fitting birth-death processes to phylogenies of extant species. Despite their wide use, the reliability of these tools is regularly questioned. It was recently demonstrated that vast 'congruent' sets of alternative diversification histories cannot be distinguished (i.e., are not identifiable) using extant phylogenies alone, reanimating the debate about the limits of phylogenetic diversification analysis. Here, we summarize what we know about the identifiability of the birth-death process and how identifiability issues can be addressed. We conclude that extant phylogenies, when combined with appropriate prior hypotheses and regularization techniques, can still tell us a lot about past diversification dynamics.

Global diversity dynamics in the fossil record are regionally heterogeneous

Global diversity patterns in the fossil record comprise a mosaic of regional trends, underpinned by spatially non-random drivers and distorted by variation in sampling intensity through time and across space. Sampling-corrected diversity estimates from spatially-standardised fossil datasets retain their regional biogeographic nuances and avoid these biases, yet diversity-through-time arises from the interplay of origination and extinction, the processes that shape macroevolutionary history. Here we present a subsampling algorithm to eliminate spatial sampling bias, coupled with advanced probabilistic methods for estimating origination and extinction rates and a Bayesian method for estimating sampling-corrected diversity. We then re-examine the Late Permian to Early Jurassic marine fossil record, an interval spanning several global biotic upheavals that shaped the origins of the modern marine biosphere. We find that origination and extinction rates are regionally heterogenous even during events that manifested globally, highlighting the need for spatially explicit views of macroevolutionary processes through geological time.

Global diversity trends in the fossil record vary regionally through time and space, affecting our ability to interpret macroevolutionary history. Here, the authors propose a method to eliminate spatial sampling bias, estimate origination and extinction rates, and generate diversity estimates, applying this method to the Late Permian to Early Jurassic marine fossil record.

Fast-growing species shape the evolution of reef corals

Ecological interactions are ubiquitous on tropical coral reefs, where sessile organisms coexist in limited space. Within these high-diversity systems, reef-building scleractinian corals form an intricate interaction network. The role of biotic interactions among reef corals is well established on ecological timescales. However, its potential effect on macroevolutionary patterns remains unclear. By analysing the rich fossil record of Scleractinia, we show that reef coral biodiversity experienced marked evolutionary rate shifts in the last 3 million years, possibly driven by biotic interactions. Our models suggest that there was an overwhelming effect of staghorn corals (family Acroporidae) on the fossil diversity trajectories of other coral groups. Staghorn corals showed an unparalleled spike in diversification during the Pleistocene. But surprisingly, their expansion was linked with increases in both extinction and speciation rates in other coral families, driving a nine-fold increase in lineage turnover. These results reveal a double-edged effect of diversity dependency on reef evolution. Given their fast growth, staghorn corals may have increased extinction rates via competitive interactions, while promoting speciation through their role as ecosystem engineers. This suggests that recent widespread human-mediated reductions in staghorn coral cover, may be disrupting the key macroevolutionary processes that established modern coral reef ecosystems.

The effect of biotic interactions among reef corals on macroevolutionary patterns is unclear. Here, the authors study the rich coral fossil record, finding that reef coral diversity experienced potentially biotic interaction-driven evolutionary rate changes, and that Staghorn corals affected fossil diversity trajectories of other coral groups.

The multicausal twilight of South American native mammalian predators (Metatheria, Sparassodonta)

Sparassodonts were the apex mammalian predators of South America throughout most of the Cenozoic, diversifying into a wide array of niches including fox-like and even saber-toothed forms. Their extinction is still controversial, with different authors suggesting competition with other predators (placental carnivorans, terror birds, and carnivorous opossums), extinction of prey, and climate change as causal explanations. Here, we analyse these hypotheses using a novel approach implicating Bayesian analyses. We find that speciation and extinction rates of sparassodonts can be correlated with (i) intrinsic biotic factors such as changes in body mass and diversity of sparassodonts, (ii) extrinsic biotic factors such as potential prey diversity, and iii) extrinsic abiotic factors like the atmospheric CO₂, sea level, temperature, and uplift of the Andes. Thus, sparassodonts are a good example of a multilevel mixed model of evolution, where various factors drove the evolutionary history of this clade in a pluralistic way. There is no evidence for competition between Sparassodonta and others predators, and the effect of competition in the face of extinctions of fossil species should be tested and not assumed. Furthermore, we propose a novel approach for evaluating the fossil record when performing macroevolutionary analyses.

Distributions of extinction times from fossil ages and tree topologies: the example of mid-Permian synapsid extinctions

Given a phylogenetic tree that includes only extinct, or a mix of extinct and extant taxa, where at least some fossil data are available, we present a method to compute the distribution of the extinction time of a given set of taxa under the Fossilized-Birth-Death model. Our approach differs from the previous ones in that it takes into account (i) the possibility that the taxa or the clade considered may diversify before going extinct and (ii) the whole phylogenetic tree to estimate extinction times, whilst previous methods do not consider the diversification process and deal with each branch independently. Because of this, our method can estimate extinction times of lineages represented by a single fossil, provided that they belong to a clade that includes other fossil occurrences. We assess and compare our new approach with a standard previous one using simulated data. Results show that our method provides more accurate confidence intervals. This new approach is applied to the study of the extinction time of three Permo-Carboniferous synapsid taxa (Ophiacodontidae, Edaphosauridae, and Sphenacodontidae) that are thought to have disappeared toward the end of the Cisuralian (early Permian), or possibly shortly thereafter. The timing of extinctions of these three taxa and of their component lineages supports the idea that the biological crisis in the late Kungurian/early Roadian consisted of a progressive decline in biodiversity throughout the Kungurian.

The largest hoplophonine and a complex new hypothesis of nimravid evolution

Nimravids were the first carnivorans to evolve saberteeth, but previously portrayed as having a narrow evolutionary trajectory of increasing degrees of sabertooth specialization. Here I present a novel hypothesis about the evolution of this group, including a description of Eusmilus adelos, the largest known hoplophonine, which forces a re-evaluation of not only their relationships, but perceived paleoecology. Using a tip-dated Bayesian analysis with sophisticated evolutionary models, nimravids can now be viewed as following two paths of evolution: one led to numerous early dirk-tooth forms, including E. adelos, while the other converged on living feline morphology, tens of millions of years before its appearance in felids.

Dinosaur biodiversity declined well before the asteroid impact, influenced by ecological and environmental pressures

The question why non-avian dinosaurs went extinct 66 million years ago (Ma) remains unresolved because of the coarseness of the fossil record. A sudden extinction caused by an asteroid is the most accepted hypothesis but it is debated whether dinosaurs were in decline or not before the impact. We analyse the speciation-extinction dynamics for six key dinosaur families, and find a decline across dinosaurs, where diversification shifted to a declining-diversity pattern ~76 Ma. We investigate the influence of ecological and physical factors, and find that the decline of dinosaurs was likely driven by global climate cooling and herbivorous diversity drop. The latter is likely due to hadrosaurs outcompeting other herbivores. We also estimate that extinction risk is related to species age during the decline, suggesting a lack of evolutionary novelty or adaptation to changing environments. These results support an environmentally driven decline of non-avian dinosaurs well before the asteroid impact.

Unifying Phylogenetic Birth-Death Models in Epidemiology and Macroevolution

Birth–death stochastic processes are the foundations of many phylogenetic models and are widely used to make inferences about epidemiological and macroevolutionary dynamics. There are a large number of birth–death model variants that have been developed; these impose different assumptions about the temporal dynamics of the parameters and about the sampling process. As each of these variants was individually derived, it has been difficult to understand the relationships between them as well as their precise biological and mathematical assumptions. Without a common mathematical foundation, deriving new models is nontrivial. Here, we unify these models into a single framework, prove that many previously developed epidemiological and macroevolutionary models are all special cases of a more general model, and illustrate the connections between these variants. This unification includes both models where the process is the same for all lineages and those in which it varies across types. We also outline a straightforward procedure for deriving likelihood functions for arbitrarily complex birth–death(-sampling) models that will hopefully allow researchers to explore a wider array of scenarios than was previously possible. By rederiving existing single-type birth–death sampling models, we clarify and synthesize the range of explicit and implicit assumptions made by these models. [Birth–death processes; epidemiology; macroevolution; phylogenetics; statistical inference.]

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

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

Evolutionary drivers, morphological evolution and diversity dynamics of a surviving mammal clade: cainotherioids at the Eocene-Oligocene transition

The Eocene-Oligocene transition (EOT) represents a period of global environmental changes particularly marked in Europe and coincides with a dramatic biotic turnover. Here, using an exceptional fossil preservation, we document and analyse the diversity dynamics of a mammal clade, Cainotherioidea (Artiodactyla), that survived the EOT and radiated rapidly immediately after. We infer their diversification history from Quercy Konzentrat-Lagerstätte (south-west France) at the species level using Bayesian birth-death models. We show that cainotherioid diversity fluctuated through time, with extinction events at the EOT and in the late Oligocene, and a major speciation burst in the early Oligocene. The latter is in line with our finding that cainotherioids had a high morphological adaptability following environmental changes throughout the EOT, which probably played a key role in the survival and evolutionary success of this clade in the aftermath. Speciation is positively associated with temperature and continental fragmentation in a time-continuous way, while extinction seems to synchronize with environmental change in a punctuated way. Within-clade interactions negatively affected the cainotherioid diversification, while inter-clade competition might explain their final decline during the late Oligocene. Our results provide a detailed dynamic picture of the evolutionary history of a mammal clade in a context of global change.

Ecological and biogeographic drivers of biodiversity cannot be resolved using clade age-richness data

Estimates of evolutionary diversification rates - speciation and extinction - have been used extensively to explain global biodiversity patterns. Many studies have analyzed diversification rates derived from just two pieces of information: a clade's age and its extant species richness. This "age-richness rate" (ARR) estimator provides a convenient shortcut for comparative studies, but makes strong assumptions about the dynamics of species richness through time. Here we demonstrate that use of the ARR estimator in comparative studies is problematic on both theoretical and empirical grounds. We prove mathematically that ARR estimates are non-identifiable: there is no information in the data for a single clade that can distinguish a process with positive net diversification from one where net diversification is zero. Using paleontological time series, we demonstrate that the ARR estimator has no predictive ability for real datasets. These pathologies arise because the ARR inference procedure yields "point estimates" that have been computed under a saturated statistical model with zero degrees of freedom. Although ARR estimates remain useful in some contexts, they should be avoided for comparative studies of diversification and species richness.

Re-evaluation of the "law of constant extinction" for ruminants at different taxonomical scales

The "law of constant extinction," proposed by Van Valen, states that long and short-lived taxa have equal chances of going extinct. This pattern of age-independent extinction was originally inferred using the fossil record of several different taxa and relied on survivorship curves built from the literal reading of the fossil record. Van Valen's seminal work was mostly done at higher taxonomic levels, hence its prevalence at the species level could not be directly inferred. The surprisingly few subsequent studies done at the species level have challenged the prevalence of age-independent extinction, but those have, for the most part, failed to explicitly incorporate inherent biases of the fossil record. Using a recent Bayesian framework that accounts for several of those biases, including the fact that very short-living lineages might never make to the record itself, we showed that Ruminantia species present age-dependent extinction, where extinction probability decreases with species age. An analysis at the genus level suggested age-independent extinction but further examination suggested that the pattern might be more complex than previously reported by Van Valen. Our results indicate that different taxonomic levels may present different extinction regimes, which could justify the development of new macroevolutionary theory and methods.

The fossil record of extant elasmobranchs

Sharks and their relatives (Elasmobranchii) are highly threatened with extinction due to various anthropogenic pressures. The abundant fossil record of fossil taxa has allowed the tracing of the evolutionary history of modern elasmobranchs to at least 250 MYA; nonetheless, exactly how far back the fossil record of living taxa goes has never been collectively surveyed. In this study, the authors assess the representation and extent of the fossil record of elasmobranchs currently living in our oceans by collecting their oldest records and quantifying first appearance dates at different taxonomic levels (i.e., orders, families, genera and species), ecological traits (e.g., body size, habitat and feeding mechanism) and extinction risks (i.e., threatened, not threatened and data deficient). The results of this study confirm the robust representation of higher taxonomic ranks, with all orders, most of the families and over half of the extant genera having a fossil record. Further, they reveal that 10% of the current global species diversity is represented in the geological past. Sharks are better represented and extend deeper in time than rays and skates. While the fossil record of extant genera (e.g., the six gill sharks, Hexanchus) goes as far back as c. 190 MYA, the fossil record of extant species (e.g., the sand shark, Carcharias taurus Rafinesque 1810) extends c. 66 MYA. Although no significant differences were found in the extent of the fossil record between ecological traits, it was found that the currently threatened species have a significantly older fossil record than the not threatened species. This study demonstrate that the fossil record of extant elasmobranchs extends deep into the geologic time, especially in the case of threatened sharks. As such, the elasmobranch geological history has great potential to advance the understanding of how species currently facing extinction have responded to different stressors in the past, thereby providing a deep-time perspective to conservation.

52 million years old Eucalyptus flower sheds more than pollen grains

PREMISE

Fossils provide fundamental evidence of the evolutionary processes that crafted today's biodiversity and consequently for understanding life on Earth. We report the finding of Myrtaceidites eucalyptoides pollen grains preserved within the anthers of a 52-million-year-old Eucalyptus flower collected at Laguna del Hunco locality of Argentinean Patagonia and discuss its implications in understanding the evolutionary history of the iconic Australian genus Eucalyptus.

METHODS

Pollen grains were extracted from the flower's anthers and were then observed under light microscopy and scanning electron microscopy. The phylogenetic position of the fossil was investigated by adding pollen data to a previously published total-evidence matrix and analyzing it using parsimony.

RESULTS

We erect the species Eucalyptus xoshemium for the fossil flower. Pollen extracted from E. xoshemium belongs to the species Myrtaceidites eucalyptoides, which, until now, was only known as dispersed pollen. The numerous pollen grains recovered from the single flower allowed estimation of M. eucalyptoides' variability. Results of the phylogenetic analysis reinforce the position of this fossil within crown group Eucalyptus.

CONCLUSIONS

The discovery of these pollen grains within a Patagonian Eucalyptus fossil flower confirms the hypothesis that Myrtaceidites eucalyptoides represents fossil pollen in the Eucalyptus lineage, extends the geographic and stratigraphic fossil pollen record, and supports an earlier age for crown-group eucalypts.

Investigating Biotic Interactions in Deep Time

Recent renewed interest in using fossil data to understand how biotic interactions have shaped the evolution of life is challenging the widely held assumption that long-term climate changes are the primary drivers of biodiversity change. New approaches go beyond traditional richness and co-occurrence studies to explicitly model biotic interactions using data on fossil and modern biodiversity. Important developments in three primary areas of research include analysis of (i) macroevolutionary rates, (ii) the impacts of and recovery from extinction events, and (iii) how humans (Homo sapiens) affected interactions among non-human species. We present multiple lines of evidence for an important and measurable role of biotic interactions in shaping the evolution of communities and lineages on long timescales.

The past and future human impact on mammalian diversity

To understand the current biodiversity crisis, it is crucial to determine how humans have affected biodiversity in the past. However, the extent of human involvement in species extinctions from the Late Pleistocene onward remains contentious. Here, we apply Bayesian models to the fossil record to estimate how mammalian extinction rates have changed over the past 126,000 years, inferring specific times of rate increases. We specifically test the hypothesis of human-caused extinctions by using posterior predictive methods. We find that human population size is able to predict past extinctions with 96% accuracy. Predictors based on past climate, in contrast, perform no better than expected by chance, suggesting that climate had a negligible impact on global mammal extinctions. Based on current trends, we predict for the near future a rate escalation of unprecedented magnitude. Our results provide a comprehensive assessment of the human impact on past and predicted future extinctions of mammals.

Ancient tropical extinctions at high latitudes contributed to the latitudinal diversity gradient

Global biodiversity currently peaks at the equator and decreases toward the poles. Growing fossil evidence suggest this hump-shaped latitudinal diversity gradient (LDG) has not been persistent through time, with similar diversity across latitudes flattening out the LDG during past greenhouse periods. However, when and how diversity declined at high latitudes to generate the modern LDG remains an open question. Although diversity-loss scenarios have been proposed, they remain mostly undemonstrated. We outline the "asymmetric gradient of extinction and dispersal" framework that contextualizes previous ideas behind the LDG under a time-variable scenario. Using phylogenies and fossils of Testudines, Crocodilia, and Lepidosauria, we find that the hump-shaped LDG could be explained by (1) disproportionate extinctions of high-latitude tropical-adapted clades when climate transitioned from greenhouse to icehouse, and (2) equator-ward biotic dispersals tracking their climatic preferences when tropical biomes became restricted to the equator. Conversely, equivalent diversification rates across latitudes can account for the formation of an ancient flat LDG. The inclusion of fossils in macroevolutionary studies allows revealing time-dependent extinction rates hardly detectable from phylogenies only. This study underscores that the prevailing evolutionary processes generating the LDG during greenhouses differed from those operating during icehouses.

Fossil and phylogenetic analyses reveal recurrent periods of diversification and extinction in dictyopteran insects

Variations of speciation and extinction rates determine the fate of clades through time. Periods of high diversification and extinction (possibly mass-extinction events) can punctuate the evolutionary history of various clades, but they remain loosely defined for many biological groups, especially nonmarine invertebrates like insects. Here, we examine whether the cockroaches, mantises and termites (altogether included in Dictyoptera) have experienced episodic pulses of speciation or extinction and how these pulses may be associated with environmental fluctuations or mass extinctions. We relied on molecular phylogeny and fossil data to shed light on the times and rates at which dictyopterans diversified. The diversification of Dictyoptera has alternated between (i) periods of high diversification in the late Carboniferous, Early-Middle Triassic, Early Cretaceous and middle Palaeogene, and (ii) periods of high extinction rates particularly at the Permian-Triassic boundary, but not necessarily correlated with the major global biodiversity crises as in the mid-Cretaceous. This study advocates the importance of analyzing, when possible, both molecular phylogeny and fossil data to unveil diversification and extinction periods for a given group. The causes and consequences of extinction must be studied beyond mass-extinction events alone to gain a broader understanding of how clades wax and wane.

Selective extinction against redundant species buffers functional diversity

The extinction of species can destabilize ecological processes. A way to assess the ecological consequences of species loss is by examining changes in functional diversity. The preservation of functional diversity depends on the range of ecological roles performed by species, or functional richness, and the number of species per role, or functional redundancy. However, current knowledge is based on short timescales and an understanding of how functional diversity responds to long-term biodiversity dynamics has been limited by the availability of deep-time, trait-based data. Here, we compile an exceptional trait dataset of fossil molluscs from a 23-million-year interval in the Caribbean Sea (34 011 records, 4422 species) and develop a novel Bayesian model of multi-trait-dependent diversification to reconstruct mollusc (i) diversity dynamics, (ii) changes in functional diversity, and (iii) extinction selectivity over the last 23 Myr. Our results identify high diversification between 23–5 Mya, leading to increases in both functional richness and redundancy. Conversely, over the last three million years, a period of high extinction rates resulted in the loss of 49% of species but only 3% of functional richness. Extinction rates were significantly higher in small, functionally redundant species suggesting that competition mediated the response of species to environmental change. Taken together, our results identify long-term diversification and selective extinction against redundant species that allowed functional diversity to grow over time, ultimately buffering the ecological functions of biological communities against extinction.

A quantitative workflow for modeling diversification in material culture

Questions about the evolution of material culture are widespread in the humanities and social sciences. Statistical modeling of long-term changes in material culture is less common due to a lack of appropriate frameworks. Our goal is to close this gap and provide robust statistical methods for examining changes in the diversity of material culture. We provide an open-source and quantitative workflow for estimating rates of origination, extinction, and preservation, as well as identifying key shift points in the diversification histories of material culture. We demonstrate our approach using two distinct kinds of data: age ranges for the production of American car models, and radiocarbon dates associated with archaeological cultures of the European Neolithic. Our approach improves on existing frameworks by disentangling the relative contributions of origination and extinction to diversification. Our method also permits rigorous statistical testing of competing hypotheses to explain changes in diversity. Finally, we stress the value of a flexible approach that can be applied to data in various forms; this flexibility allows scholars to explore commonalities between forms of material culture and ask questions about the general properties of cultural change.

The dynamics of stem and crown groups

The fossil record of the origins of major groups such as animals and birds has generated considerable controversy, especially when it conflicts with timings based on molecular clock estimates. Here, we model the diversity of "stem" (basal) and "crown" (modern) members of groups using a "birth-death model," the results of which qualitatively match many large-scale patterns seen in the fossil record. Typically, the stem group diversifies rapidly until the crown group emerges, at which point its diversity collapses, followed shortly by its extinction. Mass extinctions can disturb this pattern and create long stem groups such as the dinosaurs. Crown groups are unlikely to emerge either cryptically or just before mass extinctions, in contradiction to popular hypotheses such as the "phylogenetic fuse". The patterns revealed provide an essential context for framing ecological and evolutionary explanations for how major groups originate, and strengthen our confidence in the reliability of the fossil record.

A Cretaceous peak in family-level insect diversity estimated with mark-recapture methodology

The history of insects' taxonomic diversity is poorly understood. The two most common methods for estimating taxonomic diversity in deep time yield conflicting results: the 'range through' method suggests a steady, nearly monotonic increase in family-level diversity, whereas 'shareholder quorum subsampling' suggests a highly volatile taxonomic history with family-level mass extinctions occurring repeatedly, even at the midpoints of geological periods. The only feature shared by these two diversity curves is a steep increase in standing diversity during the Early Cretaceous. This apparent diversification event occurs primarily during the Aptian, the pre-Cenozoic interval with the most described insect occurrences, raising the possibility that this feature of the diversity curves reflects preservation and sampling biases rather than insect evolution and extinction. Here, the capture-mark-recapture (CMR) approach is used to estimate insects' family-level diversity. This method accounts for the incompleteness of the insect fossil record as well as uneven sampling among time intervals. The CMR diversity curve shows extinctions at the Permian/Triassic and Cretaceous/Palaeogene boundaries but does not contain any mass extinctions within geological periods. This curve also includes a steep increase in diversity during the Aptian, which appears not to be an artefact of sampling or preservation bias because this increase still appears when time bins are standardized by the number of occurrences they contain rather than by the amount of time that they span. The Early Cretaceous increase in family-level diversity predates the rise of angiosperms by many millions of years and can be better attributed to the diversification of parasitic and especially parasitoid insect lineages.

Enforced symmetry: the necessity of symmetric waxing and waning

A fundamental question in ecology is how the success of a taxon changes through time and what drives this change. This question is commonly approached using trajectories averaged over a group of taxa. Using results from probability theory, we show analytically and using examples that averaged trajectories will be more symmetric as the number of averaged trajectories increases, even if none of the original trajectories they were derived from is symmetric. This effect is not only based on averaging, but also on the introduction of noise and the incorporation of a priori known origination and extinction times. This implies that averaged trajectories are not suitable for deriving information about the processes driving the success of taxa. In particular, symmetric waxing and waning, which is commonly observed and interpreted to be linked to a number of different paleobiological processes, does not allow drawing any conclusions about the nature of the underlying process.

Climate cooling and clade competition likely drove the decline of lamniform sharks

Understanding heterogeneity in species richness between closely related clades is a key research question in ecology and evolutionary biology. Multiple hypotheses have been proposed to interpret such diversity contrasts across the tree of life, with most studies focusing on speciation rates to explain clades' evolutionary radiations, while often neglecting extinction rates. Here we study a notorious biological model as exemplified by the sister relationships between mackerel sharks (Lamniformes, 15 extant species) and ground sharks (Carcharhiniformes, ∼290 extant species). Using a comprehensive fossil dataset, we found that the diversity dynamics of lamniforms waxed and waned following repeated cycles of radiation phases and declining phases. Radiation phases peaked up to 3 times the current diversity in the early Late Cretaceous. In the last 20 million years, the group declined to its present-day diversity. Along with a higher extinction risk for young species, we further show that this declining pattern is likely attributed to a combination of abiotic and biotic factors, with a cooling-driven extinction (negative correlation between temperature and extinction) and clade competition with some ground sharks. Competition from multiple clades successively drove the demise and replacement of mackerel sharks due to a failure to originate facing the rise of ground sharks, particularly since the Eocene. These effects came from ecologically similar carcharhiniform species inhibiting diversification of medium- and large-sized lamniforms. These results imply that the interplay between abiotic and biotic drivers had a substantial role in extinction and speciation, respectively, which determines the sequential rise and decline of marine apex predators.