Divergence time estimation using fossils as terminal taxa and the origins of Lissamphibia

Were molecular data available for extinct taxa, questions regarding the origins of many groups could be settled in short order. As this is not the case, various strategies have been proposed to combine paleontological and neontological data sets. The use of fossil dates as node age calibrations for divergence time estimation from molecular phylogenies is commonplace. In addition, simulations suggest that the addition of morphological data from extinct taxa may improve phylogenetic estimation when combined with molecular data for extant species, and some studies have merged morphological and molecular data to estimate combined evidence phylogenies containing both extinct and extant taxa. However, few, if any, studies have attempted to estimate divergence times using phylogenies containing both fossil and living taxa sampled for both molecular and morphological data. Here, I infer both the phylogeny and the time of origin for Lissamphibia and a number of stem tetrapods using Bayesian methods based on a data set containing morphological data for extinct taxa, molecular data for extant taxa, and molecular and morphological data for a subset of extant taxa. The results suggest that Lissamphibia is monophyletic, nested within Lepospondyli, and originated in the late Carboniferous at the earliest. This research illustrates potential pitfalls for the use of fossils as post hoc age constraints on internal nodes and highlights the importance of explicit phylogenetic analysis of extinct taxa. These results suggest that the application of fossils as minima or maxima on molecular phylogenies should be supplemented or supplanted by combined evidence analyses whenever possible.     

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.     

Closing the gap between rocks and clocks using total-evidence dating

Total-evidence dating (TED) allows evolutionary biologists to incorporate a wide range of dating information into a unified statistical analysis. One might expect this to improve the agreement between rocks and clocks but this is not necessarily the case. We explore the reasons for such discordance using a mammalian dataset with rich molecular, morphological and fossil information. There is strong conflict in this dataset between morphology and molecules under standard stochastic models. This causes TED to push divergence events back in time when using inadequate models or vague priors, a phenomenon we term ‘deep root attraction’ (DRA). We identify several causes of DRA. Failure to account for diversified sampling results in dramatic DRA, but this can be addressed using existing techniques. Inadequate morphological models also appear to be a major contributor to DRA. The major reason seems to be that current models do not account for dependencies among morphological characters, causing distorted topology and branch length estimates. This is particularly problematic for huge morphological datasets, which may contain large numbers of correlated characters. Finally, diversification and fossil sampling priors that do not incorporate all the available background information can contribute to DRA, but these priors can also be used to compensate for DRA. Specifically, we show that DRA in the mammalian dataset can be addressed by introducing a modest extra penalty for ghost lineages that are unobserved in the fossil record, for instance by assuming rapid diversification, rare extinction or high fossil sampling rate; any of these assumptions produces highly congruent divergence time estimates with a minimal gap between rocks and clocks. Under these conditions, fossils have a stabilizing influence on divergence time estimates and significantly increase the precision of those estimates, which are generally close to the dates suggested by palaeontologists.This article is part of the themed issue ‘Dating species divergences using rocks and clocks’.     

Divergence Times and the Evolutionary Radiation of New World Monkeys (Platyrrhini,Primates): An Analysis of Fossil and Molecular Data

The estimation of phylogenetic relationships and divergence times among a group of organisms is a fundamental first step toward understanding its biological diversification. The time of the most recent or last common ancestor (LCA) of extant platyrrhines is one of the most controversial among scholars of primate evolution. Here we use two molecular based approaches to date the initial divergence of the platyrrhine clade, Bayesian estimations under a relaxed-clock model and substitution rate plus generation time and body size, employing the fossil record and genome datasets. We also explore the robustness of our estimations with respect to changes in topology, fossil constraints and substitution rate, and discuss the implications of our findings for understanding the platyrrhine radiation. Our results suggest that fossil constraints, topology and substitution rate have an important influence on our divergence time estimates. Bayesian estimates using conservative but realistic fossil constraints suggest that the LCA of extant platyrrhines existed at ca. 29 Ma, with the 95% confidence limit for the node ranging from 27–31 Ma. The LCA of extant platyrrhine monkeys based on substitution rate corrected by generation time and body size was established between 21–29 Ma. The estimates based on the two approaches used in this study recalibrate the ages of the major platyrrhine clades and corroborate the hypothesis that they constitute very old lineages. These results can help reconcile several controversial points concerning the affinities of key early Miocene fossils that have arisen among paleontologists and molecular systematists. However, they cannot resolve the controversy of whether these fossil species truly belong to the extant lineages or to a stem platyrrhine clade. That question can only be resolved by morphology. Finally, we show that the use of different approaches and well supported fossil information gives a more robust divergence time estimate of a clade.     

