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’.     

Introduction and synthesis: Plant phylogeny and the origin of major biomes

Phylogenetic trees based upon DNA sequence data, when calibrated with a dimension of time, allow inference of: (i) the pattern of accumulation of lineages through time; (ii) the time of origin of monophyletic groups; (iii) when lineages arrived in different geographical areas; (iv) the time of origin of biome-specific morphologies. This gives a powerful new view of the history of biomes that in many cases is not provided by the incomplete plant fossil record. Dated plant phylogenies for angiosperm families such as Leguminoaceae (Fabaceae), Melastomataceae sensu stricto, Annonaceae and Rhamnaceae indicate that long-distance, transoceanic dispersal has played an important role in shaping their distributions, and that this can obscure any effect of tectonic history, previously assumed to have been the major cause of their biogeographic patterns. Dispersal from other continents has also been important in the assembly of the Amazonian rainforest flora and the Australian flora. Comparison of dated biogeographic patterns of plants and animals suggests that recent long-distance dispersal might be more prevalent in plants, which has major implications for community assembly and coevolution. Dated plant phylogenies also reveal the role of past environmental changes on the evolution of lineages in species-rich biomes, and show that recent Plio-Pleistocene diversification has contributed substantially to their current species richness. Because of the critical role of fossils and morphological characters in assigning ages to nodes in phylogenetic trees, future studies must include careful morphological consideration of fossils and their extant relatives in a phylogenetic context. Ideal study systems will be based upon DNA sequence data from multiple loci and multiple fossil calibrations. This allows cross-validation both of age estimates from different loci, and from different fossil calibrations. For a more complete view of biome history, future studies should emphasize full taxon sampling in ecologically important groups, and should focus on geographical areas for which few species-level phylogenies are available, such as tropical Africa and Asia. These studies are urgent because understanding the history of biomes can both inform conservation decisions, and help predict the effects of future environmental changes at a time when biodiversity is being impacted on an unprecedented scale.     

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.     

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.     

Phylogenetic relationships of the Cretaceous frog Beelzebufo from Madagascar and the placement of fossil constraints based on temporal and phylogenetic evidence

The placement of fossil calibrations is ideally based on the phylogenetic analysis of extinct taxa. Another source of information is the temporal variance for a given clade implied by a particular constraint when combined with other, well-supported calibrations. For example, the frog Beelzebufo ampinga from the Cretaceous of Madagascar has been hypothesized to be a crown-group member of the New World subfamily Ceratophryinae, which would support a Late Cretaceous connection with South America. However, phylogenetic analyses and molecular divergence time estimates based on other fossils do not support this placement. We derive a metric, Δt, to quantify temporal divergence among chronograms and find that errors resulting from mis-specified calibrations are localized when additional nodes throughout the tree are properly calibrated. The use of temporal information from molecular data can further assist in testing phylogenetic hypotheses regarding the placement of extinct taxa.     

Model-based multi-locus estimation of decapod phylogeny and divergence times

Phylogenetic relationships among all of the major decapod infraorders have never been estimated using molecular data, while morphological studies produce conflicting results. In the present study, the phylogenetic relationships among the decapod basal suborder Dendrobranchiata and all of the currently recognized decapod infraorders within the suborder Pleocyemata (Caridea, Stenopodidea, Achelata, Astacidea, Thalassinidea, Anomala, and Brachyura) were inferred using 16S mtDNA, 18S and 28S rRNA, and the histone H3 gene. Phylogenies were reconstructed using the model-based methods of maximum likelihood and Bayesian methods coupled with Markov Chain Monte Carlo inference. The phylogenies revealed that the seven infraorders are monophyletic, with high clade support values (bp>70; pP>0.95) under both methods. The two suborders also were recovered as monophyletic, but with weaker support (bp=70; pP=0.74). Although the nodal support values for infraordinal relationships were low (bp<50; pP<0.77) the Anomala and Brachyura were basal to the rest of the 'Reptantia' in both reconstructions and using Bayesian tree topology tests alternate morphology-based hypotheses were rejected (P<0.01). Newly developed multi-locus Bayesian and likelihood heuristic rate-smoothing methods to estimate divergence times were compared using eight fossil and geological calibrations. Estimated times revealed that the Decapoda originated earlier than 437MYA and that the radiation within the group occurred rapidly, with all of the major lineages present by 325MYA. Node time estimation under both approaches is severely affected by the number and phylogenetic distribution of the fossil calibrations chosen. For analyses incorporating fossils as fixed ages, more consistent results were obtained by using both shallow and deep or clade-related calibration points. Divergence time estimation using fossils as lower and upper limits performed well with as few as one upper limit and a single deep fossil lower limit calibration.     

