Bias in the introduction of variation explained

Bias in the introduction of variation ("arrival bias") is a theory in the domain of evolutionary biology that asserts biases in the introduction of heritable variation are reflected in the outcome of evolution. It is relevant to topics in molecular evolution, evo-devo, and self-organization. In the context of this theory, "introduction" ("origination") is a technical term for events that shift an allele frequency upward from zero (mutation is the genetic process that converts one allele to another, whereas introduction is the population genetic process that adds to the set of alleles in a population with non-zero frequencies). Formal models demonstrate that when an evolutionary process depends on introduction events, mutational and developmental biases in the generation of variation may influence the course of evolution by a first come, first served effect, so that evolution reflects the arrival of the likelier, not just the survival of the fitter.Whereas mutational explanations for evolutionary patterns are typically assumed to imply or require neutral evolution, the theory of arrival biases distinctively predicts the possibility of mutation-biased adaptation. Direct evidence for the theory comes from laboratory studies showing that adaptive changes are systematically enriched for mutationally likely types of changes.Retrospective analyses of natural cases of adaptation also provide support for the theory.This theory is notable as an example of contemporary structuralist thinking, contrasting with a classical functionalist view in which the course of evolution is determined by natural selection (see [1]).

History

The theory of biases in the introduction process as a cause of orientation or direction in evolution has been explained as the convergence of two threads. The first, from theoretical population genetics, is the explicit recognition by theoreticians (toward the end of the 20th century) that a correct treatment of evolutionary dynamics requires a rate-dependent process of introduction (origination) missing from classical treatments of evolution as a process of shifting frequencies of available alleles. This recognition is evident in the emergence of origin-fixation models that depict evolution as a 2-step process of origination and fixation (by drift or selection), with a rate specified by multiplying a rate of introduction (based on the mutation rate) with a probability of fixation (based on the fitness effect). Origin-fixation models appeared in the midst of the molecular revolution, a half-century after the origins of theoretical population genetics: they were soon widely applied in neutral models for rates and patterns of molecular evolution; their use in models of molecular adaptation was popularized in the 1990s; by 2014 they were described as a major branch of formal theory.

The second thread is a long history of attempts to establish the thesis that mutation and development exert a dispositional influence on evolution by presenting options for subsequent functional evaluation, i.e., acting in a manner that is logically prior to selection. Many evolutionary thinkers have proposed some form of this idea. In the early 20th-century, authors such as Eimer or Cope held that development constrains or channels evolution so strongly that the effect of selection is of secondary importance.[2] [3] Early geneticists such as Morgan and Punnett proposed that common parallelisms (e.g., involving melanism or albinism) may reflect mutationally likely changes.[4] Expanding on Vavilov's (1922) [5] exploration of this theme, Spurway (1949)[6] wrote that "the mutation spectrum of a group may be more important than many of its morphological or physiological features."

Similar thinking featured in the emergence of evo-devo, e.g., Alberch (1980) suggests that "in evolution, selection may decide the winner of a given game but development non-randomly defines the players" (p. 665)[7] (see also [8]). Thomson (1985),[9] reviewing multiple volumes addressing the new developmentalist thinking— a book by Raff and Kaufman (1983)[10] and conference volumes edited by Bonner (1982)[11]  and Goodwin, et al (1983) [12] —wrote that "The whole thrust of the developmentalist approach to evolution is to explore the possibility that asymmetries in the introduction of variation at the focal level of individual phenotypes, arising from the inherent properties of developing systems, constitutes a powerful source of causation in evolutionary change" (p. 222). Likewise, the paleontologists Elisabeth Vrba and Niles Eldredge summarized this new developmentalist thinking by saying that "bias in the introduction of phenotypic variation may be more important to directional phenotypic evolution than sorting by selection." [13]

However, the notion of a developmental influence on evolution was rejected by Mayr and others such as Maynard Smith ("If we are to understand evolution, we must remember that it is a process which occurs in populations, not in individuals.")[14] and Bruce Wallace ("problems concerned with the orderly development of the individual are unrelated to those of the evolution of organisms through time"),[15] as being inconsistent with accepted concepts of causation. This conflict between evo-devo and neo-Darwinism is the focus of a book-length treatment by philosopher Ron Amundson [16] (see also Scholl and Pigliucci, 2015 [17]). In the theory of evolution as shifting gene frequencies that prevailed at the time, evolutionary causes are "forces" that act as mass pressures (i.e., the aggregate effects of countless individual events) shifting allele frequencies (see Ch. 4 of), thus development did not qualify as an evolutionary cause. A widely cited 1985 commentary on "developmental constraints" advocated the importance of developmental influences, but did not anchor this claim with a theory of causation, a deficiency noted by critics, e.g., Reeve and Sherman (1993) defended the adaptationist program (against the developmentalists and the famous critique of adaptationism by Gould and Lewontin), arguing that the "developmental constraints" argument simply restates the idea that development shapes variation, without explaining how such preferences prevail against the pressure of selection. Mayr (1994)[18] insisted that developmentalist thinking was "hopelessly mixed up" because development is a proximate cause and not an evolutionary one. In this way, developmentalist thinking was received in the 1980s and 1990s as speculation without a rigorous grounding in causal theories, an attitude that persists (e.g., Lynch, 2007 [19]).

In response to these rebukes, developmentalists concluded that population genetics cannot provide a complete account of evolutionary causation:[20] instead, a dry statistical account of changes in gene frequencies from population genetics must be supplemented with a wet biological account of changes in developmental-genetic organization (called "lineage explanation" in[21]). The beliefs that (1) developmental biology was never integrated into the Modern Synthesis and (2) population genetics must be supplemented with alternative narratives of developmental causation, are now widely repeated in the evo-devo literature and are given explicitly as motivations for reform via an Extended Evolutionary Synthesis.[22]

The proposal to recognize the introduction process formally as an evolutionary cause provides a different resolution to this conflict. Under this proposal, the key to understanding the structuralist thesis of the developmental biologists was a previously missing population-genetic theory for the consequences of biases in introduction. The authors criticized classical reasoning for framing the efficacy of variational tendencies as a question of evolution by mutation pressure, i.e., the transformation of populations by recurrent mutation. They argued that, if generative biases are important, this cannot be because they out-compete selection as forces under the shifting-gene-frequencies theory, but because they act prior to selection, via introduction. Thus the theory of arrival biases proposes that the generative dispositions of a developmental-genetic system (i.e., its tendencies to respond to genetic perturbation in preferential ways) shape evolution by mediating biases in introduction. The theory, which applies to both mutational and developmental biases, addresses how such preferences can be effective in shaping the course of evolution even while strong selection is at work.

Systematic evidence for predicted effects of introduction biases first began to appear from experimental studies of adaptation in bacteria and viruses.[23] Since 2017, this support has widened to include systematic quantitative results from laboratory adaptation, and similar but less extensive results from the retrospective analysis of natural adaptations traced to the molecular level (see below). The empirical case that biases in mutation shape adaptation is considered to be established for practical purposes such as evolutionary forecasting (e.g., [24]).However, the implications of the theory have not been tested critically in regard to morphological and behavioral traits in animals and plants that are the traditional targets of evolutionary theorizing (see Ch. 9 of). Thus, the relevance of the theory to molecular adaptation has been established, but the significance for evo-devo remains unclear. The theory sometimes appears associated with calls for reform from advocates of evo-devo (e.g.,), though it has not yet appeared in textbooks or in broad treatments of challenges in evolutionary biology (e.g.,[25] [26]).

Simple model

The kind of dual causation proposed by the theory has been explained with the analogy of "Climbing Mount Probable." Imagine a robot on a rugged mountain landscape, climbing by a stochastic 2-step process of proposal and acceptance. In the proposal step, the robot reaches out with its limbs to sample various hand-holds, and in the acceptance step, the robot commits and shifts its position. If the acceptance step is biased to favor higher hand-holds, the climber will ascend. But one also may imagine a bias in the proposal step, e.g., the robot may sample more hand-holds on the left than on the right. Then the dual proposal-acceptance process will show both an upward bias due to a bias in acceptance, and a leftward bias due to a bias in proposal. If the landscape is rugged, the ascent will end on a local peak that (due to the proposal bias) will tend to be to the left of the starting point. On a perfectly smooth landscape, the climber will simply spiral to the left until the single global peak is reached.In either case, the trajectory of the climber is subject to a dual bias. These two biases are not pressures competing to determine an allele frequency: they act at different steps, along non-identical dimensions.

The dual effect predicted by the theory was demonstrated originally with a population-genetic model of a 1-step adaptive walk with 2 options, i.e., the climber faces two upward choices, one with a higher selection coefficient and the other with a higher mutation rate. A key feature of model is that neither of the alternatives is present in the initial population: they must be introduced. In simulated adaptation under this model, the population frequently reaches fixation for the mutationally favored allele, even though it is not the most fit option. The form of the model is agnostic with respect to whether the biases are mutational or developmental. Subsequent theoretical work (below) has generalized on the theory of one-step walks, and also considered longer-term adaptive walks on complex fitness landscapes. The general implication for parallel evolution is that biases in introduction may contribute strongly to parallelism. The general implication for the directionality and repeatability of adaptive walks is simply that some paths are more evolutionarily favorable due to being mutationally favorable. The general implication for the long-term predictability of outcomes, e.g., particular phenotypes, is that some phenotypes are more findable than others due to mutational effects, and such effects may strongly shape the distribution of evolved phenotypes.

The application of the theory to problems in evo-devo and self-organization relies formally on the concept of a genotype-phenotype (GP) map. The genetic code, for example, is a GP map that induces asymmetries in mutationally accessible phenotypes. Consider evolution from the Met (amino acid) phenotype encoded by the ATG (codon) genotype. A phenotypic shift from Met to Val requires an ATG to GTG mutation; a shift from Met to Leu can occur by 2 different mutations (ATG to CTG or TTG); a shift from Met to Ile can occur by 3 different mutations (to ATT, ATC, or ATA). If each type of genetic mutation has the same rate, i.e., with no mutation bias per se, the GP map induces 3 different rates of introduction of the alternative phenotypes Val, Leu and Ile. Due to this bias in introduction, evolution from Met to Ile is favored, and this is not due to a mutational bias (in the sense of a bias reflecting the mechanisms of mutagenesis), but rather an asymmetric mapping of phenotypes to mutationally accessible genotypes.

