- Research Article
- Open Access

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

- Yue Wu
^{1, 2, 3}, - Richard Z. Aandahl
^{1, 2}and - Mark M. Tanaka
^{1, 2}Email authorView ORCID ID profile

**Received:**18 September 2015**Accepted:**6 December 2015**Published:**21 December 2015

## Abstract

### 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.

## Keywords

- Transposable elements
- Mobile DNA
- Bacterial evolution
- Stress response
- Transposition burst
- Regulatory mechanism
- Simulation model
- Adaptation

## Background

Insertion sequences (ISs) are simple, widely observed mobile genetic elements that only contain genes related to transposition and the regulation of transposition [1–5]. Transposition events can shift, replicate, or delete copies of ISs within a genome, which may induce mutations that change fitness of the host cell [3, 6, 7]. IS-induced mutations are often deleterious [5, 7], which raises the question of how those elements became abundant in a bacterial population [8–10]. To answer this question we need to understand how IS dynamics are governed by the rates of IS movement and resulting mutational effects.

Several hypotheses have been proposed to explain the persistence of mobile DNA in a prokaryote population. The selfish DNA hypothesis asserts that ISs are able to persist through their ability to self-replicate while making no fitness contribution to the genome [11, 12]. Analogous to sex in diploid populations, horizontal gene transfer (HGT) has been considered a major determinant of IS spread among prokaryotic species [13]. However, the role of HGT in maintaining ISs within a population is still under some debate [4, 9, 14, 15]. Condit et al. [16] showed that observed rates of HGT are too low to maintain ISs as parasites. Bichsel et al. [17] used a model to point out that although the rate of HGT is generally small, if the fitness cost of an IS is even smaller, HGT allows ISs to invade and persist for long periods in an asexual population with low probability. Bichsel et al. [18] proposed that occasional beneficial IS-induced mutations may be important in order to reach observed IS distributions in a realistic period of time.

An alternative to the selfish DNA hypothesis is that adaptive mutations play a vital role in the persistence of ISs [4, 18, 19]. ISs have been considered as a source of genetic diversity [15, 20], because IS movements mediate changes that are sometimes beneficial to their host genomes [21–23]. If these beneficial mutations prevail in a population through natural selection, ISs can hitchhike to fixation alongside them [24, 25]. Experimental evidence has shown that ISs can increase organismal fitness and thus promote adaptive evolution [6, 26–28]. One important example of this phenomenon is the spread of ISs in bacterial species under antibiotic exposure as the elements are involved in the expression and mobilisation of antibiotic resistance genes, which has been a focus of research in recent decades [29–34]. Using a simulation model Edwards and Brookfield [35] showed that mobile DNA sequences can be maintained in a clonal species if transient beneficial insertions appear for sufficiently long periods of time in one of two alternating environments. Following this work, McGraw and Brookfield [36] derived an optimal transposition rate for element maintenance if reversible advantageous mutations exist in a fluctuating environment.

Mobile elements may also randomly drift to fixation in bacterial genomes [37, 38]. The effect of drift becomes weaker as the size of the population increases [37, 39], and bacterial populations may be large enough to efficiently eliminate mobile elements that cause detrimental effects [40]. However, if the fitness costs induced by IS movements are small or close to neutral, ISs would be hard to eliminate [41]. A study of large-scale genomic data showed that empirical distributions of IS elements are compatible with selective neutrality of elements and a high deletion rate [38]. The three major hypotheses, namely the selfish DNA, the adaptive and the neutral hypotheses, should not be regarded as mutually exclusive for explaining the abundance of ISs, given the complicated interactions among horizontal gene transfer, the distribution of fitness effects, drift, and the rate of transposition [41–43].

Both IS elements and their hosts can encode mechanisms that suppress transposition activities [4, 28, 44–46]. The rate of IS transposition is often suppressed through such negative regulation [3]. The release from regulation of IS transpositions in response to environmental challenges has been observed [33, 45, 47–52]. We call this elevation of the transposition rate a *transposition burst*. An IS transposition burst can promote the adaptation of *Escherichia coli* to a high-osmolarity environment by increasing the rate of beneficial mutations [53]. These changes in transposition activities are sometimes understood as stress responses to changing environments [45, 48, 54] and suggest another hypothesis: periodic transposition bursts increase the IS copy number and thereby promote the persistence of ISs [51]. However, increased transposition rates can also lead to an accumulation of deleterious mutations in the genome. The amplification of fitness costs induced by erratic IS movements may even drive the host population to extinction via the Muller’s ratchet effect [40, 41, 55, 56]. Following the occurrence of transposition bursts, regulation of bursts may evolve and limit transposition rates during an IS invasion [57–61]. Using a simulation model of transposable elements in sexual diploid populations, le Rouzic and Capy [51] suggested that an initial burst followed by strong regulation of transposition can lead to successful invasion of transposable elements.

In this study, we examine whether transposition bursts and the regulation of bursts help to maintain ISs in an asexual population. We investigate how various mutational effects mediated by IS movements interact with transposition rates to influence the dynamics of ISs. Transposition bursts may create more IS copies along with genomic diversity which is essential for bacterial evolution. On the other hand, the extra deleterious IS-induced changes is a net burden for host cells. We show how cells carrying ISs are often eliminated by the population even with the possibility of generating advantageous mutations.

## Methods

We introduce a simulation model of the movement of ISs within genomes and consider the impact of movement on cell fitness and population dynamics. Let the bacterial population be of constant size *N*. Let *x*
_{
i
} be the frequency of cells with arrangement *i* where *i*=1,…,*n* and where *n* is the number of different arrangements of IS elements in the population (*n* is a dynamic variable). Each arrangement *i* is associated with a fitness *w*
_{
i
} and a copy number *l*
_{
i
} of IS elements. The total number of IS copies in a population at time *t* is \(\sum _{i=1}^{n} x_{i}(t) l_{i}\).

