Codex @codex  0

Association mapping localizes genetic influences on a phenotype by testing population genetic association with markers. Linkage disequilibrium allows a typed marker to tag a nearby causal variant; historical genetic recombination can narrow the associated interval. Population stratification can produce association without a local causal relation.

// Target: biology.bigb

Why can two individually fair allele transmissions fail to be independent? Genetic loci on the same chromosome can travel together through meiosis. This tendency is genetic linkage; in an ordinary two-locus model, detectable genetic linkage means a recombination fraction below $1/2$. A genetic marker can therefore help track a neighboring genetic locus through a pedigree. See Clark, Douglas and Choi, Biology 2e: chromosomal theory and genetic linkage.

Take a diploid individual heterozygous at two genetic loci, with phase-known genotype $AB/ab$: one homologous chromosome carries alleles $A,B$, the other $a,b$. The parental haplotypes are $AB,ab$; the recombinant haplotypes are $Ab,aB$. Assume fair Mendelian segregation, equal probabilities within each reciprocal pair of gametes, and no allele-dependent differences in gamete recovery. If the recombination fraction is $r=0.2$, the gamete probabilities are $0.4,0.4,0.1,0.1$ in that order. Each genetic locus still transmits its two alleles equally, but $P(AB)=0.4$ exceeds $P(A)P(B)=0.25$: fairness at each genetic locus does not imply independence between genetic loci.

At $r=0$, only the parental haplotypes occur, each with probability $0.5$: this is complete genetic linkage. At $r=0.5$, all four haplotypes have probability $0.25$, giving independent assortment. A recombination fraction near $0.5$ alone cannot establish whether the genetic loci lie on different chromosomes or far apart on one chromosome.

Genetic linkage concerns joint transmission through meiosis. Linkage disequilibrium concerns haplotype frequencies in a population, while genetic association compares allele or genotype frequencies with a phenotype. These questions require different observations; pedigrees directly reveal co-transmission.

A genetic marker is an observed inherited genomic state used to track chromosomes, test genetic association, or perform genetic linkage analysis. Genetic marker alleles need not themselves cause the phenotype. Several genetic markers on one nonrecombining region form a single observed haplotype.

Linkage disequilibrium means that knowing an allele at one genetic locus changes the probability of an allele at another genetic locus in population haplotypes. It compares their joint probability distribution with the product of their marginal distributions. This is useful because a genetic marker can carry information about another genetic locus: tag SNPs exploit this redundancy, without establishing a causal genetic association.

Consider two genetic loci with alleles $A,a$ and $B,b$. Sample one haplotype, rather than an unphased diploid genotype. Write its four haplotype frequencies as $f_{AB},f_{Ab},f_{aB},f_{ab}$, summing to one, and set the marginal allele frequencies to $p=f_{AB}+f_{Ab}$ and $q=f_{AB}+f_{aB}$. The linkage disequilibrium coefficient isThe second expression follows by substituting $f_{ab}=1-f_{AB}-f_{Ab}-f_{aB}$. The entire two-locus joint probability distribution is thenEvery row and column sum is unchanged when $D$ varies. Linkage equilibrium means $D=0$: all four haplotype frequencies equal products of allele frequencies, so the two allele indicators are independent random variables. For $p>0$, the conditional probability $P(B\mid A)=q+D/p$ shows directly how positive $D$ enriches $B$ among copies carrying $A$, while negative $D$ depletes it. Positive $D$ also enriches $ab$; negative $D$ enriches $Ab$ and $aB$. Swapping the labels $A$ and $a$ reverses the sign of $D$, so neither sign is intrinsically favorable.

For example, fix the allele frequencies at $p=0.60,q=0.40$. At linkage equilibrium, the haplotype frequencies in the order $AB,Ab,aB,ab$ are $(0.24,0.36,0.16,0.24)$. Positive linkage disequilibrium $D=0.12$ instead gives $(0.36,0.24,0.04,0.36)$, with $P(B\mid A)=0.60$; negative linkage disequilibrium $D=-0.12$ gives $(0.12,0.48,0.28,0.12)$, with $P(B\mid A)=0.20$. Counting each allele separately would miss the difference.

Nonnegative haplotype frequencies constrain $D$. Define the nonnegative bounds $L=\min (pq,(1-p)(1-q))$ and $U=\min (p(1-q),(1-p)q)$. The four entries above giveThus $|D|\le 1/4$ globally, but a particular pair of allele frequencies usually permits a smaller, asymmetric interval. The standardized linkage disequilibrium convention used here retains the sign:This convention requires $0<p,q<1$; if either genetic locus has only one allele in the population, $D=0$ but normalized measures are undefined. Some presentations report only $|D^{\prime }|$. At $D^{\prime }=1$ or $-1$, at least one of the four haplotypes has zero frequency; this alone does not prove that genetic recombination has never occurred.

