Calibration of galliform molecular clocks using multiple fossils and genetic partitions

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

Calibration of avian molecular clocks

Molecular clocks can be calibrated using fossils within the group under study (internal calibration) or outside of the group (external calibration). Both types of calibration have their advantages and disadvantages. An internal calibration may reduce extrapolation error but may not be from the best fossil record, raising the issue of nonindependence. An external calibration may be more independent but also may have a greater extrapolation error. Here, we used the advantages of both methods by applying a sequential calibration to avian molecular clocks. We estimated a basal divergence within birds, the split between fowl (Galliformes) and ducks (Anseriformes), to be 89.8 +/- 6.97 MYA using an external calibration and 12 rate-constant nuclear genes. In turn, this time estimate was used as an internal calibration for three species-rich avian molecular data sets: mtDNA, DNA-DNA hybridization, and transferrin immunological distances. The resulting time estimates indicate that many major clades of modern birds had their origins within the Cretaceous. This supports earlier studies that identified large gaps in the avian fossil record and suggests that modern birds may have coexisted with other avian lineages for an extended period during the Cretaceous. The new time estimates are concordant with a continental breakup model for the origin of ratites.     

Evolutionary timescale of monocots determined by the fossilized birth‐death model using a large number of fossil records

Although the phylogenetic relationships between monocot orders are sufficiently understood, a timescale of their evolution is needed. Several studies on molecular clock dating are available, but their results have been biased by their calibration schemes. Recently, the fossilized birth‐death model, a type of Bayesian dating method, was proposed, and it does not require prior calibration and allows the use all available fossils. Using this model, we conducted divergence‐time estimations of monocots to explore their evolutionary timeline without calibration bias. This is the first application of this model to seed plants. The dataset contained the matK and rbcL chloroplast genes of 118 monocot genera covering all extant orders. We employed information from 247 monocot fossils, which exceeded previous dating analyses that used a maximum of 12 monocot fossils. The crown group of monocots was dated to approximately the Late Jurassic–Early Cretaceous periods, and most extant monocot orders were estimated to diverge throughout the Early Cretaceous. Our results overlapped with the divergence time of insect lineages, such as beetles and flies, suggesting an association with pollinators in early monocot evolution. In addition, we proposed three new orders based on divergence time: Orchidales separated from Asparagales and Tofieldiales and Arales separated from Aslimatales.     

A molecular timescale for galliform birds accounting for uncertainty in time estimates and heterogeneity of rates of DNA substitutions across lineages and sites

A recent molecular timescale for major lineages of the Galliformes indicated that Megapodiidae and possibly Cracidae, originated in the Cretaceous, while the remaining families originated in the Tertiary. This timescale was based on clock-like evolution in genetic and taxonomic partitions of mitochondrial ND2 and cyt b DNA sequences, and assumed that ordinal diversification of Galloanserae around 90 million years ago and imposed, whenever appropriate, minimum age constraints based on the fossil record. This approach is not ideal, as it did not account for uncertainty in estimating branch lengths and time, including the calibration time, and heterogeneity in the rate of DNA substitution among sites and in different lineages. Furthermore all the information available in the DNA sequences was not included, and may have been affected by stochastic error in individual gene partitions. Here, we present a follow-up analysis by estimating divergence times using a Bayesian framework that accounts for these possible sources of uncertainty. Our results based on combined and separate analyses of mitochondrial DNA sequences comprised of 1756 sites of 12S rDNA, ND2 and cyt b indicated that (1) Megapodiidae and Cracidae, and likely Odontophoridae, originated in the Cretaceous; (2) estimates based on concatenated genes are less affected by stochastic error among sites and less influenced by the phylogenetic signals of individual gene partitions, which are unequally distributed along the phylogenetic tree; and (3) the use of only an external molecular calibration results in lower estimation of most ingroup node ages. We also point out that galliform fossils may not be as useful for point calibrations as was previously suggested, but instead may be better employed as priors for the estimation of node ages under a Bayesian approach.     

New specimens of the early Eocene stem group galliform Paraortygoides (Gallinuloididae), with comments on the evolution of a crop in the stem lineage of Galliformes

