Bioinformatics Advance Access originally published online on August 4, 2005
Bioinformatics 2005 21(18):3674-3676; doi:10.1093/bioinformatics/bti610
| ||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
Blast2GO: a universal tool for annotation, visualization and analysis in functional genomics research
1Centro de Genómica, Instituto Valenciano de Investigaciones Agrarias Moncada, Valencia, Spain
2BET-ITACA, Universidad Politécnica de Valencia Valencia, Spain
*To whom correspondence should be addressed.
 |
Abstract |
|---|
 TOP Abstract INTRODUCTIONOBTAINING GO TERMSANNOTATION ASSIGNMENTSTATISTICSVISUALIZATIONVALIDATIONCONCLUSIONSREFERENCES |
|---|
Summary: We present here Blast2GO (B2G), a research tool designed with the main purpose of enabling Gene Ontology (GO) based data mining on sequence data for which no GO annotation is yet available. B2G joints in one application GO annotation based on similarity searches with statistical analysis and highlighted visualization on directed acyclic graphs. This tool offers a suitable platform for functional genomics research in non-model species. B2G is an intuitive and interactive desktop application that allows monitoring and comprehension of the whole annotation and analysis process.
Availability: Blast2GO is freely available via Java Web Start at http://www.blast2go.de
Supplementary material: http://www.blast2go.de -> Evaluation
Contact: aconesa{at}ivia.es; stefang{at}fis.upv.es
|
INTRODUCTION |
|---|
|
TOPAbstractINTRODUCTIONOBTAINING GO TERMSANNOTATION ASSIGNMENTSTATISTICSVISUALIZATIONVALIDATIONCONCLUSIONSREFERENCES |
|---|
One of the most important aspects in mining genomics data is to associate individual sequences and related expression information with biological function. Automatic functional annotation is an effective approach to solve this problem. Functional annotation allows categorization of genes in functional classes, which can be very useful to understand the physiological meaning of large amounts of genes and to assess functional differences between subgroups of sequences. The Gene Ontology (GO) developed at the GO Consortium (Ashburner et al., 2000) provides a suitable framework for this kind of analysis, due to the wide scope of biology covered and its directed acyclic graph (DAG) structure that enables visualization in the context of biological dependences. Different development teams have released software to analyze sequences by the use of GO. A variety of desktop and web applications are available to electronically assign GO terms to unknown sequences based on similarity (Martin et al., 2004; Groth et al., 2004; Khan et al., 2003; Zehetner, 2003) or to analyze genomic data in the context of gene annotation (Al-Shahrour et al., 2004; Doniger et al., 2003). However, when trying to perform GO-based analysis in poorly characterized organisms we encountered a number of drawbacks. In general, these tools are either not designed for high-throughput sequence annotation, are limited in their mining and visualization capabilities, or accept only gene or probe identifiers as input data, making them restrictive to annotated sequences already deposited in public databases. In order to provide a suitable solution to these limitations we have developed Blast2GO (B2G), a universal GO annotation, visualization and statistics framework that brings advanced functional analysis to the genomics research of non-model species. B2G has been design to (1) allow automatic and high throughput sequence annotation and (2) integrate functionality for annotation-based data mining. Briefly, B2G uses BLAST (Altschul et al., 1990) to find homologs to fasta formatted input sequences. The program extracts GO terms to each obtained hit by mapping to existent annotation associations. An annotation rule finally assigns GO terms to the query sequence. Annotation and functional analysis can be visualized in a graph form reconstructing the GO relationships and color-highlighting the most relevant areas (Fig. 1). B2G was conceived to be an attractive tool for research environments where genetic and/or computational resources are limited and where much work is still done in an explorative fashion. B2G is a user-friendly, easy to distribute and low maintenance tool. It allows monitoring and interaction at different steps of the analysis, and emphasizes visualization as an important component of knowledge acquisition. B2G is a Java application made available by Java Web Start. It is platform independent and has no further requirements than an Internet connection.
|
|
OBTAINING GO TERMS |
|---|
|
TOPAbstractINTRODUCTIONOBTAINING GO TERMSANNOTATION ASSIGNMENTSTATISTICSVISUALIZATIONVALIDATIONCONCLUSIONSREFERENCES |
|---|
The first step in B2G is to find sequences similar to a query set by Blast searching. Homology search can either be done at public databases (e.g. NCBI nr and est using QBlast) or custom databases (e.g. GO annotated sequence sets and single species DBs) when a local www-Blast installation is available. Blast expectation values (E-value) and hit number thresholds are provided to retrieve significant results. To avoid the danger of annotation by short matches with low E-values, an additional filter can be set to the minimal alignment length (hsp-length). Annotation, however, will ultimately be based on sequence similarity levels since similarity percentages are independent on database size and more intuitive than E-values.
In order to retrieve GO terms associated with the obtained hits, a quite straight forward mapping is made. By using Blast hit gene identifiers (gi) and gene accessions B2G retrieves all GO annotations for the hit sequences, together with their evidence codes (EC). ECs can be interpreted as an index of the trustworthiness of the GO annotation. At the end of the mapping processes, for each query sequence, a set of candidate annotations from different hits of diverse similarity levels and various annotation sources is gathered.
|
ANNOTATION ASSIGNMENT |
|---|
|
TOPAbstractINTRODUCTIONOBTAINING GO TERMSANNOTATION ASSIGNMENTSTATISTICSVISUALIZATIONVALIDATIONCONCLUSIONSREFERENCES |
|---|
Annotation is performed by applying an annotation rule (AR) to the obtained ontologies. The rule seeks to find the most specific annotations with a certain level of reliability. This process is adjustable in specificity and stringency.
For each candidate GO an annotation score (AS) is computed. The AS is composed of two additive terms. The first, direct term (DT), represents the highest hit similarity of this GO weighted by a factor corresponding to its EC. By employing ECs, B2G promotes the assignment of annotations with experimental evidence and penalizes electronic annotations or low traceability. The EC weights have been taken following recommendations of the GO Consortium and can be modified if desired. The second term (AT) of the AS provides the possibility of abstraction. This is defined as annotation to a parent node when several child nodes are present in the GO candidate collection. This term multiplies the number of total GOs unified at the node by a user defined GO weight factor that controls the possibility and strength of abstraction. Finally, the AR selects the lowest term per branch that lies over a user defined threshold. In an analytical form, DT, AT and the AR terms are defined as follows:
 |
