Growth rate-dependent global effects on gene expression in bacteria

Achieving robust synthetic tolerance in industrial E. coli through negative auto-regulation of a DsrA-Hfq module

In industrial fermentation processes, microorganisms often encounter acid stress, which significantly impact their productivity. This study focused on the acid-resistant module composed of small RNA (sRNA) DsrA and the sRNA chaperone Hfq. Our previous study had shown that this module improved the cell growth of Escherichia coli MG1655 at low pH, but failed to obtain this desired phenotype in industrial strains. Here, we performed a quantitative analysis of DsrA-Hfq module to determine the optimal expression mode. We then assessed the potential of the CymR-based negative auto-regulation (NAR) circuit for industrial application, under different media, strains and pH levels. Growth assay at pH 4.5 revealed that NAR-05D04H circuit was the best acid-resistant circuit to improve the cell growth of E. coli MG1655. This circuit was robust and worked well in the industrial lysine-producing strain E. coli SCEcL3 at a starting pH of 6.8 and without pH control, resulting in a 250 % increase in lysine titer and comparable biomass in shaking flask fermentation compared to the parent strain. This study showed the practical application of NAR circuit in regulating DsrA-Hfq module, effectively and robustly improving the acid tolerance of industrial strains, which provides a new approach for breeding industrial strains with tolerance phenotype.

Transcriptional regulation of living materials via extracellular electron transfer

Engineered living materials combine the advantages of biological and synthetic systems by leveraging genetic and metabolic programming to control material-wide properties. Here, we demonstrate that extracellular electron transfer (EET), a microbial respiration process, can serve as a tunable bridge between live cell metabolism and synthetic material properties. In this system, EET flux from Shewanella oneidensis to a copper catalyst controls hydrogel cross-linking via two distinct chemistries to form living synthetic polymer networks. We first demonstrate that synthetic biology-inspired design rules derived from fluorescence parameterization can be applied toward EET-based regulation of polymer network mechanics. We then program transcriptional Boolean logic gates to govern EET gene expression, which enables design of computational polymer networks that mechanically respond to combinations of molecular inputs. Finally, we control fibroblast morphology using EET as a bridge for programmed material properties. Our results demonstrate how rational genetic circuit design can emulate physiological behavior in engineered living materials.

How total mRNA influences cell growth

Experimental observations tracing back to the 1960s imply that ribosome quantities play a prominent role in determining a cell's growth. Nevertheless, in biologically relevant scenarios, growth can also be influenced by the levels of mRNA and RNA polymerase. Here, we construct a quantitative model of biosynthesis providing testable scenarios for these situations. The model explores a theoretically motivated regime where RNA polymerases compete for genes and ribosomes for transcripts and gives general expressions relating growth rate, mRNA concentrations, ribosome, and RNA polymerase levels. On general grounds, the model predicts how the fraction of ribosomes in the proteome depends on total mRNA concentration and inspects an underexplored regime in which the trade-off between transcript levels and ribosome abundances sets the cellular growth rate. In particular, we show that the model predicts and clarifies three important experimental observations, in budding yeast and Escherichia coli bacteria: i) that the growth-rate cost of unneeded protein expression can be affected by mRNA levels, ii) that resource optimization leads to decreasing trends in mRNA levels at slow growth, and iii) that ribosome allocation may increase, stay constant, or decrease, in response to transcription-inhibiting antibiotics. Since the data indicate that a regime of joint limitation may apply in physiological conditions and not only to perturbations, we speculate that this regime is likely self-imposed.

A library of reporters of the global regulators of gene expression in Escherichia coli

The topology of the transcription factor network (TFN) of Escherichia coli is far from uniform, with 22 global regulator (GR) proteins controlling one-third of all genes. So far, their production rates cannot be tracked by comparable fluorescent proteins. We developed a library of fluorescent reporters for 16 GRs for this purpose. Each consists of a single-copy plasmid coding for green fluorescent protein (GFP) fused to the full-length copy of the native promoter. We tracked their activity in exponential and stationary growth, as well as under weak and strong stresses. We show that the reporters have high sensitivity and specificity to all stresses tested and detect single-cell variability in transcription rates. Given the influence of GRs on the TFN, we expect that the new library will contribute to dissecting global transcriptional stress-response programs of E. coli. Moreover, the library can be invaluable in bioindustrial applications that tune those programs to, instead of cell growth, favor productivity while reducing energy consumption.IMPORTANCECells contain thousands of genes. Many genes are involved in the control of cellular activities. Some activities require a few hundred genes to run largely synchronous transcriptional programs. To achieve this, cells have evolved global regulator (GR) proteins that can influence hundreds of genes simultaneously. We have engineered a library of Escherichia coli strains to track the levels over time of these, phenotypically critical, GRs. Each strain has a single-copy plasmid coding for a fast-maturing green fluorescent protein whose transcription is controlled by a copy of the natural GR promoter. By allowing the tracking of GR levels, with sensitivity and specificity, this library should become of wide use in scientific research on bacterial gene expression (from molecular to synthetic biology) and, later, be used in applications in therapeutics and bioindustries.

Homeostasis of cytoplasmic crowding by cell wall fluidization and ribosomal counterions

In bacteria, algae, fungi, and plant cells, the wall must expand in concert with cytoplasmic biomass production, otherwise cells would experience toxic molecular crowding1,2 or lyse. But how cells achieve expansion of this complex biomaterial in coordination with biosynthesis of macromolecules in the cytoplasm remains unexplained3, although recent works have revealed that these processes are indeed coupled4,5. Here, we report a striking increase of turgor pressure with growth rate in E. coli, suggesting that the speed of cell wall expansion is controlled via turgor. Remarkably, despite this increase in turgor pressure, cellular biomass density remains constant across a wide range of growth rates. By contrast, perturbations of turgor pressure that deviate from this scaling directly alter biomass density. A mathematical model based on cell wall fluidization by cell wall endopeptidases not only explains these apparently confounding observations but makes surprising quantitative predictions that we validated experimentally. The picture that emerges is that turgor pressure is directly controlled via counterions of ribosomal RNA. Elegantly, the coupling between rRNA and turgor pressure simultaneously coordinates cell wall expansion across a wide range of growth rates and exerts homeostatic feedback control on biomass density. This mechanism may regulate cell wall biosynthesis from microbes to plants and has important implications for the mechanism of action of antibiotics6.

The effect of DNA-binding proteins on insertion sequence element transposition upstream of the bgl operon in Escherichia coli

The bglGFB operon in Escherichia coli K-12 strain BW25113, encoding the proteins necessary for the uptake and metabolism of β-glucosides, is normally not expressed. Insertion of either IS1 or IS5 upstream of the bgl promoter activates expression of the operon only when the cell is starving in the presence of a β-glucoside, drastically increasing transcription and allowing the cell to survive and grow using this carbon source. Details surrounding the exact mechanism and regulation of the IS insertional event remain unclear. In this work, the role of several DNA-binding proteins in how they affect the rate of insertion upstream of bgl are examined via mutation assays and protocols measuring transcription. Both Crp and IHF exert a positive effect on insertional Bgl+ mutations when present, active, and functional in the cell. Our results characterize IHF's effect in conjunction with other mutations, show that IHF's effect on IS insertion into bgl also affects other operons, and indicate that it may exert its effect by binding to and altering the DNA conformation of IS1 and IS5 in their native locations, rather than by directly influencing transposase gene expression. In contrast, the cAMP-CRP complex acts directly upon the bgl operon by binding upstream of the promoter, presumably altering local DNA into a conformation that enhances IS insertion.

