Skip to main content
  • Research article
  • Open access
  • Published:

Genome plasticity favours double chromosomal Tn4401b-blaKPC-2 transposon insertion in the Pseudomonas aeruginosa ST235 clone

Abstract

Background

Pseudomonas aeruginosa Sequence Type 235 is a clone that possesses an extraordinary ability to acquire mobile genetic elements and has been associated with the spread of resistance genes, including genes that encode for carbapenemases. Here, we aim to characterize the genetic platforms involved in resistance dissemination in blaKPC-2-positive P. aeruginosa ST235 in Colombia.

Results

In a prospective surveillance study of infections in adult patients attended in five ICUs in five distant cities in Colombia, 58 isolates of P. aeruginosa were recovered, of which, 27 (46.6%) were resistant to carbapenems. The molecular analysis showed that 6 (22.2%) and 4 (14.8%) isolates harboured the blaVIM and blaKPC-2 genes, respectively. The four blaKPC-2-positive isolates showed a similar PFGE pulsotype and belonged to ST235. Complete genome sequencing of a representative ST235 isolate shows a unique chromosomal contig of 7097.241 bp with eight different resistance genes identified and five transposons: a Tn6162-like with ant(2″)-Ia, two Tn402-like with ant(3″)-Ia and blaOXA-2 and two Tn4401b with blaKPC-2. All transposons were inserted into the genomic islands. Interestingly, the two Tn4401b copies harbouring blaKPC-2 were adjacently inserted into a new genomic island (PAGI-17) with traces of a replicative transposition process. This double insertion was probably driven by several structural changes within the chromosomal region containing PAGI-17 in the ST235 background.

Conclusion

This is the first report of a double Tn4401b chromosomal insertion in P. aeruginosa, just within a new genomic island (PAGI-17). This finding indicates once again the great genomic plasticity of this microorganism.

Background

Pseudomonas aeruginosa is an opportunistic pathogen that causes numerous infections in hospitalized patients, especially in the intensive care unit (ICU), some of which are difficult to treat because of their multidrug-resistant (MDR) phenotype. In P. aeruginosa carbapenem resistance has been associated mainly with two different mechanisms: first by preventing the access of the antibiotic to its active site by means of chromosomal mutations that upregulate of multidrug efflux systems or alterations in the outer membrane permeability and second by the acquisition and production of carbapenemases [1]. In the last case, the most frequent carbapenemases found in P. aeruginosa are GES, IMP, VIM, SPM and more recently KPC and NDM [1,2,3,4,5,6].

The KPC enzyme is the most frequent class A carbapenemase and has been reported in countries from five continents. KPC has been found mainly in the epidemic Klebsiella pneumoniae ST258 clone and in some Enterobacteriaceae [7]; however, it has also been identified in P. aeruginosa, first in Colombia [2], then in Trinidad and Tobago [8], the United States of America [3], Puerto Rico [9], China [10] and Brazil [11]. The blaKPC gene (that encodes the KPC enzyme) is mobilized on Tn4401, an active transposon belonging to the Tn3 family, with a size of 10 kb, that is delimited by two 39-bp inverted repeat sequences, and its replicative transposition involves a 5-bp target site duplication, without target site specificity [12]. Currently, eight different isoforms of Tn4401 have been reported (a-h), which differ by size-variable deletions upstream at the blaKPC gene, truncated tnpA or lacking tnpR [13, 14]. In P. aeruginosa, the blaKPC gene is rare (with variant 2 being the most frequent) and is mainly mobilized on both complete or truncated Tn4401b transposons, which have been reported in two different plasmid backbones (8 kb and 31.5 kb) belonging to the IncU and IncP-6 incompatibility groups, respectively [15, 16]. This blaKPC-2 variant has been identified in P. aeruginosa isolates belonging to sequence type (ST) 308, 1006, 1060 and 235 [15,16,17].

Pseudomonas aeruginosa ST235 is an infectious and virulent clone that possesses an extraordinary ability to acquire mobile genetic elements (MGE) and has been associated with the spread of multiple resistance genes, including genes that encode for carbapenemases [18,19,20]. Currently, this clone has been identified across the world, showing great disseminative success [18, 19, 21, 22]. The genomic plasticity of P. aeruginosa can be seen in ST235 strains, where the content of the resistance genes varies greatly due to the capture, permanence and rearrangement of resistance regions [18, 23, 24]. In Colombia, the ST235 clone has become prevalent since 2005 and has suffered an evolution in clinically important resistance genes by lateral gene transfer (LGT), e.g. some strains have acquired an atypical active class 1 integron with an unusual 3′ consensus sequence (tni module), which can provide a putative self-transposition ability [25], the acquisition of IS26 and the blaKPC-2 gene, which are most frequent in Enterobacteriaceae but rare in P. aeruginosa [2, 25]. In this study, after active clinical surveillance in the ICU at five distant cities, we found a broad population distribution of P. aeruginosa in Colombia; nonetheless, the ST235 clone was the only one harbouring blaKPC-2, and interestingly, those ST235 isolates harboured two copies of this gene that were located in two Tn4401b transposons in a new genomic island (GI) at the chromosome. As far as we know, this is the first report of a double chromosomal blaKPC-2 insertion in P. aeruginosa.

Methods

Strain collection, susceptibility profile and PCR detection of resistance genes

This prospective study was carried out based on the active surveillance of P. aeruginosa infections in adult patients who were attended in five ICUs in five distant cities in Colombia between October 2014 and March 2016. A total of 58 P. aeruginosa isolates were recovered. The susceptibility profiles of the isolates to meropenem, imipenem, ertapenem and doripenem, ceftazidime, gentamicin, amikacin, ciprofloxacin, trimethoprim/sulfamethoxazole, and Piperacillin/Tazobactam were determined by the automated systems VITEK® 2 or Microscan®. Additionally, the minimum inhibitory concentration (MIC) was established by the broth dilution method. The susceptibility profiles were established using the breakpoints defined by the Clinical and Laboratory Standards Institute according to the 2017 update [26]. The ATCC 9721 P. aeruginosa strain was used as control. In addition, the following resistance genes were assessed in all the isolates: blaTEM, blaSHV, blaCTX-M, blaFOX, blaACT, blaMIR blaACC, blaDHA, blaCMY, blaMOX, blaIMP, blaOXA-48, blaVIM, blaGES, blaKPC, blaNDM, blaOXA-23, blaOXA-24, blaOXA-43, blaOXA-51, blaOXA-58, aac(6′)-lb, aac(6′)-lb-cr, qnrA, B y S, mcbg, sul1, sul2 and sul3. The blaKPC-2 variant was determined by digestion with RsaI enzyme and confirmed by sequencing. Finally, the intl1, intl2 and intl3 genes were assessed.

