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.     

Species selection and random drift in macroevolution

Species selection resulting from trait‐dependent speciation and extinction is increasingly recognized as an important mechanism of phenotypic macroevolution. However, the recent bloom in statistical methods quantifying this process faces a scarcity of dynamical theory for their interpretation, notably regarding the relative contributions of deterministic versus stochastic evolutionary forces. I use simple diffusion approximations of birth‐death processes to investigate how the expected and random components of macroevolutionary change depend on phenotype‐dependent speciation and extinction rates, as can be estimated empirically. I show that the species selection coefficient for a binary trait, and selection differential for a quantitative trait, depend not only on differences in net diversification rates (speciation minus extinction), but also on differences in species turnover rates (speciation plus extinction), especially in small clades. The randomness in speciation and extinction events also produces a species‐level equivalent to random genetic drift, which is stronger for higher turnover rates. I then show how microevolutionary processes including mutation, organismic selection, and random genetic drift cause state transitions at the species level, allowing comparison of evolutionary forces across levels. A key parameter that would be needed to apply this theory is the distribution and rate of origination of new optimum phenotypes along a phylogeny.     

QUANTIFYING VARIATION IN SPECIATION AND EXTINCTION RATES WITH CLADE DATA

High‐level phylogenies are very common in evolutionary analyses, although they are often treated as incomplete data. Here, we provide statistical tools to analyze what we name “clade data,” which are the ages of clades together with their numbers of species. We develop a general approach for the statistical modeling of variation in speciation and extinction rates, including temporal variation, unknown variation, and linear and nonlinear modeling. We show how this approach can be generalized to a wide range of situations, including testing the effects of life‐history traits and environmental variables on diversification rates. We report the results of an extensive simulation study to assess the performance of some statistical tests presented here as well as of the estimators of speciation and extinction rates. These latter results suggest the possibility to estimate correctly extinction rate in the absence of fossils. An example with data on fish is presented.     

Exploring the power of Bayesian birth‐death skyline models to detect mass extinction events from phylogenies with only extant taxa

Mass extinction events (MEEs), defined as significant losses of species diversity in significantly short time periods, have attracted the attention of biologists because of their link to major environmental change. MEEs have traditionally been studied through the fossil record, but the development of birth‐death models has made it possible to detect their signature based on extant‐taxa phylogenies. Most birth‐death models consider MEEs as instantaneous events where a high proportion of species are simultaneously removed from the tree (“single pulse” approach), in contrast to the paleontological record, where MEEs have a time duration. Here, we explore the power of a Bayesian Birth‐Death Skyline (BDSKY) model to detect the signature of MEEs through changes in extinction rates under a “time‐slice” approach. In this approach, MEEs are time intervals where the extinction rate is greater than the speciation rate. Results showed BDSKY can detect and locate MEEs but that precision and accuracy depend on the phylogeny's size and MEE intensity. Comparisons of BDSKY with the single‐pulse Bayesian model, CoMET, showed a similar frequency of Type II error and neither model exhibited Type I error. However, while CoMET performed better in detecting and locating MEEs for smaller phylogenies, BDSKY showed higher accuracy in estimating extinction and speciation rates.     

Phylogenetically patterned speciation rates and extinction risks change the loss of evolutionary history during extinctions

If we are to plan conservation strategies that minimize the loss of evolutionary history through human-caused extinctions, we must understand how this loss is related to phylogenetic patterns in current extinction risks and past speciation rates. Nee & May (1997, Science 278, 692-694) showed that for a randomly evolving clade (i) a single round of random extinction removed relatively little evolutionary history, and (ii) extinction management (choosing which taxa to sacrifice) offered only marginal improvement. However, both speciation rates and extinction risks vary across lineages within real clades. We simulated evolutionary trees with phylogenetically patterned speciation rates and extinction risks (closely related lineages having similar rates and risks) and then subjected them to several biologically informed models of extinction. Increasing speciation rate variation increases the extinction-management pay-off. When extinction risks vary among lineages but are uncorrelated with speciation rates, extinction removes more history (compared with random trees), but the difference is small. When extinction risks vary and are correlated with speciation rates, history loss can dramatically increase (negative correlation) or decrease (positive correlation) with speciation rate variation. The loss of evolutionary history via human-caused extinctions may therefore be more severe, yet more manageable, than first suggested.     

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.     

Genome Size and Species Diversification

