Equilibrium properties of a multi-locus, haploid-selection, symmetric-viability model

Under haploid selection, a multi-locus, diallelic, two-niche Levene (1953) model is studied. Viability coefficients with symmetrically opposing directional selection in each niche are assumed, and with a further simplification that the most and least favored haplotype in each niche shares no alleles in common, and that the selection coefficients monotonically increase or decrease with the number of alleles shared. This model always admits a fully polymorphic symmetric equilibrium, which may or may not be stable.We show that a stable symmetric equilibrium can become unstable via either a supercritical or subcritical pitchfork bifurcation. In the supercritical bifurcation, the symmetric equilibrium bifurcates to a pair of stable fully polymorphic asymmetric equilibria; in the subcritical bifurcation, the symmetric equilibrium bifurcates to a pair of unstable fully polymorphic asymmetric equilibria, which then connect to either another pair of stable fully polymorphic asymmetric equilibria through saddle-node bifurcations, or to a pair of monomorphic equilibria through transcritical bifurcations. As many as three fully polymorphic stable equilibria can coexist, and jump bifurcations can occur between these equilibria when model parameters are varied.In our Levene model, increasing recombination can act to either increase or decrease the genetic diversity of a population. By generating more hybrid offspring from the mating of purebreds, recombination can act to increase genetic diversity provided the symmetric equilibrium remains stable. But by destabilizing the symmetric equilibrium, recombination can ultimately act to decrease genetic diversity.     

Multilocus genetics and the coevolution of quantitative traits

We develop and analyze an explicit multilocus genetic model of coevolution. We assume that interactions between two species (mutualists, competitors, or victim and exploiter) are mediated by a pair of additive quantitative traits that are also subject to direct stabilizing selection toward intermediate optima. Using a weak-selection approximation, we derive analytical results for a symmetric case with equal locus effects and no mutation, and we complement these results by numerical simulations of more general cases. We show that mutualistic and competitive interactions always result in coevolution toward a stable equilibrium with no more than one polymorphic locus per species. Victim-exploiter interactions can lead to different dynamic regimes including evolution toward stable equilibria, cycles, and chaos. At equilibrium, the victim is often characterized by a very large genetic variance, whereas the exploiter is polymorphic in no more than one locus. Compared to related one-locus or quantitative genetic models, the multilocus model exhibits two major new properties. First, the equilibrium structure is considerably more complex. We derive detailed conditions for the existence and stability of various classes of equilibria and demonstrate the possibility of multiple simultaneously stable states. Second, the genetic variances change dynamically, which in turn significantly affects the dynamics of the mean trait values. In particular, the dynamics tend to be destabilized by an increase in the number of loci.     

Evolutionary Stability for Interactions among Kin under Quantitative Inheritance

The theory of evolutionarily stable strategies (ESS) predicts the long-term evolutionary outcome of frequency-dependent selection by making a number of simplifying assumptions about the genetic basis of inheritance. I use a symmetrized multilocus model of quantitative inheritance without mutation to analyze the results of interactions between pairs of related individuals and compare the equilibria to those found by ESS analysis. It is assumed that the fitness changes due to interactions can be approximated by the exponential of a quadratic surface. The major results are the following. (1) The evolutionarily stable phenotypes found by ESS analysis are always equilibria of the model studied here. (2) When relatives interact, one of the two conditions for stability of equilibria differs between the two models; this can be accounted for by positing that the inclusive fitness function for quantitative characters is slightly different from the inclusive fitness function for characters determined by a single locus. (3) The inclusion of environmental variance will in general change the equilibrium phenotype, but the equilibria of ESS analysis are changed to the same extent by environmental variance. (4) A class of genetically polymorphic equilibria occur, which in the present model are always unstable. These results expand the range of conditions under which one can validly predict the evolution of pairwise interactions using ESS analysis.     

Selection in Complex Genetic Systems I. the Symmetric Equilibria of the Three-Locus Symmetric Viability Model

The symmetric equilibria of the three-locus symmetric viability model are determined and their stability analyzed. For tight linkage there may be four stable equilibria, each characterized by having one pair of complementary chromosomes in high frequencies, with all others low. For looser linkage the only stable symmetric equilibrium is that with complete linkage equilibrium. For intermediate recombination values both types of equilibria may be stable. A new class of equilibria with all pairwise linkage disequilibria zero, but with third order linkage disequilibrium, has been discovered. It may be stable for tight linkage.     

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.     

Sex-specific meiotic drive and selection at an imprinted locus

We present a one-locus model that breaks two symmetries of Mendelian genetics. Whereas symmetry of transmission is breached by allowing sex-specific segregation distortion, symmetry of expression is breached by allowing genomic imprinting. Simple conditions for the existence of at least one polymorphic stable equilibrium are provided. In general, population mean fitness is not maximized at polymorphic equilibria. However, mean fitness at a polymorphic equilibrium with segregation distortion may be higher than mean fitness at the corresponding equilibrium with Mendelian segregation if one (or both) of the heterozygote classes has higher fitness than both homozygote classes. In this case, mean fitness is maximized by complete, but opposite, drive in the two sexes. We undertook an extensive numerical analysis of the parameter space, finding, for the first time in this class of models, parameter sets yielding two stable polymorphic equilibria. Multiple equilibria exist both with and without genomic imprinting, although they occurred in a greater proportion of parameter sets with genomic imprinting.     