Full-Combinatorial Concentration Gradient Array with 3D Micro/Nanofluidics for Antibiotic Susceptibility Testing

Recent advancements in micro/nanofluidics have facilitated on-chip microscopy of cellular responses in a high-throughput and controlled microenvironment with the desired physicochemical properties. Despite its potential benefits to combination drug discovery, generating a complete combinatorial set of concentration gradients for multiple reagents in an array format remains challenging. The main reason is limited layouts of conventional micro/nanofluidic systems based on two-dimensional channel networks. In this paper, we present a device with three-dimensional (3D) interconnection of micro/nanochannels capable of generating a complete combinatorial set of concentration gradients for two reagents. The device was readily fabricated by laminating a pair of multilayered monolithic films containing a Christmas tree-like mixer, a cell culture chamber array, and through-holes, all within each single film. We assessed the reliable generation of a full-combinatorial concentration gradient array and validated it by using numerical analysis. We applied the proposed device to test the antibiotic susceptibility of bacterial cells in a convenient one-step manner. Furthermore, we explored the potential of the device to accommodate the arrayed complete combinatorial set for two or more drugs, while extending the capabilities of our laminated object manufacturing method for realizing 3D micro/nanofluidic systems.

Comprehensive Characterization of fucAO Operon Activation in Escherichia coli

Wildtype Escherichia coli cells cannot grow on L-1,2-propanediol, as the fucAO operon within the fucose (fuc) regulon is thought to be silent in the absence of L-fucose. Little information is available concerning the transcriptional regulation of this operon. Here, we first confirm that fucAO operon expression is highly inducible by fucose and is primarily attributable to the upstream operon promoter, while the fucO promoter within the 3'-end of fucA is weak and uninducible. Using 5'RACE, we identify the actual transcriptional start site (TSS) of the main fucAO operon promoter, refuting the originally proposed TSS. Several lines of evidence are provided showing that the fucAO locus is within a transcriptionally repressed region on the chromosome. Operon activation is dependent on FucR and Crp but not SrsR. Two Crp-cAMP binding sites previously found in the regulatory region are validated, where the upstream site plays a more critical role than the downstream site in operon activation. Furthermore, two FucR binding sites are identified, where the downstream site near the first Crp site is more important than the upstream site. Operon transcription relies on Crp-cAMP to a greater degree than on FucR. Our data strongly suggest that FucR mainly functions to facilitate the binding of Crp to its upstream site, which in turn activates the fucAO promoter by efficiently recruiting RNA polymerase.

Dynamical model of antibiotic responses linking expression of resistance genes to metabolism explains emergence of heterogeneity during drug exposures

Antibiotic responses in bacteria are highly dynamic and heterogeneous, with sudden exposure of bacterial colonies to high drug doses resulting in the coexistence of recovered and arrested cells. The dynamics of the response is determined by regulatory circuits controlling the expression of resistance genes, which are in turn modulated by the drug's action on cell growth and metabolism. Despite advances in understanding gene regulation at the molecular level, we still lack a framework to describe how feedback mechanisms resulting from the interdependence between expression of resistance and cell metabolism can amplify naturally occurring noise and create heterogeneity at the population level. To understand how this interplay affects cell survival upon exposure, we constructed a mathematical model of the dynamics of antibiotic responses that links metabolism and regulation of gene expression, based on the tetracycline resistancetetoperon inE. coli. We use this model to interpret measurements of growth and expression of resistance in microfluidic experiments, both in single cells and in biofilms. We also implemented a stochastic model of the drug response, to show that exposure to high drug levels results in large variations of recovery times and heterogeneity at the population level. We show that stochasticity is important to determine how nutrient quality affects cell survival during exposure to high drug concentrations. A quantitative description of how microbes respond to antibiotics in dynamical environments is crucial to understand population-level behaviors such as biofilms and pathogenesis.

High-throughput screening of human mercury exposure based on a low-cost naked eye-recognized biosensing platform

Whole-cell biosensors could be helpful for in situ disease diagnosis. However, their use in analyzing biological samples has been hindered by unstable responses, low signal enhancement, and growth inhibition in complex media. Here, we offered a solution by building a visual whole-cell biosensor for urinary mercury determination. With deoxyviolacein as the preferred signal for the mercury biosensor for the first time, it enabled the quantitative detection of urinary mercury with a favorable linear range from 1.57 to 100 nM. The biosensor can accurately diagnose urine mercury levels exceeding the biological exposure index with 95.8% accuracy. Thus, our study provided a biosensing platform with great potential to serve as a stable, user-friendly, and high-throughput alternative for the daily monitoring or estimating of urinary mercury.

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.

A hybrid transistor with transcriptionally controlled computation and plasticity

Organic electrochemical transistors (OECTs) are ideal devices for translating biological signals into electrical readouts and have applications in bioelectronics, biosensing, and neuromorphic computing. Despite their potential, developing programmable and modular methods for living systems to interface with OECTs has proven challenging. Here we describe hybrid OECTs containing the model electroactive bacterium Shewanella oneidensis that enable the transduction of biological computations to electrical responses. Specifically, we fabricated planar p-type OECTs and demonstrated that channel de-doping is driven by extracellular electron transfer (EET) from S. oneidensis. Leveraging this mechanistic understanding and our ability to control EET flux via transcriptional regulation, we used plasmid-based Boolean logic gates to translate biological computation into current changes within the OECT. Finally, we demonstrated EET-driven changes to OECT synaptic plasticity. This work enables fundamental EET studies and OECT-based biosensing and biocomputing systems with genetically controllable and modular design elements.

Repurposing the mammalian RNA-binding protein Musashi-1 as an allosteric translation repressor in bacteria

The RNA recognition motif (RRM) is the most common RNA-binding protein domain identified in nature. However, RRM-containing proteins are only prevalent in eukaryotic phyla, in which they play central regulatory roles. Here, we engineered an orthogonal post-transcriptional control system of gene expression in the bacterium Escherichia coli with the mammalian RNA-binding protein Musashi-1, which is a stem cell marker with neurodevelopmental role that contains two canonical RRMs. In the circuit, Musashi-1 is regulated transcriptionally and works as an allosteric translation repressor thanks to a specific interaction with the N-terminal coding region of a messenger RNA and its structural plasticity to respond to fatty acids. We fully characterized the genetic system at the population and single-cell levels showing a significant fold change in reporter expression, and the underlying molecular mechanism by assessing the in vitro binding kinetics and in vivo functionality of a series of RNA mutants. The dynamic response of the system was well recapitulated by a bottom-up mathematical model. Moreover, we applied the post-transcriptional mechanism engineered with Musashi-1 to specifically regulate a gene within an operon, implement combinatorial regulation, and reduce protein expression noise. This work illustrates how RRM-based regulation can be adapted to simple organisms, thereby adding a new regulatory layer in prokaryotes for translation control.

