A System for Gene Expression Noise Control in Yeast

Enhanced metabolic entanglement emerges during the evolution of an interkingdom microbial community

While different stages of mutualism can be observed in natural communities, the dynamics and mechanisms underlying the gradual erosion of independence of the initially autonomous organisms are not yet fully understood. In this study, by conducting the laboratory evolution on an engineered microbial community, we reproduce and molecularly track the stepwise progression towards enhanced partner entanglement. We observe that the evolution of the community both strengthens the existing metabolic interactions and leads to the emergence of de novo interdependence between partners for nitrogen metabolism, which is a common feature of natural symbiotic interactions. Selection for enhanced metabolic entanglement during the community evolution repeatedly occurred indirectly, via pleiotropies and trade-offs within cellular regulatory networks, and with no evidence of group selection. The indirect positive selection of metabolic dependencies between microbial community members, which results from the direct selection of other coupled traits in the same regulatory network, may therefore be a common but underappreciated driving force guiding the evolution of natural mutualistic communities.

Synthetic gene circuit evolution: Insights and opportunities at the mid-scale

Directed evolution focuses on optimizing single genetic components for predefined engineering goals by artificial mutagenesis and selection. In contrast, experimental evolution studies the adaptation of entire genomes in serially propagated cell populations, to provide an experimental basis for evolutionary theory. There is a relatively unexplored gap at the middle ground between these two techniques, to evolve in vivo entire synthetic gene circuits with nontrivial dynamic function instead of single parts or whole genomes. We discuss the requirements for such mid-scale evolution, with hypothetical examples for evolving synthetic gene circuits by appropriate selection and targeted shuffling of a seed set of genetic components in vivo. Implementing similar methods should aid the rapid generation, functionalization, and optimization of synthetic gene circuits in various organisms and environments, accelerating both the development of biomedical and technological applications and the understanding of principles guiding regulatory network evolution.

Deep-learning-assisted Sort-Seq enables high-throughput profiling of gene expression characteristics with high precision

Owing to the nondeterministic and nonlinear nature of gene expression, the steady-state intracellular protein abundance of a clonal population forms a distribution. The characteristics of this distribution, including expression strength and noise, are closely related to cellular behavior. However, quantitative description of these characteristics has so far relied on arrayed methods, which are time-consuming and labor-intensive. To address this issue, we propose a deep-learning-assisted Sort-Seq approach (dSort-Seq) in this work, enabling high-throughput profiling of expression properties with high precision. We demonstrated the validity of dSort-Seq for large-scale assaying of the dose-response relationships of biosensors. In addition, we comprehensively investigated the contribution of transcription and translation to noise production in Escherichia coli, from which we found that the expression noise is strongly coupled with the mean expression level. We also found that the transcriptional interference caused by overlapping RpoD-binding sites contributes to noise production, which suggested the existence of a simple and feasible noise control strategy in E. coli.

Nearly maximal information gain due to time integration in central dogma reactions

Living cells process information about their environment through the central dogma processes of transcription and translation, which drive the cellular response to stimuli. Here, we study the transfer of information from environmental input to the transcript and protein expression levels. Evaluation of both experimental and analogous simulation data reveals that transcription and translation are not two simple information channels connected in series. Instead, we demonstrate that the central dogma reactions often create a time-integrating information channel, where the translation channel receives and integrates multiple outputs from the transcription channel. This information channel model of the central dogma provides new information-theoretic selection criteria for the central dogma rate constants. Using the data for four well-studied species we show that their central dogma rate constants achieve information gain because of time integration while also keeping the loss because of stochasticity in translation relatively low (<0.5 bits).

Single-cell measurement quality in bits

Single-cell measurements have revolutionized our understanding of heterogeneity in cellular response. However, there is no universally comparable way to assess single-cell measurement quality. Here, we show how information theory can be used to assess and compare single-cell measurement quality in bits, which provides a universally comparable metric for information content. We anticipate that the experimental and theoretical approaches we show here will generally enable comparisons of quality between any single-cell measurement methods.