Results of theoretical modeling

One-step adaptive walks

As noted above, in the simplest case of the "Climbing Mount Probable" effect, one may consider a climber facing just two fixed choices: up and to the left, or up and to the right. This case is modeled using simulations by,[27] and is given a more complete treatment by [28] In general, the limiting behavior of evolution as the supply of new mutations becomes arbitrarily small, i.e., as

\muN0

, is called "origin-fixation" dynamics .[29] The origin-fixation approximation for choosing between the left and right options

j

and

k

(respectively) in the Yampolsky-Stoltzfus model is given by the following:

where

\muij

(or

\muik

) and

sij

(or

sik

) are the mutation rate and selection coefficient for the left (or right) alternative, and assuming that the probability of fixation

\pis2s

. In the Yampolsky-Stoltzfus model, this approximation is good for

\muN<1/10

.

For 1-step walks under origin-fixation conditions, the behavior given by Eqn generalizes from 2 to many alternatives. For instance, Cano, et al. (2022)[30] consider a model gene with many different beneficial mutations, and under low mutation supply, the effects of mutation bias are proportional on the spectrum of adaptive changes.

When

\muN

is not very small, different beneficial alleles may be present simultaneously, competing and slowing down adaptation, an effect known as clonal interference. Clonal interference reduces the effect of mutation bias in models of evolution in finite genetic spaces: alleles favored by mutation still tend to arrive sooner, but before they reach fixation, later-arising alleles that are more beneficial can out-compete them, enhancing the effect of fitness differences. Under the most extreme condition when all possible beneficial alleles are reliably present in a large population, the most fit allele wins deterministically and there is no room for an effect of mutation bias.[30] [31] Stated differently, when all the beneficial alleles are present and selection determines the winner, the chance of success is 1 for the most fit allele, and 0 for all other alleles. Thus, in a gene model with a finite set of beneficial mutations, the influence of mutation bias is expected to be strongest when

\muN<<1

but to fall off as

\muN

becomes large.

The influence of mutation under varying degrees of clonal interference can be quantified precisely using the regression method of Cano, et al (2022).[30] Suppose that the expected number of changes of a given class of mutational changes defined by starting and ending states is directly proportional to the product of (1) the frequency

f(c)

of the starting state and (2) the mutation rate

\mu(c,a)

raised to the power of

\beta

, that is,

Taking the logarithm of this equation gives

where

\alpha

is the logarithm of the constant of proportionality. Thus, when

\beta

is unknown, it may be estimated as the coefficient for the regression of log(counts) on log(expected counts). Simulations of a gene model (figure at right from [30]) show a range from

\beta0

under high mutation supply to

\beta1

when the mutation supply is low. While this approach was developed to assess how the mutation spectrum influenced adaptive missense changes (defined by a starting codon and an ending amino acid), the equation reflects a generic framework applicable to any mutationally defined classes of change.

Note that these considerations apply to finite genetic spaces. In an infinite genetic space, clonal interference still slows down the rate of adaptation due to competition, but it does not prevent an effect of mutation bias because there are always mutationally favored alternative alleles among the most-fit class of alleles.[32] [33]

Contribution of mutation to parallelism

In general, if there is some set of possible steps each with a probability

pi

, then the chance of parallelism is given by summing the squares,

Ppara=\sumi

2
p
i
. It follows from the definition of the variance

V

or the coefficient of variation

cv

, that (see Box 2 of or Ch. 8 of)

That is, parallelism is increased by anything that decreases the number of choices or increases the heterogeneity in their chances (as measured by

V

or

cv

). This result validates the intuition of Shull[34] that "It strains one's faith in the laws of chance to imagine that identical changes should crop out again and again if the possibilities are endless and the probabilities equal" (p. 448). To the extent that heterogeneity in

pi

reflects heterogeneity in mutational chances, mutation contributes to parallelism.

In particular, for the case of origin-fixation dynamics, each value of

pi

is a product of a mutational origin term and a fixation term, so that heterogeneity in either contributes similarly to the chances of parallelism, and it is possible to partition effects of mutation and selection in accounting for the repeatability of evolution. Under origin-fixation conditions, and assuming

\pi(s,N)2s

, it follows that [35]

where

cv(\boldsymbol{s})

and

cv(\boldsymbol{\mu})

are coefficients of variation for vectors of selection coefficients and mutation rates, respectively. Numeric examples in Box 2 of suggest that mutation sometimes contributes more to parallelism than selection, although the authors note that

n

in the denominator above confounds effects of mutation and selection in a hidden way (because, in practice,

n

reflects the set of paths that are sufficiently favored by selection and sufficiently mutationally likely to be observed).

Longer-term effects: trends, navigability, and findability

For a systematic view of long-term effects of evolution in discrete genotypic space, consider the 4 perspectives below, focusing on the influence of a mutation spectrum (characteristic for some evolving system) on various ways of defining the chances of evolution (following the treatment by):

X

with access to other nearby points in genotype space. From

X

, there are 0 or more upward steps or paths that differ in mutational favorability (as a function of the mutation spectrum) and fitness benefits. The evolvability-from-

X

is a function of this set of steps. Under simple conditions, each step has a probability

pi

and the repeatability of evolution is derived by squaring the

pi

values.

Theoretical results relating to each of these perspectives are available. For instance, in a simulation of adaptive walks of protein-coding genes in the context of an abstract NK landscape,[36] the effect of a GC-AT mutation bias is to alter the protein sequence composition in a manner qualitatively consistent with the analogy of Climbing Mount Probable (above). Each adaptive walk begins with a random sequence and ends on some local peak; the direction of the walk and the final peak depend on the mutation bias.[36] For instance, adaptive walks under a mutation bias toward GC result in proteins that have more of the amino acids with GC-rich codons (Gly, Ala, Arg, Pro), and likewise, adaptive walks under AT bias result in proteins with more of the amino acids with AT-rich codons (Phe, Tyr, Met, Ile, Asn, Lys). On a rough landscape, the initial effect is similar, but the adaptive walks are shorter.[36] That is, the mutation bias imposes a preference (on the adaptive walks) for steps, paths, and local peaks that are enriched in outcomes favored by the mutation bias. This illustrates the concept of a directional trend in which the system moves cumulatively in a particular direction along an axis of composition.

The influence of transition-transversion bias has been explored using empirical fitness landscapes for transcription factor binding sites.[37] Each landscape is based on generating thousands of different 8-nucleotide fragments and measuring how well they bind to a particular transcription factor. Each peak on each landscape is accessible by some set of paths made of steps that are nucleotide changes, each one being either a transition or a transversion. Among all possible genetic changes, the ratio of transitions to transversions is 1:2. However, the collection of paths leading to a given peak (on a given empirical landscape) has a specific transition-transversion composition that may differ from 1:2. Likewise, any evolving system has a particular transition-transversion bias in mutation. The more closely the mutation bias (of the evolving system) matches the composition bias (of the landscape), the more likely that the evolving system will find the peak.[37] Thus, for a given evolving system with its characteristic transition-transversion bias, some landscapes are more navigable than others. Navigability is maximized when the mutation bias of the evolving system matches the composition bias of the landscape.[37]

Finally, rather than organizing genotypes by fitness (in terms of peaks, upward paths, and collections of paths leading to a peak), we can organize genotypes by phenotype using a genotype–phenotype map. A given phenotype identifies a network in genotype-space including all of the genotypes with that phenotype. An evolving system may diffuse neutrally within the network of genotypes with the same phenotype, but conversions between phenotypes are assumed to be non-neutral. Each phenotypically defined network has a findability that is, as a first approximation, a function of the number of genotypes in the network.For instance, using the canonical genetic code as a genotype-phenotype map, the phenotype Leucine has 6 codons whereas Tryptophan has 1: Leucine is more findable because there are more mutational paths from non-Leucine genotypes. This idea can be applied to the way that RNA folds (considered as phenotypes) map to RNA sequences. For instance, evolutionary simulations show that the RNA folds with more sequences are more findable, and this is due to the way that they are over-sampled by mutation.[38] A similar point has been made in regard to substructures of regulatory networks [39] (see also [40]).

The above results apply, as before, to finite spaces. In infinite spaces, the set of remaining beneficial mutations to be explored is infinite and includes an infinite supply of mutationally favored and mutationally disfavored options. Therefore evolution in infinite spaces can continue forever in the mutationally favored direction with no diminution of the mutational effect that applies in the short-term, e.g., one could consider Eqn for such an infinite space. The model of Gomez, et al (2020) allows unlimited adaptation via two traits, one with a higher rate of beneficial mutation, and the other with larger selective benefits. In this model, mutation bias continues to be important in long-term evolution even when mutation supply is very high.

Distinctive implications

The theory of biases in the introduction process as a cause of orientation or direction in evolution may be contrasted with other theories that have been used by evolutionary biologists to reason about the role of variation in evolution:

\bf{G}

matrix (standing variation) is a source of dimensional but not directional asymmetry, depending on the amount of variation available along any given dimension in trait-space.[43] [44]

Relative to these theories, the theory of arrival biases has distinctive implications, some of which are supported empirically as described below, e.g., the most frequent outcome of an adaptive process such as the emergence of antibiotic resistance is not necessarily the most beneficial, but is often a moderately beneficial outcome favored by a high rate of mutational origin. Likewise, the theory implies that evolution can have directions that are not adaptive, or tendencies that are not optimal, an implication one commentator on Arthur's bookfound "disturbing".[45] This theory is defined, not by any particular problem, taxon, level of organization, or field of study, but by a mechanism defined at the level of population genetics, namely the ability of biases in introduction to impose biases on evolution. Some implications are as follows (see).

Effects do not require neutrality or high mutation rates. In contrast to the theory of evolution by mutation pressure explored (and rejected) by Haldane and others, variational dispositions under the theory of arrival biases do not depend on neutral evolution and do not require high mutation rates.

