Open Access

FunSpec: a web-based cluster interpreter for yeast

  • Mark D Robinson1,
  • Jörg Grigull1,
  • Naveed Mohammad1 and
  • Timothy R Hughes1, 2Email author
BMC Bioinformatics20023:35

DOI: 10.1186/1471-2105-3-35

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


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 [13]. 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 [6] and the integration of a wide variety of data types including Pfam domains, SWISS-PROT keywords and known pathways (e.g. DRAGON View; see [7]). 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 [2]. Similarly, the Gene Ontology (GO) Consortium organizes the same information under the categorizations: Biological Process, Molecular Function and Cellular Component [10]. Other online knowledge bases encompass computationally-derived gene properties, such as domain and motif information (e.g. SMART [11], Pfam [12]). In addition, there are now mass stores of high-throughput biological data accumulating from experimental results such as synthetic lethal [13], protein interactions by yeast two hybrid [14, 15] or mass spectrometry [16, 17], immunolocalization [18], large-scale phenotypic analysis [19], and microarray analysis (eg., [8]).

A listing of the current databases that can be evaluated with FunSpec is given in Table 1. The published experimental information that can be queried is given in Table 2. Through our own analyses of recent literature on ribosomal RNA processing, additional categories are provided to elicit enrichment of pre-ribosome components. Also, a separate category has been added to encompass genes that are essential for viability [20].
Table 1

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


Table 2

FunSpec-searchable published datasets

Data Type


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


Proteome Localization




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 [1] 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 [21]. 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.


In this section, we describe instances where FunSpec may be useful. In Figure 1, a 2-dimensional clustering of gene expression data from [8], covering 237 genes and 271 experiments is presented. The 76 genes indicated by the red bar were submitted to MIPS Functional Classification database on the FunSpec website. The categories "amino acid biosynthesis" and "amino acid metabolism" were deemed most unlikely to occur by chance (p < 10-14 in both cases). To show the specificity of evaluating clusters in this manner, the 6 genes indicated in blue (isolated by visual inspection) have the same overall expression profile as genes generally involved in amino acid metabolism, but with a few distinct differences. FunSpec identifies these as involved in "methionine metabolism." (p = 1.0 × 10-7).
Figure 1

A 2-D clustering analysis of yeast gene expression data from [8].

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 [22]) (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 [18], 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 [21]), 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 [25]).

There are other methods of determining significance (e.g. binomial distribution [26]; representation factor [27]) and other means of relating clusters to functional annotation (e.g. associated terms and phrases in Medline abstracts [6]) 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.

Authors’ Affiliations

Banting and Best Department of Medical Research, University of Toronto
Department of Medical Genetics and Microbiology, University of Toronto


  1. 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
  2. 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
  3. 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
  4. 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
  5. 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
  6. 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
  7. 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
  8. 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
  9. 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
  10. The Gene Ontology Consortium: Creating the gene ontology resource: design and implementation. Genome Res 2001, 11: 1425–1433. 10.1101/gr.180801View ArticleGoogle Scholar
  11. 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
  12. 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
  13. 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
  14. 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
  15. 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
  16. 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
  17. 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
  18. 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
  19. 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
  20. 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
  21. 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
  22. Hubert L, Arabie P: Comparing Partitions. Journal of Classification 1985, 2: 193–218.View ArticleGoogle Scholar
  23. 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
  24. 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
  25. 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
  26. 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
  27. 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
  28. 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


© Robinson et al; licensee BioMed Central Ltd. 2002

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.