Enhancing circuit stability under growth feedback with supplementary repressive regulation

The field of synthetic biology and biosystems engineering increasingly acknowledges the need for a holistic design approach that incorporates circuit-host interactions into the design process. Engineered circuits are not isolated entities but inherently entwined with the dynamic host environment. One such circuit-host interaction, 'growth feedback', results when modifications in host growth patterns influence the operation of gene circuits. The growth-mediated effects can range from growth-dependent elevation in protein/mRNA dilution rate to changes in resource reallocation within the cell, which can lead to complete functional collapse in complex circuits. To achieve robust circuit performance, synthetic biologists employ a variety of control mechanisms to stabilize and insulate circuit behavior against growth changes. Here we propose a simple strategy by incorporating one repressive edge in a growth-sensitive bistable circuit. Through both simulation and in vitro experimentation, we demonstrate how this additional repressive node stabilizes protein levels and increases the robustness of a bistable circuit in response to growth feedback. We propose the incorporation of repressive links in gene circuits as a control strategy for desensitizing gene circuits against growth fluctuations.

Context-dependent redesign of robust synthetic gene circuits

Cells provide dynamic platforms for executing exogenous genetic programs in synthetic biology, resulting in highly context-dependent circuit performance. Recent years have seen an increasing interest in understanding the intricacies of circuit-host relationships, their influence on the synthetic bioengineering workflow, and in devising strategies to alleviate undesired effects. We provide an overview of how emerging circuit-host interactions, such as growth feedback and resource competition, impact both deterministic and stochastic circuit behaviors. We also emphasize control strategies for mitigating these unwanted effects. This review summarizes the latest advances and the current state of host-aware and resource-aware design of synthetic gene circuits.

The motility-matrix production switch in Bacillus subtilis-a modeling perspective

Phenotype switching can be triggered by external stimuli and by intrinsic stochasticity. Here, we focus on the motility-matrix production switch in Bacillus subtilis. We use modeling to describe the SinR-SlrR bistable switch and its regulation by SinI and to distinguish different sources of stochasticity. Our simulations indicate that intrinsic fluctuations in the synthesis of SinI are insufficient to drive spontaneous switching and suggest that switching is triggered by upstream noise from the Spo0A phosphorelay. IMPORTANCE The switch from motility to matrix production is the first step toward biofilm formation and, thus, to multicellular behavior in Bacillus subtilis. The transition is governed by a bistable switch based on the interplay of the regulators SinR and SlrR, while SinI transmits upstream signals to that switch. Quantitative modeling can be used to study the switching dynamics. Here, we build such a model step by step to describe the dynamics of the switch and its regulation and to study how spontaneous switching is triggered by upstream noise from the Spo0A phosphorelay.

Adaptive laboratory evolution and transcriptomics-guided engineering of Escherichia coli for increased isobutanol tolerance

As a renewable energy from biomass, isobutanol is considered as a promising alternative to fossil fuels. To biotechnologically produce isobutanol, strain development using industrial microbial hosts, such as Escherichia coli, has been conducted by introducing a heterologous isobutanol synthetic pathway. However, the toxicity of produced isobutanol inhibits cell growth, thereby restricting improvements in isobutanol titer, yield, and productivity. Therefore, the development of robust microbial strains tolerant to isobutanol is required. In this study, isobutanol-tolerant mutants were isolated from two E. coli parental strains, E. coli BL21(DE3) and MG1655(DE3), through adaptive laboratory evolution (ALE) under high isobutanol concentrations. Subsequently, 16 putative genes responsible for isobutanol tolerance were identified by transcriptomic analysis. When overexpressed in E. coli, four genes (fadB, dppC, acs, and csiD) conferred isobutanol tolerance. A fermentation study with a reverse engineered isobutanol-producing E. coli JK209 strain showed that fadB or dppC overexpression improved isobutanol titers by 1.5 times, compared to the control strain. Through coupling adaptive evolution with transcriptomic analysis, new genetic targets utilizable were identified as the basis for the development of an isobutanol-tolerant strain. Thus, these new findings will be helpful not only for a fundamental understanding of microbial isobutanol tolerance but also for facilitating industrially feasible isobutanol production.

A calibration and uncertainty quantification analysis of classical, fractional and multiscale logistic models of tumour growth

BACKGROUND AND OBJECTIVE

The validation of mathematical models of tumour growth is frequently hampered by the lack of sufficient experimental data, resulting in qualitative rather than quantitative studies. Recent approaches to this problem have attempted to extract information about tumour growth by integrating multiscale experimental measurements, such as longitudinal cell counts and gene expression data. In the present study, we investigated the performance of several mathematical models of tumour growth, including classical logistic, fractional and novel multiscale models, in terms of quantifying in-vitro tumour growth in the presence and absence of therapy. We further examined the effect of genes associated with changes in chemosensitivity in cell death rates.

METHODS

The multiscale expansion of logistic growth models was performed by coupling gene expression profiles to the cell death rates. State-of-the-art Bayesian inference, likelihood maximisation and uncertainty quantification techniques allowed a thorough evaluation of model performance.

RESULTS

The results suggest that the classical single-cell population model (SCPM) was the best fit for the untreated and low-dose treatment conditions, while the multiscale model with a cell death rate symmetric with the expression profile of OCT4 (Sym-SCPM) yielded the best fit for the high-dose treatment data. Further identifiability analysis showed that the multiscale model was both structurally and practically identifiable under the condition of known OCT4 expression profiles.

CONCLUSIONS

Overall, the present study demonstrates that model performance can be improved by incorporating multiscale measurements of tumour growth when high-dose treatment is involved.

Plasticity of growth laws tunes resource allocation strategies in bacteria

Bacteria like E. coli grow at vastly different rates on different substrates, however, the precise reason for this variability is poorly understood. Different growth rates have been attributed to 'nutrient quality', a key parameter in bacterial growth laws. However, it remains unclear to what extent nutrient quality is rooted in fundamental biochemical constraints like the energy content of nutrients, the protein cost required for their uptake and catabolism, or the capacity of the plasma membrane for nutrient transporters. Here, we show that while nutrient quality is indeed reflected in protein investment in substrate-specific transporters and enzymes, this is not a fundamental limitation on growth rate, at least for certain 'poor' substrates. We show that it is possible to turn mannose, one of the 'poorest' substrates of E. coli, into one of the 'best' substrates by reengineering chromosomal promoters of the mannose transporter and metabolic enzymes required for mannose degradation. This result falls in line with previous observations of more subtle growth rate improvement for many other carbon sources. However, we show that this faster growth rate comes at the cost of diverse cellular capabilities, reflected in longer lag phases, worse starvation survival and lower motility. We show that addition of cAMP to the medium can rescue these phenotypes but imposes a corresponding growth cost. Based on these data, we propose that nutrient quality is largely a self-determined, plastic property that can be modulated by the fraction of proteomic resources devoted to a specific substrate in the much larger proteome sector of catabolically activated genes. Rather than a fundamental biochemical limitation, nutrient quality reflects resource allocation decisions that are shaped by evolution in specific ecological niches and can be quickly adapted if necessary.

