Population Genetics - College Biology

SECTION 1

Historical Context & Motivation

In the early twentieth century, a profound intellectual crisis threatened to fracture the biological sciences. Charles Darwin's theory of evolution by natural selection, published in 1859, proposed that organisms with favorable traits survive and reproduce at higher rates, gradually transforming populations over time. Meanwhile, Gregor Mendel's rediscovered laws of inheritance, which described discrete hereditary factors (what we now call genes), seemed to many biologists irreconcilable with Darwin's vision of continuous, gradual change. How could evolution proceed smoothly if inheritance was particulate? This apparent conflict set the stage for population genetics, the discipline that would ultimately unify Mendelian genetics with Darwinian evolution.

1908 • Hardy–Weinberg Principle

Godfrey Harold Hardy and Wilhelm Weinberg independently demonstrated that allele frequencies in a population remain constant across generations in the absence of evolutionary forces, providing a mathematical null model for detecting evolution.

1918–1932 • The Modern Synthesis Founders

Ronald Fisher, J.B.S. Haldane, and Sewall Wright developed rigorous mathematical frameworks showing that Mendelian inheritance was fully compatible with natural selection operating on continuous traits.

1942 • The Modern Evolutionary Synthesis

Ernst Mayr and Theodosius Dobzhansky integrated population genetics with systematics and field biology, establishing the Modern Synthesis as the central paradigm of evolutionary biology.

1968 • Neutral Theory of Molecular Evolution

Motoo Kimura proposed that the majority of molecular-level evolutionary changes are selectively neutral and driven by genetic drift rather than natural selection, challenging adaptationist assumptions and expanding the theoretical scope of population genetics.

2000s–present • Genomics Era

High-throughput sequencing enables researchers to survey allele frequencies across entire genomes, applying population genetic models to detect selection signatures, reconstruct demographic histories, and inform conservation and medical genetics.

The central question that population genetics addresses is both deceptively simple and deeply consequential: what forces cause allele frequencies to change—or remain stable—within populations over time? By framing evolution as a change in allele frequencies rather than a change in individual organisms, population genetics provides the quantitative toolkit needed to predict evolutionary trajectories, test hypotheses about adaptation and drift, and bridge the gap between micro- and macroevolutionary processes.

SECTION 2

Core Principles & Definitions

Population genetics operates on a distinct level of biological organization. Rather than studying individual organisms or single genes in isolation, it examines the gene pool—the complete set of alleles present in a breeding population at a given time. An allele frequency (sometimes called gene frequency) is the proportion of a specific allele among all copies of that gene in the population, and tracking changes in these frequencies is the fundamental metric by which we quantify evolution at the population level. The discipline is built upon several foundational ideas that connect Mendelian genetics, probability theory, and evolutionary biology into a coherent analytical framework.

Allele & Genotype Frequencies

Allele frequency (p, q) represents the relative abundance of each allele at a locus. Genotype frequency describes the proportion of individuals carrying each genotype (AA, Aa, aa). These are connected by the Hardy–Weinberg equation: p² + 2pq + q² = 1.

Hardy–Weinberg Equilibrium

A population is at Hardy–Weinberg equilibrium when allele frequencies remain constant across generations. This requires no mutation, no migration, random mating, infinite population size, and no selection—conditions that serve as a null hypothesis for detecting evolutionary change.

Evolutionary Forces

Five major mechanisms disrupt equilibrium: natural selection (differential fitness), genetic drift (random sampling), mutation (new alleles), gene flow (migration), and nonrandom mating.

Fitness & Selection Coefficient

Fitness (w) quantifies an organism's relative reproductive success. The selection coefficient (s) measures the reduction in fitness of a genotype relative to the most fit genotype (w = 1 − s), providing the basis for predicting how rapidly selection alters allele frequencies.

Effective Population Size (Nₑ)

The effective population size is the number of individuals in an idealized population that would experience the same rate of genetic drift as the actual population. It is typically smaller than the census size due to unequal sex ratios, variance in reproductive success, and population bottlenecks.

✦

KEY TAKEAWAY

Think of a population's gene pool as a bag of colored marbles. Each marble represents one allele copy. The Hardy–Weinberg principle tells us that if you draw marbles randomly and put them back (no bias, no adding new colors, no marbles leaving or entering), the proportion of each color stays the same indefinitely. Evolution is any process that changes the marble ratios—whether that's selectively removing one color (selection), accidentally losing some when marbles spill (drift), or introducing new colors from another bag (gene flow).

SECTION 3

Visualizing Allele Frequency Change

