Estimating Age-Dependent Extinction: Contrasting Evidence from Fossils and Phylogenies.

Rates of Evolution

Estimates of diversification rates are invaluable for many macroevolutionary studies. Recently, an approach called BAMM (Bayesian Analysis of Macro-evolutionary Mixtures) has become widely used for estimating diversification rates and rate shifts. At the same time, several articles have concluded that estimates of net diversification rates from the method-of-moments (MS) estimators are inaccurate. Yet, no studies have compared the ability of these two methods to accurately estimate clade diversification rates. Here, we use simulations to compare their performance. We found that BAMM yielded relatively weak relationships between true and estimated diversification rates. This occurred because BAMM underestimated the number of rates shifts across each tree, and assigned high rates to small clades with low rates. Errors in both speciation and extinction rates contributed to these errors, showing that using BAMM to estimate only speciation rates is also problematic. In contrast, the MS estimators (particularly using stem group ages), yielded stronger relationships between true and estimated diversification rates, by roughly twofold. Furthermore, the MS approach remained relatively accurate when diversification rates were heterogeneous within clades, despite the widespread assumption that it requires constant rates within clades. Overall, we caution that BAMM may be problematic for estimating diversification rates and rate shifts.

For mammals today, mountains are diverse ecosystems globally, yet the strong relationship between species richness and topographic complexity is not a persistent feature of the fossil record. Based on fossil-occurrence data, diversity and diversification rates in the intermontane western North America varied through time, increasing significantly during an interval of global warming and regional intensification of tectonic activity from 18 to 14 Ma. However, our ability to infer origination and extinction rates reliably from the fossil record is affected by variation in preservation history. To investigate the influence of preservation on estimates of diversification rates, I simulated fossil records under four alternative diversification hypotheses and six preservation scenarios. Diversification hypotheses included tectonically controlled speciation pulses, while preservation scenarios were based on common trends (e.g., increasing rock record toward the present) or derived from fossil occurrences and the continental rock record. For each scenario, I estimated origination, extinction, and diversification rates using three standard methods—per capita, three-timer, and capture–mark–recapture (CMR) metrics—and evaluated the ability of the simulated fossil records to accurately recover the underlying diversification dynamics. Despite variable and low preservation probabilities, simulated fossil records retained the signal of true rates in several of the scenarios. The three metrics did not exhibit similar behavior under each preservation scenario: while three-timer and CMR metrics produced more accurate rate estimates, per capita rates tended to better reproduce true shifts in origination rates. All metrics suffered from spurious peaks in origination and extinction rates when highly volatile preservation impacted the simulated record. Results from these simulations indicate that elevated diversification rates in relation to tectonic activity during the middle Miocene are likely to be evident in the fossil record, even if preservation in the North American fossil record was variable. Input from the past is necessary to evaluate the ultimate mechanisms underlying speciation and extinction dynamics.

There is no scientific consensus about whether and how species' evolutionary age, or the elapsed time since their origination, might affect their probability of going extinct. Different age-dependent extinction (ADE) patterns have been proposed in theoretical and empirical studies, while the existence of a consistent and universal pattern across the tree of life remains debated. If evolutionary age predicts species extinction probability, then the study of ADE should comprise the elapsed time and the ecological process acting on species from their origin to their extinction or to the present for extant species. Additionally, given that closely related species share traits associated with fitness, evolutionary proximity could generate similar ADE patterns. Considering the historical context and extinction selectivity based on evolutionary relatedness, we build on previous theoretical work to formalize the Clade Replacement Theory (CRT) as a framework that considers the ecological and evolutionary aspects of species age and extinction probability to produce testable predictions on ADE patterns. CRT's domain is the diversification dynamics of two or more clades competing for environmental space throughout time, and its propositions or derived hypotheses are as follows: (i) incumbency effects by an early arriving clade that limit the colonization and the diversification of a younger clade leading to a negative ADE scenario (younger species more prone to extinction than older ones) and (ii) an ecological shift triggered by an environmental change that imposes a new selective regime over the environmental space and leads to a positive ADE scenario (extinction probability increasing with age). From these propositions, we developed the prediction that the ADE scenario would be defined by whether an ecological shift happens or not. We discuss how the CRT could be tested with empirical data and provide examples where it could be applied. We hope this article will provide a common ground to unify results from different fields and foster new empirical tests of the mechanisms derived here while providing insights into CRT theoretical structuration.

