The causes of species richness patterns across space, time, and clades and the role of "ecological limits"

A major goal of research in ecology and evolution is to explain why species richness varies across habitats, regions, and clades. Recent reviews have argued that species richness patterns among regions and clades may be explained by "ecological limits" on diversity over time, which are said to offer an alternative explanation to those invoking speciation and extinction (diversification) and time. Further, it has been proposed that this hypothesis is best supported by failure to find a positive relationship between time (e.g., clade age) and species richness. Here, I critically review the evidence for these claims, and propose how we might better study the ecological and evolutionary origins of species richness patterns. In fact, ecological limits can only influence species richness in clades by influencing speciation and extinction, and so this new "alternative paradigm" is simply one facet of the traditional idea that ecology influences diversification. The only direct evidence for strict ecological limits on richness (i.e., constant diversity over time) is from the fossil record, but many studies cited as supporting this pattern do not, and there is evidence for increasing richness over time. Negative evidence for a relationship between clade age and richness among extant clades is not positive evidence for constant diversity over time, and many recent analyses finding no age-diversity relationship were biased to reach this conclusion. More comprehensive analyses strongly support a positive age-richness relationship. There is abundant evidence that both time and ecological influences on diversification rates are important drivers of both large-scale and small-scale species richness patterns. The major challenge for future studies is to understand the ecological and evolutionary mechanisms underpinning the relationships between time, dispersal, diversification, and species richness patterns.     

When should species richness be energy limited,and how would we know?

Energetic constraints are fundamental to ecology and evolution, and empirical relationships between species richness and estimates of available energy (i.e. resources) have led some to suggest that richness is energetically constrained. However, the mechanism linking energy with richness is rarely specified and predictions of secondary patterns consistent with energy‐constrained richness are lacking. Here, we lay out the necessary and sufficient assumptions of a causal relationship linking energy gradients to richness gradients. We then describe an eco‐evolutionary simulation model that combines spatially explicit diversification with trait evolution, resource availability and assemblage‐level carrying capacities. Our model identified patterns in richness and phylogenetic structure expected when a spatial gradient in energy availability determines the number of individuals supported in a given area. A comparison to patterns under alternative scenarios, in which fundamental assumptions behind energetic explanations were violated, revealed patterns that are useful for evaluating the importance of energetic constraints in empirical systems. We use a data set on rockfish (genus Sebastes) from the northeastern Pacific to show how empirical data can be coupled with model predictions to evaluate the role of energetic constraints in generating observed richness gradients.     

Historical processes enhance patterns of diversity along latitudinal gradients

One of the more vexing issues in ecology is how historical processes affect contemporary patterns of biodiversity. Accordingly, few models have been presented. Two corollary models (centre of origin, time-for-speciation) can be used to make quantitative predictions characterizing the tropical niche conservatism hypothesis and describe diversification as diffusion and subsequent cladogenesis of species away from the place of origin of a higher taxon in the tropics. Predictions derived from such models are: (i) species richness declines toward the periphery of the range of a higher taxon; (ii) taxa are more derived toward the periphery than the centre; (iii) ages of taxa are lower toward the periphery than the centre; and (iv) ages and measures of derivedness are less variable toward the periphery of the range of a higher taxon. I tested these predictions to better understand the formation of one of the most ubiquitous patterns of biodiversity-the latitudinal gradient in species richness. Results indicate well-supported predictions for New World leaf-nosed bats and that diversification has had strong influences on latitudinal gradients of species richness. A better understanding of how evolutionary diversification of taxa contributes to formation of patterns of species richness along environmental gradients is necessary to fully understand spatial variation in biodiversity.     

Niche conservatism as an emerging principle in ecology and conservation biology

The diversity of life is ultimately generated by evolution, and much attention has focused on the rapid evolution of ecological traits. Yet, the tendency for many ecological traits to instead remain similar over time [niche conservatism (NC)] has many consequences for the fundamental patterns and processes studied in ecology and conservation biology. Here, we describe the mounting evidence for the importance of NC to major topics in ecology (e.g. species richness, ecosystem function) and conservation (e.g. climate change, invasive species). We also review other areas where it may be important but has generally been overlooked, in both ecology (e.g. food webs, disease ecology, mutualistic interactions) and conservation (e.g. habitat modification). We summarize methods for testing for NC, and suggest that a commonly used and advocated method (involving a test for phylogenetic signal) is potentially problematic, and describe alternative approaches. We suggest that considering NC: (1) focuses attention on the within‐species processes that cause traits to be conserved over time, (2) emphasizes connections between questions and research areas that are not obviously related (e.g. invasives, global warming, tropical richness), and (3) suggests new areas for research (e.g. why are some clades largely nocturnal? why do related species share diseases?).     