|
STATISTICS |
|---|
|
TOPAbstractINTRODUCTIONOBTAINING GO TERMSANNOTATION ASSIGNMENTSTATISTICSVISUALIZATIONVALIDATIONCONCLUSIONSREFERENCES |
|---|
Once GO annotation is available through B2G (uploading an existing annotation file is also supported), the application offers the possibility of direct statistical analysis on gene function information. A common analysis is the statistical assessment of GO term enrichments in a group of interesting genes when compared with a reference group. This functionality was introduced in B2G by integrating Gossip (Blüthgen et al., 2004). Gossip computes Fisher's Exact Test applying robust FDR (false discovery rate) correction for multiple testing and returns a list of significant GO terms ranked by their corrected or one-test P-values. Furthermore B2G offers various statistical charts summarizing the results obtained at blasting, mapping or annotation. Bar or pie charts of similarity/E-value distributions, EC distributions and annotation statistics (GOs/Seqs) can be generated, saved and printed.
|
VISUALIZATION |
|---|
|
TOPAbstractINTRODUCTIONOBTAINING GO TERMSANNOTATION ASSIGNMENTSTATISTICSVISUALIZATIONVALIDATIONCONCLUSIONSREFERENCES |
|---|
Visualization is an important aspect in B2G. For each sequence, the progress in the annotation process and the final annotation step are visualized on the main application table by successive color changes. This allows the researcher to readily spot sequences that failed the initial annotation process and, if desired, modify annotation parameters for those. Furthermore, the joined biological meaning of a set of sequences can be visualized on the GO DAG by color-intensity highlighting of the most relevant nodes in a combined sequence graph. Those nodes are identified by computing a node score that takes into account the number of sequences converging at one node and penalizes by the distance to the node where each sequence was annotated. Alternatively, when an enrichment analysis is available, graph color highlighting by statistical results will show the GO-term specificity of the query subset.
|
VALIDATION |
|---|
|
TOPAbstractINTRODUCTIONOBTAINING GO TERMSANNOTATION ASSIGNMENTSTATISTICSVISUALIZATIONVALIDATIONCONCLUSIONSREFERENCES |
|---|
The performance of Blast2GO has been tested using a dataset for which annotation and functional information was available. The methodology and results of this evaluation are given as supplementary material and are available at the B2G site. Our results show that Blast2GO reaches an annotation accuracy of 65–70%, which is commonly reported in automatic GO annotation methods (Martin et al., 2004; Khan et al., 2003). More interestingly, this evaluation shows that the tool is successful in extracting relevant functional features of these sequences based on the use of the predicted annotation.
|
CONCLUSIONS |
|---|
|
TOPAbstractINTRODUCTIONOBTAINING GO TERMSANNOTATION ASSIGNMENTSTATISTICSVISUALIZATIONVALIDATIONCONCLUSIONSREFERENCES |
|---|
By joining annotation to function analysis B2G provides a powerful data mining tool ideally suited to support genomic research in non-model species. Its species-independent character and different data input fronts makes it a valuable mining resource for potentially any organism. B2G combines high-throughput analysis, statistical evaluation and biology framed visualization with a high degree of user interaction. Further developments of Blast2GO will include extension to multiple annotation types and novel statistical analysis tools.
|
Acknowledgments |
|---|
The authors thank Dr Timothy Williams for fruitful discussions and comments on the software and Nils Bluethgen for kindly providing the Gossip software and supporting integration in B2G. This work has been funded by MCyT (GEN 2001 - 4885-C05-03) and eTumour Project (FP6-2002-LIFESCIHEALTH 503094). The authors thank the INBIOMED G03/160 research thematic network financed by FIS of the Instituto de Salud Carlos III.
Conflict of Interest: none declared.
|
Footnotes |
|---|
The authors wish it to be known that, in their opinion, the first two authors should be regarded as joint First Authors. 
Received on June 27, 2005; revised on July 28, 2005; accepted on July 29, 2005
|
REFERENCES |
|---|
|
TOPAbstractINTRODUCTIONOBTAINING GO TERMSANNOTATION ASSIGNMENTSTATISTICSVISUALIZATIONVALIDATIONCONCLUSIONSREFERENCES |
|---|
Al-Shahrour, F., et al. (2004) FatiGO: a web tool for finding significant associations of Gene Ontology terms with groups of genes. Bioinformatics, 20, 578–580
Altschul, S.F., et al. (1990) Basic local alignment search tool. J. Mol. Biol., 215, 403–410[CrossRef][ISI][Medline].
Ashburner, M., et al. (2000) Gene Ontology: tool for the unification of biology. Nat. Genet., 25, 25–29[CrossRef][ISI][Medline].
Blüthgen, N., Brand, K., Cajavec, B., Swat, M., Herzel, H., Beule, D. (2004) Biological Profiling of Gene Groups utilizing Gene Ontology – A Statistical Framework. arXiv:q-bio.GN/0407034, 1, 1.
Doniger, S., et al. (2003) MAPPFinder: using Gene Ontology and GenMAPP to create a global gene-expression profile from microarray data. Genome Biol., 4, R7[CrossRef][Medline].
Groth, D., et al. (2004) GOblet: a platform for Gene Ontology annotation of anonymous sequence data. Nucleic Acids Res., 32, 313–317.
Khan, S., et al. (2003) GoFigure: automated Gene OntologyTM annotation. Bioinformatics, 19, 2484–2485
Martin, D., et al. (2004) GOtcha: a new method for prediction of protein function assessed by the annotation of seven genomes. BMC Bioinformatics, 5, 178[CrossRef][Medline].
Zehetner, G. (2003) OntoBlast function: from sequence similarities directly to potential functional annotations by ontology terms. Nucleic Acids Res., 31, 3799–3803
CiteULike 
Connotea 
Del.icio.us What's this?
This article has been cited by other articles:
|
|
 |
