Facilitation can increase the phylogenetic diversity of plant communities

With the advent of molecular phylogenies the assessment of community assembly processes has become a central topic in community ecology. These processes have focused almost exclusively on habitat filtering and competitive exclusion. Recent evidence, however, indicates that facilitation has been important in preserving biodiversity over evolutionary time, with recent lineages conserving the regeneration niches of older, distant lineages. Here we test whether, if facilitation among distant-related species has preserved the regeneration niche of plant lineages, this has increased the phylogenetic diversity of communities. By analyzing a large worldwide database of species, we showed that the regeneration niches were strongly conserved across evolutionary history. Likewise, a phylogenetic supertree of all species of three communities driven by facilitation showed that nurse species facilitated distantly related species and increased phylogenetic diversity.     

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.     

Is regional species diversity bounded or unbounded?

Two conflicting hypotheses have been proposed to explain large‐scale species diversity patterns and dynamics. The unbounded hypothesis proposes that regional diversity depends only on time and diversification rate and increases without limit. The bounded hypothesis proposes that ecological constraints place upper limits on regional diversity and that diversity is usually close to its limit. Recent evidence from the fossil record, phylogenetic analysis, biogeography, and phenotypic disparity during lineage diversification suggests that diversity is constrained by ecological processes but that it is rarely asymptotic. Niche space is often unfilled or can be more finely subdivided and still permit coexistence, and new niche space is often created before ecological limits are reached. Damped increases in diversity over time are the prevalent pattern, suggesting the need for a new ‘damped increase hypothesis'. The damped increase hypothesis predicts that diversity generally increases through time but that its rate of increase is often slowed by ecological constraints. However, slowing due to niche limitation must be distinguished from other possible mechanisms creating similar patterns. These include sampling artifacts, the inability to detect extinctions or declines in clade diversity with some methods, the distorting effects of correlated speciation‐extinction dynamics, the likelihood that opportunities for allopatric speciation will vary in space and time, and the role of undetected natural enemies in reducing host ranges and thus slowing speciation rates. The taxonomic scope of regional diversity studies must be broadened to include all ecologically similar species so that ecological constraints may be accurately inferred. The damped increase hypothesis suggests that information on evolutionary processes such as time‐for‐speciation and intrinsic diversification rates as well as ecological factors will be required to explain why regional diversity varies among times, places and taxa.     

Diversification and biogeographic patterns in four island radiations of passerine birds

Declining diversification rates over time are a well-established evolutionary pattern, often interpreted as indicating initial rapid radiation with filling of ecological niche space. Here, we test the hypothesis that island radiations may show constant net diversification rates over time, due to continued expansion into new niche space in highly dispersive taxa. We investigate diversification patterns of four passerine bird families originating from the Indo-Pacific archipelagos, and link these to biogeographic patterns to provide independent indications of niche filling. We find a declining diversification rate for only one family, the Paradisaeidae (41 species). These are almost completely restricted to New Guinea, and have on average smaller species ranges and higher levels of species richness within grid cells than the other three families. In contrast, we cannot reject constant diversification rates for Campephagidae (93 species), Oriolidae (35 species), and Pachycephalidae (53 species), groups that have independently colonized neighboring archipelagos and continents. We propose that Paradisaeidae have reached the diversity limit imposed by their restricted distribution, whereas high dispersal and colonization success across the geologically dynamic Indo-Pacific archipelagos may have sustained high speciation rates for the other three families. Alternatively, increasing extinction rates may have obscured declining speciation rates in those three phylogenies.     

How phenotypic matching based on neutral mating cues enables speciation in locally adapted populations

Maynard Smith's (American Naturalist, 1966, 100, 637) suggestion that in some cases a prerequisite for speciation is the existence of local ecological adaptations has not received much attention to date. Here, we test the hypothesis using a model like that of Maynard Smith but differing in the way animals disperse between niches. In previous studies, males disperse randomly between niches but females stay put in their natal niche. As a first step toward generalizing the model, we here analyze the case that equal proportions of the two sexes disperse between niches before breeding. Supporting Maynard Smith's (1966) hypothesis, we find that once local adaptations are established, a neutral mating cue at an independent locus can rapidly enable speciation in populations with a suitable mechanism for phenotype matching. We find that stable ecological polymorphisms are relatively insensitive to the strength of selection, but depend crucially on the extent of dispersal between niches, with a threshold of ~5% if population sizes in two niches are equal. At higher levels of dispersal, ecological differentiation is lost. These results contrast with those of earlier studies and shed light on why parapatric speciation is limited by the extent of gene flow. Our testable model provides a candidate explanation for the rapid speciation rates, diversity of appearance and occurrence of “species flocks” observed among some African cichlids and neotropical birds and may also have implications for the occurrence of punctuational change on phylogenies.     

