Skip to content
Main -> Dating -> Human mitochondrial molecular clock - Wikipedia
  • 29.01.2019
  • by Zushakar

Human mitochondrial molecular clock - Wikipedia

Mitochondrial DNA inheritance pattern

The human mitochondrial molecular clock is the rate at which mutations have been accumulating in the mitochondrial genome of hominids during the course of human evolution. The archeological record of human activity from early periods in human prehistory is relatively limited and its interpretation has been controversial. Because of the uncertainties from the archeological record, scientists have turned to molecular dating techniques in order to refine the timeline of human evolution. A major goal of scientists in the field is to develop an accurate hominid mitochondrial molecular clock which could then be used to confidently date events that occurred during the course of human evolution. Estimates of the mutation rate of human mitochondrial DNA mtDNA vary greatly depending on the available data and the method used for estimation. The two main methods of estimation, phylogeny based methods and pedigree based methods, have produced mutation rates that differ by almost an order of magnitude. Current research has been focused on resolving the high variability obtained from different rate estimates.

The number of synonymous sites,was taken as the sum over all 4-fold and 2-fold degenerate sites. We note that the low occurrence of transversions justifies the counting of 2-fold sites as synonymous sites. Only 1. The overall rate was derived from the ratio of all substitutions to synonymous substitutions. The basic rate was determined by the human-chimpanzee comparisons. This value was equal to the non-synonymous to synonymous ratio of 0. The bibliography of the Mitomap tree from www.

Mitomap tree. Phylogenetic tree of human mitochondrial DNA coding region positions sequences from Mitomap web page as of March, www. Inside the figure the legend explains the color code for the substitution classes.

Identified mis-assignments Table S2 have been underlined and the corrections maintaining the original color scheme have been added to the figure. Browse Subject Areas? Click through the PLOS taxonomy to find articles in your field. Abstract The molecular clock of mitochondrial DNA has been extensively used to date various genetic events. Introduction The idea of genes accumulating new mutations in a clock-like manner provides a versatile tool for dating genetic events [1][2] even when the clock is stochastic and its rate variable [3].

Results The Mitomap tree of human mitochondrial DNA coding region sequences Table S1Figure S1 was analyzed for time dependent changes in the relative frequencies of substitutions. Download: PPT. Figure 1. Table 1. The extent and significance of temporal changes in the mutation ratios. Figure 2. Temporal changes in the fractions of non-synonymous A and RNA plus intergenic mutations B after the removal of sub-trees defined by mutations favored by positive selection.

Table 2. Distribution of synonymous sites by the number of received mutation hits. Table 3. Distribution of synonymous sites according to 3 rates of substitution. Table 4. Comparison of the 1 rate and 3-rates models of synonymous transitions using the expected number of observed differences with the assumption of 6. Table 5. Proportion of observed transitional differences at synonymous sites between a human and a chimpanzee [56] sequences according to the sites' substitution rate.

Discussion The overall apparent substitution rate of mtDNA coding region shows no notable change for hundreds of thousands of years implying that the decrease of non-neutral mutations with time must be very slow Figure 3. Figure 3. Changes in the overall substitution rate of the coding region. Figure 4. Changes in the fraction of synonymous substitutions. Table 6. Dates of expansion of the main clades in human mtDNA tree by different calibrations.

Figure 5. Temporal changes in the relative rates of non-synonymous A and RNA plus intergenic mutations B without counting the mutations from the branches that carried only 1 or 2 individuals. Methods Revision of the Calibration of the Molecular Clock of Synonymous Positions The calibration of average synonymous transition rate was based on the accumulated transversions in rRNA and synonymous sites in between human and chimpanzee mtDNA lineages, assuming 6. The Mitomap Data Temporal changes in the proportion of non-synonymous and RNA gene variants in human mtDNA were analyzed in sequences considering variation of coding region between positions and The Sliding Window Analysis Temporal changes in the ratios of non-synonymous to synonymous and RNA to synonymous were assessed by sliding window analysis.

Generation of the Reference Dataset The empirical variation was compared to the variation of the references that otherwise resembled the empirical dataset but had uniform mutation ratios over all age groups. The Substitution Rate Over All Coding Sites The overall rate was derived from the ratio of all substitutions to synonymous substitutions.

