Interpretation of morphogen gradients by a synthetic bistable circuit

A three-node Turing gene circuit forms periodic spatial patterns in bacteria

Turing patterns are self-organizing systems that can form spots, stripes, or labyrinths. Proposed examples in tissue organization include zebrafish pigmentation, digit spacing, and many others. The theory of Turing patterns in biology has been debated because of their stringent fine-tuning requirements, where patterns only occur within a small subset of parameters. This has complicated the engineering of synthetic Turing gene circuits from first principles, although natural genetic Turing networks have been identified. Here, we engineered a synthetic genetic reaction-diffusion system where three nodes interact according to a non-classical Turing network with improved parametric robustness. The system reproducibly generated stationary, periodic, concentric stripe patterns in growing E. coli colonies. A partial differential equation model reproduced the patterns, with a Turing parameter regime obtained by fitting to experimental data. Our synthetic Turing system can contribute to nanotechnologies, such as patterned biomaterial deposition, and provide insights into developmental patterning programs. A record of this paper's transparent peer review process is included in the supplemental information.

Pattern Formation and Bistability in a Synthetic Intercellular Genetic Toggle

Differentiation within multicellular organisms is a complex process that helps to establish spatial patterning and tissue formation within the body. Often, the differentiation of cells is governed by morphogens and intercellular signaling molecules that guide the fate of each cell, frequently using toggle-like regulatory components. Synthetic biologists have long sought to recapitulate patterned differentiation with engineered cellular communities, and various methods for differentiating bacteria have been invented. Here, we couple a synthetic corepressive toggle switch with intercellular signaling pathways to create a "quorum-sensing toggle". We show that this circuit not only exhibits population-wide bistability in a well-mixed liquid environment but also generates patterns of differentiation in colonies grown on agar containing an externally supplied morphogen. If coupled to other metabolic processes, circuits such as the one described here would allow for the engineering of spatially patterned, differentiated bacteria for use in biomaterials and bioelectronics.

Two orthogonal differentiation gradients locally coordinate fruit morphogenesis

Morphogenesis requires the coordination of cellular behaviors along developmental axes. In plants, gradients of growth and differentiation are typically established along a single longitudinal primordium axis to control global organ shape. Yet, it remains unclear how these gradients are locally adjusted to regulate the formation of complex organs that consist of diverse tissue types. Here we combine quantitative live imaging at cellular resolution with genetics, and chemical treatments to understand the formation of Arabidopsis thaliana female reproductive organ (gynoecium). We show that, contrary to other aerial organs, gynoecium shape is determined by two orthogonal, time-shifted differentiation gradients. An early mediolateral gradient controls valve morphogenesis while a late, longitudinal gradient regulates style differentiation. Local, tissue-dependent action of these gradients serves to fine-tune the common developmental program governing organ morphogenesis to ensure the specialized function of the gynoecium.

Identifying LasR Quorum Sensors with Improved Signal Specificity by Mapping the Sequence-Function Landscape

Programmable intercellular signaling using components of naturally occurring quorum sensing can allow for coordinated functions to be engineered in microbial consortia. LuxR-type transcriptional regulators are widely used for this purpose and are activated by homoserine lactone (HSL) signals. However, they often suffer from imperfect molecular discrimination of structurally similar HSLs, causing misregulation within engineered consortia containing multiple HSL signals. Here, we studied one such example, the regulator LasR from Pseudomonas aeruginosa. We elucidated its sequence-function relationship for ligand specificity using targeted protein engineering and multiplexed high-throughput biosensor screening. A pooled combinatorial saturation mutagenesis library (9,486 LasR DNA sequences) was created by mutating six residues in LasR's β5 sheet with single, double, or triple amino acid substitutions. Sort-seq assays were performed in parallel using cognate and noncognate HSLs to quantify each corresponding sensor's response to each HSL signal, which identified hundreds of highly specific variants. Sensor variants identified were individually assayed and exhibited up to 60.6-fold (p = 0.0013) improved relative activation by the cognate signal compared to the wildtype. Interestingly, we uncovered prevalent mutational epistasis and previously unidentified residues contributing to signal specificity. The resulting sensors with negligible signal crosstalk could be broadly applied to engineer bacteria consortia.

