# The Diversity of Prokaryotic DDE Transposases of the Mutator Superfamily, Insertion Specificity, and Association with Conjugation Machineries

Romain Guérillot, Patricia Siguier, Edith Gourbeyre, Michael Chandler, Philippe Glaser

<figure><img src="image_1.png" /></figure>

## ► To cite this version:

Romain Guérillot, Patricia Siguier, Edith Gourbeyre, Michael Chandler, Philippe Glaser. The Diversity of Prokaryotic DDE Transposases of the Mutator Superfamily, Insertion Specificity, and Association with Conjugation Machineries. Genome Biology and Evolution, 2014, 6 (2), pp.260-272. 10.1093/gbe/evu010. hal-01324840

HAL Id: hal-01324840

https://hal.sorbonne-universite.fr/hal-01324840

Submitted on 1 Jun 2016

HAL is a multi-disciplinary open access archive for the deposit and dissemination of scientific research documents, whether they are published or not. The documents may come from teaching and research institutions in France or abroad, or from public or private research centers.

L'archive ouverte pluridisciplinaire HAL, est destinée au dépôt et à la diffusion de documents scientifiques de niveau recherche, publiés ou non, émanant des établissements d'enseignement et de recherche français ou étrangers, des laboratoires publics ou privés.

<figure><img src="image_2.png" /></figure>

Distributed under a Creative Commons Attribution 4.0 International License
# The Diversity of Prokaryotic DDE Transposases of the Mutator Superfamily, Insertion Specificity, and Association with Conjugation Machineries

Romain Guérillot<fnref n="1" />,<fnref n="2" />,<fnref n="3" />, Patricia Siguier<fnref n="4" />, Edith Gourbeyre<fnref n="4" />, Michael Chandler<fnref n="4" />, and Philippe Glaser<fnref n="1" />,<fnref n="2" />,*

Accepted: January 6, 2014

## Abstract

Transposable elements (TEs) are major components of both prokaryotic and eukaryotic genomes and play a significant role in their evolution. In this study, we have identified new prokaryotic DDE transposase families related to the eukaryotic Mutator-like transposases. These genes were retrieved by cascade PSI-Blast using as initial query the transposase of the streptococcal integrative and conjugative element (ICE) TnGBS2. By combining secondary structure predictions and protein sequence alignments, we predicted the DDE catalytic triad and the DNA-binding domain recognizing the terminal inverted repeats. Furthermore, we systematically characterized the organization and the insertion specificity of the TEs relying on these prokaryotic Mutator-like transposases (p-MULT) for their mobility. Strikingly, two distant TE families target their integration upstream σA dependent promoters. This allowed us to identify a transposase sequence signature associated with this unique insertion specificity and to show that the dissymmetry between the two inverted repeats is responsible for the orientation of the insertion. Surprisingly, while DDE transposases are generally associated with small and simple transposons such as insertion sequences (ISs), p-MULT encoding TEs show an unprecedented diversity with several families of IS, transposons, and ICEs ranging in size from 1.1 to 52 kb.

**Key words:** transposase, integrative and conjugative element, insertion sequence, evolution, genome dynamics.

## Introduction

Since their discovery by Barbara McClintock in the 1940s, transposable elements (TEs) have gradually attracted increasing interest. TEs were first thought to be potentially harmful parasitic entities and now are recognized as major contributors to genome evolution. TEs have been found in nearly all sequenced organisms where they can represent an important proportion of the host genome. A recent analysis of the annotation of 10 million protein-encoding genes in sequenced eukaryotic, archaeal, bacterial, and viral genomes and metagenomes revealed that transposases are the most abundant and the most ubiquitous genes in nature (Aziz et al. 2010).

Transposases, the enzymes catalyzing transposition of DNA segments, are classified in phylogenetically and structurally

unrelated families, and DDE transposases represent one of the major classes. DDE transposases show similar catalytic domain architectures with a conserved triad of essential amino acids (Asp, Asp, and Glu) which coordinates a divalent metal ion (Haren et al. 1999; Curcio and Derbyshire 2003). This catalytic domain is associated with a DNA-binding region generally located in the N-terminal part of the protein responsible for the recognition of the terminal inverted repeats (IRs) at the TE extremities which correctly localizes the transposase for strand cleavage and transfer in the transposition reaction. DDE transposases are present in all three domains of life: eukaryotes, eubacteria, and archaea. Indeed, the integrases of retroviruses and that of the *Escherichia coli* bacteriophage mu represent extensively studied classes of DDE transposases

<footnote n="1">¹Unité de Biologie des Bactéries pathogènes à Gram-positif, Institut Pasteur, Paris, France</footnote>

<footnote n="2">²UMR 3525, CNRS, Paris, France</footnote>

<footnote n="3">³Université Pierre et Marie Curie, Paris, France</footnote>

<footnote n="4">⁴Laboratoire de Microbiologie et Génétique Moléculaires, CNRS, Toulouse, France</footnote>

<footnote>*Corresponding author: E-mail: pglaser@pasteur.fr.</footnote>
(Montano et al. 2012). Although DDE transposases are struc-
turally linked, they show an overall low sequence conservation
leading to numerous unannotated or misannotated represen-
tatives in genome databases. In eukaryotes, comparative anal-
ysis of DDE transposases led to the identification of 19
superfamilies (Jurka et al. 2005). A recent phylogenetic anal-
ysis suggested that all eukaryotic cut-and-paste transposable
element superfamilies have a common evolutionary origin and
define three major phyla (Yuan and Wessler 2011). Among
these groups, the highly mutagenic Mutator and Mutator-like
elements (MULEs) represent a diverse family that is related to
the prokaryotic IS256 family (Eisen et al. 1994).

In prokaryotes, insertion sequences (ISs) are the simplest
and the most abundant autonomous TEs. They were defined
as TEs that only code the functions required for their mobility:
a transposase gene surrounded by IRs that define the borders
of the mobile DNA. IS have a dedicated repository database,
ISfinder (www-is.biotoul.fr, last accessed January 27, 2014),
that contains more than 4,000 carefully annotated ISs (Siguier
et al. 2006). Some prokaryotic TEs harbor different "passen-
ger genes," implicated in regulatory or accessory functions
such as antibiotic resistance genes which confer a selective
advantage to the host. However, the organization of these
transposons may be even more complex. We have recently
characterized a new family of TEs in streptococci, the TnGBS
family, encoding a DDE transposase associated with different
conjugative machineries that promote their horizontal transfer
(Brochet et al. 2009; Guerillot et al. 2013). This represents the
first family of integrative and conjugative elements (ICEs) in
which the phage-like integrase responsible for the excision
and integration of the element is substituted by a DDE trans-
posase. TnGBSs were shown to transpose specifically 15-
17 bp upstream different σA promoters (Brochet et al.
2009). Similarity searches of the public databases revealed IS
elements expressing transposases related to TnGBS transpo-
sases (Brochet et al. 2009). These enzymes were not related to
any known transposase.

