Escherichia coli can survive stress by noisy growth modulation

Drug tolerance and persistence in bacteria, fungi and cancer cells: Role of non-genetic heterogeneity

A common feature of bacterial, fungal and cancer cell populations upon treatment is the presence of tolerant and persistent cells able to survive, and sometimes grow, even in the presence of usually inhibitory or lethal drug concentrations, driven by non-genetic differences among individual cells in a population. Here we review and compare data obtained on drug survival in bacteria, fungi and cancer cells to unravel common characteristics and cellular pathways, and to point their singularities. This comparative work also allows to cross-fertilize ideas across fields. We particularly focus on the role of gene expression variability in the emergence of cell-cell non-genetic heterogeneity because it represents a possible common basic molecular process at the origin of most persistence phenomena and could be monitored and tuned to help improve therapeutic interventions.

Single cell technologies for monitoring protein secretion heterogeneity

Cell-to-cell heterogeneity presents challenges across various fields, from biomedicine to bioproduction, where precise cellular responses are vital. While single cell technologies have significantly enhanced our understanding of population heterogeneity, the predominant focus has been on monitoring intracellular compounds. Recognizing the added complexity introduced by the secretion system, in this review, we first provide a systematic overview of the distinct steps necessary for driving protein secretion. We discuss the various sources of noise acting from the synthesized preprotein to the secretory protein released based on a Gram-positive cellular system as a model. We next explore the applicability of single cell technologies for monitoring protein secretion throughout these functional stages. We also emphasize the importance of applying these single cell technologies for monitoring protein secretion during bioproduction.

Fluid flow drives phenotypic heterogeneity in bacterial growth and adhesion on surfaces

Bacteria often thrive in surface-attached communities, where they can form biofilms affording them multiple advantages. In this sessile form, fluid flow is a key component of their environments, renewing nutrients and transporting metabolic products and signaling molecules. It also controls colonization patterns and growth rates on surfaces, through bacteria transport, attachment and detachment. However, the current understanding of bacterial growth on surfaces neglects the possibility that bacteria may modulate their division behavior as a response to flow. Here, we employed single-cell imaging in microfluidic experiments to demonstrate that attached Escherichia coli cells can enter a growth arrest state while simultaneously enhancing their adhesion underflow. Despite utilizing clonal populations, we observed a non-uniform response characterized by bistable dynamics, with co-existing subpopulations of non-dividing and actively dividing bacteria. As the proportion of non-dividing bacteria increased with the applied flow rate, it resulted in a reduction in the average growth rate of bacterial populations on flow-exposed surfaces. Dividing bacteria exhibited asymmetric attachment, whereas non-dividing counterparts adhered to the surface via both cell poles. Hence, this phenotypic diversity allows bacterial colonies to combine enhanced attachment with sustained growth, although at a reduced rate, which may be a significant advantage in fluctuating flow conditions.

Stochastic Gene Expression in Proliferating Cells: Differing Noise Intensity in Single-Cell and Population Perspectives

Random fluctuations (noise) in gene expression can be studied from two complementary perspectives: following expression in a single cell over time or comparing expression between cells in a proliferating population at a given time. Here, we systematically investigated scenarios where both perspectives lead to different levels of noise in a given gene product. We first consider a stable protein, whose concentration is diluted by cellular growth, and the protein inhibits growth at high concentrations, establishing a positive feedback loop. For a stochastic model with molecular bursting of gene products, we analytically predict and contrast the steady-state distributions of protein concentration in both frameworks. Although positive feedback amplifies the noise in expression, this amplification is much higher in the population framework compared to following a single cell over time. We also study other processes that lead to different noise levels even in the absence of such dilution-based feedback. When considering randomness in the partitioning of molecules between daughters during mitosis, we find that in the single-cell perspective, the noise in protein concentration is independent of noise in the cell cycle duration. In contrast, partitioning noise is amplified in the population perspective by increasing randomness in cell-cycle time. Overall, our results show that the commonly used single-cell framework that does not account for proliferating cells can, in some cases, underestimate the noise in gene product levels. These results have important implications for studying the inter-cellular variation of different stress-related expression programs across cell types that are known to inhibit cellular growth.

Feedback between stochastic gene networks and population dynamics enables cellular decision-making

Phenotypic selection occurs when genetically identical cells are subject to different reproductive abilities due to cellular noise. Such noise arises from fluctuations in reactions synthesizing proteins and plays a crucial role in how cells make decisions and respond to stress or drugs. We propose a general stochastic agent-based model for growing populations capturing the feedback between gene expression and cell division dynamics. We devise a finite state projection approach to analyze gene expression and division distributions and infer selection from single-cell data in mother machines and lineage trees. We use the theory to quantify selection in multi-stable gene expression networks and elucidate that the trade-off between phenotypic switching and selection enables robust decision-making essential for synthetic circuits and developmental lineage decisions. Using live-cell data, we demonstrate that combining theory and inference provides quantitative insights into bet-hedging-like response to DNA damage and adaptation during antibiotic exposure in Escherichia coli.

Shaping of microbial phenotypes by trade-offs

Growth rate maximization is an important fitness strategy for microbes. However, the wide distribution of slow-growing oligotrophic microbes in ecosystems suggests that rapid growth is often not favored across ecological environments. In many circumstances, there exist trade-offs between growth and other important traits (e.g., adaptability and survival) due to physiological and proteome constraints. Investments on alternative traits could compromise growth rate and microbes need to adopt bet-hedging strategies to improve fitness in fluctuating environments. Here we review the mechanistic role of trade-offs in controlling bacterial growth and further highlight its ecological implications in driving the emergences of many important ecological phenomena such as co-existence, population heterogeneity and oligotrophic/copiotrophic lifestyles.

Synergistic and antagonistic effects of deterministic and stochastic cell-cell variations

By diversifying, cells in a clonal population can together overcome the limits of individuals. Diversity in single-cell growth rates allows the population to survive environmental stresses, such as antibiotics, and grow faster than the undiversified population. These functional cell-cell variations can arise stochastically, from noise in biochemical reactions, or deterministically, by asymmetrically distributing damaged components. While each of the mechanisms is well understood, the effect of the combined mechanisms is unclear. To evaluate the contribution of the deterministic component we developed a mathematical model by mapping the growing population to the Ising model. To analyze the combined effects of stochastic and deterministic contributions we introduced the analytical results of the Ising-mapping into an Euler-Lotka framework. Model results, confirmed by simulations and experimental data, show that deterministic cell-cell variations increase near-linearly with stress. As a consequence, we predict that the gain in population doubling time from cell-cell variations is primarily stochastic at low stress but may cross over to deterministic at higher stresses. Furthermore, we find that while the deterministic component minimizes population damage, stochastic variations antagonize this effect. Together our results may help identifying stress-tolerant pathogenic cells and thus inspire novel antibiotic strategies.

