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 consequences of gene flow for local adaptation and differentiation: a two-locus two-deme model

We consider a population subdivided into two demes connected by migration in which selection acts in opposite direction. We explore the effects of recombination and migration on the maintenance of multilocus polymorphism, on local adaptation, and on differentiation by employing a deterministic model with genic selection on two linked diallelic loci (i.e., no dominance or epistasis). For the following cases, we characterize explicitly the possible equilibrium configurations: weak, strong, highly asymmetric, and super-symmetric migration, no or weak recombination, and independent or strongly recombining loci. For independent loci (linkage equilibrium) and for completely linked loci, we derive the possible bifurcation patterns as functions of the total migration rate, assuming all other parameters are fixed but arbitrary. For these and other cases, we determine analytically the maximum migration rate below which a stable fully polymorphic equilibrium exists. In this case, differentiation and local adaptation are maintained. Their degree is quantified by a new multilocus version of $F_\mathrm{ST}$ and by the migration load, respectively. In addition, we investigate the invasion conditions of locally beneficial mutants and show that linkage to a locus that is already in migration-selection balance facilitates invasion. Hence, loci of much smaller effect can invade than predicted by one-locus theory if linkage is sufficiently tight. We study how this minimum amount of linkage admitting invasion depends on the migration pattern. This suggests the emergence of clusters of locally beneficial mutations, which may form ‘genomic islands of divergence’. Finally, the influence of linkage and two-way migration on the effective migration rate at a linked neutral locus is explored. Numerical work complements our analytical results.     

Polymorphism in the two-locus Levene model with nonepistatic directional selection

For the Levene model with soft selection in two demes, the maintenance of polymorphism at two diallelic loci is studied. Selection is nonepistatic and dominance is intermediate. Thus, there is directional selection in every deme and at every locus. We assume that selection is in opposite directions in the two demes because otherwise no polymorphism is possible. If at one locus there is no dominance, then a complete analysis of the dynamical and equilibrium properties is performed. In particular, a simple necessary and sufficient condition for the existence of an internal equilibrium and sufficient conditions for global asymptotic stability are obtained. These results are extended to deme-independent degree of dominance at one locus. A perturbation analysis establishes structural stability within the full parameter space. In the absence of genotype-environment interaction, which requires deme-independent dominance at both loci, nongeneric equilibrium behavior occurs, and the introduction of arbitrarily small genotype-environment interaction changes the equilibrium structure and may destroy stable polymorphism. The volume of the parameter space for which a (stable) two-locus polymorphism is maintained is computed numerically. It is investigated how this volume depends on the strength of selection and on the dominance relations. If the favorable allele is (partially) dominant in its deme, more than 20% of all parameter combinations lead to a globally asymptotically stable, fully polymorphic equilibrium.     

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.     

Multilocus selection in subdivided populations I. Convergence properties for weak or strong migration

The dynamics and equilibrium structure of a deterministic population-genetic model of migration and selection acting on multiple multiallelic loci is studied. A large population of diploid individuals is distributed over finitely many demes connected by migration. Generations are discrete and nonoverlapping, migration is irreducible and aperiodic, all pairwise recombination rates are positive, and selection may vary across demes. It is proved that, in the absence of selection, all trajectories converge at a geometric rate to a manifold on which global linkage equilibrium holds and allele frequencies are identical across demes. Various limiting cases are derived in which one or more of the three evolutionary forces, selection, migration, and recombination, are weak relative to the others. Two are particularly interesting. If migration and recombination are strong relative to selection, the dynamics can be conceived as a perturbation of the so-called weak-selection limit, a simple dynamical system for suitably averaged allele frequencies. Under nondegeneracy assumptions on this weak-selection limit which are generic, every equilibrium of the full dynamics is a perturbation of an equilibrium of the weak-selection limit and has the same stability properties. The number of equilibria is the same in both systems, equilibria in the full (perturbed) system are in quasi-linkage equilibrium, and differences among allele frequencies across demes are small. If migration is weak relative to recombination and epistasis is also weak, then every equilibrium is a perturbation of an equilibrium of the corresponding system without migration, has the same stability properties, and is in quasi-linkage equilibrium. In both cases, every trajectory converges to an equilibrium, thus no cycling or complicated dynamics can occur.      