Here, by a cascade iterated Blast search we have ex-
panded our vision of the diversity of p-MULT with the dis-
covery of four new families, in addition to the IS256 family.
TnGBS and related ISs transposases represent one of these
new families. Similarity and secondary structure predictions
allowed us to determine that these five families share an
RNase H fold and to identify, unambiguously, the catalytic
triad as well as the IR DNA-binding regions. By systematic
analysis of the genomic context we showed that these
transposases are responsible for the mobility of ISs both in
eubacteria and in archaea but also of ICEs previously de-
scribed in *Mycoplasma* (Dordet Frisoni et al. 2013). The
identification of a new family of Mutator-like elements shar-
ing the same insertion specificity as TnGBS upstream σA
promoters provides further insights into this unusual prop-
erty among prokaryotic transposases.

**Materials and Methods**

Cascade PSI-Blast Search of Transposases Related to
Gbs1118

The primary protein sequence of TnGBS2 transposase (Gbs1118) was used as an initial query in a PSI-Blast (Altschul et al. 1997) search against the NCBI nonredundant protein sequence database. Two rounds of PSI-Blast searches were performed without low complexity filter and with otherwise default parameters. Protein hits with an *E*-value above 0.005 and query coverage <60% were filtered out. Retained hits were then aligned using the MAFFT algorithm (Katoh and Standley 2013) with default parameters and a tree was built with Jalview using the average distances calculated with the BLOSUM62 matrix (Waterhouse et al. 2009). Based on this tree, a protein hit distantly related to the query was chosen to perform a second PSI-Blast search. New protein hits obtained by this second round of PSI-Blast were retained and a new query for subsequent rounds of PSI-Blast search was selected by applying the same filtering, alignment, and tree building method. The systematic propagation of PSI-Blast searches through distantly related homologs allowed us to overcome the query dependence and asymmetry of the classical use of PSI-Blast (Bhadra et al. 2006). In total, we performed seven rounds of PSI-Blast search with the following queries: Gbs1118 (NP_735564), Hore_07130 (YP_002508465), Krac_8686 (ZP_06969815), MAE_08640 (YP_001655878), MAGa5060 (YP_003515670), NAS141_01721 (ZP_00964747), and Calow_0284 (YP_004001685).

Transposase Family Clustering

Protein hits retrieved by cascade PSI-Blast were first com-
pared by all-against-all BlastP comparisons. Similarities with
an *E*-value lower than 10<sup>-4</sup> were retained to build a simi-
larity network using the Cytoscape software (Cline et al.
2007). We applied the force directed layout of Cytoscape
to visualize the generated similarity graph where each node
corresponds to a transposase homolog interconnected by
edges representing the BlastP results. We applied a contin-
uous mapping of the edge opacity to further weigh rela-
tionships between transposase homologs (the more opaque
edges correspond to lower BlastP *E*-values). Transposase ho-
mologs were then clustered by using the Markov cluster
algorithm (MCL) (http://micans.org/mcl/, last accessed
January 27, 2014) implemented in the clusterMaker plugin
of Cytoscape (Morris et al. 2011). We converted the edge
weight with -log(*E*-value) and applied an inflation factor (IF)
of 1.2. This inflation value was chosen as it has been shown
to be effective in clustering other well-defined IS families
(Siguier et al. 2009).
## Identification of Transposable Elements and of Their Insertion Sites

The identification of TEs was performed semiautomatically by using scripts written in Python programming language and using the Biopython module (www.biopython.org/, last accessed January 27, 2014). First, the DNA coding sequences of all transposase homologs were retrieved together with 400 bp of up- and downstream sequences. The extracted DNA sequences were then used as BlastN queries against the complete or draft genome sequences in which the TE is inserted. If the BlastN result gave more than three hits, the most likely TE boundaries were determined automatically based on the majority start and end of high-scoring segment pairs (hsp). BlastN results giving multiple hsp that align with the first base of the query likely correspond to larger TEs encoding several open reading frames (ORFs) upstream the transposase gene. For these transposases, the length of the surrounding DNA sequence was extended until the extremities of the TE were reached. For transposases present in less than three copies, the extracted DNA sequence was compared by BlastN with the genomic sequence of other isolates from the same species, if available. All TE boundaries were manually validated by the identification of IRs and direct repeats (DRs). All transposons and ISs identified in this study (supplementary table S1, Supplementary Material online) were submitted to the ISfinder database (Siguier et al. 2006).