Bet-hedging and variability in plant development: seed germination and beyond

When future conditions are unpredictable, bet-hedging strategies can be advantageous. This can involve isogenic individuals producing different phenotypes, under the same environmental conditions. Ecological studies provide evidence that variability in seed germination time has been selected for as a bet-hedging strategy. We demonstrate how variability in germination time found in Arabidopsis could function as a bet-hedging strategy in the face of unpredictable lethal stresses. Despite a body of knowledge on how the degree of seed dormancy versus germination is controlled, relatively little is known about how differences between isogenic seeds in a batch are generated. We review proposed mechanisms for generating variability in germination time and the current limitations and new possibilities for testing the model predictions. We then look beyond germination to the role of variability in seedling and adult plant growth and review new technologies for quantification of noisy gene expression dynamics. We discuss evidence for phenotypic variability in plant traits beyond germination being under genetic control and propose that variability in stress response gene expression could function as a bet-hedging strategy. We discuss open questions about how noisy gene expression could lead to between-plant heterogeneity in gene expression and phenotypes. This article is part of a discussion meeting issue 'Causes and consequences of stochastic processes in development and disease'.

RpoS and the bacterial general stress response

SUMMARYThe general stress response (GSR) is a widespread strategy developed by bacteria to adapt and respond to their changing environments. The GSR is induced by one or multiple simultaneous stresses, as well as during entry into stationary phase and leads to a global response that protects cells against multiple stresses. The alternative sigma factor RpoS is the central GSR regulator in E. coli and conserved in most γ-proteobacteria. In E. coli, RpoS is induced under conditions of nutrient deprivation and other stresses, primarily via the activation of RpoS translation and inhibition of RpoS proteolysis. This review includes recent advances in our understanding of how stresses lead to RpoS induction and a summary of the recent studies attempting to define RpoS-dependent genes and pathways.

Deep model predictive control of gene expression in thousands of single cells

Gene expression is inherently dynamic, due to complex regulation and stochastic biochemical events. However, the effects of these dynamics on cell phenotypes can be difficult to determine. Researchers have historically been limited to passive observations of natural dynamics, which can preclude studies of elusive and noisy cellular events where large amounts of data are required to reveal statistically significant effects. Here, using recent advances in the fields of machine learning and control theory, we train a deep neural network to accurately predict the response of an optogenetic system in Escherichia coli cells. We then use the network in a deep model predictive control framework to impose arbitrary and cell-specific gene expression dynamics on thousands of single cells in real time, applying the framework to generate complex time-varying patterns. We also showcase the framework's ability to link expression patterns to dynamic functional outcomes by controlling expression of the tetA antibiotic resistance gene. This study highlights how deep learning-enabled feedback control can be used to tailor distributions of gene expression dynamics with high accuracy and throughput without expert knowledge of the biological system.

Dynamic Adaptation in Extant Porins Facilitates Antibiotic Tolerance in Energetic Escherichia coli

Bacteria can tolerate antibiotics despite lacking the genetic components for resistance. The prevailing notion is that tolerance results from depleted cellular energy or cell dormancy. In contrast to this view, many cells in the tolerant population of Escherichia coli can exhibit motility - a phenomenon that requires cellular energy, specifically, the proton-motive force (PMF). As these motile-tolerant cells are challenging to isolate from the heterogeneous tolerant population, their survival mechanism is unknown. Here, we discovered that motile bacteria segregate themselves from the tolerant population under micro-confinement, owing to their unique ability to penetrate micron-sized channels. Single-cell measurements on the motile-tolerant population showed that the cells retained a high PMF, but they did not survive through active efflux alone. By utilizing growth assays, single-cell fluorescence studies, and chemotaxis assays, we showed that the cells survived by dynamically inhibiting the function of existing porins in the outer membrane. A drug transport model for porin-mediated intake and efflux pump-mediated expulsion suggested that energetic tolerant cells withstand antibiotics by constricting their porins. The novel porin adaptation we have uncovered is independent of gene expression changes and may involve electrostatic modifications within individual porins to prevent extracellular ligand entry.

New layers of regulation of the general stress response sigma factor RpoS

The general stress response (GSR) sigma factor RpoS from Escherichia coli has emerged as one of the key paradigms for study of how numerous signal inputs are accepted at multiple levels into a single pathway for regulation of gene expression output. While many studies have elucidated the key pathways controlling the production and activity of this sigma factor, recent discoveries have uncovered still more regulatory mechanisms which feed into the network. Moreover, while the regulon of this sigma factor comprises a large proportion of the E. coli genome, the downstream expression levels of all the RpoS target genes are not identically affected by RpoS upregulation but respond heterogeneously, both within and between cells. This minireview highlights the most recent developments in our understanding of RpoS regulation and expression, in particular those which influence the regulatory network at different levels from previously well-studied pathways.

Quantifying Dynamic Phenotypic Heterogeneity in Resistant Escherichia coli under Translation-Inhibiting Antibiotics

Understanding the phenotypic heterogeneity of antibiotic-resistant bacteria following treatment and the transitions between different phenotypes is crucial for developing effective infection control strategies. The study expands upon previous work by explicating chloramphenicol-induced phenotypic heterogeneities in growth rate, gene expression, and morphology of resistant Escherichia coli using time-lapse microscopy. Correlating the bacterial growth rate and cspC expression, four interchangeable phenotypic subpopulations across varying antibiotic concentrations are identified, surpassing the previously described growth rate bistability. Notably, bacterial cells exhibiting either fast or slow growth rates can concurrently harbor subpopulations characterized by high and low gene expression levels, respectively. To elucidate the mechanisms behind this enhanced heterogeneity, a concise gene expression network model is proposed and the biological significance of the four phenotypes is further explored. Additionally, by employing Hidden Markov Model fitting and integrating the non-equilibrium landscape and flux theory, the real-time data encompassing diverse bacterial traits are analyzed. This approach reveals dynamic changes and switching kinetics in different cell fates, facilitating the quantification of observable behaviors and the non-equilibrium dynamics and thermodynamics at play. The results highlight the multi-dimensional heterogeneous behaviors of antibiotic-resistant bacteria under antibiotic stress, providing new insights into the compromised antibiotic efficacy, microbial response, and associated evolution processes.