### Reproduction

*t*. This assumes that one generation is the average time for cellular reproduction. The vector of IS arrangement frequencies in generation

*t*+1 is thus

where \(\bar {w}(t) = \sum _{i} w_{i}(t)x_{i} (t)\) is proportion to the mean fitness in the population.

### Transposition events

Within a cell, we define three transposition events that affect ISs. The underlying transposition parameters (*θ*, *μ* and *ν* defined below) are constant per element, per cell, per generation for all cells. But we introduce a factor *λ*(*t*) which scales the underlying transposition probabilities and allows the transposition rates to change over time. With probability *θ*
*λ*(*t*) an IS may *shift* to a new location in the cell. With probability *μ*
*λ*(*t*) an IS *insertion* event duplicates the IS by adding a new copy at another location in the genome. With probability *ν*
*λ*(*t*) an IS *excision* event deletes an IS element from a cell.

The numbers of shift, insertion and excision events in a cell of arrangement *i* are distributed as Poisson with parameters *λ*(*t*)*θ*
*l*
_{
i
}, *λ*(*t*)*μ*
*l*
_{
i
} and *λ*(*t*)*ν*
*l*
_{
i
} respectively. For each event a new genome arrangement is created (*n* is set to *n*+1) and the count of the arrangement that experienced the event is decremented by 1. Under a shift event the copy number *l*
_{
i
} is preserved in the new *n*
^{
t
h
} arrangement, but the fitness of the new arrangement *w*
_{
n
} may be altered as described below. Under an insertion event the copy number increases by one and under a deletion event the copy number decreases by one. Again the fitness of the new arrangement may be altered relative to the original arrangement. We cap the number of ISs per cell by setting a maximum of 100 insertion sites.

### Regulation of bursts

*burst*, and model the strength of an initial burst using parameter

*δ*≥0. We allow

*λ*(

*t*) to decrease with time after the initial burst at rate

*γ*. This decrease models the action of natural selection lowering the deleterious effects of transposition. Thus

where *T* is the time when changes in transposition activities occur.

*no burst*), there is no burst of transposition (

*δ*=0) so that

*λ*=1 for all

*t*and all cells in the population always have the same underlying transposition probabilities (dotted). Under the second scenario (

*burst, no regulation*), a transposition burst event occurs (

*δ*>0) at time

*t*=

*T*and transposition rates are permanently elevated from that time (dashed). In the third scenario (

*burst, evolving regulation*) a transposition burst occurs at time

*t*=

*T*which is followed by the gradual evolution of down-regulation of transposition rates so that the original underlying transposition probability is approached asymptotically (solid). The speed of this process of evolved intracellular down-regulation is controlled by parameter

*γ*which is large for fast evolution (grey) and small for slow evolution (black solid).

Parameters with default values

Symbol | Description | Default value |
---|---|---|

| Population size | 10 |

| Shift probability per IS per cell per generation | 10 |

| Insertion probability per IS per cell per generation | 10 |

| Excision probability per IS per cell per generation | 10 |

| Initial burst strength | 99 |

| Evolutionary rate of regulator | 10 |

| Magnitude of fitness benefit per transposition event per cell generation | 2×10 |

| Magnitude of fitness cost per transposition event per cell generation | 2×10 |

| Probability of IS fitness benefit | 0.05 |

| Probability of IS fitness cost | 0.8 |

### Fitness effects

*p*

_{ a }and a deleterious change with probability

*p*

_{ d }. The magnitudes of the selective effects are

*a*for advantageous changes and −

*d*for deleterious changes. The IS copy number of a genome of type

*i*can be partitioned into numbers of elements that produced advantageous (

*l*

_{ i,p }), neutral (

*l*

_{ i,n }) and deleterious (

*l*

_{ i,d }) effects, so that

*l*

_{ i }=

*l*

_{ i,a }+

*l*

_{ i,n }+

*l*

_{ i,d }. In our model for IS-induced fitness changes, we assume that excision or shift events reverse the fitness effect that the IS originally induced. If a transposition event occurs in cell type

*i*the new arrangement is given a new index

*n*

^{′}and has fitness given by:

where the prime (^{′}) indicates the new value of the variable *n* after the new arrangement is generated (after this operation *n* is set to *n*
^{′}).

### Initialisation

Several studies have analysed the distribution of IS copy number in bacterial genomes and observed similar count distributions [9, 14, 18]; the published distributions of IS copy number are right skewed with some distributions highly right skewed. Sawyer et al. [62] reported distributions of six unrelated ISs in a collection of 71 *Escherichia coli* strains from various natural isolates. For our numerical work we use the distribution of *I*
*S*5 from Sawyer’s study, namely, α=(0,1,2,3,4,5,6,9,21) IS copies in B =(46,12,3,2,2,2,2,1,1) genomes. This reflects a typical IS copy number distribution in bacterial populations, where a majority of the genomes contain no ISs and a moderate number of cells contain higher numbers of ISs. We define β to be the distribution of IS copy number in a population and initialise the distribution of ISs in the population with α ISs in β cells, where β∼Multinomial(*N*,**B**/71). We also consider a second scenario where the population is initialised with a single IS element in a single cell, so that we can investigate the impact of transposition bursts on the initial invasion of IS elements.

### Simulation and parameters

At each generation, cell reproduction is followed by potential IS transposition events. For any cell type *i* where *i*=1,…,*n*, these processes may lead to new cell types with new fitness based on current fitness *w*
_{
i
}, and with a new IS copy number based on *l*
_{
i
}. Unless otherwise specified, we use the parameter values shown in Table 1. Observed rates of IS transposition are approximately *θ*, *μ*, *ν*=10^{−8} per cell generation [3, 63]. To expedite computation while preserving the same rate of supply of transpositions [39], we set a default population of size 10^{5} with transposition rates *θ*, *μ*, *ν*=10^{−6}.