|
  F. Al-Shahrour, J. Carbonell, P. Minguez, S. Goetz, A. Conesa, J. Tarraga, I. Medina, E. Alloza, D. Montaner, and J. Dopazo Babelomics: advanced functional profiling of transcriptomics, proteomics and genomics experiments Nucleic Acids Res., July 1, 2008; 36(suppl_2): W341 - W346. [Abstract] [Full Text] [PDF]
|
|
|
|
|
|
S. Gotz, J. M. Garcia-Gomez, J. Terol, T. D. Williams, S. H. Nagaraj, M. J. Nueda, M. Robles, M. Talon, J. Dopazo, and A. Conesa High-throughput functional annotation and data mining with the Blast2GO suite Nucleic Acids Res., June 1, 2008; 36(10): 3420 - 3435. [Abstract] [Full Text] [PDF]
|
|
|
|
 |
|
 C. Martens, K. Vandepoele, and Y. Van de Peer Whole-genome analysis reveals molecular innovations and evolutionary transitions in chromalveolate species PNAS, March 4, 2008; 105(9): 3427 - 3432. [Abstract] [Full Text] [PDF]
|
|
|
|
 |
|
 S. A. M. Martin, J. B. Taggart, P. Seear, J. E. Bron, R. Talbot, A. J. Teale, G. E. Sweeney, B. Hoyheim, D. F. Houlihan, D. R. Tocher, et al. Interferon type I and type II responses in an Atlantic salmon (Salmo salar) SHK-1 cell line by the salmon TRAITS/SGP microarray Physiol Genomics, December 19, 2007; 32(1): 33 - 44. [Abstract] [Full Text] [PDF]
