Next Article in Journal
Is Antibiotic Prophylaxis Necessary before Dental Implant Procedures in Patients with Orthopaedic Prostheses? A Systematic Review
Next Article in Special Issue
The Evaluation of Eazyplex® SuperBug CRE Assay Usefulness for the Detection of ESBLs and Carbapenemases Genes Directly from Urine Samples and Positive Blood Cultures
Previous Article in Journal
Antimicrobial Photosensitizing Material Based on Conjugated Zn(II) Porphyrins
Previous Article in Special Issue
The Antimicrobial Resistance Characteristics of Imipenem-Non-Susceptible, Imipenemase-6-Producing Escherichia coli
 
 
Search for Articles:
Title / Keyword
Author / Affiliation / Email
Journal
Article Type
 
 
Advanced Search
 
Section
Special Issue
Volume
Issue
Number
Page
 
Journals
Antibiotics
Volume 11
Issue 1
10.3390/antibiotics11010092
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

NDM-1 Introduction in Portugal through a ST11 KL105 Klebsiella pneumoniae Widespread in Europe

by
Ângela Novais
1,2,*,
Rita Veiga Ferraz
3,
Mariana Viana
3,
Paula Martins da Costa
3,4 and
Luísa Peixe
1,2,*
1
UCIBIO, Laboratório de Microbiologia, Departamento de Ciências Biológicas, Faculdade de Farmácia, Universidade do Porto, 4050-313 Porto, Portugal
2
Associate Laboratory i4HB—Institute for Health and Bioeconomy, Faculty of Pharmacy, University of Porto, 4050-313 Porto, Portugal
3
Centro Hospitalar do Tâmega e Sousa, 4564-007 Penafiel, Portugal
4
Department of Clinical Chemistry, Centro Hospitalar Universitário do Porto, 4099-001 Porto, Portugal
*
Authors to whom correspondence should be addressed.
Antibiotics 2022, 11(1), 92; https://doi.org/10.3390/antibiotics11010092
Submission received: 19 November 2021 / Revised: 5 January 2022 / Accepted: 8 January 2022 / Published: 12 January 2022
(This article belongs to the Special Issue Carbapenemase-Producing Enterobacterales)
Browse Figures
Versions Notes

Abstract

:
The changing epidemiology of carbapenem-resistant Klebsiella pneumoniae in Southern European countries is challenging for infection control, and it is critical to identify and track new genetic entities (genes, carbapenemases, clones) quickly and with high precision. We aimed to characterize the strain responsible for the first recognized outbreak by an NDM-1-producing K. pneumoniae in Portugal, and to elucidate its diffusion in an international context. NDM-1-producing multidrug-resistant K. pneumoniae isolates from hospitalized patients (2018–2019) were characterized using FTIR spectroscopy, molecular typing, whole-genome sequencing, and comparative genomics with available K. pneumoniae ST11 KL105 genomes. FT-IR spectroscopy allowed the rapid (ca. 4 h after incubation) identification of the outbreak strains as ST11 KL105, supporting outbreak control. Epidemiological information supports a community source but without linkage to endemic regions of NDM-1 producers. Whole-genome comparison with previous DHA-1-producing ST11 KL105 strains revealed the presence of different plasmid types and antibiotic resistance traits, suggesting the entry of a new strain. In fact, this ST11 KL105 clade has successfully disseminated in Europe with variable beta-lactamases, but essentially as ESBL or DHA-1 producers. We expand the distribution map of NDM-1-producing K. pneumoniae in Europe, at the expense of a successfully established ST11 KL105 K. pneumoniae clade circulating with variable plasmid backgrounds and beta-lactamases. Our work further supports the use of FT-IR as an asset to support quick infection control.

1. Introduction

The dissemination of carbapenem resistant K. pneumoniae (CRKP) is a well-known problem, but there is a great asymmetry in the geographic distribution of carbapenemase types and K. pneumoniae lineages across countries [1,2,3]. Moreover, changing epidemiology over time is a very challenging situation that affects several countries from Southern Europe [4,5,6,7]. Knowing and to understanding these dynamics is essential in order to promptly recognize new genetic entities with enhanced antibiotic resistance, virulence and/or transmission, and to allow a prompt redesign of diagnostic tools and infection control policies to curtail further dissemination [2,8]. Whole-genome sequencing has provided the greatest resolution in establishing transmission pathways at local and global levels of CRKP [1,2], and in enlightening the phylogeny and evolutionary history of some highly frequent KP clonal groups (e.g., CG258, CG307) [8,9,10]. However, it is still unable to support real-time infection control decisions.
ST11 is one of the most diffused clones, and is especially prevalent in Asia and particularly in China [11]. The heterogeneity within ST11 is well recognized in studies on this geographic area, with clades reflecting capsule recombination events (KL64, KL47, KL105) [9,12]. For this reason, because Fourier transform infrared (FT-IR) spectroscopy is a suitable tool to identify K. pneumoniae capsular (KL) types with a fast turnaround time [13,14], it is ideal to discriminate and identify ST11 variants, which are understudied in Europe.
In Portugal, resistance to carbapenems is a problem with growing proportions, especially among K. pneumoniae, of which the rates among invasive isolates recorded in national surveillance networks increased ca. 80% in just three years (2014–2017) [15]. Reporting of CPE from national laboratories has been mandatory since 2013 and there have been official regulations for the management of outbreaks and for CPE infection control since 2017 [4]. The available epidemiological and population analysis data suggest the dominance all over the country of particular genetic traits (the ST147 clone, KPC-3 carbapenemase) [16,17,18,19], with an increasing diversity of species, carbapenemases types and clones [2,4,20]. KPC-3 and to a lesser extent OXA-48 or OXA-181, which are occasionally associated with GES-5 [18,20], are the most common carbapenemase types among K. pneumoniae. NDM-1 has been sporadically described in Providencia stuartii, Morganella morganii and Proteus mirabilis [21,22] but, prior to this study, no cases of colonization and/or infection with NDM-producing K. pneumoniae were reported. NDM-1 is especially prevalent in India, the Middle East and the Balkans, and its emergence in other countries (especially in Escherichia coli or K. pneumoniae) is considered a public health priority due to the risk of subsequent (frequently polyclonal) dissemination [3,23,24].
In this study, we used FT-IR spectroscopy to support the identification and control of the first recognized outbreak of NDM-1-producing ST11 capsular-type KL105 K. pneumoniae in Portugal. Furthermore, we used whole-genome sequencing to trace the variability and diffusion of this particular ST11 K. pneumoniae clade from an international perspective.

