Effects of the shape of distribution of mutant effect in nearly neutral mutation models

Models of the theory of nearly neutral mutation incorporate a continuous distribution of mutation effects in contrast to the theory of purely neutral mutation which allows no mutations with intermediate effects. Previous studies of one such model, namely the house-of-cards mutation model, assumed normal distribution of mutation effect. Here I study the house-of-cards mutation model in random-mating finite populations using the weak-mutation approximation, paying attention to the effects of the distribution of mutant effects. The average selection coefficient, substitution rate and average heterozygosity in the equilibrium and transient states were studied mainly by computer simulation. The main findings are: (i) Very rapid decrease of the substitution rate and very slow approach to equilibrium as selection becomes stronger are characteristics of assuming normal distribution of mutant effect. If the right tail of the mutation distribution decays more rapidly than that of the normal distribution, the decrease of substitution rate becomes slower and equilibrium is achieved more quickly. (ii) The dispersion index becomes smaller or larger than 1 depending on the time and the intensity of selection, (iii) LetN be the population size. When selection is strong the ratio of 4N times the substitution rate to the average heterozygosity, which is expected to be 1 under neutrality, is larger than 1 in earlier generations but becomes less than 1 in later generations. These findings show the importance of the distribution of mutant effect and time in determination of the behaviour of various statistics frequently used in the study of molecular evolution.     

Mutation-selection equilibrium as a possible cause of an interchange chromosome polymorphism in a cultivar of rye,Secale cereale L.

Summary A persistent chromosomal polymorphism exists in a population of cultivated rye, Secale cereale (Candela et al. 1979). In order to ascertain the possible causes that maintain it, we have estimated the fitness values of structurally homozygous and heterozygous plants and the mutation rate of spontaneous interchange. The estimates of the selection coefficient against heterozygotes (s=0.15-0.40) and of the mutation rate u=6.12×10^–2 support a mutation-selection equilibrium as a likely cause of the interchange chromosome polymorphism.     

Expected Relationship Between the Silent Substitution Rate and the GC Content: Implications for the Evolution of Isochores

The relationship between the silent substitution rate (K _s) and the GC content along the genome is a focal point of the debate about the origin of the isochore structure in vertebrates. Recent estimation of the silent substitution rate showed a positive correlation between K _s and GC content, in contradiction with the predictions of both the regional mutation bias model and the selection or biased gene conversion model. The aim of this paper is to help resolve this contradiction between theoretical studies and data. We analyzed the relationship between K _s and GC content under (1) uniform mutation bias, (2) a regional mutation bias, and (3) mutation bias and selection. We report that an increase in K _s with GC content is expected under mutation bias because of either nonequilibrium of the isochore structure or an increasing mutation rate from AT toward GC nucleotides in GC-richer isochores. We show by simulations that CpG deamination tends to increase the mutation rate with GC content in a regional mutation bias model. We also demonstrate that the relationship between K _s and GC under the selectionist or biased gene conversion model is positive under weak selection if the mutation selection equilibrium GC frequency is less than 0.5. Received: 28 March 2001 / Accepted: 16 May 2001     

Linkage and the Maintenance of Heritable Variation by Mutation-Selection Balance

The role of linkage in influencing heritable variation maintained through a balance between mutation and stabilizing selection is investigated for two different models. In both cases one trait is considered and the interactions within and between loci are assumed to be additive. Contrary to most earlier investigations of this problem no a priori assumptions on the distribution of genotypic values are imposed. For a deterministic two-locus two-allele model with recombination and mutation, related to the symmetric viability model, a complete nonlinear analysis is performed. It is shown that, depending on the recombination rate, multiple stable equilibria may coexist. The equilibrium genetic and genic variances are calculated. For a polygenic trait in a finite population with a possible continuum of allelic effects a simulation study is performed. In both models the equilibrium genetic and genic variances are roughly equal to the house-of-cards prediction or its finite population counterpart as long as the recombination rate is not extremely low. However, negative linkage disequilibrium builds up. If the loci are very closely linked the equilibrium additive genetic variance is slightly lower than the house-of-cards prediction, but the genic variance is much higher. Depending on whether the parameters are in favor of the house-of-cards or the Gaussian approximation, different behavior of the genetic system occurs with respect to linkage.     

