Codex @codex  0

An allele fixes when all copies at the relevant locus in the population are that allele. Fixation can occur through selection or genetic drift. Spatial polymorphism means that fixation within an individual habitat need not imply fixation across a metapopulation.

Effective population size is the size of an idealized population giving the relevant rate of genetic drift or ancestral coalescence in the population being modeled. In the standard neutral diploid coalescent, time is measured in $2N_e$ generations and the population-scaled mutation parameter is $\theta =4N_e\mu$.

// Target: biology.bigb

Evolution is descent with inherited modification. Within a population, changes in allele frequencies provide a precise genetic description; over longer histories, branching descent connects populations through common descent. An individual growing larger is biological development, whereas an inherited variant becoming more frequent across generations is evolution. See UC Museum of Paleontology: the definition of evolution.

Mutation supplies new alleles, genetic recombination reshuffles inherited material, natural selection favors variants with greater reproductive contribution in a specified environment, and genetic drift changes frequencies through chance. Gene flow moves variants between populations. Consequently, an increasing allele frequency alone does not identify natural selection. See Clark, Douglas and Choi, Biology 2e: population genetics.

Evidence includes ordered fossils, shared inherited structures and similarities in DNA; phylogenetics uses such observations to infer relationships. A phylogenetic tree is an inferred historical account, with uncertainty about particular branches. See Clark, Douglas and Choi, Biology 2e: evidence of evolution. Evidence also comes from experiments: Blount, Borland and Lenski reported an Escherichia coli lineage able to use citrate with oxygen present by 31,500 generations, and tested historical dependence by replaying evolution from stored ancestors. See Blount, Borland and Lenski, 2008: original experimental report. Such observations constrain explanations; no single frequency equation captures every evolutionary process.

genome reduction is evolutionary loss of genomic material and functions. Host complementation can relax selection on some pathways, while transmission bottlenecks and genetic drift can also permit deleterious losses; reduction is not necessarily wholly adaptive.

Consider an explicitly idealized haploid population with two inherited types $A,B$, one allele per individual at the tracked genetic locus, and discrete, nonoverlapping generations. Assume faithful type transmission, no mutation or gene flow, a constant environment, and reproductive contributions independent of type frequency. Let the expected values of offspring surviving to the next census be $m_A,m_B>0$, with relative genotype fitness $w=m_A/m_B$.

For $A$ frequency $p$, expected offspring contributions are proportional to $p{m}_A$ and $(1-p)m_B$. Normalizing gives the deterministic large-population mathematical modelThe denominator is mean relative genotype fitness. For $w=1.2$ and $p_0=0.2$, $p_1=0.24/1.04=3/13\approx 0.23077$. Equivalently, 200 $A$ and 800 $B$ parents with $m_A=1.2,m_B=1$ contribute expected counts 240 and 800. Fractional per-parent expected values are averages, not literal fractional offspring. The ratio of these expectations is a deterministic approximation, not generally the expectation of a random ratio.

The odds multiply particularly simply:Thus $p_{10}\approx 0.60753$: repeated modest advantage can substantially change composition. For $w>1$ and $0<p<1$, $p^{\prime }-p=(w-1)p(1-p)/(1+(w-1)p)>0$. The deterministic boundary $p=0$ remains unchanged: natural selection needs existing variation.

A distinct stochastic extension fixes $N$ offspring and samples their types independently using a binomial distribution: $K\sim \mathrm{Binomial}(N,q)$, where $q=wp/(wp+1-p)$, then $p^{\prime }=K/N$. Its conditional expected value is $q$ and conditional variance is $q(1-q)/N$. With $w=1$, the conditional expected value remains $p$, yet realized frequencies wander through genetic drift. With $w>1$, genetic drift can still lower $p$ or lose $A$ entirely. Over several generations, the deterministic curve need not equal the mean stochastic trajectory, because the recurrence is nonlinear.

The Hardy-Weinberg principle describes the genotype proportions produced by independently combining two parental gametes from the same allele pool. For alleles $A,a$ with allele frequencies $p,q$, where $p+q=1$, the Hardy-Weinberg proportions areThe two ordered heterozygote outcomes contribute $pq+qp$. These proportions sum to $(p+q{)}^2=1$, and their implied $A$ allele frequency is $p^2+(2pq)/2=p$. Reshuffling copies therefore preserves their frequencies.

For stable Hardy-Weinberg equilibrium across generations, use a large population with random mating, equal allele frequencies in the parental pools, ordinary Mendelian segregation, and no natural selection, mutation, gene flow or genetic drift. Negligible genetic drift is an idealization of large population size. The one-generation random-mating equilibrium needs independent union of gametes; persistence additionally requires processes that change allele frequencies to be absent.

At $p=0.70,q=0.30$, the predicted genotype proportions are $0.49,0.42,0.09$. For 100 sampled individuals the expected values of the counts are $49,42,9$, rather than the $25,50,25$ of a specified $Aa\times Aa$ cross. The latter fixes both parental genotypes; Hardy-Weinberg proportions average over population mating.

An observed population need not have these proportions. For example, $70AA$, no $Aa$, and $30aa$ among 100 individuals still gives $p=0.70$, but has no heterozygotes. Random mating predicts $0.42$ heterozygosity in the next generation under the stated model. A finite sample can also depart by chance, so a discrepancy is a starting point for investigating inbreeding, population stratification or genotyping error, rather than proof of one cause.

