Long-term evolution of polygenic traits under frequency-dependent intraspecific competition

We analytically investigate the long-term evolution of a continuously varying quantitative character in a diploid population that is determined additively by a finite number of loci. The trait is under a mixture of frequency-dependent disruptive selection induced by intraspecific competition and frequency-independent stabilizing selection. Moreover, the trait is restricted to a finite range by constraints on the particular loci. Our investigations are based on explicit analytical results (provided by Bürger [2005. A multilocus analysis of intraspecific competition and stabilizing selection on a quantitative trait. J. Math. Biol. 50, 355-396]; Schneider [2006. A multilocus-multiallele analysis of frequency-dependent selection induced by intraspecific competition. J. Math. Biol. 52, 483-523]) on the short-term dynamics under the assumption of linkage equilibrium. We show that the population always reaches a long-term equilibrium (LTE), i.e., an equilibrium that is resistant against perturbations of mutations of sufficiently small effect. In general, several LTEs can coexist. They can be calculated explicitly, and we provide necessary and sufficient conditions for their existence. In the case that more than one LTE exists, we exemplify numerically that the evolutionary outcome depends crucially on the initial genetic architecture, on the joint distribution of mutational effects across loci, and on the particular realization of the mutation process. Therefore, long-term evolution cannot be predicted from the ecology alone. We further show that a partial order exists for the LTEs. The set of LTEs has a 'largest' element, an LTE which is reached during long-term evolution if the effects of the occurring mutant alleles are sufficiently large.     

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.     

Competitive divergence in non-random mating populations

A haploid model of frequency-dependent selection and assortative mating is introduced and analyzed for the case of a single multiallelic autosomal locus. Frequency-dependent selection is due to intraspecific competition mediated by a quantitative character under stabilizing or directional selection. Assortment is induced by the same trait. We analyze the equilibrium structure and the local stability properties of all possible equilibria. In the limit of weak selection we obtain global stability properties by finding a Lyapunov function. We provide necessary and sufficient conditions for the maintenance of polymorphism in terms of the strength of stabilizing selection, intraspecific competition and assortment. Our results also include criteria for the ability of extreme types to invade the population. Furthermore, we study the occurrence of disruptive selection and provide necessary and sufficient conditions for intraspecific divergence to occur.     

The effects of intraspecific competition and stabilizing selection on a polygenic trait

The equilibrium properties of an additive multilocus model of a quantitative trait under frequency- and density-dependent selection are investigated. Two opposing evolutionary forces are assumed to act: (i) stabilizing selection on the trait, which favors genotypes with an intermediate phenotype, and (ii) intraspecific competition mediated by that trait, which favors genotypes whose effect on the trait deviates most from that of the prevailing genotypes. Accordingly, fitnesses of genotypes have a frequency-independent component describing stabilizing selection and a frequency- and density-dependent component modeling competition. We study how the equilibrium structure, in particular, number, degree of polymorphism, and genetic variance of stable equilibria, is affected by the strength of frequency dependence, and what role the number of loci, the amount of recombination, and the demographic parameters play. To this end, we employ a statistical and numerical approach, complemented by analytical results, and explore how the equilibrium properties averaged over a large number of genetic systems with a given number of loci and average amount of recombination depend on the ecological and demographic parameters. We identify two parameter regions with a transitory region in between, in which the equilibrium properties of genetic systems are distinctively different. These regions depend on the strength of frequency dependence relative to pure stabilizing selection and on the demographic parameters, but not on the number of loci or the amount of recombination. We further study the shape of the fitness function observed at equilibrium and the extent to which the dynamics in this model are adaptive, and we present examples of equilibrium distributions of genotypic values under strong frequency dependence. Consequences for the maintenance of genetic variation, the detection of disruptive selection, and models of sympatric speciation are discussed.     

The conditions for speciation through intraspecific competition

