An antiplasmid system drives antibiotic resistance gene integration in carbapenemase-producing Escherichia coli lineages

Zongo, Pengdbamba Dieudonné; Cabanel, Nicolas; Royer, Guilhem; Depardieu, Florence; Hartmann, Alain; Naas, Thierry; Glaser, Philippe; Rosinski-Chupin, Isabelle

doi:10.1038/s41467-024-48219-y

Download PDF

Article
Open access
Published: 15 May 2024

An antiplasmid system drives antibiotic resistance gene integration in carbapenemase-producing Escherichia coli lineages

Nature Communications volume 15, Article number: 4093 (2024) Cite this article

2078 Accesses
27 Altmetric
Metrics details

Subjects

Abstract

Plasmids carrying antibiotic resistance genes (ARG) are the main mechanism of resistance dissemination in Enterobacterales. However, the fitness-resistance trade-off may result in their elimination. Chromosomal integration of ARGs preserves resistance advantage while relieving the selective pressure for keeping costly plasmids. In some bacterial lineages, such as carbapenemase producing sequence type ST38 Escherichia coli, most ARGs are chromosomally integrated. Here we reproduce by experimental evolution the mobilisation of the carbapenemase bla_OXA-48 gene from the pOXA-48 plasmid into the chromosome. We demonstrate that this integration depends on a plasmid-induced fitness cost, a mobile genetic structure embedding the ARG and a novel antiplasmid system ApsAB actively involved in pOXA-48 destabilization. We show that ApsAB targets high and low-copy number plasmids. ApsAB combines a nuclease/helicase protein and a novel type of Argonaute-like protein. It belongs to a family of defense systems broadly distributed among bacteria, which might have a strong ecological impact on plasmid diffusion.

A novel antibiotic class targeting the lipopolysaccharide transporter

Article Open access 03 January 2024

A Gram-negative-selective antibiotic that spares the gut microbiome

Article 29 May 2024

Molecular mechanisms of antibiotic resistance revisited

Article 21 November 2022

Introduction

Antibiotic resistance is of major public health concern, having been responsible for more than 1.2 million deaths in 2019¹. Mobile genetic elements, particularly plasmids carrying antibiotic resistance genes (ARG), are key contributors to antibiotic resistance spread within and between species². While often carrying features beneficial under specific adverse environments, plasmids frequently confer a fitness cost to their bacterial host^3,4. This is expected to lead to their loss over time and, under positive selection, to the integration of beneficial traits into the host chromosome as predicted by Bergstrom et al.’ mathematical model of plasmid persistence⁵. However, some plasmids may persist over a long period even in the absence of selective pressure, a phenomenon known as the plasmid paradox^6,7. Solutions to this paradox include compensatory mutations improving the fitness of plasmid-carrying bacteria^8,9,10 and high conjugation rates¹¹. On the other hand, bacteria have developed defense mechanisms to counteract invasion by parasitic DNAs, some of them able to limit plasmid diffusion^12,13,14.

Carbapenemases are enzymes responsible for the hydrolysis of carbapenems, broad-spectrum antibiotics of major clinical importance. In Western Europe, the class D enzyme OXA-48 is the most frequent carbapenemase in Enterobacterales^15,16,17. This is generally attributed to the high conjugation rate of the IncL group pOXA-48 plasmids¹⁸. These 56–76 kb-long plasmids induce variable fitness costs to their hosts¹⁹. In addition, OXA-48 is also encoded on the chromosome^15,20,21, mainly in E. coli isolates of sequence type ST38, a major cause of extra-intestinal infections²². These ST38 isolates belong to disseminated clones with the bla_OXA-48 gene embedded in a complete or partial composite transposon, Tn6237, that likely derives from a pOXA-48 plasmid.

The factors driving ARG integration into the chromosome are largely unknown. Here, to identify such factors and to decipher the biological mechanisms leading to the selection of ARG integration, we experimentally characterize the move of the carbapenemase bla_OXA-48 gene from a pOXA-48 plasmid into the chromosome. We followed the integration of bla_OXA-48 through experimental evolution and show that the emergence of ST38 lineages with integrated bla_OXA-48 depends on the fitness cost imposed by pOXA-48 and on its elimination by a novel antiplasmid system.

Results

pOXA-48 plasmids induce a fitness cost and are unstable in ST38 E. coli

To set up an experimental model to study bla_OXA-48 integration we used E. coli ST38 because of its evolutionary history characterized by pervasive ARG integrations^23,24,25,26 in addition to its clinical relevance. We chose three ST38 isolates from water or sewage that did not carry bla_OXA-48 (Supplementary Data 1) and transferred into these isolates three variants of the IncL pOXA-48 plasmid (pOXA-48_1, pOXA-48_2 and pOXA-48_3) by conjugation. The three isolates differed by their position in ST38 phylogenetic tree (Supplementary Information Supplementary Fig. 1). In particular, ST38-1 belonged to a branch of the tree also encompassing three sub-clades characterized by different bla_OXA-48 chromosomal integrations¹⁵ (Supplementary Information Supplementary Fig. 1). It is closely related to one of these sub-clades, colored in blue on the tree, sharing the same H-and 0-antigens (O2-H30). Plasmids were scarce in the three selected isolates. ST38_1 and ST38_3 contained only an IncFII and an IncFIC(FII) plasmid, respectively, whereas ST38_2 had no plasmid. Their complete genome sequence showed that, in ST38_1 and ST38_2, all ARGs were chromosomally inserted. In contrast, in ST38_3 most ARGs were carried by the IncFIC(FII) plasmid (Supplementary Data 1). We used as donor strains of the three pOXA-48 plasmids three clinical Klebsiella pneumoniae isolates. K. pneumoniae is the most common pOXA-48 plasmid-bearing species found in hospitals²⁷. The three plasmids mainly differed by the presence and position of IS1 generating or not a composite transposon, Tn6237 (Fig. 1a). This 21.5 kb transposon consists of two IS1 bracketing bla_OXA-48 and other plasmid genes. While this structure was present in pOXA-48_1, it was inactivated in pOXA48-2 by a 6.8 kb-long deletion, including the bla_OXA-48-distant IS1. In pOXA-48_3, the 4.5 kb region including bla_OXA-48 and the bla_OXA-48-proximal IS1 was inverted, leaving bla_OXA-48 outside the composite transposon (Fig. 1a). We first showed that the plasmid structure did not influence the meropenem minimum inhibitory concentration (MIC) of the transconjugant unlike the genetic background of the recipient strain with a higher MIC for ST38-1 transconjugants, 0.5 µg ml-1, compared to 0.25 µg ml-1 and 0.38 µg ml-1 for ST38_2 and ST38_3 transconjugants respectively (Supplementary Information Supplementary Table 1).

**Fig. 1: pOXA-48 plasmids induce variable fitness costs and are unstable in three ST38 strains.**

To assess the fitness cost induced by pOXA-48 plasmids we compared the maximum growth rate of the transconjugants and isogenic plasmid-free strains. We estimated the relative fitness as the ratio of doubling time of transconjugants over plasmid-free strain. At least 15% doubling time increase was induced in lysogeny broth (LB) medium by pOXA-48 plasmids in ST38_1 (Fig. 1b). In contrast, in this medium, the three plasmids induced a lower ( < 5%) fitness cost in ST38_2 and ST38_3 (Fig. 1b). Fitness cost was dependent on the growth medium as in M9 glucose minimal medium it rose up to 12% for ST38_2 transconjugants, while decreasing to 7–11% in ST38_1 transconjugants.

As pOXA-48 plasmids often conferred a fitness cost to their host, we quantified their stability in the three ST38 strains, following ten-day serial passages of the transconjugants (ca. 80 generations) in LB in the absence of antibiotic. The three pOXA-48 plasmids were gradually lost and at day 10, more than 95% of the population in ST38_1 and 40–80% in ST38_2 and ST38_3 had lost pOXA-48, irrespective of the plasmid variant (Fig. 1c).

Twenty-eight days of in vitro evolution led to frequent bla _OXA-48 chromosomal integration in ST38_1

