Species Richness and Evolutionary Niche Dynamics: A Spatial Pattern–Oriented Simulation Experiment

Elevational Goldilocks zone underlies the exceptional diversity of a large lizard radiation (Liolaemus; Liolaemidae)

Mountains are among the most biodiverse regions on the planet, and how these landforms shape diversification through the interaction of biological traits and geo-climatic dynamics is integral to understanding global biodiversity. In this study, we investigate the dual roles of climate change and mountain uplift on the evolution of a hyper-diverse radiation, Liolaemus lizards, with a spatially explicit model of diversification using a reconstruction of uplift and paleotemperature in central and southern South America. The diversification model captures a hotspot for Liolaemus around 40°S in lineages with low-dispersal ability and narrow niche breadths. Under the model, speciation rates are highest in low latitudes (<35°S) and mid elevations (~1,000 m), while extinction rates are highest at higher latitudes (>35°S) and higher elevations (>2,000 m). Temperature change through the Cenozoic explained variation in speciation and extinction rates through time and across different elevational bands. Our results point to the conditions of mid elevations being optimal for diversification (i.e., Goldilocks Zone), driven by the combination of (1) a complex topography that facilitates speciation during periods of climatic change, and (2) a relatively moderate climate that enables the persistence of ectothermic lineages and buffers species from extinction.

Potential present and future distributions of the genus Atta of Mexico

Temperature and precipitation influence insect distribution locally and drive large-scale biogeographical patterns. We used current and future climate data from the CHELSA database to create ensemble species distribution models for three Atta leaf-cutting ant species (Atta cephalotes, A. mexicana, and A. texana) found in Mexico. These models were used to estimate the potential impact of climate change on the distribution of these species in the future. Our results show that bioclimatic variables influence the distribution of each Atta species occupying a unique climatic niche: A. cephalotes is affected by temperature seasonality, A. mexicana by isothermality, and A. texana by the minimum temperature of the coldest month. Atta texana and A. mexicana are expected to decline their range by 80% and 60%, respectively, due to rising temperatures, decreased rainfall, and increased drought. Due to rising temperatures and increased humidity, Atta cephalotes is expected to expand its range by 30%. Since Atta species are important pests, our coexistence with them requires knowledge of their ecological functions and potential future distribution changes. In addition, these insects serve as bioindicators of habitat quality, and they can contribute to the local economy in rural areas since they are eaten as food for the nutritional value of the queens. In this sense, presenting a future perspective of these species' distribution is important for forest and crop management. Education programs also are necessary to raise awareness of the importance of these ants and the challenges they face because of climate change. Our results offer a perspective of climate change studies to define conservation and adaptation strategies for protecting vulnerable areas such as high-elevation remnant forests.

Attributes of host-specificity better explain the diversified wood-boring longhorn beetles in tropical SW China than plant species diversity

A long-debated question in ecology is whether the hyper-diversity of tropical plant-feeding insects is a direct consequence of high tropical plant diversity or should be attributed to increases in host plant specialization. In this study, we used Cerambycidae (the wood-boring longhorn beetles whose larval stages feed on the xylems of trees and lianas) and plants as study materials to explore which hypothesis is more favoured. Multiple analyses were used to show the differences in host specificity of Cerambycidae in tropical and subtropical forests. From these analyses, we found that the alpha diversity of beetles in tropical forests was significantly higher than that in subtropical forests but not in plants. The relationship between plants and beetles was also closer in tropical areas than in subtropical areas. Our results imply that the wood-boring longhorn beetles show higher degrees of niche conservatism and host-specificity in tropical forests than in subtropical forests. The high diversity of wood-boring longhorn beetles in tropical forests might be explained to a large extent by their more finely partitioned diet breadth.

A quixotic view of spatial bias in modelling the distribution of species and their diversity

