What affects power to estimate speciation rate shifts?

Exploring the evolutionary dynamics of myrmecophytism: Perspectives from the Southeast Asian Macaranga ant-plant symbiosis

Myrmecophytic plants utilise defensive services offered by obligate ant partners nesting in their domatia in a novel means of survival in tropical habitats. Although much is known about the ecology of myrmecophytism, there aren't enough empirical examples to demonstrate whether it substantially influences evolutionary patterns in host plant lineages. In this study, we make use of the species-rich Macaranga (Euphorbiaceae) ant-plant symbiosis distributed in the Southeast Asian Sundaland to delve into the evolutionary dynamics of myrmecophytism in host plants. We generated the most comprehensive dated phylogeny of myrmecophytic Macaranga till date using genotyping-by-sequencing (GBS). With this in hand, we traced the evolutionary history of myrmecophytism in Macaranga using parametric biogeography and ancestral state reconstruction. Diversification rate analysis methods were employed to determine if myrmecophytism enhanced diversification rates in the genus. Our results demonstrate that myrmecophytism is labile and easily lost. Ancestral state reconstruction supported a single origin of myrmecophytism in Macaranga ∼18 mya on Borneo followed by multiple losses. Diversification rate analysis methods did not yield sufficient evidence to support the hypothesis that myrmecophytism enhanced diversification rates in Macaranga; we found that topographical features on Borneo may have played a more direct role in the divergence of clades instead. Our study provides evidence that while the acquisition of domatia clearly functions as a key innovation that has enabled host plants to exploit the environment in novel ways, it may not necessarily enhance diversification rates. In fact, we hypothesise that overly specialised cases of myrmecophytism may even be an evolutionary dead end.

Chemical signaling glands are unlinked to species diversification in lizards

Sexual selection has long been thought to increase species diversification. Sexually selected traits, such as sexual signals that contribute to reproductive isolation, were thought to promote diversification. However, studies exploring links between sexually selected traits and species diversification have thus far primarily focused on visual or acoustic signals. Many animals often employ chemical signals (i.e., pheromones) for sexual communications, but large-scale analyses on the role of chemical communications in driving species diversification have been missing. Here, for the first time, we investigate whether traits associated with chemical communications-the presence of follicular epidermal glands-promote diversification across 6,672 lizard species. In most analyses, we found no strong association between the presence of follicular epidermal glands and species diversification rates, either across all lizard species or at lower phylogenetic scales. Previous studies suggest that follicular gland secretions act as species recognition signals that prevent hybridization during speciation in lizards. However, we show that geographic range overlap was no different in sibling species pairs with and without follicular epidermal glands. Together, these results imply that either follicular epidermal glands do not primarily function in sexual communications or sexually selected traits in general (here chemical communication) have a limited effect on species diversification. In our additional analysis accounting for sex-specific differences in glands, we again found no detectable effect of follicular epidermal glands on species diversification rates. Thus, our study challenges the general role of sexually selected traits in broad-scale species diversification patterns.

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.

Exceptional evolutionary lability of flower-like inflorescences (pseudanthia) in Apiaceae subfamily Apioideae

PREMISE

Pseudanthia are widespread and have long been postulated to be a key innovation responsible for some of the angiosperm radiations. The aim of our study was to analyze macroevolutionary patterns of these flower-like inflorescences and their potential correlation with diversification rates in Apiaceae subfamily Apioideae. In particular, we were interested to investigate evolvability of pseudanthia and evaluate their potential association with changes in the size of floral display.

METHODS

The framework for our analyses consisted of a time-calibrated phylogeny of 1734 representatives of Apioideae and a morphological matrix of inflorescence traits encoded for 847 species. Macroevolutionary patterns in pseudanthia were inferred using Markov models of discrete character evolution and stochastic character mapping, and a principal component analysis was used to visualize correlations in inflorescence architecture. The interdependence between net diversification rates and the occurrence of pseudocorollas was analyzed with trait-independent and trait-dependent approaches.