To determine whether the fitness of pOXA-48 transconjugants could be improved by plasmid-host coevolution, we first performed 28-day experimental evolutions of ST38_1/pOXA-48_1 transconjugants (Table 1). We hypothesized that this strain/plasmid combination might lead to chromosomal integration of bla_OXA-48 and pOXA-48_1 loss. Indeed, ST38_1 was phylogenetically close to sub-clades with chromosomal bla_OXA-48 (Supplementary Information Supplementary Fig. 1), the Tn6237 structure in pOXA-48_1 had the characteristics of a mobile composite transposon and pOXA-48_1 induced a high fitness cost in ST38_1. Given the instability of pOXA-48_1, the growth medium was supplemented with subinhibitory concentration (0.1 µg.ml^-1) of meropenem every passage or every third passage (Table 1 and Supplementary Information Supplementary Table 2). The ST38_1 plasmid-free strain was evolved as control. Five independent lineages were derived for each condition. Experimental evolution led to a rapid growth rate improvement for all transconjugant lineages under both meropenem conditions (Fig. 2a). To characterize the populations, we first performed pool-sequencing after 7, 14, 21 and 28 days of evolution. We observed a progressive decrease in read coverage of the pOXA-48_1 DNA-region outside Tn6237 (Fig. 2b, Supplementary Information Supplementary Fig. 2), suggesting the enrichment of bacteria having lost pOXA-48_1 but keeping bla_OXA-48.

Table 1 Experimental scheme of the different evolution experiments

Full size table

**Fig. 2: Experimental evolution of pOXA-48_1 transconjugants selects Tn6237 transposition and pOXA-48 loss in ST38_1.**

To further characterize putative transposition events of Tn6237, we PCR-screened meropenem-resistant colonies isolated at 28-day evolution for bla_OXA-48 and plasmid origin of replication (repA). We identified bla_OXA-48⁺ repA^- colonies in all ST38_1/pOXA-48_1 lineages (Supplementary Data 2). Two to three of these colonies per lineage were whole-genome-sequenced (WGS). Variant analysis revealed IS1 new junctions in the chromosome or in the ST38_1 IncFII plasmid. We verified by PCR and Sanger sequencing that these new junctions corresponded to Tn6237 transposition. In total, ten integrations at different positions in the chromosome and three in the IncFII plasmid were identified (Supplementary Data 3 and Fig. 2c). Analysis of pool sequencing reads revealed 44 additional new junctions, considering a threshold of 5% of the reads (Supplementary Data 4). Twenty were chosen for PCR analysis among which 16 were confirmed as Tn6237 integration. Integrations were detected already at seven days of evolution and multiple integrations at multiple positions on the chromosome and the IncFII plasmid were selected in each lineage after 28 days (Supplementary Data 4 and Fig. 2c).

To determine whether bla_OXA-48 integration was associated with a fitness improvement, we selected three colonies, one in the chromosome and two in the IncFII plasmid. The growth rate of the three clones was increased with a 10–15% reduction in doubling time compared to the original transconjugant (Fig. 2d) and the meropenem MIC decreased by 50% (0.25 vs 0.5 µg ml-1) (Supplementary Information Supplementary Table 1). Therefore, in ST38_1, bla_OXA-48 integration appears as a major mechanism to relieve the fitness cost associated with pOXA-48_1 carriage.

bla _OXA-48 integration depends on pOXA-48 structure and on recipient strain

To determine whether bla_OXA-48 integration during experimental evolution was dependent on the presence of Tn6237, we then evolved ST38_1 transconjugants carrying plasmids pOXA-48_1, pOXA-48_2 or pOXA-48_3 in LB with meropenem added every third passage (Table 1). In contrast to pOXA-48_1, no fitness improvement throughout time was observed for pOXA-48_2, while the fitness of the population was slightly improved for pOXA-48_3 transconjugants at 28-day evolution (Fig. 3a). Bulk DNA sequencing of the whole populations at day 21 and at day 28 revealed after 28-day evolution only a relative decrease in read-coverage of the pOXA-48_3 region equivalent to the one lost in pOXA-48_1 (Supplementary Information Supplementary Fig. 3a, b) and a few new IS1 junctions. Two Tn6237 integrations in the chromosome were confirmed by PCR, suggesting the reconstitution of Tn6237 in pOXA-48_3 (Supplementary Data 4). This showed that bla_OXA-48 integration was dependent on its inclusion in Tn6237.

**Fig. 3: Experimental evolution of pOXA-48_transconjugants depends on plasmid structure and plasmid-induced fitness cost.**

To determine whether the result of the experimental evolution was dependent on the ST38_1/pOXA-48_1 combination, we performed similar experimental evolution of ST38_2/pOXA-48_1 and ST38_3/pOXA-48_1 transconjugants (Table 1). For both experiments, we added meropenem every third passage in LB. For ST38_2/pOXA-48_1 we also checked that a daily addition had no influence on the final result. In all cases, no evolution of the fitness was observed (Fig. 3b) while PCR-screening of meropenem-resistant colonies at 28-day did not identify any bla_OXA-48⁺ repA^- colonies (Supplementary Data 2). This suggested that in the absence of a sufficient fitness cost, bla_OXA-48 integrations were not enriched enough to be detected. To test this hypothesis, we performed a similar experiment in minimal medium in which pOXA-48_1 induced a higher fitness cost to the ST38_2 transconjugants (Fig. 1b). A rapid improvement of growth rate was observed (Fig. 3c). We detected a few bla_OXA-48⁺ repA^- colonies, associated with a Tn6237 chromosomal integration, in only two of the five evolved lineages by PCR screening (Supplementary Data 2 and Supplementary Data 3). These colonies showed a full fitness recovery (Supplementary Information Supplementary Fig. 3c). Pool sequencing did not reveal any integration site at a frequency superior to 5%, nor a significant plasmid loss (Supplementary Information Supplementary Fig. 3d). Therefore, chromosomal integration and plasmid loss were not the main contributors to the fitness recovery observed in ST38_2 transconjugants, in contrast to ST38_1 transconjugants.

pOXA-48 plasmids are stabilized by inactivation or mutation of a novel antiplasmid system

In all evolved lineages, including those where Tn6237 transposition was selected, bacteria still retaining the plasmid could be found after 28-day evolution. To identify other possible paths for plasmid-host coadaptation, we performed WGS of three to five bla_OXA-48⁺ repA⁺ colonies of each evolved lineage. We identified sporadic mutations occurring in the three ST38 evolved transconjugants and control plasmid-free strains, including only nine in different loci of pOXA-48 (eight in pOXA-48_1 and one in pOXA-48_2) (Supplementary Data 5). Mutations potentially decreasing susceptibility to carbapenems (in ompC or envZ for instance) were encountered. Large chromosomal deletions ranging from 15 to 40 kb and encompassing mutS and rpoS were also observed in evolved ST38_1 lineages. They were independent on pOXA-48 presence and meropenem addition, as they also occurred in control experiments. mutS loss led to a large number of mutations (n = 13 to 38) likely resulting from an hypermutator phenotype. On the other hand, in 62% (86/138) of evolved transconjugants, we observed convergent evolution with mutations in two adjacent genes, F3141 or F3140, located in a ST38_1 genomic island (Fig. 4a) or its ortholog in ST38_2 (Supplementary Data 6). The most frequent mutations were IS1 insertions (n = 67 with 37 different insertion sites). We also detected insertions of other IS (n = 3), non-synonymous (n = 3, two different), non-sense (n = 2), or frameshift (n = 7, four different) mutations and complete or partial deletions (n = 4) of these genes (Supplementary Data 6). All (15/15) sequenced bla_OXA-48⁺ repA⁺ colonies from the M9-evolved ST38_2/pOXA-48_1 transconjugants carried mutations in orthologs of F3141 or F3140 and all (n = 30) from the five LB-evolved lineages shared the same IS1 insertion in F3140 suggesting that it was present but not detectable in sequencing reads of the initial culture (Supplementary Data 6). No mutation in F3141-F3140 orthologs was detected in ST38_3/pOXA-48_1 transconjugants evolved in LB.

**Fig. 4: Mutations in the F3141-F3140 operon led to pOXA-48 stabilization.**

In ST38_1, six out of the ten tested mutations led to a fitness improvement compared to the original transconjugant, which was variable and always lower than following bla_OXA-48 integration and pOXA-48_1 loss (Fig. 4b). In ST38_2 the six tested mutations in F3141-F3140 orthologs led, in M9, to a full fitness recovery (Fig. 4b). Three ST38_1 and three ST38_2 evolved transconjugants mutated in F3140 or in F3141 were tested for pOXA-48_1 plasmid stability by 10-day serial passages in the absence of meropenem (Fig. 4c). In all tested transconjugants an increase in plasmid persistence was observed compared to the original transconjugant (Fig. 4c). pOXA-48_1 was kept in more than 90% of ST38-2 mutated transconjugants. In ST38-1 the stability was more variable with 50 to 100% bacteria keeping pOXA-48_1.

