Quantitative genetic models for the balance between migration and stabilizing selection

The evolution of a quantitative trait subject to stabilizing selection and immigration, with the immigrants deviating from the local optimum, is considered under a number of different models of the underlying genetic basis of the trait. By comparing exact predictions under the infinitesimal model obtained using numerical methods with predictions of a simplified approximate model based on ignoring linkage disequilibrium, the increase in the expressed genetic variance as a result of linkage disequilibrium generated by migration is shown to be relatively small and negligible, provided that the genetic variance relative to the squared deviation of immigrants from the local optimum is sufficiently large or selection and migration is sufficiently weak. Deviation from normality is shown to be less important by comparing predictions of the infinitesimal model with a model presupposing normality. For a more realistic symmetric model, involving a finite number of loci only, no linkage and equal effects and frequencies across loci, additional changes in the genetic variance arise as a result of changes in underlying allele frequencies. Again, provided that the genetic variance relative to the squared deviation of the immigrants from the local optimum is small, the difference between the predictions of infinitesimal and the symmetric model are small unless the number of loci is very small. However, if the genetic variance relative to the squared deviation of the immigrants from the local optimum is large, or if selection and migration are strong, both linkage disequilibrium and changes in the genetic variance as a result of changes in underlying allele frequencies become important.     

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.     

LOCAL ADAPTATION WHEN COMPETITION DEPENDS ON PHENOTYPIC SIMILARITY

Recent work incorporating demographic–genetic interactions indicates the importance of population size, gene flow, and selection in influencing local adaptation. This work typically assumes that density‐dependent survival affects individuals equally, but individuals in natural population rarely compete equally. Among‐individual differences in resource use generate stronger competition between more similar phenotypes (frequency‐dependent competition) but it remains unclear how this additional form of selection changes the interactions between population size, gene flow, and local stabilizing selection. Here, we integrate migration–selection dynamics with frequency‐dependent competition. We developed a coupled demographic‐quantitative genetic model consisting of two patches connected by dispersal and subject to local stabilizing selection and competition. Our model shows that frequency‐dependent competition slightly increases local adaptation, greatly increases genetic variance within patches, and reduces the amount that migration depresses population size, despite the increased genetic variance load. The effects of frequency‐dependence depend on the strength of divergent selection, trait heritability, and when mortality occurs in the life cycle in relation to migration and reproduction. Essentially, frequency‐dependent competition reduces the density‐dependent interactions between migrants and residents, the extent to which depends on how different and common immigrants are compared to residents. Our results add new dynamics that illustrate how competition can alter the effects of gene flow and divergent selection on local adaptation and population carrying capacities.     

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.     

Pleiotropic Models of Polygenic Variation, Stabilizing Selection, and Epistasis

We show that in polymorphic populations many polygenic traits pleiotropically related to fitness are expected to be under apparent ``stabilizing selection' independently of the real selection acting on the population. This occurs, for example, if the genetic system is at a stable polymorphic equilibrium determined by selection and the nonadditive contributions of the loci to the trait value either are absent, or are random and independent of those to fitness. Stabilizing selection is also observed if the polygenic system is at an equilibrium determined by a balance between selection and mutation (or migration) when both additive and nonadditive contributions of the loci to the trait value are random and independent of those to fitness. We also compare different viability models that can maintain genetic variability at many loci with respect to their ability to account for the strong stabilizing selection on an additive trait. Let V(m) be the genetic variance supplied by mutation (or migration) each generation, V(g) be the genotypic variance maintained in the population, and n be the number of the loci influencing fitness. We demonstrate that in mutation (migration)-selection balance models the strength of apparent stabilizing selection is order V(m)/V(g). In the overdominant model and in the symmetric viability model the strength of apparent stabilizing selection is approximately 1/(2n) that of total selection on the whole phenotype. We show that a selection system that involves pairwise additive by additive epistasis in maintaining variability can lead to a lower genetic load and genetic variance in fitness (approximately 1/(2n) times) than an equivalent selection system that involves overdominance. We show that, in the epistatic model, the apparent stabilizing selection on an additive trait can be as strong as the total selection on the whole phenotype.     

A multilocus analysis of intraspecific competition and stabilizing selection on a quantitative trait

