Functional Interactions between the Carbon and Iron Utilization Regulators, Crp and Fur, in Escherichia coli

ABSTRACT In Escherichia coli, the ferric uptake regulator (Fur) controls expression of the iron regulon in response to iron availability while the cyclic AMP receptor protein (Crp) regulates expression of the carbon regulon in response to carbon availability. We here identify genes subject to significant changes in expression level in response to the loss of both Fur and Crp. Many iron transport genes and several carbon metabolic genes are subject to dual control, being repressed by the loss of Crp and activated by the loss of Fur. However, the sodB gene, encoding superoxide dismutase, and the aceBAK operon, encoding the glyoxalate shunt enzymes, show the opposite responses, being activated by the loss of Crp and repressed by the loss of Fur. Several other genes including the sdhA-D, sucA-D, and fumA genes, encoding key constituents of the Krebs cycle, proved to be repressed by the loss of both transcription factors. Finally, the loss of both Crp and Fur activated a heterogeneous group of genes under σS control encoding, for example, the cyclopropane fatty acid synthase, Cfa, the glycogen synthesis protein, GlgS, the 30S ribosomal protein, S22, and the mechanosensitive channel protein, YggB. Many genes appeared to be regulated by the two transcription factors in an apparently additive fashion, but apparent positive or negative cooperativity characterized several putative Crp/Fur interactions. Relevant published data were evaluated, putative Crp and Fur binding sites were identified, and representative results were confirmed by real-time PCR. Molecular explanations for some, but not all, of these effects are provided.

Processes that alter the rates of transcriptional initiation, elongation and termination, mRNA degradation, and mRNA translation control gene expression in all living organisms. Of these, regulation of transcriptional initiation is of greatest physiological importance in bacteria (34). About one-fourth of all Escherichia coli proteins interact functionally with nucleic acids, and about one-fourth of these are pleiotropic transcriptional regulators that control expression of multioperon regulons, primarily by influencing the transcriptional initiation step (57,72). In fact, each regulon is defined by the activity of a transcription factor that coordinately controls expression of the operons included within that regulon (28). The global regulator, responsive to select cellular stimuli, binds to specific sites in the control regions of these operons. Examples of such pleiotropic transcription factors in E. coli include Crp, a primary sensor of carbon availability (8,15,23,65), NtrBC, a sensor of nitrogen availability (32,63,77), CysB, the sensor of sulfur availability (38,39,41,42), and Fur, a dominant sensor of iron availability (30,31,50,71). Crp when complexed with cyclic AMP (cAMP) binds to a consensus sequence in the DNA and activates transcription of many genes while repressing transcription of a few others. Some genes are subject to both positive and negative regulation by the cAMP-Crp complex due to the presence of multiple promoters (36,64). Fur generally represses iron transport and iron siderophore biosynthetic genes when complexed with ferrous iron. Under iron limiting conditions, iron dissociates from Fur, and increased transcription ensues (6,30,31).
Functional interactions between these regulons are believed to coordinate the activities of the different metabolons so that the supply of one type of nutrient matches the supply of other essential types of nutrients (28). Thus, for example, multiple links between carbon and nitrogen metabolism have been identified (43,59). Moreover, data providing evidence for functional links between carbon and sulfur metabolism (62) as well as carbon and iron metabolism (18) have been presented. In the case of the iron-carbon connection, Crp, the primary transcriptional regulator of carbon metabolism, modulates transcription of the fur gene which encodes the ferric uptake regulator, Fur (18). These and other potential global interrelationships serve as the bacterial "nervous system," coordinating the various activities of the cell (28,44,66).
In this paper we use a whole-transcriptome approach to provide evidence for previously unidentified potential interactive mechanisms involving Fur and Crp. We show that many iron metabolic genes are subject to dual Crp/Fur control although our transcriptome results cannot distinguish between direct and indirect effects of the loss of these transcription factors. For many iron transport genes (class A genes), the consequence of the loss of Crp (transcriptional down regulation) is opposite to that of the loss of Fur (transcriptional up regulation). We also show that the loss of Crp and Fur regulates expression of the genes encoding the oxidative stress enzyme, superoxide dismutase, and the aceBAK operon encoding the glyoxalate shunt enzymes (class B genes) in directions opposite to those observed for the iron transport genes. These genes are directly or indirectly repressed by the loss of Fur and activated by the loss of Crp in a fur background. Moreover, other genes were identified that are either repressed (class C gene) or activated (class D genes) by the loss of both transcription factors. For example, several enzyme complexes of the Krebs cycle are directly or indirectly repressed by the loss of both Fur and Crp. The transcriptome results were confirmed in representative cases by using the real-time PCR. Putative factor binding sites were identified when present, and based in part on published data, mechanisms of transcriptional regulation could sometimes be proposed. Physiological relevance was also considered. This work serves to suggest a previously unrecognized transcriptional regulatory connection between iron and carbon metabolism. present in the RNA preparations was removed by RNase-free DNase (Stratagene). cDNAs were synthesized using a superscript first-strand synthesis kit (Invitrogen) according to the manufacturer's instructions and stored at Ϫ20°C prior to use. polA, encoding DNA polymerase I, was included as an internal control (50). Real-time PCR was carried out on a LightCycler instrument (Roche Diagnosis Corporation) following the manufacturer's recommended protocol. Primers used for the real-time PCR were as shown in the supplemental Table S1 on our website. Data analysis. Raw expression intensity values were determined from array images analyzed using the Affymetrix Microarray Suite 5.0 software. These data were then scaled, normalized to a target value of 1,000, and imported into a Microsoft Excel spreadsheet for further analysis. The intensity values in the data sets ranged from about 3 ϫ 10 Ϫ1 to 8 ϫ 10 4 . Three data sets were obtained for each of eight experimental condition: (i) wild-type (wt) LB, (ii) wt LB plus glucose, (iii) crp LB, (iv) crp LB plus glucose, (v) fur LB, (vi) fur LB plus glucose, (vii) crp fur LB, and (viii) crp fur LB plus glucose. Pairwise comparisons were performed with the three data sets for each strain and condition to determine the Pearson correlation coefficient using the following formula: For each triplicate data set, the two sets with the highest Pearson correlation coefficient were retained for further analysis. The statistical significance of n-fold changes between data sets from different strains and conditions was determined as described by Caldwell et al. (9). Briefly, for each pair of replicate data sets, a table was generated from the average of the logs of the two replicate signals (ALS) and their log ratio (LR). From this table, using a series of sliding windows with a size of 201 ALS values, the mean ALS values and the standard deviation of the corresponding LR were calculated (SD LR ).
Phenotypic studies. Phenotypic analyses were performed using 96-well Biolog microtiter plates including GN2 and Phenotypic Microarray plates PM1 and PM2 (5). These Biolog plates are designed to measure oxidation of various carbon sources by gram-negative bacteria.

