Biotic colonization of subtropical East Asian caves through time.

Caves are home to unique and fragile biotas with high levels of endemism. However, little is known about how the biotic colonization of caves has developed over time, especially in caves from middle and low latitudes. Subtropical East Asia holds the world's largest karst landform with numerous ancient caves, which harbor a high diversity of cave-dwelling organisms and are regarded as a biodiversity hotspot. Here, we assess the temporal dynamics of biotic colonization of subtropical East Asian caves through a multi-taxon analysis with representatives of green plants, animals, and fungi. We then investigate the consequences of paleonviromental changes on the colonization dynamics of these caves in combination with reconstructions of vegetation, temperature, and precipitation. We discover that 88% of cave colonization events occurred after the Oligocene-Miocene boundary, and organisms from the surrounding forest were a major source for subtropical East Asian cave biodiversity. Biotic colonization of subtropical East Asian caves during the Neogene was subject to periods of acceleration and decrease, in conjunction with large-scale, seasonal climatic changes and evolution of local forests. This study highlights the long-term evolutionary interaction between surface and cave biotas; our climate-vegetation-relict model proposed for the subtropical East Asian cave biota may help explain the evolutionary origins of other mid-latitude subterranean biotas.

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.

Earth history events shaped the evolution of uneven biodiversity across tropical moist forests.

Far from a uniform band, the biodiversity found across Earth's tropical moist forests varies widely between the high diversity of the Neotropics and Indomalaya and the relatively lower diversity of the Afrotropics. Explanations for this variation across different regions, the "pantropical diversity disparity" (PDD), remain contentious, due to difficulty teasing apart the effects of contemporary climate and paleoenvironmental history. Here, we assess the ubiquity of the PDD in over 150,000 species of terrestrial plants and vertebrates and investigate the relationship between the present-day climate and patterns of species richness. We then investigate the consequences of paleoenvironmental dynamics on the emergence of biodiversity gradients using a spatially explicit model of diversification coupled with paleoenvironmental and plate tectonic reconstructions. Contemporary climate is insufficient in explaining the PDD; instead, a simple model of diversification and temperature niche evolution coupled with paleoaridity constraints is successful in reproducing the variation in species richness and phylogenetic diversity seen repeatedly among plant and animal taxa, suggesting a prevalent role of paleoenvironmental dynamics in combination with niche conservatism. The model indicates that high biodiversity in Neotropical and Indomalayan moist forests is driven by complex macroevolutionary dynamics associated with mountain uplift. In contrast, lower diversity in Afrotropical forests is associated with lower speciation rates and higher extinction rates driven by sustained aridification over the Cenozoic. Our analyses provide a mechanistic understanding of the emergence of uneven diversity in tropical moist forests across 110 Ma of Earth's history, highlighting the importance of deep-time paleoenvironmental legacies in determining biodiversity patterns.

Diversity across organisational scale emerges through dispersal ability and speciation dynamics in tropical fish.

BACKGROUND: Biodiversity exists at different levels of organisation: e.g. genetic, individual, population, species, and community. These levels of organisation all exist within the same system, with diversity patterns emerging across organisational scales through several key processes. Despite this inherent interconnectivity, observational studies reveal that diversity patterns across levels are not consistent and the underlying mechanisms for variable continuity in diversity across levels remain elusive. To investigate these mechanisms, we apply a spatially explicit simulation model to simulate the global diversification of tropical reef fishes at both the population and species levels through emergent population-level processes. RESULTS: We find significant relationships between the population and species levels of diversity which vary depending on both the measure of diversity and the spatial partitioning considered. In turn, these population-species relationships are driven by modelled biological trait parameters, especially the divergence threshold at which populations speciate. CONCLUSIONS: To explain variation in multi-level diversity patterns, we propose a simple, yet novel, population-to-species diversity partitioning mechanism through speciation which disrupts continuous diversity patterns across organisational levels. We expect that in real-world systems this mechanism is driven by the molecular dynamics that determine genetic incompatibility, and therefore reproductive isolation between individuals. We put forward a framework in which the mechanisms underlying patterns of diversity across organisational levels are universal, and through this show how variable patterns of diversity can emerge through organisational scale.

An integrated high-resolution mapping shows congruent biodiversity patterns of Fagales and Pinales.