Ecological processes are often spatially and temporally structured, potentially leading to autocorrelation either in environmental variables or species distribution data. Because of that, spatially-biased in-situ samples or predictors might affect the outcomes of ecological models used to infer the geographic distribution of species and diversity. There is a vast heterogeneity of methods and approaches to assess and measure spatial bias; this paper aims at addressing the spatial component of data-driven biases in species distribution modelling, and to propose potential solutions to explicitly test and account for them. Our major goal is not to propose methods to remove spatial bias from the modelling procedure, which would be impossible without proper knowledge of all the processes generating it, but rather to propose alternatives to explore and handle it. In particular, we propose and describe three main strategies that may provide a fair account of spatial bias, namely: (i) how to represent spatial bias; (ii) how to simulate null models based on virtual species for testing biogeographical and species distribution hypotheses; and (iii) how to make use of spatial bias - in particular related to sampling effort - as a leverage instead of a hindrance in species distribution modelling. We link these strategies with good practice in accounting for spatial bias in species distribution modelling.

Process-explicit models reveal the structure and dynamics of biodiversity patterns

With ever-growing data availability and computational power at our disposal, we now have the capacity to use process-explicit models more widely to reveal the ecological and evolutionary mechanisms responsible for spatiotemporal patterns of biodiversity. Most research questions focused on the distribution of diversity cannot be answered experimentally, because many important environmental drivers and biological constraints operate at large spatiotemporal scales. However, we can encode proposed mechanisms into models, observe the patterns they produce in virtual environments, and validate these patterns against real-world data or theoretical expectations. This approach can advance understanding of generalizable mechanisms responsible for the distributions of organisms, communities, and ecosystems in space and time, advancing basic and applied science. We review recent developments in process-explicit models and how they have improved knowledge of the distribution and dynamics of life on Earth, enabling biodiversity to be better understood and managed through a deeper recognition of the processes that shape genetic, species, and ecosystem diversity.

Species-range-size distributions: Integrating the effects of speciation, transformation, and extinction

The species-range size distribution is a product of speciation, transformation of range-sizes, and extinction. Previous empirical studies showed that it has a left-skewed lognormal-like distribution. We developed a new mathematical framework to study species-range-size distributions, one in which allopatric speciation, transformation of range size, and the extinction process are explicitly integrated. The approach, which we call the gain-loss-allopatric speciation model, allows us to explore the effects of various speciation scenarios. Our model captures key dynamics thought to lead to known range-size distributions. We also fitted the model to empirical range-size distributions of birds, mammals, and beetles. Since geographic range dynamics are linked to speciation and extinction, our model provides predictions for the dynamics of species richness. When a species-range-size distribution initially evolves away from the range sizes at which the likelihood of speciation is low, it tends to cause diversification slowdown even in the absence of (bio)diversity dependence in speciation rate. Using the mathematical model developed here, we give a potential explanation for how observed range-size distributions emerge from range-size dynamics. Although the framework presented is minimalistic, it provides a starting point for examining hypotheses based on more complex mechanisms.

gen3sis: A general engine for eco-evolutionary simulations of the processes that shape Earth's biodiversity

Understanding the origins of biodiversity has been an aspiration since the days of early naturalists. The immense complexity of ecological, evolutionary, and spatial processes, however, has made this goal elusive to this day. Computer models serve progress in many scientific fields, but in the fields of macroecology and macroevolution, eco-evolutionary models are comparatively less developed. We present a general, spatially explicit, eco-evolutionary engine with a modular implementation that enables the modeling of multiple macroecological and macroevolutionary processes and feedbacks across representative spatiotemporally dynamic landscapes. Modeled processes can include species' abiotic tolerances, biotic interactions, dispersal, speciation, and evolution of ecological traits. Commonly observed biodiversity patterns, such as α, β, and γ diversity, species ranges, ecological traits, and phylogenies, emerge as simulations proceed. As an illustration, we examine alternative hypotheses expected to have shaped the latitudinal diversity gradient (LDG) during the Earth's Cenozoic era. Our exploratory simulations simultaneously produce multiple realistic biodiversity patterns, such as the LDG, current species richness, and range size frequencies, as well as phylogenetic metrics. The model engine is open source and available as an R package, enabling future exploration of various landscapes and biological processes, while outputs can be linked with a variety of empirical biodiversity patterns. This work represents a key toward a numeric, interdisciplinary, and mechanistic understanding of the physical and biological processes that shape Earth's biodiversity.

From theory to practice in pattern-oriented modelling: identifying and using empirical patterns in predictive models

