The authors have declared that no competing interests exist.

Processes of genetic-mathematical modeling of a population development are considered. A basic distinction in the mathematical description of a family tree and a population is shown. In a family tree alternation of generations has discrete character. In a population there is a continuous alternation of generations. The method of the differential equations is applied for the description of a population. It is shown that mutational process in a population can be described with use of a Green’s function. For radiating influence on a population the universal evolutionary law is found.

In the nature there are two opposite directed processes: preservation of species and evolutionary changes in these species.

The primary theory of evolution has been developed by Ch. Darwin. He considered driving forces of evolution a hereditary variability and natural selection. He also has shown that various breeds of pets and cultural plants are created by the man during artificial selection. However Darwin’s theory a question on essence and mechanisms of hereditary variability has left without the answer. At artificial selection useful properties of organisms spontaneously without intervention of the person are not kept and organisms gradually pass in a wild status.

In the beginning of 20-th century the mutational theory of evolution has appeared. This theory asserts that mutational forms arise abruptly in consequence there are new species

Now be of the opinion the synthetic theory of evolution which includes as elements Darwin’s, and mutational theories of evolution.

The population will consist of family trees set which periodically contact among themselves. On

For an example three family trees are shown. The dotted line shows possible ways of crossing of male and female individuals in family trees. If for a family tree it is possible to consider the time of one alternation generation approximately T≈25 - 30 years for a population this time can be any _{1}_{2}, etc. For a population continuous alternation of generations is characteristic. Therefore a family tree and a population from the point of view of mathematical genetics laws should be analyzed differently.

The purpose of given article is attempt of mathematical modeling of mutational processes evolution in a population.

Before to proceed to studying influence of the mutational factor to a population we shall consider some mathematical questions of genetic inheritance in case of this factor absence influence.

For a family tree the analysis we shall carry out on the basis of Hardy – Weinberg law ^{2} so

According to this law genotypes

