Molecular clock in neutral protein evolution.

BACKGROUND

A frequent observation in molecular evolution is that amino-acid substitution rates show an index of dispersion (that is, ratio of variance to mean) substantially larger than one. This observation has been termed the overdispersed molecular clock. On the basis of in silico protein-evolution experiments, Bastolla and coworkers recently proposed an explanation for this observation: Proteins drift in neutral space, and can temporarily get trapped in regions of substantially reduced neutrality. In these regions, substitution rates are suppressed, which results in an overall substitution process that is not Poissonian. However, the simulation method of Bastolla et al. is representative only for cases in which the product of mutation rate micro and population size Ne is small. How the substitution process behaves when micro Ne is large is not known.

RESULTS

Here, I study the behavior of the molecular clock in in silico protein evolution as a function of mutation rate and population size. I find that the index of dispersion decays with increasing micro Ne, and approaches 1 for large micro Ne. This observation can be explained with the selective pressure for mutational robustness, which is effective when micro Ne is large. This pressure keeps the population out of low-neutrality traps, and thus steadies the ticking of the molecular clock.

CONCLUSIONS

The molecular clock in neutral protein evolution can fall into two distinct regimes, a strongly overdispersed one for small micro Ne, and a mostly Poissonian one for large micro Ne. The former is relevant for the majority of organisms in the plant and animal kingdom, and the latter may be relevant for RNA viruses.