Evolutionary Radiations and Extinctions

Biodiversity indices often combine data from different species when used in monitoring programs. Heuristic properties can suggest preferred indices, but we lack objective ways to discriminate between indices with similar heuristics. Biodiversity indices can be evaluated by determining how well they reflect management objectives that a monitoring program aims to support. For example, the Convention on Biological Diversity requires reporting about extinction rates, so simple indices that reflect extinction risk would be valuable. We developed 3 biodiversity indices that are based on simple models of population viability that relate extinction risk to abundance. We based the first index on the geometric mean abundance of species and the second on a more general power mean. In a third index, we integrated the geometric mean abundance and trend. These indices require the same data as previous indices, but they also relate directly to extinction risk. Field data for butterflies and woodland plants and experimental studies of protozoan communities show that the indices correlate with local extinction rates. Applying the index based on the geometric mean to global data on changes in avian abundance suggested that the average extinction probability of birds has increased approximately 1% from 1970 to 2009.Conectando Índices para el Monitoreo de la Biodiversidad con la Teoría de Riesgo de ExtinciónResumenLos índices de biodiversidad combinan frecuentemente los datos de diferentes especies cuando se usan en los programas de monitoreo. Las propiedades heurísticas pueden sugerir índices preferidos, pero carecemos de medios objetivos para discriminar a los índices con propiedades heurísticas similares. Los índices de biodiversidad pueden evaluarse al determinar qué tan bien reflejan los objetivos de manejo que un programa de monitoreo busca apoyar. Por ejemplo, la Convención sobre la Diversidad Biológica requiere reportar las tasas de extinción, así que los índices que reflejan el riesgo de extinción serían valiosos. Desarrollamos 3 índices de biodiversidad que se basan en modelos sencillos de viabilidad de población y que relacionan el riesgo de extinción con la abundancia. Basamos el primer índice en la media geométrica de la abundancia de especies, y el segundo en una media de poder más general. En el tercer índice integramos la media geométrica y la tendencia. Estos índices requieren los mismos datos que índices previos, pero también se relacionan directamente con el riesgo de extinción. La información de campo sobre mariposas y plantas de bosque, y los estudios experimentales de comunidades protozoarias, muestran que los índices se correlacionan con las tasas locales de extinción. Al aplicar el índice basado en la media geométrica sobre los datos globales de los cambios en la abundancia de aves, sugirió que la probabilidad de extinción promedio de aves ha incrementado aproximadamente 1% desde 1970 hasta 2009.Palabras ClaveÍndice de biodiversidad, media geométrica, medida de la biodiversidad, riesgo de extinción

Rates of Evolution

Life-history traits interact in important ways. Relatively few studies, however, have explored the relationships between life-history traits in long-lived taxa such as trees. We examined patterns of energy allocation to components of reproduction and growth in three species of California oaks (Quercus spp.) using a combination of annual acorn censuses, dendrometer bands to measure radial increment, and litterfall traps. Our results are generally consistent with the hypothesis that energy invested in reproduction detracts from the amount of energy available for growth in these long-lived taxa; i.e., there are trade-offs between these traits. The relationships between reproduction and growth varied substantially among specific trait combinations and tree species, however, and in some cases were in the direction opposite that expected based on the assumption of trade-offs between them. This latter finding appears to be a consequence of the pattern of resource use across years in these long-lived trees contrasting with the expected partitioning of resource use within years in short-lived taxa. Thus, the existence and magnitude of putative trade-offs varied depending on whether the time scale considered was within or across years. Collectively, our results indicate that negative relationships between fundamental life-history traits can be important at multiple levels of modular organization and that energy invested in reproduction can have measurable consequences in terms of the amount of energy available for future reproduction and both current and future growth.

Time-calibrated molecular phylogenies provide a valuable window into the tempo and mode of species diversification, especially for the large number of groups that lack adequate fossil records. Molecular phylogenetic data frequently suggest an initial "explosive speciation" phase, leading to widespread speculation that ecological niche-filling processes might govern the dynamics of species diversification during evolutionary radiations. However, these patterns are difficult to reconcile with the fossil record. The fossil record strongly suggests that extinction rates have been high relative to speciation rates, but such elevated background extinction should erase the signal of early, rapid speciation from molecular phylogenies. For this reason, extinction rates in molecular phylogenies are frequently estimated as zero under the widely used birth-death model. Here, I construct a simple model that combines phylogenetically patterned extinction with pulsed turnover dynamics and constant diversity through time. Using approximate Bayesian methods, I show that heritable extinction can easily explain the phenomenon of explosive early diversification, even when net diversification rates are equal to zero. Several assumptions of the model are more consistent with both the fossil record and neontological data than the standard birth-death model and it may thus represent a viable alternative interpretation of phylogenetic diversification patterns. These results suggest that variation in the absolute rate of lineage turnover through time, in conjunction with phylogenetically nonrandom extinction, may underlie the apparent diversity-dependent speciation observed in molecular phylogenies.

