Abstract:BackgroundThe recent development of single-cell transcriptomics has enabled gene expression to be measured in individual cells instead of being population-averaged. Despite this considerable precision improvement, inferring regulatory networks remains challenging because stochasticity now proves to play a fundamental role in gene expression. In particular, mRNA synthesis is now acknowledged to occur in a highly bursty manner.ResultsWe propose to view the inference problem as a fitting procedure for a mechanistic gene network model that is inherently stochastic and takes not only protein, but also mRNA levels into account. We first explain how to build and simulate this network model based upon the coupling of genes that are described as piecewise-deterministic Markov processes. Our model is modular and can be used to implement various biochemical hypotheses including causal interactions between genes. However, a naive fitting procedure would be intractable. By performing a relevant approximation of the stationary distribution, we derive a tractable procedure that corresponds to a statistical hidden Markov model with interpretable parameters. This approximation turns out to be extremely close to the theoretical distribution in the case of a simple toggle-switch, and we show that it can indeed fit real single-cell data. As a first step toward inference, our approach was applied to a number of simple two-gene networks simulated in silico from the mechanistic model and satisfactorily recovered the original networks.ConclusionsOur results demonstrate that functional interactions between genes can be inferred from the distribution of a mechanistic, dynamical stochastic model that is able to describe gene expression in individual cells. This approach seems promising in relation to the current explosion of single-cell expression data.

Modelling the burden caused by gene expression: an in silico investigation into the interactions between synthetic gene circuits and their chassis cell

Prediction of cellular burden with host-circuit models

Thermodynamically-consistent, reduced models of gene regulatory networks

A coarse-grained bacterial cell model for resource-aware analysis and design of synthetic gene circuits

Two Cellular Resource Based Models Linking Growth and Parts Characteristics Aids the Study and Optimization of Synthetic Gene Circuits

Characterization and mitigation of gene expression burden in mammalian cells

Orthogonality and Burdens of Heterologous AND Gate Gene Circuits in E. coli

Synthesising Executable Gene Regulatory Networks from Single-cell Gene Expression Data

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

Model-guided gene circuit design for engineering genetically stable cell populations in diverse applications

Understanding resource competition to achieve predictable synthetic gene expression in eukaryotes

A resource dependent protein synthesis model for evaluating synthetic circuits

Interdependence of cell growth and gene expression: origins and consequences

Simplifications and approximations in a single-gene circuit modeling

A bottom-up approach to gene regulation

Internal workings of a genetic circuit

Generative modeling and latent space arithmetics predict single-cell perturbation response across cell types, studies and species

Evaluating the Contribution of Model Complexity in Predicting Robustness in Synthetic Genetic Circuits

In silico analysis of metabolic burden effects on a multicellular integral controller

A Simulation Study on Gene Expression Regulation Via Stochastic Model

Limits to the inference of gene regulation from bulk tissue expression data