Abstract It has been shown theoretically that sympatric speciation can occur if intraspecific competition is strong enough to induce disruptive selection. However, the plausibility of the involved processes is under debate, and many questions on the conditions for speciation remain unresolved. For instance, is strong disruptive selection sufficient for speciation? Which roles do genetic architecture and initial composition of the population play? How strong must assortative mating be before a population can split in two? These are some of the issues we address here. We investigate a diploid multilocus model of a quantitative trait that is under frequency‐dependent selection caused by a balance of intraspecific competition and frequency‐independent stabilizing selection. This trait also acts as mating character for assortment. It has been established previously that speciation can occur only if competition is strong enough to induce disruptive selection. We find that speciation becomes more difficult for very strong competition, because then extremely strong assortment is required. Thus, speciation is most likely for intermediate strengths of competition, where it requires strong, but not extremely strong, assortment. For this range of parameters, however, it is not obvious how assortment can evolve from low to high levels, because with moderately strong assortment less genetic variation is maintained than under weak or strong assortment sometimes none at all. In addition to the strength of frequency‐dependent competition and assortative mating, the roles of the number of loci, the distribution of allelic effects, the initial conditions, costs to being choosy, the strength of stabilizing selection, and the particular choice of the fitness function are explored. A multitude of possible evolutionary outcomes is observed, including loss of all genetic variation, splitting in two to five species, as well as very short and extremely long stable limit cycles. On the methodological side, we propose quantitative measures for deciding whether a given distribution reflects two (or more) reproductively isolated clusters.     

Additive genetic variation under intraspecific competition and stabilizing selection: a two-locus study

A diallelic two-locus model is investigated in which the loci determine the genotypic value of a quantitative trait additively. Fitness has two components: stabilizing selection on the trait and a frequency-dependent component, as induced, for instance, if the ability to utilize different food resources depends on this trait. Since intraspecific competition induces disruptive selection, this model leads to a conflict of selective forces. We study how the underlying genetics (recombination rate and allelic effects) interacts with the selective forces, and explore the resulting equilibrium structure. For the special case of equal effects, global stability results are proved. Unless the locus effects are sufficiently different, the genetic variance maintained at equilibrium displays a threshold-like dependence on the strength of competition. For loci with equal effects, the equilibrium fitnesses of genotypic values exhibit disruptive selection if and only if competition is strong enough to maintain a stable two-locus polymorphism. For unequal effects, disruptive selection can be observed for weaker competition and in the absence of a stable polymorphism.     

On a genetic model of intraspecific competition and stabilizing selection

A genetic model is investigated in which two recombining loci determine the genotypic value of a quantitative trait additively. Two opposing evolutionary forces are assumed to act: stabilizing selection on the trait, which favors genotypes with an intermediate phenotype, and intraspecific competition mediated by that trait, which favors genotypes whose effect on the trait deviates most from that of the prevailing genotypes. Accordingly, fitnesses of genotypes have a frequency-independent component describing stabilizing selection and a frequency- and density-dependent component modeling competition. We study how the underlying genetics, in particular recombination rate and relative magnitude of allelic effects, interact with the conflicting selective forces and derive the resulting, surprisingly complex equilibrium patterns. We also investigate the conditions under which disruptive selection on the phenotypes can be observed and examine how much genetic variation can be maintained in such a model. We discovered a number of unexpected phenomena. For instance, we found that with little recombination the degree of stably maintained polymorphism and the equilibrium genetic variance can decrease as the strength of competition increases relative to the strength of stabilizing selection. In addition, we found that mean fitness at the stable equilibria is usually much lower than the maximum possible mean fitness and often even lower than the fitness at other, unstable equilibria. Thus, the evolutionary dynamics in this system are almost always nonadaptive.     

Intraspecific competitive divergence and convergence under assortative mating

