Specific growth rate dependent transcriptome profiling of Escherichia coli K12 MG1655 in accelerostat cultures

Type I fimbriae subunit fimA enhances Escherichia coli biofilm formation but affects L-threonine carbon distribution

The biofilm (BF) provides favorable growth conditions to cells, which has been exploited in the field of industrial biotechnology. Based on our previous research works on type I fimbriae for the biosynthesis of L-threonine (LT) in Escherichia coli, in this study, a fimA-overexpressing strain was engineered, which improved BF formation under industrial fermentation conditions. The morphological observation and characterization of BF formation were conducted to verify the function of the subunit FimA. However, it was not suitable for repeated-batch immobilized fermentation as the LT titer was not elevated significantly. The underlying molecular mechanisms of BF formation and the LT carbon flux were explored by transcriptomic analysis. The results showed that fimA regulated E. coli BF formation but affected LT carbon distribution. This study will stimulate thoughts about how the fimbriae gene regulated biofilms and amino acid excretion and will bring some consideration and provide a reference for the development of BF-based biomanufacturing processes in E. coli.

Genome-Scale Reconstruction of Microbial Dynamic Phenotype: Successes and Challenges

This review is a part of the SI 'Genome-Scale Modeling of Microorganisms in the Real World'. The goal of GEM is the accurate prediction of the phenotype from its respective genotype under specified environmental conditions. This review focuses on the dynamic phenotype; prediction of the real-life behaviors of microorganisms, such as cell proliferation, dormancy, and mortality; balanced and unbalanced growth; steady-state and transient processes; primary and secondary metabolism; stress responses; etc. Constraint-based metabolic reconstructions were successfully started two decades ago as FBA, followed by more advanced models, but this review starts from the earlier nongenomic predecessors to show that some GEMs inherited the outdated biokinetic frameworks compromising their performances. The most essential deficiencies are: (i) an inadequate account of environmental conditions, such as various degrees of nutrients limitation and other factors shaping phenotypes; (ii) a failure to simulate the adaptive changes of MMCC (MacroMolecular Cell Composition) in response to the fluctuating environment; (iii) the misinterpretation of the SGR (Specific Growth Rate) as either a fixed constant parameter of the model or independent factor affecting the conditional expression of macromolecules; (iv) neglecting stress resistance as an important objective function; and (v) inefficient experimental verification of GEM against simple growth (constant MMCC and SGR) data. Finally, we propose several ways to improve GEMs, such as replacing the outdated Monod equation with the SCM (Synthetic Chemostat Model) that establishes the quantitative relationships between primary and secondary metabolism, growth rate and stress resistance, process kinetics, and cell composition.

Inactivation of the Pta-AckA pathway impairs fitness of Bacillus anthracis during overflow metabolism

Under conditions of glucose excess, aerobically growing bacteria predominantly direct carbon flux towards acetate fermentation, a phenomenon known as overflow metabolism or the bacterial 'Crabtree effect'. Numerous studies of the major acetate-generating pathway, the Pta-AckA, revealed its important role in bacterial fitness through the control of central metabolism to sustain balanced growth and cellular homeostasis. In this work, we highlight the contribution of the Pta-AckA pathway to fitness of the spore-forming bacterium, Bacillus anthracis We demonstrate that disruption of the Pta-AckA pathway causes a drastic growth reduction in the mutants and alters the metabolic and energy status of the cells. Our results revealed that inactivation of the Pta-AckA pathway increases the glucose consumption rate, affects intracellular ATP, NAD⁺ and NADH levels and leads to a metabolic block at the pyruvate and acetyl-CoA nodes. Consequently, accumulation of intracellular acetyl-CoA and pyruvate forces bacteria to direct carbon into the TCA and/or glyoxylate cycles as well as fatty acid and poly(3-hydroxybutyrate) (PHB) biosynthesis pathways. Notably, the presence of phosphate butyryltransferase in B. anthracis partially compensates for the loss of phosphotransacetylase activity. Furthermore, overexpression of the ptb gene not only eliminates the negative impact of the pta mutation on B. anthracis fitness, but also restores normal growth in the pta mutant of the non-butyrate-producing bacterium, Staphylococcus aureus Taken together, the results of this study demonstrate the importance of the Pta-AckA pathway for B. anthracis fitness by revealing its critical contribution to the maintenance of metabolic homeostasis during aerobic growth under conditions of carbon overflow.IMPORTANCE B. anthracis, the etiologic agent of anthrax, is a highly pathogenic, spore-forming bacterium that causes acute, life-threatening disease in both humans and livestock. A greater understanding of the metabolic determinants governing fitness of B. anthracis is essential for the development of successful therapeutic and vaccination strategies aimed at lessening the potential impact of this important biodefense pathogen. This study is the first to demonstrate the vital role of the Pta-AckA pathway in preserving energy and metabolic homeostasis in B. anthracis under conditions of carbon overflow, therefore, highlighting this pathway as a potential therapeutic target for drug discovery. Overall, the results of this study provide important insight into understanding the metabolic processes and requirements driving rapid B. anthracis proliferation during vegetative growth.

Correlated chromosomal periodicities according to the growth rate and gene expression

Linking genetic information to population fitness is crucial to understanding living organisms. Despite the abundant knowledge of the genetic contribution to growth, the overall patterns/features connecting genes, their expression, and growth remain unclear. To reveal the quantitative and direct connections, systematic growth assays of single-gene knockout Escherichia coli strains under both rich and poor nutritional conditions were performed; subsequently, the resultant growth rates were associated with the original expression levels of the knockout genes in the parental genome. Comparative analysis of growth and the transcriptome identified not only the nutritionally differentiated fitness cost genes but also a significant correlation between the growth rates of the single-gene knockout strains and the original expression levels of these knockout genes in the parental strain, regardless of the nutritional variation. In addition, the coordinated chromosomal periodicities of the wild-type transcriptome and the growth rates of the strains lacking the corresponding genes were observed. The common six-period periodicity was somehow attributed to the essential genes, although the underlying mechanism remains to be addressed. The correlated chromosomal periodicities associated with the gene expression-growth dataset were highly valuable for bacterial growth prediction and discovering the working principles governing minimal genetic information.

Continuous Cultivation as a Tool Toward the Rational Bioprocess Development With Pichia Pastoris Cell Factory