The effects of spontaneous mutation on competitive fitness in Drosophila melanogaster

Abstract We have analysed the effect of 288 generations of mutation accumulation (MA) on chromosome II competitive fitness in 21 full‐sib lines of Drosophila melanogaster and in a large control population, all derived from the same isogenic base. The rate of mean log‐fitness decline and that of increase of the between‐line variance were consistent with a low rate (λ ≈ 0.03 per gamete and generation), and moderate average fitness effect [E(s) ≈ 0.1] of deleterious mutation. Subsequently, crosses were made between pairs of MA lines, and these were maintained with effective size on the order of a few tens. In these crosses, MA recombinant chromosomes quickly recovered to about the average fitness level of control chromosomes. Thus, deleterious mutations responsible for the fitness decline were efficiently selected against in relatively small populations, confirming that their effects were larger than a few percent.     

Hardy-Weinberg equilibria in random mating populations

The structure of multiloci random mating populations is examined. Sufficient conditions for the existence of stable local Hardy-Weinberg equilibria for n loci and an arbitrary number of alleles per locus, are then derived for specified situations under the assumption of multiplicative gene action between loci. It is shown that a stable Hardy-Weinberg equilibrium can not be a local maximum of the mean fitness function with multiplicative gene action between loci. The stability of Hardy-Weinberg type border points and the condition for the increase of newly introduced genes are topics on which some n-loci results are also obtained for an arbitrary number of alleles per locus in systems that allow Hardy-Weinberg equilibria.     

Simultaneous Estimation of Additive and Mutational Genetic Variance in an Outbred Population of Drosophila serrata

How new mutations contribute to genetic variation is a key question in biology. Although the evolutionary fate of an allele is largely determined by its heterozygous effect, most estimates of mutational variance and mutational effects derive from highly inbred lines, where new mutations are present in homozygous form. In an attempt to overcome this limitation, middle-class neighborhood (MCN) experiments have been used to assess the fitness effect of new mutations in heterozygous form. However, because MCN populations harbor substantial standing genetic variance, estimates of mutational variance have not typically been available from such experiments. Here we employ a modification of the animal model to analyze data from 22 generations of Drosophila serrata bred in an MCN design. Mutational heritability, measured for eight cuticular hydrocarbons, 10 wing-shape traits, and wing size in this outbred genetic background, ranged from 0.0006 to 0.006 (with one exception), a similar range to that reported from studies employing inbred lines. Simultaneously partitioning the additive and mutational variance in the same outbred population allowed us to quantitatively test the ability of mutation-selection balance models to explain the observed levels of additive and mutational genetic variance. The Gaussian allelic approximation and house-of-cards models, which assume real stabilizing selection on single traits, both overestimated the genetic variance maintained at equilibrium, but the house-of-cards model was a closer fit to the data. This analytical approach has the potential to be broadly applied, expanding our understanding of the dynamics of genetic variance in natural populations.     

Genealogy and subpopulation differentiation under various models of population structure

 The structured coalescent is used to calculate some quantities relating to the genealogy of a pair of homologous genes and to the degree of subpopulation differentiation, under a range of models of subdivided populations and assuming the infinite alleles model of neutral mutation. The classical island and stepping-stone models of population structure are considered, as well as two less symmetric models. For each model, we calculate the Laplace transform of the distribution of the coalescence time of a pair of genes from specified locations and the corresponding mean and variance. These results are then used to calculate the values of Wright’s coefficient F _ST , its limit as the mutation rate tends to zero and the limit of its derivative with respect to the mutation rate as the mutation rate tends to zero. From this derivative it is seen that F _ST can depend strongly on the mutation rate, for example in the case of an essentially one-dimensional habitat with many subpopulations where gene flow is restricted to neighbouring subpopulations. Received: 1 October 1997 / Revised version: 15 March 1998     

