Metabolic variability in bioprocessing: implications of microbial phenotypic heterogeneity

Giving the cells what they need when they need it: Biosensor-based feeding control

"Giving the cells exactly what they need, when they need it" is the core idea behind the proposed bioprocess control strategy: operating bioprocess based on the physiological behavior of the microbial population rather than exclusive monitoring of environmental parameters. We are envisioning to achieve this through the use of genetically encoded biosensors combined with online flow cytometry (FCM) to obtain a time-dependent "physiological fingerprint" of the population. We developed a biosensor based on the glnA promoter (glnAp) and applied it for monitoring the nitrogen-related nutritional state of Escherichia coli. The functionality of the biosensor was demonstrated through multiple cultivation runs performed at various scales-from microplate to 20 L bioreactor. We also developed a fully automated bioreactor-FCM interface for on-line monitoring of the microbial population. Finally, we validated the proposed strategy by performing a fed-batch experiment where the biosensor signal is used as the actuator for a nitrogen feeding feedback control. This new generation of process control, -based on the specific needs of the cells, -opens the possibility of improving process development on a short timescale and therewith, the robustness and performance of fermentation processes.

Genomic and metabolic instability during long-term fermentation of an industrial Saccharomyces cerevisiae strain engineered for C5 sugar utilization

The genetic stability and metabolic robustness of production strains is one of the key criteria for the production of bio-based products by microbial fermentation on an industrial scale. These criteria were here explored in an industrial ethanol-producer strain of Saccharomyces cerevisiae able to co-ferment D-xylose and L-arabinose with glucose through the chromosomal integration of several copies of pivotal genes for the use of these pentose (C5) sugars. Using batch sequential cultures in a controlled bioreactor that mimics long-term fermentation in an industrial setting, this strain was found to exhibit significant fluctuations in D-xylose and L-arabinose consumption as early as the 50th generation and beyond. These fluctuations seem not related to the few low-consumption C5 sugar clones that appeared throughout the sequential batch cultures at a frequency lower than 1.5% and that were due to the reduction in the number of copies of transgenes coding for C5 sugar assimilation enzymes. Also, subpopulations enriched with low or high RAD52 expression, whose expression level was reported to be proportional to homologous recombination rate did not exhibit defect in C5-sugar assimilation, arguing that other mechanisms may be responsible for copy number variation of transgenes. Overall, this work highlighted the existence of genetic and metabolic instabilities in an industrial yeast which, although modest in our conditions, could be more deleterious in harsher industrial conditions, leading to reduced production performance.

Deep in situ microscopy for real-time analysis of mammalian cell populations in bioreactors

An in situ microscope based on pulsed transmitted light illumination via optical fiber was combined to artificial-intelligence to enable for the first time an online cell classification according to well-known cellular morphological features. A 848 192-image database generated during a lab-scale production process of antibodies was processed using a convolutional neural network approach chosen for its accurate real-time object detection capabilities. In order to induce different cell death routes, hybridomas were grown in normal or suboptimal conditions in a stirred tank reactor, in the presence of substrate limitation, medium addition, pH regulation problem or oxygen depletion. Using such an optical system made it possible to monitor real-time the evolution of different classes of animal cells, among which viable, necrotic and apoptotic cells. A class of viable cells displaying bulges in feast or famine conditions was also revealed. Considered as a breakthrough in the catalogue of process analytical tools, in situ microscopy powered by artificial-intelligence is also of great interest for research.

Monitoring Phenotype Heterogeneity at the Single-Cell Level within Bacillus Populations Producing Poly-3-hydroxybutyrate by Label-Free Super-resolution Infrared Imaging

Phenotypic heterogeneity is commonly found among bacterial cells within microbial populations due to intrinsic factors as well as equipping the organisms to respond to external perturbations. The emergence of phenotypic heterogeneity in bacterial populations, particularly in the context of using these bacteria as microbial cell factories, is a major concern for industrial bioprocessing applications. This is due to the potential impact on overall productivity by allowing the growth of subpopulations consisting of inefficient producer cells. Monitoring the spread of phenotypes across bacterial cells within the same population at the single-cell level is key to the development of robust, high-yield bioprocesses. Here, we discuss the novel development of optical photothermal infrared (O-PTIR) spectroscopy to probe phenotypic heterogeneity within Bacillus strains by monitoring the production of the bioplastic poly-3-hydroxybutyrate (PHB) at the single-cell level. Measurements obtained on single-point and in imaging mode show significant variability in the PHB content within bacterial cells, ranging from whether or not a cell produces PHB to variations in the intragranular biochemistry of PHB within bacterial cells. Our results show the ability of O-PTIR spectroscopy to probe PHB production at the single-cell level in a rapid, label-free, and semiquantitative manner. These findings highlight the potential of O-PTIR spectroscopy in single-cell microbial metabolomics as a whole-organism fingerprinting tool that can be used to monitor the dynamic of bacterial populations as well as for understanding their mechanisms for dealing with environmental stress, which is crucial for metabolic engineering research.

LTP2 hypomorphs show genotype-by-environment interaction in early seedling traits in Arabidopsis thaliana

Isogenic individuals can display seemingly stochastic phenotypic differences, limiting the accuracy of genotype-to-phenotype predictions. The extent of this phenotypic variation depends in part on genetic background, raising questions about the genes involved in controlling stochastic phenotypic variation. Focusing on early seedling traits in Arabidopsis thaliana, we found that hypomorphs of the cuticle-related gene LTP2 greatly increased variation in seedling phenotypes, including hypocotyl length, gravitropism and cuticle permeability. Many ltp2 hypocotyls were significantly shorter than wild-type hypocotyls while others resembled the wild type. Differences in epidermal properties and gene expression between ltp2 seedlings with long and short hypocotyls suggest a loss of cuticle integrity as the primary determinant of the observed phenotypic variation. We identified environmental conditions that reveal or mask the increased variation in ltp2 hypomorphs, and found that increased expression of its closest paralog LTP1 is necessary for ltp2 phenotypes. Our results illustrate how decreased expression of a single gene can generate starkly increased phenotypic variation in isogenic individuals in response to an environmental challenge.

Strengthening microbial cell factories for efficient production of bioactive molecules

High demand of bioactive molecules (food additives, antibiotics, plant growth enhancers, cosmetics, pigments and other commercial products) is the prime need for the betterment of human life where the applicability of the synthetic chemical product is on the saturation due to associated toxicity and ornamentations. It has been noticed that the discovery and productivity of such molecules in natural scenarios are limited due to low cellular yields as well as less optimized conventional methods. In this respect, microbial cell factories timely fulfilling the requirement of synthesizing bioactive molecules by improving production yield and screening more promising structural homologues of the native molecule. Where the robustness of the microbial host can be potentially achieved by taking advantage of cell engineering approaches such as tuning functional and adjustable factors, metabolic balancing, adapting cellular transcription machinery, applying high throughput OMICs tools, stability of genotype/phenotype, organelle optimizations, genome editing (CRISPER/Cas mediated system) and also by developing accurate model systems via machine-learning tools. In this article, we provide an overview from traditional to recent trends and the application of newly developed technologies, for strengthening the systemic approaches and providing future directions for enhancing the robustness of microbial cell factories to speed up the production of biomolecules for commercial purposes.

