CENTRAL EQUILIBRIA IN MULTILOCUS SYSTEMS. I. GENERALIZED NONEPISTATIC SELECTION REGIMES

A natural class of multilocus recombination processes and related measures of crossover interference

Various classifications and representations of multilocus recombination structures are delimited. A class of recombination distributions called the count–location chiasma process is parametrized by a distribution of the number of crossover events and for such crossover events by a conditional distribution of crossover locations. A number of properties and examples of this recombination structure are developed connecting orderings among the recombination mapping functions and the nature of interference.

An Alternative to the Breeder’s and Lande’s Equations

The breeder’s equation is a cornerstone of quantitative genetics, widely used in evolutionary modeling. Noting the mean phenotype in parental, selected parents, and the progeny by E(Z₀), E(Z_W), and E(Z₁), this equation relates response to selection R = E(Z₁) − E(Z₀) to the selection differential S = E(Z_W) − E(Z₀) through a simple proportionality relation R = h²S, where the heritability coefficient h² is a simple function of genotype and environment factors variance. The validity of this relation relies strongly on the normal (Gaussian) distribution of the parent genotype, which is an unobservable quantity and cannot be ascertained. In contrast, we show here that if the fitness (or selection) function is Gaussian with mean μ, an alternative, exact linear equation of the form R′ = j²S′ can be derived, regardless of the parental genotype distribution. Here R′ = E(Z₁) − μ and S′ = E(Z_W) − μ stand for the mean phenotypic lag with respect to the mean of the fitness function in the offspring and selected populations. The proportionality coefficient j² is a simple function of selection function and environment factors variance, but does not contain the genotype variance. To demonstrate this, we derive the exact functional relation between the mean phenotype in the selected and the offspring population and deduce all cases that lead to a linear relation between them. These results generalize naturally to the concept of G matrix and the multivariate Lande’s equation Δz¯=GP−1S. The linearity coefficient of the alternative equation are not changed by Gaussian selection.

Analysis of central equilibrium configurations for certain multi-locus systems in subdivided populations

The multi-locus systems expressing non-epistatic and generalized symmetric selection lend themselves to the study of the stability of certain central polymorphic equilibria. These equilibria persist when any form of migration connects demes which share a common equilibrium. The analysis of the stability of the equilibrium in the global system is tractable, thus supplementing known protection results for two alleles at one locus with stability conditions on an internal equilibrium involving an arbitrary number of loci, each with an arbitrary number of alleles. Two of the principal findings are that stability of central Hardy–Weinberg type equilibria increase with ‘more’ migration and ‘more’ recombination. As a corollary, local stability in each deme implies stability in a system with migration superimposed; but instability in each deme when isolated does not imply instability when migration is superimposed.

On the meaning of non-epistatic selection

In population genetics, the additive and multiplicative viability models are often used for the quantitative description of models in which the genetic contributions of several different loci are independent; that is, there is no epistasis. Non-epistasis may also be quantitatively defined in terms of measures of interaction used widely in statistics. Setting these measures of epistasis to zero yields alternative definitions of non-epistasis. We show here that these two definitions of non-epistasis are equivalent; that is, in the most general case of a multilocus, multiallele system, the additive and multiplicative viability models are unique solutions of the additive and multiplicative conditions, respectively, for non-epistasis.

Lewontin and Kojima meet Fisher: linkage in a symmetric model of sex determination

The effect of linkage and epistasis on the evolution of the sex-ratio is studied in a symmetric two-locus model of autosomal sex determination closely related to the symmetric viability model of R. C. Lewontin and K. Kojima. R. A. Fisher's expectation of an even sex ratio for autosomal sex determination by a single gene governs the dynamics when the loci are tightly linked. However, recombination may preclude optimization of the sex ratio just as occurs in viability selection models. Many of the evolutionary phenomena known for the symmetric viability model also occur here. In addition, we exhibit a series of new phenomena related to the presence of surfaces of even sex ratio.

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.

Two-locus autosomal sex determination: on the evolutionary genetic stability of the even sex ratio

In two-locus models of sex determination, there are two kinds of interior (polymorphic) equilibria. One class has the even sex ratio, and the other has equal allele frequencies in the two sexes. Equilibria of the second class may exhibit linkage disequilibrium. The condition for external stability of these second-class equilibria to invasion by a new allele is that the appropriately averaged sex ratio near the equilibrium be moved closer to the even sex ratio than the average among the resident genotypes. However, invasion by a new chromosome depends on the recombination fraction in a way that appears to preclude general results about the evolutionary genetic stability of the even sex ratio in this situation.

Levels of multiallelic overdominance fitness, heterozygote excess and heterozygote deficiency

Concepts and results on selection balance in multiallelic systems are described. These include a multidimensional concept of heterozygote excess and heterozygote deficiency, a hierarchy of means of assessment of heterozygote advantage, comparisons and contrasts of allelic versus gametic polymorphic states, and conditions defining stable equilibria of complementary gametic sets. The concepts are illustrated in the context of viability selection and behavioral models of kin selection and for two major categories of multilocus selection regimes.

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 2n 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.

An evolutionary reduction principle for genetic modifiers

The joint evolution of major genes under viability selection and a modifier locus that controls recombination between the major genes, mutation at the major gene, or migration between two demes is studied. The modifying locus is selectively neutral and may have an arbitrary number of alleles. For each case a class of polymorphic equilibria exists in which the frequencies of the modifying alleles are those computed by assuming that the recombination, mutation, or migration rates were viabilities and in which the major and modifier loci are not statistically associated. These are called viability-analogous Hardy-Weinberg (VAHW) equilibria. A new allele introduced near these equilibria will enter the population if its marginal average rate of recombination, mutation, or migration (whichever applies) is less than the population average prior to its introduction. Stability properties of these VAHW equilibria are also reported.

Preferential mating in symmetric multilocus systems

A class of multilocus models that incorporate both preferential mating and viability selection is studied. Symmetry in alleles is supposed, resulting in the phenotypes being dependent only on the location of heterozygous loci. Otherwise, an arbitrary number of loci, number of alleles per locus, and arbitrary recombination schemes, viability parameters, and preferential mating pattern are allowed. The conditions for stability of a central polymorphism, c(*), are indicated and interpreted. Mating and viability parameters enter as one combined quantity for each phenotypic class, which represents a generalized fitness. The effect on stability of c(*) of increasing the number of alleles per locus and the number of loci requires the formulation to be set in terms of frequency-dependent preference parameters.

Representation of Nonepistatic selection models and analysis of multilocus Hardy-Weinberg Equilibrium configurations

The paper develops conditions for the existence and the stability of central equilibria emanating from selection recombination interaction with generalized nonepistatic selection forms operating in multilocus multiallele systems. The selection structure admits a natural representation as simple sums of Kronecker products based on a common set of marginal selection components. A flexible parametrization of the recombination process is introduced leading to a canonical derivation of the transformation equations connecting gamete frequency states over successive generations. Conditions for the existence and stability of multilocus Hardy-Weinberg (H.W.) type equilibria are elaborated for the classical nonepistatic models (multiplicative and additive viability effects across loci) as well as for generalized nonepistatic selection expressions. It is established that the range of recombination distributions maintaining a stable H.W. polymorphic equilibrium is confined to loose linkage in the pure multiplicative case, but is not restricted in the additive model. In the bisexual case we ascertain for the generalized nonepistatic model the stability conditions of a common H.W polymorphism.