Conjugation of the pOXA-48 plasmids to genetically modified derivatives of ST38_1 and ST38_3 deleted for F3141-F3140 confirmed the stabilization of the plasmid in the absence of F3141-F3140 (Fig. 4d). Of note, ST38_2 was not amenable to genetic manipulation under the conditions we used. Furthermore, complementation in ST38_1∆F3141-F3140 strain with a pHV7-derivative plasmid expressing F3141-F3140 under the control of the arabinose inducible promoter led to more than 80% loss of any of the three pOXA-48 plasmids at 48 h following induction (Fig. 4e). Induced expression of the operon with a non-sense mutation in either F3141 or in F3140 did not destabilize pOXA-48_1 confirming that the two proteins were needed for plasmid loss (Fig. 4f).

To determine whether F3141-F3140 might be involved in a more global antiplasmid activity, we compared the stability of five plasmids with different replication origins and copy numbers. In addition to the IncL pOXA-48 plasmids, F3141-F3140 deletion increased the stability of p15A, ColE1 and pMB1-type plasmids, but had no effect on pSC101 and IncFII/IncFIB plasmid stability (Fig. 5a). This confirmed that F3141-F3140 corresponds to a novel antiplasmid defense system that we renamed apsAB for antiplasmid system AB.

**Fig. 5: F3141-F3140 (ApsAB) acts as an antiplasmid system active against low and high copy number plasmids.**

We found that ApsAB also reduces the conjugation frequency of the three pOXA-48 plasmids while comparing the transfer frequency with ST38_1 wild type or ST38_1 ∆apsAB as recipient strains (Fig. 5b). This suggests that ApsAB interferes with foreign DNA acquisition. However, ApsAB did not affect the transformation frequency of p15A and ColE1 plasmids (Fig. 5c). As some antiplasmid systems are also involved in antiphage defense, we tested the activity of ApsAB chromosomally expressed in MG1655 against eight different phages from E. coli (lambda, T4, P1, 186cIts, CLB_P2, LF82_P8, T5 and T7)²⁸ and did not observe any antiphage activity (Supplementary Information Supplementary Fig. 4).

Finally, to determine whether ApsAB actively eliminates plasmids, we quantified the persistence of a ColE1 plasmid over time following the induction of apsAB expression from a chromosomally integrated copy. Plasmid elimination was observed from between 2h30 and 3 h of arabinose induction and was almost complete ( > 95%) by a six-hour induction, whether or not chloramphenicol was added at 2h30 to stop cell division. These results indicate that ApsAB actively eliminates the plasmid, likely by promoting its degradation (Fig. 5d).

ApsAB defense system is necessary for the selection of bla _OXA-48 integration

We hypothesized that the elimination of pOXA-48 plasmids by ApsAB could contribute to the emergence of lineages with bla_OXA-48 inserted in the chromosome. To test this hypothesis, we deleted the apsAB operon in the original ST38_1/pOXA-48_1 transconjugant, which has the capacity to rapidly evolve towards Tn6237 integration. No significant difference in fitness was detected between the deleted strain and the original transconjugant (Fig. 6a). As expected, apsAB deletion led to pOXA-48_1 stabilization (Fig. 6b). We then performed experimental evolution in LB medium of five lineages of ST38_1∆apsAB/pOXA-48_1 (Table 1). Contrary to the original ST38_1/pOXA-48_1 (Fig. 2a), no fitness recovery was detected after 28 days of evolution of ST38_1∆apsAB/pOXA-48_1 (Fig. 6c). No bla_OXA_-₄₈⁺ repA^- colonies out of 120 tested colonies at day 28 was detected by PCR screening. Similarly, pool sequencing of the whole population at days 7, 14, 21, and 28 showed no loss of pOXA-48_1 plasmid by analyzing plasmid read-coverage and no Tn6237 transposition based on the PCR testing of the few new IS1 junctions (n = 4) (Fig. 6d). Our results, therefore, show that plasmid destabilization through the activity of ApsAB is a key factor in the emergence of lineages with chromosomally integrated bla_OXA-48.

**Fig. 6: *bla*_OXA-48 integration and pOXA-48_1 loss require the *apsAB* operon.**

ApsAB is the first characterized member of a broad family of Argonaute-like systems

Sequence similarity search against the NCBI nr public database by BLASTP revealed positive matches of ApsA and ApsB only with proteins of unknown function. Nevertheless, using similarity search based on structure predictions, we predicted in ApsA a central helicase domain with conserved residues characteristic of superfamily 2 helicase²⁹ and a C-terminal domain containing a PD-(D/E)XK- superfamily nuclease motif ³⁰ (Fig. 7a). Directed mutagenesis showed that substitutions predicted to impair ATP hydrolysis (E533A) and helicase activity (K221A) completely abolished pOXA-48 plasmids destabilization while a mutation in the predicted nuclease active site (K1435A) had a partial effect (Fig. 7b and Supplementary Information Supplementary Fig. 5). On the other hand, structural modeling of ApsB revealed a loose structural similarity with prokaryotic Argonaute proteins (pAgos) acting as nucleic acid-guided endonucleases. However, ApsB lacked typical PIWI and PAZ domains characteristic of Argonautes and was not identified among pAgos³¹ (Fig. 7c). Sequence alignment of ApsB homologs retrieved from the NCBI database identified a conserved Y[X]₃K[X]_nQG[X]_nK motif (Fig. 7d, Supplementary Information Supplementary Fig. 6 and Supplementary Data 7), reminiscent of the conserved residues Y[X]₃K[X]_nQ[X]_nK characteristic of the MID domain of long-B pAgos that interacts with the 5’-end of guide nucleic acids³². The K413A substitution in this motif completely abolished pOXA-48 plasmid destabilization, supporting the requirement of this motif for ApsAB antiplasmid activity (Fig. 7b).

**Fig. 7: ApsAB associates a protein with helicase and nuclease domains and a novel Argonaute-like protein.**

The association of an Argonaute-like protein and a protein with helicase and nuclease domains was reminiscent of the DdmDE plasmid defense system recently characterized in Vibrio cholerae¹³. No sequence similarity could be found between ApsB and DdmE. Their predicted 3D structures were different (Root-Mean-Square Deviation of atomic positions (RMSD) = 30.6) and DdmE did not contain the MID-like motif. However, searching for ApsA homologs by PSI-BLAST in databases unveiled a wide family of proteins ranging from ApsA-like to DdmD-like proteins. (Fig. 8a, Supplementary Data 8).

**Fig. 8: ApsAB belongs to a broad family of putative antiplasmid systems.**

We identified ApsAB-like systems mainly among Enterobacterales but also in other gamma-proteobacteria and some beta-proteobacteria and cyanobacteria (Fig. 8a and Supplementary Data 7, 8). In E. coli, complete or partially deleted apsAB homologs were located in at least three different families of genomic islands, some of them encoding other defense systems like Shango systems (Fig. 8b). We evaluated the antiplasmid activity of two representatives of E. coli ApsAB-like systems, showing 92/92 % or 28/25 % protein sequence identity with ST38_1 ApsA/ApsB respectively. Both systems were found to destabilize a ColE1 multicopy plasmid (Fig. 8c), indicating that the antiplasmid activity is not restricted to ApsAB from ST38 strains.

Discussion

Mathematical modeling predicts that selective pressures, such as those exerted on bacteria during antibiotic treatment, will promote the integration of beneficial plasmid genes into the chromosome⁵. In this study, by investigating chromosomal integrations of bla_OXA-48, we demonstrate that, in addition to the selective advantage and a genetic structure promoting the mobility of the ARG, integration also depends upon a bacterial host antiplasmid system (ApsAB), here discovered. ApsAB influences plasmid-induced fitness cost and plasmid persistence. It belongs to a large family widespread across bacterial species ranging from cyanobacteria to Enterobacterales (Fig. 8a, Supplementary Data 7, 8).