2. Results

2.1. Epidemiological Context and Infection Control Measures

We identified and characterized the first recognized outbreak of NDM-1-producing K. pneumoniae isolates in Portugal, involving a total of seven patients. The index case was from a patient (patient 1) with no history of previous hospitalization or travel abroad, whereas most other cases were linked to concomitant hospitalization in the surgical unit at the same time as the index patient (patients 2–6) (Table 1). None of the other patients had a record of foreign travel in the 3 months prior to strain isolation. Patient 3 had been previously hospitalized in a different hospital but screening at admission was negative. Patient 7 came from a household and no epidemiological link could be established with other patients, but this patient had also been previously hospitalized and tested negative at first admission.
The patient cohort, screening of putative carriers, contact precautions and the reinforcement of cleaning and disinfection practices prevented new cases. Regular screening of CPE carriers at admission and during hospitalization, as well as weekly audits to assess compliance to the protocols and implement corrective measures, were used to monitor for the appearance of new cases, and these are still in practice to control the spread of CPE.
In the three months prior to NDM-1 detection, the patients who received antibiotherapy were treated with beta-lactam/beta-lactamase inhibitors, cefuroxime, ceftazidime and ciprofloxacin.

2.2. Strain Relatedness

FT-IR spectroscopy was performed retrospectively on suspected isolates received at the Faculty of Pharmacy, University of Porto. After reception, our FT-IR based workflow allowed us to establish a close relationship between the first five isolates, and to communicate strain relatedness less than 4 h after bacterial culture. Comparison of the spectra with our database allowed the classification of the K-type as KL105, together with previous DHA-1 producing ST11-KL105 isolates that had been circulating in Portuguese hospitals for years (Figure S1) [25]. The four isolates identified subsequently clustered with these, confirming the K-type and strain relatedness by means of FTIR, which was further confirmed using PFGE for all studied isolates (data not shown). The genes blaNDM-1, as well as blaCTX-M-15, were detected in all isolates via PCR, and all isolates exhibited identical susceptibility profiles.

2.3. Molecular Characterization

Isolates revealed multidrug resistance profiles, being resistant to all antibiotics tested except amikacin, tigecycline and colistin. Sequenced isolates carried several antibiotic resistance genes (blaNDM-1, blaSHV-182, blaCTX-M-15, aac(6′)-Ib-cr, aac(3)-IId, aph(3′’)-Ib, aph(6)-Id, oqxAB, qnrB1, catB3, sul2, tetD and dfrA14), supporting the phenotypic data. They also carried chromosomal mutations in gyrA (S83I) and parC (S80I), and UhpT (E350Q), conferring reduced susceptibility to fluoroquinolones and fosfomycin, respectively. They were also enriched with virulence genes encoding enterobactin (entB), yersiniabactin (irp1/2, ybt10 associated with ICEKp4), iron receptors (fyuA, kfuBC) and type I and type III fimbriae (fimH, mrkABC). All strains had wzi75, which corresponds to a novel capsular type (KL105), and KLEBORATE confirmed 100% identity with the KL105 locus and the O2v2 O antigen locus.

2.4. Whole-Genome Data and Phylogenetic Analysis

The genomes of K606 and K624 strains had 5.5 Gb, a GC content of 57.2% and were assembled into 76–97 contigs (data not shown). The SNP-based phylogenetic tree was constructed on the basis of 4.87 Gb common positions covered by the core genome, corresponding to 88.5% of the genome of the reference strain. We observed a close relatedness between the available ST11 KL105 genomes from Europe, America and China (2013–2017) since isolates varied at 0–297 sites distributed over their 4.8 million base pair (bp) core genome (~medium 31SNP/Gb) (Figure 1, Supplementary Figure S2). All of them carried blaSHV-182 in the chromosome.
Two main branches were identified within the tree (Figure 1). One of these included highly related (<20 SNPs) KPC-2-producing strains from China and carried IncR and ColRNAI plasmids and a few additional antibiotic resistance genes. The second branch included isolates producing CTX-M-3, CTX-M-15 or DHA-1, occasionally together with NDM-1, most of them (n = 13/15; 87%) identified in Europe, and many from the Euscape project [2]. Most (92%) of these carried several antibiotic resistance genes, some of which were signatures of particular genetic platforms (e.g., catB3, arr-3, qnrB4). Most carried R replicons, but also FIA(HI1), FIIk and/or FIB, together with variable F replicons (FIA, FIB and/or FIIk, and colRNAI) (Figure 1). All of them consistently carried ybt10 (ICEKp4), ipr1/irp2 and fyuA2 virulence genes.
The SNP-based matrix (supplementary Figure S2) showed that the two sequenced strains differ in only six SNPs, supporting their close relatedness. However, due to the absence of an epidemiological link between both patients, it is reasonable to consider an unrecognized source. Both strains showed ~100 SNPs (20 SNP/Gb) with previous DHA-1 ST11-KL105 producers circulating in Portuguese hospitals (Whole Genome Shotgun project accession Nº QTTC00000000 and QTTD00000000). Three of the public genomes, corresponding to isolates from Romania and the USA, also carried blaNDM-1 and showed only 80–90 SNP differences with our outbreak strains.

2.5. Genetic Support for blaNDM-1 Acquisition

