Detecting shifts in diversity limits from molecular phylogenies: what can we know?

Large complete species-level molecular phylogenies can provide the most direct information about the macroevolutionary history of clades having poor fossil records. However, extinction will ultimately erode evidence of pulses of rapid speciation in the deep past. Assessment of how well, and for how long, phylogenies retain the signature of such pulses has hitherto been based on a--probably untenable--model of ongoing diversity-independent diversification. Here, we develop two new tests for changes in diversification 'rules' and evaluate their power to detect sudden increases in equilibrium diversity in clades simulated with diversity-dependent speciation and extinction rates. Pulses of diversification are only detected easily if they occurred recently and if the rate of species turnover at equilibrium is low; rates reported for fossil mammals suggest that the power to detect a doubling of species diversity falls to 50 per cent after less than 50 Myr even with a perfect phylogeny of extant species. Extinction does eventually draw a veil over past dynamics, suggesting that some questions are beyond the limits of inference, but sudden clade-wide pulses of speciation can be detected after many millions of years, even when overall diversity is constrained. Applying our methods to existing phylogenies of mammals and angiosperms identifies intervals of elevated diversification in each.     

Primate diversification inferred from phylogenies and fossils

Biodiversity arises from the balance between speciation and extinction. Fossils record the origins and disappearance of organisms, and the branching patterns of molecular phylogenies allow estimation of speciation and extinction rates, but the patterns of diversification are frequently incongruent between these two data sources. I tested two hypotheses about the diversification of primates based on ～600 fossil species and 90% complete phylogenies of living species: (1) diversification rates increased through time; (2) a significant extinction event occurred in the Oligocene. Consistent with the first hypothesis, analyses of phylogenies supported increasing speciation rates and negligible extinction rates. In contrast, fossils showed that while speciation rates increased, speciation and extinction rates tended to be nearly equal, resulting in zero net diversification. Partially supporting the second hypothesis, the fossil data recorded a clear pattern of diversity decline in the Oligocene, although diversification rates were near zero. The phylogeny supported increased extinction ～34 Ma, but also elevated extinction ～10 Ma, coinciding with diversity declines in some fossil clades. The results demonstrated that estimates of speciation and extinction ignoring fossils are insufficient to infer diversification and information on extinct lineages should be incorporated into phylogenetic analyses.     

EXTINCTION DURING EVOLUTIONARY RADIATIONS: RECONCILING THE FOSSIL RECORD WITH 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.     

EXTINCTION RATES SHOULD NOT BE ESTIMATED FROM MOLECULAR PHYLOGENIES

Molecular phylogenies contain information about the tempo and mode of species diversification through time. Because extinction leaves a characteristic signature in the shape of molecular phylogenetic trees, many studies have used data from extant taxa only to infer extinction rates. This is a promising approach for the large number of taxa for which extinction rates cannot be estimated from the fossil record. Here, I explore the consequences of violating a common assumption made by studies of extinction from phylogenetic data. I show that when diversification rates vary among lineages, simple estimators based on the birth–death process are unable to recover true extinction rates. This is problematic for phylogenetic trees with complete taxon sampling as well as for the simpler case of clades with known age and species richness. Given the ubiquity of variation in diversification rates among lineages and clades, these results suggest that extinction rates should not be estimated in the absence of fossil data.     

A model of onshore-offshore change in faunal diversity

Onshore-offshore patterns of faunal change occurred at many taxonomic scales during the Paleozoic Era, ranging from replacement of the Cambrian evolutionary fauna by the Paleozoic fauna to the environmental expansion of many orders and classes. A simple mathematical model is constructed to investigate such change. The environmental gradient across the marine shelf-slope is treated as a linear array of discrete habitats, each of which holds a set number of species, as observed in the fossil record. During any interval of time, some portion of the species in each habitat becomes extinct by background processes, with rates of extinction varying among both clades and habitats, as also observed in the record. After extinction, species are replaced from within the habitat and from immediately adjacent habitats, with proportions dependent on surviving species. This model leads to the prediction that extinction-resistant clades will always diversify at the expense of extinction-prone clades. But if extinction intensity is highest in nearshore habitats, extinction-resistant clades will expand preferentially in the onshore direction, build up diversity there, and then diversify outward toward the offshore. Thus, onshore-offshore patterns of diversification may be the expectation for faunal change quite independently of whether or not clades originate onshore. When the model is parameterized for Paleozoic trilobites and brachiopods, numerical solutions exhibit both a pattern of faunal change and a time span for diversification similar to that seen in the fossil record. They also generate structure similar to that seen in global diversification, including logistic patterns of growth, declining origination but constant extinction within clades through time, and declining overall extinction across clades through time.     

