Constructing gene co-expression networks and predicting functions of unknown genes by random matrix theory

Abstract

BackgroundLarge-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.ResultsApplication 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.ConclusionOur 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.

Fungal Communities Respond to Long-Term CO 2 Elevation by Community Reassembly

“…The RMT identifies the threshold by observing a transition point of nearest-neighbor spacing distribution of eigenvalues from Gaussian to Poisson distribution, which are two universal extreme distributions (36). The RMT-based approach is a reliable and robust tool for network construction and has been successfully applied to construct various networks, including gene regulatory networks (49)(50)(51)(52)(53), functional molecular ecological networks (36), and phylogenetic molecular ecological networks (37). Second, the same cutoff of 0.78 was applied to construct co-occurrence networks for fungal communities at aCO 2 and eCO 2 , with the purpose of comparing between different networks.…”

The RavA-ViaA Chaperone-Like System Interacts with and Modulates the Activity of the Fumarate Reductase Respiratory Complex

“…These data were available in the form of Robust Multi Array (RMA) normalized profiles [53], which enables the direct comparison of the expression profiles of different protein-encoding genes across multiple experimental conditions. The Pearson 25 correlations, used for comparing the similarity of expression profiles, were computed for all 4,125 genes present on the Affymetrix chip against both ravA and viaA. This allowed the identification of genes that exhibit the most similar expression profiles to the seed set of genes.…”

A Key Role for Apoplastic H2O2 in Norway Spruce Phenolic Metabolism

“…The resulting network was then analyzed using RmtGeneNet (Luo et al, 2007) to determine the weight threshold for which its structure deviated from a Poisson to a Gaussian distribution, as assessed by a x 2 goodness-of-fit test. This allowed for the selection of the core of the network.…”