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.     

Calibration of galliform molecular clocks using multiple fossils and genetic partitions

For more than a century, members of the traditional avian order Galliformes (i.e., pheasants, partridges, junglefowl, and relatives) have been among the most intensively studied birds, but still a comprehensive timeframe for their evolutionary history is lacking. Thanks to a number of recent cladistic interpretations for several galliform fossils, candidates now exist that can potentially be used as accurate internal calibrations for molecular clocks. Here, we describe a molecular timescale for Galliformes based on cytochrome b and ND2 using nine mostly internal fossil-based anchorpoints. Beyond application of calibrations spanning the entire evolutionary history of Galliformes, care was taken to investigate the effects of calibration choice, substitution saturation, and rate heterogeneity among lineages on divergence time estimation. Results show broad consistency in time estimation with five out of the nine total calibrations. Our divergence time estimates, based on these anchorpoints, indicate that the early history of Galliformes took place in the Cretaceous, including the origin of the basal-most megapode and perhaps cracid lineages, but that the remaining morphological diversification likely started in the earliest Tertiary. The multi-calibration/multi-genetic partition approach used here highlights the importance of understanding the genetic saturation, variation, and rate constancy spectra for the accurate calculation of divergence times by use of molecular clocks.     

Origin and rapid diversification of a tropical moss

Molecular sequences rarely evolve at a constant rate. Yet, even in instances where a clock can be assumed or approximated for a particular set of sequences, fossils or clear patterns of vicariance are rarely available to calibrate the clock. Thus, obtaining absolute timing for diversification of natural lineages can prove difficult. Unfortunately, without absolute time we cannot develop a complete understanding of important evolutionary processes, including adaptive radiations and key innovations. In the present study, the coding sequence of the nuclear gene, glyceraldehyde 3-phosphate dehydrogenase (gpd), extracted from the paleotropical moss, Mitthyridium, was found to exhibit clocklike behavior and used to reconstruct the history of 80 distinct molecular lineages that cover the full geographic range of Mitthyridium. Two separate clades endemic to two geographically distinct oceanic archipelagos were revealed by this phylogenetic analysis. This allowed the use of island age (as derived from potassium-argon dating) as a maximum age of origin of each monophyletic group, providing two independent time anchors for the clock found in gpd, the final piece needed to study absolute time. Based on results from both maximum age calibrations, which separately yielded highly consistent estimates, the ancestor of this moss group arose approximately 8 million years ago, and then diversified at the rapid rate of 0.56 +/- 0.004 new lineages per million years. Such a rate is on par with the highest diversification rates reported in the literature including rapidly radiating insular groups like the Hawaiian silversword alliance, a classic example of an adaptive radiation. Using independent sources of data, it was found that neither the age nor diversification estimates were affected by the use of molecular lineages rather than species as the operational taxonomic units. Identifying the cause for this rapid diversification requires further testing, but it appears to be related to a general shift in reproductive strategy from sexual to asexual, which may be a key innovation for this young group.     

Application of molecular clocks in ornithology revisited

Molecular clocks have seen many applications in ornithology, but many applications are uncritical. In this commentary, I point out logical inconsistencies in many uses of clocks in avian molecular systematics. I call for greater rigor in application of molecular clocks – clocks should only be used when clocklike behavior has been tested and confirmed, and when appropriate calibrations are available. Authors and reviewers should insist on such rigor to assure that systematics is indeed scientific, and not just storytelling.     

Molecular phylogeny of the Eremias persica complex of the Iranian plateau (Reptilia: Lacertidae), based on mtDNA sequences

The Persian racerunner Eremias persica Blanford, 1875 is confined to the Iranian plateau, and forms one of the most widespread but rarely studied species of the family Lacertidae. With many local populations inhabiting a variety of habitats, and exhibiting considerable morphological, genetic, and ecological variations, it represents a species complex. We analysed sequences of mitochondrial cytochrome b and 12S ribosomal RNA (rRNA) genes derived from 13 geographically distant populations belonging to the E. persica complex. Using our knowledge of palaeogeographical events, a molecular clock was calibrated to assess the major events in fragmentation, radiation, and intraspecific variation. The sequence data strongly support a basal separation of the highland populations of western Iran from those of the open steppes and deserts, occurring in the east. The subsequent radiation, fragmentation, and evolution of these major assemblages have led to several discernable geographical lineages across the wide area of the Iranian plateau. The results indicate a middle‐Miocene origin for the clade as a whole. The first split, isolating the western and eastern clades, appears to have occurred 11–10 Mya. Further fragmentations and divergence within the major clades began about 8 Mya, with an evolutionary rate of 1.6% sequence divergence per million years among the lineages in the genes studied (combined data set). Molecular and morphological data strongly support a taxonomic revision of this species complex. At least four of the discovered clades should be raised to species, and two to subspecies, rank. © 2009 The Linnean Society of London, Zoological Journal of the Linnean Society, 2010, 158 , 641–660.     