EXPLOSIVE EVOLUTIONARY RADIATIONS: DECREASING SPECIATION OR INCREASING EXTINCTION THROUGH TIME?

A common pattern in time-calibrated molecular phylogenies is a signal of rapid diversification early in the history of a radiation. Because the net rate of diversification is the difference between speciation and extinction rates, such "explosive-early" diversification could result either from temporally declining speciation rates or from increasing extinction rates through time. Distinguishing between these alternatives is challenging but important, because these processes likely result from different ecological drivers of diversification. Here we develop a method for estimating speciation and extinction rates that vary continuously through time. By applying this approach to real phylogenies with explosive-early diversification and by modeling features of lineage-accumulation curves under both declining speciation and increasing extinction scenarios, we show that a signal of explosive-early diversification in phylogenies of extant taxa cannot result from increasing extinction and can only be explained by temporally declining speciation rates. Moreover, whenever extinction rates are high, "explosive early" patterns become unobservable, because high extinction quickly erases the signature of even large declines in speciation rates. Although extinction may obscure patterns of evolutionary diversification, these results show that decreasing speciation is often distinguishable from increasing extinction in the numerous molecular phylogenies of radiations that retain a preponderance of early lineages.     

Evolutionary diversification of reef corals: a comparison of the molecular and fossil records

Understanding historical patterns of diversity dynamics is of paramount importance for many evolutionary questions. The fossil record has long been the only source of information on patterns of diversification, but the molecular record, derived from time-calibrated phylogenies, is becoming an important additional resource. Both fossil and molecular approaches have shortcomings and biases. These have been well studied for fossil data but much less so for molecular data and empirical comparisons between approaches are lacking. Here, we compare the patterns of diversification derived from fossil and molecular data in scleractinian reef coral species. We also assess the robustness of molecular diversification rates to poor taxon sampling. We find that the temporal pattern of molecular diversification rates is robust to incomplete sampling when rates are calculated per interval. The major obstacle of molecular methods is that rate estimates are distorted because diversification rates can never be negative, whereas the fossil record suffers from incomplete preservation and inconsistent taxonomy. Nevertheless, the molecular pattern of diversification is comparable to the pattern we observe in the fossil record, with the timing of major diversification pulses coinciding in each dataset. For example, both agree that the end-Triassic coral extinction was a catastrophic bottleneck in scleractinian evolution.     

Incorporating phylogenetic uncertainty on phylogeny‐based palaeontological dating and the timing of turtle diversification

Methods improving the performance of molecular dating of divergence time of clades have improved dramatically in recent years. The calibration of molecular dating using the first appearance of a clade in the fossil record is a crucial step towards inferring the minimal diversification time of various groups and the choice of extinct taxa can strongly influence the molecular dates. Here, we evaluate the uncertainty on the phylogenetic position of extinct taxa through non‐parametric bootstrapping. The recognition of phylogenetic uncertainty resulted in the definition of the Bootstrap Uncertainty Range (BUR) for the age of first appearance of a given clade. The BUR is calculated as the interval of geological time in which the diversification of a given clade can be inferred to have occurred, based on the temporal information of the fossil record and the topologies of the bootstrap trees. Divergence times based on BUR analyses were calculated for three clades of turtles: Testudines, Pleurodira and Cryptodira. This resulted in extensive uncertainty ranges of topology‐dependent minimal divergence dates for these clades.     

Density-dependent diversification in North American wood warblers