As Van Valen has demonstrated, the taxonomic survivorship curve is a valuable means of investigating extinction rates in the fossil record. He suggested that within an adaptive zone, related taxa display stochastically constant and equal extinction rates. Such a condition is evidenced by straight survivorship curves for species and higher taxa. Van Valen's methods of survivorship analysis can be improved upon and several suggestions are presented. With proper manipulation of data, it is possible to pool the information from extinct and living taxa to produce a single survivorship curve and therefore a single estimate of extinction rate. If extinction rate is constant at the species level (producing a straight survivorship curve), higher taxa in the same group should be expected to have convex survivorship curves. The constancy of extinction rates (here termed Van Valen's Law) can and should be tested rigorously. Several methods are available, of which the Total Life method of Epstein is particularly effective.

Biases in the survivorship curves of fossil taxa

The estimation of origination and extinction rates and their temporal variation is central to understanding diversity patterns and the evolutionary history of clades. The fossil record provides the only direct evidence of extinction and biodiversity changes through time and has long been used to infer the dynamics of diversity changes in deep time. The software PyRate implements a Bayesian framework to analyze fossil occurrence data to estimate the rates of preservation, origination, and extinction while incorporating several sources of uncertainty. Building upon this framework, we present a suite of methodological advances including more complex and realistic models of preservation and the first likelihood-based test to compare the fit across different models. Further, we develop a new reversible jump Markov chain Monte Carlo algorithm to estimate origination and extinction rates and their temporal variation, which provides more reliable results and includes an explicit estimation of the number and temporal placement of statistically significant rate changes. Finally, we implement a new C++ library that speeds up the analyses by orders of magnitude, therefore facilitating the application of the PyRate methods to large data sets. We demonstrate the new functionalities through extensive simulations and with the analysis of a large data set of Cenozoic marine mammals. We compare our analytical framework against two widely used alternative methods to infer origination and extinction rates, revealing that PyRate decisively outperforms them across a range of simulated data sets. Our analyses indicate that explicit statistical model testing, which is often neglected in fossil-based macroevolutionary analyses, is crucial to obtain accurate and robust results.

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.

For some higher taxa, species can be identified in the fossil record with a high degree of reliability. The great geological durations of species indicate that phyletic evolution is normally so slow that little change occurs within a lineage during 105–107 generations. Failure to recognize sibling species in the fossil record has no bearing on this conclusion because they embody virtually no morphological change. Although slowness is the rule, we have no more precise assessment of morphological rates of phyletic evolution for any major taxon. Morphological data that have been assembled to assess rates of phyletic evolution have been meager, unrepresentative, and predominantly reflective of nothing more than body size. Net selection pressures within long segments of phylogeny—even ones that embrace large amounts of evolution—are infinitesimal and seemingly unsustainable against random fluctuations. This suggests that natural selection operates in a highly episodic fashion.Rates of adaptive radiation and extinction at the species level can be estimated for many taxa and, from them, rates of speciation in adaptive radiation. Species selection should universally tend to increase rate of speciation and decrease rate of extinction, yet these rates are positively correlated in the animal world, apparently because they are linked by common controls: both rate of speciation and rate of extinction seem to increase with level of stereotypical behavior and to decrease with dispersal ability. Only a few “supertaxa” have been able to combine high rates of speciation with moderate rates of extinction.

