Fossil calibrations and molecular divergence time estimates in centrarchid fishes (Teleostei: Centrarchidae)

Molecular clock methods allow biologists to estimate divergence times, which in turn play an important role in comparative studies of many evolutionary processes. It is well known that molecular age estimates can be biased by heterogeneity in rates of molecular evolution, but less attention has been paid to the issue of potentially erroneous fossil calibrations. In this study we estimate the timing of diversification in Centrarchidae, an endemic major lineage of the diverse North American freshwater fish fauna, through a new approach to fossil calibration and molecular evolutionary model selection. Given a completely resolved multi-gene molecular phylogeny and a set of multiple fossil-inferred age estimates, we tested for potentially erroneous fossil calibrations using a recently developed fossil cross-validation. We also used fossil information to guide the selection of the optimal molecular evolutionary model with a new fossil jackknife method in a fossil-based model cross-validation. The centrarchid phylogeny resulted from a mixed-model Bayesian strategy that included 14 separate data partitions sampled from three mtDNA and four nuclear genes. Ten of the 31 interspecific nodes in the centrarchid phylogeny were assigned a minimal age estimate from the centrarchid fossil record. Our analyses identified four fossil dates that were inconsistent with the other fossils, and we removed them from the molecular dating analysis. Using fossil-based model cross-validation to determine the optimal smoothing value in penalized likelihood analysis, and six mutually consistent fossil calibrations, the age of the most recent common ancestor of Centrarchidae was 33.59 million years ago (mya). Penalized likelihood analyses of individual data partitions all converged on a very similar age estimate for this node, indicating that rate heterogeneity among data partitions is not confounding our analyses. These results place the origin of the centrarchid radiation at a time of major faunal turnover as the fossil record indicates that the most diverse lineages of the North American freshwater fish fauna originated at the Eocene-Oligocene boundary, approximately 34 mya. This time coincided with major global climate change from warm to cool temperatures and a signature of elevated lineage extinction and origination in the fossil record across the tree of life. Our analyses demonstrate the utility of fossil cross-validation to critically assess individual fossil calibration points, providing the ability to discriminate between consistent and inconsistent fossil age estimates that are used for calibrating molecular phylogenies.     

Evaluating fossil calibrations for dating phylogenies in light of rates of molecular evolution: a comparison of three approaches

Evolutionary and biogeographic studies increasingly rely on calibrated molecular clocks to date key events. Although there has been significant recent progress in development of the techniques used for molecular dating, many issues remain. In particular, controversies abound over the appropriate use and placement of fossils for calibrating molecular clocks. Several methods have been proposed for evaluating candidate fossils; however, few studies have compared the results obtained by different approaches. Moreover, no previous study has incorporated the effects of nucleotide saturation from different data types in the evaluation of candidate fossils. In order to address these issues, we compared three approaches for evaluating fossil calibrations: the single-fossil cross-validation method of Near, Meylan, and Shaffer (2005. Assessing concordance of fossil calibration points in molecular clock studies: an example using turtles. Am. Nat. 165:137-146), the empirical fossil coverage method of Marshall (2008. A simple method for bracketing absolute divergence times on molecular phylogenies using multiple fossil calibration points. Am. Nat. 171:726-742), and the Bayesian multicalibration method of Sanders and Lee (2007. Evaluating molecular clock calibrations using Bayesian analyses with soft and hard bounds. Biol. Lett. 3:275-279) and explicitly incorporate the effects of data type (nuclear vs. mitochondrial DNA) for identifying the most reliable or congruent fossil calibrations. We used advanced (Caenophidian) snakes as a case study; however, our results are applicable to any taxonomic group with multiple candidate fossils, provided appropriate taxon sampling and sufficient molecular sequence data are available. We found that data type strongly influenced which fossil calibrations were identified as outliers, regardless of which method was used. Despite the use of complex partitioned models of sequence evolution and multiple calibrations throughout the tree, saturation severely compressed basal branch lengths obtained from mitochondrial DNA compared with nuclear DNA. The effects of mitochondrial saturation were not ameliorated by analyzing a combined nuclear and mitochondrial data set. Although removing the third codon positions from the mitochondrial coding regions did not ameliorate saturation effects in the single-fossil cross-validations, it did in the Bayesian multicalibration analyses. Saturation significantly influenced the fossils that were selected as most reliable for all three methods evaluated. Our findings highlight the need to critically evaluate the fossils selected by data with different rates of nucleotide substitution and how data with different evolutionary rates affect the results of each method for evaluating fossils. Our empirical evaluation demonstrates that the advantages of using multiple independent fossil calibrations significantly outweigh any disadvantages.     