Application of an Escherichia coli triple reporter strain for at-line monitoring of single-cell physiology during L-phenylalanine production

Biotechnological production processes are sustainable approaches for the production of biobased components such as amino acids for food and feed industry. Scale-up from ideal lab-scale bioreactors to large-scale processes is often accompanied by loss in productivity. This may be related to population heterogeneities of cells originating from isogenic cultures that arise due to dynamic non-ideal conditions in the bioreactor. To better understand this phenomenon, deeper insights into single-cell physiologies in bioprocesses are mandatory before scale-up. Here, a triple reporter strain (3RP) was developed by chromosomally integrating the fluorescent proteins mEmerald, CyOFP1, and mTagBFP2 into the L-phenylalanine producing Escherichia coli strain FUS4 (pF81kan) to allow monitoring of growth, oxygen availability, and general stress response of the single cells. Functionality of the 3RP was confirmed in well-mixed lab-scale fed-batch processes with glycerol as carbon source in comparison to the strain without fluorescent proteins, leading to no difference in process performance. Fluorescence levels could successfully reflect the course of related process state variables, revealed population heterogeneities during the transition between different process phases and potentially subpopulations that exhibit superior process performance. Furthermore, indications were found for noise in gene expression as regulation strategy against environmental perturbation.

Controlling microbial co-culture based on substrate pulsing can lead to stability through differential fitness advantages

Microbial consortia are an exciting alternative for increasing the performances of bioprocesses for the production of complex metabolic products. However, the functional properties of microbial communities remain challenging to control, considering the complex interaction mechanisms occurring between co-cultured microbial species. Indeed, microbial communities are highly dynamic and can adapt to changing environmental conditions through complex mechanisms, such as phenotypic diversification. We focused on stabilizing a co-culture of Saccharomyces cerevisiae and Escherichia coli in continuous cultures. Our preliminary data pointed out that transient diauxic shifts could lead to stable co-culture by providing periodic fitness advantages to the yeast. Based on a computational toolbox called MONCKS (for MONod-type Co-culture Kinetic Simulation), we were able to predict the dynamics of diauxic shift for both species based on a cybernetic approach. This toolbox was further used to predict the frequency of diauxic shift to be applied to reach co-culture stability. These simulations were successfully reproduced experimentally in continuous bioreactors with glucose pulsing. Finally, based on a bet-hedging reporter, we observed that the yeast population exhibited an increased phenotypic diversification process in co-culture compared with mono-culture, suggesting that this mechanism could be the basis of the metabolic fitness of the yeast.

Being able to manipulate the dynamics of microbial co-cultures is a technical challenge that need to be addressed in order to get a deeper insight about how microbial communities are evolving in their ecological context, as well as for exploiting the potential offered by such communities in an applied context e.g., for setting up more robust bioprocesses relying on the use of several microbial species. In this study, we used continuous cultures of bacteria (E. coli) and yeast (S. cerevisiae) in order to demonstrate that a simple nutrient pulsing strategy can be used for adjusting the composition of the community with time. As expected, during growth on glucose, E. coli quickly outcompeted S. cerevisiae. However, when glucose is pulsed into the culture, increased metabolic fitness of the yeast was observed upon reconsumption of the main side metabolites i.e., acetate and ethanol, leading to a robust oscillating growth profile for both species. The optimal pulsing frequency was predicted based on a cybernetic version of a Monod growth model taking into account the main metabolic routes involved in the process. Considering the limited number of metabolic details needed, this cybernetic approach could be generalized to other communities.

Molecular Origins of Transcriptional Heterogeneity in Diazotrophic Klebsiella oxytoca

Phenotypic heterogeneity in clonal bacterial batch cultures has been shown for a range of bacterial systems; however, the molecular origins of such heterogeneity and its magnitude are not well understood. Under conditions of extreme low-nitrogen stress in the model diazotroph Klebsiella oxytoca, we found remarkably high heterogeneity of nifHDK gene expression, which codes for the structural genes of nitrogenase, one key enzyme of the global nitrogen cycle. This heterogeneity limited the bulk observed nitrogen-fixing capacity of the population. Using dual-probe, single-cell RNA fluorescent in situ hybridization, we correlated nifHDK expression with that of nifLA and glnK-amtB, which code for the main upstream regulatory components. Through stochastic transcription models and mutual information analysis, we revealed likely molecular origins for heterogeneity in nitrogenase expression. In the wild type and regulatory variants, we found that nifHDK transcription was inherently bursty, but we established that noise propagation through signaling was also significant. The regulatory gene glnK had the highest discernible effect on nifHDK variance, while noise from factors outside the regulatory pathway were negligible. Understanding the basis of inherent heterogeneity of nitrogenase expression and its origins can inform biotechnology strategies seeking to enhance biological nitrogen fixation. Finally, we speculate on potential benefits of diazotrophic heterogeneity in natural soil environments. IMPORTANCE Nitrogen is an essential micronutrient for both plant and animal life and naturally exists in both reactive and inert chemical forms. Modern agriculture is heavily reliant on nitrogen that has been "fixed" into a reactive form via the energetically expensive Haber-Bosch process, with significant environmental consequences. Nitrogen-fixing bacteria provide an alternative source of fixed nitrogen for use in both biotechnological and agricultural settings, but this relies on a firm understanding of how the fixation process is regulated within individual bacterial cells. We examined the cell-to-cell variability in the nitrogen-fixing behavior of Klebsiella oxytoca, a free-living bacterium. The significance of our research is in identifying not only the presence of marked variability but also the specific mechanisms that give rise to it. This understanding gives insight into both the evolutionary advantages of variable behavior as well as strategies for biotechnological applications.

Digital models in biotechnology: Towards multi-scale integration and implementation

Industrial biotechnology encompasses a large area of multi-scale and multi-disciplinary research activities. With the recent megatrend of digitalization sweeping across all industries, there is an increased focus in the biotechnology industry on developing, integrating and applying digital models to improve all aspects of industrial biotechnology. Given the rapid development of this field, we systematically classify the state-of-art modelling concepts applied at different scales in industrial biotechnology and critically discuss their current usage, advantages and limitations. Further, we critically analyzed current strategies to couple cell models with computational fluid dynamics to study the performance of industrial microorganisms in large-scale bioprocesses, which is of crucial importance for the bio-based production industries. One of the most challenging aspects in this context is gathering intracellular data under industrially relevant conditions. Towards comprehensive models, we discuss how different scale-down concepts combined with appropriate analytical tools can capture intracellular states of single cells. We finally illustrated how the efforts could be used to develop digitals models suitable for both cell factory design and process optimization at industrial scales in the future.

Development of a ribosome profiling protocol to study translation in Kluyveromyces marxianus