According to when they attained high diversity, major taxa of marine animals have been clustered into three groups, the Cambrian, Paleozoic, and Modern Faunas. Because the Cambrian Fauna was a relatively minor component of the total fauna after mid-Ordovician time, the Phanerozoic history of marine animal diversity is largely a matter of the fates of the Paleozoic and Modern Faunas. The fact that most late Cenozoic genera belong to taxa that have been radiating for tens of millions of years indicates that the post-Paleozoic increase in diversity indicated by fossil data is real, rather than an artifact of improvement of the fossil record toward the present.Assuming that ecological crowding produced the so-called Paleozoic plateau for family diversity, various workers have used the logistic equation of ecology to model marine animal diversification as damped exponential increase. Several lines of evidence indicate that this procedure is inappropriate. A plot of the diversity of marine animal genera through time provides better resolution than the plot for families and has a more jagged appearance. Generic diversity generally increased rapidly during the Paleozoic, except when set back by pulses of mass extinction. In fact, an analysis of the history of the Paleozoic Fauna during the Paleozoic Era reveals no general correlation between rate of increase for this fauna and total marine animal diversity. Furthermore, realistically scaled logistic simulations do not mimic the empirical pattern. In addition, it is difficult to imagine how some fixed limit for diversity could have persisted throughout the Paleozoic Era, when the ecological structure of the marine ecosystem was constantly changing. More fundamentally, the basic idea that competition can set a limit for marine animal diversity is incompatible with basic tenets of marine ecology: predation, disturbance, and vagaries of recruitment determine local population sizes for most marine species. Sparseness of predators probably played a larger role than weak competition in elevating rates of diversification during the initial (Ordovician) radiation of marine animals and during recoveries from mass extinctions. A plot of diversification against total diversity for these intervals yields a band of points above the one representing background intervals, and yet this band also displays no significant trend (if the two earliest intervals of the initial Ordovician are excluded as times of exceptional evolutionary innovation). Thus, a distinctive structure characterized the marine ecosystem during intervals of evolutionary radiation—one in which rates of diversification were exceptionally high and yet increases in diversity did not depress rates of diversification.Particular marine taxa exhibit background rates of origination and extinction that rank similarly when compared with those of other taxa. Rates are correlated in this way because certain heritable traits influence probability of speciation and probability of extinction in similar ways. Background rates of origination and extinction were depressed during the late Paleozoic ice age for all major marine invertebrate taxa, but remained correlated. Also, taxa with relatively high background rates of extinction experienced exceptionally heavy losses during biotic crises because background rates of extinction were intensified in a multiplicative manner; decimation of a large group of taxa of this kind in the two Permian mass extinctions established their collective identity as the Paleozoic Fauna.Characteristic rates of origination and extinction for major taxa persisted from Paleozoic into post-Paleozoic time. Because of the causal linkage between rates of origination and extinction, pulses of extinction tended to drag down overall rates of origination as well as overall rates of extinction by preferentially eliminating higher taxa having relatively high background rates of extinction. This extinction/origination ratchet depressed turnover rates for the residual Paleozoic Fauna during the Mesozoic Era. A decline of this fauna's extinction rate to approximately that of the Modern Fauna accounts for the nearly equal fractional losses experienced by the two faunas in the terminal Cretaceous mass extinction.Viewed arithmetically, the fossil record indicates slow diversification for the Modern Fauna during Paleozoic time, followed by much more rapid expansion during Mesozoic and Cenozoic time. When viewed more appropriately as depicting geometric—or exponential—increase, however, the empirical pattern exhibits no fundamental secular change: the background rate of increase for the Modern Fauna—the fauna that dominated post-Paleozoic marine diversity—simply persisted, reflecting the intrinsic origination and extinction rates of constituent taxa. Persistence of this overall background rate supports other evidence that the empirical record of diversification for marine animal life since Paleozoic time represents actual exponential increase. This enduring rate makes it unnecessary to invoke environmental change to explain the post-Paleozoic increase of marine diversity.Because of the resilience of intrinsic rates, an empirically based simulation that entails intervals of exponential increase for the Paleozoic and Modern Faunas, punctuated by mass extinctions, yields a pattern that is remarkably similar to the empirical pattern. It follows that marine animal genera and species will continue to diversify exponentially long into the future, barring disruption of the marine ecosystem by human-induced or natural environmental changes.

