Efficient translation initiation dictates codon usage at gene start

The genetic code is degenerate; thus, protein evolution does not uniquely determine the coding sequence. One of the puzzles in evolutionary genetics is therefore to uncover evolutionary driving forces that result in specific codon choice. In many bacteria, the first 5–10 codons of protein‐coding genes are often codons that are less frequently used in the rest of the genome, an effect that has been argued to arise from selection for slowed early elongation to reduce ribosome traffic jams. However, genome analysis across many species has demonstrated that the region shows reduced mRNA folding consistent with pressure for efficient translation initiation. This raises the possibility that unusual codon usage is a side effect of selection for reduced mRNA structure. Here we discriminate between these two competing hypotheses, and show that in bacteria selection favours codons that reduce mRNA folding around the translation start, regardless of whether these codons are frequent or rare. Experiments confirm that primarily mRNA structure, and not codon usage, at the beginning of genes determines the translation rate.     

Selection Intensity for Codon Bias

The patterns of nonrandom usage of synonymous codons (codon bias) in enteric bacteria were analyzed. Poisson random field (PRF) theory was used to derive the expected distribution of frequencies of nucleotides differing from the ancestral state at aligned sites in a set of DNA sequences. This distribution was applied to synonymous nucleotide polymorphisms and amino acid polymorphisms in the gnd and putP genes of Escherichia coli. For the gnd gene, the average intensity of selection against disfavored synonymous codons was estimated as approximately 7.3 X 10(-9); this value is significantly smaller than the estimated selection intensity against selectively disfavored amino acids in observed polymorphisms (2.0 X 10(-8)), but it is approximately of the same order of magnitude. The selection coefficients for optimal synonymous codons estimated from PRF theory were consistent with independent estimates based on codon usage for threonine and glycine. Across 118 genes in E. coli and Salmonella typhimurium, the distribution of estimated selection coefficients, expressed as multiples of the effective population size, has a mean and standard deviation of 0.5 +/- 0.4. No significant differences were found in the degree of codon bias between conserved positions and replacement positions, suggesting that translational misincorporation is not an important selective constraint among synonymous polymorphic codons in enteric bacteria. However, across the first 100 codons of the genes, conserved amino acids with identical codons have significantly greater codon bias than of either synonymous or nonidentical codons, suggesting that there are unique selective constraints, perhaps including mRNA secondary structures, in this part of the coding region.     

Absence of effect of varying Thr-Leu codon pairs on protein synthesis in a T7 system

The over-represented threonine-leucine (Thr-Leu) codon pair ACG CUG has been previously reported to be inhibitory to translation compared to the synonymous under-represented Thr-Leu codon pair ACC CUG, in an E. coli system in which the codon pairs were located either 3 and 4, or 6 and 7, or 9 and 10 codons downstream from the initiating codon for the message [Irwin, B., Heck, J. D., and Hatfield, G. W. (1995) J. Biol. Chem. 270, 22801-22806]. In the work reported here, these synonymous codon pairs were tested in a T7 system, with the codon pairs located either 14 and 15, or 6 and 7 codons downstream from the AUG start codon. In contrast to the reported findings in the E. coli system, there was no difference found in translation between mRNAs containing the respective codon pairs in the T7 system. The reasons for the different findings remain unclear, but presumably are a consequence of differences between the E. coli and T7 systems used to assay gene expression. Nevertheless, as a result of this work, it appears that the effect of varying codon pairs reported in the E. coli system is not due to a difference in translational step times through the respective codon pairs, as previously proposed.     

Evolution of prokaryotic DNA: intragenic and extragenic divergences observed with orthologs from three related species

This study compared orthologous gene pairs from Escherichia coli K12, E. coli O157:H7 EDL933, Salmonella typhimurium LT2, and Yersinia pestis CO92 using only homologs of equal length, and differing nucleotides were counted and mapped. The data showed very clearly how the rates of divergence change with intragenic and extragenic position. The rate of synonymous mutation is lowest near start codons and near stop codons, and, a little surprisingly, the opposite seemed to be true for nonsynonymous substitutions. Analysis outside genes reveals that nucleotide divergences occur less frequently upstream of start codons than downstream of stop codons, and a remarkable drop in divergences is seen for two of the data sets around N = 9 nucleotides upstream of start codons; that is, the Shine-Dalgarno region changes at a lower rate. The explanation is likely to be the link between expressivity and sequence complementarity to the 3' end of 16S ribosomal rRNA. The latter is highly conserved across many bacterial and archaebacterial species.     

