Genes and Development Attend a BioResearch Product Faire

Home Help [Feedback] [For Subscribers] [Archive] [Search] [Contents]
 QUICK SEARCH:   [advanced]


     


GENES & DEVELOPMENT 20:1776-1789, 2006
©2006 by Cold Spring Harbor Laboratory Press; ISSN 0890-9369/ $5.00
This Article
Right arrow Abstract Freely available
Right arrow Full Text (PDF)
Right arrow Supplemental Research Data
Right arrow Alert me when this article is cited
Right arrow Alert me if a correction is posted
Right arrow Citation Map
Services
Right arrow Email this article to a friend
Right arrow Similar articles in this journal
Right arrow Similar articles in PubMed
Right arrow Alert me to new issues of the journal
Right arrow Download to citation manager
Right arrow reprints & permissions
Citing Articles
Right arrow Citing Articles via HighWire
Right arrow Citing Articles via Google Scholar
Google Scholar
Right arrow Articles by Nonaka, G.
Right arrow Articles by Rhodius, V. A.
Right arrow Search for Related Content
PubMed
Right arrow PubMed Citation
Right arrow Articles by Nonaka, G.
Right arrow Articles by Rhodius, V. A.
Social Bookmarking
 Add to CiteULike   Add to Connotea   Add to Del.icio.us   Add to Digg   Add to Reddit   Add to Technorati  
What's this?

Regulon and promoter analysis of the E. coli heat-shock factor, {sigma}32, reveals a multifaceted cellular response to heat stress

Gen Nonaka1,3, Matthew Blankschien4, Christophe Herman1,4, Carol A. Gross1,2,5 and Virgil A. Rhodius1

1 Department of Microbiology and Immunology, University of California at San Francisco, San Francisco, California 94143, USA; 2 Department of Cell and Tissue Biology, University of California at San Francisco, San Francisco, California 94143, USA


    Abstract
 Top
 Abstract
 Results
 Discussion
 Materials and methods
 Acknowledgments
 References
 
The heat-shock response (HSR), a universal cellular response to heat, is crucial for cellular adaptation. In Escherichia coli, the HSR is mediated by the alternative {sigma} factor, {sigma}32. To determine its role, we used genome-wide expression analysis and promoter validation to identify genes directly regulated by {sigma}32 and screened ORF overexpression libraries to identify {sigma}32 inducers. We triple the number of genes validated to be transcribed by {sigma}32 and provide new insights into the cellular role of this response. Our work indicates that the response is propagated as the regulon encodes numerous global transcriptional regulators, reveals that {sigma}70 holoenzyme initiates from 12% of {sigma}32 promoters, which has important implications for global transcriptional wiring, and identifies a new role for the response in protein homeostasis, that of protecting complex proteins. Finally, this study suggests that the response protects the cell membrane and responds to its status: Fully 25% of {sigma}32 regulon members reside in the membrane and alter its functionality; moreover, a disproportionate fraction of overexpressed proteins that induce the response are membrane localized. The intimate connection of the response to the membrane rationalizes why a major regulator of the response resides in that cellular compartment.

[Keywords: Heat-shock response; {sigma}32; transcription; microarray]

Received March 13, 2006; revised version accepted April 25, 2006.


When cells are shifted from low to high temperature, synthesis of the heat-shock proteins (hsps) is rapidly and selectively induced. The heat-shock response (HSR), was first identified by Ritossa (1963)Go, who showed that exposure to heat lead to transient changes in the puffing pattern of salivary chromosomes in Drosophila; Tissieres et al. (1974)Go demonstrated that these changes reflected the transient induction of several proteins. Initially, hsp function was unclear; however, experiments in several organisms revealed that many hsps were chaperones that promote protein folding (Pelham 1986Go; Beckmann et al. 1990Go; Gaitanaris et al. 1990Go; Skowyra et al. 1990Go). These studies not only suggested that a major function of the HSR is to maintain the protein folding state of the cell, but also indicated that some of these chaperones, such as Hsp70 and Hsp90, are present in all organisms (Bardwell and Craig 1984Go, 1987Go). Thus, both the HSR and some hsps are universally conserved among organisms.

In Escherichia coli, {sigma}32, an alternative {sigma} factor, controls the HSR by directing RNA polymerase to transcribe hsps (Yamamori and Yura 1980Go; Grossman et al. 1984Go; Taylor et al. 1984Go; Cowing et al. 1985Go). Synthesis of hsps is induced upon temperature upshift and repressed upon temperature downshift (Lemaux et al. 1978Go; Yamamori et al. 1978Go; Straus et al. 1987Go, 1989Go; Taura et al. 1989Go), thereby allowing a rapid cellular response to changes in temperature. {sigma}32 is controlled by negative feedback loops controlling its activity (Straus et al. 1989Go; Blaszczak et al. 1999Go) and stability (Straus et al. 1987Go) and a feed forward loop controlling its synthesis (Yuzawa et al. 1993Go; Morita et al. 1999aGo, bGo; El-Samad et al. 2005Go). {sigma}32 regulon members mediate both feedback loops: The FtsH protease controls {sigma}32 stability (Herman et al. 1995Go; Tatsuta et al. 1998Go), and the DnaKJ and GroEL/S chaperone machines control {sigma}32 activity (Gamer et al. 1996Go; Tatsuta et al. 1998Go; Guisbert et al. 2004Go), ensuring homeostatic control.

Given the importance of the heat-shock response for cellular adaptation, it is crucial to identify the genes that are induced so that these cellular processes can be identified and investigated. Several global approaches have been utilized in E. coli to identify genes under {sigma}32 control. Both two-dimensional protein gels (Lemaux et al. 1978Go) and global transcriptional approaches have identified genes induced after transfer to high temperature (Chuang et al. 1993bGo; Richmond et al. 1999Go) or after overexpression of {sigma}32 (Zhao et al. 2005Go); however, there has been no systematic determination of whether these induced genes are directly expressed from {sigma}32-dependent promoters. The few cases where the functions of {sigma}32-dependent genes were identified indicates the importance of this determination: Hsp33, a widely conserved hsp, is a redox activated chaperone (Jakob et al. 1999Go), thereby providing the first identification of a chaperone of this type, and FtsJ was shown to be a methyltransferase whose substrate is 23S RNA, revealing an unexpected link between the HSR and RNA (Bugl et al. 2000Go).