Evaluating DFHBI-Responsive RNA Light-Up Aptamers as Fluorescent Reporters for Gene Expression

Protein-based fluorescent reporters have been widely used to characterize and localize biological processes in living cells. However, these reporters may have certain drawbacks for some applications, such as transcription-based studies or biological interactions with fast dynamics. In this context, RNA nanotechnology has emerged as a promising alternative, suggesting the use of functional RNA molecules as transcriptional fluorescent reporters. RNA-based aptamers can bind to nonfluorescent small molecules to activate their fluorescence. However, their performance as reporters of gene expression in living cells has not been fully characterized, unlike protein-based reporters. Here, we investigate the performance of three RNA light-up aptamers─F30-2xdBroccoli, tRNA-Spinach, and Tornado Broccoli─as fluorescent reporters for gene expression in Escherichia coli and compare them to a protein reporter. We examine the activation range and effect on the cell growth of RNA light-up aptamers in time-course experiments and demonstrate that these aptamers are suitable transcriptional reporters over time. Using flow cytometry, we compare the variability at the single-cell level caused by the RNA fluorescent reporters and protein-based reporters. We found that the expression of RNA light-up aptamers produced higher variability in a population than that of their protein counterpart. Finally, we compare the dynamical behavior of these RNA light-up aptamers and protein-based reporters. We observed that RNA light-up aptamers might offer faster dynamics compared to a fluorescent protein in E. coli. The implementation of these transcriptional reporters may facilitate transcription-based studies, gain further insights into transcriptional processes, and expand the implementation of RNA-based circuits in bacterial cells.

Steering and controlling evolution - from bioengineering to fighting pathogens

Control interventions steer the evolution of molecules, viruses, microorganisms or other cells towards a desired outcome. Applications range from engineering biomolecules and synthetic organisms to drug, therapy and vaccine design against pathogens and cancer. In all these instances, a control system alters the eco-evolutionary trajectory of a target system, inducing new functions or suppressing escape evolution. Here, we synthesize the objectives, mechanisms and dynamics of eco-evolutionary control in different biological systems. We discuss how the control system learns and processes information about the target system by sensing or measuring, through adaptive evolution or computational prediction of future trajectories. This information flow distinguishes pre-emptive control strategies by humans from feedback control in biotic systems. We establish a cost-benefit calculus to gauge and optimize control protocols, highlighting the fundamental link between predictability of evolution and efficacy of pre-emptive control.

A universal mechanism of biomass density homeostasis via ribosomal counterions

In all growing cells, the cell envelope must expand in concert with cytoplasmic biomass to prevent lysis or molecular crowding. The complex cell wall of microbes and plants makes this challenge especially daunting and it unclear how cells achieve this coordination. Here, we uncover a striking linear increase of cytoplasmic pressure with growth rate in E. coli. Remarkably, despite this increase in turgor pressure with growth rate, cellular biomass density was constant across a wide range of growth rates. In contrast, perturbing pressure away from this scaling directly affected biomass density. A mathematical model, in which endopeptidase-mediated cell wall fluidization enables turgor pressure to set the pace of cellular volume expansion, not only explains these confounding observations, but makes several surprising quantitative predictions that we validated experimentally. The picture that emerges is that changes in turgor pressure across growth rates are mediated by counterions of ribosomal RNA. Profoundly, the coupling between rRNA and cytoplasmic pressure simultaneously coordinates cell wall expansion across growth rates and exerts homeostatic feedback control on biomass density. Because ribosome content universally scales with growth rate in fast growing cells, this universal mechanism may control cell wall biosynthesis in microbes and plants and drive the expansion of ribosome-addicted tumors that can exert substantial mechanical forces on their environment.

Motility mediates satellite formation in confined biofilms

Bacteria have spectacular survival capabilities and can spread in many, vastly different environments. For instance, when pathogenic bacteria infect a host, they expand by proliferation and squeezing through narrow pores and elastic matrices. However, the exact role of surface structures-important for biofilm formation and motility-and matrix density in colony expansion and morphogenesis is still largely unknown. Using confocal laser-scanning microscopy, we show how satellite colonies emerge around Escherichia coli colonies embedded in semi-dense hydrogel in controlled in vitro assays. Using knock-out mutants, we tested how extra-cellular structures, (e.g., exo-polysaccharides, flagella, and fimbria) control this morphology. Moreover, we identify the extra-cellular matrix' density, where this morphology is possible. When paralleled with mathematical modelling, our results suggest that satellite formation allows bacterial communities to spread faster. We anticipate that this strategy is important to speed up expansion in various environments, while retaining the close interactions and protection provided by the community.

Controlling circuitry underlies the growth optimization of Saccharomyces cerevisiae

Microbial growth emerges from coordinated synthesis of various cellular components from limited resources. In Saccharomyces cerevisiae, cyclic AMP (cAMP)-mediated signaling is shown to orchestrate cellular metabolism; however, it remains unclear quantitatively how the controlling circuit drives resource partition and subsequently shapes biomass growth. Here we combined experiment with mathematical modeling to dissect the signaling-mediated growth optimization of S. cerevisiae. We showed that, through cAMP-mediated control, the organism achieves maximal or nearly maximal steady-state growth during the utilization of multiple tested substrates as well as under perturbations impairing glucose uptake. However, the optimal cAMP concentration varies across cases, suggesting that different modes of resource allocation are adopted for varied conditions. Under settings with nutrient alterations, S. cerevisiae tunes its cAMP level to dynamically reprogram itself to realize rapid adaptation. Moreover, to achieve growth maximization, cells employ additional regulatory systems such as the GCN2-mediated amino acid control. This study establishes a systematic understanding of global resource allocation in S. cerevisiae, providing insights into quantitative yeast physiology as well as metabolic strain engineering for biotechnological applications.

Genomic footprinting uncovers global transcription factor responses to amino acids in Escherichia coli

Our knowledge of transcriptional responses to changes in nutrient availability comes primarily from few well-studied transcription factors (TFs), often lacking an unbiased genome-wide perspective. Leveraging recent advances allowing bacterial genomic footprinting, we comprehensively mapped the genome-wide regulatory responses of Escherichia coli to exogenous leucine, methionine, alanine, and lysine. The global TF Lrp was found to individually sense three amino acids and mount three different target gene responses. Overall, 531 genes had altered RNA polymerase occupancy, and 32 TFs responded directly or indirectly to the presence of amino acids, including regulators of membrane and osmotic pressure homeostasis. About 70% of the detected TF-DNA interactions had not been reported before. We thus identified 682 previously unknown TF-binding locations, for a subset of which the involved TFs were identified by affinity purification. This comprehensive map of amino acid regulation illustrates the incompleteness of the known transcriptional regulation network, even in E. coli.

Bacterial toxin-antitoxin systems: Novel insights on toxin activation across populations and experimental shortcomings