To investigate the factors involved in bla_OXA-48 gene integration into another replicon, we performed experimental evolutions of pOXA-48s transconjugants of three E. coli ST38 strains. Subinhibitory concentration of carbapenem was added every passage or every three passages, leading to alternate selection for growth and for resistance (Table 1). Rapidly, for one of the five strain/plasmid pairs tested (ST38_1/pOXA48_1), bla_OXA-48 integrated lineages emerged and became dominant at day 28 (Fig. 2b). Tn6237 transposition from pOXA-48_1 appears as particularly efficient to move bla_OXA-48 into the chromosome or the resident F-plasmid, as in each evolved lineage multiple integration events were selected. bla_OXA-48 chromosomal integration was also observed with pOXA-48_3, but at late time points. These integration events were also due to Tn6237 transposition. We hypothesize that the Tn6237 structure has been reconstituted following homologous recombination between the two IS1999 copies surrounding bla_OXA-48 and the bla_OXA-48 proximal IS1 (Fig. 1a). Altogether the experimental evolution successfully reproduced bla_OXA-48 integrations observed in worldwide disseminated ST38 clinical lineages^20,21. Our results show that low doses of carbapenem, such as those encountered in the gut during parenteral administration of carbapenems³³, might have contributed to the emergence of bla_OXA-48 chromosomally integrated lineages.

In addition to bla_OXA-48 integration by transposition, we observed during experimental evolution convergent mutations (mainly inactivation) of the apsAB operon. ApsAB increases the fitness cost of pOXA-48 plasmids in ST38 transconjugants and destabilizes these plasmids (Fig. 4b, c). In experiments with ST38_1, the fitness gain linked to bla_OXA-48 integration and plasmid loss was high and lineages with bla_OXA-48 integration emerged rapidly. This was not the case with ST38_2 where both events led to equivalent fitness gain. This probably explains why apsAB inactivation events were generally selected at the expense of bla_OXA-48 chromosomal integration in this strain. More specifically, we showed that pOXA-48 plasmids destabilization was needed for the selection of bla_OXA-48 integration (Fig. 6d). apsAB location, embedded in a genomic island, and its distribution among ST38 isolates (Supplementary Information Supplementary Fig. 1) suggested that it has been horizontally acquired. This represents a typical case of MGE negative interactions between a genomic island and plasmids^34,35. apsAB was present in the three strains of the study and was associated to a low stability of the three pOXA-48 plasmids. However, while apsAB was frequently inactivated in colonies that retained pOXA-48 in ST38_1 and ST38_2, we did not detect inactivation events in ST38_3, suggesting that conservation of apsAB might be positively selected for an unknown reason in this genetic background. Nevertheless, apsAB deletion in this strain led to a stabilization of pOXA-48 plasmids (Fig. 4d). In addition to ApsAB, nine to eleven other defense systems were predicted in the three strains (Supplementary Data 1). None was specific to ST38_1 and we did not detect any mutation in these systems during experimental evolutions. Therefore, they probably did not influence the evolution results. apsAB is present in the four ST38 clinical lineages in which bla_OXA-48 is integrated in the chromosome (Supplementary Information Supplementary Fig. 1)¹⁵ and has likely driven the integration event leading to their dissemination. However, as ARGs integrated in the genome are frequently being reported in resistant E. coli and apsAB homologs only found in a subset of strains, other plasmid destabilizing systems are probably to be involved. Alternatively, in the absence of an active antiplasmid system, the mobilization to the chromosome might be much less frequently selected.

In recent years, a plethora of bacterial defense systems have been discovered. But, only a few of them, such as CRISPR, restriction-modification, pAgos, DdmDE and Wadjet were described to have antiplasmid activity^13,36,37. By combining different in silico search strategies, we uncovered a broad family of ApsA homologs that combine helicase and nuclease domains. Strikingly, this family was subdivided into two classes associated with two families of argonaute-like proteins, whose prototypes were ApsB and the recently described DdmE¹³. We therefore connected DdmDE- to ApsAB-like systems. However, unique to ApsB-like proteins is the conservation of a MID-like domain of long pAgos, which we showed to be essential for its activity. This suggests that ApsB is involved in a guide-dependent recognition of the target, ApsA bringing the nuclease activity, similarly to what has recently been described for a group of catalytically inactive long true pAgos³⁸. While these Argonaute systems confer immunity also via abortive infection³⁸, we do not have evidence that ApsAB kills plasmid invaded cells. Despite the absence of similarity between DdmE and ApsB, both DdmDE and ApsAB share an antiplasmid activity, like a second E. coli system, only 28/25% similar to ST38 ApsAB (Fig. 8c). This suggests that the compendium of systems we brought out represents a new diverse family of defense systems acting on plasmids. Their precise functional mechanisms, in particular the nature and specificity of the guides, remain to be characterized. These systems might also provide new tools for genomic engineering and plasmid curing, including those carrying ARGs.

To our knowledge, we provide the first evidence of an antiplasmid defense system influencing the evolutionary fate of an antibiotic resistance gene. Conversely, during our analysis of apsAB homologs, we retrieved many partially deleted systems, which might reflect the counter-selection of these systems in a context of antibiotic pressure that favors the stabilization and dissemination of ARG carrying plasmids.

Methods

Bacterial strains, growth conditions, and plasmids

Characteristics of the strains used in this work are summarized in Supplementary Data 1. All strains were submitted to Illumina sequencing and for some of them to long-read sequencing. Donor strains in conjugation experiments were three clinical K. pneumoniae isolates carrying either pOXA-48_1, pOXA-48_2 or pOXA-48_3. The organization of IS1 sequences in the three plasmids was confirmed by PCR. Recipient strains were three E. coli ST38 environmental isolates characterized by different O-types/H-types. Their position in the ST38 phylogeny is given in (Supplementary Information Supplementary Fig. 1a). ST38_1 carried an endogenous IncFII, 70.8 kb-long plasmid deprived of any ARG while ST38_3 contained a 161 kb-long IncFIC(FII) plasmid. ST38_2 has no plasmid (Supplementary Data 1). E. coli strains CNR36C9 (ST219) and CNR81D10 (ST10) encoding closely related (91%/92% amino acid (a.a.) identity) or distantly related (28%/25% a.a. identity) apsAB-like systems respectively were used for PCR-cloning of apsAB homologs. E. coli K-12 strain MG1655 was used for apsAB expression experiments from a mini-Tn7 inserted downstream glmS. E. coli strains DH5a and XL1 blue were used for cloning and MFD pir³⁹ for propagation of plasmids with RK6 origin of replication and for bacterial mating.

Bacterial growths were performed at 37 °C. Liquid cultures were performed in LB Miller or in M9 medium supplemented with glucose 0.4%, MgSO₄ 1 mM and CaCl₂ 0.1 mM with shaking. Where appropriate, antibiotics were added: meropenem (0.1 µg ml-1), apramycin (40 µg ml-1), zeocin (30 µg ml-1) and chloramphenicol (10 µg ml-1). Growth and mating experiments with MFDpir derivatives were performed in LB or LB agar complemented with 0.3 mM diaminopimelic acid (DAP, ThermoScientific). Induction from p_BAD promoter was performed in LB Miller supplemented with L-arabinose to a final concentration of 0.02 or 0.2% as indicated. E-test were performed on Mueller Hinton agar (MHA) medium. A list of plasmids used in the study and their main characteristics is provided in (Supplementary Data 9). All inserted fragments were verified by Sanger sequencing (Eurofins Genomics), and expression plasmids carrying apsAB were WGS (PlasmidSaurus or Eurofins Genomics).

Conjugation assay

Overnight precultures in LB of donors and recipient strains were diluted 1:100 into fresh LB and grown to an optical density at 600 nm (OD₆₀₀) of 0.6. After mixing at 1:1 ratio, 200 µl were spread on a filter (MILLIPORE type HAEP 0.45 µM) placed on LB agar and incubated at 37 °C overnight or for one hour. Bacteria were harvested in physiological water followed by serial dilutions and plated on LB agar containing two antibiotics: meropenem (MEM) 0.1 µg ml-1 and tetracycline (TET) 10 µg ml-1 to select transconjugants. Transconjugants were WGS and those devoid of mutations were selected for further analysis. Conjugation frequency was calculated as the ratio of transconjugants over donor after a 1h-mating followed by selection of transconjugants, donors, and recipients (MEM 0.1 µg ml-1 and TET 10 µg ml-1; MEM 0.1 µg ml-1, TET 10 µg ml-1 respectively).

Growth curves and relative fitness assessments