Phylogeny of extant and fossil Juglandaceae inferred from the integration of molecular and morphological data sets

It is widely acknowledged that integrating fossils into data sets of extant taxa is imperative for proper placement of fossils, resolution of relationships, and a better understanding of character evolution. The importance of this process has been further magnified because of the crucial role of fossils in dating divergence times. Outstanding issues remain, including appropriate methods to place fossils in phylogenetic trees, the importance of molecules versus morphology in these analyses, as well as the impact of potentially large amounts of missing data for fossil taxa. In this study we used the angiosperm clade Juglandaceae as a model for investigating methods of integrating fossils into a phylogenetic framework of extant taxa. The clade has a rich fossil record relative to low extant diversity, as well as a robust molecular phylogeny and morphological database for extant taxa. After combining fossil organ genera into composite and terminal taxa, our objectives were to (1) compare multiple methods for the integration of the fossils and extant taxa (including total evidence, molecular scaffolds, and molecular matrix representation with parsimony [MRP]); (2) explore the impact of missing data (incomplete taxa and characters) and the evidence for placing fossils on the topology; (3) simulate the phylogenetic effect of missing data by creating "artificial fossils"; and (4) place fossils and compare the impact of single and multiple fossil constraints in estimating the age of clades. Despite large and variable amounts of missing data, each of the methods provided reasonable placement of both fossils and simulated "artificial fossils" in the phylogeny previously inferred only from extant taxa. Our results clearly show that the amount of missing data in any given taxon is not by itself an operational guideline for excluding fossils from analysis. Three fossil taxa (Cruciptera simsonii, Paleoplatycarya wingii, and Platycarya americana) were placed within crown clades containing living taxa for which relationships previously had been suggested based on morphology, whereas Polyptera manningii, a mosaic taxon with equivocal affinities, was placed firmly as sister to two modern crown clades. The position of Paleooreomunnea stoneana was ambiguous with total evidence but conclusive with DNA scaffolds and MRP. There was less disturbance of relationships among extant taxa using a total evidence approach, and the DNA scaffold approach did not provide improved resolution or internal support for clades compared to total evidence, whereas weighted MRP retained comparable levels of support but lost crown clade resolution. Multiple internal minimum age constraints generally provided reasonable age estimates, but the use of single constraints provided by extinct genera tended to underestimate clade ages.     

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.     

Molecular dating of the 'Gondwanan' plant family Proteaceae is only partially congruent with the timing of the break-up of Gondwana

Aim  The flowering plant family Proteaceae is putatively of Gondwanan age, with modern and fossil lineages found on all southern continents. Here we test whether the present distribution of Proteaceae can be explained by vicariance caused by the break-up of Gondwana.
Location  Africa, especially southern Africa, Australia, New Zealand, South America, New Caledonia, New Guinea, Southeast Asia, Sulawesi, Tasmania.
Methods  We obtained chloroplast DNA sequence data from the rbc L gene, the rbc L- atp B spacer, and the atp B gene from leaf samples of forty-five genera collected from the field and from living collections. We analysed these data using Bayesian phylogenetic and molecular dating methods, with five carefully selected fossil calibration points to obtain age estimates for the nodes within the family.
Results  Four of eight trans-continental disjunctions of sister groups within our sample of the Proteaceae post-date the break-up of Gondwana. These involve independent lineages, two with an Africa-Australia disjunction, one with an Africa–South America disjunction, and one with a New Zealand–Australasia disjunction. The date of the radiation of the bird-pollinated Embothriinae corresponds approximately to the hypothesized date of origin of nectar-feeding birds in Australia.
Main conclusions  The findings suggest that disjunct distributions in Proteaceae result from both Gondwanan vicariance and transoceanic dispersal. Our results imply that ancestors of some taxa dispersed across oceans rather than rafting with Gondwanan fragments as previously thought. This finding agrees with other studies of Gondwanan plants in dating the divergence of Australian, New Zealand and New Caledonian taxa in the Eocene, consistent with the existence of a shared, ancestral Eocene flora but contrary to a vicariance scenario based on accepted geological knowledge.     

