Constructing gene co-expression networks and predicting functions of unknown genes by random matrix theory.
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.


CONCLUSIONS
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.
Luo
Feng
F
Environmental Sciences Division, Oak Ridge National Laboratory, Oak Ridge, Tennessee 37831, USA. luofeng@cs.clemson.edu
Yang
Yunfeng
Y
Zhong
Jianxin
J
Gao
Haichun
H
Khan
Latifur
L
Thompson
Dorothea K
DK
Zhou
Jizhong
J
D000465:Algorithms; D003198:Computer Simulation; D020869:Gene Expression Profiling / Q000379:methods*; D008954:Models, Biological; D005810:Multigene Family / Q000502:physiology*; D020411:Oligonucleotide Array Sequence Analysis / Q000379:methods*; D020543:Proteome / Q000378:metabolism*; D015398:Signal Transduction / Q000502:physiology*
D020543:Proteome
Bioinformatics. 2001 Jun;17(6):520-5
11395428
Nature. 2003 Oct 16;425(6959):686-91
14562095
Bioinformatics. 2002 May;18(5):735-46
12050070
J Comput Biol. 2000;7(3-4):601-20
11108481
Bioinformatics. 2002 Nov;18(11):1486-93
12424120
EMBO J. 2002 Oct 15;21(20):5539-47
12374754
Nature. 2006 Mar 30;440(7084):637-43
16554755
Nat Genet. 2005 Apr;37(4):382-90
15778709
Stat Appl Genet Mol Biol. 2005;4:Article17
16646834
Proc Natl Acad Sci U S A. 2003 Dec 23;100(26):15522-7
14673099
Phys Rev E Stat Nonlin Soft Matter Phys. 2003 Mar;67(3 Pt 1):031902
12689096
Pac Symp Biocomput. 1999;:17-28
10380182
Pac Symp Biocomput. 2000;:418-29
10902190
Phys Rev E Stat Nonlin Soft Matter Phys. 2006 Mar;73(3 Pt 1):031924
16605575
Nucleic Acids Res. 2004 Jan 1;32(Database issue):D311-4
14681421
Pac Symp Biocomput. 1998;:18-29
9697168
Science. 2002 May 3;296(5569):910-3
11988575
Phys Life Rev. 2005 Mar;2(1):65-88
20416858
Nat Rev Genet. 2004 Feb;5(2):101-13
14735121
Mol Cancer Ther. 2003 Feb;2(2):199-205
12589037
Science. 2003 Dec 12;302(5652):1956-60
14671301
Genome Res. 2003 Dec;13(12):2568-76
14656964
Proc Natl Acad Sci U S A. 2002 Oct 1;99(20):12783-8
12196633
Science. 2003 Oct 10;302(5643):240-1
14551426
Science. 2001 May 4;292(5518):929-34
11340206
Nat Genet. 2002 May;31(1):64-8
11967538
Nat Genet. 2003 Jun;34(2):166-76
12740579
Nat Biotechnol. 2003 Nov;21(11):1337-42
14555958
Science. 2003 Jul 4;301(5629):102-5
12843395
Science. 2002 Aug 30;297(5586):1551-5
12202830
Curr Opin Struct Biol. 2003 Apr;13(2):193-202
12727512
Nature. 2002 Jan 10;415(6868):141-7
11805826
Science. 2002 Sep 27;297(5590):2270-5
12351791
J Bacteriol. 2004 Dec;186(24):8385-400
15576789
Science. 2003 Sep 26;301(5641):1866-7
14512615
Phys Rev Lett. 2003 Nov 7;91(19):198104
14611622
J Bacteriol. 2004 Nov;186(22):7796-803
15516594
Nature. 1999 Dec 2;402(6761 Suppl):C47-52
10591225
Mol Biol Cell. 2000 Dec;11(12):4241-57
11102521
Proc Natl Acad Sci U S A. 2002 Apr 30;99(9):6163-8
11983907
Mol Biol Cell. 1998 Dec;9(12):3273-97
9843569
Nat Genet. 2002 Aug;31(4):370-7
12134151
Genome Res. 2004 Aug;14(8):1654-63
15289483
Nature. 2003 Mar 20;422(6929):297-302
12646919
Science. 2003 Oct 10;302(5643):249-55
12934013
J Bacteriol. 2002 Oct;184(19):5513-7
12218041
Proc Natl Acad Sci U S A. 2001 Jan 30;98(3):1176-81
11158614
Phys Rev B Condens Matter. 1993 Dec 15;48(23):16979-16985
10008297