Supporting Information. File S1.

Figure S1. Table S1. Mutation counts and coalescence ages of the analyzed clades. Table S2. Corrected misannotations of the mutations on the Mitomap tree. References 1. Zuckerkandl E, Pauling L Molecular disease, evolution, and genetic heterogeneity. In: Kasha M, Pullman B, editors. Horizons in Biochemistry. New York: Academic Press. Kumar S Molecular clocks: four decades of evolution. Nat Rev Genet 6: View Article Google Scholar 3.

Bromham L, Penny D The modern molecular clock. Nat Rev Genet 4: View Article Google Scholar 4. Cambridge, Massachusetts: Harvard University Press. Trends Genet View Article Google Scholar 6. Annu Rev Genet View Article Google Scholar 7. Am J Hum Genet View Article Google Scholar 8. View Article Google Scholar 9. Berlin Heidelberg: Springer-Verlag. Mol Biol Evol View Article Google Scholar Am J Phys Anthropol Heredity Subramanian S Temporal trails of natural selection in human mitogenomes.

Nat Genet Genetica 61 Moilanen JS, Majamaa K Phylogenetic network and physicochemical properties of nonsynonymous mutations in the protein-coding genes of human mitochondrial DNA. Science Genetics Penny D Evolutionary biology: relativity for molecular clocks. Nature PLoS Biol 6: e Trends Genet 79 Emerson BC Alarm bells for the molecular clock?

No support for Ho et al.

"Mitochondrial DNA variation in Human Origins and Disease"

Syst Biol Bandelt HJ Clock debate: when times are a-changin': time dependency of molecular rate estimates: tempest in a teacup.

Heredity 12.

Ingman M, Gyllensten U Rate variation between mitochondrial domains and adaptive evolution in humans. Hum Mol Genet Hum Mutat Trends Ecol Evol Genomics Curr Biol 18. Nat Genet J Hum Evol Cambridge Archeol J 7: 3 Mol Ecol PLoS Genet 4: e Int J Legal Med 64 Biochem Biophys Res Commun Achaz G Testing for neutrality in samples with sequencing errors.

Human mitochondrial molecular clock

Genome Res BMC Genomics 8: J Mol Evol Nat Genet 7: Cell In: Munro HN, editor. Mammalian protein metabolism.

The American Journal of Human Genetics BMC Genet 2: We estimated SEs based on the Bayesian framework described in ref. We thank Claude Bherer and Mark Lipson for helpful discussions about characterizing the uncertainty in the genetic maps.

We thank Thomas Higham for helpful discussions about the biases and reliability of radiocarbon dates. The authors declare no conflict of interest. This article contains supporting information online at www. National Center for Biotechnology InformationU. Published online May 2. Author information Copyright and License information Disclaimer.

Email: ude. Edited by Andrew G. Copyright notice. This article has been cited by other articles in PMC. Significance We report a method for dating ancient human samples that uses the recombination clock. Keywords: molecular clock, generation interval, ancient DNA, branch shortening. Abstract The study of human evolution has been revolutionized by inferences from ancient DNA analyses. Results Model and Simulations. Accounting for Uncertainty in Parameters in Real Data.

Case Studies. Open in a separate window. Robustness of Age Estimates. Historical Generation Interval in Humans. Discussion We have developed a genetic approach for dating ancient human specimens that is applicable for dating ancient non-African samples that share a history of Neanderthal admixture with extant non-Africans. Comparison with Radiocarbon Dating. Supplementary Material Supplementary File Click here to view.

Acknowledgments We thank Claude Bherer and Mark Lipson for helpful discussions about characterizing the uncertainty in the genetic maps. Footnotes The authors declare no conflict of interest. References 1. Pickrell JK, Reich D. Toward a new history and geography of human genes informed by ancient DNA.

Trends Genet. Bronk Ramsey C. Radiocarbon dating: Revolutions in understanding. Godwin H. Half-life of radiocarbon. Fu Q, et al. A revised timescale for human evolution based on ancient mitochondrial genomes. Curr Biol.

