Characterizing carbapenemase-producing Escherichia coli isolates from Spain: high genetic heterogeneity and wide geographical spread

Introduction: Carbapenemase-Producing Escherichia coli (CP-Eco) isolates, though less prevalent than other CP-Enterobacterales, have the capacity to rapidly disseminate antibiotic resistance genes (ARGs) and cause serious difficult-to-treat infections. The aim of this study is phenotypically and genotypically characterizing CP-Eco isolates collected from Spain to better understand their resistance mechanisms and population structure.

Methods: Ninety representative isolates received from 2015 to 2020 from 25 provinces and 59 hospitals Spanish hospitals were included. Antibiotic susceptibility was determined according to EUCAST guidelines and whole-genome sequencing was performed. Antibiotic resistance and virulence-associated genes, phylogeny and population structure, and carbapenemase genes-carrying plasmids were analyzed.

Results and discussion: The 90 CP-Eco isolates were highly polyclonal, where the most prevalent was ST131, detected in 14 (15.6%) of the isolates. The carbapenemase genes detected were bla_OXA-48 (45.6%), bla_VIM-1 (23.3%), bla_NDM-1 (7.8%), bla_KPC-3 (6.7%), and bla_NDM-5 (6.7%). Forty (44.4%) were resistant to 6 or more antibiotic groups and the most active antibiotics were colistin (98.9%), plazomicin (92.2%) and cefiderocol (92.2%). Four of the seven cefiderocol-resistant isolates belonged to ST167 and six harbored bla_NDM. Five of the plazomicin-resistant isolates harbored rmt. IncL plasmids were the most frequent (45.7%) and eight of these harbored bla_VIM-1. bla_OXA-48 was found in IncF plasmids in eight isolates. Metallo-β-lactamases were more frequent in isolates with resistance to six or more antibiotic groups, with their genes often present on the same plasmid/integron. ST131 isolates were associated with sat and pap virulence genes. This study highlights the genetic versatility of CP-Eco and its potential to disseminate ARGs and cause community and nosocomial infections.

1 Introduction

Antibiotic-resistant Escherichia coli is one of the principal pathogens that can cause community-acquired and nosocomial infections with increased morbidity and mortality rates (Hu et al., 2022). They are considered as one of the principal causes of Urinary Tract Infections (UTIs), bacteremia, and Intra-Abdominal Infections (IAI) (Balasubramanian et al., 2023). E. coli have a heightened ability to acquire Antibiotic-Resistance Genes (ARGs), such as the bla_CTX-M-15 Extended Spectrum β-Lactamase (ESBL), and rapidly disseminate them throughout the community (González et al., 2020). Carbapenemase-producing Escherichia coli (CP-Eco), though not as frequently isolated in the clinical setting as compared to other CP-producing Enterobacterales (such as Klebsiella pneumoniae and Enterobacter cloacae complex), is especially worrying. This is because they are increasing in prevalence (Cañada-García et al., 2022), creating a concern that they would be able to spread carbapenemase genes in the community in a similar manner as it was observed for ESBLs (González et al., 2020). Moreover, these isolates are typically resistant to several other antibiotics, making their related infections difficult to treat (Boutzoukas et al., 2023).

All main carbapenemase families have been detected in CP-Eco (Grundmann et al., 2017), in addition to a wide range of virulence determinants that negatively affect clinical outcome (Čurová et al., 2020). All this led to the declaration by the World Health Organization that CP-Eco are a critical priority issue (Tacconelli et al., 2018). Globally, antibiotic-resistant E. coli has the highest incidence of hospital-acquired infections in high-to-mid income countries, causing three to twenty-five million infections per year (Balasubramanian et al., 2023). In Europe, the median number of infections caused by CP-Eco was 2,619 in 2015, with a median of 141 attributable deaths (Cassini et al., 2019). In Spain, the incidence of CP-Eco has evolved from isolated cases in 2013 (Oteo et al., 2015), to being detected in 10 different Spanish provinces in 2019 (Cañada-García et al., 2022).

Given the importance and the widespread dissemination of CP-Eco, the aim of this study is the phenotypic and genotypic characterization of a set of such isolates obtained from different Spanish regions over a period of five years.

2 Materials and methods

2.1 Study design and bacterial collection

As part of the ongoing Antibiotic Resistance Surveillance Program, the Spanish National Center for Microbiology (CNM) receives clinical carbapenem-resistant E. coli isolates for molecular characterization from Spanish hospitals. Participation in this program is voluntary and the criteria to send the isolates to the CNM is having a Minimum Inhibitory Concentration (MIC) to meropenem of over 0.12 mg/L (Giske et al., 2017).

Of the total CP-Eco isolates received in this Surveillance Program between January 2015 and December 2020, a subset was selected for further analysis according to the following criteria: i) one isolate for each type of carbapenemase per year was included from each of the collaborating hospitals, ii) in the event that there was more than one isolate for each carbapenemase type in the same year, the selection was prioritized based on the type of infection in which they were involved according to the following order: 1) invasive infections 2) other infections 3) colonizations; iii) in all cases, the first isolate received per year was selected. These selection criteria allowed for a temporal-spatial representative sample of the different types of carbapenemases circulating in E. coli in Spain. Not all hospitals had isolates from all the five years, and from some hospitals more than one isolate per year was included since they produced different carbapenemase types. The information sent along with the isolate from the participating hospitals included collection date, age and sex of the patients, sample type, whether the isolate was a colonizer or causing infections (as determined by the attending microbiologist), and isolate origin [hospital, healthcare center, or community setting; i.e. where the infection is contracted without interaction with healthcare settings during the last 48 hours (Plachouras et al., 2018)].

2.2 DNA extraction and identification of carbapenemase genes

All the carbapenem-resistant E. coli isolates received at the CNM from 2015 to 2020 were cultured on MacConkey agar (Becton Dickinson Microbiology Systems, Cockeysville, MD, USA) and incubated overnight at 37°C. DNA was extracted using the QIAamp^® DNA Mini Kit (QIAGEN, Hilden, Germany) according to the manufacturer’s instructions. Real-time PCR assays targeting the carbapenemase gene families bla_OXA-48, bla_KPC, bla_VIM, bla_IMP and bla_NDM were performed (Ortega et al., 2016). All the isolates were stored at -80°C until used.

2.3 Antibiotic susceptibility testing

Antibiotic Susceptibility Testing (AST) was determined for the selected isolates with broth microdilution using the YDKMGN Sensititre™ Gram Negative panels (Thermo Fisher, Waltham, MA, USA) (International Organization for Standardization (ISO), 2006). E-tests using meropenem/vaborbactam and imipenem/relebactam and disc diffusion assays for cefiderocol (Liofilchem, Roseto degli Abruzzi, Italy) were performed and interpreted according to the EUCAST guidelines (European Society of Clinical Microbiology and Infectious Diseases (EUCAST), 2022). Additionally, E-tests using plazomicin (Liofilchem, Roseto degli Abruzzi, Italy) and fosfomycin (bioMérieux, Marcy l’Etoile, France) strips were performed and interpreted according to an FDA-approved breakpoint of susceptibility of ≤2mg/L for plazomycin (Cañada-García et al., 2022), and according to the EUCAST cutoff values for urinary tract infections (UTI) for fosfomycin due to lack of clinical breakpoints for other types of samples (European Society of Clinical Microbiology and Infectious Diseases (EUCAST), 2022), respectively.

Resistance to different antibiotic families was analyzed considering 14 categories: nine categories according to the classification proposed by Magiorakos et al. (Magiorakos et al., 2012) (extended spectrum cephalosporins, monobactams, carbapenems, aminoglycosides, fluoroquinolones, fosphonic acid derivatives, folic pathway inhibitors, glycylcyclines and polymyxins); and 5 additional categories: ceftazidime/avibactam (extended spectrum cephalosporin + diazabicycloctane inhibitor), imipenem/relebactam (carbapenem + diazabicycloctane inhibitor), meropenem/vaborbactam (carbapenem + boronic acid inhibitor), cefiderocol (siderophore cephalosporin) and plazomicin (extended spectrum aminoglycoside).