Ecologically driven sympatric speciation has received much attention recently. We investigate a multilocus model of a quantitative trait that is under frequency-dependent selection caused by intraspecific competition and acts as mating character for assortment. We identify the conditions that lead to the establishment of reproductively isolated clusters. This may be interpreted as evolutionary splitting or sympatric speciation. In our model, there are parameters that independently determine the strength of assortment, the costs for being choosy, and the strength of frequency-dependent natural selection. Sufficiently strong frequency dependence leads to disruptive selection on the phenotypes. The population consists of (sexual) haploid individuals. If frequency dependence is strong enough to induce disruptive selection and costs are absent or low, the result of evolution depends in a distinctive nonlinear way on the strength of assortment: under moderately strong assortment, less genetic variation is maintained than under weak or strong assortment, and sometimes there is none at all. Evolutionary splitting occurs only if frequency dependence and assortment are both strong enough and costs are low. Even then, the evolutionary outcome depends on the genetics and the initial conditions. The roles of the number of loci, of linkage, and of asymmetric selection are also explored.     

The two-locus model of Gaussian stabilizing selection

We study the equilibrium structure of a well-known two-locus model in which two diallelic loci contribute additively to a quantitative trait that is under Gaussian stabilizing selection. The population is assumed to be infinitely large, randomly mating, and having discrete generations. The two loci may have arbitrary effects on the trait, the strength of selection and the recombination rate may also be arbitrary. We find that 16 different equilibrium patterns exist, having up to 11 equilibria; up to seven interior equilibria may coexist, and up to four interior equilibria, three in negative and one in positive linkage disequilibrium, may be simultaneously stable. Also, two monomorphic and two fully polymorphic equilibria may be simultaneously stable. Therefore, the result of evolution may be highly sensitive to perturbations in the initial conditions or in the underlying genetic parameters. For the special case of equal effects, global stability results are proved. In the general case, we rely in part on numerical computations. The results are compared with previous analyses of the special case of extremely strong selection, of an approximate model that assumes linkage equilibrium, and of the much simpler quadratic optimum model.     

Optimization under frequency-dependent selection

We consider a model of frequency-dependent selection, which we refer to as the Wildcard Model. A variety of more specific models, representing quite diverse biological situations, are covered by the Wildcard Model as particular cases. Two very different particular models that are subsumed by the Wildcard Model are the game theoretically motivated two-phenotype model of Lessard [Lessard, S.,1984. Evolutionary dynamics in frequency-dependent two-phenotype models, Theor. Popul. Biol. 25, 210-234], and the model of selection on a continuous trait due to intraspecific competition of Bürger [Bürger, R., 2005. A multilocus analysis of intraspecific competition and stabilizing selection on a quantitative trait. J. Math. Biol. 50 (4), 355-396] and Schneider [Schneider, K.A., 2006. A multilocus-multiallele analysis of frequency-dependent selection induced by intraspecific competition. J. Math. Biol. 52 (4), 483-523]. Both these models have been shown in the past to have a global Lyapunov function (LF) under appropriate genetic assumptions. We show that (i) the Wildcard Model in continuous time for a single multiallelic locus, or for multiple multiallelic loci in linkage equilibrium, has a global LF, of which the Lessard and Bürger-Scheneider LF are special cases in spite of their widely different biological interpretations; (ii) the LF of the Wildcard Model can be derived from an LF previously identified for a model of density- and frequency-dependent selection due to Lotka-Volterra competition, with one locus, multiple alleles, multiple species and continuous-time dynamics [Matessi, C., Jayakar, S.D., 1981. Coevolution of species in competition: A theoretical study. Proc. Natl. Acad. Sci. USA, 78 (2, part2), 1081-1084]. We extend the LF with density and frequency dependence to the multilocus case with linkage-equilibrium dynamics. As a possible application of our results, the optimization principle we established can be used as a tool in the study of long-term evolution of various models subsumed by the Wildcard Model based on explicit short-term dynamics.     

The evolution of genetic architecture under frequency-dependent disruptive selection