Theoretically, there are reasons to believe that large genome size should favour speciation. Several major factors contributing to genome size, such as duplications and transposable element activity have been proposed to facilitate the formation of new species. However, it is also possible that small genome size promotes speciation. For example, selection for genome reduction may be resolved in different ways in incipient species, leading to incompatibilities. Mutations and chromosomal rearrangements may also be more stably inherited in smaller genomes. Here I review the following lines of empirical evidence bearing on this question: (i) Correlations between genome size and species richness of taxa are often negative. (ii) Fossil evidence in lungfish shows that the accumulation of DNA in the genomes of this group coincided with a reduction in species diversity. (iii) Estimates of speciation interval in mammals correlate positively with genome size. (iv) Genome reductions are inferred at the base of particular species radiations and genome expansions at the base of others. (v) Insect clades that have been increasing in diversity up to the present have smaller genomes than clades that have remained stable or have decreased in diversity. The general pattern emerging from these observations is that higher diversification rates are generally found in small-genome taxa. Since diversification rates are the net effect of speciation and extinction, large genomes may thus either constrain speciation rate, increase extinction rate, or both. I argue that some of the cited examples are unlikely to be explained by extinction alone.     

Comparing the rates of speciation and extinction between phylogenetic trees

Over the past decade or so it has become increasingly popular to use reconstructed evolutionary trees to investigate questions about the rates of speciation and extinction. Although the methodology of this field has grown substantially in its sophistication in recent years, here I will take a step back to present a very simple model that is designed to investigate the relatively straightforward question of whether the tempo of diversification (speciation and extinction) differs between two or more phylogenetic trees, without attempting to attribute a causal basis to this difference. It is a likelihood method, and I demonstrate that it generally shows type I error that is close to the nominal level. I also demonstrate that parameter estimates obtained with this approach are largely unbiased. As this method can be used to compare trees of unknown relationship, it will be particularly well‐suited to problems in which a difference in diversification rate between clades is suspected, but in which these clades are not particularly closely related. As diversification methods can easily take into account an incomplete sampling fraction, but missing lineages are assumed to be missing at random, this method is also appropriate for cases in which we have hypothesized a difference in the process of diversification between two or more focal clades, but in which many unsampled groups separate the few of interest. The method of this study is by no means an attempt to replace more sophisticated models in which, for instance, diversification depends on the state of an observed or unobserved discrete or continuous trait. Rather, my intention is to provide a complementary approach for circumstances in which a simpler hypothesis is warranted and of biological interest.     

Exploring the tempo of species diversification in legumes

Whatever criteria are used to measure evolutionary success – species numbers, geographic range, ecological abundance, ecological and life history diversity, background diversification rates, or the presence of rapidly evolving clades – the legume family is one of the most successful lineages of flowering plants. Despite this, we still know rather little about the dynamics of lineage and species diversification across the family through the Cenozoic, or about the underlying drivers of diversification. There have been few attempts to estimate net species diversification rates or underlying speciation and extinction rates for legume clades, to test whether among-lineage variation in diversification rates deviates from null expectations, or to locate species diversification rate shifts on specific branches of the legume phylogenetic tree. In this study, time-calibrated phylogenetic trees for a set of species-rich legume clades – Calliandra, Indigofereae, Lupinus, Mimosa and Robinieae – and for the legume family as a whole, are used to explore how we might approach these questions. These clades are analysed using recently developed maximum likelihood and Bayesian methods to detect species diversification rate shifts and test for among-lineage variation in speciation, extinction and net diversification rates. Possible explanations for rate shifts in terms of extrinsic factors and/or intrinsic trait evolution are discussed. In addition, several methodological issues and limitations associated with these analyses are highlighted emphasizing the potential to improve our understanding of the evolutionary dynamics of legume diversification by using much more densely sampled phylogenetic trees that integrate information across broad taxonomic, geographical and temporal levels.     

Statistical analysis of diversification with species traits

Testing whether some species traits have a significant effect on diversification rates is central in the assessment of macroevolutionary theories. However, we still lack a powerful method to tackle this objective. I present a new method for the statistical analysis of diversification with species traits. The required data are observations of the traits on recent species, the phylogenetic tree of these species, and reconstructions of ancestral values of the traits. Several traits, either continuous or discrete, and in some cases their interactions, can be analyzed simultaneously. The parameters are estimated by the method of maximum likelihood. The statistical significance of the effects in a model can be tested with likelihood ratio tests. A simulation study showed that past random extinction events do not affect the Type I error rate of the tests, whereas statistical power is decreased, though some power is still kept if the effect of the simulated trait on speciation is strong. The use of the method is illustrated by the analysis of published data on primates. The analysis of these data showed that the apparent overall positive relationship between body mass and species diversity is actually an artifact due to a clade-specific effect. Within each clade the effect of body mass on speciation rate was in fact negative. The present method allows to take both effects (clade and body mass) into account simultaneously.     

