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.     

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.     

The shape and temporal dynamics of phylogenetic trees arising from geographic speciation

Phylogenetic trees often depart from the expectations of stochastic models, exhibiting imbalance in diversification among lineages and slowdowns in the rate of lineage accumulation through time. Such departures have led to a widespread perception that ecological differences among species or adaptation and subsequent niche filling are required to explain patterns of diversification. However, a key element missing from models of diversification is the geographical context of speciation and extinction. In this study, we develop a spatially explicit model of geographic range evolution and cladogenesis, where speciation arises via vicariance or peripatry, and explore the effects of these processes on patterns of diversification. We compare the results with those observed in 41 reconstructed avian trees. Our model shows that nonconstant rates of speciation and extinction are emergent properties of the apportioning of geographic ranges that accompanies speciation. The dynamics of diversification exhibit wide variation, depending on the mode of speciation, tendency for range expansion, and rate of range evolution. By varying these parameters, the model is able to capture many, but not all, of the features exhibited by birth-death trees and extant bird clades. Under scenarios with relatively stable geographic ranges, strong slowdowns in diversification rates are produced, with faster rates of range dynamics leading to constant or accelerating rates of apparent diversification. A peripatric model of speciation with stable ranges also generates highly unbalanced trees typical of bird phylogenies but fails to produce realistic range size distributions among the extant species. Results most similar to those of a birth-death process are reached under a peripatric speciation scenario with highly volatile range dynamics. Taken together, our results demonstrate that considering the geographical context of speciation and extinction provides a more conservative null model of diversification and offers a very different perspective on the phylogenetic patterns expected in the absence of ecology.     

The molecular phylogenetic signature of clades in decline

Molecular phylogenies have been used to study the diversification of many clades. However, current methods for inferring diversification dynamics from molecular phylogenies ignore the possibility that clades may be decreasing in diversity, despite the fact that the fossil record shows this to be the case for many groups. Here we investigate the molecular phylogenetic signature of decreasing diversity using the most widely used statistic for inferring diversity dynamics from molecular phylogenies, the γ statistic. We show that if a clade is in decline its molecular phylogeny may show evidence of the decrease in the diversification rate that occurred between its diversification and decline phases. The ability to detect the change in diversification rate depends largely on the ratio of the speciation rates of the diversification and decline phases, the higher the ratio the stronger the signal of the change in diversification rate. Consequently, molecular phylogenies of clades in relative rapid decline do not carry a signature of their decreasing diversification. Further, the signal of the change in diversification rate, if present, declines as the diversity drop. Unfortunately, the molecular signature of clades in decline is the same as the signature produced by diversity dependent diversification. Given this similarity, and the inability of current methods to detect declining diversity, it is likely that some of the extant clades that show a decrease in diversification rate, currently interpreted as evidence for diversity dependent diversification, are in fact in decline. Unless methods can be developed that can discriminate between the different modes of diversification, specifically diversity dependent diversification and declining diversity, we will need the fossil record, or data from some other source, to distinguish between these very different diversity trajectories.     

Prolonging the past counteracts the pull of the present: protracted speciation can explain observed slowdowns in diversification

Phylogenetic trees show a remarkable slowdown in the increase of number of lineages towards the present, a phenomenon which cannot be explained by the standard birth-death model of diversification with constant speciation and extinction rates. The birth-death model instead predicts a constant or accelerating increase in the number of lineages, which has been called the pull of the present. The observed slowdown has been attributed to nonconstancy of the speciation and extinction rates due to some form of diversity dependence (i.e., species-level density dependence), but the mechanisms underlying this are still unclear. Here, we propose an alternative explanation based on the simple concept that speciation takes time to complete. We show that this idea of "protracted" speciation can be incorporated in the standard birth-death model of diversification. The protracted birth-death model predicts a realistic slowdown in the rate of increase of number of lineages in the phylogeny and provides a compelling fit to four bird phylogenies with realistic parameter values. Thus, the effect of recognizing the generally accepted fact that speciation is not an instantaneous event is significant; even if it cannot account for all the observed patterns, it certainly contributes substantially and should therefore be incorporated into future studies.     

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.     

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.     

ECOLOGICAL CAUSES OF DECELERATING DIVERSIFICATION IN CARNIVORAN MAMMALS

Clade diversification is a central topic in macroevolutionary studies. Recently, it has been shown that diversification rates appear to decelerate over time in many clades. What causes this deceleration remains unclear, but it has been proposed that competition for limited resources between sympatric, ecologically similar species slows diversification. Employing carnivoran mammals as a model system, we test this hypothesis using a comprehensive time‐calibrated phylogeny. We also explore several conceptually related explanations including limited geographic area and limited rates of niche evolution. We find that diversification slowdowns are strong in carnivorans. Surprisingly, these slowdowns are independent of geographic range overlap between related species and are also decoupled from rates of niche evolution, suggesting that slowdowns are unrelated to competition and niche filling. When controlling for the effects of clade diversity, diversification slowdowns appear independent of geographic area. There is a significant effect of clade diversity on diversification slowdowns, but simulations show that this relationship may arise as a statistical artifact (i.e., greater clade diversity increases the ability of the gamma statistic to refute constant diversification). Overall, our results emphasize the need to test hypotheses about the causes of diversification slowdowns with ecological data, rather than assuming ecological processes from phylogenies alone.     

