Synonymous codon usage in Escherichia coli: selection for translational accuracy

In many organisms, selection acts on synonymous codons to improve translation. However, the precise basis of this selection remains unclear in the majority of species. Selection could be acting to maximize the speed of elongation, to minimize the costs of proofreading, or to maximize the accuracy of translation. Using several data sets, we find evidence that codon use in Escherichia coli is biased to reduce the costs of both missense and nonsense translational errors. Highly conserved sites and genes have higher codon bias than less conserved ones, and codon bias is positively correlated to gene length and production costs, both indicating selection against missense errors. Additionally, codon bias increases along the length of genes, indicating selection against nonsense errors. Doublet mutations or replacement substitutions do not explain our observations. The correlations remain when we control for expression level and for conflicting selection pressures at the start and end of genes. Considering each amino acid by itself confirms our results. We conclude that selection on synonymous codon use in E. coli is largely due to selection for translational accuracy, to reduce the costs of both missense and nonsense errors.     

Gradients in nucleotide and codon usage along Escherichia coli genes

The usage of codons and nucleotide combinations varies along genes and systematic variation causes gradients in usage. We have studied such gradients of nucleotides and nucleotide combinations and their immediate context in Escherichia coli. To distinguish mutational and selectional effects, the genes were subdivided into three groups with different codon usage bias and the gradients of nucleotide usage were studied in each group. Some combinations that can be associated with a propensity for processivity errors show strong negative gradients that become weaker in genes with low codon bias, consistent with a selection on translational efficiency. One of the strongest gradients is for third position G, which shows a pervasive positive gradient in usage in most contexts of surrounding bases.     

Optimality of codon usage in Escherichia coli due to load minimization

The canonical genetic code is known to be highly efficient in minimizing the effects of mistranslational errors and point mutations, an ability which in term is designated "load minimization". One parameter involved in calculating the load minimizing property of the genetic code is codon usage. In most bacteria, synonymous codons are not used with equal frequencies. Different factors have been proposed to contribute to codon usage preference. It has been shown that the codon preference is correlated with the composition of the tRNA pool. Selection for translational efficiency and translational accuracy both result in such a correlation. In this work, it is shown that codon usage bias in Escherichia coli works so as to minimize the consequences of translational errors, i.e. optimized for load minimization.     

High-level periplasmic expression in Escherichia coli using a eukaryotic signal peptide: importance of codon usage at the 5' end of the coding sequence

We investigated the ability of signal peptides of eukaryotic origin (human, mouse, and yeast) to efficiently direct model proteins to the Escherichia coli periplasm. These were compared against a well-characterized prokaryotic signal peptide-OmpA. Surprisingly, eukaryotic signal peptides can work very efficiently in E. coli, but require optimization of codon usage by codon-based mutagenesis of the signal peptide coding region. Analysis of the 5' of periplasmic and cytoplasmic E. coli genes shows some codon usage differences.     

Error minimization explains the codon usage of highly expressed genes in Escherichia coli

Different organisms use synonymous codons with different preferences. Several measures have been introduced to compute the extent of codon usage bias within a gene or genome, among which the codon adaptation index (CAI) has been shown to be well correlated with mRNA levels of Escherichia coli. In this work an error adaptation index (eAI) is introduced, which estimates the level at which a gene can tolerate the effects of mistranslations. It is shown that the eAI has a strong correlation with CAI, as well as with mRNA levels, which suggests that the codons of highly expressed genes are selected so that mistranslation would have the minimum possible effect on the structure and function of the related proteins.     

Horizontal gene transfer in glycosyl hydrolases inferred from codon usage in Escherichia coli and Bacillus subtilis.

Glycosyl hydrolase (GH) genes from Escherichia coli and Bacillus subtilis were used to search for cases of horizontal gene transfer. Such an event was inferred by G + C content, codon usage analysis, and a phylogenetic congruency test. The codon usage analysis used is a procedure based on a distance derived from a Pearson linear correlation coefficient determined from a pairwise codon usage comparison. The distances are then used to generate a distance-based tree with which we can define clusters and rapidly compare codon usage. Three genes (yagH from E. coli and xynA and xynB from B. subtilis) were determined to have arrived by horizontal gene transfer and were located in E. coli CP4-6 prophage, and B. subtilis prophages 6 and 5, respectively. In this study, we demonstrate that with codon usage analysis, the proposed horizontally transferred genes can be distinguished from highly expressed genes.     

The relationship between synonymous codon usage and protein structure in Escherichia coli and Homo sapiens

The role of silent position in the codon on the protein structure is an interesting and yet unclear problem. In this paper, 563 Homo sapiens genes and 417 Escherichia coli genes coding for proteins with four different folding types have been analyzed using variance analysis, a multivariate analysis method newly used in codon usage analysis, to find the correlation between amino acid composition, synonymous codon, and protein structure in different organisms. It has been found that in E. coli, both amino acid compositions in differently folded proteins and synonymous codon usage in different gene classes coding for differently folded proteins are significantly different. It was also found that only amino acid composition is different in different protein classes in H. sapiens. There is no universal correlation between synonymous codon usage and protein structure in these two different organisms. Further analysis has shown that GC content on the second codon position can distinguish coding genes for different folded proteins in both organisms.     