The evolution of recombination in a heterogeneous environment

Most models describing the evolution of recombination have focused on the case of a single population, implicitly assuming that all individuals are equally likely to mate and that spatial heterogeneity in selection is absent. In these models, the evolution of recombination is driven by linkage disequilibria generated either by epistatic selection or drift. Models based on epistatic selection show that recombination can be favored if epistasis is negative and weak compared to directional selection and if the recombination modifier locus is tightly linked to the selected loci. In this article, we examine the joint effects of spatial heterogeneity in selection and epistasis on the evolution of recombination. In a model with two patches, each subject to different selection regimes, we consider the cases of mutation-selection and migration-selection balance as well as the spread of beneficial alleles. We find that including spatial heterogeneity extends the range of epistasis over which recombination can be favored. Indeed, recombination can be favored without epistasis, with negative and even with positive epistasis depending on environmental circumstances. The selection pressure acting on recombination-modifier loci is often much stronger with spatial heterogeneity, and even loosely linked modifiers and free linkage may evolve. In each case, predicting whether recombination is favored requires knowledge of both the type of environmental heterogeneity and epistasis, as none of these factors alone is sufficient to predict the outcome.     

Selection for recombination in structured populations

In finite populations, linkage disequilibria generated by the interaction of drift and directional selection (Hill-Robertson effect) can select for sex and recombination, even in the absence of epistasis. Previous models of this process predict very little advantage to recombination in large panmictic populations. In this article we demonstrate that substantial levels of linkage disequilibria can accumulate by drift in the presence of selection in populations of any size, provided that the population is subdivided. We quantify (i) the linkage disequilibrium produced by the interaction of drift and selection during the selective sweep of beneficial alleles at two loci in a subdivided population and (ii) the selection for recombination generated by these disequilibria. We show that, in a population subdivided into n demes of large size N, both the disequilibrium and the selection for recombination are equivalent to that expected in a single population of a size intermediate between the size of each deme (N) and the total size (nN), depending on the rate of migration among demes, m. We also show by simulations that, with small demes, the selection for recombination is stronger than both that expected in an unstructured population (m = 1 - 1/n) and that expected in a set of isolated demes (m = 0). Indeed, migration maintains polymorphisms that would otherwise be lost rapidly from small demes, while population structure maintains enough local stochasticity to generate linkage disequilibria. These effects are also strong enough to overcome the twofold cost of sex under strong selection when sex is initially rare. Overall, our results show that the stochastic theories of the evolution of sex apply to a much broader range of conditions than previously expected.     

The Evolution of the Y Chromosome with X-Y Recombination

A theoretical population genetic model is developed to explore the consequences of X-Y recombination in the evolution of sex chromosome polymorphism. The model incorporates one sex-determining locus and one locus subject to natural selection. Both loci have two alleles, and the rate of classical meiotic recombination between the loci is r. The alleles at the sex-determining locus specify whether the chromosome is X or Y, and the alleles at the selected locus are arbitrarily labeled A and a. Natural selection is modeled as a process of differential viabilities. The system can be expressed in terms of three recurrence equations, one for the frequency of A on the X-bearing gametes produced by females, one for each of the frequency of A on the X- and Y-bearing gametes produced by males. Several special cases are examined, including X chromosome dominance and symmetric selection. Unusual equilibria are found with the two sexes having very different allele frequencies at the selected locus. A significant finding is that the allowance of recombination results in a much greater opportunity for polymorphism of the Y chromosome. Tighter linkage results in a greater likelihood for equilibria with a large difference between the sex chromosomes in allele frequency.     

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.     

