Synthetic dosage-compensating miRNA circuits allow precision gene therapy for Rett syndrome.

A longstanding challenge in gene therapy is expressing a dosage-sensitive gene within a tight therapeutic window. For example, loss of MECP2 function causes Rett syndrome, while its duplication causes MECP2 duplication syndrome. Viral gene delivery methods generate variable numbers of gene copies in individual cells, creating a need for gene dosage-invariant expression systems. Here, we introduce a compact miRNA-based, incoherent feed-forward loop circuit that achieves precise control of Mecp2 expression in cells and brains, and improves outcomes in an AAV-based mouse model of Rett syndrome gene therapy. Single molecule analysis of endogenous and ectopic Mecp2 mRNA revealed precise, sustained expression across a broad range of gene dosages. Delivered systemically in a brain-targeting AAV capsid, the circuit strongly suppressed Rett behavioral symptoms for over 24 weeks, outperforming an unregulated gene therapy. These results demonstrate that synthetic miRNA-based regulatory circuits can enable precise in vivo expression to improve the safety and efficacy of gene therapy.

miRNA circuit modules for precise, tunable control of gene expression.

The ability to express transgenes at specified levels is critical for understanding cellular behaviors, and for applications in gene and cell therapy. Transfection, viral vectors, and other gene delivery methods produce varying protein expression levels, with limited quantitative control, while targeted knock-in and stable selection are inefficient and slow. Active compensation mechanisms can improve precision, but the need for additional proteins or lack of tunability have prevented their widespread use. Here, we introduce a toolkit of compact, synthetic miRNA-based circuit modules that provide precise, tunable control of transgenes across diverse cell types. These circuits, termed DIMMERs (Dosage-Invariant miRNA-Mediated Expression Regulators) use multivalent miRNA regulatory interactions within an incoherent feed-forward loop architecture to achieve nearly uniform protein expression over more than two orders of magnitude variation in underlying gene dosages or transcription rates. They also allow coarse and fine control of expression, and are portable, functioning across diverse cell types. In addition, a heuristic miRNA design algorithm enables the creation of orthogonal circuit variants that independently control multiple genes in the same cell. These circuits allowed dramatically improved CRISPR imaging, and super-resolution imaging of EGFR receptors with transient transfections. The toolbox provided here should allow precise, tunable, dosage-invariant expression for research, gene therapy, and other biotechnology applications.

Synthetic protein circuits for programmable control of mammalian cell death.

Natural cell death pathways such as apoptosis and pyroptosis play dual roles: they eliminate harmful cells and modulate the immune system by dampening or stimulating inflammation. Synthetic protein circuits capable of triggering specific death programs in target cells could similarly remove harmful cells while appropriately modulating immune responses. However, cells actively influence their death modes in response to natural signals, making it challenging to control death modes. Here, we introduce naturally inspired "synpoptosis" circuits that proteolytically regulate engineered executioner proteins and mammalian cell death. These circuits direct cell death modes, respond to combinations of protease inputs, and selectively eliminate target cells. Furthermore, synpoptosis circuits can be transmitted intercellularly, offering a foundation for engineering synthetic killer cells that induce desired death programs in target cells without self-destruction. Together, these results lay the groundwork for programmable control of mammalian cell death.

Lineage motifs as developmental modules for control of cell type proportions.

In multicellular organisms, cell types must be produced and maintained in appropriate proportions. One way this is achieved is through committed progenitor cells or extrinsic interactions that produce specific patterns of descendant cell types on lineage trees. However, cell fate commitment is probabilistic in most contexts, making it difficult to infer these dynamics and understand how they establish overall cell type proportions. Here, we introduce Lineage Motif Analysis (LMA), a method that recursively identifies statistically overrepresented patterns of cell fates on lineage trees as potential signatures of committed progenitor states or extrinsic interactions. Applying LMA to published datasets reveals spatial and temporal organization of cell fate commitment in zebrafish and rat retina and early mouse embryonic development. Comparative analysis of vertebrate species suggests that lineage motifs facilitate adaptive evolutionary variation of retinal cell type proportions. LMA thus provides insight into complex developmental processes by decomposing them into simpler underlying modules.

