Symbolic recording of signalling and cis-regulatory element activity to DNA.

Measurements of gene expression or signal transduction activity are conventionally performed using methods that require either the destruction or live imaging of a biological sample within the timeframe of interest. Here we demonstrate an alternative paradigm in which such biological activities are stably recorded to the genome. Enhancer-driven genomic recording of transcriptional activity in multiplex (ENGRAM) is based on the signal-dependent production of prime editing guide RNAs that mediate the insertion of signal-specific barcodes (symbols) into a genomically encoded recording unit. We show how this strategy can be used for multiplex recording of the cell-type-specific activities of dozens to hundreds of cis-regulatory elements with high fidelity, sensitivity and reproducibility. Leveraging signal transduction pathway-responsive cis-regulatory elements, we also demonstrate time- and concentration-dependent genomic recording of WNT, NF-κB and Tet-On activities. By coupling ENGRAM to sequential genome editing via DNA Typewriter1, we stably record information about the temporal dynamics of two orthogonal signalling pathways to genomic DNA. Finally we apply ENGRAM to integratively record the transient activity of nearly 100 transcription factor consensus motifs across daily windows spanning the differentiation of mouse embryonic stem cells into gastruloids, an in vitro model of early mammalian development. Although these are proof-of-concept experiments and much work remains to fully realize the possibilities, the symbolic recording of biological signals or states within cells, to the genome and over time, has broad potential to complement contemporary paradigms for how we make measurements in biological systems.

Multiplex profiling of developmental cis-regulatory elements with quantitative single-cell expression reporters.

The inability to scalably and precisely measure the activity of developmental cis-regulatory elements (CREs) in multicellular systems is a bottleneck in genomics. Here we develop a dual RNA cassette that decouples the detection and quantification tasks inherent to multiplex single-cell reporter assays. The resulting measurement of reporter expression is accurate over multiple orders of magnitude, with a precision approaching the limit set by Poisson counting noise. Together with RNA barcode stabilization via circularization, these scalable single-cell quantitative expression reporters provide high-contrast readouts, analogous to classic in situ assays but entirely from sequencing. Screening >200 regions of accessible chromatin in a multicellular in vitro model of early mammalian development, we identify 13 (8 previously uncharacterized) autonomous and cell-type-specific developmental CREs. We further demonstrate that chimeric CRE pairs generate cognate two-cell-type activity profiles and assess gain- and loss-of-function multicellular expression phenotypes from CRE variants with perturbed transcription factor binding sites. Single-cell quantitative expression reporters can be applied in developmental and multicellular systems to quantitatively characterize native, perturbed and synthetic CREs at scale, with high sensitivity and at single-cell resolution.

Multiplex single-cell chemical genomics reveals the kinase dependence of the response to targeted therapy.

Chemical genetic screens are a powerful tool for exploring how cancer cells' response to drugs is shaped by their mutations, yet they lack a molecular view of the contribution of individual genes to the response to exposure. Here, we present sci-Plex-Gene-by-Environment (sci-Plex-GxE), a platform for combined single-cell genetic and chemical screening at scale. We highlight the advantages of large-scale, unbiased screening by defining the contribution of each of 522 human kinases to the response of glioblastoma to different drugs designed to abrogate signaling from the receptor tyrosine kinase pathway. In total, we probed 14,121 gene-by-environment combinations across 1,052,205 single-cell transcriptomes. We identify an expression signature characteristic of compensatory adaptive signaling regulated in a MEK/MAPK-dependent manner. Further analyses aimed at preventing adaptation revealed promising combination therapies, including dual MEK and CDC7/CDK9 or nuclear factor κB (NF-κB) inhibitors, as potent means of preventing transcriptional adaptation of glioblastoma to targeted therapy.

Induction and in silico staging of human gastruloids with neural tube, segmented somites & advanced cell types.

Embryonic organoids are emerging as powerful models for studying early mammalian development. For example, stem cell-derived 'gastruloids' form elongating structures containing all three germ layers1-4. However, although elongated, human gastruloids do not morphologically resemble post-implantation embryos. Here we show that a specific, discontinuous regimen of retinoic acid (RA) robustly induces human gastruloids with embryo-like morphological structures, including a neural tube and segmented somites. Single cell RNA-seq (sc-RNA-seq) further reveals that these human 'RA-gastruloids' contain more advanced cell types than conventional gastruloids, including neural crest cells, renal progenitor cells, skeletal muscle cells, and, rarely, neural progenitor cells. We apply a new approach to computationally stage human RA-gastruloids relative to somite-resolved mouse embryos, early human embryos and other gastruloid models, and find that the developmental stage of human RA-gastruloids is comparable to that of E9.5 mouse embryos, although some cell types show greater or lesser progression. We chemically perturb WNT and BMP signaling in human RA-gastruloids and find that these signaling pathways regulate somite patterning and neural tube length, respectively, while genetic perturbation of the transcription factors PAX3 and TBX6 markedly compromises the formation of neural crest and somites/renal cells, respectively. Human RA-gastruloids complement other embryonic organoids in serving as a simple, robust and screenable model for decoding early human embryogenesis.

Retinoic acid induces human gastruloids with posterior embryo-like structures.

Gastruloids are a powerful in vitro model of early human development. However, although elongated and composed of all three germ layers, human gastruloids do not morphologically resemble post-implantation human embryos. Here we show that an early pulse of retinoic acid (RA), together with later Matrigel, robustly induces human gastruloids with posterior embryo-like morphological structures, including a neural tube flanked by segmented somites and diverse cell types, including neural crest, neural progenitors, renal progenitors and myocytes. Through in silico staging based on single-cell RNA sequencing, we find that human RA-gastruloids progress further than other human or mouse embryo models, aligning to E9.5 mouse and CS11 cynomolgus monkey embryos. We leverage chemical and genetic perturbations of RA-gastruloids to confirm that WNT and BMP signalling regulate somite formation and neural tube length in the human context, while transcription factors TBX6 and PAX3 underpin presomitic mesoderm and neural crest, respectively. Looking forward, RA-gastruloids are a robust, scalable model for decoding early human embryogenesis.