Binary State Speciation And Extinction Research Articles

High-elevation-induced decrease in soil pH weakens ecosystem multifunctionality by influencing soil microbiomes

Plant environmental stress response has become a global research hotspot, yet there is a lack of clear understanding regarding the mechanisms that maintain microbial diversity and their ecosystem services under environmental stress. In our research, we examined the effects of moderate elevation on the rhizosphere soil characteristics, microbial community composition, and ecosystem multifunctionality (EMF) within agricultural systems. Our findings revealed a notable negative correlation between EMF and elevation, indicating a decline in multifunctionality at higher elevations. Additionally, our analysis across bacterial and protistan communities showed a general decrease in microbial richness with increasing elevation. Using random forest models, pH was identified as the key environmental stressor influencing microbial communities. Furthermore, we found that microbial community diversity is negatively correlated with stability by mediating complexity. Interestingly, while pH was found to affect the complexity within bacterial networks, it did not significantly impact the ecosystem stability along the elevation gradients. Using a Binary-State Speciation and Extinction (BiSSE) model to explore the evolutionary dynamics, we found that Generalists had higher speciation rates and lower extinction rates compared to specialists, resulting in a skewed distribution towards higher net diversification for generalists under increasing environmental stress. Moreover, structural equation modeling (SEM) analysis highlighted a negative correlation between environmental stress and community diversity, but showed a positive correlation between environmental stress and degree of cooperation & competition. These interactions under environmental stress indirectly increased community stability and decreased multifunctionality. Our comprehensive study offers valuable insights into the intricate relationship among environmental factors, microbial communities, and ecosystem functions, especially in the context of varying elevation gradients. These findings contribute significantly to our understanding of how environmental stressors affect microbial diversity and ecosystem services, providing a foundation for future ecological research and management strategies in similar contexts.

Patterns of bacterial generalists and specialists in lakes and reservoirs along a latitudinal gradient

AbstractAimThe geographical distribution of bacteria is an important, but poorly understood, topic in microbial ecology. A major question is how broadly distributed generalist taxa, and limitedly distributed specialist taxa, vary across a latitudinal gradient in freshwater ecosystems. We predict that: (a) generalists and specialists exhibit latitudinal diversity gradient with different patterns; (b) their community assemblies are mainly driven by stochastic processes; and (c) generalists coexist through niche differences, while specialists coexist through fitness differences.LocationSixty lakes and reservoirs in China, ranging from 18° N to 50° N.Time PeriodJuly–August in 2012 and 2018.Taxa StudiedBacterioplankton.MethodsWe clustered zero‐radius operational taxonomic unit (zOTU) (subspecies‐like level) into operational taxonomic unit (OTU) (species‐like level) and calculated the intra‐specific variation and spatial asynchrony of OTU. We examined the relative effects of environment versus space (latitude) on bacterial community and quantified the community assembly processes of generalists and specialists, respectively. The binary‐state speciation and extinction (BiSSE) model has been used to estimate the evolutionary characteristics. The abundances of the top zOTUs within each OTU were identified to reveal their niche preferences (distribution patterns at different sites).ResultsWe found that bacterial generalists had a higher intra‐specific variation, greater asynchrony and wider niche width than the specialists. With increasing latitude, bacterial diversity decreased with abundance variations in generalists and occurrence variations in specialists, respectively. Diversification in stochastic processes was the main driving factor in both groups, though deterministic processes showed greater effects on specialists than generalists in local communities. Bacterial generalists coexisted through niche differences, whereas specialists coexisted through average fitness differences.Main ConclusionsWe demonstrated distinct latitudinal distribution patterns of bacterial generalists and specialists. These results highlight differences in biogeography of microbial generalist and specialist taxa mediated by intra‐specific variation and emphasize their distinctly vital roles mediated by differing mechanisms in freshwater ecosystems.

Lifestyle Evolution Analysis by Binary-State Speciation and Extinction (BiSSE) Model.

Phylogenetic comparative methods (PCMs) combine statistics and evolutionary models to infer the dynamics of trait evolution and diversification that underlie the observed phylogeny. While PCMs have been used to study macro-evolutionary processes and evolutionary transitions of macroorganisms, their application to microbes is still limited. With the abundance of publicly available genomic and trait character data for diverse microbes nowadays, applications of PCMs on these data can provide insights into the fundamental principles that govern microbial evolution. Here, we introduce the Binary-State Speciation and Extinction (BiSSE) model, which is a relatively simple yet powerful approach for analyzing trait evolution. We begin by explaining the theoretical background and intuition behind the BiSSE model. Then, R commands for running the BiSSE model are presented. Finally, we introduce a case study that successfully applied the BiSSE model to investigate generalist and specialist microbial lifestyle evolution.