Phylogenetic community structure metrics and null models: a review with new methods and software

Competitive exclusion and habitat filtering influence community assembly, but ecologists and evolutionary biologists have not reached consensus on how to quantify patterns that would reveal the action of these processes. Currently, at least 22 α‐diversity and 10 β‐diversity metrics of community phylogenetic structure can be combined with nine null models (eight for β‐diversity metrics), providing 278 potentially distinct approaches to test for phylogenetic clustering and overdispersion. Selecting the appropriate approach for a study is daunting. First, we describe similarities among metrics and null models across variance in phylogeny size and shape, species abundance, and species richness. Second, we develop spatially explicit, individual‐based simulations of neutral, competitive exclusion, or habitat filtering community assembly, and quantify the performance (type I and II error rates) of all 278 metric and null model combinations against each assembly process. Many α‐diversity metrics and null models are at least functionally equivalent, reducing the number of truly unique metrics to 12 and the number of unique metric + null model combinations to 72. An even smaller subset of metric and null model combinations showed robust statistical performance. For α‐diversity metrics, phylogenetic diversity and mean nearest taxon distance were best able to detect habitat filtering, while mean pairwise phylogenetic distance‐based metrics were best able to detect competitive exclusion. Overall, β‐diversity metrics tended to have greater power to detect habitat filtering and competitive exclusion than α‐diversity metrics, but had higher type 1 error in some cases. Across both α‐ and β‐diversity metrics, null model selection affected type I error rates more than metric selection. A null model that maintained species richness, and approximately maintained species occurrence frequency and abundance across sites, exhibited low type I and II error rates. This regional null model simulates neutral dispersal of individuals into local communities by sampling from a regional species pool. We introduce a flexible new R package, metricTester, to facilitate robust analyses of method performance.     

Detection of "punctuated equilibrium" by bayesian estimation of speciation and extinction rates, ancestral character states, and rates of anagenetic and cladogenetic evolution on a molecular phylogeny

Algorithms are presented to simultaneously estimate probabilities of speciation and extinction, rates of anagenetic and cladogenetic phenotypic evolution, as well as ancestral character states, from a complete ultrametric species-level phylogeny with dates assigned to all bifurcations and one or more phenotypes in three or more extant species, using Metropolis-Hastings Markov Chain Monte Carlo sampling. The algorithms also estimate missing phenotypes of extant species and numbers of speciation events that occurred on all branches of the phylogeny. The algorithms are discussed and their performance is evaluated using simulated data. That evaluation shows that precise estimation of rates of evolution of one or a few phenotypes requires large phylogenies. Estimation accuracy improves with the number of species on the phylogeny.     

SEARCHING FOR EVOLUTIONARY PATTERNS IN THE SHAPE OF A PHYLOGENETIC TREE

If all species in a clade are equally likely to speciate or become extinct, then highly symmetric and highly asymmetric phylogenetic trees are unlikely to result. Variation between species in speciation and extinction rates can cause excessive asymmetry. We developed six non-parametric statistical tests that test for nonrandom patterns of branching in any bifurcating tree. The tests are demonstrated by applying them to two published phylogenies for genera of beetles. Comparison of the power of the six statistics under a simple model of biased speciation suggests which of them may be most useful for detecting nonrandom tree shapes.     

Santa Rosalia revisited: Why are there so many species of bacteria?

The diversity of bacteria in the world is very poorly known. Usually less than one percent of the bacteria from natural communities can be grown in the laboratory. This has caused us to underestimate bacterial diversity and biased our view of bacterial communities. The tools are now available to estimate the number of bacterial species in a community and to estimate the difference between communities. Using what data are available, I have estimated that thirty grams of forest soil contains over half a million species. The species difference between related communities suggests that the number of species of bacteria may be more than a thousand million. I suppose that the explanation for such a large number of bacterial species is simply that speciation in bacteria is easy and extinction difficult, giving a rate of speciation higher than the rate of extinction, leading to an ever increasing number of species over time. The idea that speciation is easy is justified from the results of recent experimental work in bacterial evolution.     

Does coevolution promote species richness in parasitic cuckoos?