The methylotrophic yeast Pichia pastoris (Komagataella phaffii) is currently considered one of the most promising hosts for recombinant protein production (RPP) and metabolites due to the availability of several tools to efficiently regulate the recombinant expression, its ability to perform eukaryotic post-translational modifications and to secrete the product in the extracellular media. The challenge of improving the bioprocess efficiency can be faced from two main approaches: the strain engineering, which includes enhancements in the recombinant expression regulation as well as overcoming potential cell capacity bottlenecks; and the bioprocess engineering, focused on the development of rational-based efficient operational strategies. Understanding the effect of strain and operational improvements in bioprocess efficiency requires to attain a robust knowledge about the metabolic and physiological changes triggered into the cells. For this purpose, a number of studies have revealed chemostat cultures to provide a robust tool for accurate, reliable, and reproducible bioprocess characterization. It should involve the determination of key specific rates, productivities, and yields for different C and N sources, as well as optimizing media formulation and operating conditions. Furthermore, studies along the different levels of systems biology are usually performed also in chemostat cultures. Transcriptomic, proteomic and metabolic flux analysis, using different techniques like differential target gene expression, protein description and ¹³C-based metabolic flux analysis, are widely described as valued examples in the literature. In this scenario, the main advantage of a continuous operation relies on the quality of the homogeneous samples obtained under steady-state conditions, where both the metabolic and physiological status of the cells remain unaltered in an all-encompassing picture of the cell environment. This contribution aims to provide the state of the art of the different approaches that allow the design of rational strain and bioprocess engineering improvements in Pichia pastoris toward optimizing bioprocesses based on the results obtained in chemostat cultures. Interestingly, continuous cultivation is also currently emerging as an alternative operational mode in industrial biotechnology for implementing continuous process operations.

The highly conserved chromosomal periodicity of transcriptomes and the correlation of its amplitude with the growth rate in Escherichia coli

The growth rate, representing the fitness of a bacterial population, is determined by the transcriptome. Chromosomal periodicity, which is known as the periodic spatial pattern of a preferred chromosomal distance in microbial genomes, is a representative overall feature of the transcriptome; however, whether and how it is associated with the bacterial growth rate are unknown. To address these questions, we analysed a total of 213 transcriptomes of multiple Escherichia coli strains growing in an assortment of culture conditions varying in terms of temperature, nutrition level and osmotic pressure. Intriguingly, Fourier transform analyses of the transcriptome identified a common chromosomal periodicity of transcriptomes, which was independent of the variation in genomes and environments. In addition, fitting of the data to a theoretical model, we found that the amplitudes of the periodic transcriptomes were significantly correlated with the growth rates. These results indicated that the amplitude of periodic transcriptomes is a parameter representing the global pattern of gene expression in correlation with the bacterial growth rate. Thus, our study provides a novel parameter for evaluating the adaptiveness of a growing bacterial population and quantitatively predicting the growth dynamics according to the global expression pattern.

Amino acids are key substrates to Escherichia coli BW25113 for achieving high specific growth rate

Studying substrate consumption in nutrient-rich conditions is challenging because often the growth medium includes undefined components like yeast extract or peptone. For clear and consistent results, it is necessary to use defined medium, where substrate utilization can be followed. In the present work, Escherichia coli BW25113 batch growth in a medium supplemented with 20 proteinogenic amino acids and glucose was studied. Focus was on the quantitative differences in substrate consumption and proteome composition between minimal and nutrient-rich medium. In the latter, 72% of carbon used for biomass growth came from amino acids and 28% from glucose. Serine was identified as the most consumed substrate with 41% of total carbon consumption. Proteome comparison between nutrient-rich and minimal medium revealed changes in TCA cycle and acetate producing enzymes that together with extracellular metabolite data pointed to serine being consumed mainly for energy generation purposes. Serine removal from the growth medium decreased specific growth rate by 22%. In addition, proteome comparison between media revealed a large shift in amino acid synthesis and translation related proteins. Overall, this work describes in quantitative terms the batch growth carbon uptake profile and proteome allocation of E. coli BW25113 in minimal and nutrient-rich medium.

Quantitative Flow Cytometry to Understand Population Heterogeneity in Response to Changes in Substrate Availability in Escherichia coli and Saccharomyces cerevisiae Chemostats

Microbial cells in bioprocesses are usually described with averaged parameters. But in fact, single cells within populations vary greatly in characteristics such as stress resistance, especially in response to carbon source gradients. Our aim was to introduce tools to quantify population heterogeneity in bioprocesses using a combination of reporter strains, flow cytometry, and easily comprehensible parameters. We calculated mean, mode, peak width, and coefficient of variance to describe distribution characteristics and temporal shifts in fluorescence intensity. The skewness and the slope of cumulative distribution function plots illustrated differences in distribution shape. These parameters are person-independent and precise. We demonstrated this by quantifying growth-related population heterogeneity of Saccharomyces cerevisiae and Escherichia coli reporter strains in steady-state of aerobic glucose-limited chemostat cultures at different dilution rates and in response to glucose pulses. Generally, slow-growing cells showed stronger responses to glucose excess than fast-growing cells. Cell robustness, measured as membrane integrity after exposure to freeze-thaw treatment, of fast-growing cells was strongly affected in subpopulations of low membrane robustness. Glucose pulses protected subpopulations of fast-growing but not slower-growing yeast cells against membrane damage. Our parameters could successfully describe population heterogeneity, thereby revealing physiological characteristics that might have been overlooked during traditional averaged analysis.

Efficient Biofilm-Based Fermentation Strategies for L-Threonine Production by Escherichia coli

Biofilms provide cells favorable growth conditions, which have been exploited in industrial biotechnological processes. However, industrial application of the biofilm has not yet been reported in Escherichia coli, one of the most important platform strains, though the biofilm has been extensively studied for pathogenic reasons. Here, we engineered E. coli by overexpressing the fimH gene, which successfully enhanced its biofilm formation under industrial aerobic cultivation conditions. Subsequently, a biofilm-based immobilized fermentation strategy was developed. L-threonine production was increased from 10.5 to 14.1 g/L during batch fermentations and further to 17.5 g/L during continuous (repeated-batch) fermentations with enhanced productivities. Molecular basis for the enhanced biofilm formation and L-threonine biosynthesis was also studied by transcriptome analysis. This study goes beyond the conventional research focusing on pathogenic aspects of E. coli biofilm and represents a successful application case of engineered E. coli biofilm to industrial processes.

Using Bacterial Transcriptomics to Investigate Targets of Host-Bacterial Interactions in Caenorhabditis elegans

The interactions between a host and its resident microbes form complicated networks that can affect host physiology. Disentangling these host-microbe interactions can help us better understand mechanisms by which bacteria affect hosts, while also defining the integral commensal protection that host-associated microbiota offer to promote health. Here we utilize a tractable genetic model organism, Caenorhabditis elegans, to study the effects of host environments on bacterial gene expression and metabolic pathways. First, we compared the transcriptomic profiles of E. coli OP50 in vitro (on agar plates) versus in vivo (fed to C. elegans host). Our data revealed that 110 biosynthetic genes were enriched in host-associated E. coli. Several of these expressed genes code for the precursors and products needed for the synthesis of lipopolysaccharides (LPS), which are important for innate immune and stress responses, as well as pathogenicity. Secondly, we compared the transcriptomic profiles of E. coli fed to hosts with different genetic backgrounds, including the long-lived daf-2/insulin like growth factor (IGF) receptor and short lived daf-16/FOXO transcription factor mutants. We find that hosts genetics also alters bacterial metabolic pathways. Given that bacteria influence host health, this transcriptomics approach can elucidate genes mediating host aging.

A novel nicotinamide adenine dinucleotide control strategy for increasing the cell density of Haemophilus parasuis