Hardy's 1908 note made the conceptual payoff explicit: a dominant phenotype does not by itself make its allele spread. Complete dominance describes expression; changes in frequency require a separate mechanism such as natural selection. See Hardy, Mendelian Proportions in a Mixed Population, Science 28 (1908), 49–50.

Genetic recombination erodes linkage disequilibrium by reshuffling alleles between haplotypes, while preserving marginal allele frequencies in the idealized model. Let $f_{AB,t}$ denote the haplotype frequency in the gamete pool of generation $t$, let $p,q$ denote its marginal allele frequencies, and let $r\in [0,1/2]$ be the constant recombination fraction between two autosomal genetic loci. Assume random mating within one population, nonoverlapping generations, fair Mendelian segregation, equal reciprocal gamete production, the same $r$ in both sexes, and a deterministic population large enough to neglect genetic drift. Exclude natural selection, mutation and gene flow. ThenHere $p$ and $q$ remain constant: recombining changes which alleles occur together, not their total numbers. The other haplotype frequencies change by $+r{D}_t$ for $Ab,aB$ and $-r{D}_t$ for $ab$. Consequently the joint probability distribution stays within its allowed bounds and approaches the product of its marginal distributions when $r>0$.

For the original positive example $p=0.60,q=0.40,D_0=0.12$, a recombination fraction $r=0.10$ gives $D_1=0.108$ and $D_{10}\approx 0.04184$. The negative example $D_0=-0.12$ has exactly the opposite sequence: genetic recombination reduces the magnitude without reversing the sign. With $r=0$, the initial linkage disequilibrium persists in this model; even independent assortment at $r=1/2$ only halves $D$ each generation. A population association does not vanish after one generation merely because the genetic loci are unlinked.

The geometric sequence gives a generation scale: for $r>0$, the continuous interpolation reaches $|D_0|/2$ at $t_{1/2}=\log (1/2)/\log (1-r)$. Actual generations are integers; the first at or below half is $\lceil t_{1/2}\rceil$. For small $r$, $t_{1/2}\approx (\log 2)/r$. Because the denominator of allelic correlation stays fixed, $r_{\mathrm{LD},t}^2=(1-r{)}^{2t}{r}_{\mathrm{LD},0}^2$; its decay is faster than the magnitude of $D_t$.

A short derivation uses random mating to assign frequency $2f_{AB}{f}_{ab}$ to the double heterozygote $AB/ab$ and $2f_{Ab}{f}_{aB}$ to $Ab/aB$. Relative to $r=0$, each loses or gains an $AB$ gamete fraction $r/2$, respectively. Their net contribution is $r(f_{Ab}{f}_{aB}-f_{AB}{f}_{ab})=-rD$. The other genotypes produce no net change in $f_{AB}$, proving the recurrence and, by iteration, its solution.

This temporal model does not identify elapsed generations from physical separation alone. The recombination fraction depends on genetic recombination, including recombination hotspots, while population stratification, natural selection, gene flow or genetic drift can create or maintain linkage disequilibrium. Repeated sampling across time would also require accounting for sample size and haplotype estimation. A cross-sectional plot of association against distance between genetic loci is therefore a different observation from the idealized trajectory plotted here.

Natural selection involves systematic differences in expected reproductive success associated with variation in phenotype. For a continuing evolutionary response, the relevant variation must be inherited: offspring must tend to preserve the differences being favored. Actual offspring counts also fluctuate by chance, producing genetic drift. Survival matters insofar as it affects contribution to future generations. Genotype fitness is therefore contextual, depending on environment and the life cycle being counted, rather than simply measuring strength or longevity. See Clark, Douglas and Choi, Biology 2e: population evolution.

In chapter IV of his 1859 book, Darwin explained natural selection through preservation of favorable variations and greater chances of survival and reproduction. See Darwin, On the Origin of Species, first edition, chapter IV. The beautiful consequence is cumulative change without foresight: inherited differences are filtered repeatedly by present conditions. Natural selection neither generates a needed mutation on demand nor guarantees progress toward a predetermined goal.

Intracellular bacteria face a protected but highly controlled habitat: they must obtain host resources, avoid destruction, coordinate replication and reach the next host cell or host generation. Dependence evolves when access to that habitat is combined with loss or specialization of functions needed outside it. The selective pressures differ sharply between a pathogen spreading among hosts and a vertically inherited mutualist.

A first pressure is acquisition of a suitable intracellular niche. A bacterium entering by an uptake pathway may be delivered toward an acidified, degradative compartment. Chlamydia trachomatis instead replicates in an chlamydial inclusion whose interactions with host trafficking are modified by bacterial effectors. The niche supplies access to host-derived materials while keeping the replicating population spatially organized. Rickettsia prowazekii, in contrast, occupies the host cytoplasm after escape from an entry compartment. These are different solutions; it would be misleading to describe all intracellular bacteria as surviving by blocking lysosome fusion inside one generic vacuole.

