Abstract:Motivation: Cell function is regulated by gene regulatory networks (GRNs) defined by protein-mediated interaction between constituent genes. Despite advances in experimental techniques, we can still measure only a fraction of the processes that govern GRN dynamics. To infer the properties of GRNs using partial observation, unobserved sequential processes can be replaced with distributed time delays, yielding non-Markovian models. Inference methods based on the resulting model suffer from the curse of dimensionality. Results: We develop a simulation-based Bayesian MCMC method employing an approximate likelihood for the efficient and accurate inference of GRN parameters when only some of their products are observed. We illustrate our approach using a two-step activation model: an activation signal leads to the accumulation of an unobserved regulatory protein, which triggers the expression of observed fluorescent proteins. With prior information about observed fluorescent protein synthesis, our method successfully infers the dynamics of the unobserved regulatory protein. We can estimate the delay and kinetic parameters characterizing target regulation including transcription, translation, and target searching of an unobserved protein from experimental measurements of the products of its target gene. Our method is scalable and can be used to analyze non-Markovian models with hidden components. Availability and implementation: Our code is implemented in R and is freely available with a simple example data at https://github.com/Mathbiomed/SimMCMC.

Explicitly Solvable Continuous-time Inference for Partially Observed Markov Processes

The Expectation-Maximization Algorithm for Continuous-time Hidden Markov Models

Foundation Inference Models for Markov Jump Processes

CAUSAL INFERENCE FOR CONTINUOUS-TIME PROCESSES WHEN COVARIATES ARE OBSERVED ONLY AT DISCRETE TIMES

Probabilistic Model Checking for Continuous Time Markov Chains via Sequential Bayesian Inference

CTMCs with Imprecisely Timed Observations

Markov Chain Monte Carlo for Continuous-Time Switching Dynamical Systems

Likelihood free inference for Markov processes: a comparison

Markov Decision Processes with Incomplete Information and Semi-Uniform Feller Transition Probabilities

A Finite Mixture Hidden Markov Model for Intermittently Observed Disease Process with Heterogeneity and Partially Known Disease Type

Variational Inference for Continuous-Time Switching Dynamical Systems

Forward-Backward Latent State Inference for Hidden Continuous-Time semi-Markov Chains

Continuous-Time Higher Order Markov Chains: Formulation and Parameter Estimation

Entropic Matching for Expectation Propagation of Markov Jump Processes

State and parameter estimation from exact partial state observation in stochastic reaction networks

Continuous time Stochastic optimal control under discrete time partial observations

Inferring delays in partially observed gene regulation processes

Inference from gated first-passage times

Expectation propagation for continuous time stochastic processes

Probabilistic Modeling for Sequences of Sets in Continuous-Time

A joint model for multistate disease processes and random informative observation times, with applications to electronic medical records data