Density dependent selection 1: A stable feasible equilibrium may not be attainable

The joint evolution of gene frequency p, and population size N is studied. It is shown that when the genotypic fitnesses are logistic functions of the population size, sets of initial states exist which lead to bizarre behavior. Even though equilibria may be locally stable, these sets of initial conditions eventually produce negative fitnesses. Alternative models are discussed as are some general properties of the mean fitness.     

Stable Equilibria at Two Loci in Populations with Large Selfing Rates

The equilibrium structure of two-locus, two-allele models with very large selfing rates is found using perturbation techniques. For free recombination, r = 1/2, the following results hold. If the heterozygotes do not have at least an approximate 30% advantage in fitness relative to homozygotes, a stable equilibrium with all alleles present is possible only if all of the homozygote fitnesses differ at most by approximately the outcrossing rate, t, and all stable polymorphic equilibria have disequilibrium values, D, that are at most on the order of the outcrossing rate. Once the heterozygote fitnesses are above the threshold, there are stable equilibria possible with D near its maximum possible value. The results show that the observed disequilibria in highly selfed plant populations are not likely to result from selection leading to an equilibrium.     

The Alternative Fitness Sets Which Preserve Allele Trajectories: A General Treatment

A general solution is presented of the problem of specifying all alternative, generally frequency-dependent, (absolute) fitness sets which give rise to the same allele frequency changes and population dynamics as a given fitness set. The one- and two-locus cases are analyzed in detail and the method is then extended to the n-locus case. It is shown that if biological constraints can be used to specify the mean fitness of the population and the relative fitnesses of the heterozygotes, then the allele frequency trajectories determine a unique fitness set.     

Protected polymorphism in the two-locus haploid model with unpredictable fitnesses

 In an unpredictable environment, the distributions of alleles from which polymorphism can be maintained forever belong to a certain set, the C-viability kernel. Such a set is calculated in the two-locus haploid model, as well as the corresponding fitnesses at any time which make this maintenance possible. The dependence of the C-viability kernel on the set U of admissible fitnesses and on the recombination rate r is studied. Notably, the C-viability kernel varies rapidly in the neighborhood of equal fitness of AB and ab; it becomes empty when ab has a fitness below a certain function, which is delineated, of the recombination rate. The properties of the two-locus model under constraints, out of equilibrium and with unpredictable selection are thus presented. Received: 20 May 1999     

Establishment of allelic two-locus polymorphisms

An exact analysis of necessary and sufficient conditions for the establishment and protectedness of biallelic two-locus polymorphisms is developed for the classical model with constant, sexually symmetric fitnesses and random association of the successful gametes. To demonstrate application of the results to common model types, the model of symmetric viabilities depending on the degree of heterozygosity only is chosen as a paradigm. It is pointed out that a unique locally stable internal equilibrium may exist even though all marginal equilibria (including the fixation states) are locally attractive. This example is quoted as an indication of the priority that analyses of protectedness deserve over analyses of local stability or instability of internal equilibria. Further applications of broader appeal concern the role that recombination plays in protecting polymorphisms. Probably the most interesting finding is that with increasing recombination frequency the chances for protectedness of a polymorphism generally decline. Yet, if a certain hierarchic ordering of the fitnesses with respect to the degree of heterozygosity is realized, the polymorphism is protected for arbitrary amounts of recombination. If recombination is rare, heterozygote advantage is not a universal precondition for persistence of polymorphisms. This phenomenon is utilized to derive conditions under which deleterious recessive mutants can be maintained in a population.     

Selection in heterogeneous environments maintains the gene arrangement polymorphism of Drosophila pseudoobscura

