
==== Front
JAC Antimicrob Resist
JAC Antimicrob Resist
jacamr
JAC-Antimicrobial Resistance
2632-1823
Oxford University Press UK

10.1093/jacamr/dlae073
dlae073
Original Article
AcademicSubjects/MED00740
AcademicSubjects/SCI01150
Emergence of an IncX3 plasmid co-harbouring the carbapenemase genes blaNDM-5 and blaOXA-181
Zuo Hui Antimicrobial Resistance Research Center, National Institute of Infectious Diseases, Tokyo, Japan

https://orcid.org/0000-0002-9897-3905
Sugawara Yo Antimicrobial Resistance Research Center, National Institute of Infectious Diseases, Tokyo, Japan

Kondo Kohei Antimicrobial Resistance Research Center, National Institute of Infectious Diseases, Tokyo, Japan

Kayama Shizuo Antimicrobial Resistance Research Center, National Institute of Infectious Diseases, Tokyo, Japan

Kawakami Sayoko Antimicrobial Resistance Research Center, National Institute of Infectious Diseases, Tokyo, Japan

Uechi Kohei Division of Clinical Laboratory and Blood Transfusion, University of the Ryukyus Hospital, Okinawa, Japan

Nakano Ami Division of Clinical Laboratory and Blood Transfusion, University of the Ryukyus Hospital, Okinawa, Japan

Yahara Koji Antimicrobial Resistance Research Center, National Institute of Infectious Diseases, Tokyo, Japan

Sugai Motoyuki Antimicrobial Resistance Research Center, National Institute of Infectious Diseases, Tokyo, Japan

Corresponding author. E-mail: suga-yo@niid.go.jp
6 2024
13 5 2024
13 5 2024
6 3 dlae07325 9 2023
23 4 2024
© The Author(s) 2024. Published by Oxford University Press on behalf of British Society for Antimicrobial Chemotherapy.
2024
https://creativecommons.org/licenses/by/4.0/ This is an Open Access article distributed under the terms of the Creative Commons Attribution License (https://creativecommons.org/licenses/by/4.0/), which permits unrestricted reuse, distribution, and reproduction in any medium, provided the original work is properly cited.

Abstract

Background

The spread of transmissible plasmids with carbapenemase genes has contributed to a global increase in carbapenemase-producing Enterobacterales over the past two decades, with blaNDM and blaOXA among the most prevalent carbapenemase genes.

Objectives

To characterize an Escherichia coli isolate co-carrying blaNDM-5 and blaOXA-181 (JBEHAAB-19-0176) that was isolated in the Japan Antimicrobial Resistant Bacterial Surveillance in 2019–20, and to evaluate the functional advantage of carrying both genes as opposed to only one.

Methods

The whole-genome sequence of the isolate was determined using long- and short-read sequencing. Growth assay and co-culture experiments were performed for phenotypic characterization in the presence of different β-lactam antibiotics.

Results

WGS analysis showed that blaNDM-5 and blaOXA-181 were carried by the same IncX3 plasmid, pJBEHAAB-19-0176_NDM-OXA. Genetic characterization of the plasmid suggested that the plasmid emerged through the formation of a co-integrate and resolution of two typical IncX3 plasmids harbouring blaNDM-5 and blaOXA-181, which involved two recombination events at the IS3000 and IS26 sequences. When cultured in the presence of piperacillin or cefpodoxime, the growth rate of the transformant co-harbouring blaNDM-5 and blaOXA-181 was significantly higher than the transformant with only blaNDM-5. Furthermore, in co-culture where the two blaNDM-5-harbouring transformants were allowed to compete directly, the strain additionally harbouring blaOXA-181 showed a marked growth advantage.

Conclusions

The additional carriage of blaOXA-181 confers a selective advantage to bacteria in the presence of piperacillin and cefpodoxime. These findings may explain the current epidemiology of carbapenemase-producing Enterobacterales, in which bacteria carrying both blaNDM-5 and blaOXA-48-like genes have emerged independently worldwide.

Japan Agency for Medical Research and Development 10.13039/100009619 JP19fk0108061 23fk0108604 JSPS 10.13039/501100001691 JP20K07496
==== Body
pmcIntroduction

There has been a global emergence of carbapenem-resistant Enterobacterales with escalating resistance to all available β-lactam antibiotics, which WHO has highlighted as critical priority pathogens.1 The major mechanism underlying carbapenem resistance is the degradation of carbapenems via bacterial enzymes called carbapenemases. One such enzyme is NDM, which belongs to the Ambler class B β-lactamases and is encoded by blaNDM, which confers resistance not only to carbapenems but also to almost all β-lactams.2 The blaNDM gene is disseminated worldwide, predominantly in Asia,2 and is usually carried by transferable plasmids, such as IncX3-type plasmids, enabling it to spread broadly among species.3–5 The OXA-48 enzyme is also a carbapenemase of concern. It is difficult to detect because of the low-level in vitro resistance to carbapenem antibiotics it confers, posing a challenge to clinical laboratories.6,7 Genes encoding OXA-48-like enzymes, such as blaOXA-48, blaOXA-181 and blaOXA-232, are globally widespread.6,7 The blaNDM and blaOXA-48-like genes were identified as the second and third most prevalent carbapenemase genes among carbapenemase-producing Enterobacterales (CPE), respectively, in a global surveillance programme conducted between 2008 and 2014.8 Specifically, blaOXA-181 and blaNDM-5 were the first and second most prevalent carbapenemase genes, respectively, among global carbapenemase-producing Escherichia coli strains collected in 2015–17.9 Although the sole presence of blaNDM confers sufficiently high-level carbapenem resistance to bacteria, there have been many previous reports on isolates co-carrying blaNDM and blaOXA-48-like genes.10 Most harbour these two genes on two separate plasmids, whereas only a few co-harbour them on the same plasmid, such as the Acinetobacter plasmid GR59.11 Here, we report an IncX3 plasmid co-harbouring blaNDM-5 and blaOXA-181 and explore the potential advantage of carrying both of these genes through comparative phenotypic analysis with typical endemic IncX3 plasmids harbouring only one of the two genes.

