Next Article in Journal
Antimicrobial Effects of Minocycline, Tigecycline, Ciprofloxacin, and Levofloxacin against Elizabethkingia anophelis Using In Vitro Time-Kill Assays and In Vivo Zebrafish Animal Models
Next Article in Special Issue
Molecular Epidemiology of Carbapenem-Resistant Acinetobacter baumannii Isolates from Northern Africa and the Middle East
Previous Article in Journal
Rapid Molecular Tests for Detecting Respiratory Pathogens Reduced the Use of Antibiotics in Children
Previous Article in Special Issue
Whole-Genome Sequencing Reveals the Presence of the blaCTX-M-65 Gene in Extended-Spectrum β-Lactamase-Producing and Multi-Drug-Resistant Clones of Salmonella Serovar Infantis Isolated from Broiler Chicken Environments in the Galapagos Islands
 
 
Search for Articles:
Title / Keyword
Author / Affiliation / Email
Journal
Article Type
 
 
Advanced Search
 
Section
Special Issue
Volume
Issue
Number
Page
 
Journals
Antibiotics
Volume 10
Issue 3
10.3390/antibiotics10030284
antibiotics-logo
Submit to this Journal Review for this Journal Propose a Special Issue
► Article Menu

Article Menu

share Share announcement Help format_quote Cite question_answer Discuss in SciProfiles thumb_up
...
Endorse textsms
...
Comment
first_page
settings
Order Article Reprints
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Molecular Analysis of Polymyxin Resistance among Carbapenemase-Producing Klebsiella pneumoniae in Colombia

by
Elsa De La Cadena
1,*,
María Fernanda Mojica
1,2,3,
Juan Carlos García-Betancur
1,
Tobías Manuel Appel
1,
Jessica Porras
1,
Christian José Pallares
1,4,
Juan Sebastián Solano-Gutiérrez
5,
Laura J. Rojas
2,3 and
María Virginia Villegas
1,4
1
Grupo de Investigación en Resistencia Antimicrobiana y Epidemiologia Hospitalaria, Universidad El Bosque, Bogotá 110121, Colombia
2
Department of Infectious Diseases, Case Western Reserve University, Cleveland, OH 44106-7164, USA
3
Research Service, Louis Stokes Veterans Affairs Medical Center, Cleveland, OH 44106-7164, USA
4
Comité de Infecciones y Vigilancia Epidemiológica, Clínica Imbanaco, Cali 760031, Colombia
5
Departamento de Ciencias Biológicas, Escuela de Ciencias, Universidad EAFIT, Medellín 050022, Colombia
*
Author to whom correspondence should be addressed.
Antibiotics 2021, 10(3), 284; https://doi.org/10.3390/antibiotics10030284
Submission received: 31 January 2021 / Revised: 2 March 2021 / Accepted: 4 March 2021 / Published: 10 March 2021
(This article belongs to the Special Issue Global Antimicrobial Resistance: What Are the Current Emerging Threats?)
Browse Figures
Review Reports Versions Notes

Abstract

:
Polymyxin resistance in Klebsiella pneumoniae has been attributed to mutations in mgrB, phoPQ, pmrAB, and crrAB and to the presence of mcr plasmid-mediated genes. Herein, we describe the molecular characteristics of 24 polymyxin- and carbapenem-resistant K. pneumoniae isolates recovered from six Colombian cities between 2009 and 2019. Minimum inhibitory concentrations (MICs) to polymyxin were confirmed by broth microdilution, and whole-genome sequencing was performed to determine sequence type, resistome, and mutations in the genes related to polymyxin resistance, as well the presence of mcr. The results showed high-level resistance to polymyxin (MICs ≥ 4 μg/mL). blaKPC-3 was present in the majority of isolates (17/24; 71%), followed by blaKPC-2 (6/24; 25%) and blaNDM-1 (1/24; 4%). Most isolates belonged to the CG258 (17/24; 71%) and presented amino acid substitutions in PmrB (22/24; 92%) and CrrB (15/24; 63%); mutations in mgrB occurred in only five isolates (21%). Additional mutations in pmrA, crrA, and phoPQ nor any of the mcr resistance genes were identified. In conclusion, we found clonal dissemination of polymyxin and carbapenem-resistant K. pneumoniae isolates in Colombia, mainly associated with CG258 and blaKPC-3. Surveillance of this multidrug-resistant clone is warranted due to the limited therapeutic options for the treatment of carbapenem-resistant K. pneumoniae infections.

1. Introduction