Two new specimens of the fossil stem group galliform Paraortygoides messelensis Mayr 2000 (Gallinuloididae) are described from the Middle Eocene of Messel in Germany, including a complete skeleton in which the hitherto unknown skull of this species is preserved. The shorter and more protruding crista deltopectoralis of the humerus, also for the first time visible in one of the new specimens, shows gallinuloidids to be the sister taxon of all other, extinct and extant, galliform birds. Gallinuloidids distinctly differ from modern Galliformes in several other plesiomorphic osteological features, mainly of the pectoral girdle, of which the absence of a spina interna on the sternum is here reported for the first time. It is assumed that major evolutionary transformations in the stem lineage of Galliformes are related to the evolution of a large crop, which appears to have been absent in gallinuloidids. The vegetarian food component of gallinuloidids thus probably mainly consisted of soft plant matter rather than coarse material such as seeds.     

Early penguin fossils, plus mitochondrial genomes, calibrate avian evolution

Testing models of macroevolution, and especially the sufficiency of microevolutionary processes, requires good collaboration between molecular biologists and paleontologists. We report such a test for events around the Late Cretaceous by describing the earliest penguin fossils, analyzing complete mitochondrial genomes from an albatross, a petrel, and a loon, and describe the gradual decline of pterosaurs at the same time modern birds radiate. The penguin fossils comprise four naturally associated skeletons from the New Zealand Waipara Greensand, a Paleocene (early Tertiary) formation just above a well-known Cretaceous/Tertiary boundary site. The fossils, in a new genus (Waimanu), provide a lower estimate of 61-62 Ma for the divergence between penguins and other birds and thus establish a reliable calibration point for avian evolution. Combining fossil calibration points, DNA sequences, maximum likelihood, and Bayesian analysis, the penguin calibrations imply a radiation of modern (crown group) birds in the Late Cretaceous. This includes a conservative estimate that modern sea and shorebird lineages diverged at least by the Late Cretaceous about 74 +/- 3 Ma (Campanian). It is clear that modern birds from at least the latest Cretaceous lived at the same time as archaic birds including Hesperornis, Ichthyornis, and the diverse Enantiornithiformes. Pterosaurs, which also coexisted with early crown birds, show notable changes through the Late Cretaceous. There was a decrease in taxonomic diversity, and small- to medium-sized species disappeared well before the end of the Cretaceous. A simple reading of the fossil record might suggest competitive interactions with birds, but much more needs to be understood about pterosaur life histories. Additional fossils and molecular data are still required to help understand the role of biotic interactions in the evolution of Late Cretaceous birds and thus to test that the mechanisms of microevolution are sufficient to explain macroevolution.     

A crown‐group demosponge from the early Cambrian Sirius Passet Biota,North Greenland

Calibration of the divergence times of sponge lineages and understanding of their phylogenetic history are hampered by the difficulty in recognizing crown versus stem groups in the fossil record. A new specimen from the lower Cambrian (Series 2, Stage 3; approximately 515 Ma) Sirius Passet Biota of North Greenland has yielded a diagnostic spicule assemblage of the extant demosponge lineages Haploscleromorpha and/or Heteroscleromorpha. The specimen has disarticulated approximately in situ, but represents an individual sponge that possessed monaxon spicules combined with a range of slightly smaller sigma, toxa and unique spiral morphologies. The combination of spicule forms, together with their relatively large size, suggests that the sponge represents the stem lineage of Haploscleromorpha + Heteroscleromorpha. This is the first crown‐group demosponge described from the early Cambrian and provides the most reliable calibration point currently available for phylogenetic studies.     

The Paleogene fossil record of birds in Europe