The evolution of biased codon and amino acid usage in nematode genomes

Despite the degeneracy of the genetic code, whereby different codons encode the same amino acid, alternative codons and amino acids are utilized nonrandomly within and between genomes. Such biases in codon and amino acid usage have been demonstrated extensively in prokaryote genomes and likely reflect a balance between the action of mutation, selection, and genetic drift. Here, we quantify the effects of selection and mutation drift as causes of codon and amino acid-usage bias in a large collection of nematode partial genomes from 37 species spanning approximately 700 Myr of evolution, as inferred from expressed sequence tag (EST) measures of gene expression and from base composition variation. Average G + C content at silent sites among these taxa ranges from 10% to 63%, and EST counts range more than 100-fold, underlying marked differences between the identities of major codons and optimal codons for a given species as well as influencing patterns of amino acid abundance among taxa. Few species in our sample demonstrate a dominant role of selection in shaping intragenomic codon-usage biases, and these are principally free living rather than parasitic nematodes. This suggests that deviations in effective population size among species, with small effective sizes among parasites, are partly responsible for species differences in the extent to which selection shapes patterns of codon usage. Nevertheless, a consensus set of optimal codons emerges that is common to most taxa, indicating that, with some notable exceptions, selection for translational efficiency and accuracy favors similar sets of codons regardless of the major codon-usage trends defined by base compositional properties of individual nematode genomes.     

Synonymous codon ordering: a subtle but prevalent strategy of bacteria to improve translational efficiency

Background

In yeast coding sequences, once a particular codon has been used, subsequent occurrence of the same amino acid tends to use codons sharing the same tRNA. Such a phenomenon of co-tRNA codons pairing bias (CTCPB) is also found in some other eukaryotes but it is not known whether it occurs in prokaryotes.

Methodology/Principal Findings

In this study, we focused on a total of 773 bacterial genomes to investigate their synonymous codon pairing preferences. After calculating the actual frequencies of synonymous codon pairs and comparing them with their expected values, we detected an obvious pairing bias towards identical codon pairs. This seems consistent with the previously reported CTCPB phenomenon, since identical codons are certainly read by the same tRNA. However, among co-tRNA but non-identical codon pairs, only 22 were often found overrepresented, suggesting that many co-tRNA codons actually do not preferentially pair together in prokaryotes. Therefore, the previously reported co-tRNA codons pairing rule needs to be more rigorously defined. The affinity differences between a tRNA anticodon and its readable codons should be taken into account. Moreover, both within-gene-shuffling tests and phylogenetic analyses support the idea that translational selection played an important role in shaping the observed synonymous codon pairing pattern in prokaryotes.

Conclusions

Overall, a high level of synonymous codon pairing bias was detected in 73% investigated bacterial species, suggesting the synonymous codon ordering strategy has been prevalently adopted by prokaryotes to improve their translational efficiencies. The findings in this study also provide important clues to better understand the complex dynamics of translational process.     

Rare codon clusters at 5'-end influence heterologous expression of archaeal gene in Escherichia coli