We report the identification of most {sigma}32-dependent genes obtained by combining whole-genome expression analysis of genes induced after {sigma}32 overexpression with start site mapping of these genes and in vitro transcription studies, and also on the types of overproduced proteins that signal induction of {sigma}32. Our analysis has almost tripled the number of genes validated to be expressed from {sigma}32-dependent promoters, and revealed that the HSR response targets multiple processes, protects complex proteins, and targets the membrane as well as the cytoplasm. That the {sigma}32-mediated response is intimately related to membrane functionality rationalizes why the FtsH protease controlling {sigma}32 stability resides in the cytoplasmic membrane.


    Results
 Top
 Abstract
 Results
 Discussion
 Materials and methods
 Acknowledgments
 References
 

Identification of genes whose transcription increases following overexpression of {sigma}32

To identify genes regulated by {sigma}32, we compared expression from cells having a plasmid-borne IPTG-inducible rpoH gene with cells carrying the empty vector using whole-genome expression analysis. This approach is preferable to examining cells lacking {sigma}32 because such cells are very sick and grow slowly. Because accumulation of DnaKJ and GroELS damps expression due to feedback inhibition (Gamer et al. 1996Go; Guisbert et al. 2004Go), we first determined the kinetics of hsp synthesis following {sigma}32 overexpression. Hsp synthesis peaked at ~10–15 min following {sigma}32 overexpression (Fig. 1A); thus, we analyzed an expression time course ranging from 0 to 20 min. Relative mRNA levels were determined by parallel two-color hybridization to glass-slide cDNA microarrays. SAM analysis (statistical analysis of microarrays) (Tusher et al. 2001Go) of gene expression from four independent cultures at 10 min after induction indicated that 105 genes were significantly induced. We used hierarchical clustering to examine the RNA expression pattern of these genes throughout the time course of induction (Fig. 1B) and then determined whether the induced genes consisted of subgroups with distinct induction characteristics using SOM (self-organized maps). We find that induced genes show both temporal distinction (rapid and slower responders) (Fig. 1C) and quantitative differences (strong and weak responders) (Fig. 1D).


Figure 1
View larger version (37K):
[in this window]
[in a new window]
 
Figure 1. Expression profiles of {sigma}32-regulon members after overexpressing rpoH. (A) Activity of {sigma}32 following overexpression of rpoH. An exponential phase culture of strain CAG50002 (which carries an IPTG-inducible copy of rpoH) growing at 30°C in M9 complete–methionine was induced with IPTG at OD450 = 0.3 (t = 0). At various times, pulse chase analysis was used to determine the rate of synthesis of two {sigma}32-dependent hsps, DnaK, and GroEL. Data is normalized to their synthesis rates at induction (0 min). (B) Hierarchical clustering of 105 genes whose expression is significantly altered following rpoH overexpression (CAG50002) vs. wild type (CAG50001). The color chart illustrates the average expression level at each time point for each gene from three time course experiments. Red denotes increased and green denotes decreased mRNA expression in CAG50002 vs. CAG50001: Maximum intensity represents greater than fourfold change. Time in minutes after induction of rpoH in the time-course experiments is indicated at the top of the figure; genes are identified by their unique ID and name. (C,D) SOM analysis of significantly induced genes following rpoH overexpression. The expression ratios for each gene across three time courses were averaged for each time point. Genes were partitioned based on their induction kinetics (fast/slow; C) or magnitude of induction (strong/weak; D). Each line represents an average trace of expression pattern for that group of genes. Relative Expression Level indicates mean and variance normalized (C) or raw (D) log2 (rpoH overexpressed/wild type) expression ratios.

 

Identification of induced genes having a {sigma}32 promoter

First, we organized the significantly induced genes into transcription units (TUs) and, where appropriate, expanded the TUs to include other gene members that were less strongly induced based on the criterion that they are adjacent, in the same orientation, induced with similar kinetics, and have an expression pattern typical of an operon (the first gene has the highest induction, followed by decreasing induction of downstream genes). This gave 127 genes organized in 66 potential TUs. We then searched the 300-nucleotide (nt) region upstream of each TU for candidate {sigma}32-promoter sequences using MEME and BioProspector: This predicted 42 promoters upstream of 40/66 TUs (Table 1A,B), of which 14 were previously identified {sigma}32 promoters and 28 were new predictions.


View this table:
[in this window]
[in a new window]
 
Table 1. {sigma}32-regulon members

 
We tested whether the newly predicted promoters and TUs lacking known or predicted promoters had {sigma}32-dependent start sites by using 5'RACE (rapid amplification of cDNA ends) to compare the 5' ends of mRNAs from cells overexpressing {sigma}32 with those from cells lacking {sigma}32 (Fig. 2). This confirmed 18 out of 28 promoter predictions, using as a criterion a {sigma}32-dependent start site located just downstream of a promoter prediction and identified four new {sigma}32-dependent start sites upstream of the 26 TUs with no predicted promoters (Table 1A). In 11 instances, either no {sigma}32-dependent start site was identified near a prediction or the start site location was inconsistent with that of the prediction. Given that the rrnBP1 has an overlapping {sigma}32 and {sigma}70-dependent promoter (Newlands et al. 1993Go) and dual recognition promoters would score as {sigma}32 independent in the 5'RACE assay, we retested these 11 promoters in vitro for transcription by both {sigma}32 and {sigma}70 holoenzyme. This identified five dual {sigma}32/{sigma}70-dependent promoters and two additional {sigma}32-dependent promoters (Fig. 3). In toto, we identified 29 new {sigma}32 promoters (Table 1A). We also tested seven previously proposed {sigma}32 promoters located upstream of genes that were not significantly induced on our arrays and detected no {sigma}32-dependent transcripts (Table 1C). Finally, neither 5'RACE nor in vitro transcription identified functional promoters for three of the {sigma}32 promoters proposed by Zhao et al. (2005)Go (ldhA, macB, and ybbN; Table 1B) based on electrophoretic mobility gel shift by {sigma}32-holoenzyme.


Figure 2
View larger version (44K):
[in this window]
[in a new window]
 
