Skip to main content
  • Research article
  • Open access
  • Published:

A quantitative account of genomic island acquisitions in prokaryotes

Abstract

Background

Microbial genomes do not merely evolve through the slow accumulation of mutations, but also, and often more dramatically, by taking up new DNA in a process called horizontal gene transfer. These innovation leaps in the acquisition of new traits can take place via the introgression of single genes, but also through the acquisition of large gene clusters, which are termed Genomic Islands. Since only a small proportion of all the DNA diversity has been sequenced, it can be hard to find the appropriate donors for acquired genes via sequence alignments from databases. In contrast, relative oligonucleotide frequencies represent a remarkably stable genomic signature in prokaryotes, which facilitates compositional comparisons as an alignment-free alternative for phylogenetic relatedness.

In this project, we test whether Genomic Islands identified in individual bacterial genomes have a similar genomic signature, in terms of relative dinucleotide frequencies, and can therefore be expected to originate from a common donor species.

Results

When multiple Genomic Islands are present within a single genome, we find that up to 28% of these are compositionally very similar to each other, indicative of frequent recurring acquisitions from the same donor to the same acceptor.

Conclusions

This represents the first quantitative assessment of common directional transfer events in prokaryotic evolutionary history. We suggest that many of the resident Genomic Islands per prokaryotic genome originated from the same source, which may have implications with respect to their regulatory interactions, and for the elucidation of the common origins of these acquired gene clusters.

Background

The advent of whole genome sequencing has drastically altered our perspective on life's evolutionary history. Bacterial genomes are now known to be largely mosaics made up of horizontally transferred genes [14]. In fact, many bacteria that cause disease, like those that cause plague, meningitis, tetanus or cholera have only become virulent after they acquired virulence genes [58], highlighting the impact of horizontal gene transfer on human health [9]. In many cases, questions remain about the natural reservoir of these acquired genes [10, 11].

These horizontally transferred genes are not necessarily acquired one at a time. Numerous bacterial genomes show clusters of recently acquired genes that are known as Genomic Islands (GIs) [3, 1216]. Even though many GIs have unknown functions, some of these acquired gene clusters are involved in pathogenicity (the Pathogenicity Associated Islands), though several other clustered collective functions are known (metabolic islands, degradation islands et cetera) [17, 18]. We hypothesise that, when several GIs reside in a single genome, it is possible that a single donor has been responsible for multiple gene transfer events to that host.

In a previous study, we found that compositionally similar sequences can be clustered together, and a genomic acquisition account of large acquired gene clusters can be established [19]. Such alignment-free compositional analyses focus on the similarity between two sequences with respect to their relative dinucleotide frequencies. In brief, relative dinucleotide frequencies are known to be typical for a given genome, a genomic signature, and similar between related species. This parameter can be used to assess the similarity in composition between different sequences, for example in binning sequences that are thought to be derived from the same organism [2024]. Compositional analyses have been used before to detect compositionally anomalous genes [25], which could subsequently be identified as putative horizontally acquired genes. In addition, similar comparative analyses have been applied to metagenomic datasets, in which genomic fragments were assigned to their probable host based on their compositional similarities [26, 27]. However, few attempts have been made to compositionally compare clusters of acquired genes, in order to indicate common donors, analogous to assigning donors to sequences from metagenomic libraries.

Here we focus on the acquisition accounts of GIs that are identified in the genomes of a large collection of bacterial and archaeal species [28]. By comparing the compositional similarities of all GIs that reside in the same genome, for a large number of genomes, conservative estimates of the maximum number of compositionally distinct donors can be assessed. This will shed new light onto the evolutionary histories of prokaryotes, and the quantitative dynamics of recurrent horizontal gene transfer events of large gene clusters.

Methods

