Complex patterns of gene flow and convergence in the evolutionary history of the spiral-horned antelopes (Tragelaphini)

The Tragelaphini, also known as spiral-horned antelope, is a phenotypically diverse mammalian tribe comprising a single genus, Tragelaphus. The evolutionary history of this tribe has attracted the attention of taxonomists and molecular geneticists for decades because its diversity is characterised by conflicts between morphological and molecular data as well as between mitochondrial, nuclear and chromosomal DNA. These inconsistencies point to a complex history of ecological diversification, coupled by either phenotypic convergence or introgression. Therefore, to unravel the phylogenetic relationships among spiral-horned antelopes, and to further investigate the role of divergence and gene flow in trait evolution, we sequenced genomes for all nine accepted species of the genus Tragelaphus, including a genome each for the highly divergent bushbuck lineages (T. s. scriptus and T. s. sylvaticus). We successfully reconstructed the Tragelaphus species tree, providing genome-level support for the early Pliocene divergence and monophyly of the nyala (T. angasii) and lesser kudu (T. imberbis), the monophyly of the two eland species (T. oryx and T. derbianus) and, importantly, the monophyly of kéwel (T. s. scriptus) and imbabala (T. s. sylvaticus) bushbuck. We found strong evidence for gene flow in at least four of eight nodes on the species tree. Among the six phenotypic traits assessed here, only habitat type mapped onto the species tree without homoplasy, showing that trait evolution was the result of complex patterns of divergence, introgression and convergent evolution.

Deep Learning from Phylogenies for Diversification Analyses

Birth-death (BD) models are widely used in combination with species phylogenies to study past diversification dynamics. Current inference approaches typically rely on likelihood-based methods. These methods are not generalizable, as a new likelihood formula must be established each time a new model is proposed; for some models, such a formula is not even tractable. Deep learning can bring solutions in such situations, as deep neural networks can be trained to learn the relation between simulations and parameter values as a regression problem. In this paper, we adapt a recently developed deep learning method from pathogen phylodynamics to the case of diversification inference, and we extend its applicability to the case of the inference of state-dependent diversification models from phylogenies associated with trait data. We demonstrate the accuracy and time efficiency of the approach for the time-constant homogeneous BD model and the Binary-State Speciation and Extinction model. Finally, we illustrate the use of the proposed inference machinery by reanalyzing a phylogeny of primates and their associated ecological role as seed dispersers. Deep learning inference provides at least the same accuracy as likelihood-based inference while being faster by several orders of magnitude, offering a promising new inference approach for the deployment of future models in the field.

Interactions between fitness components across the life cycle constrain competitor coexistence

Numerous mechanisms can promote competitor coexistence. Yet, these mechanisms are often considered in isolation from one another. Consequently, whether multiple mechanisms shaping coexistence combine to promote or constrain species coexistence remains an open question. Here, we aim to understand how multiple mechanisms interact within and between life stages to determine frequency-dependent population growth, which has a key role stabilizing local competitor coexistence. We conducted field experiments in three lakes manipulating relative frequencies of two Enallagma damselfly species to evaluate demographic contributions of three mechanisms affecting different fitness components across the life cycle: the effect of resource competition on individual growth rate, predation shaping mortality rates, and mating harassment determining fecundity. We then used a demographic model that incorporates carry-over effects between life stages to decompose the relative effect of each fitness component generating frequency-dependent population growth. This decomposition showed that fitness components combined to increase population growth rates for one species when rare, but they combined to decrease population growth rates for the other species when rare, leading to predicted exclusion in most lakes. Because interactions between fitness components within and between life stages vary among populations, these results show that local coexistence is population specific. Moreover, we show that multiple mechanisms do not necessarily increase competitor coexistence, as they can also combine to yield exclusion. Identifying coexistence mechanisms in other systems will require greater focus on determining contributions of different fitness components across the life cycle shaping competitor coexistence in a way that captures the potential for population-level variation.

The Major Features of Macroevolution

Evolutionary dynamics operating across deep time leave footprints in the shapes of phylogenetic trees. For the last several decades, researchers have used increasingly large and robust phylogenies to study the evolutionary history of individual clades and to investigate the causes of the glaring disparities in diversity among groups. Whereas typically not the focal point of individual clade-level studies, many researchers have remarked on recurrent patterns that have been observed across many different groups and at many different time scales. Whereas previous studies have documented various such regularities in topology and branch length distributions, they have typically focused on a single pattern and used a disparate collection (oftentimes, of quite variable reliability) of trees to assess it. Here we take advantage of modern megaphylogenies and unify previous disparate observations about the shapes embedded in the Tree of Life to create a catalog of the "major features of macroevolution." By characterizing such a large swath of subtrees in a consistent way, we hope to provide a set of phenomena that process-based macroevolutionary models of diversification ought to seek to explain.

Statistical methods to identify mechanisms in studies of eco-evolutionary dynamics

While the reciprocal effects of ecological and evolutionary dynamics are increasingly recognized as an important driver for biodiversity, detection of such eco-evolutionary feedbacks, their underlying mechanisms, and their consequences remains challenging. Eco-evolutionary dynamics occur at different spatial and temporal scales and can leave signatures at different levels of organization (e.g., gene, protein, trait, community) that are often difficult to detect. Recent advances in statistical methods combined with alternative hypothesis testing provides a promising approach to identify potential eco-evolutionary drivers for observed data even in non-model systems that are not amenable to experimental manipulation. We discuss recent advances in eco-evolutionary modeling and statistical methods and discuss challenges for fitting mechanistic models to eco-evolutionary data.

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.