Synthetic symmetry breaking and programmable multicellular structure formation

During development, cells undergo symmetry breaking into differentiated subpopulations that self-organize into complex structures.1,2,3,4,5 However, few tools exist to recapitulate these behaviors in a controllable and coupled manner.6,7,8,9 Here, we engineer a stochastic recombinase genetic switch tunable by small molecules to induce programmable symmetry breaking, commitment to downstream cell fates, and morphological self-organization. Inducers determine commitment probabilities, generating tunable subpopulations as a function of inducer dosage. We use this switch to control the cell-cell adhesion properties of cells committed to each fate.10,11 We generate a wide variety of 3D morphologies from a monoclonal population and develop a computational model showing high concordance with experimental results, yielding new quantitative insights into the relationship between cell-cell adhesion strengths and downstream morphologies. We expect that programmable symmetry breaking, generating precise and tunable subpopulation ratios and coupled to structure formation, will serve as an integral component of the toolbox for complex tissue and organoid engineering.

Microbial models of development: Inspiration for engineering self-assembled synthetic multicellularity

While the field of synthetic developmental biology has traditionally focused on the study of the rich developmental processes seen in metazoan systems, an attractive alternate source of inspiration comes from microbial developmental models. Microbes face unique lifestyle challenges when forming emergent multicellular collectives. As a result, the solutions they employ can inspire the design of novel multicellular systems. In this review, we dissect the strategies employed in multicellular development by two model microbial systems: the cellular slime mold Dictyostelium discoideum and the biofilm-forming bacterium Bacillus subtilis. Both microbes face similar challenges but often have different solutions, both from metazoan systems and from each other, to create emergent multicellularity. These challenges include assembling and sustaining a critical mass of participating individuals to support development, regulating entry into development, and assigning cell fates. The mechanisms these microbial systems exploit to robustly coordinate development under a wide range of conditions offer inspiration for a new toolbox of solutions to the synthetic development community. Additionally, recreating these phenomena synthetically offers a pathway to understanding the key principles underlying how these behaviors are coordinated naturally.

The Mutual Inhibition between PLETHORAs and ARABIDOPSIS RESPONSE REGULATORs Controls Root Zonation

During organogenesis, a key step toward the development of a functional organ is the separation of cells into specific domains with different activities. Mutual inhibition of gene expression has been shown to be sufficient to establish and maintain these domains during organogenesis in several multicellular organisms. Here, we show that the mutual inhibition between the PLETHORA transcription factors (PLTs) and the ARABIDOPSIS RESPONSE REGULATORs (ARRs) transcription factors is sufficient to separate cell division and cell differentiation during root organogenesis. In particular, we show that ARR1 suppresses PLT activities and that PLTs suppress ARR1 and ARR12 by targeting their proteins for degradation via the KISS ME DEADLY 2 F-box protein. These findings reveal new important aspects of the complex process of root zonation and development.

Synthetic spatial patterning in bacteria: advances based on novel diffusible signals

Engineering multicellular patterning may help in the understanding of some fundamental laws of pattern formation and thus may contribute to the field of developmental biology. Furthermore, advanced spatial control over gene expression may revolutionize fields such as medicine, through organoid or tissue engineering. To date, foundational advances in spatial synthetic biology have often been made in prokaryotes, using artificial gene circuits. In this review, engineered patterns are classified into four levels of increasing complexity, ranging from spatial systems with no diffusible signals to systems with complex multi-diffusor interactions. This classification highlights how the field was held back by a lack of diffusible components. Consequently, we provide a summary of both previously characterized and some new potential candidate small-molecule signals that can regulate gene expression in Escherichia coli. These diffusive signals will help synthetic biologists to successfully engineer increasingly intricate, robust and tuneable spatial structures.

