Substitution Rate Calculators

0 calculators tagged with “Substitution Rate”

The substitution rate (also called evolutionary rate) is the rate at which nucleotides or amino acids are replaced over evolutionary time — measured as substitutions per site per year. It is the foundation of the molecular clock, which uses substitution rates to estimate divergence times between lineages. Rates vary enormously: RNA viruses evolve at ~10⁻³–10⁻⁴ substitutions/site/year (very fast); mitochondrial DNA ~2% per million years; nuclear genes ~10⁻⁹ per site per year. Under Kimura's neutral theory, most nucleotide substitutions are neutral — the rate equals the mutation rate for neutral sites.

All Calculators

No calculators found for this topic.

Substitution Rate Formula

r = d / (2t)

r = substitution rate (substitutions/site/year); d = number of substitutions per site between two lineages (pairwise divergence); t = time since divergence (years). Divide by 2 because divergence accumulates on both lineages. Example: two genes with 3% divergence (d = 0.03) and a fossil calibration of 15 million years since divergence: r = 0.03 / (2 × 15 × 10⁶) = 1.0 × 10⁻⁹ substitutions/site/year.

Neutral Theory and Molecular Clock

Kimura's neutral theory (1968): most mutations at the molecular level are neutral or nearly neutral — they do not affect fitness. Under neutrality: rate of substitution = rate of neutral mutation = μ_neutral. The molecular clock: if the substitution rate is constant over time, genetic divergence is proportional to time → use divergence to date evolutionary events when calibrated with fossils or biogeography.

Variation in Substitution Rates

  • Synonymous (silent) substitutions: 3rd codon position; not change amino acid; highest rate; neutral; ~1.5–5 × 10⁻⁹/site/year (mammals)
  • Non-synonymous substitutions: change amino acid; subject to selection; lower rate
  • Purifying selection: reduces non-synonymous rate at conserved positions
  • Positive selection: accelerates substitutions at positively selected sites (dN/dS > 1)

dN/dS Ratio

ω = dN/dS (Ka/Ks): dN = non-synonymous substitution rate; dS = synonymous substitution rate. ω < 1: purifying selection (most protein-coding genes). ω = 1: neutral evolution. ω > 1: positive selection (adaptive evolution).

Glossary

Substitution Rate
r = d/(2t); the rate of nucleotide or amino acid replacement per site per year; foundation of the molecular clock; varies from ~10⁻³/site/year (RNA viruses) to ~10⁻⁹/site/year (nuclear DNA).
Molecular Clock
The use of constant substitution rates to estimate divergence times from genetic differences; calibrated with fossil ages; relaxed clock models allow rate variation between lineages.
dN/dS (ω)
Ratio of non-synonymous to synonymous substitution rates; ω < 1 = purifying selection; ω = 1 = neutral; ω > 1 = positive selection; used to identify adaptively evolving genes.

Frequently Asked Questions

The substitution rate (r) = rate at which nucleotides or amino acids are replaced per site per year. Calculation: r = d / (2t), where d = divergence (substitutions per site between two sequences) and t = time since common ancestor (years). Factor of 2 because substitutions accumulate on both lineages since the split. Example: human and chimpanzee cytochrome b genes show ~12% divergence (d = 0.12); human-chimp divergence ~6 million years ago: r = 0.12/(2 × 6 × 10⁶) = 10⁻⁸ substitutions/site/year. Divergence is estimated from aligned sequences using substitution models (JC69, HKY85, GTR) that correct for multiple hits.

The molecular clock hypothesis (Zuckerkandl and Pauling, 1965) proposes that molecular sequences evolve at approximately constant rates over time. If true: genetic divergence is proportional to divergence time → sequences can be used as clocks to date evolutionary events. Support: many gene families show roughly linear accumulation of substitutions over geological time. Problems: rate variation exists between lineages (generation time effect; metabolic rate effect; effective population size effects) — the molecular clock is not strictly constant. Relaxed molecular clock models (Bayesian: BEAST software) allow rates to vary between lineages while still using divergence to estimate timing. Calibration: fossil ages or biogeographic events anchor the clock by providing known divergence times for rate calculation.

dN/dS (ω or Ka/Ks) compares the rate of non-synonymous substitutions (dN, changing amino acid) to synonymous substitutions (dS, silent at protein level). dS ≈ neutral rate; dN depends on selective constraint. ω < 1: purifying (negative) selection — most amino acid changes are harmful → selected against → dN < dS. Typical for conserved genes (histones ω ≈ 0.001). ω = 1: neutral evolution — amino acid changes are neither beneficial nor harmful. ω > 1: positive (Darwinian) selection — amino acid changes are beneficial → selected for → dN > dS. Found at sites of immune-parasite coevolution, reproductive proteins, adaptive loci. Detecting positive selection: PAML (codeml), HyPhy, dN/dS-based methods scan genomes for adaptively evolving genes.

RNA virus evolution is ~10⁶ times faster than nuclear DNA genomes: ~10⁻³–10⁻⁴ substitutions/site/year for RNA viruses vs. ~10⁻⁹/site/year for nuclear DNA. Reasons: (1) Error-prone replication: RNA-dependent RNA polymerases (RdRp) lack proofreading activity → 10⁻⁴–10⁻⁵ errors per site per replication (vs. 10⁻⁸–10⁻¹⁰ for DNA with proofreading). (2) High replication rate: RNA viruses replicate millions of times per day in a host. (3) Large population size: huge numbers of viral particles increase the rate at which beneficial mutations are fixed. (4) Short generation time: more replication cycles per year → more substitutions per year. These rapid evolution rates allow RNA viruses to rapidly evade immunity, develop drug resistance, and cause annual flu pandemics requiring updated vaccines.