Inferring the number of evolutionary events from DNA coding sequence differences

The estimation of the amount of evolutionary divergence that has taken place between two DNA coding sequences depends strongly on the degree of constraint on amino acid replacements. If amino acid replacements are relatively unconstrained, the individual nucleotide is the appropriate unit of analysis and the method of Tajima and Nei can be used. If amino acid replacements are constrained, however, this method is shown to be inapplicable. For sequences with strong amino acid constraints, a method is outlined analogous to the Tajima and Nei method using codons as the unit of analysis. Only synonymous substitutions are used. Codon usage data can be employed to estimate the necessary parameters of the calculation, or a priori models of substitution may be employed. Sequences with significant but intermediate constraints on amino acid replacements are, in principle, unanalyzable.