The frequency of the perfect genotype in a population subject to pleiotropic mutation

We consider a large population of asexual organisms characterised by a number of quantitative traits that are subject to stabilising selection. Mutation is taken to act pleiotropically, with every mutation generally changing all of the traits under selection. We focus on the equilibrium distribution of the population, where mutation and selection are in balance. It has been previously established that the equilibrium distribution of genotypic effects may be anomalous, as it may contain a singular spike--a Dirac delta function--corresponding to a non-zero proportion of the population having exactly optimal genotypic values. In the present work, we present exact results for the case where three traits are under selection. These results give the equilibrium genetic variance of the population, and the proportion of the population that have the optimal genotype. This is achieved for two different spherically symmetric distributions of mutant effects. Additionally, a simple and robust numerical approach is also presented that allows the treatment of some other mutation distributions, where there are an arbitrary number of selected traits.     

An explicit representation of the Luria–Delbrück distribution

The probability distribution of the number of mutant cells in a growing single-cell population is presented in explicit form. We use a discrete model for mutation and population growth which in the limit of large cell numbers and small mutation rates reduces to certain classical models of the Luria-Delbrück distribution. Our results hold for arbitrarily large values of the mutation rate and for cell populations of arbitrary size. We discuss the influence of cell death on fluctuation experiments and investigate a version of our model that accounts for the possibility that both daughter cells of a non-mutant cell might be mutants. An algorithm is presented for the quick calculation of the distribution. Then, we focus on the derivation of two essentially different limit laws, the first of which applies if the population size tends to infinity while the mutation rate tends to zero such that the product of mutation rate times population size converges. The second limit law emerges after a suitable rescaling of the distribution of non-mutant cells in the population and applies if the product of mutation rate times population size tends to infinity. We discuss the distribution of mutation events for arbitrary values of the mutation rate and cell populations of arbitrary size, and, finally, consider limit laws for this distribution with respect to the behavior of the product of mutation rate times population size. Thus, the present paper substantially extends results due to Lea and Coulson (1949), Bartlett (1955), Stewart et al. (1990), and others.     

Stabilizing selection,mutational bias,and the evolution of sex*

Stabilizing selection around a fixed phenotypic optimum is expected to disfavor sexual reproduction, since asexually reproducing organisms can maintain a higher fitness at equilibrium, while sex disrupts combinations of compensatory mutations. This conclusion rests on the assumption that mutational effects on phenotypic traits are unbiased, that is, mutation does not tend to push phenotypes in any particular direction. In this article, we consider a model of stabilizing selection acting on an arbitrary number of polygenic traits coded by bialellic loci, and show that mutational bias may greatly reduce the mean fitness of asexual populations compared with sexual ones in regimes where mutations have weak to moderate fitness effects. Indeed, mutation and drift tend to push the population mean phenotype away from the optimum, this effect being enhanced by the low effective population size of asexual populations. In a second part, we present results from individual‐based simulations showing that positive rates of sex are favored when mutational bias is present, while the population evolves toward complete asexuality in the absence of bias. We also present analytical (QLE) approximations for the selective forces acting on sex in terms of the effect of sex on the mean and variance in fitness among offspring.     

Quantifying natural seasonal variation in mutation parameters with mutation accumulation lines

Mutations create novel genetic variants, but their contribution to variation in fitness and other phenotypes may depend on environmental conditions. Furthermore, natural environments may be highly heterogeneous. We assessed phenotypes associated with survival and reproductive success in over 30,000 plants representing 100 mutation accumulation lines of Arabidopsis thaliana across four temporal environments at a single field site. In each of the four assays, environmental variance was substantially larger than mutational variance. For some traits, whether mutational variance was significantly varied between seasons. The founder genotype had mean trait values near the mean of the distribution of the mutation accumulation lines in all field experiments. New mutations also contributed more phenotypic variation than would be predicted, given phenotypic and sequence‐level divergence among natural populations of A. thaliana. The combination of large environmental variance with a mean effect of mutation near zero suggests that mutations could contribute substantially to standing genetic variation.     

