ACCOUNTING FOR RATE VARIATION AMONG LINEAGES IN COMPARATIVE DEMOGRAPHIC ANALYSES

Genetic analyses of contemporary populations can be used to estimate the demographic histories of species within an ecological community. Comparison of these demographic histories can shed light on community responses to past climatic events. However, species experience different rates of molecular evolution, and this presents a major obstacle to comparative demographic analyses. We address this problem by using a Bayesian relaxed‐clock method to estimate the relative evolutionary rates of 22 small mammal taxa distributed across northwestern North America. We found that estimates of the relative molecular substitution rate for each taxon were consistent across the range of sampling schemes that we compared. Using three different reference rates, we rescaled the relative rates so that they could be used to estimate absolute evolutionary timescales. Accounting for rate variation among taxa led to temporal shifts in our skyline‐plot estimates of demographic history, highlighting both uniform and idiosyncratic evolutionary responses to directional climate trends for distinct ecological subsets of the small mammal community. Our approach can be used in evolutionary analyses of populations from multiple species, including comparative demographic studies.     

Do complex population histories drive higher estimates of substitution rate in phylogenetic reconstructions?

Our curiosity about biodiversity compels us to reconstruct the evolutionary past of species. Molecular evolutionary theory now allows parameterization of mathematically sophisticated and detailed models of DNA evolution, which have resulted in a wealth of phylogenetic histories. But reconstructing how species and population histories have played out is critically dependent on the assumptions we make, such as the clock-like accumulation of genetic differences over time and the rate of accumulation of such differences. An important stumbling block in the reconstruction of evolutionary history has been the discordance in estimates of substitution rate between phylogenetic and pedigree-based studies. Ancient genetic data recovered directly from the past are intermediate in time scale between phylogenetics-based and pedigree-based calibrations of substitution rate. Recent analyses of such ancient genetic data suggest that substitution rates are closer to the higher, pedigree-based estimates. In this issue, Navascués & Emerson (2009) model genetic data from contemporary and ancient populations that deviate from a simple demographic history (including changes in population size and structure) using serial coalescent simulations. Furthermore, they show that when these data are used for calibration, we are likely to arrive at upwardly biased estimates of mutation rate.     

Calibrating molecular estimates of substitution rates and divergence times in birds

Many evolutionary studies of birds rely on the estimation of molecular divergence times and substitution rates. In order to perform such analyses, it is necessary to incorporate some form of calibration information: a known substitution rate, radiometric ages of heterochronous sequences, or inferred ages of lineage splitting events. All three of these techniques have been employed in avian molecular studies, but their usage has not been entirely satisfactory. For example, the 'traditional' avian mitochondrial substitution rate of 2% per million years is frequently adopted without acknowledgement of the associated uncertainty. Similarly, fossil and biogeographic information is almost always converted into an errorless calibration point. In both cases, the resulting estimates of divergence times and substitution rates will be artificially precise, which has a considerable impact on hypothesis testing. In addition, using such a simplistic approach to calibration discards much of the information offered by the fossil record. A number of more sophisticated calibration methods have recently been introduced, culminating in the development of probability distribution-based calibrations. In this article, I discuss the use of this new class of methods and offer guidelines for choosing a calibration technique.     

Analyses of evolutionary dynamics in viruses are hindered by a time-dependent bias in rate estimates

Time-scales of viral evolution and emergence have been studied widely, but are often poorly understood. Molecular analyses of viral evolutionary time-scales generally rely on estimates of rates of nucleotide substitution, which vary by several orders of magnitude depending on the timeframe of measurement. We analysed data from all major groups of viruses and found a strong negative relationship between estimates of nucleotide substitution rate and evolutionary timescale. Strikingly, this relationship was upheld both within and among diverse groups of viruses. A detailed case study of primate lentiviruses revealed that the combined effects of sequence saturation and purifying selection can explain this time-dependent pattern of rate variation. Therefore, our analyses show that studies of evolutionary time-scales in viruses require a reconsideration of substitution rates as a dynamic, rather than as a static, feature of molecular evolution. Improved modelling of viral evolutionary rates has the potential to change our understanding of virus origins.     

