Impact of meropenem on Klebsiella pneumoniae metabolism

Claudio Foschi; Melissa Salvo; Luca Laghi; Chenglin Zhu; Simone Ambretti; Antonella Marangoni; Maria Carla Re

doi:10.1371/journal.pone.0207478

Abstract

The aim of this study was to analyze the metabolome of several Klebsiella pneumoniae strains characterized by different resistance patterns. A total of 59 bacterial strains (27 carbapenemase-negative and 32 carbapenemase-positive) were included and their metabolic features were assessed in basal conditions. Moreover, 8 isolates (4 wild-type and 4 KPC-producers) were randomly selected to evaluate the impact of sub-lethal concentrations of meropenem on bacterial metabolism. The metabolomic analysis was performed by ¹H-NMR spectroscopy both on filtered supernatants and cell lysates. A total of 40 and 20 molecules were quantified in the intracellular and the extracellular metabolome, respectively. While in basal conditions only five metabolites showed significant differences between carbapenemase-positive and negative strains, the use of meropenem had a profound impact on the whole bacterial metabolism. In the intracellular compartment, a reduction of different overflow metabolites and organic acids (e.g. formate, acetate, isobutyrate) was noticed, whereas, in the extracellular metabolome, the levels of several organic acids (e.g. succinate, acetate, formate, lactate) and amino acids (aspartate, threonine, lysine, alanine) were modified by meropenem stimulation. Interestingly, carbapenemase-positive and negative strains reacted differently to meropenem in terms of number and type of perturbed metabolites. In wild-type strains, meropenem had great impact on the metabolic pathways related to methane metabolism and alanine, aspartate and glutamate metabolism, whereas in KPC-producers the effect was predominant on pyruvate metabolism. The knowledge about the bacterial metabolic profiles could help to set up innovative diagnostic methods and new antimicrobial strategies to fight the global crisis against carbapenemase-positive K. pneumoniae.

Citation: Foschi C, Salvo M, Laghi L, Zhu C, Ambretti S, Marangoni A, et al. (2018) Impact of meropenem on Klebsiella pneumoniae metabolism. PLoS ONE 13(11): e0207478. https://doi.org/10.1371/journal.pone.0207478

Editor: John J. Kilbane, Illinois Institute of Technology Paul V Galvin Library, UNITED STATES

Received: August 21, 2018; Accepted: October 31, 2018; Published: November 15, 2018

Copyright: © 2018 Foschi et al. This is an open access article distributed under the terms of the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original author and source are credited.

Data Availability: All relevant data are within the manuscript and its Supporting Information files.

Funding: This work was supported by ‘Fondazione del Monte di Bologna e Ravenna’. CZ gratefully acknowledges financial support from Chinese Scholarship Council [grant number 201606910076]. The funders had no role in study design, data collection and analysis, or preparation of the manuscript.

Competing interests: The authors have declared that no competing interests exist.

Introduction

The global spread of carbapenemase-producing Enterobacteriaceae (CPE) is of great concern to health services worldwide [1, 2]. In particular, multi-drug resistant Klebsiella pneumoniae strains, harboring KPC enzymes, have been causing epidemics of international proportions [3, 4]. Healthcare-associated infections caused by CPE represent an alarming and dramatic problem for different reasons. Firstly, most carbapenemase-encoding genes are located on transferable genetic elements that are often associated with other antibiotic resistance genes, thus leading to their rapid transfer and to the spread of multi-drug resistant superbugs [5].

Moreover, the morbidity associated to CPE infections is usually high, with a relevant clinical and economic impact. Indeed, the mortality rate due to CPE infections is usually about 20–30% and can reach 70% in case of bacteremia or pneumonia in critically ill patients [6].

Finally, the therapeutic options for CPE infections are few and limited to old and toxic drugs, thus leading to the onset and spread of new resistance mechanisms (e.g. colistin resistance) [7].

In the last years, several approaches have been proposed to fight the global burden of CPE. On the one hand, screening and surveillance hospital protocols, as well as strict infection control measures (e.g. hand hygiene, patient isolation, personal protection equipment), have been adopted. On the other hand, rapid laboratory assays for the identification of CPE and new antimicrobials acting against KPC-producing Enterobacteriaceae (e.g. ceftazidime/avibactam) have been introduced in the diagnostic and clinical practice, respectively [8–12].

Nevertheless, the global crisis against multi-drug resistant CPE is still ongoing and constitutes a major public health challenge [13, 14]. For this reason, rapid and reliable diagnostic approaches, as well as new antimicrobial drugs, are urgently needed for a better management of CPE infections.

In this context, information about the metabolic profiles of CPE in basal conditions and under ‘antibiotic stress’ could represent an intriguing approach to obtain useful information to set up new diagnostic methods and to develop new antimicrobial strategies.

Proton-based nuclear magnetic resonance (¹H-NMR) spectroscopy proved to be a fast and reliable analytical technique to detect and quantify low molecular weight metabolites in different matrices and biological fluids [15–17]. Moreover, it has been previously used to study the intracellular and extracellular metabolic profiles of different Gram-positive (e.g. lactobacilli and Staphylococcus aureus) and Gram-negative bacteria (e.g. Acinetobacter baumannii, Pseudomonas aeruginosa) [18–21].

It has been shown that data about the bacterial metabolome can be useful to get new insights on the microbial physiology and pathogenesis, to find correlations with the taxonomy and to distinguish bacterial strains on the basis of different growth patterns (e.g. planktonic vs biofilm) or based on various antimicrobial resistance phenotypes (e.g. methicillin-resistant vs methicillin-susceptible Staphylococcus aureus) [18–23].

The aim of this study was to characterize and compare the metabolic profiles of several carbapenemase-positive and carbapenemase-negative K. pneumoniae strains, by means of ¹H-NMR spectroscopy, both in basal conditions and after the exposure to sub-lethal concentrations of meropenem. Major potential applications of this work concern the identification of specific metabolic biomarkers able to differentiate K. pneumoniae strains with different resistance patterns and the detection of critical metabolic pathways, to be used as potential targets of novel antimicrobial drugs.

Materials and methods

Bacterial strains

K. pneumoniae strains were isolated from surveillance and clinical samples submitted to the Microbiology Unit of S. Orsola-Malpighi University Hospital of Bologna (Italy) for routine diagnostic procedures.

The bacterial identification at the species level was obtained by means of a matrix-assisted laser desorption/ionization time-of-flight mass spectrometry (MALDI-TOF MS), using the VITEK MS system (bioMérieux), whereas antimicrobial susceptibility testing (AST) was performed using Vitek AST-N 201 cards on a Vitek 2 automated instrument (bioMérieux).

EUCAST criteria (www.eucast.org) were used for the interpretation of results and for the definition of antimicrobial susceptibility category. In the case of a strain resistant to a third-generation cephalosporins (i.e. cefotaxime and/or ceftazidime), the production of extended-spectrum β-lactamase enzymes was confirmed by double-disk synergy test [24]. Moreover, in the case of a minimum inhibitory concentration (MIC) of meropenem > 0.25 mg/L, a combination disk testing (disks of meropenem ± various inhibitors) was performed, in order to confirm the positivity for carbapenemase and to define the type of enzyme produced [25]. Temocillin susceptibility testing was used for the detection of OXA-48-like enzymes [26].

Starting from overnight grown bacterial colonies, the presence of a carbapenemase was further investigated by a genotypic method, based on the real-time amplification of the major carbapenemase genes (bla_KPC, bla_IMP, bla_VIM, bla_NDM and bla_OXA-48), using a commercially available assay (Xpert Carba-R, Cepheid).

For all the isolates, the exact MIC values of meropenem were calculated by means of a broth micro-dilution assay, following CLSI guidelines [27].

Sample preparation

For each strain, bacteria were grown overnight in Mueller-Hinton (MH; Sigma-Aldrich) agar at 37°C in aerobic atmosphere. Afterwards, 400 μl of a 2 McFarland (McF) turbidity suspension (corresponding to 0.5 optical density [OD] at 550 nm) were inoculated into 5 mL of MH broth and incubated under shaking for 6–8 hours at 37°C. Each bacterial suspension was then adjusted to the same turbidity value (2.8 McF, corresponding approximately to 8 × 10⁸ colony forming units/mL and 0.7 OD at 550 nm) by adding MH broth. One mL of the bacterial suspension was centrifuged at 5,000 × g for 20 minutes at 4°C in order to separate the cell pellet from the supernatant. Then, supernatants were filtered through a 0.2 μm membrane filter to obtain cell free supernatants and stored at -80°C.

Cell pellets were then washed twice in sterile saline and internal metabolites were extracted with cold methanol, following the protocol described by Park et al., with slight modifications [28]. Briefly, the biomass pellets were suspended in 250 μL of 100% cold methanol (−20°C) and 250 μL of sterile saline, thoroughly vortexed and frozen at −80°C for 30 minutes. Afterwards, the cell pellets were allowed to thaw on ice for 30 minutes. The freeze–thaw cycles were repeated three times to cause the metabolites to leak from the cells. The suspensions were then centrifuged at 13,000 × g for 10 min. The supernatants were collected and stored on dry ice, whereas 500-μL aliquots of 100% cold methanol (250 μL) and sterile saline (250 μL) were added to the pellets. The entire above procedure was repeated, then the two aliquots were combined. Finally, the solutions were evaporated under reduced pressure and the samples were stored at −80°C.