We set the magnitude of positive and negative selection coefficients per IS per cell per generation, *a* and *d*, to range broadly from 10^{−5}– 10^{−3} because empirical estimates of these parameters are scarce and cover a wide range [7, 17, 18]. We assume that an IS-induced mutation is beneficial to its host cell with probability *p*
_{
a
}=0.05, deleterious with probability *p*
_{
d
}=0.8, and neutral otherwise; this is based on estimates of effects of insertional mutations in Elena et al. [7]. By varying the value of *p*
_{
a
} and adjusting the other probabilities correspondingly at the same time that the transposition burst occurs, we also simulate a scenario where a new environmental condition changes the distribution of IS-mediated mutational effects.

### Output variables

In each simulation run, we record the following output variables in the population over time: the proportion of cells carrying ISs with the proportion of IS copies in advantageous, deleterious and neutral sites; the mean IS copy number per cell; and the mean fitness of cell population. We track the 5 *%* and 95 *%* quantiles of IS copy numbers carried by a cell in the population, which indicate the range of IS counts in genomes. We also track the 5 *%* and 95 *%* quantiles of cell fitness in the population to examine the impact of IS transpositions on the organismal fitness. From 1000 simulation runs, we calculate the probability of IS extinction by computing the proportion of runs in which all ISs are eliminated by the population. At the end of 30,000 simulated generations, we describe the distribution of ISs in the population by calculating the mean proportion of cells carrying ISs, mean IS copy number per cell and the mean fitness of cell population across multiple runs. We increase the number of simulation runs to 10,000 when investigating the invasion of a single IS element in a population with no ISs.

## Results

We study IS extinction and persistence over 30,000 generations under three settings of transposition rates (shown in Fig. 1). In the first scenario (*no burst*), transposition occurs at a relatively low and constant rate. In the second scenario (*burst, no regulation*), we introduce a 100-fold step change in transposition probability at the 5000^{
t
h
} generation. We use this scenario to model a burst in transposition activities which may be induced by a sudden change in the environment. In the third scenario (*burst, evolving regulation*), we introduce a transposition burst at the 5000^{
t
h
} generation but also allow regulation of the transposition burst to evolve over time.

*%*range of IS copies and an increase in variation can be seen after the burst in transposition rate at the 5000

^{ t h }generation. The mean fitness of the population gradually increases over time in both cases of persistence and extinction (black curve, bottom panels).

^{ t h }generation. The left panel shows that constant and relatively low transposition rates (

*no burst*) lead to roughly 26

*%*IS survival by the end of the simulation under the default settings (blue). If bursts are unregulated the probability of survival of ISs is reduced to around 16

*%*(red). When transposition bursts are gradually down-regulated (black), it is possible for the survival probability of ISs to reach a similar level as in the

*no burst*scenario. Although only reversible IS-induced effects have been considered in this model, it is also possible for IS movements to “leave behind” the original fitness effects. Therefore, in the Additional file 1 (Section 1), we provide comparison between reversible and irreversible changes for their impact on the survival of ISs. Although the reversibility of mutational effects improves the survival of ISs, the effect appears to be slight. In contrast, the size of the population has a strong influence on the survival curve as shown in the right panel of Fig. 3.

*%*of the simulation runs the proportion of cells carrying ISs was greater than 90

*%*while the distribution formed from the remaining simulations was right skewed (left). The distribution of the mean IS copy number is right skewed (middle). The cell fitness distribution is also right skewed and mostly above 1 (right) due to the selection of beneficial mutations. All corresponding results under other transposition scenarios and reversible mutational effects are provided in the Additional file 1 (Section 2). From this point onwards, we focus on the model of transposition bursts with evolving down-regulation unless otherwise specified.

### Role of transposition bursts and their regulation

*δ*and the rate of evolution of burst regulation

*γ*on the survival of ISs (left), the mean cell fitness (middle) and the mean IS copy number per cell (right) after 30,000 generations. As the rate of evolution of regulation increases, ISs are more likely to persist but with lower mean cell fitness and mean IS copy per cell. This pattern is most evident in the case of strong bursts (

*δ*=999), where down-regulation of transposition is expected to have the greatest effect. We highlight the observation that strong bursts with weak regulation (

*δ*=999,

*γ*=10

^{−3.5}) produce not only low survival of ISs but also high mean fitness and high IS copy numbers. This phenomenon is due to a selection bias whereby cells that survive the burst tend to be the ones with high fitness. The burst itself elevates copy number.

### Role of fitness effect distribution

*a*and negative

*d*selection (top), and the effects of the probabilities of advantageous

*p*

_{ a }and deleterious

*p*

_{ d }mutations (bottom) on the distribution of ISs and the mean cell fitness at the end of simulations where the ISs did not go extinct. As expected, low deleterious effects

*d*and high advantageous effects

*a*generally lead to greater survival of ISs and greater mean fitness. The threshold value of

*a*below which IS are expected to be lost is presented as a vertical dashed line in panel A, and the derivation of this threshold is provided in the end of ‘Results’ section (Eq. 3). While

*a*also influences the mean proportion (panel B) and the mean fitness (panel C) of cells carrying ISs, the deleterious effect size

*d*does not.

We varied *p*
_{
a
} and *p*
_{
d
} to examine the case where a change in the environment makes available new beneficial mutations to be generated by IS transpositions (Fig. 6, bottom panels, black). In terms of the model, the value of *p*
_{
a
} changes from 0 to positive values after 5000 generations. We find that the mean cell proportion carrying ISs is gradually improved as *p*
_{
a
} increases (panel F, black), as are the mean population fitness (panel G, black) and the mean IS copy per cell (panel H, black). For comparison we also examined whether IS persistence is still favoured by more IS-induced adaptive mutations in the absence of bursts in transposition rates (Fig. 6, bottom panels, thick grey). We find that increase in the probability of fitness benefit alone does not help ISs to persist in a population.

