Skip to main content

The SeqWord Genome Browser: an online tool for the identification and visualization of atypical regions of bacterial genomes through oligonucleotide usage

Abstract

Background

Data mining in large DNA sequences is a major challenge in microbial genomics and bioinformatics. Oligonucleotide usage (OU) patterns provide a wealth of information for large scale sequence analysis and visualization. The purpose of this research was to make OU statistical analysis available as a novel web-based tool for functional genomics and annotation. The tool is also available as a downloadable package.

Results

The SeqWord Genome Browser (SWGB) was developed to visualize the natural compositional variation of DNA sequences. The applet is also used for identification of divergent genomic regions both in annotated sequences of bacterial chromosomes, plasmids, phages and viruses, and in raw DNA sequences prior to annotation by comparing local and global OU patterns. The applet allows fast and reliable identification of clusters of horizontally transferred genomic islands, large multi-domain genes and genes for ribosomal RNA. Within the majority of genomic fragments (also termed genomic core sequence), regions enriched with housekeeping genes, ribosomal proteins and the regions rich in pseudogenes or genetic vestiges may be contrasted.

Conclusion

The SWGB applet presents a range of comprehensive OU statistical parameters calculated for a range of bacterial species, plasmids and phages. It is available on the Internet at http://www.bi.up.ac.za/SeqWord/mhhapplet.php.

Background

The study of genome OU signatures has a long history dating back to early publications by Karlin et al. who focused mainly on dinucleotide compositional biases and their evolutionary implications [13]. Statistical approaches of OU comparison were further advanced by Deschavanne et al., who applied chaos game algorithms [4]; and by Pride et al., who extended the analysis to tetranucleotides using Markov Chain Model simulations [5]. Later, a number of practical tools for phylogenetic comparison of bacterial genomes [4, 6, 7], identification of horizontally transferred genomic islands [813] and assignment of unknown genomic sequences [14, 15] based on OU statistics became publicly available. These approaches exploited the notion that genomic OU composition was less variable within genomes rather than between them, regardless of which genomic regions had been taken into consideration [16]. A general belief was that if a significant compositional difference was discovered in genomic fragments relative to the core genome, these loci most likely can be assigned to horizontally transferred genetic elements (transposons, prophages or integrated plasmids). This approach was criticized by several researchers [17, 18] who pointed out that codon bias and base composition are poor indicators of horizontal gene transfer. Therefore, there is a need for more informative parameters which also take into account higher order DNA variation. An overview of the current OU statistical methods based on di-, tetra- and hexanucleotides has been published recently. The conclusion of the review was that all methods were context dependent and, though being efficient and powerful, none of them were superior in all applications [19]. Thus, the major motivation of our work was to develop more flexible and informative algorithms seamlessly integrating di- to heptanucleotides OU analysis for reliable identification of divergent genomic regions.

Recently we have introduced the concept of OU patterns into the literature [20]. Each OU pattern is characterized by a number of OU statistical parameters namely, local pattern deviation (D), pattern skew (PS), relative variance (RV) and others (see Methods section). Novelties of the developed algorithms relative to other existing methods include the following: i) distances between patterns of different word length (from di- through to heptanucleotides) calculated for the same sequences are comparable; i.e. one may use longer word patterns to perform a large scale analysis and then switch to shorter word patterns for a more detailed view; ii) OU patterns calculated for sequences of different lengths are comparable provided that the length of the sequence is longer than the corresponding thresholds (specified in the Methods section); iii) alterations of OU patterns may be analyzed by different non-redundant parameters (D, PS and RV with different schemes of normalization by frequencies of shorter constituent words). Superimposition of these OU characteristics allows better discrimination of divergent genomic regions relative to other contemporary approaches [21].

Implementation

Calculation of OU statistical parameters has been described previously [20, 21]. OU pattern was denoted as a matrix of deviations Δ[ξ 1...ξ N]of observed from expected counts for all possible words of length N:

Δ[ξ 1...ξ N]= (C[ξ 1...ξ N]|obs- C[ξ 1...ξ N]|e)/C[ξ 1...ξ N]|0 (1)

where ξn is any nucleotide A, T, G or C in the N-long word; C[ξ 1...ξ N]|obsis the observed count of the word [ξ1...ξN]; C[ξ 1...ξ N]|eis the expected count and C[ξ 1...ξ N]|0is a standard count estimated from the assumption of an equal distribution of words in the sequence: (C[ξ 1...ξ N]|0= L seq × 4-N).

Expected counts of words C[ξ 1...ξ N]|ewere calculated in accordance with the applied normalization scheme. Thus, C[ξ 1...ξ N]|e= C[ξ 1...ξ N]|0if OU is not normalized, or C[ξ 1...ξ N]|e= C[ξ 1...ξ N]|nif OU is normalized by empirical frequencies of all shorter words of the length n. The expected count of a word C[ξ 1...ξ N]|eof length N in a L seq long sequence normalized by frequencies of n-mers (n <N) was calculated as follows:

C [ ξ 1 ... ξ l w ] | n = L s e q × F [ ξ 1 ... ξ n ] × i = 2 N n + 1 ( F [ ξ i ... ξ i + n 1 ] ξ i + n ξ A , T , G , C F [ ξ i ... ξ i + n ] ξ ) MathType@MTEF@5@5@+=feaagaart1ev2aaatCvAUfKttLearuWrP9MDH5MBPbIqV92AaeXatLxBI9gBaebbnrfifHhDYfgasaacPC6xNi=xI8qiVKYPFjYdHaVhbbf9v8qqaqFr0xc9vqFj0dXdbba91qpepeI8k8fiI+fsY=rqGqVepae9pg0db9vqaiVgFr0xfr=xfr=xc9adbaqaaeGaciGaaiaabeqaaeqabiWaaaGcbaGaem4qam0aaSbaaSqaaiabcUfaBjabe67a4naaBaaameaacqaIXaqmaeqaaSGaeiOla4IaeiOla4IaeiOla4IaeqOVdG3aaSbaaWqaaiabdYgaSnaaBaaabaGaem4DaChabeaaaeqaaSGaeiyxa0LaeiiFaWNaemOBa4gabeaakiabg2da9iabdYeamnaaBaaaleaacqWGZbWCcqWGLbqzcqWGXbqCaeqaaOGaey41aqRaemOray0aaSbaaSqaaiabcUfaBjabe67a4naaBaaameaacqaIXaqmaeqaaSGaeiOla4IaeiOla4IaeiOla4IaeqOVdG3aaSbaaWqaaiabd6gaUbqabaWccqGGDbqxaeqaaOGaey41aq7aaebCaeaadaqadaqcfayaamaalaaabaGaemOray0aaSbaaeaacqGGBbWwcqaH+oaEdaWgaaqaaiabdMgaPbqabaGaeiOla4IaeiOla4IaeiOla4IaeqOVdG3aaSbaaeaacqWGPbqAcqGHRaWkcqWGUbGBcqGHsislcqaIXaqmaeqaaiabc2faDjabe67a4naaBaaabaGaemyAaKMaey4kaSIaemOBa4gabeaaaeqaaaqaamaaqahabaGaemOray0aaSbaaeaacqGGBbWwcqaH+oaEdaWgaaqaaiabdMgaPbqabaGaeiOla4IaeiOla4IaeiOla4IaeqOVdG3aaSbaaeaacqWGPbqAcqGHRaWkcqWGUbGBaeqaaiabc2faDjabe67a4bqabaaabaGaeqOVdGhabaGaemyqaeKaeiilaWIaemivaqLaeiilaWIaem4raCKaeiilaWIaem4qameacqGHris5aaaaaOGaayjkaiaawMcaaaWcbaGaemyAaKMaeyypa0JaeGOmaidabaGaemOta4KaeyOeI0IaemOBa4Maey4kaSIaeGymaedaniabg+Givdaaaa@9433@
(2)

where the F[ξ 1...ξ n]values are the observed frequencies of the particular word of length n in the sequence and ξ is any nucleotide A, T, G or C. For example, expected count of a word ATGC in a sequence of L seq nucleotides normalized by frequencies of trinucleotides is:

C A T G C = L s e q × F A T G × F T G C F T G A + F T G T + F T G G + F T G C MathType@MTEF@5@5@+=feaagaart1ev2aaatCvAUfKttLearuWrP9MDH5MBPbIqV92AaeXatLxBI9gBaebbnrfifHhDYfgasaacPC6xNi=xI8qiVKYPFjYdHaVhbbf9v8qqaqFr0xc9vqFj0dXdbba91qpepeI8k8fiI+fsY=rqGqVepae9pg0db9vqaiVgFr0xfr=xfr=xc9adbaqaaeGaciGaaiaabeqaaeqabiWaaaGcbaGaem4qam0aaSbaaSqaaiabdgeabjabdsfaujabdEeahjabdoeadbqabaGccqGH9aqpcqWGmbatdaWgaaWcbaGaem4CamNaemyzauMaemyCaehabeaakiabgEna0kabdAeagnaaBaaaleaacqWGbbqqcqWGubavcqWGhbWraeqaaOGaey41aqBcfa4aaSaaaeaacqWGgbGrdaWgaaqaaiabdsfaujabdEeahjabdoeadbqabaaabaGaemOray0aaSbaaeaacqWGubavcqWGhbWrcqWGbbqqaeqaaiabgUcaRiabdAeagnaaBaaabaGaemivaqLaem4raCKaemivaqfabeaacqGHRaWkcqWGgbGrdaWgaaqaaiabdsfaujabdEeahjabdEeahbqabaGaey4kaSIaemOray0aaSbaaeaacqWGubavcqWGhbWrcqWGdbWqaeqaaaaaaaa@5B4B@
(3)

Two approaches of normalization have been exploited where the F values were calculated for the complete sequence of a chromosome, plasmid, etc (generalized normalization) or for a given sliding window (local normalization). The normalization by equation 2 allows identification of words, frequencies of which cannot be predicted exactly by frequencies of shorter constituent words.

The distance D between two patterns was calculated as the sum of absolute distances between ranks of identical words (w, in a total 4Ndifferent words) after ordering of words by Δ[ξ 1...ξ N]values (see equation 1) in patterns i and j as follows:

D ( % ) = 100 × w 4 N | r a n k w , i r a n k w , j | D min D max D min MathType@MTEF@5@5@+=feaagaart1ev2aaatCvAUfKttLearuWrP9MDH5MBPbIqV92AaeXatLxBI9gBaebbnrfifHhDYfgasaacPC6xNi=xI8qiVKYPFjYdHaVhbbf9v8qqaqFr0xc9vqFj0dXdbba91qpepeI8k8fiI+fsY=rqGqVepae9pg0db9vqaiVgFr0xfr=xfr=xc9adbaqaaeGaciGaaiaabeqaaeqabiWaaaGcbaGaemiraqKaeiikaGIaeiyjauIaeiykaKIaeyypa0JaeGymaeJaeGimaaJaeGimaaJaey41aqBcfa4aaSaaaeaadaaeWbqaamaaemaabaGaemOCaiNaemyyaeMaemOBa4Maem4AaS2aaSbaaeaacqWG3bWDcqGGSaalcqWGPbqAaeqaaiabgkHiTiabdkhaYjabdggaHjabd6gaUjabdUgaRnaaBaaabaGaem4DaCNaeiilaWIaemOAaOgabeaaaiaawEa7caGLiWoacqGHsislcqWGebardaWgaaqaaiGbc2gaTjabcMgaPjabc6gaUbqabaaabaGaem4DaChabaGaeGinaqZaaWbaaeqabaGaemOta4eaaaGaeyyeIuoaaeaacqWGebardaWgaaqaaiGbc2gaTjabcggaHjabcIha4bqabaGaeyOeI0Iaemiraq0aaSbaaeaacyGGTbqBcqGGPbqAcqGGUbGBaeqaaaaaaaa@64A6@
(4)

Application of ranks instead of relative oligonucleotide frequency statistics made the comparison of OU patterns less biased to the sequence length provided that the sequences are longer than the limits of 0.3, 1.2, 5, 18.5, 74 and 295 kbp for di-, tri-, tetra-, penta-, hexa- and heptanucleotides, respectively [20].