The ability to manipulate trafficking and obtain metabolites becomes valuable. The Chlamydia trachomatis genome revealed a compact but specialized repertoire that includes transport and secretion-related functions, consistent with extensive host dependence.^[ref] genome annotation provides hypotheses about substrate use, not proof of actual import or enzyme activity. Test a candidate transporter's contribution using substrate uptake and intracellular growth, with genetic rescue and host-cell viability controls. A reduction in bacterial DNA can result from fewer entering bacteria, slower division or increased killing, so these stages should be measured separately.

Resource availability is selective rather than unlimited. A host cytoplasm may contain nutrients that are inaccessible across an chlamydial inclusion membrane, and an activated host can restrict metals or amino acids. Dependence on host supply favors uptake systems and scavenging, but makes the bacterium vulnerable to host nutrient restriction. The Rickettsia prowazekii genome showed extensive biosynthetic losses while retaining components of the citric acid cycle and respiration.^[ref] Thus an intracellular lifestyle does not mean universal loss of energy metabolism or reliance on one identical ATP-acquisition strategy. Imported metabolites and bacterial energy production can coexist.

Another pressure is coordination of the infectious and replicative phases. Chlamydia alternates between an elementary body, adapted for transmission and entry, and a reticulate body, adapted for growth inside the host. A durable infectious form can spend time outside a cell even when productive replication requires intracellular conditions. Developmental regulation therefore solves the conflict between surviving transfer and exploiting a nutrient-rich intracellular niche. Measure infectious progeny as well as genome copies: replication without successful differentiation need not yield transmission to another host cell.

Host defense creates additional pressures. Recognition of bacterial components can trigger inflammatory signaling, antimicrobial responses or death of the host cell before bacterial replication is complete. Pathogens can alter exposure, modify host signaling or delay destructive responses; the useful strategy depends on the niche and infection stage. Host-cell death may be harmful early yet assist release later, so “preventing apoptosis” cannot be treated as a universal permanent objective. To assign a bacterial effector's role, compare entry, replication, host response and release with an otherwise matched mutant, then restore the normal gene. Host toxicity caused by overexpressing an isolated effector is insufficient evidence for its physiological function during infection.

Some cytosolic species use the host cytoskeleton for movement and spread. Biochemical work on Rickettsia parkeri Sca2 demonstrated a formin-like actin assembly activity that helps explain the unbranched actin tails associated with spotted-fever rickettsiae.^[ref] Such movement can improve access to neighboring cells and reduce extracellular exposure. The result should not be generalized to every rickettsial species; different lineages have different motility repertoires. Inhibition of host actin also affects many host processes, so a bacterial separation-of-function mutant is more informative than a broad actin poison alone.

Long-term intracellular residence can produce genome reduction. Functions redundantly supplied by the host may become less strongly selected, permitting loss. Restricted transmission, small effective population size and limited recombination can also allow mildly deleterious changes to accumulate through genetic drift. Not every lost gene should be narrated as an adaptive improvement in efficiency. Compare closely related lineages, pseudogenes and host dependence; distinguish a dispensable ancestral pathway from a costly defect tolerated because the population cannot efficiently remove it.

Buchnera aphidicola illustrates a different outcome from pathogenic manipulation. It lives in aphid bacteriocytes and contributes to nutritional complementation in a host feeding on a restricted diet. genome sequencing showed strong reduction together with retention of substantial amino-acid biosynthetic capacity relevant to the association.^[ref] A vertically transmitted symbiont depends on host reproduction, which can favor cooperation and reliable transmission rather than maximal short-term host damage. Host and symbiont exchange still creates conflicts over resource allocation; stable association does not mean every trait benefits both partners equally.

Intracellular dependence should be tested operationally. Specialized host-cell-free culture allowed Coxiella burnetii, previously treated as obligatorily intracellular in culture, to grow under defined environmental and nutritional conditions.^[ref] It remains a highly specialized intracellular organism in its normal infection setting, but this experiment shows why failure in a standard laboratory medium is not proof that a living host is chemically indispensable. The main examples above concern established intracellular dependencies; Coxiella is a useful boundary case, not evidence that all such bacteria can be grown axenically.

Successful adaptation combines niche access, selective resource acquisition, control of host interactions and a transmission strategy. Over evolutionary time, those specializations can become irreversible dependencies through gene loss and population constraints, while the particular balance between host damage and cooperation follows the bacterium's route to the next generation.

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.

A point mutation changes an individual nucleotide position, introducing heritable sequence variation. Replication errors and spontaneous chemical damage are possible origins; mutation does not occur because a future change would be useful. Natural selection changes the frequencies of variants through their effects on reproduction, while genetic drift causes random frequency changes, especially in small populations. A sequence difference fixed between lineages is a substitution; many new mutations disappear and never become substitutions. Under a simple neutral diploid model, $2N_e\mu$ new neutral variants per generation times fixation probability $1/(2N_e)$ gives substitution rate $\mu$, motivating a molecular clock for suitable neutral sites.

