Changing Population Size in McDonald-Kreitman Style Analyses: Artifactual Correlations and Adaptive Evolution between Humans and Chimpanzees

It is known that methods to estimate the rate of adaptive evolution, which are based on the McDonald-Kreitman test, can be biased by changes in effective population size. Here, we demonstrate theoretically that changes in population size can also generate an artifactual correlation between the rate of adaptive evolution and any factor that is correlated to the strength of selection acting against deleterious mutations. In this context, we have investigated whether several site-level factors influence the rate of adaptive evolution in the divergence of humans and chimpanzees, two species that have been inferred to have undergone population size contraction since they diverged. We find that the rate of adaptive evolution, relative to the rate of mutation, is higher for more exposed amino acids, lower for amino acid pairs that are more dissimilar in terms of their polarity, volume, and lower for amino acid pairs that are subject to stronger purifying selection, as measured by the ratio of the numbers of nonsynonymous to synonymous polymorphisms (p(N)/p(S)). All of these correlations are opposite to the artifactual correlations expected under contracting population size. We therefore conclude that these correlations are genuine.

Automated summary

In the divergence of humans and chimpanzees, the rate of adaptive evolution is found to be negatively correlated with the dissimilarity of amino acid pairs in terms of polarity and volume, as well as the strength of purifying selection. These correlations are not affected by population size contraction.