The blaNDM-1 gene was identified in pAN_K624, an 84 Kb IncR+IncFIA(HI1) multirreplicon plasmid with high identity (100%) and coverage (92%) to the pAR_0109 plasmid (Genbank acc. Number CP032212), which is of an unknown origin and which carries the typical genetic environment of blaCTX-M-15. blaNDM-1 is located in a Tn125-like multidrug resistance region (MRR) comprising deltaISAba125, blaNDM-1, bleMBL, trpF, dsbC, cutA, groES, groEL, psp, qnrB1 and SDR family reductase, 99% identical to that of the reference plasmid pNDM-MAR (Genbank acc. Number JN420336). The MRR region (~8.5Kb) containing blaNDM-1 was flanked by an IS3000 element that might have mediated its insertion into the backbone of a pAR_0109-like plasmid between a Tn5403-like transposon and pspF (Figure 2).
Additionally, these strains also carried a 149 Kb IncFIB(K) plasmid with high identity and coverage (>99%) with the ~160 Kb p002SK2_A plasmid (Genbank acc. Number CP025516), a non-antibiotic resistant plasmid identified in a K. pneumoniae ST147 isolate from wastewater in Switzerland. There were no traces of the DHA-1-producing pKPS30 IncR plasmid previously identified in ST11-KL105 isolates circulating in Portugal [25], nor of the IncC plasmid previously detected with blaNDM-1 in species other than K. pneumoniae [21]. Likewise, none of the other NDM-1-producing ST11 KL105 public strains carried a similar plasmid replicon content, suggesting a new plasmid as the vehicle of blaNDM-1.

3. Discussion