Haemophilus parasuis is the causative agent of Glässer's disease and is a major source of economic losses in the swine industry each year. To enhance the production of an inactivated vaccine against H. parasuis, the availability of nicotinamide adenine dinucleotide (NAD) must be carefully controlled to ensure a sufficiently high cell density of H. parasuis. In the present study, the real-time viable cell density of H. parasuis was calculated based on the capacitance of the culture. By assessing the relationship between capacitance and viable cell density/NAD concentration, the NAD supply rate could be adjusted in real time to maintain the NAD concentration at a set value based on the linear relationship between capacitance and NAD consumption. The linear relationship between cell density and addition of NAD indicated that 7.138 × 10⁹ NAD molecules were required to satisfy per cell growth. Five types of NAD supply strategy were used to maintain different NAD concentration for H. parasuis cultivation, and the results revealed that the highest viable cell density (8.57, OD₆₀₀ ) and cell count (1.57 × 10¹⁰ CFU/mL) were obtained with strategy III (NAD concentration maintained at 30 mg/L), which were 1.46- and 1.45- times more, respectively, than cultures with using NAD supply strategy I (NAD concentration maintained at 10 mg/L). An extremely high cell density of H. parasuis was achieved using this NAD supply strategy, and the results demonstrated a convenient and reliable method for determining the real-time viable cell density relative to NAD concentration. Moreover, this method provides a theoretical foundation and an efficient approach for high cell density cultivation of other auxotroph bacteria.

Metabolic engineering for improving l-tryptophan production in Escherichia coli

l-Tryptophan is an important aromatic amino acid that is used widely in the food, chemical, and pharmaceutical industries. Compared with the traditional synthetic methods, production of l-tryptophan by microbes is environmentally friendly and has low production costs, and feed stocks are renewable. With the development of metabolic engineering, highly efficient production of l-tryptophan in Escherichia coli has been achieved by eliminating negative regulation factors, improving the intracellular level of precursors, engineering of transport systems and overexpression of rate-limiting enzymes. However, challenges remain for l-tryptophan biosynthesis to be cost-competitive. In this review, successful and applicable strategies derived from metabolic engineering for increasing l-tryptophan accumulation in E. coli are summarized. In addition, perspectives for further efficient production of l-tryptophan are discussed.

Metabolic modeling and response surface analysis of an Escherichia coli strain engineered for shikimic acid production

Background

Classic metabolic engineering strategies often induce significant flux imbalances to microbial metabolism, causing undesirable outcomes such as suboptimal conversion of substrates to products. Several mathematical frameworks have been developed to understand the physiological and metabolic state of production strains and to identify genetic modification targets for improved bioproduct formation. In this work, a modeling approach was applied to describe the physiological behavior and the metabolic fluxes of a shikimic acid overproducing Escherichia coli strain lacking the major glucose transport system, grown on complex media.

Results

The obtained flux distributions indicate the presence of high fluxes through the pentose phosphate and Entner-Doudoroff pathways, which could limit the availability of erythrose-4-phosphate for shikimic acid production even with high flux redirection through the pentose phosphate pathway. In addition, highly active glyoxylate shunt fluxes and a pyruvate/acetate cycle are indicators of overflow glycolytic metabolism in the tested conditions. The analysis of the combined physiological and flux response surfaces, enabled zone allocation for different physiological outputs within variant substrate conditions. This information was then used for an improved fed-batch process designed to preserve the metabolic conditions that were found to enhance shikimic acid productivity. This resulted in a 40% increase in the shikimic acid titer (60 g/L) and 70% increase in volumetric productivity (2.45 gSA/L*h), while preserving yields, compared to the batch process.

Conclusions

The combination of dynamic metabolic modeling and experimental parameter response surfaces was a successful approach to understand and predict the behavior of a shikimic acid producing strain under variable substrate concentrations. Response surfaces were useful for allocating different physiological behavior zones with different preferential product outcomes. Both model sets provided information that could be applied to enhance shikimic acid production on an engineered shikimic acid overproducing Escherichia coli strain.

Electronic supplementary material

The online version of this article (10.1186/s12918-018-0632-4) contains supplementary material, which is available to authorized users.

Avoiding amino acid depletion in a complex medium results in improved Escherichia coli BW25113 growth

We studied Escherichia coli BW25113 growth in a complex medium with emphasis on amino acid consumption. The aim was to profile amino acid utilization in acid-hydrolysed casein and a defined nutrient-rich medium and based on these measurements modify the medium for better growth performance. Amino acid depletions in both media caused apparent biomass growth stops that prolonged growth duration. Obtained amino acid consumption values enabled a new defined medium to be formulated, where no growth stops were observed, the specific growth rate was constant, and the provided substrates were fully utilized. Similarly, we modified the acid-hydrolysed casein medium by adding pure amino acids that removed the apparent biomass growth stops. Key to our results was the combination of growth medium analysis and process monitoring data, specifically oxygen partial pressure and produced carbon dioxide that were used to track growth changes. Our findings showed the deficiencies of the nutrient-rich medium and how rational medium design, based on consumption values, removed these shortcomings. The resulting balanced medium gives a high specific growth rate and is suitable for studying E. coli physiology at fast growth.

Gene Expression Order Attributed to Genome Reduction and the Steady Cellular State in Escherichia coli

Transcriptomes not only reflect the growth status but also link to the genome in bacteria. To investigate if and how genome or cellular state changes contribute to the gene expression order, the growth profile-associated transcriptomes of an assortment of genetically differentiated Escherichia coli either exponentially growing under varied conditions or in response to environmental disturbance were analyzed. A total of 168 microarray data sets representing 56 transcriptome variations, were categorized by genome size (full length or reduced) and cellular state (steady or unsteady). At the genome-wide level, the power-law distribution of gene expression was found to be significantly disturbed by the genome size but not the cellular state. At the regulatory network level, more networks with improved coordination of growth rates were observed in genome reduction than at the steady state. At the single-gene level, both genome reduction and steady state increased the correlation of gene expression to growth rate, but the enriched gene categories with improved correlations were different. These findings not only illustrate the order of gene expression attributed to genome reduction and steady cellular state but also indicate that the accessory sequences acquired during genome evolution largely participated in the coordination of transcriptomes to growth fitness.

Influence of Glucose Availability and CRP Acetylation on the Genome-Wide Transcriptional Response of Escherichia coli: Assessment by an Optimized Factorial Microarray Analysis

Background: While in eukaryotes acetylation/deacetylation regulation exerts multiple pleiotropic effects, in Escherichia coli it seems to be more limited and less known. Hence, we aimed to progress in the characterization of this regulation by dealing with three convergent aspects: the effector enzymes involved, the master regulator CRP, and the dependence on glucose availability.

Methods: The transcriptional response of E. coli BW25113 was analyzed across 14 relevant scenarios. These conditions arise when the wild type and four isogenic mutants (defective in deacetylase CobB, defective in N(ε)-lysine acetyl transferase PatZ, Q- and R-type mutants of protein CRP) are studied under three levels of glucose availability (glucose-limited chemostat and glucose-excess or glucose-exhausted in batch culture). The Q-type emulates a permanent stage of CRP^acetylated, whereas the R-type emulates a permanent stage of CRP^deacetylated. The data were analyzed by an optimized factorial microarray method (Q-GDEMAR).