Gene expression tradeoffs determine bacterial survival and adaptation to antibiotic stress

To optimize their fitness, cells face the crucial task of efficiently responding to various stresses. This necessitates striking a balance between conserving resources for survival and allocating resources for growth and division. The fundamental principles governing these tradeoffs is an outstanding challenge in the physics of living systems. In this study, we introduce a coarse-grained theoretical framework for bacterial physiology that establishes a connection between the physiological state of cells and their survival outcomes in dynamic environments, particularly in the context of antibiotic exposure. Predicting bacterial survival responses to varying antibiotic doses proves challenging due to the profound influence of the physiological state on critical parameters, such as the Minimum Inhibitory Concentration (MIC) and killing rates, even within an isogenic cell population. Our proposed theoretical model bridges the gap by linking extracellular antibiotic concentration and nutrient quality to intracellular damage accumulation and gene expression. This framework allows us to predict and explain the control of cellular growth rate, death rate, MIC and survival fraction in a wide range of time-varying environments. Surprisingly, our model reveals that cell death is rarely due to antibiotic levels being above the maximum physiological limit, but instead survival is limited by the inability to alter gene expression sufficiently quickly to transition to a less susceptible physiological state. Moreover, bacteria tend to overexpress stress response genes at the expense of reduced growth, conferring greater protection against further antibiotic exposure. This strategy is in contrast to those employed in different nutrient environments, in which bacteria allocate resources to maximize growth rate. This highlights an important tradeoff between the cellular capacity for growth and the ability to survive antibiotic exposure.

Chaos in a bacterial stress response

Cellular responses to environmental changes are often highly heterogeneous and exhibit seemingly random dynamics. The astonishing insight of chaos theory is that such unpredictable patterns can, in principle, arise without the need for any random processes, i.e., purely deterministically without noise. However, while chaos is well understood in mathematics and physics, its role in cell biology remains unclear because the complexity and noisiness of biological systems make testing difficult. Here, we show that chaos explains the heterogeneous response of Escherichia coli cells to oxidative stress. We developed a theoretical model of the gene expression dynamics and demonstrate that chaotic behavior arises from rapid molecular feedbacks that are coupled with cell growth dynamics and cell-cell interactions. Based on theoretical predictions, we then designed single-cell experiments to show we can shift gene expression from periodic oscillations to chaos on demand. Our work suggests that chaotic gene regulation can be employed by cell populations to generate strong and variable responses to changing environments.

Design and validation of a dual-fluorescence reporter system to monitor bacterial gene expression in the gut environment

Fluorescence-based reporter systems are valuable tools for studying gene expression dynamics in living cells. However, available strategies to follow gene expression in bacteria within their natural ecosystem that can be typically rich and complex are scarce. In this work, we designed a plasmid-based tool ensuring both the identification of a strain of interest in complex environments and the monitoring of gene expression through the combination of two distinct fluorescent proteins as reporter genes. The tool was validated in Escherichia coli to monitor the expression of eut genes involved in the catabolism of ethanolamine. We demonstrated that the constructed reporter strain gradually responds with a bimodal output to increasing ethanolamine concentrations during in vitro cultures. The reporter strain was next inoculated to mice, and flow cytometry was used to detect the reporter strain among the dense microbiota of intestinal samples and to analyze specifically the expression of eut genes. This novel dual-fluorescent reporter system would be helpful to evaluate transcriptional processes in bacteria within complex environments. KEY POINTS: • A reporter tool was developed to monitor bacterial gene expression in complex environments. • Ethanolamine utilization (eut) genes are expressed by commensal E. coli in the mouse gut. • Expression of eut genes follows a bimodal distribution.

Fitness cost associated with cell phenotypic switching drives population diversification dynamics and controllability

Isogenic cell populations can cope with stress conditions by switching to alternative phenotypes. Even if it can lead to increased fitness in a natural context, this feature is typically unwanted for a range of applications (e.g., bioproduction, synthetic biology, and biomedicine) where it tends to make cellular response unpredictable. However, little is known about the diversification profiles that can be adopted by a cell population. Here, we characterize the diversification dynamics for various systems (bacteria and yeast) and for different phenotypes (utilization of alternative carbon sources, general stress response and more complex development patterns). Our results suggest that the diversification dynamics and the fitness cost associated with cell switching are coupled. To quantify the contribution of the switching cost on population dynamics, we design a stochastic model that let us reproduce the dynamics observed experimentally and identify three diversification regimes, i.e., constrained (at low switching cost), dispersed (at medium and high switching cost), and bursty (for very high switching cost). Furthermore, we use a cell-machine interface called Segregostat to demonstrate that different levels of control can be applied to these diversification regimes, enabling applications involving more precise cellular responses.

System-wide optimization of an orthogonal translation system with enhanced biological tolerance

Over the past two decades, synthetic biological systems have revolutionized the study of cellular physiology. The ability to site-specifically incorporate biologically relevant non-standard amino acids using orthogonal translation systems (OTSs) has proven particularly useful, providing unparalleled access to cellular mechanisms modulated by post-translational modifications, such as protein phosphorylation. However, despite significant advances in OTS design and function, the systems-level biology of OTS development and utilization remains underexplored. In this study, we employ a phosphoserine OTS (pSerOTS) as a model to systematically investigate global interactions between OTS components and the cellular environment, aiming to improve OTS performance. Based on this analysis, we design OTS variants to enhance orthogonality by minimizing host process interactions and reducing stress response activation. Our findings advance understanding of system-wide OTS:host interactions, enabling informed design practices that circumvent deleterious interactions with host physiology while improving OTS performance and stability. Furthermore, our study emphasizes the importance of establishing a pipeline for systematically profiling OTS:host interactions to enhance orthogonality and mitigate mechanisms underlying OTS-mediated host toxicity.

The σB alternative sigma factor circuit modulates noise to generate different types of pulsing dynamics

Single-cell approaches are revealing a high degree of heterogeneity, or noise, in gene expression in isogenic bacteria. How gene circuits modulate this noise in gene expression to generate robust output dynamics is unclear. Here we use the Bacillus subtilis alternative sigma factor σB as a model system for understanding the role of noise in generating circuit output dynamics. σB controls the general stress response in B. subtilis and is activated by a range of energy and environmental stresses. Recent single-cell studies have revealed that the circuit can generate two distinct outputs, stochastic pulsing and a single pulse response, but the conditions under which each response is generated are under debate. We implement a stochastic mathematical model of the σB circuit to investigate this and find that the system's core circuit can generate both response types. This is despite one response (stochastic pulsing) being stochastic in nature, and the other (single response pulse) being deterministic. We demonstrate that the main determinant for whichever response is generated is the degree with which the input pathway activates the core circuit, although the noise properties of the input pathway also biases the system towards one or the other type of output. Thus, our work shows how stochastic modelling can reveal the mechanisms behind non-intuitive gene circuit output dynamics.