Overnight precultures in LB or M9 of the three ST38 E. coli plasmid-free isolates and their isogenic transconjugants were inoculated without and with antibiotic (MEM 0.1 µg ml-1) respectively. 96-well plates were inoculated with 100 µL of precultures diluted to 5*10⁵ colony forming unit (CFU) per milliliter in LB or M9 medium and incubated at 37 °C with shaking for 6 (LB) or 16 h (M9) in an automatic plate reader (Tecan infinite M Nano under i-control 2.0.10). Three independent experiments were carried out for each strain. For each independent biological replicate, the doubling time was obtained by calculating the mean from three to five technical replicates. The doubling time was calculated using the formula: G = ln (2)/µ_max. µ_max is the maximum growth rate, corresponding to the slope of the curve at exponential phase. The relative doubling time was calculated using the formula: W ₌ G_{transconjugant}/G_{plasmid-free.}

Experimental evolution

Overnight precultures in LB or in M9 of transconjugants and plasmid-free strains were diluted 1:200 into 10 ml fresh LB medium or into complemented M9 medium with and without MEM 0.1 µg ml-1 respectively and incubated at 37 °C with shaking (220 r.p.m, INFORS HT Minitron). For each experiment, five independent biological replicate cultures were evolved. Serial transfers were achieved for 28 days using 1:200 dilutions every day into 10 ml of fresh medium (ca. 8 generations per day). For the evolution of transconjugant lineages, MEM was added every day or every three days at 0.1 µg ml-1 as indicated. Each seven days, whole populations were collected and frozen at −80 °C, and bacterial pellets for DNA sequencing were obtained by centrifuging 1 ml of culture. Relative fitness of whole population through experimental evolution was monitored by growth curves. PCR using primers (Sigma-Aldrich) targeting bla_OXA-48 and repA (Supplementary Data 9) were performed on isolated colonies selected on MEM 0.1 µg.ml^-1 to identify potential integration of bla_OXA-48 (Tn6237) and pOXA-48 loss at day 7, day 14, day 21 and day 28.

Plasmid stability assay

Overnight precultures of pOXA-48 transconjugants in LB with MEM 0.1 µg ml-1 were set as day 0. Plasmid stability was assessed by serial passages using 1:200 dilutions into 10 ml fresh LB without MEM for ten days. Bacteria were collected and frozen at −80 °C at days 0, 5 and 10 and CFU were determined by plating 100 µl diluted cultures on LB agar with and without MEM 0.05 µg ml-1. Colonies growing on MEM were used as a proxy for plasmid-carrying bacteria as Tn6237 integration was estimated as a rare event under these conditions. Automatic bacterial colony counting was done with scan4000 (INTERSCIENCE) and plasmid stability was calculated as the ratio of plasmid-carrying bacteria over total population. Due to dilution and plating biaises this might occasionally result in ratios slightly superior to 100%. The stability of other plasmids was similarly tested after introduction by electroporation (pBbS8c, pBbE8c, pACYt, pUC19) or conjugation (pKPC, using CNR146C9 as donor) with adequate antibiotic selection (Supplementary Data 9).

Plasmid constructions for expression of wild-type and mutated copies of apsAB

F3141-F3140 operon, F3140 (apsB) and F3141 (apsA) were cloned under the control of a p_BAD inducible promotor in a pHV7 vector coding for apramycin resistance. All growing steps of pHV7-F3141-F3140 transformants were performed in the presence of glucose 0.2 % to repress apsAB induction. Two mutated versions of this operon where a stop codon was introduced in F3141 or F3140 sequence were generated during this cloning and the corresponding plasmids were used for activity testing. For site-directed mutagenesis, selected codons were modified in the pHV7-F3141-F3140 plasmid by using the Q5 Site-Directed Mutagenesis kit (New England Biolabs) according to manufacturer’s recommendations. pHV7 plasmid mutants were WGS (Eurofins). For chromosomal expression of apsAB (F3141-F3140) in MG1655, apsAB or apsAB homologs were cloned under the control of a p_BAD promoter in a miniTn7 transposon (miniTnapsAB) using primers described in Supplementary Data 9. miniTnapsAB was integrated downstream glmS, a neutral chromosomal position of E. coli K12 MG1655 strain following triparental mating as previously described⁴⁰. Antiplasmid activity of these constructs was analyzed by electroporating the ColE1 plasmid pBbE8c. An MG1655 strain harboring an empty miniTn7 (miniTnzeoR) integrated at glmS was used as control. When stated, chloramphenicol was added at 10 µg.ml^-1 to arrest bacterial growth 2.5 h after addition of arabinose 0.2%.

F3140-3141 (apsAB) chromosomal deletion and complementation

Complete deletion of the F3140-3141 (apsAB) operon was obtained by λ Red recombination⁴¹ by using an in-lab p15red vector carrying the lambda recombinase under the control of an inducible pBAD promotor (Supplementary Data 9). Primers (Sigma-Aldrich) used for PCR amplification of the ZeoR cassette between F3141-F3140 homology sequences are shown in Supplementary Data 9. zeoR marker was eliminated by using an in-lab p15Flip plasmid encoding the Flippase. SacB counterselection was used to eliminate recombineering plasmids. F3141-F3140 (apsAB) operon deletion were confirmed by PCR. The absence of mutation in the ST38_1∆apsAB clone used for experimental evolution was determined by WGS. For complementation, recipient cells were electroporated with the plasmids pHV7-empty, pHV7-F3141, pHV7-F3140, pHV7-F3141-F3140 and pHV7-F3141-F3140-mutants (Supplementary Data 9). pHV7-derivatives carrying strains were incubated 24 h in LB, apramycin 50 µg ml-1, 0.4% glucose to repress p_BAD activity (t₀); two serial passages using 1:200 dilutions were then performed in LB, apramycin 50 µg ml-1 and 0.02% arabinose to induce F3141-F3140 transcription. Cultures were diluted and plated on LB agar with or without MEM 0.05 µg ml-1 and pOXA-48 plasmids stability assessed as the ratio of resistant colonies over whole population.

Carbapenem susceptibility testing

MEM Minimal Inhibitory concentration (MIC) was determined by Etest (Biomerieux). The plates were inoculated by flooding 2 ml of bacterial culture (10⁶ bacteria ml-1), spread by a gentle rocking motion, and excessive liquid was removed leaving 0.5 ml of culture (+/− 10%). The flooding method was selected over swab streaking as it provides a more accurate reading of the Etest result.

Bacteriophages plaque assays

Phage plaque assays have been performed as previously described and by using the same phage collection²⁸ (bacteriophages lambda, T4, P1, 186cIts, CLB_P2, LF82_P8, AL505_P2, and T5, Supplementary Table 1). Phages were obtained as active cultures. The preys, E. coli K12 MG1655 strain and its isogenic derivatives chromosomally encoding the miniTnapsAB or the miniTnzeoR as control were grown overnight. Overnight cultures were diluted to 1:20 in the presence of arabinose 0.2% in LB medium and incubated at 37 °C for 3 h to induce apsAB expression. Bacterial lawns were prepared by mixing 200 μL of the induced culture with 100 μL of CaCl₂ 1 M and 20 ml of LB + 0.5% agar and poured onto 12 × 12 cm square plates of LB containing 0.2% arabinose. High-titer ( > 10⁸ pfu ml-1) stocks of phages lambda, T4, P1, 186cIts, CLB_P2, LF82_P8, AL505_P2, and T5 serially diluted were spotted on each plate and incubated at 37 °C overnight except for phage T7 incubated overnight at room temperature.

Whole genome sequencing and mutation identification

ST38_1, ST38_2 and ST38_3 strains were fully sequenced and used as reference sequences by combining Illumina sequencing and PacBio sequencing (ST38_1, ST38_2) or Oxford Nanopore technology (ST38_3). DNA was extracted at the exponential phase with Qiagen Puregene Yeast/Bact kit B. PacBio sequencing libraries were prepared with NANOBIND CBB KIT RT PacBio. ST38_3 Nanopore sequencing library was prepared by using Native Barcoding Kit 24 V14 (ref SQK-NBD114.24), and sequencing was performed with flowcell R10.4.1 (ref FLO-MIN114) on MinION Mk1C device with Guppy 5.0.14 as a basecaller. Long-read PacBio sequences were assembled with hybridSPAdes v 3.15.5 for hybrid assembly of short and long reads⁴². Hybrid assembly of short and long-read Nanopore sequences were performed by using Canu 2.2⁴³ and Circlator 1.5.5⁴⁴.

