Temporal patterns of fruit fly (Drosophila) evolution revealed by mutation clocks

Drosophila melanogaster has been a canonical model organism to study genetics, development, behavior, physiology, evolution, and population genetics for nearly a century. Despite this emphasis and the completion of its nuclear genome sequence, the timing of major speciation events leading to the origin of this fruit fly remain elusive because of the paucity of extensive fossil records and biogeographic data. Use of molecular clocks as an alternative has been fraught with non-clock-like accumulation of nucleotide and amino-acid substitutions. Here we present a novel methodology in which genomic mutation distances are used to overcome these limitations and to make use of all available gene sequence data for constructing a fruit fly molecular time scale. Our analysis of 2977 pairwise sequence comparisons from 176 nuclear genes reveals a long-term fruit fly mutation clock ticking at a rate of 11.1 mutations per kilobase pair per Myr. Genomic mutation clock-based timings of the landmark speciation events leading to the evolution of D. melanogaster show that it shared most recent common ancestry 5.4 MYA with D. simulans, 12.6 MYA with D. erecta+D. orena, 12.8 MYA with D. yakuba+D. teisseri, 35.6 MYA with the takahashii subgroup, 41.3 MYA with the montium subgroup, 44.2 MYA with the ananassae subgroup, 54.9 MYA with the obscura group, 62.2 MYA with the willistoni group, and 62.9 MYA with the subgenus Drosophila. These and other estimates are compatible with those known from limited biogeographic and fossil records. The inferred temporal pattern of fruit fly evolution shows correspondence with the cooling patterns of paleoclimate changes and habitat fragmentation in the Cenozoic.     

Estimation of primate speciation dates using local molecular clocks

Protein-coding genes of the mitochondrial genomes from 31 mammalian species were analyzed to estimate the speciation dates within primates and also between rats and mice. Three calibration points were used based on paleontological data: one at 20-25 MYA for the hominoid/cercopithecoid divergence, one at 53-57 MYA for the cetacean/artiodactyl divergence, and the third at 110-130 MYA for the metatherian/eutherian divergence. Both the nucleotide and the amino acid sequences were analyzed, producing conflicting results. The global molecular clock was clearly violated for both the nucleotide and the amino acid data. Models of local clocks were implemented using maximum likelihood, allowing different evolutionary rates for some lineages while assuming rate constancy in others. Surprisingly, the highly divergent third codon positions appeared to contain phylogenetic information and produced more sensible estimates of primate divergence dates than did the amino acid sequences. Estimated dates varied considerably depending on the data type, the calibration point, and the substitution model but differed little among the four tree topologies used. We conclude that the calibration derived from the primate fossil record is too recent to be reliable; we also point out a number of problems in date estimation when the molecular clock does not hold. Despite these obstacles, we derived estimates of primate divergence dates that were well supported by the data and were generally consistent with the paleontological record. Estimation of the mouse-rat divergence date, however, was problematic.     

An insect molecular clock dates the origin of the insects and accords with palaeontological and biogeographic landmarks

A unified understanding of >390 Myr of insect evolution requires insight into their origin. Molecular clocks are widely applied for evolutionary dating, but clocks for the class Insecta have remained elusive. We now define a robust nucleotide and amino acid mitochondrial molecular clock encompassing five insect orders, including the Blattaria (cockroaches), Orthoptera (crickets and locusts), Hemiptera (true bugs), Diptera, and Lepidoptera (butterflies and moths). Calibration of the clock using one of the earliest, most extensive fossil records for insects (the early ancestors of extant Blattaria) was congruent with all available insect fossils, with biogeographic history, with the Cambrian explosion, and with independent dating estimates from Lepidopteran families. In addition, dates obtained from both nucleotide and amino acid clocks were congruent with each other. Of particular interest to vector biology is the early date of the emergence of triatomine bugs (99.8-93.5 MYA), coincident with the formation of the South American continent during the breakup of Gondwanaland. More generally, we reveal the insects arising from a common ancestor with the Anostraca (fairy shrimps) at around the Silurian-Ordovician boundary (434.2-421.1 MYA) coinciding with the earliest plant megafossil. We explore Tilyard's theory proposing that the terrestrial transition of the aquatic arthropod ancestor to the insects is associated with a particular plant group (early vascular plants). The major output of the study is a comprehensive series of dates for deep-branching points within insect evolution that can act as calibration points for further dating studies within insect families and genera.     