Graduated biases can have graduated effects. In contrast to what is implied by the language of "constraints" or "limits" employed in historic appeals to internal sources of direction in evolution, the theory of arrival biases is not deterministic and does not require an absolute distinction between possible and impossible forms. Instead, the theory is probabilistic, and graduated biases can have graduated effects.

Regime-dependency with regard to population genetics. Under the theory, variation biases do not have a guaranteed effect independent of the details of population-genetics. The influence of mutation biases reaches a maximum (proportional influence) under origin-fixation conditions and can disappear almost entirely under high levels of mutation supply.

Parity in fixation biases and origination biases (under limiting conditions). In classical neo-Darwinian thinking, selection governs and shapes evolution, whereas variation plays a passive role of supplying materials. By contrast, under limiting origin-fixation conditions, the theory of arrival biases establishes a condition of parity such that (for instance) a 2-fold bias in fixation and a 2-fold bias in introduction both have the same 2-fold effect on the chances of evolution.[29]

Generality with regard to sources of variational bias. In the evolutionary literature, mutation biases, developmental "constraints", and self-organization in the sense of findability are all treated as separate topics. Under the theory of arrival biases, these are all manifestations of the same kind of population-genetic mechanism in which biases in the introduction of variants imposes biases on evolution. Any short-term bias is either a mutational bias in the sense of a difference in rates for two fully specified genotypic conversions, or it can be treated as a scheme of differential phenotypic aggregation over genotypes (see Box 2 of).

In addition to these direct implications, some more sophisticated or indirect implications have emerged in the literature.

Non-causal associations induced by mutation and selection. Due to a dual dependence on mutation and selection, the distribution of adaptive changes may show non-causal associations of mutation rates and selection coefficients, somewhat akin to Berkson's paradox, as suggested in Ch. 8 of. and developed in more detail by Gitschlag, et al (2023).[46]

Conditions for composition and decomposition of causes. Under limiting origin-fixation conditions, the chances of evolution reflect two factors multiplied together, representing biases in introduction and biases in fixation, as in Eqn . Thus, conditions exist under which it is possible to quantify and directly compare the dispositional influences of mutation and selection. This approach has already been used in a few empirical arguments addressed below.[47] [48] [49] (Box 2 of).

Biased depletion of the spectrum of beneficial mutations. In any case of a system adapting via mutation and selection, there is some set of possible beneficial mutations, characterized by a distribution of selection coefficients and mutation rates. As adaptation occurs in a mutation-biased manner, this spectrum of possible beneficial mutations is depleted in a biased way. The theory for this depletion [50] is relevant to experimental work showing that "shifts in mutation spectra enhance access to beneficial mutations".[51] That is, the experimentally observed favorability of shifts in mutation spectra depends on a pattern of biased depletion of beneficial mutations that is itself a sign of mutation-biased adaptation.

Evidence

Evidence for the theory has been summarized recently, e.g., Gomez, et al (2020) [33] present a table listing 8 different studies providing evidence of effect of mutation bias on adaptation, and Ch. 9 of Mutation, Randomness and Evolution is devoted to empirical support for the theory (see also [32]). Biases in introduction are expected to influence evolution whether neutral or adaptive, but an effect on neutral evolution is not considered intuitively surprising or controversial, and so is not given much attention. Instead, accounts of evidence focus on mutation-biased adaptation, because this highlights how predictions of the theory clash with the classical conception of mutation as a weak pressure easily overcome by selection, per the "opposing pressures" argument of Fisher and Haldane.

Direct evidence of causation under controlled conditions

Direct evidence that the spectrum of mutation shapes the spectrum of adaptive changes comes from studies that manipulate the mutation spectrum directly. In one study,[52] resistance to cefotaxime was evolved repeatedly, using 3 strains of E. coli with different mutation spectra: wild-type, mutH and mutT. The spectrum of resistance mutations among the evolved strains showed the same patterns of spontaneous mutations as the parental strains. Specifically, the

A:TC:G

transversions favored by mutT (left block of bars) are highly enriched among resistant isolates from mutT parents (blue in the accompanying figure), and likewise, the resistant strains from mutH parents (red) tend to have the nucleotide transition mutations favored by mutH (center block of bars). Thus, changing the mutation spectrum changes the spectrum of adaptive changes in a corresponding manner.

Another study[53] showed that the AR2 strain of P. fluorescens adapted to the loss of motility overwhelmingly (> 95% of the time) by one specific change, an A289C change in the ntrB gene, while the Pf0-2x strain adapted via diverse changes in several genes. The pattern in AR2 derivatives was traced to a mutational hotspot. Because the hotspot behavior was associated mainly with synonymous differences between the two strains, the experimenters were able to use genetic engineering to remove the hotspot from AR2 and add it to Pf0-2x, without changing the encoded amino acid sequence. This reversed the qualitative pattern of outcomes, so that the modified AR2 (engineered to remove the hotspot) adapted via diverse changes, while the modified Pf0-2x with the engineered hotspot adapted via the A289C change 80% of the time.

Graduated effects

A different use of available evidence is to focus on the idea of graduated effects, which distinguishes the theory of arrival biases from the intuitive notion of "constraints" or "limits" on possible forms. In particular, one may set aside the dramatic effects associated with hotspots and mutator alleles, and consider the effects of ordinary quantitative biases in nucleotide mutations. A number of studies have established that modest several-fold biases in mutation can have a several-fold effect on evolution, and some studies indicate a roughly proportional relation between mutation rates and the chances of an adaptive change.[54] [55] [56] [57] [58]

For instance, Sackman, et al (2017) studied parallel evolution in 4 related bacteriophages. In each case, they adapted 20 cultures in parallel, then sequenced a sample of the adapted culture to identify causative mutations. The results showed a strong preference for nucleotide transitions, 29:5 for paths (white vs. black bars in the figure at right) and 74:6 for events.

In a study of resistance to Rifampicin in Pseudomonas aeruginosa, MacLean, et al (2010)[59] measured selection coefficients and frequency of evolution for 35 resistance mutations in the rpoB (RNA polymerase) gene, and reported mutation rates for 11 of these. The mutation rates vary over a 30-fold range. The frequency with which a resistant variant appears in the set of 284 replicate cultures correlates strongly and roughly linearly with the mutation rate (figure at right). This is not explained by a correlation between selection coefficients and mutation rates, which are not correlated (see Ch. 9 of).

As explained above, the influence of the mutation spectrum on the spectrum of adaptive changes can be captured in a single parameter

\beta

, defined as a coefficient of binomial regression of observed counts on the expected counts from a mutational model.[30] Based on theoretical considerations, expected values of

\beta

range from 0 (no influence) to 1 (proportional influence). This method was applied by Cano, et al.[30] to 3 large data sets of adaptive changes, comparing a model based on independent measures of the mutation spectrum with adaptive changes previously identified in studies of (1) clinical antibiotic resistance of Mycobacterium tuberculosis, (2) laboratory adaptation in E. coli, and (3) laboratory adaptation of the yeast Saccharomyces cerevisiae to environmental stress. In each case,

\beta1

, indicating a roughly proportional influence of mutation bias. The authors report that this is not just due to the influence of transition-transversion bias, because

\beta1

applies both to transition-transversion bias and to the other aspects of the nucleotide mutation spectrum.[30]

Scope of applicability

A final use of available evidence is to consider the range of natural conditions under which the theory may be relevant. Whereas laboratory studies can be used to establish causation and assess effect-sizes, they do not provide direct guidance to where the theory applies in nature. As noted in, most studies of adaptation do not include a genetic analysis that identifies specific mutations, and in the rare cases in which an attempt is made to identify causative mutations, the results typically implicate only very small numbers of changes that are subject to questions of interpretation.[60] Therefore, the strategy followed in key studies[61] has been to focus on trusted cases of adaptation in which the proposed functional effects of putative adaptive mutations have been verified using techniques of genetics.

Payne, et al looked for an effect of transition bias among causative mutations for antibiotic resistance in clinically identified strains of Mycobacterium tuberculosis, which exhibits a strong mutation bias toward nucleotide transitions.[62] They compared the observed transition-transversion ratio to the 1:2 null expectation under the absence of mutation bias. Using two different curated databases, they found transition:transversion ratios of1755:1020 and 1771:900, i.e., enrichments 3.4-fold and 3.9-fold over the null, respectively. They also took advantage of the special case of Met-to-Ile replacements, which can take place by 1 transition (ATG to ATA) or 2 different transversions (ATG to ATT or ATC). This 1:2 ratio of possibilities again represents a null expectation for effects independent of mutation bias. In fact, the mutations in resistant isolates have transition-transversion ratios of 88:49 and 96:39 (for the 2 datasets), i.e., 3.6-fold and 4.9-fold above null expectations. This result cannot be due to selection at the amino acid level, because the changes are all Met to Ile. The significance of this result is not that mutation bias only works when the options are selectively indistinguishable: instead, the lesson is that the bias toward nucleotide transitions is roughly 4-fold both for the Met-to-Ile case, and for amino-acid-changing substitutions generally.

A much broader taxonomic scope is implicated in a meta-analysis of published studies of parallel adaptation in nature. In this study, the authors curated a data set covering 10 published cases of parallel adaptation traced to the molecular level, including well known cases involving spectral tuning, resistance to natural toxins such as cardiac glycosides[63] and tetrodotoxin, [64] [65] foregut fermentation, and so on. The results shown below (table) indicate a transition-transversion ratio of 132:99, a 2.7-fold enrichment relative to the null expectation of 1:2 (the ratio of paths, which is less sensitive to extreme values, is 27:28, a 2-fold enrichment). Thus, this study shows that a bias toward transitions is observed in well known cases of parallel adaptation in diverse taxa, including animals and plants.

