One model fits all: Combining inference and simulation of gene regulatory networks

PLoS Comput Biol. 2023 Mar 27;19(3):e1010962. doi: 10.1371/journal.pcbi.1010962. eCollection 2023 Mar.

Abstract

The rise of single-cell data highlights the need for a nondeterministic view of gene expression, while offering new opportunities regarding gene regulatory network inference. We recently introduced two strategies that specifically exploit time-course data, where single-cell profiling is performed after a stimulus: HARISSA, a mechanistic network model with a highly efficient simulation procedure, and CARDAMOM, a scalable inference method seen as model calibration. Here, we combine the two approaches and show that the same model driven by transcriptional bursting can be used simultaneously as an inference tool, to reconstruct biologically relevant networks, and as a simulation tool, to generate realistic transcriptional profiles emerging from gene interactions. We verify that CARDAMOM quantitatively reconstructs causal links when the data is simulated from HARISSA, and demonstrate its performance on experimental data collected on in vitro differentiating mouse embryonic stem cells. Overall, this integrated strategy largely overcomes the limitations of disconnected inference and simulation.

Publication types

  • Research Support, Non-U.S. Gov't

MeSH terms

  • Algorithms*
  • Animals
  • Computer Simulation
  • Gene Expression Profiling / methods
  • Gene Regulatory Networks* / genetics
  • Mice

Grants and funding

This work was supported by funding from French agency ANR (SingleStatOmics; ANR-18-CE45-0023-03) to OG. The funders had no role in study design, data collection and analysis, decision to publish, or preparation of the manuscript.