Bees diversified in the age of eudicots

Reliable estimates on the ages of the major bee clades are needed to further understand the evolutionary history of bees and their close association with flowering plants. Divergence times have been estimated for a few groups of bees, but no study has yet provided estimates for all major bee lineages. To date the origin of bees and their major clades, we first perform a phylogenetic analysis of bees including representatives from every extant family, subfamily and almost all tribes, using sequence data from seven genes. We then use this phylogeny to place 14 time calibration points based on information from the fossil record for an uncorrelated relaxed clock divergence time analysis taking into account uncertainties in phylogenetic relationships and the fossil record. We explore the effect of placing a hard upper age bound near the root of the tree and the effect of different topologies on our divergence time estimates. We estimate that crown bees originated approximately 123 Ma (million years ago) (113–132 Ma), concurrently with the origin or diversification of the eudicots, a group comprising 75 per cent of angiosperm species. All of the major bee clades are estimated to have originated during the Middle to Late Cretaceous, which is when angiosperms became the dominant group of land plants.     

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.     

Combining fossil and molecular data to date the diversification of New World Primates

Recent methodological advances in molecular dating associated with the growing availability of sequence data have prompted the study of the evolution of New World Anthropoidea in recent years. Motivated by questions regarding historical biogeography or the mode of evolution, these works aimed to obtain a clearer scenario of Platyrrhini origins and diversification. Although some consensus was found, disputed issues, especially those relating to the evolutionary affinities of fossil taxa, remain. The use of fossil taxa for divergence time analysis is traditionally restricted to the provision of calibration priors. However, new analytical approaches have been developed that incorporate fossils as terminals and, thus, directly assign ages to the fossil tips. In this study, we conducted a combined analysis of molecular and morphological data, including fossils, to derive the timescale of New World anthropoids. Differently from previous studies that conducted total‐evidence analysis of molecules and morphology, our approach investigated the morphological clock alone. Our results corroborate the hypothesis that living platyrrhines diversified in the last 20 Ma and that Miocene Patagonian fossils compose an independent evolutionary radiation that diversified in the late Oligocene. When compared to the node ages inferred from the molecular timescale, the inclusion of fossils augmented the precision of the estimates for nodes constrained by the fossil tips. We show that morphological data can be analysed using the same methodological framework applied in relaxed molecular clock studies.     

Early divergence dates of demosponges based on mitogenomics and evaluated fossil calibrations

Demosponges are among the most primitive biomineralized metazoans to appear first in the fossil record with hard skeletons; their confirmed earliest fossils are from the lower Cambrian rocks about 520 Ma, with putative demosponge biomarkers reported from 713 to 635 Ma sediments. In this study, we use mitogenomic data to approach the early divergence timescale of demosponges using relaxed molecular clock techniques and likelihood-evaluated fossil calibration strategies. We found that among various molecular dating models, the correlated rate model yielded time estimates of demosponges in this analysis which is most congruent with the fossil appearance dates of demosponges. Our dating analyses show that crown groups of Demospongiae appeared at about 704 (674–741) Ma, and the silicification in demosponges (divergence of spicular sponges) began about 633 (616–648) Ma indicating a gap of over 100 million years between the origin of silicification and their first unequivocal appearance of siliceous spicules in the fossil record (520–525 Ma); demosponges with tetraxon-type spicules (Tetractinellida) are dated here at about 514 (498–530) Ma, an estimate comparable with the earliest tetraxial megasclere fossil records (510–520 Ma, Ordian Age, middle Cambrian).     

The impact of fossil stratigraphic ranges on tip-calibration,and the accuracy and precision of divergence time estimates