Combinatorial expression motifs in signaling pathways.

In animal cells, molecular pathways often comprise families of variant components, such as ligands or receptors. These pathway components are differentially expressed by different cell types, potentially tailoring pathway function to cell context. However, it has remained unclear how pathway expression profiles are distributed across cell types and whether similar profiles can occur in dissimilar cell types. Here, using single-cell gene expression datasets, we identified pathway expression motifs, defined as recurrent expression profiles that are broadly distributed across diverse cell types. Motifs appeared in core pathways, including TGF-ß, Notch, Wnt, and the SRSF splice factors, and involved combinatorial co-expression of multiple components. Motif usage was weakly correlated between pathways in adult cell types and during dynamic developmental transitions. Together, these results suggest a mosaic view of cell type organization, in which different cell types operate many of the same pathways in distinct modes.

Principles of Computation by Competitive Protein Dimerization Networks.

Many biological signaling pathways employ proteins that competitively dimerize in diverse combinations. These dimerization networks can perform biochemical computations, in which the concentrations of monomers (inputs) determine the concentrations of dimers (outputs). Despite their prevalence, little is known about the range of input-output computations that dimerization networks can perform (their "expressivity") and how it depends on network size and connectivity. Using a systematic computational approach, we demonstrate that even small dimerization networks (3-6 monomers) are expressive, performing diverse multi-input computations. Further, dimerization networks are versatile, performing different computations when their protein components are expressed at different levels, such as in different cell types. Remarkably, individual networks with random interaction affinities, when large enough (≥8 proteins), can perform nearly all (~90%) potential one-input network computations merely by tuning their monomer expression levels. Thus, even the simple process of competitive dimerization provides a powerful architecture for multi-input, cell-type-specific signal processing.

Diversity in Notch ligand-receptor signaling interactions.

The Notch signaling pathway uses families of ligands and receptors to transmit signals to nearby cells. These components are expressed in diverse combinations in different cell types, interact in a many-to-many fashion, both within the same cell (in cis) and between cells (in trans), and their interactions are modulated by Fringe glycosyltransferases. A fundamental question is how the strength of Notch signaling depends on which pathway components are expressed, at what levels, and in which cells. Here, we used a quantitative, bottom-up, cell-based approach to systematically characterize trans-activation, cis-inhibition, and cis-activation signaling efficiencies across a range of ligand and Fringe expression levels in two mammalian cell types. Each ligand (Dll1, Dll4, Jag1, and Jag2) and receptor variant (Notch1 and Notch2) analyzed here exhibited a unique profile of interactions, Fringe-dependence, and signaling outcomes. All four ligands were able to bind receptors in cis and in trans, and all ligands trans-activated both receptors except for Jag1, which failed to activate Notch1. Cis-interactions were predominantly inhibitory, with the exception of the Dll1- and Dll4-Notch2 pairs, which exhibited cis-activation stronger than trans-activation. Lfng strengthened Delta-mediated trans-activation and weakened Jagged-mediated trans-activation for both receptors. Finally, cis-ligands showed diverse cis-inhibition strengths, which depended on the identity of the trans-ligand as well as the receptor. The map of receptor-ligand-Fringe interaction outcomes revealed here should help guide rational perturbation and control of the Notch pathway.

Tuning Methylation-Dependent Silencing Dynamics by Synthetic Modulation of CpG Density.

Methylation of cytosines in CG dinucleotides (CpGs) within promoters has been shown to lead to gene silencing in mammals in natural contexts. Recently, engineered recruitment of methyltransferases (DNMTs) at specific loci was shown to be sufficient to silence synthetic and endogenous gene expression through this mechanism. A critical parameter for DNA methylation-based silencing is the distribution of CpGs within the target promoter. However, how the number or density of CpGs in the target promoter affects the dynamics of silencing by DNMT recruitment has remained unclear. Here, we constructed a library of promoters with systematically varying CpG content, and analyzed the rate of silencing in response to recruitment of DNMT. We observed a tight correlation between silencing rate and CpG content. Further, methylation-specific analysis revealed a constant accumulation rate of methylation at the promoter after DNMT recruitment. We identified a single CpG site between TATA box and transcription start site (TSS) that accounted for a substantial part of the difference in silencing rates between promoters with differing CpG content, indicating that certain residues play disproportionate roles in controlling silencing. Together, these results provide a library of promoters for synthetic epigenetic and gene regulation applications, as well as insights into the regulatory link between CpG content and silencing rate.