Extant primitively segmented spiders have recently diversified from an ancient lineage

Living fossils are lineages that have retained plesiomorphic traits through long time periods. It is expected that such lineages have both originated and diversified long ago. Such expectations have recently been challenged in some textbook examples of living fossils, notably in extant cycads and coelacanths. Using a phylogenetic approach, we tested the patterns of the origin and diversification of liphistiid spiders, a clade of spiders considered to be living fossils due to their retention of arachnid plesiomorphies and their exclusive grouping in Mesothelae, an ancient clade sister to all modern spiders. Facilitated by original sampling throughout their Asian range, we here provide the phylogenetic framework necessary for reconstructing liphistiid biogeographic history. All phylogenetic analyses support the monophyly of Liphistiidae and of eight genera. As the fossil evidence supports a Carboniferous Euramerican origin of Mesothelae, our dating analyses postulate a long eastward over-land dispersal towards the Asian origin of Liphistiidae during the Palaeogene (39–58 Ma). Contrary to expectations, diversification within extant liphistiid genera is relatively recent, in the Neogene and Late Palaeogene (4–24 Ma). While no over-water dispersal events are needed to explain their evolutionary history, the history of liphistiid spiders has the potential to play prominently in vicariant biogeographic studies.     

Phylogenomic Resolution of Sea Spider Diversification through Integration of Multiple Data Classes

Despite significant advances in invertebrate phylogenomics over the past decade, the higher-level phylogeny of Pycnogonida (sea spiders) remains elusive. Due to the inaccessibility of some small-bodied lineages, few phylogenetic studies have sampled all sea spider families. Previous efforts based on a handful of genes have yielded unstable tree topologies. Here, we inferred the relationships of 89 sea spider species using targeted capture of the mitochondrial genome, 56 conserved exons, 101 ultraconserved elements, and 3 nuclear ribosomal genes. We inferred molecular divergence times by integrating morphological data for fossil species to calibrate 15 nodes in the arthropod tree of life. This integration of data classes resolved the basal topology of sea spiders with high support. The enigmatic family Austrodecidae was resolved as the sister group to the remaining Pycnogonida and the small-bodied family Rhynchothoracidae as the sister group of the robust-bodied family Pycnogonidae. Molecular divergence time estimation recovered a basal divergence of crown group sea spiders in the Ordovician. Comparison of diversification dynamics with other marine invertebrate taxa that originated in the Paleozoic suggests that sea spiders and some crustacean groups exhibit resilience to mass extinction episodes, relative to mollusk and echinoderm lineages.     

Relaxed molecular clocks for dating historical plant dispersal events

Age estimation from molecular sequences has emerged as a powerful tool for inferring when a plant lineage arrived in a particular area. Knowing the tenure of lineages within a region is key to understanding the evolution of traits, the evolution of biotic interactions, and the evolution of floras. New analytical methods model change in substitution rates along individual branches of a phylogenetic tree by combining molecular data with time constraints, usually from fossils. These "relaxed clock" approaches can be applied to several gene regions that need not all have the same substitution rates, and they can also incorporate multiple simultaneous fossil calibrations. Since 1995, at least 100 plant biogeographic studies have used molecular-clock dating, and about a fifth has used relaxed clocks. Many of these report evidence of long-distance dispersal. Meta-analyses of studies from the same geographic region can identify directional biases because of prevailing wind or water currents and the relative position and size of landmasses.     

Diversity-dependence brings molecular phylogenies closer to agreement with the fossil record