For Illumina sequencing, DNA was extracted from stationary phase cultures by using the Qiagen Blood and Tissue DNeasy kit, libraries were prepared with the NEBNext Ultra II FS DNA Library Prep Kit and sequencing was performed with NovaSeq6000 or NextSeq500 sequencing platforms. Short-read Illumina sequences were assembled with SPAdes 3.15.5⁴⁵ or aligned to the reference sequences by using Breseq 0.35.7⁴⁶ to identify SNPs, deletions, insertions and recombination events. IGV 2.11.9⁴⁷ was used to visually confirm mutation events. Illumina-reads of the whole population of the evolved lineages collected at day 7, day 14, day 21, and day 28 were analyzed for new junction evidence and coverage distribution by using Breseq 0.35.7⁴⁶ with the -p option for pool sequencing, a polymorphism frequency cutoff 2.5% option with at least 10 polymorphic reads. Plasmid coverage inside and outside the Tn6237 sequence was determined as the number of reads on three 15 kb-long regions, one in Tn6237, two outside, containing no IS, by using BAM files and the GRanges function of the GenomicAlignments package 1.22.1 in RStudio (R 3.6.3). The ratio was calculated as the ratio of the average of the number of reads mapping on the two regions outside Tn6237 to the number of reads mapping on Tn6237. IS1 new junctions detected by Breseq were tested as potential Tn6237 integration sites by PCR using a primer located near the new IS1 insertion site and a primer in bla_OXA-48 (Supplementary Data 9). PCR products were Sanger-sequenced to confirm the insertion site.

Phylogenetic analysis of E. coli ST38 and sequence annotation