The following diagram illustrates how the five major evolutionary forces act upon allele frequencies in a population's gene pool, starting from Hardy–Weinberg equilibrium at the center and showing how each mechanism drives the system away from that equilibrium state. Understanding these forces visually helps clarify why real populations almost never conform perfectly to Hardy–Weinberg expectations.

Five evolutionary forces radiate outward from Hardy–Weinberg equilibrium. Natural selection and genetic drift are the primary agents of allele frequency change, while mutation, gene flow, and nonrandom mating each contribute in distinct ways.

Notice that Hardy–Weinberg equilibrium occupies the center of the diagram: it is the default state that persists only when none of the five forces are acting. In reality, every natural population is subject to at least some of these forces simultaneously, which is why allele frequencies are nearly always changing—sometimes imperceptibly slowly, sometimes quite rapidly. The power of the Hardy–Weinberg model lies not in its literal accuracy as a description of nature, but in its utility as a benchmark against which deviations can be measured and attributed to specific evolutionary mechanisms.

SECTION 4

Mathematical Framework

The quantitative backbone of population genetics rests on a few elegant equations that connect allele frequencies to genotype frequencies and predict how evolutionary forces alter them over time. These equations derive from basic probability theory applied to diploid organisms with Mendelian inheritance, and they assume a single autosomal locus with two alleles as the simplest illustrative case.

Hardy–Weinberg Equations

ALLELE FREQUENCY CONSTRAINT

p + q = 1

Where p = frequency of the dominant allele (A), and q = frequency of the recessive allele (a). For a biallelic locus, these two frequencies must sum to 1.

HARDY–WEINBERG GENOTYPE FREQUENCIES

p² + 2pq + q² = 1

p² = frequency of homozygous dominant (AA); 2pq = frequency of heterozygotes (Aa); q² = frequency of homozygous recessive (aa). This equation is simply the binomial expansion of (p + q)², reflecting the random union of gametes.

Selection Model

CHANGE IN ALLELE FREQUENCY UNDER SELECTION

Δq = −spq² / w̄

This simplified form (for selection against a recessive homozygote) describes the change in frequency of allele q per generation. Here, s is the selection coefficient (fitness cost to aa), p = 1 − q is the frequency of the dominant allele, and w̄ is the mean fitness of the population (w̄ = p² + 2pq + q²(1 − s)). When s > 0, the recessive allele declines in frequency each generation, but increasingly slowly as q becomes small because the allele hides in heterozygotes.

Genetic Drift

VARIANCE IN ALLELE FREQUENCY DUE TO DRIFT

Var(Δq) = pq / (2Nₑ)

The variance in allele frequency change per generation is inversely proportional to the effective population size (Nₑ). In small populations, Var(Δq) is large, meaning substantial random fluctuations can occur; in large populations, drift becomes negligible and selection dominates. The ratio Nₑs determines whether selection (Nₑs >> 1) or drift (Nₑs << 1) governs the fate of an allele.

ℹ

When Does Drift Override Selection?

A key threshold in population genetics is Nₑs ≈ 1. When Nₑs >> 1, natural selection efficiently drives allele frequency changes and drift is a minor perturbation. When Nₑs << 1, the selective advantage or disadvantage of an allele is too small relative to the noise of random sampling, and the allele behaves as if it were selectively neutral. This is why Kimura's neutral theory predicts that many molecular polymorphisms are maintained by drift alone.

SECTION 5

Genetic Drift: Population Size & Stochasticity

One of the most important insights from population genetics is that evolution is not purely deterministic. Genetic drift—the random fluctuation of allele frequencies due to finite sampling of gametes each generation—can cause alleles to be lost or fixed entirely by chance, independent of their effects on fitness. The magnitude of drift is inversely proportional to population size: small populations experience dramatic stochastic swings, while large populations are buffered against random change. The diagram below contrasts the trajectories of a neutral allele (starting at q = 0.5) in populations of different sizes, illustrating how drift leads to fixation or loss far more rapidly in smaller populations.

Two replicate simulations (solid and dashed lines) are shown for each population size. In the N = 20 population (red), alleles rapidly drift to fixation (q = 1.0) or loss (q = 0.0). The N = 200 population (amber) shows moderate fluctuations. The N = 10,000 population (cyan) remains nearly stable at q = 0.5, demonstrating the buffering effect of large population size.

Two special cases of extreme drift deserve attention. A bottleneck effect occurs when a population's size is drastically reduced by a catastrophic event—such as a natural disaster, disease outbreak, or habitat destruction—causing a random subset of alleles to survive. The resulting population may have dramatically different allele frequencies than the original, regardless of which alleles were adaptive. Similarly, the founder effect arises when a small group of individuals colonizes a new habitat, carrying only a fraction of the original gene pool. Both phenomena can lead to the fixation of otherwise rare alleles and the loss of genetic diversity, with lasting consequences for the population's evolutionary potential and susceptibility to inbreeding depression.