Recent application of time-varying birth-death models to molecular phylogenies suggests that a decreasing diversification rate can only be observed if there was a decreasing speciation rate coupled with extremely low or no extinction. However, from a paleontological perspective, zero extinction rates during evolutionary radiations seem unlikely. Here, with a more comprehensive set of computer simulations, we show that substantial extinction can occur without erasing the signal of decreasing diversification rate in a molecular phylogeny. We also find, in agreement with the previous work, that a decrease in diversification rate cannot be observed in a molecular phylogeny with an increasing extinction rate alone. Further, we find that the ability to observe decreasing diversification rates in molecular phylogenies is controlled (in part) by the ratio of the initial speciation rate (Lambda) to the extinction rate (Mu) at equilibrium (the LiMe ratio), and not by their absolute values. Here we show in principle, how estimates of initial speciation rates may be calculated using both the fossil record and the shape of lineage through time plots derived from molecular phylogenies. This is important because the fossil record provides more reliable estimates of equilibrium extinction rates than initial speciation rates.

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.

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.

Recent analyses have suggested that extinction and origination rates exhibit long-range correlations, implying that the fossil record may be controlled by self-organized criticality or other scale-free internal dynamics of the biosphere. Here we directly test for correlations in the fossil record by calculating the autocorrelation of extinction [corrected] and origination rates through time. Our results show that extinction rates are uncorrelated beyond the average duration of a stratigraphic interval. Thus, they lack the long-range correlations predicted by the self-organized criticality hypothesis. In contrast, origination rates show strong autocorrelations due to long-term trends. After detrending, origination rates generally show weak positive correlations at lags of 5-10 million years (Myr) and weak negative correlations at lags of 10-30 Myr, consistent with aperiodic oscillations around their long-term trends. We hypothesize that origination rates are more correlated than extinction rates because originations of new taxa create new ecological niches and new evolutionary pathways for reaching them, thus creating conditions that favour further diversification.

Short-term variations in rates of taxonomic extinction and origination in the fossil record may be the result of true changes in rates of turnover, variable rates of fossil preservation, or some combination of the two. Here, positive extinction and origination rate excursions among Phanerozoic marine animal genera are reexpressed in terms of the amount of normal, background time they represent. In addition to providing a background-adjusted calibration of rate intensities, this reexpression determines the durations of sampling gaps that would be required to explain entirely all observed rate excursions as sampling artifacts. This possibility is explored by analyzing a new compilation of the timing and duration of sedimentary hiatuses in North America. Hiatuses spanning more than approximately one million years (Myr) in the marine sedimentary rock record have a mean duration of 73 Myr. There are two major hiatus types—those that form in response to long-duration tectonic cycles and that bound the major Sloss-scale sequences (mean duration 200 Myr), and those that form in response to shorter-duration changes in sediment accommodation space and that occur within major Sloss-scale sequences (mean duration less than 23 Myr). The latter are approximately exponentially distributed and have a mean duration that is comparable to the mean duration of intervening sedimentary rock packages. Although sedimentary hiatuses are generally long enough in duration to accommodate the hypothesis that short-term variations in rates of genus origination and extinction are artifacts of sampling failures at major unconformities (“Unconformity Bias” hypothesis), the observed evolutionary rates are not correlated with hiatus durations. Moreover, hiatuses that follow the major mass extinctions are not long in comparison to most other non-mass extinction intervals. These results do not support the hypothesis that hiatuses at major unconformities alone have artificially clustered genus first and last occurrences, thereby causing many of the documented statistical similarities between the temporal structure of the sedimentary rock record and macroevolutionary patterns. Instead, environmental changes related to the expansion and contraction of marine environments may have been the primary forcers of both biological turnover and the spatio-temporal pattern of sediment accumulation. Fully testing this “Common Cause” hypothesis requires quantifying and overcoming lingering taxonomic, biostratigraphic, facies, and numerous other biases that are both inherent in geologic data and imposed by imperfect knowledge of the fossil record.

Comments on: Periodicity in the extinction rate and possible astronomical causes – comment on mass extinctions over the last 500 myr: an astronomical cause? (Erlykin <i>et al</i>.)