Tuning methylation-dependent silencing dynamics by synthetic modulation of CpG density.

Methylation of cytosines in CG dinucleotides (CpGs) within promoters has been shown to lead to gene silencing in mammals in natural contexts. Recently, engineered recruitment of methyltransferases (DNMTs) at specific loci was shown to be sufficient to silence synthetic and endogenous gene expression through this mechanism. A critical parameter for DNA methylation-based silencing is the distribution of CpGs within the target promoter. However, how the number or density of CpGs in the target promoter affects the dynamics of silencing by DNMT recruitment has remained unclear. Here we constructed a library of promoters with systematically varying CpG content, and analyzed the rate of silencing in response to recruitment of DNMT. We observed a tight correlation between silencing rate and CpG content. Further, methylation-specific analysis revealed a constant accumulation rate of methylation at the promoter after DNMT recruitment. We identified a single CpG site between TATA box and transcription start site (TSS) that accounted for a substantial part of the difference in silencing rates between promoters with differing CpG content, indicating that certain residues play disproportionate roles in controlling silencing. Together, these results provide a library of promoters for synthetic epigenetic and gene regulation applications, as well as insights into the regulatory link between CpG content and silencing rate.

Reconstitution of morphogen shuttling circuits.

Developing tissues form spatial patterns by establishing concentration gradients of diffusible signaling proteins called morphogens. The bone morphogenetic protein (BMP) morphogen pathway uses a family of extracellular modulators to reshape signaling gradients by actively "shuttling" ligands to different locations. It has remained unclear what circuits are sufficient to enable shuttling, what other patterns they can generate, and whether shuttling is evolutionarily conserved. Here, using a synthetic, bottom-up approach, we compared the spatiotemporal dynamics of different extracellular circuits. Three proteins-Chordin, Twsg, and the BMP-1 protease-successfully displaced gradients by shuttling ligands away from the site of production. A mathematical model explained the different spatial dynamics of this and other circuits. Last, combining mammalian and Drosophila components in the same system suggests that shuttling is a conserved capability. Together, these results reveal principles through which extracellular circuits control the spatiotemporal dynamics of morphogen signaling.

Engineering RNA export for measurement and manipulation of living cells.

A system for programmable export of RNA molecules from living cells would enable both non-destructive monitoring of cell dynamics and engineering of cells capable of delivering executable RNA programs to other cells. We developed genetically encoded cellular RNA exporters, inspired by viruses, that efficiently package and secrete cargo RNA molecules from mammalian cells within protective nanoparticles. Exporting and sequencing RNA barcodes enabled non-destructive monitoring of cell population dynamics with clonal resolution. Further, by incorporating fusogens into the nanoparticles, we demonstrated the delivery, expression, and functional activity of exported mRNA in recipient cells. We term these systems COURIER (controlled output and uptake of RNA for interrogation, expression, and regulation). COURIER enables measurement of cell dynamics and establishes a foundation for hybrid cell and gene therapies based on cell-to-cell delivery of RNA.

Lineage motifs: developmental modules for control of cell type proportions.

In multicellular organisms, cell types must be produced and maintained in appropriate proportions. One way this is achieved is through committed progenitor cells that produce specific sets of descendant cell types. However, cell fate commitment is probabilistic in most contexts, making it difficult to infer progenitor states and understand how they establish overall cell type proportions. Here, we introduce Lineage Motif Analysis (LMA), a method that recursively identifies statistically overrepresented patterns of cell fates on lineage trees as potential signatures of committed progenitor states. Applying LMA to published datasets reveals spatial and temporal organization of cell fate commitment in zebrafish and rat retina and early mouse embryo development. Comparative analysis of vertebrate species suggests that lineage motifs facilitate adaptive evolutionary variation of retinal cell type proportions. LMA thus provides insight into complex developmental processes by decomposing them into simpler underlying modules.