Evidence from both molecular phylogenies and the fossil record suggests that rates of species diversification often decline through time during evolutionary radiations. One proposed explanation for this pattern is ecological opportunity, whereby an initial abundance of resources and lack of potential competitors facilitate rapid diversification. This model predicts density-dependent declines in diversification rates, but has not been formally tested in any species-level radiation. Here we develop a new conceptual framework that distinguishes density dependence from alternative processes that also produce temporally declining diversification, and we demonstrate this approach using a new phylogeny of North American Dendroica wood warblers. We show that explosive lineage accumulation early in the history of this avian radiation is best explained by a density-dependent diversification process. Our results suggest that the tempo of wood warbler diversification was mediated by ecological interactions among species and that lineage and ecological diversification in this group are coupled, as predicted under the ecological opportunity model.     

Density-dependent cladogenesis in birds

A characteristic signature of adaptive radiation is a slowing of the rate of speciation toward the present. On the basis of molecular phylogenies, studies of single clades have frequently found evidence for a slowdown in diversification rate and have interpreted this as evidence for density dependent speciation. However, we demonstrated via simulation that large clades are expected to show stronger slowdowns than small clades, even if the probability of speciation and extinction remains constant through time. This is a consequence of exponential growth: clades, which, by chance, diversify at above the average rate early in their history, will tend to be large. They will also tend to regress back to the average diversification rate later on, and therefore show a slowdown. We conducted a meta-analysis of the distribution of speciation events through time, focusing on sequence-based phylogenies for 45 clades of birds. Thirteen of the 23 clades (57%) that include more than 20 species show significant slowdowns. The high frequency of slowdowns observed in large clades is even more extreme than expected under a purely stochastic constant-rate model, but is consistent with the adaptive radiation model. Taken together, our data strongly support a model of density-dependent speciation in birds, whereby speciation slows as ecological opportunities and geographical space place limits on clade growth.     

Phylogenetic estimates of diversification rate are affected by molecular rate variation

Molecular phylogenies are increasingly being used to investigate the patterns and mechanisms of macroevolution. In particular, node heights in a phylogeny can be used to detect changes in rates of diversification over time. Such analyses rest on the assumption that node heights in a phylogeny represent the timing of diversification events, which in turn rests on the assumption that evolutionary time can be accurately predicted from DNA sequence divergence. But there are many influences on the rate of molecular evolution, which might also influence node heights in molecular phylogenies, and thus affect estimates of diversification rate. In particular, a growing number of studies have revealed an association between the net diversification rate estimated from phylogenies and the rate of molecular evolution. Such an association might, by influencing the relative position of node heights, systematically bias estimates of diversification time. We simulated the evolution of DNA sequences under several scenarios where rates of diversification and molecular evolution vary through time, including models where diversification and molecular evolutionary rates are linked. We show that commonly used methods, including metric‐based, likelihood and Bayesian approaches, can have a low power to identify changes in diversification rate when molecular substitution rates vary. Furthermore, the association between the rates of speciation and molecular evolution rate can cause the signature of a slowdown or speedup in speciation rates to be lost or misidentified. These results suggest that the multiple sources of variation in molecular evolutionary rates need to be considered when inferring macroevolutionary processes from phylogenies.     

BOOM AND BUST: ANCIENT AND RECENT DIVERSIFICATION IN BICHIRS (POLYPTERIDAE: ACTINOPTERYGII), A RELICTUAL LINEAGE OF RAY‐FINNED FISHES

Understanding the history that underlies patterns of species richness across the Tree of Life requires an investigation of the mechanisms that not only generate young species‐rich clades, but also those that maintain species‐poor lineages over long stretches of evolutionary time. However, diversification dynamics that underlie ancient species‐poor lineages are often hidden due to a lack of fossil evidence. Using information from the fossil record and time calibrated molecular phylogenies, we investigate the history of lineage diversification in Polypteridae, which is the sister lineage of all other ray‐finned fishes (Actinopterygii). Despite originating at least 390 million years (Myr) ago, molecular timetrees support a Neogene origin for the living polypterid species. Our analyses demonstrate polypterids are exceptionally species depauperate with a stem lineage duration that exceeds 380 million years (Ma) and is significantly longer than the stem lineage durations observed in other ray‐finned fish lineages. Analyses of the fossil record show an early Late Cretaceous (100.5–83.6 Ma) peak in polypterid genus richness, followed by 60 Ma of low richness. The Neogene species radiation and evidence for high‐diversity intervals in the geological past suggest a “boom and bust” pattern of diversification that contrasts with common perceptions of relative evolutionary stasis in so‐called “living fossils.”     