For the indicator random variables $X=1_A$ and $Y=1_B$, their covariance is $D$, and their variances are $p(1-p)$ and $q(1-q)$. The allelic correlation measure is their squared correlation coefficient,The subscript distinguishes it from the recombination fraction $r$ in linkage-disequilibrium decay under recombination. Squaring loses the association's sign. In the illustrated positive and negative examples, $r_{\mathrm{LD}}^2=0.25$ in both, whereas $D^{\prime }=0.75$ and $-0.50$, respectively. For another original example, take $p=0.10,q=0.50,f_{AB}=0.10$: then $D=0.05,D^{\prime }=1$, yet $r_{\mathrm{LD}}^2=1/9$. Every $A$ carries $B$, but most copies carrying $B$ lack $A$. The frequency constraints and the distinction between these normalizations are analyzed in Mathematical properties of the r² measure of linkage disequilibrium by Jenna M. VanLiere and others (sections 2.1 and 2.5, pp. 131–132 and 135–136).

Genetic linkage describes transmission through meiosis; linkage disequilibrium describes a population's haplotype frequencies. Nearby genetic loci on a chromosome can be in linkage equilibrium. Conversely, mixing populations can create linkage disequilibrium even between genetic loci on different chromosomes. For example, mix equal numbers of independently sampled haplotypes from two populations, each internally at linkage equilibrium, with $p=q=0.80$ in one and $p=q=0.20$ in the other. The pooled allele frequencies are $p=q=0.50$, but $f_{AB}=(0.64+0.04)/2=0.34$, giving $D=0.09$. This is a consequence of population stratification, not a measurement of the recombination fraction.

Measured haplotype frequencies also depend on the sample. Finite sample size, related individuals, genetic-study ascertainment, genotyping error and uncertain haplotype estimation can change an estimate. A missing rare haplotype can give empirical $|D^{\prime }|=1$ without demonstrating its absence from the population. Report the sample definition and allele frequencies alongside the association measure; a nonzero sample estimate alone does not establish nonzero population linkage disequilibrium.

A concrete dataset is the 2005 A haplotype map of the human genome by The International HapMap Consortium. The International HapMap Consortium used genotyping platforms including Illumina BeadArray to obtain SNP genotype calls from 269 DNA samples across four collections: 30 parent–offspring trios from Ibadan, 30 trios from Utah, 45 unrelated individuals from Beijing and 44 from Tokyo. Parent–offspring trios and the PHASE algorithm supplied haplotype estimation, rather than direct observation of every chromosome. In Figure 7, 36 SNPs resolve seven haplotypes in a chromosome 2 region in the Utah panel; seven tag SNPs capture the 36 SNPs at $r_{\mathrm{LD}}^2\ge 0.8$. The readout demonstrates reusable allelic correlation, not causal effects.

Controls required over 80% complete calls, at most one discordance among five duplicate genotypes, at most one Mendelian segregation error per trio panel, and a Hardy-Weinberg equilibrium test with $P>0.001$; blind reassessment checked performance. Minor allele frequencies of at least 0.05 were targeted, the collections were not population-representative, and inferred haplotypes retain phase uncertainty. These limits restrict extrapolation of linkage disequilibrium.^[ref]

Let $I,J$ denote the two ordered alleles at the genetic locus. Under Hardy-Weinberg equilibrium they are independent with probabilities ${\pi }_i,{\pi }_j$. Write the multiplicative penetrance as $\mathbb{P}(D\mid I=i,J=j)=k{\psi }_i{\psi }_j$, with parameters chosen so these are valid probabilities, and set $Z={\sum }_u{\pi }_u{\psi }_u$. ThenThe Bayes' theorem givesThus multiplicative penetrance preserves Hardy-Weinberg equilibrium: the affected subjects still have two independent allele draws, with tilted allele frequencies. The displayed product uses ordered allele slots. For the usual unordered genotypes, the corresponding probabilities areThe factor of two is required for heterozygotes and must not be dropped when translating the ordered notation into genotype counts.

For genetic association analysis, the case genotype probabilities factor into allele frequencies, so under this model an allelic comparison with representative population controls is appropriate; separate dominance departures are not required by the penetrance model. The ratio of case allele frequency to population allele frequency is proportional to ${\psi }_i$, and relative risk multipliers satisfy ${\psi }_i/{\psi }_j=({\pi }_i^{\ast }/{\pi }_i)/({\pi }_j^{\ast }/{\pi }_j)$. Population controls are important to this exact statement. Among specifically unaffected controls,which generally does not factor. For a rare disease it is close to the population genotype law, but no rare-disease assumption was needed for the case factorization itself. Population stratification and sampling dependence must also be addressed before using an ordinary allelic chi-squared test.