Figure 2. 5'RACE identified 22 {sigma}32-dependent transcription starts. To identify {sigma}32-specific transcripts, mRNA from rpoH+ and rpoH strains (CAG50002 and CAG50003, respectively) was 5'-labeled with an RNA oligo, reverse-transcribed, amplified by PCR, and then visualized by 7.5% PAGE (see Materials and Methods). Bands present in rpoH+ but not in rpoH reactions were regarded as {sigma}32 specific. Two known {sigma}32 promoters, dnaK and ibpA, were tested; 22 new {sigma}32-dependent promoters were identified from the 54 newly identified {sigma}32-induced TUs. Displayed are duplicate examples of two new (ibpA, ileS) and two known (dnaK, ibpA) {sigma}32-specific promoters.

 


Figure 3
View larger version (29K):
[in this window]
[in a new window]
 
Figure 3. In vitro transcription assays identified seven new promoters transcribed by {sigma}32. Multiround in vitro transcription assays were performed to test 11 predicted and six previously documented promoters. Each promoter template was tested with RNAP containing either {sigma}70 or {sigma}32: {sigma}32-dependent transcripts were obtained from seven promoters, five of which also generated {sigma}70-dependent transcripts (yadF, ybeD, glnS, yceJ, and xerD).

 

Characteristics of {sigma}32 promoters

The 29 newly identified {sigma}32 promoters, 20 confirmed promoters, and two previously validated promoters upstream of sequences not present on our arrays (Table 1C) comprise a total of 51 promoters that drive the expression of 49 TUs. The sequence logos of the conserved sequence motifs upstream of the 50 chromosomal {sigma}32 promoters (note repE is on the F factor) together with their information content are displayed in Figure 4A. Virtually all of the total information content of these promoters (18.3 bits) is from the conserved –10 (CCCCATWT) and –35 (TTGAAA) core motifs, with very little contribution from flanking sequences. Figure 4, B and C display histograms of the distance distributions of the promoter elements from each other. Most promoters prefer a 7-nt spacing between –10 and +1 and a 13–14-nt spacing between –10 and –35. However, five promoters contained very short spacers (Fig. 4B). Four of these have very poor –35 sequences, raising the possibility that at these promoters, {sigma}32 makes nonspecific interactions with sequences further upstream at the optimal spacing. Interestingly, four of these promoters contain A/T-rich sequences in the spacer region, suggesting that the intervening DNA may distort upon binding of RNAP. Finally, 66% of all {sigma}32 promoters are located within 100 nt of the downstream gene translation start (Fig. 4D), as is the case for known {sigma}70 promoters (Burden et al. 2005Go).


Figure 4
View larger version (23K):
[in this window]
[in a new window]
 
