Central equilibria in multilocus systems. I. Generalized nonepistatic selection regimes

The generalized nonepistatic selection regime encompasses combinations of multiplicative and neutral viability effects distributed across a set of loci. These subsume, in particular, mixtures of the classical modes of multiplicative and additive fitness evaluations for multilocus traits. Exact analytic conditions for existence and stability of a multilocus Hardy-Weinberg (H-W) polymorphic equilibrium configuration are ascertained. It is established that the central H-W polymorphism is stable only if the component loci are "over-dominant" and sufficient recombination is in force. The H-W central equilibrium is never stable for tight linkage whenever some multiplicative selection effects are contributed by at least two of the loci involved. In the case of additive selection expression and individual overdominant loci, the H-W polymorphism is stable independently of the level of recombination. In the context of "natural" recombination schemes, "more recombination" enhances the stability of the H-W polymorphic equilibrium.     

Central Equilibria in Multilocus Systems. II. Bisexual Generalized Nonepistatic Selection Models

This paper is a continuation of the paper "Central Equilibria in Multilocus Systems I," concentrating on existence and stability properties accruing to central H-W type equilibria in multilocus bisexual systems acted on by generalized nonepistatic selection forces coupled to recombination events. The stability conditions are discussed and interpreted in three perspectives, and the influence of sexual differences in linkage relationships together with sex-dependent selection is appraised. In this case we deduce that the stability conditions of the H-W polymorphism in the bisexual model coincide exactly with the conditions for the corresponding monoecious model, provided that the recombination distribution imposed is that of the arithmetic mean of the male and female recombination distributions. A second concern has the same recombination distribution for both sexes, but contrasting selection regimes between sexes. It is then established that, with respect to discerning the relevance of the H-W equilibrium, there is an equivalent monoecious selection regime which is an appropriate "weighted combination" of the male and female selection forms. Finally, in the case where the selection and recombination structures are both sex dependent, a hierarchy of comparisons is elaborated, seeking to unravel the nature of selection-recombination interaction for monoecious versus diocecious systems.     

Maintenance of multilocus variability under strong stabilizing selection

We present exact conditions for stability of monomorphic equilibria in a general multilocus multiallele system and of specific polymorphic equilibria in general one- and two-locus multiallele systems. We show how these exact results on one- and two-locus systems can be used in approximate analysis of polymorphic equilibria in multilocus systems under selection strong relative to recombination. We determine conditions for existence and stability of polymorphic equilibria in specific models of quadratic stabilizing selection on additive polygenic traits.     

The generalized multiplicative model for viability selection at multiple loci

Selection due to differential viability is studied in an n-locus two-allele model using a set indexation that allows the simplicity of the one-locus two-allele model to be carried to multi-locus models. The existence condition is analyzed for polymorphic equilibria with linkage equilibrium: Robbins' equilibria. The local stability condition is given for the Robbins' equilibria on the boundaries in the generalized non-epistatic selection regimes of Karlin and Liberman (1979). These generalized non-epistatic regimes include the additive selection model, the multiplicative selection model and the multiplicative interaction model, and their symmetric versions cover all the symmetric viability models.Research supported by grant no. 11-7805 from the Danish Natural Science Research Council, by NIH grant GM 28016, by a fellowship from the Research Foundation of Aarhus University, and by a visiting fellowship from the University of New England, N.S.W.     

Global convergence properties in multilocus viability selection models: the additive model and the Hardy-Weinberg law

A natural coordinate system is introduced for the analysis of the global stability of the Hardy-Weinberg (HW) polymorphism under the general multilocus additive viability model. A global convergence criterion is developed and used to prove that the HW polymorphism is globally stable when each of the loci is diallelic, provided the loci are overdominant and the multilocus recombination is positive. As a corollary the multilocus Hardy-Weinberg law for neutral selection is derived.Research supported in part by NIH grants GM 39907-01, GM 10452-26 and NSF Grant DMS 86-06244Research supported in part by a US-Israel Binational Science Foundation grant 85-00021 and NIH grant GM 28016     

Equilibria in an epistatic viability model under arbitrary strength of selection

A class of viability models that generalize the standard additive model for the case of pairwise additive by additive epistatic interactions is considered. Conditions for existence and stability of steady states in the corresponding two-locus model are analyzed. Using regular perturbation techniques, the case when selection is weaker than recombination and the case when selection is stronger than recombination are investigated. The results derived are used to make conclusions on the dependence of population characteristics on the relation between the strength of selection and the recombination rate.     

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.     

Evolution and polymorphism in the multilocus Levene model with no or weak epistasis

