Consider an allele couple A1 et A2 of frequencies pt and qt respectively in generation t; consider u as the direct mutation rate of A1 towards A2 in each generation and v the reverse mutation rate of A2 towards A1.
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 A1 A1 which contribute to the production of 2N gametes in the next generation.
If u represents the mutation rate of A1 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 pi and qi frequencies in generation i, but also on the mutation rates of u and v.
The frequency of an allele A1, one generation after mutation (pi + 1), corresponds to the frequency of an allele A1 in the preceding generation (frequency pi) which would not have muted (probability equal to 1-u), or to an allele A2 (frequency qi) which would have muted in A1 (probability v).
pi + 1 = (1-u) pi + v qi
qi + 1 = (1-v) qi + u pi
II-1. 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 light-coloured bark of trees. Since the 19th century, this species has been often studied; it is represented by two morphs, one light-coloured (typica) and another dark-coloured, 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.
II-1.1. Genetic determinism and allele frequencies
|Phenotypes||[Melanic form]||[Light form]|
|Genotypes||CC ou Cc||cc|
In the middle of the 19th 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.
II-1.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:
q1 = q0 - vq0
The c allele frequency will be the same as the preceding generation (q0) reduced by the muted allele frequency (vq0). 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:
qt = q 0 (1-v)nWhere n represents the generation number.
v = 1 - (qt/q0)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 = (1-v) q0 + u p0If we pose at this level of equilibrium state
q1 = q0 - v q0 + u p0
Δq = the difference between 2 successive generations
Δq = qn+1 - q n
Following the sign Δ q we know if:
Δq >0 ----> qn+1 > q n ----> A2 increase
Δq < 0 ----> qn+1 < qn ----> A2 decrease
Δq = 0 ----> qn+1 = qn = qe Equilibrium state
qe = qe - v qe + u peThe equilibrium state is function of the mutation rate u and v.
u pe = v qe
u pe = v ( 1 – pe)
u pe + v pe = v
pe = v / (u + v)
qe = u / (u + v)
qe ε] 0 , 1 [
If u = v ---> qe = 1/2
If u = 10-5 and v = 10-6 ---> qe = 0.91
What ever the mutation rates are, we can have any equilibrium state.
To pass of q0 = 0.51 to qn = 0.71 and qe = 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:
pn+1 = pe + ( p0 - pe ) ( 1-u-v)n
qn+1 = qe + ( q0 - qe ) ( 1-u-v)n
Study by simulation:
The proposed program permits to simulate the evolution of the frequency q through the generations, from q0 to qe, for different values of mutation rate parameters u and v and of q0.
Translation: Maureen Labarussias
|Institut de Recherche sur la Biologie de l'Insecte, IRBI - CNRS - ESA 6035, Av. Monge, F-37200 Tours, France|
|© Atlas of Genetics and Cytogenetics in Oncology and Haematology||indexed on : Tue Mar 14 13:57:40 CET 2017|
For comments and suggestions or contributions, please contact us