A comparative study of mutations in Escherichia coli and Salmonella typhimurium shows that codon conservation is strongly correlated with codon usage

Escherichia coli and Salmonella typhimurium are closely related species of enteric bacteria, having diverged from 120 to 160 million years ago, according to the estimate of Ochman & Wilson (1987. J. Mol. Evol.26, 74-86). In order to study base substitution mutations in the genomes of these bacteria, we have compared pairs of genes for the same product in the two species, and have selected a sample in which the protein length is the same in both E. coli and S. typhimurium. From the alignment of these gene pairs, we observe that frequently used codons are more conserved than infrequently used codons, i.e., the apparent mutation rate is higher for rare codons than for popular codons.     

Design, synthesis and expression of a human interleukin-2 gene incorporating the codon usage bias found in highly expressed Escherichia coli genes.

A synthetic gene encoding human interleukin-2 (IL-2) was designed such that the codon usage bias resembled that found in highly expressed Escherichia coli genes. The percentage of preferred codons was increased from 43% in the native cDNA sequence to 85% in the synthetic sequence. The cDNA and synthetic IL-2 genes were placed under the control of the trc promoter and expressed in E. coli JM101. While Northern blot analysis of IL-2 mRNA from each genetic construct demonstrated equivalent message half-lives, immunoblot and bioactivity analyses showed the synthetic gene to direct the synthesis of up to 16 times more IL-2 than the native cDNA sequence.     

High-level expression of Falcipain-2 in Escherichia coli by codon optimization and auto-induction

Falcipain-2, the major cysteine hemoglobinase from the human malaria parasite Plasmodium falciparum, is critical for parasite development and is considered a promising chemotherapeutic target. In order to facilitate the high-throughput screening of Falcipain-2 inhibitors from natural sources, we developed an economic and highly-productive overexpression system in Escherichia coli using a codon-optimized proFalcipain-2 construct. Very high expression levels (35-55% of total host proteins) were observed when proFalcaipain-2 expression was induced with 1mM isopropyl-1-thio-β-D-galactopyranoside (IPTG) in several E. coli strains, with the highest level observed for BL21(DE3). A lower expression (~40% of total host proteins) was observed when BL21(DE3) was grown in ZYM-5052 auto-induction medium, containing 0.2% lactose as inducer. However, the culture grew to notably higher cellular density, increasing ~1.5 times the overall yield of the system when compared with conventional IPTG-induction. Although several conditions were modified to achieve the expression of soluble and active Falcipain-2, the enzyme was mainly obtained in the form of insoluble aggregates. After purification and refolding, ~50 mg of active enzyme were obtained per liter of culture at low cost using a regular incubator shaker, and recombinant Falcipain-2 exhibited structural and functional characteristics very similar to the natural counterpart. Due to its versatility and simplicity, this strategy can be straightforwardly adapted to other proteins from Plasmodium species or any other organism with an AT-rich genome.     

Influence of the codon following the AUG initiation codon on the expression of a modified lacZ gene in Escherichia coli.

In a lacZ expression vector (pMC1403Plac), all 64 codons were introduced immediately 3' from the AUG initiation codon. The expression of the second codon variants was measured by immunoprecipitation of the plasmid-coded fusion proteins. A 15-fold difference in expression was found among the codon variants. No distinct correlation could be made with the level of tRNA corresponding to the codons and large differences were observed between synonymous codons that use the same tRNA. Therefore the effect of the second codon is likely to be due to the influence of its composing nucleotides, presumably on the structure of the ribosomal binding site. An analysis of the known sequences of a large number of Escherichia coli genes shows that the use of codons in the second position deviates strongly from the overall codon usage in E. coli. It is proposed that codon selection at the second position is not based on requirements of the gene product (a protein) but is determined by factors governing gene regulation at the initiation step of translation.     

Control of basal-level codon misreading in Escherichia coli

Basal-level misreading of asparagine codons was examined in a number of Escherichia coli strains. Lysine substitutions were measured by quantitating the amount of charge heterogeneity in MS2 coat protein. In most strains the heterogeneity was consistent with misreading of AAU codons at a frequency of 3-6 X 10(-3). Strains with streptomycin resistance mutations (rpsL) have reduced levels of misreading. There is no significant difference in the frequency of basal-level errors in stringent (relA+) and relaxed (relA) strains, even during starvation for amino acids unrelated to the substitution being studied.     

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.     

UGA suppression by normal tRNA Trp in Escherichia coli: codon context effects.

The nucleotide sequences at the 3' side of in-phase UGA termination codons in mRNAs of various prokaryotic genes were re-examined. An adenine (A) residue is found to be adjacent to the 3' side of UGA in mRNAs which code for readthrough proteins by the suppression of UGA by normal Escherichia coli tRNA Trp. It is suggested that the nature of the nucleotide following a UGA codon determines whether the UGA signals inefficiently or efficiently the termination of polypeptide chain synthesis: an A residue at this position permits the UGA readthrough process.