![](corehtml/pmc/pmcents/thinsp.gif)
![](corehtml/pmc/pmcents/thinsp.gif)
|
Copyright Xuhua Xia. This is an open-access article distributed under the terms of the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original author and source are credited. The +4G Site in Kozak Consensus Is Not Related to the Efficiency of Translation Initiation Emmanouil Dermitzakis, Academic Editor The Wellcome Trust Sanger Institute, United Kingdom * To whom correspondence should be addressed. E-mail: xxia/at/uottawa.ca Conceived and designed the experiments: XX. Performed the experiments: XX. Analyzed the data: XX. Contributed reagents/materials/analysis tools: XX. Wrote the paper: XX. Received September 7, 2006; Accepted January 5, 2007. | |||||
Abstract The optimal context for translation initiation in mammalian species is GCCRCCaugG (where R ![]() ![]() | |||||
While translation initiation in prokaryotes is mediated by base-paring between the Shine-Dalgarno sequence at the 5-UTR on the mRNA and the anti-Shine-Dalgarno sequence at the 3′-end of the 16S rRNA [1], [2], translation initiation in eukaryotes is mediated by the Kozak consensus [3]–[6]. The optimal context for translation initiation in mammalian species is GCCRCCaugG (where R The interpretation of +4G has been controversial. It has been suggested that +4G may have little to do with initiation site recognition, but is constrained by the requirement for particular type of amino acid residue at the N-terminus of the protein [9]. One piece of supporting evidence came from a detailed study of an influenza virus NS cDNA derivative [10] which showed that both +4 and +5 sites were important and changes at these sites reduced protein production. In contrast, the +6 site (the third codon position of the second codon) is less important. A simple explanation of this result is that changes at the +4 and +5 sites alter the amino acid, whereas those at the +6 site may not. Recent studies, especially those involving the removal of the initiator methionine (Met) and myristoylation, revived the alternative explanation of amino acid constraint for the presence of +4G in protein-coding genes. First, amino-terminal modifications of nascent peptides occur in nearly all proteins in both prokaryotes and eukaryotes, and the removal of the initiator Met, which occurs soon after the amino terminus of the growing polypeptide chain emerges from the ribosome, is not only an important amino-terminal modification in itself, but also required for further amino-terminal modifications. The efficiency of removing the initiator Met depends heavily on the penultimate (the second) amino acid, with the cleavage occurring most efficiently when the penultimate amino acid is small [11]. Alanine (Ala) and glycine (Gly) happen to be the two smallest amino acids and both are coded by G-starting codons, i.e., Ala by the GCN (where N stands for any nucleotide) and Gly by the GGN codons. The need for removing the initiator Met in proteins implies the presence of many Ala and Gly at the penultimate amino acid position and consequently many +4G due to the GCN and GGN codons coding for Ala and Gly, respectively. Another factor contributing to the prevalence of +4G, but independent of the efficiency of translation initiation, is the myristoylation process. For example, in Coxsackievirus B3, the initiation codon is flanked by both −3R and +4G, and viral mutants with a mutation from +4G to +4C is not viable [12]. This may seem to confirm what one would expect based on the necessity of the Kozak consensus for efficient translation initiation in highly expressed genes. However, it turns out that the +4G is required in Coxsackievirus B3 not because it is essential for translation initiation, but because it is needed for coding Gly (coded by GGN). The Gly at the amino terminus, after the removal of the initiator methionine, is needed to attach to a myristoyl (C14H28O2) fatty acid side chain, and myristoylation occurs only on a Gly residue [13]. Myristoylation may involve many proteins, and are implicated in protein subcellular relocalization [13], apoptosis [14], [15], signal transduction [16], [17], and the virulence and colonization of pathogens [12], [18]–[21]. The need for myristoylation in proteins would contribute to the presence of +4G in CDSs. We thus have two alternative hypotheses for the presence of +4G in protein-coding genes. The conventional translation initiation hypothesis argues that the presence of +4G is necessary for highly expressed proteins, with two predictions. First, the selection favoring +4G should drive the increased usage of amino acids coded by GNN codons (e.g., Ala coded by GCN, Asp by GAY, Glu by GAR, Gly by GGN, and Val by GUN) at the penultimate amino acid site. Second, the +4G should be more prevalent in highly expressed than in lowly expressed genes. In contrast, the amino acid constraint hypothesis, based on the amino-terminal modification involving the removal of the initiator Met and myristoylation, has two different predictions. First, not all GNN codons should have increased usage, but only GCN coding Ala and GGN coding Gly should have increased usage. Second, highly expressed genes may need more efficient N-terminal processing and may consequently need more GCN and GGN codons. This may increase the frequency of +4G in highly expressed genes relative to lowly expressed genes. | |||||
Differential use of GNN codons at penultimate site Results from 34169 human coding sequences (CDSs) do not support the translation initiation hypothesis for the presence of +4G. While the five amino acids coded by GNN codons (Ala, Asp, Gly, Glu, Val) account for a majority (64.24%) of the amino acids at the penultimate site (which implies that nucleotide G is the consensus nucleotide at the +4 site), there is no consistent overuse of amino acids coded by GNN codons at the second amino acid site relative to other sites (Fig. 1). This pattern also holds for mouse genes (data not shown).
The expected number of codons (Fig. 1) at the penultimate site is calculated as follows. The total number of codons at non-penultimate sites is 16347992 (excluding the initiation and termination codon). Designate the number of codon XYZ at non-penultimate site as NXYZ. If codon usage at penultimate sites is the same as the rest of the genes, then the expected frequency of codon XYZ at the penultimate site is simply NXYZ*34169/16347992. Only alanine (GCN) codons deviate dramatically from the expected value (Fig. 1).GUN codons (coding for valine) is in fact underused at the penultimate site than at other sites (Fig. 1). Thus, there is no general increase in GNN codon usage at the penultimate site. Differences in +4G frequencies in highly and lowly expressed genes The translation initiation hypothesis also predicts that highly expressed genes should be more likely to have +4G than lowly expressed genes. Ideally we should have genes with different protein expression for testing the prediction. However, there is now substantial evidence suggesting a strong correlation between mRNA level and protein production, not only in Saccharomyces cerevisiae [22]–[25], but also in mammalian species [26]. We used published SAGE (serial analysis of gene expression) data to characterize gene expression because comparative studies have also demonstrated a much higher reproducibility of SAGE (serial analysis of gene expression) experiments than microarray experiments in characterizing mRNA levels [27].To check any possible differences in amino acid usage at the penultimate site and at other sites between highly and lowly expressed genes, I used the 987 unique SAGE tags that were found ubiquitously in human tissues [28]. The main reason for using ubiquitously expressed genes is that the relationship between mRNA level and protein abundance is generally weak in cell-specific genes [29]. These 987 unique tags were matched against the 34169 human CDSs. One gene (ASNA1) matched 3 tags, 16 genes matched 2 tags and 987 genes matched exactly 1 tag (The number happens to be the same as the number of unique tags, but it is accidental). For those 17 multiple-match genes (MMGs), it is difficult to assign gene expression values. For example, if a gene matches two tags, one with 10 copies/cell and another with 100 copies/cell, there is no unequivocal way of assigning an expression value to the gene. For this reason, only 987 single-match genes (SMG) are used and the data set will be referred to as SMG data set. The SMG data set has the problem involving multiple-match tags. For example, if a tag has n copies per cell and matches two genes, say SMG1 and SMG2, it is impossible to know if the n copies/cell of the tag is contributed by SMG 1 only, or SMG 2 only, or by both. In order to assign expression values unequivocally to genes, we also compiled a more limited data set, with on 168 genes that match only single-match tags (SMTs), i.e., each gene matches only one tag which matches only one gene. These genes are designated as SSGs (to reflect the fact that they are from SMG-SMT gene-tag pairs) and their expression values range from 11 copies/cell (gene LRFN4 matching GGGGGGCUGC, excluding the leading 4-bp NlaIII anchoring enzyme site) to 4374 copies/cell (gene GRIN2C matching GGUGACCACG). This small data set will be referred to as the SSG data set. The 168 genes in the SSG data set were divided into a high-expression (HE) group, including 83 SSGs with expression level of at least 50 copies/cell, and a low expression (LE) group, including 85 SSGs with expression level less than 50 copies/cell. The proportion of +4G is 43.37% in the HE group and 49.40% in the LE group. We also contrasted 30 most highly expressed SSGs (expression level at least 114 copies/cell) with 30 least expressed SSGs (expression level equal to 27 copies/cell or less). The proportion of +4G is 43.33% for the former and 46.67% for the latter. Thus, there is no indication that highly expressed genes are more likely to have +4G than lowly expressed genes. The difference is in fact in the opposite direction. This result does not support the prediction of the translation initiation hypothesis. We then categorized the 168 SSGs into five groups according to the codon at the penultimate site (GCN for alanine, GGN for glycine, GAN for aspartate and glutamate, GUN for valine, and genes without +4G designated as NonG), and compared their expression values by one-way analysis of variance (ANOVA). Gene expression differs significantly among the five groups (F We have performed similar ANOVA for the SMG data set and found the same result, i.e., genes with GGN codons at their penultimate sites have higher expression value than the other four groups, with the average gene expression being 178.51 for GCN genes, 263.21 for GGN genes, 163.79 for GAN genes, 175.52 for GUN genes and 145.14 for NonG genes. There is no other significance difference among the five groups. Thus, we may conclude from analyzing the two SAGE data sets that (1) there is no consistent pattern that GNN codons are overused in highly expressed genes, and (2) genes with GGN (glycine) codons at their penultimate site tend to be more highly expressed than other genes. The result is inconsistent with the translation initiation hypothesis but not incompatible with the amino acid constraint hypothesis. An alternative index of gene expression is codon adaptation index, or CAI [31], which has been shown to correlate well with published gene expression in terms of mRNA level and protein abundance [22], [23], [32]. Note that CAI is computed with a codon usage table from a set of reference genes known to be highly expressed. I used the reference set in the Ehum.cut file that is distributed with EMBOSS [33]. However, the cai program in EMBOSS is biased because it does not exclude codon families with a single codon, e.g., AUG coding methionine and UGG coding tryptophan in the standard genetic code (see Materials and Methods for details). I used DAMBE [34,35, version 4.5.10] to calculate CAI values. We focus on two groups of genes, the high-CAI group with CAI>0.8 and the low-CAI group with CAI<0.7. Overall, high-CAI genes tend to have more genes with +4G than low-CAI genes (chi-square test, X2
The use of CAI value [31] as an index of gene expression in human has been controversial [36], [37]. While it is well established that codon usage in bacterial species and vertebrate mitochondria is strongly constrained by the relative tRNA abundance and that there is selection pressure favoring codon-anticodon adaptation [38]–[41], there is only limited evidence for eukaryotes [42]. One additional piece of evidence supporting codon-anticodon adaptation is that the codon frequencies of the 34169 annotated human CDSs are positively correlated with the copy number of their cognate tRNA genes found at The Genomic tRNA Database (http://lowelab.ucsc.edu/GtRNAdb/), compiled with the tRNAscan-SE program [43]. For example, from the 505 human tRNA genes decoding the regular set of 20 amino acids, one can obtain their cognate codon frequencies from their anticodons. These tRNA-derived cognate codon frequencies correlate positively with the codon frequencies of the 168 genes in the SSG data set (Pearson r | |||||
The original study documenting the importance of +4G [44] does not constitute a sufficient proof that +4G is important in translation initiation. The study was based on the production of proinsulin from preproinsulin. The latter has a signal peptide at its amino terminus. The signal peptide is removed during translation, generating proinsulin. When +4G is mutated to other nucleotide, the production of proinsulin is reduced. This reduced proinsulin production was assumed to be caused by reduced efficiency in translation initiation due to the mutation of +4G to other nucleotides. However, one should note that altering +4G also alter the amino acid sequence in the signal peptide and may consequently affect the removal of the signal peptide, leading to reduced production of proinsulin. Thus, the result is compatible with the amino acid constraint hypothesis for the presence of +4G. One may think that the translation initiation hypothesis is partially correct because, after all, two of the five amino acids with G-starting codons (especially alanine) show increased usage. This is wrong. The increased usage of alanine exists not only in eukaryotes, but also in prokaryotes [45] that do not use the scanning mechanism for translation initiation and consequently do not need the +4G. The overuse of small amino acids at the penultimate amino acid site in both prokaryotes and eukaryotes is better explained by the necessity for removing the initiator Met. While the results above are inconsistent with the predictions of the translation initiation hypothesis, they generally appear to support the amino acid constraint hypothesis. First, the latter predicted the overuse of Ala and Gly at the second codon position (to facilitate the removal of the initiator Met and myristoylation), and Ala and Gly (especially Ala) are indeed overused (Table 1). Second, all differences between highly expressed and lowly expressed genes involve GCN and GGN (coding for alanine and glycine, respectively) codons. Both amino acid constraint hypothesis and the translation initiation hypothesis have difficulties in explain certain observations. For example, Kozak [6] found that +4G generally enhances translation initiation, but does not when it occurs in a GUN codon (coding for valine). An associated finding in this paper is that GUN is underused at the penultimate site (Fig. 1). While such findings are difficult to explain by the translation initiation hypothesis, they are also difficult for the amino acid constraint hypothesis unless valine at the penultimate site reduces the efficiency of initiator Met cleavage. Previous studies on prokaryotes and eukaryotes [46] and on the yeast, Saccharomyces cerevisiae [11], suggest that cleavage of initiator Met occurs with valine at the penultimate site, but a recent study on Escherichia coli [47] demonstrates that peptides with Val at the penultimate site dramatically reduces the efficiency of initiator Met cleavage relative to other amino acids such as Ala, Cys, Gly, Pro, or Ser in this position. Further studies on initiator Met cleavage on mammalian species are needed before one can reach a solid conclusion. In summary, we conclude that the presence of +4G is poorly explained by the translation initiation hypothesis that claims the necessity of +4G for efficient translation initiation, but well explained by the alternative amino acid constraint hypothesis that claims the necessity of Ala and Gly at the second amino acid position in many proteins (for the removal of the initiator Met or myristoylation) as the cause of the prevalence of +4G because Ala and Gly happen to be coded by GCN and GGN codons. The necessity of +4G for efficient translation initiation appears to be a misconception that has existed in the molecular biology textbooks for too long. | |||||
I retrieved the rna.gbk.gz file at ftp://ftp.ncbi.nih.gov/genomes/H_sapiens/RNA/, dated Sept. 3, 2006, and extracted all 34169 annotated coding sequences (CDSs) for evaluating the translation initiation hypothesis and the amino acid constraint hypothesis. CDS extraction and computation of codon adaptation index [31], as well as the analysis of codon usage at the second codon were carried out by using DAMBE (Xia 2001; Xia and Xie 2001). CDSs that are not multiples of three or not terminated with a stop codon are excluded. Because the translation initiation hypothesis predicts that highly expressed genes should be more likely to have +4G than lowly expressed genes, we have used genes of different expression levels to check this prediction. Gene expression level is measured in two ways in this study, one by SAGE (serial analysis of gene expression) data and one by using codon adaptation index. The use of SAGE data instead of available microarray data is mainly because of the much higher reproducibility of SAGE experiments than microarray experiments [27]. SAGE data were retrieved from http://www.nature.com/ng/journal/v23/n4/extref/ng1299-387b-S1.pdf which listed, among others, 987 unique tags that are ubiquitously expressed in different human tissues, together with their abundance in copies/cell. We searched these tags against the 34169 human CDSs for exact matches. To facilitate presentation, we define multiple-match genes (MMGs) as those CDSs each matching multiple tags, single-match genes (SMGs) as those each matching a single tag, multiple-match tags (MMTs) as those tags each matching multiple CDSs, and single-match tags (SMTs) as those each matching a single CDS. It is difficult to assess the expression level of an MMG because different tags it matches have different copies/cell values. The presence of MMTs causes an even more serious problem. For example, when a tag with 100 copies/cell matches two genes, it is impossible to know if the 100 copies are contributed by only one gene or by both. It would be methodologically wrong to assign both genes an expression value of 100 copies/cell. For this reason, we have compiled two data sets, one including all SMGs and the other including only SMGs that match SMTs (i.e., a SMG that matches a MMT is not included). CAI for a gene is computed from (1) the codon frequencies of the gene and (2) a codon usage table from a set of reference genes known to be highly expressed, according to the following equation [31]:
It is important to exclude codon families with a single codon. Note that for such codons (e.g., AUG and UGG in the standard genetic code), their corresponding wi value will always be 1 regardless of codon usage bias of the gene. If a gene happens to use a high proportion of methionine and tryptophan, then it will have a high CAI value even if the codon usage is not at all biased. The cai program in EMBOSS [33] does not exclude codon families with a single codon because the CAI values from that program are the same as those I computed without excluding the AUG and UGG codons. I used DAMBE [34,35, version 4.5.10] which excludes AUG and UGG in computing CAI. The reason for using CAI as an index of gene expression instead of taking advantage of the availability of gene expression data is that, in higher eukaryotes such as human, many genes are highly expressed only at specific time and in specific tissues. For this reason, a gene with no detectable expression in a specific study, which typically involves few time points and few tissues, should not be taken as a low expressed gene. However, the availability of the gene expression data does vindicate the use of CAI as a general measure of gene expression [32]. Because CAI is based on codon frequencies of a gene with respect to the codon usage of a reference set of genes known to be highly expressed, short sequences with few codons may produce unreliable CAI values. For this reason, two separate analyses were performed, one with all CDSs and the other excluding CDSs shorter than 300 bp. The two sets of results are almost identical because short CDSs constitute only a small fraction. | |||||
Acknowledgments I thank S. Aris-Brosou, M. Carullo, S. Khalouei, P. Ma, J. Mennigen, A. Meyer, Z. Song, H. Xiong, and X. Yao for comments and discussion. An anonymous reviewer and M. Dermitzakis provided many helpful suggestions. | |||||
Footnotes Competing Interests: The author has declared that no competing interests exist. Funding: This project is supported by the Discovery, Strategic Research, and Research Tools and Instruments Grants from the Natural Science and Engineering Research Council of Canada. | |||||
References 1. Shine, J; Dalgarno, L. Terminal-sequence analysis of bacterial ribosomal RNA. Correlation between the 3′-terminal-polypyrimidine sequence of 16-S RNA and translational specificity of the ribosome. Eur J Biochem. 1975;57:221–230. [PubMed] 2. Shine, J; Dalgarno, L. Determinant of cistron specificity in bacterial ribosomes. Nature. 1975;254:34–38. [PubMed] 3. Kozak, M. Initiation of translation in prokaryotes and eukaryotes. Gene. 1999;234:187–208. [PubMed] 4. Kozak, M. Effects of long 5′ leader sequences on initiation by eukaryotic ribosomes in vitro. Gene Expr. 1991;1:117–125. [PubMed] 5. Kozak, M. Possible role of flanking nucleotides in recognition of the AUG initiator codon by eukaryotic ribosomes. Nucleic Acids Res. 1981;9:5233–5252. [PubMed] 6. Kozak, M. Recognition of AUG and alternative initiator codons is augmented by G in position +4 but is not generally affected by the nucleotides in positions +5 and +6. Embo J. 1997;16:2482–2492. [PubMed] 7. Kozak, M. Not every polymorphism close to the AUG codon can be explained by invoking context effects on initiation of translation. Blood. 2003;101:1202–1203; author reply 1203. [PubMed] 8. Kozak, M. A progress report on translational control in eukaryotes. Sci STKE. 2001;2001:PE1. [PubMed] 9. Cigan, AM; Donahue, TF. Sequence and structural features associated with translational initiator regions in yeast–a review. Gene. 1987;59:1–18. [PubMed] 10. Grunert, S; Jackson, RJ. The immediate downstream codon strongly influences the efficiency of utilization of eukaryotic translation initiation codons. Embo J. 1994;13:3618–3630. [PubMed] 11. Moerschell, RP; Hosokawa, Y; Tsunasawa, S; Sherman, F. The specificities of yeast methionine aminopeptidase and acetylation of amino-terminal methionine in vivo. Processing of altered iso-1-cytochromes c created by oligonucleotide transformation. J Biol Chem. 1990;265:19638–19643. [PubMed] 12. Harkins, S; Cornell, CT; Whitton, JL. Analysis of translational initiation in coxsackievirus B3 suggests an alternative explanation for the high frequency of R+4 in the eukaryotic consensus motif. J Virol. 2005;79:987–996. [PubMed] 13. Farazi, TA; Waksman, G; Gordon, JI. The Biology and Enzymology of Protein N-Myristoylation. J Biol Chem. 2001;276:39501–39504. [PubMed] 14. Vilas, GL; Corvi, MM; Plummer, GJ; Seime, AM; Lambkin, GR, et al. Posttranslational myristoylation of caspase-activated p21-activated protein kinase 2 (PAK2) potentiates late apoptotic events. Proc Natl Acad Sci U S A. 2006;103:6542–6547. [PubMed] 15. Sakurai, N; Utsumi, T. Posttranslational N-myristoylation is required for the anti-apoptotic activity of human tGelsolin, the C-terminal caspase cleavage product of human gelsolin. J Biol Chem. 2006;281:14288–14295. [PubMed] 16. Rowe, DC; McGettrick, AF; Latz, E; Monks, BG; Gay, NJ, et al. The myristoylation of TRIF-related adaptor molecule is essential for Toll-like receptor 4 signal transduction. Proc Natl Acad Sci U S A. 2006;103:6299–6304. [PubMed] 17. de Vries, JS; Andriotis, VM; Wu, AJ; Rathjen, JP. Tomato Pto encodes a functional N-myristoylation motif that is required for signal transduction in Nicotiana benthamiana. Plant J. 2006;45:31–45. [PubMed] 18. Robert-Seilaniantz, A; Shan, L; Zhou, JM; Tang, X. The pseudomonas syringae pv. tomato DC3000 Type III effector HopF2 has a putative myristoylation site required for its avirulence and virulence functions. Mol Plant Microbe Interact. 2006;19:130–138. [PubMed] 19. Bentham, M; Mazaleyrat, S; Harris, M. Role of myristoylation and N-terminal basic residues in membrane association of the human immunodeficiency virus type 1 Nef protein. J Gen Virol. 2006;87:563–571. [PubMed] 20. Breuer, S; Gerlach, H; Kolaric, B; Urbanke, C; Opitz, N, et al. Biochemical indication for myristoylation-dependent conformational changes in HIV-1 Nef. Biochemistry. 2006;45:2339–2349. [PubMed] 21. Provitera, P; El-Maghrabi, R; Scarlata, S. The effect of HIV-1 Gag myristoylation on membrane binding. Biophys Chem. 2006;119:23–32. [PubMed] 22. Futcher, B; Latter, GI; Monardo, P; McLaughlin, CS; Garrels, JI. A sampling of the yeast proteome. Mol Cell Biol. 1999;19:7357–7368. [PubMed] 23. Gygi, SP; Rochon, Y; Franza, BR; Aebersold, R. Correlation between protein and mRNA abundance in yeast. Mol Cell Biol. 1999;19:1720–1730. [PubMed] 24. Ideker, T; Thorsson, V; Ranish, JA; Christmas, R; Buhler, J, et al. Integrated genomic and proteomic analyses of a systematically perturbed metabolic network. Science. 2001;292:929–934. [PubMed] 25. Griffin, TJ; Gygi, SP; Ideker, T; Rist, B; Eng, J, et al. Complementary profiling of gene expression at the transcriptome and proteome levels in Saccharomyces cerevisiae. Mol Cell Proteomics. 2002;1:323–333. [PubMed] 26. Tian, Q; Stepaniants, SB; Mao, M; Weng, L; Feetham, MC, et al. Integrated Genomic and Proteomic Analyses of Gene Expression in Mammalian Cells. Mol Cell Proteomics. 2004;3:960–969. [PubMed] 27. Dinel, S; Bolduc, C; Belleau, P; Boivin, A; Yoshioka, M, et al. Reproducibility, bioinformatic analysis and power of the SAGE method to evaluate changes in transcriptome. Nucleic Acids Res. 2005;33:e26. [PubMed] 28. Velculescu, VE; Madden, SL; Zhang, L; Lash, AE; Yu, J, et al. Analysis of human transcriptomes. Nat Genet. 1999;23:387–388. [PubMed] 29. Chen, G; Gharib, TG; Huang, CC; Taylor, JM; Misek, DE, et al. Discordant protein and mRNA expression in lung adenocarcinomas. Mol Cell Proteomics. 2002;1:304–313. [PubMed] 30. Zar, JH. Biostatistical Aanalysis. Upper Saddle River, New Jersey: Prentice Hall; 1999. 31. 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] 32. Coghlan, A; Wolfe, KH. Relationship of codon bias to mRNA concentration and protein length in Saccharomyces cerevisiae. Yeast. 2000;16:1131–1145. [PubMed] 33. Rice, P; Longden, I; Bleasby, A. EMBOSS: the European Molecular Biology Open Software Suite. Trends Genet. 2000;16:276–277. [PubMed] 34. Xia, X; Xie, Z. DAMBE: Software package for data analysis in molecular biology and evolution. Journal of Heredity. 2001;92:371–373. [PubMed] 35. Xia, X. Data analysis in molecular biology and evolution. Boston: Kluwer Academic Publishers; 2001. p. 277 p. 36. Duret, L. Evolution of synonymous codon usage in metazoans. Current Opinion in Genetics & Development. 2002;12:640–649. [PubMed] 37. Semon, M; Lobry, JR; Duret, L. No evidence for tissue-specific adaptation of synonymous codon usage in humans. Mol Biol Evol. 2006;23:523–529. [PubMed] 38. Gouy, M; Gautier, C. Codon usage in bacteria: correlation with gene expressivity. Nucleic Acids Res. 1982;10:7055–7064. [PubMed] 39. 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. [PubMed] 40. Xia, X. Mutation and selection on the anticodon of tRNA genes in vertebrate mitochondrial genomes. Gene. 2005;345:13–20. [PubMed] 41. Xia, X. How optimized is the translational machinery in Escherichia coli, Salmonella typhimurium and Saccharomyces cerevisiae? Genetics. 1998;149:37–44. [PubMed] 42. Plotkin, JB; Robins, H; Levine, AJ. Tissue-specific codon usage and the expression of human genes. Proc Natl Acad Sci U S A. 2004;101:12588–12591. [PubMed] 43. Lowe, TM; Eddy, SR. tRNAscan-SE: a program for improved detection of transfer RNA genes in genomic sequence. Nucleic Acids Res. 1997;25:955–964. [PubMed] 44. Kozak, M. Point mutations define a sequence flanking the AUG initiator codon that modulates translation by eukaryotic ribosomes. Cell. 1986;44:283–292. [PubMed] 45. Tats, A; Remm, M; Tenson, T. Highly expressed proteins have an increased frequency of alanine in the second amino acid position. BMC Genomics. 2006;7:28. [PubMed] 46. Flinta, C; Persson, B; Jornvall, H; von Heijne, G. Sequence determinants of cytosolic N-terminal protein processing. Eur J Biochem. 1986;154:193–196. [PubMed] 47. Frottin, F; Martinez, A; Peynot, P; Mitra, S; Holz, RC, et al. The Proteomics of N-terminal Methionine Cleavage. Mol Cell Proteomics. 2006;5:2336–2349. [PubMed] | |||||