Kluyveromyces marxianus is an interesting and important yeast because of particular traits such as thermotolerance and rapid growth, and for applications in food and industrial biotechnology. For both understanding its biology and developing bioprocesses, it is important to understand how K. marxianus responds and adapts to changing environments. For this, a full suite of omics tools to measure and compare global patterns of gene expression and protein synthesis is needed. We report here the development of a ribosome profiling method for K. marxianus, which allows codon resolution of translation on a genome-wide scale by deep sequencing of ribosome locations on mRNAs. To aid in the analysis and sharing of ribosome profiling data, we added the K. marxianus genome as well as transcriptome and ribosome profiling data to the publicly accessible GWIPS-viz and Trips-Viz browsers. Users are able to upload custom ribosome profiling and RNA-Seq data to both browsers, therefore allowing easy analysis and sharing of data. We also provide a set of step-by-step protocols for the experimental and bioinformatic methods that we developed.

New insights into phenotypic heterogeneity for the distinct lipid accumulation of Schizochytrium sp. H016

BACKGROUND

Schizochytrium sp. is a marine heterotrophic protist and an important sustainable resource for high value-added docosahexaenoic acid in the future. The production of different phenotypes during the continuous subculture of Schizochytrium sp. results in a serious reduction in lipid yield and complicates the used of this strain in scientific research and industrial production. Hence, obtaining an improved understanding of the phenotypic differences and molecular mechanisms underlying the cell-to-cell heterogeneity of Schizochytrium sp. is necessary.

RESULTS

After continuous culture passage, Schizochytrium sp. H016 differentiated into two subpopulations with different morphologies and showed decreased capacity for lipid production. The presence of cell subpopulations with degraded lipid droplets led to a substantial decrease in overall lipid yield. Here, a rapid screening strategy based on fluorescence-activated cell sorting was proposed to classify and isolate subpopulations quickly in accordance with their lipid-producing capability. The final biomass and lipid yield of the subpopulation with high cell lipid content (i.e., H016-H) were 38.83 and 17.22 g/L, respectively, which were 2.07- and 5.38-fold higher than those of the subpopulation with low lipid content (i.e., H016-L), respectively. Subsequently, time‑resolved transcriptome analysis was performed to elucidate the mechanism of phenotypic heterogeneity in different subpopulations. Results showed that the expression of genes related to the cell cycle and lipid degradation was significantly upregulated in H016-L, whereas the metabolic pathways related to fatty acid synthesis and glyceride accumulation were remarkably upregulated in H016-H.

CONCLUSION

This study innovatively used flow cytometry combined with transcriptome technology to provide new insights into the phenotypic heterogeneity of different cell subpopulations of Schizochytrium sp. Furthermore, these results lay a strong foundation for guiding the breeding of oleaginous microorganisms with high lipid contents.

Xylanase production by Thermobacillus xylanilyticus is impaired by population diversification but can be mitigated based on the management of cheating behavior

Background

The microbial production of hemicellulasic cocktails is still a challenge for the biorefineries sector and agro-waste valorization. In this work, the production of hemicellulolytic enzymes by Thermobacillus xylanilyticus has been considered. This microorganism is of interest since it is able to produce an original set of thermostable hemicellulolytic enzymes, notably a xylanase GH11, Tx-xyn11. However, cell-to-cell heterogeneity impairs the production capability of the whole microbial population.

Results

Sequential cultivations of the strain on xylan as a carbon source has been considered in order to highlight and better understand this cell-to-cell heterogeneity. Successive cultivations pointed out a fast decrease of xylanase activity (loss of ~ 75%) and Tx-xyn11 gene expression after 23.5 generations. During serial cultivations on xylan, flow cytometry analyses pointed out that two subpopulations, differing at their light-scattering properties, were present. An increase of the recurrence of the subpopulation exhibiting low forward scatter (FSC) signal was correlated with a progressive loss of xylanase activity over several generations. Cell sorting and direct observation of the sorted subpopulations revealed that the low-FSC subpopulation was not sporulating, whereas the high-FSC subpopulation contained cells at the onset of the sporulation stage. The subpopulation differences (growth and xylanase activity) were assessed during independent growth. The low-FSC subpopulation exhibited a lag phase of 10 h of cultivation (and xylanase activities from 0.15 ± 0.21 to 3.89 ± 0.14 IU/mL along the cultivation) and the high-FSC subpopulation exhibited a lag phase of 5 h (and xylanase activities from 0.52 ± 0.00 to 4.43 ± 0.61 over subcultivations). Serial cultivations on glucose, followed by a switch to xylan led to a ~ 1.5-fold to ~ 15-fold improvement of xylanase activity, suggesting that alternating cultivation conditions could lead to an efficient population management strategy for the production of xylanase.

Conclusions

Taken altogether, the data from this study point out that a cheating behavior is responsible for the progressive reduction in xylanase activity during serial cultivations of T. xylanilyticus. Alternating cultivation conditions between glucose and xylan could be used as an efficient strategy for promoting population stability and higher enzymatic productivity from this bacterium.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12934-022-01762-z.

Advances in automated real-time flow cytometry for monitoring of bioreactor processes

Flow cytometry and its technological possibilities have greatly advanced in the past decade as analysis tool for single cell properties and population distributions of different cell types in bioreactors. Along the way, some solutions for automated real-time flow cytometry (ART-FCM) were developed for monitoring of bioreactor processes without operator interference over extended periods with variable sampling frequency. However, there is still great potential for ART-FCM to evolve and possibly become a standard application in bioprocess monitoring and process control. This review first addresses different components of an ART-FCM, including the sampling device, the sample-processing unit, the unit for sample delivery to the flow cytometer and the settings for measurement of pre-processed samples. Also, available algorithms are presented for automated data analysis of multi-parameter fluorescence datasets derived from ART-FCM experiments. Furthermore, challenges are discussed for integration of fluorescence-activated cell sorting into an ART-FCM setup for isolation and separation of interesting subpopulations that can be further characterized by for instance omics-methods. As the application of ART-FCM is especially of interest for bioreactor process monitoring, including investigation of population heterogeneity and automated process control, a summary of already existing setups for these purposes is given. Additionally, the general future potential of ART-FCM is addressed.

Inducible Population Quality Control of Engineered Bacillus subtilis for Improved N-Acetylneuraminic Acid Biosynthesis

Biosynthesis by microorganisms using renewable feedstocks is an important approach for realizing sustainable chemical manufacturing. However, cell-to-cell variation in biosynthesis capability during fermentation restricts the robustness and efficiency of bioproduction, hampering the industrialization of biosynthesis. Herein, we developed an inducible population quality control system (iPopQC) for dynamically modulating the producing and nonproducing subpopulations of engineered Bacillus subtilis, which was constructed via inducible promoter- and metabolite-responsive biosensor-based genetic circuit for regulating essential genes. Moreover, iPopQC achieved a 1.97-fold increase in N-acetylneuraminic acid (NeuAc) titer by enriching producing cell subpopulation during cultivation, representing 52% higher than that of previous PopQC. Strains with double-output iPopQC cocoupling the expression of double essential genes with NeuAc production improved production robustness further, retaining NeuAc production throughout 96 h of fermentation, upon which the strains cocoupling one essential gene expression with NeuAc production abolished the production ability.