Pleistocene evolutionary history of the Clouded Apollo (Parnassius mnemosyne): genetic signatures of climate cycles and a 'time-dependent' mitochondrial substitution rate

Genetic data are currently providing a large amount of new information on past distribution of species and are contributing to a new vision of Pleistocene ice ages. Nonetheless, an increasing number of studies on the 'time dependency' of mutation rates suggest that date assessments for evolutionary events of the Pleistocene might be overestimated. We analysed mitochondrial (mt) DNA (COI) sequence variation in 225 Parnassius mnemosyne individuals sampled across central and eastern Europe in order to assess (i) the existence of genetic signatures of Pleistocene climate shifts; and (ii) the timescale of demographic and evolutionary events. Our analyses reveal a phylogeographical pattern markedly influenced by the Pleistocene/Holocene climate shifts. Eastern Alpine and Balkan populations display comparatively high mtDNA diversity, suggesting multiple glacial refugia. On the other hand, three widely distributed and spatially segregated lineages occupy most of northern and eastern Europe, indicating postglacial recolonization from different refugial areas. We show that a conventional 'phylogenetic' substitution rate cannot account for the present distribution of genetic variation in this species, and we combine phylogeographical pattern and palaeoecological information in order to determine a suitable intraspecific rate through a Bayesian coalescent approach. We argue that our calibrated 'time-dependent' rate (0.096 substitutions/ million years), offers the most convincing time frame for the evolutionary events inferred from sequence data. When scaled by the new rate, estimates of divergence between Balkan and Alpine lineages point to c. 19 000 years before present (last glacial maximum), and parameters of demographic expansion for northern lineages are consistent with postglacial warming (5-11 000 years before present).     

A Bayesian evaluation of human mitochondrial substitution rates

Accurate estimates of mitochondrial substitution rates are central to molecular studies of human evolution, but meaningful comparisons of published studies are problematic because of the wide range of methodologies and data sets employed. These differences are nowhere more pronounced than among rates estimated from phylogenies, genealogies, and pedigrees. By using a data set comprising mitochondrial genomes from 177 humans, we estimate substitution rates for various data partitions by using Bayesian phylogenetic analysis with a relaxed molecular clock. We compare the effect of multiple internal calibrations with the customary human-chimpanzee split. The analyses reveal wide variation among estimated substitution rates and divergence times made with different partitions and calibrations, with evidence of substitutional saturation, natural selection, and significant rate heterogeneity among lineages and among sites. Collectively, the results support dates for migration out of Africa and the common mitochondrial ancestor of humans that are considerably more recent than most previous estimates. Our results also demonstrate that human mitochondrial genomes exhibit a number of molecular evolutionary complexities that necessitate the use of sophisticated analytical models for genetic analyses.     

Elevated substitution rates estimated from ancient DNA sequences

Ancient DNA sequences are able to offer valuable insights into molecular evolutionary processes, which are not directly accessible via modern DNA. They are particularly suitable for the estimation of substitution rates because their ages provide calibrating information in phylogenetic analyses, circumventing the difficult task of choosing independent calibration points. The substitution rates obtained from such datasets have typically been high, falling between the rates estimated from pedigrees and species phylogenies. Many of these estimates have been made using a Bayesian phylogenetic method that explicitly accommodates heterochronous data. Stimulated by recent criticism of this method, we present a comprehensive simulation study that validates its performance. For datasets of moderate size, it produces accurate estimates of rates, while appearing robust to assumptions about demographic history. We then analyse a large collection of 749 ancient and 727 modern DNA sequences from 19 species of animals, plants and bacteria. Our new estimates confirm that the substitution rates estimated from ancient DNA sequences are elevated above long-term phylogenetic levels.     

Bayesian relaxed clock estimation of divergence times in foraminifera

