Quantitative relationship between synonymous codon usage bias and GC composition across unicellular genomes
© Wan et al; licensee BioMed Central Ltd. 2004
Received: 28 January 2004
Accepted: 28 June 2004
Published: 28 June 2004
Codon usage bias has been widely reported to correlate with GC composition. However, the quantitative relationship between codon usage bias and GC composition across species has not been reported.
Based on an informatics method (SCUO) we developed previously using Shannon informational theory and maximum entropy theory, we investigated the quantitative relationship between codon usage bias and GC composition. The regression based on 70 bacterial and 16 archaeal genomes showed that in bacteria, SCUO = -2.06 * GC3 + 2.05*(GC3)2 + 0.65, r = 0.91, and that in archaea, SCUO = -1.79 * GC3 + 1.85*(GC3)2 + 0.56, r = 0.89. We developed an analytical model to quantify synonymous codon usage bias by GC compositions based on SCUO. The parameters within this model were inferred by inspecting the relationship between codon usage bias and GC composition across 70 bacterial and 16 archaeal genomes. We further simplified this relationship using only GC3. This simple model was supported by computational simulation.
The synonymous codon usage bias could be simply expressed as 1+ (p/2)log2(p/2) + ((1-p)/2)log2((l-p)/2), where p = GC3. The software we developed for measuring SCUO (codonO) is available at http://digbio.missouri.edu/~wanx/cu/codonO.
All amino acids except Met and Trp are encoded by more than one codon. DNA sequence data from diverse organisms have shown that synonymous codons for any amino acid are not used with equal frequency, even though choices among codons should be equivalent in terms of protein sequences [1–6]. Previous codon usage analyses showed that codon usage bias is very complicated and is associated with various biological factors, such as gene expression level [7–10], gene length [11–13], gene translation initiation signal , protein amino acid composition [6, 15], protein structure [16, 17], tRNA abundance [18–21], mutation frequency and patterns [22, 23], and GC composition [24–27]. In this paper, we further explore the relationship between codon usage and GC composition.
GC composition may be described at three levels: 1) Overall GC content. The overall genome GC content in living organisms varies from 25–75% . However, within a single gene, the overall GC content is 7–95%. 2) Local GC composition. Local GC composition is defined based on the positions on the genetic codons. GC1 is the GC composition at the first site of codons, GC2 the GC composition at the second site of codons, and GC3 the GC composition at the third site of codons. 3) The ratio of G/C or A/T within a single strand of DNA. Based on Watson-Crick base pairing rules, the overall GC content is the same between the plus strand and minus strand of the DNA sequence . However, within a single strand, the ratio of G to C and A to T may not be 1. Bacterial genomes were found to be relatively enriched in G over C and T over A, and slightly depleted in G+C, in their weakly evolution-selected positions (intergenic regions and third codon positions) in the leading strand compared with the lagging strand . Although the overall GC content for a genome is reduced by AT-rich intergenic regions of the genome, the gene GC composition is tightly correlated with the genome's overall GC composition. The higher the overall GC content bias, the higher the local GC composition (GC1, GC2, and GC3).
It has been generally accepted that genome GC content is correlated with amino acid usage and codon usage . A very low or very high GC composition is associated with a large codon usage bias. Recently, through a simple mutational model, Knight et al.  showed that it was the GC composition that drives codon and amino-acid usage although both mutation and selection play important roles. By using the corresponding analysis of codon usage over 32 bacterial and 8 archaeal genomes, Lynn et al.  further showed that codon usage bias was affected by GC composition and environment (e.g., temperature).
However, there has been no systematic, quantitative evaluation of the relationship between codon usage bias and GC composition. In addition, the theoretical basis underlying the relationship between codon usage bias and GC composition has not been illustrated. In this paper, we applied an informatics method , which is based on Shannon informatics theory and the entropy theory, to explore the relationship between GC composition and synonymous codon usage bias among 70 bacterial and 16 archaeal genomes. We presented an analytical model to quantify the non-linear relationship between GC3 and a measurement of codon usage bias (synonymous codon usage order, or SCUO), which reveals that GC3 is the key factor driving synonymous codon usage and that this mechanism is independent of species. These results were supported by our simulations. Our results also showed that the asymmetric distribution of G over C and A over T at the third codon position may increase codon usage bias. The underlying mechanisms behind GC composition and codon usage bias are discussed.
We recently developed an informatics method  to provide an estimate for the orderliness of synonymous codon usage (SCUO) and the amount of synonymous codon usage bias. This method was based on the Shannon informatics theory and the entropy theory and allows the comparison of codon usage bias within and across genomes.
To calculate SCUO, we created a codon table for the amino acids that have more than one codon, indexed in an arbitrary way, so that we could unambiguously refer to the j-th (degenerate) codon of amino acid i, 1 ≤ i ≤ 18. In mycoplasmas, Trp was also included into the codon table since a standard stop codon TGA encodes Trp in this specific species so that 1 ≤ i ≤ 19. For simplicity, the following description of the method is only based on the standard genetic codon table although the actual SCUO computation considered special cases for different organisms.
Let n i represent the number of degenerate codons for amino acid i, so 1 ≤ j ≤ n i : for example, 1 ≤ j ≤ 6 for leucine, 1 ≤ j ≤ 2 for tyrosine, etc. For each sequence, let x ij represent the number of times that synonymous codon j of amino acid i is present, 1 ≤ i ≤ 18, 1 ≤ j ≤ n i . Normalizing the x ij by their sum over j gives the frequency of the j-th degenerate codon for amino acid i in each sequence.
According to the information theory, we define the entropy H ij of the i-th amino acid of the j-th codon in each sequence by
H ij = -p ij log p ij 2
Summing over the codons representing amino acid i gives the entropy of the i-th amino acid in the each sequence
If the synonymous codons for the i-th amino acid were used at random, one would expect a uniform distribution of them as representatives for the i-th amino acid. Thus, the maximum entropy for the i-th amino acid in each sequence is
If only one of the synonymous codons is used for the i-th amino acid, i.e., the usage of the synonymous codons is biased to the extreme, then the i-th amino acid in each sequence has the minimum entropy:
Unlike Shannon's definition of information, Gatlin  and Layzer  define the information as the difference between the maximum entropy and the actual entropy as an index of orderliness. The greater the information, the more ordered the sequence will be . In our case, this information measures the nonrandomness in synonymous codon usage and therefore describes the degree of orderliness for synonymous codon usage for the i-th amino acid in each sequence.
Let O i be the normalized difference between the maximum entropy and the observed entropy for the i-th amino acid in each sequence, i.e.,
Obviously, 0 ≤ O i ≤ 1. When the synonymous codon usage for the i-th amino acid is random, O i = 0. When this usage is biased to the extreme, 0 i = 1. Thus, O i can be thought of as a measure of the bias in synonymous codon usage for the i-th amino acid in each sequence. We designate the statistics O i as the synonymous codon usage order (SCUO) for the i-th amino acid in each sequence.
Let F i be the composition ratio of the i-th amino acid in each sequence:
Then the average SCUO for each sequence can be represented as
A software package called codonO was written using the C programming language to calculate SCUO for each open reading frame (ORF). This program is available at http://digbio.missouri.edu/~wanx/cu/codonO/. The 86 unicellular genomes explored in this paper can be found at our Web site http://digbio.missouri.edu/~wanx/cu/genomelist.htm.
Theoretical model between GC3 and codon usage
From the standard genetic code for synonymous codons, each amino acid was encoded by xy U/C and/or xy A/G. The third site of each codon should be purine (R, either U or C) or pyrimidine (Y, either A or G). For those amino acids with double genetic codons, they can only be coded by xy U/C or xy A/G. We define = A, = U, and vice versa. We assume the probability of G+C at the GC3 is p. Then the probability of A+T at the GC3 is 1-p. Let α be P(G), i.e., the probability of G, and β be P(A), i.e., the probability of A. Hence, 0 < α ≤ p, 0 < β ≤ 1 - p. If we apply Equation 2, we will obtain entropy for each amino acid - α log2α if the third site is G or -(p - α)log2(p - α) if the third position is C (we use base 2 through our theoretical model). The entropy for each amino acid will be - β log2β if the third site is A or -(1-p-β)log2(1-p-β) if the third position is T. Let us use Gln as an example, which is coded by CGA/G. If the sequence was coded only by this single amino acid, α = p and β = 1-p. Therefore, by applying Equation 3, we can calculate the overall entropy for synonymous codons, -p log2p - (1-p)log2(1-p). In this case, the maximum entropy will be log22 = 1. Thus, the overall codon bias for Gln will be 1 + p log2p + (1-p)log2(1-p). For all of the amino acids with 2 genetic codons, we can use the same equation:
SCUO n2 = 1 + p log2p + (1-p)log2(1-p) 10
We can apply the similar deducing process for amino acids with 3, 4 or 6 genetic codes. Based on our definition of SCUO (Equation 9), the overall codon usage bias will be the sum of each amino acid. Thus, we obtain this first-order approximation for the overall codon bias:
SCUO = 1 + μ*SCUOn2 + ν*SCUOn3 + ω*SCUOn4 + γ*SCUOn6 11
Adjusted factors in archaea and bacteria
To simplify the analysis process, we simplify the analytical model as a binary selection model. Among the 18 synonymous codons, only Ile has three genetic codons. All other synonymous codons are encoded in pairs with xyZ and xy . Assuming no additional codon bias results from a biased distribution between A and T or between C and G, i.e., P(A) = P(T), P(G) = P(C), the synonymous codons bias will be simplified to be a binary selection model:
SCUO = 1+ (p/2)log2(p/2) + ((l-p)/2)log2((l-p)/2) 12
In this case, the codon bias is solely determined by GC3, p. This simplified estimation between codon usage bias and GC3, which will reflect the lowest boundary of the codon usage bias for a specific GC3.
GC composition affects SCUO within a single genome
Relationship between GC composition and SCUO across archaeal and bacterial genomes
We plotted the simplified binary selection model in the GC3 vs. SCUO of both archaea and bacteria. The binary selection model closely fit the lowest boundary of the correlation between GC3 and SCUO. The codon usage bias was generally above the plotted curve because in most cases, P(A) is not equal P(T) and P(G) is not equal to P(C) in the natural world. The distribution of Ile may also increase the codon usage bias.
Simulation of codon usage bias and GC composition
Impact of asymmetric distribution of G/C and A/T on codon usage bias
In this paper, we applied the informatics method  to measure SCUO and GC composition within and across unicellular genomes. We first explored the relationship between SCUO and GC composition in E. coli. Our results demonstrated that the relationship between codon usage bias and GC3 formed a "U" shape, which may reflect a direct mathematic consequence of GC3 over codon usage bias. It is interesting that there was a weak horn in the center of "U" shape (Fig. 1), which is similar to previous reports . The genes in the central horns were discussed in Karlin et al. . We can detect some weak horns in some other genomes, such as Lactococcus lactis. However, the central horns were not shown within some other genomes we explored (data not shown). The comparison between GC composition and SCUO confirmed previous reports that GC3 was the most important factor in codon bias among GC, GC1, GC2, and GC3.
The comparison of codon usage and GC composition across 16 archaeal and 70 bacterial genomes showed that a similar, non-linear relationship existed between GC3 and SCUO across the unicellular cells. This also confirmed that GC3 is the main driving force for the codon usage bias. These results were simulated using pseudo ORFs.
To understand the theoretical basis for the relationship between GC composition and codon usage bias, we derived a binary selection model of GC3 based on the Shannon informatics theory and the entropy concept. The non-linear relationship agreed very well with the theoretical equations (Figs. 3 and 4).
Similar to our results with unicellular organisms, the nonlinear relationship between GC3 and codon usage bias has been reported in human and mouse . This indicates that the relationship is independent of species. However, compared to genes in the same species, it is not known whether the features associated with codon usage bias, such as gene expression level, protein structure, etc., are associated with codon usage bias across species. It may be interesting to investigate the relationships between those reported gene features [7–27] and codon usage bias of those horizontally transferred genes in their new adapted hosts. The methodology described here provides a simple way (SCUO) to compare the codon usage bias across species.
GC3 were demonstrated to be tightly associated with cell functional significance [38, 39]. Sequence analysis of human receptor tyrosine kinase genes demonstrated that functionally important transmembrane hydrophobic amino acids are specified by codons containing a higher GC frequency at the third bases than are transmembrane neutral amino acids.
We have demonstrated that the asymmetric distribution of G over C and T over A at the third codon positions increased both the codon usage bias values and ranges. The asymmetric distribution of mutational pressure at the third codon positions provides more flexible selection ability during the environmental adaptation process. The asymmetric distribution has been observed in many bacterial genomes, which were found to be relatively enriched in G over C and T over A at GC3 in the leading strand compared with the lagging strand . The effect of asymmetric mutation pressure on the amino-acid composition of proteins has been reported elsewhere [40–42]. Thus, the asymmetric distribution of mutational pressure contributes to the codon usage bias besides the GC3, although we did not include such an effect in our simplified binary selection model.
It is interesting that the GC3 is not asymmetrically distributed (Fig. 2,3). One possible reason might be that the GC content of protein coding region is overall higher than the non-coding region. The other possible reason might be the length of the distribution. Within a single genome, we found that the shorter sequences have a wider distribution of GC3s (data not shown). Within our simulation, we assumed the even distribution of the gene length between 200 and 1000. This might result in the symmetric distribution GC3 in the simulation.
In summary, we developed a simple binary selection model that mimicked the quantitative relationship between the codon usage bias and GC composition in the unicellular organisms, which was supported by systematically characterizing the relationship between codon usage bias and GC composition among 86 unicellular organisms. Our simulation results support this finding by demonstrating that the asymmetric distribution of mutational pressure at the third codon positions has an impact on codon usage bias.
The bacteria and archaea genomic sequences and annotations were downloaded from ftp://ftp.ncbi.nlm.nih.gov/genomes/Bacteria/ in August, 2002.
To simulate the relationship between codon usage bias and GC composition, we generated 15,514 pseudo ORFs with a random size of 200–1000 codons. For each ORF, the nucleotides (A,T,G, or C) for three positions (1s, 2s, and 3s) were generated separately using the following procedure: 1) generate nucleotides for the third codon position randomly with a GC content between 0 and 1; 2) generate nucleotides for the first codon position with a GC content between 0.42*GC3 + 0.18 and 0.42*GC3 + 0.46, which represents a similar linear relationship between GC3 and GC1 in actual genomes, as shown in the Results section; 3) generate nucleotides for the second codon position with a GC content between 0.27*GC3 + 0.11 and 0.27*GC3 + 0.39, which represents a similar linear relationship between GC3 and GC2 in actual genomes; 4) randomly generate A/T and G/C with equal or random probability between A and T, as well as between G and C, for each position based on the associated GC composition; and 5) combine the generated nucleotides as an ORF, and discard any pseudo ORF containing a stop codon. We calculated SCUO for each pseudo ORF using codonO.
Lists of abbreviations used
Synonymous Codon Usage Order
The GC composition for the first site of each codon
The GC composition for the second site of each codon
The GC composition for the third site of each codon
The authors wish to acknowledge Dorothea Thompson for her critical review and Lynn Kszos for editorial assistance. We also appreciate two anonymous reviewers for their invaluable comments and suggestions. The research of XFW and JZ was supported by the US DOE Office of Science as part of its Biological and Environmental Research Programs in Genome To Life and Microbial Genome Programs. DX's work was funded by the US Department of Energy's Genomes to Life program http://www.doegenomestolife.org under project, "Carbon Sequestration in Synechococcus Sp.: From Molecular Machines to Hierarchical Modeling" http://www.genomes-to-life.org. Oak Ridge National Laboratory is managed by the University of Tennessee-Battelle LLC for the Department of Energy under contract DE-AC05-00OR22725.
- Grantham R, Gautier C, Gouy M: Codon frequencies in 119 individual genes confirm consistent choices of degenerate bases according to genome type. Nucleic Acids Res. 1980, 8: 1893-1912.PubMed CentralView ArticlePubMed
- Aota S, Ikemura T: Diversity in G + C content at the third position of codons in vertebrate genes and its cause. Nucleic Acids Res. 1986, 14: 6345-6355.PubMed CentralView ArticlePubMed
- Murray EE, Lotzer J, Eberle M: Codon usage in plant genes. Nucleic Acids Res. 1989, 17: 477-198.PubMed CentralView ArticlePubMed
- Sharp PM, Cowe E, Higgins DG, Shields DC, Wolfe KH, Wright F: Codon usage patterns in Escherichia coli, Bacillus subtilis, Saccharomyces cerevisiae, Schizosaccharomyces pombe, Drosophila melanogaster and Homo sapiens; a review of the considerable within-species diversity. Nucleic Acids Res. 1988, 16: 8207-8211.PubMed CentralView ArticlePubMed
- Shields DC, Sharp PM, Higgins DG, Wright F: " Silent" sites in Drosophila genes are not neutral: evidence of selection among synonymous codons. Mol Biol Evol. 1988, 5: 704-716.PubMed
- D'Onofrio G, Mouchiroud D, Aissani B, Gautier C, Bernardi G: Correlations between the compositional properties of human genes, codon usage, and amino acid composition of proteins. J Mol Evol. 1991, 32: 504-510.View ArticlePubMed
- Gouy M, Gautier C: Codon usage in bacteria: correlation with gene expressivity. Nucleic Acids Res. 1982, 10: 7055-7074.PubMed CentralView ArticlePubMed
- Sharp PM, Li WH: An evolutionary perspective on synonymous codon usage in unicellular organisms. J Mol Evol. 1986, 24: 28-38.View ArticlePubMed
- Sharp PM, Tuohy TM, Mosurski KR: Codon usage in yeast: cluster analysis clearly differentiates highly and lowly expressed genes. Nucleic Acids Res. 1986, 14: 5125-5143.PubMed CentralView ArticlePubMed
- Sharp PM, Li WH: The codon Adaptation Index – a measure of directional synonymous codon usage bias, and its potential applications. Nucleic Acids Res. 1987, 15: 1281-1295.PubMed CentralView ArticlePubMed
- Bains W: Codon distribution in vertebrate genes may be used to predict gene length. J Mol Biol. 1987, 197: 379-388.View ArticlePubMed
- Eyre-Walker A: Synonymous codon bias is related to gene length in Escherichia coli: selection for translational accuracy?. Mol Biol Evol. 1996, 13: 864-872.View ArticlePubMed
- Wan X, Xu D, Zhou J: A new informatics method for measuring synonymous codon usage bias. In Intelligent engineering systems through artificial neural networks. Edited by: Dagli et al. 2003, ASME Press, New York, NY, 13: 1101-1018.
- Ma J, Campbell A, Karlin S: Correlations between Shine-Dalgarno sequences and gene features such as predicted expression levels and operon structures. J Bacteriol. 2002, 184: 5733-5745. 10.1128/JB.184.20.5733-5745.2002.PubMed CentralView ArticlePubMed
- Lobry JR, Gautier C: Hydrophobicity, expressivity and aromaticity are the major trends of amino-acid usage in 999 Escherichia coli chromosome-encoded genes. Nucleic Acids Res. 1994, 22: 3174-3180.PubMed CentralView ArticlePubMed
- Xie T, Ding D, Tao X, Dafu D: The relationship between synonymous codon usage and protein structure. FEBS Lett. 1998, 434: 93-96. 10.1016/S0014-5793(98)00955-7.View ArticlePubMed
- D'Onofrio G, Ghosh TC, Bernardi G: The base composition of the genes is correlated with the secondary structures of the encoded proteins. Gene. 2002, 300: 179-187. 10.1016/S0378-1119(02)01045-4.View ArticlePubMed
- Ikemura T: 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. J Mol Biol. 1981, 151: 389-409.View ArticlePubMed
- Ikemura T: Correlation between the abundance of yeast transfer RNAs and the occurrence of the respective codons in protein genes. Differences in synonymous codon choice patterns of yeast and Escherichia coli with reference to the abundance of isoaccepting transfer RNAs. J MolBiol. 1982, 158: 573-597.View Article
- Bulmer M: Coevolution of codon usage and transfer RNA abundance. Nature. 1987, 325: 728-730. 10.1038/325728a0.View ArticlePubMed
- 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.View ArticlePubMed
- Sueoka N: Directional mutation pressure, selective constraints, and genetic equilibria. J Mol Evol. 1999, 34: 95-114.
- Sueoka N: Two aspects of DNA base composition: G+C content and translation-coupled deviation from intra-strand rule of A = T and G = C. J Mol Evol. 1999, 49: 49-62.View ArticlePubMed
- Bernardi G: Compositional constraints and genome evolution. J Mol Evol. 1986, 24: 1-11.View ArticlePubMed
- Karlin S, Mrazek J: What drives codon choices in human genes?. J Mol Biol. 1996, 262: 459-72. 10.1006/jmbi.1996.0528.View ArticlePubMed
- Antezana MA, Kreitman M: The nonrandom location of synonymous codons suggests that reading frame-independent forces have patterned codon preferences. J Mol Evol. 1999, 49: 36-43.View ArticlePubMed
- Sueoka N, Kawanishi Y: DNA G+C content of the third codon position and codon usage biases of human genes. Gene. 2000, 261: 53-62. 10.1016/S0378-1119(00)00480-7.View ArticlePubMed
- Lobry JR, Sueoka N: Asymmetric directional mutation pressures in bacteria. Genome Biol. 2002, 3: RESEARCH0058-10.1186/gb-2002-3-10-research0058.PubMed CentralView ArticlePubMed
- Watson JD, Crick FHC: A structure for deoxyribose nucleic acid. Nature. 1953, 171: 737-738.View ArticlePubMed
- Ermolaeva MD: Synonymous codon usage in bacteria. Curr Issues Mol Biol. 2001, 3: 91-97.PubMed
- Knight RD, Freeland SJ, Landweber LF: A simple model based on mutation and selection explains trends in codon and amino-acid usage and GC composition within and across genomes. Genome Biol. 2001, 2: RESEARCH0010-PubMed CentralPubMed
- Lynn DJ, Singer GA, Hickey DA: Synonymous codon usage is subject to selection in thermophilic bacteria. Nucleic Acids Res. 2002, 30: 4272-277. 10.1093/nar/gkf546.PubMed CentralView ArticlePubMed
- Gatlin LL:Information Theory and the Living System. 1972, Columbia University Press,
- Layzer D: Information in cosmology, physics and biology. Int J Quantum Chem. 1977, 12 (suppl 1): 185-195.
- Brooks DR, Wiley EO:Evolution as entropy: toward a unified theory of biology. 1988, The University of Chicago Press, Chicago, 2,
- Karlin S, Mrazek J, Campbell AM: Codon usages in different gene classes of the Escherichia coli genome. Mol Microbiol. 1998, 29: 1341-1355. 10.1046/j.1365-2958.1998.01008.x.View ArticlePubMed
- Zeeberg B: Shannon information theoretic computation of synonymous codon usage biases in coding regions of human and mouse genomes. Genome Res. 2002, 12: 944-955. 10.1101/gr.213402.PubMed CentralView ArticlePubMed
- Smith NGC, Eyre-Walker A: Why are translationally sub-optimal synonymous codons used in Esherichia coli?. J Mol Evol. 2001, 53: 225-236. 10.1007/s002390010212.View ArticlePubMed
- Epstein RJ, Lin K, Tan TW: A functional significance for codon third bases. Gene. 2000, 245: 291-298. 10.1016/S0378-1119(00)00042-1.View ArticlePubMed
- Lafay B, Lloyd AT, McLean MJ, Devine KM, Sharp PM, Wolfe KH: Proteome composition and codon usage in spirochaete: species-specific and DNA strand- specific mutational biases. Nucleic Acids Res. 1999, 27: 1642-1649. 10.1093/nar/27.7.1642.PubMed CentralView ArticlePubMed
- Rocha EPC, Danchin A, Viari A: Universal replication biases in bacteria. Mol Microbiol. 1999, 32: 11-16. 10.1046/j.1365-2958.1999.01334.x.View ArticlePubMed
- Mackiewicz P, Gierlik A, Kowalczuk M, Dudek MR, Cebrat S: How does replication-associated mutational pressure influence amino acid composition of proteins?. Genome Res. 1999, 9: 409-416.PubMed Central
This article is published under license to BioMed Central Ltd. This is an Open Access article: verbatim copying and redistribution of this article are permitted in all media for any purpose, provided this notice is preserved along with the article's original URL.