Insertion Sequence IS26 Reorganizes Plasmids in Clinically Isolated Multidrug-Resistant Bacteria by Replicative Transposition Susu He,a Alison Burgess Hickman,a Alessandro M. Varani,b Patricia Siguier,c Michael Chandler,c John P. Dekker,d Fred Dydaa Laboratory of Molecular Biology, National Institute of Diabetes and Digestive and Kidney Diseases, National Institutes of Health, Bethesda, Maryland, USAa; Departamento de Tecnologia, Faculdade de Ciências Agrárias e Veterinárias de Jaboticabal, Universidade Estadual Paulista, Jaboticabal, São Paulo, Brazilb; Laboratoire de Microbiologie et Génétique Moléculaires, Centre national de la recherche scientifique, Toulouse, Francec; Department of Laboratory Medicine, Clinical Center, Microbiology Service, National Institutes of Health, Bethesda, Maryland, USAd ABSTRACT Carbapenemase-producing Enterobacteriaceae (CPE), which are resistant to most or all known antibiotics, consti- tute a global threat to public health. Transposable elements are often associated with antibiotic resistance determinants, suggest- ing a role in the emergence of resistance. One insertion sequence, IS26, is frequently associated with resistance determinants, but its role remains unclear. We have analyzed the genomic contexts of 70 IS26 copies in several clinical and surveillance CPE iso- lates from the National Institutes of Health Clinical Center. We used target site duplications and their patterns as guides and found that a large fraction of plasmid reorganizations result from IS26 replicative transpositions, including replicon fusions, DNA inversions, and deletions. Replicative transposition could also be inferred for transposon Tn4401, which harbors the car- bapenemase blaKPC gene. Thus, replicative transposition is important in the ongoing reorganization of plasmids carrying multidrug-resistant determinants, an observation that carries substantial clinical and epidemiological implications for under- standing how such extreme drug resistance phenotypes evolve. IMPORTANCE Although IS26 is frequently reported to reside in resistance plasmids of clinical isolates, the characteristic hall- mark of transposition, target site duplication (TSD), is generally not observed, raising questions about the mode of transposition for IS26. The previous observation of cointegrate formation during transposition implies that IS26 transposes via a replicative mechanism. The other possible outcome of replicative transposition is DNA inversion or deletion, when transposition occurs intramolecularly, and this would also generate a specific TSD pattern that might also serve as supporting evidence for the trans- position mechanism. The numerous examples we present here demonstrate that replicative transposition, used by many mobile elements (including IS26 and Tn4401), is prevalent in the plasmids of clinical isolates and results in significant plasmid reorgani- zation. This study also provides a method to trace the evolution of resistance plasmids based on TSD patterns. Received 5 May 2015 Accepted 11 May 2015 Published 9 June 2015 Citation He S, Hickman AB, Varani AM, Siguier P, Chandler M, Dekker JP, Dyda F. 2015. Insertion Sequence IS26 reorganizes plasmids in clinically isolated multidrug-resistant bacteria by replicative transposition. mBio 6(3):e00762-15. doi:10.1128/mBio.00762-15. Editor Julian E. Davies, University of British Columbia Copyright © 2015 He et al. This is an open-access article distributed under the terms of the Creative Commons Attribution-Noncommercial-ShareAlike 3.0 Unported license, which permits unrestricted noncommercial use, distribution, and reproduction in any medium, provided the original author and source are credited. Address correspondence to Fred Dyda, fred.dyda@nih.gov. This article is a direct contribution from a Fellow of the American Academy of Microbiology. Members of the family Enterobacteriaceae are major causes of severe and often lethal hospital-acquired infections. In the past two decades, carbapenemase-producing Enterobacteriaceae (CPE) have appeared that are resistant to most or all clinically available antibiotics, including carbapenems, which are often con- sidered the antibiotics of last resort (1). It has been known for over 40 years that antibiotic resistance determinants are often associ- ated with mobile genetic elements (2), a phenomenon that has been reinforced by genome sequencing of clinical isolates. For instance, blaCTX, the coding gene of an extended-spectrum �-lactamase (ESBL), is often located downstream of insertion se- quence (IS) ISEcp1 (3–6), and among carbapenemase-coding genes, blaKPC is typically carried by the Tn3 family transposon Tn4401 (7, 8), whereas blaNDM is frequently located downstream of ISAba125 (9, 10). Mobile elements are also associated with formation of resis- tance gene clusters in which different determinants that give rise to a multidrug-resistant phenotype are found in close proximity. For example, specific recombination platforms, integrons (11), can recruit cassettes carrying additional resistance genes via site- specific recombination (12–14). These are, in turn, often associ- ated with transposons and ISCRs (insertion sequences with a com- mon region), a class of mobile elements that has been implicated in capturing genes or gene fragments (15). Mobile elements can also integrate control elements such as promoters upstream of resistance determinants, leading to increased expression (16–18). The horizontal dissemination of resistance between bacteria can occur via conjugative plasmids and integrative conjugative ele- ments (ICEs) (19, 20). ISs are the simplest autonomous mobile elements in bacterial genomes. They are typically composed of short terminal se- quences, often arranged as inverted repeats (terminal inverted re- peats [TIRs]) at their ends and at least one open reading frame (ORF) that encodes the transposase, the only protein essential for RESEARCH ARTICLE crossmark May/June 2015 Volume 6 Issue 3 e00762-15 ® mbio.asm.org 1 on June 17, 2019 by guest http://m bio.asm .org/ D ow nloaded from http://creativecommons.org/licenses/by-nc-sa/3.0/ http://crossmark.crossref.org/dialog/?doi=10.1128/mBio.00762-15&domain=pdf&date_stamp=2015-6-9 mbio.asm.org http://mbio.asm.org/ mobility (21). The TIRs serve as specific recognition sites for the transposase, the enzyme that carries out the strand cleavage and transfer reactions which confer mobility. While the average IS DNA content of bacterial chromosomes is below 3%, in plas- mids larger than 20 kb it is between 5 and 15%, with one extreme case exceeding 40% (22). The 4,240 different currently identified ISs can be divided into 30 families based on their overall genetic organization and specific genetic signatures (https://www-is .biotoul.fr/) (23). ISs within one family are assumed to use the same or a similar transposition mechanism to one another, al- though only a subset of IS families has been investigated experi- mentally, and the mechanisms of many remain poorly under- stood. A vast majority of presently identified ISs encode transposases whose catalytic domain is homologous to those of RNase H and members of the retroviral integrase superfamily (24, 25), the so- called DDE transposases. Some ISs transpose by a cut-and-paste mechanism by excising from one genetic location and integrating into another (21). Others transpose by using a copy-in (26) or replicative transposition (21) pathway that requires extensive DNA replication and typically generates products known as cointegrates, which have a second copy of the IS at the target site while the original copy is left intact (27, 28). During replicative transposition, a branched or forked “Sha- piro intermediate” is formed between the donor IS and the target DNA (Fig. 1a) (28). When the IS and the target site are in two different replicons (intermolecular transposition) (Fig. 1a), sub- sequent DNA replication at the intermediate branch fuses the two replicons, duplicating both the transposable element and a short nucleotide sequence (target site duplication [TSD]) flanking the insertion site. The TSD becomes two new direct repeats (Fig. 1a, red arrows). Cointegrates can be subsequently resolved by recom- bination between the two directly repeated ISs, either by a dedi- cated resolvase or a host-mediated process, to regenerate the two replicons, each with a single IS copy, and the TSDs will flank the copy in the target replicon. When the donor IS and target sites are present in the same replicon (intramolecular transposition) (Fig. 1b), there are two possible and very different outcomes. After the initial DNA cleav- age steps, the liberated 3=-OH groups at the two IS ends can attack the target DNA in two different orientations, i.e., the 3=-OH of the “top” strand of the IS can attack either the “top” or “bottom” strand at the target site (27). Attack in one orientation (Fig. 1b, cis) results in the deletion of the DNA between the IS and the target site, leaving only one copy of the IS. A circular product is also created that contains a second IS copy, but this will survive only if it contains an origin of repli- cation. Neither of the resulting IS copies is flanked by TSDs: these are distributed between the two products of deletion. Attack in the other orientation (Fig. 1b, trans) results in IS duplication in an inverted configuration, and the DNA segment between the original IS and the target site is inverted; in addition, the target site is duplicated. However, as the TSD is at one end of the original IS and the opposite end of the copy, again, neither of the two resulting IS copies is flanked by TSDs. In principle, cis and trans attacks are equally likely, given the pseudo twofold sym- metry of double-stranded DNA (dsDNA) and the general lack of target sequence specificity of ISs. Two identical IS copies within one replicon can generate a composite transposon in which cleavage at two outer ends can mobilize both ISs and the DNA segment between them. If the IS uses a replicative transposition pathway, an alternative parsimo- nious way could be used to mobilize the composite transposon. Transposition of one of the ISs can generate a cointegrate with three IS copies (cointegrate, Fig. 1c). Resolution of this cointegrate between two IS copies (arrows, Fig. 1c) will result in the transfer of the composite transposon into the target replicon (21). This path- way has been shown experimentally in Escherichia coli for IS26 (38). In 2011, the National Institutes of Health Clinical Center (NIHCC) experienced an outbreak of carbapenem-resistant Kleb- siella pneumoniae infections initially involving 18 patients with a high fatality rate. Sequencing of all isolates recovered during the 2011 outbreak demonstrated that the outbreak was clonal (29). The NIHCC Hospital Epidemiology and Microbiology Services now routinely performs surveillance cultures of patients and the hospital environment, and other CPE isolates have since been re- covered. Here, we have analyzed CPE genomes from the NIHCC collec- tion isolated between 2011 and the end of 2013, focusing largely on IS content and particularly on replicative transposons. Al- though the possible outcomes of replicative transposition have been described (27, 28), no systematic investigation on the genomic level has been undertaken to understand its impact in clinically relevant circumstances. Our analysis of the distinct pat- terns of sequences that flank these mobile elements allows us to propose mechanisms and pathways that explain several changes due to transposition events, and we demonstrate that IS26 in par- ticular actively remodels resistance plasmids by both inter- and intramolecular replicative transposition. RESULTS IS annotations of the genomes of carbapenem-resistant Kleb- siella pneumoniae isolates from the 2011 NIH outbreak. Since the 2011 outbreak, the NIHCC has routinely cultured CPE isolates from patients and from the clinical environment. In the subse- quent 2 years, a number of CPE isolates, including Klebsiella pneu- moniae, Klebsiella oxytoca, Enterobacter cloacae complex, Citrobac- ter freundii complex, Escherichia coli, and Pantoea spp. have been isolated from these clinical and surveillance cultures. The ge- nomes of some of these strains have been sequenced at a very high coverage using Pacific Biosciences RS II single-molecule real-time (SMRT) technology which generates 14- to 40-kb-long reads (30). We therefore took advantage of these sequence data to investigate the diversity of ISs present in clinical isolates by first performing IS annotation of one representative isolate, KPNIH1 (29), from the clonal NIHCC K. pneumoniae outbreak. The assembled genome of this isolate contains one chromosome and three plasmids of �15 kb. Carbapenemase resistance is conferred by blaKPC which is carried by Tn4401 in plasmid pKpQIL, as previously ob- served in isolates from other hospitals (31, 32). Whole-genome IS annotation of KPNIH1 using ISsaga (http:// issaga.biotoul.fr/) (33) revealed 29 different ISs belonging to 12 different IS families (Fig. 2). Some of these are present either as a single copy or only in a partial form. In contrast, the ISKpn1, ISKpn26, IS903B, and IS26 representatives are present in numer- ous copies, suggesting that they are currently active. The locations of these ISs vary: all ISKpn1 and most ISKpn26 copies are located in the chromosome, whereas most IS26 copies are located in plas- mids. He et al. 2 ® mbio.asm.org May/June 2015 Volume 6 Issue 3 e00762-15 on June 17, 2019 by guest http://m bio.asm .org/ D ow nloaded from https://www-is.biotoul.fr https://www-is.biotoul.fr http://issaga.biotoul.fr http://issaga.biotoul.fr mbio.asm.org http://mbio.asm.org/ donor cointegrate target resolution a b a b a b a b ab ab a b cis trans 1 2 0 1 2 1 2 0 1 2 target 0 target 0 Shapiro- intermediate 1 2 0 0 donor 1 0 0 2 0 0 donor 1 2 1 2 1 2 0 0 0 1 0 0 2 1 2 0 0 1 0 2 0 a b cointegrates resolution target 0 donor1 2 4 3 c Composite transposon Target / TSD Flanking sequence IS 3’OH 1 2 3 4 0 0 4 3 2 1 0 0 1 4 0 2 0 3 left right donor target FIG 1 Schema of replicative transposition by an insertion sequence. (a) Intermolecular transposition. Cleavage at both TIRs of the IS results in nicks on both strands, generating 3=-OH groups which attack the target site, leading to the formation of a Shapiro intermediate. DNA replication generates the cointegrate containing a duplication of the IS and the target site. The cointegrate can be subsequently resolved into a plasmid identical to the original donor plasmid and a modified target plasmid carrying a copy of the IS flanked by TSDs arranged as direct repeats. (b) Intramolecular transposition. When an IS targets a target site in the same replicon, cleavages at both TIRs generate 3=-OH groups that can either attack the target site on the same strand (cis) or the opposite strand (trans). In the cis pathway, DNA between the IS and target site (dashed lines) becomes circularized and contains one IS copy and target site. In the trans pathway, DNA between IS and target site is instead inverted (“a b” becomes “b a”), bracketed by the original IS and a new copy in an inverted orientation. The target site is also duplicated but in inverted orientation, and each TSD is associated with one IS copy. (c) Intermolecular transposition of a composite transposon. A composite transposon consisting of two ISs (i.e., the DNA between TSDs 1 and 4 in the donor plasmid; purple) forms a cointegrate when either of the two flanking IS copies (left or right) transposes into a target plasmid (dashed lines), in which both the IS and the target site are duplicated. Resolution via homologous recombination yields a target plasmid carrying the whole composite transposon. IS, dark yellow rectangle; left and right TIRs, open and filled triangles; red arrows, targets for transposition; black arrows, potential TSDs from previous transposition events; gray oval or rectangle, origin of replication; filled green circle, 3=-OH. Arrow- heads indicate orientation, and different numbers represent different sequences. IS26 Reorganizes Plasmids by Replicative Transposition May/June 2015 Volume 6 Issue 3 e00762-15 ® mbio.asm.org 3 on June 17, 2019 by guest http://m bio.asm .org/ D ow nloaded from mbio.asm.org http://mbio.asm.org/ The four IS families with multiple copies present in KPNIH1 are all reported to generate TSDs (ISfinder tool, at https://www -is.biotoul.fr/). TSDs are genomic signatures of transposition events carried out by DDE transposases, and their length is a char- acteristic property of the transposon, reflecting the particular structure of the transpososome (the molecular assembly com- posed of transposase and transposon ends that catalyzes transpo- sition) (34). Among the IS families with multiple representatives, most chromosome-located ISs had directly abutted TSDs consis- tent with ISfinder records, yet those located on plasmids rarely had them (see Table S1 in the supplemental material). Interestingly, none of the six full IS26 copies, including two found on the chro- mosome, had directly flanking TSDs. Lack of TSDs has also been recently reported in another U.S. isolate of K. pneumoniae, ATCC BAA-2146, which carries a blaNDM carbapenemase gene and in which none of the 11 genomic IS26 copies has TSDs (35). Al- though the lack of TSDs was attributed by those authors to ho- mologous recombination following IS26 transposition—and this is certainly possible—we also suspect that, since IS26 transposes via a cointegrating mechanism, its movement might in itself gen- erate copies without flanking TSDs. Transposition mechanism of IS26. The evidence that IS26 transposes by using replicative transposition and cointegrate for- mation is compelling. IS26, a member of the IS6 insertion se- quence family, is 820 bp long and has 14 bp TIRs (36). It is a component of Tn2680, in which two IS26 copies bracket a 5-kb DNA segment, including a kanamycin resistance cassette on the Proteus vulgaris R plasmid Rtsl (36). Tn2680 transposition was observed from Rtsl into the bacteriophage P1 genome in Esche- richia coli, indicating that the two IS26 copies were responsible for mobility (37). Movement of IS26 from bacteriophage P1 to a pBR322 plasmid resulted in a cointegrate between the entire cir- cular prophage genome and the target plasmid and generated an additional IS26 copy. An 8-bp TSD was observed at the junction of IS26 and the pBR322 insertion site. RecA was required to resolve the cointegrate to generate pBR322 containing the insertion (36, 38). Another study involving IS15, the product of insertion of an IS26-like element into itself, also demonstrated that a cointegrate is the transposition product and that resolution requires RecA- initiated homologous recombination (39). IS26 in the genomes of Enterobacteriaceae strains isolated in the NIHCC from 2011 to 2013. When we extended our ISsaga analysis of NIHCC clinical strains by annotating eight additional fully assembled genomes (for a total of 36 replicons [see Table S2 in the supplemental material]), IS26 was observed to be the most predominant IS in the plasmids of these strains without any ap- parent species specificity. Although these eight genomes con- tained about 70 IS26 copies, only 6 were chromosomal. Further- more, not a single IS26 copy exhibited flanking 8-bp TSDs (Table 1). On the other hand, we observed several identical flank- ing sequences at a single end of two different IS26 copies (high- lighted with the same color label in Table 1). These were identified not only in IS copies within the same replicon but also in different replicons. These identical sequences are therefore “tracers” and allow us to follow the replicative transposition of IS26 and its consequences in the form of DNA rearrangements within these plasmids. Intermolecular transposition events in the analyzed NIH Clinical Center strains. When we examined the pairs of identical 8-bp flanking sequences of different IS26 copies in their genomic contexts, it became clear that these sequences provide evidence for replicative transposition. For example, we detected the hallmark of replicative transposition— cointegrate formation—when we compared pAAC154-a50 from the KPNIH1 strain and pKPN-819 from strain KPNIH24 from a patient hospitalized in NIHCC in July 2012. pAAC154-a50 is a small 15-kb plasmid that carries an ISSwi1 transposon derivative of the Tn3 family, ISSwi1-m2, which appears to be a descendant of ISSwi1 (Fig. 3a). ISSwi1 encodes a 0 2 4 6 8 10 12 14 IS 1R IS 1X 3 IS K pn 14 IS K pn 1 IS K pn 18 IS 14 00 IS Ec l1 IS 13 84 IS K pn 26 IS 90 3B IS Ps p1 IS A s4 IS 26 IS 61 00 IS 21 IS K pn 7 IS Sb a1 4 IS Ec 22 IS Ec 8 IS K pn 28 IS 60 9 IS K pn 6 IS A s1 2 IS K pn 25 IS Ec 38 IS Sw i1 IS Sh es 11 IS A cs p1 Tn 44 01 co py n um be r plasmid partial complete chromosome partial complete IS1IS family IS3 IS5 IS6 IS21 IS66 ISL3 Tn3 IS 41 8 IS 20 0/ IS 60 5 IS 11 82 IS A s1 FIG 2 IS distribution in Klebsiella pneumoniae strain KPNIH1. IS copy numbers in the three plasmids and single chromosome are shown in red and blue, respectively. Pale color, partial IS; dark color, complete IS. IS family classifications (in green) were determined using ISfinder (https://www-is.biotoul.fr/). He et al. 4 ® mbio.asm.org May/June 2015 Volume 6 Issue 3 e00762-15 on June 17, 2019 by guest http://m bio.asm .org/ D ow nloaded from https://www-is.biotoul.fr/ https://www-is.biotoul.fr/ https://www-is.biotoul.fr/ mbio.asm.org http://mbio.asm.org/ transposase (tnpA), a resolvase (tnpR), and a TEM beta-lactamase (blaTEM-1) and is present in plasmids pZM-3 from Salmonella en- terica serovar Wien (40) and pNJST258N4 from K. pneumoniae NJST258_1 (41). ISSwi1-m2 most likely arose from another ISSwi1 derivative, ISSwi1-m1 (also called Tn1331 [42–44]), which is present in plasmid pNJST258N5 of K. pneumoniae NJST258_1 due to an IS26 insertion into blaTEM-1 of parental ISSwi1 and ad- jacent DNA deletion. We noted that pAAC154-a50 is highly similar to a portion of pKPN-819 (Fig. 3b; see also Fig. S1a and b for further details). The indication that IS26 is involved is that the sequence in pKPN-819 homologous to pAAC154-a50 is bordered by two IS26 copies in the same orientation and is flanked by the same 8 bp at one side of each of the copies (labeled “0” in Fig. 3b). The most parsimonious explanation for the homology is that an IS26 in the ancestor of pKPN-819 inserted into the target sequence “0” in pAAC154-a50 to generate an additional IS26 copy and TSDs (Fig. 3b). Therefore, pKPN-819 appears to be an IS26-mediated cointegrate of two plasmids. There are three IS26 copies in pKPN-819, and it appears that two of these form an active composite transposon (Fig. 3b, high- lighted in purple) able to insert into the K. pneumoniae genome. The IS26 copies bracket a 5-kb DNA segment, and the proposed composite transposon can be perfectly aligned to a chromosomal locus in KPNIH24 (see Fig. S1c in the supplemental material). Furthermore, the 5-kb DNA segment is duplicated at the chromo- TABLE 1 The 8-bp flanking sequences of IS26 in nine genomes of Enterobacteriaceaea Strain Replicon Le Right Strain Replicon Le Right KPNIH1 chromosome GGTTTTCG TTCTACGG KONIH1 pKOX-137 GTCCGGAG CTCGAAAA GAGATTGG GGTTTTCG pKOX-86d GAGATTGG TAAAATAA pKpQIL CAAATTCG AGGCAGGC ATTGCTGA ATCTTTCT CGCCGATG pKPC-727 AAGCAGAC AGAAAATA pAAC154-a50 TTGGAAAC CTGATAAA AGAAAATA AGCTGCGG pKPN-498 CATCAGAG GGATTGAA GTTTCAGG GGAAGGAA AATTTATG ATTGTTTT GGTCTTAA ECNIH2 chromosome CGGCAGAG CCGCTGTA pKEC-39c GTGGACTG CGGTAGAT KPNIH24 chromosome GGTTCGCG GGTTTTCG CGGTAGAT GTTTAACC GGTTCGCG CTGATAAA GAGATTGG TAAAATAA GAGATTGG CTGATAAA ATTGCTGA TTCTACGG GGTTTTCG TTCTACGG GCTGCGCT CCCTGTAG pKPC-484 ATTGTTTT GGATTGAA GAAAAAAC ATATCCTG CAAATTCG AGGCAGGC TTTTAAGC TTTCGCCC TTGGAAAC CTGATAAA GCTTCCAC GGGAGGAC CGATCAAC CATTTATG pKPC_272 GCTGTTGC ATACTGTG pKPN-819 GGTTCGCG ACAATGTG AGAGTAGA GATAGATC TTGGAAAC CTGATAAA CCGGCGTT GGTTCGCG ECNIH3 pENT-576 ATATTCCT CGAATGAC CGAATGAC AAAAAAAC KPNIH27 pKEC-dc3 GTGGACTG CGGTAGAT CATCGCCG TTCGCCCG CGGTAGAT GTTTAACC pENT-8a4 CTCATTGT GGCTCTAC GAGATTGG TAAAATAA AGAGTAGA AGTGAGTT ATTGCTGA TTCTACGG GGCTCTAC AAGCTGGC GCTGCGCT CCCTGTAG CCGGCGTT GCGCGACG GAAAAAAC ATATCCTG CTAAACAG CGAGCCCC TTTTAAGC TTTCGCCC AAAGAAAT TTCTGAAA GCTTCCAC GGGAGGAC pKPC_47e GTGGACTG ATATCCTG pKPN-068 ATTGTGAT GCTCCATT CGCTGGAC CCGGCGTT ATTGTGAT CAATTTCT CCCTGTTG ECR091 pENT-08e CTCATTGT GGCTCTAC GCTCCATT GTTCAACG pKPC-47e GTGGACTG ATATCCTG pKPN-262 CAAAAAAC AAGGCATT CGCTGGAC CTTACGCC TAAGTCAG pKPN-b0b GAAAAGCA CATTACGT CFNIH1 chromosome GTCCAGAC CTTTCCTG TATTAGAA ATTTTTCC pKEC-a3c GTGGACTG CGGTAGAT ATTTTTCC CACCGGAG CGGTAGAT GTTTAACC GAGATTGG TAAAATAA KPR0928 pKPN-294 TTGGAAAC CTGATAAA ATTGCTGA TTCTACGG pKpQIL-531 CAAATTCG AGGCAGGC GCTGCGCT CCCTGTAG CGCCGATG GAAAAAAC ATATCCTG TTTTAAGC TTTCGCCC GCTTCCAC GGGAGGAC a Identical sequences are shown in the same color pattern. A blank implies a partial IS26 copy. IS26 Reorganizes Plasmids by Replicative Transposition May/June 2015 Volume 6 Issue 3 e00762-15 ® mbio.asm.org 5 on June 17, 2019 by guest http://m bio.asm .org/ D ow nloaded from mbio.asm.org http://mbio.asm.org/ pAAC154-a50 0 pKPN-819 1 0 3 4 Composite transposition KPNIH1_chromosome Ancestor 4 integron Ancestor 1 Ancestor 2 Unequal cross-over between sister chromosomes Ancestor 3 segregants 65 7 6 6 7 765 4 4 0 0 KPNIH24_chromosome 1 2 Ancestor 3 4 2 0 5 7 65 4 0 7 65 7 7 0 Composite transposition 5kb IS26 ISSwi1-m2 Target / TSD Integron Composite transposon 1 Composite transposon 2 ISSwi1-m1 (Tn1331) tnpAtnpRaacA6blaTEM-1 aadA2blaOXA-9 tnpR_Cter pNJST258N5 pAAC154-a50 ISSwi1-m2 tnpAtnpRaacA6 blaTEM-1_Cter aadA2_Nter IS26 ISSwi1 tnpAtnpRblaTEM-1 pZM-3 pNJST258N4 IRL IRR Resistance genes ISSwi1 and derivatives Transposase (tnpA) or resolvase (tnpR) + three additional resistance genes IS26 insertion + DNA deletion IS26 2 0 pKEC-dc3 pKOX-86d 0 1 6 0 0 6 1 2 Ancestor Intermolecular Intramolecular (cis) 7 7 3 4 4 3 IS26 Target/TSD Composite transposon blaKPC-2 cis 4 27kb 74kb FIG 3 Examples of intermolecular transposition events observed within plasmids from NIHCC isolates. (a) Alignment of three ISSwi1 derivatives from different plasmids. ISSwi1 is shown in light beige, and IS26 is shown in dark yellow. (b) Intermolecular transposition of IS26, leading to replicon fusion and composite transposition. The composite IS26-IS26 transposon is highlighted with a purple background. A second composite transposon derived from ancestor 2 is highlighted with a green background. (c) IS26-mediated plasmid rearrangements. Transposition of a composite IS26-IS26 transposon (background in blue) into an ancestor plasmid to generate pKEC-dc3 is shown. pKOX-86d is derived from an intramolecular IS26 transposition (cis refers to the insertion orientation in Fig. 1b) in this ancestor. The sequences of the different flanking sequences 0 to 7 are listed in Fig. S2 in the supplemental material. Plasmid names are shown in black. Hypothetical ancestor plasmids are labeled in red. ISs are shown as rectangles (different colors represent different ISs), with triangles representing the TIRs. Full-length genes are labeled in black, and truncated genes are labeled in blue. Left TIR, open triangle; right TIR, closed triangle; integron, gray rectangle; transposase or resolvase (tnpA and tnpR), dark red arrows; resistance genes, dark orange arrows. The different flanking 8-bp sequences of IS26 are labeled with different numbers, and their sequences are listed in figures in the supplemental material. Arrowheads indicate the orientation. These symbols and color patterns are applied in all subsequent figures. He et al. 6 ® mbio.asm.org May/June 2015 Volume 6 Issue 3 e00762-15 on June 17, 2019 by guest http://m bio.asm .org/ D ow nloaded from mbio.asm.org http://mbio.asm.org/ somal locus as tandem repeats separated by an IS26 and bracketed by two additional IS26 copies. This remarkable DNA structure suggests that the chromosomal locus is likely the result of IS26 composite transposition (Fig. 1c) followed by nonreciprocal re- combination at the IS26 copies (unequal crossover [45, 46]) (see Fig. S1d). Further supporting evidence for this proposed pathway is the absence of the tandem repeat in the chromosomal locus of KPNIH1. Instead, we found an IS26 flanked by octanucleotides 6 and 7 (Fig. 3b) that flank the tandem repeats in KPNIH24. There- fore, one possibility (Fig. 3b) is that the KPNIH1 and KPNIH24 chromosomes are two segregants of an unequal crossover event that occurred at the chromosomal IS26 composite transposon and resulted in a duplication and a deletion segregant. Alternatively, recombination between two IS26s in KPNIH24 could have given rise to the pattern of IS26s seen in KPNIH1. If the composite transposon in KPNIH24 arose via the pathway in Fig. 1c, TSDs would be expected at the target replicon, whereas the flanking sequences 6 and 7 in KPNIH24 are different. How- ever, sequence 7 is observed elsewhere in the KPNIH24 and KPNIH1 chromosomes as the flanking 8-bp sequence of another IS26 upstream of an integron. This sequence pattern thus suggests a second composite transposon, embracing the original compos- ite transposon much like nested Russian dolls, transposed prior to unequal crossover. There have also been several IS26-mediated plasmid rear- rangements in the pKEC plasmid series from the NIHCC CPE isolates (Fig. 3c). Three isolates—K. pneumoniae KPNIH27, C. freundii CFNIH1, and E. cloacae ECNIH2— carry highly simi- lar pKEC series plasmids (see Fig. S2a in the supplemental mate- rial). In all three, we found a 27-kb DNA segment bracketed by two IS26 copies in the same orientation (Fig. 3c, highlighted in blue; for simplicity, only pKEC-dc3 from KPNIH27 is shown) and with identical 8-bp TSDs arranged as direct repeats at one side of each copy (“0”). This 27-kb DNA segment encodes many conju- gative transfer and plasmid stability proteins, the highly mutagenic polymerase UmuDC, as well as anti-restriction- modification proteins. This entire region is absent in the otherwise-highly homologous pKOX-86d from K. oxytocia strain KONH1 (Fig. 3c, bottom left), but the “0” 8-bp sequence is pres- ent. This suggests either intermolecular IS26 transposition result- ing in replicon fusion (pathway in Fig. 1a, if the 27-kb region contained a replication origin) or transposition by an IS26-IS26 composite transposon followed by resolution (pathway in Fig. 1c). Another difference between these two plasmids is a 74-kb region in pKEC-dc3—missing in pKOX-86d—that is next to an IS26 copy (Fig. 3c; see also Fig. S2b in the supplemental material). This is probably due to IS26-mediated deletion in the ancestral plasmid upon intramolecular transposition in cis. Intramolecular transposition leading to DNA inversions. Matched pairs of 8-bp sequences flanking different IS26 copies within the same plasmid provide evidence for intramolecular rep- licative transposition in trans, leading to DNA inversion; four ex- amples are shown in Fig. 4. Plasmid pKPN-068 from KPNIH27 (Fig. 4a, right) carries multiple antibiotic resistance genes for two beta-lactamases (blaTEM-1 and blaSHV-12), two aminoglycoside modification enzymes (ant2 and cata1), and a quinolone resis- tance gene (qnrB). This 80-kb plasmid has seven ISs, including four IS26 copies, one copy of ISEcp1 (an IS1380 family member), an ISAcsp1-like IS (a Tn3 family member), and one ISSwi1 copy that is disrupted into three segments (pA, pB, and pC). This dis- ruption pattern is likely the result of ISAcsp1 insertion and IS26- mediated intramolecular transposition giving rise to inversion. The first disruption is clearly made by insertion of the ISAcsp1-like IS (shown in green) into the ISSwi1 transposase gene (tnpA) (see Fig. S3a in the supplemental material) as indicated by the flanking 5-bp TSDs characteristic of ISAcsp1 (shown as unfilled arrows in pKPN-068). The second disruption is an inversion of pA relative to the rest of ISSwi1 and its separation from pB by ~30 kb. In pKPN-068, two copies of IS26 adjacent to the pA and pB segments are in opposite orientation and are flanked by identical 8-bp re- peats inverted relative to each other (“0” in red). This arrange- ment is consistent with intramolecular transposition in trans of a single IS26 into a “0” target site in an assumed hypothetical ances- tor at the junction of the pB and pA segments of ISSwi1, resulting in the inversion of the intervening DNA. The hypothetical ancestor plasmid shown on the left in Fig. 4a would have carried three copies of IS26, two with iden- tical 8-bp flanking DNAs (“1”) and in the same orientation. This was probably a consequence of intermolecular transposi- tion of an IS26-IS26 composite transposon (grey background) carrying the blaSHV-2, qnrB, and cata1 resistance genes, followed by resolution (pathway in Fig. 1c). These steps would have brought two antibiotic resistance determinants, blaTEM-1 and ant2, previously 30 kb apart, into the immediate neighborhood of blaSHV-12, thereby augmenting an existing resistance cluster. This illustrates an IS26-mediated mechanism of resistance cluster for- mation different from and independent of integron- or ISCR- mediated rearrangements. We also found a pattern of diagnostic repeats within the plas- mid pKPC-727 from K. oxytoca KONIH1 (Fig. 4b), where two IS26 copies located next to a split conjugal transfer operon include two 8-bp flanking sequences (“0” in red) in an inverted configu- ration. This is consistent with IS26-mediated transposition result- ing in intramolecular inversion, leading to the remodeling of the operon. We found a similar pattern of IS26 elements and flanking repeats in plasmid pENT-576 of E. cloacae ECNIH3 (Fig. 4c), where a 16-kb DNA region that included ABC transporter and urease genes was probably inverted due to intramolecular IS26 transposition into an intergenic region. Examples of IS26-mediated intramolecular plasmid rearrange- ments extend beyond the NIHCC isolates. We detected similar IS26-related patterns in a series of pKpQIL-like carbapenemase blaKPC3-carrying plasmids embedded in strains of K. pneumoniae, E. coli, and Enterobacter aerogenes isolates from New York and New Jersey hospitals (47). Those authors noted an inversion bracketed by two IS26 copies in pKpQIL-234, but only one of these IS26 copies is present in pKpQIL-10 (Fig. 4d; see also Fig. S4 in the supplemental material). This was interpreted as an insertion of a second IS26 followed by subsequent DNA inversion between two IS26 elements. However, it seems more likely that pKpQIL- 234 is a product of a one-step intramolecular transposition in pKpQIL-10 by IS26 into the target “0” octanucleotide. Those au- thors also noted another rearrangement in pKpQIL-Ec (Fig. 4d, bottom), a deletion of 8 kb of DNA starting from the right flank of ISKpn14 (an IS1 family element). We suggest that this is consistent with intramolecular replicative transposition in cis of ISKpn14 resulting in DNA deletion, as IS1 family members also appear to use replicative transposition (48). ISSwi1 is a hot spot for transposition. As shown in Fig. 1b, one possible outcome of intramolecular transposition is DNA seg- IS26 Reorganizes Plasmids by Replicative Transposition May/June 2015 Volume 6 Issue 3 e00762-15 ® mbio.asm.org 7 on June 17, 2019 by guest http://m bio.asm .org/ D ow nloaded from mbio.asm.org http://mbio.asm.org/ ment deletion without generating telltale abutted TSDs. Deletion events are difficult to identify unless replicons both before and after the transposition event have been detected. Although this is rare, deletion can be inferred when a partial gene is left behind. As over half the 8-bp flanking sequences in Table 1 are unique, we believe some of these are due to transposition-mediated deletion. Figure 5a depicts the many ways that the 16 copies of the ISSwi1 locus in different plasmids in the analyzed NIHCC genomes have pKPN-068 30kb pC pB pA cata1 qnrB 1 4 5 3 2 1 30kb pC pB pA 1 4 0 qnrB cata13 2 1 5 Resistance genes ISSwi1 IS26 ISAcsp1 ISEcp1 Ancestora Composite transposon Inverted region Ancestor 2 123kb 0 pA pB pKPC-727 2 0 23kb 1 pA pB 0 trans Target / TSD of IS26 for intramolecular transposition (trans) Other target / TSD of IS26 Conjugal transfer operon Ancestor 2 1 16kb 0 pENT-576 2 0 16kb 10 trans ISPa38_partial b c pKpQIL-Ec pKpQIL-10 22kb pKpQIL-234 2 0 22kb trans d 1 2 0 1 0 ISKpn14 Target / TSD of other ISs 1 2 trans cis cis ISEcp1 ISEcp1 FIG 4 Intramolecular transposition (trans) resulting in DNA inversion. Cases of deduced intramolecular transposition events within hospital isolated strains are shown. Inverted regions are highlighted with a pink background. The cis and trans notations indicate the insertion orientation, as shown in Fig. 1b. Different ISs are shown in different colors. (a) Intramolecular transposition leading to pKPN-068. pA, pB, and pC represent the disrupted segments of ISSwi1. Further details, including sequences of the target and the flanking sequences 1 to 5 can be found in Fig. S3a in the supplemental material. (b) Intramolecular transposition leading to pKPC-727. pA and pB represent the partial DNAs constituting the conjugal transfer operon; further details, including the sequences of the target and the flanking 1 and 2 sequences can be found in Fig. S3b. (c) Intramolecular transposition leading to pENT-576. The sequences of target and flanking sequences 1 and 2 can be found in Fig. S3c. (d) Intramolecular transposition events among a series of pKpQIL plasmids isolated from New York and New Jersey hospitals. Further details, including sequences 1 and 2, can be found in Fig. S4 in the supplemental material. He et al. 8 ® mbio.asm.org May/June 2015 Volume 6 Issue 3 e00762-15 on June 17, 2019 by guest http://m bio.asm .org/ D ow nloaded from mbio.asm.org http://mbio.asm.org/ been disrupted. ISSwi1 is also found in hospital-isolated strains ATCC BAA-2146 (35), NJST258_1, and NJST258_2 (41) and, again, most of the copies are disrupted. The disruption pattern is usually different among the different copies but, remarkably, it is always caused by other transposable elements, such as IS26 and Tn4401 (interestingly, the three isoforms of Tn4401 we observed [a, b, and d] give rise to the highest blaKPC expression among five reported Tn4401 isoforms [49]). Thus, it appears that ISSwi1 might be a hot spot targeted by these elements. The mechanism responsible for this targeting is not clear, but as ISSwi1 carries one blaTEM-1 gene and Tn4401 carries carbapenemase gene blaKPC, the result is often the clustering of these genes. When the disruptions are caused by IS26, deletions are often observed at the antibiotic resistance gene locus. We suspect that these observations indicate ongoing plasmid streamlining that continues even after the acqui- sition of resistance determinants. Tn4401 is partially deleted by IS26 in pKPC-47e and the pKEC series plasmids (Fig. 5b) and is flanked by two IS26 copies, sug- gesting the possibility of a functional composite transposon. If this were the case, this composite transposon could be a new vehicle for blaKPC dissemination. As there are many more copies of IS26 than of Tn4401, we speculate that this composite element may have higher transposition activity than the original Tn4401 and could be an example of IS26-mediated mobility enhancement of a critical resistance determinant. Intramolecular transposition by Tn4401 creates an antibi- otic resistance gene cluster. Within the NIHCC strains, replica- tive transposition is not restricted to IS26; Tn4401 also appears to tnpRblaTEM-1 tnpA pNJST258C1 pHgISEcp1 pKPC-272 pKPC-f91pKPN-068 ISAcsp1 Tn4401b tnpRblaTEM-1 tnpAtnpR_p aacA6aadA2blaOXA-9 IS26IS26 IS26 tnpR tnpAaacA6 pNJST258N4 pNJST258N5 pAAC154-a50 pKPC-484 pKPN-819 pKPN-294 pKpQIL-6e6 pKPC-484 pKpQIL-531 pNJST258N2 Tn4401d IS26 pKPC-727 Tn4401b IS26 Tn4401b pNJST258C2 blaTEM-1-p aadA2-p pKEC-dc3 pKEC-39c pKEC-a3c pKPC-47e b Tn4401b-p ISSwi1-pIS26 IS26-p Tn4401b-p IS26 ISAcsp1 IS26 ISSwi1-p Left part deletion after insertion Right part deletion after insertion Simple insertion Intramolecular transposition (trans) Inserted DNA to generate ISSwi1-m1 Deleted DNA by IS26 insertion to generate ISSwi1-m2 IS26 FIG 5 ISSwi1 derivatives and IS26-mediated gene deletions. (a) ISSwi1 serves as an insertion hot spot for IS26 and other ISs. Plasmids isolated from hospitals other than the NIHCC are boxed by rectangles in gray, and the detailed genetic context of different ISSwi1 derivatives in these plasmids is shown in Fig. S5 in the supplemental material. Three ISSwi1 derivatives and the corresponding plasmid names are shown. (Top) Intact ISSwi1; (middle) ISSwi1-m1; (bottom) ISSwi1-m2. Inserted DNA to generate ISSwi1-m1 is highlighted in light orange. Deleted DNA by IS26 insertion to generate ISSwi1-m2 is highlighted in gray. Arrows show the other disruptions, by either IS26 or other ISs in the indicated plasmids; note that the IS26 insertions often disrupt the antibiotic resistance gene locus (aminoglycoside modification genes or blaTEM-1). Complete filled arrows indicate simple insertions, and half arrows indicate that one part of ISSwi1 was deleted after IS insertion. Red empty arrows indicate disruption was caused by intramolecular transposition in trans leading to DNA inversion. (b) Tn4401 truncations by IS26s and other ISs. IS26 Reorganizes Plasmids by Replicative Transposition May/June 2015 Volume 6 Issue 3 e00762-15 ® mbio.asm.org 9 on June 17, 2019 by guest http://m bio.asm .org/ D ow nloaded from mbio.asm.org http://mbio.asm.org/ play a role in plasmid reorganization. Although Tn3 family ele- ments such as Tn4401 are known to use replicative transposition (27), they are believed to be constrained, as they often display transposition immunity (50), a regulatory mechanism preventing the insertion of a second copy into a replicon already containing a copy of the same transposable element. Consistent with this, most replicons analyzed here carry only a single copy of Tn4401. pKPC- 727 from K. oxytoca KONIH1 is an exception (Fig. 6a, right). It carries two inverted Tn4401 copies, 73 kb apart. Neither copy has a characteristic 5-bp abutted TSD (51). Interestingly, the two Tn4401 copies are located next to two partial ISSwi1 elements (pA and pB) which are inverted and flanked by the same 5-bp sequence (“0” in red), an expected length for Tn3 family transposons and present in an inverted orientation. The DNA structures in pKPC- 727 are thus consistent with a single step of Tn4401 intramolecular transposition in trans. Tn4401 in the hypothetical ancestor would have had perfect pentanucleotide TSDs (“1”), indicating a con- ventional intermolecular transposition event as shown in Fig. 1a. Although this intramolecular transposition event disrupted blaTEM-1, it duplicated blaKPC-2 and simultaneously brought sev- eral other antibiotic resistance determinants close to one blaKPC-2 copy, thereby forming a resistance cluster. Tn4401 also appears to have carried out intramolecular trans- position within a plasmid from a K. pneumoniae isolate S9 from a New York City Hospital (Fig. 6b) (8). S9 carries two Tn4401 ele- ments in inverted orientation and without detectable abutted TSDs. As group II intron fragments were detected at both ends of Tn4401, those authors hypothesized that Tn4401 duplication was due to self-splicing by the intron. However, we suggest that the duplication is much more likely due to the intramolecular trans- position of Tn4401 into a target site within the group II intron. DISCUSSION We have carried out a global analysis of bacterial insertion se- quences in the genomes of nine CPE isolates from the NIH Clin- ical Center collected from the 2011 outbreak to the end of 2013. As these genomes have been fully sequenced with long reads and with high coverage (30), we were able to precisely analyze their mobile elements using ISsaga. Analyzing multiple isolates was key to fa- cilitating the identification of DNA rearrangements within the bacterial genomes. Although the chromosomes of these strains appear relatively stable, the plasmids where most of the resistance determinants reside are dynamic and display a large degree of interstrain variation. The presence of a highly active IS26 is con- sistent with previous reports on other clinical isolates (52–54), but here we have been able to trace its impact by taking into account its a ISSwi1-m1 IS26 Inverted region Target / TSD of Tn4401 for intramolecular transposition (trans) Other TSD of Tn4401 Ancestor S9 1 0 5kb 2 pA pB 01 0 5kb pA pB 2 b Group II intron-encoded reverse transcriptase/maturase trans Tn4401 Antibiotic resistance gene 73kb pA pB 1 73kb pB pA 1 0 1 1 Ancestor pKPC-727 Tn4401b Tn4401b Tn4401b Tn4401b Tn4401b Tn4401b FIG 6 Examples of intramolecular transposition of Tn4401. (a) Deduced intramolecular transposition event leading to pKPC-727. pA and pB represent DNA segments which partially reconstitute ISSwi1-m1. Reversed regions are highlighted with a pink background. The notation trans indicates the insertion orientation shown in Fig. 1b. Further details, including the sequences 0 and 1 are shown in Fig. S6 in the supplemental material. (b) Deduced intramolecular transposition event leading to a plasmid from K. pneumoniae S9. Details, including the sequences 0, 1, and 2 are shown in Fig. S7 in the supplemental material. He et al. 10 ® mbio.asm.org May/June 2015 Volume 6 Issue 3 e00762-15 on June 17, 2019 by guest http://m bio.asm .org/ D ow nloaded from mbio.asm.org http://mbio.asm.org/ replicative transposition mechanism and using its expected TSD patterns as guides. In several cases, our findings allow us to trace not only precisely how resistance plasmids are reorganized by this mobile element, but also the temporal relationship between re- lated plasmids. High activity of IS26 in the analyzed strains is suggested by its high copy number. Although target site duplication is usually con- sidered a fingerprint of transposition, within the genomes ana- lyzed we did not find a single IS26 copy flanked by a matching pair of TSDs. However, the lack of flanking TSDs cannot, by itself, be interpreted as lack of recent IS26 activity or be arbitrarily assumed to indicate an unconventional transposition mechanism, as inter- molecular replicative transposition without resolution does not generate abutted TSDs flanking a single IS copy. Also, the various possible consequences of intramolecular transposition must be taken into account. Nevertheless, the lack of flanking TSDs does imply that either intermolecular transposition by IS26 occurs at a very low frequency or that the resolution step of this process is very inefficient. As resolution is believed to be dependent on RecA-mediated homologous recombination (39), it seems un- likely that this host machinery is limiting. It is clear that intramolecular transposition in trans that results in DNA inversion occurs fairly often, leading to large DNA rear- rangements. This pathway generates not only additional IS26 cop- ies but also an inverted DNA segment bracketed between two copies of IS26 which in principle can act as a composite trans- poson, mobilizing the DNA segment between them in subsequent transposition events. Recently, Hall and coworkers investigated IS26 transposition in RecA-deficient E. coli (55) and reported that the presence of an IS26 in a target plasmid greatly increased the frequency of addi- tional IS26 insertions into that plasmid. Using a PCR-based assay, they found replicon fusions between IS26-containing donor and target plasmids, but the resulting cointegrate did not contain an additional copy of IS26 as would be expected for replicative trans- position. Despite the absence of RecA, the observed cointegrate was structurally equivalent to the recombination product between the two IS26 copies in the donor and target plasmids. Although we cannot exclude RecA-independent mechanisms of template switching (56), the transposase of IS26 might be responsible for this homology searching and pairing. One possibility is that two IS ends from different IS copies in separate replicons form a synapse in the same transpososome (see Fig. S8 in the supplemental ma- terial), a process that resembles transposition of a composite transposon which involves TIRs from two different IS copies lo- cated in the same replicon. Strand exchange followed by replica- tion would thus result in the observed replicon fusion. Such in trans synapsis thus allows the possibility of “self-targeted” inser- tion. A similar mechanism has been invoked to explain “targeted” insertion of IS3 and IS30 family members into TIRs (57–59). This process would be expected to be independent of RecA but, in the case of the IS3 family member IS911, it is known to depend on the RecG helicase (60). However, in these cases, the strand exchange reaction is expected to produce tandem copies of the element, due to their different transposition mechanisms. IS911 is known to use a copy-and-paste mechanism (21), whereas IS30 has been pro- posed to use both direct insertion and cointegrate formation (59). It would be certainly interesting to investigate experimentally whether such events occur and how often, or if IS26 in some cir- cumstances transposes via mechanisms different from replicative transposition. As we observed, intramolecular transposition can lead to an increase in the IS26 copy number within a genome (trans [Fig. 1b]), but if the opposite strand of the target site is used (cis [Fig. 1b]), the outcome is DNA deletion. More than half the 70 IS26 copies in the nine genomes we analyzed have flanking 8-bp sequences not associated with any of the other IS26 copies. These may well be the result of intramolecular transposition, implying that it is relatively common, as might be expected if there is no strand preference during the strand transfer step of transposition. Intramolecular transposition accompanied by deletion is clearly a useful mechanism for keeping bacterial genome size under con- trol during the onslaught of exogenous DNA acquired, for in- stance, through horizontal transfer. As the horizontal acquisition of resistance determinants is clearly important, it is possible that IS26-mediated genome trimming plays a significant role in this process, enabling the maintenance of a reasonable genome size and therefore viability (of course, we do not exclude the possibility that other ISs or transposons can accomplish the same thing). Deletions could also eliminate certain antigens to allow escape from the host immunity system, as seen in the Bordetellae, where IS activity resulted in a genome size reduction concomitant with gene inactivation (61). We have also seen that cis-mediated dis- ruptions and deletions can occur within antibiotic resistance genes (Fig. 5). Most of these involve aminoglycoside modification enzymes or narrow-spectrum beta-lactamases. Perhaps this is yet another sign of genome optimization to increase fitness by elimi- nating “narrow-spectrum” resistance genes no longer needed once a broader-spectrum resistance gene in the same class has been acquired (for example, elimination of an ESBL gene once a carbapenemase gene is acquired). In addition to playing a role in the organization of resistance genes, we observed that transposition of IS26 remodeled some plasmid-borne genes involved in conjugational transfer and plas- mid stability, for instance. As this may have conflicting effects on plasmid maintenance and stability, some mechanisms would be needed to rectify conflicts. Furthermore, replicative transposition could lead to plasmid fusion or plasmid integration into the chro- mosome in the absence of homologous sequences. For example, this has been observed in Yersinia pseudotuberculosis due to the activity of the IS6 family member ISYps1 (62). Replicative transposition is not restricted to the IS6 family, and other presumed replicative mobile elements, such as Tn4401 (Tn3 family), IS903B (IS5 family), and ISKpn14 (IS1 family), also exist (27, 48, 63). In contrast to the IS26 examples identified here, we found several copies of these family members flanked by direct repeats of the expected lengths, suggesting intermolecular trans- position followed by resolution. While Tn3 family members such as Tn4401 often encode their own resolvase (64), IS903B and ISKpn14 do not (see the ISfinder database). Perhaps resolution of cointegrates generated by these elements relies on RecA-mediated homologous recombination, as suggested for IS26, or these IS el- ements may be able to transpose by simple insertion (65, 66). Results of our analysis suggest that IS26 prefers to transpose within plasmids rather than into the chromosome. This is not a general feature of ISs, as there are many other ISs present at very high copy number in the chromosomes analyzed here. For exam- ple, up to 17 copies of ISKpn26 were observed in certain chromo- somes (e.g., KONIH1 [data not shown]). As the chromosomal IS IS26 Reorganizes Plasmids by Replicative Transposition May/June 2015 Volume 6 Issue 3 e00762-15 ® mbio.asm.org 11 on June 17, 2019 by guest http://m bio.asm .org/ D ow nloaded from mbio.asm.org http://mbio.asm.org/ copies tend to have conserved flanking TSDs more frequently than those in the plasmids, it seems that most chromosomal transposi- tion events employ simple transposition mechanisms, such as cut- and-paste or copy out-and-paste (26), which might be better tol- erated, as they cause less overall disruption than replicative transposition. Although there are several chromosomal IS26 cop- ies within the nine analyzed NIHCC strains, these are exclusively located in a prophage-carrying region and may have transposed into the phage genome before its integration into the host chro- mosome. We also believe that selection bias might be involved, which would constrain the IS26 copy number in the chromosome so as to avoid potentially harmful deletions. However, we cannot exclude the possibility of replicative transposition within chromo- somes: IS903B may have deleted a 15-kb region from the KPNIH1 chromosome (region 2716602 to 2731997 in CP008827) that is rich in transporter genes and whose loss results in the absence of this region in NJST258_1 (data not shown). This is consistent with the reported intramolecular replicative transposition activity of IS903B (63). The availability of whole-genome sequencing carried out with new tools, such as Pacific Biosciences RS II SMRT technology, makes possible the unambiguous identification of mobile ele- ments, and we expect many more clinical isolates to be sequenced. Combining these sequence data with IS annotation tools like ISsaga provides robust and reliable data from which to derive mechanistic information about the activity of mobile elements in their native hosts and therefore to understand the evolution of their host genomes. This in turn may open up a new avenue to trace IS dissemination and host adaptation. In particular, follow- ing transposition events of mobile elements such as IS26 is a way to understand the evolution of resistance plasmids, revealing changes that may optimize the viability of the pathogen in a hos- pital environment. The ability to follow these events clearly estab- lishes directionality as the rearrangements generated by replica- tive transposition events are not reversible. It is not obvious why IS26, among the large variety of mobile genetic elements, gained such prevalence in the resistance plas- mids in the examined enterobacterial strains. One possibility is a founder effect, in which the first IS26 element inserted in proxim- ity to resistance genes, which allowed it a foothold with these par- ticular genes. Alternatively, there might be some currently un- known host factor involved in IS26 mobility, making it highly active and perhaps influencing the choice toward intramolecular pathways. By extending these investigations to other clinically iso- lated pathogens, we hope to uncover common features of mobile genetic element-driven resistance plasmid reorganizations and to understand their effects in a mechanistic framework. MATERIALS AND METHODS Genomic sequences of the nine NIHCC isolates as well as other clinical isolates analyzed in this paper are from the NCBI public database (see Table S2 in the supplemental material for a list of the analyzed genomes). IS annotation for these genomes was carried out with the ISsaga program (insertion sequence semi-automatic genome annotation; at http://issaga- .biotoul.fr/ISsaga/issaga_index.php) (33). SnapGene 2.5 was used to visu- alize the annotation results. Mauve 2.3.1 was used to perform comparative genome alignments. SUPPLEMENTAL MATERIAL Supplemental material for this article may be found at http://mbio.asm.org/ lookup/suppl/doi:10.1128/mBio.00762-15/-/DCSupplemental. Figure S1, EPS file, 2 MB. Figure S2, PDF file, 0.8 MB. Figure S3, EPS file, 1.5 MB. Figure S4, EPS file, 2.3 MB. Figure S5, EPS file, 1.7 MB. Figure S6, EPS file, 1.3 MB. Figure S7, EPS file, 1.8 MB. Figure S8, EPS file, 0.9 MB. Table S1, DOCX file, 0.02 MB. Table S2, DOCX file, 0.02 MB. ACKNOWLEDGMENTS This work was partially supported by the Intramural Program of the Na- tional Institute of Diabetes and Digestive and Kidney Diseases (S.H., A.B.H., and F.D.) and the NIH Clinical Center (J.P.D.). S.H., A.B.H., M.C., J.P.D., and F.D. contributed to project design and cowrote the paper. S.H. carried out genome IS annotations and data anal- ysis. P.S. performed the IS annotations of the NJST258_1 and NJST258_2 genomes. A.M.V. adapted the ISsaga annotation tool to allow batch ge- nome analysis. We declare no competing financial interests. REFERENCES 1. McKenna M. 2013. Antibiotic resistance: the last resort. Nature 499: 394 –396. http://dx.doi.org/10.1038/499394a. 2. Hedges RW, Jacob AE. 1974. Transposition of ampicillin resistance from RP4 to other replicons. Mol Gen Genet 132:31– 40. http://dx.doi.org/ 10.1007/BF00268228. 3. Kiiru J, Butaye P, Goddeeris BM, Kariuki S. 2013. Analysis for preva- lence and physical linkages amongst integrons, ISEcp1, ISCR1, Tn21 and Tn7 encountered in Escherichia coli strains from hospitalized and non- hospitalized patients in Kenya during a 19-year period (1992–2011). BMC Microbiol 13:109. http://dx.doi.org/10.1186/1471-2180-13-109. 4. Dhanji H, Doumith M, Hope R, Livermore DM, Woodford N. 2011. ISEcp1-mediated transposition of linked blaCTX-M-3 and blaTEM-1b from the IncI1 plasmid pEK204 found in clinical isolates of Escherichia coli from Belfast, UK. J Antimicrob Chemother 66:2263–2265. http://dx.doi.org/ 10.1093/jac/dkr310. 5. Wang Y, Song C, Duan G, Zhu J, Yang H, Xi Y, Fan Q. 2013. Trans- position of ISEcp1 modulates blaCTX-M-55-mediated Shigella flexneri resis- tance to cefalothin. Int J Antimicrob Agents 42:507–512. http:// dx.doi.org/10.1016/j.ijantimicag.2013.08.009. 6. Karim A, Poirel L, Nagarajan S, Nordmann P. 2001. Plasmid-mediated extended-spectrum beta-lactamase (CTX-M-3 like) from India and gene association with insertion sequence ISEcp1. FEMS Microbiol Lett 201: 237–241. http://dx.doi.org/10.1016/S0378-1097(01)00276-2. 7. Naas T, Cuzon G, Villegas MV, Lartigue MF, Quinn JP, Nordmann P. 2008. Genetic structures at the origin of acquisition of the beta-lactamase blaKPC gene. Antimicrob Agents Chemother 52:1257–1263. http:// dx.doi.org/10.1128/AAC.01451-07. 8. Gootz TD, Lescoe MK, Dib-Hajj F, Dougherty BA, He W, Della-Latta P, Huard RC. 2009. Genetic organization of transposase regions sur- rounding blaKPC carbapenemase genes on plasmids from Klebsiella strains isolated in a New York City hospital. Antimicrob Agents Chemother 53: 1998 –2004. http://dx.doi.org/10.1128/AAC.01355-08. 9. Partridge SR, Iredell JR. 2012. Genetic contexts of blaNDM-1. Antimicrob Agents Chemother 56:6065–6067. http://dx.doi.org/10.1128/AAC.00117-12. 10. Pfeifer Y, Wilharm G, Zander E, Wichelhaus TA, Göttig S, Hunfeld KP, Seifert H, Witte W, Higgins PG. 2011. Molecular characterization of blaNDM-1 in an Acinetobacter baumannii strain isolated in Germany in 2007. J Antimicrob Chemother 66:1998 –2001. http://dx.doi.org/10.1093/ jac/dkr256. 11. Mazel D. 2006. Integrons: agents of bacterial evolution. Nat Rev Micro- biol 4:608 – 620. http://dx.doi.org/10.1038/nrmicro1462. 12. Stokes HW, Hall RM. 1989. A novel family of potentially mobile DNA elements encoding site-specific gene-integration functions: integrons. Mol Microbiol 3:1669 –1683. http://dx.doi.org/10.1111/j.1365 -2958.1989.tb00153.x. 13. Roy Chowdhury P, Ingold A, Vanegas N, Martínez E, Merlino J, Merkier AK, Castro M, González Rocha G, Borthagaray G, Centrón D, Bello He et al. 12 ® mbio.asm.org May/June 2015 Volume 6 Issue 3 e00762-15 on June 17, 2019 by guest http://m bio.asm .org/ D ow nloaded from http://issaga.biotoul.fr/ISsaga/issaga_index.php http://issaga.biotoul.fr/ISsaga/issaga_index.php http://mbio.asm.org/lookup/suppl/doi:10.1128/mBio.00762-15/-/DCSupplemental http://mbio.asm.org/lookup/suppl/doi:10.1128/mBio.00762-15/-/DCSupplemental http://dx.doi.org/10.1038/499394a http://dx.doi.org/10.1007/BF00268228 http://dx.doi.org/10.1007/BF00268228 http://dx.doi.org/10.1186/1471-2180-13-109 http://dx.doi.org/10.1093/jac/dkr310 http://dx.doi.org/10.1093/jac/dkr310 http://dx.doi.org/10.1016/j.ijantimicag.2013.08.009 http://dx.doi.org/10.1016/j.ijantimicag.2013.08.009 http://dx.doi.org/10.1016/S0378-1097(01)00276-2 http://dx.doi.org/10.1128/AAC.01451-07 http://dx.doi.org/10.1128/AAC.01451-07 http://dx.doi.org/10.1128/AAC.01355-08 http://dx.doi.org/10.1128/AAC.00117-12 http://dx.doi.org/10.1093/jac/dkr256 http://dx.doi.org/10.1093/jac/dkr256 http://dx.doi.org/10.1038/nrmicro1462 http://dx.doi.org/10.1111/j.1365-2958.1989.tb00153.x http://dx.doi.org/10.1111/j.1365-2958.1989.tb00153.x mbio.asm.org http://mbio.asm.org/ Toledo H, Márquez CM, Stokes HW. 2011. Dissemination of multiple drug resistance genes by class 1 integrons in Klebsiella pneumoniae isolates from four countries: a comparative study. Antimicrob Agents Chemother 55: 3140–3149. http://dx.doi.org/10.1128/AAC.01529-10. 14. Partridge SR, Recchia GD, Scaramuzzi C, Collis CM, Stokes HW, Hall RM. 2000. Definition of the attI1 site of class 1 integrons. Microbiology 146:2855–2864. 15. Toleman MA, Bennett PM, Walsh TR. 2006. ISCR elements: novel gene- capturing systems of the 21st century? Microbiol Mol Biol Rev 70: 296 –316. http://dx.doi.org/10.1128/MMBR.00048-05. 16. Siguier P, Gourbeyre E, Chandler M. 2014. Bacterial insertion sequences: their genomic impact and diversity. FEMS Microbiol Rev 38:865– 891. http://dx.doi.org/10.1111/1574-6976.12067. 17. Olliver A, Vallé M, Chaslus-Dancla E, Cloeckaert A. 2005. Overexpres- sion of the multidrug efflux operon acrEF by insertional activation with IS1 or IS10 elements in Salmonella enterica serovar Typhimurium DT204 acrB mutants selected with fluoroquinolones. Antimicrob Agents Che- mother 49:289 –301. http://dx.doi.org/10.1128/AAC.49.1.289-301.2005. 18. Jellen-Ritter AS, Kern WV. 2001. Enhanced expression of the multidrug efflux pumps AcrAB and AcrEF associated with insertion element trans- position in Escherichia coli mutants selected with a fluoroquinolone. An- timicrob Agents Chemother 45:1467–1472. http://dx.doi.org/10.1128/ AAC.45.5.1467-1472.2001. 19. Toleman MA, Walsh TR. 2011. Combinatorial events of insertion se- quences and ICE in gram-negative bacteria. FEMS Microbiol Rev 35: 912–935. http://dx.doi.org/10.1111/j.1574-6976.2011.00294.x. 20. Schubert S, Dufke S, Sorsa J, Heesemann J. 2004. A novel integrative and conjugative element (ICE) of Escherichia coli: the putative progenitor of the Yersinia high-pathogenicity island. Mol Microbiol 51:837– 848. http:// dx.doi.org/10.1046/j.1365-2958.2003.03870.x. 21. Mahillon J, Chandler M. 1998. Insertion sequences. Microbiol Mol Biol Rev 62:725–774. 22. Siguier P, Filée J, Chandler M. 2006. Insertion sequences in prokaryotic genomes. Curr Opin Microbiol 9:526 –531. http://dx.doi.org/10.1016/ j.mib.2006.08.005. 23. Siguier P, Perochon J, Lestrade L, Mahillon J, Chandler M. 2006. ISfinder: the reference centre for bacterial insertion sequences. Nucleic Acids Res 34:D32–D36. http://dx.doi.org/10.1093/nar/gkj014. 24. Dyda F, Hickman AB, Jenkins TM, Engelman A, Craigie R, Davies DR. 1994. Crystal structure of the catalytic domain of HIV-1 integrase: simi- larity to other polynucleotidyl transferases. Science 266:1981–1986. http:// dx.doi.org/10.1126/science.7801124. 25. 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. http://dx.doi.org/10.3109/10409230903505596. 26. Curcio MJ, Derbyshire KM. 2003. The outs and ins of transposition: from Mu to kangaroo. Nat Rev Mol Cell Biol 4:865– 877. http://dx.doi.org/ 10.1038/nrm1241. 27. Arthur A, Sherratt D. 1979. Dissection of the transposition process: a transposon-encoded site-specific recombination system. Mol Gen Genet 175:267–274. http://dx.doi.org/10.1007/BF00397226. 28. Shapiro JA. 1979. Molecular model for the transposition and replication of bacteriophage Mu and other transposable elements. Proc Natl Acad Sci U S A 76:1933–1937. http://dx.doi.org/10.1073/pnas.76.4.1933. 29. Snitkin ES, Zelazny AM, Thomas PJ, Stock F, NISC Comparative Sequencing Program Group, Henderson DK, Palmore TN, Segre JA. 2012. Tracking a hospital outbreak of carbapenem-resistant Klebsiella pneumoniae with whole-genome sequencing. Sci Transl Med 4:148ra116. http://dx.doi.org/10.1126/scitranslmed.3004129. 30. Conlan S, Thomas PJ, Deming C, Park M, Lau AF, Dekker JP, Snitkin ES, Clark TA, Luong K, Song Y, Tsai YC, Boitano M, Dayal J, Brooks SY, Schmidt B, Young AC, Thomas JW, Bouffard GG, Blakesley RW, NISC Comparative Sequencing Program, Mullikin JC, Korlach J, Hen- derson DK, Frank KM, Palmore TN, Segre JA. 2014. Single-molecule sequencing to track plasmid diversity of hospital-associated carbapenemase-producing Enterobacteriaceae. Sci Transl Med 6:254ra126. http://dx.doi.org/10.1126/scitranslmed.3009845. 31. Chen L, Chavda KD, Melano RG, Jacobs MR, Koll B, Hong T, Rojtman AD, Levi MH, Bonomo RA, Kreiswirth BN. 2014. Comparative genomic analysis of KPC-encoding pKpQIL-like plasmids and their distribution in New Jersey and New York hospitals. Antimicrob Agents Chemother 58: 2871–2877. http://dx.doi.org/10.1128/AAC.00120-14. 32. Leavitt A, Chmelnitsky I, Carmeli Y, Navon-Venezia S. 2010. Complete nucleotide sequence of KPC-3-encoding plasmid pKpQIL in the epidemic Klebsiella pneumoniae sequence type 258. Antimicrob Agents Chemother 54:4493– 4496. http://dx.doi.org/10.1128/AAC.00175-10. 33. Varani AM, Siguier P, Gourbeyre E, Charneau V, Chandler M. 2011. ISsaga is an ensemble of web-based methods for high throughput identi- fication and semi-automatic annotation of insertion sequences in pro- karyotic genomes. Genome Biol 12:R30. http://dx.doi.org/10.1186/gb -2011-12-3-r30. 34. Dyda F, Chandler M, Hickman AB. 2012. The emerging diversity of transpososome architectures. Q Rev Biophys 45:493–521. http:// dx.doi.org/10.1017/S0033583512000145. 35. Hudson CM, Bent ZW, Meagher RJ, Williams KP. 2014. Resistance determinants and mobile genetic elements of an NDM-1-encoding Kleb- siella pneumoniae strain. PLoS ONE 9:e99209. http://dx.doi.org/10.1371/ journal.pone.0099209. 36. Mollet B, Iida S, Shepherd J, Arber W. 1983. Nucleotide sequence of IS26, a new prokaryotic mobile genetic element. Nucleic Acids Res 11: 6319 – 6330. http://dx.doi.org/10.1093/nar/11.18.6319. 37. Iida S, Meyer J, Linder P, Goto N, Nakaya R, Reif HJ, Arber W. 1982. The kanamycin resistance transposon Tn2680 derived from the R plasmid Rts1 and carried by phage P1Km has flanking 0.8-kb-long direct repeats. Plasmid 8:187–198. http://dx.doi.org/10.1016/0147-619X(82)90056-7. 38. Iida S, Mollet B, Meyer J, Arber W. 1984. Functional characterization of the prokaryotic mobile genetic element IS26. Mol Gen Genet 198:84 – 89. http://dx.doi.org/10.1007/BF00328705. 39. Trieu-Cuot P, Courvalin P. 1985. Transposition behavior of IS15 and its progenitor IS15-�: are cointegrates exclusive end products? Plasmid 14: 80 – 89. http://dx.doi.org/10.1016/0147-619X(85)90034-4. 40. Villa L, García-Fernández A, Fortini D, Carattoli A. 2010. Replicon sequence typing of IncF plasmids carrying virulence and resistance deter- minants. J Antimicrob Chemother 65:2518 –2529. http://dx.doi.org/ 10.1093/jac/dkq347. 41. DeLeo FR, Chen L, Porcella SF, Martens CA, Kobayashi SD, Porter AR, Chavda KD, Jacobs MR, Mathema B, Olsen RJ, Bonomo RA, Musser JM, Kreiswirth BN. 2014. Molecular dissection of the evolution of carbapenem-resistant multilocus sequence type 258 Klebsiella pneu- moniae. Proc Natl Acad Sci U S A 111:4988 – 4993. http://dx.doi.org/ 10.1073/pnas.1321364111. 42. Tolmasky ME, Crosa JH. 1987. Tn1331, a novel multiresistance trans- poson encoding resistance to amikacin and ampicillin in Klebsiella pneu- moniae. Antimicrob Agents Chemother 31:1955–1960. http://dx.doi.org/ 10.1128/AAC.31.12.1955. 43. Tolmasky ME, Crosa JH. 1993. Genetic organization of antibiotic resis- tance genes (aac(6=)-Ib, aadA, and oxa9) in the multiresistance transposon Tn1331. Plasmid 29:31– 40. http://dx.doi.org/10.1006/plas.1993.1004. 44. Chen L, Chavda KD, Al Laham N, Melano RG, Jacobs MR, Bonomo RA, Kreiswirth BN. 2013. Complete nucleotide sequence of a blaKPC- harboring IncI2 plasmid and its dissemination in New Jersey and New York hospitals. Antimicrob Agents Chemother 57:5019 –5025. http:// dx.doi.org/10.1128/AAC.01397-13. 45. Chandler M, Galas DJ. 1983. IS1-mediated tandem duplication of plas- mid pBR322. Dependence on recA and on DNA polymerase I. J Mol Biol 165:183–190. http://dx.doi.org/10.1016/S0022-2836(83)80249-6. 46. Smith GP. 1976. Evolution of repeated DNA sequences by unequal cross- over. Science 191:528 –535. http://dx.doi.org/10.1126/science.1251186. 47. Billard-Pomares T, Fouteau S, Jacquet ME, Roche D, Barbe V, Cas- tellanos M, Bouet JY, Cruveiller S, Médigue C, Blanco J, Clermont O, Denamur E, Branger C. 2014. Characterization of a P1-like bacteriophage carrying an SHV-2 extended-spectrum beta-lactamase from an Esche- richia coli strain. Antimicrob Agents Chemother 58:6550 – 6557. http:// dx.doi.org/10.1128/AAC.03183-14. 48. Turlan C, Chandler M. 1995. IS1-mediated intramolecular rearrangements: formation of excised transposon circles and replicative deletions. EMBO J 14:5410 –5421. 49. Naas T, Cuzon G, Truong HV, Nordmann P. 2012. Role of ISKpn7 and deletions in blaKPC gene expression. Antimicrob Agents Chemother 56: 4753– 4759. http://dx.doi.org/10.1128/AAC.00334-12. 50. Lambin M, Nicolas E, Oger CA, Nguyen N, Prozzi D, Hallet B. 2012. Separate structural and functional domains of Tn4430 transposase con- tribute to target immunity. Mol Microbiol 83:805– 820. http://dx.doi.org/ 10.1111/j.1365-2958.2012.07967.x. 51. Cuzon G, Naas T, Nordmann P. 2011. Functional characterization of Tn4401, a Tn3-based transposon involved in blaKPC gene mobilization. IS26 Reorganizes Plasmids by Replicative Transposition May/June 2015 Volume 6 Issue 3 e00762-15 ® mbio.asm.org 13 on June 17, 2019 by guest http://m bio.asm .org/ D ow nloaded from http://dx.doi.org/10.1128/AAC.01529-10 http://dx.doi.org/10.1128/MMBR.00048-05 http://dx.doi.org/10.1111/1574-6976.12067 http://dx.doi.org/10.1128/AAC.49.1.289-301.2005 http://dx.doi.org/10.1128/AAC.45.5.1467-1472.2001 http://dx.doi.org/10.1128/AAC.45.5.1467-1472.2001 http://dx.doi.org/10.1111/j.1574-6976.2011.00294.x http://dx.doi.org/10.1046/j.1365-2958.2003.03870.x http://dx.doi.org/10.1046/j.1365-2958.2003.03870.x http://dx.doi.org/10.1016/j.mib.2006.08.005 http://dx.doi.org/10.1016/j.mib.2006.08.005 http://dx.doi.org/10.1093/nar/gkj014 http://dx.doi.org/10.1126/science.7801124 http://dx.doi.org/10.1126/science.7801124 http://dx.doi.org/10.3109/10409230903505596 http://dx.doi.org/10.1038/nrm1241 http://dx.doi.org/10.1038/nrm1241 http://dx.doi.org/10.1007/BF00397226 http://dx.doi.org/10.1073/pnas.76.4.1933 http://dx.doi.org/10.1126/scitranslmed.3004129 http://dx.doi.org/10.1126/scitranslmed.3009845 http://dx.doi.org/10.1128/AAC.00120-14 http://dx.doi.org/10.1128/AAC.00175-10 http://dx.doi.org/10.1186/gb-2011-12-3-r30 http://dx.doi.org/10.1186/gb-2011-12-3-r30 http://dx.doi.org/10.1017/S0033583512000145 http://dx.doi.org/10.1017/S0033583512000145 http://dx.doi.org/10.1371/journal.pone.0099209 http://dx.doi.org/10.1371/journal.pone.0099209 http://dx.doi.org/10.1093/nar/11.18.6319 http://dx.doi.org/10.1016/0147-619X(82)90056-7 http://dx.doi.org/10.1007/BF00328705 http://dx.doi.org/10.1016/0147-619X(85)90034-4 http://dx.doi.org/10.1093/jac/dkq347 http://dx.doi.org/10.1093/jac/dkq347 http://dx.doi.org/10.1073/pnas.1321364111 http://dx.doi.org/10.1073/pnas.1321364111 http://dx.doi.org/10.1128/AAC.31.12.1955 http://dx.doi.org/10.1128/AAC.31.12.1955 http://dx.doi.org/10.1006/plas.1993.1004 http://dx.doi.org/10.1128/AAC.01397-13 http://dx.doi.org/10.1128/AAC.01397-13 http://dx.doi.org/10.1016/S0022-2836(83)80249-6 http://dx.doi.org/10.1126/science.1251186 http://dx.doi.org/10.1128/AAC.03183-14 http://dx.doi.org/10.1128/AAC.03183-14 http://dx.doi.org/10.1128/AAC.00334-12 http://dx.doi.org/10.1111/j.1365-2958.2012.07967.x http://dx.doi.org/10.1111/j.1365-2958.2012.07967.x mbio.asm.org http://mbio.asm.org/ Antimicrob Agents Chemother 55:5370 –5373. http://dx.doi.org/10.1128/ AAC.05202-11. 52. Bertini A, Poirel L, Bernabeu S, Fortini D, Villa L, Nordmann P, Carattoli A. 2007. Multicopy blaOXA-58 gene as a source of high-level resistance to car- bapenems in Acinetobacter baumannii. Antimicrob Agents Chemother 51: 2324–2328. http://dx.doi.org/10.1128/AAC.01502-06. 53. Lee SY, Park YJ, Yu JK, Jung S, Kim Y, Jeong SH, Arakawa Y. 2012. Prevalence of acquired fosfomycin resistance among extended-spectrum beta-lactamase-producing Escherichia coli and Klebsiella pneumoniae clin- ical isolates in Korea and IS26-composite transposon surrounding fosA3. J Antimicrob Chemother 67:2843–2847. http://dx.doi.org/10.1093/jac/ dks319. 54. Miriagou V, Carattoli A, Tzelepi E, Villa L, Tzouvelekis LS. 2005. IS26-associated In4-type integrons forming multiresistance loci in enter- obacterial plasmids. Antimicrob Agents Chemother 49:3541–3543. http:// dx.doi.org/10.1128/AAC.49.8.3541-3543.2005. 55. Harmer CJ, Moran RA, Hall RM. 2014. Movement of IS26-associated antibiotic resistance genes occurs via a translocatable unit that includes a single IS26 and preferentially inserts adjacent to another IS26. mBio 5(5): e01801-01814. http://dx.doi.org/10.1128/mBio.01801-14. 56. Dutra BE, Sutera VA, Jr, Lovett ST. 2007. RecA-independent recombi- nation is efficient but limited by exonucleases. Proc Natl Acad Sci U S A 104:216 –221. http://dx.doi.org/10.1073/pnas.0608293104. 57. Loot C, Turlan C, Chandler M. 2004. Host processing of branched DNA intermediates is involved in targeted transposition of IS911. Mol Micro- biol 51:385–393. http://dx.doi.org/10.1046/j.1365-2958.2003.03850.x. 58. Loot C, Turlan C, Rousseau P, Ton-Hoang B, Chandler M. 2002. A target specificity switch in IS911 transposition: the role of the OrfA pro- tein. EMBO J 21:4172– 4182. http://dx.doi.org/10.1093/emboj/cdf403. 59. Olasz F, Stalder R, Arber W. 1993. Formation of the tandem repeat (IS30)2 and its role in IS30-mediated transpositional DNA rearrange- ments. Mol Gen Genet 239:177–187. 60. Turlan C, Loot C, Chandler M. 2004. IS911 partial transposition prod- ucts and their processing by the Escherichia coli RecG helicase. Mol Micro- biol 53:1021–1033. http://dx.doi.org/10.1111/j.1365-2958.2004.04165.x. 61. Parkhill J, Sebaihia M, Preston A, Murphy LD, Thomson N, Harris DE, Holden MTG, Churcher CM, Bentley SD, Mungall KL, Cerdeño- Tárraga AM, Temple L, James K, Harris B, Quail MA, Achtman M, Atkin R, Baker S, Basham D, Bason N, Cherevach I, Chillingworth T, Collins M, Cronin A, Davis P, Doggett J, Feltwell T, Goble A, Hamlin N, Hauser H, Holroyd S, Jagels K, Leather S, Moule S, Norberczak H, O’Neil S, Ormond D, Price C, Rabbinowitsch E, Rutter S, Sanders M, Saunders D, Seeger K, Sharp S, Simmonds M, Skelton J, Squares R, Squares S, Stevens K, Unwin L, Whitehead S, Barrell BG, Maskell DJ. 2003. Comparative analysis of the genome sequences of Bordetella pertus- sis, Bordetella parapertussis and Bordetella bronchiseptica. Nat Genet 35: 32– 40. http://dx.doi.org/10.1038/ng1227. 62. Lesic B, Zouine M, Ducos-Galand M, Huon C, Rosso ML, Prévost MC, Mazel D, Carniel E. 2012. A natural system of chromosome transfer in Yersinia pseudotuberculosis. PLoS Genet 8:e1002529. http://dx.doi.org/ 10.1371/journal.pgen.1002529. 63. Weinert TA, Schaus NA, Grindley ND. 1983. Insertion sequence dupli- cation in transpositional recombination. Science 222:755–765. http:// dx.doi.org/10.1126/science.6314502. 64. Blake DG, Boocock MR, Sherratt DJ, Stark WM. 1995. Cooperative binding of Tn3 resolvase monomers to a functionally asymmetric binding site. Curr Biol 5:1036 –1046. http://dx.doi.org/10.1016/S0960-9822 (95)00208-9. 65. Galas DJ, Chandler M. 1982. Structure and stability of Tn9-mediated cointegrates. Evidence for two pathways of transposition. J Mol Biol 154: 245–272. http://dx.doi.org/10.1016/0022-2836(82)90063-8. 66. Grindley NDF, Joyce CM. 1981. Analysis of the structure and function of the kanamycin-resistance transposon Tn903. Cold Spring Harb Symp Quant Biol 45:125–133. http://dx.doi.org/10.1101/SQB.1981.045.01.021. He et al. 14 ® mbio.asm.org May/June 2015 Volume 6 Issue 3 e00762-15 on June 17, 2019 by guest http://m bio.asm .org/ D ow nloaded from http://dx.doi.org/10.1128/AAC.05202-11 http://dx.doi.org/10.1128/AAC.05202-11 http://dx.doi.org/10.1128/AAC.01502-06 http://dx.doi.org/10.1093/jac/dks319 http://dx.doi.org/10.1093/jac/dks319 http://dx.doi.org/10.1128/AAC.49.8.3541-3543.2005 http://dx.doi.org/10.1128/AAC.49.8.3541-3543.2005 http://dx.doi.org/10.1128/mBio.01801-14 http://dx.doi.org/10.1073/pnas.0608293104 http://dx.doi.org/10.1046/j.1365-2958.2003.03850.x http://dx.doi.org/10.1093/emboj/cdf403 http://dx.doi.org/10.1111/j.1365-2958.2004.04165.x http://dx.doi.org/10.1038/ng1227 http://dx.doi.org/10.1371/journal.pgen.1002529 http://dx.doi.org/10.1371/journal.pgen.1002529 http://dx.doi.org/10.1126/science.6314502 http://dx.doi.org/10.1126/science.6314502 http://dx.doi.org/10.1016/S0960-9822(95)00208-9 http://dx.doi.org/10.1016/S0960-9822(95)00208-9 http://dx.doi.org/10.1016/0022-2836(82)90063-8 http://dx.doi.org/10.1101/SQB.1981.045.01.021 mbio.asm.org http://mbio.asm.org/ RESULTS IS annotations of the genomes of carbapenem-resistant Klebsiella pneumoniae isolates from the 2011 NIH outbreak. Transposition mechanism of IS26. IS26 in the genomes of Enterobacteriaceae strains isolated in the NIHCC from 2011 to 2013. Intermolecular transposition events in the analyzed NIH Clinical Center strains. Intramolecular transposition leading to DNA inversions. ISSwi1 is a hot spot for transposition. Intramolecular transposition by Tn4401 creates an antibiotic resistance gene cluster. DISCUSSION MATERIALS AND METHODS SUPPLEMENTAL MATERIAL ACKNOWLEDGMENTS REFERENCES