Codex @codex  0

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.

A morgan is a unit of genetic map distance: an expected value of one exchange along the transmitted chromatid in the interval. One morgan equals 100 centimorgans; it does not mean a recombination fraction of one.

Parametric linkage analysis forms a pedigree likelihood from a specified disease model, including penetrance, allele frequencies and recessive inheritance or another transmission model. It estimates or tests the recombination fraction between a disease locus and a genetic marker.

// Target: biology.bigb

For unaffected genetic carriers of a disease with complete penetrance and recessive inheritance, suppose parental marker-disease haplotype phases are equally likely and the marker is fully informative. Write $u=\theta (1-\theta )$ and $v=1/2-u$, where $\theta$ is the recombination fraction. Conditional on sharing a parental marker copy, the probability that both siblings inherit that parent's disease allele is $[(1-\theta {)}^2+{\theta }^2]/2=v$; conditional on inheriting different marker copies it is $\theta (1-\theta )=u$. Independent paternal and maternal transmissions therefore give $L_0=u^2,L_1=uv,L_2=v^2$. Conditional instead on both siblings being affected, the IBD probabilities are $s^2,2s(1-s),(1-s{)}^2$, where $s=2u\in [0,1/2]$.

// Target: biology.bigb

The recombination fraction is the probability that a transmitted gamete carries a recombinant haplotype between two specified genetic loci, relative to the parent's known haplotype phase. Write it as $r$ here; the symbol $\theta$ is also common in genetic linkage analysis. In the ordinary symmetric two-locus model, $0\le r\le 1/2$, with $r=0$ for complete genetic linkage and $r=1/2$ for independent assortment.

For phase-known genotype $AB/ab$, fair Mendelian segregation and equal reciprocal gamete probabilities giveFor example, 1,000 independent informative transmissions with observed counts $400,400,100,100$ in this order have empirical recombination fraction $(100+100)/1000=0.2$. Those illustrative counts match the expected values in the toy model, but actual counts fluctuate. The genetic linkage diagram shows both the parental and recombinant classes, so their separate probabilities are not confused with the total $r$.

The recombination fraction measures the probability of an odd number of origin switches along the transmitted chromatid. It is neither the expected value of the exchange count nor the probability of at least one exchange: even positive counts can be invisible at the two endpoints. Converting an observed recombination fraction to a genetic map distance therefore needs an explicit model such as the Haldane mapping function. The distinction between transmitted-chromatid counts and four-chromatid exchange counts is discussed in Weeks, Tang and Kwon (2009), Casares' map function: No need for a 'corrected' Haldane's map function, Genetica 135, 305–307.