The two-locus model of Gaussian stabilizing selection

We study the equilibrium structure of a well-known two-locus model in which two diallelic loci contribute additively to a quantitative trait that is under Gaussian stabilizing selection. The population is assumed to be infinitely large, randomly mating, and having discrete generations. The two loci may have arbitrary effects on the trait, the strength of selection and the recombination rate may also be arbitrary. We find that 16 different equilibrium patterns exist, having up to 11 equilibria; up to seven interior equilibria may coexist, and up to four interior equilibria, three in negative and one in positive linkage disequilibrium, may be simultaneously stable. Also, two monomorphic and two fully polymorphic equilibria may be simultaneously stable. Therefore, the result of evolution may be highly sensitive to perturbations in the initial conditions or in the underlying genetic parameters. For the special case of equal effects, global stability results are proved. In the general case, we rely in part on numerical computations. The results are compared with previous analyses of the special case of extremely strong selection, of an approximate model that assumes linkage equilibrium, and of the much simpler quadratic optimum model.     

Polymorphic equilibria with assortative mating and selection in subdivided populations

Two modes of assortative mating, partial selfing and assorting by phenotypic classes, are studied in a subdivided population. Differential viability is allowed and the selection intensities and assorting tendencies are permitted to vary among the habitats. There exists a symmetric polymorphism; we delimit its level of heterozygosity and stability nature (dependent on the selection intensities and assorting propensities). This complements studies of the fixation states and thereby provides further insight into the global equilibrium structure in subdivided populations. Circumstances are given where the fixation states and symmetric polymorphism comprise the global equilibrium structure. Examples are also given where migration engenders stable polymorphic equilibria and stable polymorphic equilibrium cycles which are absent in single demes without migration.     

Maintenance of Genetic Variability under Strong Stabilizing Selection: A Two-Locus Model

We study a two locus model with additive contributions to the phenotype to explore the relationship between stabilizing selection and recombination. We show that if the double heterozygote has the optimum phenotype and the contributions of the loci to the trait are different, then any symmetric stabilizing selection fitness function can maintain genetic variability provided selection is sufficiently strong relative to linkage. We present results of a detailed analysis of the quadratic fitness function which show that selection need not be extremely strong relative to recombination for the polymorphic equilibria to be stable. At these polymorphic equilibria the mean value of the trait, in general, is not equal to the optimum phenotype, there exists a large level of negative linkage disequilibrium which ``hides' additive genetic variance, and different equilibria can be stable simultaneously. We analyze dependence of different characteristics of these equilibria on the location of optimum phenotype, on the difference in allelic effect, and on the strength of selection relative to recombination. Our overall result that stabilizing selection does not necessarily eliminate genetic variability is compatible with some experimental results where the lines subject to strong stabilizing selection did not have significant reductions in genetic variability.     

Population genetics of modifiers of meiotic drive III. Equilibrium analysis of a general model for the genetic control of segregation distortion

Prout, Bungaard and Bryant (1973, Theor. Popul. Biol. 4, 446–465) presented the first formal treatment of a model of meiotic drive involving a modifier locus which controls the intensity of drive. They studied the equilibrium behavior in the simplest model where it is assumed that drive is maximal when not suppressed. In that case there is one polymorphic equilibrium at which there is linkage disequilibrium. The equilibrium solutions in the general model of meiotic drive proposed by Prout, et al. are given in this paper together with a stability analysis. It is shown that up to three polymorphic equilibria may exist, two of which are in linkage disequilibrium and one in linkage equilibrium. These equilibria exhibit behavior qualitatively opposite to what is widely accepted as the usual for two locus systems and which is not seem in the simple case originally treated. The polymorphic equilibria with linkage disequilibrium may be stable for loose linkage and not for tight while that with linkage equilibrium is stable in an interval of relatively tight linkage values.     

Equilibrium structure and stability in a frequency-dependent,two-population diploid model

We investigate the equilibrium structure for an evolutionary genetic model in discrete time involving two monoecious populations subject to intraspecific and interspecific random pairwise interactions. A characterization for local stability of an equilibrium is found, related to the proximity of this equilibrium with evolutionarily stable strategies (ESS). This extends to a multi-population framework a principle initially proposed for single populations, which states that the mean population strategy at a locally stable equilibrium is as close as possible to an ESS.     

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.     

The Maintenance of Genetic Variability in Two-Locus Models of Stabilizing Selection