A SPECIATIONAL HISTORY OF “LIVING FOSSILS”: MOLECULAR EVOLUTIONARY PATTERNS IN HORSESHOE CRABS

Horseshoe crabs' exceptional morphological conservatism over the past 150 My has led to their reputation as “living fossils,” but also has served to obscure phylogenetic relationships within the complex. Here we employ nucleotide sequences from two mitochondrial genes to assess molecular evolutionary rates and patterns among all extant horseshoe crab species. The American species Limulus polyphemus proved to be the sister taxon to a clade composed of the Asiatic species Tachypleus gigas, T. tridentatus, and Carcinoscorpius rotundicauda, whose relationships inter se were not resolved definitively. Both absolute and relative rate tests suggest a moderate slowdown in sequence evolution in horseshoe crabs. Nonetheless, dates of the lineage separations remain uncertain primarily because of reservations about molecular-clock calibrations resulting from large rate variances at examined loci across Arthropods and other animal lineages, as inferred in this and prior studies. Thus, ironically, separation dates as estimated by molecular evidence in general may remain most insecure in taxonomic groups for which such information is needed most—those lacking strong biogeographic or fossil benchmarks for internal-clock calibrations. In any event, the current results show that large numbers of molecular characters distinguish even these most morphologically conservative of organisms. Furthermore, comparisons against previously published mitochondrial sequence data in the morphologically dynamic hermit crab–king crab complex demonstrates that striking heterogeneity in levels of morphotypic differentiation can characterize Arthropod lineages at similar magnitudes of molecular divergence.     

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.     

Multiple Miocene Melastomataceae dispersal between Madagascar, Africa and India

Melastomataceae sensu stricto (excluding Memecylaceae) comprise some 3000 species in the neotropics, 1000 in Asia, 240 in Africa, and 230 in Madagascar. Previous family-wide morphological and DNA analyses have shown that the Madagascan species belong to at least three unrelated lineages, which were hypothesized to have arrived by trans-oceanic dispersal. An alternative hypothesis posits that the ancestors of Madagascan, as well as Indian, Melastomataceae arrived from Africa in the Late Cretaceous. This study tests these hypotheses in a Bayesian framework, using three combined sequence datasets analysed under a relaxed clock and simultaneously calibrated with fossils, some not previously used. The new fossil calibration comes from a re-dated possibly Middle or Upper Eocene Brazilian fossil of Melastomeae. Tectonic events were also tentatively used as constraints because of concerns that some of the family's fossils are difficult to assign to nodes in the phylogeny. Regardless of how the data were calibrated, the estimated divergence times of Madagascan and Indian lineages were too young for Cretaceous explanations to hold. This was true even of the oldest ages within the 95% credibility interval around each estimate. Madagascar's Melastomeae appear to have arrived from Africa during the Miocene. Medinilla, with some 70 species in Madagascar and two in Africa, too, arrived during the Miocene, but from Asia. Gravesia, with 100 species in Madagascar and four in east and west Africa, also appears to date to the Miocene, but its monophyly has not been tested. The study afforded an opportunity to compare divergence time estimates obtained earlier with strict clocks and single calibrations, with estimates based on relaxed clocks and different multiple calibrations and taxon sampling.     

Multiple fossil calibrations,nuclear loci and mitochondrial genomes provide new insight into biogeography and divergence timing for true seals (Phocidae,Pinnipedia)

Aim To better understand the historical biogeography of the true seals, Phocidae, by combining nuclear DNA (nDNA) and mitochondrial DNA (mtDNA) in a divergence time analysis using multiple fossil calibrations. Location Arctic, Antarctic, Pacific and Atlantic Oceans, Lake Baikal, Caspian Sea. Methods Fifteen nuclear genes totalling 8935 bp plus near‐complete mitochondrial genome sequences were used in a Bayesian divergence time analysis, incorporating eight soft‐bound fossil calibrations across the phylogeny. All species of true seals were included, plus the walrus, three otariids and seven carnivore outgroups. The majority of the nuclear sequences and four phocid mitochondrial genomes (plus three non‐phocid mitochondrial genomes) were newly generated for this study using DNA extracted from tissue samples; other sequences were obtained from GenBank. Results Using multiple nuclear genes and multiple fossil calibrations resulted in most divergence time estimations within Phocidae being much more recent than predicted by other molecular studies incorporating only mtDNA and using a single calibration point. A new phylogenetic hypothesis was recovered for the Antarctic seals. Main conclusions Incorporating multiple nuclear genes and fossil calibrations had a profound effect on the estimated divergence times. Most estimated divergences within Phocinae (Arctic seals) correspond to Arctic oceanic events and all occur within the last 12 Myr, a time when the Arctic and Atlantic oceans were freely exchanging and perennial Arctic sea ice existed, indicating that the Arctic seals may have had a longer association with ice than previously thought. The Monachinae (‘southern’ seals) split from the Phocinae c. 15 Ma on the eastern US coast. Several early trans‐Atlantic dispersals possibly occurred, leaving no living descendants, as divergence estimates suggest that the Monachus (monk seal) species divergences occurred in the western Atlantic c. 6 Ma, with the Mediterranean monk seal ancestor dispersing afterwards. The tribes Lobodontini (Antarctic seals) and Miroungini (elephant seals) are also estimated to have diverged in the eastern Atlantic c. 7 Ma and a single Lobodontini dispersal to Antarctica occurred shortly afterwards. Many of the newly estimated dates are used to infer how extinct lineages/taxa are allied with their living relatives.     

