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.
Full author list & citation
Jean-Benoît Lalanne, Samuel G. Regalado, Silvia Domcke, Diego Calderon, Beth K. Martin, Xiaoyi Li, Tony Li, Chase C. Suiter, Choli Lee, Cole Trapnell, Jay Shendure. Multiplex profiling of developmental cis-regulatory elements with quantitative single-cell expression reporters. 2024-05-09. https://doi.org/10.1038/s41592-024-02260-3
A piggyBac-integrated library of developmental CREs and promoter controls was measured by bulk MPRA during mouse embryoid-body differentiation. The table aggregates mBC-level GSE217679 counts per CRE, replicate and time point and includes the corresponding published activity summary when it matched.
Five promoter constructs (noP, minP, UBCp, Pgk1p and EEF1A1p) were tested with barcode-linked reporters in HEK293T, HepG2 and K562. The table retains matched author-provided bulk RNA/DNA activity and single-cell reporter activity at the mBC/oBC level.
The study used scQers in mouse embryoid bodies to test paired CREs in alternative orientations, TFBS-perturbed variants of developmental CREs, and literature-selected CRE controls. The table combines the author-provided cluster-level activity summaries into one long-form table.
Twenty promoter libraries spanning five promoters and four cHS4/U6-oBC reporter architectures were bottlenecked to estimated 3k, 10k or 30k starting clones in K562 and HEK293. The table reports the source per-barcode activity distribution and its interquartile fold range.
The study screened 204 accessible-chromatin developmental CREs from 23 loci, together with promoter controls, in mouse embryoid bodies. This table contains aggregate scQer activity and specificity with tested mm10 coordinates and cognate ATAC annotations where available.
Circular U6-Tornado barcodes and truncated linear barcodes were integrated into K562 cells with piggyBac and compared by bulk RNA/DNA barcode counting. The table provides per-barcode normalized activity for the two barcode architectures in duplicate biological replicates.
Pooled piggyBac reporter libraries varied U6/oBC presence, CRE position, promoter configuration and reporter architecture during mouse embryoid-body differentiation. The table contains construct-level winsorized RNA/DNA activity across days 0, 4, 18, 20 and 22.
Original supplemental and deposited inputs retained for this study. Download files individually or together as a ZIP; nested folders are preserved. Source reuse terms apply, and sequencing reads may be omitted.