Transiently beneficial insertions could maintain mobile DNA sequences in variable environments

The dynamic network of IS30 transposition pathways

The E. coli element IS30 has adopted the copy-out-paste-in transposition mechanism that is prevalent in a number of IS-families. As an initial step, IS30 forms free circular transposition intermediates like IS minicircles or tandem IS-dimers by joining the inverted repeats of a single element or two, sometimes distantly positioned IS copies, respectively. Then, the active IR-IR junction of these intermediates reacts with the target DNA, which generates insertions, deletions, inversions or cointegrates. The element shows dual target specificity as it can insert into hot spot sequences or next to its inverted repeats. In this study the pathways of rearrangements of transposition-derived cointegrate-like structures were examined. The results showed that the probability of further rearrangements in these structures depends on whether the IS elements are flanked by hot spot sequences or take part in an IR-IR junction. The variability of the deriving products increases with the number of simultaneously available IRs and IR-IR joints in the cointegrates or the chromosome. Under certain conditions, the parental structures whose transposition formed the cointegrates are restored and persist among the rearranged products. Based on these findings, a novel dynamic model has been proposed for IS30, which possibly fits to other elements that have adopted the same transposition mechanism. The model integrates the known transposition pathways and the downstream rearrangements occurring after the formation of different cointegrate-like structures into a complex network. Important feature of this network is the presence of “feedback loops” and reversible transposition rearrangements that can explain how IS30 generates variability and preserves the original genetic constitution in the bacterial population, which contributes to the adaptability and evolution of host bacteria.

How sequence populations persist inside bacterial genomes

Compared to their eukaryotic counterparts, bacterial genomes are small and contain extremely tightly packed genes. Repetitive sequences are rare but not completely absent. One of the most common repeat families is REPINs. REPINs can replicate in the host genome and form populations that persist for millions of years. Here, we model the interactions of these intragenomic sequence populations with the bacterial host. We first confirm well-established results, in the presence and absence of horizontal gene transfer (hgt) sequence populations either expand until they drive the host to extinction or the sequence population gets purged from the genome. We then show that a sequence population can be stably maintained, when each individual sequence provides a benefit that decreases with increasing sequence population size. Maintaining a sequence population of stable size also requires the replication of the sequence population to be costly to the host, otherwise the sequence population size will increase indefinitely. Surprisingly, in regimes with high hgt rates, the benefit conferred by the sequence population does not have to exceed the damage it causes to its host. Our analyses provide a plausible scenario for the persistence of sequence populations in bacterial genomes. We also hypothesize a limited biologically relevant parameter range for the provided benefit, which can be tested in future experiments.

Evolutionary dynamics of transposable elements in bdelloid rotifers

Transposable elements (TEs) are selfish genomic parasites whose ability to spread autonomously is facilitated by sexual reproduction in their hosts. If hosts become obligately asexual, TE frequencies and dynamics are predicted to change dramatically, but the long-term outcome is unclear. Here, we test current theory using whole-genome sequence data from eight species of bdelloid rotifers, a class of invertebrates in which males are thus far unknown. Contrary to expectations, we find a variety of active TEs in bdelloid genomes, at an overall frequency within the range seen in sexual species. We find no evidence that TEs are spread by cryptic recombination or restrained by unusual DNA repair mechanisms. Instead, we find that that TE content evolves relatively slowly in bdelloids and that gene families involved in RNAi-mediated TE suppression have undergone significant expansion, which might mitigate the deleterious effects of active TEs and compensate for the consequences of long-term asexuality.

Dynamics of bacterial insertion sequences: can transposition bursts help the elements persist?

Background

Currently there is no satisfactory explanation for why bacterial insertion sequences (ISs) widely occur across prokaryotes despite being mostly harmful to their host genomes. Rates of horizontal gene transfer are likely to be too low to maintain ISs within a population. IS-induced beneficial mutations may be important for both prevalence of ISs and microbial adaptation to changing environments but may be too rare to sustain IS elements in the long run. Environmental stress can induce elevated rates of IS transposition activities; such episodes are known as ‘transposition bursts’. By examining how selective forces and transposition events interact to influence IS dynamics, this study asks whether transposition bursts can lead to IS persistence.

