Diversity of preferred nucleotide sequences around the translation initiation codon in eukaryote genomes

Understanding regulatory mechanisms of protein synthesis in eukaryotes is essential for the accurate annotation of genome sequences. Kozak reported that the nucleotide sequence GCCGCC(A/G)CCAUGG (AUG is the initiation codon) was frequently observed in vertebrate genes and that this 'consensus' sequence enhanced translation initiation. However, later studies using invertebrate, fungal and plant genes reported different 'consensus' sequences. In this study, we conducted extensive comparative analyses of nucleotide sequences around the initiation codon by using genomic data from 47 eukaryote species including animals, fungi, plants and protists. The analyses revealed that preferred nucleotide sequences are quite diverse among different species, but differences between patterns of nucleotide bias roughly reflect the evolutionary relationships of the species. We also found strong biases of A/G at position -3, A/C at position -2 and C at position +5 that were commonly observed in all species examined. Genes with higher expression levels showed stronger signals, suggesting that these nucleotides are responsible for the regulation of translation initiation. The diversity of preferred nucleotide sequences around the initiation codon might be explained by differences in relative contributions from two distinct patterns, GCCGCCAUG and AAAAAAAUG, which implies the presence of multiple molecular mechanisms for controlling translation initiation.     

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.     

Context sequences of translation initiation codon in plants

In this survey of 5074 plant genes for their AUG context sequences, purines are present at the _3 and +4 positions in about 80% of the sequences. Although this observation is similar to the vertebrate consensus sequence, the number of plant mRNAs with purines at the _3 position is lower and at the +4 position is higher than reported for vertebrate mRNAs. Higher plants have an AC-rich consensus sequence, caA(A/C)aAUGGCg as a context of translation initiator codon. Between the two major groups of angiosperms, the context of the AUG codon in dicot mRNAs is aaA(A/C)aAUGGCu which is similar to the higher-plant consensus but monocot mRNAs have c(a/c)(A/G)(A/C)cAUGGCG as a consensus which exhibits an overall similarity with the vertebrate consensus. The experimental evidence regarding the importance of the AUG context in plants is discussed.     

Translation initiation AUG context varies with codon usage bias and gene length in Drosophila melanogaster

The relationship between the codon usage bias and the sequence context surrounding the AUG translation initiation codon was examined in 1100 Drosophila melanogaster mRNA sequences. The codon usage bias measured by the "codon adaptation index" (CAI), and the effectiveness of the AUG context for translation initiation assessed by the "AUG context adaptation index" (AUGCAI), showed a significant positive relationship (correlation coefficient: r = 0.34, p <0.0001), indicating that these two factors are evolutionally under a similar natural selection constraint at the translational level. The importance of each position of the AUG context in relation to codon usage bias was examined, and the preference for the nucleotide at the -13, -12, -11, -10, -7, -6, -5, -4, -3, -2, and -1 positions showed a significant positive correlation to the codon usage bias, suggesting the action of natural selection on these very specific positions of the Drosophila genome. The relationship between AUGCAI value and gene length was also examined, and a significant negative relationship was found (r = -0.15, p <0.0001), suggesting a general tendency of higher expressivity of shorter genes, and of lower expressivity of longer genes in D. melanogaster.     

鸡基因组差异表达基因翻译起始密码子上下游序列的比较研究

翻译起始调控是基因表达调控的一个关键步骤之一。本文以鸡为研究材料,比较研究了鸡基因组高表达基因和低表达基因翻译起始密码子上下游的碱基序列差异,旨在寻找影响鸡基因表达水平的特异性调控位点。全部3 020个单剪接基因完整的mRNA序列及有详细注释的5＇UTRs序列从Ensembl下载。编写计算机程序,读取每个基因mRNA起始密码子上下游各位点的碱基。研究发现,起始密码子上游-3、-2位点可能是鸡基因组基因表达起始密码子正确识别的关键位点。起始密码子上下游的碱基组成分析发现,高表达基因和低表达基因起始密码子的上游均倾向使用（G＋C）,高表达基因的使用偏倚尤为强烈。序列差异比较发现,高表达基因在-9、-6、-3、＋4位点显著偏向G,在-1、-2、-4、-5位点显著偏向C。低表达基因起始密码子上游使用A、U的频率显著高于低表达基因。在-19位点强烈偏向A,在＋1、＋11、＋14位点强烈偏向U。     

Cooperative effects by the initiation codon and its flanking regions on translation initiation