Accurate and precise estimation of divergence times during the Neo-Proterozoic is necessary to understand the speciation dynamic of early Eukaryotes. However such deep divergences are difficult to date, as the molecular clock is seriously violated. Recent improvements in Bayesian molecular dating techniques allow the relaxation of the molecular clock hypothesis as well as incorporation of multiple and flexible fossil calibrations. Divergence times can then be estimated even when the evolutionary rate varies among lineages and even when the fossil calibrations involve substantial uncertainties. In this paper, we used a Bayesian method to estimate divergence times in Foraminifera, a group of unicellular eukaryotes, known for their excellent fossil record but also for the high evolutionary rates of their genomes. Based on multigene data we reconstructed the phylogeny of Foraminifera and dated their origin and the major radiation events. Our estimates suggest that Foraminifera emerged during the Cryogenian (650-920 Ma, Neo-Proterozoic), with a mean time around 770 Ma, about 220 Myr before the first appearance of reliable foraminiferal fossils in sediments (545 Ma). Most dates are in agreement with the fossil record, but in general our results suggest earlier origins of foraminiferal orders. We found that the posterior time estimates were robust to specifications of the prior. Our results highlight inter-species variations of evolutionary rates in Foraminifera. Their effect was partially overcome by using the partitioned Bayesian analysis to accommodate rate heterogeneity among data partitions and using the relaxed molecular clock to account for changing evolutionary rates. However, more coding genes appear necessary to obtain more precise estimates of divergence times and to resolve the conflicts between fossil and molecular date estimates.     

Expansion dating: calibrating molecular clocks in marine species from expansions onto the Sunda Shelf Following the Last Glacial Maximum

The rate of change in DNA is an important parameter for understanding molecular evolution and hence for inferences drawn from studies of phylogeography and phylogenetics. Most rate calibrations for mitochondrial coding regions in marine species have been made from divergence dating for fossils and vicariant events older than 1-2 My and are typically 0.5-2% per lineage per million years. Recently, calibrations made with ancient DNA (aDNA) from younger dates have yielded faster rates, suggesting that estimates of the molecular rate of change depend on the time of calibration, decaying from the instantaneous mutation rate to the phylogenetic substitution rate. aDNA methods for recent calibrations are not available for most marine taxa so instead we use radiometric dates for sea-level rise onto the Sunda Shelf following the Last Glacial Maximum (starting ～18,000 years ago), which led to massive population expansions for marine species. Instead of divergence dating, we use a two-epoch coalescent model of logistic population growth preceded by a constant population size to infer a time in mutational units for the beginning of these expansion events. This model compares favorably to simpler coalescent models of constant population size, and exponential or logistic growth, and is far more precise than estimates from the mismatch distribution. Mean rates estimated with this method for mitochondrial coding genes in three invertebrate species are elevated in comparison to older calibration points (2.3-6.6% per lineage per million years), lending additional support to the hypothesis of calibration time dependency for molecular rates.     

Evaluating fossil calibrations for dating phylogenies in light of rates of molecular evolution: a comparison of three approaches

Evolutionary and biogeographic studies increasingly rely on calibrated molecular clocks to date key events. Although there has been significant recent progress in development of the techniques used for molecular dating, many issues remain. In particular, controversies abound over the appropriate use and placement of fossils for calibrating molecular clocks. Several methods have been proposed for evaluating candidate fossils; however, few studies have compared the results obtained by different approaches. Moreover, no previous study has incorporated the effects of nucleotide saturation from different data types in the evaluation of candidate fossils. In order to address these issues, we compared three approaches for evaluating fossil calibrations: the single-fossil cross-validation method of Near, Meylan, and Shaffer (2005. Assessing concordance of fossil calibration points in molecular clock studies: an example using turtles. Am. Nat. 165:137-146), the empirical fossil coverage method of Marshall (2008. A simple method for bracketing absolute divergence times on molecular phylogenies using multiple fossil calibration points. Am. Nat. 171:726-742), and the Bayesian multicalibration method of Sanders and Lee (2007. Evaluating molecular clock calibrations using Bayesian analyses with soft and hard bounds. Biol. Lett. 3:275-279) and explicitly incorporate the effects of data type (nuclear vs. mitochondrial DNA) for identifying the most reliable or congruent fossil calibrations. We used advanced (Caenophidian) snakes as a case study; however, our results are applicable to any taxonomic group with multiple candidate fossils, provided appropriate taxon sampling and sufficient molecular sequence data are available. We found that data type strongly influenced which fossil calibrations were identified as outliers, regardless of which method was used. Despite the use of complex partitioned models of sequence evolution and multiple calibrations throughout the tree, saturation severely compressed basal branch lengths obtained from mitochondrial DNA compared with nuclear DNA. The effects of mitochondrial saturation were not ameliorated by analyzing a combined nuclear and mitochondrial data set. Although removing the third codon positions from the mitochondrial coding regions did not ameliorate saturation effects in the single-fossil cross-validations, it did in the Bayesian multicalibration analyses. Saturation significantly influenced the fossils that were selected as most reliable for all three methods evaluated. Our findings highlight the need to critically evaluate the fossils selected by data with different rates of nucleotide substitution and how data with different evolutionary rates affect the results of each method for evaluating fossils. Our empirical evaluation demonstrates that the advantages of using multiple independent fossil calibrations significantly outweigh any disadvantages.     