To infer a phylogenetic tree, compare homologous, preferably orthologous sequences and align corresponding sites. A nucleotide substitution model converts observed sequence differences into estimated evolutionary branch lengths, allowing for multiple changes at one site. For example, in the equal-frequency, equal-rate Jukes-Cantor model, a site changes to any alternative with rate $\mu /3$, where $\mu$ is its total substitution rate. The probability of matching its initial state after elapsed time $t$ is $1/4+(3/4)e^{-4\mu t/3}$. Thus for an evolutionary separation with expected substitutions per site $d$, the observed mismatch fraction is $p=(3/4)(1-e^{-4d/3})$ andFor small $p$, $d\simeq p$; near saturation $p=3/4$, large evolutionary distances become difficult to distinguish. The correction does not rescue fully saturated data. Real analyses may allow unequal nucleotide frequencies, transition biases and variation of rates between sites.

Tree topology and branch lengths can be estimated using maximum likelihood estimation or Bayesian phylogenetic inference, with support assessed from resampling or posterior uncertainty. A gene tree need not equal the species tree, because incomplete lineage sorting, genetic recombination and horizontal gene transfer can produce different histories. Rooting with an appropriate outgroup establishes ancestral direction; pairwise differences alone do not supply a root.

Branch lengths measured in substitutions per site are not yet calendar times. Under a strict molecular clock with rate $\mu$ per site per year, a branch of duration $t$ has length $\mu t$; two contemporaneously sampled descendants separated at time $t$ have total path length $d=2\mu t$, hence $\boxed{t=d/(2\mu )}$. Calibration comes from independently dated constraints such as fossil calibration, geological information or sufficiently informative sampling dates. A dated common ancestor is not automatically identical in age to the earliest observed fossil. A relaxed molecular clock permits branch-specific rates while imposing a statistical model on them.

Without an external temporal calibration or an independently known rate, multiplying all times by a constant and dividing all rates by that constant leaves expected sequence divergence unchanged. Absolute dates are therefore unidentifiable from sequence divergence alone. A useful time-calibrated phylogenetic tree reports uncertainty from sequence sampling, substitution models, rate variation, topology and calibration bounds, rather than displaying an unjustified exact date for each ancestor.

Sequence divergence alone does not identify its cause. A difference can arise through mutation and spread through genetic drift, natural selection, or both. The useful comparison is with an explicit neutral expectation that accounts for shared ancestry, mutation rate and demographic history, rather than treating every conspicuous genomic difference as an adaptation.

First establish orthology, align homologous sequence accurately, and sample several individuals from each species. An apparent fixed difference in one pair of genomes may actually be segregating variation, an assembly error or comparison of paralogues. An outgroup can help polarize changes, although ancestral polymorphism and introgression can make different genomic regions have different histories. Estimate population size changes, migration and relatedness from many suitably chosen loci before interpreting an exceptional region.

For coding sequence, compare nonsynonymous and synonymous substitution rates per eligible site through the nonsynonymous-to-synonymous substitution rate ratio. A ratio greater than one in an appropriate branch or site model supports adaptive amino-acid substitution; a low average ratio indicates pervasive constraint but can hide positive selection at a few sites. A rise toward one can reflect relaxed purifying selection rather than adaptation. Synonymous sites are useful controls, not guaranteed neutral standards: codon usage, splicing regulation and mutation-context effects can matter.

Within-species polymorphism adds an independent timescale. In the McDonald-Kreitman test, let $D_N,D_S$ be nonsynonymous and synonymous fixed differences and $P_N,P_S$ their polymorphic counterparts. A simple neutral model predicts matching ratios. An excess $D_N/D_S>P_N/P_S$ is consistent with adaptive fixation, and the idealized estimatorestimates the adaptive fraction of nonsynonymous divergence. John H. McDonald and Martin Kreitman introduced this logic with the Drosophila Adh locus. Slightly deleterious variants inflate $P_N$, demographic changes affect which variants persist, and low counts make the ratio unstable; the model and uncertainty must therefore be examined rather than interpreting every positive estimate literally.^[ref]

Other signatures are complementary. A recent selective sweep can reduce linked diversity and distort the site frequency spectrum and linkage disequilibrium. Excess differentiation measured by FST can suggest geographically divergent selection, but can also result from reduced within-population diversity or population history. Background selection can mimic some diversity reductions. Comparisons with matched recombination and mutation environments, demographic simulations and several independent populations help discriminate these explanations. Long-maintained variation can instead leave signatures of balancing selection.

Noncoding differences require their own functional tests. A conserved regulatory element with an accelerated lineage-specific substitution rate is a candidate for altered function, but altered mutation bias or GC-biased gene conversion can also produce rapid sequence change. Reporter assays, allele replacement and measurements of expression, phenotype and fitness connect a candidate to a mechanism. Demonstrating an effect on expression does not by itself establish an advantage in the historical environment.

The strongest inference combines an appropriate population-genetic null, agreement among independent signatures, and a functional explanation connected to reproductive success. Conversely, failure to reject neutrality means that the evidence does not discriminate the models at the available resolution; it does not prove that a difference has no effect.

Visible divergence and reproductive incompatibility measure different outcomes of evolution. A few changes in a gene regulatory network can transform a conspicuous phenotype, whereas differences in the machinery of meiosis or biological development can prevent reproduction without greatly changing adult appearance. Neither membership of the same genus nor morphological similarity specifies divergence time or the extent of reproductive isolation.

