|
|
|
Vol. 11, Issue 1, 28-42, January 2001
LETTER
| ||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
| |
ABSTRACT |
|---|
|
|
|---|
The responses of Haemophilus influenzae to DNA gyrase inhibitors were analyzed at the transcriptional and the translational level. High-density microarrays based on the genomic sequence were used to monitor the expression levels of >80% of the genes in this bacterium. In parallel the proteins were analyzed by two-dimensional electrophoresis. DNA gyrase inhibitors of two different functional classes were used. Novobiocin, as a representative of one class, inhibits the ATPase activity of the enzyme, thereby indirectly changing the degree of DNA supercoiling. Ciprofloxacin, a representative of the second class, obstructs supercoiling by inhibiting the DNA cleavage-resealing reaction. Our results clearly show that different responses can be observed. Treatment with the ATPase inhibitor Novobiocin changed the expression rates of many genes, reflecting the fact that the initiation of transcription for many genes is sensitive to DNA supercoiling. Ciprofloxacin mainly stimulated the expression of DNA repair systems as a response to the DNA damage caused by the stable ternary complexes. In addition, changed expression levels were also observed for some genes coding for proteins either annotated as "unknown function" or "hypothetical" or for proteins not directly involved in DNA topology or repair.
[The sequence data described in this paper have been submitted to the EMBL data library under accession nos. AJ297131 and AL135960.]
| |
INTRODUCTION |
|---|
|
|
|---|
In the last few years, the total genomic sequences of many
prokaryotes were determined. In parallel, techniques were developed that allow the monitoring of the expression levels of thousands of
genes simultaneously. One of these techniques,
developed by Affymetrix, is based on the principle of photolithography
and conventional oligonucleotide synthesis, which allows the synthesis of short oligonucleotides in high-density arrays directly on a solid
surface (Fodor et al. 1991
, 1993
; Chee et al. 1996
; Lockhart et al.
1996
; Wodicka et al. 1997
). Total genomic DNA sequences are used to
select sets of unique oligonucleotides to represent each open reading
frame (ORF). To further increase the sensitivity and specificity of
detection, a mismatch partner, which is identical except for a single
base difference at the central position, is synthesized for each
perfect-match oligonucleotide. These mismatch oligonucleotides serve as
internal controls for the specificity of the probes.
For the gene expression analysis described here, a high-density
microarray containing selected oligonucleotides for ~2000 genes from
the bacterium Streptococcus pneumoniae and for ~1800 genes
from Haemophilus influenzae was used (Fleischmann et al. 1995
). In addition, the microarray contains many control genes, sequence information from intergenic regions, and genes coding for
ribosomal and transfer RNA. A set of 25-mer oligonucleotides for a
specific gene usually includes 25 probe pairs (a probe pair consists of
the perfect match and the corresponding mismatch oligonucleotide) and
at least 20 probe pairs for very short genes. This microarray was used
to simultaneously determine the changes in RNA levels for all the genes
transcribed by H. influenzae following the addition of DNA
gyrase inhibitors.
Regulation at the transcriptional level is only one possibility for a cell to respond to changing growth conditions. Other regulatory mechanisms act at the level of mRNA translation. Therefore, two-dimensional polyacrylamide gel electrophoresis (2D-PAGE), another tool for expression analysis that has been facilitated by the availability of whole-genome sequences and new developments in mass spectrometry, was performed in parallel. This allowed the comparison of the changes at the protein level with those of the transcriptional pattern.
DNA gyrase (E.C. 5.99.1.3.), a prokaryotic topoisomerase II enzyme
essential for viability, consists of two subunits, A and B, the active
enzyme being an A2B2 tetrameric complex (for
reviews, see Reece and Maxwell 1991
; Menzel and Gellert 1994
; Luttinger 1995
; Roca 1995
; Sharma and Mondragon 1995
). The enzyme has no direct
mammalian counterpart and is the only enzyme known to be able to
introduce negative supercoils into DNA by using the energy derived from
ATP hydrolysis. A key step in this supercoiling reaction is the DNA
gyrase mediated cleavage of DNA. It has been shown that the class of
subunit A inhibitors, the quinolones and the pyrimido[1,6-a]benzimidazoles, interrupt the cleavage
and resealing cycle at the cleavage step (for reviews, see Hooper and
Wolfson 1991
; Hubschwerlen et al. 1992
; Hooper 1993
; Gmuender et al.
1995
, 1997
). However, there is evidence that quinolones acting in vivo have effects beyond the inhibition of DNA gyrase. They induce the
formation of a stable ternary complex consisting of the enzyme, DNA,
and the inhibitor, resulting in DNA damage, which in turn blocks
replication and transcription (for reviews, see Drlica and Zhao 1997
;
Maxwell 1999
). As a consequence, the expression of DNA repair systems,
mainly the SOS system, is induced (Piddock and Wise 1987
; Walters et
al. 1989
; Piddock et al. 1990
). Another class of DNA gyrase inhibitors,
the cyclothialidines and the coumarins, bind to the ATP binding site
located in the subunit B, thereby inhibiting the supercoiling activity
of the enzyme but leaving the DNA otherwise intact (Contreras and
Maxwell 1992
; Ali et al. 1993
; Goetschi et al. 1993
; Maxwell 1993
;
Nakada et al. 1994
; Ali et al. 1995
; Nakada et al. 1995
; Gormley et al.
1996
; Lewis et al. 1996
; Oram et al. 1996
; Tingey and Maxwell 1996
;
Kampranis et al. 1999
). The initiation of transcription of many genes
is sensitive to DNA supercoiling, often exhibiting an optimum with respect to the degree of supercoiling (Jovanovich and Lebowitz 1987
;
Steck et al. 1993
; Wang and Lynch 1993
).
The goals of this study were (1) to investigate if both classes of antibiotics, although inhibiting the same enzyme but through different mechanisms, induce different mechanism-related expression and translational patterns; (2) to cross-validate the two technologies; and (3) to test to what extent the use of both technologies in conjunction enhances the power of expression analysis.
| |
RESULTS AND DISCUSSION |
|---|
|
|
|---|
Methodological Studies
Sensitivity
The Affymetrix chip analysis detected transcripts for typically 70%-85% (in some experiments, even 90%) of the genes represented on the chips, indicating that most of the genes are transcribed independently of the growth conditions and that the sensitivity of the system is high enough to detect even low abundance transcripts. Cross-hybridization with S. pneumoniae oligonucleotides present on the same chip was negligible in all experiments. Two-dimensional gel electrophoresis followed by the detection of radioactive spots by using a Phosphorimager and the PDQuest program reproducibly detected 560 spots, which corresponds to ~30% of all theoretical gene products. Of these spots, 274 can be assigned to the corresponding gene by comparison with a 2D protein map (Langen et al. 2000Reproducibility: Transcriptional Imaging
To estimate the variation between expression patterns derived from RNA isolated from the same culture or from different cultures but grown under the same conditions, we hybridized reverse transcribed RNA from control cultures to the microarrays and normalized and analyzed the results by using the Affymetrix GeneChip software. RNA isolated from the same culture but reverse transcribed and hybridized independently shows a highly reproducible gene expression pattern and almost identical results after hybridization (Table 1; Fig. 1). About 90% of the ratios between the average differences of two normalized experiments lay between 0.66 and 1.5 with only a few transcripts being identified as increased or decreased. No significant changes, that is fold changes of the corresponding gene transcripts >2 or <-2, could be observed.
|
|
Reproducibility: Protein Quantification
To estimate the gel-to-gel reproducibility, we resolved a protein extract from a control culture on two 2D gels in parallel and matched the resulting gel images. A sample from a second culture obtained under the same conditions was also analyzed and the gel image was compared to the parallel gels. Of the 560 spots detected, the calculated change factor was >2 or <-2 for 185 spots (32%) when gels from the same sample were compared and 220 (39%) when comparing different samples. A comparison for the 274 spots for which an identity has been assigned resulted in values of 68 (25%) for the comparison of two gels from the same sample and 140 (50%) for the comparison of two different samples. These figures indicate that the reproducibility for the annotated spots is higher, most likely reflecting the stronger intensity of these spots, which makes their quantification more accurate. The reproducibility of the 2D-gel-based protein quantification is clearly inferior to that observed for transcriptional imaging by using Affymetrix chips. To take account of these findings, we chose to carry out five parallel experiments and to analyze two on Affymetrix chips and all five on 2D gels. Only parallel samples were compared to avoid misinterpretations resulting from culture-to-culture variability.Quantitative Comparison between Total mRNA Levels and Protein Synthesis Levels
In previous studies comparing mRNA and protein levels, samples were not prepared in parallel so that only an approximation of the reproducibility of mRNA quantification and quantification of proteins could be made (Anderson and Seilhamer 1997Comparison between Differential RNA and Protein Synthesis
To estimate the degree of consistency between the results obtained by RNA and protein quantification, we calculated the average fold changes and compared for the Novobiocin treatment at the higher concentration. Genes were selected that showed significant changes by using either one of both detection methods (143 genes). The values were grouped into induced (fold change >2), repressed (fold change <-2), and unchanged genes/proteins. In 55% of all cases, the results fell into the same category. About 40% of all changes were detected by only one of the technologies. In most of these cases, however, this finding is due to the thresholds set for categorization, and the sign of the fold changes correlated in both techniques. In a small number of cases (3.5%), a clear discrepancy was found. In most of these cases, the proteins were represented in more than one spot and only one of these showed results conflicting with the mRNA measurements. The spots may therefore correspond to degradation products or to covalent reaction intermediates of enzymes. Only one protein, phosphoglycerate kinase (HI0525), was clearly detected as repressed by 2D-gel comparison and induced by transcriptional imaging. Statistical analysis of the value pairs was performed and yielded correlation coefficients between 0.02 and 0.52 for the three time points, indicating a very weak degree of relatedness between the obtained values. These observations suggest that, in H. influenzae, the gene expression changes observed on the mRNA and protein level may be qualitatively similar, but that the magnitude of the change detected differs significantly between the two technologies.Proteins Represented by Multiple Spots
Forty-three proteins were represented by more than one spot on the 2D gels. These protein isoforms could reflect protein posttranslational modifications, degradation products, covalent reaction intermediates, or artifacts such as acrylamide adducts. Many of these proteins are enzymes known to form covalent intermediates with the respective reaction educt (glyceraldehyde-3-phosphate dehydrogenase, threonine synthase, malonyl transacylase, fructose-1,6,-bisphosphate aldolase, transaldolase). Not surprisingly, in most of these cases the expression of only one of these spots followed the mRNA levels, whereas the intensity of the other one, presumably the covalent reaction intermediate, remained unchanged. Three proteins present as isoforms were tRNA synthetases (Asp, Gly, Lys), suggestive of a covalent reaction intermediate or of proteolytic cleavage. Other proteins have not been reported to undergo covalent modification or to act by an unknown mechanism (e.g., hslUV, aspartase, phosphoenolpyruvate carboxykinase). The study of these protein isoforms, although not performed in detail in this work, could yield important information about protein modifications that affect protein function.The Effect of DNA Gyrase Inhibitors on the Gene Expression Pattern
We addressed the question of whether two different classes of DNA
gyrase inhibitors would induce different responses at the level of
transcription or of translation or of both. As described earlier, a key
step in the DNA supercoiling reaction is the DNA gyrase-mediated
double-stranded cleavage and resealing of DNA. It has been shown that
the two classes of subunit A inhibitors, the quinolones and the
pyrimido[1,6-a]benzimidazoles, interrupt this cleavage
and resealing cycle at the cleavage step (Hooper and Wolfson 1991
;
Hubschwerlen et al. 1992
; Hooper 1993
; Gmuender et al. 1995
). Two other
classes of gyrase inhibitors, cyclothialidines and coumarins, bind to
the ATP binding site located in the subunit B and inhibit the ATPase
activity (Ali et al. 1993
; Maxwell 1993
; Ali et al. 1995
; Nakada et al.
1995
; Lewis et al. 1996
; Oram et al. 1996
). These compounds therefore
inactivate the enzyme without introducing DNA strand breaks. Novobiocin
(a coumarin) and Ciprofloxacin (a quinolone) were chosen as
well-characterized DNA gyrase inhibitors representing these two
functional classes. To elucidate the concentration- and time-dependence
of the response, we used low inhibitor concentrations, approximating
the minimal inhibitory concentration (MIC, values derived from
conventional agar plate techniques), and 10-fold higher concentrations.
Cells were collected after 10, 30, and 60 min. Both inhibitors were
added at an OD of 0.4 (mid log phase). Control cells grown from the
same cultures in the absence of inhibitors were collected at the same
time. RNA from treated and from control cells was isolated, reverse
transcribed, and hybridized to microarrays. The expression patterns
from the inhibitor treated cells were compared with the corresponding
control cells after normalization by using the Affymetrix GeneChip
software. At the time of harvesting, an aliquot was removed from the
culture and the bacteria were pulse labeled with radioactive
methionine. The protein extract from this culture was used for the
quantification of protein synthesis rates by 2D-PAGE followed by
computerized spot quantification.
Novobiocin
Novobiocin is not a bactericidal DNA gyrase inhibitor but indirectly influences the optimal supercoiling. The degree of supercoiling can influence promoter activity (Jovanovich and Lebowitz 1987
|
|
|
Ciprofloxacin
Quinolones are bactericidal and provoke the so-called RecA (SOS) DNA repair system (Piddock and Wise 1987
|
|
|
Expression Changes Specific for Novobiocin or Ciprofloxacin
The response to Ciprofloxacin shows some clear differences when compared with the Novobiocin-induced response. As shown in Table 8, the induction of DNA repair systems distinguished the response of the cells to Ciprofloxacin from that to Novobiocin. Furthermore, the onset of the response was delayed for Ciprofloxacin and immediate for Novobiocin. The induction of these proteins could therefore represent a molecular marker to distinguish between the response to quinolones compared with other DNA gyrase inhibitors. Novobiocin had little effect at a low concentration. At the higher concentration, however, it affected many more genes than did Ciprofloxacin, suggesting that it exerts a stronger action on the cell at 10×MIC than does Ciprofloxacin. It is therefore difficult to determine genes that are specifically affected by Novobiocin. Detailed studies of concentration-dependent responses could clarify this issue.
|
Common Effects
To examine whether the commonly affected genes code for proteins belonging to a specific functional group, they were classified as increased, decreased, or not changed at the respective highest antibiotic concentration and at the 60-min time point and were ordered in functional groups according to the scheme proposed by Fleischmann et al. (1995)
|
Conclusions
In summary, the results show that the high-density microarrays yield highly reproducible results and that the main difficulties for reproducible analysis of low abundance transcripts lie in biological sampling. Our results show that if conditions are kept as reproducible as possible, most of the transcripts can be detected even when present at low concentrations. Our studies analyzing the number of transcripts called present after hybridization showed that in all experiments, up to 85% of transcribed genes could be detected. Transcripts showing very low intensity values, only a few percent above background fluorescence, could also be reproducibly identified as present, although the uncertainties of the signals in this fluorescence range is higher. Not surprisingly, the sensitivity and reproducibility of the expression analysis by using oligonucleotide chip technology was clearly better than expression analysis by using 2D-PAGE followed by computerized image analysis. The dynamic ranges were 103 for transcript analysis and 105 for the quantification of proteins. However, although qualitatively similar, there are some quantitative differences in the response detected by protein quantification compared with mRNA quantification. This highlights the importance of combining both technologies to obtain important information on the level (transcriptional or translational) at which the regulatory mechanisms act. Moreover, although relatively rare in bacteria, posttranslational modifications constitute an important additional level of regulation and can only be studied by proteome investigations. The detection of proteins present as multiple spots underscores this point.
Expression analysis by using the bacterial microarray system or 2D gels can be used to profile the effect of an inhibitor on a cell, but the main initial challenge is to discover the appropriate concentration and time point. In the described experiments a low concentration, around the MIC values, and a rather short incubation period that is within minutes, seem to provide the best results for the detection of the genes that are primarily affected. Because of the overwhelming amount of data it is more difficult to analyze the responses after incubation with a higher inhibitor concentration and/or after longer incubation times but, on the other hand, these results may help us analyze and understand more complex response patterns. Profiling the response of a selected inhibitor class may also give indications for a classification of an unknown inhibitor because its profile can be compared with those from known inhibitors. The example of Ciprofloxacin and Novobiocin illustrates that the response to an antibiotic can yield important information as to its mode of action. Both compounds induce the expression of DNA gyrase and negatively affect topoisomerase I expression. Ciprofloxacin, in addition, induces the SOS response. The two different modes of action are thus clearly reflected in the cellular response. These genes may therefore be useful indicators for gyrase inhibition by either mechanism. More experiments with different gyrase inhibitors would be required to substantiate these findings. Expression analysis will prove to be an invaluable tool not only for the study of disease processes but also for the characterization of novel pharmaceuticals.
| |
METHODS |
|---|
|
|
|---|
Cell Growth
H. influenzae Rd KW 20 was used as the model organism for
these experiments. Bacterial cultures were grown in minimal medium with
a reduced methionine concentration (0.6 µM) to an OD600
of 0.4 (Barcak et al. 1991
). The cultures were divided into aliquots of
300 mL and antibiotics were added to a final concentration of 30 and
300 ng/mL for Ciprofloxacin, and 12.5 and 125 ng/mL for Novobiocin,
respectively. Controls without antibiotic were grown in parallel. After
10, 30, and 60 min, aliquots were taken for metabolic labeling and RNA
extraction. To the first aliquot (1mL), 0.7 MBq of
L-[35S]Methionine (>37 TBq/mmol, Amersham
Radiochemicals) were added and incubation was continued for 2 min. The
cells were then rapidly chilled on ice, harvested by centrifugation,
and frozen at
20°C. For RNA preparation, a 35-mL aliquot was
spun in a chilled centrifuge and the cell pellet was snap-frozen in
liquid nitrogen and kept at
80°C. Five individual experiments
were performed, whereby aliquots for RNA extraction were only collected
for two experiments.
RNA Extraction, Preparation, Array Hybridization, and Scanning
Bacterial RNA was isolated, labeled, and hybridized to the chips essentially as described (de Saizieu et al. ). Before fragmentation of the biotin-labeled cDNA, an additional purification step was performed by using Chromaspin-100 columns (Clontech) and the fragmented cDNA was centrifuged quickly through 0.22-µM filter units (Millipore).
2D-PAGE
The cell pellets were washed once in PBS buffer (Life Technologies). The cells were then lysed by resuspension in sample buffer containing 8 M urea, 4% CHAPS, 40 mM of Tris base (Fluka), 65 mM of 1,4-dithioerythritol (Merck), and 2% ampholytes (Resolyte 3-10, BDH). The extracts were centrifuged at 100,000g and the supernatant recovered. The amount of incorporated radioactivity was determined in a Model 2500 TR liquid scintillation counter (Packard Instrument Co.).
Aliquots of the protein extracts containing 4 × 106 cpm
of radioactivity were loaded onto Immobiline 3-10 nonlinear pH
gradient strips (Pharmacia) at the basic end and resolved according to the manufacturer's recommendations. The strips were equilibrated as
described (Sanchez et al. 1995
) and loaded onto 1-mm thick vertical
12% polyacrylamide slab gels. After electrophoresis, the gels were
dried on 3 MM Whatman filter paper and exposed to PhosphorImager
screens (Molecular Dynamics). Images were analyzed by using PDQuest
software (BioRad). Parallel samples (one sample per time point and
concentration and their corresponding controls) were run on parallel
gels (same batch of strips, same isoeletric focussing run, same batch
of gels for SDS-PAGE, and same SDS-PAGE run). Only pairs of gels that
had been obtained under identical conditions were considered for analysis.
The 2D gels were matched and the data were stored in an Oracle database. The spot intensities were normalized so that the sum of all the spot intensities was equal for all gels. The spot intensity data were exported to Microsoft Excel for further analysis.
Data Handling and Analysis
For transcriptional imaging, the hybridization intensities were
processed by using Affymetrix GeneChip software. Pairwise comparisons
of hybridization intensities were performed and the results were
exported to Microsoft Access for further analysis. The parameters that
were used for analysis were partly adopted from and partly derived from
the output of the GeneChip software (Table 10).
|
The fold change values were averaged for all experiments performed in
duplicate or triplicate. Results were considered significant if the
averaged fold change was >2 or <-2, the standard deviation of
the averaged fold change <0.25, the difference call "D" or "I", the purity factor >0.9, and the sort score >1 or <
1.
The spot intensities from the 2D gel experiments were exported to Excel and the significance of the results was estimated by using the t-test for paired samples. When the obtained P values were lower than .05, the changes were considered significant. For the calculation of average fold change values, the values were first converted to change factor format. The averages were then reconverted to fold change format to improve comparability and clarity.
The publication costs of this article were defrayed in part by payment of page charges. This article must therefore be hereby marked "advertisement" in accordance with 18 USC section 1734 solely to indicate this fact.
| |
FOOTNOTES |
|---|
1 Corresponding author.
E-MAIL hans.gmuender{at}roche.com; FAX 41-61-688 27 29.
Article and publication are at www.genome.org/cgi/doi/10.1101/gr.157701.
| |
REFERENCES |
|---|
|
|
|---|
New aspects.
Drugs
3:
8-14.
A possible new target for antibacterial agents.
Biochem. Biophys. Res. Commun.
250:
1-4[CrossRef][Medline].Received August 1, 2000; accepted in revised form October 12, 2000.
This article has been cited by other articles:
![]() |
D. A. Ostrov, J. A. Hernandez Prada, P. E. Corsino, K. A. Finton, N. Le, and T. C. Rowe Discovery of Novel DNA Gyrase Inhibitors by High-Throughput Virtual Screening Antimicrob. Agents Chemother., October 1, 2007; 51(10): 3688 - 3698. [Abstract] [Full Text] [PDF] |
||||
![]() |
P. D. Rogers, T. T. Liu, K. S. Barker, G. M. Hilliard, B. K. English, J. Thornton, E. Swiatlo, and L. S. McDaniel Gene expression profiling of the response of Streptococcus pneumoniae to penicillin J. Antimicrob. Chemother., April 1, 2007; 59(4): 616 - 626. [Abstract] [Full Text] [PDF] |
||||
![]() |
N. G. Coldham, L. P. Randall, L. J. V. Piddock, and M. J. Woodward Effect of fluoroquinolone exposure on the proteome of Salmonella enterica serovar Typhimurium J. Antimicrob. Chemother., December 1, 2006; 58(6): 1145 - 1153. [Abstract] [Full Text] [PDF] |
||||
![]() |
C. Freiberg, N. Brunner, L. Macko, and H. P. Fischer Discovering Antibiotic Efficacy Biomarkers: TOWARD MECHANISM-SPECIFIC HIGH CONTENT COMPOUND SCREENING Mol. Cell. Proteomics, December 1, 2006; 5(12): 2326 - 2335. [Abstract] [Full Text] [PDF] |
||||
![]() |
M. A. Campos, P. Morey, and J. A. Bengoechea Quinolones sensitize gram-negative bacteria to antimicrobial peptides. Antimicrob. Agents Chemother., July 1, 2006; 50(7): 2361 - 2367. [Abstract] [Full Text] [PDF] |
||||
![]() |
E. Marrer, A. T. Satoh, M. M. Johnson, L. J. V. Piddock, and M. G. P. Page Global Transcriptome Analysis of the Responses of a Fluoroquinolone-Resistant Streptococcus pneumoniae Mutant and Its Parent to Ciprofloxacin Antimicrob. Agents Chemother., January 1, 2006; 50(1): 269 - 278. [Abstract] [Full Text] [PDF] |
||||
![]() |
M. D. Brazas and R. E. W. Hancock Ciprofloxacin Induction of a Susceptibility Determinant in Pseudomonas aeruginosa Antimicrob. Agents Chemother., August 1, 2005; 49(8): 3222 - 3227. [Abstract] [Full Text] [PDF] |
||||
|
|