Heterozygote Advantage

Heterozygote advantage (also termed overdominance or superdominance) is a mode of natural selection in which diploid individuals carrying two distinct alleles at a specific locus (A₁A₂) achieve higher lifetime Darwinian fitness (w₁₂) than individuals homozygous for either allele (w₁₂ > w₁₁ and w₁₂ > w₂₂). Within the mathematical framework of population genetics, first solved by Ronald Fisher in 1922, heterozygote advantage represents one of the primary mechanisms of balancing selection: instead of driving one allele to fixation and eliminating the alternative allele, overdominant selection actively protects both alleles from loss and converges toward a globally stable internal polymorphic equilibrium (p̂ ∈ (0, 1)).

Because Mendelian segregation during meiosis continuously regenerates less-fit homozygotes (A₁A₁ and A₂A₂) from matings between superior heterozygotes (A₁A₂ × A₁A₂ → 1/4 A₁A₁ + 1/2 A₁A₂ + 1/4 A₂A₂), a population maintained at overdominant equilibrium incurs a permanent reduction in mean fitness known as segregation load. Empirically, heterozygote advantage explains why certain recessive lethal or debilitating genetic disorders—most famously the sickle-cell hemoglobin allele (HbS) in malarial regions—persist at high frequencies across human and wild animal populations.

Interactive Simulation: Set relative genotype fitnesses so that w₁₂ > w₁₁ and w₁₂ > w₂₂ to observe trajectories converging to the overdominant equilibrium in the Allele Drift Simulator.

Selection Coefficients (s, t) and Allele Frequency Dynamics

Consider a diploid autosomal locus with two alleles, A₁ and A₂, occurring at frequencies p and q = 1 − p in a randomly mating population (Hardy-Weinberg principle). Normalizing the fitness of the superior heterozygote A₁A₂ to unity (w₁₂ = 1), the reduced fitnesses of the two homozygotes are parameterized by positive selection coefficients s and t (where 0 < s ≤ 1 and 0 < t ≤ 1):

• Genotype A₁A₁ (frequency p²): Relative fitness w₁₁ = 1 − s
• Genotype A₁A₂ (frequency 2pq): Relative fitness w₁₂ = 1
• Genotype A₂A₂ (frequency q²): Relative fitness w₂₂ = 1 − t

The marginal fitness of allele A₁ is w₁ = pw₁₁ + qw₁₂ = 1 − sp, and the marginal fitness of allele A₂ is w₂ = pw₁₂ + qw₂₂ = 1 − tq. Mean population fitness (w̄) across all zygotes is:

w̄ = p²(1 − s) + 2pq(1) + q²(1 − t) = 1 − sp² − tq²

Applying Wright's general selection equation, the single-generation change in the allele frequency of A₁ is:

Δp = [ pq · (w₁ − w₂) ] / w̄ = [ pq · (tq − sp) ] / (1 − sp² − tq²)

Equilibrium Frequencies (p̂, q̂) and Global Stability

Setting Δp = 0 reveals two trivial boundary equilibria at p = 0 (fixation of A₂) and p = 1 (fixation of A₁), which are both unstable whenever s > 0 and t > 0, plus a unique internal polymorphic equilibrium where the marginal fitnesses of the two alleles are equal (w₁ = w₂ ⇒ tq̂ = sp̂). Substituting q̂ = 1 − p̂ yields the equilibrium allele frequencies:

p̂ = t / (s + t) and q̂ = s / (s + t)

Because Δp > 0 whenever 0 < p < p̂ (as tq > sp) and Δp < 0 whenever p̂ < p < 1 (as tq < sp), this internal equilibrium is globally stable in an infinite population: any initial polymorphism converges toward p̂. Furthermore, in accordance with Fisher's fundamental theorem of natural selection, p̂ maximizes mean population fitness (dw̄/dp = −2sp + 2t(1 − p) = 0 at p = p̂). Substituting p̂ and q̂ into the expression for w̄ gives the maximum equilibrium mean fitness (w̄_max) and the resulting segregation load (L_seg = 1 − w̄_max):

w̄_max = 1 − [ s · t / (s + t) ] ⇒ L_seg = ( s · t ) / (s + t)

Even at optimal equilibrium, the population loses a fraction st / (s + t) of its reproductive output each generation due to the production of homozygotes.

