Comparative analysis of the base biases at the gene terminal portions in seven eukaryote genomes

Adenine nucleotides have been found to appear preferentially in the regions after the initiation codons or before the termination codons of bacterial genes. Our previous experiments showed that AAA and AAT, the two most frequent second codons in Escherichia coli, significantly enhance translation efficiency. To determine whether such a characteristic feature of base frequencies exists in eukaryote genes, we performed a comparative analysis of the base biases at the gene terminal portions using the proteomes of seven eukaryotes. Here we show that the base appearance at the codon third positions of gene terminal regions is highly biased in eukaryote genomes, although the codon third positions are almost free from amino acid preference. The bias changes depending on its position in a gene, and is characteristic of each species. We also found that bias is most outstanding at the second codon, the codon after the initiation codon. NCN is preferred in every genome; in particular, GCG is strongly favored in human and plant genes. The presence of the bias implies that the base sequences at the second codon affect translation efficiency in eukaryotes as well as bacteria.     

An ab initio molecular orbital study on the sequence-dependency of DNA conformation: an evaluation of intra- and inter-strand stacking interaction energy

The various nearest neighbor stacking interaction energies of stacked base pairs in the DNA double helix are calculated for both A- and B-type conformations using an ab initio molecular orbital method. It is demonstrated that the sequence-dependent conformational preference for A- or B-type results from the stacking interaction. In particular, the base sequence showing the highest preference for an A-type conformation is revealed as GC/GC, and the one with the next highest preference, AT/AT; for a B-type conformation, the respective sequences are CG/CG and CA/TG. The overall conformation of a DNA fragment is not determined by these particular sequences only but is influenced by all base pair steps. An intrinsically favorable conformation is predicted from the constituent stacking interaction.     

tRNA properties help shape codon pair preferences in open reading frames

Translation elongation is an accurate and rapid process, dependent upon efficient juxtaposition of tRNAs in the ribosomal A- and P-sites. Here, we sought evidence of A- and P-site tRNA interaction by examining bias in codon pair choice within open reading frames from a range of genomes. Three distinct and marked effects were revealed once codon and dipeptide biases had been subtracted. First, in the majority of genomes, codon pair preference is primarily determined by a tetranucleotide combination of the third nucleotide of the P-site codon, and all 3 nt of the A-site codon. Second, pairs of rare codons are generally under-used in eukaryotes, but over-used in prokaryotes. Third, the analysis revealed a highly significant effect of tRNA-mediated selection on codon pairing in unicellular eukaryotes, Bacillus subtilis, and the gamma proteobacteria. This was evident because in these organisms, synonymous codons decoded in the A-site by the same tRNA exhibit significantly similar P-site pairing preferences. Codon pair preference is thus influenced by the identity of A-site tRNAs, in combination with the P-site codon third nucleotide. Multivariate analysis identified conserved nucleotide positions within A-site tRNA sequences that modulate codon pair preferences. Structural features that regulate tRNA geometry within the ribosome may govern genomic codon pair patterns, driving enhanced translational fidelity and/or rate.     

Clonorchis sinensis: codon usage in nuclear genes

Codon usage in Clonorchis sinensis was analyzed using 12,515 codons from 38 coding sequences. Total GC content was 49.83%, and GC1, GC2 and GC3 contents were 56.32%, 43.15% and 50.00%, respectively. The effective number of codons converged at 51-53 codons. When plotted against total GC content or GC3, codon usage was distributed in relation to GC3 biases. Relative synonymous codon usage for each codon revealed a single major trend, which was highly correlated with GC content at the third position when codons began with A or U at the first two positions. In codons beginning with G or C base at the first two positions, the G or C base rarely occurred at the third position. These results suggest that codon usage is shaped by a bias towards G or C at the third base, and that this is affected by the first and second bases.     

Programmed translational -1 frameshifting on hexanucleotide motifs and the wobble properties of tRNAs