Results

We show through a simulation model that ISs are gradually eliminated from a population even if IS transpositions occasionally cause advantageous mutations. With beneficial mutations, transposition bursts create variation in IS copy numbers and improve cell fitness on average. However, these benefits are not usually sufficient to overcome the negative selection against the elements, and transposition bursts amplify the mean fitness effect which, if negative, simply accelerates the extinction of ISs. If down regulation of transposition occurs, IS extinctions are reduced while ISs still generate variation amongst bacterial genomes.

Conclusions

Transposition bursts do not help ISs persist in a bacterial population in the long run because most burst-induced mutations are deleterious and therefore not favoured by natural selection. However, bursts do create more genetic variation through which occasional advantageous mutations can help organisms adapt. Regulation of IS transposition bursts and stronger positive selection of the elements interact to slow down the burst-induced extinction of ISs.

Electronic supplementary material

The online version of this article (doi:10.1186/s12862-015-0560-5) contains supplementary material, which is available to authorized users.

Extinction probabilities and stationary distributions of mobile genetic elements in prokaryotes: The birth-death-diversification model

Theoretical approaches are essential to our understanding of the complex dynamics of mobile genetic elements (MGEs) within genomes. Recently, the birth-death-diversification model was developed to describe the dynamics of mobile promoters (MPs), a particular class of MGEs in prokaryotes. A unique feature of this model is that genetic diversification of elements was included. To explore the implications of diversification on the longterm fate of MGE lineages, in this contribution we analyze the extinction probabilities, extinction times and equilibrium solutions of the birth-death-diversification model. We find that diversification increases both the survival and growth rate of MGE families, but the strength of this effect depends on the rate of horizontal gene transfer (HGT). We also find that the distribution of MGE families per genome is not necessarily monotonically decreasing, as observed for MPs, but may have a peak in the distribution that is related to the HGT rate. For MPs specifically, we find that new families have a high extinction probability, and predict that the number of MPs is increasing, albeit at a very slow rate. Additionally, we develop an extension of the birth-death-diversification model which allows MGEs in different regions of the genome, for example coding and non-coding, to be described by different rates. This extension may offer a potential explanation as to why the majority of MPs are located in non-promoter regions of the genome.

Birth, death, and diversification of mobile promoters in prokaryotes

A previous study of prokaryotic genomes identified large reservoirs of putative mobile promoters (PMPs), that is, homologous promoter sequences associated with nonhomologous coding sequences. Here we extend this data set to identify the full complement of mobile promoters in sequenced prokaryotic genomes. The expanded search identifies nearly 40,000 PMP sequences, 90% of which occur in noncoding regions of the genome. To gain further insight from this data set, we develop a birth-death-diversification model for mobile genetic elements subject to sequence diversification; applying the model to PMPs we are able to quantify the relative importance of duplication, loss, horizontal gene transfer (HGT), and diversification to the maintenance of the PMP reservoir. The model predicts low rates of HGT relative to the duplication and loss of PMP copies, rapid dynamics of PMP families, and a pool of PMPs that exist as a single copy in a genome at any given time, despite their mobility. We report evidence of these "singletons" at high frequencies in prokaryotic genomes. We also demonstrate that including selection, either for or against PMPs, was not necessary to describe the observed data.

Genomic parasites or symbionts? Modeling the effects of environmental pressure on transposition activity in asexual populations

Transposable elements are DNA segments capable of persisting in host genomes by self-replication in spite of deleterious mutagenic effects. The theoretical dynamics of these elements within genomes has been studied extensively, and population genetic models predict that they can invade and maintain as a result of both intra-genomic and inter-individual selection in sexual species. In asexuals, the success of selfish DNA is more difficult to explain. However, most theoretical work assumes constant environment. Here, we analyze the impact of environmental change on the dynamics of transposition activity when horizontal DNA exchange is absent, based on a stochastic computational model of transposable element proliferation. We argue that repeated changes in the phenotypic optimum in a multidimensional fitness landscape may induce explosive bursts of transposition activity associated with faster adaptation. However, long-term maintenance of transposition activity is unlikely. This could contribute to the significant variation in the transposable element copy number among closely related species.