RESULTS

Pseudanthia evolved in 10 major clades of Apioideae with at least 36 independent origins and 46 reversals. The morphospace analysis recovered differences in color and compactness between floral and hyperfloral pseudanthia. A correlation between pseudocorollas and size of inflorescence was also strongly supported. Contrary to our predictions, pseudanthia are not responsible for variation in diversification rates identified in this subfamily.

CONCLUSIONS

Our results suggest that pseudocorollas evolve as an answer to the trade-off between enlargement of floral display and costs associated with production of additional flowers. The high evolvability and architectural differences in apioid pseudanthia may be explained on the basis of adaptive wandering and evolutionary developmental biology.

Evolution of diverse host infection mechanisms delineates an adaptive radiation of lampsiline freshwater mussels centered on their larval ecology

North American watersheds contain a high diversity of freshwater mussels (Unionoida). During the long-lived, benthic phase of their life cycle, up to 40 species can co-occur in a single riffle and there is typically little evidence for major differences in their feeding ecology or microhabitat partitioning. In contrast, their brief parasitic larval phase involves the infection of a wide diversity of fish hosts and female mussels have evolved a spectrum of adaptations for infecting host fish with their offspring. Many species use a passive broadcast strategy: placing high numbers of larvae in the water column and relying on chance encounters with potential hosts. Many other species, including most members of the Lampsilini, have a proactive strategy that entails the use of prey-mimetic lures to change the behavior of the hosts, i.e., eliciting a feeding response through which they become infected. Two main lure types are collectively produced: mantle tissue lures (on the female’s body) and brood lures, containing infective larvae, that are released into the external environment. In this study, we used a phylogenomic approach (ddRAD-seq) to place the diversity of infection strategies used by 54 North American lampsiline mussels into an evolutionary context. Ancestral state reconstruction recovered evidence for the early evolution of mantle lures in this clade, with brood lures and broadcast infection strategies both being independently derived twice. The most common infection strategy, occurring in our largest ingroup clade, is a mixed one in which mimetic mantle lures are apparently the predominant infection mechanism, but gravid females also release simple, non-mimetic brood lures at the end of the season. This mixed infection strategy clade shows some evidence of an increase in diversification rate and most members use centrarchids (Micropterus & Lepomis spp.) as their predominant fish hosts. Broad linkage between infection strategies and predominant fish host genera is also seen in other lampsiline clades: worm-like mantle lures of Toxolasma spp. with sunfish (Lepomis spp.); insect larvae-like brood lures (Ptychobranchus spp.), or mantle lures (Medionidus spp., Obovaria spp.), or mantle lures combined with host capture (Epioblasma spp.) with a spectrum of darter (Etheostoma & Percina spp.) and sculpin (Cottus spp.) hosts, and tethered brood lures (Hamiota spp.) with bass (Micropterus spp.). Our phylogenetic results confirm that discrete lampsiline mussel clades exhibit considerable specialization in the primary fish host clades their larvae parasitize, and in the host infection strategies they employ to do so. They are also consistent with the hypothesis that larval resource partitioning of fish hosts is an important factor in maintaining species diversity in mussel assemblages. We conclude that, taking their larval ecology and host-infection mechanisms into account, lampsiline mussels may be legitimately viewed as an adaptive radiation.

Molecules and fossils tell distinct yet complementary stories of mammal diversification