Chromosomal rearrangements may play an important role in how populations adapt to a local environment. The gene arrangement polymorphism on the third chromosome of Drosophila pseudoobscura is a model system to help determine the role that inversions play in the evolution of this species. The gene arrangements are the likely target of strong selection because they form classical clines across diverse geographic habitats, they cycle in frequency over seasons, and they form stable equilibria in population cages. A numerical approach was developed to estimate the fitness sets for 15 gene arrangement karyotypes in six niches based on a model of selection-migration balance. Gene arrangement frequencies in the six different niches were able to reach a stable meta-population equilibrium that matched the observed gene arrangement frequencies when recursions used the estimated fitnesses with a variety of initial inversion frequencies. These analyses show that a complex pattern of selection is operating in the six niches to maintain the D. pseudoobscura gene arrangement polymorphism. Models of local adaptation predict that the new inversion mutations were able to invade populations because they held combinations of two to 13 local adaptation loci together.     

On the number of stable equilibria and the simultaneous stability of fixation and polymorphism in two-locus models

It is shown that in simple symmetric two-locus, two-allele constant fitness models the bound of four simultaneously stable equilibria previously accepted for general two-locus, two-allele models is exceeded. Situations with five and six stable equilibria are exhibited. These involve four chromosomal fixations and either one or two polymorphic stable equilibria.     

Examples of the effect of genetic variation on competing species

An ordinary differential equation model for two competing populations with genetic variation in one population is presented. The degree of frequency dependence needed to produce various configurations of stable equilibria is discussed. For example, if the fitnesses are frequency independent then there may exist stable polymorphism although the genetically varying population becomes extinct in each fixation plane. Stable polymorphism where the genetically invariant population becomes extinct in each fixation plane requires frequency dependence in the fitness of the genetically invariant population.     

Polymorphism and differential selection for the sexes.

Various genetic models with different fitnesses for the sexes are investigated. Only a limited set of fitness values will result in a stable polymorphism, and the rate of approach to these equilibrium frequencies is extremely slow. These results indicate that there are problems associated with the interpretation of some human genetic traits by such models.     

Disequilibrium in Two-Locus Mutation-Selection Balance Models

Equilibrium behavior of two-locus mutation-selection balance models is analyzed using perturbation techniques. The classical result of Haldane for one locus is shown to carry over to two loci, if fitnesses are replaced by marginal fitnesses. If the fitness of the double heterozygote is smaller than would be produced by a multiplicative model, as in additive or quantitative fitness models, the disequilibrium is negative--an excess of gametes with one rare allele. In this case the disequilibrium can be as large as one-half its maximum value possible, if the recombination rate is small, not greater than the strength of selection. If the fitness of the double heterozygote is larger than would be produced by a multiplicative model, the disequilibrium is positive, and is very small relative to its maximum value possible, even if the recombination rate is zero.     

The maintenance (or not) of polygenic variation by soft selection in heterogeneous environments

On the basis of single-locus models, spatial heterogeneity of the environment coupled with strong population regulation within each habitat (soft selection) is considered an important mechanism maintaining genetic variation. We studied the capacity of soft selection to maintain polygenic variation for a trait determined by several additive loci, selected in opposite directions in two habitats connected by dispersal. We found three main types of stable equilibria. Extreme equilibria are characterized by extreme specialization to one habitat and loss of polymorphism. They are analogous to monomorphic equilibria in singe-locus models and are favored by similar factors: high dispersal, weak selection, and low marginal average fitness of intermediate genotypes. At the remaining two types of equilibria the population mean is intermediate but variance is very different. At fully polymorphic equilibria all loci are polymorphic, whereas at low-variance equilibria at most one locus remains polymorphic. For most parameters only one type of equilibrium is stable; the transition between the domains of fully polymorphic and low-variance equilibria is typically sharp. Low-variance equilibria are favored by high marginal average fitness of intermediate genotypes, in contrast to single-locus models, in which marginal overdominance is particularly favorable for maintenance of polymorphism. The capacity of soft selection to maintain polygenic variation is thus more limited than extrapolation from single-locus models would suggest, in particular if dispersal is high and selection weak. This is because in a polygenic model, variance can evolve independently of the mean, whereas in the single-locus two-allele case, selection for an intermediate mean automatically leads to maintenance of polymorphism.     