The branching times of molecular phylogenies allow us to infer speciation and extinction dynamics even when fossils are absent. Troublingly, phylogenetic approaches usually return estimates of zero extinction, conflicting with fossil evidence. Phylogenies and fossils do agree, however, that there are often limits to diversity. Here, we present a general approach to evaluate the likelihood of a phylogeny under a model that accommodates diversity-dependence and extinction. We find, by likelihood maximization, that extinction is estimated most precisely if the rate of increase in the number of lineages in the phylogeny saturates towards the present or first decreases and then increases. We demonstrate the utility and limits of our approach by applying it to the phylogenies for two cases where a fossil record exists (Cetacea and Cenozoic macroperforate planktonic foraminifera) and to three radiations lacking fossil evidence (Dendroica, Plethodon and Heliconius). We propose that the diversity-dependence model with extinction be used as the standard model for macro-evolutionary dynamics because of its biological realism and flexibility.     

A phylogeny of the fossil and extant zeiform-like fishes, Upper Cretaceous to Recent, with comments on the putative zeomorph clade (Acanthomorpha)

Tyler, J. C. & Santini, F. (2005). A phylogeny of the fossil and extant zeiform-like fishes, Upper Cretaceous to Recent, with comments on the putative zeomorph clade (Acanthomorpha). — Zoological Scripta , ** , ***–***.
A phylogenetic hypothesis based on 107 morphological characters is proposed for a data set of 43 taxa. Thirty-three are extant and belong to the orders Zeiformes (20 taxa), Caproiformes (2), Tetraodontiformes (2), Beryciformes (3), Stephanoberyciformes (3) and Perciformes (3). Ten are fossil taxa previously assigned to the Zeiformes (3), Caproiformes (1), Tetraodontiformes (2), Perciformes (1), and to two extinct Eocene families, the Sorbinipercidae (2) and the Zorzinichthyidae (1). This analysis indicates the existence of a previously undocumented clade formed by the families Sorbinipercidae + Zorzinichthyidae that may be related to the tetraodontiforms. It also shows that two uppermost Palaeocene species, Archaeozeus skamolensis and Protozeus kuehnei , sequentially represent the two most basal lineages of zeiforms, whereas the most ancient known zeiform, the Upper Cretaceous Cretazeus rinaldii , belongs within the clade of extant species in a polytomy with many other zeiform lineages. A reduced data set of 25 mostly zeiform taxa, after the removal of most outgroups, shows at least weak support for Cretazeus being nested deeply within the extant zeiforms; such a placement would indicate that at least six lineages of zeiforms were present during the Upper Cretaceous, and survived the Cretaceous/Tertiary (K/T) extinction to radiate in Cenozoic seas.     

Assessing concordance of fossil calibration points in molecular clock studies: an example using turtles

Although still controversial, estimation of divergence times using molecular data has emerged as a powerful tool to examine the tempo and mode of evolutionary change. Two primary obstacles in improving the accuracy of molecular dating are heterogeneity in DNA substitution rates and accuracy of the fossil record as calibration points. Recent methodological advances have provided powerful methods that estimate relative divergence times in the face of heterogeneity of nucleotide substitution rates among lineages. However, relatively little attention has focused on the accuracy of fossil calibration points that allow one to translate relative divergence times into absolute time. We present a new cross-validation method that identifies inconsistent fossils when multiple fossil calibrations are available for a clade and apply our method to a molecular phylogeny of living turtles with fossil calibration times for 17 of the 22 internal nodes in the tree. Our cross-validation procedure identified seven inconsistent fossils. Using the consistent fossils as calibration points, we found that despite their overall antiquity as a lineage, the most species-rich clades of turtles diversified well within the Cenozoic. Many of the truly ancient lineages of turtles are currently represented by a few, often endangered species that deserve high priority as conservation targets.     

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.     

Molecular phylogenetics of the spider infraorder Mygalomorphae using nuclear rRNA genes (18S and 28S): conflict and agreement with the current system of classification

