The polymyxin B-induced transcriptomic response of a clinical, multidrug-resistant Klebsiella pneumoniae involves multiple regulatory elements and intracellular targets
- Pablo Ivan Pereira Ramos1, 2,
- Márlon Grégori Flores Custódio1,
- Guadalupe del Rosario Quispe Saji1,
- Thiago Cardoso1,
- Gisele Lucchetti da Silva1,
- Graziela Braun3,
- Willames M. B. S. Martins3,
- Raquel Girardello3,
- Ana Tereza Ribeiro de Vasconcelos1,
- Elmer Fernández4,
- Ana Cristina Gales3Email author and
- Marisa Fabiana Nicolás1Email author
© The Author(s). 2016
Published: 25 October 2016
The emergence of multidrug-resistant Klebsiella pneumoniae is a major public health concern. Many K. pneumoniae infections can only be treated when resorting to last-line drugs such as polymyxin B (PB). However, resistance to this antibiotic is also observed, although insufficient information is described on its mode of action as well as the mechanisms used by resistant bacteria to evade its effects. We aimed to study PB resistance and the influence of abiotic stresses in a clinical K. pneumoniae strain using whole transcriptome profiling.
We sequenced 12 cDNA libraries of K. pneumoniae Kp13 bacteria, from two biological replicates of the original strain Kp13 (Kp13) and five derivative strains: induced high-level PB resistance in acidic pH (Kp13pH), magnesium deprivation (Kp13Mg), high concentrations of calcium (Kp13Ca) and iron (Kp13Fe), and a control condition with PB (Kp13PolB). Our results show the involvement of multiple regulatory loci that differentially respond to each condition as well as a shared gene expression response elicited by PB treatment, and indicate the participation of two-regulatory components such as ArcA-ArcB, which could be involved in re-routing the K. pneumoniae metabolism following PB treatment. Modules of co-expressed genes could be determined, which correlated to growth in acid stress and PB exposure. We hypothesize that polymyxin B induces metabolic shifts in K. pneumoniae that could relate to surviving against the action of this antibiotic.
We obtained whole transcriptome data for K. pneumoniae under different environmental conditions and PB treatment. Our results supports the notion that the K. pneumoniae response to PB exposure goes beyond damaged membrane reconstruction and involves recruitment of multiple gene modules and intracellular targets.
Antibiotic discovery represented one of the most significant events in human health improvement. The antibiotic ‘golden age’, when many classes were identified and commercialized, extended from 1950 to 1960 . Parallel to antibiotic availability, there were ever increasing numbers of bacteria that could overcome treatment and displayed a resistant phenotype. Extreme cases involve bacteria that display resistance against all clinically available agents . Although the emergence of resistance pre-dates the antibiotic era, widespread multi-drug resistance (MDR) is now a major concern that limit therapeutic options in the treatment of many bacterial infections.
Polymyxins, including polymyxin B and colistin, are cationic antimicrobial lipopeptides (CAMPs) that have been available for decades, but due to their nephrotoxicity and the lack of pharmacodynamics and pharmacokinetics information, have had its use greatly diminished [3, 4]. The mode of action of these drugs is also not completely elucidated, but generally involves interaction with lipopolysaccharides (LPS) located in the outer membrane of Gram-negative bacteria, competing with the calcium and magnetic ions that stabilize LPS, allowing for drug uptake to the cell interior . It has been reported that PB activity may involve inhibition of vital respiratory enzymes located in the inner membrane, such as NADH-quinone oxidoreductase (NDH-2) , thereby hampering bacterial respiration.
Recently, polymyxins have re-emerged, now as ‘last resort’ drugs, used in patients with difficult-to-treat infections caused by MDR Gram-negative bacteria including K. pneumoniae, Pseudomonas aeruginosa and Acinetobacter baumannii. Resistance to polymyxins, however, has also been described in many settings [6, 7].
In order to better understand the mechanisms employed by MDR bacteria to evade the action of polymyxins, we applied whole transcriptomic profiling of a clinical K. pneumoniae strain with induced, high-level resistance to polymyxin B. We previously reported the complete genome of this strain , and we identified that it has a truncated mgrB gene, consistent with reports that mutations in this gene represent one of the various mechanisms for acquired colistin resistance [9–11]. Despite Kp13 being already resistant to polymyxin B, we aimed to induce additional adaptive responses by growth in high polymyxin B concentrations while also probing in vitro the effects of diverse environmental stimuli. Variations in abiotic stimuli govern many transcriptional regulatory responses, including in two-component regulatory systems (TCRSs) that may play a role in drug resistance. For instance, varying concentrations of divalent cations (such as Mg2+ and Ca2+) alter the expression of important global regulators such as the PhoP-PhoQ TCRS, which regulates over 40 genes and has been implicated in colistin resistance . On the other hand, acid pH and iron supplementation were shown to induce colistin resistance in A. baumannii probably involving alterations in regulatory circuits .
While genome analyses provide a complete picture of the genes present (and potentially expressed), only through the use of high throughput expression profiling techniques such as RNA sequencing (RNA-seq), as performed in this report, can we find consistent expression patterns that can aid in pinpointing the intracellular targets and metabolic processes related to polymyxin B mode of action and resistance.
Bacterial strain, higher-level polymyxin B resistance induction and growth conditions
Klebsiella pneumoniae subsp. pneumoniae Kp13 (hereafter referred as Kp13) was isolated in 2009 from a patient during an event of nosocomial outbreak due to KPC-2-producing K. pneumoniae bacteria in an intensive care unit from the Hospital Universitário (UEL), South Brazil. The complete genome of this strain is comprised by one chromosome and six plasmids and was completely sequenced by our group in a previous work . Kp13 harbors multiple virulence and resistance determinants in its genome, and presents resistant to many antibiotics [8, 14].
Kp13 was originally resistant to polymyxin B (PB)  at a minimum inhibitory concentration (MIC) of 32 μg mL−1 (EUCAST breakpoint ). We have induced an increased, high-level resistance to this antibiotic by growing the bacteria in solid Luria-Bertani medium (LB, Oxoid, Basingstoke, England) in the presence of crescent polymyxin B (Sigma-Aldrich, St. Louis, MO, USA) concentrations and passaging the bacteria in serial dilutions of PB beginning with a concentration of 8 μg mL−1 up to 64 μg mL−1. Before and after the induction of resistance, PB MICs were confirmed by CLSI broth microdilutions.
Strains used in RNA sequencing experiments and respective growth conditions
Concentration (mg L−1)
MICb (μg mL−1)
Original strain, no polymyxin B supplementation
LB supplemented with 4 μg mL−1 of polymyxin B
PB + high calcium
LB supplemented with 4 μg mL−1 of polymyxin B plus high [Ca2+]
PB + no magnesium
LB supplemented with 4 μg mL−1 of polymyxin B without Mg2+
PB + high iron
LB supplemented with 4 μg mL−1 of polymyxin B plus high [Fe2+]
PB + low pH
LB supplemented with 4 μg mL−1 of polymyxin B in acid pH
Total RNA extraction and Illumina sequencing
For each condition, bacterial cultures were grown until final-log phase and cells were all harvested at OD 600 nm. Total RNA extraction was performed using RNeasy Mini Kit (Qiagen) with DNAse treatment (Qiagen). Enrichment for mRNA was performed using the MICROBExpress Bacterial mRNA Purification (Ambion), and rRNA were removed using Ribo-Zero kit (Epicentre). Sequencing of two biological replicates per condition was performed in an Illumina HiSeq 2500 (Fasteris, Geneva, Switzerland) using a paired-end strategy.
Expression validation by real time PCR
Primers sequence used in qRT-PCR experiments
5′ GCTGAACGTCCAACTGAAAG 3′
5′ GGAGGATTGCTGTTCGAGC 3′
5′ TGCCGGATGAAGACGGACTA 3′
5′ AGGGAGATCACCTGTGAGGC 3′
5′ GCTGATCCAGCGTCTCGATC 3′
5′ CAACAGCACCTGCTGGTAGC 3′
5′ CAGCTCGTGTCGTGAGATGT 3′
5′ CGTAAGGGCCATGATGACTT 3′
Growth curve and analysis of pH variance
Strains Kp13PolB and Kp13 (non-induced) were subcultivated in the absence or presence of polymyxin B to evaluate differences in bacterial fitness. The growth of bacterial strains was monitored until stationary phase of growth was reached. The strains were cultured in the absence or presence of polymyxin B sulfate 4 μg mL−1 in LB broth overnight. Three milliliters of each culture were diluted in 447 mL of fresh LB broth in absence or presence of polymyxin B, and incubated under constant shaking at 37 °C. Ten milliliters of each bacterial suspension was removed every hour. One milliliter was transferred to a cuvette for turbidity measure at 600 nm wavelength, while the remaining 9 mL were used to quantify the pH variation during the bacteria growth.
Raw read preprocessing and quality control (QC) were performed with FastQC , which allows assessment of diverse quality metrics of the sequencing data. Removal of adapters and low-quality reads were performed with Skewer . We applied the normalization method Removal of Unwanted Variation (RUV) in order to properly correct for batch effects in experimental data obtained through RNA sequencing . For this approach we used the MLST genes of Klebsiella, available at the Institut Pasteur MLST database website , as control genes among the cDNA libraries data . Reads passing QC were aligned to the K. pneumoniae Kp13 reference genome (NCBI BioProject accession PRJNA78291) using Top-Hat2 , and SAMtools  was used to parse alignments and QC of the mappings. Gene counts were obtained using HTSeq  and differential gene expression was evaluated with edgeR , which models expression data as negative binomial distributed. The PB condition (Table 1) was the reference condition for comparisons performed in the RNA-seq analyses unless otherwise noted in the text.
Co-expression network construction
Gene co-expression networks were created using weighted gene co-expression network analysis (WGCNA) package in R . Variance-stabilized count data were used as a measure of gene expression. The biweight midcorrelation algorithm implemented in the bicor function in WGCNA was used to calculate expression correlation between each pair of genes. This correlation measure is similar to Pearson’s correlation with the advantage of being more robust to outliers . The correlation matrix was transformed into a weighted adjacency matrix representing connection strengths using a soft-thresholding approach by applying a power transformation with an appropriate β parameter, f(x) = x β. β is chosen such that the resulting adjacency matrix leads to a weighted co-expression network that displays approximately scale-free topology. In this framework, a value of β = 14 was used in our analyses, which resulted in a scale-free model fitting with R 2 >0.80. Once the adjacency matrix was constructed, we derived from it a topological overlap matrix (TOM), which takes gene expression connectivity into account. 1-TOM was used as dissimilarity measure for hierarchical clustering and module detection. Module assignment was determined using the dynamic tree cut algorithm within WGCNA . Module eigengenes were tested for association with the different experimental conditions and those displaying significant correlation, defined as Pearson's |r| ≥ 0.6; p-value <0.05, were more thoroughly studied by means of Gene Ontology (GO) enrichment analysis.
GO enrichment analysis
Genes belonging to each module of interest were submitted to ClueGO v. 2.2.4  for Cytoscape v.3.3.0 , which performs Gene Ontology enrichment analysis and cluster detection. Only biological process terms were considered, and GO terms were grouped if their kappa score was ≥0.4. The statistical test chosen to identify enriched biological process terms was the right-sided hypergeometric test, considering as significant terms with adjusted p-value <0.05 and comprising at least three genes. The Benjamini-Hochberg procedure was applied to control the false discovery rate. To study both the more general biological processes associated to each module as well as the more specific terms, the “GO tree interval” parameter was set to either (min = 2, max = 5) or (min = 7, max = 15), respectively.
Results and discussion
Evaluation of the global transcriptional response of K. pneumoniae to increased polymyxin B resistance
Firstly, we investigated the general transcriptional response of K. pneumoniae Kp13 in view of increased resistance to polymyxin B (PB). Prior to high-level experimental induction of resistance, Kp13 was already PB-resistant (MIC 32 μg mL−1) and had a truncated mgrB gene due to transposase insertion, a mechanism already described in the literature as providing PB resistance [9, 11]. We were interested in detecting additional mechanisms of resistance to this antibiotic as well as understand the effect of PB treatment combined with diverse abiotic stresses in the global gene expression response of K. pneumoniae.
Among the shared PB-dependent up-regulated genes are various transcriptional factors (TFs). These TFs (and corresponding locus) include the OmpR regulator (KP13_00688) as well as the argR gene (KP13_31778) and genes encoding for a putative HTH-transcriptional regulator (KP13_03008; KP13_00251). Also detected was arcB (KP13_03225), part of the TCRS arcAB which has been shown to be highly pleiotropic, and modulates the expression of genes encoding for proteins with membrane modification functions as well as TCA cycle enzymes depending on oxygen levels or redox state [30, 31]. Other TFs detected included rpoH (KP13_08477), rpoN (KP13_03232). rpoH, encoding for the RNA polymerase σ32 factor is involved in stress response, mobilizing genes that constitute the heat-shock stimulon, such as those encoding for chaperones and proteases , while rpoN codes the σ54 factor, involved in the regulation of nitrogen utilization genes. Among the upregulated genes we also highlight wzzE (KP13_13130) and the ORF KP13_00225, encoding for a lipopolysaccharide biosynthesis protein and a putative inner membrane protein, respectively. Overall, the combined evidence show that the main up-regulated functions in PB treatment are those related to membrane functionality and maintenance. This is in line with the primary mode of action reported for this antibiotic, which permeabilizes the Gram-negative bacterial OM via direct interaction with the LPS .
Examination of down-regulated genes revealed that both the glycerol uptake (glpF and glpT) and metabolic pathways related to anaerobic catabolism of glycerol (dhaD, dhaK, dhaM, glpA, glpB, glpC, glpD, glpQ and plsY) are repressed by exposure to polymyxin B independently of environmental stimuli. Interestingly, a secondary mode of action for PB has been described and involves inhibitory activity against NDH-2, a vital respiratory enzyme in the bacterial inner membrane . This would lead to the activation of genes involved in the response to decreased oxidative power within the cell, including arcB, as we have observed. Besides repressing TCA genes it is also known that the ArcA-ArcB TCRS acts repressing genes encoding primary dehydrogenases such as glpD , which would explain our findings.
Respiratory/fermentative pathways and hydroxyl radical scavenging form part of the response of K. pneumoniae to polymyxin B and abiotic stimuli
Relative expression of quinone biosynthesis and oxidoreductase genes belonging to the K. pneumoniae Kp13 respiratory chain
Primary dehydrogenases (DH):
Formate dehydrogenase-N subunit alpha (FDH-N subunit alpha)
Formate dehydrogenase-N subunit beta (FDH-N subunit beta)
Formate dehydrogenase-N subunit (FDH-N subunit gamma)
Formate hydrogenlyase H
Formate hydrogenlyase regulatory protein hycA
Formate hydrogenlyase subunit 2 to 7
0.9, 1.6, 0.9, 0.8, 0.9, 1.2
Formate hydrogenlyase maturation protein HycH
NADH-quinone oxidoreductase (NDH-1)
0.1, 0, 0.2, 0, 0, −0.2, 0, −0.3, −0.1, −0.2, −0.1, −0.1, −0.1
NADH-quinone oxidoreductase (NDH-2)
Glycerol-3-P DH N (subunit A)
Glycerol-3-P DH N (subunit B)
Glycerol-3-P DH N (subunit C)
D -Lactate DH
L -Lactate DH
lctD or lldD
D -Amino acid DH
0.3, 0.1, 0.1, 0, −0.1
Naphthoquinone menaquinone (MK)
Quinol oxidase bo 3
Quinol oxidase bd
Quinol oxidase III (Cyx)
Nitrate reductase A
0.2, 0.5, 0.6, 0.3
Nitrate reductase Z
0.6, 0.3, 0.7, 0.5
Nitrate reductase, periplasmic
0.8, 0.5, −0.2
Thus, our results show that polymyxin B affects the expression of several genes encoding for respiratory chain enzymes, which may lead to a shift in alternative respiratory and fermentation pathways for bacterial energy metabolism and to reoxidize the NADH produced by glycolysis. Therefore, we investigated the response of genes encoding for enzymes in these alternative pathways, both by the effect of PB with or without abiotic stimuli. In Additional file 3 we depict a model, according to the data obtained, of carbon flux in Kp13 under fermentative conditions by stresses imposed by PB (Additional file 3, A) or PB plus abiotic stimuli (high calcium/iron, magnesium depletion or low pH) (Additional file 3, B). Consequently to stress imposed by PB without abiotic stimuli the expression of genes encoding for the three enzymatic components of 2-oxoglutarate dehydrogenase complex were slightly repressed (sucA: log2FC = −0.1 with FDR≤0.05, sucB: log2FC = −0.3 with FDR≤0.01 and lpdA: log2FC = −0.2 with FDR≤0.05). This repression was more pronounced with stress of PB in combination with abiotic stimuli, in particular in PB+high Ca (sucA: log2FC = −1.4, sucB: log2FC = −0.8 and lpdA: log2FC = −1.6, all with FDR≤0.01) and PB+low pH (sucA: log2FC = −2.4, sucB: log2FC = −1.6 and lpdA: log2FC = −1.6, all with FDR≤0.01) (Table 3). One explanation for this repression may relate to the observation that arcA-arcB showed high expression level following the pressure of PB without abiotic stimuli, and this TCRS represses the expression of the genes involved in the TCA cycle, such as the sucABCD system . Under low expression of the genes encoding for the 2-oxoglutarate dehydrogenase complex the citric acid cycle may be interrupted, leading to a shift from aerobic pathways to the fermentation route . Upon PB stress it seems that lactate and acetolactate can be produced actively because ldhA and lvG showed increased expression following PB treatment, almost independently of abiotic stimuli (Additional file 3, A; Table 3). Conversely, the influence of abiotic changes together with PB induces up-regulation of genes involved in the production of ethanol, acetate and ATP, formate, butane 2,3-ediol and 1,3-propanediol (Additional file 3, B; Table 3). The synthesis of acid by-products would be observable as a lowering in the pH of the medium. To test this hypothesis, we performed an experiment of K. pneumoniae Kp13 growth under polymyxin B exposure (analogous to the PB condition in Table 1) and absence, concurrently measuring the medium pH at hourly intervals (see Additional file 4). We observed that, while there is a general trend of acidification of the culture pH, this lowering is more pronounced in bacteria grown under PB stress (Additional file 4), suggesting that PB exposure induces metabolic shifts in K. pneumoniae that lead to synthesis of acidic by-products.
The metabolism of glycerol by Kp13 under the pressure of PB (Additional file 3, C) and PB plus abiotic stimuli (Additional file 3, D) were also studied. Notably, it seems that all genes of this pathway showed decreased expression following PB treatment, whereas PB in combination with abiotic stimuli induces up-regulation of all genes for production of 1,3-propanediol.
In this sense, the shift in synthesis of metabolic by-products might reflect a decrease in NADH oxidation via the respiratory chain, which is compensated by the NADH consuming lactate or acetolactate (under PB stress), and ethanol or 1,3-propanediol (under PB+abiotic stimuli stresses) production.
On the other hand, disruption in the electron transport chain induces the production of superoxide (HO2 .-) that is converted in hydrogen peroxide (H2O2), which in turn participate in the Fenton reaction and additionally induce synthesis of hydroxyl radicals [38, 39]. Large concentrations of superoxide and hydrogen peroxide trigger bacterial antioxidant defenses, but small amounts of hydroxyl radicals (OH−) do not elicit adequate defense mechanisms and could quickly become fatal [38, 40]. In Gram-negative bacteria, the intracellular superoxide dismutases Mn-SodA and Fe-SodB, catalyze the dismutation of metabolic sources of superoxide to hydrogen peroxide and oxygen . Then, in the presence of transition iron ions (Fe2+), the hydrogen peroxide can be converted to hydroxyl radicals via Fenton chemistry: Fe2+ + H2O2 → Fe3+ + OH. + OH−. In bacteria, the scavenging enzymes that can resolve the hydrogen peroxide production are glutathione peroxidases and catalases [39, 40]. For that, glutathione peroxidases by oxidation of glutathione (GSH): 2GSH + H2O2 → GSSG + 2H2O and catalase via the reaction 2H2O2 → O2 + 2H2O .
In our RNA-seq analysis, we observed that sodA (KP13_00587) was up-regulated in PB without the influence of abiotic changes, whereas sodB (KP13_05143) appeared slightly overexpressed in high iron and magnesium deprivation conditions (Additional file 2). Genes encoding for a glutathione peroxidase (KP13_31955) and Catalase HPII (katE, KP13_05384) were overexpressed in polymyxin B independently of abiotic changes. While katG (KP13_04579) encoding for a catalase-peroxidase bufunctional enzyme appeared overexpressed in high calcium/iron conditions (Additional file 2).
One of the mechanisms used by bacteria to cope with the oxidative stress induced by Fenton products involve the intracellular regulation of iron metabolism mediated by bacterioferritin (Bfr) or ferritin-like proteins. These proteins sequester intracellular Fe2+ ions and stores it in the form of Fe3+, where one hydrogen peroxide oxidizes two Fe2+ ions, thus preventing hydroxyl radical production by the Fenton reaction . Here, we identified a gene encoding for ferritin-like protein (KP13_13049) overexpressed in PB considering all abiotic changes, with a specially consistent up-regulation in PB with high iron condition (log2FC = 1.6) (Additional file 2). It has been reported that the production of hydroxyl radicals induced by polymyxin treatment is concurrent with rapid killing of Gram-negative A. baumannii .
Additionally, Fenton oxidants cause oxidation of DNA and free nucleotides, which can induce mismatched pairing, mutagenesis, and DNA breaks. We detected significant elevated relative expression of genes encoding for DNA damage repair proteins in polymyxin B with acid pH (recF - KP13_00050 - log2FC = 2.8; recO - KP13_00815 - log2FC = 0.7, recR - KP13_03619 - log2FC = 2.3; uvrA - KP13_00395 - log2FC = 2.3; uvrB - KP13_03032 - log2FC = 0.2; uvrC - KP13_01656 - log2FC = 1.2) and high iron condition (recN - KP13_02431 - log2FC = 0.2) (see Additional file 2). These observations were confirmed by the detection of co-expressed gene modules related to DNA damage repair processes, which are presented in the last section of the results.
Also, we identified that gene expression level of tisB (KP13_31738, Additional file 2), which codes a protein part of the SOS-response regulon controlled by LexA, was increased during PB pressure with high calcium/iron condition (log2FC>1.4) (Additional file 2). In E. coli, it has been shown that this protein is toxic since its overproduction cause slow or stop bacterial growth, probably allowing DNA repair before cells continue to grow .
With regard to the defense mechanisms against acidic pH conditions, we identified a gene (KP13_00169) encoding for a chloride channel/voltage gated protein up-regulated in PB with all abiotic changes (log2FC>1.4), but having considerable higher overexpression in acid pH condition (log2FC = 3.5) (Additional file 2). This protein belongs to the chloride channel (TC 2.A.49) family and probably acts as an electrical shunt for an outwardly-directed proton pump that is linked to amino acid decarboxylation, as part of the acid resistance response .
On the other hand, we also identified PB-dependent overexpression of a few genes encoding for proteins that are activated in response to cold shock, such as deaD (KP13_01093, cold-shock DEAD box protein A with highest log2FC>2.9 in PB with acid pH condition), cspA (KP13_00250, cold shock protein CspA with highest log2FC>4.8 in PB with acid pH condition) (Additional file 2).
We performed quantitative real-time PCR (qRT-PCR) to independently validate the expression of some of the studied genes (phoP, pmrB and arcB). These genes were chosen based on their involvement in the response to PB exposure, as discussed, for instance, with arcB. With few exceptions, the three genes displayed similar expression trends in the various experimental conditions tested, thus validating the digital gene expression measurements provided by the RNA-seq experiments (Additional file 5).
Environmental stimuli elicit differential responses in the expression of regulatory systems
On the other hand, differential activation of several regulators in specific abiotic stimuli were detected in our experiments, including the phoP-phoQ system which was found more activated in low pH and high calcium/iron conditions (Fig. 2; Additional file 5); zraR-zraS is a zinc-responsive TCRS involved in tolerance to metals and its differential activation occurred mostly in high calcium and iron conditions. This system has been regarded as important for pathogenic Salmonella during infection in pig . basR (or pmrA) and basS (or pmrB), which have been implicated in PB resistance, appear overexpressed in low pH and high iron conditions (Fig. 2), in line with previous observations . In high iron condition, however, there were opposite responses comparing the RNA-seq data with the qRT-PCR, where in the latter pmrB appeared repressed. cpxA, part of the CpxA-CpxR TCRS, appeared slightly overexpressed in acid pH and high iron/calcium concentrations (Fig. 2). This system is involved in envelope stress response, adhesion, motility and pathogenesis, and modulates the expression of a large number of genes . Furthermore, cpxA-cpxR has also been implicated in oxidative stress tolerance in K. pneumoniae .
We have also identified a strong (log2FC≥2) activation in acid pH and high calcium/iron conditions of evgS, which form part of a TCRS along with evgA. EvgA-EvgS system has been reported as modulator of the expression of several genes including multidrug-related, mostly in response to osmolarity changes in the medium. Thus, the finding that the sensor component evgS was more activated in the low pH condition was not unexpected. The analysis of genes previously reported to be under regulation of this system in E. coli  allowed the identification of orthologs in K. pneumoniae Kp13 that display similar expression response in the low pH condition. For instance, we identified in Kp13 an ortholog of yfdX (55 % identity to the E. coli sequence, BLASTP; Kp13 locus KP13_31877), that is strongly activated in low pH condition (log2FC = 2.0), in a similar fashion to the cognate gene in E. coli . These collateral results are of value to better explore the uncharacterized fraction of K. pneumoniae genes and represent a complementary result attained by this transcriptomic study in better delineating conserved regulatory responses in clinically important bacteria.
Co-expression network analysis reveals functional signatures related to polymyxin B, iron limitation and acidic pH
Module brown, also related to low pH, was enriched in terms related to ion transmembrane transport (BH-adjusted p-value = 0.2), which includes genes amtB (encoding for an ammonia channel) and the chloride channel KP13_00169 (previously discussed), mechanisms which would allow, respectively, for the entrance of NH3 (that can react with H+ yielding NH4 + and increasing pH) and the pumping out of H+ species . Indeed, amtB expression is higher when in low pH condition (log2FC = 1.1). In this module was also found arnB, arnC, arnD, arnE, arnF and arnT genes, encoding for proteins that participate in the addition of 4-amino-4-deoxy-L-arabinose to lipid A, contributing to a more positive charge on LPS, which in turn decreases PB affinity to the membrane . ArnA, encoded by gene arnA (found in the pink module), is also involved in this process, and was strongly up-regulated in acid pH (log2FC = 3.1; low pH vs. PB condition). The arn genes showed highest overexpression in PB with low pH. Genes related to stress and cold shock response (cspA, deaD), previously discussed, were found up-regulated in acid pH and also belong to this module.
The yellow module, having 350 associated genes, presented negative correlation to the low pH condition (Pearson’s r = −0.84; p-value = 7 × 10−4). Overall, it appears that this module comprises genes mostly related to respiration and oxidative cell processes. Functional categories identified in this module included aerobic respiration (BH-adjusted p-value = 3.7 × 10−5; genes sdhA, sdhB, sdhC, sdhD encoding for succinate DH). These genes were previously discussed in the context of PB presence, appearing slightly activated in PB exposure. When considering acid stress, however, these genes were strongly down-regulated, with log2FC values between −2.9 (sdhC) and −3.5 (shdA) (comparing low pH to PB). Generation of precursor metabolites and energy (BH-adjusted p-value = 2.1 × 10−5) was also a significantly enriched term, and genes encoding for ubiquinol oxidase (cyoAB), components of 2-oxoglutarate dehydrogenase complex (sucAB), succinyl-CoA ligase (sucC) could be related to these metabolic processes. All of these genes were found down-regulated when in acid pH. Also present were small heat shock proteins coded in the ibpA/ibpB genes as well as universal stress proteins A/B (uspA/B). uspA is reported to be highly transcribed in response to growth stasis, starvation or stress agents . Down-regulation of uspA (log2FC = −1.8) and uspB (log2FC = −3.5) in acid pH were identified in our analyses. The combined evidence argues for a diminished oxidative activity combined with stress response factors when in acid pH.
Another module of interest was the tan, which contains 76 genes with associated functional terms including response to antibiotic (BH-adjusted p-value = 0.002), DNA repair (BH-adjusted p-value = 0.08), DNA topological change (BH-adjusted p-value = 0.02). The involved genes include uvrB, dnaB, gyrA, parC, macA. Also, beta-lactamases bla SHV-12 and bla SHV-110 and genes encoding for fimbrial protein (KP13_04093), chaperone genes (hscAB) were found in this module. Although not directly related to PB resistance, it is interesting to note the activation of bla genes in PB with acid pH. While module tan correlated strongly only to the acid condition (Fig. 4), the finding that it groups resistance, virulence and DNA repair genes warrants further study on the common regulatory mechanisms that might control their expression, as disruption of these genes would affect many important correlates of infection and resistance.
Module darkgreen, which groups 38 genes, comprised several GO terms related to translation (adjusted p-value = 2.2 × 10−13) and peptide metabolic processses (adjusted p-value = 1.4 × 10−9); genes involved in this module included elongation factor (fusA, KP13_00742), several ribosomal subunit proteins (rplA, rplC, rplD, rplX, rplR, rplE, rplO, rplI, rplN, rplL, rpsA) and a protein translocase (secA, KP13_01897). That these ribosomal genes group into the same module is consistent with previous reports indicating their cotranscription in other bacteria , thus supporting the co-expression based analytical approach that we undertook. The negative correlation of this module with high iron and PB exposure hints an effect that protein synthesis might be inhibited at these conditions. Most interestingly is the possibility that polymyxin B could be exherting an inhibitory activity within the cell interior and effectively blocking these cellular processes, an observation that is in line with literature evidence for other CAMPs [61, 62]. This is also in agreement with our data that PB not only disrupts membrane functions, but also targets several central cellular processes such as carbon metabolic pathways (previously discussed) and provides further evidence to the idea that membrane disruption per se does not justify the killing of bacterial cells by these drugs, which probably involves additional intracellular targets.
It is worth noting that the genes coding for capsule biosynthesis (cps Kp13, which we previously studied in ) did not cluster in any single module. This could be because capsule production involves diverse biological processes, ranging from oligosaccharide chain biosynthesis, the polymerization of this chain and transport through the inner to the external membrane . As such, each of these steps should be regulated independently, hindering the clustering of these genes based in co-expression patterns alone. However, there are several evidences that capsule overproduction is also correlated with resistance to polymyxins [64, 65]. This phenomenon was observed in our RNA-seq experiments, but was only consistent for almost the whole cps Kp13 gene cluster in PB with low pH condition (Additional file 2 and Fig. 4).
The study of co-expressed modules and their association to the different experimental conditions shed light on important aspects underlying the K. pneumoniae response to PB and to abiotic stresses. Our results are consistent with the proposition that multiple metabolic aspects are affected in the course of PB exposure, and bacteria display a coordinated response which involves shifts in oxidative mechanisms, respiratory control, amino acid usage and transformation as well as regulatory components.
The present report provides a comprehensive gene expression dataset related to polymyxin B treatment under diverse abiotic stresses using K. pneumoniae Kp13 as model. We emphasize the functional transcriptomic aspects which are governed by PB treatment and the regulatory responses that may relate to exposure and resistance to this antibiotic, pinpointing intracellular targets that appear linked to the response elicited by PB exposure, such as the ArcA-ArcB system, previously unlinked to polymyxin response. Here, we highlight that polymyxin B exposure induced a metabolic shift into fermentative growth in Klebsiella pneumoniae Kp13.
Since PB resistance is an increasing trend in MDR bacteria, the study of its emergence by combination of expression profiling data from diverse bacteria, expanding on the present dataset, may aid on the search for common traits of resistance mechanisms linked to PB. A better understanding of these mechanisms will support the search of alternative or combination therapies to control these hard-to-treat infections.
PIPR is a recipient of a PhD fellowship from CAPES. We thank Eliane C. Vespero and Marsileni Pelisson for providing the Kp13 strain. We also thank Dr. Gabriela Bindea for preparing custom GO mapping files for use in CloudGO.
About this supplement
This article has been published as part of BMC Genomics Volume 17 Supplement 8: Selected articles from the Sixth International Conference of the Iberoamerican Society for Bioinformatics on Bioinformatics and Computational Biology for Innovative Genomics. The full contents of the supplement are available online at https://bmcgenomics.biomedcentral.com/articles/supplements/volume-17-supplement-8.
Fundação de Amparo à Pesquisa do Estado do Rio de Janeiro (FAPERJ) and Fundação de Amparo à Pesquisa do Estado de São Paulo (FAPESP) supported this work through grants awarded to MFN (process no. E-26/110.315/2014) and ACG (process no. 2010/12891-9), respectively. The publication charges for this article were funded by Coordenação de Aperfeiçoamento de Pessoal de Nível Superior (CAPES) (process no. 23038.010041/2013-13) through grant awarded to ATRV. The funders of the study had no role in the design, data analysis, interpretation or writing of this report.
Availability of data and materials
Expression values (fold-changes relative to the PB condition) are available in Additional file 2, and raw data is available from the corresponding authors upon request.
ACG, MFN conceived and coordinated the project. EF, GLS, GRQS, MGFC, PIPR, TC performed data analysis and wrote bioinformatics scripts. GB, RG, WMBSM carried out experimental procedures. ACG supervised experimental procedures. ATRV supervised the sequencing experiments. MFN and PIPR drafted the paper with input from the other authors. All authors read and approved the final manuscript.
The authors declare that they have no competing interests.
Consent for publication
Ethics approval and consent to participate
Open AccessThis article is distributed under the terms of the Creative Commons Attribution 4.0 International License (http://creativecommons.org/licenses/by/4.0/), which permits unrestricted use, distribution, and reproduction in any medium, provided you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated.
- Brown ED, Wright GD. Antibacterial drug discovery in the resistance era. Nature. 2016;529:336–43.View ArticlePubMedGoogle Scholar
- Davies J, Davies D. Origins and evolution of antibiotic resistance. Microbiol Mol Biol Rev. 2010;74:417–33.View ArticlePubMedPubMed CentralGoogle Scholar
- Zavascki AP, Goldani LZ, Li J, Nation RL. Polymyxin B for the treatment of multidrug-resistant pathogens: a critical review. J Antimicrob Chemother. 2007;60:1206–15.View ArticlePubMedGoogle Scholar
- Li J, Nation RL, Turnidge JD, Milne RW, Coulthard K, Rayner CR, Paterson DL. Colistin: the re-emerging antibiotic for multidrug-resistant Gram-negative bacterial infections. Lancet Infect Dis. 2006;6:589–601.View ArticlePubMedGoogle Scholar
- Deris ZZ, Akter J, Sivanesan S, Roberts KD, Thompson PE, Nation RL, Li J, Velkov T. A secondary mode of action of polymyxins against Gram-negative bacteria involves the inhibition of NADH-quinone oxidoreductase activity. The Journal of Antibiotics. 2014;67:147–51.View ArticlePubMedGoogle Scholar
- Lee J, Patel G, Huprikar S, Calfee DP, Jenkins SG. Decreased susceptibility to polymyxin B during treatment for carbapenem-resistant Klebsiella pneumoniae infection. J Clin Microbiol. 2009;47:1611–2.View ArticlePubMedPubMed CentralGoogle Scholar
- Gales AC, Jones RN, Sader HS. Contemporary activity of colistin and polymyxin B against a worldwide collection of Gram-negative pathogens: results from the SENTRY Antimicrobial Surveillance Program (2006–09). J Antimicrob Chemother. 2011;66:2070–4.View ArticlePubMedGoogle Scholar
- Ramos PIP, Picão RC, Almeida LGP, Lima NCB, Girardello R, Vivan ACP, Xavier DE, Barcellos FG, Pelisson M, Vespero EC, Médigue C, Vasconcelos ATR, Gales AC, Nicolás MF. Comparative analysis of the complete genome of KPC-2-producing Klebsiella pneumoniae Kp13 reveals remarkable genome plasticity and a wide repertoire of virulence and resistance mechanisms. BMC Genomics. 2014;15:54.View ArticlePubMedPubMed CentralGoogle Scholar
- Poirel L, Jayol A, Bontron S, Villegas M-V, Ozdamar M, Turkoglu S, Nordmann P. The mgrB gene as a key target for acquired resistance to colistin in Klebsiella pneumoniae. J Antimicrob Chemother. 2015;70:75–80.View ArticlePubMedGoogle Scholar
- Cannatelli A, Giani T, D’Andrea MM, Di Pilato V, Arena F, Conte V, Tryfinopoulou K, Vatopoulos A, Rossolini GM. MgrB inactivation is a common mechanism of colistin resistance in KPC-producing Klebsiella pneumoniae of clinical origin. Antimicrob Agents Chemother. 2014;58:5696–703.View ArticlePubMedPubMed CentralGoogle Scholar
- Yu Z, Qin W, Lin J, Fang S, Qiu J. Antibacterial mechanisms of polymyxin and bacterial resistance. Biomed Res Int. 2015;2015:1–11.Google Scholar
- Jayol A, Nordmann P, Brink A, Poirel L. Heteroresistance to colistin in Klebsiella pneumoniae associated with alterations in the PhoPQ regulatory system. Antimicrob Agents Chemother. 2015;59:2780–4.View ArticlePubMedPubMed CentralGoogle Scholar
- Adams MD, Nickel GC, Bajaksouzian S, Lavender H, Murthy AR, Jacobs MR, Bonomo RA. Resistance to colistin in Acinetobacter baumannii associated with mutations in the PmrAB two-component system. Antimicrob Agents Chemother. 2009;53:3628–34.View ArticlePubMedPubMed CentralGoogle Scholar
- Ramos PIP, Picão RC, Vespero EC, Pelisson M, Zuleta LFG, Almeida LGP, Gerber AL, Vasconcelos ATR, Gales AC, Nicolás MF. Pyrosequencing-based analysis reveals a novel capsular gene cluster in a KPC-producing Klebsiella pneumoniae clinical isolate identified in Brazil. BMC Microbiol. 2012;12:173.View ArticlePubMedPubMed CentralGoogle Scholar
- The European Committee on Antimicrobial Susceptibility Testing. Breakpoint tables for interpretation of MICs and zone diameters. Version 3.0, January 2013Google Scholar
- Andrews SC. FastQC: a quality control tool for high throughput sequence data. 2010 [http://www.bioinformatics.babraham.ac.uk/projects/fastqc/].Google Scholar
- Jiang H, Lei R, Ding S-W, Zhu S. Skewer: a fast and accurate adapter trimmer for next-generation sequencing paired-end reads. BMC Bioinformatics. 2014;15:182.View ArticlePubMedPubMed CentralGoogle Scholar
- Peixoto L, Risso D, Poplawski SG, Wimmer ME, Speed P, Wood MA, Abel T. How data analysis affects power reproducibility and biological insight of RNA-seq studies in complex datasets. Nucleic Acids Res. 2015;43:1–11.View ArticleGoogle Scholar
- Institut Pasteur MLST database [http://bigsdb.web.pasteur.fr/klebsiella/klebsiella.html]
- Risso D, Ngai J, Speed TP, Dudoit S. Normalization of RNA-seq data using factor analysis of control genes or samples. Nat Biotechnol. 2014;32:896–902.View ArticlePubMedPubMed CentralGoogle Scholar
- Kim D, Pertea G, Trapnell C, Pimentel H, Kelley R, Salzberg SL. TopHat2: accurate alignment of transcriptomes in the presence of insertions, deletions and gene fusions. Genome Biol. 2013;14:R36.View ArticlePubMedPubMed CentralGoogle Scholar
- Li H, Handsaker B, Wysoker A, Fennell T, Ruan J, Homer N, Marth G, Abecasis G, Durbin R. The Sequence Alignment/Map format and SAMtools. Bioinformatics. 2009;25:2078–9.View ArticlePubMedPubMed CentralGoogle Scholar
- Anders S, Pyl PT, Huber W. HTSeq–a Python framework to work with high-throughput sequencing data. Bioinformatics. 2015;31:166–9.View ArticlePubMedGoogle Scholar
- Robinson MD, McCarthy DJ, Smyth GK. edgeR: a Bioconductor package for differential expression analysis of digital gene expression data. Bioinformatics (Oxford, England). 2010;26:139–40.View ArticleGoogle Scholar
- Langfelder P, Horvath S. WGCNA: an R package for weighted correlation network analysis. BMC Bioinformatics. 2008;9:559.View ArticlePubMedPubMed CentralGoogle Scholar
- Langfelder P, Horvath S. Fast R functions for robust correlations and hierarchical clustering. J Stat Softw. 2012;46Google Scholar
- Bindea G, Mlecnik B, Hackl H, Charoentong P, Tosolini M, Kirilovsky A, Fridman W-H, Pages F, Trajanoski Z, Galon J. ClueGO: a Cytoscape plug-in to decipher functionally grouped gene ontology and pathway annotation networks. Bioinformatics. 2009;25:1091–3.View ArticlePubMedPubMed CentralGoogle Scholar
- Shannon P, Markiel A, Ozier O, Baliga NS, Wang JT, Ramage D, Amin N, Schwikowski B, Ideker T. Cytoscape: a software environment for integrated models of biomolecular interaction networks. Genome Res. 2003;13:2498–504.View ArticlePubMedPubMed CentralGoogle Scholar
- Bruchmann S, Muthukumarasamy U, Pohl S, Preusse M, Bielecka A, Nicolai T, Hamann I, Hillert R, Kola A, Gastmeier P, Eckweiler D, Häussler S. Deep transcriptome profiling of clinical K lebsiella pneumoniae isolates reveals strain and sequence type-specific adaptation. Environ Microbiol. 2015;17:4690–710.View ArticlePubMedGoogle Scholar
- Zhou L, Lei X-H, Bochner BR, Wanner BL. Phenotype microArray analysis of Escherichia coli K-12 mutants with deletions of all two-component systems. J Bacteriol. 2003;185:4956–72.View ArticlePubMedPubMed CentralGoogle Scholar
- Shimizu K. Metabolic regulation of a bacterial cell system with emphasis on Escherichia coli metabolism. ISRN Biochemistry. 2013;2013:1–47.View ArticleGoogle Scholar
- Ramirez-Santos J. Conserved regulatory elements of the promoter sequence of the gene rpoH of enteric bacteria. Nucleic Acids Res. 2001;29:380–6.View ArticlePubMedPubMed CentralGoogle Scholar
- Deris ZZ, Swarbrick JD, Roberts KD, Azad MAK, Akter J, Horne AS, Nation RL, Rogers KL, Thompson PE, Velkov T, Li J. Probing the penetration of antimicrobial polymyxin lipopeptides into gram-negative bacteria. Bioconjug Chem. 2014;25:750–60.View ArticlePubMedPubMed CentralGoogle Scholar
- Unden G, Bongaerts J. Alternative respiratory pathways of Escherichia coli: energetics and transcriptional regulation in response to electron acceptors. Biochim Biophys Acta. 1997;1320:217–34.View ArticlePubMedGoogle Scholar
- Spiro S, Guest JR. Adaptive responses to oxygen limitation in Escherichia coli. Trends Biochem Sci. 1991;16:310–4.View ArticlePubMedGoogle Scholar
- Sampson TR, Liu X, Schroeder MR, Kraft CS, Burd EM, Weiss DS. Rapid killing of Acinetobacter baumannii by polymyxins is mediated by a hydroxyl radical death pathway. Antimicrob Agents Chemother. 2012;56:5642–9.View ArticlePubMedPubMed CentralGoogle Scholar
- Tran QH, Bongaerts J, Vlad D, Unden G. Requirement for the proton-pumping NADH dehydrogenase I of Escherichia coli in respiration of NADH to fumarate and its bioenergetic implications. Eur J Biochem. 1997;244:155–60.View ArticlePubMedGoogle Scholar
- Belenky P, Ye JD, Porter CBM, Cohen NR, Lobritz MA, Ferrante T, Jain S, Korry BJ, Schwarz EG, Walker GC, Collins JJ. Bactericidal antibiotics induce toxic metabolic perturbations that lead to cellular damage. Cell Rep. 2015;13:968–80.View ArticlePubMedPubMed CentralGoogle Scholar
- Imlay JA, Chin SM, Linn S. Toxic DNA damage by hydrogen peroxide through the Fenton reaction in vivo and in vitro. Science (New York, NY). 1988;240:640–2.View ArticleGoogle Scholar
- Vatansever F, de Melo WCMA, Avci P, Vecchio D, Sadasivam M, Gupta A, Chandran R, Karimi M, Parizotto NA, Yin R, Tegos GP, Hamblin MR. Antimicrobial strategies centered around reactive oxygen species – bactericidal antibiotics, photodynamic therapy, and beyond. FEMS Microbiol Rev. 2013;37:955–89.View ArticlePubMedPubMed CentralGoogle Scholar
- Broxton CN, Culotta VC. SOD enzymes and microbial pathogens: surviving the oxidative storm of infection. PLoS Pathog. 2016;12:e1005295.View ArticlePubMedPubMed CentralGoogle Scholar
- Betteridge DJ. What is oxidative stress? Metab Clin Exp. 2000;49(2 Suppl 1):3–8.View ArticlePubMedGoogle Scholar
- Nair S, Finkel SE. Dps protects cells against multiple stresses during stationary phase. J Bacteriol. 2004;186:4192–8.View ArticlePubMedPubMed CentralGoogle Scholar
- Fozo EM, Kawano M, Fontaine F, Kaya Y, Mendieta KS, Jones KL, Ocampo A, Rudd KE, Storz G. Repression of small toxic protein synthesis by the Sib and OhsC small RNAs. Mol Microbiol. 2008;70:1076–93.View ArticlePubMedPubMed CentralGoogle Scholar
- Iyer R, Iverson TM, Accardi A, Miller C. A biological role for prokaryotic ClC chloride channels. Nature. 2002;419:715–8.View ArticlePubMedGoogle Scholar
- Lamarche MG, Wanner BL, Crépin S, Harel J. The phosphate regulon and bacterial virulence: a regulatory network connecting phosphate homeostasis and pathogenesis. FEMS Microbiol Rev. 2008;32:461–73.View ArticlePubMedGoogle Scholar
- Goh E-B, Bledsoe PJ, Chen L-L, Gyaneshwar P, Stewart V, Igo MM. Hierarchical control of anaerobic gene expression in Escherichia coli K-12: the nitrate-responsive NarX-NarL regulatory system represses synthesis of the fumarate-responsive DcuS-DcuR regulatory system. J Bacteriol. 2005;187:4890–9.View ArticlePubMedPubMed CentralGoogle Scholar
- Chen Y-T, Liao T-L, Wu K-M, Lauderdale T-L, Yan J-J, Huang I-W, Lu M-C, Lai Y-C, Liu Y-M, Shu H-Y, Wang J-T, Su I-J, Tsai S-F. Genomic diversity of citrate fermentation in Klebsiella pneumoniae. BMC Microbiol. 2009;9:168.View ArticlePubMedPubMed CentralGoogle Scholar
- Albright LM, Huala E, Ausubel FM. Prokaryotic signal transduction mediated by sensor and regulator protein pairs. Annu Rev Genet. 1989;23:311–36.View ArticlePubMedGoogle Scholar
- Huang Y, Leming CL, Suyemoto M, Altier C. Genome-wide screen of salmonella genes expressed during infection in pigs, using in vivo expression technology. Appl Environ Microbiol. 2007;73:7522–30.View ArticlePubMedPubMed CentralGoogle Scholar
- Falagas ME, Rafailidis PI, Matthaiou DK. Resistance to polymyxins: mechanisms, frequency and treatment options. Drug Resist Updat. 2010;13:132–8.View ArticlePubMedGoogle Scholar
- Srinivasan VB, Vaidyanathan V, Mondal A, Rajamohan G. Role of the Two component signal transduction system CpxAR in conferring cefepime and chloramphenicol resistance in klebsiella pneumoniae NTUH-K2044. PLoS One. 2012;7:e33777.View ArticlePubMedPubMed CentralGoogle Scholar
- Nishino K, Inazumi Y, Yamaguchi A. Global analysis of genes regulated by EvgA of the two-component regulatory system in Escherichia coli. J Bacteriol. 2003;185:2667–72.View ArticlePubMedPubMed CentralGoogle Scholar
- Krulwich TA, Sachs G, Padan E. Molecular aspects of bacterial pH sensing and homeostasis. Nat Rev Microbiol. 2011;9:330–43.View ArticlePubMedPubMed CentralGoogle Scholar
- Koponen J, Laakso K, Koskenniemi K, Kankainen M, Savijoki K, Nyman TA, de Vos WM, Tynkkynen S, Kalkkinen N, Varmanen P. Effect of acid stress on protein expression and phosphorylation in Lactobacillus rhamnosus GG. J Proteome. 2012;75:1357–74.View ArticleGoogle Scholar
- Lund P, Tramonti A, De Biase D. Coping with low pH: molecular strategies in neutralophilic bacteria. FEMS Microbiol Rev. 2014;38:1091–125.View ArticlePubMedGoogle Scholar
- Padilla E, Llobet E, Doménech-Sánchez A, Martínez-Martínez L, Bengoechea JA, Albertí S. Klebsiella pneumoniae AcrAB efflux pump contributes to antimicrobial resistance and virulence. Antimicrob Agents Chemother. 2010;54:177–83.View ArticlePubMedGoogle Scholar
- Gatzeva-Topalova PZ, May AP, Sousa MC. Structure and mechanism of ArnA: conformational change implies ordered dehydrogenase mechanism in key enzyme for polymyxin resistance. Structure. 2005;13:929–42.View ArticlePubMedPubMed CentralGoogle Scholar
- Sousa MC, McKay DB. Structure of the universal stress protein of Haemophilus influenzae. Structure. 2001;9:1135–41.View ArticlePubMedGoogle Scholar
- Shaik YB, Grogan S, Davey M, Sebastian S, Goswami S, Szmigielski B, Genco CA. Expression of the iron-activated nspA and secY genes in Neisseria meningitidis group B by fur-dependent and -independent mechanisms. J Bacteriol. 2007;189:663–9.View ArticlePubMedGoogle Scholar
- Spindler EC, Hale JDF, Giddings TH, Hancock REW, Gill RT. Deciphering the mode of action of the synthetic antimicrobial peptide Bac8c. Antimicrob Agents Chemother. 2011;55:1706–16.View ArticlePubMedPubMed CentralGoogle Scholar
- Podda E, Benincasa M, Pacor S, Micali F, Mattiuzzo M, Gennaro R, Scocchi M. Dual mode of action of Bac7, a proline-rich antibacterial peptide. Biochim Biophys Acta Gen Subj. 2006;1760:1732–40.View ArticleGoogle Scholar
- Whitfield C. Biosynthesis and assembly of capsular polysaccharides in Escherichia coli. Annu Rev Biochem. 2006;75:39–68.View ArticlePubMedGoogle Scholar
- Campos MA, Vargas MA, Regueiro V, Llompart CM, Albertí S, Bengoechea JA. Capsule polysaccharide mediates bacterial resistance to antimicrobial peptides. Infect Immun. 2004;72:7107–14.View ArticlePubMedPubMed CentralGoogle Scholar
- Llobet E, Tomás JM, Bengoechea JA. Capsule polysaccharide is a bacterial decoy for antimicrobial peptides. Microbiology. 2008;154(Pt 12):3877–86.View ArticlePubMedGoogle Scholar
- Provasi Cardoso J, Cayô R, Girardello R, Gales AC. Diversity of mechanisms conferring resistance to β-lactams among OXA-23–producing Acinetobacter baumannii clones. Diagn Microbiol Infect Dis. 2016;85:90–7.View ArticleGoogle Scholar
- Krzywinski M, Schein J, Birol I, Connors J, Gascoyne R, Horsman D, Jones SJ, Marra MA. Circos: an information aesthetic for comparative genomics. Genome Res. 2009;19:1639–45.View ArticlePubMedPubMed CentralGoogle Scholar
- Bardou P, Mariette J, Escudié F, Djemiel C, Klopp C. jvenn: an interactive Venn diagram viewer. BMC Bioinformatics. 2014;15:293.View ArticlePubMedPubMed CentralGoogle Scholar