The molecular clock provides the only viable means of establishing realistic evolutionary timescales but it remains unclear how best to calibrate divergence time analyses. Calibrations can be applied to the tips and/or to the nodes of a phylogeny. Tip-calibration is an attractive approach since it allows fossil species to be included alongside extant relatives in molecular clock analyses. However, most fossil species are known from multiple stratigraphical horizons and it remains unclear how such age ranges should be interpreted to codify tip-calibrations. We use simulations and empirical data to explore the impact on precision and accuracy of different approaches to informing tip-calibrations. In particular, we focus on the effect of using tip-calibrations defined using the oldest vs youngest stratigraphic occurrences, the full stratigraphical range, as well as confidence intervals on these data points. The results of our simulations show that using different calibration approaches leads to different divergence-time estimates and demonstrate that concentrating tip-calibrations near the root of the dated phylogeny improves both precision and accuracy of estimated divergence times. Finally, our results indicate that the highest levels of accuracy and precision are achieved when fossil tips are calibrated based on the fossil occurrence from which the morphological data were derived. These trends were corroborated by analysis of an empirical dataset for Ursidae. Overall, we conclude that tip-dating analyses should, in particular, employ tip calibrations close to the root of the tree and they should be calibrated based on the age of the fossil used to inform the morphological data used in Total Evidence Dating.     

Molecular phylogenetic dating of asterid flowering plants shows early Cretaceous diversification

We present a phylogenetic dating of asterids, based on a 111-taxon tree representing all major groups and orders and 83 of the 102 families of asterids, with an underlying data set comprising six chloroplast DNA markers totaling 9914 positions. Phylogenetic dating was done with semiparametric rate smoothing by penalized likelihood. Confidence intervals were calculated by bootstrapping. Six reference fossils were used for calibration. To explore the effects of various sources of error, we repeated the analyses with alternative dating methods (nonparametric rate smoothing and the Langley-Fitch clock-based method), alternative tree topologies, reduced taxon sampling (22 of the 111 taxa deleted), partitioning the data into three genes and three noncoding regions, and calibrating with single reference fossils. The analyses with alternative topologies, reduced taxon sampling, and coding versus noncoding sequences all yielded small or in some cases no deviations. The choice of method influenced the age estimates of a few nodes considerably. Calibration with reference fossils is a critical issue, and use of single reference fossils yielded different results depending on the fossil. The bootstrap confidence intervals were generally small. Our results show that asterids and their major subgroups euasterids, campanulids, and lamiids diversified during the Early Cretaceous. Cornales, Ericales, and Aquifoliales also have crown node ages from the Early Cretaceous. Dipsacales and Solanales are from the Mid-Cretaceous, the other orders of core campanulids and core lamiids from the Late Cretaceous. The considerable diversity exhibited by asterids almost from their first appearance in the fossil record also supports an origin and first phase of diversification in the Early Cretaceous.     

Molecular and Paleontological Evidence for a Post-Cretaceous Origin of Rodents

The timing of the origin and diversification of rodents remains controversial, due to conflicting results from molecular clocks and paleontological data. The fossil record tends to support an early Cenozoic origin of crown-group rodents. In contrast, most molecular studies place the origin and initial diversification of crown-Rodentia deep in the Cretaceous, although some molecular analyses have recovered estimated divergence times that are more compatible with the fossil record. Here we attempt to resolve this conflict by carrying out a molecular clock investigation based on a nine-gene sequence dataset and a novel set of seven fossil constraints, including two new rodent records (the earliest known representatives of Cardiocraniinae and Dipodinae). Our results indicate that rodents originated around 61.7–62.4 Ma, shortly after the Cretaceous/Paleogene (K/Pg) boundary, and diversified at the intraordinal level around 57.7–58.9 Ma. These estimates are broadly consistent with the paleontological record, but challenge previous molecular studies that place the origin and early diversification of rodents in the Cretaceous. This study demonstrates that, with reliable fossil constraints, the incompatibility between paleontological and molecular estimates of rodent divergence times can be eliminated using currently available tools and genetic markers. Similar conflicts between molecular and paleontological evidence bedevil attempts to establish the origination times of other placental groups. The example of the present study suggests that more reliable fossil calibration points may represent the key to resolving these controversies.     

