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.     

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.     

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.     

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.     

A theoretical and numerical assessment of genetic variability

The equilibrium behavior of one-locus viability selection models is studied numerically. The selection schemes include randomly chosen viabilities, viabilities chosen to measure a hypothetical distance between the alleles making up the genotype and viabilities that obey various allelic dominance relations. From 3 to 8 alleles are considered. Among the key conclusions are (1) equilibria that are most polymorphic do not usually have the highest mean fitness, (2) the more structure there is in the choice of the viability model, the greater is the level of polymorphism at equilibrium, and (3) for the numbers of alleles chosen here, the equilibrium reached by iteration from the centroid of the allele frequency simplex is the best predictor of the equilibrium attainable from randomly chosen starting vectors. Preliminary evidence shows that this is not the case for 16 alleles.     

Deviations of Genotypic Structures from Hardy-Weinberg Proportions under Random Mating and Differential Selection between the Sexes

Population genetic models, such as differential viability selection between the sexes and differential multiplicative fecundity contributions of the sexes, are considered for a single multiallelic locus. These selection models usually produce deviations of the zygotic genotype frequencies from Hardy-Weinberg proportions. The deviations are investigated (with special emphasis put on equilibrium states) to quantify the effect of selective asymmetry in the two sexes. For many selection regimes, the present results demonstrate a strong affinity of zygotic genotype frequencies for Hardy-Weinberg proportions after two generations, at the latest. It is shown that the deviations of genotypic equilibria from the corresponding Hardy-Weinberg proportions can be expressed and estimated by means of selection components of only that sex with the lower selection intensity. This corresponds to the well-known fact that viability selection acting in only one sex yields Hardy-Weinberg equilibria.     

DYNAMICS OF SEXUAL SELECTION IN DIPLOID POPULATIONS

We describe results for a diploid, two-locus model for the evolution of a female mating preference directed at an attractive male trait that is subject to viability and/or fertility selection. Using computer simulation, we studied a large, random sample of parameter values, assuming additivity of alleles at the preference locus and partial dominance at the trait locus. Simulation results were classifiable into nine types of parameter sets, each differing in equilibria, evolutionary trajectories, and rates of evolution. For many parameters, evolutionary trajectories converged on curves within the allelic frequency plane and subsequently evolved along the curves toward fixation. Neutrally stable curves of equilibria did not occur in Fisherian models that assume only viability and sexual selection unless there is complete dominance at the trait locus. The Fisherian models also exhibited oscillation of allelic frequencies and unique polymorphic equilibria. “Sexy son” models in which attractive males had reduced fertility were much less likely to lead to increase in traits and preferences than were the Fisherian models. However, if less fertile males had increased viability, trait polymorphisms and fixation of rare “sexy” alleles occurred. In general, the behavior of the diploid model was much more complex than that of analogous haploid or polygenic models.     

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.     

Population models of genomic imprinting. I. Differential viability in the sexes and the analogy with genetic dominance

Many single-locus, two-allele selection models of genomic imprinting have been shown to reduce formally to one-locus Mendelian models with a modified parameter for genetic dominance. One exception is the model where selection at the imprinted locus affects the sexes differently. We present two models of maternal inactivation with differential viability in the sexes, one with complete inactivation, and the other with a partial penetrance for inactivation. We show that, provided dominance relations at the imprintable locus are the same in both sexes, a globally stable polymorphism exists for a range of viabilities that is independent of the penetrance of imprinting. The conditions for a polymorphism are the same as in previous models with differential viability in the sexes but without imprinting and in a model of the paternal X-inactivation system in marsupials. The model with incomplete inactivation is used to illustrate the analogy between imprinting and dominance by comparing equilibrium bifurcation plots for fixed values of dominance and penetrance. We also derive a single expression for the dominance parameter that leaves the frequency and stability of equilibria unchanged for all levels of inactivation. Although an imprinting model with sex differences does not formally reduce to a nonimprinting scheme, close theoretical parallels clearly exist.     

The effect of variable environments on polymorphism at loci with several alleles I. A symmetric model

