Agnès B. Jousset1,2,3,4Isabelle Rosinski-Chupin4,5Julie Takissian3,4
Philippe Glaser4,5
Rémy A. Bonnin2,3,4*
Thierry Naas1,2,3,4*- 1Department of Bacteriology-Parasitology-Hygiene, Bicêtre Hospital, Assistance Publique-Hôpitaux de Paris, Le Kremlin-Bicêtre, France
- 2Associated French National Reference Center for Antibiotic Resistance, Le Kremlin-Bicêtre, France
- 3EA7361 “Structure, dynamic, function and expression of broad spectrum β-lactamases”, Faculty of Medicine, Paris-Sud University, Le Kremlin-Bicêtre, France
- 4Joint Research Unit Evolution and Ecology of Resistance to Antibiotics, Institut Pasteur-APHP-University Paris Sud, Paris, France
- 5CNRS, UMRS 3525, Paris, France
The diffusion of KPC-2 carbapenemase is closely related to the spread of Klebsiella pneumoniae of the clonal-group 258 and linked to IncFIIK 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 IncFIIK-IncFIB blaKPC−2-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 (blaKPC−2, blaTEM 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.
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 Italy (Bonomo et al., 2017). The blaKPC gene, coding for an Ambler class A carbapenemase is mostly associated with K. pneumoniae but has been reported to a lesser extent in other Enterobacteriacae, 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 blaKPC 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 blaKPC genes have been identified on several transferable plasmids of different incompatibility goups (e.g., IncFIIK, IncA/C, IncF, IncN, IncP, IncR, IncX, and ColE1), but IncFIIK-type plasmids are the most common and are involved in their spreading (Chen et al., 2014b). The pKpQIL plasmid was the first characterized IncFIIK2 blaKPC-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 IncFIIK1-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 IncFIIK 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 blaKPC revealing the versatile expression of this gene (Naas et al., 2012; Cheruvanky et al., 2017; Girlich et al., 2017). Measuring blaKPC RNA expression revealed that blaKPC 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 blaKPC 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 IncFIIK2-FIB-type blaKPC−2 carrying plasmid was performed. To gain more information about the complex regulation of the blaKPC 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 blaKPC−2-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, transmission, and heavy metal (arsenic, copper, and silver) resistance, and presents the replicase repA2 of IncFIB type.
Figure 1. Schematic representation of plasmid pBIC-1a main features. Gray boxes indicate regions with high nucleotide identity (>97%) with pKPN3 (NC_009649.1), pKpQIL (GU595196.1), and pKpQIL-IT (NC_019155.1) plasmids. Genes from the forward strand are represented above the line, and genes from the reverse strand below. Antimicrobial resistance genes are indicated in red, the replicase genes in black, the heavy metal resistance operons in blue, the transfer operon in white, the genes involved in genetic mobile elements in yellow, and other plasmid-located genes in green. Δ indicates a truncated gene. Tnp: transposase. As for pKpQIL-IT, pBIC-1a carries an aminoglycoside resistance gene on a putative composite transposon-like element, IS26-aph(3′)-Ia-ΔIS26, located downstream of Tn4401a. Unlike in pKpQIL-IT, the IS26 located downstream of aph(3′)-Ia is truncated. In addition, an inversion of 3,872 nt fragment is observed as compared to pKpQIL-IT.
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 IncFIIK2 type according to replicon sequence typing (Carattoli et al., 2014), making pBIC1a a multi-replicon IncFIIK2-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 blaKPC−2 gene was found as part of the class II transposon Tn4401a. Additional antimicrobial resistance genes were present: blaTEM−1, blaOXA−9 (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 Tn4401a. 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 IncFIIK2-IncFIB plasmids carrying blaKPC, 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 blaTEM−1 and blaOXA−9, 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 2. Nucleotidic alignment of blaKPC-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), alignement 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 IncFIIK, repA IncFIB, and the regulator repB IncFIIK 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 IncFIIK 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).
Table 1. Expression analysis of pBIC1a genes without imipenem exposure using RPKM values.
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 IncFIIK-type and repA2 IncFIB-type) are highlighted in green. Genes associated to plasmid transfer are boxed in red.
Tansformation 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 blaKPC−2 gene was the most expressed plasmidic gene (Table 1; Figure 3). Two promoters were previously shown to contribute to blaKPC−2 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 blaKPC−2 gene transcription (Figure S3). In addition, RNA-seq analysis revealed no clear transcription termination with reads continuously covering more than one kb downstream blaKPC−2 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 blaTEM−1 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 blaKPC 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 blaKPC-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 blaKPC, blaTEM, 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.
Imipenem-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).
Imipenem-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 FC values < 2. With two oxidation states Cu+ and Cu2+, copper is a cofactor in redox enzymes that uses dioxygen as a substrate (Rademacher and Masepohl, 2012). Bacteria synthetize 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 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).
Table 2. Differentially expressed genes carried by pBIC1a.
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 blaKPC−2 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 blaKPC−2 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).
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.
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…; Table 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 through the electron transport chain, which is dependent on the tricarboxylic acid (TCA) cycle. Hyperactivation of the electron transport chain stimulates superoxide (O2−) formation, which overwhelms superoxide dismutase defenses and leads to the oxidation of iron-sulfur clusters ([4Fe−4S]2+) employed by numerous dehydratase enzymes. The release of ferrous (Fe2+) 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).
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.
A critical parameter to suppress Fenton reaction is the scavenging of H2O2 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 H2O2 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).
The uptaken Fe2+ 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 Fe2+ 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, Fe2+ 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 σ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 susceptiblity 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 IncFIIK-IncFIB blaKPC−2-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 IncFIIK 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 blaKPC−2 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 IncFIIK-IncFIB plasmid carrying blaKPC−2 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 (OD600 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 OD600 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−ΔΔ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).
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.
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 blaKPC−2 carbapenemase gene. (A) Strand-oriented read mapping of Tn4401a using Artemis software (Carver et al., 2008). No clear transcription termination of blaKPC−2 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 Rodríguez-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/gkp148
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.02166-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-Ríos, 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-Martínez, 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 ISKpn7 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.01509-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
Keywords: carbapenemase, KPC-producing plasmids, transcriptome, RNA-seq, oxidative-stress
Citation: Jousset AB, Rosinski-Chupin I, Takissian J, Glaser P, Bonnin RA and Naas T (2018) Transcriptional Landscape of a blaKPC-2 Plasmid and Response to Imipenem Exposure in Escherichia coli TOP10. Front. Microbiol. 9:2929. doi: 10.3389/fmicb.2018.02929
Received: 19 September 2018; Accepted: 14 November 2018;
Published: 03 December 2018.
Edited by:
Raffaele Zarrilli, Department of Public Health, Università degli Studi di Napoli Federico II, ItalyReviewed by:
Vincenzo Di Pilato, Università degli Studi di Firenze, ItalyEliana De Gregorio, Università degli Studi di Napoli Federico II, Italy
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.
*Correspondence: Rémy A. Bonnin, remy.bonnin@u-psud.fr
Thierry Naas, thierry.naas@aphp.fr