Heterosis of fitness and phenotypic variance in the evolution of a diploid gene regulatory network

Heterosis describes the phenomenon, whereby a hybrid population has higher fitness than an inbred population, which has previously been explained by either Mendelian dominance or overdominance under the general assumption of a simple genotype-phenotype relationship. However, recent studies have demonstrated that genes interact through a complex gene regulatory network (GRN). Furthermore, phenotypic variance is reportedly lower for heterozygotes, and the origin of such variance-related heterosis remains elusive. Therefore, a theoretical analysis linking heterosis to GRN evolution and stochastic gene expression dynamics is required. Here, we investigated heterosis related to fitness and phenotypic variance in a system with interacting genes by numerically evolving diploid GRNs. According to the results, the heterozygote population exhibited higher fitness than the homozygote population, indicating fitness-related heterosis resulting from evolution. In addition, the heterozygote population exhibited lower noise-related phenotypic variance in expression levels than the homozygous population, implying that the heterozygote population is more robust to noise. Furthermore, the distribution of the ratio of heterozygote phenotypic variance to homozygote phenotypic variance exhibited quantitative similarity with previous experimental results. By applying dominance and differential gene expression rather than only a single gene expression model, we confirmed the correlation between heterosis and differential gene expression. We explain our results by proposing that the convex high-fitness region is evolutionarily shaped in the genetic space to gain noise robustness under genetic mixing through sexual reproduction. These results provide new insights into the effects of GRNs on variance-related heterosis and differential gene expression.

Independent control of mean and noise by convolution of gene expression distributions

Gene expression noise can reduce cellular fitness or facilitate processes such as alternative metabolism, antibiotic resistance, and differentiation. Unfortunately, efforts to study the impacts of noise have been hampered by a scaling relationship between noise and expression level from individual promoters. Here, we use theory to demonstrate that mean and noise can be controlled independently by expressing two copies of a gene from separate inducible promoters in the same cell. We engineer low and high noise inducible promoters to validate this result in Escherichia coli, and develop a model that predicts the experimental distributions. Finally, we use our method to reveal that the response of a promoter to a repressor is less sensitive with higher repressor noise and explain this result using a law from probability theory. Our approach can be applied to investigate the effects of noise on diverse biological pathways or program cellular heterogeneity for synthetic biology applications.

Synthetic biology applications of the yeast mating signal pathway

Cell fusion is a fundamental biological process that is involved in the development of most eukaryotic organisms. During the fusion process in Saccharomyces cerevisiae, cells respond to pheromones to trigger the MAPK (mitogen-activated protein kinase) cascade to initiate mating, followed by polarization, cell-wall remodeling, membrane fusion, and karyogamy. We highlight the applications of the yeast mating signal pathway in promoter engineering for tuning the expression of output genes, as well as in metabolic engineering for decoupling growth and metabolism, biosensors for sensitive detection and signal amplification, genetic circuits for programmable biological functionalities, and artificial consortia for cell-cell communication. Strategies such as exploiting rational engineering of modular circuits and optimizing the reproductive pathway to precisely maneuver physiological events have implications for scientific research and industrial development.

Comparison of bias and resolvability in single-cell and single-transcript methods

Single-cell and single-transcript measurement methods have elevated our ability to understand and engineer biological systems. However, defining and comparing performance between methods remains a challenge, in part due to the confounding effects of experimental variability. Here, we propose a generalizable framework for performing multiple methods in parallel using split samples, so that experimental variability is shared between methods. We demonstrate the utility of this framework by performing 12 different methods in parallel to measure the same underlying reference system for cellular response. We compare method performance using quantitative evaluations of bias and resolvability. We attribute differences in method performance to steps along the measurement process such as sample preparation, signal detection, and choice of measurand. Finally, we demonstrate how this framework can be used to benchmark different methods for single-transcript detection. The framework we present here provides a practical way to compare performance of any methods.