The maintenance of genetic variability at two diallelic loci under stabilizing selection is investigated. Generations are discrete and nonoverlapping; mating is random; mutation and random genetic drift are absent; selection operates only through viability differences. The determination of the genotypic values is purely additive. The fitness function has its optimum at the value of the double heterozygote and decreases monotonically and symmetrically from its optimum, but is otherwise arbitrary. The resulting fitness scheme is identical to the symmetric viability model. Linkage disequilibrium is neglected, but the results are otherwise exact. Explicit formulas are found for all the equilibria, and explicit conditions are derived fro their existence and stability. A complete classification of the six possible global convergence patterns is presented. In addition to the symmetric equilibrium (with gene frequency 1/2 at both loci), a pair of unsymmetric equilibria may exist; the latter are usually, but not always, unstable. If the ratio of the effect of the major locus to that of the minor one exceeds a critical value, both loci will be stably polymorphic. If selection is weak at the minor locus, the more rapidly the fitness function decreases near the optimum, the lower is this critical value; for rapidly decreasing fitness functions, the critical value is close to one. If the fitness function is smooth at the optimum, then a stable polymorphism exists at both loci only if selection is strong at the major locus.     

A two-locus mutation—Selection model and some of its evolutionary implications

The deterministic properties of a two-locus model with mutation and selection have been investigated. The mutation process is unidirectional, and the model is so constructed that the genetic variation at one locus is selectively neutral in the absence of a mutant allele at the other locus. All genotypes with three or four mutant alleles are deleterious, while the double heterozygotes may have the same fitness as the standard genotype. If one of the mutant alleles becomes fixed in the population, then the other locus will show a regular one-locus mutation-selection balance. Such a boundary equilibrium may be unstable or stable in the full two-locus setting. In the symmetric case, which is analyzed in details, the population will either go to one of the two boundary equilibria, or to a fully polymorphic equilibrium at which both the mutant alleles are rare. The origin of reproductive separation between two populations via the fixation of complementary deleterious mutants at different loci, and the fixation of nonfunctional alleles at duplicated loci, are two biological processes which both can be studied with the present model. In the last part of the paper we show how the results from the deterministic analysis can be used to predict how different factors will influence the rates of evolution in these systems.     

On some models of fertility selection

Additive, multiplicative and symmetric models of fertility controlled by one diallelic gene are studied. For the completely symmetric fertility system a complete equilibrium and local stability analysis is possible. Contrary to previous conjectures, asymmetric equilibria can be stable. Conditions are derived under which a multiplicative model can be regarded as equivalent to a symmetric fertility system.     

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.     

Four simultaneously stable polymorphic equilibria in two-locus two-allele models

The existence of four simultaneously stable equilibria with both loci polymorphic is shown for the Lewontin-Kojima version of the two-locus two-allele symmetric viability model, using bifurcation theory. This exceeds the previously claimed bound of two stable polymorphisms. Biological implications of the result are discussed.     

A semi-symmetric two-locus model

The two-locus symmetric viability model characterized by its invariance with respect to the exchange of alleles at each locus, is a well-studied model of classical two-locus theory. The symmetric model introduced by Lewontin and Kojima is among the few multi-locus models with epistatic interactions between loci for which a polymorphism with linkage equilibrium can be stable and this happens when recombination is sufficiently large. We show that an analogous property holds true for a different model, in which symmetry need exist at only one locus. The properties of this new semi-symmetric model are compared with those of the classical symmetric model. For tight linkage, two classes of polymorphisms are possible, depending on the magnitude of additive epistasis. The recombination rate above which linkage equilibrium becomes stable is derived analytically. As in the symmetric model, intervals of recombination in which no polymorphism is stable are possible, and stable polymorphisms can coexist with stable fixations.     

Dominance and the maintenance of polymorphism in multiallelic migration-selection models with two demes

The maintenance of genetic variation in a spatially heterogeneous environment has been one of the main research themes in theoretical population genetics. Despite considerable progress in understanding the consequences of spatially structured environments on genetic variation, many problems remain unsolved. One of them concerns the relationship between the number of demes, the degree of dominance, and the maximum number of alleles that can be maintained by selection in a subdivided population. In this work, we study the potential of maintaining genetic variation in a two-deme model with deme-independent degree of intermediate dominance, which includes absence of G×E interaction as a special case. We present a thorough numerical analysis of a two-deme three-allele model, which allows us to identify dominance and selection patterns that harbor the potential for stable triallelic equilibria. The information gained by this approach is then used to construct an example in which existence and asymptotic stability of a fully polymorphic equilibrium can be proved analytically. Noteworthy, in this example the parameter range in which three alleles can coexist is maximized for intermediate migration rates. Our results can be interpreted in a specialist-generalist context and (among others) show when two specialists can coexist with a generalist in two demes if the degree of dominance is deme independent and intermediate. The dominance relation between the generalist allele and the specialist alleles play a decisive role. We also discuss linear selection on a quantitative trait and show that G×E interaction is not necessary for the maintenance of more than two alleles in two demes.