Abstract:Abstract Background A cell exhibits a variety of responses to internal and external cues. These responses are possible, in part, due to the presence of an elaborate gene regulatory network (GRN) in every single cell. In the past 20 years, many groups worked on reconstructing the topological structure of GRNs from large-scale gene expression data using a variety of inference algorithms. Insights gained about participating players in GRNs may ultimately lead to therapeutic benefits. Mutual information (MI) is a widely used metric within this inference/reconstruction pipeline as it can detect any correlation (linear and non-linear) between any number of variables ( n -dimensions). However, the use of MI with continuous data (for example, normalized fluorescence intensity measurement of gene expression levels) is sensitive to data size, correlation strength and underlying distributions, and often requires laborious and, at times, ad hoc optimization. Results In this work, we first show that estimating MI of a bi- and tri-variate Gaussian distribution using k -nearest neighbor (kNN) MI estimation results in significant error reduction as compared to commonly used methods based on fixed binning. Second, we demonstrate that implementing the MI-based kNN Kraskov–Stoögbauer–Grassberger (KSG) algorithm leads to a significant improvement in GRN reconstruction for popular inference algorithms, such as Context Likelihood of Relatedness (CLR). Finally, through extensive in-silico benchmarking we show that a new inference algorithm CMIA (Conditional Mutual Information Augmentation), inspired by CLR, in combination with the KSG-MI estimator, outperforms commonly used methods. Conclusions Using three canonical datasets containing 15 synthetic networks, the newly developed method for GRN reconstruction—which combines CMIA, and the KSG-MI estimator—achieves an improvement of 20–35% in precision-recall measures over the current gold standard in the field. This new method will enable researchers to discover new gene interactions or better choose gene candidates for experimental validations.

A comparative study of Gaussian Graphical Model approaches for genomic data

Using Matrix of Thresholding Partial Correlation Coefficients to Infer Regulatory Network

Introducing Gaussian covariance graph models in genome-wide prediction

Modelling correlated marker effects in genome-wide prediction via Gaussian concentration graph models

Netreg: Network-Regularized Linear Models for Biological Association Studies

L1-Penalized Censored Gaussian Graphical Model

Multivariate dependence and genetic networks inference

High-Dimensional Gaussian Graphical Regression Models with Covariates

A Log-Linear Graphical Model for Inferring Genetic Networks from High-Throughput Sequencing Data

NetDiff – Bayesian model selection for differential gene regulatory network inference

Statistical Inference on High Dimensional Gaussian Graphical Regression Models

A sparse conditional Gaussian graphical model for analysis of genetical genomics data

Biological network inference using low order partial correlation

Reciprocal Graphical Models for Integrative Gene Regulatory Network Analysis

Identifying Graphical Models

Gene regulation network inference using k-nearest neighbor-based mutual information estimation: revisiting an old DREAM

Extraction and integration of genetic networks from short-profile omic datasets

Information theoretic approaches for inference of biological networks from continuous-valued data

Inferring genetic networks: An information theoretic approach

A computational framework for gene regulatory network inference that combines multiple methods and datasets

Variational inference for sparse network reconstruction from count data