Diversification dynamics in Caribbean rain frogs (Eleutherodactylus) are uncoupled from the anuran community and consistent with adaptive radiation

Adaptive radiation is proposed to play a key role in generating differences in species richness among lineages and geographical regions. Due to the importance of ecological divergence in adaptive radiation, species richness is predicted to be influenced by equilibrium diversity dynamics, although the concept continues to generate much debate. An additional important question is whether radiating clades have intrinsic biological characteristics that make them particularly prone to diversify. We tackle these questions by analysing (i) the temporal patterns of diversification of Caribbean Eleutherodactylus frogs, and (ii) assembly of the complete native anuran community of the Caribbean archipelago (197 species), testing for the presence of equilibrium dynamics and whether diversification patterns of Eleutherodactylus differ from those of the rest of the Caribbean anurans. Diversification rates follow the predicted pattern of rapid diversification early in the radiation which gradually decreases towards the present. Eleutherodactylus diversification is significantly faster than that of the Caribbean anuran community, and although equilibrium dynamics influence richness of all Caribbean anurans, Eleutherodactylus shows higher carrying capacity. Our results indicate that ecological opportunity per se is not sufficient for adaptive radiation and that diverse lineages present intrinsic characteristics that enable them to make the most of available opportunity.

Illusion of flight? Absence, evidence and the age of winged insects

The earliest fossils of winged insects (Pterygota) are mid-Carboniferous (latest Mississippian, 328–324 Mya), but estimates of their age based on fossil-calibrated molecular phylogenetic studies place their origin at 440–370 Mya during the Silurian or Devonian. This discrepancy would require that winged insects evaded fossilization for at least the first ~50 Myr of their history. Here, we examine the plausibility of such a gap in the fossil record, and possible explanations for it, based on comparisons with the fossil records of other arthropod groups, the distribution of first occurrence dates of pterygote families, phylogenetically informed simulations of the fossilization of Palaeozoic insects, and re-analysis of data presented by Misof and colleagues using updated fossil calibrations under a variety of prior probability settings. We do not find support for the mechanisms previously suggested to account for such an extended gap in the pterygote fossil record, including sampling bias, preservation bias, and body size. We suggest that inference of an early origin of Pterygota long prior to their first appearance in the fossil record is probably an analytical artefact of taxon sampling and choice of fossil calibration points, possibly compounded by heterogeneity in rates of sequence evolution or speciation, including radiations or ‘bursts’ during their early history.

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.

Cryptic diversity impacts model selection and macroevolutionary inferences in diversification analyses

Species persist in landscapes through ecological dynamics but proliferate at wider spatial scales through evolutionary mechanisms. Disentangling the contribution of each dynamic is challenging, but the increasing use of dated molecular phylogenies opened new perspectives. First, the increasing use of DNA sequences in biodiversity inventory shed light on a substantial amount of cryptic diversity in species-rich ecosystems. Second, explicit diversification models accounting for various eco-evolutionary models are now available. Integrating both advances, we explored diversification trajectories among 10 lineages of freshwater fishes in Sundaland, for which time-calibrated and taxonomically rich phylogenies are available. By fitting diversification models to dated phylogenies and incorporating DNA-based species delimitation methods, the impact of cryptic diversity on diversification model selection and related inferences is explored. Eight clades display constant speciation rate model as the most likely if cryptic diversity is accounted, but nine display a signature of diversification slowdowns when cryptic diversity is ignored. Cryptic diversification occurs during the last 5 Myr for most groups, and palaeoecological models received little support. Most cryptic lineages display restricted range distribution, supporting geographical isolation across homogeneous landscapes as the main driver of diversification. These patterns question the persistence of cryptic diversity and its role during species proliferation.

Mapping species diversification metrics in macroecology: Prospects and challenges

The intersection of macroecology and macroevolution is one of today’s most active research in biology. In the last decade, we have witnessed a steady increment of macroecological studies that use metrics attempting to capture macroevolutionary processes to explain present-day biodiversity patterns. Evolutionary explanations of current species richness gradients are fundamental for understanding how diversity accumulates in a region. Although multiple hypotheses have been proposed to explain the patterns we observe in nature, it is well-known that the present-day diversity patterns result from speciation, extinction, colonization from nearby areas, or a combination of these macroevolutionary processes. Whether these metrics capture macroevolutionary processes across space is unknown. Some tip-rate metrics calculated directly from a phylogenetic tree (e.g., mean root distance -MRD-; mean diversification rate -mDR-) seem to return very similar geographical patterns regardless of how they are estimated (e.g., using branch lengths explicitly or not). Model-based tip-rate metrics —those estimated using macroevolutionary mixtures, e.g., the BAMM approach— seem to provide better net diversification estimates than only speciation rates. We argue that the lack of appropriate estimates of extinction and dispersal rates in phylogenetic trees may strongly limit our inferences about how species richness gradients have emerged at spatial and temporal scales. Here, we present a literature review about this topic and empirical comparisons between select taxa with several of these metrics. We implemented a simple null model approach to evaluate whether mapping of these metrics deviates from a random sampling process. We show that phylogenetic metrics by themselves are relatively poor at capturing speciation, extinction, and dispersal processes across geographical gradients. Furthermore, we provide evidence of how parametric biogeographic methods can improve our inference of past events and, therefore, our conclusions about the evolutionary processes driving biodiversity patterns. We recommend that further studies include several approaches simultaneously (e.g., spatial diversification modeling, parametric biogeographic methods, simulations) to disentangle the relative role of speciation, extinction, and dispersal in the generation and maintenance of species richness gradients at regional and global scales.

Temperature-dependent evolutionary speed shapes the evolution of biodiversity patterns across tetrapod radiations

