The cytonuclear effects of facultative apomixis. I. Disequilibrium dynamics in diploid populations.

We comprehensively analyze the cytonuclear effects of generalized mixed mating, including all combinations of selfing, outcrossing, and apomixis, the asexual production of seeds. After first deriving the time-dependent solutions for nonrandom associations (disequilibria) between a diallelic cytoplasmic marker and the alleles and genotypes at a diploid nuclear locus, we delimit all possible dynamical behaviors and the conditions under which each occurs. As in standard mixed mating systems, all disequilibria ultimately decay to zero except when outcrossing is absent, in which case permanent disequilibria result if the allelic association is initially nonzero. When at least some outcrossing is present, any initial allelic association decays at a constant geometric rate, whereas genotypic disequilibria may first increase in magnitude or change sign. Although selfing and apomixis tend to retard the decay of disequilibria (or approach to equilibrium) and often to the same extent, apomixis can have a stronger effect under some conditions. We also determine the dynamics of cytonuclear disequilibria in specific examples that may be of particular interest for empirical studies of hybrid zones. The results suggest several practical guidelines for experimental design and data analysis and show how the cytonuclear disequilibrium dynamics under mating system alone furnish a quantitative baseline for null hypotheses against which to test for the presence of other evolutionary forces.

The cytonuclear effects of facultative apomixis. II. Definitions and dynamics of disequilibria in tetraploid populations.

We develop a cytonuclear framework for tetraploid populations in which a diallelic nuclear marker exhibits tetrasomic inheritance. This system requires two separate parameterizations, with six cytonuclear disequilibria (nonrandom associations) in tetraploid individuals and four in their diploid gametes. Double reduction during meiosis adds further complexity by causing gametic output to vary with the distance of the nuclear locus from the centromere. We derive and analyze dynamical solutions for the disequilibria under generalized mixed mating, with any combination of apomixis, selfing, and outcrossing, with and without double reduction. As in comparable diploid systems, all disequilibria ultimately decay to zero, unless nuclear and cytoplasmic alleles are nonrandomly associated and outcrossing is absent, in which case permanent associations result. Selfing and apomixis retard the decay of disequilibria (or approach to equilibrium), and often to the same extent. In contrast, double reduction can accelerate the loss of tetraploid cytonuclear associations, but only negligibly in hybrid zones, and this loss is never faster than in diploids. Only in the absence of allelic associations or outcrossing is the asymptotic approach to equilibrium differentially affected by apomixis and selfing or slower under tetrasomic than disomic inheritance. To facilitate empirical applications, we also examine tetraploid hybrid zone dynamics and offer practical guidelines for experimental design and data analysis, showing how the consequences of the mating system alone provide a valuable baseline for drawing evolutionary inferences from the observed patterns of cytonuclear associations.

Genetic diversity at a single locus under viability selection and facultative apomixis: equilibrium structure and deviations from Hardy-Weinberg frequencies.

We extensively analyze the maintenance of genetic variation and deviations from Hardy-Weinberg frequencies at a diallelic locus under mixed mating with apomixis and constant viability selection. Analytical proofs show that: (1) at most one polymorphic equilibrium exists, (2) polymorphism requires overdominant or underdominant selection, and (3) a simple, modified overdominance condition is sufficient to maintain genetic variation. In numerical analyses, only overdominant polymorphic equilibria are stable, and these are stable whenever they exist, which happens for approximately 78% of random fitness and mating parameters. The potential for maintaining both alleles increases with increasing apomixis or outcrossing and decreasing selfing. Simulations also indicate that equilibrium levels of heterozygosity will often be statistically indistinguishable from Hardy-Weinberg frequencies and that adults, not seeds, should usually be censused to maximize detecting deviations. Furthermore, although both censuses more often have an excess rather than a deficit of heterozygotes, analytical sign analyses of the fixation indices prove that, overall, adults are more likely to have an excess and seeds a deficit at equilibrium.