Accuracy of ancestral state reconstruction for non-neutral traits

The assumptions underpinning ancestral state reconstruction are violated in many evolutionary systems, especially for traits under directional selection. However, the accuracy of ancestral state reconstruction for non-neutral traits is poorly understood. To investigate the accuracy of ancestral state reconstruction methods, trees and binary characters were simulated under the BiSSE (Binary State Speciation and Extinction) model using a wide range of character-state-dependent rates of speciation, extinction and character-state transition. We used maximum parsimony (MP), BiSSE and two-state Markov (Mk2) models to reconstruct ancestral states. Under each method, error rates increased with node depth, true number of state transitions, and rates of state transition and extinction; exceeding 30% for the deepest 10% of nodes and highest rates of extinction and character-state transition. Where rates of character-state transition were asymmetrical, error rates were greater when the rate away from the ancestral state was largest. Preferential extinction of species with the ancestral character state also led to higher error rates. BiSSE outperformed Mk2 in all scenarios where either speciation or extinction was state dependent and outperformed MP under most conditions. MP outperformed Mk2 in most scenarios except when the rates of character-state transition and/or extinction were highly asymmetrical and the ancestral state was unfavoured.

Binary-state speciation and extinction method is conditionally robust to realistic violations of its assumptions

BackgroundPhylogenetic comparative methods allow us to test evolutionary hypotheses without the benefit of an extensive fossil record. These methods, however, make simplifying assumptions, among them that clades are always increasing or stable in diversity, an assumption we know to be false. This study simulates hypothetical clades to test whether the Binary State Speciation and Extinction (BiSSE) method can be used to correctly detect relative differences in diversification rate between ancestral and derived character states even as net diversification rates are declining overall. We simulate clades with declining but positive diversification rates, as well those in which speciation rates decline below extinction rates so that they are losing richness for part of their history. We run these analyses both with simulated symmetric and asymmetric speciation rates to test whether BiSSE can be used to detect them correctly.ResultsFor simulations with a neutral character, the fit for a BiSSE model with a neutral character is better than alternative models so long as net diversification rates remain positive. Once net diversification rates become negative, the BiSSE model with the greatest likelihood often has a non-neutral character, even though there is no such character in the simulation. BiSSE’s usefulness in detecting real asymmetry in speciation rates improves with clade age, even well after net diversification rates have become negative.ConclusionsBiSSE is most useful in analyzing clades of intermediate age, before they have reached peak diversity and gone into decline. After this point, users of BiSSE risk incorrectly inferring differential evolutionary rates when none exist. Fortunately, most studies using BiSSE and similar models focus on rapid, recent diversifications, and are less likely to encounter the biases BiSSE models are subject to for older clades. For extant groups that were once more diverse than now, however, caution should be taken in inferring past diversification patterns without fossil data.

Predicting pathogenicity behavior in Escherichia coli population through a state dependent model and TRS profiling.

The Binary State Speciation and Extinction (BiSSE) model is a branching process based model that allows the diversification rates to be controlled by a binary trait. We develop a general approach, based on the BiSSE model, for predicting pathogenicity in bacterial populations from microsatellites profiling data. A comprehensive approach for predicting pathogenicity in E. coli populations is proposed using the state-dependent branching process model combined with microsatellites TRS-PCR profiling. Additionally, we have evaluated the possibility of using the BiSSE model for estimating parameters from genetic data. We analyzed a real dataset (from 251 E. coli strains) and confirmed previous biological observations demonstrating a prevalence of some virulence traits in specific bacterial sub-groups. The method may be used to predict pathogenicity of other bacterial taxa.

Rate heterogeneity across Squamata, misleading ancestral state reconstruction and the importance of proper null model specification.

