Subterranean archipelago: mitochondrial DNA phylogeography of stygobitic isopods (Oniscidea�:�Haloniscus) from the Yilgarn region of Western Australia

DNA sequence evolution through nucleotide substitution may be assimilated to a stationary Markov process. The fundamental equations of the general model, with 12 independent substitution parameters, are used to obtain a formula which corrects the effect of multiple and parallel substitutions on the measure of evolutionary divergence between two homologous sequences. We show that only reversible models, with six independent parameters, allow the calculation of the substitution rates. Simulation experiments on DNA sequence evolution through nucleotide substitution call into question the effectiveness of the general model (and of any other more detailed description); nevertheless, the general model results are slightly superior to any of its particular cases.

Although several decades of study have revealed the ubiquity of variation of evolutionary rates among sites, reliable methods for studying rate variation were not developed until very recently. Early methods fit theoretical distributions to the numbers of changes at sites inferred by parsimony and substantially underestimate the rate variation. Recent analyses show that failure to account for rate variation can have drastic effects, leading to biased dating of speciation events, biased estimation of the transition:transversion rate ratio, and incorrect reconstruction of phylogenies.