Abstract:Background Large-scale sequencing of entire genomes has ushered in a new age in biology. One of the next grand challenges is to dissect the cellular networks consisting of many individual functional modules. Defining co-expression networks without ambiguity based on genome-wide microarray data is difficult and current methods are not robust and consistent with different data sets. This is particularly problematic for little understood organisms since not much existing biological knowledge can be exploited for determining the threshold to differentiate true correlation from random noise. Random matrix theory (RMT), which has been widely and successfully used in physics, is a powerful approach to distinguish system-specific, non-random properties embedded in complex systems from random noise. Here, we have hypothesized that the universal predictions of RMT are also applicable to biological systems and the correlation threshold can be determined by characterizing the correlation matrix of microarray profiles using random matrix theory. Results Application of random matrix theory to microarray data of S. oneidensis , E. coli , yeast, A. thaliana , Drosophila , mouse and human indicates that there is a sharp transition of nearest neighbour spacing distribution (NNSD) of correlation matrix after gradually removing certain elements insider the matrix. Testing on an in silico modular model has demonstrated that this transition can be used to determine the correlation threshold for revealing modular co-expression networks. The co-expression network derived from yeast cell cycling microarray data is supported by gene annotation. The topological properties of the resulting co-expression network agree well with the general properties of biological networks. Computational evaluations have showed that RMT approach is sensitive and robust. Furthermore, evaluation on sampled expression data of an in silico modular gene system has showed that under-sampled expressions do not affect the recovery of gene co-expression network. Moreover, the cellular roles of 215 functionally unknown genes from yeast, E. coli and S. oneidensis are predicted by the gene co-expression networks using guilt-by-association principle, many of which are supported by existing information or our experimental verification, further demonstrating the reliability of this approach for gene function prediction. Conclusion Our rigorous analysis of gene expression microarray profiles using RMT has showed that the transition of NNSD of correlation matrix of microarray profile provides a profound theoretical criterion to determine the correlation threshold for identifying gene co-expression networks.

A tensor-based Markov chain method for module identification from multiple networks

Multiple Networks Modules Identification by a Multi-Dimensional Markov Chain Method

ModularBoost: an Efficient Network Inference Algorithm Based on Module Decomposition.

Towards Identification of Human Disease Phenotype-Genotype Association via a Network-Module Based Method

One-to-Many Node Matching Between Complex Networks

Hierarchical Modular Structure Identification with Its Applications in Gene Coexpression Networks

Using Matrix of Thresholding Partial Correlation Coefficients to Infer Regulatory Network

Protein functional module identification method combining topological features and gene expression data

Comparison of Modularization Methods in Application to Different Biological Networks

MUMI: Multitask Module Identification for Biological Networks

Integrating data and knowledge to identify functional modules of genes: a multilayer approach

Co-regulated Gene Module Detection for Time Series Gene Expression Data

Constructing Gene Co-Expression Networks and Predicting Functions of Unknown Genes by Random Matrix Theory

A New Approach to Identify Functional Modules Using Random Matrix Theory

Gene expression complex networks: synthesis, identification, and analysis

A comprehensive evaluation of module detection methods for gene expression data

A Tensor-based Markov Chain Model for Heterogeneous Information Network Collective Classification

ClustEx2: Gene Module Identification Using Density-Based Network Hierarchical Clustering

Identification of functional modules using network topology and high-throughput data

Simultaneous clustering and estimation of networks in multiple graphical models

HetFCM: functional co-module discovery by heterogeneous network co-clustering