Multiple High-Throughput Analyses Monitor the Response of E. coli to Perturbations

See allHide authors and affiliations

Science  27 Apr 2007:
Vol. 316, Issue 5824, pp. 593-597
DOI: 10.1126/science.1132067


Analysis of cellular components at multiple levels of biological information can provide valuable functional insights. We performed multiple high-throughput measurements to study the response of Escherichia coli cells to genetic and environmental perturbations. Analysis of metabolic enzyme gene disruptants revealed unexpectedly small changes in messenger RNA and proteins for most disruptants. Overall, metabolite levels were also stable, reflecting the rerouting of fluxes in the metabolic network. In contrast, E. coli actively regulated enzyme levels to maintain a stable metabolic state in response to changes in growth rate. E. coli thus seems to use complementary strategies that result in a metabolic network robust against perturbations.

The complex relationships between different layers of biological information that give rise to cellular functions can only be captured by combining global measurements across these different levels (18). Moreover, most current global methods such as DNA microarray and shotgun proteomic analyses aim to interrogate the composition of intracellular components rather than to accurately determine the quantities of molecules in cells. More quantitative methods such as quantitative reverse transcription polymerase chain reaction (qRT-PCR) and metabolic flux analysis tend to be more targeted but are essential to detect small changes in response to perturbation. In the present study, we measured the global response of Escherichia coli K-12 cells to genetic and environmental perturbations at the level of protein and gene expression, and compared this with a detailed quantitative analysis of specific metabolic pathways.

To probe the global cellular response, we determined the relative quantities of gene and protein expression using DNA microarray analysis and two-dimensional differential gel electrophoresis (2D-DIGE) (9). For more detailed analysis, in addition to qRT-PCR and metabolic flux analysis, we used capillary electrophoresis time-of-flight mass spectrometry (CE-TOFMS) for metabolome analysis (10, 11) and liquid chromatography tandem mass spectrometry (LC-MS/MS) for absolute quantification of proteins (12). As molecular targets for these more detailed quantitative analyses, we selected genes, proteins, metabolites, and fluxes that make up the bulk of reactions of the central carbon metabolism of E. coli K-12, including glycolysis, the pentose phosphate pathway, and the tricarboxylic acid cycle (fig. S1). Because these pathways play a crucial role in energy generation and in the production of important macromolecular precursors, we wanted to characterize the cellular response to disruptions in order to fulfill basic metabolic needs and maintain the growth phenotype.

We collected data using these multiple high-throughput analyses for wild-type E. coli K-12 and 24 single-gene disruptants, which were selected from the Keio collection (13) and cover most viable glycolysis and pentose phosphate pathway disruptants (i.e., galM, glk, pgm, pgi, pfkA, pfkB, fbp, fbaB, gapC, gpmA, gpmB, pykA, pykF, ppsA, zwf, pgl, gnd, rpe, rpiA, rpiB, tktA, tktB, talA, and talB). The cells were grown at a single fixed dilution rate of 0.2 hours–1 in glucose-limited chemostat cultures, and wild-type cells cultured at the same specific growth rate were used as a reference sample for comparison. To allow a comparison of the effects of these genetic perturbations with the effects of environmental perturbations, wild-type cells were examined at several different dilution rates (0.1, 0.2, 0.4, 0.5, and 0.7 hours–1). In chemostat cultures, the concentration of growth-limiting substrate can be controlled by the dilution rate (14); the dilution rate was thus varied in this study from an almost glucose-starved state to a nearly unlimited glucose supply. Changes in dilution rate therefore affect growth rate and can be considered a specific type of environmental change.

To evaluate differences between samples in the various perturbation experiments, we used an expression index (EI) (12), which scales the data to the same deviations for all molecules examined, regardless of their nature (SOM text), so that changes in mRNA, protein, and metabolite can be compared more easily. Figure 1 shows the EIs derived from the targeted analysis of metabolic pathway components. For mRNA and protein, the effects of changes in growth rate are obvious and appear to be larger than the effects of most gene disruptions. In contrast, except for nucleotides, changes in metabolites did not show any distinctive patterns, which could indicate different responses in metabolite levels against various perturbations. The lower levels of nucleotides observed in some disruptants (glk, pgm, pgi, pykA, and pykF) may reflect problems in quantification owing to the chemical instability or inconsistent extraction efficiency of these molecules, and we thus provide results calculated including and excluding nucleotide measurements.