First separate barriers. premating isolation can arise through differences in courtship, habitat or breeding season even when a forced cross produces fertile offspring. postzygotic isolation concerns the survival or fertility of the offspring themselves. Thus fertile captive hybrids do not refute substantial reproductive isolation in nature: they may never encounter suitable mates, may have intermediate ecological requirements, or may be unattractive to both parental biological populations. Under the biological species concept, the relevant issue is effective gene flow among natural biological populations, not the success of one laboratory mating. Conversely, a laboratory failure can reflect husbandry rather than an intrinsic incompatibility.

Rapid visible divergence is especially plausible when sexual selection changes plumage or other signals, or natural selection changes a feeding structure. Such changes can involve a small part of the genome, leaving most interactions required for fertility compatible. In Darwin's finches, genome comparisons associated an ALX1 region with beak shape and also revealed interspecific gene flow. This is an example of conspicuous adaptive differentiation coexisting with exchange of alleles, rather than evidence that all avian hybrids are fertile.^[ref] Continued gene flow can itself help maintain compatibility through much of the genome while divergent natural selection preserves differences at selected genetic loci.

For cryptic incompatibility, consider a Dobzhansky-Muller incompatibility. Two lineages start with compatible alleles $a,b$. One evolves $A$ in the presence of $b$, while the other evolves $B$ in the presence of $a$. Each change is tolerable in its own genetic background, but the combination $A,B$ disrupts function in hybrids. This requires neither lineage to pass through a low-fertility intermediate. Such epistasis can affect spermatogenesis, chromosome pairing, or interactions between nuclear and mitochondrial genomes, while leaving the adult body plan nearly unchanged. The relevant changes can arise through natural selection, genetic drift, or conflict over inheritance; adaptive divergence in external form is not required.

A particularly informative mammalian experiment concerns PRDM9, which helps specify recombination hotspots. In a particular house-mouse cross, changing its DNA-binding region to the corresponding human sequence restored chromosome pairing and male fertility. The experiment implicates compatibility between a genetic recombination regulator and its binding sites, rather than visible morphological divergence. Crucially, this was a rescue in a specified cross and genetic background, not a demonstration that PRDM9 explains all mammalian speciation.^[ref]

Comparative work on birds lacking PRDM9 found relatively stable recombination hotspots associated with functional genomic features.^[ref] This suggests a possible difference in the routes to incompatibility: turnover of PRDM9 and its targets is unavailable in those lineages. The inference is mechanistic, however; hotspot stability by itself does not demonstrate that avian hybrids will be fertile.

Differences in karyotypes provide another mechanism. A balanced chromosomal rearrangement can have little visible effect in either parent yet cause faulty pairing or unbalanced gametes in a heterozygous hybrid. However, a different chromosome count alone does not prove sterility: the arrangement and segregation behavior matter. Differences in genomic imprinting may also disturb placental development or embryonic growth in mammalian crosses, sometimes producing direction-dependent outcomes. Parent-of-origin growth abnormalities and disrupted genomic imprinting in particular rodent hybrids provide evidence that this mechanism can matter.^[ref] These are hypotheses to test in the relevant lineages, not a universal explanation for the bird–mammal comparison.

The sex of the offspring matters. Haldane's rule predicts that, when one sex of first-generation hybrids is absent, rare or sterile, it tends to be the heterogametic sex: typically XY males in mammals and ZW females in birds.^[ref] Recessive incompatibilities on a sex chromosome can be exposed when there is no corresponding second copy; rapid evolution of reproductive functions can contribute too. Therefore a report of fertile avian males is insufficient without testing females. The rule predicts a sex bias in a given cross, not that birds as a class have weaker reproductive isolation.

To discriminate these mechanisms, compare multiple independently diverged pairs, match their estimated divergence times using phylogenetics, and avoid treating taxonomic rank as a clock. Perform reciprocal crosses with parental within-lineage controls; measure fertilization, embryonic survival, both sexes' fertility, and second-generation and backcross success. A fertile first generation can conceal hybrid breakdown after genetic recombination disrupts coadapted combinations. Pair karyotype analysis and observations of meiosis with genome mapping and carefully controlled rescue experiments. Field measurements of mating and ancestry then determine whether compatibility translates into gene flow. The contrast is explained by different rates and genetic targets of visible divergence and reproductive divergence, with sex, cross direction, ancestry and assay conditions determining the observed outcome.

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.

DNA repair preserves usable inherited information, while genetic recombination can improve the combinations in which that information is transmitted. Their evolution need not have the production or prevention of diversity as an ultimate goal. natural selection favors consequences for survival and reproduction in particular contexts; novelty is neither automatically beneficial nor something an organism must maximize.

mutation creates new sequence variants; genetic recombination usually rearranges variants already present. A chromosomal crossover can produce a new haplotype without changing the allele frequency at either of two genetic loci. Conversely, DNA repair can restore sequence integrity after damage without eliminating the inherited differences between individuals. In particular, a normal alternative allele is not an error merely because it differs from another individual's sequence. These operations address different features of inheritance.

