Testing the molecular clock: molecular and paleontological estimates of divergence times in the Echinoidea (Echinodermata)

The phylogenetic relationships of 46 echinoids, with representatives from 13 of the 14 ordinal-level clades and about 70% of extant families commonly recognized, have been established from 3 genes (3,226 alignable bases) and 119 morphological characters. Morphological and molecular estimates are similar enough to be considered suboptimal estimates of one another, and the combined data provide a tree that, when calibrated against the fossil record, provides paleontological estimates of divergence times and completeness of their fossil record. The order of branching on the cladogram largely agrees with the stratigraphic order of first occurrences and implies that their fossil record is more than 85% complete at family level and at a resolution of 5-Myr time intervals. Molecular estimates of divergence times derived from applying both molecular clock and relaxed molecular clock models are concordant with estimates based on the fossil record in up to 70% of cases, with most concordant results obtained using Sanderson's semiparametric penalized likelihood method and a logarithmic-penalty function. There are 3 regions of the tree where molecular and fossil estimates of divergence time consistently disagree. Comparison with results obtained when molecular divergence dates are estimated from the combined (morphology + gene) tree suggests that errors in phylogenetic reconstruction explain only one of these. In another region the error most likely lies with the paleontological estimates because taxa in this region are demonstrated to have a very poor fossil record. In the third case, morphological and paleontological evidence is much stronger, and the topology for this part of the molecular tree differs from that derived from the combined data. Here the cause of the mismatch is unclear but could be methodological, arising from marked inequality of molecular rates. Overall, the level of agreement reached between these different data and methodological approaches leads us to believe that careful application of likelihood and Bayesian methods to molecular data provides realistic divergence time estimates in the majority of cases (almost 80% in this specific example), thus providing a remarkably well-calibrated phylogeny of a character-rich clade of ubiquitous marine benthic invertebrates.     

Timing and tempo of primate speciation

Published molecular clocks for primates are used to estimate typical divergence times for phylogroups (1.6 Ma), species (3.3 Ma), sister species (2.7 Ma), genera (8.9 Ma) and sister genera (8.6 Ma). Significant median differences exist between major groups (infraorders and superfamilies) for various divergence times. These data are employed to estimate typical maximum duration of speciation. Typical primate values (1.1 Ma) suggest this process to be faster than is characteristic of many vertebrates. However, after considering divergence times for hybridizing congeneric and confamilial primates, this value is likely only to estimate the commencement of prezygotic isolating mechanisms, rather than the completion of reproductive isolation. Thus, speciation typically takes around 1.0 Ma to more than 4.0 Ma to occur, depending on whether prezygotic or post-zygotic isolating mechanisms are emphasized. Typical primate genus age is around 5.3 Ma, but we note differences among major groups. In light of these estimates, the classification of humans and chimpanzees is reconsidered using a molecular yardstick approach. Three taxonomic frameworks may flow from molecular analyses, all of them having major implications for understanding the evolution of humans and chimpanzees.     

The endemic gastropod fauna of Lake Titicaca: correlation between molecular evolution and hydrographic history

Lake Titicaca, situated in the Altiplano high plateau, is the only ancient lake in South America. This 2- to 3-My-old (where My is million years) water body has had a complex history that included at least five major hydrological phases during the Pleistocene. It is generally assumed that these physical events helped shape the evolutionary history of the lake's biota. Herein, we study an endemic species assemblage in Lake Titicaca, composed of members of the microgastropod genus Heleobia, to determine whether the lake has functioned as a reservoir of relic species or the site of local diversification, to evaluate congruence of the regional paleohydrology and the evolutionary history of this assemblage, and to assess whether the geographic distributions of endemic lineages are hierarchical. Our phylogenetic analyses indicate that the Titicaca/Altiplano Heleobia fauna (together with few extralimital taxa) forms a species flock. A molecular clock analysis suggests that the most recent common ancestor (MRCAs) of the Altiplano taxa evolved 0.53 (0.28-0.80) My ago and the MRCAs of the Altiplano taxa and their extralimital sister group 0.92 (0.46-1.52) My ago. The endemic species of Lake Titicaca are younger than the lake itself, implying primarily intralacustrine speciation. Moreover, the timing of evolutionary branching events and the ages of two precursors of Lake Titicaca, lakes Cabana and Ballivián, is congruent. Although Lake Titicaca appears to have been the principal site of speciation for the regional Heleobia fauna, the contemporary spatial patterns of endemism have been masked by immigration and/or emigration events of local riverine taxa, which we attribute to the unstable hydrographic history of the Altiplano. Thus, a hierarchical distribution of endemism is not evident, but instead there is a single genetic break between two regional clades. We also discuss our findings in relation to studies of other regional biota and suggest that salinity tolerance was the most likely limiting factor in the evolution of Altiplano species flocks.     