Consequences of genetic linkage for the maintenance of sexually antagonistic polymorphism in hermaphrodites

When selection differs between males and females, pleiotropic effects among genes expressed by both sexes can result in sexually antagonistic selection (SA), where beneficial alleles for one sex are deleterious for the other. For hermaphrodites, alleles with opposing fitness effects through each sex function represent analogous genetic constraints on fitness. Recent theory based on single‐locus models predicts that the maintenance of SA genetic variation should be greatly reduced in partially selfing populations. However, selfing also reduces the effective rate of recombination, which should facilitate selection on linked allelic combinations and expand opportunities for balancing selection in a multilocus context. Here, I develop a two‐locus model of SA selection for simultaneous hermaphrodites, and explore the joint influence of linkage, self‐fertilization, and dominance on the maintainance of SA polymorphism. I find that the effective reduction in recombination caused by selfing significantly expands the parameter space where SA polymorphism can be maintained relative to single‐locus models. In particular, linkage facilitates the invasion of male‐beneficial alleles, partially compensating for the “female‐bias” in the net direction of selection created by selfing. I discuss the implications of accounting for linkage among SA loci for the maintenance of SA genetic variation and mixed mating systems in hermaphrodites.     

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.     

The effects of linkage and gene flow on local adaptation: a two-locus continent-island model

Population subdivision and migration are generally considered to be important causes of linkage disequilibrium (LD). We explore the combined effects of recombination and gene flow on the amount of LD, the maintenance of polymorphism, and the degree of local adaptation in a subdivided population by analyzing a diploid, deterministic continent–island model with genic selection on two linked loci (i.e., no dominance or epistasis). For this simple model, we characterize explicitly all possible equilibrium configurations. Simple and intuitive approximations for many quantities of interest are obtained in limiting cases, such as weak migration, weak selection, weak or strong recombination. For instance, we derive explicit expressions for the measures and r² (the squared correlation in allelic state) of LD. They depend in qualitatively different ways on the migration rate. Remarkably high values of r² are maintained between weakly linked loci, especially if gene flow is low. We determine how the maximum amount of gene flow that admits preservation of the locally adapted haplotype, hence of polymorphism at both loci, depends on recombination rate and selection coefficients. We also investigate the evolution of differentiation by examining the invasion of beneficial mutants of small effect that are linked to an already present, locally adapted allele. Mutants of much smaller effect can invade successfully than predicted by naive single-locus theory provided they are at least weakly linked. Finally, the influence of linkage on the degree of local adaptation, the migration load, and the effective migration rate at a neutral locus is explored. We discuss possible consequences for the evolution of genetic architecture, in particular, for the emergence of clusters of tightly linked, slightly beneficial mutations and the evolution of recombination and chromosome inversions.     

Effects of Interference Between Selected Loci on the Mutation Load,Inbreeding Depression,and Heterosis

A classical prediction from single-locus models is that inbreeding increases the efficiency of selection against partially recessive deleterious alleles (purging), thereby decreasing the mutation load and level of inbreeding depression. However, previous multilocus simulation studies found that increasing the rate of self-fertilization of individuals may not lead to purging and argued that selective interference among loci causes this effect. In this article, I derive simple analytical approximations for the mutation load and inbreeding depression, taking into account the effects of interference between pairs of loci. I consider two classical scenarios of nonrandomly mating populations: a single population undergoing partial selfing and a subdivided population with limited dispersal. In the first case, correlations in homozygosity between loci tend to reduce mean fitness and increase inbreeding depression. These effects are stronger when deleterious alleles are more recessive, but only weakly depend on the strength of selection against deleterious alleles and on recombination rates. In subdivided populations, interference increases inbreeding depression within demes, but decreases heterosis between demes. Comparisons with multilocus, individual-based simulations show that these analytical approximations are accurate as long as the effects of interference stay moderate, but fail for high deleterious mutation rates and low dominance coefficients of deleterious alleles.     