The purine-rich Shine-Dalgarno (SD) sequence located a few bases upstream of the mRNA initiation codon supports translation initiation by complementary binding to the anti-SD in the 16S rRNA, close to its 3' end. AUG is the canonical initiation codon but the weaker UUG and GUG codons are also used for a minority of genes. The codon sequence of the downstream region (DR), including the +2 codon immediately following the initiation codon, is also important for initiation efficiency. We have studied the interplay between these three initiation determinants on gene expression in growing Escherichia coli. One optimal SD sequence (SD(+)) and one lacking any apparent complementarity to the anti-SD in 16S rRNA (SD(-)) were analyzed. The SD(+) and DR sequences affected initiation in a synergistic manner and large differences in the effects were found. The gene expression level associated with the most efficient of these DRs together with SD(-) was comparable to that of other DRs together with SD(+). The otherwise weak initiation codon UUG, but not GUG, was comparable with AUG in strength, if placed in the context of two of the DRs. The +2 codon was one, but not the only, determinant for this unexpectedly high efficiency of UUG.     

Alpha-thalassemia due to the deletion of nucleotides -2 and -3 preceding the AUG initiation codon affects translation efficiency both in vitro and in vivo.

We previously hypothesized that a 2 nucleotide deletion, causing a A-greater than C change at position -3 preceding the ATG initiation codon of alpha globin gene, reduced translation efficiency of alpha globin mRNA and was responsible for a form of alpha + thalassemia displayed by an Algerian patient. We presently show that this deletion leads to a 30-45% reduction in translation efficiency of synthetic alpha globin mRNA in rabbit reticulocyte lysate. In other experiments, we constructed alpha/G gamma hybrid globin genes in which the 3' end of normal or mutated alpha globin genes downstream to the ATG initiation codon was substituted by the 3' part of a G gamma globin gene. COS cells transfected with either of these 2 hybrid genes were shown to synthesize a similar amount of alpha/G gamma hybrid mRNAs but 50% less G gamma globin when transfected with the alpha/G gamma hybrid gene carrying the deletion. These results definitively establish that the 2 nucleotide deletion reduces translation efficiency by 30-50%. This contrasts with the 93% reduction induced by a similar A-greater than C change at position -3 in the different nucleotide context preceding the ATG codon of the rat preproinsulin gene.     

The +4G site in Kozak consensus is not related to the efficiency of translation initiation

The optimal context for translation initiation in mammalian species is GCCRCCaugG (where R = purine and "aug" is the initiation codon), with the -3R and +4G being particularly important. The presence of +4G has been interpreted as necessary for efficient translation initiation. Accumulated experimental and bioinformatic evidence has suggested an alternative explanation based on amino acid constraint on the second codon, i.e., amino acid Ala or Gly are needed as the second amino acid in the nascent peptide for the cleavage of the initiator Met, and the consequent overuse of Ala and Gly codons (GCN and GGN) leads to the +4G consensus. I performed a critical test of these alternative hypotheses on +4G based on 34169 human protein-coding genes and published gene expression data. The result shows that the prevalence of +4G is not related to translation initiation. Among the five G-starting codons, only alanine codons (GCN), and glycine codons (GGN) to a much smaller extent, are overrepresented at the second codon, whereas the other three codons are not overrepresented. While highly expressed genes have more +4G than lowly expressed genes, the difference is caused by GCN and GGN codons at the second codon. These results are inconsistent with +4G being needed for efficient translation initiation, but consistent with the proposal of amino acid constraint hypothesis.     

A role for initiation codon context in chloroplast translation

To study the role of initiation codon context in chloroplast protein synthesis, we mutated the three nucleotides immediately upstream of the initiation codon (the -1 triplet) of two chloroplast genes in the alga Chlamydomonas reinhardtii. In prokaryotes, the -1 triplet has been proposed to base pair with either the 530 loop of 16S rRNA or the extended anticodon of fMet-tRNA. We found that in vivo, none of the chloroplast mutations affected mRNA stability. However, certain mutations did cause a temperature-sensitive decrease in translation and a more dramatic decrease at room temperature when combined with an AUU initiation codon. These mutations disrupt the proposed extended base pairing interaction with the fMet-tRNA anticodon loop, suggesting that this interaction may be important in vivo. Mutations that would still permit base pairing with the 530 loop of the 16S rRNA also had a negative effect on translation, suggesting that this interaction does not occur in vivo. Extended base pairing surrounding the initiation codon may be part of a mechanism to compensate for the lack of a classic Shine-Dalgarno rRNA interaction in the translation of some chloroplast mRNAs.     