The estimation of diversification rates using phylogenetic data has attracted a lot of attention in the past decade. In this context, the analysis of incomplete phylogenies (e.g. phylogenies resolved at the family level but unresolved at the species level) has remained difficult. I present here a likelihood-based method to combine partly resolved phylogenies with taxonomic (species-richness) data to estimate speciation and extinction rates. This method is based on fitting a birth-and-death model to both phylogenetic and taxonomic data. Some examples of the method are presented with data on birds and on mammals. The method is compared with existing approaches that deal with incomplete phylogenies. Some applications and generalizations of the approach introduced in this paper are further discussed.

Are polyploids really evolutionary dead-ends (again)? A critical reappraisal of Mayrose etal. ().

The integration of Geographic Information System (GIS) methodology within a phylogenetic and statistical framework provides a background against which to evaluate the relationship between biogeographic changes and evolution in the fossil record. A case study based on patterns in Middle and Late Devonian phyllocarids (Crustacea) illustrates the usefulness of this integrated approach. Using a combined approach enhances determination of rates of biodiversity change and the relationship between biogeographic and evolutionary changes. Because the interaction between speciation and extinction rates fundamentally determines biodiversity dynamics, and speciation and extinction rates are influenced by the geographic ranges of component taxa, the relationship between biogeography and evolution is important. Furthermore, GIS makes it possible to quantify paleobiogeographic ranges.Phylogenetic biogeography resolved patterns of both vicariance and geodispersal and revealed that range expansions were more abundant than range contractions in Devonian phyllocarids. In addition, statistical tests on GIS-constrained species ranges and evolutionary-rate data revealed a relationship between increasing species' ranges and increases in both speciation and extinction rates. Extinction rate, however, increased more rapidly than speciation rate in the phyllocarids. The pattern of extinction rate increasing faster than speciation rate in the phyllocarids may illuminate aspects of the Late Devonian biodiversity crisis in particular, and protracted biodiversity crises in general.

Studies in insular environments have often documented a positive association of extinction risk and evolutionary uniqueness (i.e., how distant a species is from its closest living relative). However, the cause of this association is unclear. One explanation is that species threatened with extinction are evolutionarily unique because they are old, implying that extinction risk increases with time since speciation (age-dependent extinction). An alternative explanation is that such threatened species are last survivors of clades that have undergone an elevated extinction rate, and that their uniqueness results from the extinction of their close relatives. Distinguishing between these explanations is difficult but important, since they imply different biological processes determining extinction patterns. Here, we designed a simulation approach to distinguish between these alternatives using living species, and applied it to 12 insular radiations that show a positive association between extinction risk and evolutionary uniqueness. We also tested the sensitivity of results to underlying assumptions and variable extinction rates. Despite differences among the radiations considered, age-dependent extinction was supported as best explaining the majority of the empirical cases. Biological processes driving characteristic changes in abundance with species duration (age-dependency) may merit further investigation.

Author SummaryWhy are there more species in the tropics? This question has fascinated ecologists and evolutionary biologists for decades, generating hundreds of hypotheses, yet basic questions remain: Are rates of speciation higher in the tropics? Are rates of extinction higher in temperate regions? Do the tropics act as a source of diversity for temperate regions? We estimated rates of speciation, extinction, and range expansion associated with mammals living in tropical and temperate regions, using an almost complete mammalian phylogeny. Contrary to what has been suggested before for this class of vertebrates, we found that diversification rates are strikingly consistent with diversity patterns, with latitudinal peaks in species richness being associated with high speciation rates, low extinction rates, or both, depending on the mammalian order (rodents, bats, primates, etc.). We also found evidence for an asymmetry in range expansion, with more expansion “out of” than “into” the tropics. Taken together, these results suggest that tropical regions are not only a reservoir of biodiversity, but also the main place where biodiversity is generated.