Branching patterns in phylogenies cannot distinguish diversity-dependent diversification from time-dependent diversification

Analysis of intrinsic evolutionary factors leading to microendemic distributions in New Caledonian leaf beetles

Microendemicity, or the condition of some species having local ranges, is a relatively common pattern in nature. However, the factors that lead to this pattern are still largely unknown. Most studies addressing this issue tend to focus on extrinsic factors associated with microendemic distributions, such as environmental conditions, hypothesising a posteriori about underlying potential speciation mechanisms, linked or not to these conditions. Here, we use a multi-faceted approach mostly focusing on intrinsic factors instead, namely diversification dynamics and speciation modes in two endemic sibling genera of leaf beetles with microendemic distributions, Taophila and Tricholapita, in a microendemicity hotspot, New Caledonia. Results suggest that the diversification rate in this lineage slowed down through most of the Neogene and consistently with a protracted speciation model possibly combined with several ecological and environmental factors potentially adding rate-slowing effects through time. In turn, species accumulated following successive allopatric speciation cycles, possibly powered by marked geological and climatic changes in the region in the last 25 million years, with daughter species ranges uncorrelated with the time of speciation. In this case, microendemicity seems to reflect a mature state for the system, rather than a temporary condition for recent species, as suggested for many microendemic organisms.

Slow and steady wins the race: Diversification rate is independent from body size and lifestyle in Malagasy skinks (Squamata: Scincidae: Scincinae)

Most of the unique and diverse vertebrate fauna that inhabits Madagascar derives from in situ diversification from colonisers that reached this continental island through overseas dispersal. The endemic Malagasy Scincinae lizards are amongst the most species-rich squamate groups on the island. They colonised all bioclimatic zones and display many ecomorphological adaptations to a fossorial (burrowing) lifestyle. Here we propose a new phylogenetic hypothesis for their diversification based on the largest taxon sampling so far compiled for this group. We estimated divergence times and investigated several aspects of their diversification (diversification rate, body size and fossorial lifestyle evolution, and biogeography). We found that diversification rate was constant throughout most of the evolutionary history of the group, but decreased over the last 6-4 million years and independently from body size and fossorial lifestyle evolution. Fossoriality has evolved from fully quadrupedal ancestors at least five times independently, which demonstrates that even complex morphological syndromes - in this case involving traits such as limb regression, body elongation, modification of cephalic scalation, depigmentation, and eyes and ear-opening regression - can evolve repeatedly and independently given enough time and eco-evolutionary advantages. Initial diversification of the group likely occurred in forests, and the divergence of sand-swimmer genera around 20 Ma appears linked to a period of aridification. Our results show that the large phenotypic variability of Malagasy Scincinae has not influenced diversification rate and that their rich species diversity results from a constant accumulation of lineages through time. By compiling large geographic and trait-related datasets together with the computation of a new time tree for the group, our study contributes important insights on the diversification of Malagasy vertebrates.

Evolutionary models demonstrate rapid and adaptive diversification of Australo-Papuan pythons

Lineages may diversify when they encounter available ecological niches. Adaptive divergence by ecological opportunity often appears to follow the invasion of a new environment with open ecological space. This evolutionary process is hypothesized to explain the explosive diversification of numerous Australian vertebrate groups following the collision of the Eurasian and Australian plates 25 Mya. One of these groups is the pythons, which demonstrate their greatest phenotypic and ecological diversity in Australo-Papua (Australia and New Guinea). Here, using an updated and near complete time-calibrated phylogenomic hypothesis of the group, we show that following invasion of this region, pythons experienced a sudden burst of speciation rates coupled with multiple instances of accelerated phenotypic evolution in head and body shape and body size. These results are consistent with adaptive radiation theory with an initial rapid niche-filling phase and later slow-down approaching niche saturation. We discuss these findings in the context of other Australo-Papuan adaptive radiations and the importance of incorporating adaptive diversification systems that are not extraordinarily species-rich but ecomorphologically diverse to understand how biodiversity is generated.

Is species richness mediated by functional and genetic divergence? A global analysis in birds