PS is a particular case of D where patterns i and j were calculated for the same DNA but for direct and reversed strands, respectively. Dmax = 4N× (4N- 1)/2 and Dmin = 0 when calculating a D or, in a case of PS calculation, Dmin = 4Nif N is an odd number or Dmin = 4N- 2Nif N is an even number due to presence of palindromic words [20]. Normalization of D-values by Dmax ensures that the distances between two sequences are comparable regardless of the word length of OU patterns.

Relative variance of an OU pattern was calculated by the following equation:

R V = w 4 N Δ w 2 ( 4 N 1 ) σ 0 2 MathType@MTEF@5@5@+=feaagaart1ev2aaatCvAUfKttLearuWrP9MDH5MBPbIqV92AaeXatLxBI9gBaebbnrfifHhDYfgasaacPC6xNi=xI8qiVKYPFjYdHaVhbbf9v8qqaqFr0xc9vqFj0dXdbba91qpepeI8k8fiI+fsY=rqGqVepae9pg0db9vqaiVgFr0xfr=xfr=xc9adbaqaaeGaciGaaiaabeqaaeqabiWaaaGcbaGaemOuaiLaemOvayLaeyypa0tcfa4aaSaaaeaadaaeWbqaaiabfs5aenaaDaaabaGaem4DaChabaGaeGOmaidaaaqaaiabdEha3bqaaiabisda0maaCaaabeqaaiabd6eaobaaaiabggHiLdaabaWaaeWaaeaacqaI0aandaahaaqabeaacqWGobGtaaGaeyOeI0IaeGymaedacaGLOaGaayzkaaGaeq4Wdm3aa0baaeaacqaIWaamaeaacqaIYaGmaaaaaaaa@4360@
(5)

where N is word length; Δ2 w is the square of a word w count deviation (see equation 1); and σ20 is the expected variance of the word distribution in a randomly generated sequence that depends on the sequence length and the word length:

σ 0 2 = 0.14 + 4 N L s e q MathType@MTEF@5@5@+=feaagaart1ev2aaatCvAUfKttLearuWrP9MDH5MBPbIqV92AaeXatLxBI9gBaebbnrfifHhDYfgasaacPC6xNi=xI8qiVKYPFjYdHaVhbbf9v8qqaqFr0xc9vqFj0dXdbba91qpepeI8k8fiI+fsY=rqGqVepae9pg0db9vqaiVgFr0xfr=xfr=xc9adbaqaaeGaciGaaiaabeqaaeqabiWaaaGcbaGaeq4Wdm3aa0baaSqaaiabicdaWaqaaiabikdaYaaakiabg2da9iabicdaWiabc6caUiabigdaXiabisda0iabgUcaRKqbaoaalaaabaGaeGinaqZaaWbaaeqabaGaemOta4eaaaqaaiabdYeamnaaBaaabaGaem4CamNaemyzauMaemyCaehabeaaaaaaaa@3DE7@
(6)

where L seq is sequence length, and N is word length. Normalization of OU pattern variance by σ0 makes the variances comparable regardless of the word length of OU patterns and the sequence length. The regression equation was tested on 300 randomly generated sequences with an equiprobable occurrence of all 4 nucleotides by the DataFit 7.1.44 software.

The SWGB is coded in Java to be used as an applet in a Web-browser either on the Internet or locally (the programs OligoWords in Python and SeqWord_Viewer, which respectively calculate and visualize the OU patterns for DNA sequences, are available for download from the SWGB website). SWGB should run on any platform with a Java 1.5.x runtime environment or newer.

The pre-calculated data-sets are saved in a MySQL Server 5.0 database. The size of the sliding window and the OU pattern type were applied according to the sequence length (Table 1). At the time of writing, the SeqWord database contained OU patterns pre-calculated for the sequences of 682 bacterial chromosomes belonging to 637 different organisms (strains and species), 412 plasmids, 100 bacteriophages and 39 other viruses, which were downloaded from the NCBI [44].

Table 1 Sliding window size and OU pattern types (oligomer lengths) selected for sequences of different length present in the SeqWord database.

Results

User familiarity with the abbreviations of the various OU statistical parameters is important. Different types of OU patterns were abbreviated as type_N mer. Types might be "n0" for non-normalized, or "n1" for normalized by mononucleotide frequencies. For example, the non-normalized tetranucleotide usage pattern is denoted as n0_4mer; tetranucleotide usage pattern normalized by mononucleotide content is n1_4mer etc. The genomes in the SWGB database were analyzed by the following statistical parameters: D – distance between two patterns of the same type (in this work we used distances D between local patterns calculated for overlapping genome fragments and the global genome patterns calculated for the complete sequence – the local pattern deviation); PS – pattern skew, distance between the two patterns of the direct and reverse strands of the same DNA sequence; RV and GRV – oligonucleotide usage variances normalized locally and globally, respectively, and reduced to the OU variance expected for a randomly generated sequence (see Materials section); GC-content (GC) and GC-skew (GCS) in DNA fragments.

The applet GUI and database of pre-calculated OU patterns

The SeqWord Genome Browser (SWGB) applet is available via the Internet [2224] and is mouse and menu driven. The Web-based applet is used to visualize DNA compositional variations in bacterial and viral genomes stored in the SeqWord database. Every genome in the database is represented by a set of statistical OU parameters (D, PS, GV, GRV, GC and GCS) calculated for genomic fragments, which were selected by a sliding window (sliding window length and step were set according to the total length of the sequence as demonstrated in Table 1). While in 70 to 99% of genomic fragments the OU compositional bias is similar to the complete genome OU pattern, some regions with atypical OU composition, however, are always present. Superimposition of different OU parameters allows discrimination of divergent genomic regions, as was published previously [21]. Briefly: rRNA operons are characterized by extremely high PS and low RV; giant genes with multiple repeated elements have high or moderate PS and high RV; horizontally transferred genetic elements are characterized by increased divergence between RV and GRV accompanied by high D; and genes for ribosomal proteins show a moderate increase of D, PS and RV above genomic averages. Having analyzed 1243 sequences of different microorganisms including viruses and plasmids in the SeqWord database, we confirmed that the approaches we have developed and tested previously [25] (mainly on Pseudomonas putida KT2440 chromosomal DNA) are appropriate and useful for analysis of genomic sequences of other microorganisms and viruses.