RESULTS
Identification of genes whose expression is influenced by the loss of both Fur and Crp. Analysis of transcriptome data allowed the identification of genes that showed significant changes in expression level in response to the loss of both Crp and Fur. Most of the genes presented in Table 1 met the following criteria. (i) They showed statistically significant changes in gene expression upon introduction of both the crp and fur mutations. The n-fold change was at least 4 times higher than the SD LR for the ALS of the gene when comparing wt and fur or wt and crp fur (see reference 23 for explanation of the statistical analyses of the microarray data). (ii) A large signal was observed in at least one of the mutant strains (signal intensity value, Ն1,000). (iii) The ratio for crp fur/fur was Յ0.5 or Ն2 following growth in LB broth. Table 1 also indicates whether or not putative Crp and/or Fur binding sites in the control regions of these operons could be found using the GRASP-DNA program (68). The putative binding site sequences and their locations can be found in supplemental Table S2A on our website (http://www.biology.ucsd.edu/ϳmsaier/supmat/crp_fur /index.html). Published regulatory data for these operons are summarized in Table S2B in the supplemental material. It is important to note that some of the genes in Table 1 are in operons, and not all of the genes in a specific operon always met these criteria. However, in the discussion below, all of the cistrons of each operon are usually considered.
As shown in Table 1, genes meeting these criteria can be divided into four categories: A, genes that appear to be repressed by Fur but activated by Crp in a fur background (i.e., genes involved in iron sequestration and transport); B, genes that appear to be activated by Fur but repressed by Crp in a fur background (i.e., the sodB gene and the aceBAK operon); C, a Genes tabulated gave high reproducibility for replicate values, large signals for at least one of the strains assayed, fur/wt ratios greater than 2 or less than 0.6, and crp fur/fur ratios of greater than 2 or less than 0.6 for most of the genes in an operon. Brackets to the left of the gene designations indicate genes included within a single operon. Note that transcriptome data do not by themselves allow direct and indirect effects of the loss of a transcription factor to be distinguished. Category A, genes apparently repressed by Fur and activated by Crp in a fur background; category B, genes apparently activated by Fur and repressed by Crp in a fur background; category C, genes apparently activated by Fur and activated by Crp in a fur background; category D, genes apparently repressed by both Crp and Fur. b Values are the average signal values provided as described in Materials and Methods. The abbreviations are as follows: wt, wild type; crp, the isogenic crp null mutant; fur, the isogenic fur null mutant; crp fur, the isogenic double mutant; G, grown in LB medium plus glucose. When G is not indicated, cells were grown in LB medium.
c The presence of putative binding sites for Fur and Crp is indicated by an X. The predicted binding sites are present within or close to the regulated gene in most cases (see our website). Putative sigma-70 promoters (34) were identified for all tabulated genes. The putative Crp and Fur binding sites and their positions are presented on our website (Table S2A). The occurrence of other established binding sites (72) within these control regions is also reported on our website (Table S2B).
d Two putative Crp binding sites are present upstream of pdhR.
genes that appear to be activated by Fur and activated by Crp in a fur background (i.e., the cspA cold shock regulatory gene and the ompF outer membrane porin gene as well as several genes encoding Krebs cycle enzymes); and D, genes that appear to be repressed by both Crp and Fur (e.g., the glgS glycogen biosynthetic gene and the rpsV ribosomal protein gene).
It should be noted that the transcriptome data alone do not distinguish direct from indirect effects of the loss of a transcription factor. Genes apparently repressed by Fur and activated by Crp in a fur background. Fur is known to repress many iron transport system-encoding genes under conditions of iron sufficiency (50). For representative genes of this group, the detected transcript levels increased 1.8-to 8.8-fold as a result of fur inactivation when the wild-type and mutant strains were grown in LB medium (Table 1). Simultaneous inactivation of crp and fur caused a substantial decrease in transcript level for all genes in this group relative to the fur single mutant. The crp fur/fur ratio for these genes ranged from 0.2 to 0.7 (Table 1). Putative Crp and/or Fur binding sites could usually be identified in these operons, suggesting that many of the mRNA level changes observed were direct consequences of factor binding.
The 21 genes that comprise category A are included within several different transcription units ( Table 1). Some of these will be discussed. Genes aceE and aceF encode pyruvate dehydrogenase (E1) and lipoate acetyltransferase (E2), respectively. These two proteins and lipoamide dehydrogenase (E3), encoded by gene lpd, constitute the pyruvate dehydrogenase (PDH) complex. Synthesis of the PDH complex is induced by pyruvate, repressed during anaerobic growth, and partially repressed by growth in media containing excess glucose, citric acid intermediates, or acetate (20). Genes encoding PDH subunits are part of the complex operon pdhR-aceEF-lpd where the pdhR gene encodes a transcriptional repressor. Two different transcripts originate from this group of genes, one including pdhR-aceEF-lpd and another only lpd. The larger transcript originates from a single promoter in front of the pdhR gene. Weak activation of the pdhR promoter by Crp and Fnr and corresponding putative binding sites upstream of the Ϫ35 promoter sequence have been reported (61).
The ackA and pta genes form an operon and encode acetate kinase and phosphotransacetylase, respectively. These enzymes constitute a pathway for the generation of acetate and ATP from acetyl-coenzyme A (CoA). The ackA-pta operon is part of a group of genes regulated by the two-component regulator CreCB. Expression of the ackA-pta operon is induced during growth in minimal medium containing glucose (2).
The fepA gene, encoding the outer membrane receptor for enterobactin, forms an operon with entD, encoding an enzyme involved in enterobactin biosynthesis. The genes fes, entF, and fepE form a second operon and encode two more biosynthetic enzymes and a ferric enterobactin transport protein, respectively. The fepA and fes operons, divergently transcribed from overlapping promoters, are controlled by the binding of Fur to a single binding site (35).
The entCEBA-ybdB operon, encoding enterobactin biosynthetic enzymes, is transcribed from a bidirectional promoter region shared with fepB that encodes a protein involved in enterobactin transport. Transcriptome data for fepB showed low levels of transcript expression. Therefore, it did not meet our selection criteria. Divergent promoters for the entC and fepB operons are separated by 31 bp, and for each, a functional Fur binding sequence has been identified (7).
Genes ybiL, ybiX, and ybiI are clustered and possibly form an operon. The gene ybiL, also known as fiu (ferric iron uptake), encodes an outer membrane protein having a putative TonB box (29,51,52,56). The ybiX gene encodes a protein similar to an iron-regulated hydroxylase-encoding gene from Pseudomonas aeruginosa (piuC), and ybiI encodes a probable C 4 -Zn finger protein (50,54). Sequence and functional analyses of the promoter region of fiu have revealed the presence of four overlapping Fur binding sites (52).
The gene yncE potentially encodes a protein with sequence similarity to a pyrolo-quinoline quinone containing periplasmic oxidase. A putative Fur binding site has been found in the promoter region of this gene (50).
Iron-citrate import depends on the genes of the fecABCDE operon. fecA encodes an outer membrane iron-citrate receptor, while fecBCDE encodes components of an ABC ironcitrate transporter. A separate operon, fecIR, encodes the alternative sigma factor, FecI, and an anti-sigma factor, FecR. When iron concentrations are low and iron citrate is present, the receptor protein FecA transduces a signal to the transmembrane FecR protein, which causes FecI to become active and promote transcription from the fecA promoter (6). Functional Fur binding sites have been detected in the promoter regions of both the fecABCDE and the fecIR operons (1). According to our transcriptome data, the fecI gene showed a twofold increase in expression in the fur mutant relative to the wild type, and deletion of the crp gene caused a slight decrease in fecI expression in both the wild-type and fur mutant backgrounds (data not shown). fecR gave a low signal. Real-time PCR data confirmed that the fur mutation caused a substantial increase (5ϫ) in fecI expression (data not shown). Our results are in general agreement with those of Braun et al. (6). The effects of the fur mutation on fecI gene expression could partially explain the data reported in Table 1 for the fec structural genes. Other genes in category A include genes encoding an additional outer membrane receptor for iron transport, enzymes involved in carbon metabolism, and two proteins of unknown function.
The cirA gene encodes the colicin I receptor. This protein may bind a catechol-type siderophore (14,29,56). cirA is transcribed from two overlapping promoters. Footprinting analysis revealed the presence of a Fur binding site overlapping both promoters. Transcription from cirA is significantly reduced in cya and crp mutants, and a potential Crp binding site has been identified downstream of the tandem promoters (25).
Genes apparently activated by Fur and repressed by Crp in a fur background. The second group of genes showing significant changes in expression level when crp was inactivated in a fur background includes the sodB and the aceBA genes ( Table  1). The sodB gene encodes iron-dependent superoxide dismutase (3) while the aceB and aceA genes encode malate synthase and isocitrate lyase, respectively, the enzymes of the glyoxalate shunt (12). Fur has been shown to activate the sodB gene, at least in part by promoting mRNA stabilization (21). Our results showed that in the fur mutant, the sodB transcript level decreased fivefold compared to the wild-type strain. In VOL. 187, 2005 Crp/Fur INTERACTIONS 983 the crp fur double mutant, the sodB transcript level increased 4.7-fold over that observed in the fur mutant although the crp mutation in the wild-type background had no effect on gene expression (Table 1). It is surprising that in the double mutant, sodB expression is repressed by glucose, although this is not true in the wild-type or single mutant genetic backgrounds. Regulation of sodB expression has been shown to be dependent on the small RNA encoded by the ryhB gene (45). RyhB RNA exerts posttranscriptional down regulation of sodB and other genes involved in iron metabolism, apparently by an antisense RNA mechanism (45). ryhB expression is repressed by Fur (73). If RyhB decreases the stability or causes mRNA truncation of the sodB mRNA, its effect could in part explain the transcriptional results reported in Table 1.
Other genes known to be regulated by RyhB (acnB, fumA, bfr, fln, and sdhC) also show a general trend of a decrease in transcript level in response to this small RNA. The bfr (bacterioferrin) gene shows increased expression in the crp mutant, slightly decreased expression in the fur mutant, but dramatically increased expression in the crp fur double mutant (data not shown). Particularly in this last respect, the bfr gene behaves the same as the sodB gene. Also, the acnA gene shows similar behavior although the changes in gene expression levels are substantially less dramatic (data not shown).
The aceBAK operon encodes the glyoxalate shunt enzymes required for growth of the bacterium on acetate as a carbon source (11,12). Regulation of this operon is extremely complex, as suggested by the transcription factor binding sites presented in Fig. 1. Putative overlapping Crp and Fur binding sites were identified upstream of the aceB gene, and these also overlap with two IclR binding sites ( Fig. 1; see Discussion). Expression of this operon is subject to glucose repression which is largely abolished by either the crp or the fur mutation, either of which also reduces the expression level to a low value. The double mutant shows greatly enhanced expression, comparable to that of the wild-type strain grown in the absence of glucose, and this expression is no longer subject to strong glucose repression (Table 1).
Genes apparently activated by Fur and activated by Crp in a fur genetic background. The cspA gene encodes a cold shock transcriptional regulator (22,27). Expression of this gene appeared to be activated by glucose, but also by Fur, and was further activated by Crp (Table 1). cspA showed depressed expression when either the crp or the fur mutation was introduced, and the double crp fur mutant showed a more than additive effect. Glucose activation was not abolished by either the crp mutation or the fur mutation although expression in the double mutant was largely glucose independent.
Another gene, ompF, encoding an outer membrane porin, showed apparent repression by Crp in the wild-type genetic background and activation by Fur in the same background (Table 1). However, while the crp mutation enhanced ompF gene expression in the wild-type genetic background, it reduced expression in the fur background. ompF expression is known to be subject to complex regulation by multiple factors, but regulation by Fur had not been reported previously (16,17,40,48,49,53,60). Our results suggest that Fur is required for Crp-dependent repression and that in the absence of Fur, Crp activates gene expression.
The four genes of the sdh operon (sdhCDAB), encoding the four subunits of iron-dependent succinate dehydrogenase, showed strong glucose repression, mild apparent Fur activation in a wild-type background, and strong apparent activation by Crp both in the wild-type background and in the fur genetic background. Expression of the sdh operon is known to be subject to negative control by the RyhB antisense RNA that promotes mRNA truncation at a region of complementarity in the first gene of the operon (sdhC). RyhB requires the RNA binding protein, Hfq, for this activity. It seems that RyhB provides a general mechanism for the down regulation of nonessential iron-containing genes when iron is limiting (45).
The sucABCD operon that clusters together with the sdh operon showed similar regulation in response to the loss of Crp and Fur. Genes sdhCDAB-sucABCD form a complex operon. Three promoters have been identified in this cluster, one in front of the sdhC gene, a second within sdhC, and a third in front of the sucA gene. Transcription starting from the promoter in sdhC usually terminates in the sdhB-sucA intergenic region. However, a large 10-kb transcript including sdhCDAB-sucABCD has been detected. The transcript starting in front of sucA terminates at the 3Ј end of sucD (13,75). Expression from the sdhC promoter is depressed by anaerobiosis, inclusion of glucose in the growth medium, and inactivation of either the crp or the cya gene. Binding sites for ArcA, Crp, and Fnr have been identified in this promoter region (55). The promoter upstream of the sucA gene is repressed by integration host factor (IHF) and partially activated by S . It is not repressed by glucose or anaerobiosis. Binding sites for ArcA, Fnr, and IHF have been identified in this promoter region. It is believed that the strong glucose repression (Table 1) is due to read-through from the sdh operon (13).
As for the cspA and ompF genes, Crp appeared to activate transcription of the sdh and suc genes although the activation by Fur was less pronounced. As noted above, the sdh operon is subject to regulation by RyhB, and the same could be true of the suc operon. The consequences of the loss of both Crp and Fur appear to be less than the additive effects of the loss of either one alone for both of these operons. Finally, the fumA gene, encoding fumarase, shows strong apparent activation by both Crp and Fur, but the loss of both of these transcription factors together, as for the sdh and suc operons, gave a much less than additive effect. fumA gene expression, like sdh operon expression, is negatively regulated by RyhB (45).
Genes apparently repressed by both Crp and Fur. Nine genes appeared to be strongly repressed by both Crp and Fur, and in all such cases loss of both Crp and Fur resulted in greater than additive enhancement of gene expression ( Table  1). Three of these genes encode proteins of unknown function, but the cfa gene encodes a cyclopropane fatty acid synthase (CFA synthase); the glgS gene encodes a glycogen biosynthetic protein, and rpsV encodes a 30S ribosomal protein, S22 or SRA (see below). The osmE gene, an osmotically induced gene, encodes an envelope protein of unknown function; the tktB gene encodes transketolase II, and the yggB gene encodes a mechanosensitive channel of the MscS family (58). Several of these genes are subject to S control as noted below.
During growth of E. coli, CFA synthase activity is detected throughout the growth cycle, but the highest activity is observed upon entering the stationary phase. The cfa gene is transcribed from two promoters, one requiring 70 and the other requiring S . No binding sites for transcriptional regulators have been reported for this gene. CFA synthase activity is unstable in vivo, explaining a decrease in activity when cells progress into stationary phase (74).
The glgS gene encodes a 7.9-kDA protein that stimulates glycogen synthesis. However, the biochemical function of GlgS has not been elucidated. Two promoters have been identified and transcription from both promoters requires functional cya and rpoS genes (33).
Gene rpsV (sra) encodes the ribosome-associated protein D, the abundance of which increases upon entry into stationary phase. For this reason, protein D has been renamed the SRA (stationary phase-induced ribosome-associated) protein. This gene possesses a single promoter that is partly dependent on S . Transcription from this promoter is decreased in mutants lacking S , FIS (factor for inversion stimulation), HNS (histone-like nucleoid structuring protein), cAMP, and ppGpp (37).
The osmE gene displays biphasic expression when cells are grown at elevated osmotic pressure. Induction is detected during exponential growth, but stronger induction occurs at the onset of stationary phase. Transcription of osmE starts at a single promoter. DNA supercoiling is involved in osmE expression during exponential growth in media of high osmolarity, but stationary phase induction depends on S (10).
All nine genes respond similarly to the loss of either Crp or Fur, showing about twofold enhancement of gene expression in either case. The simultaneous loss of both proteins showed synergistic enhancement, with increased expression, relative to the wild-type strain, of about 10-fold. These results suggest that both Crp and Fur somehow antagonize the repressive effects of the other transcription factor. While all but one of these genes exhibit putative Crp binding sites in their control regions, none could be shown to have a recognizable Fur binding site, suggesting an indirect mechanism, possibly involving interaction of Fur with the Crp binding site (see Discussion and Table 4).
Confirmation of transcriptome results using real-time PCR. Ten genes from various categories of Table 1 were selected for confirmation using real-time PCR. Table 2 shows the results of these analyses. The five genes analyzed from category A are fepA, entC, ybiL, fecA, and cirA. When expression levels are measured by real-time PCR, a much more dramatic increase in transcript level is observed for these genes when fur is inactivated compared with the levels measured with the microarrays. In the crp fur mutant, the transcript levels of these genes are dramatically decreased compared to those of the fur mutant, regardless of the method used to measure relative expression level. Nevertheless, relative values are in good agreement. In a similar fashion, expression levels for sodB from category B, cspA and ompF from category C, and glgS and rpsV from category D were analyzed by real-time PCR. Qualitative agree- ment when using the two methods for measuring relative expression levels was always obtained ( Table 2). The results obtained by real-time PCR are thus in good agreement with the data obtained from the transcriptome analysis presented in Table 1 although a systematic quantitative difference was consistently observed. Phenotypic studies. Table 3 summarizes the results when substrate oxidation was measured for wild-type, crp, fur, and crp fur mutant strains. In general, the fur mutation by itself showed no effect on the apparent rates of oxidation although the crp mutation alone frequently had dramatic effects. The double mutant usually showed depressed oxidation relative to the crp mutant although the opposite was sometimes observed (Table 3). For several substrates, the crp mutant strain showed positive responses while the crp fur double mutant strain showed significantly reduced responses. However, for D-mannitol and inosine the crp mutant gave significantly reduced oxidation responses while the double mutant gave more-positive responses. Since the wild-type and fur mutant strains showed fully positive responses, the results suggest that Fur and Crp may interact to regulate cellular carbon metabolic activities in previously unrecognized ways. For a more complete description of the phenotypic analyses, see Table S3 in the supplemental material.
The phenotypic data cannot be compared directly with the transcriptome data because the substrates listed in Table 3 are not subject to metabolism by enzymes and transporters encoded by genes listed in Table 1 and the conditions are not comparable. However, there is agreement in the following respects: (i) the crp mutation often depresses the rate of oxidation as expected; (ii) the fur mutation sometimes has no effect by itself, but usually modifies the phenotype of the crp mutant; and (iii) for a majority of substrates, a cooperative effect of the two mutations is observed, but (iv) for a minority of the substrates (D-mannitol and inosine), apparently antagonistic effects are observed. It seems that Fur and Crp can somehow act either cooperatively or antagonistically to produce a phenotype. Whether these effects represent direct or indirect consequences of the loss of Crp and Fur cannot be ascertained.
Effect of glucose on growth of E. coli strains. The four strains used in these studies were grown in LB media Ϯ 0.4% glucose in shake flasks under exactly the same conditions used for the microarray and real-time PCR experiments for quantitative measurement of growth rates. The growth curves are shown in Fig. 2. In the absence of glucose ( Fig. 2A), the fur mutant strain   grew almost as well as the isogenic wild-type strain. The crp mutant grew decidedly more slowly, and the crp fur double mutant grew the slowest. Addition of glucose (Fig. 2B) increased the growth rates of all strains. Thus, the wild-type and fur mutants grew equally well in the presence of glucose, significantly faster than observed for the LB grown cells. The crp mutant and the crp fur double mutant also showed substantial increases in growth rate when glucose was present. These relative growth rates were confirmed by examining colony sizes after growth on LB plus glucose agar media and by measuring growth rates in minimal liquid media (data not shown).

DISCUSSION
The comparative analyses of transcriptome data from wildtype, crp, fur, and crp fur mutant strains allowed the identification of genes that are subject to regulation resulting from the loss of both Fur and Crp. Many of these genes were known to belong to the Fur regulon while others were known to belong to the Crp regulon. Some of these genes, however, had not previously been known to belong to either regulon, and very few of them were known to belong to both.
Since all the genes examined in this report showed changes in transcript levels in the crp fur double mutant relative to the fur single mutant, it was possible that they would have Crp binding sites in or near their promoters. A search for putative Crp binding sites in the genes shown in Table 1 was performed using the GRASP-DNA program (68). This search yielded significant matches to the reported Crp binding consensus sequence in or near many of these genes (Table 1). Detailed descriptions of Crp and Fur binding sites using the GRASP-DNA program (including their positions, sequences, and scores) are reported in Table S2A and published regulatory information about these genes is summarized in Table S2B in the supplemental material. Some of these data may be relevant to the control of gene expression in response to the availability of Fur and Crp. In a few cases (e.g., the cirA gene), direct binding of both Crp and Fur to the upstream control region has been reported (25,26).
Expression levels of all of the genes listed in Table 1 are significantly altered in response to the lack of Fur and Crp. The major subset of known Fur regulon genes displaying negative regulation by Fur but positive apparent regulation by Crp encode proteins related to iron transport and chelation ( Table  1). Most of these genes are repressed by Fur in the presence of excess iron, and therefore showed increased expression in the fur mutant under the growth conditions used. These Fur-repressed, Crp-activated genes are grouped into seven transcriptional units that have been shown to be regulated directly by Fur binding to their promoter regions (fepA-entD, fes-entF-fepE, entCEBA-ybdB, ybiLXI, yncE, fecABCDE, and cirA). Expression of the aceEF operon showed strong glucose activation that was abolished by both the crp mutation and the fur mutation; while both mutations enhanced operon expression about 2.5-fold, the double mutant showed about threefolddepressed expression relative to the fur or crp single mutant. The entC operon showed mild glucose repression, mild repression by the loss of Crp, about threefold activation by the loss of Fur, and restoration of the wild-type gene expression level in the crp fur double mutant with complete loss of glucose repres-sion. The results obtained by real-time PCR were qualitatively similar, but the regulatory effects were generally larger in magnitude. Thus, we observed good qualitative agreement although there were quantitative differences. It should be pointed out that transcriptional initiation of the fecABCDE operon is dependent on both Fur clearance from the promoter region and the subsequent binding of RNA polymerase containing the FecI sigma subunit (70). Since fecI expression proved to be regulated by Fur and Crp, some of the effects reported for the fec structural genes could be indirect.
The aceBAK operon and the sodB gene (Table 1) appeared to be activated by Fur and repressed by Crp in a fur background. The aceBAK operon (76) is subject to strong glucose repression, and this repression is abolished by either the crp or the fur mutation. The crp fur double mutant showed about sixfold-enhanced expression over either of the two single mutants. Overlapping putative Crp and Fur binding sites found upstream of the aceB gene ( Fig. 1 and supplemental Table  S2A) suggest direct control of this operon by both global regulators. The fact that activation of the aceBAK operon depends on the presence of both Crp and Fur (for glucose-repressible expression) or depends on the absence of both transcription factors (for glucose-independent expression) additionally suggests that protein-protein interactions play a role. However, the aceBAK operon has previously been reported to be regulated in a complex fashion by four transcription factors: IclR, IHF, Cra, and FadR ( Fig. 1) (12). The operon-specific IclR repressor binds to three sites in the aceB promoter region. Two of these sites overlap each other and also overlap the region to which Crp and Fur may bind (Fig. 1). IclR also autoregulates its own expression (76). It is possible, therefore, that Crp and Fur regulate iclR expression and/or that IclR interacts with Crp and/or Fur as well as with the C-terminal domain of the RNA polymerase alpha-subunit (76). Glucose-independent highlevel expression in the absence of Crp and Fur may occur if IclR repression is dependent on an interaction with either Crp or Fur or both. Additionally, physical interactions between Crp and/or Fur with other regulators such as IHF, Cra, and FadR could be important.
Expression of the sodB gene (Table 1) has been shown to be indirectly regulated by Fur via potential antisense interactions involving a small RNA encoded by the ryhB gene (21,45,46,47). This regulatory effect of RyhB, however, does not appear to fully explain the results reported here. Since no good Crp or Fur binding sites could be identified upstream or within the sodB gene, sodB regulation is likely to be indirect and prove complex.
A few operons listed in Table 1 (such as sdhC and fumA) that appear to be subject to dual control by Crp and Fur are also known to be regulated by RyhB (24,45). However, to what extent this small RNA or other small RNAs that influence mRNA stability account for the results reported here has yet to be determined.
In the case of genes included in Table 1, at least one Crp binding site per operon but no Fur sites could be identified. The mechanism of regulation of these genes by Fur may be indirect and may prove to be complex, possibly involving S . A careful analysis of several known Crp and Fur binding sites revealed conserved nucleotides common to both of the binding sequences (Table 4). Therefore, it is possible that Crp binding VOL. 187, 2005 Crp/Fur INTERACTIONS 987 sites are recognized by Fur and vice versa. Since the various Crp and Fur binding sites differ in sequence, it is probable that cross-binding specificity will be site specific. A regulatory interaction between the Fur and Crp regulons had been reported previously. Transcriptional levels of fur gene expression were reported to be repressed by Fur and stimulated by Crp (18). Both Fur repression and Crp stimulation of transcription caused modest changes in fur gene expression during exponential growth (Յ2-fold). In our studies, a total lack of the Fur protein in the fur mutant eliminated the possibility that the consequences of the loss of Crp in this genetic background could be explained by regulatory effects of Crp loss on fur gene expression. Thus, the effects of the crp mutation on the expression of other target genes in a fur genetic background must be independent of Fur. The observations of De Lorenzo et al. (18) are evidently not relevant to most of the results reported here. It is possible that the observed effects depend on direct binding of both factors to the DNA, to the binding of one factor to the other, or to regulation of elements higher in the hierarchy of transcriptional regulation.
Examination of transcriptome data for genes encoding proteins involved in transcriptional initiation revealed that expression of the rpoD gene, encoding the sigma-70 subunit of RNA polymerase, was significantly lower in the crp (fivefold) and crp fur mutants (twofold) when compared to wild-type and fur transcript levels, respectively (data not shown), and these results were confirmed using real-time PCR (Z. Zhang and M. H. Saier, Jr., unpublished results). It was therefore possible that the reduced transcript levels observed for Fur-repressed genes in the crp mutant were the result of decreased availability of sigma-70. Our preliminary experiments, however, suggest that this was not the case for class A genes (Table 1). Overexpression of the rpoD gene did not enhance expression of any of the genes listed in Table 2 (Zhang and Saier, unpublished). In the case of the sodB gene, the observed increase in transcript level in the crp fur double mutant when compared to the fur single mutant might be explained assuming that sigma-70 limitation causes a lower expression level of the small RNA RyhB which destabilizes the sodB mRNA. This possibility has yet to be examined.
As noted above, many genes subject to simultaneous Fur and Crp regulation are apparently repressed by Fur and activated by Crp in a fur mutant. However, the three other possible consequences of the loss of these two transcription factors in a wild-type background indicate that gene expression can be activated and repressed, activated and activated, or repressed and repressed by Fur and Crp, respectively. Although we confirmed the transcriptome results by real-time PCR for representative genes, the detailed regulatory mechanisms have been examined in very few of the genes listed in Table 1.
Crp regulation of Fur-controlled genes and Fur regulation of Crp-controlled genes could allow integration of signals for iron and carbon sufficiency. Because iron is required for some carbon metabolic processes, particularly aerobic Krebs cycle intermediary metabolism and electron flow, but not others, such as anaerobic sugar metabolism, the signals coordinating iron and carbon metabolism are likely to be complex. This report provides a glimpse into the nature and consequences of these interactions.