INTRODUCTION

Top Abstract Introduction Materials and Methods Results Discussion References

Escherichia coli responds to oxidative stress by modifying the expression of many genes. Early studies using two-dimensional gels to analyze variations in protein expression have shown that the synthesis of more than 80 proteins is activated in response to oxidative stress (19). Some of these induced proteins were identified as possessing fundamental antioxidant functions, e.g., superoxide dismutase and catalase. The search for mutants with altered antioxidant defenses led to the isolation and characterization of pleiotropic regulators that operate as redox-regulated genetic switches (3, 20, 42, 43, 45). The best-characterized pleiotropic regulators of the antioxidant responses are the OxyR and SoxR proteins (36). Both proteins have the remarkable ability of directly transducing oxidative signals to genetic regulation. Both proteins are expressed constitutively in an inactive state and are transiently activated in cells under specific types of oxidative stress. The activation of the OxyR and SoxR proteins results in the transcriptional enhancement of sets of genes (regulons) whose products relieve the stress by eliminating oxidants and preventing or repairing oxidative damage (36).

SoxR is a member of the MerR family of metal-binding transcription factors, and it exists in solution as a homodimer, with each subunit containing a [2Fe-2S] cluster. In nonactivated SoxR, these clusters are in the reduced state and their oxidation activates SoxR as a powerful transcription factor (12, 16). Interestingly, SoxR can also be activated by nitric oxide (NO) by direct nitrosylation of the iron-sulfur clusters (11). The active (oxidized or nitrosylated) form of SoxR activates transcription of the soxS gene up to 100-fold. The soxS gene product, SoxS protein, belongs to the AraC/XylS family of DNA-binding transcription factors (3), but its activity seems to be regulated solely at the level of expression. Conventional analysis using limited proteomics and genetic approaches showed that SoxS activates the expression of 17 genes or operons. The known SoxS-activated genes are sodA (encoding Mn-superoxide dismutase), fpr (NADPH-ferredoxin oxidoreductase), micF (antisense RNA, repressor of OmpF translation), ribA (cyclic GMP hydrolase), inaA (unknown function), fldA and fldB (flavodoxins A and B), nfo (endonuclease IV), marRAB (multiple-antibiotic-resistance operon), nfsA (also called mdaA, a nitroreductase), zwf (glucose-6-phosphate dehydrogenase), fur (an iron-binding repressor of iron uptake), fumC (fumarase C), acnA (aconitase), tolC (outer membrane protein), acrAB (drug efflux pump), and rimK (a modifier of ribosomal protein S6). Activating tolC, acrAB, micF, and rimK alters the sensitivity of E. coli and Salmonella enterica to a broad range of antibiotics (9, 10, 27, 33). SoxS is also a repressor of the soxS gene (34) and thus limits its own synthesis.

The diversity of genes activated by OxyR and SoxR illustrates the variety of cellular defense mechanisms against oxidative stress. Antioxidant mechanisms include the scavenging of reactive species (sodA, ahpCF), synthesis of reducing species (acnA, zwf), repair of oxidative damage (nfo, fpr), drug efflux (acrAB, tolC), reduction of cell permeability (micF), and replacement of redox-sensitive isozymes by redox-resistant isozymes (fumC). This variety is hardly surprising, given the large number of targets for oxidative damage; virtually all biological macromolecules can be damaged by oxidants. Particularly sensitive are electron-rich moieties, such as metal centers in proteins, unsaturated bonds in phospholipids, aromatic amino acids, and the double bonds of bases in nucleic acids (40). Oxidative stress has other, indirect effects, such as depletion of reducing power by consumption of NADH and NADPH in antioxidant reactions (40).

Interestingly, the induction of the marRAB regulon also results in enhanced resistance to oxidative agents and multiple antibiotics (3, 18). The genes of the marRAB regulon overlap significantly with those of the soxRS regulon (2, 5). This overlap evidently results from the structural similarity between the MarA and SoxS proteins and their respective DNA-binding sites (29, 30). Thus, two overlapping sets of genes are modulated by different signals sensed by different regulatory circuits. While the soxS gene is under the redox-regulated, positive control of SoxR, marA is under negative control by MarR, a repressor whose DNA-binding activity is regulated by the binding of small molecules with toxic effects (1, 31, 39). Among marRAB inducers are sodium salicylate (NaSal), the naphthoquinones menadione and plumbagin, and dinitrophenol. Although NaSal does not activate the soxRS regulon and although paraquat (PQ), a superoxide-generating agent, only marginally activates the marA regulon (30), the naphthoquinones activate both regulons (33). In contrast to PQ, menadione is a natural plant product. This fact has led to the hypothesis that the early evolution of the soxRS and marA regulons was shaped by the environmental stress mediated by naphthoquinones and other noxious xenobiotics from natural sources (33).