Results: (a) By analyzing one mutant against the other, we were able to unravel the true genes that participate in the interaction between ΔcobB/ΔpatZ mutations and glucose availability; (b) Increasing stages of glucose limitation appear to be associated with the up-regulation of specific sets of target genes rather than with the loss of genes present when glucose is in excess; (c) Both CRP^deacetylated and CRP^acetylated produce extensive changes in specific subsets of genes, but their number and identity depend on the glucose availability; (d) In other sub-sets of genes, the transcriptional effect of CRP seems to be independent of its acetylation or deacetylation; (e) Some specific ontology functions can be associated with each of the different sets of genes detected herein.

Conclusions: CRP cannot be thought of only as an effector of catabolite repression, because it acts along all the glucose conditions tested (excess, limited, and exhausted), exerting both positive and negative effects through different sets of genes. Acetylation of CRP does not seem to be a binary form of regulation, as there is not a univocal relationship between its activation/inhibitory effect and its acetylation/deacetylation stage. All the combinatorial possibilities are observed. During the exponential growth phase, CRP also exerts a very significant transcriptional effect, mainly on flagellar assembly and chemotaxis (FDR = 7.2 × 10⁻⁴⁴).

Detection of growth rate-dependent product formation in miniaturized parallel fed-batch cultivations

Saccharomyces cerevisiae is a popular expression system for recombinant proteins. In most cases, production processes are performed as carbon-limited fed-batch cultures to avoid aerobic ethanol formation. Especially for constitutive expression systems, the specific product formation rate depends on the specific growth rate. The development of optimal feeding strategies strongly depends on laboratory-scale cultivations, which are time and resource consuming, especially when continuous experiments are carried out. It is therefore beneficial for accelerated process development to look at alternatives. In this study, S. cerevisiae AH22 secreting a heterologous endo-polygalacturonase (EPG) was characterized in microwell plates with an enzyme-based fed-batch medium. Through variation of the glucose release rate, different growth profiles were established and the impact on EPG secretion was analyzed. Product formation rates of 200-400 U (gx h)-1 were determined. As a reference, bioreactor experiments using the change-stat cultivation technique were performed. The growth-dependent product formation was analyzed over dilution rates of D = 0.01-0.35 with smooth change of D at a rate of 0.003 h-2. EPG production was found to be comparable with a qp of 400 U (gx h)-1 at D = 0.27 h-1. The presented results indicate that parallel miniaturized fed-batch cultures can be applied to determine product formation profiles of putative production strains. With further automation and parallelization of the concept, strain characterization can be performed in shorter time.

Precise, High-throughput Analysis of Bacterial Growth

Bacterial growth is a central concept in the development of modern microbial physiology, as well as in the investigation of cellular dynamics at the systems level. Recent studies have reported correlations between bacterial growth and genome-wide events, such as genome reduction and transcriptome reorganization. Correctly analyzing bacterial growth is crucial for understanding the growth-dependent coordination of gene functions and cellular components. Accordingly, the precise quantitative evaluation of bacterial growth in a high-throughput manner is required. Emerging technological developments offer new experimental tools that allow updates of the methods used for studying bacterial growth. The protocol introduced here employs a microplate reader with a highly optimized experimental procedure for the reproducible and precise evaluation of bacterial growth. This protocol was used to evaluate the growth of several previously described Escherichia coli strains. The main steps of the protocol are as follows: the preparation of a large number of cell stocks in small vials for repeated tests with reproducible results, the use of 96-well plates for high-throughput growth evaluation, and the manual calculation of two major parameters (i.e., maximal growth rate and population density) representing the growth dynamics. In comparison to the traditional colony-forming unit (CFU) assay, which counts the cells that are cultured in glass tubes over time on agar plates, the present method is more efficient and provides more detailed temporal records of growth changes, but has a stricter detection limit at low population densities. In summary, the described method is advantageous for the precise and reproducible high-throughput analysis of bacterial growth, which can be used to draw conceptual conclusions or to make theoretical observations.

Why Do Fast-Growing Bacteria Enter Overflow Metabolism? Testing the Membrane Real Estate Hypothesis

Bacteria and other cells show a puzzling behavior. At high growth rates, E. coli switch from respiration (which is ATP-efficient) to using fermentation for additional ATP (which is inefficient). This overflow metabolism results in a several-fold decrease in ATP produced per glucose molecule provided as food. By integrating diverse types of experimental data into a simple biophysical model, we give evidence that this onset is the result of the membrane real estate hypothesis: Fast growth drives cells to be bigger, reducing their surface-to-volume ratios. This decreases the membrane area available for respiratory proteins despite growing demand, causing increased crowding. Only when respiratory proteins reach their crowding limit does the cell activate fermentation, since fermentation allows faster ATP production per unit membrane area. Surface limitation thus creates a Pareto trade-off between membrane efficiency and ATP yield that links metabolic choice to the size and shape of a bacterial cell. By exploring the predictions that emerge from this trade-off, we show how consideration of molecular structures, energetics, rates, and equilibria can provide important insight into cellular behavior.

Coordinated Changes in Mutation and Growth Rates Induced by Genome Reduction

Genome size is determined during evolution, but it can also be altered by genetic engineering in laboratories. The systematic characterization of reduced genomes provides valuable insights into the cellular properties that are quantitatively described by the global parameters related to the dynamics of growth and mutation. In the present study, we analyzed a small collection of W3110 Escherichia coli derivatives containing either the wild-type genome or reduced genomes of various lengths to examine whether the mutation rate, a global parameter representing genomic plasticity, was affected by genome reduction. We found that the mutation rates of these cells increased with genome reduction. The correlation between genome length and mutation rate, which has been reported for the evolution of bacteria, was also identified, intriguingly, for genome reduction. Gene function enrichment analysis indicated that the deletion of many of the genes encoding membrane and transport proteins play a role in the mutation rate changes mediated by genome reduction. Furthermore, the increase in the mutation rate with genome reduction was highly associated with a decrease in the growth rate in a nutrition-dependent manner; thus, poorer media showed a larger change that was of higher significance. This negative correlation was strongly supported by experimental evidence that the serial transfer of the reduced genome improved the growth rate and reduced the mutation rate to a large extent. Taken together, the global parameters corresponding to the genome, growth, and mutation showed a coordinated relationship, which might be an essential working principle for balancing the cellular dynamics appropriate to the environment.

Genome reduction is a powerful approach for investigating the fundamental rules for living systems. Whether genetically disturbed genomes have any specific properties that are different from or similar to those of natively evolved genomes has been under investigation. In the present study, we found that Escherichia coli cells with reduced genomes showed accelerated nucleotide substitution errors (mutation rates), although these cells retained the normal DNA mismatch repair systems. Intriguingly, this finding of correlation between reduced genome size and a higher mutation rate was consistent with the reported evolution of mutation rates. Furthermore, the increased mutation rate was quantitatively associated with a decreased growth rate, indicating that the global parameters related to the genome, growth, and mutation, which represent the amount of genetic information, the efficiency of propagation, and the fidelity of replication, respectively, are dynamically coordinated.

Correlation between genome reduction and bacterial growth

