BioSAVE: Display of scored annotation within a sequence context
© Pollock and Adryan; licensee BioMed Central Ltd. 2008
Received: 17 January 2008
Accepted: 20 March 2008
Published: 20 March 2008
Visualization of sequence annotation is a common feature in many bioinformatics tools. For many applications it is desirable to restrict the display of such annotation according to a score cutoff, as biological interpretation can be difficult in the presence of the entire data. Unfortunately, many visualisation solutions are somewhat static in the way they handle such score cutoffs.
We present BioSAVE, a sequence annotation viewer with on-the-fly selection of visualisation thresholds for each feature. BioSAVE is a versatile OS X program for visual display of scored features (annotation) within a sequence context. The program reads sequence and additional supplementary annotation data (e.g., position weight matrix matches, conservation scores, structural domains) from a variety of commonly used file formats and displays them graphically. Onscreen controls then allow for live customisation of these graphics, including on-the-fly selection of visualisation thresholds for each feature.
Possible applications of the program include display of transcription factor binding sites in a genomic context or the visualisation of structural domain assignments in protein sequences and many more. The dynamic visualisation of these annotations is useful, e.g., for the determination of cutoff values of predicted features to match experimental data. Program, source code and exemplary files are freely available at the BioSAVE homepage.
Visualization of sequence annotation is a common feature in many bioinformatics tools. For example, the widely used UCSC Genome Browser  and Ensembl  web sites, and also standalone programs such as Artemis  or the Integrated Genome Browser  graphically display genetic features in their genomic context. These tools also provide means to import additional user-defined numerical annotation, e.g., for the display of experimental data along the genome. These data are then displayed either as continuous line or bar graphs along the sequence, or as boxes covering a sequence range coloured according to the numerical data, i.e., a scored annotation.
A better understanding of a data set can often be achieved when its complexity can be reduced. While all of the aforementioned frameworks are suited to display user-defined data in its entirety, they do not allow for convenient and dynamic customisation of the graphics on-the-fly. The genome browsers are similar to other previously published tools for visualising annotation, such as gff2ps , in that they produce "static" output as either a Web page or a PostScript file. The Artemis and Integrated Genome Browser programs have their strength in presenting gene annotation or tiling array data, but are not easily customisable for display of other data types such as position weight matrix matches.
Here, we present BioSAVE (from Biological Sequence Annotation Visualisation), a graphical Mac OS X application for visualisation of scored annotation in a sequence context. It is written in Objective C and makes use of OS X-specific system libraries for responsive display of annotation changes. A plug-in interface allows to invoke external scripts (e.g., written in Perl or Python) upon loading an annotation file, which is useful for converting external formats into GFF2.
Results and Discussion
BioSAVE will read sequence data in FASTA format  or raw text, annotation scores from GFF2 files  or, using a plug-in interface, any other parse-able file format that contains coordinates and scores for an annotation. Only data with scores that lie within a user-specified score range will then be displayed along the sequence. Adjusting the score thresholds will dynamically change which annotations are visible, and at which colour level. Where other tools only support a fixed score range (e.g., 0..1000), this range automatically adjusts to the data present in the user's input. In addition to the dynamic view, BioSAVE can display additional tracks along the sequence comprising a conventional bar graph visualisation. The sequence view itself can be adjusted and allows for continuous zooming from the entire sequence down to the nucleotide level. In our hands, on a 1.67 GHz PowerPC G4 with 1 GB of RAM, BioSAVE can reliably handle up to 3 Mb of sequence along with about 3,000 microarray data points for numerical annotation.
The versatility of BioSAVE extends to the display of annotation scores also in protein sequences. Only a small number of tools for handling protein sequence allow to load and display additional numerical annotation (e.g., the Utopia Toolset ). However, these programs are comparable to the static display of Artemis and the Integrated Genome Browser, and do not allow for the on-the-fly customisation found in BioSAVE.
In the following section, we discuss two applications in which BioSAVE could be used in day-to-day research.
Usage example: Transcription factor binding site matches in a DNA sequence
The binding specificity of a transcription factor (TF) or other DNA-binding protein is often represented as position weight matrix (PWM), which allows for a numerical description of the binding site rather than a fixed search string. Many theoretical frameworks and bioinformatic tools exists to identify nucleotide sequences that are matched by a PWM, e.g., Patser , Clover  and Motifscanner . While these programs are powerful in the identification of PWM matches to a sequence, they all lack an immediate and intuitive display of the results, which is key to make useful biological interpretations.
The difficulty in the interpretation of PWM matches to a sequence is often the lack of a score threshold that reliably distinguishes biologically relevant matches from closely related hits in the sequence. Thus, it is desirable to display the identified PWM matches at various score cutoffs. This helps to establish a cutoff at which the PWM match is convincing, both in terms of its similarity with the sequence, and, if known, whether biologically relevant sequences are being recognised. For Patser and Clover, the Web interfaces RSA-tools  and MotifViz  are available to visualize PWM matches in input sequences, thus making the results accessible to the bench biologist. However, by default, these interfaces display every possible match to the PWM and lack the ability to dynamically adjust score thresholds.
Usage example: Structural hidden Markov model matches in a protein sequence
BioSAVE is a highly focused, easy-to-use tool for viewing scored annotation in a sequence context. It has many advantages over other sequence annotation visualisation systems not specifically designed for this task. For example, the aforementioned Web interfaces do not allow for prompt user interaction, and other tools developed primarily for visualisation of microarray data are too inflexible to handle scored annotation of different types or sources appropriately, e.g., quick highlighting of matches to a particular PWM etc. To our knowledge, BioSAVE is the only tool for annotation visualisation that handles DNA and protein sequences alike.
This paper describes several applications in which BioSAVE is used in our day-to-day research. Through its support of standard sequence and annotation formats, as well as the versatile plug-in system (employing any script or program that outputs GFF format), BioSAVE can be used in a wide range of applications not specifically discussed here. For example, it can be used for display of mutation hotspots in DNA or protein sequences (e.g., showing frequency) or display of predicted interaction surfaces along a protein sequence (e.g., showing a p-value).
Availability and requirements
The universal binary for Mac OS X is available for download at the BioSAVE web site , along with an extensive tutorial detailing many features not previously mentioned in this paper. The Objective C source code is available from the web site under the GNU GPL.
Project name: BioSAVE
Home page: http://www.mrc-lmb.cam.ac.uk/genomes/BioSAVE
Operating system: OS X
Programming language: Objective C
Other requirements: None
License: Freely available, source code under GNU GPL.
List of abbreviations
general feature format
position weight matrix
We thank Sarah Teichmann and Madan Babu for comments on the manuscript, and Charlie Massie and Ian Mills for testing the software extensively. Richard Pollock was supported by a Trinity College research grant to Sarah Teichmann.
- Kuhn RM, Karolchik D, Zweig AS, Trumbower H, Thomas DJ, Thakkapallayil A, Sugnet CW, Stanke M, Smith KE, Siepel A, Rosenbloom KR, Rhead B, Raney BJ, Pohl A, Pedersen JS, Hsu F, Hinrichs AS, Harte RA, Diekhans M, Clawson H, Bejerano G, Barber GP, Baertsch R, Haussler D, Kent WJ: The UCSC genome browser database: update 2007. Nucleic Acids Res 2007, 35: D668–673.PubMed CentralView ArticlePubMedGoogle Scholar
- Hubbard TJ, Aken BL, Beal K, Ballester B, Caccamo M, Chen Y, Clarke L, Coates G, Cunningham F, Cutts T, Down T, Dyer SC, Fitzgerald S, Fernandez-Banet J, Graf S, Haider S, Hammond M, Herrero J, Holland R, Howe K, Howe K, Johnson N, Kahari A, Keefe D, Kokocinski F, Kulesha E, Lawson D, Longden I, Melsopp C, Megy K, Meidl P, Ouverdin B, Parker A, Prlic A, Rice S, Rios D, Schuster M, Sealy I, Severin J, Slater G, Smedley D, Spudich G, Trevanion S, Vilella A, Vogel J, White S, Wood M, Cox T, Curwen V, Durbin R, Fernandez-Suarez XM, Flicek P, Kasprzyk A, Proctor G, Searle S, Smith J, Ureta-Vidal A, Birney E: Ensembl 2007. Nucleic Acids Res 2007, 35: D610–617.PubMed CentralView ArticlePubMedGoogle Scholar
- Rutherford K, Parkhill J, Crook J, Horsnell T, Rice P, Rajandream MA, Barrell B: Artemis: sequence visualization and annotation. Bioinformatics 2000, 16: 944–945.View ArticlePubMedGoogle Scholar
- Abril JF, Guigó R: gff2ps: visualizing genomic annotations. Bioinformatics 2000, 16: 743–744.View ArticlePubMedGoogle Scholar
- McDermott P, Sinnott J, Thorne D, Pettifer S, Attwood T: An Architecture for Visualisation and Interactive Analysis of Proteins. Proceedings of the Fourth International Conference on Coordinated and Multiple Views in Exploratory Visualization, Washington, DC, USA 2006, 55–65.View ArticleGoogle Scholar
- Hertz GZ, Stormo GD: Identifying DNA and protein patterns with statistically significant alignments of multiple sequences. Bioinformatics 1999, 15: 563–577.View ArticlePubMedGoogle Scholar
- Frith MC, Fu Y, Yu L, Chen JF, Hansen U, Weng Z: Detection of functional DNA motifs via statistical over-representation. Nucleic Acids Res 2004, 32: 1372–1381.PubMed CentralView ArticlePubMedGoogle Scholar
- Down TA, Hubbard TJ: NestedMICA: sensitive inference of over-represented motifs in nucleic acid sequence. Nucleic Acids Res 2005, 33: 1445–1453.PubMed CentralView ArticlePubMedGoogle Scholar
- van Helden J: Regulatory sequence analysis tools. Nucleic Acids Res 2003, 31: 3593–3596.PubMed CentralView ArticlePubMedGoogle Scholar
- Fu Y, Frith MC, Haverty PM, Weng Z: MotifViz: an analysis and visualization tool for motif discovery. Nucleic Acids Res 2004, 32: W420–423.PubMed CentralView ArticlePubMedGoogle Scholar
- Massie CE, Adryan B, Barbosa-Morais NL, Lynch AG, Tran MG, Neal DE, Mills IG: New androgen receptor genomic targets show an interaction with the ETS1 transcription factor. EMBO Rep 2007, 8: 871–878.PubMed CentralView ArticlePubMedGoogle Scholar
- Adryan B, Woerfel G, Birch-Machin I, Gao S, Quick M, Meadows L, Russell S, White R: Genomic mapping of Suppressor of Hairy-wing binding sites in Drosophila. Genome Biol 2007, 8(8):R167.PubMed CentralView ArticlePubMedGoogle Scholar
- Siepel A, Bejerano G, Pedersen JS, Hinrichs AS, Hou M, Rosenbloom K, Clawson H, Spieth J, Hillier LW, Richards S, Weinstock GM, Wilson RK, Gibbs RA, Kent WJ, Miller W, Haussler D: Evolutionarily conserved elements in vertebrate, insect, worm, and yeast genomes. Genome Res 2005, 15: 1034–1050.PubMed CentralView ArticlePubMedGoogle Scholar
- Wilson D, Madera M, Vogel C, Chothia C, Gough J: The SUPERFAMILY database in 2007: families and functions. Nucleic Acids Res 2007, 35(Database issue):D308-D313.PubMed CentralView ArticlePubMedGoogle Scholar
- Finn RD, Mistry J, Schuster-Bockler B, Griffiths-Jones S, Hollich V, Lassmann T, Moxon S, Marshall M, Khanna A, Durbin R, Eddy SR, Sonnhammer EL, Bateman A: Pfam: clans, web tools and services. Nucleic Acids Res 2006, 34(Database issue):D247-D251.PubMed CentralView ArticlePubMedGoogle Scholar
- Blom N, Gammeltoft S, Brunak S: Sequence and structure-based prediction of eukaryotic protein phosphorylation sites. J Mol Biol 1999, 294: 1351–1362.View ArticlePubMedGoogle Scholar
This article is published under license to BioMed Central Ltd. This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/2.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.