Regulation by cyclic di-GMP attenuates dynamics and enhances robustness of bimodal curli gene activation in Escherichia coli

Curli amyloid fibers are a major constituent of the extracellular biofilm matrix formed by bacteria of the Enterobacteriaceae family. Within Escherichia coli biofilms, curli gene expression is limited to a subpopulation of bacteria, leading to heterogeneity of extracellular matrix synthesis. Here we show that bimodal activation of curli gene expression also occurs in well-mixed planktonic cultures of E. coli, resulting in all-or-none stochastic differentiation into distinct subpopulations of curli-positive and curli-negative cells at the entry into the stationary phase of growth. Stochastic curli activation in individual E. coli cells could further be observed during continuous growth in a conditioned medium in a microfluidic device, which further revealed that the curli-positive state is only metastable. In agreement with previous reports, regulation of curli gene expression by the second messenger c-di-GMP via two pairs of diguanylate cyclase and phosphodiesterase enzymes, DgcE/PdeH and DgcM/PdeR, modulates the fraction of curli-positive cells. Unexpectedly, removal of this regulatory network does not abolish the bimodality of curli gene expression, although it affects dynamics of activation and increases heterogeneity of expression levels among individual cells. Moreover, the fraction of curli-positive cells within an E. coli population shows stronger dependence on growth conditions in the absence of regulation by DgcE/PdeH and DgcM/PdeR pairs. We thus conclude that, while not required for the emergence of bimodal curli gene expression in E. coli, this c-di-GMP regulatory network attenuates the frequency and dynamics of gene activation and increases its robustness to cellular heterogeneity and environmental variation.

Phenotypic heterogeneity in the bacterial oxidative stress response is driven by cell-cell interactions

Genetically identical bacterial cells commonly display different phenotypes. This phenotypic heterogeneity is well known for stress responses, where it is often explained as bet hedging against unpredictable environmental threats. Here, we explore phenotypic heterogeneity in a major stress response of Escherichia coli and find it has a fundamentally different basis. We characterize the response of cells exposed to hydrogen peroxide (H2O2) stress in a microfluidic device under constant growth conditions. A machine-learning model reveals that phenotypic heterogeneity arises from a precise and rapid feedback between each cell and its immediate environment. Moreover, we find that the heterogeneity rests upon cell-cell interaction, whereby cells shield each other from H2O2 via their individual stress responses. Our work shows how phenotypic heterogeneity in bacterial stress responses can emerge from short-range cell-cell interactions and result in a collective phenotype that protects a large proportion of the population.

Microfluidic single-cell scale-down systems: introduction, application, and future challenges

Performance losses during the scaling-up of bioprocesses from the laboratory to the production scale are common obstacles caused by the formation of concentration gradients in bioreactors. To overcome these obstacles, so-called scale-down bioreactors are used to analyze selected large-scale conditions and are one of the most important predictive tools for the successful transfer of bioprocesses from the lab to the industrial scale. In this regard, cellular behavior is usually measured as an averaged value, neglecting possible cell-to-cell heterogeneity within the culture. In contrast, microfluidic single-cell cultivation (MSCC) systems offer the possibility of understanding cellular processes on a single-cell level. To date, most MSCC systems have a limited choice of cultivation parameters that are not representative of bioprocess-relevant environmental conditions. Herein, we critically review recent advances in MSCC that allow the cultivation and analysis of cells under dynamic (bioprocess-relevant) environmental conditions. Finally, we discuss what technological advances and efforts are needed to bridge the gap between current MSCC systems and the use of these systems as single-cell scale-down devices.

Quiescence of Escherichia coli Aerosols to Survive Mechanical Stress during High-Velocity Collection

A low cutpoint wetted wall bioaerosol sampling cyclone (LCP-WWC), with an aerosol sampling flow rate of 300 L/min at 55″ H2O pressure drop and a continuous liquid outflow rate of about 0.2 mL/min, was developed by upgrading an existing system. The laboratory strain Escherichia coli MG1655 was aerosolized using a six-jet Collison Nebulizer and collected at high velocity using the LCP-WWC for 10 min with different collection liquids. Each sample was quantitated during a 15-day archiving period after aerosolization for culturable counts (CFUs) and gene copy numbers (GCNs) using microbial plating and whole-cell quantitative polymerase chain (qPCR) reaction. The samples were analyzed for protein composition and antimicrobial resistance using protein gel electrophoresis and disc diffusion susceptibility testing. Aerosolization and collection were followed by an initial period of quiescence or dormancy. After 2 days of archiving at 4 °C and RT, the bacteria exhibited increased culturability and antibiotic resistance (ABR), especially to cell wall inhibitors (ampicillin and cephalothin). The number of resistant bacteria on Day 2 increased nearly four-times compared to the number of cells at the initial time of collection. The mechanical stress of aerosolization and high-velocity sampling likely stunned the cells triggering a response of dormancy, though with continued synthesis of vital proteins for survival. This study shows that an increase in intensity in environmental conditions surrounding airborne bacteria affects their ability to grow and their potential to develop antimicrobial resistance.

Memory Effect on the Survival of Deinococcus radiodurans after Exposure in Near Space