The binary-state speciation and extinction (BiSSE) model has been used in many instances to identify state-dependent diversification and reconstruct ancestral states. However, recent studies have shown that the standard procedure of comparing the fit of the BiSSE model to constant-rate birth-death models often inappropriately favours the BiSSE model when diversification rates vary in a state-independent fashion. The newly developed HiSSE model enables researchers to identify state-dependent diversification rates while accounting for state-independent diversification at the same time. The HiSSE model also allows researchers to test state-dependent models against appropriate state-independent null models that have the same number of parameters as the state-dependent models being tested. We reanalyse two data sets that originally used BiSSE to reconstruct ancestral states within squamate reptiles and reached surprising conclusions regarding the evolution of toepads within Gekkota and viviparity across Squamata. We used this new method to demonstrate that there are many shifts in diversification rates across squamates. We then fit various HiSSE submodels and null models to the state and phylogenetic data and reconstructed states under these models. We found that there is no single, consistent signal for state-dependent diversification associated with toepads in gekkotans or viviparity across all squamates. Our reconstructions show limited support for the recently proposed hypotheses that toepads evolved multiple times independently in Gekkota and that transitions from viviparity to oviparity are common in Squamata. Our results highlight the importance of considering an adequate pool of models and null models when estimating diversification rate parameters and reconstructing ancestral states.

Detecting Hidden Diversification Shifts in Models of Trait-Dependent Speciation and Extinction.

The distribution of diversity can vary considerably from clade to clade. Attempts to understand these patterns often employ state-dependent speciation and extinction models to determine whether the evolution of a particular novel trait has increased speciation rates and/or decreased extinction rates. It is still unclear, however, whether these models are uncovering important drivers of diversification, or whether they are simply pointing to more complex patterns involving many unmeasured and co-distributed factors. Here we describe an extension to the popular state-dependent speciation and extinction models that specifically accounts for the presence of unmeasured factors that could impact diversification rates estimated for the states of any observed trait, addressing at least one major criticism of BiSSE (Binary State Speciation and Extinction) methods. Specifically, our model, which we refer to as HiSSE (Hidden State Speciation and Extinction), assumes that related to each observed state in the model are "hidden" states that exhibit potentially distinct diversification dynamics and transition rates than the observed states in isolation. We also demonstrate how our model can be used as character-independent diversification models that allow for a complex diversification process that is independent of the evolution of a character. Under rigorous simulation tests and when applied to empirical data, we find that HiSSE performs reasonably well, and can at least detect net diversification rate differences between observed and hidden states and detect when diversification rate differences do not correlate with the observed states. We discuss the remaining issues with state-dependent speciation and extinction models in general, and the important ways in which HiSSE provides a more nuanced understanding of trait-dependent diversification.

Notes on the Statistical Power of the Binary State Speciation and Extinction (BiSSE) Model.

The Binary State Speciation and Extinction (BiSSE) method is one of the most popular tools for investigating the rates of diversification and character evolution. Yet, based on previous simulation studies, it is commonly held that the BiSSE method requires phylogenetic trees of fairly large sample sizes (>300 taxa) in order to distinguish between the different models of speciation, extinction, or transition rate asymmetry. Here, the power of the BiSSE method is reevaluated by simulating trees of both small and large sample sizes (30, 60, 90, and 300 taxa) under various asymmetry models and root state assumptions. Results show that the power of the BiSSE method can be much higher, also in trees of small sample size, for detecting differences in speciation rate asymmetry than anticipated earlier. This, however, is not a consequence of any conceptual or mathematical flaw in the method per se but rather of assumptions about the character state at the root of the simulated trees and thus the underlying macroevolutionary model, which led to biased results and conclusions in earlier power assessments. As such, these earlier simulation studies used to determine the power of BiSSE were not incorrect but biased, leading to an overestimation of type-II statistical error for detecting differences in speciation rate but not for extinction and transition rates.

Multiple independent origins of auto-pollination in tropical orchids (Bulbophyllum) in light of the hypothesis of selfing as an evolutionary dead end.