Materials and methods

Bacterial isolates and antimicrobial susceptibility testing

An E. coli isolate, JBEHAAB-19-0176, was isolated during the Japan Antimicrobial Resistant Bacterial Surveillance, focusing on Gram-negative bacteria (JARBS-GNR) that had been conducted during 2019–20, which also focused on nosocomial Enterobacterales isolates with reduced carbapenem susceptibility and/or resistance to third-generation cephalosporins.12 The current study was approved by the International Review Board of the National Institute of Infectious Diseases (approval number: 1553). The isolate was recovered from the faeces of an outpatient at the University of Ryukyus Hospital (Okinawa, Japan). Antimicrobial susceptibility testing was performed using the MicroScan WalkAway Plus system with Neg MIC EN 2J and 3.31E panels (Beckman Coulter, Brea, CA, USA). The MICs of meropenem, imipenem, ampicillin, piperacillin and cefpodoxime were further determined using a broth microdilution method, according to the CLSI guidelines.

Bioinformatics and plasmid replicon analysis

Genomic and plasmid DNA were extracted using a Monarch HMW DNA Extraction Kit for Tissues (New England Biolabs, Ipswich, MA, USA) and a QIAGEN Plasmid Mini kit (QIAGEN, Hilden, Germany), respectively, and WGS data obtained using HiSeq X (Illumina, San Diego, CA, USA) and GridION systems (Oxford Nanopore Technologies, Oxford, UK), as previously described.12 Obtained raw reads were assembled, as previously described,12 and antimicrobial resistance gene detection and plasmid Inc typing performed using ResFinder v2.113 and PlasmidFinder v1.3,14 respectively. Plasmid sequences were annotated using DFAST v1.2.015 and compared using BLAST (http://blast.ncbi.nlm.nih.gov) and Easyfig v2.2.3.16 Transposon and insertion sequences were identified using ISfinder.17

Transformation and conjugation

Plasmid DNA was introduced into E. coli HST08 cells (Takara Bio, Shiga, Japan) via electroporation using a Gene Pulser Xcell (Bio-Rad, Hercules, CA, USA). The size of the plasmids in the transformants and the presence of blaNDM-5 and blaOXA-181 on them were confirmed through S1 nuclease digestion of whole genomic DNA, followed by PFGE and Southern hybridization, as described previously.5 The conjugation assay was conducted as previously described, with slight modifications,5 using rifampicin-resistant E. coli ML490918 as the recipient. Transformants carrying IncX3 plasmids were mixed with recipient cells in a 1:10 ratio and incubated on a nitrocellulose membrane filter at 37°C for 2 h. The bacterial mixture was then suspended in brain heart infusion (BHI) broth and plated onto BHI agar containing meropenem (0.125 mg/L) and rifampicin (100 mg/L). The presence of blaNDM-5 and/or blaOXA-181 in the colonies was confirmed using colony-direct PCR. The conjugation frequency was calculated by dividing the number of cfu of the transconjugants by the number of cfu of the donor and transconjugants.

Growth curve and pairwise competition assay

Growth rates were determined according to a previously published study,19 with a few modifications. Briefly, bacterial suspensions of pJBEHAAB-19-0176_NDM-OXA and pJBBDAGF-19-0019_NDM-5 transformants were adjusted to OD600 = 0.5 and then diluted 100-fold in LB broth containing piperacillin (64 mg/L), cefpodoxime (64 mg/L), meropenem (128 mg/L), ampicillin (2048 mg/L), cefazolin (512 mg/L), cefoxitin (512 mg/L), flomoxef (256 mg/L), ceftazidime (512 mg/L), cefotaxime (512 mg/L), cefepime (128 mg/L) or imipenem (128 mg/L). Drug concentrations were set to maximize any differences in growth rates. The OD600 was measured for 36 h at 37°C with vigorous shaking (800 rpm) using a LogPhase 600 Microbiology Reader (Agilent Technologies, Santa Clara, CA, USA).

Competition assays were conducted as previously described,19 with slight modifications. Bacterial suspensions of the above two transformants were adjusted to OD600 = 0.5 and mixed at a 1:1 ratio. The mixture was then diluted 100-fold in LB broth with or without piperacillin (128 mg/L), cefpodoxime (64 mg/L) or meropenem (16 mg/L) and incubated at 37°C with shaking. The number of cells in each culture was evaluated at 24 h timepoints by spreading serially diluted cells and then culturing them on two types of BHI agar plates. One plate contained 0.25 mg/L meropenem to count the number of the two transformants, and the other plate contained 0.25 mg/L meropenem and 0.125 mg/L levofloxacin to count the number of pJBEHAAB-19-0176_NDM-OXA transformants carrying the quinolone resistance gene qnrS1. The competition indices were calculated by dividing the number of pJBEHAAB-19-0176_NDM-OXA transformants by the total number of transformants.

Plasmid construction and induction of blaOXA-181 expression

The blaOXA-181 gene was cloned into the arabinose-inducible expression vector, pBAD18-cm,20 as described previously.21 The pBAD18-cm and the one harbouring blaOXA-181 were introduced into E. coli HST08 carrying the pJBBDAGF-19-0019_NDM-5 plasmid via electroporation. Growth and competition assays were performed as described above in the presence of 0.02% arabinose and 30 mg/L chloramphenicol to induce blaOXA-181 expression and to maintain the introduced plasmids, respectively. For the competition assay, the ratio of E. coli cells with and without blaOXA-181 was determined using PCR. Each culture was spread onto BHI agar plates containing 0.25 mg/L meropenem and 30 mg/L chloramphenicol. At least 20 colonies per condition were subjected to colony-direct PCR using a primer pair targeting the upstream and downstream regions of the pBAD18-cm cloning site (5′-GATTAGCGGATCCTACCTGAC-3′ and 5′-CTTCTCTCATCCGCCAAAAC-3′). Competition indices were calculated by dividing the number of pJBBDAGF-19-0019_NDM-5 transformants positive for a longer insert (i.e. blaOXA-181) by the total number of transformants.

Nucleotide sequence accession numbers

The nucleotide sequence of strain JBEHAAB-19-0176 was deposited in DDBJ/ENA/GenBank under the BioSample accession number: SAMD00502403.

Results and discussion

An isolate co-carrying blaNDM-5 and blaOXA-181, JBEHAAB-19-0176, was obtained from the faeces of an outpatient who had returned from a trip to Bangladesh and was diagnosed with Campylobacter enteritis. These carbapenemase genes are the most common and the third most common carbapenemase genes encountered in a hospital surveillance in Bangladesh, respectively;22 however, both genes have rarely been detected in Japan.12 Therefore, it is likely that the isolate was imported from Bangladesh. This strain belonged to ST648, a globally disseminated carbapenem-resistant clone and the first reported isolate with NDM-5 in the UK in 2011.23 Complete sequences showed that the isolates had a chromosome (5 224 427 bp) and five plasmids including the IncX3 plasmid pJBEHAAB-19-0176_NDM-OXA (63 152 bp) co-harbouring blaNDM-5, blaOXA-181 and a partial ColKP3 replicon. Details of the four other plasmids, namely pJBEHAAB-19-0176_1 (151 726 bp) harbouring IncFIB and IncFII replicons, pJBEHAAB-19-0176_2 (109 582 bp) harbouring IncFIB replicons, pJBEHAAB-19-0176_4 (59 260 bp) harbouring the IncI plasmid, and pJBEHAAB-19-0176_5 (3066 bp) harbouring the Col replicon, are listed in Table 1. The strain harboured resistance genes against aminoglycosides (aadA2), sulphonamides (sul1) and trimethoprim (dfrA12) on the chromosome. The quinolone resistance gene qnrS1 was also located on the IncX3 plasmid, together with blaNDM-5 and blaOXA-181, whereas the macrolide resistance gene mph(A) was present on both the chromosome and IncFIB plasmids. Additionally, blaCMY-141 was located on an IncI (gamma) plasmid.

Table 1. Genomic information of E. coli JBEHAAB-19-0176

Replicon	Length (bp)	Inc type	Acquired antimicrobial resistance gene(s)	GenBank
accession no.	
Chromosome	5 224 427	ND	mph(A), sul1, aadA2, dfrA12, qepA4	AP028869	
pJBEHAAB-19-0176_NDM-OXA	63 152	IncX3, ColKP3	bla NDM-5, blaOXA-181, qnrS1	AP028870	
pJBEHAAB-19-0176_1	151 726	IncFIB, IncFII	mph(A), erm(B)	AP028871	
pJBEHAAB-19-0176_2	109 582	IncFIB	ND	AP028872	
pJBEHAAB-19-0176_4	59 260	IncI(gamma)	bla CMY-141	AP028873	
pJBEHAAB-19-0176_5	3066	Col	ND	AP028874	
ND, not detected.

IncX3 plasmids harbouring blaNDM-5 or blaOXA-181 have been frequently reported worldwide.3–7 We compared pJBEHAAB-19-0176_NDM-OXA with two representative IncX3 plasmids, a typical plasmid harbouring blaNDM-5 (pJBBDAGF-19-0019_NDM-5) and another harbouring blaOXA-181 (pJBCDAAC-19-0068_OXA-181) (Figure 1a), which were identified in the JBBDAGF-19-0019 and JBCDAAC-19-0068 E. coli isolates, respectively, in the JARBS-GNR surveillance.12 The three plasmids shared backbone genetic structures of the IncX3 plasmid, such as replication and conjugal transfer genes, which were 99.9% identical among the plasmids. A conjugation assay using a recipient strain E. coli ML4909 confirmed that pJBEHAAB-19-0176_NDM-OXA could be transferable with a transfer efficiency of 3.1 × 10−2, which was comparable to that of pJBBDAGF-19-0019_NDM-5 (1.4 × 10−2) and noticeably higher than that of pJBCDAAC-19-0068_OXA-181 (2.3 × 10−4). In addition to the IncX3 backbone, pJBEHAAB-19-0176_NDM-OXA had two genetic regions harbouring blaOXA-181 and blaNDM-5 (Figure 1a, i and ii), which were homologous to those found in pJBBDAGF-19-0019_NDM-5 and pJBCDAAC-19-0068_OXA-181, respectively. Among the three plasmids, two IS3000 sequences and Tn5403 were specific to pJBEHAAB-19-0176_NDM-OXA, whereas the other two plasmids had only one IS3000 sequence and lacked Tn5403. The pJBEHAAB-19-0176_NDM-OXA and pJBBDAGF-19-0019_NDM-5 plasmids had similar gene synteny with identical blaNDM-5 regions; however, the blaOXA-181 region was located in opposite orientations between the two IS3000s in pJBCDAAC-19-0068_OXA-181. There was a cluster of three mobile genetic elements, namely IS26, a putative Tn3-family transposase, and IS3000, located on the left boundary of the inverted blaOXA-181 region (Figure 1a, iii). The same cluster was found in pJBCDAAC-19-0068_OXA-181 (Figure 1a, iii’). However, the IS26 flanking sequences did not match between them, and the eight nucleotide sequences flanking the right side of IS26 of cluster iii matched that of another IS26 sequence on pJBCDAAC-19-0068_OXA-181 (flag 3 in Figure 1). These structural features suggest that pJBEHAAB-19-0176_NDM-OXA was formed through two recombination events (Figure 1b): (1) formation of a co-integrate between pJBBDAGF-19-0019_NDM-5 and pJBCDAAC-19-0068_OXA-181 via homologous recombination at IS3000 sequences; followed by (2) a second recombination event at IS26 sequences, resulting in excision of the IncX3 backbone of pJBCDAAC-19-0068_OXA-181. The Tn5403 transposon insertion into Tn2 could have occurred either after (Figure 1b, i) or before (Figure 1b, ii) these events. Both IS3000 and IS26 are known for their crucial roles in forming self-transmissible mobile genomic elements through a replicative mechanism and co-integrate formation during transposition.24–26 Our hypothesis assumes the coexistence of the two plasmids; however, it is unclear how they can coexist in a bacterial cell as they cannot be stably maintained due to plasmid incompatibility. One possible explanation is that an IncX3 plasmid carrying one gene was introduced into the bacteria along with the IncX3 plasmid carrying another gene via conjugal transfer, resulting in transient coexistence of the two plasmids. Considering the remarkably higher transfer efficiency of the plasmid with blaNDM-5 than that with blaOXA-181 found in this study (almost 100-fold), it is likely that the blaNDM-5-carrying plasmid was introduced into an organism already carrying blaOXA-181.

Figure 1. Genetic features of pJBEHAAB-19-0176_NDM-OXA and proposed models for its formation. (a) Linear comparison of the pJBEHAAB-19-0176_NDM-OXA plasmid with two other IncX3 plasmids harbouring blaNDM-5 or blaOXA-181. The coding sequences are represented by boxed arrows, and homologous regions (>99% identity) highlighted in grey shading. (b) Schematic diagram of pJBEHAAB-19-0176_NDM-OXA formation. The red line denotes an IncX3 plasmid with blaNDM-5, and the blue one denotes that with blaOXA-181. The red and orange arrows indicate IS3000, and the light blue arrows indicate IS26. Flags represent 8 bp flanking sequences of IS26. Wavy arrows show the recombination event.

Dual carriage of carbapenemase genes can contribute to high carbapenem resistance compared with that of organisms with single carriage.27 We generated the E. coli HST08 transformants carrying one of three IncX3 plasmids to investigate the phenotypic differences and potential benefits of co-carrying blaNDM-5 and blaOXA-181 compared with when either gene is carried alone. Antimicrobial susceptibility testing via a broth microdilution method showed that for the pJBEHAAB-19-0176_NDM-OXA transformants, the carbapenem MICs were identical to those of the pJBBDAGF-19-0019_NDM-5 transformants (Table 2). Interestingly, the MICs of ampicillin, piperacillin and cefpodoxime were 2-fold higher in the pJBEHAAB-19-0176_NDM-OXA transformants than in the pJBBDAGF-19-0019_NDM-5 transformants, indicating that the additional carriage of blaOXA-181 could increase the MICs of these antibiotics. Although pJBCDAAC-19-0068_OXA-181 showed 16 times lower carbapenem MICs than those of the other two transformants, it showed a >16 times higher MIC than that of the host (E. coli HST08).

Table 2. Antimicrobial susceptibility patterns of blaNDM-5- and/or blaOXA-181-harbouring isolates and transformants

Antimicrobial agents	MICs (mg/L)						
JBCDAAC-19-0068_OXA-181	JBEHAAB-19-0176_NDM-OXA	JBBDAGF-19-0019_NDM-5	HST:: pJBCDAAC-19-0068_OXA-181	HST:: pJBEHAAB-19-0176_NDM-OXA	HST:: pJBBDAGF-19-0019_NDM-5	HST08	
Meropenema	0.5	128	64	2	32	32	≤0.125	
Imipenema	2	64	64	2	32	32	≤0.125	
Ampicillina	2048	32 768	16 384	1024	8192	4096	≤4	
Piperacillina	2048	>4096	>4096	128	1024	512	1	
Cefpodoximea	2048	4096	2048	1	256	128	0.5	
Doripenem	≤0.5	>8	>8	1	>8	>8	≤0.5	
Ampicillin/sulbactam	32/16	>32/16	>32/16	32/16	>32/16	>32/16	≤4/2	
Piperacillin/tazobactam	64	>64	>64	8	64	>64	≤4	
Cefazolin	>16	>16	>16	16	>16	>16	≤1	
Cefoperazone/sulbactam	32/16	>32/16	>32/16	≤8/4	>32/16	>32/16	≤8/4	
Cefoxitin	≤8	>32	>32	8	>32	>32	≤8	
Cefmetazole	2	>32	>32	4	16	16	≤0.5	
Cefotetan	≤1	>32	>32	4	>32	>32	≤1	
Flomoxef	≤8	>32	>32	≤8	>32	>32	≤8	
Ceftazidime	16	>128	>128	≤0.5	>128	>128	≤0.5	
Ceftazidime/clavulanic acid	2/4	>32/4	>32/4	0.12/4	>32/4	>32/4	0.12/4	
Cefotaxime	>128	>128	>128	≤0.5	>128	>128	≤0.5	
Cefotaxime/clavulanic acid	>32/4	>32/4	>32/4	0.25/4	>32/4	>32/4	≤0.125/4	
Ceftriaxone	>64	>64	>64	≤0.5	>64	>64	≤0.5	
Cefepime	>32	>32	>32	≤1	32	>32	≤1	
Cefozopran	>16	>16	>16	2	>16	>16	≤1	
Aztreonam	64	64	2	≤0.5	≤0.5	≤0.5	≤0.5	
Gentamicin	≤1	2	>8	≤1	≤1	≤1	≤1	
Tobramycin	2	2	>8	≤1	≤1	≤1	≤1	
Amikacin	≤4	8	>32	≤4	≤4	≤4	≤4	
Levofloxacin	≤0.5	>8	>8	1	1	≤0.5	≤0.5	
Ciprofloxacin	0.5	>4	>4	1	0.5	≤0.25	≤0.25	
Minocycline	2	≤1	8	≤1	≤1	≤1	≤1	
Trimethoprim/sulfamethoxazole	>2/38	>2/38	>2/38	≤1/19	≤1/19	≤1/19	≤1/19	
Fosfomycin	16	16	>16	≤4	≤4	≤4	≤4	
Chloramphenicol	≤8	16	>16	≤8	≤8	≤8	≤8	
Colistin	≤1	≤1	≤1	≤1	≤1	≤1	≤1	
Tigecycline	≤0.25	≤0.25	≤0.25	≤0.25	≤0.25	≤0.25	≤0.25	
aMIC was determined using the broth microdilution method.

We further investigated the growth rates, i.e. the time to reach the exponential phase, of the two transformants through measuring changes in the OD600 over time in the presence of different antibiotics. Notably, the growth rate of the pJBEHAAB-19-0176_NDM-OXA transformant was higher than that of the pJBBDAGF-19-0019_NDM-5 transformant in the presence of piperacillin and cefpodoxime (Figure 2a and b). In contrast, both transformants showed comparable growth rates in the presence of meropenem and other β-lactams tested (Figure 2c and Figure S1, available as Supplementary data at JAC-AMR Online). The pJBEHAAB-19-0176_NDM-OXA transformants reached OD600 = 0.4 at 7.6 and 7.0 h after being inoculated in the presence of piperacillin and cefpodoxime, respectively. In contrast, the pJBBDAGF-19-0019_NDM-5 transformant took longer to reach OD600 = 0.4 (11.6 and 21.0 h, respectively) under the same conditions. These results suggest that the blaNDM-5/blaOXA-181 co-carrier had a growth advantage in the presence of these antibiotics, particularly cefpodoxime, compared with that of the blaNDM-5 single carrier. We tested this using a competition assay in which pJBEHAAB-19-0176_NDM-OXA and pJBBDAGF-19-0019_NDM-5 transformants were mixed and incubated in the presence of the above antibiotics. The bacteria were counted 24 h after inoculation (Figure 2d), and results showed that the growth of these transformants was almost comparable in the absence of any drugs and in the presence of meropenem. In contrast, the pJBEHAAB-19-0176_NDM-OXA transformant outcompeted the pJBBDAGF-19-0019_NDM-5 transformant in the presence of piperacillin and cefpodoxime, suggesting that the additional carriage of blaOXA-181 confers a competitive advantage over blaNDM-5 single carriers under these conditions.

Figure 2. The transformant carrying a plasmid co-harbouring blaNDM-5 and blaOXA-181 outcompetes one carrying a plasmid with only blaNDM-5 in the presence of piperacillin and cefpodoxime. Bacterial growth curves of pJBEHAAB-19-0176_NDM-OXA (HST::pJBEHAAB-19-0176_NDM-OXA, green) and pJBBDAGF-19-0019_NDM (HST::pJBBDAGF-19-0019_NDM-5, blue) transformants. OD600 was monitored in the presence of 64 mg/L piperacillin (a), 64 mg/L cefpodoxime (b) or 128 mg/L meropenem (c). The experiment was repeated three times. The pJBEHAAB-19-0176_NDM-OXA and pJBBDAGF-19-0019_NDM-5 transformants were mixed and cultured in the presence or absence of 128 mg/L piperacillin, 64 mg/L cefpodoxime and 16 mg/L meropenem. The bacterial counts were measured 24 h after inoculation, and the abundance ratio of pJBEHAAB-19-0176_NDM-OXA transformants calculated (d). The experiment was repeated eight times. Error bars represent standard deviations. Asterisks show significant differences (Steel test, P < 0.05).

To test this notion further, we compared the growth rate of pJBBDAGF-19-0019_NDM-5 transformants carrying an arabinose-inducible expression vector with (pBAD-OXA-181) or without (pBAD) blaOXA-181. When the experiment was performed with piperacillin and cefpodoxime, the growth rate of the transformant with pBAD-OXA-181 was faster than that with pBAD in the presence of arabinose (Figure 3a and b). Similar differences were not observed in the absence of arabinose (Figure 3d and e) and in the condition when meropenem was used instead of piperacillin/cefpodoxime (Figure 3c and f). In addition, the pJBBDAGF-19-0019_NDM-5 transformant carrying pBAD-OXA-181 outcompeted the one carrying pBAD in the presence of arabinose and piperacillin or cefpodoxime, but not meropenem (Figure 3g). Thus, the additional expression of blaOXA-181 was sufficient to confer a growth advantage in the presence of piperacillin and cefpodoxime. OXA-48 efficiently hydrolyses piperacillin, and OXA-181 shows similar substrate specificity.28,29 Therefore, it is probable that OXA-181 can hydrolyse piperacillin, and the additional presence of its encoding gene provides additional piperacillin hydrolytic activity to the bacteria. In contrast, it is unclear why additional carriage of blaOXA-181 provides a growth advantage in the presence of cefpodoxime because a subset of OXA-48 enzymes, including OXA-181, is incapable of hydrolysing extended-spectrum cephalosporins.6 This issue should be addressed in future studies.

Figure 3. Expression of blaOXA-181 confers a growth advantage to the transformant carrying blaNDM-5 in the presence of piperacillin and cefpodoxime. The OD600 of the pJBBDAGF-19-0019_NDM-5 transformant carrying the pBAD18-cm empty vector (HST::pJBBDAGF-19-0019_NDM-5::pBAD, blue) or the one harbouring blaOXA-181 (HST::pJBBDAGF-19-0019_NDM-5::pBAD-OXA, green) were monitored with (a–c) or without (d–f) 0.02% arabinose in the presence of 64 mg/L piperacillin (a, d), 64 mg/L cefpodoxime (b, e) or 128 mg/L meropenem (c, f). (g) The pJBBDAGF-19-0019_NDM-5 transformant carrying pBAD18-cm or the one harbouring blaOXA-181 (pBAD-OXA) were mixed and cultured in the presence or absence of 128 mg/L piperacillin, 64 mg/L cefpodoxime and 16 mg/L meropenem—0.02% arabinose was included in all test conditions. The abundance ratio of the transformants carrying pBAD-OXA was determined 24 h after inoculation. The experiment was repeated three times. Error bars represent standard deviations. Asterisks show significant differences (Tukey’s test, P < 0.05).

From the perspective of the molecular epidemiology of CPE, isolates co-carrying the two carbapenemase genes, blaNDM and blaOXA-48-like, have been found in various genotypes of Enterobacterales and have been increasingly reported worldwide.10,30–33 It was also observed that isolates co-carrying the two carbapenemase genes, blaNDM-1 and blaOXA-232 or blaOXA-181, emerged in three different STs of Klebsiella pneumoniae during a 4 year hospital surveillance, with a concomitant decrease in blaNDM-1-carrying isolates.34 The present study provides a clue to understanding the changing epidemiology of CPE. Piperacillin and cefpodoxime could be selective agents for the emergence of isolates co-carrying blaNDM and blaOXA-48-like genes. Meanwhile, it is conceivable that there are other advantages conferred by double carriage of these genes, which were not explored in this study. In addition, it should be noted that the changing epidemiology can result from many other factors, including characteristics of the host bacteria and plasmids carrying these genes, carriage of other antimicrobial resistance genes, and antimicrobial use.

In conclusion, we identified an IncX3 plasmid co-carrying dual carbapenemase genes in a national surveillance study using short- and long-read WGS. The plasmid was formed via the convergence of two highly disseminated plasmids in an endemic region and then imported to Japan. Due to the global spread of such problematic antimicrobial-resistant plasmids, continued surveillance and efforts to limit their further spread are warranted.

Supplementary Material

dlae073_Supplementary_Data

Acknowledgements

We thank Satoshi Nakano, Chika Arai, Shaheem Elahi, Noriko Sakamoto and Liansheng Yu for their technical assistance and helpful discussions. We also thank Yoshikazu Ishii for kindly providing a rifampicin-resistant E. coli strain.

Funding

This study was supported by the Japan Agency for Medical Research and Development (AMED) on Emerging and Re-emerging Infectious Diseases (grant numbers: JP19fk0108061 and 23fk0108604). This study was also supported by the Japan Society for the Promotion of Science (JSPS) KAKENHI (grant number: JP20K07496).

Transparency declarations

None to declare.

Supplementary data

Figure S1 is available as Supplementary data at JAC-AMR Online.
==== Refs
References

1 WHO . WHO publishes list of bacteria for which new antibiotics are urgently needed. 2017. https://www.who.int/news/item/27-02-2017-who-publishes-list-of-bacteria-for-which-new-antibiotics-are-urgently-needed.
2 Wu  W, Feng  Y, Tang  G  et al  NDM metallo-β-lactamases and their bacterial producers in health care settings. Clin Microbiol Rev  2019; 32 : e00115-18. 10.1128/CMR.00115-18 30700432
3 Zhu  W, Wang  X, Qin  J  et al  Dissemination and stability of the blaNDM-5-carrying IncX3-type plasmid among multiclonal Klebsiella pneumoniae isolates. mSphere  2020; 5 : e00917-20. 10.1128/mSphere.00917-20
4 Ho  P-L, Wang  Y, Liu  MC-J  et al  Incx3 epidemic plasmid carrying blaNDM-5 in Escherichia coli from swine in multiple geographic areas in China. Antimicrob Agents Chemother  2018; 62 : e02295-17. 10.1128/AAC.02295-17 29311058
5 Sugawara  Y, Akeda  Y, Hagiya  H  et al  Spreading patterns of NDM-producing Enterobacteriaceae in clinical and environmental settings in Yangon, Myanmar. Antimicrob Agents Chemother  2019; 63 : e01924-18. 10.1128/AAC.01924-18 30530602
6 Boyd  SE, Holmes  A, Peck  R  et al  OXA-48-like β-lactamases: global epidemiology, treatment options, and development pipeline. Antimicrob Agents Chemother  2022; 66 : e0021622. 10.1128/aac.00216-22 35856662
7 Pitout  JDD, Peirano  G, Kock  MM  et al  The global ascendency of OXA-48-type carbapenemases. Clin Microbiol Rev  2019; 33 : e00102-19. 10.1128/CMR.00102-19 31722889
8 Karlowsky  JA, Lob  SH, Kazmierczak  KM  et al  In vitro activity of imipenem against carbapenemase-positive Enterobacteriaceae isolates collected by the SMART global surveillance program from 2008 to 2014. J Clin Microbiol  2017; 55 : 1638–49. 10.1128/JCM.02316-16 28298454
9 Peirano  G, Chen  L, Nobrega  D  et al  Genomic epidemiology of global carbapenemase-producing Escherichia coli, 2015–2017. Emerg Infect Dis  2022; 28 : 924–31. 10.3201/eid2805.212535 35451367
10 Ham  DC, Mahon  G, Bhaurla  SK  et al  Gram-negative bacteria harbouring multiple carbapenemase genes, United States, 2012–2019. Emerg Infect Dis  2021; 27 : 2475–9. 10.3201/eid2709.210456 34424168
11 Li  Y, Qiu  Y, Fang  C  et al  Coexistence of blaOXA-58 and blaNDM-1 on a novel plasmid of GR59 from an Acinetobacter towneri isolate. Antimicrob Agents Chemother  2022; 66 : e0020622. 10.1128/aac.00206-22 35546112
12 Kayama  S, Yahara  K, Sugawara  Y  et al  National genomic surveillance integrating standardized quantitative susceptibility testing clarifies antimicrobial resistance in Enterobacterales. Nat Commun  2023; 14 : 8046. 10.1038/s41467-023-43516-4 38052776
13 Zankari  E, Hasman  H, Cosentino  S  et al  Identification of acquired antimicrobial resistance genes. J Antimicrob Chemother  2012; 67 : 2640–4. 10.1093/jac/dks261 22782487
14 Carattoli  A, Zankari  E, Garcia-Fernandez  A  et al  In silico detection and typing of plasmids using PlasmidFinder and plasmid multilocus sequence typing. Antimicrob Agents Chemother  2014; 58 : 3895–903. 10.1128/AAC.02412-14 24777092
15 Tanizawa  Y, Fujisawa  T, Nakamura  Y. DFAST: a flexible prokaryotic genome annotation pipeline for faster genome publication. Bioinformatics  2018; 34 : 1037–9. 10.1093/bioinformatics/btx713 29106469
16 Sullivan  MJ, Petty  NK, Beatson  SA. Easyfig: a genome comparison visualizer. Bioinformatics  2011; 27 : 1009–10. 10.1093/bioinformatics/btr039 21278367
17 Siguier  P, Perochon  J, Lestrade  L  et al  ISfinder: the reference centre for bacterial insertion sequences. Nucleic Acids Res  2006; 34 : D32–6. 10.1093/nar/gkj014 16381877
18 Ma  L, Ishii  Y, Ishiguro  M  et al  Cloning and sequencing of the gene encoding Toho-2, a class A β-lactamase preferentially inhibited by tazobactam. Antimicrob Agents Chemother  1998; 42 : 1181–6. 10.1128/AAC.42.5.1181 9593147
19 Liu  Z, Wang  Z, Lu  X  et al  Structural diversity, fitness cost, and stability of a blaNDM-1-bearing co-integrate plasmid in Klebsiella pneumoniae and Escherichia coli. Microorganisms  2021; 9 : 2435. 10.3390/microorganisms9122435 34946035
20 Guzman  LM, Belin  D, Carson  MJ  et al  Tight regulation, modulation, and high-level expression by vectors containing the arabinose PBAD promoter. J Bacteriol  1995; 177 : 4121–30. 10.1128/jb.177.14.4121-4130.1995 7608087
21 Shanmugakani  RK, Sugawara  Y, Akeda  Y  et al  blaOXA-731, a new chromosome-encoded blaOXA-48-like variant in Shewanella sp. from the aquatic environment in Myanmar. Environ Microbiol Rep  2020; 12 : 548–54. 10.1111/1758-2229.12875 32776437
22 Farzana  R, Jones  LS, Rahman  MA  et al  Genomic insights into the mechanism of carbapenem resistance dissemination in Enterobacterales from a tertiary public heath setting in South Asia. Clin Infect Dis  2023; 76 : 119–33. 10.1093/cid/ciac287 35412593
23 Hornsey  M, Phee  L, Wareham  DW. A novel variant, NDM-5, of the New Delhi metallo-β-lactamase in a multidrug-resistant Escherichia coli ST648 isolate recovered from a patient in the United Kingdom. Antimicrob Agents Chemother  2011; 55 : 5952–4. 10.1128/AAC.05108-11 21930874
24 He  S, Hickman  AB, Varani  AM  et al  Insertion sequence IS26 reorganizes plasmids in clinically isolated multidrug-resistant bacteria by replicative transposition. mBio  2015; 6 : e00762. 10.1128/mBio.00762-15 26060276
25 Sun  J, Yang  RS, Zhang  Q  et al  Co-transfer of blaNDM-5 and mcr-1 by an IncX3–X4 hybrid plasmid in Escherichia coli. Nat Microbiol  2016; 1 : 16176. 10.1038/nmicrobiol.2016.176 27668643
26 Partridge  SR, Kwong  SM, Firth  N  et al  Mobile genetic elements associated with antimicrobial resistance. Clin Microbiol Rev  2018; 31 : e00088-17. 10.1128/CMR.00088-17 30068738
27 Li  X, He  J, Jiang  Y  et al  Genetic characterisation and passage instability of a hybrid plasmid co-harbouring blaIMP-4 and blaNDM-1 reveal the contribution of insertion sequences during plasmid formation and evolution. Microbiol Spectr  2021; 9 : e0157721. 10.1128/Spectrum.01577-21 34908434
28 Dortet  L, Oueslati  S, Jeannot  K  et al  Genetic and biochemical characterisation of OXA-405, an OXA-48-type extended-spectrum β-lactamase without significant carbapenemase activity. Antimicrob Agents Chemother  2015; 59 : 3823–8. 10.1128/AAC.05058-14 25870062
29 Oueslati  S, Nordmann  P, Poirel  L. Heterogeneous hydrolytic features for OXA-48-like β-lactamases. J Antimicrob Chemother  2015; 70 : 1059–63. 10.1093/jac/dku524 25583748
30 Hopkins  KL, Ellaby  N, Ellington  MJ  et al  Diversity of carbapenemase-producing Enterobacterales in England as revealed by whole-genome sequencing of isolates referred to a national reference laboratory over a 30-month period. J Med Microbiol  2022; 71 : 001518. 10.1099/jmm.0.001518
31 Castanheira  M, Doyle  TB, Collingsworth  TD  et al  Increasing frequency of OXA-48-producing Enterobacterales worldwide and activity of ceftazidime/avibactam, meropenem/vaborbactam and comparators against these isolates. J Antimicrob Chemother  2021; 76 : 3125–34. 10.1093/jac/dkab306 34459890
32 Sader  HS, Mendes  RE, Carvalhaes  CG  et al  Changing epidemiology of carbapenemases among carbapenem-resistant Enterobacterales from United States hospitals and the activity of aztreonam-avibactam against contemporary Enterobacterales (2019–2021). Open Forum Infect Dis  2023; 10 : ofad046. 10.1093/ofid/ofad046
33 Estabrook  M, Muyldermans  A, Sahm  D  et al  Epidemiology of resistance determinants identified in meropenem-nonsusceptible Enterobacterales collected as part of a global surveillance study, 2018 to 2019. Antimicrob Agents Chemother  2023; 67 : e0140622. 10.1128/aac.01406-22 37074173
34 Sakamoto  N, Laolerd  W, Akeda  Y  et al  Temporal shifts in the predominant carbapenemase gene types among carbapenemase-producing Klebsiella pneumoniae isolated in Bangkok, Thailand, during 2013–2016. J Med Microbiol  2023; 72 : 001711. 10.1099/jmm.0.001711