A two-locus model for polymorphism for sex-linked meiotic drive modifiers with possible applications to Aedes aegypti

A two-locus model is presented which shows the possibility of maintaining a polymorphism for modifiers of sex-linked meiotic drive in the absence of fitness differences. The model is very similar to the situation actually found in some laboratory strains of the mosquito Aedes aegypti. The existence of a stable polymorphism usually requires sufficiently loose linkage between the two loci.     

The significance of over- and underdominance for the maintenance of genetic polymorphisms. I. Underdominance and stability

It is pointed out that the standard selection models in population genetics all require some form of heterozygote advantage in fitness in order to guarantee the maintenance or stability of genetic polymorphisms. Even more recent results demonstrating the existence of stable two-locus polymorphisms with marginal underdominance at both loci are based on certain epistatically acting heterosis assumptions. This raises the question as to whether heterozygote advantage in fitness is indeed a generally valid principle of maintaining polymorphisms. To avoid ambiguity in definition of heterozygote advantage (overdominance) as it appears in multiallele or multilocus systems, a one-locus-two-allele model is considered. This model allows for sexually asymmetric selection and random mating. It is shown that the model produces globally stable polymorphisms exhibiting underdominance in fitness for a considerable and biologically reasonable range of selection values. Having thus properly refuted the general validity of the common overdominance principle, a modified version is suggested which covers the classical viability selection model and its extension to arbitrary, sexually asymmetric viability and fertility selection. This modified overdominance principle is based on the notion of fractional fitnesses and relates protectedness of biallelic polymorphisms to the extent to which each genotype reproduces its own type. The fact that the model treated displays frequency dependent fitnesses which may change in ranking while approaching equilibrium is discussed in relation to problems of the evolution of overdominance and underdominance.     

The significance of over- and underdominance for the maintenance of genetic polymorphisms. II. Overdominance and instability with random mating

Part I of the present series demonstrates that globally stable polymorphic equilibria may show underdominance in Darwinian fitness. Hence, overdominance in fitness can no longer be conceived of as a necessary condition for the stability of a polymorphism. In the present paper, the question is posed as to whether overdominance is at least sufficient for this stability. A population of randomly mating individuals is considered, where selection operates uniquely through differential fecundities of particular mating types and may generate either a heterozygote excess or deficit relative to Hardy-Weinberg proportions. It turns out that both unstable central overdominance and stable central underdominance are possible and that their occurrence is strongly related to an excess or a deficiency of heterozygotes in the vicinity of the regions of instability or stability. As one consequence, the above suggested sufficiency of heterozygote superiority is not valid, even in random mating populations. Based on the results of both papers of this series, which demonstrate the inadequacy of over- and underdominance as indicators of stability or instability, a modified overdominance principle is discussed. This principle states that a biallelic polymorphism is maintained if the heterozygote is superior in its degree of "heterogamous self-replication" to the degrees of "autogamous self-replication" of the corresponding homozygotes. It is derived with the help of fractional fitnesses, and it is pointed out that certain ratios of these may be more useful for finding evolutionary constants which govern the maintenance of genetic polymorphisms than are ratios of total fitnesses.     

Evolutionary principles for general frequency-dependent two-phenotype models in sexual populations

The evolutionary dynamics in general two-sex two-phenotype frequency-dependent selection models are studied with respect to underlying multi-allele one-locus genetic systems. Two classes of equilibria come into play: genotypic equilibria, with equilibrium allelic frequencies independent of the phenotype, and phenotypic equilibria, which are characterized by equal mean phenotypic fitnesses. The exact conditions for genotypic equilibria to exist and be stable and for phenotypic equilibria to exist and be evolutionarily attractive are examined. Using adequate definitions of mean fitnesses in general contexts of frequency-dependent selection in dioecious populations, we show that two phenotypes, when they can coexist in the population, tend to balance their fitnesses as far as is allowed by the genetic system as more alleles responsible for phenotype determination are introduced into the population.     