Timing the eastern Asian-eastern North American floristic disjunction: molecular clock corroborates paleontological estimates

Sequence data of the chloroplast gene rbcL were used to estimate the time of the well-known eastern Asian-eastern North American floristic disjunction. Sequence divergence of rbcL was examined for 22 species of 11 genera (Campsis, Caulophyllum, Cornus, Decumaria, Liriodendron, Menispermum, Mitchella, Pachysandra, Penthorum, Podophyllum, and Phryma) representing a diverse array of flowering plants occurring disjunctly in eastern Asia and eastern North America. Divergence times of putative disjunct species pairs were estimated from synonymous substitutions, using rbcL molecular clocks calibrated for Cornus. Relative rate tests were performed to assess rate constancy of rbcL evolution among lineages. Corrections of estimates of divergence times for each species pair were made based on rate differences of rbcL between Cornus and other species pairs. Results of these analyses indicate that the time of divergence of species pairs examined ranges from 12.56 +/- 4.30 million years to recent (<0.31 million years), with most within the last 10 million years (in the late Miocene and Pliocene). These results suggest that the isolation of most morphologically similar disjunct species in eastern Asia and eastern North America occurred during the global climatic cooling period that took place throughout the late Tertiary and Quaternary. This estimate is closely correlated with paleontological evidence and in agreement with the hypothesis that considers the eastern Asian-eastern North American floristic disjunction to be the result of the range restriction of a once more or less continuously distributed mixed mesophytic forest of the Northern Hemisphere that occurred during the late Tertiary and Quaternary. This implies that in most taxa the disjunction may have resulted from vicariance events. However, long-distance dispersal may explain the disjunct distribution of taxa with low divergence, such as Menispermum.     

Why do species vary in their rate of molecular evolution?

Despite hopes that the processes of molecular evolution would be simple, clock-like and essentially universal, variation in the rate of molecular evolution is manifest at all levels of biological organization. Furthermore, it has become clear that rate variation has a systematic component: rate of molecular evolution can vary consistently with species body size, population dynamics, lifestyle and location. This suggests that the rate of molecular evolution should be considered part of life-history variation between species, which must be taken into account when interpreting DNA sequence differences between lineages. Uncovering the causes and correlates of rate variation may allow the development of new biologically motivated models of molecular evolution that may improve bioinformatic and phylogenetic analyses.     

Time-dependent rates of molecular evolution

For over half a century, it has been known that the rate of morphological evolution appears to vary with the time frame of measurement. Rates of microevolutionary change, measured between successive generations, were found to be far higher than rates of macroevolutionary change inferred from the fossil record. More recently, it has been suggested that rates of molecular evolution are also time dependent, with the estimated rate depending on the timescale of measurement. This followed surprising observations that estimates of mutation rates, obtained in studies of pedigrees and laboratory mutation-accumulation lines, exceeded long-term substitution rates by an order of magnitude or more. Although a range of studies have provided evidence for such a pattern, the hypothesis remains relatively contentious. Furthermore, there is ongoing discussion about the factors that can cause molecular rate estimates to be dependent on time. Here we present an overview of our current understanding of time-dependent rates. We provide a summary of the evidence for time-dependent rates in animals, bacteria and viruses. We review the various biological and methodological factors that can cause rates to be time dependent, including the effects of natural selection, calibration errors, model misspecification and other artefacts. We also describe the challenges in calibrating estimates of molecular rates, particularly on the intermediate timescales that are critical for an accurate characterization of time-dependent rates. This has important consequences for the use of molecular-clock methods to estimate timescales of recent evolutionary events.     

The genome as a life-history character: why rate of molecular evolution varies between mammal species