Assessing the quality of molecular divergence time estimates by fossil calibrations and fossil-based model selection

Estimates of species divergence times using DNA sequence data are playing an increasingly important role in studies of evolution, ecology and biogeography. Most work has centred on obtaining appropriate kinds of data and developing optimal estimation procedures, whereas somewhat less attention has focused on the calibration of divergences using fossils. Case studies with multiple fossil calibration points provide important opportunities to examine the divergence time estimation problem in new ways. We discuss two cross-validation procedures that address different aspects of inference in divergence time estimation. 'Fossil cross-validation' is a procedure used to identify the impact of different individual calibrations on overall estimation. This can identify fossils that have an exceptionally large error effect and may warrant further scrutiny. 'Fossil-based model cross-validation' is an entirely different procedure that uses fossils to identify the optimal model of molecular evolution in the context of rate smoothing or other inference methods. Both procedures were applied to two recent studies: an analysis of monocot angiosperms with eight fossil calibrations and an analysis of placental mammals with nine fossil calibrations. In each case, fossil calibrations could be ranked from most to least influential, and in one of the two studies, the fossils provided decisive evidence about the optimal molecular evolutionary model.     

Integrating fossil preservation biases in the selection of calibrations for molecular divergence time estimation

The selection of fossil data to use as calibration age priors in molecular divergence time estimates inherently links neontological methods with paleontological theory. However, few neontological studies have taken into account the possibility of a taphonomic bias in the fossil record when developing approaches to fossil calibration selection. The Sppil-Rongis effect may bias the first appearance of a lineage toward the recent causing most objective calibration selection approaches to erroneously exclude appropriate calibrations or to incorporate multiple calibrations that are too young to accurately represent the divergence times of target lineages. Using turtles as a case study, we develop a Bayesian extension to the fossil selection approach developed by Marshall (2008. A simple method for bracketing absolute divergence times on molecular phylogenies using multiple fossil calibrations points. Am. Nat. 171:726-742) that takes into account this taphonomic bias. Our method has the advantage of identifying calibrations that may bias age estimates to be too recent while incorporating uncertainty in phylogenetic parameter estimates such as tree topology and branch lengths. Additionally, this method is easily adapted to assess the consistency of potential calibrations to any one calibration in the candidate pool.     

Calibrating the Tree of Life: fossils, molecules and evolutionary timescales

Background

Molecular dating has gained ever-increasing interest since the molecular clock hypothesis was proposed in the 1960s. Molecular dating provides detailed temporal frameworks for divergence events in phylogenetic trees, allowing diverse evolutionary questions to be addressed. The key aspect of the molecular clock hypothesis, namely that differences in DNA or protein sequence between two species are proportional to the time elapsed since they diverged, was soon shown to be untenable. Other approaches were proposed to take into account rate heterogeneity among lineages, but the calibration process, by which relative times are transformed into absolute ages, has received little attention until recently. New methods have now been proposed to resolve potential sources of error associated with the calibration of phylogenetic trees, particularly those involving use of the fossil record.

Scope and Conclusions

The use of the fossil record as a source of independent information in the calibration process is the main focus of this paper; other sources of calibration information are also discussed. Particularly error-prone aspects of fossil calibration are identified, such as fossil dating, the phylogenetic placement of the fossil and the incompleteness of the fossil record. Methods proposed to tackle one or more of these potential error sources are discussed (e.g. fossil cross-validation, prior distribution of calibration points and confidence intervals on the fossil record). In conclusion, the fossil record remains the most reliable source of information for the calibration of phylogenetic trees, although associated assumptions and potential bias must be taken into account.     

Bayesian estimation of species divergence times under a molecular clock using multiple fossil calibrations with soft bounds