Four well-constrained calibration points from the vertebrate fossil record for molecular clock estimates

Recent controversy about the use of the vertebrate fossil record for external calibration of molecular clocks centers on two issues, the number of dates used for calibration and the reliability of the fossil calibration date. Viewing matters from a palaeontological perspective, we propose three qualitative, phylogenetic criteria that can be used within a comparative framework for the selection of well-constrained calibration dates from the vertebrate fossil record. On the basis of these criteria, we identify three highly suitable new fossil calibration dates for molecular clock estimates: the lungfish-tetrapod split (between 419 and 408 Ma), the bird-crocodile split (between 251 and 243 Ma), and the alligator-caiman split (between 71 and 66 Ma). Together with our previously suggested bird-lizard split (between 252 and 257 Ma), these four fossil dates span a range of very different ages. They are, in our opinion, more suitable for molecular clock calibration than the traditionally used mammal-bird split, which is less well constrained. We plea for closer interactions between paleontologists and molecular biologists in studying the timescale of vertebrate evolution.     

Evolution and biogeography of primates: a new model based on molecular phylogenetics,vicariance and plate tectonics

Heads, M. Evolution and biogeography of primates: a new model based on molecular phylogenetics, vicariance and plate tectonics. —Zoologica Scripta, 39, 107–127. The ages of the oldest fossils suggest an origin for primates in the Paleocene (～56 Ma). Fossil‐calibrated molecular clock dates give Cretaceous dates (～80–116 Ma). Both these estimates are minimum dates although they are often ‘transmogrified’ and treated as maximum or absolute dates. Oldest fossils can underestimate ages by tens of millions of years and instead of calibrating the time‐course of evolution with a scanty fossil record, the geographical boundaries of the main molecular clades of primates are calibrated here with radiometrically dated tectonic events. This indicates that primates originated when a globally widespread ancestor (early Archonta) differentiated into a northern group (Plesiadapiformes, extinct), a southern group (Primates), and two south‐east Asian groups (Dermoptera and Scandentia). The division occurred with the breakup of Pangea in the Early Jurassic and the opening of the central Atlantic (～185 Ma). Within primates, the strepsirrhines and haplorhines diverged with volcanism and buckling on the Lebombo Monocline, a volcanic rifted margin in south‐east Africa (Early Jurassic, ～180 Ma). Within strepsirrhines, lorises and galagos (Africa and Asia) and lemurs (Madagascar) diverged with the formation of the Mozambique Channel (Middle Jurassic, ～160 Ma). Within haplorhines, Old World monkeys and New World monkeys diverged with the opening of the Atlantic (Early Cretaceous, ～130 Ma). The main aspects of primate distribution are interpreted as the result of plate tectonics, phylogeny and vicariance, with some subsequent range expansion leading to secondary overlap. Long‐distance, trans‐oceanic dispersal events are not necessary. The primate ancestral complex was already widespread globally when sea‐floor spreading, strike‐slip rifting and orogeny fractured and deformed distributions through the Jurassic and Cretaceous, leading to the origin of the modern clades. The model suggests that the topology of the phylogenetic tree reflects a sequence of differentiation in a widespread ancestor rather than a series of dispersal events.     

Revisiting taxonomy, morphological evolution, and fossil calibration strategies in Chloranthaceae

