*


Consider an allele couple A_{1} et A_{2} of frequencies p_{t} and q_{t} respectively in generation t; consider u as the direct mutation rate of A_{1} towards A_{2} in each generation and v the reverse mutation rate of A_{2} towards A_{1}.
We define the mutation rate as the probability for a mutation to appear per gamete and per generation.
For example, suppose a population only composed of individuals of the genotype A_{1} A_{1} which contribute to the production of 2N gametes in the next generation.
If u represents the mutation rate of A_{1} locus towards all the other possible alleles in a diploid population of N individuals, as we have 2N concerned genes, the number of new mutants that we get are 2N x u, in a given generation. In the population, on 2N genes, only one gene among them (1/2N) will be fixed. Thus, the probability for a mutant allele to appear in the generation with a fixed gene, is 2N x u x 1/2N = u.
Therefore, the probability for one allele to be substitute by a muted allele is equal to the mutation rate, per gamete and per generation, in a given locus.
II Recurrence law
In the presence of mutations, the evolution of allele frequencies depends on p_{i} and q_{i} frequencies in generation i, but also on the mutation rates of u and v.
The frequency of an allele A_{1}, one generation after mutation (p_{i} + 1), corresponds to the frequency of an allele A_{1} in the preceding generation (frequency p_{i}) which would not have muted (probability equal to 1u), or to an allele A_{2} (frequency q_{i}) which would have muted in A_{1} (probability v).
p_{i + 1} = (1u) p_{i} + v q_{i}
q_{i + 1} = (1v) q_{i} + u p_{i}
II1. Example: The peppered moth (with the kind authorisation of Pr. Georges Périquet, University of Tours)
The peppered moth, Biston betularia (Order: Lepidoptera) is common in Northern Europe. The individuals fly during the night and in the day lay down against the lightcoloured bark of trees. Since the 19th century, this species has been often studied; it is represented by two morphs, one lightcoloured (typica) and another darkcoloured, or melanic (carbonaria).
The genetic determinism of this coloration is monogenetic and autosomic. The allele carbonaria (C) is dominant on the allele typica (c). The allele frequencies C and c will be named p and q respectively.II1.1. Genetic determinism and allele frequencies
Phenotypes  [Melanic form]  [Light form] 
Genotypes  CC ou Cc  cc 
In the middle of the 19^{th} century, among the British populations, the typica form was largely majority. It is only in 1848 that an individual carbonaria capture had been reported in Manchester area. This form frequency had strongly increased and other melanic individuals, afterwards, had been captured in others industrial areas of England. The frequency increase of the carbonaria form was very rapid in Manchester region; it reached 98 %, in 1895. In less than 50 years, this form had become the great majority into the region. As the variations always occurred in the way of a carbonaria form increase of frequency, it could not be due to a random phenomenon.
Two hypothesis are possible; the effect of recurrent mutations, or the action by selection.II1.2. Two hypothesis
The hypothesis of recurrent mutations:
We could think that the carbonaria allele was obtained frequently enough by mutation and thus replace the typica allele. Therefore, we can calculate the mutation rate v from c towards C, which would corresponds to an evolution as rapid as the one observed in the nature.
Under mutations effect alone, the recurrence relation for the allele frequencies c, between two successive generations, can be written as:
q_{1} = q_{0}  vq_{0}
The c allele frequency will be the same as the preceding generation (q_{0}) reduced by the muted allele frequency (vq_{0}). Here, We consider that the allele C does not mute in c, which correspond to the most rapid possible evolution under the mutation effect.
As the process is repeated with the course of the generations, we deduce:
q_{t} = q _{0} (1v)^{n}Where n represents the generation number.and
v = 1  (q_{t}/q_{0})^{1/n}
If we consider that the melanic form frequency, in the Manchester region, was 0,01 in 1848 and 0.98 in 1895, and knowing that the peppered moth gives one generation per year in this region, it is possible, therefore, to calculate the v value which better express the rapidity of the observed evolution.
Thus, in 47 generations, we obtain:
v = 1  (0,141/0,995)^{1/47} = 0,041.
Then, the calculate value obtained was about 4 %. That means, to obtain a rapid evolution, as fast as what we can observe in the nature, about 4 gametes per 100 at each generation should be muted from typica to carbonaria. In fact, it is an incompatible rate mutation in respect with the spontaneous mutation rate observed in reality (about 10^{5} to 10^{6}). Furthermore, we have to explain the brutal change of mutation rate since 1848!
If the rate of 4 5.10^{2} really exists, it could be explained by the existence of Transposons (doubtful) or by the effect of other factors of evolution such as the selection.
III Equilibrium state
The equilibrium state in which performed direct and reverse mutations.
q _{1 }= (1v) q_{0} + u p_{0}If we pose at this level of equilibrium stateq_{1 }= q_{0 } v q_{0} + u p_{0}
Δ_{q} = the difference between 2 successive generations
Δ_{q} = q_{n+1}  q_{ n}
Following the sign Δ_{ q} we know if:
Δ_{q} >0 > q_{n+1} > q _{n} > A_{2} increase
Δ_{q} < 0 > q_{n+1} < q_{n} > A_{2} decrease
Δ_{q } = 0 > q_{n+1} = q_{n } = q_{e} Equilibrium state
q_{e} = q_{e}  v q_{e} + u p_{e}The equilibrium state is function of the mutation rate u and v.u p_{e} = v q_{e}
u p_{e }= v ( 1 – p_{e})
u p_{e} + v p_{e} = v
p_{e} = v / (u + v)
and
q_{e} = u / (u + v)
q_{e} ε] 0 , 1 [
Example:
If u = v > q_{e} = 1/2If u = 10^{5} and v = 10^{6} > q_{e} = 0.91
What ever the mutation rates are, we can have any equilibrium state.
To pass of q_{0} = 0.51 to q_{n} = 0.71 and q_{e} = 0. 91
Mutations can produce equilibrium states, but they will be reached very slowly through the time. The mutation process does not have major effect on the genetic structure of populations; the variation of allele frequencies is very low through the time. This evolution of allele frequencies varies according the equation:
p_{n+1} = p_{e} + ( p_{0}  p_{e} ) ( 1uv)^{n}
q_{n+1} = q_{e} + ( q_{0}  q_{e} ) ( 1uv)^{n}
Study by simulation:
The proposed program permits to simulate the evolution of the frequency q through the generations, from q_{0} to q_{e}, for different values of mutation rate parameters u and v and of q_{0}.
Translation: Maureen Labarussias
Contributor(s) 
Written  200503  R Kalmes 
Institut de Recherche sur la Biologie de l'Insecte, IRBI  CNRS  ESA 6035, Av. Monge, F37200 Tours, France 
© Atlas of Genetics and Cytogenetics in Oncology and Haematology  indexed on : Fri Jun 30 11:25:43 CEST 2017 
For comments and suggestions or contributions, please contact us