The documentation of biodiversity distribution through species range identification is crucial for macroecology, biogeography, conservation, and restoration. However, for plants, species range maps remain scarce and often inaccurate. We present a novel approach to map species ranges at a global scale, integrating polygon mapping and species distribution modelling (SDM). We develop a polygon mapping algorithm by considering distances and nestedness of occurrences. We further apply an SDM approach considering multiple modelling algorithms, complexity levels, and pseudo-absence selections to map the species at a high spatial resolution and intersect it with the generated polygons. We use this approach to construct range maps for all 1957 species of Fagales and Pinales with data compilated from multiple sources. We construct high-resolution global species richness maps of these important plant clades, and document diversity hotspots for both clades in southern and south-western China, Central America, and Borneo. We validate the approach with two representative genera, Quercus and Pinus, using previously published coarser range maps, and find good agreement. By efficiently producing high-resolution range maps, our mapping approach offers a new tool in the field of macroecology for studying global species distribution patterns and supporting ongoing conservation efforts.

Species ecology explains the spatial components of genetic diversity in tropical reef fishes.

Generating genomic data for 19 tropical reef fish species of the Western Indian Ocean, we investigate how species ecology influences genetic diversity patterns from local to regional scales. We distinguish between the α, ß and γ components of genetic diversity, which we subsequently link to six ecological traits. We find that the α and γ components of genetic diversity are strongly correlated so that species with a high total regional genetic diversity display systematically high local diversity. The α and γ diversity components are negatively associated with species abundance recorded using underwater visual surveys and positively associated with body size. Pelagic larval duration is found to be negatively related to genetic ß diversity supporting its role as a dispersal trait in marine fishes. Deviation from the neutral theory of molecular evolution motivates further effort to understand the processes shaping genetic diversity and ultimately the diversification of the exceptional diversity of tropical reef fishes.

A Minimal Model for the Latitudinal Diversity Gradient Suggests a Dominant Role for Ecological Limits.

The latitudinal diversity gradient (LDG) is one of Earth's most iconic biodiversity patterns and still one of the most debated. Explanations for the LDG are often categorized into three broad pathways in which the diversity gradient is created by (1) differential diversification rates, (2) differential carrying capacities (ecological limits), or (3) differential time to accumulate species across latitude. Support for these pathways has, however, been mostly verbally expressed. Here, we present a minimal model to clarify the essential assumptions of the three pathways and explore the sensitivity of diversity dynamics to these pathways. We find that an LDG arises most easily from a gradient in ecological limits compared with a gradient in the time for species accumulation or diversification rate in most modeled scenarios. Differential diversification rates create a stronger LDG than ecological limits only when speciation and dispersal rates are low, but then the predicted LDG seems weaker than the observed LDG. Moreover, range dynamics may reduce an LDG created by a gradient in diversification rates or time for species accumulation, but they cannot reduce an LDG induced by differential ecological limits. We conclude that our simple model provides a null prediction for the effectiveness of the three LDG pathways and can thus aid discussions about the causal mechanisms underlying the LDG or motivate more complex models to confirm or falsify our findings.

Estimating Age-Dependent Extinction: Contrasting Evidence from Fossils and Phylogenies.

The estimation of diversification rates is one of the most vividly debated topics in modern systematics, with considerable controversy surrounding the power of phylogenetic and fossil-based approaches in estimating extinction. Van Valen's seminal work from 1973 proposed the "Law of constant extinction," which states that the probability of extinction of taxa is not dependent on their age. This assumption of age-independent extinction has prevailed for decades with its assessment based on survivorship curves, which, however, do not directly account for the incompleteness of the fossil record, and have rarely been applied at the species level. Here, we present a Bayesian framework to estimate extinction rates from the fossil record accounting for age-dependent extinction (ADE). Our approach, unlike previous implementations, explicitly models unobserved species and accounts for the effects of fossil preservation on the observed longevity of sampled lineages. We assess the performance and robustness of our method through extensive simulations and apply it to a fossil data set of terrestrial Carnivora spanning the past 40 myr. We find strong evidence of ADE, as we detect the extinction rate to be highest in young species and declining with increasing species age. For comparison, we apply a recently developed analogous ADE model to a dated phylogeny of extant Carnivora. Although the phylogeny-based analysis also infers ADE, it indicates that the extinction rate, instead, increases with increasing taxon age. The estimated mean species longevity also differs substantially, with the fossil-based analyses estimating 2.0 myr, in contrast to 9.8 myr derived from the phylogeny-based inference. Scrutinizing these discrepancies, we find that both fossil and phylogeny-based ADE models are prone to high error rates when speciation and extinction rates increase or decrease through time. However, analyses of simulated and empirical data show that fossil-based inferences are more robust. This study shows that an accurate estimation of ADE from incomplete fossil data is possible when the effects of preservation are jointly modeled, thus allowing for a reassessment of Van Valen's model as a general rule in macroevolution.