Biodiversity varies predictably with environmental energy around the globe, but the underlaying mechanisms remain incompletely understood. The evolutionary speed hypothesis predicts that environmental kinetic energy shapes variation in speciation rates through temperature- or life history-dependent rates of evolution. To test whether variation in evolutionary speed can explain the relationship between energy and biodiversity in birds, mammals, amphibians, and reptiles, we simulated diversification over 65 million years of geological and climatic change with a spatially explicit eco-evolutionary simulation model. We modelled four distinct evolutionary scenarios in which speciation-completion rates were dependent on temperature (M1), life history (M2), temperature and life history (M3), or were independent of temperature and life-history (M0). To assess the agreement between simulated and empirical data, we performed model selection by fitting supervised machine learning models to multidimensional biodiversity patterns. We show that a model with temperature-dependent rates of speciation (M1) consistently had the strongest support. In contrast to statistical inferences, which showed no general relationships between temperature and speciation rates in tetrapods, we demonstrate how process-based modelling can disentangle the causes behind empirical biodiversity patterns. Our study highlights how environmental energy has played a fundamental role in the evolution of biodiversity over deep time.

The scaling of diversification rates with age is likely explained by sampling bias

Numerous phylogenetic studies reported the existence of a pervasive scaling relationship between the ages of extant eukaryotic clades and their estimated diversification rates. The causes of this age-rate-scaling (ARS), whether biological and/or artifactual, remain unresolved. Here we fit diversification models to thousands of eukaryotic time-calibrated phylogenies to explore multiple potential causes of the ARS including parameter non-identifiability, model inadequacy, biases in taxonomic practice, and an important and ubiquitous form of sampling bias-preferentially analyzing larger extant clades. We distinguish between two mechanism by which such sampling biases can cause an ARS: First, by favoring clades that happen to be unusually large merely by chance (i.e., due to the stochastic nature of the cladogenic process), thus leading to rate overestimation, and second, by favoring clades that have truly higher diversification rates. We find that, of the proposed explanations, only sampling biases are likely to contribute to the observed ARS. We develop methods for fully correcting for sampling bias mechanism 1, and find that despite these corrections a substantial ARS remains. We then confirm using simulations that preferring trees with truly higher rates (mechanism 2) likely explains this residual ARS. Since we do not have a completely unbiased sample of clades, including extinct ones, for phylogenetic analyses, it is difficult to demonstrate unambiguously that sampling biases are the sole cause of the ARS. Sampling biases are, however, a parsimonious and plausible explanation for this widely observed macroevolutionary pattern, and this has implications for how we interpret the distribution of diversification rate estimates in extant clades.

The enrichment paradox in adaptive radiations: Emergence of predators hinders diversification in resource rich environments

Adaptive radiations are known for rapid niche diversification in response to ecological opportunity. While most resources usually exist prior to adaptive radiation, novel niches associated with novel resources can be created as a clade diversifies. For example, in African lake cichlid radiations some species prey upon other species of the clade (intraclade consumers). Using a trait-based eco-evolutionary model, we investigate the evolution of intraclade consumers in adaptive radiations and the effect of this novel trophic interaction on the diversification process of the radiating clade. We find that the evolutionary emergence of intraclade consumers halts the diversification processes of other ecomorphs as a result of increased top-down control of density. Because high productivity enables earlier evolution of intraclade consumers, highly productive environments come to harbour less species-rich radiations than comparable radiations in less productive environments. Our results reveal how macroevolutionary and community patterns can emerge from ecological and microevolutionary processes.

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.

Multiple paths to morphological diversification during the origin of amniotes

Early terrestrial vertebrates (amniotes) provide a classic example of diversification following adaptive zone invasion. The initial terrestrialization of vertebrates was closely followed by dietary diversification, as evidenced by a proliferation of craniomandibular and dental adaptations. However, morphological evolution of early amniotes has received limited study, in analyses with restricted taxonomic scope, leaving substantial questions about the dynamics of this important terrestrial radiation. We use novel analyses of discrete characters to quantify variation in evolutionary rates and constraints during diversification of the amniote feeding apparatus. We find evidence for an early burst, comprising high rates of anatomical change that decelerated through time, giving way to a background of saturated morphological evolution. Subsequent expansions of phenotypic diversity were not associated with increased evolutionary rates. Instead, variation in the mode of evolution became important, with groups representing independent origins of herbivory evolving distinctive, group-specific morphologies and thereby exploring novel character-state spaces. Our findings indicate the importance of plant-animal interactions in structuring the earliest radiation of amniotes and demonstrate the importance of variation in modes of phenotypic divergence during a major evolutionary radiation.

Modelling the tempo and mode of lineage dispersal

Lineage dispersal is a basic macroevolutionary process shaping the distribution of biodiversity. Probabilistic approaches in biogeography, epidemiology, and macroecology often model dispersal as a background process to explain extant or infer past distributions. We propose framing questions around the mode, timing, rate, and direction of lineage dispersal itself, from a lineage- or geography-centric perspective. We review available methods for modelling lineage dispersal. Likelihood- and simulation-based approaches to modelling dispersal have made progress in accounting for the variation of lineage dispersal over space, time, and branches of a phylogeny and its interaction with diversification. Methodological improvements, guided by a focus on model adequacy, will lead to more realistic models that can answer fundamental questions about the tempo and mode of lineage dispersal.

Molecular assessment of dietary niche partitioning in an endemic island radiation of tropical mammals