The alarming rise in hard-to-treat bacterial infections is of great concern to human health. Thus, the identification of molecular mechanisms that enable the survival and growth of pathogens is of utmost urgency for the development of more efficient antimicrobial therapies. In challenging environments, such as presence of antibiotics, or during host infection, metabolic adjustments are essential for microorganism survival and competitiveness. Toxin-antitoxin systems (TASs) consisting of a toxin with metabolic modulating activity and a cognate antitoxin that antagonizes that toxin are important elements in the arsenal of bacterial stress defense. However, the exact physiological function of TA systems is highly debatable and with the exception of stabilization of mobile genetic elements and phage inhibition, other proposed biological functions lack a broad consensus. This review aims at gaining new insights into the physiological effects of TASs in bacteria and exploring the experimental shortcomings that lead to discrepant results in TAS research. Distinct control mechanisms ensure that only subsets of cells within isogenic cultures transiently develop moderate levels of toxin activity. As a result, TASs cause phenotypic growth heterogeneity rather than cell stasis in the entire population. It is this feature that allows bacteria to thrive in diverse environments through the creation of subpopulations with different metabolic rates and stress tolerance programs.

Dynamical model of antibiotic responses linking expression of resistance to metabolism explains emergence of heterogeneity during drug exposures

Antibiotic responses in bacteria are highly dynamic and heterogeneous, with sudden exposure of bacterial colonies to high drug doses resulting in the coexistence of recovered and arrested cells. The dynamics of the response is determined by regulatory circuits controlling the expression of resistance genes, which are in turn modulated by the drug's action on cell growth and metabolism. Despite advances in understanding gene regulation at the molecular level, we still lack a framework to describe how feedback mechanisms resulting from the interdependence between expression of resistance and cell metabolism can amplify naturally occurring noise and create heterogeneity at the population level. To understand how this interplay affects cell survival upon exposure, we constructed a mathematical model of the dynamics of antibiotic responses that links metabolism and regulation of gene expression, based on the tetracycline resistance tet operon in E. coli. We use this model to interpret measurements of growth and expression of resistance in microfluidic experiments, both in single cells and in biofilms. We also implemented a stochastic model of the drug response, to show that exposure to high drug levels results in large variations of recovery times and heterogeneity at the population level. We show that stochasticity is important to determine how nutrient quality affects cell survival during exposure to high drug concentrations. A quantitative description of how microbes respond to antibiotics in dynamical environments is crucial to understand population-level behaviors such as biofilms and pathogenesis.

Plasticity of growth laws tunes resource allocation strategies in bacteria

Bacteria like E. coli grow at vastly different rates on different substrates, however, the precise reason for this variability is poorly understood. Different growth rates have been attributed to 'nutrient quality', a key parameter in bacterial growth laws. However, it remains unclear to what extent nutrient quality is rooted in fundamental biochemical constraints like the energy content of nutrients, the protein cost required for their uptake and catabolism, or the capacity of the plasma membrane for nutrient transporters. Here, we show that while nutrient quality is indeed reflected in protein investment in substrate-specific transporters and enzymes, this is not a fundamental limitation on growth rate. We show that it is possible to turn mannose, one of the 'poorest' substrates of E. coli, into one of the 'best' substrates by reengineering chromosomal promoters of the mannose transporter and metabolic enzymes required for mannose degradation. However, we show that this faster growth rate comes at the cost of diverse cellular capabilities, reflected in longer lag phases, worse starvation survival and lower motility. We show that addition of cAMP to the medium can rescue these phenotypes but imposes a corresponding growth cost. Based on these data, we propose that nutrient quality is largely a self-determined, plastic property that can be modulated by the fraction of proteomic resources devoted to a specific substrate in the much larger proteome sector of catabolically activated genes. Rather than a fundamental biochemical limitation, nutrient quality reflects resource allocation decisions that are shaped by evolution in specific ecological niches and can be quickly adapted if necessary.

Unbalanced response to growth variations reshapes the cell fate decision landscape

The global regulation of cell growth rate on gene expression perturbs the performance of gene networks, which would impose complex variations on the cell-fate decision landscape. Here we use a simple synthetic circuit of mutual repression that allows a bistable landscape to examine how such global regulation would affect the stability of phenotypic landscape and the accompanying dynamics of cell-fate determination. We show that the landscape experiences a growth-rate-induced bifurcation between monostability and bistability. Theoretical and experimental analyses reveal that this bifurcating deformation of landscape arises from the unbalanced response of gene expression to growth variations. The path of growth transition across the bifurcation would reshape cell-fate decisions. These results demonstrate the importance of growth regulation on cell-fate determination processes, regardless of specific molecular signaling or regulation.

Long-term history dependence of growth rates of E. coli after nutrient shifts

According to a widely accepted paradigm of microbiology, steady-state growth rates are determined solely by current growth conditions1-3 and adaptations between growth states are rapid, as recently recapitulated by simple resource allocation models4. However, even in microbes overlapping regulatory networks can yield multi-stability or long-term cellular memory. Species like Listeria monocytogenes5 and Bacillus subtilis "distinguish" distinct histories for the commitment to sporulation6, but it is unclear if these states can persist over many generations. Remarkably, studying carbon co-utilization of Escherichia coli, we found that growth rates on combinations of carbon sources can depend critically on the previous growth condition. Growing in identical conditions, we observed differences in growth rates of up to 25% and we did not observe convergence of growth rates over 15 generations. We observed this phenomenon occurs across combinations of different phosphotransferase (PTS) substrates with various gluconeogenic carbon sources and found it to depend on the transcription factor Mlc.

A Hybrid Transistor with Transcriptionally Controlled Computation and Plasticity

Organic electrochemical transistors (OECTs) are ideal devices for translating biological signals into electrical readouts and have applications in bioelectronics, biosensing, and neuromorphic computing. Despite their potential, developing programmable and modular methods for living systems to interface with OECTs has proven challenging. Here we describe hybrid OECTs containing the model electroactive bacterium Shewanella oneidensis that enable the transduction of biological computations to electrical responses. Specifically, we fabricated planar p-type OECTs and demonstrated that channel de-doping is driven by extracellular electron transfer (EET) from S. oneidensis. Leveraging this mechanistic understanding and our ability to control EET flux via transcriptional regulation, we used plasmid-based Boolean logic gates to translate biological computation into current changes within the OECT. Finally, we demonstrated EET-driven changes to OECT synaptic plasticity. This work enables fundamental EET studies and OECT-based biosensing and biocomputing systems with genetically controllable and modular design elements.

Bye bye, linearity, bye: quantification of the mean for linear CRNs in a random environment

Molecular reactions within a cell are inherently stochastic, and cells often differ in morphological properties or interact with a heterogeneous environment. Consequently, cell populations exhibit heterogeneity both due to these intrinsic and extrinsic causes. Although state-of-the-art studies that focus on dissecting this heterogeneity use single-cell measurements, the bulk data that shows only the mean expression levels is still in routine use. The fingerprint of the heterogeneity is present also in bulk data, despite being hidden from direct measurement. In particular, this heterogeneity can affect the mean expression levels via bimolecular interactions with low-abundant environment species. We make this statement rigorous for the class of linear reaction systems that are embedded in a discrete state Markov environment. The analytic expression that we provide for the stationary mean depends on the reaction rate constants of the linear subsystem, as well as the generator and stationary distribution of the Markov environment. We demonstrate the effect of the environment on the stationary mean. Namely, we show how the heterogeneous case deviates from the quasi-steady state (Q.SS) case when the embedded system is fast compared to the environment.