At the time of metabolomic analysis, all the samples were gently thawed. Both the filtered supernatants and the cell lysates were analyzed by ¹H-NMR to examine the extracellular and intracellular metabolome, respectively, as described below.

The same protocol was used to investigate the metabolic profile of 8 randomly selected K. pneumoniae strains (4 wild-type isolates and 4 KPC-producers) under ‘antibiotic-stress’. In particular, bacteria were allowed to grow for 6–8 hours in antibiotic-free MH broth and, in parallel, in MH broth added from the beginning of the incubation with sub-lethal concentrations (1/8 of the MIC value) of meropenem (Sigma-Aldrich).

As preliminary investigation, an evaluation of the bacterial growth curves in absence and in presence of various sub-MIC levels of meropenem (1/2, 1/4 and 1/8 of the MIC values) was performed by the measurement of the optical density values at different time points, at 550 nm. For a value of 1/8 of the MIC, the bacterial growth in presence of meropenem was comparable to the cultures without antibiotic (S1 Fig). Therefore, 1/8 of the MIC was chosen for antibiotic exposure assays, as the highest sub-lethal dose of meropenem that did not affect the bacterial growth.

After a sample size calculation based on previously published data [29], the 8 strains used for antibiotic-exposure experiments were selected randomly, considering that the metabolic profiles of the 59 K. pneumoniae isolates in basal conditions were very similar to each other (see Results section). To assess whether the randomly selected strains were representative of the corresponding groups, a rPCA model was calculated on all the strains and all the molecules of the extracellular compartment (S2 Fig). A screen test allowed to observe that three components offered the best compromise between variance explained and model complexity [30]. Along each of the three components, representing as much as the 74.8% of the total variance, the samples selected never exceeded 1.5 times the interquartile range of the corresponding group and could be therefore considered representative of the corresponding groups [31].

Metabolomic analysis

The metabolomic investigation was performed by means of ¹H-NMR spectroscopy, starting from the filtered supernatants (‘external metabolome’) and the cell lysates (‘internal metabolome’).

For each K. pneumoniae strain, 700 μl of cell-free supernatant and 100 μl of cellular lysate were added to 100 μl of a D₂O solution of 3-(trimethylsilyl)-propionic-2,2,3,3-d4 acid sodium salt (TSP) 10 mM set to pH 7.0 ± 0.02 by means of a 1M phosphate buffer, containing also NaN₃ 2 mM to avoid bacteria proliferation. Cellular lysate samples were also added with 600 μl of deionized water, to meet the sample volume specifications of the NMR probe.

¹H-NMR spectra were recorded at 298 K with an AVANCE III spectrometer (Bruker, Milan, Italy) operating at a frequency of 600.13 MHz, equipped with Topspin software (Ver. 3.5) [17, 32]. The spectra were formed by 32K points, spanning a spectral width of 11.97 ppm, registered for 256 scans. Each spectrum was then processed (alignment, baseline adjustment) following the procedure previously described [17, 33, 34].

The signals were assigned by comparing their chemical shift and multiplicity with the Human Metabolome Database [35] and the compounds library (Ver. 10) of Chenomx software (Chenomx Inc., Canada, Ver. 8.3). Quantification of the molecules was performed in the first sample acquired by employing the added TSP as an internal standard. In the case of cell-free supernatant, the first sample was represented by the growth medium. In order to compensate for differences in solids content, any other sample was then normalized to such sample by means of probabilistic quotient normalization [36]. Integration of the signals was performed for each molecule by means of rectangular integration.

Metabolic pathway analysis

MetaboAnalyst 3.0 (www.metaboanalyst.ca), a web server designed to permit comprehensive metabolomic data analysis, was used to evaluate the impact of sub-lethal meropenem concentration on the bacterial metabolic pathways [37].