The computational capabilities of many-to-many protein interaction networks.

Many biological circuits comprise sets of protein variants that interact with one another in a many-to-many, or promiscuous, fashion. These architectures can provide powerful computational capabilities that are especially critical in multicellular organisms. Understanding the principles of biochemical computations in these circuits could allow more precise control of cellular behaviors. However, these systems are inherently difficult to analyze, due to their large number of interacting molecular components, partial redundancies, and cell context dependence. Here, we discuss recent experimental and theoretical advances that are beginning to reveal how promiscuous circuits compute, what roles those computations play in natural biological contexts, and how promiscuous architectures can be applied for the design of synthetic multicellular behaviors.

Alternative splicing of GSDMB modulates killer lymphocyte-triggered pyroptosis.

Granzyme A from killer lymphocytes cleaves gasdermin B (GSDMB) and triggers pyroptosis in targeted human tumor cells, eliciting antitumor immunity. However, GSDMB has a controversial role in pyroptosis and has been linked to both anti- and protumor functions. Here, we found that GSDMB splicing variants are functionally distinct. Cleaved N-terminal (NT) fragments of GSDMB isoforms 3 and 4 caused pyroptosis, but isoforms 1, 2, and 5 did not. The nonfunctional isoforms have a deleted or modified exon 6 and therefore lack a stable belt motif. The belt likely contributes to the insertion of oligomeric GSDMB-NTs into the membrane. Consistently, noncytotoxic GSDMB-NTs blocked pyroptosis caused by cytotoxic GSDMB-NTs in a dominant-negative manner. Upon natural killer (NK) cell attack, GSDMB3-expressing cells died by pyroptosis, whereas GSDMB4-expressing cells died by mixed pyroptosis and apoptosis, and GSDMB1/2-expressing cells died only by apoptosis. GSDMB4 partially resisted NK cell-triggered cleavage, suggesting that only GSDMB3 is fully functional. GSDMB1-3 were the most abundant isoforms in the tested tumor cell lines and were similarly induced by interferon-γ and the chemotherapy drug methotrexate. Expression of cytotoxic GSDMB3/4 isoforms, but not GSDMB1/2 isoforms that are frequently up-regulated in tumors, was associated with better outcomes in bladder and cervical cancers, suggesting that GSDMB3/4-mediated pyroptosis was protective in those tumors. Our study indicates that tumors may block and evade killer cell-triggered pyroptosis by generating noncytotoxic GSDMB isoforms. Therefore, therapeutics that favor the production of cytotoxic GSDMB isoforms by alternative splicing may improve antitumor immunity.

The sound of silence: Transgene silencing in mammalian cell engineering.

To elucidate principles operating in native biological systems and to develop novel biotechnologies, synthetic biology aims to build and integrate synthetic gene circuits within native transcriptional networks. The utility of synthetic gene circuits for cell engineering relies on the ability to control the expression of all constituent transgene components. Transgene silencing, defined as the loss of expression over time, persists as an obstacle for engineering primary cells and stem cells with transgenic cargos. In this review, we highlight the challenge that transgene silencing poses to the robust engineering of mammalian cells, outline potential molecular mechanisms of silencing, and present approaches for preventing transgene silencing. We conclude with a perspective identifying future research directions for improving the performance of synthetic gene circuits.

Periodic spatial patterning with a single morphogen.

During multicellular development, periodic spatial patterning systems generate repetitive structures, such as digits, vertebrae, and teeth. Turing patterning provides a foundational paradigm for understanding such systems. The simplest Turing systems are believed to require at least two morphogens to generate periodic patterns. Here, using mathematical modeling, we show that a simpler circuit, including only a single diffusible morphogen, is sufficient to generate long-range, spatially periodic patterns that propagate outward from transient initiating perturbations and remain stable after the perturbation is removed. Furthermore, an additional bistable intracellular feedback or operation on a growing cell lattice can make patterning robust to noise. Together, these results show that a single morphogen can be sufficient for robust spatial pattern formation and should provide a foundation for engineering pattern formation in the emerging field of synthetic developmental biology.

Embryo model completes gastrulation to neurulation and organogenesis.