Phylogenetic relationship among Pseudomonas aeruginosa isolates

The genetic relatedness between the isolates was determined by genome macro-restriction using the SpeI enzyme (Promega) and PFGE separation according to the protocol reported by Herschleb et al [27]. The dendrogram was obtained using the GelCompar II program (Applied Maths NV) with a tolerance position of 1.5% and a Dice coefficient of 1.0%. Two pulsotypes were considered different when they had < 80% similarity. The sequence type (ST) to some representative isolates of the most frequent pulsotypes was determined using the protocol reported by Curran et al [28].

Whole genome sequencing of a ST235 bla KPC-2-positive Pseudomonas aeruginosa isolate

The complete genome sequence of the ST235 blaKPC-2-positive 24pae112 isolate was obtained using the PacBio RS II platform (Pacific Biosciences, USA) and assembled through the previously reported procedure [29]. In brief, the total DNA was extracted, and a BluePippin (Sage Science) 20-kb size-selected library was prepared and then sequenced using one SMRT cell on the PacBio RS II platform. The sequencing process of the 24Pae112 isolate yielded 136,958 reads and a total of 1,221,367,395 bases (coverage of 172X). Sequencing reads were de novo assembled using the HGAP 3 protocol with default parameters. The assembly was visually inspected and manually verified using BWA-MEM (Burrows-Wheeler Aligner with maximal exact matches) [30] and Tablet v1.15.09.01 [31]. Misassembled terminal repeat overlap sequences were identified with Gepard (Genome Pair Rapid Dotter) [32] and trimmed manually. The whole genome was annotated using Prokka v1.11 [33], and the interesting regions were manually confirmed using BLASTn and BLASTp and edited in Artemis [34].

Phylogenetic analysis and comparative genomics of ST235 isolates

To perform the phylogenetic analysis and comparative genomics, the 43 publicly available P. aeruginosa ST235 genomes at the time were included (41 drafts and 2 complete genomes available at the NCBI genomes repository in January 2017), plus the genomes of the 24pae112 and 9Ps50 isolates (Additional file 1: Table S1). This 9Ps50 isolate was previously identified in 2006 by our group, causing an infection in institution C (included in this study), and did not harbour blaKPC-2 [25]. To build the phylogenetic tree based on the core-genome SNPs, partially assembled genomes were annotated using Prokka v1.11 [33], and an alignment of the concatenated core genes (genes present in all genomes with ≥90% of nucleotide identity) was created with Roary [35] using PRANK [36]. Poorly aligned positions and divergent regions were eliminated using Gblocks [37]. The phylogenetic tree was created using RAxML version 8.2.9 [38] running 1000 bootstrap replicates under the generalized time reversible model (GTRCAT). Finally, the consensus tree was plotted using Dendroscope [39]. Branch lengths are expressed in changes/nucleotide position (scale bar). In addition, an accessory-genome-based phylogenetic tree was built (presence/absence of genes) to assess the genetic diversity of the foreign DNA into the ST235 clone using the programs mentioned above.

A genomic comparative analysis was performed to determine the genetic variability in all P. aeruginosa ST235 strains using BLASTn [40] and visualized with BRIG. Genomic and pathogenicity islands (GI and PI, respectively), prophages, insertion sequence (IS) elements and transposons were identified using IslandViewer [41], PHASTER [42] and ISfinder [43]. The genetic platforms identified were visualized and illustrated by Artemis [34], ACT [44] and EasyFig [45].

Resistome identification