Feature	Bottleneck Effect	Founder Effect
Cause	Population crash (disaster, disease, habitat loss)	Colonization of new area by small group
Mechanism	Random survival of a subset of original population	Sampling of alleles carried by founders
Genetic outcome	Reduced diversity; shifted allele frequencies	Reduced diversity; some alleles over- or underrepresented
Example	Northern elephant seals (hunted to ~20 individuals in 1890s)	Amish populations with elevated frequency of Ellis–van Creveld syndrome

Comparison of two important drift-amplifying events

SECTION 6

Worked Example: Hardy–Weinberg Analysis

Suppose you are studying a population of wildflowers in which petal color is controlled by a single autosomal locus with two alleles. The C^R allele produces red petals and is completely dominant over the C^W allele, which produces white petals. You survey 500 plants and find that 80 have white flowers. Determine the allele and genotype frequencies, and predict the number of heterozygous carriers.

Hardy–Weinberg Allele and Genotype Frequency Calculation

1
Step 1 — Identify the recessive phenotype frequency

White-flowered plants are homozygous recessive (C^WC^W). Their frequency is 80/500 = 0.16. Under Hardy–Weinberg assumptions, this frequency equals q².

q² = 0.16
2
Step 2 — Calculate q (frequency of $C^W$ allele)

Taking the square root of q²: q = √0.16 = 0.4. This means the C^W allele constitutes 40% of all alleles at this locus.

q = 0.4
3
Step 3 — Calculate p (frequency of $C^R$ allele)

Since p + q = 1, we find p = 1 − 0.4 = 0.6. The dominant allele frequency is 60%.

p = 0.6
4
Step 4 — Determine all genotype frequencies

Applying the Hardy–Weinberg equation: homozygous dominant (C^RC^R) = p² = (0.6)² = 0.36; heterozygous (C^RC^W) = 2pq = 2(0.6)(0.4) = 0.48; homozygous recessive (C^WC^W) = q² = 0.16. Verification: 0.36 + 0.48 + 0.16 = 1.00 ✓

p² = 0.36, 2pq = 0.48, q² = 0.16
5
Step 5 — Predict the number of heterozygous carriers

The heterozygote frequency is 0.48, so in a population of 500: 0.48 × 500 = 240 plants are expected to be heterozygous carriers. These plants appear phenotypically red (because C^R is dominant) but carry one copy of the C^W allele. Notably, this is the largest genotype class—the 'hidden' genetic variation that selection cannot directly see when the allele is recessive.

240 heterozygous carriers out of 500 plants

SECTION 7

Comparing Evolutionary Forces

Each of the five evolutionary forces has distinct properties regarding directionality, predictability, and its effect on genetic variation within and between populations. Understanding these differences is essential for correctly interpreting population genetic data—for example, distinguishing whether an unusual allele frequency pattern reflects strong directional selection or a recent bottleneck event. The table below synthesizes the key properties of each force.

Evolutionary Force	Directional?	Effect on Variation	Depends on Pop. Size?
Natural Selection	Yes — favors specific alleles based on fitness	Can increase (balancing) or decrease (directional, disruptive) variation	Partially — efficacy depends on Nₑs
Genetic Drift	No — random fluctuations, no preferred direction	Decreases within-population variation; increases between-population divergence	Yes — inversely proportional to Nₑ
Mutation	Weakly — introduces new alleles at very low rates	Increases variation (the ultimate source of all genetic novelty)	No — mutation rate is per-gene, per-generation
Gene Flow	Yes — moves alleles from source to recipient population	Increases within-population variation; decreases between-population divergence	Partially — affected by migration rate (m)
Nonrandom Mating	No — alters genotype frequencies, not allele frequencies directly	Inbreeding increases homozygosity; assortative mating shifts genotype ratios	Indirectly — inbreeding effects amplified in small populations

Properties of the five major evolutionary forces in population genetics

✦

KEY TAKEAWAY

Think of population genetics as a weather forecasting system for evolution. Selection is like prevailing wind patterns—it has a consistent direction determined by the environment. Drift is like turbulence—random perturbations that matter most when the system is small (a leaf in a wind tunnel versus a leaf in the open sky). Gene flow is like weather fronts moving between regions, homogenizing conditions across areas. Understanding which 'weather system' dominates in a given population is the core skill of a population geneticist.

SECTION 8