Relationship between epithelial organization and morphogen interpretation

Despite molecular noise and genetic differences between individuals, developmental outcomes are remarkably constant. Decades of research has focused on the underlying mechanisms that ensure this precision and robustness. Recent quantifications of chemical gradients and epithelial cell shapes provide novel insights into the basis of precise development. In this review, we argue that these two aspects may be linked in epithelial morphogenesis.

Synthetic cell-based materials extract positional information from morphogen gradients

Biomaterials composed of synthetic cells have the potential to adapt and differentiate guided by physicochemical environmental cues. Inspired by biological systems in development, which extract positional information (PI) from morphogen gradients in the presence of uncertainties, we here investigate how well synthetic cells can determine their position within a multicellular structure. To calculate PI, we created and analyzed a large number of synthetic cellular assemblies composed of emulsion droplets connected via lipid bilayer membranes. These droplets contained cell-free feedback gene circuits that responded to gradients of a genetic inducer acting as a morphogen. PI is found to be limited by gene expression noise and affected by the temporal evolution of the morphogen gradient and the cell-free expression system itself. The generation of PI can be rationalized by computational modeling of the system. We scale our approach using three-dimensional printing and demonstrate morphogen-based differentiation in larger tissue-like assemblies.

Modular, robust and extendible multicellular circuit design in yeast

Division of labor between cells is ubiquitous in biology but the use of multi-cellular consortia for engineering applications is only beginning to be explored. A significant advantage of multi-cellular circuits is their potential to be modular with respect to composition but this claim has not yet been extensively tested using experiments and quantitative modeling. Here, we construct a library of 24 yeast strains capable of sending, receiving or responding to three molecular signals, characterize them experimentally and build quantitative models of their input-output relationships. We then compose these strains into two- and three-strain cascades as well as a four-strain bistable switch and show that experimentally measured consortia dynamics can be predicted from the models of the constituent parts. To further explore the achievable range of behaviors, we perform a fully automated computational search over all two-, three- and four-strain consortia to identify combinations that realize target behaviors including logic gates, band-pass filters and time pulses. Strain combinations that are predicted to map onto a target behavior are further computationally optimized and then experimentally tested. Experiments closely track computational predictions. The high reliability of these model descriptions further strengthens the feasibility and highlights the potential for distributed computing in synthetic biology.

Implications of diffusion and time-varying morphogen gradients for the dynamic positioning and precision of bistable gene expression boundaries

The earliest models for how morphogen gradients guide embryonic patterning failed to account for experimental observations of temporal refinement in gene expression domains. Following theoretical and experimental work in this area, dynamic positional information has emerged as a conceptual framework to discuss how cells process spatiotemporal inputs into downstream patterns. Here, we show that diffusion determines the mathematical means by which bistable gene expression boundaries shift over time, and therefore how cells interpret positional information conferred from morphogen concentration. First, we introduce a metric for assessing reproducibility in boundary placement or precision in systems where gene products do not diffuse, but where morphogen concentrations are permitted to change in time. We show that the dynamics of the gradient affect the sensitivity of the final pattern to variation in initial conditions, with slower gradients reducing the sensitivity. Second, we allow gene products to diffuse and consider gene expression boundaries as propagating wavefronts with velocity modulated by local morphogen concentration. We harness this perspective to approximate a PDE model as an ODE that captures the position of the boundary in time, and demonstrate the approach with a preexisting model for Hunchback patterning in fruit fly embryos. We then propose a design that employs antiparallel morphogen gradients to achieve accurate boundary placement that is robust to scaling. Throughout our work we draw attention to tradeoffs among initial conditions, boundary positioning, and the relative timescales of network and gradient evolution. We conclude by suggesting that mathematical theory should serve to clarify not just our quantitative, but also our intuitive understanding of patterning processes.