In an open applet window, the user has the ability to choose from an ever growing list of available sequences (Fig. 1). The user also has the option of restricting the list to display only bacterial chromosomes, plasmids, phages, viruses or all sequences by selecting the corresponding filter button. Users have to select a genome in the list and click the 'Display in the Applet' button to retrieve the pre-calculated data. All OU parameters calculated for a given genome may be exported to a local text file by using the 'Export' function from the applet's 'File' menu. Later, instead of again having to connect to the database, users may open and view their local files (previously exported from the applet or calculated by the OligoWords program, see below) via the 'Open' function in the 'File' menu.

Figure 1
figure 1

General view of the web-based SWGB with a list of genomes present in the database and an enclosed Java applet for data visualization. To show OU statistical parameters for a selected genome, click the 'Display in the Applet' button. Click a filter button to order genomes by the corresponding category and use the interactive letters at the top to scroll the list to a sequence of interest.

The SWGB is basically comprised of two views, denoted by the 'Gene Map' and 'Diagram' tabs. The applet is instrumental for visualization of natural variation in DNA sequences by the interactive diagrams on the 'Gene Map' and 'Diagram' tabs. Users may save the current diagram in JPG format by using the 'Save picture' function in the 'File' menu.

The 'Gene Map' tab offers a simple view of an entire genome at a glance and gives users access to a number of important pre-calculated OU statistics superimposed on the gene map (Fig. 2). Displays for each of the statistical parameters can be toggled on/off by checking items in the 'Hide Elements' menu. By merely mousing over any region on the plot, a message displaying detailed information for the pointed curve will be shown in the 'Message' bar. Clicking a gene on the map displays a dialog with the annotation details (Fig. 2).

Figure 2
figure 2

Identification of divergent genomic regions on the 'Gene Map' view. Superimposition of different OU parameters such as GC (black line), GCS (pink), PS (green), D (blue), GRV (upper brown line) and RV (lower brown line) allows discrimination of divergent genomic regions. In this example a part of the chromosome of Pseudomonas putida KT2440 (127–774 kbp) is displayed in the applet window. A genomic fragment was highlighted using the function 'Select region' and a giant gene, PP0168, was selected by 'Select gene'. A pop-up window 'Gene Details' was opened by double-clicking the gene on the map. Genes are indicated by red and grey (for hypotheticals) bars. The black horizontal line separates genes by their direction of translation.

The 'Zoom' function is straight-forward and allows users to control the amount of data viewed in the plot area. Clicking the 'Enter' button after setting the desired zoom value will then redraw the map. A 'Zoom into region' function under the 'Tools' drop-down menu allows users to zoom into exact genomic regions by merely entering their desired co-ordinates into the pop-up dialog box. The 'Tools' → 'Select region' menu item allows highlighting of selected regions without zooming. Use the option 'Clear ...' in the 'Tools' menu to undo zooming or highlighting. To locate a genomic region by gene, click the button 'Select Gene'. In the pop-up dialog box one may order the gene list by gene names, functionality or coordinates, then select a gene in the list and click 'OK'. When a gene annotation is not available, the values of the locus coordinates are used as a gene name. The applet window will be scrolled to the selected gene highlighted on the map (see Fig. 2).

The 'Diagram' tab allows flexible filtering of the underlying data based on the criteria chosen by users. Although the underlying data is pre-calculated, the user may, by simply changing selected parameters, generate very different images which give different insights into the natural genomic variation. To start with, the 'Diagram' view offers a bar chart or a dot-plot presentation of the pre-calculated data. To view a bar chart of the distribution statistics for a given OU parameter, select the desired parameters from the X or Y-axis drop-downs and click 'Enter'. The number of bars displayed can be adjusted using the '# Bars' selector.

On the dot-plot diagram, each genomic fragment (selected by the sliding window) is represented by a dot with X and Y coordinates that correspond to values of OU parameters chosen from X and Y drop-down lists, respectively. The Z axis parameter may be set as well. In this case, the dots are coloured by values of OU parameters selected for the Z axis, and the colour range is displayed on the vertical colour bar on the left of the plot area (Fig. 3).

Figure 3
figure 3

The 'Diagram' view. To draw a diagram, first select corresponding OU parameters using the dropdown lists and click the 'Enter' button. In this example n1_4mer:RV, n1_4mer:GRV and n0_4mer:D were selected for the X, Y and Z axes, respectively. Every dot on the dot-plot corresponds to a genomic fragment selected by the sliding window. Dots are spread and coloured in accordance with their values of the selected statistical OU parameters. Information for each dot may be found by one of the following methods: i) information for a dot under the mouse pointed by the mouse is shown in the 'Message' bar; ii) double clicking a dot returns us to the 'Gene map' tab with the corresponding genomic fragment highlighted; iii) framing the dots and clicking the 'Get' button opens a new applet window with the information about all selected regions. In this example the genomic regions of Salmonella typhimurium LT2 (NC_003197) that correspond to horizontally transferred genetic elements were selected (see discussion in the text).

Having set up the dot-plot, users will be able to identify divergent genomic regions (see next section). To retrieve annotations of genomic fragments corresponding to a group of dots, frame the dots of interest by clicking and dragging over the desired area. A selector frame then appears around the dots (Fig. 3). Clicking the 'Get' button displays the selected genomic fragments with their coordinates and gene annotations. Furthermore, identification and isolation of specific genomic regions may be improved significantly by filtering dots by OU parameters. The simplest way of filtering is by the third (Z axis) parameter. One may select an area on the colour bar to exclude all dots from the plot lying outside of the selected colour range (see an example in help files on-line). The hidden dots will not be selected by the 'Get' button. A more sophisticated way to filter genomic regions is provided by the 'Filter' button. An example will be discussed below.