Unravelling why species richness shows such dramatic spatial variation is an ongoing challenge. Common to many theories is that increasing species richness (e.g. with latitude) requires a compensatory trade-off on an axis of species' ecology. Spatial variation in species richness may also affect genetic diversity if large numbers of coexisting, related species result in smaller population sizes.Here, we test whether increasing species richness results in differential occupation of morphospace by the constituent species, or decreases species' genetic diversity. We test for two potential mechanisms of morphological accommodation: denser packing in ecomorphological space, and expansion of the space. We then test whether species differ in their nucleotide diversity depending on allopatry or sympatry with relatives, indicative of potential genetic consequences of coexistence that would reduce genetic diversity in sympatry. We ask these questions in a spatially explicit framework, using a global database of avian functional trait measurements in combination with >120,000 sequences downloaded from GenBank.We find that higher species richness within families is not systematically correlated with either packing in morphological space or overdispersion but, at the Class level, we find a general positive relationship between packing and species richness, but that points sampled in the tropics have comparatively greater packing than temperate ones relative to their species richness. We find limited evidence that geographical co-occurrence with closely related species or tropical distributions decreases nucleotide diversity of nuclear genes; however, this requires further analysis.Our results suggest that avian families can accumulate species regionally with minimal tradeoffs or cost, implying that external biotic factors do not limit species richness. Read the free Plain Language Summary for this article on the Journal blog.

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.

Causes and Consequences of Apparent Timescaling Across All Estimated Evolutionary Rates

Evolutionary rates play a central role in connecting micro- and macroevolution. All evolutionary rate estimates, including rates of molecular evolution, trait evolution, and lineage diversification, share a similar scaling pattern with time: The highest rates are those measured over the shortest time interval. This creates a disconnect between micro- and macroevolution, although the pattern is the opposite of what some might expect: Patterns of change over short timescales predict that evolution has tremendous potential to create variation and that potential is barely tapped by macroevolution. In this review, we discuss this shared scaling pattern across evolutionary rates. We break down possible explanations for scaling into two categories, estimation error and model misspecification, and discuss how both apply to each type of rate. We also discuss the consequences of this ubiquitous pattern, which can lead to unexpected results when comparing ratesover different timescales. Finally, after addressing purely statistical concerns, we explore a few possibilities for a shared unifying explanation across the three types of rates that results from a failure to fully understand and account for how biological processes scale over time.

Why extinction estimates from extant phylogenies are so often zero

Time-calibrated phylogenies of extant species ("extant timetrees") are widely used to estimate historical speciation and extinction rates by fitting stochastic birth-death models.¹ These approaches have long been controversial, as many phylogenetic studies report zero extinction in many taxa, contradicting the high extinction rates seen in the fossil record and the fact that the majority of species ever to have existed are now extinct.^2-9 To date, the causes of this discrepancy remain unresolved. Here, we provide a novel and simple explanation for these "zero-inflated" extinction estimates, based on the recent discovery that there exist many alternative "congruent" diversification scenarios that cannot be distinguished based solely on extant timetrees.¹⁰ Due to such congruencies, estimation methods tend to converge to some scenario congruent to (i.e., statistically indistinguishable from) the true diversification scenario, but not necessarily to the true diversification scenario itself. This congruent scenario may exhibit negative extinction rates, a biologically meaningless but mathematically feasible situation, in which case estimators will tend to stick to the boundary of zero extinction. Based on this explanation, we make multiple testable predictions, which we confirm using analyses of simulated trees and 121 empirical trees. In contrast to other proposed mechanisms for erroneous extinction rate estimates,^5,11-14 our proposed mechanism specifically explains the zero inflation of previous extinction rate estimates in the absence of detectable model violations, even for large trees. Not only do our results likely resolve a long-standing mystery in phylogenetics, they demonstrate that model congruencies can have severe consequences in practice.

Empirical and Methodological Challenges to the Model-Based Inference of Diversification Rates in Extinct Clades

Changes in speciation and extinction rates are key to the dynamics of clade diversification, but attempts to infer them from phylogenies of extant species face challenges. Methods capable of synthesizing information from extant and fossil species have yielded novel insights into diversification rate variation through time, but little is known about their behavior when analyzing entirely extinct clades. Here, we use empirical and simulated data to assess how two popular methods, PyRate and Fossil BAMM, perform in this setting. We inferred the first tip-dated trees for ornithischian dinosaurs, and combined them with fossil occurrence data to test whether the clade underwent an end-Cretaceous decline. We then simulated phylogenies and fossil records under empirical constraints to determine whether macroevolutionary and preservation rates can be teased apart under paleobiologically realistic conditions. We obtained discordant inferences about ornithischian macroevolution including a long-term speciation rate decline (BAMM), mostly flat rates with a steep diversification drop (PyRate) or without one (BAMM), and episodes of implausibly accelerated speciation and extinction (PyRate). Simulations revealed little to no conflation between speciation and preservation, but yielded spuriously correlated speciation and extinction estimates while time-smearing tree-wide shifts (BAMM) or overestimating their number (PyRate). Our results indicate that the small phylogenetic datasets available to vertebrate paleontologists and the assumptions made by current model-based methods combine to yield potentially unreliable inferences about the diversification of extinct clades. We provide guidelines for interpreting the results of the existing approaches in light of their limitations, and suggest how the latter may be mitigated.