HOW DOES ECOLOGICAL OPPORTUNITY INFLUENCE RATES OF SPECIATION,EXTINCTION, AND MORPHOLOGICAL DIVERSIFICATION IN NEW WORLD RATSNAKES (TRIBE LAMPROPELTINI)?

Ecological adaptive radiation theory predicts an increase in both morphological and specific diversification when organisms colonize new environments. Accordingly, bursts of morphological diversification, characterized by low within‐subclade morphological disparity, may be associated with these increases in speciation rates. Conversely, increasing species density, reduction in available habitat, or increasing extinction rates are expected to cause rates of diversification to decline. We test these hypotheses by examining the tempo and mode of speciation in the lampropeltinine snakes, a morphologically variable group that colonized the New World ～24 million years ago and radiated throughout the Miocene. We show that specific diversification increased early in the history of the group, and that most morphological variation is partitioned among, rather than within subclades. These patterns provide further evidence for the hypothesis that morphological variation tends to be strongly partitioned among lineages when clades undergo early bursts of species diversification. A reduction in speciation rates may be indicative of density dependent effects due to a saturation of available ecological opportunity, rather than increases in extinction rates at the onset of the Pleistocene/Pliocene glacial cycles. This evidence runs counter to the general Pleistocene species pump model.     

The shape of mammalian phylogeny: patterns, processes and scales

Mammalian phylogeny is far too asymmetric for all contemporaneous lineages to have had equal chances of diversifying. We consider this asymmetry or imbalance from four perspectives. First, we infer a minimal set of 'regime changes'-points at which net diversification rate has changed-identifying 15 significant radiations and 12 clades that may be 'downshifts'. We next show that mammalian phylogeny is similar in shape to a large set of published phylogenies of other vertebrate, arthropod and plant groups, suggesting that many clades may diversify under a largely shared set of 'rules'. Third, we simulate six simple macroevolutionary models, showing that those where speciation slows down as geographical or niche space is filled, produce more realistic phylogenies than do models involving key innovations. Lastly, an analysis of the spatial scaling of imbalance shows that the phylogeny of species within an assemblage, ecoregion or larger area always tends to be more unbalanced than expected from the phylogeny of species at the next more inclusive spatial scale. We conclude with a verbal model of mammalian macroevolution, which emphasizes the importance to diversification of accessing new regions of geographical or niche space.     

Ecological opportunity alters the timing and shape of adaptive radiation

The uneven distribution of diversity is a conspicuous phenomenon across the tree of life. Ecological opportunity is a prominent catalyst of adaptive radiation and therefore may alter patterns of diversification. We evaluated the distribution of shifts in diversification rates across the cichlid phylogeny and the distribution of major clades across phylogenetic space. We also tested if ecological opportunity influenced these patterns. Colonization‐associated ecological opportunity altered the tempo and mode of diversification during the adaptive radiation of cichlid fishes. Clades that arose following colonization events diversified faster than other clades. Speciation rate shifts were nonrandomly distributed across the phylogeny such that they were disproportionally concentrated around nodes that corresponded with colonization events (i.e., of continents, river basins, or lakes). Young clades tend to expand faster than older clades; however, colonization‐associated ecological opportunity accentuated this pattern. There was an interaction between clade age and ecological opportunity that explained the trajectory of clades through phylogenetic space over time. Our results indicate that ecological opportunities afforded by continental and ecosystem‐scale colonization events explain the dramatic speciation rate heterogeneity and phylogenetic imbalance that arose during the evolutionary history of cichlid fishes.     

Rapid speciation and ecological divergence in the American seven-spined gobies (Gobiidae,Gobiosomatini) inferred from a molecular phylogeny

Abstract.— The American seven-spined gobies (Gobiidae, Gobiosomatini) are highly diverse both in morphology and ecology with many endemics in the Caribbean region. We have reconstructed a molecular phylogeny of 54 Gobio-somatini taxa (65 individuals) based on a 1646-bp region that includes the mitochondrial 12S rRNA, tRNA-Val, and 16S rRNA genes. Our results support the monophyly of the seven-spined gobies and are in agreement with the existence of two major groups within the tribe, the Gobiosoma group and the Microgobius group. However, they reject the monophyly of some of the Gobiosomatini genera. We use the molecular phylogeny to study the dynamics of speciation in the Gobiosomatini by testing for departures from the constant speciation rate model. We observe a burst of speciation in the early evolutionary history of the group and a subsequent slowdown. Our results show a split among clades into coastal-estuarian, deep ocean, and tropical reef habitats. Major habitat shifts account for the early significant acceleration in lineage splitting and speciation rate and the initial divergence of the main Gobiosomatini clades. We found that subsequent diversification is triggered by behavior and niche specializations at least in the reef-associated clades. Overall, our results confirm that the diversity of Gobiosomatini has arisen during episodes of adaptive radiation, and emphasize the importance of ecology in marine speciation.     