Although a complete profile of the cellular responses to oxidative stress has been lacking, the availability of the complete sequence of the E. coli genome now provides important tools for analyzing gene expression. The analysis of the variations of the "transcriptome," or transcriptional profiling, has already yielded abundant biological information in several organisms, including E. coli (4, 37, 44), Bacillus subtilis (14), Caulobacter crescentus (25), and Saccharomyces cerevisiae (15, 23).

Here we have determined genome-wide transcriptional profiles of E. coli cells exposed to the superoxide-producing agent PQ or to NaSal. We have also begun to dissect regulatory hierarchies within these responses by expressing, in the absence of stress, individual transcription factors that respond to environmental insult.

	MATERIALS AND METHODS

Top Abstract Introduction Materials and Methods Results Discussion References

Bacterial strains and plasmids. E. coli strain GC4468 (K-12 lacU169 rpsL) (19) was used as the wild type in all experiments. The GC4468 soxRS derivative DJ901 (20) transformed with pJP105 (placZ::soxS lacI Ap^r) (35) was used in the SoxS expression experiment. The GC4468 marRAB derivative MB106 was constructed for this study by transduction of the kanamycin-tagged 1.24-kb deletion spanning the marRAB operon from strain AG100/Kan (28) and was transformed with pMB102 (placZ::marA lacI Ap). Plasmid pMB102 was constructed by replacing the robA gene in plasmid pMB101 (6) with a PCR product containing the marA gene. The sequence of the cloned marA gene was confirmed at the Molecular Biology Core Facility of the Dana-Farber Cancer Center.

Culture growth and RNA isolation. Overnight cultures were diluted 1/100 in 15 ml of Luria broth contained in 125-ml Erlenmeyer flasks and were grown at 37°C and 250 rpm to an optical density at 600 nm (OD₆₀₀) of 0.5. At this point, the cultures were either left untreated or exposed to the different inducers for 45 min. The cells were harvested by centrifugation of 1.5-ml aliquots in a microcentrifuge for 30 s. An RNAeasy kit (Qiagen) was used to lyse the cells and extract total RNA. The RNA preparation obtained from the RNAeasy columns was extracted with phenol (pH 5.0) at 60°C, followed by extraction with a mixture of acidic phenol and chloroform and finally with chloroform alone. The RNA was precipitated with 90% ethanol and 3 M sodium acetate overnight at 20°C and was recovered by centrifugation for 10 min at 4°C. The RNA pellet was washed twice with cold 70% ethanol, dried in a Speed-Vac, and resuspended in diethylpyrocarbonate-treated water. The RNA was quantified by light absorption at 260/280 nm. Samples of the RNA were separated by electrophoresis in agarose gels and stained with ethidium bromide to verify the purity and integrity of the RNA preparation, using the 16S and 23S ribosomal bands as indicators. The absence of chromosomal DNA was verified by using the RNA preparations as the template in PCR amplifications. Briefly, ~100 ng of total RNA or ~20 ng of chromosomal DNA was incubated in PCR cocktails containing sodB-specific primers. While the reactions containing chromosomal DNA consistently yielded a DNA fragment of the expected size, the reactions containing RNA preparations always failed to yield a product.

cDNA synthesis. The RNA preparations from each sample were used as the template for cDNA synthesis by employing a commercial set of 4,290 open reading frame (ORF)-specific oligonucleotides (Sigma-Genosys) as primers. A total of 1 µg of total RNA was incubated in reaction buffer with the ORF-specific oligonucleotides, dTTP, dATP, and dGTP, at 90°C for 2 min and was then cooled to 42°C at a rate of 2.4°C/min. Two hundred units of avian myeloma virus reverse transcriptase (Sigma-Genosys) was added to each reaction, together with 20 µCi of [³³P]dCTP (3,000 Ci/mmol). The mixture was incubated at 42°C for 2 h. Labeled products were purified on gel filtration columns (Sephadex G-50).