Much of the extant polymorphism has been attributed to spatial and temporal variation among selection regimes. Analysis of models entailing two alleles at a single locus has demonstrated that polymorphism may result from variation among selection regimes which prescribe monomorphism if constant. This relationship is studied in the context of several alleles at a locus.One result which is not valid with only two alleles is that variation among selection regimes which specify polymorphic equilibria may lead to a stable monomorphic equilibrium. The analyses of temporal variation and total panmixia spatial variation among environments show that temporal variation allows the simultaneous stability of equilibrium configurations which cannot be simultaneously stable under total panmixia spatial variation (hard or soft selection). The principle that polymorphism is more readily maintained with spatial than temporal variation is invalidated.Supported in part by Purdue Research Foundation and National Science Foundation (USA) grant MCS-8002227     

Polymorphism in multiallelic migration–selection models with dominance

The evolution of the gene frequencies at a single multiallelic locus under the joint action of migration and viability selection with dominance is investigated. The monoecious, diploid population is subdivided into finitely many panmictic colonies that exchange adult migrants independently of genotype. Underdominance and overdominance are excluded. If the degree of dominance is deme independent for every pair of alleles, then under the Levene model, the qualitative evolution of the gene frequencies (i.e., the existence and stability of the equilibria) is the same as without dominance. In particular: (i) the number of demes is a generic upper bound on the number of alleles present at equilibrium; (ii) there exists exactly one stable equilibrium, and it is globally attracting; and (iii) if there exists an internal equilibrium, it is globally asymptotically stable. Analytic examples demonstrate that if either the Levene model does not apply or the degree of dominance is deme dependent, then the above results can fail. A complete global analysis of weak migration and weak selection on a recessive allele in two demes is presented.     

Analysis of some nonrandom mating models

In this paper a few asymmetric models are presented taking account of the effects of assortative mating on an autosomal trait controlled by a single locus possibly with multiple alleles. The models are developed by specifying the intensities for preference mating for various phenotypes. The analysis is confined to the case in which preference is exercised by the individuals of one sex only. It is assumed that males possess unlimited fertility.The dynamics of the population and its equilibrium distribution are discussed. The gene frequency usually changes with time and equilibrium distribution in most cases depends only on the assortment parameters. Expressions are obtained giving the additive and dominance components of variance, and covariances for relatives for populations in equilibrium for some of the models.     

Gene-Culture Coevolution and Sex-Ratios: The Effects of Infanticide, Sex-Selective Abortion, Sex Selection, and Sex-Biased Parental Investment on the Evolution of Sex Ratios

The evolutionary consequences of culturally transmitted practices that cause differential mortality between the sexes, thereby distorting the sex ratio (e.g., female infanticide and sex-selective abortion), are explored using dynamic models of gene-culture coevolution. We investigate how a preference for the sex of offspring may affect the selection of genes distorting the primary sex ratio. Sex-dependent differences in mortality have been predicted to select for a male- or female-biased primary sex ratio, to have no effect, or to favor either under different circumstances. We find that when a mating pair′s behavior modifies mortality rates in favor of one sex, but does not change the number of offspring produced in the mating, the primary sex ratio will evolve a bias against the favored sex However, when the total number of offspring of a mating pair is significantly seduced as a consequence of their prejudice, the primary sex ratio will evolve to favor the preferred sex. These results hold irrespective of whether the sex ratio is distorted by the mother′s, the father′s or the individual′s own autosomal genes. The use of dynamic models of gene-culture coevolution allows us to explore the evolution of alleles which distort the sex ratio, as well as the final equilibrium states of the system. Gene-culture interactions can provide equilibria different from those in purely genetic systems, slow the approach to these equilibria by orders of magnitude, and move the primary (PSR) and the adult sex ratio (ASR) away from any stable equilibrium for hundreds of generations.     

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.     

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.     

Protected polymorphism and evolutionary stability in pleiotropic models with trait-specific dominance