This study enlarges the distribution map of NDM-1 carbapenemase in Europe and raises the alarm for its emergence in our country in a multidrug-resistant ST11-KL105 K. pneumoniae lineage, that has previously caused the long-term dissemination of DHA-1 in Portuguese hospitals [2,25]. Simultaneously, it increases the pool of carbapenemases circulating among high-risk K. pneumoniae lineages in our country, which is worrisome considering the possibility of further transmission of blaNDM-1 to other already circulating high-risk K. pneumoniae clones or even to other Enterobacterales [4,17]. In fact, other NDM-1-producing K. pneumoniae ST11-KL105 strains with identical PFGE patterns have been more recently identified in other hospitals in our area, suggesting further dissemination and wider expansion (data not shown).
It is of interest to highlight that none of the patients had a record of hospitalization or travel abroad in the 3 months prior to strain isolation, and the absence of an epidemiological link for the index case suggests community acquisition from an unrecognized source. A chain of nosocomial dissemination could be established, and our data further reinforce the importance of contact screening of CPE using fast and reliable methods such as FT-IR, together with quick identification of carbapenemases, in order to optimize infection control, to detect new emergences and prevent subsequent spread [4].
The quick identification of the outbreak using FT-IR spectroscopy (ca. 4 h after incubation) allowed the prompt establishment of infection control measures for outbreak control, and whole-genome sequencing supported the classification of the outbreak strain as ST11 KL105. Our comparative genomic analysis of available ST11 KL105 genomes revealed that this lineage has been found circulating in several European countries since at least 2006, in most cases carrying blaDHA-1, blaCTX-M-3 or blaCTX-M-15 (Figure 1). This is in agreement with data from the Euscape project [2], in which ca. 62% of ST11 K. pneumoniae isolates were non-carbapenemase producers. Here, we additionally demonstrate that these Euscape isolates belong to the same ST11 clade expressing KL105, which differs from variants prevalent in China and eastern European countries [9,12,26].
Our comparative genomic analysis shows that NDM-1 was not acquired by ST11 KL105 strains previously circulating in Portugal [25] but through the introduction of a new ST11 KL105 strain, carrying blaNDM-1 in a different genetic context. In fact, comparison of NDM-1- and DHA-1-producing ST11 strains revealed variability of beta-lactamases (blaNDM-1, blaCTX-M-15 vs blaDHA-1), antibiotic resistance traits (qnrB1, catB3, dfrA14 vs qnrB4, arr-3, catA1, dfrA12) and plasmids (IncR+FIA(HI1) + IncFIB(K) vs IncR). The pAN_K624 plasmid (IncR+FIA(HI1)) encoding NDM-1 was also different from the IncC platform identified previously in other Enterobacterales in Portugal [21], supporting a new entry of the blaNDM-1 gene and the circulation of different platforms encoding NDM-1. Furthermore, our data on the plasmid genetic context suggest the introduction of blaNDM-1 into a multidrug resistance plasmid that is different from those carrying blaNDM-1 in other ST11 KL105 strains from Europe, suggesting a new platform. Also, it differs from the pKPX-1 plasmid, associated with massive NDM-1 spread in Poland [26]. In the ST11 KL105 clade, the acquisition of blaNDM-1 seems to be sporadic, in some cases over DHA-1 or ESBL genetic backgrounds (Figure 1, Supplementary Table S1), contributing to high genetic plasticity and plasmid shuffling in this lineage, extending previous observations [27].
Our work also highlights the need to further deepen the population analysis of K. pneumoniae clones to obtain maximum discrimination for the efficient tracking of outbreak strains. Indeed, there is high genetic variability within some of the most frequent clonal groups, including CG11 [2,27,28]. Despite the high capsular diversity identified among the most common clonal groups, many studies have demonstrated an extraordinary selection of lineages with specific capsular types within CG258 [28] and CG11 [9,12], suggesting that the capsule type is a good evolutionary marker of clinically relevant K. pneumoniae lineages [13]. Considering this evolutionary scenario, FT-IR based typing based on a representative spectral database is a reliable tool to differentiate these lineages, providing information that, together with epidemiological data, might be of great value in supporting the identification of outbreaks and real-time infection control decisions at a quick and low-cost rate [13,14,29]. Whole-genome sequencing requires a higher cost and time but it can be used downstream to provide full resolution, allowing a thorough understanding of the dynamics of the spread of antimicrobial resistance [8,10,27].
This ST11-KL105 K. pneumoniae lineage consistently carries blaSHV-182, encoding a natural SHV-type beta-lactamase, according to the BLBD beta-lactamase database (http://bldb.eu accessed on 11 November 2019) [30]. In all the genomes analyzed, blaSHV-182 was chromosomally located and might be potentially used as a marker for this specific lineage.

4. Materials and Methods

4.1. Setting

The outbreak occurred at a 500-bed hospital that serves a population of about 520,000 in the north of Portugal, together with another hospital unit at a distance of 30 km. It offers healthcare specialties within medical, surgical and emergency (adults and pediatric), as well as outpatient attendance. Control of the dissemination of CPE bacteria is implemented according to official national recommendations and includes active screening and the reinforcement of infection control procedures. At this hospital, rectal screening of CPE is performed (i) at admission when transferring from another hospital with hospitalizations >48 h or from a household or long-term care facility; (ii) when there are previous (12 months) hospitalizations, (iii) every 7 days of hospitalization. At admission, genes encoding specific carbapenemases were detected using a rapid PCR kit (Xpert CarbaR), whereas subsequent screenings on hospitalized patients were performed using a cultural method.

4.2. Index Case

The first NDM-1-producing K. pneumoniae was isolated in the surgical ward in the bile of a 71-year-old patient with a cholangiocarcinoma on 27 November 2018. In the absence of clinical symptoms it was considered possible drain colonization, and CPE rectal screening performed the day after revealed the carriage of a NDM-1-producing isolate. This patient had been hospitalized since 26 October 2018 but there were no criteria for screening at admission, thus this information is missing.
These isolates were resistant to third-generation cephalosporins, carbapenems, aztreonam, beta-lactam/beta-lactamase inhibitors, fluoroquinolones, trimethoprim-sulfamethoxazole, nitrofurantoin and several aminoglycosides; intermediate to fosfomycin; and susceptible only to amikacin, tigecycline and colistin.

4.3. Additional Cases

Six additional patients were identified with NDM-1-encoding K. pneumoniae. All patients in the same nursery were considered exposed and were screened for rectal carriage independently on the result at admission (if applicable). Three of them carried a K. pneumoniae isolate with the same susceptibility profile. All patients with positive screening results were transferred to cohort areas and subjected to isolation and contact precautions. Four additional isolates were subsequently identified in the rectal swabs or caused a urinary infection in three outpatients that were admitted between December 2018 and March 2019 from the emergency room. Two of them were identified as carriers of NDM-producing K. pneumoniae at admission, whereas the last one was negative for carriage at admission during a previous hospitalization but that missed subsequent screenings.

4.4. Bacterial Identification and Carbapenemase Detection

In total, 9 NDM-1-encoding K. pneumoniae isolates were identified in seven patients. These isolates were recovered from feces (n = 9), urine (n = 1) and bile (n = 1) (Table 1). Carbapenemase genes were primarily detected via molecular biology using Xpert CarbaR (Cepheid). Bacterial identification and preliminary susceptibility testing were performed using a VITEK 2 system (bioMérieux) or broth microdilution (for colistin). Available clinical and epidemiological data, as well as information on recognized risk factors for infection by CPE, were recorded and analyzed.

4.5. Antibiotic Resistance Phenotypes and Genotypes

Antibiotic susceptibility profiles to 25 antibiotics (all beta-lactams, beta-lactam/beta-lactamase combinations, fluoroquinolones, aminoglycosides, tetracyclines, fosfomycin, trimethoprim-sulfamethoxazole, chloramphenicol) were determined via disk diffusion and interpreted following EUCAST guidelines (www.eucast.org 6 January 2020). Carbapenemase production was confirmed using the Blue-Carba test [31] and the type of carbapenemase was identified via PCR directed to the most frequent carbapenemase gene families (blaKPC, blaOXA-48, blaIMP; blaVIM, blaNDM) and further sequencing [17,32]. The presence of blaCTX-M-15 was confirmed using PCR in all isolates.

4.6. Identification of Outbreak Strains

We used Fourier transform infrared (FTIR) spectroscopy to assess isolates’ relationships according to our previously established workflow and models based on the comparison of bacterial discriminatory spectra. The analysis was performed retrospectively on suspected isolates received in early December 2018 and March 2019. Upon reception, isolates were re-grown in standardized culture conditions (37 °C/18 h), after which a colony was directly deposited on the ATR accessory of our FT-IR instrument (Perkin-Elmer Frontier) and air-dried. Three spectra per strain were acquired (technical replicates) from 600 cm−1 to 4000 cm−1, and the spectral region of interest (accumulating polysaccharide vibrations) was compared with spectra included in our in-house K. pneumoniae database and machine-learning model [13]. The K-type and clone of the suspected strains were then predicted with an accuracy of 100%, because of previously established correlations in internationally circulating lineages [13]. These FT-IR-based assignments were confirmed by means of PCR and sequencing of the wzi gene [33], as well as pulsed-field gel electrophoresis.

4.7. Whole Genome and Phylogenetic Analysis

Two NDM-1-producing isolates (the first and the last outbreak strain with no epidemiological link) were selected for whole genome sequencing. The DNA was extracted using the Wizard Genomic DNA purification kit (Promega Corporation, Madison, WI, USA) according to the manufacturer’s instructions and the concentration was determined with a Qubit 3.0 Fluorometer (Invitrogen, Thermo Fisher Scientific, Waltham, MA, USA). DNA was sequenced using an Illumina Novaseq 6000 (2× 300 bp pair-ended runs, ~6 Gb genome, coverage 100×), reads were trimmed de novo using Trimommatic 0.39 to remove adapters (http://www.usadellab.org/cms/?page=trimmomatic accessed on 5 November 2019) and assembled using SPAdes version 3.9.0 (cab.spbu.ru/software/spades/ accessed on 5 November 2019). The quality of the reads and the assembly were assessed using FastQC v0.11.9 (http://www.bioinformatics.babraham.ac.uk/projects/fastqc/ accessed on 5 November 2019) and Quast (http://cab.cc.spbu.ru/quast/ accessed on 5 November 2019). Contigs were further annotated with RAST (https://rast.nmpdr.org accessed on 6 November 2019). This Whole Genome Shotgun project has been deposited at DDBJ/ENA/GenBank under the BioProject accession number PRJNA484888, with Biosamples SAMN16522901 and SAMN16522902.
Provided with a snapshot of the ST11 clade distribution in previous papers [9,12], we decided to focus on understanding the phylogeny and geographic distribution of the ST11 KL105 clade since these isolates were poorly represented in those series. Twenty-eight ST11-KL105 K. pneumoniae genomes available on GenBank database (https://www.ncbi.nlm.nih.gov/genbank accessed on 18 November 2019) were included for comparison (Supplementary Table S1). Antibiotic resistance and replicon content were extracted from Center for Genomic Epidemiology tools (ResFinder, PlasmidFinder at http://www.genomicepidemiology.org accessed on 5 November 2019) or CARD (https://card.mcmaster.ca accessed on 5 November 2019). Plasmid contigs were identified using mlplasmids [34]. Plasmid assembly and partial reconstruction was obtained, mapping ST11 contigs on reference plasmid sequences and BLAST (https://blast.ncbi.nlm.nih.gov/Blast.cgi accessed on 5 November 2019). Geneious version 9.1.8 was used to represent the blaNDM-1 genetic environment.
The virulence gene content and wzi, capsular (KL) and O antigen loci were assessed using the K. pneumoniae Institut Pasteur database (https://bigsdb.pasteur.fr/cgi-bin/bigsdb/bigsdb.pl?db=pubmlst_klebsiella_seqdef accessed on 5 November 2019) and KLEBORATE (https://github.com/katholt/Kleborate accessed on 5 November 2019). The core genome of the two sequenced isolates (K606 and K624), the two previously described DHA-1-producing ST11-KL105 strains from Portugal [25] and those of 26 other ST11-KL105 strains deposited in the GenBank database were extracted and used to construct a SNP-based phylogenetic tree using the CSI phylogeny tool, using default parameters (https://cge.cbs.dtu.dk/services/CSIPhylogeny/ accessed on 6 November 2019) and the maximum-likelihood method to calculate bootstrap values. The closest genome corresponding to strain KLPN_19 (a ST11 strain with KL107 GCF_003861305.1) was used as a reference, and the tree was further represented with Interactive Tree of Life (iTOL, https://itol.embl.de accessed on 7 November 2019).

5. Conclusions

Data from the first recognized outbreak of NDM-1-producing K. pneumoniae in our country raise the alarm for the potential silent dissemination of NDM-1 through unnoticed carriage in the community, which needs to be monitored to prevent further spread. This report extends the distribution map of NDM-1 carbapenemase in Europe, which emerged through the re-introduction of ST11 KL105 K. pneumoniae clade into our country and a new plasmid backbone carrying blaNDM-1. We provide the first comparative analysis of the ST11 KL105 K. pneumoniae clade and show that it is successfully established in European countries carrying a variable pool of ESBL, DHA-1 or NDM-1-encoding plasmids. Finally, we further demonstrate the usefulness of the speed and the accuracy of FT-IR spectroscopy to support outbreak identification and infection control.

Supplementary Materials

The following are available https://www.mdpi.com/article/10.3390/antibiotics11010092/s1, Table S1: List of K. pneumoniae ST11-KL105 genomes deposited on NCBI public databases used for phylogenetic analysis. Figure S1: Projection of FT-IR spectra from the outbreak K. pneumoniae isolates (gray dots) in the in-house partial-least squares discriminant model (PLSDA) [13]. Figure S2: SNP matrix obtained from core-genome DNA comparison between K. pneumoniae ST11 KL105 genomes, supporting phylogenetic tree.

Author Contributions

Â.N. designed the study, acquired and interpreted molecular and whole-genome sequencing data and wrote the manuscript. R.V.F., M.V. and P.M.d.C. performed the characterization of isolates and collected epidemiological information. L.P. revised the manuscript for scientific content. All authors have read and agreed to the published version of the manuscript.

Funding

This work is financed by national funds from FCT—Fundação para a Ciência e a Tecnologia, I.P., in the scope of the project UIDP/04378/2020 and UIDB/04378/2020 of the Research Unit on Applied Molecular Biosciences—UCIBIO and the project LA/P/0140/2020 of the Associate Laboratory Institute for Health and Bioeconomy—i4HB. Additionally, we would like to acknowledge the support of HEALTH-UNORTE I&D&I project (NORTE-01-0145- FEDER-000039) co-financed by the European Regional Development Fund (ERDF), through the NORTE 2020 (Programa Operacional Regional do Norte 2014/2020). Ângela Novais is supported by national funds through FCT, I.P., in the context of the transitional norm (DL57/2016/CP1346/CT0032).

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

This Whole Genome Shotgun project has been deposited at DDBJ/ENA/GenBank under the BioProject accession number PRJNA484888, with Biosamples SAMN16522901 and SAMN16522902.

Acknowledgments

We thank Mariana Silva Nunes for the help in initial isolate processing and organization of information. We are also grateful to João Pedro Rueda Furlan for his support on Geneious.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Wyres, K.L.; Lam, M.M.C.; Holt, K.E. Population genomics of Klebsiella pneumoniae. Nat. Rev. Microbiol. 2020, 18, 344–359. [Google Scholar] [CrossRef] [PubMed]
  2. David, S.; Reuter, S.; Harris, S.R.; Glasner, C.; Feltwell, T.; Argimon, S.; Abudahab, K.; Goater, R.; Giani, T.; Errico, G.; et al. Epidemic of carbapenem-resistant Klebsiella pneumoniae in Europe is driven by nosocomial spread. Nat. Microbiol. 2019, 4, 1919–1929. [Google Scholar] [CrossRef]
  3. Lee, C.R.; Lee, J.H.; Park, K.S.; Kim, Y.B.; Jeong, B.C.; Lee, S.H. Global Dissemination of Carbapenemase-Producing Klebsiella pneumoniae: Epidemiology, Genetic Context, Treatment Options, and Detection Methods. Front. Microbiol. 2016, 13, 895. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  4. Guerra, A.M.; Lira, A.; Lameirão, A.; Selaru, A.; Abreu, G.; Lopes, P.; Mota, M.; Novais, Â.; Peixe, L. Multiplicity of carbapenemase-producers three years after a KPC-3-producing Klebsiella pneumoniae ST147-K64 hospital outbreak. Antibiotics 2020, 13, 806. [Google Scholar] [CrossRef]
  5. Bonnin, R.A.; Jousset, A.B.; Chiarelli, A.; Emeraud, C.; Glaser, P.; Naas, T.; Dortet, L. Emergence of New Non-Clonal Group 258 High-Risk Clones among Klebsiella pneumoniae Carbapenemase-Producing K. pneumoniae Isolates, France. Emerg. Infect. Dis. 2020, 26, 1212–1220. [Google Scholar] [CrossRef] [PubMed]
  6. Hernández-García, M.; Pérez-Viso, B.; Carmen Turrientes, M.; Díaz-Agero, C.; López-Fresneña, N.; Bonten, M.; Malhotra-Kumar, S.; Ruiz-Garbajosa, P.; Cantón, R. Characterization of carbapenemase-producing Enterobacteriaceae from colonized patients in a university hospital in Madrid, Spain, during the R-GNOSIS project depicts increased clonal diversity over time with maintenance of high-risk clones. J. Antimicrob. Chemother. 2018, 73, 3039–3043. [Google Scholar] [CrossRef]
  7. Bonura, C.; Giuffrè, M.; Aleo, A.; Fasciana, T.; Di Bernardo, F.; Stampone, T.; Giammanco, A.; The MDR-GN Working Group; Palma, D.M.; Mammina, C. An Update of the Evolving Epidemic of blaKPC Carrying Klebsiella pneumoniae in Sicily, Italy, 2014: Emergence of Multiple Non-ST258 Clones. PLoS ONE 2015, 10, e0132936. [Google Scholar] [CrossRef] [Green Version]
  8. Chen, L.; Mathema, B.; Pitout, J.D.; DeLeo, F.R.; Kreiswirth, B.N. Epidemic Klebsiella pneumoniae ST258 is a hybrid strain. mBio 2014, 5, e01355-14. [Google Scholar] [CrossRef] [Green Version]
  9. Dong, N.; Zhang, R.; Liu, L.; Li, R.; Lin, D.; Chan, E.W.; Chen, S. Genome analysis of clinical multilocus sequence Type 11 Klebsiella pneumoniae from China. Microb. Genom. 2018, 4, e000149. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  10. Wyres, K.L.; Hawkey, J.; Hetland, M.A.K.; Fostervold, A.; Wick, R.R.; Judd, L.M.; Hamidian, M.; Howden, B.P.; Löhr, I.H.; Holt, K.E. Emergence and rapid global dissemination of CTX-M-15-associated Klebsiella pneumoniae strain ST307. J. Antimicrob. Chemother. 2019, 74, 577–581. [Google Scholar] [CrossRef] [Green Version]
  11. Zhang, R.; Liu, L.; Zhou, H.; Chan, E.W.; Li, J.; Fang, Y.; Li, Y.; Liao, K.; Chen, S. Nationwide surveillance of clinical carbapenem-resistant Enterobacteriaceae (CRE) strains in China. EBioMedicine 2017, 19, 98–106. [Google Scholar] [CrossRef] [Green Version]
  12. Zhou, K.; Xiao, T.; David, S.; Wang, Q.; Zhou, Y.; Guo, L.; Aanensen, D.; Holt, K.E.; Thomson, N.T.; Grundmann, H.; et al. Novel Subclone of Carbapenem-Resistant Klebsiella pneumoniae Sequence Type 11 with Enhanced Virulence and Transmissibility, China. Emerg. Infect. Dis. 2020, 26, 289–297. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  13. Rodrigues, C.; Sousa, C.; Lopes, J.A.; Novais, Â.; Peixe, L. A Front Line on Klebsiella pneumoniae Capsular Polysaccharide Knowledge: Fourier Transform Infrared Spectroscopy as an Accurate and Fast Typing Tool. mSystems 2020, 5, e00386-19. [Google Scholar] [CrossRef] [Green Version]
  14. Silva, L.; Rodrigues, C.; Lira, A.; Leão, M.; Mota, M.; Lopes, P.; Novais, Â.; Peixe, L. Fourier transform infrared (FT-IR) spectroscopy typing: A real-time analysis of an outbreak by carbapenem-resistant Klebsiella pneumoniae. Eur. J. Clin. Microbiol. Infect. Dis. 2020, 39, 2471–2475. [Google Scholar] [CrossRef]
  15. Direção-Geral da Saúde. Infeções e Resistências aos Antimicrobianos: Relatório Anual do Programa Prioritário 2018. Available online: https://www.dgs.pt/portal-da-estatistica-da-saude/diretorio-de-informacao/diretorio-de-informacao/por-serie-1003038-pdf.aspx?v=%3d%3dDwAAAB%2bLCAAAAAAABAArySzItzVUy81MsTU1MDAFAHzFEfkPAAAA (accessed on 5 November 2019).
  16. Rodrigues, C.; Machado, E.; Ramos, H.; Peixe, L.; Novais, Â. Expansion of ESBL-producing Klebsiella pneumoniae in hospitalized patients: A successful story of international clones (ST15, ST147, ST336) and epidemic plasmids (IncR, IncFIIK). Int. J. Med. Microbiol. 2014, 304, 1100–11088. [Google Scholar] [CrossRef]
  17. Rodrigues, C.; Bavlovič, J.; Machado, E.; Amorim, J.; Peixe, L.; Novais, Â. KPC-3-Producing Klebsiella pneumoniae in Portugal Linked to Previously Circulating Non-CG258 Lineages and Uncommon Genetic Platforms (Tn4401d-IncFIA and Tn4401d-IncN). Front. Microbiol. 2016, 7, 1000. [Google Scholar] [CrossRef] [Green Version]
  18. Aires-de-Sousa, M.; Ortiz de la Rosa, J.M.; Gonçalves, M.L.; Pereira, A.L.; Nordmann, P.; Poirel, L. Epidemiology of Carbapenemase-Producing Klebsiella pneumoniae in a Hospital, Portugal. Emerg. Infect. Dis. 2019, 25, 1632–1638. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  19. Perdigão, J.; Modesto, A.; Pereira, A.L.; Neto, O.; Matos, V.; Godinho, A.; Phelan, J.; Charleston, J.; Sadar, A.; de Sessions, P.F.; et al. Whole-genome sequencing resolves a polyclonal outbreak by extended-spectrum beta-lactam and carbapenem-resistant Klebsiella pneumoniae in a Portuguese tertiary-care hospital. Microb. Genom. 2019, 7, 000349. [Google Scholar] [CrossRef]
  20. Lopes, E.; Saavedra, M.J.; Costa, E.; de Lencastre, H.; Poirel, L.; Aires-de-Sousa, M. Epidemiology of carbapenemase-producing Klebsiella pneumoniae in northern Portugal: Predominance of KPC-2 and OXA-48. J. Glob. Antimicrob. Resist. 2020, 22, 349–353. [Google Scholar] [CrossRef] [PubMed]
  21. Aires-de-Sousa, M.; Ortiz de la Rosa, J.M.; Goncalves, M.L.; Costa, A.; Nordmann, P.; Poirel, L. Occurrence of NDM-1-producing Morganella morganii and Proteus mirabilis in a single patient in Portugal: Probable in vivo transfer by conjugation. J. Antimicrob. Chemother. 2020, 75, 903–906. [Google Scholar] [CrossRef]
  22. Manageiro, V.; Sampaio, D.A.; Pereira, P.; Rodrigues, P.; Veira, L.; Palos, C.; Caniça, M. Draft Genome Sequence of the First NDM-1-Producing Providencia stuartii Strain Isolated in Portugal. Genome Announc. 2015, 3, e01077-15. [Google Scholar] [CrossRef] [Green Version]
  23. Wailan, A.M.; Paterson, D.L. The spread and acquisition of NDM-1: A multifactorial problem. Expert. Rev. Anti. Infect. Ther. 2014, 12, 91–115. [Google Scholar] [CrossRef] [PubMed]
  24. Politi, L.; Gartzonika, K.; Spanakis, N.; Zarkotou, O.; Poulou, A.; Skoura, L.; Vrioni, G.; Tsakris, A. Emergence of NDM-1-producing Klebsiella pneumoniae in Greece: Evidence of a widespread clonal outbreak. J. Antimicrob. Chemother. 2019, 74, 2197–2202. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  25. Ribeiro, T.G.; Novais, Â.; Rodrigues, C.; Nascimento, R.; Freitas, F.; Machado, E.; Peixe, L. Dynamics of clonal and plasmid backgrounds of Enterobacteriaceae producing acquired AmpC in Portuguese clinical settings over time. Int. J. Antimicrob. Agents 2019, 53, 650–656. [Google Scholar] [CrossRef]
  26. Izdebski, R.; Sitkiewicz, M.; Urbanowicz, P.; Krawczyk, M.; Gniadkowski, M. Genomic background of the Klebsiella pneumoniae NDM-1 outbreak in Poland, 2012-18. J. Antimicrob. Chemother. 2020, 75, 3156–3162. [Google Scholar] [CrossRef] [PubMed]
  27. van Dorp, L.; Wang, Q.; Shaw, L.P.; Acman, M.; Brynildsrud, O.B.; Eldholm, V.; Wang, R.; Gao, H.; Yin, Y.; Chen, H.; et al. Rapid phenotypic evolution in multidrug-resistant Klebsiella pneumoniae hospital outbreak strains. Microb. Genom. 2019, 5, e000263. [Google Scholar] [CrossRef]
  28. Bowers, J.R.; Kitchel, B.; Driebe, E.M.; MacCannell, D.R.; Roe, C.; Lemmer, D.; de Man, T.; Rasheed, J.K.; Engelthaler, D.M.; Keim, P.; et al. Genomic Analysis of the Emergence and Rapid Global Dissemination of the Clonal Group 258 Klebsiella pneumoniae Pandemic. PLoS ONE 2015, 10, e0133727. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  29. Novais, Â.; Freitas, A.R.; Rodrigues, C.; Peixe, L. Fourier transform infrared spectroscopy: Unlocking fundamentals and prospects for bacterial strain typing. Eur. J. Clin. Microbiol. Infect. Dis. 2019, 38, 427–448. [Google Scholar] [CrossRef]
  30. Naas, T.; Oueslati, S.; Bonnin, R.A.; Dabos, M.L.; Zavala, A.; Dortet, L.; Retailleau, P.; Iorga, B.I. Beta-lactamase database (BLDB)—structure and function. J. Enzyme Inhib. Med. Chem. 2017, 32, 917–919. [Google Scholar] [CrossRef]
  31. Pires, J.; Novais, A.; Peixe, L. Blue-carba, an easy biochemical test for detection of diverse carbapenemase producers directly from bacterial cultures. J. Clin. Microbiol. 2013, 51, 4281–4283. [Google Scholar] [CrossRef] [Green Version]
  32. Bogaerts, P.; Rezende de Castro, R.; de Mendonça, R.; Huang, T.D.; Denis, O.; Glupczynski, Y. Validation of carbapebemase and extended-spectrum β-lactamase multiplex endpoint PCR assays according to ISO 15189. J. Antimicrob. Chemother. 2013, 68, 1576–1582. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  33. Brisse, S.; Passet, V.; Haugaard, A.B.; Babosan, A.; Kassis-Chikhani, N.; Struve, C.; Decré, D. wzi gene sequencing, a rapid method to determine the capsular type of Klebsiella strains. J. Clin. Microbiol. 2013, 51, 4073–4078. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  34. Arredondo-Alonso, S.; Rogers, M.R.C.; Braat, J.C.; Verschuuren, T.D.; Top, J.; Corander, J.; Willems, R.J.L.; Schürch, A.C. mlplasmids: A user-friendly tool to predict plasmid- and chromosome-derived sequences for single species. Microb. Genom. 2018, 4, e000224. [Google Scholar] [CrossRef] [PubMed] [Green Version]
Antibiotics 11 00092 g001 550
Figure 1. SNP-based phylogenetic tree of ST11-KL105 genomes and heatmap of antimicrobial resistance genotype. Genomes sequenced from this study (shown in bold) were compared with those from 28 other ST11-KL105 strains deposited in public databases (accessed on 5 November 2019). Phylogenetic tree was obtained using the SNP method of CSI phylogeny with default parameters and the reference strain KLPN_19, the closest ST11 with KL107 (accession number GCF_003861305.1). Country and year of isolation are shown, as represented in metadata, as well as their content of acquired beta-lactamases (ESBL, acquired AmpC and carbapenemases), other antibiotic resistance genes and replicons. Numbers in the left indicate bootstrap values >0.85, calculated based on the maximum-likelihood method.
Figure 1. SNP-based phylogenetic tree of ST11-KL105 genomes and heatmap of antimicrobial resistance genotype. Genomes sequenced from this study (shown in bold) were compared with those from 28 other ST11-KL105 strains deposited in public databases (accessed on 5 November 2019). Phylogenetic tree was obtained using the SNP method of CSI phylogeny with default parameters and the reference strain KLPN_19, the closest ST11 with KL107 (accession number GCF_003861305.1). Country and year of isolation are shown, as represented in metadata, as well as their content of acquired beta-lactamases (ESBL, acquired AmpC and carbapenemases), other antibiotic resistance genes and replicons. Numbers in the left indicate bootstrap values >0.85, calculated based on the maximum-likelihood method.
Antibiotics 11 00092 g001
Antibiotics 11 00092 g002 550
Figure 2. Genetic environment of blaNDM-1 in plasmid pAN_K624, represented using Geneious version 9.1.8. The location of the insertion of the MRR region into the pAR_0109-like plasmid backbone is shown.
Figure 2. Genetic environment of blaNDM-1 in plasmid pAN_K624, represented using Geneious version 9.1.8. The location of the insertion of the MRR region into the pAR_0109-like plasmid backbone is shown.
Antibiotics 11 00092 g002
Table
Table 1. Epidemiological and clinical data of patients colonized/infected by NDM-1-producing K. pneumoniae ST11-KL105.
Table 1. Epidemiological and clinical data of patients colonized/infected by NDM-1-producing K. pneumoniae ST11-KL105.
Patient
N° 		Sex/Age Isolates Date of Isolation Hospital
Unit 		Period of Hospitalization Patient’s Origin Sample Type Pathology Colonization at Admission Antibiotherapy (Previous 3 Months) Previous Hospitalization (Date, Unit)
1 M/71 K607
K606 *		27/11/18
28/11/18 		Surgery 26–30 October 2018 Hospital C (bile)
C (RS) 		Cholangiocarcinoma Not tested Piperacillin-Tazobactam None
2 M/79 K608 28/11/18 Surgery 15 November–3 December 2018 Hospital C (RS) Acute pancreatitis Not tested None None
3 M/55 K609 30/11/18 Surgery 24 October–5 December 2018 Hospital C (RS) Trauma - Amoxacillin-clavulanic acid 15 September–22 October 2018 (different hospital)
4 M/66 K610 30/11/18 Surgery 23–30 November 2018 Hospital C (RS) Adenocarcinoma Not tested None None
5 M/66 K611 14/12/19 Emergency 14 December 2018–24 January 2019 Outpatient C (RS) Colon carcinoma + Ciprofloxacin 10–28 November 2018 (Surgery)
6 M/68 K612 25/12/19 Emergency 25 December 2018–11 January 2019 (Medicine) Outpatient C (RS) Soft tissue infection + Tigecycline,
Ceftazidime 		31 October–9 December 2018 (Surgery)
7 M/87 K624 *
K625 		25/03/19 Emergency 25 March–18 April 2019 Household I (urine)
C (RS) 		Pneumonia - Cefuroxime 28 January–4 March 2019
* Isolates selected for whole-genome sequencing. M = male; C = colonization; I = infection; RS = rectal swab.
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Novais, Â.; Ferraz, R.V.; Viana, M.; da Costa, P.M.; Peixe, L. NDM-1 Introduction in Portugal through a ST11 KL105 Klebsiella pneumoniae Widespread in Europe. Antibiotics 2022, 11, 92. https://doi.org/10.3390/antibiotics11010092

AMA Style

Novais Â, Ferraz RV, Viana M, da Costa PM, Peixe L. NDM-1 Introduction in Portugal through a ST11 KL105 Klebsiella pneumoniae Widespread in Europe. Antibiotics. 2022; 11(1):92. https://doi.org/10.3390/antibiotics11010092

Chicago/Turabian Style

Novais, Ângela, Rita Veiga Ferraz, Mariana Viana, Paula Martins da Costa, and Luísa Peixe. 2022. "NDM-1 Introduction in Portugal through a ST11 KL105 Klebsiella pneumoniae Widespread in Europe" Antibiotics 11, no. 1: 92. https://doi.org/10.3390/antibiotics11010092

APA Style

Novais, Â., Ferraz, R. V., Viana, M., da Costa, P. M., & Peixe, L. (2022). NDM-1 Introduction in Portugal through a ST11 KL105 Klebsiella pneumoniae Widespread in Europe. Antibiotics, 11(1), 92. https://doi.org/10.3390/antibiotics11010092

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