Comparative evolutionary analysis of chalcone synthase and alcohol dehydrogenase loci in Arabidopsis, Arabis, and related genera (Brassicaceae)

We analyzed sequence variation for chalcone synthase (CHS:) and alcohol dehydrogenase (ADH:) loci in 28 species in the genera ARABIDOPSIS: and ARABIS: and related taxa from tribe Arabideae. CHS: was single-copy in nearly all taxa examined, while ADH: duplications were found in several species. Phylogenies constructed from both loci confirmed that the closest relatives of Arabidopsis thaliana include Arabidopsis lyrata, Arabidopsis petraea, and Arabidopsis halleri (formerly in the genus CARDAMINOPSIS:). Slightly more distant are the North American n = 7 Arabis (BOECHERA:) species. The genus ARABIS: is polyphyletic-some unrelated species appear within this taxonomic classification, which has little phylogenetic meaning. Fossil pollen data were used to compute a synonymous substitution rate of 1.5 x 10 substitutions per site per year for both CHS: and Adh. Arabidopsis thaliana diverged from its nearest relatives about 5 MYA, and from Brassica roughly 24 MYA. Independent molecular and fossil data from several sources all provide similar estimates of evolutionary timescale in the Brassicaceae.     

Unraveling the evolutionary radiation of the thoracican barnacles using molecular and morphological evidence: a comparison of several divergence time estimation approaches

The Thoracica includes the ordinary barnacles found along the sea shore and is the most diverse and well-studied superorder of Cirripedia. However, although the literature abounds with scenarios explaining the evolution of these barnacles, very few studies have attempted to test these hypotheses in a phylogenetic context. The few attempts at phylogenetic analyses have suffered from a lack of phylogenetic signal and small numbers of taxa. We collected DNA sequences from the nuclear 18S, 28S, and histone H3 genes and the mitochondrial 12S and 16S genes (4,871 bp total) and data for 37 adult and 53 larval morphological characters from 43 taxa representing all the extant thoracican suborders (except the monospecific Brachylepadomorpha). Four Rhizocephala (highly modified parasitic barnacles) taxa and a Rhizocephala + Acrothoracica (burrowing barnacles) hypothetical ancestor were used as the outgroup for the molecular and morphological analyses, respectively. We analyzed these data separately and combined using maximum likelihood (ML) under "hill-climbing" and genetic algorithm heuristic searches, maximum parsimony procedures, and Bayesian inference coupled with Markov chain Monte Carlo techniques under mixed and homogeneous models of nucleotide substitution. The resulting phylogenetic trees answered key questions in barnacle evolution. The four-plated Iblomorpha were shown as the most primitive thoracican, and the plateless Heteralepadomorpha were placed as the sister group of the Lepadomorpha. These relationships suggest for the first time in an invertebrate that exoskeleton biomineralization may have evolved from phosphatic to calcitic. Sessilia (nonpedunculate) barnacles were depicted as monophyletic and appear to have evolved from a stalked (pedunculate) multiplated (5+) scalpelloidlike ancestor rather than a five-plated lepadomorphan ancestor. The Balanomorpha (symmetric sessile barnacles) appear to have the following relationship: (Chthamaloidea(Coronuloidea(Tetraclitoidea, Balanoidea))). Thoracican divergence times were estimated under ML-based local clock, Bayesian, and penalized likelihood approaches using an 18S data set and three calibration points: Heteralepadomorpha = 530 million years ago (MYA), Scalpellomorpha = 340 MYA, and Verrucomorpha = 120 MYA. Estimated dates varied considerably within and between approaches depending on the calibration point. Highly parameterized local clock models that assume independent rates (r > or = 15) for confamilial or congeneric species generated the most congruent estimates among calibrations and agreed more closely with the barnacle fossil record. Reasonable estimates were also obtained under the Bayesian procedure of Kishino et al. (2001, Mol. Biol. Evol. 18:352-361) but using multiple calibrations. Most of the dates estimated under the Bayesian procedure of Aris-Brosou and Yang (2002, Syst. Biol. 51:703-714) and the penalized likelihood method using single and/or multiple calibrations were inconsistent among calibrations and did not fit the fossil record.     