Reducing phenotypic instabilities of a microbial population during continuous cultivation based on cell switching dynamics

Predicting the fate of individual cells among a microbial population (i.e., growth and gene expression) remains a challenge, especially when this population is exposed to very dynamic environmental conditions, such as those encountered during continuous cultivation. Indeed, the dynamic nature of a continuous cultivation process implies the potential diversification of the microbial population resulting in genotypic and phenotypic heterogeneity. The present work focused on the induction of the arabinose operon in Escherichia coli as a model system to study this diversification process in continuous cultivations. As a preliminary step, the green fluorescent protein (GFP) level triggered by an arabinose-inducible ParaBAD promoter was tracked by flow cytometry in chemostat cultivations with glucose-arabinose co-feeding. For a wide range of glucose-arabinose co-feeding concentrations in the chemostats, the simultaneous occurrence of GFP positive and negative subpopulation was observed. In the second set of experiments, continuous cultivation was performed by adding glucose continuously and arabinose based on the capability of individual cells to switch from low GFP to high GFP expression states, performed with a technology setup called segregostat. In the segregostat cultivation mode, on-line flow cytometry analysis was used for adjusting the arabinose/glucose transitions based on the phenotypic switching profiles of the microbial population. This strategy allowed finding an appropriate arabinose pulsing frequency, leading to prolonged maintenance of the induction level with a limited increase in the phenotypic diversity for more than 60 generations. The results suggest that the steady forcing of individual cells into a given phenotypic trajectory may not be the best strategy for controlling cell populations. Instead, allowing individual cells to switch periodically around a predefined threshold seems to be a more robust strategy leading to oscillations, but within a predictable cell population behavior range.

A yeast FRET biosensor enlightens cAMP signaling

The cAMP-PKA signaling cascade in budding yeast regulates adaptation to changing environments. We developed yEPAC, a FRET-based biosensor for cAMP measurements in yeast. We used this sensor with flow cytometry for high-throughput single cell-level quantification during dynamic changes in response to sudden nutrient transitions. We found that the characteristic cAMP peak differentiates between different carbon source transitions and is rather homogenous among single cells, especially for transitions to glucose. The peaks are mediated by a combination of extracellular sensing and intracellular metabolism. Moreover, the cAMP peak follows the Weber-Fechner law; its height scales with the relative, and not the absolute, change in glucose. Last, our results suggest that the cAMP peak height conveys information about prospective growth rates. In conclusion, our yEPAC-sensor makes possible new avenues for understanding yeast physiology, signaling, and metabolic adaptation.

A synthetic RNA-based biosensor for fructose-1,6-bisphosphate that reports glycolytic flux

RNA-based sensors for intracellular metabolites are a promising solution to the emerging issue of metabolic heterogeneity. However, their development, i.e., the conversion of an aptamer into an in vivo-functional intracellular metabolite sensor, still harbors challenges. Here, we accomplished this for the glycolytic flux-signaling metabolite, fructose-1,6-bisphosphate (FBP). Starting from in vitro selection of an aptamer, we constructed device libraries with a hammerhead ribozyme as actuator. Using high-throughput screening in yeast with fluorescence-activated cell sorting (FACS), next-generation sequencing, and genetic-environmental perturbations to modulate the intracellular FBP levels, we identified a sensor that generates ratiometric fluorescent readout. An abrogated response in sensor mutants and occurrence of two sensor conformations-revealed by RNA structural probing-indicated in vivo riboswitching activity. Microscopy showed that the sensor can differentiate cells with different glycolytic fluxes within yeast populations, opening research avenues into metabolic heterogeneity. We demonstrate the possibility to generate RNA-based sensors for intracellular metabolites for which no natural metabolite-binding RNA element exits.

Dynamic control in metabolic engineering: Theories, tools, and applications

Metabolic engineering has allowed the production of a diverse number of valuable chemicals using microbial organisms. Many biological challenges for improving bio-production exist which limit performance and slow the commercialization of metabolically engineered systems. Dynamic metabolic engineering is a rapidly developing field that seeks to address these challenges through the design of genetically encoded metabolic control systems which allow cells to autonomously adjust their flux in response to their external and internal metabolic state. This review first discusses theoretical works which provide mechanistic insights and design choices for dynamic control systems including two-stage, continuous, and population behavior control strategies. Next, we summarize molecular mechanisms for various sensors and actuators which enable dynamic metabolic control in microbial systems. Finally, important applications of dynamic control to the production of several metabolite products are highlighted, including fatty acids, aromatics, and terpene compounds. Altogether, this review provides a comprehensive overview of the progress, advances, and prospects in the design of dynamic control systems for improved titer, rate, and yield metrics in metabolic engineering.

Potential of Integrating Model-Based Design of Experiments Approaches and Process Analytical Technologies for Bioprocess Scale-Down

Typically, bioprocesses on an industrial scale are dynamic systems with a certain degree of variability, system inhomogeneities, and even population heterogeneities. Therefore, the scaling of such processes from laboratory to industrial scale and vice versa is not a trivial task. Traditional scale-down methodologies consider several technical parameters, so that systems on the laboratory scale tend to qualitatively reflect large-scale effects, but not the dynamic situation in an industrial bioreactor over the entire process, from the perspective of a cell. Supported by the enormous increase in computing power, the latest scientific focus is on the application of dynamic models, in combination with computational fluid dynamics to quantitatively describe cell behavior. These models allow the description of possible cellular lifelines which in turn can be used to derive a regime analysis for scale-down experiments. However, the approaches described so far, which were for a very few process examples, are very labor- and time-intensive and cannot be validated easily. In parallel, alternatives have been developed based on the description of the industrial process with hybrid process models, which describe a process mechanistically as far as possible in order to determine the essential process parameters with their respective variances. On-line analytical methods allow the characterization of population heterogeneity directly in the process. This detailed information from the industrial process can be used in laboratory screening systems to select relevant conditions in which the cell and process related parameters reflect the situation in the industrial scale. In our opinion, these technologies, which are available in research for modeling biological systems, in combination with process analytical techniques are so far developed that they can be implemented in industrial routines for faster development of new processes and optimization of existing ones.

A Monte Carlo method to estimate cell population heterogeneity from cell snapshot data