DNA sequences evolve at different rates in different species. This rate variation has been most closely examined in mammals, revealing a large number of characteristics that can shape the rate of molecular evolution. Many of these traits are part of the mammalian life-history continuum: species with small body size, rapid generation turnover, high fecundity and short lifespans tend to have faster rates of molecular evolution. In addition, rate of molecular evolution in mammals might be influenced by behaviour (such as mating system), ecological factors (such as range restriction) and evolutionary history (such as diversification rate). I discuss the evidence for these patterns of rate variation, and the possible explanations of these correlations. I also consider the impact of these systematic patterns of rate variation on the reliability of the molecular date estimates that have been used to suggest a Cretaceous radiation of modern mammals, before the final extinction of the dinosaurs.     

The local-clock permutation test: a simple test to compare rates of molecular evolution on phylogenetic trees

Rates of molecular evolution vary substantially between lineages, and a growing effort is directed at uncovering the causes and consequences of this variation. Comparing local-clocks (rates of molecular evolution estimated from different sets of branches of a phylogenetic tree) is a common tool in this research effort. Here, I show that a commonly used test (the Likelihood Ratio Test, LRT) will not be statistically valid for comparing local-clocks in most cases. Instead, I propose the local-clock permutation test (LCPT), a simple test that can be used to test the significance of differences between local-clocks. The LCPT could also be used to test for differences between any parameter that can be assigned to individual branches on a phylogenetic tree. Using simulated data, I show that the LCPT has good power to detect differences between local-clocks.     

Calibrating the avian molecular clock

Molecular clocks are widely used to date phylogenetic events, yet evidence supporting the rate constancy of molecular clocks through time and across taxonomic lineages is weak. Here, we present 90 candidate avian clock calibrations obtained from fossils and biogeographical events. Cross-validation techniques were used to identify and discard 16 inconsistent calibration points. Molecular evolution occurred in an approximately clock-like manner through time for the remaining 74 calibrations of the mitochondrial gene, cytochrome b . A molecular rate of approximately 2.1% (± 0.1%, 95% confidence interval) was maintained over a 12-million-year interval and across most of 12 taxonomic orders. Minor but significant variance in rates occurred across lineages but was not explained by differences in generation time, body size or latitudinal distribution as previously suggested.     

Sociality and the rate of molecular evolution

The molecular clock does not tick at a uniform rate in all taxa but may be influenced by species characteristics. Eusocial species (those with reproductive division of labor) have been predicted to have faster rates of molecular evolution than their nonsocial relatives because of greatly reduced effective population size; if most individuals in a population are nonreproductive and only one or few queens produce all the offspring, then eusocial animals could have much lower effective population sizes than their solitary relatives, which should increase the rate of substitution of "nearly neutral" mutations. An earlier study reported faster rates in eusocial honeybees and vespid wasps but failed to correct for phylogenetic nonindependence or to distinguish between potential causes of rate variation. Because sociality has evolved independently in many different lineages, it is possible to conduct a more wide-ranging study to test the generality of the relationship. We have conducted a comparative analysis of 25 phylogenetically independent pairs of social lineages and their nonsocial relatives, including bees, wasps, ants, termites, shrimps, and mole rats, using a range of available DNA sequences (mitochondrial and nuclear DNA coding for proteins and RNAs, and nontranslated sequences). By including a wide range of social taxa, we were able to test whether there is a general influence of sociality on rates of molecular evolution and to test specific predictions of the hypothesis: (1) that social species have faster rates because they have reduced effective population sizes; (2) that mitochondrial genes would show a greater effect of sociality than nuclear genes; and (3) that rates of molecular evolution should be correlated with the degree of sociality. We find no consistent pattern in rates of molecular evolution between social and nonsocial lineages and no evidence that mitochondrial genes show faster rates in social taxa. However, we show that the most highly eusocial Hymenoptera do have faster rates than their nonsocial relatives. We also find that social parasites (that utilize the workers from related species to produce their own offspring) have faster rates than their social relatives, which is consistent with an effect of lower effective population size on rate of molecular evolution. Our results illustrate the importance of allowing for phylogenetic nonindependence when conducting investigations of determinants of variation in rate of molecular evolution.     

Biogeographic calibrations for the molecular clock

