Experiment 04Molecular EvolutionDivergence Dating
Molecular Clock & Divergence Dating: Estimating Speciation Times & $dN/dS$ ($omega$) Selection Ratios
Duration: 4 Hours • Computational Molecular Lab
🎯 Experimental Aim & Learning Objectives
To calibrate a linear molecular clock using fossil divergence benchmarks and protein/nucleotide sequence substitution rates ($mu$); estimate divergence dates for un-dated phylogenetic nodes; and calculate codon non-synonymous ($dN$) to synonymous ($dS$) substitution rate ratios ($omega = dN/dS$) to distinguish purifying, neutral, and positive Darwinian selection.
🔬 Theoretical Foundation & Molecular Clock Equations
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.
Figure 1.4: Linear molecular clock divergence dating ($T = d / 2mu$) and selective constraint regimes defined by the non-synonymous to synonymous substitution ratio ($dN/dS$ or $omega$).
🧪 Computational Tools & Benchmark Data
- 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).
📋 Step-by-Step Computational Protocol
- 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$).
📊 Representative Molecular Clock & Selection Output
| 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.
{reading_nav}
