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.

Discriminate the falsely predicted protein-coding genes in Aeropyrum Pernix K1 genome based on graphical representation

Identify protein-coding genes in the genomes of Aeropyrum pernix K1 and Chlorobium tepidum TLS.

Discovering putative peptides encoded from non-coding RNAs in ribosome profiling data of Arabidopsis thaliana.

Gene Prediction by the Noise-Assisted MEMD and Wavelet Transform for Identifying the Protein Coding Regions

Re-prediction of protein-coding genes in the genome of Amsacta moorei entomopoxvirus.

Gene Recognition Based on Nucleotide Distribution of ORFs in a Hyper-Thermophilic Crenarchaeon, Aeropyrum Pernix K1.

Theoretical Prediction and Experimental Verification of Protein-Coding Genes in Plant Pathogen Genome Agrobacterium tumefaciens Strain C58

Flawed machine-learning confounds coding sequence annotation

Accurate Prediction of Human Essential Genes Using Only Nucleotide Composition and Association Information

Gene Re-annotation in Genome of the ExtremophilePyrobaculum Aerophilumby Using Bioinformatics Methods

Re-Annotation of Protein-Coding Genes in the Genome of Saccharomyces Cerevisiae Based on Support Vector Machines

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

Prediction of Protein-Coding Small ORFs in Multi-Species Using Integrated Sequence-Derived Features and the Random Forest Model.

Recognizing Short Coding Sequences of Prokaryotic Genome Using a Novel Iteratively Adaptive Sparse Partial Least Squares Algorithm.

Analysis and Prediction of Exon, Intron, Intergenic Region and Splice Sites for A. Thaliana and C. Elegans Genomes

Inferring an Organism-Specific Optimal Threshold for Predicting Protein Coding Regions in Eukaryotes Based on a Bootstrapping Algorithm

Comprehensive evaluation of protein-coding sORFs prediction based on a random sequence strategy

The shrinking human protein coding complement: are there now fewer than 20,000 genes?

A Novel Approach to Identify Protein Coding Domains by Sampling Binary Profiles from Genome

Iesgene-Zcpseknc: Identify Essential Genes Based on Z Curve Pseudo $k$ -Tuple Nucleotide Composition.

GeneAlign: a Coding Exon Prediction Tool Based on Phylogenetical Comparisons.