Niche width impacts vertebrate diversification

Aim The size of the climatic niche of a species is a major factor determining its distribution and evolution. In particular, it has been proposed that niche width should be associated with the rate of species diversification. Here, we test whether species niche width affects the speciation and extinction rates of three main clades of vertebrates: amphibians, mammals and birds. Location Global. Methods We obtained the time‐calibrated phylogenies, IUCN conservation status, species distribution maps and climatic data for 2340 species of amphibians, 4563 species of mammals and 9823 species of birds. We computed the niche width for each species as the mean annual temperature across the species range. We estimated speciation, extinction and transition rates associated with lineages with either narrow (specialist) or wide (generalist) niches using phylogeny‐based birth–death models. We also tested if current conservation status was correlated with the niche width of species. Results We found higher net diversification rates in specialist species than in generalist species. This result was explained by both higher speciation rates (for the three taxonomic groups) and lower extinction rates (for mammals and birds only) in specialist than in generalist species. In contrast, current specialist species tended to be more threatened than generalist species. Main conclusions Our diversification analysis shows that the width of the climatic niche is strongly associated with diversification rates and may thus be a crucial factor for understanding the emergence of diversity patterns in vertebrates. The striking difference between our diversification results and current conservation status suggests that the current extinction process may be different from extinction rates estimated from the whole history of the group.     

Positive correlation between diversification rates and phenotypic evolvability can mimic punctuated equilibrium on molecular phylogenies

The hypothesis of punctuated equilibrium proposes that most phenotypic evolution occurs in rapid bursts associated with speciation events. Several methods have been developed that can infer punctuated equilibrium from molecular phylogenies in the absence of paleontological data. These methods essentially test whether the variance in phenotypes among extant species is better explained by evolutionary time since common ancestry or by the number of estimated speciation events separating taxa. However, apparent "punctuational" trait change can be recovered on molecular phylogenies if the rate of phenotypic evolution is correlated with the rate of speciation. Strong support for punctuational models can arise even if the underlying mode of trait evolution is strictly gradual, so long as rates of speciation and trait evolution covary across the branches of phylogenetic trees, and provided that lineages vary in their rate of speciation. Species selection for accelerated rates of ecological or phenotypic divergence can potentially lead to the perception that most trait divergence occurs in association with speciation events.     

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.     

Modeling colonization rates over time: Generating null models and testing model adequacy in phylogenetic analyses of species assemblages*

A central theme connecting macroevolutionary processes to macroecological patterns is the shaping of regional biodiversity over time through speciation, extinction, migration, and range shifts. The use of phylogenies to explore the dynamics of diversification due to variation in speciation and extinction rates has been well-developed and there are established methods for inferring speciation times from phylogenies and generating its null distributions (as represented by node heights on molecular phylogenies). But inferring colonization events from phylogenies is more challenging. Unlike speciation events, represented by nodes, colonization events could occur at any point along a branch connecting species in the assemblage to the regional pool. We account for uncertainty in identification of colonization lineages and timing of colonization events by using an efficient analytical solution to inferring the distribution of colonization times from an assemblage phylogeny. Using the same solution, we efficiently derive the null distribution of colonization times, which provides us with a general approach to testing the adequacy of a model to describe colonization events into the assemblage. We illustrate this approach by demonstrating how the movement of squamate lineages into Madagascar has been uneven over time, peaking in the early Cenozoic when ocean conditions favored colonization.     

Testing macro-evolutionary models using incomplete molecular phylogenies.

Phylogenies reconstructed from gene sequences can be used to investigate the tempo and mode of species diversification. Here we develop and use new statistical methods to infer past patterns of speciation and extinction from molecular phylogenies. Specifically, we test the null hypothesis that per-lineage speciation and extinction rates have remained constant through time. Rejection of this hypothesis may provide evidence for evolutionary events such as adaptive radiations or key adaptations. In contrast to previous approaches, our methods are robust to incomplete taxon sampling and are conservative with respect to extinction. Using simulation we investigate, first, the adverse effects of failing to take incomplete sampling into account and, second, the power and reliability of our tests. When applied to published phylogenies our tests suggest that, in some cases, speciation rates have decreased through time.     