Reconstructing the tempo at which biodiversity arose is a fundamental goal of evolutionary biologists, yet the relative merits of evolutionary-rate estimates are debated based on whether they are derived from the fossil record or time-calibrated phylogenies (timetrees) of living species. Extinct lineages unsampled in timetrees are known to "pull" speciation rates downward, but the temporal scale at which this bias matters is unclear. To investigate this problem, we compare mammalian diversification-rate signatures in a credible set of molecular timetrees (n = 5,911 species, ∼70% from DNA) to those in fossil genus durations (n = 5,320). We use fossil extinction rates to correct or "push" the timetree-based (pulled) speciation-rate estimates, finding a surge of speciation during the Paleocene (∼66-56 million years ago, Ma) between the Cretaceous-Paleogene (K-Pg) boundary and the Paleocene-Eocene Thermal Maximum (PETM). However, about two-thirds of the K-Pg-to-PETM originating taxa did not leave modern descendants, indicating that this rate signature is likely undetectable from extant lineages alone. For groups without substantial fossil records, thankfully all is not lost. Pushed and pulled speciation rates converge starting ∼10 Ma and are equal at the present day when recent evolutionary processes can be estimated without bias using species-specific "tip" rates of speciation. Clade-wide moments of tip rates also enable enriched inference, as the skewness of tip rates is shown to approximate a clade's extent of past diversification-rate shifts. Molecular timetrees need fossil-correction to address deep-time questions, but they are sufficient for shallower time questions where extinctions are fewer.

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.

Exploring diversification drivers in golden orbweavers

Heterogeneity in species diversity is driven by the dynamics of speciation and extinction, potentially influenced by organismal and environmental factors. Here, we explore macroevolutionary trends on a phylogeny of golden orbweavers (spider family Nephilidae). Our initial inference detects heterogeneity in speciation and extinction, with accelerated extinction rates in the extremely sexually size dimorphic Nephila and accelerated speciation in Herennia, a lineage defined by highly derived, arboricolous webs, and pronounced island endemism. We evaluate potential drivers of this heterogeneity that relate to organisms and their environment. Primarily, we test two continuous organismal factors for correlation with diversification in nephilids: phenotypic extremeness (female and male body length, and sexual size dimorphism as their ratio) and dispersal propensity (through range sizes as a proxy). We predict a bell-shaped relationship between factor values and speciation, with intermediate phenotypes exhibiting highest diversification rates. Analyses using SSE-class models fail to support our two predictions, suggesting that phenotypic extremeness and dispersal propensity cannot explain patterns of nephilid diversification. Furthermore, two environmental factors (tropical versus subtropical and island versus continental species distribution) indicate only marginal support for higher speciation in the tropics. Although our results may be affected by methodological limitations imposed by a relatively small phylogeny, it seems that the tested organismal and environmental factors play little to no role in nephilid diversification. In the phylogeny of golden orbweavers, the recent hypothesis of universal diversification dynamics may be the simplest explanation of macroevolutionary patterns.

Geo-Climatic Changes and Apomixis as Major Drivers of Diversification in the Mediterranean Sea Lavenders (Limonium Mill.)

The Mediterranean realm, comprising the Mediterranean and Macaronesian regions, has long been recognized as one of the world's biodiversity hotspots, owing to its remarkable species richness and endemism. Several hypotheses on biotic and abiotic drivers of species diversification in the region have been often proposed but rarely tested in an explicit phylogenetic framework. Here, we investigate the impact of both species-intrinsic and -extrinsic factors on diversification in the species-rich, cosmopolitan Limonium, an angiosperm genus with center of diversity in the Mediterranean. First, we infer and time-calibrate the largest Limonium phylogeny to date. We then estimate ancestral ranges and diversification dynamics at both global and regional scales. At the global scale, we test whether the identified shifts in diversification rates are linked to specific geological and/or climatic events in the Mediterranean area and/or asexual reproduction (apomixis). Our results support a late Paleogene origin in the proto-Mediterranean area for Limonium, followed by extensive in situ diversification in the Mediterranean region during the late Miocene, Pliocene, and Pleistocene. We found significant increases of diversification rates in the "Mediterranean lineage" associated with the Messinian Salinity Crisis, onset of Mediterranean climate, Plio-Pleistocene sea-level fluctuations, and apomixis. Additionally, the Euro-Mediterranean area acted as the major source of species dispersals to the surrounding areas. At the regional scale, we infer the biogeographic origins of insular endemics in the oceanic archipelagos of Macaronesia, and test whether woodiness in the Canarian Nobiles clade is a derived trait linked to insular life and a biotic driver of diversification. We find that Limonium species diversity on the Canary Islands and Cape Verde archipelagos is the product of multiple colonization events followed by in situ diversification, and that woodiness of the Canarian endemics is indeed a derived trait but is not associated with a significant shift to higher diversification rates. Our study expands knowledge on how the interaction between abiotic and biotic drivers shape the uneven distribution of species diversity across taxonomic and geographical scales.