Molecular estimates of evolutionary timescales have an important role in a range of biological studies. Such estimates can be made using methods based on molecular clocks, including models that are able to account for rate variation across lineages. All clock models share a dependence on calibrations, which enable estimates to be given in absolute time units. There are many available methods for incorporating fossil calibrations, but geological and climatic data can also provide useful calibrations for molecular clocks. However, a number of strong assumptions need to be made when using these biogeographic calibrations, leading to wide variation in their reliability and precision. In this review, we describe the nature of biogeographic calibrations and the assumptions that they involve. We present an overview of the different geological and climatic events that can provide informative calibrations, and explain how such temporal information can be incorporated into dating analyses.     

Studies of primate protein variation and evolution: Microelectrophoretic detection

Genetic variation at 16 protein and enzyme loci in Cercopithecus aethiops and several other primate species has been surveyed, using cellulose acetate microelectrophoresis. Resolution of several standard variant proteins is comparable to that achieved on starch gel or polyacrylamide gel. Although both intraspecific and interspecific variation was observed for some loci, the data generally support the concept that extracellular proteins are more likely to be polymorphic within a species, while intracellular proteins generally vary between species, if at all. These methodologies are particularly appropriate for screening multiple-locus variation in large numbers of samples; their relevance to studies of molecular evolution and evaluation of theories of kin selection is discussed.This research was supported in part by California State Agricultural Experiment Station Funds to the University of California, Berkeley, and the Wenner Gren Foundation.     

Time‐dependent estimates of molecular evolutionary rates: evidence and causes

We are writing in response to a recent critique by Emerson & Hickerson ( 2015 ), who challenge the evidence of a time‐dependent bias in molecular rate estimates. This bias takes the form of a negative relationship between inferred evolutionary rates and the ages of the calibrations on which these estimates are based. Here, we present a summary of the evidence obtained from a broad range of taxa that supports a time‐dependent bias in rate estimates, with a consideration of the potential causes of these observed trends. We also describe recent progress in improving the reliability of evolutionary rate estimation and respond to the concerns raised by Emerson & Hickerson ( 2015 ) about the validity of rates estimated from time‐structured sequence data. In doing so, we hope to dispel some misconceptions and to highlight several research directions that will improve our understanding of time‐dependent biases in rate estimates.     

A mitogenomic timescale for birds detects variable phylogenetic rates of molecular evolution and refutes the standard molecular clock

Current understanding of the diversification of birds is hindered by their incomplete fossil record and uncertainty in phylogenetic relationships and phylogenetic rates of molecular evolution. Here we performed the first comprehensive analysis of mitogenomic data of 48 vertebrates, including 35 birds, to derive a Bayesian timescale for avian evolution and to estimate rates of DNA evolution. Our approach used multiple fossil time constraints scattered throughout the phylogenetic tree and accounts for uncertainties in time constraints, branch lengths, and heterogeneity of rates of DNA evolution. We estimated that the major vertebrate lineages originated in the Permian; the 95% credible intervals of our estimated ages of the origin of archosaurs (258 MYA), the amniote-amphibian split (356 MYA), and the archosaur-lizard divergence (278 MYA) bracket estimates from the fossil record. The origin of modern orders of birds was estimated to have occurred throughout the Cretaceous beginning about 139 MYA, arguing against a cataclysmic extinction of lineages at the Cretaceous/Tertiary boundary. We identified fossils that are useful as time constraints within vertebrates. Our timescale reveals that rates of molecular evolution vary across genes and among taxa through time, thereby refuting the widely used mitogenomic or cytochrome b molecular clock in birds. Moreover, the 5-Myr divergence time assumed between 2 genera of geese (Branta and Anser) to originally calibrate the standard mitochondrial clock rate of 0.01 substitutions per site per lineage per Myr (s/s/l/Myr) in birds was shown to be underestimated by about 9.5 Myr. Phylogenetic rates in birds vary between 0.0009 and 0.012 s/s/l/Myr, indicating that many phylogenetic splits among avian taxa also have been underestimated and need to be revised. We found no support for the hypothesis that the molecular clock in birds "ticks" according to a constant rate of substitution per unit of mass-specific metabolic energy rather than per unit of time, as recently suggested. Our analysis advances knowledge of rates of DNA evolution across birds and other vertebrates and will, therefore, aid comparative biology studies that seek to infer the origin and timing of major adaptive shifts in vertebrates.     

Choosing the right molecular genetic markers for studying biodiversity: from molecular evolution to practical aspects