Genome reduction by removing dispensable genomic sequences in bacteria is commonly used in both fundamental and applied studies to determine the minimal genetic requirements for a living system or to develop highly efficient bioreactors. Nevertheless, whether and how the accumulative loss of dispensable genomic sequences disturbs bacterial growth remains unclear. To investigate the relationship between genome reduction and growth, a series of Escherichia coli strains carrying genomes reduced in a stepwise manner were used. Intensive growth analyses revealed that the accumulation of multiple genomic deletions caused decreases in the exponential growth rate and the saturated cell density in a deletion-length-dependent manner as well as gradual changes in the patterns of growth dynamics, regardless of the growth media. Accordingly, a perspective growth model linking genome evolution to genome engineering was proposed. This study provides the first demonstration of a quantitative connection between genomic sequence and bacterial growth, indicating that growth rate is potentially associated with dispensable genomic sequences.

The Escherichia coli transcriptome linked to growth fitness

A series of Escherichia coli strains with varied genomic sequences were subjected to high-density microarray analyses to elucidate the fitness-correlated transcriptomes. Fitness, which is commonly evaluated by the growth rate during the exponential phase, is not only determined by the genome but is also linked to growth conditions, e.g., temperature. We previously reported genetic and environmental contributions to E. coli transcriptomes and evolutionary transcriptome changes in thermal adaptation. Here, we describe experimental details on how to prepare microarray samples that truly represent the growth fitness of the E. coli cells. A step-by-step record of sample preparation procedures that correspond to growing cells and transcriptome data sets that are deposited at the GEO database (GSE33212, GSE52770, GSE61739) are also provided for reference.

Proteome reallocation in Escherichia coli with increasing specific growth rate

Cells usually respond to changing growth conditions with a change in the specific growth rate (μ) and adjustment of their proteome to adapt and maintain metabolic efficiency. Description of the principles behind proteome resource allocation is important for understanding metabolic regulation in response to changing μ. Thus, we analysed the proteome resource allocation dynamics of Escherichia coli into different metabolic processes in response to changing μ. E. coli was grown on minimal and defined rich media in steady state continuous cultures at different μ and characterised combining two LC-MS/MS-based proteomics methods: stable isotope labelling by amino acids in cell culture (SILAC) and intensity based label-free absolute quantification. We detected slowly growing cells investing more proteome resources in energy generation and carbohydrate transport and metabolism whereas for achieving faster growth cells needed to devote most resources to translation and processes closely related to the protein synthesis pipeline. Furthermore, down-regulation of energy generation and carbohydrate metabolism proteins with faster growth displayed very similar expression dynamics with the global transcriptional regulator CRP (cyclic AMP receptor protein), pointing to a dominant protein resource allocating role of this protein. Our data also suggest that acetate overflow may be the result of global proteome resource optimisation as cells saved proteome resources by switching from fully respiratory to respiro-fermentative growth. The presented results give a quantitative overview of how E. coli adjusts its proteome to achieve faster growth and in future could contribute to the design of more efficient cell factories through proteome optimisation.

Bacterial transcriptome reorganization in thermal adaptive evolution

Background

Evolution optimizes a living system at both the genome and transcriptome levels. Few studies have investigated transcriptome evolution, whereas many studies have explored genome evolution in experimentally evolved cells. However, a comprehensive understanding of evolutionary mechanisms requires knowledge of how evolution shapes gene expression. Here, we analyzed Escherichia coli strains acquired during long-term thermal adaptive evolution.

Results

Evolved and ancestor Escherichia coli cells were exponentially grown under normal and high temperatures for subsequent transcriptome analysis. We found that both the ancestor and evolved cells had comparable magnitudes of transcriptional change in response to heat shock, although the evolutionary progression of their expression patterns during exponential growth was different at either normal or high temperatures. We also identified inverse transcriptional changes that were mediated by differences in growth temperatures and genotypes, as well as negative epistasis between genotype—and heat shock-induced transcriptional changes. Principal component analysis revealed that transcriptome evolution neither approached the responsive state at the high temperature nor returned to the steady state at the regular temperature. We propose that the molecular mechanisms of thermal adaptive evolution involve the optimization of steady-state transcriptomes at high temperatures without disturbing the heat shock response.

Conclusions

Our results suggest that transcriptome evolution works to maintain steady-state gene expression during constrained differentiation at various evolutionary stages, while also maintaining responsiveness to environmental stimuli and transcriptome homeostasis.

Electronic supplementary material

The online version of this article (doi:10.1186/s12864-015-1999-x) contains supplementary material, which is available to authorized users.

Effect of Global Regulators RpoS and Cyclic-AMP/CRP on the Catabolome and Transcriptome of Escherichia coli K12 during Carbon- and Energy-Limited Growth

For heterotrophic microbes, limited availability of carbon and energy sources is one of the major nutritional factors restricting the rate of growth in most ecosystems. Physiological adaptation to this hunger state requires metabolic versatility which usually involves expression of a wide range of different catabolic pathways and of high-affinity carbon transporters; together, this allows for simultaneous utilization of mixtures of carbonaceous compounds at low concentrations. In Escherichia coli the stationary phase sigma factor RpoS and the signal molecule cAMP are the major players in the regulation of transcription under such conditions; however, their interaction is still not fully understood. Therefore, during growth of E. coli in carbon-limited chemostat culture at different dilution rates, the transcriptomes, expression of periplasmic proteins and catabolomes of strains lacking one of these global regulators, either rpoS or adenylate cyclase (cya), were compared to those of the wild-type strain. The inability to synthesize cAMP exerted a strong negative influence on the expression of alternative carbon source uptake and degradation systems. In contrast, absence of RpoS increased the transcription of genes belonging to high-affinity uptake systems and central metabolism, presumably due to reduced competition of σ^D with σ^S. Phenotypical analysis confirmed this observation: The ability to respire alternative carbon substrates and to express periplasmic high-affinity binding proteins was eliminated in cya and crp mutants, while these properties were not affected in the rpoS mutant. As expected, transcription of numerous stress defence genes was negatively affected by the rpoS knock-out mutation. Interestingly, several genes of the RpoS stress response regulon were also down-regulated in the cAMP-negative strain indicating a coordinated global regulation. The results demonstrate that cAMP is crucial for catabolic flexibility during slow, carbon-limited growth, whereas RpoS is primarily involved in the regulation of stress response systems necessary for the survival of this bacterium under hunger conditions.

Global coordination in adaptation to gene rewiring

Gene rewiring is a common evolutionary phenomenon in nature that may lead to extinction for living organisms. Recent studies on synthetic biology demonstrate that cells can survive genetic rewiring. This survival (adaptation) is often linked to the stochastic expression of rewired genes with random transcriptional changes. However, the probability of adaptation and the underlying common principles are not clear. We performed a systematic survey of an assortment of gene-rewired Escherichia coli strains to address these questions. Three different cell fates, designated good survivors, poor survivors and failures, were observed when the strains starved. Large fluctuations in the expression of the rewired gene were commonly observed with increasing cell size, but these changes were insufficient for adaptation. Cooperative reorganizations in the corresponding operon and genome-wide gene expression largely contributed to the final success. Transcriptome reorganizations that generally showed high-dimensional dynamic changes were restricted within a one-dimensional trajectory for adaptation to gene rewiring, indicating a general path directed toward cellular plasticity for a successful cell fate. This finding of global coordination supports a mechanism of stochastic adaptation and provides novel insights into the design and application of complex genetic or metabolic networks.