Fig. 1.

Overall representation of changes in the targeted cellular components. The heat map shows the EI values of intracellular components that were detected in more than half the samples. RF, reference sample (wild-type cells cultured at a specific growth rate of 0.2 hours–1); GR, wild-type cells cultured at the indicated specific growth rates; KO, single-gene knockout mutants cultured at a specific growth rate of 0.2 hours–1. The numbers of individual components shown are metabolites, 130; proteins, 57; and mRNA transcripts, 85. The same type of representation, including all component names, is shown in fig. S3, and all EI values are shown in table S8.

To capture the magnitudes of the effects seen in Fig. 1, we calculated an average EI, the average expression index (AEI) (12), for each type of component (mRNA, protein, and metabolite) and compared them between samples (Fig. 2 and tables S10 and S11). One-way analysis of variance (one-way ANOVA) and post hoc tests (Games-Howell test) were performed (12) forall comparisons to test for statistically significant differences (P < 0.05) (SOM text).

Fig. 2.

AEI values for mRNA, protein, and metabolite levels. RF, GR, and KO are defined in the legend of Fig. 1. (A and B) The AEI values of a specific sample are represented as a heat map. The three panels in each box show three quantitative values corresponding to the AEI for mRNAs as determined by qRT-PCR (left panel), proteins as determined by LC-MS/MS (middle panel), and metabolites as determined by CE-TOFMS (right panel), respectively. (A) AEI values for GR. The labels in each row of the heat map correspond to either a specific growth rate (hours–1) or the two reference samples (RF02 and RF03). (B) AEI values for all gene disruptants displayed at their respective locations on the metabolic pathway map. Abbreviations for metabolites are as shown in table S1. (C) Quantitative measurements obtained by targeted analysis. (D) Semiquantitative cell-wide measurements. The numerical values are provided in table S9. In (C) and (D), numbers 1, 2, 3, and 4 correspond to specific growth rates of 0.1, 0.4, 0.5, and 0.7 hours–1, and numbers 5, 6, and 7 correspond to rpe, pgi, and pgm disruptants, respectively.

At higher growth rates, the AEI for mRNAs and proteins gradually increased in both targeted and global quantification methods, suggesting that E. coli actively regulates global gene and protein levels to meet increasing metabolic demands (Fig. 2, A, C, and D). In contrast, the AEI values for metabolites did not change significantly with growth rate. This relative stability in metabolite level may be a consequence of the observed regulation of enzyme levels. Closer examination of particular components revealed the induction and repression of mRNAs and specific enzymes after a change in specific growth rate (figs. S5 and S6). The levels of proteins involved in glucose uptake (PtsH and PtsI) were induced considerably when the specific growth rate increased. This is to be expected, because substrate uptake is the first step in fulfilling nutrient demands. The observed increase in the expression of phosphoenolpyruvate carboxylase (Ppc), decrease in the expression of phosphoenolpyruvate carboxykinase (PckA), and repression of glyoxylate shunt enzymes are also logical responses to enhance the tricarboxylic acid cycle flux (2), which corresponds to energy production under aerobic conditions and the components of which are important precursors in biomass production.

