Continuous selective models with mutation and migration

The continuous selective model formulated previously for a single locus with multiple alleles in a monoecious population is extended to include mutation and migration. Somatic and germ line genotypic frequencies are distinguished, and the alternative hypotheses of constant mutation rates and age-independent mutation frequencies are analyzed in detail for arbitrary selection and mating schemes. With any mating pattern, if there is no selection, the equilibrium allelic frequencies are shown to be unaffected by the generalizations introduced in this paper. If, in addition, mating is at random, the equilibrium genotypic frequencies are proved to be in Hardy-Weinberg proportions. For both models, the nature of the approach to equilibrium is discussed. Migration is treated in the island model.     

Separable models for age-structured population genetics

This paper is concerned with the applications of nonlinear age-dependent dynamics to population genetics. Age-structured models are formulated for a single autosomal locus with an arbitrary number of alleles. The following cases are considered: a) haploid populations with selection and mutation; b) monoecious diploid populations with or without mutation reproducing by self-fertilization or by two types of random mating. The diploid models do not deal with selection. For these cases the genic and genotypic frequencies evolve towards time-persistent forms, whether the total population size tends towards exponential growth or not.     

The effect of positive assortative mating at one locus on a second linked locus

Summary Considerations proceed from a model of positive assortative mating based on genotype at one locus, with an arbitrary number of alleles, assuming no selection, mutation, or migration, hypothetically infinite population size, and discrete non-overlapping generations. From these conditions, inferences are made about the genotypic structure at a linked locus, as well as about the corresponding 2-locus gametic structure.The following main results are presented: in the course of the generations, the genotypic structure at the second locus and the 2-locus gametic structure always tend to a limit responsive to the initial conditions concerning the joint genotypic structure at the two loci and the degree of assortativity and linkage. A complete, analytical representation of the limits is given. In particular, if assortative mating is only partial and at the same time linkage is not complete, a population is not able to maintain a permanent deviation of the gametic structure from linkage equilibrium, and thus the genotypic structure at the second locus tends to Hardy-Weinberg proportions. On the other hand, if initial linkage disequilibrium is combined with partial assortative mating and complete linkage (or with complete assortative mating and unlinked loci) the population maintains this disequilibrium and thus the genotypic structure at the second locus need not tend to Hardy-Weinberg proportions. It turns out that the conditions not only of complete linkage, but also of unlinked loci together with complete assortativity, imply no change in gametic structure from the initial structure.In order to demonstrate the influence of several parameters on the speed of convergence to and the magnitude of the respective limits, several graphs are included.     

A multi-locus continuous-time selection model

Summary A continuous time selection model is formulated for a diploid monoecious population with multiple alleles at each of an arbitrary number of loci, incorporating differential fertility and mortality as well as arbitrary mating and age structure. The model is simplified in the case of age-independence and for the case of a stable age distribution. The age-independent model is examined in detail for the special case of multiple alleles at each of two loci. This model is analyzed under the assumptions of random mating and additive fertilities, with close attention given to the behavior of the system with respect to Hardy-Weinberg proportions and linkage equilibrium.M. M. was supported by a U.S. Public Health Service training grant (Grant No. GM780).     

The multiple-locus symmetric fertility model

Selection due to variation in the fecundity among matings of genotypes with respect to many loci each with two alleles is studied. The fitness of a mating depends only on the genotypic distinction between homozygote and heterozygote at each locus in the two individuals, and differences among loci are allowed. This symmetric fertility model is therefore a generalization of the multiple-locus symmetric viability model. The phenomena seen in the two-locus symmetric fertility model generalize—e.g., the possibility of joint stability of equilibria with linkage equilibrium and with linkage disequilibrium, and the existence of different types of totally polymorphic equilibria with the gametic proportions in linkage equilibrium. The central equilibrium with genotypic frequencies in Hardy-Weinberg proportions and gametic frequencies in Robbins proportions exists for all symmetric fertility models. For some symmetric fertility regimes additional equilibria exist with gametic frequencies in linkage equilibrium and with genotypic frequencies in Hardy-Weinberg proportions at all except one locus. These equilibria may exist in the dioecious symmetric viability model, and then they will be locally stable. For free recombination the stable equilibria show linkage equilibrium, but several of these with different numbers of polymorphic loci may be stable simultaneously.     

Deviations of Genotypic Structures from Hardy-Weinberg Proportions under Random Mating and Differential Selection between the Sexes