Age-dependent speciation can explain the shape of empirical phylogenies.

Tens of thousands of phylogenetic trees, describing the evolutionary relationships between hundreds of thousands of taxa, are readily obtainable from various databases. From such trees, inferences can be made about the underlying macroevolutionary processes, yet remarkably these processes are still poorly understood. Simple and widely used evolutionary null models are problematic: Empirical trees show very different imbalance between the sizes of the daughter clades of ancestral taxa compared to what models predict. Obtaining a simple evolutionary model that is both biologically plausible and produces the imbalance seen in empirical trees is a challenging problem, to which none of the existing models provide a satisfying answer. Here we propose a simple, biologically plausible macroevolutionary model in which the rate of speciation decreases with species age, whereas extinction rates can vary quite generally. We show that this model provides a remarkable fit to the thousands of trees stored in the online database TreeBase. The biological motivation for the identified age-dependent speciation process may be that recently evolved taxa often colonize new regions or niches and may initially experience little competition. These new taxa are thus more likely to give rise to further new taxa than a taxon that has remained largely unchanged and is, therefore, well adapted to its niche. We show that age-dependent speciation may also be the result of different within-species populations following the same laws of lineage splitting to produce new species. As the fit of our model to the tree database shows, this simple biological motivation provides an explanation for a long standing problem in macroevolution.

Dispersal-diversity feedbacks and their consequences for macroecological patterns.

Dispersal is a key process in ecology and evolution. While the effects of dispersal on diversity are broadly acknowledged, our understanding of the influence of diversity on dispersal remains limited. This arises from the dynamic, context-dependent, nonlinear and ubiquitous nature of dispersal. Diversity outcomes, such as competition, mutualism, parasitism and trophic interactions can feed back on dispersal, thereby influencing biodiversity patterns at several spatio-temporal scales. Here, we shed light on the dispersal-diversity causal links by discussing how dispersal-diversity ecological and evolutionary feedbacks can impact macroecological patterns. We highlight the importance of dispersal-diversity feedbacks for advancing our understanding of macro-eco-evolutionary patterns and their challenges, such as establishing a unified framework for dispersal terminology and methodologies across various disciplines and scales. This article is part of the theme issue 'Diversity-dependence of dispersal: interspecific interactions determine spatial dynamics'.

Freshwater fish diversity in the western Amazon basin shaped by Andean uplift since the Late Cretaceous.

South America is home to the highest freshwater fish biodiversity on Earth, and the hotspot of species richness is located in the western Amazon basin. The location of this hotspot is enigmatic, as it is inconsistent with the pattern observed in river systems across the world of increasing species richness towards a river's mouth. Here we investigate the role of river capture events caused by Andean mountain building and repeated episodes of flooding in western Amazonia in shaping the modern-day richness pattern of freshwater fishes in South America, and in Amazonia in particular. To this end, we combine a reconstruction of river networks since 80 Ma with a mechanistic model simulating dispersal, allopatric speciation and extinction over the dynamic landscape of rivers and lakes. We show that Andean mountain building and consequent numerous small river capture events in western Amazonia caused freshwater habitats to be highly dynamic, leading to high diversification rates and exceptional richness. The history of marine incursions and lakes, including the Miocene Pebas mega-wetland system in western Amazonia, played a secondary role.

Climatic effects on niche evolution in a passerine bird clade depend on paleoclimate reconstruction method.

Climatic niches describe the climatic conditions in which species can persist. Shifts in climatic niches have been observed to coincide with major climatic change, suggesting that species adapt to new conditions. We test the relationship between rates of climatic niche evolution and paleoclimatic conditions through time for 65 Old-World flycatcher species (Aves: Muscicapidae). We combine niche quantification for all species with dated phylogenies to infer past changes in the rates of niche evolution for temperature and precipitation niches. Paleoclimatic conditions were inferred independently using two datasets: a paleoelevation reconstruction and the mammal fossil record. We find changes in climatic niches through time, but no or weak support for a relationship between niche evolution rates and rates of paleoclimatic change for both temperature and precipitation niche and for both reconstruction methods. In contrast, the inferred relationship between climatic conditions and niche evolution rates depends on paleoclimatic reconstruction method: rates of temperature niche evolution are significantly negatively related to absolute temperatures inferred using the paleoelevation model but not those reconstructed from the fossil record. We suggest that paleoclimatic change might be a weak driver of climatic niche evolution in birds and highlight the need for greater integration of different paleoclimate reconstructions.