Proteins from hyperthermophilic microorganisms are attractive candidates for novel biocatalysts because of their high resistance to temperature extremes. However, archaeal genes are usually poorly expressed in Escherichia coli because of differences in codon usage. Genes from the thermoacidophilic archaea Sulfolobus solfataricus and Thermoplasma acidophilum contain high proportions of rare codons for arginine, isoleucine, and leucine, which are recognized by the tRNAs encoded by the argU, ileY, and leuW genes, respectively, and which are rarely used in E. coli. To examine the effects of these rare codons on heterologous expression, we expressed the Sso_gnaD and Tac_gnaD genes from S. solfataricus and T. acidophilum, respectively, in E. coli. The Sso_gnaD product was expressed at very low levels when the open reading frame (ORF) was cloned in pRSET and expressed in E. coli BL21(DE3), and was expressed at much higher levels in the E. coli BL21(DE3)-CodonPlus RIL strain, which contains extra copies of the argU, ileY, and leuW tRNA genes. In contrast, Tac_gnaD was expressed at similar levels in both E. coli strains. Comparison of the Sso_gnaD and Tac_gnaD gene sequences revealed that the 5'-end of the Sso_gnaD sequence was rich in AGA(arg) and ATA(Ile) codons. These codons were replaced with the codons commonly used in E. coli by polymerase chain reaction-mediated site-directed mutagenesis. The results of expression studies showed that a non-tandem repeat of rare codons is critical in the observed interference in heterologous expression of this gene. We concluded that the level of heterologous expression of Sso_gnaD in E. coli was limited by the clustering of the rare codons in the ORF, rather than on the rare codon frequency.     

Activity of Chi Recombinational Hotspots in SALMONELLA TYPHIMURIUM

Chi sites have previously been shown to stimulate homologous recombination by the Escherichia coli RecBC pathway. To test the activity of Chi in another organism, bacteriophage lambda crosses were carried out in Salmonella typhimurium strains bearing the E. coli lambda receptor protein. Chi is active in these crosses in S. typhimurium, but is less active than in the same crosses carried out in E. coli. The lower Chi activity in S. typhimurium appears to be intrinsic to the S. typhimurium RecBC enzyme, since the Chi activity in E. coli-S. typhimurium hybrids depends on the species of origin of their RecBC enzyme. For these studies we constructed and F' factor and a pBR322-derived plasmid carrying the thyA+ recC+ recB+ argA+ region of the S. typhimurium chromosome.     

Preferential usage of some minor codons in bacteria

In many bacterial species, such as Deinococcus radiodurans, Haemophilus influenzae, and Methanobacterium thermoautotrophicum, some minor codons are preferentially used near the initiation codon. Among these codons, there are some minor codons that have strong preference for the initiation site in the high codon adaptation index (CAI) group (comprising of highly expressed genes) rather than in the low CAI genes group (comprising of low expressing genes). In the present study, codon usage in the initiation site and in the rest of the gene was systematically compared in the 27 complete bacterial genomes and Saccharomyces cerevisiae genome. Furthermore, we classified genes into two groups according to the CAI values and conducted the same analysis for each of the two groups. Our results suggest a role for some minor codons in the initiation site of the regulating translation system in many bacteria. We have summarized codons that are preferentially used in the initiation site and probably play a role in regulating genes expression in these organisms.     

Efficient expression of gene variants that harbour AGA codons next to the initiation codon

In an effort to improve the knowledge about the rules which direct the effect of the early ORF sequences on translation efficiency, we have analyzed the effect of pairs of the six arginine codons at the second and third positions on the expression of lacZ variants. Whereas the pairs of identical AGA or AGG codons were favorable for the gene expression, identical pairs of each of the four CGN codons were very inefficient. This result was unexpected because tandems of AGA or AGG codons located in more internal gene positions provoke deficient expression whilst internally located CGU and CGC are the most abundant and efficiently translated arginine codons. The mixed combinations of AGA and each of the CGN codons usually resulted in efficient rates of lacZ expression independently of the peptidyl-tRNA propensity to dissociate from the ribosome. Thus, the variant harboring the pair of AGA codons was expressed as efficiently as the variant carrying a pair of AAA codons in the same positions, a configuration reported as one of the most common and efficient for gene expression. We explain these results assuming that the presence of adenines in these early positions enhance gene expression. As expected, specific mRNA levels correlated with the intensity of lacZ expression for each variant. However, the induction of lacZ AGA AGA gene in pth cells accumulated peptidyl-tRNA^Arg4 as well as a short 5′-proximal lacZ mRNA fragment suggesting ribosome stalling due to depletion of aminoacylated-tRNA^Arg4.     

Lactic acid bacteria as prime candidates for codon optimization