Out-of-equilibrium gene expression fluctuations in the presence of extrinsic noise

Cell-to-cell variability in protein concentrations is strongly affected by extrinsic noise, especially for highly expressed genes. Extrinsic noise can be due to fluctuations of several possible cellular factors connected to cell physiology and to the level of key enzymes in the expression process. However, how to identify the predominant sources of extrinsic noise in a biological system is still an open question. This work considers a general stochastic model of gene expression with extrinsic noise represented as fluctuations of the different model rates, and focuses on the out-of-equilibrium expression dynamics. Combining analytical calculations with stochastic simulations, we characterize how extrinsic noise shapes the protein variability during gene activation or inactivation, depending on the prevailing source of extrinsic variability, on its intensity and timescale. In particular, we show that qualitatively different noise profiles can be identified depending on which are the fluctuating parameters. This indicates an experimentally accessible way to pinpoint the dominant sources of extrinsic noise using time-coarse experiments.

Growth Feedback Confers Cooperativity in Resource-Competing Synthetic Gene Circuits

Modularity is a key concept in designing synthetic gene circuits, as it allows for constructing complex molecular systems using well-characterized building blocks. One of the major challenges in this field is that these modular components often do not function as expected when assembled into larger circuits. One of the major issues is caused by resource competition, where multiple genes in the circuit compete for the same limited cellular resources, such as transcription factors and ribosomes. In addition, the mutual inhibition between synthetic gene circuits and cell growth results in growth feedback that significantly impacts its host-circuit dynamics. However, the complexity of the gene circuit dynamics under intertwined resource competition and growth feedback is not fully understood. This study developed a theoretical framework to examine the dynamics of synthetic gene circuits by considering both growth feedback and resource competition. Our results suggest a cooperative behavior between resource-competing gene circuits under growth feedback. Cooperation or competition is non-monotonically determined by the metabolic burden threshold. These two diverse effects could lead to the activation or deactivation of one circuit by the other. Lastly, the cooperativity mediated by growth feedback can attenuate the winner-takes-all resource competition. These findings show that coupling growth feedback and resource competition plays a crucial role in the dynamics of the host-circuit system, and understanding its effects helps control unexpected gene expression behaviors.

Evolutionary barriers to horizontal gene transfer in macrophage-associated Salmonella

Horizontal gene transfer (HGT) is a powerful evolutionary force facilitating bacterial adaptation and emergence of novel phenotypes. Several factors, including environmental ones, are predicted to restrict HGT, but we lack systematic and experimental data supporting these predictions. Here, we address this gap by measuring the relative fitness of 44 genes horizontally transferred from Escherichia coli to Salmonella enterica in infection-relevant environments. We estimated the distribution of fitness effects in each environment and identified that dosage-dependent effects across different environments are a significant barrier to HGT. The majority of genes were found to be deleterious. We also found longer genes had stronger negative fitness consequences than shorter ones, showing that gene length was negatively associated with HGT. Furthermore, fitness effects of transferred genes were found to be environmentally dependent. In summary, a substantial fraction of transferred genes had a significant fitness cost on the recipient, with both gene characteristics and the environment acting as evolutionary barriers to HGT.

Effects of growth feedback on gene circuits: A dynamical understanding

The successful integration of engineered gene circuits into host cells remains a significant challenge in synthetic biology due to circuit-host interactions, such as growth feedback, where the circuit influences cell growth and vice versa. Understanding the dynamics of circuit failures and identifying topologies resilient to growth feedback are crucial for both fundamental and applied research. Utilizing transcriptional regulation circuits with adaptation as a paradigm, we systematically study 435 distinct topological structures and uncover six categories of failures. Three dynamical mechanisms of circuit failures are identified: continuous deformation of the response curve, strengthened or induced oscillations, and sudden switching to coexisting attractors. Our extensive computations also uncover a scaling law between a circuit robustness measure and the strength of growth feedback. Despite the negative effects of growth feedback on the majority of circuit topologies, we identify a few circuits that maintain optimal performance as designed, a feature important for applications.

Shaping bacterial gene expression by physiological and proteome allocation constraints

Networks of molecular regulators are often the primary objects of focus in the study of gene regulation, with the machinery of protein synthesis tacitly relegated to the background. Shifting focus to the constraints imposed by the allocation of protein synthesis flux reveals surprising ways in which the actions of molecular regulators are shaped by physiological demands. Using carbon catabolite repression as a case study, we describe how physiological constraints are sensed through metabolic fluxes and how flux-controlled regulation gives rise to simple empirical relations between protein levels and the rate of cell growth.

Emergent statistical laws in single-cell transcriptomic data

Large-scale data on single-cell gene expression have the potential to unravel the specific transcriptional programs of different cell types. The structure of these expression datasets suggests a similarity with several other complex systems that can be analogously described through the statistics of their basic building blocks. Transcriptomes of single cells are collections of messenger RNA abundances transcribed from a common set of genes just as books are different collections of words from a shared vocabulary, genomes of different species are specific compositions of genes belonging to evolutionary families, and ecological niches can be described by their species abundances. Following this analogy, we identify several emergent statistical laws in single-cell transcriptomic data closely similar to regularities found in linguistics, ecology, or genomics. A simple mathematical framework can be used to analyze the relations between different laws and the possible mechanisms behind their ubiquity. Importantly, treatable statistical models can be useful tools in transcriptomics to disentangle the actual biological variability from general statistical effects present in most component systems and from the consequences of the sampling process inherent to the experimental technique.

Improved Bacterial Single-Cell RNA-Seq through Automated MATQ-Seq and Cas9-Based Removal of rRNA Reads

Bulk RNA sequencing technologies have provided invaluable insights into host and bacterial gene expression and associated regulatory networks. Nevertheless, the majority of these approaches report average expression across cell populations, hiding the true underlying expression patterns that are often heterogeneous in nature. Due to technical advances, single-cell transcriptomics in bacteria has recently become reality, allowing exploration of these heterogeneous populations, which are often the result of environmental changes and stressors. In this work, we have improved our previously published bacterial single-cell RNA sequencing (scRNA-seq) protocol that is based on multiple annealing and deoxycytidine (dC) tailing-based quantitative scRNA-seq (MATQ-seq), achieving a higher throughput through the integration of automation. We also selected a more efficient reverse transcriptase, which led to reduced cell loss and higher workflow robustness. Moreover, we successfully implemented a Cas9-based rRNA depletion protocol into the MATQ-seq workflow. Applying our improved protocol on a large set of single Salmonella cells sampled over different growth conditions revealed improved gene coverage and a higher gene detection limit compared to our original protocol and allowed us to detect the expression of small regulatory RNAs, such as GcvB or CsrB at a single-cell level. In addition, we confirmed previously described phenotypic heterogeneity in Salmonella in regard to expression of pathogenicity-associated genes. Overall, the low percentage of cell loss and high gene detection limit makes the improved MATQ-seq protocol particularly well suited for studies with limited input material, such as analysis of small bacterial populations in host niches or intracellular bacteria. IMPORTANCE Gene expression heterogeneity among isogenic bacteria is linked to clinically relevant scenarios, like biofilm formation and antibiotic tolerance. The recent development of bacterial single-cell RNA sequencing (scRNA-seq) enables the study of cell-to-cell variability in bacterial populations and the mechanisms underlying these phenomena. Here, we report a scRNA-seq workflow based on MATQ-seq with increased robustness, reduced cell loss, and improved transcript capture rate and gene coverage. Use of a more efficient reverse transcriptase and the integration of an rRNA depletion step, which can be adapted to other bacterial single-cell workflows, was instrumental for these improvements. Applying the protocol to the foodborne pathogen Salmonella, we confirmed transcriptional heterogeneity across and within different growth phases and demonstrated that our workflow captures small regulatory RNAs at a single-cell level. Due to low cell loss and high transcript capture rates, this protocol is uniquely suited for experimental settings in which the starting material is limited, such as infected tissues.

