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.     

Computer Simulation 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 an 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, H ⁿ _k max D ⁿ _kof 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.     

Response of the MGE 412 pattern to truncation selection of a quantitative trait in an isogenic line of drosophila after severe heat shock (SHS)

Positive and negative selection on the total length of two fragments of an interrupted longitudinal wing vein in an isogenic line of Drosophila melanogaster was accompanied by changes in the genomic localization pattern of MGE 412. Strong truncation selection was conducted in the population of effective size Ne = 160 for 50 generations. Twenty-six out of 35 polymorphic HHS-induced segments of MGE localization behaved as independent copies and markers, whereas 9 segments proved to be selective. The second group included "hot" segments of HHS transposition induction (43B, 97E, etc.). Thus, final consensus patterns of induced MGE transpositions have a random and an adaptive component in generation 50 of positive and negative selection. Selective patterns probably include modifier MGEs, which generate induced genetic regulatory variation of polygenes controlling the selected quantitative trait in the isogenic line after HHS.     

Mobile genetic elements and quantitative characters in Drosophila: facts and hypotheses]

This review is dedicated to the comparison of the facts obtained and the proposed hypotheses, to the critical analysis of the situation arisen, and to the estimation of key propositions of the concept developed. The main point is that mobile genetic elements (MGEs) participate directly in expression, variability, selection and evolution of different quantitative characters. Genetic and selection data are considered, and hypotheses of random fixation, marker effect and direct participation of MGE patterns in expression and selection of quantitative characters are discussed. The consequences of temperature treatment are considered and hypotheses of masked selection and temperature induction of transpositions are discussed. The marker effects are shown to be non-sufficient to explain the properties of quantitative character radius incompletus system. The MGE patterns are important components of genetical system of determination of a quantitative character. MGEs modify, enhance the expression of neighbouring polygenes. Temperature effects could be explained by the influence of stress temperature treatment through the system of heat shock response on the capacity of MGEs to transcribe and transpose. The system of diversed MGE patterns in drosophila chromosomes could be believed to be universal genomic system of "soft" modification of the polygenic control of any limiting quantitative characters.     

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.     

Conflicting selection alters the trajectory of molecular evolution in a tripartite bacteria–plasmid–phage interaction

Bacteria engage in a complex network of ecological interactions, which includes mobile genetic elements (MGEs) such as phages and plasmids. These elements play a key role in microbial communities as vectors of horizontal gene transfer but can also be important sources of selection for their bacterial hosts. In natural communities, bacteria are likely to encounter multiple MGEs simultaneously and conflicting selection among MGEs could alter the bacterial evolutionary response to each MGE. Here, we test the effect of interactions with multiple MGEs on bacterial molecular evolution in the tripartite interaction between the bacterium, Pseudomonas fluorescens, the lytic bacteriophage, SBW25φ2, and conjugative plasmid, pQBR103, using genome sequencing of experimentally evolved bacteria. We show that individually, both plasmids and phages impose selection leading to bacterial evolutionary responses that are distinct from bacterial populations evolving without MGEs, but that together, plasmids and phages impose conflicting selection on bacteria, constraining the evolutionary responses observed in pairwise interactions. Our findings highlight the likely difficulties of predicting evolutionary responses to multiple selective pressures from the observed evolutionary responses to each selective pressure alone. Understanding evolution in complex microbial communities comprising many species and MGEs will require that we go beyond studies of pairwise interactions.     

Prolongation of MGE 412 transposition induction after gamma-irradiation in an isogenic line of Drosophila melanogaster]]

The dose dependence of the rate of gamma-induced transpositions and consequent dynamics of the MGE 412 pattern after gamma-irradiation were investigated in isogenic line 49 in generations F1, F12, F140, and F170. It was shown that the results on dose dependence of transpositions was very similar with the corresponding results of the classic works by Timofeeff-Ressovsky et al. (1935). It is suggested that the transcribed copies of retrotransposon 412 "cure" gamma-radiation-induced double-strand DNA breaks. The phenomenon of prolongation of MGE transposition induction during early generations after treatment was shown. In this period (F1-F12), the maximum transposition rate (lambda approximately equal to 2 x 10(-2) events per MGE copy, per haploid genome, per generation) and the maximum number of heterozygous MGE copies were achieved. In the late generations (F140 and F170), the reduced induction level (lambda approximately 10(-3) was established. In the population of effective size Ne = 2000 individuals, this corresponds to the state when lambda > 1/4Ne, i.e., when the transposition flow prevails over the MGE copy loss by genetic drift. These data together with some indirect evidence argue for the hypothesis that the spontaneous transposition rate is proportional to the average number of heterozygous MGE copies per diploid genome.     

Analysis of MGE Dm412 localization pattern in Drosophila melanogaster lines undergoing long-term selection for a quantitative character]