In contrast to the changes observed in wild-type cells cultured at various growth rates, most gene disruptants showed only subtle changes in AEIs for both mRNA and protein of central carbon metabolism enzymes, which fell within the range of variations observed in wild-type samples at the same specific growth rate (Fig. 2C). The AEI values for these targeted analyses of mRNAs and proteins in all disruptants were smaller than the AEI values observed for wild-type cells at a specific growth rate of 0.7 hours–1, and the difference was significant in 23 disruptants for mRNAs and in 16 disruptants for proteins. Similar results were obtained for the AEI values representing the global analysis of mRNA and protein expression (Fig. 2D). These findings suggest that E. coli does not appreciably respond to the loss of a single central carbon metabolism enzyme by regulating the abundance of other compensatory enzymes. This was especially apparent when examining changes in the expression of genes encoding isozymes of the respective disrupted enzymes (fig. S5). Isozymes that might be expected to compensate for the loss of an important enzyme were not up-regulated, although many were induced at higher growth rates. The two exceptions to this general pattern were the induction of pfkB in the pfkA disruptant and the very considerable induction of rpiB in the rpiA disruptant. However, in both cases, we found that cells accumulated additional mutations that are known to induce expression of the isozymes (15, 16). Therefore, the induction of these genes was not solely dependent on direct regulation (SOM text). It seems apparent that, except for these genes, enzymatic capacity is already expressed in the cell through structural redundancy (presence of isozymes and/or alternative routes).

The AEI values for metabolites in three disruptants (pgm, pgi, and rpe) were higher than those of wild-type cells at high specific growth rates (0.5 and 0.7 hours–1) (Fig. 2C). This phenomenon does not appear to be caused by the problems associated with nucleotide measurements in some disruptants, because similar results were obtained even when nucleotides were excluded from the AEI calculation (Fig. 2C). In the rpe disruptant, which had a particularly high AEI for metabolites, the levels of a few linked metabolites were increased to a large extent (Fig. 3A, rpe). However, these local changes were not responsible for the large metabolite AEI for this particular disruptant. Also, an overall stability in metabolite levels, i.e., low AEI, was maintained in most disruptants, even if the local metabolic network was perturbed (Fig. 2B and Fig. 3B, talB). This stability of metabolite levels is likely achieved through structural robustness of the metabolic network itself (17, 18), even when considerable rerouting of fluxes occurs. For example, in the zwf disruptant, the overall flux within the pentose phosphate pathway was reversed to compensate for the blocked route for the production of essential molecules, but, interestingly, this task can be achieved without large changes in enzyme levels in the pathways involved (Fig. 3C and table S5).

Fig. 3.

Perturbations in the local metabolic network caused by enzyme gene disruptions. The EI values for proteins (circled) and metabolites (boxed) in rpe (A), talB (B), and zwf (C) disruptants are shown mapped onto the local metabolic network. The disrupted enzymes are highlighted with an X, and the specific reaction eliminated by the disruption is shown with a dashed line. Reactions indicated with red arrows are running in the direction opposite to that in the wild-type. Colorless symbols represent components that were not measurable, and components not detected in the specific mutant are shown with an asterisk.

Our results reveal that the intracellular metabolic network of E. coli is remarkably stable and robust in the face of various types of disturbances. Highlighting the functional importance of this global stability for efficient cell growth, E. coli can actively respond to changes in the concentration of growth-limiting substrate by regulating the level of enzyme expression to maximize growth rate, which is reflected in the observed stability of metabolite levels. Using this strategy, E. coli can react quickly and effectively to this environmental change, but the strategy might be costly, because the cell must provide additional systems (such as sensor proteins, signal mediators, and transcriptional regulators) to detect and react appropriately to each specific perturbation. Because active response to the environment at the level of gene expressions is a widespread phenotype of E. coli, this might be a common strategy that the cell can use when it faces environmental changes. This strategy contrasts with the finding that E. coli does not appear to react significantly to the disruption of most single metabolic genes by regulating multiple other mRNA or protein levels. In this case, structural redundancy in the metabolic network itself provides the necessary robustness (19, 20), allowing the levels of most metabolites to remain constant, although some localized perturbations are inevitable. This strategy is likely to be less costly than the one mentioned above, because it requires no specific sensors for each mutation, while providing resistance against many of these mutations. Even if this strategy appears insufficient in the face of some mutations, E. coli may survive by accumulating additional mutations, as we observed for pfkA and rpiA disruptants. Using multiple strategies may thus enable E. coli to maintain a stable state when exposed to various types of perturbations. Additional information is available at the project Web site (21).

Supporting Online Material

Materials and Methods

SOM Text

Figs. S1 to S7

Tables S1 to S12


References and Notes

View Abstract

Stay Connected to Science

Editor's Blog

Navigate This Article