Table of Contents Author Guidelines Submit a Manuscript
Advances in Bioinformatics
Volume 2009, Article ID 316936, 11 pages
http://dx.doi.org/10.1155/2009/316936
Research Article

Synonymous Codon Usage Analysis of Thirty Two Mycobacteriophage Genomes

Tuberculosis Research Centre Indian Council of Medical Research, Chennai 600 031, India

Received 30 July 2009; Accepted 27 October 2009

Academic Editor: Huixiao Hong

Copyright © 2009 Sameer Hassan et al. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

Abstract

Synonymous codon usage of protein coding genes of thirty two completely sequenced mycobacteriophage genomes was studied using multivariate statistical analysis. One of the major factors influencing codon usage is identified to be compositional bias. Codons ending with either C or G are preferred in highly expressed genes among which C ending codons are highly preferred over G ending codons. A strong negative correlation between effective number of codons (Nc) and GC3s content was also observed, showing that the codon usage was effected by gene nucleotide composition. Translational selection is also identified to play a role in shaping the codon usage operative at the level of translational accuracy. High level of heterogeneity is seen among and between the genomes. Length of genes is also identified to influence the codon usage in 11 out of 32 phage genomes. Mycobacteriophage Cooper is identified to be the highly biased genome with better translation efficiency comparing well with the host specific tRNA genes.

1. Introduction

The genetic code uses 64 codons to represent the 20 standard amino acids and the translation termination signal. Each codon is recognized by a subset of a cell’s transfer ribonucleotide acid molecules (tRNAs), and with the exception of a few codons that have been reassigned in some lineages [1, 2] the genetic code is remarkably conserved, although it is still in a state of evolution [3].

Codons can be grouped into 20 disjoint families and each family in the universal genetic code contains between 1 and 6 codons. The usage of alternate synonymous codons in an organism is understood to be nonrandom. Grantham et al. [4, 5] proposed that each genome has a particular codon usage signature that reflects particular evolutionary forces acting within that genome. The synonymous codons and amino acids are not used at equal frequencies both within and between organisms [57]; the patterns of codon usage vary considerably among organisms, and also among genes from the same genome [8].

Several factors such as directional mutational bias [911], translational selection [1115], secondary structure of proteins [1621], replicational and transcriptional selection [22, 23], and environmental factors [24, 25] have been reported to influence the codon usage in various organisms. In contrast, amino acid usage has been shown to be influenced by factors such as hydrophobicity, aromaticity, cysteine residue (Cys) content, and mean molecular weight (MMW) [24, 2630]. Compositional constraints and translational selection are thought to be the two main factors for the codon usage variation among the genes in and across genomes. Compositional bias shapes the codon usage variation among the genes in the extremely AT or GC rich unicellular organisms [3133]. Analysis of codon usage pattern can provide a basis for understanding the relevant mechanism for biased usage of synonymous codons [34] and has both practical and theoretical importance in understanding the basics of molecular biology [7, 3539].

Bacteriophages generally use the translational machinery of their hosts to synthesize both their structural and regulatory proteins. This indicates that the amount of codon usage in the protein coding genes in the phages and their bacterial hosts should be similar.

Mycobacteriophages have the potential to be used in diagnosis of tuberculosis and as molecular tools to study mycobacteria. Understanding the codon usage pattern of these phages should guide in the selection of appropriate ones for such purposes. The present venture is to study and understand codon usage patterns of all the mycobacteriophages so far sequenced.

Codon usage analysis was previously done for fourteen phages by Sau et al. [40]. Eighteen more mycobacteriophage genomes were subsequently sequenced and became available in Genbank. In the present work we have analysed all these 32 phage genomes to study and compare their codon usage pattern. Other codon usage indices that affect the genomes of these phages are also studied.

2. Materials and Methods

2.1. Sequences

The complete genome sequences of 32 mycobacteriophages were downloaded from GenBank. Genes having more than 100 codons with proper start and stop codons and without any intermediate stop codon were selected for the current study.

2.2. Analysis