Variation is characteristic of all living systems. Laboratory techniques such as flow cytometry can probe individual cells, and, after decades of experimentation, it is clear that even members of genetically identical cell populations can exhibit differences. To understand whether variation is biologically meaningful, it is essential to discern its source. Mathematical models of biological systems are tools that can be used to investigate causes of cell-to-cell variation. From mathematical analysis and simulation of these models, biological hypotheses can be posed and investigated, then parameter inference can determine which of these is compatible with experimental data. Data from laboratory experiments often consist of "snapshots" representing distributions of cellular properties at different points in time, rather than individual cell trajectories. These data are not straightforward to fit using hierarchical Bayesian methods, which require the number of cell population clusters to be chosen a priori. Nor are they amenable to standard nonlinear mixed effect methods, since a single observation per cell is typically too few to estimate parameter variability. Here, we introduce a computational sampling method named "Contour Monte Carlo" (CMC) for estimating mathematical model parameters from snapshot distributions, which is straightforward to implement and does not require that cells be assigned to predefined categories. The CMC algorithm fits to snapshot probability distributions rather than raw data, which means its computational burden does not, like existing approaches, increase with the number of cells observed. Our method is appropriate for underdetermined systems, where there are fewer distinct types of observations than parameters to be determined, and where observed variation is mostly due to variability in cellular processes rather than experimental measurement error. This may be the case for many systems due to continued improvements in resolution of laboratory techniques. In this paper, we apply our method to quantify cellular variation for three biological systems of interest and provide Julia code enabling others to use this method.

Microbial metabolic noise

From the time a cell was first placed under the microscope, it became apparent that identifying two clonal cells that "look" identical is extremely challenging. Since then, cell-to-cell differences in shape, size, and protein content have been carefully examined, informing us of the ultimate limits that hinder two cells from occupying an identical phenotypic state. Here, we present recent experimental and computational evidence that similar limits emerge also in cellular metabolism. These limits pertain to stochastic metabolic dynamics and, thus, cell-to-cell metabolic variability, including the resulting adapting benefits. We review these phenomena with a focus on microbial metabolism and conclude with a brief outlook on the potential relationship between metabolic noise and adaptive evolution. This article is categorized under: Metabolic Diseases > Computational Models Metabolic Diseases > Biomedical Engineering.

Metabolic Engineering Strategies in Diatoms Reveal Unique Phenotypes and Genetic Configurations With Implications for Algal Genetics and Synthetic Biology

Diatoms are photosynthetic microeukaryotes that dominate phytoplankton populations and have increasing applicability in biotechnology. Uncovering their complex biology and elevating strains to commercial standards depends heavily on robust genetic engineering tools. However, engineering microalgal genomes predominantly relies on random integration of transgenes into nuclear DNA, often resulting in detrimental “position-effects” such as transgene silencing, integration into transcriptionally-inactive regions, and endogenous sequence disruption. With the recent development of extrachromosomal transgene expression via independent episomes, it is timely to investigate both strategies at the phenotypic and genomic level. Here, we engineered the model diatom Phaeodactylum tricornutum to produce the high-value heterologous monoterpenoid geraniol, which, besides applications as fragrance and insect repellent, is a key intermediate of high-value pharmaceuticals. Using high-throughput phenotyping we confirmed the suitability of episomes for synthetic biology applications and identified superior geraniol-yielding strains following random integration. We used third generation long-read sequencing technology to generate a complete analysis of all transgene integration events including their genomic locations and arrangements associated with high-performing strains at a genome-wide scale with subchromosomal detail, never before reported in any microalga. This revealed very large, highly concatenated insertion islands, offering profound implications on diatom functional genetics and next generation genome editing technologies, and is key for developing more precise genome engineering approaches in diatoms, including possible genomic safe harbour locations to support high transgene expression for targeted integration approaches. Furthermore, we have demonstrated that exogenous DNA is not integrated inadvertently into the nuclear genome of extrachromosomal-expression clones, an important characterisation of this novel engineering approach that paves the road to synthetic biology applications.

Statistical mechanics of interacting metabolic networks

We cast the metabolism of interacting cells within a statistical mechanics framework considering both the actual phenotypic capacities of each cell and its interaction with its neighbors. Reaction fluxes will be the components of high-dimensional spin vectors, whose values will be constrained by the stochiometry and the energy requirements of the metabolism. Within this picture, finding the phenotypic states of the population turns out to be equivalent to searching for the equilibrium states of a disordered spin model. We provide a general solution of this problem for arbitrary metabolic networks and interactions. We apply this solution to a simplified model of metabolism and to a complex metabolic network, the central core of Escherichia coli, and demonstrate that the combination of selective pressure and interactions defines a complex phenotypic space. We also present numerical results for cells fixed in a grid. These results reproduce the qualitative picture discussed for the mean-field model. Cells may specialize in producing or consuming metabolites complementing each other, and this is described by an equilibrium phase space with multiple minima, like in a spin-glass model.

An Improved ATP FRET Sensor For Yeast Shows Heterogeneity During Nutrient Transitions

Adenosine 5-triphosphate (ATP) is the main free energy carrier in metabolism. In budding yeast, shifts to glucose-rich conditions cause dynamic changes in ATP levels, but it is unclear how heterogeneous these dynamics are at a single-cell level. Furthermore, pH also changes and affects readout of fluorescence-based biosensors for single-cell measurements. To measure ATP changes reliably in single yeast cells, we developed yAT1.03, an adapted version of the AT1.03 ATP biosensor, that is pH-insensitive. We show that pregrowth conditions largely affect ATP dynamics during transitions. Moreover, single-cell analyses showed a large variety in ATP responses, which implies large differences of glycolytic startup between individual cells. We found three clusters of dynamic responses, and we show that a small subpopulation of wild-type cells reached an imbalanced state during glycolytic startup, characterized by low ATP levels. These results confirm the need for new tools to study dynamic responses of individual cells in dynamic environments.

Quantification of the dynamics of population heterogeneities in CHO cultures with stably integrated fluorescent markers

Cell population heterogeneities and their changes in mammalian cell culture processes are still not well characterized. In this study, the formation and dynamics of cell population heterogeneities were investigated with flow cytometry and stably integrated fluorescent markers based on the lentiviral gene ontology (LeGO) vector system. To achieve this, antibody-producing CHO cells were transduced with different LeGO vectors to stably express single or multiple fluorescent proteins. This enables the tracking of the transduced populations and is discussed in two case studies from the field of bioprocess engineering: In case study I, cells were co-transduced to express red, green, and blue fluorescent proteins and the development of sub-populations and expression heterogeneities were investigated in high passage cultivations (total 130 days). The formation of a fast-growing and more productive population was observed with a simultaneous increase in cell density and product titer. In case study II, different preculture growth phases and their influence on the population dynamics were investigated in mixed batch cultures with flow cytometry (offline and automated). Four cell line derivatives, each expressing a different fluorescent protein, were generated and cultivated for different time intervals, corresponding to different growth phases. Mixed cultures were inoculated from them, and changes in the composition of the cell populations were observed during the first 48 h of cultivation with reduced process productivity. In summary, we showed how the dynamics of population heterogeneities can be characterized. This represents a novel approach to investigate the dynamics of cell population heterogeneities under near-physiological conditions with changing productivity in mammalian cell culture processes.

Development and characterization of Escherichia coli triple reporter strains for investigation of population heterogeneity in bioprocesses

Background

Today there is an increasing demand for high yielding robust and cost efficient biotechnological production processes. Although cells in these processes originate from isogenic cultures, heterogeneity induced by intrinsic and extrinsic influences is omnipresent. To increase understanding of this mechanistically poorly understood phenomenon, advanced tools that provide insights into single cell physiology are needed.

