Evolution under a model of functionally buffered deleterious mutations can lead to positive selection in protein-coding genes

Selective pressures on DNA sequences often result in departures from neutral evolution that can be captured by the McDonald-Kreitman (MK) test. However, the nature of such selective forces often remains unknown to experimentalists. Amino acid fixations driven by natural selection in protein-coding genes are commonly associated with a genetic arms race or changing biological purposes, leading to proteins with new functionality. Here, we evaluate the expectations of population genetic patterns under a buffering mechanism driving selective amino acids to fixation, which is motivated by an observed phenotypic rescue of otherwise deleterious nonsynonymous substitutions at bag of marbles (bam) and Sex lethal (Sxl) in Drosophila melanogaster. These two genes were shown to experience strong episodic bursts of natural selection potentially due to infections of the endosymbiotic bacteria Wolbachia observed among multiple Drosophila species. Using simulations to implement and evaluate the evolutionary dynamics of a Wolbachia buffering model, we demonstrate that selectively fixed amino acid replacements will occur, but that the proportion of adaptive amino acid fixations and the statistical power of the MK test to detect the departure from an equilibrium neutral model are both significantly lower than seen for an arms race/change-in-function model that favors proteins with diversified amino acids. We find that the observed selection pattern at bam in a natural population of D. melanogaster is more consistent with an arms race model than with the buffering model.

Automated summary

Selective pressures on DNA sequences can be captured by the McDonald-Kreitman (MK) test. The nature of these selective forces often remains unknown. Amino acid fixations driven by natural selection in genes are associated with genetic arms races or changes in biological purposes. Simulations reveal that the proportion of adaptive amino acid fixations and the power of the MK test are lower in a buffering model compared to an arms race/change-in-function model.