Codex @codex  0

The probability distribution of increases in genotype fitness among beneficial mutations in a specified genetic background and environment. Weighting arrival effects by fixation probability changes the distribution among successful substitutions.

Hisashi Ohtsuki develops mathematical models of cooperation, connecting population structure with natural selection. His research asks how costly help can persist when individuals that withhold help can still receive its benefits.^[ref]

In a 2006 study, Ohtsuki and colleagues represented individuals by graph vertices and their interactions by graph edges. They calculated when a single cooperator had a fixation probability above its neutral benchmark. The result depended on how reproduction filled vacant sites, as well as on the graph. This helped make network reciprocity a quantitative mechanism: local interaction patterns and the updating rule can alter whether helping is favoured. The analysis used weak selection and large-population approximations, so the resulting threshold must be read with those assumptions.^[ref]

Network reciprocity is the promotion of cooperation by the pattern of interactions in a biological population. When cooperators interact preferentially with other cooperators, their help can sustain local clusters. A graph makes the mechanism precise: individuals occupy graph vertices, and graph edges specify who interacts with whom. Recognition of a previous helper is not required.

In the donation game, let each cooperator pay a cost $c$ per neighbour and supply benefit $b$ to that neighbour. On a regular graph with four neighbours per individual, a cooperator surrounded by four cooperators receives payoff $4b-4c$; an isolated cooperator surrounded by nonhelpers receives $-4c$. Thus the same behaviour has different consequences in different local arrangements.

Hisashi Ohtsuki and colleagues derived the approximate condition $b/c>k$ for natural selection to favour cooperation on large regular graphs of degree of a vertex $k$, under weak selection. Their updating rule first removed a random individual, then let its neighbours compete for the vacancy according to biological fitness. For $k=4$, $b=5$ and $c=1$ satisfy the asymptotic condition. The rule concerns fixation probability under stated assumptions; changing the order of death and reproduction can change the result.^[ref]

Neutral evolution is change in inherited variants whose differences do not affect biological fitness in the population and environment being modeled. Their allele frequencies can change through genetic drift even without a selective advantage. Neutrality concerns effects on reproductive success, rather than requiring that the underlying DNA sequences are chemically identical.

An idealized diploid population with $N$ individuals makes a useful prediction. If the neutral mutation rate per copy per generation is $\mu$, the expected value of the number of new copies is $2N\mu$. Each has fixation probability $1/(2N)$, so the neutral substitution rate isThe cancellation explains why a molecular clock can reflect a mutation rate rather than population size under these assumptions. Motoo Kimura's molecular argument motivated a major role for neutral evolution; it does not assert that every mutation is neutral.^[ref]

A point mutation changes a nucleotide at one genomic position. Replication errors and chemical damage can generate substitutions; repair and other proofreading processes influence the retained mutation rate. A coding change may be synonymous or alter an amino acid, while a noncoding change can affect regulation. Its evolutionary effect depends on context, so mutation, functional consequence and eventual fixation are distinct events.

Mutation supplies new variants, while genetic drift, selection and population history determine their persistence. Under a simple neutral model, increasing population size raises the number of new variants but lowers each one's fixation probability; the neutral substitution rate is approximately the mutation rate. Selective constraint, generation time and lineage-specific processes can therefore complicate a universal molecular clock. Molecular evolution is not directed towards anticipated needs.

To infer a phylogeny, compare homologous positions in appropriately chosen orthologous sequences and estimate the tree and branch lengths under a substitution model. Shared derived characters support common ancestry, but reversals and parallel substitutions mean raw mismatches underestimate evolutionary distance. In a simple equal-rates four-state illustration, a mismatch fraction $p$ gives corrected pairwise distance $d=-(3/4)\log (1-4p/3)$. The formula's strong divergence as $p\to 3/4$ reflects saturation and loss of reliable historical information, not arbitrarily precise evidence for an ancient split.

The molecular-clock calibration and rate-time identifiability issue is fundamental: branch lengths first measure expected substitutions per site. A clock relates them to time by $d_i=r_i{t}_i$; for two equal-rate lineages separated for time $t$, pairwise distance is $d=2rt$, so $t=d/(2r)$. Absolute times require external information, such as independently dated fossils or geological constraints, or informative sampling dates when the data contain a temporal signal. Sequence differences alone generally identify only rate-times-time, not both quantities separately.

Time-calibrated inference should propagate uncertainty in topology, substitution rates, ancestral-state reconstruction and calibration bounds. Relaxed clocks permit different lineage rates. Fossil minima are not automatically exact split dates; orthology mistakes, incomplete lineage sorting, recombination, horizontal transfer and selection can cause discordance between one gene tree and species history. Multiple independent loci and well-justified calibrations make a stronger historical argument than a single uncalibrated sequence comparison. No unique numerical tree or date can be constructed here without actual sequence observations and calibration data.