Investigating the timing of origin and evolutionary processes shaping regional species diversity: Insights from simulated data and neotropical butterfly diversification rates

Different diversification scenarios have been proposed to explain the origin of extant biodiversity. However, most existing meta‐analyses of time‐calibrated phylogenies rely on approaches that do not quantitatively test alternative diversification processes. Here, I highlight the shortcomings of using species divergence ranks, which is a method widely used in meta‐analyses. Divergence ranks consist of categorizing cladogenetic events to certain periods of time, typically to either Pleistocene or to pre‐Pleistocene ages. This approach has been claimed to shed light on the origin of most extant species and the timing and dynamics of diversification in any biogeographical region. However, interpretations drawn from such method often confound two fundamental questions in macroevolutionary studies, tempo (timing of evolutionary rate shifts) and mode (“how” and “why” of speciation). By using simulated phylogenies under four diversification scenarios, constant‐rate, diversity‐dependence, high extinction, and high speciation rates in the Pleistocene, I showed that interpretations based on species divergence ranks might have been seriously misleading. Future meta‐analyses of dated phylogenies need to be aware of the impacts of incomplete taxonomic sampling, tree topology, and divergence time uncertainties, as well as they might be benefited by including quantitative tests of alternative diversification models that acknowledge extinction and diversity dependence.     

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.     

Climate-driven diversity dynamics in plants and plant-feeding insects

The origin of species-rich insect-plant food webs has traditionally been explained by diversifying antagonistic coevolution between plant defences and herbivore counter-defences. However, recent studies combining paleoclimatic reconstructions with time-calibrated phylogenies suggest that variation in global climate determines the distribution, abundance and diversity of plant clades and, hence, indirectly influences the balance between speciation and extinction in associated herbivore groups. Extant insect communities tend to be richest on common plant species that have many close relatives. This could be explained either by climate-driven diffuse cospeciation between plants and insects, or by elevated speciation and reduced extinction in herbivore lineages associated with expanding host taxa (resources). Progress in paleovegetation reconstructions in combination with the rapidly increasing availability of fossil-calibrated phylogenies provide means to discern between these alternative hypotheses. In particular, the 'Diffuse cospeciation' scenario predicts closely matching main diversification periods in plants and in the insects that feed upon them, while the 'Resource abundance-dependent diversification' hypothesis predicts that both positive and negative responses of insect diversity are lagged in relation to host-plant availability. The dramatic Cenozoic changes in global climate provide multiple possibilities for studying the mechanisms by which climatic shifts may drive diversity dynamics in plants and insect herbivores.     

Wood decomposition is more strongly controlled by temperature than by tree species and decomposer diversity in highly species rich subtropical forests

While the number of studies on the role of biodiversity on ecosystem functioning is steadily increasing, a key component of biogeochemical cycling in forests, dead wood decay, has been largely neglected. It remains widely unknown whether and how dead wood decay is affected by diversity loss in forests. We studied the hierarchical effects of tree species diversity on wood decay rates in a subtropical forest landscape in southeast China via its influence on fungal OTU richness and invertebrate diversity using piecewise structural equation models. The experiment was conducted in natural forest plots that span a wide gradient of tree species diversity embedded in a heterogeneous topography. To account for interactions between macro‐invertebrates and fungi, that potentially modify the influence of tree biodiversity and climate on dead wood decay, we compared a macro‐invertebrate exclusion treatment with a control treatment that allowed access to all types of decomposers. Diversity effects of trees on wood decay rates were mostly negative and mediated by the diversity of macro‐invertebrates. However, the effects of tree species diversity or fungal OTU richness and macro‐invertebrate diversity on wood decay rates were comparatively weak. Temperature affected decay rates positively and had the strongest influence in all treatments. While the exclusion of macro‐invertebrates did not lead to a reduction of wood decay rates, our results suggest that they may however have a mediating role in the process. In the presence of invertebrates the predictability of wood decay rates was higher and we observed a tendency of a stronger temperature control. Our results suggest that there is evidence for diversity effects on wood decomposition, but the temperature control is still more important. Thus, an increase in mean annual temperature will increase carbon and nutrient turnover through wood decomposition in subtropical forest irrespective of biotic composition.     