Sickle-Cell Hemoglobin (HbS) and Plasmodium falciparum Malaria

The canonical molecular demonstration of heterozygote advantage in humans was established by Anthony C. Allison in 1954 across East Africa. The sickle-cell allele (HbS) arises from a single nucleotide point mutation (A → T) in the sixth codon of the HBB β-globin gene on chromosome 11, substituting hydrophobic valine for hydrophilic glutamic acid (Glu6Val).

In individuals homozygous for the mutant allele (HbS/HbS), deoxygenated hemoglobin polymerizes into rigid fibers that deform erythrocytes into a sickle shape, causing vaso-occlusive crises, hemolytic anemia, and a pre-modern fitness deficit of t ≈ 0.80–0.90. Conversely, in holoendemic malarial regions, wild-type homozygotes (HbA/HbA) suffer childhood mortality and severe cerebral or placental malaria caused by the protozoan parasite Plasmodium falciparum, incurring a selective disadvantage of s ≈ 0.10–0.18. Heterozygous carriers (HbA/HbS, sickle-cell trait) experience minimal clinical anemia because normal HbA inhibits polymer elongation under physiological oxygen tensions, yet parasitized HbA/HbS red blood cells undergo premature sickling under the low-pH, low-O₂ intracellular environment induced by the parasite. These infected cells are rapidly cleared by macrophages in the spleen and display impaired parasite cytoadherence via PfEMP1 surface proteins, conferring an ~80–90% protection against life-threatening malaria. Substituting s = 0.15 and t = 0.85 yields the exact equilibrium frequency q̂ = 0.15 / (0.15 + 0.85) = 0.15 observed across equatorial West and Central Africa.

Other Human, Animal, and Agricultural Examples

• Hemoglobinopathies and Enzymopathies: Parallel malarial heterozygote advantage maintains Hemoglobin C (HbC), α- and β-thalassemias, Southeast Asian ovalocytosis (SLC4A1), and X-linked glucose-6-phosphate dehydrogenase (G6PD) deficiency across tropical populations.
• Major Histocompatibility Complex (MHC / HLA): Vertebrate MHC Class I and Class II glycoproteins bind intracellular and extracellular pathogen peptides for presentation to T lymphocytes. Because MHC expression is codominant, a heterozygous individual (MHC¹/MHC²) expresses two distinct sets of peptide-binding grooves (divergent allele advantage, Wakeland et al., 1990), enabling recognition of a broader repertoire of viral and bacterial epitopes than either homozygote—verified empirically in human HIV-1 progression cohorts and salmonid infectious hematopoietic necrosis virus challenges.
• Warfarin Resistance in Brown Rats (Rattus norvegicus): Anticoagulant rodenticides such as warfarin inhibit vitamin K epoxide reductase (VKORC1). In areas treated with warfarin, wild-type homozygotes (RR) die from hemorrhage (s = 0.37), whereas homozygous resistant mutants (SS) suffer a 20-fold increase in dietary vitamin K requirements and die from spontaneous vitamin K deficiency (t = 0.63); heterozygotes (RS) survive both warfarin exposure and normal dietary vitamin K levels, stabilizing the mutant allele at q̂ ≈ 0.37.
• Single-Gene Overdominance in Tomato Yield (SFT): In domesticated tomatoes (Solanum lycopersicum), Zachary Lippman and Dani Zamir (2010) demonstrated that plants heterozygous for a loss-of-function mutation in the flowering gene SINGLE FLOWER TRUSS (SFT/sft) exhibit up to a 60% increase in fruit yield compared to either homozygote by optimizing the balance between vegetative sympodial shoots and reproductive inflorescences.

Distinguishing True Overdominance from Pseudo-Overdominance

In genomic scans and experimental breeding, two distinct multi-locus phenomena can mimic single-locus heterozygote advantage without any single locus being intrinsically overdominant. Associative overdominance (or pseudo-overdominance), analyzed by Tomoko Ohta (1971), occurs when two tightly linked loci in repulsion-phase genetic linkage disequilibrium carry completely or partially recessive deleterious alleles (A₁b₂ on chromosome 1 and a₂B₁ on chromosome 2). Here, both chromosomal homozygotes (A₁b₂/A₁b₂ and a₂B₁/a₂B₁) express a recessive genetic defect, whereas the repulsion double heterozygote (A₁b₂/a₂B₁) complements both defects (w₁₂ = 1 > w₁₁, w₂₂). Pseudo-overdominance is common inside chromosomal inversions and centromeric low-recombination regions, where genetic recombination cannot easily generate the superior coupling-phase chromosome (A₁B₁).