Subpopulation-proteomics reveal growth rate, but not cell cycling, as a major impact on protein composition in Pseudomonas putida KT2440

Population heterogeneity occurring in industrial microbial bioprocesses is regarded as a putative effector causing performance loss in large scale. While the existence of subpopulations is a commonly accepted fact, their appearance and impact on process performance still remains rather unclear. During cell cycling, distinct subpopulations differing in cell division state and DNA content appear which contribute individually to the efficiency of the bioprocess. To identify stressed or impaired subpopulations, we analyzed the interplay of growth rate, cell cycle and phenotypic profile of subpopulations by using flow cytometry and cell sorting in conjunction with mass spectrometry based global proteomics. Adjusting distinct growth rates in chemostats with the model strain Pseudomonas putida KT2440, cells were differentiated by DNA content reflecting different cell cycle stages. The proteome of separated subpopulations at given growth rates was found to be highly similar, while different growth rates caused major changes of the protein inventory with respect to e.g. carbon storage, motility, lipid metabolism and the translational machinery. In conclusion, cells in various cell cycle stages at the same growth rate were found to have similar to identical proteome profiles showing no significant population heterogeneity on the proteome level. In contrast, the growth rate clearly determines the protein composition and therefore the metabolic strategy of the cells.

High-level production of L-threonine by recombinant Escherichia coli with combined feeding strategies

The process of L-threonine production using Escherichia coli TRFC was investigated, and the result showed that there was a large amount of acetic acid in the broth. The effects of acetic acid, which is a known inhibitory metabolite in E. coli cultivation, on L-threonine production by recombinant E. coli TRFC were evaluated, and the result indicated that the growth of E. coli TRFC and L-threonine formation were significantly inhibited in the presence of acetic acid. Two combined feeding strategies were applied to L-threonine fed-batch fermentation in order to investigate the effects of the feeding strategy on L-threonine fermentation. The results showed that using the combined feeding strategy of pseudo-exponential feeding and glucose-stat feeding resulted in high cell density (36.67 g L⁻¹) and L-threonine production (124.57 g L⁻¹) as well as low accumulation of by-products. This work provides a useful approach for large-scale production of L-threonine.

Co-Expression of Threonine Dehydratase and Acetolactate Synthase in Escherichia Coli for L-Isoleucine Production

Metabolic engineering ofCorynebacterium glutamicumhas sought to divert carbon into L-isoleucine. However, the fermentation period of this strain is long. TheC.glutamicumYILW strain (LeuL, AHVr, SGr, Leu-MEr) was previously derived by repeated compound mutagenesis which could accumulate 20.2 g/L L-isoleucine in a 5-L jar fermentor. Overexpression of the threonine dehydratase gene (ilvA) fromCorynebacterium glutamicumYILW and coexpression of threonine dehydratase and acetolactate synthase (ilvBN) from it were employed to divert carbon flux toward L-isoleucine. The strainE. coliTRFC with the expression ofilvA could accumulate L-isoleucine of 6.8 g/L without accumulation of any L-threonine by fed-batch fermentation in a 5-L jar fermentor. However, the production of L-isoleucine by the strainE.coliTRFC with the co-expression ofilvA andilvBN was decreased by 19.1%, and the production of L-valine was increased by 40% compared with that ofE. coliTRFC with the expression ofilvA.

A central role for carbon-overflow pathways in the modulation of bacterial cell death

Similar to developmental programs in eukaryotes, the death of a subpopulation of cells is thought to benefit bacterial biofilm development. However mechanisms that mediate a tight control over cell death are not clearly understood at the population level. Here we reveal that CidR dependent pyruvate oxidase (CidC) and α-acetolactate synthase/decarboxylase (AlsSD) overflow metabolic pathways, which are active during staphylococcal biofilm development, modulate cell death to achieve optimal biofilm biomass. Whereas acetate derived from CidC activity potentiates cell death in cells by a mechanism dependent on intracellular acidification and respiratory inhibition, AlsSD activity effectively counters CidC action by diverting carbon flux towards neutral rather than acidic byproducts and consuming intracellular protons in the process. Furthermore, the physiological features that accompany metabolic activation of cell death bears remarkable similarities to hallmarks of eukaryotic programmed cell death, including the generation of reactive oxygen species and DNA damage. Finally, we demonstrate that the metabolic modulation of cell death not only affects biofilm development but also biofilm-dependent disease outcomes. Given the ubiquity of such carbon overflow pathways in diverse bacterial species, we propose that the metabolic control of cell death may be a fundamental feature of prokaryotic development.

Many bacterial species including the pathogen Staphylococcus aureus are capable of adhering to surfaces and forming complex communities called biofilms. This mode of growth can be particularly challenging from an infection control standpoint, as they are often refractory to antibiotics and host immune system. Although developmental processes underlying biofilm formation are not entirely clear, recent evidence suggests that cell death of a subpopulation is crucial for its maturation. In this study we provide insight regarding the metabolic pathways that control cell death and demonstrate that acetate, a by-product of glucose catabolism, potentiates a form of cell death that exhibits physiological and biochemical hallmarks of apoptosis in eukaryotic organisms. Finally, we demonstrate that altering the ability of metabolic pathways that regulate acetate mediated cell death in S. aureus affects the outcome of biofilm-related diseases, such as infective endocarditis.

The metabolic potential of Escherichia coli BL21 in defined and rich medium

BACKGROUND

The proteome reflects the available cellular machinery to deal with nutrients and environmental challenges. The most common E. coli strain BL21 growing in different, commonly employed media was evaluated using a detailed quantitative proteome analysis.

RESULTS

The presence of preformed biomass precursor molecules in rich media such as Luria Bertani supported rapid growth concomitant to acetate formation and apparently unbalanced abundances of central metabolic pathway enzymes, e.g. high levels of lower glycolytic pathway enzymes as well as pyruvate dehydrogenase, and low levels of TCA cycle and high levels of the acetate forming enzymes Pta and AckA. The proteome of cells growing exponentially in glucose-supplemented mineral salt medium was dominated by enzymes of amino acid synthesis pathways, contained more balanced abundances of central metabolic pathway enzymes, and a lower portion of ribosomal and other translational proteins. Entry into stationary phase led to a reconstruction of the bacterial proteome by increasing e.g. the portion of proteins required for scavenging rare nutrients and general cell protection. This proteomic reconstruction during entry into stationary phase was more noticeable in cells growing in rich medium as they have a greater reservoir of recyclable proteins from the translational machinery.

CONCLUSIONS

The proteomic comparison of cells growing exponentially in different media reflected the antagonistic and competitive regulation of central metabolic pathways through the global transcriptional regulators Cra, Crp, and ArcA. For example, the proteome of cells growing exponentially in rich medium was consistent with a dominating role of phosphorylated ArcA most likely a result from limitations in reoxidizing reduced quinones in the respiratory chain under these growth conditions. The proteomic alterations of exponentially growing cells into stationary phase cells were consistent with stringent-like and stationary phase responses and a dominating control through DksA-ppGpp and RpoS.