Results

Two Escherichia coli triple reporter strains have been designed based on the industrially relevant production host E. coli BL21(DE3) and a modified version thereof, E. coli T7E2. The strains carry three different fluorescence proteins chromosomally integrated. Single cell growth is followed with EmeraldGFP (EmGFP)-expression together with the ribosomal promoter rrnB. General stress response of single cells is monitored by expression of sigma factor rpoS with mStrawberry, whereas expression of the nar-operon together with TagRFP657 gives information about oxygen limitation of single cells. First, the strains were characterized in batch operated stirred-tank bioreactors in comparison to wildtype E. coli BL21(DE3). Afterwards, applicability of the triple reporter strains for investigation of population heterogeneity in bioprocesses was demonstrated in continuous processes in stirred-tank bioreactors at different growth rates and in response to glucose and oxygen perturbation simulating gradients on industrial scale. Population and single cell level physiology was monitored evaluating general physiology and flow cytometry analysis of fluorescence distributions of the triple reporter strains. Although both triple reporter strains reflected physiological changes that were expected based on the expression characteristics of the marker proteins, the triple reporter strain based on E. coli T7E2 showed higher sensitivity in response to environmental changes. For both strains, noise in gene expression was observed during transition from phases of non-growth to growth. Apparently, under some process conditions, e.g. the stationary phase in batch cultures, the fluorescence response of EmGFP and mStrawberry is preserved, whereas TagRFP657 showed a distinct response.

Conclusions

Single cell growth, general stress response and oxygen limitation of single cells could be followed using the two triple reporter strains developed in this study. They represent valuable tools to study population heterogeneity in bioprocesses significantly increasing the level of information compared to the use of single reporter strains.

Microbioreactors for Process Development and Cell-Based Screening Studies

Microbioreactors (MBRs) have emerged as potent cultivation devices enabling automated small-scale experiments in parallel while enhancing their cost efficiency. The widespread use of MBRs has contributed to recent advances in industrial and pharmaceutical biotechnology, and they have proved to be indispensable tools in the development of many modern bioprocesses. Being predominantly applied in early stage process development, they open up new fields of research and enhance the efficacy of biotechnological product development. Their reduced reaction volume is associated with numerous inherent advantages - particularly the possibility for enabling parallel screening operations that facilitate high-throughput cultivations with reduced sample consumption (or the use of rare and expensive educts). As a result, multiple variables can be examined in a shorter time and with a lower expense. This leads to a simultaneous acceleration of research and process development along with decreased costs.MBRs range from simple miniaturized cultivations vessels (i.e., in the milliliter scale with limited possibilities for process control) to highly complex and automated small-scale microreactors with integrated sensors that allow for comprehensive screenings in very short time or a precise reflection of large-scale cultivation conditions. Progressive developments and improvements in manufacturing and automation techniques are already helping researchers to make use of the advantages that MBRs offer. This overview of current MBR systems surveys the diverse application for microbial and mammalian cell cultivations that have been developed in recent years.

Quantitation and Comparison of Phenotypic Heterogeneity Among Single Cells of Monoclonal Microbial Populations

Phenotypic heterogeneity within microbial populations arises even when the cells are exposed to putatively constant and homogeneous conditions. The outcome of this phenomenon can affect the whole function of the population, resulting in, for example, new "adapted" metabolic strategies and impacting its fitness at given environmental conditions. Accounting for phenotypic heterogeneity becomes thus necessary, due to its relevance in medical and applied microbiology as well as in environmental processes. Still, a comprehensive evaluation of this phenomenon requires a common and unique method of quantitation, which allows for the comparison between different studies carried out with different approaches. Consequently, in this study, two widely applicable indices for quantitation of heterogeneity were developed. The heterogeneity coefficient (HC) is valid when the population follows unimodal activity, while the differentiation tendency index (DTI) accounts for heterogeneity implying outbreak of subpopulations and multimodal activity. We demonstrated the applicability of HC and DTI for heterogeneity quantitation on stable isotope probing with nanoscale secondary ion mass spectrometry (SIP-nanoSIMS), flow cytometry, and optical microscopy datasets. The HC was found to provide a more accurate and precise measure of heterogeneity, being at the same time consistent with the coefficient of variation (CV) applied so far. The DTI is able to describe the differentiation in single-cell activity within monoclonal populations resolving subpopulations with low cell abundance, individual cells with similar phenotypic features (e.g., isotopic content close to natural abundance, as detected with nanoSIMS). The developed quantitation approach allows for a better understanding on the impact and the implications of phenotypic heterogeneity in environmental, medical and applied microbiology, microbial ecology, cell biology, and biotechnology.

Coupled metabolic-hydrodynamic modeling enabling rational scale-up of industrial bioprocesses

Metabolomics aims to address what and how regulatory mechanisms are coordinated to achieve flux optimality, different metabolic objectives as well as appropriate adaptations to dynamic nutrient availability. Recent decades have witnessed that the integration of metabolomics and fluxomics within the goal of synthetic biology has arrived at generating the desired bioproducts with improved bioconversion efficiency. Absolute metabolite quantification by isotope dilution mass spectrometry represents a functional readout of cellular biochemistry and contributes to the establishment of metabolic (structured) models required in systems metabolic engineering. In industrial practices, population heterogeneity arising from fluctuating nutrient availability frequently leads to performance losses, that is reduced commercial metrics (titer, rate, and yield). Hence, the development of more stable producers and more predictable bioprocesses can benefit from a quantitative understanding of spatial and temporal cell-to-cell heterogeneity within industrial bioprocesses. Quantitative metabolomics analysis and metabolic modeling applied in computational fluid dynamics (CFD)-assisted scale-down simulators that mimic industrial heterogeneity such as fluctuations in nutrients, dissolved gases, and other stresses can procure informative clues for coping with issues during bioprocessing scale-up. In previous studies, only limited insights into the hydrodynamic conditions inside the industrial-scale bioreactor have been obtained, which makes case-by-case scale-up far from straightforward. Tracking the flow paths of cells circulating in large-scale bioreactors is a highly valuable tool for evaluating cellular performance in production tanks. The "lifelines" or "trajectories" of cells in industrial-scale bioreactors can be captured using Euler-Lagrange CFD simulation. This novel methodology can be further coupled with metabolic (structured) models to provide not only a statistical analysis of cell lifelines triggered by the environmental fluctuations but also a global assessment of the metabolic response to heterogeneity inside an industrial bioreactor. For the future, the industrial design should be dependent on the computational framework, and this integration work will allow bioprocess scale-up to the industrial scale with an end in mind.

The effect of acetate on population heterogeneity in different cellular characteristics of Escherichia coli in aerobic batch cultures