|
|
|
|
 |
|
 J. K. Hane, R. G.T. Lowe, P. S. Solomon, K.-C. Tan, C. L. Schoch, J. W. Spatafora, P. W. Crous, C. Kodira, B. W. Birren, J. E. Galagan, et al. Dothideomycete Plant Interactions Illuminated by Genome Sequencing and EST Analysis of the Wheat Pathogen Stagonospora nodorum PLANT CELL, November 1, 2007; 19(11): 3347 - 3368. [Abstract] [Full Text] [PDF]
|
|
|
|
 |
|
 K. A. Duthie, L. C. Osborne, L. J. Foster, and N. Abraham Proteomics Analysis of Interleukin (IL)-7-induced Signaling Effectors Shows Selective Changes in IL-7R{alpha}449F Knock-in T Cell Progenitors Mol. Cell. Proteomics, October 1, 2007; 6(10): 1700 - 1710. [Abstract] [Full Text] [PDF]
|
|
|
|
 |
|
 J. D. Hackett, H. S. Yoon, S. Li, A. Reyes-Prieto, S. E. Rummele, and D. Bhattacharya Phylogenomic Analysis Supports the Monophyly of Cryptophytes and Haptophytes and the Association of Rhizaria with Chromalveolates Mol. Biol. Evol., August 1, 2007; 24(8): 1702 - 1713. [Abstract] [Full Text] [PDF]
|
|
|
|
 |
|
 M. J. Nueda, A. Conesa, J. A. Westerhuis, H. C. J. Hoefsloot, A. K. Smilde, M. Talon, and A. Ferrer Discovering gene expression patterns in time course microarray experiments by ANOVA SCA Bioinformatics, July 15, 2007; 23(14): 1792 - 1800. [Abstract] [Full Text] [PDF]
|
|
|
|
|
|
A. Labarga, F. Valentin, M. Anderson, and R. Lopez Web Services at the European Bioinformatics Institute Nucleic Acids Res., July 13, 2007; 35(suppl_2): W6 - W11. [Abstract] [Full Text] [PDF]
|
|
|
|
|
|
S. H. Nagaraj, N. Deshpande, R. B. Gasser, and S. Ranganathan ESTExplorer: an expressed sequence tag (EST) assembly and annotation platform Nucleic Acids Res., July 13, 2007; 35(suppl_2): W143 - W147. [Abstract] [Full Text] [PDF]
|
|
|
|
|
|
G. Lopez, A. Valencia, and M. L. Tress firestar--prediction of functionally important residues using structural templates and alignment reliability Nucleic Acids Res., July 13, 2007; 35(suppl_2): W573 - W577. [Abstract] [Full Text] [PDF]
|
|
|
|
|
|
J. Adachi, C. Kumar, Y. Zhang, and M. Mann In-depth Analysis of the Adipocyte Proteome by Mass Spectrometry and Bioinformatics Mol. Cell. Proteomics, July 1, 2007; 6(7): 1257 - 1273. [Abstract] [Full Text] [PDF]
|
|
|
|
|
|
B. Macek, I. Mijakovic, J. V. Olsen, F. Gnad, C. Kumar, P. R. Jensen, and M. Mann The Serine/Threonine/Tyrosine Phosphoproteome of the Model Bacterium Bacillus subtilis Mol. Cell. Proteomics, April 1, 2007; 6(4): 697 - 707. [Abstract] [Full Text] [PDF]
|
|
|
|
 |
|
 S. H. Nagaraj, R. B. Gasser, and S. Ranganathan A hitchhiker's guide to expressed sequence tag (EST) analysis Brief Bioinform, January 1, 2007; 8(1): 6 - 21. [Abstract] [Full Text] [PDF]
|
|
|
|
 |
|
 D. S. Durica, D. Kupfer, F. Najar, H. Lai, Y. Tang, K. Griffin, P. M. Hopkins, and B. Roe EST library sequencing of genes expressed during early limb regeneration in the fiddler crab and transcriptional responses to ecdysteroid exposure in limb bud explants Integr. Comp. Biol., December 1, 2006; 46(6): 948 - 964. [Abstract] [Full Text] [PDF]
|
|
|
|
|
|
Q. W. T. Chan, C. G. Howes, and L. J. Foster Quantitative Comparison of Caste Differences in Honeybee Hemolymph Mol. Cell. Proteomics, December 1, 2006; 5(12): 2252 - 2262. [Abstract] [Full Text] [PDF]
|
|
| ||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||