The equilibrium structure of an additive, diallelic multilocus model of a quantitative trait under frequency- and density-dependent selection is derived. The trait is under stabilizing selection and mediates intraspecific competition as induced, for instance, by differential resource utilization. It is assumed that stabilizing selection is weak, but the strength of competition may be arbitrary relative to it. Density dependence is caused by population regulation, which may be of a very general kind. The number and effects of loci are arbitrary, and stabilizing selection is not necessarily symmetric with respect to the range of phenotypic values. All previously studied models of intraspecific competition for a continuum of resources known to the author reduce to a special case of the present model if overall selection is weak. Therefore, in this case our results are applicable as approximations to all these models. Our central result is the (nearly) complete characterization of the equilibrium and stability structure in terms of all parameters. It is derived under the sole assumption that selection is weak enough relative to recombination to ignore linkage disequilibrium. In particular, necessary and sufficient conditions on the strength of competition relative to stabilizing selection are found that ensure the maintenance of multilocus polymorphism and the occurrence of disruptive selection. In this case, explicit formulas for the number of polymorphic loci at equilibrium, the allele frequencies, the genetic variance, and the strength of disruptive selection are obtained. For two loci, the effects of linkage are investigated analytically; for several loci, they are studied numerically.     

Maintenance of Genetic Variability under Strong Stabilizing Selection: A Two-Locus Model

We study a two locus model with additive contributions to the phenotype to explore the relationship between stabilizing selection and recombination. We show that if the double heterozygote has the optimum phenotype and the contributions of the loci to the trait are different, then any symmetric stabilizing selection fitness function can maintain genetic variability provided selection is sufficiently strong relative to linkage. We present results of a detailed analysis of the quadratic fitness function which show that selection need not be extremely strong relative to recombination for the polymorphic equilibria to be stable. At these polymorphic equilibria the mean value of the trait, in general, is not equal to the optimum phenotype, there exists a large level of negative linkage disequilibrium which ``hides' additive genetic variance, and different equilibria can be stable simultaneously. We analyze dependence of different characteristics of these equilibria on the location of optimum phenotype, on the difference in allelic effect, and on the strength of selection relative to recombination. Our overall result that stabilizing selection does not necessarily eliminate genetic variability is compatible with some experimental results where the lines subject to strong stabilizing selection did not have significant reductions in genetic variability.     

The Selection Limit Due to the Conflict between Truncation and Stabilizing Selection with Mutation

Long-term selection response could slow down from a decline in genetic variance or in selection differential or both. A model of conflict between truncation and stabilizing selection in infinite population size is analysed in terms of the reduction in selection differential. Under the assumption of a normal phenotypic distribution, the limit to selection is found to be a function of kappa, the intensity of truncation selection, omega 2, a measure of the intensity of stabilizing selection, and sigma 2, the phenotypic variance of the character. The maintenance of genetic variation at this limit is also analyzed in terms of mutation-selection balance by the use of the "House-of-cards" approximation. It is found that truncation selection can substantially reduce the equilibrium genetic variance below that when only stabilizing selection is acting, and the proportional reduction in variance is greatest when the selection is very weak. When truncation selection is strong, any further increase in the strength of selection has little further influence on the variance. It appears that this mutation-selection balance is insufficient to account for the high levels of genetic variation observed in many long-term selection experiments.     

Quantitative genetic variability maintained by mutation-stabilizing selection balance: sampling variation and response to subsequent directional selection

A model of genetic variation of a quantitative character subject to the simultaneous effects of mutation, selection and drift is investigated. Predictions are obtained for the variance of the genetic variance among independent lines at equilibrium with stabilizing selection. These indicate that the coefficient of variation of the genetic variance among lines is relatively insensitive to the strength of stabilizing selection on the character. The effects on the genetic variance of a change of mode of selection from stabilizing to directional selection are investigated. This is intended to model directional selection of a character in a sample of individuals from a natural or long-established cage population. The pattern of change of variance from directional selection is strongly influenced by the strengths of selection at individual loci in relation to effective population size before and after the change of regime. Patterns of change of variance and selection responses from Monte Carlo simulation are compared to selection responses observed in experiments. These indicate that changes in variance with directional selection are not very different from those due to drift alone in the experiments, and do not necessarily give information on the presence of stabilizing selection or its strength.     

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.     

Evolution at a multiallelic locus under migration and uniform selection

The semilinear parabolic system that describes the evolution of the gene frequencies in the diffusion approximation for migration and selection at a multiallelic locus is investigated. The population occupies a finite habitat of arbitrary dimensionality and shape. The drift and diffusion coefficients may depend on position, but the selection coefficients do not. It is established that if p is a uniform equilibrium point under pure selection, then p is a migration-selection equilibrium, and that generically the introduction of migration does not change the stability of p. It is also proved that if p is a uniform, globally asymptotically stable, internal equilibrium point under pure selection, then the gene frequencies converge to p when both migration and selection are present. Thus, in this case, after a sufficiently long time, there is no genetic indication of the spatial distribution of the population.     