Metaproteomics reveals abundant transposase expression in mutualistic endosymbionts

Transposases, enzymes that catalyze the movement of mobile genetic elements, are the most abundant genes in nature. While many bacteria encode an abundance of transposases in their genomes, the current paradigm is that the expression of transposase genes is tightly regulated and generally low due to its severe mutagenic effects. In the current study, we detected the highest number of transposase proteins ever reported in bacteria, in symbionts of the gutless marine worm Olavius algarvensis with metaproteomics. At least 26 different transposases from 12 different families were detected, and genomic and proteomic analyses suggest that many of these are active. This high expression of transposases indicates that the mechanisms for their tight regulation have been disabled or no longer exist.

Relaxed natural selection alone does not permit transposable element expansion within 4,000 generations in Escherichia coli

Insertion sequences (ISs) are transposable genetic elements in bacterial genomes. IS elements are common among bacteria but are generally rare within free-living species, probably because of the negative fitness effects they have on their hosts. Conversely, ISs frequently proliferate in intracellular symbionts and pathogens that recently transitioned from a free-living lifestyle. IS elements can profoundly influence the genomic evolution of their bacterial hosts, although it is unknown why they often expand in intracellular bacteria. We designed a laboratory evolution experiment with Escherichia coli K-12 to test the hypotheses that IS elements often expand in intracellular bacteria because of relaxed natural selection due to (1) their generally small effective population sizes (N (e)) and thus enhanced genetic drift, and (2) their nutrient rich environment, which makes many biosynthetic genes unnecessary and thus selectively neutral territory for IS insertion. We propagated 12 populations under four experimental conditions: large N (e) versus small N (e), and nutrient rich medium versus minimal medium. We found that relaxed selection over 4,000 generations was not sufficient to permit IS element expansion in any experimental population, thus leading us to hypothesize that IS expansion in intracellular symbionts may often be spurred by enhanced transposition rates, possibly due to environmental stress, coupled with relaxed natural selection.

Mobile DNA can drive lineage extinction in prokaryotic populations

Natural selection ultimately acts on genes and other DNA sequences. Adaptations that are good for the gene can have adverse effects at higher levels of organization, including the individual or the population. Mobile genetic elements illustrate this principle well, because they can self-replicate within a genome at a cost to their host. As they are costly and can be transmitted horizontally, mobile elements can be seen as genomic parasites. It has been suggested that mobile elements may cause the extinction of their host populations. In organisms with very large populations, such as most bacteria, individual selection is highly effective in purging genomes of deleterious elements, suggesting that extinction is unlikely. Here we investigate the conditions under which mobile DNA can drive bacterial lineages to extinction. We use a range of epidemiological and ecological models to show that harmful mobile DNA can invade, and drive populations to extinction, provided their transmission rate is high and that mobile element-induced mortality is not too high. Population extinction becomes more likely when there are more elements in the population. Even if elements are costly, extinction can still occur because of the combined effect of horizontal gene transfer, a mortality induced by mobile elements. Our study highlights the potential of mobile DNA to be selected at the population level, as well as at the individual level.

Correlations between bacterial ecology and mobile DNA

Several factors can affect the density of mobile DNA in bacterial genomes including rates of exposure to novel gene pools, recombination, and reductive evolution. These traits are difficult to measure across a broad range of bacterial species, but the ecological niches occupied by an organism provide some indication of the relative magnitude of these forces. Here, by analyzing 384 bacterial genomes assigned to three ecological categories (obligate intracellular, facultative intracellular, and extracellular), we address two, related questions: How does the density of mobile DNA vary across the Bacteria? And is there a statistically supported relationship between ecological niche and mobile element gene density? We report three findings. First, the fraction of mobile element genes in bacterial genomes ranges from 0 to 21% and decreases significantly: facultative intracellular > extracellular > obligate intracellular bacteria. Results further show that the obligate intracellular bacteria that host switch have a higher mobile DNA gene density than the obligate intracellular bacteria that are vertically transmitted. Second, while bacteria from the three ecological niches differ in their average mobile DNA contents, the ranges of mobile DNA found in each category overlap a surprising extent, suggesting bacteria with different lifestyles can tolerate similar amounts of mobile DNA. Third, mobile DNA gene densities increase with genome size across the entire dataset, and the significance of this correlation is dependent on the obligate intracellular bacteria. Further, mobile DNA gene densities do not correlate with evolutionary relationships in a 16S rDNA phylogeny. These findings statistically support a compelling link between mobile element evolution and bacterial ecology.