The Paleogene (Paleocene-Oligocene) fossil record of birds in Europe is reviewed and recent and fossil taxa are placed into a phylogenetic framework, based on published cladistic analyses. The pre-Oligocene European avifauna is characterized by the complete absence of passeriform birds, which today are the most diverse and abundant avian taxon. Representatives of small non-passeriform perching birds thus probably had similar ecological niches before the Oligocene to those filled by modern passerines. The occurrence of passerines towards the Lower Oligocene appears to have had a major impact on these birds, and the surviving crown-group members of many small arboreal Eocene taxa show highly specialized feeding strategies not found or rare in passeriform birds. It is detailed that no crown-group members of modern 'families' are known from pre-Oligocene deposits of Europe, or anywhere else. The phylogenetic position of Paleogene birds thus indicates that diversification of the crown-groups of modern avian 'families' did not take place before the Oligocene, irrespective of their relative position within Neornithes (crown-group birds). The Paleogene fossil record of birds does not even support crown-group diversification of Galliformes, one of the most basal taxa of neognathous birds, before the Oligocene, and recent molecular studies that dated diversification of galliform crown-group taxa into the Middle Cretaceous are shown to be based on an incorrect interpretation of the fossil taxa used for molecular clock calibrations. Several taxa that occur in the Paleogene of Europe have a very different distribution than their closest extant relatives. The modern survivors of these Paleogene lineages are not evenly distributed over the continents, and especially the great number of taxa that are today restricted to South and Central America is noteworthy. The occurrence of stem-lineage representatives of many taxa that today have a restricted Southern Hemisphere distribution conflicts with recent hypotheses on a Cretaceous vicariant origin of these taxa, which were deduced from the geographical distribution of the basal crown-group members.     

Error in Estimation of Rate and Time Inferred from the Early Amniote Fossil Record and Avian Molecular Clocks

The best reconstructions of the history of life will use both molecular time estimates and fossil data. Errors in molecular rate estimation typically are unaccounted for and no attempts have been made to quantify this uncertainty comprehensively. Here, focus is primarily on fossil calibration error because this error is least well understood and nearly universally disregarded. Our quantification of errors in the synapsid–diapsid calibration illustrates that although some error can derive from geological dating of sedimentary rocks, the absence of good stem fossils makes phylogenetic error the most critical. We therefore propose the use of calibration ages that are based on the first undisputed synapsid and diapsid. This approach yields minimum age estimates and standard errors of 306.1±8.5 MYR for the divergence leading to birds and mammals. Because this upper bound overlaps with the recent use of 310 MYR, we do not support the notion that several metazoan divergence times are significantly overestimated because of serious miscalibration (sensu Lee 1999). However, the propagation of relevant errors reduces the statistical significance of the pre-K–T boundary diversification of many bird lineages despite retaining similar point time estimates. Our results demand renewed investigation into suitable loci and fossil calibrations for constructing evolutionary timescales.[Reviewing Editor: Martin Kreitman]     

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.     

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

Fossil calibration of molecular divergence infers a moderate mutation rate and recent radiations for pinus

Silent mutation rate estimates for Pinus vary 50-fold, ranging from angiosperm-like to among the slowest reported for plants. These differences either reflect extraordinary genomic processes or inconsistent fossil calibration, and they have important consequences for population and biogeographical inferences. Here we estimate mutation rates from 4 Pinus species that represent the major lineages using 11 nuclear and 4 chloroplast loci. Calibration was tested at the divergence of Pinus subgenera with the oldest leaf fossil from subg. Strobus (Eocene; 45 MYA) or a recently published subg. Strobus wood fossil (Cretaceous; 85 MYA). These calibrations place the origin of Pinus 190-102 MYA and give absolute silent rate estimates of 0.70-1.31x10(-9) and 0.22-0.42x10(-9).site-1.year-1 for the nuclear and chloroplast genomes, respectively. These rates are approximately 4- to 20-fold slower than angiosperms, but unlike many previous estimates, they are more consistent with the high per-generation deleterious mutation rates observed in pines. Chronograms from nuclear and chloroplast genomes show that the divergence of subgenera accounts for about half of the time since Pinus diverged from Picea, with subsequent radiations occurring more recently. By extending the sampling to encompass the phylogenetic diversity of Pinus, we predict that most extant subsections diverged during the Miocene. Moreover, subsect. Australes, Ponderosae, and Contortae, containing over 50 extant species, radiated within a 5 Myr time span starting as recently as 18 MYA. An Eocene divergence of pine subgenera (using leaf fossils) does not conflict with fossil-based estimates of the Pinus-Picea split, but a Cretaceous divergence using wood fossils accommodates Oligocene fossils that may represent modern subsections. Because homoplasy and polarity of character states have not been tested for fossil pine assignments, the choice of fossil and calibration node represents a significant source of uncertainty. Based on several lines of evidence (including agreement with ages inferred using calibrations outside of Pinus), we conclude that the 85 MYA calibration at the divergence of pine subgenera provides a reasonable lower bound and that further refinements in age and mutation rate estimates will require a synthetic examination of pine fossil history.     