Island radiations represent unique evolutionary histories in unique ecological contexts. These radiations provide opportunities to investigate ecological diversification in groups that typically exhibit niche partitioning among their constituents, including partitioning of food resources. DNA metabarcoding produces finer levels of diet identification than traditional methods, allowing us to examine dietary niche partitioning in communities or clades in which species share superficially similar diets. Here, we use DNA metabarcoding to investigate dietary niche partitioning in an endemic radiation of mammals in the Philippines. Our data reveal niche partitioning as well as phylogenetically-uncorrelated adaptive evolution in this small mammal community. Because 70% of the focal species belong to the tribe Chrotomyini, an endemic Philippine radiation of murid rodents that feed extensively on earthworms, this study sheds light on dietary adaptation and its role in the co-occurrence of closely related species. Our results reveal fine-scale resource partitioning within this community; our data provide compelling evidence for niche partitioning of diet that was masked by previous diet categories and will help in further dissecting the model adaptive radiation of endemic small mammals on Luzon. This study reinforces the notion that DNA metabarcoding can be a valuable tool for investigating both ecological relationships and evolutionary ecology at the community and phylogenetic level, respectively.

Evolutionary Rate Variation Among Lineages in Gene Trees has a Negative Impact on Species-Tree Inference

Phylogenetic analyses of genomic data provide a powerful means of reconstructing the evolutionary relationships among organisms, yet such analyses are often hindered by conflicting phylogenetic signals among loci. Identifying the signals that are most influential to species-tree estimation can help to inform the choice of data for phylogenomic analysis. We investigated this in an analysis of 30 phylogenomic data sets. For each data set, we examined the association between several branch-length characteristics of gene trees and the distance between these gene trees and the corresponding species trees. We found that the distance of each gene tree to the species tree inferred from the full data set was positively associated with variation in root-to-tip distances and negatively associated with mean branch support. However, no such associations were found for gene-tree length, a measure of the overall substitution rate at each locus. We further explored the usefulness of the best-performing branch-based characteristics for selecting loci for phylogenomic analyses. We found that loci that yield gene trees with high variation in root-to-tip distances have a disproportionately distant signal of tree topology compared with the complete data sets. These results suggest that rate variation across lineages should be taken into consideration when exploring and even selecting loci for phylogenomic analysis.[Branch support; data filtering; nucleotide substitution model; phylogenomics; substitution rate; summary coalescent methods.]

Extraordinary diversification of the “bristle flies” (Diptera: Tachinidae) and its underlying causes

The family Tachinidae (“bristle flies”) is the most diverse and ecologically important group of insect parasitoids outside the parasitic wasps. It is among the most species rich families of flies (Diptera) and has experienced a recent adaptive radiation across the globe. We make use of a molecular phylogeny of the family to examine its rapid radiation and explore the traits of tachinid lineages that may have contributed to variation in their diversification. We apply a range of diversification analyses to assess the consistency and robustness of effects. We find that the Tachinidae are among the most rapidly diversifying families of animals. Six to eight clades of bristle flies, distributed across the phylogeny, exhibit strong evidence of accelerated diversification. Our results suggest that the use of holometabolous insect larvae, and specifically caterpillars (Lepidoptera), as hosts, is associated with increased diversification rates. However, these effects were inconsistent across analyses. We detected little influence of oviposition strategy (egg type) or host feeding habit, and we recovered evidence that unmeasured “hidden” traits may explain greater variance in diversification. We evaluated the strengths and weaknesses of different Maximum Likelihood and Bayesian approaches for analysing diversification and the potential for extrinsic factors, such as geography, to influence patterns of richness and diversification. In general, we conclude that although certain traits may provide opportunities for diversification, whether this is capitalized on may depend on additional traits and/or historical contingency.

Impact of Pleistocene Eustatic Fluctuations on Evolutionary Dynamics in Southeast Asian Biodiversity Hotspots

Pleistocene Climatic Fluctuations (PCF) are frequently highlighted as important evolutionary engines that triggered cycles of biome expansion and contraction. While there is ample evidence of the impact of PCF on biodiversity of continental biomes, the consequences in insular systems depend on the geology of the islands and the ecology of the taxa inhabiting them. The idiosyncratic aspects of insular systems are exemplified by the islands of the Sunda Shelf in Southeast Asia (Sundaland), where PCF-induced eustatic fluctuations had complex interactions with the geology of the region, resulting in high species diversity and endemism. Emergent land in Southeast Asia varied drastically with sea level fluctuations during the Pleistocene. Climate-induced fluctuations in sea level caused temporary connections between insular and continental biodiversity hotspots in Southeast Asia. These exposed lands likely had freshwater drainage systems that extended between modern islands: the Paleoriver Hypothesis. Built upon the assumption that aquatic organisms are among the most suitable models to trace ancient river boundaries and fluctuations of landmass coverage, the present study aims to examine the evolutionary consequences of PCF on the dispersal of freshwater biodiversity in Southeast Asia. Time-calibrated phylogenies of DNA-delimited species were inferred for six species-rich freshwater fish genera in Southeast Asia (Clarias, Channa, Glyptothorax, Hemirhamphodon, Dermogenys, Nomorhamphus). The results highlight rampant cryptic diversity and the temporal localization of most speciation events during the Pleistocene, with 88% of speciation events occurring during this period. Diversification analyses indicate that sea level-dependent diversification models poorly account for species proliferation patterns for all clades excepting Channa. Ancestral area estimations point to Borneo as the most likely origin for most lineages, with two waves of dispersal to Sumatra and Java during the last 5 Myr. Speciation events are more frequently associated with boundaries of the paleoriver watersheds, with 60%, than islands boundaries, with 40%. In total, one-third of speciation events are inferred to have occured within paleorivers on a single island, suggesting that habitat heterogeneity and factors other than allopatry between islands substantially affected diversification of Sundaland fishes. Our results suggest that species proliferation in Sundaland is not wholly reliant on Pleistocene sea-level fluctuations isolating populations on different islands.