Bayesian models of episodic evolution support a late precambrian explosive diversification of the Metazoa

Multicellular animals, or Metazoa, appear in the fossil records between 575 and 509 million years ago (MYA). At odds with paleontological evidence, molecular estimates of basal metazoan divergences have been consistently older than 700 MYA. However, those date estimates were based on the molecular clock hypothesis, which is almost always violated. To relax this hypothesis, we have implemented a Bayesian approach to describe the change of evolutionary rate over time. Analysis of 22 genes from the nuclear and the mitochondrial genomes under the molecular clock assumption produced old date estimates, similar to those from previous studies. However, by allowing rates to vary in time and by taking small species-sampling fractions into account, we obtained much younger estimates, broadly consistent with the fossil records. In particular, the date of protostome-deuterostome divergence was on average 582 +/- 112 MYA. These results were found to be robust to specification of the model of rate change. The clock assumption thus had a dramatic effect on date estimation. However, our results appeared sensitive to the prior model of cladogenesis, although the oldest estimates (791 +/- 246 MYA) were obtained under a suboptimal model. Bayes posterior estimates of evolutionary rates indicated at least one major burst of molecular evolution at the end of the Precambrian when protostomes and deuterostomes diverged. We stress the importance of assumptions about rates on date estimation and suggest that the large discrepancies between the molecular and fossil dates of metazoan divergences might partly be due to biases in molecular date estimation.     

Primate molecular divergence dates

With genomic data, alignments can be assembled that greatly increase the number of informative sites for analysis of molecular divergence dates. Here, we present an estimate of the molecular divergence dates for all of the major primate groups. These date estimates are based on a Bayesian analysis of approximately 59.8 kbp of genomic data from 13 primates and 6 mammalian outgroups, using a range of paleontologically supported calibration estimates. Results support a Cretaceous last common ancestor of extant primates (approximately 77 mya), an Eocene divergence between platyrrhine and catarrhine primates (approximately 43 mya), an Oligocene origin of apes and Old World monkeys (approximately 31 mya), and an early Miocene (approximately 18 mya) divergence of Asian and African great apes. These dates are examined in the context of other molecular clock studies.     

Proceedings of the SMBE Tri-National Young Investigators' Workshop 2005. Southern hemisphere springtails: could any have survived glaciation of Antarctica?

Throughout the Southern Hemisphere many terrestrial taxa have circum-Antarctic distributions. This pattern is generally attributed to ongoing dispersal (by wind, water, or migrating birds) or relict Gondwanan distributions. Few of these terrestrial taxa have extant representatives in Antarctica, but such taxa would contribute to our understanding of the evolutionary origins of the continental Antarctic fauna. Either these taxa have survived the harsh climate cooling in Antarctica over the last 23 Myr (Gondwanan/vicariance origin) or they have dispersed there more recently (<2 MYA). In this context, we examined mtDNA (COI) sequence variation in Cryptopygus and related extant Antarctic and subantarctic terrestrial springtails (Collembola). Sequence divergence was estimated under a maximum likelihood model (general time reversible+I+Gamma) between individuals from subantarctic islands, Australia, New Zealand, Patagonia, Antarctic Peninsula, and continental Antarctica. Recent dispersal/colonization (<2 MYA) of Cryptopygus species was inferred between some subantarctic islands, and there was a close association between estimated times of divergences based on a molecular clock and proposed geological ages of islands. Most lineages generally grouped according to geographic proximity or by inferred dispersal/colonization pathways. In contrast, the deep divergences found for the four endemic Antarctic species indicate that they represent a continuous chain of descent dating from the break up of Gondwana to the present. We suggest that the diversification of these springtail species (21-11 MYA) in ice-free glacial refugia throughout the Trans-Antarctic Mountains was caused by the glaciation of the Antarctic continent during the middle to late Miocene.     