Acetate as the major by-product in industrial-scale bioprocesses with Escherichia coli is found to decrease process efficiency as well as to be toxic to cells, which has several effects like a significant induction of cellular stress responses. However, the underlying phenomena are poorly explored. Therefore, we studied time-resolved population heterogeneity of the E. coli growth reporter strain MG1655/pGS20PrrnBGFPAAV expressing destabilized green fluorescent protein during batch growth on acetate and glucose as sole carbon sources. Additionally, we applied five fluorescent stains targeting different cellular properties (viability as well as metabolic and respiratory activity). Quantitative analysis of flow cytometry data verified that bacterial populations in the bioreactor are more heterogeneous in growth as well as stronger metabolically challenged during growth on acetate as sole carbon source, compared to growth on glucose or acetate after diauxic shift. Interestingly, with acetate as sole carbon source, significant subpopulations were found with some cells that seem to be more robust than the rest of the population. In conclusion, following batch cultures population heterogeneity was evident in all measured parameters. Our approach enabled a deeper study of heterogeneity during growth on the favored substrate glucose as well as on the toxic by-product acetate. Using a combination of activity fluorescent dyes proved to be an accurate and fast alternative as well as a supplement to the use of a reporter strain. However, the choice of combination of stains should be well considered depending on which population traits to aim for.

Growth-dependent recombinant product formation kinetics can be reproduced through engineering of glucose transport and is prone to phenotypic heterogeneity

Background

Escherichia coli W3110 and a group of six isogenic derivatives, each displaying distinct specific rates of glucose consumption were characterized to determine levels of GFP production and population heterogeneity. These strains have single or combinatory deletions in genes encoding phosphoenolpyruvate:sugar phosphotransferase system (PTS) permeases as PtsG and ManX, as well as common components EI, Hpr protein and EIIA, also the non-PTS Mgl galactose/glucose ABC transporter. They have been transformed for expressing GFP based on a lac-based expression vector, which is subject to bistability.

Results

These strains displayed specific glucose consumption and growth rates ranging from 1.75 to 0.45 g/g h and 0.54 to 0.16 h⁻¹, respectively. The rate of acetate production was strongly reduced in all mutant strains when compared with W3110/pV21. In bioreactor cultures, wild type W3110/pV21 produced 50.51 mg/L GFP, whereas strains WG/pV21 with inactive PTS IICBGlc and WGM/pV21 with the additional inactivation of PTS IIABMan showed the highest titers of GFP, corresponding to 342 and 438 mg/L, respectively. Moreover, we showed experimentally that bistable expression systems, as lac-based ones, induce strong phenotypic segregation among microbial populations.

Conclusions

We have demonstrated that reduction on glucose consumption rate in E. coli leads to an improvement of GFP production. Furthermore, from the perspective of phenotypic heterogeneity, we observed in this case that heterogeneous systems are also the ones leading to the highest performance. This observation suggests reconsidering the generally accepted proposition stating that phenotypic heterogeneity is generally unwanted in bioprocess applications.

Electronic supplementary material

The online version of this article (10.1186/s12934-019-1073-5) contains supplementary material, which is available to authorized users.

Maximum entropy and population heterogeneity in continuous cell cultures

Continuous cultures of mammalian cells are complex systems displaying hallmark phenomena of nonlinear dynamics, such as multi-stability, hysteresis, as well as sharp transitions between different metabolic states. In this context mathematical models may suggest control strategies to steer the system towards desired states. Although even clonal populations are known to exhibit cell-to-cell variability, most of the currently studied models assume that the population is homogeneous. To overcome this limitation, we use the maximum entropy principle to model the phenotypic distribution of cells in a chemostat as a function of the dilution rate. We consider the coupling between cell metabolism and extracellular variables describing the state of the bioreactor and take into account the impact of toxic byproduct accumulation on cell viability. We present a formal solution for the stationary state of the chemostat and show how to apply it in two examples. First, a simplified model of cell metabolism where the exact solution is tractable, and then a genome-scale metabolic network of the Chinese hamster ovary (CHO) cell line. Along the way we discuss several consequences of heterogeneity, such as: qualitative changes in the dynamical landscape of the system, increasing concentrations of byproducts that vanish in the homogeneous case, and larger population sizes.

Bioreactor-Scale Strategies for the Production of Recombinant Protein in the Yeast Yarrowia lipolytica

Recombinant protein production represents a multibillion-dollar market. Therefore, it constitutes an important research field both in academia and industry. The use of yeast as a cell factory presents several advantages such as ease of genetic manipulation, growth at high cell density, and the possibility of post-translational modifications. Yarrowia lipolytica is considered as one of the most attractive hosts due to its ability to metabolize raw substrate, to express genes at a high level, and to secrete protein in large amounts. In recent years, several reviews have been dedicated to genetic tools developed for this purpose. Though the construction of efficient cell factories for recombinant protein synthesis is important, the development of an efficient process for recombinant protein production in a bioreactor constitutes an equally vital aspect. Indeed, a sports car cannot drive fast on a gravel road. The aim of this review is to provide a comprehensive snapshot of process tools to consider for recombinant protein production in bioreactor using Y. lipolytica as a cell factory, in order to facilitate the decision-making for future strain and process engineering.

Estimation methods for heterogeneous cell population models in systems biology

Heterogeneity among individual cells is a characteristic and relevant feature of living systems. A range of experimental techniques to investigate this heterogeneity is available, and multiple modelling frameworks have been developed to describe and simulate the dynamics of heterogeneous populations. Measurement data are used to adjust computational models, which results in parameter and state estimation problems. Methods to solve these estimation problems need to take the specific properties of data and models into account. The aim of this review is to give an overview on the state of the art in estimation methods for heterogeneous cell population data and models. The focus is on models based on the population balance equation, but stochastic and individual-based models are also discussed. It starts with a brief discussion of common experimental approaches and types of measurement data that can be obtained in this context. The second part describes computational modelling frameworks for heterogeneous populations and the types of estimation problems occurring for these models. The third part starts with a discussion of observability and identifiability properties, after which the computational methods to solve the various estimation problems are described.

Pulsatile inputs achieve tunable attenuation of gene expression variability and graded multi-gene regulation

Many natural transcription factors are regulated in a pulsatile fashion, but it remains unknown whether synthetic gene expression systems can benefit from such dynamic regulation. Here we find, using a fast-acting, optogenetic transcription factor in Saccharomyces cerevisiae, that dynamic pulsatile signals reduce cell-to-cell variability in gene expression. We then show that by encoding such signals into a single input, expression mean and variability can be independently tuned. Further, we construct a light-responsive promoter library and demonstrate how pulsatile signaling also enables graded multi-gene regulation at fixed expression ratios, despite differences in promoter dose-response characteristics. Pulsatile regulation can thus lead to beneficial functional behaviors in synthetic biological systems, which previously required laborious optimization of genetic parts or the construction of synthetic gene networks.

Dynamic metabolic control: towards precision engineering of metabolism