Tectonic blocks and molecular clocks

Evolutionary timescales have mainly used fossils for calibrating molecular clocks, though fossils only really provide minimum clade age constraints. In their place, phylogenetic trees can be calibrated by precisely dated geological events that have shaped biogeography. However, tectonic episodes are protracted, their role in vicariance is rarely justified, the biogeography of living clades and their antecedents may differ, and the impact of such events is contingent on ecology. Biogeographic calibrations are no panacea for the shortcomings of fossil calibrations, but their associated uncertainties can be accommodated. We provide examples of how biogeographic calibrations based on geological data can be established for the fragmentation of the Pangaean supercontinent: (i) for the uplift of the Isthmus of Panama, (ii) the separation of New Zealand from Gondwana, and (iii) for the opening of the Atlantic Ocean. Biogeographic and fossil calibrations are complementary, not competing, approaches to constraining molecular clock analyses, providing alternative constraints on the age of clades that are vital to avoiding circularity in investigating the role of biogeographic mechanisms in shaping modern biodiversity.This article is part of the themed issue ‘Dating species divergences using rocks and clocks’.     

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.     

GENOMIC INSIGHTS INTO EVOLUTIONARY RELATIONSHIPS AMONG HETEROKONT LINEAGES EMPHASIZING THE PHAEOPHYCEAE1

Heterokonts comprise a large and diverse group of organisms unified by the heterokont biflagellate condition. Monophyly of many of these lineages is well established, but evolutionary relationships among the various lineages remain elusive. Among these lineages, the brown algae (Phaeophyceae) are a monophyletic, taxonomically diverse, and ecologically critical group common to marine environments. Despite their biological and scientific importance, consensus regarding brown algal phylogeny and taxonomic relationships is missing. Our long‐term research goal is to produce a well‐resolved taxon‐rich phylogeny of the class to assess evolutionary patterns and taxonomic relationships among brown algal lineages and their relationship to other closely related heterokont groups. To accomplish this goal and augment existing loci for phaeophycean‐wide systematic studies, we generated expressed sequence tags (ESTs) from several major brown algal lineages and from the heterokont lineage representing the closest sister group to brown algae. To date, we have successfully constructed cDNA libraries for two lineages (Choristocarpus tenellus Zanardini and Schizocladia ischiensis E. C. Henry, Okuda et H. Kawai) and in the library test phase obtained up to 1,600 ESTs per organism. Annotation results showed a gene discovery rate of 45%–50% for each library revealing 500–700 unique genes from each organism. We have identified several potential genes for phylogenetic inference and used these loci for preliminary molecular clock analyses. Our molecular clock analysis suggests that the basal divergence in brown algae occurred around the time of the pennate‐centric diatom divergence. Here we report this analysis and other uses of ESTs in brown algal phylogenomics and the utility of these data for resolving the phylogeny of this group.     

Biogeographic history of geminate cirrhitoids (Perciformes: Cirrhitoidea) with east–west allopatric distributions across southern Australia, based on molecular data

The biogeographic history of three cirrhitoid species pairs with east–west allopatric distributions across southern Australia was examined by determining levels of mitochondrial DNA sequence divergence and applying molecular clock calibrations. Similar levels of genetic divergence were observed for Aplodactylus Valenciennes and Goniistius Gill species pairs, but these were more than twice that observed for a Nemadactylus Richardson pair. Molecular clock calibrations suggested divergences occurred during the late Miocene and mid Pliocene, respectively. Given evidence of high dispersal capabilities, the habitat and climatic barriers of the Australian south coast appear too small to have facilitated speciation of the cirrhitoids examined. A mechanism is proposed by which ancestral cirrhitoids were vicariantly isolated into east and west coast populations during periods of climate change. Although Aplodactylus and Goniistius divergences occurred during the same period, separate vicariant events across the Australian north and south coasts are invoked.     