Hybridization. The Panorama gene arrays (Sigma-Genosys) are positively charged nylon membranes onto which ORF-specific PCR products have been spotted in duplicate. The arrays represent the complete set of known and predicted ORFs as deduced from the complete genomic sequence of the E. coli K-12 strain MG1655 (7). For all experiments, a pair of membranes was used: one was hybridized with cDNA synthesized from the untreated cell culture, and the other was hybridized with the cDNA synthesized from the treated cell culture. Hybridization and washing steps were carried out following the manufacturer's instructions. The nylon filters were prehybridized at 65°C for 1 h in cylindrical tubes containing 5 ml of hybridization solution (5× SSPE [1× SSPE is 0.18 M NaCl, 10 mM sodium phosphate, and 1 mM EDTA {pH 7.7}], 2% sodium dodecyl sulfate, 1× Denhardt's reagent, and 100 µg of sonicated salmon testes DNA/ml). The whole cDNA preparation was first denatured at 100°C in 3 ml of hybridization solution, and the prehybridization solution was discarded from the tubes and replaced by the mix containing the labeled cDNA. The filters were hybridized for ~16 h at 65°C in a rotary oven. After this, the filters were rinsed with washing solution (0.5× SSPE, 0.2% sodium dodecyl sulfate) twice at room temperature and three more times at 60°C. The filters were then air dried and wrapped in clear plastic film. Two series of experiments were performed using two different pairs of membranes. For each experimental series, the hybridized membranes were stripped and reprobed up to four times, following the procedures recommended by the manufacturer.

Array imaging and analysis. Phosphorimaging screens were exposed to the hybridized filters for 48 to 72 h at room temperature. The screens were scanned in a Storm 840 PhosphorImager (Molecular Dynamics) at a 50-µm resolution. The resulting files were analyzed by determination of pixel density using Arrayvision software, which determined the intensity of each duplicate spot, measured in arbitrary units. The background signal was determined for each filter by averaging the intensities of 42 spots that did not contain DNA. This average background was then subtracted from the intensity at each DNA spot, and the corrected intensity of each spot was expressed as a percentage of the sum of all intensities. This treatment allowed comparison between filters independently of total hybridization intensity. The corrected intensities of duplicate spots were averaged, and the untreated and treated intensities from two independent experiments were averaged. The expression ratio for each gene was calculated as the treated/untreated intensities. Thus, an expression ratio of 1 indicated an invariable level of a transcript, whereas expression ratios larger or smaller than 1 indicated up- or down-regulated levels of mRNA, respectively.

