Estimating Substitution Rates from Molecular Data Using the Coalescent

A coalescent model is used to estimate the rate at which neutral substitutions occur in a DNA sequence, without the necessity for an independent estimate of divergence times. Given a random sample of molecular sequences from a finite population, the distribution of the time to a common ancestor can be obtained from the coalescent model. With this principle, summary statistics are developed that use the distribution of molecular diversity within the sample to estimate the relative magnitude of nucleotide substitution rates. If, in addition, the effective population size is known, absolute substitution rates can also be estimated. These techniques are illustrated by estimating the transition rates that underlie the evolution of the first 360 nucleotides of the mitochondrial control region in an Amerindian tribal population.