Comparison of codon usage and tRNAs in mitochondrial genomes of Candida species

To gain insight into the nature of the mitochondrial genomes (mtDNA) of different Candida species, the synonymous codon usage bias of mitochondrial protein coding genes and the tRNAs in C. albicans, C. parapsilosis, C. stellata, C. glabrata and the closely related yeast Saccharomyces cerevisiae were analyzed. Common features of the mtDNA in Candida species are a strong A+T pressure on protein coding genes, and insufficient mitochondrial tRNA species are encoded to perform protein synthesis. The wobble site of the anticodon is always U for the NNR (NNA and NNG) codon families, which are dominated by A-ending codons, and always G for the NNY (NNC and NNU) codon families, which is dominated by U-ending codons, and always U for the NNN (NNA, NNU, NNC and NNG) codon families, which are dominated by A-ending codons and U-ending codons. Patterns of synonymous codon usage of Candida species can be classified into three groups: (1) optimal codon-anticodon usage, Glu, Lys, Leu (translated by anti-codon UAA), Gln, Arg (translated by anti-codon UCU) and Trp are containing NNR codons. NNA, whose corresponding tRNA is encoded in the mtDNA, is used preferentially. (2) Non-optimal codon-anticodon usage, Cys, Asp, Phe, His, Asn, Ser (translated by anti-codon GCU) and Tyr are containing NNY codons. The NNU codon, whose corresponding tRNA is not encoded in the mtDNA, is used preferentially. (3) Combined codon-anticodon usage, Ala, Gly, Leu (translated by anti-codon UAG), Pro, Ser (translated by anti-codon UGA), Thr and Val are containing NNN codons. NNA (tRNA encoded in the mtDNA) and NNU (tRNA not encoded in the mtDNA) are used preferentially. In conclusion, we propose that in Candida species, codons containing A or U at third position are used preferentially, regardless of whether corresponding tRNAs are encoded in the mtDNA. These results might be useful in understanding the common features of the mtDNA in Candida species and patterns of synonymous codon usage.     

Translation efficiencies of synonymous codons are not always correlated with codon usage in tobacco chloroplasts

Codon usage in chloroplasts is different from that in prokaryotic and eukaryotic nuclear genomes. However, no experimental approach has been made to analyse the translation efficiency of individual codons in chloroplasts. We devised an in vitro assay for translation efficiencies using synthetic mRNAs, and measured the translation efficiencies of five synonymous codon groups in tobacco chloroplasts. Among four alanine codons (GCN, where N is U, C, A or G), GCU was the most efficient for translation, whereas the chloroplast genome lacks tRNA genes corresponding to GCU. Phenylalanine and tyrosine are each encoded by two codons (UUU/C and UAU/C, respectively). Phenylalanine UUC and tyrosine UAC were translated more than twice as efficiently than UUU and UAU, respectively, contrary to their codon usage, whereas translation efficiencies of synonymous codons for alanine, aspartic acid and asparagine were parallel to their codon usage. These observations indicate that translation efficiencies of individual codons are not always correlated with codon usage in vitro in chloroplasts. This raises an important issue for foreign gene expression in chloroplasts.     

Accurate translation of the genetic code depends on tRNA modified nucleosides

Transfer RNA molecules translate the genetic code by recognizing cognate mRNA codons during protein synthesis. The anticodon wobble at position 34 and the nucleotide immediately 3' to the anticodon triplet at position 37 display a large diversity of modified nucleosides in the tRNAs of all organisms. We show that tRNA species translating 2-fold degenerate codons require a modified U(34) to enable recognition of their cognate codons ending in A or G but restrict reading of noncognate or near-cognate codons ending in U and C that specify a different amino acid. In particular, the nucleoside modifications 2-thiouridine at position 34 (s(2)U(34)), 5-methylaminomethyluridine at position 34 (mnm(5)U(34)), and 6-threonylcarbamoyladenosine at position 37 (t(6)A(37)) were essential for Watson-Crick (AAA) and wobble (AAG) cognate codon recognition by tRNA(UUU)(Lys) at the ribosomal aminoacyl and peptidyl sites but did not enable the recognition of the asparagine codons (AAU and AAC). We conclude that modified nucleosides evolved to modulate an anticodon domain structure necessary for many tRNA species to accurately translate the genetic code.     

The preferential codon usages in variable and constant regions of immunoglobulin genes are quite distinct from each other.

