Correlation of codon biases and potential secondary structures with mRNA translation efficiency in unicellular organisms

Gene expression is known to correlate with the degree of codon bias in many unicellular organisms. However, such a correlation is not observed in some organisms. It was demonstrated that inverted complementary repeats within coding DNA sequences (ORFs) should be considered for proper estimation of the translation efficiency because they can form secondary structures that obstruct ribosome movement. A program was developed for estimating the potential expression of ORFs in unicellular organisms on the basis of their genome sequences. The program computes the elongation efficiency index (EEI) and takes into account three key factors: codon bias, the average number of inverted complementary repeats, and the free energies of potential stem-loop structures formed by these repeats. The influence of these factors on translation was numerically estimated. Their optimal ratio was computed for each organism. EEIs of 384 unicellular organisms (351 bacteria, 28 archaea, and 5 eukaryotes) were computed using the annotated genomes available from GenBank. Five potential evolutionary strategies of translational optimization were determined in the organisms studied. A considerable difference in preferential translational strategies was observed between bacteria and archaea. Significant correlations between EEIs and gene expression levels were shown for two species (Saccharomyces cerevisiae and Helicobacter pylori), using the available microarray data. The method allows the numerical estimation of the translation efficiency of an ORF and optimization of the nucleotide composition of heterologous genes in specified unicellular organisms. The program is available at http://wwwmgs.bionet.nsc.ru/mgs/programs/eei-calculator.     

Correlation of codon biases and potential secondary structures with mRNA translation efficiency in unicellular organisms

Gene expression is known to correlate with degree of codon bias in many unicellular organisms. However, such correlation is absent in some organisms. Recently we demonstrated that inverted complementary repeats within coding DNA sequence must be considered for proper estimation of translation efficiency, since they may form secondary structures that obstruct ribosome movement. We have developed a program for estimation of potential coding DNA sequence expression in defined unicellular organism using its genome sequence. The program computes elongation efficiency index. Computation is based on estimation of coding DNA sequence elongation efficiency, taking into account three key factors: codon bias, average number of inverted complementary repeats, and free energy of potential stem-loop structures formed by the repeats. The influence of these factors on translation is numerically estimated. An optimal proportion of these factors is computed for each organism individually. Quantitative translational characteristics of 384 unicellular organisms (351 bacteria, 28 archaea, 5 eukaryota) have been computed using their annotated genomes from NCBI GenBank. Five potential evolutionary strategies of translational optimization have been determined among studied organisms. A considerable difference of preferred translational strategies between Bacteria and Archaea has been revealed. Significant correlations between elongation efficiency index and gene expression levels have been shown for two organisms (S. cerevisiae and H. pylori) using available microarray data. The proposed method allows to estimate numerically the coding DNA sequence translation efficiency and to optimize nucleotide composition of heterologous genes in unicellular organisms. Availability: http://www.mgs.bionet.nsc.ru/mgs/programs/eei-calculator/.     

Nucleotide composition-based prediction of gene expression efficacy

A correlation between efficacy of gene expression and nucleotide composition of its protein-coding sequences was studied. It was found that measures based exclusively on codon frequencies are analogous to codon adaptation index [1] and do not adequately reveal this correlation. A more general measure is suggested, i.e., elongation efficacy index, which takes into account both the codon frequencies and the level of local mRNA complementarity. Utilization of this measure facilitated adequate recognition of highly expressed genes in 18 unicellular organisms, i.e., 14 eubacteria, three archaebacteria, and yeast.     

Dissecting the contributions of GC content and codon usage to gene expression in the model alga Chlamydomonas reinhardtii