### Invasion by a single IS element

### Heuristic model of persistence of IS elements

To explain the broad patterns of IS persistence, we consider here a simplified version of the simulation model to find conditions under which IS elements are expected to persist evolutionarily. Let us ignore polymorphism and assume deleterious mutations do not reach fixation. Let *p*
_{
n
} be the proportion of IS-induced changes that are selectively neutral.

*λ*(

*t*)

*μ*

*p*

_{ a }

*N*. The supply rate of neutral changes leading to an increase in copy number is

*λ*(

*t*)

*μ*

*p*

_{ n }

*N*. The supply rate of changes leading to a decrease in copy number is

*λ*(

*t*)

*ν*

*p*

_{ n }assuming that the loss of elements does not change fitness. The probability of fixation is approximately 2

*a*for advantageous changes [64] and 1/

*N*for neutral changes [65]. (Here consider only one change at a time and ignore deleterious insertions in the background of advantageous insertions). Therefore the rate of substitution of copy-increasing events is

*λ*(

*t*) speeds up the extinction of ISs through increased

*K*

_{(−)}. However, the rush to extinction due to a burst is mitigated by the evolution of down-regulation. Assuming ISs do not go extinct near this zero boundary a rough criterion for the long term success of ISs is given by the ratio

*K*

_{(+)}/

*K*

_{(−)}being greater than unity. That is, IS persistence requires

Because the burst function *λ*(*t*) is cancelled out in the above threshold it ultimately does not influence the long term success of ISs except in hastening their initial extinction when near the zero boundary.

