# Transcriptional Landscape of a bla KPC-2 Plasmid and Response to Imipenem Exposure in Escherichia coli TOP10

Agnès Jousset, Isabelle Rosinski-Chupin, Julie Takissian, Philippe Glaser, Rémy Bonnin, Thierry Naas

► To cite this version:

Agnès Jousset, Isabelle Rosinski-Chupin, Julie Takissian, Philippe Glaser, Rémy Bonnin, et al. Transcriptional Landscape of a bla KPC-2 Plasmid and Response to Imipenem Exposure in Escherichia coli TOP10. Frontiers in Microbiology, 2018, 9, pp.2929. 10.3389/fmicb.2018.02929. pasteur-02565586

HAL Id: pasteur-02565586

https://pasteur.hal.science/pasteur-02565586

Submitted on 6 May 2020

HAL is a multi-disciplinary open access
archive for the deposit and dissemination of sci-
entific research documents, whether they are pub-
lished or not. The documents may come from
teaching and research institutions in France or
abroad, or from public or private research centers.

L'archive ouverte pluridisciplinaire HAL, est
destinée au dépôt et à la diffusion de documents
scientifiques de niveau recherche, publiés ou non,
émanant des établissements d'enseignement et de
recherche français ou étrangers, des laboratoires
publics ou privés.

<figure><img src="image_1.png" /></figure>

Distributed under a Creative Commons Attribution 4.0 International License
<figure><img src="image_2.png" /></figure>

Check for
updates

# Transcriptional Landscape of a *bla*<sub>KPC-2</sub> Plasmid and Response to Imipenem Exposure in *Escherichia coli* TOP10

Agnès B. Jousset<sup>1,2,3,4</sup>, Isabelle Rosinski-Chupin<sup>4,5</sup>, Julie Takissian<sup>3,4</sup>, Philippe Glaser<sup>4,5</sup>, Rémy A. Bonnin<sup>2,3,4*</sup> and Thierry Naas<sup>1,2,3,4*</sup>

## OPEN ACCESS

**Edited by:**
Raffaele Zarrilli,
Department of Public Health,
Università degli Studi di Napoli
Federico II, Italy

**Reviewed by:**
Vincenzo Di Pilato,
Università degli Studi di Firenze, Italy
Eliana De Gregorio,
Università degli Studi di Napoli
Federico II, Italy

*Correspondence:
Rémy A. Bonnin
remy.bonnin@u-psud.fr
Thierry Naas
thierry.naas@aphp.fr

**Specialty section:**
This article was submitted to *Antimicrobials, Resistance and Chemotherapy*, a section of the journal *Frontiers in Microbiology*

Received: 19 September 2018

Accepted: 14 November 2018

Published: 03 December 2018

**Citation:**
Jousset AB, Rosinski-Chupin I, Takissian J, Glaser P, Bonnin RA and Naas T (2018) Transcriptional Landscape of a *bla*<sub>KPC-2</sub> Plasmid and Response to Imipenem Exposure in *Escherichia coli* TOP10. *Front. Microbiol.* 9:2929. doi: 10.3389/fmicb.2018.02929

