Combining tree-based and dynamical systems for the inference of gene regulatory networks

Reconstructing the topology of gene regulatory networks (GRNs) from time series of gene expression data remains an important open problem in computational systems biology. Existing GRN inference algorithms face one of two limitations: model-free methods are scalable but suffer from a lack of interpretability, and cannot in general be used for out of sample predictions. On the other hand, model-based methods focus on identifying a dynamical model of the system; these are clearly interpretable and can be used for predictions, however they rely on strong assumptions and are typically very demanding computationally.

Here, we propose a new hybrid approach for GRN inference, called Jump3, exploiting time series of expression data. Jump3 is based on a formal on/off model of gene expression, but uses a non-parametric procedure based on decision trees (called “jump trees”) to reconstruct the GRN topology, allowing the inference of networks of hundreds of genes. We show the good performance of Jump3 on in silico and synthetic networks, and applied the approach to identify regulatory interactions activated in the presence of interferon gamma.

Our MATLAB implementation of Jump3 is available at http://www.montefiore.ulg.ac.be/~huynh-thu/software.html.