Dated Phylogenies of the Sister Genera Macaranga and Mallotus (Euphorbiaceae): Congruence in Historical Biogeographic Patterns?

Molecular phylogenies and estimates of divergence times within the sister genera Macaranga and Mallotus were estimated using Bayesian relaxed clock analyses of two generic data sets, one per genus. Both data sets were based on different molecular markers and largely different samples. Per genus three calibration points were utilised. The basal calibration point (crown node of all taxa used) was taken from literature and used for both taxa. The other three calibrations were based on fossils of which two were used per genus. We compared patterns of dispersal and diversification in Macaranga and Mallotus using ancestral area reconstruction in RASP (S-DIVA option) and contrasted our results with biogeographical and geological records to assess accuracy of inferred age estimates. A check of the fossil calibration point showed that the Japanese fossil, used for dating the divergence of Mallotus, probably had to be attached to a lower node, the stem node of all pioneer species, but even then the divergence time was still younger than the estimated age of the fossil. The African (only used in the Macaranga data set) and New Zealand fossils (used for both genera) seemed reliably placed. Our results are in line with existing geological data and the presence of stepping stones that provided dispersal pathways from Borneo to New Guinea-Australia, from Borneo to mainland Asia and additionally at least once to Africa and Madagascar via land and back to India via Indian Ocean island chains. The two genera show congruence in dispersal patterns, which corroborate divergence time estimates, although the overall mode and tempo of dispersal and diversification differ significantly as shown by distribution patterns of extant species.     

Morphology, fossils, divergence timing, and the phylogenetic relationships of Gavialis

Although morphological data have historically favored a basal position for the Indian gharial (Gavialis gangeticus) within Crocodylia and a Mesozoic divergence between Gavialis and all other crocodylians, several recent molecular data sets have argued for a sister-group relationship between Gavialis and the Indonesian false gharial (Tomistoma schlegelii) and a divergence between them no earlier than the Late Tertiary. Fossils were added to a matrix of 164 discrete morphological characters and subjected to parsimony analysis. When morphology was analyzed alone, Gavialis was the sister taxon of all other extant crocodylians whether or not fossil ingroup taxa were included, and a sister-group relationship between Gavialis and Tomistoma was significantly less parsimonious. In combination with published sequence and restriction site fragment data, Gavialis was the sister taxon of all other living crocodylians, but the position of Tomistoma depended on the inclusion of fossil ingroup taxa; with or without fossils, preferred morphological and molecular topologies were not significantly different. Fossils closer to Gavialis than to Tomistoma can be recognized in the Late Cretaceous, and fossil relatives of Tomistoma are known from the basal Eocene, strongly indicating a divergence long before the Late Tertiary. Comparison of minimum divergence time from the fossil record with different measures of molecular distance indicates evolutionary rate heterogeneity within Crocodylia. Fossils strongly contradict a post-Oligocene divergence between Gavialis and any other living crocodylian, but the phylogenetic placement of Gavialis is best viewed as unresolved.     

The impact of fossils and taxon sampling on ancient molecular dating analyses

