"The theory does not deny the role of natural selection in determining the course of adaptive evolution, but it assumes that only a minute fraction of DNA changes in evolution are adaptive in nature" — Motoo Kimura, The Neutral Theory of Molecular Evolution, 19831
"we have reached this peculiar situation where the analysis of amino acid substitution rates has both led to a clear rejection of a major claim of the neutral theory—that most observable changes are neutral—and at the same time indisputably demonstrated the importance of drift in protein evolution." — Nicolas Galtier, Genome Biology and Evolution, 20242
Summary
The neutral theory of molecular evolution says that most of the DNA differences between species were not chosen by natural selection. They are mutations that made no difference to survival, and they spread through a population by chance, the process called genetic drift. Motoo Kimura proposed it in 1968, and Jack King and Thomas Jukes independently the next year.3 4 This article asks whether it is true, and what follows for evolutionary biology and for the young-earth and design arguments that use it. The formulas behind it are in Population Genetics; Kimura's answer to Haldane's cost of selection is in Haldane's Dilemma.
- The theory is a claim about most changes, not all of them. Kimura's own statement is that "only a minute fraction of DNA changes in evolution are adaptive," and that the theory "does not deny the role of natural selection in determining the course of adaptive evolution."1 Its modern defenders add that most mutations that do matter are harmful, and selection removes them.5
- The founding argument was wrong, and both sides now say so. Kimura's 1968 case was that DNA changes too fast for selection to pay for it. He got that rate by spreading the protein rate over all four billion letters of a mammal genome. Andrew Kern and Matthew Hahn, who oppose the theory, write that this overstates the protein-coding sites "by two orders of magnitude,"6 and Masatoshi Nei, who defended it, made the same correction in 2005.7
- Across the whole human genome, the mainstream's own numbers agree with the theory. About 15 million single-letter changes separate the human line from the human–chimpanzee ancestor.8 The best recent estimate of the adaptive protein changes among them is about 4,000, with a range of 2,900 to 5,100.9 That is about 0.03% of the total, and the most generous figures in the literature raise it to about 0.4%. This site's Functional DNA article says over 99.98% of mammal differences spread by chance. The mainstream count gives about 99.97% at its central estimate and 99.6% at its most generous, so that page's figure is right in substance and slightly high in its last digit.8
- Inside the proteins of species with large populations, the theory's central claim has failed. About half the amino-acid changes between fruit-fly species look adaptive,10 so do more than half in enteric bacteria11 and wild mice,12 and across 44 animal species the median is 0.65.13 Humans sit at 0 to about 25%,14 15 and most plants near zero.16 Nicolas Galtier's 2024 review calls this "a clear rejection of a major claim of the neutral theory," and in the same sentence says drift's importance in protein evolution has been "indisputably demonstrated."2
- The nearly neutral version passes a test it set itself. Tomoko Ohta predicted in 1973 that slightly harmful mutations would fix more often in small populations.17 The adaptive share does track population size, mainly because large populations weed out more of the slightly harmful changes.18 Relative to silent change, protein change runs at about 0.20 on the human and chimpanzee branches against 0.13 on the mouse and rat branches, and the genome authors put the excess down to relaxed constraint, not adaptation.19
- Lewontin's paradox is real and unsolved, and it is not a junk-DNA question. Under strict neutrality, genetic diversity should rise with population size. Across 172 animal taxa, population sizes span about twelve orders of magnitude and diversity about three.20 Selection at linked sites removes some diversity, about 17% in humans,21 but by its defenders' own arithmetic not nearly enough.22 The leading explanation is that large populations are still far from equilibrium after past growth.22 A young earth is the extreme case of that, so the paradox is no evidence against one, nor for one.
- The 2018 challenge is not a junk-DNA argument. Kern and Hahn say outright that the large nonfunctional fraction of eukaryotic genomes is not in dispute.6 Their case rests on proteins and on linked selection, and the reply from seven senior population geneticists defends the theory on the same ground.5
- But the genome-wide form of the theory does lean on junk. The defenders answer the fruit-fly figures by pointing out that they concern "only the small fraction of the genome that codes for proteins."5 Ohta built her theory on the premise that "most genomic DNA of higher organisms can freely accumulate base substitutions."23 Kimura's showcase for the theory's top speed was pseudogenes, which he regarded "as indeed dead."1 Today the premise is measured, not assumed: estimates of the human genome under purifying selection run from 3.5% to about 26%.24 25 26 27 If far more of it had a function that depends on its exact letters, the 15 million human-branch changes could not mostly be neutral, and no evolutionary account could easily explain them, as the mainstream's Dan Graur concedes.28
- Young-earth writers mostly use the theory rather than refute it. John Sanford's genetic entropy argument needs Ohta's slightly harmful mutations to drift to fixation, and Robert Carter calls it an idea that "leverages the beliefs and language of the neutralist position."29 That premise is measured: a quarter to a third of new amino-acid mutations in humans are too weakly harmful for selection to see,30 and Kimura himself wrote that such mutations "accumulate continuously in every species," though he judged the loss tiny and easily offset.31 How much of the genome lies in that zone is the real dispute. Carter dates Adam and Eve with the theory's rate result,32 and Walter ReMine concedes that "the vast majority of substitutions would be neutral."33 Their claim that the theory requires junk DNA is right about its genome-wide form and wrong as a statement of the theory, which allows function without sequence-specific constraint.1 On a young earth, neutral substitution at the measured rate changes about three sites in a million per lineage,8 so differences between created kinds are not fixed mutations at all; within kinds, drift sorts created variety fast enough only in small populations.
Verdict. Most DNA differences between mammal species are neutral or nearly so, by the mainstream's own count, and that conclusion rests mainly on the measured finding that most of the genome is not under constraint. Inside the proteins of species with large populations the original theory has failed. For the design argument, the theory is a real concession that selection did not build most differences, but it becomes evidence against evolution only if a large, sequence-specific function is shown in the DNA that differs, and the data within living humans do not yet show one. For the genetic-entropy argument the theory supplies the premise rather than the refutation: Kimura conceded that effectively neutral mutations are slightly harmful and accumulate, and the open question is how much of the genome they touch.
Part 1: What the neutral theory claims
Kimura's 1968 argument
Motoo Kimura was a Japanese population geneticist who had spent the 1950s and 1960s working out the mathematics of chance in small populations. In February 1968 he published a three-page letter in Nature whose subtitle states the whole idea: "Calculating the rate of evolution in terms of nucleotide substitutions seems to give a value so high that many of the mutations involved must be neutral ones."3
His reasoning ran in three steps. First, he averaged the rates at which three proteins had changed across mammals: haemoglobin, cytochrome c and triosephosphate dehydrogenase. Second, he spread that rate over the whole genome, which he took to be four billion DNA letters, and found that a mammal lineage swaps one letter about every two years. Third, he compared that with J. B. S. Haldane's 1957 estimate that selection can afford to replace one gene about every 300 generations. The gap was several hundredfold, so he concluded that the replacements could not all have been paid for by selection: "the very high rate of nucleotide substitution which I have calculated can only be reconciled with the limit set by the substitutional load by assuming that most mutations produced by nucleotide replacement are almost neutral in natural selection."3
The same letter contains the theory's most useful result. For a mutation that makes no difference to fitness, the chance that it eventually spreads to the whole population equals its starting frequency. The number of new neutral mutations each generation is proportional to population size, and each one's chance of spreading is inversely proportional to it, so population size cancels out. In Kimura's words, "new alleles may be produced at the same rate per individual as they are substituted in the population in evolution."3 The rate at which neutral changes accumulate between species equals the mutation rate. That is the basis of the molecular clock, which has its own article. The derivation is in Population Genetics.
King and Jukes, 1969
A year later two American biochemists, Jack Lester King and Thomas Jukes, reached the same conclusion from the chemistry of proteins and published it in Science as "Non-Darwinian Evolution." Their summary line is the one most quoted: "Natural selection is the editor, rather than the composer, of the genetic message."4 They made two observations that matter for this article. At the cytochrome c gene, "at least 90 percent of all substitutional mutations ... are harmful and are rejected by natural selection,"4 so the theory never said that most mutations are harmless, only that most of the ones that survive to become differences are. And they argued that "not much more than 1 percent of mammalian DNA codes for proteins,"4 so most of the genome had no protein to protect.
The theory in its mature form
Kimura's 1983 book, The Neutral Theory of Molecular Evolution, states the claim with its limits:
The neutral theory asserts that the great majority of evolutionary changes at the molecular level, as revealed by comparative studies of protein and DNA sequences, are caused not by Darwinian selection but by random drift of selectively neutral or nearly neutral mutants. The theory does not deny the role of natural selection in determining the course of adaptive evolution, but it assumes that only a minute fraction of DNA changes in evolution are adaptive in nature, while the great majority of phenotypically silent molecular substitutions exert no significant influence on survival and reproduction and drift randomly through the species.1
Two things follow from that paragraph. The theory is about molecules, not about eyes or wings; Kimura left the shaping of bodies to selection. And "neutral" does not mean "unconstrained." Kimura wrote that "the existence of selective constraint itself does not contradict the neutral theory."1 A gene can be under heavy constraint, with most of its possible mutations rejected, and still have most of its accepted changes be neutral ones.
The book's main evidence was a pattern Kimura stated as a rule: "Functionally less important molecules or parts of molecules evolve (in terms of mutant substitutions) faster than more important ones."1 The surface of haemoglobin changes faster than its oxygen pocket. The middle piece of proinsulin, which is cut away before insulin works, changes several times faster than insulin itself. Synonymous sites, where a letter change does not alter the amino acid, change faster still. And the fastest of all were pseudogenes, broken copies of genes, which Kimura regarded "as indeed dead."1 Kimura argued that a selection-driven view gives no reason to expect this order, while the neutral view requires it: the less a part matters, the larger the share of its mutations that do no harm, and so the more of them drift through.1 The data matched that order.
Ohta's nearly neutral theory
Tomoko Ohta, Kimura's colleague at the National Institute of Genetics in Japan, extended the theory in 1973 to mutations that are slightly harmful: "The concept of neutral mutant substitution in the population by random genetic drift can be extended to include random fixation of very slightly deleterious mutations which have more chance of being selected against than of being selected for. If this class of mutant substitution is important, we can predict that the evolution is rapid in small populations or at the time of speciation."17
The reason is that selection's grip depends on population size. A mutation that lowers fitness by one part in a hundred thousand is weeded out efficiently in a population of a million, but in a population of ten thousand it behaves almost as if it were neutral and can drift to fixation. The threshold, the "drift barrier," is roughly one over twice the effective population size; Population Genetics explains why it is a soft shoulder rather than a wall. Ohta's 1992 review put her thesis this way: "I propose that a substantial fraction of mutant substitutions at the molecular level are caused by random fixation of very slightly deleterious mutations."23 She divided the genome in two. "The fundamental idea is that most genomic DNA of higher organisms can freely accumulate base substitutions, i.e. most new mutations are selectively neutral. On the other hand, amino acid substitutions are more likely to be influenced by natural selection, i.e. many of them may be regarded as nearly neutral, or very slightly deleterious."23 That two-tier picture comes back in Part 5.
The modern statement
The fullest recent statement is a 2019 reply by seven senior population geneticists, Jeffrey Jensen, Bret Payseur, Wolfgang Stephan, Charles Aquadro, Michael Lynch, Deborah Charlesworth and Brian Charlesworth. They define the theory as the claim that "most de novo mutations are either sufficiently deleterious in their effects on fitness that they have little chance of becoming fixed in the population, or are under such weak selection that they may become fixed as a result of genetic drift."5 Their five-point "modern version" opens with two claims: "A large fraction of the genome of organisms studied to date is subject to mutations that are effectively neutral with respect to their fitness effects, and hence evolve under genetic drift,"5 and "The great majority of newly arising mutations that do affect fitness (i.e., non-neutral mutations) are deleterious, and the predominant mode of natural selection is purifying in nature, removing these deleterious mutations from populations."5
So the theory makes five checkable predictions:
- Neutral changes accumulate between species at the mutation rate, whatever the population size.
- Parts of the genome under less constraint change faster.
- Within a species, neutral genetic diversity should scale with population size times mutation rate.
- Polymorphism within species and divergence between species are two stages of one process, so the two should be in proportion.
- On the nearly neutral version, slightly harmful changes accumulate faster in species with small populations.
Prediction 1 is the molecular clock, and it holds only roughly. If neutral substitutions arrived at random, like the clicks of a Geiger counter, the variance in their number across lineages would equal the mean, an index of dispersion of 1. Trevor Bedford, Ilan Wapinski and Daniel Hartl measured mean indexes of 2.1 across 3,081 yeast proteins, 4.2 across 7,174 fruit-fly proteins and 6.5 across 8,065 mammal proteins.34 The clock is least regular where populations are smallest, and the authors add that this "does not necessarily imply positive selection or adaptive evolution; overdispersion is compatible with nearly neutral evolution occurring in a heterogeneous fashion."34 Rates measured over a few generations also run faster than rates measured over millions of years. One of the mainstream's proposed causes is the slightly harmful class of mutations described in Part 3, counted in families before selection removes them. Molecular Clocks works through both departures, and through the young-earth clocks of Robert Carter and Nathaniel Jeanson that rest on the short-term rate. Prediction 4 is what the McDonald–Kreitman test checks, in Part 3. Prediction 3 is Lewontin's paradox, in Part 4.
The founding argument has not survived, on either side's account
Kimura's 1968 calculation needed every one of four billion letters to be a possible target of selection. Kern and Hahn, the theory's leading critics, write that "Kimura overestimates the number of protein-coding sites in the genome by two orders of magnitude (he uses 4×10⁹ bp, whereas the real number is closer to 3×10⁷ bp in the largest genomes)," and that fixing this number alone removes the conflict with Haldane.6 Masatoshi Nei, who defended both Haldane's cost and the neutral theory, made the same correction in 2005: counted per gene, the substitution rate falls far below Haldane's limit, and counted per amino acid it is "about two times higher," not hundreds of times.7 Haldane's Dilemma has the full history. The theory's second original pillar, that proteins change at a steady rate in every lineage, Kern and Hahn call "an illusion" once more proteins were sequenced.6 Whether the theory stands today depends on the later evidence, which is the subject of the rest of this article.
Part 2: The 2018 challenge and the reply
Kern and Hahn: "overwhelmingly rejected"
For the theory's fiftieth anniversary, Molecular Biology and Evolution asked Andrew Kern of Rutgers University and Matthew Hahn of Indiana University how it had fared. Their answer, in the paper's second sentence, was "In a word, poorly."6 Their abstract: "We argue that the neutral theory was supported by unreliable theoretical and empirical evidence from the beginning, and that in light of modern, genome-scale data, we can firmly reject its universality."6
They began by stating what they were not disputing. The neutral theory, they wrote, "is not simply a statement about the presence of neutral mutations, nor about the large fraction of eukaryotic genomes that are nonfunctional—neither of these assertions would be contested by current competing hypotheses."6 What they rejected was a stronger thesis: "that differences between species are due to neutral substitutions (not adaptive evolution), and that polymorphisms within species are not only neutral but also have dynamics dominated by mutation-drift equilibrium."6
Their case has four parts.
- The original evidence failed. Kimura's cost argument used the wrong number of sites, and steady protein rates turned out to be uneven, as Part 1 describes.
- Many amino-acid changes look adaptive. McDonald–Kreitman studies in two fruit-fly species found "∼50% of all amino acid substitutions have been fixed by positive selection," and similar results followed in other species.6 They concede the human exception: "In humans, though the overall fraction of amino acid substitutions fixed by positive selection is estimated to be zero," while genes that deal with pathogens show strong signs of adaptation.6
- Selection reaches almost everywhere through linkage. Across 40 genomes, regions where chromosomes recombine less carry less diversity. Under strict neutrality recombination should not matter, so they read the pattern as the footprint of selection on nearby sites: "these results imply that almost no loci are free from the effects of selection, in any organism."6 Because "An equilibrium between mutation and drift is a central tenet of the neutral theory," they conclude that "current data appear to be fundamentally incompatible with the neutral theory."6
- Adaptation is not limited by mutation. It draws on standing variation, on genes borrowed from other species, and on many genes of small effect at once, none of which leaves the classic signature of a single sweep.
Their verdict: "we believe that the neutral theory has been overwhelmingly rejected."6
Jensen and six colleagues: "abundantly clear" that it is correct
The reply came the following January, in Evolution, from seven population geneticists, among them Michael Lynch and Brian and Deborah Charlesworth, two of the three authors who first described background selection in 1993.35 They called Kern and Hahn's position "probably a minority view" and wrote that "it is now abundantly clear that the foundational ideas presented five decades ago by Kimura and Ohta are indeed correct."5
Their answers, point by point:
- The fruit-fly numbers concern a small part of the genome. "These inferred frequencies of adaptive substitutions mostly concern only the small fraction of the genome that codes for proteins (e.g., <2% of the human genome; see Lander et al. 2001)."5 Focusing on fast-evolving genes and then rejecting neutrality is, they write, circular.5
- Most linked selection is purifying, not adaptive. "In an explicit comparison between models of widespread purifying selection on weakly deleterious alleles versus recurrent positive selection on beneficial alleles, Lohmueller et al. (2011) found a much better fit of the former to the observed pattern in humans."5 Background selection, the loss of neutral variation that happens to sit beside harmful mutations being weeded out, explains the recombination pattern without adaptation.
- Linkage does not change the rate of neutral change between species. Neither background selection nor sweeps "affect the probability of fixation of neutral mutations ..., which determines the rate of neutral sequence evolution."5 Both lower the effective population size of the affected region, which makes more mutations behave as neutral, not fewer.
- Without a neutral baseline, stories multiply. "In the absence of an appropriate null model accounting for these processes that are common to the genome as a whole, inappropriate adaptive story-telling will be likely to proliferate."5
Where the dispute stands
The two papers mostly disagree about what the theory says. Kern and Hahn test the strong form: between-species differences are neutral, and diversity sits at mutation–drift equilibrium. Jensen and colleagues defend a weaker form: much of the genome evolves by drift, and most selection removes harmful mutations. Both can be true at once. The fruit-fly protein figures falsify the strong form for fruit-fly proteins. The noncoding bulk of a mammal genome fits the weaker form.
Nicolas Galtier's 2024 review of fifty years of the debate gives the cleanest verdict from inside the field. The theory's two founding observations, abundant polymorphism and roughly steady protein evolution, "are still very much alive and important aspects of current molecular evolutionary thinking."2 But on proteins, the field has "both led to a clear rejection of a major claim of the neutral theory—that most observable changes are neutral—and at the same time indisputably demonstrated the importance of drift in protein evolution."2 His diagnosis of the founders: "Enthused by their discovery, they slightly overstated the prevalence of neutrality."2 The question is now one of proportion: "The issue, however, is no longer phrased in terms of corroborating versus challenging the neutral theory, but rather in terms of quantifying the relative contribution of neutral and selective processes to the observed patterns."2
None of the three papers argues that most of the genome is functional. That matters for how creationists have read the 2018 exchange, and it is taken up in Part 5.
Part 3: How much of the change is adaptive? The measurements
The McDonald–Kreitman test
In 1991 John McDonald and Martin Kreitman found a way to count adaptive changes directly. Their reasoning: "If the observed substitutions are neutral, the ratio of replacement to synonymous fixed differences between species should be the same as the ratio of replacement to synonymous polymorphisms within species."36 Silent changes serve as the neutral yardstick. If protein-changing differences between species outnumber what the variation within species predicts, the excess is the share that selection drove to fixation. That share is called alpha (α). At the Adh gene of fruit flies they found "more fixed replacement differences between species than expected."36
The test has known biases, and they matter for the numbers below. Slightly harmful protein changes linger as rare variants inside a species but seldom fix, so they inflate the within-species count and push α down, sometimes below zero. Population growth or shrinkage pushes it either way. Modern methods filter out rare variants or model the distribution of fitness effects (taken from Beneficial and Harmful Mutations as an input here), and different methods give different answers on the same data.15
The numbers by taxon
| Group | Adaptive share of substitutions (α) | Source |
|---|---|---|
| Fruit flies, D. simulans and D. yakuba | 45% of amino-acid changes | Smith and Eyre-Walker 200210 |
| Fruit flies | about 50% of amino-acid changes, about 20% of intron changes | Eyre-Walker and Keightley 200937 |
| Fruit flies, noncoding | about 20% intronic and intergenic, 60% UTRs | Andolfatto 200538 |
| 44 animal species pairs | median 0.65, range 0.18 to 0.93 | Galtier 201618 13 |
| Wild mice, Mus musculus castaneus | 57% (0.30 to 0.76) | Halligan et al. 201012 |
| House mice by subspecies | about 40% in castaneus, 13% and 12% in the two smaller subspecies | Phifer-Rixey et al. 201239 |
| Mice, noncoding | 19% in UTRs and conserved noncoding elements | Halligan et al. 201340 |
| Nine plant species pairs | mean −0.18 (SE 0.47), little evidence except sunflowers | Gossmann et al. 201016 |
| Enteric bacteria, E. coli and Salmonella | more than 50% (56%, lower bound 45%); 34% with Salmonella data | Charlesworth and Eyre-Walker 200611 |
| Humans | 10–20% | Boyko et al. 200814 |
| Humans | 13.5% (9.6% to 17%) | Uricchio, Petrov and Enard 20199 |
| Humans | at least 20% under their preferred model; −0.41 by the plain test | Zhen et al. 202115 15 |
Estimates across species run "from nearly zero ... up to roughly 90%."2 The plain conclusion is that in most animals with large populations, about half or more of the protein changes between related species were driven by selection. Galtier, after measuring 44 species pairs, wrote that "Observing a high α therefore seems to be the rule, and a low α the exception, making humans and apes a peculiar, intriguing case."18 For proteins in those species, the neutral theory's founding claim is false.
Humans are the exception, and the count is in the thousands
Human estimates run from zero to about a quarter, depending on method. Kern and Hahn summarise the human figure as "estimated to be zero," citing Adam Boyko and colleagues among others,6 but Boyko's own abstract gives 10–20%.14 Ying Zhen and colleagues show why the methods disagree: the plain test gives −0.41 for human against chimpanzee, and the figure rises to 0.25 once the sizes of the outgroup and ancestral populations are modelled.15
The best recent count is from Lawrence Uricchio, Dmitri Petrov and David Enard. They counted 29,925 protein-changing differences on the human branch since the split from chimpanzees and estimated α at 0.135.9 That is about 4,040 adaptive protein changes, 2,900 to 5,100 on their interval. Out of about 15 million single-letter changes on the human branch, it is 0.027%; at Boyko's top figure of 20% it is still only 0.04%.8 Haldane's Dilemma compares these counts with the reproductive budget for paying for them.
So for humans, both halves of the theory hold: most protein changes look neutral or nearly so, and the genome as a whole is overwhelmingly so. The Haldane's Dilemma article found that the popular replies to young-earth critics concede fewer adaptive changes than this literature now estimates, so the literature's count is the one to use, on either side.
Population size is the pattern, and it is the nearly neutral prediction
The adaptive share is not random across species. It tracks population size. Galtier found α "positively correlated to species effective population size," but traced most of the correlation to the other side of the ratio: it "appears to be primarily driven by a decreased rate of nearly-neutral amino-acid substitution because of more efficient purifying selection in large populations."18 Within one species of mouse, the subspecies with the largest population has α near 40%, the others about 12% and 13%, and "the higher estimate of a for M. m. castaneus appears to reflect not only more adaptive fixations but also more effective purifying selection."39 Plants mostly show little adaptation, and Gossmann and colleagues suggest it "may be because many plant species have modest effective population sizes."16
This is what Ohta predicted in 1973. Slightly harmful changes fix more readily in small populations, so small populations have more non-adaptive protein changes and a lower α. The chimpanzee genome paper shows the same effect directly. The protein-changing rate relative to the silent rate is about 0.20 on the human and chimpanzee branches against about 0.13 on the mouse and rat branches, so about 35% of the hominid amino-acid changes are in excess of the rodent pattern. The authors conclude that relaxed constraint, not extra adaptation, "seems to be the principal explanation."19 That is why Galtier could write that the same data rejected a major claim of the neutral theory and "indisputably demonstrated the importance of drift."2 The strict theory failed; the nearly neutral theory passed a test it set in advance.
The slightly harmful class, measured
The nearly neutral theory needs a large class of mutations whose harm is too small for selection to see. That class has been measured, and it is real. Adam Eyre-Walker and Peter Keightley's 2007 review starts from a definition: "it seems unlikely that any mutation is truly neutral in the sense that it has no effect on fitness." A mutation counts as "effectively neutral" when the product of the effective population size and its selection coefficient, Nes, "is much less than 1," so that its fate "is largely determined by random genetic drift."41 For hominids, with effective sizes of 10,000 to 30,000, they estimated that fewer than 30% of amino-acid-changing mutations are effectively neutral, which "implies that at least 70% are deleterious."41 Fits to human variation agree. Adam Boyko and colleagues put 27–29% of new amino-acid mutations below an effect of 0.01% and 30–42% between 0.01% and 1%.14 Bernard Kim, Christian Huber and Kirk Lohmueller, with samples of over 1,000 Europeans, found 34–37% below 0.01% and about a quarter below 0.001%, and concluded that "nearly neutral forces play a larger role in human evolution than previously thought."30 Beneficial and Harmful Mutations compares these fits and their limits: each depends on a model of population history, and most have no beneficial class, so slightly helpful mutations are counted as neutral.
Kimura saw what this class implies. In the 1979 paper that set out his distribution of mutation effects, he wrote: "Under the present model, effectively neutral, but, in fact, very slightly deleterious mutants accumulate continuously in every species." He put the loss at about one part in ten million of fitness per generation, from 10,000 protein-coding genes, and called the threat to a species "a moot point," because "this will easily be taken care of by adaptive gene substitutions that must occur from time to time (say once every few hundred generations)."31 In his 1983 book he called the accumulation "essentially equivalent to the deterioration of the environment," which adaptive substitutions must answer "to save the species from extinction."1 Alexey Kondrashov extended the arithmetic to the whole genome in 1995. Mutations whose harm lies between about one over the genome size and one over four times the effective size are "neutral enough to accumulate almost freely, but are still deleterious enough to make an impact at the level of the whole genome." If a tenth of all sites were in that range, an average person would carry about 100 lethal equivalents. Unlike Kimura, Kondrashov wrote that "This paradox cannot be resolved by invoking beneficial mutations or environmental fluctuations," and he weighed other ways out, including soft selection and interactions that make many small harms count for more together than apart.42
This is the premise of the young-earth genetic-entropy argument, stated by the founder of the theory and by a leading mainstream theorist in their own words. It does not settle the size of the effect. Kimura's own figure of one part in ten million per generation compounds to about 0.002% over 240 generations and about 2% over 200,000.43 The dispute is over how many sites outside the protein-coding genes sit in the zone, and the constraint measurements of Part 5 cannot answer it. A site whose harm sits at the upper edge of Kondrashov's range still changes at 58% of the neutral rate, and one ten times weaker at 95%, so a conservation scan, which looks for sites that change far more slowly than neutral ones, sees little of the range.43 Christian Huber, Bernard Kim and Kirk Lohmueller tested a standard conservation score by simulation in 2020. Once sites gain and lose function over evolutionary time, as most human noncoding constraint seems to, "There is no improvement in statistical power with increasing tree size for very weak selection," and "Commonly used GERP score thresholds miss over half of the noncoding sites in the human genome where mutations are under purifying selection."44 44 Genetic Entropy prices that question and the mainstream's answers to it.
dN/dS: most protein change is rejected
The ratio of protein-changing to silent differences, written dN/dS or Ka/Ks, measures how much of the mutational supply survives. For human against chimpanzee it is 0.23, which the genome authors read as "77% of amino acid alterations in hominid genes are sufficiently deleterious as to be eliminated by natural selection."19 Only 585 of 13,454 genes, 4.4%, showed a ratio above 1, the classic sign of adaptation, and the authors' simulations expected at least 263 such genes "simply by chance."19 Between human and mouse the median ratio across 12,845 genes is 0.115.45 A scan of about 16,500 human genes against five other mammals found significant evidence of positive selection in 400, with the authors warning that the counts reflect "differences in power more than differences in the prevalence of selection."46 This is King and Jukes's pattern at genome scale: most changes to working proteins are thrown out, and most of the ones kept look neutral or slightly harmful.
Noncoding DNA
Adaptive change is not confined to proteins. In fruit flies, Peter Andolfatto estimated that about 20% of intron and intergenic divergence and 60% of UTR divergence was driven by positive selection, and concluded "that a large fraction of the non-translated genome is functionally important."38 Fruit-fly genomes are compact, with far less noncoding DNA per gene than ours. In wild mice, α in conserved noncoding elements is lower than in proteins, 0.19 against 0.32–0.38, but those elements are so much larger a target that they hold at least four times as many adaptive changes as the proteins do.40 In humans, Eyre-Walker and Keightley found "little evidence for adaptive substitutions" in protein-coding and noncoding data alike, while warning that the true protein figure "could be as high as 40%."37 Noncoding adaptation is real in species with large populations and compact genomes, and poorly measured in ours.
Part 4: Lewontin's paradox and selection at linked sites
The paradox
If most variation within a species is neutral, the amount of it should be set by a balance between new mutations and their loss by drift. For neutral DNA the expected diversity is about four times the effective population size times the mutation rate. Mutation rates differ only about tenfold across animals,20 while population sizes differ by many orders of magnitude. So an insect species numbering in the billions should be vastly more diverse than an elephant species.
Richard Lewontin, a Harvard geneticist whose own protein surveys in the 1960s had helped show how variable natural populations are, pointed out in 1974 that this does not happen. The heterozygosities then known lay "between 0.056 and 0.184," which on the neutral formula meant "we are required to believe that higher organisms including man, mouse, Drosophila and the horseshoe crab all have population sizes within a factor of 4 of each other."47 His conclusion: "Despite the unchallengeable agreement between the heterozygosity observed in any given case and the predictions of neutral theory, the application of neutral theory to every case of heterozygosity turns out to lead to an absurd result."47 This is Lewontin's paradox of variation.
DNA sequencing has widened the range without solving the problem. Across 167 species, Ellen Leffler and colleagues found diversity "ranges from 0.01% per base pair in Lynx lynx to 8.01% in Ciona savignyi, a span of only 800-fold."48 Vince Buffalo's 2021 compilation of 172 animal taxa puts diversity across about three orders of magnitude against population sizes that vary over about twelve.20
The paradox began as an argument against the neutral theory. John Maynard Smith and John Haigh, who first described genetic hitchhiking in 1974, wrote that "the uniformity of H between species is powerful evidence against the view that the observed polymorphisms are in the main selectively neutral," and called their own hitchhiking model "a last ditch attempt to save the neutral mutation theory."49
First explanation: selection at linked sites
When selection acts on one site, it changes the fate of neutral variants sitting next to it on the chromosome. A favourable mutation that sweeps to fixation drags its neighbours with it and wipes out their variation; this is hitchhiking. A harmful mutation that selection removes takes its neighbours with it too; this is background selection. Both effects are strongest where recombination is rare, so they predict lower diversity in low-recombination regions. That correlation is found in fruit flies50 and in most of the 40 genomes Russell Corbett-Detig and colleagues examined, where the effect was larger in species with larger populations. They wrote that "natural selection may provide an explanation for this longstanding paradox of population genetics."51
How much diversity does linked selection remove?
| Species | Estimate of neutral diversity removed | Source |
|---|---|---|
| Humans | about 6% (B ≈ 0.94) | Charlesworth 201252 |
| Humans | 19–26% on the autosomes | McVicker et al. 200953 |
| Humans | about 17% on average, sweeps adding "little or no" effect | Murphy et al. 202221 |
| Drosophila melanogaster | median B 0.56 on autosomes, so about 44% | Comeron 201454 |
| Drosophila melanogaster | 77–89% genome-wide, a model fit | Elyashiv et al. 201655 |
The largest figure has a weakness its authors state: the fitted harmful-mutation rates behind it "far exceed direct estimates of the total mutation rate (by 4-9-fold in exons and UTRs".55 And even taken at face value it does not solve the paradox. Brian Charlesworth and Jeffrey Jensen, defenders of the neutral theory, worked through it in 2022: an 80% reduction means the neutral coalescence time is "about 5-fold larger than the value indicated by current diversity levels; this is far smaller than the factor of 100 or more suggested by the probable current population size of D. melanogaster ..., and so is insufficient to resolve LP for this species."22 Buffalo's verdict across all 172 taxa is the same: "even assuming strong and abundant selection, models of linked selection are unlikely to explain the observed relationship between diversity and census sizes across species."20
In humans, linked selection reaches almost every variant. Florian Pouyet and colleagues found that background selection and a recombination-linked repair bias "together affect as much as 95% of the variants of our genome."56 That is linkage to a small constrained fraction, not a sign that 95% of sites are themselves functional.
Second explanation: populations far from equilibrium
The neutral formula gives the diversity a population reaches once mutation and drift have balanced. Reaching that balance takes about four times the effective size in generations, which for a large population is an enormous time. Charlesworth and Jensen calculate that for a population that grew to an effective size of a billion, "the time for π to approach ¾ would be of the order of 10⁹ generations."22 Populations that expanded long ago may still be far below their equilibrium diversity, and "This scenario thus seems an excellent candidate for resolving LP."22 Buffalo reached the same place from the other side: since linked selection fails, "this supports the hypothesis the diversity across species are shaped primarily by past demographic fluctuations."20 Life history fits this reading. Across 76 animal species, Jonathan Romiguier and colleagues found diversity "accurately predicted by key species traits related to parental investment," not by geographic range.57
What the paradox does to the neutral theory
It falsifies the strict version Kern and Hahn attacked, in which diversity sits at mutation–drift equilibrium. Nobody defends that version now, and Galtier notes that the demographic explanations are "compatible with the neutral theory."2 It does not touch the theory's claim about divergence between species, because linked selection does not change how often a neutral mutation fixes.5 The irony is that a datum first raised against neutralism is now mostly explained by neutral processes plus a modest share of selection, most of it purifying.
What it means for a young earth
A young earth is the most extreme non-equilibrium history there is. At the measured human rate, reaching mutation–drift balance takes 40,000 to 80,000 generations, 1 to 2.3 million years;8 for a large insect population it takes far more generations. On a young earth no species is anywhere near equilibrium, so diversity would not be expected to track population size. The paradox is therefore not evidence against a young earth, and the mainstream's leading explanation for it is the same kind of explanation a young-earth model would give.
It is not evidence for a young earth either. An expansion five million years ago does the same job in Charlesworth and Jensen's arithmetic.22 And a young-earth model still has to say where the existing diversity came from. Mutation at the measured rate over about 240 human generations adds only about three new differences per million sites per lineage,8 while living humans differ at about one site in a thousand, the figure behind the familiar human effective size of 10,000 to 20,000 explained in Population Genetics. The young-earth answer is created diversity, which is argued in Adam, Eve, and Human Genetic Diversity.
Part 5: Does the theory depend on junk DNA?
What the 2018 critics did not say
It is easy to read the 2018 exchange as "the neutral theory failed because the genome turned out to be functional." That is not what Kern and Hahn argued. They set junk DNA aside at the start: the theory "is not simply a statement ... about the large fraction of eukaryotic genomes that are nonfunctional—neither of these assertions would be contested by current competing hypotheses."6 The strongest critics of the theory accept that most of a large genome does nothing. Their case is about adaptive change inside proteins and about selection reaching neutral sites through linkage. Nobody in the professional dispute argued that the genome is mostly functional, and an article that cites Kern and Hahn for that is citing them for the opposite of what they wrote.
Where the theory leans on a nonfunctional genome
Three of the theory's supports do lean on most of the genome being free of constraint.
- The genome-wide claim. Faced with fruit-fly proteins where half the changes look adaptive, the theory's defenders answered that those figures "mostly concern only the small fraction of the genome that codes for proteins."5 The claim that most genomic differences are neutral is carried by the 98% that is not protein-coding. If that 98% were under selection, the answer would fail.
- Ohta's explanation of the generation-time effect. DNA divergence in general runs faster in short-lived animals, while protein divergence does not. Ohta explained this with two tiers: "most genomic DNA of higher organisms can freely accumulate base substitutions," while amino-acid changes are often slightly harmful.23 The first tier is a no-constraint premise for most of the genome.
- Kimura's top speed. The fastest-changing DNA Kimura knew was in pseudogenes, where "base substitutions at the first and the second codon positions in the pseudogenes occur just as frequently as those at the third position."1 That is what a dead gene should look like, and he used it as the theory's best case. But he answered the suggestion that pseudogenes might have "an unknown but important function" by opinion: "In my opinion, it is much more appropriate to regard pseudogenes as indeed dead."1 This site's Pseudogenes article finds that a few dozen of about 14,700 human pseudogenes have a measured function, and most have never been tested one by one.
Where it does not
The protein evidence does not lean on junk at all. The human-chimpanzee ratio of protein-changing to silent changes is 0.23, which the chimpanzee genome authors read as "77% of amino acid alterations in hominid genes are sufficiently deleterious as to be eliminated by natural selection."19 The McDonald–Kreitman estimates in Part 3, the hominid–murid contrast and the pattern that less important parts of a protein change faster are all measured inside genes whose function no one disputes. Whatever the noncoding genome turns out to do, these results stand.
How much is measured to be under constraint
The mainstream estimates of how much of the human genome is held in place by purifying selection are these:
| Estimate | Share of the human genome | Method | Source |
|---|---|---|---|
| Sullivan et al. 2023 | 3.5% of single bases | per-base conservation, 240 mammals | 27 |
| Huber et al. 2020 | 4.51% of the noncoding genome | population-genetic model of conservation scores, with turnover | 44 |
| Rands et al. 2014 | 8.2% (7.1–9.2%) | indel deficit, with turnover | 24 |
| Zoonomia 2023 | 10.7%, "a minimum" | conservation per element, 240 mammals | 26 |
| Ponting and Hardison 2011 | 10–15% | review, with turnover | 58 |
| Kellis et al. 2014 | 12–15%, plus 4–11% human-specific | union of methods | 25 |
Every one of these is far below a majority. The cross-species estimates assume common descent and treat ancient transposon copies as the neutral yardstick, so a design proponent can fairly say they measure shared conservation, not function. This site's Functional DNA article makes that objection. But the human-only part of the Kellis estimate, 4–11%, rests on variation within living humans and assumes no common ancestry.
Selection does touch almost the whole genome indirectly. Florian Pouyet and colleagues found that background selection and a recombination-linked repair bias "together affect as much as 95% of the variants of our genome."56 That is the reach of linkage to a small constrained fraction, not a finding that 95% of sites are functional.
What changes if much more of the genome is functional
There are two very different ways the genome could be "more functional."
Function that does not depend on the exact letters. Dan Graur, one of the most outspoken defenders of junk DNA, allows a category he calls "indifferent DNA": "DNA sites that are functional, but show no evidence of selection against point mutations," such as "spacers and flanking elements whose presence is required but whose sequence is not im[portant]."28 If most of the extra function is of this kind, the neutral theory is untouched. A letter change in a spacer is neutral even though the spacer does a job. Much of the biochemical activity the ENCODE project reported could fall here, and nothing in the neutral theory forbids it.
Function that depends on the exact letters. This is the claim that matters. The Functional DNA article estimates that at least 30% of human DNA is sequence-specific. This site's GWAS article reran that arithmetic and got 2% to 30% depending on inputs, and 6% to 10% with heritability inputs; but take 30% as the test case. If the 15 million human-line changes fell at random, about 4.5 million of them landed in DNA whose letters matter.8 They cannot mostly be neutral, by the definition of "sequence-specific." Each one must then be one of three things.
- Adaptive. Then selection paid for millions of substitutions. Haldane's Dilemma puts the budget at hundreds to a few thousand on the human line. No one on either side proposes millions.
- Slightly harmful, carried by drift. Then each had an effect below the drift barrier, about 2.5 × 10⁻⁵ for a human effective size of 20,000.8 Function that weak is "sequence-specific" in only the faintest sense, and millions of such fixations are a fitness decline, which is the Genetic Entropy argument. This is close to John Sanford's own picture: "While most sites are probably nearly-neutral, very few, if any, should be absolutely neutral."59
- Never fixed by any population process. This is the young-earth answer: the differences between human and chimpanzee were created.
Or the changes did not fall at random, and the 30% shows a deficit of differences, which is exactly what a conservation measurement looks for. The measurements find a deficit in about 3.5% to 15% of the genome across mammals, and 4% to 11% more within humans. So the data do not show a sequence-specific 30% that evolution would have to explain. These scans see function strong enough for selection to slow change. They cannot rule out the very slightly harmful class of item 2, which changes at more than half the neutral rate (Part 3). If such a fraction were shown, the professional literature agrees on what would follow. Graur wrote that "there exists a misconception among functional genomicists that the evolutionary process can produce a genome that is mostly functional."28
A larger sequence-specific fraction would also raise the number of harmful new mutations per newborn. At the pedigree mutation rate a child carries about 77 new point mutations, so 3.5% sequence-specific DNA means about 2.7 harmful ones a child, and 30% means about 23.8 That is the genetic-load argument, which Genetic Entropy prices against the mainstream's answers, and it is not re-argued here. And more constrained DNA would mean more background selection, so lower neutral diversity everywhere: that pushes the data toward Kern and Hahn's side of the dispute, not toward a functional genome that is also neutral.
The junction, argued both ways
For the design side. The junction is real, and the mainstream says so. The genome-wide form of the neutral theory is not independent evidence against widespread function, because two of its supports, Ohta's free-accumulating bulk and Kimura's dead pseudogenes, assumed the answer. If sequence-specific function were shown in a large part of the genome, the neutral account of the 15 million human-line differences would fail, and so would every other evolutionary account of them, as Graur concedes.
For the mainstream. The no-constraint premise is no longer only an assumption. It is measured, by conservation across mammals and by variation within humans, and the measurements agree on a minority. Pseudogenes passed the neutral theory's own test: all three codon positions change at the same rate, as they should in dead DNA, while in working genes the third position changes fastest. The theory's sharpest critics concede the nonfunctional bulk. And the theory is not the reason the bulk is thought nonfunctional; it is one of several lines that agree.
What would settle it. A sequence-specific function shows up as missing variation in living humans, whatever one believes about ancestry. The test a design proponent needs to pass is a within-human depletion of variation across far more than the 4–11% now measured, in DNA outside known genes and regulatory elements. The large population sequencing projects now under way can run that test.
Part 6: What follows for both sides
What young-earth and design writers have said
Young-earth writers have said a good deal about the neutral theory, and most of it uses the theory rather than attacks it.
- They have Kimura's motive right. Robert Carter of Creation Ministries International wrote that "Kimura saw Haldane's Dilemma as a serious problem and listed it as his main reason for proposing the Neutral Theory of evolution."60 John Sanford and colleagues: "Kimura agreed with Haldane's assessment of the problem, and this motivated him to develop his neutral model of evolution."61 Walter ReMine: "These arguments were never 'for' the neutral theory. They are arguments against selective evolution, severely limiting its speed."62 Kimura's 1968 paper bears this out, as Part 1 shows. What they leave out is that the cost argument rested on the four-billion-site denominator both sides now reject.
- They link it to junk DNA. Carter in 2009: "Junk DNA is a necessary mathematical extrapolation. It was invented to solve a theoretical evolutionary dilemma," and "If 'Haldane's Dilemma' was correct, the majority of DNA must be non-functional."63 Jeffrey Tomkins and Jerry Bergman say new findings of function "present serious obstacles to the evolutionary necessity of a large fraction of the genome being 'junk'."64 Sanford, in Genetic Entropy, says Kimura's theory "led him to believe that most mutations are unselectable, and therefore most genetic sequences must be irrelevant,"59 and that "As the functional genome expands, the likelihood of neutral mutations shrinks."59
- ReMine concedes the central point. He writes that his book "has an entire chapter on neutral evolution, and argues that the vast majority of substitutions would be neutral, due to genetic drift, not selection."33 His dispute is about whether a neutral substitution has a cost, not about whether most substitutions are neutral.
- Genetic entropy is built on the nearly neutral theory. Sanford writes that "most sites are probably nearly-neutral" while doubting that any are "absolutely neutral."59 Carter describes genetic entropy as "an idea which leverages the beliefs and language of the neutralist position to make a strong case against it."29 The argument needs Ohta to be right: slightly harmful mutations must drift to fixation in populations as small as ours. Part 3 shows that this is the part of the theory the data support best, and that Kimura himself wrote that effectively neutral mutants are "in fact, very slightly deleterious" and "accumulate continuously in every species."31 Sanford overstates how the mainstream has treated them. He writes that in the standard account "everything in the 'near-neutral box' is redefined as being completely neutral, and is thereby dismissed."59 Kimura did not dismiss the box. He priced it, at one part in ten million of fitness per generation for the protein-coding genes, and answered it with occasional adaptive substitutions. Kondrashov priced the genome-wide case in 1995 and rejected that answer.42 So the question between Sanford and the mainstream is not whether the class exists, but how much of the genome lies in it and what offsets it.
- Some young-earth arguments need the neutral rate result. Carter uses the equality of mutation rate and substitution rate to date the human ancestors with pedigree mutation rates: "since most mutations are selectively neutral, genealogical mutation rates should be directly applicable to estimates of the timing of Y Chromosome Adam and Mitochondrial Eve."32 Elsewhere he writes that "The attempted rescue mechanism of neutral evolution theory, and the associated belief in 'junk DNA', failed as an escape mechanism."65 Both cannot be fully right unless he says which part of the theory failed; he does not.
- Nobody answered 2018 or Lewontin. No creationist reply to Kern and Hahn, and no creationist treatment of Lewontin's paradox, was found in searches of the main young-earth and design sites.
Is "neutral" the same as "junk"?
Partly. The creationist step from neutral theory to junk DNA is right about where the theory's genome-wide claim gets its weight, as Part 5 shows. It is wrong as a statement of what the theory requires. Kimura wrote that "the existence of selective constraint itself does not contradict the neutral theory,"1 and his 1979 model built constraint into the theory explicitly.31 A site can have a function and still accept most of its changes neutrally, and a stretch of DNA can do a job without its exact letters mattering, as Graur's "indifferent DNA" shows. What the theory cannot absorb is a large share of the genome whose exact letters matter strongly enough for selection to see. That is a narrower claim, and it is the one a design argument has to prove.
What follows for the mainstream
The theory is the mainstream's own statement that natural selection did not build most of the DNA differences between humans and chimpanzees. On the best current count, selection fixed about 4,000 protein changes on the human branch, with an unmeasured number in regulatory DNA, out of about 15 million. Popular accounts that present the human–chimpanzee difference count as a record of adaptation are describing something the field's own theory denies.
It also means the theory cannot be used as independent evidence that the genome is mostly junk. Its genome-wide success was partly assumed from the start, in Ohta's freely changing bulk and Kimura's dead pseudogenes. The independent evidence is the constraint measurements, and that is where the argument belongs.
What follows for the young-earth side
The concession is real. The figure on this site's Functional DNA page, that over 99.98% of mammal differences spread by chance, is close to the mainstream count: about 99.97% of human-branch changes were not adaptive at Uricchio's central estimate, and 99.6% on the most generous figures.8 A young-earth writer can cite the neutral theory for the claim that selection did not build most DNA differences, and should cite Uricchio's 4,000, not Kimura's 1968 arithmetic.
The inference needs care. "A highly neutral evolutionary process could create nothing better," as that page argues, applies to the DNA that changed: the 4% that differs between humans and chimpanzees. On common descent the shared 96% was inherited, not built in six million years, and the neutral theory says nothing about how it arose. The mainstream also proposes a neutral route to complexity, "constructive neutral evolution," a "unidirectional evolutionary ratchet leading to complexity";66 it is a scenario, not a measured rate, but it means the claim that neutral processes build nothing is contested.
Between created kinds, the neutral theory does not apply on a young earth. At the measured rate, neutral substitution changes about three sites in a million per lineage in 6,000 years.8 So none of the differences between created kinds can be fixed mutations, neutral or otherwise. Young-earth models say they were created, and Nathaniel Jeanson predicts that "the majority — if not the vast majority — of nuclear DNA differences among species within families (for vertebrates) ... are functional participants in the biology of each species, not functionally neutral."67 That is a clear, testable claim. The McDonald–Kreitman estimates for close relatives give a mixed answer: in fruit flies about half of protein changes and a fifth of intron changes look adaptive, while between mouse and rat, which share a family, the adaptive share is 19% to 38% in proteins and conserved elements and the rest of the genome is treated as neutral.40 Neither is a direct test, because "adaptive" and "functional" are not the same thing. The test the prediction calls for, measuring selection on the within-family differences themselves, has not been run.
Within kinds, drift has to work fast. Jeanson and Jason Lisle propose that created variety was sorted after the Flood, and that "genetic drift and other processes likely have diluted this concentrated heterozygosity to lower levels."68 Drift on that timetable is the neutral theory's own mechanism, and it has a speed limit. A neutral allele starting at half the population takes on average 4Ne ln 2, about 2.8 Ne generations, to fix, when it does.69 In 4,500 years that allows an effective size of up to about 1,600 for an animal that breeds every year, and about 160 for one that breeds every ten years.8 Small post-Flood populations meet that limit; large ones do not, unless selection rather than drift did the sorting. The Creation Model of Rapid Speciation article takes up how the models handle it.
Objections and responses
"Kern and Hahn showed the neutral theory is dead, so this article is out of date."
They argued that its universality is rejected, and on proteins in species with large populations they are right, as Part 3 shows. But they concede the human figure for proteins is near zero,6 and they set junk DNA aside as uncontested.6 Seven senior population geneticists replied that their view is "probably a minority view,"5 and the 2024 review from inside the field treats the question as one of proportions, not of life or death.2 For the human genome as a whole, nothing in the 2018 paper moves the conclusion.
"Selection reaches 95% of human variants, so nothing is neutral."
That figure measures linkage to selected sites, not selection on the sites themselves.56 Background selection lowers diversity near constrained DNA but does not change the rate at which neutral mutations fix between species.5 And Murphy and colleagues found that in humans sweeps had "little or no effect," so most of that reach is the removal of harmful mutations, which is the neutral theory's own second clause.21
"The 4,000 adaptive changes is an underestimate."
It may be. Eyre-Walker and Keightley warn that the human protein figure "could be as high as 40%,"37 and adaptive change in human regulatory DNA has hardly been measured. Take the generous case: 40% of the 29,925 human-branch protein changes, plus four times as many in noncoding DNA, the ratio Halligan and colleagues found in wild mice. That is about 60,000 adaptive changes, 0.4% of the human branch.8 The conclusion that most of the 15 million changes were not adaptive survives any figure the literature supports.
"Conservation only measures shared DNA under common descent. Under common design it means nothing."
For the cross-species estimates, that is a fair point, and this site's Functional DNA article makes it. But constraint can also be measured within one species, from which variants are missing in living people, and that needs no assumption about ancestry. The human-only estimate is 4% to 11% beyond protein-coding DNA.25 A design argument that relies on a much larger sequence-specific fraction has to show it in that kind of data.
"The neutral theory proves evolution cannot build function, so function refutes evolution."
The theory says selection did not fix most of the differences between species. It says nothing about the shared DNA, which on common descent was inherited, not built in the last six million years. The mainstream also proposes that neutral processes can build complexity, though this "constructive neutral evolution" is a scenario rather than a measured rate.66 The argument has force only if a large share of the DNA that differs turns out to be sequence-specific, which is the test in Part 5.
"Kimura himself said most of the genome must be irrelevant."
Sanford reads him that way.59 Kimura's own words are more careful: "the existence of selective constraint itself does not contradict the neutral theory."1 The fair version of the charge is narrower and correct: Kimura's showcase evidence, the fast-changing pseudogenes, rested on his opinion that they were "indeed dead."1
"Effectively neutral mutations are harmless by definition, so genetic entropy misreads the theory."
Not by the theory's own definition. Eyre-Walker and Keightley call neutrality "operational rather than functional": it depends on whether selection can act in a given population, "not solely on the effect of the mutation on fitness."41 Kimura wrote that effectively neutral mutants are "in fact, very slightly deleterious" and "accumulate continuously in every species."31 The theory and genetic entropy agree that the class exists and accumulates. They disagree about its size and about what offsets it: Kimura put the loss at one part in ten million per generation for protein-coding genes and trusted adaptive substitutions to cover it, while Kondrashov's genome-wide version is a paradox he could not resolve with beneficial mutations.42 Genetic Entropy prices that dispute.
"Synonymous mutations are mostly harmful (the 2022 yeast study), so the neutral yardstick is broken."
The yeast result has not held in human cells, and in apes about a fifth to two fifths of synonymous sites show weak constraint. Synonymous Mutations and Codon Usage works through it and finds the bias in dN/dS is modest, a factor of 1.26 to 1.57, and in the direction that makes protein change look less constrained, not more. It does not overturn the numbers here.
"The neutral theory can fit anything by adjusting population size, so it explains nothing."
Lewontin made a version of this complaint in 1974. But the theory made risky predictions, and some of them failed: protein evolution in flies is not mostly neutral, and diversity does not track population size at equilibrium. Others passed: pseudogenes change equally at all three codon positions,1 and slightly harmful changes are commoner in small populations, as Ohta predicted.19 A theory that can fail and sometimes does is testable.
"A young earth cannot be squared with any of this."
Part 4 and Part 6 price the two places it bears. Lewontin's paradox is what a recent, non-equilibrium history predicts, so it is no evidence against a young earth. But drift sorts created variety slowly in large populations: a neutral allele at half frequency needs about 2.8 Ne generations to fix,69 which fits only small populations in 4,500 years.
Issues and unanswered questions
- Adaptive change in human noncoding DNA is barely measured. Every human α in this article is for proteins. Regulatory DNA is most of the constrained genome, and the mouse data suggest it holds most adaptive changes. The human-branch count could rise several-fold without changing the conclusion that most change was not adaptive.
- The methods disagree. For human against chimpanzee, the same data give −0.41 by the plain McDonald–Kreitman test and 0.25 under a fuller model.15 Gossmann and colleagues in 2012 found the adaptive rate rising with population size; Galtier in 2016 found it flat. These are not settled.
- Lewontin's paradox is open. Linked selection falls short;20 non-equilibrium history is the leading candidate but has not been shown species by species. No young-earth writer was found to have addressed it.
- Jeanson's prediction is untested. Whether most nuclear differences within a created family are "functional participants" or neutral could be tested by measuring selection on those differences in, for example, mouse and rat or cat species. No one has done it.
- Drift's speed inside young-earth models. The fixation-time arithmetic in Part 6 assumes a constant effective size. Real post-Flood populations grew from a few founders; a model with growth and selection would give a better answer, and none has been published.
- ReMine's claim that neutral substitutions carry a cost is a definitional dispute with Brookfield and others that this article does not settle. A simulation counting the lost reproduction per neutral fixation would settle it.
- A citation defect in the 2018 paper. Kern and Hahn cite Boyko and colleagues for a human figure of zero; Boyko's abstract gives 10–20%.14 Li and colleagues 2008, their other source, was not checked.
- Carter's two positions on the neutral theory, that it failed as an escape from Haldane's dilemma and that its rate result dates Adam and Eve, need reconciling by their author.
- How much of the genome is very slightly harmful when changed is unmeasured. Conservation scans cannot see harm near the drift barrier, because such sites still change at more than half the neutral rate (Part 3). So the 3.5% to 26% constraint figures bound sequence-specific function that selection can act on, not the class the genetic-entropy argument needs. No current method bounds that class directly.
Verdict
Most DNA differences fixed between mammal species are neutral or nearly so, and the mainstream's own count puts the adaptive share on the human branch at a few hundredths of a percent; even the most generous reading leaves it under half a percent. That conclusion rests mainly on measurements showing that most of the genome is not under constraint, so the theory's genome-wide form does depend on most of the genome lacking sequence-specific function, as young-earth critics say. Inside proteins of species with large populations, the theory's original claim has failed, and the 2018 challenge is right about that while conceding the rest. For the design argument, the neutral theory is a real concession that selection did not build most differences, but it becomes evidence against evolution only if a large, sequence-specific function is shown in the DNA that differs, and the within-human data do not yet show one. For the genetic-entropy argument, the theory supplies the premise: its founder wrote that effectively neutral mutations are "in fact, very slightly deleterious" and accumulate in every species, and the dispute is over how many sites they touch and what offsets the loss.31
Sources
- 1.
Kimura, Motoo. The Neutral Theory of Molecular Evolution. Cambridge University Press. 1983. Kimura's book-length statement of the theory and its evidence fifteen years on: the mathematics of drift, the history of the selectionist debate, and the comparative data on protein, synonymous and pseudogene evolution. It is the fullest statement of the theory by its founder and the source most later critics, including creationists, quote. Mirror: Archive.org
- 1a.
Section Preface: “The neutral theory asserts that the great majority of evolutionary changes at the molecular level, as revealed by comparative studies of protein and DNA sequences, are caused not by Darwinian selection but by random drift of selectively neutral or nearly neutral mutants. The theory does not deny the role of natural selection in determining the course of adaptive evolution, but it assumes that only a minute fraction of DNA changes in evolution are adaptive in nature, while the great majority of phenotypically silent molecular substitutions exert no significant influence on survival and reproduction and drift randomly through the species.” ↩ ↩ ↩
- 1b.
Page 103: “(ii) Functionally less important molecules or parts of molecules evolve (in terms of mutant substitutions) faster than more important ones.” ↩
- 1c.
Page 103: “the probability of a mutation not being harmful and therefore selectively neutral is larger if the mutation occurs ina functionally less important molecule or a part of a molecule, and thus has a higher chance of being fixed in the population by random genetic drift.” Page 160: “It is known that this part evolves at about the rate 2.4 x 10⁻⁹ per amino acid site per year which is several times as fast as that of insulin.” Page 313: “We have also seen that in hemoglobins the surface part evolves about 10 times as fast as the heme pocket”. "ina" is the scan's text for "in a". The page 160 passage is about the middle segment of proinsulin, which is cut out when insulin is made. ↩
- 1d.
Page 160: “selective constraint (negative selection) may still be acting to some extent even in hypervariable regions, but the existence of selective constraint itself does not contradict the neutral theory as I have already emphasized in previous places in this book.” ↩ ↩ ↩ ↩
- 1e.
Page 178: “It is sometimes called a ‘dead gene’, although some have objected to this expression (see for example, Vanin et al., 1980) saying that a pseudogene may have an unknown but important function such as controlling the expression of normal genes. In my opinion, it is much more appropriate to regard pseudogenes as indeed dead; they have been liberated from the constraint of negative selection after they were created by gene duplication followed by loss of”. ↩ ↩ ↩ ↩
- 1f.
Page 179: “Unlike the ‘conservative’ mode of change that characterizes the evolution of many normal genes, base substitutions at the first and the second codon positions in the pseudogenes occur just as frequently as those at the third position.” ↩ ↩
- 1g.
Page 115: “We note also that accumulation of very slightly deleterious mutants by random drift is essentially equivalent to the deterioration of the environment, and definitely adaptive gene substitutions must occur from time to time to save the species from extinction.” ↩
- 1a.
- 2.
Galtier, Nicolas. “Half a century of controversy: the neutralist/selectionist debate in molecular evolution.” Genome Biology and Evolution 16(2):evae003. 2024. A short review of fifty years of the neutralist–selectionist debate in the two journals of the Society for Molecular Biology and Evolution, by a Montpellier population geneticist who has himself measured adaptive protein evolution in 44 animal species. It is the most balanced verdict from inside the field.
- 2a.
Section Introduction: “Interestingly, the 2 cornerstones of the neutral theory are still very much alive and important aspects of current molecular evolutionary thinking.” ↩
- 2b.
Section Within-Species Polymorphism and Lewontin's Paradox: “The issue, however, is no longer phrased in terms of corroborating versus challenging the neutral theory, but rather in terms of quantifying the relative contribution of neutral and selective processes to the observed patterns.” ↩ ↩
- 2c.
Section Substitution Rates and Protein Adaptation: “This led to estimates of the fraction of adaptive amino acid substitution that varied from nearly zero (Gossmann et al. 2010) up to roughly 90% (Galtier 2016)—an impressively wide range, likely in part explained by differences in the number and nature of analyzed genes, as well as in methodological choices.” ↩
- 2d.
Section Substitution Rates and Protein Adaptation: “Therefore, we have reached this peculiar situation where the analysis of amino acid substitution rates has both led to a clear rejection of a major claim of the neutral theory—that most observable changes are neutral—and at the same time indisputably demonstrated the importance of drift in protein evolution.” ↩ ↩ ↩ ↩
- 2e.
Section Why a Controversy?: “Enthused by their discovery, they slightly overstated the prevalence of neutrality.” ↩
- 2f.
Section Within-Species Polymorphism and Lewontin's Paradox: “These hypotheses account for the narrow range of the level of genetic polymorphism observed across species, while being compatible with the neutral theory.” ↩
- 2a.
- 3.
Kimura, Motoo. “Evolutionary rate at the molecular level.” Nature 217:624–626. Feb 17, 1968. The three-page letter that founded the neutral theory. Kimura extrapolated the amino-acid substitution rates of three proteins to the whole mammal genome, found about one nucleotide substitution every two years, judged that too fast for Haldane's cost of selection, and concluded that most substitutions must be nearly neutral. It also states that neutral changes accumulate at the mutation rate.
- 3a.
Page 624: “Calculating the rate of evolution in terms of nucleotide substitutions seems to give a value so high that many of the mutations involved must be neutral ones.” ↩ ↩
- 3b.
Page 625: “Thus the very high rate of nucleotide substitution which I have calculated can only be reconciled with the limit set by the substitutional load by assuming that most mutations produced by nucleotide replacement are almost neutral in natural selection.” ↩
- 3c.
Page 625: “This means that new alleles may be produced at the same rate per individual as they are substituted in the population in evolution.” ↩
- 3a.
- 4.
King, Jack Lester, and Thomas H. Jukes. “Non-Darwinian evolution.” Science 164(3881):788–798. May 16, 1969. The second founding paper of the neutral theory, reached independently of Kimura from protein chemistry: amino-acid composition, the genetic code, and the substitution patterns of a few well-sequenced proteins. It argues that most protein evolution is the drift of neutral changes and that most mammalian DNA codes for no protein.
- 4a.
Page 788: “Natural selection is the editor, rather than the composer, of the genetic message.” ↩ ↩
- 4b.
Page 790: “it is probable that not much more than 1 percent of mammalian DNA codes for proteins.” ↩
- 4c.
Page 792: “one can conclude that at least 90 percent of all substitutional mutations at the cytochrome-c locus are harmful and are rejected by natural selection.” ↩
- 4a.
- 5.
Jensen, Jeffrey D. et al. “The importance of the Neutral Theory in 1968 and 50 years on: a response to Kern and Hahn 2018.” Evolution 73(1):111–114. Jan 2019. A reply to Kern and Hahn by seven senior population geneticists, including two of the authors who first described background selection. It argues that the 2018 paper used too narrow a definition of the theory and leaned on protein data from a small part of the genome, and it sets out a five-point modern version of the theory.
- 5a.
Paragraph Abstract: “Although this is probably a minority view, it is important to evaluate such claims carefully in the context of current knowledge, as inaccuracies can sometimes morph into an accepted narrative for those not familiar with the underlying science. We here critically examine and ultimately reject Kern and Hahn’s arguments and assessment, and instead propose that it is now abundantly clear that the foundational ideas presented five decades ago by Kimura and Ohta are indeed correct.” ↩ ↩
- 5b.
Paragraph 1: “The Neutral Theory of Molecular Evolution asserts that most de novo mutations are either sufficiently deleterious in their effects on fitness that they have little chance of becoming fixed in the population, or are under such weak selection that they may become fixed as a result of genetic drift”. ↩
- 5c.
Paragraph 4: “Apart from the inherent uncertainty in these estimates (discussed by Fay 2011), it is misleading to use them to make the general claim that the Neutral Theory is insufficient to explain genome-wide patterns of variation and evolution; these inferred frequencies of adaptive substitutions mostly concern only the small fraction of the genome that codes for proteins (e.g., <2% of the human genome; see Lander et al. 2001).” ↩ ↩ ↩ ↩
- 5d.
Paragraph 4: “The circularity involved in ignoring the vast majority of neutral or nearly neutral substitutions across the genome, and then rejecting a significant role for neutrality, hardly justifies the need for the “selection theory of molecular evolution” advocated by Hahn (2008).” ↩
- 5e.
Paragraph 5: “In an explicit comparison between models of widespread purifying selection on weakly deleterious alleles versus recurrent positive selection on beneficial alleles, Lohmueller et al. (2011) found a much better fit of the former to the observed pattern in humans”. ↩
- 5f.
Paragraph 7: “Regardless of the precise interplay of the two forms of hitchhiking, background selection and selective sweeps, in shaping patterns of variation, it is important to note that neither affect the probability of fixation of neutral mutations (Birky and Walsh 1988), which determines the rate of neutral sequence evolution.” ↩ ↩ ↩
- 5g.
Section Modern version, point 1: “A large fraction of the genome of organisms studied to date is subject to mutations that are effectively neutral with respect to their fitness effects, and hence evolve under genetic drift.” ↩
- 5h.
Section Modern version, point 2: “The great majority of newly arising mutations that do affect fitness (i.e., non-neutral mutations) are deleterious, and the predominant mode of natural selection is purifying in nature, removing these deleterious mutations from populations.” ↩ ↩
- 5i.
Section Modern version, point 5: “In the absence of an appropriate null model accounting for these processes that are common to the genome as a whole, inappropriate adaptive story-telling will be likely to proliferate.” ↩
- 5a.
- 6.
Kern, Andrew D., and Matthew W. Hahn. “The neutral theory in light of natural selection.” Molecular Biology and Evolution 35(6):1366–1371. May 2, 2018. A perspective written for the theory's fiftieth anniversary by two population geneticists, at Rutgers and Indiana, who argue that natural selection dominates molecular evolution. It reviews the original evidence, McDonald–Kreitman estimates across species, and the correlation between recombination and diversity in 40 genomes, and concludes the theory has been rejected. It is the main modern challenge to the theory. Mirror: PMC
- 6a.
Page 1366: “On the 50th anniversary of the neutral theory of molecular evolution, we have been charged with the task of asking: how has the neutral theory fared in light of adaptive variation within and between species? In a word, poorly.” ↩
- 6b.
Page 1366: “We argue that the neutral theory was supported by unreliable theoretical and empirical evidence from the beginning, and that in light of modern, genome-scale data, we can firmly reject its universality.” ↩
- 6c.
Page 1366: “It is not simply a statement about the presence of neutral mutations, nor about the large fraction of eukaryotic genomes that are nonfunctional—neither of these assertions would be contested by current competing hypotheses.” ↩ ↩ ↩ ↩
- 6d.
Page 1366: “The neutral theory instead posits a positive thesis about nature: that differences between species are due to neutral substitutions (not adaptive evolution), and that polymorphisms within species are not only neutral but also have dynamics dominated by mutation-drift equilibrium.” ↩
- 6e.
Page 1366: “For instance, Kimura overestimates the number of protein-coding sites in the genome by two orders of magnitude (he uses 4×10⁹ bp, whereas the real number is closer to 3×10⁷ bp in the largest genomes). Even given the data on rates Kimura had available, fixing this number alone would remove the conflict between the rate of protein evolution at the nucleotide level and Haldane’s (1957) upper limit of 300 years per adaptive substitution.” Haldane's limit was 300 generations, not 300 years. ↩ ↩
- 6f.
Page 1367: “Accumulating evidence from additional proteins, coupled with better analyses (Langley and Fitch 1974), soon showed that this constancy was an illusion”. ↩
- 6g.
Page 1367: “The first such studies were carried out in Drosophila melanogaster and D. simulans, finding that ∼50% of all amino acid substitutions have been fixed by positive selection”. ↩
- 6h.
Page 1367: “In humans, though the overall fraction of amino acid substitutions fixed by positive selection is estimated to be zero (Boyko et al. 2008; Li et al. 2008), careful functional characterization of individual proteins has revealed that a large fraction of all genes which interact with pathogens show pervasive evidence for positive selection”. ↩ ↩ ↩
- 6i.
Page 1368: “First, these results imply that almost no loci are free from the effects of selection, in any organism.” "These results" are the correlations between recombination rate and diversity in 40 genomes plotted in their Figure 1, from data of Corbett-Detig and colleagues 2015. ↩
- 6j.
Page 1368: “An equilibrium between mutation and drift is a central tenet of the neutral theory (Kimura and Ohta 1971); therefore, current data appear to be fundamentally incompatible with the neutral theory.” ↩
- 6k.
Page 1369: “As a consequence, we believe that the neutral theory has been overwhelmingly rejected, and that as a field we must continue to develop alternate theories of molecular evolution.” ↩
- 6a.
- 7.
Nei, Masatoshi. “Selectionism and neutralism in molecular evolution.” Molecular Biology and Evolution 22(12):2318–2342. 2005. A long historical review by one of the founders of molecular population genetics, who defends both Haldane's cost of selection and the neutral theory. It recomputes Kimura's 1968 argument with genes and amino acids as the unit of selection instead of every DNA letter.
- 7a.
Section Cost of Natural Selection: “However, Kimura's paper had at least two deficiencies. First, in the computation of the cost of natural selection, he assumed that all nucleotides in the genome were subject to natural selection. In practice, the unit of selection should be a gene or an amino acid because noncoding regions of DNA are largely irrelevant to the evolution of proteins and organisms.” “In contrast, if we consider an amino acid as the unit of selection and each gene encodes on average 450 amino acid sites (Zhang 2000), the average rate of amino acid substitution will be one substitution every 636 (=2.86 × 10⁵/450) years. This rate is about two times higher than Haldane's upper limit.” ↩ ↩
- 7a.
- 8.
This article's arithmetic, in research/pricing.ts, which rebuilds every figure from its published inputs and checks itself. Human-branch single-letter changes: 35 million differences between the two genomes19, times the fixed share (1.06% of 1.23%), halved for the two branches, gives about 15.1 million. Adaptive protein changes: 29,925 human-branch protein-changing substitutions times α = 0.135 gives 4,040 (2,873 to 5,087 on the 0.096 to 0.17 interval); with Boyko's 10–20%, 2,993 to 5,985. Shares of 15.1 million: 0.027% at 0.135, 0.040% at 20%. Generous case: 40% of 29,925 is 11,970, plus four times as many in noncoding DNA, gives 59,850, or 0.40%. So the share not adaptive is about 99.97% at α = 0.135 and 99.60% in the generous case. New point mutations per child: 2 × 3.1 × 10⁹ sites × 1.25 × 10⁻⁸ = 77.5; times a sequence-specific fraction of 3.5% gives 2.7, of 30% gives 23.3 (scale only). If 30% of the genome were sequence-specific and changes fell at random, 4.5 million of the 15.1 million would land in it. Neutral change in g generations is g × μ: 240 generations give 3.0 × 10⁻⁶ per site. Time to mutation–drift balance, about 4Ne generations: 40,000 to 80,000 for Ne = 10,000 to 20,000, or 1.0 to 2.3 million years at 25 to 29 years a generation. Drift barrier 1/(2Ne): 2.5 × 10⁻⁵ at Ne = 20,000. ↩ ↩ ↩ ↩ ↩ ↩ ↩ ↩ ↩ ↩ ↩ ↩ ↩
- 9.
Uricchio, Lawrence H. et al. “Exploiting selection at linked sites to infer the rate and strength of adaptation.” Nature Ecology & Evolution 3(6):977–984. 2019. A statistical method that separates weakly and strongly beneficial changes, applied to protein-coding DNA of 661 people of African ancestry from the 1000 Genomes Project, with the human-specific changes since the chimpanzee split. It estimates that 13.5% of human protein-changing substitutions were adaptive, most of them weakly so, and it is the source of this article's count of about 4,000.
- 9a.
Section Human adaptation rate is shaped by linked selection: “We estimate that the total α̂ = 0.135, nearly twice the estimate obtained with the same dataset using the original aMK approach”. Table 1: “29925”. The source prints α with a hat, for an estimate. Table 1, whole-genome row: 29,925 nonsynonymous fixations, α 0.135 with a 95% interval of 0.096 to 0.17. The methods describe these as "human-specific fixed substitutions in coding sequences since divergence with chimpanzees." ↩ ↩ ↩
- 9a.
- 10.
Smith, Nick G. C., and Adam Eyre-Walker. “Adaptive protein evolution in Drosophila.” Nature 415:1022–1024. 2002. A genome-scale McDonald–Kreitman estimate for fruit flies, using the divergence of Drosophila simulans and D. yakuba. It is the first estimate that about half of the protein changes between fruit-fly species were adaptive.
- 11.
Charlesworth, Jane, and Adam Eyre-Walker. “The rate of adaptive evolution in enteric bacteria.” Molecular Biology and Evolution 23(7):1348–1356. 2006. A McDonald–Kreitman analysis of 410 genes shared by six Escherichia coli and six Salmonella enterica genomes. It found more than half of the protein changes between the two lineages adaptive, while warning that the estimate rises with the frequency cut-off and does not settle.
- 11a.
Page 1348: “We estimate that more than 50% of amino acid substitutions in this set of genes have been fixed by positive selection between the E. coli and S. enterica lineages.” Page 1352: “this is 56%, with a lower confidence interval (CI) of 45%. Using polymorphism data from S. enterica, it is 34%, with a lower CI of 14%.” ↩ ↩
- 11a.
- 12.
Halligan, Daniel L. et al. “Evidence for pervasive adaptive protein evolution in wild mice.” PLoS Genetics 6(1):e1000825. 2010. Protein-coding variation in 15 wild house mice of the subspecies Mus musculus castaneus from India, compared with another mouse species and the rat. It found a majority of protein changes adaptive in this large-population mammal, against low figures for hominids.
- 13.
Our summary of Galtier's per-species table:70 44 species pairs, median α 0.65, minimum 0.18, maximum 0.93; 31 pairs at or above 0.5. The twelve mammal pairs have a median of 0.36, and the human pair, with chimpanzee as outgroup, 0.24. ↩ ↩
- 14.
Boyko, Adam R. et al. “Assessing the evolutionary impact of amino acid mutations in the human genome.” PLoS Genetics 4(5):e1000083. May 30, 2008. A distribution of fitness effects fitted to 47,576 human variants in 11,404 genes, used to predict how many human–chimpanzee protein differences were fixed by selection. It is one of the two papers Kern and Hahn cite for a human figure of zero, though its own estimate is 10–20%.
- 14a.
Page 1: “Our results are consistent with 10–20% of amino acid differences between humans and chimpanzees having been fixed by positive selection with the remainder of differences being neutral or nearly neutral.” ↩ ↩ ↩ ↩
- 14b.
Page 1: “we predict 27–29% of amino acid changing (nonsynonymous) mutations are neutral or nearly neutral (|s|<0.01%), 30–42% are moderately deleterious (0.01%<|s|<1%), and nearly all the remainder are highly deleterious or lethal (|s|>1%).” ↩
- 14a.
- 15.
Zhen, Ying et al. “Greater strength of selection and higher proportion of beneficial amino acid changing mutations in humans compared with mice and Drosophila melanogaster.” Genome Research 31(1):110–120. 2021. A comparison of methods for estimating the adaptive share in humans, mice and fruit flies, which shows that the plain McDonald–Kreitman test is thrown off by the outgroup and ancestral population sizes. Its preferred model gives humans an adaptive share of at least 20%, and its table shows how far the answer moves with the method.
- 15a.
Page 110: “estimate that at least 20% of nonsynonymous substitutions between humans and an outgroup species were fixed by positive selection.” ↩
- 15b.
Page 111, table 1: Human–chimpanzee: plain MK test with all variants −0.41, with variants under 20% frequency removed −0.01, Simple model 0.11, Complex model 0.25, DFE-alpha 0.24. D. melanogaster–D. simulans: −0.13, 0.49, 0.53, 0.60, 0.71. Mouse–rat: 0.25, 0.40, 0.45, 0.41, 0.51. ↩ ↩ ↩ ↩ ↩
- 15a.
- 16.
Gossmann, Toni I. et al. “Genome wide analyses reveal little evidence for adaptive evolution in many plant species.” Molecular Biology and Evolution 27(8):1822–1832. 2010. More than 2,400 genes from nine pairs of plant species analysed with the Eyre-Walker and Keightley method. It found little sign of adaptive protein evolution in any plant except sunflowers.
- 16a.
Page 1822: “we find little evidence of adaptive amino acid substitution in any comparison except sunflowers. This may be because many plant species have modest effective population sizes.” ↩ ↩
- 16b.
Page 1826: “The mean estimate of α across all ten species-wide comparisons is −0.18 (standard error [SE] = 0.47). This is not significantly different from zero”. ↩
- 16a.
- 17.
Ohta, Tomoko. “Slightly deleterious mutant substitutions in evolution.” Nature 246:96–98. Nov 9, 1973. The letter that introduced the nearly neutral theory. Ohta proposed that slightly harmful mutations, not only strictly neutral ones, drift to fixation, and that they do so more often in small populations, which explains why protein clocks tick per year rather than per generation.
- 17a.
Page 96: “The concept of neutral mutant substitution in the population by random genetic drift can be extended to include random fixation of very slightly deleterious mutations which have more chance of being selected against than of being selected for3,4. If this class of mutant substitution is important, we can predict that the evolution is rapid in small populations or at the time of speciation5.” The numerals 3,4 and 5 are reference marks, left out of the quotation in the article. ↩ ↩
- 17a.
- 18.
Galtier, Nicolas. “Adaptive protein evolution in animals and the effective population size hypothesis.” PLoS Genetics 12(1):e1005774. Jan 11, 2016. A new, uniform data set of protein-coding genes from 44 pairs of non-model animal species, from sea urchins and earthworms to penguins and primates, analysed with one method. It found the adaptive share above one half in most animals and low in primates, and traced its link to population size mainly to fewer nearly neutral changes in large populations.
- 18a.
Page 1: “The proportion of adaptive amino-acid substitution is found to be positively correlated to species effective population size. This relationship, however, appears to be primarily driven by a decreased rate of nearly-neutral amino-acid substitution because of more efficient purifying selection in large populations.” ↩ ↩
- 18b.
Page 12: “Observing a high α therefore seems to be the rule, and a low α the exception, making humans and apes a peculiar, intriguing case.” ↩
- 18c.
Page 14: “We report a wide range of estimated α (minimum: 0.18; maximum: 0.93) and ωa (minimum: 0.023; maximum: 0.15) among animals, without obvious connection to species biology or taxonomy.” ↩
- 18a.
- 19.
Chimpanzee Sequencing and Analysis Consortium. “Initial sequence of the chimpanzee genome and comparison with the human genome.” Nature 437:69–87. Sep 1, 2005. The first comparison of the whole human and chimpanzee genomes. It counted the differences, measured the ratio of protein-changing to silent change in 13,454 matched genes, and compared it with mice and rats. Its protein results are the standard reference for how much of human–chimpanzee protein divergence is neutral, harmful or adaptive.
- 19a.
Page 77: “The K A/K S ratio is markedly lower for murids than for hominids (q murid < 0.13 compared with q hominid < 0.20) (Fig. 9). This implies that there is an ,35% excess of the amino-acid-changing mutations in the two hominids, relative to the two murids. Excess amino acid divergence may be explained by either increased adaptive evolution or relaxation of evolutionary constraints. As shown in the next section, the latter seems to be the principal explanation.” The text layer renders the printed ω as q, ≈ as < and ∼ as a comma; the printed page reads ω murid ≈ 0.13, ω hominid ≈ 0.20 and ∼35%. ↩ ↩ ↩
- 19b.
Page 76: “Under the assumption that synonymous mutations are selectively neutral, the results imply that 77% of amino acid alterations in hominid genes are sufficiently deleterious as to be eliminated by natural selection.” ↩ ↩
- 19c.
Page 77: “A total of 585 of the 13,454 human–chimpanzee orthologues (4.4%) have observed K A/K I . 1”. “Simulations show that estimates of K A/K I . 1 would be expected to occur simply by chance in at least 263 cases if purifying selection is allowed to act non-uniformly across genes”. ". 1" is the text layer's rendering of "> 1". ↩
- 19d.
Page 69: “constituting approximately thirty-five million single-nucleotide changes, five million insertion/deletion events, and various chromosomal rearrangements.” “Single-nucleotide substitutions occur at a mean rate of 1.23% between copies of the human and chimpanzee genome, with 1.06% or less corresponding to fixed divergence between the species.” ↩
- 19a.
- 20.
Buffalo, Vince. “Quantifying the relationship between genetic diversity and population size suggests natural selection cannot explain Lewontin's Paradox.” eLife 10:e67509. Aug 19, 2021. A comparison of genetic diversity in 172 animal taxa with new estimates of their census sizes, followed by a test of whether models of linked selection could squeeze diversity into its observed narrow range. It concludes that they cannot, and points to past changes in population size instead.
- 20a.
Page 1: “germline mutation rates only differ 10-fold”. ↩
- 20b.
Page 1: “However, I show that even assuming strong and abundant selection, models of linked selection are unlikely to explain the observed relationship between diversity and census sizes across species.” ↩ ↩
- 20c.
Page 6, Figure 2 caption: “Pairwise diversity (data from Leffler et al., 2012, Corbett-Detig et al., 2015, and Romiguier et al., 2014), which varies over three orders of magnitude, shows a weak relationship with approximate population size, which varies over 12 orders of magnitude.” ↩ ↩
- 20d.
Page 14: “Given that I find that models of linked selection are incapable of explaining the observed relationship between Nc and p, this supports the hypothesis the diversity across species are shaped primarily by past demographic fluctuations.” "p" is the text layer's rendering of π. ↩
- 20a.
- 21.
Murphy, David A. et al. “Broad-scale variation in human genetic diversity levels is predicted by purifying selection on coding and non-coding elements.” eLife 11:e76065. 2022. A joint model of background selection and sweeps fitted to human genomes from the 1000 Genomes Project, controlling for variation in mutation rate. Background selection alone explained about 60% of the variation in diversity at the megabase scale, and sweeps added nothing detectable.
- 22.
Charlesworth, Brian, and Jeffrey D. Jensen. “How can we resolve Lewontin's paradox?” Genome Biology and Evolution 14(7):evac096. 2022. A review of the candidate explanations for Lewontin's paradox by two defenders of the neutral theory, with worked numbers for each. It judges linked selection insufficient on its own and favours populations that are still far below their equilibrium diversity after growing.
- 22a.
Page 8: “if the step change had involved a change in Ne to 10⁹ from a much lower value that is compatible with the modest π values seen in most species, the time for π to approach ¾ would be of the order of 10⁹ generations.” The same page gives a step change 5 million years ago, at ten generations a year, as producing π = 0.036, typical of invertebrates. ↩ ↩
- 22b.
Page 8: “The above calculations show that populations which have expanded from much smaller numbers in even the quite distant past may still be far from their very high equilibrium diversities. This scenario thus seems an excellent candidate for resolving LP.” ↩ ↩
- 22c.
Page 12: “This estimate would imply that Tn, the purely neutral value of the coalescent time, is about 5-fold larger than the value indicated by current diversity levels; this is far smaller than the factor of 100 or more suggested by the probable current population size of D. melanogaster (Karasov et al. 2010; Buffalo 2021), and so is insufficient to resolve LP for this species.” ↩ ↩
- 22a.
- 23.
Ohta, Tomoko. “The nearly neutral theory of molecular evolution.” Annual Review of Ecology and Systematics 23:263–286. 1992. Ohta's own review of the nearly neutral theory after twenty years: its population-size predictions, the generation-time effect, and the data on proteins and genomic DNA. It states the two-tier view of the genome that this article examines in Part 5.
- 23a.
Page 272: “In other words, I propose that a substantial fraction of mutant substitutions at the molecular level are caused by random fixation of very slightly deleterious mutations.” ↩
- 23b.
Page 272: “The fundamental idea is that most genomic DNA of higher organisms can freely accumulate base substitutions (80), i.e. most new mutations are selectively neutral. On the other hand, amino acid substitutions are more likely to be influenced by natural selection, i.e. many of them may be regarded as nearly neutral, or very slightly deleterious.” "(80)" is a reference mark, left out of the quotation in the article. ↩ ↩ ↩
- 23a.
- 24.
Rands, Chris M. et al. “8.2% of the human genome is constrained: variation in rates of turnover across functional element classes in the human lineage.” PLoS Genetics 10(7):e1004525. Jul 24, 2014. An estimate of the constrained fraction of the human genome from the shortfall of insertions and deletions across many mammal comparisons, allowing for constrained sequence to be gained and lost over time. It finds that most constrained noncoding sequence in humans is not shared with mouse.
- 25.
Kellis, Manolis et al. “Defining functional DNA elements in the human genome.” PNAS 111(17):6131–6138. Apr 29, 2014. A perspective by leading ENCODE scientists on the genetic, evolutionary and biochemical ways of defining function, written after the 2012 controversy. It gives the upper range of the constraint estimates the mainstream accepts.
- 25a.
Section Evolutionary approach: “estimates that incorporate alternate references, shape-based constraint (88), evolutionary turnover (89), or lineage-specific constraint (90) each suggests roughly two to three times more constraint than previously (12–15%), and their union might be even larger”. “human population studies suggest that an additional 4–11% of the genome may be under lineage-specific constraint after specifically excluding protein-coding regions”. ↩ ↩ ↩
- 25a.
- 26.
Christmas, Matthew J. et al. “Evolutionary constraint and innovation across hundreds of placental mammals.” Science 380(6643):eabn3943. Apr 28, 2023. The Zoonomia Consortium's estimate of how much of the human genome is conserved across 240 placental mammals, measured against neutrally evolving repeats. Its authors call the figure a minimum. Mirror: PMC
- 27.
Sullivan, Patrick F. et al. “Leveraging base-pair mammalian constraint to understand genetic variation and human disease.” Science 380(6643):eabn2937. 2023. A Zoonomia Consortium paper that scores every base of the human genome for conservation across 240 placental mammals and tests the scores against human variation and disease. It gives the share of individual bases that are significantly constrained. Mirror: PMC
- 28.
Graur, Dan et al. “On the immortality of television sets: "function" in the human genome according to the evolution-free gospel of ENCODE.” Genome Biology and Evolution 5(3):578–590. 2013. A polemical critique of the ENCODE project's claim that 80% of the human genome is functional, by an evolutionary biologist who defends junk DNA. It argues from the selected-effect meaning of function and from mutational load.
- 28a.
Section Junk DNA Is Dead! Long Live Junk DNA: “Indifferent DNA refers to DNA sites that are functional, but show no evidence of selection against point mutations. Deletion of these sites, however, is deleterious, and is subject to purifying selection. Examples of indifferent DNA are spacers and flanking elements whose presence is required but whose sequence is not im”. ↩
- 28b.
Section Junk DNA Is Dead! Long Live Junk DNA: “Second, there exists a misconception among functional genomicists that the evolutionary process can produce a genome that is mostly functional.” ↩ ↩
- 28a.
- 29.
Carter, Robert. “Mutations are more harmful than we thought: silencing the 'wobble' in the codon table.” Creation Ministries International. Aug 23, 2022. A popular article reading a 2022 yeast study of synonymous mutations as support for genetic entropy and against the neutral theory.
- 30.
Kim, Bernard Y. et al. “Inference of the distribution of selection coefficients for new nonsynonymous mutations using large samples.” Genetics 206(1):345–361. May 2017. A fit of the distribution of fitness effects of new amino-acid mutations to the allele frequencies of about 1,300 Europeans in each of two exome data sets and 432 in the 1000 Genomes Project. It finds more weakly harmful and fewer strongly harmful mutations than earlier fits from smaller samples. The model has no beneficial class, and the authors note that slightly helpful mutations would be counted in its neutral share.
- 30a.
Page 345, abstract: “Our results suggest that nearly neutral forces play a larger role in human evolution than previously thought.” Page 351, right column: “We infer a strongly leptokurtic distribution where there are many neutral and nearly neutral mutations (i.e., 34–37% of new mutations have |s| < 10⁻⁴), as well as a class of strongly deleterious mutations (i.e., 15–22% of new mutations have |s| > 10⁻²).” Page 354, table 3: At the authors' preferred mutation rate, the share of new amino-acid mutations below |s| = 10⁻⁵ is 0.237, 0.263 and 0.242 in the three data sets. ↩ ↩
- 30a.
- 31.
Kimura, Motoo. “Model of effectively neutral mutations in which selective constraint is incorporated.” PNAS 76(7):3440–3444. 1979. Kimura's model of a continuous distribution of harmful effects, in which the share of mutations that behave as neutral depends on how strongly a site is constrained and on population size. It is the model Sanford's Genetic Entropy draws its curve of mutation effects from. ↩
- 31a.
Page 3444, left column: “Finally, there is one biological problem that we have to consider. Under the present model, effectively neutral, but, in fact, very slightly deleterious mutants accumulate continuously in every species. The selective disadvantage of such mutants (in terms of an individual's survival and reproduction—i.e., in Darwinian fitness) is likely to be of the order of 10⁻⁵ or less, but with 10⁴ loci per genome coding for various proteins and each accumulating the mutants at the rate of 10⁻⁶ per generation, the rate of loss of fitness per generation may amount to 10⁻⁷ per generation. Whether such a small rate of deterioration in fitness constitutes a threat to the survival and welfare of the species (not to the individual) is a moot point, but this will easily be taken care of by adaptive gene substitutions that must occur from time to time (say once every few hundred generations).” ↩ ↩ ↩ ↩ ↩ ↩
- 31a.
- 32.
Carter, Robert W. “Genealogical vs. phylogenetic mutation rates: answering a challenge.” Proceedings of the Ninth International Conference on Creationism, 168–180. 2023. A conference paper defending the use of pedigree mutation rates to date the most recent common ancestors of human Y chromosomes and mitochondria, against the objection that long-term rates are slower.
- 33.
ReMine, Walter James. “Response to talk.origins article CB121.” Saint Paul Science. Nov 2006. ReMine's reply to the TalkOrigins summary of his Haldane's-dilemma argument, in which he corrects how his book is described.
- 34.
Bedford, Trevor et al. “Overdispersion of the Molecular Clock Varies Between Yeast, Drosophila and Mammals.” Genetics 179(2):977–984. 2008. A genome-wide test of how irregular the protein clock is, counting amino-acid substitutions in 3,081 yeast, 7,174 fruit-fly and 8,065 mammal proteins on fixed species trees. The clock is overdispersed in all three groups, most in mammals, and the overdispersion falls as effective population size rises. The authors read the pattern as compatible with nearly neutral evolution and do not take it as a sign of adaptation.
- 34a.
Page 979, table 1: Mean index of dispersion R(t) of amino-acid substitution counts: yeast 2.0993 over 3,081 proteins, Drosophila 4.1892 over 7,174, mammals 6.4790 over 8,065. A steady random clock gives 1. ↩
- 34b.
Page 977, abstract: “we show that the extent of overdispersion shows a strong negative correlation with the effective population size of these organisms”. Page 981: “We also emphasize that the finding of overdispersion does not necessarily imply positive selection or adaptive evolution; overdispersion is compatible with nearly neutral evolution occurring in a heterogeneous fashion.” ↩
- 34a.
- 35.
Charlesworth, Brian et al. “The effect of deleterious mutations on neutral molecular variation.” Genetics 134(4):1289–1303. 1993. The theory paper that described background selection: selection removing harmful mutations also removes the neutral variation linked to them, lowering diversity most where recombination is rare. It is the mechanism Jensen and colleagues offer in place of adaptive sweeps. ↩
- 36.
McDonald, John H., and Martin Kreitman. “Adaptive protein evolution at the Adh locus in Drosophila.” Nature 351:652–654. 1991. The letter that introduced the McDonald–Kreitman test, comparing changes that alter the protein with silent ones, within and between species, at the alcohol dehydrogenase gene of three fruit-fly species. It found an excess of protein-changing differences between species and attributed them to selection.
- 36a.
Page 652: “If the observed substitutions are neutral, the ratio of replacement to synonymous fixed differences between species should be the same as the ratio of replacement to synonymous polymorphisms within species. DNA sequence data on the Adh locus (encoding alcohol dehydrogenase, EC 1.1.1.1) in three species in the Drosophila melanogaster species subgroup do not fit this expectation; instead, there are more fixed replacement differences between species than expected.” ↩ ↩
- 36a.
- 37.
Eyre-Walker, Adam, and Peter D. Keightley. “Estimating the rate of adaptive molecular evolution in the presence of slightly deleterious mutations and population size change.” Molecular Biology and Evolution 26(9):2097–2108. 2009. A method that corrects the McDonald–Kreitman test for slightly harmful mutations, applied to fruit flies and to humans compared with macaques. It gives the standard figures of about half of fly protein changes adaptive and little evidence in humans, and warns that the human figure may be biased down.
- 37a.
Page 2097: “In Drosophila genes, we estimate that approximately 50% of amino acid substitutions and approximately 20% of substitutions in introns are adaptive. In protein-coding and noncoding data sets of humans, comparison to macaque sequences reveals little evidence for adaptive substitutions. However, the true frequency of adaptive substitutions in human-coding DNA could be as high as 40%, because estimates based on current polymorphism may be strongly downwardly biased by a decrease in the effective population size along the human lineage.” ↩ ↩ ↩
- 37a.
- 38.
Andolfatto, Peter. “Adaptive evolution of non-coding DNA in Drosophila.” Nature 437:1149–1152. 2005. An extension of the McDonald–Kreitman test to fruit-fly introns, intergenic DNA and untranslated regions. It found that much of that noncoding DNA is constrained and that a substantial share of its divergence was adaptive.
- 38a.
Page 1149: “However, I also use an extension to the McDonald-Kreitman test to show that a substantial fraction of the nucleotide divergence in these regions was driven to fixation by positive selection (about 20% for most intronic and intergenic DNA, and 60% for UTRs). On the basis of these observations, I suggest that a large fraction of the non-translated genome is functionally important and subject to both purifying selection and adaptive evolution.” A reference mark after "test" is left out. ↩ ↩
- 38a.
- 39.
Phifer-Rixey, Megan et al. “Adaptive evolution and effective population size in wild house mice.” Molecular Biology and Evolution 29(10):2949–2955. 2012. The adaptive share compared across three house-mouse subspecies with different population sizes, a test of the population-size effect inside one species. Only the largest subspecies showed a clearly nonzero share.
- 39a.
Page 2949: “There was significant evidence of positive selection in M. m. castaneus, the lineage with the largest Ne, with a estimated to be approximately 40%.” “Interestingly, the higher estimate of a for M. m. castaneus appears to reflect not only more adaptive fixations but also more effective purifying selection.” The abstract gives 13% for M. m. domesticus and 12% for M. m. musculus. ↩ ↩
- 39a.
- 40.
Halligan, Daniel L. et al. “Contributions of protein-coding and regulatory change to adaptive molecular evolution in murid rodents.” PLoS Genetics 9(12):e1003995. 2013. Whole genomes of ten wild mice, used to estimate the adaptive share in protein sites, untranslated regions and conserved noncoding elements. It found the share lower in noncoding elements but the number of adaptive changes there several times larger, because those elements are a much bigger target.
- 40a.
Page 2, table 1: Fraction adaptive (α): zero-fold protein sites 0.32, two-fold 0.38, UTRs 0.19, conserved noncoding elements 0.19. Adaptive substitutions: 0.32, 0.063, 0.12 and 1.4 Mb. The abstract: "we estimate that there are at least four times as many adaptive substitutions in CNEs and UTRs as in proteins." ↩ ↩ ↩
- 40a.
- 41.
Eyre-Walker, Adam, and Peter D. Keightley. “The distribution of fitness effects of new mutations.” Nature Reviews Genetics 8(8):610–618. Aug 2007. A review of the distribution of fitness effects of new mutations, drawing on mutagenesis experiments and on the frequencies of variants in natural populations. It defines neutrality by whether selection can act in a given population, not by zero effect, and estimates the effectively neutral share of amino-acid mutations at under 30% in hominids, at most 16% in fruit flies and 2.8% in enteric bacteria. It is the source Sanford's Genetic Entropy quotes on neutrality. Mirror: Homepages.ed.ac.uk
- 41a.
Page 612, right column: “The first point to make is one of definition; it seems unlikely that any mutation is truly neutral in the sense that it has no effect on fitness. All mutations must have some effect, even if that effect is vanishingly small. However, there is a class of mutations that we can term effectively neutral. These are mutations for which Nes is much less than 1, the fate of which is largely determined by random genetic drift. As such, the definition of neutrality is operational rather than functional; it depends on whether natural selection is effective on the mutation in the population or the genomic context in which it segregates, not solely on the effect of the mutation on fitness.” Reference numerals after “drift” are left out. ↩ ↩
- 41b.
Page 613, left column, foot, into right column: “In hominids, which seem to have effective population sizes in the range of 10,000 to 30,000 … the ratio dn/ds is less than 0.3 …, and this suggests that fewer than 30% of amino-acid-changing mutations are effectively neutral.” “if we infer that at most 30% of non-synonymous mutations are neutral in humans, this implies that at least 70% are deleterious.” The omissions are reference marks. ↩
- 41a.
- 42.
Kondrashov, Alexey S. “Contamination of the genome by very slightly deleterious mutations: why have we not died 100 times over?” Journal of Theoretical Biology 175:583–594. Aug 21, 1995. A theoretical paper by a mainstream population geneticist on mutations too weakly harmful for selection to remove but numerous enough to matter across a whole genome. It calculates that if a tenth of all sites were in that range a typical individual would carry about 100 lethal equivalents, calls the result a paradox, and weighs resolutions such as soft selection and synergistic epistasis. It is the mainstream statement of the problem the genetic-entropy argument rests on.
- 42a.
Page 583, abstract: “Mutations with s within this range are neutral enough to accumulate almost freely, but are still deleterious enough to make an impact at the level of the whole genome.” “This paradox cannot be resolved by invoking beneficial mutations or environmental fluctuations. Several possible resolutions are considered, including soft selection and synergistic epistasis among very slightly deleterious mutations.” The range is selection coefficients between one over the genome size and one over four times the effective population size. ↩ ↩ ↩ ↩
- 42a.
- 43.
This article's arithmetic, in research/pricing.ts. Kimura's loss of 10⁻⁷ per generation31, compounded, is 1 − (1 − 10⁻⁷) to the power 240 ≈ 0.0024% over 240 generations, about 6,000 years, and 1 − (1 − 10⁻⁷) to the power 200,000 ≈ 1.98% over 200,000 generations, roughly the human line since the chimpanzee split on mainstream dates. A harmful mutation with S = 4Nes fixes at S / (e to the power S, minus 1) times the rate of a neutral one, from Kimura's fixation probability for a new mutant: 0.95 at S = 0.1, 0.58 at S = 1, which is the upper edge of Kondrashov's range42, and 0.0005 at S = 10. ↩ ↩
- 44.
Huber, Christian D. et al. “Population genetic models of GERP scores suggest pervasive turnover of constrained sites across mammalian evolution.” PLoS Genetics 16(5):e1008827. May 29, 2020. A simulation study of how well GERP, a widely used conservation score, detects sites under purifying selection, fitted to the human genome's GERP scores. It finds the score tracks strong selection but loses power for weak selection and when sites change function over time, estimates 4.51% of the noncoding human genome under purifying selection, and finds the usual thresholds miss over half of it. It bears on what the constraint figures can and cannot see. Mirror: PMC
- 44a.
Page 1, abstract: “We show that mutations in 4.51% of the noncoding human genome are under purifying selection and that most of this sequence has likely experienced changes in selection coefficients throughout mammalian evolution.” ↩ ↩
- 44b.
Page 8: “There is no improvement in statistical power with increasing tree size for very weak selection (Nes < -0.25).” Page 1, abstract: “Commonly used GERP score thresholds miss over half of the noncoding sites in the human genome where mutations are under purifying selection.” The first sentence describes simulations in which selected sites gain and lose function over time at the rate estimated for noncoding elements. ↩
- 44a.
- 45.
Mouse Genome Sequencing Consortium. “Initial sequencing and comparative analysis of the mouse genome.” Nature 420:520–562. Dec 5, 2002. The first mouse genome paper, which compared 12,845 matched mouse and human genes. Its ratio of protein-changing to silent change is the standard measure of how strongly mammal proteins are held in place.
- 45a.
Page 543: “the median KA /KS ratio was 0.115”. For 12,845 mouse–human one-to-one orthologues. ↩
- 45a.
- 46.
Kosiol, Carolin et al. “Patterns of positive selection in six mammalian genomes.” PLoS Genetics 4(8):e1000144. 2008. A gene-by-gene scan for positive selection across human, chimpanzee, macaque, mouse, rat and dog. It counts genes with detectable adaptation, not the share of substitutions that were adaptive, and it says its power is limited.
- 46a.
Page 1: “400 genes showed significant evidence of positive selection”. Page 4: “differences in numbers of identified PSGs almost certainly reflect differences in power more than differences in the prevalence of selection.” Out of about 16,500 human genes with high-confidence matches in at least two other species. ↩
- 46a.
- 47.
Lewontin, Richard C. The Genetic Basis of Evolutionary Change. Columbia University Press. 1974. Lewontin's book-length review of how much genetic variation natural populations hold and what theory can explain it, written after his protein-gel surveys showed variation was far higher than expected. Chapter 5, "The Paradox of Variation," is the origin of the paradox named after him. Mirror: Archive.org
- 47a.
Page 208: “The observed range of heterozygosities over all the species listed in table 22 lies in the sensitive region, between 0.056 and 0.184. This range corresponds to values of Nμ between 0.015 and 0.057. Since there is no reason to suppose that mutation rate has been specially adjusted in evolution to be the reciprocal of population size for higher organisms, we are required to believe that higher organisms including man, mouse, Drosophila and the horseshoe crab all have population sizes within a factor of 4 of each other.” ↩
- 47b.
Page 208: “Despite the unchallengeable agreement between the heterozygosity observed in any given case and the predictions of neutral theory, the application of neutral theory to every case of heterozygosity turns out to lead to an absurd result.” ↩
- 47a.
- 48.
Leffler, Ellen M. et al. “Revisiting an old riddle: what determines genetic diversity levels within species?” PLoS Biology 10(9):e1001388. 2012. A compilation of published DNA diversity for 167 species and a review of what could explain its narrow range. It is the source of the "0.01% to 8%" span most later papers quote, and it concludes that the explanation is still open.
- 48a.
Page 4: “Among the 167 species that met our criteria for the autosomes, nucleotide diversity p ranges from 0.01% per base pair in Lynx lynx to 8.01% in Ciona savignyi, a span of only 800-fold”. "p" is the text layer's rendering of π. ↩
- 48a.
- 49.
Maynard Smith, John, and John Haigh. “The hitch-hiking effect of a favourable gene.” Genetical Research 23(1):23–35. 1974. The theory paper that described genetic hitchhiking: a favourable mutation sweeping to fixation lowers variation at linked neutral sites. Its authors offered it as a possible rescue of the neutral theory from the uniformity of variation across species.
- 49a.
Page 34: “Thus it seems that the uniformity of H between species is powerful evidence against the view that the observed polymorphisms are in the main selectively neutral. The investigation in §§ 2—4 can therefore be regarded as a last ditch attempt to save the neutral mutation theory by showing that there is another process which can account for the uniformity of H between species. The attempt has been partially successful.” ↩
- 49a.
- 50.
Begun, David J., and Charles F. Aquadro. “Levels of naturally occurring DNA polymorphism correlate with recombination rates in D. melanogaster.” Nature 356:519–520. 1992. The first genome-wide evidence, from 20 gene regions in fruit flies, that diversity is lower where recombination is rarer while divergence between species is not. It is the pattern every later linked-selection study builds on.
- 50a.
Page 519: “Here we show that for 20 gene regions from across the genome, the amount of nucleotide diversity in natural populations of D. melanogaster is positively correlated with the regional rate of recombination.” ↩
- 50a.
- 51.
Corbett-Detig, Russell B. et al. “Natural selection constrains neutral diversity across a wide range of species.” PLoS Biology 13(4):e1002112. 2015. A fit of background-selection and hitchhiking models to diversity and recombination maps in 40 sexual species. It found linked selection removes more neutral diversity in species with larger ranges and smaller bodies, a proxy for larger populations, but it could not measure census sizes directly.
- 51a.
Page 1: “This implies that natural selection may provide an explanation for this longstanding paradox of population genetics.” ↩
- 51a.
- 52.
Charlesworth, Brian. “The effects of deleterious mutations on evolution at linked sites.” Genetics 190(1):5–22. 2012. A review of background selection theory by one of its authors, with rough applications to human and fruit-fly chromosomes.
- 52a.
Page 11: “a fairly trivial effect”. For a centrally located site on a typical human autosome he estimates B ≈ 0.94, the share of neutral diversity left after background selection, and compares it with an estimate that hitchhiking reduced average human diversity by 6%. ↩
- 52a.
- 53.
McVicker, Graham et al. “Widespread genomic signatures of natural selection in hominid evolution.” PLoS Genetics 5(5):e1000471. 2009. A background-selection map of the human genome, fitted to human diversity and ape divergence around conserved sequence. It gives one of the larger estimates of how much diversity selection has removed in humans, and could not separate background selection from hitchhiking.
- 53a.
Page 1: “a surprisingly large genome average diversity reduction due to selection of 19–26% on the autosomes and 12–40% on the X chromosome.” ↩
- 53a.
- 54.
Comeron, Josep M. “Background selection as baseline for nucleotide variation across the Drosophila genome.” PLoS Genetics 10(6):e1004434. 2014. A genome-wide prediction of background selection in fruit flies built from gene annotations and recombination maps, without fitting the diversity data, then compared with observed diversity.
- 54a.
Page 4: “median B of 0.559 and 0.736 for autosomes and the X chromosome, respectively”. ↩
- 54a.
- 55.
Elyashiv, Eyal et al. “A genomic map of the effects of linked selection in Drosophila.” PLoS Genetics 12(8):e1006130. 2016. A joint model of background selection and sweeps fitted to fruit-fly diversity, which gives the largest published estimate of how much neutral diversity linked selection removes, with the authors' own caveat about its fitted parameters.
- 55a.
Page 17: “Doing so indicates an average reduction of 77%-89% in neutral diversity levels genome-wide (excluding low- recombination regions for which the reduction should be greater).” ↩
- 55b.
Page 15: “The strong selection coefficient is also associated with unreasonably high estimates of the deleterious mutation rate, which far exceed direct estimates of the total mutation rate (by 4-9-fold in exons and UTRs”. ↩
- 55a.
- 56.
Pouyet, Fanny et al. “Background selection and biased gene conversion affect more than 95% of the human genome and bias demographic inferences.” eLife 7:e36317. 2018. An analysis of high-quality human genomes that sorts each variant by how strongly background selection and GC-biased gene conversion, a repair bias tied to recombination, should affect it. Almost none escapes both, which matters for any method that treats noncoding variants as a neutral baseline.
- 57.
Romiguier, Jonathan et al. “Comparative population genomics in animals uncovers the determinants of genetic diversity.” Nature 515:261–263. 2014. Transcriptome sequencing of two to ten individuals in each of 76 non-model animal species, from Galtier's group, to see what predicts a species' genetic diversity. Life history predicted it; geographic range did not.
- 57a.
Page 261: “The distribution of genetic diversity between species revealed no detectable influence of geographic range or invasive status but was accurately predicted by key species traits related to parental investment: long-lived or low-fecundity species with brooding ability were genetically less diverse than short-lived or highly fecund ones.” ↩
- 57a.
- 58.
Ponting, Chris P., and Ross C. Hardison. “What fraction of the human genome is functional?” Genome Research 21(11):1769–1776. 2011. A review of the evolutionary and experimental estimates of the functional fraction of the human genome, which argues that turnover of constrained sequence makes cross-species estimates too low.
- 58a.
Section Abstract: “Approaches that take account of turnover estimate that the steady-state value of α sel lies between 10% and 15%.” ↩
- 58a.
- 59.
Sanford, John C. Genetic Entropy. FMS Publications. 4th edition. 2014. The fourth edition of the book that founded the young-earth argument that human genomes are degenerating because slightly harmful mutations escape selection. Its argument depends on Kimura's and Ohta's nearly neutral class of mutations.
- 59a.
Page 21: “In just a few years, many geneticists have shifted from believing that less than 3% of the total genome is functional to believing that much more than 30% is functional, and most recently evidence suggests 80% or more is functional (The ENCODE Project Consortium, 2012). As the functional genome expands, the likelihood of neutral mutations shrinks.” ↩
- 59b.
Page 21: “While most sites are probably nearly-neutral, very few, if any, should be absolutely neutral.” ↩ ↩
- 59c.
Page 23: “Given this distribution of mutations, one would naturally ask, “How can theorists possibly explain evolutionary progress?” It is done as follows: everything in the “near-neutral box” is redefined as being completely neutral, and is thereby dismissed. It is then assumed that the mutations to the left of the near-neutral box can be entirely eliminated using natural selection.” ↩
- 59d.
Page 199: “Paradoxically, his theory led him to believe that most mutations are unselectable, and therefore most genetic sequences must be irrelevant, and most evolutionary change must be independent of selection. Because he was totally committed to the Primary Axiom, Kimura apparently never considered the possibility that his cost arguments could most rationally be used to argue against the Primary Axiom itself.” "Primary Axiom" is Sanford's name for the claim that mutation and selection produced the genome. ↩ ↩
- 59a.
- 60.
Carter, Robert W. “The Neutral Model of evolution and recent African origins.” Journal of Creation 23(1):70–77. 2009. A young-earth critique, by a geneticist on the staff of Creation Ministries International, of the "recent African origin" model of human ancestry. It argues that the model rests on the standard neutral model and that each of that model's assumptions is open to question.
- 60a.
Page 71: “Kimura saw Haldane’s Dilemma as a serious problem and listed it as his main reason for proposing the Neutral Theory of evolution.” ↩
- 60a.
- 61.
Sanford, John et al. “The waiting time problem in a model hominin population.” Theoretical Biology and Medical Modelling 12:18. 2015. Simulations by a team led by the Cornell plant geneticist who wrote Genetic Entropy, estimating how long a human-sized population waits for specific strings of mutations to arise and spread. Its introduction gives the young-earth account of why Kimura proposed the neutral theory.
- 61a.
Page 2: “Kimura [3] agreed with Haldane’s assessment of the problem, and this motivated him to develop his neutral model of evolution.” ↩
- 61a.
- 62.
ReMine, Walter James. The Biotic Message: Evolution versus Message Theory. St. Paul Science. 1993. A book by a creationist engineer arguing that life was designed to carry a message against evolution, best known for its modern statement of Haldane's dilemma. Chapter 9 treats the neutral theory as an attempt to escape that dilemma. Mirror: Archive.org
- 62a.
Page 237: “These arguments were never ‘for’ the neutral theory. They are arguments against selective evolution, severely limiting its speed.” ↩
- 62a.
- 63.
Carter, Robert W. “The slow, painful death of junk DNA.” Creation Ministries International. May 20, 2009. A popular article arguing that junk DNA was a theory-driven idea, required by Haldane's dilemma and the neutral theory, and that discoveries of function are refuting it.
- 63a.
“Junk DNA is a necessary mathematical extrapolation. It was invented to solve a theoretical evolutionary dilemma.” “If “Haldane’s Dilemma” was correct, the majority of DNA must be non-functional.” A reference mark after "Haldane’s Dilemma" is left out. ↩
- 63a.
- 64.
Tomkins, Jeffrey P., and Jerry Bergman. “Evolutionary molecular genetic clocks—a perpetual exercise in futility and failure.” Journal of Creation 29(2):26–35. 2015. A young-earth review arguing that molecular clocks have failed, that the clock rests on the neutral theory, and that evidence of widespread function undermines both.
- 64a.
Page 33: “Furthermore, because the molecular clock idea is directly tied to the neutral model theory of evolution, recent discoveries in full codon utility and pervasive genomewide biochemical functionality, present serious obstacles to the evolutionary necessity of a large fraction of the genome being ‘junk’.” ↩
- 64a.
- 65.
Carter, Robert W. “Robert Carter gets everything wrong?” Creation Ministries International. Jun 7, 2021. Carter's reply to a video critique of his work, including a section on Haldane's dilemma.
- 65a.
“The attempted rescue mechanism of neutral evolution theory, and the associated belief in ‘junk DNA’, failed as an escape mechanism.” ↩
- 65a.
- 66.
Lukeš, Julius et al. “How a neutral evolutionary ratchet can build cellular complexity.” IUBMB Life 63(7):528–537. 2011. A perspective by five cell and evolutionary biologists describing "constructive neutral evolution," a proposed route by which complex machinery could grow without positive selection. It offers scenarios, not measured rates.
- 67.
Jeanson, Nathaniel T. Replacing Darwin: The New Origin of Species. Master Books. 2017. A young-earth model of speciation by a Harvard-trained cell biologist at Answers in Genesis, in which differences between families were created and differences within them arose after the Flood. It makes explicit predictions about which genetic differences should be functional. Mirror: Archive.org
- 67a.
Page 231: “In other words, creationists predict that the majority — if not the vast majority — of nuclear DNA differences among species within families (for vertebrates) or within subgenera/genera (invertebrates and fungi) are functional participants in the biology of each species, not functionally neutral.” ↩
- 67a.
- 68.
Jeanson, Nathaniel T., and Jason Lisle. “On the origin of eukaryotic species' genotypic and phenotypic diversity.” Answers Research Journal 9:81–122. 2016. A young-earth model in which most present genetic diversity was created as heterozygosity in the founding pairs and then sorted after the Flood, tested against mutation rates and population growth curves.
- 68a.
Page 115: “genetic drift and other processes likely have diluted this concentrated heterozygosity to lower levels.” ↩
- 68a.
- 69.
Kimura, Motoo, and Tomoko Ohta. “The average number of generations until fixation of a mutant gene in a finite population.” Genetics 61(3):763–771. Mar 1969. The paper that derived how long a mutation takes to spread through a finite population, for neutral and selected alleles. Its equation 14 gives the mean time to fixation of a neutral allele from any starting frequency.
- 70.
Galtier, Nicolas. “S1 Table: species, sample sizes and estimates of α and ωa for 44 animal species pairs.” PLoS Genetics 12(1):e1005774. Jan 11, 2016. The per-species table behind Galtier 2016, used here for the median and counts. ↩