Bioreactor mixing efficiency modulates the activity of a prpoS::GFP reporter gene in E. coli
Microbial Cell Factoriesvolume 8, Article number: 15 (2009)
Extensive studies have shown that up-scaling of bioprocesses has a significant impact on the physiology of the microorganisms. Among the factors associated with the fluid dynamics of the bioreactor, concentration gradients induced by loss of the global mixing efficiency associated with the increasing scale is the main phenomena leading to strong physiological modifications at the level of the microbial population. These changes are not fully understood since they involve complex physiological mechanisms. In this work, we intend to investigate, at the single cell level, the expression of the rpoS gene associated with the stress response of E. coli. The cultures of the reporter strain have been performed in a small scale reactor as well as in a series of scaled-down bioreactors able to induce extracellular perturbations with increasing level of magnitude.
The rpoS level has been monitored by the aim of a transcriptional reporter gene based on the synthesis of the green fluorescent protein (GFP). It has been observed that the level of GFP increases during the transition from batch to fed-batch phase. After this initial increase, the GFP content of the cell drops, primarily due to the dilution by cell division. However, a significant drop of the GFP content has been observed if using a partitioned bioreactor, for which the mixing conditions are very bad, leading to the exposure of the cells to cyclic and stochastic extracellular fluctuations. If considering the flow cytometric profile of the cell to cell GFP content, this drop has to be attributed to the appearance of segregation at the level of the GFP content among the microbial population.
The generation of extracellular perturbations (in the present case, at the level of the sugar concentration and the dissolved oxygen level) has led to a drop at the level of the rpoS expression level. This drop has to be attributed to a segregation phenomenon in microbial population, with a major sub-population exhibiting a low expression level and a minor sub-population keeping its initial elevated expression level. The intensity of the segregation, as well as its time of appearance during the culture can be related to the bioreactor mixing efficiency.
Bioreactors are designed in order to promote an optimal environment for growth and/or metabolites synthesis. However, during scale-up, the hydrodynamic efficiency is strongly altered, potentially leading to cell exposure to heterogeneities . Microbial cells are highly sensitive to these heterogeneities and are able to respond by several physiological mechanisms involving a complex set of cellular networks (metabolome, proteome, genome), each exhibiting their own spatial and temporal organization . At this time, several cellular mechanisms associated with the bioreactor heterogeneities remain to be elucidated. These mechanisms can be observed at the lab scale by using the so-called scale-down strategy . The scale-down reactors are able to reproduce at small scale the extracellular fluctuations experienced by the cells in large-scale bioreactors. Extracellular fluctuations perceived by the cells depend on the interrelation between two fluid dynamics related processes, i.e. the global mixing efficiency of the reactor governing the spatial intensity of the concentration gradient field (for an aerobic process, mainly the glucose and the dissolved oxygen level), and the circulation paths followed by cells. To this end, a precise definition of the extracellular environment perceived by the microbial cells requires the superimposition of the two hydrodynamic mechanisms [4, 5]. At this level, difficulty comes from the fact that the circulation process is partly governed by random phenomena, giving rise to a stochastic process at the level of the heterogeneities perceived by the cells [5–7]. Several scale-down apparatus have been described in order to reproduce these heterogeneities , i.e. controlled reactors , two interconnected lab-scale stirred bioreactors [9–11] and lab-scale stirred reactor connected with a tubular section [12, 13], all these systems being based on the principle that the microbial cells must be submitted to oscillating or stochastic (partially random) fluctuations to reproduce the flow conditions encountered at large-scale. At the physiological level, the primary impact of the heterogeneities generated by the different scale-down strategies is a drop at the level of the overall biomass yield and this effect has been observed in several studies involving scale-down reactors [12, 14, 15]. Other physiological parameters that have been followed during scaling down and scaling up studies are : side products excretion , cell viability [13, 17–19] and mRNA levels . Magnitude of microbial stress is usually characterized by techniques based on membrane integrity . However, system biology has led a lot of informative material that can be used to find new parameters in order to characterize more precisely microbial stress in bioreactors . These progresses have notably allowed to get more insight about regulation mechanisms involved in the stress response of E. coli. Microbial cells have developed signal transduction systems to sense environmental state and to induce coordinated expression of genes involved in the appropriate stress response, this response being critical for the microorganism to adapt and survive [22, 23]. For bacteria, the coordination of gene expression program is mediated by small proteins, called sigma factors σ, that bind to RNA polymerases (RNAP). Sigma factors increase the affinity of RNAP for specific promoter regions, according to the class of sigma factor. E. coli presents seven classes of sigma factors. The main (or housekeeping) σ factor is involved in the transcription from a majority of the promoters. The six other σ factors, called alternate σ factors, induce the activation of more specific promoters involved in the response to specific environmental stimuli (e.g., heat shock, pH shift,...). Among these alternate factors, σS, coded by rpoS, is the master regulator of the general stress response. System biology approaches have allowed to get more insight at the level of the organization of the transcriptional network of E. coli, showing that rpoS is among the global regulator located at the top level of the hierarchical structure . The gene rpoS controls the expression of more than 50 genes , and its inactivation makes the cells more vulnerable to stress conditions . It is thus interesting to monitor the activity of this gene as a global stress reporter for microbial cells cultivated in heterogeneous reactors. We therefore used a prpoS transcriptional reporter based on the green fluorescent protein expression. Similar construction has been previously used to monitor the effect of various stresses experienced in the environment, such as the osmotic stress [26, 27]. A very recent publication involves a nar::gfp reporter to monitor oxygen availability during E. coli cultivations in lab scale reactors . In this work, we propose to extent the application of GFP sensor to monitor the effect of several bioprocess specific stresses, i.e. mainly glucose excess, limitation and starvation and oxygen exhaustion.
Global trends observed for the different reactor configurations
Culture of E. coli is generally performed in fed-batch reactor with precise control of the substrate addition in order to reach high cell densities and to avoid the excretion of by-products . Three cultivation strategies have been considered in order to alter progressively the mixing efficiency of the reactor (figure 1). At first, a reference culture in a well-mixed reactor with an exponential feed has been performed (figure 1A). In a second time, cyclic perturbations at the level of the sugar concentration and dissolved oxygen level have been generated by aim of an abrupt ON/OFF DO-feed control (figure 1B). In a third time, these perturbations have been enhanced by addition of a tubular section to the original apparatus and by performing the DO-feed control at the level of this section (figure 1C). The partitioned reactors used in this study are based on the scale-down configurations (plug-flow reactor PFR connected to a stirred tank reactor STR) used in previous works [12, 15, 18], except for the DO-controlled feed. In the case of the well mixed bioreactor with exponential feed, hydrodynamics is improved and the resulting extracellular conditions are therefore expected to be the best the growth of the microorganisms. In fact, some limited heterogeneities have to be expected, owing to the pulses generated by the feed pump, but these pulses are rapidly attenuated by the mixing performances of the reactor (typically, about a few seconds for a stirred, small scale, bioreactor). This configuration generates an overall substrate profile slightly above the limitation level (figure 1A). During scale-up, mixing efficiency drops, leading to the development of extracellular heterogeneities. This case study has been reproduced at small scale by using a bioreactor with an ON/OFF dissolved oxygen feed control. Since the regulation trend is not smoothed by a proportional-integral-differential (PID) controller, the glucose is delivered in excess when the dissolved oxygen level rises above 30% from saturation and is subsequently followed by a period of glucose limitation and starvation until the dissolved oxygen level reach again the set point. This situation is depicted at figure 1B, by substrate profile oscillating between the excess and the limitation level. However, the situation is more complex in a large-scale bioreactor, where the extracellular conditions met by the cells do not simply follow an oscillatory trend, and are function of the mixing efficiency of the bioreactor and of the circulation pattern followed by the cells. This latter can be assimilated to a stochastic process  and induces a significant amount of noise at the level of the extracellular environment met by the microorganisms. This case figure has been tentatively reproduced at small-scale by the use of a partitioned scale-down, for which glucose is injected at the level of the tubular part (figure 1C). In this case, cells are exposed to limitation and starvation at the level of the stirred bioreactor and are stochastically pumped through the tubular part where they encounter a strong starvation in the first half of the section and a strong glucose excess in the second half. The effects of the different culture conditions at the level of the general bioprocess trend are depicted at figure 2 for the DO level, and at figure 3 for the microbial growth. It can be clearly seen from figure 3 that there is a significant difference at the level of the biomass evolution between the well-mixed reactor with an exponential feed and the perturbed reactors (with the well-mixed DO-controlled reactor exhibiting a slightly higher biomass yield than the partitioned reactors). These results are in accordance with previous scale-down studies [17, 18]. In the case of the DO-controlled reactor, the oxygen level fluctuations govern the dynamics of the substrate addition (figure 2). For the well mixed-non partitioned reactor, it can be seen that, for the dissolved oxygen profile, frequency of substrate addition is relatively high. In the case of partitioned bioreactors, more spaced dissolved oxygen peaks can be distinguished, revealing a significantly lower substrate feeding frequency. In this case, glucose is punctually added in large quantities, leading to strong concentration gradient inside the partitioned reactor, and particularly throughout the tubular section. In the case of well mixed reactor performed with an exponential feeding strategy, the dissolved oxygen fluctuations can be attributed to the fact that the feeding algorithm is designed in order to operate close to the substrate and dissolved oxygen uptake capacity of the microorganism.
Extracellular fluctuations experienced in heterogeneous bioreactors induce segregation of the rpoS expression level
In addition to the global parameters described in the previous section (biomass and dissolved oxygen evolution), the fate of the microbial population has been monitored by flow cytometry. This technique has been applied for all the bioreactor configurations in order to estimate the expression level associated to a rpoS::GFPmut2 construction inserted in the cells. For each reactor configuration, flow cytometry profiles have been observed at different times during the culture. The whole set of data is too important to be inserted in the text and is thus presented in additional files (additional files 1, 2, 3 and 4, for each reactor configuration). The main information is summarized at figure 4, showing the distribution of rpoS expression level at different times and for each reactor configuration. Segregation at the level of the rpoS expression level can be clearly observed in the case of cultures performed in partitioned reactors. Indeed, two sub-populations of cells can be distinguished, the first sub-population exhibiting a low expression of rpoS and the second one exhibiting a high level of expression associated to the promoter of this gene. The same phenomena can also be observed in a lower extent for the late stages of the culture performed in the well-mixed reactor with an ON/OFF DO-controlled feed. It is thus interesting to note that the intensity as well as the time of appearance of the segregation are both correlated with the mixing efficiency of the reactor. For the less efficient reactor (i.e., the partitioned reactor with the lower recirculation flow rate), the segregation appears early, and two sub-populations are clearly visible after 11 hours of culture. For the partitioned reactor with the higher recirculation flow rate, segregation is only partially observed at 23 hours and begins to be clearly visible at 26 h. For the well-mixed reactor with the DO-feed control, the segregation appears very lately during the culture (26 h) and is fairly attenuated if compared to the previous cases.
The evolution of the GFP expression is depicted at figure 5, considering evolution of the amount of GFP positive cells and the evolution of the total fluorescence. For all the reactor configurations tested, it can be seen that the GFP is extensively synthesized between 5 and 10 hours of culture, this period including the transition from the batch phase to the fed-batch phase. This transition is generally associated with a decrease of the growth rate, a condition leading to the induction of rpoS [30, 31]. In the same context, it has been shown that ppGpp, an alarmone responsible for the induction of rpoS, is strongly synthesized after transition from bath to the fed-batch phase . After this global increase of the GFP content of the cells, fluorescence intensity tends to remain constant for the well-mixed reactor with an exponential feed, and tends to drop in the case of the heterogeneous reactors. It is interesting to note that the drop is significantly higher in the case of the partitioned bioreactors, where the extracellular perturbations experienced by the cells are very intensive. Indeed, after the beginning of the fed-batch regulation, the curves depicted at figure 5 show a progressive decrease proportional to the intensity of the heterogeneities encountered in the corresponding reactor. The drop in the amount of GFP+ cells (or more roughly, the drop of the fluorescence intensity of the microbial population) is enhanced by the segregation mechanism illustrated at figure 4. This effect is distributed among the population since in the worst case, i.e. partitioned bioreactor with the lower recirculation flow rate, there is a minor subpopulation showing a standard level of expression (standard means that in this case, the fluorescence intensities of these cells are in the range of those encountered in the well-mixed reactor). In front of these considerations, it appears that the stochasticity inherent to the extracellular noise (in our case mainly with respect to the glucose and oxygen availability) enhance the segregation phenomena.
However, several artifacts can be linked to GFP expression at the level of the heterogeneous reactors. Indeed, it has been reported that heterogeneities induce several changes at the level of the microorganism physiology (viability, membrane integrity,...) and it is not excluded that these changes have an impact at the level of protein conformation, leading to a fluorescence drop. In order to examine this hypothesis, additional experiments have been carried out with another reporter strain carrying a cyaA::GFPmut2 plasmid. The cyaA codes for adenylate cyclase, a key enzyme promoting the conversion of ATP to cAMP. As shown in figure 6, this enzyme seems to be constitutively expressed, since flow cytometry profiles show a peak of expression of high intensity with a narrow dispersion from cell to cell (the sensitivity of the photomultiplier has been lowered compared with the measurements carried out with the rpoS strain in order to overcome channel saturation as reported in the Methods section). For the experiments carried out with the cyaA reporter strains, the dissolved oxygen profile as well as the microbial growth curve is similar to that obtained in equivalent hydrodynamic conditions with the rpoS reporter strain (data not shown). The most important information depicted at figure 6 is that the evolution of the expression profile is not affected by the bioreactor hydrodynamics (the full datasets are available in additional files 5 and 6, and photographs of pcyaA::gfp cells are available in additional file 7). This result suggests that the rpoS expression level recorded in the previous experiments is not promoted by artifacts linked with the protein maturation and conformation.
The characteristic time for GFP synthesis does not allow to obtain instantaneous informations about the partitioned reactors hydrodynamics
In order to gain some insights about the GFP synthesis process and to assess the reproducibility of the results, the experiment performed in the partitioned reactor with a Qrecirc = 18 L/h has been reiterated. Microbial growth curve and dissolved oxygen profile are similar to those obtained in the previous experiment with the same reactor configuration (data not shown). Samples for flow cytometry have been taken both at the level of the stirred part and at the outlet of the tubular part. The results are shown at figure 7. The segregation phenomenon is observed again at the early stages of the fed-batch process. No significant differences can be observed at the level of the GFP expression profiles between samples collected at the level of the stirred reactor and the tubular part. This observation implies that the characteristic time related to GFP synthesis and maturation is larger than the characteristic time for the exchanges between the mixed and the tubular part of the partitioned bioreactor. In the working conditions used to run the experiment, the theoretical residence time inside the tubular part of the reactor working at a recirculation flow rate of 18 L/h is 3.33 minutes. Informations obtained in the literature indicate that the GFPmut2 variant is a fast folding protein with a synthesis and maturation time of 4 minutes  (the original GFP exhibits a maturation time of approximately 95 minutes ). It is thus normal to observe no difference of expression level between the two parts of the reactor. In this case, the GFP expression level represents the stress accumulated by the cells exposed to extracellular fluctuations. In order to point out physiological parameters giving the instantaneous state of the cells, other molecules with a higher synthesis rate, must be considered. The mRNA transcript analysis constitutes in this case a performant tool, since several techniques, such as RT-PCR , and more recently an at-line technique based on a sandwich hybridization assay , are available. However, mRNA based techniques don't allow to obtain information at the single cell level, but have to be overlooked as a good complement to the reporter gene technology.
The segregation is reversible and seems to have no significant effect on cell viability
At the late stage of the stationary phase of culture depicted at figure 7, it appears, by comparison with the expression profiles obtained in the well-mixed reactor, that the microbial population shifts from a bimodal to a standard unimodal expression profile. This result is confirmed if taking a sample at the level of the partitioned reactor and by putting the collected cells in fresh medium. As shown at figure 8, the population is shifted from a bimodal to an unimodal GFP expression profile after 2 hours of culture in a fresh medium. These results tend to indicate that the mechanisms leading to segregation are reversible and disappear when cells are cultivated without extracellular fluctuations. In the same order of idea, the viability of the cells has been checked at the level of the two subpopulations (GFP+ and GFP-). As shown at figure 9, no difference can be observed at the level of the PI intensity (FL3 channel) between the GFP+ and the GFP- populations. This suggests that cell viability is not correlated with segregation at the level of the stress gene expression. However, and owing to the fact that a lot of mechanisms are involved in rpoS regulation, more precise informations should be obtained if considering the two sub populations separately (e.g., by using cell sorting techniques).
An E. coli reporter strain has been cultivated in different reactor configurations with decreasing level of mixing efficiency. The goal of this work is to observe a physiological parameter that can be easily monitored, i.e. in our case the GFP fluorescence intensity resulting from the activation of the rpoS promoter, in order to put this parameter in relation with the bioreactor mixing efficiency. We have observed a significant decrease of the fluorescence level when mixing conditions were altered.
The following discussions will be mainly focused on possible explanations of the segregation phenomena observed at the level of the prpoS level and to propose some practical applications of the finding pointed out in this work.
When exposed to stress, E. coli is able to respond by the coordinated expression of a set of genes [22, 36]. In normal, non stress, conditions, bacteria have a housekeeping sigma factor (σ70) which governs the main transcriptional machinery. In stress conditions, bacteria have also alternate sigma factors that are able to redirect RNAP to specific promoters. E. coli has six alternate sigma factors (i.e., σS, σH, σN, σE, σF, σFecl) that respond to different cellular stresses . In fluctuating environmental conditions, these alternate sigma factors are in competition with the housekeeping sigma factor. For example, when exposed to nutrient starvation, the σS coded by the rpoS gene and the σ70 coded by the rpoD gene enter in competition [31, 37]. These genomic responses correspond to two distinct strategies developed by the cell to cope with nutrient limitation. In the first case figure, cells react by inducing a set of genes involving the global enhancement for the nutrient uptake in order to lead to an improvement of the growth. This first strategy is usually called the hunger response and is triggered by the housekeeping sigma factor σ70. In the second case figure, cells induce a stringent response involving a complete subset of genes that will lead to a stationary, dormant state. This pathway corresponds to the general stress response triggered by the σS subunit, the product of the rpoS gene that has been tracked throughout this work. When exposed to limitation, there is a constant competition of the two subunits for the core RNA polymerase. In several cases, prolonged starvation leads to the appearance of rpoS mutants lacking in the general stress response. The generally accepted explanation about the appearance of such mutants is an unbalance of the competing process with the prevalence of the σ70 at the level of the core RNA polymerases, resulting in a general reduction of the expression of the promoters associated with the general stress response. Another justification is that rpoS makes E. coli less competitive in hunger conditions, since it decreases the expression of high affinity transport pathway . In such conditions, rpoS mutation can be viewed as a selective advantage and has been defined as a specific phenotype (GASP for Growth Advantage in Stationary Phase) [39, 40].
However, the picture is by far more complex since stress response in bacteria is governed by complex regulatory networks. The rpoS regulation can be performed at different levels: transcription, mRNA turnover and proteolysis . At the transcriptional level, rpoS is induced by ppGpp, an alarmone synthesized under stress conditions. This alarmone is very important for the modulation of starvation response in fed-batch bioreactors [41, 42]. The amount of rpoS transcripts can be regulated by the presence of ppGpp, by but also by mutation at the level of the corresponding synthases (RelA for the synthesis of ppGpp in case of amino acid starvation, and SpoT for the synthesis of ppGpp in case of carbon starvation). However, additional difficulties arise when considering ppGpp as a master regulator of stress response. Indeed, ppGpp modulates competition of σ factors for core RNA polymerase and can thus induce promoters other than σS dependent . On the other hand, ppGpp is not absolutely required for the activity of σS dependent promoters, because underproduction of σ70 (housekeeping) or mutation in rpoD restored the expression of σS dependent promoters . A protein, DksA, acts as a coregulator of genes controlled by ppGpp, and can compensate for the loss of ppGpp [36, 45]. Thus, possible mutations accumulating in the genes involved in this regulation can be unnoticed while studying prpoS activity. In front of these considerations, the phenomena observed in this work can be only understood if considering the entire network or sub-network involved in the general stress response, highlighting the need for integrating system biology into biochemical engineering approaches. Indeed, system biology proposes simplified numerical tools to study the complexity of gene network. As an example, the concept of bistability involves gene networks in order to describe the global behavior of a population of cells, characterized by a bimodal distribution for a physiological characteristic among the population . Bistability is induced by noise and several models involving interactions between subset of genes (called network motifs) have been proposed to explain its occurrence [47–49]. More specifically, bistability arises when a master regulator of the gene network is involved. In our case, we have investigated the level of expression of the rpoS gene, which is a recognized master regulator of the E. coli genes network belonging to the top level of the hierarchical structure of this network . It is also recognized that bistability (or in our case segregation phenomena) is strongly enhanced when the stochasticity of the extracellular fluctuations is increased . This observation is in accordance with our results, showing that the most intensive segregation phenomena has been observed in the partitioned bioreactors for which the extracellular perturbations are strongly governed by a stochastic mechanism at the level of the passage between the stirred part and the tubular part .
A possible explanation of the segregation phenomena can be postulated on the basis of the properties of GFP. As said before, rpoS regulation occurs at three levels: transcription, mRNA stability and proteolysis. Without excluding the numerous phenomena occurring at the transcriptional level, and considering the properties of GFPmut2, proteolysis can be involved in the fluorescence drop associated with scale-down effect. Indeed, GFPmut2 is a fast folding mutant of the original protein with a high stability, and it is generally assumed that the lifetime of this protein is longer than the cell cycle . In our case, segregation phenomena, i.e. strong fluorescence decreases for a major part of the microorganisms (about 80%), occurs within one hour (see flow cytometry profiles in additional files), suggesting that fluorescence drop cannot be only attributed to GFP dilution due to cell division. The possible explanation involves proteolysis of the GFP content of the cells. This phenomenon has not been observed in the case of the pcyaA::gfp reporter strain cultivated in the same conditions. However, in this case, the GFP content of cells was significantly higher than for prpoS strain, and proteolysis can be less marked.
The results acquired during this work suggest a strong modulation of the rpoS promoter in front of the intensity, but also the frequency at which perturbations occur. Two kinds of scale-down reactor have been used throughout the study. The first one is the DO-controlled reactor inducing oscillating extracellular conditions according to the glucose and dissolved oxygen level. In this kind of reactor, all the cells are exposed alternatively to the oscillating conditions. The second one is the partitioned bioreactor which is different from the previous one. In this kind of reactor, cells are exposed to oscillating environmental conditions, but a fraction of cells can be exposed to stochastic fluctuations when crossing the tubular section. It seems that the stochastic component plays an important role by enhancing the intensity of the segregation mechanism among the population. Surprisingly, sub-optimal conditions induce a drop of prpoS activity. This observation can be explained by the induction of stress proteases that has been proposed previously as the possible mechanism leading to the diminution of GFP content of the cells. However, the segregation mechanism cannot be entirely attributed to the stochastic nature of the fluctuations, i.e. to the physical fragmentation of the reacting volume into a well-mixed and a heterogeneous zone, because the segregation mechanism has also been observed in the DO-controlled reactor. Lin et al.  have reported complex mechanisms having strong influence on protein synthesis by E. coli in function of the frequencies of the extracellular segregation. In this case, additional informations can be collected if using a method with a smaller response time than the reporter gene technology. Indeed, it has been reported that the characteristic time related to GFP is too large, even if a fast folding mutant is used and if the characteristic time associated to the exchange process in the partitioned reactors is important, to observe the instantaneous stress level of a cell in a particular flow region. No difference of GFP level has been observed between the mixed and the nonmixed part of the partitioned reactor. This fact highlights the potentialities of mRNA based techniques as a complementary tool for GFP expression studies.
Considering the different rpoS expression levels obtained in function of the hydrodynamic efficiency of the reactor, this gene can be considered as a good reporter of the physiological stress experienced by E. coli in process conditions. In practise, some specific probes are available for the on-line monitoring of the GFP fusion proteins and could be used to make some useful control of the bioreactor in direct relation with the physiological status of the cells (the initial concept of fluorescence probe that can be fit in bioreactors has been developed by Randers-Eichhorn and co-workers , and this concept has been more recently improved ). The main limitation with this approach lies on the fact that these fluorescence probes only give a signal averaging the physiological status over the whole population. However, it has been shown that there is a strong segregation of the physiological status of the cells in the heterogeneous reactors, leading to a subpopulation with a strong expression level and another one with a low expression level. In a bioreactor control perspective, it will be thus very useful to adapt an on-line flow cytometer. The efficiency of such automated apparatus have been estimated in the literature but is not yet routinely used for the on-line monitoring of bioreactors .
Previous scale-down experiments involve the observation of direct (mRNA analysis, cellular viability) or indirect (global growth rate, oxygen consumption rate,...) physiological parameters, but it is the first time that a GFP reporter strain is used in this context. The use of a prpoS::GFPmut2 reporter strain shows an increasing level of segregation when the mixing efficiency of the bioreactor is altered. Experiments carried out in different perturbed reactor (i.e., where glucose and dissolved oxygen fluctuations occur) show that a fraction of the cellular population lose its stress response capabilities in heterogeneous conditions. It appears thus that the intensity, as well as the frequency at which extracellular perturbations occur, are of importance in the repartition of gene expression level among the population. This imply that prpoS::gfp reporter is thus a promising tool for monitoring stress level encountered by microbial cells in process conditions.
Strain and cultivation media
E. coli K12 MG1655 bearing a pMS201 (4260 bp) plasmid with either the prpoS::GFPmut2 or the pcyaA::GFPmut2 transcriptional reporter and a kanamycin resistance gene. These strains comes from a cloning vector library elaborated at the Weizmann Institute of Science . The strain is maintained at -80°C in working seeds vials (2 mL) in solution with LB media and with 40% of glycerol. The precultures and cultures have been performed on a defined liquid medium containing (in g/L): K2HPO4 14.6, NaH2PO4.2H2O 3.6; Na2SO4 2; (NH4)2SO4 2.47, NH4Cl 0.5, (NH4)2-H-citrate 1, glucose 5, thiamine 0.01, kanamycin 0.1. Thiamine and kanamycin are sterilized by filtration (0.2 μm). The medium is supplemented with 3 mL/L of trace solution, 3 mL/L of a FeCl3.6H2O solution (16.7 g/L), 3 mL/L of an EDTA solution (20.1 g/L) and 2 mL/L of a MgSO4 solution (120 g/L). The trace solution contains (in g/L) : CoCl2.H2O 0.74, ZnSO4.7H2O 0.18, MnSO4.H2O 0.1, CuSO4.5H2O, CoSO4.7H2O. The precultures are performed in 600 mL of the above mentioned medium in shake flask at 37°C. The bioreactor experiments are performed initially with 10 L of the above mentioned medium. The feed solution contains 2 L of the basal medium but with all the concentrations doubled, except for the glucose which is at 500 g/L. During the culture, cell growth has been monitored by optical density at a wavelength of 600 nm. Cell dry weight has been determined on the basis of filtered samples (0.45 μm) dried during 24 hours at 105°C.
Three kinds of reactor configuration have been used and are schematized at figure 1. The stirred vessel is a standard bioreactor (vessel diameter: 0.22 m; initial working volume: 10 L) equipped with four equally spaced baffles and a single rushton disk turbine with 6 blades. Agitation rate is maintained at 500 min-1 throughout the culture. Air is injected through a ring sparger. Air flow rate is maintained at 10 L/min during the batch phase and is increased to 40 L/min during the fed-batch phase. During the fed-batch phase the reactor internal pressure is increased to 0.4 bars in order to improve oxygen solubility. In the case of the partitioned reactors, stirred reactor is connected to a tubular section having a diameter of 0.05 m and a length of 1 m. The tubular section contains 14 static mixing elements (Kenics) in order to keep a plug-flow hydrodynamics. Recirculation between the stirred vessel and the tubular section is ensured by a peristaltic pump (Watson Marlow series 323). For each assay, temperature is maintained at 37°C and at pH 7 by a digital control system (ABB). The pH is maintained constant by the addition of ammonia solution (25%). Dissolved oxygen is measured by a polarographic probe (Mettler Toledo). In the case of the exponential feed, the equation has been calculated and adjusted based on a previous fed-batch experiment. The equation has been implemented in a MatLab program driving the feed pump (Watson Marlow series 101 U/R) working by pulse on the basis of a flow rate of 15 mL/min. In the case of the DO-controlled feed, the pump activation is ensured by the ON/OFF contact option at the level of the DO transmitter (Knick oxy 2402) with a set point of 30% from saturation. In the case of well-mixed reactors, glucose is added by a special port at the top of the vessel. In the case of the partitioned reactors, glucose is added at the mid length of the tubular part, inducing the appearance of a glucose limited zone above the injection port and a glucose excess zone below the injection port.
Flow cytometry analysis
The analysis of the GFP expression level has been performed with a FACscan (Becton Dickinson) flow cytometer. Samples are taken directly from the reactor and are diluted in 900 μL of PBS and 100 μL of a cycloheximide solution (1 mg/mL) in order to stop protein synthesis. For each measurement, 30,000 cells are analyzed (GFP is excited at 488 nm and emission signals are collected by using filters at 530 nm). The measurements are repeated 3 times at different FL1 channel intensities (FL1 at 480, 550 and 620 for the rpoS strain and FL1 at 400, 480 and 550 for the cyaA strain). The results have been analyzed by the CellQuest (Becton Dickinson) software and are subsequently exported to WinMDI and MatLab for further analysis. The GFP negative fraction of cells has been determined on the basis of the initial FL1 distribution for each reactor experiment. The initial FL1 distribution is very comparable from a reactor condition to another and corresponds to a very low GFP emission (monitored by epifluorescence microscopy). Cell viability estimation is carried out by using propidium iodide (PI) at a working concentration of 5 μg/ml. Before addition of PI, cells are washed with PBS. Sample is then divided in two parts, the first being untreated and the second being stained with PI for 5 minutes at room temperature.
Hewitt CJ, Nienow AW: The scale-up of microbial batch and fed-batch fermentation processes. Advances in applied microbiology. 2007, 62: 105-135.
Deckwer WD, Jahn D, Hempel D, Zeng AP: Systems biology approaches to bioprocess development. Engineering in life sciences. 2006, 6 (5): 455-469.
Lara AR, Galindo E, Ramirez OT, Palomares LA: Living with heterogeneities in bioreactors – Understanding the effects of environmental gradients on cells. Molecular biotechnology. 2006, 34: 355-381.
Delvigne F, Destain J, Thonart P: A methodology for the design of scale-down bioreactors by the use of mixing and circulation stochastic models. Biochemical engineering journal. 2006, 28 (3): 256-268.
Lapin A, Schmid J, Reuss M: Modeling the dynamics of E. coli populations in the three-dimensional turbulent field of a stirred bioreactor – A structured-segregated approach. Chemical engineering science. 2006, 61: 4783-4797.
Delvigne F, Lejeune A, Destain J, Thonart P: Stochastic models to study the impact of mixing on a fed-batch culture of Saccharomyces cerevisiae. Biotechnology progress. 2006, 22: 259-269.
Lapin A, Müller D, Reuss M: Dynamic behavior of microbial populations in stirred bioreactors simulated with Euler-Lagrange methods : traveling along the lifelines of single cells. Industrial and engineering chemistry research. 2004, 43: 4647-4656.
Fowler JD, Dunlop EH: Effects of reactant heterogeneity and mixing on catabolite repression in cultures of Saccharomyces cerevisiae. Biotechnology and bioengineering. 1989, 33: 1039-1046.
Oosterhuis NMG, Kossen NWF, Olivier APC, Schenk ES: Scale-down and optimization studies of the gluconic acid fermentation by Gluconobacter oxydans. Biotechnology and bioengineering. 1985, 27: 711-720.
Lara AR, Leal L, Flores N, Gosset G, Bolivar F, Ramirez OT: Transcriptional and metabolic response of recombinant Escherichia coli to spatial dissolved oxygen tension gradients simulated in a scale-down system. Biotechnology and bioengineering. 2005, 93 (2): 372-385.
Sandoval-Basurto EA, Gosset G, Bolivar F, Ramirez OT: Culture of Escherichia coli under dissolved oxygen gradients simulated in a two-compartment scale-down system : metabolic response and production of recombinant protein. Biotechnology and bioengineering. 2004, 89 (4): 453-463.
Neubauer P, Häggström L, Enfors SO: Influence of substrate oscillations on acetate formation and growth yield in Escherichia coli glucose limited fed-batch cultivations. Biotechnology and bioengineering. 1995, 47: 139-146.
Hewitt CJ, Onyeaka H, Lewis G, Taylor IW, Nienow AW: A Comparison of High Cell Density Fed-Batch Fermentations Involving Both Induced and Non-Induced Recombinant Escherichia coli Under Well-Mixed Small-Scale and Simulated Poorly Mixed Large-Scale Conditions. Biotechnology and bioengineering. 2007, 96 (3): 495-505.
Delvigne F, Destain J, Thonart P: Towards a stochastic formulation of the microbial growth in relation with the bioreactor performances : case study of an E. coli fed-batch process. Biotechnology progress. 2006, 22 (4): 1114-1124.
Enfors SO, Jahic M, Rozkov A, Xu B, Hecker M, Jürgen B, Krüger E, Schweder T, Hamer G, O'Beirne D, Noisommit-Rizzi N, Reuss M, Boone L, Hewitt C, McFarlane C, Nienow A, Kovacs T, Trägardh C, Fuchs L, Revstedt J, Friberg PC, Hjertager B, Blomsten G, Skogman H, Hjort S, Hoeks F, Lin HY, Neubauer P, Lans van der R, Luyben K, Vrabel P, Manelius A: Physiological responses to mixing in large scale bioreactors. Journal of biotechnology. 2001, 85: 175-185.
Xu B, Jahic M, Blomsten G, Enfors SO: Glucose overflow metabolism and mixed-acid fermentation in aerobic large-scale fed-batch processes with Escherichia coli. Applied microbiology and biotechnology. 1999, 51: 564-571.
Hewitt CJ, Nebe-Von Caron G, Nienow AW, Mc Farlane CM: Use of multi-staining flow cytometry to characterise the physiological state of Escherichia coli W3110 in high cell density fed-batch cultures. Biotechnology and bioengineering. 1999, 63 (6): 705-711.
Hewitt CJ, Nebe-Von Caron G, Axelsson B, Mc Farlane CM, Nienow AW: Studies related to the scale-up of high-cell-density E. coli fed-batch fermentations using multiparameter flow cytometry : effect of a changing microenvironment with respect to glucose and dissolved oxygen concentration. Biotechnology and bioengineering. 2000, 70 (4): 381-390.
Sundstrom H, Wallberg F, Ledung E, Norrman B, Hewitt CJ, Enfors SO: Segregation to non-dividing cells in recombinant Escherichia coli fed-batch fermentation processes. Biotechnology letters. 2004, 26: 1533-1539.
Schweder T, Krüger E, Xu B, Jürgen B, Blomsten G, Enfors SO, Hecker M: Monitoring of genes that respond to process related stress in large-scale bioprocesses. Biotechnology and bioengineering. 1999, 65 (2): 151-159.
Booth IR: Stress and the single cell : intrapopulation diversity is a mechanism to ensure survival upon exposure to stress. International journal of food microbiology. 2002, 78: 19-30.
Abee T, Wouters JA: Microbial stress response in minimal processing. International journal of food microbiology. 1999, 50: 65-91.
Chung HJ, Bang W, Drake MA: Stress response of Escherichia coli. Comprehensive reviews in food science and food safety. 2006, 5: 52-64.
Ma HW, Buer J, Zeng AP: Hierarchical structure and modules in the Escherichia coli transcriptional regulatory network revealed by a new top-down approach. BMC bioinformatics. 2004, 5: 199-
Brown L, Gentry D, Elliott T, Cashel M: DksA Affects ppGpp Induction of RpoS at a Translational Level. Journal of bacteriology. 2002, 184 (16): 4455-4465.
Funabashi H, Mie M, Yanagida Y, Kobatake E, Aizawa M: Fluorescent monitoring of cellular physiological status depending on the accumulation of ppGpp. Biotechnology letters. 2002, 24: 269-273.
Funabashi H, Haruyama T, Mie M, Yanagida Y, Kobatake E, Aizawa M: Non-destructive monitoring of rpoS promoter activity as stress marker for evaluating cellular physiological status. Journal of biotechnology. 2002, 95 (1): 85-93.
Garcia JR, Cha HJ, Rao GG, Marten MR, Bentley WE: Microbial nar-GFP cell sensors reveal oxygen limitations in highly agitated and aerated laboratory-scale fermentors. Microbial cell factories. 2009, 8: 6-
Shiloach J, Fass R: Growing E. coli to high cell density – A historical perspective on method development. Biotechnology advances. 2005, 23: 345-357.
Notley L, Ferenci T: Induction of RpoS-dependent functions in glucose-limited continuous culture : what level of nutrient limitation induces the stationary phase of Escherichia coli. Journal of bacteriology. 1996, 178 (5): 1465-1468.
Notley-McRobb L, King T, Ferenci T: rpoS mutations and loss of general stress resistance in Escherichia coli populations as a consequence of conflict between competing stress responses. Journal of bacteriology. 2002, 184 (3): 806-811.
Andersson L, Yang S, Neubauer P, Enfors SO: Impact of plasmid presence and induction on cellular responses in fed batch cultures of Escherichia coli. Journal of biotechnology. 1996, 46: 255-263.
Zaslaver A, Bren A, Ronen M, Itzkovitz S, Kikoin I, Shavit S, Liebermeister W, Surette MG, Alon U: A comprehensive library of fluorescent transcriptional reporters for Escherichia coli. Nature methods. 2006, 3 (8): 623-628.
DeLisa MP, Li J, Rao G, Weigand WA, Bentley WE: Monitoring GFP operon fusion protein expression during high cell density cultivation of Escherichia coli using an on-line optical sensor. Biotechnology and bioengineering. 1999, 65 (1): 54-64.
Pioch D, Jürgen B, Evers S, Maurer KH, Hecker M, Schweder T: At-line monitoring of bioprocess-relevant marker genes. Engineering in life sciences. 2007, 7 (4): 373-379.
Costanzo A, Nicoloff H, Barchinger SE, Banta AB, Gourse RL, Ades SE: ppGpp and DksA likely regulate the activity of the extracytoplasmic stress factor sigma E in Escherichia coli by both direct and indirect mechanisms. Molecular microbiology. 2008, 67 (3): 61-632.
Aertsen A, Michiels CW: Stress and how bacteria cope with death and survival. Critical reviews in microbiology. 2004, 30: 263-273.
Ferenci T: Hungry bacteria – definition and properties of a nutritional state. Environmental microbiology. 2001, 3 (10): 605-611.
Fukuda T, Nakahigashi K, Inokuchi H: Viability of Escherichia coli cells under long-term cultivation in a rich nutrient medium. Genes and genetic systems. 2001, 76: 271-278.
Zinser ER, Kolter R: Escherichia coli evolution during stationary phase. Research in microbiology. 2004, 155: 328-336.
Teich A, Meyer S, Lin HY, Andersson L, Enfors SO, Neubauer P: Growth rate related concentration changes of the starvation response regulators sigma S and ppGpp in glucose-limited fed-batch and continuous cultures of Escherichia coli. Biotechnology progress. 1999, 15: 123-129.
Neubauer P, Ahman M, Törnkvist M, Larsson G, Enfors SO: Response of guanosine tetraphosphate to glucose fluctuations in fed-batch cultivations of Escherichia coli. Journal of biotechnology. 1995, 43: 195-204.
Laurie AD, Bernardo LMD, Sze CC, Skärfstad E, Szalewska-Palasz A, Nyström T, Shingler V: The role of the alarmone (p)ppGpp in sigma N competition for core RNA polymerase. The journal of biological chemistry. 2003, 278 (3): 1494-1503.
Jishage M, Kvint K, Shingler V, Nyström T: Regulation of sigma factor competition by the alarmone ppGpp. Genes and development. 2002, 16: 1260-1270.
Magnusson LU, Gummesson B, Joksimovic P, Farewell A, Nyström T: Identical, independent, and opposing roles of ppGpp and DksA in Escherichia coli. Journal of bacteriology. 2007, 189 (14): 5193-5202.
Dubnau D, Losick R: Bistability in bacteria. Molecular microbiology. 2006, 61 (3): 564-572.
Alon U: Network motifs : theory and experimental approaches. Nature. 2007, 8: 450-461.
Süel GM, Garcia-Ojalvo J, Liberman LM, Elowitz MB: An excitable gene regulatory circuit induces transient cellular differentiation. Nature. 2006, 440 (23): 545-550.
Shen-Orr SS, Milo R, Manga S, Alon U: Network motifs in the transcriptional regulation network of Escherichia coli. Nature genetics. 2002, 31: 64-68.
Kussell E, Leibler S: Phenotypic diversity, population growth and information in fluctuating environments. Science. 2005, 309: 2075-2078.
Ou J, Yamada T, Nagahisa K, Hirassawa T, Furusawa C, Yomo T, Shimizu H: Dynamic change in promoter activation during lysine biosynthesis in Escherichia coli cells. Molecular biosystems. 2008, 4: 128-134.
Lin HY, Neubauer P: Influence of controlled glucose oscillations on a fed-batch process of recombinant Escherichia coli. Journal of biotechnology. 2000, 79: 27-37.
Randers-Eichhorn L, Albano CR, Sipior J, Bentley WE, Rao G: On-line green fluorescence protein sensor with LED excitation. Biotechnology and bioengineering. 1997, 55 (6): 921-926.
Jones JJ, Bridges AM, Fosberry AP, Gardner S, Lowers RR, Newby RR, James PJ, Hall RM, Jenkins O: Potential of real-time measurement of GFP-fusion proteins. Journal of biotechnology. 2004, 109: 201-211.
Patkar A, Vijayasankaran N, Urry DW, Srienc F: Flow cytometry as a useful tool for process development : rapid evaluation of expression systems. Journal of biotechnology. 2002, 93: 217-229.
Frank Delvigne and Mathieu Boxus are both post-doctoral researchers at the Belgian "Fond de la Recherche Scientifique" (FNRS) and gratefully acknowledge the FRNS for financial support. The authors also acknowledge the CGRI that has supported this work by a European scientific exchange program (PAI Tournesol) and for the resulting fruitful discussions with the peoples of the LISBP from the INSA Toulouse (Dr. Nathalie Gorret, Dr. Stéphane Guillouet and Prof. Carole Molina-Jouve). A special thanks to Prof. R. Kettman for giving us the opportunity to work in his lab (unité de biologie cellulaire et moléculaire, FUSAGx).
The authors declare that they have no competing interests.
FD planned the experiments, performed the first half of the experiments carried out in bioreactors and wrote the final version of the manuscript. MB performed the flow cytometry experiments. SI performed the second half of the experiments carried out in bioreactors. PT supervised the work. All authors read and approved the final version of the manuscript.