In species having a strong correlation of expressivity and codon bias it has been shown that heterologous expression can be optimized by changing codons of the introduced gene towards the set of codons that the host organism naturally uses in its highly expressed genes. Even though two lactic acid bacteria are fully sequenced, there are no reports on attempts of codon optimization in the literature. In this report it is demonstrated that codons used in highly expressed genes tend to differ from the codons in lowly expressed genes, and that there is a strong correlation of codon bias and empirical expressivity (codon adaptation index) in Lactococcus lactis and Lactobacillus plantarum. This strongly suggests that codon optimization strategies could be applied to expression systems with lactic acid bacteria as producer strains. A good example of a candidate for codon optimization is the mouse interleukin-2 gene, which in its natural form has an extremely low codon adaptation index for expression in Lc. lactis.     

Analysis of synonymous codon usage bias in Chlamydia

Chlamydiae are obligate intracellular bacterial pathogens that cause ocular and sexuallytransmitted diseases,and are associated with cardiovascular diseases.The analysis of codon usage mayimprove our understanding of the evolution and pathogenesis of Chlamydia and allow reengineering of targetgenes to improve their expression for gene therapy.Here,we analyzed the codon usage of C.muridarum,C.trachomatis(here indicating biovar trachoma and LGV),C.pneumoniae,and C.psittaci using the codonusage database and the CUSP(Create a codon usage table)program of EMBOSS(The European MolecularBiology Open Software Suite).The results show that the four genomes have similar codon usage patterns,with a strong bias towards the codons with A and T at the third codon position.Compared with Homosapiens,the four chlamydial species show discordant seven or eight preferred codons.The ENC(effectivenumber of codons used in a gene)-plot reveals that the genetic heterogeneity in Chlamydia is constrained bythe G+C content,while translational selection and gene length exert relatively weaker influences.Moreover,mutational pressure appears to be the major determinant of the codon usage variation among the chlamydialgenes.In addition,we compared the codon preferences of C.trachomatis with those of E.coli,yeast,adenovirus and Homo sapiens.There are 23 codons showing distinct usage differences between C.trachomatisand E.coli,24 between C.trachomatis and adenovirus,21 between C.trachomatis and Homo sapiens,butonly six codons between C.trachomatis and yeast.Therefore,the yeast system may be more suitable for theexpression of chlamydial genes.Finally,we compared the codon preferences of C.trachomatis with those ofsix eukaryotes,eight prokaryotes and 23 viruses.There is a strong positive correlation between the differ-ences in coding GC content and the variations in codon bias(r=0.905,P<0,001).We conclude that thevariation of codon bias between C.trachomatis and other organisms is much less influenced by phylogeneticlineage and primarily determined by the extent of disparities in GC content.     

Escherichia coli translation initiation factor 3 discriminates the initiation codon in vivo

In a genetic selection designed to isolate Escherichia coli mutations that increase expression of the IS 10 transposase gene ( tnp ), we unexpectedly obtained viable mutants defective in translation initiation factor 3 (IF3). Several lines of evidence led us to conclude that transposase expression, per se , was not increased. Rather, these mutations appear to increase expression of the tnp'–'lacZ gene fusions used in this screen, by increasing translation initiation at downstream, atypical initiation codons. To test this hypothesis we undertook a systematic analysis of start codon requirements and measured the effects of IF3 mutations on initiation from various start codons. Beginning with an efficient translation initiation site, we varied the AUG start codon to all possible codons that differed from AUG by one nucleotide. These potential start codons fall into distinct classes with regard to translation efficiency in vivo : Class I codons (AUG, GUG, and UUG) support efficient translation; Class IIA codons (CUG, AUU, AUC, AUA, and ACG) support translation at levels only 1–3% that of AUG; and Class IIB codons (AGG and AAG) permit levels of translation too low for reliable quantification. Importantly, the IF3 mutations had no effect on translation from Class I codons, but they increased translation from Class II codons 3–5-fold, and this same effect was seen in other gene contexts. Therefore, IF3 is generally able to discriminate between efficient and inefficient codons in vivo , consistent with earlier in vitro observations. We discuss these observations as they relate to IF3 autoregulation and the mechanism of IF3 function.     