Klebsiella pneumoniae is a notorious nosocomial pathogen given its worldwide distribution and its capacity to develop resistance to multiple classes of antibiotics, including aminoglycosides, fluoroquinolones, extended-spectrum cephalosporins, and carbapenems [1]. Carbapenems are the preferred treatment for infections caused by K. pneumoniae isolates harboring extended-spectrum β-lactamases (ESBLs) [2,3]. Infections caused by carbapenem-resistant K. pneumoniae are difficult to treat and associated with high mortality rates, on average 41% [4]. Resistance to carbapenems has increased in K. pneumoniae, mainly due to the production of carbapenemases belonging to Ambler class A (KPC), class B (NDM and VIM), and class D (OXA-48). Polymyxins have re-emerged in clinical practice due to the limited antibiotic development pipeline and worldwide increasing prevalence of nosocomial infections caused by multidrug-resistant (MDR) Gram-negative bacteria including carbapenemase-producing K. pneumoniae. Polymyxin B and colistin (polymyxin E) have been ultimately considered as the last-resort treatment of such infections [5].
Polymyxins are cyclic polypeptides, with primary amines carrying a net polycationic charge at physiological pH and a hydrophobic fatty-acyl side chain. The amphipathic nature of polymyxins enables them to form ionic and hydrophobic molecular interactions at the bacterial outer membrane. Specifically, polymyxins initially bind to the lipid A of the lipopolysaccharide (LPS) through an ionic interaction that leads to the disruption of the outer membrane and the hydrophobic insertion of the fatty acyl side chain of polymyxin into lipid A. Subsequently, cytoplasmic membrane disruption causing cell death. Naturally, the increased use of polymyxins has led to the appearance of diverse mechanisms of polymyxin resistance, particularly during the course of therapy [6].
Resistance to polymyxin is provided by the addition of cationic groups to lipid A, either via chromosomally encoded pathways or plasmid-borne mobile colistin resistance gene (mcr) [7]. In K. pneumoniae, modification in several two-component systems such as PhoPQ, PmrAB, and CrrAB, as well as partial deletions or inactivation by different insertion sequence (IS) of the mgrB gene are currently the most commonly reported mechanisms of resistance [8]. In addition to these mechanisms, Jana et al. identified 35 nonessential chromosomal genes involved in polymyxin resistance among K. pneumoniae ST258 strains. Of those genes, collectively called “secondary resistome”, dedA encoding a putative integral membrane was shown to have a prominent role, as its deletion completely restored colistin susceptibility [9].
In Colombia, several molecular epidemiological studies have established that K. pneumoniae-producing KPC is endemic and represents a significant proportion of all clinical K. pneumoniae isolates in the country [10]. Furthermore, half of them belong to the clonal group (CG) 258 [11], which has been associated with the emergence of polymyxin resistance in K. pneumoniae in several countries [2,11,12]. Nevertheless, data of the prevalence of polymyxin-resistant K. pneumoniae in the country are scarce, mainly due to the shortcomings of the routine laboratory testing for polymyxin susceptibility [2,13]. Therefore, in this study we aimed to determine the profile and mechanism of polymyxin resistance in carbapenem-resistant K. pneumoniae isolates from six Colombian cities over a period of 10 years. Hereby, we describe the clonal dissemination of polymyxin and carbapenem-resistant K. pneumoniae isolates in several Colombian cities, mainly due to simultaneous mutations in mgrB, pmrB, and CrrB associated with the high-risk clone ST258.

2. Results

2.1. Antimicrobial Susceptibility

The susceptibility profile of the 24 isolates of K. pneumoniae included in this study is shown in Table 1. According to the standard definitions for acquired resistance proposed by Magiorakos et al. [14], all isolates exhibited multidrug-resistant phenotypes and in particular, were resistant to polymyxin with MICs ranging between 4 to ≥8 mg/L. The results also revealed high-level resistance to ceftazidime, cefepime, piperacillin-tazobactam, ertapenem, and ciprofloxacin. Ceftazidime-avibactam had the highest susceptibility percentage; as expected, only isolates producing NDM were resistant to it.

2.2. Molecular Typing and Virulome Analysis

According to the MLST analysis shown in Figure 1, the 24 isolates belonged to five different sequence types (STs). Most isolates belonged to CG258 (71%) (ST258, 16/24 and ST512, 1/24), followed by the ST14 (4/24, 17%), ST129 (2/24, 8%), and ST219 with one isolate (4%). Most of the KPC-3-K. pneumoniae producers belonged to the CG258 (16/17; 94%), while the isolate producing NDM-1 belonged to ST129 (Figure 1). Phylogenomic analysis performed using the CSI phylogeny tool grouped the different polymyxin-resistant K. pneumoniae isolates in two different clusters, coinciding with the presence of the blaKPC-3 and blaKPC-2 and in agreement with the MLST analysis, which clustered all CG258 isolates together (Figure 1). Regarding to virulence genes, KPC-3-producing K. pneumoniae belonging to ST258 and ST512 harbored multiple variants of distinct type 3 fimbriae genes, the yersiniabactin siderophore (ybt), and the genotoxic colibactin (clb), which have been associated with the type of capsular 154 (Supplementary Figure S1). In addition, the analysis of plasmid types revealed sequences associated with the incompatibility groups colRNAI, IncFIB (K), IncFII (K), IncFII (YP), IncI2, and IncR (Supplementary Figure S1).

2.3. Resistome Analysis

A wide variety of antibiotic resistance determinants to aminoglycosides, β-lactams, fluoroquinolones, macrolides, chloramphenicol, sulfonamide, fosfomycin, phenicol, tetracycline, co-trimoxazole, and rifampin were found within the genome of the isolates studied (Figure 1). Regarding quinolone resistance, analysis revealed the presence of several plasmid-mediated quinolone resistance (PMQR) genes such as qnrB, qnrS, and aac(6′)lb-cr. Additionally, mutations were found in the quinolone resistance-determining region (QRDR) of gyrA, leading to amino acid changes in Ser83Ile and Asp87Asn as well as in parC, leading to the amino acid change in Ser80Ile. The presence of oqxA and fosA was found in all isolates. Most of the isolates harbored blaKPC, being blaKPC-3 the most prevalent as it was detected in 17 of 24 isolates (71%); blaKPC-2 was carried by six isolates (25%), and only one isolate harbored blaNDM-1. On the other hand, analysis of the genetic context of the blaKPC revealed that the majority of isolates carried this gene within a Tn4401b structure (87%), followed by Tn4401a (8.7%), and not-Tn4401 (4.3%) (Figure 1).

2.4. Mechanisms of Resistance towards Polymyxins

