Score-based prediction of genomic islands in prokaryotic genomes using hidden Markov models.
BACKGROUND
Horizontal gene transfer (HGT) is considered a strong evolutionary force shaping the content of microbial genomes in a substantial manner. It is the difference in speed enabling the rapid adaptation to changing environmental demands that distinguishes HGT from gene genesis, duplications or mutations. For a precise characterization, algorithms are needed that identify transfer events with high reliability. Frequently, the transferred pieces of DNA have a considerable length, comprise several genes and are called genomic islands (GIs) or more specifically pathogenicity or symbiotic islands.


RESULTS
We have implemented the program SIGI-HMM that predicts GIs and the putative donor of each individual alien gene. It is based on the analysis of codon usage (CU) of each individual gene of a genome under study. CU of each gene is compared against a carefully selected set of CU tables representing microbial donors or highly expressed genes. Multiple tests are used to identify putatively alien genes, to predict putative donors and to mask putatively highly expressed genes. Thus, we determine the states and emission probabilities of an inhomogeneous hidden Markov model working on gene level. For the transition probabilities, we draw upon classical test theory with the intention of integrating a sensitivity controller in a consistent manner. SIGI-HMM was written in JAVA and is publicly available. It accepts as input any file created according to the EMBL-format.It generates output in the common GFF format readable for genome browsers. Benchmark tests showed that the output of SIGI-HMM is in agreement with known findings. Its predictions were both consistent with annotated GIs and with predictions generated by different methods.


CONCLUSIONS
SIGI-HMM is a sensitive tool for the identification of GIs in microbial genomes. It allows to interactively analyze genomes in detail and to generate or to test hypotheses about the origin of acquired genes.
Waack
Stephan
S
Institut für Informatik, Universität Göttingen, Lotzestr, 16-18, 37083 Göttingen, Germany.
Keller
Oliver
O
Asper
Roman
R
Brodag
Thomas
T
Damm
Carsten
C
Fricke
Wolfgang Florian
WF
Surovcik
Katharina
K
Meinicke
Peter
P
Merkl
Rainer
R
D000465:Algorithms; D001412:Bacillus subtilis / Q000235:genetics; D001482:Base Composition; D003062:Codon; D019295:Computational Biology; D004269:DNA, Bacterial / Q000235:genetics*; D048168:Escherichia coli K12 / Q000235:genetics; D022761:Gene Transfer, Horizontal; D016680:Genome, Bacterial; D044404:Genomic Islands; D008390:Markov Chains; D008957:Models, Genetic; D012984:Software; D014734:Vibrio cholerae / Q000235:genetics
D003062:Codon; D004269:DNA, Bacterial