To robustly predict the effects of disturbance and ecosystem changes on species, it is necessary to produce structurally realistic models with high predictive power and flexibility. To ensure that these models reflect the natural conditions necessary for reliable prediction, models must be informed and tested using relevant empirical observations. Pattern-oriented modelling (POM) offers a systematic framework for employing empirical patterns throughout the modelling process and has been coupled with complex systems modelling, such as in agent-based models (ABMs). However, while the production of ABMs has been rising rapidly, the explicit use of POM has not increased. Challenges with identifying patterns and an absence of specific guidelines on how to implement empirical observations may limit the accessibility of POM and lead to the production of models which lack a systematic consideration of reality. This review serves to provide guidance on how to identify and apply patterns following a POM approach in ABMs (POM-ABMs), specifically addressing: where in the ecological hierarchy can we find patterns; what kinds of patterns are useful; how should simulations and observations be compared; and when in the modelling cycle are patterns used? The guidance and examples provided herein are intended to encourage the application of POM and inspire efficient identification and implementation of patterns for both new and experienced modellers alike. Additionally, by generalising patterns found especially useful for POM-ABM development, these guidelines provide practical help for the identification of data gaps and guide the collection of observations useful for the development and verification of predictive models. Improving the accessibility and explicitness of POM could facilitate the production of robust and structurally realistic models in the ecological community, contributing to the advancement of predictive ecology at large.

The mathematical influence on global patterns of biodiversity

Although we understand how species evolve, we do not appreciate how this process has filled an empty world to create current patterns of biodiversity. Here, we conduct a numerical experiment to determine why biodiversity varies spatially on our planet. We show that spatial patterns of biodiversity are mathematically constrained and arise from the interaction between the species' ecological niches and environmental variability that propagates to the community level. Our results allow us to explain key biological observations such as (a) latitudinal biodiversity gradients (LBGs) and especially why oceanic LBGs primarily peak at midlatitudes while terrestrial LBGs generally exhibit a maximum at the equator, (b) the greater biodiversity on land even though life first evolved in the sea, (c) the greater species richness at the seabed than at the sea surface, and (d) the higher neritic (i.e., species occurring in areas with a bathymetry lower than 200 m) than oceanic (i.e., species occurring in areas with a bathymetry higher than 200 m) biodiversity. Our results suggest that a mathematical constraint originating from a fundamental ecological interaction, that is, the niche-environment interaction, fixes the number of species that can establish regionally by speciation or migration.

Climatic-niche evolution follows similar rules in plants and animals

Climatic niches are essential in determining where species can occur and how they will respond to climate change. However, it remains unclear if climatic-niche evolution is similar in plants and animals or is intrinsically different. For example, previous authors have proposed that plants have broader environmental tolerances than animals but are more sensitive to climate change. Here, we test ten predictions about climatic-niche evolution in plants and animals, using phylogenetic and climatic data for 19 plant clades and 17 vertebrate clades (2,087 species total). Surprisingly, we find that for all ten predictions, plants and animals show similar patterns. For example, in both groups, climatic niches change at similar mean rates and species have similar mean niche breadths, and niche breadths show similar relationships with latitude across groups. Our results suggest that there are general 'rules' of climatic-niche evolution that span plants and animals, despite the fundamental differences in their biology. These results may help to explain why plants and animals have similar responses to climate change and why they often have shared species richness patterns, biogeographic regions, biomes and biodiversity hotspots.

Dispersion fields reveal the compositional structure of South American vertebrate assemblages

The causes of continental patterns in species richness continue to spur heated discussion. Hypotheses based on ambient energy have dominated the debate, but are increasingly being challenged by hypotheses that model richness as the overlap of species ranges, ultimately controlled by continental range dynamics of individual species. At the heart of this controversy lies the question of whether species richness of individual grid cells is controlled by local factors, or reflects larger-scale spatial patterns in the turnover of species’ ranges. Here, we develop a new approach based on assemblage dispersion fields, formed by overlaying the geographic ranges of all species co-occurring in a grid cell. We created dispersion fields for all tetrapods of South America, and characterized the orientation and shape of dispersion fields as a vector field. The resulting maps demonstrate the existence of macro-structures in the turnover of biotic similarity at continental scale that are congruent among vertebrate classes. These structures underline the importance of continental-scale processes for species richness in individual assemblages.