Near space (20 to 100 km in altitude) is an extreme environment with high radiation and extreme cold, making it difficult for organisms to survive. However, many studies had shown that there were still microbes living in this extremely harsh environment. It was particularly important to study which factors affected the survival of microorganisms living in near space after exposure to irradiation, as this was related to many studies, such as studies of radioresistance mechanisms, panspermia hypothesis, long-distance microbial transfer, and developing extraterrestrial habitats. Survival after radiation was probably influenced by the growth condition before radiation, which is called the memory effect. In this research, we used different growth conditions to affect the growth of Deinococcus radiodurans and lyophilized bacteria in exponential phase to maintain the physiological state at this stage. Then high-altitude scientific balloon exposure experiments were carried out by using the Chinese Academy of Sciences Balloon-Borne Astrobiology Platform (CAS-BAP) at Dachaidan, Qinghai, China (37°44'N, 95°21'E). The aim was to investigate which factors influence survival after near-space exposure. The results suggested that there was a memory effect on the survival of D. radiodurans after exposure. If the differences in growth rate were caused by differences in nutrition, the survival rate and growth rate were positively correlated. Moreover, the addition of paraquat and Mn2+ during the growth phase can also increase survival. This finding may help to deepen the understanding of the mechanics of radiation protection and provide relevant evidence for many studies, such as of long-distance transfer of microorganisms in near space. IMPORTANCE Earth's near space is an extreme environment with high radiation and extreme cold. Which factors affect the survival of microbes in near space is related to many studies, such as studies of radioresistance mechanisms, panspermia hypothesis, long-distance microbial transfer, and developing extraterrestrial habitats. We performed several exposure experiments with Deinococcus radiodurans in near space to investigate which factors influence the survival rate after near-space exposure; that is, there was a relationship between survival after radiation and the growth condition before radiation. The results suggested that there was a memory effect on the survival of D. radiodurans after exposure. This finding may help to deepen the understanding of the mechanism of radiation protection and provide relevant evidence for many studies, such as of long-distance transfer of microorganisms in near space.

Navigating Environmental Transitions: the Role of Phenotypic Variation in Bacterial Responses

The ability of bacteria to respond to changes in their environment is critical to their survival, allowing them to withstand stress, form complex communities, and induce virulence responses during host infection. A remarkable feature of many of these bacterial responses is that they are often variable across individual cells, despite occurring in an isogenic population exposed to a homogeneous environmental change, a phenomenon known as phenotypic heterogeneity. Phenotypic heterogeneity can enable bet-hedging or division of labor strategies that allow bacteria to survive fluctuating conditions. Investigating the significance of phenotypic heterogeneity in environmental transitions requires dynamic, single-cell data. Technical advances in quantitative single-cell measurements, imaging, and microfluidics have led to a surge of publications on this topic. Here, we review recent discoveries on single-cell bacterial responses to environmental transitions of various origins and complexities, from simple diauxic shifts to community behaviors in biofilm formation to virulence regulation during infection. We describe how these studies firmly establish that this form of heterogeneity is prevalent and a conserved mechanism by which bacteria cope with fluctuating conditions. We end with an outline of current challenges and future directions for the field. While it remains challenging to predict how an individual bacterium will respond to a given environmental input, we anticipate that capturing the dynamics of the process will begin to resolve this and facilitate rational perturbation of environmental responses for therapeutic and bioengineering purposes.

Intercellular interaction mechanisms promote diversity in intracellular ATP concentration in Escherichia coli populations

In fluctuating environments, many microorganisms acquire phenotypic heterogeneity as a survival tactic to increase the likelihood of survival of the overall population. One example of this interindividual heterogeneity is the diversity of ATP concentration among members of Escherichia coli populations under glucose deprivation. Despite the importance of such environmentally driven phenotypic heterogeneity, how the differences in intracellular ATP concentration emerge among individual E. coli organisms is unknown. In this study, we focused on the mechanism through which individual E. coli achieve high intracellular ATP concentrations. First, we measured the ATP retained by E. coli over time when cultured at low (0.1 mM) and control (22.2 mM) concentrations of glucose and obtained the chronological change in ATP concentrations. Then, by comparing these chronological change of ATP concentrations and analyzing whether stochastic state transitions, periodic oscillations, cellular age, and intercellular communication-which have been reported as molecular biological mechanisms for generating interindividual heterogeneity-are involved, we showed that the appearance of high ATP-holding individuals observed among E. coli can be explained only by intercellular transmission. By performing metabolomic analysis of post-culture medium, we revealed a significant increase in the ATP, especially at low glucose, and that the number of E. coli that retain significantly higher ATP can be controlled by adding large amounts of ATP to the medium, even in populations cultured under control glucose concentrations. These results reveal for the first time that ATP-mediated intercellular transmission enables some individuals in E. coli populations grown at low glucose to retain large amounts of ATP.

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.

Microfluidics for long-term single-cell time-lapse microscopy: Advances and applications

Cells are inherently dynamic, whether they are responding to environmental conditions or simply at equilibrium, with biomolecules constantly being made and destroyed. Due to their small volumes, the chemical reactions inside cells are stochastic, such that genetically identical cells display heterogeneous behaviors and gene expression profiles. Studying these dynamic processes is challenging, but the development of microfluidic methods enabling the tracking of individual prokaryotic cells with microscopy over long time periods under controlled growth conditions has led to many discoveries. This review focuses on the recent developments of one such microfluidic device nicknamed the mother machine. We overview the original device design, experimental setup, and challenges associated with this platform. We then describe recent methods for analyzing experiments using automated image segmentation and tracking. We further discuss modifications to the experimental setup that allow for time-varying environmental control, replicating batch culture conditions, cell screening based on their dynamic behaviors, and to accommodate a variety of microbial species. Finally, this review highlights the discoveries enabled by this technology in diverse fields, such as cell-size control, genetic mutations, cellular aging, and synthetic biology.

An expanded whole-cell model of E. coli links cellular physiology with mechanisms of growth rate control

Growth and environmental responses are essential for living organisms to survive and adapt to constantly changing environments. In order to simulate new conditions and capture dynamic responses to environmental shifts in a developing whole-cell model of E. coli, we incorporated additional regulation, including dynamics of the global regulator guanosine tetraphosphate (ppGpp), along with dynamics of amino acid biosynthesis and translation. With the model, we show that under perturbed ppGpp conditions, small molecule feedback inhibition pathways, in addition to regulation of expression, play a role in ppGpp regulation of growth. We also found that simulations with dysregulated amino acid synthesis pathways provide average amino acid concentration predictions that are comparable to experimental results but on the single-cell level, concentrations unexpectedly show regular fluctuations. Additionally, during both an upshift and downshift in nutrient availability, the simulated cell responds similarly with a transient increase in the mRNA:rRNA ratio. This additional simulation functionality should support a variety of new applications and expansions of the E. coli Whole-Cell Modeling Project.

Alteration of DNA supercoiling serves as a trigger of short-term cold shock repressed genes of E. coli

