Summary
The more relevant retained-change family estimates give sapiens ages in the hundreds of thousands of years. Many faster estimates count variants whose long-term survival is unknown. Including archaic humans deepens the distance-based root; fossil dates introduce additional assumptions that must remain explicit. These calculations accompany the main article.
Reading the rate and age tables
The tables below compare mutations per maternal transmission in the region actually studied, then ask what age follows if that quantity is used as a constant lineage clock. They include raw detections, transmission-model estimates, near-homoplasmic or homoplasmic changes, and calibrated substitutions. Those are different measurements. A calculated age beside a raw detection rate is an illustration of the fast-clock argument, not an endorsement of that rate for dating ancestry.
Each figure sits under the stretch of the molecule its study sequenced, and spans exactly the columns it covers, so a control-region figure is never scaled onto the whole molecule. A cell spanning several stretches takes the blend of their colours. A number spanning HVR-I and HVR-II is their combined total, not that number in each region. Blank cells mean no separate measurement. Under each figure is the raw count as events / transmissions (or the model's formula), the denominator in bases, and for whole-genome cells the rate per site per year at 25-year generations. Figures marked model are converted or fitted, not counted. The printed base count is the study's denominator where known; it is not automatically the full width of a named region. HVR-I and HVR-II are short, variable parts of the larger control region. "Coding" here means everything outside the control region, including RNA genes; Arbeithuber's narrower definition excludes 88 intergenic bases. The header widths are only loosely proportional to length: drawn to scale, the coding region would take 93% of the table. Each table scrolls vertically and horizontally, with its header kept visible; study names link to their source notes.
Direct pedigree rows retain events divided by transmissions, and the Transmissions and Pedigrees columns give the study's sample. Ages use a common 25-year generation, so a pedigree age at 20 years is 0.8 times the displayed age, and at 30 years it is 1.2 times it. A study's own measured interval raises its age by that ratio: deCODE's 29.3 years by 1.172, Madrigal's 28.3 by 1.132, Connell's 26.9 by 1.076. The calibrated rows instead start with a published annual rate; their per-generation number is a 25-year display conversion and does not supply an independent clock. For a region of length L, annual rate r and generation interval g, μ = rLg. An age from a root-to-tip count B is Bg/μ = B/(rL).
Sapiens age uses Behar's published mean of 57.1 root-to-tip changes for whole mtDNA. Regional ages allocate that total using directly counted differences from the reconstructed RSRS root in a 54-record modern-human panel. That panel averages 50.241 whole-molecule differences, so its regional counts are multiplied by 57.1/50.241. For example, its control-region mean of 13.556 becomes a working root depth of 15.406. This is a profile approximation, not a regional measurement published by Behar, and it replaces an older allocation by share of pairwise differences (HVR-I 7.9, HVR-II 5.5, rest 1.8, control region 15.3, coding 41.8). Using the small panel's unscaled root counts would lower these sapiens ages by about 12%. It would not change the combined-human column.12
All-three age baseline includes sapiens, Neanderthal and Denisovan mtDNA. It takes the largest of the three between-group mean distances in the matched region and divides by twice the lineage rate, because a pairwise distance contains two branches. The full-molecule value is 360.830/2 = 180.415 changes, supplied by human-Denisovan comparisons. In the control region, Neanderthal-Denisovan comparisons supply the largest value, 52.623/2 = 26.312. The distance table marks each region's largest mean with a dagger. This is a distance-based approximation with all tips treated as contemporaneous. It is not a fitted three-group MRCA date, nor a statistical lower confidence bound. Fossil ages, repeated changes at the same site, branch-rate variation and topology need a fuller analysis. Half the average of all pooled pairs would underestimate the root because many pairs share younger ancestors.2
These are central-value comparisons. ka means thousand years and Ma means million years. A zero observed count supplies no finite age estimate. An asterisked "greater than" age instead uses the idealized one-sided 95% Poisson upper rate, −ln(0.05) divided by transmissions. That limit assumes complete detection of the stated class and independent events; it is not a bound valid under arbitrary selection, missed mutations or rate changes. Rows from the same study, and pooled rows containing earlier studies, are not independent evidence. The six rows marked "as tabulated by Howell 2003" are taken from Howell's Table 2, not from the original papers.
Region boundaries follow Connell 2022's Table 2 in revised Cambridge Reference Sequence numbering; the control region wraps through position 1, which is why HVR-I sits at the high end of the numbering and HVR-II at the low end. Exact windows were counted when recovered, and they differ from the header in places: Howell 1996's 1,194 bases (the control region plus 72 coding bases), Heyer's 673, Forster's 579, Sigurðardóttir's 673-base mask, deCODE's callable 1,092 and 15,443, and in Connell's partition HVR-II ending at position 371, which leaves 447 bases for the rest of the control region. Approximate ages remain visibly qualified where endpoints or site classes are unresolved. The common sequence alignment masks length-variable tracts and ignores ambiguous bases; a study may detect events in sites that this alignment excludes. Thus even an exact interval match is not a complete match of assay, masking and mutation class. Short control-region windows are particularly vulnerable to recurrent mutations. The calculation inputs, full precision results, regional counts, distance panels and code make those choices inspectable.
Measured directly, in living families
| Study | Whole genome1–16,569 · 16,569 bp | Count / method | Transmissions | Pedigrees | Sapiens age | All-three age baseline | Method and caveats | |||
|---|---|---|---|---|---|---|---|---|---|---|
| Control region16,024–576 · 1,122 bp | Coding region577–16,023 · 15,447 bp | |||||||||
| HVR-I16,024–16,383 · 360 bp | HVR-II57–372 · 316 bp | Rest16,384–56, 373–576 · 446 bp | ||||||||
| Howell 1996 | 0.02272 / 88 · 1,194 bp | — | DetectedNew germline mutations, both heteroplasmic (HVR-II nt 152 and 195), plus one somatic | 88 | 1 LHON pedigree, 45 members | 17 ka | 29.5 ka | Opened the pedigree-rate literature: the “ten times faster” claim began with two events in one disease pedigree. The primary sequenced nt 16024–16569 and 1–648 (about 1,194 bp), the control region plus 72 adjoining coding bases. Howell 2003 restates 2/88; Sigurðardóttir 2000 says 2/81. | ||
| Bendall 1996as tabulated by Howell 2003 | 0.01114 / 360 · 360 bp | — | — | — | DetectedNewly arising substitutions, heteroplasmic or homoplasmic, as Howell pooled them | 360 | Not given in Howell’s table | 17.7 ka | 32.8 ka | Howell 2003 Table 2: 4/360, divergence rate 0.99 (99.5% interval 0.14–3.35) per site per Myr, computed over Howell’s 1,122-bp convention at 20-year generations. HVR-1 only. The exact window was not checked; the age uses the 360-bp HVR-I window. |
| Mumm 1997as tabulated by Howell 2003 | 0.01691 / 59 · 360 bp | — | — | — | DetectedNewly arising substitutions, heteroplasmic or homoplasmic, as Howell pooled them | 59 | Not given in Howell’s table | 11.6 ka | 21.5 ka | Howell 2003 Table 2: 1/59, divergence rate 1.51 (99.5% interval 0.00–11.71) per site per Myr, computed over Howell’s 1,122-bp convention at 20-year generations. HVR-1 only; one event. The exact window was not checked; the age uses the 360-bp HVR-I window. |
| Parsons 1997 | 0.030610 / 327 · ~610 bp | — | — | DetectedSubstitutions between close relatives; heteroplasmic and homoplasmic both counted, per Howell 2003 | 327 | 134 maternal lineages of close relatives | 10.2 ka | 14.4 ka | The load-bearing young-earth citation. Ten differences across about 610 bp of HVR-I and HVR-II, not the whole molecule. The ~610 bp and the 20-year generation are reverse-engineered from the published 2.5/site/Myr; exact endpoints are unverified, so the age uses a broad HVR-I+II profile. The abstract gives the count. Howell’s Table 2 books this laboratory’s data as Parsons and Holland 1998 (10/306) plus the Forensic Science Service set (1/32): 11/338 together, rows below. | |
| Parsons 1997: Forensic Science Service setas tabulated by Howell 2003 | 0.03131 / 32 · 676 bp | — | — | DetectedNewly arising substitutions, heteroplasmic or homoplasmic, as Howell pooled them | 32 | Not given in Howell’s table | 11.1 ka | 15.6 ka | Howell 2003 Table 2: 1/32, divergence rate 2.76 (99.5% interval 0.01–20.43) per site per Myr, computed over Howell’s 1,122-bp convention at 20-year generations. HVR-1 and HVR-2. Howell’s footnote: “the analyses at the U.K. Forensic Science Service that were reported in Parsons et al. (1997).” One event. The age uses the 676-bp HVR-I+II window. | |
| Soodyall 1997as tabulated by Howell 2003 | 00 / 108 · 676 bp | — | — | DetectedNewly arising substitutions, heteroplasmic or homoplasmic, as Howell pooled them | 108 | Not given in Howell’s table | > 12.5 ka* | > 17.6 ka* | Howell 2003 Table 2: 0/108, divergence rate 0.00 (99.5% interval 0.00–4.79) per site per Myr, computed over Howell’s 1,122-bp convention at 20-year generations. HVR-1 and HVR-2. The Tristan da Cunha study, founded by a handful of women in the nineteenth century; Gibbons 1998 reports the same zero in 108. A zero count does not establish a zero rate. The age uses the 676-bp HVR-I+II window. | |
| Parsons and Holland 1998as tabulated by Howell 2003 | 0.032710 / 306 · 676 bp | — | — | DetectedNewly arising substitutions, heteroplasmic or homoplasmic, as Howell pooled them | 306 | Not given in Howell’s table | 10.6 ka | 14.9 ka | Howell 2003 Table 2: 10/306, divergence rate 2.91 (99.5% interval 0.99–6.44) per site per Myr, computed over Howell’s 1,122-bp convention at 20-year generations. HVR-1 and HVR-2; the same laboratory as Parsons 1997. The exact window was not checked; the age uses the 676-bp HVR-I+II window. | |
| Cavelier 2000: CR | 00 / 292 · 1,122 bp | — | HomoplasmicHomoplasmic substitutions only | 292 | 33 maternal lineages, 194 people | > 37.5 ka* | > 64.1 ka* | No homoplasmic substitutions detected; a zero count does not establish a zero mutation rate. Only length heteroplasmy, at three poly-C sites (65, 309, 16189). The zero result routinely left out of compilations on both sides. The limit assumes independent detectable events. Howell 2003 books the study as HVR-1 and HVR-2, 0/292. | ||
| Cavelier 2000: coding sample | — | — | — | 00 / 292 · positions 5,550–6,050 · 500 bp | HomoplasmicHomoplasmic substitutions only | Same 292 | Same | > 1.23 ka* | > 6.79 ka* | 5550–6050 is 501 inclusive positions; the paper calls it 500 bp. The zero-count limit uses that window. |
| Sigurðardóttir 2000 | 0.004263 / 705 · 673 bp | — | — | HomoplasmicBase substitutions between relatives: the three homoplasmic changes, per Howell 2003 | 705 | 26 large Icelandic pedigrees, 272 people | 81.1 ka | 114 ka | Largest pre-genomic study, about 8 times slower than Parsons. The published HVR mask is 673 bp; the fully sequenced amplicon is nt 16012–394 (~952 bp), and inclusive endpoint arithmetic does not reproduce 673. Madrigal’s Table 3 pairs Howell’s count with Sigurðardóttir’s rate. | |
| Sigurðardóttir 2000: rescored | 0.007095 / 705 · 673 bp | — | — | Detected; rescoredHomoplasmic and heteroplasmic changes, Howell’s rescoring | Same 705 | Same | 48.7 ka | 68.6 ka | Howell’s five-event rescoring of the same data, using the study’s published 673-bp HVR mask; not an independent study. | |
| Heyer 2001 | 0.007874 / 508 · 673 bp | — | — | Pedigree differencesSubstitutions between relatives; pooled with earlier studies, 18 / 1,729 = 0.0104 | 508 | Deep-rooting French-Canadian pedigrees (Quebec) | 44 ka | 62 ka | HVR-I 16024–16383 and HVR-II 58–370. Only 12 of the pooled 18 mutations sat at independent sites; three positions mutated twice. Fitted model: 40 fast sites (CI 27–54) at 274 per site per million generations. The pooled count overlaps other studies. | |
| Forster 2002: low radiation | 0.0051 / 200 · 579 bp | — | Detected heteroplasmyNew heteroplasmies only; none homoplasmic at either dose | 200 | 248 Kerala families, 988 people, across 2 to 4 generations (795 transmissions in both groups) | 68.2 ka | 86.5 ka | 16093–16390 and 35–315: parts of both HVRs plus adjoining bases. Kerala natural radioactivity, low-dose group. | ||
| Forster 2002: high radiation | 0.03722 / 595 · 579 bp | — | Detected heteroplasmyNew heteroplasmies only; none homoplasmic | 595 | Same families | 9.22 ka | 11.7 ka | A 7.4-fold rise in heteroplasmies appearing, and zero fixations. The radiation-hit positions “mutate, on average, 12 times faster during evolution.” The exposure comparison is not a universal clock. | ||
| Howell 2003: own pedigreesas tabulated by Howell 2003 | 0.005411 / 185 · 1,122 bp | — | DetectedNewly arising substitutions, heteroplasmic or homoplasmic, as Howell pooled them | 185 | Not given in Howell’s table | 71.3 ka | 122 ka | Howell 2003 Table 2: 1/185, divergence rate 0.48 (99.5% interval 0.00–3.87) per site per Myr, computed over Howell’s 1,122-bp convention at 20-year generations. Control region nt 16024–16569 and 1–576 in five LHON families. One event; included in the pooled row below. | ||
| Howell 2003: pooled CR | 0.010628 / 2,633 · pooled · 1,122 bp | — | Detected; pooledEverything newly arising, heteroplasmic or homoplasmic | 2,633 | 5 LHON pedigrees (185 transmissions) plus nine earlier data sets | 36.2 ka | 61.9 ka | The canonical pooled rate, published as a divergence rate of 0.95/site/Myr (twice the mutation rate; 99.5% interval 0.53–1.57) at 20-year generations, treating every study as 1,122 bp. Multiplying it across 16,569 sites is Carter’s error. Regions and assays differ between the pooled studies. Table 2’s ten study rows sum to 2,643 transmissions, not the printed 2,633. | ||
| Howell 2003: coding counts | — | — | — | 0.008664 / 462 · Subset | Pedigree eventsCoding-region events in partially sequenced LHON pedigrees | 462 | LHON pedigrees and Cavelier’s families | Not calculated | Not calculated | Mixed coverage: earlier LHON fragments, later coding sequences, Cavelier’s 500 bp. The 6,773 bp in Howell 1996 is not the common mask of all 462 transmissions. The raw count has no matched age. |
| Howell 2003: coding approximation | — | — | — | 0.00927 model0.06/site/Myr divergence, halved, × 15,447 × 20 y · 15,447 bp | Full-region lower estimateThe same coding events, converted to a rate | Same | Same | 112 ka | 427 ka | Published 0.06/site/Myr divergence rate (99.5% CI 0.008–0.20), halved, then converted using the study’s 20-year generation. Treating partial coverage as complete coding coverage makes this conservative; it is not a directly observed whole-coding rate. |
| Santos 2005: germinal | 0.01876 / 321 · ~610 bp | — | — | Detected; germinal classNew point substitutions, all heteroplasmic, classified somatic or germinal | 321 | Azores families; heteroplasmy in 72% of them | 16.7 ka | 23.6 ka | Six germinal events; 0.034 = 11/321 counting the five judged somatic. Counts from the abstract; exact windows are unverified, so the regional age is approximate. | |
| Santos 2005: corrected | 0.0037 modelmodel, corrected to fixation · ~610 bp | — | — | Transmission modelThe germinal events, corrected for loss before fixation | Same 321 | Same | 84.4 ka | 119 ka | First pedigree study to apply the arising-to-fixing correction; its corrected 0.2415/site/Myr lands “in the upper end” of phylogenetic values. A model-corrected rate for the measured region, not a direct homoplasmic count; it must not be assigned to 16,569 bases. | |
| Madrigal 2012: maximum | 0.02567 / 273 · 360 bp | — | — | — | Pedigree differencesSubstitutions between relatives; heteroplasmy status not stated | 273 | 19 deep-rooting Costa Rican pedigrees of at least 7 generations; 53 people typed | 7.67 ka | 14.2 ka | HVR-I, 360 bp. Generation time measured, 28.3 years over 289 meioses: at that interval multiply these ages by 1.132. The text reports 7 events over 299 meioses; its printed 71.1 × 10-6 per site per generation equals 7/273. Three of the mutated sites are known hot spots. The abstract’s 1.24 × 10-6 per site per year appears nowhere in the results. |
| Madrigal 2012: minimum | 0.009092 / 220 · 360 bp | — | — | — | Pedigree differencesSubstitutions between relatives; heteroplasmy status not stated | 220 | Same | 21.6 ka | 40.1 ka | Excluding a pedigree and hot-spot sites as the paper does leaves 2/220. Different branch and hot-spot exclusions define the two counts. |
| Guo 2013 | 00 / 26 · 3,981× depth · 16,569 bp · < 2.78 × 10-7 per site per year | HomoplasmicHomoplasmic changes; heteroplasmies followed down to 0.1% | 26 | 18 families, 44 people; mother–child pairs plus 9 sibling pairs | > 12.4 ka* | > 39.1 ka* | No discordant homoplasmic SNPs; 12 heteroplasmies shared by mother and child. A zero count does not exclude rare events, and 26 transmissions give only a weak bound. The mothers were childhood-cancer survivors with uterine irradiation, which biases the count upward. Implies a bottleneck near 200, the high outlier. | |||
| Rebolledo-Jaramillo 2014: arising | 0.143 modelα = μ0/(2 ln 99); μ0 = 51/(39 × 16,569) observed · 16,569 bp · 3.45 × 10-7 per site per year | Arising-rate modelGermline point heteroplasmies at ≥ 1% in both blood and buccal tissue, fed into a transmission model | 39 | 39 mother–child pairs, two tissues each | 9.98 ka | 31.5 ka | Fitted arising parameter α. The paper estimates μ0 as 51/(39 × 16,569 bp): the 51 are heteroplasmies observed above its 1% threshold, which include inherited variation; they are not 51 new mutations. 39% of heteroplasmies were tissue-specific, so somatic. | |||
| Rebolledo-Jaramillo 2014: retained | 0.0045 model2.7 × 10-7 per site per generation × 16,569 ≈ 0.0045; bottleneck 32.3 · 16,569 bp · 1.09 × 10-8 per site per year | Transmission modelThe same heteroplasmies, divided by the measured bottleneck | Same 39 | Same | 317 ka | 1 Ma | The paper that writes μ = α/N. Fitted μ after frequency and bottleneck modelling (IQR 8.3 × 10-8 to 8.2 × 10-7 per site per generation), not a direct count of homoplasmic changes. The authors call their estimate a lower bound. The best usable clock parameter from this study. | |||
| Rebolledo-Jaramillo 2014: D-loop retained | 0.00168 model1.5 × 10-6 per site per generation × 1,122 · 1,122 bp | — | Transmission modelD-loop heteroplasmies, modelled separately | Same 39 | Same | 229 ka | 391 ka | The paper’s separately modelled D-loop rate; its interquartile range, 4.8 × 10-7 to 4.7 × 10-6, spans tenfold. | ||
| Ding 2015: heteroplasmy gains | 0.423141 / 333 · ≥ 4% MAF · ~180× · 16,569 bp · 1.02 × 10-6 per site per year | Detected discordanceHeteroplasmic discordance between mother and child, lymphocytes only | 333 | 333 Sardinian parent–child trios | 3.37 ka | 10.7 ka | Child variants absent from the mother’s called class. Discordance between tissues sampled from mother and child is not automatically a new germline mutation; threshold, somatic origin and calling errors prevent identifying a retained rate. | |||
| Ding 2015: homoplasmy discordance | 0.18963 / 333 · 35 gains (0.105) + 28 losses · 16,569 bp · 4.57 × 10-7 per site per year | Discordant callsHomoplasmic discordance between mother and child, either direction | Same 333 | Same | Not calculated | Not calculated | The dataset behind Jeanson’s rate: 28 of his 63 discordances are variants the mother had and the child lost. The authors read the 35 child-not-mother calls as noise, and they were never re-examined at the read level. Gains and losses cannot both simply be counted as new inherited substitutions. Not a validated mutation rate. | |||
| Jeanson 2015 synthesis, for comparison | 0.15863 / 398 · 16,569 bp · 3.82 × 10-7 per site per year | Proposed mutation clockHomoplasmic discordances in Ding 2015, pooled with Guo and Rebolledo-Jaramillo (both zero), treated as homoplasmic changes | 398 | Three studies pooled | 9.02 ka | 28.5 ka | Mostly Ding’s unvalidated discordances. Predicts 4.1 homoplasmic changes in Guo’s 26 transmissions (0 seen), 54.6 in Norfolk’s 345 (2 seen) and 242 near-homoplasmic events in Wei’s 1,526 (3 seen). Shown to test the proposal, not as a measured rate. | |||
| Wei 2019: detected | 0.313477 / 1,526 · > 1% VAF · ~1,900× · 16,569 bp · 7.55 × 10-7 per site per year | Detected de novoNew heteroplasmies above 1% in the child’s blood | 1,526 | 1,526 mother–offspring pairs, whole blood | 4.57 ka | 14.4 ka | 614 maternal heteroplasmies lost against 477 gained. Offspring blood variants undetected in the mother; not all germline or destined to persist. Low-level against near-homoplasmic is a factor of 159 in one generation: a snapshot, not the eventual fate of the 477. | |||
| Wei 2019: control | 0.0858131 / 1,526 · 1,122 bp | — | Detected de novoThe same events in the D-loop | Same | Same | 4.49 ka | 7.66 ka | 131 D-loop events; the supplementary Table S2’s printed 13 is inconsistent with its rate and totals. Regional clock illustration only. | ||
| Wei 2019: coding | — | — | — | 0.227346 / 1,526 · 15,447 bp | Detected de novoThe same events outside the D-loop (genes, rRNA, tRNA) | Same | Same | 4.6 ka | 17.4 ka | 346 outside-D-loop events. Same tissue and threshold limitations. |
| Wei 2019: frequency weighted | 0.0194 model477 × 0.062 / 1,526 · 16,569 bp · 4.68 × 10-8 per site per year | Neutral illustrationThe 477, each weighted by its heteroplasmic fraction | Same | Same | 73.7 ka | 233 ka | Σ-frequency approximation. Blood is not the germline; selection and undetected variants remain unmodelled. | |||
| Wei 2019: above 80% | 0.001973 / 1,526 · 16,569 bp · 4.75 × 10-9 per site per year | Near-homoplasmicDe novo variants above 80% heteroplasmic fraction | Same | Same | 726 ka | 2.29 Ma | Two of the three events have ascertainment concerns. Above 80% is not complete homoplasmy, and no population-wide fixation estimate follows. | |||
| Zaidi 2019 | 0.00782 model4.72 × 10-7 per site per generation × 16,569 · ~11,400× · 16,569 bp · 1.89 × 10-8 per site per year | Transmission modelGermline heteroplasmies (in both blood and buccal) converted to a fixation rate through a bottleneck of 10.3 | 249 | 96 families, 345 people, 2 to 4 generations, at least 2 siblings per mother | 183 ka | 577 ka | Independent of Rebolledo-Jaramillo, with a bottleneck three times smaller; per transmission it is 1.7 times Rebolledo-Jaramillo’s retained rate. Nonsynonymous-to-synonymous occurrence 0.43. A model-inferred germline rate, not a count of observed homoplasmic transitions. | |||
| Connell 2022: all | 0.02619 / 345 · Sanger · 16,569 bp · 6.30 × 10-8 per site per year | Detected; mixed classHomoplasmic changes plus heteroplasmies above 20% MAF | 345 | 45 maternal lineages, 225 people, from the last four generations of an 11-generation island pedigree | 54.7 ka | 173 ka | Norfolk Island. Seven heteroplasmies and two homoplasmic changes; retaining or excluding heteroplasmy changes the estimand. Its control-region rate, 0.469/site/Myr, matches Howell 2003’s pooled 0.475. Population-specific pedigree. | |||
| Connell 2022: homoplasmic | 0.00582 / 345 · 16,569 bp · 1.40 × 10-8 per site per year | HomoplasmicThe two homoplasmic transitions only | Same 345 | Same | 246 ka | 778 ka | The best direct retained-change count here, but only two events, so its interval is broad. No bottleneck model anywhere in the derivation; its published 1.30 × 10-8 per site per year is 0.78 times the phylogenetic rate. At the paper’s 26.9-year interval, multiply ages by 1.076. Detection and family structure still require modelling. | |||
| Connell 2022: control, all | 0.01455 / 345 · 1,122 bp | — | Detected; mixed classThe nine events, by position | Same 345 | Same | 26.6 ka | 45.4 ka | The same nine-event study, partitioned by the positions in its Table 1. HVR-II ends at 371 here. One or zero events give weak constraints. The count-derived CR rate differs slightly from the printed Table 2 estimate. | ||
| Connell 2022: control, homoplasmic | 0.00291 / 345 · 1,122 bp | — | HomoplasmicThe two homoplasmic transitions, by position | Same 345 | Same | 133 ka | 227 ka | The same nine-event study, partitioned by the positions in its Table 1. HVR-II ends at 371 here. One or zero events give weak constraints. The count-derived CR rate differs slightly from the printed Table 2 estimate. | ||
| Connell 2022: HVR-I, all | 0.01164 / 345 · 360 bp | — | — | — | Detected; mixed classThe nine events, by position | Same 345 | Same | 17 ka | 31.4 ka | The same nine-event study, partitioned by the positions in its Table 1. HVR-II ends at 371 here. One or zero events give weak constraints. |
| Connell 2022: HVR-I, homoplasmic | 0.00291 / 345 · 360 bp | — | — | — | HomoplasmicThe two homoplasmic transitions, by position | Same 345 | Same | 67.9 ka | 126 ka | The same nine-event study, partitioned by the positions in its Table 1. HVR-II ends at 371 here. One or zero events give weak constraints. |
| Connell 2022: HVR-II, all | — | 0.00291 / 345 · 315 bp | — | — | Detected; mixed classThe nine events, by position | Same 345 | Same | 51.6 ka | 42.8 ka | The same nine-event study, partitioned by the positions in its Table 1. HVR-II ends at 371 here. One or zero events give weak constraints. |
| Connell 2022: HVR-II, homoplasmic | — | 00 / 345 · 315 bp | — | — | HomoplasmicThe two homoplasmic transitions, by position | Same 345 | Same | > 17.2 ka* | > 14.3 ka* | The same nine-event study, partitioned by the positions in its Table 1. HVR-II ends at 371 here. One or zero events give weak constraints. |
| Connell 2022: rest of CR, all | — | — | 00 / 345 · 447 bp | — | Detected; mixed classThe nine events, by position | Same 345 | Same | > 4.48 ka* | > 19.5 ka* | The same nine-event study, partitioned by the positions in its Table 1. HVR-II ends at 371 here. One or zero events give weak constraints. |
| Connell 2022: coding, all | — | — | — | 0.01164 / 345 · 15,447 bp | Detected; mixed classThe nine events, by position | Same 345 | Same | 89.9 ka | 341 ka | The same nine-event study, partitioned by the positions in its Table 1. HVR-II ends at 371 here. One or zero events give weak constraints. |
| Connell 2022: coding, homoplasmic | — | — | — | 0.00291 / 345 · 15,447 bp | HomoplasmicThe two homoplasmic transitions, by position | Same 345 | Same | 360 ka | 1.36 Ma | The same nine-event study, partitioned by the positions in its Table 1. HVR-II ends at 371 here. One or zero events give weak constraints. |
| Árnadóttir 2024: all | 0.07038,199 / 116,663 · 16,569 bp · 1.70 × 10-7 per site per year | Detected; branch inferredEvery mutation reaching 5% of reads, with one carrier at 10% or more | 116,663 | 64,806 living Icelanders in 2,548 matrilines (deCODE genealogy) | 20.3 ka | 64.2 ka | Seventy-five times larger than anything before it. Hypermutable positions dominate the count; the control region runs 8.3 times the coding rate per site. At deCODE’s 29.3-year interval, multiply ages by 1.172. A tight bottleneck is not a universal divisor converting this count into retained substitutions. | |||
| Árnadóttir 2024: control | 0.026 model2.38 × 10-5 per site × 1,092 callable · 1,092 bp | — | Detected; regionalThe same events in the control region | Same | Same | 14.4 ka | 24.6 ka | Callable CR denominator: 1,092 bp after exclusions. Published per-site parameter; the age uses the broader CR profile. Threshold and hypermutable sites limit its use as a clock. | ||
| Árnadóttir 2024: coding | — | — | — | 0.0443 model2.87 × 10-6 per site × 15,443 callable · 15,443 bp | Detected; regionalThe same events in the coding region | Same | Same | 23.5 ka | 89.2 ka | 16,535 attempted sites minus 1,092 CR sites gives 15,443. The age uses the broader coding profile; selection and detection thresholds affect the connection to a long-run clock. |
| Arbeithuber 2025: total | 0.74352 / 70 · ≥ 1% VAF · duplex sequencing, error < 10-7 · 16,569 bp · 1.79 × 10-6 per site per year | Detected in oocytesDe novo mutations at 1% or more of a single oocyte’s mtDNA copies, counted before any transmission | 70 | 80 single oocytes from 22 women aged 20–42, with blood and saliva from each donor | 1.92 ka | 6.07 ka | The fastest human rate measured and the purest arising rate: counted in the germ cell itself, before bottleneck, drift or transmission. No rise with maternal age in oocytes (P = 0.688) against a rise in blood (P = 0.0027). Mature oocytes have already experienced selection; this is not a retained-lineage rate. | |||
| Arbeithuber 2025: D-loop | 0.443 model3.95 × 10-4 per site per generation × 1,122 · 1,122 bp | — | Detected in oocytesThe same mutations in the D-loop | Same | Same | 0.869 ka | 1.48 ka | The authors’ own regional split. Surviving egg variants are not equivalent to homoplasmic descendant sequences. | ||
| Arbeithuber 2025: coding genes | — | — | — | 0.3 model1.95 × 10-5 per site per generation × 15,359 · 15,359 bp | Detected in oocytesThe same mutations in genes, rRNA and tRNA | Same | Same | 3.46 ka | 13.1 ka | Protein, rRNA and tRNA genes total 15,359 bp; the 88 bp of other non-coding sequence are reported separately and are not in this denominator. The age uses a coding-region approximation. |
Table sources: 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20
Calibrated on ancient DNA or the chimpanzee
| Study | Whole genome1–16,569 · 16,569 bp | What was compared | Calibration assumed | Authors’ age of Eve | Sapiens age | All-three age baseline | Method and caveats | |||
|---|---|---|---|---|---|---|---|---|---|---|
| Control region16,024–576 · 1,122 bp | Coding region577–16,023 · 15,447 bp | |||||||||
| HVR-I16,024–16,383 · 360 bp | HVR-II57–372 · 316 bp | Rest16,384–56, 373–576 · 446 bp | ||||||||
| Horai 1995: control | 0.002 model1,143 bp · 7.00 × 10-8 per site per year | — | 3 human and 4 ape complete genomes | Orangutan–African ape split at 13 Myr; the human–chimp 4.87 Myr is an output | 143 ka ± 18 (synonymous and D-loop); 450 ka ± 70 (nonsynonymous and RNA sites, same paper) | 196 ka | 335 ka | The 7.00 × 10-8 D-loop rate is itself derived from the 3.89 × 10-8 synonymous rate. 1,143-bp control-region convention; the regional distance is approximate. Other site classes give different dates. | ||
| Ingman 2000 | — | — | — | 0.00656 model15,426 bp · 1.70 × 10-8 per site per year | 53 human genomes and 1 chimpanzee; D-loop excluded | Human–chimp split at 5 Myr | 171.5 ka ± 50 | 159 ka | 602 ka | Uses the coding distance, not the 57.1 whole depth. |
| Mishmar 2003 | — | — | — | 0.00487 model15,447 bp · 1.26 × 10-8 per site per year | 104 human genomes, coding region | Human–chimp split at 6.5 Myr | 198 ka ± 19 | 214 ka | 813 ka | Selection and calibration uncertainty remain relevant; the paper argues for climatic selection while dating with a strict clock. |
| Kivisild 2006: synonymous | — | — | — | Not determined modelSubset · 3.50 × 10-8 per site per year | 277 genomes, synonymous transitions only | 6 Myr species split; 6.5 Myr mtDNA coalescence | 160 ka (SD 22) | Not calculated | Not calculated | The rate is per synonymous site. No matching synonymous-site count and root profile was established here, so neither a per-transmission figure nor an implied age is calculated. |
| Endicott & Ho 2008: external | — | — | — | 0.000687 model7,198 bp · 3.82 × 10-9 per site per year | 177 human and 2 chimpanzee genomes; 1st and 2nd codon positions (7,198 bp) analysed alone | Human–chimp split, lognormal, 95% within 5–7.5 Myr | 162 ka (122–213) from the concatenated 11,940 bp at 6.81 × 10-8; 310 ka (158–508) from these positions alone | Not calculated | Not calculated | Table 3’s separately analysed PC1+2 rate (95% HPD 2.44–5.06 × 10-9). Matching root counts for these positions were not established, so no age is calculated here. |
| Endicott & Ho 2008: internal | — | — | — | 0.00234 model7,198 bp · 1.30 × 10-8 per site per year | The same 177 genomes and software, calibration swapped | Haplogroup P in Sahul at 40–55 ka, plus H1 and H3 at about 18 ka. No ape | 108 ka (82–134) from the concatenated 11,940 bp at 9.66 × 10-8; 79 ka (60–106) from these positions alone | Not calculated | Not calculated | Table 3’s separately analysed PC1+2 rate (1.01–1.58 × 10-8), 3.4 times the external one. Not a direct pedigree measurement. |
| Soares 2009 | 0.0069 model16,569 bp · 1.67 × 10-8 per site per year | Tree of 2,196 complete genomes, corrected for purifying selection | Human–chimp split at 6.5 Myr plus 0.5 Myr lineage coalescence | 192.4 ka (Figure 6) | 207 ka | 654 ka | The benchmark every pedigree rate is compared against; a model output, not a count. The table uses its summary whole-molecule rate (one mutation every 3,624 years), not the full age-dependent function. The coding-region cross-check gives a very similar 1.708 × 10-8. | |||
| Soares 2009: HVS-I | 0.0015 model350 bp · 1.71 × 10-7 per site per year | — | — | — | The same tree, HVS-I 16051–16400 (350 bp) | Scaled against the synonymous clock | Same tree | 128 ka | 236 ka | The paper’s average of 1.62 × 10-7 and 1.784 × 10-7: “1 mutation every 16,677 years” in the segment. Its window differs from HVR-I 16024–16383; the age uses the HVR-I profile. |
| Fu 2013 | 0.0111 model16,569 bp · 2.67 × 10-8 per site per year | 10 dated ancient and 311 modern genomes | Radiocarbon-dated tips; Bayesian relaxed clock | 157 ka (120–197) | 129 ka | 408 ka | Dates, site corrections and shared ancestry matter. | |||
| Brotherton 2013 | 0.00994 model16,569 bp · 2.40 × 10-8 per site per year | 39 Neolithic and Bronze Age genomes | Holocene-only tip calibration | — | 144 ka | 454 ka | A published genome-rate summary rather than a raw homoplasmic count. | |||
| Rieux 2014 | 0.00888 model16,569 bp · 2.14 × 10-8 per site per year | 26 retained ancient and 320 modern genomes; damage screening | Tip, node and root calibration compared | 143 ka (112–180) | 161 ka | 508 ka | Its fitted mtEve age and the standardized depth division here are different calculations. | |||
| Fu 2014 | 0.0105 model16,569 bp · 2.53 × 10-8 per site per year | The Ust’-Ishim femur and 311 modern genomes | A single dated tip | — | 136 ka | 430 ka | The specimen date is part of the inference. | |||
| Posth 2016 | 0.0113 model16,569 bp · 2.74 × 10-8 per site per year | 66 dated ancient and 311 modern genomes | Dated tips; Bayesian skyline, strict rate | — | 126 ka | 397 ka | Shares some comparison data with earlier work. | |||
| Cabrera 2021: youngest bin | 0.0179 model16,569 bp · 4.33 × 10-8 per site per year | 423 dated and 784 modern genomes | Direct mutation counts by age bin, no tree | — | 79.6 ka | 251 ka | The same data, binned by age, give a rate that falls 2.3-fold from the youngest bin to the oldest. A bin-specific rate is not automatically a rate for the complete tree. | |||
| Cabrera 2021: oldest bin | 0.00791 model16,569 bp · 1.91 × 10-8 per site per year | 423 dated and 784 modern genomes | Direct mutation counts by age bin, no tree | — | 180 ka | 570 ka | The same data, binned by age, give a rate that falls 2.3-fold from the youngest bin to the oldest. A bin-specific rate is not automatically a rate for the complete tree. | |||
| Översti & Palo 2022: U4 | 0.0365 model16,569 bp · 8.80 × 10-8 per site per year | 42 Holocene U4 genomes, no moderns | Root-to-tip regression on dated tips | — | 39.2 ka | 124 ka | The fastest defensible number in the mainstream literature, and the best single figure available to the young-earth side. It is also 42 sequences from one twig of the tree with R² = 0.37, and its own x-intercept still puts that twig’s root at 12,723 years. Not an established universal mtDNA rate. | |||
| Project: U-subclade regression | 0.0107 model16,519 bp · 2.59 × 10-8 per site per year | 229 dated ancient genomes | Root-to-tip, per-subclade intercepts | — | 133 ka | 421 ka | An in-house descriptive calculation. Independent-error intervals omit phylogenetic covariance and date uncertainty. | |||
Table sources: 21 22 23 24 25 26 27 28 29 30 31 32 33 2
The two Endicott and Ho rows are the experiment. Same sequences, same alignment, same software, one thing changed, and the answer moves from 162,000 years to 108,000. On first-and-second codon positions alone the same swap moves it from 310,000 to 79,000, a fourfold spread from a methodological choice. Anyone who says the ape calibration does no work is wrong.25
The distances the rates have to explain
Every rate above has to explain these distances. They are counts of differences, not rates. Each panel’s first line is the whole molecule: the mean over all pairs, its maximum, and the control-region total. The second line splits the mean by region, with maxima. A dagger marks the largest between-group mean in each region; half of it is the all-three column’s depth. The Behar row is a root-to-tip depth, a different quantity from a pairwise mean, and its second line is this project’s regional allocation used by the sapiens column. The Soares column converts each mean to the age of the average pair’s common ancestor at 1.665 × 10-8 per site per year, with every tip treated as living, so it sits below the root by construction.261343536
| Panel | Whole genome1–16,569 · 16,569 bp | What is measured | Pairs | Sample | Average pair’s ancestor at Soares’s rate | Notes | |||
|---|---|---|---|---|---|---|---|---|---|
| Control region16,024–576 · 1,122 bp | Coding region577–16,023 · 15,447 bp | ||||||||
| HVR-I16,024–16,383 · 360 bp | HVR-II57–372 · 316 bp | Rest16,384–56, 373–576 · 446 bp | |||||||
| All living humans54 genomes | 59.40mean of all pairs · max 106 · control region 15.88 (max 33) | Pairwise differences within the group | 1,431 | rCRS plus Ingman’s 53 globally sampled genomes | 108 ka | Africa-heavy by design, which is why its numbers are higher than Kim and Schuster’s. Half the average pair is a pair-specific branch length, not the root-to-tip depth. | |||
| 8.19max 19 | 5.76max 16 | 1.92max 5 | 43.52max 87 | ||||||
| Behar root to tip18,843 genomes | 57.1mean · range 41–77 | Substitutions from the reconstructed root (RSRS) to each living genome | — | 18,843 complete genomes | 207 ka | The sapiens-age column’s whole-molecule anchor; median 56, SD 5.9. A 6,000-year root needs 5.74 × 10-7 per site per year: 26.8 times Rieux 2014’s rate and 34.5 times Soares 2009’s. The second line is this project’s allocation of the 57.1 by the RSRS profile of the 54 genomes (scaled by 57.1/50.241), not a figure Behar published. | |||
| 7.87allocated | 5.98allocated | 1.56allocated | 41.69allocated | ||||||
| Kim and Schuster7,098 genomes | 39.8Π · max 123 | Nucleotide diversity Π: differences between two random genomes | — | 7,098 complete genomes, 10% African | 72.1 ka | A Europe- and Asia-heavy panel. Its control-region and coding figures are Π over 1,122 and 15,447 bp; no HVR split is published. | |||
| 9.9max 32 | 29.9max 98 | ||||||||
| All Neanderthals36 records, 35 individuals | 31.49mean of all pairs · max 118 · control region 10.40 (max 29) | Pairwise differences within the group | 630 | Every near-complete Neanderthal mitogenome in GenBank at the time of the relative-rate test | 57.1 ka | About half the modern human spread. | |||
| 5.84max 20 | 3.57max 10 | 0.98max 6 | 21.09max 93 | ||||||
| All Denisovans7 records | 42.76mean of all pairs · max 86 · control region 13.76 (max 25) | Pairwise differences within the group | 21 | Seven mitogenomes | 77.5 ka | Seven records only, so the maximum means little. Wider than the Neanderthal spread despite the smaller sample. Denisova 19 and 21 have identical mtDNA. | |||
| 8.33max 14 | 4.67max 9 | 0.76max 2 | 29.00max 61 | ||||||
| Human with Neanderthalbetween groups | 199.35mean of all pairs · max 225 · control region 35.06 (max 45) | Mean pairwise differences between the two groups | 1,944 | The two panels above | 361 ka | Spread over 6,000 years at 25-year generations, it needs 0.415 retained changes per lineage per generation, with tips treated as contemporaneous. | |||
| 20.50max 30 | 9.03max 16 | 5.52max 9 | 164.29max 182 | ||||||
| Human with Denisovanbetween groups | 360.83†mean of all pairs · max 391 · control region 44.42 (max 56) | Mean pairwise differences between the two groups | 378 | The two panels above | 654 ka | Spread over 6,000 years at 25-year generations, it needs 0.752 retained changes per lineage per generation, with tips treated as contemporaneous. | |||
| 25.44max 32 | 8.70max 16 | 10.27max 14 | 316.42†max 341 | ||||||
| Neanderthal with Denisovanbetween groups | 349.79mean of all pairs · max 378 · control region 52.62† (max 57) | Mean pairwise differences between the two groups | 252 | The two panels above | 634 ka | Spread over 6,000 years at 25-year generations, it needs 0.729 retained changes per lineage per generation, with tips treated as contemporaneous. | |||
| 29.13†max 34 | 9.92†max 12 | 13.58†max 17 | 297.17max 322 | ||||||
| Human with Sima de los Huesosbetween groups | 279.33mean of all pairs · max 296 · control region 20.15 (max 27) | Mean pairwise differences between the groups | 54 | 54 living humans against the one Sima sequence | 506 ka | One sequence only. Not used by the all-three column, which takes the three groups with several sequences each. | |||
| 7.17max 12 | 7.98max 13 | 5.00max 7 | 259.19max 269 | ||||||
| All three pooled97 genomes | 154.17mean of all pairs · max 391 · control region 27.44 (max 57) | Pairwise differences across all 97 | 4,656 | The three panels above | 279 ka | Uncorrected pairwise figures, not model split dates. Half the pooled average underestimates the root, because many pairs share younger ancestors. | |||
| 15.55max 34 | 7.29max 16 | 4.60max 17 | 126.73max 341 | ||||||
Accounting for a fossil's age
A fossil has had less time to accumulate mutations than a living person. With tip ages a and b years before present, a pairwise distance D, and a lineage rate R in changes per region per year, the constant-clock calculation is:
Common-ancestor age = (D / R + a + b) / 2
The correction needs the dates of the actual sequences being compared. An age assigned to one Denisovan cannot simply be added to the average distance across seven Denisovan records. Their ages and contexts differ, and Denisova 19 and 21 have identical mtDNA and may represent one individual or maternal relatives.37
One check can be made without assigning dates to the whole panel. Denisova 19's sequence differs from the 54 modern sequences by an average of 351.111 whole-molecule substitutions under the stated mask. Brown and colleagues place its depositional layer at 187–217 ka in their 95.4% interval. Substituting that layer-age range produces the following human-Denisova-19 estimates. These are one dated-pair comparison, not the fitted root of all three groups.372
| Rate or proposal | Human–Denisova-19 ancestor | Ancestor older than fossil? |
|---|---|---|
| Rebolledo-Jaramillo 2014: retained | 1069–1084 ka | Yes throughout the layer interval |
| Zaidi 2019 | 655–670 ka | Yes throughout the layer interval |
| Connell 2022: homoplasmic | 851–866 ka | Yes throughout the layer interval |
| Rieux 2014 | 588–603 ka | Yes throughout the layer interval |
| Jeanson 2015 synthesis, for comparison | 121–136 ka | No: incompatible plug-in result |
| Árnadóttir 2024: all | 156–171 ka | No: incompatible plug-in result |
| Arbeithuber 2025: total | 99–114 ka | No: incompatible plug-in result |
The three fastest rows yield a calculated ancestor younger than the fossil’s assigned layer age. That is an invalid ancestor-descendant ordering for the plug-in calculation, exposing a failure of that combination of rate and dating assumptions. It does not show that an observed low-frequency mutation count was wrong. A young-earth account can dispute the geological chronology; if it does, the sequence distances alone do not supply the missing years. Conversely, accepting the layer chronology while using those raw counts as a constant substitution clock fails this check.
What these comparisons establish
Rebolledo-Jaramillo's retained-change model gives about 317 ka for sapiens and a 1 Ma all-three baseline; Zaidi's model gives about 183 ka and 577 ka. Connell's two homoplasmic changes give about 246 ka and 778 ka at 25 years per generation. At Connell's measured 26.9-year interval, both ages rise by 7.6%. These are the more relevant family-study quantities for a persistent-lineage clock, but the model-based results are not direct counts of fixation, and Connell's count is only two.131718
Uncertainty is especially large for rare events. Using Connell's own annual homoplasmic rate interval, 0.003–0.047 changes/site/Myr, gives a standardized sapiens interval of roughly 73 ka to 1.15 Ma. That propagates the paper's rate interval only; it omits root-depth and model uncertainty. A point estimate such as 246 ka must not conceal that width.18
Calibrated ancient-DNA estimates generally give a similar order of magnitude. Faster results remain in the table: Översti and Palo's U4 regression gives 8.8 × 10-8 per site per year, and Cabrera's youngest bin gives 4.33 × 10-8. They deserve attention, but a fast rate in one clade or age bin is not automatically valid throughout the human tree. No row by itself settles the question; the useful test is whether one mutation-and-transmission model can explain these different observations together.3233
Superseded: the per-site-per-year table, kept for auditing
This is the older article's second table, which normalised each study to a rate per site per year and then to a 25-year generation. The tables above replace it, and it is kept only so its numbers can be audited. Three things differ from the tables above. Its per-site-per-year figures carry each study's own generation interval and length convention, so they do not divide back to the counts. Its two Wei 2019 rows (0.26 and 0.0016 per 25-year generation, from 6.3 × 10-7 and 4.0 × 10-9 per site per year) disagree with Wei's own per-transmission counts, 477/1,526 = 0.313 and 3/1,526 = 0.00197; use the counts. And its Rebolledo-Jaramillo row is the retained-rate model, not a count.
| Study | Region sequenced | Design | What was counted | per site per year | per 25-yr generation |
|---|---|---|---|---|---|
| Howell 19963 | control region, 1,122 bp | 1 LHON pedigree, 88 transmissions | new heteroplasmies | 1.01 × 10-6 | 0.028 |
| Parsons 19975 | HV1+HV2, ~610 bp | 134 lineages, 327 transmissions | both, on Howell's re-scoring | 2.5 × 10-6 | 0.038 |
| Cavelier 20006 | control region + 500 bp coding | 33 lineages, 292 transmissions | fixed only — none observed | < 2.1 × 10-7 | < 0.006 |
| Sigurðardóttir 20007 | HV1+HV2, ~660 bp | 26 Icelandic pedigrees, 705 transmissions | substitutions between relatives | 3.2 × 10-7 | 0.0053 |
| Heyer 20018 | control region, 673 bp | French-Canadian deep pedigrees, 508 transmissions | substitutions between relatives | 3.9 × 10-7 | 0.0066 |
| Forster 20029 | control region, 579 bp | Kerala, 200 low-radiation transmissions | new heteroplasmies only | 4.3 × 10-7 | 0.0063 |
| Howell 20034 | control region, 1,122 bp | nine studies pooled, 2,633 transmissions | everything newly arising | 4.72 × 10-7 | 0.0132 |
| Santos 200510 | HVI+HVII, ~610 bp | Azores, 321 transmissions | germinal heteroplasmies, model-corrected | 2.42 × 10-7 | 0.0037 |
| Madrigal 201211 | HVR-I only, 360 bp | Costa Rica, 19 pedigrees, 289 meioses | substitutions between relatives | 0.89–2.51 × 10-6 | 0.008–0.023 |
| Guo 201312 | whole genome, 3,981x depth | 26 mother–child transmissions | fixed only — none observed | < 2.3 × 10-7 | < 0.096 |
| Rebolledo-Jaramillo 201413 | whole genome, ~20,000x | 39 mother–child pairs, two tissues | heteroplasmies, converted to fixation | 1.3 × 10-8 | 0.0056 |
| Ding 201514 | whole genome, ~180x | 333 Sardinian trios, lymphocytes | homoplasmic + heteroplasmic discordance | 8.5 × 10-7 | 0.35 |
| Wei 2019, arising16 | whole genome, ~1,900x | 1,526 transmissions | new heteroplasmies above 1% | 6.3 × 10-7 | 0.26 |
| Wei 2019, near-fixed in the child16 | whole genome, ~1,900x | the same 1,526 transmissions | de novo variants reaching >80% | 4.0 × 10-9 | 0.0016 |
| Zaidi 201917 | whole genome, ~11,400x | 96 families, 249 transmissions | germline heteroplasmy, converted to fixation | 1.57 × 10-8 | 0.0065 |
| Connell 2022, all18 | whole genome, Sanger | Norfolk Island, 345 transmissions | fixed + heteroplasmic | 5.85 × 10-8 | 0.024 |
| Connell 2022, fixed18 | whole genome, Sanger | the same 345 transmissions | fixed only — 2 observed | 1.30 × 10-8 | 0.0054 |
| Árnadóttir 202419 | whole genome | deCODE, 116,663 transmissions | every mutation reaching 5% of reads | 1.41 × 10-7 | 0.059 |
| Arbeithuber 202520 | whole genome, duplex | 80 single oocytes, 70 mother-to-oocyte transmissions | new mutations at 1% or more, inside the egg | 1.8 × 10-6 | 0.74 |
| Jeanson 2015, for comparison15 | whole genome | three studies pooled, 398 transmissions | treated as fixed | 3.8 × 10-7 | 0.158 |
Sources
- 1.
Behar, Doron M. et al. “A 'Copernican' Reassessment of the Human Mitochondrial DNA Tree from its Root.” American Journal of Human Genetics 90(4):675–684. 2012. The root lies between L0 and L1'2'3'4'5'6, the paper says they are. “separated from each other by 14 coding and four control-region mutations.” A reassessment of the human mitochondrial DNA family tree that proposes a new reference sequence rooted on Neanderthal mitogenomes, the Reconstructed Sapiens Reference Sequence. The abstract says the phylogeny was refined by adding information from 8,216 modern mitogenomes, while the mean distance below is reported over a larger set of 18,843, so the two numbers refer to different things. The mean distance from the reconstructed root to a modern tip is 57.1 substitutions (standard deviation 5.9, range 41-77), and the paper also reports statistically significant molecular-clock violations in M and L2, which is a reason for caution about every published interval on either side. The article uses it for the 57.1-substitution mean depth of the tree, its 41-77 range, and the bifurcation at the root. ↩ ↩ ↩
- 2.
Region-matched calculations for the rate tables: rate inputs and results, full-precision results, the generator and its independent checks. Each age divides a region's root-to-tip depth by the study's rate per transmission at 25-year generations. The sapiens depth is Behar's 57.1 for the whole molecule, allocated to smaller windows by the root-to-tip profile of 54 living genomes measured from the reconstructed root and scaled by 57.1/50.241 (counts, code). The all-three baseline is half the largest between-group mean distance in the same window. The U-subclade regression row comes from a separate script and its results. ↩ ↩ ↩ ↩
- 3.
Howell, Neil et al. “How rapidly does the human mitochondrial genome evolve?” American Journal of Human Genetics 59(3):501–509. 1996. Page 501: “triplasmic descendants with three mitochondrial genotypes”. Both mutations produced these descendants, which are heteroplasmic (carry more than one mitochondrial type), not fixed. A pedigree study that opened the pedigree-rate literature. It sequenced the D-loop (control region) of 45 members of a large matrilineal Leber hereditary optic neuropathy pedigree and found two germline mutations in 88 transmissions, normalised here as 2/88 over 1,122 bp of control region. The article uses it as one row of its table of pedigree rates. ↩ ↩
- 4.
Howell, Neil et al. “The pedigree rate of sequence divergence in the human mitochondrial genome: there is a difference between phylogenetic and pedigree rates.” American Journal of Human Genetics 72(3):659–670. 2003. It counts arising, not fixing, in the paper's own words: “All newly arising substitutions, both heteroplasmic and homoplasmic, were included in our divergence-rate estimates.” A pooled analysis of nine control-region pedigree studies, counting 28 mutations in 2,633 transmissions and counting every newly arising change, whether it sat in some of a person's mitochondrial copies or in all of them. Its figure of 0.95 mutations per base pair per million years is a divergence rate, twice the mutation rate, which is verifiable from the paper's own table; halving it and keeping it inside the 1,122 base pairs of the control region gives the 4.72 × 10-7 tabulated in the article. It is the canonical pooled control-region figure, and the article relies on it. Multiplying the published figure by 16,569 sites and doubling it, as Carter 2025 endnote 10 does, inflates the answer roughly thirty-fold. ↩ ↩
- 5.
Parsons, Thomas J. et al. “A high observed substitution rate in the human mitochondrial DNA control region.” Nature Genetics 15(4):363–368. Apr 1997. A genetic study that sequenced 610 base pairs of two hypervariable segments (HV1 and HV2) of mitochondrial DNA directly in 327 mother-to-child transmissions, not by RFLP analysis as one online rebuttal says. It found a substitution rate of 2.5 per site per million years, about 20 times the rate predicted by phylogenetic analyses. The article uses it because its authors worked out the age that this fast rate gives, about 6,500 years, and called it clearly incompatible with the known age of modern humans. ↩ ↩
- 6.
Cavelier, Lucia et al. “MtDNA substitution rate and segregation of heteroplasmy in coding and noncoding regions.” Human Genetics 107(1):45–50. 2000. The study reports: “No homoplasmic nucleotide substitutions were detected in a total of 292 transmissions.” The forgotten null result, and the first study to count fixed changes specifically. The article uses it as the no-fixed-changes row of its table of pedigree rates. ↩ ↩
- 7.
Sigurðardóttir, Sigrún et al. “The mutation rate in the human mtDNA control region.” American Journal of Human Genetics 66(5):1599–1609. 2000. A pedigree study of the control region of mtDNA in 272 people from 26 large Icelandic family trees, covering 705 mother-to-child transmissions, which found 3 substitutions. Howell 2003 re-scored the same dataset as 5 by counting newly arising heteroplasmies as well, which shows the arising-versus-fixing distinction inside a single dataset. The article uses it as one row of its table of pedigree rates. ↩ ↩
- 8.
Heyer, Evelyne et al. “Phylogenetic and familial estimates of mitochondrial substitution rates: study of control region mutations in deep-rooting pedigrees.” American Journal of Human Genetics 69(5):1113–1126. 2001. The authors write that hypermutable sites. “are expected to be of no or little informativity for the evaluation of mutational distances at the phylogenetic time scale.” A pedigree study of control region mutations in deep-rooting French-Canadian family trees, which found 4 substitutions in 508 maternal transmissions over 673 bp. Its more important result is the pooled one: of 18 substitutions in 1,729 transmissions, only 12 were at independent sites, because three positions had each mutated twice. The article uses it as the direct empirical rebuttal to any argument that extrapolates a control-region pedigree count backwards. ↩ ↩
- 9.
Forster, Lucy et al. “Natural radioactivity and human mitochondrial DNA mutations.” PNAS 99(21):13950–13954. 2002. A pedigree study of 248 families in Kerala, covering 579 bp of the control region, compared with families on low natural radiation. It counts new heteroplasmies explicitly and exclusively: 22 in 595 high-radiation transmissions against 1 in 200 low-radiation transmissions, and none of the 23 was fixed. The article uses the low-radiation arm as its figure, and uses the study as its test of an environmental insult: a severe natural exposure raised arising mutations without fixing any. ↩ ↩
- 10.
Santos, Cristina et al. “Understanding differences between phylogenetic and pedigree-derived mtDNA mutation rate: a model using families from the Azores Islands (Portugal).” Molecular Biology and Evolution 22(6):1490–1505. 2005. The authors write: “Eleven new point substitutions were found, all of them in heteroplasmy, from which five appear to be somatic mutations and six can be considered germinal.” The last-named author is Aluja, María Pilar. A pedigree study of Azores families, built as a model to explain why phylogenetic and pedigree mutation rates differ. It is the first study to correct an arising count into a fixing rate explicitly, and zero reached homoplasmy. The tabulated rate is the corrected figure, a model output; the raw data are 6 germinal mutations in 321 transmissions. The article uses it as its model-corrected Azores row. ↩ ↩
- 11.
Madrigal, Lorena et al. “High mitochondrial mutation rates estimated from deep-rooting Costa Rican pedigrees.” American Journal of Physical Anthropology 148(3):327–333. 2012. A pedigree study of 19 deep-rooting genealogies from Costa Rica, covering 289 meioses (transmissions from parent to child), that found 7 mutations in hypervariable region I, the fast-changing part of the control region. The paper uses a measured female generation interval of 28.3 years, which is worth noting against the 15 years the young-earth calculation requires; its abstract headline rate of 1.24 × 10-6 per site per year could not be reproduced from its own Results or Table 2, which give 2.51 × 10-6 (maximum) and 0.89 × 10-6 (minimum). The article quotes its authors' caution that the fast and slow estimates should both be considered. ↩ ↩
- 12.
Guo, Yan et al. “Very low-level heteroplasmy mtDNA variations are inherited in humans.” Journal of Genetics and Genomics 40(12):607–615. 2013. Their finding, in their words: “Homoplasmic SNPs were inherited without exception within each family.” A sequencing study of 26 mother-child pairs at a median depth of 3,981-fold, more than twenty times deeper than the Sardinian data the young-earth rate is built on. It found zero fixed changes in 26 transmissions, where Jeanson's rate predicts 4.1; the article uses it as a test of that rate. ↩ ↩
- 13.
Rebolledo-Jaramillo, Boris et al. “Maternal age effect and severe germ-line bottleneck in the inheritance of human mitochondrial DNA.” PNAS 111(43):15474–15479. 2014. A study of 39 mother-child pairs, using blood and buccal (cheek) tissue sequenced at 20,000-fold coverage with a 1% detection threshold, that asks whether a mother's age changes the mitochondrial DNA she passes on. It estimated the germline bottleneck at 30-35 mitochondrial genomes. The article uses it for the bottleneck formula that separates mutations that arise from mutations that become fixed, and for its maternal-age numbers. ↩ ↩ ↩
- 14.
Ding, Jun et al. “Assessing Mitochondrial DNA Variation and Copy Number in Lymphocytes of ~2,000 Sardinians Using Tailored Sequencing Analysis Tools.” PLoS Genetics 11(7):e1005306. 2015. A study that sequenced lymphocytes (white blood cells) from about 2,000 Sardinians at roughly 180-fold coverage and reported how often mothers and children share mitochondrial variants, which is concordance, not a mutation rate. It is the study Jeanson's whole-genome rate rests on, and the article uses it to show that his rate repurposes a concordance figure as a mutation count. ↩ ↩
- 15.
Jeanson, Nathaniel T. “A Young-Earth Creation Human Mitochondrial DNA 'Clock': Whole Mitochondrial Genome Mutation Rate Confirms D-Loop Results.” Answers Research Journal 8:375–378. 2015. A young-earth creationist article that builds a whole-genome mutation rate by pooling three whole-genome parent-child studies into one count of transmissions and mutations. Its Supplemental Table 1, with formulas intact, reproduces every figure in Table 1 exactly when recomputed from its own inputs, so there is no arithmetic error in it. The article uses it as the whole-genome step of Jeanson's clock argument. Mirror: Assets.answersresearchjournal.org ↩ ↩
- 16.
Wei, Wei et al. “Germline selection shapes human mitochondrial DNA diversity.” Science 364(6442):eaau6520. 2019. A genetic study that sequenced 1,526 mother-offspring pairs and counted the new heteroplasmies in each child. Supplementary Table S2 lists 13 de novo D-loop variants, but the rate printed on the same row requires 131; every other row is exactly consistent with n = 1,526, and this one implies n = 151. It is a dropped digit. The article uses it as its main direct comparison of mutations appearing with mutations staying. Mirror: Repository.cam.ac.uk ↩ ↩ ↩
- 17.
Zaidi, Arslan A. et al. “Bottleneck and selection in the germline and maternal age influence transmission of mitochondrial DNA in human pedigrees.” PNAS 116(50):25172–25178. 2019. A genetic study of heteroplasmy, the mix of mutant and normal mitochondrial DNA inside one person, in 96 multigenerational families, with two tissues per person and a mean sequencing depth of about 11,400-fold. It estimated a germline bottleneck of 10.3 segregating units (95% HPD 8.5–12.3) and a germline rate of 4.72 × 10-7 per site per generation. The heteroplasmy count is a measurement; the bottleneck correction that turns it into a fixation rate is a model. The article uses it for the one measured lever on germline mutation, maternal age, and for the bottleneck that maternal age tightens. ↩ ↩ ↩
- 18.
Connell, Jennifer R. et al. “Pedigree derived mutation rate across the entire mitochondrial genome of the Norfolk Island population.” Scientific Reports 12:6827. 2022. A pedigree study that sequenced the complete mitochondrial genome of 225 individuals from 45 maternal lineages in the documented Norfolk Island genealogy, covering 345 mother-child transmissions, and found nine changes, seven of them heteroplasmic (a variant present in only some of a person's mitochondrial copies) and two fixed changes. Counting all nine gives 0.058 mutations per site per million years; excluding the seven heteroplasmic changes gives 0.013, that is 1.30 × 10-8 per site per year, against the phylogenetic 1.665 × 10-8. Jeanson's rate of 0.158 per generation predicts 54.6 fixed changes in 345 transmissions, and the probability of observing 2 or fewer is about 3 × 10-21; this is the cleanest direct measurement of the whole-genome fixing rate in a deep human pedigree and the most important source in the article. ↩ ↩ ↩ ↩ ↩
- 19.
Árnadóttir, Erla R. et al. “The rate and nature of mitochondrial DNA mutations in human pedigrees.” Cell 187(15):3904–3918.e8. 2024. A study from deCODE Genetics/Amgen in Reykjavik of 64,806 Icelanders from 2,548 matrilines, with 116,663 mother-child transmissions and 8,199 mutations, which the article calls the largest such study ever run. The article gives the rate as 0.0703 mutations per genome per transmission. The article uses it as its direct measurement of the pedigree rate, divided into the measured depth of the human mitochondrial tree, which gives roughly 24,000 years. Mirror: Decode.com ↩ ↩
- 20.
Arbeithuber, Barbara et al. “Allele frequency selection and no age-related increase in human oocyte mitochondrial mutations.” Science Advances 11(32):eadw4954. Aug 6, 2025. The oocytes came from 28 eggs of eight women aged 20-29, 25 from six aged 30-34, and 27 from eight aged 35-42. The age result, from the abstract: “with age, mutations increased in blood and saliva but not in oocytes … mtDNA in human oocytes is protected against accumulation of mutations with aging and having functional consequences.” The mixed model gives Z-test P = 2.7 × 10-3 for blood, 5.7 × 10-2 for saliva and 0.688 for oocytes; the permutation test on old versus young oocytes gives P = 0.550 and a 0.99-fold decrease. The paper's count is 52 de novo mutations at minor allele frequency >= 1% over. “70 mother-to-oocyte transmissions … the germline mutation rate to be 3.95 × 10-4 … and 1.95 × 10-5 … mutations per site per generation for the D-loop and coding region, respectively.” The per-genome figure quoted in this article is computed from those two published rates over the regions they cover: 3.95 × 10-4 × 1,122 = 0.443 for the D-loop, and 1.95 × 10-5 × 15,359 = 0.300 for what the paper calls the coding region (protein-coding 11,341 bp, rRNA 2,513 bp and tRNA 1,505 bp). The 88 bp of non-coding sequence outside the D-loop is reported separately and omitted here. Together they give 0.743 mutations per genome per generation over the 16,481 bp rated, or 4.5 × 10-5 per site per generation and 1.8 × 10-6 per site per year at 25 years. The implied Eve age divides Behar's 57.1 root-to-tip by that figure. A 2025 study that applied duplex sequencing, whose error rate is below 10-7, to 80 single human eggs (oocytes) from 22 women aged 20 to 42, with blood and saliva from the same donors. The article uses it as the only direct measurement of the quantity a young-earth rescue would need, and it points the wrong way for that rescue. Mirror: PMC ↩ ↩
- 21.
Horai, Satoshi et al. “Recent African origin of modern humans revealed by complete sequences of hominoid mitochondrial DNAs.” PNAS 92(2):532–536. 1995. A research paper that sequenced the complete mitochondrial DNA of three humans, three African apes and an orangutan to estimate how fast hominoid mtDNA changes and when modern humans shared a common ancestor. The calibration is two steps removed from anything datable: the orangutan and African ape split is fixed at 13 million years, and the human-chimpanzee split of 4.87 ± 0.23 million years (Myr) comes out of that. The paper's own data give two irreconcilable dates for the human ancestor: 143,000 ± 18,000 years from synonymous sites and the D-loop, and 450,000 ± 70,000 from nonsynonymous and RNA-gene sites. Only the first is ever quoted. The article uses it in its table of rate estimates, where it gives both dates. ↩
- 22.
Ingman, Max et al. “Mitochondrial genome variation and the origin of modern humans.” Nature 408(6813):708–713. 2000. “has not evolved at a constant rate across all human lineages”. “consequently less suitable for dating evolutionary events”. A genetic study of 53 complete mitochondrial genomes, the whole molecule, from people around the world. The mean pairwise difference over the whole molecule was 61.1 for all humans (n=53), 38.5 for non-Africans (n=32) and 76.7 for Africans (n=21), and the authors left out the D-loop, a control region, because it does not evolve at a constant rate. The article uses it for the rate that depends on an assumed 5-million-year human-chimpanzee split, and for its whole-genome sequences, which first re-established the African root. ↩
- 23.
Mishmar, Dan et al. “Natural selection shaped regional mtDNA variation in humans.” PNAS 100(1):171–176. 2003. The last-named author is Wallace, Douglas C. A study of 104 complete mtDNA genomes from all global regions, arguing that natural selection, with climate as one influence, shaped regional mtDNA variation. The measured quantity is a ratio: human tree depth to the human-chimpanzee ancestor, 24.88 against 818.05 × 10-4 substitutions per nucleotide, or 1 : 32.8, and the 198 ka date comes from fixing the denominator at 6.5 Myr. Its separate claim of climate-driven selection on ATP6 rests on counts as small as 3 versus 7 nonsynonymous changes and was not reproduced by Kivisild 2006. The article uses it for its 198 ka row, the date that the human-chimpanzee calibration gives. ↩
- 24.
Kivisild, Toomas et al. “The role of selection in the evolution of human mitochondrial genomes.” Genetics 172(1):373–387. 2006. A 2006 study of complete coding-region sequences from 277 people on five continents, arranged as a family tree of mitochondrial lineages. Of the 1,788 mutations in the tree, 1,758 occurred at 1,292 variable coding-region sites, and the surplus of amino-acid-changing mutations is a general feature of the young branches, so the ratio of amino-acid-changing to silent variation declines with the age of the branch. It did not replicate the climate-based regional differences in amino-acid-changing mutations that Ruiz-Pesini et al. 2004 claimed, although Mishmar et al. 2003 is among the studies it cites for the claim that natural selection has significantly shaped mtDNA evolution. Its coalescent date for the tree is 160,000 years (SD 22,000). The article reads the observed decline as variants carried by healthy living people still being removed, and uses it for the claim that mitochondrial sites are not neutral. ↩
- 25.
Endicott, Phillip, and Simon Y. W. Ho. “A Bayesian evaluation of human mitochondrial substitution rates.” American Journal of Human Genetics 82(4):895–902. 2008. A Bayesian phylogenetic analysis of 177 human and two chimpanzee mitochondrial genomes that dated the human mitochondrial ancestor twice on the same data, once with a human-chimpanzee calibration and once with dates of human settlement inside the tree. The swap moved the answer by half, which is why the article treats every ape-calibrated date as conditional on its calibration. ↩ ↩
- 26.
Soares, Pedro et al. “Correcting for Purifying Selection: An Improved Human Mitochondrial Molecular Clock.” American Journal of Human Genetics 84(6):740–759. 2009. “one mutation every 3624 years”. “a severe shortage of such uncontroversial episodes”. A study that recalibrated the human mitochondrial clock for the whole molecule, using more than 2,000 complete mitochondrial genomes. Its rate of 1.665 × 10-8 substitutions per nucleotide per year is calibrated on an assumed 6.5–7 million-year human-chimpanzee split, so it is a model output rather than a measurement, and the young-earth circularity critique is correct about this paper specifically. The paper also concedes that there are few uncontroversial episodes available for archaeological calibration. The article uses it for the phylogenetic rate that the young-earth argument has to overturn, and for the size of the selection correction, which the paper's own curve puts at 1.50 times at 0–1 thousand years falling to 1.13 times at 50–100 thousand years. ↩ ↩
- 27.
Fu, Qiaomei et al. “A revised timescale for human evolution based on ancient mitochondrial genomes.” Current Biology 23(7):553–559. 2013. A dating study that used ten radiocarbon-dated ancient mitochondrial genomes, spanning 40,000 years, as tip calibrations (samples of known age that anchor the clock), with no ape calibration anywhere in it. It gives a whole-genome rate of 2.67 × 10-8 per site per year (95% HPD, the range the model gives the rate with 95% probability, 2.16-3.16) and a mitochondrial TMRCA, the time back to the most recent common ancestor, of 157,000 years (120,000-197,000); the same data analysed by linear regression give 1.92 × 10-8, a 39% difference, and the regression's own confidence interval spans a 2.3-fold range on nine points. The article uses it for its calibrated rate in the comparison table and for the Cro-Magnon 1 redating quoted below. ↩
- 28.
Brotherton, Paul et al. “Neolithic mitochondrial haplogroup H genomes and the genetic origins of Europeans.” Nature Communications 4:1764. 2013. The last-named authors are Cooper, Alan and the Genographic Consortium. A study of 39 complete haplogroup H genomes from Central European skeletons dated 5500-2000 BC. Every calibration point sits inside the dendrochronologically anchored part of the radiocarbon curve, with no Pleistocene sample, no ape split and no Out-of-Africa node, and a date-randomisation test confirmed the dates carry enough signal. The authors note their 2.4 × 10-8 is 1.45x the ape-calibrated 1.66 × 10-8. The article uses it as its Holocene-only calibration row. ↩
- 29.
Rieux, Adrien et al. “Improved calibration of the human mitochondrial clock using ancient genomes.” Molecular Biology and Evolution 31(10):2780–2792. 2014. A recalibration of the human mitochondrial clock using 30 high-quality ancient whole mitochondrial genomes and 320 contemporary ones, of which 26 ancient genomes were kept after damage screening. It gives a whole-genome rate of 2.143 × 10-8 per site per year and a mitochondrial time to the most recent common ancestor of 143,000 years (112,000-180,000). One error to note: Table S3 lists the six Blätterhöhle samples' ages in calibrated BC while labelling the column "Age (CalBP)," so each is exactly 2,000 years too young. The article uses it for the tip-calibrated clock and its tests of rate constancy, damage and calibration. ↩
- 30.
Fu, Qiaomei et al. “Genome sequence of a 45,000-year-old modern human from western Siberia.” Nature 514(7523):445–449. 2014. A genome study of one femur from Ust'-Ishim in western Siberia, dated by two independent radiocarbon measurements of 41,400 ± 1,300 and 41,400 ± 1,400 BP, which calibrate to 46,880-43,210 cal BP. Its mitochondrial rate is 2.53 × 10-8 per site per year (95% HPD 1.76-3.23), and the abstract's rounded interval of 1.8 to 3.2 × 10-8 is not a point estimate; the paper's own supplement gives the rate of Fu 2013 as 2.47 × 10-8, where the published value is 2.67. The article uses it as the single dated tip of its rate table and as its answer to the claim that fossil DNA is too degraded to use. ↩
- 31.
Posth, Cosimo et al. “Pleistocene mitochondrial genomes suggest a single major dispersal of non-Africans and a Late Glacial population turnover in Europe.” Current Biology 26(6):827–833. 2016. The last-named author is Krause, Johannes. A study that recovered 35 new mitochondrial genomes from Late Pleistocene and early Holocene European hunter-gatherers. It used 66 radiocarbon-dated ancient genomes as tip calibrations alongside 311 modern ones, spanning 3,600 to 46,880 cal BP, with roughly a third of the calibration points Holocene. Four model combinations give 2.62, 2.67, 2.74 and 2.77 × 10-8; the preferred Bayesian-skyline strict-rate model gives 2.74 × 10-8 (2.44-3.01 HPD). The article uses it as its table row for that model. ↩
- 32.
Cabrera, Vicente M. “Human molecular evolutionary rate, time dependency and transient polymorphism effects viewed through ancient and modern mitochondrial DNA genomes.” Scientific Reports 11:5036. 2021. A direct count of mutation differences between 423 radiocarbon-dated ancient genomes and 784 modern ones, binned into ten age periods, with no phylogenetic model at all. Mean differences rise with age, almost monotonically (one small dip, from 1.63 to 1.51, at the fourth period): 0.63 at 1,119 BP, 3.28 at 8,504, 11.07 at 29,638, 14.39 at 40,160, which is arguably the cleanest single table in the subject. Flagged discrepancy: the rates in its Table 2 cannot be reproduced from Table 1 by the paper's own stated formula, and reverse-engineering implies mutually inconsistent sequence lengths that are never stated. The direction and monotonicity are not in doubt; the exact rate figures are not reproducible from published numbers, and are tabulated here with that caveat. The article uses its youngest and oldest bins as two rows of its rate table. ↩ ↩
- 33.
Översti, Sanni, and Jukka U. Palo. “Variation in the substitution rates among the human mitochondrial haplogroup U sublineages.” Genome Biology and Evolution 14(7):evac097. 2022. A dating study of 234 ancient human mitochondrial genomes from haplogroup U, a branch of the maternal family tree, each dated by radiocarbon (calibrated years before present, cal BP). Of these, 212 are Holocene and 84 are younger than 4,500 cal BP. The authors ran twelve root-to-tip regressions (tests of how many mutations build up with time) on non-clock trees run in TempEst, all with positive slopes and R² from 0.10 to 0.84. Their haplogroup U4 regression on 42 Holocene-only samples gives 8.80 × 10-8 per site per year, the fastest defensible figure in the literature and the best number available to the young-earth side; even it implies a U4 coalescence beyond 12,000 years, for one small twig of the tree. The article uses it for the largest lineage-to-lineage rate difference in human mitochondrial DNA, about twofold between U2/U5a and U5b. ↩ ↩
- 34.
Kim, Hie Lim, and Stephan C. Schuster. “Poor Man's 1000 Genome Project: Recent Human Population Expansion Confounds the Detection of Disease Alleles in 7,098 Complete Mitochondrial Genomes.” Frontiers in Genetics 4. 2013. A survey of 7,098 complete human mitochondrial genomes that reports how much people's mitochondrial DNA differs on average. The mean pairwise difference across the whole molecule is 39.8 worldwide, with a maximum of 123; within Africa it is 64.0 (maximum 123), among non-Africans 35.4 (maximum 71), in Asians 30.1 and in Europeans 30.6. The control region differs by 9.9 over 1,122 bp and the rest of the molecule by 29.9 over 15,447 bp, a 4.6-fold difference in per-site diversity that the young-earth clock, whose D-loop and whole-genome rates differ by only 1.13-fold, does not reconcile. The paper is number 13 in volume 4. Jeanson's stated non-African targets of 39.8 and 39.3 are in fact the global figures. The article uses it as the observed comparator for every clock calculation in the article. ↩ ↩
- 35.
Mean and maximum pairwise substitutions by region for 54 living humans (the revised Cambridge Reference Sequence plus Ingman's 53 genomes), 36 Neanderthal records representing 35 individuals, and 7 Denisovan records, all on reference coordinates with the length-variable tracts 303–315, 16,184–16,193 and 8,281–8,289 left out. The code reads the pairwise matrices and the sequence inventory and writes the panels. ↩ ↩
- 36.
Mean and maximum pairwise substitutions between groups (living humans, Neanderthals, Denisovans and the one Sima de los Huesos sequence), from the same matrices and code as the within-group panels. Every tip is treated as living today, so these are distances, not split dates; a fossil's age shortens its branch and is handled separately. ↩ ↩
- 37.
Brown, Samantha et al. “The earliest Denisovans and their cultural adaptation.” Nature Ecology & Evolution 6:28–35. 2022. A dating and ancient-DNA study of the oldest Denisovan remains from Denisova Cave. It places the layer holding Denisova 19 at 187,000 to 217,000 years (95.4% interval) and finds that Denisova 19 and 21 carry identical mitochondrial DNA. These are ages of the layer, not direct dates of the individuals; the companion uses them for its one dated-pair check. ↩ ↩