Distant horizontal gene transfer is rare for multiple families of prokaryotic insertion sequences

Horizontal gene transfer in prokaryotes is rampant on short and intermediate evolutionary time scales. It poses a fundamental problem to our ability to reconstruct the evolutionary tree of life. Is it also frequent over long evolutionary distances? To address this question, we analyzed the evolution of 2,091 insertion sequences from all 20 major families in 438 completely sequenced prokaryotic genomes. Specifically, we mapped insertion sequence occurrence on a 16S rDNA tree of the genomes we analyzed, and we also constructed phylogenetic trees of the insertion sequence transposase coding sequences. We found only 30 cases of likely horizontal transfer among distantly related prokaryotic clades. Most of these horizontal transfer events are ancient. Only seven events are recent. Almost all of these transfer events occur between pairs of human pathogens or commensals. If true also for other, non-mobile DNA, the rarity of distant horizontal transfer increases the odds of reliable phylogenetic inference from sequence data.

Frequency-Dependent Selection in a Periodic Environment

We examine the action of natural selection in a periodically changing environment where two competing strains are specialists respectively for each environmental state. When the relative fitness of the strains is subject to a very general class of frequency-dependent selection, we show that coexistence rather than extinction is the likely outcome. This coexistence may be a stable periodic equilibrium, stable limit cycles of varying lengths, or be deterministically chaotic. Our model is applicable to the population dynamics commonly found in many types of viruses.

The fate of transposable elements in asexual populations

Sexual reproduction and recombination are important for maintaining a stable copy number of transposable elements (TEs). In sexual populations, elements can be contained by purifying selection against host carriers with higher element copy numbers; however, in the absence of sex and recombination, asexual populations could be driven to extinction by an unchecked proliferation of TEs. Here we provide a theoretical framework for analyzing TE dynamics under asexual reproduction. Analytic results show that, in an infinite asexual population, an equilibrium in copy number is achieved if no element excision is possible, but that all TEs are eliminated if there is some excision. In a finite population, computer simulations demonstrate that small populations are driven to extinction by a Muller's ratchet-like process of element accumulation, but that large populations can be cured of vertically transmitted TEs, even with excision rates well below transposition rates. These results may have important consequences for newly arisen asexual lineages and may account for the lack of deleterious retrotransposons in the putatively ancient asexual bdelloid rotifers.

The interaction between mobile DNAs and their hosts in a fluctuating environment

The interaction between mobile DNA sequences and their hosts raises important questions in the context of hosts which reproduce clonally with only rare horizontal transmission between clones. The activity of some mobile DNAs as reversible mutators of genes raises the possibility that, in a fluctuating environment, cells may gain an advantage if they have mobile DNAs which mutate genes whose inactivation is favoured in one of the environments that the population encounters. Here we analyse a model of this process and ask what would be the optimal rate of transposition in a population whose elements are maintained by this mechanism. We also examine the impact of horizontal transfer on such a population. With movement of elements between cells, we can imagine elements with differing rates of transposition and host cells with differing rates of transposition. We find that evolution in the population of elements favours a rapid rate of transposition, and evolution of the host cells favours cells in which this rapid rate of element-dependent transposition results in an optimal rate of transposition per cell. However, when horizontal transfer rates are high, some unexpected features of the model are observed. In particular, a polymorphism between cell types (some with an optimal rate of transposition and some with no transposition at all from endogenous elements) can be stably maintained. We consider the possible biological predictions of this analysis.

Retrotransposon Sequence Variation in Four Asexual Plant Species