The pattern of codon utilization in the variable and constant regions of immunoglobulin genes are compared. It is shown that, in these regions, codon utilizations are quite distinct from one another: For most degenerate codons, there is a selective bias that prefers C and/or G ending codons to U and/or A ending codons in the constant region compared with the bias in the variable region. This would strongly suggest that, in immunoglobulin genes, the bias in code word usage is determined by other factors than those concerning with the translational mechanism such as tRNA availability and codon-anticodon interaction. A possibility is also suggested that this differance of code word usage between them is due to the existence of secondary structure in the constant region but not in the variable region.     

Novel anticodon composition of transfer RNAs in Micrococcus luteus, a bacterium with a high genomic G + C content. Correlation with codon usage.

The number and relative amount of isoacceptor tRNAs for each amino acid in Micrococcus luteus, a Gram-positive bacterium with high genomic G + C content, have been determined by sequencing their anticodon loop and its adjacent regions and by selective labelling of tRNAs. Thirty-one tRNA species with 29 different anticodon sequences have been detected. All the tRNAs have G or C at the anticodon first position except for tRNA(ICGArg) and tRNA(NGASer), in response to the abundant usage of NNC and NNG codons. No tRNA with the anticodon UNN capable of translating codon NNA has been detected, in accordance with a very low or zero usage of NNA codons. The relative amount of isoacceptor tRNAs for an amino acid determined by selective labelling strongly correlates with usage of the corresponding codons. On the basis of these and other observations in this and other eubacterial species, we conclude that the relative amount and anticodon composition of isoacceptor tRNA species are flexible, and their changes are mainly adaptive phenomena that have been primarily affected by codon usage, which in turn is affected by directional mutation pressure.     

The genetic code as a periodic table: algebraic aspects

The systematics of indices of physico-chemical properties of codons and amino acids across the genetic code are examined. Using a simple numerical labelling scheme for nucleic acid bases, A=(-1,0), C=(0,-1), G=(0,1), U=(1,0), data can be fitted as low order polynomials of the six coordinates in the 64-dimensional codon weight space. The work confirms and extends the recent studies by Siemion et al. (1995. BioSystems 36, 231-238) of the conformational parameters. Fundamental patterns in the data such as codon periodicities, and related harmonics and reflection symmetries, are here associated with the structure of the set of basis monomials chosen for fitting. Results are plotted using the Siemion one-step mutation ring scheme, and variants thereof. The connections between the present work, and recent studies of the genetic code structure using dynamical symmetry algebras, are pointed out.     

Comparative analysis of codon usage bias in Crenarchaea and Euryarchaea genome reveals differential preference of synonymous codons to encode highly expressed ribosomal and RNA polymerase proteins

The present study was undertaken to investigate the pattern of optimal codon usage in Archaea. Comparative analysis was executed to understand the pattern of codon usage bias between the high expression genes (HEG) and the whole genomes in two Archaeal phyla, Crenarchaea and Euryarchaea. The G + C% of the HEG was found to be less in comparison to the genome G + C% in Crenarchaea, whereas reverse was the case in Euryarchaea. The preponderance of U/A ending codons that code for HEG in Crenarchaea was in sharp contrast to the C/G ended ones in Euryarchaea. The analysis revealed prevalence of U-ending codons even within the WWY (nucleotide ambiguity code) families in Crenarchaea vis-à-vis Euryarchaea, bacteria and Eukarya. No plausible interpretation of the observed disparity could be made either in the context of tRNA gene composition or genome G + C%. The results in this study attested that the preferential biasness for codons in HEG of Crenarchaea might be different from Euryarchaea. The main highlights are (i) varied CUB in the HEG and in the whole genomes in Euryarchaea and Crenarchaea. (ii) Crenarchaea was found to have some unusual optimal codons (OCs) compared to other organisms. (iii) G + C% (and GC₃) of the HEG were different from the genome G + C% in the two phyla. (iv) Genome G + C% and tRNA gene number failed to explain CUB in Crenarchaea. (v) Translational selection is possibly responsible for A + T rich OCs in Crenarchaea.     

Codon sextets with leading role of serine create “ideal” symmetry classification scheme of the genetic code