What explains patterns of species richness? The relative importance of climatic‐niche evolution,morphological evolution,and ecological limits in salamanders

A major goal of evolutionary biology and ecology is to understand why species richness varies among clades. Previous studies have suggested that variation in richness among clades might be related to variation in rates of morphological evolution among clades (e.g., body size and shape). Other studies have suggested that richness patterns might be related to variation in rates of climatic‐niche evolution. However, few studies, if any, have tested the relative importance of these variables in explaining patterns of richness among clades. Here, we test their relative importance among major clades of Plethodontidae, the most species‐rich family of salamanders. Earlier studies have suggested that climatic‐niche evolution explains patterns of diversification among plethodontid clades, whereas rates of morphological evolution do not. A subsequent study stated that rates of morphological evolution instead explained patterns of species richness among plethodontid clades (along with “ecological limits” on richness of clades, leading to saturation of clades with species, given limited resources). However, they did not consider climatic‐niche evolution. Using phylogenetic multiple regression, we show that rates of climatic‐niche evolution explain most variation in richness among plethodontid clades, whereas rates of morphological evolution do not. We find little evidence that ecological limits explain patterns of richness among plethodontid clades. We also test whether rates of morphological and climatic‐niche evolution are correlated, and find that they are not. Overall, our results help explain richness patterns in a major amphibian group and provide possibly the first test of the relative importance of climatic niches and morphological evolution in explaining diversity patterns.     

Spatial diversity patterns of Pristimantis frogs in the Tropical Andes

Although biodiversity gradients have been widely documented, the factors governing broad‐scale patterns in species richness are still a source of intense debate and interest in ecology, evolution, and conservation biology. Here, we tested whether spatial hypotheses (species–area effect, topographic heterogeneity, mid‐domain null model, and latitudinal effect) explain the pattern of diversity observed along the altitudinal gradient of Andean rain frogs of the genus Pristimantis. We compiled a gamma‐diversity database of 378 species of Pristimantis from the tropical Andes, specifically from Colombia to Bolivia, using records collected above 500 m.a.s.l. Analyses were performed at three spatial levels: Tropical Andes as a whole, split in its two main domains (Northern and Central Andes), and split in its 11 main mountain ranges. Species richness, area, and topographic heterogeneity were calculated for each 500‐m‐width elevational band. Spatial hypotheses were tested using linear regression models. We examined the fit of the observed diversity to the mid‐domain hypothesis using randomizations. The species richness of Pristimantis showed a hump‐shaped pattern across most of the altitudinal gradients of the Tropical Andes. There was high variability in the relationship between area and species richness along the Tropical Andes. Correcting for area effects had little impact in the shape of the empirical pattern of biodiversity curves. Mid‐domain models produced similar gradients in species richness relative to empirical gradients, but the fit varied among mountain ranges. The effect of topographic heterogeneity on species richness varied among mountain ranges. There was a significant negative relationship between latitude and species richness. Our findings suggest that spatial processes partially explain the richness patterns of Pristimantis frogs along the Tropical Andes. Explaining the current patterns of biodiversity in this hot spot may require further studies on other possible underlying mechanisms (e.g., historical, biotic, or climatic hypotheses) to elucidate the factors that limit the ranges of species along this elevational gradient.     

Richness gradients of stream invertebrates across the USA: taxonomy- and trait-based approaches

Large-scale diversity patterns in relationship to environmental factors at multiple spatial scales have been well-studied for many taxonomic groups; however, freshwater ecosystems remain understudied. Biodiversity is now widely recognized to encompass many more factors than just species numbers, particularly the inclusion of functional attributes. In this study, we examined richness patterns of stream invertebrate genera and their biological traits (“functional” richness) across 364 sites in the contiguous USA. In particular, we focused on the relationship between taxonomy- and trait-based richness to test for functional redundancy in stream communities. Further, we obtained environmental data to model the relative importance of local and watershed-scale environmental factors and residual spatial (latitude, longitude) influences on taxonomy- and trait-based richness. Trait richness increased linearly with genus richness (slope ≪ 1), although this appears to be an artifact of the restricted range of genus richness in our study (32 genera maximum). Furthermore, trait richness was significantly lower than expected under random community assembly. In contrast, the Ephemeroptera, Plecoptera, and Trichoptera (EPT) genera exhibited a saturating pattern between trait and genus richness and trait richness was no different from random. Our study indicates that there is functional redundancy among stream invertebrate genera, likely as a result of harsh habitat filters limiting trait diversity. Environmental factors (including spatially structured environmental factors) were always more important than spatial factors (latitude, longitude) in structuring richness despite strong longitudinal patterns of all richness measures (these differences were only significant for EPT genera). Finally, we found no significant difference in the relative importance of local and watershed scale environmental factors for taxonomy- and trait-based richness.     