Programmed -1 ribosomal frameshifting, involving tRNA re-pairing from an AAG codon to an AAA codon, has been reported to occur at the sequences CGA AAG and CAA AAG. In this study, using the recoding region of insertion sequence IS3, we have investigated the influence on frameshifting in Escherichia coli of the first codon of this type of motif by changing it to all other NNA codons. Two classes of NNA codons were distinguished, depending on whether they favor or limit frameshifting. Their degree of shiftiness is correlated with wobble propensity, and base 34 modification, of their decoding tRNAs. A more flexible anticodon loop very likely makes the tRNAs with extended wobble more prone to liberate the third codon base, A, for re-pairing of tRNALys in the -1 frame.     

Wobble position modified nucleosides evolved to select transfer RNA codon recognition: a modified-wobble hypothesis.

While recognized that some wobble exists in the base pairing of the first base of the tRNA anticodon with the third of the codon, specific base modifications have evolved to select particular codons. This modified-wobble theory would be exemplified by a single codon recognition imposed on the anticodon by modification of the tRNA wobble position nucleoside.     

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.     

Codon catalog usage and the genome hypothesis.

Frequencies for each of the 61 amino acid codons have been determined in every published mRNA sequence of 50 or more codons. The frequencies are shown for each kind of genome and for each individual gene. A surprising consistency of choices exists among genes of the same or similar genomes. Thus each genome, or kind of genome, appears to possess a "system" for choosing between codons. Frameshift genes, however, have widely different choice strategies from normal genes. Our work indicates that the main factors distinguishing between mRNA sequences relate to choices among degenerate bases. These systematic third base choices can therefore be used to establish a new kind of genetic distance, which reflects differences in coding strategy. The choice patterns we find seem compatible with the idea that the genome and not the individual gene is the unit of selection. Each gene in a genome tends to conform to its species' usage of the codon catalog; this is our genome hypothesis.     

Comparison of base composition and codon usage in insect mitochondrial genomes

Insects, the most biodiverse taxonomic group, have high AT content in their mitochondrial genomes. Although codon usage tends to be AT-rich, base composition and codon usage of mitochondrial genomes may vary among taxa. Thus, we compare base composition and codon usage patterns of 49 insect mitochondrial genomes. For protein coding genes, AT content is as high as 80% in the Hymenoptera and Lepidoptera and as low as 72% in the Orthopotera. The AT content is high at positions 1 and 3, but A content is low at position 2. A close correlation occurs between codon usage and tRNA abundance in nuclear genomes. Optimal codons can pair well with the antr codons of the most abundant tRNAs. One tRNA gene translates a synonymous codon family in vertebrate mitochondrial genomes and these tRNA anticodons can pair with optimal codons. However, optimal codons cannot pair with anticodons in mtDNA ofCochiiomyia hominivorax (Dipteral: CaLliphoridae). Ten optimal codons cannot pair with tRNA anticodons in all 49 insect mitochondrial genomes; non-optimal codon-anticodon usage is common and codon usage is not influenced by tRNA abundance.     

Codon usage tabulated from international DNA sequence databases: status for the year 2000

The frequencies of each of the 257 468 complete protein coding sequences (CDSs) have been compiled from the taxonomical divisions of the GenBank DNA sequence database. The sum of the codons used by 8792 organisms has also been calculated. The data files can be obtained from the anonymous ftp sites of DDBJ, Kazusa and EBI. A list of the codon usage of genes and the sum of the codons used by each organism can be obtained through the web site http://www.kazusa.or.jp/codon/ . The present study also reports recent developments on the WWW site. The new web interface provides data in the CodonFrequency-compatible format as well as in the traditional table format. The use of the database is facilitated by keyword based search analysis and the availability of codon usage tables for selected genes from each species. These new tools will provide users with the ability to further analyze for variations in codon usage among different genomes.     

原核表达中优化起始密码下游序列的软件设计与实现

