Evolutionary Biology Laboratory Manual: Population Genetics, Molecular Clocks & Morphometrics
Experimental protocols in Hardy-Weinberg equilibrium verification and Chi-square (χ²) goodness-of-fit testing, microevolutionary simulations of natural selection against recessive lethals and overdominant balancing selection, stochastic Wright-Fisher genetic drift and effective population size ($N_e$) estimation, molecular clock calibration and $dN/dS$ (ω) selection profiling, bivariate allometric cranial morphometrics, and Dobzhansky-Muller speciation incompatibility modeling aligned with CSIR-NET Unit 11 and GATE EY syllabi.
Population Sampling Ethics, Fossil Handling & Simulation Rigor Protocols
Evolutionary biology bench protocols combine live model organisms (e.g., Drosophila melanogaster), comparative osteological/fossil casts, and Monte Carlo computational simulations:
- Live Organism Biosafety (Drosophila): When scoring phenotypic markers (white eye, vestigial wing, ebony body), anesthetize flies using FlyNap (triethylamine) or carbon dioxide ($CO_2$) pads. Never expose flies to excessive anesthetic that causes mortality. Autoclave and freeze morgue bottles prior to disposal.
- Osteological & Fossil Specimen Handling: Handle delicate vertebrate cranial casts and fossil hominin reproductions on cushioned felt mats. Use digital vernier calipers with smooth nylon/Teflon jaws to prevent scratching cranial bone landmarks (bregma, lambda, nasion, prosthion).
- Human Population Genetic Ethics: When collecting human phenotypic data (PTC tasting sensitivity, ABO blood grouping, Darwin’s ear tubercle), ensure informed voluntary consent, maintain strict anonymity of personal health data, and autoclave all lancets and alcohol swabs immediately in biohazard sharps containers.
- In Silico Simulation Verification: In Monte Carlo Wright-Fisher simulations, run a minimum of 1,000 replicate trials per parameter combination to establish robust probability distributions of allele fixation and heterozygosity decay.
Testing Hardy-Weinberg Equilibrium & Chi-Square ($chi^2$) Goodness-of-Fit in Natural and Model Populations
The Hardy-Weinberg Principle forms the null model of evolutionary biology. In an infinitely large, panmictic (randomly mating) diploid population free of mutation, migration, natural selection, and genetic drift, allele and genotype frequencies remain constant across generations:
- Binomial Expansion: For a locus with two alleles $A$ (frequency $p$) and $a$ (frequency $q$), where $p + q = 1$:
$$(p + q)^2 = p^2 + 2pq + q^2 = 1$$
where $p^2 = text{Freq}(AA)$, $2pq = text{Freq}(Aa)$, and $q^2 = text{Freq}(aa)$. - Allele Frequency Calculation (Direct Gene Counting):
$$p = frac{2 cdot N_{AA} + N_{Aa}}{2N}, quad q = frac{2 cdot N_{aa} + N_{Aa}}{2N}$$ - Pearson’s Chi-Square ($chi^2$) Goodness-of-Fit:
$$chi^2 = sum frac{(O – E)^2}{E} = frac{(N_{AA} – N p^2)^2}{N p^2} + frac{(N_{Aa} – 2Npq)^2}{2Npq} + frac{(N_{aa} – N q^2)^2}{N q^2}$$
The degrees of freedom is: $text{df} = k – 1 – m = 3 – 1 – 1 = 1$, where $k=3$ genotypic classes and $m=1$ estimated independent parameter ($p$). - Decision Rule: The critical value at $alpha = 0.05$ with $text{df} = 1$ is $chi^2_{text{critical}} = 3.841$. If $chi^2 > 3.841$, the null hypothesis of equilibrium is rejected, demonstrating the action of evolutionary forces or non-random mating.
- Phenotyping Assay: Phenylthiocarbamide (PTC) taste test paper strips (impregnated with $5text{ mu g}$ PTC) and control untreated paper strips; sanitized disposable forceps.
- Model Alternative: Two-color bead simulation kits (400 red beads representing allele $A$, 600 white beads representing allele $a$, total $N = 500$ diploid individuals).
- Human Phenotyping (PTC Tasting): Place a control paper strip on the tongue, followed by a PTC paper strip. Score individuals as Tasters (dominant phenotype, $TT$ or $Tt$, experiencing intensely bitter taste) or Non-Tasters (recessive phenotype, $tt$, tasting only cellulose). Record total sample size ($N$).
- Model Sampling (Direct Counting): In a model population with codominance (or bead drawing), draw 100 pairs of beads randomly without replacement to form 100 diploid genotypes: $AA$, $Aa$, and $aa$.
- Frequency Derivation: Calculate observed counts $N_{AA}$, $N_{Aa}$, and $N_{aa}$. Compute allele frequencies $p$ and $q$ by direct gene counting.
- Expected Genotypic Numbers: Compute expected numbers: $E(AA) = N cdot p^2$, $E(Aa) = N cdot 2pq$, and $E(aa) = N cdot q^2$.
- Chi-Square Testing: Calculate $chi^2 = sum frac{(O – E)^2}{E}$. Compare with $chi^2_{0.05, 1} = 3.841$ and state biological conclusion.
| Genotype Class | Observed ($O$) | Expected Proportion | Expected ($E$) | $(O – E)$ | $(O – E)^2 / E$ |
|---|---|---|---|---|---|
| $AA$ (Homozygous Dominant) | 140 | $p^2 = (0.60)^2 = 0.360$ | 144.0 | -4.0 | 0.111 |
| $Aa$ (Heterozygote) | 200 | $2pq = 2(0.6)(0.4) = 0.480$ | 192.0 | +8.0 | 0.333 |
| $aa$ (Homozygous Recessive) | 60 | $q^2 = (0.40)^2 = 0.160$ | 64.0 | -4.0 | 0.250 |
| Total ($N = 400$) | 400 | 1.000 | 400.0 | 0.0 | $chi^2 = 0.694$ |
Statistical Inference: Calculated $chi^2 = 0.694 < chi^2_{0.05, 1} = 3.841$ ($p = 0.405$). We fail to reject the null hypothesis; the population is in Hardy-Weinberg equilibrium at this locus.
Simulating Microevolutionary Forces: Natural Selection against Recessive Lethals vs Heterozygote Advantage (Sickle-Cell Model)
Natural selection acts as a deterministic evolutionary force altering allele frequencies based on differential reproductive success (relative fitness, $w$):
- Selection Against Recessive Lethal ($w_{AA} = 1, w_{Aa} = 1, w_{aa} = 0$): All homozygous recessive zygotes die prior to reproductive maturity ($s = 1.0$). The recessive allele frequency $q_t$ after $t$ generations decays as:
$$q_t = frac{q_0}{1 + t cdot q_0}$$
The number of generations required to reduce the frequency from $q_0$ to $q_t$ is: $t = frac{1}{q_t} – frac{1}{q_0}$. As $q$ becomes very small, almost all recessive alleles reside in heterozygous carriers ($Aa$), making complete elimination by natural selection virtually impossible. - Overdominant Balancing Selection (Heterozygote Advantage): When the heterozygote exhibits higher fitness than either homozygote ($w_{Aa} > w_{AA}$ and $w_{Aa} > w_{aa}$):
• $w_{AA} = 1 – s_1$, $w_{Aa} = 1$, $w_{aa} = 1 – s_2$
The population converges onto a stable internal equilibrium allele frequency ($hat{q}$):
$$hat{q} = frac{s_1}{s_1 + s_2}, quad hat{p} = frac{s_2}{s_1 + s_2}$$
In Sub-Saharan Africa, where $HbA/HbA$ individuals suffer mortality from Plasmodium falciparum malaria ($s_1 approx 0.12$) and $HbS/HbS$ individuals suffer severe sickle-cell anemia ($s_2 approx 0.85$), balancing selection actively maintains the lethal $HbS$ allele at $hat{q} approx 0.124$ (12.4%).
- Benchtop Simulation Kit: Large sorting containers with 500 red beads (allele $A$) and 500 green beads (allele $a$). Genotype collection cups labeled $AA$, $Aa$, and $aa$ (lethal discard bin).
- Digital Simulator: Populus 5.5 / R script running deterministic recurrence equations.
- Generation 0 Initialization: Mix 500 red ($A$) and 500 green ($a$) beads. $p_0 = q_0 = 0.50$. Draw 100 random pairs to form 100 individuals.
- Recessive Lethal Selection Step: Completely discard all $aa$ pairs (100% mortality, $s = 1$). Retain all $AA$ and $Aa$ pairs.
- Gene Pool Replenishment: Count surviving alleles in $AA$ and $Aa$ survivors. Calculate new $q_1$. Replenish the bead container to maintain total pool proportion. Repeat for generations $t = 1, 2, 3, 4, 5$.
- Overdominance Simulation: Reset pool. Impose selective penalties: discard a fraction $s_1 = 0.20$ of $AA$ pairs and $s_2 = 0.80$ of $aa$ pairs. Keep 100% of $Aa$ heterozygotes. Recompute frequencies across 10 cycles until $hat{q}$ stabilizes.
| Generation ($t$) | Recessive Lethal $q_t$ (Theoretical) | Recessive Lethal $q_t$ (Simulated) | Balancing Selection $q_t$ ($s_1=0.2, s_2=0.8$) |
|---|---|---|---|
| 0 | 0.500 | 0.500 | 0.500 |
| 1 | 0.333 | 0.328 | 0.320 |
| 2 | 0.250 | 0.245 | 0.245 |
| 5 | 0.143 | 0.138 | 0.202 (Equilibrium → 0.200) |
| 10 | 0.083 | 0.080 | 0.200 (Stable $hat{q}$) |
Evolutionary Conclusion: Purging efficiency plummets as $q$ decreases because $99%$ of lethal alleles reside in phenotypically normal heterozygotes. In contrast, overdominance drives the population inexorably to the stable attractor point $hat{q} = frac{0.2}{0.2 + 0.8} = 0.200$, maintaining genetic diversity permanently.
Genetic Drift, Founder Effects & Effective Population Size ($N_e$) in Finite Demes: Wright-Fisher Simulation
In finite populations, chance sampling of gametes produces random genetic drift, leading to non-adaptive evolutionary change independent of natural selection:
- Probability of Fixation & Fixation Time: The probability that a neutral allele will eventually become fixed equals its current frequency ($P(text{fixation}) = p_0$). In a diploid population of size $N$, the mean time to fixation for a newly arisen neutral mutation is:
$$bar{t}_{text{fixation}} = 4N text{ generations}$$
whereas a lost allele takes on average $bar{t}_{text{loss}} = 2 left(frac{N}{p}right) ln(2N)$ generations. - Decay of Heterozygosity: Drift continuously erodes genetic variation at a rate inversely proportional to population size:
$$H_t = H_0 left(1 – frac{1}{2N_e}right)^t$$ - Wright’s Effective Population Size ($N_e$): The size of an idealized Fisher-Wright population that would experience the same rate of genetic drift or inbreeding as the actual population:
• Unequal Sex Ratio: $N_e = frac{4 N_m N_f}{N_m + N_f}$ (A population of 2 males and 98 females has an $N_e approx 7.84$, not 100).
• Demographic Bottlenecks: Over $t$ generations with fluctuating sizes $N_1, N_2, dots, N_t$, $N_e$ is governed by the harmonic mean:
$$frac{1}{N_e} = frac{1}{t} sum_{i=1}^t frac{1}{N_i}$$
A single bottleneck generation severely depresses long-term genetic variation.
- Physical Deme Modeling: Multi-well culture plates with 20 beads per well (10 diploid individuals) representing isolated island populations (founder populations).
- Computational Tool: R / Python script running 10,000 Wright-Fisher binomial draws:
rbinom(1, 2*N, p) / (2*N).
- Small Deme Simulation ($N = 10$): Initialize 10 replicate populations with $N = 10$ individuals ($2N = 20$ alleles) at $p_0 = 0.50$. In each generation, draw 20 alleles with replacement according to current $p$. Track $p$ until every replicate either fixes ($p = 1.0$) or is lost ($p = 0.0$).
- Large Deme Comparison ($N = 500$): Repeat with $2N = 1000$ alleles over 50 generations. Observe the dampening of stochastic fluctuations.
- Founder Effect Simulation: From a large mainland pool ($p = 0.50$), draw 2 individuals ($N = 2$, 4 alleles) to colonize an empty island. Repeat 20 times. Record the dramatic disparity in initial island allele frequencies ($p = 0.0, 0.25, 0.50, 0.75, 1.0$).
- Effective Size Calculations: Compute $N_e$ for an endangered wildlife herd with 5 breeding males and 45 breeding females. Calculate $N_e$ for a population passing through sizes $1000 rightarrow 20 rightarrow 1000$ over 3 generations.
Demographic Bottleneck Impact:
Generations: $N_1 = 1000$, $N_2 = 25$ (Catastrophic bottleneck), $N_3 = 1000$.
$$frac{1}{N_e} = frac{1}{3} left(frac{1}{1000} + frac{1}{25} + frac{1}{1000}right) = frac{1}{3} (0.001 + 0.040 + 0.001) = frac{0.042}{3} = 0.014$$
$$N_e = frac{1}{0.014} approx mathbf{71.4text{ individuals}}$$
Conservation Implication: Despite a census recovery back to 1000 animals in generation 3, the effective genetic population size is permanently reduced to ~71, causing enduring loss of heterozygosity and founder effect inbreeding depression.
Molecular Clock & Divergence Dating: Estimating Speciation Times & $dN/dS$ ($omega$) Selection Ratios
Zuckerkandl and Pauling (1962) observed that amino acid replacements accumulate at an approximately constant rate over geological time, creating a molecular clock:
- Linear Molecular Clock: The genetic distance ($d$) between two lineages diverged for time $T$ is:
$$d = 2 cdot mu cdot T implies T = frac{d}{2mu}$$
where $mu$ is the evolutionary substitution rate per site per million years, and the factor of 2 accounts for independent divergence down both sister lineages from the common ancestor. - Neutral Theory Explanation (Kimura): For neutral mutations, the rate of molecular substitution ($k$) exactly equals the neutral mutation rate ($mu_{text{neutral}}$), completely independent of population size: $k = 2N mu cdot left(frac{1}{2N}right) = mu$.
- Codon Selection Ratio ($omega = dN / dS$):
• $dS$ (Synonymous substitutions per site): Silent mutations that do not alter amino acids; serves as a proxy for the baseline neutral mutation rate.
• $dN$ (Non-synonymous substitutions per site): Missense mutations altering the amino acid.
• $omega = frac{dN}{dS} < 1.0$: Purifying (negative) selection eliminating deleterious changes.
• $omega approx 1.0$: Neutral evolution.
• $omega > 1.0$: Positive (Darwinian diversifying) selection actively favoring adaptive innovations.
- Software: MEGA 11 (Timetree & Codon Selection Analysis), PAML 4.9 (codeml package), and HyPhy suite.
- Sequence Dataset: Primate mitochondrial Cytochrome c Oxidase Subunit I (COI) coding sequences and primate fossil calibrations (Human-Chimp split: 6-8 MYA, Catarrhini-Platyrrhini split: 35 MYA).
- Codon Alignment: Align coding nucleotide sequences of COI using translation-aware aligner (MUSCLE in MEGA). Verify no frameshifting stop codons occur.
- Evolutionary Distance: Compute evolutionary distance $d$ using the Tamura-Nei model with gamma-distributed rate variation.
- Clock Calibration: Using the 35 MYA fossil divergence point between Old World and New World monkeys ($d = 0.140$), calibrate $mu$:
$mu = frac{d}{2T} = frac{0.140}{2 times 35} = 0.0020text{ substitutions/site/MYA}$. - Node Dating: Calculate divergence time for the Human-Chimpanzee split ($d = 0.026$):
$T = frac{0.026}{2 times 0.0020} = 6.5text{ MYA}$. - Codon-by-Codon $dN/dS$ Analysis: Run
Selection → Codon-based Test of Positive Selectionin MEGA using the Nei-Gojobori method. Identify specific codon sites under positive selection ($omega > 1, p < 0.05$).
| Lineage Pair | Pairwise Distance ($d$) | Estimated Divergence Time ($T$) | Fossil Benchmark |
|---|---|---|---|
| Human – Chimpanzee | 0.026 | 6.5 MYA | 6.0 – 7.5 MYA (Sahelanthropus) |
| Human – Gorilla | 0.038 | 9.5 MYA | 8.5 – 10.5 MYA (Chororapithecus) |
| Human – Orangutan | 0.062 | 15.5 MYA | 14.0 – 16.0 MYA (Sivapithecus) |
Selective Constraint Assessment: Mean $dN = 0.004$, mean $dS = 0.052 implies omega = frac{0.004}{0.052} = 0.077 ll 1.0$. The mitochondrial COI gene is subjected to intense purifying selection ($>92%$ of amino acid alterations eliminated) to maintain critical electron transport chain bioenergetic efficiency.
Evolutionary Morphometrics & Allometric Growth: Bivariate Allometry ($y = b cdot x^a$) & Cranial Indexes
Morphological evolution frequently proceeds through heterochrony and differential growth rates among body components:
- Huxley’s Allometric Power Law (1932): Describes the non-linear relationship between trait size ($y$) and body/skull size ($x$):
$$y = b cdot x^a iff log_{10}(y) = log_{10}(b) + a cdot log_{10}(x)$$
where $a$ is the allometric growth ratio (scaling exponent) and $b$ is the initial growth index. - Scaling Classifications:
• Isometry ($a = 1.0$): Trait grows in direct geometric proportion to overall size.
• Positive Allometry / Hyperallometry ($a > 1.0$): Trait becomes disproportionately larger as body size increases (e.g., weapon traits, Irish elk antlers, sexual display ornaments).
• Negative Allometry / Hypoallometry ($a < 1.0$): Trait becomes disproportionately smaller with increasing body size (e.g., mammalian brain-to-body mass scaling, $a approx 0.67 – 0.75$). - Retzius Cephalic / Cranial Index:
$$text{Cranial Index} = left(frac{text{Maximum Cranial Breadth (Euryon – Euryon)}}{text{Maximum Cranial Length (Glabella – Opisthocranion)}}right) times 100$$
• Dolichocephalic ($< 75$): Long, narrow skull (e.g., early hominins, Homo erectus).
• Mesocephalic ($75 – 80$): Intermediate skull shape.
• Brachycephalic ($> 80$): Short, broad, globular skull (e.g., modern Homo sapiens).
- Instruments: Digital spreading calipers (0-300 mm, 0.01 mm precision), sliding vernier calipers, coordinate measuring grid, and osteometric boards.
- Specimens: Comparative series of adult primate/hominin cranial casts: Chimpanzee (Pan troglodytes), Australopithecus afarensis, Homo erectus, and modern Homo sapiens.
- Cranial Landmark Identification: Locate standard anthropological landmarks: Glabella (most anterior mid-sagittal point on frontal bone), Opisthocranion (most posterior point on occiput), Euryon (most lateral points on parietal bones), Basion (anterior margin of foramen magnum), and Prosthion.
- Linear Caliper Measurements: Measure Maximum Cranial Length ($L$, Glabella-Opisthocranion), Maximum Cranial Breadth ($B$, Euryon-Euryon), Facial Breadth (Bizygomatic diameter), and Cranial Vault Height ($H$, Basion-Bregma).
- Data Transformation: Log-transform measurements: $log_{10}(L)$ and $log_{10}(B)$.
- Regression Analysis: Perform Ordinary Least Squares (OLS) regression: $log_{10}(B) = log_{10}(b) + a log_{10}(L)$. Determine slope $a$. Test $H_0: a = 1.0$ using a Student’s $t$-test.
- Index Calculations: Calculate the Cranial Index and Facial Index for each taxon.
| Taxon | Cranial Length (mm) | Cranial Breadth (mm) | Cranial Index | Cranial Morphology Classification |
|---|---|---|---|---|
| Pan troglodytes (Chimp) | 132.5 | 94.2 | 71.1 | Dolichocephalic (Elongated) |
| Australopithecus afarensis | 145.0 | 105.8 | 73.0 | Dolichocephalic |
| Homo erectus (Peking) | 194.0 | 142.5 | 73.5 | Dolichocephalic (Platycephalic) |
| Homo sapiens (Modern) | 182.0 | 144.8 | 79.6 | Mesocephalic / Globular |
Allometric Regression Analysis: Regression slope $a = 1.34 pm 0.12$ ($p < 0.01$). The cranial vault breadth displays positive hyperallometry relative to facial length across hominin phylogeny, documenting the evolutionary neurocranial expansion and globularization characteristic of encephalization in the genus Homo.
Simulating Reproductive Isolation & Speciation Modes: Allopatric Vicariance vs Dobzhansky-Muller Incompatibilities (DMI)
Speciation occurs when gene flow between diverging demes is permanently interrupted by the evolution of reproductive isolation barriers:
- Allopatric Speciation Modes:
• Vicariance: Extrinsic geological fracturing (mountain uplift, continental drift) subdivides a continuous population into two large allopatric demes.
• Peripatric (Founder Effect): A small peripheral colony buds off, undergoing rapid genetic drift and founder-induced peak shifts. - Darwin’s Dilemma of Hybrid Sterility: Darwin noted that natural selection cannot directly favor sterility or inviability in hybrids, as sterility confers zero fitness. How do reproductive barriers arise if intermediate states are unfit?
- Dobzhansky-Muller Incompatibility (DMI) Resolution:
1. The ancestral population has genotype $A_0A_0 B_0B_0$.
2. Geographic isolation divides the species into Deme 1 and Deme 2.
3. In Deme 1, a new allele $A_1$ arises and fixes ($A_1A_1 B_0B_0$). $A_1$ is fully compatible with $B_0$.
4. In Deme 2, a new allele $B_1$ arises and fixes ($A_0A_0 B_1B_1$). $B_1$ is fully compatible with $A_0$.
5. Neither deme ever experiences a drop in fitness during divergence.
6. When secondary contact occurs, hybridization produces the $F_1$ genotype $A_1A_0 B_1B_0$. The untested combination of $A_1$ and $B_1$ causes severe deleterious negative epistasis, resulting in hybrid inviability or sterility. - Haldane’s Rule (1922): “When in the $F_1$ offspring of two different animal races one sex is absent, rare, or sterile, that sex is the heterozygous [heterogametic, $XY$ or $ZW$] sex.” Explained by the dominance theory of DMIs: recessive X-linked incompatibilities are fully expressed in the heterogametic sex ($XY$ males in mammals/Drosophila; $ZW$ females in birds/Lepidoptera).
- Model System:Drosophila melanogaster × Drosophila simulans inter-specific hybrid lethality system ($Hmr$ • $Lhr$ incompatibility).
- In Silico Simulation: Python / R simulation modeling 2-locus epistatic fitness matrices under divergent selection and gene flow.
- DMI Epistatic Fitness Matrix Construction: Define fitnesses across 9 two-locus genotypes. Set $w(A_0A_0 B_0B_0) = 1.0$, $w(A_1A_1 B_0B_0) = 1.0$, $w(A_0A_0 B_1B_1) = 1.0$. Set the incompatible interaction $w(A_1A_0 B_1B_0) = 0.0$ (complete $F_1$ hybrid lethality).
- Gene Flow Migration Simulation ($m$): Introduce migration rate $m$ ($0.001$ to $0.1$) between demes. Determine critical migration threshold $m_{text{crit}}$ above which gene flow swamps divergence and halts speciation.
- Simulating the Drosophila Hmr-Lhr Lethality: Model the interaction between the D. melanogaster X-linked Hybrid male rescue ($Hmr$) gene and the D. simulans autosomal Lethal hybrid rescue ($Lhr$) gene. Observe the exact molecular mechanism of Haldane’s rule: male hybrids ($XY$) die during larval-pupal transition due to mitotic catastrophe and heterochromatin decompaction.
- Secondary Contact & Reinforcement: Simulate post-zygotic inviability causing natural selection to favor pre-zygotic assortative mating (reinforcement / Wallace effect). Plot pre-zygotic isolation index across 100 generations of sympatry.
| Cross Configuration | $F_1$ Female Viability | $F_1$ Male Viability | Phenotypic Outcome & Haldane’s Rule Check |
|---|---|---|---|
| D. melanogaster ♀ × D. simulans ♂ | 100% Viable (Sterile) | 0% (100% Lethal) | Haldane’s Rule Confirmed: Heterogametic $XY$ males completely inviable. |
| D. simulans ♀ × D. melanogaster ♂ | 0% (100% Lethal) | 100% Viable (Sterile) | Reciprocal lethal effect governed by maternal cytoplasm / nucleocytoplasmic clash. |
| Rescue Cross: melanogaster ($Hmr^-$) × simulans | 100% Viable (Sterile) | 100% Viable (Sterile) | Knockout of $Hmr$ completely suppresses epistatic lethality, rescuing male hybrids. |
Speciation Conclusion: Post-zygotic reproductive isolation evolves as an incidental by-product of independent genomic divergence. The Dobzhansky-Muller mechanism elegantly resolves Darwin’s conundrum by demonstrating how hybrid sterility arises without natural selection ever having to select for inviability.
Comprehensive Practical Viva-Voce Examination & Evolutionary Assessment Bank
Essential questions testing understanding of Hardy-Weinberg testing, selection kinetics, Wright-Fisher drift, molecular clock calibration, allometric scaling, and Dobzhansky-Muller speciation genetics:
Ans: 1. Infinitely large population size; 2. Completely random mating (panmixia); 3. No mutation; 4. No migration (gene flow); 5. No differential natural selection. Under these conditions, allele and genotype frequencies remain constant indefinitely.
Ans: Degrees of freedom df = (number of genotypic classes) – (number of alleles) = 3 – 2 = 1. Since p + q = 1, estimating one allele frequency automatically constrains the other.
Ans: As the recessive allele becomes rare (q << 0.01), over 99% of recessive alleles reside in heterozygous carriers (Aa), where they are shielded from phenotypic exposure and negative selection by the dominant allele.
Ans: Overdominance (heterozygote advantage) occurs when the heterozygote has higher relative fitness than either homozygote (w_Aa > w_AA and w_Aa > w_aa). A classic example is sickle-cell trait (HbA/HbS), which confers selective resistance against severe malaria while avoiding fatal sickle cell anemia.
Ans: q_hat = s_1 / (s_1 + s_2), where s_1 is the selection coefficient against the homozygous dominant genotype and s_2 is the selection coefficient against the homozygous recessive genotype.
Ans: Genetic drift is the non-directional random change in allele frequencies across generations due to sampling error in gamete transmission. Its magnitude is inversely proportional to population size (Var(Delta q) = p·q / (2N)); small demes experience rapid fixation or loss.
Ans: Census size N is the total head count of living individuals. Effective population size Ne is the number of breeding individuals in an idealized Wright-Fisher population that would experience the same rate of genetic drift or inbreeding.
Ans: Ne = (4 · Nm · Nf) / (Nm + Nf). If there are few males (e.g., 2 males and 98 females), all offspring inherit 50% of their genes from the 2 males, causing a genetic bottleneck with Ne ~ 7.8 individuals.
Ans: Ne across generations is calculated by the harmonic mean: 1/Ne = (1/t) · sum(1/N_i). Because harmonic means are dominated by small values, a single bottleneck generation (e.g., N = 10) permanently depresses long-term effective size.
Ans: Neutral amino acid or nucleotide substitutions accumulate in diverging lineages at an approximately constant rate over geological time, allowing sequence divergence (d) to date the time of common ancestry (T = d / 2mu).
Ans: Kimura postulated that the vast majority of evolutionary changes at the molecular level are caused by random genetic drift of selectively neutral or nearly neutral mutant alleles, rather than by positive Darwinian natural selection.
Ans: An omega ratio > 1 indicates that non-synonymous substitution rate exceeds synonymous substitution rate, which is the definitive molecular signature of positive (diversifying) Darwinian selection actively driving adaptive protein innovations.
Ans: Purifying (negative) selection. Deleterious non-synonymous amino acid changes are actively purged to preserve critical tertiary protein structure and enzyme function (typical for core histones and metabolic enzymes).
Ans: y = b · x^a, which in logarithmic form is: log(y) = log(b) + a · log(x), where y is organ/trait size, x is reference body size, a is the allometric growth exponent, and b is the initial growth coefficient.
Ans: Isometry: a = 1.0 (proportionate geometric scaling). Hyperallometry (positive allometry): a > 1.0 (trait grows disproportionately faster than body, e.g., antlers). Hypoallometry (negative allometry): a < 1.0 (trait grows disproportionately slower, e.g., brain volume).
Ans: Dolichocephaly (a long, narrow cranial vault), characteristic of non-human primates and early archaic hominins such as Homo erectus.
Ans: A two-locus epistatic model where two isolated populations independently fix compatible mutations (A_1 in Deme 1, B_1 in Deme 2). When they interbreed, the novel combination A_1 and B_1 causes deleterious negative epistasis in F_1 hybrids, producing sterility or inviability without either parental deme ever crossing a fitness valley.
Ans: When in the F_1 offspring of two different animal species one sex is absent, rare, or sterile, that sex is the heterogametic sex (XY in male mammals/Drosophila, ZW in female birds/butterflies).
Ans: The Dominance Theory: Recessive X-linked incompatibility alleles are masked by a second X chromosome in the homogametic sex (XX), but are fully expressed and unmasked in the hemizygous heterogametic sex (XY), causing inviability or sterility.
Ans: The evolution of pre-zygotic isolating mechanisms (mating calls, courtship dances, pheromones) driven by natural selection to prevent wasteful hybridization when two diverged populations with post-zygotic inviability come into secondary contact.
Explore Theoretical Evolutionary Biology Modules
Connect population genetic models, molecular clocks, and allometric scaling with natural selection, sexual selection, macroevolution, and hominin phylogeny modules.