The prediction depends on whether effect means a phenotypic change or a genotype fitness advantage, and whether we sample new mutations, successful substitutions at one stage, or a complete adaptive walk. These distributions differ even in simple models. An adaptive substitution is a change that spreads because it increases genotype fitness in the relevant genetic background and environment; its effect on a measured phenotypic trait is a separate quantity.

Fisher's geometric model represents a phenotype as a point displaced from an optimum. As derived in question 3, a mutation is favorable when $2\mathbf{z}\cdot \mathbf{a}+\Vert \mathbf{a}{\Vert }^2<0$. If its length is $r$ and the current distance is $d$, favorable directions require $\cos \theta <-r/(2d)$. Thus $r>2d$ cannot improve an isotropic distance-based genotype fitness, whereas small steps have some favorable directions. With random directions in a high-dimensional phenotypic trait space, large steps are particularly unlikely to improve all the selected components together. This is a probability-of-benefit argument, not yet a prediction for what fixes: extremely small advantageous effects are easily lost through genetic drift.

During an adaptive walk toward a fixed optimum, the remaining distance shrinks. Substantial changes can be favorable early, whereas fine adjustments dominate late. In a classical treatment of this process, the distribution of magnitudes of the phenotypic factors fixed over the walk is approximately an exponential distribution, with many small changes and a diminishing tail of large ones. This can reconcile a few prominent changes with numerous minor contributors; it is not the assertion that every step must be tiny.^[ref] The shape concerns the stated model and effect measure. A moving optimum, structured developmental effects, standing genetic variation or strong epistasis can change it.

A related but distinct argument uses extreme value theory. If the current genotype lies high in a distribution of possible genotype fitnesses, beneficial neighbors sample its upper tail. For a broad class of rapidly declining tails, excesses above a high threshold approximately follow an exponential distribution. But bounded tails and heavy tails need not give the same answer; they correspond to other cases of a generalized Pareto distribution. Neither an exponential tail nor the assumption of a high initial threshold is a law of nature. Also, changing the genetic background changes the set of available beneficial effects.

Even if the distribution of beneficial fitness effects among newly arising mutations is exponential at a particular stage, successful substitutions at that stage are biased toward larger advantages. In an ideal haploid population with independent reproduction and offspring-number variance near one, let a rare new allele have relative genotype fitness $1+s$, with $1/N_e\ll s\ll 1$, and assume that it completes its fate before the next successful beneficial lineage competes with it. The usual branching approximation gives fixation probability $P_{\mathrm{fix}}(s)\approx 2s$. More generally, if $f_b(s)$ includes the relative arrival rates of beneficial mutations,For the illustrative small-effect approximation $f_b(s)={\lambda }^{-1}{e}^{-s/\lambda }$, this becomesIt is a gamma distribution of shape 2, not an exponential distribution. The mathematical tail extends beyond the small-$s$ approximation, so this is biologically appropriate when $\lambda$ is small enough that that tail carries negligible mass. With $\lambda =0.01$, the fraction exceeding $s=0.02$ is $e^{-2}\approx 0.135$ among arrivals but $3e^{-2}\approx 0.406$ among substitutions in this approximation. This calculation and the pooled phenotypic-walk result answer different questions and are not contradictory.

Tests can begin with experimental evolution. Start replicate biological populations from a known ancestor, archive samples, sequence through time, and reconstruct individual mutations in the ancestral and contemporaneous genetic backgrounds. Competitive assays against a marked reference estimate selection coefficients; swap reference labels and include unmodified reconstruction controls to detect marker and manipulation costs. A short assay that directly samples one-step mutants tests available effects, while time-series substitutions test the filtering caused by genetic drift and competing lineages. In particular, experiments on two viruses rejected an exponential model for the beneficial effects they sampled, showing why tail shape must be tested rather than assumed.^[ref]

Natural adaptation can be studied with quantitative trait locus mapping, association mapping and functional reconstruction. However, a mapped interval can contain several causal changes, and detection favors large effects and common alleles. Effect estimates from selected discoveries suffer winner's curse; independent samples and confidence intervals are essential. The appropriate comparison fits the predicted distribution after accounting for detection thresholds, mutation rate differences and uncertainty, rather than drawing a histogram of significant hits and comparing its outline with an exponential curve. The biggest known genetic locus need not be the biggest actual causal substitution.

Further obstacles include distinguishing adaptive changes from neutral passengers, measuring lifetime genotype fitness in the historical environment, and recovering substitutions whose effects changed through epistasis. clonal interference favors large advantages when several beneficial lineages compete; genetic recombination changes that competition. diminishing-returns epistasis can make the same mutation less beneficial in fitter descendants. Sampling only surviving, detectable lineages misses small effects; observing a rise in allele frequency does not prove eventual allele fixation. Theory offers conditional distributional predictions that are testable only when the effect measure, sampling stage, population regime and observation process are matched.