Fast and slow implementations of relaxed-clock methods show similar patterns of accuracy in estimating divergence times.
Ontology highlight
ABSTRACT:
SUBMITTER: Battistuzzi FU
PROVIDER: S-EPMC3203548 | biostudies-literature | 2011 Sep
REPOSITORIES: biostudies-literature
Publications
Fast and slow implementations of relaxed-clock methods show similar patterns of accuracy in estimating divergence times.
Battistuzzi Fabia U FU Billing-Ross Paul P Paliwal Aditya A Kumar Sudhir S
Molecular biology and evolution 20110415 9
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 relax ...[more]
Similar Datasets
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:Accurately dating when the first bilaterally symmetrical animals arose is crucial to our understanding of early animal evolution. The earliest unequivocally bilaterian fossils are approximately 555 million years old. In contrast, molecular-clock analyses calibrated by using the fossil record of vertebrates estimate that vertebrates split from dipterans (Drosophila) approximately 900 million years ago (Ma). Nonetheless, comparative genomic analyses suggest that a significant rate difference exists between vertebrates and dipterans, because the percentage difference between the genomes of mosquito and fly is greater than between fish and mouse, even though the vertebrate divergence is almost twice that of the dipteran. Here we show that the dipteran rate of molecular evolution is similar to other invertebrate taxa (echinoderms and bivalve molluscs) but not to vertebrates, which significantly decreased their rate of molecular evolution with respect to invertebrates. Using a data set consisting of the concatenation of seven different amino acid sequences from 23 ingroup taxa (giving a total of 11 different invertebrate calibration points scattered throughout the bilaterian tree and across the Phanerozoic), we estimate that the last common ancestor of bilaterians arose somewhere between 573 and 656 Ma, depending on the value assigned to the parameter scaling molecular substitution rate heterogeneity. These results are in accord with the known fossil record and support the view that the Cambrian explosion reflects, in part, the diversification of bilaterian phyla.
| S-EPMC404080 | biostudies-literature
Project description:Fast and slow earthquakes are two modes of energy release by the slip in tectonic fault rupture. Although fast and slow slips were observed in the laboratory stick-slip experiments, due to the sampling rate limitation, the details of the fault thickness variation were poorly understood. Especially, why a single fault would show different modes of slip remains elusive. Herein, we report on ring shear experiments with an ultrahigh sampling rate (10 MHz) that illuminate the different physical processes between fast and slow slip events. We show that the duration of slips ranged from dozens to hundreds of milliseconds. Fast slip events are characterized by continuous large-amplitude AE (acoustic emission) and somewhat intricate variation of the sample thickness: A short compaction pulse during the rapid release of stress is followed by dilation and vibrations of the sample thickness. As the slip ends, the thickness of the sample first recovers by slow compaction and then dilates again before nucleation of the following slip event. In contrast, during slow slip events, the shear stress reduction is accompanied by intermittent bursts of low-amplitude AE and sample dilation. We observed the detailed thickness variation during slips and found that dilation occurs during both fast and slow slips, which is consistent with natural observations of coseismic dilatation. This study may be used to reveal the mechanism of fault slips during fast and slow earthquakes, which explain the potential effect of fast and slow slips on stress redistribution and structural rearrangement in faults.
| S-EPMC10666130 | biostudies-literature
Project description:Molecular dating of species divergences has become an important means to add a temporal dimension to the Tree of Life. Increasingly larger datasets encompassing greater taxonomic diversity are becoming available to generate molecular timetrees by using sophisticated methods that model rate variation among lineages. However, the practical application of these methods is challenging because of the exorbitant calculation times required by current methods for contemporary data sizes, the difficulty in correctly modeling the rate heterogeneity in highly diverse taxonomic groups, and the lack of reliable clock calibrations and their uncertainty distributions for most groups of species. Here, we present a method that estimates relative times of divergences for all branching points (nodes) in very large phylogenetic trees without assuming a specific model for lineage rate variation or specifying any clock calibrations. The method (RelTime) performed better than existing methods when applied to very large computer simulated datasets where evolutionary rates were varied extensively among lineages by following autocorrelated and uncorrelated models. On average, RelTime completed calculations 1,000 times faster than the fastest Bayesian method, with even greater speed difference for larger number of sequences. This speed and accuracy will enable molecular dating analysis of very large datasets. Relative time estimates will be useful for determining the relative ordering and spacing of speciation events, identifying lineages with significantly slower or faster evolutionary rates, diagnosing the effect of selected calibrations on absolute divergence times, and estimating absolute times of divergence when highly reliable calibration points are available.
| S-EPMC3511068 | biostudies-literature
Project description:Understanding causes of biological diversity may be greatly enhanced by knowledge of divergence times. Strict and relaxed clock models are used in Bayesian estimation of divergence times. We examined whether: i) strict clock models are generally more appropriate in shallow phylogenies where rate variation is expected to be low, ii) the likelihood ratio test of the clock (LRT) reliably informs which model is appropriate for dating divergence times. Strict and relaxed models were used to analyse sequences simulated under different levels of rate variation. Published shallow phylogenies (Black bass, Primate-sucking lice, Podarcis lizards, Gallotiinae lizards, and Caprinae mammals) were also analysed to determine natural levels of rate variation relative to the performance of the different models.Strict clock analyses performed well on data simulated under the independent rates model when the standard deviation of log rate on branches, ?, was low (? 0.1), but were inappropriate when ?>0.1 (95% of rates fall within 0.0082-0.0121 subs/site/Ma when ? = 0.1, for a mean rate of 0.01). The independent rates relaxed clock model performed well at all levels of rate variation, although posterior intervals on times were significantly wider than for the strict clock. The strict clock is therefore superior when rate variation is low. The performance of a correlated rates relaxed clock model was similar to the strict clock. Increased numbers of independent loci led to slightly narrower posteriors under the relaxed clock while older root ages provided proportionately narrower posteriors. The LRT had low power for ? = 0.01-0.1, but high power for ? = 0.5-2.0. Posterior means of ?2 were useful for assessing rate variation in published datasets. Estimates of natural levels of rate variation ranged from 0.05-3.38 for different partitions. Differences in divergence times between relaxed and strict clock analyses were greater in two datasets with higher ?2 for one or more partitions, supporting the simulation results.The strict clock can be superior for trees with shallow roots because of low levels of rate variation between branches. The LRT allows robust assessment of suitability of the clock model as does examination of posteriors on ?2.
| S-EPMC3205074 | biostudies-literature
Project description:Accurate estimation of divergence times of soil bacteria that form nitrogen-fixing associations with most leguminous plants is challenging because of a limited fossil record and complexities associated with molecular clocks and phylogenetic diversity of root nodule bacteria, collectively called rhizobia. To overcome the lack of fossil record in bacteria, divergence times of host legumes were used to calibrate molecular clocks and perform phylogenetic analyses in rhizobia. The 16S rRNA gene and intergenic spacer region remain among the favored molecular markers to reconstruct the timescale of rhizobia. We evaluate the performance of the random local clock model and the classical uncorrelated lognormal relaxed clock model, in combination with four tree models (coalescent constant size, birth-death, birth-death incomplete sampling, and Yule processes) on rhizobial divergence time estimates. Bayes factor tests based on the marginal likelihoods estimated from the stepping-stone sampling analyses strongly favored the random local clock model in combination with Yule process. Our results on the divergence time estimation from 16S rRNA gene and intergenic spacer region sequences are compatible with age estimates based on the conserved core genes but significantly older than those obtained from symbiotic genes, such as nodIJ genes. This difference may be due to the accelerated evolutionary rates of symbiotic genes compared to those of other genomic regions not directly implicated in nodulation processes.
| S-EPMC4856229 | biostudies-literature
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:BackgroundMatrices of morphological characters are frequently used for dating species divergence times in systematics. In some studies, morphological and molecular character data from living taxa are combined, whereas others use morphological characters from extinct taxa as well. We investigated whether morphological data produce time estimates that are concordant with molecular data. If true, it will justify the use of morphological characters alongside molecular data in divergence time inference.ResultsWe systematically analyzed three empirical datasets from different species groups to test the concordance of species divergence dates inferred using molecular and discrete morphological data from extant taxa as test cases. We found a high correlation between their divergence time estimates, despite a poor linear relationship between branch lengths for morphological and molecular data mapped onto the same phylogeny. This was because node-to-tip distances showed a much higher correlation than branch lengths due to an averaging effect over multiple branches. We found that nodes with a large number of taxa often benefit from such averaging. However, considerable discordance between time estimates from molecules and morphology may still occur as some intermediate nodes may show large time differences between these two types of data.ConclusionsOur findings suggest that node- and tip-calibration approaches may be better suited for nodes with many taxa. Nevertheless, we highlight the importance of evaluating the concordance of intrinsic time structure in morphological and molecular data before any dating analysis using combined datasets.
| S-EPMC8117668 | biostudies-literature
Project description:Higher N170 amplitudes to words and to faces were recently reported for faster readers of German. Since the shallow German orthography allows phonological recoding of single letters, the reported speed advantages might have their origin in especially well-developed visual processing skills of faster readers. In contrast to German, adult readers of Hebrew are forced to process letter chunks up to whole words. This dependence on more complex visual processing might have created ceiling effects for this skill. Therefore, the current study examined whether also in the deep Hebrew orthography visual processing skills as reflected by N170 amplitudes explain reading speed differences. Forty university students, native speakers of Hebrew without reading impairments, accomplished a lexical decision task (i.e., deciding whether a visually presented stimulus represents a real or a pseudo word) and a face decision task (i.e., deciding whether a face was presented complete or with missing facial features) while their electroencephalogram was recorded from 64 scalp positions. In both tasks stronger event related potentials (ERPs) were observed for faster readers in time windows at about 200 ms. Unlike in previous studies, ERP waveforms in relevant time windows did not correspond to N170 scalp topographies. The results support the notion of visual processing ability as an orthography independent marker of reading proficiency, which advances our understanding about regular and impaired reading development.
| S-EPMC4117504 | biostudies-literature
Project description:Background and aimsAlthough dioecy, which characterizes only 6 % of angiosperm species, has been considered an evolutionary dead end, recent studies have demonstrated that this is not necessarily the case. Moraceae (40 genera, 1100 spp., including Ficus, 750 spp.) are particularly diverse in breeding systems (including monoecy, gynodioecy, androdioecy and dioecy) and thus represent a model clade to study macroevolution of dioecy.MethodsAncestral breeding systems of Ficus and Moraceae were inferred. To do so, a new dated phylogenetic tree of Ficus and Moraceae was first reconstructed by combining a revised 12-fossil calibration set and a densely sampled molecular data set of eight markers and 320 species. Breeding system evolution was then reconstructed using both parsimony and model-based (maximum likelihood and Bayesian) approaches with this new time scale.Key resultsThe crown group ages of Ficus and Moraceae were estimated in the Eocene (40.6-55.9 Ma) and Late Cretaceous (73.2-84.7 Ma), respectively. Strong support was found for ancestral dioecy in Moraceae. Although the ancestral state of Ficus remained particularly sensitive to model selection, the results show that monoecy and gynodioecy evolved from dioecy in Moraceae, and suggest that gynodioecy probably evolved from monoecy in Ficus.ConclusionsDioecy was found not to be an evolutionary dead end in Moraceae. This study provides a new time scale for the phylogeny and a new framework of breeding system evolution in Ficus and Moraceae.
| S-EPMC6344110 | biostudies-literature