Research Article | Open Access
Codon Usage Patterns in Corynebacterium glutamicum: Mutational Bias, Natural Selection and Amino Acid Conservation
The alternative synonymous codons in Corynebacterium glutamicum, a well-known bacterium used in industry for the production of amino acid, have been investigated by multivariate analysis. As C. glutamicum is a GC-rich organism, G and C are expected to predominate at the third position of codons. Indeed, overall codon usage analyses have indicated that C and/or G ending codons are predominant in this organism. Through multivariate statistical analysis, apart from mutational selection, we identified three other trends of codon usage variation among the genes. Firstly, the majority of highly expressed genes are scattered towards the positive end of the first axis, whereas the majority of lowly expressed genes are clustered towards the other end of the first axis. Furthermore, the distinct difference in the two sets of genes was that the C ending codons are predominate in putatively highly expressed genes, suggesting that the C ending codons are translationally optimal in this organism. Secondly, the majority of the putatively highly expressed genes have a tendency to locate on the leading strand, which indicates that replicational and transciptional selection might be invoked. Thirdly, highly expressed genes are more conserved than lowly expressed genes by synonymous and nonsynonymous substitutions among orthologous genes fromthe genomes of C. glutamicum and C. diphtheriae. We also analyzed other factors such as the length of genes and hydrophobicity that might influence codon usage and found their contributions to be weak.
It is well established that the codon usage patterns are generally not used with equal frequency. Grantham et al. firstly explained the phenomena of unequal usage and proposed the “genome hypothesis", stating that the biases are species specific , and multivariate analysis methods were used to analyze codon usage and amino acid composition [2–4]. As more and more complete genome sequences of diverse species are investigated, researchers have found that biased usage of synonymous codons may result from various factors. Some unicellular species have extremely biased compositions, where compositional constraints are the main factors in determining the codon usage variation among genes [5–7]. In contrast, both translational selection and compositional constraint operate on the codon usage variation in other organisms [8–14]. Moreover, in several bacteria, the replication and translational selection is responsible for the codon usage variation among genes [15–18]. In organisms, such as Escherichia coli , Drosophila melanogaster , and Caenorhabditis elegans , the frequency of codon usage is directly proportional to the corresponding tRNA population and the preferred codons in highly expressed genes are recognized by the most abundant tRNAs. Meanwhile, it has been reported that amino acid conservation and hydrophobicity are the main factors shaping codon usage among the genes in Mycobacteria [21, 22]. Other factors may also influence the synonymous codon usage, such as protein secondary structure [23–26], mRNA folding stability [27, 28], gene function [29, 30], and gene length [31–33].
Corynebacterium glutamicum ATCC 13032, used industrially for the production of amino acids, is an aerobic, gram-positive rod-shaped bacteria capable of growing on a variety of sugars or organic acids . In this study, we used the available complete genome sequence of this organism and analyzed its codon usage, aiming to understand the genetic organization of the C. glutamicum genome. Our results show that mutational bias, natural selection, and amino acid conservation are the main factors driving codon usage patterns in C. glutamicum genes.
2. Materials and Methods
2.1. Genome Sequence Data
The complete genome sequences and coding sequences of C. glutamicum and C. diphtheriae were obtained from the NCBI ftp site (ftp://ftp.ncbi.nih.gov/genomes/Bacteria/). To minimize sampling errors, only genes of at least 100 codons in length with correct initiation and termination codons were used in further analysis.
2.2. Multivariate Analysis of Codon Usage
The COA (codon usage correspondence analysis, plots the codon usage data in a multidimensional space of 59 axes, excluding Met, Trp, and termination codons, identifies the axes which represent the most prominent factors contributing to the variation among genes), GC3s (the frequency of at the third synonymously variable coding position, excluding Met, Trp, and termination codons), ENC (the “effective number of codons"; a measure of the bias in codon usage of genes, usually highly expressed genes display lower values compared with lowly expressed ones), RSCU (the “relative synonymous codon usage"; a value greater than 1.0 indicates that the corresponding codon is more frequently used than expected, whereas the reverse is true for RSCU values less than 1.0), CAI (the “codon adaptation index"; high values mean higher codon usage bias and higher expressed level), Fop (the “frequency of optional codons"), GRAVY index of hydrophobicity, and A3s, G3s, C3s, and T3s (the composition of each individual base A, G, C, and T at the third synonymous codon positions) were performed using the program CodonW1.42 (http://codonw.sourceforge.net/).The CAI was calculated taking the codon usage of the ribosomal proteins as a reference. Other statistical analyses were performed with the SPSS statistical software version 11.0.
2.3. Locating Genes Situated on the Leading and Lagging Strands of Replication
Asymmetrical mutational bias between the two complementary strands may contribute to variations in codon usage. To locate the genes on the leading or lagging strand of replication, the sites of origin and termination were determined by using the oriloc program (ftp://pbil.univlyon1.fr/pub/logiciel/oriloc/oriloc.c) and GC skew (/) was determined using the GC Skewing program (http://www.genomicsplace.com/gc_skew/gc_skew.html) by taking a 24 kb window size and a step size of 3 kb to locate the leading and lagging strands.
2.4. Orthologous Gene Pairs and Analysis
Orthologous genes were identified by the reciprocal best blast hit approach as those pairs displaying value of 60% identity, an E-value of and overlapped by at least 60% of the length of the longest protein, with at least 100 amino acids in length using the local BLASTP program (http://www.ncbi.nlm.nih.gov/BLAST/download.shtml). The protein sequences of 1525 orthologous gene pairs were aligned using the MUSCLE program (http://www.drive5.com/muscle); then the aligned protein sequences were used to generate the corresponding codon alignment. The Ka (the number of synonymous substitutions per site) and Ks (the number of nonsynonymous substitutions per site) for each pair of aligned sequences were estimated using the PAML version 4.3 package (http://abacus.gene.ucl.ac.uk/software/paml.html) with runmode and CodonFreq . Only those pairs of sequences having Ks values below 1.0 were considered in further analysis and the final dataset was comprised of 437 gene pairs.
3.1. Overall Codon Usage
As shown in Figure 1(a), the genome of C. glutamicum is biased towards high contents ranging from 40% to 68% with an average of 54.7% and a standard deviation of 3.7%. With the exception of small regions, its genome shows little variation around the mean value. Due to composition constraints, G and C are expected to predominate at the third position of codons. Indeed, the codon usage indicated that C ending codons are predominant overall (data not shown). In order to understand the codon usage variation among different genes, ENC and GC3s values were calculated (Figure 1(b)). ENC values vary from 24.46 to 61.00 with a mean of 46.9 and standard deviation of 7.55%. The heterogeneity of codon usage was further confirmed from the GC3s values ranging from 28% to 87% with a mean of 57.18% and standard deviation of 8.3%. Wright suggested that plotting ENC against GC3s values could be used to effectively explore codon usage variation among genes . If GC3s are the only determination of the codon usage variation among genes, then the values of ENC would fall on the continuous curve. The GC3s versus ENC plot reveals that only a small proportion of points lie on the expected curve (Figure 1(b)), which indicates that apart from the effect of compositional constraints, there might be some additional factors driving codon usage variation among the genes.
3.2. Gene Expression and Codon Usage Bias
In order to investigate the other possible trends in shaping codon usage variation among the genes in C. glutamicum, we subjected the data to multivariate statistical analysis. Figure 1(c) shows the position of genes along the first two axes. At the positive end of the first axis, it comprises of putatively highly expressed genes, such as ribosomal proteins, translation elongation factors, while the majority of putatively lowly expressed genes are scattered towards the other extreme. A more important result emerged when the genes were sorted according to their respective CAI values, and the highest positions were displayed not only by the genes encoding ribosomal protein but also by almost the same genes along the extreme of the first axis. Table 1 shows the first axis accounts for 20.33%, compared with 10.5% of the second axis and this value of the first axis is high and much larger than that of the second axis, indicating a primary trend in codon usage across genes. Furthermore, there are positive correlations between the first axis and CAI (, ), with Fop (, ), with GC3s (, ), and especially with C3s (, ). Those results suggest that gene expression may be the main factor shaping the codon usage in this organism, the first axis is associated with expression levels, and highly expressed genes have higher () content, especially C content at their synonymous third codon position than lowly expressed genes.
**Correlation is significant at the 0.01 level.|
To investigate the differences between highly and lowly expressed genes, we compared the codon usage of genes that locate the two extremes of the first axis (Table 2). Chi square tests were performed taking as the significant criterion. We found that there were 22 coding codons (corresponding to 18 amino acids) that are more highly used in putatively highly expressed genes than putatively lowly expressed genes. Among the 20 codons, there are 14 C ending codons and 3 G ending codons, which demonstrate that the presumed highly expressed genes tend to be C3-rich.
|: the number of codons; AA: amino acid.|
*Codon with significantly () higher frequencies in highly expressed genes.
High: codons in highly expressed genes; Low: codons in lowly expressed genes.
3.3. Replicational and Transcriptional Selection and Codon Usage
Recent reports of several bacterial strains show that codon usage bias is mainly governed by transcriptional and translational selection [19, 20, 36, 37]. After the origin versus termination and leading versus lagging strands were determined, we located the genes on the leading or lagging strands of replication and found that the proportion of genes located on the leading strands increases with CAI, from about 55% for low CAI genes (0.35) to 67% for high CAI genes (0.65) in the organism. For the putatively highly expressed ribosomal proteins, the proportion of genes on the leading strands reaches 84% (44/52) (Table 3). This observation is consistent with previous research results that essential genes are enriched to a greater extent than nonessential genes in the leading strand .
3.4. Gene Conservation and Codon Usage
The rate of synonymous substitutions has been reported to be nonuniform among different genes in the same species . When we calculated the Ka and Ks between the orthologous genes from C. glutamicum and C. diphtheriae, several results were determined. Firstly, there is a negative correlation between the Ka and CAI value (, ), comparative with Ks and CAI values with and (Figure 2). When the genes are sorted according to the respective Ks, the genes displaying the lowest values are those presumed highly expressed genes, such as ribosomal protein and translation elongation factors., Taken together, this indicates that highly expressed genes have diverged less at the synonymous position than lowly expressed genes. Secondly, the Ka and Ks are correlated with , . Thirdly, there is significant correlation between Ks and Fop (, ), which indicated that the genes diverging less are the ones displaying highest frequencies of optional codon usage.
Finally, we also investigated the relationship between codon usage and gene length (, ), codon usage, and hydrophobicity (, ), suggesting that their contributions to the codon usage variation are weak.
Among prokaryotes, it is generally accepted that the preferences of synonymous codons can be explained as the result of mutational bias and natural selection acting at the level of translation. In C. glutamicum, the composition bias towards GC constraint indicates that these bases are predominant at the third codon positions across all genes. Indeed, the putatively highly expressed genes show an increment of several codons, most of which are C-ending triplets. Ikemura showed that there is a match between these codons and the most abundant tRNAs . In Escherichia coli , Drosophila melanogaster , and Caenorhabditis elegans , highly expressed genes have a strong selective preference for codons with a high concentration for the corresponding acceptor tRNA molecule; the preferred codons are those best recognized by the most abundant tRNAs. This trend has been interpreted as the coadaptation between amino acid composition of protein and tRNA-pools to enhance the translational efficiency. Remarkably, in this study, there is a strong positive correlation (, ) between the Fop in each gene and respective CAI value. This strongly suggests that translational selection influenced the codon usage of C. glutamicum and the “optional codons" were more frequent in highly expressed genes.
As more prokaryotic genomes are analyzed, it becomes evident that codon usage is rather dependent on mutational bias and natural selection. For example, the complex pattern of codon usage in Chlamydia trachomatis is inferred to be the result of strand-specific mutation, natural selection, the hydropathy level of each protein and amino acid conservation . In this study, we present evidence suggesting that, apart from mutational bias and natural selection, strand-specific and amino acid conservation also contribute to the codon usage of C. glutamicum. Strand bias also dominates codon usage in other symbiotic or parasitic bacteria, such as Rickettsia prowazekii, Borrelia burgdorferi, and Lawsonia intracellularis [15, 40, 41]. We found a distribution bias of genes (particularly for those with a high CAI) on the leading strands in C. glutamicum. This is usually interpreted as the result of “replicational selection”, by which presence on the leading strand would permit the avoidance of collision between polymerases when replication and transcription occur at the same time .
It was reported that the codon usage is more biased for amino acid that are more conserved between species [42, 43]; natural selection has a larger contribution than mutation to the observed correlation between evolutionary rates and gene expression level in Chlamydomonas . A correlation between Ks and Fop was also identified. This correlation with Ka might be explained in many ways. Akashi argued that the selection for translation accuracy maintains a high frequency of preferred codons for highly conserved amino acids . Two additional hypotheses for this pattern are a possible mechanistic bias in mutation and the fact that synonymous sites are also subject to some degree of selection . The latter scenario could mean either selection on codon usage, or that synonymous substitutions might not always be silent or evolutionary responses to adaptations . A similar interaction between the level of expression, the level of codon bias, and gene conservation was demonstrated in Mycobacterium .
In summary, this study has shown that the codon usage variation among the genes of C. glutamicum is influenced by mutational bias, translational selection, and amino acid conservation. As more complete prokaryotic genomes are being studied, different factors shaping the pattern of codon usage might be found.
This work is supported by the National Natural Science Foundation of China (30571009). G. Liu and J. Wu contributed equally to this work.
- 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.
- C. Medigue, T. Rouxel, P. Vigier, A. Henaut, and A. Danchin, “Evidence for horizontal gene transfer in Escherichia coli speciation,” Journal of Molecular Biology, vol. 222, no. 4, pp. 851–856, 1991.
- G. Pascal, C. Médigue, and A. Danchin, “Universal biases in protein composition of model prokaryotes,” Proteins, vol. 60, no. 1, pp. 27–35, 2005.
- G. Pascal, C. Médigue, and A. Danchin, “Persistent biases in the amino acid composition of prokaryotic proteins,” BioEssays, vol. 28, no. 7, pp. 726–738, 2006.
- T. Ohama, A. Muto, and S. Osawa, “Role of GC-biased mutation pressure on synonymous codon choice in Micrococcus luteus, a bacterium with a high genomic GC-content,” Nucleic Acids Research, vol. 18, no. 6, pp. 1565–1569, 1990.
- S. G. E. Andersson and P. M. Sharp, “Codon usage in the Mycobacterium tuberculosis complex,” Microbiology, vol. 142, no. 4, pp. 915–925, 1996.
- S. G. E. Andersson and P. M. Sharp, “Codon usage and base composition in Rickettsia prowazekii,” Journal of Molecular Evolution, vol. 42, no. 5, pp. 525–536, 1996.
- M. Malumbres, J. A. Gil, and J. F. Martin, “Codon preference in corynebacteria,” Gene, vol. 134, no. 1, pp. 15–24, 1993.
- T. C. Ghosh, S. K. Gupta, and S. Majumdar, “Studies on codon usage in Entamoeba histolytica,” International Journal for Parasitology, vol. 30, no. 6, pp. 715–722, 2000.
- H. Romero, A. Zavala, and H. Musto, “Compositional pressure and translational selection determine codon usage in the extremely GC-poor unicellular eukaryote Entamoeba histolytica,” Gene, vol. 242, no. 1-2, pp. 307–311, 2000.
- H. Musto, H. Romero, and A. Zavala, “Translational selection is operative for synonymous codon usage in Clostridium perfringens and Clostridium acetobutylicum,” Microbiology, vol. 149, no. 4, pp. 855–863, 2003.
- S. Das, S. Paul, S. Chatterjee, and C. Dutta, “Codon and amino acid usage in two major human pathogens of genus Bartonella—optimization between replicational-transcriptional selection, translational control and cost minimization,” DNA Research, vol. 12, no. 2, pp. 91–102, 2005.
- K. Sau, S. Sau, S. C. Mandal, and T. C. Ghosh, “Factors influencing the synonymous codon and amino acid usage bias in AT-rich Pseudomonas aeruginosa phage PhiKZ,” Acta Biochimica et Biophysica Sinica, vol. 37, no. 9, pp. 625–633, 2005.
- K. Sau, S. K. Gupta, S. Sau, S. C. Mandal, and T. C. Ghosh, “Factors influencing synonymous codon and amino acid usage biases in Mimivirus,” BioSystems, vol. 85, no. 2, pp. 107–113, 2006.
- 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.
- B. Lafay, A. T. Lloyd, M. J. McLean, K. M. Devine, P. M. Sharp, and K. H. Wolfe, “Proteome composition and codon usage in spirochaetes: species-specific and DNA strand-specific mutational biases,” Nucleic Acids Research, vol. 27, no. 7, pp. 1642–1649, 1999.
- H. Romero, A. Zavala, and H. Musto, “Codon usage in Chlamydia trachomatis is the result of strand-specific mutational biases and a complex pattern of selective forces,” Nucleic Acids Research, vol. 28, no. 10, pp. 2084–2090, 2000.
- N. Stoletzki and A. Eyre-Walker, “Synonymous codon usage in Escherichia coli: selection for translational accuracy,” Molecular Biology and Evolution, vol. 24, no. 2, pp. 374–381, 2007.
- E. N. Moriyama and J. R. Powell, “Codon usage bias and tRNA abundance in Drosophila,” Journal of Molecular Evolution, vol. 45, no. 5, pp. 514–523, 1997.
- L. Duret, “tRNA gene number and codon usage in the C. elegans genome are co-adapted for optimal translation of highly expressed genes,” Trends in Genetics, vol. 16, no. 7, pp. 287–289, 2000.
- A. B. de Miranda, F. Alvarez-Valin, K. Jabbari, W. M. Degrave, and G. Bernardi, “Gene expression, amino acid conservation, and hydrophobicity are the main factors shaping codon preferences in Mycobacterium tuberculosis and Mycobacterium leprae,” Journal of Molecular Evolution, vol. 50, no. 1, pp. 45–55, 2000.
- T. Zhou, X. Sun, and Z. Lu, “Synonymous codon usage in environmental chlamydia UWE25 reflects an evolutional divergence from pathogenic chlamydiae,” Gene, vol. 368, no. 1-2, pp. 117–125, 2006.
- M. Oresic and D. Shalloway, “Specific correlations between relative synonymous codon usage and protein secondary structure,” Journal of Molecular Biology, vol. 281, no. 1, pp. 31–48, 1998.
- T. Xie and D. Ding, “The relationship between synonymous codon usage and protein structure,” FEBS Letters, vol. 434, no. 1-2, pp. 93–96, 1998.
- S. K. Gupta, S. Majumdar, T. K. Bhattacharya, and T. C. Ghosh, “Studies on the relationships between the synonymous codon usage and protein secondary structural units,” Biochemical and Biophysical Research Communications, vol. 269, no. 3, pp. 692–696, 2000.
- W. Gu, T. Zhou, J. Ma, X. Sun, and Z. Lu, “The relationship between synonymous codon usage and protein structure in Escherichia coli and Homo sapiens,” BioSystems, vol. 73, no. 2, pp. 89–97, 2004.
- J. V. Chamary, J. L. Parmley, and L. D. Hurst, “Hearing silence: non-neutral evolution at synonymous sites in mammals,” Nature Reviews Genetics, vol. 7, no. 2, pp. 98–108, 2006.
- B. Kahali, S. Basak, and T. C. Ghosh, “Reinvestigating the codon and amino acid usage of S. cerevisiae genome: a new insight from protein secondary structure analysis,” Biochemical and Biophysical Research Communications, vol. 354, no. 3, pp. 693–699, 2007.
- R. J. Epstein, K. Lin, and T. W. Tan, “A functional significance for codon third bases,” Gene, vol. 245, no. 2, pp. 291–298, 2000.
- A. Fuglsang, “Strong associations between gene function and codon usage,” APMIS, vol. 111, no. 9, pp. 843–847, 2003.
- 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.
- L. Duret and D. Mouchiroud, “Expression pattern and, surprisingly, gene length shape codon usage in Caenorhabditis, Drosophila, and Arabidopsis,” Proceedings of the National Academy of Sciences of the United States of America, vol. 96, no. 8, pp. 4482–4487, 1999.
- G. Marais and L. Duret, “Synonymous codon usage, accuracy of translation, and gene length in Caenorhabditis elegans,” Journal of Molecular Evolution, vol. 52, no. 3, pp. 275–280, 2001.
- J. Kalinowski, B. Bathe, D. Bartels et al., “The complete Corynebacterium glutamicum ATCC 13032 genome sequence and its impact on the production of L-aspartate-derived amino acids and vitamins,” Journal of Biotechnology, vol. 104, no. 1–3, pp. 5–25, 2003.
- A. Fuglsang, “The ‘effective number of codons’ revisited,” Biochemical and Biophysical Research Communications, vol. 317, no. 3, pp. 957–964, 2004.
- T. Ikemura, “Correlation between the abundance of Escherichia coli transfer RNAs and the occurrence of the respective codons in its protein genes: a proposal for a synonymous codon choice that is optimal for the E. coli translational system,” Journal of Molecular Biology, vol. 151, no. 3, pp. 389–409, 1981.
- R. Banerjee and D. Roy, “Codon usage and gene expression pattern of Stenotrophomonas maltophilia R551-3 for pathogenic mode of living,” Biochemical and Biophysical Research Communications, vol. 390, no. 2, pp. 177–181, 2009.
- E. P. C. Rocha and A. Danchin, “Essentiality, not expressiveness, drives gene-strand bias in bacteria,” Nature Genetics, vol. 34, no. 4, pp. 377–378, 2003.
- Y. Kawahara and T. Imanishi, “A genome-wide survey of changes in protein evolutionary rates across four closely related species of Saccharomyces sensu stricto group,” BMC Evolutionary Biology, vol. 7, article 9, 2007.
- J. J. Davis and G. J. Olsen, “Modal Codon Usage: assessing the typical codon usage of a genome,” Molecular Biology and Evolution, vol. 27, no. 4, pp. 800–810, 2010.
- F.-B. Guo and J.-B. Yuan, “Codon usages of genes on chromosome, and surprisingly, genes in plasmid are primarily affected by strand-specific mutational biases in lawsonia intracellularis,” DNA Research, vol. 16, no. 2, pp. 91–104, 2009.
- A. Ticher and D. Graur, “Nucleic acid composition, codon usage, and the rate of synonymous substitution in protein-coding genes,” Journal of Molecular Evolution, vol. 28, no. 4, pp. 286–298, 1989.
- H. Akashi, “Synonymous codon usage in Drosophila melanogaster: natural selection and translational accuracy,” Genetics, vol. 136, no. 3, pp. 927–935, 1994.
- C. E. Popescu, T. Borza, J. P. Bielawski, and R. W. Lee, “Evolutionary rates and expression level in chlamydomonas,” Genetics, vol. 172, no. 3, pp. 1567–1576, 2006.
- I. K. Jordan, I. B. Rogozin, Y. I. Wolf, and E. V. Koonin, “Essential genes are more evolutionarily conserved than are nonessential genes in bacteria,” Genome Research, vol. 12, no. 6, pp. 962–968, 2002.
- D. A. Drummond and C. O. Wilke, “The evolutionary consequences of erroneous protein synthesis,” Nature Reviews Genetics, vol. 10, no. 10, pp. 715–724, 2009.
Copyright © 2010 Guiming Liu 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.