Analysis of fluorescent reporters indicates heterogeneity in glucose uptake and utilization in clonal bacterial populations
BMC Microbiology volume 13, Article number: 258 (2013)
In this study, we aimed at investigating heterogeneity in the expression of metabolic genes in clonal populations of Escherichia coli growing on glucose as the sole carbon source. Different metabolic phenotypes can arise in these clonal populations through variation in the expression of glucose transporters and metabolic enzymes. First, we focused on the glucose transporters PtsG and MglBAC to analyze the diversity of glucose uptake strategies. Second, we analyzed phenotypic variation in the expression of genes involved in gluconeogenesis and acetate scavenging (as acetate is formed and excreted during bacterial growth on glucose), which can reveal, for instance, phenotypic subpopulations that cross-feed through the exchange of acetate. In these experiments, E. coli MG1655 strains containing different transcriptional GFP reporters were grown in chemostats and reporter expression was measured with flow cytometry.
Our results suggest heterogeneous expression of metabolic genes in bacterial clonal populations grown in glucose environments. The two glucose transport systems exhibited different level of heterogeneity. The majority of the bacterial cells expressed the reporters for both glucose transporters MglBAC and PtsG and a small fraction of cells only expressed the reporter for Mgl. At a low dilution rate, signals from transcriptional reporters for acetyl-CoA synthetase Acs and phosphoenolpyruvate carboxykinase Pck indicated that almost all cells expressed the genes that are part of acetate utilization and the gluconeogenesis pathway, respectively. Possible co-existence of two phenotypic subpopulations differing in acs expression occurred at the threshold of the switch to overflow metabolism. The overflow metabolism results in the production of acetate and has been previously reported to occur at intermediate dilution rates in chemostats with high concentration of glucose in the feed.
Analysis of the heterogeneous expression of reporters for genes involved in glucose and acetate metabolism raises new question whether different metabolic phenotypes are expressed in clonal populations growing in continuous cultures fed on glucose as the initially sole carbon source.
This study focuses on the analysis of reporters for the expression of metabolic genes as a first step towards the analysis of phenotypic variation in metabolism in clonal populations of Escherichia coli. Our aim was to explore whether different systems that transport glucose exhibit different level of heterogeneity. We were also interested in whether certain conditions promote heterogeneity further downstream, in metabolic reactions. We therefore investigated if bacterial populations growing on solely glucose form phenotypically distinct subpopulations that express different parts of metabolic pathways. The rationale for these analyses was that, even under constant and homogeneous conditions, single cells can show marked differences in phenotypic traits [1, 2], including the expression of different transporters and metabolic enzymes. Such phenotypic variation can arise through a number of cellular processes; one well-studied phenomenon is ‘stochastic gene expression’ , i.e. the fact that many cellular processes are inherently variable, and that this can lead to substantial phenotypic variation that is produced independently of genetic or environmental differences [1, 4, 5].
Generally, variation in gene expression can have functional consequences and provide adaptive benefits. In situations in which the environment changes rapidly, genotypes that produce higher levels of phenotypic variation among individuals can have a higher probability to thrive [6–8]. In this study, we focus on cases in which variation in gene expression might potentially provide a different benefit. In some scenarios, it might be advantageous for cells to specialize in their metabolic function , for example due to inefficiencies or trade-offs  that arise from performing different metabolic functions within the same cell. In such cases, we might expect that individual cells within a population will either perform one function or the other, but not both. To test for instances in which we find metabolic specialization, we analyzed gene expression as a proxy for how glucose and acetate metabolism differs between single cells in clonal populations grown in glucose environments.
Previous studies have established that E. coli can employ different transport systems to take up a given carbon source from the environment. The redundancy in glucose (Glc) uptake has, in particular, been widely studied. E. coli can use five different permeases for glucose, which belong to three protein families: MglBAC is an ABC (ATP-binding cassette) transporter; GalP is a MFS (major facilitator superfamily) transporter; and PtsG/Crr, ManXYZ and NagE are parts of PTS (phosphotransferase system) [11–13]. Population-based studies have shown that the expression of a specific glucose transporter highly depends on the bacterial growth rate and the concentration of glucose in the environment [11, 12]. PtsG/Crr is the only glucose-specific PTS permease (Glc-PTS) and transcription of ptsG is induced solely by glucose . MglBAC is an uptake system that is induced by glucose and galactose, whereas GalP exhibits a wider range of specificity as it can transport different carbon sources. MglBAC and PtsG/Crr are the uptake systems that engage in most of the glucose transport in E. coli in different glucose environments [11, 12, 14–16]. The Mgl system has the leading role in glucose uptake in carbon-limited chemostat cultures. mgl expression increases rapidly at the onset to glucose depletion, particularly in glucose-limited chemostat cultures [15, 16]. MglBAC additionally allows bacteria to utilize glucose in micromolar concentrations. It is the most highly expressed transporter under glucose limitation  due to its high affinity for glucose , but PTS also transports glucose with similar micromolar affinity [12, 17, 18]. Regarding dependence of activity of glucose transporters on bacterial growth rate, at intermediate growth rates Mgl has the leading role in glucose uptake, although PtsG is active as well . Regulation of expression and activity of transporters PtsG/Crr and MglBAC is substantially different. Different groups of sigma factors, activators and repressors are responsible for regulation of their transcription, including a small RNA that additionally controls degradation of the ptsG transcript [12, 14, 19]. Furthermore, PtsG/Crr takes up and concomitantly phosphorylates glucose in an ATP-independent fashion, whereas glucose transported via ATP-dependent uptake system MglBAC is subsequently phosphorylated by a different enzyme .
Glucose is metabolized via central metabolism, which is the source of energy and biomass building blocks. First, the glycolytic enzymes break down glucose to pyruvate, which is then further metabolized to acetyl-CoA that can enter the citric acid cycle . If glucose is present in the environment as a sole carbon source, cells growing at a high rate of glucose consumption perform a fast metabolism known as overflow metabolism . The cells rapidly degrade glucose to acetyl-CoA and further to acetate, and ultimately excrete acetate . Two different pathways can catalyze the excretion of acetate: Pta-AckA (phosphate acetyltransferase - acetate kinase) during the exponential phase or PoxB (pyruvate oxidase) in the stationary phase [23, 24]. Furthermore, E. coli also has the ability to grow on acetate as a sole carbon source . Acetate can freely penetrate the cell membrane  but it also has its dedicated uptake system ActP (acetate permease) that is co-transcribed with acs encoding for acetyl-CoA synthetase . Bacteria utilize acetate by using the low affinity Pta-AckA pathway when acetate is present in high concentrations in the millimolar range. Acetyl-CoA synthetase Acs takes over acetate uptake at low concentrations of acetate in the micromolar range [21, 26]. However, the growth rate when growing solely on acetate is low: for example, the maximal growth rate on acetate is almost five times lower than on a concentration of glucose with the equivalent number of carbon atoms . In batch cultures with glucose as the sole provided carbon source, E. coli populations start to grow on the excreted acetate when glucose is depleted .
As mentioned above, acetate appears as an intermediate in reactions of glucose metabolism, and it can as well serve as a carbon source. This raises the question whether clonal bacterial populations growing in constant glucose-feed conditions form two physiologically different groups: one group that excretes acetate and has acs down-regulated and a second group that expresses acs and utilizes acetate. Such a situation would correspond to phenotypic cross-feeding. The term cross-feeding describes a metabolic interaction where the complete degradation of a substrate is partitioned between two types. One type utilizes a nutrient from the environment (e.g. glucose) and excretes the metabolized product (e.g. acetate) that is afterwards used as the primary nutrient source for the second type. Previous studies have only focused on cross-feeding between different genotypes within bacterial populations, which can spontaneously evolve in experimental microbial populations growing on glucose as the sole carbon source [28, 29]. In this study, we hypothesized that cross-feeding could also arise within an isogenic bacterial population, based on the emergence of phenotypic subpopulations with different expression of metabolic genes. Acetate cross-feeding subpopulations could potentially occur in glucose-fed clonal populations and scavenge acetate that is excreted by other cells.
Results and discussion
Different levels of phenotypic variation between different glucose transporters
Our focus was on quantifying heterogeneity in the expression of genes involved in the uptake and utilization of glucose and its metabolic intermediate acetate. We used a plasmid-based reporter system  in which fluorescence from promoter-gfp fusion constructs serves as an indirect measurement of transcription. In our recent work , we showed that signals from such plasmid-based fluorescent reporters were significantly correlated with directly measured levels of mRNA as well as with measurements of translational reporters , although the latter association was weaker. Analyses of the fluorescence of promoter-gfp reporters therefore provide partial (but not complete) information about the actual expression of a gene. We also established  that using this plasmid-based reporter system  gives comparable results of mean and variation of expression to reporter systems integrated into the chromosome.
We first investigated variation in the expression of reporters for the transporters PtsG and MglBAC, which are the most prominent glucose uptake systems in E. coli[12, 15, 16]. The aim was to test whether these glucose transporters exhibit different levels of heterogeneity in gene expression. The expression of ptsG and mglB reporters was measured in media supplemented solely with glucose (see Methods; the results are shown in Table 1, Table 2 and Additional file 1: File S1). The mean expression of PmglB-gfp was higher than PptsG-gfp in all tested glucose growth conditions (Table 1), which is consistent with previous reports that MglBAC is the most highly expressed glucose transporter at intermediate growth rates . Next, we aimed to assess whether different glucose transport strategies exhibit different levels of variation in gene expression. We used two different approaches to do so. First, we computed the coefficient of variation (CV, the ratio between the standard deviation and the mean) for each measurement of GFP fluorescence. As control, we used the reporter for rpsM, which encodes the ribosomal protein S13, previously shown to exhibit a low degree of variation in the expression between clonal cells . The ptsG reporter showed higher CVs than the mglB reporter in all glucose-feed environments (Table 2, Additional file 1: File S1), and also higher CVs than the PrpsM-gfp control (Figure 1, Table 2). However, CVs alone are not a reliable indicator for the level of heterogeneity in gene expression, since it has been previously demonstrated that CVs are dependent on the mean expression level . This relationship also manifests in our dataset in all tested growth conditions (presented in the next section of Results and Discussion).
We thus used a second measure for variation in gene expression: the fraction of cells in a clonal population that expressed the transcriptional reporter above background levels. We subtracted the background fluorescence (log10 value of 1.3; see Methods) from the measurements of expression of PptsG-gfp and PmglB-gfp, for all growth conditions that we tested. Expression of PmglB-gfp was above background in 90.1-99.8% of the cells within a population (one measurement for each environmental condition presented in Table 3; Additional file 1: File S1), depending on the growth conditions. This implies that the vast majority of cells transcribe mglBAC regardless of the carbon sources present in the media or the growth rate. Considering only cultures grown on glucose, 96.8-99.8% of the population expressed the mglB reporter above background. In the same conditions, the fraction of cells that did not express PptsG-gfp was in two cases above 5%. For instance, 8.6% of the cells in the population that was grown in the chemostats cultures  at D = 0.15 h-1 with 0.56 mM Glc did not express PptsG-gfp. It is conceivable that a subfraction of the cells that do not express PptsG-gfp is metabolically inactive. To test this, we compared the fraction of cells that does not express PptsG-gfp with the fraction of cells that does not express the ribosomal reporter PrpsM-gfp, measured under the same conditions. The ribosomal reporter indicated that only around 0.5% of the population did not transcribe the ribosomal protein (Table 3), i.e. those were probably dead or not actively dividing cells. This indirectly implies that most of the cells that did not express PptsG-gfp may be metabolically active and should thus engage in another glucose uptake strategy.
Overall, these results suggest that the promoter for mglBAC is expressed above background in a higher fraction of the population than the promoter for ptsG, and differences in ptsG expression between genetically identical cells could be an indication of glucose uptake heterogeneity within clonal populations.
Next, we used direct measurements of uptake to analyze the activity of the glucose-PTS transporter and to compare the transporter activity with the expression of PptsG-gfp. 2-NBDG, 2-[N-(7-nitrobenz-2-oxa-1,3-diazol-4-yl)amino]-2-deoxy-D-glucose, is a fluorescent D-glucose analog, and has been used to study the dynamics of glucose uptake via the phosphotransferase system (PTS) in single cells of E. coli[18, 34]. Since 2-NBDG is exclusively taken up via Glc-PTS, cells will fluoresce only if their PTS system is active and the glucose analog is transported inside the cell. As this assay uses a glucose analog that cannot be metabolized, the results can be interpreted only in the context of the activity of the transport system and not as a general measure of metabolic activity of a cell. Our data indicate that not all cells use the PTS system to take up glucose from the media (Figure 2, medium supplemented with 0.56 mM Glc). How do the rest of the cells take up glucose – do they maybe employ alternative carbon sources? There are two possibilities. First, cells might use Mgl or another glucose transporters. Second, it is possible that the cells use excreted acetate as (an additional) carbon source. We also found that even if the PptsG-gfp reporter strain fluoresces, it does not necessarily mean that PTS is actively transporting glucose (Figure 2). This became evident in control experiments where we grew cells in medium containing acetate or arabinose as the sole carbon source. Around 80% of the gated population growing in acetate (around 60% growing in arabinose) expressed the ptsG reporter above the background level, without any glucose present to induce the expression or to be transported (Additional file 1: File S1). Furthermore, in these conditions the PptsG-gfp reporter showed a high degree of variation in expression (Figure 2).
Transcriptional reporters for glucose transporters can only provide limited insights into the actual metabolic state of cells. Several recent papers have discussed discrepancies between transcriptional reporters and metabolic fluxes in specific parts of metabolic pathways [35, 36]. As a consequence, we need to be cautious when using data from transcriptional reporters to make inferences about the actual physiology of cells. Additional experiments could provide complementary insights, for instance the analysis of sugar transporter synthesis or activity, together with analysis of sugar assimilation at the single-cell level .
Variation in the expression of glucose transporters across environments
We next investigated how the variation in expression of reporters for different glucose transporters changes across different environments. We first compared the results of this study with the results from a genome-wide study of promoter-mediated phenotypic variation . Mean and variation of the expression of ptsG, mglB and rpsM reporters are shown in Figure 3 (plotted are mean values of replicates in different conditions). When power regression lines were fitted across different expression data from the same environment, all lines showed the same trend, namely that the CV of log fluorescence values decreased with mean log GFP expression (Figure 3). Our analysis suggests some general rules: variation in the expression from these three promoters was lowest in batch cultures supplemented with glucose, or glucose plus acetate, and highest in batch or chemostats cultures with acetate as a sole carbon source.
Our results also show that, while mean and variation of reporter expression in a given environment follow a clear and simple pattern (Figure 3), measurements of mean and variation in expression of a particular reporter in different environments deviate substantially from this pattern (Additional file 2: Figure S1). The specifics of these deviations were dependent on the reporter we analyzed: ptsG showed a negative association between mean expression and the variation of expression across environments, while mglB showed a positive association. We speculate that these differences between ptsG and mglB could be a consequence of distinctive regulatory features of the glucose transporters [12–15, 17, 19], different affinity towards transported sugar [12, 17], and possibly different growth rate dependencies .
Variation in the expression of genes involved in glucose and acetate utilization
Besides exhibiting heterogeneity in uptake of glucose, cells could show phenotypic variation in the expression of metabolic genes involved in utilization of glucose and acetate. In particular, we were interested in gene expression patterns that could indicate variation between cells in the consumption of acetate; in our system, acetate can come from two different sources – from the same cell or taken up from the environment where it is excreted by other cells. As discussed in the Background, the presence of cells that take up acetate produced by other cells would be indicative of phenotypic cross-feeding in clonal populations. To investigate this, we constructed a Pacs-gfp reporter to measure the expression of the gene encoding for acetyl-CoA synthetase Acs. Generally, rapid increase in acs transcription occurs when bacterial cultures are inoculated into medium containing solely acetate as a carbon source . The promoter Pacs controls the acs-yjcH-actP operon, and hence also controls transcription of the acetate permease ActP . Therefore, differential regulation of acs can also indicate altered expression of the acetate transporter and regulation of the uptake of external acetate. However, uptake via ActP is not the only acetate uptake strategy, since acetate can freely diffuse into cells . The expression of acs is down-regulated when bacteria excrete acetate  and up-regulated when bacteria utilize acetate . Accordingly, we detected increased expression of the acs reporter when bacteria were grown only on acetate in comparison to growth on glucose (Figure 4, Additional file 1: File S1). Moreover, the expression of the acs reporter was reduced when the concentration of glucose in the chemostat feed was increased (Figure 4). This is consistent with previous reports that have shown that high concentrations of glucose lead to an increase in the intracellular concentration of acetate , resulting in down-regulation of the acs operon. Furthermore, we observed that variation in Pacs-gfp expression was actually lower in glucose-limited chemostats than in acetate-limited continuous cultures (Table 2) or in glucose-acetate continuous cultures (Additional file 3: Text S1, Additional file 4: Figure S2).
Results from previous studies suggest that under the conditions used here – glucose as the only carbon source, and low dilution rates – the reactions of glyoxylate shunt and gluconeogenesis should be active, which would allow utilization of simple carbon sources such as acetate when glucose is not available . According to population-based studies on bacteria grown on glucose, the shunt operates at the dilution rates from 0.05–0.2 h-1, allowing metabolism of acetyl-CoA to succinate. The reactions of the citric acid cycle are not engaged, and this prevents carbon loss in the form of CO2[33, 41]. When acetate is used as a sole carbon source, the expression of the phosphoenolpyruvate (PEP) carboxykinase gene pck (a gluconeogenesis enzyme) is up-regulated [40, 42], indicating synthesis of glucose from non-carbohydrate precursors such as acetate . pck is also up-regulated in chemostats containing glucose as a carbon source that are run at low dilution rates .
Our experiments at the single-cell level largely support these previous population-based studies. In the following paragraph, we will discuss in more details the gene expression phenotypes that we observed in clonal populations grown in mini-chemostats at low dilution rate of D = 0.15 h-1, and with glucose as the sole carbon source at a feed concentration of 0.56 mM Glc. These are the conditions in which the majority of the cells expressed both glucose transporters mglB and ptsG, whereas some cells only expressed mglB (Figure 1, Table 3). The fraction of cells that did not express the ribosomal reporter was below 1% (Table 3), and these were the cells that presumably did not grow and divide. The residual concentration of glucose in the mini-chemostats after five volume changes (theoretical steady-state concentration ) was 1.95 ± 0.13 μM, measured by ion chromatography (our experimental setup did not allow us to accurately measure concentration of acetate).
We found that, under these conditions, almost all cells expressed the acs reporter above background level (Figure 4). This may indicate that they either recover cytoplasmic acetate or take up acetate excreted by others. Recovering cytoplasmic acetate can be consequence of homeostasis of acetyl-AMP and acetylphosphate . In chemostats run under such conditions, acetate is usually not detected [43–45], however it might be possible that scarce amounts of acetate are excreted and immediately taken up by an acetate cross-feeding subpopulation. It has been argued that the production of acetate is independent of the growth rate and that the growing bacteria can simultaneously produce and utilize acetate [45, 46]. The expression of the pck reporter also indicates that most of the cells possibly engaged in the reactions of gluconeogenesis (Additional file 5: Figure S3).
Previous studies provided evidence that transcriptional regulation does indeed have a significant impact on the direction of the metabolic flux through the pyruvate/acetyl-CoA node . Transcriptional control at this branching point allows flux to proceed via overflow metabolism, citric acid cycle and/or PEP-glyoxylate cycle . Results presented in another paper indicate that alterations of fluxes through the glyoxylate shunt and the citric acid cycle were associated with changes in the expression of these genes . Therefore, transcriptional reporters for acetate metabolism (the acs reporter) and PEP-glyoxylate pathway (the pck reporter) may indeed be indicative of the fluxes through those pathways.
Switching to overflow metabolism and bimodal expression of the acsreporter
It has been shown that the excretion of acetate (overflow metabolism) occurs in chemostat populations at a dilution rate of about 0.3 h-1[22, 44]. Increasing the concentration of glucose in the chemostat feed results in intensified production of acetate . Our results support the existence of overflow metabolism at D = 0.3 h-1 in chemostats with high concentrations (5.6 mM) of glucose in the feed. Under these conditions, decreased expression of acs and pck reporters indicated that assimilation of acetate was reduced and gluconeogenesis was shut down (Figure 5). However, not all replicate cultures showed consistent patterns in the expression of transcriptional reporters. The expression of the reporters for mglB and acs was not consistent between different experiments, in contrast to the measurements for rpsM, ptsG and pck (Figure 5). This suggests that not all replicate cultures switched to the overflow metabolism, possibly due to the fact that the mini-chemostats were operated at the threshold of the expected switch to overflow metabolism.
Several replicates showed bimodal patterns of the expression of Pacs-gfp in well-mixed chemostat cultures. This is consistent with the idea that within clonal populations two phenotypically different subpopulations existed (Figure 5) – a first group of cells that presumably scavenged acetate and expressed the acs reporter and a second group that excreted acetate and thus down-regulated expression of acs. According to this scenario, the first subpopulation performed metabolic reactions indicative of carbon source limitation whereas the expression profiles of metabolic genes in the second subpopulation did not reflect glucose-limited conditions. These results potentially support the existence of phenotypic subpopulations that engage in acetate cross-feeding, as hypothesized above. However, it is also possible that both phenotypic subpopulations utilize glucose as the primary carbon source (since the expression of the pck reporter was only slightly above background, Figure 5) and the first subpopulation additionally recovers cytoplasmic acetate to increase intracellular levels of acetyl-AMP and acetylphosphate . Future experiments with advanced continuous cultivation methods (e.g. accelerostat cultivation as described in ) would be valuable for further refining the environmental conditions where these two metabolic strategies co-exist.
Many studies refer to glucose-limited chemostats as “simple conditions”, e.g. [28, 29, 48]. Even though glucose serves as a sole carbon source in these experiments, the metabolic regimes of the populations of E. coli are far away from “simple” . Each cell within the bacterial population can take up glucose via five different transporters and metabolize it according to its needs for biomass building blocks and energy. Glucose is broken down to metabolic intermediates including acetate; acetate can be recovered in the central metabolic pathway, or it can be excreted and potentially then scavenged by other cells. Our results show that single cells within clonal population differ in their gene expression patterns and thus potentially in metabolic phenotypes when only glucose is supplied in the feed. This variation can arise through 1) different expression of glucose transporters (PtsG/Crr, MglBAC, etc.) between individual cells, 2) differences in utilization of acetate recovered within the cells and potentially, 3) uptake of excreted acetate. While the fluorescent reporter systems that we used in this study can provide insights into cell-to-cell variation in gene expression, and thus be used for the generation of hypotheses about metabolic variation, addition experiments are needed to test these hypotheses and to get direct insights into metabolic activities of cells growing on glucose.
E. coli is commonly used for the production of recombinant proteins and other valuable products, and the corresponding cultures are usually grown at high growth rates. High consumption of glucose is often associated with the excretion of acetate that inhibits recombinant protein production [44, 45]. The findings presented here can provide a better understanding of the strategies involved in metabolizing glucose (as the only carbon-source component of the medium) and acetate that is subsequently produced during glucose utilization, and thus contribute to the development of new strategies for improving growth of industrial strains.
All E.coli K-12 MG1655  strains with reporter plasmids used in this study are listed in Table 4. The strain containing the plasmid with the reporter Pacs-gfp was constructed as follows. A 858 bp-long intergenic region (comprising the region between acs and nrfA and the parts of the open reading frames) was amplified from the MG1655 chromosome using the primers Fwd_Pacs_XhoI 5’-CCGCTCGAGTAAGCTGAAGATACGGCGTGC-3’ and Rev_Pacs_BamHI 5’-CGGGATCCCCATCGGCATATAAATCGCCACC-3’ (italic parts of sequences are the restriction sites). The construct was cloned via XhoI/BamHI restriction into the plasmid containing the PptsG-gfp reporter  (thus swapping the existing ptsG promoter) and transformed into MG1655.
The cultures were grown in M9 minimal medium (Sigma-Aldrich) containing 47.76 mM Na2HPO4, 23.6 mM KH2PO4, 8.56 mM NaCl and 20.2 mM NH4Cl. 1 mL of 1 M MgSO4 (Fluka) and 100 μL of 1 M CaCl2 (Sigma-Aldrich) were added to 1 L of minimal medium. D(+)-glucose (Sigma) and/or sodium acetate (Fluka) were used as carbon source(s) and added to the desired concentration. The concentration of kanamycin sulfate (Sigma) was 50 μg/mL.
Cultivation in the chemostats
Frozen clones were first streaked on LB agar (Sigma-Aldrich) plates to obtain single colonies. The agar plates contained 50 μg/mL of kanamycin for reporter strains . A single colony was inoculated overnight in defined minimal medium (total 4 mL). 1 mL of these precultures was used to inoculate each mini-chemostat (total 5.5 mL) . The minimal speed of the inflow pump corresponding to a dilution rate of D = 0.14 h-1 was increased in 2 or 3 steps until a dilution rate of D = 0.15 h-1 was reached after 24 h (using the peristaltic pump IPC-N from Ismatec, IDEX Health & Science, Germany). The airflow was maintained with the outflow pump (model IP from Ismatec, IDEX Health & Science, Germany) at 20 mL per minute with filter-sterilized water-saturated air . Continuous formation of air-bubbles as well as small magnetic stirrer bars within the mini-chemostats ensured sufficient mixing of the bacterial cultures. The chemostats were harvested after 5 volume changes (one volume change every 6.67 hours) at the final dilution rate, i.e. after reaching the steady state  (Additional file 6: Figure S4).
For the experiments performed at D = 0.3 h-1 the total run-time was adjusted to the same number of volume changes as obtained with the experiments performed at D = 0.15 h-1.
Frozen clones were first streaked on LB agar plates (containing kanamycin when needed). A single colony was inoculated overnight in defined minimal medium (total 4 mL). The overnight cultures were diluted 200-fold into 4 mL of minimal medium and grown for 2 hours before measured in the flow cytometer.
We analyzed GFP fluorescence as a proxy for gene expression. For the strains grown in mini-chemostats, the GFP fluorescence was measured after 5 volume changes, which are required to reach steady state  (Additional file 6: Figure S4) but short enough to minimize the probability of mutations in the promoter region. GFP fluorescence was measured in the early exponential phase for the samples grown in the batch cultures.
All measurements were performed 2–5 times, as independent replicates coming from different overnight cultures. (For analysis of overflow metabolism we measured up to 20 replicates.)
We used the PAS-III flow cytometer (Partec, Muenster, Germany) equipped with 488 nm excitation laser. The green fluorescence was measured at 520 nm (FL1 channel), with the following settings – FSC (forward scatter): 270 V, SSC (side scatter): 210 V, FL1: 600 V, speed = 3, trigger on SSC.
For each strain and all growth conditions, raw data were processed using FlowJo software version 8.8.7 (Tree Star, Inc.), and gated on 10,000-12,000 cells by using the autogating tool in the densest area of the pseudo-color plots of SSC vs. FSC. These gated cells were then used for the subsequent analysis. For analysis of the negative controls (strains with the promoterless plasmid pUA66 or wild-type MG1655) no gating was applied.
The cells were considered not to express a reporter when their fluorescence values were below the background fluorescence. The background fluorescence was defined as the mean value of the 99th percentile of fluorescence intensities (Additional file 1: File S1) of the strain with the promoterless plasmid pUA66 (no gating applied) measured in various environments.
The fluorescence values for the cells within the gated populations were log10 transformed for the analysis, and thus we computed mean log expression and CV (coefficient of variation, the ratio between standard deviation and mean) of log expression.
Influence of data filtering on the results
We restricted our analysis to the fraction of cells that were in similar physiological activity and size [31, 51, 52]. The cells were gated within a narrow range of defined flow cytometry parameters. We analyzed how the number of cells in the gated fraction influences the computation of mean and CV. One sample (the measurement of the strain harboring PmglB-gfp in the chemostats cultures at D = 0.15 h-1, with 5.6 mM Glc feed) was, therefore, gated 24 times (Additional file 7: Figure S5) while varying cell number in the range 5,000-20,000 cells.
E.coli K-12 MG1655  and the PptsG-gfp strain from the plasmid library  were used for these experiments. The strains were grown in the mini-chemostats  with minimal media supplemented with a sole carbon source (0.56 mM sodium acetate, 0.56 mM L-arabinose (Sigma-Aldrich), 0.56 mM D-glucose or 5.6 mM D-glucose). After 5 volume changes at D = 0.15 h-1, cells were harvested. Fluorescence was measured with the flow cytometer, as described above. PptsG-gfp fluorescence was measured immediately upon harvesting. MG1655 samples were incubated with 10 μM 2-NBDG (Molecular Probes, Life Technologies) for 5 minutes according to , and their fluorescence was measured directly afterwards.
We analyzed glucose concentration by ion chromatography using Dionex DX-500 system with CarboPack PA10 carbohydrate column. The eluent was 200 mM NaOH, and the calibration curves were obtained by measuring glucose solutions of known concentration.
The data were analyzed in SPSS statistical software version 19 and Microsoft Excel version 14.3.
Green fluorescent protein
Coefficient of variation
Major facilitator superfamily
Davidson CJ, Surette MG: Individuality in Bacteria. Annu Rev Genet. 2008, 42: 253-268. 10.1146/annurev.genet.42.110807.091601.
Veening JW, Smits WK, Kuipers OP: Bistability, epigenetics, and bet-hedging in bacteria. Annu Rev Microbiol. 2008, 62: 193-210. 10.1146/annurev.micro.62.081307.163002.
Elowitz MB, Levine AJ, Siggia ED, Swain PS: Stochastic gene expression in a single cell. Science. 2002, 297: 1183-1186. 10.1126/science.1070919.
Raser JM, O'Shea EK: Noise in gene expression: Origins, consequences, and control. Science. 2005, 309: 2010-2013. 10.1126/science.1105891.
Raj A, van Oudenaarden A: Nature, nurture, or chance: stochastic gene expression and its consequences. Cell. 2008, 135: 216-226. 10.1016/j.cell.2008.09.050.
Kussell E, Leibler S: Phenotypic diversity, population growth, and information in fluctuating environments. Science. 2005, 309: 2075-2078. 10.1126/science.1114383.
Acar M, Mettetal JT, van Oudenaarden A: Stochastic switching as a survival strategy in fluctuating environments. Nat Genet. 2008, 40: 471-475. 10.1038/ng.110.
Arnoldini M, Mostowy R, Bonhoeffer S, Ackermann M: Evolution of stress response in the face of unreliable environmental signals. PLOS Comput Biol. 2012, 8 (8): e1002627-10.1371/journal.pcbi.1002627.
Johnson DR, Goldschmidt F, Lilja EE, Ackermann M: Metabolic specialization and the assembly of microbial communities. ISME J. 2012, 6: 1985-1991. 10.1038/ismej.2012.46.
Molenaar D, van Berlo R, de Ridder D, Teusink B: Shifts in growth strategies reflect tradeoffs in cellular economics. Mol Syst Biol. 2009, 5: 323-
Ferenci T: Adaptation to life at micromolar nutrient levels. FEMS Microbiol Rev. 1996, 18: 301-317. 10.1111/j.1574-6976.1996.tb00246.x.
Jahreis K, Pimentel-Schmitt EF, Bruckner R, Titgemeyer F: Ins and outs of glucose transport systems in eubacteria. FEMS Microbiol Rev. 2008, 32: 891-907. 10.1111/j.1574-6976.2008.00125.x.
Keseler IM, Collado-Vides J, Santos-Zavaleta A, Peralta-Gil M, Gama-Castro S, et al: EcoCyc: a comprehensive database of Escherichia coli biology. Nucleic Acids Res. 2011, 39: D583-D590. 10.1093/nar/gkq1143.
Plumbridge J: Expression of ptsG, the gene for the major glucose PTS transporter in Escherichia coli, is repressed by Mlc and induced by growth on glucose. Mol Microbiol. 1998, 29: 1053-1063. 10.1046/j.1365-2958.1998.00991.x.
Ferenci T: Hungry bacteria—definition and properties of a nutritional state. Environ Microbiol. 2001, 3: 605-611. 10.1046/j.1462-2920.2001.00238.x.
Steinsiek S, Bettenbrock K: Glucose transport in Escherichia coli mutant strains with defects in sugar transport systems. J Bacteriol. 2012, 194: 5897-5908. 10.1128/JB.01502-12.
Seeto S, Notley-McRobb L, Ferenci T: The multifactorial influences of RpoS, Mlc and cAMP on ptsG expression under glucose-limited and anaerobic conditions. Res Microbiol. 2004, 155: 211-215. 10.1016/j.resmic.2003.11.011.
Natarajan A, Srienc F: Dynamics of glucose uptake by single Escherichia coli cells. Metab Eng. 1999, 1: 320-333. 10.1006/mben.1999.0125.
Gama-Castro S, Salgado H, Peralta-Gil M, Santos-Zavaleta A, Muniz-Rascado L, et al: RegulonDB version 7.0: transcriptional regulation of Escherichia coli K-12 integrated within genetic sensory response units (Gensor Units). Nucleic Acids Res. 2011, 39: D98-D105. 10.1093/nar/gkq1110.
Madigan MT, Martinko JM, Stahl DA, Clark DP: Brock biology of microorganisms. 2012, San Francisco: Pearson Education Inc., 13
Wolfe AJ: The acetate switch. Microbiol Mol Biol Rev. 2005, 69: 12-50. 10.1128/MMBR.69.1.12-50.2005.
El-Mansi EM, Holms WH: Control of carbon flux to acetate excretion during growth of Escherichia coli in batch and continuous cultures. J Gen Microbiol. 1989, 135: 2875-2883.
Dittrich CR, Bennett GN, San KY: Characterization of the acetate-producing pathways in Escherichia coli. Biotechnol Prog. 2005, 21: 1062-1067.
Castano-Cerezo S, Pastor JM, Renilla S, Bernal V, Iborra JL, et al: An insight into the role of phosphotransacetylase (pta) and the acetate/acetyl-CoA node in Escherichia coli. Microb Cell Fact. 2009, 8: 54-10.1186/1475-2859-8-54.
Gimenez R, Nunez MF, Badia J, Aguilar J, Baldoma L: The gene yjcG, cotranscribed with the gene acs, encodes an acetate permease in Escherichia coli. J Bacteriol. 2003, 185: 6448-6455. 10.1128/JB.185.21.6448-6455.2003.
Kumari S, Beatty CM, Browning DF, Busby SJ, Simel EJ, et al: Regulation of acetyl coenzyme A synthetase in Escherichia coli. J Bacteriol. 2000, 182: 4173-4179. 10.1128/JB.182.15.4173-4179.2000.
Liu M, Durfee T, Cabrera JE, Zhao K, Jin DJ, et al: Global transcriptional programs reveal a carbon source foraging strategy by Escherichia coli. J Biol Chem. 2005, 280: 15921-15927. 10.1074/jbc.M414050200.
Rosenzweig RF, Sharp RR, Treves DS, Adams J: Microbial evolution in a simple unstructured environment: genetic differentiation in Escherichia coli. Genetics. 1994, 137: 903-917.
Treves DS, Manning S, Adams J: Repeated evolution of an acetate-crossfeeding polymorphism in long-term populations of Escherichia coli. Mol Biol Evol. 1998, 15: 789-797. 10.1093/oxfordjournals.molbev.a025984.
Zaslaver A, Bren A, Ronen M, Itzkovitz S, Kikoin I, et al: A comprehensive library of fluorescent transcriptional reporters for Escherichia coli. Nat Methods. 2006, 3: 623-628. 10.1038/nmeth895.
Silander OK, Nikolic N, Zaslaver A, Bren A, Kikoin I, et al: A genome-wide analysis of promoter-mediated phenotypic noise in Escherichia coli. PLOS Genet. 2012, 8: e1002443-10.1371/journal.pgen.1002443.
Taniguchi Y, Choi PJ, Li G-W, Chen H, Babu M, et al: Quantifying E. coli proteome and transcriptome with single-molecule sensitivity in single cells. Science. 2010, 329: 533-538. 10.1126/science.1188308.
Nanchen A, Schicker A, Sauer U: Nonlinear dependency of intracellular fluxes on growth rate in miniaturized continuous cultures of Escherichia coli. Appl Environ Microbiol. 2006, 72: 1164-1172. 10.1128/AEM.72.2.1164-1172.2006.
Natarajan A, Srienc F: Glucose uptake rates of single E. coli cells grown in glucose-limited chemostat cultures. J Microbiol Meth. 2000, 42: 87-96. 10.1016/S0167-7012(00)00180-9.
van Rijsewijk BRB H, Nanchen A, Nallet S, Kleijn RJ, Sauer U: Large-scale 13C-flux analysis reveals distinct transcriptional control of respiratory and fermentative metabolism in Escherichia coli. Mol Syst Biol. 2011, 7: 477-
Kochanowski K, Sauer U, Chubukov V: Somewhat in control—the role of transcription in regulating microbial metabolic fluxes. Curr Opin Biotech. 2013, in press
Musat N, Foster R, Vagner T, Adam B, Kuypers MMM: Detecting metabolic activities in single cells, with emphasis on nanoSIMS. FEMS Microbiol Rev. 2012, 36: 486-511. 10.1111/j.1574-6976.2011.00303.x.
Cases I, de Lorenzo V: Expression systems and physiological control of promoter activity in bacteria. Curr Opin Microbiol. 1998, 1 (3): 303-310. 10.1016/S1369-5274(98)80034-9.
Veit A, Polen T, Wendisch V: Global gene expression analysis of glucose overflow metabolism in Escherichia coli and reduction of aerobic acetate formation. Appl Microbiol Biotechnol. 2007, 74: 406-421. 10.1007/s00253-006-0680-3.
Oh M-K, Rohlin L, Kao KC, Liao JC: Global expression profiling of acetate-grown Escherichia coli. J Biol Chem. 2002, 277 (15): 13175-13183. 10.1074/jbc.M110809200.
Sauer U, Eikmanns BJ: The PEP-pyruvate-oxaloacetate node as the switch point for carbon flux distribution in bacteria. FEMS Microbiol Rev. 2005, 29: 765-794. 10.1016/j.femsre.2004.11.002.
Peng L, Shimizu K: Global metabolic regulation analysis for Escherichia coli K12 based on protein expression by 2-dimensional electrophoresis and enzyme activity measurement. Appl Microbiol Biotechnol. 2003, 61: 163-178.
Fischer E, Sauer U: A novel metabolic cycle catalyzes glucose oxidation and anaplerosis in hungry Escherichia coli. J Biol Chem. 2003, 278: 46446-46451. 10.1074/jbc.M307968200.
Valgepea K, Adamberg K, Nahku R, Lahtvee PJ, Arike L, et al: Systems biology approach reveals that overflow metabolism of acetate in Escherichia coli is triggered by carbon catabolite repression of acetyl-CoA synthetase. BMC Syst Biol. 2010, 4: 166-10.1186/1752-0509-4-166.
Renilla S, Bernal V, Fuhrer T, Castano-Cerezo S, Pastor JM, et al: Acetate scavenging activity in Escherichia coli: interplay of acetyl-CoA synthetase and the PEP-glyoxylate cycle in chemostat cultures. Appl Microbiol Biotechnol. 2012, 93: 2109-2124. 10.1007/s00253-011-3536-4.
Shin S, Chang D, Pan JG: Acetate consumption activity directly determines the level of acetate accumulation during Escherichia coli W3110 growth. J Microbiol Biotechnol. 2009, 19: 1127-1134.
Fong SS, Nanchen A, Palsson BO, Sauer U: Latent pathway activation and increased pathway capacity enable Escherichia coli adaptation to loss of key metabolic enzymes. J Biol Chem. 2006, 281: 8024-8033. 10.1074/jbc.M510016200.
Kinnersley MA, Holben WE, Rosenzweig F: E unibus plurum: Genomic analysis of an experimentally evolved polymorphism in Escherichia coli. PLOS Genet. 2009, 5: e1000713-10.1371/journal.pgen.1000713.
Notley-McRobb L, Ferenci T: The generation of multiple co-existing mal-regulatory mutations through polygenic evolution in glucose-limited populations of Escherichia coli. Environ Microbiol. 1999, 1: 45-52. 10.1046/j.1462-2920.1999.00003.x.
Blattner FR, Plunkett G, Bloch CA, Perna NT, Burland V, et al: The complete genome sequence of Escherichia coli K-12. Science. 1997, 277: 1453-1462. 10.1126/science.277.5331.1453.
Tsuru S, Ichinose J, Kashiwagi A, Ying BW, Kaneko K, et al: Noisy cell growth rate leads to fluctuating protein concentration in bacteria. Phys Biol. 2009, 6: 036015-10.1088/1478-3975/6/3/036015.
Freed NE, Silander OK, Stecher B, Böhm A, Hardt WD, et al: A simple screen to identify promoters conferring high levels of phenotypic noise. PLOS Genet. 2008, 4: e1000307-10.1371/journal.pgen.1000307.
We thank Tobias Bergmiller for great discussions, Karl Kochanowski for useful comments, Olin K Silander for critical feedback on the manuscript, Luise Wolf for technical help with data analysis and Thomas Fleischmann for technical assistance with ion chromatography.
NN and MA were supported by the Swiss National Science Foundation grant 31003A_130735.
The authors declare that they have no competing interests.
Conceived and designed the experiments: NN MA. Performed the experiments: NN TB. Analyzed the data: NN TB MA. Wrote the manuscript: NN MA. All authors read and approved the final manuscript.
Electronic supplementary material
Additional file 2: Figure S1: Variation in the expression of ptsG, mglB and rpsM reporters across different environments. The CV of log expression of PptsG-gfp (green), PmglB-gfp (blue) and PrpsM-gfp (red) was plotted against the mean log expression. Power regression was fitted to each dataset corresponding to the expression of the same reporter across different environments. The individual curves of variation in the expression of ptsG and rpsM reporters showed negative associations between the mean expression and the variation of expression across environments, whereas the mglB reporter showed a positive association. (TIFF 145 KB)
Additional file 4: Figure S2: Reporter expression in mixed-substrate environments. Expression of ptsG, mglB and acs reporters was measured in chemostats (D = 0.15 h-1) in mixed-substrate environments supplemented with 0.28 mM Glc and 0.28 mM Ac (green), or 2.8 mM Glc and 2.8 mM Ac (blue). The distributions were plotted together with the measurements of the reporter expression in the environments with only glucose in the feed (0.56 mM Glc – orange, and 5.6 mM Glc – red). The fluorescence of the promoterless strain is presented in black. (TIFF 544 KB)
Additional file 5: Figure S3: Expression of the pck reporter in different chemostat and batch conditions. Ppck-gfp fluorescence (indication of flux to gluconeogenesis) was measured in bacterial populations grown in chemostats (D = 0.15 h-1) and batch environments supplied with minimal media supplemented with only D-glucose, only sodium acetate or D-glucose plus sodium acetate. Again, background fluorescence is the fluorescence of the promoterless strain, depicted in black. The expression of the pck reporter was decreased in the exponential phase in glucose batch cultures in comparison to carbon-limited chemostats. (TIFF 757 KB)
Additional file 6: Figure S4: Changes in gfp expression prior of reaching theoretical steady-state. Pacs-gfp fluorescence was measured for five independent replicates growing on different concentration of glucose in the feed. At time point of 0 hours, chemostat experiments were started at a minimal dilution rate of D = 0.14 h-1. After 24 hours, dilution rates were increased to D = 0.15 h-1. The fluorescence plots show gfp distribution in bacterial populations without gating, together with fluorescence of the promoterless strain depicted in black. All independent replicates showed reproducible measurements of GFP fluorescence after 3.6 volume turnovers at D = 0.15 h-1. (TIFF 2 MB)
Additional file 7: Figure S5: Influence of size of the gate on the mean and CV. The strain carrying PmglB-gfp was grown in chemostats (at D = 0.15 h-1, with 5.6 mM Glc) and analyzed with flow cytometry. A) For subsequent analysis, the cells were gated using the autogating tool (FlowJo, Tree Star, Inc.) in the densest area of the pseudo-color plots of SSC vs. FSC. B) The gating was performed 24 times to capture between 5,000-20,000 cells, and the resulting distributions of GFP fluorescence were plotted. This yielded mean log expression of 2.69 ± 0.005 (mean ± standard deviation) and CV was 0.13 ± 0.0014. This suggests that the results for mean expression and CV deviated less than 1% when gate size was varying 4-fold. Our gate size varied maximally 1.2-fold when analyzing 10,000-12,000 cells, therefore the slight differences in the gate size should minimally influence the computation of mean and CV. (TIFF 681 KB)
About this article
Cite this article
Nikolic, N., Barner, T. & Ackermann, M. Analysis of fluorescent reporters indicates heterogeneity in glucose uptake and utilization in clonal bacterial populations. BMC Microbiol 13, 258 (2013). https://doi.org/10.1186/1471-2180-13-258