Enhancement of translation initiation by A/T-rich sequences downstream of the initiation codon in Escherichia coli

The region located downstream of the initiation codon constitutes part of the translation initiation signal, significantly affecting the level of protein expression in E. coli. In order to determine its influence on translation initiation, we inserted random 12-base sequences downstream of the initiation codon of the lacZ gene. A total of 119 random clones showing higher beta-galactosidase activities than the control lacZ gene were isolated and subsequently sequenced. Analysis of these clones revealed that their insertion sequences are strikingly rich in A and T, but poor in G, with no consensus sequences among them. Toeprinting experiments and polysome profile analysis confirmed that the A/T-rich sequences enhance translation at the level of initiation. Collectively, the present data demonstrate that A/T richness of the region following the initiation codon plays a significant role in E. coli gene expression.     

Codon bias at the 3'-side of the initiation codon is correlated with translation initiation efficiency in Escherichia coli

The codon that follows the AUG initiation triplet (+2 codon) affects gene expression in Escherichia coli. We have extended this analysis using two model genes lacking any apparent Shine-Dalgarno sequence. Depending on the identity of the +2 codon a difference in gene expression up to 20-fold could be obtained. The effects did not correlate with the levels of intracellular pools of cognate tRNA for the +2 codon, with putative secondary mRNA structures, or with mRNA stability. However, most +2 iso-codons that were decoded by the same species of tRNA gave pairwise similar effects, suggesting that the effect on gene expression was associated with the decoding tRNA. High adenine content of the +2 codon was associated with high gene expression. Of the fourteen +2 codons that mediated the highest efficiency, all except two had an adenine as the first base of the codon. Analysis of the 3540 E. coli genes from the TransTerm database revealed that codons associated with high gene expression in the two expression systems are over-represented at the +2 position in natural genes. Codons that are associated with low gene expression are under-represented. The data suggest that evolution has favored codons at the +2 position that give high translation initiation.     

Correlation between sequence conservation of the 5' untranslated region and codon usage bias in Mus musculus genes

The codon adaptation index (CAI) values of all protein-coding sequences of the full-length cDNA libraries of Mus musculus were computed based on the RIKEN mouse full-length cDNA library. We have also computed the extent of consensus in flanking sequences of the initiator ATG codon based on the 'relative entropy' values of respective nucleotide positions (from -20 to +12 bp relative to the initiator ATG codon) for each group of genes classified by CAI values. With regard to the two nucleotides positions (-3 and +4) known to be highly conserved in Kozak's consensus sequence, a clear correlation between CAI values and relative entropy values was observed at position -3 but this was not significant at position +4, although a significant correlation was found at position -1 of the consensus sequence. Further, although no correlation was observed at any additional positions, relative entropy values were very high at positions -4, -6, and -8 in genes with high CAI values. These findings suggest that the extent of conservation in the flanking sequence of the initiator ATG codon including Kozak's consensus sequence was an important factor in modulation of the translation efficiency as well as synonymous codon usage bias particularly in highly expressed genes.     

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.     

Stringency of start codon selection modulates autoregulation of translation initiation factor eIF5

An AUG in an optimal nucleotide context is the preferred translation initiation site in eukaryotic cells. Interactions among translation initiation factors, including eIF1 and eIF5, govern start codon selection. Experiments described here showed that high intracellular eIF5 levels reduced the stringency of start codon selection in human cells. In contrast, high intracellular eIF1 levels increased stringency. High levels of eIF5 induced translation of inhibitory upstream open reading frames (uORFs) in eIF5 mRNA that initiate with AUG codons in conserved poor contexts. This resulted in reduced translation from the downstream eIF5 start codon, indicating that eIF5 autoregulates its own synthesis. As with eIF1, which is also autoregulated through translation initiation, features contributing to eIF5 autoregulation show deep evolutionary conservation. The results obtained provide the basis for a model in which auto- and cross-regulation of eIF5 and eIF1 translation establish a regulatory feedback loop that would stabilize the stringency of start codon selection.     

Effect of codon adaptation on codon-level and gene-level translation efficiency in vivo

Background

There is a significant difference between synonymous codon usage in many organisms, and it is known that codons used more frequently generally showed efficient decoding rate. At the gene level, however, there are conflicting reports on the existence of a correlation between codon adaptation and translation efficiency, even in the same organism.

Results