Environmental niche conservatism explains the accumulation of species richness in Mediterranean‐hotspot plant genera

The causes of exceptionally high plant diversity in Mediterranean‐climate biodiversity hotspots are not fully understood. We asked whether a mechanism similar to the tropical niche conservatism hypothesis could explain the diversity of four large genera (Protea, Moraea, Banksia, and Hakea) with distributions within and adjacent to the Greater Cape Floristic Region (South Africa) or the Southwest Floristic Region (Australia). Using phylogenetic and spatial data we estimated the environmental niche of each species, and reconstructed the mode and dynamics of niche evolution, and the geographic history, of each genus. For three genera, there were strong positive relationships between the diversity of clades within a region and their inferred length of occupation of that region. Within genera, there was evidence for strong evolutionary constraint on niche axes associated with climatic seasonality and aridity, with different niche optima for hotspot and nonhotspot clades. Evolutionary transitions away from hotspots were associated with increases in niche breadth and elevated rates of niche evolution. Our results point to a process of “hotspot niche conservatism” whereby the accumulation of plant diversity in Mediterranean‐type ecosystems results from longer time for speciation, with dispersal away from hotspots limited by narrow and phylogenetically conserved environmental niches.     

Phylogenies support out‐of‐equilibrium models of biodiversity

There is a long tradition in ecology of studying models of biodiversity at equilibrium. These models, including the influential Neutral Theory of Biodiversity, have been successful at predicting major macroecological patterns, such as species abundance distributions. But they have failed to predict macroevolutionary patterns, such as those captured in phylogenetic trees. Here, we develop a model of biodiversity in which all individuals have identical demographic rates, metacommunity size is allowed to vary stochastically according to population dynamics, and speciation arises naturally from the accumulation of point mutations. We show that this model generates phylogenies matching those observed in nature if the metacommunity is out of equilibrium. We develop a likelihood inference framework that allows fitting our model to empirical phylogenies, and apply this framework to various mammalian families. Our results corroborate the hypothesis that biodiversity dynamics are out of equilibrium.     

Ancient tropical extinctions at high latitudes contributed to the latitudinal diversity gradient*

Global biodiversity currently peaks at the equator and decreases toward the poles. Growing fossil evidence suggest this hump-shaped latitudinal diversity gradient (LDG) has not been persistent through time, with similar diversity across latitudes flattening out the LDG during past greenhouse periods. However, when and how diversity declined at high latitudes to generate the modern LDG remains an open question. Although diversity-loss scenarios have been proposed, they remain mostly undemonstrated. We outline the “asymmetric gradient of extinction and dispersal” framework that contextualizes previous ideas behind the LDG under a time-variable scenario. Using phylogenies and fossils of Testudines, Crocodilia, and Lepidosauria, we find that the hump-shaped LDG could be explained by (1) disproportionate extinctions of high-latitude tropical-adapted clades when climate transitioned from greenhouse to icehouse, and (2) equator-ward biotic dispersals tracking their climatic preferences when tropical biomes became restricted to the equator. Conversely, equivalent diversification rates across latitudes can account for the formation of an ancient flat LDG. The inclusion of fossils in macroevolutionary studies allows revealing time-dependent extinction rates hardly detectable from phylogenies only. This study underscores that the prevailing evolutionary processes generating the LDG during greenhouses differed from those operating during icehouses.     

Primary controls on species richness in higher taxa

The disparity in species richness across the tree of life is one of the most striking and pervasive features of biological diversity. Some groups are exceptionally diverse, whereas many other groups are species poor. Differences in diversity among groups are frequently assumed to result from primary control by differential rates of net diversification. However, a major alternative explanation is that ecological and other factors exert primary control on clade diversity, such that apparent variation in net diversification rates is a secondary consequence of ecological limits on clade growth. Here, I consider a likelihood framework for distinguishing between these competing hypotheses. I incorporate hierarchical modeling to explicitly relax assumptions about the constancy of diversification rates across clades, and I propose several statistics for a posteriori evaluation of model adequacy. I apply the framework to a recent dated phylogeny of ants. My results reject the hypothesis that net diversification rates exert primary control on species richness in this group and demonstrate that clade diversity is better explained by total time-integrated speciation. These results further suggest that it may not possible to estimate meaningful speciation and extinction rates from higher-level phylogenies of extant taxa only.     