The standard classification scheme of the genetic code is organized for alphabetic ordering of nucleotides. Here we introduce the new, “ideal” classification scheme in compact form, for the first time generated by codon sextets encoding Ser, Arg and Leu amino acids. The new scheme creates the known purine/pyrimidine, codon–anticodon, and amino/keto type symmetries and a novel A + U rich/C + G rich symmetry. This scheme is built from “leading” and “nonleading” groups of 32 codons each. In the ensuing 4 × 16 scheme, based on trinucleotide quadruplets, Ser has a central role as initial generator. Six codons encoding Ser and six encoding Arg extend continuously along a linear array in the “leading” group, and together with four of six Leu codons uniquely define construction of the “leading” group. The remaining two Leu codons enable construction of the “nonleading” group. The “ideal” genetic code suggests the evolution of genetic code with serine as an initiator.     

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.     

Codon usage in Homo sapiens: evidence for a coding pattern on the non-coding strand and evolutionary implications of dinucleotide discrimination

This study reports the analysis of codon usage in 35 complete Homo sapiens genes. Both codon frequency and inter-codon interference exhibit patterns of evolutionary interest. There is a significant positive correlation between the frequency with which a given codon is used and the frequency with which its complement is used. Since the frequency of appearance of the complementary codon on the coding strand is equal to the frequency of appearance of the original codon on the non-coding strand, in the same phase, the non-coding strand is found to resemble the coding strand in triplet composition. The same effect has been observed in Escherichia coli. This preference for the use of certain complementary triplets as codons suggests that the evolution of the use of the genetic code depended to some extent upon the double-stranded nature of the coding material. In addition, the effect of discrimination against the use of two dinucleotides, CpG and UpA, is observed in codon usage and also in adjacent codon interference. Codons beginning with G, or A, are unlikely to be preceded by codons ending in C, or U, respectively. Consideration of codon assignment in the genetic code together with the observed CpG infrequency suggests that the evolution of the code may have been influenced by conditions in which the use of CpG dinucleotides was unfavorable. The infrequent use of UpA dinucleotides can be explained as the result of frameshift mutation during gene evolution.     

Effect of Correlated tRNA Abundances on Translation Errors and Evolution of Codon Usage Bias

Despite the fact that tRNA abundances are thought to play a major role in determining translation error rates, their distribution across the genetic code and the resulting implications have received little attention. In general, studies of codon usage bias (CUB) assume that codons with higher tRNA abundance have lower missense error rates. Using a model of protein translation based on tRNA competition and intra-ribosomal kinetics, we show that this assumption can be violated when tRNA abundances are positively correlated across the genetic code. Examining the distribution of tRNA abundances across 73 bacterial genomes from 20 different genera, we find a consistent positive correlation between tRNA abundances across the genetic code. This work challenges one of the fundamental assumptions made in over 30 years of research on CUB that codons with higher tRNA abundances have lower missense error rates and that missense errors are the primary selective force responsible for CUB.     

Codon usage in higher plants, green algae, and cyanobacteria

Codon usage is the selective and nonrandom use of synonymous codons by an organism to encode the amino acids in the genes for its proteins. During the last few years, a large number of plant genes have been cloned and sequenced, which now permits a meaningful comparison of codon usage in higher plants, algae, and cyanobacteria. For the nuclear and organellar genes of these organisms, a small set of preferred codons are used for encoding proteins. Codon usage is different for each genome type with the variation mainly occurring in choices between codons ending in cytidine (C) or guanosine (G) versus those ending in adenosine (A) or uridine (U). For organellar genomes, chloroplastic and mitochrondrial proteins are encoded mainly with codons ending in A or U. In most cyanobacteria and the nuclei of green algae, proteins are encoded preferentially with codons ending in C or G. Although only a few nuclear genes of higher plants have been sequenced, a clear distinction between Magnoliopsida (dicot) and Liliopsida (monocot) codon usage is evident. Dicot genes use a set of 44 preferred codons with a slight preference for codons ending in A or U. Monocot codon usage is more restricted with an average of 38 codons preferred, which are predominantly those ending in C or G. But two classes of genes can be recognized in monocots. One set of monocot genes uses codons similar to those in dicots, while the other genes are highly biased toward codons ending in C or G with a pattern similar to nuclear genes of green algae. Codon usage is discussed in relation to evolution of plants and prospects for intergenic transfer of particular genes.     

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.     

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 and tRNA content in unicellular and multicellular organisms