In Situ Biomanufacturing of Small Molecules in the Mammalian Gut by Probiotic Saccharomyces boulardii

Saccharomyces boulardii is a probiotic yeast that exhibits rapid growth at 37 °C, is easy to transform, and can produce therapeutic proteins in the gut. To establish its ability to produce small molecules encoded by multigene pathways, we measured the amount and variance in protein expression enabled by promoters, terminators, selective markers, and copy number control elements. We next demonstrated efficient (>95%) CRISPR-mediated genome editing in this strain, allowing us to probe engineered gene expression across different genomic sites. We leveraged these strategies to assemble pathways enabling a wide range of vitamin precursor (β-carotene) and drug (violacein) titers. We found that S. boulardii colonizes germ-free mice stably for over 30 days and competes for niche space with commensal microbes, exhibiting short (1-2 day) gut residence times in conventional and antibiotic-treated mice. Using these tools, we enabled β-carotene synthesis (194 μg total) in the germ-free mouse gut over 14 days, estimating that the total mass of additional β-carotene recovered in feces was 56-fold higher than the β-carotene present in the initial probiotic dose. This work quantifies heterologous small molecule production titers by S. boulardii living in the mammalian gut and provides a set of tools for modulating these titers.

Orthogonal control of mean and variability of endogenous genes in a human cell line

Stochastic fluctuations at the transcriptional level contribute to isogenic cell-to-cell heterogeneity in mammalian cell populations. However, we still have no clear understanding of the repercussions of this heterogeneity, given the lack of tools to independently control mean expression and variability of a gene. Here, we engineer a synthetic circuit to modulate mean expression and heterogeneity of transgenes and endogenous human genes. The circuit, a Tunable Noise Rheostat (TuNR), consists of a transcriptional cascade of two inducible transcriptional activators, where the output mean and variance can be modulated by two orthogonal small molecule inputs. In this fashion, different combinations of the inputs can achieve the same mean but with different population variability. With TuNR, we achieve low basal expression, over 1000-fold expression of a transgene product, and up to 7-fold induction of the endogenous gene NGFR. Importantly, for the same mean expression level, we are able to establish varying degrees of heterogeneity in expression within an isogenic population, thereby decoupling gene expression noise from its mean. TuNR is therefore a modular tool that can be used in mammalian cells to enable direct interrogation of the implications of cell-to-cell variability.

Noise and the molecular processes underlying cell fate decision-making

Cell fate decision-making events involve the interplay of many molecular processes, ranging from signal transduction to genetic regulation, as well as a set of molecular and physiological feedback loops. Each aspect offers a rich field of investigation in its own right, but to understand the whole process, even in simple terms, we need to consider them together. Here we attempt to characterise this process by focussing on the roles of noise during cell fate decisions. We use a range of recent results to develop a view of the sequence of events by which a cell progresses from a pluripotent or multipotent to a differentiated state: chromatin organisation, transcription factor stoichiometry, and cellular signalling all change during this progression, and all shape cellular variability, which becomes maximal at the transition state.

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.

Orthogonal tuning of gene expression noise using CRISPR-Cas

The control of gene expression noise is important for improving drug treatment and the performance of synthetic biological systems. Previous work has tuned gene expression noise by changing the rate of transcription initiation, mRNA degradation, and mRNA translation. However, these methods are invasive: they require changes to the target genetic components. Here, we create an orthogonal system based on CRISPR-dCas9 to tune gene expression noise. Specifically, we modulate the gene expression noise of a reporter gene in Escherichia coli by incorporating CRISPR activation and repression (CRISPRar) simultaneously in a single cell. The CRISPRar uses a single dCas9 that recognizes two different single guide RNAs (sgRNA). We build a library of sgRNA variants with different expression activation and repression strengths. We find that expression noise and mean of a reporter gene can be tuned independently by CRISPRar. Our results suggest that the expression noise is tuned by the competition between two sgRNAs that modulate the binding of RNA polymerase to promoters. The CRISPRar may change how we tune expression noise at the genomic level. Our work has broad impacts on the study of gene functions, phenotypical heterogeneity, and genetic circuit control.