Alan Robertson's Paradox in Finite Populations (N_e)

While deterministic equations predict that any overdominant locus (s, t > 0) maintains polymorphism indefinitely, Alan Robertson (1962) proved using diffusion theory that in a finite population of effective population size N_e, asymmetric overdominance can actually accelerate allele fixation compared to neutrality. If the equilibrium frequency q̂ = s / (s + t) lies outside the interval (0.20, 0.80) (for instance, s = 0.01, t = 0.20 ⇒ q̂ ≈ 0.048) and the composite selection parameter 2N_e(s + t) is not sufficiently large, random genetic drift pushes the rare allele across the boundary q = 0 faster than weak selection against the common homozygote can restore it. Thus, long-term retention of overdominant polymorphism in finite populations requires both strong selection (N_es ≫ 1, N_et ≫ 1) and relatively symmetric selection coefficients (0.2 < p̂ < 0.8).

How can an allele that is lethal in homozygotes (t = 1) remain common in a population?

When an allele is a recessive lethal (t = 1, so w₂₂ = 0), its equilibrium frequency under heterozygote advantage simplifies to q̂ = s / (s + 1), where s is the fitness disadvantage of the wild-type homozygote relative to the heterozygote. Even a moderate 10% to 15% survival advantage for heterozygotes (s = 0.10–0.15)—such as protection against lethal childhood malaria in HbA/HbS carriers—stabilizes the lethal allele at q̂ = 0.10 / 1.10 ≈ 9.1% to q̂ = 0.15 / 1.15 ≈ 13.0%, meaning roughly 18% to 23% of the adult population will be heterozygous carriers (2p̂q̂).

Why doesn't natural selection eliminate the segregation load of overdominance?

Segregation load (L_seg = st / [s + t]) is an inescapable consequence of diploid Mendelian segregation at a single locus: two alleles residing on homologous chromosomes must separate during anaphase I of meiosis. However, over long evolutionary timescales, gene duplication via unequal crossing-over can resolve segregation load. If a duplicated chromosome arises carrying both functional alleles in tandem (A₁-A₂ on the same chromatid), every gamete transmits both variants simultaneously, creating fixed permanent heterozygosity across all offspring with zero segregation load—a mechanism documented in tandemly duplicated esterase resistance genes in Culex pipiens mosquitoes and vertebrate MHC gene families.

Does Cystic Fibrosis (CFTR ΔF508) persist in European populations due to heterozygote advantage?

The three-base-pair deletion ΔF508 in the CFTR chloride channel gene occurs at a carrier frequency of ~3–4% (q ≈ 0.02) in populations of European descent, far above simple mutation-selection balance (q̂ = √[μ/t] ≈ 0.001). Several heterozygote advantage hypotheses propose that CFTR/+ heterozygotes experienced reduced secretory fluid loss during historic epidemics of secretory diarrheas (such as cholera or enterotoxigenic Escherichia coli), reduced uptake of Salmonella enterica serovar Typhi (typhoid fever, which uses CFTR as an epithelial entry receptor; Pier et al., 1998), or increased resistance to tuberculosis, though historical demographic bottlenecks are also debated as contributing factors.

1. Allison, A. C. (1954). "Protection afforded by sickle-cell trait against subtertian malarial infection." British Medical Journal, 1(4857), 290–294. doi:10.1136/bmj.1.4857.290
2. Fisher, R. A. (1922). "On the dominance ratio." Proceedings of the Royal Society of Edinburgh, 42, 321–341. doi:10.1017/S0370164600023993
3. Robertson, A. (1962). "Selection for heterozygotes in small populations." Genetics, 47(9), 1291–1300. doi:10.1093/genetics/47.9.1291
4. Krieger, U., Lippman, Z. B., & Zamir, D. (2010). "The flowering gene SINGLE FLOWER TRUSS drives heterosis for yield in tomato." Nature Genetics, 42(5), 459–463. doi:10.1038/ng.550