Cold shock adaptability is a key survival skill of gut bacteria of warm-blooded animals. Escherichia coli cold shock responses are controlled by a complex multi-gene, timely-ordered transcriptional program. We investigated its underlying mechanisms. Having identified short-term, cold shock repressed genes, we show that their responsiveness is unrelated to their transcription factors or global regulators, while their single-cell protein numbers' variability increases after cold shock. We hypothesized that some cold shock repressed genes could be triggered by high propensity for transcription locking due to changes in DNA supercoiling (likely due to DNA relaxation caused by an overall reduction in negative supercoiling). Concomitantly, we found that nearly half of cold shock repressed genes are also highly responsive to gyrase inhibition (albeit most genes responsive to gyrase inhibition are not cold shock responsive). Further, their response strengths to cold shock and gyrase inhibition correlate. Meanwhile, under cold shock, nucleoid density increases, and gyrases and nucleoid become more colocalized. Moreover, the cellular energy decreases, which may hinder positive supercoils resolution. Overall, we conclude that sensitivity to diminished negative supercoiling is a core feature of E. coli's short-term, cold shock transcriptional program, and could be used to regulate the temperature sensitivity of synthetic circuits.

The transcription factor network of E. coli steers global responses to shifts in RNAP concentration

The robustness and sensitivity of gene networks to environmental changes is critical for cell survival. How gene networks produce specific, chronologically ordered responses to genome-wide perturbations, while robustly maintaining homeostasis, remains an open question. We analysed if short- and mid-term genome-wide responses to shifts in RNA polymerase (RNAP) concentration are influenced by the known topology and logic of the transcription factor network (TFN) of Escherichia coli. We found that, at the gene cohort level, the magnitude of the single-gene, mid-term transcriptional responses to changes in RNAP concentration can be explained by the absolute difference between the gene's numbers of activating and repressing input transcription factors (TFs). Interestingly, this difference is strongly positively correlated with the number of input TFs of the gene. Meanwhile, short-term responses showed only weak influence from the TFN. Our results suggest that the global topological traits of the TFN of E. coli shape which gene cohorts respond to genome-wide stresses.

Nutrient Gradients Mediate Complex Colony-Level Antibiotic Responses in Structured Microbial Populations

Antibiotic treatments often fail to eliminate bacterial populations due to heterogeneity in how individual cells respond to the drug. In structured bacterial populations such as biofilms, bacterial metabolism and environmental transport processes lead to an emergent phenotypic structure and self-generated nutrient gradients toward the interior of the colony, which can affect cell growth, gene expression and susceptibility to the drug. Even in single cells, survival depends on a dynamic interplay between the drug's action and the expression of resistance genes. How expression of resistance is coordinated across populations in the presence of such spatiotemporal environmental coupling remains elusive. Using a custom microfluidic device, we observe the response of spatially extended microcolonies of tetracycline-resistant E. coli to precisely defined dynamic drug regimens. We find an intricate interplay between drug-induced changes in cell growth and growth-dependent expression of resistance genes, resulting in the redistribution of metabolites and the reorganization of growth patterns. This dynamic environmental feedback affects the regulation of drug resistance differently across the colony, generating dynamic phenotypic structures that maintain colony growth during exposure to high drug concentrations and increase population-level resistance to subsequent exposures. A mathematical model linking metabolism and the regulation of gene expression is able to capture the main features of spatiotemporal colony dynamics. Uncovering the fundamental principles that govern collective mechanisms of antibiotic resistance in spatially extended populations will allow the design of optimal drug regimens to counteract them.

RNase E and HupB dynamics foster mycobacterial cell homeostasis and fitness

RNA turnover is a primary source of gene expression variation, in turn promoting cellular adaptation. Mycobacteria leverage reversible mRNA stabilization to endure hostile conditions. Although RNase E is essential for RNA turnover in several species, its role in mycobacterial single-cell physiology and functional phenotypic diversification remains unexplored. Here, by integrating live-single-cell and quantitative-mass-spectrometry approaches, we show that RNase E forms dynamic foci, which are associated with cellular homeostasis and fate, and we discover a versatile molecular interactome. We show a likely interaction between RNase E and the nucleoid-associated protein HupB, which is particularly pronounced during drug treatment and infection, where phenotypic diversity increases. Disruption of RNase E expression affects HupB levels, impairing Mycobacterium tuberculosis growth homeostasis during treatment, intracellular replication, and host spread. Our work lays the foundation for targeting the RNase E and its partner HupB, aiming to undermine M. tuberculosis cellular balance, diversification capacity, and persistence.

•

Single mycobacterial cells exhibit phenotypic variation in RNase E expression

•

RNase E is implicated in the maintenance of mycobacterial cell growth homeostasis

•

RNase E and HupB show a functional interplay in single mycobacterial cells

•

RNase E-HupB disruption impairs Mycobacterium tuberculosis fate under drug and in macrophages

Dynamic gene expression and growth underlie cell-to-cell heterogeneity in Escherichia coli stress response

Individual bacteria that share identical genomes and growth environments can display substantial cell-to-cell differences in expression of stress-response genes and single-cell growth rates. This phenotypic heterogeneity can impact the survival of single cells facing sudden stress. However, the windows of time that cells spend in vulnerable or tolerant states are often unknown. We quantify the temporal expression of a suite of stress-response reporters, while simultaneously monitoring growth. We observe pulsatile expression across genes with a range of stress-response functions, finding that single-cell growth rates are often anticorrelated with reporter levels. These dynamic phenotypic differences have a concrete link to function, in which individual cells undergoing a pulse of elevated expression and slow growth are predisposed to survive antibiotic exposure.

Cell-to-cell heterogeneity in gene expression and growth can have critical functional consequences, such as determining whether individual bacteria survive or die following stress. Although phenotypic variability is well documented, the dynamics that underlie it are often unknown. This information is important because dramatically different outcomes can arise from gradual versus rapid changes in expression and growth. Using single-cell time-lapse microscopy, we measured the temporal expression of a suite of stress-response reporters in Escherichia coli, while simultaneously monitoring growth rate. In conditions without stress, we found several examples of pulsatile expression. Single-cell growth rates were often anticorrelated with reporter levels, with changes in growth preceding changes in expression. These dynamics have functional consequences, which we demonstrate by measuring survival after challenging cells with the antibiotic ciprofloxacin. Our results suggest that fluctuations in both gene expression and growth dynamics in stress-response networks have direct consequences on survival.

Sequence-dependent model of genes with dual σ factor preference