Unrepaired damage can block DNA replication or transcription, break a chromosome, or generate a harmful mutation. DNA repair thus has immediate benefits to a biological cell and can reduce the future burden of deleterious changes in descendants. In multicellular organisms, somatic mutation can compromise tissue function, while germline mutation changes what offspring inherit. The relevant balance includes the biochemical costs of DNA repair, the probability and severity of damage, and the possibility of DNA repair errors. natural selection for increased fidelity can become weak relative to genetic drift when an additional improvement is extremely small; perfect fidelity need not evolve. Some error-prone damage responses can also improve immediate survival even if surviving descendants carry more mutations.

The biochemical relationship is complementary. homologous recombination repairs a DNA double-strand break by using a related intact sequence as a template. In meiosis, programmed breaks can be repaired using a homologous chromosome, sometimes producing chromosomal crossovers and sometimes noncrossover products. Biochemical identification of Spo11 linked a component of the meiotic machinery to programmed break formation, demonstrating that meiotic diversification can begin with a controlled DNA-damage-and-repair process.^[ref] Crossovers also help establish the physical connections needed for reliable segregation in many meiotic systems. This immediate transmission function can favor the machinery independently of any later benefit from a novel haplotype. DNA repair with sister chromatids in another context can preserve sequence without the same exchange between parental homologues.

There are also population-level benefits of reshuffling. Suppose two favorable alleles, $A$ and $B$, occur in different lineages as $Ab$ and $aB$. With no exchange, their lineages compete and either favorable allele can be lost when the other lineage prevails. genetic recombination can create $AB$, allowing selection to combine their benefits; it can also separate a favorable allele from linked harmful changes. This is the core of the Fisher-Muller hypothesis and relates to Hill-Robertson interference, the reduction in the efficacy of natural selection caused by associations among selected genetic loci. Recreating low-load genotypes can also slow Muller's ratchet in finite biological populations.

The bookkeeping illustrates why this is not simply increasing all diversity. Let $x_{AB},x_{Ab},x_{aB},x_{ab}$ be haplotype frequencies and $D=x_{AB}{x}_{ab}-x_{Ab}{x}_{aB}$ be linkage disequilibrium. Under the standard random-mating model with recombination fraction $r$ and no other evolutionary forces,For $x_{Ab}=x_{aB}=1/2$ and the other two zero, $D=-1/4$. At $r=0.1$, the next frequencies are $x_{AB}^{\prime }=x_{ab}^{\prime }=0.025$ and $x_{Ab}^{\prime }=x_{aB}^{\prime }=0.475$. Each favorable allele still has frequency $1/2$: the process creates both the favorable combination and its unfavorable counterpart. Subsequent natural selection determines their fates. If the original associations are already advantageous, breaking them can instead impose a recombination load.

In controlled yeast experimental evolution, sexual and asexual treatments differed in their rates and genomic routes of adaptation, with the sexual treatment more efficiently separating beneficial from deleterious changes.^[ref] A mechanistic test should control the costs and bottlenecks of the sexual cycle, compare otherwise matched genetic backgrounds, and use reconstructed changes to distinguish genetic linkage from intrinsically different effects of a mutation. Such results establish a benefit under the tested population regime, not an unconditional advantage of sex in every environment. genetic recombination can disrupt favorable epistasis, and the costs of mating are distinct from those of crossover formation.

For an explanation of the evolution of the machinery itself, the association with its bearer matters. A recombination modifier must on average increase the representation of its own descendants, perhaps through improved segregation or through association with combinations that selection favors. Merely increasing long-run adaptability of the species is insufficient. Similarly, a mutator allele can spread in an asexual lineage by hitchhiking with a beneficial mutation it helped generate despite producing many harmful ones; exchange can weaken that association. The advantages of fidelity and of exchange consequently depend on genetic linkage, population size, environmental change and the current supply of beneficial variants.

DNA repair can regulate which exchanges are permitted as well. DNA mismatch repair corrects copying errors and can reject exchange between sufficiently diverged sequences. In bacteria, disabling mismatch-repair functions greatly relaxed a barrier to genetic recombination between divergent lineages.^[ref] An experiment separating these roles should compare DNA repair of copying errors with exchange between templates of controlled divergence, using a restored DNA repair gene as a rescue. DNA repair and genetic recombination can therefore be parts of the same evolved solution: preserve functional sequence, transmit chromosomes reliably, and reshuffle existing alleles when doing so improves reproductive success.

Let $u$ be the allele frequency of $A$, with random mating and sufficiently large population size to neglect genetic drift. The Hardy-Weinberg proportions before selection are $u^2,2u(1-u),(1-u{)}^2$. Write the three genotype fitnesses as $w_{AA},w_{Aa},w_{aa}$, with $w_{Aa}$ smaller than either homozygote value. After one viability-selection generation,ConsequentlyPut $a=w_{AA}-w_{Aa}>0$, $b=w_{aa}-w_{Aa}>0$. The underdominant allele-frequency dynamics has stable fixation states zero and one, separated by the unstable thresholdBelow it $A$ decreases to allele fixation of $a$; above it $A$ increases to fixation. Exactly at the threshold the deterministic frequency remains balanced, but any perturbation selects a side. This is underdominance, or heterozygote disadvantage, rather than balancing selection from heterozygote advantage.

In a continuous-time weak-selection or selection-rate convention, the corresponding equation is $\dot{u}=u(1-u)[(a+b)u-b]$. Equal homozygote fitness gives $\dot{u}=su(1-u)(2u-1)$ for a positive selection strength $s$. The next part uses this explicitly stated continuous-time convention; $s$ and migration then have the same time units.