Explaining large-scale patterns of vertebrate diversity

The major clades of vertebrates differ dramatically in their current species richness, from 2 to more than 32 000 species each, but the causes of this variation remain poorly understood. For example, a previous study noted that vertebrate clades differ in their diversification rates, but did not explain why they differ. Using a time-calibrated phylogeny and phylogenetic comparative methods, I show that most variation in diversification rates among 12 major vertebrate clades has a simple ecological explanation: predominantly terrestrial clades (i.e. birds, mammals, and lizards and snakes) have higher net diversification rates than predominantly aquatic clades (i.e. amphibians, crocodilians, turtles and all fish clades). These differences in diversification rates are then strongly related to patterns of species richness. Habitat may be more important than other potential explanations for richness patterns in vertebrates (such as climate and metabolic rates) and may also help explain patterns of species richness in many other groups of organisms.     

Intraspecific phenotypic variation in a fish predator affects multitrophic lake metacommunity structure

Contemporary insights from evolutionary ecology suggest that population divergence in ecologically important traits within predators can generate diversifying ecological selection on local community structure. Many studies acknowledging these effects of intraspecific variation assume that local populations are situated in communities that are unconnected to similar communities within a shared region. Recent work from metacommunity ecology suggests that species dispersal among communities can also influence species diversity and composition but can depend upon the relative importance of the local environment. Here, we study the relative effects of intraspecific phenotypic variation in a fish predator and spatial processes related to plankton species dispersal on multitrophic lake plankton metacommunity structure. Intraspecific diversification in foraging traits and residence time of the planktivorous fish alewife (Alosa pseudoharengus) among coastal lakes yields lake metacommunities supporting three lake types which differ in the phenotype and incidence of alewife: lakes with anadromous, landlocked, or no alewives. In coastal lakes, plankton community composition was attributed to dispersal versus local environmental predictors, including intraspecific variation in alewives. Local and beta diversity of zooplankton and phytoplankton was additionally measured in response to intraspecific variation in alewives. Zooplankton communities were structured by species sorting, with a strong influence of intraspecific variation in A. pseudoharengus. Intraspecific variation altered zooplankton species richness and beta diversity, where lake communities with landlocked alewives exhibited intermediate richness between lakes with anadromous alewives and without alewives, and greater community similarity. Phytoplankton diversity, in contrast, was highest in lakes with landlocked alewives. The results indicate that plankton dispersal in the region supplied a migrant pool that was strongly structured by intraspecific variation in alewives. This is one of the first studies to demonstrate that intraspecific phenotypic variation in a predator can maintain contrasting patterns of multitrophic diversity in metacommunities.     

Regional variation in the historical components of global avian species richness

Aim Using a global data base of the distribution of extant bird species, we examine the evidence for spatial variation in the evolutionary origins of contemporary avian diversity. In particular, we assess the possible role of the timing of mountain uplift in promoting diversification in different regions. Location Global. Methods We mapped the distribution of avian richness at four taxonomic levels on an equal‐area 1° grid. We examined the relationships between richness at successive taxonomic levels (e.g. species richness vs. genus richness). We mapped the residuals from linear regressions of these relationships to identify areas that are exceptional in the number of lower taxa relative to the number of higher taxa. We use generalized least squares models to test the influence of elevation range and temperature on lower‐taxon richness relative to higher‐taxon richness. Results Peaks of species richness in the Neotropics were congruent with patterns of generic richness, whilst peaks in Australia and the Himalayas were congruent with patterns of both genus and family richness. Hotspots in the Afrotropics did not reflect higher‐taxon patterns. Regional differences in the relationship between richness at successive taxonomic levels revealed variation in patterns of taxon co‐occurrence. Species and genus co‐occurrence was positively associated with elevational range across much of the world. Taxon occurrence in the Neotropics was associated with a positive interaction between elevational range and temperature. Conclusions These results demonstrate that contemporary patterns of richness show different associations with higher‐taxon richness in different regions, which implies that the timing of historical effects on these contemporary patterns varies across regions. We suggest that this is due to dispersal limitation and phylogenetic constraints on physiological tolerance limits promoting diversification. We speculate that diversification rates respond to long‐term changes in the Earth's topography, and that the role of tropical mountain ranges is implicated as a correlate of contemporary diversity, and a source of diversification across avian evolutionary history.     