Interchromosomal Biased Gene Conversion, Mutation and Selection in a Multigene Family

A mathematical model of the effects of interchromosomal biased gene conversion, mutation and natural selection on a multigene family is developed and analyzed. The model assumes two allelic states at each of n loci. The effects of genetic drift are ignored. The model is developed under the assumption of no recombination, but the analysis shows that, at equilibrium, there is no linkage disequilibrium, which implies that the conclusions are valid for arbitrary recombination among loci. At equilibrium, the balance between mutation, gene conversion and selection depends on the ratio of the mutation rates to the quantity [s + g(2α - 1)/ n], where s is the increment or decrement in relative fitness with each additional copy of one of the alleles, g is the conversion rate, and α is a measure of the bias in favor of one of the alleles. When this quantity is large relative to the mutation rates, the allele that has the net advantage, combining the effects of selection and conversion, will be nearly fixed in the multigene family. A comparison of these results with those from a comparable model of intrachromosomal biased conversion shows that biased interchromosomal conversion leads to approximately the same equilibrium copy number as does intrachromosomal conversion of the same strength. Interchromosomal conversion is much more effective in causing the substitution of one allele by another. The relative frequencies of interchromosomal and intrachromosomal conversion is indicated by the extent of the linkage disequilibrium among the loci in a multigene family.     

A Study on a Nearly Neutral Mutation Model in Finite Populations

As a nearly neutral mutation model, the house-of-cards model is studied in finite populations using computer simulations. The distribution of the mutant effect is assumed to be normal. The behavior is mainly determined by the product of the population size, N, and the standard deviation, sigma, of the distribution of the mutant effect. If 4N sigma is large compared to one, a few advantageous mutants are quickly fixed in early generations. Then most mutation becomes deleterious and very slow increase of the average selection coefficient follows. It takes very long for the population to reach the equilibrium state. Substitutions of alleles occur very infrequently in the later stage. If 4N sigma is the order of one or less, the behavior is qualitatively similar to that of the strict neutral case. Gradual increase of the average selection coefficient occurs and in generations of several times the inverse of the mutation rate the population almost reaches the equilibrium state. Both advantageous and neutral (including slightly deleterious) mutations are fixed. Except in the early stage, an increase of the standard deviation of the distribution of the mutant effect decreases the average heterozygosity. The substitution rate is reduced as 4N sigma is increased. Three tests of neutrality, one using the relationship between the average and the variance of heterozygosity, another using the relationship between the average heterozygosity and the average number of substitutions and Watterson's homozygosity test are applied to the consequences of the present model. It is found that deviation from the neutral expectation becomes apparent only when 4N sigma is more than two. Also a simple approximation for the model is developed which works well when the mutation rate is very small.     

Convergence of multilocus systems under weak epistasis or weak selection

 The convergence of multilocus systems under viability selection with constant fitnesses is investigated. Generations are discrete and nonoverlapping; the monoecious population mates at random. The number of multiallelic loci, the linkage map, dominance, and epistasis are arbitrary. It is proved that if epistasis or selection is sufficiently weak (and satisfies a certain nondegeneracy assumption whose genericity we establish), then there is always convergence to some equilibrium point. In particular, cycling cannot occur. The behavior of the mean fitness and some other aspects of the dynamics are also analyzed. Received: 15 November 1997 / Revised version: 25 May 1998     

The Quantitative Genetic Consequences of Pleiotropy under Stabilizing and Directional Selection