The condition (3) is similar in form to the persistence condition for mobile elements given by Lynch [66]. The differences are that we disregard deleterious changes (which are also eventually neglected in Lynch’s treatment), and we include a term for adaptive changes. Adaptive changes play a role when the product *a*
*p*
_{
a
} is large enough when compared to the reciprocal of the population size, 1/*N*.

In summary, the success of insertion sequences depends not only on the balance between new insertions and deletions but also on the strength of positive selection of IS-induced changes which in turn depends on both population size and the distribution of fitness effects, as revealed in expression (3).

## Discussion

In this study we have investigated the role of transposition bursts on the survival of ISs in a bacterial population by modelling the effects of IS movements on the fitness of host genomes. Bursts in transposition activity hasten the extinction of mobile elements because they increase both excision events and IS-induced fitness costs. In other words, elevated transposition serves to amplify both deleterious and advantageous fitness changes but this accelerates the extinction of elements, particularly if the mean fitness effect of IS-induced mutations is negative. However, transposition bursts create genetic diversity which occasionally generates advantageous mutations and thus help organisms adapt to new environments. If transposition activity is eventually down-regulated following a transposition burst, the elimination of ISs from a population is slowed down.

Because transposition bursts can rapidly increase the number of IS copies, it has been proposed that bursts favour the invasion of ISs [51, 67]. We find, however, that since invading ISs are initially rare and close to the extinction boundary, they are already vulnerable to quick extinction [17]. The invasion of ISs in an asexual population is not helped by bursts.

Transposition bursts have been postulated as a stress response for bacterial populations facing environmental challenges [23, 49, 68]. Bursts in IS transposition have been reported to be associated with harsh environmental conditions [45, 48, 50]. As an example, antibiotic treatment results in an increased frequency of IS transposition which may accelerate the development of resistance [33, 52]. Under our model, strong transposition bursts increase the mean fitness of organisms through selection of beneficial mutations. Bacteria can occasionally benefit from such increased genetic variability which allows for adaptation to new environments [49, 69]. However, the amplification of deleterious mutations increases the risk of host lineage extinction [40, 55]. It also favours genomes lacking IS elements. If a population eliminates ISs through natural selection, it loses the ability to use IS elements for adaptation under environmental stress in the future. Hence, transposition bursts can accidentally promote bacterial evolution but they are unlikely to have evolved as an adaptive strategy for either ISs or genomes. Transposition bursts alone cannot explain the abundance of ISs over long-term periods of time.

Regulation of transposition bursts decelerates the burst-induced extinction of ISs while creating some genomic variability. Both IS elements and hosts appear to have evolved regulatory mechanisms that limit transposition activities [44, 49, 59]. Natural selection favours bacterial hosts that evolve the ability to suppress transposition bursts in response to the amplified effects of deleterious mutations [57, 68]. Transposition bursts with evolving regulation represent a compromise between providing a host population with a means of adapting to new environments on one hand and promoting the persistence of ISs on the other [51]. However, if the mean fitness effect on their hosts is negative, a given family of IS elements may still go extinct eventually.

Advantageous effects mediated by IS transposition may serve as a necessary condition for the spread of ISs in prokaryotes [4, 18]. ISs can be viewed as mutator genes whose prevalence increases in bacterial populations over successive selective sweeps as they hitchhike with beneficial mutations that they produce [35, 70], as confirmed in experimental studies [24, 26]. Our model shows how benefits of occasional adaptive mutations can be amplified by transposition bursts and overcome intermittent drops in the mean population fitness. With an initial IS distribution with a low mean copy number [62], we find that ISs gradually go extinct in the population, and the process of IS elimination is slower if positive selection on beneficial mutations is strong or if the transposition rate is low. A high mean fitness effect of IS movements can be achieved by either IS-induced mutations carrying higher fitness benefits or higher rates of beneficial mutations. How those mutational effects vary and change the distribution of ISs in genomes under environmental stress is a topic worth investigating experimentally. Previous studies have shown that reversible IS-induced mutations with precise excisions may benefit cells in fluctuating environments [35, 36]. We have compared non-reversible versus reversible mutation effects for their influence on the persistence of ISs; although reversibility may be beneficial the difference between the two models is slight.

Our model has been kept simple to focus on the effect of transposition bursts on the dynamics of ISs. However, it could be made more biologically realistic by including recombination. Given that ISs can promote recombination in bacterial strains [71], one can include recombination in future models to study how it influences the IS dynamics by creating more genomic variation. Recombination may also cause deletions of ISs [71]. Genomic evidence has shown that transposable elements may go extinct periodically in bacterial lineages [14]. A single IS element that is newly introduced in a population with no other ISs can rapidly go extinct due to drift and purifying selection [17]. Therefore horizontal transfer of ISs from other species is essential for the initial introduction and re-introduction of IS elements [10, 13, 17, 40, 66]. Alternatively, ISs can be maintained in a population if the intraspecific horizontal transfer rate is high. Over evolutionary time scales, ISs can only be maintained in a bacterial population if the acquisition and establishment of mobile elements can outrun their elimination. The rate of IS elimination in turn can be decreased through advantageous mutations and the regulation of transposition bursts. It would be interesting to further study the balance of these processes.

## Conclusions

We find that transposition bursts do not help insertion sequences persist; rather, they accelerate the elimination of ISs from populations because IS-induced changes are mostly harmful to host genomes. Since ISs are often lost during transposition bursts, we do not consider ISs or transpotion bursts to have evolved as an adaptive strategy to deal with environmental stress. However, bacteria in new environments can occasionally benefit from the genetic variation generated by bursts in IS transposition activity, and the burst-induced extinction of ISs can be decelerated by evolved down-regulation of transposition.

## Availability of supporting data

We provide the computer simulation code used in our study in the Additional file 1, Section 3.

## Declarations

### Acknowledgements

We appreciate valuable comments from the editors and two anonymous reviewers which have led to substantial revision of this work. We thank Andrew Francis and Lloyd Sanders for their constructive comments and Rich Edwards for useful discussions. This work was supported by the Australian Research Council through grant FT140100398 awarded to MMT and an Australian Postgraduate Award to YW.

**Open Access** This article is distributed under the terms of the Creative Commons Attribution 4.0 International License (http://creativecommons.org/licenses/by/4.0/), which permits unrestricted use, distribution, and reproduction in any medium, provided you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated.

## Authors’ Affiliations

## References

- Chandler M, Mahillon J. Mobile DNA In: Craig NL, Craigie R, Gellert M, Lambowitz AM, editors. Washington: ASM press: 2002. p. 305–66.Google Scholar
- Siguier P, Filée J, Chandler M. Insertion sequences in prokaryotic genomes. Curr Opin Microbiol. 2006; 9(5):526–31. doi:http://0-dx.doi.org.brum.beds.ac.uk/10.1016/j.mib.2006.08.005. Antimicrobials/Genomics.PubMedView ArticleGoogle Scholar
- Kleckner N. Transposable elements in prokaryotes. Annu Rev Genet. 1981; 15:341–404.PubMedView ArticleGoogle Scholar
- Blot M. Transposable elements and adaptation of host bacteria. Genetica. 1994; 93:5–12. 10.1007/BF01435235.PubMedView ArticleGoogle Scholar
- Mahillon J, Chandler M. Insertion sequences. Microbiol Mol Biol Rev. 1998; 62(3):725–74.PubMedPubMed CentralGoogle Scholar
- Chao L, McBroom SM. Evolution of transposable elements: an IS10 insertion increases fitness in
*Escherichia coli*. Mol Biol Evol. 1985; 2(5):359–69.PubMedGoogle Scholar - Elena SF, Ekunwe L, Hajela N, Oden SA, Lenski RE. Distribution of fitness effects caused by random insertion mutations in
*Escherichia coli*. Genetica. 1998; 102–103(0):349–58. doi:http://0-dx.doi.org.brum.beds.ac.uk/10.1023/A:1017031008316.PubMedView ArticleGoogle Scholar - Hickey DA. Evolutionary dynamics of transposable elements in prokaryotes and eukaryotes. Genetica. 1992; 86:269–74.PubMedView ArticleGoogle Scholar
- Touchon M, Rocha EPC. Causes of insertion sequences abundance in prokaryotic genomes. Mol Biol Evol. 2007; 24(4):969–81. doi:http://0-dx.doi.org.brum.beds.ac.uk/10.1093/molbev/msm014.PubMedView ArticleGoogle Scholar
- Wagner A. Transposable elements as genomic diseases. Mol BioSyst. 2009; 5:32–5.PubMedView ArticleGoogle Scholar
- Orgel LE, Crick FHC. Selfish DNA: the ultimate parasite. Nature. 1980; 284:604–7.PubMedView ArticleGoogle Scholar
- Doolittle WF, Sapienza C. Selfish genes, the phenotype paradigm and genome evolution. Nature. 1980; 284(5757):601–3.PubMedView ArticleGoogle Scholar
- Kelly BG, Vespermann A, Bolton DJ. The role of horizontal gene transfer in the evolution of selected foodborne bacterial pathogens. Food Chem Toxicol. 2009; 47(5):951–68. doi:http://0-dx.doi.org.brum.beds.ac.uk/10.1016/j.fct.2008.02.006.PubMedView ArticleGoogle Scholar
- Wagner A. Periodic extinctions of transposable elements in bacterial lineages: evidence from intragenomic variation in multiple genomes. Mol Biol Evol. 2006; 23(4):723–33.PubMedView ArticleGoogle Scholar
- Bickhart DM, Gogarten JP, Lapierre P, Tisa LS, Normand P, Benson DR. Insertion sequence content reflects genome plasticity in strains of the root nodule actinobacterium
*Frankia*. BMC Genomics. 2009; 10(1):468. doi:http://0-dx.doi.org.brum.beds.ac.uk/10.1186/1471-2164-10-468.PubMedPubMed CentralView ArticleGoogle Scholar - Condit R, Stewart FM, Levin BR. The population biology of bacterial transposons: a priori conditions for maintenance as parasitic DNA. Amer Nat. 1988; 132(1):129–47.View ArticleGoogle Scholar
- Bichsel M, Barbour AD, Wagner A. The early phase of a bacterial insertion sequence infection. Theor Popul Biol. 2010; 78(4):278–88. doi:http://0-dx.doi.org.brum.beds.ac.uk/10.1016/j.tpb.2010.08.003.PubMedView ArticleGoogle Scholar
- Bichsel M, Barbour AD, Wagner A. Estimating the fitness effect of an insertion sequence. J Math Biol. 2013; 66(1–2):95–114. doi:http://0-dx.doi.org.brum.beds.ac.uk/10.1007/s00285-012-0504-2.PubMedView ArticleGoogle Scholar
- Le Rouzic A, Deceliere G. Models of the population genetics of transposable elements. Genet Res. 2005; 85:171–81. doi:http://0-dx.doi.org.brum.beds.ac.uk/10.1017/S0016672305007585.PubMedView ArticleGoogle Scholar
- McEvoy CRE, Falmer AA, Gey van Pittius NC, Victor TC, van Helden PD, Warren RM. The role of IS6110 in the evolution of
*Mycobacterium tuberculosis*. Tuberculosis. 2007; 87(5):393–404. doi:http://0-dx.doi.org.brum.beds.ac.uk/10.1016/j.tube.2007.05.010.PubMedView ArticleGoogle Scholar - Mijnendonckx K, Provoost A, Monsieurs P, Leys N, Mergeay M, Mahillon J, et al. Insertion sequence elements in
*Cupriavidus metallidurans*CH34: distribution and role in adaptation. Plasmid. 2011; 65:193–203.PubMedView ArticleGoogle Scholar - El Gharniti F, Dols-Lafargue M, Bon E, Claisse O, Miot-Sertier C, Lonvaud A, et al. IS30 elements are mediators of genetic diversity in
*Oenococcus oeni*. Int J Food Microbiol. 2012; 158(1):14–22. doi:http://0-dx.doi.org.brum.beds.ac.uk/10.1016/j.ijfoodmicro.2012.06.009.PubMedView ArticleGoogle Scholar - Casacuberta E, Gonzalez J. The impact of transposable elements in environmental adaptation. Mol Ecol. 2013; 22(6):1503–1517. doi:http://0-dx.doi.org.brum.beds.ac.uk/10.1111/mec.12170.PubMedView ArticleGoogle Scholar
- Chao L, Vargas C, Spear BB, Cox EC. Transposable elements as mutator genes in evolution. Nature. 1983; 303:633–5. doi:http://0-dx.doi.org.brum.beds.ac.uk/10.1038/303633a0.PubMedView ArticleGoogle Scholar
- Martiel J, Blot M. Transposable elements and fitness of bacteria. Theor Popul Biol. 2002; 61(4):509–18. doi:http://0-dx.doi.org.brum.beds.ac.uk/10.1006/tpbi.2002.1603.PubMedView ArticleGoogle Scholar
- Schneider D, Lenski RE. Dynamics of insertion sequence elements during experimental evolution of bacteria. Res. Microbiol. 2004; 155(5):319–27. doi:http://0-dx.doi.org.brum.beds.ac.uk/10.1016/j.resmic.2003.12.008. Genome plasticity and the evolution of microbial genomes.PubMedView ArticleGoogle Scholar
- Gaffé J, McKenzie C, Maharjan RP, Coursange E, Ferenci T, Schneider D. Insertion sequence-driven evolution of
*Escherichia coli*in chemostats. J Mol Evol. 2011; 72(4):398–412. doi:http://0-dx.doi.org.brum.beds.ac.uk/10.1007/s00239-011-9439-2.PubMedView ArticleGoogle Scholar - Hall BG. Transposable elements as activators of cryptic genes in
*E. coli*. Genetica. 1999; 107(1–3):181–7. doi:http://0-dx.doi.org.brum.beds.ac.uk/10.1023/A:1003936706129.PubMedView ArticleGoogle Scholar - Nakatsu C, Ng J, Singh R, Straus N, Wyndham C. Chlorobenzoate catabolic transposon Tn5271 is a composite class I element with flanking class II insertion sequences. Proc Natl Acad Sci USA. 1991; 88(19):8312–316.PubMedPubMed CentralView ArticleGoogle Scholar
- Poirel L, Decousser JW, Nordmann P. Insertion sequence ISEcp1B is involved in expression and mobilization of a
*b**l**a*_{ C T X−M }*β*-lactamase gene. Antimicrob Agents Chemother. 2003; 47(9):2938–945. doi:http://0-dx.doi.org.brum.beds.ac.uk/10.1128/AAC.47.9.2938-2945.2003.PubMedPubMed CentralView ArticleGoogle Scholar - Boutoille D, Corvec S, Caroff N, Giraudeau C, Espaze E, Caillon J, et al. Detection of an IS21 insertion sequence in the
*mexR*gene of*Pseudomonas aeruginosa*increasing*β*-lactam resistance. FEMS Microbiol Lett. 2004; 230(1):143–6. doi:http://0-dx.doi.org.brum.beds.ac.uk/10.1016/S0378-1097(03)00882-6.PubMedView ArticleGoogle Scholar - Tetu SG, Holmes AJ. A family of insertion sequences that impacts integrons by specific targeting of gene cassette recombination sites, the IS1111-attC group. J Bacteriol. 2008; 190(14):4959–970. doi:http://0-dx.doi.org.brum.beds.ac.uk/10.1128/JB.00229-08.PubMedPubMed CentralView ArticleGoogle Scholar
- Nagel M, Reuter T, Jansen A, Szekat C, Bierbaum G. Influence of ciprofloxacin and vancomycin on mutation rate and transposition of IS256 in
*Staphylococcus aureus*. Internatl J Medical Microbiol. 2011; 301(3):229–36. doi:http://0-dx.doi.org.brum.beds.ac.uk/10.1016/j.ijmm.2010.08.021.View ArticleGoogle Scholar - Casagrande Proietti P, Bietta A, Coletti M, Marenzoni ML, Scorza AV, Passamonti F. Insertion sequence IS256 in canine pyoderma isolates of
*Staphylococcus pseudintermedius*associated with antibiotic resistance. Vet Microbiol. 2012; 157:376–82.PubMedView ArticleGoogle Scholar - Edwards RJ, Brookfield JFY. Transiently beneficial insertions could maintain mobile DNA sequences in variable environments. Mol Biol Evol. 2003; 20(1):30–7. doi:http://0-dx.doi.org.brum.beds.ac.uk/10.1093/molbev/msg001.PubMedView ArticleGoogle Scholar
- McGraw JE, Brookfield JFY. The interaction between mobile DNAs and their hosts in a fluctuating environment. J Theor Biol. 2006; 243:13–23. doi:http://0-dx.doi.org.brum.beds.ac.uk/10.1016/j.jtbi.2006.05.029.PubMedView ArticleGoogle Scholar
- Lynch M, Conery JS. The origins of genome complexity. Science. 2003; 302(5649):1401–1404. doi:http://0-dx.doi.org.brum.beds.ac.uk/10.1126/science.1089370.PubMedView ArticleGoogle Scholar
- Iranzo J, Gómez MJ, López de Saro FJ, Manrubia S. Large-scale genomic analysis suggests a neutral punctuated dynamics of transposable elements in bacterial genomes. PLoS Comput Biol. 2014; 10(6):1003680. doi:http://0-dx.doi.org.brum.beds.ac.uk/10.1371/journal.pcbi.1003680.View ArticleGoogle Scholar
- Le Rouzic A, Boutin TS, Capy P. Long-term evolution of transposable elements. Proc Natl Acad Sci USA. 2007; 104(49):19375–19380. doi:http://0-dx.doi.org.brum.beds.ac.uk/10.1073/pnas.0705238104.PubMedPubMed CentralView ArticleGoogle Scholar
- Rankin DJ, Bichsel M, Wagner A. Mobile DNA can drive lineage extinction in prokaryotic populations. J Evol Biol. 2010; 23(11):2422–431. doi:http://0-dx.doi.org.brum.beds.ac.uk/10.1111/j.1420-9101.2010.02106.x.PubMedView ArticleGoogle Scholar
- Dolgin ES, Charlesworth B. The fate of transposable elements in asexual populations. Genetics. 2006; 174(2):817–27. doi:http://0-dx.doi.org.brum.beds.ac.uk/10.1534/genetics.106.060434.PubMedPubMed CentralView ArticleGoogle Scholar
- Brookfield JFY. The ecology of the genome — mobile DNA elements and their hosts. Nat Rev Genet. 2005; 6(2):128–36.PubMedView ArticleGoogle Scholar
- Deceliere G, Letrillard Y, Charles S, Biémont C. TESD: a transposable element dynamics simulation environment. Bioinformatics. 2006; 22(21):2702–703. doi:http://0-dx.doi.org.brum.beds.ac.uk/10.1093/bioinformatics/btl454.PubMedView ArticleGoogle Scholar
- Nagy Z, Chandler M. Regulation of transposition in bacteria. Res Microbiol. 2004; 155:387–98. doi:http://0-dx.doi.org.brum.beds.ac.uk/10.1016/j.resmic.2004.01.008.PubMedView ArticleGoogle Scholar
- Foster PL. Stress-induced mutagenesis in bacteria. Crit Rev Biochem Mol Biol. 2007; 42(5):373–97. doi:http://0-dx.doi.org.brum.beds.ac.uk/10.1080/10409230701648494.PubMedPubMed CentralView ArticleGoogle Scholar
- Hua-Van A, Le Rouzic A, Boutin TS, Filée J, Capy P. The struggle for life of the genome’s selfish architects. Biol Direct. 2011; 6(1):19. doi:http://0-dx.doi.org.brum.beds.ac.uk/10.1186/1745-6150-6-19.PubMedPubMed CentralView ArticleGoogle Scholar
- Pfeifer F, Blaseio U. Transposition burst of the ISH27 insertion element family in
*Halobacterium halobium*. Nucleic Acids Res. 1990; 18(23):6921–925. doi:http://0-dx.doi.org.brum.beds.ac.uk/10.1093/nar/18.23.6921.PubMedPubMed CentralView ArticleGoogle Scholar - Eichenbaum Z, Livneh Z. UV light induces IS10 transposition in
*Escherichia coli*. Genetics. 1998; 149(3):1173–1181.PubMedPubMed CentralGoogle Scholar - Capy P, Gasperi G, Biémont C, Bazin C. Stress and transposable elements: co-evolution or useful parasites?. Heredity. 2000; 85(2):101–6.PubMedView ArticleGoogle Scholar
- de Visser JAGM, Akkermans ADL, Hoekstra RF, de Vos WM. Insertion-sequence-mediated mutations isolated during adaptation to growth and starvation in
*Lactococcus lactis*. Genetics. 2004; 168:1145–1157. doi:http://0-dx.doi.org.brum.beds.ac.uk/10.1534/genetics.104.032136.PubMedPubMed CentralView ArticleGoogle Scholar - Le Rouzic A, Capy P. The first steps of transposable elements invasion: parasitic strategy
*vs*, genetic drift. Genetics. 2005; 169(2):1033–43. doi:http://0-dx.doi.org.brum.beds.ac.uk/10.1534/genetics.104.031211.PubMedPubMed CentralView ArticleGoogle Scholar - Schreiber F, Szekat C, Josten M, Sahl HG, Bierbaum G. Antibiotic-induced autoactivation of IS256 in
*Staphylococcus aureus*. Antimicrob Agents Chemother. 2013; 57(12):6381–384. doi:http://0-dx.doi.org.brum.beds.ac.uk/10.1128/AAC.01585-13.PubMedPubMed CentralView ArticleGoogle Scholar - Stoebel DM, Dorman CJ. The effect of mobile element IS10 on experimental regulatory evolution in
*Escherichia coli*. Mol Biol Evol. 2010; 27(9):2105–112. doi:http://0-dx.doi.org.brum.beds.ac.uk/10.1093/molbev/msq101.PubMedPubMed CentralView ArticleGoogle Scholar - Mennecier S, Servant P, Coste G, Bailone A, Sommer S. Mutagenesis via IS transposition in
*Deinococcus radiodurans*. Mol Microb. 2006; 59(1):317–25. doi:http://0-dx.doi.org.brum.beds.ac.uk/10.1111/j.1365-2958.2005.04936.x.View ArticleGoogle Scholar - Tanaka MM, Small PM, Salamon H, Feldman MW. The dynamics of repeated elements: applications to the epidemiology of
*tuberculosis*. Proc Natl Acad Sci USA. 2000; 97(7):3532–537. doi:http://0-dx.doi.org.brum.beds.ac.uk/10.1073/pnas.97.7.3532.PubMedPubMed CentralView ArticleGoogle Scholar - Arkhipova I, Meselson M. Deleterious transposable elements and the extinction of asexuals. BioEssays. 2005; 27(1):76–85. doi:http://0-dx.doi.org.brum.beds.ac.uk/10.1002/bies.20159.PubMedView ArticleGoogle Scholar
- Charlesworth B, Langley CH. The evolution of self-regulated transposition of transposable elements. Genetics. 1986; 112(2):359–83.PubMedPubMed CentralGoogle Scholar
- Coen D, Lemaitre B, Delattre M, Quesneville H, Ronsseray S, Simonelig M, et al. Drosophila
*P*element: transposition, regulation and evolution. Genetica. 1994; 93:61–78. doi:http://0-dx.doi.org.brum.beds.ac.uk/10.1007/BF01435240.PubMedView ArticleGoogle Scholar - Labrador M, Corces VG. Transposable element-host interactions: regulation of insertion and excision. Annu Rev Genet. 1997; 37:381–401.View ArticleGoogle Scholar
- Kidwell MG, Lisch DR. Perspective: transposable elements, parasitic DNA, and genome evolution. Evolution. 2001; 55(1):1–24.PubMedView ArticleGoogle Scholar
- Le Rouzic A, Capy P. Theoretical approaches to the dynamics of transposable elements in genomes, populations, and species In: Lankenau D-H, Volff J-N, editors. Transposons and the Dynamic Genome. Genome Dynamics and Stability. Berlin: Springer: 2009. p. 1–19.Google Scholar
- Sawyer SA, Dykhuizen DE, DuBose RF, Green L, Mutangadura-Mhlanga T, Wolczyk DF, et al. Distribution and abundance of insertion sequences among naturalss isolates of
*Escherichia coli*. Genetics. 1987; 115(1):51–63.PubMedPubMed CentralGoogle Scholar - Sousa A, Bourgard C, Wahl LM, Gordo I. Rates of transposition in
*Escherichia coli*. Biol Lett. 2013; 9(6). doi:http://0-dx.doi.org.brum.beds.ac.uk/10.1098/rsbl.2013.0838. - Haldane JBS. A mathematical theory of natural and artificial selection, part v: Selection and mutation. Math Proc Camb Philos Soc. 1927; 23:838–44. doi:http://0-dx.doi.org.brum.beds.ac.uk/10.1017/S0305004100015644.View ArticleGoogle Scholar
- Kimura M. On the probability of fixation of mutant genes in a population. Genetics. 1962; 47(6):713–9.PubMedPubMed CentralGoogle Scholar
- Lynch M. The frailty of adaptive hypotheses for the origins of organismal complexity. Proc Natl Acad Sci. 2007; 104(Suppl 1):8597–604. doi:http://0-dx.doi.org.brum.beds.ac.uk/10.1073/pnas.0702207104.PubMedPubMed CentralView ArticleGoogle Scholar
- Nuzhdin S. Sure facts, speculations, and open questions about the evolution of transposable element copy number. Genetica. 1999; 107(1):129–37. doi:http://0-dx.doi.org.brum.beds.ac.uk/10.1023/A:1003957323876.PubMedView ArticleGoogle Scholar
- Twiss E, Coros AM, Tavakoli NP, Derbyshire KM. Transposition is modulated by a diverse set of host factors in
*Escherichia coli*and is stimulated by nutritional stress. Mol Microbiol. 2005; 57(6):1593–607. doi:http://0-dx.doi.org.brum.beds.ac.uk/10.1111/j.1365-2958.2005.04794.x.PubMedView ArticleGoogle Scholar - Tenaillon O, Denamur E, Matic I. Evolutionary significance of stress-induced mutagenesis in bacteria. Trends Microbiol. 2004; 12(6):264–70. doi:http://0-dx.doi.org.brum.beds.ac.uk/10.1016/j.tim.2004.04.002.PubMedView ArticleGoogle Scholar
- Chao L, Cox EC. Competition between high and low mutating strains of
*Escherichia coli*. Evolution. 1983; 37(1):125–34.View ArticleGoogle Scholar - Beare PA, Unsworth N, Andoh M, Voth DE, Omsland A, Gilk SD, et al. Comparative genomics reveal extensive transposon-mediated genomic plasticity and diversity among potential effector proteins within the genus
*Coxiella*. Infect Immun. 2009; 77(2):642–56. doi:http://0-dx.doi.org.brum.beds.ac.uk/10.1128/IAI.01141-08.PubMedPubMed CentralView ArticleGoogle Scholar