Detecting local diversity‐dependence in diversification

Whether there are ecological limits to species diversification is a hotly debated topic. Molecular phylogenies show slowdowns in lineage accumulation, suggesting that speciation rates decline with increasing diversity. A maximum‐likelihood (ML) method to detect diversity‐dependent (DD) diversification from phylogenetic branching times exists, but it assumes that diversity‐dependence is a global phenomenon and therefore ignores that the underlying species interactions are mostly local, and not all species in the phylogeny co‐occur locally. Here, we explore whether this ML method based on the nonspatial diversity‐dependence model can detect local diversity‐dependence, by applying it to phylogenies, simulated with a spatial stochastic model of local DD speciation, extinction, and dispersal between two local communities. We find that type I errors (falsely detecting diversity‐dependence) are low, and the power to detect diversity‐dependence is high when dispersal rates are not too low. Interestingly, when dispersal is high the power to detect diversity‐dependence is even higher than in the nonspatial model. Moreover, estimates of intrinsic speciation rate, extinction rate, and ecological limit strongly depend on dispersal rate. We conclude that the nonspatial DD approach can be used to detect diversity‐dependence in clades of species that live in not too disconnected areas, but parameter estimates must be interpreted cautiously.     

Macroevolution: dynamics of diversity

The fossil record typically exhibits very dynamic patterns of innovation, diversification and extinction. In contrast, molecular phylogenies suggest smoother patterns of evolutionary change. Several new studies reconcile this difference and reveal more about the mechanisms behind macroevolutionary change.     

Macroevolutionary Dynamics and Historical Biogeography of Primate Diversification Inferred from a Species Supermatrix

Phylogenetic relationships, divergence times, and patterns of biogeographic descent among primate species are both complex and contentious. Here, we generate a robust molecular phylogeny for 70 primate genera and 367 primate species based on a concatenation of 69 nuclear gene segments and ten mitochondrial gene sequences, most of which were extracted from GenBank. Relaxed clock analyses of divergence times with 14 fossil-calibrated nodes suggest that living Primates last shared a common ancestor 71–63 Ma, and that divergences within both Strepsirrhini and Haplorhini are entirely post-Cretaceous. These results are consistent with the hypothesis that the Cretaceous-Paleogene mass extinction of non-avian dinosaurs played an important role in the diversification of placental mammals. Previous queries into primate historical biogeography have suggested Africa, Asia, Europe, or North America as the ancestral area of crown primates, but were based on methods that were coopted from phylogeny reconstruction. By contrast, we analyzed our molecular phylogeny with two methods that were developed explicitly for ancestral area reconstruction, and find support for the hypothesis that the most recent common ancestor of living Primates resided in Asia. Analyses of primate macroevolutionary dynamics provide support for a diversification rate increase in the late Miocene, possibly in response to elevated global mean temperatures, and are consistent with the fossil record. By contrast, diversification analyses failed to detect evidence for rate-shift changes near the Eocene-Oligocene boundary even though the fossil record provides clear evidence for a major turnover event (“Grande Coupure”) at this time. Our results highlight the power and limitations of inferring diversification dynamics from molecular phylogenies, as well as the sensitivity of diversification analyses to different species concepts.     

Testing for unequal rates of morphological diversification in the absence of a detailed phylogeny: a case study from Characiform fishes

This study develops the random phylogenies rate test (RAPRATE), a likelihood method that simulates morphological evolution along randomly generated phylogenies, and uses it to determine whether a considerable difference in morphological diversity between two sister clades of South American fishes should be taken as evidence of differing rates of morphological change or lineage turnover. Despite identical ages of origin, similar species richness, and sympatric geographic distributions, the morphological and ecological diversity of the superfamily Anostomoidea exceeds that of the Curimatoidea. The test shows with 90% confidence (using variance among species as the measure of morphological diversity) or 99% confidence (using volume of occupied morphospace) that the rate of morphological change per unit time in the Anostomoidea likely exceeded that of the Curimatoidea. Variation in the rate of lineage turnover (speciation and extinction rates) is not found to affect greatly the morphological diversity of simulated clades and is not a likely explanation of the observed difference in morphological diversity in this case study. Though a 17% or greater delay in the onset of diversification in the Curimatoidea remains a possible alternative explanation of unequal morphological diversification, further simulations suggest that two clades drawn from the possible treespace of the Anostomoidea and Curimatoidea will rarely differ so greatly in the onset of diversification. Several uniquely derived morphological and ecological features of the Anostomoidea and Curimatoidea may have accelerated or decelerated their rate of morphological change, including a marked lengthening of the quadrate that may have relaxed structural constraints on the evolution of the anostomoid jaw.     