The number and complexity of molecular dating studies has increased over the past decade. Along with a broadening acceptance of their utility has come significant controversy over the methods and models that are appropriate, as well as the accuracy of the estimates yielded by molecular clock analyses. Radically different age estimates have been published for the same divergences from analyses of different datasets with different fossil constraints obtained with different methods, and the underlying explanation for these differences is often unclear. Here we utilize two previously published datasets to examine the effect of fossil calibrations and taxon sampling on the age estimates for two deep eukaryote divergences in an attempt to discern the relative impact of these factors. Penalized likelihood, non-parametric rate smoothing, and Bayesian methods were utilized to generate age estimates for the origin of the Metazoa from a 7-gene dataset and for the divergence of Eukaryotes from a 129-gene dataset. From these analyses, it is clear that the fossil calibrations chosen and the method for applying constraints to these nodes have a large impact on age estimates, while the degree of taxon sampling within a dataset is less important in terms of the resulting age estimates. Concerns and recommendations for addressing these two factors when initiating a dating analysis are discussed.     

Back to Gondwanaland: can ancient vicariance explain (some) Indian Ocean disjunct plant distributions?

Oceans, or other wide expanses of inhospitable environment, interrupt present day distributions of many plant groups. Using molecular dating techniques, generally incorporating fossil evidence, we can estimate when such distributions originated. Numerous dating analyses have recently precipitated a paradigm shift in the general explanations for the phenomenon, away from older geological causes, such as continental drift, in favour of more recent, long-distance dispersal (LDD). For example, the ‘Gondwanan vicariance’ scenario has been dismissed in various studies of Indian Ocean disjunct distributions. We used the gentian tribe Exaceae to reassess this scenario using molecular dating with minimum (fossil), maximum (geological), secondary (from wider analyses) and hypothesis-driven age constraints. Our results indicate that ancient vicariance cannot be ruled out as an explanation for the early origins of Exaceae across Africa, Madagascar and the Indian subcontinent unless a strong assumption is made about the maximum age of Gentianales. However, both the Gondwanan scenario and the available evidence suggest that there were also several, more recent, intercontinental dispersals during the diversification of the group.     

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.     

Incorporating phylogenetic uncertainty on phylogeny‐based palaeontological dating and the timing of turtle diversification

Methods improving the performance of molecular dating of divergence time of clades have improved dramatically in recent years. The calibration of molecular dating using the first appearance of a clade in the fossil record is a crucial step towards inferring the minimal diversification time of various groups and the choice of extinct taxa can strongly influence the molecular dates. Here, we evaluate the uncertainty on the phylogenetic position of extinct taxa through non‐parametric bootstrapping. The recognition of phylogenetic uncertainty resulted in the definition of the Bootstrap Uncertainty Range (BUR) for the age of first appearance of a given clade. The BUR is calculated as the interval of geological time in which the diversification of a given clade can be inferred to have occurred, based on the temporal information of the fossil record and the topologies of the bootstrap trees. Divergence times based on BUR analyses were calculated for three clades of turtles: Testudines, Pleurodira and Cryptodira. This resulted in extensive uncertainty ranges of topology‐dependent minimal divergence dates for these clades.     

Systematics and evolution of the Pan‐Alcidae (Aves,Charadriiformes)

Puffins, auks and their allies in the wing‐propelled diving seabird clade Pan‐Alcidae (Charadriiformes) have been proposed to be key pelagic indicators of faunal shifts in Northern Hemisphere oceans. However, most previous phylogenetic analyses of the clade have focused only on the 23 extant alcid species. Here we undertake a combined phylogenetic analysis of all previously published molecular sequence data (～ 12 kb) and morphological data (n = 353 characters) with dense species level sampling that also includes 28 extinct taxa. We present a new estimate of the patterns of diversification in the clade based on divergence time estimates that include a previously vetted set of twelve fossil calibrations. The resultant time trees are also used in the evaluation of previously hypothesized paleoclimatic drivers of pan‐alcid evolution. Our divergence dating results estimate the split of Alcidae from its sister taxon Stercorariidae during the late Eocene (～ 35 Ma), an evolutionary hypothesis for clade origination that agrees with the fossil record and that does not require the inference of extensive ghost lineages. The extant dovekie Alle alle is identified as the sole extant member of a clade including four extinct Miocene species. Furthermore, whereas an Uria + Alle clade has been previously recovered from molecular analyses, the extinct diversity of closely related Miocepphus species yields morphological support for this clade. Our results suggest that extant alcid diversity is a function of Miocene diversification and differential extinction at the Pliocene–Pleistocene boundary. The relative timing of the Middle Miocene climatic optimum and the Pliocene–Pleistocene climatic transition and major diversification and extinction events in Pan‐Alcidae, respectively, are consistent with a potential link between major paleoclimatic events and pan‐alcid cladogenesis.     

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.     

