Assessing Rapid Relaxed-Clock Methods for Phylogenomic Dating.
Ontology highlight
ABSTRACT:
SUBMITTER: Barba-Montoya J
PROVIDER: S-EPMC8633771 | biostudies-literature |
REPOSITORIES: biostudies-literature
Similar Datasets
Project description:Relaxed clock models enable estimation of molecular substitution rates across lineages and are widely used in phylogenetics for dating evolutionary divergence times. Under the (uncorrelated) relaxed clock model, tree branches are associated with molecular substitution rates which are independently and identically distributed. In this article we delved into the internal complexities of the relaxed clock model in order to develop efficient MCMC operators for Bayesian phylogenetic inference. We compared three substitution rate parameterisations, introduced an adaptive operator which learns the weights of other operators during MCMC, and we explored how relaxed clock model estimation can benefit from two cutting-edge proposal kernels: the AVMVN and Bactrian kernels. This work has produced an operator scheme that is up to 65 times more efficient at exploring continuous relaxed clock parameters compared with previous setups, depending on the dataset. Finally, we explored variants of the standard narrow exchange operator which are specifically designed for the relaxed clock model. In the most extreme case, this new operator traversed tree space 40% more efficiently than narrow exchange. The methodologies introduced are adaptive and highly effective on short as well as long alignments. The results are available via the open source optimised relaxed clock (ORC) package for BEAST 2 under a GNU licence (https://github.com/jordandouglas/ORC).
| S-EPMC7880504 | biostudies-literature
Project description:Phylogenetic dating is one of the most powerful and commonly used methods of drawing epidemiological interpretations from pathogen genomic data. Building such trees requires considering a molecular clock model which represents the rate at which substitutions accumulate on genomes. When the molecular clock rate is constant throughout the tree then the clock is said to be strict, but this is often not an acceptable assumption. Alternatively, relaxed clock models consider variations in the clock rate, often based on a distribution of rates for each branch. However, we show here that the distributions of rates across branches in commonly used relaxed clock models are incompatible with the biological expectation that the sum of the numbers of substitutions on two neighboring branches should be distributed as the substitution number on a single branch of equivalent length. We call this expectation the additivity property. We further show how assumptions of commonly used relaxed clock models can lead to estimates of evolutionary rates and dates with low precision and biased confidence intervals. We therefore propose a new additive relaxed clock model where the additivity property is satisfied. We illustrate the use of our new additive relaxed clock model on a range of simulated and real data sets, and we show that using this new model leads to more accurate estimates of mean evolutionary rates and ancestral dates.
| S-EPMC8480190 | biostudies-literature
Project description:Evolutionary timescales can be inferred from molecular sequence data using a Bayesian phylogenetic approach. In these methods, the molecular clock is often calibrated using fossil data. The uncertainty in these fossil calibrations is important because it determines the limiting posterior distribution for divergence-time estimates as the sequence length tends to infinity. Here, we investigate how the accuracy and precision of Bayesian divergence-time estimates improve with the increased clock-partitioning of genome-scale data into clock-subsets. We focus on a data set comprising plastome-scale sequences of 52 angiosperm taxa. There was little difference among the Bayesian date estimates whether we chose clock-subsets based on patterns of among-lineage rate heterogeneity or relative rates across genes, or by random assignment. Increasing the degree of clock-partitioning usually led to an improvement in the precision of divergence-time estimates, but this increase was asymptotic to a limit presumably imposed by fossil calibrations. Our clock-partitioning approaches yielded highly precise age estimates for several key nodes in the angiosperm phylogeny. For example, when partitioning the data into 20 clock-subsets based on patterns of among-lineage rate heterogeneity, we inferred crown angiosperms to have arisen 198-178 Ma. This demonstrates that judicious clock-partitioning can improve the precision of molecular dating based on phylogenomic data, but the meaning of this increased precision should be considered critically.
| S-EPMC5647803 | biostudies-literature
Project description:In phylogenetics, the unrooted model of phylogeny and the strict molecular clock model are two extremes of a continuum. Despite their dominance in phylogenetic inference, it is evident that both are biologically unrealistic and that the real evolutionary process lies between these two extremes. Fortunately, intermediate models employing relaxed molecular clocks have been described. These models open the gate to a new field of "relaxed phylogenetics." Here we introduce a new approach to performing relaxed phylogenetic analysis. We describe how it can be used to estimate phylogenies and divergence times in the face of uncertainty in evolutionary rates and calibration times. Our approach also provides a means for measuring the clocklikeness of datasets and comparing this measure between different genes and phylogenies. We find no significant rate autocorrelation among branches in three large datasets, suggesting that autocorrelated models are not necessarily suitable for these data. In addition, we place these datasets on the continuum of clocklikeness between a strict molecular clock and the alternative unrooted extreme. Finally, we present analyses of 102 bacterial, 106 yeast, 61 plant, 99 metazoan, and 500 primate alignments. From these we conclude that our method is phylogenetically more accurate and precise than the traditional unrooted model while adding the ability to infer a timescale to evolution.
| S-EPMC1395354 | biostudies-literature
Project description:Over recent years, several alternative relaxed clock models have been proposed in the context of Bayesian dating. These models fall in two distinct categories: uncorrelated and autocorrelated across branches. The choice between these two classes of relaxed clocks is still an open question. More fundamentally, the true process of rate variation may have both long-term trends and short-term fluctuations, suggesting that more sophisticated clock models unfolding over multiple time scales should ultimately be developed. Here, a mixed relaxed clock model is introduced, which can be mechanistically interpreted as a rate variation process undergoing short-term fluctuations on the top of Brownian long-term trends. Statistically, this mixed clock represents an alternative solution to the problem of choosing between autocorrelated and uncorrelated relaxed clocks, by proposing instead to combine their respective merits. Fitting this model on a dataset of 105 placental mammals, using both node-dating and tip-dating approaches, suggests that the two pure clocks, Brownian and white noise, are rejected in favour of a mixed model with approximately equal contributions for its uncorrelated and autocorrelated components. The tip-dating analysis is particularly sensitive to the choice of the relaxed clock model. In this context, the classical pure Brownian relaxed clock appears to be overly rigid, leading to biases in divergence time estimation. By contrast, the use of a mixed clock leads to more recent and more reasonable estimates for the crown ages of placental orders and superorders. Altogether, the mixed clock introduced here represents a first step towards empirically more adequate models of the patterns of rate variation across phylogenetic trees.This article is part of the themed issue 'Dating species divergences using rocks and clocks'.
| S-EPMC4920333 | biostudies-literature
Project description:The rapid expansion of sequence data and the development of statistical approaches that embrace varying evolutionary rates among lineages have encouraged many more investigators to use DNA and protein data to time species divergences. Here, we report results from a systematic evaluation, by means of computer simulation, of the performance of two frequently used relaxed-clock methods for estimating these times and their credibility intervals (CrIs). These relaxed-clock methods allow rates to vary in a phylogeny randomly over lineages (e.g., BEAST software) and in autocorrelated fashion (e.g., MultiDivTime software). We applied these methods for analyzing sequence data sets simulated using naturally derived parameters (evolutionary rates, sequence lengths, and base substitution patterns) and assuming that clock calibrations are known without error. We find that the estimated times are, on average, close to the true times as long as the assumed model of lineage rate changes matches the actual model. The 95% CrIs also contain the true time for >or=95% of the simulated data sets. However, the use of incorrect lineage rate model reduces this frequency to 83%, indicating that the relaxed-clock methods are not robust to the violation of underlying lineage rate model. Because these rate models are rarely known a priori and are difficult to detect empirically, we suggest building composite CrIs using CrIs produced from MultiDivTime and BEAST analysis. These composite CrIs are found to contain the true time for >or=97% data sets. Our analyses also verify the usefulness of the common practice of interpreting the congruence of times inferred from different methods as a reflection of the accuracy of time estimates. Overall, our results show that simple strategies can be used to enhance our ability to estimate times and their CrIs when using the relaxed-clock methods.
| S-EPMC2877995 | biostudies-literature
Project description:Information about the geological timeframe during which animals radiated into their major subclades is crucial to understanding early animal ecology and evolution. Unfortunately, the pre-Cambrian fossil record is sparse and its interpretation controversial. Relaxed molecular-clock methods provide an alternative means of estimating the timing of cladogenesis deep in the metazoan tree of life. So far, thorough molecular clock studies focusing specifically on Metazoa as a whole have been based on relatively small datasets or incomplete representation of the main non-bilaterian lineages (such as sponges and ctenophores), which are fundamental for understanding early metazoan evolution. Here, we use a previously published phylogenomic dataset that includes a fair sampling of all relevant groups to estimate the timing of early animal evolution with Bayesian relaxed-clock methods. According to our results, all non-bilaterian phyla, as well as total-group Bilateria, evolved in an ancient radiation during a geologically relatively short time span, before the onset of long-term global glaciations ("Snowball Earth"; ~720-635?Ma). Importantly, this result appears robust to alterations of a number of important analytical variables, such as models of among-lineage rate variation and sets of fossil calibrations used.
| S-EPMC5472626 | biostudies-literature
Project description:Phylogenetic analyses are using increasingly larger data sets for estimating divergence times. With this increase in data sizes, the computation time required is becoming a bottleneck in evolutionary investigations. Our recent study of two relaxed-clock programs (BEAST and MultiDivTime [MDT]) showed their usefulness in time estimation; however, they place a significant computational time burden on biologists even for moderately small data sets. Here, we report speed and accuracy of another relaxed-clock program (MCMCTree, MC2T). We find it to be much faster than both MDT and BEAST while producing comparable time estimates. These results will encourage the analysis of larger data sets as well as the evaluation of the robustness of estimated times to changes in the model of evolutionary rates and clock calibrations.
| S-EPMC3203548 | biostudies-literature
Project description:Arecaceae is a species-rich clade of Arecales, while also being regarded as a morphologically diverse angiosperm family with numerous species having significant economic, medicinal, and ornamental value. Although in-depth studies focused on the chloroplast structure of Arecaceae, as well as inferring phylogenetic relationships using gene fragments, have been reported in recent years, a comprehensive analysis of the chloroplast structure of Arecaceae is still needed. Here we perform a comprehensive analysis of the structural features of the chloroplast genome of Arecaceae, compare the variability of gene sequences, infer phylogenetic relationships, estimate species divergence times, and reconstruct ancestral morphological traits. In this study, 74 chloroplast genomes of Arecaceae were obtained, covering five subfamilies. The results show that all chloroplast genomes possess a typical tetrad structure ranging in size between 153,806-160,122 bp, with a total of 130-137 genes, including 76-82 protein-coding genes, 29-32 tRNA genes, and 4 rRNA genes. Additionally, the total GC content was between 36.9-37.7%. Analysis of the SC/IR boundary indicated that the IR region underwent expansion or contraction. Phylogenetic relationships indicate that all five subfamilies in Arecaceae are monophyletic and that Ceroxyloideae and Arecoideae are sister groups (BS/PP = 100/1). The results of molecular dating indicate that the age of the crown group of Arecaceae is likely to be 96.60 [84.90-107.60] Ma, while the age of the stem group is 102.40 [93.44-111.17] Ma. Reconstruction of ancestral traits indicate that the ancestral characteristics of the family include monoecious plants, one seed, six stamens, and a smooth pericarp.
| S-EPMC9552784 | biostudies-literature
Project description:We present molecular dating analyses for land plants that incorporate 33 fossil calibrations, permit rates of molecular evolution to be uncorrelated across the tree, and take into account uncertainties in phylogenetic relationships and the fossil record. We attached a prior probability to each fossil-based minimum age, and explored the effects of relying on the first appearance of tricolpate pollen grains as a lower bound for the age of eudicots. Many of our divergence-time estimates for major clades coincide well with both the known fossil record and with previous estimates. However, our estimates for the origin of crown-clade angiosperms, which center on the Late Triassic, are considerably older than the unequivocal fossil record of flowering plants or than the molecular dates presented in recent studies. Nevertheless, we argue that our older estimates should be taken into account in studying the causes and consequences of the angiosperm radiation in relation to other major events, including the diversification of holometabolous insects. Although the methods used here do help to correct for lineage-specific heterogeneity in rates of molecular evolution (associated, for example, with evolutionary shifts in life history), we remain concerned that some such effects (e.g., the early radiation of herbaceous clades within angiosperms) may still be biasing our inferences.
| S-EPMC2851901 | biostudies-other