We implement a Bayesian Markov chain Monte Carlo algorithm for estimating species divergence times that uses heterogeneous data from multiple gene loci and accommodates multiple fossil calibration nodes. A birth-death process with species sampling is used to specify a prior for divergence times, which allows easy assessment of the effects of that prior on posterior time estimates. We propose a new approach for specifying calibration points on the phylogeny, which allows the use of arbitrary and flexible statistical distributions to describe uncertainties in fossil dates. In particular, we use soft bounds, so that the probability that the true divergence time is outside the bounds is small but nonzero. A strict molecular clock is assumed in the current implementation, although this assumption may be relaxed. We apply our new algorithm to two data sets concerning divergences of several primate species, to examine the effects of the substitution model and of the prior for divergence times on Bayesian time estimation. We also conduct computer simulation to examine the differences between soft and hard bounds. We demonstrate that divergence time estimation is intrinsically hampered by uncertainties in fossil calibrations, and the error in Bayesian time estimates will not go to zero with increased amounts of sequence data. Our analyses of both real and simulated data demonstrate potentially large differences between divergence time estimates obtained using soft versus hard bounds and a general superiority of soft bounds. Our main findings are as follows. (1) When the fossils are consistent with each other and with the molecular data, and the posterior time estimates are well within the prior bounds, soft and hard bounds produce similar results. (2) When the fossils are in conflict with each other or with the molecules, soft and hard bounds behave very differently; soft bounds allow sequence data to correct poor calibrations, while poor hard bounds are impossible to overcome by any amount of data. (3) Soft bounds eliminate the need for "safe" but unrealistically high upper bounds, which may bias posterior time estimates. (4) Soft bounds allow more reliable assessment of estimation errors, while hard bounds generate misleadingly high precisions when fossils and molecules are in conflict.     

Can fast early rates reconcile molecular dates with the Cambrian explosion?

Molecular dates consistently place the divergence of major metazoan lineages in the Precambrian, leading to the suggestion that the 'Cambrian explosion' is an artefact of preservation which left earlier forms unrecorded in the fossil record. While criticisms of molecular analyses for failing to deal with variation in the rate of molecular evolution adequately have been countered by analyses which allow both site-to-site and lineage-specific rate variation, no analysis to date has allowed the rates to vary temporally. If the rates of molecular evolution were much higher early in the metazoan radiation, molecular dates could consistently overestimate the divergence times of lineages. Here, we use a new method which uses multiple calibration dates and an empirically determined range of possible substitution rates to place bounds on the basal date of divergence of lineages in order to ask whether faster rates of molecular evolution early in the metazoan radiation could possibly account for the discrepancy between molecular and palaeontological date estimates. We find that allowing basal (interphylum) lineages the fastest observed substitution rate brings the minimum possible divergence date (586 million years ago) to the Vendian period, just before the first multicellular animal fossils, but excludes divergence of the major metazoan lineages in a Cambrian explosion.     

Calibration choice, rate smoothing, and the pattern of tetrapod diversification according to the long nuclear gene RAG-1

A phylogeny of tetrapods is inferred from nearly complete sequences of the nuclear RAG-1 gene sampled across 88 taxa encompassing all major clades, analyzed via parsimony and Bayesian methods. The phylogeny provides support for Lissamphibia, Theria, Lepidosauria, a turtle-archosaur clade, as well as most traditionally accepted groupings. This tree allows simultaneous molecular clock dating for all tetrapod groups using a set of well-corroborated calibrations. Relaxed clock (PLRS) methods, using the amniote = 315 Mya (million years ago) calibration or a set of consistent calibrations, recovers reasonable divergence dates for most groups. However, the analysis systematically underestimates divergence dates within archosaurs. The bird-crocodile split, robustly documented in the fossil record as being around approximately 245 Mya, is estimated at only approximately 190 Mya, and dates for other divergences within archosaurs are similarly underestimated. Archosaurs, and particulary turtles have slow apparent rates possibly confounding rate modeling, and inclusion of calibrations within archosaurs (despite their high deviances) not only improves divergence estimates within archosaurs, but also across other groups. Notably, the monotreme-therian split ( approximately 210 Mya) matches the fossil record; the squamate radiation ( approximately 190 Mya) is younger than suggested by some recent molecular studies and inconsistent with identification of approximately 220 and approximately 165 Myo (million-year-old) fossils as acrodont iguanians and approximately 95 Myo fossils colubroid snakes; the bird-lizard (reptile) split is considerably older than fossil estimates (< or = 285 Mya); and Sphenodon is a remarkable phylogenetic relic, being the sole survivor of a lineage more than a quarter of a billion years old. Comparison with other molecular clock studies of tetrapod divergences suggests that the common practice of enforcing most calibrations as minima, with a single liberal maximal constraint, will systematically overestimate divergence dates. Similarly, saturation of mitochondrial DNA sequences, and the resultant greater compression of basal branches means that using only external deep calibrations will also lead to inflated age estimates within the focal ingroup.     

