Journal of Bacteriology, September 2002, p . 4881-4890, Vol . 184, No . 17

Genome-Wide Analysis of the Stationary-Phase Sigma Factor (Sigma-H) Regulon of Bacillus subtilis

Robert A . Britton,¹ Patrick Eichenberger,² Jose Eduardo Gonzalez-Pastor,² Paul Fawcett,^2, Rita Monson,¹ Richard Losick,² and Alan D . Grossman^1*

Department of Biology, Massachusetts Institute of Technology, Cambridge, Massachusetts 02139 ,¹ Department of Molecular and Cellular Biology, Harvard University, Cambridge, Massachusetts 02138²

Received 5 April 2002/ Accepted 11 June 2002

Sigma-H, the sigH (spo0H) gene product, directs the transcription of several genes that function in the transition from exponential growth to stationary phase, including the initiation of spore formation and entry into the state of genetic competence (1, 7, 12, 20) . Sigma-H is required at an early stage of sporulation and directly activates transcription of several sporulation genes including spo0A, spo0F, kinA, spo0M, spoVG, and spoVS and the spoIIA operon (2, 4, 28, 41, 42, 46, 54, 59, 61) . Sigma-H also directs the transcription of several members of the phr family of genes, which encode secreted peptide pheromones (31, 37) . Each Phr peptide likely inhibits the activity of a corresponding Rap phosphatase that modulates entry into genetic competence, sporulation, and perhaps other processes (32, 39) .

Several of the genes that are transcribed by a sigma-H-recognized promoter have additional promoters that are recognized by other sigma factors . For example, spo0A (sporulation response regulator), ftsA (cell division), dnaG (DNA replication), sigA (encoding sigma-A, the major sigma factor), and citG (tricarboxylic acid cycle) are transcribed under sigma-H control but are also transcribed from sigma-A-dependent promoters .

In addition to genes that are under the direct control of sigma-H, there are many genes whose transcription is indirectly influenced by sigma-H . For example, during sporulation, sigma-H stimulates transcription of the master regulator of sporulation, spo0A, from a sigma-H-recognized promoter, Ps (41) . Expression from Ps is essential for efficient sporulation (49) . The product of spo0A (Spo0A), in turn, activates or represses a large number of genes, many of which are transcribed by RNA polymerase containing sigma-A . Thus, sigma-H contributes to gene expression during sporulation both directly and indirectly .

The regulation of sigma-H itself is complex . Transcription of sigH is controlled directly by the transcriptional repressor AbrB and indirectly by the phosphorylated form of Spo0A (Spo0AP), which represses abrB, and by sigma-H itself (6, 40, 52, 57) . spo0A gene expression is driven in part by sigma-H (see above), which results in increased levels of Spo0AP (under the appropriate conditions) . Increased levels of Spo0AP result in more repression of abrB and therefore increased levels of sigH transcription, thereby setting up a self-reinforcing cycle . Sigma-H activity is also controlled at the posttranscriptional, translational, and posttranslational levels and responds to a variety of external conditions including pH, carbon source, and availability of amino acids (3, 9, 10, 17, 34) . The precise mechanisms of regulation are not completely understood .

We used a combination of DNA microarray analysis and a bioinformatics approach to identify genes of the sigma-H regulon . We performed two types of DNA microarray experiments, one comparing RNAs from wild-type cells to RNAs from a sigH-null mutant and the other identifying RNAs that were induced by the overexpression of sigma-H during growth . We found that this combined microarray approach, along with the use of a hidden Markov model (HMM) database of possible sigma-H promoters, was successful both in confirming the identification of genes previously known to be directly controlled by sigma-H and in assigning many additional genes to the sigma-H regulon . Our results indicate that, in addition to its previously known role in sporulation and genetic competence, sigma-H controls many genes that are involved in helping cells adapt to conditions of nutrient depletion . The products of these genes are involved in a variety of processes, including transport, cell wall metabolism, proteolysis, and cytochrome biogenesis . In combination with the several regulatory genes that it controls, sigma-H contributes significantly to the intertwined networks that influence physiological decisions during entry into stationary phase .

To induce ectopic expression of sigH, we constructed a strain with sigH under the control of the LacI-repressible, isopropyl-ß-D-thiogalactopyranoside (IPTG)-inducible promoter Pspank-hy . pDR111 (gift from David Rudner and Federico Gueiros Filho, Harvard University) is a derivative of the Pspac-hy plasmid pJQ43 (43) that contains an additional lacO binding site to achieve better repression in the absence of the inducer IPTG . An EcoRI/BamHI fragment of pDR111 was cloned into pDG1727 (21) to generate plasmid pPE30 . A 603-bp fragment, including the ribosome binding site and the N-terminal coding sequence of sigH, was amplified by PCR with primers PE193 and PE195 (sequences available on request) and chromosomal DNA from PY79 as template . This PCR fragment was digested with SalI and SphI and cloned into pPE30 to yield pPE31 . This plasmid was integrated into PY79 by single-crossover recombination to generate strain EG232 . To analyze the transcriptional profile of a strain in which sigH was overexpressed, cells were grown in Luria-Bertani medium at 37°C to mid-exponential phase (optical density at 600 nm = 0.5), at which time the culture was split in two and IPTG was added to one culture . Samples were taken for analysis immediately after and 15, 30, and 60 min after the addition of IPTG and compared to the same time points from the parallel culture without IPTG .