We propose a model to analyze a quantitative trait under frequency-dependent disruptive selection. Selection on the trait is a combination of stabilizing selection and intraspecific competition, where competition is maximal between individuals with equal phenotypes. In addition, there is a density-dependent component induced by population regulation. The trait is determined additively by a number of biallelic loci, which can have different effects on the trait value. In contrast to most previous models, we assume that the allelic effects at the loci can evolve due to epistatic interactions with the genetic background. Using a modifier approach, we derive analytical results under the assumption of weak selection and constant population size, and we investigate the full model by numerical simulations. We find that frequency-dependent disruptive selection favors the evolution of a highly asymmetric genetic architecture, where most of the genetic variation is concentrated on a small number of loci. We show that the evolution of genetic architecture can be understood in terms of the ecological niches created by competition. The phenotypic distribution of a population with an adapted genetic architecture closely matches this niche structure. Thus, evolution of the genetic architecture seems to be a plausible way for populations to adapt to regimes of frequency-dependent disruptive selection. As such, it should be seen as a potential evolutionary pathway to discrete polymorphisms and as a potential alternative to other evolutionary responses, such as the evolution of sexual dimorphism or assortative mating.     

Determinants of the strength of disruptive and/or divergent selection arising from resource competition

We investigate how the intensity of competition for resources affects the strength of disruptive selection on a resource acquisition trait. This is done by analyzing several consumer–resource models in which consumers use a linear array of resources. We show that disruptive selection can be diminished under both strong and weak competition, making disruptive selection a unimodal function of the strength of competition. Weak selection under strong competition arises when competition causes the extinction (for self-reproducing resources) or depletion (for abiotic resources) of the most rapidly caught resources. Weak selection under weak competition is a consequence of minimal effects of consumers on resources. The precise relationship between intensity of competition and strength of disruptive selection is sensitive to the shape of the consumer's resource utilization curve and the nature of resource growth. The most strongly unimodal competition–selection relationships result from utilization curves with long tails. Our results show that a simple comparison of the width of the resource abundance distribution and the consumer's utilization function is not sufficient to determine whether selection is disruptive. The results may explain some contradictory experimental findings regarding the effect of consumer mortality on the strength of disruptive selection.     

Evolution of dominance under frequency-dependent intraspecific competition

A population-genetic analysis is performed of a two-locus two-allele model, in which the primary locus has a major effect on a quantitative trait that is under frequency-dependent disruptive selection caused by intraspecific competition for a continuum of resources. The modifier locus determines the degree of dominance at the trait level. We establish the conditions when a modifier allele can invade and when it becomes fixed if sufficiently frequent. In general, these are not equivalent because an unstable internal equilibrium may exist and the condition for successful invasion of the modifier is more restrictive than that for eventual fixation from already high frequency. However, successful invasion implies global fixation, i.e., fixation from any initial condition. Modifiers of large effect can become fixed, and also invade, in a wider parameter range than modifiers of small effect. We also study modifiers with a direct, frequency-independent deleterious fitness effect. We show that they can invade if they induce a sufficiently high level of dominance and if disruptive selection on the ecological trait is strong enough. For deleterious modifiers, successful invasion no longer implies global fixation because they can become stuck at an intermediate frequency due to a stable internal equilibrium. Although the conditions for invasion and for fixation if sufficiently frequent are independent of the linkage relation between the two loci, the rate of spread depends strongly on it. The present study provides further support to the view that evolution of dominance may be an efficient mechanism to remove unfit heterozygotes that are maintained by balancing selection. It also demonstrates that an invasion analysis of mutants of very small effect is insufficient to obtain a full understanding of the evolutionary dynamics under frequency-dependent selection.     

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.     

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

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

The strength of the selection barrier between populations

Genetic differences among populations exposed to selection form barriers against genetic exchange by mortality among hybrids. The strength of such a selection barrier, with which one (recipient) population reacts against immigration from another (donor) population, may be measured as the cumulative mean fitness of hybrids and their descendants relative to the fitness of the recipient population. Previous work analysed a case of weak selection with pairwise epistatic interactions by assuming small genetic distance between two populations in contact. The present study allows large genetic difference between the donor and recipient populations and considers weak multilocus selection with arbitrary epistatic interactions between two or more linked loci. An approximate analytical expression for the barrier strength is obtained as an expansion in which the strength of selection plays the role of a small parameter. It is shown that allele frequencies and gametic linkage disequilibria contribute in different ways to the strength of the selection barrier.     

Models of sexual selection on a quantitative genetic trait when preference is acquired by sexual imprinting