Speciation in North American black basses,Micropterus (Actinopterygii: Centrarchidae)

Abstract.— The Pleistocene Epoch has been frequently cited as a period of intense speciation for a significant portion of temperate continental biotas. To critically assess the role of Pleistocene glaciations on the evolution of the freshwater fish clade Micropterus , we use a phylogenetic analysis of complete gene sequences from two mitochondrial genes (cytochrome b and ND2), and a fossil calibration of the molecular clock to estimate ages of speciation events and rates of diversification. The absence of substantial morphological and ecological divergence together with endemism of five of the eight species in North American tributaries of the Gulf of Mexico may be interpreted as the result of a recent Pleistocene origin for these species. Speciation dates in Micropterus range from 1.01 ± 0.32 to 11.17 ± 1.02 million years ago. Only one speciation event is dated to the Pleistocene, and rates of diversification are not significantly variable in Micropterus. The premise that the Pleistocene was an exceptional period of speciation in Micropterus is not supported. Instead, a Gulf Coast allopatric speciation model is proposed, and predicts periods of dynamic speciation driven by sea level fluctuations in the Late Miocene and Pliocene. The Pleistocene, however, was a period of significant intraspecific mitochondrial lineage diversification. The application of the Gulf Coast allopatric speciation model to the remaining aquatic fauna of the Gulf of Mexico coast in North America will rely on robust phylogenetic hypotheses and accurate age estimations of speciation events.     

Time dependency of molecular rate estimates and systematic overestimation of recent divergence times

Studies of molecular evolutionary rates have yielded a wide range of rate estimates for various genes and taxa. Recent studies based on population-level and pedigree data have produced remarkably high estimates of mutation rate, which strongly contrast with substitution rates inferred in phylogenetic (species-level) studies. Using Bayesian analysis with a relaxed-clock model, we estimated rates for three groups of mitochondrial data: avian protein-coding genes, primate protein-coding genes, and primate d-loop sequences. In all three cases, we found a measurable transition between the high, short-term (< 1-2 Myr) mutation rate and the low, long-term substitution rate. The relationship between the age of the calibration and the rate of change can be described by a vertically translated exponential decay curve, which may be used for correcting molecular date estimates. The phylogenetic substitution rates in mitochondria are approximately 0.5% per million years for avian protein-coding sequences and 1.5% per million years for primate protein-coding and d-loop sequences. Further analyses showed that purifying selection offers the most convincing explanation for the observed relationship between the estimated rate and the depth of the calibration. We rule out the possibility that it is a spurious result arising from sequence errors, and find it unlikely that the apparent decline in rates over time is caused by mutational saturation. Using a rate curve estimated from the d-loop data, several dates for last common ancestors were calculated: modern humans and Neandertals (354 ka; 222-705 ka), Neandertals (108 ka; 70-156 ka), and modern humans (76 ka; 47-110 ka). If the rate curve for a particular taxonomic group can be accurately estimated, it can be a useful tool for correcting divergence date estimates by taking the rate decay into account. Our results show that it is invalid to extrapolate molecular rates of change across different evolutionary timescales, which has important consequences for studies of populations, domestication, conservation genetics, and human evolution.     

Exploring patterns and extent of bias in estimating divergence time from mitochondrial DNA sequence data in a particular lineage: a case study of salamanders (order Caudata)