Habitat preference and diversification rates in a speciose lineage of diving beetles

The long-term geological stability of aquatic habitats has been demonstrated to be a determinant in the evolution of macroinvertebrate fauna, with species in running (lotic) waters having lower dispersal abilities, smaller ranges and higher gene flow between populations than species in standing (lentic) environments. Lotic species have been hypothesized to be more specialised, but the diversification dynamics of both habitat types have not been studied in detail. Using a speciose lineage of water beetles we test here whether diversification rates are related to the habitat preference of the species and its consequences on turnover, which we expect to be higher for lotic taxa. Moreover, we tested whether life in lotic environments is acting as an evolutionary dead-end as it is considered an ecological specialisation. We built a comprehensive molecular phylogeny with 473 terminals representing 421 of the 689 known species of the tribe Hydroporini (Coleoptera, Dytiscidae), using a combination of sequences from four mitochondrial and two nuclear genes plus 69 mitogenomes obtained with NGS. We found a general pattern of gradual acceleration of diversification rate with time, with 2-3 significant diversification shifts. However, habitat is not the main factor driving diversification in Hydroporini based on SecSSE analyses. The most recent common ancestor of Hydroporini was reconstructed as a lotic species, with multiple shifts to lentic environments. Most frequent transitions were estimated from lentic and lotic habitats to the category "both", followed by transitions from lotic to lentic and lentic to lotic respectively, although with very similar rates. Contrary to expectations, we found little evidence for differences in diversification dynamics between habitats, with lotic environments clearly not acting as evolutionary dead-ends in Hydroporini.

The Dynamics of Bird Diversity in the New World

Three prominent explanations have been proposed to explain the dramatic differences in species richness across regions and elevations, (i) time for speciation, (ii) diversification rates, and (iii) ecological limits. But the relative importance of these explanations and, especially, their interplay and possible synthesis remain largely elusive. Integrating diversification analyses, null models, and geographic information systems, I study avian richness across regions and elevations of the New World. My results reveal that even though the three explanations are differentially important (with ecological limits playing the dominant role), each contributes uniquely to the formation of richness gradients. Further, my results reveal the likely interplay between the explanations. They indicate that ecological limits hinder the diversification process, such that the accumulation of species within a region gradually slows down over time. Yet, it does not seem to converge toward a hard ceiling on regional richness. Instead, species-rich regions show suppressed, but continued, diversification, coupled with signatures of possible competition (esp. Neotropical lowlands). Conversely, species-poor, newly-colonized regions show fast diversification and weak to no signs of competition (esp. Nearctic highlands). These results held across five families of birds, across grid cells, biomes, and elevations. Together, my findings begin to illuminate the rich, yet highly consistent, interplay of the mechanisms that together shape richness gradients in the New World, including the most species-rich biodiversity hotspots on the planet, the Andes and the Amazon. [Biogeography; community; competition; macroevolution; phylogenetics; richness gradient.]

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

One of the primary goals of macroevolutionary biology has been to explain general trends in long‐term diversity patterns, including whether such patterns correspond to an upscaling of processes occurring at lower scales. Reconstructed phylogenies often show decelerated lineage accumulation over time. This pattern has often been interpreted as the result of diversity‐dependent (DD) diversification, where the accumulation of species causes diversification to decrease through niche filling. However, other processes can also produce such a slowdown, including time dependence without diversity dependence. To test whether phylogenetic branching patterns can be used to distinguish these two mechanisms, we formulated a time‐dependent, but diversity‐independent model that matches the expected diversity through time of a DD model. We simulated phylogenies under each model and studied how well likelihood methods could recover the true diversification mode. Standard model selection criteria always recovered diversity dependence, even when it was not present. We correct for this bias by using a bootstrap method and find that neither model is decisively supported. This implies that the branching pattern of reconstructed trees contains insufficient information to detect the presence or absence of diversity dependence. We advocate that tests encompassing additional data, for example, traits or range distributions, are needed to evaluate how diversity drives macroevolutionary trends.

Skeletal marine animal biodiversity is built by families with long macroevolutionary lag times

The clade dynamics of marine animals have changed markedly over the Phanerozoic. Long-term diversification is associated with decreasing origination and extinction rates, and with increasing taxon longevity. Here we use the diversification trajectories of skeletal non-colonial marine families to infer the mechanisms that generated these trends. Suggested mechanisms behind these trends include stochastic extinction of taxa with high evolutionary volatility and selection for traits that buffer against extinction. We find an increasing predominance of Phanerozoic families with long lag times between first appearance and peak diversity, over those with 'early burst' diversification trajectories. Long-lag families persisted for longer and had lower evolutionary volatilities, higher genus-level occupancies and genera with larger niche breadths than early burst families. However, they do not preferentially show ecological modes known to protect against extinction. We interpret the rise of the long-lag families as reflecting an intensification of ecosystem-level mechanisms supporting both long-term coexistence and transient dynamics, which increased the capacity of marine ecosystems to accommodate highly diverse communities.

Diversification in evolutionary arenas-Assessment and synthesis

