The two-locus model of Gaussian stabilizing selection

We study the equilibrium structure of a well-known two-locus model in which two diallelic loci contribute additively to a quantitative trait that is under Gaussian stabilizing selection. The population is assumed to be infinitely large, randomly mating, and having discrete generations. The two loci may have arbitrary effects on the trait, the strength of selection and the recombination rate may also be arbitrary. We find that 16 different equilibrium patterns exist, having up to 11 equilibria; up to seven interior equilibria may coexist, and up to four interior equilibria, three in negative and one in positive linkage disequilibrium, may be simultaneously stable. Also, two monomorphic and two fully polymorphic equilibria may be simultaneously stable. Therefore, the result of evolution may be highly sensitive to perturbations in the initial conditions or in the underlying genetic parameters. For the special case of equal effects, global stability results are proved. In the general case, we rely in part on numerical computations. The results are compared with previous analyses of the special case of extremely strong selection, of an approximate model that assumes linkage equilibrium, and of the much simpler quadratic optimum model.