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.

Inferring stochastic regulatory networks from perturbations of the non-equilibrium steady state

Inferring Regulatory Networks by Combining Perturbation Screens and Steady State Gene Expression Profiles

Stochastic Stability of Genetic Networks with Disturbance Attenuation.

Data-driven approximation for extracting the transition dynamics of a genetic regulatory network with non-Gaussian Lévy noise

Learning a nonlinear dynamical system model of gene regulation: A perturbed steady-state approach

Gene Regulatory Network Inference from Perturbed Time-Series Expression Data Via Ordered Dynamical Expansion of Non-Steady State Actors.

Gene Regulatory Network Inference from Multifactorial Perturbation Data

Gene network reconstruction from transcriptional dynamics under kinetic model uncertainty: a case for the second derivative

Identifying Latent Dynamic Components in Biological Systems

Modeling Stochasticity and Variability in Gene Regulatory Networks

Stochastic Modeling of Biophysical Responses to Perturbation

A Posterior Probability Approach for Gene Regulatory Network Inference in Genetic Perturbation Data

Inference and Uncertainty Quantification of Stochastic Gene Expression via Synthetic Models

Dynamic probabilistic threshold networks to infer signaling pathways from time-course perturbation data

Network Inference Using Steady State Data and Goldbeter-Koshland Kinetics

Stochastic fluctuations can reveal the feedback signs of gene regulatory networks at the single-molecule level

Extending the linear-noise approximation to biochemical systems influenced by intrinsic noise and slow lognormally distributed extrinsic noise

Inference of gene regulatory networks based on nonlinear ordinary differential equations

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

Inferring genetic networks: An information theoretic approach

Effects of Small Random Perturbations in the Extended Glass–Kauffman Model of Gene Regulatory Networks