From Microbial Communities to Distributed Computing Systems

A distributed biological system can be defined as a system whose components are located in different subpopulations, which communicate and coordinate their actions through interpopulation messages and interactions. We see that distributed systems are pervasive in nature, performing computation across all scales, from microbial communities to a flock of birds. We often observe that information processing within communities exhibits a complexity far greater than any single organism. Synthetic biology is an area of research which aims to design and build synthetic biological machines from biological parts to perform a defined function, in a manner similar to the engineering disciplines. However, the field has reached a bottleneck in the complexity of the genetic networks that we can implement using monocultures, facing constraints from metabolic burden and genetic interference. This makes building distributed biological systems an attractive prospect for synthetic biology that would alleviate these constraints and allow us to expand the applications of our systems into areas including complex biosensing and diagnostic tools, bioprocess control and the monitoring of industrial processes. In this review we will discuss the fundamental limitations we face when engineering functionality with a monoculture, and the key areas where distributed systems can provide an advantage. We cite evidence from natural systems that support arguments in favor of distributed systems to overcome the limitations of monocultures. Following this we conduct a comprehensive overview of the synthetic communities that have been built to date, and the components that have been used. The potential computational capabilities of communities are discussed, along with some of the applications that these will be useful for. We discuss some of the challenges with building co-cultures, including the problem of competitive exclusion and maintenance of desired community composition. Finally, we assess computational frameworks currently available to aide in the design of microbial communities and identify areas where we lack the necessary tools.

Design of a MAPK signalling cascade balances energetic cost versus accuracy of information transmission

Cellular processes are inherently noisy, and the selection for accurate responses in presence of noise has likely shaped signalling networks. Here, we investigate the trade-off between accuracy of information transmission and its energetic cost for a mitogen-activated protein kinase (MAPK) signalling cascade. Our analysis of the pheromone response pathway of budding yeast suggests that dose-dependent induction of the negative transcriptional feedbacks in this network maximizes the information per unit energetic cost, rather than the information transmission capacity itself. We further demonstrate that futile cycling of MAPK phosphorylation and dephosphorylation has a measurable effect on growth fitness, with energy dissipation within the signalling cascade thus likely being subject to evolutionary selection. Considering optimization of accuracy versus the energetic cost of information processing, a concept well established in physics and engineering, may thus offer a general framework to understand the regulatory design of cellular signalling systems.

Cellular signalling networks provide information to the cell, but the trade-off between accuracy of information transfer and energetic cost of doing so has not been assessed. Here, the authors investigate a MAPK signalling cascade in budding yeast and find that information is maximised per unit energetic cost.

Application of combinatorial optimization strategies in synthetic biology

In the first wave of synthetic biology, genetic elements, combined into simple circuits, are used to control individual cellular functions. In the second wave of synthetic biology, the simple circuits, combined into complex circuits, form systems-level functions. However, efforts to construct complex circuits are often impeded by our limited knowledge of the optimal combination of individual circuits. For example, a fundamental question in most metabolic engineering projects is the optimal level of enzymes for maximizing the output. To address this point, combinatorial optimization approaches have been established, allowing automatic optimization without prior knowledge of the best combination of expression levels of individual genes. This review focuses on current combinatorial optimization methods and emerging technologies facilitating their applications.

Modeling the tunability of the dual-feedback genetic oscillator