Population genetic models, such as differential viability selection between the sexes and differential multiplicative fecundity contributions of the sexes, are considered for a single multiallelic locus. These selection models usually produce deviations of the zygotic genotype frequencies from Hardy-Weinberg proportions. The deviations are investigated (with special emphasis put on equilibrium states) to quantify the effect of selective asymmetry in the two sexes. For many selection regimes, the present results demonstrate a strong affinity of zygotic genotype frequencies for Hardy-Weinberg proportions after two generations, at the latest. It is shown that the deviations of genotypic equilibria from the corresponding Hardy-Weinberg proportions can be expressed and estimated by means of selection components of only that sex with the lower selection intensity. This corresponds to the well-known fact that viability selection acting in only one sex yields Hardy-Weinberg equilibria.     

Modifiers of mutation rate: a general reduction principle

A deterministic two-locus population genetic model with random mating is studied. The first locus, with two alleles, is subject to mutation and arbitrary viability selection. The second locus, with an arbitrary number of alleles, controls the mutation at the first locus. A class of viability-analogous Hardy-Weinberg equilibria is analyzed in which the selected gene and the modifier locus are in linkage equilibrium. It is shown that at these equilibria a reduction principle for the success of new mutation-modifying alleles is valid. A new allele at the modifier locus succeeds if its marginal average mutation rate is less than the mean mutation rate of the resident modifier allele evaluated at the equilibrium. Internal stability properties of these equilibria are also described.     

The effect of positive assortative mating at one locus on a second linked locus

Summary A model for positive assortative mating based on genotype for one locus is employed to investigate the effect of this mating system on the genotypic structure of a second linked locus as well as on the joint genotypic structure of these two loci. It is shown that the second locus does not attain a precise positive assortative mating structure, but yet it shares a property that is characteristic of positive assortative mating, namely an increase in the frequency of homozygotes over that typically found in panmictic structures. Given any arbitrary genotypic structure for the parental population, the resulting offspring generation possesses a structure at the second locus that does not depend on the recombination frequency, while the joint structure of course does. In case assortative mating as well as linkage are not complete, there exists a unique joint equilibrium state for the two loci, which is characterized by complete stochastic independence between the two loci as well as by Hardy-Weinberg proportions at the second locus. For the second locus alone, Hardy-Weinberg equilibrium is realized if and only if gametic linkage equilibrium and an additionally specified condition are realized.     

The Evolution of One- and Two-Locus Systems

Assuming age-independent fertilities and mortalities and random mating, continuous-time models for a monoecious population are investigated for weak selection. A single locus with multiple alleles and two alleles at each of two loci are considered. A slow-selection analysis of diallelic and multiallelic two-locus models with discrete nonoverlapping generations is also presented. The selective differences may be functions of genotypic frequencies, but their rate of change due to their explicit dependence on time (if any) must be at most of the second order in s, (i.e., O( s²)), where s is the intensity of natural selection. Then, after several generations have elapsed, in the continuous time models the time-derivative of the deviations from Hardy-Weinberg proportions is of O(s²), and in the two-locus models the rate of change of the linkage disequilibrium is of O(s²). It follows that, if the rate of change of the genotypic fitnesses is smaller than second order in s (i.e., o(s²)), then to O(s²) the rate of change of the mean fitness of the population is equal to the genic variance. For a fixed value of s, however, no matter how small, the genic variance may occasionally be smaller in absolute value than the (possibly negative) lower order terms in the change in fitness, and hence the mean fitness may decrease. This happens if the allelic frequencies are changing extremely slowly, and hence occurs often very close to equilibrium. Some new expressions are derived for the change in mean fitness. It is shown that, with an error of O( s), the genotypic frequencies evolve as if the population were in Hardy-Weinberg proportions and linkage equilibrium. Thus, at least for the deterministic behavior of one and two loci, deviations from random combination appear to have very little evolutionary significance.     

The diffusion model for migration and selection in a dioecious population

The diffusion approximation is derived for migration and selection at a multiallelic locus in a dioecious population subdivided into a lattice of panmictic colonies. Generations are discrete and nonoverlapping; autosomal and X-linked loci are analyzed. The relation between juvenile and adult subpopulation numbers is very general and includes both soft and hard selection; the zygotic sex ratio is the same in every colony. All the results hold for both adult and juvenile migration. If ploidy-weighted average selection, drift, and diffusion coefficients are used, then the ploidy-weighted average allelic frequencies satisfy the corresponding partial differential equation for a monoecious population. The boundary conditions and the unidimensional transition conditions for coincident discontinuities in the carrying capacity and migration rate extend identically. The previous unidimensional formulation and analysis of symmetric, nearest-neighbor migration of a monoecious population across a geographical barrier is generalized to symmetric migration of arbitrary finite range, and the transition conditions are shown to hold for a dioecious population. Thus, the entire theory of clines and of the wave of advance of favorable alleles is applicable to dioecious populations.This work was supported by National Science Foundation grant BSR-9006285     

