A SIMULATION STUDY OF MULTILOCUS CLINES

Fisher's method of junctions is used to investigate the degree of association between selected alleles in a cline, in the limit where there is divergence between very many genes. A computer model is used to simulate one of a pair of infinite demes that exchange individuals each generation. Selection is on haploids; it is additive and is equivalent to heterozygote disadvantage. Recombination is uniform over a single chromosome. A "critical value" of selection exists at equilibrium, below which loci act independently and above which they act in association (Barton 1983). Starting with secondary contact, simulation results contrast markedly with the equilibrium solution. The "critical value" is not apparent in the simulated clines, even after many generations. Rather, loci remain associated to some extent under all degrees of selection. The simulation is consistent with the equilibrium analysis in all other respects, and therefore indicates that under weak selection the approach to equilibrium is very slow. This is borne out by further numerical calculations. The slow approach to equilibrium enables us to estimate the time since contact between two demes under idealized conditions. Extending this work toward natural hybrid zones is discussed.