The use of molecular genetic markers (MGMs) has become widespread among evolutionary biologists, and the methods of analysis of genetic data improve rapidly, yet an organized framework in which scientists can work is lacking. Elements of molecular evolution are summarized to explain the origin of variation at the DNA level, its measures, and the relationships linking genetic variability to the biological parameters of the studied organisms. MGM are defined by two components: the DNA region(s) screened, and the technique used to reveal its variation. Criteria of choice belong to three categories: (1) the level of variability, (2) the nature of the information (e.g. dominance vs. codominance, ploidy, ... ) which must be determined according to the biological question and (3) some practical criteria which mainly depend on the equipment of the laboratory and experience of the scientist. A three-step procedure is proposed for drawing up MGMs suitable to answer given biological questions, and compiled data are organized to guide the choice at each step: (1) choice, determined by the biological question, of the level of variability and of the criteria of the nature of information, (2) choice of the DNA region and (3) choice of the technique.     

Molecular‐clock methods for estimating evolutionary rates and timescales

The molecular clock presents a means of estimating evolutionary rates and timescales using genetic data. These estimates can lead to important insights into evolutionary processes and mechanisms, as well as providing a framework for further biological analyses. To deal with rate variation among genes and among lineages, a diverse range of molecular‐clock methods have been developed. These methods have been implemented in various software packages and differ in their statistical properties, ability to handle different models of rate variation, capacity to incorporate various forms of calibrating information and tractability for analysing large data sets. Choosing a suitable molecular‐clock model can be a challenging exercise, but a number of model‐selection techniques are available. In this review, we describe the different forms of evolutionary rate heterogeneity and explain how they can be accommodated in molecular‐clock analyses. We provide an outline of the various clock methods and models that are available, including the strict clock, local clocks, discrete clocks and relaxed clocks. Techniques for calibration and clock‐model selection are also described, along with methods for handling multilocus data sets. We conclude our review with some comments about the future of molecular clocks.     

Patterns of molecular evolution and diversification in a biodiversity hotspot: the California Floristic Province

The California Floristic Province harbours more endemic plant and animal taxa and more identifiable subspecies than any other area of comparable size in North America. We present evidence that physical historical processes have resulted in congruent patterns of genetic diversity over the past 2-10 million years. Using a molecular clock approach we show that diversification and establishment of spatial genetic structure across six taxonomic groups coincide with the putative age of California's mountain ranges and aridification in the region. Our results demonstrate the importance of geographical barriers and climatological events to species diversification and the overall geographical structure of biodiversity. These results should facilitate conservation efforts in this biodiversity hotspot for taxa whose population genetic structure is still unknown and may suggest the potential utility of this approach in regional conservation planning efforts.     

Evolutionary processes in the genus Coreocarpus: insights from molecular phylogenetics

A molecular phylogenetic study of the plant genus Coreocarpus was conducted using nuclear (ITS) and plastid (rpl16 intron) DNA sequences, with phylogenies of the nuclear and plastid sequences highly congruent in defining a monophyletic group of six species (core Coreocarpus), although three other species often placed within the genus were excluded. Relationships within the genus are largely but not totally concordant with prior biosystematic studies. Despite strong molecular support, no morphological characters uniting the six species of core Coreocarpus have been identified; retention of plesiomorphic characters and the genetic lability of characters are two probable factors contributing to lack of consistent defining characters. The age of the core Coreocarpus is estimated at 1 million years because the basal species is endemic to a volcanic island that emerged in the past million years. Mapping the results of earlier breeding studies on the molecular phylogeny showed that use of cross-compatibility as a criterion for species delimitation would result in the recognition of paraphyletic species. Prior field, morphological, and biosystematic studies provided no indication of past hybridization in the evolution of Coreocarpus, and species in the genus appeared to be well defined morphologically. However, three instances of incongruence were observed. Two of these were between the nuclear and plastid partitions, and the third was between the morphological species assignment of one accession and the molecular data. If hybridization accounts for incongruence between the nuclear and plastid data, it occurred between species that now appear to be cross-incompatible and allopatric. The incongruence between morphological species assignment and the molecular data could be the result of parallel fixation of characters that have a simple genetic basis. This study suggests that the evolutionary history of Coreocarpus is much more complex than indicated from prior biosystematic investigations and that biosystematic and molecular phylogenetic studies may complement each other for elucidating the evolution and phylogeny of a group.