To determine the resistance gene content variability in the P. aeruginosa ST235 clone, the antibiotic resistance genes were identified in the genome of all the sequenced ST235 isolates using ARIBA (https://github.com/sanger-pathogens/ariba/wiki), with the ResFinder [46], CARD (The Comprehensive Antibiotic Resistance Database) [47] and ARG-ANNOT (Antibiotic Resistance Gene-ANNOTation) [48] databases. The results were visualized through Phandango software.

Results

In the surveillance study, we identified 58 P. aeruginosa that caused ICU infections among adult patients. With respect to the resistance behaviour, among all 58 isolates, 27 (46.6%) were resistant to at least one carbapenem antibiotic (imipenem, meropenem, ertapenem or doripenem). Of these 27 carbapenem-resistant P. aeruginosa, six (22.2%) and four (14.8%) isolates contained the blaVIM and blaKPC-2 genes, respectively. Multilocus sequence typing shows that the six blaVIM-positive isolates belonged to four different ST’s: ST111 (two isolates) previously reported in other Colombian cities [17]; ST244 (one isolate) reported in Europe [49, 50], Asia [21] and Africa [51]; ST1978 (two isolates), which was not reported previously; and ST235 (one isolate). Interestingly, all blaKPC-2-positive isolates only belonged to ST235. Analysis of the meropenem MIC showed a heterogeneous behaviour among the carbapenem-resistant isolates. Nine isolates reached an MIC value of 1024 μg/ml; among them were all the blaKPC-2-positive ST235 isolates and two blaVIM-positive ST111 isolates. Curiously, three isolates had an MIC of 1024 μg/ml but did not contain any of the assessed carbapenemases genes (Fig. 1).

Fig. 1
figure 1

Behaviour of the minimum inhibitory concentration (MIC) to meropenem of the 27 carbapenem-resistant isolates. Yellow and red points indicate the VIM-producing and KPC-producing isolates, respectively. ST: Sequence type

The PFGE analysis revealed 34 pulsotypes in the 58 P. aeruginosa isolates, showing a polyclonal behaviour. The six VIM-producing isolates were found in four unrelated PFGE pulsotypes, each of them with a different ST, suggesting that diverse clones have acquired this carbapenemase (Additional file 2: Figure S1). In contrast, in the four KPC-2-producing isolates (24Pae65, 24Pae112, 24Pae153 and 24Pae250), two very similar PFGE pulsotypes (a unique band of difference, > 95% of similarity) were found (Additional file 2: Figure S1).

Genomic comparative analysis of the ST235 isolates

In order to go deeper into the genetic platforms mobilizing the blaKPC-2 and other resistance genes, the complete genome sequence of the 24Pae112 isolate was established. After PacBio sequencing and assembly process, the 24Pae112 genome had a size of 7097,245 bp, in a unique contig (suggesting that this isolate does not contains plasmids), with a 66.0% G + C content, comprising 6634 predicted genes (coding density of 88.2%), which include 328 (4.9%) membrane transport proteins, 224 (3.4%) cell wall synthesis proteins, 202 (3.0%) virulence factors (pathogenicity and defence), 97 (1.5%) proteins of mobile genetic elements (prophages, insertion sequences and transposons), 63 (0.9%) tRNAs and 12 (0.2%) rRNAs. The antimicrobial resistome to the 24Pae112 chromosome evaluated using ARIBA included resistance genes to aminoglycosides (ant(2″)-Ia, ant(3″)-Ia and aph(3′)-IIb), beta-lactams (blaOXA-2 and blaKPC-2), sulphonamides (sul1), chloramphenicol (catB7) and fosfomycin (fosA). In addition, 8 prophages, 7 genomic islands (GI) and the different transposons within them were identified in the genome (Table 1). Recently, 40 and three ST235 genomes (complete or draft) were reported at the public NCBI genome database from the different countries, and we performed a phylogenetic analysis based on the core genome SNPs to determine the relationship between them and the Colombian ST235 isolates (Fig. 2). This analysis reveals that ST235 can be divided in two different genetic clades (Clade 1 and 2), showing the circulation of two genetic “sublineages” within the ST235 clone. No phylogenetic correlation with the origin country or identification date could be demonstrated because several isolates do not cluster according to their geographic origin or the date when they were identified, e.g. of the nine isolates identified in the USA, seven and two were grouped in clade 2 and 1, respectively; the six Colombian isolates were grouped in each of the clades. Meanwhile, the 24Pae112 and 9Ps50 isolates were grouped in different clades. These results suggest that two genetic “sublineages” are circulating within the ST235 clone. Interestingly, the accessory-genome-based tree of the ST235 clone contained more branches than the core-genome-based tree, showing a relatively wide genetic variability in the accessory genome of the ST235 isolates. The gene content on the accessory genome of the 24Pae112 isolate was closely related to the AZPAE14726 isolate that was identified in the USA in 2012, and very distant to the other Colombian ST235 isolates suggesting that the 24Pae112 isolate acquired a very different accessory DNA, perhaps by differential LGT events.

Table 1 Genomic islands and prophages present in the Pseudomonas aeruginosa 24Pae112 genome
Fig. 2
figure 2

Maximum likehood phylogenetic tree based on core (a) and accessory (b) genome of the 45 ST235 isolates. The isolates are named as follow: strain name, acronym of the origin country in bold (see Additional file 1: Table S1) and date of identification. The strains that start as AZPAE1 were simplified in the figure as A1. The branches into the tree were established when the sum of all internal nodes were ≤ 2, every branch is show in different colours. The bar scale (1.0 and 0.001) indicate the SNP number

Pseudomonas aeruginosa ST235 resistome

The P. aeruginosa ST235 clone is recognized for its ability to acquire and spread resistance genes worldwide. We determined the resistance genes repertoire in the 45 ST235 genomes included in the study. Multiple resistance genes to different families of antibiotics were found (Fig. 3). Of note, five different genes related to carbapenem resistance were found distributed in 8 countries: blaGES, blaKPC-2, blaVIM-11, blaIMP-1, and blaVIM-13 in 10 (22.2%), 4 (8.9%), 4 (8.9%), 3 (6.7%) and 1 (2.2%) isolates, respectively. Nonetheless, the gene blaKPC-2 (the unique variant detected in this clone) was exclusively found in the Colombian ST235 isolates, showing that the acquisition of this gene by the ST235 clone may be a particular phenomenon in our country. In addition, 22 gene variants related to aminoglycoside resistance were identified. The qnr, mcbg and aac(6′)-Ib-cr genes related to plasmid-mediated quinolone resistance were not found.

Fig. 3
figure 3

Resistome of Pseudomonas aeruginosa ST235 clone. Resistance genes identified in the strains belonging to ST235. The blue and red squares represent the presence of the resistance genes and the isolates harbouring the blaKPC gene, respectively. The isolates are named as follow: strain name, acronym of the origin country in bold (see Additional file 1: Table S1) and date of identification. The orange circle indicates the position of the 24Pae112 isolate

Diverse genetic platforms transport resistance genes in 24Pae112 isolate

To gain a better understanding of the complex process of the acquisition of resistance genes by the 24Pae112 isolate, the genetic platforms transporting these genes were identified, and the genome localization was compared with the remaining ST235 genomes and visualized with BRIG (Figs. 4 and 5). Five transposons were identified to harbour resistance genes: the Tn6162-like transposon with a size of 34,220 bp, which was located in P. aeruginosa genomic island 1 (GI-1). This transposon harbours a classical class 1 integron with the gene cassette array aadB-gcuD (Fig. 5a). Sample 24Pae112 also harbours two distant copies of the Tn402-like transposon, with the first one in the prophage F10 and the second one in the PAGI-5. Curiously, both transposons contained a class 1 integron structure, but with different variable regions, harbouring only one gene cassette, aadA1 and blaOXA-2, respectively (Fig. 5b and c). In addition, two blaKPC-2-harbouring Tn4401b transposons were also identified (see below). It was not possible to determine the genetic platforms mobilizing the aph(3′)-IIb, catB7 and fosA genes.

Fig. 4
figure 4

Genomic comparison of Pseudomonas aeruginosa ST235 isolates. Matched alignment between 24Pae112 genome and the 40 and four remain ST235 genomes visualized using BRIG, the order of the strains in the rings can be seen in Additional file 1: Table S1. The colours of every ring indicate the continent of the origin and the colour intensity represents the identity percentage. The outer ring shows the location of resistance genes and their genetic platforms of mobilization (when was identified), genomic islands and prophages in 24Pae112

Fig. 5
figure 5

Organization of the different genetic platform harbouring resistance genes within the 24Pae112 isolate. a Tn6162-like transposon inserted into the genomic island 1. b First copy of the Tn402-like transposon inserted into the prophage F10. c Second copy of the Tn402-like transposon within the PAGI-5. d Two copies of the Tn4401 transposon into the new genomic island PAGI-17. IR: Inverted repeats

Two bla KPC-2-harbouring Tn4401b are inserted into the 24Pae112 chromosome

Surprisingly, the genome sequence analysis of 24Pae112 reveals the presence of two complete and intact copies of blaKPC-2-harbouring Tn4401 (isoform b) transposons inserted into the chromosome within a new genomic island (Figs. 5d and 6). Comparative genomic analysis allowed us to establish that this same genomic island had already been described by Rau et al. in the P. aeruginosa DK2 clone and named generically as GI4 [52]; however, the sequence and structure of this island is different with respect to the 16 PAGI currently reported (including PAGI-4, data not shown) [53,54,55,56,57]. Consequently, here, we renamed this genomic island as Pseudomonas aeruginosa genomic island 17 (PAGI-17), according to the widely accepted nomenclature system used for genomic islands in this microorganism [53,54,55,56,57]. This double Tn4401b chromosomal insertion was experimentally confirmed by PCR using specific primers (data not shown). The two Tn4401b copies were identical, with a size of 10,006 bp, in an inverted orientation, for a separate DNA fragment of 5877 bp, which harbours the parB, yhbS, chrA, padR and arsC genes. The first Tn4401b copy was inserted within the acr3 gene, which encodes the Acr3 protein related to arsenic resistance [58], with a 5-bp target site duplication (CTGCT). The second Tn4401b copy was inserted within the acrB gene, which encodes the transporter protein AcrB, related to acriflavine, antibiotics, disinfectants, dyes, and detergent resistance; AcrB gene is located upstream of the acrA gene which encodes the RND multidrug efflux membrane permease protein AcrA. The second insertion generated a 5-bp distinct target site duplication (GTCCA) (Fig. 7). This confirms the lack of target site specificity of this transposon, as has previously been reported [12].

Fig. 6
figure 6

Comparison of the different genetics platform mobilizing blaKPC-2 gene in Pseudomonas aeruginosa. The pPA-2 (KC609322.1) and pCOL-1 (KC609323.1) plasmids were identified in Colombia; the p10265-KPC plasmid (KU578314.1) was identified in China

Fig. 7
figure 7

Insertion sites analysis of the two blaKPC-2-harbouring Tn4401b transposons and PAGI-17. The PAGI-17 was inserted at tRNAGly gene as it is showed in the PAGI-17-negative PAO1 reference strain (NC_002516.2). The two couples of Tn4401b direct repeats and the three direct repeats of the PAGI-17 are showed. The GenBank accession number for DK2, LESlike1, LESlike4, LESlike5 and LESlike7 genomes are NC_018080.1, NZ_CP006984.1, NZ_CP006985.1, NZ_CP006980.1 and NZ_CP006981.1; respectively

Although the information about PAGI-17 is very limited, to unravel the genetic mechanisms that possibly allow the Tn4401b chromosomal insertion, a more detailed analysis of this genomic island was carried out. Comparing the chromosome sequence of the reference strain P. aeruginosa PAO1 (devoid of PAGI-17) with the 24Pae112 genome, we determined that in 24Pae112, the PAGI-17 was inserted at the tRNAGly gene, in the positions 2,724,254 to 2,841,225 (to the 24Pae112 genome). The PAGI-17 insertion produced the duplication of a 16-bp sequence (CTCGTTTCCCGCTCCA) located at the 3′-end of the tRNAGly gene (Fig. 7). The duplication did not alter the sequence of this gene. The total size of PAGI-17 in 24Pae112 was 116,971 bp, with a G + C content of 64.1% and 113 open reading frames. Comparing PAGI-17 from 24Pae112 with PAGI-17 from P. aeruginosa DK2, we identified an additional 16-bp identical sequence (CTCGTTTCCCGCTCCA, similar to that generated by the PAGI-17 insertion in 24Pae112) located downstream of the second duplicated sequence. This new duplication delimits a fragment of a 2264-bp, which is unique in the 24Pae112 PAGI-17 and harbours the fimA gene that encodes the type 1 fimbria component protein, an important virulence factor that mediates mannose-sensitive adherence to the eukaryotic cell.

Analysis of the chromosomal environment surrounding the Tn4401b insertion sites revealed a 24,134-bp specific DNA fragment, which is rarely found in P. aeruginosa since it was only found in PAGI-17 in five unrelated P. aeruginosa strains: DK2-CF333(ST386), SCV20265 (ST299), PAER4 (ST260), 8380(ST2619) and BAMCPAO7 (ST313). None of the ST235 isolates contained this region (Fig. 7). While the complete or truncated PAGI-17 is present in several P. aeruginosa strains, only six strains (including the 24Pae112) possess this 24,134-bp-size DNA fragment. Interestingly, a genomic comparative analysis carried out by Rau et al. on 45 isolates belonging to the DK2 lineage and recovered from patients suffering from cystic fibrosis over the course of 35 years (from 1973 to 2008) found that this 24,134-bp DNA fragment was acquired early in the DK2 lineage. Specifically, the DK2-CF114 isolate (the first DK2-lineage isolate recovered at 1973) did not contain exactly this fragment, whereas the remaining 44 isolates incorporated it into their chromosomes [52]. This result suggests that this DNA portion was possibly acquired by LGT in P. aeruginosa a long time ago in a few genetic lineages, including the ST235 clone. The acquisition of this DNA may favour the incorporation and permanence of the two Tn4401b-blaKPC-2 transposons.

Discussion

The problem of resistance to carbapenems has become highly relevant for the health care authorities, as this class of antibiotics has frequently been used as the last-line-of-defence agents to treat serious bacterial infections. Additionally, carbapenem resistance is a much bigger problem in non-fermentering bacteria, such as P. aeruginosa and Acinetobacter baumanni, because these bacteria possess frequent resistance (intrinsic or acquired) to other families of antibiotics, and the lack of available treatment options against some MDR strains is worrying. In our study, we found that almost half (46.6%) of the P. aeruginosa isolates were already resistant to carbapenems, which is alarming given that in Colombia the number of MDR P. aeruginosa infections in hospitalized patients has risen [59]. High percentages of carbapenem resistance have been reported in other countries, such as Brazil (43.9%), Spain (88.0%) and Iran (96.0%) [60,61,62]. Additionally, among all the carbapenem resistant isolates, almost 40% are produced by carbapenemase production (VIM and KPC) in different clones (ST111, ST244, ST1978 and ST235). We identified four blaKPC-2-positive and multidrug resistant isolates belonging to ST235 (only susceptible to colistin), but these do not contain the blaKPC-2–harbouring plasmids that were previously reported [15, 16]. For these reasons, we performed complete genome sequencing of the blaKPC-2-harbouring P. aeruginosa isolates to establish the different genetic platforms moving the resistance genes. This analysis led us to find a new feature in the dynamic blaKPC-2 gene mobilization through Tn4401 in this bacterium.

The blaKPC-2 gene in P. aeruginosa was first reported in Colombia in 2007 [2] but has already been found in several countries, including the United States of America and China [3, 10, 11]. Previous studies in P. aeruginosa revealed that the blaKPC-2 gene is mobilized into only one copy of Tn4401 (complete or truncated forms) within the two plasmid backbones. However, comparison of the genetic structures surrounding Tn4401b in the 24Pae112 isolate indicate to us that its incorporation was not produced by the chromosome integration of any two previously identified Inc-P or Inc-U plasmids in Colombia (Fig. 6) [15]. Conversely, the Tn4401b acquisition occurred by transposition at two different, but nearby sites. As far as we know, this is the first report of a double chromosomal blaKPC-2–harbouring Tn4401 insertion in P. aeruginosa. Even though this duplication is a rare genetic event, it was previously reported within a 70.6-kb plasmid from a Klebsiella pneumoniae S9 strain identified in New York [63]. Although the Tn4401 duplication could provide higher production of the KPC enzyme, it is necessary to determine whether it provides any advantage in the carbapenem hydrolysis and consequently an increase in the carbapenem MICs.

The origin of Tn4401b in the 24Pae112 isolate is unknown, although it may be related to the process of intercellular mobility (the MGE may be excised from one location and to undergo replicative transposition before to be reintegrated elsewhere in the genome) [64] from a plasmid. The 24Pae112 isolate probably acquired a blaKPC-2–harbouring plasmid from a Klebsiella pneumoniae isolate that favoured Tn4401b transposition to the chromosome, given that the KPC-producing Klebsiella pneumoniae clone is already endemic in Colombia. This Tn4401b mobilization from the plasmid to the chromosome has already been reported in Klebsiella pneumoniae [63].

Our comparative genomic analysis shows a relatively high variability in the accessory genome of ST235 isolates, suggesting great genome plasticity in this clone. This genomic flexibility could contribute to the fact that the ST235 clone can survive in the most hostile environments. The variability in the accessory genome is produced for the differential acquisition of diverse MGE, such as prophages, genomic islands and transposons (in some cases moving resistance genes), and for changes within them produced by genetic recombination events, duplications and deletions mediated by transposases and recombinases [64]. In the case of the 24Pae112 isolate, this was probably necessary for several genetic events to occur before the double Tn4401 insertion was done. First, the acquisition of PAGI-17 by different genetics lineages, including an ST235 common ancestor because it is in almost all the ST235 isolates. Second, the insertion of a 24-kb fragment within the 5′ region of PAGI-17, this is not a genetic phenomenon exclusively associated with the ST235 clone, given that it was only found in 24Pae112. However, it is also a rare event in P. aeruginosa because it was only found in five isolates belonging to different genetic lineages (ST386, ST299, ST260, ST2619, and ST313). Third, the acquisition and maintenance of a MGE (plasmid) harbouring a Tn4401b transposon, and fourth, the transposition of the Tn4401 from MGE to the chromosome within a site that did not affect the essential genes for the bacteria, such as PAGI-17.

From the clinical point of view, the great genomic plasticity of P. aeruginosa and its extraordinary ability to incorporate several copies of resistance genes generates an elevated production of B-lactam antibiotic-degrading enzymes; this limits the carbapenem use and pose the urgent need to find different treatment strategies, such as the use of combined antibiotics and the development of new therapeutic technologies. Accordingly, knowing in detail the mechanisms of bacterial resistance is an important factor to contain and prevent this great problem. For other hand, the 24Pae112 isolate was recovered from a remitted patient, who had a healthcare-associated bloodstream infection, suggesting that this P. aeruginosa clone could be already circulating in other regions of the country. Understanding the dynamics of the transfer of elements that confer resistance is one of the important factors to design the best practices to prevent the spread of multidrug resistant bacteria and healthcare-associated infections.

Conclusion

Our findings together with those previously reported in the ST235 clone show the extraordinary genomic plasticity of Pseudomonas aeruginosa, quality that allows it to acquire, keep and modify foreign DNA into its genome. This DNA acquisition contributes to more genetic armament to face new hostile environmental conditions, such as high concentrations of antibiotics. Our findings also indicate the essential role that genomic islands play in the constant and successful evolution of P. aeruginosa. This has been demonstrated here with the PAGI-17, but it has also been demonstrated for PAGI-1 and -2 [22] because they serve as both catchers and natural “warehouses” of foreign DNA storage, where this bacterium can keep virulence and resistance genes with few secondary effects on its genome stability.

Abbreviations

CLSI:

Clinical and Laboratory Standards Institute

GI:

Genomic islands

ICU:

Intensive care unit

IS:

Insertion sequence

KPC:

Klebsiella pneumoniae carbapenemase

LGT:

Lateral gene transfer

MDR:

Multidrug-resistant

MGE:

Mobile genetics elements

MIC:

Minimum inhibitory concentration

MLST:

Multilocus sequence typing

ORF:

Open reading frame

PAGI:

Pseudomonas aeruginosa genomic islands

PCR:

Polymerase chain reaction

PFGE:

Pulsed-field gel electrophoresis

SNP:

Single nucleotide polymorphism

ST:

Sequence type

VIM:

Verona Integron-encoded Metallo-β-lactamase

WGS:

Whole-genome sequencing

References

  1. Poole K. Pseudomonas aeruginosa: resistance to the max. Front Microbiol. 2011;2:65.

    Article  CAS  Google Scholar 

  2. Villegas MV, Lolans K, Correa A, Kattan JN, Lopez JA, Quinn JP, et al. First identification of Pseudomonas aeruginosa isolates producing a KPC-type carbapenem-hydrolyzing beta-lactamase. Antimicrob Agents Chemother. 2007;51:1553–5.

    Article  CAS  Google Scholar 

  3. Poirel L, Nordmann P, Lagrutta E, Cleary T, Munoz-Price LS. Emergence of KPC-producing Pseudomonas aeruginosa in the United States. Antimicrob Agents Chemother. 2010;54:3072.

    Article  CAS  Google Scholar 

  4. Paul D, Dhar D, Maurya AP, Mishra S, Sharma GD, Chakravarty A, et al. Occurrence of co-existing bla VIM-2 and bla NDM-1 in clinical isolates of Pseudomonas aeruginosa from India. Ann Clin Microbiol Antimicrob. 2016;15:31.

    Article  Google Scholar 

  5. Rahman M, Prasad KN, Pathak A, Pati BK, Singh A, Ovejero CM, et al. RmtC and RmtF 16S rRNA methyltransferase in NDM-1-producing Pseudomonas aeruginosa. Emerg Infect Dis. 2015;21:2059–62.

    Article  CAS  Google Scholar 

  6. Khan AU, Maryam L, Zarrilli R. Structure, genetics and worldwide spread of New Delhi Metallo-beta-lactamase (NDM): a threat to public health. BMC Microbiol. 2017;17:101.

    Article  Google Scholar 

  7. Nordmann P, Dortet L, Poirel L. Carbapenem resistance in Enterobacteriaceae: here is the storm! Trends Mol Med. 2012;18:263–72.

    Article  CAS  Google Scholar 

  8. Akpaka PE, Swanston WH, Ihemere HN, Correa A, Torres JA, Tafur JD, et al. Emergence of KPC-producing Pseudomonas aeruginosa in Trinidad and Tobago. J Clin Microbiol. 2009;47:2670–1.

    Article  CAS  Google Scholar 

  9. Robledo IE, Aquino EE, Vazquez GJ. Detection of the KPC gene in Escherichia coli, Klebsiella pneumoniae, Pseudomonas aeruginosa, and Acinetobacter baumannii during a PCR-based nosocomial surveillance study in Puerto Rico. Antimicrob Agents Chemother. 2011;55:2968–70.

    Article  CAS  Google Scholar 

  10. Ge C, Wei Z, Jiang Y, Shen P, Yu Y, Li L. Identification of KPC-2-producing Pseudomonas aeruginosa isolates in China. J Antimicrob Chemother. 2011;66:1184–6.

    Article  CAS  Google Scholar 

  11. Jacome PR, Alves LR, Cabral AB, Lopes AC, Maciel MA. First report of KPC-producing Pseudomonas aeruginosa in Brazil. Antimicrob Agents Chemother. 2012;56:4990.

    Article  CAS  Google Scholar 

  12. Cuzon G, Naas T, Nordmann P. Functional characterization of Tn4401, a Tn3-based transposon involved in blaKPC gene mobilization. Antimicrob Agents Chemother. 2011;55:5370–3.

    Article  CAS  Google Scholar 

  13. Naas T, Cuzon G, Truong HV, Nordmann P. Role of ISKpn7 and deletions in blaKPC gene expression. Antimicrob Agents Chemother. 2012;56:4753–9.

    Article  CAS  Google Scholar 

  14. Cheruvanky A, Stoesser N, Sheppard AE, Crook DW, Hoffman PS, Weddle E, et al. Enhanced Klebsiella pneumoniae Carbapenemase expression from a novel Tn4401 deletion. Antimicrob Agents Chemother. 2017;61.

  15. Naas T, Bonnin RA, Cuzon G, Villegas MV, Nordmann P. Complete sequence of two KPC-harbouring plasmids from Pseudomonas aeruginosa. J Antimicrob Chemother. 2013;68:1757–62.

    Article  CAS  Google Scholar 

  16. Dai X, Zhou D, Xiong W, Feng J, Luo W, Luo G, et al. The IncP-6 plasmid p10265-KPC from Pseudomonas aeruginosa carries a novel DeltaISEc33-associated bla KPC-2 gene cluster. Front Microbiol. 2016;7:310.

    PubMed  PubMed Central  Google Scholar 

  17. Vanegas JM, Cienfuegos AV, Ocampo AM, Lopez L, del Corral H, Roncancio G, et al. Similar frequencies of Pseudomonas aeruginosa isolates producing KPC and VIM carbapenemases in diverse genetic clones at tertiary-care hospitals in Medellin, Colombia. J Clin Microbiol. 2014;52:3978–86.

    Article  Google Scholar 

  18. Kos VN, Deraspe M, McLaughlin RE, Whiteaker JD, Roy PH, Alm RA, et al. The resistome of Pseudomonas aeruginosa in relationship to phenotypic susceptibility. Antimicrob Agents Chemother. 2015;59:427–36.

    Article  Google Scholar 

  19. Cho HH, Kwon KC, Sung JY, Koo SH. Prevalence and genetic analysis of multidrug-resistant Pseudomonas aeruginosa ST235 isolated from a hospital in Korea, 2008-2012. Ann Clin Lab Sci. 2013;43:414–9.

    CAS  PubMed  Google Scholar 

  20. Kim MJ, Bae IK, Jeong SH, Kim SH, Song JH, Choi JY, et al. Dissemination of metallo-beta-lactamase-producing Pseudomonas aeruginosa of sequence type 235 in Asian countries. J Antimicrob Chemother. 2013;68:2820–4.

    Article  CAS  Google Scholar 

  21. Hong JS, Yoon EJ, Lee H, Jeong SH, Lee K. Clonal dissemination of Pseudomonas aeruginosa sequence type 235 isolates carrying blaIMP-6 and emergence of blaGES-24 and blaIMP-10 on novel Genomic Islands PAGI-15 and -16 in South Korea. Antimicrob Agents Chemother. 2016;60:7216–23.

    Article  CAS  Google Scholar 

  22. Roy Chowdhury P, Scott M, Worden P, Huntington P, Hudson B, Karagiannis T, et al. Genomic islands 1 and 2 play key roles in the evolution of extensively drug-resistant ST235 isolates of Pseudomonas aeruginosa. Open Biol. 2016;6.

  23. Liakopoulos A, Mavroidi A, Katsifas EA, Theodosiou A, Karagouni AD, Miriagou V, et al. Carbapenemase-producing Pseudomonas aeruginosa from Central Greece: molecular epidemiology and genetic analysis of class I integrons. BMC Infect Dis. 2013;13:505.

    Article  Google Scholar 

  24. Edelstein MV, Skleenova EN, Shevchenko OV, D'Souza JW, Tapalski DV, Azizov IS, et al. Spread of extensively resistant VIM-2-positive ST235 Pseudomonas aeruginosa in Belarus, Kazakhstan, and Russia: a longitudinal epidemiological and clinical study. Lancet Infect Dis. 2013;13:867–76.

    Article  Google Scholar 

  25. Martinez E, Perez JE, Buelvas F, Tovar C, Vanegas N, Stokes HW. Establishment and multi drug resistance evolution of ST235 Pseudomonas aeruginosa strains in the intensive care unit of a Colombian hospital. Res Microbiol. 2014;165:852–6.

    Article  CAS  Google Scholar 

  26. Institute. CaLS: M100-S25 performance standards for antimicrobial susceptibility testing; twenty-fifth informational supplement. In: Wayne, PA 19087 USA: 2015; 2015. p. 237.

    Google Scholar 

  27. Herschleb J, Ananiev G, Schwartz DC. Pulsed-field gel electrophoresis. Nat Protocols. 2007;2:677–84.

    Article  CAS  Google Scholar 

  28. Curran B, Jonas D, Grundmann H, Pitt T, Dowson CG. Development of a multilocus sequence typing scheme for the opportunistic pathogen Pseudomonas aeruginosa. J Clin Microbiol. 2004;42:5644–9.

    Article  CAS  Google Scholar 

  29. Marquez-Ortiz RA, Haggerty L, Olarte N, Duarte C, Garza-Ramos U, Silva-Sanchez J, et al. Genomic epidemiology of NDM-1-encoding plasmids in Latin American clinical isolates reveals insights into the evolution of multidrug resistance. Genome Biol Evol. 2017;9:1725–41.

    Article  CAS  Google Scholar 

  30. Li H, Durbin R. Fast and accurate long-read alignment with burrows-wheeler transform. Bioinformatics. 2010;26:589–95.

    Article  Google Scholar 

  31. Milne I, Stephen G, Bayer M, Cock PJ, Pritchard L, Cardle L, et al. Using tablet for visual exploration of second-generation sequencing data. Brief Bioinform. 2013;14:193–202.

    Article  CAS  Google Scholar 

  32. Krumsiek J, Arnold R, Rattei T. Gepard: a rapid and sensitive tool for creating dotplots on genome scale. Bioinformatics. 2007;23:1026–8.

    Article  CAS  Google Scholar 

  33. Seemann T. Prokka: rapid prokaryotic genome annotation. Bioinformatics. 2014;30(14):2068–9.

    Article  CAS  Google Scholar 

  34. Rutherford K, Parkhill J, Crook J, Horsnell T, Rice P, Rajandream MA, et al. Artemis: sequence visualization and annotation. Bioinformatics. 2000;16:944–5.

    Article  CAS  Google Scholar 

  35. Page AJ, Cummins CA, Hunt M, Wong VK, Reuter S, Holden MT, et al. Roary: rapid large-scale prokaryote pan genome analysis. Bioinformatics. 2015;31:3691–3.

    Article  CAS  Google Scholar 

  36. Loytynoja A. Phylogeny-aware alignment with PRANK. Methods Mol Biol. 2014;1079:155–70.

    Article  Google Scholar 

  37. Talavera G, Castresana J. Improvement of phylogenies after removing divergent and ambiguously aligned blocks from protein sequence alignments. Syst Biol. 2007;56:564–77.

    Article  CAS  Google Scholar 

  38. Stamatakis A. RAxML version 8: a tool for phylogenetic analysis and post-analysis of large phylogenies. Bioinformatics. 2014;30:1312–3.

    Article  CAS  Google Scholar 

  39. Huson DH, Scornavacca C. Dendroscope 3: an interactive tool for rooted phylogenetic trees and networks. Syst Biol. 2012;61:1061–7.

    Article  Google Scholar 

  40. Altschul SF, Gish W, Miller W, Myers EW, Lipman DJ. Basic local alignment search tool. J Mol Biol. 1990;215:403–10.

    Article  CAS  Google Scholar 

  41. Langille MG, Brinkman FS. IslandViewer: an integrated interface for computational identification and visualization of genomic islands. Bioinformatics. 2009;25:664–5.

    Article  CAS  Google Scholar 

  42. Arndt D, Grant JR, Marcu A, Sajed T, Pon A, Liang Y, et al. PHASTER: a better, faster version of the PHAST phage search tool. Nucleic Acids Res. 2016;44:W16–21.

    Article  CAS  Google Scholar 

  43. Siguier P, Varani A, Perochon J, Chandler M. Exploring bacterial insertion sequences with ISfinder: objectives, uses, and future developments. Methods Mol Biol. 2012;859:91–103.

    Article  CAS  Google Scholar 

  44. Carver TJ, Rutherford KM, Berriman M, Rajandream MA, Barrell BG, Parkhill J. ACT: the Artemis comparison tool. Bioinformatics. 2005;21:3422–3.

    Article  CAS  Google Scholar 

  45. Sullivan MJ, Petty NK, Beatson SA. Easyfig: a genome comparison visualizer. Bioinformatics. 2011;27:1009–10.

    Article  CAS  Google Scholar 

  46. Zankari E, Hasman H, Cosentino S, Vestergaard M, Rasmussen S, Lund O, et al. Identification of acquired antimicrobial resistance genes. J Antimicrob Chemother. 2012;67:2640–4.

    Article  CAS  Google Scholar 

  47. McArthur AG, Waglechner N, Nizam F, Yan A, Azad MA, Baylay AJ, et al. The comprehensive antibiotic resistance database. Antimicrob Agents Chemother. 2013;57:3348–57.

    Article  CAS  Google Scholar 

  48. Gupta SK, Padmanabhan BR, Diene SM, Lopez-Rojas R, Kempf M, Landraud L, et al. ARG-ANNOT, a new bioinformatic tool to discover antibiotic resistance genes in bacterial genomes. Antimicrob Agents Chemother. 2014;58:212–20.

    Article  Google Scholar 

  49. Ruiz-Roldan L, Belles A, Bueno J et al. Pseudomonas aeruginosa isolates from Spanish children: occurrence in faecal samples, antimicrobial resistance, virulence, and molecular typing. Biomed Res Int 2018;2018:8060178.

  50. Rojo-Bezares B, Cavalie L, Dubois D, Oswald E, Torres C, Saenz Y. Characterization of carbapenem resistance mechanisms and integrons in Pseudomonas aeruginosa strains from blood samples in a French hospital. J Med Microbiol. 2016;65:311–9.

    Article  CAS  Google Scholar 

  51. Moyo S, Haldorsen B, Aboud S, Blomberg B, Maselle SY, Sundsfjord A, et al. Identification of VIM-2-producing Pseudomonas aeruginosa from Tanzania is associated with sequence types 244 and 640 and the location of blaVIM-2 in a TniC integron. Antimicrob Agents Chemother. 2015;59:682–5.

    Article  Google Scholar 

  52. Rau MH, Marvig RL, Ehrlich GD, Molin S, Jelsbak L. Deletion and acquisition of genomic content during early stage adaptation of Pseudomonas aeruginosa to a human host environment. Environ Microbiol. 2012;14:2200–11.

    Article  CAS  Google Scholar 

  53. Liang X, Pham XQ, Olson MV, Lory S. Identification of a genomic island present in the majority of pathogenic isolates of Pseudomonas aeruginosa. J Bacteriol. 2001;183:843–53.

    Article  CAS  Google Scholar 

  54. Larbig KD, Christmann A, Johann A, Klockgether J, Hartsch T, Merkl R, et al. Gene islands integrated into tRNA(Gly) genes confer genome diversity on a Pseudomonas aeruginosa clone. J Bacteriol. 2002;184:6665–80.

    Article  CAS  Google Scholar 

  55. Klockgether J, Reva O, Larbig K, Tummler B. Sequence analysis of the mobile genome island pKLC102 of Pseudomonas aeruginosa C. J Bacteriol. 2004;186:518–34.

    Article  CAS  Google Scholar 

  56. Battle SE, Meyer F, Rello J, Kung VL, Hauser AR. Hybrid pathogenicity island PAGI-5 contributes to the highly virulent phenotype of a Pseudomonas aeruginosa isolate in mammals. J Bacteriol. 2008;190:7130–40.

    Article  CAS  Google Scholar 

  57. Battle SE, Rello J, Hauser AR. Genomic islands of Pseudomonas aeruginosa. FEMS Microbiol Lett. 2009;290:70–8.

    Article  CAS  Google Scholar 

  58. Aaltonen EK, Silow M. Transmembrane topology of the Acr3 family arsenite transporter from Bacillus subtilis. Biochim Biophys Acta. 2008;1778:963–73.

    Article  CAS  Google Scholar 

  59. Villa LM, Cortes JA, Leal AL, Meneses A, Melendez MP, Grebo. Resistance to antibiotics in Pseudomonas aeruginosa in Colombian hospitals. Rev Chil Infectol. 2013;30:605–10.

    Article  Google Scholar 

  60. Rossi Goncalves I, Dantas RC, Ferreira ML, Batistao DW, Gontijo-Filho PP, Ribas RM. Carbapenem-resistant Pseudomonas aeruginosa: association with virulence genes and biofilm formation. Braz J Microbiol. 2016.

  61. Estepa V, Rojo-Bezares B, Azcona-Gutierrez JM, Olarte I, Torres C, Saenz Y. Characterisation of carbapenem-resistance mechanisms in clinical Pseudomonas aeruginosa isolates recovered in a Spanish hospital. Enferm Infecc Microbiol Clin. 2016.

  62. Radan M, Moniri R, Khorshidi A, Gilasi H, Norouzi Z, Beigi F, et al. Emerging Carbapenem-resistant Pseudomonas aeruginosa isolates carrying blaIMP among burn patients in Isfahan, Iran. Arch Trauma Res. 2016;5:e33664.

    Article  Google Scholar 

  63. Gootz TD, Lescoe MK, Dib-Hajj F, Dougherty BA, He W, Della-Latta P, et al. Genetic organization of transposase regions surrounding blaKPC carbapenemase genes on plasmids from Klebsiella strains isolated in a New York City hospital. Antimicrob Agents Chemother. 2009;53:1998–2004.

    Article  CAS  Google Scholar 

  64. Frost LS, Leplae R, Summers AO, Toussaint A. Mobile genetic elements: the agents of open source evolution. Nat Rev Microbiol. 2005;3:722–32.

    Article  CAS  Google Scholar 

Download references

Acknowledgments

We gratefully acknowledge to the participant institutions, Hospital Universitario San Jorge E.S.E., Hospital El Tunal E.S.E., Hospital Universitario Departamental de Nariño, Hospital San Jerónimo and Clínica Estrios for allowing us to conduct this study. We gratefully acknowledge to Vice Chancellery for Research of El Bosque University (especially to Dr. Miguel Otero for your invaluable support).

Availability of the data and materials

The annotated, complete 24pae112 genome and the trimmed and filtered MiSeq sequencing reads for 9Ps50 genome sequenced in this study are available in the DDBJ/EMBL/GenBank public databases under the BioProject number PRJNA416504, BioSample number SAMN07975642, and GenBank accession number NZ_CP029605.

Funding

This study was funded by the Departamento Administrativo de Ciencia, Tecnología e Innovación, Colciencias (grant number 1308–657-41107 607–2014 and partially by grant number 1308–569-33463) and Vice Chancellery for Research of El Bosque University (grant number PCI63–2014), University of Cartagena and University of Sinu. DA was supported by Colciencias through the grant number 1308–657-41107.

The funders had no role in study design, data collection and interpretation, or the decision to submit the work for publication.

Author information

Authors and Affiliations

Authors

Contributions

JEP designed research; JIMO, NMOE, NR, CT, HFS, YMGG and CELU collected the isolates; performed acquisition, analysis, and interpretation of data in each institution, DA, BCC, RAMO, and ZLCR performed research; RAMO, DA and JC performed the bioinformatics analysis; DA, RAMO, NVG and JEP interpreted the data and DA, RAMO and JEP wrote the paper. All authors read and approved the final manuscript.

Corresponding author

Correspondence to Javier Escobar-Pérez.

Ethics declarations

Ethics approval and consent to participate

The study and protocol were approved by the Institutional Ethics Committees of all participant institutions (El Bosque University, University of Cartagena and University of Sinu). The protocol was presented to all participants and written and signed informed consent was obtained. In addition, information and records of all participants included in the study were processed anonymously.

Consent for publication

Not applicable.

Competing interests

The authors declare that they have no competing interests.

Publisher’s Note

Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Additional files

Additional file 1:

Table S1. Characteristics of the P. aeruginosa ST235 genomes included in the phylogenetic tree and the genomic comparative analysis. (XLSX 16 kb)

Additional file 2:

Figure S1. Genetic relationship by PFGE of the carbapenem-resistant P. aeruginosa isolates identified in the study. The black line corresponds to a genetic relationship of 80% of similarity. Numbers in the colour boxes indicate to the main PFGE pulsotype found. ST: Sequence type. The abbreviations in the susceptibility profile: TZP: piperacillin/tazobactam, CAZ: ceftazidime, MEM: meropenem, DOR: doripenem, IMP: imipenem, STX: trimethoprim/sulfamethoxazole, CIP: ciprofloxacin, AK: amikacin, GEN: gentamicin y COL: colistin. The isolate marked with asterisk (*) was selected to whole genome sequencing. (TIF 445 kb)

Rights and permissions

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

Reprints and permissions

About this article

Check for updates. Verify currency and authenticity via CrossMark

Cite this article

Abril, D., Marquez-Ortiz, R.A., Castro-Cardozo, B. et al. Genome plasticity favours double chromosomal Tn4401b-blaKPC-2 transposon insertion in the Pseudomonas aeruginosa ST235 clone. BMC Microbiol 19, 45 (2019). https://doi.org/10.1186/s12866-019-1418-6

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: https://doi.org/10.1186/s12866-019-1418-6

Keywords