在原核表达中影响外源基因表达效率的因素有很多,这其中翻译起始效率起了非常重要的作用。翻译起始效率又主要受SD序列、SD序列与起始密码子之间的间距、DB(DownstreamBox)序列、mRNA翻译起始区(TIR)的二级结构和稀有密码子等因素的影响。主要针对DB序列和5′端稀有密码子的优化设计了软件。通过计算机对序列进行分析比对后,按照匹配碱基数、匹配位置、密码子使用频率平均值的顺序进行排序,给出一些优化序列,并给出了软件算法。     

On the Evolution of the Standard Genetic Code: Vestiges of Critical Scale Invariance from the RNA World in Current Prokaryote Genomes

Herein two genetic codes from which the primeval RNA code could have originated the standard genetic code (SGC) are derived. One of them, called extended RNA code type I, consists of all codons of the type RNY (purine-any base-pyrimidine) plus codons obtained by considering the RNA code but in the second (NYR type) and third (YRN type) reading frames. The extended RNA code type II, comprises all codons of the type RNY plus codons that arise from transversions of the RNA code in the first (YNY type) and third (RNR) nucleotide bases. In order to test if putative nucleotide sequences in the RNA World and in both extended RNA codes, share the same scaling and statistical properties to those encountered in current prokaryotes, we used the genomes of four Eubacteria and three Archaeas. For each prokaryote, we obtained their respective genomes obeying the RNA code or the extended RNA codes types I and II. In each case, we estimated the scaling properties of triplet sequences via a renormalization group approach, and we calculated the frequency distributions of distances for each codon. Remarkably, the scaling properties of the distance series of some codons from the RNA code and most codons from both extended RNA codes turned out to be identical or very close to the scaling properties of codons of the SGC. To test for the robustness of these results, we show, via computer simulation experiments, that random mutations of current genomes, at the rates of 10⁻¹⁰ per site per year during three billions of years, were not enough for destroying the observed patterns. Therefore, we conclude that most current prokaryotes may still contain relics of the primeval RNA World and that both extended RNA codes may well represent two plausible evolutionary paths between the RNA code and the current SGC.     

Codon and base biases after the initiation codon of the open reading frames in the Escherichia coli genome and their influence on the translation efficiency

Nucleotide sequences around the boundaries of all open reading frames in the Escherichia coli whole genome were analyzed. Characteristic base biases were observed after the initiation codon and before the termination codon. We examined the effect of the base sequence after the initiation codon on the translation efficiency, by introducing mutations after the initiation codon of the E. coli dihydrofolate reductase (DHFR) gene, considering codon and base biases, and using in vitro and in vivo translation systems. In both assay systems, the two most frequent second codons, AAA and AAU, enhanced the translation efficiency compared with the wild type, whereas the effects of lower frequency codons were not significant. Experiments using 16S rRNA variants with mutations in the putative complementary sequence to the region downstream of the initiation codon showed that the translation efficiency of none of the DHFR mutants was affected. These results demonstrate that the statistically most frequent sequences for the second codon enhance translation efficiency, and this effect seems to be independent of base pairing between mRNA and 16S rRNA.     

Potential secondary structure at translation-initiation sites.

Since translational start codons also occur internally, more-complex features within mRNA must determine initiation. We compare the potential secondary structure of 123 prokaryotic mRNA start regions to that of regions coding for internal methionines. The latter display an unexpectedly-uniform, almost-periodic pattern of pairing potential. In contrast, sequences 5' to start codons have little self-pairing, and do not pair extensively with the proximal coding region. Pairing potential surrounding start codons was found to be less than half of that found near internal AUGs. In groups of random sequences where the distribution of nucleotides at each position, or of trinucleotides at each in-frame codon position, matched the observed natural distribution, there was no periodicity in the pairing potential of the internal sequences. Randomized internal sequences had less pairing: the ratio of pairing intensity between internals and starts was reduced from 2.0 to 1.6 by randomization. We propose that the transition from the relatively-unstructured start domains to the highly-structured internal sequences may be an important determinant of translational start-site recognition.     

Contextual constraints on synonymous codon choice