The 'Mark' button enables genomic fragments to be selected by their coordinates and highlighted on the dot-plot. Click the 'Mark' button to open a dialog and enter coordinates of one or multiple fragments (Fig. 4). Co-ordinates of each fragment must be added to the list by clicking the 'Add' button. Close the dialog by clicking 'OK'. The corresponding dots on the dot-plot will be highlighted as shown in Fig. 4.

Figure 4
figure 4

Identification of divergent genomic regions by plotting and highlighting. In this example the genome of Synechococcus sp. WH8102 was analysed. The parameters n0_4mer:D, n1_4mer:RV and n0_4mer:PS were selected for the X, Y and Z axes, respectively. The genomic regions covering the giant gene for the surface protein SwmB [29] were highlighted by entering the coordinates of this gene into the 'Mark loci by coordinates' dialog. The genomic regions enriched with i) housekeeping genes; ii) genes for ribosomal proteins; iii) vestigial genetic elements (comprising pseudogenes, transposons, prophages and IS-elements) are indicated.

Identification of divergent genomic islands

Several routines have been developed to identify the horizontally transferred genomic islands, genes for ribosomal RNA and proteins, non-functional pseudogenes and genes of other functional categories. All these routines are described in detail with illustrations in supplementary web-pages (use the 'Help' link in the applet window).

The approach to identify inserts of foreign genomic elements by OU statistical parameters has been described recently [21]. While several algorithms allow identification of horizontally transferred genomic islands [813], the multiple oligomer parameters used in the SWGB even allows tentative attribution of genomic fragments (and, given the right scale, genes or gene clusters) to different functional classes using only a FASTA sequence as input. However, the emphasis of the SWGB is not primarily its annotation capability, but its ability to display the natural internal variability of genome sequences. We use Pseudomonas putida KT2440, a known mosaic genome with 105 genomic islands above 4000 bp in length [26] as an example. Many of these features can be visualized at a glance using the SWGB without any in depth analysis (see Fig. 2). On the 'Diagram' view the parameters n1_4mer:RV, n1_4mer:GRV and n0_4mer:D were selected for the X, Y and Z axes, respectively, as we showed previously (see Fig. 3). Plotting local relative oligomer variance (RV) against global relative variance (GRV) basically shows the effect of normalization by global mononucleotide content. The core genome is then represented on the dot plot as the positive linear correlation line where RV ≈ GRV (Fig. 3). In other words, these fragments exhibit such compositional closeness to the core genome that normalizing by local mononucleotide content does not have a different effect compared to normalizing by global content. These genomic fragments also exhibit a low distance from the genomic average; and are therefore coloured blue. Scattered dots lying peripheral to the expected strong linear correlation do not belong to the core genome and also have a higher distance from the genomic average and are hence coloured green. Using the filter settings recommended in Fig. 5, twenty one fragments were found to be genomic islands (note that while border values of OU parameters are not the same for different genomes, the grading notches of the sliders represent relative values that allows identification of homologous regions in many different genomes). For a number of reasons, many more islands were found in a similar analysis by Weinel et al. [26]. Firstly, the sliding window size of 8 kbp means many of the 4 kbp features from their analysis were not identified automatically. Furthermore, they were looking for all compositionally atypical regions, whereas here we restrict ourselves to horizontally transferred regions.

Figure 5
figure 5

Filtering genomic regions by multiple parameters. Click the 'Filter' button to open a dialog as shown in the figure. Setting up border values of multiple OU statistical parameters allows more precise localisation of regions of interest.

A known 40 kbp bacteriophage insertion [2586000–2626000] is, surprisingly, not among the genomic fragments selected in the SWGB using this filter. Although the prophage is still perceptible on the 'Gene Map' view (see a figure in the supplementary help web-pages), the OU parameters of the region do not differ markedly enough from the core sequence to be isolated automatically as a horizontally transferred region.

As the SWGB uses parameters that are based on comparison of local fragments to the global genomic average, strains with abundant insertions of homogenous DNA can confound this form of analysis. One example is the Methanosarcina acetivorans C2A genome which is composed of an estimated 25% of putatively horizontally acquired DNA, one of the highest amounts discovered to date [11]. As a result of these insertions, the genomic signature has been strongly influenced, resulting in a large amount of scatter and a poorly defined core genome on the plots. On the other hand, this type of analysis allows estimation of genome stability in a simple, multi parameter view (see the Vibrio cholerae N16961-O1-eltor example in the online help files). To conclude, filtering provides a convenient way to automatically isolate divergent genomic regions of interest. However, some regions may erroneously remain undetected due to possible amelioration of older inserts [27] or a higher level of noise in unstable genomes. However, many problematic genomic fragments can in some cases be easily attributed to functional gene categories using the SWGB 'Diagram' window (see Fig. 2).

Methodologies for discovering long modular genes have already been discussed in a previous publication [28]. Briefly, long genes display a particular tetranucleotide usage and can be discovered by plotting n0_4mer:D (X axis) versus n1_4mer:RV (Y axis). The positively linear correlated outlier fragments (towards the top right of the image) are often fragments of long genes with their characteristic repeats. An example using the gene encoding the 1.12 megadalton cell surface protein of Synechococcus sp. WH8102 [29] marked on the dot-plot is shown in Fig. 4. Ribosomal RNA operons (but not genes for ribosomal proteins) are characterized by extremely high pattern skew and a large distance from the core genome (Fig. 2). Thus, there is a tendency to find many genomic fragments containing rRNA genes coloured dark brown to red in the bottom right section of the 'Diagram' tab. The annotation for rRNA operons is not present in the database; therefore, these are seen in the 'Gene Map' tab as unannotated areas with high pattern skew (Fig. 2). Ribosomal proteins tend to be increasingly present at a slightly greater than average RV and above average D (see Fig. 2), which is in agreement with observations that highly expressed genes for ribosomal proteins have a highly specific codon usage compared to housekeeping genes of the organism [30]. The majority of genomic fragments form a cluster characterized by average and higher than average RV, stable OU patterns (low D) and low PS. These tend to be the core, or bulk genes and genomic regions with their typical tetranucleotide usage. Some other core sequence fragments spread from this area toward lower RV and less specific OU patterns (higher D and PS) – these are all characteristics of an unstable or randomly generated sequence [20]. These regions were found to be enriched with many hypothetical genes, prophages and transposons. (The data is not shown but is easily verified with any genome using the 'Get' button. Consider, for example, this area in the pseudogene rich Mycobacterium leprae TN or Methanosarcina acetivorans C2A genomes [11, 31], and the relatively homogenous Alcanivorax borkumensis SK2 genome [32].) These regions were thus categorized as rich in 'vestigial' genes in contrast to the core genome regions rich in housekeeping genes (Fig. 4).

