Bayesian inference of gene regulatory networks at stochastic steady state

Gene regulatory networks (GRNs) form the regulatory backbone that coordinates gene expression. The architecture of GRNs shapes their function and constrains the biochemical pathways through which information flows. Inferring the structure of regulatory interactions is thus essential for understanding biological systems and designing targeted therapies. Despite substantial progress in GRN inference, most approaches-from statistical methods to deep learning-do not take into account fundamental biochemical processes that drive regulatory dynamics. To address this shortcoming, here, we present a novel Bayesian inference approach based on using the chemical Langevin equation as a model of gene expression dynamics at stochastic equilibrium. Interactions in GRNs are sparse, and we thus use a regularized horseshoe prior enabling selective shrinkage of unsupported interactions while identifying strong regulatory edges. We evaluate our method using synthetic gene expression data, allowing for benchmarking against a known ground truth. Our approach allows us to infer kinetic parameters, identify network structure and infer regulatory cycles without the need to observe transient dynamics. This Bayesian alternative to current methods thus provides both biological interpretability and structural identifiability in GRN inference.