Gene annotation. The expression values and ratios for each gene were transferred to Excel spreadsheets for statistical analysis and integration into updated, annotated databases (38), accessible online (http://genprotec.mbl.edu).

Web access. The complete data sets for all experiments are available online (http://www.hsph.harvard.edu/demplelab/genomics).

Gene probe synthesis and Northern blot analysis. The RNA samples used in the gene array experiments were also used in Northern blot experiments as an independent way to assess the quantitative validity of transcriptional profiling. Probes for specific genes were generated by PCR amplification using chromosomal DNA from strain GC4468 as template and gene-specific primers (ORF-mers) obtained from Sigma-Genosys. Typically, PCR amplifications were carried out in 30 cycles of annealing at 60°C (45 s), elongation at 72°C (1 min), and denaturation at 94°C (30 s). The PCR products were resolved by electrophoresis in 1.25% agarose gels, recovered by excision from the gel, and purified using Qiaquick DNA-binding microspin columns (Qiagen). The DNA fragments were labeled by nick translation using a Klenow DNA polymerase fragment, random hexamers (Gibco BRL), and [³²P]dCTP (3,000 Ci/mmol) plus unlabeled dATP, dGTP, and dTTP. The labeled probe was purified by gel filtration in Sephadex G-25 columns (Pharmacia). For the Northern blot experiments, 2 µg of total RNA was run per lane in 1.25% agarose gels containing formaldehyde and was transferred to Nytran membranes using a Turboblotter setup (Schleicher & Schuell). The RNA was cross-linked to the membrane by UV irradiation, and the membranes were then hybridized at 65°C with radioactively labeled DNA fragments in cylindrical tubes using QuickHyb solution (Stratagene). The membranes were washed according to the instructions from the manufacturer. X-ray films were exposed to the membranes at 70°C and were developed using a Fuji automatic developer. The radioactive signals were measured using an Applied Biosystems phosphorimager.

	RESULTS

Top Abstract Introduction Materials and Methods Results Discussion References

A transcriptional profile of E. coli exposed to PQ. The genomic transcriptional profile of E. coli growing on Luria broth under superoxide stress was determined by comparison of cultures in which one sample was left untreated while the other was exposed to 250 µM PQ for 45 min. This concentration of PQ induced the genes of the soxRS regulon (35) but failed to inhibit growth rate in any substantial way during the 45-min exposure. We expected that these exposure conditions would prevent the activation of "general stress" pathways that generally accompanies growth inhibition (22) and would thus maximize the chance to identify specific, oxidative stress-responsive genes. Figure 1 shows the growth of the untreated and PQ-treated cultures and demonstrates the lack of effect of the chosen PQ concentration on growth during the 45-min exposure.

View larger version (12K): [in this window] [in a new window]   FIG. 1.   Growth of E. coli strain GC4468. OD600 as a function of time for strain GC4468 was measured. The cultures were either left untreated (empty circles) or treated with 250 µM PQ (filled circles). PQ was added to log-phase cultures (empty arrow), and cells were harvested after 45 min (filled arrow). For detailed culture conditions, see Materials and Methods.

View larger version (17K): [in this window] [in a new window]   FIG. 2.   Scatter plot of expression levels for the E. coli genome: untreated and PQ-treated cells. The percent total intensity for each gene represented in the arrays is plotted on a log scale. A relatively small number of genes had associated signal intensities that were below average background levels. The subtraction of the background value from these signals resulted in a negative corrected value and therefore a negative percentage of total signal. These genes were given an arbitrary value of 0.000001% of total intensity.

Validation by Northern blotting of expression ratios from transcriptional profiling. In order to test independently the values for the expression ratios obtained from the gene arrays, we performed Northern blot analysis for 21 genes. The same RNA samples used as templates in one of the PQ-treatment experiments were run on an agarose gel, transferred to nylon membranes, and hybridized with labeled gene probes (see Materials and Methods). The intensity of each signal was measured by phosphorimaging, and the expression ratio (treated/untreated intensities) was calculated for each gene. Figure 3A shows the comparison of the expression ratios obtained from gene array experiments and conventional Northern blots. The 21 genes selected for the comparison included 11 genes that scored as activated in the gene array experiment, 7 genes that scored as unmodified (expression ratios between 0.5 and 2), and 3 genes that scored as down-regulated. In general, the results of the two methods were similar across at least 2 orders of magnitude. However, the correlation between the Northern blot and the array values was weaker for genes with basal (or down-regulated) expression levels that were close to background. Figure 3B shows three examples of genes (nfo, uraA, and yaiA) with comparatively large differences between the expression ratios obtained by transcriptional profiling and Northern blotting. In all three cases, the basal or down-regulated level of mRNA was extremely low as estimated by Northern blotting. In contrast, genes with detectable basal levels (sodA, gshB, and ptsG) show good correlation in their expression ratios measured by the two methods (Fig. 3A).

View larger version (15K): [in this window] [in a new window]   FIG. 3.   Comparison of expression ratios from transcriptional profiling and Northern blots. (A) Total RNA was extracted from untreated or PQ-treated cultures (250 µM for 45 min). Aliquots from these preparations were used as a template for cDNA synthesis and hybridization with gene arrays or were run in agarose gels, transferred to Nytran membranes, and probed with labeled gene-specific PCR fragments. The expression ratio of 21 genes is shown for Northern blotting (horizontal axis) or cDNA synthesis and hybridization to gene arrays (vertical axis). The genes tested were cyoD, cysK, dnaE, gshB, inaA, lpxC, nadE, nfo, nupC, pdhR, ptsG, pyrB, rpsS, sodA, speE, uraA, yadJ, ybjC, yhiM, and zwf. See Materials and Methods for detailed protocols. (B) Northern blots of the nfo, yaiA, uraA, sodA, gshB, and ptsG genes. , absence of PQ; +, presence of PQ. A replicate gel was run and stained with ethidium bromide (EtBr), revealing the 16S and 23S rRNA, which serves as loading control.

A search for novel SoxS-regulated genes. Although transcriptional profiling provides information about relative RNA levels, it does not establish regulatory hierarchies among genes. In an effort to begin dissecting the regulatory cascade within the response to PQ, we expressed SoxS protein in the absence of oxidative stress. The plasmid pJP105 expresses the SoxS protein from an isopropyl--D-thiogalactopyranoside (IPTG)-regulated promoter and has been used to identify SoxS-regulated operon fusions (35). The activity of the SoxS protein is regulated exclusively by its intracellular concentration, and thus the artificial induction of soxS expression from pJP105 should be sufficient to modulate all genes of the soxRS regulon. The E. coli strain DJ901 (soxRS) was transformed to Ap^r with plasmid pJP105, and cultures were either left untreated or were treated with IPTG. The labeled cDNA from these cultures was hybridized to Panorama gene arrays, and the blots were analyzed using the same stringency conditions applied to the PQ experiment. Again, the vast majority of genes was unaffected, with only 95 genes modulated beyond the statistical threshold: 37 genes were activated and 58 genes were down-regulated. The activated genes included 11 out of the 16 known SoxS-activated genes represented in the gene array (data not shown). The indirect down-regulation of ompF was also registered. Of the 37 SoxS-activated genes revealed by transcriptional profiling, 14 were also registered in the PQ exposure experiments. The genes that were activated both by PQ and by expression of SoxS are indicated in Table 1 in boldface. The complete results from the SoxS-expression experiment can be obtained online (http://www.harvard.hsph.edu/demplelab/genomics).

A transcriptional profile of the E. coli response to NaSal. NaSal is a natural product with a role in signal transduction in the response to plant infection (8). In the millimolar concentration range, NaSal dissipates the proton gradient across the inner membrane, chelates iron, inhibits growth, and induces the heat shock and marA regulons (41), the latter by binding to the repressor of the marRAB operon (31, 39). In order to characterize the global response to NaSal and to identify novel regulatory overlaps with the response to PQ, the transcriptional profile of growing cultures of E. coli treated with NaSal was determined.

View larger version (13K): [in this window] [in a new window]   FIG. 4.   Growth of E. coli GC4468 exposed to NaSal. OD600 as a function of time for strain GC4468 was determined. The cultures were either left untreated (empty circles) or treated with 5 mM NaSal (filled circles). NaSal was added to log-phase cultures (empty arrow), and cells were harvested after 45 min (filled arrow). For detailed culture conditions, see Materials and Methods.

View larger version (18K): [in this window] [in a new window]   FIG. 5.   Scatter plot of expression levels for the E. coli genome: untreated and NaSal-treated cells. The percent total intensity for each gene represented in the arrays is plotted on a log scale. A relatively small number of genes had associated signal intensities that were below average background levels. The subtraction of the background value from these signals resulted in a negative corrected value and therefore a negative percentage of total signal. These genes were given an arbitrary value of 0.000001% of total intensity.

A search for MarA-regulated genes. To begin dissecting the regulatory cascade of the response to NaSal, we expressed MarA in the absence of exogenous toxic agents. The plasmid pMB102 harbors the marA gene under the control of an IPTG-regulated promoter. As with SoxS, the activity of the MarA protein is regulated exclusively by its intracellular concentration (2). Thus, the artificial induction of marA from pMB102 should be sufficient to modulate all MarA target genes. Cultures of an E. coli marRAB strain containing plasmid pMB102 were grown and either left untreated or treated with IPTG. After 45 min, total RNA was extracted from both cultures and analyzed by hybridization to gene arrays. As for the previous experiments, background-corrected expression values were determined and used to calculate the expression ratios for each gene and the same significance thresholds were applied. In total, 88 genes were modulated by MarA expression; 67 genes were activated and 21 genes were down-regulated. The activated genes include 21 of the 62 postulated MarA-regulated genes (5) (data not shown). An additional 19 genes had similar expression trends as reported but fell below our statistical threshold. Those genes activated both by exposure to NaSal and by expression of MarA, 20 in total, are indicated in boldface in Table 2. The complete results from the MarA-expression experiment can be obtained online (http://www.hsph.harvard.edu/demplelab/genomics).

The overlap between the PQ and NaSal stimulons. Sixteen genes were modulated with the same trend by both PQ and NaSal (Table 3). Some genes were known from standard genetic studies to be regulatory targets for SoxS and MarA (fumC, inaA, marA, and sodA). In addition, the PQ- or SoxS-responsive genes gatABD, gltA, nfnB, and ybjC were recently shown to be activated by constitutive expression of MarA (5). The ybjC gene was also activated by expression of SoxS in our experiments (Table 1). This study identified five additional genes activated by both PQ and NaSal. These genes code for products involved in arginine transport (artP), cysteine synthesis (cysK), protection of DNA from iron-mediated oxidative damage (dps), and salvage of nucleotides (deoB). One other gene, activated by both PQ and NaSal, has no known function or extensive homology to any other gene of known function (b1452). Only one gene was down-regulated by both treatments: pyrB, which codes for aspartate transcarbamylase, involved in pyrimidine biosynthesis.

	DISCUSSION

Top Abstract Introduction Materials and Methods Results Discussion References

We have used a functional genomics approach to identify novel genes that respond to oxidative stress. In all the experiments, our results not only bear the hallmarks of the cellular responses that we intended to evoke but also provide new insights in the physiology of oxidative stress. The resulting tally of the genes involved in these responses should be considered an underestimate for reasons that are inherent to the method used. First, commercially available gene arrays are a fixed platform that do not admit modifications. The arrays used excluded untranslated RNAs, which are clearly involved in responses to changing redox conditions (46-48). Second, transcriptional profiling reveals only comparative, steady-state levels of mRNAs, without any information about posttranscriptional processes or actual protein expression. Third, despite the almost complete coverage of the genome, transcriptional profiling experiments consistently fail to detect changes in genes known to be modulated by the stimuli of interest. For example, the use of gene arrays to analyze the heat shock response revealed only 23 of the 51 known genes (37). Finally, the results can vary significantly between transcriptional profiling experiments. We have tried to address this problem by averaging replicate experiments and by establishing statistical thresholds for the expression ratios beyond 99% confidence. For those experiments involving induction of the SoxS and MarA transcription factors, we performed a single experiment. However, we included in our analysis only those genes that were also modulated in the duplicate experiment involving the corresponding global inducer.

PQ activates genes involved in pathways that reconstitute NADH and NADPH pools. PQ is reduced intracellularly at the expense of NADPH in a reaction catalyzed by at least three oxidoreductases (26). PQ reduced by one electron is oxidized by O₂ to form superoxide, resulting in a redox cycle that produces a flux of superoxide. Thus, the cell faces a double threat under exposure to PQ and to other redox-cycling agents: the deleterious effects of superoxide itself and the decreased level of NADPH that limits biosynthetic capabilities. Equilibration of NADPH with NADH would generalize this limitation of cellular reducing power.

Evidence for adaptation and repair pathways under oxidative stress. Challenge of growing E. coli cells with PQ induces the genes coding for nine ribosomal proteins (Table 1). This increase in ribosomal building blocks was not predicted, in that the bacterial growth rate was not significantly affected during the 45-min treatment with PQ. Another gene coding for a translational regulator, fmt, is also induced by PQ. The product of this gene, methionyl-tRNA formyltransferase, is an important factor in translational initiation (21). One possibility is that increased translational capacity counterbalances a faster turnover of proteins due to oxidative damage and increased degradation. In this scenario, enhanced synthesis would be required to maintain the high growth rate. The possibility of an increased metabolic rate under superoxide stress is consistent with the activation of genes coding for products involved in crucial anabolic and catabolic pathways. These genes include nuoI and nuoK coding for subunits of NADPH dehydrogenase, the first electron acceptor of the respiratory chain. Interestingly, 11 out of the 12 genes coding for NADPH dehydrogenase subunits showed some degree of activation in both PQ exposure experiments, albeit below the statistical threshold levels.

Common themes in the responses to PQ and NaSal. The list of genes activated by treatment of E. coli with PQ or NaSal (Table 3) is a first approximation to the common solutions to the physiological challenges posed by these two compounds. In addition to the common genes listed in Table 3, many other genes activated by the individual stresses have comparable functions. For example, the induction of genes involved in sugar transport occurs in both situations, albeit with specificity for different sugars. Exposure to both PQ and NaSal results in activation of the galactitol (gat) operon, while exposure to PQ activates the additional glucose and maltose transport genes (Table 1) and exposure to NaSal activates genes for the transport of sorbitol and mannose (Table 2).

From statistical significance to biological relevance. Transcriptional profiling can measure changes in only the steady-state levels of mRNA, and the significance of these changes can be estimated by statistical analysis. However, the biological relevance of the detected variations in mRNA levels must be substantiated by genetic and biochemical analysis of the functions. We have constructed deletion mutants of the genes cysK and b2962 and measured the sensitivity of these strains to PQ and other oxidants. Preliminary results show that both deletion mutants are hypersensitive to PQ but not to H₂O₂ or to tert-butyl hydroperoxide (data not shown). These results suggest that even relatively small increases in gene expression, as in the cases of cysK and b2962, can point to genes with important roles in defense against oxidants.