Bayesian relaxed clock estimation of divergence times in foraminifera

Accurate and precise estimation of divergence times during the Neo-Proterozoic is necessary to understand the speciation dynamic of early Eukaryotes. However such deep divergences are difficult to date, as the molecular clock is seriously violated. Recent improvements in Bayesian molecular dating techniques allow the relaxation of the molecular clock hypothesis as well as incorporation of multiple and flexible fossil calibrations. Divergence times can then be estimated even when the evolutionary rate varies among lineages and even when the fossil calibrations involve substantial uncertainties. In this paper, we used a Bayesian method to estimate divergence times in Foraminifera, a group of unicellular eukaryotes, known for their excellent fossil record but also for the high evolutionary rates of their genomes. Based on multigene data we reconstructed the phylogeny of Foraminifera and dated their origin and the major radiation events. Our estimates suggest that Foraminifera emerged during the Cryogenian (650-920 Ma, Neo-Proterozoic), with a mean time around 770 Ma, about 220 Myr before the first appearance of reliable foraminiferal fossils in sediments (545 Ma). Most dates are in agreement with the fossil record, but in general our results suggest earlier origins of foraminiferal orders. We found that the posterior time estimates were robust to specifications of the prior. Our results highlight inter-species variations of evolutionary rates in Foraminifera. Their effect was partially overcome by using the partitioned Bayesian analysis to accommodate rate heterogeneity among data partitions and using the relaxed molecular clock to account for changing evolutionary rates. However, more coding genes appear necessary to obtain more precise estimates of divergence times and to resolve the conflicts between fossil and molecular date estimates.     

Are pollen fossils useful for calibrating relaxed molecular clock dating of phylogenies? A comparative study using Myrtaceae

The identification and application of reliable fossil calibrations represents a key component of many molecular studies of evolutionary timescales. In studies of plants, most paleontological calibrations are associated with macrofossils. However, the pollen record can also inform age calibrations if fossils matching extant pollen groups are found. Recent work has shown that pollen of the myrtle family, Myrtaceae, can be classified into a number of morphological groups that are synapomorphic with molecular groups. By assembling a data matrix of pollen morphological characters from extant and fossil Myrtaceae, we were able to measure the fit of 26 pollen fossils to a molecular phylogenetic tree using parsimony optimisation of characters. We identified eight Myrtaceidites fossils as appropriate for calibration based on the most parsimonious placements of these fossils on the tree. These fossils were used to inform age constraints in a Bayesian phylogenetic analysis of a sequence alignment comprising two sequences from the chloroplast genome (matK and ndhF) and one nuclear locus (ITS), sampled from 106 taxa representing 80 genera. Three additional analyses were calibrated by placing pollen fossils using geographic and morphological information (eight calibrations), macrofossils (five calibrations), and macrofossils and pollen fossils in combination (12 calibrations). The addition of new fossil pollen calibrations led to older crown ages than have previously been found for tribes such as Eucalypteae and Myrteae. Estimates of rate variation among lineages were affected by the choice of calibrations, suggesting that the use of multiple calibrations can improve estimates of rate heterogeneity among lineages. This study illustrates the potential of including pollen-based calibrations in molecular studies of divergence times.     

Dated Phylogenies of the Sister Genera Macaranga and Mallotus (Euphorbiaceae): Congruence in Historical Biogeographic Patterns?