Additive Uncorrelated Relaxed Clock Models for the Dating of Genomic Epidemiology Phylogenies

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.     

Mitochondrial DNA evolution at a turtle's pace: evidence for low genetic variability and reduced microevolutionary rate in the Testudines.

Evidence is compiled suggesting a slowdown in mean microevolutionary rate for turtle mitochondrial DNA (mtDNA). Within each of six species or species complexes of Testudines, representing six genera and three taxonomic families, sequence divergence estimates derived from restriction assays are consistently lower than expectations based on either (a) the dates of particular geographic barriers with which significant mtDNA genetic clades appear associated or (b) the magnitudes of sequence divergence between mtDNA clades in nonturtle species that otherwise exhibit striking phylogeographic concordance with the genetic partitions in turtles. Magnitudes of the inferred rate slowdowns average eightfold relative to the "conventional" mtDNA clock calibration of 2%/Myr sequence divergence between higher animal lineages. Reasons for the postulated deceleration remain unknown, but two intriguing correlates are (a) the exceptionally long generation length most turtles and (b) turtles' low metabolic rate. Both factors have been suspected of influencing evolutionary rates in the DNA sequences of some other vertebrate groups. Uncertainities about the dates of cladogenetic events in these Testudines leave room for alternatives to the slowdown interpretation, but consistency in the direction of the inferred pattern, across several turtle species and evolutionary settings, suggests the need for caution in acceptance of a universal mtDNA-clock calibration for higher animals.     

Maximum-likelihood divergence date estimates based on rRNA gene sequences suggest two scenarios of Trypanosoma cruzi intraspecific evolution.

The phylogenetic relationships of Trypanosoma cruzi strains were inferred using maximum-likelihood from complete 18S rDNA sequences and D7-24Salpha rDNA regions from 20 representative strains of T. cruzi. For this we sequenced the 18S rDNA of 14 strains and the D7-24Salpha rDNA of four strains and aligned them to previously published sequences. Phylogenies inferred from these data sets identified four groups, named Riboclades 1, 2, 3, and 4, and a basal dichotomy that separated Riboclade 1 from Riboclades 2, 3, and 4. Substitution models and other parameters were optimized by hierarchical likelihood tests, and our analysis of the 18S rDNA molecular clock by the likelihood ratio test suggests that a taxa subset encompassing all 2,150 positions in the alignment supports rate constancy among lineages. The present analysis supports the notion that divergence dates of T. cruzi Riboclades can be estimated from 18S rDNA sequences and therefore, we present alternative evolutionary scenarios based on two different views of T. cruzi intraspecific divergence. The first assumes a faster evolutionary rate, which suggests that the divergence between T. cruzi I and II and the extant strains occurred in the Tertiary period (37-18 MYA). The other, which supports the hypothesis that the divergence between T. cruzi I and II occurred in the Cretaceous period (144-65 MYA) and the divergence of the extant strains occurred in the Tertiary period of the Cenozoic era (65-1.8 MYA), is consistent with our previously proposed hypothesis of divergence by geographical isolation and mammalian host coevolution.     

Integration of Bayesian molecular clock methods and fossil-based soft bounds reveals early Cenozoic origin of African lacertid lizards

Background  

Although current molecular clock methods offer greater flexibility in modelling evolutionary events, calibration of the clock with dates from the fossil record is still problematic for many groups. Here we implement several new approaches in molecular dating to estimate the evolutionary ages of Lacertidae, an Old World family of lizards with a poor fossil record and uncertain phylogeny. Four different models of rate variation are tested in a new program for Bayesian phylogenetic analysis called TreeTime, based on a combination of mitochondrial and nuclear gene sequences. We incorporate paleontological uncertainty into divergence estimates by expressing multiple calibration dates as a range of probabilistic distributions. We also test the reliability of our proposed calibrations by exploring effects of individual priors on posterior estimates.     

Local molecular clocks in three nuclear genes: divergence times for rodents and other mammals and incompatibility among fossil calibrations