Escherichia coli uses σ factors to quickly control large gene cohorts during stress conditions. While most of its genes respond to a single σ factor, approximately 5% of them have dual σ factor preference. The most common are those responsive to both σ⁷⁰, which controls housekeeping genes, and σ³⁸, which activates genes during stationary growth and stresses. Using RNA-seq and flow-cytometry measurements, we show that 'σ⁷⁰⁺³⁸ genes' are nearly as upregulated in stationary growth as 'σ³⁸ genes'. Moreover, we find a clear quantitative relationship between their promoter sequence and their response strength to changes in σ³⁸ levels. We then propose and validate a sequence dependent model of σ⁷⁰⁺³⁸ genes, with dual sensitivity to σ³⁸ and σ⁷⁰, that is applicable in the exponential and stationary growth phases, as well in the transient period in between. We further propose a general model, applicable to other stresses and σ factor combinations. Given this, promoters controlling σ⁷⁰⁺³⁸ genes (and variants) could become important building blocks of synthetic circuits with predictable, sequence-dependent sensitivity to transitions between the exponential and stationary growth phases.

Toward Multiplexed Optogenetic Circuits

Owing to its ubiquity and easy availability in nature, light has been widely employed to control complex cellular behaviors. Light-sensitive proteins are the foundation to such diverse and multilevel adaptive regulations in a large range of organisms. Due to their remarkable properties and potential applications in engineered systems, exploration and engineering of natural light-sensitive proteins have significantly contributed to expand optogenetic toolboxes with tailor-made performances in synthetic genetic circuits. Progressively, more complex systems have been designed in which multiple photoreceptors, each sensing its dedicated wavelength, are combined to simultaneously coordinate cellular responses in a single cell. In this review, we highlight recent works and challenges on multiplexed optogenetic circuits in natural and engineered systems for a dynamic regulation breakthrough in biotechnological applications.

Genome-wide gene expression noise in Escherichia coli is condition-dependent and determined by propagation of noise through the regulatory network

Although it is well appreciated that gene expression is inherently noisy and that transcriptional noise is encoded in a promoter’s sequence, little is known about the extent to which noise levels of individual promoters vary across growth conditions. Using flow cytometry, we here quantify transcriptional noise in Escherichia coli genome-wide across 8 growth conditions and find that noise levels systematically decrease with growth rate, with a condition-dependent lower bound on noise. Whereas constitutive promoters consistently exhibit low noise in all conditions, regulated promoters are both more noisy on average and more variable in noise across conditions. Moreover, individual promoters show highly distinct variation in noise across conditions. We show that a simple model of noise propagation from regulators to their targets can explain a significant fraction of the variation in relative noise levels and identifies TFs that most contribute to both condition-specific and condition-independent noise propagation. In addition, analysis of the genome-wide correlation structure of various gene properties shows that gene regulation, expression noise, and noise plasticity are all positively correlated genome-wide and vary independently of variations in absolute expression, codon bias, and evolutionary rate. Together, our results show that while absolute expression noise tends to decrease with growth rate, relative noise levels of genes are highly condition-dependent and determined by the propagation of noise through the gene regulatory network.

Genome-wide flow cytometry measurements reveal that gene expression noise in bacteria is highly condition-dependent; while absolute noise levels of all genes decrease with growth-rate, theoretical modeling shows that the relative noise levels of different genes are determined by the propagation of noise through the gene regulatory network (GRN). Thus GRN structure controls not only mean expression but also noise levels.

Feedback regulation and coordination of the main metabolism for bacterial growth and metabolic engineering for amino acid fermentation

Living organisms such as bacteria are often exposed to continuous changes in the nutrient availability in nature. Therefore, bacteria must constantly monitor the environmental condition, and adjust the metabolism quickly adapting to the change in the growth condition. For this, bacteria must orchestrate (coordinate and integrate) the complex and dynamically changing information on the environmental condition. In particular, the central carbon metabolism (CCM), monomer synthesis, and macromolecular synthesis must be coordinately regulated for the efficient growth. It is a grand challenge in bioscience, biotechnology, and synthetic biology to understand how living organisms coordinate the metabolic regulation systems. Here, we consider the integrated sensing of carbon sources by the phosphotransferase system (PTS), and the feed-forward/feedback regulation systems incorporated in the CCM in relation to the pool sizes of flux-sensing metabolites and αketoacids. We also consider the metabolic regulation of amino acid biosynthesis (as well as purine and pyrimidine biosyntheses) paying attention to the feedback control systems consisting of (fast) enzyme level regulation with (slow) transcriptional regulation. The metabolic engineering for the efficient amino acid production by bacteria such as Escherichia coli and Corynebacterium glutamicum is also discussed (in relation to the regulation mechanisms). The amino acid synthesis is important for determining the rate of ribosome biosynthesis. Thus, the growth rate control (growth law) is further discussed on the relationship between (p)ppGpp level and the ribosomal protein synthesis.

Evolutionary trade-offs between growth and survival: The delicate balance between reproductive success and longevity in bacteria

All living cells strive to allocate cellular resources in a way that promotes maximal evolutionary fitness. While there are many competing demands for resources the main decision making process centres on whether to proceed with growth and reproduction or to "hunker down" and invest in protection and survival (or to strike an optimal balance between these two processes). The transcriptional programme active at any given time largely determines which of these competing processes is dominant. At the top of the regulatory hierarchy are the sigma factors that commandeer the transcriptional machinery and determine which set of promoters are active at any given time. The regulatory inputs controlling their activity are therefore often highly complex, with multiple layers of regulation, allowing relevant environmental information to produce the most beneficial response. The tension between growth and survival is also evident in the developmental programme necessary to promote biofilm formation, which is typically associated with low growth rates and enhanced long-term survival. Nucleotide second messengers and energy pools (ATP/ADP levels) play critical roles in determining the fate of individual cells. Regulatory small RNAs frequently play important roles in the decision making processes too. In this review we discuss the trade-off that exists between reproduction and persistence in bacteria and discuss some of the recent advances in this fascinating field.

Tunable phenotypic variability through an autoregulatory alternative sigma factor circuit

Genetically identical individuals in bacterial populations can display significant phenotypic variability. This variability can be functional, for example by allowing a fraction of stress prepared cells to survive an otherwise lethal stress. The optimal fraction of stress prepared cells depends on environmental conditions. However, how bacterial populations modulate their level of phenotypic variability remains unclear. Here we show that the alternative sigma factor σV circuit in Bacillus subtilis generates functional phenotypic variability that can be tuned by stress level, environmental history and genetic perturbations. Using single-cell time-lapse microscopy and microfluidics, we find the fraction of cells that immediately activate σV under lysozyme stress depends on stress level and on a transcriptional memory of previous stress. Iteration between model and experiment reveals that this tunability can be explained by the autoregulatory feedback structure of the sigV operon. As predicted by the model, genetic perturbations to the operon also modulate the response variability. The conserved sigma-anti-sigma autoregulation motif is thus a simple mechanism for bacterial populations to modulate their heterogeneity based on their environment.