It must be stressed that with an average length of genes being around 1 kbp and overlapping sliding windows of 8 kbp, one cannot expect precise separation of housekeeping and vestigial genes by the method described above. However, when analyzing an unknown DNA sequence prior to annotation, it may be helpful to identify genomic regions enriched with a higher proportion of these so called housekeeping genes and other regions rich in vestigial genes. These tentative results should be verified with other complementary algorithms such as BLAST, gene finding and annotation techniques.

The most important feature of the supplemented software available from the SWGB web-server for download is the ability to quickly and easily analyze a novel sequence on a local computer. The command-line Python program OligoWords is first used to analyse FASTA or GenBank formatted sequences. The program is available for download [33] in several packages as precompiled executable files and as Python source code. The command-line interface of the OligoWords program is shown in Fig. 6. Parameters such as oligomer length and window size can all be set depending on the sequence length and desired resolution (see Table 1 for suggestions). Since the SWGB is implemented as a Java applet, it can be run within a web browser locally. The HTML-embedded applet is available for download from the same FTP site [33] (select SeqWord_Viewer.zip). The output file from OligoWords is read into the SWGB via the 'Open' function of the 'File' menu, and the complete functionality of the online system is then available. For example, a new sequence can be analysed for ribosomal gene clusters, putative horizontally transferred elements or other regions of atypical DNA structure prior to the lengthy annotation step. A complete description of how to run the SWGB and OligoWords locally is presented in the online help files.

Figure 6
figure 6

Command-line interface of the OligoWords program. To change the setting for the current run, type the option's letter and enter a new value as prompted. Users may change: T) the set of statistical OU parameters to be calculated for every local pattern; L) length of the sliding window; S) step of the sliding window; I) the name of the input folder that contains FASTA and/or GenBank files with source DNA sequences; and O) the name of the output folder where the result files will be stored.

Conclusion

The SWGB applies novel OU statistics to visualize and discern divergent genomic regions. It has been extensively tested in practice for large scale genome analysis [32, 34], and for identification and comparison of horizontally transferred genomic islands [35]. The applet is linked to a database of pre-calculated OU patterns of bacterial genomes (1243 complete sequences, including bacterial chromosomes, plasmids and some viruses were available at the time of manuscript submission, however, new sequences are regularly being added). The SWGB allows tentative annotation of the various divergent regions and provides overviews for use in comparative genomics. Users may download the command line version of the OligoWords program to analyze their own sequences. A packaged version of the SWGB allows users to view and manipulate their OligoWords results locally using a compatible web-browser.

Although there are several readily available tools for DNA compositional analysis, genomic island identification and large scale genome analysis [3642], the SWGB surpasses previous approaches in making use of a wider range of parameters which allow identification of divergent genomic regions and even visual tentative attribution of these DNA fragments to various categories. We have found superimposition of these parameters to be more informative than a simple GC average or a relative OU frequency deviation since they allow discrimination of divergent genomic regions (large modular genes, ribosomal RNAs, ribosomal protein clusters and the horizontally transferred genomic islands, see Fig. 2) all of which are characterized by an alternative OU composition relative to the core sequence. In addition, our approach provides some insight into the physicochemical state of the analysed DNA and the stability/state of flux of a genome as tetranucleotides exert a strong structural signal [20, 21]. Consideration of flux inferred oligonucleotide usage is particularly interesting when comparing, for example, multiple replicating units of the same strain. Using the simple analysis described here, the second chromosome of Vibrio cholerae N16961-O1-eltor was demonstrated to be far less conserved than the first, with differences in mononucleotide content and distance from core genomic values implying a more heterogenous chromosome consistent with its role as a gene capture system [43].

Furthermore, no single oligonucleotide word size has been found to be optimal for all purposes, such as finding conserved or horizontally transferred DNA, plasmid host comparisons or testing distant homology [19]. The SWGB crucially provides the opportunity to analyse DNA sequences with various oligomer lengths and normalisation schemes. For example, genomic regions of particular interest may be multiply analysed with progressively smaller oligomer sizes to provide more detailed information on oligomer usage in individual genes.

Availability and requirements

The SWGB applet is freely available to any researcher wishing to use it for non-commercial purposes via the Internet [2224]. It has been tested on openSUSE 10.2, Gentoo Linux 2.6, Fedora Core 5 and Microsoft Windows XP workstations using Microsoft Internet Explorer 6.0, Maxthon 1.5.9, Mozilla Firefox 2.0, Mozilla SeaMonkey 1.1.1, Safari 3.0.4 for Mac, Konqueror 3.5.5 and 3.5.7, and Opera 9.10 browsers with Java 1.5. At the time of manuscript submission, a problem likely related to the local firewall was encountered with the Firefox browser on SUSE 10.2 and some other browsers (see 'Compatibility' link on the SWGB front page). The problem will be tackled in later releases of the SWGB. Feedback from users (addressed to the corresponding author) is very much appreciated.

Abbreviations

SWGB:

SeqWord Genome Browser

OU:

oligonucleotide usage

D:

distance between two oligonucleotide usage patterns

PS:

oligonucleotide usage pattern skew

RV:

relative variance of the oligonucleotide usage

GRV:

globally normalized relative variance of the oligonucleotide usage

GC:

guanine + cytosine content

GCS:

guanine versus cytosine skew in DNA strands.