When alleles have pleiotropic effects on a number of quantitative traits, the degree of dominance between a pair of alleles can be different for each trait. Such trait-specific dominance has been studied previously in models for the maintenance of genetic variation by antagonistic effects of an allele on two fitness components. By generalizing these models to an arbitrary number of fitness components or other phenotypic traits with different degrees of dominance, I show that genetic polymorphism is generally impossible without antagonistic fitness effects of different traits and without trait-specific dominance. I also investigate dominance and pleiotropy from a more long-term evolutionary perspective, allowing for the study of general ecological scenarios, and I discuss the effects of trait-specific dominance on evolutionary stability criteria. When selection is mainly directional and only trait-specific dominance and antagonism cause the emergence of polymorphism, then these polymorphisms can be overtaken by single mutants again, such that they are probably short-lived on an evolutionary time scale. Near evolutionarily singular points where directional selection is absent, trait-specific dominance and overdominance facilitate the emergence of polymorphism and cause evolutionary divergence in some cases. An important outcome of these models is that trait-specific dominance allows for the emergence of genetic polymorphisms without a selective disadvantage for heterozygotes. This removes the scope for the evolution of assortative mate choice and affects dominance modification. Sympatric speciation by disruptive ecological selection requires this heterozygote disadvantage in order to evolve, and therefore it becomes less plausible if the emergence of genetic polymorphism usually occurs via trait-specific dominance and antagonistic effects.     

The resolution of sexual antagonism by gene duplication

Disruptive selection between males and females can generate sexual antagonism, where alleles improving fitness in one sex reduce fitness in the other. This type of genetic conflict arises because males and females carry nearly identical sets of genes: opposing selection, followed by genetic mixing during reproduction, generates a population genetic "tug-of-war" that constrains adaptation in either sex. Recent verbal models suggest that gene duplication and sex-specific cooption of paralogs might resolve sexual antagonism and facilitate evolutionary divergence between the sexes. However, this intuitive proximal solution for sexual dimorphism potentially belies a complex interaction between mutation, genetic drift, and positive selection during duplicate fixation and sex-specific paralog differentiation. The interaction of these processes--within the explicit context of duplication and sexual antagonism--has yet to be formally described by population genetics theory. Here, we develop and analyze models of gene duplication and sex-specific differentiation between paralogs. We show that sexual antagonism can favor the fixation and maintenance of gene duplicates, eventually leading to the evolution of sexually dimorphic genetic architectures for male and female traits. The timescale for these evolutionary transitions is sensitive to a suite of genetic and demographic variables, including allelic dominance, recombination, sex linkage, and population size. Interestingly, we find that female-beneficial duplicates preferentially accumulate on the X chromosome, whereas male-beneficial duplicates are biased toward autosomes, independent of the dominance parameters of sexually antagonistic alleles. Although this result differs from previous models of sexual antagonism, it is consistent with several findings from the empirical genomics literature.     

A Simulation Study of Truncation Selection for a Quantitative Trait Opposed by Natural Selection

A quantitative character controlled at one locus with two alleles was submitted to artificial (mass) selection and to three modes of opposing natural selection (directional selection, overdominance and underdominance) in a large random-mating population. The selection response and the limits of the selective process were studied by deterministic simulation. The lifetime of the process was generally between 20 and 100 generations and did not appear to depend on the mode of natural selection. However, depending on the values of the parameters (initial gene frequency, selection intensity, ratio of the effect of the gene to the environmental standard deviation, fitness values) the following outcomes of selection were observed: fixation of the allele favored by artificial selection, stable nontrivial equilibrium, unstable equilibrium and loss of the allele favored by artificial selection. Finally, the results of the simulation were compared to the results of selection experiments.     

Sexual selection and fertility

Genetic models are analyzed in which sexual selection is combined with fertility selection. In these models, the sexual selection acts on males, the fertility selection on either males, females or both sexes. The phenotypes thus selected may be determined either by dominant and recessive alleles or by each homozygous and heterozygous genotype. Polymorphisms of dominant and recessive phenotypes can be maintained in equilibrium by a balance between sexual and fertility selection. Generally fertility selection has a greater effect than viability selection in determining the point of equilibrium. The dominant phenotype is maintained at a lower frequency when at a fertility disadvantage than when at a viability disadvantage. When about 20% or more of the females mate preferentially, the models show that equilibria will be established at very different frequencies depending on whether fertility selection acts on males, females or both sexes. These results, applied to data of preferential mating of melanic two-spot ladybirds, predict differences in fertility which can be use to test the models. Symmetric models of preferences for each genotype also give rise to polymorphisms if the heterozygotes obtain an overall advantage.