- Open Access
CBDB: The codon bias database
BMC Bioinformatics volume 13, Article number: 62 (2012)
In many genomes, a clear preference in the usage of particular codons exists. The mechanisms that induce codon biases remain an open question; studies have attributed codon usage to translational selection, mutational bias and drift. Furthermore, correlations between codon usage within host genomes and their viral pathogens have been observed for a myriad of host-virus systems. As such, numerous studies have investigated codon usage and codon bias in an effort to better understand how species evolve. Numerous metrics have been developed to identify biases in codon usage. In addition, a few data repositories of codon bias data are available, differing in the metrics reported as well as the number and taxonomy of strains examined.
We have created a new web resource called the Codon Bias Database (CBDB) which provides information regarding the codon bias within the set of highly expressed genes for 300+ bacterial genomes. CBDB was developed to provide a resource for researchers investigating codon bias in bacteria, facilitating comparisons between strains and species. Furthermore, the site was created to serve those studying adaptation in phage; the genera selected for this first release of CBDB all have sequenced, annotated bacteriophages. The annotations and sequences for the highly expressed gene set are available for each strain in addition to the strain’s codon bias measurements.
Comparing species and strains provides a comprehensive look at how codon usage has been shaped over evolutionary time and can elucidate the putative mechanisms behind it. The Codon Bias Database provides a centralized repository of look-up tables and codon usage bias measures for a wide variety of genera, species and strains. Through our analysis of the variation in codon usage within the strains presently available, we find that most members of a genus have a codon composition most similar to other members of its genus, although not necessarily other members of its species.
Why does codon bias exist?
The redundancy within the genetic code accommodates a variable number of codons to encode for the same amino acid. Codon usage biases have been found to exist, ranging from relatively neutral to extremely strong [1–6]. Debate within the scientific community continues as to exactly why codon biases exist. Theories based on translational selection, mutational biases, and drift have all been found to contribute to codon biases [7–11]. Sequence-based analysis revealed that within organisms having a biased genome, the most frequently occurring codons often reflect the most abundant transfer RNA (tRNA) available [12–15]. Furthermore, through direct molecular manipulation protein-throughput was increased by designing coding sequences to utilize the most abundant tRNAs [16–18]. These bioinformatic and experimental studies suggest that translational selection may be the primary factor shaping codon usage.
Correspondences between codon usage and tRNA availability may signify selection for translational accuracy, selection for translational efficiency, or both [19–30]. Theories for selection for translational accuracy assume that the codon with the highest tRNA abundance has a lower missense error rate than its synonymous codons, considering tRNA gene copy numbers of both cognate and near-cognate tRNA abundances (e.g. [11, 31]). If codon bias exists as a result of selection to maximize protein-throughput, codon usage often reflects available tRNAs. Because codon bias is often strongest within highly expressed genes and there exists a correspondence between generation time and the strength of codon bias within these highly expressed genes , it is believed that selection for translational efficiency plays some role in determining a species’ codon usage. Moreover, recent studies have shown that the usage of particular synonymous codons can also impact protein folding or misfolding [32–34].
As the abundances of individual tRNAs vary from one species to the next, so too do the preferences in codons. Translational selection seems most prevalent in viral species. Correspondence in virus-host codon usage has been observed in DNA-based, RNA-based, and retro-transcribing viruses (e.g. [35–39]). This observation is not isolated to eukaryote-infecting viruses; bacteriophages also exhibit codon biases similar to their host species (e.g. [40–42]). Given that viruses are heavily dependent upon their host for biosynthesis, utilizing the most prevalent tRNAs within the host would likely give the virus a translational advantage.
Examining codon bias
In order to quantify the bias present within a species, a number of metrics have been developed. Calculating the codon usage frequencies within the genome can reveal biases, in particular when comparing these usage profiles to expected usage patterns and/or the usage profiles of synonymous codons, e.g. the frequency of preferred codons (FOP) metric , the synonymous codon usage order (SCUO) metric , and the %MinMax algorithm . In order to compare the strength of biases within and between species, the relative synonymous codon usage (RSCU) and the geometric mean of the RSCU values, the codon adaptation index (CAI), were developed in which codon usage frequencies within highly expressed genes (HEGs) are specifically examined [6, 23, 45]. Numerous extensions of the CAI metric have been proposed, e.g. the self-consistent codon index (SCCI)  and relative codon adaption index (rCAI) . Rather than looking at the codons themselves to ascertain biases, a second approach exists in which biases are assessed relative to individual tRNA abundances. The tRNA adaptation index (tAI), also inspired by the CAI metric, considers not only the gene copy number of the tRNA with the perfectly matched anticodon but also those tRNAs which can bind imperfectly . The “local tAI” metric takes a similar approach, however the tAI measure is averaged for sliding windows across the gene rather than for the whole gene sequence . In many species the codon usage in HEGs matches tRNA abundance [12, 48].
Given the information encoded within the usage of codons, a number of tools and databases have been developed for analyzing the codon and tRNA content of genic and genomic sequences. CAIcal , CodonExplorer , and CodonW  calculate CAI values for user input sequences while E-CAI  calculates the expected CAI values by generating random sequences with a G+C content and amino acid composition similar to the user input sequence. The application CodonO  can analyze individual genomes or compare genomes using the SCUO metric. JCat  and GCUA  both calculate the RSCU values of user input sequences relative to a reference organism’s usage profile. The RSCU values for many species are contained in the Codon Usage Database ; this collection, however, does not appear to have been updated since June 15, 2007. The Microbial Genome Codon Usage Database  and Prokaryotes Codon Usage Database  have lookup tables of codon counts and frequencies for over 500 and 800 species, respectively (as of June 2011). The most comprehensive database of codon statistics for individual microbial strains can be found in the Codon Usage Bias Database (CUB-DB) ; individual links guide the user to individual strain values for many of the aforementioned metrics as well as two metrics developed by the database’s author .
The codon bias database
We have developed a new web resource called the Codon Bias Database (CBDB) which lists RSCU, normalized RSCU, and frequency biases (FB) values for 300+ and counting bacterial strains. The genera selected for this first release of CBDB all have sequenced, annotated bacteriophages, thus providing a reference for those studying adaptation in phage. Following the metric proposed by Paul Sharp and collaborators [6, 23, 45], analysis is performed for the set of highly expressed genes (HEGs) and these gene sets have been manually curated. CBDB is organized in such a way to easily accommodate codon usage comparisons between strains and species.
Construction and content
In this initial release, over 300 strains belonging to 17 genera were selected. The FASTA format (*.fna) and protein coding gene (*.ptt) files were downloaded from the NCBI FTP site . Code developed in C++ generated *.ptt files containing only the annotations of the HEGs which includes 40 genes previously used for analysis of codon bias ; these genes encode for the translation elongation factor Tu (tufA), Ts (tsf) and G (fusA) as well as 37 ribosomal proteins (rplA rplF rplI rplT, and rpsB rpsT). Each HEG file was then manually inspected to account for variations in naming conventions. The HEG *.ptt files, which include information regarding the location of the gene within the genome, the gene name, and information about the protein product, are available for download for each strain.
Using code developed in C++, the frequency of each codon within each strain’s set of HEG sequences was calculated. These frequencies were then used to calculate three metrics representing the codon usage patterns: the relative synonymous codon usage (RSCU), normalized relative synonymous codon usage (NRSCU) and frequency bias (HEG FB) values. The RSCU metric, introduced by Sharp et al. , is the observed frequency of the codon divided by the frequency expected under the assumption of equal usage of synonymous codons. The commonly used metric of quantifying codon bias, the CAI value, is derived by referencing the set of RSCU values for a species’ HEG sequence set. RSCU values can vary from zero to the number of synonymous codons available for a particular amino acid. The NRSCU value scales the RSCU values such that each codon’s value is between zero and one. Thus, each amino acid is weighted equally while still retaining the variance in usage between synonymous codons. The last metric, the HEG FB metric, is the relative frequency of the codon within the set of all codons present in the HEG sequences. The FB values also reflect the variation in synonymous codon usage, but in contrast to the NRSCU metric it also captures biases in amino acid usage within the HEG set. These metrics are calculated as follows:
where x ij is the number of occurrences of the jth codon for the ith amino acid and n i is the number of synonymous codons which encode for the ith amino acid. The user can download a table containing the results of these calculations for each strain available through CBDB as well as a table containing these values for all strains of the selected genus.
The nucleotide and protein sequences of the HEGs for each strain are also available for download from CBDB in FASTA format.
Utility & discussion
The motivation behind creating the CBDB site was two-fold. Firstly, this site was designed to provide a resource of codon frequency look-up tables for researchers who are investigating translational selection and codon bias in bacteria. Researchers can assess if a nonsynonymous mutation observed is for a codon which is more/less preferred. Secondly, a number of groups, including our own, utilize the phage-host system to explore viral evolution. In particular, CBDB is a resource for recognizing selection within evolving species.
While CUB-DB  includes much of the same data included in CBDB, we wanted to present the data in such a way that one could easily compare the variation in biases between strains and species within a genus. Secondly, CBDB also contains the FASTA format sequences of the highly expressed genes as well as the gene annotation information, which could then be subsequently analyzed by any of the aforementioned codon analysis software tools. Thus, researchers can go to one repository for this data. In contrast with the resource of the Highly Expressed Genes Database or HEG-DB , CBDB includes only the set of HEG identified by Sharp et al.  and the calculation of CAI values. Furthermore, this initial release of CBDB contains more species than is available through HEG-DB .
Figure1 shows a screen shot of the CBDB website, showing the results for the first strain listed for the genus Acinetobacter. Each strain’s name and NCBI RefSeq number, which includes a hyperlink to the NCBI genome record, is listed followed by links to download: the data table (format *.xls), the tab-delimited gene annotation file for the HEG (format *.ptt), the FASTA format amino acid sequences of the HEG set, and the FASTA format nucleotide sequences of the HEG set. Below these links is the data table containing the amino acid abbreviation for each codon and each codon’s RSCU, NRSCU, and HEG FB values. At the top of each genus’ page is a link to download all of the data tables in a single MS Excel document. Also included is a link to a single MS Excel file which contains the RSCU values for each codon for each strain and an interactive graph to compare biases between strains. The list of all of the available genera appears in the navigation on the left pane of each page and all of the species available for a selected genus. Moreover, the navigation pane includes a link describing the codon bias metrics used and associated references.
While in this release of the database only 300+ species are available, additional species will be added on a regular basis; all of the annotated bacterial species available through NCBI will be included over the coming months. Furthermore, we are in the process of developing and making publicly available through the site functionality for conducting inter- and intra-species as well as phage-host comparisons of codon usage profiles.
Exploring codon bias with CBDB
Previous analysis of the variation in codon usage between bacterial species selected just a few representatives of a genus and species . We were interested to see how codon usage varied between species in the same genus as well as between different strains of the same species. For each of the three metrics included in CBDB, we calculated the correlation in codon usage for each pair of strains. A distance matrix was then computed as (1-r)/2 such that the instance of a pair being anticorrelated has a distance of 1 and a perfectly correlated pair has a distance of 0. The FITCH application of the PHYLIP package  was used to derive a tree to visualize similarities/dissimilarities in codon usage. Figure2A shows the tree derived when the NRSCU values were compared for all strains. (Clades were collapsed when the subtree contained only one genus. Shigella was found to group both with Escherichia and Salmonella. Chlamydia and Chlamydophila were also collapsed into single nodes as some strains were closer to strains of the other species than of their own species.) From Figure2A, we also identified two branches including Pseudomonas (maroon circle) and Salmonella & Shigella (orange square) species. The leftmost Salmonella & Shigella branch, separating the pseudomonads is a single species – Salmonella enterica subsp. arizonae serovar 62:z4,z23:-- str. RSK2980 (GenBank: NC_010067). Also from this visualization of codon usage bias we noticed that the 34 Bacillus strains are interspersed amongst other species in the tree. The two B. lincheniformis and single B. clausii strain show a more similar codon usage pattern to the Vibrio fisheri sequences than they do to the other bacilli, as is shown in Figure2B. Looking at the Bacillus genus data file available from CBDB which includes statistics for all of the strains, one can see variation in usage between species.
The genera included in this first release vary in the strength of their bias. For instance, as has been previously documented , the Chlamydiae phylum does not exhibit a significant bias. The eight Chlamydia and Chlamydophila species examined here do not show a strong bias. ANOVAs were performed for all three statistics revealing that the variation between the species is not statistically significant. The Vibrio species, however, are strongly biased  and exhibit differences in codon usage between species. For instance, Figure3 shows the biases for the 11 Vibrio species examined for the six leucine codons. Referring to the Genomic tRNA database , one can find that the non-cholerae species have more TAG-anticodon tRNAs (blue) than any other Leucine tRNAs. This thus can explain the preference within these species for CTA codon. The V. cholerae species also have more TAG-anticodon tRNAs than any other Leucine tRNAs, but in contrast to the non-cholerae species they also have CAG-anticodon tRNAs (green) .
As recent research has found, the use of particular codons can improve translational accuracy and efficiency in addition to serving as a signal for co-translational protein folding [11,19-34]. Comparing usage between species and strains can expose variation in usage between closely related organisms. Furthermore, this site serves as a resource for studying the mechanisms shaping codon usage within bacteriophages. The Codon Bias Database (CBDB) provides a centralized repository of look-up tables and codon usage bias measures for a wide variety of genera, species and strains, facilitating comparisons in codon usage between closely related species such as the one presented here.
Availability and requirements
The Codon Bias Database (CBDB) is freely available without restriction at http://www.cbdb.info. This website has been tested with browsers of Safari, Internet Explorer and Firefox.
AH and JS collected the data as well as participated in the design of the site and writing of the manuscript. CP conceived of the study and site and developed the code used for generating data and analyses in addition to participating in the design of the site and writing of the manuscript. All authors read and approved the final manuscript.
Sueoka N: On the genetic basis of variation and heterogeneity of DNA base composition. Proc Natl Acad Sci USA 1962, 34: 95–114.
Bennetzen JL, Hall BD: Codon selection in yeast. J Biol Chem 1982, 257: 3026–3031.
Gautier C: Compositional bias in DNA. Curr Opin Genet Dev 2000, 10: 656–661. 10.1016/S0959-437X(00)00144-1
dos Reis M, Savva R, Wernisch L: Solving the riddle of codon usage preferences: a test for translational selection. Nucleic Acids Res 2004, 32: 5036–5044. 10.1093/nar/gkh834
Urrutia AO, Hurst LD: Codon usage bias covaries with expression breadth and the rate of synonymous evolution in humans, but this is not evidence for selection. Genetics 2001, 159: 1191–1199.
Sharp PM, Bailes E, Grocock RJ, Peden JF, Sockett RE: Variation in the strength of selected codon usage bias among bacteria. Nucleic Acids Res 2005, 33: 1141–1153. 10.1093/nar/gki242
Bulmer M: The selection-mutation-drift theory of synonymous codon usage. Genetics 1991, 129: 897–907.
Palidwor GA, Perkins TJ, Xia X: A general model of codon bias due to GC mutational bias. PLoS One 2010, 5: e13431. 10.1371/journal.pone.0013431
Supek F, Škunca N, Repar J, Vlahoviček K, Šmuc T: Translational selection is ubiquitous in prokaryotes. PLoS Genet 2010, 6: e1001004. 10.1371/journal.pgen.1001004
Hershberg R, Petrov DA: Selection on codon bias. Annu Rev Genet 2008, 42: 287–299. 10.1146/annurev.genet.42.110807.091442
Shah P, Gilchrist MA: Explaining complex codon usage patterns with selection for translational efficiency, mutation bias, and genetic drift. Proc Natl Acad Sci U S A 2011, 108: 10231–10236. 10.1073/pnas.1016719108
Bulmer M: Coevolution of codon usage and transfer RNA abundance. Nature 1987, 325: 728–730. 10.1038/325728a0
Ikemura T: Correlation between the abundance of transfer RNAs and the occurrence of the respective codons in its protein genes: A proposal for a synonymous codon choice that is optimal for translational system. J Mol Biol 1981, 151: 389–409. 10.1016/0022-2836(81)90003-6
Ikemura T: Codon usage and tRNA content in unicellular and multicellular organisms. Mol Biol Evol 1985, 2: 13–34.
Kanaya S, Yamada Y, Kudo Y, Ikemura T: Studies of codon usage and tRNA genes of 18 unicellular organisms and quantification of Bacillus subtilis tRNAs: gene expression level and species-specific diversity of codon usage based on multivariate analysis. Gene 1999, 238: 143–155. 10.1016/S0378-1119(99)00225-5
Burgess-Brown NA, Sharma S, Sobott F, Loenarz C, Oppermann U, Gileadi O: Codon optimaization can improve expression of human genes in Escherichia coli: A multi-gene study. Protein Expr Purif 2008, 59: 94–102. 10.1016/j.pep.2008.01.008
Maertens B, Spriestersbach A, von Groll U, Roth U, Kubicek J, Gerrits M, Graf M, Liss M, Daubert D, Wagner R, Schäfer F: Gene optimization mechanisms: a multi-gene study reveals a high success rate of full-length human proteins expressed in Escherichia coli. Protein Sci 2010, 19: 1312–1326. 10.1002/pro.408
Welch M, Govindarajan S, Ness JE, Villalobos A, Gurney A, Minshull J, Gustafsson C: Design parameters to control synthetic gene expression in Escherichia coli. PLoS One 2009, 4: e7002. 10.1371/journal.pone.0007002
Akashi H: Synonymous codon usage in Drosophila melanogaster: natural selection and translational accuracy. Genetics 1994, 136: 927–935.
Marais G, Duret L: Synonymous codon usage, accuracy of translation, and gene length in Caenorhabditis elegans. J Mol Evol 2001, 52: 275–280.
Stoletzki N, Eyre-Walker A: Synonymous codon usage in Escherichia coli: selection for translational accuracy. Mol Biol Evol 2007, 24: 374–381.
Zhou T, Weems M, Wilke CO: Translationally optimal codons associate with structurally sensitive sites in proteins. Mol Biol Evol 2009, 26: 1571–1580. 10.1093/molbev/msp070
Sharp PM, Li W-H: The codon adaptation index: a measure of directional synonymous codon usage, and its potential applications. Nucleic Acids Res 1987, 15: 1281–1295. 10.1093/nar/15.3.1281
dos Reis M, Wernisch L, Savva R: Unexpected correlations between gene expression and codon usage bias from microarray data for the whole Escherichia coli K-12 genome. Nucleic Acids Res 2003, 22: 2437–2446.
Brockmann R, Beyer A, Heinisch JJ, Wilhelm T: Posttranscriptional expression regulation: what determines translation rates? PLoS Comput Biol 2007, 3: e57. 10.1371/journal.pcbi.0030057
Tuller T, Carmi A, Vestsigian K, Navon S, Dorfan Y, Zaborske J, Pan T, Dahan O, Furman I, Pilpel Y: An evolutionary conserved mechanism for controlling the efficiency of protein translation. Cell 2010, 141: 344–354. 10.1016/j.cell.2010.03.031
Zhang Z, Zhou L, Hu L, Zhu Y, Xu H, Liu Y, Chen X, Yi X, Kong X, Hurst LD: Nonsense-mediated decay targets have multiple sequence-related features that can inhibit translation. Mol Syst Biol 2010, 6: 442.
Gingold H, Pilpel Y: Determinants of translation efficiency and accuracy. Mol Syst Biol 2011, 7: 481.
Plotkin JB, Kudla G: Synonymous but not the same: the causes and consequences of codon bias. Nat Rev Genet 2011, 12: 32–42. 10.1038/nrg2899
Kudla G, Murray AW, Tollervey D, Plotkin JB: Coding-sequence determinants of gene expression in Escherichia coli. Science 2009, 324: 255–258. 10.1126/science.1170160
Gilchrist MA, Shah P, Zaretzki R: Measuring and detecting molecular adaptation in codon usage against nonsense errors during protein translation. Genetics 2009, 183: 1493–1505. 10.1534/genetics.109.108209
Tsai CJ, Sauna ZE, Kimchi-Sarfaty C, Ambudkar SV, Gottesman MM, Nussinov R: Synonymous mutations and ribosome stalling can lead to altered folding pathways and distinct minima. J Mol Biol 2008, 383: 281–291. 10.1016/j.jmb.2008.08.012
Marin M: Folding at the rhythm of the rare codon beat. Biotechnol J 2008, 3: 1047–1057. 10.1002/biot.200800089
Ivanov IG, Saraffova AA, Abouhaidar MG: Unusual effect of clusters of rare arginine (AGG) codons on the expression of human interferon alpha 1 gene in Escherichia coli. Int J Biochem Cell Biol 1997, 29: 659–666. 10.1016/S1357-2725(96)00161-6
Karlin S, Doerfler W, Cardon LR: Why is CpG suppressed in the genomes of virtually all small eukaryotic viruses but not in those of large eukaryotic viruses? J Virol 1994, 68: 2889–2897.
Shackelton LA, Parrish CR, Holmes EC: Evolutionary basis of codon usage and nucleotide composition bias in vertebrate DNA viruses. J Mol Evol 2006, 62: 551–563. 10.1007/s00239-005-0221-1
Greenbaum BD, Levine AJ, Bhanot G, Rabadan R: Patterns of evolution and host gene mimicry in influenza and other RNA viruses. PLoS Pathog 2008, 4: e1000079. 10.1371/journal.ppat.1000079
Schubert AM, Putonti C: Evolution of the sequence composition of Flaviviruses. Infect Genet Evol 2010, 10: 129–136. 10.1016/j.meegid.2009.11.004
Jenkins GM, Holmes EC: The extent of codon usage bias in human RNA viruses and its evolutionary origin. Vir Res 2003, 92: 1–7. 10.1016/S0168-1702(02)00309-X
Kunisawa T, Kanaya S, Kutter E: Comparison of synonymous codon distribution patterns of bacteriophage and host genomes. DNA Res 1998, 5: 319–326. 10.1093/dnares/5.6.319
Lucks JB, Nelson DR, Kudla GR, Plotkin JB: Genome landscapes and bacteriophage codon usage. PLoS Comput Biol 2008, 4: e1000001. 10.1371/journal.pcbi.1000001
Carbone A: Codon bias is a major factor explaining phage evolution in translationally biased hosts. J Mol Evol 2008, 66: 210–223. 10.1007/s00239-008-9068-6
McLachlan AD, Staden R, Boswell DR: A method for measuring the non-random bias of a codon usage table. Nucleic Acids Res 1984, 12: 9567–9575. 10.1093/nar/12.24.9567
Clarke TF, Clark PL: Rare codons cluster. PLoS One 2008, 3: e3412. 10.1371/journal.pone.0003412
Sharp PM, Tuohy TMF, Mosurski KR: Codon usage in yeast: Cluster analysis clearly differentiates highly and lowly expressed genes. Nucleic Acids Res 1986, 14: 5125–5143. 10.1093/nar/14.13.5125
Carbone A, Zinovyev A, Képès F: Codon Adaptation Index as a measure of dominating codon bias. Bioinformatics 2003, 19: 2005–2015. 10.1093/bioinformatics/btg272
Lee S, Weon S, Lee S, Kang C: Relative codon adaptation index, a sensitive measure of codon usage bias. Evol Bioinform Online 2010, 6: 47–55.
Sharp PM, Emery LR, Zeng K: Forces that influence the evolution of codon bias. Phil Trans R Soc B 2010, 365: 1203–1212. 10.1098/rstb.2009.0305
Puigbò P, Bravo IG, Garcia-Vallvé S: E-CAI: a novel server to estimate an expected value of Codon Adaptation Index (eCAI). BMC Bioinforma 2008, 9: 65. 10.1186/1471-2105-9-65
Harnady M, Wilson SA, Zaneveld J, Sueoka N, Knight R: CodonExplorer: an online tool for analyzing codon usage and sequence composition, scaling from genes to genomes. Bioinformatics 2009, 25: 1331–1332. 10.1093/bioinformatics/btp141
Peden J: CodonW. http://codonw.sourceforge.net/
Puigbò P, Bravo IG, Garcia-Vallvé S: CAIcal: a combined set of tools to access codon usage adaptation. Biol Direct 2008, 3: 38. 10.1186/1745-6150-3-38
Angellotti MC, Bhuiyan SB, Chen G, Wan XF: CodonO: codon usage bias analysis within and across genomes. Nucleic Acids Res 2007, 35: W132-W136. 10.1093/nar/gkm392
Grote A, Hiller K, Scheer M, Münch R, Nörtemann B, Hempel DC, Jahn D: JCat: a novel tool to adapt codon usage of a target gene to its potential expression host. Nucleic Acids Res 2005, 33: W526-W531. 10.1093/nar/gki376
Fuhrmann M, Hausherr A, Ferbitz L, Schödl T, Heitzer M, Hegemann P: Monitoring dynamic expression of nuclear genes in Chlamydomonas reinhardtii by using a synthetic luciferase reporter gene. Plant Mol Biol 2004, 55: 869–881.
Nakamura Y, Gojobori T, Ikemura T: Codon usage tabulated from international DNA sequence databases: status for the year 2000. Nucleic Acids Res 2000, 28: 292. 10.1093/nar/28.1.292
MGCUD:Microbial Genome Codon Usage Database http://bioinformatics.forsyth.org/mgcud/
Zhu W, Lomsadze A, Borodovsky M: Ab initio gene identification in metagenomic sequences. Nucleic Acids Res 2010, 38: e132. 10.1093/nar/gkq275
Raiford D: CUB-DB: Codon Usage Database. http://web-dev.cs.umt.edu/cub_db/
NCBI FTP ftp://ftp.ncbi.nlm.nih.gov/
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.
Felsenstein J: PHYLIP. http://evolution.genetics.washington.edu/phylip.html
Chan PP, Lowe TM: GtRNAdb: A database of transfer RNA genes detected in genomic sequence. Nucleic Acids Res 2009, 37: D93-D97. 10.1093/nar/gkn787
This study was partially funded by a research support grant from Loyola University Chicago (CP).
The authors declare that they have no competing interests.