The evolution of a quantitative genetic trait under stabilizing viability selection and sexual selection is modeled for a polygynous species in which female mating preferences are acquired by sexual imprinting on the parents and by exposure to the surviving population at large. Stabilizing viability selection acts equally on both sexes in the case of a sexually monomorphic trait and on males only in the case of a dimorphic trait. A genetically fixed sensory or perceptual bias defines the origin of the scale on which the trait is measured, and the possibility is incorporated that female preferences may deviate asymmetrically from the familiar-either toward or away from this origin. When viability selection is strong relative to sexual selection, the models predict that the mean trait value will evolve to the viability optimum. With intermediate ratios of the strength of viability to sexual selection, a stable equilibrium can occur on either side of this viability optimum, depending on the direction of asymmetry in female preferences. When viability selection is relatively weak and certain other conditions are also satisfied, runaway selection is predicted.     

An analytically tractable model for competitive speciation

Several recent models have shown that frequency-dependent disruptive selection created by intraspecific competition can lead to the evolution of assortative mating and, thus, to competitive sympatric speciation. However, since most of these results rely on limited numerical analyses, their generality has been debated. Here, we consider one of the standard models (the so-called Roughgarden model) with a simplified genetics where the selected trait is determined by a single diallelic locus. This model is sufficiently complex to maintain key properties of the general multilocus case but simple enough to allow for comprehensive analytical treatment by means of invasion fitness arguments. Depending on (1) the strength and (2) the shape of stabilizing selection, (3) the strength and (4) the shape of pairwise competition, (5) the shape of the mating function, and (6) whether assortative mating leads to sexual selection, we find five different evolutionary regimes. In one of these regimes, complete reproductive isolation can evolve through arbitrarily small steps in the strength of assortative mating. Our approach provides a mechanistic understanding of several phenomena that have been found in previous models. The results demonstrate how even in a simple model, the evolutionary outcome depends in a complex way on ecological and genetic parameters.     

Multilocus selection in subdivided populations II. Maintenance of polymorphism under weak or strong migration

The potential of maintaining multilocus polymorphism by migration-selection balance is studied. A large population of diploid individuals is distributed over finitely many demes connected by migration. Generations are discrete and nonoverlapping, selection may vary across demes, and loci are multiallelic. It is shown that if migration and recombination are strong relative to selection, then with weak or no epistasis and intermediate dominance at every locus and in every deme, arbitrarily many alleles can be maintained at arbitrarily many loci at a stable equilibrium. If migration is weak relative to selection and recombination, then with weak or no epistasis and intermediate dominance at every locus and in every deme, as many alleles as there are demes can be maintained at arbitrarily many loci at equilibrium. In both cases open sets of such parameter combinations are constructed, thus the results are robust with respect to small, but arbitrary, perturbations in the parameters. For weak migration, the number of demes is, in fact, a generic upper bound to the number of alleles that can be maintained at any locus. Thus, several scenarios are identified under which multilocus polymorphism can be maintained by migration-selection balance when this is impossible in a panmictic population.      

Gene duplications contribute to the overrepresentation of interactions between proteins of a similar age

Background

Disruptive selection has been documented in a growing number of natural populations. Yet, its prevalence within individual systems remains unclear. Furthermore, few studies have sought to identify the ecological factors that promote disruptive selection in the wild. To address these issues, we surveyed 15 populations of Mexican spadefoot toad tadpoles, Spea multiplicata, and measured the prevalence of disruptive selection acting on resource-use phenotypes. We also evaluated the relationship between the strength of disruptive selection and the intensity of intraspecific competition??an ecological agent hypothesized to be an important driver of disruptive selection.

Results

Disruptive selection was the predominant mode of quadratic selection across all populations. However, a directional component of selection favoring an extreme ecomorph??a distinctive carnivore morph??was also common. Disruptive selection was strongest in populations experiencing the most intense intraspecific competition, whereas stabilizing selection was only found in populations experiencing relatively weak intraspecific competition.

Conclusions

Disruptive selection can be common in natural populations. Intraspecific competition for resources may be a key driver of such selection.