DNA microarray construction. Our microarrays consist of >99% (4,074 of 4,106) of the annotated protein coding genes of the B . subtilis genome . Primers (B . subtilis ORFmers) to 4,100 of the 4,106 genes of the B . subtilis genome were purchased from Sigma-Genosys . Each gene was amplified by PCR from genomic DNA from B . subtilis strain JH642 by using Hot Start Taq Master mix (Qiagen) . A second PCR was performed with a dilution of the first PCR as a template . Of the 4,100 genes, 180 were not successfully amplified with the B . subtilis ORFmer primers . We constructed new primer sets for these 180 genes, and all but 26 were successfully amplified for spotting on the microarrays . PCR products were purified with QIAquick 96-well PCR purification kits (Qiagen) . The DNA was dried down and resuspended in 3x SSC (1x SSC is 0.15 M NaCl plus 0.015 M sodium citrate)-0.1% Sarkosyl for spotting . DNA was spotted on CMT-GAPS slides (Corning) with the Affymetrix 417 arrayer . In addition to the B . subtilis genes, four control sequences from Escherichia coli (ybaS, yfiF, yciC, and ygjU) were amplified and spotted on the array . These were chosen because there are no homologous sequences in B . subtilis . Each is represented 15 to 20 times throughout the array . Slides were prepared for hybridization (postprocessed) as previously described (13) . The quality of each set of printed arrays was checked by hybridizing to the arrays genomic DNA isolated from B . subtilis strains JH642 and PY79 that had been labeled with cyanine 3 (Cy3) or Cy5 . Briefly, genomic DNA was digested with HpaII and labeled with Cy3-dUTP or Cy5-dUTP with random primers and the Klenow fragment of DNA polymerase I . Spots that did not show proper hybridization in the control experiment were excluded from further analysis .

Sample preparation and RNA isolation. Samples of cells taken for RNA isolation were immediately mixed with an equal volume of methanol (-20°C) and kept at room temperature for 1 to 2 min . Cells were centrifuged for 5 min, and cell pellets were stored at -80°C . RNA was isolated by using RNeasy RNA isolation kits (Qiagen) or by a hot acid-phenol isolation procedure (15) . RNA prepared by either method gave similar results . RNA was treated with DNase on Qiagen columns as described by the manufacturer . The quality of the RNA was checked by visualizing the integrity of the 23S and 16S rRNA bands on an agarose gel .

Labeling and hybridization conditions. Labeled cDNA was generated from RNA samples by direct incorporation of Cy3- or Cy5-labeled dUTP into cDNA . Differentially labeled samples from two different conditions (strains) were mixed and hybridized to the DNA microarrays, and each experiment was done at least three times (see below) . For labeling, RNA (10 to 50 µg) was incubated with 1 µg of random hexamers and E . coli control RNA at 70°C for 10 min and then placed on ice for 2 min . A labeling mix containing 2x reverse transcription buffer (Life Technologies), 5 mM MgCl₂, 20 mM dithiothreitol, deoxynucleoside triphosphates (1 mM dATP, 1 mM dGTP, 1 mM dCTP, and 0.4 mM dTTP), and either Cy3-dUTP or Cy5-dUTP (Perkin-Elmer Life Sciences) was added to the RNA-primer mixture and incubated at 25°C for 5 min . Superscript II reverse transcriptase (300 U) (Life Technologies) was added, and the mixture was incubated at 25°C for 10 min and then at 42°C for 70 min . The reaction was stopped by heating the reaction mixture to 70°C for 15 min . RNA was digested by adding RNase A and RNase H and incubating the mixture at 37°C for 30 min . Unincorporated nucleotides were removed by using QiaQuick purification spin columns (Qiagen) or DyeEx spin columns (Qiagen) . Labeled cDNA was dried and resuspended in hybridization buffer (25 mM HEPES [pH 8.0], 1 mM EDTA, 0.8 µg of yeast tRNA/µl, 3x SSC, 0.2% sodium dodecyl sulfate) . Hybridizations were performed as described previously (13) . Slides were scanned on a GenePix 4000B scanner (Axon Instruments, Inc.) . E . coli control RNA corresponded to the genes ybaS, yfiF, yciC, and ygjU . These genes were amplified by PCR from E . coli with an upstream primer that contained a promoter recognized by T7 RNA polymerase . RNA was made by in vitro transcription with T7 RNA polymerase with the PCR products as template .

Data analysis . (i) Image analysis and normalization. Images were processed and analyzed with GenePix 3.0 software (Axon Instruments, Inc.) . To be considered a valid signal, 40% of the pixels in a spot had to be at least 1 standard deviation above the local background in at least one of the channels . Spots not making this cutoff were excluded from further analysis . Because many of the genes in the analysis are expressed under only one condition, we had to assign a value to the spots that did not contain a significant signal in one channel . This value was the lowest signal in a channel that met our criterion of being at least 1 standard deviation above background . Background signal was not subtracted from the signal intensity of the spots . Once spots with significant signals were identified, the two channels were normalized by making the total signal in each of the channels equal .

(ii) Determination of outliers. Genes whose expression differed significantly between the two conditions being compared were determined by two independent methods . We analyzed our data by a method similar to the previously described iterative outlier analysis (35) . Each time point was the average of at least three independent experiments (independently grown and prepared samples) . In at least one hybridization, the fluorophores were swapped to help decrease bias introduced by the dyes . In cases where multiple hybridizations were done from the same RNA sample, the data were averaged and treated as a single value for the experiment . The ratios from the independent samples were log₂ transformed, and then the data for each individual spot were averaged between the replicate experiments . We then calculated the geometric mean and standard deviation of the entire population . Any spot that had a ratio that was more than 2.5 standard deviations away from the mean was considered an outlier . Outliers were then removed from the population, and the means and standard deviations were recalculated . Once again, any spot more than 2.5 standard deviations away from the mean was considered an outlier . This process was repeated until few or no outliers were detected . In these experiments generally three iterations were needed to identify all outliers in the population .

