Previous Article | Next Article ![]()
Journal of Bacteriology, November 2005, p. 7655-7666, Vol. 187, No. 22
0021-9193/05/$08.00+0 doi:10.1128/JB.187.22.7655-7666.2005
Copyright © 2005, American Society for Microbiology. All Rights Reserved.
Department of Chemistry, Williams College, Williamstown, Massachusetts 01267,1 Department of Biology, Massachusetts Institute of Technology, Cambridge, Massachusetts2
Received 17 June 2005/ Accepted 30 August 2005
|
|
|---|
|
|
|---|
By contrast with E. coli, where more than 30 SOS genes have been identified (7, 8), only 5 B. subtilis SOS genes have been shown to meet both SOS gene criteria thus far: recA, lexA, uvrB (formerly dinA), dinB, and dinC (also called tagC) (4, 9, 15, 25). The E. coli uvrB gene encodes part of the UvrABC endonuclease, which catalyzes nucleotide excision repair of a variety of DNA lesions (41). Both the E. coli uvrA and uvrB genes are damage-inducible SOS genes with LexA binding sites overlapping their promoters (8); the E. coli uvrC gene is not damage inducible (32), and LexA does not bind to its promoter region in vitro (12). Homologs of all three uvr genes are present in B. subtilis, and genetic and biochemical evidence indicate that the uvrA, uvrB, and uvrC genes are involved in excision repair (14, 23). The functions of the dinB and dinC genes are unknown. Because it is adjacent to the tag operon (which codes for enzymes involved in teichoic acid synthesis), the dinC gene has been named tagC; however, there is no evidence for its involvement in teichoic acid synthesis. We will refer to it as dinC in this report.
We report here the identification of 28 additional SOS genes in B. subtilis. They all have LexA binding sites, or SOS boxes, within their promoter regions, and they are induced by at least two distinct DNA-damaging treatments in RecA+, but not recA null, cells. The consensus operator sequence for the binding of a B. subtilis LexA dimer, 5'-CGAACN4GTTCG-3', has been characterized by mutational analyses and DNA binding studies (4, 5, 27, 53). A study of LexA binding to recA operator mutants suggests the requirements for site-specific LexA binding summarized in Fig. 1 (E. S. Groban, N. Au, M. B. Johnson, P. Banky, P. G. Burnett, G. L. Calderon, E. C. Dwyer, S. N. Fuller, B. Gebre, L. M. King, I. N. Sheren, L. D. Von Mutius, T. M. O'Gara, and C. M. Lovett, submitted for publication). According to the study, the thermodynamically preferred half site sequence for LexA binding is 5'-CGAACAT-3'; certain substitutions do not reduce binding affinity significantly, while others (labeled destabilizing replacements) abolish binding altogether. Guided by these binding requirements, we searched the genome for sites within putative promoter regions that could potentially bind LexA. We assessed binding activity using mobility shift assays and we identified genes that show RecA-dependent induction by DNA damage using genomic microarrays.
![]() View larger version (15K): [in a new window] |
FIG. 1. Sequence requirements for LexA binding. The preferred half site sequence based on a thermodynamic analysis of LexA binding to recA operator mutants. Base substitutions labeled as destabilizing abolish LexA binding to the recA operator (Groban et al., submitted).
|
|
|
|---|
Preparation of promoter regions for mobility shift assays.
DNA containing putative SOS operators was prepared by PCR amplification of B. subtilis YB886 (59) DNA (10 ng/ml) using synthetic oligonucleotide primers (2 µM) with a Peltier PTC-200 thermal cycler (MJ Research). Samples of amplified DNA were electrophoresed alongside DNA samples of known concentration; gels were analyzed by densitometry with an Alpha Innotech imaging system, and the concentration of amplified DNA was interpolated from DNA standard curves. The promoter regions prepared by PCR amplification were radiolabeled with [
-32P]ATP using T4 kinase. Radiolabeled DNA was purified by electrophoresis on an 8% nondenaturing polyacrylamide gel.
Mobility shift assays. For competition experiments, purified LexA was incubated with radiolabeled recA promoter DNA (5 to 10 nM) and a 5- to 50-fold molar excess of competitor DNA for 30 min at 25°C in mobility shift buffer, which consisted of 12 mM HEPES-NaOH (pH 7.9), 4 mM Tris-Cl (pH 7.9), 12% glycerol, 60 mM KCl, 1 mM EDTA, 1 mM dithiothreitol (DTT), 2 µg poly(dI-dC) · poly(dI-dC), and 0.3 mg/ml bovine serum albumin. This incubation mixture (10 µl) was loaded on a 4% (acrylamide:bisacrylamide ratio of 80:1) nondenaturing polyacrylamide gel, and electrophoresis was begun immediately. The buffer within the gel and the running buffer were both 25 mM Tris-Cl (pH 8.5), 250 mM glycine, and 1 mM EDTA. Samples were electrophoresed, and the dried gel was subjected to densitometric analysis using a Bio-Rad Molecular Imager FX phosphorimaging system.
For determination of the apparent binding constant, Kd, purified LexA was incubated with radiolabeled promoter DNA (5 to 10 nM) for 30 min at 25°C in mobility shift buffer, and 10 µl of this incubation mixture was subjected to electrophoresis and phosphorimaging analysis as described above. For Kd determinations using promoter fragments as competitive inhibitors of LexA binding to the recA operator, purified LexA was incubated with radiolabeled recA promoter DNA (5 to 10 nM) and competitor DNA (50 nM) for 30 min at 25°C in mobility shift buffer and analyzed as described above.
Microarray analyses. DNA microarrays were prepared using PCR products from >99% of the annotated B. subtilis open reading frames spotted onto Corning GAPS II slides, essentially as described previously (16). Prior to hybridization with biological samples, arrays were prehybridized for at least 45 min at 42°C in 1% bovine serum albumin, 5x SSC (1x SSC is 0.15 M NaCl plus 0.015 M sodium citrate), 0.1% sodium dodecyl sulfate (SDS), washed in water, and dried (16).
Cultures (25 ml) of strain YB886 (metB5 trpC2 xin-1 SPß amyE sigB) or YB3000 (YB886 recA260) were grown in defined minimal medium to an optical density at 600 nm of 0.3, treated with 1 µg/ml mitomycin C (MC) or 25 mJ/m2 UV radiation, incubated for 1 h at 37°C, and harvested by centrifugation. For UV treatment, cells were pelleted, resuspended in SMS minimal medium (16a), transferred to sterile petri plates, and exposed to a 254-nm germicidal lamp to a dose of 25 mJ/m2. Survival at this dose was between 20% and 50%. Untreated cells were handled similarly except they were not UV irradiated. Culture samples were immediately mixed with an equal volume of methanol (prechilled to 20°C). Samples were then spun to pellet the cells, the supernatant was discarded, and cell pellets were frozen at 80°C until further use. RNA was extracted using QIAGEN RNeasy kits combined with on-column DNase treatment according to the manufacturer's recommendations. RNA was then precipitated with ethanol and LiCl and resuspended in RNase-free water to a desired concentration, usually
1.0 mg/ml. The quality of RNA was checked on agarose gels by visualizing the integrity of the 23S and 16S rRNA.
To generate cDNA, RNA from the different experimental conditions was reverse transcribed in the presence of amino-allyl-dUTP, followed by coupling to Cy5 for all experimental samples or Cy3 for all reference RNAs. Reference RNA was made by pooling RNA samples from various strains grown under conditions similar to the experimental conditions. For reverse transcriptase reactions, 10 µg RNA template was mixed with 2.5 µg random hexamers (in 18 µl) and incubated at 70°C for 10 min and on ice for 5 min. Reverse transcription reactions were then started by the addition of a cocktail resulting in a final mix of RNA template, random hexamer primers, and 300 U SuperScript II reverse transcriptase, 1x reverse transcription buffer, 10 mM DTT, and deoxyribonucleoside triphosphates (0.5 mM each dATP, dCTP, and dGTP; 0.1 mM dTTP; 0.4 mM aminoallyl-dUTP) in a final volume of 30 µl. The labeling reaction mixtures were incubated at 25°C for 10 min and at 42°C for 70 min and then shifted to 70°C for 15 min to stop the reactions. RNA in the reaction mixtures was degraded by adding NaOH (33 mM final concentration) and incubating at 70°C for 10 min. HCl (33 mM) was added to each reaction mixture to neutralize the pH.
Reaction mixtures were purified with QIAGEN MinElute kits and eluted in 10-µl volumes, and 0.5 µl of 1 M NaHCO3 (pH 9.0) was added to adjust the pH for the coupling reactions. To couple the fluorescent dyes to cDNA, 1 µl freshly dissolved Cy3 or Cy5 dye (Amersham) was added to cDNA and incubated for 1 h in the dark, mixing every 15 min. Reactions were quenched by incubation with 1.4 M hydroxylamine for 15 min. Each experimental (Cy5-labeled) sample was mixed with an aliquot of reference RNA (labeled with Cy3), and mixed samples were purified with QIAGEN MinElute kits. The labeled samples were mixed with 10 µg salmon sperm DNA and 0.8 µg yeast tRNA, and the volume was adjusted to 14 µl. The samples were heated to 100°C for 5 min, spun down, mixed with 2x hybridization buffer (0.05% SDS, 5x SSC, 25% formamide final concentration) and hybridized to DNA on a microarray for at least 16 h at 42°C. Following hybridization, arrays were washed with 1x SSC-0.2% SDS for 5 min at 42°C, followed by a 5-min wash with 0.1x SSC-0.2% SDS at room temperature, and a final 5-min wash in 0.1x SSC at room temperature. Arrays were spun to remove extra liquid and dried with nitrogen gas.
Arrays were scanned and analyzed with GenePix 3.0 software (Axon Instruments, Inc.). The signal intensity for each spot (gene) [(Cy5/Cy3)g] was normalized to the total signal intensity [(Cy5/Cy3)t] on the array (essentially the sum of all the spots). Normalized ratios of experimental RNA abundance over reference RNA was obtained for each spot for which 80% of the pixels had intensities at least 1 standard deviation above background. In an average experiment, 96% of all genes gave such a signal. For genes of interest, we verified that the intensities of 80% of the pixels in the spot were at least 2 standard deviations above background. To compare two conditions, the normalized signals from condition A (e.g., cells treated with MC) were divided by the corresponding signals from condition B (cells not treated with MC): [(Cy5/Cy3)g/(Cy5/Cy3)t] under condition A/[(Cy5/Cy3)g/(Cy5/Cy3)t] under condition B.
Since all experiments were done at least in triplicate, we report the average ratio of ratios from all repetitions. Statistical analysis of microarrays (48) was used for all experiments. The input data contained the replicate ratios of experimental over reference samples in log2 format. Thus, for each gene, at least six values were input for calculation of significancethree from one condition and three from the other. For all experiments, an effect on a gene is considered statistically significant if there is less than 1% probability that this change occurred by chance (i.e., a false discovery rate of less than 1.0%).
|
|
|---|
|
View this table: [in a new window] |
TABLE 1. B. subtilis genes with LexA binding sites in their putative promoter regions
|
|
View this table: [in a new window] |
TABLE 2. Canonical SOS boxes within promoter regions that are not bound by LexA
|
A promoter sequence.
![]() View larger version (54K): [in a new window] |
FIG. 2. Binding of B. subtilis LexA to potential SOS promoters. Mobility shift assays were conducted with purified LexA, radiolabeled recA promoter DNA (5 to 10 nM), and a 5- to 50-fold molar excess of the indicated promoter DNA as described in Materials and Methods. The lower and upper bands correspond to unbound and LexA-bound recA promoter DNA, respectively. Lanes with no LexA protein or competitor DNA added are indicated.
|
![]() View larger version (18K): [in a new window] |
FIG. 3. Binding of B. subtilis LexA to the recA promoter. Graphical analyses of mobility shift titration of 32P-labeled recA promoter (10 nM) incubated with increasing concentrations of LexA as described in Materials and Methods.
|
![]() View larger version (78K): [in a new window] |
FIG. 4. Binding of B. subtilis LexA to the lexA and yqjW promoters. Mobility shift assays were conducted with purified LexA (0 to 96 nM) and radiolabeled lexA (12 nM) or yqjW (12 nM) promoter DNA as described in Materials and Methods.
|
B. subtilis SOS operator consensus sequence. Among the operator sequences listed in Table 1, there are 33 distinct sites or 66 half sites. Alignment of these distinct half site sequences gives the consensus operator sequence, 5'-CGAACATATGTTCG-3' (bottom of Table 1), which expands the previously determined consensus sequence by four internal base pairs. This sequence has also been shown to be the thermodynamically preferred sequence for LexA binding (Groban et al., submitted). Although sequences outside the operator may contribute to binding and all mismatches are not equal, there is a correlation between binding affinity and the number of mismatches relative to the consensus sequence. The SOS operators with the highest LexA affinity generally differ from the expanded consensus sequence by 2 or 3 base pairs, unless they are adjacent to another site that could contribute to binding through cooperative interactions. Lower affinity binding sites have four or five mismatches, as do the sequences that are not bound by LexA. There are no sequences in the entire genome with less than two mismatches relative to the expanded consensus sequence.
Identification of six additional LexA binding sites in the B. subtilis genome. To search for other LexA binding sites that may be located elsewhere in the B. subtilis chromosome, we searched the entire genome for the expanded consensus sequence allowing for up to five nondestabilizing mismatches. In this search we also allowed for one mismatch in both of the GAAC sections to test the possibility that we missed some binding sites in our initial search (which allowed for only one GAAC mismatch). In addition to the sequences shown in Table 1, we found 18 canonical sites. Six of these sites (listed in the middle of Table 2) are located within putative promoter regions; they were not identified in our initial search because they all contain one mismatch in both of the 5'-GAAC-3' sites. Consistent with the prediction, none of these sites displaced LexA from the recA operator when added in 50-fold molar excess over the recA site. Of the remaining 12 canonical sites located outside putative promoter regions, only 4 of these sites, located upstream of the yonT gene and inside the fliP, ylbC, and yqkA genes, bind specifically to LexA (Table 3).
|
View this table: [in a new window] |
TABLE 3. Canonical SOS boxes located outside putative promoter regions
|
Microarray analysis of DNA damage-inducible genes. According to the SOS system model, any genes with LexA binding sites overlapping their promoters in a way that inhibits RNA polymerase binding should be induced by DNA damage. We used genomic microarrays to identify genes that are induced by mitomycin C and UV radiation in a RecA+ strain, but not in a recA null strain. UV radiation and MC are known inducers of the SOS response that generate the inducing signal differently (29, 42), and our microarray analysis shows that they cause substantially different changes in overall gene expression (data not shown). By definition, an SOS gene should be induced by both treatments but only in cells containing a functional RecA protein.
Microarray experiments were done, in triplicate, on samples from RecA+ cells and recA null cells treated with either MC or UV radiation as described in Materials and Methods. During the 60-minute treatment time, there was no visible effect on growth. We found 37 genes that met the SOS induction criteria; that is, they were induced by both MC and UV in wild-type cells, but not in a recA mutant (Table 4). Thirty-three of these genes, corresponding to 18 operons, are preceded by LexA binding sites (Table 1). All but four genes (yqjY, yqjZ, uvrC, and yhaN) in Table 4 exhibited statistically significant induction, corresponding to a 99% confidence level, by both treatments. Although induction of the yqjY, yqjZ, uvrC, and yhaN genes did not meet the 99% confidence level, we presume they are transcribed from promoters containing LexA binding sites, and they show slight, but reproducible, RecA-dependent induction by both treatments. The uvrC gene has a LexA binding site overlapping its promoter, and the other three genes are apparently part of damage-inducible operons in which the upstream genes are induced. The lower induction levels for the yqjY, yqjZ, and yhaN genes are consistent with the reduced induction we observed for the downstream genes in all the putative operons containing more than two genes. For example, ynzC and yhjB induction levels by both treatments are about 20% of the yneA and yhjD levels, respectively. The induction levels of yhaN and yqjY or yqjZ by both treatments are about 50% of the yhaO and yqjW levels, respectively.
|
View this table: [in a new window] |
TABLE 4. B. subtilis genes induced by UV and MC in wild-type cells, but not in recA null mutants
|
Four of the genes listed in Table 4licA, ymaC, ymaD, and yvsGdo not have LexA binding sites. Of these, licA and the ymaCD operon have no upstream sequences resembling an SOS operator and LexA does not bind specifically to their promoter regions. The yvsG promoter contains a canonical SOS operator sequence (Table 2), but LexA does not bind specifically to DNA containing this sequence. As with the other genes listed in Table 2, no binding was detected at LexA concentrations below 400 nM and a supershift was observed at higher LexA concentrations (data not shown).
The 20 genes listed in Table 5 contain LexA binding sites in their putative promoter regions, but they did not meet our induction criteria. In every case there was no statistically significant RecA-dependent induction following one or both of the treatments. Many of them show low RecA-dependent induction but not enough to support including them in our list of SOS genes. We also found five genes that showed statistically significant RecA-dependent induction by MC, but not by UV treatment (Table 6). Of these, the yraH gene and the ycgMNO operon have canonical SOS operator sites in their promoter regions (Table 2). Both sites have seven mismatches relative to the expanded consensus, and LexA does not bind to either site. The yolC gene, which has a LexA binding site overlapping its putative promoter, is the only gene that showed RecA-dependent induction by UV, but not by MC (Table 5).
|
View this table: [in a new window] |
TABLE 5. B. subtilis genes that contain LexA binding sites within their promoter regions, but are not significantly induced by MC or UV in wild-type cells
|
|
View this table: [in a new window] |
TABLE 6. B. subtilis genes that exhibit RecA-dependent induction by MC, but not by UV
|
|
|
|---|
![]() View larger version (16K): [in a new window] |
FIG. 5. Genetic map locations of B. subtilis SOS genes. Primary (black) and secondary (gray) SOS genes are indicated, with arrows depicting the direction of transcription.
|
Although the excision repair proteins have been highly conserved in E. coli and B. subtilis, there are differences in the organization and regulation of the corresponding uvr genes. Unlike in E. coli, where the uvrA and uvrB genes are separated on the chromosome, the B. subtilis genes are contiguous and presumably transcribed from the same promoter; correspondingly, they are induced to about the same degree following DNA damageabout 8- to 12-fold, which is about twice the induction level of the E. coli genes (37). Unlike in E. coli, the B. subtilis uvrC gene has a LexA binding site and is marginally induced by MC and UV. The E. coli uvrD gene is induced five- to sevenfold by DNA damage (37) compared with about twofold for pcrA; the uvrD gene is also not part of a multigene operon like its B. subtilis counterpart. The inducible gene following pcrA in that putative operon, ligA, codes for DNA ligase, which shares 49% amino acid identity with E. coli DNA ligase; the E. coli gene has not been shown to be a damage-inducible gene.
Recombinational repair. The main proteins involved in recombinational repair are RecA and the Ruv proteins. The contiguous B. subtilis ruvA and ruvB genes apparently comprise an operon as they do in E. coli. The RuvA and RuvB proteins have also been highly conserved in the two bacteria, and genetic evidence indicates a similar enzymatic role (3). Prior to resolution of the Holliday junction in homologous recombination (and recombinational repair), the E. coli RuvA and RuvB proteins are thought to act together, presumably with RecA, to promote ATP-dependent branch migration (36). In E. coli the Holliday junction is cleaved by the RuvC protein, which is not SOS regulated (51, 52). There is no RuvC homolog in B. subtilis; however, the B. subtilis RecU protein, which is found in gram-positive, but not gram-negative, bacteria has a similar Holliday junction resolvase activity (1).
The sequences of other SOS genes suggest possible involvements in recombinational repair. The product of the yneB gene, part of the yneA operon, has significant homology with several Bacillus species site-specific recombinases of the resolvase family. The product of the ynzC gene, also part of the yneA operon, shares 31% identity with part of the B. subtilis GyrB protein. The putative parEC operon codes for the subunits of topoisomerase IV, a type II topoisomerase that is essential for segregation of replicated chromosomes in B. subtilis (19); the products of the parC and parE genes are homologous with the GyrA and GyrB subunits of DNA gyrase.
Error-prone repair (translesion DNA synthesis). The yqjW gene product shares 26% identity with E. coli UmuC protein and the N-terminal two-thirds of YgjW is 33% identical with E. coli DinP (also called DinB). The two E. coli proteins belong to the Y superfamily of DNA polymerases, which can replicate over various DNA lesions (11, 21). The better characterized is the product of the umuC gene, part of the damage-inducible umuDC operon, whose products catalyze translesion DNA synthesis during the SOS response (i.e., error-prone repair) (46). The UmuD protein is cleaved by activated RecA to produce UmuD', which associates as a dimer with UmuC to form the functional UmuCD'2 enzyme. No homolog of E. coli UmuD exists in B. subtilis.
Deletion of the B. subtilis yqjW gene decreases UV-induced mutagenesis (45), suggesting that the product of the yqjW gene and/or downstream members of the putative operon are involved in error-prone repair. yqjW is the first gene in a putative operon also containing the yqjX, yqjY, and yqjZ genes. The functions of the yqjX, yqjY, and yqjZ gene products are unknown, and none of them have any homology with E. coli UmuD. The yqjX gene product shares limited homology with the products of the B. subtilis yolD and yozL genes, which both have SOS boxes but did not meet our induction criteria. (It is noteworthy that the product of the uvrX gene, which follows yolD in a putative operon, is homologous with YqjW and other Y-family DNA polymerases.) The yqjY gene codes for a protein with a GCN5-related N-acetyltransferase (GNAT) domain that is homologous with other Bacillus acetyltransferases as well as the product of the E. coli yfiQ gene. The yqjZ gene product shares 59% identity with a conserved bacterial protein involved in polyketide biosynthesis, and it is homologous with the C-terminal domains of P. aeruginosa CTP synthase and E. coli DNA photolyase (17).
A role for DnaE in error-prone repair has been suggested by its ability to bypass certain DNA lesions and by the loss of UV-induced mutagenesis when DnaE is depleted (22). Although we did not find significant induction of the dnaE gene in our microarray analyses, it probably is an SOS gene; it has an upstream LexA binding site, and it was recently shown to be induced about threefold by both MC and nalidixic acid (22). The corresponding protein is an essential DNA polymerase that lacks 3'
5' proofreading exonuclease activity.
Filamentation. Inhibition of cell division, or filamentation, is a property that has long been associated with the SOS response. In E. coli the product of the sulA gene inhibits cell division by binding to the major component of the cell division machinery, FtsZ (6, 34). There is no sulA homolog in B. subtilis, but there is evidence that the B. subtilis yneA gene product plays a similar role in the inhibition of B. subtilis cell division. Studies of SOS-induced filamentation in yneA and lexA mutants suggest that the YneA protein suppresses cell division by inhibiting FtsZ ring formation (20). yneA is the first gene of the yneAB-ynzC operon.
Other putative DNA repair functions. Some of the SOS genes code for proteins homologous with DNA repair or modification enzymes. The yhaZ gene codes for a protein that shares 40 to 50% identity with DNA alkylation repair enzymes of a variety of bacterial species. The ydiO and ydiP gene products have some homology with each other and are similar to cytosine-specific methyltransferases from a variety of bacteria; recent evidence suggests that these proteins are responsible for the modification of BsuM restriction sites (35). The yhaO gene codes for a protein that is homologous with DNA repair endonucleases, the C-terminal end of the yhaN gene product has a DNA repair-associated ATPase domain, and yhaM has a metal-dependent phosphohydrolase domain.
Of the remaining SOS gene products, only threeAprX, CwlD, and LicAhave known activities, but their roles in the SOS response are not clear. The aprX gene codes for a subtilisin-like protease that appears to be a member of a new family of proteases (49). The gene has two putative
A promoters, but transcriptional analysis showed that only the downstream promoter is used and only during stationary phase (47). The aprX SOS box overlaps the upstream promoter sequence, which could explain why no transcription was observed from this promoter. Our induction results are consistent with repression of both promoters during exponential growth; DNA damage would induce transcription from the upstream promoter through LexA cleavage, and the downstream promoter may be induced by another mechanism during stationary phase.
The cwlD gene codes for N-acetylmuramoyl-L-alanine amidase, an enzyme that hydrolyzes a linkage in the cell wall (10). It is located immediately downstream of the ybaK gene, and together the two constitute an operon; the function of the ybaK gene product is unknown. The licA gene, the third gene of the licCBAH operon codes for an enzyme IIA component of the lichenan phosphotransferase system (47). It is not clear why licA is induced by DNA damage and not the other genes in the operon, or how the phosphotransferase system could be involved in the SOS response.
Are there more SOS genes? We have identified 33 primary SOS genes, but we do not rule out the possibility that other genes listed in Table 1 are also part of the B. subtilis SOS regulon. Indeed, our inability to detect significant dnaE induction suggests that we may have missed other inducible genes using our microarray conditions. However, assuming we have identified all of the LexA binding sites, we can put an upper limit of 54 on the number of primary SOS genes in YB886 cells. Although it is possible that our search missed some LexA binding sites, the results from our less stringent search of the entire genome argues that we probably did not. Besides the genes regulated by LexA, the induction of prophage genes has long been associated with the SOS response. Because the strains we used in this study are noninducible for prophage PBSX and have been cured of prophage SPß, we did not detect the induction of the corresponding bacteriophage genes.
A recent study of the B. subtilis oxidative stress response suggests that the SOS regulon may overlap with other regulons comprising the B. subtilis adaptational network. Macroarray analysis of B. subtilis cells exposed to hydrogen peroxide showed the induction of several SOS operons: recA, lexA, uvrBA, uvrC, dinB, dinC, yhaONM, yhaZ, yneAB-ynzC, and the secondary SOS genes ymaCD (33). However, in a similar study of B. subtilis cells exposed to the same concentration of hydrogen peroxide and analyzed using microarrays, no SOS gene induction was reported; however, there was significant down-regulation of prophage PBSX genes, contrary to what would be expected if the SOS response were activated (18). Both analyses showed significant induction of the perR and sigB regulons by hydrogen peroxide, and the macroarray analysis also showed a stringent response. Neither the perR nor sigB regulon is induced by MC or UV, and there is no stringent response following either treatment (A. I. Goranov, E. Kuester-Schoeck, J. D. Wang, and A. D. Grossman, unpublished results). Thus, there is no conclusive evidence for overlap between the oxidative stress response and the SOS response.
Differential binding and expression of B. subtilis SOS genes. The binding constants for LexA binding to B. subtilis SOS operators range from 2 nM to over 100 nM, similar to the range reported for E. coli (43). There is also a wide range of induction levels, ranging from about 2-fold to 175-fold, although the very high induction levels for dinB and dinC are probably due, in part, to low basal levels of expression. Genes that are needed during normal growth, such as recA and lexA, have higher basal levels of expression and lower induction levels. Nevertheless, there is some correlation between binding affinity and induction level. As expected, genes with higher affinity binding sites are generally induced to a greater extent, but the induction level also depends on the strength of the promoter, the position of the operator relative to the promoter, stability of the mRNA, and interactions with other molecules. Our results indicate that the position of the operator may be particularly important. For example, the operators for the dinB and dinC genes, which have high-affinity binding sites and show the greatest induction, are both located in the same position relative to the putative promoters. In general, the operators of genes that showed the greatest level of induction are located between 14 and 45 base pairs upstream of the 3' end of the 10 region of the promoter; genes with operators downstream of the 10 region were moderately induced; and genes with operators more than 50 base pairs upstream of the 10 region were not significantly induced regardless of binding affinity.
In summary, we have demonstrated a very effective and generally applicable approach for identifying specific DNA binding sites that regulate genes scattered throughout the genome. Using the results from a thermodynamic analysis of LexA binding to recA operator mutants, we identified 40 potential SOS boxes and we showed that 33 of them are bound specifically by LexA. A less stringent search did not reveal any additional sites within promoter regions, although we found six sites located outside promoters. Without the operator binding study as a guide, a search for a 14-mer allowing for up to five mismatches relative to the SOS operator consensus sequence yields over 18,000 sequences in the B. subtilis genome.
|
|
|---|
This article has been cited by other articles:
| |||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
Copyright © 2009 by the American Society for Microbiology. For an alternate route to Journals.ASM.org, visit: http://intl-journals.asm.org | More Info»