Harnessing paleo-environmental modeling and genetic data to predict intraspecific genetic structure.

Spatially explicit simulations of gene flow within complex landscapes could help forecast the responses of populations to global and anthropological changes. Simulating how past climate change shaped intraspecific genetic variation can provide a validation of models in anticipation of their use to predict future changes. We review simulation models that provide inferences on population genetic structure. Existing simulation models generally integrate complex demographic and genetic processes but are less focused on the landscape dynamics. In contrast to previous approaches integrating detailed demographic and genetic processes and only secondarily landscape dynamics, we present a model based on parsimonious biological mechanisms combining habitat suitability and cellular processes, applicable to complex landscapes. The simulation model takes as input (a) the species dispersal capacities as the main biological parameter, (b) the species habitat suitability, and (c) the landscape structure, modulating dispersal. Our model emphasizes the role of landscape features and their temporal dynamics in generating genetic differentiation among populations within species. We illustrate our model on caribou/reindeer populations sampled across the entire species distribution range in the Northern Hemisphere. We show that simulations over the past 21 kyr predict a population genetic structure that matches empirical data. This approach looking at the impact of historical landscape dynamics on intraspecific structure can be used to forecast population structure under climate change scenarios and evaluate how species range shifts might induce erosion of genetic variation within species.

The Latitudinal Diversity Gradient: Novel Understanding through Mechanistic Eco-evolutionary Models.

The latitudinal diversity gradient (LDG) is one of the most widely studied patterns in ecology, yet no consensus has been reached about its underlying causes. We argue that the reasons for this are the verbal nature of existing hypotheses, the failure to mechanistically link interacting ecological and evolutionary processes to the LDG, and the fact that empirical patterns are often consistent with multiple explanations. To address this issue, we synthesize current LDG hypotheses, uncovering their eco-evolutionary mechanisms, hidden assumptions, and commonalities. Furthermore, we propose mechanistic eco-evolutionary modeling and an inferential approach that makes use of geographic, phylogenetic, and trait-based patterns to assess the relative importance of different processes for generating the LDG.

TreeSimGM: Simulating phylogenetic trees under general Bellman-Harris models with lineage-specific shifts of speciation and extinction in R.

Understanding macroevolutionary processes using phylogenetic trees is a challenging and complex process that draws on mathematics, computer science and biology. Given the development of complex mathematical models and the growing computational processing power, simulation tools are becoming increasingly popular.In order to simulate phylogenetic trees, most evolutionary biologists are forced to build their own algorithms or use existing tools built on different platforms and/or as standalone programmes. The absence of a simulation tool accommodating for user-chosen model specifications limits, amongst others, model testing and pipelining with approximate Bayesian computation methods or other subsequent statistical analysis.We introduce "TreeSimGM," an r-package simulation tool for phylogenetic trees under a general Bellman and Harris model. This package allows the user to specify any desired probability distribution for the waiting times until speciation and extinction (e.g. age-dependent speciation/extinction). Upon speciation, the user can specify whether one descendant species corresponds to the ancestor species inheriting its age or whether both descendant species are new species of age 0. Moreover, it is possible to scale the waiting time to speciation/extinction for newly formed species. Thus, "TreeSimGM" not only allows the user to simulate stochastic phylogenetic trees assuming several popular existing models, such as the Yule model, the constant-rate birth-death model, and proportional to distinguishable arrangement models, but it also allows the user to formulate new models for exploration. A short explanation of the supported models and a few examples of how to use our package are presented here.As an r-package, "TreeSimGM" allows flexible and powerful stochastic phylogenetic tree simulations. Moreover, it facilitates the pipelining of outputs or inputs with other functions in r. "TreeSimGM" contributes to the tools available to the r community in the fields of ecology and evolution, is freely available under the GPL-2 licence and can be downloaded at https://cran.r-project.org/web/packages/TreeSimGM.

Evaluating alternative explanations for an association of extinction risk and evolutionary uniqueness in multiple insular lineages.

