Rates of species origination and extinction can vary over time during evolutionary radiations, and it is possible to reconstruct the history of diversification using molecular phylogenies of extant taxa only. Maximum likelihood methods provide a useful framework for inferring temporal variation in diversification rates. LASER is a package for the R programming environment that implements maximum likelihood methods based on the birth-death process to test whether diversification rates have changed over time. LASER contrasts the likelihood of phylogenetic data under models where diversification rates have changed over time to alternative models where rates have remained constant over time. Major strengths of the package include the ability to detect temporal increases in diversification rates and the inference of diversification parameters under multiple rate-variable models of diversification. The program and associated documentation are freely available from the R package archive at http://cran.r-project.org.
Time-calibrated phylogenies of living species have been widely used to study the tempo and mode of species diversification. However, it is increasingly clear that inferences about species diversification-extinction rates in particular-can be unreliable in the absence of paleontological data. We introduce a general framework based on the fossilized birth-death process for studying speciation-extinction dynamics on phylogenies of extant and extinct species. The model assumes that phylogenies can be modeled as a mixture of distinct evolutionary rate regimes and that a hierarchical Poisson process governs the number of such rate regimes across a tree. We implemented the model in BAMM, a computational framework that uses reversible jump Markov chain Monte Carlo to simulate a posterior distribution of macroevolutionary rate regimes conditional on the branching times and topology of a phylogeny. The implementation, we describe can be applied to paleontological phylogenies, neontological phylogenies, and to phylogenies that include both extant and extinct taxa. We evaluate performance of the model on data sets simulated under a range of diversification scenarios. We find that speciation rates are reliably inferred in the absence of paleontological data. However, the inclusion of fossil observations substantially increases the accuracy of extinction rate estimates. We demonstrate that inferences are relatively robust to at least some violations of model assumptions, including heterogeneity in preservation rates and misspecification of the number of occurrences in paleontological data sets.
For many species, both local abundance and regional occupancy are highest near the centre of their geographic distributions. One hypothesis for this pattern is that niche suitability declines with increasing distance from a species geographic centre, such that populations near range margins are characterized by reduced density and increased patchiness. In these smaller edge populations, genetic drift is more powerful, leading to the loss of genetic diversity. This simple verbal model has been formalized as the central-marginal hypothesis, which predicts that core populations should have greater genetic diversity than edge populations. Here, we tested the central-marginal hypothesis using a genomic data set of 25 species-level taxa of Australian scincid lizards in the genera Ctenotus and Lerista. A majority of taxa in our data set showed range-wide patterns of genetic variation consistent with central-marginal hypothesis, and eight of 25 taxa showed significantly greater genetic diversity in the centre of their range. We then explored biological, historical, and methodological factors that might predict which taxa support the central-marginal hypothesis. We found that taxa with the strongest evidence for range expansion were the least likely to follow predictions of the central-marginal hypothesis. The majority of these taxa had range expansions that originated at the range edge, which led to a gradient of decreasing genetic diversity from the range edge to the core, contrary to the central-marginal hypothesis.
Advances in the generation, retrieval, and analysis of phylogenetic data have enabled researchers to create phylogenies that contain many thousands of taxa. These "macrophylogenies"-large trees that typically derive from megaphylogeny, supermatrix, or supertree approaches-provide researchers with an unprecedented ability to conduct evolutionary analyses across broad phylogenetic scales. Many studies have now used these phylogenies to explore the dynamics of speciation, extinction, and phenotypic evolution across large swaths of the tree of life. These trees are characterized by substantial phylogenetic uncertainty on multiple levels, and the stability of macroevolutionary inferences from these data sets has not been rigorously explored. As a case study, we tested whether five recently published phylogenies for squamate reptiles-each consisting of more than 4000 species-yield congruent inferences about the processes that underlie variation in species richness across replicate evolutionary radiations of Australian snakes and lizards. We find discordance across the five focal phylogenies with respect to clade age and several diversification rate metrics, and in the effects of clade age on species richness. We also find that crown clade ages reported in the literature on these Australian groups are in conflict with all of the large phylogenies examined. Macrophylogenies offer an unprecedented opportunity to address evolutionary and ecological questions at broad phylogenetic scales, but accurately representing the uncertainty that is inherent to such analyses remains a critical challenge to our field. [Australia; macroevolution; macrophylogeny; squamates; time calibration.].
Summary Understanding the dynamics of speciation, extinction and phenotypic evolution is a central challenge in evolutionary biology. Here, we present BAMM tools, an r package for the analysis and visualization of macroevolutionary dynamics on phylogenetic trees. BAMM tools is a companion package to BAMM , an open‐source program for reversible‐jump MCMC analyses of diversification and trait evolution. Functions in BAMM tools operate directly on output from the BAMM program. The package is oriented towards reconstructing and visualizing changes in evolutionary rates through time and across clades in a Bayesian statistical framework. BAMM tools enables users to extract credible sets of diversification shifts and to identify diversification histories with the maximum a posteriori probability. Users can compare the fit of alternative diversification models using Bayes factors and by directly comparing model posterior probabilities. By providing a robust framework for quantifying uncertainty in macroevolutionary dynamics, BAMM tools will facilitate inference on the complex mixture of processes that have shaped the distribution of species and phenotypes across the tree of life.
A bstract Rates of character evolution in macroevolutionary datasets are typically estimated by maximizing the likelihood function of a continuous-time Markov chain (CTMC) model of character evolution over all possible histories of character state change, a technique known as maximum average likelihood. An alternative approach is to estimate ancestral character states independently of rates using parsimony and to then condition likelihood-based estimates of transition rates on the resulting ancestor-descendant reconstructions. We use maximum parsimony reconstructions of possible pathways of evolution to implement this alternative approach for single-character datasets simulated on empirical phylogenies using a two-state CTMC. We find that transition rates estimated using parsimonious ancestor-descendant reconstructions have lower mean squared error than transition rates estimated by maximum average likelihood. Although we use a binary state character for exposition, the approach remains valid for an arbitrary number of states. Finally, we show how this method can be used to rapidly and easily detect phylogenetic variation in tempo and mode of character evolution with two empirical examples from squamates. These results highlight the mutually informative roles of parsimony and likelihood when testing hypotheses of character evolution in macroevolution.