The subpopulation composed of the mixture of Drosophila isogenic lines with interrupted wing radial vein (mutation radius incompletus, ri) was subjected to long-term selection in different directions for increase or decrease in expression of the ri gene. As a result, the lines with contrasting different values of mean character phenotype were developed. The isogenic lines of mean character phenotype were developed. The isogenic lines and F2 from their crosses with selected lines were analysed by the pattern of copia-like MGE DM412 localization. The isogenic lines were shown to have individual pattern, the selected lines differing strongly from them. Selection led to the loss of Dm412 localization sites during negative selection, while positive selection results both in loss and acquisition of sites. Correlation between the phenotype of the quantitative character and the pattern of MGE Dm412 was revealed.     

Horizontal DNA Transfer Mechanisms of Bacteria as Weapons of Intragenomic Conflict

Horizontal DNA transfer (HDT) is a pervasive mechanism of diversification in many microbial species, but its primary evolutionary role remains controversial. Much recent research has emphasised the adaptive benefit of acquiring novel DNA, but here we argue instead that intragenomic conflict provides a coherent framework for understanding the evolutionary origins of HDT. To test this hypothesis, we developed a mathematical model of a clonally descended bacterial population undergoing HDT through transmission of mobile genetic elements (MGEs) and genetic transformation. Including the known bias of transformation toward the acquisition of shorter alleles into the model suggested it could be an effective means of counteracting the spread of MGEs. Both constitutive and transient competence for transformation were found to provide an effective defence against parasitic MGEs; transient competence could also be effective at permitting the selective spread of MGEs conferring a benefit on their host bacterium. The coordination of transient competence with cell–cell killing, observed in multiple species, was found to result in synergistic blocking of MGE transmission through releasing genomic DNA for homologous recombination while simultaneously reducing horizontal MGE spread by lowering the local cell density. To evaluate the feasibility of the functions suggested by the modelling analysis, we analysed genomic data from longitudinal sampling of individuals carrying Streptococcus pneumoniae. This revealed the frequent within-host coexistence of clonally descended cells that differed in their MGE infection status, a necessary condition for the proposed mechanism to operate. Additionally, we found multiple examples of MGEs inhibiting transformation through integrative disruption of genes encoding the competence machinery across many species, providing evidence of an ongoing “arms race.” Reduced rates of transformation have also been observed in cells infected by MGEs that reduce the concentration of extracellular DNA through secretion of DNases. Simulations predicted that either mechanism of limiting transformation would benefit individual MGEs, but also that this tactic’s effectiveness was limited by competition with other MGEs coinfecting the same cell. A further observed behaviour we hypothesised to reduce elimination by transformation was MGE activation when cells become competent. Our model predicted that this response was effective at counteracting transformation independently of competing MGEs. Therefore, this framework is able to explain both common properties of MGEs, and the seemingly paradoxical bacterial behaviours of transformation and cell–cell killing within clonally related populations, as the consequences of intragenomic conflict between self-replicating chromosomes and parasitic MGEs. The antagonistic nature of the different mechanisms of HDT over short timescales means their contribution to bacterial evolution is likely to be substantially greater than previously appreciated.     

Expression of the quantitative trait radius incompletus, temperature effects and localization of mobile genetic elements in Drosophila. II. Mobile genetic elements Dm-412

Two "selection" sub-populations (ris- and ris+), as well as two "temperature" ones (ric113 and ric149) were earlier developed from the control ric sub-population with interrupted vein of the fly wing. All five sub-populations were investigated for hybridization of MGE Dm-412 with drosophila polytene chromosomes in situ. The tree of similarity of MGE Dm-412 hybridization patterns was built by the methods of matrix clusterization. The sub-populations with the most resembling expressions of characters (ris- and ric113, ris+ and ric149) were found to be also most similar in patterns of MGE localization and their changes. Nonrandomness of these changes was shown, the similarity of patterns being demonstrated to be mainly the result of the changes. There is evidence that such effects cannot be accounted for by genetic drift and independent stochastic changes in MGE localization.     

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 Simulation of Concerted Evolution of Patterns of Polygenes, Transposable Elements, and Origin Identity Labels: Negative Selection

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.     