The effect of subdivision on variation at multi-allelic loci under balancing selection

Simulations are used to investigate the expected pattern of variation at loci under different forms of multi-allelic balancing selection in a finite island model of a subdivided population. The objective is to evaluate the effect of restricted migration among demes on the distribution of polymorphism at the selected loci at equilibrium, and to compare the results with those expected for a neutral locus. The results show that the expected number of alleles maintained, and numbers of nucleotide differences between alleles, are relatively insensitive to the migration rate, and differentiation remains low even under very restricted migration. However, nucleotide divergence between copies of functionally identical alleles increases sharply when migration decreases. These results are discussed in relation to published surveys of allelic diversity in MHC and plant self-incompatibility systems, and to the possibility of inferring ancient population genetic events and processes. In addition, it is shown that, for sporophytic self-incompatibility systems, it is not necessarily true in a subdivided population that recessive alleles are more frequent than dominant ones.     

The average number of sites separating DNA sequences drawn from a subdivided population

The "infinite sites" model in the absence of recombination is examined in a subdivided population in which there is arbitrary migration among demes. It is shown that, if the migration matrix is symmetric and irreducible, the average number of sites that differ in two alleles chosen from the same deme depends only on an effective size of the whole population and not on either the elements of the migration matrix or the size of each deme separately. If there are n demes all of size N, the average number of sites that differ in two alleles chosen from the same deme is 4nN mu, where mu is the average mutation rate per site. This is the same value as for two alleles drawn from a panmictic population of size nN. The average number of sites that differ in alleles drawn from the same and from different demes can provide some information about the degree of population subdivision, as is illustrated by using the data of Kreitman and Aquadé (1986, Proc. Nat. Acad. Sci. U.S.A., 83, 3562) on Drosophila melanogaster.     

Spatial heterogeneity in the strength of selection against deleterious alleles and the mutation load

According to current estimates of genomic deleterious mutation rates (which are often of the order 0.1-1) the mutation load (defined as a reduction in the average fitness of a population due to the presence of deleterious alleles) may be important in many populations. In this paper, I use multilocus simulations to explore the effect of spatial heterogeneity in the strength of selection against deleterious alleles on the mutation load (for example, it has been suggested that stressful environments may increase the strength of selection). These simulations show contrasted results: in some situations, spatial heterogeneity may greatly reduce the mutation load, due to the fact that migrants coming from demes under stronger selection carry relatively few deleterious alleles, and benefit from a strong advantage within demes under weaker selection (where individuals carry many more deleterious alleles); in other situations, however, deleterious alleles accumulate within demes under stronger selection, due to migration pressure from demes under weaker selection, leading to fitness erosion within those demes. This second situation is more frequent when the productivity of the different demes is proportional to their mean fitness. The effect of spatial heterogeneity is greatly reduced, however, when the response to environmental differences is inconsistent across loci.     

Molecular mechanisms of dominance evolution in Müllerian mimicry

Natural selection acting on dominance between adaptive alleles at polymorphic loci can be sufficiently strong for dominance to evolve. However, the molecular mechanisms underlying such evolution are generally unknown. Here, using Müllerian mimicry as a case‐study for adaptive morphological variation, we present a theoretical analysis of the invasion of dominance modifiers altering gene expression through different molecular mechanisms. Toxic species involved in Müllerian mimicry exhibit warning coloration, and converge morphologically with other toxic species of the local community, due to positive frequency‐dependent selection acting on these colorations. Polymorphism in warning coloration may be maintained by migration–selection balance with fine scale spatial heterogeneity. We modeled a dominance modifier locus altering the expression of the warning coloration locus, targeting one or several alleles, acting in cis or trans, and either enhancing or repressing expression. We confirmed that dominance could evolve when balanced polymorphism was maintained at the color locus. Dominance evolution could result from modifiers enhancing one allele specifically, irrespective of their linkage with the targeted locus. Nonspecific enhancers could also persist in populations, at frequencies tightly depending on their linkage with the targeted locus. Altogether, our results identify which mechanisms of expression alteration could lead to dominance evolution in polymorphic mimicry.     