Molecular phylogenies and estimates of divergence times within the sister genera Macaranga and Mallotus were estimated using Bayesian relaxed clock analyses of two generic data sets, one per genus. Both data sets were based on different molecular markers and largely different samples. Per genus three calibration points were utilised. The basal calibration point (crown node of all taxa used) was taken from literature and used for both taxa. The other three calibrations were based on fossils of which two were used per genus. We compared patterns of dispersal and diversification in Macaranga and Mallotus using ancestral area reconstruction in RASP (S-DIVA option) and contrasted our results with biogeographical and geological records to assess accuracy of inferred age estimates. A check of the fossil calibration point showed that the Japanese fossil, used for dating the divergence of Mallotus, probably had to be attached to a lower node, the stem node of all pioneer species, but even then the divergence time was still younger than the estimated age of the fossil. The African (only used in the Macaranga data set) and New Zealand fossils (used for both genera) seemed reliably placed. Our results are in line with existing geological data and the presence of stepping stones that provided dispersal pathways from Borneo to New Guinea-Australia, from Borneo to mainland Asia and additionally at least once to Africa and Madagascar via land and back to India via Indian Ocean island chains. The two genera show congruence in dispersal patterns, which corroborate divergence time estimates, although the overall mode and tempo of dispersal and diversification differ significantly as shown by distribution patterns of extant species.     

A simple method for bracketing absolute divergence times on molecular phylogenies using multiple fossil calibration points

A central challenge facing the temporal calibration of molecular phylogenies is finding a quantitative method for estimating maximum age constraints on lineage divergence times. Here, I provide such a method. This method requires an ultrametric tree generated without reference to the fossil record. Exploiting the fact that the relative branch lengths on the ultrametric tree are proportional to time, this method identifies the lineage with the greatest proportion of its true temporal range covered by the fossil record. The oldest fossil of this calibration lineage is used as the minimum age constraint. The maximum age constraint is obtained by adding a confidence interval onto the end point of the calibration lineage, thus making it possible to bracket the true divergence times of all lineages on the tree. The approach can also identify fossils that have been grossly misdated or misassigned to the phylogeny. The method assumes that the relative branch lengths on the ultrametric tree are accurate and that fossilization is random. The effect of violations of these assumptions is assessed. This method is simple to use and is illustrated with a reanalysis of Near et al.'s turtle data.     

Comparison of likelihood and Bayesian methods for estimating divergence times using multiple gene Loci and calibration points, with application to a radiation of cute-looking mouse lemur species

Divergence time and substitution rate are seriously confounded in phylogenetic analysis, making it difficult to estimate divergence times when the molecular clock (rate constancy among lineages) is violated. This problem can be alleviated to some extent by analyzing multiple gene loci simultaneously and by using multiple calibration points. While different genes may have different patterns of evolutionary rate change, they share the same divergence times. Indeed, the fact that each gene may violate the molecular clock differently leads to the advantage of simultaneous analysis of multiple loci. Multiple calibration points provide the means for characterizing the local evolutionary rates on the phylogeny. In this paper, we extend previous likelihood models of local molecular clock for estimating species divergence times to accommodate multiple calibration points and multiple genes. Heterogeneity among different genes in evolutionary rate and in substitution process is accounted for by the models. We apply the likelihood models to analyze two mitochondrial protein-coding genes, cytochrome oxidase II and cytochrome b, to estimate divergence times of Malagasy mouse lemurs and related outgroups. The likelihood method is compared with the Bayes method of Thorne et al. (1998, Mol. Biol. Evol. 15:1647-1657), which uses a probabilistic model to describe the change in evolutionary rate over time and uses the Markov chain Monte Carlo procedure to derive the posterior distribution of rates and times. Our likelihood implementation has the drawbacks of failing to accommodate uncertainties in fossil calibrations and of requiring the researcher to classify branches on the tree into different rate groups. Both problems are avoided in the Bayes method. Despite the differences in the two methods, however, data partitions and model assumptions had the greatest impact on date estimation. The three codon positions have very different substitution rates and evolutionary dynamics, and assumptions in the substitution model affect date estimation in both likelihood and Bayes analyses. The results demonstrate that the separate analysis is unreliable, with dates variable among codon positions and between methods, and that the combined analysis is much more reliable. When the three codon positions were analyzed simultaneously under the most realistic models using all available calibration information, the two methods produced similar results. The divergence of the mouse lemurs is dated to be around 7-10 million years ago, indicating a surprisingly early species radiation for such a morphologically uniform group of primates.     

Fossil calibration of molecular divergence infers a moderate mutation rate and recent radiations for pinus