Finite metapopulation models with density-dependent migration and stochastic local dynamics

The effects of small density-dependent migration on the dynamics of a metapopulation are studied in a model with stochastic local dynamics. We use a diffusion approximation to study how changes in the migration rate and habitat occupancy affect the rates of local colonization and extinction. If the emigration rate increases or if the immigration rate decreases with local population size, a positive expected rate of change in habitat occupancy is found for a greater range of habitat occupancies than when the migration is density-independent. In contrast, the reverse patterns of density dependence in respective emigration and immigration reduce the range of habitat occupancies where the metapopulation will be viable. This occurs because density-dependent migration strongly influences both the establishment and rescue effects in the local dynamics of metapopulations.     

Spatial Autocorrelation of Genotypes under Directional Selection

The spatial distributions of genetic variation under selection-mutation equilibrium within populations that have limited dispersal are investigated. The results show that directional selection with moderate strength rapidly reduces the amount of genetic structure and spatial autocorrelations far below that predicted for selectively neutral loci. For the latter, homozygotes are spatially clustered into separate areas or patches, each consisting of several hundred homozygotes. When selection is added the patches of the deleterious homozygotes are much smaller, in the range of 25 to 50 individuals. Selection also reduces temporal correlations. Also investigated are the effects of random replacement processes, such as mutation, immigration, and long-distance migration, on spatial and temporal correlations. The detection of natural selection through spatial pattern analysis is discussed, and applied to data from populations of the morning glory, Ipomoea purpurea.     

Deleterious Mutations, Apparent Stabilizing Selection and the Maintenance of Quantitative Variation

Apparent stabilizing selection on a quantitative trait that is not causally connected to fitness can result from the pleiotropic effects of unconditionally deleterious mutations, because as N. Barton noted, "...individuals with extreme values of the trait will tend to carry more deleterious alleles...." We use a simple model to investigate the dependence of this apparent selection on the genomic deleterious mutation rate, U; the equilibrium distribution of K, the number of deleterious mutations per genome; and the parameters describing directional selection against deleterious mutations. Unlike previous analyses, we allow for epistatic selection against deleterious alleles. For various selection functions and realistic parameter values, the distribution of K, the distribution of breeding values for a pleiotropically affected trait, and the apparent stabilizing selection function are all nearly Gaussian. The additive genetic variance for the quantitative trait is kQa2, where k is the average number of deleterious mutations per genome, Q is the proportion of deleterious mutations that affect the trait, and a2 is the variance of pleiotropic effects for individual mutations that do affect the trait. In contrast, when the trait is measured in units of its additive standard deviation, the apparent fitness function is essentially independent of Q and a2; and beta, the intensity of selection, measured as the ratio of additive genetic variance to the "variance" of the fitness curve, is very close to s = U/k, the selection coefficient against individual deleterious mutations at equilibrium. Therefore, this model predicts appreciable apparent stabilizing selection if s exceeds about 0.03, which is consistent with various data. However, the model also predicts that beta must equal Vm/VG, the ratio of new additive variance for the trait introduced each generation by mutation to the standing additive variance. Most, although not all, estimates of this ratio imply apparent stabilizing selection weaker than generally observed. A qualitative argument suggests that even when direct selection is responsible for most of the selection observed on a character, it may be essentially irrelevant to the maintenance of variation for the character by mutation-selection balance. Simple experiments can indicate the fraction of observed stabilizing selection attributable to the pleiotropic effects of deleterious mutations.     

A genetic model of interpopulation variation and covariation of quantitative characters

Evolutionary consequences of natural selection, migration, genotype-environment interaction, and random genetic drift on interpopulation variation and covariation of quantitative characters are analysed in terms of a selection model that partitions natural selection into directional and stabilizing components. Without migration, interpopulation variation and covariation depend mainly on the pattern and intensities of selection among populations and the harmonic mean of effective population sizes. Both transient and equilibrium covariance structures are formulated with suitable approximations. Migration reduces the differentiation among populations, but its effect is less with genotype-environment interaction. In some special cases of genotype-environment interaction, the equilibrium interpopulation variation and covariation is independent of migration.     

GENE FLOW FROM DOMESTICATED SPECIES TO WILD RELATIVES: MIGRATION LOAD IN A MODEL OF MULTIVARIATE SELECTION