The potential for sexually antagonistic polymorphism in different genome regions

Sex differences in the fitness effects of alleles at a single locus (intralocus sexual antagonism, or SA) have several evolutionary consequences. Among the consequences of SA, polymorphisms at genes partially linked to the sex-determining region of the sex chromosome pair potentially drive the evolution of suppressed recombination between the sex chromosomes. Understanding the conditions under which SA polymorphism can exist at such pseudo-autosomal (or PAR) loci should increase understanding of the evolution of recombination between sex chromosome pairs, and can help predict when we may expect potentially empirically detectable allele frequency differences between the sexes. Models so far published have concluded that PAR genes can maintain SA polymorphisms over a wider range of selection coefficients than autosomal ones, but have used restrictive assumptions. We expand the modeling of SA alleles at a single locus with the full range of degrees of linkage to the male-specific region, to include strong or weak selection and the possibility of different dominance coefficients in the two sexes. We confirm the previous major conclusion that SA polymorphisms are generally maintained in a larger region of parameter space if the locus is in the PAR than if it is autosomal.     

A chip off the old block: a model for the evolution of genomic imprinting via selection for parental similarity

A consequence of genomic imprinting is that offspring are more similar to one parent than to the other, depending on which parent's genes are inactivated in those offspring. We hypothesize that genomic imprinting may have evolved at some loci because of selection to be similar to the parent of one sex or the other. We construct and analyze an evolutionary-genetic model of a two-locus two-deme system, in which one locus codes for a character under local selection and the second locus is a potential cis-acting modifier of imprinting. A proportion of males only migrate between demes every generation, and prebreeding males are less fit, on average, than females. We examine the conditions in which an imprinting modifier allele can invade a population fixed for a nonimprinting modifier allele and vice versa. We find that the conditions under which the imprinting modifier invades are biologically restrictive (high migration rates and high values of recombination between the two loci) and thus this hypothesis is unlikely to explain the evolution of imprinting. Our modeling also shows that, as with several other hypotheses, polymorphism of imprinting status may evolve under certain circumstances, a feature not predicted by verbal accounts.     

Selective sweeps in multilocus models of quantitative traits

We study the trajectory of an allele that affects a polygenic trait selected toward a phenotypic optimum. Furthermore, conditioning on this trajectory we analyze the effect of the selected mutation on linked neutral variation. We examine the well-characterized two-locus two-allele model but we also provide results for diallelic models with up to eight loci. First, when the optimum phenotype is that of the double heterozygote in a two-locus model, and there is no dominance or epistasis of effects on the trait, the trajectories of selected mutations rarely reach fixation; instead, a polymorphic equilibrium at both loci is approached. Whether a polymorphic equilibrium is reached (rather than fixation at both loci) depends on the intensity of selection and the relative distances to the optimum of the homozygotes at each locus. Furthermore, if both loci have similar effects on the trait, fixation of an allele at a given locus is less likely when it starts at low frequency and the other locus is polymorphic (with alleles at intermediate frequencies). Weaker selection increases the probability of fixation of the studied allele, as the polymorphic equilibrium is less stable in this case. When we do not require the double heterozygote to be at the optimum we find that the polymorphic equilibrium is more difficult to reach, and fixation becomes more likely. Second, increasing the number of loci decreases the probability of fixation, because adaptation to the optimum is possible by various combinations of alleles. Summaries of the genealogy (height, total length, and imbalance) and of sequence polymorphism (number of polymorphisms, frequency spectrum, and haplotype structure) next to a selected locus depend on the frequency that the selected mutation approaches at equilibrium. We conclude that multilocus response to selection may in some cases prevent selective sweeps from being completed, as described in previous studies, but that conditions causing this to happen strongly depend on the genetic architecture of the trait, and that fixation of selected mutations is likely in many instances.