![]()
|
|
||||||||
J. Biol. Chem., Vol. 282, Issue 15, 11230-11237, April 13, 2007
| ||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||





1
From the
Department of Molecular Biology and Biotechnology, The University of Sheffield, Sheffield S10 2TN, United Kingdom and the
Department of Computer Science, The University of Sheffield, Sheffield S1 4DP, United Kingdom
Received for publication, January 25, 2007 , and in revised form, February 15, 2007.
| ABSTRACT |
|---|
|
|
|---|
| INTRODUCTION |
|---|
|
|
|---|
Adaptation to changes in O2 availability is regulated at the level of transcription by two well characterized systems, FNR and ArcBA (6-9). FNR is a direct O2 sensor (10), whereas the ArcBA two-component system senses O2 availability indirectly by monitoring the redox state of the quinone pool (11). FNR and ArcBA cooperate to regulate some operons, such as cydAB (encoding cytochrome bd) and focA-pflB (encoding PFL) (5, 12). Under aerobic conditions expression of cydAB is minimal, but as O2 availability decreases, ArcA binds at the cydAB promoter to facilitate maximal expression under micro-aerobic conditions (5). As O2 availability decreases further, FNR binds the cydAB promoter to repress expression (5). Similarly, under aerobic conditions expression of pflB is minimal, but as O2 availability decreases, ArcA and FNR bind to the focA-pflB promoter region to activate expression (12). Thus, the transcriptional outputs are dependent on the relative amounts of active transcription factors (in these cases FNR and ArcA) and their regulatory strengths. It is difficult to determine regulatory activities of transcription factors experimentally, but they can be inferred using models where information on network connectivity and time-resolved transcript profiling is available. Such models promise to provide new insight into the responses of regulatory networks in changing environments (13).
In the present work, changes in the E. coli transcriptome were monitored as aerobic chemostat cultures adapted to micro-aerobic conditions. It is shown that, as well as predictable changes in transcripts encoding proteins involved in central metabolic processes, the bacteria respond to changes in metal ion physiology and to cell envelope stress during the transition to micro-aerobic growth. Modeling transcriptome dynamics revealed that the responses to lower O2 availability and cell envelope stress involved fast and slow reacting components.
| EXPERIMENTAL PROCEDURES |
|---|
|
|
|---|
For competitive index experiments approximately equal numbers of parent (W3110) and yfiD (JRG4403; 15) strains were mixed and used to inoculate (1:100) 10 ml of Evansdefined medium (14) with glucose (30 mM) as the carbon source. The cultures were incubated at 37 °C under 20, 7, or 0% O2 atmospheres for 16 h. The number of colony forming units (cfu) on L agar and L agar containing kanamycin (30 mg l-1) were then measured. Competitive index (CI) values were calculated as follows: mutantanaerobic/parentanaerobic: mutantaerobic/parentaerobic and mutantmicroaerobic/parentmicroaerobic: mutantaerobic/parentaerobic.
The effects of mutations in key regulators on gene expression were determined using promoter-lacZ fusions (16) of cultures grown under the indicated conditions in L broth (17) according to Miller (18). Oxygen transfer to cultures was measured as described (19). The phenotype of the ompW mutant was investigated using Biolog phenotype microarrays (Biolog, Hayward CA).
Transcript ProfilingTranscript profiling and data analysis were as described by Partridge et al. (20). Briefly culture samples (10 ml) were flushed from the chemostat directly into RNA Protect (Qiagen) before total RNA purification using a RNeasy Midi kit (Qiagen). Equal quantities of RNA from control (aerobic steady state) and experimental (5, 10, 15, and 60 min lowered exposure to air) cultures were labeled using Cy3-dCTP or Cy5-dCTP. For each E. coli K12 OciChip (Ocimum Biosolutions, product number 2140-000000), dye swap experiments were done for each pairing with two biological replicates (separate chemostat runs).
Data were analyzed using Imagene, version 5.1 and Genesight version 4 (Biodiscovery Inc) as previously described (20). Genes showing greater than 2-fold changes in transcript abundance at one or more time points with a p value of
0.05 (as determined by a Student's t test) were deemed to be differentially regulated.
Real Time PCRReactions and control samples (using the gene-specific primers listed in supplemental Table S2) were assembled and analyzed as previously described (20).
Transcript Mapping by 5'-RACE PCRRNA was isolated from anaerobic cultures using Qiagen RNeasy mini kits according to the manufacturer's instructions. The ompW and ftnB transcript starts were identified using 2 µg of E. coli RNA per RACE reaction according to the manufacturer's instructions (Roche Applied Science) with the oligonucleotides listed in supplemental Table S2.
Modeling Transcription Factor ActivitiesThe response of the regulatory network during the transition from aerobic to micro-aerobic conditions was analyzed using a probabilistic model (13). The model adopts a log-linear approximation to the kinetics of the regulatory system where the transcription profiles are modeled as sparse linear combinations of the transcription factor activity profiles. Assumptions are incorporated as prior probability distributions on the active transcription factor profiles and the gene-specific intensity with which they regulate their targets. These are, essentially, a continuity assumption on the activity profiles (encapsulated in the Markov chain structure of the prior distribution) and equal prior probability of repression and activation. The sparsity of the model is provided by the structure of the regulatory network. The gene expression data is used to infer transcription factor activity profiles and the gene specific regulatory intensities using Bayes' theorem. A more detailed description is provided under supplemental materials.
| RESULTS AND DISCUSSION |
|---|
|
|
|---|
210 µM dissolved O2) to microaerobic conditions (O2 transfer rate to the culture of 12 ± 0.5 mmol min-1 liter-1). The choice of glucose as the carbon source will ultimately allow the full range of metabolic modes (aerobic respiration, anaerobic respiration and fermentation) to be investigated without the need to change the carbon source, but will affect the responses of transcripts such as ansB, that are co-dependent on the cyclic-AMP receptor protein for expression. The dilution rate chosen (0.2 h-1) was lower than the µmax under micro-aerobic conditions, thus allowing the cultures the potential to respond to changes in O2 availability without changes in growth rate. Nevertheless, any temporary changes in growth rate during adaptation to the new conditions forms an integral component of the response and will be reflected in the transcriptome. Samples from the aerobic steady-state culture were rapidly flushed into RNA protect (Qiagen) for isolation of total RNA. Switching off the air that was sparging the culture perturbed the aerobic steady-state. Samples were obtained for RNA preparation 5, 10, 15, and 60 min after the air supply to the culture was lowered. The corresponding labeled cDNAs were synthesized and used to probe E. coli K12 arrays (Ocimum Biosolutions). For initial analysis of the data a significant change in transcript abundance was defined as
2-fold (p
0.05) at one or more of the time points analyzed. The full dataset is available (GEO accession number GSE6644
[NCBI GEO]
).
|
|
The phenotype of the yfiD mutant was analyzed in competitive index experiments, to determine whether the changes in transcript abundance were reflected in culture physiology. An index of <1.0 indicates that the mutant is impaired compared with the parent. For micro-aerobic conditions, the value was 0.12 ± 0.03; for anaerobic conditions the value was 0.65 ± 0.11 (mean ± S.D., n = 6). Thus, growth of the yfiD mutant was significantly impaired under micro-aerobic conditions compared with the parent, suggesting that the changes in transcription are physiologically relevant.
FNR and ArcA are the major regulators responsible for adaptation of central metabolism to changes in O2 availability (Table 1). To gain further insight into this regulatory network, a probabilistic model was used to infer FNR and ArcA activity profiles during the transition from aerobic to micro-aerobic conditions. The activity profiles show that FNR (Fig. 1B, blue line) responds rapidly to the change in culture aeration, whereas ArcA (Fig. 1B, red line) responds more slowly. Furthermore, it is inferred from the model that FNR activity peaks after 10 min and then significantly decreases as time proceeds, while the activity of ArcA increases up to 15 min and then stabilizes. The delay in the ArcBA response compared with FNR might reflect the signals perceived by these two transcription factors; FNR is a direct O2 sensor and responds rapidly, whereas ArcBA monitors O2 indirectly, via the redox state of the quinone pool, and responds more slowly. Moreover, arcA expression is regulated by FNR (29).
Interruption of the O2 Supply to Aerobic Cultures Changes Metal Ion PhysiologyMetal ion homeostasis is important for maintaining viability because many metal ions are not only essential, but also toxic. Switching E. coli cultures from aerobic to micro-aerobic growth conditions increased the abundance of transcripts involved in copper management (Table 2). Under aerobic conditions the CopA/CueO system is the major contributor to copper tolerance (30). CopA removes copper from the cytoplasm to the periplasm where toxic Cu(I) is oxidized to the less toxic Cu(II) by the multi-copper oxidase, CueO, in a copper- and O2-dependent reaction. Expression of copA and cueO is regulated by CueR, but only copA responded to lower O2 availability, exhibiting maximal abundance 10 min into the switch to micro-aerobic conditions (Table 2). Thus, copA transcript abundance during the transition was similar to the inferred activity of FNR (see above). Analysis of the copA promoter revealed a potential FNR site (cTGATGCAAATCgA), with only two mismatches (shown in lowercase) compared with the consensus. This site is centered at -93.5 relative to the copA transcript start (31), suggesting that copA has a class I FNR-dependent promoter, which was confirmed by in vivo transcription using a single copy copA-lacZ fusion that responded to O2 availability in an FNR-dependent manner (supplemental Table S1). Thus, the model correctly predicted that FNR regulates copA.
|
In contrast to CopA, the CusCFBA transporter moves copper ions directly to the external milieu (30). Expression of the cusCFBA operon is controlled by the Cu(I)-responsive two-component system CusSR (31) and the cusSR transcripts are present in increased abundance (1.7-3.8-fold) after decreasing the supply of O2 (Table 2). It has previously been reported that exposure to copper causes a significant (
3-fold) increase in the amount of CusR protein in the bacterial cell (31). The activation of CusSR-regulated genes suggests that restricting the O2 supply leads to increased concentrations of Cu(I) in the periplasm. This could be caused by lower activity of the O2-dependent multicopper oxidase CueO, changes in the redox state of the periplasm, and/or lower demand for copper accompanying the switch from the copper-containing cytochrome bo3 to the copper-independent cytochrome bd under micro-aerobic conditions. The enhanced abundance of the cusCFBA transcript detected by transcript profiling was confirmed by RT-PCR of the cusF mRNA (supplemental Fig. S1). Thus, it would seem that expression of cusCFAB is enhanced under micro-aerobic conditions in the presence of only 5 µM copper (the copper concentration used in the growth medium) (Table 2), as well as under conditions of extreme copper stress (30). In a previous study the cusCFAB operon was assigned to a group of genes that are repressed by FNR, but show no O2-dependent changes in gene expression (6). A putative FNR site was identified 151.5 base pairs upstream of the transcript start (6), but it has only weak identity to, and is located further upstream than, known functional FNR sites (33). Consistent with this analysis, a cus-lacZ fusion was not regulated by FNR (6) and cus transcript abundance does not match the activity profile of FNR. Therefore, it is suggested that the increased abundance of the cusCFBA operon transcript(s) observed here, and perhaps previously (6), is mediated by the CusSR two-component system to compensate for the impairment of CueO function at low O2 concentrations, and possibly to manage the portion of the copper pool that would be sequestered in cytochrome bo3 under fully aerobic conditions.
As the concentration of O2 decreases, the solubility of iron increases. This is reflected in the enhanced abundance the feoAB transcripts encoding the ferrous ion transporter FeoAB (Table 2). Immediately downstream of feoB is an open reading frame, yhgG. There are no nucleotides between the feoB stop codon and the yhgG start codon. Thus, it is likely that there is a feoAB-yhgG operon. This is supported by the transcript profiles. The feo operon transcripts are most abundant after 10 min, and the dynamics of feo transcript abundance matched the predicted changes in FNR activity (see above). Expression of feoAB is repressed by Fur under conditions of iron sufficiency, and activated by FNR under anaerobic conditions (34). Therefore, the responses of the feo transcripts are most likely FNR-driven, allowing the increased availability of Fe(II) under conditions of lower aeration to be exploited by the bacteria.
In addition to enhanced abundance of transcripts involved in Fe(II) acquisition (feoAB-yhgG), the ftnB transcript, encoding a ferritin-like iron-storage protein, was also enhanced (Table 2). Increased abundance of the latter transcript was also observed in cultures exposed to external copper (31). Inspection of the ftnB promoter region revealed the presence of a consensus FNR binding site (TTGATTAATATCAA), suggesting that ftnB, like feoAB-yhgG might be activated by FNR under micro-aerobic conditions. However, although expression of a plasmid-based ftnB-lacZ fusion was shown to be O2 responsive, the effect of O2 was independent of FNR (not shown). Nevertheless, gel retardation analysis indicated that FNR was able to bind at the ftnB promoter region (not shown), and so the apparent lack of FNR-mediated regulation is not caused by impaired recognition of the FNR site. Transcript mapping by 5'-RACE PCR showed the FNR site was centered at -136.5 relative to the transcript start, and is thus located further upstream than any known functional FNR site (33). Thus, it is likely that FNR alone is unable to act so far upstream of RNA polymerase, and other transcription factors, responding to particular metabolic or environmental signals, might be required to allow FNR to function at this promoter. Consistent with this, ftnB transcript abundance did not correlate well with inferred FNR activity (Fig. 1B).
In summary, the changes in transcript profiles indicate that during the transition from aerobic to micro-aerobic conditions, metal ion physiology is perturbed, and that FNR plays an important role in managing the uptake and efflux of the redoxactive metals copper and iron.
Lower O2 Availability Enhances the Abundance of Transcripts Associated with Cell Envelope StressTen minutes after switching to micro-aerobic conditions, transcripts controlled by the cell envelope stress regulators CpxR or RpoE (35-37) increased in abundance (Table 3 and Fig. 2A). Modeling the activity profiles of CpxR and RpoE suggested that RpoE responds rapidly to decreased O2 availability, exhibiting maximal activity after 15 min, and lower activity after 60 min (Fig. 2B). In contrast, CpxR did not respond initially (0-5 min), but then its activity increased and stabilized after 15 min (Fig. 2B). The profiles of the transcripts in Table 3 did not match the inferred activities of any of the individual regulators analyzed above. This suggested that these transcripts are controlled by multiple transcription factors with differential effects during the transition. For example, a potential CpxR binding motif (GcAAAAAAATGTAAA) has been identified in the yfiD-ung intergenic region (38) and the activity of the FNR- and ArcA-regulated yfiD promoter was
1.5-fold lower in a cpxR mutant under micro-aerobic and anaerobic conditions (supplemental Table S1).
|
|
-barrel cell envelope protein that is probably involved in the transport of small molecules across the outer membrane, as well as being the receptor for colicin S4 (39). Biolog phenotyping of the ompW mutant showed metabolic defects when dipeptides (Gly-Asp, Gly-Glu, Gly-Pro) were supplied as carbon source, suggesting OmpW has a role in transporting these compounds. Homologs of E. coli OmpW have been implicated in bacterial adaptation to a range of stresses, including salinity, temperature, and O2 availability (40, 41). Inspection of the ompW promoter region revealed the presence of potential FNR (TTGATTTAAATCAc), CpxR (GTAAAGCGGAGgtAA), and ArcA (AGTTAATatT; gGTgAATTAT) binding sites (lowercase letters indicate deviations from the consensuses). The transcript start for ompW was located by 5'-RACE PCR placing the potential sites for FNR at -132.5, ArcA at -70, and -167 and CpxR at -224. In vivo transcription assays using an ompW-lacZ fusion plasmid in parent, fnr, cpxR, and arcA strains showed that ompW expression was activated by FNR, and repressed by ArcA and CpxR, under micro-aerobic conditions (supplemental Table S1). Comparison of the Aerobic to Micro-aerobic Transition to the Anaerobic to Aerobic TransitionAlthough the transition described here (aerobic to micro-aerobic conditions) is not an exact reversal of the anaerobic to aerobic transition that we described previously (20), some informative comparisons can be made.
First, comparing the responses of central metabolic transcripts suggests the identity of the key control points. Thus, here we observed decreased abundance of the PDHC transcript, increased abundance of the PFL transcript, and decreased abundance of the succinate dehydrogenase/succinyl-CoA synthetase transcript (Table 1). In the previous transition from anaerobic to aerobic conditions, these same transcripts also responded, but in the opposite direction (20). This common response confirms the importance of pyruvate metabolism and regulation of the sdh-suc operon in adapting central metabolism when O2 availability changes (26). Previously, introduction of O2 into anaerobic cultures caused the abundances of transcripts encoding anaerobic reductases and the micro-aerobic quinol oxidase (Cyd) to decrease (20). Here, upon lowering O2 availability, the opposite effects were observed. Moreover, analysis of the previous anaerobic-aerobic data to predict the activities of FNR and ArcA using the probabilistic model described above yielded activity profiles resembling those shown Fig. 1B, except that the activities are decreased during the transition rather than increased. Thus, FNR was again predicted to respond rapidly and overshoot, whereas ArcA responds more slowly with no overshoot (not shown).
Second, one of the transcriptional responses observed during the transition from anaerobic to aerobic conditions was the increased abundance of methionine biosynthetic transcripts (20). This response was transient and related to inhibition of MetE (cobalamin-independent methionine synthase) by oxidative stress (42). Interestingly, switching cultures from aerobic to micro-aerobic conditions also changed the abundances of the methionine biosynthetic transcripts, but in the opposite direction, i.e. they were less abundant under micro-aerobic conditions (not shown), consistent with an increase in the pool of active MetE when the O2 supply is restricted.
Third, the switch from anaerobic to aerobic conditions was accompanied by a transient increase in transcripts associated with the degradation of putrescine (20). Here, the switch from aerobic to micro-aerobic conditions caused a transient increase in the abundance of the speB transcript, encoding agmatine ureohydrolase, a putrescine biosynthetic enzyme. The speB transcript was most abundant (3.3-fold increase) 5 min after the air supply to the cultures was restricted, but after 60 min had almost returned to the levels present aerobically. Assuming that the increase in speB transcript abundance results in increased SpeB activity, then the rate of putrescine biosynthesis should be increased under micro-aerobic conditions. Thus, this response is consistent with the proposal that modulating intracellular putrescine concentrations is important in adaptation to changes in O2 availability, possibly through altering DNA topology (20).
In contrast to the common features of the two transitions described above, the copper management and cell envelope stress responses only seem to be associated with the switch to micro-aerobic conditions. This suggests that there are effective systems to manage these aspects in the aerobic and anaerobic steady states. Similarly, the transient OxyR-mediated oxidative stress response appears to be only associated with the anaerobic to aerobic transition (20). Thus, regardless of the direction of the transition (i.e. anaerobic to aerobic, or aerobic to microaerobic) central metabolism, methionine biosynthesis, and putrescine turnover appear to form part of a common core response to changes in O2 availability, but this common core is supported by additional responses to the specific challenges posed by a particular transition.
ConclusionsThe ability of bacteria to respond rapidly and effectively to environmental perturbation is a distinguishing and vital aspect of their physiology. The work described here has revealed that transfer of E. coli from aerobic to micro-aerobic growth conditions causes changes in gene expression associated with central metabolism, metal ion physiology and cell envelope stress. Probabilistic modeling of the regulator activities driving transcriptome dynamics during the transition from aerobic to micro-aerobic conditions revealed fast and slow reacting components working together to mediate adaptation. Thus, the fast-reacting direct O2 sensor, FNR, and the slow-reacting indirect O2 sensor, ArcA, combine to modulate transcription of central metabolic genes. It is suggested that systems composed of rapid and slow reacting components might be a feature of other regulatory networks in which an environmental signal is perceived directly and indirectly by two or more functionally related transcription factors controlling overlapping regulons. Such systems may have evolved to allow rapid adaptation (via the fast reacting regulator), while limiting the response to short term environmental changes by mediating a full response (via both the fast and slow reacting regulators) only when the altered conditions persist.
| FOOTNOTES |
|---|
The on-line version of this article (available at http://www.jbc.org) contains supplemental Tables S1-S3, Fig. S1, and methods. ![]()
1 To whom correspondence should be addressed: Dept. of Molecular Biology and Biotechnology, The University of Sheffield, Sheffield S10 2TN, UK. Tel.: 44-114-222-4403; Fax: 44-114-222-2800; E-mail: jeff.green{at}sheffield.ac.uk.
2 The abbreviations used are: PDHC, pyruvate dehydrogenase complex; PFL, pyruvate formate-lyase; Arc, aerobic respiratory control; cfu, colony forming unit; FNR, regulator of fumarate and nitrate reduction; RACE, rapid amplification of cDNA ends. ![]()
| ACKNOWLEDGMENTS |
|---|
| REFERENCES |
|---|
|
|
|---|
This article has been cited by other articles:
![]() |
D. M. Doughty, E. G. Kurth, L. A. Sayavedra-Soto, D. J. Arp, and P. J. Bottomley Evidence for Involvement of Copper Ions and Redox State in Regulation of Butane Monooxygenase in Pseudomonas butanovora J. Bacteriol., April 15, 2008; 190(8): 2933 - 2938. [Abstract] [Full Text] [PDF] |
||||
![]() |
H. Arai, J. H. Roh, and S. Kaplan Transcriptome Dynamics during the Transition from Anaerobic Photosynthesis to Aerobic Respiration in Rhodobacter sphaeroides 2.4.1 J. Bacteriol., January 1, 2008; 190(1): 286 - 299. [Abstract] [Full Text] [PDF] |
||||
![]() |
A. K. Schmid, D. J. Reiss, A. Kaur, M. Pan, N. King, P. T. Van, L. Hohmann, D. B. Martin, and N. S. Baliga The anatomy of microbial cell state transitions in response to oxygen Genome Res., October 1, 2007; 17(10): 1399 - 1413. [Abstract] [Full Text] [PDF] |
||||
![]() |
M. C. Sanchez-Sutil, N. Gomez-Santos, A. Moraleda-Munoz, L. O. Martins, J. Perez, and J. Munoz-Dorado Differential Expression of the Three Multicopper Oxidases from Myxococcus xanthus J. Bacteriol., July 1, 2007; 189(13): 4887 - 4898. [Abstract] [Full Text] [PDF] |
||||
| ||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
| HOME | HELP | FEEDBACK | SUBSCRIPTIONS | ARCHIVE | SEARCH | TABLE OF CONTENTS |
| All ASBMB Journals | Molecular and Cellular Proteomics |
| Journal of Lipid Research | ASBMB Today |