ABC: software for interactive browsing of genomic multiple sequence alignment data
© Cooper et al; licensee BioMed Central Ltd. 2004
Received: 22 October 2004
Accepted: 08 December 2004
Published: 08 December 2004
Alignment and comparison of related genome sequences is a powerful method to identify regions likely to contain functional elements. Such analyses are data intensive, requiring the inclusion of genomic multiple sequence alignments, sequence annotations, and scores describing regional attributes of columns in the alignment. Visualization and browsing of results can be difficult, and there are currently limited software options for performing this task.
The Application for Browsing Constraints (ABC) is interactive Java software for intuitive and efficient exploration of multiple sequence alignments and data typically associated with alignments. It is used to move quickly from a summary view of the entire alignment via arbitrary levels of resolution to individual alignment columns. It allows for the simultaneous display of quantitative data, (e.g., sequence similarity or evolutionary rates) and annotation data (e.g. the locations of genes, repeats, and constrained elements). It can be used to facilitate basic comparative sequence tasks, such as export of data in plain-text formats, visualization of phylogenetic trees, and generation of alignment summary graphics.
The ABC is a lightweight, stand-alone, and flexible graphical user interface for browsing genomic multiple sequence alignments of specific loci, up to hundreds of kilobases or a few megabases in length. It is coded in Java for cross-platform use and the program and source code are freely available under the General Public License. Documentation and a sample data set are also available http://mendel.stanford.edu/sidowlab/downloads.html.
Functional elements in a genome accumulate inter-specific substitutions more slowly than neutral DNA throughout evolution . Therefore, comparing orthologous genomic sequences from related species is useful for the identification of elements that play important roles in the biology of an organism [2–7]. While the statistical and computational methods for extracting comparative information are variable, the types of data involved are generally quite similar. First, a multiple sequence alignment is necessary. Second, a vector of quantitative scores is produced that describes the similarity of the nucleotides observed in small windows, or individual columns, of the alignment; percent identity is the metric used by the popular program VISTA , while a variety of other scoring methods also exist [9–12]. Third, annotations are generated that highlight regions of the alignment that are under constraint or meet some other quantitative threshold. Fourth, annotations of features like transcripts, promoters, coding exons, and repeats provide functional context. Finally, the phylogenetic tree that relates the aligned sequences is important for both performing comparative analyses and for interpreting their results.
Simultaneous visualization of complex data such as these is of utmost importance both for experimentalists and for computational biologists. Several options currently exist for such visualization, but there are a variety of characteristics that distinguish the ABC from them. VISTA, for example, generates a static image that is not interactive . Other popular browsers such as phylo-VISTA  and PipMaker [14, 15] require the use of a particular alignment program and scoring scheme. Also, the ABC is not suited for genome-wide visualization. Other tools exist for this and are quite useful for the browsing of very large genomic intervals and major evolutionary events such as genomic rearrangements [16–19]. However, these programs are generally part of larger, more complex interfaces that are not necessarily ideal for targeted analysis of an individual alignment or genomic locus. Finally, we note that the ABC allows many annotation types and colors, is not web-based and can be used on a local machine as intensively as necessary, and the source code is open and freely available allowing users to modify and add features if desired.
Results and discussion
By default, the ABC displays graphical summaries of the quantitative information associated with the alignment. Scores are summarized regionally in consecutive non-overlapping windows. The size of these windows depends on the resolution, defined as the number of alignment columns summarized per pixel. The ABC has three distinct display modes, chosen automatically depending on the density of the information. At very low resolution, a histogram is displayed that plots the number of data points in each window that are at or below a specified value; note that regions containing many low scores will stand out as peaks in the histogram, as demonstrated by the clear association of peaks and the location of exons (Figure 1, top panel). At intermediate resolutions, a 'wiggly plot' is displayed, in which the average score for each regional window is plotted; in this case, regions containing many low scores will appear as valleys in the plot (Figure 1; middle panels). Finally, at very high resolution, the user may view the sequence data directly, along with the sequence names and a tree relating the sequences (Figure 1; lower panel).
A mobile and scalable zoom window allows for exploration of the summary views (Figure 1; upper panels). The user may drag and resize this rectangle, and when a desired region is selected a more detailed view can be obtained. This region will be expanded immediately below the parent display, with the resolution, score plot, and annotation adjusted accordingly (Figure 1; compare the start and stop coordinates of the black rectangle in the top panel to the start and stop coordinates of the entire panel immediately below). At all resolutions, annotation tracks are displayed immediately above the score/sequence display (Figure 1; all panels). An arbitrary number of tracks may be displayed, but the bottom two tracks are reserved for displaying information about transcripts with exons and introns. Colors for features can be specified individually using standard RGB notation. Other key features of the ABC include:
Mouse-over highlighting to reveal annotations, scores, coordinates, etc
Exporting of sequence, score, and annotation data
Searching sequence data for particular nucleotide strings
GoTo feature to quickly bring up a desired region
The ABC is flexible in that is has the ability to visualize diverse quantitative information and it has the capacity to display an arbitrary number of annotation types. It does not have a built-in scoring function; all data needs to be generated and formatted prior to being displayed in the ABC. While this may seem to be a drawback, it is in fact the intended function for an interface that has no preconceptions about the methodology that generated the data. Finally, the ABC is interactive, allowing the user to zoom in quickly from summary views of the comparative data to individual alignment columns. Zoom levels remain in the display, allowing the user to keep a birds-eye view of a large genomic region while focusing at much higher resolution on a small section within it.
The ABC is stand-alone alignment browsing software that is relatively easy to use and customize. While it was not designed as a genome-wide browser, it is well-suited for tasks associated with comparative sequence analysis: exploration of alignments of individual genomic loci; analyzing the relationship between known biological features and quantitative comparative data; visualizing results for researchers who develop and test methods for comparative sequence analysis; isolating sequence elements in a genomic locus for downstream applications like motif-discovery or primer design; generating graphics that characterize a multiple sequence alignment or region of an alignment; and potentially more applications that we have not yet considered. In our own research, for example, we have used it to display SNPs between different mouse strains in the context of a comparative alignment (not shown). This flexibility should be beneficial to researchers whose primary interest is comparative sequence analysis, but should also be valuable to those who use comparative analyses in support of other types of projects, such as experimental characterization of constrained elements. We also note that this flexibility distinguishes the ABC from other browsers that require built-in or specific types of score data and/or the use of a particular alignment program. The software is written in Java for cross-platform support, and the source code is freely available under the General Public License (GPL).
Availability and requirements
Project name: Application for Browsing Constraints (ABC)
Project home page: http://mendel.stanford.edu/sidowlab/downloads.html
Operating system: Platform independent
Programming language: Java
Other requirements: Java 1.4 or later
Application for Browsing Constraints
Graphical User Interface
General Public License
Single Nucleotide Polymorphism
University of California, Santa Cruz
GMC is a Howard Hughes Medical Institute Pre-doctoral fellow. AS acknowledges support from NIH/NHGRI. We thank two anonymous reviewers for their comments.
- Kimura M: The neutral theory of molecular evolution. Cambridge [Cambridgeshire] ; New York, Cambridge University Press; 1983:xv, 367.View ArticleGoogle Scholar
- Boffelli D, Nobrega MA, Rubin EM: Comparative genomics at the vertebrate extremes. Nat Rev Genet 2004, 5: 456–465. 10.1038/nrg1350View ArticlePubMedGoogle Scholar
- Cooper GM, Sidow A: Genomic regulatory regions: insights from comparative sequence analysis. Curr Opin Genet Dev 2003, 13: 604–610. 10.1016/j.gde.2003.10.001View ArticlePubMedGoogle Scholar
- Dermitzakis ET, Reymond A, Scamuffa N, Ucla C, Kirkness E, Rossier C, Antonarakis SE: Evolutionary discrimination of mammalian conserved non-genic sequences (CNGs). Science 2003, 302: 1033–1035. 10.1126/science.1087047View ArticlePubMedGoogle Scholar
- Göttgens B, Barton LM, Chapman MA, Sinclair AM, Knudsen B, Grafham D, Gilbert JG, Rogers J, Bentley DR, Green AR: Transcriptional regulation of the stem cell leukemia gene (SCL)--comparative analysis of five vertebrate SCL loci. Genome Res 2002, 12: 749–759. 10.1101/gr.45502PubMed CentralView ArticlePubMedGoogle Scholar
- Hardison RC: Comparative genomics. PLoS Biol 2003, 1: E58. 10.1371/journal.pbio.0000058PubMed CentralView ArticlePubMedGoogle Scholar
- Sidow A: Sequence first. Ask questions later. Cell 2002, 111: 13. 10.1016/S0092-8674(02)01003-6View ArticlePubMedGoogle Scholar
- Mayor C, Brudno M, Schwartz JR, Poliakov A, Rubin EM, Frazer KA, Pachter LS, Dubchak I: VISTA : visualizing global DNA sequence alignments of arbitrary length. Bioinformatics 2000, 16: 1046–1047. 10.1093/bioinformatics/16.11.1046View ArticlePubMedGoogle Scholar
- Sumiyama K, Kim CB, Ruddle FH: An efficient cis-element discovery method using multiple sequence comparisons based on evolutionary relationships. Genomics 2001, 71: 260–262. 10.1006/geno.2000.6422View ArticlePubMedGoogle Scholar
- Boffelli D, McAuliffe J, Ovcharenko D, Lewis KD, Ovcharenko I, Pachter L, Rubin EM: Phylogenetic shadowing of primate sequences to find functional regions of the human genome. Science 2003, 299: 1391–1394. 10.1126/science.1081331View ArticlePubMedGoogle Scholar
- Margulies EH, Blanchette M, Haussler D, Green ED: Identification and characterization of multi-species conserved sequences. Genome Res 2003, 13: 2507–2518. 10.1101/gr.1602203PubMed CentralView ArticlePubMedGoogle Scholar
- Cooper GM, Brudno M, Stone EA, Dubchak I, Batzoglou S, Sidow A: Characterization of evolutionary rates and constraints in three Mammalian genomes. Genome Res 2004, 14: 539–548. 10.1101/gr.2034704PubMed CentralView ArticlePubMedGoogle Scholar
- Shah N, Couronne O, Pennacchio LA, Brudno M, Batzoglou S, Bethel EW, Rubin EM, Hamann B, Dubchak I: Phylo-VISTA: interactive visualization of multiple DNA sequence alignments. Bioinformatics 2004, 20: 636–643. 10.1093/bioinformatics/btg459View ArticlePubMedGoogle Scholar
- Schwartz S, Elnitski L, Li M, Weirauch M, Riemer C, Smit A, Green ED, Hardison RC, Miller W: MultiPipMaker and supporting tools: Alignments and analysis of multiple genomic DNA sequences. Nucleic Acids Res 2003, 31: 3518–3524. 10.1093/nar/gkg579PubMed CentralView ArticlePubMedGoogle Scholar
- Elnitski L, Riemer C, Petrykowska H, Florea L, Schwartz S, Miller W, Hardison R: PipTools: a computational toolkit to annotate and analyze pairwise comparisons of genomic sequences. Genomics 2002, 80: 681–690. 10.1006/geno.2002.7018View ArticlePubMedGoogle Scholar
- Kalafus KJ, Jackson AR, Milosavljevic A: Pash: efficient genome-scale sequence anchoring by Positional Hashing. Genome Res 2004, 14: 672–678. 10.1101/gr.1963804PubMed CentralView ArticlePubMedGoogle Scholar
- Kent WJ, Sugnet CW, Furey TS, Roskin KM, Pringle TH, Zahler AM, Haussler D: The human genome browser at UCSC. Genome Res 2002, 12: 996–1006. 10.1101/gr.229102. Article published online before print in May 2002PubMed CentralView ArticlePubMedGoogle Scholar
- Kozik A, Kochetkova E, Michelmore R: GenomePixelizer--a visualization program for comparative genomics within and between species. Bioinformatics 2002, 18: 335–336. 10.1093/bioinformatics/18.2.335View ArticlePubMedGoogle Scholar
- Chakrabarti K, Pachter L: Visualization of multiple genome annotations and alignments with the K-BROWSER. Genome Res 2004, 14: 716–720. 10.1101/gr.1957004PubMed CentralView ArticlePubMedGoogle Scholar
- Cooper GM, Brudno M, Program NC, Green ED, Batzoglou S, Sidow A: Quantitative Estimates of Sequence Divergence for Comparative Analyses of Mammalian Genomes. Genome Res 2003, 13: 813–820. 10.1101/gr.1064503PubMed CentralView ArticlePubMedGoogle Scholar
- Smit AFA, Hubley R, Green P: RepeatMasker Open-3.0. <http://wwwrepeatmaskerorg> 1996.Google Scholar
- Brudno M, Do CB, Cooper GM, Kim MF, Davydov E, Green ED, Sidow A, Batzoglou S: LAGAN and Multi-LAGAN: efficient tools for large-scale multiple alignment of genomic DNA. Genome Res 2003, 13: 721–731. 10.1101/gr.926603PubMed CentralView 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.