Sequence analysis of tyrosine recombinases allows annotation of mobile genetic elements in prokaryotic genomes

Mobile genetic elements (MGEs) sequester and mobilize antibiotic resistance genes across bacterial genomes. Efficient and reliable identification of such elements is necessary to follow resistance spreading. However, automated tools for MGE identification are missing. Tyrosine recombinase (YR) proteins drive MGE mobilization and could provide markers for MGE detection, but they constitute a diverse family also involved in housekeeping functions. Here, we conducted a comprehensive survey of YRs from bacterial, archaeal, and phage genomes and developed a sequence‐based classification system that dissects the characteristics of MGE‐borne YRs. We revealed that MGE‐related YRs evolved from non‐mobile YRs by acquisition of a regulatory arm‐binding domain that is essential for their mobility function. Based on these results, we further identified numerous unknown MGEs. This work provides a resource for comparative analysis and functional annotation of YRs and aids the development of computational tools for MGE annotation. Additionally, we reveal how YRs adapted to drive gene transfer across species and provide a tool to better characterize antibiotic resistance dissemination.     

Mobile genetic elements in plant sex evolution

The most significant theories of the appearance and maintenance of sex are presented. However, in the overwhelming majority of existing theories, the problem of sex, which is the central problem of evolutionary biology, is considered primarily through the prism of reproductive features of living organisms, whereas the issue of molecular driving forces of sexual reproduction is restricted to the possible role of mobile genetic elements (MGEs) in the appearance of sexual reproduction. The structural and functional significance of MGEs in the genomic organization of plants is illustrated. It is shown that MGEs could act as important molecular drivers of sex evolution in plants. The involvement of MGEs in the formation of sex chromosomes and possible participation in seeds-without-sex reproduction (apomixis) is discussed. Thus, the hypothesis on the active MGE participation in sex evolution is in good agreement with the modern views on pathways and directions of sex evolution in plants.     

Expression of the quantitative trait radius incompletus in Drosophila and localization of mobile elements MDG1 and copia

A series of subpopulations earlier obtained were studied for hybridization of mobile genetic elements (MGE). The subpopulations examined were two "selected" (ris- and ris+), two "temperature" (ri(c113) and ri(c149)) and the control (ric). The method of in situ hybridization with polytene chromosomes of larval salivary glands was used to determine the patterns of MGE localization for all subpopulations. The patterns obtained appeared to be quite different from that of mdg-2. The trees of similarity for subpopulations according to the patterns of every MGE localization were built by conventional clustering methods. These trees were topologically similar to each other and to mdg-2. Distinction spectra of patterns of four daughter subpopulations, in comparison with the control one, were shown for each of these MGE to be independent and individual. However, there are some common regularities among copia-like MGE-mdg-1, copia, mdg-2 and, probably, mdg-3, namely: non-random property of the majority of changes, the similarity of patterns for subpopulations with similar phenotypes etc. So, Drosophila genome can be conceived as a complex system of patterns of different MGE localization, capable of common or independent mass transpositions after external stress action.     

Computer Simulation of the Joint Evolution of the Patterns of Polygenes, Transposable Elements, and Origin Identity Labels: Stabilizing Selection

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 0.7. The inbreeding coefficient reaches plateau (at 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.     

Airborne and indigenous microbiomes co-drive the rebound of antibiotic resistome during compost storage

Composting is widely used to reduce the abundance of antibiotic resistance genes (ARGs) in solid waste. While ARG dynamics have been extensively investigated during composting, the fate and abundance of residual ARGs during the storage remain unexplored. Here, we tested experimentally how ARG and mobile genetic element (MGE) abundances change during compost storage using metagenomics, quantitative PCR and direct culturing. We found that 43.8% of ARGs and 39.9% of MGEs quickly recovered already during the first week of storage. This rebound effect was mainly driven by the regrowth of indigenous, antibiotic-resistant bacteria that survived the composting. Bacterial transmission from the surrounding air had a much smaller effect, being most evident as MGE rebound during the later stages of storage. While hyperthermophilic composting was more efficient at reducing the relative abundance of ARGs and MGEs, relatively greater ARG rebound was observed during the storage of hyperthermophilic compost, exceeding the initial levels of untreated sewage sludge. Our study reveals that residual ARGs and MGEs left in the treated compost can quickly rebound during the storage via airborne introduction and regrowth of surviving bacteria, highlighting the need to develop better storage strategies to prevent the rebound of ARGs and MGEs after composting.