References

  1. Karlin S, Burge C: Dinucleotide relative abundance extremes: a genomic signature. Trends Genet 1995, 11(7):283–290. 10.1016/S0168-9525(00)89076-9

    Article  CAS  PubMed  Google Scholar 

  2. Karlin S, Mrázek J, Campbell A: Compositional biases of bacterial genomes and evolutionary implications. J Bacteriol 1997, 179(12):3899–3913.

    PubMed Central  CAS  PubMed  Google Scholar 

  3. Karlin S: Global dinucleotide signatures and analysis of genomic heterogeneity. Curr Opin Microbiol 1998, 1: 598–610. 10.1016/S1369-5274(98)80095-7

    Article  CAS  PubMed  Google Scholar 

  4. Deschavanne PJ, Giron A, Vilain J, Fagot G, Fertil B: Genomic signature: characterization and classification of species assessed by chaos game representation of sequences. Mol Biol Evol 1999, 16(10):1391–1399.

    Article  CAS  PubMed  Google Scholar 

  5. Pride DT, Meinersmann RJ, Wassenaar TM, Blaser MJ: Evolutionary implications of microbial genome tetanucleotide frequency biases. Genome Res 2003, 13: 145–155. 10.1101/gr.335003

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  6. Coenye T, Vandamme P: Use of the genomic signatures in bacterial classification and identification. System Appl Microbiol 2004, 27: 175–185. 10.1078/072320204322881790

    Article  CAS  Google Scholar 

  7. van Passel MW, Bart A, Luyf AC, van Kampen AH, Ende A: The reach of the genome signature in prokaryotes. BMC Evol Biol 2006, 6: 84. 10.1186/1471-2148-6-84

    Article  PubMed Central  PubMed  Google Scholar 

  8. Mrázek J, Karlin S: Detecting alien genes in bacterial genomes. Ann NY Acad Sci 1999, 870: 314–329. 10.1111/j.1749-6632.1999.tb08893.x

    Article  PubMed  Google Scholar 

  9. Azad RK, Lawrence JG: Use of artificial genomes in assessing methods for atypical gene detection. PLoS Comput Biol 2005, 1: e56. 10.1371/journal.pcbi.0010056

    Article  PubMed Central  PubMed  Google Scholar 

  10. Becq J, Gutierrez MC, Rosas-Magallanes V, Rauzier J, Gicquel B, Neyrolles O, Deschavanne P: Contribution of horizontally acquired genomic islands to the evolution of tubercle bacilli. Mol Biol Evol 2008, 24: 1861–1871. 10.1093/molbev/msm111

    Article  Google Scholar 

  11. Dufraigne C, Fertil B, Lespinats S, Giron A, Deschavanne P: Detection and characterization of horizontal transfers in prokaryotes using genomic signature. Nucleic Acids Res 2005, 33: e6. 10.1093/nar/gni004

    Article  PubMed Central  PubMed  Google Scholar 

  12. Nakamura Y, Itoh T, Matsuda H, Gojobori T: Biased biological functions of horizontally transferred genes in prokaryotic genomes. Nat Genet 2004, 36: 760–766. 10.1038/ng1381

    Article  CAS  PubMed  Google Scholar 

  13. Pride DT, Blaser MJ: Identification of horizontally acquired elements in Helicobacter pylori and other prokaryotes using oligonucleotide difference analysis. Genome Let 2002, 1: 2–15. 10.1166/gl.2002.003

    Article  CAS  Google Scholar 

  14. Abe T, Kanaya S, Kinouchi M, Ichiba Y, Kozuki T, Ikemura T: Informatics for unveiling hidden genome signatures. Genome Res 2003, 13: 693–702. 10.1101/gr.634603

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  15. Teeling H, Meyerdierks A, Bauer M, Amann R, Glöckner FO: Application of tetranucleotide frequencies for the assignment of genomic fragments. Environ Microbiol 2004, 6: 938–947. 10.1111/j.1462-2920.2004.00624.x

    Article  CAS  PubMed  Google Scholar 

  16. Jernigan RW, Baran RH: Pervasive properties of the genomic signature. BMC Genomics 2002, 3: 23. 10.1186/1471-2164-3-23

    Article  PubMed Central  PubMed  Google Scholar 

  17. Koski LB, Morton RA, Golding GB: Codon bias and base composition are poor indicators of horizontally transferred genes. Mol Biol Evol 2001, 18: 404–412.

    Article  CAS  PubMed  Google Scholar 

  18. Wang B: Limitations of compositional approach to identify horizontally transferred genes. J Mol Evol 2001, 53: 244–250. 10.1007/s002390010214

    Article  CAS  PubMed  Google Scholar 

  19. Bohlin J, Skjerve E, Ussery D: Reliability and applications of statistical methods based on oligonucleotide frequencies in bacterial and archaeal genomes. BMC Genomics 2008, 9: 104. 10.1186/1471-2164-9-104

    Article  PubMed Central  PubMed  Google Scholar 

  20. Reva ON, Tümmler B: Global features of sequences of bacterial chromosomes, plasmids and phages revealed by analysis of oligonucleotide usage patterns. BMC Bioinformatics 2004, 5: 90. 10.1186/1471-2105-5-90

    Article  PubMed Central  PubMed  Google Scholar 

  21. Reva ON, Tümmler B: Differentiation of regions with atypical oligonucleotide composition in bacterial genomes. BMC Bioinformatics 2005, 6: 251. 10.1186/1471-2105-6-251

    Article  PubMed Central  PubMed  Google Scholar 

  22. SWGB mirror site at the University of Pretoria in South Africa[http://www.bi.up.ac.za/SeqWord/mhhapplet.php]

  23. SWGB mirror site at Hannover Medical School in Germany[http://genomics1.mh-hannover.de/seqword/genomebrowser/mhhapplet.php]

  24. SWGB mirror site at Penn State University in the USA[http://seqword.bx.psu.edu/mhhapplet.php]

  25. Reva O, Tümmler B: Oligonucleotide usage signatures of the Pseudomonas putida KT2440 genome. In Pseudomonas: Genomics and Molecular Biology. Volume Chapter 3. Edited by: Cornelis P. Caister Academic Press; 2008:43–64.

    Google Scholar 

  26. Weinel C, Nelson KE, Tümmler B: Global features of the Pseudomonas putida KT2440 genome sequence. Environ Microbiol 2002, 4: 809–818. 10.1046/j.1462-2920.2002.00331.x

    Article  CAS  PubMed  Google Scholar 

  27. Lawrence JG, Ochman H: Amelioration of bacterial genomes: rates of change and exchange. J Mol Evol 1997, 44: 383–397. 10.1007/PL00006158

    Article  CAS  PubMed  Google Scholar 

  28. Reva O, Tümmler B: Think big – giant genes in bacteria. Environ Microbiol 2008, 10: 768–777. 10.1111/j.1462-2920.2007.01500.x

    Article  CAS  PubMed  Google Scholar 

  29. McCarren J, Brahamsha B: SwmB, a 1.12-megadalton protein that is required for nonflagellar swimming motility in Synechococcus . J Bacteriol 2007, 189: 1158–1162. 10.1128/JB.01500-06

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  30. Puigbò P, Romeu A, Garcia-Vallvé S: HEG-DB: a database of predicted highly expressed genes in prokaryotic complete genomes under translational selection. Nucleic Acids Res 2008, 36: D524-D527. 10.1093/nar/gkm831

    Article  PubMed Central  PubMed  Google Scholar 

  31. Lawrence JG, Hendrix RW, Casjens S: Where are the pseudogenes in bacterial genomes? Trends Microbiol 2001, 9: 535–540. 10.1016/S0966-842X(01)02198-9

    Article  CAS  PubMed  Google Scholar 

  32. Reva ON, Hallin PF, Willenbrock H, Sicheritz-Ponten T, Tümmler B, Ussery DW: Global features of the Alcanivorax borkumensis SK2 genome. Environ Microbiol 2008, 10: 614–625. 10.1111/j.1462-2920.2007.01483.x

    Article  CAS  PubMed  Google Scholar 

  33. SWGB FTP site[ftp://milliways.bi.up.ac.za/SeqWord/GenomeBrowser/]

  34. Chen XH, Koumoutsi A, Scholz R, Eisenreich A, Schneider K, Heinemeyer I, Morgenstern B, Voss B, Hess WR, Reva O, Junge H, Voigt B, Jungblut PR, Vater J, Süssmuth R, Liesegang H, Strittmatter A, Gottschalk G, Borriss R: Comparative analysis of the complete genome sequence of the plant growth-promoting bacterium Bacillus amyloliquefaciens FZB42. Nat Biotechnol 2007, 25: 1007–1014. 10.1038/nbt1325

    Article  CAS  PubMed  Google Scholar 

  35. Klockgether J, Reva O, Tümmler B: Spread of genomic islands between clinical and environmental isolates. In Prokaryotic Diversity: Mechanisms and Significance. Edited by: Logan NA, Lapping-Scott HM, Oyston PCF. Cambridge: Cambridge University Press; 2006:187–200.

    Chapter  Google Scholar 

  36. TETRA[http://www.megx.net/tetra_new/index.html]

  37. GENESTYLE[http://genstyle.imed.jussieu.fr/]

  38. δρ-WEB and Compare_Islands[http://deltarho.amc.nl/cgi-bin/bin/index.cgi]

  39. The Horizontal Gene Transfer Database (HGT-DB)[http://www.tinet.org/~debb/HGT/]

  40. IslandPath[http://www.pathogenomics.sfu.ca/islandpath]

  41. Artemis[http://www.sanger.ac.uk/Software/Artemis/]

  42. Waack S, Keller O, Asper R, Brodag T, Damm C, Fricke WF, Surovcik K, Meinicke P, Merkl R: Score-based prediction of genomic islands in prokaryotic genomes using hidden Markov models. BMC Bioinformatics 2006, 7: 142. 10.1186/1471-2105-7-142

    Article  PubMed Central  PubMed  Google Scholar 

  43. Heidelberg JF, Eisen JA, Nelson WC, Clayton RA, Gwinn ML, Dodson RJ, Haft DH, Hickey EK, Peterson JD, Umayam L, Gill SR, Nelson KE, Read TD, Tettelin H, Richardson D, Ermolaeva MD, Vamathevan J, Bass S, Qin H, Dragoi I, Sellers P, McDonald L, Utterback T, Fleishmann RD, Nierman WC, White O, Salzberg SL, Smith HO, Colwell RR, Mekalanos JJ, Venter JC, Fraser CM: DNA sequence of both chromosomes of the cholera pathogen Vibrio cholerae . Nature 2000, 406: 477–483. 10.1038/35020000

    Article  CAS  PubMed  Google Scholar 

  44. NCBI, Complete Microbial Genomes[http://www.ncbi.nlm.nih.gov/genomes/lproks.cgi]

Download references

Acknowledgements

This work was funded by the National Bioinformatics Network of South Africa http://www.nbn.ac.za. CFD is supported by the Deutsche Forschungsgemeinschaft grant DFG 653/3. The authors are thankful to Prof. Anton Nekrutenko from Penn State University (Pennsylvania, USA) for hosting the SeqWord applet at http://seqword.bx.psu.edu/.

Author information

Authors and Affiliations

Authors

Corresponding author

Correspondence to Oleg N Reva.

Additional information

Authors' contributions

HG participated in development of the MySQL database and the SWGB web-site, Python programming, applet testing; ASR contributed by Java and PHP programming; CFD participated in development and support of the SWGB web-site, Help files and the applet testing; BT participated in development and testing of the OU statistical algorithms; ONR participated in development and testing of the OU statistical algorithms, supervision of the SeqWord project and Python programming.

Authors’ original submitted files for images

Rights and permissions

Open Access This article is published under license to BioMed Central Ltd. This is an Open Access article is distributed under the terms of the Creative Commons Attribution License ( https://creativecommons.org/licenses/by/2.0 ), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

Reprints and permissions

About this article

Cite this article

Ganesan, H., Rakitianskaia, A.S., Davenport, C.F. et al. The SeqWord Genome Browser: an online tool for the identification and visualization of atypical regions of bacterial genomes through oligonucleotide usage. BMC Bioinformatics 9, 333 (2008). https://doi.org/10.1186/1471-2105-9-333

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: https://doi.org/10.1186/1471-2105-9-333

Keywords