The early wasp plucks the flower: disparate extant diversity of sawfly superfamilies (Hymenoptera: ‘Symphyta’) may reflect asynchronous switching to angiosperm hosts

The insect order Hymenoptera originated during the Permian nearly 300 Mya. Ancestrally herbivorous hymenopteran lineages today make up the paraphyletic suborder ‘Symphyta’, which encompasses c. 8200 species with very diverse host-plant associations. We use phylogeny-based statistical analyses to explore the drivers of diversity dynamics within the ‘Symphyta’, with a particular focus on the hypothesis that diversification of herbivorous insects has been driven by the explosive radiation of angiosperms during and after the Cretaceous. Our ancestral-state estimates reveal that the first symphytans fed on gymnosperms, and that shifts onto angiosperms and pteridophytes – and back – have occurred at different time intervals in different groups. Trait-dependent analyses indicate that average net diversification rates do not differ between symphytan lineages feeding on angiosperms, gymnosperms or pteridophytes, but trait-independent models show that the highest diversification rates are found in a few angiosperm-feeding lineages that may have been favoured by the radiations of their host taxa during the Cenozoic. Intriguingly, lineages-through-time plots show signs of an early Cretaceous mass extinction, with a recovery starting first in angiosperm-associated clades. Hence, the oft-invoked assumption of herbivore diversification driven by the rise of flowering plants may overlook a Cretaceous global turnover in insect herbivore communities during the rapid displacement of gymnosperm- and pteridophyte-dominated floras by angiosperms.

Molecular phylogeny and historical biogeography of Parnara butterflies (Lepidoptera: Hesperiidae)

The butterfly genus Parnara (Hesperiinae: Baorini), of which some are major pests of economic crops (e.g., rice, wild rice stems and sugarcane), currently consists of 10 species and several subspecies and has a highly disjunct distribution in Australia, Africa, and Asia. We determined the systematic relationships and biogeographical history of the genus by reconstructing the phylogeny based on eight genes and 101 specimens representing all 10 recognized species. Four species delimitation methods (ABGD, bPTP, GMYC and BPP) were also employed to assess the taxonomic status of each species. Based on these results and analyses, we recognize 11 extant species in the genus. The status of the taxon P. naso poutieri (Boisduval, 1833) from Madagascar is revised as a distinct species, Parnara poutieri (Boisduval, 1833) stat. rev. The subspecies P. guttata mangala (Moore, 1866) syn. nov. is synonymized with P. guttata guttata (Bremer & Grey, 1853), while P. bada (Moore, 1878) is provisionally treated as a complex of two species, namely P. bada and P. apostata (Snellen, 1886). The monophyly of Parnara is strongly supported, with the following relationships: P. amalia + ((P. monasi + (P. poutieri + P. naso)) + ((P. kawazoei + P. bada complex) + (P. ganga + (P. ogasawarensis + (P. guttata + P. batta))))). Divergence time and ancestral range estimates indicate that the common ancestor of Parnara originated in an implausible area of Australia, Africa, and Oriental region in the mid-Oligocene and then differentiated in the late Miocene-late Pliocene. Dispersal and range expansion have played an important role in diversification of the genus in Asia and Afica. Relatively stable geotectonic plates at the time when most extant lineages appeared during the late Miocene-early Pliocene might have been the factor responsible for the relatively constant low dynamic rate of diversification within the group.

Clade-age-dependent diversification under high species turnover shapes species richness disparities among tropical rainforest lineages of Bulbophyllum (Orchidaceae)