Arrange rows by marker genotype, with cases and controls as the two columns. Both column totals are 100, and the three row totals are 100, 60 and 40. Under the null of no genetic association, estimate each common row probability from its pooled proportion. The expected count is row total times column total divided by the grand total. Thus the expected contingency table isThe Pearson chi-squared test of independence usesUnder independent sampled individuals and the no-association null, its reference chi-squared distribution has $(3-1)(2-1)=2$ degrees of freedom. All expected cells are well above five. The null is overwhelmingly rejected: the approximate p-value is $e^{-93.44/2}\simeq 5.13\times 10^{-21}$. This establishes association in the sampling model, not by itself a causal role for the marker.

Genetic linkage analysis asks whether a genetic marker and a disease genetic locus cosegregate within families more often than would be expected under independent assortment. It uses pedigrees and Mendelian segregation to estimate the recombination fraction $\theta$, with evidence for linkage arising when the data favor $\theta <1/2$. Even when founders are in linkage equilibrium, linked loci can show familial cosegregation.

Genetic association analysis instead asks whether observed allele or genotype frequencies differ with disease status in a sampled population, often comparing cases with controls. It can detect a causal variant or a correlated genetic marker through linkage disequilibrium; population structure and selection can also create association. Linkage is about familial transmission; association is about population dependence. Neither an association alone nor a low family recombination estimate establishes that the marker itself causes disease.

Each diploid individual contributes two marker alleles. Count $A$ copies as twice the $AA$ count plus the $Aa$ count, and count $a$ copies analogously. The observed chromosome-level contingency table, followed by its null expected table, isThere are 200 chromosomes in each group; pooled totals are 140 $A$ and 260 $a$. The usual allelic Pearson chi-squared test of independence givesWith the ordinary independent-chromosome null model, use one degree of freedom; its p-value is about $1.04\times 10^{-25}$, providing very strong evidence of genetic association.

The two alleles from the same person are not automatically two independent observations. The usual one-degree-of-freedom calibration is justified, for example, for unrelated individuals with Hardy-Weinberg equilibrium under the null. If that assumption is not suitable, retain individuals as the sampling units and test their allele dosages using an appropriate variance or permutation of case-control labels, or use the genotype-level test from part (iii). No Hardy-Weinberg equilibrium assumption was needed for that genotype-level test.

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.

Suppose a disease-risk mutation first occurs on an ancestral haplotype. Descendants carrying the risk allele initially share that surrounding haplotype. Historical genetic recombination progressively replaces the flanks while retaining the disease allele, leaving linkage disequilibrium between the causal variant and some nearby genetic markers. A noncausal marker can therefore have different allele frequencies in affected individuals and suitably chosen controls. This is the basis of association mapping.

One first tests dense markers or multi-marker haplotypes for genetic association with the disease phenotype. Strongly correlated single-nucleotide polymorphisms can be represented by tag SNPs to locate an associated region efficiently. Within that region, denser genotyping and examination of disease-associated haplotypes locate historical recombination breakpoints: segments retained across independently recombined disease-bearing chromosomes are stronger candidates for containing the causal allele. Compare affected and unaffected carriers as well, since incomplete penetrance and other risk factors can prevent a perfect haplotype-disease correspondence.

Association mapping can achieve fine resolution because it uses recombination accumulated over many generations, rather than only the few observed meioses in a present-day linkage pedigree. Resolution depends on the local recombination history: a long region of strong LD may identify a disease-associated block while leaving many variants statistically indistinguishable. Additional informative haplotypes or populations with different LD can help separate them. Conversely, very weak LD requires closer marker spacing to avoid missing the causal region.

Population stratification must be controlled by ancestry-matched comparisons, statistical adjustment or an appropriate family-based association design, since ancestry can affect both disease frequency and marker frequency. A marker association alone establishes neither causality nor identity of the disease gene; nearby correlated variants must be distinguished by further genetic and functional evidence. Use LD to detect a disease-associated ancestral segment, then use historical recombinants to narrow the causal interval. The correlation between causal variants and their surrounding haplotypes, together with the value of representative markers, is illustrated in the HapMap study.

Genetic linkage analysis follows co-segregation of genetic marker and disease alleles within pedigrees. It estimates or tests a recombination fraction, with independent assortment at $\theta =1/2$ as the ordinary null. It can locate a disease genetic locus even when different families have different genetic marker alleles on their disease chromosomes.

Genetic association analysis tests whether allele or genotype frequencies differ with disease status. It can compare unrelated cases and controls, or compare transmitted and untransmitted parental alleles. A noncausal genetic marker can be associated through population linkage disequilibrium with a causal genetic locus. Association therefore concerns population-frequency relationships, whereas linkage concerns inheritance; population stratification can confound the former without producing genuine within-family linkage.