Evolution and the maintenance of polymorphism under the multilocus Levene model with soft selection are studied. The number of loci and alleles, the number of demes, the linkage map, and the degree of dominance are arbitrary, but epistasis is absent or weak. We prove that, without epistasis and under mild, generic conditions, every trajectory converges to a stationary point in linkage equilibrium. Consequently, the equilibrium and stability structure can be determined by investigating the much simpler gene-frequency dynamics on the linkage-equilibrium manifold. For a haploid species an analogous result is shown. For weak epistasis, global convergence to quasi-linkage equilibrium is established. As an application, the maintenance of multilocus polymorphism is explored if the degree of dominance is intermediate at every locus and epistasis is absent or weak. If there are at least two demes, then arbitrarily many multiallelic loci can be maintained polymorphic at a globally asymptotically stable equilibrium. Because this holds for an open set of parameters, such equilibria are structurally stable. If the degree of dominance is not only intermediate but also deme independent, and loci are diallelic, an open set of parameters yielding an internal equilibrium exists only if the number of loci is strictly less than the number of demes. Otherwise, a fully polymorphic equilibrium exists only nongenerically, and if it exists, it consists of a manifold of equilibria. Its dimension is determined. In the absence of genotype-by-environment interaction, however, a manifold of equilibria occurs for an open set of parameters. In this case, the equilibrium structure is not robust to small deviations from no genotype-by-environment interaction. In a quantitative-genetic setting, the assumptions of no epistasis and intermediate dominance are equivalent to assuming that in every deme directional selection acts on a trait that is determined additively, i.e., by nonepistatic loci with dominance. Some of our results are exemplified in this quantitative-genetic context.     

A multilocus-multiallele analysis of frequency-dependent selection induced by intraspecific competition

The study of the mechanisms that maintain genetic variation has a long history in population genetics. We analyze a multilocus-multiallele model of frequency- and density-dependent selection in a large randomly mating population. The number of loci and the number of alleles per locus are arbitrary. The n loci are assumed to contribute additively to a quantitative character under stabilizing or directional selection as well as under frequency-dependent selection caused by intraspecific competition. We assume the strength of stabilizing selection to be weak, whereas the strength of frequency dependence may be arbitrary. Density-dependence is induced by population regulation. Our main result is a characterization of the equilibrium structure and its stability properties in terms of all parameters. It turns out that no equilibrium exists with more than two alleles segregating per locus. We give necessary and sufficient conditions on the strength of frequency dependence to ensure the maintenance of multilocus polymorphism. We also give explicit formulas on the number of polymorphic loci maintained at equilibrium. These results are based on the assumption that selection is sufficiently weak compared with recombination, so that linkage equilibrium can be assumed. If additionally the population size is assumed to be constant, we prove that the dynamics of the model form a generalized gradient system. For the model in its general form we are able to derive necessary and sufficient conditions for the stability of the monomorphic equilibria. Furthermore, we briefly analyze a special symmetric two-locus two-allele model for a constant population size but allowing for linkage disequilibrium. Finally, we analyze a single diallelic locus with dominance to illustrate the complications that can occur if the assumption of additivity is relaxed.     

Random temporal variation in selection intensities: One-locus two-allele model

Summary The paper formulates a one locus two allele diploid model influenced by temporal random selection intensities. The concept of stochastic local stability for equilibria states is formulated. Necessary and sufficient conditions are derived for the stochastic local stability of the fixation states for non-dominant and dominant traits. A general model, where a fixed polymorphic equilibrium is maintained is investigated. The complete local evolutionary picture is derived and some cases of global convergence are established.Supported in part by N. I. H. Grant GM 10452-09 and N. S. F. Grant 7507129.     

A multilocus analysis of intraspecific competition and stabilizing selection on a quantitative trait

The equilibrium structure of an additive, diallelic multilocus model of a quantitative trait under frequency- and density-dependent selection is derived. The trait is under stabilizing selection and mediates intraspecific competition as induced, for instance, by differential resource utilization. It is assumed that stabilizing selection is weak, but the strength of competition may be arbitrary relative to it. Density dependence is caused by population regulation, which may be of a very general kind. The number and effects of loci are arbitrary, and stabilizing selection is not necessarily symmetric with respect to the range of phenotypic values. All previously studied models of intraspecific competition for a continuum of resources known to the author reduce to a special case of the present model if overall selection is weak. Therefore, in this case our results are applicable as approximations to all these models. Our central result is the (nearly) complete characterization of the equilibrium and stability structure in terms of all parameters. It is derived under the sole assumption that selection is weak enough relative to recombination to ignore linkage disequilibrium. In particular, necessary and sufficient conditions on the strength of competition relative to stabilizing selection are found that ensure the maintenance of multilocus polymorphism and the occurrence of disruptive selection. In this case, explicit formulas for the number of polymorphic loci at equilibrium, the allele frequencies, the genetic variance, and the strength of disruptive selection are obtained. For two loci, the effects of linkage are investigated analytically; for several loci, they are studied numerically.     

Quantitative variability and multilocus polymorphism under epistatic selection