Silent mutation rate estimates for Pinus vary 50-fold, ranging from angiosperm-like to among the slowest reported for plants. These differences either reflect extraordinary genomic processes or inconsistent fossil calibration, and they have important consequences for population and biogeographical inferences. Here we estimate mutation rates from 4 Pinus species that represent the major lineages using 11 nuclear and 4 chloroplast loci. Calibration was tested at the divergence of Pinus subgenera with the oldest leaf fossil from subg. Strobus (Eocene; 45 MYA) or a recently published subg. Strobus wood fossil (Cretaceous; 85 MYA). These calibrations place the origin of Pinus 190-102 MYA and give absolute silent rate estimates of 0.70-1.31x10(-9) and 0.22-0.42x10(-9).site-1.year-1 for the nuclear and chloroplast genomes, respectively. These rates are approximately 4- to 20-fold slower than angiosperms, but unlike many previous estimates, they are more consistent with the high per-generation deleterious mutation rates observed in pines. Chronograms from nuclear and chloroplast genomes show that the divergence of subgenera accounts for about half of the time since Pinus diverged from Picea, with subsequent radiations occurring more recently. By extending the sampling to encompass the phylogenetic diversity of Pinus, we predict that most extant subsections diverged during the Miocene. Moreover, subsect. Australes, Ponderosae, and Contortae, containing over 50 extant species, radiated within a 5 Myr time span starting as recently as 18 MYA. An Eocene divergence of pine subgenera (using leaf fossils) does not conflict with fossil-based estimates of the Pinus-Picea split, but a Cretaceous divergence using wood fossils accommodates Oligocene fossils that may represent modern subsections. Because homoplasy and polarity of character states have not been tested for fossil pine assignments, the choice of fossil and calibration node represents a significant source of uncertainty. Based on several lines of evidence (including agreement with ages inferred using calibrations outside of Pinus), we conclude that the 85 MYA calibration at the divergence of pine subgenera provides a reasonable lower bound and that further refinements in age and mutation rate estimates will require a synthetic examination of pine fossil history.     

Estimating divergence dates from molecular sequences

The ability to date the time of divergence between lineages using molecular data provides the opportunity to answer many important questions in evolutionary biology. However, molecular dating techniques have previously been criticized for failing to adequately account for variation in the rate of molecular evolution. We present a maximum- likelihood approach to estimating divergence times that deals explicitly with the problem of rate variation. This method has many advantages over previous approaches including the following: (1) a rate constancy test excludes data for which rate heterogeneity is detected; (2) date estimates are generated with confidence intervals that allow the explicit testing of hypotheses regarding divergence times; and (3) a range of sequences and fossil dates are used, removing the reliance on a single calculated calibration rate. We present tests of the accuracy of our method, which show it to be robust to the effects of some modes of rate variation. In addition, we test the effect of substitution model and length of sequence on the accuracy of the dating technique. We believe that the method presented here offers solutions to many of the problems facing molecular dating and provides a platform for future improvements to such analyses.      

Assessing calibration uncertainty in molecular dating: the assignment of fossils to alternative calibration points

Although recent methodological advances have allowed the incorporation of rate variation in molecular dating analyses, the calibration procedure, performed mainly through fossils, remains resistant to improvements. One source of uncertainty pertains to the assignment of fossils to specific nodes in a phylogeny, especially when alternative possibilities exist that can be equally justified on morphological grounds. Here we expand on a recently developed fossil cross-validation method to evaluate whether alternative nodal assignments of multiple fossils produce calibration sets that differ in their internal consistency. We use an enlarged Crypteroniaceae-centered phylogeny of Myrtales, six fossils, and 72 combinations of calibration points, termed calibration sets, to identify (i) the fossil assignments that produce the most internally consistent calibration sets and (ii) the mean ages, derived from these calibration sets, for the split of the Southeast Asian Crypteroniaceae from their West Gondwanan sister clade (node X). We found that a correlation exists between s values, devised to measure the consistency among the calibration points of a calibration set (Near and Sanderson, 2004), and nodal distances among calibration points. By ranking all sets according to the percent deviation of s from the regression line with nodal distance, we identified the sets with the highest level of corrected calibration-set consistency. These sets generated lower standard deviations associated with the ages of node X than sets characterized by lower corrected consistency. The three calibration sets with the highest corrected consistencies produced mean age estimates for node X of 79.70, 79.14, and 78.15 My. These timeframes are most compatible with the hypothesis that the Crypteroniaceae stem lineage dispersed from Africa to the Deccan plate as it drifted northward during the Late Cretaceous.     