The origin and diversification of eukaryotes: problems with molecular phylogenetics and molecular clock estimation

Determining the relationships among and divergence times for the major eukaryotic lineages remains one of the most important and controversial outstanding problems in evolutionary biology. The sequencing and phylogenetic analyses of ribosomal RNA (rRNA) genes led to the first nearly comprehensive phylogenies of eukaryotes in the late 1980s, and supported a view where cellular complexity was acquired during the divergence of extant unicellular eukaryote lineages. More recently, however, refinements in analytical methods coupled with the availability of many additional genes for phylogenetic analysis showed that much of the deep structure of early rRNA trees was artefactual. Recent phylogenetic analyses of a multiple genes and the discovery of important molecular and ultrastructural phylogenetic characters have resolved eukaryotic diversity into six major hypothetical groups. Yet relationships among these groups remain poorly understood because of saturation of sequence changes on the billion-year time-scale, possible rapid radiations of major lineages, phylogenetic artefacts and endosymbiotic or lateral gene transfer among eukaryotes. Estimating the divergence dates between the major eukaryote lineages using molecular analyses is even more difficult than phylogenetic estimation. Error in such analyses comes from a myriad of sources including: (i) calibration fossil dates, (ii) the assumed phylogenetic tree, (iii) the nucleotide or amino acid substitution model, (iv) substitution number (branch length) estimates, (v) the model of how rates of evolution change over the tree, (vi) error inherent in the time estimates for a given model and (vii) how multiple gene data are treated. By reanalysing datasets from recently published molecular clock studies, we show that when errors from these various sources are properly accounted for, the confidence intervals on inferred dates can be very large. Furthermore, estimated dates of divergence vary hugely depending on the methods used and their assumptions. Accurate dating of divergence times among the major eukaryote lineages will require a robust tree of eukaryotes, a much richer Proterozoic fossil record of microbial eukaryotes assignable to extant groups for calibration, more sophisticated relaxed molecular clock methods and many more genes sampled from the full diversity of microbial eukaryotes.     

The evolution of hominoid cranial diversity: A quantitative genetic approach

Hominoid cranial evolution is characterized by substantial phenotypic diversity, yet the cause of this variability has rarely been explored. Quantitative genetic techniques for investigating evolutionary processes underlying morphological divergence are dependent on the availability of good ancestral models, a problem in hominoids where the fossil record is fragmentary and poorly understood. Here, we use a maximum likelihood approach based on a Brownian motion model of evolutionary change to estimate nested hypothetical ancestral forms from 15 extant hominoid taxa. These ancestors were then used to calculate rates of evolution along each branch of a phylogenetic tree using Lande's generalized genetic distance. Our results show that hominoid cranial evolution is characterized by strong stabilizing selection. Only two instances of directional selection were detected; the divergence of Homo from its last common ancestor with Pan, and the divergence of the lesser apes from their last common ancestor with the great apes. In these two cases, selection gradients reconstructed to identify the specific traits undergoing selection indicated that selection on basicranial flexion, cranial vault expansion, and facial retraction characterizes the divergence of Homo, whereas the divergence of the lesser apes was defined by selection on neurocranial size reduction.     

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.