- Methodology article
- Open Access
FunSpec: a web-based cluster interpreter for yeast
© Robinson et al; licensee BioMed Central Ltd. 2002
Received: 8 September 2002
Accepted: 13 November 2002
Published: 13 November 2002
For effective exposition of biological information, especially with regard to analysis of large-scale data types, researchers need immediate access to multiple categorical knowledge bases and need summary information presented to them on collections of genes, as opposed to the typical one gene at a time.
We present here a web-based tool (FunSpec) for statistical evaluation of groups of genes and proteins (e.g. co-regulated genes, protein complexes, genetic interactors) with respect to existing annotations (e.g. functional roles, biochemical properties, localization). FunSpec is available online at http://funspec.med.utoronto.ca
FunSpec is helpful for interpretation of any data type that generates groups of related genes and proteins, such as gene expression clustering and protein complexes, and is useful for predictive methods employing "guilt-by-association."
A large proportion of the online resources compiling information about gene and protein functions present information one-gene-at-a-time [1–3]. Interpretation of large-scale data types, however, often requires a rapid means to access and present information on large numbers of genes simultaneously [4, 5]. Recent tools that provide annotation of gene expression data involve text mining of biomedical literature  and the integration of a wide variety of data types including Pfam domains, SWISS-PROT keywords and known pathways (e.g. DRAGON View; see ). We have developed an online resource called FunSpec (for Functional Specification), which uses information from public databases to evaluate lists of yeast genes to determine if they are enriched for particular attributes, using a well-accepted statistical model. The data tables underlying FunSpec are updated regularly from major public-domain knowledge sources, and from our analyses of large-scale data sets. Hence, FunSpec is a convenient, current, and comprehensive resource for any application requiring interpretation of groups of genes, including clustering of microarray expression data, lists of protein IDs from mass spectrometry, and genetic interactions.
Online database and published data repositories
Yeast, perhaps the most thoroughly-studied eukaryote, already has a number of comprehensive online databases containing functional annotation information for close to two-thirds of all known or predicted open reading frames (ORFs) in the genome, and numerous large-scale experimental results covering practically every gene/protein (eg., [8, 9]). In many cases there are links between the data sources, though the majority of these are at the individual gene level. For example, the Comprehensive Yeast Genome Database at the Munich Information Center for Protein Sequences (MIPS) lists five catalogues encompassing a curation of several thousand publications and spanning functional classification, known protein complexes, protein classes, mutant phenotypes and subcellular localization . Similarly, the Gene Ontology (GO) Consortium organizes the same information under the categorizations: Biological Process, Molecular Function and Cellular Component . Other online knowledge bases encompass computationally-derived gene properties, such as domain and motif information (e.g. SMART , Pfam ). In addition, there are now mass stores of high-throughput biological data accumulating from experimental results such as synthetic lethal , protein interactions by yeast two hybrid [14, 15] or mass spectrometry [16, 17], immunolocalization , large-scale phenotypic analysis , and microarray analysis (eg., ).
FunSpec-searchable public databases
Munich Information Center for Protein Sequences (MIPS)
Gene Ontology (GO)
Simple Modular Architecture Research Tool
Protein families database of alignments and HMMs
FunSpec-searchable published datasets
Protein-Protein Interactions via Yeast Two Hybrid Method
Protein-Protein Interactions via Yeast Two Hybrid Method (core + full datasets)
Synthetic Lethal Interactions via Synthetic Genetic Array analysis
Protein Complex Identification using TAP tagging and mass spectrometry
Protein Complex Identification using Flag tagging and mass spectrometry
Interpretation of clusters
Upon collection of data resulting in clusters of genes/proteins having similar properties, responses, or other associations, often a researcher may wish to summarize and interpret the current state of knowledge of genes in such a cluster. A cluster of genes can be said to be "functionally enriched" for an attribute if the proportion of genes within the cluster known to have that attribute exceeds the number that could reasonably be expected from random chance. For example, 215 of the 6267 genes in yeast  are known to function in ribosome biogenesis. If, in a clustering analysis of gene expression profiles, one finds a cluster of 100 genes to be co-regulated and of those, 60 are known to function in ribosome biogenesis, then intuitively the cluster is enriched for ribosome biogenesis (60% in the cluster versus 3% in the genome). The degree of functional enrichment for a given cluster and category can be quantitatively assessed by the hypergeometric distribution . For each category, the probability (p-value) of observing such an overlap by chance is calculated as:
where G is the size of the genome, C is the number of genes in the genome having that attribute, n is the size of the query cluster, of which k are known to possess the attribute.
If this probability is sufficiently low for a given category, then a list of genes (e.g., a cluster) is said to be enriched for that attribute (see examples below). FunSpec inputs a list of genes and computes the hypergeometric P-value in many knowledge sources (MIPS, GO, SMART and Pfam domains, published protein complexes, 2-hybrid interactions, and subcellular localizations). A Bonferroni correction may be applied to compensate for multiple testing over many categories of a knowledge base.
Components of purified protein complexes can also be interpreted using FunSpec. It is noteworthy that functional annotations, most of which are based on phenotypic rather than biochemical data, generally correlate with large protein complexes in the literature (rand index = 0.94, adjusted rand index = 0.15; rand index measures the similarity of two partitions where 0 represents no similarity and 1 is complete correspondence; the adjusted rand index is a much more conservative measure ) (data not shown). This underscores the utility of interpreting one data type with information from another [23, 24]. Subcellular localization correlates positively with membership in protein complexes, as would be expected to enable physical association (among all protein complexes in the MIPS database, 17% are enriched (P < 0.01) for a specific subcellular compartment in the Kumar experimental data , in comparison to 1% for random lists of proteins with the same size distribution) (data not shown). For newly-identified complexes, FunSpec provides not only information about localization and potential function, but also a measure of confidence regarding the biological validity of the complex.
The current FunSpec implementation utilizes the hypergeometric distribution (following ), which assumes independent tests (i.e. each gene has a single classification). In reality, the MIPS and GO categories are hierarchical. We are currently considering other statistical models to account for the hierarchical nature of the categorization, as well as alternatives to the Bonferroni correction (eg. False Discovery Rate ).
There are other methods of determining significance (e.g. binomial distribution ; representation factor ) and other means of relating clusters to functional annotation (e.g. associated terms and phrases in Medline abstracts ) but the hypergeometric distribution has proven to be the simplest to use. The advantages of FunSpec over current tools is the integration of many relevant and recent data sources and a convenient website enabling the user to explore special features of gene groupings quickly and exhaustively.
In the near future, as more functional annotation information becomes available for other organisms such as mouse and human, the FunSpec resource will be updated to accommodate the new knowledge bases on these organisms, presumably including a great deal of new experimental data. The underlying assembly of categorical annotations is useful on its own and can be downloaded from the FunSpec web site as a text file.
FunSpec takes clusters (or any collection of genes or proteins) and identifies attributes (based on prior knowledge) for which the cluster is enriched. In our own research we use FunSpec routinely for rapid interpretation of all clusters in full-genome clustering analyses, as well as for "guilt-by-association" predictive methodologies (e.g. [5, 8, 28]). We and others have also found it useful as a filter for prioritizing screening hits, interpreting genetic interactions, and selecting thresholds/boundaries for membership in protein complexes (data not shown).
The FunSpec backend is implemented in C++. Upon receiving the information from the web page, the relevant databases are queried, the statistics are calculated and the results are sent back to the user. At present, the databases are stored in flat text files.
The authors thank members of the Best Institute (University of Toronto) for testing the web implementation, and Gary Bader, Mike Tyers, and Andrew Emili for critical evaluation of this manuscript. This work was supported by Genome Canada, the CIHR, and the University of Toronto Connaught Foundation.
- Cherry JM, Adler C, Ball C, Chervitz SA, Dwight SS, Hester ET, Jia Y, Juvik G, Roe T, Schroeder M, et al.: SGD: Saccharomyces Genome Database. Nucleic Acids Res 1998, 26: 73–79. 10.1093/nar/26.1.73PubMed CentralView ArticlePubMedGoogle Scholar
- Mewes HW, Frishman D, Gruber C, Geier B, Haase D, Kaps A, Lemcke K, Mannhaupt G, Pfeiffer F, Schuller C, et al.: MIPS: a database for genomes and protein sequences. Nucleic Acids Res 2000, 28: 37–40. 10.1093/nar/28.1.37PubMed CentralView ArticlePubMedGoogle Scholar
- Costanzo MC, Crawford ME, Hirschman JE, Kranz JE, Olsen P, Robertson LS, Skrzypek MS, Braun BR, Hopkins KL, Kondu P, et al.: YPD, PombePD and WormPD: model organism volumes of the BioKnowledge library, an integrated resource for protein information. Nucleic Acids Res 2001, 29: 75–79. 10.1093/nar/29.1.75PubMed CentralView ArticlePubMedGoogle Scholar
- Greenbaum D, Luscombe NM, Jansen R, Qian J, Gerstein M: Interrelating different types of genomic data, from proteome to secretome: 'oming in on function. Genome Res 2001, 11: 1463–1468. 10.1101/gr.207401View ArticlePubMedGoogle Scholar
- Wu LF, Hughes TR, Davierwala AP, Robinson MD, Stoughton R, Altschuler SJ: Large-scale prediction of Saccharomyces cerevisiae gene function using overlapping transcriptional clusters. Nat Genet 2002, 31: 255–265. 10.1038/ng906View ArticlePubMedGoogle Scholar
- Blaschke C, Oliveros JC, Valencia A: Mining functional information associated with expression arrays. Funct Integr Genomics 2001, 1: 256–268. 10.1007/s101420000036View ArticlePubMedGoogle Scholar
- Bouton CM, Pevsner J: DRAGON View: information visualization for annotated microarray data. Bioinformatics 2002, 18: 323–324. 10.1093/bioinformatics/18.2.323View ArticlePubMedGoogle Scholar
- Hughes TR, Marton MJ, Jones AR, Roberts CJ, Stoughton R, Armour CD, Bennett HA, Coffey E, Dai H, He YD, et al.: Functional discovery via a compendium of expression profiles. Cell 2000, 102: 109–126.View ArticlePubMedGoogle Scholar
- Kumar A, Cheung KH, Ross-Macdonald P, Coelho PS, Miller P, Snyder M: TRIPLES: a database of gene function in Saccharomyces cerevisiae. Nucleic Acids Res 2000, 28: 81–84. 10.1093/nar/28.1.81PubMed CentralView ArticlePubMedGoogle Scholar
- The Gene Ontology Consortium: Creating the gene ontology resource: design and implementation. Genome Res 2001, 11: 1425–1433. 10.1101/gr.180801View ArticleGoogle Scholar
- Schultz J, Milpetz F, Bork P, Ponting CP: SMART, a simple modular architecture research tool: identification of signaling domains. Proc Natl Acad Sci U S A 1998, 95: 5857–5864. 10.1073/pnas.95.11.5857PubMed CentralView ArticlePubMedGoogle Scholar
- Bateman A, Birney E, Cerruti L, Durbin R, Etwiller L, Eddy SR, Griffiths-Jones S, Howe KL, Marshall M, Sonnhammer EL: The Pfam protein families database. Nucleic Acids Res 2002, 30: 276–280. 10.1093/nar/30.1.276PubMed CentralView ArticlePubMedGoogle Scholar
- Tong AH, Evangelista M, Parsons AB, Xu H, Bader GD, Page N, Robinson M, Raghibizadeh S, Hogue CW, Bussey H, et al.: Systematic genetic analysis with ordered arrays of yeast deletion mutants. Science 2001, 294: 2364–2368. 10.1126/science.1065810View ArticlePubMedGoogle Scholar
- Uetz P, Giot L, Cagney G, Mansfield TA, Judson RS, Knight JR, Lockshon D, Narayan V, Srinivasan M, Pochart P, et al.: A comprehensive analysis of protein-protein interactions in Saccharomyces cerevisiae. Nature 2000, 403: 623–627. 10.1038/35001009View ArticlePubMedGoogle Scholar
- Ito T, Chiba T, Ozawa R, Yoshida M, Hattori M, Sakaki Y: A comprehensive two-hybrid analysis to explore the yeast protein interactome. Proc Natl Acad Sci U S A 2001, 98: 4569–4574. 10.1073/pnas.061034498PubMed CentralView ArticlePubMedGoogle Scholar
- Gavin AC, Bosche M, Krause R, Grandi P, Marzioch M, Bauer A, Schultz J, Rick JM, Michon AM, Cruciat CM, et al.: Functional organization of the yeast proteome by systematic analysis of protein complexes. Nature 2002, 415: 141–147. 10.1038/415141aView ArticlePubMedGoogle Scholar
- Ho Y, Gruhler A, Heilbut A, Bader GD, Moore L, Adams SL, Millar A, Taylor P, Bennett K, Boutilier K, et al.: Systematic identification of protein complexes in Saccharomyces cerevisiae by mass spectrometry. Nature 2002, 415: 180–183. 10.1038/415180aView ArticlePubMedGoogle Scholar
- Kumar A, Agarwal S, Heyman JA, Matson S, Heidtman M, Piccirillo S, Umansky L, Drawid A, Jansen R, Liu Y, et al.: Subcellular localization of the yeast proteome. Genes Dev 2002, 16: 707–719. 10.1101/gad.970902PubMed CentralView ArticlePubMedGoogle Scholar
- Ross-Macdonald P, Coelho PS, Roemer T, Agarwal S, Kumar A, Jansen R, Cheung KH, Sheehan A, Symoniatis D, Umansky L, et al.: Large-scale analysis of the yeast genome by transposon tagging and gene disruption. Nature 1999, 402: 413–418. 10.1038/46558View ArticlePubMedGoogle Scholar
- Giaever G, Chu AM, Ni L, Connelly C, Riles L, Veronneau S, Dow S, Lucau-Danila A, Anderson K, Andre B, et al.: Functional profiling of the Saccharomyces cerevisiae genome. Nature 2002, 418: 387–391. 10.1038/nature00935View ArticlePubMedGoogle Scholar
- Tavazoie S, Hughes JD, Campbell MJ, Cho RJ, Church GM: Systematic determination of genetic network architecture. Nat Genet 1999, 22: 281–285. 10.1038/10343View ArticlePubMedGoogle Scholar
- Hubert L, Arabie P: Comparing Partitions. Journal of Classification 1985, 2: 193–218.View ArticleGoogle Scholar
- Ge H, Liu Z, Church GM, Vidal M: Correlation between transcriptome and interactome mapping data from Saccharomyces cerevisiae. Nat Genet 2001, 29: 482–486. 10.1038/ng776View ArticlePubMedGoogle Scholar
- Kemmeren P, van Berkum NL, Vilo J, Bijma T, Donders R, Brazma A, Holstege FC: Protein interaction verification and functional annotation by integrated analysis of genome-scale data. Mol Cell 2002, 9: 1133–1143.View ArticlePubMedGoogle Scholar
- Benjamini Y, Hochberg Y: Controlling the false discovery rate: a practical and powerful approach to multiple testing. Journal of the Royal Statistical Society, Series B 1995, 57: 289–300.Google Scholar
- Cho RJ, Huang M, Campbell MJ, Dong H, Steinmetz L, Sapinoso L, Hampton G, Elledge SJ, Davis RW, Lockhart DJ: Transcriptional regulation and function during the human cell cycle. Nat Genet 2001, 27: 48–54.PubMedGoogle Scholar
- Kim SK, Lund J, Kiraly M, Duke K, Jiang M, Stuart JM, Eizinger A, Wylie BN, Davidson GS: A gene expression map for Caenorhabditis elegans. Science 2001, 293: 2087–2092. 10.1126/science.1061603View ArticlePubMedGoogle Scholar
- Marcotte EM, Pellegrini M, Thompson MJ, Yeates TO, Eisenberg D: A combined algorithm for genome-wide prediction of protein function. Nature 1999, 402: 83–86. 10.1038/47048View ArticlePubMedGoogle Scholar
This article is published under license to BioMed Central Ltd. This is an Open Access article: verbatim copying and redistribution of this article are permitted in all media for any purpose, provided this notice is preserved along with the article's original URL.