Reconstructing the chronology of mammalian evolution is a debated issue between molecule- and fossil-based inferences. A methodological limitation of molecules is the evolutionary rate variation among lineages, precluding the application of the global molecular clock. We considered 2422 first and second codon positions of the combined ADRA2B, IRBP, and vWF nuclear genes for a well-documented set of placentals including an extensive sampling of rodents. Using seven independent calibration points and a maximum-likelihood framework, we evaluated whether molecular and paleontological estimates of mammalian divergence dates may be reconciled by the local molecular clocks approach, allowing local constancy of substitution rates with variations at larger phylogenetic scales. To handle the difficulty of choosing among all possible rate assignments for various lineages, local molecular clocks were based on the results of branch-length and two-cluster tests. Extensive lineage-specific variation of evolutionary rates was detected, even among rodents. Cross-calibrations indicated some incompatibilities between divergence dates based on different paleontological references. To decrease the impact of a single calibration point, estimates derived from independent calibrations displaying only slight reciprocal incompatibility were averaged. The divergence dates inferred for the split between mice and rats (approximately 13-19 Myr) was younger than previously published molecular estimates. The most recent common ancestors of rodents, primates and rodents, boreoeutherians, and placentals were estimated to be, respectively, approximately 60, 70, 75, and 78 Myr old. Global clocks, local clocks, and quartet dating analyses suggested a Late Cretaceous origin of the crown placental clades followed by a Tertiary radiation of some placental orders like rodents.     

Ancient dates or accelerated rates? Morphological clocks and the antiquity of placental mammals

Analyses of a comprehensive morphological character matrix of mammals using ‘relaxed’ clock models (which simultaneously estimate topology, divergence dates and evolutionary rates), either alone or in combination with an 8.5 kb nuclear sequence dataset, retrieve implausibly ancient, Late Jurassic–Early Cretaceous estimates for the initial diversification of Placentalia (crown-group Eutheria). These dates are much older than all recent molecular and palaeontological estimates. They are recovered using two very different clock models, and regardless of whether the tree topology is freely estimated or constrained using scaffolds to match the current consensus placental phylogeny. This raises the possibility that divergence dates have been overestimated in previous analyses that have applied such clock models to morphological and total evidence datasets. Enforcing additional age constraints on selected internal divergences results in only a slight reduction of the age of Placentalia. Constraining Placentalia to less than 93.8 Ma, congruent with recent molecular estimates, does not require major changes in morphological or molecular evolutionary rates. Even constraining Placentalia to less than 66 Ma to match the ‘explosive’ palaeontological model results in only a 10- to 20-fold increase in maximum evolutionary rate for morphology, and fivefold for molecules. The large discrepancies between clock- and fossil-based estimates for divergence dates might therefore be attributable to relatively small changes in evolutionary rates through time, although other explanations (such as overly simplistic models of morphological evolution) need to be investigated. Conversely, dates inferred using relaxed clock models (especially with discrete morphological data and MrBayes) should be treated cautiously, as relatively minor deviations in rate patterns can generate large effects on estimated divergence dates.     

Catarrhine primate divergence dates estimated from complete mitochondrial genomes: concordance with fossil and nuclear DNA evidence