Embryonic stem (ES) cells can undergo many aspects of mammalian embryogenesis in vitro1-5, but their developmental potential is substantially extended by interactions with extraembryonic stem cells, including trophoblast stem (TS) cells, extraembryonic endoderm stem (XEN) cells and inducible XEN (iXEN) cells6-11. Here we assembled stem cell-derived embryos in vitro from mouse ES cells, TS cells and iXEN cells and showed that they recapitulate the development of whole natural mouse embryo in utero up to day 8.5 post-fertilization. Our embryo model displays headfolds with defined forebrain and midbrain regions and develops a beating heart-like structure, a trunk comprising a neural tube and somites, a tail bud containing neuromesodermal progenitors, a gut tube, and primordial germ cells. This complete embryo model develops within an extraembryonic yolk sac that initiates blood island development. Notably, we demonstrate that the neurulating embryo model assembled from Pax6-knockout ES cells aggregated with wild-type TS cells and iXEN cells recapitulates the ventral domain expansion of the neural tube that occurs in natural, ubiquitous Pax6-knockout embryos. Thus, these complete embryoids are a powerful in vitro model for dissecting the roles of diverse cell lineages and genes in development. Our results demonstrate the self-organization ability of ES cells and two types of extraembryonic stem cells to reconstitute mammalian development through and beyond gastrulation to neurulation and early organogenesis.

Dynamics and functional roles of splicing factor autoregulation.

Non-core spliceosome components are essential, conserved regulators of alternative splicing. They provide concentration-dependent control of diverse pre-mRNAs. Many splicing factors direct unproductive splicing of their own pre-mRNAs through negative autoregulation. However, the impact of such feedback loops on splicing dynamics at the single-cell level remains unclear. Here, we developed a system to quantitatively analyze negative autoregulatory splicing dynamics by splicing factor SRSF1 in response to perturbations in single HEK293 cells. We show that negative autoregulatory splicing provides critical functions for gene regulation, establishing a ceiling of SRSF1 protein concentration, reducing cell-cell heterogeneity in SRSF1 levels, and buffering variation in transcription. Most important, it adapts SRSF1 splicing activity to variations in demand from other pre-mRNA substrates. A minimal mathematical model of autoregulatory splicing explains these experimentally observed features and provides values for effective biochemical parameters. These results reveal the unique functional roles that splicing negative autoregulation plays in homeostatically regulating transcriptional programs.

The context-dependent, combinatorial logic of BMP signaling.

Cell-cell communication systems typically comprise families of ligand and receptor variants that function together in combinations. Pathway activation depends on the complex way in which ligands are presented extracellularly and receptors are expressed by the signal-receiving cell. To understand the combinatorial logic of such a system, we systematically measured pairwise bone morphogenetic protein (BMP) ligand interactions in cells with varying receptor expression. Ligands could be classified into equivalence groups based on their profile of positive and negative synergies with other ligands. These groups varied with receptor expression, explaining how ligands can functionally replace each other in one context but not another. Context-dependent combinatorial interactions could be explained by a biochemical model based on the competitive formation of alternative signaling complexes with distinct activities. Together, these results provide insights into the roles of BMP combinations in developmental and therapeutic contexts and establish a framework for analyzing other combinatorial, context-dependent signaling systems.

Ligand-receptor promiscuity enables cellular addressing.

In multicellular organisms, secreted ligands selectively activate, or "address," specific target cell populations to control cell fate decision-making and other processes. Key cell-cell communication pathways use multiple promiscuously interacting ligands and receptors, provoking the question of how addressing specificity can emerge from molecular promiscuity. To investigate this issue, we developed a general mathematical modeling framework based on the bone morphogenetic protein (BMP) pathway architecture. We find that promiscuously interacting ligand-receptor systems allow a small number of ligands, acting in combinations, to address a larger number of individual cell types, defined by their receptor expression profiles. Promiscuous systems outperform seemingly more specific one-to-one signaling architectures in addressing capability. Combinatorial addressing extends to groups of cell types, is robust to receptor expression noise, grows more powerful with increases in the number of receptor variants, and is maximized by specific biochemical parameter relationships. Together, these results identify design principles governing cellular addressing by ligand combinations.