We study multilocus polymorphism under selection, using a class of fitness functions that account for additive, dominant, and pairwise additive-by-additive epistatic interactions. The dynamic equations are derived in terms of allele frequencies and disequilibria, using the notions of marginal systems and marginal fitnesses, without any approximations. Stationary values of allele frequencies and pairwise disequilibria under weak selection are calculated by regular perturbation techniques. We derive conditions for existence and stability of the multilocus polymorphic states. Using these results, we then analyze a number of models describing stabilizing selection on additive characters, with some other factors, and determine the conditions under which genetic quantitative variability is maintained.     

General analysis of frequency-dependent sexual selection at a multi-allelic locus

A general model is analyzed in which arbitrarily frequency-dependent selection acts on one sex of a diploid population with several alleles at one locus, as a result of viability or mating-success differences. The existence of boundary and polymorphic equilibria is examined, and conditions for local stability, internal and external, are obtained. The status of Hardy-Weinberg approximations in studying stability and approach to equilibria is also considered. The general principles are then applied to two specific models: one where genotypes fall into two phenotypic classes; and one with a hierarchy of dominance where viability and sexual selection are opposed. In the latter case it is found that, of all the equilibria present, there is one and only one which could possibly be stable: the existence of a unique globally stable equilibrium might then be inferred.     

Long-term evolution of multilocus traits

 We analyze monomorphic equilibria of long-term evolution for one or two continuous traits, controlled by an arbitrary number of autosomal loci and subject to constant viability selection. It turns out that fitness maximization always obtains at long term equilibria, but in the case of two traits, linkage determines the precise nature of the fitness measure that is maximized. We then consider local convergence to long term equilibria, for two multilocus traits subject to either constant or frequency dependent selection. From a model of long-term dynamics near an equilibrium we derive a criterion of local long-term stability for 2-dimensional equilibria. It turns out that mutation can be a decisive factor for stability. Received 26 January 1994; received in revised form 26 September 1994     

Epistasis in the multiple locus symmetric viability model

The n-locus two-allele symmetric viability model is considered in terms of the parameters measuring the additive epistasis in fitness. The dynamics is analysed using a simple linear transformation of the gametic frequencies, and then the recurrence equations depend on the epistatic parameters and Geiringer's recombination distribution only. The model exhibits an equilibrium, the central equilibrium, where the 2 ⁿ gametes are equally frequent. The transformation simplifies the stability analysis of the central point, and provides the stability conditions in terms of the existence conditions of other equilibria. For total negative epistasis (all epistatic parameters are negative) the central point is stable for all recombination distributions. For free recombination either a central point (segregating one, two, ... or n loci) or the n-locus fixation states are stable. For no recombination and some epistatic parameters positive the central point is unstable and several boundary equilibria may be locally stable. The sign structure of the additive epistasis is therefore an important determinant of the dynamics of the n-locus symmetric viability model. The non-symmetric multiple locus models previously analysed are dynamically related, and they all have an epistatic sign structure that resembles that of the multiplicative viability model. A non-symmetric model with total negative epistasis which share dynamical properties with the similar symmetric model is suggested.Supported in part by NIH grant GM 28016, and by grant 81-5458 from the Danish Natural Science Research Council     

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.     

Dynamics of Genetic Variability in Two-Locus Models of Stabilizing Selection

We study a two locus model, with additive contributions to the phenotype, to explore the dynamics of different phenotypic characteristics under stabilizing selection and recombination. We demonstrate that the interaction of selection and recombination results in constraints on the mode of phenotypic evolution. Let V(g) be the genic variance of the trait and C(L) be the contribution of linkage disequilibrium to the genotypic variance. We demonstrate that, independent of the initial conditions, the dynamics of the system on the plane (V(g), C(L)) are typically characterized by a quick approach to a straight line with slow evolution along this line afterward. We analyze how the mode and the rate of phenotypic evolution depend on the strength of selection relative to recombination, on the form of fitness function, and the difference in allelic effect. We argue that if selection is not extremely weak relative to recombination, linkage disequilibrium generated by stabilizing selection influences the dynamics significantly. We demonstrate that under these conditions, which are plausible in nature and certainly the case in artificial stabilizing selection experiments, the model can have a polymorphic equilibrium with positive linkage disequilibrium that is stable simultaneously with monomorphic equilibria.     

Preferential mating in symmetric multilocus systems: stability conditions of the central equilibrium

Several multilocus models that incorporate both preferential mating and viability selection are studied. Specifically, a class of symmetric heterozygosity models are considered that assign individuals to phenotypic classes according to which loci are in heterozygous state regardless of the actual allelic content. Otherwise, an arbitrary number of loci, number of alleles per locus, and arbitrary recombination scheme, viability parameters and preferential mating pattern based on phenotypes are allowed. The conditions for the stability of a central polymorphism are indicated and interpreted. The effects of viability and preference selection may be summarized in a single quantity for each phenotypic class, a generalized fitness. Preferential assortative mating alone can produce stability for a central polymorphism as in the case of viability selection when sexual attractiveness or general fitness increases with higher levels of heterozygosity. The situation is more complex with sexual selection.     

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.     

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.