MICRAT: a novel algorithm for inferring gene regulatory networks using time series gene expression data

Abstract

BackgroundReconstruction of gene regulatory networks (GRNs), also known as reverse engineering of GRNs, aims to infer the potential regulation relationships between genes. With the development of biotechnology, such as gene chip microarray and RNA-sequencing, the high-throughput data generated provide us with more opportunities to infer the gene-gene interaction relationships using gene expression data and hence understand the underlying mechanism of biological processes. Gene regulatory networks are known to exhibit a multiplicity of interaction mechanisms which include functional and non-functional, and linear and non-linear relationships. Meanwhile, the regulatory interactions between genes and gene products are not spontaneous since various processes involved in producing fully functional and measurable concentrations of transcriptional factors/proteins lead to a delay in gene regulation. Many different approaches for reconstructing GRNs have been proposed, but the existing GRN inference approaches such as probabilistic Boolean networks and dynamic Bayesian networks have various limitations and relatively low accuracy. Inferring GRNs from time series microarray data or RNA-sequencing data remains a very challenging inverse problem due to its nonlinearity, high dimensionality, sparse and noisy data, and significant computational cost, which motivates us to develop more effective inference methods.ResultsWe developed a novel algorithm, MICRAT (Maximal Information coefficient with Conditional Relative Average entropy and Time-series mutual information), for inferring GRNs from time series gene expression data. Maximal information coefficient (MIC) is an effective measure of dependence for two-variable relationships. It captures a wide range of associations, both functional and non-functional, and thus has good performance on measuring the dependence between two genes. Our approach mainly includes two procedures. Firstly, it employs maximal information coefficient for constructing an undirected graph to represent the underlying relationships between genes. Secondly, it directs the edges in the undirected graph for inferring regulators and their targets. In this procedure, the conditional relative average entropies of each pair of nodes (or genes) are employed to indicate the directions of edges. Since the time delay might exist in the expression of regulators and target genes, time series mutual information is combined to cooperatively direct the edges for inferring the potential regulators and their targets. We evaluated the performance of MICRAT by applying it to synthetic datasets as well as real gene expression data and compare with other GRN inference methods. We inferred five 10-gene and five 100-gene networks from the DREAM4 challenge that were generated using the gene expression simulator GeneNetWeaver (GNW). MICRAT was also used to reconstruct GRNs on real gene expression data including part of the DNA-damaged response pathway (SOS DNA repair network) and experimental dataset in E. Coli. The results showe...

Learning causal networks using inducible transcription factors and transcriptome‐wide time series

“…Integration of multiple ‘omic technologies combined with time‐series measurements can help identify direct functional interactions to elucidate GRNs, as was done in a recent study that combined RNA‐seq, NET‐seq, and ChIP‐seq to identify a core regulon for Hsf1 in yeast (Solís et al , ). Finally, there is a growing literature of computational methods for reconstructing GRNs from high‐throughput data (Pe'er et al , ; Markowetz et al , ; Stolovitzky et al , ; Marbach et al , , ,b; Yip et al , ; Yang et al , ). The Dialogue on Reverse Engineering Assessment and Methods (DREAM) project, which is organized around annual challenges, provides a framework to benchmark network inference methods (Marbach et al , ).…”

Genetic Complexity of the Human Genome in Health and Disease: Basic Concepts

“…Mathematical tools [10][11][12][13][14]. Studies on the proper execution of the gene expression program for each cell as well as on structural changes of the genome, reveal complex processes, which are demonstrated by a massive body of data regarding genetic information and networks.…”

Gat2Get: A Novel Approach to Infer Gene Regulatory Network from Gene Activity using Dynamic Bayesian Network learning

“…The Gene Ontology (GO) is a biological repository specific to genes. It is designed to encapsulate the known relationships between biological terms and all genes that are related to these terms [3]. An expanded version of GO's proposed framework for the field of gene regulation is the Gene Regulation Ontology (GRO).…”