Understanding how and why rates of evolutionary diversification vary is a key issue in evolutionary biology, ecology, and biogeography. Evolutionary rates are the net result of interacting processes summarized under concepts such as adaptive radiation and evolutionary stasis. Here, we review the central concepts in the evolutionary diversification literature and synthesize these into a simple, general framework for studying rates of diversification and quantifying their underlying dynamics, which can be applied across clades and regions, and across spatial and temporal scales. Our framework describes the diversification rate (d) as a function of the abiotic environment (a), the biotic environment (b), and clade-specific phenotypes or traits (c); thus, d ~ a,b,c. We refer to the four components (a-d) and their interactions collectively as the "Evolutionary Arena." We outline analytical approaches to this framework and present a case study on conifers, for which we parameterize the general model. We also discuss three conceptual examples: the Lupinus radiation in the Andes in the context of emerging ecological opportunity and fluctuating connectivity due to climatic oscillations; oceanic island radiations in the context of island formation and erosion; and biotically driven radiations of the Mediterranean orchid genus Ophrys. The results of the conifer case study are consistent with the long-standing scenario that low competition and high rates of niche evolution promote diversification. The conceptual examples illustrate how using the synthetic Evolutionary Arena framework helps to identify and structure future directions for research on evolutionary radiations. In this way, the Evolutionary Arena framework promotes a more general understanding of variation in evolutionary rates by making quantitative results comparable between case studies, thereby allowing new syntheses of evolutionary and ecological processes to emerge.

Impacts of spatial and environmental differentiation on early Palaeozoic marine biodiversity

The unprecedented diversifications in the fossil record of the early Palaeozoic (541-419 million years ago) increased both within-sample (α) and global (γ) diversity, generating considerable ecological complexity. Faunal difference (β diversity), including spatial heterogeneity, is thought to have played a major role in early Palaeozoic marine diversification, although α diversity is the major determinant of γ diversity through the Phanerozoic. Drivers for this Phanerozoic shift from β to α diversity are not yet resolved. Here, we evaluate the impacts of environmental and faunal heterogeneity on diversity patterns using a global spatial grid. We present early Palaeozoic genus-level α, β and γ diversity curves for molluscs, brachiopods, trilobites and echinoderms and compare them with measures of spatial lithological heterogeneity, which is our proxy for environmental heterogeneity. We find that α and β diversity are associated with increased lithological heterogeneity, and that β diversity declines over time while α increases. We suggest that the enhanced dispersal of marine taxa from the Middle Ordovician onwards facilitated increases in α diversity by encouraging the occupation of narrow niches and increasing the prevalence of transient species, simultaneously reducing spatial β diversity. This may have contributed to a shift from β to α diversity as the major determinant of γ diversity increase over this critical evolutionary interval.

Beyond Reproductive Isolation: Demographic Controls on the Speciation Process

Studies of speciation typically investigate the evolution of reproductive isolation between populations, but several other processes can serve as key steps limiting the formation of species. In particular, the probability of successful speciation can be influenced by factors that affect the frequency with which population isolates form as well as their persistence through time. We suggest that population isolation and persistence have an inherently spatial dimension that can be profitably studied using a conceptual framework drawn from metapopulation ecology. We discuss models of speciation that incorporate demographic processes and highlight the need for a broader application of phylogenetic comparative approaches to evaluate the general importance of population isolation, persistence, and reproductive isolation in speciation. We review diverse and nontraditional data sources that can be leveraged to study isolation and persistence in a comparative framework. This incorporation of spatial demographic information facilitates the integration of perspectives on speciation across disciplines and timescales.

Repeated evolution of terrestrial lineages in a continental lizard radiation

The "early-burst" model of adaptive radiation predicts an early increase in phenotypic disparity concurrent with lineage diversification. Although most studies report a lack of this coupled pattern, the underlying processes are not identified. The continental radiation of Hemidactylus geckos from Peninsular India includes morphologically diverse species that occupy various microhabitats. This radiation began diversifying ~36 Mya with an early increase in lineage diversification. Here, we test the "early-burst" hypothesis by investigating the presence of ecomorphs and examining the pattern of morphological diversification in a phylogenetic framework. Two ecomorphs-terrestrial and scansorial species-that vary significantly in body size and toepad size were identified. Unlike the prediction of the "early-burst" model, we find that disparity in toepad morphology accumulated more recently ~14 Mya and fit the Ornstein-Ulhenbeck model. Ancestral state reconstruction of the two ecomorphs demonstrates that terrestrial lineages evolved independently at least five times from scansorial ancestors, with the earliest diversification in terrestrial lineages 19-12 Mya. Our study demonstrates a delayed increase in morphological disparity as a result of the evolution of terrestrial ecomorphs. The diversification of terrestrial lineages is concurrent with the establishment of open habitat and grasslands in Peninsular India, suggesting that the appearance of this novel resource led to the adaptive diversification.

Understanding the effect of competition during evolutionary radiations: an integrated model of phenotypic and species diversification

Competition can drive macroevolutionary change, for example during adaptive radiations. However, we still lack a clear understanding of how it shapes diversification processes and patterns. To better understand the macroevolutionary consequences of competition, as well as the signal left on phylogenetic data, we developed a model linking trait evolution and species diversification in an ecological context. We find four main results: first, competition spurs trait diversity but not necessarily species richness; second, competition produces slowdowns in species diversification even in the absence of explicit ecological limits, but not in phenotypic diversification even in the presence of such limits; third, early burst patterns do not provide a reliable way of testing for adaptive radiations; and fourth, looking for phylogenetic signal in trait data and support for phenotypic models incorporating competition is a better alternative. Our results clarify the macroevolutionary consequences of competition and could help design more powerful tests of adaptive radiations in nature.

Widespread Wolbachia infection in an insular radiation of damselflies (Odonata, Coenagrionidae)