Accurate divergence date estimates improve scenarios of primate evolutionary history and aid in interpretation of the natural history of disease-causing agents. While molecule-based estimates of divergence dates of taxa within the superfamily Hominoidea (apes and humans) are common in the literature, few such estimates are available for the Cercopithecoidea (Old World monkeys), the sister taxon of the hominoids in the primate infraorder Catarrhini. To help fill this gap, we have sequenced the entire mitochondrial DNA (mtDNA) genomes from a representative of three cercopithecoid tribes, Cercopithecini (Chlorocebus aethiops), Colobini (Colobus guereza), and Presbytini (Trachypithecus obscurus), and analyzed these new data together with other catarrhine mtDNA genomes available in public databases. Molecular divergence date estimates are dependent on calibration points gleaned from the paleontological record. We defined criteria for the selection of good calibration points and identified three points meeting these criteria: Homo-Pan, 6.0 Ma; Pongo-hominines, 14.0 Ma; hominoid/cercopithecoid, 23.0 Ma. Because a uniform molecular clock does not fit the catarrhine mtDNA data, we estimated divergence dates using a penalized likelihood and a Bayesian method, both of which take into account the effects of rate differences on lineages, phylogenetic tree structure, and multiple calibration points. The penalized likelihood method applied to the coding regions of the mtDNA genome yielded the following divergence date estimates, with approximate 95% confidence intervals: cercopithecine-colobine, 16.2 (14.4-17.9) Ma; colobin-presbytin, 10.9 (9.6-12.3) Ma; cercopithecin-papionin, 11.6 (10.3-12.9) Ma; and Macaca-Papio, 9.8 (8.6-10.9) Ma. Within the hominoids, the following dates were inferred: hylobatid-hominid, 16.8 (15.0-18.5) Ma; Gorilla-Homo+Pan, 8.1 (7.1-9.0) Ma; Pongo pygmaeus pygmaeus-P. p. abelii, 4.1 (3.5-4.7) Ma; and Pan troglodytes-P. paniscus, 2.4 (2.0-2.7) Ma. These dates were similar to those found using penalized likelihood on other subsets of the data, but slightly younger than several of the Bayesian estimates.     

Resolving the phylogenetic relationships and evolutionary history of the Drosophila virilis group using multilocus data

The Drosophila virilis group is one of the major lineages of Drosophila previously recognised and it has been used as a model for different types of studies. It comprises 13 species whose phylogenetic relationships are not well resolved. In the present study, six nuclear genes (Adh, fused, Gpdh, NonA, CG9631 and CG7219) and the mitochondrial ribosomal RNA genes (12S-16S) have been used to estimate the evolutionary tree of the group using different methods of phylogenetic reconstruction. Different competing evolutionary hypotheses have also been compared using the Approximately Unbiased test to further evaluate the robustness of the inferred trees. Results are, in general, consistent with previous studies in recovering the four major lineages of the group (D. virilis phylad, Drosophila montana subphylad, Drosophila kanekoi subphylad and Drosophila littoralis subphylad), although D. kanekoi, D. littoralis and Drosophila ezoana are here inferred to be more closely related to the D. virilis phylad than to the D. montana subphylad. The age of the crown group, estimated with a Bayesian method that assumes a relaxed molecular clock, is placed in the late Miocene (～ 10 Mya). The oldest lineages also appeared during this period (～ 7.5 to ～ 8.9 Mya), while the ages of the basal nodes of the montana subphylad and the virilis phylad are located in the early Pliocene (～ 4.9 and ～ 4.1 Mya). Major cladogenesis events correlate to geological and palaeoclimatic occurrences that most likely affected the freshwater and deciduous forests where these species are found. The inferred biogeographical history of the group, based on the statistical dispersal-vicariance analysis, indicates that the last common ancestor of the group had a Holarctic distribution from which the North American and the Eurasian lineages evolved as a result of a vicariant event.     

Molecular phylogeny and evolution of Sorex shrews (Soricidae: insectivora) inferred from mitochondrial DNA sequence data

Shrews of the genus Sorex are characterized by a Holarctic distribution, and relationships among extant taxa have never been fully resolved. Phylogenies have been proposed based on morphological, karyological, and biochemical comparisons, but these analyses often produced controversial and contradictory results. Phylogenetic analyses of partial mitochondrial cytochrome b gene sequences (1011 bp) were used to examine the relationships among 27 Sorex species. The molecular data suggest that Sorex comprises two major monophyletic lineages, one restricted mostly to the New World and one with a primarily Palearctic distribution. Furthermore, several sister-species relationships are revealed by the analysis. Based on the split between the Soricinae and Crocidurinae subfamilies, we used a 95% confidence interval for both the calibration of a molecular clock and the subsequent calculation of major diversification events within the genus Sorex. Our analysis does not support an unambiguous acceleration of the molecular clock in shrews, the estimated rate being similar to other estimates of mammalian mitochondrial clocks. In addition, the data presented here indicate that estimates from the fossil record greatly underestimate divergence dates among Sorex taxa.     

Playing chicken (Gallus gallus): methodological inconsistencies of molecular divergence date estimates due to secondary calibration points