Environmental correlates of leguminosae species richness in Mexico: Quantifying the contributions of energy and environmental seasonality

Explaining species richness patterns is a central issue in ecology, but a general explanation remains elusive. Environmental conditions have been proposed to be important drivers of these patterns, but we still need to better understand the relative contribution of environmental factors. Here, we aim at testing two environmental hypotheses for richness gradients: energy availability and environmental seasonality using diversity patterns of the family Leguminosae across Mexico. We compiled a data base of 502 species and 32,962 records. After dividing Mexico into 100 × 100 km grid cells, we constructed a map of variation in species richness that accounts for heterogeneity in sampling effort. We found the cells with the highest species richness of legumes are in the Neotropical region of Pacific coastal and southern Mexico, where the legume family dominates the tropical rain forests and seasonally dry tropical forests. Regression models show that energy and seasonality predictors can explain 25% and 49% of the variation in richness, respectively. Spatial autocorrelation analysis showed that richness has a strong spatial structure, but that most of this structure disappears when both energy and seasonality are used to account for richness gradient. Our study demonstrates multiple environmental conditions contribute complementarily to explain diversity gradients. Moreover, it shows that in some regions, environmental seasonality can be more important than energy availability, contradicting studies at coarser spatial scales. More basic taxonomic and floristic work is needed to help describe patterns of diversity for many groups to allow for testing the underlying mechanisms responsible for diversity gradients. Abstract in Spanish is available with online material.     

Clade age and not diversification rate explains species richness among animal taxa

Animal taxa show remarkable variability in species richness across phylogenetic groups. Most explanations for this disparity postulate that taxa with more species have phenotypes or ecologies that cause higher diversification rates (i.e., higher speciation rates or lower extinction rates). Here we show that clade longevity, and not diversification rate, has primarily shaped patterns of species richness across major animal clades: more diverse taxa are older and thus have had more time to accumulate species. Diversification rates calculated from 163 species-level molecular phylogenies were highly consistent within and among three major animal phyla (Arthropoda, Chordata, Mollusca) and did not correlate with species richness. Clades with higher estimated diversification rates were younger, but species numbers increased with increasing clade age. A fossil-based data set also revealed a strong, positive relationship between total extant species richness and crown group age across the orders of insects and vertebrates. These findings do not negate the importance of ecology or phenotype in influencing diversification rates, but they do show that clade longevity is the dominant signal in major animal biodiversity patterns. Thus, some key innovations may have acted through fostering clade longevity and not by heightening diversification rate.     

Faster diversification on land than sea helps explain global biodiversity patterns among habitats and animal phyla

