Abstract
The degeneracy of codons allows a multitude of possible sequences to code for the same protein. Hidden within the particular choice of sequence for each organism are over 100 previously undiscovered biologically significant, short oligonucleotides (length, 2 to 7 nucleotides). We present an information-theoretic algorithm that finds these novel signals. Applying this algorithm to the 209 sequenced bacterial genomes in the NCBI database, we determine a set of oligonucleotides for each bacterium which uniquely characterizes the organism. Some of these signals have known biological functions, like restriction enzyme binding sites, but most are new. An accompanying scoring algorithm is introduced that accurately (92%) places sequences of 100 kb with their correct species among the choice of hundreds. This algorithm also does far better than previous methods at relating phage genomes to their bacterial hosts, suggesting that the lists of oligonucleotides are "genomic fingerprints" that encode information about the effects of the cellular environment on DNA sequence. Our approach provides a novel basis for phylogeny and is potentially ideally suited for classifying the short DNA fragments obtained by environmental shotgun sequencing. The methods developed here can be readily extended to other problems in bioinformatics.
MeSH Terms
Algorithms
Bacteria/classification,genetics
Bacteriophages/genetics
Computational Biology
DNA Fingerprinting
DNA, Bacterial/genetics
DNA, Viral/genetics
Databases, Genetic
Genome, Bacterial
Genome, Viral
Genomics
Oligodeoxyribonucleotides/genetics
Phylogeny
Sequence Homology, Nucleic Acid
Chemicals
DNA, Bacterial
DNA, Viral
Oligodeoxyribonucleotides
Authors & Affiliations
4 authors, click to expand affiliations / ORCID
Robins Harlan
Institute for Advanced Study, Natural Sciences, Einstein Drive, Princeton, NJ 08540, USA.
[email protected]
Krasnitz Michael
Barak Hagar
Levine Arnold J
References (11)
11 references, click to expand
-
Genome-wide analysis of mRNA decay in resting and activated primary human T lymphocytes.
Nucleic Acids Res. 2002 Dec 15;30(24):5529-38
PMID: 12490721
-
Decay rates of human mRNAs: correlation with functional characteristics and sequence attributes.
Genome Res. 2003 Aug;13(8):1863-72
PMID: 12902380
-
Isochores and tissue-specificity.
Nucleic Acids Res. 2003 Sep 1;31(17):5212-20
PMID: 12930973
-
The relationship between palindrome avoidance and intragenic codon usage variations: a Monte Carlo study.
Biochem Biophys Res Commun. 2004 Apr 9;316(3):755-62
PMID: 15033465
-
Environmental genome shotgun sequencing of the Sargasso Sea.
Science. 2004 Apr 2;304(5667):66-74
PMID: 15001713
-
Similarities and dissimilarities of phage genomes.
Proc Natl Acad Sci U S A. 1996 Jun 11;93(12):5854-9
PMID: 8650182
-
Maximum entropy modeling of short sequence motifs with applications to RNA splicing signals.
J Comput Biol. 2004;11(2-3):377-94
PMID: 15285897
-
Tissue-specific codon usage and the expression of human genes.
Proc Natl Acad Sci U S A. 2004 Aug 24;101(34):12588-91
PMID: 15314228
-
Correlation between the abundance of Escherichia coli transfer RNAs and the occurrence of the respective codons in its protein genes.
J Mol Biol. 1981 Feb 15;146(1):1-21
PMID: 6167728
-
Computational DNA sequence analysis.
Annu Rev Microbiol. 1994;48:619-54
PMID: 7826021
-
RESCUE-ESE identifies candidate exonic splicing enhancers in vertebrate exons.
Nucleic Acids Res. 2004 Jul 1;32(Web Server issue):W187-90
PMID: 15215377