Amino acid sequences from other chromosomal genes found in these K. pneumoniae isolates associated with resistance to polymyxin were compared to those available in the GenBank, and the changes found are listed in Figure 2. Mutations in phoPQ and crrA genes were not detected. We found changes in mgrB gene in five isolates (21%), including the presence of an IS-interrupting the mgrB gene in three of the 24 isolates (12.5%). Insertions of the IS element within the gene were found at nucleotides 74 and 99. Moreover, in two isolates (8%), mgrB was prematurely terminated due to a single-nucleotide change at position 88 causing the introduction of a stop codon, which generates a shorter protein of only 29 amino acids (Figure 2).
On the other hand, sequence analysis of the pmrB genes revealed several non-synonymous mutations that yield the following amino acid substitutions in 22 polyxymin resistant isolates: T157P, N105K, R256G, and T246A (Figure 2). According to the prediction of the SMART tool (http://smart.embl-heidelberg.de/) (accessed on 28 July 2020), T157 is located in the histidine kinase A (HisKA) phosphoacceptor domain, N105 and R256 in the HATPase_c domain, and T246 does not appear to fall within any conserved domain. In contrast, one-point mutation in the pmrA was found in only one isolate. In regard to CrrB, we identified three amino acid substitutions in this protein (C68S, G183V, and Q296L) among 15 of the 24 polymyxin-resistant isolates examined. Structurally, C68 is in the transmembrane region, G183 in the histidine kinase A (HisKA) phosphoacceptor domain and Q296 in the HATPase_c domain. The crrB gene was absent from nine isolates (Figure 2). Interestingly, mcr genes were not detected in any of the isolates.

3. Discussion

Several studies have shown that mutations identified within the mgrB gene and its insertional inactivation mediated by mobile IS are the most prevalent mechanisms involved in the development of polymyxin resistance in K. pneumoniae [15]. Our results confirmed these findings and showed that polymyxin resistance was associated to loss of function in mgrB and mutations in the pmrB and crrB genes. Major loss of function mutations in mgrB (premature stop codon, insertional inactivation, or deletion) were detected in five (21%) of the polymyxin-resistant strains (Figure 2). In three isolates (12.5%), the inactivation of mgrB gene was identified by insertion of IsKpn26 and two with premature stop codons. Jaidane et al. showed that colistin resistance in K. pneumoniae was directlylinked to mutations identified within the MgrB protein [16]; however, we did not find any amino acid substitutions in MgrB. Other studies have shown that the insertion of IS5 element into the mgrB gene was involved in resistance to polymyxins [15,17]. A large number of different ISs have been associated with transposition into mgrB including IS5, ISL3, IS5-like, IsKpn13, ISKpn14, and ISKpn26 with the IS5 family being the most commonly reported [8,18,19]. Surprisingly, in the isolates analyzed in this study, an IS5-like sequence was inserted at the same location in two of the isolates (UEB_12 and UEB_13). Insertional inactivation of mgrB occurred at nucleotide 74 and was caused by an IS that shares identity with ISKpn26 and belongs to the IS5 family [12,15]. The repeated observation of this insertional inactivation could indicate that is a hot spot for insertions. Likewise, Novovic et al. described one isolate with premature stop codon (TAG) due to a C-to-T change at position 88, which generates a truncated MgrB protein of 29 amino [20], just like the one found in the isolates UEB_14 and UEB_18. Truncations identified at positions 26, 28, and 30 of mgrB have also been previously described [21,22].
Sequence analysis of pmrB and crrB genes known to be involved in LPS synthesis revealed several polymorphisms in nucleotide sequences compared to the wild-type sequence of K. pneumoniae [18,23]. Notably, these polymorphisms were linked to specific STs. Within strains belonging to ST258, amino acid changes were observed in T157P, T246A, and R256G in PmrB, and C68S and Q296L in CrrB. Most belonging to the CG258 carried the double PmrB variant T246A R256G; in addition to these, the isolate identified as UEB_2 had also the T157P substitution. Interestingly, Jaidane et al. reported the T246A PmrB variant in K. pneumoniae ST111 isolates and the T246A R256G double PmrB variant in isolates belonging to other STs (392 and 247) collected in Tunisia [16]. Likewise, Jayol et al. described the PmrB variant T157P, which contributes to polymyxin resistance, in isolates recovered from different geographic regions and belonging to different STs. However, the Colombian isolate reported in that study carrying that mutation belonged to the CG258 [24]. On the other hand, the N105K substitution in PmrB was exclusively found among isolates belonging to the ST14. A different substitution in the same position (N105S) had been previously identified in K. pneumoniae isolates from Taiwan and South Korea [25], although molecular evidence about its role in conferring polymyxin resistance is not yet available. On the other hand, mutations in the crrB gene have also been associated with polymyxin resistance. Specifically, those leading to the Q10L, L94M, and N141Y amino acid substitutions have been shown to confer high-level resistance to polymyxin [18,26]. None of these crrB mutations were found in our analyzed isolates, and more studies are needed to elucidate the role of the C68S, G183V, and Q296L substitutions in the function of CrrB and their association with polymyxin resistance.
Different studies carried out in France with Colombian isolates mostly including ST258 isolates carrying blaKPC-2 and blaKPC-3 found mutations in CrrB (G183V), particularly in isolates producing KPC-3 [18]. Additionally, five isolates from the CG258 producing KPC-2 showed the insertion of different IS within the mgrB gene (IS5-like, IsKpn13, IsKpn14, Is10R) [12], and three isolates, two of them belonging to ST258 and producing KPC-3, presented a mutation in pmrB (T157P) [24]. In line with the previously reported molecular epidemiology of KPC-producing K. pneumoniae isolates in Colombia [10,27,28], 17 of the 24 isolates (71%) belonged to the epidemic CG258, and 16/17 (94%) harbored the blaKPC-3. As shown in Figure 1, WGS phylogenomic analysis effectively clustered the polymyxin-resistant K. pneumoniae isolates in agreement with both the presence of the blaKPC-3 and blaKPC-2 and the clustering obtained by MLST analysis, demonstrating a clonal pattern. Consistently with the clonal spread reported in this study, large outbreaks of polymyxin-resistant, KPC-3-harboring K. pneumoniae from ST258 and ST512 (SLV of ST258) have been reported in Israel, Italy, Greece, and the United States [10,29,30].
Interestingly, isolates found within the two main clusters have diverse geographical origins and additionally were isolated in a broad time range (from 2009 to 2019). Nevertheless, the close genetic relatedness of isolates collected from the same geographical region that belong to the same ST group and carry the same blaKPC gene, suggests that they also tend to share a particular polymyxin-resistance profile. Indeed, the inactivation of mgrB by ISKpn26 was only observed in isolates from one hospital in Medellin (M2) belonging to the ST258 (Figure 2). Cienfuegos-Gallet et al. reported similar results in K. pneumoniae isolates from this Colombian city [31]. In that study, they analyzed 32 carbapenem and colistin-resistant K. pneumoniae isolates, of which 27 belonged to ST512. They reported the presence of ISKpn25 within the mgrB gene in a cluster of 19 isolates from the same clonal group isolated from the same hospital. Likewise, additional mutations in phoPQ, pmrAB, and crrAB were not identified [31].
Of note, in accordance with the secondary resistome associated with colistin resistance in K. pneumoniae reported by Jana et al., [9] we found that all our polymyxin resistant isolates carried wild type dedA. Lastly, mcr genes were not detected suggesting that although this gene has been found in clinical Enterobacterales isolates, including K. pneumoniae [32], its prevalence among this species remains low. Therefore, resistance to polymyxin in isolates of K. pneumoniae in Colombia appears to be mainly restricted to chromosomal mutations, as previously described by Berglund et al. [33].

4. Materials and Methods

4.1. Bacterial Isolates and Antimicrobial Susceptibility

A total of 24 nonduplicate polymyxin and carbapenem resistant K. pneumoniae isolates from nine hospitals from six different Colombian cities were used in this study. These clinical isolates were collected through the Colombian Bacterial Resistance Surveillance Network between 2009 and 2019. Species identification was performed using the VITEK2 system (bioMerieux, Marcy l’Etoile, France). The minimum inhibitory concentration (MIC) of ceftazidime, cefepime, piperacillin-tazobactam, ertapenem, imipenem, meropenem, doripenem, ceftazidime-avibactam, ciprofloxacin, and tigecycline was determined by broth microdilution assay using customized Sensititre plates (Sensititre™; Thermo Fisher, Waltham, MA, USA) and interpreted according to the Clinical and Laboratory Standards Institute guidelines (CLSI) [34]. Resistance was determined according to CLSI guidelines except for tigecycline for which the European Committee on Antimicrobial Susceptibility Testing (EUCAST) was used (http://www.eucast.org) (accessed on 16 February 2020). Polymyxin resistance was interpreted according to CLSI (MIC > 2 mg/L), based on the macrodilution method with polymyxin sulfate (Sigma-Aldrich, St. Louis, MO, USA) and cation-adjusted Mueller-Hinton broth according to CLSI M7 guidelines [35]. Escherichia coli ATCC 25922 (polymyxin susceptible) and an E. coli isolate 6770 harboring the mcr-1 gene previously characterized (GenBank: MVPG00000000.1) [32] and with an expected MIC > 4 g/mL) were used as control strains for antimicrobial susceptibility testing (Table 1).

4.2. Molecular Typing

The presence of carbapenemases was confirmed by qPCR as published by Correa et al. [36]. A total of 24 polymyxin-resistant and carbapenemase-producing K. pneumoniae isolates were whole-genome sequenced (WGS) using the Illumina Nextera XT library prep kit and the Illumina MiSeq sequencing platform (Illumina, San Diego, CA, USA). Sequence annotation was done in PATRIC, the bacterial bioinformatics database and analysis resource (http://www.patricbrc.org) (accessed on 20 February 2020) [37]. Genome assemblies were uploaded to the Center for Genomic Epidemiology (CGE; https://cge.cbs.dtu.dk) (accessed on 25 February 2020) to identify antibiotic-resistance genes, plasmidic profiles, and STs, using the ResFinder, Plasmid Finder database [38,39], and the MLST (multilocus sequence typing) 1.8 server, respectively (the MLST scheme included: rpoB, gapA, mdh, pgi, phoE, infB, and tonB) [40]. Virulence genes and capsular typing via the wzi gene were performed in silico using the typing scheme at BIGSdb (http://bigsdb.pasteur.fr/) (accessed on 9 March 2020) [41]. To identify the genetic environment of the blaKPC gene, sequences were mapped into the Tn4401 transposon sequence (GenBank accession number KT378597.1), which was used as reference. The consensus sequences from the mapping were compared with tnpR, tnpA, istB, istA, and blaKPC by BLAST analysis.
Insertion sequences (IS) and mutations in mgrB gene or in the two-component regulatory systems pmrAB, phoPQ, and crrAB were searched within these isolates. The insertional sequence was identified using the ISfinder database (http://www-is.biotoul.fr) (accessed on 12 May 2020) [42], and the resulting sequences were compared to those of the polymyxin-susceptible strain K. pneumoniae (GenBank accession number CP000647). Whole-genome, single-nucleotide-polymorphism (SNP) phylogeny among isolates was inferred using the CSI Phylogeny version 1.4 (https://cge.cbs.dtu.dk/services/CSIPhylogeny/) (accessed on 21 March 2020) [43] to determine the epidemiological relationship between the isolates. To infer phylogenies using this SNP alignment, we used all the default values suggested by the developers, including the FastTree method (PMID 19377059), which infers phylogenies from big sequence-datasets based on the ‘‘minimum-evolution’’ principle instead of distance matrixes and uses significantly lower computational times. FastTree provides local support values based on the Shimodaira Hasegawa (SH) test (DOI: https://doi.org/10.1093/oxfordjournals.molbev.a026201) (accessed on 10 September 2020) with 1,000 bootstrap replicates to estimate the confidence in the given split. The graphical representation of the resulting phylogeny was obtained using the iTol interactive online tool (https://itol.embl.de/) (accessed on 21 September 2020) [44]. Variants were called in relation to the reference genome of the K. pneumoniae strain 30660/NJST258_1 (NCBI accession CP006923.1). Sequencing data of the polymyxin-resistant and carbapenemase-producing K. pneumoniae isolates from this study are available in NCBI BioProjects: PRJNA657895.

5. Conclusions

Herein, we report that the endemic CG258 represented by ST258 and ST512 (17 of 24 strains) exhibit high-level resistance to polymyxin. Since isolates belonging to this clonal group are one of the main drivers of the expansion of blaKPC in the country, the massive use of polymyxins to treat infected patients with carbapenemase-producing K. pneumoniae may lead to selective pressure and increase resistance. Active surveillance is warranted to prevent further spread of this clone in the country and to track the appearance of more polymyxin resistance determinants, including mcr-1.

Supplementary Materials

The following are available online at https://www.mdpi.com/2079-6382/10/3/284/s1, Figure S1: Virulence factors and plasmid replicons in isolates of K. pneumoniae.

Author Contributions

Conceptualization, M.F.M., E.D.L.C., and M.V.V.; methodology, E.D.L.C., M.F.M., J.C.G.-B., L.J.R., and J.S.S.-G.; validation, M.F.M., T.M.A., J.C.G.-B., J.P., C.J.P., and M.V.V.; formal analysis, E.D.L.C. and M.F.M.; investigation, E.D.L.C., M.F.M., and T.M.A.; data curation, E.D.L.C., M.F.M., J.S.S.-G., L.J.R., J.P., and J.C.G.-B.; writing—original draft preparation, E.D.L.C.; writing—review and editing, M.F.M., E.D.L.C., T.M.A., J.C.G.-B., J.P., and M.V.V.; supervision, E.D.L.C., M.F.M., J.C.G.-B., and M.V.V.; project administration, E.D.L.C.; and funding acquisition, E.D.L.C. and M.V.V. All authors have read and agreed to the published version of the manuscript.

Funding

This research was funded by Universidad El Bosque, grant number PCI 2018–10106.

Institutional Review Board Statement

Ethics approval was not required for this study.

Data Availability Statement

Sequencing data of the polymyxin-resistant and carbapenemase-producing K. pneumoniae isolates from this study are available in NCBI BioProject: PRJNA657895.

Acknowledgments

We would like to thank the Centro Internacional de Entrenamiento e Investigaciones Médicas (CIDEIM-Cali, Colombia) for kindly providing some of the strains used in this study.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Van Duin, D.; Perez, F.; Rudin, S.D.; Cober, E.; Hanrahan, J.; Ziegler, J.; Webber, R.; Fox, J.; Mason, P.; Richter, S.S.; et al. Surveillance of carbapenem-resistant Klebsiella pneumoniae: Tracking molecular epidemiology and outcomes through a regional network. Antimicrob. Agents Chemother. 2014, 58, 4035–4041. [Google Scholar] [CrossRef] [Green Version]
  2. Rojas, L.J.; Salim, M.; Cober, E.; Richter, S.S.; Perez, F.; Salata, R.A.; Kalayjian, R.C.; Watkins, R.R.; Marshall, S.; Rudin, S.D.; et al. Colistin Resistance in Car-bapenem-Resistant Klebsiella pneumoniae: Laboratory Detection and Impact on Mortality. Clin. Infect. Dis. 2017, 64, 711–718. [Google Scholar] [PubMed] [Green Version]
  3. Queenan, A.M.; Bush, K. Carbapenemases: The Versatile β-Lactamases. Clin. Microbiol. Rev. 2007, 20, 440–458. [Google Scholar] [CrossRef] [Green Version]
  4. Xu, L.; Sun, X.; Ma, X. Systematic review and meta-analysis of mortality of patients infected with carbapenem-resistant Klebsiella pneumoniae. Ann. Clin. Microbiol. Antimicrob. 2017, 16, 1–12. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  5. Zavascki, A.P.; Goldani, L.Z.; Li, J.; Nation, R.L. Polymyxin B for the treatment of multidrug-resistant pathogens: A critical review. J. Antimicrob. Chemother. 2007, 60, 1206–1215. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  6. Biswas, S.; Brunel, J.M.; Dubus, J.C.; Reynaud-Gaubert, M.; Rolain, J.M. Colistin: An update on the antibiotic of the 21st century. Expert Rev. Anti-Infect. Ther. 2012, 10, 917–934. [Google Scholar] [CrossRef] [PubMed]
  7. Jeannot, K.; Bolard, A.; Plésiat, P. Resistance to polymyxins in Gram-negative organisms. Int. J. Antimicrob. Agents 2017, 49, 526–535. [Google Scholar] [CrossRef] [PubMed]
  8. Cannatelli, A.; Giani, T.; D’Andrea, M.M.; Di Pilato, V.; Arena, F.; Conte, V.; Tryfinopoulou, K.; Vatopoulos, A.; Rossolini, G.M. MgrB Inactivation Is a Common Mechanism of Colistin Resistance in KPC-Producing Klebsiella pneumoniae of Clinical Origin. Antimicrob. Agents Chemother. 2014, 58, 5696–5703. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  9. Jana, B.; Cain, A.K.; Doerrler, W.T.; Boinett, C.J.; Fookes, M.C.; Parkhill, J.; Guardabassi, L. The secondary resistome of multidrug-resistant Klebsiella pneumoniae. Sci. Rep. 2017, 7, srep42483. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  10. Rojas, L.J.; Weinstock, G.M.; De La Cadena, E.; Diaz, L.; Rios, R.; Hanson, B.M.; Brown, J.S.; Vats, P.; Phillips, D.S.; Nguyen, H.; et al. An Analysis of the Epidemic of Klebsiella pneumoniae Carbapenemase-Producing K. pneumoniae: Convergence of Two Evolutionary Mechanisms Creates the “Perfect Storm”. J. Infect. Dis. 2018, 217, 82–92. [Google Scholar] [CrossRef] [Green Version]
  11. Poirel, L.; Jayol, A.; Bontron, S.; Villegas, M.-V.; Ozdamar, M.; Türkoglu, S.; Nordmann, P. The mgrB gene as a key target for acquired resistance to colistin in Klebsiella pneumoniae. J. Antimicrob. Chemother. 2014, 70, 75–80. [Google Scholar] [CrossRef] [PubMed]
  12. Bogdanovich, T.; Adams-Haduch, J.M.; Tian, G.-B.; Nguyen, M.H.; Kwak, E.J.; Muto, C.A.; Doi, Y. Colistin-Resistant, Klebsiella pneumoniae Carbapenemase (KPC)-Producing Klebsiella pneumoniae Belonging to the International Epidemic Clone ST258. Clin. Infect. Dis. 2011, 53, 373–376. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  13. Behera, B.; Mathur, P.; Das, A.; Kapil, A.; Gupta, B.; Bhoi, S.; Farooque, K.; Sharma, V.; Misra, M. Evaluation of susceptibility testing methods for polymyxin. Int. J. Infect. Dis. 2010, 14, e596–e601. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  14. Magiorakos, A.-P.; Srinivasan, A.; Carey, R.B.; Carmeli, Y.; Falagas, M.E.; Giske, C.G.; Harbarth, S.; Hindler, J.F.; Kahlmeter, G.; Olsson-Liljequist, B.; et al. Multi-drug-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, 268–281. [Google Scholar]
  15. Cannatelli, A.; D’Andrea, M.M.; Giani, T.; Di Pilato, V.; Arena, F.; Ambretti, S.; Gaibani, P.; Rossolini, G.M. In vivo emergence of colistin resistance in Klebsiella pneumoniae producing KPC-type carbapenemases mediated by insertional inactivation of the PhoQ/PhoP mgrB regulator. Antimicrob. Agents Chemother. 2013, 57, 5521–5526. [Google Scholar] [CrossRef] [Green Version]
  16. Jaidane, N.; Bonnin, R.A.; Mansour, W.; Girlich, D.; Creton, E.; Cotellon, G.; Chaouch, C.; Boujaafar, N.; Bouallegue, O.; Naas, T. Genomic Insights into Colistin-Resistant Klebsiella pneumoniae from a Tunisian Teaching Hospital. Antimicrob. Agents Chemother. 2017, 62. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  17. Jayol, A.; Poirel, L.; Villegas, M.-V.; Nordmann, P. Modulation of mgrB gene expression as a source of colistin resistance in Klebsiella oxytoca. Int. J. Antimicrob. Agents 2015, 46, 108–110. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  18. Jayol, A.; Nordmann, P.; Brink, A.; Villegas, M.-V.; Dubois, V.; Poirel, L. High-Level Resistance to Colistin Mediated by Various Mutations in the crrB Gene among Carbapenemase-Producing Klebsiella pneumoniae. Antimicrob. Agents Chemother. 2017, 61, 1–4. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  19. Kumar, A.; Biswas, L.; Omgy, N.; Mohan, K.; Vinod, V.; Sajeev, A.; Nair, P.; Singh, S.; Biswas, R. Colistin resistance due to inser-tional inactivation of the mgrB in klebsiella pneumoniae of clinical origin: First report from India. Rev. Esp. Quimioter. 2018, 31, 406–410. [Google Scholar]
  20. Novović, K.; Trudić, A.; Brkić, S.; Vasiljević, Z.; Kojić, M.; Medić, D.; Ćirković, I.; Jovčić, B. Molecular Epidemiology of Colistin-Resistant, Carbapenemase-Producing Klebsiella pneumoniae in Serbia from 2013 to 2016. Antimicrob. Agents Chemother. 2017, 61, e02550-16. [Google Scholar] [CrossRef] [Green Version]
  21. Olaitan, A.O.; Diene, S.M.; Kempf, M.; Berrazeg, M.; Bakour, S.; Gupta, S.K.; Thongmalayvong, B.; Akkhavong, K.; Somphavong, S.; Paboriboune, P.; et al. Worldwide emergence of colistin resistance in Klebsiella pneumoniae from healthy humans and pa-tients in Lao PDR, Thailand, Israel, Nigeria and France owing to inactivation of the PhoP/PhoQ regulator mgrB: An epidemi-ological and molecular study. Int. J. Antimicrob. Agents 2014, 44, 500–507. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  22. Pitt, M.E.; Elliot, A.G.; Cao, M.D.; Ganesamoorthy, D.; Karaiskos, I.; Giamarellou, H.; Abboud, C.S.; Blaskovich, M.A.T.; Cooper, M.A.; Coin, L.J.M. Multifactorial chromosomal variants regulate polymyxin resistance in extensively drug-resistant Klebsiella pneumoniae. Microb. Genom. 2018, 4, e000158. [Google Scholar] [CrossRef] [PubMed]
  23. Cheng, Y.-H.; Lin, T.-L.; Pan, Y.-J.; Wang, Y.-P.; Lin, Y.-T.; Wang, J.-T. Colistin Resistance Mechanisms in Klebsiella pneumoniae Strains from Taiwan. Antimicrob. Agents Chemother. 2015, 59, 2909–2913. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  24. Jayol, A.; Poirel, L.; Brink, A.; Villegas, M.-V.; Yilmaz, M.; Nordmann, P. Resistance to Colistin Associated with a Single Amino Acid Change in Protein PmrB among Klebsiella pneumoniae Isolates of Worldwide Origin. Antimicrob. Agents Chemother. 2014, 58, 4762–4766. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  25. Cheong, H.S.; Kim, S.Y.; Seo, J.; Wi, Y.M.; Peck, K.R.; Ko, K.S. Colistin Resistance and Extensive Genetic Variations in PmrAB and PhoPQ in Klebsiella Pneumoniae Isolates from South Korea. Curr. Microbiol. 2020, 77, 2307–2311. [Google Scholar] [CrossRef]
  26. Jayol, A.; Poirel, L.; Dortet, L.; Nordmann, P. 2016. National survey of colistin resistance among carbapenemase-producing Enterobacteriaceae and outbreak caused by colistin-resistant OXA-48-producing Klebsiella pneumoniae, France, 2014. Eurosurveillance 2016, 21, 30339. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  27. Mojica, M.F.; Correa, A.; Vargas, D.A.; Maya, J.J.; Montealegre, M.C.; Rojas, L.J.; Ruiz, S.J.; Quinn, J.P.; Villegas, M.V.; Colombian Nosocomial Bacterial Resistance Study Group. Molecular correlates of the spread of KPC-producing Enterobacteriaceae in Colombia. Int. J. Antimicrob. Agents 2012, 40, 277–279. [Google Scholar] [CrossRef] [PubMed]
  28. Ocampo, A.M.; Vargas, C.A.; Sierra, P.M.; Cienfuegos, A.V.; Jiménez, J.N. Caracterización molecular de un brote de Klebsiella pneumoniae resistente a carbapenémicos en un hospital de alto nivel de complejidad de Medellín, Colombia. Biomédica 2015, 35, 496–504. [Google Scholar] [CrossRef] [Green Version]
  29. Lalaoui, R.; Bakour, S.; Livnat, K.; Assous, M.V.; Diene, S.M.; Rolain, J.-M. Spread of Carbapenem and Colistin-Resistant Klebsiella pneumoniae ST512 Clinical Isolates in Israel: A Cause for Vigilance. Microb. Drug Resist. 2019, 25, 63–71. [Google Scholar] [CrossRef]
  30. Giani, T.; Pini, B.; Arena, F.; Conte, V.; Bracco, S.; Migliavacca, R.; Pantosti, A.; Pagani, L.; Luzzaro, F.; Rossolini, G.M. Epidemic diffusion of KPC carbapenemase-producing Klebsiella pneumoniae in Italy: Results of the first countrywide survey, 15 May to 30 June 2011. Eurosurveillance 2013, 18, 1–9. [Google Scholar]
  31. Cienfuegos-Gallet, A.V.; Chen, L.; Kreiswirth, B.N.; Jiménez, J.N. Colistin Resistance in Carbapenem-Resistant Klebsiella pneumoniae Mediated by Chromosomal Integration of Plasmid DNA. Antimicrob. Agents Chemother. 2017, 61, e00404-17. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  32. Saavedra, S.Y.; Diaz, L.; Wiesner, M.; Correa, A.; Arévalo, S.A.; Reyes, J.; Hidalgo, A.M.; De La Cadena, E.; Perenguez, M.; Montaño, L.A.; et al. Genomic and Molecular Characterization of Clinical Isolates of Enterobacteriaceae Harboring mcr-1 in Colombia, 2002 to 2016. Antimicrob. Agents Chemother. 2017, 61, e00841-17. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  33. Berglund, B. Acquired Resistance to Colistin via Chromosomal And Plasmid-Mediated Mechanisms in Klebsiella pneumoniae. Infect. Microbes. Dis. 2019, 1, 10–19. [Google Scholar] [CrossRef]
  34. CLSI. Performance Standards for Antimicrobial Susceptibility Testing, 28th ed.; CLSI Supplement M100; Clinical and Laboratory Standards Institute: Wayne, PA, USA, 2018. [Google Scholar]
  35. CLSI. Methods for Dilution Antimicrobial Susceptibility Tests for Bacteria That Grow Aerobically; Approved Standard-Thenth Edition; CLSI Document M07-A10; Clinical and Laboratory Standards Institute: Wayne, PA, USA, 2015. [Google Scholar]
  36. Correa, A.; Del Campo, R.; Perenguez, M.; Blanco, V.M.; Rodríguez-Baños, M.; Perez, F.; Maya, J.J.; Rojas, L.; Cantón, R.; Arias, C.A.; et al. Dissemination of High-Risk Clones of Extensively Drug-Resistant Pseudomonas aeruginosa in Colombia. Antimicrob. Agents Chemother. 2015, 59, 2421–2425. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  37. Davis, J.J.; Wattam, A.R.; Aziz, R.K.; Brettin, T.; Butler, R.; Butler, R.M.; Chlenski, P.; Conrad, N.; Dickerman, A.; Dietrich, E.M.; et al. The PATRIC Bioinformatics Resource Center: Expanding data and analysis capabilities. Nucleic Acids Res. 2019, 48, D606–D612. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  38. Zankari, E.; Hasman, H.; Cosentino, S.; Vestergaard, M.; Rasmussen, S.; Lund, O.; Aarestrup, F.M.; Larsen, M.V. Identification of acquired antimicrobial resistance genes. J. Antimicrob. Chemother. 2012, 67, 2640–2644. [Google Scholar] [CrossRef]
  39. Carattoli, A.; Zankari, E.; García-Fernández, A.; Larsen, M.V.; Lund, O.; Villa, L.; Aarestrup, F.M.; Hasman, H. In SilicoDetection and Typing of Plasmids using PlasmidFinder and Plasmid Multilocus Sequence Typing. Antimicrob. Agents Chemother. 2014, 58, 3895–3903. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  40. Larsen, M.V.; Cosentino, S.; Rasmussen, S.; Friis, C.; Hasman, H.; Marvig, R.L.; Jelsbak, L.; Sicheritz-Pontén, T.; Ussery, D.W.; Aarestrup, F.M.; et al. Multilocus sequence typing of total-genome-sequenced bacteria. J. Clin. Microbiol. 2012, 50, 1355–1361. [Google Scholar] [CrossRef] [Green Version]
  41. Chan, M.-S.; Maiden, M.C.J.; Spratt, B.G. Database-driven Multi Locus Sequence Typing (MLST) of bacterial pathogens. Bioinformatics 2001, 17, 1077–1083. [Google Scholar] [CrossRef]
  42. Siguier, P.; Perochon, J.; Lestrade, L.; Mahillon, J.; Chandler, M. ISfinder: The reference centre for bacterial insertion sequenc-es. Nucleic Acids Res. 2006, 34, 32–36. [Google Scholar] [CrossRef] [Green Version]
  43. Kaas, R.S.; Leekitcharoenphon, P.; Aarestrup, F.M.; Lund, O. Solving the Problem of Comparing Whole Bacterial Genomes across Different Sequencing Platforms. PLoS ONE 2014, 9, e104984. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  44. Letunic, I.; Bork, P. Interactive Tree Of Life (iTOL) v4: Recent updates and new developments. Nucleic Acids Res. 2019, 47, W256–W259. [Google Scholar] [CrossRef] [PubMed] [Green Version]
Antibiotics 10 00284 g001 550
Figure 1. Genomic characteristics of the 24 K. pneumoniae isolates included in this study. C1: Cali, M1: institution 1 Medellin, M2: institution 2 Medellin, B1: Institution 1 Bogotá, B2: institution 2 Bogotá, B3: institution 3 Bogotá, N1: Neiva, I1: Ibague, H1: Pasto. Black indicates presence, * not-Tn4401.
Figure 1. Genomic characteristics of the 24 K. pneumoniae isolates included in this study. C1: Cali, M1: institution 1 Medellin, M2: institution 2 Medellin, B1: Institution 1 Bogotá, B2: institution 2 Bogotá, B3: institution 3 Bogotá, N1: Neiva, I1: Ibague, H1: Pasto. Black indicates presence, * not-Tn4401.
Antibiotics 10 00284 g001
Antibiotics 10 00284 g002 550
Figure 2. Characteristics of polymyxin resistant K. pneumoniae isolates included in this study. C1: Cali, M1: institution 1 Medellin, M2: institution 2 Medellin, B1: Institution 1 Bogotá, B2: institution 2 Bogotá, B3: institution 3 Bogotá, N1: Neiva, I1: Ibague, H1: Pasto. Mutations previously reported associated with polymyxin resistance are highlighted in red. Q30STOP: Premature stop codon at 30th amino acid. The interactive tree of life (iTOL) tool was used for the phylogenetic tree display and annotation. Bootstrap values shown in Supplemental Figure S2.
Figure 2. Characteristics of polymyxin resistant K. pneumoniae isolates included in this study. C1: Cali, M1: institution 1 Medellin, M2: institution 2 Medellin, B1: Institution 1 Bogotá, B2: institution 2 Bogotá, B3: institution 3 Bogotá, N1: Neiva, I1: Ibague, H1: Pasto. Mutations previously reported associated with polymyxin resistance are highlighted in red. Q30STOP: Premature stop codon at 30th amino acid. The interactive tree of life (iTOL) tool was used for the phylogenetic tree display and annotation. Bootstrap values shown in Supplemental Figure S2.
Antibiotics 10 00284 g002
Table
Table 1. Characteristics and minimum inhibitory concentrations of polymyxin-resistant K. pneumoniae (mg/L).
Table 1. Characteristics and minimum inhibitory concentrations of polymyxin-resistant K. pneumoniae (mg/L).
IDHOSPITALYEARSOURCEROOMPOLCZACAZFEPTZPETPIMIMEMDORCIPTGC
UEB_1M12009BLOODHOSP44/4>32>64>64/4>32>32>32>16>24
UEB_2C12010SKIN ER8≤1/4>32>64>64/4>3232>32>16>2≤0.5
UEB_3C12010BLOODICU84/4>32>64>128/4>3232>32>16>2≤0.5
UEB_4C12010SKIN ICU>84/4>32>64>128/4>32>32>32>16>21
UEB_5B12010OTHERHOSP4≤1/4>3264>128/42111>2≤0.5
UEB_6M12011GIICU>82/4>3216>128/4164441≤0.5
UEB_7M12011SECRETIONICU84/4>3216>128/4>3216>32>1611
UEB_8M12011GIICU>82/4>32>64>128/4>321632>1611
UEB_9M12011GIICU8≤1/4>3216128/4211≤0.511
UEB_10C12011GIICU>8≤1/4>3264>128/43241616>22
UEB_11B22011BLOODICU>84/4>33>64>128/4>32>32>32>16>22
UEB_12M22011RESP TRACTICU>84/4>32>64>128/4>32>32>32>16>21
UEB_13M22011BLOODICU>84/4>3264>128/4>32>32>32>16>2≤0.5
UEB_14M22011GIICU>82/4>3264>128/4>32>32>32>16>21
UEB_15M22012URINEHOSP>84/4>3264>128/4>32>32>32>16>21
UEB_16C12014BLOODHOSP>8≤1/4>32>64>128/432>32>32>16>24
UEB_17H12015RESP TRACTICU>82/4>3264>128/4>32>32>32>1612
UEB_18N12017BLOODER>8≤1/4>3232>128/44212≤0.254
UEB_19H12017URINEHOSP>8≤1/4>3264>128/4>323232>16>24
UEB_20I12017URINEHOSP>8≤1/4>32>64>128/4328161612
UEB_21B32019RESP TRACTHOSP>84/4>32>64>128/4>32>32>32>16>21
UEB_22B32019PERIT LIQHOSP>82/4>32>64>128/4>32>32>32>16>21
UEB_23B32019ASCIT LIQHOSP>84/4>32>64>128/4>32>32>32>16>21
UEB_24I12016RESP TRACTHOSP>8>64/4>32>64>128/4>32>3232>16>22
POL: polymyxin, CZA: ceftazidime-avibactam, CAZ: ceftazidime, FEP: cefepime, TZP: piperacillin-tazobactam, ETP: ertapenem, IMI: imipenem, MEM: meropenem, DOR: doripenem, CIP: ciprofloxacin, and TGC: tygecicline. C1: Cali, M1: institution 1 Medellin, M2: institution 2 Medellin, B1: Institution 1 Bogotá, B2: institution 2 Bogotá, B3: institution 3 Bogotá, N1: Neiva, I1: Ibague, H1: Pasto. ICU: intensive care unit, ER: emergency room, and HOSP: hospitalization room. GI: gastrointestinal, SKIN: Skin and Soft Tissue, RESP TRACT: respiratory Tract, PERIT LIQU: peritoneal liquid, ASCITIC LIQ: ascitic liquid.
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

De La Cadena, E.; Mojica, M.F.; García-Betancur, J.C.; Appel, T.M.; Porras, J.; Pallares, C.J.; Solano-Gutiérrez, J.S.; Rojas, L.J.; Villegas, M.V. Molecular Analysis of Polymyxin Resistance among Carbapenemase-Producing Klebsiella pneumoniae in Colombia. Antibiotics 2021, 10, 284. https://doi.org/10.3390/antibiotics10030284

AMA Style

De La Cadena E, Mojica MF, García-Betancur JC, Appel TM, Porras J, Pallares CJ, Solano-Gutiérrez JS, Rojas LJ, Villegas MV. Molecular Analysis of Polymyxin Resistance among Carbapenemase-Producing Klebsiella pneumoniae in Colombia. Antibiotics. 2021; 10(3):284. https://doi.org/10.3390/antibiotics10030284

Chicago/Turabian Style

De La Cadena, Elsa, María Fernanda Mojica, Juan Carlos García-Betancur, Tobías Manuel Appel, Jessica Porras, Christian José Pallares, Juan Sebastián Solano-Gutiérrez, Laura J. Rojas, and María Virginia Villegas. 2021. "Molecular Analysis of Polymyxin Resistance among Carbapenemase-Producing Klebsiella pneumoniae in Colombia" Antibiotics 10, no. 3: 284. https://doi.org/10.3390/antibiotics10030284

APA Style

De La Cadena, E., Mojica, M. F., García-Betancur, J. C., Appel, T. M., Porras, J., Pallares, C. J., Solano-Gutiérrez, J. S., Rojas, L. J., & Villegas, M. V. (2021). Molecular Analysis of Polymyxin Resistance among Carbapenemase-Producing Klebsiella pneumoniae in Colombia. Antibiotics, 10(3), 284. https://doi.org/10.3390/antibiotics10030284

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop