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.     

Estimating diversification rates for higher taxa: BAMM can give problematic estimates of rates and rate shifts

Estimates of diversification rates are invaluable for many macroevolutionary studies. Recently, an approach called BAMM (Bayesian Analysis of Macro‐evolutionary Mixtures) has become widely used for estimating diversification rates and rate shifts. At the same time, several articles have concluded that estimates of net diversification rates from the method‐of‐moments (MS) estimators are inaccurate. Yet, no studies have compared the ability of these two methods to accurately estimate clade diversification rates. Here, we use simulations to compare their performance. We found that BAMM yielded relatively weak relationships between true and estimated diversification rates. This occurred because BAMM underestimated the number of rates shifts across each tree, and assigned high rates to small clades with low rates. Errors in both speciation and extinction rates contributed to these errors, showing that using BAMM to estimate only speciation rates is also problematic. In contrast, the MS estimators (particularly using stem group ages), yielded stronger relationships between true and estimated diversification rates, by roughly twofold. Furthermore, the MS approach remained relatively accurate when diversification rates were heterogeneous within clades, despite the widespread assumption that it requires constant rates within clades. Overall, we caution that BAMM may be problematic for estimating diversification rates and rate shifts.     

In and out of refugia: historical patterns of diversity and demography in the North American Caesar's mushroom species complex

Some of the effects of past climate dynamics on plant and animal diversity make‐up have been relatively well studied, but to less extent in fungi. Pleistocene refugia are thought to harbour high biological diversity (i.e. phylogenetic lineages and genetic diversity), mainly as a product of increased reproductive isolation and allele conservation. In addition, high extinction rates and genetic erosion are expected in previously glaciated regions. Some of the consequences of past climate dynamics might involve changes in range and population size that can result in divergence and incipient or cryptic speciation. Many of these dynamic processes and patterns can be inferred through phylogenetic and coalescent methods. In this study, we first delimit species within a group of closely related edible ectomycorrhizal Amanita from North America (the American Caesar's mushrooms species complex) using multilocus coalescent‐based approaches; and then address questions related to effects of Pleistocene climate change on the diversity and genetics of the group. Our study includes extensive geographical sampling throughout the distribution range, and DNA sequences from three nuclear protein‐coding genes. Results reveal cryptic diversity and high speciation rates in refugia. Population sizes and expansions seem to be larger at midrange latitudes (Mexican highlands and SE USA). Range shifts are proportional to population size expansions, which were overall more common during the Pleistocene. This study documents responses to past climate change in fungi and also highlights the applicability of the multispecies coalescent in comparative phylogeographical analyses and diversity assessments that include ancestral species.     

Faster Speciation and Reduced Extinction in the Tropics Contribute to the Mammalian Latitudinal Diversity Gradient

The increase in species richness from the poles to the tropics, referred to as the latitudinal diversity gradient, is one of the most ubiquitous biodiversity patterns in the natural world. Although understanding how rates of speciation and extinction vary with latitude is central to explaining this pattern, such analyses have been impeded by the difficulty of estimating diversification rates associated with specific geographic locations. Here, we use a powerful phylogenetic approach and a nearly complete phylogeny of mammals to estimate speciation, extinction, and dispersal rates associated with the tropical and temperate biomes. Overall, speciation rates are higher, and extinction rates lower, in the tropics than in temperate regions. The diversity of the eight most species-rich mammalian orders (covering 92% of all mammals) peaks in the tropics, except that of the Lagomorpha (hares, rabbits, and pikas) reaching a maxima in northern-temperate regions. Latitudinal patterns in diversification rates are strikingly consistent with these diversity patterns, with peaks in species richness associated with low extinction rates (Primates and Lagomorpha), high speciation rates (Diprotodontia, Artiodactyla, and Soricomorpha), or both (Chiroptera and Rodentia). Rates of range expansion were typically higher from the tropics to the temperate regions than in the other direction, supporting the “out of the tropics” hypothesis whereby species originate in the tropics and disperse into higher latitudes. Overall, these results suggest that differences in diversification rates have played a major role in shaping the modern latitudinal diversity gradient in mammals, and illustrate the usefulness of recently developed phylogenetic approaches for understanding this famous yet mysterious pattern.     

Quantitative traits and diversification

Quantitative traits have long been hypothesized to affect speciation and extinction rates. For example, smaller body size or increased specialization may be associated with increased rates of diversification. Here, I present a phylogenetic likelihood-based method (quantitative state speciation and extinction [QuaSSE]) that can be used to test such hypotheses using extant character distributions. This approach assumes that diversification follows a birth-death process where speciation and extinction rates may vary with one or more traits that evolve under a diffusion model. Speciation and extinction rates may be arbitrary functions of the character state, allowing much flexibility in testing models of trait-dependent diversification. I test the approach using simulated phylogenies and show that a known relationship between speciation and a quantitative character could be recovered in up to 80% of the cases on large trees (500 species). Consistent with other approaches, detecting shifts in diversification due to differences in extinction rates was harder than when due to differences in speciation rates. Finally, I demonstrate the application of QuaSSE to investigate the correlation between body size and diversification in primates, concluding that clade-specific differences in diversification may be more important than size-dependent diversification in shaping the patterns of diversity within this group.     

Minimal effects of latitude on present-day speciation rates in New World birds

The tropics contain far greater numbers of species than temperate regions, suggesting that rates of species formation might differ systematically between tropical and non-tropical areas. We tested this hypothesis by reconstructing the history of speciation in New World (NW) land birds using BAMM, a Bayesian framework for modelling complex evolutionary dynamics on phylogenetic trees. We estimated marginal distributions of present-day speciation rates for each of 2571 species of birds. The present-day rate of speciation varies approximately 30-fold across NW birds, but there is no difference in the rate distributions for tropical and temperate taxa. Using macroevolutionary cohort analysis, we demonstrate that clades with high tropical membership do not produce species more rapidly than temperate clades. For nearly any value of present-day speciation rate, there are far more species in the tropics than the temperate zone. Any effects of latitude on speciation rate are marginal in comparison to the dramatic variation in rates among clades.     

THE ROLE OF CLIMATIC TOLERANCES AND SEED TRAITS IN REDUCED EXTINCTION RATES OF TEMPERATE POLYGONACEAE

The latitudinal diversity gradient (LDG) is one of the most striking and consistent biodiversity patterns across taxonomic groups. We investigate the species richness gradient in the buckwheat family, Polygonaceae, which exhibits a reverse LDG and is, thus, decoupled from dominant gradients of energy and environmental stability that increase toward the tropics and confound mechanistic interpretations. We test competing age and evolutionary diversification hypotheses, which may explain the diversification of this plant family over the past 70 million years. Our analyses show that the age hypothesis, which posits that clade richness is positively correlated with the ecological and evolutionary time since clade origin, fails to explain the richness gradient observed in Polygonaceae. However, an evolutionary diversification hypothesis is highly supported, with diversification rates being 3.5 times higher in temperate clades compared to tropical clades. We demonstrate that differences in rates of speciation, migration, and molecular evolution insufficiently explain the observed patterns of differential diversification rates. We suggest that reduced extinction rates in temperate clades may be associated with adaptive responses to selection, through which seed morphology and climatic tolerances potentially act to minimize risk in temporally variable environments. Further study is needed to understand causal pathways among these traits and factors correlated with latitude.     

Origins of marine patterns of biodiversity: some correlates and applications

Abstract: Marine shelf diversity patterns correlate with macroecological features of basic importance that may play causal roles in macroevolution. We have investigated the global diversity pattern of living Bivalvia, which is dominated by the latitudinal diversity gradient (LDG), maintained by high tropical origination rates. Generic‐level lineages expand poleward, chiefly through speciation, so that species richness within provinces and globally is positively correlated with generic geographical ranges. A gradient in diversity accommodation progressively lowers both immigration and speciation rates in higher latitudes. The LDG correlates with seasonality of trophic resources but not with area; tropical provinces are not diverse because they are large but because they are tropical. A similar dynamic evidently underlays Jurassic and Carboniferous LDGs. Larval developmental modes correlate with the LDG and thus with resource seasonality, with tropical dominance of planktotrophs offset by increasing nonplanktotrophy to poleward. The acquisition of planktotrophy in several early Palaeozoic clades indicates a change in macroecological relationships during Cambrian and Ordovician radiations.     

EVOLUTIONARY AND BIOGEOGRAPHIC ORIGINS OF HIGH TROPICAL DIVERSITY IN OLD WORLD FROGS (RANIDAE)

Differences in species richness between regions are ultimately explained by patterns of speciation, extinction, and biogeographic dispersal. Yet, few studies have considered the role of all three processes in generating the high biodiversity of tropical regions. A recent study of a speciose group of predominately New World frogs (Hylidae) showed that their low diversity in temperate regions was associated with relatively recent colonization of these regions, rather than latitudinal differences in diversification rates (rates of speciation–extinction). Here, we perform parallel analyses on the most species-rich group of Old World frogs (Ranidae; ∼1300 species) to determine if similar processes drive the latitudinal diversity gradient. We estimate a time-calibrated phylogeny for 390 ranid species and use this phylogeny to analyze patterns of biogeography and diversification rates. As in hylids, we find a strong relationship between the timing of colonization of each region and its current diversity, with recent colonization of temperate regions from tropical regions. Diversification rates are similar in tropical and temperate clades, suggesting that neither accelerated tropical speciation rates nor greater temperate extinction rates explain high tropical diversity in this group. Instead, these results show the importance of historical biogeography in explaining high species richness in both the New World and Old World tropics.     

FiSSE: A simple nonparametric test for the effects of a binary character on lineage diversification rates

It is widely assumed that phenotypic traits can influence rates of speciation and extinction, and several statistical approaches have been used to test for correlations between character states and lineage diversification. Recent work suggests that model‐based tests of state‐dependent speciation and extinction are sensitive to model inadequacy and phylogenetic pseudoreplication. We describe a simple nonparametric statistical test (“FiSSE”) to assess the effects of a binary character on lineage diversification rates. The method involves computing a test statistic that compares the distributions of branch lengths for lineages with and without a character state of interest. The value of the test statistic is compared to a null distribution generated by simulating character histories on the observed phylogeny. Our tests show that FiSSE can reliably infer trait‐dependent speciation on phylogenies of several hundred tips. The method has low power to detect trait‐dependent extinction but can infer state‐dependent differences in speciation even when net diversification rates are constant. We assemble a range of macroevolutionary scenarios that are problematic for likelihood‐based methods, and we find that FiSSE does not show similarly elevated false positive rates. We suggest that nonparametric statistical approaches, such as FiSSE, provide an important complement to formal process‐based models for trait‐dependent diversification.     

Evolution and the latitudinal diversity gradient: speciation, extinction and biogeography

A latitudinal gradient in biodiversity has existed since before the time of the dinosaurs, yet how and why this gradient arose remains unresolved. Here we review two major hypotheses for the origin of the latitudinal diversity gradient. The time and area hypothesis holds that tropical climates are older and historically larger, allowing more opportunity for diversification. This hypothesis is supported by observations that temperate taxa are often younger than, and nested within, tropical taxa, and that diversity is positively correlated with the age and area of geographical regions. The diversification rate hypothesis holds that tropical regions diversify faster due to higher rates of speciation (caused by increased opportunities for the evolution of reproductive isolation, or faster molecular evolution, or the increased importance of biotic interactions), or due to lower extinction rates. There is phylogenetic evidence for higher rates of diversification in tropical clades, and palaeontological data demonstrate higher rates of origination for tropical taxa, but mixed evidence for latitudinal differences in extinction rates. Studies of latitudinal variation in incipient speciation also suggest faster speciation in the tropics. Distinguishing the roles of history, speciation and extinction in the origin of the latitudinal gradient represents a major challenge to future research.     

The great American biotic interchange and diversification history in Dynastes beetles (Scarabaeidae; Dynastinae)

Biotic interchange between geographic regions can promote rapid diversification. However, what are the important factors that determine the rate of diversification (e.g., trait‐dependent diversification) vary between study systems. The evolutionary history of Dynastes beetles, which can be found in both North and South Americas and exhibit two different altitudinal preferences (highland and lowland) is tested for the effects of biotic interchange between continents and different ecological preferences on the rate of species diversification. Additionally, the hypotheses of geological time‐dependent and lineage specific diversification rates are also tested. Results from this study indicate that in Dynastes beetles a pre‐landbridge dispersal hypothesis from South to North America is preferred and that the speciation rates estimated using BAMM are similar between lineages of different geographic origins and different altitudinal preferences (i.e., diversification rate is not trait‐dependent). On the other hand, my result from marcoevolutionary cohort analysis based on BAMM outputs suggests that the rate of speciation in Dynastes beetles is, instead of trait‐dependent (geographic and ecological), lineage specific. Furthermore, a steadily increasing speciation rate can be found in Pliocene and Pleistocene, which implies that geological and climatic events, i.e., colonizing North America, habitat reformation in the Amazonia, and forest contraction in Pleistocene, may have together shaped the current biodiversity pattern in Dynastes beetles.     

Large-scale phylogenetic analyses reveal the causes of high tropical amphibian diversity

Many groups show higher species richness in tropical regions but the underlying causes remain unclear. Despite many competing hypotheses to explain latitudinal diversity gradients, only three processes can directly change species richness across regions: speciation, extinction and dispersal. These processes can be addressed most powerfully using large-scale phylogenetic approaches, but most previous studies have focused on small groups and recent time scales, or did not separate speciation and extinction rates. We investigate the origins of high tropical diversity in amphibians, applying new phylogenetic comparative methods to a tree of 2871 species. Our results show that high tropical diversity is explained by higher speciation in the tropics, higher extinction in temperate regions and limited dispersal out of the tropics compared with colonization of the tropics from temperate regions. These patterns are strongly associated with climate-related variables such as temperature, precipitation and ecosystem energy. Results from models of diversity dependence in speciation rate suggest that temperate clades may have lower carrying capacities and may be more saturated (closer to carrying capacity) than tropical clades. Furthermore, we estimate strikingly low tropical extinction rates over geological time scales, in stark contrast to the dramatic losses of diversity occurring in tropical regions presently.     

Can the tropical conservatism hypothesis explain temperate species richness patterns? An inverse latitudinal biodiversity gradient in the New World snake tribe Lampropeltini

Aim  A latitudinal gradient in species richness, defined as a decrease in biodiversity away from the equator, is one of the oldest known patterns in ecology and evolutionary biology. However, there are also many known cases of increasing poleward diversity, forming inverse latitudinal biodiversity gradients. As only three processes (speciation, extinction and dispersal) can directly affect species richness in areas, similar factors may be responsible for both classical (high tropical diversity) and inverse (high temperate diversity) gradients. Thus, a modified explanation for differential species richness which accounts for both patterns would be preferable to one which only explains high tropical biodiversity.
Location  The New World.
Methods  We test several proposed ecological, temporal, evolutionary and spatial explanations for latitudinal diversity gradients in the New World snake tribe Lampropeltini, which exhibits its highest biodiversity in temperate regions.
Results  We find that an extratropical peak in species richness is not explained by latitudinal variation in diversification rate, the mid-domain effect, or Rapoport's rule. Rather, earlier colonization and longer duration in the temperate zones allowing more time for speciation to increase biodiversity, phylogenetic niche conservatism limiting tropical dispersal and the expansion of the temperate zones in the Tertiary better explain inverse diversity gradients in this group.
Main conclusions  Our conclusions are the inverse of the predictions made by the tropical conservatism hypothesis to explain higher biodiversity near the equator. Therefore, we suggest that the processes invoked are not intrinsic to the tropics but are dependent on historical biogeography to determine the distribution of species richness, which we refer to as the 'biogeographical conservatism hypothesis'.     

WHAT CAN MULTIPLE PHYLOGENIES SAY ABOUT THE LATITUDINAL DIVERSITY GRADIENT? A NEW LOOK AT THE TROPICAL CONSERVATISM,OUT OF THE TROPICS,AND DIVERSIFICATION RATE HYPOTHESES

We reviewed published phylogenies and selected 111 phylogenetic studies representing mammals, birds, insects, and flowering plants. We then mapped the latitudinal range of all taxa to test the relative importance of the tropical conservatism, out of the tropics, and diversification rate hypotheses in generating latitudinal diversity gradients. Most clades originated in the tropics, with diversity peaking in the zone of origin. Transitions of lineages between latitudinal zones occurred at 16–22% of the tree nodes. The most common type of transition was range expansions of tropical lineages to encompass also temperate latitudes. Thus, adaptation to new climatic conditions may not represent a major obstacle for many clades. These results contradict predictions of the tropical conservatism hypothesis (i.e., few clades colonizing extratropical latitudes), but support the out‐of‐the‐tropics model (i.e., tropical originations and subsequent latitudinal range expansions). Our results suggest no difference in diversification between tropical and temperate sister lineages; thus, diversity of tropical clades was not explained by higher diversification rates in this zone. Moreover, lineages with latitudinal stasis diversified more compared to sister lineages entering a new latitudinal zone. This preserved preexisting diversity differences between latitudinal zones and can be considered a new mechanism for why diversity tends to peak in the zone of origin.     

Evolutionary and ecological causes of the latitudinal diversity gradient in hylid frogs: treefrog trees unearth the roots of high tropical diversity

Why are there more species in the tropics than in temperate regions? In recent years, this long-standing question has been addressed primarily by seeking environmental correlates of diversity. But to understand the ultimate causes of diversity patterns, we must also examine the evolutionary and biogeographic processes that directly change species numbers (i.e., speciation, extinction, and dispersal). With this perspective, we dissect the latitudinal diversity gradient in hylid frogs. We reconstruct a phylogeny for 124 hylid species, estimate divergence times and diversification rates for major clades, reconstruct biogeographic changes, and use ecological niche modeling to identify climatic variables that potentially limit dispersal. We find that hylids originated in tropical South America and spread to temperate regions only recently (leaving limited time for speciation). There is a strong relationship between the species richness of each region and when that region was colonized but not between the latitudinal positions of clades and their rates of diversification. Temperature seasonality seemingly limits dispersal of many tropical clades into temperate regions and shows significant phylogenetic conservatism. Overall, our study illustrates how two general principles (niche conservatism and the time-for-speciation effect) may help explain the latitudinal diversity gradient as well as many other diversity patterns across taxa and regions.     

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.     

Into and out of the tropics: global diversification patterns in a hyperdiverse clade of ectomycorrhizal fungi

Ectomycorrhizal (ECM) fungi, symbiotic mutualists of many dominant tree and shrub species, exhibit a biogeographic pattern counter to the established latitudinal diversity gradient of most macroflora and fauna. However, an evolutionary basis for this pattern has not been explicitly tested in a diverse lineage. In this study, we reconstructed a mega‐phylogeny of a cosmopolitan and hyperdiverse genus of ECM fungi, Russula, sampling from annotated collections and utilizing publically available sequences deposited in GenBank. Metadata from molecular operational taxonomic unit cluster sets were examined to infer the distribution and plant association of the genus. This allowed us to test for differences in patterns of diversification between tropical and extratropical taxa, as well as how their associations with different plant lineages may be a driver of diversification. Results show that Russula is most species‐rich at temperate latitudes and ancestral state reconstruction shows that the genus initially diversified in temperate areas. Migration into and out of the tropics characterizes the early evolution of the genus, and these transitions have been frequent since this time. We propose the ‘generalized diversification rate’ hypothesis to explain the reversed latitudinal diversity gradient pattern in Russula as we detect a higher net diversification rate in extratropical lineages. Patterns of diversification with plant associates support host switching and host expansion as driving diversification, with a higher diversification rate in lineages associated with Pinaceae and frequent transitions to association with angiosperms.     

Global patterns of diversification and species richness in amphibians

Geographic patterns of species richness ultimately arise through the processes of speciation, extinction, and dispersal, but relatively few studies consider evolutionary and biogeographic processes in explaining these diversity patterns. One explanation for high tropical species richness is that many species-rich clades originated in tropical regions and spread to temperate regions infrequently and more recently, leaving little time for species richness to accumulate there (assuming similar rates of diversification in temperate and tropical regions). However, the major clades of anurans (frogs) and salamanders may offer a compelling counterexample. Most salamander families are predominately temperate in distribution, but the one primarily tropical clade (Bolitoglossinae) contains nearly half of all salamander species. Similarly, most basal clades of anurans are predominately temperate, but one largely tropical clade (Neobatrachia) contains approximately 96% of anurans. In this article, I examine patterns of diversification in frogs and salamanders and their relationship to large-scale patterns of species richness in amphibians. I find that diversification rates in both frogs and salamanders increase significantly with decreasing latitude. These results may shed light on both the evolutionary causes of the latitudinal diversity gradient and the dramatic but poorly explained disparities in the diversity of living amphibian clades.     

Detecting Trait-Dependent Diversification Under Diversification Slowdowns

Testing whether a certain biological trait significantly affects clade diversification is central to macroevolutionary research. To this end, many scientists use constant-rate estimators (CR estimators) of diversification. However, it has never been examined whether these estimators report meaningful relationships between traits and diversification even when the diversification itself decelerates over time. In this study, I simulate trait-driven diversification concurrently with diversification slowdowns. Then, I test whether CR estimators manage to uncover the simulated relationships. Results suggest that CR estimators are robust against violation of rate constancy and successfully detect trait-dependent diversification in spite of diversification declines. Interestingly, correct results were recovered whether clade age correlated with clade diversity or not. Further comparison of CR estimators with QuaSSE suggested that QuaSSE performs better under constant diversification, but tends to report spuriously significant outcomes when diversification decelerates (=elevated Type I error). Given that diversification slowdowns have been recently reported for a wide range of taxa, these findings may be of particular relevance for future diversification studies.