Figure 4. Core promoter motifs of 50 {sigma}32 promoters. Fifty of the 51 validated {sigma}32-dependent promoters were used to derive sequence logos of the core motifs (see Table 1A,B) (the repE promoter was excluded due to the large number of multiple start points that made it difficult to confidently identify the upstream –10 and –35 motifs). The sequences were initially aligned by their start sites and WCONSENSUS was used to search small windows upstream for the conserved –10 and –35 motifs (see Materials and Methods). (A) Sequence logos (http://weblogo.berkeley.edu; Crooks et al. 2004Go) of the –35, –10, and +1 start site motifs. Note that only 43 promoters were used to derive the start site motif. Seven promoters were excluded that had either multiple starts with no clearly preferred position (b0473 P1 and P2, b1057, b3179, and b3400) or no defined starts (b0209 and b2954). The information content (Iseq) of each conserved sequence window is indicated. (B) Frequency distribution of distance between –35 (TTGAAA) and –10 (CCCATAT) motifs. (C) Frequency distribution of distance between –10 (CCCATAT) motif and +1 (A/G) start. (D) Histogram of the 50 {sigma}32-promoter start sites upstream of the gene translation start.

 
There are several A/T-rich motifs directly upstream of the –35 core motif reminiscent of the A/T-rich UP-elements that bind the {alpha} C-terminal domains ({alpha}CTDs) of RNAP (Ross et al. 1993Go; Estrem et al. 1998Go). We termed these motifs "Complete," "Proximal," and "Distal" by analogy to their UP element counterparts (Fig. 5). About 40% of all {sigma}32 promoters contain such motifs, with about half having a distal motif only and the remainder having a complete motif or a combination of complete, distal, and proximal motifs (Fig. 6). We suggest that these A/T-rich {sigma}32-promoter motifs are likely to be {alpha}CTD-binding sequences; in fact, the upstream sequences of the rrnB P1 promoter, which contain all three motifs (Fig. 6), has already been shown to stimulate {sigma}32-dependent transcription (Newlands et al. 1993Go). Interestingly, neither {sigma}32 (Fig. 6) nor {sigma}S promoters (Typas and Hengge 2005Go) have solely proximal sites, although such sites are common in {sigma}70 promoters. This is probably because the {sigma}70 residues that contact {alpha}CTD and mediate enhanced transcription from the proximal site are not conserved in {sigma}S (Typas and Hengge 2005Go) and are poorly conserved in {sigma}32 (Supplementary Fig. 1).


Figure 5
View larger version (27K):
[in this window]
[in a new window]
 
Figure 5. Sequence logos of A/T-rich motifs upstream of the {sigma}32 promoters. WCONSENSUS was used to search sequences upstream of the –35 element of 50 {sigma}32 promoters (excluding repE). Distal, Proximal, and Complete motifs were identified from the search windows –46 to –60, –36 to –51, and –36 to –60, respectively, by aligning the promoters with respect to their –35 elements and assigning the first "T" of the –35 motif as position –35.

 


Figure 6
View larger version (21K):
[in this window]
[in a new window]
 
Figure 6. Location of A/T-rich motifs in {sigma}32 promoters. The location and arrangement of the Distal, Proximal, and Complete A/T-rich motifs identified in Figure 5 are shown for the 50 {sigma}32 promoters.

 

Functional analysis of the {sigma}32regulon members

We examined the location and functional classification of validated members of the {sigma}32 regulon (Table 2A,B). A significant fraction (~25%) of known regulon members reside in the inner membrane, raising the possibility that one role of the response may be to monitor and preserve the membrane during stress. In addition to chaperones and proteases, the regulon is involved in maintaining the integrity of cellular DNA; the connection of the {sigma}32 regulon to the RNA state is more extensive than previously realized; and the regulon is intimately connected to central metabolism and to transport of substrates. Finally, {sigma}32 is a master regulator, altering expression of many transcription factors as well as activity of the transcriptional apparatus itself. Thus, the {sigma}32-mediated hsr has a general role in maintaining cellular homeostasis upon exposure to heat stress.


View this table:
[in this window]
[in a new window]
 
Table 2. Location and functional classification of {sigma}32-regulon members

 

Identification of endogenous proteins that induce the {sigma}32 HSR

A global analysis of proteins that induce a response when overexpressed is a powerful method for determining the signals normally sensed by the response. For example, overexpressed porins induce the {sigma}E envelope stress response (Mecsas et al. 1993Go); later work demonstrated that unassembled porins initiate the signal transduction pathway that induces {sigma}E (Walsh et al. 2003Go). We used the ASKA library, which contains every ORF of E. coli K-12 expressed from an IPTG-inducible promoter on separate plasmids (see Table 4, below) to search for inducers of the {sigma}32 HSR. The plasmids were pooled, introduced into a reporter strain that expresses the lacZ operon from the {sigma}32-dependent htpG promoter, and transformants screened for higher than average {sigma}32 activity by plating on triphenyltetrazolium lactose indicator plates and selecting white colonies (high-lactose fermenters). Candidates were confirmed by restreaking on MacConkey lactose plates (selecting red colonies) and then on Xgal and IPTG (selecting dark blue colonies). Thirty dark-blue colonies were selected at random, the ASKA plasmids extracted and sequenced to determine which ORF was responsible for inducing the HSR (Table 3), and the level of induction of the HSR measured by overexpressing each of the candidates in strains grown in liquid culture (Supplementary Fig. 2). Twenty-five candidates induced the HSR, and surprisingly, ~60% of these were membrane proteins; this is a significant overrepresentation, given that only 20% of all E. coli proteins are located in the membrane (Serres et al. 2004Go). Not all overexpressed membrane proteins induced the HSR (e.g., LacY) (Supplementary Fig. 2).


View this table:
[in this window]
[in a new window]
 
Table 3. Inducers of the heat-shock response

 

    Discussion
 Top
 Abstract
 Results
 Discussion
 Materials and methods
 Acknowledgments
 References
 
The ability of cells to maintain homeostasis in response to temperature stress is essential for viability. In E. coli, {sigma}32-controlled genes govern the immediate response to temperature stress. In this report we used global approaches to provide complete or nearly complete identification of the {sigma}32 regulon and to identify endogenous signals of the response. Our analysis, together with data in the literature, indicates that the multifaceted {sigma}32 stress response not only maintains protein homeostasis but protects DNA and RNA and is further propagated because {sigma}32 controls the expression of other global regulators. Finally, our results reveal an unanticipated connection with the cellular membrane: {sigma}32 senses imbalances in membrane proteins and {sigma}32-controlled genes contribute to membrane homeostasis.


The {sigma}32 regulon

We identified 29 additional {sigma}32 promoters and 57 additional members of the {sigma}32 regulon. At present, the regulon is comprised of 49 TUs regulated by 51 {sigma}32-dependent promoters, which together encode on the chromosome 89 ORFs and one rRNA, and one ORF on the F factor (Table 1). We suggest that our analysis has identified most, if not all, {sigma}32-dependent TUs present on the arrays (~96% of the genome). This assertion is based on our success in validating previously proposed {sigma}32-dependent TUs. Of the 25 of the previously proposed TUs, 23 were present on our arrays; of these, 16 showed rapid induction after overexpression of {sigma}32. The seven uninduced TUs also failed to generate {sigma}32-dependent transcripts either in vitro or in vivo and contained poorly conserved {sigma}32 promoters (Table 1C). We therefore suggest that these seven promoters either require additional regulators or are not recognized by {sigma}32. Thus, every bona fide previously proposed {sigma}32 TUs present on our array was validated in this analysis.

The consensus promoter sequence derived from the 50 validated {sigma}32 promoters indicates the following highly conserved bases: STTGAAA-N11–12-GNCCCCATWT (Fig. 4; note the optimal spacer length [13–14 bp] is from the core-10 motif, CCCCATWT). The importance of these conserved bases is indicated by their close correspondence to the "functional" {sigma}32-promoter sequence, CTTGA-14 bp-GNCCCCATNT, derived by mutating each position of the core groE {sigma}32-promoter sequence to every other nucleotide (Wang and deHaseth 2003Go). There are subtle differences in the less-preferred bases between our aligned promoter profile and the derived "functional" {sigma}32-promoter sequence, suggesting that the importance of individual bases may vary according to the context of the local sequence.

We identified five overlapping {sigma}32/{sigma}70 promoters in addition to the previously reported rrnBP1 {sigma}32/{sigma}70 dual promoter (Newlands et al. 1993Go). Given the similarity between the {sigma}32 and {sigma}70 core promoter sequences ({sigma}70 TTGACA-N17-TATAAT; {sigma}32 TTGAAA-N14-CCCCAT WT) it is likely that both holoenzymes are recognizing the same core promoter sequences. However, as 22 of our {sigma}32 promoters generated only {sigma}32-specific transcripts as detected by 5'RACE, it is unlikely that {sigma}70 holoenzyme directly transcribes many more {sigma}32 promoters, (Table 1A; Fig. 2). However, in Pseudomonas putida, both {sigma}S and {sigma}32 holoenzymes recognize the Pm promoter (Dominguez-Cuevas et al. 2005Go), raising the possibility for promoter overlap between {sigma}S and {sigma}32 in E. coli.

Many significantly induced genes did not have a detectable {sigma}32 promoter. Interestingly, of the 47 rapidly induced genes (Fig. 1C), 43 were under {sigma}32 control whether or not they were strongly induced, whereas most induced genes without an identified {sigma}32 promoter (18/25) were only slowly and weakly induced (Table 1; Supplementary Fig. 3). Since {sigma}32 controls a number of transcriptional regulators (see below), these transcription factors may be responsible for delayed expression of genes lacking a {sigma}32 promoter. These data suggest that even in bacteria where events are very rapid, careful temporal analysis may be able to resolve cascades of induction.


{sigma}32-mediated temperature adaptation has many cellular targets

A primary role of the {sigma}32 regulon in protein homeostasis was suggested by the fact that >50% of previously known regulon members were either chaperones or proteases (Table 1A). Our results suggest that the regulon also contributes to protein homeostasis by maintaining an adequate supply of complex proteins that contain cofactors, Fe++-S centers, or lipoyl modification. Since the extent to which damaged complex proteins can be repaired is not well known, thermolabile proteins of this type may need to be synthesized at higher rates during temperature stress. The regulon encodes two Fe++-S center proteins (SdaA, IspH), a predicted Fe++-S center protein (YggW), and an IscA homolog (GntY) (Giel et al. 2006Go) likely to be involved in Fe++-S cluster biogenesis. The regulon also encodes lipoyl protein ligase (LipB) and enzymes participating in the biosynthesis of the cofactors riboflavin (RibE), thiamin (ThiL), folate (FolP), and probably biotin (YafE). If cofactor supply is tightly coupled to demand, these pathways may need to be induced during temperature stress. The particular enzymes in the regulon are likely to be rate-limiting and highly regulated, either because they are at the junction of two pathways (RibE and FolP) or because they utilize ATP (ThiL). Thus, altering the amount of these enzymes is likely to influence the amount of cofactor available to the cell.

{sigma}32-regulon members protect both DNA and RNA (Table 2B). Maintaining genomic integrity is essential for survival of the organism and high temperatures may increase both misincorporation and chromosome damage. The regulon encodes enzymes that mediate mismatch and excision repair (MutL, MutM), general recombination functions that resolve double-strand breaks (RecA, RecJ), and a site-specific recombination system to resolve chromosome dimers (XerD). This response also covalently modifies both classes of stable RNAs, possibly improving their thermal resistance. {sigma}32-regulon members have previously been shown to methylate 23S RNA (FtsJ) (Bugl et al. 2000Go) and to modify tRNA by transferring {Delta}3-isopentyl-PP to A37 of tRNA (MiaA; just adjacent to and 3' of the anticodon) (Tsui et al. 1996Go). A new regulon member, IspH, also produces {Delta}3-isopentyl-PP; this tRNA modification stabilizes codon–anticodon pairing (Bjork and Hagervall 2005Go) and cells lacking this modification are temperature sensitive (Tsui et al. 1996Go). Finally, the regulon encodes the RNA chaperone Hfq (Table 2B) (Chuang and Blattner 1993Go; Tsui et al. 1996Go), which is required for the proper function of small RNAs.

Several regulon functions may promote transcription and translation at high temperatures. Two newly discovered regulon members, HepA and NusB, join TopA (Qi et al. 1996Go) as general effectors of transcription, and all three may respond to altered supercoiling at high temperature. Previous work indicated that thermal induction of TopA, an RNA polymerase-associated topoisomerase (Cheng et al. 2003Go) is necessary for normal transcription patterns at high temperature and for thermal resistance (Qi et al. 1996Go). HepA (RapA) is also an RNA polymerase-associated protein (Sukhodolets et al. 2001Go) specifically required for RNA polymerase recycling of tightly compacted negatively supercoiled DNA. Induction of NusB, essential for antitermination of rRNA transcription (Zellars and Squires 1999Go; Torres et al. 2004Go) may reflect altered requirements for rRNA antitermination when DNA supercoiling is altered, possibly to decrease R-loop formation (Drolet 2006Go). Finally, a general connection of the regulon to translation was previously established by Korber et al. (2000)Go, who showed that regulon member Hsp15 (YfrH) binds free 50S ribosomal subunits with nascent chains and may promote ribosome recycling. This state is more prevalent at high temperatures (Korber et al. 2000Go), possibly because translation often terminates prematurely at lethal temperatures (VanBogelen and Neidhardt 1990Go).

{sigma}32 is also a master regulator controlling expression of seven transcriptional factors (including RpoD) (Table 2B) in addition to functions that alter RNA polymerase action (described above). The proportion of the {sigma}32 regulon devoted to transcriptional factors (~8%) is similar to that for {sigma}S, previously thought to be the only alternative E. coli {sigma} with a propagated response (Weber et al. 2005Go). The transcriptional regulators under {sigma}32 control mostly sense extracellular conditions: Mlc senses glucose, Pho P/Q senses low Mg++ and is induced in acid stress, ZntR responds to excess Zn++, and NarP induction requires nitrite/nitrate under anaerobic conditions. Intriguingly, many of the genes controlled by these regulators are membrane-localized transporters and other membrane proteins, pointing to the connection of this response to the status of the cytoplasmic membrane.

A comparison of genes directly regulated by HSF (heat-shock factor; controlling the eukaryotic counterpart response) with the {sigma}32-mediated heat-shock response indicates that only chaperones are encoded by both regulons (Hahn et al. 2004Go; Yamamoto et al. 2005Go). However, both regulons are of similar size and there is definite overlap in functional classes induced; for example, proteolysis, small molecule transport, additional transcription factors, proteins that modify RNA polymerase, and DNA related proteins. A broad response to heat stress appears to be universal among organisms.


The cytoplasmic membrane is a target of the {sigma}32-mediated hsr

The cytoplasmic membrane maintains cellular integrity and coordinates processes that traverse cellular compartments including secretion and processing of envelope proteins, energy generation, sensing the extracytoplasmic environment, and synthesis and/or transfer of the building blocks for phospholipids, LPS, and peptidoglycan. Previously identified {sigma}32 regulon members FtsH (Herman et al. 1995Go), HflX, HflC (Chuang and Blattner 1993Go; Tsui et al. 1996Go), and HtpX (Kornitzer et al. 1991Go) play a vital role in membrane quality control; FtsH, together with its regulators HflX and HflC, degrades and dislocates unassembled membrane proteins, and the function of HtpX overlaps with that of FtsH. Our new data indicates that {sigma}32-regulon members have important additional roles in membrane homeostasis (Table 2):

  1. Increasing the potential to make fatty acids by up-regulating expression of carbonic anhydrase (YadF), which converts C02 to bicarbonate (Merlin et al. 2003Go). The first committed step of fatty acid biosynthesis is a major consumer of bicarbonate (Merlin et al. 2003Go).
  2. Altering transport properties of the membrane by up-regulating several transporters (Table 2B).
  3. Maintaining the potential for disulfide bond formation and isomerization by increasing menaquinone production (MenA), which is necessary to reoxidize DsbB under anaerobic conditions (Bader et al. 1999Go) and by up-regulating the disufide bond isomerase DsbC.
  4. Maintaining lipoprotein maturation by up-regulating both signal peptidase II (IspA), which cleaves the proprotein, and CutE, the enzyme that activates lipoproteins by adding palmitate to their N-terminal cysteine (Wu 1996Go). Because there are >100 lipoproteins, the precise physiological role of this regulation is unclear.
  5. Altering the composition of the lipid bilayer by up-regulating PgpA, an enzyme involved in synthesis of cardiolipin, a minority phospholipid in the inner membrane whose abundance increases during stationary phase (Hiraoka et al. 1993Go). As cells induce a variety of substances during stationary phase to prevent cell lysis, cardiolipin may impart additional stability to the lipid bilayer.
  6. Altering sensing capacity of the membrane by inducing the expression of two membrane-localized histidine kinases (CreC and PhoQ) and their cytoplasmic response regulators.

 These results strongly implicate {sigma}32-regulon members in maintaining membrane functionality and homeostasis. This connection rationalizes why the {sigma}E-envelope stress response regulates {sigma}32 in all organisms examined thus far, making this one of the most conserved features of the response (Rhodius et al. 2006Go). Because the inner membrane controls the flow of building blocks and proteins to the envelope, induction of {sigma}32 is essential to maintaining envelope homeostasis.

It is provocative that genome-wide analysis indicated that when overproduced, very few cytoplasmic proteins and a disproportionate number of membrane proteins induce the {sigma}32 response. Whether this simply reflects an increased propensity for accumulated membrane proteins to misfold or whether a novel mechanism is involved, the fact that {sigma}32 is particularly sensitive to imbalances in membrane proteins provides a rationale for utilizing FtsH to degrade {sigma}32. Because FtsH carries out both membrane quality control and degradation of {sigma}32, these two processes can be coordinated, thereby providing the cell with a means to assess membrane homeostasis. Interestingly, the overexpressed membrane proteins that induce the {sigma}32-mediated response also induce the response mediated by the two-component CpxRA system (C. Herman, unpubl.). Moreover, at least two proteins, HtpX and DsbC, are jointly induced by both CpxRA and {sigma}32. Additional evidence that CpxRA senses membrane status has recently been presented by Ito (Shimohata et al. 2002Go). Thus, {sigma}32 and CpxRA may jointly protect the cell membrane during times of stress.


    Materials and methods
 Top
 Abstract
 Results
 Discussion
 Materials and methods
 Acknowledgments
 References
 

Medium, strains, and plasmids

M9 medium (Sambrook et al. 1989Go) was supplemented with 0.2% glucose, 1 mM MgSO4, vitamins, all amino acids (40 µg/mL), ampicillin (100 µg/mL), and/or kanamycin (30 µg/mL) as required and is referred to as M9 complete. Bacterial strains and plasmids used in this study are listed in Table 4.


View this table:
[in this window]
[in a new window]
 
Table 4. Bacterial strains and plasmids used in this study

 

Pulse chase

A 750-µL aliquot of cells was removed from a culture growing exponentially in M9 complete–methionine at 30°C, labeled for 1 min with 80 µC [35S]methionine, chased for 1 min with excess unlabeled methionine, precipitated with 100 µL of ice-cold 50% TCA (15 min on ice), centrifuged, and resuspended in 50 µL of 2% SDS and 50 mM Tris (pH 7.5). Equal counts were loaded in each lane of a 7.5% polyacrylamide gel and bands corresponding to GroEL and DnaK proteins were quantified using the program ImageQuant version 1.2.


Time course microarray experiments

Cells were inoculated to OD450 = 0.02 from fresh overnight cultures of CAG50002 (Table 4; vector has IPTG-inducible rpoH) or CAG50001 (Table 4; control vector) into 100 mL of M9 complete in 500-mL conical flasks and grown aerobically at 30°C in a gyratory water bath (model G76 from New Brunswick Scientific) shaking at 240 rpm. At OD450 = 0.3, cultures were induced with IPTG (1 mM final concentration). Ten-milliliter samples for microarray analysis were removed immediately prior to and at 2.5, 5, 10, and 20 min after induction, added to ice-cold 5% water–saturated phenol in ethanol, and centrifuged at 6600g. Cell pellets were flash frozen in liquid N2 and stored at –80°C. Preparation of labeled probes and microarray procedures were performed exactly as described in Rhodius et al. (2006)Go and Supplemental Material. Raw and normalized microarray expression data are available on the NCBI Gene Expression Onmibus (GEO) Web site (http://www.ncbi.nih.gov/geo) under the accession code GSE4321 [NCBI GEO] .

Additional Materials and Methods are provided in Supplemental Material.


    Acknowledgments
 Top
 Abstract
 Results
 Discussion
 Materials and methods
 Acknowledgments
 References
 
We thank Chi Zen Lu for supplying purified RNAP core, {sigma}32, and {sigma}70 proteins, H. Mori from the Japanese E. coli consortium for the generous gift of the ASKA plasmid ORF library, and P. Kiley for advice about FeS proteins and assembly of complex proteins. This work was supported by Ajinomoto Co., Inc., and National Institutes of Health (NIH) Grants GM57755 and GM36278 (to C.A.G.).


    Footnotes
 
3 Present addresses:

Ajinomoto Co., Inc., 15-1, Kyobashi 1-chome, Chuo-ku, Tokyo 104-8315, Japan; Back

4 Department of Molecular and Human Genetics, Baylor College of Medicine, One Baylor Plaza, Houston, TX 77030, USA. Back

5 Corresponding author.

E-MAIL cgross{at}cgl.ucsf.edu; FAX (415) 514-4080. Back

Supplemental material is available at http://www.genesdev.org.

Article is online at http://www.genesdev.org/cgi/doi/10.1101/gad.1428206


    References
 Top
 Abstract
 Results
 Discussion
 Materials and methods
 Acknowledgments
 References
 
Bader M., Muse W., Ballou D.P., Gassner C., Bardwell J.C. 1999. Oxidative protein folding is driven by the electron transport system. Cell 98: 217–227.[CrossRef][Medline]

Bardwell J.C. and Craig E.A. 1984. Major heat shock gene of Drosophila and the Escherichia coli heat-inducible dnaK gene are homologous. Proc. Natl. Acad. Sci. 81: 848–852.[Abstract/Free Full Text]

Bardwell J.C. and Craig E.A. 1987. Eukaryotic Mr 83,000 heat shock protein has a homologue in Escherichia coli.. Proc. Natl. Acad. Sci. 84: 5177–5181.[Abstract/Free Full Text]

Beckmann R.P., Mizzen L.E., Welch W.J. 1990. Interaction of Hsp 70 with newly synthesized proteins: Implications for protein folding and assembly. Science 248: 850–854.[Abstract/Free Full Text]

Biran D., Brot N., Weissbach H., Ron E.Z. 1995. Heat shock-dependent transcriptional activation of the metA gene of Escherichia coli.. J. Bacteriol. 177: 1374–1379.[Abstract/Free Full Text]

Bjork G. and Hagervall T. 2005. Transfer RNA Modification. In Escherichia coli and Salmonella: Cellular and molecular biology. In Module 4.6.2 (ed. Curtis III R.) . . ASM Press, Washington, D.C.

Blaszczak A., Georgopoulos C., Liberek K. 1999. On the mechanism of FtsH-dependent degradation of the {sigma} 32 transcriptional regulator of Escherichia coli and the role of the Dnak chaperone machine. Mol. Microbiol. 31: 157–166.[CrossRef][Medline]

Bugl H., Fauman E.B., Staker B.L., Zheng F., Kushner S.R., Saper M.A., Bardwell J.C., Jakob U. 2000. RNA methylation under heat shock control. Mol. Cell 6: 349–360.[CrossRef][Medline]

Burden S., Lin Y.X., Zhang R. 2005. Improving promoter prediction for the NNPP2.2 algorithm: A case study using Escherichia coli DNA sequences. Bioinformatics 21: 601–607.[Abstract/Free Full Text]

Charpentier B. and Branlant C. 1994. The Escherichia coli gapA gene is transcribed by the vegetative RNA polymerase holoenzyme E {sigma} 70 and by the heat shock RNA polymerase E {sigma} 32. J. Bacteriol. 176: 830–839.[Abstract/Free Full Text]

Cheng B., Zhu C.X., Ji C., Ahumada A., Tse-Dinh Y.C. 2003. Direct interaction between Escherichia coli RNA polymerase and the zinc ribbon domains of DNA topoisomerase I. J. Biol. Chem. 278: 30705–30710.[Abstract/Free Full Text]

Chin D.T., Goff S.A., Webster T., Smith T., Goldberg A.L. 1988. Sequence of the lon gene in Escherichia coli. A heat-shock gene which encodes the ATP-dependent protease La. J. Biol. Chem. 263: 11718–11728.[Abstract/Free Full Text]

Chuang S.E. and Blattner F.R. 1993. Characterization of twenty-six new heat shock genes of Escherichia coli.. J. Bacteriol. 175: 5242–5252.[Abstract/Free Full Text]

Chuang S.E., Burland V., Plunkett III G., Daniels D.L., Blattner F.R. 1993a. Sequence analysis of four new heat-shock genes constituting the hslTS/ibpAB and hslVU operons in Escherichia coli.. Gene 134: 1–6.[CrossRef][Medline]

Chuang S.E., Daniels D.L., Blattner F.R. 1993b. Global regulation of gene expression in Escherichia coli.. J. Bacteriol. 175: 2026–2036.[Abstract/Free Full Text]

Cowing D.W., Bardwell J.C., Craig E.A., Woolford C., Hendrix R.W., Gross C.A. 1985. Consensus sequence for Escherichia coli heat shock gene promoters. Proc. Natl. Acad. Sci. 82: 2679–2683.[Abstract/Free Full Text]

Crooks G.E., Hon G., Chandonia J.M., Brenner S.E. 2004. WebLogo: A sequence logo generator. Genome Res. 14: 1188–1190.[Abstract/Free Full Text]

Dartigalongue C. and Raina S. 1998. A new heat-shock gene, ppiD, encodes a peptidyl-prolyl isomerase required for folding of outer membrane proteins in Escherichia coli.. EMBO J. 17: 3968–3980.[CrossRef][Medline]

Dominguez-Cuevas P., Marin P., Ramos J.L., Marques S. 2005. RNA polymerase holoenzymes can share a single transcription start site for the Pm promoter. Critical nucleotides in the –7 to –18 region are needed to select between RNA polymerase with {sigma}38 or {sigma}32. J. Biol. Chem. 280: 41315–41323.[Abstract/Free Full Text]

Drolet M. 2006. Growth inhibition mediated by excess negative supercoiling: The interplay between transcription elongation, R-loop formation and DNA topology. Mol. Microbiol. 59: 723–730.[CrossRef][Medline]

El-Samad H., Kurata H., Doyle J.C., Gross C.A., Khammash M. 2005. Surviving heat shock: Control strategies for robustness and performance. Proc. Natl. Acad. Sci. 102: 2736–2741.[Abstract/Free Full Text]

Estrem S.T., Gaal T., Ross W., Gourse R.L. 1998. Identification of an UP element consensus sequence for bacterial promoters. Proc. Natl. Acad. Sci. 95: 9761–9766.[Abstract/Free Full Text]

Gaitanaris G.A., Papavassiliou A.G., Rubock P., Silverstein S.J., Gottesman M.E. 1990. Renaturation of denatured {lambda} repressor requires heat shock proteins. Cell 61: 1013–1020.[CrossRef][Medline]

Gamer J., Multhaup G., Tomoyasu T., McCarty J.S., Rudiger S., Schonfeld H.J., Schirra C., Bujard H., Bukau B. 1996. A cycle of binding and release of the DnaK, DnaJ and GrpE chaperones regulates activity of the Escherichia coli heat shock transcription factor {sigma}32. EMBO J. 15: 607–617.[Medline]

Giel J., Rodionov D., Liu M., Blattner F.R., Kiley P.J. 2006. IscR-dependent gene expression links iron–sulfur assembly to the control of O2-regulated genes in Escehrichia coli.. Mol. Microbiol. 60: 1058–1075.[CrossRef][Medline]

Grossman A.D., Erickson J.W., Gross C.A. 1984. The htpR gene product of E. coli is a {sigma} factor for heat-shock promoters. Cell 38: 383–390.[CrossRef][Medline]

Guisbert E., Herman C., Lu C.Z., Gross C.A. 2004. A chaperone network controls the heat shock response in E. coli.. Genes & Dev. 18: 2812–2821.[Abstract/Free Full Text]

Guyer M.S., Reed R.R., Steitz J.A., Low K.B. 1981. Identification of a sex-factor-affinity site in E. coli as {gamma} {delta}. Cold Spring Harb. Symp. Quant. Biol. 45: 135–140.[Medline]

Hahn J.S., Hu Z., Thiele D.J., Iyer V.R. 2004. Genome-wide analysis of the biology of stress responses through heat shock transcription factor. Mol. Cell. Biol. 24: 5249–5256.[Abstract/Free Full Text]

Herman C., Thevenet D., D’Ari R., Bouloc P. 1995. Degradation of {sigma} 32, the heat shock regulator in Escherichia coli, is governed by HflB. Proc. Natl. Acad. Sci. 92: 3516–3520.[Abstract/Free Full Text]

Hiraoka S., Matsuzaki H., Shibuya I. 1993. Active increase in cardiolipin synthesis in the stationary growth phase and its physiological significance in Escherichia coli.. FEBS Lett. 336: 221–224.[CrossRef][Medline]

Jakob U., Muse W., Eser M., Bardwell J.C. 1999. Chaperone activity with a redox switch. Cell 96: 341–352.[CrossRef][Medline]

Jensen K.F. 1993. The Escherichia coli K-12 ‘wild types’ W3110 and MG1655 have an rph frameshift mutation that leads to pyrimidine starvation due to low pyrE expression levels. J. Bacteriol. 175: 3401–3407.[Abstract/Free Full Text]

Kanemori M., Mori H., Yura T. 1994. Effects of reduced levels of GroE chaperones on protein metabolism: Enhanced synthesis of heat shock proteins during steady-state growth of Escherichia coli.. J. Bacteriol. 176: 4235–4242.[Abstract/Free Full Text]

Kitagawa M., Wada C., Yoshioka S., Yura T. 1991. Expression of ClpB, an analog of the ATP-dependent protease regulatory subunit in Escherichia coli, is controlled by a heat shock {sigma} factor ({sigma} 32). J. Bacteriol. 173: 4247–4253.[Abstract/Free Full Text]

Kitagawa M., Ara T., Arifuzzaman M., Ioka-Nakamichi T., Inamoto E., Toyonaga H., Mori H. 2005. Complete set of ORF clones of Escherichia coli ASKA library (A complete set of E. coli K-12 ORF archive): Unique resources for biological research. DNA Res. 12: 291–299.[Abstract/Free Full Text]

Korber P., Stahl J.M., Nierhaus K.H., Bardwell J.C. 2000. Hsp15: A ribosome-associated heat shock protein. EMBO J. 19: 741–748.[CrossRef][Medline]

Kornitzer D., Teff D., Altuvia S., Oppenheim A.B. 1991. Isolation, characterization, and sequence of an Escherichia coli heat shock gene, htpX. J. Bacteriol. 173: 2944–2953.[Abstract/Free Full Text]

Lemaux P.G., Herendeen S.L., Bloch P.L., Neidhardt F.C. 1978. Transient rates of synthesis of individual polypeptides in E. coli following temperature shifts. Cell 13: 427–434.[CrossRef][Medline]

Lesley S.A., Jovanovich S.B., Tse-Dinh Y.C., Burgess R.R. 1990. Identification of a heat shock promoter in the topA gene of Escherichia coli.. J. Bacteriol. 172: 6871–6874.[Abstract/Free Full Text]

Li C., Tao Y.P., Simon L.D. 2000. Expression of different-size transcripts from the clpP-clpX operon of Escherichia coli during carbon deprivation. J. Bacteriol. 182: 6630–6637.[Abstract/Free Full Text]

Liberek K., Galitski T.P., Zylicz M., Georgopoulos C. 1992. The DnaK chaperone modulates the heat shock response of Escherichia coli by binding to the {sigma}32 transcription factor. Proc. Natl. Acad. Sci. 89: 3516–3520.[Abstract/Free Full Text]

Lipinska B., King J., Ang D., Georgopoulos C. 1988. Sequence analysis and transcriptional regulation of the Escherichia coli grpE gene, encoding a heat shock protein. Nucleic Acids Res. 16: 7545–7562.[Abstract/Free Full Text]

Mecsas J., Rouviere P.E., Erickson J.W., Donohue T.J., Gross C.A. 1993. The activity of {sigma} E, an Escherichia coli heat-inducible {sigma}-factor, is modulated by expression of outer membrane proteins. Genes & Dev. 7: 2618–2628.[Abstract/Free Full Text]

Merlin C., Masters M., McAteer S., Coulson A. 2003. Why is carbonic anhydrase essential to Escherichia coli? J. Bacteriol. 185: 6415–6424.[Abstract/Free Full Text]

Missiakas D. and Raina S. 1997. Signal transduction pathways in response to protein misfolding in the extracytoplasmic compartments of E. coli: Role of two new phosphoprotein phosphatases PrpA and PrpB. EMBO J. 16: 1670–1685.[CrossRef][Medline]

Missiakas D., Georgopoulos C., Raina S. 1993. The Escherichia coli heat shock gene htpY: Mutational analysis, cloning, sequencing, and transcriptional regulation. J. Bacteriol. 175: 2613–2624.[Abstract/Free Full Text]

Morita M., Kanemori M., Yanagi H., Yura T. 1999a. Heat-induced synthesis of {sigma}32 in Escherichia coli: Structural and functional dissection of rpoH mRNA secondary structure. J. Bacteriol. 181: 401–410.[Abstract/Free Full Text]

Morita M.T., Tanaka Y., Kodama T.S., Kyogoku Y., Yanagi H., Yura T. 1999b. Translational induction of heat shock transcription factor {sigma}32: Evidence for a built-in RNA thermosensor. Genes & Dev. 13: 655–665.[Abstract/Free Full Text]

Newlands J.T., Gaal T., Mecsas J., Gourse R.L. 1993. Transcription of the Escherichia coli rrnB P1 promoter by the heat shock RNA polymerase (E {sigma} 32) in vitro. J. Bacteriol. 175: 661–668.[Abstract/Free Full Text]

Pelham H.R. 1986. Speculations on the functions of the major heat shock and glucose-regulated proteins. Cell 46: 959–961.[CrossRef][Medline]

Qi H., Menzel R., Tse-Dinh Y.C. 1996. Effect of the deletion of the {sigma} 32-dependent promoter (P1) of the Escherichia coli topoisomerase I gene on thermotolerance. Mol. Microbiol. 21: 703–711.[CrossRef][Medline]

Raina S. and Georgopoulos C. 1990. A new Escherichia coli heat shock gene, htrC, whose product is essential for viability only at high temperatures. J. Bacteriol. 172: 3417–3426.