The independence of two phenotypic characters affected by both pleiotropic and nonpleiotropic mutations is investigated using a generalization of M. Slatkin's stepwise mutation model of 1987. The model is used to determine whether predictions of either the multivariate normal model introduced in 1980 by R. Lande or the house-of-cards model introduced in 1985 by M. Turelli can be regarded as typical of models that are intermediate between them. We found that, under stabilizing selection, the variance of one character at equilibrium may depend on the strength of stabilizing selection on the other character (as in the house-of-cards model) or not (as in the multivariate normal model) depending on the types of mutations that can occur. Similarly, under directional selection, the genetic covariance between two characters may increase substantially (as in the house-of-cards model) or not (as in the multivariate normal model) depending on the kinds of mutations that are assumed to occur. Hence, even for the simple model we consider, neither the house-of-cards nor the multivariate normal model can be used to make predictions, making it unlikely that either could be used to draw general conclusions about more complex and realistic models.     

Adaptive Topography in Fertility-Viability Selection Models: The Haplodiploid Case

A linear combination of partial changes of mean fitnesses from one generation to the next one is shown to be approximately equal to the additive genetic variance in fitness after enough generations and away from equilibrium in random mating haplodiploid populations under arbitrary weak frequency-dependent selection on sex-differentiated viability of individuals and sex-differentiated fertility of matings controlled at a single multiallelic locus. The result can be applied to X-linked locus models in diploid populations. The result is used to deduce approximate adaptive topographies far frequency-independent selection models in the cases of nonsex-differentiated fertilities and multiplicative sex-differentiated fertilities and for kin selection models in family-structured populations under the assumptions of single insemination and multiple insemination of females. Multiple insemination creates frequency-dependent selection regimes.     

The Evolution of One- and Two-Locus Systems

Assuming age-independent fertilities and mortalities and random mating, continuous-time models for a monoecious population are investigated for weak selection. A single locus with multiple alleles and two alleles at each of two loci are considered. A slow-selection analysis of diallelic and multiallelic two-locus models with discrete nonoverlapping generations is also presented. The selective differences may be functions of genotypic frequencies, but their rate of change due to their explicit dependence on time (if any) must be at most of the second order in s, (i.e., O( s²)), where s is the intensity of natural selection. Then, after several generations have elapsed, in the continuous time models the time-derivative of the deviations from Hardy-Weinberg proportions is of O(s²), and in the two-locus models the rate of change of the linkage disequilibrium is of O(s²). It follows that, if the rate of change of the genotypic fitnesses is smaller than second order in s (i.e., o(s²)), then to O(s²) the rate of change of the mean fitness of the population is equal to the genic variance. For a fixed value of s, however, no matter how small, the genic variance may occasionally be smaller in absolute value than the (possibly negative) lower order terms in the change in fitness, and hence the mean fitness may decrease. This happens if the allelic frequencies are changing extremely slowly, and hence occurs often very close to equilibrium. Some new expressions are derived for the change in mean fitness. It is shown that, with an error of O( s), the genotypic frequencies evolve as if the population were in Hardy-Weinberg proportions and linkage equilibrium. Thus, at least for the deterministic behavior of one and two loci, deviations from random combination appear to have very little evolutionary significance.     

Asymptotic method for peristaltic transport

The purpose of this paper is to justify an asymptotic method developed for the study of peristaltic transport in a tube of arbitrary cross section. Within the framework of long wave approximation, the three-dimensional nonlinear Navier-Stokes equations are reduced to a sequence of two-dimensional linear boundary value problems of Laplace and biharmonic operators. It is shown that, if a Reynolds number is less than some constant, the solution of the approximate equations is indeed an asymptotic approximation to the exact solution of the problem as the ratio of the maximum radius of the tube to the wave length of the peristaltic motion of the wall tends to zero, and the error estimates are expressed inL ₂ norms. Furthermore, under the same condition the exact solution is shown to be unique and stable under arbitrary perturbation of spatially periodic disturbance. Application of the stability condition to peristaltic transport in a tube of circular cross section is given.