Conditions for the Existence of Clines

A very general partial differential equation in space and time satisfied by the gene frequency in a monoecious population distributed continuously over an arbitrary habitat is derived. The treatment is restricted to a single diallelic locus in the absence of mutation and random drift, and it is supposed that time is continuous, births and deaths occur at random, and migration is independent of genotype. With the further assumptions that migration is isotropic and homogeneous, the population density is constant and uniform (as permitted by the population regulation mechanism included in the formulation), and Hardy-Weinberg proportions obtain locally, this partial differential equation reduces to the simplest multidimensional generalization of the classical Fisher-Haldane cline model. The efficacy of migration and selection in maintaining genetic variability at equilibrium in this model is investigated by deducing conditions for the existence of clines under various circumstances. The effects of the degree of dominance, a neutral belt between the regions where a particular allele is advantageous and deleterious, finiteness of the habitat, and habitat dimensionality are evaluated. Provided at least one of the alleles is favored only in a finite region, excluding the special case in which its total effective selective coefficient is zero, if conditions for supporting a cline are too unfavorable because migration is too strong, selection is too weak, or both, a cline cannot exist at all. Thus, unless there is overdominance, the population must be monomorphic. It is possible for a cline which can barely exist under the prevailing ecological circumstances to show a large amount of variation in gene frequency.     

MATING PATTERN AND FITNESS-COMPONENT ANALYSIS ASSOCIATED WITH INVERSION POLYMORPHISM IN A NATURAL POPULATION OF DROSOPHILA BUZZATII

Direct studies of mating success or mating pattern associated with Mendelian factors rarely have been carried out in nature. From the samples taken for the standard analyses of selection components, it is not usually possible to obtain the mating table, and only directional selection for male mating success can be detected. Both processes, mating pattern and differential mating probability, together with other fitness components, have been investigated for the inversion polymorphism of a natural population of the cactophilic species Drosophila buzzatii. Two independent samples of adult flies were collected: nonmating or single individuals (base population) and mating pairs (mating population). All individuals were karyotyped for the second and fourth chromosomes. A sequence of models with increasing simplicity was fitted to the data to test null hypotheses of no selection and random union of gametes and karyotypes. The main results were (1) no deviations from random mating were found; (2) differential mating probability was nonsignificant in both sexes; (3) inversion and karyotypic frequencies did not differ between sexes; and (4) karyotypic frequencies did not depart from Hardy-Weinberg expectations. These results are discussed in light of complementary evidence showing the need for interpreting with caution no-effect hypotheses such as the ones tested here. The use of complementary selective tests in these studies is suggested.     

Geographical invariance in population genetics

Various quantities of evolutionary interest are shown to be independent of geographical structure. A diploid, monoecious population is subdivided into a finite number of panmictic colonies that exchange migrants. The migration pattern is fixed and ergodic, but otherwise arbitrary. Generations are discrete and non-overlapping; the analysis is restricted to a single locus. Previous results are generalized in the neutral multiallelic case. With selection, it is assumed that there are only two alleles, dominance is absent, selection has the same intensity in all demes, migration does not change the subpopulation numbers, and all evolutionary forces are weak. A diffusion approximation is established for the gene frequencies, and the invariance of the fixation probability and of the moments of the conditional and unconditional total heterozygosities before absorption is demonstrated by a martingale argument.     

ASSORTATIVE MATE CHOICE AND DOMINANCE MODIFICATION: ALTERNATIVE WAYS OF REMOVING HETEROZYGOTE DISADVANTAGE

In genetic polymorphisms of two alleles, heterozygous individuals may contribute to the next generation on average more or fewer descendants than the homozygotes. Two different evolutionary responses that remove a disadvantageous heterozygote phenotype from the population are the evolution of strictly assortative mate choice, and that of a modifier making one of the two alleles completely dominant. We derive invasion fitness of mutants introducing dominance or assortative mate choice in a randomly mating population with a genetic polymorphism for an ecological trait. Mutations with small effects as well as mutants introducing complete dominance or perfect assorting are considered. Using adaptive dynamics techniques, we are able to calculate the ratio of fitness gradients for the effects of a dominance modifier and a mate choice locus, near evolutionary branching points. With equal resident allele frequencies, selection for mate choice is always stronger. Dominance is more strongly selected than assortative mating when the resident (common) alleles have very unequal frequencies at equilibrium. With female mate choice the difference in frequencies where dominance is more strongly selected is smaller than when mutants of both sexes can choose without costs. A symmetric resource-competition model illustrates the results.     

The robustness of neutral models of geographical variation