Escherichia coli achieves faster growth by increasing catalytic and translation rates of proteins

Regulation levels of the gene expression cascade controlling protein levels and metabolic fluxes for cells to achieve faster growth have not been elaborated in acceptable detail. Furthermore, there is need for specific growth rate (μ) dependent absolute quantitative transcriptome and proteome data to understand the molecular relationships for enabling cells to modify μ. We address these questions, for the first time, by presenting regulatory strategies for more efficient metabolism of Escherichia coli at higher μ by statistical covariance analysis of genome-wide intracellular mRNA and protein concentrations coupled to metabolic flux analysis in the steady state range of μ = 0.11-0.49 h(-1). Our analyses show dominating post-transcriptional control of protein abundances and post-translational control of flux rates. On average, E. coli achieved five-times faster growth through 3.7-fold increase of apparent catalytic rates of enzymes (kapp) and 2.5-fold increased translation rates, demonstrating the relevance of post-translational regulation for increasing flux throughput. Interestingly, pathways carrying the highest flux showed both high protein abundance and kapp values. Furthermore, co-regulation analysis of enzymatic capacities revealed tightly coupled regulatory dependencies of protein synthesis and RNA precursor synthesis, substrate utilization, biosynthetic and energy generation pathways carrying the highest flux. We also observed metabolic pathway and COG specific protein and metabolic flux control levels, protein expression costs and genome-wide principles for translation efficiency and transcription unit polarity. This work contributes to the much needed quantitative understanding of coordinated gene expression regulation and metabolic flux control. Our findings will also advance modeling and metabolic engineering of industrial strains.

Dual role of transcription and transcript stability in the regulation of gene expression in Escherichia coli cells cultured on glucose at different growth rates

Microorganisms extensively reorganize gene expression to adjust growth rate to changes in growth conditions. At the genomic scale, we measured the contribution of both transcription and transcript stability to regulating messenger RNA (mRNA) concentration in Escherichia coli. Transcriptional control was the dominant regulatory process. Between growth rates of 0.10 and 0.63 h(-1), there was a generic increase in the bulk mRNA transcription. However, many transcripts became less stable and the median mRNA half-life decreased from 4.2 to 2.8 min. This is the first evidence that mRNA turnover is slower at extremely low-growth rates. The destabilization of many, but not all, transcripts at high-growth rate correlated with transcriptional upregulation of genes encoding the mRNA degradation machinery. We identified five classes of growth-rate regulation ranging from mainly transcriptional to mainly degradational. In general, differential stability within polycistronic messages encoded by operons does not appear to be affected by growth rate. We show here that the substantial reorganization of gene expression involving downregulation of tricarboxylic acid cycle genes and acetyl-CoA synthetase at high-growth rates is controlled mainly by transcript stability. Overall, our results demonstrate that the control of transcript stability has an important role in fine-tuning mRNA concentration during changes in growth rate.

Dynamic transcriptional response of Escherichia coli to inclusion body formation

Escherichia coli is used intensively for recombinant protein production, but one key challenge with recombinant E. coli is the tendency of recombinant proteins to misfold and aggregate into insoluble inclusion bodies (IBs). IBs contain high concentrations of inactive recombinant protein that require recovery steps to salvage a functional recombinant protein. Currently, no universally effective method exists to prevent IB formation in recombinant E. coli. In this study, DNA microarrays were used to compare the E. coli gene expression response dynamics to soluble and insoluble recombinant protein production. As expected and previously reported, the classical heat-shock genes had increased expression due to IB formation, including protein folding chaperones and proteases. Gene expression levels for protein synthesis-related and energy-synthesis pathways were also increased. Many transmembrane transporter and corresponding catabolic pathways genes had decreased expression for substrates not present in the culture medium. Additionally, putative genes represented over one-third of the genes identified to have significant expression changes due to IB formation, indicating many important cellular responses to IB formation still need to be characterized. Interestingly, cells grown in 3% ethanol had significantly reduced gene expression responses due to IB formation. Taken together, these results indicate that IB formation is complex, stimulates the heat-shock response, increases protein and energy synthesis needs, and streamlines transport and catabolic processes, while ethanol diminished all of these responses.

Growth rate-coordinated transcriptome reorganization in bacteria

Background

Cell growth rate reflects an organism’s physiological state and largely relies on the ability of gene expression to respond to the environment. The relationship between cellular growth rate and gene expression remains unknown.

Results

Growth rate-coordinated changes in gene expression were discovered by analyzing exponentially growing Escherichia coli cells cultured under multiple defined environments, in which osmotic pressure, temperature and starvation status were varied. Gene expression analyses showed that all 3,740 genes in the genome could be simply divided into three clusters (C1, C2 and C3), which were accompanied by a generic trend in the growth rate that was coordinated with transcriptional changes. The direction of transcriptional change in C1 indicated environmental specificity, whereas those in C2 and C3 were correlated negatively and positively with growth rates, respectively. The three clusters exhibited differentiated gene functions and gene regulation task division.

Conclusions

We identified three gene clusters, exhibiting differential gene functions and distinct directions in their correlations with growth rates. Reverses in the direction of the growth rate correlated transcriptional changes and the distinguished duties of the three clusters indicated how transcriptome homeostasis is maintained to balance the total expression cost for sustaining life in new habitats.

Genome-scale models of metabolism and gene expression extend and refine growth phenotype prediction

Growth is a fundamental process of life. Growth requirements are well-characterized experimentally for many microbes; however, we lack a unified model for cellular growth. Such a model must be predictive of events at the molecular scale and capable of explaining the high-level behavior of the cell as a whole. Here, we construct an ME-Model for Escherichia coli--a genome-scale model that seamlessly integrates metabolic and gene product expression pathways. The model computes ~80% of the functional proteome (by mass), which is used by the cell to support growth under a given condition. Metabolism and gene expression are interdependent processes that affect and constrain each other. We formalize these constraints and apply the principle of growth optimization to enable the accurate prediction of multi-scale phenotypes, ranging from coarse-grained (growth rate, nutrient uptake, by-product secretion) to fine-grained (metabolic fluxes, gene expression levels). Our results unify many existing principles developed to describe bacterial growth.

Genetic engineering of Escherichia coli to enhance production of L-tryptophan

