Coevolution of robustness, epistasis, and recombination favors asexual reproduction

Abstract

The prevalence of sexual reproduction remains one of the most perplexing phenomena in evolutionary biology. The deterministic mutation hypothesis postulates that sexual reproduction will be advantageous under synergistic epistasis, a condition in which mutations cause a greater reduction in fitness when combined than would be expected from their individual effects. The inverse condition, antagonistic epistasis, correspondingly is predicted to favor asexual reproduction. To assess this hypothesis, we introduce a finite population evolutionary process that combines a recombination modifier formalism with a gene-regulatory network model. We demonstrate that when reproductive mode and epistasis are allowed to coevolve, asexual reproduction outcompetes sexual reproduction. In addition, no correlation is found between the level of synergistic epistasis and the fixation time of the asexual mode. However, a significant correlation is found between the level of antagonistic epistasis and asexual mode fixation time. This asymmetry can be explained by the greater reduction in fitness imposed by sexual reproduction as compared with asexual reproduction. Our findings present evidence and suggest plausible explanations that challenge both the deterministic mutation hypothesis and recent arguments asserting the importance of emergent synergistic epistasis in the maintenance of sexual reproduction.

Fig. 1.

Comparison of sexual and asexual initial phases for cases where the modifier allele (r) reducing recombination fixes (dominant mode). (a and b) After the sexual initial phase, epistasis and fixation time are uncorrelated (r = −0.040, R² = 0.002) (a), whereas for the asexual initial phase, they are negatively correlated (r = −0.490, R² = 0.240) (b). (c and d) In both cases, epistasis and inclusive robustness (see Methods) are negatively correlated: sexual initial phase (r = −0.320, R² = 0.102) (c) and asexual initial phase (r = −0.466, R² = 0.217) (d).