Strong mitochondrial DNA support for a Cretaceous origin of modern avian lineages

Background  

Determining an absolute timescale for avian evolutionary history has proven contentious. The two sources of information available, paleontological data and inference from extant molecular genetic sequences (colloquially, 'rocks' and 'clocks'), have appeared irreconcilable; the fossil record supports a Cenozoic origin for most modern lineages, whereas molecular genetic estimates suggest that these same lineages originated deep within the Cretaceous and survived the K-Pg (Cretaceous-Paleogene; formerly Cretaceous-Tertiary or K-T) mass-extinction event. These two sources of data therefore appear to support fundamentally different models of avian evolution. The paradox has been speculated to reflect deficiencies in the fossil record, unrecognized biases in the treatment of genetic data or both. Here we attempt to explore uncertainty and limit bias entering into molecular divergence time estimates through: (i) improved taxon (n = 135) and character (n = 4594 bp mtDNA) sampling; (ii) inclusion of multiple cladistically tested internal fossil calibration points (n = 18); (iii) correction for lineage-specific rate heterogeneity using a variety of methods (n = 5); (iv) accommodation of uncertainty in tree topology; and (v) testing for possible effects of episodic evolution.     

Neither phylogenomic nor palaeontological data support a Palaeogene origin of placental mammals

O''Leary et al. (O''Leary et al. 2013 Science 339, 662–667. (doi:10.1126/science.1229237)) performed a fossil-only dating analysis of mammals, concluding that the ancestor of placentals post-dated the Cretaceous–Palaeogene boundary, contradicting previous palaeontological and molecular studies that placed the ancestor in the Cretaceous. They incorrectly used fossil ages as species divergence times for crown groups, while in fact the former should merely form minimum-age bounds for the latter. Statistical analyses of the fossil record have shown that crown groups are significantly older than the oldest ingroup fossil, so that fossils do not directly reflect the true ages of clades. Here, we analyse a 20 million nucleotide genome-scale alignment in conjunction with a probabilistic interpretation of the fossil ages from O''Leary et al. Our combined analysis of fossils and molecules demonstrates that Placentalia originated in the Cretaceous.     

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.     

Exploring the impact of fossil constraints on the divergence time estimates of derived liverworts

In this study, we evaluate the impact of fossil assignments and different models of calibration on divergence time estimates carried out as Bayesian analyses. Estimated ages from preceding studies and liverwort inclusions from Baltic amber are used as constraints on a molecular phylogeny of Cephaloziineae (Jungermanniopsida) obtained from sequences of two chloroplast coding regions: rbcL and psbA. In total, the comparison of 12 different analyses demonstrates that an increased reliability of the chronograms is linked to the number of fossils assigned and to the accuracy of their assignments. Inclusion of fossil constraints leads to older ages of most crown groups, but has no influence on lineage through time plots suggesting a nearly constant accumulation of diversity since the origin of Cephaloziineae in the early to Middle Jurassic. Our results provide a note of caution regarding the interpretation of chronograms derived from DNA sequence variation of extant species based on a single calibration point and/or low accuracy of the assignment of fossils to nodes in the phylogeny.     

Multiple lineages of lice pass through the K-Pg boundary

For modern lineages of birds and mammals, few fossils have been found that predate the Cretaceous-Palaeogene (K-Pg) boundary. However, molecular studies using fossil calibrations have shown that many of these lineages existed at that time. Both birds and mammals are parasitized by obligate ectoparasitic lice (Insecta: Phthiraptera), which have shared a long coevolutionary history with their hosts. Evaluating whether many lineages of lice passed through the K-Pg boundary would provide insight into the radiation of their hosts. Using molecular dating techniques, we demonstrate that the major louse suborders began to radiate before the K-Pg boundary. These data lend support to a Cretaceous diversification of many modern bird and mammal lineages.