Mygalomorph spiders, which include the tarantulas, trapdoor spiders, and their kin, represent one of three main spider lineages. Mygalomorphs are currently classified into 15 families, comprising roughly 2500 species and 300 genera. The few published phylogenies of mygalomorph relationships are based exclusively on morphological data and reveal areas of both conflict and congruence, suggesting the need for additional phylogenetic research utilizing new character systems. As part of a larger combined evidence study of global mygalomorph relationships, we have gathered approximately 3.7 kb of rRNA data (18S and 28S) for a sample of 80 genera, representing all 15 mygalomorph families. Taxon sampling was particularly intensive across families that are questionable in composition-Cyrtaucheniidae and Nemesiidae. The following primary results are supported by both Bayesian and parsimony analyses of combined matrices representing multiple 28S alignments: (1) the Atypoidea, a clade that includes the families Atypidae, Antrodiaetidae, and Mecicobothriidae, is recovered as a basal lineage sister to all other mygalomorphs, (2) diplurids and hexathelids form a paraphyletic grade at the base of the non-atypoid clade, but neither family is monophyletic in any of our analyses, (3) a clade consisting of all sampled nemesiids, Microstigmata and the cyrtaucheniid genera Kiama, Acontius, and Fufius is consistently recovered, (4) other sampled cyrtaucheniids are fragmented across three separate clades, including a monophyletic North American Euctenizinae and a South African clade, (5) of the Domiothelina, only idiopids are consistently recovered as monophyletic; ctenizids are polyphyletic and migids are only weakly supported. The Domiothelina is not monophyletic. The molecular results we present are consistent with more recent hypotheses of mygalomorph relationship; however, additional work remains before mygalomorph classification can be formally reassessed with confidence-increased taxonomic sampling and the inclusion of additional character systems (more genes and morphology) are required.     

Explosive radiation of Malpighiales supports a mid-cretaceous origin of modern tropical rain forests

Fossil data have been interpreted as indicating that Late Cretaceous tropical forests were open and dry adapted and that modern closed-canopy rain forest did not originate until after the Cretaceous-Tertiary (K/T) boundary. However, some mid-Cretaceous leaf floras have been interpreted as rain forest. Molecular divergence-time estimates within the clade Malpighiales, which constitute a large percentage of species in the shaded, shrub, and small tree layer in tropical rain forests worldwide, provide new tests of these hypotheses. We estimate that all 28 major lineages (i.e., traditionally recognized families) within this clade originated in tropical rain forest well before the Tertiary, mostly during the Albian and Cenomanian (112-94 Ma). Their rapid rise in the mid-Cretaceous may have resulted from the origin of adaptations to survive and reproduce under a closed forest canopy. This pattern may also be paralleled by other similarly diverse lineages and supports fossil indications that closed-canopy tropical rain forests existed well before the K/T boundary. This case illustrates that dated phylogenies can provide an important new source of evidence bearing on the timing of major environmental changes, which may be especially useful when fossil evidence is limited or controversial.     

Taking into account phylogenetic and divergence‐time uncertainty in a parametric biogeographical analysis of the Northern Hemisphere plant clade Caprifolieae

Aim To examine the biogeographical history of the angiosperm clade Caprifolieae (Caprifoliaceae) using parametric biogeographical reconstruction methods. The existing parametric method was extended to evaluate biogeographical reconstructions over the distribution of dated phylogenies. This method provides a framework for reconstructing large‐scale biogeography with parametric methods, while accounting for uncertainty in the phylogenetic relationships and divergence time. Location Asia, Europe and North America. Methods The biogeographical history of the major lineages of Caprifolieae was reconstructed over the posterior distribution of dated trees generated from Bayesian divergence‐time analyses. Results from a model with no geological constraints were compared with those from one where movement is disallowed across the North Atlantic after the Eocene. The most plausible scenarios were segregated at each node to test whether particular scenarios were reconstructed for particular divergence times. The parametric biogeographical method was also extended to estimate connectivity between areas so that the probability of dispersal between the major areas of the Northern Hemisphere could be explored. Results Phylogenetic results for Caprifolieae agreed with previous estimates using smaller sampling, but uncertainty remained despite efforts to resolve the relationships of the four genera within this clade using multiple markers. In addition to topological uncertainty, there were few fossils available for calibrations, resulting in large confidence intervals for divergence times. Divergence‐time analyses put the diversification of Caprifolieae at between 36 and 51 Ma and showed that Caprifolieae probably originated in Asia, with multiple movements into Europe and western and eastern North America. Main conclusions Newly developed parametric methods for biogeographical reconstruction incorporate more data and better models. Here, the parametric biogeographical reconstruction method has been extended to allow for topological and divergence‐time uncertainty. The analyses of Caprifolieae demonstrated that biogeographical hypotheses can be explored even where there are large confidence intervals on divergence times and uncertainty in topology. These results add to the growing evidence that Asia was an important source of Northern Hemisphere diversity throughout the Cenozoic.     

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.     

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.