Transposable elements (TEs) can be viewed as genetic parasites that persist in populations due to their capacity for increase in copy number and the inefficacy of selection against them. A corollary of this hypothesis is that TEs are more likely to spread within sexual populations and be eliminated or inactivated within asexual populations. While previous work with animals has shown that asexual taxa may contain less TE diversity than sexual taxa, comparable work with plants has been lacking. Here we report the results of a study of Ty1/copia, Ty3/gypsy, and LINE-like retroelement diversity in four asexual plant species. Retroelement-like sequences, with a high degree of conservation both within and between species, were isolated from all four species. The sequences correspond to several previously annotated retroelement subfamilies. They also exhibit a pattern of nucleotide substitution characterized by an excess of synonymous substitutions, suggestive of a history of purifying selection. These findings were compared with retroelement sequence evolution in sexual plant taxa. One likely explanation for the discovery of conserved TE sequences in the genomes of these asexual taxa is simply that asexuality within these taxa evolved relatively recently, such that the loss and breakdown of TEs is not yet detectable through analysis of sequence diversity. This explanation is examined by conducting stochastic simulation of TE evolution and by using published information to infer rough estimates of the ages of asexual taxa.

Periodic Extinctions of Transposable Elements in Bacterial Lineages: Evidence from Intragenomic Variation in Multiple Genomes

Most previous work on the evolution of mobile DNA was limited by incomplete sequence information. Whole genome sequences allow us to overcome this limitation. I study the nucleotide diversity of prominent members of five insertion sequence families whose transposition activity is encoded by a single transposase gene. Eighteen among 376 completely sequenced bacterial genomes and plasmids carry between 3 and 20 copies of a given insertion sequence. I show that these copies generally show very low DNA divergence. Specifically, more than 68% of the transposase genes are identical within a genome. The average number of amino acid replacement substitutions at amino acid replacement sites is Ka = 0.013, that at silent sites is Ks = 0.1. This low intragenomic diversity stands in stark contrast to a much higher divergence of the same insertion sequences among distantly related genomes. Gene conversion among protein-coding genes is unlikely to account for this lack of diversity. The relation between transposition frequencies and silent substitution rates suggests that most insertion sequences in a typical genome are evolutionarily young and have been recently acquired. They may undergo periodic extinction in bacterial lineages. By implication, they are detrimental to their host in the long run. This is also suggested by the highly skewed and patchy distribution of insertion sequences among genomes. In sum, one can think of insertion sequences as slow-acting infectious diseases of cell lineages.

Mobile genetic elements and sexual reproduction

Transposable elements (TE) are prominent components of most eukaryotic genomes. In addition to their possible participation in the origin of sexual reproduction in eukaryotes, they may be also involved in its maintenance as important contributors to the deleterious mutation load. Comparative analyses of transposon content in the genomes of sexually reproducing and anciently asexual species may help to understand the contribution of different TE classes to the deleterious load. The apparent absence of deleterious retrotransposons from the genomes of ancient asexuals is in agreement with the hypothesis that they may play a special role in the maintenance of sexual reproduction and in early extinction for which most species are destined upon the abandonment of sex.

Evolutionary forces generating sequence homogeneity and heterogeneity within retrotransposon families

Genome projects allow us to sample copies of a retrotransposon sequence family residing in a host genome. The variation in DNA sequence between these individual copies will reflect the evolutionary process that has spread the sequences through the genome. Here I review quantitatively the expected diversity of elements belonging to a transposable genetic element family. I use a simple neutral model for replicative mobile DNAs such as retrotransposons to predict the extent of sequence variability between members of a single family of transposable elements, both within and between species. The effects of horizontal transfer are also explored. I also consider the impact on these distributions of an increase in transposition rate arising from a mutational change in copy of the sequence. In addition, I consider the question of the interaction between retrotransposons and their hosts, and the causes of the abundance of transposable elements in the genomes that they occupy.

Modeling the amplification dynamics of human Alu retrotransposons