The efficiency of gene expression in all organisms depends on the nucleotide composition of the coding region. GC content and codon usage are the two key sequence features known to influence gene expression, but the underlying molecular mechanisms are not entirely clear. Here we have determined the relative contributions of GC content and codon usage to the efficiency of nuclear gene expression in the unicellular green alga Chlamydomonas reinhardtii. By comparing gene variants that encode an identical amino acid sequence but differ in their GC content and/or codon usage, we show that codon usage is the key factor determining translational efficiency and, surprisingly, also mRNA stability. By contrast, unfavorable GC content affects gene expression at the level of the chromatin structure by triggering heterochromatinization. We further show that mutant algal strains that permit high‐level transgene expression are less susceptible to epigenetic transgene suppression and do not establish a repressive chromatin structure at the transgenic locus. Our data disentangle the relationship between GC content and codon usage, and suggest simple strategies to overcome the transgene expression problem in Chlamydomonas.     

Evolutionary rates and expression level in Chlamydomonas

In many biological systems, especially bacteria and unicellular eukaryotes, rates of synonymous and nonsynonymous nucleotide divergence are negatively correlated with the level of gene expression, a phenomenon that has been attributed to natural selection. Surprisingly, this relationship has not been examined in many important groups, including the unicellular model organism Chlamydomonas reinhardtii. Prior to this study, comparative data on protein-coding sequences from C. reinhardtii and its close noninterfertile relative C. incerta were very limited. We compiled and analyzed protein-coding sequences for 67 nuclear genes from these taxa; the sequences were mostly obtained from the C. reinhardtii EST database and our C. incerta EST data. Compositional and synonymous codon usage biases varied among genes within each species but were highly correlated between the orthologous genes of the two species. Relative rates of synonymous and nonsynonymous substitution across genes varied widely and showed a strong negative correlation with the level of gene expression estimated by the codon adaptation index. Our comparative analysis of substitution rates in introns of lowly and highly expressed genes suggests that natural selection has a larger contribution than mutation to the observed correlation between evolutionary rates and gene expression level in Chlamydomonas.     

Translation efficiency in humans: tissue specificity,global optimization and differences between developmental stages

Various studies in unicellular and multicellular organisms have shown that codon bias plays a significant role in translation efficiency (TE) by co-adaptation to the tRNA pool. Yet, in humans and other mammals the role of codon bias is still an open question, with contradictory results from different studies. Here we address this question, performing a large-scale tissue-specific analysis of TE in humans, using the tRNA Adaptation Index (tAI) as a direct measure for TE. We find tAI to significantly correlate with expression levels both in tissue-specific and in global expression measures, testifying to the TE of human tissues. Interestingly, we find significantly higher correlations in adult tissues as opposed to fetal tissues, suggesting that the tRNA pool is more adjusted to the adult period. Optimization based analysis suggests that the tRNA pool—codon bias co-adaptation is globally (and not tissue-specific) driven. Additionally, we find that tAI correlates with several measures related to the protein functionally importance, including gene essentiality. Using inferred tissue-specific tRNA pools lead to similar results and shows that tissue-specific genes are more adapted to their tRNA pool than other genes and that related sets of functional gene groups are translated efficiently in each tissue. Similar results are obtained for other mammals. Taken together, these results demonstrate the role of codon bias in TE in humans, and pave the way for future studies of tissue-specific TE in multicellular organisms.     

Codon usage-mediated inhibition of HIV-1 gag expression in mammalian cells occurs independently of translation

Codon usage is considered one of the critical factors that limit the expression rate of heterologous genes. Impaired translation efficiency, specifically insufficient amount of corresponding tRNAs and changed startcodon context, are believed to account for the low translation initiation and elongation rates during the protein biosynthesis in unicellular organisms. Translational efficiency is probably not the primary factor influencing codon usage diversity in mammalian cells. However, the other possible mechanisms preventing expression of genes with low-usage such as mRNA stability, processing and nucleocytoplasmic transport, are not adequately explored. In our work, we addressed the question of whether codon usage differences affect exclusively translational efficiency of mammalian gene products. We demonstrated that the CMV-induced expression of gag-reporter in human H1299 cell line was influenced by the nucleotide composition of the mRNA, and the limitation of gag expression appeared to be inversely related to the level of codon optimization. However, cytoplasmic expression of the gag-reporter driven by vaccinia virus/T7 RNA polymerase hybrid system rescued its expression independently of HIV-1 gag mRNA nucleotide content. We concluded that impaired HIV-1 gag expression may be caused by translation-independent mechanisms, which probably play a major role in codon usage-mediated defects in heterologous gene expression in mammalian cells.     