An ABA-GA bistable switch can account for natural variation in the variability of Arabidopsis seed germination time

Genetically identical plants growing in the same conditions can display heterogeneous phenotypes. Here we use Arabidopsis seed germination time as a model system to examine phenotypic variability and its underlying mechanisms. We show extensive variation in seed germination time variability between Arabidopsis accessions and use a multiparent recombinant inbred population to identify two genetic loci involved in this trait. Both loci include genes implicated in modulating abscisic acid (ABA) sensitivity. Mutually antagonistic regulation between ABA, which represses germination, and gibberellic acid (GA), which promotes germination, underlies the decision to germinate and can act as a bistable switch. A simple stochastic model of the ABA-GA network shows that modulating ABA sensitivity can generate the range of germination time distributions we observe experimentally. We validate the model by testing its predictions on the effects of exogenous hormone addition. Our work provides a foundation for understanding the mechanism and functional role of phenotypic variability in germination time.

Pulses and delays, anticipation and memory: seeing bacterial stress responses from a single-cell perspective

Stress responses are crucial for bacteria to survive harmful conditions that they encounter in the environment. Although gene regulatory mechanisms underlying stress responses in bacteria have been thoroughly characterised for decades, recent advances in imaging technologies helped to uncover previously hidden dynamics and heterogeneity that become visible at the single-cell level. Despite the diversity of stress response mechanisms, certain dynamic regulatory features are frequently seen in single cells, such as pulses, delays, stress anticipation and memory effects. Often, these dynamics are highly variable across cells. While any individual cell may not achieve an optimal stress response, phenotypic diversity can provide a benefit at the population level. In this review, we highlight microscopy studies that offer novel insights into how bacteria sense stress, regulate protective mechanisms, cope with response delays and prepare for future environmental challenges. These studies showcase developments in the single-cell imaging toolbox including gene expression reporters, FRET, super-resolution microscopy and single-molecule tracking, as well as microfluidic techniques to manipulate cells and create defined stress conditions.

Bacterial phenotypic heterogeneity in DNA repair and mutagenesis

Genetically identical cells frequently exhibit striking heterogeneity in various phenotypic traits such as their morphology, growth rate, or gene expression. Such non-genetic diversity can help clonal bacterial populations overcome transient environmental challenges without compromising genome stability, while genetic change is required for long-term heritable adaptation. At the heart of the balance between genome stability and plasticity are the DNA repair pathways that shield DNA from lesions and reverse errors arising from the imperfect DNA replication machinery. In principle, phenotypic heterogeneity in the expression and activity of DNA repair pathways can modulate mutation rates in single cells and thus be a source of heritable genetic diversity, effectively reversing the genotype-to-phenotype dogma. Long-standing evidence for mutation rate heterogeneity comes from genetics experiments on cell populations, which are now complemented by direct measurements on individual living cells. These measurements are increasingly performed using fluorescence microscopy with a temporal and spatial resolution that enables localising, tracking, and counting proteins with single-molecule sensitivity. In this review, we discuss which molecular processes lead to phenotypic heterogeneity in DNA repair and consider the potential consequences on genome stability and dynamics in bacteria. We further inspect these concepts in the context of DNA damage and mutation induced by antibiotics.

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.

Functional roles of microbial cell-to-cell heterogeneity and emerging technologies for analysis and control

Clonal cell populations often display significant cell-to-cell phenotypic heterogeneity, even when maintained under constant external conditions. This variability can result from the inherently stochastic nature of transcription and translation processes, which leads to varying numbers of transcripts and proteins per cell. Here, we showcase studies that reveal links between stochastic cellular events and biological functions in isogenic microbial populations. Then, we highlight emerging tools from engineering, computation, and synthetic and molecular biology that enable precise measurement, control, and analysis of gene expression noise in microorganisms. The capabilities offered by this sophisticated toolbox will shape future directions in the field and generate insight into the behavior of living systems at the single-cell level.

Function, Evolution, and Composition of the RpoS Regulon in Escherichia coli

For many bacteria, successful growth and survival depends on efficient adaptation to rapidly changing conditions. In Escherichia coli, the RpoS alternative sigma factor plays a central role in the adaptation to many suboptimal growth conditions by controlling the expression of many genes that protect the cell from stress and help the cell scavenge nutrients. Neither RpoS or the genes it controls are essential for growth and, as a result, the composition of the regulon and the nature of RpoS control in E. coli strains can be variable. RpoS controls many genetic systems, including those affecting pathogenesis, phenotypic traits including metabolic pathways and biofilm formation, and the expression of genes needed to survive nutrient deprivation. In this review, I review the origin of RpoS and assess recent transcriptomic and proteomic studies to identify features of the RpoS regulon in specific clades of E. coli to identify core functions of the regulon and to identify more specialized potential roles for the regulon in E. coli subgroups.

Quantifying plasmid dynamics using single-cell microfluidics and image bioinformatics

Multicopy plasmids play an important role in bacterial ecology and evolution by accelerating the rate of adaptation and providing a platform for rapid gene amplification and evolutionary rescue. Despite the relevance of plasmids in bacterial evolutionary dynamics, evaluating the population-level consequences of randomly segregating and replicating plasmids in individual cells remains a challenging problem, both in theory and experimentally. In recent years, technological advances in fluorescence microscopy and microfluidics have allowed studying temporal changes in gene expression by quantifying the fluorescent intensity of individual cells under controlled environmental conditions. In this paper, we will describe the manufacture, experimental setup, and data analysis pipeline of different microfluidic systems that can be used to study plasmid dynamics, both in single-cells and in populations. To illustrate the benefits and limitations of microfluidics to study multicopy plasmid dynamics, we will use an experimental model system consisting on Escherichia coli K12 carrying non-conjugative, multicopy plasmids (19 copies per cell, in average) encoding different fluorescent markers and β-lactam resistance genes. First, we will use an image-based flow cytometer to estimate changes in the allele distribution of a heterogeneous population under different selection regimes. Then we will use a mothermachine microfluidic device to obtain time-series of fluorescent intensity of individual cells to argue that plasmid segregation and replication dynamics are inherently stochastic processes. Finally, using a microchemostat, we track thousands of cells in time to reconstruct bacterial lineages and evaluate the allele frequency distributions that emerge in response to a range of selective pressures.