Solution

ID: past-exam-of-the-biological-physics-course-of-the-university-of-oxford/2021/6/b/solution

A point mutation changes an individual nucleotide position, introducing heritable sequence variation. Replication errors and spontaneous chemical damage are possible origins; mutation does not occur because a future change would be useful. Natural selection changes the frequencies of variants through their effects on reproduction, while genetic drift causes random frequency changes, especially in small populations. A sequence difference fixed between lineages is a substitution; many new mutations disappear and never become substitutions. Under a simple neutral diploid model, new neutral variants per generation times fixation probability gives substitution rate , motivating a molecular clock for suitable neutral sites.
To infer a phylogenetic tree, compare homologous, preferably orthologous sequences and align corresponding sites. A nucleotide substitution model converts observed sequence differences into estimated evolutionary branch lengths, allowing for multiple changes at one site. For example, in the equal-frequency, equal-rate Jukes-Cantor model, a site changes to any alternative with rate , where is its total substitution rate. The probability of matching its initial state after elapsed time is . Thus for an evolutionary separation with expected substitutions per site , the observed mismatch fraction is and
For small , ; near saturation , large evolutionary distances become difficult to distinguish. The correction does not rescue fully saturated data. Real analyses may allow unequal nucleotide frequencies, transition biases and variation of rates between sites.
Tree topology and branch lengths can be estimated using maximum likelihood estimation or Bayesian phylogenetic inference, with support assessed from resampling or posterior uncertainty. A gene tree need not equal the species tree, because incomplete lineage sorting, genetic recombination and horizontal gene transfer can produce different histories. Rooting with an appropriate outgroup establishes ancestral direction; pairwise differences alone do not supply a root.
Branch lengths measured in substitutions per site are not yet calendar times. Under a strict molecular clock with rate per site per year, a branch of duration has length ; two contemporaneously sampled descendants separated at time have total path length , hence . Calibration comes from independently dated constraints such as fossil calibration, geological information or sufficiently informative sampling dates. A dated common ancestor is not automatically identical in age to the earliest observed fossil. A relaxed molecular clock permits branch-specific rates while imposing a statistical model on them.
Without an external temporal calibration or an independently known rate, multiplying all times by a constant and dividing all rates by that constant leaves expected sequence divergence unchanged. Absolute dates are therefore unidentifiable from sequence divergence alone. A useful time-calibrated phylogenetic tree reports uncertainty from sequence sampling, substitution models, rate variation, topology and calibration bounds, rather than displaying an unjustified exact date for each ancestor.

New to topics? Read the docs here!