We recalibrate the molecular clock of human mtDNA as years The molecular clock of mitochondrial DNA has been extensively used to date various genetic events . Such an approach was taken in a recent study [16]. A generalized leastsquares method was applied in fitting a model to mtDNA sequence Although there is some uncertainty in the clock, this dating may pose a. Mitochondrial dna dating method - Find single man in the US with mutual relations. Looking for love in all the wrong places? Now, try the right place. Register.

Stadler T, Yang Z. Dating phylogenies with sequentially sampled tips. Syst Biol. Meyer M, et al. A high-coverage genome sequence from an archaic Denisovan individual. Determinants of mutation rate variation in the human germline. Annu Rev Genomics Hum Genet. Hinch AG, et al. The landscape of recombination in African Americans.

The date of interbreeding between Neandertals and modern humans. PLoS Genet. Green RE, et al. A draft sequence of the Neandertal genome. Selection and reduced population size cannot explain higher amounts of Neandertal ancestry in East Asian than in European human populations.

Am J Hum Genet. Vernot B, Akey JM. Resurrecting surviving Neandertal lineages from modern human genomes. Chakraborty R, Weiss KM. Admixture as a tool for finding linked genes and detecting that difference from allelic association between loci. Hellenthal G, et al. A genetic atlas of human admixture history. Moorjani P, et al. Loh P-R, et al. Inferring admixture histories of human populations using linkage disequilibrium. Genome sequence of a 45,year-old modern human from western Siberia.

The complete genome sequence of a Neanderthal from the Altai Mountains. Kong A, et al. Fine-scale recombination rate differences between sexes, populations and individuals. Common and low-frequency variants associated with genome-wide recombination rate. Nat Genet. Fenner JN. Cross-cultural estimation of the human generation interval for use in genetics-based population divergence studies. Am J Phys Anthropol. A populationwide coalescent analysis of Icelandic matrilineal and patrilineal genealogies: Evidence for a faster evolutionary rate of mtDNA lineages than Y chromosomes.

Amster G, Sella G.

Fossil evidence has been frequently used to estimate a date for the MRCA of two Attempts at calculating the human mitochondrial DNA (mtDNA) Using a previously published contamination estimation method [11], four. By this method, scientists have calculated that the common ancestor of all this places the date of the common ancestor at around 6, years. The HGDP-CEPH mtDNA sequences are part of a wider sequencing the PS and SS methods as recently recommended by Baele et al. Uncertainty around the dating was modeled in BEAST using a.

Life history effects on the molecular clock of autosomes and sex chromosomes. Sun JX, et al. A direct characterization of human mutation based on microsatellites. The Genomes Project Consortium A map of human genome variation from population-scale sequencing.

Nature Higham T, et al. The timing and spatiotemporal patterning of Neanderthal disappearance. Rasmussen M, et al. The genome of a Late Pleistocene human from a Clovis burial site in western Montana. Using tip calibration, we also estimated the coalescence dates of various nodes of interests in the tree.

We obtained a value of 4. This estimation may appear too young when compared with the dates that are generally derived from the fossil record. We estimated a split time between Homo neanderthalensis and H. Considering the widely held view that H. The Boxgrove tibia from Sussex, Englandattributed to H.

This estimation rather places a conservative upper bound of 93 ka for the time of the last major gene exchange between non-African and sub-Saharan African populations. As pointed out by Fu, Mittnik, et al. Finally, our results also allowed us to check whether the coalescence dates of some major haplogroups associated with human migrations table 2 are consistent with the archaeological evidence supplementary appendix S However, in the case of the Canary Islands, Remote Oceania, New Zealand, and the Americas, the estimated coalescence times were systematically older than the archaeological evidence.

Potential explanations for such discrepancies include ancestral polymorphism in the founding population or complex demographic histories involving multiples waves of colonists. N ote. In conclusion, our results demonstrate that the recent availability of ancient high-quality mtDNA genomes offers a powerful tool to robustly date past evolutionary events of our own species.

Using the age of ancient sequences leads to far more reproducible inferences and allows circumventing the large number of assumptions behind node and root calibration, which in turn should lead to an improvement of the estimation of human mitochondrial substitution rates. It should be possible to obtain increasingly narrow and precise substitution rate estimates by including additional ancient genomes in the analyses, as they will become available.

