- Research article
- Open Access
Genome-wide dynamics of a bacterial response to antibiotics that target the cell envelope
BMC Genomics volume 12, Article number: 226 (2011)
A decline in the discovery of new antibacterial drugs, coupled with a persistent rise in the occurrence of drug-resistant bacteria, has highlighted antibiotics as a diminishing resource. The future development of new drugs with novel antibacterial activities requires a detailed understanding of adaptive responses to existing compounds. This study uses Streptomyces coelicolor A3(2) as a model system to determine the genome-wide transcriptional response following exposure to three antibiotics (vancomycin, moenomycin A and bacitracin) that target distinct stages of cell wall biosynthesis.
A generalised response to all three antibiotics was identified which involves activation of transcription of the cell envelope stress sigma factor σE, together with elements of the stringent response, and of the heat, osmotic and oxidative stress regulons. Attenuation of this system by deletion of genes encoding the osmotic stress sigma factor σB or the ppGpp synthetase RelA reduced resistance to both vancomycin and bacitracin. Many antibiotic-specific transcriptional changes were identified, representing cellular processes potentially important for tolerance to each antibiotic. Sensitivity studies using mutants constructed on the basis of the transcriptome profiling confirmed a role for several such genes in antibiotic resistance, validating the usefulness of the approach.
Antibiotic inhibition of bacterial cell wall biosynthesis induces both common and compound-specific transcriptional responses. Both can be exploited to increase antibiotic susceptibility. Regulatory networks known to govern responses to environmental and nutritional stresses are also at the core of the common antibiotic response, and likely help cells survive until any specific resistance mechanisms are fully functional.
The bacterial cell wall is a key target for antibiotic discovery; it is crucial for cell growth, and provides a physical protective barrier between the cell and its environment. Antibiotics that inhibit bacterial cell wall biosynthesis, such as penicillin and vancomycin, are extremely important in the clinical treatment of infectious diseases. Free living bacteria resident in soil and marine environments produce a vast array of chemically diverse biologically active molecules, and the actinomycete Streptomyces spp. in particular are a rich source of antibiotics [1, 2]. Antibiotics do not kill the organisms that produce them since they have co-evolved systems that make them resistant or tolerant to their effects, but it is when similar systems develop in the target pathogenic bacteria allowing them to survive antibiotic treatments that major problems arise. The increase in the number of cases of methicillin-resistant Staphylococcus aureus (MRSA), and the emergence of vancomycin-resistant MRSA in hospital-acquired infections are two such examples [3–6].
Understanding how antibiotics can fail to be effective against bacteria has the potential both to provide novel insights into their mode of action, and to suggest new targets for anti-infective therapy. In this study we use Streptomyces coelicolor A3(2) as a model system to rigorously analyse the dynamic transcriptional response to sub-lethal concentrations of three antibiotics (vancomycin, bacitracin and moenomycin A) that target distinct stages of cell wall biosynthesis. Streptomyces species produce about 70% of known antibiotics and are likely to be the ultimate source of most clinically relevant antibiotic resistance genes. Consequently they possess many genes involved in sensing and responding to extracellular antibiotics, and provide an excellent opportunity to understand the processes involved. S. coelicolor exhibits significant resistance to vancomycin (minimum inhibitory concentration (MIC) ~80 μg/ml), bacitracin (MIC ~45 μg/ml) and moenomycin ((MIC ~300 μg/ml), and learning how the presence of these compounds, and the damage to the cell wall that they cause, is communicated to the bacterial chromosome, and how the subsequent reprogramming of gene expression acts to counteract this damage, can provide new ideas for drug discovery or for prolonging the therapeutic usefulness of existing compounds. Indeed, a biochemical analysis of the VanRS two-component sensor system of S. coelicolor recently defined exactly how bacteria sense vancomycin, a mechanism that triggers resistance to the only antibiotic in widespread use for the treatment of MRSA .
The development of techniques for studying gene transcription on a genome-wide scale allows characterisation of the response to antibiotic exposure, and several recent studies have used DNA microarrays to analyse the effect of antimicrobial compounds on bacterial strains [8–13]. However, only a minority have focused on antibiotics active against the cell wall, and these have tended to use either lethal drug concentrations, or to analyse only one time point following treatment. In this work we have studied three time points after treatment, and use the antibiotics at concentrations well below their MIC. Sub-lethal concentrations of antibiotics have been reported to be optimal for data quality in this type of study , perhaps because bacteriocidal levels seem ultimately to result in the same mechanism of cellular death, thereby obscuring the link between the transcriptional response to the antibiotic and its mode of action .
Vancomycin, a glycopeptide antibiotic, inhibits the transpeptidase reaction that occurs as one of the final steps in peptidoglycan biosynthesis not by binding to the enzyme but to the enzyme substrate, the D-alanyl-D-alanine terminus of the pentapeptide of the peptidoglycan precursor lipid II . Moenomycin A is a member of the phosphoglycolipid family of antibiotics, the only natural products known to target directly the transglycosylase enzyme that catalyzes polymerization of the peptidoglycan sugar backbone prior to transpeptidase-mediated cross-linking . Bacitracin is a branched cyclic dodecylpeptide antibiotic that inhibits cycling of lipid II by binding to the lipid II carrier undecaprenol pyrophosphate [17, 18]. By characterising the transcriptional response over a period of 90 min following treatment with each of these antibiotics we identified expression changes common to two or more compounds, or unique to each antibiotic. These genes are likely to be involved in adaptation or resistance to the antibiotics, the former representing a generalised response to cell wall damage, and the latter the antibiotic-specific responses. Antibiotic susceptibility studies using mutant strains selected for analysis on the basis of the microarray results confirmed a role for several genes in contributing to antibiotic resistance, validating the utility of this approach.
Exposure of S. coelicolor to antibiotics that inhibit different stages of peptidoglycan biosynthesis induced changes in transcription in up to 27% of the genome
Transcriptome analysis identified 2094 genes whose expression was significantly different at the 1% probability level, and more than 1.5-fold up- or down-regulated, as a result of exposure of exponentially growing S. coelicolor cells to vancomycin, bacitracin or moenomycin A (Figure 1, Additional files 1, 2, 3 and 4). This corresponds to 27% of the 7,825 genes predicted to make up the genome , and includes 1162 genes that were repressed at one or more time points after antibiotic addition, and 1062 genes whose expression was up-regulated. 130 genes were common to both the repressed and induced lists. For each antibiotic treatment, the majority of the genes identified as induced or repressed were present in at least two of the three time points analysed (Figure 1A and 1C). Strikingly, only a minority (0.06%) of the 1147 genes (68/1147) changing in expression in response to bacitracin treatment responded uniquely to this compound; the vast majority were similarly affected by vancomycin or by all three drugs (see Figure 1B). In contrast, 42% (454/1079) of the genes significantly down-regulated following addition of vancomycin responded only to this compound, and 31% (235/744) were also uniquely induced. Moenomycin treatment notably down-regulated the expression of far fewer genes than either vancomycin or bacitracin, and 16% of these (30/184) responded uniquely. Moenomycin however up-regulated a similar number of genes to either vancomycin or bacitracin, approximately 60% (364/616) of which were also induced by one or more of the other antibiotics, while 40% (252/616) were uniquely induced. Interestingly, a core set of 243 genes were up-regulated in response to any one of the compounds and 118 were similarly down-regulated, representing a common response to the antibiotic treatments.
To determine the similarities and differences in the effects of each drug, the differentially expressed sets of genes illustrated in Figure 1 were subjected to gene ontology (GO) analysis (Additional files 5, 6, 7, 8, 9, 10 and 11), and compared with in-house curated sets of functionally related genes (Tables 1 and 2, Additional file 12). The latter approach was included since the computationally generated GOA annotation available for the S. coelicolor proteome (see Methods) has not been manually curated and is therefore incomplete. The results of these analyses are summarised in the following sections.
Common changes in expression of known stress-response genes following antibiotic treatment suggest a general environmental stress response system in S. coelicolor
i) The σE cell envelope stress response system is induced by all three drugs
σE is an extracytoplasmic function sigma factor required for normal cell envelope integrity in S. coelicolor, and its expression is induced by a wide-variety of cell wall damaging agents . Reassuringly, transcription of the four-gene operon containing the σE gene was significantly induced following treatment with each of the antibiotics (Table 1, Figure 2). In addition, similarity analysis indicated that the expression of 158 genes was closely correlated with that of σE (Pearson correlation >0.9; Additional file 13). This group is likely to include genes directly dependent on σE for their transcription, and indeed contains σhrdD (SCO3202) previously characterised to be dependent on σE for expression . Interestingly, expression of the only aminoacyl-tRNA synthetase gene (SC03397) to be up-regulated in response to any drug correlates with σE transcription (similarity coefficient = 0.93); of the remaining 27 tRNA synthetase genes annotated in the genome sequence, transcription of 21 was significantly repressed following antibiotic treatment. Mutation in SCO3397, encoding a lysyl-tRNA synthetase, markedly reduced the MIC towards both vancomycin and bacitracin (Table 3) indicating that it plays an important role in resistance to these compounds.
ii) Effects caused by exposure to cell wall antibiotics show significant similarity to those resulting from synthesis of the stringent factor ppGpp
The stringent factor ppGpp is an intracellular signalling molecule that is synthesised in S. coelicolor in response to amino acid limitation, and which co-ordinates an adaptive response to the changing nutritional conditions that includes down-regulation of processes associated with active growth, and up-regulation of systems that may help it survive the nutritional stress . Comparison of the lists of genes that were down-regulated in response to antibiotic treatment with the genes shown to be repressed following controlled induction of ppGpp synthesis  revealed a significant commonality for each of the compounds (Table 2). The 118 genes down-regulated in response to all three antibiotics were also significantly enriched for ppGpp-repressed genes, containing 25/189 (p-value < 1E-10; Additional file 12). These included 9 genes from conservons 5, 6 and 13. Conservons are conserved operons believed to encode membrane-associated signalling complexes , and more than half of the 51 cvn genes (from cvns 1, 2, 4, 5, 6, 10, 12 and 13) were repressed by bacitracin or vancomycin, while 11 (from cvns 5, 6, and 13) were repressed by moenomycin. In addition, ppGpp-induced genes were significantly over-represented in each of the vancomycin-, bacitracin-, and moenomycin-induced lists (Table 1), and in the 243 genes commonly induced by all of the antibiotics (Additional file 12).
To investigate the possibility that an increase in ppGpp synthesis helps mediate the observed responses to the antibiotic treatments, intracellular nucleotide pools were measured in cultures following exposure to vancomycin. No increase in ppGpp synthesis was observed at 15, 30 or 60 min following addition of vancomycin, but ATP and GTP pools were reduced on average by 23% and 35% after 15 min, respectively (data not shown). MIC tests indicated that a relA mutant strain that is incapable of ppGpp synthesis  is more susceptible to both bacitracin and vancomycin (Table 3).
iii) Antibiotic treatment induces changes in expression of many genes from the heat, osmotic and oxidative stress regulons
Recent studies by Bucca et al.  and Lee et al.  characterised the heat-shock regulon and σB-dependent salt-stress regulon in S. ceolicolor, identifying 119 and 87 genes, respectively. Comparison of the genes significantly differently expressed in response to antibiotic treatment in this study indicated a significant overlap between the sets of data (Tables 1 and 2). The 744 genes up-regulated in response to vancomycin contained 51 of the 119 genes reported by  to be up-regulated in response to heat shock, an over-representation with a significance p-value of <1 E-33. Almost half of the σB-dependent salt-stress regulon genes (39/87) were similarly present (p-value <1 E-33), while the 616 and 481 genes induced by moenomycin and bacitracin, respectively, contained 42/87 and 33/87 of the salt stress genes, and 43/119 and 27/119 of the heat shock genes, respectively. Indeed, of the 243 genes commonly up-regulated in response to all three antibiotics, 26 belonged to the heat shock regulon and 29 to the σB-dependent salt-stress regulon. A σB-deletion mutant showed a lowered MIC for both bacitracin and vancomycin when compared to the parent strain (Table 3), indicating an important role for σB in responding to cell envelope stress.
Genes from the well-characterised σR regulon responsible for disulphide bond homeostasis in response to oxidative stress  were significantly over-represented in both the vancomycin- and bacitracin-repressed gene lists, suggesting a decrease in σR-mediated transcriptional activity after antibiotic treatment.
iv) Vancomycin, bacitracin and moenomycin up-regulate ectoine biosynthesis and genes for gas vesicle production
Ectoine is a pyrimidine carboxylic acid derivative that serves as a compatible solute in several bacterial species, protecting against salt and also temperature stress [27–32]. Exogenous addition of ectoine and 5-hydroxyectoine to S. coelicolor cultures confers protection against salt and heat stress, and a cluster of genes (SCO1864-1867) homologous to the ectABCD gene cluster required for ectoine biosynthesis in other bacteria was identified in the S. coelicolor genome sequence . GO terms for ectoine biosynthesis and production were significantly over-represented in the lists of genes induced following exposure to each of the three antibiotics and the ectABCD homologues were up-regulated at all three time points after their addition (see Additional files 1,2 and 3, Additional file 14 Figure 1g). Genes required for gas vesicle biosynthesis were similarly induced, but only those from one of the two gvp gene loci, gvp1 (Table 1, Additional files 1, 2 and 3, Additional file 14 Figure 1h). A strong induction of gvp1 gene expression has also been reported following exposure to high concentrations of salt, and it has been proposed that gas vesicles play a role in the response to hyperosmotic stress .
In MIC tests, supplementation of the agar media with 1 mM ectoine provided no significant protection against the effects of either bacitracin or vancomycin (data not shown).
Changes in expression of genes with functions related to production and maintenance of the cell wall show antibiotic-specificity
i) Genes required for cell wall peptidoglycan biosynthesis
Bacterial cell wall peptidoglycan biosynthesis can be separated into reactions occurring intracellularly to generate the lipid II precursor subunits, and those taking place extracellularly to join the subunits together into a three dimensional structure [35–42] (Figure 3A). Transcription of many of the genes for the cytoplasmic steps required for precursor biosynthesis, particularly the mur operon (SCO2083-2089) and the femX gene (SCO3904), were significantly repressed following treatment with either bacitracin or vancomycin, but were unchanged following moenomycin addition (Table 2, Figure 3B). In contrast, genes encoding penicillin-binding proteins (PBPs) required for the extracellular reactions in peptidoglycan production were over-represented in the genes induced by all three antibiotics (Table 1, Figure 3B). The secreted PBPs SCO1875, SCO2897 and SCO5039 were significantly induced following treatment with each compound, while SCO2608 was induced by vancomycin and moenomycin but not by bacitracin. In MIC tests, deletion of SCO1875 markedly increased susceptibility to both bacitracin and vancomycin (Table 3) indicating that this PBP is essential for optimum resistance to these compounds, and suggesting that up-regulation of its transcription is an important response to the antibiotic treatments. Deletion of SCO2897 also increased susceptibility to both antibiotics (lowered the MIC), while mutation in SCO2608 and SCO5039 increased sensitivity to vancomycin only (Table 3). Vancomycin treatment alone down-regulated expression of SCO2090, SCO3901 and SCO4013, and was the only antibiotic to up-regulate transcription of SCO3156 and SCO3847. Consistent with this, deletion of SCO3847 slightly lowered the MIC for vancomycin but showed no effect with bacitracin (Table 3). The PBP encoded by SCO3157 was up-regulated by both vancomycin and bacitracin but not by moenomycin.
The previously characterised van gene cluster responsible for remodelling peptidoglcan precursor biosynthesis to produce subunits terminating in D-alanyl-D-lactate (rather than D-alanyl-D-alanine) was, as expected, induced strongly only by vancomycin, and not by either of the other two antibiotics [43, 44](Figure 3B). vanH, encoding D-lactate dehydrogenase, exhibited the largest change in expression observed for any gene: a 300-fold increase 90 min after vancomycin addition.
ii) Genes encoding cell wall remodelling enzymes and sortase substrates
In addition to peptidoglycan biosynthesis, hydrolysis of the cell wall is an important process for growth allowing remodelling of cell shape to enable expansion and differentiation. Haiser et al.  identified 57 putative cell wall hydrolase genes in the S. coelicolor genome in an in silico search, and confirmed that four could efficiently degrade S. coelicolor cell walls using biochemical approaches. In this study, genes encoding cell wall hydrolases were significantly over-represented in the lists of genes up-regulated in response to bacitracin or moenomycin, but not vancomycin (see Table 1). In total, expression of 23 of the 57 genes were significantly affected following treatment with one or more of the cell wall antibiotics, only a minority (four) of these however responded in a similar way to all three compounds (Additional file 14 Figure 1e). Seven genes were up-regulated only in response to moenomycin, and four of these were actually repressed following treatment with at least one of the other two antibiotics. In total 13 of the 23 genes were down-regulated in response to either bacitracin or vancomycin, while moenomycin treatment did not negatively affect transcription of any. Interestingly, of the seven members of the resuscitation-promoting factor (Rpf) subgroup of hydrolases, five were significantly induced in response to moenomycin. Only rpfA (SCO3097) was induced by all three antibiotics, with expression of the SCO3098 homologue induced by bacitracin and moenomycin, but significantly repressed following vancomycin treatment (Additional file 14 Figure 1e). rpfA and SCO3098 encode proteins with putative LysM peptidoglycan-binding domains, and RpfA has been shown to hydrolyse S. coelicolor peptidoglycan . Of the three other rpf homologues induced by meonomycin, SCO5029 was repressed following bacitracin and vancomycin treatment, while SCO0974 expression was down-regulated by bacitracin. Transcription of SCO3150 was not significantly affected by either vancomycin or bacitracin. The Nlp60-family cell wall hydrolase encoded by SCO4108 was induced in response to both bacitracin and moenomycin, but repressed by vancomycin. Vancomycin and bacitracin strongly induced both SCO5487 (amidase family) and SCO5660 (carboxypeptidase family), while moenomycin up-regulated only the former of these two genes.
Septum synthesis and localisation is a specialised aspect of cell wall biosynthesis important for cell division and sporulation in S. coelicolor. SsgA-like proteins are intimately involved in developmental cell division in streptomycetes, and although their precise function has yet to be fully characterised they are believed to control the fate of peptidoglycan during sporulation [46, 47]. SsgA-like proteins were significantly over-represented in the vancomycin- and bacitracin-induced gene lists (see Table 1), and ssgR (SCO3925), the regulator activating transcription of ssgA (SCO3926), was induced 2-to 6-fold at all times following treatment with vancomycin (Additional file 14 Figure 1o). Expression of ssgD (SCO6722) was significantly up-regulated by all three antibiotics, while ssgE (SCO3158) was induced by vancomycin and moenomycin only.
Sortase enzymes are transpeptidases that covalently link proteins to peptidoglycan in many Gram-positive bacteria . Substrate proteins are recognisable from the presence of amino acid sequence motifs, most notably by a conserved C-terminal LP(A)XTG sequence. Genes encoding 17 substrate proteins were identified in S. coelicolor, and presumably play a role in cell wall structure and function. For example, three are known to encode chaplin proteins that assemble as a hydrophobic coat during aerial mycelium formation when cultured on solid agar media . Genes encoding sortase substrates were significantly over-represented in the lists of genes repressed by each of the drugs (see Table 2), with six genes significantly altered in expression following treatment with one or more of the antibiotics (Additional file 14 Figure 1n). Four (SCO2492, SCO2682, SCO3176 and SCO5650) were down-regulated in comparison to the control experiment by all three compounds, while SCO2270 was repressed only by bacitracin and vancomycin. SCO1023, encoding a protein of unknown function, was induced in response to all three antibiotics. As expected, expression of the chaplin genes was not significantly altered since cells were grown in submerged liquid culture where aerial mycelium formation does not take place.
The biosynthesis of antibiotics is induced in response to cell wall specific drugs
In the gene ontology analysis (Additional files 5, 6,7, 8, 9, 10 and 11), the GO term for acyl carrier activity was prominent in the lists from the term-for-term and MGSA analysis of genes induced by vancomycin or bacitracin in the 90 min sample, with as many as 8 genes associated with this term present in the lists of up-regulated genes. These correspond to four genes each from the biosynthetic clusters responsible for producing the antibiotics CDA and undecylprodiginine (Red), and indicate that activation of antibiotic production is part of the cellular response to exposure to both vancomycin and bacitracin. Indeed, 39 of the 47 CDA cluster probe sets (corresponding to 32 of the 40 genes) were up-regulated after addition of vancomycin, along with all 23 Red cluster probe sets (22 of the 22 genes; Table 1, Additional files 14 Figures 1d and 1i). The corresponding figures for bacitracin treatment were 35 and 9 probe sets, respectively, while moenomycin induced 11 probe sets from the CDA cluster. Activation of CDA synthesis is therefore a common response to treatment with all of the antibiotics. The expression of the CDA and Red biosynthesis genes, and that of the only other antibiotic biosynthetic gene cluster identified in the genome, the actinorhodin (Act) cluster, is ultimately dependent on activation by regulatory proteins encoded by each cluster [51–53]. The changes in expression of these activator genes brought about by treatment with the cell wall antibiotics is summarised in Table 4, and indicates that vancomycin and bacitracin in particular had a strong inducing effect on transcription of all four activators. Although an increase in expression of the actinorhodin biosynthesis genes was not observed in the 90 min in liquid culture following addition of any of the antibiotics, vancomycin treatment induced production of the blue pigmented antibiotic in a paper disc assay using growth on a solid agar medium (Figure 4). This response was markedly reduced, but not abolished, in an M570 mutant defective in ppGpp synthesis, indicating that the signalling molecule is not essential for this response to vancomycin.
Sigma factor gene expression is a major response to exposure to cell wall antibiotics
S. coelicolor RNA polymerase can interact with any one of 64 different sigma factors, each specifying recognition of a different cognate consensus promoter sequence . Changes in the relative abundance of sigma factors within a cell can therefore significantly influence the pattern of gene transcription. Sigma factor genes were significantly over-represented in the lists of genes induced by each antibiotic (Table 1) indicating that alteration in sigma factor gene expression is an important response to treatment with antibiotics targeting cell wall biosynthesis. This is also reinforced by the prominence of genes known to be regulated by σE, σB and σR in the differentially expressed gene lists (see Tables 1 and 2), and by the observed increase in susceptibility of σE and σB mutant strains towards bacitracin and vancomycin (see Table 3). Eight sigma factor genes, including σE (SCO3356) and σB (SCO0600), were up-regulated following treatment with all three compounds. Of these, the genes encoding the putative extra-cytoplasmic function (ECF) sigma factors encoded by SCO4005, SCO4908 and SCO5147 were particularly strongly induced (5- to 7-fold by vancomycin; Additional file 14 Figure 1l). A mutant strain in which SCO4005 was disrupted by insertion of a transposon showed increased susceptibility to bacitracin, but was not significantly altered in its resistance to vancomycin (Table 3). However, a similar mutant in which SCO5147 had been disrupted showed no increase in susceptibility to either compound and was slightly more resistant to vancomycin than the parental strain (Table 3). Transcription of the alternative sigma factor σhrdD (SCO3202) known to be dependent on σE was also up-regulated by all three antibiotics, while the expression of σbldN (SCO3323), a sigma factor important for S. coelicolor differentiation , was repressed by all three. Vancomycin alone repressed transcription of four other sigma factor genes (SCO3068, SCO3626, SCO5216 (σR), and SCO5243); two of these, including σR, were up-regulated after treatment with at least one of the other two antibiotics. Consistent with this, genes from the σR-regulon were significantly over-represented in both the vancomycin-repressed and moenomycin-induced gene lists (see Tables 1 and 2).
The zinc-responsive zur regulon is induced exclusively by vancomycin
Zinc plays an indispensable role in cellular biochemistry as a catalytic or structural cofactor for a wide variety of metalloproteins. It can however be toxic if accumulated to excess, and its intracellular concentration is controlled within safe limits by homeostatic regulatory mechanisms. In bacteria this is largely achieved by tight control of import and export mechanisms via the zinc-responsive regulatory protein Zur [56, 57], and the Zur regulon in S. coelicolor has recently been characterised [58–60]. Interestingly in this work the majority of the genes in the regulon, including the znuABC operon (SCO2505-2507) encoding the high-affinity Zn2+ import system, were transiently up-regulated following addition of vancomycin but not by either of the other two antibiotics (Table 1, Additional file 14 Figure 1p). Expression of znuA (SCO2505) was induced 12-fold 30 min after vancomycin addition, decreasing to about 2-fold up-regulated by 90 min (see Additional file 3).
Genes required for branched chain amino acid and biotin biosynthesis are induced by moenomycin
GO terms associated with amino acid biosynthesis were over-represented in the genes induced in response to moenomycin treatment (see Additional files 5, 6, 7, 8, 9 and 10). Comparison of the significantly differently expressed genes with the groups of functionally related genes in Tables 1 and 2 supported this observation, indicating that moenomycin induced genes associated with branched chain amino acid (leucine, isoleucine and valine), lysine and methionine biosynthesis. Indeed, the ilvD (SCO3345), leuB (SCO5522) and leuC (SCO5553) genes for the biosynthesis of the branched chain amino acids leucine, valine and isoleucine were up to 7-fold up-regulated by moenomycin treatment (Additional file 14 Figure 1c). In addition, genes associated with the transport of amino acids were over-represented in both the moenomycin up-regulated genes, and those down-regulated by the antibiotic (Tables 1 and 2). Interestingly, the down-regulated genes included four from the putative branched-chain amino acid import operon SCO2008-2012 (Additional file 14 Figure 1a).
GO analysis of the 252 genes specifically up-regulated by moenomycin in Figure 1B indicated a significant over-representation not only of amino acid biosynthesis genes, but also of biotin biosynthesis. The bioB-D operon (SCO1244-1246) was approximately 2-fold up-regulated in the 90 min moenomycin-treated sample. Biotin is a cofactor implicated in catabolism of fatty acids and branched chain amino acids, and up-regulation of its production may indicate an increased use of branched chain amino acids as an energy source in the moenomycin treated cells. Indeed, expression of accD1 (SCO2776) encoding the biotin-dependent acetyl-CoA carboxylase enzyme, a key enzyme in branched chain amino acid degradation, was 5-fold up-regulated 90 min after moenomycin addition (see Additional file 3).
Genes encoding two-component signal transduction systems are specifically induced by bacitracin
GO analysis of the 39 genes whose transcription was induced only in response to bacitracin (see Figure 1B) indicated genes associated with two-component system signal transduction as the most over-represented process (see Additional file 11), with ten of the genes encoding response regulator or sensor kinase proteins. These included two clusters of three genes encoding a response regulator-kinase-kinase (SCO4596-4598), and a regulator-regulator-kinase (SCO3638, SCO3640 and SCO3641). This implies that they represent multi-component signal transduction systems that are important for coordinating a response to bacitracin (but not moenomycin or vancomycin). Other two-component system genes uniquely up-regulated in response to bacitracin were the SCO5683-5684 regulator-kinase pair, and the kinases encoded by SCO2359 and SCO7649. Mutant strains deleted for SCO3638-3641, SCO4596-4598 or SCO5683-5684 did not however show any significant change in their resistance to bacitracin (data not shown).
Identification of novel ABC transport systems important for bacitracin resistance
SCO3089-3090 and SCO3010-3011 encode putative ABC transport proteins that show a high degree of homology: SCO3089 is 75% identical to SCO3111 at the amino acid level, while SCO3090 shows 38% identity over its full-length to SCO3110 (Figure 5A). Transcription of all four genes was strikingly up-regulated 30-fold or higher in response to each of the three compounds (Figure 5B, Additional files 1, 2 and 3), suggesting an important role for these transporters in the response to the antibiotic treatments. To test this further, mutant strains carrying deletions in one or both of the transport systems were constructed and analysed for their susceptibility to vancomycin and bacitracin (Table 3 and Figure 5C). While the single and double mutant strains showed at most only a marginal increase in susceptibility to vancomycin, the double mutant strain H1004 exhibited markedly increased susceptibility to bacitracin, reducing the MIC to below 10 μg/ml compared to 40-45 μg/ml for the parental strain. Both single mutant strains also showed an increased susceptibility to bacitracin, but MIC values (15-20 μg/ml for H1003 (ΔSCO3089-3090) and 30-35 μg/ml for H1002 (ΔSCO3110-3111)) were not as low as for the double mutant.
Stress response genes contribute towards antibiotic resistance in S. coelicolor
Vancomycin, bacitracin and moenomycin target distinctly different stages of cell wall biosynthesis (see Figure 3A), yet exposure of S. coelicolor to each antibiotic resulted in changes in expression of many of the same genes (see Figure 1). Previous studies on the role of the ECF (e xtrac ytoplasmic f unction) sigma factor σE in maintaining cell wall homeostasis in S. coelicolor identified a key role in responding to antibiotics that damage the cell wall . This study reveals that the transcriptional response to vancomycin, moenomycin and bacitracin extends far beyond the involvement of σE, incorporating elements of the heat shock response, the stringent response, and the response to osmotic and oxidative stress. Whether this reflects different, condition-specific regulatory cascades influencing the same genes, or triggering of the same regulatory cascade by different stimuli is unclear. The failure of vancomycin to increase synthesis of the stringent factor ppGpp while modulating expression of many stringently controlled genes suggests the former, but simultaneous induction of transcription of σB and the σB regulon, and of σE and the σE-dependent gene SCO3202, by all three antibiotics would be indicative of the latter. A combination of the two appears likely. These stress responses were not necessary for antibiotic resistance, but were required for maximal resistance. A σB-deletion mutant strain and a mutant defective in ppGpp synthesis both exhibited increased susceptibility towards bacitracin and vancomycin. S. coelicolor carries an inducible system that specifically confers high-level (MIC ~80 μg/ml) resistance to vancomycin, and functions by reprogramming cell wall precursor biosynthesis [43, 44]. Although the van genes involved are among the most highly induced in response to vancomycin treatment at all time points tested (see Figure 3B, Additional file 1), it is unlikely that cells become instantly fully resistant to the effects of the antibiotic and expression of the stress response genes may facilitate survival during this transition period. A transient slowing in growth was observed following treatment with both vancomycin and bacitracin, consistent with the existence of such a "stress" phase (data not shown).
S. coelicolor is susceptible to bacitracin (MIC 40-45 μg/ml) but highly resistant to moenomycin (MIC >300 μg/ml). All of the mutants tested were also resistant to 300 μg/ml moenomycin. It is perhaps therefore not surprising that moenomycin repressed the transcription of many fewer genes than either bacitracin or vancomycin (see Figure 1A). Interestingly however, moenomycin activated the transcription of a comparable number of genes to both vancomycin and bacitracin, including genes from the heat shock regulon, the stringent response, and the osmotic and oxidative stress regulons. Treatment with 10 μg/ml moenomycin, although only a fraction of the MIC, clearly induced a physiological response, perhaps reflecting induction of mechanisms for cell wall homeostasis well below inhibitory concentrations of the antibiotic or even a signalling role for the compound (see below).
Identification of new genes with roles in antibiotic resistance
Genes whose transcription was activated in response to two or more of the compounds may be components of a general response to treatment with antibiotics inhibiting cell wall biosynthesis. Inactivation of two putative ABC transporters, SCO3089-3090 and SCO3110-3111, both markedly up-regulated after treatment with all three antibiotics, significantly increased susceptibility to bacitracin. The integral membrane proteins encoded by these two systems are similar to lysophospholipase L1 and transport system components associated with lipoprotein release (http://streptomyces.org.uk), possibly suggesting a role in remodelling the cell envelope. Interestingly, orthologous and syntenous systems also occur in the genomes of Streptomyces clavuligerus, Streptomyces avermitilis and Streptomyces griseus, suggesting conserved evolution and a common biological function. The ECF sigma factor encoded by SCO4005 also contributes to bacitracin resistance; identifying the SCO4005 regulon could provide new insights into how bacterial cells cope with the antibiotic.
Of the other genes inactivated and analysed for their effect on antibiotic resistance, only disruption of the PBP encoded by SCO1875 and of the putative lysyl-tRNA synthetase encoded by SCO3397 markedly increased susceptibility to both bacitracin and vancomycin. Transcription of both SCO1875 and SCO3397 is markedly up-regulated following treatment with each of the antibiotics and is therefore an important part of the antibiotic resistance response. Navratna et al.  reported that deletion of the gene encoding PBP4 in Staphylococcus aureus increased susceptibility to vancomycin and was associated with a reduction in cross-linked muropeptide. PBP5 in Enterococcus faecalis is required for resistance to cephalosporins . None of these three PBPs are essential for growth in their respective organisms, yet modulation of their activities appears to be an important component of the response to compounds that inhibit peptidoglycan biosynthesis. The response to vancomycin treatment included the up-regulation of transcription of eight genes encoding putative PBPs, including two only induced by this antibiotic. These observations, together with three PBPs genes uniquely repressed following vancomycin addition are discussed later. Of the 28 genes predicted to encode tRNA synthetases in the S. coelicolor genome, only SCO3397 was up-regulated following drug treatment, the majority (21) typically being down-regulated. In addition to their central role in protein biosynthesis, aminoacyl-tRNAs are involved in porphyrin biosynthesis, N-terminal protein modification, peptidoglycan biosynthesis and aminoacyl-phosphotidylglycerol biosynthesis . SCO3397 is required for maximal resistance to bacitracin and vancomycin in S. coelicolor, and homologous MprF proteins in Staphylococcus aureus and Listeria monocytogenes have also been implicated in resistance to antimicrobial peptides [64, 65]. Although it is not known whether SCO3397 possesses phospholipid lysylination activity similar to MprF, this protein class may represent a useful target for counteracting resistance to antibiotics active against the cell envelope.
Genes responding uniquely to individual drug treatments offer new insights into activity and resistance
While there were many similarities in the changes in gene expression following treatment with vancomycin, moenomycin or bacitracin, genes responding uniquely to only one of the antibiotics were also identified (see Figure 1B). Given the different activities of the three compounds, these results potentially provide new insights into the mode of action of, and mechanisms of resistance towards, the individual antibiotics.
In total, transcription of 235 genes was uniquely activated by vancomycin, while 454 were uniquely repressed (see Figure 1B). As expected, the genes uniquely up-regulated included the set of seven van genes previously characterised as encoding an inducible resistance system specific to this drug [43, 44]. Surprisingly, they also included 19/22 genes from the Zur regulon required for zinc homeostasis [58–60] that were transiently but strongly up-regulated in the 30 min sample. The implication is that vancomycin addition caused a decrease in intracellular zinc availability that was subsequently corrected via the Zur-dependent zinc starvation response. How vancomycin causes zinc limitation is not known, but could result from changes in the cell envelope that favour zinc sequestration or perhaps via chelation of zinc by the antibiotic itself (vancomycin interacts with Cu2+ ions in vitro[66, 67]). Transcription of five genes encoding PBPs and therefore with putative roles in extracellular peptidoglycan biosynthesis were uniquely significantly affected by vancomycin treatment (SCO3156 and SCO3847 were up-regulated while SCO2090, SCO3901 and SCO4013 were repressed). It is interesting to speculate that at least some of these changes may be associated with the reprogrammed peptidoglycan biosynthesis producing precursors terminating in D-alanyl-D-lactate and perhaps requiring PBPs to match the changes in substrates. SCO5034-5035 encoding a putative ABC transporter was 3.5-fold or more up-regulated at all time points in response to vancomycin, but not moenomycin or bacitracin (see Additional file 14, Figure 1b). The products of these genes are homologous to the transporters encoded by SCO3089-3090 and SCO3110-3111 (about 20% identity between the permease proteins, and 40% between the ATP-binding subunits) that are strongly induced in response to all three antibiotics. Functional redundancy associated with vancomycin-dependent expression may explain why the double mutant strain H1004 (ΔSCO3089-3090 + ΔSCO3110-3111) was sensitive to bacitracin but not to vancomycin. Interestingly, the putative regulatory gene SCO5036 immediately upstream of SCO5034-5035 was also induced by vancomycin, and was >15-fold up-regulated at all times following antibiotic treatment. From their chromosomal arrangement these genes appear to form a single transcription unit directed from a promoter upstream of SCO5036.
Genes exclusively induced by bacitracin are particularly interesting because none correspond to previously known bacitracin resistance genes in other organisms, such as bacA, bcrA, and bcrC[68–72]. This suggests that S. coelicolor may possess a novel mechanism for bacitracin resistance. Treatment with bacitracin activated the transcription of many two-component regulatory systems, including several whose strong induction was unique to this antibiotic. In fact, about 25% of the genes uniquely up-regulated in response to bacitracin treatment encoded response regulators or sensor histidine kinases. However, mutant strains carrying deletions in individual sensory systems showed no change in susceptibility to bacitracin, indicating, in contrast to vancomycin, the existence of multiple, functionally redundant sensory systems. The bcrA and bcrC homologues induced by bacitracin in this work were also similarly up-regulated by at least one of the other two antibiotics, and are therefore likely to be involved in a generalised antibiotic response. Transcription of the bacA homologue SCO7047 was not significantly affected by any of the antibiotics.
iii) Moenomycin A
Moenomycin A exposure was the only treatment that activated transcription of genes required for amino acid biosynthesis; vancomycin and bacitracin tended to repress their expression. Genes encoding enzymes with roles in branched chain amino acid biosynthesis were particularly strongly induced, while those required for their import were simultaneously repressed. While this may reflect inhibition by moenomycin of the uptake of these amino acids, it could also indicate an increase in catabolism of branched chain amino acids. The latter is supported by a simultaneous and moenomycin-specific up-regulation in the catabolic enzyme accD1 and in genes required for the biosynthesis of biotin, an important cofactor for catabolic processes. Branched chain amino acids are important precursors in secondary metabolite biosynthesis, and their availability can affect secondary metabolite production in S. coelicolor. Preferential use of these precursors for energy yielding catabolic processes in the moenomycin-treated cells could potentially explain the lack of activation by moenomycin of the secondary metabolic pathways for CDA and Red.
Induction of antibiotic biosynthesis following antibiotic treatment: chemical warfare in the soil?
S. coelicolor M600 produces three antibiotics, actinorhodin (Act), undecylprodiginine (Red) and the calcium-dependent antibiotic CDA. Strikingly, expression of the regulatory genes required for activating transcription of all three biosynthetic gene clusters was up-regulated within 30-60 min of exposure to both vancomycin and bacitracin, and was accompanied by a marked increase in transcription of the Red and CDA biosynthesis genes. Moenomycin also induced the CDA biosynthetic gene cluster. Production of the pigmented antibiotics Act and Red in response to vancomycin treatment was readily detectable in an agar disc assay (see Figure 4). S. coelicolor shares its native soil environment with a host of other microbial strains capable of secreting their own antibiotic metabolites into their immediate surroundings, and it appears to have evolved to respond by activating production of its own bioactive metabolites. Whether it does so as a form of retaliatory chemical warfare or as a means of inter-species communication is open to speculation.
Antibiotics as signalling molecules?
The production of metabolites with antibiotic activity is a widespread trait of soil microorganisms, and must confer some evolutionary advantage. While some of these compounds may indeed be used to inhibit competing species, they may also act as signalling molecules, particularly at sub-inhibitory concentrations [74–77], where they can modulate the transcriptional profiles of other bacteria [78–80]. The significant transcriptional response of S. coelicolor to treatment with moenomycin at levels corresponding to <5% of the MIC may represent another such example.
This work has used the novel approach of identifying common and unique transcriptional changes following the treatment of bacteria with three antibiotics that share the same key cellular target (cell wall biosynthesis) but differ at the molecular level in their mode of action. The characterised responses provide a wealth of information pertinent to generalised and antibiotic-specific cellular processes important for resistance and adaptation to the compounds used. Regulatory networks known to govern responses to environmental and nutritional stresses were shown to be at the heart of the common antibiotic response, and likely help cells survive until any specific resistance mechanisms are fully functional. Attenuation of these networks by mutation increased susceptibility to antibiotic treatment. Genes with putative roles in the extracellular synthesis and remodelling of peptidoglycan exhibited antibiotic-specific transcriptional changes, as did genes encoding two-component sensor systems and sigma factors. Antibiotic susceptibility studies using mutants constructed on the basis of the transcriptome profiling confirmed a role for several genes in antibiotic resistance, validating the usefulness of the approach. This included genes encoding two novel ABC transport systems required for bacitracin resistance, and a penicillin-binding protein and lysyl-tRNA synthetase that contribute to both bacitracin and vancomycin resistance. These represent new leads for further improving our understanding of how bacteria resist antimicrobial treatments at the molecular level, knowledge which may well be vital in enabling the future development of more effective antibiotics.
Materials and methods
Bacterial strains and culture conditions
S. coelicolor M600, a prototrophic plasmid-free derivative of S. coelicolor A3(2), was used for the microarray study. Spores of S. coelicolor M600 were germinated and grown to mid-log phase (OD600 nm approx. 0.5) in NMMP as described previously [19, 43, 44]. Biological triplicate cultures were treated at this point by adding a sub-lethal concentration (10 μg/ml) of vancomycin, bacitracin or moenomycin A and samples taken 0, 30, 60, 90 minutes after addition of the antibiotic. A negative control which received no antibiotic treatment was also performed. All of the antibiotics used for this study except moenomycin A were purchased from Sigma-Aldrich. Purified moenomycin A was a kind gift from Professor Peter Walzel in the University of Leipzig. All strains used are listed in Table 3.
Paper disc assays were used to study the induction of pigmented antibiotics by vancomycin. Spore lawns of M600 and M570 were spread on SMMS agar plates , and paper discs impregnated with vancomycin (or blank) were immediately applied to the surface. Cultures were incubated for 3 days at 30°C.
RNA isolation and DNA microarray analysis
RNA was isolated from liquid cultures according to Hesketh et al. (2007) . Purified total RNA (10 μg) was processed into labelled and fragmented cDNA for hybridisation to Streptomyces diS_div712a Affymetrix GeneChip arrays as previously described . Following scanning of the arrays, the data quality was verified using a variety of tools including the 'affyPLM', 'affy' and 'simpleaffy' packages for the statistical computing environment R , quality control methods available within GeneSpring 11 (Agilent), and data from report files generated in the Affymetrix Genechip Operating Software. One replicate from the bacitracin experiment failed quality control and was excluded from further analysis, leaving duplicate data for this antibiotic treatment. All microarray data are available from ArrayExpress (http://www.ebi.ac.uk/arrayexpress/) under accession number E-MEXP-3032.
To identify differentially expressed genes, array data were first imported into GeneSpring 11 with normalisation using the Robust Multichip Average (RMA) algorithm of Irizarry et al. . The data were then filtered to remove genes deemed to be expressed at a level below reliable quantification by determining those with a raw signal value below a defined background cut-off value of 100 in all samples. The results were further filtered to remove genes deemed to be unresponsive to the conditions under test in the experiment by identifying those with normalised expression values between 0.667-1.5 (1.5-fold change limit) in all conditions. The filtered data (2789 genes) were then subjected to two-way ANOVA to identify genes significantly altered under the experimental conditions, contrasting the expression values in each antibiotic-treated sample with the corresponding time point from the untreated control. This was performed using the parametric test option with a false discovery rate of P < 0.01. P-values were corrected using the Benjamini and Hochberg false discovery rate multiple testing correction procedure, and were computed asymptotically. Details of the statistical calculations used in the software can be accessed through the manufacturer's manual.
Gene ontology (GO) analysis was performed on lists of differentially expressed genes to identify over-represented biological functions or processes using the S. coelicolor GOA annotation from EMBL-EBI (ftp://ftp.ebi.ac.uk/pub/databases/GO/goa/proteomes/84.S_coelicolor.goa) as listed on 17/12/2009. GO analyses were realised using the Ontologizer tool (http://compbio.charite.de/index.php/ontologizer2.html), performing term-for-term hypergeometric testing with the Benjamini and Hochberg false discovery rate correction, or a Bayesian modelling approach using the model-based gene set analysis (MGSA) option . Statistical comparison of lists of differentially expressed genes to in-house curated lists of functionally related genes was performed using the 'find similar entity lists' tool of Genespring 11. The results (summarised in Tables 1 and 2) determined whether the number of genes from a particular functional class were over-represented in the data, and therefore biologically meaningful, or could have occurred by chance. Absence of a functional class in Tables 1 and 2 therefore does not mean that all genes in that class are absent in the data, but rather that the number of genes present could have arisen simply by chance at the 5% probability level used. Quality threshold (QT) clustering was performed in Genespring 7.3.
Construction of mutant strains for determining their antibiotic susceptibility
Targeted gene deletions were constructed using the PCR-directed mutagenesis procedure of Gust et al.  or by using transposon mutated cosmids from the library of Bishop et al. . Briefly, for PCR-directed replacement of genes with an antibiotic resistance cassette (e.g. aac(3)IV for apramycin resistance, hyg for hygromycin resistance), cosmid DNA in E. coli was targeted with a disruption cassette created by PCR using gene-specific primers and pIJ773 (for aac(3)IV) or pIJ10007 (for hyg) as template. Candidate antibiotic resistant clones were verified by PCR and restriction digestion, and conjugated into S. coelicolor M600. Double crossover integrants were selected as apramycin/hygromycin resistant, kanamycin sensitive clones, and verified by PCR. For transposon mutants, mutated cosmids from the library were conjugated into S. coelicolor M600 and double crossover integrants were selected as apramycin resistant, kanamycin sensitive clones.
Assay for determining the minimum inhibitory concentration (MIC) of the antibiotics
Minimum inhibitory concentration (MIC) values were evaluated using a method similar to that of Andrews (2001) . McFarland turbidity standards of 0 (1.0 ml distilled water), 0.5 (5.0 μl 1.175% BaCl2, 995 μl 1% H2SO4), 1 (10.0 μl 1.175% BaCl2, 990 μl 1% H2SO4), and 2 (20 μl 1.175% BaCl2, 980 μl 1% H2SO4) were prepared. S. coelicolor spores from 20% glycerol stocks were added to 1.0 ml distilled water until the turbidity of the spore suspension was between McFarland standards 1 and 2. Spore suspensions were stored at 4°C and used within 1 week, during which time they showed no loss of viability. 5.0 μl aliquots of the spore suspensions were spotted onto a series of MMCGT agar plates  containing the antibiotic under test at a range of concentrations (vancomycin 30-150 μg/ml in 10 μg/ml increments; bacitraicin 10-45 μg/ml in 5 μg/ml increments; moenomycin 0-300 μg/ml in 20 and 50 μg/ml increments). The spots were allowed to dry, and the plates were incubated at 30°C. MIC values were determined by visual inspection of growth over the range of antibiotic concentrations after 48 hrs incubation.
Analysis of intracellular nucleotide pools
Quantification of intracellular nucleotides was performed by HPLC as described in Hesketh et al. . Duplicate cultures of S. coelicolor grown to mid-exponential phase (OD600 nm approx. 0.5) and treated with 10 μg/ml vancomycin were extracted and analysed at 0, 15 30 and 60 min following addition of the antibiotic.
analysis of variance
model-based gene set analysis
minimum inhibitory concentration
minimal medium supplemented with casamino acids, glucose, tiger milk
methicillin resistant Staphylococcus aureus
Bentley SD, Chater KF, Cerdeno-Tarraga AM, Challis GL, Thomson NR, James KD, Harris DE, Quail MA, Kieser H, Harper D, Bateman A, Brown S, Chandra G, Chen CW, Collins M, Cronin A, Fraser A, Goble A, Hidalgo J, Hornsby T, Howarth S, Huang CH, Kieser T, Larke L, Murphy L, Oliver K, O'Neil S, Rabbinowitsch E, Rajandream MA, Rutherford K, et al: Complete genome sequence of the model actinomycete Streptomyces coelicolor A3(2). Nature. 2002, 417: 141-147. 10.1038/417141a.
Ikeda H, Ishikawa J, Hanamoto A, Shinose M, Kikuchi H, Shiba T, Sakaki Y, Hattori M, Omura S: Complete genome sequence and comparative analysis of the industrial microoganism Streptomyces avermitilis. Nat Biotechnol. 2003, 21: 526-531. 10.1038/nbt820.
Chang S, Sievert DM, Jefferey C, Hageman JC, Boulton ML, Tenover FC, Downes FP, Shah S, Rudrick JT, Pupp GR, Brown WJ, Cardo D, Fridkin SK: Infection with vancomycin-resistant Staphylococcus aureus containing the vanA resistance gene. N Engl J Med. 2003, 348: 1342-1347. 10.1056/NEJMoa025025.
Pearson H: 'Superbug' hurdles key drug barrier. Nature. 2002, 418: 469-470.
Tenover FC, Weigel LM, Appelbaum PC, McDougal LK, Chaitram J, McAllister S, Clark N, Killgore G, O'Hara CM, Jevitt L, Patel JB, Bozdogan B: Vancomycin-resistant Staphylococcus aureus isolate from a patient in Pennsylvania. Antimicrob Agents Chemother. 2004, 48: 275-280. 10.1128/AAC.48.1.275-280.2004.
Weigel LM, Clewell DB, Gill SR, Clark NC, McDougal LK, Flannagan SE, Kolonay JF, Shetty J, Killogore GE, Tenover FC: Genetic analysis of a high-level vancomycin-resistant isolate of Staphylococcus aureus. Science. 2003, 28: 1567-1571.
Koteva K, Hong H-J, Wang XD, Nazi I, Hughes D, Naldrett MJ, Buttner MJ, Wright GD: A vancomycin photoprobe identifies the histidine kinase VanSsc as a vancomycin receptor. Nature Chem Biol. 2010, 6: 327-329.
Kohanski MA, Dwyer DJ, Hayete B, Lawrence CA, Collins JJ: A common mechanism of cellular death induced by bactericidal antibiotics. Cell. 2007, 130: 797-810. 10.1016/j.cell.2007.06.049.
McCallum N, Spehar G, Bischoff M, Berger-Bächi B: Strain dependence of the cell wall-damage induced stimulon in Staphylococcus aureus. Biochim Biophys Acta. 2006, 1760: 1475-1481.
O'Neill AJ, Lindsay JA, Gould K, Hinds J, Chopra I: Transcriptional signature following inhibition of early-stage cell wall biosynthesis in Staphylococcus aureus. Antimicrob Agents Chemother. 2009, 53: 1701-1704. 10.1128/AAC.01309-08.
Pietiäinen M, François P, Hayyryläinen H-L, Tangomo M, Sass V, Sahl HG, Schrenzel J, Kontinen VP: Transcriptome analysis of the responses of Staphylococcus aureus to antimicrobial peptides and characterization of the roles of vraDE and vraSR in antimicrobial resistance. BMC Genomics. 2009, 10: 429-443. 10.1186/1471-2164-10-429.
Provvedi R, Boldrin F, Falciani F, Palù G, Manganelli R: Global transcriptional response to vancomycin in Mycobacterium tuberculosis. Microbiology. 2009, 155: 1093-1102. 10.1099/mic.0.024802-0.
Utaida S, Dunman PM, Murphy ME, Projan SJ, Singh VK, Jayaswal RK, Wilkinson BJ: Genome-wide transcriptional profiling of the response of Staphylococcus aureus to cell wall active antibiotics reveals a cell wall-stress stimulon. Microbiology. 2003, 149: 2719-2732. 10.1099/mic.0.26426-0.
Hutter B, Schaab C, Alktrecht S, Borgmann M, Brunner NA, Freiberg C, Ziegelbauer K, Rock CO, Ivanov I, Loferer H: Prediction of mechanisms of action of antibacterial compounds by gene expression profiling. Antimicrob Agents Chemother. 2004, 48: 2838-2844. 10.1128/AAC.48.8.2838-2844.2004.
Bugg TDH, Wright GD, Dutka-Malen S, Arthur M, Courvalin P, Walsh CT: Molecular basis for vancomycin resistance in Enterococcus faecium BM4147: biosynthesis of a depsipeptide peptidoglycan precursor by vancomycin resistance proteins VanH and VanA. Biochemistry. 1991, 30: 10408-10415. 10.1021/bi00107a007.
van Heijenoort J: Formation of the glycan chains in the synthesis of bacterial peptidoglycan. Glycobiology. 2001, 11: 25R-36R. 10.1093/glycob/11.3.25R.
Azevedo EC, Rios EM, Fukushima K, Campos-Takaki GM: Bacitracin production by a new strain of Bacillus subtilis. Extraction, purification and characterization. Appl Biochem Biotechnol. 1993, 42: 1-7. 10.1007/BF02788897.
Ishihara H, Takoh M, Nishibayashi R, Sato A: Distribution and variation of bacitracin synthetase gene sequences in laboratory stock strains of Bacillus licheniformis. Curr Microbiol. 2002, 45: 18-23. 10.1007/s00284-001-0041-5.
Hong H-J, Paget MSB, Buttner MJ: A signal transduction system in Streptomyces coelicolor that activated expression of a putative cell wall glycan operon in response to vancomycin and other cell wall-specific antibiotics. Mol Microbiol. 2002, 44: 1199-1211. 10.1046/j.1365-2958.2002.02960.x.
Paget MSB, Chamberlin L, Atrih A, Foster SJ, Buttner MJ: Evidence that the extracytoplasmic function sigma factor σE is required for normal cell wall structure in Streptomyces coelicolor A3(2). J Bacteriol. 1999, 181: 204-211.
Hesketh A, Chen WJ, Ryding J, Chang S, Bibb M: The global role of ppGpp synthesis in morphological differentiation and antibiotic production in Streptomyces coelicolor A3(2). Genome Biol. 2007, 8: R161-10.1186/gb-2007-8-8-r161.
Komatsu M, Takano H, Hiratsuka T, Ishigaki Y, Shimada K, Beppu T, Ueda K: Proteins encoded by the conservon of Streptomyces coelicolor A3(2) comprise a membrane-associated heterocomplex that resembles eukaryotic G protein-coupled regulatory system. Mol Microbiol. 2006, 62: 1534-1546. 10.1111/j.1365-2958.2006.05461.x.
Chakraburtty R, Bibb J: The ppGpp synthetase gene (relA) of Streptomyces coelicolor A3(2) plays a conditional role in antibiotic production and morphological differentiation. J Bacteriol. 1997, 179: 5854-5861.
Bucca G, Laing E, Mersinias V, Allenby N, Hurd D, Holdstock J, Brenner V, Harrison M, Smith CP: Development and application of versatile high density microarrays for genome-wide analysis of Streptomyces coelicolor: characterization of the HspR regulon. Genome Biol. 2009, 10: R5-10.1186/gb-2009-10-1-r5.
Lee EJ, Cho YH, Kim HS, Ahn BE, Roe J-H: Regulation of σB by an anti- and an anti-anti-sigma factor in Streptomyces coelicolor in response to osmotic stress. J Bacteriol. 2004, 186: 8490-8498. 10.1128/JB.186.24.8490-8498.2004.
Paget MS, Molle V, Cohen G, Aharonowitz Y, Buttner MJ: Defining the disulphide stress response in Streptomyces coelicolor A3(2): identification of the σR regulon. Mol Microbiol. 2001, 42: 1007-1020. 10.1046/j.1365-2958.2001.02675.x.
Bursy J, Pierik AJ, Pica N, Bremer E: Osmotically induced synthesis of the compatible solute hycroxyectoine is mediated by an evolutionarily conserved ectoine hydroxylase. J Biol Chem. 2007, 282: 31147-31155. 10.1074/jbc.M704023200.
da Costa MS, Santos H, Galinski EA: An overview of the role and diversity of compatible solutes in bacteria and archaea. Adv Biochem Eng Biotechnol. 1998, 61: 117-153.
Galinski EA, Trüper HG: Microbial behaviour in salt-stressed ecosystems. FEMS Microbiol Rev. 1994, 15: 95-108. 10.1111/j.1574-6976.1994.tb00128.x.
García-Estepa R, Argandoña M, Reina-Bueno M, Capote N, Iglesias-Guerra F, Nieto JJ, Vargas C: The ectD gene, which is involved in the synthesis of the compatible solute hydroxyectoine, is essential for thermoprotection of the halophilic bacterium Chromohalobacter salexigens. J Bacteriol. 2006, 188: 3774-3784. 10.1128/JB.00136-06.
Kuhlmann AU, Bremer E: Osmotically regulated synthesis of the compatible solute ectoine in Bacillus pasteurii and related Bacillus spp. Appl Environ Microbiol. 2002, 68: 772-783. 10.1128/AEM.68.2.772-783.2002.
Kuhlmann AU, Bursy J, Gimpel S, Hoffmann T, Bremer E: Synthesis of the compatible solute ectoine in Virgibacillus pantothenticus is triggered by high salinity and low growth temperature. Appl Environ Microbiol. 2008, 74: 4560-4563. 10.1128/AEM.00492-08.
Bursy J, Kuhlmann AU, Pittelkow M, Hartmann H, Jebbar M, Pierik AJ, Bremer E: Synthesis and uptake of the compatible solutes ectoine and 5-hydroxyectoine by Streptomyces coelicolor A3(2) in response to salt and heat stresses. Appl Environ Microbiol. 2008, 74: 7286-7296. 10.1128/AEM.00768-08.
Lee E-J, Karoonuthaisiri N, Kim HS, Park JH, Cha CJ, Kao CM, Roe JH: A master regulator σB governs osmotic and oxidative response as well as differentiation via a network of sigma factors in Streptomyces coelicolor. Mol Microbiol. 2005, 57: 1252-1264. 10.1111/j.1365-2958.2005.04761.x.
Bugg TDH: Bacterial peptidoglycan biosynthesis and its inhibition. Comprehensive Natural Products Chemistry. Edited by: Meth-Cohn O, Barton D, Nakanishi K. 1999, Elsevier Science
Foster SJ, Popham DL: Structure and synthesis of cell wall, spore cortex, teichoic acid, S-layer, and capsules. Bacillus subtilis and its relatives: from Genes to Cells. Edited by: Sonenshein L, Losick R, Hoch JA. 2001, American Soceity for Microbiology Press
Goffin C, Ghuysen JM: Multimodular penicillin-binding proteins: An enigmatic family of orthologs and paralogs. Microbiol Mol Biol Review. 1998, 62: 1079-1093.
Ishino F, Matsuhashi M: Peptidoglycan synthetic enzyme activities of highly purified penicillin-binding protein 3 in Escherichia coli: a septum-forming reaction sequence. Biochem Biophys Res Commun. 1981, 101: 905-911. 10.1016/0006-291X(81)91835-0.
McPherson DC, Driks A, Popham DL: Two class A high-molecular weight penicillin-binding proteins of Bacitllus subtilis play redundant roles in sporulation. J Bacteriol. 2001, 183: 6046-6053. 10.1128/JB.183.20.6046-6053.2001.
Murray T, Popham DL, Pearson CB, Hand AR, Setlow P: Analysis of outgrowth of Bacillus subtilis spores lacking penicillin-binding protein 2a. J Bacteriol. 1998, 180: 6493-6502.
Murray T, Popham DL, Setlow P: Bacillus subtilis cells lacking penicillin-binding protein 1 required increased levels of divalent cations for growth. J Bacteriol. 1998, 180: 4555-4563.
Ruiz N: Bioinformatics identification of MurJ (MviN) as the peptidoglycan lipid II flippase in Escherichia coli. Proc Natl Acad Sci. 2008, 105: 15553-15557. 10.1073/pnas.0808352105.
Hong H-J, Hutchings MI, Neu JM, Wright GD, Buttner MJ: Characterization of an inducible vancomycin resistance system in Streptomyces coelicolor reveals a novel gene (vanK) required for drug resistance. Mol Microbiol. 2004, 52: 1107-1121. 10.1111/j.1365-2958.2004.04032.x.
Hong H-J, Hutchings MI, Hill L, Buttner MJ: The role of the novel Fem protein VanK in vancomycin resistance in Streptomyces coelicolor. J Biol Chem. 2005, 280: 13055-13061.
Haiser HJ, Yousef MR, Elliot MA: Cell wall hydrolases affect germination, vegetative growth, and sporulation in Streptomyces coelicolor. J Bacteriol. 2009, 191: 6501-6512. 10.1128/JB.00767-09.
Noens EEE, Mersinias V, Traag BA, Smith CP, Koerten HK, van Wezel GP: SsgA-like proteins determine the fate of peptidoglycan during sporulation of Streptomyces coelicolor. Mol Microbiol. 2005, 58: 929-944. 10.1111/j.1365-2958.2005.04883.x.
Traag BA, van Wezel GP: The SsgA-like proteins in actinomycetes: small proteins up to a big task. Antonie van Leeuwenh.
Cossart P, Jonquières R: Sortase, a universal target for therapeutic agents against Gram-positive bacteria?. Proc Natl Acad Sci. 2000, 97: 5013-5015. 10.1073/pnas.97.10.5013.
Boekhorst J, de Been MWHJ, Ckeerebezem M, Siezen RJ: Genome-wide detection and analysis of cell wall-bound proteins with LPxTG-like sorting motifs. J Bacteriol. 2005, 184: 4928-4934.
Elliot MA, Karoonuthaisiri N, Huang J, Bibb MJ, Cohen SN, Kao CM, Buttner MJ: The chaplins: a family of hycrophobic cell-surpface proteins involved in aerial mycelium formation in Streptomyces coelicolor. Genes Dev. 2003, 17: 1727-1740. 10.1101/gad.264403.
Gramajo HC, Takano E, Bibb MJ: Stationary-phase production of the antibiotic actinorhodin in Streptomyces coelicolor A3(2) is transcriptionally regulated. Mol Microbiol. 1993, 7: 837-845. 10.1111/j.1365-2958.1993.tb01174.x.
Ryding NJ, Anderson TB, Champness WC: Regulation of the Streptomyces coelicolor calcium-dependent antibiotic by absA, encoding a cluster-linked two-component system. J Bacteriol. 2002, 184: 794-805. 10.1128/JB.184.3.794-805.2002.
Takano E, Gramajo HC, Strauch E, Andres N, White J, Bibb MM: Transcriptional regulation of the redD transcriptional activator gene accounts for growth-phase-dependent production of the antibiotic undecylprodigiosin in Streptomyces coelicolor A3(2). Mol Microbiol. 1992, 6: 2797-2804. 10.1111/j.1365-2958.1992.tb01459.x.
Kang JG, Hahn MY, Ishihama A, Roe JH: Identification of sigma factors for growth phase-related promoter selectivity of RNA polymerases from Streptomyces coelicolor A3(2). Nucleic Acids Res. 1997, 25: 2566-2573. 10.1093/nar/25.13.2566.
Bibb MJ, Molle V, Buttner MJ: σBldN, an extracytoplasmic function RNA polymerase sigma factor required for aerial mycelium formation in Streptomyces coelicolor A3(2). J Bacteriol. 2002, 182: 4606-4616.
Gaballa A, Wang T, Ye RW, Helmann JD: Functional analysis of the Bacillus subtilis Zur regulon. J Bacteriol. 2002, 184: 6508-6514. 10.1128/JB.184.23.6508-6514.2002.
Patzer SI, Hantke K: The ZnuABC high-affinity zinc uptake system and its regulator Zur in Escherichia coli. Mol Microbiol. 1998, 28: 1199-1210. 10.1046/j.1365-2958.1998.00883.x.
Kallifidas D, Pascoe B, Owen GA, Strain-Damerell CM, Hong H-J, Paget MSB: The zinc-responsive regulator Zur controls expression of the Coelibactin gene cluster in Streptomyces coelicolor. J Bacteriol. 2009, 192: 608-611.
Owen GA, Pascoe B, Kallifidas D, Paget MS: Zinc-responsive regulation of alternative ribosomal protein genes in Streptomyces coelicolor involves zur and σR. J Bacteriol. 2007, 189: 4078-4086. 10.1128/JB.01901-06.
Shin JH, Oh SY, Kim S-J, Roe J-H: The zinc-responsive regulator Zur controls a zinc uptake system and some ribosomal proteins in Streptomcyes coelicolor A3(2). J Bacteriol. 2007, 189: 4070-4077. 10.1128/JB.01851-06.
Navratna V, Nadig S, Sood V, Prasad K, Arakere G, Gopal B: Molecular basis for the role of Staphylococcus aureus penicillin binding protein 4 in antimicrobial resistance. J Bacteriol. 2010, 192: 134-144. 10.1128/JB.00822-09.
Arbeloa A, Segal H, Hugonnet JE, Josseaume N, Dubost L, Brouard JP, Gutmann L, Mengin-Lecreulx D, Arthur M: Role of class A penicillin-binding proteins in PBP5-mediated β-lactam resistance in Enterococcus faecalis. J Bacteriol. 2004, 186: 1221-1228. 10.1128/JB.186.5.1221-1228.2004.
RajBhandary UL, Söll D: Aminoacy-tRNAs, the bacterial cell envelope, and antibiotics. Proc Natl Acad Sci. 2008, 105: 5285-5286. 10.1073/pnas.0801193105.
Nishi H, Komatsuzawa H, Fujiwara T, McCallum N, Sugai M: Reduced content of lysyl-phosphatidylglycerol in the cytoplasmic membrane affects susceptibility to moenomycin, as well as vancomycin, gentamicin, and antimicrobial peptides in Staphylococcus aureus. Antimicrob Agents Chemother. 2004, 48: 4800-4807. 10.1128/AAC.48.12.4800-4807.2004.
Gottilieb CT, Thomsen LE, Ingmer H, Mygind PH, Kristensen H-H, Gram L: Antimicrobial peptides effectively kill a broad spectrum of Listeria monocytogenes and Staphylococcus aureus strains independently of origin, sub-type, or virulence factor expression. BMC Microbiol. 2006, 8: 205-
Kucharczyk M, Brzezowska M, Maciag A, Lis T, Jezowska-Bojczuk M: Structural features of the Cu(2+)-vancomycin complex. J Inorg Biochem. 2008, 102: 936-942. 10.1016/j.jinorgbio.2007.12.014.
Swiatek M, Valensin D, Migliorini C, Gaggelli E, Valensin G, Jezowska-Bojczuk M: Unusual binding ability of vancomycin towards Cu2+ ions. Dalton Trans. 2005, 7: 3808-3813.
Bernard R, Ghachi ME, Mengin-Lecreulx D, Chippaux M, Denizot F: BcrC from Bacillus subtilis acts as an undecaprenyl pyrophosphate phosphatase in bacitracin resistance. J BIol Chem. 2005, 280: 28852-28857. 10.1074/jbc.M413750200.
Cain BD, Norton PJ, Eubanks W, Nick HS, Allen CM: Amplication of the bacA gene confers bacitracin resistance to Escherichia coli. J Bacteriol. 1993, 175: 3784-3789.
Cao M, Helmann JD: Regulation of the Bacillus subtilis bcrC bacitracin resistance gene by two extracytoplasmic function σ factors. J Bcteriol. 2002, 184: 6123-6129. 10.1128/JB.184.22.6123-6129.2002.
Chalker AF, Ingraham KA, Swayne Lunsford R, Bryant AP, Bryant J, Wallis NG, Broskey JP, Pearson SC, Homes DJ: The bacA gene, which determines bacitracin susceptibility in Streptococcus pneumonia and Staphylococcus aureus, is also required for virulence. Microbiology. 2000, 146: 1547-1553.
Ohki R, Tateno K, Okada Y, Okajima H, Asai K, Sadaie Y, Murata M, Aiso T: A bacitracin-resistant Bacillus subtilis gene encodes a homologue of the membrane-spanning subunit of the Bacillus licheniformis ABC transporter. J Bacteriol. 2003, 185: 51-59. 10.1128/JB.185.1.51-59.2003.
Potter CA, Baumberg S: End-product control of enzymes of branched-chain amino acid biosynthesis in Streptomyces coelicolor. Microbiology. 1996, 142: 1945-1952. 10.1099/13500872-142-8-1945.
Davies J: Are antibiotics naturally antibiotics?. J Ind Microbiol Biotechnol. 2006, 33: 496-499. 10.1007/s10295-006-0112-5.
Davies J: Small molecules: The lexicon of biodiversity. J Biotechnol. 2007, 129: 3-5. 10.1016/j.jbiotec.2006.11.023.
Yim G, Wang HH, Davies J: The truth about antibiotics. Int J Med Microbiol. 2006, 296: 163-170. 10.1016/j.ijmm.2006.01.039.
Yim G, Wang HH, Davies J: Antibiotics as signalling molecules. Phil Trans R Soc. 2007, 362: 1195-1200. 10.1098/rstb.2007.2044.
Goh E-B, Yim G, Tsui W, McClure J, Surette MG, Davies J: Transcriptional modulation of bacterial gene expression by subinhibitory concentrations of antibiotics. Proc Natl Acad Sci. 2002, 99: 17025-17030. 10.1073/pnas.252607699.
Mesak LR, Miao V, Davies J: Effects of subinhibitory concentrations of antibiotics on SOS and DNA repair gene expression in Staphylococcus aureus. Antimicrob Agents Chemother. 2008, 52: 3394-3397. 10.1128/AAC.01599-07.
Yim G, de la Cruz F, Spiegelman GB, Davies J: Transcription modulation of Salmonella enteric serovar typhimurium promoters by sub-MIC levels of rifampin. J Bacteriol. 2006, 188: 7988-7991. 10.1128/JB.00791-06.
Kieser T, Bibb MJ, Buttner MJ, Chater KF, Hopwood DA: Practical Streptomyces Genetics. 2000, John Innes Foundation
The R project for statistical computing. [http://www.r-project.org/]
Irizarry RA, Bolstad BM, Collin F, Cope LM, Hobbs B, Speed TP: Summaries of Affymetrix GeneChip probe level data. Nucleic Acids Res. 2003, 31: e15-10.1093/nar/gng015.
Bauer S, Gagneur J, Robinson PN: Going Bayesian: model-based gene set analysis of genome-scale data. Nucleic Acids Res. 2010, 38: 3523-3532. 10.1093/nar/gkq045.
Gust B, Challis GL, Fowler K, Kieser T, Chater KF: PCR-targeted Streptomyces gene replacement identifies a protein domain needed for biosynthesis of the sesquiterpene soil odor geosmin. Proc Natl Acad Sci. 2003, 100: 1541-1546. 10.1073/pnas.0337542100.
Bishop A, Fielding S, Dyson P, Herron P: Systematic Insertional Mutagenesis of a Streptomycete Genome: A Link Between Osmoadaptation and Antibiotic Production. Genome Res. 2004, 14: 893-900. 10.1101/gr.1710304.
Andrews J: Determination of minimum inhibitory concentrations. J Antimicrob Chemother. 2001, 48: 5-16.
Rajaram S, Oono Y: NeatMap: non-clustering heat map alternatives in R. BMC Bioinformatics. 2010, 11: 45-10.1186/1471-2105-11-45.
HJH was supported by the Royal Society and the Medical Research Council. AH and MB were supported by the BBSRC. We thank Peter Welzel in the University of Leipzig for his kind gift of Moenomycin A for the microarray study. We also thank Mark Buttner and Gabriela Kelemen for kindly supplying strains to test for antibiotic resistance, and are also grateful to both for helpful discussions.
The authors declare that they have no competing interests.
HJH conceived and designed the study, performed the microarray experiments, assisted with data analysis and interpretation, and lead the writing of the manuscript. AH assisted with the microarray experiments, performed the microarray data analysis and interpretation, drafted the manuscript, and contributed to the construction and phenotypic characterisation of mutant strains. MB participated in the experimental design, and contributed to the writing of the manuscript. CH, JM, GN and NT all contributed to the construction and phenotypic characterisation of mutant strains. All authors read and approved the final manuscript.
Electronic supplementary material
Additional file 1:Summary of transcriptome analysis of genes significantly differently expressed in response to vancomycin treatment. (XLSX 2 MB)
Additional file 2:Summary of transcriptome analysis of genes significantly differently expressed in response to bacitracin treatment. (XLSX 1 MB)
Additional file 3:Summary of transcriptome analysis of genes significantly differently expressed in response to moenomycin treatment. (XLSX 1 MB)
Additional file 4:Normalised expression data for the differentially expressed gene groups illustrated in the venn diagram of Figure1B. (XLSX 788 KB)
Additional file 5:Term for term GO analysis of the genes significantly differently expressed in response to vancomycin. (XLSX 589 KB)
Additional file 6:Term for term GO analysis of the genes significantly differently expressed in response to bacitracin. (XLSX 431 KB)
Additional file 7:Term for term GO analysis of the genes significantly differently expressed in response to moenomycin. (XLSX 355 KB)
Additional file 8:MGSA GO analysis of the genes significantly differently expressed in response to vancomycin. (XLSX 465 KB)
Additional file 9:MGSA GO analysis of the genes significantly differently expressed in response to bacitracin. (XLSX 256 KB)
Additional file 10:MGSA GO analysis of the genes significantly differently expressed in response to moenomycin. (XLSX 277 KB)
Additional file 11:Term for term GO analysis of the differentially expressed gene groups illustrated in the venn diagram of Figure1B. (XLSX 71 KB)
Additional file 12:Analysis of the differentially expressed gene groups illustrated in the venn diagram of Figure1B, looking for significant similarity with in-house curated lists of functionally related genes. (XLSX 26 KB)
Additional file 13:Genes whose expression following treatment with all three antibiotics is closely correlated (Pearson correlation > 0. 9) with transcription of the cell wall stress sigma factor σE (SCO3356). (XLSX 64 KB)
Additional file 14:Heatmaps summarising the expression of genes significantly differently expressed in response to drug treatment, grouped according to related function. (PDF 1 MB)
About this article
Cite this article
Hesketh, A., Hill, C., Mokhtar, J. et al. Genome-wide dynamics of a bacterial response to antibiotics that target the cell envelope. BMC Genomics 12, 226 (2011). https://doi.org/10.1186/1471-2164-12-226
- Minimum Inhibitory Concentration
- Branch Chain Amino Acid
- Sigma Factor