To resolve this issue, we cultured Escherichia coli under conditions designed to maintain constant levels of mRNA and protein and subjected the cells to ribosome profiling (RP) and mRNA-seq analyses. We showed that the RP results correlated more closely with protein levels generated under similar culture conditions than with the mRNA abundance from the mRNA-seq. Our result indicated that RP/mRNA ratio could be used as a measure of translation efficiency at gene level. On the other hand, the RP data showed that codon-specific ribosome density at the decoding site negatively correlated with codon usage, consistent with the hypothesis that preferred codons display lower ribosome densities due to their faster decoding rate. However, highly codon-adapted genes showed higher ribosome densities at the gene level, indicating that the efficiency of translation initiation, rather than higher elongation efficiency of preferred codons, exerted a greater effect on ribosome density and thus translation efficiency.

Conclusions

These findings indicate that evolutionary pressure on highly expressed genes influenced both codon bias and translation initiation efficiency and therefore explains contradictory findings that codon usage bias correlates with translation efficiency of native genes, but not with the artificially created gene pool, which was not subjected to evolution pressure.

Electronic supplementary material

The online version of this article (doi:10.1186/1471-2164-15-1115) contains supplementary material, which is available to authorized users.     

Translation of chloroplast-encoded mRNA: potential initiation and termination signals.

A survey of 196 protein-coding chloroplast DNA sequences demonstrated the preference for AUG and UAA codons for initiation and termination of translation, respectively. As in prokaryotes at every nucleotide position from -25 to +25 (AUG is +1 to +3) and for 25 nucleotides 5' and 3' to the termination codon an A or U is predominant, except for C at +5 and G at +22. A Shine-Dalgarno (SD) sequence (GGAGG or tri- or tetranucleotide variant) was found within 100 bp 5' to the AUG codon in 92% of the genes. In 40% of these cases, the location of the SD sequence was similar to that of the consensus for prokaryotes (-12 to -7 5' to AUG), presumed to be optimal for translation initiation. A SD sequence could not be located in 6% of the chloroplast sequences. We propose that mRNA secondary structures may be required for the relocation of a distal SD sequences to within the optimal region (-12 to -7) for initiation of translation. We further suggest that termination at UGA codons in chloroplast genes may occur by a mechanism, involving 16S rRNA secondary structure, which has been proposed for UGA termination in E. coli.     

The initiation codon determines the efficiency but not the site of translation initiation in Chlamydomonas chloroplasts.

To study translation initiation in Chlamydomonas chloroplasts, we mutated the initiation codon AUG to AUU, ACG, ACC, ACU, and UUC in the chloroplast petA gene, which encodes cytochrome f of the cytochrome b6/f complex. Cytochrome f accumulated to detectable levels in all mutant strains except the one with a UUC codon, but only the mutant with an AUU codon grew well at 24 degrees C under conditions that require photosynthesis. Because no cytochrome f was detectable in the UUC mutant and because each mutant that accumulated cytochrome f did so at a different level, we concluded that any residual translation probably initiates at the mutant codon. As a further demonstration that alternative initiation sites are not used in vivo, we introduced in-frame UAA stop codons immediately downstream or upstream or in place of the initiation codon. Stop codons at or downstream of the initiation codon prevented accumulation of cytochrome f, whereas the one immediately upstream of the initiation codon had no effect on the accumulation of cytochrome f. These results suggest that an AUG codon is not required to specify the site of translation initiation in chloroplasts but that the efficiency of translation initiation depends on the identity of the initiation codon.     

Codon usage and protein length-dependent feedback from translation elongation regulates translation initiation and elongation speed

Essential cellular functions require efficient production of many large proteins but synthesis of large proteins encounters many obstacles in cells. Translational control is mostly known to be regulated at the initiation step. Whether translation elongation process can feedback to regulate initiation efficiency is unclear. Codon usage bias, a universal feature of all genomes, plays an important role in determining gene expression levels. Here, we discovered that there is a conserved but codon usage-dependent genome-wide negative correlation between protein abundance and CDS length. The codon usage effects on protein expression and ribosome flux on mRNAs are influenced by CDS length; optimal codon usage preferentially promotes production of large proteins. Translation of mRNAs with long CDS and non-optimal codon usage preferentially induces phosphorylation of initiation factor eIF2α, which inhibits translation initiation efficiency. Deletion of the eIF2α kinase CPC-3 (GCN2 homolog) in Neurospora preferentially up-regulates large proteins encoded by non-optimal codons. Surprisingly, CPC-3 also inhibits translation elongation rate in a codon usage and CDS length-dependent manner, resulting in slow elongation rates for long CDS mRNAs. Together, these results revealed a codon usage and CDS length-dependent feedback mechanism from translation elongation to regulate both translation initiation and elongation kinetics.