Array data were also analyzed with the Rosetta Resolver application Axon error model (Rosetta Biosoftware) . The lists of outliers from the two analysis methods were compared, and only those genes that were considered significantly changed in both were considered further . The range of ratios of the outliers was from a high of 12 to a low of 1.6 .

HMM analysis. The sigma-H promoter sequence was modeled by using the HMMER 2.1.1 suite of software packages to create a series of closely related HMMs (http://hmmer.wustl.edu) . Known sigma-H promoter sequences from citG, ftsAZ, kinA, sigA-P3, spo0A, spo0F, spoIIA, spoVG, spoVS, ureaABC, dnaG/sigA, minC, and ytxG were manually aligned and used to hand specify models with the HMMBUILD module . The various models differed only in their accommodation of alternate direction of transcription and variation in the size of the spacer region . In all cases, a background nucleotide distribution consistent with that of B . subtilis was assigned to portions of the promoter sequences corresponding to the spacer region, while regions corresponding to the -10 and -35 boxes were assigned to an HMM "match" state . To enable detection of multiple occurrences of promoters within the genome, all HMMs were required to match globally with respect to the model but locally with respect to the B . subtilis genome . By using the HMMSEARCH module, B . subtilis genome release 14.2 was searched with each HMM, and the positions of hits were compiled and correlated with nearby reading frames .

We found a total of 433 genes that had significantly different levels of RNA from at least one time point in the sigH mutant compared to the sigH⁺ strain . Of the genes altered, 245 were dependent on sigma-H for expression (that is, they were more highly expressed in wild-type cells) and 188 had higher expression levels in the sigH mutant . Together, over 10% of the genes in the B . subtilis genome were altered in the sigH mutant (discussed further below), demonstrating the important role of sigma-H in cellular physiology . Graphs comparing the relative abundance of RNA for each gene between the two strains are shown in Fig . 1A . Points that fall on or near the line with a slope of 1 are the majority of genes whose expression is not significantly affected in the sigH-null mutant . Points that are significantly off the line (outliers) are different between the two strains . As expected, most of the effects of the sigH mutation are seen an hour into stationary phase . We found 312 genes that had a significant change at this time point, 199 that were more highly expressed in the wild type, and 113 that had increased expression in the sigH mutant . Many of the genes dependent on sigma-H had previously been characterized as dependent on Spo0A (15) .

 
Overexpression of sigH during vegetative growth. To identify many of the genes that are directly activated by sigma-H, we analyzed gene expression in response to overexpression of sigH during exponential growth in Luria-Bertani medium . Despite the complex posttranscriptional regulation of sigma-H, this strategy was quite useful . The sigH gene was placed under the control of the LacI-repressible, IPTG-inducible promoter Pspank-hy (Materials and Methods) . Cells were grown to mid-exponential phase, at which time the cultures were split in two and sigH was induced by addition of IPTG to one of the cultures . Samples of cells for RNA extraction were taken just after addition of the inducer, IPTG, and 15, 30, and 60 min after induction and were compared to the same time points in the parallel cultures without IPTG . We found 160 genes that were significantly changed in response to overexpression of sigH; 110 of these had higher levels of expression at at least one of the time points after induction .

Data from a representative induction experiment are presented in Fig . 1B . Immediately after the addition of inducer, a single outlier is detected, sigH (Fig . 1B) . This indicates how rapidly transcription from the promoter Pspank-hy is induced in response to addition of IPTG . Fifteen minutes after induction, most of the genes that are significantly changed are more highly expressed in the cells overproducing sigma-H, consistent with the function of sigma-H as an activator of transcription (Fig . 1B) . That this approach is successful is underscored by the fact that many of these genes are known to be regulated by sigma-H, such as dnaG and citG . There is a more even distribution of genes that are up and down regulated 60 min after the overexpression of sigma-H (Fig . 1B) . This is likely because we are observing some indirect effects on gene expression caused by increased levels of sigma-H and/or cells are beginning to proceed from exponential growth to stationary phase at this time . Interestingly, we observe known sigma-H promoters that additionally require the transcription factor Spo0A for activation (i.e., spoIIA operon) at 60 min after induction and not 15 min, most likely because Spo0A is more active at the transition to stationary phase, which is occurring about 60 min after addition of IPTG . Additional transcription factors that regulate gene expression at the transition from exponential growth to stationary phase such as AbrB and CodY are likely involved in the increase in genes regulated at the 60-min time point .

Identifying potential sigma-H binding sites in the genome by using an HMM approach. HMMs of the promoter sequence recognized by RNA polymerase containing sigma-H were created (Materials and Methods) by using sequences from experimentally determined sigma-H promoters reported in the literature . Promoter sequences from the genes citG, ftsAZ, kinA, sigA-P3, spo0A, spo0F, spoIIA, spoVG, spoVS, ureABC, dnaG, minC, and ytxG were used to create HMMs of the sigma-H promoter . These models were used to search the B . subtilis genome to create a database of possible sigma-H promoters . These potential promoters were then compared to the microarray data described above . Genes that showed dependence on sigma-H for expression and had an HMM-predicted promoter were considered strong candidates for direct regulation by sigma-H .

We found that the models gave us a large number of false positives . This was probably because there were relatively few promoters known to be recognized directly by RNA polymerase containing sigma-H that could be used to make the HMMs and because the consensus promoter sequence is AT rich, as is the B . subtilis genome . Thus, without the transcriptional profiling data it would have been difficult to accurately predict sigma-H promoters based on sequence alone . Because we knew that the models were not ideal, we considered the possibility that some sigma-H promoters would be missed . Therefore, we additionally used a sequence based on the consensus to search all genes that were dependent on sigma-H for expression that did not already have a predicted HMM . The two approaches were then combined to give a final list of potential sigma-H promoters (discussed below) .

Comparison of the microarray data with known sigma-H-regulated genes. To validate our experimental approach, we analyzed the behavior of known sigma-H-regulated genes in the two types of DNA microarray experiments . The experiments reliably detected genes previously known to be activated by sigma-H . We were aware of 23 known sigma-H-regulated promoters driving the expression of 33 genes (Table 1 and Fig . 2 ) . Of these 23 promoters we found sigma-H-dependent gene expression for 18 of the promoters in at least one of the microarray experiments (sigH⁺ versus sigH mutant or sigH overexpression) . Classical sigma-H-regulated genes such as citG, spoVG, and kinA showed the correct regulation in both DNA microarray experiments .

 
Five promoters known to be regulated by sigma-H were not identified in our microarray analysis, and biological or technical reasons account for the discrepancy with previous reports . Three of the genes not found in our experiments were phr genes (phrG, phrI, and phrK), which are reading frames of under 200 nucleotides . Although we successfully detected several small open reading frames on our arrays (other phr genes, spoVG, and spoVS), DNA fragments under 200 bp tended to give lower signals with higher variability . Thus, it is not surprising that some of the small genes were missed (47) . One of the other two genes (spo0A) would likely not have been detected under the experimental conditions used . The sigma-H-dependent promoter that drives expression of spo0A is a secondary promoter and has a less-than-twofold effect on the overall expression of spo0A (41), and therefore, we did not expect to find this gene in our experiments . The last gene, spo0M, is expressed at a low level during sporulation (23), and we did not detect a significant signal for this gene on our arrays .

The analysis of how the known sigma-H genes behaved in the microarray experiments was important in interpreting the data from the full-genome arrays . First, some genes that appeared to be dependent on sigma-H for expression are likely detected because the promoter for the downstream gene resides in their coding region . Examples of this are the promoters for the phr genes that reside in the coding sequence of the rap gene located immediately upstream . In these cases only a portion of the rap gene is transcribed and would result in a positive signal on the microarrays although the entire gene is not being expressed . Also, we found that promoters that require an additional transcription factor for activity (such as the requirement of Spo0A for the expression of the spoIIA operon) were detected only at time points when the transcription factor is active . Thus, in the sigH overexpression experiment, genes that were expressed at the 60-min time point, as cells were entering stationary phase, are excellent candidates for those that are regulated by additional transcription factors .

In some cases it is possible that a gene that appears to be independent of sigma-H (in the sigH⁺ versus sigH mutant experiment) is actually regulated by sigma-H . For example, ytxG is known to be under the control of sigma-H and is also directly controlled by the general stress response sigma factor sigma-B (56) . We found that sigma-B activity was increased at the T_-1 time point in the sigH mutant background, consistent with previous findings that expression of sigma-B-dependent genes is increased in sigH mutants (45, 56) . RNA from ytxG was higher in the sigH mutant due to the fact that the promoter controlled by sigma-B was more active in the absence of sigma-H (56) . We detected this regulation in the two types of microarray experiments: ytxG was induced by overexpression of sigma-H but showed a higher level of expression in a sigH mutant strain (Table 1) . This example shows that one cannot necessarily assume that a gene that is regulated by sigma-H will behave as expected in the sigH⁺ versus sigH mutant experiment and demonstrates the importance of performing both types of microarray experiments . yvyD is also known to be controlled by both sigma-H and sigma-B (11) and had a pattern of expression similar to that of ytxG (Table 1) .

Identification of genes that are strong candidates to be regulated directly by sigma-H. In addition to the previously described genes activated by sigma-H, we found 26 operons containing 54 genes that showed dependence on sigma-H for expression in the microarray experiments and had a potential sigma-H promoter within 200 bp of the start codon (except yojLM, for which the promoter was 292 bp upstream) (Table 2 and Fig . 2) . In 23 of the operons there was a promoter predicted by an HMM . Promoters for the other three operons were uncovered by using the pattern search algorithm on the SubtiList website (http://genolist.pasteur.fr/SubtiList/) . In our analysis of genes that were differentially affected in the microarrays, we inferred direct regulation by sigma-H if a good sigma-H binding site was located upstream of the gene . It seems likely that most of the other genes are indirectly regulated by sigma-H . Many of the known or putative functions of the 54 genes that may be directly controlled by sigma-H highlight the critical role of sigma-H in adaptation to nutrient deprivation . Several genes under direct sigma-H control appear to be involved in adaptation to nutrient deficiency . Many of the proteins in this class are known or predicted to be secreted and could be used to modify the extracellular environment (Table 2) . Proteins such as Vpr (extracellular serine protease) (50) could be used to scavenge for food in the extracellular environment by degrading proteins . In addition to Vpr there are other secreted proteases that are indirectly regulated by sigma-H, AprE and NprE, which also could be used to digest extracellular material . In addition there is a putative secreted nuclease, YhcR, that could be used to degrade nucleic acid that could also be used for food . (A recent study of E . coli has shown that DNA can be used as a sole source of carbon for the cell [16].) Another group of genes that could provide alternative nutrients to the cell are transporters . For example, gltP encodes a glutamate transporter and is predicted to be directly regulated by sigma-H . The yhaQ gene product has similarity to ABC transporter ATP binding proteins, likely has a role in transport, and appears to be directly regulated by sigma-H . Lastly, genes such as ccdA and the qcr operon encode proteins that are involved in the synthesis of cytochrome c and the cytochrome bc complex, respectively . These proteins are involved in the electron transport chain and may be up regulated by sigma-H in response to nutrient-limiting conditions in an attempt to generate energy . We also find the expression of resABC, which is also required for cytochrome c synthesis, to be indirectly regulated by sigma-H . Expression of resABC was previously shown to be induced upon entry into stationary phase by a putative sigma-A promoter (53) .

 
Expression of the ccdA operon was previously known to coincide with the time that sigma-H is fully active (33) . The ccdA operon transcript has been mapped by primer extension, and the authors indicated that sigma-A was likely to drive transcription of the ccdA operon . We have identified a potential sigma-H promoter that overlaps the putative sigma-A promoter . Thus, sigma-H could be responsible for additional regulation of ccdA . ccdA mutant strains are deficient in sporulation at a very late stage (48), similar to what is observed with spoVS, a gene controlled by sigma-H .

The best-known role of sigma-H is to activate sporulation . Many of the known sigma-H-controlled genes are involved in the signal transduction pathway involved in the initiation of sporulation (kinA, spo0F, and spo0A) and in the early stages of sporulation (spoIIAA, spoIIAB, and sigF) . We find that an additional histidine kinase implicated in controlling the initiation of sporulation, kinE, is controlled by sigma-H (27) . We also searched the MICADO (http://locus.jouy.inra.fr/cgi-bin/genmic/madbase/progs/madbase.operl) and JAFAN (http://bacillus.genome.ad.jp) databases to determine if any of the newly identified direct sigma-H targets display a sporulation defect when mutated . Two genes, ymaH and yoeA, are reported to be defective in endospore formation when mutated .

Recently, it was shown that natural isolates of B . subtilis can form multicellular aerial structures that may be important for the dispersal of spores and that sigH is required for the formation of these structures (5) . Two genes required for the formation of these aerial structures, yveQ and yveR, are thought to be part of a large operon responsible for exopolysaccharide production . We find that genes that are members of the yve operon that is involved in fruiting body formation, yveKLMNOPQRST-yvfABCDEF, are dependent on sigma-H .

Three additional groups of newly identified sigma-H-controlled genes are worth mentioning: transcription factors, cell wall binding proteins/autolysins, and proteins involved in detoxification . Sigma-H was already known to regulate three transcription factors, sigma-A, Spo0A, and sigma-F . Two additional putative transcription factors show sigma-H-dependent gene expression, ykoM (MarR family transcriptional regulator) and yttP (TetR/AcrR family) . How these putative transcriptional regulators contribute to gene expression will be an interesting avenue of future investigation .

The cell wall binding proteins affected by sigma-H are of interest because sigH mutants are unable to form the asymmetric septum during sporulation . With the exception of ftsA and ftsZ, we do not find that the expression of any known cell division genes is affected . We do find three genes that show dependence on sigma-H, yojL (similar to major autolysins lytE and lytF), yrvJ (similar to N-acetylmuramoyl-L-alanine amidase), and yuxL (similar to acylaminoacyl-peptidase), which are likely to be involved in the modification of the cell wall and possibly in the formation of the asymmetric septum . Interestingly two of the previously identified sigma-H-controlled genes, dacC and spoVG, may also be involved in cell wall modification (36, 38) . Alternatively, these gene products may be involved in generating nutrients for the cell by digesting cell wall material . Conversely, the major autolysins lytC, lytD, and lytF (yhdD) appear to be indirectly controlled by sigma-H and are up regulated in the sigH mutant presumably via activation of sigma-D in the sigH mutant cells (see below) .

We found four genes that are predicted to be involved in adapting to changing environmental conditions . They are yoeA (similar to multidrug efflux), yojM (similar to superoxide dismutase), ywfF (similar to efflux protein), and bsaA (putative glutathione peroxidase) . These proteins likely provide protective properties to the cell and could contribute to the effects of sigma-H on survival in stationary phase at high and low pH and in the presence of ethanol (18) .

Global changes in cell physiology. As mentioned above, we found that the expression of over 10% of the genes is significantly different between sigH⁺ and sigH mutant cells . The majority of these genes are almost certainly not under the direct control of sigma-H . Many of the genes whose expression is altered are involved in metabolism, transport of nutrients, and antibiotic production . The large numbers of genes that appear to be indirectly controlled by sigma-H are listed in supplementary material (http://mcb.harvard.edu/losick) . Many of these indirect effects could be due to increased association of alternate sigma factors, especially sigma-B and sigma-D, with core RNA polymerase in the absence of sigma-H .

Many of the changes that we observed in the sigH (spo0H) mutant were also found in spo0A mutants (15) . The regulation of three transcription factors, Spo0A, AbrB, and sigma-H, is intertwined, and mutations in one gene affect the others (see the introduction) . AbrB represses transcription of sigH, and this repression is relieved via Spo0AP repressing transcription of abrB . Thus, a spo0A mutant lacks full sigma-H activity . In turn, sigma-H is partially responsible for the activation of transcription of spo0A and expression of genes needed for (kinA and spo0F) or influencing (phr genes) phosphorylation of Spo0A . Thus, in a sigma-H mutant Spo0A activity is greatly reduced . The result is that mutations in either gene have significant effects on the activity of the other, and it is not surprising that their transcriptional profiles are similar .

An example of how a sigH mutation can have profound effects on a large number of genes via its role in activating Spo0A is highlighted by the expression of genes under the control of sigma-D . Sigma-D regulates expression of genes involved in motility, chemotaxis, and autolysin production (25) . We found that expression of many of these genes was increased in the sigH mutant . This effect is almost certainly due to the role of sigma-H in activation of Spo0A; a spo0A mutation also causes increased expression of genes dependent on sigma-D via regulation of sinI and sinR (15, 44) .

Although many of the same genes are affected in spo0A and sigH mutants, the genes directly regulated by Spo0A and sigma-H are mostly different . The only promoters that have binding sites for both Spo0A and sigma-H are involved in the initial stages of sporulation, spo0F, spoIIA (sigF), and spo0A itself .

Conclusion. Global analysis of gene expression has provided insights into the role of sigma-H in the transition from exponential growth to stationary phase and sporulation . We have identified 54 genes expressed by 26 newly identified putative sigma-H promoters . This brings the total regulon of sigma-H to potentially 49 promoters controlling the expression of at least 87 genes . Many of these are controlled by additional transcription factors, and several of these genes encode bona fide or putative transcription factors . In addition, there are many other genes that are indirectly affected by sigma-H . Identifying most of the genes on a particular pathway controlled by a specific transcription factor is a critical step in beginning to dissect how the various regulatory pathways are interconnected to form complex networks .

This work was supported in part by Public Health Service grant GM50895 (A.D.G.), by the Merck/MIT Collaborative program (A.D.G.), and by NIH grant GM18568 (R.L.) . R.A.B . was supported in part by a postdoctoral fellowship from the NIH, P.E . was a postdoctoral fellow of the Human Frontier Science Program and of the Swiss National Science Foundation, and J.E.G.-P . was supported by the Ministerio de Educacion y Ciencia Postdoctoral Fellowship (Spain) .

R.A.B., P.E., and J.E.G.-P . contributed equally to the work .

Present address: Department of Biochemistry, Stanford University, Stanford, CA 94035 .

1. Albano, M., J . Hahn, and D . Dubnau. 1987 . Expression of competence genes in Bacillus subtilis . J . Bacteriol . 169:3110-3117.
2. Antoniewski, C., B . Savelli, and P . Stragier. 1990 . The spoIIJ gene, which regulates early developmental steps in Bacillus subtilis, belongs to a class of environmentally responsive genes . J . Bacteriol . 172:86-93.
3. Asai, K., F . Kawamura, H . Yoshikawa, and H . Takahashi. 1995 . Expression of kinA and accumulation of ^H at the onset of sporulation in Bacillus subtilis . J . Bacteriol . 177:6679-6683.
4. Bai, U., M . Lewandoski, E . Dubnau, and I . Smith. 1990 . Temporal regulation of the Bacillus subtilis early sporulation gene spo0F . J . Bacteriol . 172:5432-5439.
5. Branda, S . S., J . E . Gonzalez-Pastor, S . Ben-Yehuda, R . Losick, and R . Kolter. 2001 . Fruiting body formation by Bacillus subtilis . Proc . Natl . Acad . Sci . USA 98:11621-11626.
6. Burkholder, W . F., and A . D . Grossman. 2000 . Regulation of the initiation of endospore formation in Bacillus subtilis, p . 151-166 . In Y . V . Brun and L . J . Shimkets (ed.), Prokaryotic development . ASM Press, Washington, D.C.
7. Carter, H . L., III, and C . P . Moran, Jr. 1986 . New RNA polymerase sigma factor under spo0 control in Bacillus subtilis . Proc . Natl . Acad . Sci . USA 83:9438-9442.
8. Carter, H . L., III, L . F . Wang, R . H . Doi, and C . P . Moran, Jr. 1988 . rpoD operon promoter used by ^H-RNA polymerase in Bacillus subtilis . J . Bacteriol . 170:1617-1621.
9. Cosby, W . M., and P . Zuber. 1997 . Regulation of Bacillus subtilis ^H (Spo0H) and AbrB in response to changes in external pH . J . Bacteriol . 179:6778-6787.
10. Dixon, L . G., S . Seredick, M . Richer, and G . B . Spiegelman. 2001 . Developmental gene expression in Bacillus subtilis crsA47 mutants reveals glucose-activated control of the gene for the minor sigma factor ^H . J . Bacteriol . 183:4814-4822.
11. Drzewiecki, K., C . Eymann, G . Mittenhuber, and M . Hecker. 1998 . The yvyD gene of Bacillus subtilis is under dual control of ^B and ^H . J . Bacteriol . 180:6674-6680.
12. Dubnau, E., J . Weir, G . Nair, L . Carter III, C . Moran, Jr., and I . Smith. 1988 . Bacillus sporulation gene spo0H codes for ³⁰ (^H) . J . Bacteriol . 170:1054-1062.
13. Eisen, M . B., and P . O . Brown. 1999 . DNA arrays for analysis of gene expression . Methods Enzymol . 303:179-205.
14. Errington, J. 1996 . Determination of cell fate in Bacillus subtilis . Trends Genet . 12:31-34.
15. Fawcett, P., P . Eichenberger, R . Losick, and P . Youngman. 2000 . The transcriptional profile of early to middle sporulation in Bacillus subtilis . Proc . Natl . Acad . Sci . USA 97:8063-8068.
16. Finkel, S . E., and R . Kolter. 2001 . DNA as a nutrient: novel role for bacterial competence gene homologs . J . Bacteriol . 183:6288-6293.
17. Frisby, D., and P . Zuber. 1994 . Mutations in pts cause catabolite-resistant sporulation and altered regulation of spo0H in Bacillus subtilis . J . Bacteriol . 176:2587-2595.
18. Gaidenko, T . A., and C . W . Price. 1998 . General stress transcription factor ^B and sporulation transcription factor ^H each contribute to survival of Bacillus subtilis under extreme growth conditions . J . Bacteriol . 180:3730-3733.
19. Gonzy-Treboul, G., C . Karmazyn-Campelli, and P . Stragier. 1992 . Developmental regulation of transcription of the Bacillus subtilis ftsAZ operon . J . Mol . Biol . 224:967-979.
20. Grossman, A . D. 1995 . Genetic networks controlling the initiation of sporulation and the development of genetic competence in Bacillus subtilis . Annu . Rev . Genet . 29:477-508.
21. Guerout-Fleury, A . M., K . Shazand, N . Frandsen, and P . Stragier. 1995 . Antibiotic-resistance cassettes for Bacillus subtilis . Gene 167:335-336.
22. Haldenwang, W . G. 1995 . The sigma factors of Bacillus subtilis . Microbiol . Rev . 59:1-30.
23. Han, W . D., S . Kawamoto, Y . Hosoya, M . Fujita, Y . Sadaie, K . Suzuki, Y . Ohashi, F . Kawamura, and K . Ochi. 1998 . A novel sporulation-control gene (spo0M) of Bacillus subtilis with a ^H-regulated promoter . Gene 217:31-40.
24. Healy, J., J . Weir, I . Smith, and R . Losick. 1991 . Post-transcriptional control of a sporulation regulatory gene encoding transcription factor ^H in Bacillus subtilis . Mol . Microbiol . 5:477-487.
25. Helmann, J . D., and C . P . Moran, Jr. 2001 . RNA polymerase and sigma factors, p . 289-312 . In A . L . Sonenshein, J . A . Hoch, and R . Losick (ed.), Bacillus subtilis and its closest relatives: from genes to cells . ASM Press, Washington, D.C.
26. Ishikawa, S., Y . Hara, R . Ohnishi, and J . Sekiguchi. 1998 . Regulation of a new cell wall hydrolase gene, cwlF, which affects cell separation in Bacillus subtilis . J . Bacteriol . 180:2549-2555.
27. Jiang, M., W . Shao, M . Perego, and J . A . Hoch. 2000 . Multiple histidine kinases regulate entry into stationary phase and sporulation in Bacillus subtilis . Mol . Microbiol . 38:535-542.
28. Johnson, W . C., C . P . Moran, Jr., and R . Losick. 1983 . Two RNA polymerase sigma factors from Bacillus subtilis discriminate between overlapping promoters for a developmentally regulated gene . Nature 302:800-804.
29. Kiel, J . A., J . M . Boels, G . Beldman, and G . Venema. 1991 . Molecular cloning and nucleotide sequence of the glycogen branching enzyme gene (glgB) from Bacillus stearothermophilus and expression in Escherichia coli and Bacillus subtilis . Mol . Gen . Genet . 230:136-144.
30. Kunst, F., N . Ogasawara, I . Moszer, A . M . Albertini, G . Alloni, V . Azevedo, M . G . Bertero, P . Bessieres, A . Bolotin, S . Borchert, R . Borriss, L . Boursier, A . Brans, M . Braun, S . C . Brignell, S . Bron, S . Brouillet, C . V . Bruschi, B . Caldwell, V . Capuano, N . M . Carter, S . K . Choi, J . J . Codani, I . F . Connerton, A . Danchin, et al. 1997 . The complete genome sequence of the gram-positive bacterium Bacillus subtilis . Nature 390:249-256.
31. Lazazzera, B . A., I . G . Kurtser, R . S . McQuade, and A . D . Grossman. 1999 . An autoregulatory circuit affecting peptide signaling in Bacillus subtilis . J . Bacteriol . 181:5193-5200.
32. Lazazzera, B . A., J . M . Solomon, and A . D . Grossman. 1997 . An exported peptide functions intracellularly to contribute to cell density signaling in B . subtilis . Cell 89:917-925.
33. Le Brun, N . E., J . Bengtsson, and L . Hederstedt. 2000 . Genes required for cytochrome c synthesis in Bacillus subtilis . Mol . Microbiol . 36:638-650.
34. Liu, J., W . M . Cosby, and P . Zuber. 1999 . Role of Lon and ClpX in the post-translational regulation of a sigma subunit of RNA polymerase required for cellular differentiation in Bacillus subtilis . Mol . Microbiol . 33:415-428.
35. Loos, A., C . Glanemann, L . B . Willis, X . M . O'Brien, P . A . Lessard, R . Gerstmeir, S . Guillouet, and A . J . Sinskey. 2001 . Development and validation of Corynebacterium DNA microarrays . Appl . Environ . Microbiol . 67:2310-2318.
36. Matsuno, K., and A . L . Sonenshein. 1999 . Role of SpoVG in asymmetric septation in Bacillus subtilis . J . Bacteriol . 181:3392-3401.
37. McQuade, R . S., N . Comella, and A . D . Grossman. 2001 . Control of a family of phosphatase regulatory genes (phr) by the alternate sigma factor sigma-H of Bacillus subtilis . J . Bacteriol . 183:4905-4909.
38. Pedersen, L . B., T . Murray, D . L . Popham, and P . Setlow. 1998 . Characterization of dacC, which encodes a new low-molecular-weight penicillin-binding protein in Bacillus subtilis . J . Bacteriol . 180:4967-4973.
39. Perego, M., and J . A . Brannigan. 2001 . Pentapeptide regulation of aspartyl-phosphate phosphatases . Peptides 22:1541-1547.
40. Perego, M., G . B . Spiegelman, and J . A . Hoch. 1988 . Structure of the gene for the transition state regulator, abrB: regulator synthesis is controlled by the spo0A sporulation gene in Bacillus subtilis . Mol . Microbiol . 2:689-699.
41. Predich, M., G . Nair, and I . Smith. 1992 . Bacillus subtilis early sporulation genes kinA, spo0F, and spo0A are transcribed by the RNA polymerase containing ^H . J . Bacteriol . 174:2771-2778.
42. Price, V . A., I . M . Feavers, and A . Moir. 1989 . Role of ^H in expression of the fumarase gene (citG) in vegetative cells of Bacillus subtilis 168 . J . Bacteriol . 171:5933-5939.
43. Quisel, J . D., W . F . Burkholder, and A . D . Grossman. 2001 . In vivo effects of sporulation kinases on mutant Spo0A proteins in Bacillus subtilis . J . Bacteriol . 183:6573-6578.
44. Rashid, M . H., and J . Sekiguchi. 1996 . flaD (sinR) mutations affect SigD-dependent functions at multiple points in Bacillus subtilis . J . Bacteriol . 178:6640-6643.
45. Ray, C., M . Igo, W . Shafer, R . Losick, and C . P . Moran, Jr. 1988 . Suppression of ctc promoter mutations in Bacillus subtilis . J . Bacteriol . 170:900-907.
46. Resnekov, O., A . Driks, and R . Losick. 1995 . Identification and characterization of sporulation gene spoVS from Bacillus subtilis . J . Bacteriol . 177:5628-5635.
47. Richmond, C . S., J . D . Glasner, R . Mau, H . Jin, and F . R . Blattner. 1999 . Genome-wide expression profiling in Escherichia coli K-12 . Nucleic Acids Res . 27:3821-3835.
48. Schiott, T., and L . Hederstedt. 2000 . Efficient spore synthesis in Bacillus subtilis depends on the CcdA protein . J . Bacteriol . 182:2845-2854.
49. Siranosian, K . J., and A . D . Grossman. 1994 . Activation of spo0A transcription by ^H is necessary for sporulation but not for competence in Bacillus subtilis . J . Bacteriol . 176:3812-3815.
50. Sloma, A., G . A . Rufo, Jr., K . A . Theriault, M . Dwyer, S . W . Wilson, and J . Pero. 1991 . Cloning and characterization of the gene for an additional extracellular serine protease of Bacillus subtilis . J . Bacteriol . 173:6889-6895.
51. Stragier, P., and R . Losick. 1996 . Molecular genetics of sporulation in Bacillus subtilis . Annu . Rev . Genet . 30:297-341.
52. Strauch, M . A. 1995 . Delineation of AbrB-binding sites on the Bacillus subtilis spo0H, kinB, ftsAZ, and pbpE promoters and use of a derived homology to identify a previously unsuspected binding site in the bsuB1 methylase promoter . J . Bacteriol . 177:6999-7002.
53. Sun, G., E . Sharkova, R . Chesnut, S . Birkey, M . F . Duggan, A . Sorokin, P . Pujic, S . D . Ehrlich, and F . M . Hulett. 1996 . Regulators of aerobic and anaerobic respiration in Bacillus subtilis . J . Bacteriol . 178:1374-1385.
54. Tatti, K . M., H . L . Carter III, A . Moir, and C . P . Moran, Jr. 1989 . Sigma H-directed transcription of citG in Bacillus subtilis . J . Bacteriol . 171:5928-5932.
55. van Dijl, J . M., A . Bolhuis, H . Tjalsma, J . D . H . Jongbloed, A . De Jong, and S . Bron. 2001 . Protein transport pathways in Bacillus subtilis: a genome-based road map, p . 337-355 . In A . L . Sonenshein, J . A . Hoch, and R . Losick (ed.), Bacillus subtilis and its closest relatives: from genes to cells . ASM Press, Washington, D.C.
56. Varon, D., M . S . Brody, and C . W . Price. 1996 . Bacillus subtilis operon under the dual control of the general stress transcription factor ^B and the sporulation transcription factor ^H . Mol . Microbiol . 20:339-350.
57. Weir, J., M . Predich, E . Dubnau, G . Nair, and I . Smith. 1991 . Regulation of spo0H, a gene coding for the Bacillus subtilis ^H factor . J . Bacteriol . 173:521-529.
58. Wray, L . V., Jr., A . E . Ferson, and S . H . Fisher. 1997 . Expression of the Bacillus subtilis ureABC operon is controlled by multiple regulatory factors including CodY, GlnR, TnrA, and Spo0H . J . Bacteriol . 179:5494-5501.
59. Wu, J . J., P . J . Piggot, K . M . Tatti, and C . P . Moran, Jr. 1991 . Transcription of the Bacillus subtilis spoIIA locus . Gene 101:113-116.
60. Youngman, P., J . B . Perkins, and R . Losick. 1984 . Construction of a cloning site near one end of Tn917 into which foreign DNA may be inserted without affecting transposition in Bacillus subtilis or expression of the transposon-borne erm gene . Plasmid 12:1-9.
61. Zuber, P., and R . Losick. 1987 . Role of AbrB in Spo0A- and Spo0B-dependent utilization of a sporulation promoter in Bacillus subtilis . J . Bacteriol . 169:2223-2230.

Free Online Full-text Article

What Is Bioengineering?, What Is Fermentation?, What Is Genetics?, What Is Staphylococcus Aureus?, What Is MIC?, s, Microbes, a, Microorganism, i, Microbe, c, Microbiology, r, Bacterium, n, Streptococcal, i, Biofilms, r, Bacteria, i, Bacillus subtilis, r, Enterobacters, i, Escherichia coli, r, Microbial, c, Bacillus subtilis, i, Haemophilus, a, Bacteria, e, Escherichia coli, e, Haemophilus, r, Phage, n, Escherichia coli, e, S. cerevisiae, e, Functional genomics, i, Biofilms, o, Antibiotics, e, Lactobacilli, e, Anaerobe, i, Staphylococcus aureus




 

© 2005 Transgalactic Ltd (manufacturer of Bioscreen C software) | Privacy Statement | P.O. Box 1393, 00101 Helsinki, Finland, phone: +358 9 85172920, fax: +358 9 8749481, e-mail: microbiology@bionewsonline.com
 

Last modified: May 25, 2005