Ecologists continue to debate whether local species assemblages result from habitat filtering or from turnover among the regional species pool. Here the authors develop a “dispersion field” method to mapping species range overlaps, showing that regional turnover processes are key to local assembly.

A Grinnellian Niche Perspective on Species-Area Relationships

In this work, Grinnellian niche theory (a body of theory about geographic distributions of species in terms of noninteracting niche variables) is used to demonstrate that species-area relationships emerge with both size of environmental space and size of geographic area. As environmental space increases, more species' fundamental niches are included, thus increasing the number of species capable of living in the corresponding region. This idea is made operational by proposing a size measure for multidimensional environmental space and approximating fundamental niches with minimum volume ellipsoids. This framework allows estimating a presence-absence matrix based on the distribution of fundamental niches in environmental space, from which many biodiversity measures can be calculated, such as beta diversity. I establish that Whittaker's equation for beta diversity is equivalent to MacArthur's formula relating species numbers and niche breadth; this latter equation provides a mechanism for the species-niche space relationship. I illustrate the theoretical results via exploration of niches of the terrestrial mammals of North America (north of Panama). Each world region has a unique structure of its environmental space, and the position of fundamental niches in niche space is different for different clades; therefore, species-area relationships depend on the clades involved and the region of focus, mostly as a function of MacArthur's niche beta diversity. Analyzing species-area relationships from the perspective of niche position in environmental space is novel, shifting emphasis from demographic processes to historical, geographic, and climatic factors; moreover, the Grinnellian approach is based on available data and is computationally feasible.

Reconstructing Ecological Niche Evolution When Niches Are Incompletely Characterized

Evolutionary dynamics of abiotic ecological niches across phylogenetic history can shed light on large-scale biogeographic patterns, macroevolutionary rate shifts, and the relative ability of lineages to respond to global change. An unresolved question is how best to represent and reconstruct evolution of these complex traits at coarse spatial scales through time. Studies have approached this question by integrating phylogenetic comparative methods with niche estimates inferred from correlative and other models. However, methods for estimating niches often produce incomplete characterizations, as they are inferred from present-day distributions that may be limited in full expression of the fundamental ecological niche by biotic interactions, dispersal limitations, and the existing set of environmental conditions. Here, we test whether incomplete niche characterizations inherent in most estimates of species' niches bias phylogenetic reconstructions of niche evolution, using simulations of virtual species with known niches. Results establish that incompletely characterized niches inflate estimates of evolutionary change and lead to error in ancestral state reconstructions. Our analyses also provide a potential mechanism to explain the frequent observation that maximum thermal tolerances are more conserved than minimum thermal tolerances: populations and species experience more spatial variation in minimum temperature than in maximum temperature across their distributions and, consequently, may experience stronger diversifying selection for cold tolerance.

Modeling the ecology and evolution of biodiversity: Biogeographical cradles, museums, and graves

Individual processes shaping geographical patterns of biodiversity are increasingly understood, but their complex interactions on broad spatial and temporal scales remain beyond the reach of analytical models and traditional experiments. To meet this challenge, we built a spatially explicit, mechanistic simulation model implementing adaptation, range shifts, fragmentation, speciation, dispersal, competition, and extinction, driven by modeled climates of the past 800,000 years in South America. Experimental topographic smoothing confirmed the impact of climate heterogeneity on diversification. The simulations identified regions and episodes of speciation (cradles), persistence (museums), and extinction (graves). Although the simulations had no target pattern and were not parameterized with empirical data, emerging richness maps closely resembled contemporary maps for major taxa, confirming powerful roles for evolution and diversification driven by topography and climate.

Comparing spatial diversification and meta-population models in the Indo-Australian Archipelago