Explaining bathymetric diversity patterns in marine benthic invertebrates and demersal fishes: physiological contributions to adaptation of life at depth

Bathymetric biodiversity patterns of marine benthic invertebrates and demersal fishes have been identified in the extant fauna of the deep continental margins. Depth zonation is widespread and evident through a transition between shelf and slope fauna from the shelf break to 1000 m, and a transition between slope and abyssal fauna from 2000 to 3000 m; these transitions are characterised by high species turnover. A unimodal pattern of diversity with depth peaks between 1000 and 3000 m, despite the relatively low area represented by these depths. Zonation is thought to result from the colonisation of the deep sea by shallow‐water organisms following multiple mass extinction events throughout the Phanerozoic. The effects of low temperature and high pressure act across hierarchical levels of biological organisation and appear sufficient to limit the distributions of such shallow‐water species. Hydrostatic pressures of bathyal depths have consistently been identified experimentally as the maximum tolerated by shallow‐water and upper bathyal benthic invertebrates at in situ temperatures, and adaptation appears required for passage to deeper water in both benthic invertebrates and demersal fishes. Together, this suggests that a hyperbaric and thermal physiological bottleneck at bathyal depths contributes to bathymetric zonation. The peak of the unimodal diversity–depth pattern typically occurs at these depths even though the area represented by these depths is relatively low. Although it is recognised that, over long evolutionary time scales, shallow‐water diversity patterns are driven by speciation, little consideration has been given to the potential implications for species distribution patterns with depth. Molecular and morphological evidence indicates that cool bathyal waters are the primary site of adaptive radiation in the deep sea, and we hypothesise that bathymetric variation in speciation rates could drive the unimodal diversity–depth pattern over time. Thermal effects on metabolic‐rate‐dependent mutation and on generation times have been proposed to drive differences in speciation rates, which result in modern latitudinal biodiversity patterns over time. Clearly, this thermal mechanism alone cannot explain bathymetric patterns since temperature generally decreases with depth. We hypothesise that demonstrated physiological effects of high hydrostatic pressure and low temperature at bathyal depths, acting on shallow‐water taxa invading the deep sea, may invoke a stress–evolution mechanism by increasing mutagenic activity in germ cells, by inactivating canalisation during embryonic or larval development, by releasing hidden variation or mutagenic activity, or by activating or releasing transposable elements in larvae or adults. In this scenario, increased variation at a physiological bottleneck at bathyal depths results in elevated speciation rate. Adaptation that increases tolerance to high hydrostatic pressure and low temperature allows colonisation of abyssal depths and reduces the stress–evolution response, consequently returning speciation of deeper taxa to the background rate. Over time this mechanism could contribute to the unimodal diversity–depth pattern.     

Phylogenetic evidence of a rapid radiation of pleurocarpous mosses (Bryophyta)

Abstract Pleurocarpous mosses, characterized by lateral female gametangia and highly branched, interwoven stems, comprise three orders and some 5000 species, or almost half of all moss diversity. Recent phylogenetic analyses resolve the Ptychomniales as sister to the Hypnales plus Hookeriales. Species richness is highly asymmetric with approximately 100 Ptychomniales, 750 Hookeriales, and 4400 Hypnales. Chloroplast DNA (cpDNA) sequences were obtained to compare partitioning of molecular diversity among the orders with estimates of species richness, and to test the hypothesis that either the Hookeriales or Hypnales underwent a period (or periods) of exceptionally rapid diversification. Levels of biodiversity were quantified using explicitly historical "phylogenetic diversity" and non-historical estimates of standing sequence diversity. Diversification rates were visualized using lineage-through-time (LTT) plots, and statistical tests of alternative diversification models were performed using the methods of Paradis (1997). The effects of incomplete sampling on the shape of LTT plots and performance of statistical tests were investigated using simulated phylogenies with incomplete sampling. Despite a much larger number of accepted species, the Hypnales contain lower levels of (cpDNA) biodiversity than their sister group, the Hookeriales, based on all molecular measures. Simulations confirm previous results that incomplete sampling yields diversification patterns that appear to reflect a decreasing rate through time, even when the true phylogenies were simulated with constant rates. Comparisons between simulated results and empirical data indicate that a constant rate of diversification cannot be rejected for the Hookeriales. The Hypnales, however, appear to have undergone a period of exceptionally rapid diversification for the earliest 20% of their history.