Advances in metabolic engineering have led to the synthesis of a wide variety of valuable chemicals in microorganisms. The key to commercializing these processes is the improvement of titer, productivity, yield, and robustness. Traditional approaches to enhancing production use the “push–pull-block” strategy that modulates enzyme expression under static control. However, strains are often optimized for specific laboratory set-up and are sensitive to environmental fluctuations. Exposure to sub-optimal growth conditions during large-scale fermentation often reduces their production capacity. Moreover, static control of engineered pathways may imbalance cofactors or cause the accumulation of toxic intermediates, which imposes burden on the host and results in decreased production. To overcome these problems, the last decade has witnessed the emergence of a new technology that uses synthetic regulation to control heterologous pathways dynamically, in ways akin to regulatory networks found in nature. Here, we review natural metabolic control strategies and recent developments in how they inspire the engineering of dynamically regulated pathways. We further discuss the challenges of designing and engineering dynamic control and highlight how model-based design can provide a powerful formalism to engineer dynamic control circuits, which together with the tools of synthetic biology, can work to enhance microbial production.

Links between evolutionary processes and phenotypic robustness in microbes

The costs and benefits of random phenotypic heterogeneity in microbes have been vigorously debated and experimental tested for decades; yet, this conversation is largely independent from discussion of phenotypic robustness in other disciplines. In this review I connect microbial examples of stochasticity with studies on the ecological and population-genetic consequences of phenotypic variability. These topics illustrate the complexity of selection pressures on phenotypic robustness and provide inspiration that this complexity can be parsed with theoretical advances and the experimental power of microbial systems.

Bimodality of gene expression from yeast promoter can be instigated by DNA context, inducing conditions and strain background

Bimodality in gene expression is thought to provide a high phenotypic heterogeneity that can be favourable for adaptation or unfavourable notably in industrial processes that require stable and homogeneous properties. Whether this property is produced or suppressed in different conditions has been understudied. Here we identified tens of Saccharomyces cerevisiae genomic fragments conferring bimodal yEGFP expression on centromeric plasmid and studied some of these promoters in different DNA contexts, inducing conditions or strain backgrounds. First, we observed that the bimodal behaviour identified on plasmid is generally suppressed at the genomic level. Second, an inducible promoter such as the copper-regulated CUP1 promoter can produce bimodal expression in a time- and dose-dependent fashion. For a given copper sulphate concentration, a constant proportion of the subpopulation is induced and only the induction level of this subpopulation changed with induction duration, while for a same induction time, higher copper sulphate concentrations induced more cells at higher levels. Third, we showed that bimodality conferred by the CUP1 promoter in expression profile is strain background dependent, revealing epistasis in the generation of bimodality. The influence of these parameters on bimodality has to be taken into account when considering transgene expression for industrial microbial productions.

Efficient simulation of intrinsic, extrinsic and external noise in biochemical systems

Motivation

Biological cells operate in a noisy regime influenced by intrinsic, extrinsic and external noise, which leads to large differences of individual cell states. Stochastic effects must be taken into account to characterize biochemical kinetics accurately. Since the exact solution of the chemical master equation, which governs the underlying stochastic process, cannot be derived for most biochemical systems, approximate methods are used to obtain a solution.

Results

In this study, a method to efficiently simulate the various sources of noise simultaneously is proposed and benchmarked on several examples. The method relies on the combination of the sigma point approach to describe extrinsic and external variability and the τ-leaping algorithm to account for the stochasticity due to probabilistic reactions. The comparison of our method to extensive Monte Carlo calculations demonstrates an immense computational advantage while losing an acceptable amount of accuracy. Additionally, the application to parameter optimization problems in stochastic biochemical reaction networks is shown, which is rarely applied due to its huge computational burden. To give further insight, a MATLAB script is provided including the proposed method applied to a simple toy example of gene expression.

Availability and implementation

MATLAB code is available at Bioinformatics online.

Supplementary information

Supplementary data are available at Bioinformatics online.

Accelerating bacterial growth detection and antimicrobial susceptibility assessment in integrated picoliter droplet platform

There remains an urgent need for rapid diagnostic methods that can evaluate antibiotic resistance for pathogenic bacteria in order to deliver targeted antibiotic treatments. Toward this end, we present a rapid and integrated single-cell biosensing platform, termed dropFAST, for bacterial growth detection and antimicrobial susceptibility assessment. DropFAST utilizes a rapid resazurin-based fluorescent growth assay coupled with stochastic confinement of bacteria in 20 pL droplets to detect signal from growing bacteria after 1h incubation, equivalent to 2-3 bacterial replications. Full integration of droplet generation, incubation, and detection into a single, uninterrupted stream also renders this platform uniquely suitable for in-line bacterial phenotypic growth assessment. To illustrate the concept of rapid digital antimicrobial susceptibility assessment, we employ the dropFAST platform to evaluate the antibacterial effect of gentamicin on E. coli growth.

Analysis of population structures of the microalga Acutodesmus obliquus during lipid production using multi-dimensional single-cell analysis

Microalgae bear a great potential to produce lipids for biodiesel, feed, or even food applications. To understand the still not well-known single-cell dynamics during lipid production in microalgae, a novel single-cell analytical technology was applied to study a well-established model experiment. Multidimensional single-cell dynamics were investigated with a non-supervised image analysis technique that utilizes data from epi-fluorescence microscopy. Reliability of this technique was successfully proven via reference analysis. The technique developed was used to determine cell size, chlorophyll amount, neutral lipid amount, and deriving properties on a single-cellular level in cultures of the biotechnologically promising alga Acutodesmus obliquus. The results illustrated a high correlation between cell size and chlorophyll amount, but a very low and dynamic correlation between cell size, lipid amount, and lipid density. During growth conditions under nitrogen starvation, cells with low chlorophyll content tend to start the lipid production first and the cell suspension differentiated in two subpopulations with significantly different lipid contents. Such quantitative characterization of single-cell dynamics of lipid synthesizing algae was done for the first time and the potential of such simple technology is highly relevant to other biotechnological applications and to deeper investigate the process of microalgal lipid accumulation.

Population heterogeneity in microbial bioprocesses: origin, analysis, mechanisms, and future perspectives

Population heterogeneity is omnipresent in all bioprocesses even in homogenous environments. Its origin, however, is only so well understood that potential strategies like bet-hedging, noise in gene expression and division of labour that lead to population heterogeneity can be derived from experimental studies simulating the dynamics in industrial scale bioprocesses. This review aims at summarizing the current state of the different parts of single cell studies in bioprocesses. This includes setups to visualize different phenotypes of single cells, computational approaches connecting single cell physiology with environmental influence and special cultivation setups like scale-down reactors that have been proven to be useful to simulate large-scale conditions. A step in between investigation of populations and single cells is studying subpopulations with distinct properties that differ from the rest of the population with sub-omics methods which are also presented here. Moreover, the current knowledge about population heterogeneity in bioprocesses is summarized for relevant industrial production hosts and mixed cultures, as they provide the unique opportunity to distribute metabolic burden and optimize production processes in a way that is impossible in traditional monocultures. In the end, approaches to explain the underlying mechanism of population heterogeneity and the evidences found to support each hypothesis are presented. For instance, population heterogeneity serving as a bet-hedging strategy that is used as coordinated action against bioprocess-related stresses while at the same time spreading the risk between individual cells as it ensures the survival of least a part of the population in any environment the cells encounter.