Extinction in the Anthropocene.

Evolutionary radiations have been extensively studied especially in the fossil record and in the context of postcrisis recoveries. The concept of adaptive radiation that emerges from this very broad topic explicitly involves the effect of adaptation driven by ecological opportunity and is considered to be of the foremost importance. It is essential to be able to detect adaptive radiation because it points up factors that predispose a clade to radiate. Adaptive radiation has received much attention in recent decades based mostly on studies dealing with recent clades, but data from the fossil record are still scarce. This study begins to fill this gap with the example ofLowerJurassic ammonoids (throughc. 8 Myr of history). A survey of several clades, using both taxonomic and disparity‐based approaches, shows that they diversified successively through time, but not systematically, in terms of species numbers and morphological variety. Some clades seem to have exhibited adaptive radiation and to have become rapidly extinct. One clade (which engendered nearly all post–LowerJurassic ammonoids) has a fossil record that begins with low diversity and disparity but is superseded by a sustained radiation pattern. The results are discussed in the light of theModernSynthesis and its continuation into anExtendedEvolutionarySynthesis.

The incompleteness of the fossil record hinders the inference of evolutionary rates and patterns. Here, we derive relationships among true taxonomic durations, preservation probability, and observed taxonomic ranges. We use these relationships to estimate original distributions of taxonomic durations, preservation probability, and completeness (proportion of taxa preserved), given only the observed ranges. No data on occurrences within the ranges of taxa are required. When preservation is random and the original distribution of durations is exponential, the inference of durations, preservability, and completeness is exact. However, reasonable approximations are possible given non-exponential duration distributions and temporal and taxonomic variation in preservability. Thus, the approaches we describe have great potential in studies of taphonomy, evolutionary rates and patterns, and genealogy. Analyses of Upper Cambrian-Lower Ordovician trilobite species, Paleozoic crinoid genera, Jurassic bivalve species, and Cenozoic mammal species yield the following results: (1) The preservation probability inferred from stratigraphic ranges alone agrees with that inferred from the analysis of stratigraphic gaps when data on the latter are available. (2) Whereas median durations based on simple tabulations of observed ranges are biased by stratigraphic resolution, our estimates of median duration, extinction rate, and completeness are not biased.(3) The shorter geologic ranges of mammalian species relative to those of bivalves cannot be attributed to a difference in preservation potential. However, we cannot rule out the contribution of taxonomic practice to this difference. (4) In the groups studied, completeness (proportion of species [trilobites, bivalves, mammals] or genera [crinoids] preserved) ranges from 60% to 90%. The higher estimates of completeness at smaller geographic scales support previous suggestions that the incompleteness of the fossil record reflects loss of fossiliferous rock more than failure of species to enter the fossil record in the first place.

Because both descriptions of species and modern human-driven extinctions started around the same time (i.e., eighteenth century), a logical expectation is that a large proportion of species may have gone extinct without ever having been recorded. Despite this evident and widely recognized assumption, the loss of undescribed species has never been estimated. We quantified this loss for several taxonomic groups and regions for which undescribed species extinctions are likely to have occurred. Across a wide range of taxonomic groups, we applied known extinction rates computed from recorded species losses to assumed exponential decay in the proportion of species remaining undiscovered. Because all previous modeling attempts to project total species richness implicitly assumed that undescribed species extinctions could be neglected, we also evaluated the effect of neglecting them. Finally, because we assumed constant description and extinction probabilities, we applied our model to simulated data that did not conform to this assumption. Actual species losses were severely underestimated by considering only known species extinctions. According to our estimates, the proportion of undiscovered extinct species over all extinctions ranged from 0.15 to 0.59, depending on the taxonomic group and the region considered. This means that recent extinctions may be up to twice as large as the number recorded. When species differed in their extinction or description probabilities, our model underestimated extinctions of undescribed species by up to 20%.

Recent Fourier analyses of fossil extinction data have indicated that the power spectrum of extinction during the Phanerozoic may take the form of 1/f noise, a result which, it has been suggested, could be indicative of the presence of `critical dynamics' in the processes giving rise to extinction. In this paper we examine extinction power spectra in some detail, using family-level data from two widely available compilations. We find that although the average form of the power spectrum roughly obeys the 1/f law, the spectrum can be represented more accurately by dividing it into two regimes: a low-frequency one which is well fit by an exponential, and a high-frequency one in which it follows a power law with a 1/f2 form. We give explanations for the occurrence of each of these behaviours and for the position of the crossover between them.