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.     

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.     

Two loci with two alleles: Linkage equilibrium and linkage disequilibrium can be simultaneously stable

A class of two-locus two-allele viability matrices is exhibited which has the property that, for a large range of recombination values, both linkage equilibrium and linkage disequilibrium are stable. The model is specified by five viabilities; the classical schemes previously analyzed involve at most four selection parameters.     

A numerical solution to the equilibria of the two-locus two-allele selection model.

Examination of the equilibria of the standard two-locus two-allele selection model leads to the construction of a polynomial with coefficients derived from selective values in the genotypic fitness matrix. This polynomial can be partially factored algebraically and numerical techniques are available to extract the roots of the remainder. Each root provides a possible value of the disequilibrium coefficient and the gametic frequencies at equilibrium, and these can be readily checked for stability.     

Linkage and the Maintenance of Heritable Variation by Mutation-Selection Balance

The role of linkage in influencing heritable variation maintained through a balance between mutation and stabilizing selection is investigated for two different models. In both cases one trait is considered and the interactions within and between loci are assumed to be additive. Contrary to most earlier investigations of this problem no a priori assumptions on the distribution of genotypic values are imposed. For a deterministic two-locus two-allele model with recombination and mutation, related to the symmetric viability model, a complete nonlinear analysis is performed. It is shown that, depending on the recombination rate, multiple stable equilibria may coexist. The equilibrium genetic and genic variances are calculated. For a polygenic trait in a finite population with a possible continuum of allelic effects a simulation study is performed. In both models the equilibrium genetic and genic variances are roughly equal to the house-of-cards prediction or its finite population counterpart as long as the recombination rate is not extremely low. However, negative linkage disequilibrium builds up. If the loci are very closely linked the equilibrium additive genetic variance is slightly lower than the house-of-cards prediction, but the genic variance is much higher. Depending on whether the parameters are in favor of the house-of-cards or the Gaussian approximation, different behavior of the genetic system occurs with respect to linkage.     

MULTIPLE EQUILIBRIA AND MAINTENANCE OF ADDITIVE GENETIC VARIANCE IN A MODEL OF PLEIOTROPY

We describe a multilocus model that incorporates pleiotropic stabilizing selection on a large number of characters. We find many different stable equilibria with different levels of polymorphism and additive genetic variability. The results lend support to Wright's concept of a complex adaptive surface with many peaks of different heights. The model assumes that alleles contribute additively to the characters. We analyze the multilocus model by first considering a two-locus model. The two-locus model depends critically on having loci of different effect and on having the optimum phenotype not be that of a completely heterozygous individual. The effects of different loci need to differ only by less than a factor of two. For the multilocus, multicharacter model, we assume that completely heterozygous individuals do not have the optimum phenotype. By restricting attention to a two-allele model, we also assume that there are no alleles that can affect all characters in all possible combinations of directions.     

A complete enumeration and classification of two-locus disease models

There are 512 two-locus, two-allele, two-phenotype, fully penetrant disease models. Using the permutation between two alleles, between two loci, and between being affected and unaffected, one model can be considered to be equivalent to another model under the corresponding permutation. These permutations greatly reduce the number of two-locus models in the analysis of complex diseases. This paper determines the number of nonredundant two-locus models (which can be 102, 100, 96, 51, 50, or 58, depending on which permutations are used, and depending on whether zero-locus and single-locus models are excluded). Whenever possible, these nonredundant two-locus models are classified by their property. Besides the familiar features of multiplicative models (logical AND), heterogeneity models (logical OR), and threshold models, new classifications are added or expanded: modifying-effect models, logical XOR models, interference and negative interference models (neither dominant nor recessive), conditionally dominant/recessive models, missing lethal genotype models, and highly symmetric models. The following aspects of two-locus models are studied: the marginal penetrance tables at both loci, the expected joint identity-by-descent (IBD) probabilities, and the correlation between marginal IBD probabilities at the two loci. These studies are useful for linkage analyses using single-locus models while the underlying disease model is two-locus, and for correlation analyses using the linkage signals at different locations obtained by a single-locus model.     

Monotonic Change of the Mean Phenotype in Two-Locus Models

It is shown that the mean phenotype monotonically approaches the optimum in a class of symmetric, two-locus, two-allele models with stabilizing selection. In this model, mean fitness does not change monotonically. Thus, Fisher's fundamental theorem does not hold, even though another quantity of evolutionary interest, the mean phenotype, can be shown to change monotonically. Using this quantity, it is proven that global stability results for this model.     

Phenogenetic drift and the evolution of genotype-phenotype relationships

Maintenance of a stable two-locus polymorphism is analyzed statistically by fitting a logistic regression with a quadratic function of genotypic fitnesses to the probability for a fitness set to maintain a polymorphism. The regression is fitted using a data set containing information on stable equilibria maintained by 32,00 randomly generated fitness sets with three recombination values (0. 005, 0.05, 0.5). Fitted logistic regressions discriminate with 88 to 90% accuracy between fitness sets maintaining and not maintaining a stable internal equilibrium, which implies the existence of a fitness structure (balance of fitnesses) maintaining a two-locus polymorphism. Aspects of the balance of fitnesses revealed by logistic regressions are discussed. It is demonstrated that logistic regression also discriminates between types of a stable polymorphism: globally stable polymorphism, several simultaneously stable polymorphisms, and stable equilibria in addition to a polymorphic one, which implies that different balances of fitnesses are responsible for the maintenance of different types of polymorphism.     

