- Open Access
Genome Environment Browser (GEB): a dynamic browser for visualising high-throughput experimental data in the context of genome features
BMC Bioinformatics volume 9, Article number: 501 (2008)
There is accumulating evidence that the milieu of repeat elements and other non-genic sequence features at a given chromosomal locus, here defined as the genome environment, can play an important role in regulating chromosomal processes such as transcription, replication and recombination. The availability of whole-genome sequences has allowed us to annotate the genome environment of any locus in detail. The development of genome wide experimental analyses of gene expression, chromatin modification and chromatin proteins means that it is now possible to identify potential links between chromosomal processes and the underlying genome environment. There is a need for novel bioinformatic tools that facilitate these studies.
We developed the Genome Environment Browser (GEB) in order to visualise the integration of experimental data from large scale high throughput analyses with repeat sequence features that define the local genome environment. The browser has incorporated dynamic scales adjustable in real-time, which enables scanning of large regions of the genome as well as detailed investigation of local regions on the same page without the need to load new pages. The interface also accommodates a 2-dimensional display of repetitive features which vary substantially in size, such as LINE-1 repeats. Specific queries for preliminary quantitative analysis of genome features can also be formulated, results of which can be exported for further analysis.
The Genome Environment Browser is a versatile program which can be easily adapted for displaying all types of genome data with known genomic coordinates. It is currently available at http://web.bioinformatics.ic.ac.uk/geb/.
Common repetitive DNA elements, which include satellite DNA, long interspersed repeats (LINE), short interspersed repeat (SINE) and long terminal repeat (LTR) elements, comprise 37% of the rodent and 42% of the human genome sequence respectively [1, 2]. By comparison, exons of genes comprise only approximately 2% of sequence. These common repeat elements, together with other features such as CpG islands , scaffold-attachment regions (SARs) , and transcription factor binding sites, shape the genome environment in which a gene resides. There is accumulating evidence that the genome environment can be important for the regulation of gene expression. For example, SARs play a role in regulating MHC Class I gene expression in humans , LTR retrotransposons influence developmentally regulated expression of genes in mouse oocytes and preimplantation embryos , and LINE-1 (L1) elements modulate transcription of human genes .
With the DNA sequence data generated from genome projects, we can now paint a fuller picture of a gene's environment in silico. Added to this, the development of high throughput DNA sequence-based experimental strategies such as whole-genome gene expression microarrays and ChIP-on-chip/ChIP sequencing means that it is now possible to look for correlations between underlying sequence features, the transcriptome, and epigenetic features such as DNA methylation, covalent histone modification and chromatin protein distribution. Importantly, novel bioinformatics and software tools are needed, both to analyse the large datasets generated by such studies and to facilitate elucidation of previously unappreciated relationships between underlying sequence features, gene expression and epigenetic modification. Here we describe development of the Genome Environment Browser, a novel tool to aid visualisation and analysis of genome wide data in the context of underlying genomic features.
GEB is designed as a set of software components that automatically build a core database of genomic feature data from the Ensembl database for any available species, using the Ensembl Perl API, with the features to be retrieved defined in a configuration file. The settings for the local storage database and Ensembl connection are also stored in the configuration file so once initialized the software automatically builds the GEB data without the requirement for further user input. For repeat features, such as LINEs, individual classes of the repeat can be defined to be stored separately to view as an individual track in the GEB viewer. We have used this feature for the display of LINE L1 elements. The data is stored in a standard relational database, specifically MySQL . Alternatively we provide pre-built databases of the latest Ensembl builds for human, mouse and rat on our web site. These can be used as the basis of a core GEB installation to which users' own data can be added. Further scripts are provided for the storage of non-Ensembl features and microarray data, both expression and ChIP-chip. These scripts require the data to be in a tab delimited format, which can be created for example by parsing genomic annotation software output or from an Excel spreadsheet for microarray data. We have used this feature for the LINE L1 components (UTRs and ORFs) and CpG island predictions within our custom annotations. We found the CpG island Ensembl predictions to be conservative so for our predictions we chose to use the EMBOSS newcpgreport program , the output of which was parsed to produce a tab delimited file as required. To facilitate the ease of adding data to GEB, including the core database, a Java graphical user interface is provided to manage the uploading of the database itself and addition/deletion of other feature data, including microarray data. This program can be used in place of the Perl scripts but it is platform-dependent and is only suitable for Windows users. See additional file 1 for the full GEB installation guide.
A separate, platform-independent, standalone Java front-end application provides a graphical interface for navigating the database and visualisation of the genome features. The Java application was created using the GenoViz SDK , a Java toolkit providing custom classes for the visualisation of genomic annotation data.
Results and discussion
Genomic feature annotation and browsing
We have been investigating a proposed role for the L1 retrotransposon in facilitating X-chromosome inactivation, the process whereby a single X chromosome is epigenetically inactivated in cells of female mammals . L1s (a major family of LINEs) are non-LTR retrotransposons 6–7 kb in length that collectively account for approximately a fifth of the mammalian genome . To facilitate this study we have developed a novel software tool, the Genome Environment Browser (GEB), for visualising epigenetic and transcriptomic features in the context of the genome environment.
Currently GEB is configured to display Human and Mouse genomes. While many public genome browsers are designed to be gene-oriented, GEB focuses on both genes and non-genic elements in their vicinity. In the default version of the browser, CpG islands, LINEs (including families LINE-1, LINE-2 and LINE-3), L1s, SINEs and LTRs are annotated in detail in the context of known Ensembl genes and non-coding genes (mainly functional RNAs). Because of our interest in L1 elements we have included annotation of the functional components, the 5' UTR, ORF1, ORF2 and 3' UTR (Figure 1A), and additionally full length (FL) L1s. The distribution and density of all these features can be visualised at two levels: first a histogram display for a panoramic, whole-chromosome view (Figure 1B), and second, a physical map display for more detailed analysis of regions between 500 bp and 25 Mb (Figure 1C). See additional file 2 for the GEB user guide.
The histogram display is composed of a parallel set of histograms, each plotting the copy number/density of a genome feature across the entire chromosome (typically per Mb). It is powerful in capturing the modularity of genome features across a chromosome. For example, the previously reported positive correlation between CpGI and genes  and the reciprocal relationship between L1 and SINEs [14, 15] are reproduced graphically (Figure 1B). Specific regions with interesting patterns, e.g. a gene-rich region adjacent to a gene desert, can be easily identified and precisely selected on the histogram display to be examined in greater detail in the physical map display.
In the physical map display, unlike some public genome browsers, all genome features (including interspersed repeats) are displayed on their respective strands. A very important feature of GEB is the capability to display genomic features in two-dimensional multi-colour graphics reminiscent of a dot-plot, as illustrated by our custom annotation of L1 elements. The length of homology of an L1 element to each functional component in the FL-L1 consensus sequence is plotted on the vertical axis, whereas the physical length of the L1 element is indicated on the horizontal axis as in public genome browsers (Figure 1C). Relatively rare FL-L1s appear as a long, continuous diagonal arrangement of UTRs and ORFs, and are clearly distinguished from the bulk of truncated L1s (mostly remnants of the 3' UTR), shown as shortened lines or bars.
Another key feature of GEB is the rapid transition between high and low resolutions in the physical map display. The "condensed" view of the selected region is initially displayed at 10–200 kb resolution (depending on the size of the region), and can be zoomed in on the same screen flexibly to a detailed 2 bp resolution (and zoomed back out again), abolishing the need to refresh or load a new page. This allows uninterrupted visual scanning of genome features and pattern searching across very large genomic distances, up to 25 Mb. In essence, the physical map display provides a zoomed-in view of the spatial relationship between features across the region of interest. Local enrichment or depletion of any type of feature can be readily visualised.
Data retrieval and preliminary quantitative data analysis tools
Detailed textual annotation information on all features can be explored at a mouse click in physical map display. For the retrieval of other types of annotation data (e.g. gene ontology, protein-related annotation, homology to genes in other species), we have incorporated a functionality in GEB which brings up relevant Ensembl pages (ContigView and GeneView respectively) in the user's default web browser, depending on whether a genomic region or a gene has been selected in GEB histogram/physical map displays. Additionally, a simple quantitative sequence analysis tool is available, which calculates the copy number and/or percentage sequence representation of any type of annotated feature across a genome region of any size. For a more detailed analysis, the calculation can also be done in defined windows (e.g. every 100 kb) across a region (Figure 2).
Displaying custom experimental data in GEB
The GEB histogram and physical map display framework can be applied to one or more sets of custom genome/chromosome-wide experimental data in the context of genome environment, for instance, measurements of gene expression, location of transcription factor binding sites/DNA sequence motifs, ChIP-on-chip/ChIP-seq mapping of histone tail and CpG methylation modifications. Gene expression array data can be presented in two complementary ways in the histogram display (Figure 3A). First, the number of differentially expressed genes (DEGs) (up- or downregulated) per Mb interval can be plotted alongside other genome features of interest. Second, taking into account some regions are inherently more gene-rich than others, GEB also allows plotting in each window the proportion of DEGs normalised to the underlying gene density. Specific regions of interest can be explored in detail, again in the physical map display, with an added functionality where genes are colour-coded to reflect their expression status (downregulated, upregulated, or no change/data unavailable). Additionally, the cut-offs/thresholds defining differential expression can be altered in real-time using an incorporated slide-bar, allowing the user to observe how the pattern of DEG distribution changes (Figure 3B).
GEB's flexible and efficient browsing capabilities are pivotal in the visualisation of any patterns embedded in densely-tiled probes over long genomic distances, as in the case of ChIP-on-chip data. The interface for displaying data from tiling arrays shares many of the features described above for gene expression data, with two specific features in the physical map display accommodating for the higher probe density. A "glyphs" track has been designed for appreciating global patterns, where all or a subset of probes can be colour-coded according to their ChIP-enrichment levels/status and displayed linearly (as if the probes were individual genes) (Figure 4A). To discern local patterns, the exact ChIP-enrichment levels for each probe can also be plotted in the "graphs" track with an adjustable scale for the Y-axis (Figure 4B).
Comparison of GEB with other genome browsers
Relative to major public genome browsers (e.g. Ensembl Genome Browser , UCSC Genome Browser  and NCBI Map Viewer ), GEB provides novel visualisation tools to analyse patterns of genome features and experimental data along an entire chromosome, and additionally the capability to move rapidly to a scalable detailed view with a dynamic range between 25 megabases and 2 basepairs resolution. Thus, users can move seamlessly through data covering a large genomic region, for example zooming in and out from megabase scale gene clusters to specific genes, without the need to change to a new page. Additionally GEB has relatively detailed annotation of interspersed repeats, including directionality, and provides simple tools to quantify the density/frequency of a given genomic feature within user defined limits. The latter are important for preliminary quantitative data analysis in cases where interesting patterns have been visualised.
GEB also has unique features relative to more recently developed track-based public browsers for integrating genome and experimental data, such as Argo Genome Browser  and Integrative Genomics Viewer , notably the detailed annotation of dispersed repeats and unique display options of GEB, the use of the anti-parallel "top and bottom" strand representation (as used in Ensembl) for detecting potential strand-bias in the distribution of selected features (e.g. gene clusters), the interactive visualisation tools for chromosome-wide feature density/distribution, and the quantitative data analysis tools.
GEB offers a number of unique and novel features and as such complements existing tools for use in studies relating whole genome scale experimental data with underlying genome features. The software is currently available at http://web.bioinformatics.ic.ac.uk/geb/. The program is easily adaptable for viewing all kinds of data in the context of genome environment from different sequenced genomes, provided the chromosomal start and end positions for each feature are known, e.g. data from ChIP-sequencing and chromosome conformation capture. The two-dimensional display used for L1 data in our example is potentially useful for other repetitive features which vary substantially in size. A complete set of GEB documentation as well as files for installing and running GEB is provided in additional file 3.
Availability and requirements
Project name: GEB
Project home page: http://web.bioinformatics.ic.ac.uk/geb/
Operating system(s): Unix/Linux for data retrieval and processing, platform independent for MySQL database and Java visualisation software
Programming language: Perl and Java
Other requirements: Java 1.6+, Perl, including Ensembl Perl API, and MySQL
License: GNU GPL
Any restrictions to use by non-academics: None
Genome Environment Browser
long interspersed nuclear element
short interspersed nuclear element
long terminal repeat
open reading frame
differentially expressed genes.
Lander ES, Linton LM, Birren B, Nusbaum C, Zody MC, Baldwin J, Devon K, Dewar K, Doyle M, FitzHugh W, et al.: Initial sequencing and analysis of the human genome. Nature 2001, 409(6822):860–921. 10.1038/35057062
Waterston RH, Lindblad-Toh K, Birney E, Rogers J, Abril JF, Agarwal P, Agarwala R, Ainscough R, Alexandersson M, An P, et al.: Initial sequencing and comparative analysis of the mouse genome. Nature 2002, 420(6915):520–562. 10.1038/nature01262
Gardiner-Garden M, Frommer M: CpG islands in vertebrate genomes. J Mol Biol 1987, 196(2):261–282. 10.1016/0022-2836(87)90689-9
Mirkovitch J, Mirault ME, Laemmli UK: Organization of the higher-order chromatin loop: specific DNA attachment sites on nuclear scaffold. Cell 1984, 39(1):223–232. 10.1016/0092-8674(84)90208-3
Kumar PP, Bischof O, Purbey PK, Notani D, Urlaub H, Dejean A, Galande S: Functional interaction between PML and SATB1 regulates chromatin-loop architecture and transcription of the MHC class I locus. Nat Cell Biol 2007, 9(1):45–56. 10.1038/ncb1516
Peaston AE, Evsikov AV, Graber JH, de Vries WN, Holbrook AE, Solter D, Knowles BB: Retrotransposons regulate host genes in mouse oocytes and preimplantation embryos. Dev Cell 2004, 7(4):597–606. 10.1016/j.devcel.2004.09.004
Han JS, Szak ST, Boeke JD: Transcriptional disruption by the L1 retrotransposon and implications for mammalian transcriptomes. Nature 2004, 429(6989):268–274. 10.1038/nature02536
Rice P, Longden I, Bleasby A: EMBOSS: the European Molecular Biology Open Software Suite. Trends Genet 2000, 16(6):276–277. 10.1016/S0168-9525(00)02024-2
Lyon MF: X-chromosome inactivation: a repeat hypothesis. Cytogenet Cell Genet 1998, 80(1–4):133–137. 10.1159/000014969
Ostertag EM, Kazazian HH Jr: Biology of mammalian L1 retrotransposons. Annu Rev Genet 2001, 35: 501–538. 10.1146/annurev.genet.35.102401.091032
Antequera F, Bird A: CpG islands as genomic footprints of promoters that are associated with replication origins. Curr Biol 1999, 9(17):R661–667. 10.1016/S0960-9822(99)80418-7
Boyle AL, Ballard SG, Ward DC: Differential distribution of long and short interspersed element sequences in the mouse genome: chromosome karyotyping by fluorescence in situ hybridization. Proc Natl Acad Sci USA 1990, 87(19):7757–7761. 10.1073/pnas.87.19.7757
Korenberg JR, Rykowski MC: Human genome organization: Alu, lines, and the molecular structure of metaphase chromosome bands. Cell 1988, 53(3):391–400. 10.1016/0092-8674(88)90159-6
Ensembl Genome Browser[http://www.ensembl.org]
UCSC Genome Browser[http://genome.ucsc.edu]
NCBI Map Viewer[http://www.ncbi.nlm.nih.gov/mapview/]
Argo Genome Browser[http://www.broad.mit.edu/annotation/argo]
Integrative Genomics Viewer[http://broad.mit.edu/igv/index.html]
This project was funded in part by the Medical Research Council, UK. YAT was supported by an MRC Pre-doctoral fellowship.
DH designed the software, wrote the source code and drafted the manuscript. YAT conceived of the study, designed the software and drafted the manuscript. NB and TBN conceived of the study, participated in its design and coordination, and edited the manuscript. SB coordinated the study. All authors read and approved the final manuscript.
Derek Huntley, Y Amy Tang contributed equally to this work.
Electronic supplementary material
Additional file 3: An archive of files for the GEB Java display, GEB set up and documentation. This includes all the files mentioned in the manuscript and GEB installation guide, e.g. the Java display file (GEB.jar), the sample configuration file (GEB.ini) and the Perl scripts/modules required to build the dedicated GEB database. Additional file 3 is a zipped archive file. We recommend using gunzip on the Mac (OSX) and Unix/Linux, and WinZip or WinRAR on Windows to open the file. (GZ 3 MB)
About this article
Cite this article
Huntley, D., Tang, Y.A., Nesterova, T.B. et al. Genome Environment Browser (GEB): a dynamic browser for visualising high-throughput experimental data in the context of genome features. BMC Bioinformatics 9, 501 (2008). https://doi.org/10.1186/1471-2105-9-501
- Long Terminal Repeat
- Genome Feature
- Long Terminal Repeat Retrotransposons
- Genome Environment