Transitions and transversions involved in natural adaptive parallelisms (data from Table 3 of)
PhenotypeTaxonTargetTi counts (paths)Ti sumTv counts (paths)Tv sum
Insecticide resistanceInsectaRdl, Kdr, Ace2, 2, 5, 2, 3 (5)149, 2, 4 (3)15
Tetrodotoxin resistanceVertebrataNa channels2, 6, 3 (3)112, 2, 2, 3, 3 (5)12
Glycoside resistanceMetazoa+/+-ATPase4, 4, 2, 2 (4)127, 2, 2, 4 (4)15
Herbicide resistancePoaceaeACCase5, 2 (2)77, 2, 4, 5 (4)18
Altitude adaptationAves-hemoglobin4, 13 (2)172, 3, 2 (3)7
Trichromatic visionVertebrataOpsins2, 5 (2)76, 4, 2 (3)12
EcholocationMammaliaPrestin2, 2, 2 (3)63, 2 (2)5
Growth in RitonavirHIV1Protease25, 7, 9 (3)414 (1)4
Foregut fermentationVertebrataRibonucleases2, 4, 4 (3)10NA (0)0
Benzimidazole resistanceAscomycota-tubulin7 (1)75, 6 (2)11
(28)132(27)99

Finally, Storz, et al. analyzed changes in hemoglobin affinity associated with altitude adaptation in birds. Specifically, they studied the effect of CpG bias, an enhanced mutation rate at CpG sites due to effects of cytosine methylation on damage and repair, found widely in mammals and birds.[66] They assembled a data set consisting of 35 matched pairs of high- and low-altitude bird species. In each case, hemoglobins were evaluated for functional differences resulting in a higher oxygen affinity in the high-altitude species. The changes in affinity plausibly linked to adaptation implicated 10 different paths found a total of 22 different times. Six of the 10 paths involved CpG mutations, whereas only 1 would be expected by chance; and 10 out of 22 events involved CpG mutations, whereas only 2 would be expected by chance (both differences were significant). This enrichment of mutationally-likely genetic changes supports the theory of arrival biases and provides further evidence that predictable effects of mutation bias are important for understanding adaptation in nature.

Context in evolutionary thinking

The theory of arrival bias has been described as a cross-cutting theory[67] [68] [69] because it proposes a causal grounding (in population genetics) for diverse kinds of pre-existing claims for which a causal grounding is either unknown or mis-specified,

The context for applying the theory is illustrated in this figure (right). On the left are details of mutation and development that are responsible for tendencies in the generation of variation (varigenesis), i.e., tendencies prior to selection or drift. On the right are observable evolutionary patterns that might possibly be explained by these tendencies. The green arrow is some theory— the theory of arrival biases or some alternative theory— that specifies conditions of a cause-effect relationship linking variational tendencies to evolutionary tendencies. To apply a theory in this context is to generate evolutionary hypotheses or explanations that appeal to the internal details of mutation and development (on the left) to account for evolutionary patterns (on the right) via the conditions of causation specified by the theory. For instance, Darwin's comment that the laws of variation "bear no relation" to the structures built by selection would suggest that there are no conditions under which the internal details on the left account for the patterns on the right. The other theories all suggest that variational tendencies may influence evolution under some conditions. For instance, the theory of mutation pressure applies when mutation rates are high and unopposed by selection, thus it has a limited range of applications. The theory of evolutionary quantitative genetics can be applied very broadly to the evolution of quantitative characters, but the theory (as developed so far) does not suggest that mutation biases will have much impact.By contrast, the theory of arrival bias might apply broadly, and allows for a strong role for variational tendencies in shaping evolutionary tendencies.

Late arrival and non-obviousness

Though it seems intuitively obvious today, the theory did not emerge formally until 2001, e.g., as noted above, population geneticists did not propose the theory in the 1980s to answer an evo-devo challenge that literally called for recognizing biases in the introduction of variation. This late emergence has been attributed to a "blind spot" due to multiple factors, including a tradition of verbal arguments that minimize the role of mutation, a tendency to associate causation with processes that shift frequencies of variants rather than processes that create variants, and a formal argument from population genetics that doesn't extend to evolution from new mutations.

Specifically, when Haldane[74] and Fisher[75]   asked if tendencies of mutation could influence evolution, they framed this as a matter of the efficacy of mutation pressure (below), concluding that, because mutation rates are low, mutation is a weak force, only important in the special case of abnormally high mutation rates unopposed by selection. Their conception of evolutionary causation was modeled on selection, which operates by shifting frequencies of available alleles, and so they treated recurrent mutation in the same way. Their conclusion is correct for the case of evolution from standing variation.

More generally, in Modern Synthesis thinking, evolution was assumed to follow from a short-term process of shifting frequencies of available alleles.[76] [77] In this process, mutation is typically unimportant except when the focus is on low-frequency alleles maintained by deleterious mutation pressure (see Population genetics: Mutation), e.g., Edwards (1977) [78] addressed theoretical population genetics without considering mutation at all; Lewontin (1974) [79] stated that "There is virtually no qualitative or gross quantitative conclusion about the genetic structure of populations in deterministic theory that is sensitive to small values of migration, or any that depends on mutation rates" (p. 267).

The Haldane-Fisher "opposing pressures" argument was used repeatedly by leading thinkers to reject structuralist or internalist thinking (examples in or Ch. 6 of), e.g., Fisher (1930) stated that “The whole group of theories which ascribe to hypothetical physiological mechanisms, controlling the occurrence of mutations, a power of directing the course of evolution, must be set aside, once the blending theory of inheritance is abandoned." Seventy years later Gould (2002), citing Fisher (1930), wrote that “Since orthogenesis can only operate when mutation pressure becomes high enough to act as an agent of evolutionary change, empirical data on low mutation rates sound the death-knell of internalism.” (p. 510)[80]  

In this way, arguments from population genetics were used to reject, rather than support, speculative claims about the role of variational tendencies. The flaw in the Haldane-Fisher argument, pointed out in, is that it treats mutation only as a pressure on frequencies of existing alleles, not as a cause of the origin of new alleles. When alleles relevant to the outcome of evolution are absent initially, biases in introduction can impose strong biases on the outcome. Thus, the late appearance of this theory sheds light on how closely Modern Synthesis thinking was tied to the assumption of standing variation,[81] and to the forces theory. These commitments continue to echo in contemporary sources, e.g., in a US white paper endorsed by SSE, SMBE, ASN, ESA and other relevant professional societies, Futuyma, et al (2001)[82] state, as fact, that evolution is shifting gene frequencies, identifying the main causes of "evolution" (so defined) as selection and drift (figure).

However, toward the end of the 20th century, theoreticians began to note that long-term dynamics depend on events of mutational introduction not covered in classical theory.[83] [84] [85] In the recent literature, the assumption of evolution from standing variation is only rarely made explicit, e.g., [86] .[87] More commonly, evolution from standing variation is presented as an option to be considered together with evolution from new mutations. [88] [29] [89]

Relevance to contemporary issues

Parallelism and predictability. The application of the theory to parallelism is addressed above. The tendency for particular outcomes to recur in evolution is not merely a function of selection, but also reflects biases in introduction due to differential accessibility by mutation (or, for the case of phenotypes, by mutation and altered development). Recent reviews on prediction apply the theory to the role of mutation biases in contributing to the repeatability of evolution.

Partitioning causal responsibility for patterns to mutation and selection. In the case of origin-fixation dynamics, evolutionary dispositions can be attributed to a combination of mutation and selection, and it is possible, in principle, to untangle these contributions, as in an analysis of regulatory vs. structural effects in evolution or patterns of amino acid replacement in protein evolution.

Evo-devo, GP maps, and findability. The application of the theory of arrival biases to development and phenotypes is mediated by the concept of a genotype-phenotype map. A simple example of a biased induced by a GP map is shown at right. An evolving system diffuses neutrally within the genotypic network for its phenotype, and may occasionally jump to another phenotype. From the starting network of genotypes encoding phenotype P0, there are mutations leading to genotypic networks for P1 and P2. However, the number of mutations leading from the starting network to P2 is 4 times higher illustrating the idea that, for a given developmental-genetic system, some phenotypes are more mutationally accessible. This is not the same thing as a mutation bias per se (an asymmetry caused by the details of mutagenesis), but it can have the same effect in a population-genetic model. In this case, if all mutations happen at the same rate, the total rate of mutational introduction of P2 is 4 times higher than for P1: this bias can be mapped to the Yampolsky-Stoltzfus model and would have the same implications as a 4-fold mutation bias.

For short-term evolution, what matters is the distribution of immediately mutationally accessible phenotypes. In long-term evolution, however, one may expect two different effects that can be explained by the figure at right, after Fig. 4 of Fontana (2002).[90] The networks show the genotypes that map to 3 different phenotypes, P0, P1 and P2. Over time, a system may diffuse neutrally among different genotypes with the same phenotypes. Rarely a jump from one phenotype to another may occur. In the short term, evolution depends only on what is immediately accessible from a given point in genotype space. In the medium-term, evolution depends on the accessibility of alternative phenotype networks, relative to the starting network, e.g., starting with P0, P2 is twice as accessible as P1, even though P1 and P2 have the same number of genotypes. In the long-term, what matters is the total findability of a phenotype from all other phenotypes, which (as a first approximation) is a matter of the number of genotypes (and more precisely is a matter of the total surface area of the network accessible to other high-fitness phenotypes). In this case, P0 is more findable than P1 and P2 because it has twice as many genotypes. The variational bias toward more numerous phenotypes is called "phenotype bias" by Dingle, et al (2022) (see also).

This effect of findability is the formal basis for empirical and theoretical arguments in studies of the findability of regulatory network motifs or RNA fold families. In the figure at right, Dingle, et al. (2022) present evidence of a striking tendency for the most common RNA folds in nature to match the folds most widely distributed in sequence space.

Role in broad appeals for reform

The theory of arrival biases, proposed in 2001, appears in several subsequent appeals for reform, relative to the neo-Darwinian view of the Modern Synthesis. Per Arthur,[91] it is part of a developmentalist approach to evolution that emphasizes the internal organizing effects of "developmental reprogramming" on variation. In a different framing per, the efficacy of arrival biases undermines the historic commitment of theoreticians to viewing evolution as a process of shifting gene frequencies in an abundant gene pool, dominated by mass-action forces, and is part of a larger movement (beginning during the molecular revolution) away from the neo-Darwinism of the Modern Synthesis and towards a version of mutationism grounded in population genetics. The theory also has been invoked in the literature of the Extended Evolutionary Synthesis under the heading of Developmental bias.[92]

Distinction from other theories of mutational effects

The theory of arrival biases focuses on a kind of population-genetic causation linking intrinsic generative biases acting prior to selection with predictable evolutionary tendencies. It is distinct from other ideas that lack the same focus on causation, on intrinsic biases, or on the introduction process.