In this context, ancient isolates from geographic regions which are not represented yet, such as Africa and Australia, would be particularly helpful, as these would allow fine calibration of further clades in the human mitochondrial genome. From a more general point of view, the growing availability of ancient sequences due to sequencing technology improvements should allow reliable tip-calibrated phylogenetic rate and divergence time estimates to be obtained in many species for which internal nodes split times information are presently not available.

Our data set composed of cAMH and 30 ancient human complete mitochondrial genomes, consisting of both new and publicly available sequences. A total of new samples were obtained from two different sources. First, samples were selected by randomly choosing two individuals from each of the 51 populations of the HGDP-CEPH human genome diversity cell line panel Cann et al. Second, we randomly selected two individuals from each of the 21 Native American and one Siberian populations that had previously been genotyped at autosomal microsatellites Wang et al.

Details on the molecular and data processing are given in supplementary appendix S1Supplementary Material online. A total of public cAMH sequences were selected from GenBank to complete the geographic coverage and the haplogroup spectrum of the cAMH sequences generated in this project supplementary appendix S2Supplementary Material online. A chimpanzee sequence accession number HM was used as outgroup.

This resulted in a 16,bp aligned sequences matrix in which each nucleotide has been annotated by matching to the Cambridge reference annotation file using an in-house R script R Core Team It is well documented that substantial chemical modifications of nucleotide bases can be introduced postmortem as a result of DNA damage Paabo ; Sawyer et al. This feature has been used to distinguish between ancient damaged sequences and putative modern contaminants Green et al.

The new method refines the mtDNA calculation by taking into account the process of natural selection - which researchers realised was. We report a method for dating ancient human samples that uses the Keywords: molecular clock, generation interval, ancient DNA, A revised timescale for human evolution based on ancient mitochondrial genomes. Dating of the human-ape splitting by a molecular clock of mitochondrial DNA. Hasegawa M, Kishino H, Yano T. A new statistical method for estimating.

Here, we took advantage of the random nature of DNA damage to test for sequence quality in the ancient samples. Moreover, given that deamination errors are expected to happen at random, the probability of observing the same error in multiple sequences is low. The optimal partitioning scheme and the best-fit nucleotide substitution model for each partition of the mtDNA molecule was estimated using the software PartitionFinder Lanfear et al. We defined the following seven groups of nucleotides: 1 HVS1, 2 HVS2, 3, 4, 5 protein coding positions at 1st, 2nd, and 3rd codon, 6 tRNAs, 7 rRNAs and used PartitionFinder to analyze every scheme that includes those seven groups in any possible combination.

We also used the same algorithm to independently select the best model of nucleotide evolution to apply to the whole molecule. In all analyses on the partitioned data, substitution and clock models were unlinked, whereas tree topology was assumed to be the same between schemes.

We also performed separate phylogenetic inferences on the whole-mtDNA molecule. We compared constant size, logistic growth, and exponential growth models.

We specified a model where the probability that any transition sites of the alignment remained undamaged is assumed to decay exponentially with sample age. Model choice was based on the Bayes factors calculated from the marginal likelihoods, the latter being computed using both the PS and SS methods as recently recommended by Baele et al. Substitution models were based on the best fit from Partition Finder, and rate variation among sites both including and excluding invariant sites was modeled with a discrete gamma distribution with four rate categories.

An independent constant population size model was applied to the Neanderthal clade as recommended by Briggs et al. For each analysis, we ran four independent chains in which samples were drawn every 5, MCMC steps from a total of 50, steps, after a discarded burn-in of 5, steps. Parameter estimation was based on the samples combined from the different chains. The best supported tree was estimated from the combined samples using the maximum clade credibility method implemented in TreeAnnotator.

We compared the effect of different tree calibration strategies on the estimation of substitution rates and divergence times. Priors were placed on tips, internal nodes, and root age. Tip age was specified using reliable dates available for all high-quality ancient complete human mtDNA sequences available at the time.

Uncertainty around the dating was modeled in BEAST using a normal prior with standard deviation equal to the standard error of the radiocarbon date supplementary appendix S3Supplementary Material online. We first used date-randomization tests to determine whether the temporal and genetic information contained in the ancient sequences were sufficient for a thorough estimation of substitution rates.