In the practice of molecular dating, substitution saturation will bias the results if not properly modeled. Date estimates based on commonly used mitochondrial DNA sequences likely suffer from this problem because of their high substitution rate. Nevertheless, the patterns and extent of such expected bias remain unknown for many major evolutionary lineages, which often differ in ages, available calibrations, and substitution rates of their mitochondrial genome. In this case study of salamanders, we used estimates based on multiple nuclear exons to assess the effects of saturation on dating divergences using mitochondrial genome sequences on a timescale of ~200-300 My. The results indicated that, due to saturation for older divergences and in the absence of younger effective calibration points, dates derived from the mitochondrial data were considerably overestimated and systematically biased toward the calibration point for the ingroup root. The overestimate might be as great as 3-10 times (about 20 My) older than actual divergence dates for recent splitting events and 40 My older for events that are more ancient. For deep divergences, dates estimated were strongly compressed together. Furthermore, excluding the third codon positions of protein-coding genes or only using the RNA genes or second codon positions did not considerably improve the performance. In the order Caudata, slowly evolving markers such as nuclear exons are preferred for dating a phylogeny covering a relatively wide time span. Dates estimated from these markers can be used as secondary calibrations for dating recent events based on rapidly evolving markers for which mitochondrial DNA sequences are attractive candidates due to their short coalescent time. In other groups, similar evaluation should be performed to facilitate the choice of markers for molecular dating and making inferences from the results.     

Reconciling deep calibration and demographic history: bayesian inference of post glacial colonization patterns in Carcinus aestuarii (Nardo, 1847) and C. maenas (Linnaeus, 1758)

A precise inference of past demographic histories including dating of demographic events using bayesian methods can only be achieved with the use of appropriate molecular rates and evolutionary models. Using a set of 596 mitochondrial cytochrome c oxidase I (COI) sequences of two sister species of European green crabs of the genus Carcinus (C. maenas and C. aestuarii), our study shows how chronologies of past evolutionary events change significantly with the application of revised molecular rates that incorporate biogeographic events for calibration and appropriate demographic priors. A clear signal of demographic expansion was found for both species, dated between 10,000 and 20,000 years ago, which places the expansions events in a time frame following the Last Glacial Maximum (LGM). In the case of C. aestuarii, a population expansion was only inferred for the Adriatic-Ionian, suggestive of a colonization event following the flooding of the Adriatic Sea (18,000 years ago). For C. maenas, the demographic expansion inferred for the continental populations of West and North Europe might result from a northward recolonization from a southern refugium when the ice sheet retreated after the LGM. Collectively, our results highlight the importance of using adequate calibrations and demographic priors in order to avoid considerable overestimates of evolutionary time scales.     

A protocol for temporal calibration of general area cladograms

Aim To describe a protocol for incorporating a temporal dimension into historical biogeographical analysis, while maintaining the essential independence of all datasets, involving the generation of general area cladograms. Location Global. Methods General area cladograms (GACs) are a reconstruction of the evolutionary history of a set of areas and unrelated clades within those areas. Nodes on a GAC correspond to speciation events in a group of taxa; general nodes are those at which multiple unrelated clades speciate. We undertake temporal calibration of GACs using molecular clock estimates of splitting events between extant taxa as well as first appearance data from the fossil record. We present two examples based on re‐analysis of previously published data: first, a temporally calibrated GAC generated from secondary Brooks parsimony analysis (BPA) of six extant bird clades from the south‐west of North America using molecular clock estimates of divergence times; and second, an analysis of African Neogene mammals based on a phylogenetic analysis for comparing trees (PACT) analysis. Results A hypothetical example demonstrates how temporal calibration reveals potentially critical information about the timing of both unique and general events, while also illustrating instances of incongruence between dates generated from molecular clock estimates and fossils. For the African Neogene mammal dataset, our analysis reveals that most mammal clades underwent geodispersal associated with the Neogene climatic optimum (c. 16 Ma) and vicariant speciation in central Africa correlated with increased aridity and cooler temperatures around 2.5 Ma. Main conclusions Temporally calibrated GACs are valuable tools for assessing whether coordinated patterns of speciation are associated with large‐scale climatic or tectonic phenomena.     