The approximation of diploid migration by gametic dispersion is studied. The monoecious, diploid population is subdivided into panmictic colonies that exchange migrants. Generations are discrete and nonoverlapping; the analysis is restricted to a single locus in the absence of selection; every allele mutates to a new allele at the same rate u. Diploid-migration models without self-fertilization and with selfing at the “random” rate (equal to the reciprocal of the deme size in each deme) are investigated; in the gametic-dispersion models, selfing occurs at the random rate. It is shown for the unbounded stepping-stone model in one and two dimensions, the circular stepping-stone model, and the island model that the probabilitities of identity in state at equilibrium for diploid migration are close to those for gametic dispersion if the mutation rate is small or the deme size is large. Explicit error bounds are presented in all the above cases. It is also proved that if the number of demes is finite and the migration matrix is arbitrary but time independent and ergodic, then in the strong-migration approximation the equilibrium and the ultimate rate and pattern of convergence of both diploid-dispersion models are close to the corresponding gametic-dispersion formulae. For the strong-migration approximation at equilibrium, migration must dominate both mutation and random drift; for the convergence results, it suffices that migration dominate random drift. All the results apply to a dioecious population if the migration pattern and mutation rate are sex independent.     

Approaches to the Hardy-Weinberg manifold

Let fertilities and death rates be additive, let fertilities be positive, and let mating be random in the Nagylaki-Crow continuous model of selection at a multiallelic locus in a monoecious population. Then polymorphisms are in Hardy-Weinberg proportions. If some fertilities vanish, there is an example of a diallelic polymorphism that is not in Hardy-Weinberg proportions. If the fertilities are larger, in one sense or another, than the difference between any two death rates, then convergence to the Hardy-Weinberg manifold is shown. If, in addition, the Malthusian parameters are constant, and only a finite number of equilibria exist, then global convergence to equilibria is proved.     

The evolution of intratetrad mating rates

Intratetrad mating, the fusion of gametes formed in a single meiosis, has unusual consequences for genetic diversity, especially in genome regions linked to mating type loci. Here we investigate the fate of modifier alleles that alter the rate of intratetrad mating, under models of heterozygote advantage and of genetic load resulting from recurrent mutation. In both cases, intratetrad mating is favored if the recombination rate between the selected locus and mating type is less than the frequency of lethal recessive alleles at that locus in the population. Positive feedback often accelerates the invasion of modifiers to the intratetrad mating rate. Recombination rate and intratetrad mating rate exert indirect selection on one another, resulting in a cascading decline in outcrossing, even in the absence of any cost of sex. However, under recurrent mutation, alleles for obligate intratetrad mating invade only very slowly, perhaps explaining why outcrossing can persist at low frequencies in a largely intratetrad mating population.     

Polymorphism with selection and genotype-dependent migration

For a single autosomal locus with multiple alleles both an island and a multiple-niche model with discrete nonoverlapping generations are formulated for the maintenance of genetic variability. Both models incorporate viability selection in an arbitrary way and allow for genotypic differences in the pertinent migration structure. Random drift is ignored, and mating is at random. A global analysis is given for the island model in the neutral case. For a subdivided population, conditions are derived for the existence of a protected polymorphism, and the model is examined in some special two-niche cases. Of particular consideration is the loss of neutral alleles due solely to population regulation and genotype-dependent migration, and the possible existence of equilibrium clines without selection.M. M. was supported by USPHS Pre-doctoral training grant No. GM 7197 to the University of Chicago; this work represents part of the author's Doctoral dissertation.     

Population trajectories for single locus additive fecundity selection and related selection models

A selection model which comprises models of additive fecundities as well as models of viability, fecundity, or differential mating selection acting only in one sex, is investigated for an autosomal gene locus in a population reproducing in nonoverlapping generations. The recurrence equations and basic properties of the genotypic population trajectories and equilibrium points are formulated for the multiallelic case. For the diallelic case, the trajectory development is discussed in more detail, and it is proven that every population trajectory converges to a Hardy-Weinberg equilibrium point.     

Hardy-Weinberg equilibria in random mating populations

The structure of multiloci random mating populations is examined. Sufficient conditions for the existence of stable local Hardy-Weinberg equilibria for n loci and an arbitrary number of alleles per locus, are then derived for specified situations under the assumption of multiplicative gene action between loci. It is shown that a stable Hardy-Weinberg equilibrium can not be a local maximum of the mean fitness function with multiplicative gene action between loci. The stability of Hardy-Weinberg type border points and the condition for the increase of newly introduced genes are topics on which some n-loci results are also obtained for an arbitrary number of alleles per locus in systems that allow Hardy-Weinberg equilibria.