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