Reconstructing the processes that have shaped the emergence of biodiversity gradients is critical to understand the dynamics of diversification of life on Earth. Islands have traditionally been used as model systems to unravel the processes shaping biological diversity. MacArthur and Wilson's island biogeographic model predicts diversity to be based on dynamic interactions between colonization and extinction rates, while treating islands themselves as geologically static entities. The current spatial configuration of islands should influence meta-population dynamics, but long-term geological changes within archipelagos are also expected to have shaped island biodiversity, in part by driving diversification. Here, we compare two mechanistic models providing inferences on species richness at a biogeographic scale: a mechanistic spatial-temporal model of species diversification and a spatial meta-population model. While the meta-population model operates over a static landscape, the diversification model is driven by changes in the size and spatial configuration of islands through time. We compare the inferences of both models to floristic diversity patterns among land patches of the Indo-Australian Archipelago. Simulation results from the diversification model better matched observed diversity than a meta-population model constrained only by the contemporary landscape. The diversification model suggests that the dynamic re-positioning of islands promoting land disconnection and reconnection induced an accumulation of particularly high species diversity on Borneo, which is central within the island network. By contrast, the meta-population model predicts a higher diversity on the mainlands, which is less compatible with empirical data. Our analyses highlight that, by comparing models with contrasting assumptions, we can pinpoint the processes that are most compatible with extant biodiversity patterns.

Biodiversity Models: What If Unsaturation Is the Rule?

Improving biodiversity predictions is essential if we are to meet the challenges posed by global change. As knowledge is key to feed models, we need to evaluate how debated theory can affect models. An important ongoing debate is whether environmental constraints limit the number of species that can coexist in a community (saturation), with recent findings suggesting that species richness in many communities might be unsaturated. Here, we propose that biodiversity models could address this issue by accounting for a duality: considering communities as unsaturated but where species composition is constrained by different scale-dependent biodiversity drivers. We identify a variety of promising advances for incorporating this duality into commonly applied biodiversity modelling approaches and improving their spatial predictions.

The majority of biodiversity modelling approaches do not explicitly address the question of saturation.

Theoretical and methodological implications of saturation or unsaturation in biodiversity modelling.

Addressing saturation or unsaturation is vital to produce more reliable conservation strategies.

Integrative community modelling frameworks may be the way forward.

Environmental niche divergence among three dune shrub sister species with parapatric distributions

BACKGROUND AND AIMS

The geographical distributions of species are constrained by their ecological requirements. The aim of this work was to analyse the effects of environmental conditions, historical events and biogeographical constraints on the diversification of the three species of the western Mediterranean shrub genus Stauracanthus , which have a parapatric distribution in the Iberian Peninsula.

METHODS

Ecological niche factor analysis and generalized linear models were used to measure the response of all Stauracanthus species to the environmental gradients and map their potential distributions in the Iberian Peninsula. The bioclimatic niche overlap between the three species was determined by using Schoener's index. The genetic differentiation of the Iberian and northern African populations of Stauracanthus species was characterized with GenalEx. The effects on genetic distances of the most important environmental drivers were assessed through Mantel tests and non-metric multidimensional scaling.

KEY RESULTS

The three Stauracanthus species show remarkably similar responses to climatic conditions. This supports the idea that all members of this recently diversified clade retain common adaptations to climate and consequently high levels of climatic niche overlap. This contrasts with the diverse edaphic requirements of Stauracanthus species. The populations of the S. genistoides-spectabilis clade grow on Miocene and Pliocene fine-textured sedimentary soils, whereas S. boivinii , the more genetically distant species, occurs on older and more coarse-textured sedimentary substrates. These patterns of diversification are largely consistent with a stochastic process of geographical range expansion and fragmentation coupled with niche evolution in the context of spatially complex environmental fluctuations.

CONCLUSIONS

: The combined analysis of the distribution, realized environmental niche and phylogeographical relationships of parapatric species proposed in this work allows integration of the biogeographical, ecological and evolutionary processes driving the evolution of species adaptations and how they determine their current geographical ranges.

The impact of Cenozoic cooling on assemblage diversity in planktonic foraminifera

The Cenozoic planktonic foraminifera (PF) (calcareous zooplankton) have arguably the most detailed fossil record of any group. The quality of this record allows models of environmental controls on macroecology, developed for Recent assemblages, to be tested on intervals with profoundly different climatic conditions. These analyses shed light on the role of long-term global cooling in establishing the modern latitudinal diversity gradient (LDG)--one of the most powerful generalizations in biogeography and macroecology. Here, we test the transferability of environment-diversity models developed for modern PF assemblages to the Eocene epoch (approx. 56-34 Ma), a time of pronounced global warmth. Environmental variables from global climate models are combined with Recent environment-diversity models to predict Eocene richness gradients, which are then compared with observed patterns. The results indicate the modern LDG--lower richness towards the poles--developed through the Eocene. Three possible causes are suggested for the mismatch between statistical model predictions and data in the Early Eocene: the environmental estimates are inaccurate, the statistical model misses a relevant variable, or the intercorrelations among facets of diversity--e.g. richness, evenness, functional diversity--have changed over geological time. By the Late Eocene, environment-diversity relationships were much more similar to those found today.