The diffusion of KPC-2 carbapenemase is closely related to the spread of *Klebsiella pneumoniae* of the clonal-group 258 and linked to IncFII<sub>K</sub> plasmids. Little is known about the biology of multi-drug resistant plasmids and the reasons of their successful dissemination. Using *E. coli* TOP10 strain harboring a multi-replicon IncFII<sub>K</sub>-IncFIB *bla*<sub>KPC-2</sub>-gene carrying plasmid pBIC1a from *K. pneumoniae* ST-258 clinical isolate BIC-1, we aimed to identify basal gene expression and the effects of imipenem exposure using whole transcriptome approach by RNA sequencing (RNA-Seq). Independently of the antibiotic pressure, most of the plasmid-backbone genes were expressed at low levels. The most expressed pBIC1a genes were involved in antibiotic resistance (*bla*<sub>KPC-2</sub>, *bla*<sub>TEM</sub> and *aph*(3')-I), in plasmid replication and conjugation, or associated to mobile elements. After antibiotic exposure, 34% of *E. coli* (pBIC1a) genome was differentially expressed. Induction of oxidative stress response was evidenced, with numerous upregulated genes of the *SoxRS/OxyR* oxydative stress regulons, the Fur regulon (for iron uptake machinery), and *IscR* regulon (for iron sulfur cluster synthesis). Nine genes carried by pBIC1a were up-regulated, including the murein DD-endopeptidase *mepM* and the copper resistance operon. Despite the presence of a carbapenemase, we observed a major impact on *E. coli* (pBIC1a) whole transcriptome after imipenem exposure, but no effect on the level of transcription of antimicrobial resistance genes. We describe adaptive responses of *E. coli* to imipenem-induced stress, and identified plasmid-encoded genes that could be involved in resistance to stressful environments.

**Keywords:** carbapenemase, KPC-producing plasmids, transcriptome, RNA-seq, oxidative-stress

## INTRODUCTION

*Klebsiella pneumoniae* is a prominent opportunistic pathogen for hospital- and community-acquired infections (Navon-Venezia et al., 2017). The increasing incidence of KPC-producing *K. pneumoniae* (KPC-Kp) in health care facilities is a cause of global concern, especially in countries where the carbapenemase KPC is endemic, e.g., the United States, Israel, Greece, and

<footnote n="1"><sup>1</sup> Department of Bacteriology-Parasitology-Hygiene, Bicêtre Hospital, Assistance Publique-Hôpitaux de Paris, Le Kremlin-Bicêtre, France, <sup>2</sup> Associated French National Reference Center for Antibiotic Resistance, Le Kremlin-Bicêtre, France, <sup>3</sup> EA7361 “Structure, dynamic, function and expression of broad spectrum β-lactamases”, Faculty of Medicine, Paris-Sud University, Le Kremlin-Bicêtre, France, <sup>4</sup> Joint Research Unit Evolution and Ecology of Resistance to Antibiotics, Institut Pasteur-APHP-University Paris Sud, Paris, France, <sup>5</sup> CNRS, UMRS 3525, Paris, France</footnote>
Italy (Bonomo et al., 2017). The *bla*<sub>KPC</sub> gene, coding for an Ambler class A carbapenemase is mostly associated with *K. pneumoniae* but has been reported to a lesser extent in other Enterobacteriaceae, in *Pseudomonas* spp. and in *Acinetobacter baumannii* (Cuzon et al., 2011, 2013). KPC confers resistance or decreased susceptibility to almost all β-lactams, and KPC-producing isolates are often resistant to many other non-β-lactam drugs because of the co-occurrence of *bla*<sub>KPC</sub> gene with resistance genes to other classes of antibiotics. This multidrug resistance leaves only limited therapeutic options for antimicrobial treatment, and thereby results in high mortality rates (Lee and Burgess, 2012).

The *bla*<sub>KPC</sub> genes have been identified on several transferable plasmids of different incompatibility groups (e.g., IncFII<sub>K</sub>, IncA/C, IncF, IncN, IncP, IncR, IncX, and ColE1), but IncFII<sub>K</sub>-type plasmids are the most common and are involved in their spreading (Chen et al., 2014b). The pKpQIL plasmid was the first characterized IncFII<sub>K2</sub> *bla*<sub>KPC</sub>-carrying plasmid in Israel in 2006, and is strongly associated to the *K. pneumoniae* ST258 pandemic clone (Chen et al., 2013). Subsequently, pKpQIL-like plasmids seemed to have played a major role in KPC dissemination with several reports worldwide (Baraniak et al., 2011; García-Fernández et al., 2012; Hidalgo-Grass et al., 2012; Chen et al., 2013, 2014a; Papagiannitsis et al., 2016; Doumith et al., 2017). Complete sequencing of KPC-associated plasmids revealed their plasticity with multiple rearrangements that occurred between pKpQIL-like plasmids and pKPN3, a non-KPC-encoding IncFII<sub>K1</sub>-FIB plasmid described in *K. pneumoniae* MGH78578 (also known as ATCC 700721; Chen et al., 2013; Jousset et al., 2018). However, little is known about their basic biology, and due to the presence of multiple antimicrobial resistance genes on IncFII<sub>K</sub> plasmids, their broad host range, and their rapid dissemination, there is a need to understand the mechanisms underlying regulation and maintenance of these plasmids.

KPC is mostly associated to *K. pneumoniae* of a single clonal group (CG) designated CG258 containing 43 different MLST single-locus variants with ST258, ST11, ST512, and ST340 being the more predominant ones (Deleo et al., 2014; Rojas et al., 2017, 2018). KPC has also been identified in other Enterobacteriaceae such as *Escherichia coli*, *Citrobacter* sp, *Enterobacter* sp, *Serratia marcescens*, *Proteus*. *Mirabilis*, and *Morganella morganii*. KPC-producing *E. coli* isolates usually have low MICs for carbapenems and can remain susceptible according to clinical breakpoints (Landman et al., 2010). Regardless the level of carbapenem resistance, options for treatment of infections due to KPC-producing Enterobacteriaceae are limited and current clinical evidence for treatment guidelines are still lacking. Even if novel β-lactamase inhibitors (such as avibactam) are increasingly used with success (van Duin et al., 2017), carbapenem-based regimens remain a therapeutic option in many countries. Imipenem and meropenem remain useful mainly in combination with other classes of antibiotics or even another carbapenem, such as ertapenem (Bonomo et al., 2017). Pharmacokinetic parameters of the used molecule and Minimal Inhibitory Concentration value of the KPC-producing strain are also crucial in patient outcome (Tumbarello et al., 2012; Daikos et al., 2014). Tumbarello et al.

showed that patients infected by carbapenemase-producing Enterobacteriaceae with imipenem MIC values ≥4 µg/mL had worse outcomes than patients whose isolates had MIC ≤ 2 µg/mL. On the contrary, therapeutic failure with use of carbapenems against KPC-2-producing isolates possessing low MICs were also reported (Weisenberg et al., 2009; Daikos and Markogiannakis, 2011). Previous studies were made to understand the transcriptional regulation of *bla*<sub>KPC</sub> revealing the versatile expression of this gene (Naas et al., 2012; Cheruvanky et al., 2017; Girlich et al., 2017). Measuring *bla*<sub>KPC</sub> RNA expression revealed that *bla*<sub>KPC</sub> gene copy number and the presence of different promoters cannot be strictly correlated with MICs determination and are not sufficient to explain the level of resistance regarding carbapenems (Roth et al., 2011; Naas et al., 2012). These data triggered our interest in the regulation of *bla*<sub>KPC</sub> gene and in the entire plasmid transcriptome, in particular under antibiotic exposure.

Effect of β-lactams on *E. coli* has been characterized using transcriptomic approaches. In addition to interacting with their targets, bactericidal antibiotics induce parallel changes in bacterial metabolism that promote the formation of reactive oxygen species, which play a role in cell death (Kohanski et al., 2007; Dwyer et al., 2014). β-lactams were also shown by Miller and colleagues to trigger the SOS response *via* the activation of the DpiAB two-component system (Miller, 2004). However, the dynamic underlying oxidative and antibiotic-induced SOS stress response activation needs further exploration. In most of these studies, sub-inhibitory doses of ampicillin were used on fully susceptible *E. coli* isolates expressing no β-lactamase.

To give new insights into the biology of multidrug resistance successful plasmids, the transcriptome of *E. coli* harboring pBIC1a, an IncFII<sub>K2</sub>-FIB-type *bla*<sub>KPC-2</sub> carrying plasmid was performed. To gain more information about the complex regulation of the *bla*<sub>KPC</sub> gene and all the transcripts of the cell, the transcriptome of *E. coli* harboring pBIC1a was performed with or without imipenem exposure. Moreover, to mimic physiopathological conditions that occur during a carbapenem-based treatment, we chose to analyse the response of *E. coli* (pBIC1a) to high doses of imipenem (10 times the MIC value of the strain). We described adaptive responses of *E. coli* (pBIC1a) to imipenem-induced stress with disruption of several metabolic pathways.

# RESULTS AND DISCUSSION

## Mosaic Structure of pBIC1a Plasmid

Plasmid pBIC1a is a *bla*<sub>KPC-2</sub>-encoding plasmid of 170,415 bp in-size from a ST258 *K. pneumoniae* BIC-1 (Jousset et al., 2018). This isolate was recovered from a contaminated endoscope in France in 2009 (Naas et al., 2010). This plasmid exhibited an average GC content of 53% and 193 predicted open reading frames (ORF). The overall structure of pBIC1a was highly similar (83% query coverage and 99% nucleotide identity by Blast) to the non-KPC carrying plasmid pKPN3 (Figure 1). Plasmids pBIC1a and pKPN3 shares a common region of 120-kb. This region contains genes responsible for plasmid replication, maintenance,
<figure><img src="image_3.png" /></figure>

transmission, and heavy metal (arsenic, copper, and silver) resistance, and presents the replicase *repA2* of IncFIB type.

The rest of pBIC1a plasmid shared homology with pKpQIL-like plasmids such as pKpQIL (27% query coverage and 99.9% nucleotide identity) and pKpQIL-IT (36% query coverage and 99.3% nucleotide identity; **Figure 1**). This region contains the end of the *tra* operon, another replicase *repA* of IncFII<sub>K2</sub> type according to replicon sequence typing (Carattoli et al., 2014), making pBIC1a a multi-replicon IncFII<sub>K2</sub>-IncFIB plasmid. The multi-replicon status can expand plasmid's host range replication (Villa et al., 2010). Like most of the pKpQIL-like plasmids, the variant *bla*<sub>KPC-2</sub> gene was found as part of the class II transposon *Tn*4401a. Additional antimicrobial resistance genes were present: *bla*<sub>TEM-1</sub>, *bla*<sub>OXA-9</sub> (disrupted by a frameshift mutation), *aph(3')-I* and partial *aadA1*. Furthermore, when compared to original pKpQIL, pBIC1a had an additional truncated resistance gene carried on a composite transposon-like element, IS26-Δ*aphA1* -Δ*IS26*, located downstream of *Tn*4401a. This element has been described in pKpQIL-IT, but differed by an inversion of a 3,872 nt fragment that occurred downstream (**Figure 1**).

Therefore, pBIC1a seems to be the result of several transpositions and recombination events between pKpQIL, pKpQIL-IT, and pKPN3 plasmids (**Figure 2**). Differently rearranged pKPN3 plasmids have been described in KPC-Kp from the USA, Greece and Italy (Chen et al., 2013; Wright et al., 2014; Papagiannitsis et al., 2016; Doumith et al., 2017). pIT-O6C07 (LT009688.1), pBK32179 (JX430448) and pGMI16-005_01 (NZ_CP028181.1) are three fully sequenced IncFII<sub>K2</sub>-IncFIB plasmids carrying *bla*<sub>KPC</sub>, resulting from recombination between pKPN3 and pKpQIL-like plasmids (**Figure 2**). pIT-O6C07 is the closest to pBIC1a with 97% query coverage and 99% nucleotide identity (Papagiannitsis et al., 2016). The two plasmids only differed by a 4,695 nt region containing *bla*<sub>TEM-1</sub> and *bla*<sub>OXA-9</sub>, that is missing in pIT-O6C07. Alignment of the four KPC-producing plasmids deriving from pKPN3 indicated that recombination did not occur at the same spot precisely (**Figure 2**). Recombination events between pKpQIL-like and pKPN3-like plasmids seem common, likely due to the presence of highly homologous regions and their frequent coexistence in *K. pneumoniae*.
<figure><img src="image_4.png" /></figure>

**FIGURE 2** | Nucleotidic alignment of *bla*~KPC~-carrying plasmids deriving from pKPN3. pKPN3 (NC_009649.1) is a plasmid from *K. pneumoniae* MGH78578 that carries no β-lactamase. pIT-O6C07 (LT009688.1), pBK32179 (JX430448), pGMI16-005_01 (NZ_CP028181.1) are KPC-producing plasmids deriving from recombination between pKPN3 and pKpQIL-like plasmids. From the inside to the outside, the different circles represent: pBIC1a GC content (black) and open reading frames annotation (red), pBIC1a GC skew (purple and green), alignment with pIT-O6C07 (pink), alignment with pBK32179 (red), alignment with pGMI16-005 (blue), alignment with pKPN3 (green). This circular representation was performed using BRIG.

## Transcriptome Analysis of *E. coli* (pBIC1a)

*E. coli* TOP10, transformed with the KPC-carrying plasmid pBIC1a, was grown at mid-log phase in rich media (Brain Heart Infusion (BHI) liquid media). Two conditions in biological triplicates were studied: liquid exponential growth as control and liquid exponential growth supplemented with imipenem at final concentration of 5 µg/ml (10 times the MIC value of TOP10-pBIC1a for imipenem, which was determined at 0.5 µg/ml using E-Test). This concentration was chosen because it can be achieved in human serum during the treatment of a bacteriemia due to a KPC-producing Enterobacteriaceae after intravenous administration of 1,000 mg of imipenem-cilastatin (Singlas, 1990; Signs et al., 1992). Imipenem exposure was limited to 10 min to evaluate the immediate cellular response and limit bacterial cell death. The viability of bacterial cells was verified by time-kill analysis and no effect on survival was observed after 10 min exposure to imipenem (Figure S1).

After 10 min of exposure, extraction of total RNA was performed on each replicate. Libraries were obtained from ribo-depleted RNA. After sequencing, a mean of 27.9 ± 3.3 million reads per individual library was obtained. Reads were aligned on *E. coli* TOP10 (CP000948.1) and pBIC1a (CP022574.1) DNA sequences, with 97.25% ± 0.78% of the reads mapping on non-ribosomal regions (Table S1). The principal component analysis (PCA) and the clustering heatmap plots revealed that the samples clustered to their biological replicates (Figure S2).

## Transcriptional Landscape of pBIC1a in the Absence of Imipenem

In BHI medium and exponential growth phase, the replicases *repA* IncFII~K~, *repA* IncFIB, and the regulator *repB* IncFII~K~ genes were expressed with reads per kilobase per million mapped reads (RPKM) values of 127, 224, and 643 respectively. These genes were used for subsequent comparisons as they are involved in
plasmid replication and should be expressed consistently during exponential growth due to continuous cell division (Lang et al., 2012).

Alignment of reads on pBIC1a sequence was performed and heatmap visualization was used to identify highly or weakly transcribed regions using IGV software (Robinson et al., 2011).

Most of the plasmid genes were expressed at low levels or not expressed during exponential phase in BHI medium, with only 34 of predicted CDS (17%) having RKPM values above the *repA* IncFII<sub>K</sub> gene RPKM value (Table 1; Figure 3). These results are comparable with those found in a previous transcriptomic analysis of a multidrug resistance-encoding plasmid where most of the backbone was transcriptionally silent in *E. coli* (Lang et al., 2012). Among the plasmid scaffold genes, *higA* antitoxin gene, part of HigA/HigB type II toxin/antitoxin system was highly expressed suggesting a role in pBIC1a maintenance in *E. coli*. Most genes involved in plasmid transfer function were expressed at low levels, except the genes encoding TraM (relaxosome protein), TraA (pilin precursor), TraL (pilus assembly protein), TraT (surface exclusion protein), and TraS (entry exclusion protein; Figure 3).

Transformation of pBIC1a into the *E. coli* background only resulted into a slight increase in the MIC for imipenem (0.5 µg/ml compared to 0.25 µg/ml in *E. coli* TOP10 wild-type strain). However, our RNA-seq experiments showed that *bla*<sub>KPC-2</sub> gene was the most expressed plasmidic gene (Table 1; Figure 3). Two promoters were previously shown to contribute to *bla*<sub>KPC-2</sub> gene expression (Naas et al., 2012). Although our RNA-sequencing protocol was not designed to precisely map transcription start sites, the high coverage by sequencing reads upstream the p1 promoter sequence strongly suggests promoter p2 as the main promoter for *bla*<sub>KPC-2</sub> gene transcription (Figure S3). In addition, RNA-seq analysis revealed no clear transcription termination with reads continuously covering more than one kb downstream *bla*<sub>KPC-2</sub> stop codon (Figure S3). In agreement with this observation, no rho-independent termination site could be predicted by bioinformatic tools (http://rna.igmors.u-psud.fr/toolbox/arnold/). This transcription covered in antisense the transposase gene of ISKpn6. It might therefore act as a repressor for ISKpn6 transposition as an antisense RNA, and/or play a role in the regulation of KPC through mRNA stability as observed for other 3'-UTR (Ruiz de los Mozos et al., 2013).

Other resistance genes carried by the plasmid, namely *bla*<sub>TEM-1</sub> and *aph(3′)-Ia* genes were also highly expressed (Table 1; Figure 3). The high expression of *aph(3′)-Ia* gene is in accordance with the high level resistance to kanamycin observed (MIC > 32 µg/ml). The association of *bla*<sub>KPC</sub> gene with other antibiotic resistance determinants also highly expressed provides a very simple scenario for a carbapenemase to spread as a hitchhiker gene, especially in the absence of carbapenem selection.

Different operons involved in metal resistance were identified on pBIC1a: the *cus* (*cusCFBA*) and *cop* (*copABCDR*) operons involved in copper and other cations efflux, and the *ars* operon (*arsRDABC*) in arsenical resistance. Here, only the regulators genes *cusR* and *arsR* were highly transcribed. CusR is part of

the two component system CusRS that activates the transcription of *cusCFBA* genes in the presence of Cu/Ag (Xiao et al., 2017). pKPN3-like plasmids carrying the highly similar sil or cus operons are known to confer resistance to high silver concentration in *K. pneumoniae* and in *Enterobacter cloacae* (Finley et al., 2015). On the contrary, ArsR acts as a repressor for the arsenical resistance operon (Ren et al., 2017).

Several mobile elements and insertion sequences (IS) were identified on pBIC1a including one copy of ISKpn6, ISKpn7, ISKpn26-like, ISKpn28, ISKpn32, IS1, IS903, and ISEcl1 (Figure 1). In addition to these elements, several fragments of IS26 were identified. These IS were not or weakly expressed that is in accordance with the tight regulation of ISs (Nagy and Chandler, 2004). An interesting feature of ISKpn31 was the expression of its passenger gene. Whereas, the transposase gene was almost not expressed, its passenger gene (pBIC_0005) belonged to the top expressed genes of the plasmid. This encoded protein is named p019 and has been identified on matrix-assisted laser desorption ionization-time of flight mass spectrometry (MALDI-TOF MS) spectra of *K. pneumoniae* harboring pKpQIL-like plasmids (Lau et al., 2014; Partridge, 2014). The presence of a peptide signal indicates a probable export but this protein does not share homology with protein with known function.

Recently, Buckner et al. studied the expression of two *bla*<sub>KPC</sub>-producing pKpQIL-like plasmids in *K. pneumoniae*, and revealed that 60% of plasmid-encoded genes were expressed. These genes included genes involved in replication, transmission, stability, recombination and toxin-antitoxin systems (Buckner et al., 2018). Comparison of pKpQIL-UK, pKpQIL-D2, and pBIC1a transcripts revealed common expressed genes, including *bla*<sub>KPC</sub>, *bla*<sub>TEM</sub>, *repA*, *traALMTS*, *finO*, *tnp* gene coding for Tn3 transposase. There are growing evidence for the existence of “core” expressed genes associated to pKpQIL-like plasmids, expressed in both *K. pneumoniae* and *E. coli* genetic backgrounds.

## *Impenem-Induced Response in E. coli (pBIC1a): Global Changes*

After imipenem exposure, we observed a major impact on *E. coli* (pBIC1a) whole transcriptome with 1,563 RNAs being differentially expressed (with adjusted *p*-values ≤0.05) among a total of 4,550 (34%). 35% (1,535 over 4,357) of chromosomally-encoded genes and 14% (28 over 193) of pBIC1a genes were affected by imipenem exposure. 765 RNAs were up-regulated and 798 were down-regulated (Table S2). Among differentially expressed genes, 154 RNAs had a Fold-Change (FC) >2 and 329 RNAs had a FC <0.5. Differential expression of eight selected genes (four up-regulated and four down-regulated genes) was confirmed by RT-qPCR (Figure S4).

## *Impenem-Induced Response at the Plasmid pBIC1a Level*

Effect of imipenem on all pBIC1a transcripts is indicated in Table S3. After imipenem addition in culture medium, nine predicted CDS of pBIC1a were up-regulated but none with a FC >2 (Table 2). Transcription of copper resistance operon (*copABDRS/pcoE*) was induced after imipenem exposure, with
<table><caption>TABLE 1 | Expression analysis of pBIC1a genes without imipenem exposure using RPKM values.</caption><thead><tr><th>pBIC1a</th><th>Gene name</th><th>Orientation</th><th>Annotation</th><th>RPKM value</th></tr></thead><tbody><tr><td>pBIC_00001</td><td></td><td>+</td><td>Aldehyde dehydrogenase</td><td>35</td></tr><tr><td>pBIC_00002</td><td></td><td>+</td><td>Putative resolvase</td><td>12</td></tr><tr><td>pBIC_00003</td><td></td><td>+</td><td>Hypothetical protein</td><td>23</td></tr><tr><td>pBIC_00004</td><td></td><td>+</td><td>Hypothetical protein</td><td>144</td></tr><tr><td>pBIC_00005</td><td></td><td>-</td><td>p019 of ISKpn31</td><td>1,840</td></tr><tr><td>pBIC_00006</td><td></td><td>-</td><td>Transposase</td><td>20</td></tr><tr><td>pBIC_00007</td><td></td><td>+</td><td>Mobile element protein</td><td>8</td></tr><tr><td>pBIC_00008</td><td></td><td>+</td><td>Tnp-ISKpn6</td><td>8</td></tr><tr><td>pBIC_00009</td><td>bla<sub>KPC-2</sub></td><td>-</td><td>KPC-2</td><td>2,706</td></tr><tr><td>pBIC_00010</td><td></td><td>-</td><td>Hypothetical protein</td><td>294</td></tr><tr><td>pBIC_00011</td><td>istA</td><td>-</td><td>ISKpn7</td><td>0</td></tr><tr><td>pBIC_00012</td><td>istB</td><td>-</td><td>ISKpn7</td><td>1</td></tr><tr><td>pBIC_00013</td><td>TnpA</td><td>-</td><td>Transposase</td><td>4</td></tr><tr><td>pBIC_00014</td><td>TnpR</td><td>+</td><td>Resolvase</td><td>13</td></tr><tr><td>pBIC_00015</td><td></td><td>+</td><td>Hypothetical protein</td><td>2,261</td></tr><tr><td>pBIC_00016</td><td></td><td>+</td><td>Hypothetical protein</td><td>76</td></tr><tr><td>pBIC_00017</td><td></td><td>+</td><td>Transposase</td><td>65</td></tr><tr><td>pBIC_00018</td><td></td><td>-</td><td>Hypothetical protein</td><td>220</td></tr><tr><td>pBIC_00019</td><td>finO</td><td>-</td><td>IncF plasmid conjugative transfer fertility inhibition protein FinO</td><td>328</td></tr><tr><td>pBIC_00020</td><td>repA</td><td>-</td><td>DNA replication protein (IncFIIk)</td><td>127</td></tr><tr><td>pBIC_00021</td><td>repB</td><td>-</td><td>Replication regulatory protein</td><td>643</td></tr><tr><td>pBIC_00022</td><td>nucC</td><td>-</td><td>Endonuclease</td><td>13</td></tr><tr><td>pBIC_00023</td><td></td><td>-</td><td>Hypothetical protein</td><td>35</td></tr><tr><td>pBIC_00024</td><td>yihA</td><td>-</td><td>Hypothetical protein</td><td>13</td></tr><tr><td>pBIC_00025</td><td></td><td>-</td><td>Hypothetical protein</td><td>14</td></tr><tr><td>pBIC_00026</td><td></td><td>-</td><td>Hypothetical protein</td><td>11</td></tr><tr><td>pBIC_00027</td><td>finO</td><td>-</td><td>IncF plasmid conjugative transfer fertility inhibition protein FinO</td><td>72</td></tr><tr><td>pBIC_00028</td><td>traX</td><td>-</td><td>IncF plasmid conjugative transfer pilin acetylase TraX</td><td>14</td></tr><tr><td>pBIC_00029</td><td>tral</td><td>-</td><td>IncF plasmid conjugative transfer DNA-nicking and unwinding protein Tral</td><td>10</td></tr><tr><td>pBIC_00030</td><td>traD</td><td>-</td><td>IncF plasmid conjugative transfer protein TraD</td><td>7</td></tr><tr><td>pBIC_00031</td><td></td><td>-</td><td>Hypothetical protein</td><td>8</td></tr><tr><td>pBIC_00032</td><td>traT</td><td>-</td><td>IncF plasmid conjugative transfer surface exclusion protein TraT</td><td>282</td></tr><tr><td>pBIC_00033</td><td>traS</td><td>-</td><td>IncF plasmid conjugative transfer surface exclusion protein TraS</td><td>231</td></tr><tr><td>pBIC_00034</td><td>traG</td><td>-</td><td>IncF plasmid conjugative transfer protein TraG</td><td>56</td></tr><tr><td>pBIC_00035</td><td>traH</td><td>-</td><td>IncF plasmid conjugative transfer pilus assembly protein TraH</td><td>6</td></tr><tr><td>pBIC_00036</td><td></td><td>-</td><td>Hypothetical protein</td><td>4</td></tr><tr><td>pBIC_00037</td><td>trbB</td><td>-</td><td>IncF plasmid conjugative transfer protein TrbB</td><td>7</td></tr><tr><td>pBIC_00038</td><td>traQ</td><td>-</td><td>IncF plasmid conjugative transfer protein TraQ</td><td>6</td></tr><tr><td>pBIC_00039</td><td>traF</td><td>-</td><td>IncF plasmid conjugative transfer pilus assembly protein TraF</td><td>6</td></tr><tr><td>pBIC_00040</td><td></td><td>-</td><td>Hypothetical protein</td><td>19</td></tr><tr><td>pBIC_00041</td><td></td><td>-</td><td>Hypothetical protein</td><td>14</td></tr><tr><td>pBIC_00042</td><td>trbE</td><td>-</td><td>IncF plasmid conjugative transfer protein TrbE</td><td>5</td></tr><tr><td>pBIC_00043</td><td>traN</td><td>-</td><td>IncF plasmid conjugative transfer protein TraN</td><td>10</td></tr><tr><td>pBIC_00044</td><td>trbC</td><td>-</td><td>IncF plasmid conjugative transfer protein TrbC</td><td>7</td></tr><tr><td>pBIC_00045</td><td>traU</td><td>-</td><td>IncF plasmid conjugative transfer pilus assembly protein TraU</td><td>10</td></tr><tr><td>pBIC_00046</td><td>traW</td><td>-</td><td>IncF plasmid conjugative transfer pilus assembly protein TraW</td><td>8</td></tr><tr><td>pBIC_00047</td><td>trbl</td><td>-</td><td>IncF plasmid conjugative transfer protein Trbl</td><td>9</td></tr><tr><td>pBIC_00048</td><td>traC</td><td>-</td><td>IncF plasmid conjugative transfer pilus assembly protein TraC</td><td>11</td></tr><tr><td>pBIC_00049</td><td></td><td>-</td><td>Hypothetical protein</td><td>12</td></tr><tr><td>pBIC_00050</td><td></td><td>-</td><td>Hypothetical protein</td><td>11</td></tr><tr><td>pBIC_00051</td><td></td><td>-</td><td>Hypothetical protein</td><td>9</td></tr></tbody></table>

(Continued)
<table><caption>TABLE 1 | Continued</caption><thead><tr><th>pBIC1a</th><th>Gene name</th><th>Orientation</th><th>Annotation</th><th>RPKM value</th></tr></thead><tbody><tr><td>pBIC_00052</td><td></td><td>-</td><td>Hypothetical protein</td><td>7</td></tr><tr><td>pBIC_00053</td><td></td><td>-</td><td>Hypothetical protein</td><td>14</td></tr><tr><td>pBIC_00054</td><td>traV</td><td>-</td><td>IncF plasmid conjugative transfer pilus assembly protein TraV</td><td>22</td></tr><tr><td>pBIC_00055</td><td>traP</td><td>-</td><td>IncF plasmid conjugative transfer protein TraP</td><td>22</td></tr><tr><td>pBIC_00056</td><td>traB</td><td>-</td><td>IncF plasmid conjugative transfer pilus assembly protein TraB</td><td>21</td></tr><tr><td>pBIC_00057</td><td>traK</td><td>-</td><td>IncF plasmid conjugative transfer pilus assembly protein TraK</td><td>25</td></tr><tr><td>pBIC_00058</td><td>traE</td><td>-</td><td>IncF plasmid conjugative transfer pilus assembly protein TraE</td><td>36</td></tr><tr><td>pBIC_00059</td><td>traL</td><td>-</td><td>IncF plasmid conjugative transfer pilus assembly protein TraL</td><td>433</td></tr><tr><td>pBIC_00060</td><td>traA</td><td>-</td><td>IncF plasmid conjugative transfer pilin protein TraA</td><td>656</td></tr><tr><td>pBIC_00061</td><td>traY</td><td>-</td><td>IncF plasmid conjugative transfer regulator TraY</td><td>65</td></tr><tr><td>pBIC_00062</td><td>traJ</td><td>-</td><td>IncF plasmid conjugative transfer regulator TraJ</td><td>34</td></tr><tr><td>pBIC_00063</td><td>traM</td><td>-</td><td>IncF plasmid conjugative transfer mating signal transduction protein TraM</td><td>429</td></tr><tr><td>pBIC_00064</td><td></td><td>+</td><td>Lytic transglycosylase</td><td>295</td></tr><tr><td>pBIC_00065</td><td></td><td>-</td><td>Hypothetical protein</td><td>16</td></tr><tr><td>pBIC_00066</td><td></td><td>-</td><td>Unnamed protein product</td><td>20</td></tr><tr><td>pBIC_00067</td><td></td><td>+</td><td>Hypothetical protein</td><td>5</td></tr><tr><td>pBIC_00068</td><td>higA</td><td>-</td><td>Antitoxine HigA</td><td>223</td></tr><tr><td>pBIC_00069</td><td></td><td>-</td><td>Hypothetical protein phage-related</td><td>400</td></tr><tr><td>pBIC_00070</td><td></td><td>-</td><td>Hypothetical protein</td><td>17</td></tr><tr><td>pBIC_00071</td><td></td><td>-</td><td>Hypothetical protein</td><td>14</td></tr><tr><td>pBIC_00072</td><td></td><td>-</td><td>Hypothetical protein</td><td>23</td></tr><tr><td>pBIC_00073</td><td></td><td>-</td><td>Hypothetical protein</td><td>163</td></tr><tr><td>pBIC_00074</td><td></td><td>-</td><td>Hypothetical protein</td><td>144</td></tr><tr><td>pBIC_00075</td><td></td><td>-</td><td>Hypothetical protein</td><td>15</td></tr><tr><td>pBIC_00076</td><td></td><td>+</td><td>Mobile element protein</td><td>22</td></tr><tr><td>pBIC_00077</td><td></td><td>-</td><td>Hypothetical protein</td><td>19</td></tr><tr><td>pBIC_00078</td><td></td><td>-</td><td>Hypothetical protein</td><td>2</td></tr><tr><td>pBIC_00079</td><td>psiA</td><td>-</td><td>PsiA protein</td><td>2</td></tr><tr><td>pBIC_00080</td><td>psiB</td><td>-</td><td>PsiB protein</td><td>3</td></tr><tr><td>pBIC_00081</td><td></td><td>-</td><td>Hypothetical protein</td><td>4</td></tr><tr><td>pBIC_00082</td><td>parB</td><td>-</td><td>ParB</td><td>5</td></tr><tr><td>pBIC_00083</td><td>ssb</td><td>-</td><td>Single-stranded DNA-binding protein</td><td>7</td></tr><tr><td>pBIC_00084</td><td></td><td>+</td><td>Hypothetical protein</td><td>5</td></tr><tr><td>pBIC_00085</td><td></td><td>-</td><td>SAM-dependent methyltransferase</td><td>3</td></tr><tr><td>pBIC_00086</td><td></td><td>-</td><td>Hypothetical protein</td><td>3</td></tr><tr><td>pBIC_00087</td><td>ydaB</td><td>-</td><td>YdaB</td><td>2</td></tr><tr><td>pBIC_00088</td><td></td><td>-</td><td>Hypothetical protein</td><td>11</td></tr><tr><td>pBIC_00089</td><td></td><td>-</td><td>Hypothetical protein</td><td>53</td></tr><tr><td>pBIC_00090</td><td></td><td>-</td><td>DNA polymerase III theta subunit</td><td>52</td></tr><tr><td>pBIC_00091</td><td></td><td>-</td><td>Antirestriction protein klcA</td><td>6</td></tr><tr><td>pBIC_00092</td><td></td><td>+</td><td>Hypothetical protein</td><td>9</td></tr><tr><td>pBIC_00093</td><td></td><td>-</td><td>Hypothetical protein</td><td>47</td></tr><tr><td>pBIC_00094</td><td></td><td>-</td><td>Retron-type RNA-directed DNA polymerase; Ontology_term</td><td>49</td></tr><tr><td>pBIC_00095</td><td>umuD</td><td>+</td><td>Error-prone repair protein UmuD</td><td>17</td></tr><tr><td>pBIC_00096</td><td>umuC</td><td>+</td><td>Error-prone lesion bypass DNA polymerase V (UmuC)</td><td>11</td></tr><tr><td>pBIC_00097</td><td></td><td>-</td><td>Mobile element protein</td><td>15</td></tr><tr><td>pBIC_00098</td><td></td><td>-</td><td>Mobile element protein</td><td>26</td></tr><tr><td>pBIC_00099</td><td></td><td>+</td><td>Hypothetical protein</td><td>2</td></tr><tr><td>pBIC_00100</td><td></td><td>-</td><td>Hypothetical protein</td><td>13</td></tr><tr><td>pBIC_00101</td><td>sopB/parB</td><td>-</td><td>Chromosome (plasmid) partitioning protein ParB</td><td>55</td></tr><tr><td>pBIC_00102</td><td>sopA/parA</td><td>-</td><td>Chromosome (plasmid) partitioning protein ParA</td><td>53</td></tr></tbody></table>

(Continued)
<table><caption>TABLE 1 | Continued</caption><thead><tr><th>pBIC1a</th><th>Gene name</th><th>Orientation</th><th>Annotation</th><th>RPKM value</th></tr></thead><tbody><tr><td>pBIC_00103</td><td>repA2</td><td>+</td><td>Replicase (IncFI)</td><td>224</td></tr><tr><td>pBIC_00104</td><td>-</td><td>-</td><td>Resolvase/Recombinase</td><td>36</td></tr><tr><td>pBIC_00105</td><td>-</td><td>+</td><td>Hypothetical protein</td><td>3</td></tr><tr><td>pBIC_00106</td><td>-</td><td>+</td><td>Sensor domain-containing diguanylate cyclase</td><td>43</td></tr><tr><td>pBIC_00107</td><td>-</td><td>+</td><td>Mobile element protein (IS903B)</td><td>2</td></tr><tr><td>pBIC_00108</td><td>-</td><td>+</td><td>Hypothetical protein</td><td>36</td></tr><tr><td>pBIC_00109</td><td>-</td><td>+</td><td>Putative membrane protein YjcC</td><td>12</td></tr><tr><td>pBIC_00110</td><td>-</td><td>-</td><td>Hypothetical protein</td><td>4</td></tr><tr><td>pBIC_00111</td><td>-</td><td>-</td><td>Hypothetical protein</td><td>16</td></tr><tr><td>pBIC_00112</td><td>-</td><td>-</td><td>Secondary glycine betaine transporter BetU</td><td>54</td></tr><tr><td>pBIC_00113</td><td>-</td><td>-</td><td>Mobile element protein (IS1F)</td><td>10</td></tr><tr><td>pBIC_00114</td><td>-</td><td>-</td><td>Mobile element protein (IS1F)</td><td>18</td></tr><tr><td>pBIC_00115</td><td>-</td><td>-</td><td>Mobile element protein (ISEc11)</td><td>53</td></tr><tr><td>pBIC_00116</td><td>-</td><td>-</td><td>Mobile element protein (ISEc11)</td><td>23</td></tr><tr><td>pBIC_00117</td><td>-</td><td>+</td><td>Mobile element protein</td><td>21</td></tr><tr><td>pBIC_00118</td><td>-</td><td>+</td><td>Mobile element protein</td><td>4</td></tr><tr><td>pBIC_00119</td><td>-</td><td>-</td><td>Silver-binding protein</td><td>107</td></tr><tr><td>pBIC_00120</td><td>cusS</td><td>-</td><td>Osmosensitive K+ channel histidine kinase</td><td>48</td></tr><tr><td>pBIC_00121</td><td>cusR</td><td>-</td><td>Copper-sensing two-component system response regulator CusR</td><td>139</td></tr><tr><td>pBIC_00122</td><td>cusC</td><td>+</td><td>Cation efflux system protein CusC precursor</td><td>5</td></tr><tr><td>pBIC_00123</td><td>cusF</td><td>+</td><td>Cation efflux system protein CusF precursor</td><td>17</td></tr><tr><td>pBIC_00124</td><td>cusB</td><td>+</td><td>Cobalt/zinc/cadmium efflux RND transporter membrane fusion protein CzcB family</td><td>10</td></tr><tr><td>pBIC_00125</td><td>cusA</td><td>+</td><td>Cobalt-zinc-cadmium resistance protein CzcAB Cation efflux system protein CusA</td><td>7</td></tr><tr><td>pBIC_00126</td><td>copG</td><td>+</td><td>CopG protein</td><td>27</td></tr><tr><td>pBIC_00127</td><td>-</td><td>+</td><td>Lead cadmium zinc and mercury transporting ATPase B Copper-translocating P-type ATPase</td><td>14</td></tr><tr><td>pBIC_00128</td><td>-</td><td>+</td><td>hypothetical protein</td><td>42</td></tr><tr><td>pBIC_00129</td><td>mepM</td><td>-</td><td>Cell wall endopeptidase-family M23/M37</td><td>49</td></tr><tr><td>pBIC_00130</td><td>-</td><td>-</td><td>Copper-binding protein PcoE</td><td>22</td></tr><tr><td>pBIC_00131</td><td>-</td><td>+</td><td>Multicopper oxidase</td><td>51</td></tr><tr><td>pBIC_00132</td><td>copB</td><td>+</td><td>Copper resistance protein B</td><td>32</td></tr><tr><td>pBIC_00133</td><td>copC</td><td>+</td><td>Copper resistance protein CopC</td><td>53</td></tr><tr><td>pBIC_00134</td><td>copD</td><td>+</td><td>Copper resistance protein D</td><td>46</td></tr><tr><td>pBIC_00135</td><td>-</td><td>+</td><td>DNA-binding heavy metal response regulator</td><td>57</td></tr><tr><td>pBIC_00136</td><td>-</td><td>+</td><td>Heavy metal sensor histidine kinase</td><td>33</td></tr><tr><td>pBIC_00137</td><td>-</td><td>+</td><td>Probable copper-binding protein</td><td>30</td></tr><tr><td>pBIC_00138</td><td>-</td><td>-</td><td>Hypothetical protein</td><td>25</td></tr><tr><td>pBIC_00139</td><td>-</td><td>-</td><td>Mobile element protein</td><td>18</td></tr><tr><td>pBIC_00140</td><td>-</td><td>+</td><td>Aquaporin Z</td><td>13</td></tr><tr><td>pBIC_00141</td><td>-</td><td>-</td><td>N-acetyltransferase</td><td>48</td></tr><tr><td>pBIC_00142</td><td>arsR</td><td>+</td><td>Arsenical resistance operon repressor</td><td>21</td></tr><tr><td>pBIC_00143</td><td>arsH</td><td>+</td><td>Arsenic resistance protein ArsH</td><td>39</td></tr><tr><td>pBIC_00144</td><td>-</td><td>+</td><td>Hypothetical protein</td><td>47</td></tr><tr><td>pBIC_00145</td><td>arsD</td><td>+</td><td>Arsenical resistance operon trans-acting repressor ArsD</td><td>9</td></tr><tr><td>pBIC_00146</td><td>-</td><td>+</td><td>Arsenical pump-driving ATPase</td><td>4</td></tr><tr><td>pBIC_00147</td><td>-</td><td>+</td><td>Arsenical pump-driving ATPase</td><td>4</td></tr><tr><td>pBIC_00148</td><td>-</td><td>+</td><td>Hypothetical protein</td><td>12</td></tr><tr><td>pBIC_00149</td><td>-</td><td>-</td><td>Arsenate reductase</td><td>91</td></tr><tr><td>pBIC_00150</td><td>-</td><td>-</td><td>Arsenic efflux pump protein</td><td>18</td></tr></tbody></table>

(Continued)
<table><caption>TABLE 1 | Continued</caption><thead><tr><th>pBIC1a</th><th>Gene name</th><th>Orientation</th><th>Annotation</th><th>RPKM value</th></tr></thead><tbody><tr><td>pBIC_00151</td><td></td><td>-</td><td>Arsenical pump-driving ATPAs</td><td>48</td></tr><tr><td>pBIC_00152</td><td>arsD</td><td>-</td><td>Arsenical resistance operon trans-acting repressor ArsD</td><td>77</td></tr><tr><td>pBIC_00153</td><td>arsR</td><td>-</td><td>Arsenical resistance operon repressor</td><td>163</td></tr><tr><td>pBIC_00154</td><td></td><td>+</td><td>Hypothetical protein</td><td>77</td></tr><tr><td>pBIC_00155</td><td>ydeA</td><td>+</td><td>YdeA protein</td><td>85</td></tr><tr><td>pBIC_00156</td><td></td><td>+</td><td>DNA binding protein</td><td>47</td></tr><tr><td>pBIC_00157</td><td></td><td>+</td><td>Hypothetical protein</td><td>37</td></tr><tr><td>pBIC_00158</td><td></td><td>+</td><td>Haemolysin expression modulating protein</td><td>87</td></tr><tr><td>pBIC_00159</td><td></td><td>+</td><td>Hypothetical protein</td><td>191</td></tr><tr><td>pBIC_00160</td><td></td><td>+</td><td>Hypothetical protein</td><td>535</td></tr><tr><td>pBIC_00161</td><td></td><td>+</td><td>Hypothetical protein</td><td>269</td></tr><tr><td>pBIC_00162</td><td></td><td>+</td><td>N-acetyltransferase</td><td>333</td></tr><tr><td>pBIC_00163</td><td></td><td>+</td><td>Transposase</td><td>28</td></tr><tr><td>pBIC_00164</td><td></td><td>-</td><td>Hypothetical protein</td><td>33</td></tr><tr><td>pBIC_00165</td><td></td><td>-</td><td>Mobile element protein</td><td>79</td></tr><tr><td>pBIC_00166</td><td></td><td>-</td><td>Mobile element protein</td><td>28</td></tr><tr><td>pBIC_00167</td><td></td><td>-</td><td>Zn-dependent protease</td><td>57</td></tr><tr><td>pBIC_00168</td><td></td><td>-</td><td>Diguanylate cyclase</td><td>170</td></tr><tr><td>pBIC_00169</td><td></td><td>-</td><td>Hypothetical protein</td><td>210</td></tr><tr><td>pBIC_00170</td><td>clp</td><td>-</td><td>ATP-dependent Clp protease ATP-binding subunit ClpA</td><td>409</td></tr><tr><td>pBIC_00171</td><td>hsp</td><td>-</td><td>Small HspC2 heat shock protein</td><td>795</td></tr><tr><td>pBIC_00172</td><td></td><td>-</td><td>DNA binding protein</td><td>591</td></tr><tr><td>pBIC_00173</td><td></td><td>-</td><td>Transposase</td><td>14</td></tr><tr><td>pBIC_00174</td><td></td><td>-</td><td>Hypothetical protein</td><td>11</td></tr><tr><td>pBIC_00175</td><td></td><td>+</td><td>Mobile element protein (ISKpn26)</td><td>2</td></tr><tr><td>pBIC_00176</td><td></td><td>+</td><td>Mobile element protein (ISKpn26)</td><td>15</td></tr><tr><td>pBIC_00177</td><td></td><td>+</td><td>Hypothetical protein</td><td>7</td></tr><tr><td>pBIC_00178</td><td></td><td>-</td><td>Mobile element protein</td><td>1</td></tr><tr><td>pBIC_00179</td><td>cbbR</td><td>-</td><td>RuBisCO operon transcriptional regulator CbbR</td><td>2</td></tr><tr><td>pBIC_00180</td><td></td><td>-</td><td>Phosphonate dehydrogenase (NAD-dependent phosphite dehydrogenase)</td><td>4</td></tr><tr><td>pBIC_00181</td><td>phnE</td><td>-</td><td>Phosphonate ABC transporter permease protein phnE</td><td>5</td></tr><tr><td>pBIC_00182</td><td></td><td>-</td><td>Phosphonate ABC transporter phosphate-binding periplasmic component</td><td>2</td></tr><tr><td>pBIC_00183</td><td>phnC1</td><td>-</td><td>Phosphonate ABC transporter ATP-binding protein</td><td>3</td></tr><tr><td>pBIC_00184</td><td></td><td>+</td><td>Mobile element protein</td><td>4</td></tr><tr><td>pBIC_00185</td><td></td><td>-</td><td>Mobile element protein (IS26)</td><td>6</td></tr><tr><td>pBIC_00186</td><td></td><td>+</td><td>TnpA transposase (IS26)</td><td>8</td></tr><tr><td>pBIC_00187</td><td></td><td>+</td><td>Mobile element protein (partial IS26)</td><td>22</td></tr><tr><td>pBIC_00188</td><td>ΔblaOXA-9</td><td>+</td><td>Beta-lactamase class D (partial OXA-9)</td><td>6</td></tr><tr><td>pBIC_00189</td><td>ΔblaOXA-9</td><td>+</td><td>Beta-lactamase class D (partial OXA-9)</td><td>6</td></tr><tr><td>pBIC_00190</td><td></td><td>+</td><td>Mobile element protein (partial Tn3)</td><td>241</td></tr><tr><td>pBIC_00191</td><td>blaTEM-1</td><td>+</td><td>Beta-lactamase TEM-1</td><td>175</td></tr><tr><td>pBIC_00192</td><td></td><td>-</td><td>Mobile element protein (partial IS26)</td><td>13</td></tr><tr><td>pBIC_00193</td><td>aph(3')-I</td><td>+</td><td>Aminoglycoside phosphotransferase</td><td>655</td></tr></tbody></table>

Genes with RPKM values above the RPKM value of IncFII<sub>k</sub> replicase repA are highlighted in gray.

FC values <2. With two oxidation states Cu⁺ and Cu²⁺, copper is a cofactor in redox enzymes that uses dioxygen as a substrate (Rademacher and Masepohl, 2012). Bacteria synthesize various cuproenzymes, which play important roles in cellular processes such as energy transduction, iron mobilization and oxidative stress (Rademacher and Masepohl, 2012). In excess, copper

becomes toxic as it interacts with free proteinogenic thiol groups, destabilizes iron-sulfur cofactors, and possibly leads to formation of ROS (Chillappagari et al., 2010; Rademacher and Masepohl, 2012). Activation of the copper resistance operon could act as a defense factor to limit toxicity induced by high copper intracellular concentration, hence ROS production (Bondarczuk
<figure><img src="image_5.png" /><figcaption>**FIGURE 3** | Expression of pBIC1a genes in the absence of antibiotic. The coverage (expressed as the average density of reads over 500 nt fragments) was computed on each strand of pBIC1a sequence and visualized by using a heatmap representation (IGV software). Open Reading Frames (ORFs) annotated on the strand (+ or sens) are indicated in red boxes and ORFs annotated on the strand (− or rev) are indicated in blue boxes. Known annotation of the most expressed genes are mentioned on a supplementary line for both strands. Both replicases (repA IncFII<sub>K</sub>-type and repA2 IncFIB-type) are highlighted in green. Genes associated to plasmid transfer are boxed in red.</figcaption></figure>

and Piotrowska-Seget, 2013). The chromosomally-encoded small RybA (also known as MntS) could be involved in this positive regulation. It was up-regulated after imipenem exposure (FC = 3.5 *p*-value = 4.14e-39). RybA transcripts are stabilized under peroxide stress and members of the CusR regulon are downregulated in a Δ*rybA* mutant, suggesting that RybA could be a positive regulator of genes involved in copper detoxification (Gerstle et al., 2012).

Among the membrane associated regulated genes, a plasmid-encoded murein DD-endopeptidase called *mepM* (pBIC_00129) was up-regulated (FC = 1.37 *p*-value = 1.9e-04). Search for conserved domains using the NCBI Database identified MepM as a member of LytM domain containing proteins, that include numerous autolysins involved in peptidoglycan remodeling (Uehara et al., 2009). Bacteria mutated for such enzymes form longer chains due to a default in cell separation (Visweswaran et al., 2013). On the contrary, overexpression of LysM-associated proteins lead to increased autolysis (Uehara et al., 2009; Visweswaran et al., 2013). This small activation of *mepM* could be an indirect consequence of the interaction of imipenem with its Penicillin Binding Proteins (PBP), reflecting an early activation of bacterial lysis (Uehara et al., 2009), though not observed through time-kill experiments (Figure S1).

The chromosomal copy of *mepM* (also known as *yebA*) of *E. coli* MG1655 was previously shown to be up-regulated under oxidative stress (Seo et al., 2015). Here, we only observed the up-regulation of the plasmid-encoded endopeptidase. The two proteins differed by their size and molecular weight, with plasmid-encoded MepM being a 27.08 kDa protein vs. 49.04 kDa for the chromosomal copy. They shared only 41% of amino-acid identity, suggesting that both enzymes might have different cellular functions.

We specifically looked at the level of *bla*<sub>KPC-2</sub> gene expression that was not modified in the presence of imipenem (FC = 1.12, *p*-value = 0.14). There is no genetic argument for an inducible expression of *bla*<sub>KPC-2</sub> and our results were in agreement with a constitutive expression (Naas et al., 2012).

Finally, 19 predicted CDS of pBIC1a were down-regulated with *p*-values <0.05, and four of them with FC values <0.5 (Table 2). Six down-regulated genes coded for transposases from various IS. Down-regulation of two helix-turn-helix domain proteins (pBIC_00098 and pBIC_00172) was also observed with *p*-values of 5.31E-07 and 1.51E-08, respectively. Such proteins may act as global regulators for gene expression (Aravind et al., 2005).

## Imipenem-Induced Response at the Chromosomal Level

Most of the transcriptomic variations induced by imipenem exposure affected *E. coli* chromosomal genes, with 756 genes (17.4%) being up-regulated and 779 genes (17.9%) down-regulated over 4,357 genes.

To focused on the most significant changes, as a first step, we analyzed global changes on genes with FC values >2 (154 genes) or <0.5 (325 genes). To identify common metabolic pathways, genes were analyzed for gene ontology (GO) term enrichment using PANTHER web application. Enriched GO terms from up-regulated genes (adjusted *p*-values ranging from 2.23e-4 to 0.053) and from down-regulated genes (adjusted *p*-values ranging from 2.22e-10 to 0.047) were further explored separately on the REVIGO web application to visualize relationships and redundancy among GO terms (Supek et al., 2011; **Figure 4**). An activation of the aerobic respiration was observed with enrichment in the pathway of electron transport chain (*p*-value = 0.014). Expression of genes involved in cluster Fe-S assembly (*p*-value = 0.016) and in ion transport (*p*-value = 0.015) was also affected. Besides, downregulated genes were mostly associated to anaerobic respiration (*p*-value 2.22e-10) and generation of precursors of metabolites and energy (*p*-value 2.55e-10). Up and down-regulated genes were enriched in oxydo-reduction process (*p*-values 0.015 and 1.82e-07 respectively), suggesting that oxidative stress might have occurred in bacterial cells after imipenem exposure (Figure 4).
<table><caption>**TABLE 2** | Differentially expressed genes carried by pBIC1a.</caption><thead><tr><th>Id</th><th>Name</th><th>Counts (without imipenem)</th><th>Counts (with imipenem)</th><th>Fold-Change</th><th>log2(Fold-Change)</th><th>Adjusted p-value</th><th>Annotation</th></tr></thead><tbody><tr><td colspan="8"><strong>DOWN-REGULATED GENES</strong></td></tr><tr><td>pBIC_00006</td><td>pBIC_00006</td><td>110</td><td>77</td><td>0.697</td><td>-0.522</td><td>0.018900674</td><td>Transposase (ISKpn31)</td></tr><tr><td>pBIC_00007</td><td>pBIC_00007</td><td>247</td><td>110</td><td>0.45</td><td>-1.152</td><td>5E-12</td><td>Transposase DDE domain protein</td></tr><tr><td>pBIC_00008</td><td>pBIC_00008</td><td>256</td><td>104</td><td>0.412</td><td>-1.281</td><td>1.20E-13</td><td>Transposase (ISKpn6)</td></tr><tr><td>pBIC_00034</td><td>pBIC_00034</td><td>4,125</td><td>3,052</td><td>0.741</td><td>-0.432</td><td>3.15E-05</td><td>Conjugal transfert protein TraG</td></tr><tr><td>pBIC_00064</td><td>pBIC_00064</td><td>3,745</td><td>3,173</td><td>0.848</td><td>-0.238</td><td>0.047283777</td><td>Lytic Transglycosylase</td></tr><tr><td>pBIC_00068</td><td>higA</td><td>4,004</td><td>3,293</td><td>0.823</td><td>-0.28</td><td>0.013712564</td><td>Antitoxin HigA</td></tr><tr><td>pBIC_00069</td><td>pBIC_00069</td><td>3,600</td><td>2,784</td><td>0.774</td><td>-0.369</td><td>0.000846589</td><td>Hypothetical protein</td></tr><tr><td>pBIC_00097</td><td>pBIC_00097</td><td>323</td><td>204</td><td>0.635</td><td>-0.655</td><td>6.44E-06</td><td>Putative transposase</td></tr><tr><td>pBIC_00098</td><td>pBIC_00098</td><td>190</td><td>92</td><td>0.491</td><td>-1.027</td><td>5.31E-07</td><td>Helix-turn-helix domain protein</td></tr><tr><td>pBIC_00108</td><td>pBIC_00108</td><td>119</td><td>50</td><td>0.43</td><td>-1.218</td><td>4.94E-08</td><td>Hypothetical protein</td></tr><tr><td>pBIC_00109</td><td>yjcC</td><td>450</td><td>300</td><td>0.669</td><td>-0.579</td><td>7.65E-06</td><td>Putative membrane protein YjcC</td></tr><tr><td>pBIC_00168</td><td>pBIC_00168</td><td>3,282</td><td>2,429</td><td>0.741</td><td>-0.433</td><td>0.000196756</td><td>Diguanylate cyclase</td></tr><tr><td>pBIC_00169</td><td>pBIC_00169</td><td>1,063</td><td>840</td><td>0.791</td><td>-0.338</td><td>0.027066196</td><td>Hypothetical protein</td></tr><tr><td>pBIC_00170</td><td>clpC</td><td>29,463</td><td>22,147</td><td>0.753</td><td>-0.41</td><td>0.000216524</td><td>ATP-dependent Clp protease ATP-binding subunit ClpC</td></tr><tr><td>pBIC_00171</td><td>hspA</td><td>11,647</td><td>7,775</td><td>0.669</td><td>-0.58</td><td>4.35E-07</td><td>Heat shock protein</td></tr><tr><td>pBIC_00172</td><td>pBIC_00172</td><td>4,274</td><td>2,662</td><td>0.624</td><td>-0.68</td><td>1.51E-08</td><td>Helix-turn-helix domain protein</td></tr><tr><td>pBIC_00176</td><td>pBIC_00176</td><td>379</td><td>266</td><td>0.704</td><td>-0.507</td><td>0.001793481</td><td>Transposase (ISKpn26)</td></tr><tr><td>pBIC_00177</td><td>pBIC_00177</td><td>108</td><td>75</td><td>0.695</td><td>-0.525</td><td>0.029045775</td><td>Hypothetical protein</td></tr><tr><td>pBIC_00186</td><td>pBIC_00186</td><td>217</td><td>165</td><td>0.764</td><td>-0.388</td><td>0.040167906</td><td>Transposase (IS26)</td></tr><tr><td colspan="8"><strong>UP-REGULATED GENES</strong></td></tr><tr><td>pBIC_00010</td><td>pBIC_00010</td><td>1,807</td><td>2,107</td><td>1.165</td><td>0.22</td><td>0.048646696</td><td>Hypothetical protein</td></tr><tr><td>pBIC_00107</td><td>pBIC_00107</td><td>48</td><td>93</td><td>1.911</td><td>0.935</td><td>0.003155501</td><td>Transposase</td></tr><tr><td>pBIC_00129</td><td>mepM</td><td>937</td><td>1,288</td><td>1.374</td><td>0.458</td><td>1.94E-04</td><td>Murein DD-endopeptidase MepM</td></tr><tr><td>pBIC_00131</td><td>copA</td><td>2,404</td><td>3,018</td><td>1.253</td><td>0.326</td><td>0.003442412</td><td>Copper resistance protein A precursor</td></tr><tr><td>pBIC_00132</td><td>copB</td><td>725</td><td>939</td><td>1.291</td><td>0.369</td><td>0.035026611</td><td>Copper resistance protein B precursor</td></tr><tr><td>pBIC_00134</td><td>copD</td><td>1,107</td><td>1,373</td><td>1.239</td><td>0.309</td><td>0.008349118</td><td>Copper resistance protein D</td></tr><tr><td>pBIC_00135</td><td>copR</td><td>1,001</td><td>1,188</td><td>1.185</td><td>0.245</td><td>0.023543192</td><td>Transcriptional activator protein CopR</td></tr><tr><td>pBIC_00136</td><td>copS</td><td>1,205</td><td>1,491</td><td>1.237</td><td>0.307</td><td>0.009481833</td><td>Sensor kinase CopS</td></tr><tr><td>pBIC_00137</td><td>pcoE</td><td>336</td><td>423</td><td>1.255</td><td>0.328</td><td>0.030486108</td><td>Putative copper-binding protein PcoE precursor</td></tr></tbody></table>

Genes with FC values <0.5 are highlighted in gray.

## Imipenem Induces the Activation of Several Regulons Linked to Oxidative Stress in E. coli (pBIC1a)

We searched for genes that could be related to regulons induced by oxidative stress conditions in E. coli (Zheng et al., 2001; Blanchard et al., 2007; Seo et al., 2015). Activation of 9 genes that belonged to the SoxRS regulon was observed e.g., the superoxide dismutase sodA (FC = 2.1, adjusted *p*-value = 4.6e-20), the fumarate hydratase fumC (FC = 2.2, adjusted *p*-value = 8.3e-10), the aconitate hydratase acnA (FC = 1.9, adjusted *p*-value = 2.4e-13), the paraquat-inducible membrane protein A pqiA (FC = 1.4, adjusted *p*-value = 3.8e-06; Table S4). SoxR is a transcription factor containing [2Fe-2S] clusters, which is the sensor enhancing the expression of SoxS in response to oxidative stress (Lushchak, 2011). Then, SoxS regulates genes implicated in superoxide scavenging, carbon metabolism, DNA repair and xenobiotic efflux (Blanchard et al., 2007). Our transcriptomic

data revealed differential gene expression of other regulons linked to oxidative stress response in E. coli, such as the OxyR regulon (grxA, trxC, ...), the IscR regulon for Iron-Sulfur cluster binding and assembly (iscR, iscA, iscU, iscS, bfd, fdx...) and the Fur regulon for iron homeostasis and uptake machinery (fep, fepD, fepG, fecR, efeB, efeU, yncD...; Tab le S4).

## Activation of the Tricarboxylic Acid Cycle

Induction of oxidative stress after exposure to β-lactams has been observed in previous studies (Kaldalu et al., 2004; Belenky et al., 2015). Kohanski et al. (2007) and Dwyer et al. (2014) proposed a model explaining how bactericidal antibiotics were able to generate reactive oxygen species (ROS) in bacterial cells in addition to their main mode of action (Kohanski et al., 2007; Dwyer et al., 2014). According to these studies, interactions between β-lactams and Penicillin-Binding Proteins (PBPs) stimulate the oxidation of NADH
A Up-regulated genes

<figure><img src="image_6.png" /></figure>

B Down-regulated genes

<figure><img src="image_7.png" /></figure>

**FIGURE 4** | Biological process GO terms enrichment analysis in response to imipenem in *E. coli* (pBIC-1a). Scatterplots showing the non-redundant up (A) and down-regulated (B) GO terms significantly enriched after imipenem exposure. GO enrichment was determined by PANTHER web application and GO terms significantly enriched with adjusted *p*-values under 0.05 were summarized by using REVIGO. Only chromosomal genes were analyzed. The axes have no intrinsic meaning, similar GO terms remain close together in the plot. The bubble size indicates the frequency of the GO term. (A) upregulated genes revealed enrichment in energy production processes (aerobic respiration and electron transport chain), in ion transport and in cluster Fe-S assembly. (B) Downregulated genes were mostly associated to anaerobic respiration, generation of precursors of metabolites and energy. Both up and down-regulated genes summarization revealed an enrichment in oxidoreduction process.

through the electron transport chain, which is dependent on the tricarboxylic acid (TCA) cycle. Hyperactivation of the electron transport chain stimulates superoxide (O<sup>2</sup>-) formation, which overwhelms superoxide dismutase defenses and leads to the oxidation of iron-sulfur clusters ([4Fe-4S]<sup>2+</sup>) employed by numerous dehydratase enzymes. The release of ferrous (Fe<sup>2+</sup>) irons leads to hydroxyl radicals (OH*) via the Fenton reaction, which damage DNA, lipids and proteins and therefore contributes to antibiotic-induced cell death.

Activation of the TCA cycle and new iron-sulfur cluster synthesis have been shown to play critical roles in the initiation of the oxidative stress response after antibiotic exposure (Kohanski et al., 2010). Interestingly, our transcriptomic data revealed an overexpression of almost all TCA enzymes (Table S2; Figure 5).

A critical parameter to suppress Fenton reaction is the scavenging of H<sub>2</sub>O<sub>2</sub> generated in bacterial cells by sodA. This can occur through the activation of peroxidases like KatG, KatE, AhpC, AhpF, or the protein complex Hcp-Hcr being a putative hydroxylamine reductase or a peroxidase (Wolfe et al., 2002; Almeida et al., 2006). Here, only the two copies of ahpF were slightly up-regulated (with FC values <2). Of note, katG and hcp-hcr were even down-regulated. Our results suggest that no H<sub>2</sub>O<sub>2</sub> detoxification was highlighted by RNA-Seq after 10 min of imipenem exposure. Microarray data collected by Dwyer et al.

also observed down-regulation of *katG* and *hcp-hcr* in *E. coli* MG1655 exposed to ampicillin (Dwyer et al., 2014). Nevertheless, in this study, overexpression of *E. coli* *katG* and *ahpCF* partially protected the bacteria from β-lactams lethality, confirming the protective role these peroxidases.

## Genes Related to Fe-S Clusters Biogenesis and Iron Recruiting

In addition, we observed up-regulation of numerous Fe-S cluster binding proteins and Fe-S cluster assembly proteins (*iscS*, *iscU*, *iscA*, *yadR*, and *gntY*) indicating an activation of the Fe-S cluster biogenesis. Consistent with this, genes encoding proteins related to iron uptake and storage were regulated following imipenem exposure. Several iron permeases were up-regulated such as *efeU* (ferrous iron permease) (validated by RT-qPCR), *efeB* (exogenous heme iron acquisition), *fecI* (minor sigma factor that initiates transcription of ferric citrate transport genes), *fecR*, *fepB*, *fepD*, *fegG* (ferric enterobactin ABC transporter complex), and *yncD* (outer membrane protein involved in iron transport). In parallel, *bfd* gene encoding iron storage protein bacterioferritin-associated ferredoxin was up-regulated but not its partner coding the bacterioferritin Bfr. Bfd and Bfr can store iron as ferric iron to protect cells against ROS resulting from ferrous iron overload (Andrews et al., 2003).
<figure><img src="image_8.png" /></figure>

**FIGURE 5** | Metabolic response of *E. coli* pBIC-1a to imipenem exposure. Activation of the electron transport chain pathway dependent on the tricarboxylic acid (TCA) cycle was observed, and is thought to generate ROS according to the model proposed by Kohanski et al. and Dwyer et al. (Kohanski et al., 2007; Dwyer et al., 2014). We observed an up-regulation of most of the TCA cycle enzymes, of many genes involved in Iron-Sulfur cluster biogenesis and assembly, and of several metal ions permeases that are assumed to fuel the Fenton reaction. Detoxification process occured to limit protein damages and ROS accumulation with up-regulation of superoxide dismutase *sodA*. Several genes involved in bacterial cell wall homeostasis were disturbed. Activation of the plasmidic murein endopeptidase *mepM* might modify bacterial peptidoglycan and reflect an early activation of bacterial lysis. Transcriptional regulators are written in bold. ↑ means up-regulated and ↓ means down-regulated. LPS, lipopolysaccharide; OM, outer membrane; IM, inner membrane. Reactions of the TCA cycle: 1. Pyruvate dehydrogenase, 2. Citrate synthase, 3. Aconitase, 4. Isocitrate dehydrogenase, 5. α-ketoglutarate dehydrogenase, 6. SuccinylCoA thiokinase, 7. Succinate dehydrogenase, 8. Fumarate reductase, 9. Fumarase, 10. Malate dehydrogenase.

The uptaken $Fe^{2+}$ could also be reintegrated into Fe-S clusters to stabilize dehydratases and therefore be a metabolic feedback to counteract antibiotic toxicity (Py and Barras, 2010). YtfE is a protein known to recruit and integrate $Fe^{2+}$ to repair Fe-S cluster damaged proteins (Py and Barras, 2010). Here, *ytfE* was not induced in our conditions (FC = 0.816, *p*-value = 0.16), suggesting that no suppression of Fenton reaction occurred through *ytfE*. Therefore, $Fe^{2+}$ uptake could directly fuel Fenton reaction to produce the highly reactive HO\* and contribute to maintain oxidative stress.

## Genes Involved in Cell Wall

Imipenem is an antibiotic targeting bacterial cell wall biogenesis. Therefore, we looked at cell envelope components that could

be implicated in reducing entry of antibiotic or in resistance to osmotic stress.

Among the osmotic stress-inducible genes, *ybaY*, *bdm*, and *osmY* were up-regulated. YbaY is a lipoprotein associated to various stress responses like the oxidative stress response (Cheung et al., 2003), the RpoS-induced stress (Weber et al., 2005), and the osmotic stress response (Asakura and Kobayashi, 2009). Bdm is involved in flagellar biosynthesis, and its overexpression, regulated by the RcsB transcriptional activator, enhances biofilm formation (Sim et al., 2010; Kim et al., 2015). Here, RcsA/RcsB regulators that belong to the multicomponent RcsF/RcsC/RcsD/RcsA-RcsB phosphorelay system were indeed up-regulated upon imipenem addition (Gottesman and Stout, 1991; Majdalani and Gottesman, 2005). Any perturbation of the peptidoglycan activates the Rcs phosphorelay, and RcsB activity
enables bacteria to survive in the presence of the antibiotics (Laubacher and Ades, 2008), possibly due to a modification of bacterial cell wall through synthesis of the colanic acid capsular exopolysaccharide. Modification of expression of other genes belonging to the Rcs regulon and involved in flagellar biosynthesis was observed, like the downregulation of *flhC/flhD*.

Few genes associated to membrane components were down-regulated, like *ompW* coding for an outer-membrane porin. Expression of *ompW* was validated by RT-qPCR (Figure S4). OmpW has been implicated in the adaptation to stresses in various species, but its biological function is not yet fully characterized. *ompW* expression was decreased in *E. coli* in response to oxidative stress by an unknown mechanism (Blanchard et al., 2007). Recently, Xia et al. revealed a role in the carbon and energy metabolism in *E. coli* (Xiao et al., 2016). Expression of *ompW* is modulated by various regulators including RybB, a small RNA dependent on the alternative sigma factor $\sigma^E$, and the transcriptional regulator FNR (Fumarate and Nitrate Reduction) particularly important during anaerobic transition (Xiao et al., 2016). Here, both *rybB* and *fnr* were overexpressed after imipenem addition (FC = 1.546; *p*-value = 0.023 and FC = 2.38; *p*-value = 1.1e-46 respectively) and thus probably involved in *ompW* regulation.

In parallel to *ompW* down-regulation, several efflux pumps were up-regulated, *mdfA* (also known as *cmr*) *mdtK* and *sugE*, with FC values <2. These proteins are part of multi-drugs efflux pumps and could help bacteria to get rid of carbapenems and/or oxidative compounds.

## Small RNAs

Bacterial small RNAs are major actors for bacterial adaptation to various stresses and some of them are involved in regulatory circuits controlling antibiotic resistance (Felden and Cattoir, 2018). For instance, sRNA can control antimicrobial intracellular concentration by regulating the expression of several outer membrane proteins and efflux pumps. They can also regulate the expression of bacterial component like the LPS or enzymes involved in cell wall biosynthesis (Felden and Cattoir, 2018).

Following imipenem exposure, 14 sRNAs were differentially expressed (Table S5): 11 were up-regulated and 3 were down-regulated.

IsrB, also known as AzuC, was the most induced sRNA (FC = 10.334; *p*-value = 4.26e-80). No target has been clearly identified to date, but hydrogen peroxide treatment could also induce AzuC expression 4-fold in a previous study in *E. coli*, suggesting a role in oxidative stress response (Hemm et al., 2010).

*cyaR* (also known as *ryeE*) (FC = 1.93; *p*-values = 1.33e-12) and *omrA* (also known as *rygA*) (FC = 1.88; *p*-value = 6.3e-5) are known to decrease the expression of outer membrane proteins OmpX and OmpT, respectively (Vogel and Papenfort, 2006; Rau et al., 2015). However, it is not clear whether this down-regulation occurs at the translational level or through the decay of the mRNA targets (Vogel and Papenfort, 2006). Here, we could not observe a down-expression of *ompX* and *ompT* suggesting that base-pairing of both sRNAs with their targets does not affect RNA stability but would essentially interfere with the translation of the outer membrane proteins. OmpT

belongs to aspartyl proteases, involved in virulence by protecting cells from cationic antimicrobial peptides, activation of the anticoagulation pathway, and inflammatory responses (Hwang et al., 2007). OmpX mediates adhesion to eukaryotic cells and negatively regulates porins involved in β-lactam uptake like OmpK35 and OmpK36 in *Enterobacter aerogenes* (Dupont et al., 2004; Kim et al., 2010). Its expression also increases under high osmolarity conditions (Dupont et al., 2004). Inhibition of *ompX* translation through interaction with CyaR might indirectly increase carbapenem influx in bacterial cells.

Another target of OmrA is the transcriptional regulator CsgD that controls biofilm formation (Mika and Hengge, 2014). In biofilms, bacteria survive as persister cells and become more tolerant to antibiotics. By downregulating *csgD* and therefore biofilm formation, the expression of OmrA could increase susceptibility to antimicrobial agents.

*csrB* and *csrC* were less expressed in the presence of imipenem. These two sRNA regulate the carbon storage regulator *csrA*, which activates the expression of several metabolic enzymes involved in energy conversion (Weilbacher et al., 2003). One of the targeted enzymes is the Acetyl CoA synthase (FC = 1.34, *p*-value = 0.006), which when induced may contribute to fuel the TCA cycle by a secondary pathway.

RybA (also known as MntS) was also an up-regulated small RNA during imipenem exposure. As previously mentioned, RybA transcripts could be implicated in the up-regulation of the plasmid copper resistance operon (see section Imipenem-induced response at the plasmid pBIC1a level).

## Conclusions

Taken together, our work showed that imipenem induced an oxidative stress response in *E. coli* TOP10(pBIC1a) and perturbed several regulatory networks and cellular compounds (Figure 5). In particular, treatment with high concentrations resulted in the activation of the TCA cycle, the electron transport chain pathway and iron metabolism. Few enzymes were induced supposedly to limit the damages caused by oxidative stress, such as detoxification by superoxide dismutase SodA, and protein damage repair by thioredoxins and glutaredoxins (Figure 5). Whereas, the effect on chromosomal genes was highly significant, the impact of imipenem exposure on the transcription of plasmid genes was more limited. However, we could identify few differentially expressed genes that could be involved in resistance to these stressful conditions.

β-lactams are also known to trigger a SOS response in Gram-negative bacilli (Michel, 2005) but activation of the SOS response was not evidenced here. SOS response can be induced when bacteria are exposed to subinhibitory concentrations of antibiotics. In this study, high doses of imipenem were used instead to reproduce the effect of antibiotic concentrations reached in human serum. Another explanation might be that SOS response depends on the PBPs targeted by the β-lactam (Miller, 2004). After exposure to ampicillin, PBP3 inhibition causes filamentation and is known to stimulate the DpiAB two-component system, which activates the SOS response (Miller, 2004). At the opposite, imipenem preferentially targets PBP2
(Davies et al., 2008). By targeting preferentially some PBP, β-lactams seem to induce different stress responses.

Dissemination of KPC-encoding plasmids in Gram-negative bacteria poses a serious threat to medical treatment and patient management. To better understand the plasmid biology of a prototypical successful IncFII<sub>K</sub>-IncFIB *bla*<sub>KPC-2</sub>-carrying plasmid, a transcriptomic analysis was performed in *E. coli* TOP10, which constituted a well-characterized and annotated model. The choice of the genetic background is an important matter and has probably an impact on the expression of plasmid genes. Given the strong association between IncFII<sub>K</sub> plasmids and *K. pneumoniae* CG258, further studies using *K. pneumoniae* background are needed to understand the interaction between successful plasmids and successful clones. Here, genes involved in antibiotic resistance and basic plasmid function such as replication were highly expressed in broth culture. Carbapenem exposure did not influence the level of transcription of antimicrobial resistance genes. However, despite high transcription of *bla*<sub>KPC-2</sub> gene, we observed a major impact of imipenem on chromosomal genes and identified adaptive responses of *E. coli* to imipenem-induced oxidative stress.

# METHODS

## Bacterial Transformation

*K. pneumoniae* BIC-1 is a KPC-producing clinical strain isolated at Bicêtre Hospital in 2009 (Naas et al., 2010). This strain has been entirely sequenced using PacBio and Illumina sequencing (Jousset et al., 2018). The natural IncFII<sub>K</sub>-IncFIB plasmid carrying *bla*<sub>KPC-2</sub> gene, named pBIC1a was extracted from the BIC-1 strain by using the Kieser extraction method and subsequently analyzed by electrophoresis on a 0.7% agarose gel (Kieser, 1984). Plasmid pBIC1a was then introduced by electroporation into electrocompetent *E. coli* TOP10 (Invitrogen, Eragny, France) using a Gene Pulser II (Bio-Rad Laboratories. Marnes-la-Coquette, France). The transformant *E. coli* TOP10-pBIC1a was selected on plates supplemented with ticarcilline 50 mg/L (Sigma-Aldrich). The presence of pBIC1a was verified by PCR and by plasmid extraction followed by electrophoresis on a 0.7% agarose gel.

## pBIC1a Sequence Analysis

Annotation of pBIC1a was performed using PROKKA (Seemann, 2014). Mobile element proteins and insertion sequences were annotated with ISFinder (https://www.is-biotoul.fr/index.php; Siguier et al., 2006). Sequences of plasmids pKPN3 (NC_009649.1), pKpQIL (GU595196.1), pKpQIL-IT (NC_019155.1), pITO6-C07 (LT009688.1), pBK32179 (JX430448) and pGMI16-005-01 (NZ_CP028181.1) were used for genomic comparison. Plasmid alignment was performed using BRIG (Alikhan et al., 2011) and sequence similarity was studied using BLAST software (https://blast.ncbi.nlm.nih.gov/Blast.cgi). Allelic variant of the replicase genes were analyzed by replicon sequence typing using Plasmid Finder (https://cge.cbs.dtu.dk/services/PlasmidFinder/; Carattoli et al., 2014).

## Time-Kill Experiments

Time-kill experiments were performed to evaluate bacterial cell death after 10 min of imipenem exposure. Briefly, 10 ml of a log-phase culture (OD<sub>600</sub> at 0.5) of *E. coli* TOP10-pBIC1a was divided in two tubes of 5 ml. A final concentration of imipenem at 5 µg/ml was added in one tube. After 10 min, 10-fold serial dilutions were performed and 100 µl of diluted cultures were plated on agar plates without antibiotic. All these steps were performed in triplicates. The next day, colony-forming units (CFU) were counted.

## RNA Isolation and Purification

An overnight culture of *E. coli* TOP10-pBIC1a was diluted 100-fold in fresh BHI broth and allowed to grow to an OD<sub>600</sub> of 0.5. The culture was then split and imipenem (Imipenem-Cilastatin, Mylan, USA) was added to one of the cultures at a final concentration of 5 µg/ml. Ten min later, 2 ml of each culture were taken and mixed with RNA protect (Qiagen, Courtaboeuf, France) according to the manufacturer's recommendations. Cells were immediately pelleted by a 10 min centrifugation at 5,000 g. Biological triplicates were performed.

RNA was extracted using GeneJet RNA extraction kit (ThermoScientific, ThermoFisher Scientific, Villebon sur Yvette, France). DNAse treatment was performed with the dsDNAse kit (ThermoScientific) according to the manufacturer's protocol. RNA concentration was measured using Qubit RNA BR assay kit (Invitrogen. ThermoFisher Scientific), and RNA quality was evaluated using the Agilent 2100 Bioanalyzer RNA 6000 Nano Kit (Agilent Technologies, Santa Clara, CA). RIN >9 were obtained for all the RNA samples. Verification of complete removal of any contaminating DNA was performed via PCR amplification of *recA* housekeeping gene. A total of 2.5 µg of total RNA was treated with the Ribo-Zero rRNA Removal kit according to the manufacturer's instructions (Illumina, CA, USA) to remove 16S and 23S rRNAs. rRNA removal efficiency was then analyzed via the Bioanalyzer RNA 6000 Nano Kit. Ribodepleted RNA concentration was measured using Qubit™ RNA HS assay kit (Invitrogen).

## RNA Library and Sequencing

cDNA libraries were prepared using the NEBNext® Ultra Directional RNA Library Prep Kit for Illumina protocol (New England Biolabs, MA, USA). The quality of the cDNA was validated using the Agilent 2100 Bioanalyzer DNA1000 kit (Agilent Technologies) and quantity was determined with the Qubit dsDNA BR Assay (ThermoFisher Scientific). Sequencing of the library was performed with an Illumina HiSeq 2500 in single-read mode with 50 cycles.

## RT-qPCR

Eight up and down-regulated genes were selected and expression levels were analyzed by quantitative RT-PCR to validate the RNA-Seq data. The primers used are listed in Table S6. An input of 100 ng of total RNA, with no contaminating DNA was used. One-step RT-qPCR was performed using One-step RT-PCR kit (Qiagen) according to the manufacturer's protocol, and the reactions were carried out in 96-well plates with CFX96 Real-time PCR System
(BioRad, USA). All qRT-PCRs included an initial denaturation step of 30 s at 95°C, 35 cycles of denaturation (95°C/1s), annealing (52°C/5s), and extension (72°C/7s). Expression of each gene was normalized using *rpoB*, *idnT* and *yiaJ* genes as previously described (Zhou et al., 2011). Expression of these genes was not influenced by the presence of imipenem according to the RNA-Seq data. All reactions were carried out in triplicates. The relative abundance of gene transcripts among the imipenem treated group was calculated using the $2^{-\Delta\Delta Ct}$ method (Livak and Schmittgen, 2001).

## RNA-Seq Data Analysis and Differential Expression Analysis

Reads generated from strand-specific RNA-seq experiments were aligned to the genome of *E. coli* TOP10 genome (CP000948.1) and to that of pBIC1a sequence (CP022574.1) by using the software Bowtie (version 0.12.7) (Langmead et al., 2009). Reads that mapped in more than four different positions on the genome were discarded i.e., reads corresponding to rRNA. RNA-seq data were analyzed as described (Rosinski-Chupin et al., 2015) by using Rsamtools (version 1.26.2), GenomicAlignments (version 1.10.1). GenomicFeatures (version 1.26.4) and DESeq2 (version 1.14.1) and SARTools (Varet et al., 2016) in R 3.3.1. Only adjusted *p*-values were used and were obtained using the Benjamini-Hochberg correction for false discovery rate (Benjamini and Hochberg, 1995). Read count data were visually assessed using the Artemis genome viewer (Carver et al., 2008). The coverage (expressed as the average density of reads over 500 nucleotide fragments) was computed on each strand of pBICa sequence and visualized by using a heatmap representation (IGV software) (Robinson et al., 2011). The gene expression values were quantified in terms of reads per million (RPKM) defined as the total number of reads mapping to the feature divided by feature length (in kbp) normalized by the total number of reads (in millions) (Mortazavi et al., 2008).

The significantly up- and downregulated genes were analyzed for GO term enrichment separately with PANTHER (http://pantherdb.org/webservices/go/overrep.jsp; Mi et al., 2017) using batch upload, and the significantly enriched terms further explored on the REVIGO web application (http://revigo.irb.hr/) to identify and visualize relationships among the GO terms (Supek et al., 2011). The default settings were used on REVIGO with a list of results of medium size, and *E. coli* GO terms database used as reference.

The functions of genes and gene regulons were annotated using the following web databases: EcoCyc (http://biocyc.org/ECOLI/organism-summary; Keseler et al., 2017) and RegulonDB v8.2 (http://regulondb.ccg.unam.mx/index.jsp). Regulation for genes in the SoxRS, OxyR, IscR and Fur regulons were identified as annotated in RegulonDB v8.2 (Gama-Castro et al., 2016).

## ACCESSION NUMBERS

Plasmid pBIC1a has been deposited in NCBI database under accession number CP022574.1, and RNA-Seq data in Array Express under accession number E-MTAB-7190.

## AUTHOR CONTRIBUTIONS

AJ was in charge of the study design, library preparations, RNA-seq experiments, data analysis, and writing of the manuscript. IRC performed statistical analysis, data analysis and the proofreading of the manuscript. JT realized qRT-PCR and Time kill experiments. RB worked on study design, data analysis and the proofreading of the manuscript. PG and TN were in charge of the study design, and the proofreading of the manuscript.

## FUNDING

This work was supported by the Assistance Publique—Hôpitaux de Paris, by a grant from the Université Paris-Sud (EA7361), by the LabEx LERMIT supported by a grant from the French National Research Agency (ANR-10-LABX-33), and by a project of ANR LabEx IBEID. This work was also funded in part by a grant from Joint Program Initiative on Antimicrobial Resistance (ANR-14-JAMR-0002).

## SUPPLEMENTARY MATERIAL

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

**Figure S1** | Time-kill analysis after imipenem exposure. Bacteria were diluted using 10-fold serial dilutions and Colony-Forming Units (CFUs) were counted before (T0) and 10 min after imipenem addition (T10). The percentage of surviving bacteria is represented. The experiments were performed in triplicates.

**Figure S2** | Differential expression of *E. coli* (pBIC1a) genes in presence of imipenem (imi). Condition without antibiotic was used as control (Ct). 479 significantly differentially expressed genes were found (*p*-value <0.05). (A) Plots of the principal component analysis (PCA). All experiments were performed with biological triplicates. (B) Volcano plot with 325 genes with diminished expression and 154 genes with increased expression in presence of imipenem are specified.

**Figure S3** | Transcriptional landscape of *Tn4401a* bracketing *bla*<sub>KPC-2</sub> carbapenemase gene. (A) Strand-oriented read mapping of *Tn4401a* using Artemis software (Carver et al., 2008). No clear transcription termination of *bla*<sub>KPC-2</sub> transcripts was observed, and this transcription covered in antisense the transposase gene of *ISKpn6*. CDS and their orientations are indicated by colored arrows in the lower part. The six reading frames are indicated with vertical bars representing stop codon. Green reads represent multiple identical reads merged and black reads represent unique read mapped. (B) Mapping of the reads indicates promoter p2 as the main promoter. Promoters are indicated by black boxes and named according to the literature (Girlich et al., 2017). Inverted repeat right (IRR) of *ISKpn7* is indicated by horizontal black arrow.

**Figure S4** | RT-qPCR validation of eight representative genes with significantly differential expression after imipenem-induced stress. RT-qPCR data are represented as mean ± SD from three biological replicates.

**Table S1** | Alignment of raw data on *E. coli* TOP10(pBIC1a) sequence.

**Table S2** | Differentially expressed genes after imipenem addition in broth culture (5µg/ml).

**Table S3** | Effet of imipenem on pBIC1a genes.

**Table S4** | SoxRS/OxyRS/IscR/Fur-regulated gene expression changes.

**Table S5** | Differentially expressed regulatory RNAs.

**Table S6** | Primers used for RT-PCR.
REFERENCES

Alikhan, N.-F., Petty, N. K., Ben Zakour, N. L., and Beatson, S. A. (2011). BLAST ring image generator (BRIG): simple prokaryote genome comparisons. BMC Genomics 12:402. doi: 10.1186/1471-2164-12-402

Almeida, C. C., Romão, C. V., Lindley, P. F., Teixeira, M., and Saraiva, L. M. (2006). The role of the hybrid cluster protein in oxidative stress defense. *J. Biol. Chem.* 281, 32445-32450. doi: 10.1074/jbc.M605888200

Andrews, S. C., Robinson, A. K., and Rodriguez-Quiñones, F. (2003). Bacterial iron homeostasis. *FEMS Microbiol. Rev.* 27, 215-237. doi: 10.1016/S0168-6445(03)00055-X

Aravind, L., Anantharaman, V., Balaji, S., Babu, M. M., and Iyer, L. M. (2005). The many faces of the helix-turn-helix domain: transcription regulation and beyond. *FEMS Microbiol. Rev.* 29, 231-262. doi: 10.1016/j.femsre.2004.12.008

Asakura, Y., and Kobayashi, I. (2009). From damaged genome to cell surface: transcriptome changes during bacterial cell death triggered by loss of a restriction-modification gene complex. *Nucleic Acids Res.* 37, 3021-3031. doi: 10.1093/nar/gk148

Baraniak, A., Grabowska, A., Izdebski, R., Fiett, J., Herda, M., Bojarska, K., et al. (2011). Molecular characteristics of KPC-producing Enterobacteriaceae at the early stage of their dissemination in Poland, 2008–2009. *Antimicrob. Agents Chemother.* 55, 5493–5499. doi: 10.1128/AAC.05118-11

Belenky, P., Ye, J. D., Porter, C. B. M., Cohen, N. R., Lobritz, M. A., Ferrante, T., et al. (2015). Bactericidal antibiotics induce toxic metabolic perturbations that lead to cellular damage. *Cell Rep.* 13, 968-980. doi: 10.1016/j.celrep.2015.09.059

Benjamini, Y., and Hochberg, Y. (1995). Controlling the false discovery rate: a practical and powerful approach to multiple testing. *J. R. Stat. Soc. Ser. B Stat. Methodol.* 57, 289–300.

Blanchard, J. L., Wholey, W.-Y., Conlon, E. M., and Pomposiello, P. J. (2007). Rapid changes in gene expression dynamics in response to superoxide reveal SoxRS-dependent and independent transcriptional networks. *PLOS ONE* 2:e1186. doi: 10.1371/journal.pone.0001186

Bondarczuk, K., and Piotrowska-Seget, Z. (2013). Molecular basis of active copper resistance mechanisms in Gram-negative bacteria. *Cell Biol. Toxicol.* 29, 397-405. doi: 10.1007/s10565-013-9262-1

Bonomo, R. A., Burd, E. M., Conly, J., Limbago, B. M., Poirel, L., Segre, J. A., et al. (2017). Carbapenemase-producing organisms: a global scourge. *Clin. Infect. Dis.* 66, 1290-1297. doi: 10.1093/cid/cix893

Buckner, M. M. C., Saw, H. T. H., Osagie, R. N., McNally, A., Ricci, V., Wand, M. E., et al. (2018). Clinically relevant plasmid-host interactions indicate that transcriptional and not genomic modifications ameliorate fitness costs of *Klebsiella pneumoniae* carbapenemase-carrying plasmids. *MBio* 9:e02303-17. doi: 10.1128/MBio.02303-17

Carattoli, A., Zankari, E., García-Fernández, A., Voldby Larsen, M., Lund, O., Villa, L., et al. (2014). *In silico* detection and typing of plasmids using PlasmidFinder and plasmid multilocus sequence typing. *Antimicrob. Agents Chemother.* 58, 3895-3903. doi: 10.1128/AAC.02412-14

Carver, T., Berriman, M., Tivey, A., Patel, C., Böhme, U., Barrell, B. G., et al. (2008). Artemis and ACT: viewing, annotating and comparing sequences stored in a relational database. *Bioinformatics* 24, 2672-2676. doi: 10.1093/bioinformatics/btn529

Chen, L., Chavda, K. D., Melano, R. G., Jacobs, M. R., Koll, B., Hong, T., et al. (2014a). Comparative genomic analysis of KPC-encoding pKpQIL-like plasmids and their distribution in New Jersey and New York Hospitals. *Antimicrob. Agents Chemother.* 58, 2871-2877. doi: 10.1128/AAC.00120-14

Chen, L., Chavda, K. D., Melano, R. G., Jacobs, M. R., Levi, M. H., Bonomo, R. A., et al. (2013). Complete sequence of a bla(KPC-2)-harboring IncFII(K1) plasmid from a *Klebsiella pneumoniae* sequence type 258 strain. *Antimicrob. Agents Chemother.* 57, 1542-1545. doi: 10.1128/AAC.02332-12

Chen, L., Mathema, B., Chavda, K. D., DeLeo, F. R., Bonomo, R. A., and Kreiswirth, B. N. (2014b). Carbapenemase-producing Klebsiella pneumoniae: molecular and genetic decoding. *Trends Microbiol.* 22, 686-696. doi: 10.1016/j.tim.2014.09.003

Cheruvanky, A., Stoesser, N., Sheppard, A. E., Crook, D. W., Hoffman, P. S., Weddle, E., et al. (2017). Enhanced Klebsiella pneumoniae carbapenemase expression from a novel Tn4401 deletion. *Antimicrob. Agents Chemother.* 61:e00025-17. doi: 10.1128/AAC.00025-17

Cheung, K. J., Badarinarayana, V., Selinger, D. W., Janse, D., and Church, G. M. (2003). A microarray-based antibiotic screen identifies a regulatory role for supercoiling in the osmotic stress response of *Escherichia coli*. *Genome Res.* 13, 206-215. doi: 10.1101/gr.401003

Chillappagari, S., Seubert, A., Trip, H., Kuipers, O. P., Marahiel, M. A., and Miethke, M. (2010). Copper stress affects iron homeostasis by destabilizing iron-sulfur cluster formation in *Bacillus subtilis*. *J. Bacteriol.* 192, 2512-2524. doi: 10.1128/JB.00058-10

Cuzon, G., Naas, T., Correa, A., Quinn, J. P., Villegas, M.-V., and Nordmann, P. (2013). Dissemination of the KPC-2 carbapenemase in non-Klebsiella pneumoniae enterobacterial isolates from Colombia. *Int. J. Antimicrob. Agents* 42, 59-62. doi: 10.1016/j.ijantimicag.2013.04.002

Cuzon, G., Naas, T., Villegas, M.-V., Correa, A., Quinn, J. P., and Nordmann, P. (2011). Wide dissemination of pseudomonas aeruginosa producing lactamase blaKPC-2 gene in Colombia. *Antimicrob. Agents Chemother.* 55, 5350-5353. doi: 10.1128/AAC.00297-11

Daikos, G. L., and Markogiannakis, A. (2011). Carbapenemase-producing Klebsiella pneumoniae: (when) might we still consider treating with carbapenems? *Clin. Microbiol. Infect.* 17, 1135-1141. doi: 10.1111/j.1469-0691.2011.03553.x

Daikos, G. L., Tsaousi, S., Tzouvelekis, L. S., Anyfantis, I., Psichogiou, M., Argyropoulou, A., et al. (2014). Carbapenemase-producing Klebsiella pneumoniae bloodstream infections: lowering mortality by antibiotic combination schemes and the role of carbapenems. *Antimicrob. Agents Chemother.* 58, 2322-2328. doi: 10.1128/AAC.02 166-13

Davies, T. A., Shang, W., Bush, K., and Flamm, R. K. (2008). Affinity of doripenem and comparators to penicillin-binding proteins in *Escherichia coli* and Pseudomonas aeruginosa. *Antimicrob. Agents Chemother.* 52, 1510-1512. doi: 10.1128/AAC.01529-07

Deleo, F. R., Chen, L., Porcella, S. F., Martens, C. A., Kobayashi, S. D., Porter, A. R., et al. (2014). Molecular dissection of the evolution of carbapenem-resistant multilocus sequence type 258 Klebsiella pneumoniae. *Proc. Natl. Acad. Sci. U.S.A.* 111, 4988-4993. doi: 10.1073/pnas.1321364111

Doumith, M., Findlay, J., Hirani, H., Hopkins, K. L., Livermore, D. M., Dodgson, A., et al. (2017). Major role of pKpQIL-like plasmids in the early dissemination of KPC-type carbapenemases in the UK. *J. Antimicrob. Chemother.* 72, 2241-2248. doi: 10.1093/jac/dkx141

Dupont, M., Dé, E., Chollet, R., Chevalier, J., and Pagès, J.-M. (2004). Enterobacter aerogenes OmpX, a cation-selective channel mar- and osmo-regulated. *FEBS Lett.* 569, 27-30. doi: 10.1016/j.febslet.2004.05.047

Dwyer, D. J., Belenky, P. A., Yang, J. H., MacDonald, I. C., Martell, J. D., Takahashi, N., et al. (2014). Antibiotics induce redox-related physiological alterations as part of their lethality. *Proc. Natl. Acad. Sci. U.S.A.* 111, E2100-E2109. doi: 10.1073/pnas.1401876111

Felden, B., and Cattoir, V. (2018). Bacterial adaptation to antibiotics through regulatory RNAs. *Antimicrob. Agents Chemother.* 62:e02503-17. doi: 10.1128/AAC.02503-17

Finley, P. J., Norton, R., Austin, C., Mitchell, A., Zank, S., and Durham, P. (2015). Unprecedented silver resistance in clinically isolated enterobacteriaceae: major implications for burn and wound management. *Antimicrob. Agents Chemother.* 59, 4734-4741. doi: 10.1128/AAC.00026-15

Gama-Castro, S., Salgado, H., Santos-Zavaleta, A., Ledezma-Tejeida, D., Muñiz-Rascado, L., García-Sotelo, J. S., et al. (2016). RegulonDB version 9.0: high-level integration of gene regulation, coexpression, motif clustering and beyond. *Nucleic Acids Res.* 44, D133-D143. doi: 10.1093/nar/gkv1156

García-Fernández, A., Villa, L., Carta, C., Venditti, C., Giordano, A., Venditti, M., et al. (2012). Klebsiella pneumoniae ST258 producing KPC-3 identified in italy carries novel plasmids and OmpK36/OmpK35 porin variants. *Antimicrob. Agents Chemother.* 56, 2143-2145. doi: 10.1128/AAC.05308-11

Gerstle, K., Klätschke, K., Hahn, U., and Piganeau, N. (2012). The small RNA RybA regulates key-genes in the biosynthesis of aromatic amino acids under peroxide stress in *E. coli*. *RNA Biol.* 9, 458-468. doi: 10.4161/rna. 19065

Girlich, D., Bonnin, R. A., Jousset, A., and Naas, T. (2017). Promoter characterization and expression of the blaKPC-2 gene in *Escherichia coli*, Pseudomonas aeruginosa and Acinetobacter baumannii. *J. Antimicrob. Chemother.* 72, 1597-1601. doi: 10.1093/jac/dkx044
Gottesman, S., and Stout, V. (1991). Regulation of capsular polysaccharide synthesis in *Escherichia coli* K12. *Mol. Microbiol.* **5**, 1599–1606.

Hemm, M. R., Paul, B. J., Miranda-Rios, J., Zhang, A., Soltanzad, N., and Storz, G. (2010). Small stress response proteins in *Escherichia coli*: proteins missed by classical proteomic studies. *J. Bacteriol.* **192**, 46–58. doi: 10.1128/JB.00872-09

Hidalgo-Grass, C., Warburg, G., Temper, V., Benenson, S., Moses, A. E., Block, C., et al. (2012). KPC-9, a novel carbapenemase from clinical specimens in Israel. *Antimicrob. Agents Chemother.* **56**, 6057–6059. doi: 10.1128/AAC.01156-12

Hwang, B.-Y., Varadarajan, N., Li, H., Rodriguez, S., Iverson, B. L., and Georgiou, G. (2007). Substrate specificity of the *Escherichia coli* outer membrane protease OmpP. *J. Bacteriol.* **189**, 522–530. doi: 10.1128/JB.01493-06

Jousset, A. B., Bonnin, R. A., Rosinski-Chupin, I., Girlich, D., Cuzon, G., Cabanel, N., et al. (2018). 4.5 years within-patient evolution of a colistin resistant KPC-producing Klebsiella pneumoniae ST258. *Clin. Infect. Dis.* **67**, 1388–1394. doi: 10.1093/cid/ciy293

Kaldalu, N., Mei, R., and Lewis, K. (2004). Killing by ampicillin and ofloxacin induces overlapping changes in *Escherichia coli* transcription profile. *Antimicrob. Agents Chemother.* **48**, 890–896. doi: 10.1128/AAC.48.3.890-896.2004

Keseler, I. M., Mackie, A., Santos-Zavaleta, A., Billington, R., Bonavides-Martinez, C., Caspi, R., et al. (2017). The EcoCyc database: reflecting new knowledge about *Escherichia coli* K-12. *Nucleic Acids Res.* **45**, D543–D550. doi: 10.1093/nar/gkw1003

Kieser, T. (1984). Factors affecting the isolation of CCC DNA from *Streptomyces lividans* and *Escherichia coli*. *Plasmid* **12**, 19–36. doi: 10.1016/0147-619X(84)90063-5

Kim, J.-S., Kim, Y. J., Seo, S., Seong, M.-J., and Lee, K. (2015). Functional role of bdm during flagella biogenesis in *Escherichia coli*. *Curr. Microbiol.* **70**, 369–373. doi: 10.1007/s00284-014-0729-y

Kim, K., Kim, K.-P., Choi, J., Lim, J.-A., Lee, J., Hwang, S., et al. (2010). Outer membrane proteins A (OmpA) and X (OmpX) are essential for basolateral invasion of *Cronobacter sakazakii*. *Appl. Environ. Microbiol.* **76**, 5188–5198. doi: 10.1128/AEM.02498-09

Kohanski, M. A., Dwyer, D. J., and Collins, J. J. (2010). How antibiotics kill bacteria: from targets to networks. *Nat. Rev. Microbiol.* **8**, 423–435. doi: 10.1038/nrmicro2333

Kohanski, M. A., Dwyer, D. J., Hayete, B., Lawrence, C. A., and Collins, J. J. (2007). A common mechanism of cellular death induced by bactericidal antibiotics. *Cell* **130**, 797–810. doi: 10.1016/j.cell.2007.06.049

Landman, D., Urban, C., Bäcker, M., Kelly, P., Shah, N., Babu, E., et al. (2010). Susceptibility profiles, molecular epidemiology, and detection of KPC-producing *Escherichia coli* isolates from the New York City vicinity. *J. Clin. Microbiol.* **48**, 4604–4607. doi: 10.1128/JCM.01143-10

Lang, K. S., Danzeisen, J. L., Xu, W., and Johnson, T. J. (2012). Transcriptome mapping of pAR060302, a blaCMY-2-positive broad-host-range IncA/C plasmid. *Appl. Environ. Microbiol.* **78**, 3379–3386. doi: 10.1128/AEM.07199-11

Langmead, B., Trapnell, C., Pop, M., and Salzberg, S. L. (2009). Ultrafast and memory-efficient alignment of short DNA sequences to the human genome. *Genome Biol.* **10**:R25. doi: 10.1186/gb-2009-10-3-r25

Lau, A. F., Wang, H., Weingarten, R. A., Drake, S. K., Suffredini, A. F., Garfield, M. K., et al. (2014). A rapid matrix-assisted laser desorption ionization-time of flight mass spectrometry-based method for single-plasmid tracking in an outbreak of carbapenem-resistant Enterobacteriaceae. *J. Clin. Microbiol.* **52**, 2804–2812. doi: 10.1128/JCM.00694-14

Laubacher, M. E., and Ades, S. E. (2008). The Rcs phosphorelay is a cell envelope stress response activated by peptidoglycan stress and contributes to intrinsic antibiotic resistance. *J. Bacteriol.* **190**, 2065–2074. doi: 10.1128/JB.01740-07

Lee, G. C., and Burgess, D. S. (2012). Treatment of Klebsiella Pneumoniae carbapenemase (KPC) infections: a review of published case series and case reports. *Ann. Clin. Microbiol. Antimicrob.* **11**:32. doi: 10.1186/1476-0711-11-32

Livak, K. J., and Schmittgen, T. D. (2001). Analysis of relative gene expression data using real-time quantitative PCR and the 2(-Delta Delta C(T)) method. *Methods* **25**, 402–408. doi: 10.1006/meth.2001.1262

Lushchak, V. I. (2011). Adaptive response to oxidative stress: bacteria, fungi, plants and animals. *Comp. Biochem. Physiol. C Toxicol. Pharmacol.* **153**, 175–190. doi: 10.1016/j.cbpc.2010.10.004

Majdalani, N., and Gottesman, S. (2005). The Rcs phosphorelay: a complex signal transduction system. *Annu. Rev. Microbiol.* **59**, 379–405. doi: 10.1146/annurev.micro.59.050405.101230

Mi, H., Huang, X., Muruganujan, A., Tang, H., Mills, C., Kang, D., et al. (2017). PANTHER version 11: expanded annotation data from gene ontology and reactome pathways, and data analysis tool enhancements. *Nucleic Acids Res.* **45**, D183–D189. doi: 10.1093/nar/gkw1138

Michel, B. (2005). After 30 years of study, the bacterial SOS response still surprises us. *PLoS Biol.* **3**:e255. doi: 10.1371/journal.pbio.0030255

Mika, F., and Hengge, R. (2014). Small RNAs in the control of RpoS, CsgD, and biofilm architecture of *Escherichia coli*. *RNA Biol.* **11**, 494–507. doi: 10.4161/rna.28867

Miller, C. (2004). SOS response induction by lactams and bacterial defense against antibiotic lethality. *Science* **305**, 1629–1631. doi: 10.1126/science.1101630

Mortazavi, A., Williams, B. A., McCue, K., Schaeffer, L., and Wold, B. (2008). Mapping and quantifying mammalian transcriptomes by RNA-Seq. *Nat. Methods* **5**, 621–628. doi: 10.1038/nmeth.1226

Naas, T., Cuzon, G., Babics, A., Fortineau, N., Boytchev, I., Gayral, F., et al. (2010). Endoscopy-associated transmission of carbapenem-resistant Klebsiella pneumoniae producing KPC-2 β-lactamase. *J. Antimicrob. Chemother.* **65**, 1305–1306. doi: 10.1093/jac/dkq117

Naas, T., Cuzon, G., Truong, H.-V., and Nordmann, P. (2012). Role of ISKp7n and deletions in blaKPC gene expression. *Antimicrob. Agents Chemother.* **56**, 4753–4759. doi: 10.1128/AAC.00334-12

Nagy, Z., and Chandler, M. (2004). Regulation of transposition in bacteria. *Res. Microbiol.* **155**, 387–398. doi: 10.1016/j.resmic.2004.01.008

Navon-Venezia, S., Kondratyeva, K., and Carattoli, A. (2017). Klebsiella pneumoniae: a major worldwide source and shuttle for antibiotic resistance. *FEMS Microbiol. Rev.* **41**, 252–275. doi: 10.1093/femsre/fux013

Papagiannitsis, C. C., Pilato, V. D., Giani, T., Giakkoupi, P., Riccobono, E., Landini, G., et al. (2016). Characterization of KPC-encoding plasmids from two endemic settings, Greece and Italy. *J. Antimicrob. Chemother.* **71**, 2824–2830. doi: 10.1093/jac/dkw227

Partridge, S. R. (2014). Tn4401 carrying blaKPC is inserted within another insertion in pKpQIL and related plasmids. *J. Clin. Microbiol.* **52**, 4448–4449. doi: 10.1128/JCM.02426-14

Py, B., and Barras, F. (2010). Building Fe-S proteins: bacterial strategies. *Nat. Rev. Microbiol.* **8**, 436–446. doi: 10.1038/nrmicro2356

Rademacher, C., and Masepohl, B. (2012). Copper-responsive gene regulation in bacteria. *Microbiology* **158**, 2451–2464. doi: 10.1099/mic.0.058487-0

Rau, M. H., Bojanovič, K., Nielsen, A. T., and Long, K. S. (2015). Differential expression of small RNAs under chemical stress and fed-batch fermentation in *E. coli*. *BMC Genomics* **16**:1051. doi: 10.1186/s12864-015-2231-8

Ren, S., Li, Q., Xie, L., and Xie, J. (2017). Molecular mechanisms underlying the function diversity of ArsR family metalloregulator. *Crit. Rev. Eukaryot. Gene Expr.* **27**, 19–35. doi: 10.1615/CritRevEukaryotGeneExpr.2016018476

Robinson, J. T., Thorvaldsdóttir, H., Winckler, W., Guttman, M., Lander, E. S., Getz, G., et al. (2011). Integrative genomics viewer. *Nat. Biotechnol.* **29**, 24–26. doi: 10.1038/nbt.1754

Rojas, L. J., Salim, M., Cober, E., Richter, S. S., Perez, F., Salata, R. A., et al. (2017). Colistin resistance in carbapenem-resistant Klebsiella pneumoniae: laboratory detection and impact on mortality. *Clin. Infect. Dis.* **64**, 711–718. doi: 10.1093/cid/ciw805

Rojas, L. J., Weinstock, G. M., De La Cadena, E., Diaz, L., Rios, R., Hanson, B. M., et al. (2018). An analysis of the epidemic of Klebsiella pneumoniae carbapenemase-producing K. pneumoniae: convergence of two evolutionary mechanisms creates the “Perfect Storm.” *J. Infect. Dis.* **217**, 82–92. doi: 10.1093/infdis/jix524

Rosinski-Chupin, I., Sauvage, E., Sismeiro, O., Villain, A., Da Cunha, V., Caliot, M.-E., et al. (2015). Single nucleotide resolution RNA-seq uncovers new regulatory mechanisms in the opportunistic pathogen *Streptococcus agalactiae*. *BMC Genomics* **16**:419. doi: 10.1186/s12864-015-1583-4

Roth, A. L., Kurpiel, P. M., Lister, P. D., and Hanson, N. D. (2011). *bla* KPC RNA expression correlates with two transcriptional start sites but not always with gene copy number in four genera of gram-negative pathogens. *Antimicrob. Agents Chemother.* **55**, 3936–3938. doi: 10.1128/AAC.01109-10
Ruiz de los Mozos, I., Vergara-Irigaray, M., Segura, V., Villanueva, M., Bitarte, N., Saramago, M., et al. (2013). Base pairing interaction between 5'- and 3'-UTRs controls icaR mRNA translation in *Staphylococcus aureus*. *PLoS Genet.* 9:e1004001. doi: 10.1371/journal.pgen.1004001

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

Seo, S. W., Kim, D., Szubin, R., and Palsson, B. O. (2015). Genome-wide reconstruction of OxyR and SoxRS transcriptional regulatory networks under oxidative stress in *Escherichia coli* K-12 MG1655. *Cell Rep.* 12, 1289–1299. doi: 10.1016/j.celrep.2015.07.043

Signs, S. A., Tan, J. S., Salstrom, S. J., and File, T. M. (1992). Pharmacokinetics of imipenem in serum and skin window fluid in healthy adults after intramuscular or intravenous administration. *Antimicrob. Agents Chemother.* 36, 1400–1403.

Siguier, P., Perochon, J., Lestrade, L., Mahillon, J., and Chandler, M. (2006). ISfinder: the reference centre for bacterial insertion sequences. *Nucleic Acids Res.* 34, D32–D36. doi: 10.1093/nar/gkj014

Sim, S.-H., Yeom, J.-H., Shin, C., Song, W.-S., Shin, E., Kim, H.-M., et al. (2010). Escherichia coli ribonuclease III activity is downregulated by osmotic stress: consequences for the degradation of bdm mRNA in biofilm formation. *Mol. Microbiol.* 75, 413–425. doi: 10.1111/j.1365-2958.2009.06986.x

Singlas, E. (1990). Clinical pharmacokinetics of an imipenem-cilastatin combination. *Presse Med.* 19, 598–605.

Supek, F., Bošnjak, M., Škunca, N., and Šmuc, T. (2011). REVIGO summarizes and visualizes long lists of gene ontology terms. *PLoS ONE* 6:e21800. doi: 10.1371/journal.pone.0021800

Tumbarello, M., Viale, P., Viscoli, C., Trecarichi, E. M., Tumietto, F., Marchese, A., et al. (2012). Predictors of mortality in bloodstream infections caused by Klebsiella pneumoniae carbapenemase-producing *K. pneumoniae*: importance of combination therapy. *Clin. Infect. Dis.* 55, 943–950. doi: 10.1093/cid/cis588

Uehara, T., Dinh, T., and Bernhardt, T. G. (2009). LytM-domain factors are required for daughter cell separation and rapid ampicillin-induced lysis in *Escherichia coli*. *J. Bacteriol.* 191, 5094–5107. doi: 10.1128/JB.00505-09

van Duin, D., Lok, J. J., Earley, M., Cober, E., Richter, S. S., Perez, F., et al. (2017). Colistin vs. Ceftazidime-avibactam in the treatment of infections due to carbapenem-resistant enterobacteriaceae. *Clin. Infect. Dis.* 66, 163–171. doi: 10.1093/cid/cix783

Varet, H., Brillet-Guéguen, L., Coppée, J.-Y., and Dillies, M.-A. (2016). SARTools: a DESeq2- and EdgeR-based R pipeline for comprehensive differential analysis of RNA-SEQ data. *PLOS ONE* 11:e0157022. doi: 10.1371/journal.pone.0157022

Villa, L., García-Fernández, A., Fortini, D., and Carattoli, A. (2010). Replicon sequence typing of IncF plasmids carrying virulence and resistance determinants. *J. Antimicrob. Chemother.* 65, 2518–2529. doi: 10.1093/jac/dkq347

Visweswaran, G. R. R., Steen, A., Leenhouts, K., Szeliga, M., Ruban, B., Hesseling-Meinders, A., et al. (2013). AcmD, a homolog of the major autolysin AcmA of Lactococcus lactis, binds to the cell wall and contributes to cell separation and autolysis. *PLoS ONE* 8:e72167. doi: 10.1371/journal.pone.0072167

Vogel, J., and Papenfort, K. (2006). Small non-coding RNAs and the bacterial outer membrane. *Curr. Opin. Microbiol.* 9, 605–611. doi: 10.1016/j.mib.2006.10.006

Weber, H., Polen, T., Heuveling, J., Wendisch, V. F., and Hengge, R. (2005). Genome-wide analysis of the general stress response network in *Escherichia coli*: σS-dependent genes, promoters, and sigma factor selectivity. *J. Bacteriol.* 187, 1591–1603. doi: 10.1128/JB.187.5.1591-1603. 2005

Weilbacher, T., Suzuki, K., Dubey, A. K., Wang, X., Gudapaty, S., Morozov, I., et al. (2003). A novel sRNA component of the carbon storage regulatory system of *Escherichia coli*. *Mol. Microbiol.* 48, 657–670. doi: 10.1046/j.1365-2958.2003.03459.x

Weisenberg, S. A., Morgan, D. J., Espinal-Witter, R., and Larone, D. H. (2009). Clinical outcomes of patients with Klebsiella pneumoniae carbapenemase-producing *K. pneumoniae* after treatment with imipenem or meropenem. *Diagn. Microbiol. Infect. Dis.* 64, 233–235. doi: 10.1016/j.diagmicrobio.2009.02.004

Wolfe, M. T., Heo, J., Garavelli, J. S., and Ludden, P. W. (2002). Hydroxylamine reductase activity of the hybrid cluster protein from *Escherichia coli*. *J. Bacteriol.* 184, 5898–5902. doi: 10.1128/JB.184.21.5898-5902.2002

Wright, M. S., Perez, F., Brinkac, L., Jacobs, M. R., Kaye, K., Cober, E., et al. (2014). Population structure of KPC-producing Klebsiella pneumoniae isolates from midwestern U.S. Hospitals. *Antimicrob. Agents Chemother.* 58, 4961–4965. doi: 10.1128/AAC.00125-14

Xiao, M., Lai, Y., Sun, J., Chen, G., and Yan, A. (2016). Transcriptional regulation of the outer membrane porin gene ompW reveals its physiological role during the transition from the aerobic to the anaerobic lifestyle of *Escherichia coli*. *Front. Microbiol.* 7:799. doi: 10.3389/fmicb.2016.00799

Xiao, M., Zhu, X., Fan, F., Xu, H., Tang, J., Qin, Y., et al. (2017). Osmotolerance in *Escherichia coli* is improved by activation of copper efflux genes or supplementation with sulfur-containing amino acids. *Appl. Environ. Microbiol.* 83:e03050-16. doi: 10.1128/AEM.03050-16

Zheng, M., Wang, X., Templeton, L. J., Smulski, D. R., LaRossa, R. A., and Storz, G. (2001). DNA microarray-mediated transcriptional profiling of the *Escherichia coli* response to hydrogen peroxide. *J. Bacteriol.* 183, 4562–4570. doi: 10.1128/JB.183.15.4562-4570.2001

Zhou, K., Zhou, L., Lim, Q. 'En, Zou, R., Stephanopoulos, G., and Too, H.-P. (2011). Novel reference genes for quantifying transcriptional responses of *Escherichia coli* to protein overexpression by quantitative PCR. *BMC Mol. Biol.* 12:18. doi: 10.1186/1471-2199-12-18

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

Copyright © 2018 Jousset, Rosinski-Chupin, Takissian, Glaser, Bonnin and Naas. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.