Insertion sites and insertion specificity were analyzed upon extracting 300 bp sequences on both sides of the validated TEs after filtering identical insertions. These regions were scanned for putative σA promoters using the PPP software (http://bioinformatics.biol.rug.nl/websoftware/ppp/, last accessed January 27, 2014). Hidden Markov models of the lactococcal σA dependent RNA polymerase binding site, allowing a 15–19 bp distance between the canonical –35 and –10 promoter elements, were constructed using alignments of known σA binding sites (Zomer et al. 2007).

## Phylogeny and Transposase Sequences Analysis

For phylogenetic reconstruction, transposase sequences with BlastP similarity lower than 98% and representative of the diversity of TEs were retained. The transposase sequences were aligned using MAFFT version 6 with the E-INS-I method (Katoh and Standley 2013) and manually checked using Jalview (Waterhouse et al. 2009). Aligned positions with more than 60% of gaps were removed before constructing the tree. Phylogenetic relationships were inferred by Maximum likelihood (ML) using MEGA5 (Tamura et al. 2011). Prior to ML analysis, the best protein substitution model of Jones-Taylor-Thornton (JTT) was selected according to the Akaike information criterion given by the ProtTest software (Darriba et al. 2011). Branch support was

determined by 100 bootstrap replications. The level of conservation in protein sequence alignments was plotted using the plotcon application (http://emboss.sourceforge.net/, last accessed January 27, 2014). Secondary structure predictions were performed using the jpred3 server http://www.compbio.dundee.ac.uk/jpred (last accessed January 27, 2014) (Cole et al. 2008).

## Results

### Discovery of New Families of Mutator-Related Transposases

Transposases from *TnGBSs* and related ISs show a PFAM *nve* retroviral integrase domain with a low score (Brochet et al. 2008). However, we did not retrieve known transposase sequences by BlastP search at NCBI or at the ISfinder Web site. BlastP shows a low sensitivity. To retrieve more distantly related protein sequences, we performed a cascade PSI-Blast search in the nonredundant protein database. After applying the filters described in the Materials and Methods section, we retained 731 protein hits. Interestingly, although 70% of these (517) are currently described as hypothetical proteins, 104 were annotated as Mutator-like transposases, 23 as transposases of the IS256 family, 7 as *ISH6* transposases, and 80 as transposases of unknown families. This iterative search suggests that, contrary to our first analysis, *TnGBS* transposases are distantly related to known transposase families.

We then built a similarity graph to visualize the relatedness between protein sequences and to decipher the overall relationships of putative and characterized transposases. In this network, protein sequences are represented as nodes that are connected by edges weighted according to their BlastP E-Value (fig. 1A). All hits form an interlinked network in agreement with the overall relatedness of all protein hits observed in the course of the PSI-Blast analysis. In particular, it shows the relatedness of *TnGBS* transposases with the Mutator transposase superfamily and transposase of the previously identified IS256 and *ISH6* families. By using the MCL, the protein sequences of the similarity network were clustered in five groups, each defining a family of transposases that we named p-MULT 1–5 for prokaryotic Mutator-like transposase. Proteins previously identified as transposases of the IS256 and *ISH6* families are members of two different clusters, p-MULT 1 and p-MULT 2. *TnGBS* transposases are members of a large cluster encompassing 320 proteins (p-MULT 3) that are closely linked to two others clusters of 186 (p-MULT 4) and 31 (p-MULT 5) proteins, respectively (fig. 1). Strickingly, except for the putative p-MULT 5 transposases that are only connected to *TnGBS* transposases (p-MULT 3), the four other clusters are interconnected. The phylogenetic tree constructed with representatives of each family confirms this clustering (fig. 2). Our analysis extends a previous study reporting that
A

<figure><img src="image_3.png" /></figure>

B

<figure><img src="image_4.png" /></figure>

p-MULT 1

<figure><img src="image_5.png" /></figure>

p-MULT 2

<figure><img src="image_6.png" /></figure>

p-MULT 3

<figure><img src="image_7.png" /></figure>

p-MULT 4

<figure><img src="image_8.png" /></figure>

p-MULT 5

**Fig. 1.**—Similarity network and MCL clustering of prokaryotic Mutator-like transposases. (A) Similarity network of TnGBS transposases related protein sequences weighted according to all-against-all BlastP *E*-values. Each node represents a protein sequence obtained by the cascade PSI-Blast search. More opaque edges correspond to greater similarity according to the BlastP *E*-value. Each node was colored according to the MCL clustering of the network. (B) Similarity network of the five p-MULT families defined by MCL clustering of the all-against-all similarity graph of figure 1A.

ISH6 transposases are distantly related to the IS256 family (Filee et al. 2007). The discovery of three additional putative transposase families shows that, together with IS256 and ISH6 transposases, Mutator-like transposases in prokaryotes are much more diverse and widespread than previously thought.

## p-MULTs Are Encoded by Diverse Types of Mobile Elements: ISs, Transposons, and ICEs That Share Transposition Features

To systematically identify the TEs that encode transposases of the five p-MULT families, we analyzed the DNA regions on both sides of the transposase genes for IRs and direct repeats (DRs). In total, we accurately identified 424 TEs in addition to the 58 TnGBS related ICE previously described (supplementary table S1, Supplementary Material online). The 109 TEs encoding a p-MULT 1 transposase have the genetic organization of IS (fig. 3A). Based on BLAST analysis performed on the ISfinder database, they all belong to the IS256 family. These ISs are widely distributed among bacterial phyla (Proteobacteria, Firmicutes, Chlamydiae, Actinobacteria, and Deferribactere) and are also present in the archaeal phylum Euryarchaeota.

Similarly, the 11 TEs that encode ISH6 related transposases (p-MULT 2) are ISs (fig. 3B). Six are new representatives of this small group. Interestingly, although the ISH6 group was first identified in archaea of the Euryarchaeota phylum (Filee et al. 2007), we identified multiple copies of one IS belonging to this family in three different uncultured *Desulfobacterium* strains that are members of the proteobacterial phylum (ISDesp5, fig. 2).

In the three other families, we identified more complex TEs. TEs encoding p-MULT 3 transposases include the TnGBS ICEs and 168 ISs that form a new IS family that we named ISLre2 (fig. 3C). Unlike TnGBS that are restricted to the streptococcal genus, ISLre2 ISs were found in a broad variety of Firmicute species, and eight are present in multiple copies in a *Fusobacterium* and two *Synergistetes* strains, respectively, belonging to two distantly related phyla.

The 115 TEs encoding p-MULT 4 transposases were found in numerous phyla: Proteobacteria, Cyanobacteria, Nitrosopirae, Bacteroidetes, Actinobacteria, Planctomycetes, and Chloroflexi. According to the transposase phylogeny, three distinct groups of ISs or simple transposons, that we named ISAzba1, ISMich2, and ISKra4, are clearly distinguishable.
Among the *ISAzba1* group, 12 TEs forming a monophyletic branch harbor one to three ORFs in addition to the transposase gene (figs. 2 and 3D). They share an ORF similar to one of unknown function in *Agrobacterium rhizogenes* plasmid

pRiA4 (Endoh et al. 1990). Among these TEs, five carry a serine recombinase gene and one a tyrosine recombinase gene, respectively (supplementary table S1, Supplementary Material online). As shown for the *Tn3* family, these site-

<figure><img src="image_9.png" /><figcaption>**Fig. 2.**—Phylogenetic tree of prokaryotic Mutator-like transposases. Each p-MULT clade is colored according to figure 1. p-MULT 1 and p-MULT 2 transposases are encoded by ISs of the IS256 and *ISH6* families, respectively. p-MULT 3 transposases are encoded by both the *TnGBS* family and the *ISLre2* family. p-MULT 4 encoded by both transposons and by ISs form three different lineages: *ISAzba1*, *ISMich2*, and *ISKra4*. Transposons of the *ISAzba1* group encoding a pRiA4_Orf3-like protein are indicated by blue dots. IS of the *ISMich2* group with a predicted –1 frameshift in the transposase gene are indicated by pink dots. TE names are indicated at the extremity of the tree branches. TEs with a predicted σA promoter at a distance of 13–17 bp from the IR-genome junctions in more than 20% of their insertion sites (supplementary table S2, Supplementary Material online) are indicated by small black dots.</figcaption></figure>
<figure><img src="image_10.png" /><figcaption>**Fig. 3.**—Diversity of transposable elements encoding transposases of the five p-MULT families. Representation of the gene maps and structural/organizational diversity of the five TEs families: (A) p-MULT 1, IS256 family; (B) p-MULT 2, ISH6 family; (C) p-MULT 3, TnGBS/ISLre2 family; (D) p-MULT 4, ISAzba1, ISMich2, and ISKra4 families; (E) p-MULT 5 Mycoplasma ICE family. Arrows represent genes. Predicted functions of the gene products are indicated according to a color code shown at the bottom of the figure. Putative origins of replication are represented by yellow triangles.</figcaption></figure>

specific recombinases might be involved in the resolution of cointegrates generated by the transposition (Kostriken et al. 1981). Similarly, four *ISKra4* elements encode different genes of unknown function in addition to their transposases (fig. 3D and supplementary table S1, Supplementary Material online). *ISLdr1* from *Legionella drancourtii* LLAP12 encodes a group II intron putative reverse transcriptase (fig. 3D). Conversely, all *ISMich2* TEs are IS. However, in this group, the transposases are encoded by two ORFs, and a -1 frameshift is sufficient to restore the expression of a full-length transposase. As this feature is conserved in the *ISMich2* group, it likely reflects a mode of regulation of transposition by programmed translational frameshifting as described in the IS1 or IS3 families (Chandler and Fayet 1993; Nagy and Chandler 2004).

The 31 TEs encoding p-MULT 5 transposases were identified in the genomes of *Mycoplasma* species. We precisely characterized the boundaries of 21 elements. Surprisingly, they are all of large size, ranging from 7 to 37 kb (fig. 3E). Eight of these are present in two or three copies. Except for the 7-kb-long element, they all encode homologs of type IV secretion systems proteins responsible for the mobility of

conjugative elements (fig. 3E). Thirteen of these elements were previously predicted as ICE despite the absence of an identifiable integrase gene. The protein responsible for their integration was not known (Marenda et al. 2006). These ICEs share little sequence conservation and have diverse organizations, as illustrated by ten representatives depicted in figure 3E. Strikingly, the most unifying feature of these elements is the conservation of a p-MULT 5 putative transposase gene upstream of one of the two IR. This strongly suggests a role for this transposase in *Mycoplasma* ICE mobility. This prediction was recently experimentally demonstrated for ICEA of *Mycoplasma agalactiae* 5632 (Dordet Frisoni et al. 2013). ICEA is transferable by conjugation and its excision and integration involve a p-MULT 5 transposase encoded by CDS22. These ICEs, like the TnGBSs, therefore rely on a DDE transposase of the Mutator family and not a tyrosine or serine recombinase for their mobility. Based on the identification of the putative transposase gene, we identified five new ICEs in three *Mycoplasma hyopneumoniae* strains and one in *Mycoplasma capricolum* subsp. *capricolum* strain ATCC 27343.
Despite their genetic diversity, TEs encoding p-MULT transposases share several features. We identified 18- to 39-bp-long IRs at their extremities with a conserved terminal cytosine residue (supplementary fig. S1A, Supplementary Material online). Like eukaryotic MULEs, the insertion of these prokaryotic TEs generates DRs of 8 or 9 bp (supplementary table S1, Supplementary Material online). We have shown that TnGBSs transpose by production of an extrachromosomal circular form, which acts as a substrate for a plasmid-like replication and conjugative transfer. Similar circular forms of ISLre2 in *Lactobacillus reuteri* JCM 1112 have been detected by polymerase chain reaction (data not shown). IS256 family ISs also transpose via a circular intermediate (Loessner et al. 2002) as do *Mycoplasma* ICEA and ICEF-I. In these circular forms, the terminal IRs are separated by 6–10 bp sequences derived from one of its two flanking DNA sequences (Calcutt et al. 2002; Marenda et al. 2006; Brochet et al. 2009; Dordet Frisoni et al. 2013; Guerillot et al. 2013). These data suggest that the Mutator-like transposases of the five families catalyze transposition using a similar mechanism involving the formation of a circular intermediate.

### Insertion Specificity for Upstream Promoter Regions Is Shared by the p-MULT 3 Family and One Lineage of p-MULT 4

We performed a systematic analysis of the insertion specificity of the five p-MULT families by extracting cognate genomic DNA sequences next to the IR-right and left of each TE (supplementary table S2, Supplementary Material online). In total, we obtained 2,833 IR-genomic DNA junctions. TnGBS and the related ISs are preferentially inserted 15–17 bases upstream the –35 region of $\sigma_A$ promoters (Brochet et al. 2009; Guerillot et al. 2013). To determine whether other p-MULT families show a similar insertion specificity, we first searched for putative $\sigma_A$ promoter sequences on both sides of the TEs (supplementary table S2, Supplementary Material online). We then analyzed the position-specific enrichment of promoter detection relative to the end of the IRs. The result of this analysis is depicted in figure 4 for the insertion sites of TEs encoding p-MULT 1, 3, and 4. For the p-MULT 3 family, a strong relative increase of promoter detection is observed at a distance of 16 bp from IR-right (fig. 4A), whereas no enrichment at any specific position was detected for TEs encoding p-MULT 1 transposases (fig. 4C). This confirms the oriented insertion at a fixed distance from $\sigma_A$ promoters catalyzed by p-MULT 3. More interestingly, we detected a similar but lower signal for TEs encoding p-MULT 4 transposases (fig. 4B). Therefore, some p-MULT 4 transposases share the p-MULT 3 insertion specificity for upstream promoter regions (supplementary tables S1 and S2, Supplementary Material online). However, $\sigma_A$ promoters were predicted on both sides of the element. These 13 ISs (ISKra2, 4, 5 and 6; ISCasp2 and 3, ISLfe1 and 2; ISHhy1; ISAcce1; ISTvi1; ISUncu20; and

<figure><img src="image_11.png" /></figure>

<figure><img src="image_12.png" /></figure>

<figure><img src="image_13.png" /><figcaption>**Fig. 4.**—Relative position of putative $\sigma_A$ promoters identified in the DNA region surrounding TE insertion sites. The DNA sequences on both sides of insertions were extracted and scanned for putative $\sigma_A$ promoters using the PPP software (Zomer et al. 2007). The histogram represents the ratio of the predicted $\sigma_A$ promoters at a given position from the insertion site to the total number of $\sigma_A$ promoters predicted at a maximum distance of 30 bp. Only the results obtained for p-MULT 3 (A), p-MULT 4 (B), and p-MULT 1 (C) for which a sufficient number of promoters were predicted are represented (4,207, 523, and 1,145 $\sigma_A$ promoters, respectively). The abscise numbers correspond to the position of the predicted –35 sequence of $\sigma_A$ promoters relative to the insertion site. Negative values correspond to the IRI-genome junctions and positive values to IRr-genome junctions.</figcaption></figure>

$ISDesp4$ are characterized by an almost perfect complementarity (83–100%) between the right and left IRs. Likewise, ISLbu1 encoding a p-MULT 3 transposase of *Leptotrichia buccalis* shows perfectly complementary IR of 24 bp and was the only TE of this family found inserted in both orientations with respect to $\sigma_A$ promoters (supplementary table S2, Supplementary Material online). These observations suggest that the orientation of TE insertions relies on a differential recognition of the two IRs by the transposase during integration.

We did not observe any conserved DNA motif in the DNA flanking the insertions of TEs encoding p-MULT 1 or p-MULT 5 transposases, in agreement with previous observations
showing that insertion of IS256 and *M. agalactiae* ICEA is likely random (Ziebuhr et al. 1999; Dordet Frisoni et al. 2013). For the TEs encoding p-MULT 2 transposases, we observed a conservation of the flanking region among the 39 IS-genome junctions extracted. The consensus sequence of the DR corresponds to the AT-rich motif AANATNTT (supplementary fig. S1B, Supplementary Material online). Remarkably, we observed this sequence also at the insertion sites of the distantly related IS identified in the uncultured *Desulfobacterium* sp. (ISDesp5). The conservation of this particular targeting further supports the grouping of these bacterial ISs as new members of the ISH6 family.

## Conservation of a Mutator-Like Catalytic Domain

The catalytic triad of DDE transposases consists of two aspartyl (D) residues and a glutamyl (E) residue, located in a conserved core that forms a characteristic RNase H-like fold of mixed

α-helices and β-strands (see supplementary fig. S2A, Supplementary Material online) (Hickman et al. 2010). The first D residue is located in β1, the second D residue is in or just after β4, and the third D/E residue in or just before α4. These three catalytic residues were experimentally confirmed in the IS256 transposase (Loessner et al. 2002). To identify the catalytic residues in the five p-MULT families, we combined sequence alignments and secondary structure predictions. First, we observed that the sequence of the transposases from the five p-MULT families align perfectly at the three catalytic residues of the IS256 transposase (fig. 5). Second, secondary structure modeling of representatives of the five transposase families unraveled an RNase H-like fold with the expected positioning of the conserved DDE residues (supplementary fig. S2A, Supplementary Material online), despite divergences in the regions between these putative catalytic residues. Compared with a typical RNase H fold, some DDE transposase catalytic domains are characterized by the

<figure><img src="image_14.png" /><figcaption>**Fig. 5.**—Alignment of the protein domains encompassing the catalytic DDE residues in p-MULT. Transposase sequences were aligned by the MAFFT alignment software (Katoh and Standley 2013) and visualized using Jalview (Waterhouse et al. 2009). The alignment was filtered for redundancy to subsequently retain a subset of transposases for each p-MULT family representative of their diversity. Only regions surrounding the predicted DDE residues and the C/D(2)H motif were kept in the alignment. Numbers given in parentheses correspond to the distance in aa residues between the different motifs. Transposases accession numbers are indicated on the left.</figcaption></figure>
presence of a β-strand or a α-helical insert between the second D residue and the E residue (Hickman et al. 2010; Yuan and Wessler 2011). For the five p-MULT families, a 99- to 138-aa-long α-helical insert was predicted between the catalytic residues D2 and E, like previously shown in eukaryotic Mutator transposases (Hua-Van and Capy 2008; Yuan and Wessler 2011). Altogether these results allowed us to predict with a high confidence the three catalytic residues in the transposases of the five p-MULT families (fig. 5).

In addition to the predicted catalytic DDE residues, a specific signature (C/D(2)H) is conserved in all retrieved homologs, 11–19 aa downstream D2 (fig. 5). This motif is positioned in the α-helical insert located after the predicted strand-β5 in the five p-MULT families (supplementary fig. S2A, Supplementary Material online). Although the functional role of this motif is unknown, it has also been identified with a similar relative position in the predicted α-helical insert of eukaryotic Mutator-like transposases (Yuan and Wessler 2011).

## Prediction of an N-Terminal DNA-Binding Domain Implicated in IR Recognition

Two additional domains are conserved in the N-terminal part of the alignment of the p-MULT sequences (domain N1 and N2 in fig. 6 showing the similarity along the aligned sequences). Domain N2, situated upstream the catalytic

**Fig. 6.**—Level of conservation along the alignment of p-MULT sequences. The graphical representation of the similarity scores along the aligned p-MULT sequences was plotted using the plotcon software of the EMBOSS package (www.ebi.ac.uk/Tools/emboss/, last accessed January 29, 2014). The similarity was calculated by moving a window of 15 aa residues along the aligned sequences. The similarity score at each position corresponds to the average of all the possible pairwise scores at that position. The pairwise scores are taken from the BLOSUM62 matrix. The average of the similarity values at each position within the window was plotted.

domain, corresponds to the DNA-binding region shown to recognize the terminal IRs of IS256 (fig. 7) (Hennig and Ziebuhr 2010). For the five transposase families, an identical secondary structure of two α helices of similar length separated by a possible turn formed by three residues were predicted in this region, and several conserved residues are identically positioned in this structure (supplementary fig. S2B, Supplementary Material online). Therefore, this region probably represents the IR binding region for the five families.

However, in 5 of the 31 putative transposases of *Mycoplasma* ICE, we also found a more divergent domain at the same position matching the HTH_23 PFAM domain (supplementary table S1, Supplementary Material online). This domain is present in transcription regulators but has also been identified in the DDE transposase of the IS630 (equivalent to the mariner family) and IS30 families (Nagy et al. 2004). This observation suggests a replacement of the N2 domain in these five transposases.

## Insertion Specificity Upstream Promoters Is Associated with a Specific Transposase Sequence Signature

p-MULT 4 transposases belong to three different phylogenetic lineages. We took advantage of the fact that all p-MULT 4 transposases that catalyze integration of their cognate element upstream of putative σA promoters belong to the *ISKra4* group (fig. 2) to search for particular motifs associated with this insertion specificity. We compared the similarity of p-MULT 3 with p-MULT 4 transposases of the *ISKra4* group or of the *ISAzba1* and *ISMich2* groups (fig. 8) and identified a single region located between the conserved domains N1 and N2 which is differentially conserved. This region contains a conserved aspartyl residue (fig. 8). Comparison with transposases from the other p-MULT families showed that this motif is conserved only between the *TnGBS*, *ISLre2*, and *ISKra4* transposases that catalyze insertion upstream σA promoters. This strongly suggests that the motif is involved in insertion specificity.

## Discussion

Mutator and Mutator-like transposases are one of the major superfamilies of transposases in eukaryotes (Yuan and Wessler 2011). They are encoded by diverse TEs present in most eukaryotic lineages, including mammals, plant, fungi, and amoeba (Hua-Van and Capy 2008). The name Mutator originates from the ability of active copies of this element to induce mutations corresponding to diverse recombination events, which were first described in plants (Bennetzen 1984; Jiang et al. 2011) and later in other organisms (Amyotte et al. 2012). Among bacterial TEs, the IS256 family was shown to be related to the eukaryotic Mutator-like elements (Mahillon and Chandler 1998; Hua-Van and Capy 2008). In this study, we further expanded the Mutator transposase superfamilies in prokaryotes by the discovery of
<figure><img src="image_15.png" /><figcaption>**Fig. 7.**—Alignment of the predicted N-terminal DNA binding domain implicated in p-MULT IR recognition. Transposases identified in this study were aligned by the MAFFT alignment software (Katoh and Standley 2013) and visualized using Jalview (Waterhouse, et al. 2009). Only the predicted N2 domain encompassing the minimum IR-binding domain identified in the IS256 transposase (Hennig and Ziebuhr 2010) was retained in the alignment. The alignment was filtered for redundancy in order to keep a subset of transposases representative of the diversity of p-MULT 1, 2, 3, and 4 families. Transposases accession numbers are indicated on the left.</figcaption></figure>

four additional families related to IS256, defining five p-MULT families (fig. 1). The majority of prokaryotic DDE transposases are associated with IS. Although the IS256 (p-MULT 1) and ISH6 (p-MULT 2) families contain only ISs, the three other Mutator-like families display a wide variety of different organizations (fig. 3). TEs encoding p-MULT 3 transposases include ISs and the diverse streptococcal ICEs of the *TnGBS* family (Guerillot et al. 2013). More interestingly, all TEs encoding p-MULT 5 transposases are *Mycoplasma* ICEs previously described or identified in this study (Marenda et al. 2006).

The only other former example of an association between a DDE transposase and a conjugation machinery is ICE6013 (Han et al. 2009; Smyth and Robinson 2009). The combination of transposition with conjugation implies recombination constraints linked to the physical separation of donor and recipient molecules. The expansion of ICE families relying on different IS related DDE transposases highlights that

transposition via a circular intermediate overcomes these constraints by generating a molecular substrate compatible with the conjugative transfer. As this mode of transposition is common to several widespread families of IS, such as IS1, IS3, IS21, and IS30 (Polard et al. 1992; Turlan and Chandler 1995; Kallastu et al. 1998; Kiss and Olasz 1999; Berger and Haas 2001), the association of transposition with conjugative transfers of DNA might be underestimated. Alternatively, the transposition process catalyzed by Mutator-like transposases might be particularly adapted to conjugative transfer explaining why they are associated with two broad families of ICEs.

TnGBSs were shown to replicate both in the donor strain following circularization and in the recipient strain upon their insertion in the chromosome. This replication is dependent on a plasmid-like replicase and promotes the transfer of the ICE (Guerillot et al. 2013). *Mycoplasma* ICEs show several features suggesting a transient replication. First, in ICEC(27343),
<figure><img src="image_16.png" /><figcaption>**FIG. 8.**—Transposase sequence signature associated with insertion specificity upstream of σA promoters. Like in figure 6, the similarity score along the aligned putative transposase sequences was plotted from the MAFFT alignment using the plotcon software of the EMBOSS package. The similarity was calculated by moving a window of 10 aa residues along the aligned sequences and using the BLOSUM62 matrix. Two similarity plots were superimposed. The gray and black curves correspond to the similarity calculated on the alignment of p-MULT 3 plus *ISKra4* transposases (p-MULT 4) and of p-MULT 3 plus *ISAzba1* and *ISMich2* transposases (p-MULT 4), respectively. The regions showing a higher conservation with *ISKra4* transposases than with *ISAzba1* and *ISMich2* transposases are colored in gray. The protein motif of the region specifically conserved between p-MULT 3 and p-MULT 4 transposases of the *ISKra4* group, and putatively involved in the targeting of σA promoters, was generated by using WebLogo (Crooks et al. 2004).</figcaption></figure>

ICEM(95010), and ICEM(GM12), a sequence of 70–286 nt located between the putative tranposase gene and the terminal IR shows more than 83% identity with the DNA region containing the putative single strand origin of replication of the plasmid pMmc-95010 (Thiaucourt et al. 2011). Second, we identified in ICEC(27343), ICEM(GM12), and ICEA-III(5632) a *parA* homolog (fig. 3E). ParA proteins are implicated in the segregation of replicating plasmids (Lutkenhaus 2012). Therefore, a transient replication might be a common feature of ICEs relying for their mobility on a Mutator-like DDE transposase.

No three-dimensional structure of a Mutator-like transposase is presently available. Nevertheless, secondary structure predictions showed that, as with other DDE transposases, the five p-MULT families show an RNase H fold organization (supplementary fig. S2, Supplementary Material online). This analysis also revealed features specifically shared between the eukaryotic Mutator-like transposases and their prokaryotic counterparts, such as a conserved C/D(2)H signature a few amino acids after the second aspartyl residue of the catalytic triad and a long α-helical insert between the second aspartyl

and the glutamyl residues (fig. 5). Shared functional features further underscore the relationships between eukaryotic and prokaryotic Mutator-like transposases. Interestingly, as for IS256, *TnGBS*, and *Mycoplasma* ICE, circular forms have been observed in eukaryotic MULE, like Mu1 and Mu1.7 of the maize (Sundaresan and Freling 1987) and in the α3 MULE of the yeast *Kluyveromyces lactis* (Barsoum et al. 2010). Thus, this mode of transposition might be a unifying feature of the Mutator superfamily of transposases.

The dramatic increase of genomic data provides opportunities to decipher the insertion specificity of transposable elements by comparing multiple insertion sites. We have characterized the diversity of insertion specificity among the five p-MULT families. The insertions of IS256 (p-MULT 1) and of *Mycoplasma* ICE (p-MULT 5) appeared to be random. We have previously shown that *TnGBS* ICEs insert specifically upstream σA promoters in a conserved orientation. We show here that this property is shared by both the *ISLre2* family encoding p-MULT 3 transposases and several members of the *ISKra4* group (p-MULT 4 family). By comparing transposase sequences from these two lineages, we identified a
conserved motif predicted to be involved in this atypical insertion specificity among prokaryotic TEs (fig. 8). By analogy with the integrase of the yeast retrotransposon Ty3, which interacts with the transcription factors TFIIB and TFIIC (Kirchner et al. 1995), we proposed that the *TnGBS* transposase interacts with a subunit of the RNA polymerase initiation complex (Brochet et al. 2009). The location of this motif just upstream the domain N2 interacting with the two IRs is compatible with such a model. It has been suggested that transposase-mediated circularization of IS256 preferentially starts with a sequence-specific first-strand cleavage at the left IS terminus (Hennig and Ziebuhr 2010). Similarly, the asymmetric transposition of IS911 was shown to be a result of differential recognition of IRr and IRI by the transposase (Rousseau et al. 2008). Based on the analysis of the orientation of insertion of the different TEs that target promoter regions, we propose that the asymmetric recognition of the two IRs is responsible for the specific orientation of most TEs encoding p-MULT 3 transposases with the IRr next to the targeted promoter sequence.

All the retrieved members of the *ISH6* (p-MULT 2) family except one were identified in archaea. Promoters in archaea are more similar to eukaryotic Pol II dependent promoters with an AT-rich TATA box-like element (Palmer and Daniels 1995). Interestingly, *ISH6* preferentially targets an AT-rich motif (AANATNTT) that is duplicated upon transposition (supplementary table S2 and fig. S1, Supplementary Material online). Thus, this insertion specificity might also lead to a preferential insertion of these ISs in promoter or intergenic regions. Interestingly, Pack-MULEs that are nonautonomous MULEs carrying fragments of cellular genes have been shown to preferentially insert into the 5' end of genes (Jiang et al. 2011). Therefore, targeting promoter regions and avoiding transposition into genes seems to be a shared strategy among Mutator-like elements to limit the fitness cost on the host cell.

In conclusion, transposable elements encoding Mutator-like transposases are much more widespread and diverse in prokaryotes than previously thought. As in eukaryotes, they represent one major superfamily of transposable elements in prokaryotes. The late discovery of the expansion of this group was probably the result of the low protein sequence conservation that was only revealed by using an extensive cascade PSI-Blast search. The comparative analysis of these elements showed both unifying features in terms of the predicted structure and transposition mechanism, but also differences in terms of insertion specificity and of organization.

## Supplementary Material

Supplementary figures S1 and S2 and tables S1 and S2 are available at *Genome Biology and Evolution* online (http://www.gbe.oxfordjournals.org/).

## Acknowledgments

The authors thank Alexandre Almeida and Pierre-Emmanuel Douarre for critical reading of the manuscript and Stéphane Descorps-Declere for his help in bioinformatics. They also thank Carmen Buchrieser, Christine Citti, Patrick Trieu-Cuot, Violette Da Cunha, and Isabelle Rosinski-Chupin for fruitful discussions. This work was supported by the French National Research Agency (grants 2010-PATH-004-02) and the LabEx project IBEID.

## Literature Cited

Altschul SF, et al. 1997. Gapped BLAST and PSI-BLAST: a new generation of protein database search programs. Nucleic Acids Res. 25: 3389-3402.

Amyotte SG, et al. 2012. Transposable elements in phytopathogenic *Verticillium* spp.: insights into genome evolution and inter- and intra-specific diversification. BMC Genomics 13:314.

Aziz RK, Breitbart M, Edwards RA. 2010. Transposases are the most abundant, most ubiquitous genes in nature. Nucleic Acids Res. 38: 4207-4217.

Barsoum E, Martinez P, Astrom SU. 2010. Alpha3, a transposable element that promotes host sexual reproduction. Genes Dev. 24:33-44.

Bennetzen JL. 1984. Transposable element Mu1 is found in multiple copies only in Robertson's Mutator maize lines. J Mol Appl Genet. 2: 519-524.

Berger B, Haas D. 2001. Transposase and cointegrase: specialized transposition proteins of the bacterial insertion sequence IS21 and related elements. Cell Mol Life Sci. 58:403-419.

Bhadra R, et al. 2006. Cascade PSI-BLAST web server: a remote homology search tool for relating protein domains. Nucleic Acids Res. 34: W143-W146.

Brochet M, Couve E, Glaser P, Guedon G, Payot S. 2008. Integrative conjugative elements and related elements are major contributors to the genome diversity of *Streptococcus agalactiae*. J Bacteriol. 190: 6913-6917.

Brochet M, et al. 2009. Atypical association of DDE transposition with conjugation specifies a new family of mobile elements. Mol Microbiol. 71:948-959.

Calcutt MJ, Lewis MS, Wise KS. 2002. Molecular genetic analysis of ICEF, an integrative conjugal element that is present as a repetitive sequence in the chromosome of *Mycoplasma fermentans* PG18. J Bacteriol. 184: 6929-6941.

Chandler M, Fayet O. 1993. Translational frameshifting in the control of transposition in bacteria. Mol Microbiol. 7:497-503.

Cline MS, et al. 2007. Integration of biological networks and gene expression data using Cytoscape. Nat Protoc. 2:2366-2382.

Cole C, Barber JD, Barton GJ. 2008. The Jpred 3 secondary structure prediction server. Nucleic Acids Res. 36:W197-W201.

Crooks GE, Hon G, Chandonia JM, Brenner SE. 2004. WebLogo: a sequence logo generator. Genome Res. 14:1188-1190.

Curcio MJ, Derbyshire KM. 2003. The outs and ins of transposition: from Mu to kangaroo. Nat Rev Mol Cell Biol. 4:865-877.

Darriba D, Taboada GL, Doallo R, Posada D. 2011. ProtTest 3: fast selection of best-fit models of protein evolution. Bioinformatics 27:1164-1165.

Dordet Frisoni E, et al. 2013. ICEA of *Mycoplasma agalactiae*: a new family of self-transmissible integrative elements that confers conjugative properties to the recipient strain. Mol Microbiol. 89:1226-1239.

Eisen JA, Benito MI, Walbot V. 1994. Sequence similarity of putative transposases links the maize Mutator autonomous element and a group of bacterial insertion sequences. Nucleic Acids Res. 22:2634-2636.
Endoh H, Hirayama T, Aoyama T, Oka A. 1990. Characterization of the *virA* gene of the agropine-type plasmid pRiA4 of *Agrobacterium rhizogenes*. FEBS Lett. 271:28-32.

Filee J, Siguier P, Chandler M. 2007. Insertion sequence diversity in archaea. Microbiol Mol Biol Rev. 71:121-157.

Guerillot R, Da Cunha V, Sauvage E, Bouchier C, Glaser P. 2013. Modular evolution of *TnGBSs*, a new family of ICEs associating IS transposition, plasmid replication and conjugation for their spreading. J Bacteriol. 195:1979-1990.

Han X, et al. 2009. Identification of a novel variant of staphylococcal cassette chromosome mec, type II.5, and its truncated form by insertion of putative conjugative transposon *Tn6012*. Antimicrob Agents Chemother. 53:2616-2619.

Haren L, Ton-Hoang B, Chandler M. 1999. Integrating DNA: transposases and retroviral integrases. Annu Rev Microbiol. 53:245-281.

Hennig S, Ziebuhr W. 2010. Characterization of the transposase encoded by IS256, the prototype of a major family of bacterial insertion sequence elements. J Bacteriol. 192:4153-4163.

Hickman AB, Chandler M, Dyda F. 2010. Integrating prokaryotes and eukaryotes: DNA transposases in light of structure. Crit Rev Biochem Mol Biol. 45:50-69.

Hua-Van A, Capy P. 2008. Analysis of the DDE motif in the Mutator superfamily. J Mol Evol. 67:670-681.

Jiang N, Ferguson AA, Slotkin RK, Lisch D. 2011. Pack-Mutator-like transposable elements (Pack-MULEs) induce directional modification of genes through biased insertion and DNA acquisition. Proc Natl Acad Sci U S A. 108:1537-1542.

Jurka J, et al. 2005. Repbase Update, a database of eukaryotic repetitive elements. Cytogenet Genome Res. 110:462-467.

Kallastu A, Horak R, Kivisaar M. 1998. Identification and characterization of *IS1411*, a new insertion sequence which causes transcriptional activation of the phenol degradation genes in *Pseudomonas putida*. J Bacteriol. 180:5306-5312.

Katoh K, Standley DM. 2013. MAFFT Multiple Sequence Alignment Software Version 7: improvements in performance and usability. Mol Biol Evol. 30:772-780.

Kirchner J, Connolly CM, Sandmeyer SB. 1995. Requirement of RNA polymerase III transcription factors for in vitro position-specific integration of a retroviruslike element. Science 267:1488-1491.

Kiss J, Olasz F. 1999. Formation and transposition of the covalently closed *IS30* circle: the relation between tandem dimers and monomeric circles. Mol Microbiol. 34:37-52.

Kostriken R, Morita C, Heffron F. 1981. Transposon *Tn3* encodes a site-specific recombination system: identification of essential sequences, genes, and actual site of recombination. Proc Natl Acad Sci U S A. 78: 4041-4045.

Loessner I, Dietrich K, Dittrich D, Hacker J, Ziebuhr W. 2002. Transposase-dependent formation of circular IS256 derivatives in *Staphylococcus epidermidis* and *Staphylococcus aureus*. J Bacteriol. 184:4709-4714.

Lutkenhaus J. 2012. The ParA/MinD family puts things in their place. Trends Microbiol. 20:411-418.

Mahillon J, Chandler M. 1998. Insertion sequences. Microbiol Mol Biol Rev. 62:725-774.

Marenda M, et al. 2006. A new integrative conjugative element occurs in *Mycoplasma agalactiae* as chromosomal and free circular forms. J Bacteriol. 188:4137-4141.

Montano SP, Pigli YZ, Rice PA. 2012. The Mu transpososome structure sheds light on DDE recombinase evolution. Nature 491:413-417.

Morris JH, et al. 2011. clusterMaker: a multi-algorithm clustering plugin for Cytoscape. BMC Bioinformatics 12:436.

Nagy Z, Chandler M. 2004. Regulation of transposition in bacteria. Res Microbiol. 155: 387-398.

Nagy Z, Szabo M, Chandler M, Olasz F. 2004. Analysis of the N-terminal DNA binding domain of the *IS30* transposase. Mol Microbiol. 54: 478-488.

Palmer JR, Daniels CJ. 1995. *In vivo* definition of an archaeal promoter. J Bacteriol. 177:1844-1849.

Polard P, Prere MF, Fayet O, Chandler M. 1992. Transposase-induced excision and circularization of the bacterial insertion sequence *IS911*. EMBO J. 11:5079-5090.

Rousseau P, Loot C, Turlan C, Nolivos S, Chandler M. 2008. Bias between the left and right inverted repeats during *IS911* targeted insertion. J Bacteriol. 190:6111-6118.

Siguier P, Filee J, Chandler M. 2006. Insertion sequences in prokaryotic genomes. Curr Opin Microbiol. 9:526-531.

Siguier P, Gagnevin L, Chandler M. 2009. The new *IS1595* family, its relation to *IS1* and the frontier between insertion sequences and transposons. Res Microbiol. 160:232-241.

Smyth DS, Robinson DA. 2009. Integrative and sequence characteristics of a novel genetic element, *ICE6013*, in *Staphylococcus aureus*. J Bacteriol. 191:5964-5975.

Sundaresan V, Freeling M. 1987. An extrachromosomal form of the Mu transposons of maize. Proc Natl Acad Sci U S A. 84:4924-4928.

Tamura K, et al. 2011. MEGA5: molecular evolutionary genetics analysis using maximum likelihood, evolutionary distance, and maximum parsimony methods. Mol Biol Evol. 28:2731-2739.

Thiaucourt F, et al. 2011. *Mycoplasma mycoides*, from "mycoides Small Colony" to "Capri". A microevolutionary perspective. BMC Genomics 12:114.

Turlan C, Chandler M. 1995. *IS1*-mediated intramolecular rearrangements: formation of excised transposon circles and replicative deletions. EMBO J. 14:5410-5421.

Waterhouse AM, Procter JB, Martin DM, Clamp M, Barton GJ. 2009. Jalview Version 2—a multiple sequence alignment editor and analysis workbench. Bioinformatics 25:1189-1191.

Yuan YW, Wessler SR. 2011. The catalytic domain of all eukaryotic cut-and-paste transposase superfamilies. Proc Natl Acad Sci U S A. 108: 7884-7889.

Ziebuhr W, et al. 1999. A novel mechanism of phase variation of virulence in *Staphylococcus epidermidis*: evidence for control of the polysaccharide intercellular adhesin synthesis by alternating insertion and excision of the insertion sequence element IS256. Mol Microbiol. 32:345-356.

Zomer AL, Buist G, Larsen R, Kok J, Kuipers OP. 2007. Time-resolved determination of the CcpA regulon of *Lactococcus lactis* subsp. cremoris MG1363. J Bacteriol. 189:1366-1381.

Associate editor: Tal Dagan