Automated discovery of functional generality of human gene expression programs.
An important research problem in computational biology is the identification of expression programs, sets of co-expressed genes orchestrating normal or pathological processes, and the characterization of the functional breadth of these programs. The use of human expression data compendia for discovery of such programs presents several challenges including cellular inhomogeneity within samples, genetic and environmental variation across samples, uncertainty in the numbers of programs and sample populations, and temporal behavior. We developed GeneProgram, a new unsupervised computational framework based on Hierarchical Dirichlet Processes that addresses each of the above challenges. GeneProgram uses expression data to simultaneously organize tissues into groups and genes into overlapping programs with consistent temporal behavior, to produce maps of expression programs, which are sorted by generality scores that exploit the automatically learned groupings. Using synthetic and real gene expression data, we showed that GeneProgram outperformed several popular expression analysis methods. We applied GeneProgram to a compendium of 62 short time-series gene expression datasets exploring the responses of human cells to infectious agents and immune-modulating molecules. GeneProgram produced a map of 104 expression programs, a substantial number of which were significantly enriched for genes involved in key signaling pathways and/or bound by NF-kappaB transcription factors in genome-wide experiments. Further, GeneProgram discovered expression programs that appear to implicate surprising signaling pathways or receptor types in the response to infection, including Wnt signaling and neurotransmitter receptors. We believe the discovered map of expression programs involved in the response to infection will be useful for guiding future biological experiments; genes from programs with low generality scores might serve as new drug targets that exhibit minimal "cross-talk," and genes from high generality programs may maintain common physiological responses that go awry in disease states. Further, our method is multipurpose, and can be applied readily to novel compendia of biological data.
Gerber
Georg K
GK
Department of Computer Science and Electrical Engineering, Massachusetts Institute of Technology, Cambridge, Massachusetts, United States of America.
Dowell
Robin D
RD
Jaakkola
Tommi S
TS
Gifford
David K
DK
D000465:Algorithms; D001185:Artificial Intelligence; D003198:Computer Simulation; D015870:Gene Expression / Q000502:physiology*; D020869:Gene Expression Profiling / Q000379:methods*; D008954:Models, Biological; D005810:Multigene Family / Q000502:physiology*; D020411:Oligonucleotide Array Sequence Analysis / Q000379:methods; D010363:Pattern Recognition, Automated / Q000379:methods; D020543:Proteome / Q000378:metabolism*; D015398:Signal Transduction / Q000502:physiology*; D012984:Software
D020543:Proteome
Bioinformatics. 2005 Jun;21 Suppl 1:i144-51
15961451
J Comput Biol. 2003;10(3-4):341-56
12935332
Nucleic Acids Res. 2006 Jan 1;34(Database issue):D354-7
16381885
Proc Natl Acad Sci U S A. 2003 Mar 18;100(6):3351-6
12631705
Nucleic Acids Res. 2006 Jan 1;34(Database issue):D173-80
16381840
Bioinformatics. 2005 Aug 1;21(15):3286-93
15919724
Nature. 2005 Sep 15;437(7057):417-21
16163358
Proc Natl Acad Sci U S A. 1998 Dec 8;95(25):14863-8
9843981
Proc Natl Acad Sci U S A. 2002 Jan 22;99(2):972-7
11805339
BMC Bioinformatics. 2007 Jan 10;8:9
17214880
J Comput Biol. 2005 Sep;12(7):909-27
16201912
Nat Genet. 2004 Oct;36(10):1090-8
15448693
Proc Natl Acad Sci U S A. 2004 Mar 23;101(12):4164-9
15016911
Nature. 1999 Oct 21;401(6755):788-91
10548103
Bioinformatics. 2001 Jun;17(6):495-508
11395426
BMC Bioinformatics. 2005 Sep 21;6:232
16176576
Genome Res. 2003 Jul;13(7):1706-18
12840046
Proc Natl Acad Sci U S A. 2005 Apr 12;102(15):5310-1
15811940
Proc Int Conf Intell Syst Mol Biol. 2000;8:93-103
10977070
Bioinformatics. 2003 Oct;19 Suppl 2:ii196-205
14534190
Proc Natl Acad Sci U S A. 2004 Apr 6;101 Suppl 1:5228-35
14872004
Bioinformatics. 2004 May 22;20(8):1222-32
14871871
Proc Natl Acad Sci U S A. 2004 Apr 20;101(16):6062-7
15075390
Proc Natl Acad Sci U S A. 2006 Jan 17;103(3):696-701
16407142
Blood. 2000 Oct 1;96(7):2323-8
11001878
IEEE/ACM Trans Comput Biol Bioinform. 2004 Jan-Mar;1(1):24-45
17048406
J Allergy Clin Immunol. 2007 Feb;119(2):472-81
17125827
Proc Natl Acad Sci U S A. 2000 Aug 29;97(18):10101-6
10963673
BMC Bioinformatics. 2006 Feb 17;7:78
16503973
Pac Symp Biocomput. 2001;:422-33
11262961
Nat Rev Microbiol. 2005 Apr;3(4):281-94
15806094
Theor Biol Med Model. 2006 Apr 05;3:18
16597340
Lancet. 2004 Jan 17;363(9404):203-9
14738793
Proc Natl Acad Sci U S A. 2006 Apr 4;103(14):5454-9
16569699
Blood. 2006 Aug 1;108(3):965-73
16601243
Nucleic Acids Res. 2006 Jan 1;34(Database issue):D322-6
16381878
Genome Biol. 2005;6(3):R22
15774023
BMC Bioinformatics. 2006 Apr 05;7:191
16597342
Proc Natl Acad Sci U S A. 2002 Feb 5;99(3):1503-8
11805289
J Immunol. 2003 May 15;170(10):5203-9
12734368
Proc Natl Acad Sci U S A. 2006 Apr 11;103(15):5899-904
16595631
Bioinformatics. 2002;18 Suppl 1:S136-44
12169541
Nat Rev Mol Cell Biol. 2004 May;5(5):392-401
15122352
Cell Signal. 2007 Feb;19(2):251-60
16996249
Science. 2001 Oct 26;294(5543):870-5
11679675
J Mol Neurosci. 2006;30(1-2):107-8
17192650
J Comput Biol. 2005 Jul-Aug;12(6):882-93
16108723
Curr Opin Genet Dev. 2001 Oct;11(5):547-53
11532397
Bioinformatics. 2004 Nov 1;20(16):2493-503
15130923
Proc Natl Acad Sci U S A. 2002 Nov 12;99(23):15136-41
12411577