Simultaneous Stability of D=0 and D not equal0 for Multiplicative Viabilities at Two Loci

The two-locus, two-allele multiplicative viability model is investigated. It is shown that the well-known region of recombination values for which D = 0 is locally stable does not preclude the local stability of an equilibrium with D ≠ 0. This is shown numerically and is true for every case investigated in which both loci are overdominant and the viabilities not symmetric.     

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.     

Simulation of ‘hitch-hiking’ genealogies

An ancestral influence graph is derived, an analogue of the coalescent and a composite of Griffiths' (1991) two-locus ancestral graph and Krone and Neuhauser's (1997) ancestral selection graph. This generalizes their use of branching-coalescing random graphs so as to incorporate both selection and recombination into gene genealogies. Qualitative understanding of a ‘hitch-hiking’ effect on genealogies is pursued via diagrammatic representation of the genealogical process in a two-locus, two-allele haploid model. Extending the simulation technique of Griffiths and Tavaré (1996), computational estimation of expected times to the most recent common ancestor of samples of n genes under recombination and selection in two-locus, two-allele haploid and diploid models are presented. Such times are conditional on sample configuration. Monte Carlo simulations show that ‘hitch-hiking’ is a subtle effect that alters the conditional expected depth of the genealogy at the linked neutral locus depending on a mutation-selection-recombination balance. Received: 21 July 2000 / Published online: 5 December 2000     

Disequilibrium, Selection, and Recombination: Limits in Two-Locus, Two-Allele Models

All possible combinations of equilibria and fitnesses in two-locus, two-allele, deterministic, discrete-generation selection models are enumerated. This knowledge is used to obtain limits (which can be calculated to arbitrary precision) to the relationships among disequilibrium, selection and recombination for fixed values of allele frequencies. In all cases, the inequality|rD| < s/10 holds, where r is recombination and D is disequilibrium, and all selection coefficients lie between 1 - s and 1 + s times that of the double heterozygote. Linear programming techniques are used to calculate the minimum strength of selection needed to explain several observed nonzero values of D reported in the literature. One conclusion is that the failure to observe nonzero values of D is not surprising.     

The role of recombination and selection in the modifier theory of sex-ratio distortion

The equilibrium configurations for a two-locus multialle model of sex-linked meiotic drive are studied with regard to the recombination fraction:limit cycles can occur in the case of small recombination while stable equilibrium points associated with linkage equilibrium can exist for an intermediate range of recombination values depending on the equilibrium sex ratio, linkage disequilibrium at nearby equilibrium points taking turn with loser linkage. The evolutionary dynamics in two-locus sex-ratio distortion systems is enlightened: while equilibria with a sex ratio closer to 1/2 are more likely to be stable with respect to perturbations on the frequencies of sex-ratio distorters that are represented at equilibrium, such equilibria are also more vulnerable to the invasion of mutant distorters when there is some degree of linkage with the sex-determining locus. For X-linked multimodifier systems of sex-ratio distortion, differential fertilities and viabilities are incorporated and a maximum principle is suggested.     

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.     

Multilocus Population Genetics with Weak Epistasis. I. Equilibrium Properties of Two-Locus Two-Allele Models

Using perturbation techniques, I determine the equilibrium of two-locus two-allele models with overdominance and weak epistasis. To lowest order, the allele frequencies, the mean fitness and the covariance between heterokaryotic and homokaryotic flies arising in the Sturtevant experimental design are independent of the recombination rate, r. The disequilibrium varies as one divided by the recombination rate, in contrast to neutral models. Although the disequilibrium generated by weak epistasis is small, too small to be experimentally detected, it can be large enough to have biological importance.     

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.     

Conditions for the existence of stationary densities for some two-dimensional diffusion processes with applications in population biology

Conditions are derived that we conjecture are necessary and sufficient for the existence of stationary densities for a class of two-dimensional diffusion processes. The derivation of the conditions rests on the assumption that a two-dimensional stationary density (which can be viewed as a stable “internal equilibrium”) exists if and only if all “boundary equilibria” are unstable in the sense that small perturbations lead to moving away from the boundaries with high probability. For the models considered, the boundary equilibria are one-dimensional stationary densities and equilibrium points. To demonstrate the usefulness of the conditions, three random environment models are analyzed: a three-allele selection model, a two-species competition model, and a two-locus selection model. Several of the results obtained have been verified by alternate methods.     

Limits to the Relationship among Recombination, Disequilibrium and Epistasis in Two-Locus Models

I determine limits to the equilibrium relationship among epistasis, recombination and disequilibrium in two-locus, two-allele models using linear programming techniques. I show that when allele frequencies are one-half at each locus, the symmetric model is the fitness pattern that generates the most disequilibrium for the smallest level of epistasis. When allele frequencies deviate from one-half much larger levels of epistasis are required to generate similar levels of disequilibrium. I determine the level of epistasis required to generate observed significant levels of disequilibrium in natural populations. The overall implication is that disequilibrium will be large at equilibrium only between strongly interacting, closely linked loci.     

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.