Evolution by mutation pressure

In classic sources, evolution by "mutation pressure" means the mass transformation of a population by mutational conversion, as in Wilson and Bossert (1971, p. 42[93] [94]). The general assessment of this theory, following Haldane (1932) [95]   and Fisher (1930),[75] is that evolution by mutation pressure is implausible because it requires high mutation rates unopposed by selection. Kimura argued even more pessimistically that transformation by mutation pressure would take so long that it can be ignored for practical purposes.[96] Nevertheless, later empirical and theoretical work showed that the theory can be valuable in cases such as the loss of a complex trait encoded by many loci, e.g., loss of sporulation in experimental populations of B. subtilis, a case in which the mutation rate for loss of the trait was estimated as an unusually high value,

\mu=0.003

.[97]

Thus, the theory of mutation pressure and the theory of arrival biases both depict ways for the process of mutation to be an important influence, but they focus on different modes of causation: influencing either the fixation process (mutation pressure) or the introduction process (arrival bias). The effectiveness of mutational tendencies via these two modes is completely different, e.g., only the mutation pressure theory relies on high mutation rates unopposed by selection.

Evolution along genetic lines of least resistance

Evolutionary quantitative genetics, the body of theory that focuses on highly polygenic quantitative traits, makes a particular prediction about mutational effects that has some empirical support. In the standard theory for a set of quantitative traits, the standing variation is represented by a

G

matrix of variances and covariances, which depends (in a complex way) on mutational input represented by an

M

matrix. Phenotypic divergence will tend to be aligned (in phenotype space) with the dimension of greatest variation,

gmax

, and this predicted effect of standing variation has been seen repeatedly.[98] [99]    This effect (explained more fully in Developmental bias) is called adaptation "along genetic lines of least resistance" and could be re-stated (with variation in a positive role) as adaptation along lines of maximal variational fuel. When divergence also aligns with

M

, this suggests that mutational variability shapes divergence, but this circumstantial correlation has other interpretations and is not taken as dispositive evidence.[100]

The use of mutation bias in the sense of an asymmetric effect on trait means is not part of the standard framework. When mutation bias is included in models of a single quantitative trait under stabilizing selection, the result is a small displacement from the optimal value.[101] [102]  [103]  

Thus, models in evolutionary quantitative genetics are focusing on a different kind of problem, so that there is no simple translation between (for instance) effects of

M

and effects of biases in introduction.

Mutational contingency

Evolutionary explanations have often relied on a paradigm of "equilibrium explanation" (ch. 5 of [104]) in which outcomes are explained by appeal to what is selectively optimal, without regard to history or details of process(as explained in [105])However, attention has focused in recent decades on the idea of "contingency", i.e., the idea that the outcome of evolution cannot be explained as the predictable or predefined endpoint of a deterministic process, but takes some path that cannot be predicted easily, or can only be predicted by knowing details of the starting conditions and the subsequent dynamics. "Mutational" contingency refers to cases in which an event of evolution is associated distinctively with a particular mutation, or a mutational hotspot, e.g.,.[106] in the sense that the evolutionary change would not have happened in the observed manner if the distinctive mutation had not occurred in the manner inferred.

This notion differs from the theory of biases in the introduction process because it is an explanatory concept (rather than a mechanism), applied in idiographic explanations, i.e., explaining one-off events (token events),[107] The theory of biases in the introduction process is a theory of general causation: the result of successfully applying the theory is to assign, not a token explanation, but a general explanation like this pattern in which

AB

happens more often than

AC

is caused by a bias in introduction due to the higher chance of the mutational-developmental conversion

AB

.

Developmental constraints (developmental bias)

The concept of "constraint" is fraught.[108] Green and Jones (2016)[109] argue that evolutionary biologists use it as a flexible explanatory concept rather than as a way to refer to a specific causal theory, i.e., a constraint is a factor with some limiting influence that makes it predictive, even if the causal basis of this influence is unclear.

A simple notion of developmental constraint is that some phenotypic forms are not observed, due to being impossible (or at least very difficult) to generate developmentally [110]  , e.g., centipedes with an even number of leg-bearing segments.[111] That is, constraint is an explanation for the non-existence of phenotypes based on a variational effect (absence), within a paradigm focused on accounting for patterns of phenotype existence.

Other references to "constraint" imply graduated differences rather than the absolute difference between possible and impossible forms, e.g.,.[112] Whereas the effectiveness of absolute biases does not require a special causal theory (because a developmentally impossible form is an evolutionarily impossible form), the idea of graduated biases prompts questions of causation, due to the conflict with the classic Haldane-Fisher "opposing pressures" argument, which holds that mere variational tendencies are ineffectual because mutation rates are small. The seminal "developmental constraints" paper by Maynard Smith, et al. (1985) [112] noted this issue without providing a solution. Advocates of "constraint" were criticized for failing to provide a mechanism.[105] This is the issue that Yampolsky and Stoltzfus [27] sought to remedy.

Nevertheless, the theory of arrival biases cannot be easily mapped to the concept of "constraint" due to the latter being used widely as a synonym for "factor". In the evo-devo literature, the term "constraint" is increasingly replaced with references to developmental bias. However, the concept of developmental bias is often associated with some idea of facilitated variation or evolvability, whereas the theory of arrival biases is only about the population-genetic consequences of arbitrary biases in the generation of variation.

Facilitated variation, evolvability, and directed mutation

The theory of arrival biases does not require or imply facilitated variation or directed mutation and is not by itself a theory of the evolution of evolvability. The population-genetic models used to illustrate the theory, and the empirical cases invoked in support of the theory, focus on the effects of different forms of mutation bias, where the bias is always relative to some dimension other than fitness, e.g., transition-transversion bias, CpG bias, or the asymmetry of two traits with different mutabilities. That is, the theory does not assume that biases are beneficial with respect to fitness, and it does not propose that mutation somehow contributes to adaptedness separately from the effect of selection (contra [113]). In fact, many models illustrate the efficacy of arrival biases by focusing on a case where the most mutationally favored outcomes are not the most fit options, as in the original Yampolsky-Stoltzfus model, where one choice has a higher mutation rate but a smaller fitness benefit, and the other has a higher fitness benefit but a smaller mutation rate. The theory assumes neither that mutationally favored outcomes are more fit, nor that they are less fit.

See also