Oscillatory gene circuits are ubiquitous to biology and are involved in fundamental processes of cell cycle, circadian rhythms, and developmental systems. The synthesis of small, non-natural oscillatory genetic circuits has been increasingly used to test the fundamental principles of genetic network dynamics. While the "repressilator" was used to first demonstrate the proof of principle, a more recently developed dual-feedback, fast, tunable genetic oscillator has demonstrated a greater degree of robustness and control over oscillatory behavior by combining positive- and negative-feedback loops. This oscillator, combining lacI (negative-) and araC (positive-) feedback loops, was, however, modeled using multiple layers of differential equations to capture the molecular complexity of regulation, in order to explain the experimentally measured oscillations. In the search for design principles of such minimal oscillatory circuits, we have developed a reduced model of this dual-feedback loop oscillator consisting of just six differential equations, two of which are delay differential equations. The delay term is optimized, as the only free parameter, to fit the experimental dynamics of the oscillator period and amplitude tunability by the two inducers isopropyl β-D-1-thiogalactopyranoside (IPTG) and arabinose. We proceed to use our reduced and experimentally validated model to redesign the network by comparing the effect of asymmetry in gene expression at the level of (a) DNA copy numbers and the rates of (b) mRNA translation and (c) degradation, since experimental and theoretical work had predicted a need for an asymmetry in the copy numbers of activator (araC) and repressor (lacI) genes encoded on plasmids. We confirm that the minimal period of the oscillator is sensitive to DNA copy number asymmetry, and can demonstrate that while the asymmetry in the translation rate has an identical effect as the plasmid copy numbers, modulating the asymmetry in mRNA degradation can improve the tunability of the period and amplitude of the oscillator. Thus, our model predicts control at the level of translation can be used to redesign such networks, for improved tunability, while at the same time making the network robust to replication "noise" and the effects of the host cell cycle. Thus, our model predicts experimentally testable principles to redesign a potentially more robust oscillatory genetic network.

Gene expression noise in a complex artificial toxin expression system

Gene expression is an intrinsically stochastic process. Fluctuations in transcription and translation lead to cell-to-cell variations in mRNA and protein levels affecting cellular function and cell fate. Here, using fluorescence time-lapse microscopy, we quantify noise dynamics in an artificial operon in Escherichia coli, which is based on the native operon of ColicinE2, a toxin. In the natural system, toxin expression is controlled by a complex regulatory network; upon induction of the bacterial SOS response, ColicinE2 is produced (cea gene) and released (cel gene) by cell lysis. Using this ColicinE2-based operon, we demonstrate that upon induction of the SOS response noise of cells expressing the operon is significantly lower for the (mainly) transcriptionally regulated gene cea compared to the additionally post-transcriptionally regulated gene cel. Likewise, we find that mutations affecting the transcriptional regulation by the repressor LexA do not significantly alter the population noise, whereas specific mutations to post-transcriptionally regulating units, strongly influence noise levels of both genes. Furthermore, our data indicate that global factors, such as the plasmid copy number of the operon encoding plasmid, affect gene expression noise of the entire operon. Taken together, our results provide insights on how noise in a native toxin-producing operon is controlled and underline the importance of post-transcriptional regulation for noise control in this system.

Controlling cell-to-cell variability with synthetic gene circuits

Cell-to-cell variability originating, for example, from the intrinsic stochasticity of gene expression, presents challenges for designing synthetic gene circuits that perform robustly. Conversely, synthetic biology approaches are instrumental in uncovering mechanisms underlying variability in natural systems. With a focus on reducing noise in individual genes, the field has established a broad synthetic toolset. This includes noise control by engineering of transcription and translation mechanisms either individually, or in combination to achieve independent regulation of mean expression and its variability. Synthetic feedback circuits use these components to establish more robust operation in closed-loop, either by drawing on, but also by extending traditional engineering concepts. In this perspective, we argue that major conceptual advances will require new theory of control adapted to biology, extensions from single genes to networks, more systematic considerations of origins of variability other than intrinsic noise, and an exploration of how noise shaping, instead of noise reduction, could establish new synthetic functions or help understanding natural functions.