Use a globally coupled underdominant metapopulation with equal habitat sizes:Symmetry preserves a common frequency $x$ in the initially $AA$ habitats and $y$ in the initially $aa$ habitats. With $q=m/N$, the equations reduce exactly toThe maximum positive selection rate is $f(u_m)=s/(6\sqrt{3})$ at $u_m=(3+\sqrt{3})/6$. Comparing it with the worst outward migration flux gives a sufficient barrier condition $\sigma \max (q,1-q)<f(u_m)/u_m$. A slightly stronger optimized condition follows from maximizing the per-capita restoring term $f(u)/u=s(1-u)(2u-1)$ over $u>1/2$: its maximum is $s/8$ at $u=3/4$.

Indeed, on $x=3/4$, with $0\le y\le 1/4$, $\dot{x}\ge 3s/32-(3/4)\sigma (1-q)>0$. On $y=1/4$, with $3/4\le x\le 1$, $\dot{y}\le -3s/32+(3/4)\sigma q<0$. The other two edges point inward as well. The invariant-rectangle criterion for underdominant coexistence therefore gives the explicit robust conditionsStarting from the specified pure habitats, the populations then remain in $x\ge 3/4$, $y\le 1/4$. Their global frequency lies between $3q/4$ and $q+(1-q)/4$, so neither allele fixes. This is a sufficient condition, not a claim that the bound is sharp. If $m=0$ or $m=N$, the population is already fixed; if $\sigma =0$, any mixture of the two pure habitat types persists.

The exactly balanced case illustrates why robustness and the prescribed trajectory differ. When $q=1/2$, symmetry forces $y=1-x$ andFor $\sigma <s/2$ its heterogeneous limiting state isHere $f^{\prime }(x_{\ast })=f^{\prime }(y_{\ast })=-s+3\sigma$, so the full Jacobian matrix has diagonal entries $f^{\prime }(x_{\ast })-\sigma /2$ and off-diagonal entries $\sigma /2$. Its two full-system eigenvalues are $-s+3\sigma$ and $-s+2\sigma$: the symmetric underdominant two-cluster equilibrium is robustly attracting only for $\sigma <s/3$. For larger migration the exactly symmetric initial condition can still avoid fixation, eventually approaching the spatially uniform frequency $1/2$ when $\sigma \ge s/2$, but the symmetry-broken mean-frequency mode is unstable. A finite population with genetic drift need not maintain deterministic coexistence indefinitely.

Assume a closed population with even size $N$, unrelated founders, monogamous mating, nonoverlapping generations and exactly two offspring per couple. Sex labels are exchangeable for the autosomal calculation, with enough individuals of each sex to form the stipulated pairs. A randomly chosen distinct offspring pair are siblings with probability $1/(N-1)$: among the $\left(\genfrac{}{}{0ex}{}{N}{2}\right)$ pairs, exactly $N/2$ are within-family pairs. This fixed-family-size fact, rather than independent sampling of all parents, determines the recurrence.

First give the elementary chromosome-sharing idealization consistent with taking self-relatedness as one and ordinary sibling relatedness as $1/2$. In this model each individual carries one ancestral label, passed through one randomly chosen parent. Let $r_G$ be mean sharing between distinct individuals, and let $m_G$ be the mean sharing between mates. Siblings share their parental label with probability $(1+m_G)/2$. The average relatedness between parents from different couples isTo see the latter count, the ordered pairs from different couples are all ordered distinct pairs minus the $N$ ordered mating pairs. Each pair of parents from different couples contributes equally to the four possible cross-family offspring relationships. Combining sibling and nonsibling offspring givesWith random mating, $m_G=r_G$ in expectation, so $r_0=0$ givesThe mean sharing tends to one, on a generation scale of order $2(N-1)$. It is not linear indefinitely: repeated ancestry saturates the probability. This is an exact result for the stated one-label approximation, not an exact diploid formula after substantial inbreeding has developed.

Now prohibit full-sibling matings, with $N\ge 4$ and representative sampling of allowed nonsibling mating pairs. The following closed recurrences assume that opposite-sex mating pairs have the same mean autosomal sharing as all nonsibling pairs; a sex-specific or assortative scheme requires separate mating means. In the exchangeable one-label model write $u_G$ for mean relatedness of a nonsibling pair, so $m_G=u_G$. The two required recurrences areThe second relation includes sharing through grandparents and more distant ancestors; setting $u_G=0$ forever would be wrong. For $v_G=1-r_G$, $w_G=1-u_G$, the linear recurrence matrix isIts eigenvalues areboth of magnitude less than one for $N\ge 4$. Thus $r_G\to 1$ even when siblings never mate. Explicitly, $1-r_G=A{\lambda }_{+}^G+(1-A){\lambda }_{-}^G$, where $A=[1-1/(2(N-1))-{\lambda }_{-}]/({\lambda }_{+}-{\lambda }_{-})$. Sibling avoidance prevents immediate close-relative matings, but it does not make the remaining population unrelated. In a fixed-two-offspring model it need not reduce average pairwise sharing at every generation: a smaller mating mean actually increases the displayed increment in $r_G$.

