Mutation-Selection Balance Calculator

Calculate equilibrium allele frequency under mutation-selection balance.

Typical: 10⁻⁵ - 10⁻⁶
s against homozygous aa
h=0: fully recessive, h=0.5: additive
Usually negligible
Please check your inputs and try again.

Why Recessive and Dominant Differ So Much

At equilibrium, new mutations exactly replace those removed by selection. The equilibrium frequency depends critically on whether selection can see the allele in heterozygotes.

ModeEquilibrium frequencyReason
Recessive (h = 0)q = √(μ/s)Hidden in heterozygotes
Dominant (h = 1)q = μ/sExposed every generation
Partial (0 < h < 1)q ≈ μ/(hs)Partly exposed

The square root makes an enormous difference. At μ = 10−5 and s = 1, a recessive allele sits at frequency 0.0032 while a dominant one sits at 10−5 — over three hundred times rarer. Recessive alleles hide from selection in heterozygotes, so far more copies persist.

Note how small a dominance coefficient is needed to change the picture. Even h = 0.02 makes the allele behave far more like a dominant, because selection acting on the many heterozygotes outweighs selection on the rare homozygotes.

Why Some Diseases Are Commoner Than Predicted

Cystic fibrosis carrier frequency in European populations is around 1 in 25, which is higher than mutation–selection balance alone predicts. Several explanations are proposed:

ExplanationMechanism
Heterozygote advantageProposed resistance to cholera or typhoid — debated
Founder effectsAncestral population bottleneck
Genetic driftSmall historical population sizes
Relaxed selectionModern medicine reduces s

The heterozygote advantage hypothesis remains unproven for CF, though it is well established for sickle cell and malaria. Where it does operate, the equilibrium is set by balancing selection rather than mutation–selection balance, and frequencies can be far higher.

Genetic Load

The mutational load is the fitness reduction a population carries from deleterious alleles at equilibrium. A striking result is that load depends on mutation rate, not on selection coefficient — a strongly deleterious allele is rarer but each occurrence costs more, and the two effects cancel exactly.

Worked Examples

Example 1: CF (recessive): μ=1e-4, s=1
q_eq=√(1e-4/1)=0.01
Result: ~1/100 frequency, 1/50 carriers
Actual CF carrier freq ~1/25 — other factors
Example 2: Dominant lethal: μ=1e-5, s=1, h=1
q_eq=1e-5/(1×1)=1e-5
Result: New mutations only, no inheritance
achondroplasia: almost all de novo mutations
Example 3: The square root effect
μ = 10−5, s = 1, recessive versus dominant
Result: q = 0.0032 versus 10−5
Over three hundred times more common when recessive, because heterozygotes shield the allele from selection entirely.
Example 4: Why partial dominance matters
h = 0 versus h = 0.02
Result: Frequency drops roughly sixfold
Selection on the many heterozygotes outweighs selection on rare homozygotes, so a tiny dominance coefficient shifts the equilibrium substantially.

Common Mistakes

⚠️
Using q = μ/s for recessive alleles

Recessives are hidden in heterozygotes, so equilibrium frequency is the square root of μ/s — far higher than the dominant expression gives.

⚠️
Assuming h = 0 for real recessives

Most so-called recessive alleles have small but non-zero dominance. Even h = 0.02 shifts behaviour substantially toward the dominant case.

⚠️
Attributing every elevated frequency to heterozygote advantage

Founder effects, drift and relaxed selection also raise frequencies. Heterozygote advantage is established for sickle cell but unproven for cystic fibrosis.

⚠️
Expecting load to fall with stronger selection

Mutational load depends on mutation rate, not selection coefficient. Stronger selection makes alleles rarer but each costs proportionally more.

Frequently Asked Questions

Equilibrium concept?
Mutation adds deleterious alleles each generation; selection removes them. At equilibrium: input rate = removal rate. Recessive: many carriers shield allele from selection → higher equilibrium freq. Dominant: selection efficient → very low freq.
Genetic load?
Total fitness reduction due to deleterious alleles. For recessive: load ≈ 2μ (independent of s for recessive). For dominant: load ≈ 2μ×s. Humans carry ~250 loss-of-function variants on average. Most are recessive and heterozygous.
Why are recessive deleterious alleles more common?
Because heterozygotes are unaffected, so selection cannot remove most copies. Equilibrium frequency is √(μ/s) rather than μ/s.
Why is cystic fibrosis more common than predicted?
Possibly heterozygote advantage, though this remains unproven for CF. Founder effects and genetic drift in ancestral populations are also plausible contributors.
What is genetic load?
The fitness reduction a population carries from deleterious alleles at equilibrium. It depends on mutation rate rather than selection coefficient.
How much dominance is needed to matter?
Very little. Even h = 0.02 shifts equilibrium substantially, because selection acts on the far more numerous heterozygotes.