Dynamics of correlated genetic systems. II. Simulation studies of chromosomal segments under selection

Abstract

The dynamics of chromosomal segments under selection are investigated by comparing experimental data to simulations of simple models of selection. The simulations assume 93 loci distributed evenly along an entire chromosome. The two issues addressed in this paper concern rates of decay of linkage disequilibria for chromosomes under selection and rates of gene frequency change after perturbation of gametic frequencies to states near the edge of the gametic frequency simplex. The findings are: (1) If reasonable values of inbreeding depression are assumed, linkage disequilibria decays to zero but at a rate nearly twice that expected from neutral theory. Experimental results also show accelerated decay rates. The acceleration of decay seems to be a simple consequence of the increased heterozygosity produced by selection. It is, therefore, argued that massive linkage disequilibria, of the kind found by Franklin and Lewontin (1970) in their simulations, are unlikely to characterize the genetic structure of natural populations of random mating organisms. (2) It is possible to distinguish between two time-honored models of multilocus selection, known as the symmetric overdominant and classical models, on the basis of gene frequency change near the edge of the gametic frequency simplex assuming linkage disequilibria is intense. (3) Examination of experimental data from perturbation experiments shows that neither of these elementary models adequately account for observed rates of gene frequency change, although the symmetric overdominant model does provide the better fit. Instead the experimental data suggest a markedly nonuniform distribution of selective effects along the chromosome. The data also suggest that these selective effects combine in markedly nonadditive ways in determining joint fitness.