Stable genetic polymorphism in heterogeneous environments: balance between asymmetrical dispersal and selection in the acorn barnacle

Elucidating the processes responsible for maintaining polymorphism at ecologically relevant genes is intimately related to understanding the interplay between selection imposed by habitat heterogeneity and a species' capacity for dispersal in the face of environmental constraints. In this paper, we used a model-based approach to solve equilibria of balanced polymorphism, given values of fitness and larval dispersal among different habitats in the acorn barnacle Semibalanus balanoides from the Gulf of St Lawrence. Our results showed that allele frequencies observed at both MPI* and GPI* loci represented stable equilibria, given empirical estimates of fitness values, and that considerably more larvae dispersed from one region (north) to the other (south) than vice versa. Dispersal conditions were predicted to be similar for the maintenance of polymorphism at both loci. Moreover, the values of asymmetrical dispersal required by the model to reach stable equilibria were compatible with empirical estimates of larval dispersal and oceanic circulation documented in this system. Overall, this study illustrated the usefulness of a modified and computable version of Bulmer's model (1972) in order to test hypotheses of balanced polymorphism resulting from interactions between spatial selection and asymmetrical dispersal.     

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.     

On the mutual maintenance of a polymorphism in two nonepistatic loci with partial selfing

In a two-locus system of multiplicative fitnesses with a mixed selfing and random mating, the maintenance of a polymorphism in each locus is not independent of its maintenance in the other one.     

CYCLICAL ENVIRONMENTAL CHANGES AS A FACTOR MAINTAINING GENETIC POLYMORPHISM. 2. DIPLOID SELECTION FOR AN ADDITIVE TRAIT

The subject of this paper is polymorphism maintenance due to stabilizing selection with a moving optimum. It was shown that in case of two-locus additive control of the selected trait, global polymorphism is possible only when the geometric mean fitnesses of double homozygotes averaged over the period are lower than that of the single heterozygotes and of the double heterozygote (with a multiplier [1 – r]^p, which depends on recombination rate r and period length p). But local stability of polymorphism cannot be excluded even if geometric mean fitnesses of all double homozygotes are higher than that of all heterozygotes. We proved, that for logarithmically convex fitness functions, cyclical changes of the optimum cannot help in polymorphism maintenance in case of additive control of the selected trait by two equal loci. However, within the same class of fitness functions, nonequal gene action and/or dominance effect for one or both loci may lead to local polymorphism stability with large enough polymorphism attracting domain. The higher the intensity of selection and closer the linkage between selected loci the larger is this domain. Note that even simple cyclical selection could result in two forms of polymorphic limiting behavior: (a) usually expected forced cycle with a period equal to that of environmental changes; and (b) “supercycles,” nondumping auto-oscillations with a period comprising of hundreds of forced oscillation periods.     

Frequency-Dependent Selection: The High Potential for Permanent Genetic Variation in the Diallelic, Pairwise Interaction Model

A detailed analytic and numerical study is made of the potential for permanent genetic variation in frequency-dependent models based on pairwise interactions among genotypes at a single diallelic locus. The full equilibrium structure and qualitative gene-frequency dynamics are derived analytically for a symmetric model, in which pairwise fitnesses are chiefly determined by the genetic similarity of the individuals involved. This is supplemented by an extensive numerical investigation of the general model, the symmetric model, and nine other special cases. Together the results show that there is a high potential for permanent genetic diversity in the pairwise interaction model, and provide insight into the extent to which various forms of genotypic interactions enhance or reduce this potential. Technically, although two stable polymorphic equilibria are possible, the increased likelihood of maintaining both alleles, and the poor performance of protected polymorphism conditions as a measure of this likelihood, are primarily due to a greater variety and frequency of equilibrium patterns with one stable polymorphic equilibrium, in conjunction with a disproportionately large domain of attraction for stable internal equilibria.