Choices of synonymous codons in unicellular organisms are here reviewed, and differences in synonymous codon usages between Escherichia coli and the yeast Saccharomyces cerevisiae are attributed to differences in the actual populations of isoaccepting tRNAs. There exists a strong positive correlation between codon usage and tRNA content in both organisms, and the extent of this correlation relates to the protein production levels of individual genes. Codon-choice patterns are believed to have been well conserved during the course of evolution. Examination of silent substitutions and tRNA populations in Enterobacteriaceae revealed that the evolutionary constraint imposed by tRNA content on codon usage decelerated rather than accelerated the silent-substitution rate, at least insofar as pairs of taxonomically related organisms were examined. Codon-choice patterns of multicellular organisms are briefly reviewed, and diversity in G+C percentage at the third position of codons in vertebrate genes--as well as a possible causative factor in the production of this diversity--is discussed.      

Silent substitutions predictably alter translation elongation rates and protein folding efficiencies

Genetic code redundancy allows most amino acids to be encoded by multiple codons that are non-randomly distributed along coding sequences. An accepted theory explaining the biological significance of such non-uniform codon selection is that codons are translated at different speeds. Thus, varying codon placement along a message may confer variable rates of polypeptide emergence from the ribosome, which may influence the capacity to fold toward the native state. Previous studies report conflicting results regarding whether certain codons correlate with particular structural or folding properties of the encoded protein. This is partly due to different criteria traditionally utilized for predicting translation speeds of codons, including their usage frequencies and the concentration of tRNA species capable of decoding them, which do not always correlate. Here, we developed a metric to predict organism-specific relative translation rates of codons based on the availability of tRNA decoding mechanisms: Watson-Crick, non-Watson-Crick or both types of interactions. We determine translation rates of messages by pulse-chase analyses in living Escherichia coli cells and show that sequence engineering based on these concepts predictably modulates translation rates in a manner that is superior to codon usage frequency, which occur during the elongation phase, and significantly impacts folding of the encoded polypeptide. Finally, we demonstrate that sequence harmonization based on expression host tRNA pools, designed to mimic ribosome movement of the original organism, can significantly increase the folding of the encoded polypeptide. These results illuminate how genetic code degeneracy may function to specify properties beyond amino acid encoding, including folding.     

Complete DNA Sequence of the Mitochondrial Genome of the Black Chiton, Katharina Tunicata

The DNA sequence of the 15,532-base pair (bp) mitochondrial DNA (mtDNA) of the chiton Katharina tunicata has been determined. The 37 genes typical of metazoan mtDNA are present: 13 for protein subunits involved in oxidative phosphorylation, 2 for rRNAs and 22 for tRNAs. The gene arrangement resembles those of arthropods much more than that of another mollusc, the bivalve Mytilus edulis. Most genes abut directly or overlap, and abbreviated stop codons are inferred for four genes. Four junctions between adjacent pairs of protein genes lack intervening tRNA genes; however, at each of these junctions there is a sequence immediately adjacent to the start codon of the downstream gene that is capable of forming a stem-and-loop structure. Analysis of the tRNA gene sequences suggests that the D arm is unpaired in tRNA(ser(AGN)), which is typical of metazoan mtDNAs, and also in tRNA(ser(UCN)), a condition found previously only in nematode mtDNAs. There are two additional sequences in Katharina mtDNA that can be folded into structures resembling tRNAs; whether these are functional genes is unknown. All possible codons except the stop codons TAA and TAG are used in the protein-encoding genes, and Katharina mtDNA appears to use the same variation of the mitochondrial genetic code that is used in Drosophila and Mytilus. Translation initiates at the codons ATG, ATA and GTG. A + T richness appears to have affected codon usage patterns and, perhaps, the amino acid composition of the encoded proteins. A 142-bp non-coding region between tRNA(glu) and CO3 contains a 72-bp tract of alternating A and T.     

Possible evolutionary steps in the genetic code

In mammalian mitochondrial codes, fourfold codons wobble-pair with UNN anticodons so that U wobbles with U, C, A and G. Twofold pyrimidine-terminated codons pair with GNN and twofold purine-terminated codons pair with UNN. These properties enable a prediction to be made for evolution of the universal genetic code. It was postulated (1) that an archetypal code of 16 quartets coded for 15 amino acids. If this code used UNN anticodons, then duplication of tRNA genes, followed by mutations in the anticodons and aminoacylation sites, would give rise to the present universal code.