The main findings of a study of the evolution of modifier gene frequencies in models of deterministic population genetics are presented. A wide variety of random mating systems are subject to selection with modifiers operating, in different cases, on mutation rates, dominance, migration between subpopulations, and linkage between other loci. In all these instances the modifier frequencies evolve in such a way as to maximize the mean fitness of the population at equilibrium. This is remarkable since, except for the dominance modifier, the modifier genes are selectively neutral in the sense that they do not affect the fitness of their individual carriers. In nonrandom mating systems the mean fitness concept is not well defined, and there does not appear to be such a simple principle governing the evolution of modifier frequencies. In assortative mating systems modifiers favoring reduced assortment propensities tend to increase. In contrast, for selfing-outcrossing systems modifiers favoring increased selfing tend to increase.