The number of self-incompatibility alleles in a finite, subdivided population

The actual and effective number of gametophytic self-incompatibility alleles maintained at mutation-drift-selection equilibrium in a finite population subdivided as in the island model is investigated by stochastic simulations. The existing theory founded by Wright predicts that for a given population size the number of alleles maintained increases monotonically with decreasing migration as is the case for neutral alleles. The simulation results here show that this is not true. At migration rates above Nm = 0.01-0.1, the actual and effective number of alleles is lower than for an undivided population with the same number of individuals, and, contrary to Wright's theoretical expectation, the number of alleles is not much higher than for an undivided population unless Nm < 0.001. The same pattern is observed in a model where the alleles display symmetrical overdominant selection. This broadens the applicability of the results to include proposed models for the major histocompatibility (MHC) loci. For a subdivided population over a large range of migration rates, it appears that the number of self-incompatibility alleles (or MHC-alleles) observed can provide a rough estimate of the total number of individuals in the population but it underestimates the neutral effective size of the subdivided population.