Phylogeny, historical biogeography, and patterns of diversification for Pinus (Pinaceae): phylogenetic tests of fossil-based hypotheses

Pines comprise one of the largest coniferous genera, are distributed throughout the Northern Hemisphere, and have an abundant fossil record. Distributions of fossils have been used to derive a three-step hypothesis of early pine evolution, which postulates a Mesozoic origin for the genus, east-west expansions across Laurasia, and retraction into Eocene refugia. Here, we present phylogenetic tests of this hypothesis using chloroplast sequence data from four loci for 83 pine species. We used the fossil-based hypothesis to derive null expectations concerning monophyly of taxonomic groups, dates of cladogenesis, and patterns of diversification. Phylogenetic analyses using several algorithms subsequently provided rigorous tests of these expectations. Our inferred phylogenies illustrated broad congruence with taxonomic groups, but highlighted consistent problems within subgenus Strobus. Estimated minimum dates of divergence derived from relaxed clock methods were largely consistent with the fossil record and yielded a date for the ingroup node of Pinus of 128+/-4 mya, depending upon the calibration used for subgenus Pinus. Ancestral area reconstructions showed Pinus to have most likely originated in Eurasia. Major clades differed in biogeographic patterns, but were consistent with the fossil-based hypothesis. We found weak support, however, for a change in diversification rate in the Eocene as interpretations of fossil distributions would have predicted.     

Understanding angiosperm diversification using small and large phylogenetic trees

How will the emerging possibility of inferring ultra-large phylogenies influence our ability to identify shifts in diversification rate? For several large angiosperm clades (Angiospermae, Monocotyledonae, Orchidaceae, Poaceae, Eudicotyledonae, Fabaceae, and Asteraceae), we explore this issue by contrasting two approaches: (1) using small backbone trees with an inferred number of extant species assigned to each terminal clade and (2) using a mega-phylogeny of 55473 seed plant species represented in GenBank. The mega-phylogeny approach assumes that the sample of species in GenBank is at least roughly proportional to the actual species diversity of different lineages, as appears to be the case for many major angiosperm lineages. Using both approaches, we found that diversification rate shifts are not directly associated with the major named clades examined here, with the sole exception of Fabaceae in the GenBank mega-phylogeny. These agreements are encouraging and may support a generality about angiosperm evolution: major shifts in diversification may not be directly associated with major named clades, but rather with clades that are nested not far within these groups. An alternative explanation is that there have been increased extinction rates in early-diverging lineages within these clades. Based on our mega-phylogeny, the shifts in diversification appear to be distributed quite evenly throughout the angiosperms. Mega-phylogenetic studies of diversification hold great promise for revealing new patterns, but we will need to focus more attention on properly specifying null expectation.     

Trees of unusual size: biased inference of early bursts from large molecular phylogenies

An early burst of speciation followed by a subsequent slowdown in the rate of diversification is commonly inferred from molecular phylogenies. This pattern is consistent with some verbal theory of ecological opportunity and adaptive radiations. One often-overlooked source of bias in these studies is that of sampling at the level of whole clades, as researchers tend to choose large, speciose clades to study. In this paper, we investigate the performance of common methods across the distribution of clade sizes that can be generated by a constant-rate birth-death process. Clades which are larger than expected for a given constant-rate branching process tend to show a pattern of an early burst even when both speciation and extinction rates are constant through time. All methods evaluated were susceptible to detecting this false signature when extinction was low. Under moderate extinction, both the [Formula: see text]-statistic and diversity-dependent models did not detect such a slowdown but only because the signature of a slowdown was masked by subsequent extinction. Some models which estimate time-varying speciation rates are able to detect early bursts under higher extinction rates, but are extremely prone to sampling bias. We suggest that examining clades in isolation may result in spurious inferences that rates of diversification have changed through time.