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.

Gene expression model inference from snapshot RNA data using Bayesian non-parametrics

Accelerated Bayesian inference of gene expression models from snapshots of single-cell transcripts

Non-parametric Bayesian modelling of digital gene expression data

Inferring gene regulatory networks from single-cell data: a mechanistic approach

A deep generative model for gene expression profiles from single-cell RNA sequencing

A deep generative model for single-cell RNA sequencing with application to detecting differentially expressed genes

Bayesian inference of biochemical kinetic parameters using the linear noise approximation

Reverse engineering gene regulatory networks using approximate Bayesian computation

Approximate inference of gene regulatory network models from RNA-Seq time series data

A Bayesian framework for the inference of gene regulatory networks from time and pseudo‐time series data

A single-cell gene expression language model

Modeling dependent gene expression

Bayesian Inference of RNA Velocity from Multi-Lineage Single-Cell Data

Integrated Pipelines for Inferring Gene Regulatory Networks from Single-Cell Data

Reconstructing Gene Regulatory Networks with Bayesian Networks by Combining Expression Data with Multiple Sources of Prior Knowledge

Gene network inference from incomplete expression data: transcriptional control of hematopoietic commitment

Deciphering regulatory architectures from synthetic single-cell expression patterns

A Bayesian graphical modeling approach to microRNA regulatory network inference

An empirical Bayes approach to estimating dynamic models of co-regulated gene expression

Estimation of Isoform Expression in Rna-Seq Data Using A Hierarchical Bayesian Model

A Unified Probabilistic Modeling Framework for Eukaryotic Transcription Based on Nascent RNA Sequencing Data