Chtoranthaceae is one of the earliest diverging angiosperm families and is comprised of approximately 75 species in four genera (Chloranthus,Sarcandra,Ascarina,and Hedyosmurn).This family has received considerable attention because of its seemingly primitive morphology,disjunct tropical distribution in Asia and America,and extensive fossil record from the Early Cretaceous.In the present study,we reconstructed the phylogeny of Chloranthaceae based on a combined dataset of three plastid DNA regions and 56 species.We then estimated divergence times in the family using two relaxed molecular clock methods (BEAST and penalized likelihood).We focused on testing the influence of fossil taxa in calibrating the molecular phylogeny,and on assessing the current taxonomy of the family in light of the phylogenetic results.Our results indicate that most intrageneric divisions within Ascarina and Hedyosmum are not monophyletic.The results from the dating analysis suggest that the Hedyosmum-like fossil Asteropollis represents a stem lineage of Hedyosmum,as has been suggested previously from morphological analyses.In contrast,our results indicate that the Chloranthus-like fossil Chloranthistemon,previously suggested on morphological grounds to be a stem relative of Chloranthus,may,instead,belong to the branch leading to the clade Chloranthus + Sarcandra.The median crown ages of Chloranthus,Sarcandra,Ascarina,and Hedyosmum estimated in the BEAST analysis were 26.3,9.5,31.0 and 45.8 million years ago (Ma),respectively,whereas the divergence between Chloranthus and Sarcandra,the splitting of Ascarina with the former two genera,and Hedyosmum separating from the three genera were estimated to 63.8,95.7 and 111.1 Ma.The present study sheds further light on the temporal evolution of Chloranthaceae and exemplifies how molecular dating analyses may be used to explore alternative phylogenetic placements of fossil taxa.     

Angiosperm divergence times: the effect of genes, codon positions, and time constraints

An understanding of the evolution of modern terrestrial ecosystems requires an understanding of the dynamics associated with angiosperm evolution, including the timing of their origin and diversification into their extraordinary present-day diversity. Molecular estimates of angiosperm age have varied widely, and many substantially predate the Early Cretaceous fossil appearance of the group. In this study, the effect of different genes, codon positions, and chronological constraints on node ages are examined on divergence time estimates across seed plants, with a special focus on angiosperms. Penalized likelihood was used to estimate divergence times on a phylogenetic hypothesis for seed plants derived from Bayesian analysis, with branch lengths estimated with maximum likelihood. The plastid genes atpB, psaA, psbB, and rbcL were used individually and in combination, using first and second, third, and the three codon positions, including and excluding age constraints on 20 nodes derived from a critical examination of the land-plant fossil record. The optimal level of rate smoothing according to each unconstrained and constrained dataset was obtained with penalized likelihood. Tests for a molecular clock revealed significantly unclocklike rates in all datasets. Addition of fossil constraints resulted in even greater departures from constancy. Consistently with significant deviations from a clock, estimated optimal smoothing values were low, but a strict correlation between rate heterogeneity and optimal smoothing value was not found. Age estimates for nodes across the phylogeny varied, sometimes substantially, with gene and codon position. Nevertheless, estimates based on the four concatenated genes are very similar to the mean of the four individual gene estimates. For any given node, unconstrained age estimates are more variable than constrained estimates and are frequently younger than well-substantiated fossil members of the clade. Constrained estimates of ages of clades are older than unconstrained estimates and oldest fossil representatives, sometimes substantially so. Angiosperm age estimates decreased as rate smoothing increased. Whereas the range of unconstrained angiosperm age estimates spans the fossil age of the clade, the range of constrained estimates is narrower (and older) than the earliest angiosperm fossils. Results unambiguously indicate the relevance of constraints in reducing the variability of ages derived from different partitions of the data and diminishing the effect of the smoothing parameter. Constrained optimizations of divergence times and substitution rates across the phylogeny suggest appreciably different evolutionary dynamics for angiosperms and for gymnosperms. Whereas the gymnosperm crown group originated shortly after the origin of seed plants, a long time elapsed before the origin of crown group angiosperms. Although absolute age estimates of angiosperms and angiosperm clades are older than their earliest fossils, the estimated pace of phylogenetic diversification largely agrees with the rapid appearance of angiosperm lineages in stratigraphic sequences.     

Molecular and fossil evidence place the origin of cichlid fishes long after Gondwanan rifting