To contextualize the three ST38 strains used in this work we performed a phylogenetic analysis using 1907 sequences retrieved from public databases (April 2020): 1248 assembled genome sequences from Enterobase [https://enterobase.warwick.ac.uk/], 149 assembled genomes from the NCBI and 510 sequences retrieved as reads and assembled with SPAdes 3.12.0⁴⁵ (Supplementary Data 10). QUAST 2.2⁴⁸ was used to assess the assembly quality and contigs shorter than 500 bp were filtered out for the phylogenetic analysis. A core genome alignment was generated with Parsnp 1.5.4⁴⁹, by using a finished genome sequence as reference. Maximum-Likelihood (ML) trees were generated with RAxML 8.2.12⁵⁰ using GTRGAMMA after removing regions of recombination with Gubbins⁵¹. The ST38 single locus variant ST963 strain CNRC6O47 was used as outgroup to root the phylogenetic tree. Trees were visualized and annotated using ITOL⁵² [https://itol.embl.de/].

Genomes annotation was performed by using Prokka 1.14.5⁵³. Resistome and plasmidome were characterized by using ABRicate 1.0.1 on the ResFinder db⁵⁴ (minimum coverage 60%, minimum identity 95%) and PlasmidFinder 2.1.1 (minimum coverage 60%, minimum identity 95%)⁵⁵. Identification of defense systems was performed by using DefenseFinder [https://defensefinder.mdmlab.fr/].

In silico characterization of ApsAB antiplasmid systems

Search for apsAB homologs were performed by using PSI-BLAST with three iterations on the recently introduced NCBI clustered nr database [https://blast.ncbi.nlm.nih.gov/]. This database is composed of representative sequences of clusters, that groups NCBI nr sequences sharing 90% identity and 90% length to other members of the cluster. Only sequences longer than 1100 a.a. residues (ApsA) and 500 a.a. residues (ApsB) were kept for PSI-BLAST iterations. CDS located downstream of apsA homologs were retrieved by using GCsnap 1.0.17⁵⁶. Protein sequences were aligned by using MuscleW 3.8.31 (default options) under Jalview 2.11.3⁵⁷ and a distance tree was created by Neighbor-Joining method using a BLOSUM62 matrix. Remote homology detection was also performed with the HHpred server [https://toolkit.tuebingen.mpg.de/tools/hhpred] using PDB_mmCIF70_18_jun, SCOPe70_2_08, CATH_S40_v4.3 and UniProt-SwissProt-viral70_3_nov_2021 as target databases accessed in September 2023. Structural modelling was performed with Alphafold2⁵⁸ implemented in Neurosnap [https://neurosnap.ai/] and the resulting models were used as templates for similarity searches with Foldseek [https://search.foldseek.com/search]. Structure annotation was performed with Pymol 2.5.5 (The PyMOL Molecular Graphics System, Version 2.5.5 Schrödinger, LLC.). Alpha-fold or icn3d pdb models of 20 ApsB-like or DdmE-like proteins (Supplementary Data 8) were recovered from Uniprot [https://www.uniprot.org/] or ncbi [https://www.ncbi.nlm.nih.gov/Structure/icn3d/] and aligned to ApsB structure model under Pymol 2.5.5 and the Root-mean-square deviation of atomic positions (RMSD) was used as a proxy to estimate protein structure similarity. The genomic environment of apsAB homologs from 12 E. coli strains belonging to different STs was compared by using the web version of Clinker [https://cagecat.bioinformatics.nl/tools/clinker].

Statistical tests

All graphs were generated with R (version 3.6.3) by using the R packages (tidyverse, forcats, ggplot2, ggpubr, rstatix, broom). Normality of data were assessed by using the Shapiro test. All the statistical analyzes were performed with a pairwise two sample t.test and p-value were corrected following Benjamini-Hochberg correction (FDR).

Reporting summary

Further information on research design is available in the Nature Portfolio Reporting Summary linked to this article.

Data availability

The data that support the findings of this study are provided within the manuscript and the associated supplementary materials. Full details and links to the publicly available databases [https://enterobase.warwick.ac.uk/], [https://defensefinder.mdmlab.fr/], [https://www.uniprot.org/], [https://www.ncbi.nlm.nih.gov/Structure/icn3d/], [https://www.ncbi.nlm.nih.gov/] used in bioinformatic analyzes are provided in the methods and their associated references. Complete genome sequences of ST38_1, ST38_2, and ST38_3 and of the three pOXA-48 plasmids have been deposited at DDBJ/EMBL/GenBank (BioProject PRJEB71895). The corresponding Genome_ID are provided in Supplementary Data 1 (for ST38 strains) and 9 (for pOXA-48s). Illumina reads for the three K. pneumoniae used as donors of pOXA-48 have also been deposited at DDBJ/EMBL/GenBank (BioProject PRJEB71895), and their accession numbers are given in Supplementary Data 1. Sequence data from individual colonies and pools are available from the corresponding author upon request. Source data are provided in this paper.

References

Murray, C. J. L. et al. Global burden of bacterial antimicrobial resistance in 2019: a systematic analysis. Lancet 399, 629–655 (2022).
Article CAS Google Scholar
Rozwandowicz, M. et al. Plasmids carrying antimicrobial resistance genes in Enterobacteriaceae. J. Antimicrob. Chemother. 73, 1121–1137 (2018).
Article CAS PubMed Google Scholar
Johnson, T. J. et al. In Vivo Transmission of an IncA/C Plasmid in Escherichia coli Depends on Tetracycline Concentration, and Acquisition of the Plasmid Results in a Variable Cost of Fitness. Appl. Environ. Microbiol. 81, 3561–3570 (2015).
Article ADS CAS PubMed PubMed Central Google Scholar
Starikova, I. et al. Fitness costs of various mobile genetic elements in Enterococcus faecium and Enterococcus faecalis. J. Antimicrob. Chemother. 68, 2755–2765 (2013).
Article CAS PubMed PubMed Central Google Scholar
Bergstrom, C. T., Lipsitch, M. & Levin, B. R. Natural selection, infectious transfer and the existence conditions for bacterial plasmids. Genetics 155, 1505–1519 (2000).
Article CAS PubMed PubMed Central Google Scholar
Brockhurst, M. A. & Harrison, E. Ecological and evolutionary solutions to the plasmid paradox. Trends Microbiol. 30, 534–543 (2022).
Article CAS PubMed Google Scholar
San Millan, A., Heilbron, K. & MacLean, R. C. Positive epistasis between co-infecting plasmids promotes plasmid survival in bacterial populations. ISME J. 8, 601–612 (2014).
Article CAS PubMed Google Scholar
Dahlberg, C. & Chao, L. Amelioration of the cost of conjugative plasmid carriage in Eschericha coli K12. Genetics 165, 1641–1649 (2003).
Article CAS PubMed PubMed Central Google Scholar
Harrison, E., Guymer, D., Spiers, A. J., Paterson, S. & Brockhurst, M. A. Parallel Compensatory Evolution Stabilizes Plasmids across the Parasitism-Mutualism Continuum. Curr. Biol. 25, 2034–2039 (2015).
Article CAS PubMed Google Scholar
Loftie-Eaton, W. et al. Compensatory mutations improve general permissiveness to antibiotic resistance plasmids. Nat. Ecol. Evol. 1, 1354–1363 (2017).
Article PubMed PubMed Central Google Scholar
Lopatkin, A. J. et al. Persistence and reversal of plasmid-mediated antibiotic resistance. Nat. Commun. 8, 1689 (2017).
Article ADS PubMed PubMed Central Google Scholar
Boyle, T. A. & Hatoum-Aslan, A. Recurring and emerging themes in prokaryotic innate immunity. Curr. Opin. Microbiol. 73, 102324 (2023).
Article CAS PubMed Google Scholar
Jaskólska, M., Adams, D. W. & Blokesch, M. Two defence systems eliminate plasmids from seventh pandemic Vibrio cholerae. Nature 604, 323–329 (2022).
Article ADS PubMed PubMed Central Google Scholar
Mayo-Muñoz, D., Pinilla-Redondo, R., Birkholz, N. & Fineran, P. C. A host of armor: Prokaryotic immune strategies against mobile genetic elements. Cell Rep. 42, 112672 (2023).
Article PubMed Google Scholar
Patiño-Navarrete, R. et al. Specificities and Commonalities of Carbapenemase-Producing Escherichia coli Isolated in France from 2012 to 2015. mSystems 7, e01169-21 (2022).
Article PubMed PubMed Central Google Scholar
Pitout, J. D. D., Peirano, G., Kock, M. M., Strydom, K.-A. & Matsumura, Y. The Global Ascendency of OXA-48-Type Carbapenemases. Clin. Microbiol. Rev. 33, e00102–e00119 (2019).
Article PubMed PubMed Central Google Scholar
Wielders, C. C. H. et al. Epidemiology of carbapenem-resistant and carbapenemase-producing Enterobacterales in the Netherlands 2017–2019. Antimicrob. Resist. Infect. Control 11, 57 (2022).
Article PubMed PubMed Central Google Scholar
Poirel, L., Bonnin, R. A. & Nordmann, P. Genetic features of the widespread plasmid coding for the Carbapenemase OXA-48. Antimicrob. Agents Chemother. 56, 559–562 (2012).
Article CAS PubMed PubMed Central Google Scholar
Alonso-del Valle, A. et al. Variability of plasmid fitness effects contributes to plasmid persistence in bacterial communities. Nat. Commun. 12, 2653 (2021).
Article ADS CAS PubMed PubMed Central Google Scholar
Hendrickx, A. P. A. et al. bla OXA-48-like genome architecture among carbapenemase-producing Escherichia coli and Klebsiella pneumoniae in the Netherlands. Microb. Genomics 7, 000512 (2021).
Turton, J. F. et al. Clonal expansion of Escherichia coli ST38 carrying a chromosomally integrated OXA-48 carbapenemase gene. J. Med. Microbiol. 65, 538–546 (2016).
Article PubMed Google Scholar
Fonseca, E. L., Morgado, S. M., Caldart, R. V. & Vicente, A. C. Global Genomic Epidemiology of Escherichia coli (ExPEC) ST38 Lineage Revealed a Virulome Associated with Human Infections. Microorganisms 10, 2482 (2022).
Article CAS PubMed PubMed Central Google Scholar
Emeraud, C. et al. Emergence and Polyclonal Dissemination of OXA-244–Producing Escherichia coli, France. Emerg. Infect. Dis. 27, 1206–1210 (2021).
Article CAS PubMed PubMed Central Google Scholar
Falgenhauer, L. et al. Cross-border emergence of clonal lineages of ST38 Escherichia coli producing the OXA-48-like carbapenemase OXA-244 in Germany and Switzerland. Int. J. Antimicrob. Agents 56, 106157 (2020).
Article CAS PubMed Google Scholar
Guenther, S. et al. Chromosomally encoded ESBL genes in Escherichia coli of ST38 from Mongolian wild birds. J. Antimicrob. Chemother. 72, 1310–1313 (2017).
Article CAS PubMed Google Scholar
Shawa, M. et al. Novel chromosomal insertions of ISEcp1-blaCTX-M-15 and diverse antimicrobial resistance genes in Zambian clinical isolates of Enterobacter cloacae and Escherichia coli. Antimicrob. Resist. Infect. Control 10, 79 (2021).
Article PubMed PubMed Central Google Scholar
León-Sampedro, R. et al. Pervasive transmission of a carbapenem resistance plasmid in the gut microbiota of hospitalized patients. Nat. Microbiol. 6, 606–616 (2021).
Article PubMed PubMed Central Google Scholar
Rousset, F. et al. Phages and their satellites encode hotspots of antiviral systems. Cell Host Microbe 30, 740–753.e5 (2022).
Article CAS PubMed PubMed Central Google Scholar
White, M. F. Structure, function and evolution of the XPD family of iron-sulfur-containing 5’–>3’ DNA helicases. Biochem. Soc. Trans. 37, 547–551 (2009).
Article CAS PubMed Google Scholar
Knizewski, L., Kinch, L. N., Grishin, N. V., Rychlewski, L. & Ginalski, K. Realm of PD-(D/E)XK nuclease superfamily revisited: detection of novel families with modified transitive meta profile searches. BMC Struct. Biol. 7, 40 (2007).
Article PubMed PubMed Central Google Scholar
Lisitskaya, L., Aravin, A. A. & Kulbachinskiy, A. DNA interference and beyond: structure and functions of prokaryotic Argonaute proteins. Nat. Commun. 9, 5165 (2018).
Article ADS PubMed PubMed Central Google Scholar
Miyoshi, T., Ito, K., Murakami, R. & Uchiumi, T. Structural basis for the recognition of guide RNA and target DNA heteroduplex by Argonaute. Nat. Commun. 7, 11846 (2016).
Article ADS CAS PubMed PubMed Central Google Scholar
Kuang, H., Yang, Y., Luo, H. & Lv, X. The impact of three carbapenems at a single-day dose on intestinal colonization resistance against carbapenem-resistant Klebsiella pneumoniae. mSphere 8, e00479–23 (2023).
Article PubMed PubMed Central Google Scholar
Dionisio, F., Zilhão, R. & Gama, J. A. Interactions between plasmids and other mobile genetic elements affect their transmission and persistence. Plasmid 102, 29–36 (2019).
Article CAS PubMed Google Scholar
San Millan, A., Toll-Riera, M., Qi, Q. & MacLean, R. C. Interactions between horizontally acquired genes create a fitness cost in Pseudomonas aeruginosa. Nat. Commun. 6, 6845 (2015).
Article ADS CAS PubMed Google Scholar
Georjon, H. & Bernheim, A. The highly diverse antiphage defence systems of bacteria. Nat. Rev. Microbiol. 21, 686–700 (2023).
Article CAS PubMed Google Scholar
Millman, A. et al. An expanded arsenal of immune systems that protect bacteria from phages. Cell Host Microbe 30, 1556–1569.e5 (2022).
Article CAS PubMed Google Scholar
Song, X. et al. Catalytically inactive long prokaryotic Argonaute systems employ distinct effectors to confer immunity via abortive infection. Nat. Commun. 14, 6970 (2023).
Article ADS CAS PubMed PubMed Central Google Scholar
Jackson, S. A., Fellows, B. J. & Fineran, P. C. Complete Genome Sequences of the Escherichia coli Donor Strains ST18 and MFDpir. Microbiol. Resour. Announc. 9, e01014–e01020 (2020).
Article CAS PubMed PubMed Central Google Scholar
Choi, K.-H. et al. A Tn7-based broad-range bacterial cloning and expression system. Nat. Methods 2, 443–448 (2005).
Article CAS PubMed Google Scholar
Datsenko, K. A. & Wanner, B. L. One-step inactivation of chromosomal genes in Escherichia coli K-12 using PCR products. Proc. Natl Acad. Sci. 97, 6640–6645 (2000).
Article ADS CAS PubMed PubMed Central Google Scholar
Antipov, D., Korobeynikov, A., McLean, J. S. & Pevzner, P. A. hybrid SPA des: an algorithm for hybrid assembly of short and long reads. Bioinformatics 32, 1009–1015 (2016).
Article CAS PubMed Google Scholar
Koren, S. et al. Canu: scalable and accurate long-read assembly via adaptive k -mer weighting and repeat separation. Genome Res. 27, 722–736 (2017).
Article CAS PubMed PubMed Central Google Scholar
Hunt, M. et al. Circlator: automated circularization of genome assemblies using long sequencing reads. Genome Biol. 16, 294 (2015).
Article PubMed PubMed Central Google Scholar
Prjibelski, A., Antipov, D., Meleshko, D., Lapidus, A. & Korobeynikov, A. Using SPAdes De Novo Assembler. Curr. Protoc. Bioinforma. 70, e102 (2020).
Deatherage, D. E. & Barrick, J. E. Identification of Mutations in Laboratory-Evolved Microbes from Next-Generation Sequencing Data Using breseq. in Engineering and Analyzing Multicellular Systems (eds. Sun, L. & Shou, W.) vol. 1151 165–188 (Springer New York, New York, NY, 2014).
Robinson, J. T. et al. Integrative genomics viewer. Nat. Biotechnol. 29, 24–26 (2011).
Article CAS PubMed PubMed Central Google Scholar
Gurevich, A., Saveliev, V., Vyahhi, N. & Tesler, G. QUAST: quality assessment tool for genome assemblies. Bioinformatics 29, 1072–1075 (2013).
Article CAS PubMed PubMed Central Google Scholar
McKenna, A. et al. The Genome Analysis Toolkit: A MapReduce framework for analyzing next-generation DNA sequencing data. Genome Res. 20, 1297–1303 (2010).
Article CAS PubMed PubMed Central Google Scholar
Stamatakis, A. RAxML version 8: a tool for phylogenetic analysis and post-analysis of large phylogenies. Bioinformatics 30, 1312–1313 (2014).
Article CAS PubMed PubMed Central Google Scholar
Croucher, N. J. et al. Rapid phylogenetic analysis of large samples of recombinant bacterial whole genome sequences using Gubbins. Nucleic Acids Res. 43, e15–e15 (2015).
Article PubMed Google Scholar
Letunic, I. & Bork, P. Interactive Tree Of Life (iTOL) v4: recent updates and new developments. Nucleic Acids Res. 47, W256–W259 (2019).
Article CAS PubMed PubMed Central Google Scholar
Seemann, T. Prokka: rapid prokaryotic genome annotation. Bioinformatics 30, 2068–2069 (2014).
Article CAS PubMed Google Scholar
Zankari, E. et al. Identification of acquired antimicrobial resistance genes. J. Antimicrob. Chemother. 67, 2640–2644 (2012).
Article CAS PubMed PubMed Central Google Scholar
Carattoli, A. et al. In silico detection and typing of plasmids using plasmidfinder and plasmid multilocus sequence typing. Antimicrob. Agents Chemother. 58, 3895–3903 (2014).
Article PubMed PubMed Central Google Scholar
Pereira, J. GCsnap: interactive snapshots for the comparison of protein-coding genomic contexts. J. Mol. Biol. 433, 166943 (2021).
Article CAS PubMed Google Scholar
Waterhouse, A. M., Procter, J. B., Martin, D. M. A., Clamp, M. & Barton, G. J. Jalview Version 2—a multiple sequence alignment editor and analysis workbench. Bioinformatics 25, 1189–1191 (2009).
Article CAS PubMed PubMed Central Google Scholar
Jumper, J. et al. Highly accurate protein structure prediction with AlphaFold. Nature 596, 583–589 (2021).
Article ADS CAS PubMed PubMed Central Google Scholar
Sullivan, M. J., Petty, N. K. & Beatson, S. A. Easyfig: a genome comparison visualizer. Bioinformatics 27, 1009–1010 (2011).
Article CAS PubMed PubMed Central Google Scholar

Download references

Acknowledgements

We thank Alexandre Almeida for his critical reading of the manuscript, David Vallenet for his help in the structure analysis of ApsAB, and Jean-Marc Ghigo for providing E. coli strain MFDpir. Sequencing was performed at the Biomics Platform, C2RT, Institut Pasteur, Paris, France, supported by France Génomique (ANR-10-INBS-09) and IBISA. This work was supported by “Investissement d’Avenir” program, LABEX IBEID (Grant ANR-10-LABX-62-IBEID) and SEQ2DIAG (ANR-20-PAMR-0010). PDZ received a fellowship from CIOSPB (Center national de l’Information, de l’Orientation Scolaire et Professionnelle, et des Bourses) of the Ministry of Higher Education, Scientific Research and Innovation of Burkina Faso and from the LABEX IBEID.

Author information

These authors contributed equally: Nicolas Cabanel, Guilhem Royer.
These authors jointly supervised this work: Philippe Glaser, Isabelle Rosinski-Chupin.

Authors and Affiliations

Ecology and Evolution of Antibiotic Resistance Unit, Institut Pasteur, Paris, France
Pengdbamba Dieudonné Zongo, Nicolas Cabanel, Guilhem Royer, Philippe Glaser & Isabelle Rosinski-Chupin
Sorbonne Université, Paris, France
Pengdbamba Dieudonné Zongo
Université Paris Cité, Paris, France
Pengdbamba Dieudonné Zongo, Nicolas Cabanel, Guilhem Royer, Florence Depardieu, Philippe Glaser & Isabelle Rosinski-Chupin
Synthetic Biology Unit, Institut Pasteur, Paris, France
Florence Depardieu
UMR AgroEcologie 1347, INRAe, Université Bourgogne Franche-Comté, Dijon, France
Alain Hartmann
Team ReSIST, INSERM UMR 1184, Paris-Saclay University, Le Kremlin-Bicêtre, France
Thierry Naas
Department of Bacteriology-Hygiene, Bicêtre Hospital, APHP, Le Kremlin-Bicêtre, France
Thierry Naas
Associated French National Reference Center for Antibiotic Resistance, Bicêtre Hospital, Le Kremlin-Bicêtre, France
Thierry Naas

Authors

Pengdbamba Dieudonné Zongo
View author publications
You can also search for this author in PubMed Google Scholar
Nicolas Cabanel
View author publications
You can also search for this author in PubMed Google Scholar
Guilhem Royer
View author publications
You can also search for this author in PubMed Google Scholar
Florence Depardieu
View author publications
You can also search for this author in PubMed Google Scholar
Alain Hartmann
View author publications
You can also search for this author in PubMed Google Scholar
Thierry Naas
View author publications
You can also search for this author in PubMed Google Scholar
Philippe Glaser
View author publications
You can also search for this author in PubMed Google Scholar
Isabelle Rosinski-Chupin
View author publications
You can also search for this author in PubMed Google Scholar

Contributions

I.R.C., P.G. and P.D.Z. designed the study. P.D.Z., I.R.C., NC, F.D. and G.R. performed the experiments. P.D.Z., P.G. and I.R.C. analyzed the data. F.D., T.N. and A.H. provided materials. I.R.C., P.G. and P.D.Z. wrote the manuscript with input from G.R. All authors read and approved the final manuscript.

Corresponding author

Correspondence to Isabelle Rosinski-Chupin.

Ethics declarations

Competing interests

All authors declare no competing interests.

Peer review

Peer review information

Nature Communications thanks the anonymous reviewers for their contribution to the peer review of this work. A peer review file is available.

Additional information

Publisher’s note Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Supplementary information

Supplementary Information

Peer Review File

Description of Additional Supplementary Files

Supplementary Data 1

Supplementary Data 2

Supplementary Data 3

Supplementary Data 4

Supplementary Data 5

Supplementary Data 6

Supplementary Data 7

Supplementary Data 8

Supplementary Data 9

Supplementary Data 10

Reporting Summary

Source data

Source Data

Rights and permissions

Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article’s Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article’s Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecommons.org/licenses/by/4.0/.

Reprints and permissions

About this article

Cite this article

Zongo, P.D., Cabanel, N., Royer, G. et al. An antiplasmid system drives antibiotic resistance gene integration in carbapenemase-producing Escherichia coli lineages. Nat Commun 15, 4093 (2024). https://doi.org/10.1038/s41467-024-48219-y

Download citation

Received: 08 February 2024
Accepted: 24 April 2024
Published: 15 May 2024
DOI: https://doi.org/10.1038/s41467-024-48219-y

Comments

By submitting a comment you agree to abide by our Terms and Community Guidelines. If you find something abusive or that does not comply with our terms or guidelines please flag it as inappropriate.