Patterns of diversification of Afrotropical Otiteselline fig wasps: phylogenetic study reveals a double radiation across host figs and conservatism of host association

We studied the phylogenetic relationships of Otiteselline fig waSPS associated with Ficus in the Afrotropical region using rDNA sequences. African fig species usually host two species of Otiteselline fig waSPS. Phylogenetic analyses reveal that this pattern of association results from the radiation of two clades of waSPS superimposed on the fig system. Within each clade, wasp species generally cluster according to their host classification. The phylogenies of the two clades are also generally more congruent than expected by chance. Together these results suggest that Otiteselline wasp speciation is largely constrained by the diversification of their hosts. Finally, we show a difference in ovipositor length between the two Otiteselline species coexisting in the same Ficus species, which probably corresponds to ecological differences. The diversification of ecological niches within the fig is probably, with cospeciation, one of the key factors explaining the diversification and maintenance of species of parasites of the fig/pollinator system.     

Diversification rates have declined in the Malagasy herpetofauna

The evolutionary origins of Madagascar''s biodiversity remain mysterious despite the fact that relative to land area, there is no other place with consistently high levels of species richness and endemism across a range of taxonomic levels. Most efforts to explain diversification on the island have focused on geographical models of speciation, but recent studies have begun to address the island''s accumulation of species through time, although with conflicting results. Prevailing hypotheses for diversification on the island involve either constant diversification rates or scenarios where rates decline through time. Using relative-time-calibrated phylogenies for seven endemic vertebrate clades and a model-fitting framework, I find evidence that diversification rates have declined through time on Madagascar. I show that diversification rates have clearly declined throughout the history of each clade, and models invoking diversity-dependent reductions to diversification rates best explain the diversification histories for each clade. These results are consistent with the ecological theory of adaptive radiation, and, coupled with ancillary observations about ecomorphological and life-history evolution, strongly suggest that adaptive radiation was an important formative process for one of the most species-rich regions on the Earth. These results cast the Malagasy biota in a new light and provide macroevolutionary justification for conservation initiatives.     

Diversity versus disparity and the radiation of modern cetaceans

Modern whales are frequently described as an adaptive radiation spurred by either the evolution of various key innovations (such as baleen or echolocation) or ecological opportunity following the demise of archaic whales. Recent analyses of diversification rate shifts on molecular phylogenies raise doubts about this interpretation since they find no evidence of increased speciation rates during the early evolution of modern taxa. However, one of the central predictions of ecological adaptive radiation is rapid phenotypic diversification, and the tempo of phenotypic evolution has yet to be quantified in cetaceans. Using a time-calibrated molecular phylogeny of extant cetaceans and a morphological dataset on size, we find evidence that cetacean lineages partitioned size niches early in the evolutionary history of neocetes and that changes in cetacean size are consistent with shifts in dietary strategy. We conclude that the signature of adaptive radiations may be retained within morphological traits even after equilibrium diversity has been reached and high extinction or fluctuations in net diversification have erased any signature of an early burst of diversification in the structure of the phylogeny.     

ESTIMATING THE DURATION OF SPECIATION FROM PHYLOGENIES

Speciation is not instantaneous but takes time. The protracted birth–death diversification model incorporates this fact and predicts the often observed slowdown of lineage accumulation toward the present. The mathematical complexity of the protracted speciation model has barred estimation of its parameters until recently a method to compute the likelihood of phylogenetic branching times under this model was outlined (Lambert et al. 2014 ). Here, we implement this method and study using simulated phylogenies of extant species how well we can estimate the model parameters (rate of initiation of speciation, rate of extinction of incipient and good species, and rate of completion of speciation) as well as the duration of speciation, which is a combination of the aforementioned parameters. We illustrate our approach by applying it to a primate phylogeny. The simulations show that phylogenies often do not contain enough information to provide unbiased estimates of the speciation‐initiation rate and the extinction rate, but the duration of speciation can be estimated without much bias. The estimate of the duration of speciation for the primate clade is consistent with literature estimates. We conclude that phylogenies combined with the protracted speciation model provide a promising way to estimate the duration of speciation.     

Evolutionary diversification of clades of squamate reptiles