Effects of Mitochondrial DNA Rate Variation on Reconstruction of Pleistocene Demographic History in a Social Avian Species,Pomatostomus superciliosus

Mitochondrial sequence data is often used to reconstruct the demographic history of Pleistocene populations in an effort to understand how species have responded to past climate change events. However, departures from neutral equilibrium conditions can confound evolutionary inference in species with structured populations or those that have experienced periods of population expansion or decline. Selection can affect patterns of mitochondrial DNA variation and variable mutation rates among mitochondrial genes can compromise inferences drawn from single markers. We investigated the contribution of these factors to patterns of mitochondrial variation and estimates of time to most recent common ancestor (T_MRCA) for two clades in a co-operatively breeding avian species, the white-browed babbler Pomatostomus superciliosus. Both the protein-coding ND3 gene and hypervariable domain I control region sequences showed departures from neutral expectations within the superciliosus clade, and a two-fold difference in T_MRCA estimates. Bayesian phylogenetic analysis provided evidence of departure from a strict clock model of molecular evolution in domain I, leading to an over-estimation of T_MRCA for the superciliosus clade at this marker. Our results suggest mitochondrial studies that attempt to reconstruct Pleistocene demographic histories should rigorously evaluate data for departures from neutral equilibrium expectations, including variation in evolutionary rates across multiple markers. Failure to do so can lead to serious errors in the estimation of evolutionary parameters and subsequent demographic inferences concerning the role of climate as a driver of evolutionary change. These effects may be especially pronounced in species with complex social structures occupying heterogeneous environments. We propose that environmentally driven differences in social structure may explain observed differences in evolutionary rate of domain I sequences, resulting from longer than expected retention times for matriarchal lineages in the superciliosus clade.     

Rapid miocene-pliocene dispersal and evolution of Mediterranean rajid fauna as inferred by mitochondrial gene variation

Rajidae (colloquially known as skates and rays) experienced multiple and parallel adaptive radiations allowing high species diversity and great differences of species composition between regional faunas. Nevertheless, they show considerable conservation of bio-ecological, morphological and reproductive traits. The evolutionary history and dispersal of North-east Atlantic and Mediterranean rajid fauna were investigated throughout the sequence analysis of the control region and 16S rDNA mitochondrial genes. Molecular estimates of divergence times indicated recent origin and rapid dispersal of the present species. Compared with the ancient origin of the family (Late Cretaceous), the present species diversity arose in a relatively narrow time-window (12 Myr) from Middle Miocene to Early Pleistocene, likely by speciation processes related to dramatic geological and climatic events in the Mediterranean. Nucleotide substitution rates and phylogenetic relationships indicated Mediterranean endemic skates derived from sister species with wider distribution during Late Pliocene-Pleistocene. Skate phylogeny and systematics obtained using mitochondrial gene variation were largely consistent with those based on morpho-anatomical data.     

Pleistocene climatic oscillations and the speciation history of an alpine endemic and a widespread arctic-alpine plant

Even in cases in which geographic isolation appears to have driven the speciation of regional endemics, range shifts during the Pleistocene climatic oscillations may also have influenced their evolutionary history. Elucidating speciation history can provide novel insights into evolutionary dynamics following climatic oscillations. We demonstrated a sister relationship between the Japanese alpine endemic Cardamine nipponica and the currently allopatric, widespread arctic-alpine Cardamine?bellidifolia (Brassicaceae) based on internal transcribed spacer (ITS) sequences and 10 other nuclear genes. Speciation history was inferred using demographic parameters under the isolation with migration model. The estimated demographic parameters showed that the population size of C. nipponica was similar to that of C. bellidifolia and that gene flow occurred exclusively from C. nipponica to C. bellidifolia after speciation. The inferred speciation history, which included gene flow, suggests that geographic barriers between the peripheral C. nipponica and the widespread C. bellidifolia were reduced during the Pleistocene. The asymmetric introgression implies that genetic isolation may have been involved in the speciation of C. nipponica. Our results suggest that even currently allopatric species may not have diverged solely under geographic isolation, and that their evolutionary history may have been influenced by Pleistocene range dynamics.     

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.