Codon usage bias covaries with expression breadth and the rate of synonymous evolution in humans, but this is not evidence for selection.

In numerous species, from bacteria to Drosophila, evidence suggests that selection acts even on synonymous codon usage: codon bias is greater in more abundantly expressed genes, the rate of synonymous evolution is lower in genes with greater codon bias, and there is consistency between genes in the same species in which codons are preferred. In contrast, in mammals, while nonequal use of alternative codons is observed, the bias is attributed to the background variance in nucleotide concentrations, reflected in the similar nucleotide composition of flanking noncoding and exonic third sites. However, a systematic examination of the covariants of codon usage controlling for background nucleotide content has yet to be performed. Here we present a new method to measure codon bias that corrects for background nucleotide content and apply this to 2396 human genes. Nearly all (99%) exhibit a higher amount of codon bias than expected by chance. The patterns associated with selectively driven codon bias are weakly recovered: Broadly expressed genes have a higher level of bias than do tissue-specific genes, the bias is higher for genes with lower rates of synonymous substitutions, and certain codons are repeatedly preferred. However, while these patterns are suggestive, the first two patterns appear to be methodological artifacts. The last pattern reflects in part biases in usage of nucleotide pairs. We conclude that we find no evidence for selection on codon usage in humans.     

Functional significance of an evolutionarily conserved alanine (GCA) resume codon in tmRNA in Escherichia coli

Occasionally, ribosomes stall on mRNAs prior to the completion of the polypeptide chain. In Escherichia coli and other eubacteria, tmRNA-mediated trans-translation is a major mechanism that recycles the stalled ribosomes. The tmRNA possesses a tRNA-like domain and a short mRNA region encoding a short peptide (ANDENYALAA in E. coli) followed by a termination codon. The first amino acid (Ala) of this peptide encoded by the resume codon (GCN) is highly conserved in tmRNAs in different species. However, reasons for the high evolutionary conservation of the resume codon identity have remained unclear. In this study, we show that changing the E. coli tmRNA resume codon to other efficiently translatable codons retains efficient functioning of the tmRNA. However, when the resume codon was replaced with the low-usage codons, its function was adversely affected. Interestingly, expression of tRNAs decoding the low-usage codon from plasmid-borne gene copies restored efficient utilization of tmRNA. We discuss why in E. coli, the GCA (Ala) is one of the best codons and why all codons in the short mRNA of the tmRNA are decoded by the abundant tRNAs.     

Chemotactic transducer proteins of Escherichia coli exhibit homology with methyl-accepting proteins from distantly related bacteria.

Transducers are transmembrane, methyl-accepting proteins central to the chemotactic systems of the enteric bacteria Escherichia coli and Salmonella typhimurium. Methyl-accepting proteins have been reported in a number of species in addition to these enteric bacteria. Those species include Bacillus subtilis and Spirochaeta aurantia, representatives of groups that diverged from ancestral enteric bacteria and from each other very early in bacterial evolution. An antiserum that reacts with all transducers of E. coli precipitated specifically methyl-accepting proteins from B. subtilis and S. aurantia, indicating that these proteins share antigenic determinants with transducers of E. coli. In addition, analysis of tryptic peptides by high-pressure liquid chromatography revealed similarities in the regions of methyl-accepting sites for proteins from all three species. These observations imply that structural features have been preserved in the three species from transducers contained in a common ancestor of eubacteria. It is thus reasonable to predict that other flagellated, chemotactic bacteria will be found to contain methyl-accepting proteins homologous to transducers of enteric bacteria.     

Consequences of stop codon reassignment on protein evolution in ciliates with alternative genetic codes