Notes and References

  1. Book: R. Novick . 2023 . Structure and Function . Cambridge University Press, Cambridge . 9781009013888 .
  2. I. Popov . 2009 . The problem of constraints on variation, from Darwin to the present . Ludus Vitalis . 17 . 201–220.
  3. M. A. Ulett . 2014 . Making the case for orthogenesis: The popularization of definitely directed evolution (1890–1926) . Studies in History and Philosophy of Science Part C: Studies in History and Philosophy of Biological and Biomedical Sciences . 45 . 124–132 . 10.1016/j.shpsc.2013.11.009. 24368232 .
  4. A. Stoltzfus . K. Cable . 2014 . Mendelian-Mutationism: The Forgotten Evolutionary Synthesis . J Hist Biol . 47 . 4 . 501–46 . 10.1007/s10739-014-9383-2. 24811736 . 254552058 . free .
  5. N. I. Vavilov . 1922 . The Law of Homologous Series in Variation . Journal of Heredity . 12 . 47–89 .
  6. H. Spurway . 1949 . Remarks on Vavilov's law of homologous variation . La Ricerca Scientifica . 18 . 18–24 . A given species, family, order or class can mutate more readily towards certain phenotypes than others. . . From the evolutionary point of view, the mutation spectrum of a group may be more important than many of its morphological or physiological features. They [sic] determine its possiblities of evolution. Most of these are not exploited. .
  7. P. Alberch . 1980 . Ontogenesis and Morphological Diversification . American Zoologist . 20 . 4 . 653–667 . 10.1093/icb/20.4.653 . free .
  8. G. Oster . P. Alberch . 1982 . Evolution and Bifurcation of Developmental Programs . Evolution . 36 . 3 . 444–459 . 10.1111/j.1558-5646.1982.tb05066.x . 28568040 . 773428 . in evolution, selection may decide the winner of a given game but development non-randomly defines the players (p. 665) .
  9. K. S. Thomson . 1985 . Essay review: the relationship between development and evolution . Oxford Surveys in Evolutionary Biology . 2 . 220–233.
  10. Book: R. A. Raff . T. C. Kaufman . 1983 . Embryos, Genes, and Evolution . MacMillan, New York .
  11. Book: J. T. Bonner . 1982 . Evolution and Development: Report of the Dahlem Workshop on Evolution and Development Berlin 1981, May 10–15 . Springer-Verlag, New York .
  12. Book: B. C. Goodwin, N. Holder and C. C. Wylie . 1983 . Development and evolution . Cambridge University Press, Cambridge .
  13. E. S. Vrba . N. Eldredge . 1984 . Individuals, hierarchies and processes: towards a more complete evolutionary theory . Paleobiology . 10 . 2 . 146–171 . 10.1017/S0094837300008149 . 1984Pbio...10..146V . 86346683 . "Developmental biologists variously stress: (1) how indirect any genetic control is during certain stages of epigenesis; (2) that the system determines by downward causation which genomic constituents are stored in unexpressed form versus those which are expressed in the phenotype; (3) that bias in the introduction of phenotypic variation may be more important to directional phenotypic evolution than sorting by selection.".
  14. Book: J. Maynard Smith . Development and Evolution. B. C. Goodwin . N. Holder . C. C. Wylie . 1983 . Evolution and Development: Report of the Dahlem Workshop on Evolution and Development Berlin 1981, May 10–15 . Springer-Verlag, New York . This brings me to my final point. If we are to understand evolution, we must remember that it is a process which occurs in populations, not in individuals. Individual animals may dig, swim, climb or gallop, and they also develop, but they do not evolve. To attempt an explanation of evolution in terms of the development of individuals is to commit precisely that error of misplaced reductionism of which geneticists are sometimes accused.[p. 45].
  15. Book: B. Wallace . 1986 . Can Embryologists Contribute to an Understanding of Evolutionary Mechanisms? . W. Bechtel . Integrating Scientific Disciplines . Martinus Nijhoff Publishers, Dordrecht . 149–163 . "problems concerned with the orderly development of the individual are unrelated to those of the evolution of organisms through time".
  16. Book: R. Amundson . 2005 . The Changing Role of the Embryo in Evolution . Cambridge University Press, Cambridge, UK . The evolutionary irrelevance of developmental biology was argued on the basis of philosophical, methododological, and sometimes even historic grounds during the mid-twentieth to the late-twentieth century ... I do not know how to refute the irrelevancy arguments of the neo-Darwinians. It is not yet clear how this dilemma will be resolved..
  17. R. Scholl and M. Pigliucci . 2015 . The proximate–ultimate distinction and evolutionary developmental biology: causal irrelevance versus explanatory abstraction . Biology & Philosophy . 30 . 5 . 653–670 . 10.1007/s10539-014-9427-1. 254277341 .
  18. E. Mayr . 1994 . Response to John Beatty . Biology and Philosophy . 9 . 3 . 357–358 . 10.1007/BF00857941 . 84453520 . what amazes me is, in how much of the current literature that [proximate-ultimate] distinction is still ignored or confused. I must have read in the last two years four or five papers and one book on development and evolution. Now development, the decoding of the genetic program, is clearly a matter of proximate causations. Evolution, equally clearly, is a matter of evolutionary causations. And yet, in all these papers and that book the two kinds of causations were hopelessly mixed up.
  19. M. Lynch . 2007 . The frailty of adaptive hypotheses for the origins of organismal complexity . Proc Natl Acad Sci U S A . 104 . Suppl 1 . 8597–604. 10.1073/pnas.0702207104 . 17494740 . 1876435 . free .
  20. A. S. Wilkins . 1998 . Evolutionary Developmental Biology: Where is it going? . BioEssays . 20 . 10 . 783–784. 10.1002/(SICI)1521-1878(199810)20:10<783::AID-BIES1>3.0.CO;2-R . 85320436 .
  21. B. Calcott . 2009 . Lineage Explanations: Explaining How Biological Mechanisms Change . The British Journal for the Philosophy of Science . 60 . 51–78. 10.1093/bjps/axn047 . 1885/50236 . free .
  22. Book: L. Chiu . 2022 . Extended Evolutionary Synthesis. A Review of the Latest Scientific Research . John Templeton Foundation, West Conshohocken, Pennsylvania, USA . 10.15868/socialsector.40950 . 252252438 .
  23. D. R. Rokyta . P. Joyce . S. B. Caudle . H. A. Wichman . 2005 . An empirical test of the mutational landscape model of adaptation using a single-stranded DNA virus . Nat Genet . 37 . 4 . 441–4. 10.1038/ng1535 . 15778707 . 20296781 .
  24. M. T. Wortel . D. Agashe . S. F. Bailey . C. Bank . K. Bisschop . T. Blankers . J. Cairns . E. S. Colizzi . D. Cusseddu . M. M. Desai . B. van Dijk . M. Egas . J. Ellers . A. T. Groot . D. G. Heckel . M. L. Johnson . K. Kraaijeveld . J. Krug . L. Laan, M. Lassig . P. A. Lind . J. Meijer . L. M. Noble . S. Okasha . P. B. Rainey . D. E. Rozen . S. Shitut . S. J. Tans . O. Tenaillon . H. Teotonio . J. de Visser . M. E. Visser . R. M. A. Vroomans . G. D. A. Werner . B. Wertheim . P. S. Pennings . 3 . 2023 . Towards evolutionary predictions: Current promises and challenges . Evolutionary Applications . 16 . 1 . 3–21 . 10.1111/eva.13513 . 36699126 . 9850016 . In a wide range of taxa, mutation bias explains a non-negligible proportion of cases of parallel genetic evolution (Stern & Orgogozo 2008, Bailey et al. 2017, 2018, Stoltzfus & McCandlish 2017). For instance, an elegant study on adaptation to high altitudes in birds found parallel evolution, in part due to mutation bias at CpG sites (Storz et al. 2019). When highly beneficial mutations are under-sampled due to the existing mutational bias, other smaller-effect but more frequent mutations may fix instead. Such a pattern was observed in replicated evolving populations of bacteriophage (Sackman et al. 2017) where the mutation with the largest fitness effect was not the one that reached fixation most often, because its mutation rate was lower than that of other mutations with smaller fitness effects..
  25. J. B. Losos . S. J. Arnold . G. Bejerano . E. D. Brodie, 3rd . D. Hibbett . H. E. Hoekstra . D. P. Mindell . A. Monteiro . C. Moritz . H. A. Orr . D. A. Petrov . S. S. Renner . R. E. Ricklefs . P. S. Soltis . T. L. Turner . 3 . 2013 . Evolutionary biology for the 21st century . PLOS Biol . 11 . 1 . e1001466 . 10.1371/journal.pbio.1001466. 23319892 . 3539946 . free .
  26. Book: Jonathan B. Losos . David A. Baum . Douglas J. Futuyma . Hopi E. Hoekstra . Richard E. Lenski . Allen J. Moore . Catherine L. Peichel . Dolph Schluter . Michael J. Whitlock . 2014 . The Princeton Guide to Evolution . Princeton University Press .
  27. L. Y. Yampolsky and A. Stoltzfus . 2001 . Bias in the introduction of variation as an orienting factor in evolution . Evol Dev . 3 . 2 . 73–83 . 10.1046/j.1525-142x.2001.003002073.x. 11341676 . 26956345 .
  28. A. A. Soares . L. Wardil . L. B. Klaczko . R. Dickman . 2021 . Hidden role of mutations in the evolutionary process . Phys Rev E . 104 . 4 . 044413 . 10.1103/PhysRevE.104.044413. 34781575 . 2021PhRvE.104d4413S . 239954422 .
  29. D. M. McCandlish . A. Stoltzfus . 2014 . Modeling evolution using the probability of fixation: history and implications . Quarterly Review of Biology . 89 . 3 . 225–252 . 10.1086/677571 . 25195318 . 19619966 .
  30. A. V. Cano . H. Rozhonova . A. Stoltzfus . D. M. McCandlish . J. L. Payne . 2022-02-10 . Mutation bias shapes the spectrum of adaptive substitutions . Proc Natl Acad Sci U S A . 119 . 7 . 10.1073/pnas.2119720119. 35145034 . 8851560 . 2022PNAS..11919720C .
  31. S. F. Bailey . F. Blanquart . T. Bataillon . R. Kassen . 2017 . What drives parallel evolution?: How population size and mutational variation contribute to repeated evolution . BioEssays . 39 . 1 . 1–9 . 10.1002/bies.201600176. 27859467 . 193595 . free .
  32. A. V. Cano . B. Gitschlag . H. Rozhonova . A. Stoltzfus . D. M. McCandlish . J. L. Payne . 2023 . Proc. R. Soc. B . 378 . 20220055 . Mutation bias and the predictability of evolution . 1877 . 10.1098/rstb.2022.0055. 37004719 . 10067271 .
  33. K. Gomez . J. Bertram . J. Masel . 2020 . Mutation bias can shape adaptation in large asexual populations experiencing clonal interference . Proc. R. Soc. B . 287 . 1937 . 20201503 . 10.1098/rspb.2020.1503. 33081612 . 7661309 .
  34. A. F. Shull . 1935 . Weismann and Haeckel: One Hundred Years . Science . 81 . 2106 . 443–451. 10.1126/science.81.2106.443 . 17818758 . 1935Sci....81..443S .
  35. L. M. Chevin, G. Martin and T. Lenormand . 2010 . Fisher's model and the genomics of adaptation: restricted pleiotropy, heterogenous mutation, and parallel evolution . Evolution . 64 . 11 . 3213–31 . 10.1111/j.1558-5646.2010.01058.x. 20662921 . 205782829 . free .
  36. A. Stoltzfus . 2006 . Mutation-Biased Adaptation in a Protein NK Model . Mol Biol Evol . 23 . 10 . 1852–1862 . 10.1093/molbev/msl064 . 16857856 .
  37. Cano . Alejandro . Payne . Joshua L. . Mutation bias interacts with composition bias to influence adaptive evolution . PLOS Computational Biology . 2020 . 16 . 9 . e1008296 . 10.1371/journal.pcbi.1008296. 32986712 . 7571706 . 2020PLSCB..16E8296C . free .
  38. Schaper . S. . Louis . A.A. . The arrival of the frequent: how bias in genotype-phenotype maps can steer populations to local optima . PLOS ONE . 2014 . 9 . 2 . e86635 . 10.1371/journal.pone.0086635 . 24505262. 3914804 . 2014PLoSO...986635S . free .
  39. K. Xiong . M. Gerstein . J. Masel . 2021 . Differences in evolutionary accessibility determine which equally effective regulatory motif evolves to generate pulses . Genetics . 219 . 3 . 10.1093/genetics/iyab140. 34740240 . 8570775 .
  40. J. Garson, L. Wang and S. Sarkar . 2003 . How Development May Direct Evolution . Biology and Philosophy . 18 . 2 . 353–370 . 10.1023/A:1023996321257 . 83110476 .
  41. Book: E. Mayr . 1980 . Some Thoughts on the History of the Evolutionary Synthesis . E. Mayr and W. Provine . The Evolutionary Synthesis: Perspectives on the Unification of Biology . Harvard University Press . 1–48 .
  42. Book: C. Darwin . 1868 . Variation of Animals and Plants under Domestication . Murray, London . Throughout this chapter and elsewhere I have spoken of selection as the paramount power, yet its action absolutely depends on what we in our ignorance call spontaneous or accidental variability. Let an architect be compelled to build an edifice with uncut stones, fallen from a precipice. The shape of each fragment may be called accidental; yet the shape of each has been determined by the force of gravity, the nature of the rock, and the slope of the precipice, – events and circumstances, all of which depend on natural laws; but there is no relation between these laws and the purpose for which each fragment is used by the builder. In the same manner the variations of each creature are determined by fixed and immutable laws; but these bear no relation to the living structure which is slowly built up through the power of selection, whether this be natural or artificial selection [emphasis added by editor].
  43. R. Lande . 1979 . Quantitative Genetic Analysi of Multivariate Evolution, Applied to Brain:Body Size Allometry . Evolution . 33 . 1Part2 . 402–416. 10.1111/j.1558-5646.1979.tb04694.x . 28568194 . 17132606 .
  44. R. Lande and S. J. Arnold . 1983 . The Measurement of Selection on Correlated Characters . Evolution . 37 . 6 . 1210–1226. 10.2307/2408842 . 2408842 . 28556011 .
  45. J. Brookfield . 2005 . Where does animal diversity come from? . Current Biology . 15 . 22 . R908–R910. 10.1016/j.cub.2005.11.007 . 10714428 . free .
  46. B. L. Gitschlag . A. V. Cano . J. L. Payne . D. M. McCandlish . A. Stoltzfus . 2023 . Mutation and selection induce correlations between selection coefficients and mutation rates . The American Naturalist . 202 . 4 . 534–557 . 10.1086/726014. 37792926 . 20.500.11850/639862 . 256902349 . free .
  47. M. A. Streisfeld and M. D. Rausher . 2011 . Population genetics, pleiotropy, and the preferential fixation of mutations during adaptive evolution . Evolution . 65 . 3 . 629–42 . 10.1111/j.1558-5646.2010.01165.x. 21054357 . 33188741 . free .
  48. A. Farlow, E. Meduri and C. Schlotterer . 2011 . DNA double-strand break repair and the evolution of intron density . Trends in Genetics . 27 . 1 . 1–6 . 10.1016/j.tig.2010.10.004. 21106271 . 3020277 .
  49. A. Stoltzfus and L. Y. Yampolsky . 2009 . Climbing Mount Probable: mutation as a cause of nonrandomness in evolution . J Hered . 100 . 5 . 637–47 . 10.1093/jhered/esp048 . 19625453 . free .
  50. M. Tuffaha . S. Varakunan . D. Castellano . R. N. Gutenkunst . L. M. Wahl . 2022 . Shifts in mutation bias promote mutators by altering the distribution of fitness effects . bioRxiv . 2022.09.27.509708 . 10.1101/2022.09.27.509708. 252684737 .
  51. M. Sane . G. D. Diwan . B. A. Bhat . L. M. Wahl . D. Agashe . 2023 . Shifts in mutation spectra enhance access to beneficial mutations . Proc Natl Acad Sci U S A . 120 . 22 . e2207355120 . 10.1073/pnas.2207355120. 37216547 . 10235995 . 2023PNAS..12007355S . free .
  52. A. Couce, A. Rodríguez-Rojas and J. Blázquez . 2015 . Bypass of genetic constraints during mutator evolution to antibiotic resistance . Proceedings of the Royal Society of London B: Biological Sciences . 282.
  53. J. S. Horton . L. M. Flanagan . R. W. Jackson . N. K. Priest . T. B. Taylor . 2021 . A mutational hotspot that determines highly repeatable evolution can be built and broken by silent genetic changes . Nat Commun . 12 . 1 . 6092 . 10.1038/s41467-021-26286-9. 34667151 . 8526746 . 2021NatCo..12.6092H .
  54. A. Stoltzfus . D. M. McCandlish . 2017 . Mutational biases influence parallel adaptation . Mol Biol Evol . 34 . 9 . 2163–2172 . 10.1093/molbev/msx180. 28645195 . 5850294 .
  55. A. M. Sackman . L. W. McGee . A. J. Morrison . J. Pierce . J. Anisman . H. Hamilton . S. Sanderbeck . C. Newman . D. R. Rokyta . 2017 . Mutation-Driven Parallel Evolution During Viral Adaptation . Mol Biol Evol . 34 . 12 . 3243–3253 . 10.1093/molbev/msx257. 29029274 . 5850295 .
  56. J. Y. Pham . C. B. Ogbunugafor . A. N. Nguyen Ba . D. L. Hartl . 2019 . Experimental evolution for niche breadth in bacteriophage T4 highlights the importance of structural genes . MicrobiologyOpen . 9 . 2 . e968 . 10.1002/mbo3.968. 31778298 . 7002106 .
  57. J. L. Payne . F. Menardo . A. Trauner . S. Borrell . S. M. Gygli . C. Loiseau . S. Gagneux . A. R. Hall . 2019 . Transition bias influences the evolution of antibiotic resistance in Mycobacterium tuberculosis . PLOS Biol . 17 . 5 . e3000265 . 10.1371/journal.pbio.3000265. 31083647 . 6532934 . free .
  58. A. Stoltzfus . D. M. McCandlish . 2015 . Mutation-biased adaptation in Andean house wrens . Proc Natl Acad Sci U S A . 112 . 45 . 13753–4 . 10.1073/pnas.1518490112. 26489650 . 4653193 . 2015PNAS..11213753S . free .
  59. R. C. MacLean . G. G. Perron . A. Gardner . 2010 . Diminishing returns from beneficial mutations and pervasive epistasis shape the fitness landscape for rifampicin resistance in Pseudomonas aeruginosa . Genetics . 186 . 4 . 1345–54 . 10.1534/genetics.110.123083. 20876562 . 2998316 .
  60. D. Coombes . J. W. B. Moir . A. M. Poole . T. F. Cooper . R. C. J. Dobson . 2019 . The fitness challenge of studying molecular adaptation . Biochem Soc Trans . 47 . 5 . 1533–1542 . 10.1042/BST20180626. 31642877 . 204849713 .
  61. J. F. Storz . C. Natarajan . A. V. Signore . C. C. Witt . D. M. McCandlish . A. Stoltzfus . 2019 . The role of mutation bias in adaptive molecular evolution: insights from convergent changes in protein function . Philos Trans R Soc Lond B Biol Sci . 374 . 1777 . 20180238 . 10.1098/rstb.2018.0238. 31154983 . 6560279 .
  62. R. Hershberg and D. A. Petrov . 2010 . Evidence that mutation is universally biased towards AT in bacteria . PLOS Genet . 6 . 9 . e1001115 . 10.1371/journal.pgen.1001115. 20838599 . 2936535 . free .
  63. B. Ujvari . N. R. Casewell . K. Sunagar . K. Arbuckle . W. Wuster . N. Lo . D. O'Meally . C. Beckmann . G. F. King . E. Deplazes . T. Madsen . 2015 . Widespread convergence in toxin resistance by predictable molecular evolution . Proc Natl Acad Sci U S A . 112 . 38 . 11911–6 . 10.1073/pnas.1511706112. 26372961 . 4586833 . 2015PNAS..11211911U . free .
  64. C. R. Feldman . E. D. Brodie Jr. . E. D. Brodie, 3rd . M. E. Pfrender . 2012 . Constraint shapes convergence in tetrodotoxin-resistant sodium channels of snakes . Proceedings of the National Academy of Sciences of the United States of America . 109 . 12 . 4556–61 . 10.1073/pnas.1113468109. 22392995 . 3311348 . free .
  65. M. C. Jost . D. M. Hillis . Y. Lu . J. W. Kyle . H. A. Fozzard . H. H. Zakon . 2008 . Toxin-resistant sodium channels: parallel adaptive evolution across a complete gene family . Molecular Biology and Evolution . 25 . 6 . 1016–24 . 10.1093/molbev/msn025. 18258611 . 2877999 .
  66. L. Smeds, A. Qvarnstrom and H. Ellegren . 2016 . Direct estimate of the rate of germline mutation in a bird . Genome Res . 26 . 9 . 1211–8 . 10.1101/gr.204669.116. 27412854 . 5052036 .
  67. A. Stoltzfus . 2006 . Mutationism and the Dual Causation of Evolutionary Change . Evol Dev . 8 . 3 . 304–317 . 10.1111/j.1525-142X.2006.00101.x . 16686641 . 10469049 .
  68. Book: A. Stoltzfus . Understanding bias in the introduction of variation as an evolutionary cause . Uller . T. . Laland . K.N. . Evolutionary Causation: Biological and Philosophical Reflections . 2019 . MIT Press . Cambridge, MA.
  69. Book: A. Stoltzfus . 2021 . Mutation, Randomness and Evolution . Oxford, Oxford .
  70. N. Sueoka . 1988 . Directional mutation pressure and neutral molecular evolution . Proc Natl Acad Sci U S A . 85 . 8 . 2653–7. 10.1073/pnas.85.8.2653 . 3357886 . 280056 . 1988PNAS...85.2653S . free .
  71. K. Dingle . F. Ghaddar . P. Sulc . A. A. Louis . 2022 . Phenotype Bias Determines How Natural RNA Structures Occupy the Morphospace of All Possible Shapes . Mol Biol Evol . 39 . 10.1093/molbev/msab280. 34542628 . 8763027 .
  72. Book: J. L. King . 1971 . The Influence of the Genetic Code on Protein Evolution . E. Schoffeniels . Biochemical Evolution and the Origin of Life . North-Holland Publishing Company . 3–13 . The frequency of each amino acid is proportional to the probability of its emergence through mutation..
  73. Book: S. A. Kauffman . 1993 . The origins of order: Self-organization and evolution . Oxford University Press, New York .
  74. J. B. S. Haldane . 1927 . A mathematical theory of natural and artificial selection. V. Selection and mutation. . Proc. Camb. Phil. Soc. . 26 . 220–230 . 10.1017/S0305004100015644. 86716613 .
  75. Book: R. A. Fisher . 1930 . The Genetical Theory of Natural Selection . Oxford University Press, London .
  76. Book: T. Dobzhansky . 1937 . Genetics and the Origin of Species . Columbia University Press, New York.
  77. Book: W. B. Provine . 1971 . The Origins of Theoretical Population Genetics . University of Chicago Press, Chicago .
  78. Book: A. W. F. Edwards . 1977 . Foundations of Mathematical Genetics . Cambridge University Press, New York . All genes will be assumed stable, and mutation will not be taken into account [p. 3].
  79. Book: R. C. Lewontin . 1974 . The Genetic Basis of Evolutionary Change . Columbia University Press, New York . Deterministic theory depends critically on the parameters of selection. There is virtually no qualitative or gross quantitative conclusion about the genetic structure of populations in deterministic theory that is sensitive to small values of migration, or any that depends on mutation rates.[p. 267].
  80. Book: S. J. Gould . 2002 . The Structure of Evolutionary Theory . Harvard University Press, Cambridge, Massachusetts .
  81. J. Beatty . 2022 . The synthesis and the two scenarios . Evolution . 76 . S1 . 6–14 . 10.1111/evo.14423. 34951705 . 245445537 .
  82. D. J. Futuyma . 2001 . Evolution, Science, and Society . The American Naturalist . 158 . S1–S46 . 10.1086/509090. 18707297 . 41606175 .
  83. G. Yedid and G. Bell . 2002 . Macroevolution simulated with autonomously replicating computer programs . Nature . 420 . 6917 . 810–2 . 10.1038/nature01151 . 12490948 . 2002Natur.420..810Y . 607532 . "The process of adaptation occurs on two timescales. In the short term, natural selection merely sorts the variation already present in a population, whereas in the longer term genotypes quite different from any that were initially present evolve through the cumulation of new mutations. The first process is described by the mathematical theory of population genetics. However, this theory begins by defining a fixed set of genotypes and cannot provide a satisfactory analysis of the second process because it does not permit any genuinely new type to arise." .
  84. D. L. Hartl and C. H. Taubes . 1998 . Towards a theory of evolutionary adaptation . Genetica . 103 . 525–533 . 10.1023/A:1017071901530 . 2057562 . Almost every theoretical model in population genetics can be classified into one of two major types. In one type of model, mutations with stipulated selective effects are assumed to be present in the population as an initial condition ... The second major type of models does allow mutations to occur at random intervals of time, but the mutations are assumed to be selectively neutral or nearly neutral..
  85. Book: I. Eshel and M. W. Feldman . 2001 . 161–190 . Optimality and Evolutionary Stability under Short-term and Long-term Selection . S.H. Orzack and E. Sober . Adaptationism and optimality . Cambridge University Press, Cambridge . "Since the time of Fisher, an implicit working assumption in the quantitative study of evolutionary dynamics is that qualitative laws governing long-term evolution can be extrapolated from results obtained for the short-term process. We maintain that this extrapolation is not accurate. The two processes are qualitatively different from each other. One obvious difference is that the mutational changes that cause the transition between one set of genotype frequencies, close to its short-term equilibrium, and the next set are random even in large populations." .
  86. Book: S. C. Stearns and R. F. Hoekstra . 2005 . Evolution: an introduction . Oxford University Press, New York . the rate of evolutionary change is not limited by the occurrence of new favorable mutations (p. 99) .
  87. R. Dawkins . 2007 . Review: The Edge of Evolution . International Herald Tribune . 29 June . If correct, Behe's calculations would at a stroke confound generations of mathematical geneticists, who have repeatedly shown that evolutionary rates are not limited by mutation. Single-handedly, Behe is taking on Ronald Fisher, Sewall Wright, J.B.S. Haldane, Theodosius Dobzhansky, Richard Lewontin, John Maynard Smith and hundreds of their talented co-workers and intellectual descendants. Notwithstanding the inconvenient existence of dogs, cabbages and pouter pigeons, the entire corpus of mathematical genetics, from 1930 to today, is flat wrong. Michael Behe, the disowned biochemist of Lehigh University, is the only one who has done his sums right. You think? The best way to find out is for Behe to submit a mathematical paper to The Journal of Theoretical Biology, say, or The American Naturalist, whose editors would send it to qualified referees. .
  88. J. Hermisson and P. S. Pennings . 2005 . Soft sweeps: molecular population genetics of adaptation from standing genetic variation . Genetics . 169 . 4 . 2335–52 . 10.1534/genetics.104.036947 . 15716498 . 1449620 .
  89. D. Houle, B. Morikawa and M. Lynch . 1996 . Comparing mutational variabilities . Genetics . 143 . 3 . 1467–83. 10.1093/genetics/143.3.1467 . 8807316 . 1207413 .
  90. W. Fontana . 2002 . Modelling 'evo-devo' with RNA . BioEssays . 24 . 12 . 1164–77 . 10.1002/bies.10190. 12447981 . free .
  91. Book: W. Arthur . 2004 . Biased Embryos and Evolution . Cambridge University Press, Cambridge .
  92. T. Uller . A. P. Moczek . R. A. Watson . P. M. Brakefield . K. N. Laland . 2018 . Developmental Bias and Evolution: A Regulatory Network Perspective . Genetics . 209 . 4 . 949–966 . 10.1534/genetics.118.300995. 30049818 . 6063245 .
  93. Book: E. O. Wilson and W. H. Bosssert . 1971 . A Primer of Population Biology . Sinauer Sunderland, MA .
  94. An online demonstration is included in the educational tools from population geneticist John McDonald at https://udel.edu/~mcdonald/evolmutpress.html)
  95. Book: J. B. S. Haldane . 1932 . The Causes of Evolution . Longmans, Green and Co., New York .
  96. M. Kimura . 1980 . Average time until fixation of a mutant allele in a finite population under continued mutation pressure: Studies by analytical, numerical, and pseudo-sampling methods . Proc Natl Acad Sci U S A . 77 . 1 . 522–526 . 10.1073/pnas.77.1.522 . 16592764 . 348304 . 1980PNAS...77..522K . free .
  97. H. Maughan . J. Masel . C. W. Birky Jr. . W. L. Nicholson . 2007 . The roles of mutation accumulation and selection in loss of sporulation in experimental populations of Bacillus subtilis . Genetics . 177 . 2 . 937–48 . 10.1534/genetics.107.075663. 17720926 . 2034656 .
  98. D. Schluter . 1996 . Adaptive Radiation Along Genetic Lines of Least Resistance . Evolution . 50 . 5 . 1766–1774 . 10.2307/2410734 . 2410734 . 28565589 .
  99. K. McGuigan . 2006 . Studying phenotypic evolution using multivariate quantitative genetics . Mol Ecol . 15 . 4 . 883–96 . 10.1111/j.1365-294X.2006.02809.x. 16599954 . 32969263 . free .
  100. D. Houle . G. H. Bolstad . K. van der Linde . T. F. Hansen . 2017 . Mutation predicts 40 million years of fly wing evolution . Nature . 548 . 7668 . 447–450 . 10.1038/nature23473. 28792935 . 2017Natur.548..447H . 90225702 . 11250/2453538 . free .
  101. D. Waxman and J. R. Peck . 2003 . The anomalous effects of biased mutation . Genetics . 164 . 4 . 1615–26 . 10.1093/genetics/164.4.1615 . 12930765 . 1462658 .
  102. X. S. Zhang and W. G. Hill . 2008 . The anomalous effects of biased mutation revisited: mean-optimum deviation and apparent directional selection under stabilizing selection . Genetics . 179 . 2 . 1135–41 . 10.1534/genetics.107.083428. 18558659 . 2429868 .
  103. B. Charlesworth . 2013 . Stabilizing selection, purifying selection, and mutational bias in finite populations . Genetics . 194 . 4 . 955–71 . 10.1534/genetics.113.151555. 23709636 . 3730922 .
  104. Book: E. Sober . 1984 . The Nature of Selection: Evolutionary Theory in Philosophical Focus . MIT Press, Cambridge, Mass. .
  105. H. K. Reeve and P. W. Sherman . 1993 . Adaptation and the Goals of Evolutionary Research . Quarterly Review of Biology . 68 . 1–32 . 10.1086/417909 . 84402115 . Some analogies may be helpful here. In explaining why certain traits predominate in nature rather than conceivable others, evolutionary biologists act a little like physicists who study stable states in thermodynamic systems or astronomers who study spatial arrangements of stars and galaxies' In each case, the researcher examines the products of historical processes that unfold over centuries due to the interactions of particular entities. The products of these processes are explicable, sometimes even predictable, without a detailed knowledge of each entity's individual history. This is so because, for a given set of conditions, a large number of distinct individual histories tend to converge on a small number of stable states. For example, the observation that boulders collect in mountain valleys can be explained without knowing each boulder's precise trajectory..
  106. T. Lenormand, D. Roze and F. Rousset . 2009 . Stochasticity in evolution . Trends Ecol Evol . 24 . 3 . 157–65. 10.1016/j.tree.2008.09.014 . 19178980 .
  107. Book: J. L. King . 2017 . Wonderful Life Revisited: Chance and Contingency in the Ediacaran-Cambrian Radiation . G. Ramsey and C. H. Pence . Chance in Evolution . University of Chicago Press .
  108. J. Antonovics and P. H. van Tienderen . 1991 . Ontoecoenophyloconstraints? The Chaos of Constraint Terminoloy . Trends Ecol. Evol. . 6 . 5 . 166–168. 10.1016/0169-5347(91)90059-7 . 21232448 .
  109. S. Green and N. Jones . 2016 . Constraint-Based Reasoning for Search and Explanation: Strategies for Understanding Variation and Patterns in Biology . Dialectica . 70 . 3 . 343–374 . 10.1111/1746-8361.12145.
  110. M. E. Olson . 2012 . The developmental renaissance in adaptationism . Trends Ecol Evol . 27 . 5 . 278–87 . 10.1016/j.tree.2011.12.005. 22326724 .
  111. W. Arthur and M. Farrow . 1999 . The pattern of variation in centipede segment number as an example of developmental constraint in evolution . J Theor Biol . 200 . 2 . 183–91. 10.1006/jtbi.1999.0986 . 10504284 . 1999JThBi.200..183A .
  112. J. Maynard Smith . R. Burian . S. Kauffman . P. Alberch . J. Campbell . B. Goodwin . R. Lande . D. Raup . L. Wolpert . 1985 . Developmental Constraints and Evolution . Quart. Rev. Biol. . 60 . 265–287 . 10.1086/414425 . 85201850 . Since the classic work of Fisher (1930) and Haldane (1932) established the weakness of directional mutation as compared to selection, it has been generally held that directional bias in variation will not produce evolutionary change in the face of opposing selection. This position deserves reexamination. .
  113. E. I. Svensson and D. Berger . 2019 . The Role of Mutation Bias in Adaptive Evolution . Trends Ecol Evol . 34 . 5 . 422–434 . 10.1016/j.tree.2019.01.015. 31003616 . 125066709 .