Cichlid fishes are a key model system in the study of adaptive radiation, speciation and evolutionary developmental biology. More than 1600 cichlid species inhabit freshwater and marginal marine environments across several southern landmasses. This distributional pattern, combined with parallels between cichlid phylogeny and sequences of Mesozoic continental rifting, has led to the widely accepted hypothesis that cichlids are an ancient group whose major biogeographic patterns arose from Gondwanan vicariance. Although the Early Cretaceous (ca 135 Ma) divergence of living cichlids demanded by the vicariance model now represents a key calibration for teleost molecular clocks, this putative split pre-dates the oldest cichlid fossils by nearly 90 Myr. Here, we provide independent palaeontological and relaxed-molecular-clock estimates for the time of cichlid origin that collectively reject the antiquity of the group required by the Gondwanan vicariance scenario. The distribution of cichlid fossil horizons, the age of stratigraphically consistent outgroup lineages to cichlids and relaxed-clock analysis of a DNA sequence dataset consisting of 10 nuclear genes all deliver overlapping estimates for crown cichlid origin centred on the Palaeocene (ca 65–57 Ma), substantially post-dating the tectonic fragmentation of Gondwana. Our results provide a revised macroevolutionary time scale for cichlids, imply a role for dispersal in generating the observed geographical distribution of this important model clade and add to a growing debate that questions the dominance of the vicariance paradigm of historical biogeography.     

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.     

Dating an impressive Neotropical radiation: Molecular time estimates for the Sigmodontinae (Rodentia) provide insights into its historical biogeography

With about 400 living species and 82 genera, rodents of the subfamily Sigmodontinae comprise one of the most diverse and more broadly distributed Neotropical mammalian clades. There has been much debate on the origin of the lineage or the lineages of sigmodontines that entered South America, the timing of entrance and different aspects of further diversification within South America. The ages of divergence of the main lineages and the crown age of the subfamily were estimated by using sequences of the interphotoreceptor retinoid binding protein and cytochrome b genes for a dense sigmodontine and muroid sampling. Bayesian inference using three fossil calibration points and a relaxed molecular clock estimated a middle Miocene origin for Sigmodontinae (～12 Ma), with most tribes diversifying throughout the Late Miocene (6.9–9.4 Ma). These estimates together results of analyses of ancestral area reconstructions suggest a distribution for the most recent common ancestor of Sigmodontinae in Central-South America and a South American distribution for the most recent common ancestor of Oryzomyalia.     

Integrating fossils in a molecular‐based phylogeny and testing them as calibration points for divergence time estimates in Menispermaceae

Abstract The phylogeny of extant Menispermaceae (Ranunculales) is reconstructed based on DNA sequences of two chloroplast genes (rbcL and atpB) from 94 species belonging to 56 genera. Fossilized endocarps represent 34 genera. The positions of these are inferred using 30 morphological characters and the molecular phylogeny as a backbone constraint. Nine of the thirteen nodes that are each dated by a fossil are used as calibration points for the estimates of molecular divergence times. BEAST is used to estimate stem age (121.2 Myr) and crown age (105.4 Myr) for Menispermaceae. This method does not require an input tree topology and can also account for rate heterogeneity among lineages. The sensitivity of these estimates to fossil constraints is then evaluated by a cross‐validation procedure. The estimated origin for Menispermaceae is dated to the mid‐Jurassic if the customary maximum age of 125 Myr for eudicots is not implemented. All constraints when used alone failed to estimate node ages in some parts of the tree. Fossils from the Palaeocene and Eocene impose strict constraints. Likewise, the use of Prototinomiscium as a dating constraint for Menispermaceae appears to be a conservative approach.     

Of teeth and trees: A fossil tip‐dating approach to infer divergence times of extinct and extant squaliform sharks