Reducing the accumulation of acetate in Escherichia coli cultures can decrease carbon efflux as by-products and reduce acetate toxicity, and therefore enable high cell density cultivation. The concentration of intracellular amino acids can be decreased by genetic modifications of the corresponding amino acid transport systems. This can increase the levels of amino acids in the fermentation broth by decreasing the feedback inhibition on the corresponding biosynthetic pathways. Here, the effects of genetic manipulation of phosphate acetyltransferase (pta), high affinity tryptophan transporter (mtr) and aromatic amino acid exporter (yddG) on L-tryptophan production were investigated. The pta mutants accumulated less acetate and showed higher capacity for producing L-tryptophan as compared with the parental strain. The strains lacking mtr, or overexpressed yddG, or with the both mtr knockout and yddG overexpression, accumulated lower concentrations of intracellular tryptophan but higher production of extracellular L-tryptophan. In the L-tryptohan fed-batch fermentation of an E. coli derived from TRTH0709/pMEL03 having deletion of pta-mtr and overexpression of yddG in a 30-L fermentor, the maximum concentration of L-tryptophan (48.68 g/L) was obtained, which represented a 15.96 % increase as compared with the parental strain. Acetate accumulated to a concentration of 0.95 g/L. The intracellular concentration of L-tryptophan was low, and the glucose conversion rate reached a high level of 21.87 %, which was increased by 15.53 % as compared with the parent strain.

Strategy for pH control and pH feedback-controlled substrate feeding for high-level production of L-tryptophan by Escherichia coli

Optimum production of L-tryptophan by Escherichia coli depends on pH. Here, we established conditions for optimizing the production of L-tryptophan. The optimum pH range was 6.5-7.2, and pH was controlled using a three-stage strategy [pH 6.5 (0-12 h), pH 6.8 (12-24 h), and pH 7.2 (24-38 h)]. Specifically, ammonium hydroxide was used to adjust pH during the initial 24 h, and potassium hydroxide and ammonium hydroxide (1:2, v/v) were used to adjust pH during 24-38 h. Under these conditions, NH4 (+) and K(+) concentrations were kept below the threshold for inhibiting L-tryptophan production. Optimization was also accomplished using ratios (v/v) of glucose to alkali solutions equal to 4:1 (5-24 h) and 6:1 (24-38 h). The concentration of glucose and the pH were controlled by adjusting the pH automatically. Applying a pH-feedback feeding method, the steady-state concentration of glucose was maintained at approximately 0.2 ± 0.02 g/l, and acetic acid accumulated to a concentration of 1.15 ± 0.03 g/l, and the plasmid stability was 98 ± 0.5 %. The final, optimized concentration of L-tryptophan was 43.65 ± 0.29 g/l from 52.43 ± 0.38 g/l dry cell weight.

Comparison and applications of label-free absolute proteome quantification methods on Escherichia coli

Three different label-free proteome quantification methods--APEX, emPAI and iBAQ--were evaluated to measure proteome-wide protein concentrations in the cell. All the methods were applied to a sample from Escherichia coli chemostat culture. A Pearson squared correlation of approximately 0.6 among the three quantification methods was demonstrated. Importantly, the sum of quantified proteins by iBAQ and emPAI corresponded with the Lowry total protein quantification, demonstrating applicability of label-free methods for an accurate calculation of protein concentrations at the proteome level. The iBAQ method showed the best correlation between biological replicates, a normal distribution among all protein abundances, and the lowest variation among ribosomal protein abundances, which are expected to have equal amounts. Absolute quantitative proteome data enabled us to evaluate metabolic cost for protein synthesis and apparent catalytic activities of enzymes by integration with flux analysis. All the methods demonstrated similar ATP costs for protein synthesis for different cellular processes and that costs for expressing biomass synthesis related proteins were higher than those for energy generation. Importantly, catalytic activities of energy metabolism enzymes were an order or two higher than those of monomer synthesis. Interestingly, a staircase-like protein expression was demonstrated for most of the transcription units.

Acetate scavenging activity in Escherichia coli: interplay of acetyl-CoA synthetase and the PEP-glyoxylate cycle in chemostat cultures

Impairment of acetate production in Escherichia coli is crucial for the performance of many biotechnological processes. Aerobic production of acetate (or acetate overflow) results from changes in the expression of central metabolism genes. Acetyl-CoA synthetase scavenges extracellular acetate in glucose-limited cultures. Once converted to acetyl-CoA, it can be catabolized by the tricarboxylic acid cycle or the glyoxylate pathway. In this work, we assessed the significance of these pathways on acetate overflow during glucose excess and limitation. Gene expression, enzyme activities, and metabolic fluxes were studied in E. coli knock-out mutants related to the glyoxylate pathway operon and its regulators. The relevance of post-translational regulation by AceK-mediated phosphorylation of isocitrate dehydrogenase for pathway functionality was underlined. In chemostat cultures performed at increasing dilution rates, acetate overflow occurs when growing over a threshold glucose uptake rate. This threshold was not affected in a glyoxylate-pathway-deficient strain (lacking isocitrate lyase, the first enzyme of the pathway), indicating that it is not relevant for acetate overflow. In carbon-limited chemostat cultures, gluconeogenesis (maeB, sfcA, and pck), the glyoxylate operon and, especially, acetyl-CoA synthetase are upregulated. A mutant in acs (encoding acetyl-CoA synthetase) produced acetate at all dilution rates. This work demonstrates that, in E. coli, acetate production occurs at all dilution rates and that overflow is the result of unbalanced synthesis and scavenging activities. The over-expression of acetyl-CoA synthetase by cAMP-CRP-dependent induction limits this phenomenon in cultures consuming glucose at low rate, ensuring the recycling of the acetyl-CoA and acetyl-phosphate pools, although establishing an energy-dissipating substrate cycle.

Decrease of energy spilling in Escherichia coli continuous cultures with rising specific growth rate and carbon wasting

Background

Growth substrates, aerobic/anaerobic conditions, specific growth rate (μ) etc. strongly influence Escherichia coli cell physiology in terms of cell size, biomass composition, gene and protein expression. To understand the regulation behind these different phenotype properties, it is useful to know carbon flux patterns in the metabolic network which are generally calculated by metabolic flux analysis (MFA). However, rarely is biomass composition determined and carbon balance carefully measured in the same experiments which could possibly lead to distorted MFA results and questionable conclusions. Therefore, we carried out both detailed carbon balance and biomass composition analysis in the same experiments for more accurate quantitative analysis of metabolism and MFA.

Results

We applied advanced continuous cultivation methods (A-stat and D-stat) to continuously monitor E. coli K-12 MG1655 flux and energy metabolism dynamic responses to change of μ and glucose-acetate co-utilisation. Surprisingly, a 36% reduction of ATP spilling was detected with increasing μ and carbon wasting to non-CO₂ by-products under constant biomass yield. The apparent discrepancy between constant biomass yield and decline of ATP spilling could be explained by the rise of carbon wasting from 3 to 11% in the carbon balance which was revealed by the discovered novel excretion profile of E. coli pyrimidine pathway intermediates carbamoyl-phosphate, dihydroorotate and orotate. We found that carbon wasting patterns are dependent not only on μ, but also on glucose-acetate co-utilisation capability. Accumulation of these compounds was coupled to the two-phase acetate accumulation profile. Acetate overflow was observed in parallel with the reduction of TCA cycle and glycolysis fluxes, and induction of pentose phosphate pathway.

Conclusions

It can be concluded that acetate metabolism is one of the major regulating factors of central carbon metabolism. More importantly, our model calculations with actual biomass composition and detailed carbon balance analysis in steady state conditions with -omics data comparison demonstrate the importance of a comprehensive systems biology approach for more advanced understanding of metabolism and carbon re-routing mechanisms potentially leading to more successful metabolic engineering.