- Research article
- Open Access
Genome-based characterization of two Colombian clinical Providencia rettgeri isolates co-harboring NDM-1, VIM-2, and other β-lactamases
BMC Microbiology volume 20, Article number: 345 (2020)
Providencia rettgeri is a nosocomial pathogen associated with urinary tract infections and related to Healthcare-Associated Infection (HAI). In recent years isolates producing New Delhi Metallo-β-lactamase (NDM) and other β-lactamases have been reported that reduce the efficiency of clinical antimicrobial treatments. In this study, we analyzed antibiotic resistance, the presence of resistance genes and the clonal relationship of two P. rettgeri isolates obtained from male patients admitted to the same hospital in Bogotá – Colombia, 2015.
Antibiotic susceptibility profile evaluated by the Kirby-Bauer method revealed that both isolates were resistant to third-generation carbapenems and cephalosporins.
Whole-genome sequencing (Illumina HiSeq) followed by SPAdes assembling, Prokka annotation in combination with an in-house Python program and resistance gene detection by ResFinder identified the same six β-lactamase genes in both isolates: blaNDM-1, blaVIM-2, blaCTX-M-15, blaOXA-10, blaCMY-2 and blaTEM-1. Additionally, various resistance genes associated with antibiotic target alteration (arnA, PmrE, PmrF, LpxA, LpxC, gyrB, folP, murA, rpoB, rpsL, tet34) were found and four efflux pumps (RosAB, EmrD, mdtH and cmlA).
The additional resistance to gentamicin in one of the two isolates could be explained by a detected SNP in CpxA (Cys191Arg) which is involved in the stress response of the bacterial envelope.
Genome BLAST comparison using CGView, the ANI value (99.99%) and the pangenome (using Roary) phylogenetic tree (same clade, small distance) showed high similarity between the isolates. The rMLST analysis indicated that both isolates were typed as rST-61,696, same as the RB151 isolate previously isolated in Bucaramanga, Colombia, 2013, and the FDAARGOS_330 isolate isolated in the USA, 2015.
We report the coexistence of the carbapenemase genes blaNDM-1, and blaVIM-2, together with the β-lactamase genes blaCTX-M-15, blaOXA-10, blaCMY-2 and blaTEM-1, in P. rettgeri isolates from two patients in Colombia. Whole-genome sequence analysis indicated a circulation of P. rettgeri rST-61,696 strains in America that needs to be investigated further.
Providencia rettgeri is an opportunistic pathogen mainly associated with urinary tract infections ; it can be a causative agent of bacteremia , traveler’s diarrhea , eye infections  and wounds  and in a low incidence, of meningitis  and pyelonephritis  . In recent years, P. rettgeri has taken on major importance due to the emergence of carrier isolates of the New Delhi metallo-β-lactamase NDM-1, a carbapenemase that hydrolyzes all β-lactams for clinical use except aztreonam [8, 9]. NDM-1 was first described in 2009, in a Klebsiella pneumoniae isolated in a Swedish hospital from a patient who had been previously hospitalized in India . Since then, P. rettgeri isolates producing NDM-1 have been reported in different parts of the world, including Asia (Nepal, Israel, India, Korea, China, Pakistan) [11,12,13,14,15,16] South America (Brazil, Argentina, Colombia, Ecuador) [17,18,19,20,21] North America (Canada, Mexico) [22, 23] South Africa (Nigeria)  and South East Europe (Bulgaria) .
The first report of P. rettgeri producing NDM-1 in Colombia was published in 2015. The isolates were obtained in 2012–2013, in the departments of Santander and Cundinamarca, and in the capital city, Bogotá . In 2016, the Colombian Antimicrobial Resistance Surveillance Program in Healthcare-Associated Infections (HAI) reported P. rettgeri as one of the most commonly sent of the order Enterobacterales for confirmation of carbapenemases , and a study conducted by the National Institute of Health (INS) of Colombia in 2018 (unpublished data) P. rettgeri was determined to be the second most common carrier of the blaNDM-1 gene among Enterobacteriaceae in this country.
The presence of blaNDM-1 in a bacterial isolate is a challenge for clinical therapeutics. However, it is a greater challenge when it is co-harboring with other resistance determinants because it results in a greater limitation of the treatment options [27,28,29]. Worldwide, few reports of P. rettgeri notifying coexistence of blaNDM-1 with other resistance determinants have been published. In 2014, the first coexistence of metallo-β-lactamase blaNDM-1 was reported along with armA encoding a 16S rRNA methylase in clinical isolates of P. rettgeri in Nepal . Recently, the coexistence of blaNDM-1, along with other resistance determinants such as blaPER-1-type extended-spectrum β-lactamases (ESBLs) in Korea , and class D and A β-lactamases blaOXA-48 and blaTEM-1, respectively, in Turkey  have been described. Similarly, other coexistence in P. rettgeri that do not include the blaNDM-1 gene have been described. In 2007, the co-production of PER-1, VIM-2 and ArmA in Providencia spp. isolates in Korea was reported ; and in Nigeria, an isolation of P. rettgeri was reported in 2011 carrying genes blaOXA-10, blaVEB-1, blaCMY-4 and blaTEM-1 .
In this study, we used whole-genome sequencing (WGS) for characterization of two P. rettgeri isolates that carried blaNDM-1 and blaVIM-2 carbapenemases, recovered in March and December 2015 at a fourth-level health care institution in Colombia.
Amplification and sequencing of 16S rRNA
The 16S rRNAs obtained by PCR and sequenced for both isolates (1490 pb each of them), were aligned against the NCBI DNA nr database using blastn. The best hits of isolates GMR-RA257 and GMR-RA1153, corresponded to the seqref 16S RNA of P. rettgeri strain RB-151 (coverage 100%, value E 0.0, identity 100%) and P. rettgeri strain FDAARGOS_330 (coverage 100%, value E 0.0, identity 99%), which confirmed their species identity as Providencia rettgeri.
Antibiotic resistance of P. rettgeri isolates
The isolates GMR-RA257 and GMR-RA1153 presented a similar susceptibility profile. Both isolates were resistant to carbapenems (imipinem, meropenem), second- (cefoxitin) and third-generation cephalosporins (ceftazidime, cefotaxime, cefepime), amikacin and trimethoprim/sulfamethoxazole and were susceptible to aztreonam and piperacillin/tazobactam.
Some differences in the phenotypic profile of susceptibility were observed: GMR-RA257 showed an intermediate susceptibility to gentamicin, while GMR-RA1153 showed intermediate susceptibility to ertapenem. The phenotypic susceptibility profile of the P. rettgeri isolates showed resistance to three or more families of antibiotics, and therefore, both were classified as multi-drug resistant (MDR) according to standardized international terminology .
The imipenem-EDTA/SMA double-disk synergy test confirmed the production of enzymes of the metallo-β-lactamase (MBL) type. However, the modified Hodge test (MHT) was negative and the 3-Aminophenylboronic acid (APB)-based disk/microdilution test did not detect the production of serine-type β-lactamases.
Assembly and annotation of GMR-RA257 and GMR-RA1153 genomic sequences
Whole genome sequencing (WGS) of GMR-RA257 and GMR-RA1153 yielded 20,13 and 5,93 millions of reads respectively. These results were obtained using Illumina HiSeq 2500, and quality control was provided using FastQC, obtaining a high-quality score of Q30 (results not shown). The G + C content for these isolates was 40.5%, considered common for this species [9, 34, 35]. Detailed sequencing characteristics of GMR-RA257 and GMR-RA1153 can be found in Table S1. The estimated size of GMR-RA257 and GMR-RA1153 were 4,84 Mb for both isolates, with a greater coverage at a depth of 744X and 216X, respectively. A total of 71 tRNAs, 7 rRNAs and 4452 coding sequence (CDS) were annotated in both genomes. Those values are in the range of the previously reported genomes of P. rettgeri [9, 36]. The whole-genome sequences of P. rettgeri GMR-RA257 and GMR-RA1153 were deposited in the DDBJ/ENA/GenBank databases under the accession numbers VRPG00000000 and VRPH00000000 respectively. Details of the assembly and annotation for Providencia rettgeri isolates GMR-RA257 and GMR-RA1153 can be found in the supplementary Table S1.
The GMR-RA1153 and GMR-RA257 antimicrobial resistomes were obtained from their GFF files produced by Prokka v1.13. The resistomes included genes three different antibiotic resistance categories: antibiotic degradation or modification; antibiotic target alteration; and efflux pumps (Table 1). There was no difference in the type and the numbers of resistance genes present in both isolates: they showed a 100% identity with genes previously reported in the GenBank. Only the CpxA gene (histidine kinase sensor of the two-component system Cpx-TCS, involved in the stress response of the bacterial envelope) of GMR-RA1153 had a 99.93% identity with genes previously reported in GenBank and its alignment with the CpxA gene of GMR-RA257 showed a change in one amino acid (Cys191Arg).
The further results showed that plasmid-mediated genes associated with resistance to β-lactam antibiotics were the most abundant in the antibiotic degradation or modification category. This includes MBL genes encoding blaNDM-1 and blaVIM-2, the ESBL gene blaCTX-M-15, the AmpC gene blaCMY-2, and other β-lactamases genes like blaOXA-10, and blaTEM-1. Further commonly detected genes were associated with aminoglycoside resistance (ant(2″)-Ia, aph(3′)-VI, aac(6′)-Il, and aadA1), and sulfonamide resistance (sul1 and sul2) (Table 1).
The arnA, PmrE, PmrF, LpxA, and LpxC genes, associated with resistance to peptide antibiotics (bacitracin, colistin and polymyxin B) were the most abundant in the antibiotic target alteration category. Additionally, six different efflux pumps were identified, belonging to the major facilitator superfamily (MFS) and represented by pumps of the resistance-nodulation-cell division (RND) family (Table 1).
The carbapenemase genes blaNDM-1 and blaVIM-2 were found in two different contigs. Within the blaNDM-1 contig (7794 pb) was present the IS91 family transposase gene that could increase the blaNDM-1 expression level [41, 42]. This contig had 100% identity and 100% coverage with bigger plasmids present in many other bacteria like in Escherichia coli (GI: MN604268.1), Acinetobacter baumannii (GI: CP027528.1) and P. rettgeri strain RB151(GI: CP017672.1). The blaVIM-2 contig (2011 pb) had 100% identity and coverage within 68 and 79% with bigger plasmids present in other bacteria like in Pseudomona aeruginosa (GI: KR337992.1), Acinetobacter baumannii (GI: AF324464.1) and Pseudomona putida (GI: AF327064.1). Within the small blaVIM-2 contig was present the gene class 1 integron integrase IntI1 associated with antimicrobial resistance gene cassette. The plasmid location of the blaVIM-2 and the other plasmid resistance-associated genes (Table 1) was supported too by the high deep of coverage of their contigs in comparison with the mean coverage of the genome contigs, and the plasmid contig identification made using PlasmiSPAdes (SPAdes v3.8.0)  and PlasFlow 1.1 .
Ribosomal multilocus sequence typing (rMLST)
The rMLST was done using the PubMLST website https://pubmlst.org/rmlst/ for the genomes of the GMR-RA257 and GMR-RA1153 study and 16 Providencia rettgeri genomes available in the GenBank database. Ten of the 18 P. rettgeri genomes evaluated were assigned to seven different rST types, as follows: rST-61,696 was the predominant rST corresponding to the genomes of the study GMR-RA257, GMR-RA1153 and two genomes: RB151 isolated in Bucaramanga, Colombia, in 2013 and FDAARGOS_330 isolated in Washington D.C., the United States in 2015, both of them recovered from urine samples. The rST-37,417, rST-37,410, rST-63,059, rST-63,070, rST-37,423 and rST-37,424 were assigned to an individual reference’s genomes (Table 2).
Whole genome comparisons
The pangenome analysis carried out by Roary identified a total of 12, 539 genes, including 1, 217 (~ 10%) for the core genome, in the 18 genomes used in this work. The 11,322 non-central genes were divided into 6980 accessory genes and 4342 unique genes.
A phylogenetic tree was built from the core and accessory genes presence/absence information produced by Roary, using the isolates and the GenBank genomes above described, along with the rST sequence type, the country and year of collection (Fig. 1). Five important clades were observed, where the different rST types were grouped as follows: rST-37,417 (clade A); rST-37,410, rST-63,059, rST-63,070 and two genomes with no rST known (clade B); rST-37,423, rST-37,424, rST-61,696 and one genome with no rST known (clade C); two genomes with no rST known (clade D); and three genomes with no rST known (clade E) (Fig. 1). The phylogenetic tree shows a possible clonal relationship between GMR-RA257 and GMR-RA1153 which were assigned to clade C; both isolates were recovered in Bogotá and were isolated from two patients in the same fourth level hospital in two different times (March and December, respectively). This may indicate that a clone could be circulating within this hospital during that period of time.
The average nucleotide identity (ANI) was calculated using JSpecies. JSpecies provided a 99.99% ANI between GMR-RA257 and GMR-RA1153. Likewise, the ANI value between these genomes and the Colombian reference genome RB151 provided a value ANI of 99.67%, these results showed a close relationship between the Colombian genomes.
In addition to that, a circular comparison, between the genomic sequences of GMR-RA257 and GMR-RA1153 with the Colombian genome RB151, was performed using CGView Comparison Tool (Fig. 2). The results showed that the three Colombian samples were very similar; most of the compared genomic regions of GMR-RA257 and GMR-RA1153 showed a high genomic similarity (78.1%) with RB151 (Fig. 2).
As result of the Resfinder analysis of the 18 genomes of P. rettgeri, the genomes pret_2032 and CCBH11880 have the highest number of genomic elements associated with resistance to different categories of antibiotics (Fig. 3). GMR-RA1153 and GMR-RA257 also have several resistance genes associated with eight different categories of antibiotics (results similar to those found in the Prokka annotation), but when compared with the other genomes, they showed the highest number of β-lactamases genes associated with resistance to β-lactam antibiotics (Fig. 3). The comparison of the GMR-RA257 and GMR-RA1153 genomes with the genome of P. rettgeri RB151, isolated two years earlier in Colombia, revealed that the three isolates share six different genes associated with resistance: blaNDM-1, blaTEM-1, sul1, sul2, dfrA6, tet . However, in GMR-RA257 and GMR-RA1153 additional genes were detected that were not present in RB151: blaOXA-10, blaCMY-2, blaCTX-M-15, blaVIM-2, qnrD1, aadA1, aph (3′)-VI, aac (6′)-Il, cmlA1, ARR-2. The detailed results of all the antibiotic resistance genes in each class of antibiotics of each reference genome and collection strains are included in additional Table S2.
We focused our analysis on identifying non-synonymous single nucleotide polymorphisms (nsSNPs) that result in amino acid changes within the functional gene sequences in GMR-RA257 and GMR-RA1153. We established 108 non-synonymous single nucleotide variants (nsSNVs) for GMR-RA257 and GMR-RA1153 using SnpEff. From those, 107 variants were located in genes encoding hypothetical proteins. Only one variant (Cys191Arg) present in the antibiotic resistance gene CpxA could explain the increased resistance of GMR-RA1153 to the aminoglycoside antibiotic gentamicin (Table 1).
In this study, we used WGS data to provide information on two isolates of Providencia rettgeri. We detected an uncommon accumulation of various β-lactamases genes (blaCTX-M-15, blaCMY-2, blaOXA-10, and blaTEM-1) genes including carbapenemase genes blaNDM-1 and blaVIM-2. To the best of our knowledge this combination of β-lactamases has not been reported before in Providencia spp. The isolates showed a similar resistance profile, were resistant to most of the antibiotics evaluated, except aztreonam and piperacillin/tazobactam.
Reports in America [17,18,19, 21,22,23] as well as reports in Asia [11,12,13,14,15,16] and in South Africa , commonly associate P. rettgeri with high resistance rates to β-lactamic antibiotics, a characteristic that was conserved in Colombian isolates GMR-RA257 and GMR-RA1153. This resistance trait in P. rettgeri is frequently associated with the production of metallo-β-lactamase (blaNDM-1). GMR-RA257 and GMR-RA1153 were susceptible to aztreonam, similarly as previously reported in Israel, Argentina and in a former Colombian report [12, 19, 20]. It is interesting that the presence of blaCTX-M-15 gene in both isolates did not result in resistance to aztreonam, possibly due to a low expression of this gene. This agrees with different studies where many CTX-M-type ESBL-producing isolates that would previously have been reported as resistant to cephalosporins and aztreonam were now reported as susceptible to these drugs [41, 42, 44]. In isolates GMR-RA257 and GMR-RA1153 the blaCTX-M-15 was located in a small contig in both isolated (2013 bp and 1924 bp) flanked upstream by a cupin-fold metalloprotein WbuC gene and downstream, by a fragment of transposase gene. Both blaCTX-M-15 contigs had 92% coverage and 100% identity with bigger plasmids present in different bacteria like E. coli (GI: MN158990.1) and K. pneunoniae (GI: CP050364.1), and therefore, they could be part of a plasmid ensembled by SPAdes in many contigs, not located into the same scaffold. However, due to the small contig size no assessment of regulatory elements of blaCTX-M-15 could be made for the two study isolates.
The rMLST analysis of GMR-RA257 and GMR-RA1153 and 16 previously reported genomes, allows to identify seven different rST for ten genomes analyzed. GMR-RA257 and GMR-RA1153 belong to rST-61,696, which is the most predominant rST among the genomes used in this work. An interesting finding was that the genomes GMR-RA257 and GMR-RA1153 obtained in Bogotá in the same hospital, shared the same rST-61,696 and the same clade C in the pangenome phylogenetic tree (Fig. 1), with the Colombian genome RB151 previously reported in Bucaramanga, Colombia, 2013 . This indicate the presence of rST-61,696 in two different geographic locations of the country and the possible wider circulation of this rST in Colombia. Likewise, the FDAARGOS_330 collected in 2015, in the Children’s National Hospital (USA), showed to be closely related to the Colombian genomes GMR-RA257, GMR-RA1153 and RB151 (Fig. 1) and was also typed as rST-61,696. This suggest a wider presence of P. rettgeri rST61696 in America.
The pangenome phylogenetic analysis suggests that GMR-RA257 and GMR-RA1153 were clonal considering that they were clustered very closely (distance < 0.01) in the tree (Fig. 1). It was also supported because they came from the same hospital, the ANI was 99.99% and the high identity between them showed by the BLAST comparative circular visualization (Fig. 2). Furthermore, both isolates carried the same resistance genes (Table 1). There were few differences between them and the Colombian P. rettgeri RB151 previously reported . These differences were principally in plasmid contigs since they were in the region (4.6–4.9 Mb) where the plasmid contig of the reference genome P. rettgeri RB151 was placed (Fig. 2).
As result of the SNP analysis we detected 108 variants of a single non-synonymous nucleotide (nsSNVs). Among those 107 nsSNVs were located in genes encoding hypothetical proteins and only one variant (Cys191Arg) of the GMR-RA1153 genome was located in the CpxA protein. CpxA is a transmembrane histidine kinase sensor of the two-component system Cpx-TCS [45, 46], involved in the stress response of the bacterial envelope by hostile factors such as physical stress (osmolarity), chemical stress (ethanol, pH, detergents) and misfolded proteins [45, 47]. Several studies have implicated the Cpx system in multidrug resistance (MDR), in particular to β-lactam antibiotics (imipenem, cefepime, ceftriaxone, ceftazidime, cefotaxime), chloramphenicol  and aminoglycosides . The mutation (Cys191Arg) in cpxA is located in the protein region 185–222, where have been reported mutations affecting its activity , which could explain the resistance profile differences: resistance and intermediate resistance to gentamicin and intermediate and resistance to ertapenem for GMR-RA1153 and GMR-RA257 respectively. The Kirby-Bauer disk diffusion susceptibility test results interpreted using the 2015 CLSI (Clinical and Laboratory Standards Institute) guidelines  showed that this mutation had no influence on resistance to amikacin because both isolates were resistant. Probably, some of the other SNPs could contribute to this variation, but more research will be necessary for understanding how the SNPs located in unknown proteins or in non-coding regions influence the resistance.
The isolates GMR-RA257 and GMR-RA1153 showed an uncommon high number of β-lactamase genes. The ESBL and AmpC β-lactamase genes blaCTX-M-15 and blaCMY-2 have not been reported previously in P. rettgeri isolates and only one study has described the presence of carbapenemase NDM-1 together with an ESBL as PER-1 . Furthermore, there is not a previously report of NDM-1 together with another ESBL such as CTX-M-15.
Apart from the high number of β-lactamase genes in both study isolates further resistance genes arnA, PmrE, PmrF, LpxA, LpxC, gyrB, folP, murA, rpoB, rpsL, and tet34 genes were found that have been reported previously in Enterobacterales species but no in P. rettgeri. Additionally, four efflux pumps of the MFS family: RosAB (rosA, rosB), EmrD, mdtH and cmlA, were not previously reported in P. rettgeri.
GMR-RA257 and GMR-RA1153 were resistant to the carbapenems imipenem and meropenem; this was highly probable caused by the production of MBL types NDM-1 and VIM-2 [51, 52] and the expression of blaOXA-10 encoding an enzyme that is currently associated with weak hydrolysis of carbapenem, and whose inclusion is proposed within the group of carbapenem-hydrolyzing class D β-lactamases (CHDLs) [53, 54]. In a study in Nigeria in 2011 a carbapenem susceptible P. rettgeri isolate was reported that carried blaOXA-10, blaVEB-1 and blaCMY-4 genes with no presence of the blaNDM-1 gene . These results show that the presence of MBL as blaNDM-1 and blaVIM-2 contribute to carbapenem resistance, which is a serious problem for clinicians, as these are considered last-resort antibiotics. Unfortunately, antimicrobials such as colistin and tigecycline –which are now used as antibiotics to treat multidrug-resistant microorganisms producing carbapenemase – are not a treatment option for P. rettgeri isolates, due to the natural resistance of the latter to the said antibiotics [56, 57]. However, a possible treatment option is a combination: Avibactam, non-β-lactam inhibitor with activity against serine β-lactamase classes A, C and D ; and aztreonam, resistant to the hydrolysis of MBL. In vitro studies demonstrated good results in the values of minimum inhibitory concentration (MIC) for Enterobacterales isolates that produced MBL and serine β-lactamases such as ESBL and / or AmpC [43, 59]. Aztreonam-avibactam is not yet approved by the U.S. Food and Drug Administration (FDA), but clinicians can administer this combination [43, 59, 60].
Whole genome analysis of the two clinical P. rettgeri isolates GMB-RA257 and GMB-RA1153 revealed that both isolates harbored an uncommon combination of various β-lactamase genes: blaNDM-1, blaVIM-2, blaCTX-M-15, blaCMY-2, blaOXA-10, and blaTEM-1, that have not been reported previously in this microorganism.
The comparative genomic analyses revealed a close relationship of the isolates; the assigned rST-61,696 was found previously in Colombia and the USA.
The differences of both isolates in gentamicin susceptibility could be explained by a Cys191Arg substitution in the CpxA protein.
Our finding shows the need of detailed and continuous molecular epidemiological surveillance of P. rettgeri to assess the regional and worldwide spread of clones and to detect early the acquisition of resistance genes that limit dramatically limit the treatment options for infections with this pathogen.
The Microbiology Group of the National Institute of Health (INS) of Colombia received the isolates GMR-RA257 and GMR-RA1153 that exhibited a decreased susceptibility to carbapenems, as part of the Healthcare-Associated Infection (HAI) Surveillance Program. The isolates were recovered from urine samples obtained from two patients admitted to a 4th level hospital in Bogotá – Colombia, in different months of 2015. GMR-RA1153 was isolated from a patient who was in an intensive care unit (ICU), and GMR-RA257 without data from the collection site. The two isolates were identified as Providencia rettgeri, using the automatic VITEK® 2 system (bioMérieux).
An amplification of the 16S rRNA gene (~ 1.5 kb) was made to confirm the genus and species of each isolate. For this, the universal primers Forward 27F 5 ‘AGAGTTTGATCMTGGCTCAG 3’ and Reverse 1492R 5 ‘TACGGYTACCTTGTTACGACTT 3’ were used. Purified PCR products were sequenced, using the BigDye methodology and following the manufacturer’s protocol (Applied Biosystems, Foster City, United States) at the sequencing service of the Institute of Genetics at Universidad Nacional de Colombia. Chromatograms were analyzed using Chromas program (Technelysium Pty Ltd.). The 16S rRNA sequencing data were aligned and compared with those available in the NCBI (National Center for Biotechnology Information) GenBank database using the BLASTn algorithm.
The antibiotic susceptibility profile was determined using the Kirby-Bauer disk diffusion susceptibility test. The following disks were tested: cefepime (FEP), cefoxitin (FOX), cefotaxime (CTX), ceftazidime (CAZ), ertapenem (ERT), imipenem (IMI), meropenem (MEM), aztreonam (ATM), ciprofloxacin (CIP), amikacin (AK), piperacillin/tazobactam (TZP) and trimethoprim/sulfamethoxazole (SXT). The susceptibility results were interpreted using the 2015 CLSI (Clinical and Laboratory Standards Institute) guidelines .
Confirmation of carbapenemase production
Phenotypic detection of carbapenemases was performed using modified Hodge test (MHT) , and two double-disc synergy tests with specific inhibitors: ethylenediaminetetraacetic acid / mercaptoacetic acid (EDTA / SMA) and phenyl boronic acid (APB). For the MHT, the reference strain Escherichia coli ATCC 25922 was used as the indicator organism, strain K. pneumoniae BAA 1705 as positive control and strain K. pneumoniae BAA 1706 as negative control. For the EDTA / SMA test the positive control used was K. pneumoniae BAA 2146 and for the APB test was the strain K. pneumoniae BAA 1705. The synergy of the antibiotic imipenem towards EDTA/SMA was interpreted as a positive result for the presence of class B carbapenemase. Similarly, the halo distortion of the antibiotic imipenem towards APB was interpreted as a positive result for the presence of class A carbapenemase.
Whole genome sequencing
The isolates’ genomic DNA was obtained by using two commercial kits, according to the manufacturer’s recommendations: QIAamp® DNA mini kit (Qiagen) for GMR-RA1153, and PureLink® Genomic DNA Kit (Thermo Fisher Scientific / Invitrogen) for GMR-RA257. DNA quantification was estimated by the PicoGreen method, using Vector 3 fluorometry. The final Paired-End sequencing was done with the Illumina HiSeq 2500 platform (the TruSeq Nano DNA kit (Illumina) was used for library construction). The quality of the raw paired-end reads was assessed by FastQC v0.11.7 , and the trimming of poor-quality bases (<Q30 and < 1000 bp) was done using the Sickle v1.33 tool . Subsequently, the reads were de novo assembled using SPAdes v3.8.0 . The annotation process was done using Prokka v1.13 , which was enriched with the following resistance databases: CARD , Resfam , GIPSy  and VFDB  and we complete the analysis with the Resfinder database (https://cge.cbs.dtu.dk/services/ResFinder/).
Providencia rettgeri genome sequence data
The genomes GMR-RA257 and GMR-RA1153 were compared with 16 genomes of Providencia rettgeri available in GenBank: 14 reference genomes (RB151, H1736, CCBH11880, FDAARGOS_330, MR4, 729–12, UBA5024, AR_0082, TUM9994, Dmel1, DSM1131, PR1, Pret_2032, PR002); and 2 collection strains: (NCTC11801, NCTC74801). An important factor to take into account is that these genomes were used without discrimination of chromosome and plasmid sequences.
Ribosomal multilocus sequence typing (rMLST)
The Ribosomal Multi-Locus Sequence Typing Scheme described by Jolley et al.  and hosted on the PubMLST website https://pubmlst.org/rmlst/ was used for the typing of the isolates of this work. Ribosomal Multilocus Sequence Typing (rMLST) is an approach that indexes the variation of the 53 genes encoding the bacterial ribosome protein subunits (rps, rpl, rpm genes) as a means of integrating microbial typing. The allelic variation rMLST is cataloged using the BIGSdb platform : scalable, open source database that allows to store sequences and allelic definitions for defined loci  For the study, the rMLST scheme first identified the nucleotide sequences of the genes (rps, rpl, rpm) of the annotated genomes; then, defined the unique alleles for each of these genes based on their nucleotide sequence; and finally, assigned each genome a customized ribosomal sequence type (rST) number defined by the allelic profile of the numerical combination of the 53 genes.
Whole genome comparisons
Using the Providencia Rettgeri genome sequence data described above the pangenome was inferred with Roary version 3.7.0 . Roary produced a gene presence/absence matrix, a multi-FASTA alignment of core genes using PRANK version 0.140603  and a phylogenetic tree based on the presence and absence of accessory genes among taxa using FastTree version 2.1.9 .
Visualization of circular genome BLASTN-based comparisons between Colombian reference genome RB151 and the genomes GMR-RA257 and GMR-RA1153, in order to show a broad genome visualization of the identity of the coding, was done using CGView Comparison Tool .
The average nucleotide identity (ANI) approach was used  to determine the genetic relationship between GMR-RA257 and GMR-RA1153. ANI was calculated within the JSpecies software . GMR-RA257 and GMR-RA1153 were aligned with Mauve version 2.4.0 . From this alignment, a file with the SNPs was obtained using Mauve 2.4.0 and converted to VCF format using an in-house Perl script. Finally, the SNPs were annotated with the snpEff program (http://snpeff.sourceforge.net).
The detection of these resistance genes for Providencia Rettgeri genome sequence data was performed using the program Resfinder (https://cge.cbs.dtu.dk/services/ResFinder/), using the fasta files of the 18 genomes of P. rettgeri (GMR-RA257 and GMR-RA1153, 14 reference genomes and 2 collection strains).
Availability of data and materials
Sequence data of this project have been deposited in the DDBJ/EMBL/GenBank of the National Center for Biotechnology Information (NCBI) under the accession number VRPG00000000 and VRPH00000000.
Average Nucleotide Identity
Phenyl Boronic Acid
- EDTA / SMA:
Ethylenediaminetetraacetic acid / Mercaptoacetic acid
Intensive Care Unit
National Institute of Health
Modified Hodge test
Ribosomal Multilocus Sequence Typing
Ribosomal Sequence Type
O’Hara CM, Brenner FW, Miller JM. Classification, identification, and clinical significance of Proteus, Providencia, and Morganella. Clin Microbiol Rev. 2000;13(4):534–46.
Wie S-H. Clinical significance of Providencia bacteremia or bacteriuria. Korean J Intern Med. 2015;30(2):167–9.
Yoh M, Matsuyama J, Ohnishi M, Takagi K, Miyagi H, Mori K, et al. Importance of Providencia species as a major cause of travellers’ diarrhoea. J Med Microbiol. 2005;54(11):1077–82.
Koreishi AF, Schechter BA, Karp CL. Ocular infections caused by Providencia rettgeri. Ophthalmology. 2006;113(8):1463–6.
Washington MA, Barnhill J, Griffin JM. A case of wound infection with Providencia rettgeri and coincident gout in a patient from Guam. Hawaii J Med Public Health. 2015;74(11):375–7.
Maiti TK, Pandey P, Singh VK. Providencia rettgeri: an unusual cause of central nervous system infections. Am J Med Sci. 2013;346(2):158–9.
Lee G, Hong JH. Xanthogranulomatous pyelonephritis with nephrocutaneous fistula due to Providencia rettgeri infection. J Med Microbiol. 2011;60(7):1050–2.
Marquez-Ortiz RA, Haggerty L, Sim EM, Duarte C, Castro-Cardozo BE, Beltran M, et al. First Complete Providencia rettgeri Genome Sequence, the NDM-1-Producing Clinical Strain RB151. Genome Announc. 2017;5(3):e01472–16. https://doi.org/10.1128/genomeA.01472-16.
Olaitan AO, Diene SM, Assous MV, Rolain J-M. Genomic plasticity of multidrug-resistant NDM-1 positive clinical isolate of Providencia rettgeri. Genome Biol Evol. 2015;8(3):723–8.
Yong D, Toleman MA, Giske CG, Cho HS, Sundman K, Lee K, et al. Characterization of a new Metallo-β-lactamase gene, blaNDM-1, and a novel erythromycin esterase gene carried on a unique genetic structure in Klebsiella pneumoniae sequence type 14 from India. Antimicrob Agents Chemother. 2009;53(12):5046–54.
Tada T, Miyoshi-Akiyama T, Dahal RK, Sah MK, Ohara H, Shimada K, et al. NDM-1 Metallo-β-lactamase and ArmA 16S rRNA methylase producing Providencia rettgeri clinical isolates in Nepal. BMC Infect Dis. 2014;14:56.
Gefen-Halevi S, Hindiyeh MY, Ben-David D, Smollan G, Gal-Mor O, Azar R, et al. Isolation of genetically unrelated blaNDM-1 positive Providencia rettgeri strains in Israel. J Clin Microbiol. 2013;51(5):1642–3.
Lascols C, Hackel M, Marshall SH, Hujer AM, Bouchillon S, Badal R, et al. Increasing prevalence and dissemination of NDM-1 metallo-β-lactamase in India: data from the SMART study (2009). J Antimicrob Chemother. 2011;66(9):1992–7.
Shin S, Jeong SH, Lee H, Hong JS, Park M-J, Song W. Emergence of multidrug-resistant Providencia rettgeri isolates co-producing NDM-1 carbapenemase and PER-1 extended-spectrum β-lactamase causing a first outbreak in Korea. Ann Clin Microbiol Antimicrob. 2018;17(1):20.
Zhou G, Guo S, Luo Y, Ye L, Song Y, Sun G, et al. NDM-1–producing strains, family Enterobacteriaceae, in hospital, Beijing, China. Emerg Infect Dis. 2014;20(2):340–2.
Perry JD, Naqvi SH, Mirza IA, Alizai SA, Hussain A, Ghirardi S, et al. Prevalence of faecal carriage of Enterobacteriaceae with NDM-1 carbapenemase at military hospitals in Pakistan, and evaluation of two chromogenic media. J Antimicrob Chemother. 2011;66(10):2288–94.
Carvalho-Assef APD, Pereira PS, Albano RM, Berião GC, Chagas TPG, Timm LN, et al. Isolation of NDM-producing Providencia rettgeri in Brazil. J Antimicrob Chemother. 2013;68(12):2956–7.
do Carmo Junior NV, Filho HF, Costa DA GE, AJW C, de O GD, JJD F. First report of a NDM-producing Providencia rettgeri strain in the state of São Paulo. Braz J Infect Dis Off Publ Braz Soc Infect Dis. 2015;19(6):675–6.
Pasteran F, Meo A, Gomez S, Derdoy L, Albronoz E, Faccone D, et al. Emergence of genetically related NDM-1-producing Providencia rettgeri strains in Argentina. J Glob Antimicrob Resist. 2014;2(4):344–5.
Saavedra-Rojas S-Y, Duarte-Valderrama C, González-de-Arias M-N, Ovalle-Guerro MV. Emergencia de Providencia rettgeri NDM-1 en dos departamentos de Colombia, 2012-2013. Enfermedades Infecc Microbiol Clínica. 2017;35(6):354–8.
Zurita J, Parra H, Gestal MC, McDermott J, Barba P. First case of NDM-1-producing Providencia rettgeri in Ecuador. J Glob Antimicrob Resist. 2015;3(4):302–3.
Kus JV, Tadros M, Simor A, Low DE, McGeer AJ, Willey BM, et al. New Delhi metallo-β-lactamase-1: local acquisition in Ontario, Canada, and challenges in detection. CMAJ Can Med Assoc J. 2011;183(11):1257–61.
Barrios H, Garza-Ramos U, Reyna-Flores F, Sanchez-Perez A, Rojas-Moreno T, Garza-Gonzalez E, et al. Isolation of carbapenem-resistant NDM-1-positive Providencia rettgeri in Mexico. J Antimicrob Chemother. 2013;68(8):1934–6.
Tshisevhe VS, Lekalakala MR, Tshuma N, Janse van Rensburg S, Mbelle N. Outbreak of carbapenem-resistant Providencia rettgeri in a tertiary hospital. South Afr Med J Suid-Afr Tydskr Vir Geneeskd. 2016;107(1):31–3.
Pfeifer Y, Trifonova A, Pietsch M, Brunner M, Todorova I, Gergova I, et al. Clonal transmission of gram-negative Bacteria with Carbapenemases NDM-1, VIM-1, and OXA-23/72 in a Bulgarian hospital. Microb Drug Resist Larchmt N. 2017;23(3):301–7.
Ovalle MV, Rojas SYS, Valderrama CD, Durán MA. Title of subordinate document. In: Resultados del Programa de Vigilancia por Laboratorio de Resistencia antimicrobiana en Infecciones Asociadas a la Atención en Salud (IAAS) 2016; 2016. https://www.ins.gov.co/buscador-eventos/Informacin%20de%20laboratorio/Informe%20Vigilancia%20por%20Laboratorio%20Resistencia%20Antimicrobiana%20y%20Whonet%20IAAS%202016.pdf. Accessed 10 October 2019.
Oteo J, Miró E, Pérez-Vázquez M, Navarro F. Evolution of carbapenemase-producing Enterobacteriaceae at the global and national level: what should be expected in the future? Enferm Infecc Microbiol Clin. 2014;32(Suppl 4):17–23.
Blakiston M, Heffernan H, Roberts S, Freeman J. The clear and present danger of carbapenemase-producing Enterobacteriaceae (CPE) in New Zealand: time for a national response plan. N Z Med J. 2017;130(1454):72–9.
Dhillon RH-P, Clark J. ESBLs: a clear and present danger? Crit Care Res Pract. 2012;2012:625170.
Otlu B, Yakupoğulları Y, Gürsoy NC, Duman Y, Bayındır Y, Tekerekoğlu MS, et al. Co-production of OXA-48 and NDM-1 Carbapenemases in Providencia rettgeri: the first report. Mikrobiyol Bul. 2018;52(3):300–7.
Lee H-W, Kang H-Y, Shin K-S, Kim J. Multidrug-resistant Providencia isolates carrying blaPERR-1, blaVIM-2, and armA. J Microbiol Seoul Korea. 2007;45(3):272–4.
Aibinu IE, Pfeifer Y, Ogunsola F, Odugbemi T, Koenig W, Ghebremedhin B. Emergence of β-lactamases OXA-10, VEB-1 and CMY in Providencia spp. from Nigeria. J Antimicrob Chemother. 2011;66(8):1931–2.
Magiorakos A-P, Srinivasan A, Carey RB, Carmeli Y, Falagas ME, Giske CG, et al. Multidrug-resistant, extensively drug-resistant and pandrug-resistant bacteria: an international expert proposal for interim standard definitions for acquired resistance. Clin Microbiol Infect Off Publ Eur Soc Clin Microbiol Infect Dis. 2012;18(3):268–81.
Galac MR, Lazzaro BP. Comparative genomics of bacteria in the genus Providencia isolated from wild Drosophila melanogaster. BMC Genomics. 2012;13:612.
Pereira PS, Albano RM, Asensi MD, Carvalho-Assef APD. Draft genome sequences of three NDM-1-producing Enterobacteriaceae species isolated from Brazil. Mem Inst Oswaldo Cruz. 2015;110(4):580–2.
Wang G-H, Brucker RM. Genome sequence of Providencia rettgeri NVIT03, isolated from Nasonia vitripennis. Microbiol Resour Announc. 2019;8(3):e01157–18.
Nawaz M, Sung K, Kweon O, Khan S, Nawaz S, Steele R. Characterisation of novel mutations involved in quinolone resistance in Escherichia coli isolated from imported shrimp. Int J Antimicrob Agents. 2015;45(5):471–6.
O’Regan E, Quinn T, Pagès J-M, McCusker M, Piddock L, Fanning S. Multiple regulatory pathways associated with high-level ciprofloxacin and multidrug resistance in salmonella enterica Serovar Enteritidis: involvement of ramA and other global regulators. Antimicrob Agents Chemother. 2009;53(3):1080–7.
Bankevich A, Nurk S, Antipov D, Gurevich AA, Dvorkin M, Kulikov AS, et al. SPAdes: a new genome assembly algorithm and its applications to single-cell sequencing. J Comput Biol J Comput Mol Cell Biol. 2012;19(5):455–77.
Krawczyk PS, Lipinski L, Dziembowski A. PlasFlow: predicting plasmid sequences in metagenomic data using genome signatures. Nucleic Acids Res. 2018;46(6):e35.
Pecora N, Zhao X, Nudel K, Hoffmann M, Li N, Onderdonk A, Yokoe D, Brown E, Allard M, Bry L. Diverse vectors and mechanisms spread New Delhi Metallo-β-lactamases among Carbapenem-resistant Enterobacteriaceae in the greater Boston area. Antimicrob Agents Chemother. 2019;63(2):e02040–18.
Partridge S, Kwong S, Firth N, Jensenc S. Mobile genetic elements associated with antimicrobial resistance. Clin Microbiol Rev. 2018;31:e00088–17.
Wu W, Feng Y, Tang G, Qiao F, McNally A, Zong Z. NDM Metallo-β-lactamases and their bacterial producers in health care settings. Clin Microbiol Rev. 2019;32(2):e00115–8.
Jolley KA, Maiden MCJ. BIGSdb: scalable analysis of bacterial genome variation at the population level. BMC Bioinformatics. 2010;11:595.
Srinivasan VB, Vaidyanathan V, Mondal A, Rajamohan G. Role of the Two Component Signal Transduction System CpxAR in Conferring Cefepime and Chloramphenicol Resistance in Klebsiella pneumoniae NTUH-K2044. PLoS ONE. 2012;7(4):e33777.
DiGiuseppe PA, Silhavy TJ. Signal detection and target gene induction by the CpxRA two-component system. J Bacteriol. 2003;185(8):2432–40.
MacRitchie DM, Buelow DR, Price NL, Raivio TL. Two-component signaling and gram negative envelope stress response systems. Adv Exp Med Biol. 2008;631:80–110.
Hirakawa H, Nishino K, Hirata T, Yamaguchi A. Comprehensive studies of drug resistance mediated by overexpression of response regulators of two-component signal transduction Systems in Escherichia coli. J Bacteriol. 2003;185(6):1851–6.
Mechaly AE, Sassoon N, Betton JM, Alzari PM. Segmental helical motions and dynamical asymmetry modulate histidine kinase autophosphorylation. PLoS Biol. 2014;12(1):e1001776.
Clinical and Laboratory Standards Institute. 2015. Performance standards for antimicrobial susceptibility testing, M100-S25. Wayne, PA.
Poirel L, Naas T, Nordmann P. Diversity, epidemiology, and genetics of class D β-lactamases. Antimicrob Agents Chemother. 2010;54(1):24–38.
Antunes NT, Lamoureaux TL, Toth M, Stewart NK, Frase H, Vakulenko SB. Class D β-lactamases: are they all carbapenemases? Antimicrob Agents Chemother. 2014;58(4):2119–25.
Ye Y, Xu L, Han Y, Chen Z, Liu C, Ming L. Mechanism for carbapenem resistance of clinical Enterobacteriaceae isolates. Exp Ther Med. 2018;15(1):1143–9.
Khan AU, Maryam L, Zarrilli R. Structure, genetics and worldwide spread of New Delhi Metallo-β-lactamase (NDM): a threat to public health. BMC Microbiol. 2017;17(1):101.
Stock I. Infectious diseases caused by carbapenemase-producing Enterobacteriaceae--a particular challenge for antibacterial therapy. Med Monatsschr Pharm. 2014;37(5):162–72 quiz 173-4.
Rose WE, Rybak MJ. Tigecycline: first of a new class of antimicrobial agents. Pharmacotherapy. 2006;26(8):1099–110.
Choi HK, Kim YK, Kim HY, Park JE, Uh Y. Clinical and microbiological features of Providencia bacteremia: experience at a tertiary care hospital. Korean J Intern Med. 2015;30(2):219–25.
Bush K, Bradford PA. Interplay between β-lactamases and new β-lactamase inhibitors. Nat Rev Microbiol. 2019;17(5):295.
Karlowsky JA, Kazmierczak KM, de Jonge BLM, Hackel MA, Sahm DF, Bradford PA. In Vitro Activity of Aztreonam-Avibactam against Enterobacteriaceae and Pseudomonas aeruginosa Isolated by Clinical Laboratories in 40 Countries from 2012 to 2015. Antimicrob Agents Chemother. 2017;61(9):e00472–17. https://doi.org/10.1128/AAC.00472-17.
Emeraud C, Escaut L, Boucly A, Fortineau N, Bonnin RA, Naas T, et al. Aztreonam plus Clavulanate, Tazobactam, or Avibactam for treatment of infections caused by Metallo-β-lactamase-producing gram-negative Bacteria. Antimicrob Agents Chemother. 2019;63(5):e00010–9.
Andrews S. Title of subordinate document. In: FastQC A Quality Control tool for High Throughput Sequence Data; 2013. http://www.bioinformatics.babraham.ac.uk/projects/fastqc/. Accessed 10 October 2019.
Joshi N, Fass J. Title of subordinate document. In: Sickle: A sliding-window, adaptive, quality-based trimming tool for FastQ files; 2011. http://github.com/najoshi/sickle. Accessed 10 October 2019.
Seemann T. Prokka: rapid prokaryotic genome annotation. Bioinforma Oxf Engl. 2014;30(14):2068–9.
McArthur AG, Waglechner N, Nizam F, Yan A, Azad MA, Baylay AJ, et al. The comprehensive antibiotic resistance database. Antimicrob Agents Chemother. 2013;57(7):3348–57.
Gibson MK, Forsberg KJ, Dantas G. Improved annotation of antibiotic resistance determinants reveals microbial resistomes cluster by ecology. ISME J. 2015;9(1):207–16.
Soares SC, Geyik H, Ramos RTJ, de Sá PHCG, Barbosa EGV, Baumbach J, et al. GIPSy: genomic island prediction software. J Biotechnol. 2016;232:2–11.
Chen L, Zheng D, Liu B, Yang J, Jin Q. VFDB 2016: hierarchical and refined dataset for big data analysis—10 years on. Nucleic Acids Res. 2016;44(Database issue):D694–7.
Jolley KA, Bliss CM, Bennett JS, Bratcher HB, Brehony C, Colles FM, et al. Ribosomal multilocus sequence typing: universal characterization of bacteria from domain to strain. Microbiol Read Engl. 2012;158(Pt 4):1005–15.
Jolley KA, Bray JE, Maiden MCJ. Open-access bacterial population genomics: BIGSdb software, the PubMLST.org website and their applications. Wellcome Open Res. 2018. https://doi.org/10.12688/wellcomeopenres.14826.1.
Page AJ, Cummins CA, Hunt M, Wong VK, Reuter S, Holden MTG, et al. Roary: rapid large-scale prokaryote pan genome analysis. Bioinforma Oxf Engl. 2015;31(22):3691–3.
Löytynoja A. Phylogeny-aware alignment with PRANK. Methods Mol Biol Clifton NJ. 2014;1079:155–70.
Price MN, Dehal PS, Arkin AP. FastTree: computing large minimum evolution trees with profiles instead of a distance matrix. Mol Biol Evol. 2009;26(7):1641–50.
Grant JR, Arantes AS, Stothard P. Comparing thousands of circular genomes using the CGView comparison tool. BMC Genomics. 2012;13:202.
Figueras MJ, Beaz-Hidalgo R, Hossain MJ, Liles MR. Taxonomic affiliation of new genomes should be verified using average nucleotide identity and multilocus phylogenetic analysis. Genome Announc. 2014;2(6):e00927–14.
Richter M, Rosselló-Móra R. Shifting the genomic gold standard for the prokaryotic species definition. Proc Natl Acad Sci U S A. 2009;106(45):19126–31.
Darling ACE, Mau B, Blattner FR, Perna NT. Mauve: multiple alignment of conserved genomic sequence with rearrangements. Genome Res. 2004;14(7):1394–403.
We would like to thank the Microbiology group of the National Institute of Health of Colombia, the Molecular Epidemiology Laboratory and the Bioinformatics Group of the Universidad Nacional de Colombia and the Research Division of the Universidad Nacional de Colombia at Bogotá-Colombia
The Universidad Nacional de Colombia at Bogotá (Project 35030) financially supported this work and the design of the study, analysis, and interpretation of data, and manuscript writing. The Instituto Nacional de Salud de Colombia supported the collection, analysis, and interpretation of data. The Departamento Administrativo de Ciencia y Tecnología, COLCIENCIAS under the Convocatoria 656, 2014, Es Tiempo de Volver (FP44842–155-2015) financially supported the VR postdoctoral fellow.
Ethics approval and consent to participate
The isolates GMR-RA257 and GMR-RA1153 analyzed in this manuscript were recovered as part of the national surveillance program by the antimicrobial resistance Healthcare-associated infections (HAIs) laboratory of the National Institute of Health of Colombia. Due that the isolates are part of a national public health surveillance program, they do not require informed consent (decrees 2323 and 3518 of 2006 of the Social Protection Ministry of Colombia).
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.
About this article
Cite this article
Piza-Buitrago, A., Rincón, V., Donato, J. et al. Genome-based characterization of two Colombian clinical Providencia rettgeri isolates co-harboring NDM-1, VIM-2, and other β-lactamases. BMC Microbiol 20, 345 (2020). https://doi.org/10.1186/s12866-020-02030-z
- Whole-genome sequencing
- bla CTX-M-15
- bla CMY-2
- bla OXA-10
- And blaTEM-1