Wolbachia is one of the most common endosymbionts found infecting arthropods. Theory predicts symbionts like Wolbachia will be more common in species radiations, as host shift events occur with greatest frequency between closely related species. Further, the presence of Wolbachia itself may engender reproductive isolation, and promote speciation of their hosts. Here we screened 178 individuals belonging to 30 species of the damselfly genera Nesobasis and Melanesobasis - species radiations endemic to the Fiji archipelago in the South Pacific - for Wolbachia, using multilocus sequence typing to characterize bacterial strains. Incidence of Wolbachia was 71% in Nesobasis and 40% in Melanesobasis, and prevalence was also high, with an average of 88% in the Nesobasis species screened. We identified a total of 25 Wolbachia strains, belonging to supergroups A, B and F, with some epidemic strains present in multiple species. The occurrence of Wolbachia in both males and females, and the similar global prevalence found in both sexes rules out any strong effect of Wolbachia on the primary sex-ratio, but are compatible with the phenotype of cytoplasmic incompatibility. Nesobasis has higher species richness than most endemic island damselfly genera, and we discuss the potential for endosymbiont-mediated speciation within this group.

Darwin review: angiosperm phylogeny and evolutionary radiations

Darwin's dual interests in evolution and plants formed the basis of evolutionary botany, a field that developed following his publications on both topics. Here, we review his many contributions to plant biology—from the evolutionary origins of angiosperms to plant reproduction, carnivory, and movement—and note that he expected one day there would be a ‘true’ genealogical tree for plants. This view fuelled the field of plant phylogenetics. With perhaps nearly 400 000 species, the angiosperms have diversified rapidly since their origin in the Early Cretaceous, often through what appear to be rapid radiations. We describe these evolutionary patterns, evaluate possible drivers of radiations, consider how new approaches to studies of diversification can contribute to our understanding of angiosperm diversity, and suggest new directions for further insight into plant evolution.

Speciation Rate Is Independent of the Rate of Evolution of Morphological Size, Shape, and Absolute Morphological Specialization in a Large Clade of Birds

Whether ecological differences between species evolve in parallel with lineage diversification is a fundamental issue in evolutionary biology. These processes might be connected if conditions that favor the proliferation of species, such as release from competitors, facilitate the evolution of novel ecological relationships. Despite this, phylogenetic studies do not consistently identify such a connection. Conversely, if higher diversity caused species to become increasingly specialized ecologically, then lineage diversification might become dissociated from ecological diversification. In this analysis, we ask whether the rate of lineage diversification in a large clade of birds is correlated with morphological specialization and with rates of morphological evolution. We find that morphological variation is related to species richness within clades but that rates of morphological evolution are decoupled from the rate of lineage diversification. Additionally, morphological specialization within lineages is independent of the rate at which lineages diversify, with the results apparently robust against false negative inference. This dissociation is likely a consequence of the major ecomorphological differences between avian clades arising early in their evolutionary history, with comparatively little variation added subsequently, while avian diversification has been driven predominantly by geographic isolation and sexual selection. Accordingly, biodiversity appears to be limited by the extent to which taxa can subdivide exploited regions of ecological space and not just overall ecological opportunity.

Analysis of temporal diversification of African Cyprinidae (Teleostei, Cypriniformes)

Recent evidence that freshwater fishes diversify faster than marine fishes signifies that the evolutionary history of biodiversity in freshwater system is of particular interest. Here, the evolutionary diversification events of African Cyprinidae, a freshwater fish family with wide geographic distribution, were reconstructed and analysed. The overall diversification rate of African Cyprinidae is 0.08 species per million year (when extinction rate is very high, i.e., ε = 0.9) and 0.11 species per million year (when ε = 0). This overall rate is lower than the rate reported for African Cichlids, suggesting that African freshwaters might be less conducive for a rapid diversification of Cyprinidae. However, the observed diversification events of African Cyprinidae occurred in the last 10 million years. The temporal pattern of these events follows a non-constant episodic birth-death model (Bayes Factor > 28) and the rate-constant model never outperformed any of the non-constant models tested. The fact that most diversification events occurred in the last 10 million years supports the pattern reported for Cyprinidae in other continent, e.g., Asia, perhaps pointing to concomitant diversification globally. However, the diversification events coincided with major geologic and paleo-climatic events in Africa, suggesting that geological and climatic events may have mediated the diversification patterns of Cyprinidae on the continent.

Signatures of Microevolutionary Processes in Phylogenetic Patterns

Phylogenetic trees are representations of evolutionary relationships among species and contain signatures of the processes responsible for the speciation events they display. Inferring processes from tree properties, however, is challenging. To address this problem, we analyzed a spatially-explicit model of speciation where genome size and mating range can be controlled. We simulated parapatric and sympatric (narrow and wide mating range, respectively) radiations and constructed their phylogenetic trees, computing structural properties such as tree balance and speed of diversification. We showed that parapatric and sympatric speciation are well separated by these structural tree properties. Balanced trees with constant rates of diversification only originate in sympatry and genome size affected both the balance and the speed of diversification of the simulated trees. Comparison with empirical data showed that most of the evolutionary radiations considered to have developed in parapatry or sympatry are in good agreement with model predictions. Even though additional forces other than spatial restriction of gene flow, genome size, and genetic incompatibilities, do play a role in the evolution of species formation, the microevolutionary processes modeled here capture signatures of the diversification pattern of evolutionary radiations, regarding the symmetry and speed of diversification of lineages.

The dimensionality of niche space allows bounded and unbounded processes to jointly influence diversification