We analysed the diversification of squamate reptiles (7488 species) based on a new molecular phylogeny, and compared the results to similar estimates for passerine birds (5712 species). The number of species in each of 36 squamate lineages showed no evidence of phylogenetic conservatism. Compared with a random speciation-extinction process with parameters estimated from the size distribution of clades, the alethinophidian snakes (2600 species) were larger than expected and 13 clades, each having fewer than 20 species, were smaller than expected, indicating rate heterogeneity. From a lineage-through-time plot, we estimated that a provisional rate of lineage extinction (0.66 per Myr) was 94% of the rate of lineage splitting (0.70 per Myr). Diversification in squamate lineages was independent of their stem age, but strongly related to the area of the region within which they occur. Tropical vs. temperate latitude exerted a marginally significant influence on species richness. In comparison with passerine birds, squamates share several clade features, including: (1) independence of species richness and age; (2) lack of phylogenetic signal with respect to clade size; (3) general absence of exceptionally large clades; (4) over-representation of small clades; (5) influence of region size on clade size; and (6) similar rates of speciation and extinction. The evidence for both groups suggests that clade size has achieved long-term equilibrium, suggesting negative feedback of species richness on the rate of diversification.     

Rapid allopatric speciation in logperch darters (Percidae: Percina)

Theory predicts that clades diversifying via sympatric speciation will exhibit high diversification rates. However, the expected rate of diversification in clades characterized by allopatric speciation is less clear. Previous studies have documented significantly higher speciation rates in freshwater fish clades diversifying via sympatric versus allopatric modes, leading to suggestions that the geographic pattern of speciation can be inferred solely from knowledge of the diversification rate. We tested this prediction using an example from darters, a clade of approximately 200 species of freshwater fishes endemic to eastern North America. A resolved phylogeny was generated using mitochondrial DNA gene sequences for logperches, a monophyletic group of darters composed of 10 recognized species. Divergence times among logperch species were estimated using a fossil calibrated molecular clock in centrarchid fishes, and diversification rates in logperches were estimated using several methods. Speciation events in logperches are recent, extending from 4.20 +/- 1.06 million years ago (mya) to 0.42 +/- 0.22 mya, with most speciation events occurring in the Pleistocene. Diversification rates are high in logperches, at some nodes exceeding rates reported for well-studied adaptive radiations such as Hawaiian silverswords. The geographic pattern of speciation in logperches was investigated by examining the relationship between degree of sympatry and the absolute age of the contrast, with the result that diversification in logperches appears allopatric. The very high diversification rate observed in the logperch phylogeny is more similar to freshwater fish clades thought to represent examples of sympatric speciation than to clades representing allopatric speciation. These results demonstrate that the geographic mode of speciation for a clade cannot be inferred from the diversification rate. The empirical observation of high diversification rates in logperches demonstrates that allopatric speciation can occur rapidly.     

Renewed diversification is associated with new ecological opportunity in the Neotropical turtle ants

Ecological opportunity, defined as access to new resources free from competitors, is thought to be a catalyst for the process of adaptive radiation. Much of what we know about ecological opportunity, and the larger process of adaptive radiation, is derived from vertebrate diversification on islands. Here, we examine lineage diversification in the turtle ants (Cephalotes), a species‐rich group of ants that has diversified throughout the Neotropics. We show that crown group turtle ants originated during the Eocene (around 46 mya), coincident with global warming and the origin of many other clades. We also show a marked lineage‐wide slowdown in diversification rates in the Miocene. Contrasting this overall pattern, a species group associated with the young and seasonally harsh Chacoan biogeographic region underwent a recent burst of diversification. Subsequent analyses also indicated that there is significant phylogenetic clustering within the Chacoan region and that speciation rates are highest there. Together, these findings suggest that recent ecological opportunity, from successful colonization of novel habitat, may have facilitated renewed turtle ant diversification. Our findings highlight a central role of ecological opportunity within a successful continental radiation.     

Assessing the causes of diversification slowdowns: temperature‐dependent and diversity‐dependent models receive equivalent support

Diversification rates vary over time, yet the factors driving these variations remain unclear. Temporal declines in speciation rates have often been interpreted as the effect of ecological limits, competition, and diversity dependence, emphasising the role of biotic factors. Abiotic factors, such as climate change, are also supposed to have affected diversification rates over geological time scales, yet direct tests of these presumed effects have mainly been limited to few clades well represented in the fossil record. If warmer climatic periods have sustained faster speciation, this could explain slowdowns in speciation during the Cenozoic climate cooling. Here, we apply state‐of‐the art diversity‐dependent and temperature‐dependent phylogenetic models of diversification to 218 tetrapod families, along with constant rate and time‐dependent models. We confirm the prevalence of diversification slowdowns, and find as much support for temperature‐dependent than diversity‐dependent models. These results call for a better integration of these two processes in studies of diversification dynamics.