Analyses were carried out as described previously [19], with a few modifications for scaling up the analyses. In brief, GIs were obtained from IslandViewer [28] at http://www.pathogenomics.sfu.ca/islandviewer/download.php, selecting only those species/genera having genome sizes > 800 kb and the GIs with sizes > 10 kb. Smaller genomes are thought to be mostly devoid of GIs since they often represent intracellular symbionts, whereas the 10 kb cut-off is based on previous publications concerning sizes of GIs [29]. IslandViewer is a computational tool that integrates different genomic island prediction software suits; two sequence composition prediction methods (IslandPick [30], SIGI-HMM [31]) and a comparative GI prediction method (IslandPath-DIMOB [13]). Regions that are identified with IslandViewer are annotated as putative genomic islands, and included in our GI set [28]. This database may not cover all large acquired gene clusters, but does allow for large-scale compositional analyses.

With this collection of GIs, the compositional relatedness of each GI was subsequently compared with its respective genome. This was done by comparing the composition of the GI, with the compositions of all genomic fragments of the same size [25, 32]. Next, all GIs residing in the same genome were compared with each other by calculating the average dinucleotide relative abundance difference, or genomic dissimilarity (δ*) [20], after which we cluster all compositionally similar GIs per genome, based on their genomic dissimilarity values using δρ-Web and Compare_Islands [19, 25], respectively. More information on these methods can be found at the website http://deltarho.amc.nl. In brief, distance matrices of GI comparisons per genome revealed the GIs that have a lower genomic dissimilarity than a conservative threshold sequence to its host genome. With respect to these conservative thresholds of relatedness, we included in each comparison a chromosomal fragment of 15 kb with a very low relative dissimilarity with its genome. The relative dissimilarity signifies the dissimilarity between a query sequence and the rest of the genome; the dissimilarity of the query is relative to the collection of non-overlapping genomic fragments of identical size as the query. The relative dissimilarity is expressed as a percentage of genomic fragments with a lower genomic dissimilarity than the query sequence. A relative dissimilarity of 95% signifies that 95% of all non-overlapping genomic fragments of identical size as the query is more similar to the genome than the query is. The threshold sequences are based on progressively lower relative dissimilarity values. Core Islands CI-25, CI-10, CI-5 and CI-0 represent the four threshold sequences with relative dissimilarity values of 25%, 10%, 5% and 0%. In other words, a Core Island from a specific genome is compositionally very similar to its host. In order to test that different GIs from the same genome originate from a same donor species, the GIs need to be compositionally more similar to each other than the Core Island is to its host genome. Thus, if GIs meet these similarity thresholds, we score these GIs as clusters that have a compositionally similar background, and therefore likely a common origin. In some instances, compositionally similar GIs are not clustered together due to a high compositional similarity between a GI from a predicted cluster with a GI outside of that specific cluster. The similarity threshold between the unclustered GI and one of the other GIs in the cluster is not met. These problems in an unambiguous interpretation of the clustering are categorized as 'conflicts', and subsequently all GIs from that genome are excluded from the cluster analyses in order to reduce potential misclassifications. An example of an analysis with a clustering conflict is given in Additional File 1, which gives the compositional distance matrix of six GIs from the genome of Clostridium botulinum Ba4_657 (NC_012658), relative to the CI-25 threshold sequence. These conflicts are removed from the analyses in an attempt to obtain a conservative dataset with few ambiguities.