Evolution of climatic niche specialization: a phylogenetic analysis in amphibians

The evolution of climatic niche specialization has important implications for many topics in ecology, evolution and conservation. The climatic niche reflects the set of temperature and precipitation conditions where a species can occur. Thus, specialization to a limited set of climatic conditions can be important for understanding patterns of biogeography, species richness, community structure, allopatric speciation, spread of invasive species and responses to climate change. Nevertheless, the factors that determine climatic niche width (level of specialization) remain poorly explored. Here, we test whether species that occur in more extreme climates are more highly specialized for those conditions, and whether there are trade-offs between niche widths on different climatic niche axes (e.g. do species that tolerate a broad range of temperatures tolerate only a limited range of precipitation regimes?). We test these hypotheses in amphibians, using phylogenetic comparative methods and global-scale datasets, including 2712 species with both climatic and phylogenetic data. Our results do not support either hypothesis. Rather than finding narrower niches in more extreme environments, niches tend to be narrower on one end of a climatic gradient but wider on the other. We also find that temperature and precipitation niche breadths are positively related, rather than showing trade-offs. Finally, our results suggest that most amphibian species occur in relatively warm and dry environments and have relatively narrow climatic niche widths on both of these axes. Thus, they may be especially imperilled by anthropogenic climate change.

Spatial predictions at the community level: from current approaches to future frameworks

A fundamental goal of ecological research is to understand and model how processes generate patterns so that if conditions change, changes in the patterns can be predicted. Different approaches have been proposed for modelling species assemblage, but their use to predict spatial patterns of species richness and other community attributes over a range of spatial and temporal scales remains challenging. Different methods emphasize different processes of structuring communities and different goals. In this review, we focus on models that were developed for generating spatially explicit predictions of communities, with a particular focus on species richness, composition, relative abundance and related attributes. We first briefly describe the concepts and theories that span the different drivers of species assembly. A combination of abiotic processes and biotic mechanisms are thought to influence the community assembly process. In this review, we describe four categories of drivers: (i) historical and evolutionary, (ii) environmental, (iii) biotic, and (iv) stochastic. We discuss the different modelling approaches proposed or applied at the community level and examine them from different standpoints, i.e. the theoretical bases, the drivers included, the source data, and the expected outputs, with special emphasis on conservation needs under climate change. We also highlight the most promising novelties, possible shortcomings, and potential extensions of existing methods. Finally, we present new approaches to model and predict species assemblages by reviewing promising 'integrative frameworks' and views that seek to incorporate all drivers of community assembly into a unique modelling workflow. We discuss the strengths and weaknesses of these new solutions and how they may hasten progress in community-level modelling.

Estimating and mapping ecological processes influencing microbial community assembly

Ecological community assembly is governed by a combination of (i) selection resulting from among-taxa differences in performance; (ii) dispersal resulting from organismal movement; and (iii) ecological drift resulting from stochastic changes in population sizes. The relative importance and nature of these processes can vary across environments. Selection can be homogeneous or variable, and while dispersal is a rate, we conceptualize extreme dispersal rates as two categories; dispersal limitation results from limited exchange of organisms among communities, and homogenizing dispersal results from high levels of organism exchange. To estimate the influence and spatial variation of each process we extend a recently developed statistical framework, use a simulation model to evaluate the accuracy of the extended framework, and use the framework to examine subsurface microbial communities over two geologic formations. For each subsurface community we estimate the degree to which it is influenced by homogeneous selection, variable selection, dispersal limitation, and homogenizing dispersal. Our analyses revealed that the relative influences of these ecological processes vary substantially across communities even within a geologic formation. We further identify environmental and spatial features associated with each ecological process, which allowed mapping of spatial variation in ecological-process-influences. The resulting maps provide a new lens through which ecological systems can be understood; in the subsurface system investigated here they revealed that the influence of variable selection was associated with the rate at which redox conditions change with subsurface depth.

