Efficient inference of dynamic gene regulatory networks using discrete penalty

TL;DR

This paper introduces a scalable method for inferring dynamic gene regulatory networks from high-dimensional data using a discrete $ ext{L}_0$ penalty, improving accuracy and interpretability over traditional methods.

Contribution

It presents a novel framework employing a discrete $ ext{L}_0$ penalty for unbiased, sparse, and scalable inference of dynamic gene regulatory networks, accommodating complex sample structures.

Findings

Successfully validated on synthetic benchmarks.

Reconstructed gene networks in glioblastoma data.

Mapped network rewiring along hypoxia gradients.

Abstract

Gene regulatory networks (GRNs) orchestrate cellular decision making and survival strategies. Inferring the structure of these networks from high-dimensional transcriptomics data is a central challenge in systems biology. Traditional approaches to GRN inference, such as the graphical lasso and its joint extensions, rely on $\ell_1$ penalty to induce sparsity but can bias network recovery and require extensive hyperparameter tuning. Here, we present a scalable framework for the joint inference of dynamic GRNs using a discrete $\ell_0$ penalty, enabling direct and unbiased control over network sparsity. Leveraging recent algorithmic advances, we efficiently solve the resulting mixed-integer optimization problem for populations structured as arbitrary tree hypergraphs, accommodating both continuous and categorical distinctions among biological samples. After validating our method on synthetic benchmarks, we apply it to single-cell and spatial transcriptomics data from glioblastoma (GBM) patient tumors. Our approach reconstructs gene networks across tumor clusters, maps network rewiring along hypoxia gradients, and reveals niche-specific differences between primary and recurrent tumors. By providing a robust and interpretable tool for GRN inference in complex tissues, our work facilitates high-resolution dissection of tumor heterogeneity and adaptation, with broad applicability to emerging large-scale transcriptomic datasets.