(^{2}: (^{2} …...(1)

The Hardy - Weinberg equilibrium state is indifferent. On the autosomal inheritance it is obvious. Really, using distribution of genotypes (1) it is possible to receive, for example, the frequency of a recessive allele

In the following (

For the linked to a sex genome the analysis is necessary more complex.

Let's assume that alleles _{m} and in the women _{f} . For recessive allele _{m} and _{f}.

At crossing in the first generation there is a ratio of genotypes in the women according to product (_{f} + _{f}) (_{m} +_{m}):

Using the distribution of genotypes (3) it is possible to find a frequency of allele

…..(4)

In the deduction (4) the following obvious ratios _{mn} = 1-_{mn}, and _{fn} = 1-_{fn} are used.

Thus, the frequency of allele

For convenience of the formula (5) further analysis we shall write down it with displacement on one generation back:

At absence of mutagen effect the frequency of allele_{m}_{(n-1)} = _{f}_{(n-2)}. Using the given condition, from (6) we shall find:

The solution of the finite-differential equation (7) we search as _{fn}_{ = }^{n} where in this caseb is constant. Substituting this solution in the formula (7) we have:

Let's divide the equation (8) on ^{n-2}:

b - 2b^{2} + 1= 0 ……..(9)

We find two roots of the characteristic quadratic (9):

_{1}=1, and _{2}=1/2, …..(10)

Hence the general solution of the finite-differential equation (7) looks like:

_{fn}_{1} + C_{2}(-1/2)^{n} …….(11)

The formula (11) has been found by a different way in

Constants of integration c_{1} also c_{2} we shall find on the basis of the initial conditions: at n=0 , _{fn}=_{f0}, and at n = 1 according to (4) _{fn} = _{f1 }= _{mo} + _{f0} / 2. Thus:

C_{1} = 2q_{f0} + _{m0} /3, and C_{2} = _{f0} + _{m0} / 3. ….(12)

Therefore the solution (11) finally looks like:

.......(13)

On

The initial conditions for calculation are accepted _{f0} = 0.6, and _{m0} = 0.2. The _{fn} = 2_{f0} + _{m0} / 3 = 0.467.

Hardy - Weinberg law in the form before considered concerns to a separate family tree. Implicitly this law includes time because alternation of generations occurs through certain time. Thus Hardy - Weinberg law has the discrete character on time. The population will consist of family trees crossed among themselves and lives in continuous time. Alternation of generations of the family trees set results to generations are replaced actually according to a continuous time scale.

In it the basic distinction of the Hardy - Weinberg law for a population and a family tree. In earlier used form (1) Hardy - Weinberg law is written down for a family tree.

Let's pass to a continuous scale of the generations

Let's find out, whether there is the differential equation having the characteristic equation similar (8).

For this purpose we shall consider the differential equation:

…….(14)

where Ƞ is constant. In this formula _{f}_{(n-1)} there is recessive allele frequency at the woman in n - 1 generation. Differentiation in (14) goes on dimensionless time

Let's transform the equation (14) to finite-differential form:

where size Δ

Uniting similar members and multiplying the equation (15) on - 2 we shall find:

Let's try to identify the equations (16) and (7) for alternation of the recessive allele frequencies in generations of a family tree. For this purpose it is necessary to accept:

2(2 - ȠΔ

Wonderful feature of the algebraic equations (17) is that they have the identical solution:

ȠΔ

It’s means that the finite-differential equation (7) and the differential equation (14) can have the same characteristic equation. Taking into account (18) the equation (14) can be copied as:

The equation (19) can be integrated once:

where _{1} is constant of integration.

Further integrating the equation (20) by method of a separation of variables:

…….(21)

Where _{fo} there is initial frequency of examined recessive allele at the woman we shall find ^{fn}^{}= q_{f}(n) - the time dependence of it allele on frequency:

Identifying the solution (22) with the solution (13) we shall find:

and

……..(23)

Where _{m0} there is initial frequency of examined recessive allele at the man.

As it was expected the formulas (23) do not contradict each other. Hence, the solution (22) can be written down as:

The formula (24) is correct for frequency of female allele

Comparing (13) and (24) for even generations we have 2^{-n}=e-3/2Δ^{n}or Δ

…...(25)

that is identical to the formula (13) for even generations

Taking into account (18), and Δ

........(26)

As generations in a population are replaced continuously in the formula (26) there is no necessity to specify number of generation. The formula (26) represents Hardy - Weinberg law for a population. We shall note that Hardy - Weinberg law in form (26) as well as (1) reflects indifferent balance of a genome since to it satisfies the solution _{f} = const.

On

Let's consider influence of some mutational factor on a population. It can be abrupt change of the environment characteristics: increase of ultra-violet or radioactive radiation, emission of harmful products of a natural (volcanic) or artificial origin, distribution of pesticides, etc. All these factors change of the alleles frequency in a population, i.e. cause genic modification.

It is possible to result as an example decrease in frequency of the recessive alleles of the genes determining the concentration of a catalase in blood in 15.5 times approximately for 2 generations for people living in Hiroshima and Nagasaki in comparison with other people of Japan. Herewith the frequency of the dominant alleles practically has not changed

The characteristic of a system response on some influence generally is Green's function of this system

Let's consider the differential equation (26) reflecting the Hardy – Weinberg law for female allele linked with a sex for a population on which some mutational factor operates influencing only on a recessive allele

……..(27)

Where D (n) there is function of time describing influence of the mutational factor on a population.

Let's search the solution of the equation (27) as:

where _{f0} there is initial (before the beginning of influence on a population) frequency of recessive allele in the woman, G(n, n^{|} ) - Green's function for a population, n^{|} - the parameter of integration – time (moment) of mutational influence the response on which is investigated in time (moment)^{/} .

Let's substitute the formula (28) in the equation (27):

Let's accept the following initial conditions reflecting that initial frequency of a recessive allele does not depend on time .

In this case the equation (29) will be transformed to a kind:

Where ẟ (n - n^{/}) there is Dirac’s function.

In (30) we used known property of Dirac’s ẟ - function

where Dirac’s ẟ - function has the following determination ẟ (n - n^{/}) = 0 at n ≠ ^{|}, and ẟ (n - ^{|}) = ∞ at n = n^{/}.

Thus from the equation (30) we find the equation for Green's function of a population:

.........(32)

The linear differential equation of the second order (32) easily is solved by introduction of a new variable.

Let's note that Green's function of a population does not depend on a kind of influence. Green's function of a population is a characteristic of a population but of not the influence kind on it. It shows how the population reacts to conditional instant mutational influence as ẟ -function in the initial moment of time n - n^{|} .

Further assuming Z=^{n }^{ln 2}. For function

......(33)

which solution is:

......(34)

In a finding (34) we are used the evenness of the ẟ-function and property (31). Hence there is. Z=^{-(}^{n-n|}^{ )ln}^{2}. Green's function of a population find by a condition:

Formula (35) can be presented as:

…….(36)

where Δn = n-n^{|} there is time interval from the moment of influence n^{|} .

The sign a minus specifies that the population always resists to the external influence.

On

Let's substitute Green's function (36) in the formula (28):

where

there is parameter of a mutation, and integration is carried out on all dimensionless time of the mutational factor influence from 0 up to n^{|}_{max}. We shall notice that as against the equation for Green's function (32) in the formula (37) time of influence not instant.

Thus the law (37) of a recessive allele frequency change at women in a population in response to mutational influence carries universal on time^{|}_{max}.

For example if mutational influence is carried out with the help of gradually radioactive element decay, D(n^{|})=D(0)2^{-}^{α}^{n|}_{1/2} , and _{1/2} there is half-life of a radioactive element. Hence the parameter of a mutation is equal:

The elementary formula for parameter of a mutation occurred if the half-life of a radioactive element is equal to average time of generation alternation, i.e. ^{|}_{max}. For example to this condition for the person satisfy half-life of such widespread radioactive elements, as caesium Cs^{137}_{55} - 30.2 years, strontium Sr^{90}_{38} - 28.7 years

If the radioactive doze is proportional to time of influence the linear dependence of parameter of a mutation on time of influence well correlates with the experimental linear curve “doze – effect” for

The populations consist of family trees set which individuals are crossed among themselves. At mathematical research of the family trees genotypes distribution it is necessary to use Hardy - Weinberg law in the discrete form, and at research of a population in the form of the differential equation of the second order.

The characteristic of a population reaction on mutational influence is Green's function. It shows how the population reacts to instant mutational influence as Dirac’s ẟ-function.

The response of the population genic system on mutational influence is levelling of this influence due to crossing individuals. The population comes back to an initial level for 4-5 generations.

The law of a recessive allele frequency change in a population in reply to mutational influence has universal character on time.