For a fully diploid formulation, self-kinship also changes with inbreeding and cannot remain artificially fixed. Let ${\Phi }_G$ be mean kinship coefficient between distinct individuals, $F_G$ mean inbreeding coefficient, and ${\Phi }_{m,G}$ mean kinship of mates. A sibling pair has kinship $(1+F_G)/4+{\Phi }_{m,G}/2$, because parental self-kinship is $(1+F_G)/2$. Repeating the pair count yields finite-population chromosomal coancestry:Initially ${\Phi }_0=F_0=0$. For random mating take ${\Phi }_{m,G}={\Phi }_G$. An explicit solution is the affine matrix recurrenceIts characteristic polynomial is ${\lambda }^2-[1-1/(2(N-1))]\lambda -1/[4(N-1)]$; its two roots have magnitude less than one. Hence ${\Phi }_G,F_G\to 1$. The normalized relationship $2{\Phi }_G/(1+F_G)$ tends to one, while unnormalized additive relationship $2{\Phi }_G$ tends to two, demonstrating why that coefficient cannot always be read literally as a probability.

With sibling avoidance replace ${\Phi }_{m,G}$ by the nonsibling mean $U_G$ and addThese equations compute every generation's kinship and inbreeding, including delayed consanguinity through cousins. Under random allowed mating in a connected closed population, backward ancestral gene-copy lineages have a positive chance to meet in each sufficiently long block of generations; repeated blocks force eventual coalescence. Thus sibling avoidance alone does not prevent long-term identity by descent and genetic drift. Maintaining distinct isolated breeding groups, introducing immigrants, or changing offspring-number variance gives a different model. The elementary exponential formula answers the fixed-self-sharing approximation; the diploid recurrences give the consistent inbreeding-aware answer.

Linkage disequilibrium means that alleles at different genetic loci are associated in population haplotypes. For two biallelic loci,The first quantity measures departure from independent allele frequencies; the second is the squared allelic correlation. Standardized linkage disequilibrium $D^{\prime }$ instead compares $D$ to its possible range at the given frequencies. High $|D^{\prime }|$ can coexist with low $r^2$, so it need not make one genetic marker a good predictor of another. Genetic linkage is a different notion: it describes within-family transmission and the recombination fraction.

The chromosome-21 study of Patil and colleagues, Science 2001 found consecutive single-nucleotide polymorphisms organized into regions with comparatively few common haplotypes; typically three haplotypes accounted for over 80% of the sampled chromosomes within a block. Such haplotype blocks can be represented efficiently by tag SNPs. Strong association within blocks and weaker association across boundaries give a useful description, though block boundaries depend on allele frequencies, population history, sampling and the operational definition of a block. They do not divide chromosomes into regions where genetic recombination is respectively impossible and certain.

The ancestry explanation uses an ancestral recombination graph. Going backwards, ancestral lineages coalesce, while a genetic recombination event splits one lineage's ancestral material into contributions from two parents. At a fixed genomic position, retaining only its ancestry gives a phylogenetic tree with the appropriate coalescent law. Nearby positions often retain much of the same tree; mutations on shared branches then produce correlated allele states. Breakpoints can change the local tree, weakening that shared history.

Under a simple random-mating two-locus model, linkage-disequilibrium decay under recombination gives $D_{t+1}=(1-r)D_t$. Actual populations also undergo genetic drift, population bottlenecks, mutation and migration, which can create or preserve association. Recombination rates themselves vary along chromosomes: McVean and colleagues, Science 2004 inferred strong local variation from sequence data. These observations motivate an ancestry model with spatially varying recombination, rather than treating physical distance as the only determinant of linkage disequilibrium.

For fine mapping, a typed genetic marker can show genetic association because it is correlated with a nearby causal genetic variant. Linkage-disequilibrium marker tagging makes a sparse first-stage scan possible, but several variants in one block can produce nearly indistinguishable association patterns. Historical recombination and diverse local genealogies can narrow the interval when denser genotyping is available.

An ancestral recombination graph turns that intuition into a likelihood calculation. One can place a proposed disease-origin event on a local ancestral branch and integrate over possible graphs, mutation histories and unobserved haplotypes. The local phylogenetic tree then relates disease-bearing chromosomes and predicts which markers should remain associated. In a Bayesian inference, uncertainty in ancestry, mutation location and demographic parameters enters the posterior distribution rather than being replaced by one guessed tree. Whole-graph inference is computationally demanding, and approximate or pairwise methods often supply practical estimates.

The useful distinction is between finding an associated region and identifying its causal variant. Strong linkage disequilibrium improves marker coverage but can limit resolution inside the associated region. Population stratification, ascertainment, rare variants and misspecified population history can also mislead an association analysis. Combining fine mapping with independent biological evidence is therefore stronger than interpreting an association peak as proof of causation.

Five consequences of the ancestral recombination graph are:

Shared ancestry explains correlation; recombination makes ancestry a mosaic; mutation on the resulting branches converts genealogy into observed variation. These explanations use the neutral ARG, with a variable recombination map where stated. Selection, population subdivision or changing population size can also be studied genealogically, but require extending the event model rather than being consequences of the constant-size panmictic neutral process alone.