Connections to Advanced Evolutionary Theory

The principles of population genetics introduced here form the foundation for several advanced areas of modern evolutionary biology. As you move beyond introductory models, the mathematics becomes richer and the biological complexity deepens considerably. The Hardy–Weinberg framework extends naturally into multi-locus models, coalescent theory, and quantitative genetics—each of which addresses limitations of the single-locus, two-allele models we have explored. Understanding where introductory population genetics ends and these advanced frameworks begin helps you appreciate both the power and the boundaries of the tools developed in this lesson.

Introductory Concept	Advanced Extension	Key Addition
Hardy–Weinberg equilibrium (single locus, 2 alleles)	Multi-locus models & linkage disequilibrium	Non-independent assortment of alleles at linked loci; recombination rates affect how selection at one locus influences nearby loci
Genetic drift as variance in allele frequency	Coalescent theory	Models the genealogy of alleles backward in time; allows inference of population history from DNA sequence data
Selection on single genes	Quantitative genetics	Extends to polygenic traits influenced by many loci; uses heritability and the breeder's equation (R = h²S)
Gene flow between two populations	Landscape genetics & F-statistics	Uses Wright's F_ST to quantify population structure and gene flow across complex spatial landscapes
Neutral versus selected alleles	Molecular evolution & phylogenomics	Tests for selection using dN/dS ratios, McDonald–Kreitman tests, and genome-wide scans for selective sweeps

From introductory to advanced population genetics

Modern population genetics is increasingly intertwined with genomics and bioinformatics. Genome-wide association studies (GWAS) use population genetic principles to map disease-associated variants in humans, while conservation geneticists apply effective population size estimates and heterozygosity metrics to assess the viability of endangered species. As you encounter these applications in upper-division courses, you will find that the Hardy–Weinberg equation, drift models, and selection coefficients introduced here are not merely textbook abstractions—they are the working tools of researchers studying evolution in real time.

SECTION 9

Practice Problems

Question 1: CONCEPTUAL

A population geneticist observes that allele frequencies at a particular locus have not changed over 50 generations despite the population experiencing significant environmental fluctuations. However, the population is known to be large (N > 100,000), no migration occurs, mating is random, and mutation rates are negligible. Explain why allele frequencies might remain stable even though the environment changed.

Show answer

If the environmental fluctuations did not alter the relative fitness of genotypes at this locus—for example, if the locus controls a trait unrelated to the changing environmental variable—then natural selection would not act differentially on the alleles. With a large population size (minimizing drift), no migration, random mating, and negligible mutation, all five Hardy–Weinberg conditions are approximately met, and allele frequencies remain stable. This illustrates the important principle that not every locus is under selection in every environment; many loci may be selectively neutral with respect to a given ecological change, consistent with Kimura's neutral theory.

Question 2: BASIC CALCULATION

In a population of 1,000 individuals, 90 express the recessive phenotype for a single-locus, two-allele trait. Assuming Hardy–Weinberg equilibrium, calculate the frequencies of both alleles (p and q) and the expected number of heterozygous individuals in the population.

Show answer

The frequency of the recessive phenotype equals q² = 90/1000 = 0.09. Taking the square root, q = √0.09 = 0.3, which is the frequency of the recessive allele. Since p + q = 1, p = 1 − 0.3 = 0.7. The heterozygote frequency is 2pq = 2(0.7)(0.3) = 0.42, so the expected number of heterozygous individuals is 0.42 × 1,000 = 420. Note that heterozygotes constitute the largest genotype class, which is always the case when allele frequencies are between approximately 0.3 and 0.7.

Question 3: INTERMEDIATE

A recessive lethal allele (a) has a current frequency of q = 0.02 in a large population. If homozygous recessive individuals (aa) die before reproduction (s = 1 against aa), but heterozygotes (Aa) and homozygous dominants (AA) have equal fitness, calculate: (a) the mean fitness of the population (w̄), and (b) the change in allele frequency (Δq) in one generation. Explain why complete lethality of aa does not immediately eliminate the a allele.

Show answer

With p = 0.98 and q = 0.02, the mean fitness is w̄ = p² + 2pq + q²(1 − s) = (0.98)² + 2(0.98)(0.02) + (0.02)²(0) = 0.9604 + 0.0392 + 0 = 0.9996. The change in allele frequency Δq = −spq² / w̄ = −(1)(0.98)(0.0004)/0.9996 ≈ −0.000392. This tiny change (from q = 0.02 to approximately 0.01961) demonstrates that selection against a recessive allele becomes extremely inefficient at low frequencies because the vast majority of recessive alleles are sheltered in heterozygotes (2pq = 0.0392, or about 3.92% of the population) where they are phenotypically invisible to selection. Only q² = 0.04% of individuals express the lethal phenotype.

