Abstract:For the vast majority of genes in sequenced genomes, there is limited understanding of how they are regulated. Without such knowledge, it is not possible to perform a quantitative theory-experiment dialogue on how such genes give rise to physiological and evolutionary adaptation. One category of high-throughput experiments used to understand the sequence-phenotype relationship of the transcriptome is massively parallel reporter assays (MPRAs). However, to improve the versatility and scalability of MPRA pipelines, we need a "theory of the experiment" to help us better understand the impact of various biological and experimental parameters on the interpretation of experimental data. These parameters include binding site copy number, where a large number of specific binding sites may titrate away transcription factors, as well as the presence of overlapping binding sites, which may affect analysis of the degree of mutual dependence between mutations in the regulatory region and expression levels. To that end, in this paper we create tens of thousands of synthetic single-cell gene expression outputs using both equilibrium and out-of-equilibrium models. These models make it possible to imitate the summary statistics (information footprints and expression shift matrices) used to characterize the output of MPRAs and from this summary statistic to infer the underlying regulatory architecture. Specifically, we use a more refined implementation of the so-called thermodynamic models in which the binding energies of each sequence variant are derived from energy matrices. Our simulations reveal important effects of the parameters on MPRA data and we demonstrate our ability to optimize MPRA experimental designs with the goal of generating thermodynamic models of the transcriptome with base-pair specificity. Further, this approach makes it possible to carefully examine the mapping between mutations in binding sites and their corresponding expression profiles, a tool useful not only for better designing MPRAs, but also for exploring regulatory evolution.

Genome-scale High-Resolution Mapping of Activating and Repressive Nucleotides in Regulatory Regions

Single-cell Multi-Scale Footprinting Reveals the Modular Organization of DNA Regulatory Elements

Systematic Dissection and Optimization of Inducible Enhancers in Human Cells Using a Massively Parallel Reporter Assay

Chromatin Potential Identified by Shared Single Cell Profiling of RNA and Chromatin

Systematic dissection of regulatory motifs in 2000 predicted human enhancers using a massively parallel reporter assay.

High-resolution genome-wide functional dissection of transcriptional regulatory regions and nucleotides in human

MPRAnalyze: Statistical Framework for Massively Parallel Reporter Assays

ATAC-STARR-seq reveals transcription factor-bound activators and silencers within chromatin-accessible regions of the human genome

Massively parallel reporter assay: a novel technique for analyzing the regulation of gene expression.

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

Characterization of enhancer activity in early human neurodevelopment using Massively Parallel Reporter Assay (MPRA) and forebrain organoids

A genome-wide nucleosome-resolution map of promoter-centered interactions in human cells corroborates the enhancer-promoter looping model

Combinatorial Patterning of Chromatin Regulators Uncovered by Genome-wide Location Analysis in Human Cells

Deciphering regulatory architectures from synthetic single-cell expression patterns

Mapping and Analysis of Chromatin State Dynamics in Nine Human Cell Types

Systematic Dissection of Genomic Features Determining Transcription Factor Binding and Enhancer Function

A Tiling-Deletion-Based Genetic Screen for Cis-Regulatory Element Identification in Mammalian Cells

Sensitive dissection of a genomic regulatory landscape using bulk and targeted single-cell activation

Massively parallel reporter assays and mouse transgenic assays provide complementary information about neuronal enhancer activity

Coupling high-throughput mapping with proteomics analysis delineates cis -regulatory elements at high resolution

Decoding biology with massively parallel reporter assays and machine learning