Studies in insular environments have often documented a positive association of extinction risk and evolutionary uniqueness (i.e., how distant a species is from its closest living relative). However, the cause of this association is unclear. One explanation is that species threatened with extinction are evolutionarily unique because they are old, implying that extinction risk increases with time since speciation (age-dependent extinction). An alternative explanation is that such threatened species are last survivors of clades that have undergone an elevated extinction rate, and that their uniqueness results from the extinction of their close relatives. Distinguishing between these explanations is difficult but important, since they imply different biological processes determining extinction patterns. Here, we designed a simulation approach to distinguish between these alternatives using living species, and applied it to 12 insular radiations that show a positive association between extinction risk and evolutionary uniqueness. We also tested the sensitivity of results to underlying assumptions and variable extinction rates. Despite differences among the radiations considered, age-dependent extinction was supported as best explaining the majority of the empirical cases. Biological processes driving characteristic changes in abundance with species duration (age-dependency) may merit further investigation.

Fauna de besouros bioluminescentes (Coleoptera: Elateroidea: Lampyridae; Phengodidae, Elateridae) nos municípios de Campinas, Sorocaba-Votorantim e Rio Claro-Limeira (SP, Brasil): biodiversidade e influência da urbanização / Bioluminescent beetles (Coleoptera: Elateroidea: Lampyridae, Phengodidae, Elateridae) in the municipalities of Campinas, Sorocaba-Votorantim and Rio Claro-Limeira (SP, Brazil): biodiversity and influence of urban sprawl

O Brasil abriga a maior biodiversidade de besouros bioluminescentes do mundo. Apesar disto, levantamentos faunísticos regionais e estudos sobre o efeito da urbanização na biodiversidade de vaga-lumes são inexistentes. Este tipo de estudo pode ser especialmente importante para selecionar novos bioindicadores de ambientes noturnos. A biodiversidade de espécies bioluminescentes de Elateroidea em três grandes áreas urbanas no Estado de São Paulo: os municípios de Campinas, Sorocaba-Votorantim e Rio Claro-Limeira foi catalogada. Em Campinas, a ocorrência de espécies tem sido monitorada durante os últimos 20 anos. Foi registrada a ocorrência de 26 espécies em Campinas, 21 em Sorocaba e 19 em Rio Claro, sendo estas encontradas principalmente em áreas de Mata Atlântica, crescimentos secundários, brejos e campos. Apesar de vaga-lumes lampirídeos serem encontrados em todos os habitats, elaterídeos e fengodídeos ocorreram encontrados preferencialmente em áreas florestadas. A expansão urbana, em especial a iluminação artificial tem impacto evidente na biodiversidade deste grupo de insetos. Os lampirídeos Aspisoma lineatum Gyllendal, 1817 (Lampyrinae: Cratomorphini) e Bicellonychia lividipennis Motschulsky, 1854 (Photurinae) são as espécies mais frequentes em areas urbanizadas, ao passo que várias outras espécies de lampirídeos, elaterídeos e fengodídeos são as mais incomuns nestas áreas. Estes estudos oferecem subsídios para o uso deste grupo de vaga-lumes como potenciais bioindicadores noturnos.

Brazil hosts the richest biodiversity of bioluminescent beetles in the world. However, no surveys on the regional biodiversity have been made. Furthermore, studies on the effect of urbanization on bioluminescent beetle's biodiversity are still missing. Such kind of studies could be especially helpful to select new night environment bioindicators. Bioluminescent species of the Elateroidea superfamily occurring in three large and expanding urban areas in São Paulo State: Campinas, Sorocaba-Votorantim and Rio Claro-Limeira municipalities, were catalogued. In Campinas, the occurrence of bioluminescent beetles in different sites has been monitored during the past 20 years. There were 26 species in Campinas, 21 in Sorocaba-Votorantim and 19 in Rio Claro-Limeira, occurring mainly in conserved remnants of Atlantic Rain-forest, secondary growths, marshes and open fields. Whereas fireflies (Lampyridae) are found in almost all habitats, click beetles (Elateridae) and railroadworms (Phengodidae) occur predominantly in woody environments. Urban sprawl has an evident impact on the biodiversity of this selected group of beetles, and special attention should be given to the effect of artificial night lighting in the occurrence of these beetles. The lampyrids Aspisoma lineatum Gyllenhal, 1817 (Lampyrinae: Cratomorphini) and Bicellonychia lividipennis Motschulsky 1854 (Photurinae) are the most common species in urbanized areas, whereas many other species of lampyrids, elaterids and phengodids are rarely found in such environments. These studies offer insights on the potential use of a selected group of luminescent beetles found in this region as nocturnal environment bioindicators.