Domesticated species frequently spread their genes into populations of wild relatives through interbreeding. The domestication process often involves artificial selection for economically desirable traits. This can lead to an indirect response in unknown correlated traits and a reduction in fitness of domesticated individuals in the wild. Previous models for the effect of gene flow from domesticated species to wild relatives have assumed that evolution occurs in one dimension. Here, I develop a quantitative genetic model for the balance between migration and multivariate stabilizing selection. Different forms of correlational selection consistent with a given observed ratio between average fitness of domesticated and wild individuals offsets the phenotypic means at migration–selection balance away from predictions based on simpler one-dimensional models. For almost all parameter values, correlational selection leads to a reduction in the migration load. For ridge selection, this reduction arises because the distance the immigrants deviates from the local optimum in effect is reduced. For realistic parameter values, however, the effect of correlational selection on the load is small, suggesting that simpler one-dimensional models may still be adequate in terms of predicting mean population fitness and viability.     

Bottlenecks in large populations: the effect of immigration on population viability

We model a large population that is subject to successive short bottlenecks, in order to investigate the impact of different extents of immigration on the change in genetic load and on viability. A first simple genetic model uncovers the opposite effects of immigration on fitness according to the type of deleterious mutations considered: immigration increases fitness if the genetic load is comprised of mildly deleterious mutations, whereas it decreases fitness if it is comprised of lethals. When considering both types of mutations and adding explicit stochastic demographic considerations, in which bottlenecks are engendered by random catastrophes, the global impact of immigration on viability is dependent upon a balance between its opposite effects on the two components of the genetic load and on demographic stochasticity. In this context, immigration tends to increase the probability of extinction if occurring preferentially when population density is high, while it decreases extinction if occurring preferentially towards low-density populations.     

Geographical Variation in a Quantitative Character

A model for the evolution of the local averages of a quantitative character under migration, selection, and random genetic drift in a subdivided population is formulated and investigated. Generations are discrete and nonoverlapping; the monoecious, diploid population mates at random in each deme. All three evolutionary forces are weak, but the migration pattern and the local population numbers are otherwise arbitrary. The character is determined by purely additive gene action and a stochastically independent environment; its distribution is Gaussian with a constant variance; and it is under Gaussian stabilizing selection with the same parameters in every deme. Linkage disequilibrium is neglected. Most of the results concern the covariances of the local averages. For a finite number of demes, explicit formulas are derived for (i) the asymptotic rate and pattern of convergence to equilibrium, (ii) the variance of a suitably weighted average of the local averages, and (iii) the equilibrium covariances when selection and random drift are much weaker than migration. Essentially complete analyses of equilibrium and convergence are presented for random outbreeding and site homing, the Levene and island models, the circular habitat and the unbounded linear stepping-stone model in the diffusion approximation, and the exact unbounded stepping-stone model in one and two dimensions.     

Environmental variation, fluctuating selection and genetic drift in subdivided populations

Although there have many studies of the population genetical consequences of environmental variation, little is known about the combined effects of genetic drift and fluctuating selection in structured populations. Here we use diffusion theory to investigate the effects of temporally and spatially varying selection on a population of haploid individuals subdivided into a large number of demes. Using a perturbation method for processes with multiple time scales, we show that as the number of demes tends to infinity, the overall frequency converges to a diffusion process that is also the diffusion approximation for a finite, panmictic population subject to temporally fluctuating selection. We find that the coefficients of this process have a complicated dependence on deme size and migration rate, and that changes in these demographic parameters can determine both the balance between the dispersive and stabilizing effects of environmental variation and whether selection favors alleles with lower or higher fitness variance.     

Isolation by Distance in a Quantitative Trait

Random genetic drift in a quantitative character is modeled for a population with a continuous spatial distribution in an infinite habitat of one or two dimensions. The analysis extends Wright's concept of neighborhood size to spatially autocorrelated sampling variation in the expected phenotype at different locations. Weak stabilizing selection is assumed to operate toward the same optimum phenotype in every locality, and the distribution of dispersal distances from parent to offspring is a (radially) symmetric function. The equilibrium pattern of geographic variation in the expected local phenotype depends on the neighborhood size, the genetic variance within neighborhoods, and the strength of selection, but is nearly independent of the form of the dispersal function. With all else equal, geographic variance is smaller in a two-dimensional habitat than in one dimension, and the covariance between expected local phenotypes decreases more rapidly with the distance separating them in two dimensions than in one. The equilibrium geographic variance is less than the phenotypic variance within localities, unless the neighborhood size is small and selection is extremely weak, especially in two dimensions. Nevertheless, dispersal of geographic variance created by random genetic drift is an important mechanism maintaining genetic variance within local populations. For a Gaussian dispersal function it is shown that, even with a small neighborhood size, a population in a two-dimensional habitat can maintain within neighborhoods most of the genetic variance that would occur in an infinite panmictic population.