The program is based on several databases, including KEGG (http://www.genome.jp/kegg/) and the Human Metabolome Database (http://www.hmdb.ca/) and identifies the pathways that are mostly perturbed.

The final list of altered metabolites that were significantly influenced by drug treatment, was analyzed by Metabolomics Pathway Analysis (MetPA) software in MetaboAnalyst. The built-in prokaryote (Escherichia coli K-12 MG1655) pathway library was used as reference. The following pathways analysis algorithms were employed: the hypergeometric test for over-representation analysis and the relative-betweenness centrality for pathway topology analysis [38, 39].

The output of this program marks a metabolic pathway as significant if a significantly higher number of compounds involved in that pathway is present in the input list than would be expected by random chance. In the graphical outputs (metabolome view), each dot represents one metabolic pathway. The circle size and color shade of each circle are based on pathway impact value and P-value (red being the most significant), respectively [38, 39].

Statistical analysis

All statistical analyses were performed by using R computational language (www.r-project.org) and GraphPad Prism version 5.02 for Windows (GraphPad software).

Differences in the extracellular/intracellular metabolome composition between bacterial groups were assessed by means of Kruskal-Wallis test, two-tailed paired or unpaired Wilcoxon test or t-test, where appropriate (univariate analysis).

A robust principal component analysis (rPCA) (multivariate analysis) [40] was also performed to obtain an overview of the differences between non-stimulated and meropenem-stressed bacterial groups, on the basis of the significant metabolites found by means of the univariate analysis. The scoreplot—the representation of the samples from the new point of view—helps to visually inspect the trends underlying the samples, while the loadingplot (barplot) evidences the molecules mainly driving the trends. Statistical significance was determined at P < 0.05.

Results

Bacterial strains

A total of 59 K. pneumoniae strains were included in the study. On the basis of AST and phenotypic test results, K. pneumoniae isolates were categorized in carbapenemase-positive (n = 32) and carbapenemase-negative strains (n = 27). In particular, carbapenemase-positive strains included 24 KPC-producers, 4 strains positive for metallo-β-lactamases (MBL; specifically 3 NDM and 1 VIM) and 4 isolates with OXA-48 enzymes. The class of carbapenemase detected by phenotypic tests (combination disk testing) always agreed with the genotypic detection (real-time PCR) of carbapenemase genes.

The group of carbapenemase-negative strains included 16 wild-type K. pneumoniae strains and 11 ESBL-producers.

No significant differences were noticed with respect to bacterial growth between the different groups of strains.

The full list of bacterial strains with the source of isolation, the type of enzyme detected and the exact MIC values of meropenem are detailed in S1 Table.

Metabolomic analysis

Consistent with previous reports on similar matrices [18, 41], a total of 40 and 20 molecules were detected and quantified by ¹H-NMR analysis in the extracellular and intracellular metabolome, respectively. These metabolites mainly belong to the families of amino acids, organic acids and alcohols (S2 Table).

In order to highlight significant differences in the metabolic profiles of K. pneumoniae strains, different sets of analyses have been carried out, as follows: (i) in basal conditions, a comparison between carbapenemase-positive (n = 32) and carbapenemase-negative strains (n = 27), as well as between the three groups of strains harboring different types of carbapenemase (24 KPC, 4 MBL and 4 OXA-48), and between wild-type (n = 16) and ESBL-producing strains (n = 11); (ii) an analysis of a group of 8 randomly selected strains in antibiotic-free conditions and after meropenem exposure (1/8 of the MIC values), irrespective of the carbapenemase production; (iii) a comparison between non-stimulated and antibiotic-stressed bacteria, stratified by the carbapenemase production. In particular, for theses latter analyses, 4 KPC-producers (Kp08, Kp12, Kp15, Kp53) and 4 wild-type strains (Kp17, Kp25, Kp26, Kp61) were considered.

The results of these three different sets of analyses are detailed in the following paragraphs.

(i) Metabolome of K. pneumoniae in basal conditions

Globally, the intracellular and extracellular metabolic profiles of the 59 K. pneumoniae strains in basal conditions were rather similar. Indeed, in antibiotic-free conditions, we found 2 metabolites in the intracellular metabolome (betaine and isobutyrate) and 3 molecules in the external metabolome (phenylalanine, sarcosine, valine), respectively, whose concentration differed significantly between carbapenemase-positive and carbapenemase-negative K. pneumoniae strains (Table 1). Similarly, in basal conditions, only 4 molecules showed significantly different concentrations in the comparison between KPC, MBL and OXA-48 positive strains (Table 2). Finally, two metabolites of the intracellular compartment (glutamate and acetate) were found to be significantly different when the metabolic profiles of wild-type and ESBL-positive strains were compared.

Download:

Table 1. List of metabolites whose concentration (mmol/L; mean ± SD) differed significantly between carbapenemase-positive (+) and carbapenemase-negative (-) K. pneumoniae strains in basal conditions.

https://doi.org/10.1371/journal.pone.0207478.t001

Download:

Table 2. List of metabolites whose concentration (mmol/L; mean ± SD) differed significantly between KPC, MBL and OXA-48 K. pneumoniae strains in antibiotic-free conditions.

https://doi.org/10.1371/journal.pone.0207478.t002

(ii) Metabolome of K. pneumoniae strains under ‘meropenem stress’

A total of 9 and 22 metabolites in the internal and external metabolome (S3 and S4 Tables), respectively, showed significant variations in antibiotic-stressed bacteria compared to non-stimulated strains, thus indicating that meropenem had globally a profound impact on the bacterial metabolism.

In the intracellular compartment, the most significant changes regarded the reduction of different overflow metabolites and organic acids (e.g. formate, acetate, isobutyrate), with the exception of succinate, showing higher levels after meropenem exposure.

When the extracellular metabolome was analyzed, it was possible to highlight that the levels of several metabolites, mostly organic acids and amino acids, were modified by the meropenem stimulation. As examples, aspartate, threonine, lysine and 2-hydroxyisobutyrate were secreted in significantly larger amounts, whereas a lower production of urocanate, succinate, acetate, ethanol, propionate and alanine was detected. Moreover, the antibiotic stimulation led to a higher consumption of formate, lactate and maltose by the bacterial cells.

In order to investigate underlying trends in the metabolome subspace represented by these molecules, their concentrations were employed as the base for a rPCA model. In the intracellular metabolome, PC1, explaining the 91.9% of the total variance, effectively accounted for the stress connected to the presence of a meropenem, with non-stimulated strains appearing at low PC1 scores and clearly separated from antibiotic-stressed K. pneumoniae isolates (P<0.001) (Fig 1). The correlation between the concentration of each molecule and its importance over PC1 served as a handy way to rank the molecules in order of importance in determining such grouping. As pointed out in Fig 1, higher concentrations of succinate and lower concentrations of isocaproate and glycerol were the most significant changes after the stimulation with sub-lethal concentration of meropenem in the internal metabolome.

Download:

Fig 1.

rPCA model calculated on the space constituted by the concentration of the intracellular (upper panel) and extracellular (lower panel) molecules that significantly differed between non-stressed and meropenem-stressed bacterial strains. In the scoreplot (left), non-stressed strains and meropenem-stressed strains are represented with black square and red dots respectively, with lines connecting each strain to the median of its group. In the barplot (right), describing the correlation between the concentration of each molecule and its importance over PC1, dark gray bars highlight statistically significant correlations (P<0.05). In the barplot, bars pointing to the left denote metabolites that were more abundant in non-stressed cells, while bars pointing to the right denote molecules that were more abundant in meropenem-stressed bacteria.

https://doi.org/10.1371/journal.pone.0207478.g001

Similarly, in the external metabolome, non-stimulated strains were clearly separated from meropenem-stressed isolates along PC1, explaining the 92.1% of the total variance (Fig 1).

Higher concentrations of hydroxyisobutirate, threonine and aspartate and lower concentrations of succinate, acetate and citraconate were the main characteristics of the stimulated K. pneumoniae strains in contrast to non-stimulated bacteria (Fig 1).

Considering globally the data about the bacterial cell lysates and the supernatants, the analysis with MetaboAnalyst allowed to identify the specific metabolic pathways most affected by drug exposure.

As shown in Fig 2, the impact of meropenem on K. pneumoniae was dramatic: six metabolic pathways showed significant P-value following exposure (top portion with y-axis higher or equal to 5: glycine, serine and threonine metabolism; aminoacyl-tRNA biosynthesis; cyanoamino acid metabolism; methane metabolism; valine, leucine and isoleucine biosynthesis; pyruvate metabolism), and four metabolic pathways showed great pathway impact (x-axis larger or equal to 0.4: glycine, serine and threonine metabolism; methane metabolism; pyruvate metabolism; alanine, aspartate and glutamate metabolism).

Download:

Fig 2. Metabolome view of the metabolic pathways altered after meropenem exposure in the group of the 8 selected K. pneumoniae strains.

Each dot represents a unique metabolic pathway, with the dot size corresponding to the pathway impact score and the dot color (red being the most significant) corresponding to the–log(P) value.

https://doi.org/10.1371/journal.pone.0207478.g002

(iii) Analysis of the metabolic response to meropenem of wild-type and KPC-producing strains

As shown in Tables 3 and 4, the presence of sub-lethal concentration of meropenem led to different metabolic variations in the comparison between carbapenemase-negative (4 wild-type) and carbapenemase-positive strains (4 KPC-producers), in terms of both number and type of perturbed metabolites.

Indeed, in the group of wild type isolates, a total of 21 molecules was altered by the drug exposure, whereas 26 metabolites showed significant changes in the group of KPC-producers after meropenem stimulation. Among those molecules, 14 were altered in the two groups in the same direction (increase/reduction), while several others were perturbed only in one of the two groups (e.g. acetate in the intracellular metabolome of wild-type strains; urocanate and cadaverine in the external metabolome of KPC-positive isolates).

Considering the altered metabolites, an analysis of the most affected metabolic pathways was performed. In wild-type K. pneumoniae (Fig 3), three metabolic pathways were significantly altered after the drug exposure (alanine, aspartate and glutamate metabolism; aminoacyl-tRNA biosynthesis; pyruvate metabolism), and two showed great pathway impact (methane metabolism; alanine, aspartate and glutamate metabolism).

Download:

Fig 3.

Metabolome view of the metabolic pathways altered after meropenem exposure in the group of wild-type (left) and KPC-producing (right) K. pneumoniae strains. Each dot represents a unique metabolic pathway, with the dot size corresponding to the pathway impact score and the dot color (red being the most significant) corresponding to the–log(P) value.

https://doi.org/10.1371/journal.pone.0207478.g003

Download:

Table 3. List of intracellular/extracellular metabolites that significantly differed after the exposure with sub-lethal concentrations of meropenem in the group of carbapenemase-negative (wild-type) K. pneumoniae isolates.

The concentration of metabolites is expressed as mmol/L (mean ± SD).

https://doi.org/10.1371/journal.pone.0207478.t003

Download:

Table 4. List of intracellular/extracellular metabolites that significantly differed after the exposure with sub-lethal concentrations of meropenem in the group of carbapenemase-positive (KPC) K. pneumoniae isolates.

The concentration of metabolites is expressed as mmol/L (mean ± SD).

https://doi.org/10.1371/journal.pone.0207478.t004

Conversely, in KPC-producing strains (Fig 3), five metabolic pathways showed significant P-values following meropenem stimulation (glycine, serine and threonine metabolism; alanine, aspartate and glutamate metabolism; butanoate metabolism; C5-branched dibasic acid metabolism; citrate cycle or tricarboxylic acid cycle) and one showed great pathway impact (pyruvate metabolism).

Discussion

To the best of our knowledge, this is the first report about the analysis of the metabolic features of several K. pneumoniae strains characterized by different resistance patterns before and after the exposure to meropenem. Previous works have investigated the volatile metabolome of K. pneumoniae in human blood or during in vitro growth, but none has examined the effect of meropenem on these bacteria, divided according to carbapenemase production [42–44].

At first, we looked for differences in the metabolic profiles between carbapenemase-positive and carbapenemase-negative K. pneumoniae strains in the absence of a meropenem exposure, as well as between the strains harboring different types of carbapenemase. In basal conditions, only few metabolites showed significant differences between these groups of isolates. This result agrees with a previous report about S. aureus, demonstrating that without exposure to antibiotics, the concentrations of only 5 metabolites varied significantly in comparison of isogenic pairs of methicillin-resistant (MRSA) and methicillin-susceptible (MSSA) strains [29].

Afterwards, we evaluated the global effect of sub-lethal concentrations of meropenem on the bacterial metabolome. To this purpose, we randomly selected 4 wild-type K. pneumoniae isolates and 4 strains positive for KPC enzymes, representing the commonest types of carbapenemase in many parts of the world [45].

Irrespective of the resistance pattern, the presence of meropenem induced profound perturbations in the bacterial metabolome. Indeed, the exposure to meropenem led to significant variations in the levels of different amino acids and organic acids with a great impact on several metabolic pathways including pyruvate metabolism and alanine, aspartate and glutamate metabolism.

It has been shown that β-lactams, besides their primary mechanism of action (i.e. inhibition of cell-wall synthesis), cause subsequent metabolic changes that occur downstream of the interaction of the antibiotic with their target, playing an important role in bacterial lethality [46, 47].

Effectively, β-lactams display a strong perturbations on bacterial metabolic profiles with the increase or decrease in the concentrations of several metabolites, suggesting that antibiotic treatment have broad and complex effects on metabolism and do not simply quench all metabolic activities [48].

In agreement with our findings, it has been previously shown, that the exposure to lethal or sub-lethal concentrations of β-lactams induces significant changes in metabolites involved in several pathways, like glycolysis, acetate metabolism, tricarboxylic acid (TCA) cycle, NAD and amino acids metabolism [48–50].

For example, in S. aureus, lethal doses of ampicillin induce a higher consumption of glucose, glutamine, valine, leucine and tyrosine, as well as a higher secretion of aspartate, glutamate, isobutyrate, acetate and pyruvate [49].

Similarly, in Escherichia coli, ampicillin leads to an increase of amino acids metabolites at 60 and 90 minutes post-treatment, as well as to higher concentrations of metabolites involved in central carbon metabolism, as citrate and succinate [48].

Comparing these data with our results, we noticed that some modifications showed similar trends (e.g. higher levels of pyruvate and aspartate), whereas others were opposed to those described for E. coli and S. aureus (e.g. variations of succinate and glutamate).

It should be considered that the bacterial metabolic response to an antibiotic is not only specific for each type of microorganism, but also depends on many other factors, as the type and dose (lethal/sub-lethal concentrations) of antimicrobial and the time of the analysis following the drug treatment.

Probably, the most interesting data concerned the assessment of the bacterial response to meropenem, stratified by the capability of carbapenemase production (wild-type strains vs KPC-producers).

It has been reported that KPC-producing K. pneumoniae strains differs from carbapenemase-negative due to the presence of the bla_KPC gene, carried on large plasmids and typically present in a Tn3-based transposon. This gene encodes for the KPC enzyme, responsible for high levels of penicillin and cephalosporin resistance and variable carbapenem resistance and it is usually accompanied by other drug resistance determinants (e.g. aminoglycosides, quinolones, trimethoprim and tetracyclines) [44, 51].

Hence, we hypothesized that even the exposure to sub-lethal levels of meropenem would cause different modifications in the metabolome between carbapenemase-negative and positive isolates.

Indeed, the bacterial response to meropenem differed significantly between wild-type and KPC-positive strains, in terms of type of altered metabolites and impaired metabolic pathways.

It has been previously shown that, after the exposure to sub-lethal doses of methicillin, MRSA and MSSA strains have different metabolic reactions, regarding the type, the number and the trend of change of metabolites, with great impact on different metabolic pathways. Globally, the analysis suggested that the metabolic activities of MSSA strains were more susceptible to the perturbation of methicillin exposure compared to the MRSA strains [29, 52].

Conversely, in this work, we noticed that carbapenemase-positive K. pneumoniae strains were not more resistant to the metabolic changes caused by meropenem compared to wild-type isolates. Rather, even though in presence of common metabolic alterations, some variations were specific for KPC-producers and others acted as biomarkers for wild-type isolates.

We are fully aware that some limitations could have affected our results. Among all, the relatively low sensitivity of ¹H-NMR spectroscopy could have led to miss significant molecules in the metabolome of K. pneumoniae strains with a partial representation of the metabolic pathways perturbed by the meropenem stimulation. As an order of magnitude, ¹H-NMR has been shown to detect metabolites with concentrations of at least 1 to 5 μM, higher than other metabolomics platforms (e.g. mass spectrometry) [53–55]. On the other hand, it is worth emphasizing some important strengths of the metabolomic analysis by ¹H-NMR: (i) it is an intrinsically quantitative technique, (ii) the experimental protocol requires a very limited sample manipulation, allowing to process a large number of samples simultaneously [15, 56].

Conclusions

In conclusion, the metabolomic analysis of the bacterial isolates by ¹H-NMR allowed to detect and quantify a list of molecules that significantly differed between carbapenemase-positive and carbapenemase-negative strains, both in antibiotic-free conditions and under ‘meropenem stress’.

The data obtained could serve to set up innovative diagnostic strategies based on metabolic biomarkers, as well as to elucidate the action of β-lactams on the bacterial metabolome. The knowledge about these cellular responses are essential for improving existing therapies and identifying novel therapeutic approaches [46, 47]. In this way, it could be possible to open up significant perspectives for the control of the spread of CPE, still representing a worrisome emergency for public health structures and healthcare systems worldwide [57, 58].

Future perspectives will include the analysis of a larger panel of K. pneumoniae strains, evaluating the effect of different types of antimicrobials and various combinations of drugs on the bacterial metabolome. Moreover, the study of the effect of lethal concentrations of beta-lactams will allow the identification of the earliest and most abundant metabolites associated with the bacterial death. In this way, innovative approaches for the evaluation of the effectiveness of several antimicrobials against carbapenemase-positive K. pneumoniae could be possible.

Finally, it will be crucial to evaluate the metabolic profiles of K. pneumoniae strains getting deeper insights into their genetic background (e.g. sequence types by MLST or distinction between different classes of KPC enzymes), in order to understand if the observed data are limited to specific clones or, if they are common to bacterial isolates with different genetic origins.

Supporting information

S1 Fig. Bacterial growth curves at different sub-lethal levels of meropenem MIC.

Panel A: carbapenemase-positive (KPC-producers) strains; Panel B: carbapenemase-negative (wt) strains. The growth was determined by the measurement of the optical density at 550 nm (OD₅₅₀). For each time point, the mean value ± SEM (standard error of the mean) of the OD related to the strains selected for meropenem experiments is shown.

https://doi.org/10.1371/journal.pone.0207478.s001

(TIF)

S2 Fig.

(A) Scoreplot of an rPCA model calculated on the space constituted by the concentration of the 40 molecules identified in the extracellular metabolome. Black (Carbapenemase-negative; C-) and red (Carbapenemase-positive; C+) lines connect each strain to the median of its group, while circles evidence the 8 samples considered to investigate the effect of meropenem. (B-D) The position of the samples along PC1, 2 and 3 is summarized as boxplots. The dashed lines evidence the inter-quartile distance multiplied by 1.5.

https://doi.org/10.1371/journal.pone.0207478.s002

(TIFF)

S1 Table. List of Klebsiella pneumoniae strains included in the study.

wt: wild-type (no acquired resistance mechanisms); ESBL: extended-spectrum β-lactamase; MBL: metallo-β-lactamase; NDM: New Delhi metallo-β-lactamase; VIM: Verona Integron-encoded metallo-β-lactamase.

https://doi.org/10.1371/journal.pone.0207478.s003

(DOCX)

S2 Table. List of the molecules found in the intracellular and extracellular metabolome of Klebsiella pneumoniae strains.

https://doi.org/10.1371/journal.pone.0207478.s004

(DOCX)

S3 Table. List of intracellular molecules that significantly differed after meropenem exposure in the group of selected K. pneumoniae isolates, irrespective of the carbapenemase production.

The concentration of metabolites is expressed as mmol/L (mean ± SD).

https://doi.org/10.1371/journal.pone.0207478.s005

(DOCX)

S4 Table. List of extracellular molecules that significantly differed after meropenem exposure in the group of selected K. pneumoniae isolates, irrespective of their carbapenemase production.

The concentration of metabolites is expressed as mmol/L (mean ± SD).

https://doi.org/10.1371/journal.pone.0207478.s006

(DOCX)

Acknowledgments

We would like to thank Dr. Maria Vittoria Tamburini for the excellent technical assistance during the selection of bacterial strains.

References

1. Tzouvelekis LS, Markogiannakis A, Psichogiou M, Tassios PT, Daikos GL. Carbapenemases in Klebsiella pneumoniae and other Enterobacteriaceae: an evolving crisis of global dimensions. Clin Microbiol Rev. 2012;25: 682–707. pmid:23034326
- View Article
- PubMed/NCBI
- Google Scholar
2. European Centre for Disease Prevention and Control. Antimicrobial resistance surveillance in Europe 2014. Annual Report of the European Antimicrobial Resistance Surveillance Network (EARS-Net). Stockholm: ECDC; 2015.
3. Giani T, Pini B, Arena F, Conte V, Bracco S, Migliavacca R, et al. Epidemic diffusion of KPC carbapenemase-producing Klebsiella pneumoniae in Italy: results of the first countrywide survey, 15 May to 30 June 2011. Euro Surveill. 2013; 18: pii = 20489.
- View Article
- Google Scholar
4. Grundmann H, Livermore DM, Giske CG, Canton R, Rossolini GM, Campos J, et al. Carbapenem-non-susceptible Enterobacteriaceae in Europe: conclusions from a meeting of national experts. Euro Surveill. 2010;15: pii = 19711.
- View Article
- Google Scholar
5. Diene SM, Rolain JM. Carbapenemase genes and genetic platforms in Gram-negative bacilli: Enterobacteriaceae, Pseudomonas and Acinetobacter species. Clin Microbiol Infect. 2014;20: 831–838. pmid:24766097
- View Article
- PubMed/NCBI
- Google Scholar
6. Daikos GL, Tsaousi S, Tzouvelekis LS, Anyfantis I, Psichogiou M, Argyropoulou A, et al. Carbapenemase-producing Klebsiella pneumoniae bloodstream infections: lowering mortality by antibiotic combination schemes and the role of carbapenems. Antimicrob Agents Chemother. 2014;58: 2322–2328. pmid:24514083
- View Article
- PubMed/NCBI
- Google Scholar
7. Monaco M, Giani T, Raffone M, Arena F, Garcia-Fernandez A, Pollini S, et al. Colistin resistance superimposed to endemic carbapenem-resistant Klebsiella pneumoniae: a rapidly evolving problem in Italy, November 2013 to April 2014. Euro Surveill. 2014;19: pii = 20939.
- View Article
- Google Scholar
8. Akova M, Daikos GL, Tzouvelekis L, Carmeli Y. Interventional strategies and current clinical experience with carbapenemase-producing Gram-negative bacteria. Clin Microbiol Infect. 2012;18: 439–448. pmid:22507111
- View Article
- PubMed/NCBI
- Google Scholar
9. De Rosa FG, Corcione S, Cavallo R, Di Perri G, Bassetti M. Critical issues for Klebsiella pneumoniae KPC-carbapenemase producing K. pneumoniae infections: a critical agenda. Future Microbiol. 2015;10: 283–294. pmid:25689539
- View Article
- PubMed/NCBI
- Google Scholar
10. Findlay J, Hopkins KL, Meunier D, Woodford N. Evaluation of three commercial assays for rapid detection of genes encoding clinically relevant carbapenemases in cultured bacteria. J Antimicrob Chemother. 2015;70: 1338–1342. pmid:25630646
- View Article
- PubMed/NCBI
- Google Scholar
11. Dortet L, Agathine A, Naas T, Cuzon G, Poirel L, Nordmann P. Evaluation of the RAPIDEC® CARBA NP, the Rapid CARB Screen® and the Carba NP test for biochemical detection of carbapenemase-producing Enterobacteriaceae. J Antimicrob Chemother. 2015;70: 3014–3022. pmid:26260131
- View Article
- PubMed/NCBI
- Google Scholar
12. Sharma R, Eun Park T, Moy S. Ceftazidime-Avibactam: A Novel Cephalosporin/β-Lactamase Inhibitor Combination for the Treatment of Resistant Gram-Negative Organisms. Clin Ther. 2016;pii: S0149-2918(16)30041-8.
- View Article
- Google Scholar
13. Martens E, Demain AL. The antibiotic resistance crisis, with a focus on the United States. J Antibiot (Tokyo). 2017;70: 520–526.
- View Article
- Google Scholar
14. Tzouvelekis LS, Markogiannakis A, Psichogiou M, Tassios PT, Daikos GL. Carbapenemases in Klebsiella pneumoniae and other Enterobacteriaceae: an evolving crisis of global dimensions. Clin Microbiol Rev. 2012;25: 682–707. pmid:23034326
- View Article
- PubMed/NCBI
- Google Scholar
15. Laghi L, Picone G, Capozzi F. Nuclear magnetic resonance for foodomics beyond food analysis. TrAC Trends in Analytical Chemistry. 2014;59: 93–102.
- View Article
- Google Scholar
16. Parolin C, Foschi C, Laghi L, Zhu C, Banzola N, Gaspari V, et al. Insights Into Vaginal Bacterial Communities and Metabolic Profiles of Chlamydia trachomatis Infection: Positioning Between Eubiosis and Dysbiosis. Front Microbiol. 2018;9: 600. pmid:29643849
- View Article
- PubMed/NCBI
- Google Scholar
17. Foschi C, Laghi L, D'Antuono A, Gaspari V, Zhu C, Dellarosa N, et al. Urine metabolome in women with Chlamydia trachomatis infection. PLoS One. 2018;13: e0194827. pmid:29566085
- View Article
- PubMed/NCBI
- Google Scholar
18. Foschi C, Laghi L, Parolin C, Giordani B, Compri M, Cevenini R, et al. Novel approaches for the taxonomic and metabolic characterization of lactobacilli: Integration of 16S rRNA gene sequencing with MALDI-TOF MS and 1H-NMR. PLoS One. 2017;12: e0172483. pmid:28207855
- View Article
- PubMed/NCBI
- Google Scholar
19. Ammons MC, Tripet BP, Carlson RP, Kirker KR, Gross MA, Stanisich JJ, et al. Quantitative NMR metabolite profiling of methicillin-resistant and methicillin-susceptible Staphylococcus aureus discriminates between biofilm and planktonic phenotypes. J Proteome Res. 2014;13: 2973–2985. pmid:24809402
- View Article
- PubMed/NCBI
- Google Scholar
20. Yeom J, Shin JH, Yang JY, Kim J, Hwang GS. (1)H NMR-based metabolite profiling of planktonic and biofilm cells in Acinetobacter baumannii 1656–2. PLoS One. 2013;8: e57730. pmid:23483923
- View Article
- PubMed/NCBI
- Google Scholar
21. Chen T, Sheng J, Fu Y, Li M, Wang J, Jia AQ. 1H NMR-Based Global Metabolic Studies of Pseudomonas aeruginosa upon Exposure of the Quorum Sensing Inhibitor Resveratrol. J Proteome Res. 2017;16: 824–830. pmid:28094526
- View Article
- PubMed/NCBI
- Google Scholar
22. Antti H, Fahlgren A, Näsström E, Kouremenos K, Sundén-Cullberg J, Guo Y, et al. Metabolic profiling for detection of Staphylococcus aureus infection and antibiotic resistance. PLoS One. 2013;8: e56971. pmid:23451124
- View Article
- PubMed/NCBI
- Google Scholar
23. Liebeke M, Dörries K, Zühlke D, Bernhardt J, Fuchs S, Pané-Farré J, et al. A metabolomics and proteomics study of the adaptation of Staphylococcus aureus to glucose starvation. Mol Biosyst. 2011;7: 1241–1253. pmid:21327190
- View Article
- PubMed/NCBI
- Google Scholar
24. Paterson DL, Bonomo RA. Extended-spectrum beta-lactamases: a clinical update. Clin Microbiol Rev. 2005;18: 657–686. pmid:16223952
- View Article
- PubMed/NCBI
- Google Scholar
25. Ambretti S, Gaibani P, Berlingeri A, Cordovana M, Tamburini MV, Bua G, et al. Evaluation of phenotypic and genotypic approaches for the detection of class A and class B carbapenemases in Enterobacteriaceae. Microb Drug Resist. 2013;19: 212–215. pmid:23336577
- View Article
- PubMed/NCBI
- Google Scholar
26. van Dijk K, Voets GM, Scharringa J, Voskuil S, Fluit AC, Rottier WC, et al. A disc diffusion assay for detection of class A, B and OXA-48 carbapenemases in Enterobacteriaceae using phenyl boronic acid, dipicolinic acid and temocillin. Clin Microbiol Infect. 2014;20: 345–349. pmid:23927659
- View Article
- PubMed/NCBI
- Google Scholar
27. CLSI. Performance Standards for Antimicrobial Susceptibility Testing. 27th ed. CLSI supplement M100. Wayne, PA: Clinical and Laboratory Standards Institute; 2017.
28. Park C, Yun S, Lee SY, Park K, Lee J. Metabolic profiling of Klebsiella oxytoca: evaluation of methods for extraction of intracellular metabolites using UPLC/Q-TOF-MS. Appl Biochem Biotechnol. 2012;167: 425–438. pmid:22555499
- View Article
- PubMed/NCBI
- Google Scholar
29. Schelli K, Rutowski J, Roubidoux J, Zhu J. Staphylococcus aureus methicillin resistance detected by HPLC-MS/MS targeted metabolic profiling. J Chromatogr B Analyt Technol Biomed Life Sci. 2017;1047: 124–130. pmid:27316783
- View Article
- PubMed/NCBI
- Google Scholar
30. Ruscio J, Roche B. Determining the number of factors to retain in an exploratory factor analysis using comparison data of known factorial structure. Psychol Assess. 2012;24: 282–292. pmid:21966933
- View Article
- PubMed/NCBI
- Google Scholar
31. Ghasemi A, Zahediasl S. Normality tests for statistical analysis: a guide for non-statisticians. Int J Endocrinol Metab. 2012;10: 486–489. pmid:23843808
- View Article
- PubMed/NCBI
- Google Scholar
32. Parolin C, Marangoni A, Laghi L, Foschi C, Ñahui Palomino RA, Calonghi N, et al. Isolation of Vaginal Lactobacilli and Characterization of Anti-Candida Activity. PLoS One. 2015;10: e0131220. pmid:26098675
- View Article
- PubMed/NCBI
- Google Scholar
33. Kneen MA, Annegarn HJ. Algorithm for fitting XRF, SEM and PIXE X-ray spectra backgrounds. Nucl Instrum Methods Phys Res B. 1996;109: 209–213.
- View Article
- Google Scholar
34. Liland KH, Almøy T, Mevik BH. Optimal choice of baseline correction for multivariate calibration of spectra. Appl Spectrosc. 2010;64: 1007–1016. pmid:20828437
- View Article
- PubMed/NCBI
- Google Scholar
35. Wishart DS, Tzur D, Knox C, Eisner R, Guo AC, Young N, et al. HMDB: the human metabolome database. Nucleic Acids Res. 2007;35: D521–D526. pmid:17202168
- View Article
- PubMed/NCBI
- Google Scholar
36. Dieterle F, Ross A, Schlotterbeck G, Senn H. Probabilistic quotient normalization as robust method to account for dilution of complex biological mixtures. Application in 1H NMR metabonomics. Anal Chem. 2006;78: 4281–4290. pmid:16808434
- View Article
- PubMed/NCBI
- Google Scholar
37. Xia J, Sinelnikov IV, Han B, Wishart DS. MetaboAnalyst 3.0—making metabolomics more meaningful. Nucleic Acids Res. 2015;43: W251–W257. pmid:25897128
- View Article
- PubMed/NCBI
- Google Scholar
38. Xia J, Wishart DS. MetPA: a web-based metabolomics tool for pathway analysis and visualization. Bioinformatics. 2010;26: 2342–2344. pmid:20628077
- View Article
- PubMed/NCBI
- Google Scholar
39. Cambiaghi A, Ferrario M, Masseroli M. Analysis of metabolomic data: tools, current strategies and future challenges for omics data integration. Brief Bioinform. 2017;18: 498–510. pmid:27075479
- View Article
- PubMed/NCBI
- Google Scholar
40. Hubert M, Rousseeuw PJ, Vanden Branden K. ROBPCA: A new approach to robust principal component analysis. Technometrics. 2005;47: 64–79.
- View Article
- Google Scholar
41. Meyer H, Weidmann H, Lalk M. Methodological approaches to help unravel the intracellular metabolome of Bacillus subtilis. Microb Cell Fact. 2013;12: 69. pmid:23844891
- View Article
- PubMed/NCBI
- Google Scholar
42. Rees CA, Smolinska A, Hill JE. The volatile metabolome of Klebsiella pneumoniae in human blood. J Breath Res. 2016;10: 027101. pmid:27163334
- View Article
- PubMed/NCBI
- Google Scholar
43. Rees CA, Franchina FA, Nordick KV, Kim PJ, Hill JE. Expanding the Klebsiella pneumoniae volatile metabolome using advanced analytical instrumentation for the detection of novel metabolites. J Appl Microbiol. 2017;122: 785–795. pmid:27930839
- View Article
- PubMed/NCBI
- Google Scholar
44. Lee S, Kim B, Yang J, Jeong D, Park S, Shin SH, et al. Comparative whole genome transcriptome and metabolome analyses of five Klebsiella pneumonia strains. Bioprocess Biosyst Eng. 2015;38: 2201–2219. pmid:26296382
- View Article
- PubMed/NCBI
- Google Scholar
45. Porreca AM, Sullivan KV, Gallagher JC. The epidemiology, evolution, and treatment of KPC-producing organisms. Curr Infect Dis Rep. 2018;20: 13. pmid:29730830
- View Article
- PubMed/NCBI
- Google Scholar
46. Dwyer DJ, Belenky PA, Yang JH, MacDonald IC, Martell JD, Takahashi N, et al. Antibiotics induce redox-related physiological alterations as part of their lethality. Proc Natl Acad Sci U S A. 2014;111: E2100E2109.
- View Article
- Google Scholar
47. Dwyer DJ, Collins JJ, Walker GC. Unraveling the physiological complexities of antibiotic lethality. Annu Rev Pharmacol Toxicol 2015;55: 313–332. pmid:25251995
- View Article
- PubMed/NCBI
- Google Scholar
48. Belenky P, Ye JD, Porter CB, Cohen NR, Lobritz MA, Ferrante T, et al. Bactericidal Antibiotics Induce Toxic Metabolic Perturbations that Lead to Cellular Damage. Cell Rep. 2015;13: 968–980. pmid:26565910
- View Article
- PubMed/NCBI
- Google Scholar
49. Dörries K, Schlueter R, Lalk M. Impact of antibiotics with various target sites on the metabolome of Staphylococcus aureus. Antimicrob Agents Chemother. 2014;58: 7151–7163. pmid:25224006
- View Article
- PubMed/NCBI
- Google Scholar
50. Keaton MA, Rosato RR, Plata KB, Singh CR, Rosato AE. Exposure of clinical MRSA heterogeneous strains to β-lactams redirects metabolism to optimize energy production through the TCA cycle. PLoS One. 2013;8: e71025. pmid:23940684
- View Article
- PubMed/NCBI
- Google Scholar
51. Nordmann P, Naas T, Poirel L. Global spread of Carbapenemase-producing Enterobacteriaceae. Emerg Infect Dis. 2011;17: 1791–1798 pmid:22000347
- View Article
- PubMed/NCBI
- Google Scholar
52. Schelli K, Zhong F, Zhu J. Comparative metabolomics revealing Staphylococcus aureus metabolic response to different antibiotics. Microb Biotechnol. 2017;10: 1764–1774. pmid:28815967
- View Article
- PubMed/NCBI
- Google Scholar
53. Smolinska A, Blanchet L, Buydens LM, Wijmenga SS. NMR and pattern recognition methods in metabolomics: from data acquisition to biomarker discovery: a review. Anal Chim Acta. 2012;750: 82–97. pmid:23062430
- View Article
- PubMed/NCBI
- Google Scholar
54. Hollywood K, Brison DR, Goodacre R. Metabolomics: current technologies and future trends. Proteomics. 2006;6: 4716–4723. pmid:16888765
- View Article
- PubMed/NCBI
- Google Scholar
55. Zhang B, Powers R. Analysis of bacterial biofilms using NMR-based metabolomics. Future Med Chem. 2012;4: 1273–1306. pmid:22800371
- View Article
- PubMed/NCBI
- Google Scholar
56. Trimigno A, Marincola FC, Dellarosa N, Picone G, Laghi L. Definition of food quality by NMR-based foodomics. Current Opinion in Food Science. 2015;4: 99–104.
- View Article
- Google Scholar
57. Carmeli Y, Akova M, Cornaglia G, Daikos GL, Garau J, Harbarth S, et al. Controlling the spread of carbapenemase-producing Gram-negatives: therapeutic approach and infection control. Clin Microbiol Infect. 2010;16: 102–111. pmid:20085604
- View Article
- PubMed/NCBI
- Google Scholar
58. Centers for Disease Control and Prevention (CDC). 2009. Guidance for control of infections with carbapenem-resistant or carbapenemase-producing Enterobacteriaceae in acute care facilities. MMWR Morb Mortal Wkly Rep. 58:256–260. pmid:19300408
- View Article
- PubMed/NCBI
- Google Scholar

[ref1] 1. Tzouvelekis LS, Markogiannakis A, Psichogiou M, Tassios PT, Daikos GL. Carbapenemases in Klebsiella pneumoniae and other Enterobacteriaceae: an evolving crisis of global dimensions. Clin Microbiol Rev. 2012;25: 682–707. pmid:23034326
View Article
PubMed/NCBI
Google Scholar

[2] View Article

[3] PubMed/NCBI

[4] Google Scholar

[ref2] 2. European Centre for Disease Prevention and Control. Antimicrobial resistance surveillance in Europe 2014. Annual Report of the European Antimicrobial Resistance Surveillance Network (EARS-Net). Stockholm: ECDC; 2015.

[ref3] 3. Giani T, Pini B, Arena F, Conte V, Bracco S, Migliavacca R, et al. Epidemic diffusion of KPC carbapenemase-producing Klebsiella pneumoniae in Italy: results of the first countrywide survey, 15 May to 30 June 2011. Euro Surveill. 2013; 18: pii = 20489.
View Article
Google Scholar

[7] View Article

[8] Google Scholar

[ref4] 4. Grundmann H, Livermore DM, Giske CG, Canton R, Rossolini GM, Campos J, et al. Carbapenem-non-susceptible Enterobacteriaceae in Europe: conclusions from a meeting of national experts. Euro Surveill. 2010;15: pii = 19711.
View Article
Google Scholar

[10] View Article

[11] Google Scholar

[ref5] 5. Diene SM, Rolain JM. Carbapenemase genes and genetic platforms in Gram-negative bacilli: Enterobacteriaceae, Pseudomonas and Acinetobacter species. Clin Microbiol Infect. 2014;20: 831–838. pmid:24766097
View Article
PubMed/NCBI
Google Scholar

[13] View Article

[14] PubMed/NCBI

[15] Google Scholar

[ref6] 6. Daikos GL, Tsaousi S, Tzouvelekis LS, Anyfantis I, Psichogiou M, Argyropoulou A, et al. Carbapenemase-producing Klebsiella pneumoniae bloodstream infections: lowering mortality by antibiotic combination schemes and the role of carbapenems. Antimicrob Agents Chemother. 2014;58: 2322–2328. pmid:24514083
View Article
PubMed/NCBI
Google Scholar

[17] View Article

[18] PubMed/NCBI

[19] Google Scholar

[ref7] 7. Monaco M, Giani T, Raffone M, Arena F, Garcia-Fernandez A, Pollini S, et al. Colistin resistance superimposed to endemic carbapenem-resistant Klebsiella pneumoniae: a rapidly evolving problem in Italy, November 2013 to April 2014. Euro Surveill. 2014;19: pii = 20939.
View Article
Google Scholar

[21] View Article

[22] Google Scholar

[ref8] 8. Akova M, Daikos GL, Tzouvelekis L, Carmeli Y. Interventional strategies and current clinical experience with carbapenemase-producing Gram-negative bacteria. Clin Microbiol Infect. 2012;18: 439–448. pmid:22507111
View Article
PubMed/NCBI
Google Scholar

[24] View Article

[25] PubMed/NCBI

[26] Google Scholar

[ref9] 9. De Rosa FG, Corcione S, Cavallo R, Di Perri G, Bassetti M. Critical issues for Klebsiella pneumoniae KPC-carbapenemase producing K. pneumoniae infections: a critical agenda. Future Microbiol. 2015;10: 283–294. pmid:25689539
View Article
PubMed/NCBI
Google Scholar

[28] View Article

[29] PubMed/NCBI

[30] Google Scholar

[ref10] 10. Findlay J, Hopkins KL, Meunier D, Woodford N. Evaluation of three commercial assays for rapid detection of genes encoding clinically relevant carbapenemases in cultured bacteria. J Antimicrob Chemother. 2015;70: 1338–1342. pmid:25630646
View Article
PubMed/NCBI
Google Scholar

[32] View Article

[33] PubMed/NCBI

[34] Google Scholar

[ref11] 11. Dortet L, Agathine A, Naas T, Cuzon G, Poirel L, Nordmann P. Evaluation of the RAPIDEC® CARBA NP, the Rapid CARB Screen® and the Carba NP test for biochemical detection of carbapenemase-producing Enterobacteriaceae. J Antimicrob Chemother. 2015;70: 3014–3022. pmid:26260131
View Article
PubMed/NCBI
Google Scholar

[36] View Article

[37] PubMed/NCBI

[38] Google Scholar

[ref12] 12. Sharma R, Eun Park T, Moy S. Ceftazidime-Avibactam: A Novel Cephalosporin/β-Lactamase Inhibitor Combination for the Treatment of Resistant Gram-Negative Organisms. Clin Ther. 2016;pii: S0149-2918(16)30041-8.
View Article
Google Scholar

[40] View Article

[41] Google Scholar

[ref13] 13. Martens E, Demain AL. The antibiotic resistance crisis, with a focus on the United States. J Antibiot (Tokyo). 2017;70: 520–526.
View Article
Google Scholar

[43] View Article

[44] Google Scholar

[ref14] 14. Tzouvelekis LS, Markogiannakis A, Psichogiou M, Tassios PT, Daikos GL. Carbapenemases in Klebsiella pneumoniae and other Enterobacteriaceae: an evolving crisis of global dimensions. Clin Microbiol Rev. 2012;25: 682–707. pmid:23034326
View Article
PubMed/NCBI
Google Scholar

[46] View Article

[47] PubMed/NCBI

[48] Google Scholar

[ref15] 15. Laghi L, Picone G, Capozzi F. Nuclear magnetic resonance for foodomics beyond food analysis. TrAC Trends in Analytical Chemistry. 2014;59: 93–102.
View Article
Google Scholar

[50] View Article

[51] Google Scholar

[ref16] 16. Parolin C, Foschi C, Laghi L, Zhu C, Banzola N, Gaspari V, et al. Insights Into Vaginal Bacterial Communities and Metabolic Profiles of Chlamydia trachomatis Infection: Positioning Between Eubiosis and Dysbiosis. Front Microbiol. 2018;9: 600. pmid:29643849
View Article
PubMed/NCBI
Google Scholar

[53] View Article

[54] PubMed/NCBI

[55] Google Scholar

[ref17] 17. Foschi C, Laghi L, D'Antuono A, Gaspari V, Zhu C, Dellarosa N, et al. Urine metabolome in women with Chlamydia trachomatis infection. PLoS One. 2018;13: e0194827. pmid:29566085
View Article
PubMed/NCBI
Google Scholar

[57] View Article

[58] PubMed/NCBI

[59] Google Scholar

[ref18] 18. Foschi C, Laghi L, Parolin C, Giordani B, Compri M, Cevenini R, et al. Novel approaches for the taxonomic and metabolic characterization of lactobacilli: Integration of 16S rRNA gene sequencing with MALDI-TOF MS and 1H-NMR. PLoS One. 2017;12: e0172483. pmid:28207855
View Article
PubMed/NCBI
Google Scholar

[61] View Article

[62] PubMed/NCBI

[63] Google Scholar

[ref19] 19. Ammons MC, Tripet BP, Carlson RP, Kirker KR, Gross MA, Stanisich JJ, et al. Quantitative NMR metabolite profiling of methicillin-resistant and methicillin-susceptible Staphylococcus aureus discriminates between biofilm and planktonic phenotypes. J Proteome Res. 2014;13: 2973–2985. pmid:24809402
View Article
PubMed/NCBI
Google Scholar

[65] View Article

[66] PubMed/NCBI

[67] Google Scholar

[ref20] 20. Yeom J, Shin JH, Yang JY, Kim J, Hwang GS. (1)H NMR-based metabolite profiling of planktonic and biofilm cells in Acinetobacter baumannii 1656–2. PLoS One. 2013;8: e57730. pmid:23483923
View Article
PubMed/NCBI
Google Scholar

[69] View Article

[70] PubMed/NCBI

[71] Google Scholar

[ref21] 21. Chen T, Sheng J, Fu Y, Li M, Wang J, Jia AQ. 1H NMR-Based Global Metabolic Studies of Pseudomonas aeruginosa upon Exposure of the Quorum Sensing Inhibitor Resveratrol. J Proteome Res. 2017;16: 824–830. pmid:28094526
View Article
PubMed/NCBI
Google Scholar

[73] View Article

[74] PubMed/NCBI

[75] Google Scholar

[ref22] 22. Antti H, Fahlgren A, Näsström E, Kouremenos K, Sundén-Cullberg J, Guo Y, et al. Metabolic profiling for detection of Staphylococcus aureus infection and antibiotic resistance. PLoS One. 2013;8: e56971. pmid:23451124
View Article
PubMed/NCBI
Google Scholar

[77] View Article

[78] PubMed/NCBI

[79] Google Scholar

[ref23] 23. Liebeke M, Dörries K, Zühlke D, Bernhardt J, Fuchs S, Pané-Farré J, et al. A metabolomics and proteomics study of the adaptation of Staphylococcus aureus to glucose starvation. Mol Biosyst. 2011;7: 1241–1253. pmid:21327190
View Article
PubMed/NCBI
Google Scholar

[81] View Article

[82] PubMed/NCBI

[83] Google Scholar

[ref24] 24. Paterson DL, Bonomo RA. Extended-spectrum beta-lactamases: a clinical update. Clin Microbiol Rev. 2005;18: 657–686. pmid:16223952
View Article
PubMed/NCBI
Google Scholar

[85] View Article

[86] PubMed/NCBI

[87] Google Scholar

[ref25] 25. Ambretti S, Gaibani P, Berlingeri A, Cordovana M, Tamburini MV, Bua G, et al. Evaluation of phenotypic and genotypic approaches for the detection of class A and class B carbapenemases in Enterobacteriaceae. Microb Drug Resist. 2013;19: 212–215. pmid:23336577
View Article
PubMed/NCBI
Google Scholar

[89] View Article

[90] PubMed/NCBI

[91] Google Scholar

[ref26] 26. van Dijk K, Voets GM, Scharringa J, Voskuil S, Fluit AC, Rottier WC, et al. A disc diffusion assay for detection of class A, B and OXA-48 carbapenemases in Enterobacteriaceae using phenyl boronic acid, dipicolinic acid and temocillin. Clin Microbiol Infect. 2014;20: 345–349. pmid:23927659
View Article
PubMed/NCBI
Google Scholar

[93] View Article

[94] PubMed/NCBI

[95] Google Scholar

[ref27] 27. CLSI. Performance Standards for Antimicrobial Susceptibility Testing. 27th ed. CLSI supplement M100. Wayne, PA: Clinical and Laboratory Standards Institute; 2017.

[ref28] 28. Park C, Yun S, Lee SY, Park K, Lee J. Metabolic profiling of Klebsiella oxytoca: evaluation of methods for extraction of intracellular metabolites using UPLC/Q-TOF-MS. Appl Biochem Biotechnol. 2012;167: 425–438. pmid:22555499
View Article
PubMed/NCBI
Google Scholar

[98] View Article

[99] PubMed/NCBI

[100] Google Scholar

[ref29] 29. Schelli K, Rutowski J, Roubidoux J, Zhu J. Staphylococcus aureus methicillin resistance detected by HPLC-MS/MS targeted metabolic profiling. J Chromatogr B Analyt Technol Biomed Life Sci. 2017;1047: 124–130. pmid:27316783
View Article
PubMed/NCBI
Google Scholar

[102] View Article

[103] PubMed/NCBI

[104] Google Scholar

[ref30] 30. Ruscio J, Roche B. Determining the number of factors to retain in an exploratory factor analysis using comparison data of known factorial structure. Psychol Assess. 2012;24: 282–292. pmid:21966933
View Article
PubMed/NCBI
Google Scholar

[106] View Article

[107] PubMed/NCBI

[108] Google Scholar

[ref31] 31. Ghasemi A, Zahediasl S. Normality tests for statistical analysis: a guide for non-statisticians. Int J Endocrinol Metab. 2012;10: 486–489. pmid:23843808
View Article
PubMed/NCBI
Google Scholar

[110] View Article

[111] PubMed/NCBI

[112] Google Scholar

[ref32] 32. Parolin C, Marangoni A, Laghi L, Foschi C, Ñahui Palomino RA, Calonghi N, et al. Isolation of Vaginal Lactobacilli and Characterization of Anti-Candida Activity. PLoS One. 2015;10: e0131220. pmid:26098675
View Article
PubMed/NCBI
Google Scholar

[114] View Article

[115] PubMed/NCBI

[116] Google Scholar

[ref33] 33. Kneen MA, Annegarn HJ. Algorithm for fitting XRF, SEM and PIXE X-ray spectra backgrounds. Nucl Instrum Methods Phys Res B. 1996;109: 209–213.
View Article
Google Scholar

[118] View Article

[119] Google Scholar

[ref34] 34. Liland KH, Almøy T, Mevik BH. Optimal choice of baseline correction for multivariate calibration of spectra. Appl Spectrosc. 2010;64: 1007–1016. pmid:20828437
View Article
PubMed/NCBI
Google Scholar

[121] View Article

[122] PubMed/NCBI

[123] Google Scholar

[ref35] 35. Wishart DS, Tzur D, Knox C, Eisner R, Guo AC, Young N, et al. HMDB: the human metabolome database. Nucleic Acids Res. 2007;35: D521–D526. pmid:17202168
View Article
PubMed/NCBI
Google Scholar

[125] View Article

[126] PubMed/NCBI

[127] Google Scholar

[ref36] 36. Dieterle F, Ross A, Schlotterbeck G, Senn H. Probabilistic quotient normalization as robust method to account for dilution of complex biological mixtures. Application in 1H NMR metabonomics. Anal Chem. 2006;78: 4281–4290. pmid:16808434
View Article
PubMed/NCBI
Google Scholar

[129] View Article

[130] PubMed/NCBI

[131] Google Scholar

[ref37] 37. Xia J, Sinelnikov IV, Han B, Wishart DS. MetaboAnalyst 3.0—making metabolomics more meaningful. Nucleic Acids Res. 2015;43: W251–W257. pmid:25897128
View Article
PubMed/NCBI
Google Scholar

[133] View Article

[134] PubMed/NCBI

[135] Google Scholar

[ref38] 38. Xia J, Wishart DS. MetPA: a web-based metabolomics tool for pathway analysis and visualization. Bioinformatics. 2010;26: 2342–2344. pmid:20628077
View Article
PubMed/NCBI
Google Scholar

[137] View Article

[138] PubMed/NCBI

[139] Google Scholar

[ref39] 39. Cambiaghi A, Ferrario M, Masseroli M. Analysis of metabolomic data: tools, current strategies and future challenges for omics data integration. Brief Bioinform. 2017;18: 498–510. pmid:27075479
View Article
PubMed/NCBI
Google Scholar

[141] View Article

[142] PubMed/NCBI

[143] Google Scholar

[ref40] 40. Hubert M, Rousseeuw PJ, Vanden Branden K. ROBPCA: A new approach to robust principal component analysis. Technometrics. 2005;47: 64–79.
View Article
Google Scholar

[145] View Article

[146] Google Scholar

[ref41] 41. Meyer H, Weidmann H, Lalk M. Methodological approaches to help unravel the intracellular metabolome of Bacillus subtilis. Microb Cell Fact. 2013;12: 69. pmid:23844891
View Article
PubMed/NCBI
Google Scholar

[148] View Article

[149] PubMed/NCBI

[150] Google Scholar

[ref42] 42. Rees CA, Smolinska A, Hill JE. The volatile metabolome of Klebsiella pneumoniae in human blood. J Breath Res. 2016;10: 027101. pmid:27163334
View Article
PubMed/NCBI
Google Scholar

[152] View Article

[153] PubMed/NCBI

[154] Google Scholar

[ref43] 43. Rees CA, Franchina FA, Nordick KV, Kim PJ, Hill JE. Expanding the Klebsiella pneumoniae volatile metabolome using advanced analytical instrumentation for the detection of novel metabolites. J Appl Microbiol. 2017;122: 785–795. pmid:27930839
View Article
PubMed/NCBI
Google Scholar

[156] View Article

[157] PubMed/NCBI

[158] Google Scholar

[ref44] 44. Lee S, Kim B, Yang J, Jeong D, Park S, Shin SH, et al. Comparative whole genome transcriptome and metabolome analyses of five Klebsiella pneumonia strains. Bioprocess Biosyst Eng. 2015;38: 2201–2219. pmid:26296382
View Article
PubMed/NCBI
Google Scholar

[160] View Article

[161] PubMed/NCBI

[162] Google Scholar

[ref45] 45. Porreca AM, Sullivan KV, Gallagher JC. The epidemiology, evolution, and treatment of KPC-producing organisms. Curr Infect Dis Rep. 2018;20: 13. pmid:29730830
View Article
PubMed/NCBI
Google Scholar

[164] View Article

[165] PubMed/NCBI

[166] Google Scholar

[ref46] 46. Dwyer DJ, Belenky PA, Yang JH, MacDonald IC, Martell JD, Takahashi N, et al. Antibiotics induce redox-related physiological alterations as part of their lethality. Proc Natl Acad Sci U S A. 2014;111: E2100E2109.
View Article
Google Scholar

[168] View Article

[169] Google Scholar

[ref47] 47. Dwyer DJ, Collins JJ, Walker GC. Unraveling the physiological complexities of antibiotic lethality. Annu Rev Pharmacol Toxicol 2015;55: 313–332. pmid:25251995
View Article
PubMed/NCBI
Google Scholar

[171] View Article

[172] PubMed/NCBI

[173] Google Scholar

[ref48] 48. Belenky P, Ye JD, Porter CB, Cohen NR, Lobritz MA, Ferrante T, et al. Bactericidal Antibiotics Induce Toxic Metabolic Perturbations that Lead to Cellular Damage. Cell Rep. 2015;13: 968–980. pmid:26565910
View Article
PubMed/NCBI
Google Scholar

[175] View Article

[176] PubMed/NCBI

[177] Google Scholar

[ref49] 49. Dörries K, Schlueter R, Lalk M. Impact of antibiotics with various target sites on the metabolome of Staphylococcus aureus. Antimicrob Agents Chemother. 2014;58: 7151–7163. pmid:25224006
View Article
PubMed/NCBI
Google Scholar

[179] View Article

[180] PubMed/NCBI

[181] Google Scholar

[ref50] 50. Keaton MA, Rosato RR, Plata KB, Singh CR, Rosato AE. Exposure of clinical MRSA heterogeneous strains to β-lactams redirects metabolism to optimize energy production through the TCA cycle. PLoS One. 2013;8: e71025. pmid:23940684
View Article
PubMed/NCBI
Google Scholar

[183] View Article

[184] PubMed/NCBI

[185] Google Scholar

[ref51] 51. Nordmann P, Naas T, Poirel L. Global spread of Carbapenemase-producing Enterobacteriaceae. Emerg Infect Dis. 2011;17: 1791–1798 pmid:22000347
View Article
PubMed/NCBI
Google Scholar

[187] View Article

[188] PubMed/NCBI

[189] Google Scholar

[ref52] 52. Schelli K, Zhong F, Zhu J. Comparative metabolomics revealing Staphylococcus aureus metabolic response to different antibiotics. Microb Biotechnol. 2017;10: 1764–1774. pmid:28815967
View Article
PubMed/NCBI
Google Scholar

[191] View Article

[192] PubMed/NCBI

[193] Google Scholar

[ref53] 53. Smolinska A, Blanchet L, Buydens LM, Wijmenga SS. NMR and pattern recognition methods in metabolomics: from data acquisition to biomarker discovery: a review. Anal Chim Acta. 2012;750: 82–97. pmid:23062430
View Article
PubMed/NCBI
Google Scholar

[195] View Article

[196] PubMed/NCBI

[197] Google Scholar

[ref54] 54. Hollywood K, Brison DR, Goodacre R. Metabolomics: current technologies and future trends. Proteomics. 2006;6: 4716–4723. pmid:16888765
View Article
PubMed/NCBI
Google Scholar

[199] View Article

[200] PubMed/NCBI

[201] Google Scholar

[ref55] 55. Zhang B, Powers R. Analysis of bacterial biofilms using NMR-based metabolomics. Future Med Chem. 2012;4: 1273–1306. pmid:22800371
View Article
PubMed/NCBI
Google Scholar

[203] View Article

[204] PubMed/NCBI

[205] Google Scholar

[ref56] 56. Trimigno A, Marincola FC, Dellarosa N, Picone G, Laghi L. Definition of food quality by NMR-based foodomics. Current Opinion in Food Science. 2015;4: 99–104.
View Article
Google Scholar

[207] View Article

[208] Google Scholar

[ref57] 57. Carmeli Y, Akova M, Cornaglia G, Daikos GL, Garau J, Harbarth S, et al. Controlling the spread of carbapenemase-producing Gram-negatives: therapeutic approach and infection control. Clin Microbiol Infect. 2010;16: 102–111. pmid:20085604
View Article
PubMed/NCBI
Google Scholar

[210] View Article

[211] PubMed/NCBI

[212] Google Scholar

[ref58] 58. Centers for Disease Control and Prevention (CDC). 2009. Guidance for control of infections with carbapenem-resistant or carbapenemase-producing Enterobacteriaceae in acute care facilities. MMWR Morb Mortal Wkly Rep. 58:256–260. pmid:19300408
View Article
PubMed/NCBI
Google Scholar

[214] View Article

[215] PubMed/NCBI

[216] Google Scholar

Figures

Abstract

Introduction

Materials and methods

Bacterial strains

Sample preparation

Metabolomic analysis

Metabolic pathway analysis

Statistical analysis

Results

Bacterial strains

Metabolomic analysis

(i) Metabolome of K. pneumoniae in basal conditions

(ii) Metabolome of K. pneumoniae strains under ‘meropenem stress’

(iii) Analysis of the metabolic response to meropenem of wild-type and KPC-producing strains

Discussion

Conclusions

Supporting information

S1 Fig. Bacterial growth curves at different sub-lethal levels of meropenem MIC.

S2 Fig.

S1 Table. List of Klebsiella pneumoniae strains included in the study.

S2 Table. List of the molecules found in the intracellular and extracellular metabolome of Klebsiella pneumoniae strains.

S3 Table. List of intracellular molecules that significantly differed after meropenem exposure in the group of selected K. pneumoniae isolates, irrespective of the carbapenemase production.

S4 Table. List of extracellular molecules that significantly differed after meropenem exposure in the group of selected K. pneumoniae isolates, irrespective of their carbapenemase production.

Acknowledgments

References