Numbers of codons (Ncs), Relative Synonymous Codon Usage (RSCU), and GC composition at every position of codons were calculated for each gene. The analysis was carried out by GCUA [41], CODONW 1.4.2 (http://codonw.sourceforge.net/).

(a)Nc, the “effective number of codons” used in a gene measures the bias away from equal usage of codons within synonymous groups [19]. Nc can take values from 20 to 61, when only one codon or all synonyms in equal frequencies were used per amino acid, respectively. Nc appears to be a good measure of general codon usage bias [19, 42]. The sequences in which Nc values are 30 are highly expressed while those with 55 are poorly expressed genes [12, 43].(b)Relative Synonymous Codon Usage. Relative synonymous codon usage (RSCU) is defined as the ratio of the observed frequency of codons to the expected frequency if all the synonymous codons for those amino acids are used equally [21]. RSCU is used to observe the synonymous codon usage variation among the genes. (c) Base composition. The frequency of A, T, G, C, and GC at first, second, and third positions of synonymously variable sense codons which can potentially vary from 0 to 1 was calculated. The variation of GC3s among genes was characterized by its standard deviation.
2.3. Statistical Methods

Correspondence analysis (CA) is used to study the codon usage variation between genes in different organisms in which the data are plotted in a multidimensional space of 59 axes excluding those of Met, Trp, and stop codons [19]. For understanding the codon usage variation of mycobacteriophages chosen for the current study, RSCU values are used for CA in order to minimize the amino acid composition. To investigate the difference between high and low expressed genes, we have compared the codon usage variation between 10% of the genes located at the extreme right of axis 1 and 10% of the genes located at the extreme left of the axis 1 produced by CA using RSCU. To estimate the codon usage variation between these two sets of genes we have performed Chi square tests taking as significant criterion.

The Pearson correlation coefficient and linear regression were calculated to identify the indices that influence the codon usage variation in mycobacteriophages using SPSS version 10.0. The levels of statistical significance were defined as or .

3. Results and Discussion

3.1. Overall Codon Usage Analysis in Mycobacteriophages

The RSCU values of 32 mycobacteriophage genomes show that G- and/or C-ending codons are predominantly used (Table 1), in which 13 are C-ending and 6 are G-ending codons. This was expected, as these phages have a high genomic content. However, from the overall RSCU values, it can be assumed that compositional constraint is the only factor responsible for shaping the codon usage variation among the genes in these genomes. Although the overall RSCU values could unveil the codon usage pattern for the genomes, it may hide the codon usage variation among different genes in a genome.

tab1
Table 1: Relative synonymous codon usage for all the genes of 32 mycobacteriophages was calculated. AA and N are the amino acid and number of codons, respectively.
3.2. Codon Usage Variation in 32 Mycobacteriophages

The codon usage bias in the coding regions of 32 completely sequenced mycobacteriophages of varying G C content has been investigated. The average values of the effective numbers of codon (Nc) in different mycobacteriophages varied from 31.44 to 47.96 in mycobacteriophage Cooper and mycobacteriophage Barnyard, respectively. Nucleotide usage pattern in third codon position of all the mycobacteriophages showed high codon usage variation (Table 2). The average GC3s values for individual genomes varied from 65.84 to 89.35 in mycobacteriophage Barnyard and mycobacteriophage Cooper, respectively. In addition, there are marked intragenomic variations in Nc and GC3s values with standard deviation of 3.5 in both the indices. There seems to be a considerable heterogeneity in compositional bias and codon usage pattern within and among the genome of these phages. Of the 32 mycobacteriophages, the genome of Cooper is identified to have the lowest Nc and the highest GC3s values while Barnyard has the highest Nc and the lowest GC3s values indicating that highly GC rich genomes are more biased than poor GC rich genomes.

tab2
Table 2: Nc and GC3s values for 32 mycobacteriophages with standard deviation within brackets.

In unicellular organisms, a strong correlation between gene expressivity and the extent of codon usage bias is reported for Escherichia coli and Saccharomyces cerevisiae and phages of Staphycoccus aureus and mycobacteria [13, 40, 4448]. Our analysis reveals that the genome of mycobacteriophage Cooper is highly biased than other 31 mycobacteriophage genomes. Based on the comparison of the highly represented codons of cooper and the copy number of host specific tRNA, the data indicate that the putatively highly expressed genes of this phage have better translational efficiency comparatively and may be expressed rapidly by the host translation machinery.

Table 3 represents the base composition for the 32 completely sequenced mycobacteriophages. GC composition at the third codon position is always higher than the second and first second codon positions observed in other GC-rich genomes [911]. It is also identified that there is major variation in GC3s content among the genomes studied with no major variation in GC1s and GC2s. This suggests that GC3s has a major role to play than GC1s and GC2s and is tightly associated with the codon usage bias of these genomes.

tab3
Table 3: GC % and base composition in the third codon position for 32 mycobacteriophage genomes.
3.3. Synonymous Codon Usage in Different Mycobacteriophages

A plot of Nc versus GC3s (Nc plot) has been widely used to study the codon usage variation among genes in different organisms [19]. It was demonstrated that the comparison of the actual distribution of the genes, with the expected distribution under no selection, indicates that the codon usage bias of the genes has influences other than the compositional bias. In contrast, if GC3s were the only determinants of the codon usage variation among the genes, then the values of Nc would fall on the continuous curve between Nc and GC3s.

It is evident from Nc plot for the mycobacteriophages studied that most of the genes fall within a restricted cloud, at GC3s between 0.65 and 0.93, and Nc values 28 and 47 (Figure 1). Nc values for these genes lie just below the expected curve, indicating that these genes have additional codon usage bias apart from compositional bias. The rest of the genes have higher Nc values and lower GC3s values, mostly lying on and close to the expected curve. Consequently, the Nc values of these genes are substantially higher relative to expected values. However, strong influence of compositional constraints on codon usages bias in all the phages analyzed could be understood from the presence of significant negative correlation between GC3s and Nc (; ).

316936.fig.001
Figure 1: Nc plot of thirty two mycobacteriophages. The genes for individual phages are represented by different colors.
3.4. Differential Base Usage in Third Codon Position

The correlation of the frequencies of four bases in the third position against Nc values of different genes of these 32 mycobacteriophages has been estimated (Table 4). As there is no information of gene expression level of mycobacteriophages so far, we have considered the highly biased genes having low Nc values as those highly expressed and vice versa.

tab4
Table 4: Correlation coefficient of Nc with C3s, T3s, A3s, and G3s base composition.

In all mycobacteriophages analyzed, except Che9d, the frequency of C at the third codon position increases with decreasing Nc values, whereas frequencies of T and A increase with Nc. However, the frequency of G is not influenced in the third codon position excluding for few phages such as Bxb1, Barnyard, Corndog, Plot, PBI1, and Bethelhem. Thus the influence of mutational bias of these phages is reflected in the choice of bases in the third position. However, this is expected since the optimal codons are, in general, chosen in accordance with the mutational bias of these phages. In other words, it is due to the translational selection that the mutational bias appears to be more prominent in the third codon position of highly expressed genes [18].

3.5. Effect of Translational Selection on the Synonymous Codon Usage in Mycobacteriophage Genomes

Some of the earlier reports showed that synonymous codon usage in the highly expressed genes of diverse array of organisms is influenced by cellular tRNA abundance [44, 4952]. Kanaya et al. [52, 53] have reported that the cellular tRNA abundance in several organisms are directly proportional to their tRNA copy number.

Of the 32 mycobacteriophages analyzed, 10 phages 244 (2 tRNA genes), Bxz1 (26 tRNA genes), Omega (2 tRNA genes), Cjw1 (1 tRNA gene), Wildcat (21 tRNA genes), D29 (5 tRNA genes), L5 (3 tRNA genes), Che12 (3 tRNA genes), Catera (26 tRNA genes) encode tRNA genes for few of the amino acids. Both Bxz1 and Catera are identified to encode large number of tRNA genes (26 tRNA genes) and Wildcat encoding for 21 tRNA genes. Of the 10 phage genomes encoding tRNA genes, excluding 244 and Omega, eight carry tRNA genes for the overrepresented codons in highly and lowly expressed genes.

To see whether the synonymous codon usage of putatively highly expressed genes of these mycobacteriophage genomes is positively correlated with the host tRNA abundance, the number of over represented synonymous codons in such genes was determined by comparing with that of the putatively lowly expressed genes. It was found that among the 22 overrepresented synonymous codons in highly expressed genes, 21 codons are recognized by the M.tuberculosis specific tRNAs (data not shown). Based on the above analysis, the data indicate that the putatively highly expressed genes of these phages have translational efficiency.

3.6. Relationship between Codon Bias and Gene Length

Selection for translational accuracy is predicted to have a positive correlation between codon bias and gene length [20]. Previously, the relationship between gene length and synonymous codon usage bias has been reported for Drosophila melanogaster, Escherichia coli, Saccharomyces cerevisiae, Pseudomonas aeruginosa, and Yersinia pestis [11, 15]. From the plot drawn with gene length against Nc (Figure 2), it is understood that shorter genes have a much wider variance in Nc values, and vice versa for longer genes. Lower Nc values in longer genes may be due to the direct effect of translation time on fitness or to the extra energy cost of proofreading associated with longer translating time. A significant correlation was identified in 11 phages revealing that gene length influences codon usage of these genomes (Table 5). Similar results were also reported for S. pneumoniae, P. aeruginosa and SARS coronavirus [11, 16]. Eyre-Walker [20] has reported that the selection for fidelity in protein translation is likely to be greater in longer genes because the cost of producing a protein is proportional to its length. Therefore selection of translational accuracy predicts a positive correlation between codon usage bias and gene length. And this selection may be stronger at constrained codons coding for evolutionarily conserved amino acids than the nonconserved amino acid. As the codon bias is lower in longer genes than shorter ones, further analysis in finding these constrained codons will help us in understanding whether it is the same in all genes irrespective of their length.

tab5
Table 5: Correlation coefficient of gene length with Nc and GC3s values.
316936.fig.002
Figure 2: Plot of Nc versus Gene length for all mycobacteriophage genomes.
3.7. Correspondence Analysis Using RSCU Values

In order to determine the factors that influence variations in codon usage among the genes of mycobacteriophage genomes, correspondence analysis was conducted on the RSCU values of its genes. Only the distributions of the genes along the first two major axes were shown, as these accounted for 13.63% and 6.89% of the total variation (Figure 3).

316936.fig.003
Figure 3: Correspondence analysis of Relative Synonymous Codon Usage values of mycobacteriophages (32 genomes).

The first major axis is negatively correlated with G3s () and C3s () but correlated positively with A3s () and T3s (). Interestingly, high degree of positive correlation exists between position of genes along the first axis with Nc () (Figure 4) and high degree of negative correlation with GC3s () (Figure 5). These findings suggest that highly biased genes, those with G- and C-ending codons, are clustered on the negative side, whereas the codons ending in A and T predominate on the positive side of the first major axis.

316936.fig.004
Figure 4: Scatter plot of mycobacteriophages and Nc values.
316936.fig.005
Figure 5: Scatter plot of mycobacteriophages and GC3s values.

Additionally, significant negative correlation is observed with Nc against GC3s and GC. Highly expressed genes tend to use “C” or “G” at the synonymous positions compared with lowly expressed genes. It is also studied that C-ending codons are preferred over G-ending codons in highly expressed genes. Preference of C-ending codons in the highly expressed genes might be related to the translational efficiency of the genes as it has been reported that RNY (R-purine, N- any nucleotide base, and Y-pyrimidine) codons are more advantageous for translation [54]. Thus, compositional mutation bias possibly plays an important role in shaping the genome of these phages.

The genomes of Llij, PMC, Wildcat, TM4, Che8, Tweety, and Che9d phages showed no significant correlation between first major axis and GC3s. Whereas, phages such as 244, Bxb1, Bxz2, Che9c, Rosebush, Omega, Halo, Barnyard, Bxz1, Cjw1, Corndog, Orion, Plot, Qyrzula, and Giles show strong negative correlation with GC3s. The primary trend in codon usage variation in these phages can be attributed to the presence of putatively foreign genes acquired through horizontal gene transfer with unusually A T rich codon usage. However, in phages such as D29, L5, PBI1, PG1, Cooper, Che12, Catera, Bethelhem, and U2, axis 1 coordinates are significantly positively correlated with GC3s values (Table 6). Moreover, when G3s and C3s are considered separately, the correlation coefficient exhibited by the positions of genes along the first axis with C3s is significantly larger than that with G3s (Table 6), indicating that the contribution of C3s to the interspecies variation in overall GC3s content is greater than that of G3s.

tab6
Table 6: Correlation of Axis1 with other codon usage indices.

Table 7 shows RSCU values for each codon for the two groups of genes. The asterisk represents the codons whose occurrences are significantly higher in the genes situated on the extreme left side of axis 1, compared to the genes present on the extreme right of the first major axis. It is important to note that out of 22 codons that are statistically overrepresented in genes located on the extreme left side of axis 1 there is 16 C-ending codons and 5-G ending codons. UGA is the most frequent stop codon among highly and lowly expressed genes. Similar pattern is also seen in Mycobacterium tuberculosis genome, where the highly expressed genes prefer codons ending with “C” and “G” [18].

tab7
Table 7: Relative Synonymous Codon Usage for the highly and lowly expressed genes.

4. Conclusion

Compositional bias and translational forces had been reported to play a major role in shaping the codon usage of 14 mycobacteriophages. Our observations corroborate with the earlier report with respect to all the 32 mycobacteriophages. Gene length has a minor role in the selection of codon usage of 11 out of 32 mycobacteriophage genes analyzed. High level of heterogeneity is seen within and among the mycobacteriophage genomes. Cooper is identified to be the highly biased genome with better translation efficiency comparing well with the host specific tRNA genes.

Acknowledgments

S. Hassan acknowledges ICMR—Biomedical Informatics Project and Tuberculosis Research Centre for the funding provided. The authors acknowledge the Department of Biotechnology, New Delhi, for their financial support of this initiative under the collaborative study titled “Establishment of National Database on Tuberculosis.” They thank Dr. P. Venkatesan, Gomathi Sivaramakrishnan, Azger Dusthackeer, and Balaji Subramanyam of Tuberculosis Research Centre for their help in preparation of the manuscript.

References

  1. S. Osawa and T. H. Jukes, “Codon reassignment (codon capture) in evolution,” Journal of Molecular Evolution, vol. 28, no. 4, pp. 271–278, 1989. View at Publisher · View at Google Scholar
  2. S. Osawa, A. Muto, T. H. Jukes, and T. Ohama, “Evolutionary changes in the genetic code,” Proceedings of the Royal Society B, vol. 241, no. 1300, pp. 19–28, 1990. View at Google Scholar
  3. S. Osawa, T. H. Jukes, K. Watanabe, and A. Muto, “Recent evidence for evolution of the genetic code,” Microbiological Reviews, vol. 56, no. 1, pp. 229–264, 1992. View at Google Scholar
  4. R. Grantham, C. Gautier, M. Gouy, R. Mercier, and A. Pavé, “Codon catalog usage and the genome hypothesis,” Nucleic Acids Research, vol. 8, no. 1, pp. r49–r62, 1980. View at Google Scholar
  5. R. Grantham, C. Gautier, and M. Gouy, “Codon frequencies in 119 individual genes confirm consistent choices of degenerate bases according to genome type,” Nucleic Acids Research, vol. 8, no. 9, pp. 1893–1912, 1980. View at Google Scholar
  6. C. E. Martin and S. Scheinbach, “Expression of proteins encoded by foreign genes in Saccharomyces cerevisiae,” Biotechnology Advances, vol. 7, no. 2, pp. 155–185, 1989. View at Google Scholar
  7. A. T. Lloyd and P. M. Sharp, “Evolution of codon usage patterns: the extent and nature of divergence between Candida albicans and Saccharomyces cerevisiae,” Nucleic Acids Research, vol. 20, no. 20, pp. 5289–5295, 1992. View at Google Scholar
  8. R. J. Grocock and P. M. Sharp, “Synonymous codon usage in Cryptosporidium parvum: identification of two distinct trends among genes,” International Journal for Parasitology, vol. 31, no. 4, pp. 402–412, 2001. View at Publisher · View at Google Scholar
  9. G. D'Onofrio and G. Bernardi, “A universal compositional correlation among codon positions,” Gene, vol. 110, no. 1, pp. 81–88, 1992. View at Publisher · View at Google Scholar
  10. S. Majumdar, S. K. Gupta, V. S. Sundararajan, and T. C. Ghosh, “Compositional correlation studies among the three different codon positions in 12 bacterial genomes,” Biochemical and Biophysical Research Communications, vol. 266, no. 1, pp. 66–71, 1999. View at Publisher · View at Google Scholar · View at PubMed
  11. S. K. Gupta and T. C. Ghosh, “Gene expressivity is the main factor in dictating the codon usage variation among the genes in Pseudomonas aeruginosa,” Gene, vol. 273, no. 1, pp. 63–70, 2001. View at Publisher · View at Google Scholar
  12. P. M. Sharp and E. Cowe, “Synonymous codon usage in Saccharomyces cerevisiae,” Yeast, vol. 7, no. 7, pp. 657–678, 1991. View at Publisher · View at Google Scholar · View at PubMed
  13. F. Alvarez, C. Robello, and M. Vignali, “Evolution of codon usage and base contents in kinetoplastid protozoans,” Molecular Biology and Evolution, vol. 11, no. 5, pp. 790–802, 1994. View at Google Scholar
  14. E. N. Moriyama and J. R. Powell, “Gene length and codon usage bias in Drosophila melanogaster, Saccharomyces cerevisiae and Escherichia coli,” Nucleic Acids Research, vol. 26, no. 13, pp. 3188–3193, 1998. View at Publisher · View at Google Scholar
  15. Z.-C. Hou and N. Yang, “Factors affecting codon usage in Yersinia pestis,” Sheng Wu Hua Xue Yu Sheng Wu Wu Li Xue Bao Shanghai, vol. 35, no. 6, pp. 580–586, 2003. View at Google Scholar
  16. W. Gu, T. Zhou, J. Ma, X. Sun, and Z. Lu, “Analysis of synonymous codon usage in SARS Coronavirus and other viruses in the Nidovirales,” Virus Research, vol. 101, no. 2, pp. 155–161, 2004. View at Publisher · View at Google Scholar · View at PubMed
  17. Z.-C. Hou and N. Yang, “Analysis of factors shaping S. pneumoniae codon usage,” Yi Chuan Xue Bao, vol. 29, no. 8, pp. 747–752, 2002. View at Google Scholar
  18. A. Pan, C. Dutta, and J. Das, “Codon usage in highly expressed genes of Haemophillus influenzae and Mycobacterium tuberculosis: translational selection versus mutational bias,” Gene, vol. 215, pp. 405–413, 1998. View at Google Scholar
  19. F. Wright, “The ‘effective number of codons’ used in a gene,” Gene, vol. 87, no. 1, pp. 23–29, 1990. View at Publisher · View at Google Scholar
  20. A. Eyre-Walker, “Synonymous codon bias is related to gene length in Escherichia coli: selection for translational accuracy?” Molecular Biology and Evolution, vol. 13, no. 6, pp. 864–872, 1996. View at Google Scholar
  21. P. M. Sharp and W. H. Li, “The codon Adaptation Index—a measure of directional synonymous codon usage bias, and its potential applications,” Nucleic Acids Research, vol. 15, no. 3, pp. 1281–1295, 1987. View at Google Scholar
  22. A. Ranjan, A. S. Vidyarthi, and R. Poddar, “Evaluation of codon bias perspectives in phage therapy of Mycobacterium tuberculosis by multivariate analysis,” In Silico Biology, vol. 7, no. 4-5, article 0030, pp. 423–431, 2007. View at Google Scholar
  23. S. K. Gupta, T. K. Bhattacharyya, and T. C. Ghosh, “Synonymous codon usage in Lactococcus lactis: mutational bias versus translational selection,” Journal of Biomolecular Structure & Dynamics, vol. 21, no. 4, pp. 527–536, 2004. View at Google Scholar
  24. S. Basak, T. Banerjee, S. K. Gupta, and T. C. Ghosh, “Investigation on the causes of codon and amino acid usages variation between thermophilic Aquifex aeolicus and mesophilic Bacillus subtilis,” Journal of Biomolecular Structure & Dynamics, vol. 22, no. 2, pp. 205–214, 2004. View at Google Scholar
  25. D. J. Lynn, G. A. C. Singer, and D. A. Hickey, “Synonymous codon usage is subject to selection in thermophilic bacteria,” Nucleic Acids Research, vol. 30, no. 19, pp. 4272–4277, 2002. View at Google Scholar
  26. J. R. Lobry and C. Gautier, “Hydrophobicity, expressivity and aromaticity are the major trends of amino-acid usage in 999 Escherichia coli chromosome-encoded genes,” Nucleic Acids Research, vol. 22, no. 15, pp. 3174–3180, 1994. View at Google Scholar
  27. B. Garat and H. Musto, “Trends of amino acid usage in the proteins from the unicellular parasite Giardia lamblia,” Biochemical and Biophysical Research Communications, vol. 279, no. 3, pp. 996–1000, 2000. View at Publisher · View at Google Scholar · View at PubMed
  28. A. Zavala, H. Naya, H. Romero, and H. Musto, “Trends in codon and amino acid usage in Thermotoga maritima,” Journal of Molecular Evolution, vol. 54, no. 5, pp. 563–568, 2002. View at Publisher · View at Google Scholar · View at PubMed
  29. T. Banerjee, S. Basak, S. K. Gupta, and T. C. Ghosh, “Evolutionary forces in shaping the codon and amino acid usages in Bochmannia floridanus,” Journal of Biomolecular Structure & Dynamics, vol. 22, no. 1, pp. 13–23, 2004. View at Google Scholar
  30. H. Naya, A. Zavala, H. Romero, H. Rodríguez-Maseda, and H. Musto, “Correspondence analysis of amino acid usage within the family Bacillaceae,” Biochemical and Biophysical Research Communications, vol. 325, no. 4, pp. 1252–1257, 2004. View at Publisher · View at Google Scholar · View at PubMed
  31. S. Ohkubo, A. Muto, Y. Kawauchi, F. Yamao, and S. Osawa, “The ribosomal protein gene cluster of Mycoplasma capricolum,” Molecular & General Genetics, vol. 210, no. 2, pp. 314–322, 1987. View at Publisher · View at Google Scholar
  32. F. Wright and M. J. Bibb, “Codon usage in the G+C-rich Streptomyces genome,” Gene, vol. 113, no. 1, pp. 55–65, 1992. View at Publisher · View at Google Scholar
  33. S. K. Gupta, T. K. Bhattacharyya, and T. C. Ghosh, “Compositional correlation and codon usage studies in Buchnera aphidicola,” Indian Journal of Biochemistry & Biophysics, vol. 39, no. 1, pp. 35–48, 2002. View at Google Scholar
  34. H. Lü, W.-M. Zhao, Y. Zheng, H. Wang, M. Qi, and X.-P. Yu, “Analysis of synonymous codon usage bias in Chlamydia,” Acta Biochimica et Biophysica Sinica, vol. 37, no. 1, pp. 1–10, 2005. View at Google Scholar
  35. J. W. Fickett, “Recognition of protein coding regions in DNA sequences,” Nucleic Acids Research, vol. 10, no. 17, pp. 5303–5318, 1982. View at Google Scholar
  36. H. Chiapello, E. Ollivier, C. Landes-Devauchelle, P. Nitschke, and J.-L. Risler, “Codon usage as a tool to predict the cellular location of eukaryotic ribosomal proteins and aminoacyl-tRNA synthetases,” Nucleic Acids Research, vol. 27, no. 14, pp. 2848–2851, 1999. View at Publisher · View at Google Scholar
  37. A. Martín, J. Bertranpetit, J. L. Oliver, and J. R. Medina, “Variation in G+C-content and codon choice: differences among synonymous codon groups in vertebrate genes,” Nucleic Acids Research, vol. 17, no. 15, pp. 6181–6189, 1989. View at Google Scholar
  38. E. N. Moriyama and D. L. Hartl, “Codon usage bias and base composition of nuclear genes in Drosophila,” Genetics, vol. 134, no. 3, pp. 847–858, 1993. View at Google Scholar
  39. J. O. McInerney, “Replicational and transcriptional selection on codon usage in Borrelia burgdorferi,” Proceedings of the National Academy of Sciences of the United States of America, vol. 95, no. 18, pp. 10698–10703, 1998. View at Google Scholar
  40. K. Sahu, S. K. Gupta, S. Sau, and T. C. Ghosh, “Comparative analysis of the base composition and codon usages in fourteen mycobacteriophage genomes,” Journal of Biomolecular Structure & Dynamics, vol. 23, no. 1, pp. 63–71, 2005. View at Google Scholar
  41. J. O. McInerney, “GCUA: general codon usage analysis,” Bioinformatics, vol. 14, no. 4, pp. 372–373, 1998. View at Google Scholar
  42. J. M. Comeron and M. Aguadé, “An evaluation of measures of synonymous codon usage bias,” Journal of Molecular Evolution, vol. 47, no. 3, pp. 268–274, 1998. View at Publisher · View at Google Scholar
  43. K. Sau, S. K. Gupta, S. Sau, and T. C. Ghosh, “Synonymous codon usage bias in 16 Staphylococcus aureus phages: implication in phage therapy,” Virus Research, vol. 113, no. 2, pp. 123–131, 2005. View at Publisher · View at Google Scholar · View at PubMed
  44. R. Grantham, C. Gautier, M. Gouy, M. Jacobzone, and R. Mercier, “Codon catalog usage is a genome strategy modulated for gene expressivity,” Nucleic Acids Research, vol. 9, no. 1, pp. r43–r74, 1981. View at Google Scholar
  45. J. L. Bennetzen and B. D. Hall, “Codon selection in yeast,” Journal of Biological Chemistry, vol. 257, no. 6, pp. 3026–3031, 1982. View at Google Scholar
  46. M. Gouy and C. Gautier, “Codon usage in bacteria: correlation with gene expressivity,” Nucleic Acids Research, vol. 10, no. 22, pp. 7055–7074, 1982. View at Google Scholar
  47. P. M. Sharp and W.-H. Li, “An evolutionary perspective on synonymous codon usage in unicellular organisms,” Journal of Molecular Evolution, vol. 24, no. 1-2, pp. 28–38, 1986. View at Publisher · View at Google Scholar
  48. P. M. Sharp and K. M. Devine, “Codon usage and gene expression level in Dictyostelium discoideum: highly expressed genes do ‘prefer’ optimal codons,” Nucleic Acids Research, vol. 17, no. 13, pp. 5029–5039, 1989. View at Google Scholar
  49. P. M. Sharp, M. S. Rogers, and D. J. McConnell, “Selection pressures on codon usage in the complete genome of bacteriophage T7,” Journal of Molecular Evolution, vol. 21, no. 2, pp. 150–160, 1985. View at Publisher · View at Google Scholar
  50. T. Ikemura, Transfer RNA in Protein Synthesis, CRC Press, London, UK, 1992.
  51. J. Zhou, W. J. Liu, S. W. Peng, X. Y. Sun, and I. Frazer, “Papillomavirus capsid protein expression level depends on the match between codon usage and tRNA availability,” Journal of Virology, vol. 73, no. 6, pp. 4972–4982, 1999. View at Google Scholar
  52. S. Kanaya, Y. Yamada, Y. Kudo, and T. Ikemura, “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, vol. 238, no. 1, pp. 143–155, 1999. View at Publisher · View at Google Scholar
  53. S. Kanaya, Y. Yamada, M. Kinouchi, Y. Kudo, and T. Ikemura, “Codon usage and tRNA genes in eukaryotes: correlation of codon usage diversity with translation efficiency and with CG-dinucleotide usage as assessed by multivariate analysis,” Journal of Molecular Evolution, vol. 53, no. 4-5, pp. 290–298, 2001. View at Publisher · View at Google Scholar · View at PubMed
  54. J. C. Shepherd, “Method to determine the reading frame of a protein from the purine/pyrimidine genome sequence and its possible evolutionary justification,” Proceedings of the National Academy of Sciences of the United States of America, vol. 78, no. 3, pp. 1596–1600, 1981. View at Google Scholar