2.4 Whole-genome sequencing and read assembly

DNA was extracted from the selected isolates as described in section 2.2. Paired-end (2x150) libraries were then prepared using the Nextera DNA Flex Library Preparation Kit and sequenced using Illumina NextSeq 550 (Illumina Inc., San Diego, CA, United States) according to the manufacturer’s instructions. Quality of the reads was assessed using FASTQC (version 0.11.9), followed by de novo assembly using Unicycler (version 0.4.8) (Wick et al., 2017). The quality of the assembly was assessed using QUAST (version 5.2.0) and annotation was performed using Prokka (version 1.14.6) (Seemann, 2014). The obtained assemblies were deposited in the European Nucleotide Archive under the accession number PRJEB70795.

2.5 Phylogenetic analyses and diversity

Ridom SeqSphere+ (version 8.3.1; Ridom, Münsten, Germany) was used to perform a core-genome Multi-Locus Sequence-Typing analysis (cgMLST) using a built-in scheme for E. coli containing 2,515 core genes, and to construct a minimum spanning tree based on allelic differences. ARIBA (version 2.6.2) (Hunt et al., 2017) was used to determine STs in accordance with the University of Warwick scheme. Diversity of the samples was calculated using a Simple Diversity Index [SDI (Gastmeier et al., 2006)]. Analysis of the fimH gene was done to subclassify ST131isolates using the FimTyper tool [CGE server; https://cge.cbs.dtu.dk; (accessed 23/02/2024)].

2.6 Antibiotic resistance genes, virulence-associated genes, and plasmids

ARGs were analyzed by ARIBA (version 2.14.6) (Hunt et al., 2017) using the CARD database (https://card.mcmaster.ca; (accessed 23/02/2024)) and ResFinder (CGE server) with ID thresholds of 100% for β-lactamase variants and 98% for other resistance genes. In isolates resistant to cefiderocol, the sequences of the genes encoding PBP3 and the siderophore CirA, previously related to resistance to this antibiotic (Wang et al., 2022), were analyzed with ARIBA.

The VirulenceFinder tool (CGE server) (accessed 23/02/2024) was used to detect virulence-associated genes (five toxin-coding genes, four adhesin-coding genes, two siderophore-coding genes, one outer membrane protein-coding gene and one invasive protein-coding gene). PlasmidID [https://github.com/BU-ISCIII/plasmidID; (accessed 23/02/2024)] was used to map the reads against a curated plasmid database, perform de novo plasmid assemblies, and determine the presence of resistance and replicon genes (Pérez-Vazquez et al., 2019).

2.7 Statistical analyses

The differences in prevalence of the ARGs, virulence-related genes, STs, and sample characteristics were evaluated using the Fisher’s exact test using GraphPad Prism (version 3.02; GraphPad Software, Inc., San Diego, CA, USA). P values of less than 0.05 were considered statistically significant.

2.8 Data availability

All the data generated for this study is available in the manuscript and its related Supplementary Material. The sequenced raw reads have been deposited in the open-access European Nucleotide Archive database under the accession number PRJEB70795.

3 Results

3.1 Characteristics of all the CP-Eco isolates received from 2015 to 2020

From January 2015 until December 2020, a total of 593 CP-E. coli isolates were confirmed as carbapenems producers by the Spanish CNM Reference Lab after screening for MIC to meropenem of over 0.12 mg/L by the collaborating hospitals. Supplementary Table S1 shows the characteristics and origin of the isolates, as well as the information received from the participating hospitals (total and per year). The isolates were received from 25 Spanish provinces and 59 hospitals. Two-hundred-ninety-nine (50.4%) of the isolates were from females and the average patient age was 77.8 years. Three-hundred-one (50.8%) isolates were colonizers, and the rest were causing infections (24 were unspecified). The most common infection was UTI (175; 29.5%), followed by wound/ulcer infections (32; 5.4%), bacteremia (25; 4.2%), and IAI (17; 2.9%). Four-hundred-forty-three (74.7%) were healthcare-associated cases, 391 from hospitals and 52 from other healthcare centers, and 126 (21.2%) were from community settings; 24 were unspecified.

Using real-time PCR carbapenemase genes were detected in all 593 isolates where bla_OXA-48-like was detected in 485 (81.8%), bla_VIM-like in 71 (12%), bla_NDM-like in 17 (2.9%), and bla_KPC-like in 17 (2.9%) isolates. Three isolates harbored two carbapenemases, one of each of the following combinations bla_OXA-48 and bla_KPC, bla_OXA-48 and bla_NDM, and bla_VIM and bla_KPC. In the community setting, the most frequent carbapenemase was bla_{OXA-48-family}, detected in 109 (86.5%) isolates. The relative frequency of these genes per year was maintained throughout the study period and there was no statistically relevant correlation between the types of carbapenemase gene, sample type, location of isolation, infection, colonization, and type of carbapenemase harbored by the isolates (P > 0.05).

3.2 Characteristics and population structure of the 90 selected CP-Eco isolates

Based on the criteria set in section 2.1, 90 isolates from 59 hospitals across 25 Spanish provinces were selected. Supplementary Table S2 lists all the characteristics of these isolates. Sixty (66.6%) isolates were implicated in clinical infections, as determined by the doctor who treated the patient, and mainly were obtained from urine (27; 45%), blood (15; 25%) and abscesses and wounds (12; 20%). The remaining 30 (33.3%) isolates were considered colonizing and came mainly from rectal swabs and feces (26; 86.7%)

Forty-two (46.7%) of the isolates were obtained from the hospital setting, the majority of which harbored bla_OXA-48 (24; 57.1%), followed by bla_VIM-1 (8; 19%), bla_NDM (3 bla_NDM-1 and 3 bla_NDM-5; 14.3% in total), bla_KPC (1 bla_KPC-2 and 1 bla_KPC-3; 4.8% in total), bla_OXA-162 (1; 2.4%), and both bla_KPC-3 and bla_OXA-48 (1; 2.4%). Twelve isolates were obtained from healthcare centers where 3 harbored bla_OXA-48, 3 bla_VIM-1, 2 bla_NDM (1 bla_NDM-1 and 1 bla_NDM-7), 2 bla_KPC (1 bla_KPC-2 and 1 bla_KPC-3), 1 (8.3%) bla_OXA-162, and 1 both bla_KPC-3 and bla_OXA-48. From the community setting, 18 (20%) of the isolates were collected where 8 had bla_OXA-48, 4 bla_VIM-1, 3 bla_NDM (1 bla_NDM-1 and 2 bla_NDM-5), and 3 bla_KPC (2 bla_KPC-2 and 1 bla_KPC-3). The origin of 18 isolates could not be determined.

The carbapenemase genes detected among the infecting isolates were bla_OXA-48 (30; 50%), bla_VIM-1 (13; 21.7%), bla_NDM (5 bla_NDM-1 and 3 bla_NDM-5; 13.3% in total), bla_KPC (3 bla_KPC-2 and 3 bla_KPC-3; 10% in total), one bla_OXA-162 (1; 1.7%), and both bla_KPC and bla_OXA-48 (2; 3.3%). Among the colonizing isolates, the carbapenemases genes observed were bla_OXA-48 (11; 36.7%), bla_VIM-1 (8; 26.7%), bla_NDM (6; 20%) (3 bla_NDM-5, 2 bla_NDM-1, and 1 bla_NDM-7), bla_KPC (4; 13.3%) (3 bla_KPC-2 and 1 bla_KPC-3), and bla_OXA-162 (1; 3.3%). No statistically significant associations were found between the sample type, location of isolation, infection, colonization, and type of carbapenemase harbored by the isolates (P > 0.05)

MLST analyses showed that the 90 isolates belonged to 49 different STs (SDI=54.4) and only ST131 (14 isolates) and ST88 (6 isolates) were detected in more than five isolates. Population diversity was similar for isolates harboring bla_OXA-48 (SDI=65.1) and bla_NDM (SDI=64.3), and was lower than those harboring bla_KPC (SDI=71.0) and bla_VIM (SDI=81.0). Ridom SeqSphere+ was used to construct a maximum-likelihood phylogenetic tree using the gene-by-gene approach and allelic distance from cgMLST. Figure 1 shows a great genetic diversity of the CP-Eco studied alone and in relation to the types of carbapenemases they had. The average allelic distances in pairwise comparison of all the CP-Eco isolates was 2,120 ± 555, and only four main clusters with ≥ 4 isolates were identified. The first cluster included 14 (15.6%) CP-Eco isolates that belonged to ST131, had an average of 105 allelic differences (range:13-660), were collected from 7 provinces, and all 4 carbapenemase gene families were detected in them (including two with both bla_OXA-48 and bla_KPC-like). The second-largest cluster was formed by 6 (6.7%) isolates belonging to ST88, had an average of 34 allelic differences (range:20-58), were collected from 5 provinces and 5 isolates harbored bla_OXA-48. Cluster 3 contained 4 (4.4%) isolates belonging to ST167 that all harbored bla_NDM, were collected from 3 provinces and had an average of 142 allelic differences (range:79-264). Cluster 4 contained 4 (4.4%) isolates pertaining to ST1193 collected from 2 provinces with an average allelic difference of 24 (range:20-32); two of them harbored bla_OXA-48 and the other two bla_NDM. The rest of the isolates were grouped into STs with a frequency of 3 or less isolates. Applying a ≤ 5 allele differences threshold to identify possible transmission events, only one cluster of two OXA-48-producing ST127 isolates were detected.

Figure 1

Figure 1 Minimum spanning tree showing the population structure of the 90 isolates analyzed through core-genome multi-locus sequence typing (scheme of 2,515 genes). The distance between the genes are not up to scale, black lines represent differences higher than 50% of the alleles considered in the scheme (>1257 alleles), blue lines differences between 10-50% (252 to 1257) and red lines differences lower than 10% (<252). Sequence types with four or more isolates are circled and amplified. The blue shading highlights the only group (containing two ST127 isolates) that have a genetic distance of ≤ 5 alleles (this group has also been amplified).

Further typing of the fimH gene was conducted on the predominant ST131. Among the 14 ST131 isolates, 7 (50%) had fimH30 indicative of belonging to clade C, 4 (28.6%) had fimH41 (clade A), 1 (7.1%) had fimH22 (clade B), 1 (7.1%) had fimH298 (subtype of clade B), and 1 (7.1%) had fimH1196 (no specific clade). bla_CTX-M-15 and bla_KPC genes were more frequent in ST131 clade C isolates (71.4% and 57.1%, respectively) than in non-clade C isolates (14.3% and 28.6%, respectively) but these differences were not statistically significant (P > 0.05).

3.3 Characterization of the acquired antibiotic resistance genes

Table 1 shows the detected ARGs where an average of 5.6 ARGs were detected (range:1–13; excluding chromosomal genes and mutations). bla_OXA-48-like was detected in 43 (47.8%) isolates, bla_VIM-1 in 21 (23.3%), bla_NDM-like in 14 (15.6%), bla_KPC-like in 10 (11.1%), and both bla_OXA-48 and bla_KPC-like in 2 (2.2%) isolates. Isolates harboring bla_VIM-1 (7.2; range:3-11) and bla_NDM-like (6.7; range:1-12) had more ARGs than those isolates harboring bla_OXA-48 (5; range:1-13) and bla_KPC-like (4.6; range:1-12). Twelve isolates harbored only one ARG, six of them had bla_OXA-48-like, five had bla_KPC-like, and one had bla_VIM-1. Isolates belonging to ST88 had an average of 9.8 ARGs (range:7-13). bla_OXA-48 was most frequently encountered in ST88 isolates (5; 11.6%), bla_NDM-like in ST167 (4; 28.6%), bla_KPC-like in ST131 (6; 54.6%) and bla_VIM-like in ST131 isolates (3; 14.3%). bla_OXA-162 was detected in 2 isolates belonging to ST5968.

Table 1

Table 1 Genes and mutations associated with antibiotic resistance detected through whole-genome sequencing.

Acquired ESBL genes were detected in 27 (30%) isolates where bla_CTX-M-15 was predominant (n=13; 48.1%). bla_PER-2 was detected in a bla_NDM-1-harbouring isolate belonging to ST12. Regarding acquired AmpC genes, bla_CMY was detected in four isolates, three of them belonged to ST167 and all four harbored bla_NDM-like (Supplementary Table S2).

The non-β-lactamase acquired ARGs present in at least 25% of the CP-Eco isolates were the sulfonamide-resistance-encoding sul (n=64, 71.1%; mainly sul1); trimethoprim-resistance-encoding dfrA (n=56, 62.2%; mainly dfrA17); tetracycline-resistance-encoding tet (n=38, 42.2%; mainly tetA); aminoglycoside-resistance-encoding aac(6’) (n=31, 34.4%; mainly aac(6’)-Ib) and aadA (n = 52, 57.8%; mainly aadA1); chloramphenicol-resistance-encoding cat (n=25, 27.8%; mainly catA1); and quinolone-resistance-encoding qnr (n=23, 25.6%; mainly qnrS1) (Table 1; Supplementary Table S2). Twenty-five isolates had chromosomal mutations associated with resistance to fosfomycin and 2 additional isolates harbored the fosA7.5 and fosA3 acquired genes. These genes were located in chromosomal and plasmid contigs, respectively. The aminoglycoside methyltransferase gene rmt was detected in five isolates belonging to different STs, where 4 harbored bla_NDM-1 and 1 bla_NDM-5.

3.4 Characterization of the plasmids carrying the carbapenem resistance genes

The median size of the detected plasmids was 63,598 bp with a median coverage of 98% and a median similarity of 100% with the respective entries in Genbank. IncL plasmids were the most frequently detected (n=42; 46.7%). Thirty-two of these harboured bla_OXA-48, eight harboured bla_VIM-1, and two harboured bla_OXA-162. The second most frequent plasmid type was IncFII detected in 12 (13.3%) isolates, eight of which harbored bla_OXA-48 and four harbored bla_NDM. Five of the IncF-harboring isolates also belonged to clade C of ST131. IncN was detected in five (5.3%) isolates where four harbored bla_VIM-1 and one harbored bla_KPC-2. IncP6 was detected in four isolates harboring bla_KPC-2, and IncFIB was detected in four of the five isolates harboring bla_KPC-3. Other plasmid types were detected with a frequency of two or less and it was not possible to fully reconstruct the carbapenemase-harboring plasmids in 20 isolates (Supplementary Table S2).

No other ARG was detected on the same plasmid harboring bla_OXA-48-like nor for those harboring bla_KPC-like, except for one isolate harboring bla_KPC-2 that also harbored a class 1 integron with the sequence Int1-aac(6’)-Ib-arr3-dfrA27-aadA16-qacEΔ1-sul1. bla_VIM-1 in IncL plasmids was found in a class 1 integron with the sequence Int1-bla_VIM-1-aac(6’)-Ib-dfrB1-aadA1-catB2-qacEΔ1-sul1. In four other isolates, IncN plasmids harbored class 1 integron with bla_VIM-1 but with a different combination of genes. Six isolates harbored bla_NDM-5 (n=4), bla_NDM-1 (n=1) or bla_NDM-7 (n=1) in complex class 1 integrons with different combinations of aminoglycoside-resistance genes, including rmtC. Another isolate had bla_NDM-1 and rmtC on the same plasmid, but they did not form part of an integron (Supplementary Table S2).

3.5 Antibiotic resistance profiles

Table 2 shows the AST results that were analyzed according to the categories described in section 2.3. All isolates had MICs>0.12mg/L to meropenem, harbored at least one gene of carbapenem resistance, and were resistant to at least one carbapenem (Table 2). The two isolates producing two carbapenemases presented relatively low MICs of 4 mg/L to imipenem (both) and 4 and 1 mg/L to meropenem (each). Antibiotics with in-vitro activity against more that 90% of the isolates were colistin, plazomicin, cefiderocol, and meropenem/vaborbactam. The mcr genes and known mutations in the PhoPQ and PmrAB encoding genes related to colistin resistance were not detected in the colistin-resistant isolate.

Table 2

Table 2 Antibiotic susceptibility profiles of the 90 selected isolates determined through broth microdilution and interpreted according to the EUCAST guidelines.

Forty isolates (44.4%) were resistant to more than six out of the 14 antibiotic categories considered. These isolates harbored bla_VIM-1 (n=15; 37.5%), bla_NDM-like (n=13; 32.5%), bla_OXA-48-like (n=7; 17.5%), bla_KPC-like (n=4; 10%), and both bla_KPC-3 and bla_OXA-48 (n=1; 2.4%). Nine isolates (7 NDM-producers, 1 VIM-producer, and 1 OXA-48-producer) were only susceptible to three or four antibiotic categories. The most active agents against these isolates were colistin (9/9) and tigecycline (6/9).

Table 3 shows the antibiotic resistance rates of the CP-Eco isolates, grouped by the type of carbapenemase gene harbored. Isolates harboring bla_OXA-48 had lower resistance rates towards meropenem and imipenem as compared to non-bla_OXA-48-harboring isolates (4.7% and 7% versus 21.3% and 46.2%, respectively; P = 0.0003). OXA-48 and ESBL co-producers (n=15) were significantly more resistant to ceftazidime (46.7%), cefotaxime (100%), and aztreonam (100%) than OXA-48 producers without ESBL (n=26, with 3.8%, 15.4%, and 3.8% resistance rates, respectively (P<0.05)). No other β-lactamases combinations could be statistically studied in relation to antibiotic resistance profiles due to the very low number of isolates for each combination. All the isolates resistant to imipenem/relebactam (n=29), ceftazidime/avibactam (n=26) or meropenem/vaborbactam (n=8) were metallo-β-lactamase-(MBL)-producers, except for one OXA-48-producer that was resistant to meropenem/vaborbactam and imipenem/relebactam, one OXA-48-producer resistant to imipenem/relebactam and ceftazidime/avibactam, and two isolates resistant to ceftazidime/avibactam producers of KPC-3/OXA-48 and OXA-48 carbepenemases, respectively. Three of these isolates presented insertions in the ompC gene in comparison with ertapenem-resistant isolates from this same study but susceptible to the rest of the carbapenems and to 3rd generation cephalosporins and their combinations. Insertions detected were E87_N88insD in one OXA-48 producer, and T225_A226insGLNGYG in one OXA-48- and one KPC-3/OXA-48-producing isolates. However, these insertions do not generate a premature stop codon. ompF mutations were not detected.

Table 3

Table 3 Number and percentages of the CP-Eco isolates that are resistant to the tested antibiotics grouped by the carbapenemase gene that they harboured.

Six of the seven isolates resistant to cefiderocol were NDM-producers (4 bla_NDM-5 and 2 bla_NDM-1) and one had bla_VIM-1. Four NDM-producers belonged to ST167 and had an insertion in the pbp3 gene (Y334insYRIN), in addition to two other mutations in the same gene (E349K and I532L). Three of these isolates additionally had a frameshift mutation in the siderophore gene cirA (V89fs). Another NDM-producer (ST131) had a truncation in cirA (W243trunc), and the final NDM-producer belonged to ST405, had the Y334insYRIN insertion in the pbp3 gene, and a frameshift mutation in cirA (T508fs). The VIM-producer belonged to ST23 and had a frameshift mutation in cirA (T190fs). The median inhibition zone for NDM-producers was significantly lower compared to non-NDM-producers (20.5 versus 29.5; P = 0.00001).

There were six isolates that were resistant to all aminoglycosides, including plazomicin, and that carried rmt methyltransferases (3 rmtC, 1 rmtF1 and 1 rmtB1), with five of them additionally having bla_NDM-like. All 16 amikacin-resistant isolates harbored aac(6’)-Ib and/or methyltransferases genes, while the 46 tobramycin-resistant isolates had aac(3)-II, aac(6’)-Ib, ant(2’’)-Ia or rmtB1 genes (Supplementary Table S2). Two NDM-producers were resistant to aztreonam with no other β-lactamase gene detected, while all other MBL-producers resistant to this antibiotic harbored ESBL and/or AmpC genes (Supplementary Table S2).

Finally, 82 isolates (91.1%) had fosfomycin MICs of ≤ 8 mg/L and were considered susceptible (using EUCAST criteria for UTIs). Of those, 21 isolates (25.6%) had chromosomal mutations associated with resistance to fosfomycin. Of the 8 isolates resistant to this antibiotic, two harbored fosA genes, three had the chromosomal mutations ptsI_V25I/uhpT_E350Q, one had the chromosomal mutation uhpT_E350Q, and in the last two no known genetic mutation or gene was identified. There was a significant correlation between harboring mutations or genes associated with fosfomycin resistance and phenotypic resistance to this antibiotic (P = 0.003).

3.6 Characterization of the virulence-associated genes

Table 4 shows the detected virulence-associated genes (13 virulence genes considered), and the number of infection-causing isolates harboring them. The isolates had a median of two virulence-associated genes (range=0-9) where adhesin-coding genes were the most prevalent (68.9%). Thirty-eight isolates (42.2%) had four or more virulence-associated genes, 19 of which were causing infections. The most common infection in this group was UTI (9), followed by bacteremia (5), IAIs (2), wound infections (2), and abscess (1) (Supplementary Table S2). Isolates with 4 or more virulence-associated genes included three isolates belonging to ST131 and two belonging to ST88; and had similar averages of ARGs and susceptibility rates to the different antibiotic categories. Nineteen of these 34 isolates had bla_OXA-48, 8 had bla_VIM-1, 3 had bla_NDM-5, 2 had bla_NDM-1, 1 had bla_KPC-3, and 1 had both bla_OXA-48 and bla_KPC-3. Of the isolates pertaining to ST131, all had fimH and over half of them harbored pap and ompA. Over half of those pertaining to ST88 harbored fimH, iroN, and ompA.

Table 4

Table 4 Virulence-associated genes detected in the all 90 carbapenemase-producing E. coli isolates, in addition to the two most frequent sequence types and the number of isolates harboring these genes that were causing infections.

The adhesin-coding pap and the toxin-coding sat genes were found to be more frequent among isolates pertaining to ST131 (P = 0.02 and 0.0004, respectively). Particularly, five of the 9 isolates positive for pap and five out of the six positive for sat belonged to clade C. Additionally, the adhesion-coding fimH gene had almost double the frequency among isolates harboring bla_OXA-48 as compared to isolates harboring other carbapenemase genes (P = 0.04). There were no statistically significant differences in the virulence-associated genes between bloodstream isolates and isolates from other sources.

4 Discussion

In this study, 90 representative CP-Eco isolates obtained from 25 Spanish provinces from 2015 to 2020 were characterized. The inclusion criteria were designed to achieve the greatest possible representation of geographic areas, clones, and carbapenemases, as well as to avoid overrepresentation due to outbreaks and/or endemics. They do not reflect the real prevalence of the carbapenemase genes; an approach to this more general view is provided in section 3.1, which mirrored the epidemiological situation of CP-Eco in Spain (Cañada-García et al., 2022). Notably, CP-Eco isolates were found in 25 Spanish provinces, and major STs, such as ST131, showed inter-provincial spread. The studied CP-Eco isolates were highly polyclonal, unlike other Enterobacterales, such as Klebsiella pneumoniae that tend to show less genetic diversity (Oteo et al., 2015; Lázaro-Perona et al., 2022).

The high-risk ST131 was the most prevalent in this collection; this ST also represented 36% of CP-Eco isolates from 16 countries between 2008 and 2013 (Patiño-Navarrete et al., 2020). In Spain, ST131 represented 16.5% of CP-Eco isolates collected between 2012-2014 (Ortega et al., 2016). Furthermore, isolates belonging to this ST harbored all the carbapenemase genes detected in this study, on different plasmids and in different combinations with other acquired ARGs. This demonstrates its potential to act as a vector for disseminating resistance genes in different scenarios. Half of the ST131 isolates caused infections, and all of them were community-acquired. Moreover, a significant association was found between ST131 and the virulence genes sat (toxin-coding) and pap (adhesin-coding), highlighting its potential to cause life-threatening community-acquired infections.

The ST167/NDM-5 clone, observed in this study, had been previously described (Poirel et al., 2022b) and had shown to spread rapidly across European borders (Linkevicius et al., 2023). These isolates were also resistant to cefiderocol and detected in three different provinces. Of note, 6/7 of the cefiderocol-resistant isolates were also NDM-producers, an association that has been previously drawn, especially among ST167 isolates having the same mutations in pbp3 as shown in our study and a premature stop codon in the siderophore gene cirA (Wang et al., 2022). Additionally, NDM-producers were found to have tighter inhibition zones towards cefiderocol compared to non-NDM-producers, in line with recent findings where NDM-producers had elevated MICs towards this antibiotic (Poirel et al., 2022a). This hints at an interplay between NDM and cefiderocol resistance that warrants further investigation.

Another finding in our study is the detection of a relatively high rate of MBL-producers in isolates recovered from the community setting. The high average age of these patients could suggest an origin in nursing homes that are known reservoirs of CP-producing bacteria (Palacios-Baena et al., 2016) and whose flow of admissions-discharges from and to the community is difficult to trace. Moreover, bla_VIM-1 was detected in small IncL plasmids similar to the highly transmissible ones that carried bla_OXA-48 (Boyd et al., 2022). This combination was detected in isolates belonging to different STs, as previously described in K. pneumoniae (Cañada-García et al., 2023), and was predominantly found in healthcare-associated settings. This raises concern regarding the fast dissemination of bla_VIM-1, especially in Spain where the prevalence of VIM-producers is higher than other countries (Johnston et al., 2021; Cañada-García et al., 2022). Of note, bla_OXA-48 was in IncFII plasmids in eight isolates, an observation that was reported in another study though bla_OXA-48 was part of a different transposon (Moussa et al., 2020).

The CP-Eco isolates harboring bla_OXA-48 were significantly more susceptible to imipenem and meropenem compared to CP-Eco isolates harboring other carbapenemases, in accordance with the known weaker hydrolytic activity of this enzyme towards carbapenems (López-Camacho et al., 2019). Moreover, all OXA-48-like-producers and all but one KPC-producer did not harbor any other ARG on their plasmids, unlike other Enterobacterales such as K. pneumoniae where it is typical to find other ARGs with bla_KPC (Cañada-García et al., 2022; Cañada-García et al., 2023). In contrast, 30-40% of MBL-producers were resistant to more than six antibiotic categories as defined in this study. All the “new” antibiotics and antibiotic/inhibitor combinations, except imipenem/relebactam, showed over 75% in-vitro activity towards the CP-Eco isolates. Resistance to imipenem/relebactam, meropenem/vaborbactam or ceftzidime/avibactam was mainly due to MBL production (Galani et al., 2019). Three of the only four isolates resistant to at least one of these combinations had insertions in the OmpC porin gene, in comparison with sequences from OXA-48-producing isolates of this study only resistant to ertapenem but susceptible to the rest of the carbapenems and 3rd generation cephalosporins. The presence of mutations in the main porins of Enterobacterales has previously been related to resistance to the new beta-lactamase inhibitors combinations (Navarro-Carrera et al., 2022). VIM-producers appeared to be more resistant towards imipenem/relebactam than imipenem alone, which has been previously reported in another study (Hernández-García et al., 2022). However, this was mainly caused by the EUCAST cutoff values (European Society of Clinical Microbiology and Infectious Diseases (EUCAST), 2022) used to determine susceptibility since six of the eight VIM-1 producing isolates resistant to imipenem/relebactam but susceptible to imipenem had MICs of 4 mg/L for both antibiotics. There were also seven NDM-producers that were susceptible to meropenem/vaborbactam according to the EUCAST guidelines. However, these isolates had elevated MICs, and most would have been considered resistant in case the CLSI guidelines were followed (Jean et al., 2019).

Six isolates were resistant to plazomicin, five of them harbored bla_NDM-like and were the only isolates in which rmt methyltransferases were detected. Furthermore, unlike other Enterobacterales in which rmt methyltransferase gene seems to disseminate via clonal expansion (Arca-Suárez et al., 2022), in our collection this gene was detected in isolates from different STs and were present in different plasmids.

Fosfomycin is one of the main antibiotics used to treat community-acquired uncomplicated cystitis caused by E. coli in Spain (Sojo-Dorado et al., 2022). In this study, most of the isolates that had mutations associated with resistance to this antibiotic harbored chromosomal mutations, and only two isolates had fos genes, similarly to what was previously described (Falagas et al., 2019). Moreover, the fosA7 genes detected in the chromosome of one of the isolates has also been described to be able to be present there (Pan et al., 2021). The correlation between resistance to fosfomycin and the presence of resistance mechanisms against this antibiotic must be taken with caution since cut-off points for UTI have been used (the only ones established by EUCAST), of which there were only 26 isolates in this collection. Moreover, although a good agreement was reported between using E-test strips and disk diffusion assays (Hirsch et al., 2015; Goer et al., 2022), the EUCAST does not explicitly recommend the use of E-tests for the determination of cutoff values for fosfomycin. Nevertheless, these isolates could be displaying a heterorresistant phenotype that appeared to be susceptible in our testing methods and criteria (Campos et al., 2020).

In conclusion, our data shows a wide range of genetically diverse CP-E. coli isolates present across several Spanish geographical areas, with a heterogeneous set of ARGs, plasmids, and virulence-associated genes. These findings highlight the potential of this bacteria to acquire and disseminate ARGs, and to cause life-threatening, difficult-to-treat infections in different settings, including the community. The spread of CP-Eco in the community, like what already happened with CTX-M-15-producing E. coli, would pose a threat of great impact to the general health of the population. The integration of WGS in the surveillance and diagnosis of CP-Eco, in line with ECDC recommendations, can help control and minimize this impact.

Data availability statement

All the data generated for this study is available in the manuscript and its related Supplementary Material. The sequenced raw reads have been deposited in the open-access European Nucleotide Archive database under the accession number PRJEB70795.

Author contributions

ED: Writing – review & editing, Writing – original draft, Validation, Methodology, Investigation, Data curation. LG-M: Writing – review & editing, Writing – original draft, Methodology, Investigation, Data curation. JC-G: Writing – review & editing, Writing – original draft, Visualization, Software, Investigation, Data curation. ER: Writing – review & editing, Writing – original draft, Investigation, Formal Analysis, Data curation. AS-G: Validation, Formal Analysis, Writing – review & editing, Writing – original draft, Resources. IS-R: Formal Analysis, Writing – review & editing, Writing – original draft, Validation, Resources. LL-U: Validation, Formal Analysis, Writing – review & editing, Writing – original draft, Resources. PI: Validation, Formal Analysis, Writing – review & editing, Writing – original draft, Resources. AG-P: Validation, Resources, Formal Analysis, Writing – review & editing, Writing – original draft. JS: Writing – review & editing, Writing – original draft, Methodology, Investigation. DV-M: Writing – review & editing, Writing – original draft, Software, Investigation, Data curation. IA-G: Writing – review & editing, Writing – original draft, Investigation, Data curation. VB: Writing – review & editing, Writing – original draft, Methodology, Investigation. NL: Writing – review & editing, Writing – original draft, Methodology, Investigation. SG-C: Writing – review & editing, Writing – original draft, Investigation, Formal Analysis, Data curation. EC: Validation, Writing – review & editing, Writing – original draft, Resources, Formal Analysis. BA: Writing – review & editing, Writing – original draft, Validation, Supervision, Resources, Project administration, Formal Analysis. JO-I: Writing – review & editing, Writing – original draft, Validation, Supervision, Resources, Project administration, Funding acquisition, Formal Analysis. MP-V: Writing – review & editing, Writing – original draft, Visualization, Supervision, Software, Resources, Formal Analysis.

Other members of the Spanish Eco-Carba study group are:

Verónica Casquero and Olga Valiente (Centro Nacional de Microbiología, Instituto de Salud Carlos III, Madrid), Almudena Alhambra Mosquera (Hospital Universitario HM Montepríncipe, Madrid), Alia Eworo Ndongo (Hospital de El Escorial, Madrid), Susana Hernando Real (Hospital General de Segovia), Luis Moisés Ruiz-Velasco (Hospital Central De La Cruz Roja San José y Santa Adela, Madrid), José Leiva (Clínica Universitaria de Navarra), Nieves Balado (Hospital Xeral de Vigo, Pontevedra), Adriana Ortega (Fundación Hospital de Alcorcón, Madrid. CIBERINFEC), Mar Olga Pérez Moreno (Hospital de Tortosa Verge de La Cinta, Tortosa, Tarragona), Ana Bordes (Hospital De Gran Canaria Dr. Negrín, Las Palmas de Gran Canaria), Cristobal del Rosario Quintana (Hospital Universitario Insular de Gran Canaria), María Eugenia Portillo (Hospital Universitario de Navarra), Caridad Sainz de Baranda (Complejo Hospitalario Universitario de Albacete), Gloria Trujillo (Hospital San Joan de Deu (Althaia), Barcelona), Begoña Palop (Hospital Regional de Málaga), Carmen Aldea-Mansilla (Complejo Asistencial Universitario de Soria), Juan Cuadros (Hospital Universitario Príncipe de Asturias, Madrid), Yolanda Gil (Hospital de Móstoles, Madrid), Soledad Illescas Fernández-Bermejo (Hospital General Universitario de Ciudad Real), Ana Ramos (Hospital Universitario de Torrejón, Madrid), Salvador Giner (Hospital Universitario La Fe, Valencia), Antonio Casabella Pernas (Corporación Sanitaria Parc Taulí, Barcelona), M. Pilar Ortega Lafont (Complejo Asistencial Universitario de Burgos), María Huertas Vaquero (Hospital General La Mancha Centro, Ciudad Real), Isabel Antolín (Hospital de Medina del Campo, Valladolid), M^a de los Ángeles Pallarés (Hospital Montecelo, Pontevedra), Beatriz Iglesias (Hospital San Agustín, Asturias), Frederic Gómez-Bertomeu (Hospital Universitario Joan XXIII, Tarragona; CIBERINFEC), Ana Isabel López-Calleja (Hospital Universitario Miguel Servet, Zaragoza), Pilar Zamarrón (Hospital Virgen de La Salud, Toledo)

Funding

The author(s) declare financial support was received for the research, authorship, and/or publication of this article. This research was supported by CIBER—Consorcio Centro de Investigación Biomédica en Red (CIBERINFEC and CIBERES), Instituto de Salud Carlos III, Ministerio de Ciencia e Innovación and Unión Europea-NextGenerationEU; and by the Antibiotic Resistance Surveillance Programs of the National Center for Microbiology, Instituto de Salud Carlos III. This research was also supported by Personalized and precision medicine grant from the Instituto de Salud Carlos III (MePRAM Project, PMP22/00092), Instituto de Salud Carlos III, Ministerio de Ciencia e Innovación. Funded by NextGenerationEU funds from the European Union that finance the actions of the The Resilience and Recovery Facility.

Acknowledgments

We thank the Genomics Unit of the Centro Nacional de Microbiología for performing the DNA sequencing.

Conflict of interest

The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.

Publisher’s note

All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.

Supplementary material

The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fcimb.2024.1390966/full#supplementary-material

References

Arca-Suárez, J., Rodiño-Janeiro, B. K., Pérez, A., Guijarro-Sánchez, P., Vázquez-Ucha, J. C., Cruz, F., et al. (2022). Emergence of 16S rRNA methyltransferases among carbapenemase-producing Enterobacterales in Spain studied by whole-genome sequencing. Int. J. Antimicrob. Agents 59, 106456. doi: 10.1016/j.ijantimicag.2021.106456

PubMed Abstract | CrossRef Full Text | Google Scholar

Balasubramanian, R., Van Boeckel, T. P., Carmeli, Y., Cosgrove, S., Laxminarayan, R. (2023). Global incidence in hospital-associated infections resistant to antibiotics: An analysis of point prevalence surveys from 99 countries. PloS Med. 20, e1004178. doi: 10.1371/journal.pmed.1004178

PubMed Abstract | CrossRef Full Text | Google Scholar

Boutzoukas, A. E., Komarow, L., Chen, L., Hanson, B., Kanj, S. S., Liu, Z., et al. (2023). International epidemiology of carbapenemase-producing escherichia coli. Clin. Infect. Dis. 77, 499–509. doi: 10.1093/cid/ciad288

PubMed Abstract | CrossRef Full Text | Google Scholar

Boyd, S. E., Holmes, A., Peck, R., Livermore, D. M., Hope, W. (2022). OXA-48-like β-lactamases: global epidemiology, treatment options, and development pipeline. Antimicrob. Agents Chemother. 66, e0021622. doi: 10.1128/aac.00216-22

PubMed Abstract | CrossRef Full Text | Google Scholar

Campos, A. C. C., Andrade, N. L., Couto, N., Mutters, N., de Vos, M., de, P., et al. (2020). Characterization of fosfomycin heteroresistance among multidrug-resistant Escherichia coli isolates from hospitalized patients in Rio de Janeiro, Brazil. J. Global Antimicrobial Resistance 22, 584–593. doi: 10.1016/j.jgar.2020.04.026

CrossRef Full Text | Google Scholar

Cañada-García, J. E., Moure, Z., Sola-Campoy, P. J., Delgado-Valverde, M., Cano, M. E., Gijón, D., et al. (2022). CARB-ES-19 Multicenter Study of Carbapenemase-Producing Klebsiella pneumoniae and Escherichia coli From All Spanish Provinces Reveals Interregional Spread of High-Risk Clones Such as ST307/OXA-48 and ST512/KPC-3. Front. Microbiol. 13. doi: 10.3389/fmicb.2022.918362

CrossRef Full Text | Google Scholar

Cañada-García, J. E., Ramírez de Arellano, E., Jiménez-Orellana, M., Viedma, E., Sánchez, A., Alhambra, A., et al. (2023). Carbapenemase-producing klebsiella pneumoniae in COVID-19 intensive care patients: identification of IncL-VIM-1 plasmid in previously non-predominant sequence types. Antibiotics 12, 107. doi: 10.3390/antibiotics12010107

PubMed Abstract | CrossRef Full Text | Google Scholar

Cassini, A., Högberg, L. D., Plachouras, D., Quattrocchi, A., Hoxha, A., Simonsen, G. S., et al. (2019). Attributable deaths and disability-adjusted life-years caused by infections with antibiotic-resistant bacteria in the EU and the European Economic Area in 2015: a population-level modelling analysis. Lancet Infect. Dis. 19, 56–66. doi: 10.1016/S1473-3099(18)30605-4

PubMed Abstract | CrossRef Full Text | Google Scholar

Čurová, K., Slebodníková, R., Kmeťová, M., Hrabovský, V., Maruniak, M., Liptáková, E., et al. (2020). Virulence, phylogenetic background and antimicrobial resistance in Escherichia coli associated with extraintestinal infections. J. Infect. Public Health 13, 1537–1543. doi: 10.1016/j.jiph.2020.06.032

PubMed Abstract | CrossRef Full Text | Google Scholar

European Society of Clinical Microbiology and Infectious Diseases (EUCAST) (2022) European Committee on Antimicrobial Susceptibility Testing. Breakpoint tables for interpretation of MICs and zone diameters. Version 12.0, 2022. Available online at: https://www.eucast.org/ast_of_bacteria/previous_versions_of_documents (Accessed 23/02/2024).

Google Scholar

Falagas, M. E., Athanasaki, F., Voulgaris, G. L., Triarides, N. A., Vardakas, K. Z. (2019). Resistance to fosfomycin: mechanisms, frequency and clinical consequences. Int. J. Antimicrob. Agents 53, 22–28. doi: 10.1016/j.ijantimicag.2018.09.013

PubMed Abstract | CrossRef Full Text | Google Scholar

Galani, I., Souli, M., Nafplioti, K., Adamou, P., Karaiskos, I., Giamarellou, H., et al. (2019). In vitro activity of imipenem-relebactam against non-MBL carbapenemase-producing Klebsiella pneumoniae isolated in Greek hospitals in 2015–2016. Eur. J. Clin. Microbiol. Infect. Dis. 38, 1143–1150. doi: 10.1007/s10096-019-03517-y

PubMed Abstract | CrossRef Full Text | Google Scholar

Gastmeier, P., Schwab, F., Bärwolff, S., Rüden, H., Grundmann, H. (2006). Correlation between the genetic diversity of nosocomial pathogens and their survival time in intensive care units. J. Hosp. Infection 62, 181–186. doi: 10.1016/j.jhin.2005.08.010

CrossRef Full Text | Google Scholar

Giske, C. G., Martinez-Martinez, L., Cantón Spain, R., Stefani, S., Skov, R., Glupczynski, Y., et al. (2017) EUCAST guidelines for detection of resistance mechanisms and specific resistances of clinical and/or epidemiological importance 2017. Available online at: https://www.eucast.org/fileadmin/src/media/PDFs/EUCAST_files/Resistance_mechanisms/EUCAST_detection_of_resistance_mechanisms_170711.pdf (Accessed 23/02/2024).

Google Scholar

Goer, A., Blanchard, L. S., Van Belkum, A., Loftus, K. J., Armstrong, T. P., Gatermann, S. G., et al. (2022). Multicenter evaluation of the novel ETEST fosfomycin for antimicrobial susceptibility testing of enterobacterales, enterococcus faecalis, and staphylococcus species. J. Clin. Microbiol. 60, e0002122. doi: 10.1128/jcm.00021-22

PubMed Abstract | CrossRef Full Text | Google Scholar

González, D., Gallagher, E., Zúñiga, T., Leiva, J., Vitas, A. I. (2020). Prevalence and characterization of β-lactamase-producing Enterobacteriaceae in healthy human carriers. Int. Microbiol. 23, 171–177. doi: 10.1007/s10123-019-00087-z

PubMed Abstract | CrossRef Full Text | Google Scholar

Grundmann, H., Glasner, C., Albiger, B., Aanensen, D. M., Tomlinson, C. T., Andrasević, A. T., et al. (2017). Occurrence of carbapenemase-producing Klebsiella pneumoniae and Escherichia coli in the European survey of carbapenemase-producing Enterobacteriaceae (EuSCAPE): a prospective, multinational study. Lancet Infect. Dis. 17, 153–163. doi: 10.1016/S1473-3099(16)30257-2

PubMed Abstract | CrossRef Full Text | Google Scholar

Hernández-García, M., García-Castillo, M., Bou, G., Cercenado, E., Delgado-Valverde, M., Oliver, A., et al. (2022). Imipenem-relebactam susceptibility in enterobacterales isolates recovered from ICU patients from Spain and Portugal (SUPERIOR and STEP studies). Microbiol. Spectr. 10, e0292722. doi: 10.1128/spectrum.02927-22

PubMed Abstract | CrossRef Full Text | Google Scholar

Hirsch, E. B., Raux, B. R., Zucchi, P. C., Kim, Y., McCoy, C., Kirby, J. E., et al. (2015). Activity of fosfomycin and comparison of several susceptibility testing methods against contemporary urine isolates. Int. J. Antimicrob. Agents 46, 642–647. doi: 10.1016/j.ijantimicag.2015.08.012

PubMed Abstract | CrossRef Full Text | Google Scholar

Hu, F., Yuan, L., Yang, Y., Xu, Y., Huang, Y., Hu, Y., et al. (2022). A multicenter investigation of 2,773 cases of bloodstream infections based on China antimicrobial surveillance network (CHINET). Front. Cell Infect. Microbiol. 12. doi: 10.3389/fcimb.2022.1075185

CrossRef Full Text | Google Scholar

Hunt, M., Mather, A. E., Sánchez-Busó, L., Page, A. J., Parkhill, J., Keane, J. A., et al. (2017). ARIBA: rapid antimicrobial resistance genotyping directly from sequencing reads. Microb. Genom 3, e000131. doi: 10.1099/mgen.0.000131

PubMed Abstract | CrossRef Full Text | Google Scholar

ISO Central Secretariat. (2006). ch. de la Voie Creuse, Case postale 56, CH-1211 Genève 20, Switzerland. Available at: https://www.iso.org/files/live/sites/isoorg/files/archive/pdf/en/annual_report_2006.pdf.

Google Scholar

Jean, S., Garrett, S., Anglade, C., Bridon, L., Davies, L., Garner, O. B., et al. (2019). Multicenter clinical evaluation of etest meropenem-vaborbactam (bioMérieux) for susceptibility testing of enterobacterales (Enterobacteriaceae) and pseudomonas aeruginosa. J. Clin. Microbiol. 58, e01205–e01219. doi: 10.1128/JCM.01205-19

PubMed Abstract | CrossRef Full Text | Google Scholar

Johnston, B. D., Thuras, P., Porter, S. B., Anacker, M., VonBank, B., Vagnone, P. S., et al. (2021). Global molecular epidemiology of carbapenem-resistant Escherichia coli (2002–2017). Eur. J. Clin. Microbiol. Infect. Dis. doi: 10.1007/s10096-021-04310-6

CrossRef Full Text | Google Scholar

Lázaro-Perona, F., Dahdouh, E., Sotillo, A., Pérez-Blanco, V., Villa, J., Viedma, E., et al. (2022). Dissemination of a single ST11 clone of OXA-48-producing Klebsiella pneumoniae within a large polyclonal hospital outbreak determined by genomic sequencing. Microb. Genom 8, 000808. doi: 10.1099/mgen.0.000808

PubMed Abstract | CrossRef Full Text | Google Scholar

Linkevicius, M., Bonnin, R. A., Alm, E., Svartström, O., Apfalter, P., Hartl, R., et al. (2023). Rapid cross-border emergence of NDM-5-producing Escherichia coli in the European Union/European Economic Area, 2012 to June 2022. Eurosurveillance 28, 2300209. doi: 10.2807/1560-7917.ES.2023.28.19.2300209

PubMed Abstract | CrossRef Full Text | Google Scholar

López-Camacho, E., Paño-Pardo, J. R., Sotillo, A., Elías-López, C., Martínez-Martínez, L., Gómez-Gil, R., et al. (2019). Meropenem heteroresistance in clinical isolates of OXA-48–producing Klebsiella pneumoniae. Diagn. Microbiol. Infect. Dis. 93, 162–166. doi: 10.1016/j.diagmicrobio.2018.09.008

PubMed Abstract | CrossRef Full Text | Google Scholar

Magiorakos, A. P., Srinivasan, A., Carey, R. B., Carmeli, Y., Falagas, M. E., Giske, C. G., et al. (2012). Multidrug-resistant, extensively drug-resistant and pandrug-resistant bacteria: an international expert proposal for interim standard definitions for acquired resistance. Clin. Microbiol. Infection 18, 268–281. doi: 10.1111/j.1469-0691.2011.03570.x

CrossRef Full Text | Google Scholar

Moussa, J., Panossian, B., Nassour, E., Salloum, T., Abboud, E., Tokajian, S. (2020). Detailed characterization of an IncFII plasmid carrying bla_OXA-48 from Lebanon. J. Antimicrobial Chemotherapy 75, 2462–2465. doi: 10.1093/jac/dkaa181

CrossRef Full Text | Google Scholar

Navarro-Carrera, P., García-Clemente, P., Gómez-Gil, R., Cáceres-Sánchez, R., Bloise, I., Dahdouh, E., et al. (2022). Resistance to meropenem-vaborbactam in clinical isolates of Enterobacter roggenkampii without previous exposure to vaborbactam. Clin. Microbiol. Infection 28, 1662–1664. doi: 10.1016/j.cmi.2022.09.005

CrossRef Full Text | Google Scholar

Ortega, A., Sáez, D., Bautista, V., Fernández-Romero, S., Lara, N., Aracil, B., et al. (2016). Carbapenemase-producing Escherichia coli is becoming more prevalent in Spain mainly because of the polyclonal dissemination of OXA-48. J. Antimicrobial Chemotherapy 71, 2131–2138. doi: 10.1093/jac/dkw148

CrossRef Full Text | Google Scholar

Oteo, J., Ortega, A., Bartolomé, R., Bou, G., Conejo, C., Fernández-Martínez, M., et al. (2015). Prospective multicenter study of carbapenemase-producing enterobacteriaceae from 83 hospitals in Spain reveals high in vitro susceptibility to colistin and meropenem. Antimicrob. Agents Chemother. 59, 3406–3412. doi: 10.1128/AAC.00086-15

PubMed Abstract | CrossRef Full Text | Google Scholar

Palacios-Baena, Z. R., Oteo, J., Conejo, C., Larrosa, M. N., Bou, G., Fernández-Martínez, M., et al. (2016). Comprehensive clinical and epidemiological assessment of colonisation and infection due to carbapenemase-producing Enterobacteriaceae in Spain. J. Infection 72, 152–160. doi: 10.1016/j.jinf.2015.10.008

CrossRef Full Text | Google Scholar

Pan, Y., Hu, B., Bai, X., Yang, X., Cao, L., Liu, Q., et al. (2021). Antimicrobial resistance of Non-O157 shiga toxin-producing escherichia coli isolated from humans and domestic animals. Antibiotics (Basel) 10, 74. doi: 10.3390/antibiotics10010074

PubMed Abstract | CrossRef Full Text | Google Scholar

Patiño-Navarrete, R., Rosinski-Chupin, I., Cabanel, N., Gauthier, L., Takissian, J., Madec, J. Y., et al. (2020). Stepwise evolution and convergent recombination underlie the global dissemination of carbapenemase-producing Escherichia coli. Genome Med. 12, 10. doi: 10.1186/s13073-019-0699-6

PubMed Abstract | CrossRef Full Text | Google Scholar

Pérez-Vazquez, M., Oteo-Iglesias, J., Sola-Campoy, P. J., Carrizo-Manzoni, H., Bautista, V., Lara, N., et al. (2019). Spanish antibiotic resistance surveillance program collaborating group. Characterization of carbapenemase-producing klebsiella oxytoca in Spain, 2016–2017. Antimicrob. Agents Chemother. 63, e02529–e02518. doi: 10.1128/AAC.02529-18

PubMed Abstract | CrossRef Full Text | Google Scholar

Plachouras, D., Lepape, A., Suetens, C. (2018). ECDC definitions and methods for the surveillance of healthcare-associated infections in intensive care units. Intensive Care Med. 44, 2216–2218. doi: 10.1007/s00134-018-5113-0

PubMed Abstract | CrossRef Full Text | Google Scholar

Poirel, L., Ortiz de la Rosa, J.-M., Sadek, M., Nordmann, P. (2022a). Impact of Acquired Broad-Spectrum β-Lactamases on Susceptibility to Cefiderocol and Newly Developed β-Lactam/β-Lactamase Inhibitor Combinations in Escherichia coli and Pseudomonas aeruginosa. Antimicrob. Agents Chemother. 66, e0003922. doi: 10.1128/aac.00039-22

PubMed Abstract | CrossRef Full Text | Google Scholar

Poirel, L., Ortiz de la Rosa, J. M., Sakaoglu, Z., Kusaksizoglu, A., Sadek, M., Nordmann, P. (2022b). NDM-35-producing ST167 escherichia coli highly resistant to β-lactams including cefiderocol. Antimicrob. Agents Chemother. 66, e0031122. doi: 10.1128/aac.00311-22

PubMed Abstract | CrossRef Full Text | Google Scholar

Seemann, T. (2014). Prokka: rapid prokaryotic genome annotation. Bioinformatics 30, 2068–2069. doi: 10.1093/bioinformatics/btu153

PubMed Abstract | CrossRef Full Text | Google Scholar

Sojo-Dorado, J., López-Hernández, I., Rosso-Fernandez, C., Morales, I. M., Palacios-Baena, Z. R., Hernández-Torres, A., et al. (2022). Effectiveness of fosfomycin for the treatment of multidrug-resistant escherichia coli bacteremic urinary tract infections. JAMA Netw. Open 5, e2137277. doi: 10.1001/jamanetworkopen.2021.37277

PubMed Abstract | CrossRef Full Text | Google Scholar

Tacconelli, E., Carrara, E., Savoldi, A., Harbarth, S., Mendelson, M., Monnet, D. L., et al. (2018). Discovery, research, and development of new antibiotics: the WHO priority list of antibiotic-resistant bacteria and tuberculosis. Lancet Infect. Dis. 18, 318–327. doi: 10.1016/S1473-3099(17)30753-3

PubMed Abstract | CrossRef Full Text | Google Scholar

Wang, Q., Jin, L., Sun, S., Yin, Y., Wang, R., Chen, F., et al. (2022). Occurrence of High Levels of Cefiderocol Resistance in Carbapenem-Resistant Escherichia coli before Its Approval in China: a Report from China CRE-Network. Microbiol. Spectr. 10, e0267021. doi: 10.1128/spectrum.02670-21

PubMed Abstract | CrossRef Full Text | Google Scholar

Wick, R. R., Judd, L. M., Gorrie, C. L., Holt, K. E. (2017). Unicycler: Resolving bacterial genome assemblies from short and long sequencing reads. PloS Comput. Biol. 13, e1005595. doi: 10.1371/journal.pcbi.1005595

PubMed Abstract | CrossRef Full Text | Google Scholar