Dating the arthropod tree based on large-scale transcriptome data

Molecular sequences do not only allow the reconstruction of phylogenetic relationships among species, but also provide information on the approximate divergence times. Whereas the fossil record dates the origin of most multicellular animal phyla during the Cambrian explosion less than 540 million years ago(mya), molecular clock calculations usually suggest much older dates. Here we used a large multiple sequence alignment derived from Expressed Sequence Tags and genomes comprising 129genes (37,476 amino acid positions) and 117 taxa, including 101 arthropods. We obtained consistent divergence time estimates applying relaxed Bayesian clock models with different priors and multiple calibration points. While the influence of substitution rates, missing data, and model priors were negligible, the clock model had significant effect. A log-normal autocorrelated model was selected on basis of cross-validation. We calculated that arthropods emerged ~600 mya. Onychophorans (velvet worms) and euarthropods split ~590 mya, Pancrustacea and Myriochelata ~560 mya, Myriapoda and Chelicerata ~555 mya, and 'Crustacea' and Hexapoda ~510 mya. Endopterygote insects appeared ~390 mya. These dates are considerably younger than most previous molecular clock estimates and in better agreement with the fossil record. Nevertheless, a Precambrian origin of arthropods and other metazoan phyla is still supported. Our results also demonstrate the applicability of large datasets of random nuclear sequences for approximating the timing of multicellular animal evolution.     

Calibrating molecular estimates of substitution rates and divergence times in birds

Many evolutionary studies of birds rely on the estimation of molecular divergence times and substitution rates. In order to perform such analyses, it is necessary to incorporate some form of calibration information: a known substitution rate, radiometric ages of heterochronous sequences, or inferred ages of lineage splitting events. All three of these techniques have been employed in avian molecular studies, but their usage has not been entirely satisfactory. For example, the 'traditional' avian mitochondrial substitution rate of 2% per million years is frequently adopted without acknowledgement of the associated uncertainty. Similarly, fossil and biogeographic information is almost always converted into an errorless calibration point. In both cases, the resulting estimates of divergence times and substitution rates will be artificially precise, which has a considerable impact on hypothesis testing. In addition, using such a simplistic approach to calibration discards much of the information offered by the fossil record. A number of more sophisticated calibration methods have recently been introduced, culminating in the development of probability distribution-based calibrations. In this article, I discuss the use of this new class of methods and offer guidelines for choosing a calibration technique.     

Phylogenetics of Perissodactyla and Tests of the Molecular Clock

Two mitochondrial genes, the protein-coding cytochrome c oxidase subunit II (COII) gene and a portion of the 12S rRNA gene, were used for phylogenetic investigation of the mammalian order Perissodactyla. The primary objective of the study was to utilize the extensive fossil record of perissodactyls for calibrating molecular clocks and comparing estimates of divergence times using both genes and two fossil calibration points. Secondary objectives included clarification of previously unresolved relationships within Tapiridae and comparison of the results of separate and combined analyses of two genes. Analyses included several perissodactyl lineages representing all three families (Tapiridae, Equidae, and Rhinocerotidae), most extant genera, all four species of tapirs, two to four species of rhinoceros, and two species of Equus. The application of a relatively recent fossil calibration point and a relatively ancient calibration point produced greatly different estimates of evolutionary rates and divergence times for both genes, even though a relative rates test did not find significant rate differences among taxa. A likelihood-ratio test, however, rejected a molecular clock for both genes. Neither calibration point produced estimates of divergence times consistent with paleontological evidence over a range of perissodactyl radiations. The combined analysis of both genes produces a well-resolved phylogeny with Perissodactyla that conforms to traditional views of interfamilial relationships and supports monophyly of neotropical tapirs. Combining the data sets increases support for most nodes but decreases the support for a neotropical tapir clade because the COII and 12S rRNA data sets are in conflict for tapir relationships. Received: 6 January 1999 / Accepted: 2 August 1999