Question 4: APPLIED

Cystic fibrosis (CF) is an autosomal recessive disorder caused by mutations in the CFTR gene. Among individuals of Northern European descent, approximately 1 in 2,500 newborns is affected. Using Hardy–Weinberg analysis: (a) estimate the carrier frequency, (b) calculate the probability that two randomly chosen unaffected individuals are both carriers, and (c) discuss one biological reason the CF allele might be maintained at a frequency higher than expected from mutation-selection balance alone.

Show answer

(a) The frequency of affected individuals q² = 1/2500 = 0.0004, so q = √0.0004 = 0.02 and p = 0.98. The carrier (heterozygote) frequency is 2pq = 2(0.98)(0.02) = 0.0392, or approximately 1 in 25 individuals. (b) The probability that a randomly chosen unaffected person is a carrier equals 2pq/(1 − q²) ≈ 0.0392/0.9996 ≈ 0.0392. The probability that two randomly chosen unaffected individuals are both carriers is (0.0392)² ≈ 0.00154, or about 1 in 650 couples. (c) The heterozygote advantage hypothesis proposes that CF carriers have increased resistance to certain gastrointestinal infections (such as cholera or typhoid fever), giving them a selective advantage in environments where these diseases are prevalent. This form of balancing selection—analogous to sickle cell trait and malaria—could maintain the CF allele at a higher frequency than predicted by mutation-selection balance alone.

Question 5: CRITICAL THINKING

You are studying two island populations of a lizard species. Island A has Nₑ = 50 and Island B has Nₑ = 5,000. Both populations start with the same allele frequency (q = 0.5) for a neutral locus. After 100 generations, you genotype both populations. (a) Which population is more likely to have lost the allele entirely (q = 0 or q = 1)? Justify quantitatively using the drift variance formula. (b) If a mildly beneficial mutation (s = 0.001) arises on each island, on which island is selection more likely to fix the allele? Use the Nₑs criterion in your argument. (c) Synthesize your answers to explain the 'paradox of small populations' for conservation biology.

Show answer

(a) Island A (Nₑ = 50) is far more likely to have lost the allele. The variance in allele frequency change per generation due to drift is Var(Δq) = pq/(2Nₑ). For Island A: Var = (0.5)(0.5)/(2 × 50) = 0.0025 per generation. For Island B: Var = (0.5)(0.5)/(2 × 5000) = 0.000025. The standard deviation of change is 10× larger for Island A (0.05 vs. 0.005), meaning after 100 generations, Island A has had ample stochastic opportunity to drift to fixation or loss. (b) For Island A, Nₑs = 50 × 0.001 = 0.05, which is much less than 1, meaning the mutation behaves as effectively neutral—its selective advantage is overwhelmed by drift. For Island B, Nₑs = 5000 × 0.001 = 5, well above 1, so selection can efficiently promote fixation. (c) The 'paradox of small populations' is that they simultaneously suffer more genetic drift (eroding diversity and fixing deleterious alleles) and are less able to respond to beneficial mutations via selection. Small populations thus experience a double jeopardy: they lose existing adaptive variation through drift while being unable to incorporate new adaptive variation through selection. This is why maintaining effective population sizes above critical thresholds is a central priority in conservation genetics.

SUMMARY

Summary

Population genetics is the quantitative study of how allele frequencies change within populations over time, providing the mathematical foundation for evolutionary biology. The Hardy–Weinberg equilibrium (p² + 2pq + q² = 1) serves as the null model: allele frequencies remain constant across generations when no evolutionary forces act. Five mechanisms disrupt this equilibrium—natural selection (differential fitness alters allele frequencies directionally), genetic drift (random sampling causes stochastic fluctuations, especially potent in small populations), mutation (introduces new alleles as the raw material of evolution), gene flow (migration homogenizes allele frequencies between populations), and nonrandom mating (alters genotype proportions without directly changing allele frequencies).

Key quantitative tools include the selection coefficient (s) and fitness (w) for modeling selection, the drift variance formula Var(Δq) = pq/(2Nₑ) for quantifying stochastic effects, and the critical threshold Nₑs ≈ 1 that determines whether selection or drift governs allele fate. Special drift events—bottlenecks and founder effects—can dramatically reshape genetic diversity. These foundational concepts connect directly to advanced topics including coalescent theory, quantitative genetics, landscape genetics, and genomic analyses of natural selection, making population genetics an indispensable toolkit for understanding how evolution operates at every scale.