Geographic range shifts do not erase the historic signal of speciation in mammals

Many evolutionary analyses assume that the positions of species geographic ranges are sufficiently phylogenetically conserved that current ranges reflect ancestral ranges and retain the historic signal of speciation. The validity of this assumption has been challenged, because there is evidence that ranges can shift rapidly and extensively. Here I test the assumption of range conservatism using simulations and empirical tests of phylogenetic signal in geographic positions of ranges within mammal orders, families, and genera. In most taxa, range positions show strong phylogenetic signal, quantified using Pagel's λ, Mantel tests, and a novel method to measure phylogenetic signal near the tips of a phylogeny. Taxa with highly labile range positions are exceptions to the general pattern and include very young groups such as Sciurus that may still be in the early, rapid-expansion phase of adaptive radiation. In two orders containing many species with large distributions (Artiodactyla and Carnivora), temporal patterns of range evolution are consistent with large instantaneous shifts in range position associated with allopatric speciation. In most other taxa, range evolution is better described by models that allow ranges to evolve along branches of the phylogeny. The results point to a common pattern of phylogenetically conserved ranges where the current position of species ranges reflects their position at the time of speciation, modified by gradual drift of range boundaries through time.

Rapoport's rule revisited: geographical distributions of human languages

One of the most well studied ecological patterns is Rapoport's rule, which posits that the geographical extent of species ranges increases at higher latitudes. However, studies to date have been limited in their geographic scope and results have been equivocal. In turn, much debate exists over potential links between Rapoport's rule and latitudinal patterns in species richness. Humans collectively speak nearly 7000 different languages, which are spread unevenly across the globe, with loci in the tropics. Causes of this skewed distribution have received only limited study. We analyze the extent of Rapoport's rule in human languages at a global scale and within each region of the globe separately. We test the relationship between Rapoport's rule and the richness of languages spoken in different regions. We also explore the frequency distribution of language-range sizes. The language-range area distribution is strongly right-skewed, with 87% of languages having range areas less than 10,000 km(2), and only nine languages with range areas over 1,000,000 km(2). At a global scale, language-range extents and areas are positively correlated with latitude. At a global scale and in five of the six regions examined, language-range extent and language-range area are strongly correlated with language richness. Our results point to group boundary formation as a critical mediator of the relationship between Rapoport's rule and diversity patterns. Where strong group boundaries limit range overlap, as is the case with human languages, and range sizes increase with latitude, latitudinal richness gradients may result.

Stronger tests of mechanisms underlying geographic gradients of biodiversity: insights from the dimensionality of biodiversity

Inference involving diversity gradients typically is gathered by mechanistic tests involving single dimensions of biodiversity such as species richness. Nonetheless, because traits such as geographic range size, trophic status or phenotypic characteristics are tied to a particular species, mechanistic effects driving broad diversity patterns should manifest across numerous dimensions of biodiversity. We develop an approach of stronger inference based on numerous dimensions of biodiversity and apply it to evaluate one such putative mechanism: the mid-domain effect (MDE). Species composition of 10,000-km(2) grid cells was determined by overlaying geographic range maps of 133 noctilionoid bat taxa. We determined empirical diversity gradients in the Neotropics by calculating species richness and three indices each of phylogenetic, functional and phenetic diversity for each grid cell. We also created 1,000 simulated gradients of each examined metric of biodiversity based on a MDE model to estimate patterns expected if species distributions were randomly placed within the Neotropics. For each simulation run, we regressed the observed gradient onto the MDE-expected gradient. If a MDE drives empirical gradients, then coefficients of determination from such an analysis should be high, the intercept no different from zero and the slope no different than unity. Species richness gradients predicted by the MDE fit empirical patterns. The MDE produced strong spatially structured gradients of taxonomic, phylogenetic, functional and phenetic diversity. Nonetheless, expected values generated from the MDE for most dimensions of biodiversity exhibited poor fit to most empirical patterns. The MDE cannot account for most empirical patterns of biodiversity. Fuller understanding of latitudinal gradients will come from simultaneous examination of relative effects of random, environmental and historical mechanisms to better understand distribution and abundance of the current biota.

Diversity and niche evolution along aridity gradients in north american lizards (phrynosomatidae)