Fossil tip‐dating allows for the inclusion of morphological data in divergence time estimates based on both extant and extinct taxa. Neoselachii have a cartilaginous skeleton, which is less prone to fossilization compared to skeletons of Osteichthyans. Therefore, the majority of the neoselachian fossil record is comprised of single teeth, which fossilize more easily. Neoselachian teeth can be found in large numbers as they are continuously replaced. Tooth morphologies are of major importance on multiple taxonomic levels for identification of shark and ray taxa. Here, we review dental morphological characters of squalomorph sharks and test these for their phylogenetic signal. Subsequently, we combine DNA sequence data (concatenated exon sequences) with dental morphological characters from 85 fossil and extant taxa to simultaneously infer the phylogeny and re‐estimate divergence times using information of 61 fossil tip‐dates as well as eight node age calibrations of squalomorph sharks. Our findings show that the phylogenetic placement of fossil taxa is mostly in accordance with their previous taxonomic allocation. An exception is the phylogenetic placement of the extinct genus ?Protospinax , which remains unclear. We conclude that the high number of fossil taxa as well as the comprehensive DNA sequence data for extant taxa may compensate for the limited number of morphological characters identifiable on teeth, serving as a backbone for reliably estimating the phylogeny of both extinct and extant taxa. In general, tip‐dating mostly estimates older node ages compared to previous studies based on calibrated molecular clocks.     

Dating the origin of the major lineages of Branchiopoda

Despite the well-established phylogeny and good fossil record of branchiopods, a consistent macro-evolutionary timescale for the group remains elusive. This study focuses on the early branchiopod divergence dates where fossil record is extremely fragmentary or missing. On the basis of a large genomic dataset and carefully evaluated fossil calibration points, we assess the quality of the branchiopod fossil record by calibrating the tree against well-established first occurrences, providing paleontological estimates of divergence times and completeness of their fossil record. The maximum age constraints were set using a quantitative approach of Marshall (2008). We tested the alternative placements of Yicaris and Wujicaris in the referred arthropod tree via the likelihood checkpoints method. Divergence dates were calculated using Bayesian relaxed molecular clock and penalized likelihood methods. Our results show that the stem group of Branchiopoda is rooted in the late Neoproterozoic (563 ± 7 Ma); the crown-Branchiopoda diverged during middle Cambrian to Early Ordovician (478–512 Ma), likely representing the origin of the freshwater biota; the Phyllopoda clade diverged during Ordovician (448–480 Ma) and Diplostraca during Late Ordovician to early Silurian (430–457 Ma). By evaluating the congruence between the observed times of appearance of clade in the fossil record and the results derived from molecular data, we found that the uncorrelated rate model gave more congruent results for shallower divergence events whereas the auto-correlated rate model gives more congruent results for deeper events.     

Assessing the sensitivity of divergence time estimates to locus sampling,calibration points,and model priors in a RAD-seq phylogeny of Carex section Schoenoxiphium

Restriction site-associated DNA sequencing (RAD-seq) and related methods have become relatively common approaches to resolve species-level phylogeny. It is not clear, however, whether RAD-seq data matrices are well suited to relaxed clock inference of divergence times, given the size of the matrices and the abundance of missing data. We investigated the sensitivity of Bayesian relaxed clock estimates of divergence times to alternative analytical decisions on an empirical RAD-seq phylogenetic matrix. We explored the relative contribution of secondary calibration strategies, amount of missing data, and the data partition analyzed to overall variance in divergence times inferred using BEAST MCMC analyses of Carex section Schoenoxiphium (Cyperaceae)—a recent radiation for which we have nearly complete species sampling of RAD-seq data. The crown node for Schoenoxiphium was estimated to be 15.22 (9.56–21.18) Ma using a single calibration point and low missing data, 11.93 (8.07–16.03) Ma using multiple calibration points and low missing data, and 8.34 (5.41–11.22) using multiple calibrations but high missing data. We found that using matrices with more than half of the individuals with missing data inferred younger mean ages for all nodes. Moreover, we have found that our molecular clock estimates are sensitive to the positions of the calibration(s) in our phylogenetic tree (using matrices with low missing data), especially when only a single calibration was applied to estimate divergence times. These results argue for sensitivity analyses and caution in interpreting divergence time estimates from RAD-seq data.     

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.