There are two prominent and competing hypotheses that disagree about the effect of competition on diversification processes. The first, the bounded hypothesis, suggests that species diversity is limited (bounded) by competition between species for finite ecological niche space. The second, the unbounded hypothesis, proposes that innovations associated with evolution render competition unimportant over macroevolutionary timescales. Here we use phylogenetically structured niche modelling to show that processes consistent with both of these diversification models drive species accumulation in conifers. In agreement with the bounded hypothesis, niche competition constrained diversification, and in line with the unbounded hypothesis, niche evolution and partitioning promoted diversification. We then analyse niche traits to show that these diversification enhancing and inhibiting processes can occur simultaneously on different niche dimensions. Together these results suggest a new hypothesis for lineage diversification based on the multi-dimensional nature of ecological niches that can accommodate both bounded and unbounded evolutionary processes.

Predator macroevolution drives trophic cascades and ecosystem functioning

Biologists now recognize that ecology can drive evolution, and that evolution in turn produces ecological patterns. I extend this thinking to include longer time scales, suggesting that macroevolutionary transitions can create phenotypic differences among species, which then have predictable impacts on species interactions, community assembly and ecosystem functioning. Repeated speciation can exacerbate these patterns by creating communities with similar phenotypes and hence ecological impacts. Here, I use several experiments to test these ideas in dragonfly larvae that occupy ponds with fish, ponds without fish, or both. I show that macroevolutionary transitions between habitats cause fishless pond species to be more active relative to fish pond specialists, reducing prey abundance, shifting prey community composition and creating stronger trophic cascades. These effects scale up to the community level with predictable consequences for ecosystem multi-functioning. I suggest that macroevolutionary history can have predictable impacts on phenotypic traits, with consequences for interacting species and ecosystems.

Thermal niche evolution across replicated Anolis lizard adaptive radiations

Elucidating how ecological and evolutionary mechanisms interact to produce and maintain biodiversity is a fundamental problem in evolutionary ecology. Here, we focus on how physiological evolution affects performance and species coexistence along the thermal niche axis in replicated radiations of Anolis lizards best known for resource partitioning based on morphological divergence. We find repeated divergence in thermal physiology within these radiations, and that this divergence significantly affects performance within natural thermal environments. Morphologically similar species that co-occur invariably differ in their thermal physiology, providing evidence that physiological divergence facilitates species coexistence within anole communities. Despite repeated divergence, phylogenetic comparative analyses indicate that physiological traits have evolved more slowly than key morphological traits related to the structural niche. Phylogenetic analyses also reveal that physiological divergence is correlated with divergence in broad-scale habitat climatic features commonly used to estimate thermal niche evolution, but that the latter incompletely predicts variation in the former. We provide comprehensive evidence for repeated adaptive evolution of physiological divergence within Anolis adaptive radiations, including the complementary roles of physiological and morphological divergence in promoting community-level diversity. We recommend greater integration of performance-based traits into analyses of climatic niche evolution, as they facilitate a more complete understanding of the phenotypic and ecological consequences of climatic divergence.

The more-individuals hypothesis revisited: the role of community abundance in species richness regulation and the productivity-diversity relationship

Species richness increases with energy availability, yet there is little consensus as to the exact processes driving this species-energy relationship. The most straightforward explanation is the more-individuals hypothesis (MIH). It states that higher energy availability promotes a higher total number of individuals in a community, which consequently increases species richness by allowing for a greater number of species with viable populations. Empirical support for the MIH is mixed, partially due to the lack of proper formalisation of the MIH and consequent confusion as to its exact predictions. Here, we review the evidence of the MIH and evaluate the reliability of various predictions that have been tested. There is only limited evidence that spatial variation in species richness is driven by variation in the total number of individuals. There are also problems with measures of energy availability, with scale-dependence, and with the direction of causality, as the total number of individuals may sometimes itself be driven by the number of species. However, even in such a case the total number of individuals may be involved in diversity regulation. We propose a formal theory that encompasses these processes, clarifying how the different factors affecting diversity dynamics can be disentangled.

Anatomy of a Neotropical insect radiation

Background

Much evolutionary theory predicts that diversity arises via both adaptive radiation (diversification driven by selection against niche-overlap within communities) and divergence of geographically isolated populations. We focus on tropical fruit flies (Blepharoneura, Tephritidae) that reveal unexpected patterns of niche-overlap within local communities. Throughout the Neotropics, multiple sympatric non-interbreeding populations often share the same highly specialized patterns of host use (e.g., flies are specialists on flowers of a single gender of a single species of host plants). Lineage through time (LTT) plots can help distinguish patterns of diversification consistent with ecologically limited adaptive radiation from those predicted by ecologically neutral theories. Here, we use a time-calibrated phylogeny of Blepharoneura to test the hypothesis that patterns of Blepharoneura diversification are consistent with an “ecologically neutral” model of diversification that predicts that diversification is primarily a function of time and space.

Results

The Blepharoneura phylogeny showed more cladogenic divergence associated with geography than with shifts in host-use. Shifts in host-use were associated with ~ 20% of recent splits (< 3 Ma), but > 60% of older splits (> 3 Ma). In the overall tree, gamma statistic and maximum likelihood model fitting showed no evidence of diversification rate changes though there was a weak signature of slowing diversification rate in one of the component clades.

Conclusions

Overall patterns of Blepharoneura diversity are inconsistent with a traditional explanation of adaptive radiation involving decreases in diversification rates associated with niche-overlap. Sister lineages usually use the same host-species and host-parts, and multiple non-interbreeding sympatric populations regularly co-occur on the same hosts. We suggest that most lineage origins (phylogenetic splits) occur in allopatry, usually without shifts in host-use, and that subsequent dispersal results in assembly of communities composed of multiple sympatric non-interbreeding populations of flies that share the same hosts.

Electronic supplementary material

The online version of this article (10.1186/s12862-018-1146-9) contains supplementary material, which is available to authorized users.