BackgroundThe transition from outcrossing to selfing has long been portrayed as an ‘evolutionary dead end’ because, first, reversals are unlikely and, second, selfing lineages suffer from higher rates of extinction owing to a reduced potential for adaptation and the accumulation of deleterious mutations. We tested these two predictions in a clade of Madagascan Bulbophyllum orchids (30 spp.), including eight species where auto-pollinating morphs (i.e., selfers, without a ‘rostellum’) co-exist with their pollinator-dependent conspecifics (i.e., outcrossers, possessing a rostellum). Specifically, we addressed this issue on the basis of a time-calibrated phylogeny by means of ancestral character reconstructions and within the state-dependent evolution framework of BiSSE (Binary State Speciation and Extinction), which allowed jointly estimating rates of transition, speciation, and extinction between outcrossing and selfing.ResultsThe eight species capable of selfing occurred in scattered positions across the phylogeny, with two likely originating in the Pliocene (ca. 4.4–3.1 Ma), one in the Early Pleistocene (ca. 2.4 Ma), and five since the mid-Pleistocene (ca. ≤ 1.3 Ma). We infer that this scattered phylogenetic distribution of selfing is best described by models including up to eight independent outcrossing-to-selfing transitions and very low rates of speciation (and either moderate or zero rates of extinction) associated with selfing.ConclusionsThe frequent and irreversible outcrossing-to-selfing transitions in Madagascan Bulbophyllum are clearly congruent with the first prediction of the dead end hypothesis. The inability of our study to conclusively reject or support the likewise predicted higher extinction rate in selfing lineages might be explained by a combination of methodological limitations (low statistical power of our BiSSE approach to reliably estimate extinction in small-sized trees) and evolutionary processes (insufficient time elapsed for selfers to go extinct). We suggest that, in these tropical orchids, a simple genetic basis of selfing (via loss of the ‘rostellum’) is needed to explain the strikingly recurrent transitions to selfing, perhaps reflecting rapid response to parallel and novel selective environments over Late Quaternary (≤ 1.3 Ma) time scales.Electronic supplementary materialThe online version of this article (doi:10.1186/s12862-015-0471-5) contains supplementary material, which is available to authorized users.

Lost and found: The third molars of Callimico goeldii and the evolution of the callitrichine postcanine dentition

This study tests the hypothesis that the third molars of Callimico goeldii represent a reversal in evolutionary tooth loss within the Callitrichinae. Loss of third molars is part of a suite of unusual characters that has been used to unite marmosets and tamarins in a clade to the exclusion of Callimico. However, molecular phylogenetic studies provide consistent support for the hypothesis that marmosets are more closely related to Callimico than to tamarins, raising the possibility that some or all of the features shared by marmosets and tamarins are homoplastic. Here, I use the binary-state speciation and extinction (BiSSE) model and a sample of 249 extant primate species to demonstrate that, given the shape of the primate phylogenetic tree and the distribution of character states in extant taxa, models in which M3 loss is constrained to be irreversible are much less likely than models in which reversals are allowed to occur. This result provides support for the idea that the last common ancestor of Callimico and marmosets was characterized by the two-molared phenotype. The M3s of Callimico therefore appear to be secondarily derived rather than plesiomorphic. This conclusion may also apply to the other apparently plesiomorphic traits found in Callimico. Hypotheses regarding the re-evolution of M3 in the callitrichine clade and the origin and maintenance of the two-molared phenotype are discussed.

As old as the mountains: the radiations of the Ericaceae.

Mountains are often more species-rich than lowlands. This could be the result of migration from lowlands to mountains, of a greater survival rate in mountains, or of a higher diversification rate in mountains. We investigated this question in the globally distributed family Ericaceae, which includes c. 4426 species ranging from sea level to >5000m. We predict that the interaction of low specific leaf area (SLA) and montane habitats is correlated with increased diversification rates. A molecular phylogeny of Ericaceae based on rbcL and matK sequence data was built and dated with 18 fossil calibrations and divergence time estimates. We identified radiations using bamm and correlates of diversification rate changes using binary-state speciation and extinction (BiSSE) and multiple-state speciation and extinction (MuSSE) analyses. Analyses revealed six largely montane radiations. Lineages in mountains diversified faster than nonmountain lineages (higher speciation rate, but no difference in extinction rate), and lineages with low SLA diversified faster than high-SLA lineages. Further, habitat and trait had a positive interactive effect on diversification. Our results suggest that the species richness in mountains is the result of increased speciation rather than reduced extinction or increased immigration. Increased speciation in Ericaceae was facilitated by low SLA.

Repeated evolution of tricellular (and bicellular) pollen