From the total data set composed of sequencesages of the sequences were randomly shuffled ten times, and date-randomized data sets were reanalyzed with BEAST. Point-calibrated tips analyses i.

To prevent the chain from becoming stuck on unrealistic inflated values, flat priors uniform distributions were applied for all other internal nodes and root ages 0. Tip dating analyses were performed while considering aAMHs intraspecific calibration and archaic humans extra-specific calibration either separately or simultaneously.

Finally, to test whether the inferred substitution rates are dependent on a single-dated aAMH sequence, we performed BEAST analyses using each ancient mtDNA sequence independently as separate tip calibrations. Internal node calibration was conducted by specifying priors on the ages of some specific nodes in the tree. To avoid as much as possible any subjective bias, we enforced the following three criteria: 1 the anthropological evidence must correspond to a discrete expansion or migration, 2 at least one candidate haplotype has to be associated with the demographic event, and 3 the age of the event needs to be based on solid archaeological evidence.

We could identify ten different events which fitted all our criteria: The colonization of Sahul, Europe, Japan, America, Sardinia, Remote Oceania, Canary Islands, Madagascar, New Zealand, and the post glacial re-expansion of populations from the Franco-Iberian refuge areas see supplementary appendix S4Supplementary Material online, for details.

Mitochondrial dna dating method

Analyses were run considering all ten internal nodes either simultaneously or independently. These values were chosen for consistency with previous studies and information from the fossil record Benton and Donoghue Equal amounts of computational time number of MCMC iterations and independent chains as previously given were employed for the different calibration scenarios.

We also selected an equal number of posterior samples to estimate parameters and the marginal likelihood. Secondly, we investigated whether the calibration scenario had an influence on the tree topology. To do so, it is necessary to control for the MCMC error in the best supported topology that can be inferred among independent BEAST runs performed assuming the same calibration approach.

Topological distances between pairs of trees calculated on a sample of ten independent trees inferred per calibration scenario were generated using the Penny and Hendy statistic Penny and Hendywhich is defined as twice the number of internal branches defining different bipartitions of the tips.

To be comparable between all scenarios, trees were pruned to contain AMH tips only. We tested for deviations from a strict molecular clock by considering the standard deviation of the uncorrelated log-normal relaxed clock the parameter ucld. If this parameter is small close to 0there is little variation in rates among branches, and the data are clock like.

If this parameter takes values greater than 1 then the standard deviation in branch rates is greater than the mean rate, and the data exhibit substantial rate heterogeneity among lineages.

We also empirically investigated the effect of the chimpanzee and archaic human sequences on rate heterogeneity among lineages. To do so, we ran two analyses using alignments excluding those sequences one excluding the chimpanzee only and one excluding both the chimpanzee and Neanderthal sequences and performed model comparison through the evaluation of Bayes factors, as described previously.

The authors thank Mark Stoneking for generating the sequence data in his laboratory. They thank Mattias Jakobsson, Morten Rasmussen, Eske Willerslev, Maanasa Raghavan, and Michael Knapp for providing us access to sequence data and associated information in advance of their publication.

Oxford University Press is a department of the University of Oxford. It furthers the University's objective of excellence in research, scholarship, and education by publishing worldwide. Sign In or Create an Account. Sign In. Advanced Search. Article Navigation.

Close mobile search navigation Article Navigation. Volume Article Contents. Materials and Methods. Oxford Academic. Google Scholar. Anders Eriksson. Mingkun Li. Benjamin Sobkowiak. Lucy A. Vera Warmuth. Andres Ruiz-Linares. Andrea Manica. Associate editor: Beth Shapiro. The sequences reported in this article have been deposited in the GenBank under the accession nos.

The human mitochondrial molecular clock is the rate at which mutations have been . Pedigree methods estimate the mutation rate by comparing the mtDNA sequences of a . Methods/parameters for estimating date of mitochondrial Eve.

Cite Citation. Permissions Icon Permissions.


2 thoughts on “Human mitochondrial molecular clock - Wikipedia

Leave a Reply

Your email address will not be published. Required fields are marked *

Back to top