For any given taxonomic divergence event, one may find in the literature a wide range of time estimates. Many factors contribute to the variation in molecular date estimates for the same evolutionary event. High on the list is the choice of calibration points for converting genetic distances into evolutionary rates and, subsequently, into dates of divergence. In this study, we investigate one critical source of error in estimating divergence times, i.e. the use of secondary calibration points, which are divergence time estimates that have been derived from one molecular dataset on the basis of a primary external calibration point, and which are used again independently of the original external calibration point on a second dataset. Unless particular care is exercised, this practice leads to internal inconsistencies, and the inferred dates of divergence are by necessity unreliable. We present a consistency test for assessing the reliability of divergence time estimates based on secondary calibration points. As a case study, we examine recent estimates of divergence times among phyla and kingdoms based on multiple nuclear protein-coding genes, and show that they fail the consistency test.     

A protocol for temporal calibration of general area cladograms

Aim To describe a protocol for incorporating a temporal dimension into historical biogeographical analysis, while maintaining the essential independence of all datasets, involving the generation of general area cladograms. Location Global. Methods General area cladograms (GACs) are a reconstruction of the evolutionary history of a set of areas and unrelated clades within those areas. Nodes on a GAC correspond to speciation events in a group of taxa; general nodes are those at which multiple unrelated clades speciate. We undertake temporal calibration of GACs using molecular clock estimates of splitting events between extant taxa as well as first appearance data from the fossil record. We present two examples based on re‐analysis of previously published data: first, a temporally calibrated GAC generated from secondary Brooks parsimony analysis (BPA) of six extant bird clades from the south‐west of North America using molecular clock estimates of divergence times; and second, an analysis of African Neogene mammals based on a phylogenetic analysis for comparing trees (PACT) analysis. Results A hypothetical example demonstrates how temporal calibration reveals potentially critical information about the timing of both unique and general events, while also illustrating instances of incongruence between dates generated from molecular clock estimates and fossils. For the African Neogene mammal dataset, our analysis reveals that most mammal clades underwent geodispersal associated with the Neogene climatic optimum (c. 16 Ma) and vicariant speciation in central Africa correlated with increased aridity and cooler temperatures around 2.5 Ma. Main conclusions Temporally calibrated GACs are valuable tools for assessing whether coordinated patterns of speciation are associated with large‐scale climatic or tectonic phenomena.     

Evolution of the Diatoms (Bacillariophyta): IV. A Reconstruction of Their Age from Small Subunit rRNA Coding Regions and the Fossil Record

Small subunit ribosomal RNA (ssu rRNA) coding regions from 30 diatoms, 3 oomycetes, and 6 pelagophytes were used to construct linearized trees, maximum-likelihood trees, and neighbor-joining trees inferred from both unweighted and weighted distances. Stochastic accumulation of sequence substitutions among the diatoms was assessed with relative rate tests. Pennate diatoms evolved relatively slowly but within the limits set by a stochastic model; centric diatoms exceeded those limits. A rate distribution test was devised to identify those taxa showing an aberrant distribution of base substitutions within the ssu rRNA coding region. First appearance dates of diatom taxa from the fossil record were regressed against their corresponding branch lengths to infer the average and earliest possible age for the origin of the diatoms, the pennate diatoms, and the centric diatom order Thalassiosirales. Our most lenient age estimate (based on the median-evolving diatom taxon in the maximum-likelihood tree or on the average branch length in a linearized tree) suggests that their average age is approximately 164–166 Ma, which is close to their earliest fossil record. Both calculations suggest that it is unlikely that diatoms existed prior to 238–266 Ma. Rate variation among the diatoms' ssu rRNA coding regions and uncertainties associated with the origin of extant taxa in the fossil record contribute significantly to the variation in age estimates obtained. Different evolutionary models and the exclusion of fast or slow evolving taxa did not significantly affect age estimates; however, the inclusion of aberrantly fast evolving taxa did. Our molecular clock calibrations indicate that the rRNA coding regions in the diatoms are evolving at approximately 1% per 18 to 26 Ma, which is the fastest substitution rate reported in any pro- or eukaryotic group of organisms to date.