Selection-recombination-mutation dynamics: Gradient, limit cycle, and closed invariant curve

In this paper, the replicator dynamics of the two-locus two-allele system under weak mutation and weak selection is investigated in a generation-wise non-overlapping unstructured population of individuals mating at random. Our main finding is that the dynamics is gradient-like when the point mutations at the two loci are independent. This is in stark contrast to the case of one-locus multi-allele where the existence gradient behaviour is contingent on a specific relationship between the mutation rates. When the mutations are not independent in the two-locus two-allele system, there is the possibility of non-convergent outcomes, like asymptotically stable oscillations, through either the Hopf bifurcation or the Neimark–Sacker bifurcation depending on the strength of the weak selection. The results can be straightforwardly extended for multi-locus two-allele systems.