We have studied the statistical constraints on synonymous codon choice to evaluate various proposals regarding the origin of the bias in synonymous codon usage observed by Fiers et al. (1975), Air et al. (1976), Grantham et al. (1980) and others. We have determined the statistical dependence of the degenerate third base on either of its nearest neighbors in mitochondrial, prokaryotic, and eukaryotic coding sequences. We noted an increasing dependence of the third base on its nearest neighbors in moving from mitochrondria to prokaryotes to eukaryotes.A statistical model assuming random equiprobable selection of synonymous codons was found grossly adequate for the mitochondria, but totally indequate for prokaryotes and eukaryotes. A model assuming selection of synonymous codons reflecting a genomic strategy, i.e. the genome hypothesis of Grantham et al. (1980), gave a good approximation of the mitochondrial sequences. A statistical model which exactly maintains codon frequency, but allows the position of corresponding synonymous codons to vary was only grossly adequate for prokaryotes and totally inadequate for eukaryotes. The results of these simulations are consistent with the measures on experimental sequences and suggest that a “frequency constraint” model such as that of Grantham et al. (1980) may be an adequate explanation of the codon usage in mitochondria. However, in addition to this frequency constraint, there may be constraints on synonymous codon choice in prokaryotes due to codon context. Furthermore, any proposal to explain codon usage in eukaryotes must involve a constraint on the context of a codon in the sequence.     

mRNAs have greater negative folding free energies than shuffled or codon choice randomized sequences

An examination of 51 mRNA sequences in GenBank has revealed that calculated mRNA folding is more stable than expected by chance. Free energy minimization calculations of native mRNA sequences are more negative than randomized mRNA sequences with the same base composition and length. Randomization of the coding region of genes yields folding free energies of less negative magnitude than the original native mRNA sequence. Randomization of codon choice, while still preserving original base composition, also results in less stable mRNAs. This suggests that a bias in the selection of codons favors the potential formation of mRNA structures which contribute to folding stability.     

Factors that influence selection of coding resumption sites in translational bypassing: minimal conventional peptidyl-tRNA:mRNA pairing can suffice

This study investigates bypassing initiated from codons immediately 5' of a stop codon. The mRNA slips and is scanned by the peptidyl-tRNA for a suitable landing site, and standard decoding resumes at the next 3' codon. This work shows that landing sites with potentially strong base pairing between the peptidyl-tRNA anticodon and mRNA are preferred, but sites with little or no potential for Watson-Crick or wobble base pairing can also be utilized. These results have implications for re-pairing in ribosomal frameshifting. Shine-Dalgarno sequences in the mRNA can alter the distribution of landing sites observed. The bacteriophage T4 gene 60 nascent peptide, known to influence take-off in its native context, imposes stringent P-site pairing requirements, thereby limiting the number of suitable landing sites.     

An Extensive Study of Mutation and Selection on the Wobble Nucleotide in tRNA Anticodons in Fungal Mitochondrial Genomes

Two alternative hypotheses aim to predict the wobble nucleotide of tRNA anticodons in mitochondrion. The codon-anticodon adaptation hypothesis predicts that the wobble nucleotide of tRNA anticodon should evolve toward maximizing the Watson-Crick base pairing with the most frequently used codon within each synonymous codon family. In contrast, the wobble versatility hypothesis argues that the nucleotide at the wobble site should be occupied by a nucleotide most versatile in wobble pairing, i.e., the wobble site of the tRNA anticodon should be G for NNY codon families and U for NNR and NNN codon families (where Y stands for C or U, R for A or G, and N for any nucleotide). We examined codon usage and anticodon wobble sites in 36 fungal genomes to evaluate these two alternative hypotheses and identify exceptional cases that deserve new explanations. While the wobble versatility hypothesis is generally supported, there are interesting exceptions involving tRNA(Arg) translating the CGN codon family, tRNA(Trp) translating the UGR codon family, and tRNA(Met) translating the AUR codon family. Our results suggest that the potential to suppress stop codons, the historical inertia, and the conflict between translation initiation and elongation can all contribute to determining the wobble nucleotide of tRNA anticodons.