Intron Length and Codon Usage

The correlation was shown between the length of introns and the codon usage of the coding sequences of the corresponding genes, which in some cases can be related to the level of gene expression. The link is positive in the unicellular organisms, i.e., genes with the longer introns show the higher bias of codon usage. It is most pronounced in baker's yeast, where it is definitely related to the level of gene expression—genes with the higher level of expression have the longer introns. The correlation is inverted in multicellular organisms as compared to unicellular ones. Some organisms, however, do not show the link. The presence or absence of the link does not seem to be related to the GC percent of the coding sequences. Received: 7 December 1999 / Accepted: 10 May 2000     

A role for tRNA modifications in genome structure and codon usage

Transfer RNA (tRNA) gene content is a differentiating feature of genomes that contributes to the efficiency of the translational apparatus, but the principles shaping tRNA gene copy number and codon composition are poorly understood. Here, we report that the emergence of two specific tRNA modifications shaped the structure and composition of all extant genomes. Through the analysis of more than 500 genomes, we identify two kingdom-specific tRNA modifications as major contributors that separated archaeal, bacterial, and eukaryal genomes in terms of their tRNA gene composition. We show that, contrary to prior observations, genomic codon usage and tRNA gene frequencies correlate in all kingdoms if these two modifications are taken into account and that presence or absence of these modifications explains patterns of gene expression observed in previous studies. Finally, we experimentally demonstrate that human gene expression levels correlate well with genomic codon composition if these identified modifications are considered.     

Efficient manipulations of synonymous mutations for controlling translation rate: an analytical approach