Inferring clocks when lacking rocks: the variable rates of molecular evolution in bacteria

Background  

Because bacteria do not have a robust fossil record, attempts to infer the timing of events in their evolutionary history requires comparisons of molecular sequences. This use of molecular clocks is based on the assumptions that substitution rates for homologous genes or sites are fairly constant through time and across taxa. Violation of these conditions can lead to erroneous inferences and result in estimates that are off by orders of magnitude. In this study, we examine the consistency of substitution rates among a set of conserved genes in diverse bacterial lineages, and address the questions regarding the validity of molecular dating.     

Dating phototrophic microbial lineages with reticulate gene histories

Phototrophic bacteria are among the most biogeochemically significant organisms on Earth and are physiologically related through the use of reaction centers to collect photons for energy metabolism. However, the major phototrophic lineages are not closely related to one another in bacterial phylogeny, and the origins of their respective photosynthetic machinery remain obscured by time and low sequence similarity. To better understand the co‐evolution of Cyanobacteria and other ancient anoxygenic phototrophic lineages with respect to geologic time, we designed and implemented a variety of molecular clocks that use horizontal gene transfer (HGT) as additional, relative constraints. These HGT constraints improve the precision of phototroph divergence date estimates and indicate that stem green non‐sulfur bacteria are likely the oldest phototrophic lineage. Concurrently, crown Cyanobacteria age estimates ranged from 2.2 Ga to 2.7 Ga, with stem Cyanobacteria diverging ~2.8 Ga. These estimates provide a several hundred Ma window for oxygenic photosynthesis to evolve prior to the Great Oxidation Event (GOE) ~2.3 Ga. In all models, crown green sulfur bacteria diversify after the loss of the banded iron formations from the sedimentary record (~1.8 Ga) and may indicate the expansion of the lineage into a new ecological niche following the GOE. Our date estimates also provide a timeline to investigate the temporal feasibility of different photosystem HGT events between phototrophic lineages. Using this approach, we infer that stem Cyanobacteria are unlikely to be the recipient of an HGT of photosystem I proteins from green sulfur bacteria but could still have been either the HGT donor or the recipient of photosystem II proteins with green non‐sulfur bacteria, prior to the GOE. Together, these results indicate that HGT‐constrained molecular clocks are useful tools for the evaluation of various geological and evolutionary hypotheses, using the evolutionary histories of both genes and organismal lineages.     

Bayesian random local clocks,or one rate to rule them all

Background  

Relaxed molecular clock models allow divergence time dating and "relaxed phylogenetic" inference, in which a time tree is estimated in the face of unequal rates across lineages. We present a new method for relaxing the assumption of a strict molecular clock using Markov chain Monte Carlo to implement Bayesian modeling averaging over random local molecular clocks. The new method approaches the problem of rate variation among lineages by proposing a series of local molecular clocks, each extending over a subregion of the full phylogeny. Each branch in a phylogeny (subtending a clade) is a possible location for a change of rate from one local clock to a new one. Thus, including both the global molecular clock and the unconstrained model results, there are a total of 2^2n-2 possible rate models available for averaging with 1, 2, ..., 2n - 2 different rate categories.     

Of mice and men: Fossil-based divergence dates and molecular “clocks”

The magnitude of immunological differences between species has been used to estimate the time of divergence of lineages, especially among primates. Calibration of “clocks” based on immunological differences has relied heavily on data from Pakistan suggest that the lineages leading to Mus and Rattus diverged between 14 and 8 million years ago rather than 30 or more million years ago as suggested by techniques invoking constant rates of molecular evolution. Molecular analyses of pairs of mammal species which appear well documented paleontologically should be undertaken in order to expand our understanding of the tempo and mode of molecular evolution.     

Molecular clocks do not support the Cambrian explosion

The fossil record has long supported the view that most animal phyla originated during a brief period approximately 520 MYA known as the Cambrian explosion. However, molecular data analyses over the past 3 decades have found deeper divergences among animals (approximately 800 to 1,200 MYA), with and without the assumption of a global molecular clock. Recently, two studies have instead reported time estimates apparently consistent with the fossil record. Here, we demonstrate that methodological problems in these studies cast doubt on the accuracy and interpretations of the results obtained. In the study by Peterson et al., young time estimates were obtained because fossil calibrations were used as maximum limits rather than as minimum limits, and not because invertebrate calibrations were used. In the study by Aris-Brosou and Yang, young time estimates were obtained because of problems with rate models and other methods specific to the study, and not because Bayesian methods were used. This also led to many anomalous findings in their study, including a primate-rodent divergence at 320 MYA. With these results aside, molecular clocks continue to support a long period of animal evolution before the Cambrian explosion of fossils.