Deserts occupy approximately 12% of the Earth's land surface, and are thought to have species poor but highly specialized biotas. However, few studies have examined the evolutionary origins of desert biotas and of diversity patterns along aridity gradients. Further, it is unclear if species occurring in more extreme conditions on a given niche axis (i.e., precipitation) are more specialized for those conditions (i.e., have narrower niche breadths). We address these questions here using a time-calibrated phylogeny and climatic data for 117 species of phrynosomatid lizards. Phrynosomatids are the most species-rich family of lizards in North America, and are found from deserts to rainforests. Surprisingly, we find that phrynosomatids have higher richness in more arid environments. This pattern occurs seemingly because they have been present in more arid habitats longer (~55 million years), and lineages in mesic environments are recently derived from more arid-dwelling ancestors. We find little support for the hypothesis that species in more extreme environments are more specialized. Instead, many desert-dwelling species are broadly distributed, and species in the most mesic environments have the broadest niche breadths. In summary, phrynosomatids offer a counterexample to the idea that arid regions are inhabited by a small number of recent and highly specialized lineages.

Phylogenetic fields of species: cross-species patterns of phylogenetic structure and geographical coexistence

Differential coexistence among species underlies geographical patterns of biodiversity. Understanding such patterns has relied either on ecological or historical approaches applied separately. Recently, macroecology and community phylogenetics have tried to integrate both ecological and historical approaches. However, macroecology is mostly non-phylogenetic, whereas community phylogenetics is largely focused on local scales. Here, we propose a conceptual framework to link macroecology and community phylogenetics by exploring the evolutionary context of large-scale species coexistence, introducing the phylogenetic field concept. This is defined as the phylogenetic structure of species co-occurrence within a focal species' geographical range. We developed concepts and methods for analysing phylogenetic fields and applied them to study coexistence patterns of the bat family Phyllostomidae. Our analyses showed that phyllostomid bats coexist mostly with closely related species, revealing a north-south gradient from overdispersed to clustered phylogenetic fields. Patterns at different phylogenetic levels (i.e. all species versus close relatives only) presented the same gradient. Results support the tropical niche conservatism hypothesis, potentially mediated by higher speciation rates in the region of origin coupled with shared environmental preferences among species. The phylogenetic field approach enables species-based community phylogenetics, instead of those that are site-based, allowing the description of historical processes at more appropriate macroecological and biogeographic scales.

How is the rate of climatic-niche evolution related to climatic-niche breadth?

The rate of climatic-niche evolution is important to many research areas in ecology, evolution, and conservation biology, including responses of species to global climate change, spread of invasive species, speciation, biogeography, and patterns of species richness. Previous studies have implied that clades with higher rates of climatic-niche evolution among species should have species with narrower niche breadths, but there is also evidence suggesting the opposite pattern. However, the relationships between rate and breadth have not been explicitly analyzed. Here, we examine the relationships between the rate of climatic-niche evolution and climatic-niche breadth using phylogenetic and climatic data for 250 species in the salamander family Plethodontidae, a group showing considerable variation in both rates of climatic-niche evolution and climatic-niche breadths. Contrary to some expectations, we find no general relationship between climatic-niche breadth and the rate of climatic-niche evolution. Climatic-niche breadths for some ecologically important climatic variables considered separately (temperature seasonality and annual precipitation) do show significant relationships with the rate of climatic-niche evolution, but rates are faster in clades in which species have broader (not narrower) niche breadths. In summary, our results show that narrower niche breadths are not necessarily associated with faster rates of niche evolution.

Are range-size distributions consistent with species-level heritability?

The concept of species-level heritability is widely contested. Because it is most likely to apply to emergent, species-level traits, one of the central discussions has focused on the potential heritability of geographic range size. However, a central argument against range-size heritability has been that it is not compatible with the observed shape of present-day species range-size distributions (SRDs), a claim that has never been tested. To assess this claim, we used forward simulation of range-size evolution in clades with varying degrees of range-size heritability, and compared the output of three different models to the range-size distribution of the South American avifauna. Although there were differences among the models, a moderate-to-high degree of range-size heritability consistently leads to SRDs that were similar to empirical data. These results suggest that range-size heritability can generate realistic SRDs, and may play an important role in shaping observed patterns of range sizes.