The effects of assortative mating and migration on cytonuclear associations in hybrid zones.

We examine the influence of nonrandom mating and immigration on the evolutionary dynamics of cytonuclear associations in hybrid zones. Recursion equations for allelic and genotypic cytonuclear disequilibria were generated under models of (1) migration alone, assuming hybrid zone matings are random with respect to cytonuclear genotype; and (2) migration in conjunction with refined epistatic mating, in which females of the pure parental species preferentially mate with conspecific males. Major results are as follows: (a) even the slightest migration removes the dependency of the final outcome on initial conditions, producing a unique equilibrium in which both pure parental genotypes are maintained in the hybrid zone; (b) in contrast to nuclear genes, the dynamics of cytoplasmic allele frequencies appear robust to changes in the assumed mating system, yet are particularly sensitive to gene flow; (c) continued immigration can generate permanent cytonuclear disequilibria, whether mating is random or assortative; and (d) the order of population censusing (before versus after reproduction by immigrants) can have a dramatic effect on the magnitude but not the pattern of cytonuclear disequilibria. Using the maximum likelihood method, the parameter space of migration rates and assortative mating rates was examined for best fit to observed cytonuclear disequilibria data in a hybrid population of Hyla tree frogs. An epistatic mating model with a total immigration rate of about 32% per generation produces equilibrium gene frequencies and cytonuclear disequilibria consistent with the empirical observations.