Negative selection and computer models of the joint evolution of the patterns of polygenes,transposable elements,and origin identity labels

Computer simulation of the population dynamics of the genomic patterns of polygenes, transposable elements (TEs), and origin identity labels (OILs) in the course of negative selection for an additive quantitative trait has been performed. It was demonstrated that active polygene alleles disappear very rapidly, whereas the patterns of TEs and OILs continue their evolution determined by strict selective inbreeding and gene drift. Dendrograms of the patterns of polygenes, TEs, and OILs were constructed for all generations. It was demonstrated that the final consensus pattern of OILs consists of the fragments of the original patterns, which contain neither active polygene alleles nor modifier or marker TEs. Neutral TE copies were present in the final pattern, as should be expected in the case of gene drift. Inbreeding coefficient increased steadily but by generation 100 reached values higher than 0.9. All other parameters and initial conditions being the same, the responses to negative and positive selections were asymmetric.     

Computer system for simulating population dynamic patterns of polygenes and mobile genetic elements upon truncation selection for a quantitative trait

A computer system was developed for simulation of population dynamics of interacting polygene patterns and mobile genetic elements (MGEs) under selection for a quantitative trait. The system is stochastic (Monte Carlo) and takes into account the main sources of random change in the patterns (recombinations, transpositions, excisions), genetic drift, and determined trends of selection and other genetic processes in a finite population. Using this model, it is possible to analyze the dynamics of many population parameters that cannot be experimentally estimated: frequencies of polygenic alleles, proportions of adaptive and random fixations, average heterozygosities of polygenes and MGEs, coefficient of inbreeding, heritability, etc. In addition, the model can be used to test various hypotheses on polygene-MGE interaction.     

Computer modeling of joint selection response of MGE patterns and a polygenic system]

Using computer simulation, selection response of three genome patterns--polygenes, mobile genetic elements (MGEs), and labels of identity by origin (LIOs)--were studied. In each generation of selection, variability of each pattern type was described by on UPGMA tree. Stringent positive truncation selection on an additive polygenic trait and recombination between segments of the genetic map were considered. MGEs were classified into three groups: modifiers (enhancers) of the polygenic expression, markers, and independent copies. It was shown that at generations 30 to 40, 95-96% and 70-80% of respectively enforced and non-enforced active polygenic alleles were fixed (2-3% and 16-17% lost). In all generations, Hkn < or = Dkn of the length of the maximal route along the tree. At the same time, modifier MGEs were fixed for 85-88% (lost for 11-12%); marker MGEs, for 60-70 (lost for 21-25%); and independent copies, for 30-40 (lost for 50-60%). The behavior of independent MGE copies was generally consistent with the predictions of the genetic drift theory, modifier MGEs behaved similarly to the modified polygenes, and marker MGEs exhibited intermediate properties. The LIO patterns showed rapid homozygotization: their variability dropped dramatically between generations 10 and 30. In F50, the final consensus pattern of polygenes included 16 out of 18 enforced and 18 out of 21 non-enforced polygenic alleles. The fixation/loss ratios were 16:3 for modifier MGEs, 15:6 for marker MGEs, and 25:28 (with 7 polymorphic) for independent copies. The LIO consensus pattern contained 13 out of 100 original markers, which formed 26 fragments of one to ten map segments in size; 21 fragments contained active polygenic alleles, and 14 of them had also modifier MGEs. Recombinational shuffling of patterns was not completed. In the course of selection, active polygenic alleles take along adjacent segments, including those containing modifier MGEs and markers. These constitute the conservative part of all consensus patterns while the remaining segments are random.     

Population dynamics of the additive polygenic system under truncation selection]

Common features of the equations describing dynamics of the additive polygenic system under truncation selection are summarized. A combination of parameters playing the role of the effective selective pressure on the ith polygenic locus was revealed. The product of mean relative fitnesses of the individual polygenic loci, [formula: see text], was shown to play the role of relative mean fitness of the polygenic population. This value depends on the measurable parameters of the character distribution in the population: [formula: see text]. It was shown that under the constant population number during truncation selection, the characteristic of the best genotype increases, [formula: see text]; which is also a product of the frequencies of preferable genotypes at individual polygenic loci. This value plays the role of the proportion of the number of the best ("champion") genotype in the population. In fact, this is the champion genotype polygene consensus pattern frequency, which a priori indicates the possibility of the champion pattern fixation. The analogue of Haldane's dilemma for the polygenic system which restrict the number of polygenes simultaneously subjected to adaptive evolution [formula: see text] was obtained for the case of constant effective population number (Ne = const).     

Stabilizing selection and computer models of the joint evolution of patterns of polygenes,transposable elements,and origin identity labels

A computer model of the populations dynamics of the patterns of polygenes, transposable elements (TEs), and origin identity labels (OILs) in the course of stabilizing selection for an additive quantitative trait (with the target value being 0.4 of the maximum) was analyzed. It was demonstrated that the final plateaus of the trait value and the frequencies of the active values of polygenes are reached rapidly, namely, within five to seven generations (the effective selection period). The inbreeding coefficient during this period also grows rapidly and then gradually increases eventually reaching approximately 0.7. The inbreeding coefficient reaches plateau (at approximately 1.0) only in generations 300-350, which suggests the effect of gene drift. Dendrograms of the patterns of polygenes, TEs, and OILs were constructed for all generations. By generation 100 of selection, the final patterns of TEs and OILs were not formed completely. Fixations and losses, especially those of the OIL pattern, were delayed. In general, however, the population heterogeneity with respect to the patterns studied does not stabilize. This heterogeneity decreases the case of stabilizing selection, although more slowly than in the cases of positive and negative selections.     

Truncation family selection and nonsystematic inbreeding leads to a rapid fixation of a pattern of mobile genetic elements in a computer model

A computer simulation model of the population dynamics of a polygenic system and a pattern of mobile genetic elements (MGEs) under directional truncation selection for a quantitative trait was developed. Modifier MGEs were shown to be rapidly and adaptively fixed (or lost) together with the modified polygenes. Marker MGEs and independent MGE copies were fixed and lost just as rapidly but in a random manner. Using specific marking of initial haploid genomes and direct computing of the mean proportion of identical encounters at each locus in each generation, it was shown that the mean nonselective inbreeding coefficient F(n) dramatically increases in the course of selection, reaching values 0.7-0.9 in 15-20 generations. As a result, adaptive homozygotization of polygenes and modifier MGEs and random homozygotization of marker MGEs, independent MGE copies, and all other genes of the genome occurs. These results confirm the hypothesis on the "champion" polygene pattern advanced earlier to explain the data of selection experiments.