The Slowdown of Growth Rate Controls the Single-Cell Distribution of Biofilm Matrix Production via an SinI-SinR-SlrR Network

In Bacillus subtilis, master regulator Spo0A controls several cell-differentiation pathways. Under moderate starvation, phosphorylated Spo0A (Spo0A~P) induces biofilm formation by indirectly activating genes controlling matrix production in a subpopulation of cells via an SinI-SinR-SlrR network. Under severe starvation, Spo0A~P induces sporulation by directly and indirectly regulating sporulation gene expression. However, what determines the heterogeneity of individual cell fates is not fully understood. In particular, it is still unclear why, despite being controlled by a single master regulator, biofilm matrix production and sporulation seem mutually exclusive on a single-cell level. In this work, with mathematical modeling, we showed that the fluctuations in the growth rate and the intrinsic noise amplified by the bistability in the SinI-SinR-SlrR network could explain the single-cell distribution of matrix production. Moreover, we predicted an incoherent feed-forward loop; the decrease in the cellular growth rate first activates matrix production by increasing in Spo0A phosphorylation level but then represses it via changing the relative concentrations of SinR and SlrR. Experimental data provide evidence to support model predictions. In particular, we demonstrate how the degree to which matrix production and sporulation appear mutually exclusive is affected by genetic perturbations. IMPORTANCE The mechanisms of cell-fate decisions are fundamental to our understanding of multicellular organisms and bacterial communities. However, even for the best-studied model systems we still lack a complete picture of how phenotypic heterogeneity of genetically identical cells is controlled. Here, using B. subtilis as a model system, we employ a combination of mathematical modeling and experiments to explain the population-level dynamics and single-cell level heterogeneity of matrix gene expression. The results demonstrate how the two cell fates, biofilm matrix production and sporulation, can appear mutually exclusive without explicitly inhibiting one another. Such a mechanism could be used in a wide range of other biological systems.

Development of Nanoparticle Adaptation Phenomena in Acinetobacter baumannii: Physiological Change and Defense Response

The present work describes the evolution of a resistance phenotype to a multitargeting antimicrobial agent, namely, silver nanoparticles (nanosilver; NAg), in the globally prevalent bacterial pathogen Acinetobacter baumannii. The Gram-negative bacterium has recently been listed as a critical priority pathogen requiring novel treatment options by the World Health Organization. Through prolonged exposure to the important antimicrobial nanoparticle, the bacterium developed mutations in genes that encode the protein subunits of organelle structures that are involved in cell-to-surface attachment as well as in a cell envelope capsular polysaccharide synthesis-related gene. These mutations are potentially correlated with stable physiological changes in the biofilm growth behavior and with an evident protective effect against oxidative stress, most likely as a feature of toxicity defense. We further report a different adaptation response of A. baumannii to the cationic form of silver (Ag+). The bacterium developed a tolerance phenotype to Ag+, which was correlated with an indicative surge in respiratory activity and changes in cell morphology, of which these are reported characteristics of tolerant bacterial populations. The findings regarding adaptation phenomena to NAg highlight the risks of the long-term use of the nanoparticle on a priority pathogen. The findings urge the implementation of strategies to overcome bacterial NAg adaptation, to better elucidate the toxicity mechanisms of the nanoparticle, and preserve the efficacy of the potent alternative antimicrobial agent in this era of antimicrobial resistance. IMPORTANCE Several recent studies have reported on the development of bacterial resistance to broad-spectrum antimicrobial silver nanoparticles (nanosilver; NAg). NAg is currently one of the most important alternative antimicrobial agents. However, no studies have yet established whether Acinetobacter baumannii, a globally prevalent nosocomial pathogen, can develop resistance to the nanoparticle. The study herein describes how a model strain of A. baumannii with no inherent silver resistance determinants developed resistance to NAg, following prolonged exposure. The stable physiological changes are correlated with mutations detected in the bacterium genome. These mutations render the bacterium capable of proliferating at a toxic NAg concentration. It was also found that A. baumannii developed a "slower-to-kill" tolerance trait to Ag+, which highlights the unique antimicrobial activities between the nanoparticulate and the ionic forms of silver. Despite the proven efficacy of NAg, the observation of NAg resistance in A. baumannii emphasises the potential risks of the repeated overuse of this agent on a priority pathogen.

A Systems Biology Approach To Disentangle the Direct and Indirect Effects of Global Transcription Factors on Gene Expression in Escherichia coli

Delineating the pleiotropic effects of global transcriptional factors (TFs) is critical for understanding the system-wide regulatory response in a particular environment. Currently, with the availability of genome-wide TF binding and gene expression data for Escherichia coli, several gene targets can be assigned to the global TFs, albeit inconsistently. Here, using a systematic integrated approach with emphasis on metabolism, we characterized and quantified the direct effects as well as the growth rate-mediated indirect effects of global TFs using deletion mutants of FNR, ArcA, and IHF regulators (focal TFs) under glucose fermentative conditions. This categorization enabled us to disentangle the dense connections seen within the transcriptional regulatory network (TRN) and determine the exact nature of focal TF-driven epistatic interactions with other global and pathway-specific local regulators (iTFs). We extended our analysis to combinatorial deletions of these focal TFs to determine their cross talk effects as well as conserved patterns of regulatory interactions. Moreover, we predicted with high confidence several novel metabolite-iTF interactions using inferred iTF activity changes arising from the allosteric effects of the intracellular metabolites perturbed as a result of the absence of focal TFs. Further, using compendium level computational analyses, we revealed not only the coexpressed genes regulated by these focal TFs but also the coordination of the direct and indirect target expression in the context of the economy of intracellular metabolites. Overall, this study leverages the fundamentals of TF-driven regulation, which could serve as a better template for deciphering mechanisms underlying complex phenotypes. IMPORTANCE Understanding the pleiotropic effects of global TFs on gene expression and their relevance underlying a specific response in a particular environment has been challenging. Here, we distinguish the TF-driven direct effects and growth rate-mediated indirect effects on gene expression using single- and double-deletion mutants of FNR, ArcA, and IHF regulators under anaerobic glucose fermentation. Such dissection assists us in unraveling the precise nature of interactions existing between the focal TF(s) and several other TFs, including those altered by allosteric effects of intracellular metabolites. We were able to recapitulate the previously known metabolite-TF interactions and predict novel interactions with high confidence. Furthermore, we determined that the direct and indirect gene expression have a strong connection with each other when analyzed using the coexpressed- or coregulated-gene approach. Deciphering such regulatory patterns explicitly from the metabolism point of view would be valuable in understanding other unpredicted complex regulation existing in nature.