Tetrahymena thermophila and Paramecium tetraurelia are ciliates that reassign TAA and TAG from stop codons to glutamine codons. Because of the lack of full genome sequences, few studies have concentrated on analyzing the effects of codon reassignment in protein evolution. We used the recently sequenced genome of these species to analyze the patterns of amino acid substitution in ciliates that reassign the code. We show that, as expected, the codon reassignment has a large impact on amino acid substitutions in closely related proteins; however, contrary to expectations, these effects also hold for very diverged proteins. Previous studies have used amino acid substitution data to calculate the minimization of the genetic code; our results show that because of the lasting influence of the code in the patterns of substitution, such studies are tautological. These different substitution patterns might affect alignment of ciliate proteins, as alignment programs use scoring matrices based on substitution patterns of organisms that use the standard code. We also show that glutamine is used more frequently in ciliates than in other species, as often as expected based on the presence of the 2 new reassigned codons, indicating that the frequencies of amino acids in proteomes is mostly determined by neutral processes based on their number of codons.     

Optimal codon identities in bacteria: implications from the conflicting results of two different methods

A correlation method was recently adopted to identify selection-favored 'optimal' codons from 675 bacterial genomes. Surprisingly, the identities of these optimal codons were found to track the bacterial GC content, leading to a conclusion that selection would generally shape the codon usages to the same direction as the overall mutation does. Raising several concerns, here we report a thorough comparative study on 203 well-selected bacterial species, which strongly suggest that the previous conclusion is likely an illusion. Firstly, the previous study did not preclude species that are suffering weak or no selection pressures on their codon usages. For these species, as showed in this study, the optimal codon identities are prone to be incorrect and follow GC content. Secondly, the previous study only adopted the correlation method, without considering another method to test the reliability of inferred optimal codons. Actually by definition, optimal codons can also be identified by simply comparing codon usages between high- and low-expression genes. After using both methods to identify optimal codons for the selected species, we obtained highly conflicting results, suggesting at least one method is misleading. Further we found a critical problem of correlation method at the step of calculating gene bias level. Due to a failure of accurately defining the background mutation, the problem would result in wrong optimal codon identities. In other words, partial mutational effects on codon choices were mistakenly regarded as selective influences, leading to incorrect and biased optimal codon identities. Finally, considering the translational dynamics, optimal codons identified by comparison method can be well-explained by tRNA compositions, whereas optimal codons identified by correlation method can not be. For all above reasons, we conclude that real optimal codons actually do not track the genomic GC content, and correlation method is misleading in identifying optimal codons and better be avoided.     

Patterns of codon usage bias in three dicot and four monocot plant species

Codon usage in nuclear genes of four monocot and three dicot species was analyzed to find general patterns in codon choice of plant species. Codon bias was correlated with GC content at the third codon position. GC contents were higher in monocot species than in dicot species at all codon positions. The high GC contents of monocot species might be the result of relatively strong mutational bias that occurred in the lineage of the Poaceae species. In both dicot and monocot species, the effective number of codons (ENCs) for most genes was similar to that for the expected ENCs based on the GC content at the third codon positions. G and C ending codons were detected as the "preferred" codons in monocot species, as in Drosophila. Also, many "preferred" codons are the same in dicot species. Pyrimidine (C and T) is used more frequently than purine (G and A) in four-fold degenerate codon groups.     

Distinct paths to stop codon reassignment by the variant-code organisms Tetrahymena and Euplotes

The reassignment of stop codons is common among many ciliate species. For example, Tetrahymena species recognize only UGA as a stop codon, while Euplotes species recognize only UAA and UAG as stop codons. Recent studies have shown that domain 1 of the translation termination factor eRF1 mediates stop codon recognition. While it is commonly assumed that changes in domain 1 of ciliate eRF1s are responsible for altered stop codon recognition, this has never been demonstrated in vivo. To carry out such an analysis, we made hybrid proteins that contained eRF1 domain 1 from either Tetrahymena thermophila or Euplotes octocarinatus fused to eRF1 domains 2 and 3 from Saccharomyces cerevisiae. We found that the Tetrahymena hybrid eRF1 efficiently terminated at all three stop codons when expressed in yeast cells, indicating that domain 1 is not the sole determinant of stop codon recognition in Tetrahymena species. In contrast, the Euplotes hybrid facilitated efficient translation termination at UAA and UAG codons but not at the UGA codon. Together, these results indicate that while domain 1 facilitates stop codon recognition, other factors can influence this process. Our findings also indicate that these two ciliate species used distinct approaches to diverge from the universal genetic code.