Male gametophytes of seed plants are sexually immature at the time they are dispersed as pollen, but approximately 30% of flowering plants have tricellular pollen containing fully formed sperm at anthesis. The classic study of Brewbaker (1967: American Journal of Botany 54: 1069-1083) provided a powerful confirmation of the long-standing hypothesis that tricellular pollen had many parallel and irreversible origins within angiosperms. We readdressed the main questions of that study with modern comparative phylogenetic methods. We used our own and more recent reports to greatly expand the Brewbaker data set. We modeled trait evolution for 2511 species on a time-calibrated angiosperm phylogeny using (1) Binary State Speciation and Extinction (BiSSE), which accounts for the effect of species diversification rates on character transition rates and, (2) the hidden rates model (HRM), which incorporates variation in transition rates across a phylogeny. Seventy percent of species had bicellular pollen. BiSSE found a 1.9-fold higher bicellular to tricellular transition rate than in the reverse direction, and bicellular lineages had a 1.8-fold higher diversification rate than tricellular lineages. HRM found heterogeneity in evolutionary rates, with bidirectional transition rates in three of four rate classes. The tricellular condition is not irreversible. Pollen cell numbers are maintained at intermediate frequencies because lower net diversification rates of tricellular lineages are counterbalanced by slower state shifts to the bicellular condition. That tricellular lineages diversify slowly and give rise to bicellular lineages slowly reflects a linkage between the evolution of sporophyte lifestyles and the developmental lability of male gametophytes.

Diversitree: comparative phylogenetic analyses of diversification in R

Summary1. The R package ‘diversitree’ contains a number of classical and contemporary comparative phylogenetic methods. Key included methods are BiSSE (binary state speciation and extinction), MuSSE (a multistate extension of BiSSE), and QuaSSE (quantitative state speciation and extinction). Diversitree also includes methods for analysing trait evolution and estimating speciation/extinction rates independently.2. In this note, I describe the features and demonstrate use of the package, using a new method, MuSSE (multistate speciation and extinction), to examine the joint effects of two traits on speciation.3. Using simulations, I found that MuSSE could reliably detect that a binary trait that affected speciation rates when simultaneously accounting for additional thats that had no effect on speciation rates.4. Diversitree is an open source and available on the Comprehensive R Archive Network (cran). A tutorial and worked examples can be downloaded from http://www.zoology.ubc.ca/prog/diversitree.

ON THE REVERSIBILITY OF MANDIBULAR SYMPHYSEAL FUSION

Experimental and comparative studies suggest that a major determinant of increased ossification of the mandibular symphysis is elevated masticatory stress related to a mechanically challenging diet. However, the morphology of this joint tracks variation in dietary properties in only some mammalian clades. Extant anthropoid primates are a notable exception: synostosis is ubiquitous in this speciose group, despite its great age and diverse array of feeding adaptations. One possible explanation for this pattern is that, once synostosis evolves, reversion to a lesser degree of fusion is unlikely or even constrained. If correct, this has important implications for functional and phylogenetic analyses of the mammalian feeding apparatus. To test this hypothesis, we generated a molecular tree for 76 vespertilionoid and noctilionoid chiropterans using Bayesian phylogenetic analysis and examined character evolution using parsimony and likelihood ancestral-state reconstructions along with the binary state speciation and extinction (BiSSE) model. Results indicate that reversals have occurred within Vespertilionoidea. In contrast, noctilionoids exhibit an anthropoid-like pattern, which suggests that more detailed comparisons of the functional and developmental bases for fusion in these bat clades may provide insight into why fusion is maintained in some lineages but not in others. Potential functional and developmental explanations for the lack of reversal are discussed.

LOSS OF SEXUAL RECOMBINATION AND SEGREGATION IS ASSOCIATED WITH INCREASED DIVERSIFICATION IN EVENING PRIMROSES

The loss of sexual recombination and segregation in asexual organisms has been portrayed as an irreversible process that commits asexually reproducing lineages to reduced diversification. We test this hypothesis by estimating rates of speciation, extinction, and transition between sexuality and functional asexuality in the evening primroses. Specifically, we estimate these rates using the recently developed BiSSE (Binary State Speciation and Extinction) phylogenetic comparative method, which employs maximum likelihood and Bayesian techniques. We infer that net diversification rates (speciation minus extinction) in functionally asexual evening primrose lineages are roughly eight times faster than diversification rates in sexual lineages, largely due to higher speciation rates in asexual lineages. We further reject the hypothesis that a loss of recombination and segregation is irreversible because the transition rate from functional asexuality to sexuality is significantly greater than zero and in fact exceeded the reverse rate. These results provide the first empirical evidence in support of the alternative theoretical prediction that asexual populations should instead diversify more rapidly than sexual populations because they are free from the homogenizing effects of sexual recombination and segregation. Although asexual reproduction may often constrain adaptive evolution, our results show that the loss of recombination and segregation need not be an evolutionary dead end in terms of diversification of lineages.