Growth rate-associated transcriptome reorganization in response to genomic, environmental, and evolutionary interruptions

The genomic, environmental, and evolutionary interruptions caused the changes in bacterial growth, which were stringently associated with changes in gene expression. The growth and gene expression changes remained unclear in response to these interruptions that occurred combinative. As a pilot study, whether and how bacterial growth was affected by the individual and dual interruptions of genome reduction, environmental stress, and adaptive evolution were investigated. Growth assay showed that the presence of the environmental stressors, i.e., threonine and chloramphenicol, significantly decreased the growth rate of the wild-type Escherichia coli, whereas not that of the reduced genome. It indicated a canceling effect in bacterial growth due to the dual interruption of the genomic and environmental changes. Experimental evolution of the reduced genome released the canceling effect by improving growth fitness. Intriguingly, the transcriptome architecture maintained a homeostatic chromosomal periodicity regardless of the genomic, environmental, and evolutionary interruptions. Negative epistasis in transcriptome reorganization was commonly observed in response to the dual interruptions, which might contribute to the canceling effect. It was supported by the changes in the numbers of differentially expressed genes (DEGs) and the enriched regulons and functions. Gene network analysis newly constructed 11 gene modules, one out of which was correlated to the growth rate. Enrichment of DEGs in these modules successfully categorized them into three types, i.e., conserved, responsive, and epistatic. Taken together, homeostasis in transcriptome architecture was essential to being alive, and it might be attributed to the negative epistasis in transcriptome reorganization and the functional differentiation in gene modules. The present study directly connected bacterial growth fitness with transcriptome reorganization and provided a global view of how microorganisms responded to genomic, environmental, and evolutionary interruptions for survival from wild nature.

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.

Negatively Competitive Incoherent Feedforward Loops Mitigate Winner-Take-All Resource Competition

The effects of host resource limitations on the function of synthetic gene circuits have gained significant attention over the past years. Hosts, having evolved resource capacities optimal for their own genome, have been repeatedly demonstrated to suffer from the added burden of synthetic genetic programs, which may in return pose deleterious effects on the circuit's function. Three resource controller archetypes have been proposed previously to mitigate resource distribution problems in dynamic circuits: the local controller, the global controller, and a "negatively competitive" regulatory (NCR) controller that utilizes synthetic competition to combat resource competition. The dynamics of negative feedback forms of these controllers have been previously investigated, and here we extend the analysis of these resource allocation strategies to the incoherent feedforward loop (iFFL) topology. We demonstrate that the three iFFL controllers can attenuate Winner-Take-All resource competition between two bistable switches. We uncover that the parameters associated with the synthetic competition in the NCR iFFL controller are paramount to its increased efficacy over the local controller type, while the global controllers demonstrate to be relatively ineffectual. Interestingly, unlike the negative feedback counterpart topologies, iFFL controllers exhibit a unique coupling of switch activation thresholds which we term the "coactivation threshold shift" effect. Finally, we demonstrate that a nearly fully orthogonal set of bistable switches could be achieved by pairing an NCR controller with an appropriate level of controller resource consumption.

Principles of gene regulation quantitatively connect DNA to RNA and proteins in bacteria

Protein concentrations are set by a complex interplay between gene-specific regulatory processes and systemic factors, including cell volume and shared gene expression machineries. Elucidating this interplay is crucial for discerning and designing gene regulatory systems. We quantitatively characterized gene-specific and systemic factors that affect transcription and translation genome-wide for Escherichia coli across many conditions. The results revealed two design principles that make regulation of gene expression insulated from concentrations of shared machineries: RNA polymerase activity is fine-tuned to match translational output, and translational characteristics are similar across most messenger RNAs (mRNAs). Consequently, in bacteria, protein concentration is set primarily at the promoter level. A simple mathematical formula relates promoter activities and protein concentrations across growth conditions, enabling quantitative inference of gene regulation from omics data.

Transcriptional Fluctuations Govern the Serum-Dependent Cell Cycle Duration Heterogeneities in Mammalian Cells

Mammalian cells exhibit a high degree of intercellular variability in cell cycle period and phase durations. However, the factors orchestrating the cell cycle duration heterogeneities remain unclear. Herein, by combining cell cycle network-based mathematical models with live single-cell imaging studies under varied serum conditions, we demonstrate that fluctuating transcription rates of cell cycle regulatory genes across cell lineages and during cell cycle progression in mammalian cells majorly govern the robust correlation patterns of cell cycle period and phase durations among sister, cousin, and mother-daughter lineage pairs. However, for the overall cellular population, alteration in the serum level modulates the fluctuation and correlation patterns of cell cycle period and phase durations in a correlated manner. These heterogeneities at the population level can be fine-tuned under limited serum conditions by perturbing the cell cycle network using a p38-signaling inhibitor without affecting the robust lineage-level correlations. Overall, our approach identifies transcriptional fluctuations as the key controlling factor for the cell cycle duration heterogeneities and predicts ways to reduce cell-to-cell variabilities by perturbing the cell cycle network regulations.

Robust replication initiation from coupled homeostatic mechanisms

The bacterium Escherichia coli initiates replication once per cell cycle at a precise volume per origin and adds an on average constant volume between successive initiation events, independent of the initiation size. Yet, a molecular model that can explain these observations has been lacking. Experiments indicate that E. coli controls replication initiation via titration and activation of the initiator protein DnaA. Here, we study by mathematical modelling how these two mechanisms interact to generate robust replication-initiation cycles. We first show that a mechanism solely based on titration generates stable replication cycles at low growth rates, but inevitably causes premature reinitiation events at higher growth rates. In this regime, the DnaA activation switch becomes essential for stable replication initiation. Conversely, while the activation switch alone yields robust rhythms at high growth rates, titration can strongly enhance the stability of the switch at low growth rates. Our analysis thus predicts that both mechanisms together drive robust replication cycles at all growth rates. In addition, it reveals how an origin-density sensor yields adder correlations.

Predicting the re-distribution of antibiotic molecules caused by inter-species interactions in microbial communities

Microbes associate in nature forming complex communities, but they are often studied in purified form. Here I show that neighbouring species enforce the re-distribution of carbon and antimicrobial molecules, predictably changing drug efficacy with respect to standard laboratory assays. A simple mathematical model, validated experimentally using pairwise competition assays, suggests that differences in drug sensitivity between the competing species causes the re-distribution of drug molecules without affecting carbon uptake. The re-distribution of drug is even when species have similar drug sensitivity, reducing drug efficacy. But when their sensitivities differ the re-distribution is uneven: The most sensitive species accumulates more drug molecules, increasing efficacy against it. Drug efficacy tests relying on samples with multiple species are considered unreliable and unpredictable, but study demonstrates that efficacy in these cases can be qualitatively predicted. It also suggests that living in communities can be beneficial even when all species compete for a single carbon source, as the relationship between cell density and drug required to inhibit their growth may be more complex than previously thought.