Retrotransposons have had a considerable impact on the overall architecture of the human genome. Currently, there are three lineages of retrotransposons (Alu, L1, and SVA) that are believed to be actively replicating in humans. While estimates of their copy number, sequence diversity, and levels of insertion polymorphism can readily be obtained from existing genomic sequence data and population sampling, a detailed understanding of the temporal pattern of retrotransposon amplification remains elusive. Here we pose the question of whether, using genomic sequence and population frequency data from extant taxa, one can adequately reconstruct historical amplification patterns. To this end, we developed a computer simulation that incorporates several known aspects of primate Alu retrotransposon biology and accommodates sampling effects resulting from the methods by which mobile elements are typically discovered and characterized. By modeling a number of amplification scenarios and comparing simulation-generated expectations to empirical data gathered from existing Alu subfamilies, we were able to statistically reject a number of amplification scenarios for individual subfamilies, including that of a rapid expansion or explosion of Alu amplification at the time of human–chimpanzee divergence.

Nearly 50% of the human genome is composed of mobile elements. While much of this sequence consists of inactive “fossil” elements that are no longer actively moving or generating new copies, three families are currently proliferating in human genomes. Among these, the Alu lineage has reached a copy number of over 1 million and alone accounts for approximately 10% of the genome. While considerable evidence has been gathered concerning the underlying biological mechanisms of Alu mobilization and proliferation, a detailed understanding of Alu amplification history is currently lacking. Researchers are aware, for example, that several thousand Alu elements have inserted within the human genome since the divergence of humans and chimpanzees, but how those insertions were distributed over this ~6-million-year time period is currently unknown. In this work, the authors introduce a simulation framework that seeks to incorporate both sequence diversity and empirically gathered population data from human Alu elements, in order to provide a better understanding of the last several million years of human Alu evolution. The results suggest that a rapid explosion of Alu amplification at the time of the human–chimpanzee divergence is unlikely. Therefore, it is improbable that an increase in Alu retrotransposition activity was involved in the speciation of humans and chimpanzees.

The first steps of transposable elements invasion: parasitic strategy vs. genetic drift

Transposable elements are often considered as selfish DNA sequences able to invade the genome of their host species. Their evolutive dynamics are complex, due to the interaction between their intrinsic amplification capacity, selection at the host level, transposition regulation, and genetic drift. Here, we propose modeling the first steps of TE invasion, i.e., just after a horizontal transfer, when a single copy is present in the genome of one individual. If the element has a constant transposition rate, it will disappear in most cases: the elements with low-transposition rate are frequently lost through genetic drift, while those with high-transposition rate may amplify, leading to the sterility of their host. Elements whose transposition rate is regulated are able to successfully invade the populations, thanks to an initial transposition burst followed by a strong limitation of their activity. Self-regulation or hybrid dysgenesis may thus represent some genome-invasion parasitic strategies.

The ecology of the genome - mobile DNA elements and their hosts

The genomes of multicellular eukaryotes provide information that determines the phenotype. However, not all sequences in the genome are required for this purpose. Other sequences are often selfish in their actions and interact in complex ways. Here, an analogy is developed between the components of the genome, including mobile DNA elements, and an ecological community. Unlike ecological communities, however, the slow rates at which genomes change allow us to reconstruct patterns of interaction that stretch back tens or hundreds of millions of years.

Evolutionary capacitance as a general feature of complex gene networks

An evolutionary capacitor buffers genotypic variation under normal conditions, thereby promoting the accumulation of hidden polymorphism. But it occasionally fails, thereby revealing this variation phenotypically. The principal example of an evolutionary capacitor is Hsp90, a molecular chaperone that targets an important set of signal transduction proteins. Experiments in Drosophila and Arabidopsis have demonstrated three key properties of Hsp90: (1) it suppresses phenotypic variation under normal conditions and releases this variation when functionally compromised; (2) its function is overwhelmed by environmental stress; and (3) it exerts pleiotropic effects on key developmental processes. But whether these properties necessarily make Hsp90 a significant and unique facilitator of adaptation is unclear. Here we use numerical simulations of complex gene networks, as well as genome-scale expression data from yeast single-gene deletion strains, to present a mechanism that extends the scope of evolutionary capacitance beyond the action of Hsp90 alone. We illustrate that most, and perhaps all, genes reveal phenotypic variation when functionally compromised, and that the availability of loss-of-function mutations accelerates adaptation to a new optimum phenotype. However, this effect does not require the mutations to be conditional on the environment. Thus, there might exist a large class of evolutionary capacitors whose effects on phenotypic variation complement the systemic, environment-induced effects of Hsp90.