Background

Tropical rainforests (TRFs) harbour almost half of the world’s vascular plant species diversity while covering only about 6–7% of land. However, why species richness varies amongst the Earth’s major TRF regions remains poorly understood. Here we investigate the evolutionary processes shaping continental species richness disparities of the pantropical, epiphytic and mostly TRF-dwelling orchid mega-genus Bulbophyllum (c. 1948 spp. in total) using diversification analyses based on a time-calibrated molecular phylogeny (including c. 45–50% spp. each from Madagascar, Africa, Neotropics, and 8.4% from the Asia-Pacific region), coupled with ecological niche modelling (ENM) of geographic distributions under present and past (Last Glacial Maximum; LGM) conditions.

Results

Our results suggest an early-to-late Miocene scenario of ‘out-of-Asia-Pacific’ origin and progressive, dispersal-mediated diversification in Madagascar, Africa and the Neotropics, respectively. Species richness disparities amongst these four TRF lineages are best explained by a time-for-speciation (i.e. clade age) effect rather than differences in net diversification or diversity-dependent diversification due to present or past spatial-bioclimatic limits. For each well-sampled lineage (Madagascar, Africa, Neotropics), we inferred high rates of speciation and extinction over time (i.e. high species turnover), yet with the origin of most extant species falling into the Quaternary. In contrast to predictions of classical ‘glacial refuge’ theories, all four lineages experienced dramatic range expansions during the LGM.

Conclusions

As the Madagascan, African and Neotropical lineages display constant-rate evolution since their origin (early-to-mid-Miocene), Quaternary environmental change might be a less important cause of their high species turnover than intrinsic features generally conferring rapid population turnover in tropical orchids (e.g., epiphytism, specialization on pollinators and mycorrhizal fungi, wind dispersal). Nonetheless, climate-induced range fluctuations during the Quaternary could still have played an influential role in the origination and extinction of Bulbophyllum species in those three, if not in all four TRF regions.

Electronic supplementary material

The online version of this article (10.1186/s12862-019-1416-1) contains supplementary material, which is available to authorized users.

Live fast, diversify non-adaptively: evolutionary diversification of exceptionally short-lived annual killifishes

Background

Adaptive radiations are triggered by ecological opportunity – the access to novel niche domains with abundant available resources that facilitate the formation of new ecologically divergent species. Therefore, as new species saturate niche space, clades experience a diversity-dependent slowdown of diversification over time. At the other extreme of the radiation continuum, non-adaptively radiating lineages undergo diversification with minimal niche differentiation when ‘spatial opportunity’ (i.e. areas with suitable ‘ancestral’ ecological conditions) is available. Traditionally, most research has focused on adaptive radiations, while empirical studies on non-adaptive radiations remain lagging behind. A prolific clade of African fish with extremely short lifespan (Nothobranchius killifish), show the key evolutionary features of a candidate non-adaptive radiation – primarily allopatric species with minimal niche and phenotypic divergence. Here, we test the hypothesis that Nothobranchius killifish have non-adaptively diversified. We employ phylogenetic modelling to investigate the tempo and mode of macroevolutionary diversification of these organisms.

Results

Nothobranchius diversification has proceeded with minor niche differentiation and minimal morphological disparity among allopatric species. Additionally, we failed to identify evidence for a role of body size or biogeography in influencing diversification rates. Diversification has been homogeneous within this genus, with the only hotspot of species-richness not resulting from rapid diversification. However, species in sympatry show higher disparity, which may have been caused by character displacement among coexisting species.

Conclusions

Nothobranchius killifish have proliferated following the tempo and mode of a non-adaptive radiation. Our study confirms that this exceptionally short-lived group have diversified with minimal divergent niche adaptation, while one group of coexisting species seems to have facilitated spatial overlap among these taxa via the evolution of ecological character displacement.

Electronic supplementary material

The online version of this article (10.1186/s12862-019-1344-0) contains supplementary material, which is available to authorized users.