- Open Access
VIM-encoding IncpSTY plasmids and chromosome-borne integrative and mobilizable elements (IMEs) and integrative and conjugative elements (ICEs) in Pseudomonas
Annals of Clinical Microbiology and Antimicrobials volume 21, Article number: 10 (2022)
The carbapenem-resistance genes blaVIM are widely disseminated in Pseudomonas, and frequently harbored within class 1 integrons that reside within various mobile genetic elements (MGEs). However, there are few reports on detailed genetic dissection of blaVIM-carrying MGEs in Pseudomonas.
This study presented the complete sequences of five blaVIM-2/-4-carrying MGEs, including two plasmids, two chromosomal integrative and mobilizable elements (IMEs), and one chromosomal integrative and conjugative element (ICE) from five different Pseudomonas isolates.
The two plasmids were assigned to a novel incompatibility (Inc) group IncpSTY, which included only seven available plasmids with determined complete sequences and could be further divided into three subgroups IncpSTY-1/2/3. A detailed sequence comparison was then applied to a collection of 15 MGEs belonging to four different groups: three representative IncpSTY plasmids, two Tn6916-related IMEs, two Tn6918-related IMEs, and eight Tn6417-related ICEs and ten of these 15 MGEs were first time identified. At least 22 genes involving resistance to seven different categories of antibiotics and heavy metals were identified within these 15 MGEs, and most of these resistance genes were located within the accessory modules integrated as exogenous DNA regions into these MGEs. Especially, eleven of these 15 MGEs carried the blaVIM genes, which were located within 11 different concise class 1 integrons.
These blaVIM-carrying integrons were further integrated into the above plasmids, IMEs/ICEs with intercellular mobility. These MGEs could transfer between Pseudomonas isolates, which resulted in the accumulation and spread of blaVIM among Pseudomonas and thus was helpful for the bacteria to survival from the stress of antibiotics. Data presented here provided a deeper insight into the genetic diversification and evolution of VIM-encoding MGEs in Pseudomonas.
VIM hydrolyzes nearly all β-lactam except for aztreonam and currently consists of 76 variants (https://www.ncbi.nlm.nih.gov/pathogens/refgene) . The blaVIM genes are mostly found in Pseudomonas and Enterobacteriaceae and occur frequently within class 1 integrons with distinct gene cassette arrays (GCAs). The blaVIM-carrying integrons are always presented in various mobile genetic elements (MGEs) such as plasmids, integrative and conjugative elements (ICEs), and integrative and mobilizable elements (IMEs), enhancing the mobility and dissemination of blaVIM genes .There are a wealth of sequenced blaVIM-carrying plasmids, which can be assigned to various incompatibility (Inc) groups such as IncA, IncC, IncN, IncP and IncL/M [3,4,5,6,7]. IMEs/ICEs are frequently located in bacterial chromosomes: ICEs encode self-conjugation modules and thus are able to transfer between cells , whereas IMEs needs a helper conjugative element to achieve intercellular mobility .
Our previous studies showed the detailed genetic characteristics of two novel carbapenemase-encoding MGEs: SIM-encoding plasmid pHN39-SIM  and blaVIM-4-containing ICE Tn6413 . This follow-up study presented the complete sequences of five novel blaVIM-2/-4-carrying MGEs, including two plasmids pJ20133-VIM and p716811-VIM, and two chromosomal IMEs Tn6917 and Tn6918, and one chromosomal ICE Tn6953 from Pseudomonas isolates. These two plasmids were assigned to a novel Inc group IncpSTY. Comprehensive sequence comparisons were then applied to three representative IncpSTY plasmids and 12 chromosome-borne IMEs and ICEs (including the above five ones), providing a deeper understanding of the genetic diversification and evolution of VIM-encoding MGEs in Pseudomonas.
Materials and methods
Bacterial strains and identification
Five clinical isolates P. aeruginosa SE5443, P. putida 716811 and 159349, and P. monteilii J20133 and 918607 were collected from four different Chinese public hospitals (Additional file 1: Table S1). The 16S rRNA genes and the blaVIM genes were detected as described previously .
Genomic DNA extraction, sequencing, and sequence assembly
Bacterial genomic DNA was isolated using the UltraClean Microbial Kit (Qiagen, NW, Germany), and sequenced from a sheared DNA library with average size of 15 kb (ranged from 10 to 20 kb) on a PacBio RSII sequencer (Pacific Biosciences, CA, USA), as well as a paired-end library with an average insert size of 350 bp (ranged from 150 to 600 kb) on a HiSeq sequencer (Illumina, CA, USA) . The paired-end short Illumina reads were used to correct the long PacBio reads utilizing proovread , and then the corrected PacBio reads were assembled de novo utilizing SMARdenovo (https://github.com/ruanjue/smartdenovo).
Sequence annotation and comparison
Genome sequences were annotated by the Rapid Annotation using Subsystem Technology (RAST)  combined with BLASTP/BLASTX/BLASTN , Domain (https://www.ncbi.nlm.nih.gov/Structure/cdd/wrpsb.cgi), and RefSeq database . Annotation of resistance genes, MGEs, and other features were carried out using the online databases including CARD , ResFinder , ISfinder , INTEGRALL , and Tn Number Registry . Multiple and pairwise sequence comparisons were performed using BLASTN. Gene organization diagrams were drawn using Inkscape 1.0 (https://inkscape.org/en/).
Indicated nucleotide sequences were aligned using Clustal Omega 1.2.2 , and then maximum-likelihood phylogenetic trees were constructed from aligned sequences using MEGA X 10.1.8  with a bootstrap iteration of 1000.
Each indicated blaVIM-carrying MGE was transformed from its wild-type isolate into rifampin-resistant P. aeruginosa ATCC 27853 or PAO1, through conjugal transfer or electroporation experiments. Three milliliters of overnight cultures of each donor and recipient bacteria were mixed together, harvested, and resuspended in 80 mL of Brain Heart Infusion (BHI) broth (BD Biosciences). The mixture was spotted on a 1 cm2 hydrophilic nylon membrane filter with a 0.45 µm pore size (Millipore) that was placed on BHI agar (BD Biosciences) plate and then incubated for mating at 26 °C or 37 °C for 12 h to 18 h. Bacteria were washed from filter membrane and spotted on BHI plates, for selecting a blaVIM-carrying transconjugant. 1500 mg/mL rifampin (for ATCC 27853 or PAO1), together with 4 mg/mL meropenem (for blaVIM) was used as transconjugant selection.
Phenotypic assays and multi-locus sequence typing (MLST)
Bacterial antimicrobial susceptibility was tested by BioMérieux VITEK 2, and interpreted as per the 2020 Clinical and Laboratory Standards Institute (CLSI) guidelines . Activity of Amber class A/B/D carbapenemases was determined with a modified CarbaNP test . The sequence types (STs) of Pseudomonas isolates were identified according to the online Pseudomonas MLST schemes (https://pubmlst.org/organisms).
Identification of three STs for Pseudomonas isolates
Three different STs, namely ST129, ST17, and ST639, were identified from the five Pseudomonas isolates (Additional file 1: Table S1), whose complete genome sequences were determined in this work. P. putida 716811 and 159349 belonged to ST129 and ST17, respectively. P. aeruginosa SE5443 belonged to ST639. ST129 was a novel ST of P. putida. The two P. monteilii isolates J20133 and 918607 could not be assigned with any ST types due to lack of P. monteilii MLST scheme.
Proposal of a novel group of IncpSTY plasmids
Two blaVIM-carrying plasmids pJ20133-VIM and p716811-VIM were identified from the complete genome sequences of the two Pseudomonas isolates J20133 and 716811 (Additional file 1: Table S1). A novel IncpSTY group was proposed from a total of seven available fully sequenced single-replicon plasmids (Additional file 1: Table S2; last accessed August 10th 2020) that included the above two plasmids together with five additional ones from GenBank, because these seven plasmids harbored not only homologous repA (replication initiation protein) genes together with its iterons (Fig. 1) but also similar backbone gene organizations (Fig. 2). Since the Inc group of plasmids could be divided according to the homology of repA genes [25,26,27], a phylogenetic tree (Fig. 1) was constructed based on the repA sequences of these seven plasmids, showing that these seven plasmids could be divided into three separately clustering subgroups IncpSTY-1/2/3. As shown by pairwise comparison of repA sequences, plasmids within each subgroup showed ≥ 95% identity, while those from different subgroups displayed ≤ 95% but ≥ 82% identity (Additional file 1: Table S3). pSTY , p716811-VIM, and pHN39-SIM  were the first sequenced plasmids of IncpSTY-1, IncpSTY-2, and IncpSTY-3, respectively, and thus were identified as the references of the corresponding IncpSTY subgroups.
Comparison of three IncpSTY plasmids pSTY, pJ20133-VIM and p716811-VIM
A detailed comparative genomics analysis was applied to three representative IncpSTY plasmids including the IncpSTY-1 plasmids pSTY and pJ20133-VIM and the IncpSTY-2 plasmid p716811-VIM. The IncpSTY-3 plasmid pHN39-SIM was not included herein because it had been detailed characterized in our previous study . These three characterized plasmids varied in size from about 243 to 321 kb with predicted ORFs from 255 to 282 (Table 1). The modular structure of each plasmid was divided into the backbone, and the accessory modules that were acquired DNA regions and inserted at different sites of the backbone (Fig. 2). pSTY, pJ20133-VIM, and p716811-VIM had > 80% nucleotide identity across > 66% of their backbone sequences, while pSTY and pJ20133-VIM shared > 74% of their backbones with > 95% nucleotide identity (Additional file 1: Table S4). These three plasmids shared the key IncpSTY backbone markers: repA together with its iterons, parABC for partition, and multiple conjugal transfer genes (Fig. 2). A 47.9-kb backbone region orf426–to–orf375 was presented in pSTY, but it was deleted from pJ20133-VIM and p716811-VIM due to insertion of an 11.5-kb Tn7-family transposon remnant (see below) and a 21.6-kb orf282–to–orf237 region, respectively, at the same backbone site. The orf282 − to − orf237 region carried IS1491, mvaT (DNA-binding domain of bacterial xenogeneic silencer MvaT), five unnamed genes (encoding nucleoid-associated protein, predicted NTPase, AAA family ATPase, endonuclease, and putative transposase, respectively), and eight genes encoding hypothetical proteins.
pSTY, pJ20133-VIM, and p716811-VIM had totally different profiles of accessory modules (Table 1 and Fig. 2): (i) a 59.7-kb paa (phenylacetic acid degradation) region carrying paa locus and mer (mercury resistance) locus, Tn6734, Tn6603a, and ISPsp16 in pSTY; (ii) a 67.2-kb MDR (multi-drug resistance) region, the 11.5-kb Tn7-family transposon remnant, IS2000:P.ae.I3-like intron, ISPmo2, ISPpu22, ISPmo3, and ISPmo5 in pJ20133-VIM; and (iii) a 29.6-kb MDR region, Tn6735, orf282–to–orf237 region, IS1491, and ISPsy2 in p716811-VIM.
The integration of the 59.7-kb paa region into pSTY led to truncation of orf246–to–orf846 region (as originally observed in pJ20133-VIM) into a 138-bp orf846 remnant (Fig. 2). A presumed prototype paa region was terminally bracketed by two copies of Tn5563 ; the right copy was intact while the left copy was truncated due to insertion of Tn6481 in the 59.7-kb paa region from pSTY (Additional file 1: Fig. S1). This 59.7-kb paa region still contained a complete set of phenylacetic acid degradation genes (Additional file 1: Fig. S1). Tn6734 was an IME containing a phenol degradation gene locus (Additional file 1: Fig. S2a). Tn6603a was a Tn3-family unit transposon containing an oxidative stress defense gene osmC (Additional file 1: Fig. S2b). The 67.2-kb MDR region from pJ20133-VIM contained In58, Tn5046 and ΔTn512, a 3.6-kb Tn4662a remnant, and several IS elements (Additional file 1: Fig. S3).
According to the previous grouping scheme of Tn7-family transposons , the phylogenetic tree was constructed based on the nucleotide sequences of tnsA, which encoded the endonuclease responsible for excision of Tn7-family transposons. A phylogenetic analysis of Tn6735, the 11.5-kb Tn7-family transposon remnant, Tn6922, and Tn6921, together with additional 12 representative sequenced Tn7-family transposons from GenBank based on the tnsA genes, indicated that the first two belonged to a novel Tn6735 subfamily and an unnamable subfamily (it had no intact sequenced transposon), respectively (Fig. 3). Similar to Tn7 , both Tn6735 (carrying pyrimidine biosynthesis genes) and the 11.5-kb Tn7-family transposon remnant encoded the core transposition determinants (endonuclease TnsA, transposase TnsB, regulator TnsC, target-site selection protein TnsD, and TnsB-binding sites), and the later one still encoded an additional target-site selection protein TnsE (Additional file 1: Fig. S4).
The 29.6-kb MDR region from p716811-VIM has a primary structure of a P.p.I3-like intron with interruption by Tn6740 that carried blaVIM-2, aacA4, qnrVC1, and mer locus (Additional file 1: Fig. S5). As a derivant of the IS66-family transposon Tn5042 , Tn6740 retained the whole Tn5042 sequence with integration of a 21.2-kb region that comprised mer-containing Tn3-family unit transposon Tn5563 , qnrVC1, and blaVIM-2-carrying In1722 (see below) (Additional file 1: Fig. S5).
Collection of 12 chromosome-borne IMEs/ICEs for sequence comparison
Three blaVIM-carrying chromosomal MGEs Tn6917, Tn6918 and Tn6953 were identified from the complete genome sequences of the three Pseudomonas isolates 918607, 159349, and SE5443 (Additional file 1: Table S1). A detailed sequence comparison was applied to a total of 12 chromosomal MGEs belonging to three groups: two related IMEs Tn6916 and Tn6917; two related IMEs Tn6918 and Tn6919; and eight related ICEs Tn6417, Tn6413, Tn6953, Tn6454, Tn6455, Tn6456, Tn6457 and Tn6458 (Table 2; last accessed August 10th 2020). Tn6917 was the only blaVIM-carrying one belonging to Tn6916-related IMEs in GenBank. Tn6919 was the only one (in addition to Tn6918) belonging to Tn6918-related IMEs in GenBank, and it did not carry blaVIM. The seven Tn6417-related ICEs were selected from GenBank because Tn6417 was used as reference and the rest six ones Tn6413, Tn6954, Tn6955, Tn6956, Tn6957 and Tn6958 carried blaVIM genes. Similar to the above IncpSTY plasmids, the modular structure of each IME/ICE was divided into the backbone and the accessory modules.
Comparison of two related IMEs Tn6916 and Tn6917
Tn6916 (a 51.0-kb prototype IME initially identified in P. monteilii FDAARGOS_171 ) and Tn6917 were integrated at the same site downstream of the P. monteilii chromosomal gene dinG (ATP-dependent DNA helicase), and they shared the core backbone markers attL (attachment site at the left end), int (integrase), and attR (attachment site at the right end); in addition, the β-lactam resistance gene ampC was considered as a backbone component of these two IMEs because no associated MGEs were identified for this resistance gene (Fig. 4). Tn6916 thereby had no accessory modules, but blaVIM-2-carrying In528 (see below) was integrated at a site upstream of eamA (transcription regulator) and identified as the sole accessory module in Tn6917 (Fig. 4).
Comparison of two related IMEs Tn6918 and Tn6919
Tn6918 (a 70.49-kb prototype IME initially found in P. putida 159,349) and Tn6919 were integrated at a site downstream of the P. putida chromosomal gene tRNAThr, and they shared the core backbone markers attL, int1, int2, and attR (Fig. 5). There were at least three major modular differences in the backbone of Tn6919 relative to Tn6918: i) inversion of orf1965, ii) interruption of orf549 and orf2007 due to insertion of ΔISPpu13 and ISPen2, respectively, and iii) replacement of 3’-termimal 22.8-kb orf270–to–orf213 region by 3.4-kb orf255–to–orf282 region (Fig. 5). The ISPa122–mer region was integrated at the same site of the backbones of Tn6918 and Tn6919 (Fig. 5), and it was composed of ISPa122 and a mer-carrying ΔTn5041-like element (Additional file 1: Fig. S6). Tn6918 and Tn6919 acquired blaVIM-2-carrying In1770 (see below) and strAB region, respectively, which were inserted at the same site upstream of orf555 (Fig. 5). The strAB region contained a cryptic unit transposon Tn6920, a strAB-carrying unit transposon Tn6921, a 3.6-kb Tn4662a remnant, and several IS elements (Additional file 1: Fig. S7). Tn6920 and Tn6921 had the core transposition tniABQ− res −tniR modules with the highest sequence identities to Tn5053  and Tn5058a , respectively, and thus they belonged to the Tn5053 subfamily of Tn7 family (Fig. 3).
Comparison of eight Tn6417-related ICEs
Tn6417 was a 108.2-kb reference ICE  initially found in P. aeruginosa DHS01 . The eight related ICEs Tn6417, Tn6413, Tn6953, Tn6954 (GZAF3_GI) , Tn6955 (ICE6440) , Tn6956 (ICE6441) , Tn6957, and Tn6958  shared the core backbone markers attL/attR (these sequences showed somewhat differences but shared a core motif ‘CCTTCGCCCGCTCCA’), int, cpl (coupling protein), rlx (relaxase), and a F (TivF)-type type IV secretion system gene set (mating pair formation). However, the backbones of these eight ICEs varied in size from 67.8 to 95.1 kb and had at least four major modular differences: (i) presence of unique orf672, orf306, piL1–to–orf381 region, and orf3336–to–orf2514 region in Tn6953; (ii) presence of two highly similar regions lysR2–to–orf1068 and orf348–to–orf1188 in Tn6417 and Tn6413, respectively; (iii) absence of smc–to–orf3006 region from Tn6955; and (iv) interruption of orf2280 owing to insertion of intron E.c.I5 in Tn6954 and Tn6955 (Fig. 6). The first six ICEs were integrated into the P. aeruginosa or A. faecalis chromosomal tRNAGly gene, while the last two ones into the P. aeruginosa chromosomal gene mhpC (methyl ester carboxylesterase) and orf1203 (putative DNA-binding protein), respectively.
Each ICE carried a single accessory module, and thus there were the resulting eight different accessory modules blaVIM-2-carrying In1779 (see below), Tn6532, Tn6403, Tn6959, a 15.2-kb Tn6346-related region, Tn6960, Tn6961, and Tn6962 from the eight ICEs; the first one was inserted into ftsk (cell division protein), while all the other seven ones were integrated at a site upstream of orf582 and identified as the derivatives of Tn6346 (Fig. 7). Tn6346 was a prototype Tn3-family unit transposon originally identified in Achromobacter spp. AO22  and manifested as a hybrid of the core transposition module tnpAR–res from Tn5051 and the mer region from Tn501 (Fig. 7). The seven Tn6346 derivatives differed from Tn6346 in two major aspects: (i) insertion of IS1071 at the same position within tnpA in all these seven Tn6346 derivatives, and furthermore insertion of ISPa91 at another position within tnpA in only Tn6960; and (ii) insertion of different concise class 1 integrons In159, In127, In1365, In56, In1879, In1853, and In58 into urf2 in these seven Tn6346 derivatives, sometimes leading to truncation of surrounding regions (Fig. 7). The blaVIM genes were found the last five integrons (see below) rather than the first two (Additional file 1: Fig. S8).
Comparison of eleven bla VIM-carrying integrons
The blaVIM genes were found in the two plasmids pJ20133-VIM and p716811-VIM (Table 1) and the nine chromosomal IMEs/ICEs Tn6917, Tn6918, Tn6413, Tn6953, Tn6454, Tn6455, Tn6456, Tn6457 and Tn6458 (Table 2). The 11 corresponding local blaVIM genetic environments manifested as 11 different class 1 concise integrons, and three of them contained blaVIM-4 while the other nine carried blaVIM-2 (Fig. 8). A complete class 1 concise integron contained IRi (inverted repeat at the integrase end), 5’-conserved segment (5'-CS: intI1–attI1), GCA, 3'-CS (qacED1–sul1–orf5–orf6), tniTn402 (Tn402 core transposition module tniABQ–res–tniR), and IRt (inverted repeat at the tni end). Except for In1443Tn6413 with terminal truncations, all the other tens had the terminal IRi/IRt pairs. Insertions or truncations occurred in 5'-CS or 3'-CS or tniTn402 of all these 11 integrons, especially including In1779 that showed the longest modular structure because of integration of a mer-containing Tn3-family prototype unit transposon Tn6758. In56Tn6596 had the shortest GCA that contained only blaVIM-2, while all the other 10 integrons had the multiple-gene GCAs harboring blaVIM-2 plus one or more additional antibiotic resistance genes.
Summary of newly identified or designated MGEs
There were eight newly identified MGEs, including two IMEs Tn6917 and Tn6918, one ICE Tn6953, two unit transposons Tn6735 and Tn6740, and three integrons In1770, In1772, and In1779. Additional 13 MGEs (IMEs: Tn6734, Tn6916, Tn6919, and Tn6957; unit transposons: Tn6920, Tn6921, Tn6959, Tn6960, Tn6961, and Tn6962; IS elements: ISPmo2, ISPmo3, ISPmo5) were newly designated but with previously determined sequences. The three previously designated ICEs GZAF3_GI, ICE6440, and ICE6441 were renamed as the standard Tn designations Tn6954, Tn6955, and Tn6956, respectively.
Transferability and antimicrobial susceptibility
pJ20133-VIM, p716811-VIM, and Tn6953, which were selected to represent the IncpSTY-1 and IncpSTY-2 plasmids and the Tn6417-related ICEs, respectively, could be transferred from the relevant wild-type isolates into ATCC 27853 or PAO1 through conjugation, generating the transconjugants PAO1/pJ20133-VIM, PAO1/p716811-VIM, and ATCC 27853/Tn6953 respectively. All these wild-type and transconjugant strains had the Ambler class B carbapenemase activity (data not shown) and were resistant to imipenem and meropenem with minimal inhibit concentration (MIC) values ≥ 16 mg/L (Additional file 1: Table S5), owing to production of VIM enzyme.
The present work identifies a novel group of IncpSTY plasmids, which can be further divided into the three subgroups IncpSTY-1/2/3. A detailed genetic dissection analysis is applied to a collection of 15 Pseudomonas MGEs, including three IncpSTY plasmids, two Tn6916-related IMEs, two Tn6918-related IMEs, and eight Tn6417-related ICEs. All these IMEs/ICEs were located within the bacterial chromosomes. At least 22 genes for resistance to seven different categories of antibiotics and heavy metals are identified within these 15 MGEs (Additional file 1: Table S6). Eleven of these 15 MGEs carry the blaVIM genes, which are all located within GCAs of concise class 1 integrons. For 10 of these integrons, the blaVIM genes were presented together with other resistance genes, especially including those for resistance to β-lactams, aminoglycosides, trimethoprim, and rifampicin. These blaVIM-carrying integrons were further integrated into the above plasmids, IMEs and ICEs.
Confirmed by conjugal transfer experiments herein, IncpSTY plasmids and Tn6417-related ICEs have the ability to transfer from one cell to another cell. These two groups MGEs could transfer autonomously by utilizing self-encoded conjugation gene sets [39, 40]. However, the conjugal transfer of Tn6916- and Tn6918-related IMEs are nonautonomous due to lack of essential conjugation genes encoding Cpl and type IV secretion system, and thereby its intercellular transfer is relied on the help of other conjugative elements .
blaVIM genes are widely disseminated in Pseudomonas [41,42,43], Enterobacteriaceae [5, 44, 45], Providencia [46, 47], Achromobacter , Acinetobacter [49, 50], Aeromonas , and Morganella . Currently, at least 183 different blaVIM-carrying integrons [19, 53, 54] have been identified, and they either directly integrate into the chromosomes or reside as the inner-components of various MGEs such as plasmids, unit transposons, ICEs and IMEs [3, 55,56,57]. These results indicate that blaVIM genes have evolved to reside diverse MGEs with the intra- or inter-cellular mobility, which will promote the wide dissemination of blaVIM genes.
Until now, only seven IncpSTY plasmids are identified in six Pseudomonas isolates and one Stenotrophomonas rhizophila isolate, which come from China (n = 5), Japan (n = 1) and Germany (n = 1) (Additional file 1: Table S2). Only six Tn6916-related IMEs and two Tn6918-related IMEs are detected in eight Pseudomonas isolates, two of which are sequenced in this study. The six Tn6916-related IMEs come from China (n = 4), USA (n = 1) and Japan (n = 1), while the two Tn6918-related IMEs from China (n = 1) and Brazil (n = 1). These observations indicate that these three group MGEs are mainly prevalent in Pseudomonas isolates and meanwhile still not widely disseminated in multiple geographic regions. By contrast, Tn6417-related ICEs are frequently identified in Pseudomonas [7, 58, 59] and also occasionally in Bordetella, Alcaligenes, Klebsiella pneumoniae, Achromobacter xylosoxidans, Morganella morganii, Aeromonas caviae, Casimicrobium huifangae, and Delftia acidovorans. Tn6417-related ICEs have been already disseminated worldwide in a plenty of bacterial species.
The four groups of blaVIM-carrying MGEs characterized in this study have a wide range of hosts including Pseudomonas, and they are able to transfer across different bacterial species. Integration of blaVIM genes into these MGEs contributes to the accumulation and distribution of blaVIM genes and enhances the ability of bacteria to survive under selection pressure of carbapenems especially in hospital settings. Moreover, these four groups of MGEs display a high-level diversification in modular structures, which have complex mosaic natures and carry a large number of drug resistance genes (particularly in the MDR regions in these MGEs). This study provides a deeper insight into the genetic diversification and evolution of VIM-encoding MGEs in Pseudomonas.
Availability of data and materials
The complete sequences of plasmids pJ20133-VIM and p716811-VIM, and those of the chromosomes of strains SE5443, 918607, and 159349 were submitted to GenBank under accession numbers MN310371, MN310372, CP046405, CP043395, and CP045553 respectively.
Mobile genetic element
Integrative and conjugative element
Integrative and mobilizable element
Gene cassette array
Rapid Annotation using Subsystem Technology
Brain Heart Infusion
Clinical and Laboratory Standards Institute
Multi-locus Sequence Typing
Inverted repeat at the integrase end
Inverted repeat at the tni end
Minimal inhibit concentration
Queenan AM, Bush K. Carbapenemases: the versatile beta-lactamases. Clin Microbiol Rev. 2007;20:440–58.
Toleman MA, Vinodh H, Sekar U, et al. blaVIM-2-harboring integrons isolated in India, Russia, and the United States arise from an ancestral class 1 integron predating the formation of the 3’ conserved sequence. Antimicrob Agents Chemother. 2007;51:2636–8.
Botelho J, Grosso F, Quinteira S, et al. Two decades of blaVIM-2-producing Pseudomonas aeruginosa dissemination: an interplay between mobile genetic elements and successful clones. J Antimicrob Chemother. 2018;73:873–82.
Gaibani P, Ambretti S, Scaltriti E, et al. A novel IncA plasmid carrying blaVIM-1 in a Kluyvera cryocrescens strain. J Antimicrob Chemother. 2018;73:3206–8.
Miriagou V, Papagiannitsis CC, Kotsakis SD, et al. Sequence of pNL194, a 79.3-kilobase IncN plasmid carrying the blaVIM-1 metallo-beta-lactamase gene in Klebsiella pneumoniae. Antimicrob Agents Chemother. 2010;54:4497–502.
Peirano G, Lascols C, Hackel M, et al. Molecular epidemiology of Enterobacteriaceae that produce VIMs and IMPs from the SMART surveillance program. Diagn Microbiol Infect Dis. 2014;78:277–81.
Zeng L, Zhan Z, Hu L, et al. Genetic characterization of a blaVIM-24-carrying IncP-7β plasmid p1160-VIM and a blaVIM-4-harboring integrative and conjugative element Tn6413 from clinical Pseudomonas aeruginosa. Front Microbiol. 2019;10:213.
Botelho J, Schulenburg H. The role of integrative and conjugative elements in antibiotic resistance evolution. Trends Microbiol. 2021;29:8–18.
Guedon G, Libante V, Coluzzi C, et al. The obscure world of integrative and mobilizable elements, highly widespread elements that pirate bacterial conjugative systems. Genes (Basel). 2017;8:337.
Sun F, Zhou D, Wang Q, et al. The first report of detecting the blaSIM-2 gene and determining the complete sequence of the SIM-encoding plasmid. Clin Microbiol Infect. 2016;22:347–51.
Fu J, Zhang J, Yang L, et al. Precision methylome and in vivo methylation kinetics characterization of Klebsiella pneumoniae. GPB. 2021. https://doi.org/10.1016/j.gpb.2021.04.002.
Hackl T, Hedrich R, Schultz J, et al. Proovread: large-scale high-accuracy PacBio correction through iterative short read consensus. Bioinformatics. 2014;30:3004–11.
Brettin T, Davis JJ, Disz T, et al. RASTtk: a modular and extensible implementation of the RAST algorithm for building custom annotation pipelines and annotating batches of genomes. Sci Rep. 2015;5:8365.
Boratyn GM, Camacho C, Cooper PS, et al. BLAST: a more efficient report with usability improvements. Nucleic Acids Res. 2013;41:W29-33.
O’Leary NA, Wright MW, Brister JR, et al. Reference sequence (RefSeq) database at NCBI: current status, taxonomic expansion, and functional annotation. Nucleic Acids Res. 2016;44:D733–45.
Jia B, Raphenya AR, Alcock B, et al. CARD 2017: expansion and model-centric curation of the comprehensive antibiotic resistance database. Nucleic Acids Res. 2017;45:D566–73.
Zankari E, Hasman H, Cosentino S, et al. Identification of acquired antimicrobial resistance genes. J Antimicrob Chemother. 2012;67:2640–4.
Siguier P, Perochon J, Lestrade L, et al. ISfinder: the reference centre for bacterial insertion sequences. Nucleic Acids Res. 2006;34:D32–6.
Moura A, Soares M, Pereira C, et al. INTEGRALL: a database and search engine for integrons, integrases and gene cassettes. Bioinformatics. 2009;25:1096–8.
Roberts AP, Chandler M, Courvalin P, et al. Revised nomenclature for transposable genetic elements. Plasmid. 2008;60:167–73.
Sievers F, Higgins DG. Clustal Omega for making accurate alignments of many protein sequences. Protein Sci. 2018;27:135–45.
Kumar S, Stecher G, Li M, et al. MEGA X: molecular evolutionary genetics analysis across computing platforms. Mol Biol Evol. 2018;35:1547–9.
CLSI. Performance standards for antimicrobial susceptibility testing. Wayne: Clinical and Laboratory Standards Institute; 2020.
Sun F, Yin Z, Feng J, et al. Production of plasmid-encoding NDM-1 in clinical Raoultella ornithinolytica and Leclercia adecarboxylata from China. Front Microbiol. 2015;6:458.
Carloni E, Andreoni F, Omiccioli E, et al. Comparative analysis of the standard PCR-Based Replicon Typing (PBRT) with the commercial PBRT-KIT. Plasmid. 2017;90:10–4.
Carattoli A, Bertini A, Villa L, et al. Identification of plasmids by PCR-based replicon typing. J Microbiol Methods. 2005;63:219–28.
Partridge SR, Kwong SM, Firth N, et al. Mobile genetic elements associated with antimicrobial resistance. Clin Microbiol Rev. 2018. https://doi.org/10.1128/CMR.00088-17.
Kohler KA, Ruckert C, Schatschneider S, et al. Complete genome sequence of Pseudomonas sp. strain VLB120 a solvent tolerant, styrene degrading bacterium, isolated from forest soil. J Biotechnol. 2013;168:729–30.
Yeo CC, Tham JM, Kwong SM, et al. Tn5563, a transposon encoding putative mercuric ion transport proteins located on plasmid pRA2 of Pseudomonas alcaligenes. FEMS Microbiol Lett. 1998;165:253–60.
Peters JE, Fricker AD, Kapili BJ, et al. Heteromeric transposase elements: generators of genomic islands across diverse bacteria. Mol Microbiol. 2014;93:1084–92.
Peters JE, Craig NL. Tn7: smarter than we thought. Nat Rev Mol Cell Biol. 2001;2:806–14.
Mindlin S, Minakhin L, Petrova M, et al. Present-day mercury resistance transposons are common in bacteria preserved in permafrost grounds since the Upper Pleistocene. Res Microbiol. 2005;156:994–1004.
Sichtig H, Minogue T, Yan Y, et al. FDA-ARGOS is a database with public quality-controlled reference genomes for diagnostic use and regulatory science. Nat Commun. 2019;10:3313.
Kholodii GY, Yurieva OV, Lomovskaya OL, et al. Tn5053, a mercury resistance transposon with integron’s ends. J Mol Biol. 1993;230:1103–7.
Valot B, Rohmer L, Jacobs MA, et al. Comparative genomic analysis of two multidrug-resistant clinical isolates of ST395 epidemic strain of Pseudomonas aeruginosa obtained 12 years apart. Genome Announc. 2014;2:e00515-e614.
Al Laham N, Chavda KD, Cienfuegos-Gallet AV, et al. Genomic characterization of VIM metallo-β-Lactamase-producing Alcaligenes faecalis from Gaza. Palestine Antimicrob Agents Chemother. 2017;61:e01499-e1517.
Majewski P, Majewska P, Gutowska A, et al. Molecular characterisation of clinical pandrug-resistant Alcaligenes faecalis strain MUB14. Int J Antimicrob Agents. 2020;55:105939.
Ng SP, Davis B, Palombo EA, et al. A Tn5051-like mer-containing transposon identified in a heavy metal tolerant strain Achromobacter sp. AO22. BMC Res Notes. 2009;2:38.
Smillie C, Garcillán-Barcia MP, Francia MV, et al. Mobility of plasmids. Microbiol Mol Biol Rev. 2010;74:434–52.
Bellanger X, Payot S, Leblond-Bourget N, et al. Conjugative and mobilizable genomic islands in bacteria: evolution and diversity. FEMS Microbiol Rev. 2014;38:720–60.
Li R, Peng K, Xiao X, et al. Emergence of a multidrug resistance efflux pump with carbapenem resistance gene blaVIM-2 in a Pseudomonas putida megaplasmid of migratory bird origin. J Antimicrob Chemother. 2021;76:1455–8.
Hishinuma T, Uchida H, Tohya M, et al. Emergence and spread of VIM-type metallo-β-lactamase-producing Pseudomonas aeruginosa clinical isolates in Japan. J Glob Antimicrob Resist. 2020;23:265–8.
Ocampo-Sosa AA, Guzmán-Gómez LP, Fernández-Martínez M, et al. Isolation of VIM-2-producing Pseudomonas monteilii clinical strains disseminated in a tertiary hospital in northern Spain. Antimicrob Agents Chemother. 2015;59:1334–6.
Matsumura Y, Peirano G, Devinney R, et al. Genomic epidemiology of global VIM-producing Enterobacteriaceae. J Antimicrob Chemother. 2017;72:2249–58.
Pournaras S, Ikonomidis A, Tzouvelekis LS, et al. VIM-12, a novel plasmid-mediated metallo-beta-lactamase from Klebsiella pneumoniae that resembles a VIM-1/VIM-2 hybrid. Antimicrob Agents Chemother. 2005;49:5153–6.
Miriagou V, Tzouvelekis LS, Flevari K, et al. Providencia stuartii with VIM-1 metallo-beta-lactamase. J Antimicrob Chemother. 2007;60:183–4.
Lee HW, Kang HY, Shin KS, et al. Multidrug-resistant Providencia isolates carrying blaPER-1, blaVIM-2, and armA. J Microbiol. 2007;45:272–4.
Shin KS, Han K, Lee J, et al. Imipenem-resistant Achromobacter xylosoxidans carrying blaVIM-2-containing class 1 integron. Diagn Microbiol Infect Dis. 2005;53:215–20.
Lee MF, Peng CF, Hsu HJ, et al. Molecular characterisation of the metallo-beta-lactamase genes in imipenem-resistant Gram-negative bacteria from a university hospital in southern Taiwan. Int J Antimicrob Agents. 2008;32:475–80.
Figueiredo S, Poirel L, Papa A, et al. First identification of VIM-4 metallo-beta-lactamase in Acinetobacter spp. Clin Microbiol Infect. 2008;14:289–90.
Libisch B, Giske CG, Kovács B, et al. Identification of the first VIM metallo-beta-lactamase-producing multiresistant Aeromonas hydrophila strain. J Clin Microbiol. 2008;46:1878–80.
Tsakris A, Ikonomidis A, Spanakis N, et al. Characterization of In3Mor, a new integron carrying VIM-1 metallo-beta-lactamase and sat1 gene, from Morganella morganii. J Antimicrob Chemother. 2007;59:739–41.
Papagiannitsis CC, Verra A, Galani V, et al. Unravelling the features of success of VIM-producing ST111 and ST235 Pseudomonas aeruginosa in a Greek Hospital. Microorganisms. 2020. https://doi.org/10.3390/microorganisms8121884.
Papagiannitsis CC, Medvecky M, Chudejova K, et al. Molecular characterization of carbapenemase-producing Pseudomonas aeruginosa of czech origin and evidence for clonal spread of extensively resistant sequence type 357 expressing IMP-7 metallo-β-Lactamase. Antimicrob Agents Chemother. 2017. https://doi.org/10.1128/AAC.01811-17.
Zhao W, Hu Z. Epidemiology and genetics of VIM-type metallo-β-lactamases in Gram-negative bacilli. Future Microbiol. 2011;6:317–33.
Tseng SP, Hsueh PR, Tsai JC, et al. Tn6001, a transposon-like element containing the blaVIM-3-harboring integron In450. Antimicrob Agents Chemother. 2007;51:4187–90.
Di Pilato V, Pollini S, Rossolini GM. Tn6249, a new Tn6162 transposon derivative carrying a double-integron platform and involved with acquisition of the blaVIM-1 metallo-β-lactamase gene in Pseudomonas aeruginosa. Antimicrob Agents Chemother. 2015;59:1583–7.
Valot B, Rohmer L, Jacobs MA, et al. Comparative genomic analysis of two multidrug-resistant clinical isolates of ST395 epidemic strain of pseudomonas aeruginosa obtained 12 years apart. Genome Announc. 2014. https://doi.org/10.1128/genomeA.00515-14.
Yu T, Yang H, Li J, et al. Novel chromosome-borne accessory genetic elements carrying multiple antibiotic resistance genes in Pseudomonas aeruginosa. Front Cell Infect Microbiol. 2021;11:638087.
Pansegrau W, Lanka E, Barth PT, et al. Complete nucleotide sequence of Birmingham IncP alpha plasmids. Compilation and comparative analysis. J Mol Biol. 1994;239:623–63.
This work was supported by the National Natural Science Foundation of China (Grant No. 3214100027).
Ethics approval and consent to participate
Consent for publication
The authors declare that they have no competing interests.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Additional file 1: Fig S1.
Organization of paa region from pSTY. Fig S2. Organization of Tn6734 and Tn6603a from pSTY. Fig S3. Organization of MDR region from pJ20133-VIM, and comparison with related regions. Fig S4. Organization of two Tn7-family elements, and comparison with related region. Fig S5. Organization of MDR region from p716811-VIM, and comparison with related regions. Fig S6. Organization of ISPa122 − mer region from Tn6918 and Tn6919, and comparison with related region. Fig S7. Organization of strAB region from Tn6919, and comparison with related regions. Fig S8. Organization of In159 and In127 from Tn6532 and Tn6403 respectively, and comparison with related region. Table S1. Major features of the five Pseudomonas isolates sequenced in this study. Table S2. List of all the seven sequenced IncpSTY single-replicon plasmids. Table S3. Pairwise comparison of repAIncpSTY sequences from the seven IncpSTY plasmids using BLASTN. Table S4. Pairwise comparison of backbone sequences from the three IncpSTY plasmids using BLASTN. Table S5. Antimicrobial drug susceptibility profiles. Table S6. Resistance genes in the 15 mobile genetic elements (MGEs) analyzed.
Rights and permissions
Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecommons.org/licenses/by/4.0/. 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 in a credit line to the data.
About this article
Cite this article
Chen, F., Wang, P., Yin, Z. et al. VIM-encoding IncpSTY plasmids and chromosome-borne integrative and mobilizable elements (IMEs) and integrative and conjugative elements (ICEs) in Pseudomonas. Ann Clin Microbiol Antimicrob 21, 10 (2022). https://doi.org/10.1186/s12941-022-00502-w
- Bla VIM
- Mobile genetic elements
- Antibiotic resistance