A set of stand-alone scripts is available from the authors (at https://trac.nbic.nl/brsp200901_vanpassel/wiki), with both instructions on how to perform the analyses for GI sets automatically as well as the raw data for the analyses presented here. It iterates all calculations for each applicable GI-host and GI-GI combination within a host, allowing the user to choose different cut-off values of compositional dissimilarity (i.e., the Core Islands), as well as GI size. For individual GI/genome comparisons, Compare_Islands can be used at http://deltarho.amc.nl[19].

In order to test to test the accuracy of our clustering approach and cut-offs, we simulate a clustering fidelity by analyzing how frequently fragments from the same genome are clustered together when a pool of phylogenetically unrelated 15 kb sequences are compared. For each of the four thresholds, 100 analyses are carried out, each consisting of a set of 100 sequences; 90 originating from distinct genera, and ten non-biological randomized sequences. The 90 sequences originate from 30 genomes, with three fragments per genome, and two out of these three have a relative dissimilarity like the threshold that is being simulated. The third sequence has a relative dissimilarity of 50% with its host genome. In this simulation, the accuracy of the clustering is expressed as a percentage, which indicates how often the threshold sequences are clustered with sequences from the same host genome.

Results

First, we extracted all Genomic Islands from IslandViewer (December 2009, [28]), amounting to a total of 5447 sequences between 2.2 and 143 kb in size, originating from 339 distinct genomes. After applying the conservative criteria discussed in the Material and Methods section (genome size > 800 kb, GI size > 10 kb, monochromosomal genomes, no internal conflicts in the clustering approach using the CI-10 cut-off threshold sequence), we maintained 1787 GIs (33%) that vary in size from 10 kb to 130 kb (average 20.7 kb, Figure 1). These GIs are present in 246 genome sequences, which represent 88 species in 45 genera (Additional File 2). This means on average 7 GIs per genome, varying from 1 (in 17 genomes) to 27 (in Xanthomonas oryzae MAFF 311018) GIs per genome (Figure 2). These 1787 GIs, using the compositional threshold of CI-10, were used in our subsequent investigations, unless noted otherwise.

Figure 1
figure 1

Size distribution of 1787 Genomic Islands > 10 kb in 246 genome sequences (note the logarithmic scale on the vertical axis). The GIs are binned per 2 kb in size.

Figure 2
figure 2

Number of Genomic Islands per genome for the 246 genomes tested (Genome size > 800 kb, with GIs size > 10 kb and no conflicts).

Second, we analyzed the composition dissimilarities of these GIs with their respective host chromosomes similar to previous analyses on comparisons between plasmids and host chromosomes [33]. Of the 1787 GIs, 1394 (78%) are compositionally anomalous compared to their host genome (with a genomic dissimilarity score higher than that of 90% of the genomic fragments of equal length, Figure 3). Of these 1394 GIs, a large number (683 GIs, 49%) have a very low GC content compared to fragments of identical length from their respective host genomes (i.e., lower than 95% of identical sized fragments from their respective genome). Out of the 1787 GIs, only 11 GIs (0.6%) have a lower genomic dissimilarity with the host genome than with the threshold sequence CI-10, meaning that these GIs are compositionally extremely similar to their respective host genomes (Additional File 2).

Figure 3
figure 3

Distribution of the relative compositional similarity and GC similarity of all GIs (1787) with their respective genomes, with 1395 (78%, in red) of the GIs having a relative dissimilarity of 90%.

Third, we tested all GIs per genome for their compositional dissimilarity. In 86 genomes (of 44 species in 28 genera) we identify 134 clusters, including a total of 271 GIs (15.3% of the 1770 GIs that reside in genomes with at least 2 GIs, Figure 4, Additional File 3). The number of clustered GIs per genome varies between the minimal 2 GIs in a single cluster (in 56 genomes), to 15 GIs in a total of 7 clusters (in the EHEC strain E. coli O157H7 Sakai, Figure 5). The only clusters that contain three GIs occur in three genomes Bradyrhizobium ORS278, Escherichia coli O157H7 strain Sakai and Rhodobacter sphaeroides ATCC 17025. In Bradyrhizobium ORS278, the three GIs that are clustered together show a high compositional similarity to the host genome sequence. Two of these in fact belong to the 11 GIs that are compositionally very similar to their respective genomes, and therefore these two GIs are unlikely to represent horizontal transfer events. For E. coli O157H7 strain Sakai and R. sphaeroides ATCC 17025, we tested whether the three clustered GIs are more similar than a set of five sequences belong to the 10% most compositionally similar sequences of the genome (Tables 1 and 2, respectively). We find that for both E. coli O157H7 strain Sakai and R. sphaeroides ATCC 17025, the GIs that are clustered together in threes, are on average equally or more similar to each other than the five sequences that represent the conservative genome signature (Tables 1 and 2).

Figure 4
figure 4

Number of clustered GIs per genome.

Figure 5
figure 5

Clustering of the 24 Genomic Islands > 10 kb in Escherichia coli O157H7 strain Sakai in seven clusters and nine singletons. Below the cut-off value (red line; dissimilarity < 1.44, see Additional File 2), seven clusters are identified (six clusters with two GIs, and one with three GIs), with a total of 15 Genomic Islands (indicated with seven colored bars). The GIs and their numbers are identified in Additional File 5.

Table 1 Compositional comparison of Core Islands e1-e5 (with relative dissimilarities of 10%) of Escherichia coli O157H7 with each other (underlined), and Genomic Islands with each other (bold)
Table 2 Compositional comparison of Core Islands r1-r5 (with relative dissimilarities of 10%) of Rhodobacter sphaeroides with each other (underlined), and Genomic Islands with each other (bold)

In genomes with multiple chromosomes, we test for GIs that are compositionally very similar to each other, yet reside on a different replicon. Out of 110 GIs identified in this set of genomes, 38 GIs are assigned to a total of 19 clusters (Additional File 4). Out of these 38 GIs, ten (26%), all of them in Burkholderia genomes, are not located on the same chromosome.

When relaxing the similarity threshold by using the genome signature difference between the Core Island 25 (CI-25) and the genome, we observe only 1370 GIs in a total of 220 genome sequences that meet our criteria. With this more lenient threshold, a total of 16 GIs are now compositionally more similar to the host genome than the cut-off sequence CI-25 is to the genome. A total of 383 GIs (out of 1353 GIs that reside with at least one other GI in a genome; 28%) are now grouped together in 185 clusters, with 13 clusters containing three GIs.

In contrast, when making the composition similarity threshold substantially more conservative (i.e., using CI-5), we find 2047 GIs in a total of 260 genome sequences. Only 9 GIs are now compositionally more similar to the host genome when compared to the cut-off sequence CI-5. Still, there are 99 clusters containing 202 GIs (out of 2030 GIs that reside with at least one other GI in a genome; 10%), which show very high compositional similarity within each cluster (Table 3). Finally, using threshold CI-0, in which GIs need to be compositionally more similar to each other than the genomic fragment that has a nearly identical dinucleotide composition as the host genome, we find only 40 GIs that form 20 clusters. These 40 GIs represent only 1.8% of the total number of GIs included in this analysis.

Table 3 Overview of the characteristics of the GI analyses using decreasing similarity thresholds (for all GIs > 10 kb)

Finally, in order to validate the accuracy of the clustered GIs per threshold, we simulated cluster assignments by comparing sets of 100 15 kb fragments from random prokaryotic genomes for their compositionally most similar fragment. Each set of 100 fragments consists of three fragments per species, for 30 species of distinct genera, and includes 10 random synthetic sequences with no biological significance. For each threshold, the simulation was carried out 100 times. Accuracy is expressed as the percentage in which a sequence is found to be most similar to another sequence from the same genome. These values are used as proxies for the correct assignment of a sequence from the same genome, and range from 99.9% accurate for the strictest threshold of CI-0, to 94.8% accurate for the CI-25 threshold (Table 3).

Discussion

By comparing Genomic Islands from 339 bacterial and archaeal chromosomes, we explore the dynamics of the genome-specific acquisition accounts on a large scale. These analyses show us that in numerous cases, distinct GIs in particular genomes are remarkably similar in composition. This leads us to speculate that, using a conservative similarity cut-off, in 15.3% of the cases, multiple acquisition events of GIs have taken place from a donor with a very similar base composition as the acceptor.

For this goal, we developed a suite of scripts that allows users to customize these analyses by modifying the minimal GI length, or the similarity cut-off sequence (i.e., the genomic Core Island of each genome with a certain genomic compositional dissimilarity). By increasing the stringency to the very conservative CI-5, which means that GIs need to be more similar to each other in composition than 95% of the genomic fragments, we still find that 10% of the tested GIs can be grouped together to a total of 99 GI clusters. However, we investigate the accuracy of our assessments by simulating the clustering efficiency in a randomized sample of sequences with different thresholds. We find that the prediction accuracy according to this simulation is > 94,8% even for the least conservative compositional threshold (CI-25). This gives credibility to our findings of substantial recurrent transfer events from the same donor to the same host. In Bradyrhizobium ORS278, we find a cluster containing three GIs, two of which cannot be considered compositionally dissimilar from the genome. In this case, these GIs may have been residing for a substantial amount of time in the genome, and have ameliorated to the host's genome composition [34, 35].

This approach does not discriminate between separate introgressions of multiple compositionally very similar sequences from a common donor, and the post-acquisition intragenomic dispersal of a large Genomic Island. Technically, this is of little importance, since in both cases the host of the distinct GIs would be a similar donor. With respect to compositionally similar GIs that reside on separate chromosomes, we find that that incoming GIs seem to be indiscriminate between the replicon it integrates in, or that subsequent dispersal throughout the genome can result in a move to a different replicon.

The association of GIs with virulence factors [36] emphasize the significant role of acquired gene clusters in the evolution of numerous pathogens. Investigations into the repetitive acquisition of GIs from a common source may help identifying potential donors of these sequences through for example the association with species-specific sequence motifs such as DNA uptake sequences [37]. Also, a common origin of compositionally similar clusters may result in common regulatory modules, interactions or mobilizing capacities. For example, a study into small regulatory RNAs (sRNAs) on Genomic Islands in Salmonella typhimurium revealed that sRNAs mainly affect the expression of flanking genes [38]. If Genomic Islands disperse throughout the genome, our analysis would facilitate the identification of potential associated regulatory targets that are no longer adjacent.

Unfortunately, the forces that shape the genome signatures of prokaryotes are still unknown. It has been speculated that they could include species-specific properties such as DNA modifications, replication and repair mechanisms [20], though recently statistical support has been found for an environmental influence on the oligonucleotide compositions [39], which could mean that a similar environment could also cause similarities in genome signature. For compositional comparisons such as described here, it is of interest to understand what conditions shape the composition of DNA to which extent, in order to pinpoint potential pitfalls in grouping Genomic Islands.

Conclusions

Even when lacking sequence alignments, numerous large acquired gene clusters in sequenced genomes can be associated with each other individually via substantial compositional similarities. Our analysis suggests, backed up by simulations, that in many cases recurring horizontal gene transfer events have taken place between a donor and acceptor organism. These analyses do not only quantify these events, but also enable further investigations into the origin of these Genomic Islands, and even help analyzing possible interactions between related sequences.

References

  1. Ochman H, Jones IB: Evolutionary dynamics of full genome content in Escherichia coli. Embo J. 2000, 19 (24): 6637-6643. 10.1093/emboj/19.24.6637.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  2. Ochman H, Lawrence JG, Groisman EA: Lateral gene transfer and the nature of bacterial innovation. Nature. 2000, 405 (6784): 299-304. 10.1038/35012500.

    Article  CAS  PubMed  Google Scholar 

  3. Coleman ML, Sullivan MB, Martiny AC, Steglich C, Barry K, Delong EF, Chisholm SW: Genomic islands and the ecology and evolution of Prochlorococcus. Science. 2006, 311 (5768): 1768-1770. 10.1126/science.1122050.

    Article  CAS  PubMed  Google Scholar 

  4. van Passel MW, Marri PR, Ochman H: The emergence and fate of horizontally acquired genes in Escherichia coli. PLoS Comput Biol. 2008, 4 (4): e1000059-10.1371/journal.pcbi.1000059.

    Article  PubMed  PubMed Central  Google Scholar 

  5. Lathem WW, Price PA, Miller VL, Goldman WE: A plasminogen-activating protease specifically controls the development of primary pneumonic plague. Science. 2007, 315 (5811): 509-513. 10.1126/science.1137195.

    Article  CAS  PubMed  Google Scholar 

  6. Bille E, Zahar JR, Perrin A, Morelle S, Kriz P, Jolley KA, Maiden MC, Dervin C, Nassif X, Tinsley CR: A chromosomally integrated bacteriophage in invasive meningococci. J Exp Med. 2005, 201 (12): 1905-1913. 10.1084/jem.20050112.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  7. Bruggemann H, Baumer S, Fricke WF, Wiezer A, Liesegang H, Decker I, Herzberg C, Martinez-Arias R, Merkl R, Henne A, et al: The genome sequence of Clostridium tetani, the causative agent of tetanus disease. Proc Natl Acad Sci USA. 2003, 100 (3): 1316-1321. 10.1073/pnas.0335853100.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  8. Waldor MK, Mekalanos JJ: Lysogenic conversion by a filamentous phage encoding cholera toxin. Science. 1996, 272 (5270): 1910-1914. 10.1126/science.272.5270.1910.

    Article  CAS  PubMed  Google Scholar 

  9. Pallen MJ, Wren BW: Bacterial pathogenomics. Nature. 2007, 449 (7164): 835-842. 10.1038/nature06248.

    Article  CAS  PubMed  Google Scholar 

  10. Hsiao WW, Ung K, Aeschliman D, Bryan J, Finlay BB, Brinkman FS: Evidence of a large novel gene pool associated with prokaryotic genomic islands. PLoS Genet. 2005, 1 (5): e62-10.1371/journal.pgen.0010062.

    Article  PubMed  PubMed Central  Google Scholar 

  11. Cortez D, Forterre P, Gribaldo S: A hidden reservoir of integrative elements is the major source of recently acquired foreign genes and ORFans in archaeal and bacterial genomes. Genome Biol. 2009, 10 (6): R65-10.1186/gb-2009-10-6-r65.

    Article  PubMed  PubMed Central  Google Scholar 

  12. Hacker J, Kaper JB: Pathogenicity islands and the evolution of microbes. Annu Rev Microbiol. 2000, 54: 641-679. 10.1146/annurev.micro.54.1.641.

    Article  CAS  PubMed  Google Scholar 

  13. Hsiao W, Wan I, Jones SJ, Brinkman FS: IslandPath: aiding detection of genomic islands in prokaryotes. Bioinformatics. 2003, 19 (3): 418-420. 10.1093/bioinformatics/btg004.

    Article  CAS  PubMed  Google Scholar 

  14. Mantri Y, Williams KP: Islander: a database of integrative islands in prokaryotic genomes, the associated integrases and their DNA site specificities. Nucleic Acids Res. 2004, D55-58. 32 Database

  15. Merkl R: SIGI: score-based identification of genomic islands. BMC Bioinformatics. 2004, 5 (1): 22.-10.1186/1471-2105-5-22.

    Article  PubMed  PubMed Central  Google Scholar 

  16. Zhang R, Zhang CT: A systematic method to identify genomic islands and its applications in analyzing the genomes of Corynebacterium glutamicum and Vibrio vulnificus CMCP6 chromosome I. Bioinformatics. 2004, 20 (5): 612-622. 10.1093/bioinformatics/btg453.

    Article  CAS  PubMed  Google Scholar 

  17. Dobrindt U, Hochhut B, Hentschel U, Hacker J: Genomic islands in pathogenic and environmental microorganisms. Nat Rev Microbiol. 2004, 2 (5): 414-424. 10.1038/nrmicro884.

    Article  CAS  PubMed  Google Scholar 

  18. Langille MG, Hsiao WW, Brinkman FS: Detecting genomic islands using bioinformatics approaches. Nat Rev Microbiol. 2010, 8 (5): 373-382. 10.1038/nrmicro2350.

    Article  CAS  PubMed  Google Scholar 

  19. van Passel MW, Bart A, Thygesen HH, Luyf AC, van Kampen AH, van der Ende A: An acquisition account of genomic islands based on genome signature comparisons. BMC Genomics. 2005, 6 (1): 163-10.1186/1471-2164-6-163.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  20. 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 

  21. Karlin S: Detecting anomalous gene clusters and pathogenicity islands in diverse bacterial genomes. Trends Microbiol. 2001, 9 (7): 335-343. 10.1016/S0966-842X(01)02079-0.

    Article  CAS  PubMed  Google Scholar 

  22. Sandberg R, Branden CI, Ernberg I, Coster J: Quantifying the species-specificity in genomic signatures, synonymous codon choice, amino acid usage and G+C content. Gene. 2003, 311: 35-42.

    Article  CAS  PubMed  Google Scholar 

  23. van Passel MW, Kuramae EE, Luyf AC, Bart A, Boekhout T: The reach of the genome signature in prokaryotes. BMC Evol Biol. 2006, 6: 84-10.1186/1471-2148-6-84.

    Article  PubMed  PubMed Central  Google Scholar 

  24. 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 

  25. van Passel MW, Luyf AC, van Kampen AH, Bart A, van der Ende A: {delta}{rho}-Web, an online tool to assess composition similarity of individual nucleic acid sequences. Bioinformatics. 2005, 21 (13): 3053-3055. 10.1093/bioinformatics/bti460.

    Article  CAS  PubMed  Google Scholar 

  26. Teeling H, Waldmann J, Lombardot T, Bauer M, Glockner FO: TETRA: a web-service and a stand-alone program for the analysis and comparison of tetranucleotide usage patterns in DNA sequences. BMC Bioinformatics. 2004, 5: 163-10.1186/1471-2105-5-163.

    Article  PubMed  PubMed Central  Google Scholar 

  27. Teeling H, Meyerdierks A, Bauer M, Amann R, Glockner FO: Application of tetranucleotide frequencies for the assignment of genomic fragments. Environ Microbiol. 2004, 6 (9): 938-947. 10.1111/j.1462-2920.2004.00624.x.

    Article  CAS  PubMed  Google Scholar 

  28. Langille MG, Brinkman FS: IslandViewer: an integrated interface for computational identification and visualization of genomic islands. Bioinformatics. 2009, 25 (5): 664-665. 10.1093/bioinformatics/btp030.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  29. Hentschel U, Hacker J: Pathogenicity islands: the tip of the iceberg. Microbes Infect. 2001, 3 (7): 545-548. 10.1016/S1286-4579(01)01410-1.

    Article  CAS  PubMed  Google Scholar 

  30. Langille MG, Hsiao WW, Brinkman FS: Evaluation of genomic island predictors using a comparative genomics approach. BMC Bioinformatics. 2008, 9: 329-10.1186/1471-2105-9-329.

    Article  PubMed  PubMed Central  Google Scholar 

  31. 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  PubMed Central  Google Scholar 

  32. van Passel MW, Bart A, Waaijer RJ, Luyf AC, van Kampen AH, van der Ende A: An in vitro strategy for the selective isolation of anomalous DNA from prokaryotic genomes. Nucleic Acids Res. 2004, 32 (14): e114.-

    Article  PubMed  PubMed Central  Google Scholar 

  33. van Passel MW, Bart A, Luyf AC, van Kampen AH, van der Ende A: Compositional discordance between prokaryotic plasmids and host chromosomes. BMC Genomics. 2006, 7 (1): 26.-10.1186/1471-2164-7-26.

    Article  PubMed  PubMed Central  Google Scholar 

  34. Lawrence JG, Ochman H: Amelioration of bacterial genomes: rates of change and exchange. J Mol Evol. 1997, 44 (4): 383-397. 10.1007/PL00006158.

    Article  CAS  PubMed  Google Scholar 

  35. Marri PR, Golding GB: Gene amelioration demonstrated: the journey of nascent genes in bacteria. Genome. 2008, 51 (2): 164-168. 10.1139/G07-105.

    Article  CAS  PubMed  Google Scholar 

  36. Ho Sui SJ, Fedynak A, Hsiao WW, Langille MG, Brinkman FS: The association of virulence factors with genomic islands. PLoS One. 2009, 4 (12): e8094-10.1371/journal.pone.0008094.

    Article  PubMed  PubMed Central  Google Scholar 

  37. Sandberg R, Winberg G, Branden CI, Kaske A, Ernberg I, Coster J: Capturing whole-genome characteristics in short sequences using a naive Bayesian classifier. Genome Res. 2001, 11 (8): 1404-1409. 10.1101/gr.186401.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  38. Padalon-Brauch G, Hershberg R, Elgrably-Weiss M, Baruch K, Rosenshine I, Margalit H, Altuvia S: Small RNAs encoded within genetic islands of Salmonella typhimurium show host-induced expression and role in virulence. Nucleic Acids Res. 2008, 36 (6): 1913-1927. 10.1093/nar/gkn050.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  39. Bohlin J, Skjerve E: Examination of genome homogeneity in prokaryotes using genomic signatures. PLoS One. 2009, 4 (12): e8113-10.1371/journal.pone.0008113.

    Article  PubMed  PubMed Central  Google Scholar 

Download references

Acknowledgements and funding

MWJvP is funded by the Netherlands Organization for Scientific Research (NWO) via a VENI grant. The authors would like to thank Angela Luyf for her assistance in customizing the Delta-Rho Web and Compare_Islands web applications, and Dr. Morris Swertz for critically reading the manuscript.

Author information

Authors and Affiliations

Authors

Corresponding author

Correspondence to Mark WJ van Passel.

Additional information

Authors' contributions

MWJvP conceived the study, participated in the design of the software, analyzed the data and wrote the manuscript. TER designed the software, analyzed the data and helped draft the manuscript. Both authors read and approved the final manuscript.

Electronic supplementary material

12864_2011_3586_MOESM1_ESM.DOCX

Additional file 1: Example of a GI clustering conflict in Clostridium botulinum Ba4 657. Example of a GI clustering conflict. In Clostridium botulinum Ba4 657, six GIs larger than 10 kb are identified by IslandViewer. The CI-25 threshold sequence has a genomic dissimilarity to its genome of 30,86. GI-1 and GI-5 are compositionally more similar to each other (δ* of 22,7), as are GI-2 and GI-5 (δ* of 27,7). However, GI-1 and GI-2 are much more dissimilar (δ* of 42,1), and therefore could be considered as a clustering conflict. (DOCX 78 KB)

12864_2011_3586_MOESM2_ESM.DOCX

Additional file 2: Complete table for all genomes that contain GIs. Complete table for all genomes (> 800 kb) that contain GIs ( > 10 kb), and have no conflicts in the genome. In green are highlighted the cases where clusters of three GIs are found. The dissimilarity cut-off is expressed in the genomic dissimilarity value between the Core Island (in this case, CI-10) and the genome sequence. (DOCX 121 KB)

12864_2011_3586_MOESM3_ESM.DOCX

Additional file 3: List of genomes with clustered GIs. List of genomes with clusters in which more than 1 GI are located, using the compositional threshold of CI-10. (DOCX 90 KB)

12864_2011_3586_MOESM4_ESM.DOCX

Additional file 4: Clustered GIs in multichromosomal genomes. Clusters of GIs in genomes with multiple chromosomes using cut-off of CI-10. Highlighted in green are clustered GIs that are located on a different replicon. (DOCX 59 KB)

12864_2011_3586_MOESM5_ESM.XLSX

Additional file 5: Characteristics of all Genomic Islands in this study. Characteristics of all Genomic Islands (5447) analyzed in this study, including their number in the genome of occurrence. The 1787 GIs that comply with the criteria (Genome > 800 kb, GI size > 10 kb, no conflicts in the clustering analysis and using cut-off threshold CI-10) are separated from the rest of the GIs by a blank line. (XLSX 720 KB)

Authors’ original submitted files for images

Rights and permissions

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.

Reprints and permissions

About this article

Cite this article

Roos, T.E., van Passel, M.W. A quantitative account of genomic island acquisitions in prokaryotes. BMC Genomics 12, 427 (2011). https://doi.org/10.1186/1471-2164-12-427

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: https://doi.org/10.1186/1471-2164-12-427

Keywords