Why some lineages have diversified into larger numbers of species than others is a fundamental but still relatively poorly understood aspect of the evolutionary process. Coevolution has been recognized as a potentially important engine of speciation, but has rarely been tested in a comparative framework. We use a comparative approach based on a complete phylogeny of all living cuckoos to test whether parasite–host coevolution is associated with patterns of cuckoo species richness. There are no clear differences between parental and parasitic cuckoos in the number of species per genus. However, a cladogenesis test shows that brood parasitism is associated with both significantly higher speciation and extinction rates. Furthermore, subspecies diversification rate estimates were over twice as high in parasitic cuckoos as in parental cuckoos. Among parasitic cuckoos, there is marked variation in the severity of the detrimental effects on host fitness; chicks of some cuckoo species are raised alongside the young of the host and others are more virulent, with the cuckoo chick ejecting or killing the eggs/young of the host. We show that cuckoos with a more virulent parasitic strategy have more recognized subspecies. In addition, cuckoo species with more recognized subspecies have more hosts. These results hold after controlling for confounding geographical effects such as range size and isolation in archipelagos. Although the power of our analyses is limited by the fact that brood parasitism evolved independently only three times in cuckoos, our results suggest that coevolutionary arms races with hosts have contributed to higher speciation and extinction rates in parasitic cuckoos.     

The ecological dynamics of clade diversification and community assembly

Clades diversify in an ecological context, but most macroevolutionary models do not directly encapsulate ecological mechanisms that influence speciation and extinction. A data set of 245 chordate, arthropod, mollusk, and magnoliophyte phylogenies had a majority of clades that showed rapid lineage accumulation early with a slowing more recently, whereas a small but significant minority showed accelerated lineage accumulation in their recent histories. Previous analyses have demonstrated that macroevolutionary birth-death models can replicate the pattern of slowing lineage accumulation only by a strong decrease in speciation rate with increasing species richness and extinction rate held extremely low or absent. In contrast, the metacommunity model presented here could generate the full range of patterns seen in the real phylogenies by simply manipulating the degree of ecological differentiation of new species at the time of speciation. Specifically, the metacommunity model predicts that clades showing decelerating lineage accumulation rates are those that have diversified by ecological modes of speciation, whereas clades showing accelerating lineage accumulation rates are those that have diversified primarily by modes of speciation that generate little or no ecological diversification. A number of testable predictions that integrate data from molecular systematics, community ecology, and biogeography are also discussed.     

Latitude and rates of diversification in birds and butterflies

Central to many explanations of latitudinal diversity gradients is the idea that rates of species diversification increase towards the equator. However, there have been few explicit tests of whether or not this pattern exists. Using sister-group analyses to compare 48 clades of passerine birds and swallowtail butterflies from different latitudes, I found evidence that relative rates of diversification per unit time are indeed higher towards the equator. This pattern is explicable in terms of abiotic factors which vary continuously with latitude, and may be further enhanced by diversity-dependent speciation and extinction processes.     

Settling down of seasonal migrants promotes bird diversification

How seasonal migration originated and impacted diversification in birds remains largely unknown. Although migratory behaviour is likely to affect bird diversification, previous studies have not detected any effect. Here, we infer ancestral migratory behaviour and the effect of seasonal migration on speciation and extinction dynamics using a complete bird tree of life. Our analyses infer that sedentary behaviour is ancestral, and that migratory behaviour evolved independently multiple times during the evolutionary history of birds. Speciation of a sedentary species into two sedentary daughter species is more frequent than speciation of a migratory species into two migratory daughter species. However, migratory species often diversify by generating a sedentary daughter species in addition to the ancestral migratory one. This leads to an overall higher migratory speciation rate. Migratory species also experience lower extinction rates. Hence, although migratory species represent a minority (18.5%) of all extant birds, they have a higher net diversification rate than sedentary species. These results suggest that the evolution of seasonal migration in birds has facilitated diversification through the divergence of migratory subpopulations that become sedentary, and illustrate asymmetrical diversification as a mechanism by which diversification rates are decoupled from species richness.     

Simulating trees with a fixed number of extant species

In this paper, I develop efficient tools to simulate trees with a fixed number of extant species. The tools are provided in my open source R-package TreeSim available on CRAN. The new model presented here is a constant rate birth-death process with mass extinction and/or rate shift events at arbitrarily fixed times 1) before the present or 2) after the origin. The simulation approach for case (2) can also be used to simulate under more general models with fixed events after the origin. I use the developed simulation tools for showing that a mass extinction event cannot be distinguished from a model with constant speciation and extinction rates interrupted by a phase of stasis based on trees consisting of only extant species. However, once we distinguish between mass extinction and period of stasis based on paleontological data, fast simulations of trees with a fixed number of species allow inference of speciation and extinction rates using approximate Bayesian computation and allow for robustness analysis once maximum likelihood parameter estimations are available.     

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.