Gene translation is a central process in all living organism with important ramifications to almost every biomedical field. Previous systems evolutionary studies in the field have demonstrated that in many organisms coding sequence features undergo selection to optimize this process. In the current study, we report for the first time analytical proofs related to the various aspects of this process and its optimality. Among our results we show that coding sequences with mono- tonic increasing profiles of translation efficiency (i.e., with slower codons near the 5'UTR), mathematically optimize ribosomal allocation by minimizing the number of ribosomes needed for translating a codon per time unit. Thus, the genomic translation efficiency profile reported in previous studies for many organisms is optimal in this sense. In addition, we show that improving translation efficiency of a codon in a gene may result in a decrease in the translation rate of other genes, demonstrating that the relation between codon bias and protein translation rate is less trivial than was assumed before. Based on these observations we describe an efficient heuristic for designing coding sequences with specific translation efficiency and minimal ribosomal allocation for heterologous gene expression. We demonstrate how this heuristic can be used in biotechnology for engineering a heterologous gene before expressing it in a new host.     

A Universal Trend of Reduced mRNA Stability near the Translation-Initiation Site in Prokaryotes and Eukaryotes

Recent studies have suggested that the thermodynamic stability of mRNA secondary structure near the start codon can regulate translation efficiency in Escherichia coli, and that translation is more efficient the less stable the secondary structure. We survey the complete genomes of 340 species for signals of reduced mRNA secondary structure near the start codon. Our analysis includes bacteria, archaea, fungi, plants, insects, fishes, birds, and mammals. We find that nearly all species show evidence for reduced mRNA stability near the start codon. The reduction in stability generally increases with increasing genomic GC content. In prokaryotes, the reduction also increases with decreasing optimal growth temperature. Within genomes, there is variation in the stability among genes, and this variation correlates with gene GC content, codon bias, and gene expression level. For birds and mammals, however, we do not find a genome-wide trend of reduced mRNA stability near the start codon. Yet the most GC rich genes in these organisms do show such a signal. We conclude that reduced stability of the mRNA secondary structure near the start codon is a universal feature of all cellular life. We suggest that the origin of this reduction is selection for efficient recognition of the start codon by initiator-tRNA.     

Translational selection shapes codon usage in the GC-rich genome of Chlamydomonas reinhardtii

In unicellular species codon usage is determined by mutational biases and natural selection. Among prokaryotes, the influence of these factors is different if the genome is skewed towards AT or GC, since in AT-rich organisms translational selection is absent. On the other hand, in AT-rich unicellular eukaryotes the two factors are present. In order to understand if GC-rich genomes display a similar behavior, the case of Chlamydomonas reinhardtii was studied. Since we found that translational selection strongly influences codon usage in this species, we conclude that there is not a common pattern among unicellular organisms.     

Conservation of the relative tRNA composition in healthy and cancerous tissues

Elongation in protein translation is strongly dependent on the availability of mature transfer RNAs (tRNAs). The relative concentrations of the tRNA isoacceptors determine the translation efficiency in unicellular organisms. However, the degree of correspondence of codons and the relevant tRNA isoacceptors serves as an estimator for translation efficiency in all organisms. In this study, we focus on the translational capacity of the human proteome. We show that the correspondence between the codon usage and tRNAs can be improved by combining experimental measurements with the genomic copy number of isoacceptor groups. We show that there are technologies of tRNA measurements that are useful for our analysis. However, fragments of tRNAs do not agree with translational capacity. It was shown that there is a significant increase in the absolute levels of tRNA genes in cancerous cells in comparison to healthy cells. However, we find that the relative composition of tRNA isoacceptors in healthy, cancerous, or transformed cells remains almost identical. This result may indicate that maintaining the relative tRNA composition in cancerous cells is advantageous via its stabilizing of the effectiveness of translation.     

Predicting gene expression level from codon usage bias

The "expression measure" of a gene, E(g), is a statistic devised to predict the level of gene expression from codon usage bias. E(g) has been used extensively to analyze prokaryotic genome sequences. We discuss 2 problems with this approach. First, the formulation of E(g) is such that genes with the strongest selected codon usage bias are not likely to have the highest predicted expression levels; indeed the correlation between E(g) and expression level is weak among moderate to highly expressed genes. Second, in some species, highly expressed genes do not have unusual codon usage, and so codon usage cannot be used to predict expression levels. We outline a simple approach, first to check whether a genome shows evidence of selected codon usage bias and then to assess the strength of bias in genes as a guide to their likely expression level; we illustrate this with an analysis of Shewanella oneidensis.     

Introduction of amino acid residues at the N-terminus of the zeocin-resistance protein increases its expression in Saccharomyces cerevisiae

Nucleotide sequences proximal to the initiation codon of a gene are known to affect the expression efficiency of that gene. We screened 10-bp random sequences upstream of the initiation codon of the zeocin-resistance gene to identify sequences that could enhance its expression in Saccharomyces cerevisiae. Of the isolated sequences, 20 sequences exhibited a common feature, i.e. ATG at the position −9 through −7, which resulted in the incorporation of three amino acids at the N-terminus of the protein. The introduction of these sequences upstream of the initiation codon increased the expression levels of zeocin-resistance protein by 2.2–6.5-fold. One of these sequences increased the expression levels of three out of four human proteins, thereby suggesting that this sequence may also enhance the expression efficiency of mammalian proteins in yeast.     

Codon usage and amino acid usage influence genes expression level

Highly expressed genes in any species differ in the usage frequency of synonymous codons. The relative recurrence of an event of the favored codon pair (amino acid pairs) varies between gene and genomes due to varying gene expression and different base composition. Here we propose a new measure for predicting the gene expression level, i.e., codon plus amino bias index (CABI). Our approach is based on the relative bias of the favored codon pair inclination among the genes, illustrated by analyzing the CABI score of the Medicago truncatula genes. CABI showed strong correlation with all other widely used measures (CAI, RCBS, SCUO) for gene expression analysis. Surprisingly, CABI outperforms all other measures by showing better correlation with the wet-lab data. This emphasizes the importance of the neighboring codons of the favored codon in a synonymous group while estimating the expression level of a gene.     

Swapped green algal promoters: aphVIII-based gene constructs with Chlamydomonas flanking sequences work as dominant selectable markers in Volvox and vice versa

Production of transgenic organisms is a well-established, versatile course of action in molecular biology. Genetic engineering often requires heterologous, dominant antibiotic resistance genes that have been used as selectable markers in many species. However, as heterologous 5′ and 3′ flanking sequences often result in very low expression rates, endogenous flanking sequences, especially promoters, are mostly required and are easily obtained in model organisms, but it is much more complicated and time-consuming to get appropriate sequences from less common organisms. In this paper, we show that aminoglycoside 3′-phosphotransferase gene (aphVIII) based constructs with 3′ and 5′ untranslated flanking sequences (including promoters) from the multicellular green alga Volvox work in the unicellular green alga Chlamydomonas and flanking sequences from Chlamydomonas work in Volvox, at least if a low expression rate is compensated by an enforced high gene dosage. This strategy might be useful for all investigators that intend to transform species in which genomic sequences are not available, but sequences from related organisms exist.     

Universal function-specificity of codon usage

Synonymous codon usage has long been known as a factor that affects average expression level of proteins in fast-growing microorganisms, but neither its role in dynamic changes of expression in response to environmental changes nor selective factors shaping it in the genomes of higher eukaryotes have been fully understood. Here, we propose that codon usage is ubiquitously selected to synchronize the translation efficiency with the dynamic alteration of protein expression in response to environmental and physiological changes. Our analysis reveals that codon usage is universally correlated with gene function, suggesting its potential contribution to synchronized regulation of genes with similar functions. We directly show that coexpressed genes have similar synonymous codon usages within the genomes of human, yeast, Caenorhabditis elegans and Escherichia coli. We also demonstrate that perturbing the codon usage directly affects the level or even direction of changes in protein expression in response to environmental stimuli. Perturbing tRNA composition also has tangible phenotypic effects on the cell. By showing that codon usage is universally function-specific, our results expand, to almost all organisms, the notion that cells may need to dynamically alter their intracellular tRNA composition in order to adapt to their new environment or physiological role.     

Mutation bias is the driving force of codon usage in the Gallus gallus genome

Synonymous codons are used with different frequencies both among species and among genes within the same genome and are controlled by neutral processes (such as mutation and drift) as well as by selection. Up to now, a systematic examination of the codon usage for the chicken genome has not been performed. Here, we carried out a whole genome analysis of the chicken genome by the use of the relative synonymous codon usage (RSCU) method and identified 11 putative optimal codons, all of them ending with uracil (U), which is significantly departing from the pattern observed in other eukaryotes. Optimal codons in the chicken genome are most likely the ones corresponding to highly expressed transfer RNA (tRNAs) or tRNA gene copy numbers in the cell. Codon bias, measured as the frequency of optimal codons (Fop), is negatively correlated with the G + C content, recombination rate, but positively correlated with gene expression, protein length, gene length and intron length. The positive correlation between codon bias and protein, gene and intron length is quite different from other multi-cellular organism, as this trend has been only found in unicellular organisms. Our data displayed that regional G + C content explains a large proportion of the variance of codon bias in chicken. Stepwise selection model analyses indicate that G + C content of coding sequence is the most important factor for codon bias. It appears that variation in the G + C content of CDSs accounts for over 60% of the variation of codon bias. This study suggests that both mutation bias and selection contribute to codon bias. However, mutation bias is the driving force of the codon usage in the Gallus gallus genome. Our data also provide evidence that the negative correlation between codon bias and recombination rates in G. gallus is determined mostly by recombination-dependent mutational patterns.