Terrestrial environments occupy ~ 30% of the Earth's surface yet contain ~ 80% of all species. The causes of this dramatic biodiversity gradient have remained relatively unstudied. Here, I test the fundamental prediction that predominantly non‐marine clades have more rapid rates of diversification than marine clades, using a time‐calibrated phylogeny of animal phyla. The results strongly support this hypothesis. This pattern helps explain the higher richness of terrestrial environments and the dramatic variation in species richness among animal phyla. The results show the importance of ecology in explaining large‐scale patterns of clade richness and of diversification rates in explaining Earth's largest biodiversity patterns. The results also demonstrate remarkable niche conservatism in habitats, in some cases lasting > 800 million years. Finally, the results highlight the surprisingly high species richness of freshwater habitats, which are nearly equal to marine environments despite their much smaller area (~ 2% of Earth's surface vs. 70% for marine habitats).     

Spatial Scaling of Non-Native Fish Richness across the United States

A major goal and challenge of invasion ecology is to describe and interpret spatial and temporal patterns of species invasions. Here, we examined fish invasion patterns at four spatially structured and hierarchically nested scales across the contiguous United States (i.e., from large to small: region, basin, watershed, and sub-watershed). All spatial relationships in both richness and fraction between species groups (e.g., natives vs. exotics) were positive at large scales. However, contrary to predictions using null/neutral models, the patterns at small scales were hump-shaped (unimodal), not simply negative. The fractions of both domestic (introduced among watersheds within the USA) and foreign (introduced from abroad) exotics increased with area across scales but decreased within each scale. The foreign exotics exhibited the highest dominance (lowest evenness) and spatial variation in distribution, followed by domestic exotics and natives, although on average natives still occupy larger areas than domestic and foreign exotics. The results provide new insight into patterns and mechanisms of fish species invasions at multiple spatial scales in the United States.     

A Dated Phylogeny Complements Macroecological Analysis to Explain the Diversity Patterns in Geonoma (Arecaceae)

Integrating phylogenetic data into macroecological studies of biodiversity patterns may complement the information provided by present‐day spatial patterns. In the present study, we used range map data for all Geonoma (Arecaceae) species to assess whether Geonoma species composition forms spatially coherent floristic clusters. We then evaluated the extent to which the spatial variation in species composition reflects present‐day environmental variation vs. nonenvironmental spatial effects, as expected if the pattern reflects historical biogeography. We also examined the degree of geographic structure in the Geonoma phylogeny. Finally, we used a dated phylogeny to assess whether species richness within the floristic clusters was constrained by a specific historical biogeographic driver, namely time‐for‐diversification. A cluster analysis identified six spatially coherent floristic clusters, four of which were used to reveal a significant geographic phylogenetic structure. Variation partitioning analysis showed that 56 percent of the variation in species composition could be explained by spatial variables alone, consistent with historical factors having played a major role in generating the Geonoma diversity pattern. To test for a time‐for‐diversification effect, we correlated four different species richness measures with the diversification time of the earliest large lineage that is characteristic of each cluster. In support of this hypothesis, we found that geographic areas with higher richness contained older radiations. We conclude that current geographic diversity patterns in Geonoma reflect the present‐day climate, but to a larger extent are related to nonenvironmental spatial constraints linked to colonization time, dispersal limitation, and geological history, followed by within‐area evolutionary diversification. Abstract in Spanish is available at http://www.blackwell‐synergy.com/loi/btp .     

Global patterns of amphibian phylogenetic diversity

Aim Phylogenetic diversity can provide insight into how evolutionary processes may have shaped contemporary patterns of species richness. Here, we aim to test for the influence of phylogenetic history on global patterns of amphibian species richness, and to identify areas where macroevolutionary processes such as diversification and dispersal have left strong signatures on contemporary species richness. Location Global; equal‐area grid cells of approximately 10,000 km². Methods We generated an amphibian global supertree (6111 species) and repeated analyses with the largest available molecular phylogeny (2792 species). We combined each tree with global species distributions to map four indices of phylogenetic diversity. To investigate congruence between global spatial patterns of amphibian species richness and phylogenetic diversity, we selected Faith’s phylogenetic diversity (PD) index and the total taxonomic distinctness (TTD) index, because we found that the variance of the other two indices we examined (average taxonomic distinctness and mean root distance) strongly depended on species richness. We then identified regions with unusually high or low phylogenetic diversity given the underlying level of species richness by using the residuals from the global relationship of species richness and phylogenetic diversity. Results Phylogenetic diversity as measured by either Faith’s PD or TTD was strongly correlated with species richness globally, while the other two indices showed very different patterns. When either Faith’s PD or TTD was tested against species richness, residuals were strongly spatially structured. Areas with unusually low phylogenetic diversity for their associated species richness were mostly on islands, indicating large radiations of few lineages that have successfully colonized these archipelagos. Areas with unusually high phylogenetic diversity were located around biogeographic contact zones in Central America and southern China, and seem to have experienced high immigration or in situ diversification rates, combined with local persistence of old lineages. Main conclusions We show spatial structure in the residuals of the relationship between species richness and phylogenetic diversity, which together with the positive relationship itself indicates strong signatures of evolutionary history on contemporary global patterns of amphibian species richness. Areas with unusually low and high phylogenetic diversity for their associated richness demonstrate the importance of biogeographic barriers to dispersal, colonization and diversification processes.     

The role of spatial scale and the perception of large-scale species-richness patterns

Despite two centuries of exploration, our understanding of factors determining the distribution of life on Earth is in many ways still in its infancy. Much of the disagreement about governing processes of variation in species richness may be the result of differences in our perception of species‐richness patterns. Until recently, most studies of large‐scale species‐richness patterns assumed implicitly that patterns and mechanisms were scale invariant. Illustrated with examples and a quantitative analysis of published data on altitudinal gradients of species richness (n = 204), this review discusses how scale effects (extent and grain size) can influence our perception of patterns and processes. For example, a hump‐shaped altitudinal species‐richness pattern is the most typical (c. 50%), with a monotonic decreasing pattern (c. 25%) also frequently reported, but the relative distribution of patterns changes readily with spatial grain and extent. If we are to attribute relative impact to various factors influencing species richness and distribution and to decide at which point along a spatial and temporal continuum they act, we should not ask only how results vary as a function of scale but also search for consistent patterns in these scale effects. The review concludes with suggestions of potential routes for future analytical exploration of species‐richness patterns.     

Out of the dark: 350 million years of conservatism and evolution in diel activity patterns in vertebrates

Many animals are active only during a particular time (e.g., day vs. night), a partitioning that may have important consequences for species coexistence. An open question is the extent to which this diel activity niche is evolutionarily conserved or labile. Here, we analyze diel activity data across a phylogeny of 1914 tetrapod species. We find strong phylogenetic signal, showing that closely related species tend to share similar activity patterns. Ancestral reconstructions show that nocturnality was the most likely ancestral diel activity pattern for tetrapods and many major clades within it (e.g., amphibians, mammals). Remarkably, nocturnal activity appears to have been maintained continuously in some lineages for ～350 million years. Thus, we show that traits involved in local‐scale resource partitioning can be conserved over strikingly deep evolutionary time scales. We also demonstrate a potentially important (but often overlooked) metric of niche conservatism. Finally, we show that diurnal lineages appear to have faster speciation and diversification rates than nocturnal lineages, which may explain why there are presently more diurnal tetrapod species even though diurnality appears to have evolved more recently. Overall, our results may have implications for studies of community ecology, species richness, and the evolution of diet and communication systems.     

Global tests of regional effect on species richness of vascular plants and terrestrial vertebrates

Determining the effects of regional and contemporary factors on large-scale patterns in species richness has been a fundamental question in modern ecology and biogeography. However, few studies have examined effects of historical and regional factors on species richness at the global scale, and conclusions are often inconsistent or controversial. Here, I use a comprehensive dataset to examine regional effects on species richness of vascular plants and four taxa of terrestrial vertebrates (mammals, birds, reptiles, and amphibians) in the same set of sample units (ecoregions) in seven biogeographic realms across the globe. The same spatial scale and the same set of environmental variables, which are thought to influence large-scale patterns in species richness of vascular plants and terrestrial vertebrates, are used for all the five taxa. Species richness of each taxon is compared across biogeographic realms. Regional effect on species richness has been found for all the five taxa. Of the 90 realm-pair comparisons for the five taxa between observed richness of a region and the richness of the region predicted by the richness–environment relationship derived from the data of another region, 74 (82.2%) showed significant differences between observed and predicted species richness, indicating that a species richness–environment relationship developed for one region cannot accurately predict species richness in other regions of similar environments.     

Mutualistic mimicry enhances species diversification through spatial segregation and extension of the ecological niche space

Species richness varies among clades, yet the drivers of diversification creating this variation remain poorly understood. While abiotic factors likely drive some of the variation in species richness, ecological interactions may also contribute. Here, we examine one class of potential contributors to species richness variation that is particularly poorly understood: mutualistic interactions. We aim to elucidate large‐scale patterns of diversification mediated by mutualistic interactions using a spatially explicit population‐based model. We focus on mutualistic Müllerian mimicry between conspicuous toxic prey species, where convergence in color patterns emerges from predators' learning process. To investigate the effects of Müllerian mimicry on species diversification, we assume that some speciation events stem from shifts in ecological niches, and can also be associated with shift in mimetic color pattern. Through the emergence of spatial mosaics of mimetic color patterns, Müllerian mimicry constrains the geographical distribution of species and allows different species occupying similar ecological niches to exist simultaneously in different regions. Müllerian mimicry and the resulting spatial segregation of mimetic color patterns thus generate more balanced phylogenetic trees and increase overall species diversity. Our study sheds light on complex effects of Müllerian mimicry on ecological, spatial, and phylogenetic diversification.