<title>Background</title> <p>The wobble versatility hypothesis [##REF##6997870##1##, ####REF##14715921##2##, ##REF##15177692##3##, ##REF##6447871##4##, ##UREF##0##5##, ##REF##15716092##6####15716092##6##], abbreviated as WVH, states that the wobble site of tRNA anticodon should have G for NNY codons (where Y stands for C or U and N for any nucleotide) because G can pair with both C and U in RNA, and should have U for NNR to pair with both A and G. For NNN codon families, the wobble site should be U because U is known to be the most versatile in wobble-pairing [##UREF##1##7##, ####REF##2416594##8##, ##REF##7658467##9##, ##REF##11675591##10##, ##REF##3502716##11##, ##REF##6932013##12####6932013##12##]. In contrast, the codon-anticodon adaptation hypothesis, or CAAH for short, invokes the codon usage bias as a determining factor, i.e., the wobble site of tRNA anticodon should co-evolve with codon usage so that the nucleotide in the wobble site of tRNA anticodon should match the most abundant codon in a synonymous codon family [##REF##15716092##6##,##REF##2434856##13##, ####REF##1752426##14##, ##REF##9584084##15####9584084##15##]. The association between the major codon and the anticodon of the most abundant tRNA has been documented in <italic>Escherichia coli </italic>[##REF##6175758##16##,##REF##6760125##17##], <italic>Saccharomyces cerevisiae </italic>[##REF##7037777##18##], and other species and organelles [##REF##9584084##15##,##REF##12930740##19##, ####REF##8754221##20##, ##REF##8913770##21##, ##REF##18048402##22####18048402##22##].</p> <p>Here we develop a general hypothesis of codon-anticodon adaptation based on an analysis of wobble costs, and derive its predictions that can be tested by genomic data. The wobble cost may be viewed as reduction in decoding efficiency and accuracy because such reduction would be selected against over evolutionary time. We will refer to this new general hypothesis based on wobble cost as WCH (for wobble cost hypothesis). The two traditional hypotheses, CAAH and WVH, will be shown to be special forms of WCH.</p> <p>Following the shorthand notation of Ogle et al. [##REF##15952884##23##], I designate the translation cost through wobble base-pairing between nucleotide i at the third codon position of a codon and the nucleotide j at the wobble site of tRNA anticodon as M<sub>i3:j </sub>(where M is for wobble cost. The letter C would be more suitable to represent cost but it may confuse with the nucleotide C). We assume M<sub>i3:j </sub>= 0 if nucleotides i and j form Watson-Crick base pairing. The reason for this assumption is that C:G and A:U pairs have not been found to contribute to ribosomal stalling (which reduces translation efficiency) or amino acid mis-incorporation (which reduces translation accuracy) although almost all non-Watson-Crick base pairings have been shown to reduce translation efficiency and accuracy. We define M<sub>Y3:U </sub>as the wobble cost between a wobble U at the tRNA anticodon and a C or U at the third codon position. For simplicity, we also assume M<sub>A3:A </sub>= M<sub>C3:A </sub>= M<sub>G3:A </sub>= M<sub>A3:C </sub>= M<sub>C3:C </sub>= M<sub>U3:C </sub>= M<sub>A3:G </sub>= M<sub>G3:G </sub>= M<sub>O </sub>(where the subscript O is for "other", i.e., M<sub>O </sub>is the wobble cost between base pairs that are not Watson-Crick base pairs, not U3:G or G3:U base pairs and not Y3:U pairs). In general, M<sub>U3:G </sub>and M<sub>G3:U </sub>are expected to be smaller than M<sub>C3:U</sub>, M<sub>U3:U</sub>, and M<sub>O </sub>because G and U can base-pair in RNA, M<sub>Y3:U </sub>is expected to be smaller than M<sub>O </sub>because a wobble U at tRNA anticodon is known to be the most versatile in wobble-pairing [##UREF##1##7##, ####REF##2416594##8##, ##REF##7658467##9##, ##REF##11675591##10##, ##REF##3502716##11##, ##REF##6932013##12####6932013##12##].</p> <p>A classical study of U:G wobble pairs [##REF##724522##24##] suggests a preference for the U being at the 3' end rather than at the 5' end, which implies that U3:G wobble pair is energetically more favorable than G3:U wobble pair [##REF##17496902##25##]. Subsequent studies have shown that, while the U3:G wobble pair occurs on the ribosome, the unmodified G3:U wobble pair does not [##REF##14715921##2##,##REF##15952884##23##,##REF##10668794##26##]. These findings suggest that M<sub>U3:G </sub>may be smaller than M<sub>G3:U</sub>, although its generality is unknown.</p> <title>Two-fold NNY and NNR codon families</title> <p>First consider the NNY codon family where Y is either C or U. Designate the number of C-ending and U-ending synonymous codons by N<sub>C </sub>and N<sub>U</sub>, respectively, and the total cost of wobble pairing as M<sub>wG </sub>when the wobble site of the anticodon is G, and as M<sub>wA </sub>when the wobble nucleotide is A (we do not need to consider the case when the wobble site is U or C for NNY codon families because such cases have never been observed and because a tRNA with a wobble U or C to translate NNC and NNU codons is against physiochemical reasons). We now express the total cost M<sub>wG </sub>and M<sub>wA </sub>as</p> <p></p> <p>Note that M<sub>C3:G </sub>= M<sub>U3:A </sub>= 0 according to our definition. The dependence of M<sub>wG </sub>and M<sub>wA </sub>on the relative frequencies of N<sub>C</sub>, expressed as proportion of C-ending codons in the NNY codon family (P<sub>C</sub>), is graphically shown in Fig. ##FIG##0##1##, with M<sub>O </sub>assumed to equal 2•;M<sub>U3:G</sub>. The condition for M<sub>wG </sub>= M<sub>wA</sub>, i.e., when the wobble site of the tRNA anticodon can take either G or A without a fitness differential, is</p> <p></p> <p>In the scenario in Fig. ##FIG##0##1## with M<sub>O </sub>= 2M<sub>U3:G</sub>, the condition for M<sub>wG </sub>= M<sub>wA </sub>is P<sub>C </sub>= 1/3. Thus, when P<sub>C </sub>drifts above 1/3, it would become more beneficial (less costly) to have a G at the wobble nucleotide of tRNA anticodon. However, when PC drifts below 1/3, then a wobble A would be more advantageous and should be favored by natural selection. Naturally, if M<sub>U3:G </sub>= 0, then we should always have a wobble G at the tRNA anticodon regardless of codon usage.</p> <p>Now we consider three special cases with reference to Fig. ##FIG##0##1## and the two traditional hypotheses, CAAH and WVH. First, if N<sub>C </sub>= N<sub>U </sub>= N/2, then CAAH, which predicts that the wobble nucleotide of tRNA anticodon should form Watson-Crick base pair with the most abundant codon, has no prediction at all because the two codons are equally abundant. In contrast, if we assume M<sub>U3:G </sub>&lt; M<sub>O </sub>as is depicted in Fig. ##FIG##0##1##, we have M<sub>wG </sub>&lt; M<sub>wA</sub>, and WCH predicts that the anticodon wobble site should be a G. This is the same prediction as WVH.</p> <p>Second, if N<sub>C </sub>&gt; N<sub>U</sub>, then M<sub>wG </sub>&lt; M<sub>wA </sub>assuming M<sub>U3:G </sub>≤ M<sub>O</sub>, and WCH predicts that the anticodon wobble site should be occupied by a G. This prediction is shared by both CAAH and WVH.</p> <p>Third, when N<sub>C </sub>&lt; N<sub>U </sub>but N<sub>C </sub>&gt; N<sub>U</sub>·M<sub>U3:G</sub>/M<sub>O</sub>, then WCH will still predict a G at the wobble nucleotide of tRNA anticodon because M<sub>wG </sub>is still smaller than M<sub>wA</sub>. Take the scenario in Fig. ##FIG##0##1## for example, when N<sub>C </sub>&lt; N<sub>U </sub>but P<sub>C </sub>&gt; 1/3 (which correspond to the shaded area in Fig. ##FIG##0##1##), we have M<sub>wG </sub>&lt; M<sub>wA </sub>(Fig. ##FIG##0##1##), so natural selection should favor a G at the wobble site of tRNA anticodon. WVH happens to have the same prediction. However, CAAH will predict an A at the wobble site of tRNA anticodon if U-ending codons are more abundant than C-ending codons. This is in contrast to WVH, i.e. CAAH is inapplicable when P<sub>C </sub>is within the shaded range in Fig. ##FIG##0##1##.</p> <p>Fourth, when N<sub>C </sub>&lt;&lt; N<sub>U</sub>, especially in the extreme case when N<sub>C </sub>= 0, then Eq. (1) is reduced to M<sub>wG </sub>= N<sub>U• </sub>M<sub>U3:G </sub>and M<sub>wA </sub>= 0. Because now M<sub>wG </sub>&gt; M<sub>wA</sub>, WCH predicts an A at the anticodon wobble site. In this case, WVH would still predict a G at the wobble site because it ignores the codon frequencies (i.e., it ignores the relative magnitude of N<sub>C </sub>and N<sub>U</sub>), but CAAH would predict an A at the wobble site, which is the same prediction as WCH. Only in this particular case when N<sub>C </sub>&lt;&lt; N<sub>U </sub>can CAAH and WVH be clearly differentiated. As is depicted in Fig. ##FIG##0##1##, the chance for CAAH to be supported decreases as M<sub>U3:G </sub>decreases relative to M<sub>O</sub>. CAAH would have little chance to be unambiguously supported if M<sub>U3:G </sub>is close to 0. Thus, as long as M<sub>U3:G </sub>is small relative to M<sub>O</sub>, we should expect CAAH to be less supported than WVH. A recent study shows that CAAH indeed receives much less support than WVH when both were tested against fungal mitochondrial genomic data [##REF##18401633##27##].</p> <p>The simple Eq. (1) offers a way to assess relative costs of different wobble pairs. For example, if nature has chosen a G at the wobble site of the tRNA anticodon, then we may infer that M<sub>wG </sub>&lt; M<sub>wA</sub>. So, from Eq. (1), we have</p> <p></p> <p>Suppose we have an NNY codon family translated by a tRNA with a G at the wobble site of the anticodon. If M<sub>U3:G </sub>is large, then there would be strong selection favoring C-ending codons against U ending codons. So the ratio of N<sub>C</sub>/N<sub>U </sub>may be very large. In such a case the estimated M<sub>U3:G</sub>/M<sub>O </sub>ratio, being smaller than a very large value, is not informative. However, in different codon families, mutation pressure may allow N<sub>U </sub>to drift up relative to N<sub>C</sub>. Suppose we have three NNY codon families each translated by a tRNA with a G at the anticodon wobble site. If the N<sub>C</sub>/N<sub>U </sub>ratio for the three NNY codon families are 100/2, 80/100 and 200/20, respectively, we may infer that M<sub>U3:G</sub>/M<sub>O </sub>&lt; 80/100 (because it is the most precise). Alternatively, with many N<sub>C</sub>/N<sub>U </sub>ratios from many NNY codon families translated by tRNA with a G at its anticodon wobble site, we may compute the lower 95% confidence limit of the N<sub>C</sub>/N<sub>U </sub>ratio (LCL<sub>95.G </sub>where the subscript G indicates tRNA with a G at the anticodon wobble site) and infer that M<sub>U3:G</sub>/M<sub>C3:A </sub>&lt; LCL<sub>95.G</sub>.</p> <p>If we always have very large N<sub>C</sub>/N<sub>U </sub>ratios, we may infer that selection against U-ending codons must be strong, with little chance for mutation to elevate N<sub>U</sub>. This is a strong indication of a large M<sub>U3:G</sub>. Along the same line of reasoning, we may infer that M<sub>U3:G </sub>is very small if N<sub>U </sub>can often as large as, or even larger than, N<sub>C</sub>.</p> <p>Similarly, if nature has chosen an A at the wobble site of tRNA anticodon, then we may infer M<sub>wG </sub>&gt; M<sub>wA</sub>, so</p> <p></p> <p>We can apply exactly the same rationale for the NNR codon family leading to parallel conclusions. For example, if nature has chosen a U at the wobble site of the tRNA anticodon, then we may infer that M<sub>wU </sub>&lt; M<sub>wC</sub>, so that</p> <p></p> <p>Similarly, if M<sub>G3:U </sub>is very large, then there should be strong selection against G-ending codons in favor of A-ending codons. This will produce large N<sub>A</sub>/N<sub>G </sub>ratios. In contrast, a large N<sub>G </sub>comparable to N<sub>A </sub>indicates a very small M<sub>G3:U </sub>cost.</p> <p>Given previous studies indicating that U3:G is energetically much more favorable than G3:U [##REF##14715921##2##,##REF##15952884##23##, ####REF##724522##24##, ##REF##17496902##25##, ##REF##10668794##26####10668794##26##], we should expect M<sub>U3:G </sub>&lt; M<sub>G3:U</sub>. The reasoning above paves the way for us to test whether this is generally true among the genetically diverse fungal species.</p> <title>Three-fold AUH (Isoleucine) codons</title> <p>Designate codons ending with A, C, and U as N<sub>A</sub>, N<sub>C </sub>and N<sub>U</sub>, respectively. The wobble cost of having an A, G, C, or U at the wobble site of tRNA anticodon is</p> <p></p> <p>It is obvious that M<sub>wC </sub>is always greater than M<sub>wA</sub>, so we should never find a C at the wobble site of a tRNA<sup>Ile </sup>anticodon, i.e., we can disregard M<sub>wC</sub>. If nature has chosen G at the wobble site of tRNA<sup>Ile </sup>anticodon, then we may infer that M<sub>wG </sub>is the smallest. From M<sub>wG </sub>&lt; M<sub>wA </sub>and M<sub>wG </sub>&lt; M<sub>wU</sub>, we have</p> <p></p> <p>where the assumption is made on the basis of previous observations that U is generally the most versatile in wobble-pairing among the four nucleotides [##UREF##1##7##, ####REF##2416594##8##, ##REF##7658467##9##, ##REF##11675591##10##, ##REF##3502716##11##, ##REF##6932013##12####6932013##12##].</p> <title>Four-fold NNN codons</title> <p>Designate the number of codons ending with A, C, G, and U as N<sub>A</sub>, N<sub>C</sub>, N<sub>G</sub>, and N<sub>U</sub>, respectively. The wobble costs involving an A, C, G or U at the wobble site of tRNA anticodon are, respectively,</p> <p></p> <p>If nature has chosen a U at the wobble site of the tRNA anticodon, then Eq. (8) does not give us any simple inequality to estimate the cost ratios. However, if nature has chosen an A at the wobble site of the tRNA anticodon, then from M<sub>wA </sub>&lt; M<sub>wG </sub>and M<sub>wA </sub>&lt; M<sub>wC</sub>, we can infer</p> <p></p> <p>In what follows, we estimate M<sub>U3:G</sub>/M<sub>O</sub>, M<sub>G3:U</sub>/M<sub>O</sub>, and M<sub>U3:C</sub>/M<sub>O </sub>by using fungal mitochondrial genomic data. Cells in fungal species are generally rapid-replicating which necessitates efficient translation. Rapidly replicating unicellular organisms are theoretically expected to be under strong selection to increase the rate of biosynthesis [##REF##9584084##15##,##REF##7476129##28##] and they typically exhibit strong codon-anticodon adaptation [##UREF##2##29##]. Thus, fungal species should be ideal for evaluating evolutionary hypothesis on codon-anticodon adaptation.</p>
<title>Methods</title> <p>We retrieved 36 fungal mitochondrial genomes (Table ##TAB##0##1##) by using NCBI Entrez. Three different genetic codes are used in different fungal genomes. Among the 36 fungal mitochondrial genomes, seven genomes use translation table 3, 27 genomes use translation table 4, and two genomes use translation table 16 (Table ##TAB##0##1##). When results are similar among genomes using the same translation table, only results from a representative genome are presented. The number of codon families supporting CAAH (N<sub>CAAH</sub>) and WVH (N<sub>WVH</sub>) is compiled following the following rationale [##REF##18401633##27##]. Suppose a lysine (Lys) codon family has 20 AAA and 60 AAG codons. WVH would ignore the codon usage bias and predict a wobble U in the tRNA<sup>Lys </sup>anticodon because U can pair with both A and G, whereas CAAH would predict a wobble C in the tRNA<sup>Lys </sup>anticodon to maximize the Watson-Crick match with the more frequent G-ending codons. If the tRNA<sup>Lys </sup>anticodon is found to have a wobble U, then WVH is supported; if a wobble C is found, then CAAH is supported. If we have 60 AAA codons and 20 AAG codons and if tRNA<sup>Lys </sup>anticodon has a wobble U, then both hypotheses are supported, i.e., they are indistinguishable and are not included in Table ##TAB##0##1##. The methionine codon families are not included in Table ##TAB##0##1## but discussed in detail elsewhere [##REF##18401633##27##,##REF##17311091##30##].</p> <p>The tRNA and CDS sequences were extracted and analyzed by using DAMBE [##REF##11535656##31##,##UREF##3##32##]. The CDS-derived codon usage is also computed with DAMBE. The anticodon in almost all tRNA sequences from all species share the regular feature of being flanked by two nucleotides on either side to form a loop that is held together by a stem. For example, the anticodon loop (AC loop) of the tRNA<sup>Arg </sup>genes translating CGN codons in <italic>Epidermophyton floccosum </italic>is 28CGUG<bold>UU</bold>ACG<bold>GC</bold>CACG42, where the starting and ending numbers indicate the position of the AC loop in the tRNA sequence, with the anticodon 5'-ACG-3' (matching codon CGU) flanked by two nucleotides on either side (in bold) to form a loop that is held together by a stem made of the first and the last four nucleotides. Similarly, the other tRNA<sup>Arg </sup>translating AGR codons is 25AAAAUA<bold>CU</bold>UCU<bold>AA</bold>UAUUUU43, with the AC loop held together by a six-base stem. However, some tRNA sequences have a suspicious AC loop and DAMBE will flag them out. The AC loop is then identified by aligning the tRNA sequences against other isoaccepting tRNA sequences with a regular AC loop [##REF##15716092##6##]. Some tRNA anticodon loop has the anticodon flanked by three instead of two nucleotides. For example, the anticodon loop in tRNA<sup>Leu </sup>in the mitochondrial genome of <italic>Kluyveromyces thermotolerans </italic>is GAUAC<bold>UCU</bold>UAA<bold>GAU</bold>GUAUU, with the anticodon UAA flanked by three nucleotides (in bold) on both sides. There are a few tRNA sequences in which anticodon loop cannot be identified.</p> <p>Some mitochondrial genomes in GenBank are annotated incorrectly. For example, tRNA<sup>Pro </sup>in the mitochondrial genome of <italic>Ashbya gossypii </italic>ATCC 10895 has an anticodon of UGG matching codon CCA (the most frequently used proline codon), but the GenBank file (<ext-link ext-link-type="gen" xlink:href="NC_005789">NC_005789</ext-link>) annotated the anticodon to match codon CCU.</p> <p>A few fungal mitochondrial genomes do not have a complete set of tRNA genes. For example, the mitochondrial genomes of <italic>Hyaloraphidium curvatum </italic>and <italic>Harpochytrium sp</italic>. JEL94 have seven and eight tRNA genes, respectively, and consequently will need tRNA import from the nuclear genome. This may cause complication in analyzing codon-anticodon adaptation. However, removing such genomes does not alter the conclusions.</p> <p>Some species exhibit extreme avoidance of certain codon families. For example, <italic>Ashbya gossypii </italic>ATCC 10895 codes Arg with only AGR codons without using any CGN codons. In contrast, <italic>Hyaloraphidium curvatum </italic>codes Arg with only CGN codons without using any AGR codons. Such avoidance of certain codon families would facilitate the evolutionary loss of the associated tRNA [##REF##14708576##33##, ####REF##15781705##34##, ##REF##17541678##35####17541678##35##], although it is not always clear whether the avoidance is the cause or the consequence of the loss of the associated tRNA.</p> <p>We computed relative synonymous codon usage, or RSCU [##REF##3526280##36##], as a measure of codon usage bias within a codon family by using DAMBE [##REF##11535656##31##,##UREF##3##32##]. Some coding sequences are incomplete. For example, the cox1 CDS in <italic>Aspergillus niger </italic>is annotated as "join(&lt;19768..20614,21640..22495)". The first two nucleotides (i.e., at positions 19768 and 19769) represent a partial codon and are discarded in computing codon frequencies.</p>
<title>Results and discussion</title> <title>Wobble cost between G and U: MU3:G and MG3:U</title> <p>Recall that the two inequalities</p> <p></p> <p>are, respectively, for NNY codons translated by tRNA with a G at the wobble site of tRNA anticodon, and for NNR codons translated by tRNA with a U at the wobble site of tRNA anticodon. The observed N<sub>C</sub>/N<sub>U </sub>ratios for <italic>Allomyces macrogynus </italic>(representing fungal mitochondrial genomes with translation table 4) are much smaller than N<sub>A</sub>/N<sub>G </sub>ratios (Table ##TAB##1##2##). The smallest N<sub>C</sub>/N<sub>U </sub>value is 0.279 whereas the smallest N<sub>A</sub>/N<sub>G </sub>value is 2.372 (Table ##TAB##1##2##). We have mentioned before that, if M<sub>U3:G </sub>is very small, then a wobble G at the tRNA anticodon will not impose strong selection against U-ending codons, and N<sub>U </sub>may drift up and down with mutation relative to N<sub>C</sub>. This will lead to relatively small N<sub>C</sub>/N<sub>U </sub>ratios. From the minimum N<sub>C</sub>/N<sub>U </sub>value of 0.279, we may infer that M<sub>U3:G </sub>&lt; 0.279•M<sub>O</sub>, i.e., M<sub>U3:G </sub>is quite small relative to M<sub>O</sub>.</p> <p>The AUH codon family coding for amino acid Ile in <italic>A. macrogynus </italic>mitochondrial genome is translated by a tRNA with a GAU anticodon. According to Eq. (7), the M<sub>U3:G</sub>/Co ratio should also be smaller than the N<sub>C</sub>/N<sub>U </sub>ratio. The observed N<sub>C</sub>/N<sub>U </sub>ratio is 0.3605 (= 159/441). This is similar to the N<sub>C</sub>/N<sub>U </sub>ratio in NNY codon families (Table ##TAB##1##2##). Thus, the wobble cost of M<sub>U3:G </sub>relative to M<sub>O </sub>from the AUH codon family is similar to that derived from NNY codon families.</p> <p>The CUN codon family coding for amino acid Leu in <italic>A. macrogynus </italic>mitochondrial genome is translated by a tRNA with an AAG anticodon. Note that no A→I conversion has been observed in mitochondria [##REF##17878308##37##,##REF##7541252##38##] so Eq. (9) is applicable. According to Eq. (9), the M<sub>U3:G</sub>/Co ratio should be greater than the N<sub>C</sub>/N<sub>U </sub>ratio. The observed N<sub>C</sub>/N<sub>U </sub>ratio is 0.1905 (= 48/252). Thus we have M<sub>U3:G </sub>&gt; 0.1905•M<sub>O</sub>. This inequality, together with the previous inequality of M<sub>U3:G </sub>&lt; 0.279•M<sub>O</sub>, leads to 0.1905•M<sub>O </sub>&lt; M<sub>U3:G </sub>&lt; 0.279•M<sub>O</sub>.</p> <p>Given that M<sub>U3:G </sub>&gt; 0.1905•M<sub>O</sub>, we can infer that M<sub>wG </sub>&gt; M<sub>wA </sub>when N<sub>C</sub>/N<sub>U </sub>&lt; 0.1905, i.e., when U-ending codons are more than five times as frequent as C-ending codons. In other words, when N<sub>C</sub>/N<sub>U </sub>&lt; 0.1905, the wobble cost of having an A at the wobble site of the tRNA anticodon is smaller than that of having a G at the wobble site, and should be favored by natural selection. Among the 36 fungal mitochondrial genomes, only two genomes have a NNY codon family with a correspondent tRNA that has an A at the wobble site of its anticodon. In the mitochondrial genome of <italic>Penicillium marneffei</italic>, the tRNA translating the AAY codon family has a wobble A in its anticodon. The N<sub>C</sub>/N<sub>U </sub>ratio is 0.087 (= 41/473) which is much smaller than 0.1905. Similarly, in the mitochondrial genome of <italic>Pichia Canadensis</italic>, the tRNA translating the AGY codon family has a wobble A in its anticodon. The N<sub>C</sub>/N<sub>U </sub>ratio is 0.0083 (= 1/120) which is also much smaller than 0.1905. Thus, the prediction that a wobble A at the tRNA anticodon is advantageous over a wobble G when N<sub>C</sub>/N<sub>U </sub>&lt; 0.1905 is consistent with empirical data.</p> <p>In contrast to the small N<sub>C</sub>/N<sub>U </sub>ratios in NNY, AUH and CUN codon families, all N<sub>A</sub>/N<sub>G </sub>ratios in NNR codon families are substantially larger (Table ##TAB##1##2##). We have argued before that, if M<sub>G3:U </sub>is very small, then a U at the wobble site of tRNA anticodon would impose little selection against G-ending codons in NNR codon families, and mutation may allow N<sub>G </sub>to drift up, leading to large N<sub>G </sub>values relative to N<sub>A</sub>. However, if M<sub>G3:U </sub>is large, then G-ending codons should be strongly selected against and N<sub>G </sub>would be small relative to N<sub>A</sub>, leading to large N<sub>A</sub>/N<sub>G </sub>ratios. The much larger N<sub>A</sub>/N<sub>G </sub>ratios than N<sub>C</sub>/N<sub>U </sub>ratios (t = 5.2967, DF = 10, p = 0.0003, two-tailed test) suggest that M<sub>G3:U </sub>is much greater than M<sub>U3:G</sub>.</p> <p>There is a caveat in evaluating the relative magnitude of M<sub>U3:G</sub>/M<sub>O</sub>, and M<sub>G3:U</sub>/M<sub>O </sub>by the N<sub>C</sub>/N<sub>U </sub>and N<sub>A</sub>/N<sub>G </sub>ratios because these ratios can be affected by AT-biased mutations. The mitochondrial genome of <italic>A. macrogynus </italic>is 57472 bp, with the number of C+G being 22700 and that of A+U being 34772. If we exclude those nucleotides in coding sequences, then the numbers of C+G and A+U are 14136 and 17656, respectively. This may be considered as the background frequencies maintained by mutation bias, which leads to the expected N<sub>C</sub>/N<sub>U </sub>ratio of 0.8001 (=N<sub>C+G</sub>/N<sub>A+U</sub>) and that of N<sub>A</sub>/N<sub>G </sub>ratio of 1.2490 (=N<sub>A+U</sub>/N<sub>C+G</sub>). Thus, we expect N<sub>C</sub>/N<sub>U </sub>to be smaller than N<sub>A</sub>/N<sub>G </sub>even when there is no difference between M<sub>U3:G </sub>and M<sub>G3:U</sub>. To establish the argument that M<sub>U3:G </sub>is indeed smaller than M<sub>G3:U</sub>, we need to demonstrate that (1) there is no selection against U-ending codons in NNY codon families by showing that the observed N<sub>C</sub>/N<sub>U </sub>ratio is not greater than 0.8001, and (2) there is selection against G-ending codons in NNR codon families by showing that N<sub>A</sub>/N<sub>G </sub>is significantly greater than 1.2490. It is not enough to show that N<sub>C</sub>/N<sub>U </sub>&lt; N<sub>A</sub>/N<sub>G</sub>.</p> <p>We note that the observed N<sub>C</sub>/N<sub>U </sub>values for the seven NNY codon families in Table ##TAB##1##2## are all smaller than the expected 0.8001, suggesting no selection against U-ending codons in NNY codon families (i.e., small M<sub>U3:G</sub>). In contrast the observed N<sub>A</sub>/N<sub>G </sub>ratios for the five NNR codon families are all much greater than the expected 1.2490 (Table ##TAB##1##2##), consistent with the interpretation of selection against G-ending codons in NNR codon families (i.e., large M<sub>G3:U</sub>). This is consistent with the interpretation that M<sub>U3:G </sub>&lt; M<sub>G3:U</sub>. One can perform a χ<sup>2</sup>-test for each of the NNR codon families to see if G-ending codons are underused. The tests are all highly significant, with p &lt; 0.00001.</p> <p>The results are similar for fungal genomes with translation table 3, with the result from a representative species (<italic>Saccharomyces cerevisiae</italic>) presented in Table ##TAB##2##3##. Again the N<sub>C</sub>/N<sub>U </sub>ratios in NNY codon families are much smaller than N<sub>A</sub>/N<sub>G </sub>ratios in NNR codon families. We should note that the <italic>S. cerevisiae </italic>mitochondrial genome is much more AT-biased than the <italic>A. macrogynus </italic>mitochondrial genome, with the proportion of (G+C) in non-coding sequences being only 0.1484. The reason for the GC deficiency in yeast is not clear, but it may be caused either by mutation bias or by the low abundance of C in living cells [##REF##12044357##39##, ####REF##16475092##40##, ##REF##16687416##41####16687416##41##]. In any case, the expected N<sub>C</sub>/N<sub>U </sub>and N<sub>A</sub>/N<sub>G </sub>ratios, given the biased genomic AT content, are 0.1742 and 5.7405, respectively. We note that the observed N<sub>C</sub>/N<sub>U </sub>ratios among the NNY codon families are all smaller than the expected value of 0.1742 except for the UUY (Phe) codon family (Table ##TAB##2##3##), against suggesting little selection against U-ending codons (i.e., small M<sub>U3:G</sub>). In contrast, the N<sub>A</sub>/N<sub>G </sub>ratios are all significantly greater than the expected 5.7405 except for the AUR (Met) codon family, suggesting selection against G-ending codons (i.e., large M<sub>G3:U</sub>). The exceptional AUR (Met) codon family has a tRNA with a CAU anticodon which would favor G-ending codons and is expected to be different.</p> <p>The N<sub>C</sub>/N<sub>U </sub>ratio for the UUY (Phe) codon family is consistently greater than those of other NNY codon families (Tables ##TAB##1##2##, ##TAB##2##3##). One may suspect whether, for this particular codon family, there is a significant M<sub>U3:G</sub>. The rate of tRNA<sup>Phe </sup>with anticodon 3'-AAG-5' dissociating from the UUU codon is about twice as high as that from the fully matched UUC codon [##REF##16455492##42##]. Also, the tRNA<sup>Phe </sup>misread CUC codons more than twice as often as CUU codons [##REF##16455492##42##]. These two lines of evidence suggest that C3:G pair is much more favorable than U3:G pair, i.e., M<sub>U3:G </sub>for the UUY (Phe) codon family may indeed be substantially greater than M<sub>C3:G</sub>. Unfortunately, there has been no other similar studies on codon-anticodon pairing for other NNY codon families. One should also note that the tRNA<sup>Phe </sup>in this study comes from <italic>Escherichia coli</italic>, and the result may not be applicable to fungal species.</p> <p>The two fungal species using genetic table 16, i.e., <italic>Spizellomyces punctatus </italic>and <italic>Rhizophydium sp. 136 </italic>each have only a partial set of tRNA genes. Among NNY codon families in <italic>S. punctatus</italic>, only the GAY (coding for Asp) and UAY codon family (coding for Tyr) have an identifiable tRNA, with anticodons being GUC and GUA, respectively. The expected N<sub>C</sub>/N<sub>U </sub>and N<sub>A</sub>/N<sub>G </sub>ratios are 0.4536 and 2.2048, respectively, based on the nucleotide frequencies of non-coding sequences in the mitochondrial genome. The observed N<sub>C</sub>/N<sub>U </sub>ratios for GAY and UAY codon families are 0.1667 (= 40/240) and 0.2411 (= 81/336), respectively, suggesting no selection against U-ending codons (i.e., small M<sub>U3:G</sub>). In contrast, the N<sub>A</sub>/N<sub>G </sub>ratios are much larger, being 19.2308 (= 250/13) for the AAR codon family (coding for amino acid Lys) translated by a tRNA with a UUU anticodon and 11.1538 (= 145/13) for the CAR codon family (coding for amino acid Gln) translated by a tRNA with a UUG anticodon. These results suggest selection against G-ending codons (i.e., large M<sub>G3:U</sub>) There is no other NNR codon family with identifiable tRNAs in <italic>S. punctatus</italic>. The other species using translation table 16, <italic>Rhizophydium sp. 136</italic>, exhibits a similar pattern.</p> <p>In short, results from fungal mitochondrial genomes are consistent with no selection against U-ending codons in NNY codon families but significant selection against G-ending codons in NNR codon families, indicating that M<sub>U3:G </sub>is smaller than M<sub>G3:U</sub>. These findings corroborate previous biochemical studies demonstrating that U3:G is energetically much more favorable than G3:U [##REF##14715921##2##,##REF##15952884##23##, ####REF##724522##24##, ##REF##17496902##25##, ##REF##10668794##26####10668794##26##]. However, G3:U can be almost as good as A3:U when U is modified to xm<sup>5</sup>U [##REF##18456657##43##].</p> <p>The finding of a small M<sub>U3:G </sub>can explain puzzling observations in codon usage in fungal mitochondrial genomes. Take the tRNA<sup>Ser </sup>translating the AGY codon family in the mitochondrial genome of <italic>Ashbya gossypii </italic>ATCC 10895 for example. The genome contains 31 AGU codons and no AGC codon. CAAH would have predicted an ACU anticodon with perfect base pair with AGU codons, but the observed anticodon is GCU consistent with WVH. Such an anticodon makes sense only if M<sub>U3:G </sub>is equal to M<sub>C3:G</sub>, i.e., a wobble G in a tRNA anticodon can pair with U just as good as with C in the third codon position.</p> <p>Can this finding be generalized to nuclear genomes and the translation machinery in the cytoplasm? If M<sub>U3:G </sub>is generally small, then we have an answer to a puzzling observation that has long baffled molecular evolutionary geneticists. The UGY codon family (coding for Cys) in the nuclear genome of <italic>S. cerevisiae </italic>is translated by tRNA coded in four tRNA<sup>Cys </sup>genes all with a wobble G at the tRNA anticodon. One would have predicted that UGC codons should be strongly preferred over UGU codons to avoid the wobble cost M<sub>U3:G</sub>. However, the observed numbers of UGC and UGU codons in <italic>S. cerevisiae </italic>CDSs is 13802 and 22873, respectively, opposite to the predicted trend. The bias is even stronger in highly expressed genes. For example, in the codon usage table of highly expressed genes distributed with the EMBOSS package [[##REF##10827456##44##], Eyeastcai.cut], the numbers of UGC and UGU codons are 3 and 39, respectively. The unexpected codon usage bias (unexpected because the bias is in the wrong direction) becomes easy to understand if M<sub>U3:G </sub>is equal to M<sub>C3:G</sub>, i.e., there is no need to overuse the UGC codons to avoid M<sub>U3:G </sub>so AT-biased mutation in yeast (whose genome is highly AT-biased), which increases the frequency of U-ending codons, is not checked by counteracting selection.</p> <p>If the finding of M<sub>U3:G </sub>&lt;&lt; M<sub>G3:U </sub>is applicable to nuclear genomes, then we predict that NNY codon families need only one type of tRNA (i.e., tRNA with a wobble G at the tRNA anticodon) for translation. In contrast, NNR codon families should ideally be translated by two different types of tRNAs, one with a wobble U for NNA codons and the other with a wobble C for NNG codons (to avoid the relatively high wobble cost of M<sub>G3:U</sub>). A corollary is that, if a NNY codon family is translated by only tRNA with a wobble G and NNR codon family by tRNA with a wobble U, then codon usage bias should be smaller in the NNY codon family (in which selection against G-ending codon is weak because of the small M<sub>U3:G</sub>) than in the NNR codon family (in which selection against U-ending codon is strong because of the relatively large M<sub>G3:U</sub>). Below we test these predictions with genomic data.</p> <title>NNY codon families are translated by one type of tRNA with a wobble G and NNR codon families are translated by two types of tRNA with a wobble U and a wobble C, respectively, in fungal nuclear genomes</title> <p>We have inferred that a tRNA with a wobble G at its anticodon should be efficient not only in translating C-ending codons, but also in translating U-ending codons because of the small M<sub>U3:G</sub>. For this reason, only one type of tRNA with a wobble G should generally be sufficient in translating NNY codon families. In contrast, in NNR codon families, tRNA with a wobble U at its anticodon should be poor in translating G-ending codons because of the large M<sub>G3:U</sub>. Thus, the presence of G-ending codons in NNR codon families should favour the use of two types of tRNAs, one with a wobble U for translating NNA codons and another with a wobble C for translating NNG codons. In mitochondrial genomes with limited gene content, each codon family is generally translated by a single tRNA species. So this prediction cannot be tested. However, this prediction can be tested with nuclear genomes where the number of tRNA genes is not so limited as in mitochondrial genomes.</p> <p>The prediction is strongly supported by results from the nuclear genome of <italic>Saccharomyces cerevisiae </italic>(Table ##TAB##3##4##). All NNR codon families are translated by two types of tRNAs with anticodons matching NNA and NNG codons, respectively, whereas all NNY codon families are translated by one type of tRNA with a wobble G at the tRNA anticodon (Table ##TAB##3##4##). This is consistent with the interpretation that the inference of M<sub>U3:G </sub>&lt;&lt; M<sub>G3:U </sub>derived from fungal mitochondrial genome is also applicable to fungal nuclear genomes.</p> <p>One may propose an alternative hypothesis for the observation that NNY codon families are translated by tRNAs with a wobble G whereas NNR codon families translated by tRNAs with a wobble C and a wobble U, respectively. The wobble A in some nuclear tRNAs is known to be converted to inosine or I for short [##REF##6369251##45##, ####REF##8915538##46##, ##REF##8893855##47####8893855##47##] which may pair with A, C or U. If a tRNA translating an NNY codon family has its wobble G mutated to wobble A, then the wobble A may undergo the A→I conversion and misread NNA codons. For this reason, the wobble G→A mutation should be strongly selected against, which would explain the lack of wobble A in tRNA translating NNY codons.</p> <p>This alternative hypothesis invoking the A→I conversion, while logically sound, was dismissed quite early after the discovery [##REF##6369251##45##,##REF##8893855##47##] that the A→I conversion is quite restrictive and occurs mainly at ACG or ARV anticodons (where R is the IUB code for A or G, and V is for A, G or C). Among tRNAs translating NNY codons, only tRNA<sup>Phe </sup>has a middle A in the anticodon, and the rest do not have an R in the middle of the anticodon. This means that, for all NNY-translating tRNAs except for tRNA<sup>Phe</sup>, even if their regular wobble G mutates to A, the resulting wobble A will NOT be converted to inosine. Thus, the alternative hypothesis can only explain the avoidance of a wobble A in tRNA<sup>Phe </sup>but cannot explain the avoidance of a wobble A in other NNY-translating tRNAs.</p>
<title>Results and discussion</title> <title>Wobble cost between G and U: MU3:G and MG3:U</title> <p>Recall that the two inequalities</p> <p></p> <p>are, respectively, for NNY codons translated by tRNA with a G at the wobble site of tRNA anticodon, and for NNR codons translated by tRNA with a U at the wobble site of tRNA anticodon. The observed N<sub>C</sub>/N<sub>U </sub>ratios for <italic>Allomyces macrogynus </italic>(representing fungal mitochondrial genomes with translation table 4) are much smaller than N<sub>A</sub>/N<sub>G </sub>ratios (Table ##TAB##1##2##). The smallest N<sub>C</sub>/N<sub>U </sub>value is 0.279 whereas the smallest N<sub>A</sub>/N<sub>G </sub>value is 2.372 (Table ##TAB##1##2##). We have mentioned before that, if M<sub>U3:G </sub>is very small, then a wobble G at the tRNA anticodon will not impose strong selection against U-ending codons, and N<sub>U </sub>may drift up and down with mutation relative to N<sub>C</sub>. This will lead to relatively small N<sub>C</sub>/N<sub>U </sub>ratios. From the minimum N<sub>C</sub>/N<sub>U </sub>value of 0.279, we may infer that M<sub>U3:G </sub>&lt; 0.279•M<sub>O</sub>, i.e., M<sub>U3:G </sub>is quite small relative to M<sub>O</sub>.</p> <p>The AUH codon family coding for amino acid Ile in <italic>A. macrogynus </italic>mitochondrial genome is translated by a tRNA with a GAU anticodon. According to Eq. (7), the M<sub>U3:G</sub>/Co ratio should also be smaller than the N<sub>C</sub>/N<sub>U </sub>ratio. The observed N<sub>C</sub>/N<sub>U </sub>ratio is 0.3605 (= 159/441). This is similar to the N<sub>C</sub>/N<sub>U </sub>ratio in NNY codon families (Table ##TAB##1##2##). Thus, the wobble cost of M<sub>U3:G </sub>relative to M<sub>O </sub>from the AUH codon family is similar to that derived from NNY codon families.</p> <p>The CUN codon family coding for amino acid Leu in <italic>A. macrogynus </italic>mitochondrial genome is translated by a tRNA with an AAG anticodon. Note that no A→I conversion has been observed in mitochondria [##REF##17878308##37##,##REF##7541252##38##] so Eq. (9) is applicable. According to Eq. (9), the M<sub>U3:G</sub>/Co ratio should be greater than the N<sub>C</sub>/N<sub>U </sub>ratio. The observed N<sub>C</sub>/N<sub>U </sub>ratio is 0.1905 (= 48/252). Thus we have M<sub>U3:G </sub>&gt; 0.1905•M<sub>O</sub>. This inequality, together with the previous inequality of M<sub>U3:G </sub>&lt; 0.279•M<sub>O</sub>, leads to 0.1905•M<sub>O </sub>&lt; M<sub>U3:G </sub>&lt; 0.279•M<sub>O</sub>.</p> <p>Given that M<sub>U3:G </sub>&gt; 0.1905•M<sub>O</sub>, we can infer that M<sub>wG </sub>&gt; M<sub>wA </sub>when N<sub>C</sub>/N<sub>U </sub>&lt; 0.1905, i.e., when U-ending codons are more than five times as frequent as C-ending codons. In other words, when N<sub>C</sub>/N<sub>U </sub>&lt; 0.1905, the wobble cost of having an A at the wobble site of the tRNA anticodon is smaller than that of having a G at the wobble site, and should be favored by natural selection. Among the 36 fungal mitochondrial genomes, only two genomes have a NNY codon family with a correspondent tRNA that has an A at the wobble site of its anticodon. In the mitochondrial genome of <italic>Penicillium marneffei</italic>, the tRNA translating the AAY codon family has a wobble A in its anticodon. The N<sub>C</sub>/N<sub>U </sub>ratio is 0.087 (= 41/473) which is much smaller than 0.1905. Similarly, in the mitochondrial genome of <italic>Pichia Canadensis</italic>, the tRNA translating the AGY codon family has a wobble A in its anticodon. The N<sub>C</sub>/N<sub>U </sub>ratio is 0.0083 (= 1/120) which is also much smaller than 0.1905. Thus, the prediction that a wobble A at the tRNA anticodon is advantageous over a wobble G when N<sub>C</sub>/N<sub>U </sub>&lt; 0.1905 is consistent with empirical data.</p> <p>In contrast to the small N<sub>C</sub>/N<sub>U </sub>ratios in NNY, AUH and CUN codon families, all N<sub>A</sub>/N<sub>G </sub>ratios in NNR codon families are substantially larger (Table ##TAB##1##2##). We have argued before that, if M<sub>G3:U </sub>is very small, then a U at the wobble site of tRNA anticodon would impose little selection against G-ending codons in NNR codon families, and mutation may allow N<sub>G </sub>to drift up, leading to large N<sub>G </sub>values relative to N<sub>A</sub>. However, if M<sub>G3:U </sub>is large, then G-ending codons should be strongly selected against and N<sub>G </sub>would be small relative to N<sub>A</sub>, leading to large N<sub>A</sub>/N<sub>G </sub>ratios. The much larger N<sub>A</sub>/N<sub>G </sub>ratios than N<sub>C</sub>/N<sub>U </sub>ratios (t = 5.2967, DF = 10, p = 0.0003, two-tailed test) suggest that M<sub>G3:U </sub>is much greater than M<sub>U3:G</sub>.</p> <p>There is a caveat in evaluating the relative magnitude of M<sub>U3:G</sub>/M<sub>O</sub>, and M<sub>G3:U</sub>/M<sub>O </sub>by the N<sub>C</sub>/N<sub>U </sub>and N<sub>A</sub>/N<sub>G </sub>ratios because these ratios can be affected by AT-biased mutations. The mitochondrial genome of <italic>A. macrogynus </italic>is 57472 bp, with the number of C+G being 22700 and that of A+U being 34772. If we exclude those nucleotides in coding sequences, then the numbers of C+G and A+U are 14136 and 17656, respectively. This may be considered as the background frequencies maintained by mutation bias, which leads to the expected N<sub>C</sub>/N<sub>U </sub>ratio of 0.8001 (=N<sub>C+G</sub>/N<sub>A+U</sub>) and that of N<sub>A</sub>/N<sub>G </sub>ratio of 1.2490 (=N<sub>A+U</sub>/N<sub>C+G</sub>). Thus, we expect N<sub>C</sub>/N<sub>U </sub>to be smaller than N<sub>A</sub>/N<sub>G </sub>even when there is no difference between M<sub>U3:G </sub>and M<sub>G3:U</sub>. To establish the argument that M<sub>U3:G </sub>is indeed smaller than M<sub>G3:U</sub>, we need to demonstrate that (1) there is no selection against U-ending codons in NNY codon families by showing that the observed N<sub>C</sub>/N<sub>U </sub>ratio is not greater than 0.8001, and (2) there is selection against G-ending codons in NNR codon families by showing that N<sub>A</sub>/N<sub>G </sub>is significantly greater than 1.2490. It is not enough to show that N<sub>C</sub>/N<sub>U </sub>&lt; N<sub>A</sub>/N<sub>G</sub>.</p> <p>We note that the observed N<sub>C</sub>/N<sub>U </sub>values for the seven NNY codon families in Table ##TAB##1##2## are all smaller than the expected 0.8001, suggesting no selection against U-ending codons in NNY codon families (i.e., small M<sub>U3:G</sub>). In contrast the observed N<sub>A</sub>/N<sub>G </sub>ratios for the five NNR codon families are all much greater than the expected 1.2490 (Table ##TAB##1##2##), consistent with the interpretation of selection against G-ending codons in NNR codon families (i.e., large M<sub>G3:U</sub>). This is consistent with the interpretation that M<sub>U3:G </sub>&lt; M<sub>G3:U</sub>. One can perform a χ<sup>2</sup>-test for each of the NNR codon families to see if G-ending codons are underused. The tests are all highly significant, with p &lt; 0.00001.</p> <p>The results are similar for fungal genomes with translation table 3, with the result from a representative species (<italic>Saccharomyces cerevisiae</italic>) presented in Table ##TAB##2##3##. Again the N<sub>C</sub>/N<sub>U </sub>ratios in NNY codon families are much smaller than N<sub>A</sub>/N<sub>G </sub>ratios in NNR codon families. We should note that the <italic>S. cerevisiae </italic>mitochondrial genome is much more AT-biased than the <italic>A. macrogynus </italic>mitochondrial genome, with the proportion of (G+C) in non-coding sequences being only 0.1484. The reason for the GC deficiency in yeast is not clear, but it may be caused either by mutation bias or by the low abundance of C in living cells [##REF##12044357##39##, ####REF##16475092##40##, ##REF##16687416##41####16687416##41##]. In any case, the expected N<sub>C</sub>/N<sub>U </sub>and N<sub>A</sub>/N<sub>G </sub>ratios, given the biased genomic AT content, are 0.1742 and 5.7405, respectively. We note that the observed N<sub>C</sub>/N<sub>U </sub>ratios among the NNY codon families are all smaller than the expected value of 0.1742 except for the UUY (Phe) codon family (Table ##TAB##2##3##), against suggesting little selection against U-ending codons (i.e., small M<sub>U3:G</sub>). In contrast, the N<sub>A</sub>/N<sub>G </sub>ratios are all significantly greater than the expected 5.7405 except for the AUR (Met) codon family, suggesting selection against G-ending codons (i.e., large M<sub>G3:U</sub>). The exceptional AUR (Met) codon family has a tRNA with a CAU anticodon which would favor G-ending codons and is expected to be different.</p> <p>The N<sub>C</sub>/N<sub>U </sub>ratio for the UUY (Phe) codon family is consistently greater than those of other NNY codon families (Tables ##TAB##1##2##, ##TAB##2##3##). One may suspect whether, for this particular codon family, there is a significant M<sub>U3:G</sub>. The rate of tRNA<sup>Phe </sup>with anticodon 3'-AAG-5' dissociating from the UUU codon is about twice as high as that from the fully matched UUC codon [##REF##16455492##42##]. Also, the tRNA<sup>Phe </sup>misread CUC codons more than twice as often as CUU codons [##REF##16455492##42##]. These two lines of evidence suggest that C3:G pair is much more favorable than U3:G pair, i.e., M<sub>U3:G </sub>for the UUY (Phe) codon family may indeed be substantially greater than M<sub>C3:G</sub>. Unfortunately, there has been no other similar studies on codon-anticodon pairing for other NNY codon families. One should also note that the tRNA<sup>Phe </sup>in this study comes from <italic>Escherichia coli</italic>, and the result may not be applicable to fungal species.</p> <p>The two fungal species using genetic table 16, i.e., <italic>Spizellomyces punctatus </italic>and <italic>Rhizophydium sp. 136 </italic>each have only a partial set of tRNA genes. Among NNY codon families in <italic>S. punctatus</italic>, only the GAY (coding for Asp) and UAY codon family (coding for Tyr) have an identifiable tRNA, with anticodons being GUC and GUA, respectively. The expected N<sub>C</sub>/N<sub>U </sub>and N<sub>A</sub>/N<sub>G </sub>ratios are 0.4536 and 2.2048, respectively, based on the nucleotide frequencies of non-coding sequences in the mitochondrial genome. The observed N<sub>C</sub>/N<sub>U </sub>ratios for GAY and UAY codon families are 0.1667 (= 40/240) and 0.2411 (= 81/336), respectively, suggesting no selection against U-ending codons (i.e., small M<sub>U3:G</sub>). In contrast, the N<sub>A</sub>/N<sub>G </sub>ratios are much larger, being 19.2308 (= 250/13) for the AAR codon family (coding for amino acid Lys) translated by a tRNA with a UUU anticodon and 11.1538 (= 145/13) for the CAR codon family (coding for amino acid Gln) translated by a tRNA with a UUG anticodon. These results suggest selection against G-ending codons (i.e., large M<sub>G3:U</sub>) There is no other NNR codon family with identifiable tRNAs in <italic>S. punctatus</italic>. The other species using translation table 16, <italic>Rhizophydium sp. 136</italic>, exhibits a similar pattern.</p> <p>In short, results from fungal mitochondrial genomes are consistent with no selection against U-ending codons in NNY codon families but significant selection against G-ending codons in NNR codon families, indicating that M<sub>U3:G </sub>is smaller than M<sub>G3:U</sub>. These findings corroborate previous biochemical studies demonstrating that U3:G is energetically much more favorable than G3:U [##REF##14715921##2##,##REF##15952884##23##, ####REF##724522##24##, ##REF##17496902##25##, ##REF##10668794##26####10668794##26##]. However, G3:U can be almost as good as A3:U when U is modified to xm<sup>5</sup>U [##REF##18456657##43##].</p> <p>The finding of a small M<sub>U3:G </sub>can explain puzzling observations in codon usage in fungal mitochondrial genomes. Take the tRNA<sup>Ser </sup>translating the AGY codon family in the mitochondrial genome of <italic>Ashbya gossypii </italic>ATCC 10895 for example. The genome contains 31 AGU codons and no AGC codon. CAAH would have predicted an ACU anticodon with perfect base pair with AGU codons, but the observed anticodon is GCU consistent with WVH. Such an anticodon makes sense only if M<sub>U3:G </sub>is equal to M<sub>C3:G</sub>, i.e., a wobble G in a tRNA anticodon can pair with U just as good as with C in the third codon position.</p> <p>Can this finding be generalized to nuclear genomes and the translation machinery in the cytoplasm? If M<sub>U3:G </sub>is generally small, then we have an answer to a puzzling observation that has long baffled molecular evolutionary geneticists. The UGY codon family (coding for Cys) in the nuclear genome of <italic>S. cerevisiae </italic>is translated by tRNA coded in four tRNA<sup>Cys </sup>genes all with a wobble G at the tRNA anticodon. One would have predicted that UGC codons should be strongly preferred over UGU codons to avoid the wobble cost M<sub>U3:G</sub>. However, the observed numbers of UGC and UGU codons in <italic>S. cerevisiae </italic>CDSs is 13802 and 22873, respectively, opposite to the predicted trend. The bias is even stronger in highly expressed genes. For example, in the codon usage table of highly expressed genes distributed with the EMBOSS package [[##REF##10827456##44##], Eyeastcai.cut], the numbers of UGC and UGU codons are 3 and 39, respectively. The unexpected codon usage bias (unexpected because the bias is in the wrong direction) becomes easy to understand if M<sub>U3:G </sub>is equal to M<sub>C3:G</sub>, i.e., there is no need to overuse the UGC codons to avoid M<sub>U3:G </sub>so AT-biased mutation in yeast (whose genome is highly AT-biased), which increases the frequency of U-ending codons, is not checked by counteracting selection.</p> <p>If the finding of M<sub>U3:G </sub>&lt;&lt; M<sub>G3:U </sub>is applicable to nuclear genomes, then we predict that NNY codon families need only one type of tRNA (i.e., tRNA with a wobble G at the tRNA anticodon) for translation. In contrast, NNR codon families should ideally be translated by two different types of tRNAs, one with a wobble U for NNA codons and the other with a wobble C for NNG codons (to avoid the relatively high wobble cost of M<sub>G3:U</sub>). A corollary is that, if a NNY codon family is translated by only tRNA with a wobble G and NNR codon family by tRNA with a wobble U, then codon usage bias should be smaller in the NNY codon family (in which selection against G-ending codon is weak because of the small M<sub>U3:G</sub>) than in the NNR codon family (in which selection against U-ending codon is strong because of the relatively large M<sub>G3:U</sub>). Below we test these predictions with genomic data.</p> <title>NNY codon families are translated by one type of tRNA with a wobble G and NNR codon families are translated by two types of tRNA with a wobble U and a wobble C, respectively, in fungal nuclear genomes</title> <p>We have inferred that a tRNA with a wobble G at its anticodon should be efficient not only in translating C-ending codons, but also in translating U-ending codons because of the small M<sub>U3:G</sub>. For this reason, only one type of tRNA with a wobble G should generally be sufficient in translating NNY codon families. In contrast, in NNR codon families, tRNA with a wobble U at its anticodon should be poor in translating G-ending codons because of the large M<sub>G3:U</sub>. Thus, the presence of G-ending codons in NNR codon families should favour the use of two types of tRNAs, one with a wobble U for translating NNA codons and another with a wobble C for translating NNG codons. In mitochondrial genomes with limited gene content, each codon family is generally translated by a single tRNA species. So this prediction cannot be tested. However, this prediction can be tested with nuclear genomes where the number of tRNA genes is not so limited as in mitochondrial genomes.</p> <p>The prediction is strongly supported by results from the nuclear genome of <italic>Saccharomyces cerevisiae </italic>(Table ##TAB##3##4##). All NNR codon families are translated by two types of tRNAs with anticodons matching NNA and NNG codons, respectively, whereas all NNY codon families are translated by one type of tRNA with a wobble G at the tRNA anticodon (Table ##TAB##3##4##). This is consistent with the interpretation that the inference of M<sub>U3:G </sub>&lt;&lt; M<sub>G3:U </sub>derived from fungal mitochondrial genome is also applicable to fungal nuclear genomes.</p> <p>One may propose an alternative hypothesis for the observation that NNY codon families are translated by tRNAs with a wobble G whereas NNR codon families translated by tRNAs with a wobble C and a wobble U, respectively. The wobble A in some nuclear tRNAs is known to be converted to inosine or I for short [##REF##6369251##45##, ####REF##8915538##46##, ##REF##8893855##47####8893855##47##] which may pair with A, C or U. If a tRNA translating an NNY codon family has its wobble G mutated to wobble A, then the wobble A may undergo the A→I conversion and misread NNA codons. For this reason, the wobble G→A mutation should be strongly selected against, which would explain the lack of wobble A in tRNA translating NNY codons.</p> <p>This alternative hypothesis invoking the A→I conversion, while logically sound, was dismissed quite early after the discovery [##REF##6369251##45##,##REF##8893855##47##] that the A→I conversion is quite restrictive and occurs mainly at ACG or ARV anticodons (where R is the IUB code for A or G, and V is for A, G or C). Among tRNAs translating NNY codons, only tRNA<sup>Phe </sup>has a middle A in the anticodon, and the rest do not have an R in the middle of the anticodon. This means that, for all NNY-translating tRNAs except for tRNA<sup>Phe</sup>, even if their regular wobble G mutates to A, the resulting wobble A will NOT be converted to inosine. Thus, the alternative hypothesis can only explain the avoidance of a wobble A in tRNA<sup>Phe </sup>but cannot explain the avoidance of a wobble A in other NNY-translating tRNAs.</p>
<title>Conclusion</title> <p>In summary, our general hypothesis based on wobble costs allows the integration of the two conventional hypotheses (i.e., CAAH and WVH) on codon-anticodon wobble pairing. The integration leads to new ways of evaluating relative wobble cost of different wobble pairings. In particular, the finding that M<sub>U3:G </sub>is much smaller than M<sub>G3:U </sub>corroborates previous structural studies showing that U3:G is energetically more favourable than G3:U and leads to a better understanding of the translation efficiency mediated by codon and anticodon wobble pairing.</p>
<title>Background</title> <p>Fungal and animal mitochondrial genomes typically have one tRNA for each synonymous codon family. The codon-anticodon adaptation hypothesis predicts that the wobble nucleotide of a tRNA anticodon should evolve towards maximizing Watson-Crick base pairing with the most frequently used codon within each synonymous codon family, whereas 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 tRNA wobble nucleotide 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).</p> <title>Results</title> <p>We here integrate these two traditional hypotheses on tRNA anticodons into a unified model based on an analysis of the wobble costs associated with different wobble base pairs. This novel approach allows the relative cost of wobble pairing to be qualitatively evaluated. A comprehensive study of 36 fungal genomes suggests very different costs between two kinds of U:G wobble pairs, i.e., (1) between a G at the wobble site of a tRNA anticodon and a U at the third codon position (designated M<sub>U3:G</sub>) and (2) between a U at the wobble site of a tRNA anticodon and a G at the third codon position (designated M<sub>G3:U</sub>).</p> <title>Conclusion</title> <p>In general, M<sub>U3:G </sub>is much smaller than M<sub>G3:U</sub>, suggesting no selection against U-ending codons in NNY codon families with a wobble G in the tRNA anticodon but strong selection against G-ending codons in NNR codon families with a wobble U at the tRNA anticodon. This finding resolves several puzzling observations in fungal genomics and corroborates previous studies showing that U3:G wobble is energetically more favorable than G3:U wobble.</p>
<title>Acknowledgements</title> <p>This study is supported by the CAS/SAFEA International Partnership Program for Creative Research Teams and by NSERC's Discovery and Strategic Grants. I thank S. Aris-Brosou, P. Higgs and E. Rocha for discussion and comments and E. Westhof for references. Four anonymous reviewers provided helpful comments leading to substantial improvement of the paper.</p>
<fig id="F1" position="float"><label>Figure 1</label><caption><p>Conceptual illustration of the dependence of wobbling cost involving a G or an A at the wobble site of tRNA anticodon (M<sub>wG </sub>and M<sub>wA</sub>) on the proportion of C-ending codons (P<sub>C</sub>) in an NNY codon family, with M<sub>C3:A </sub>= 2M<sub>U3:G</sub>. The shaded area corresponds to P<sub>C </sub>smaller than 1/2 but larger than 1/3.</p></caption></fig>
<table-wrap id="T1" position="float"><label>Table 1</label><caption><p>Number of codon families unambiguously supporting the codon-anticodon adaptation hypothesis (N<sub>CAAH</sub>) and the wobble versatility hypothesis (N<sub>WVH</sub>) in each fungal species.</p></caption><table frame="hsides" rules="groups"><thead><tr><th align="left">Species</th><th align="left"><bold>Accession</bold>*</th><th align="right"><bold>Code</bold><sup>†</sup></th><th align="right">N<sub>CAAH</sub></th><th align="right">N<sub>WVH</sub></th></tr></thead><tbody><tr><td align="left"><italic>Allomyces macrogynus</italic></td><td align="left"><ext-link ext-link-type="gen" xlink:href="NC_001715">NC_001715</ext-link></td><td align="right">4</td><td align="right">1</td><td align="right">13</td></tr><tr><td align="left"><italic>Ashbya gossypii ATCC 10895</italic></td><td align="left"><ext-link ext-link-type="gen" xlink:href="NC_005789">NC_005789</ext-link></td><td align="right">3</td><td align="right">0</td><td align="right">10</td></tr><tr><td align="left"><italic>Aspergillus niger</italic></td><td align="left"><ext-link ext-link-type="gen" xlink:href="NC_007445">NC_007445</ext-link></td><td align="right">4</td><td align="right">0</td><td align="right">12</td></tr><tr><td align="left"><italic>Aspergillus tubingensis</italic></td><td align="left"><ext-link ext-link-type="gen" xlink:href="NC_007597">NC_007597</ext-link></td><td align="right">4</td><td align="right">0</td><td align="right">12</td></tr><tr><td align="left"><italic>Candida albicans SC5314</italic></td><td align="left"><ext-link ext-link-type="gen" xlink:href="NC_002653">NC_002653</ext-link></td><td align="right">4</td><td align="right">1</td><td align="right">11</td></tr><tr><td align="left"><italic>Candida glabrata</italic></td><td align="left"><ext-link ext-link-type="gen" xlink:href="NC_004691">NC_004691</ext-link></td><td align="right">3</td><td align="right">0</td><td align="right">12</td></tr><tr><td align="left"><italic>Candida metapsilosis</italic></td><td align="left"><ext-link ext-link-type="gen" xlink:href="NC_006971">NC_006971</ext-link></td><td align="right">4</td><td align="right">1</td><td align="right">15</td></tr><tr><td align="left"><italic>Candida orthopsilosis</italic></td><td align="left"><ext-link ext-link-type="gen" xlink:href="NC_006972">NC_006972</ext-link></td><td align="right">4</td><td align="right">1</td><td align="right">15</td></tr><tr><td align="left"><italic>Candida parapsilosis</italic></td><td align="left"><ext-link ext-link-type="gen" xlink:href="NC_005253">NC_005253</ext-link></td><td align="right">4</td><td align="right">1</td><td align="right">15</td></tr><tr><td align="left"><italic>Candida stellata</italic></td><td align="left"><ext-link ext-link-type="gen" xlink:href="NC_005972">NC_005972</ext-link></td><td align="right">4</td><td align="right">1</td><td align="right">13</td></tr><tr><td align="left"><italic>Moniliophthora perniciosa</italic></td><td align="left"><ext-link ext-link-type="gen" xlink:href="NC_005927">NC_005927</ext-link></td><td align="right">4</td><td align="right">0</td><td align="right">13</td></tr><tr><td align="left"><italic>Epidermophyton floccosum</italic></td><td align="left"><ext-link ext-link-type="gen" xlink:href="NC_007394">NC_007394</ext-link></td><td align="right">4</td><td align="right">2</td><td align="right">12</td></tr><tr><td align="left"><italic>Harpochytrium sp. JEL94</italic></td><td align="left"><ext-link ext-link-type="gen" xlink:href="NC_004760">NC_004760</ext-link></td><td align="right">4</td><td align="right">1</td><td align="right">3</td></tr><tr><td align="left"><italic>Harpochytrium sp. JEL105</italic></td><td align="left"><ext-link ext-link-type="gen" xlink:href="NC_004623">NC_004623</ext-link></td><td align="right">4</td><td align="right">1</td><td align="right">4</td></tr><tr><td align="left"><italic>Hyaloraphidium curvatum</italic></td><td align="left"><ext-link ext-link-type="gen" xlink:href="NC_003048">NC_003048</ext-link></td><td align="right">4</td><td align="right">1</td><td align="right">4</td></tr><tr><td align="left"><italic>Hypocrea jecorina</italic></td><td align="left"><ext-link ext-link-type="gen" xlink:href="NC_003388">NC_003388</ext-link></td><td align="right">4</td><td align="right">1</td><td align="right">13</td></tr><tr><td align="left"><italic>Kluyveromyces lactis</italic></td><td align="left"><ext-link ext-link-type="gen" xlink:href="NC_006077">NC_006077</ext-link></td><td align="right">4</td><td align="right">0</td><td align="right">13</td></tr><tr><td align="left"><italic>Kluyveromyces thermotolerans</italic></td><td align="left"><ext-link ext-link-type="gen" xlink:href="NC_006626">NC_006626</ext-link></td><td align="right">3</td><td align="right">0</td><td align="right">11</td></tr><tr><td align="left"><italic>Lecanicillium muscarium</italic></td><td align="left"><ext-link ext-link-type="gen" xlink:href="NC_004514">NC_004514</ext-link></td><td align="right">4</td><td align="right">0</td><td align="right">10</td></tr><tr><td align="left"><italic>Monoblepharella sp. JEL 15</italic></td><td align="left"><ext-link ext-link-type="gen" xlink:href="NC_004624">NC_004624</ext-link></td><td align="right">4</td><td align="right">0</td><td align="right">15</td></tr><tr><td align="left"><italic>Mortierella verticillate</italic></td><td align="left"><ext-link ext-link-type="gen" xlink:href="NC_006838">NC_006838</ext-link></td><td align="right">4</td><td align="right">1</td><td align="right">13</td></tr><tr><td align="left"><italic>Penicillium marneffei</italic></td><td align="left"><ext-link ext-link-type="gen" xlink:href="NC_005256">NC_005256</ext-link></td><td align="right">4</td><td align="right">1</td><td align="right">12</td></tr><tr><td align="left"><italic>Pichia Canadensis</italic></td><td align="left"><ext-link ext-link-type="gen" xlink:href="NC_001762">NC_001762</ext-link></td><td align="right">4</td><td align="right">2</td><td align="right">12</td></tr><tr><td align="left"><italic>Podospora anserine</italic></td><td align="left"><ext-link ext-link-type="gen" xlink:href="NC_001329">NC_001329</ext-link></td><td align="right">4</td><td align="right">0</td><td align="right">13</td></tr><tr><td align="left"><italic>Rhizophydium sp.136</italic></td><td align="left"><ext-link ext-link-type="gen" xlink:href="NC_003053">NC_003053</ext-link></td><td align="right">16</td><td align="right">0</td><td align="right">2</td></tr><tr><td align="left"><italic>Rhizopus oryzae</italic></td><td align="left"><ext-link ext-link-type="gen" xlink:href="NC_006836">NC_006836</ext-link></td><td align="right">4</td><td align="right">1</td><td align="right">14</td></tr><tr><td align="left"><italic>Saccharomyces cerevisiae</italic></td><td align="left"><ext-link ext-link-type="gen" xlink:href="NC_001224">NC_001224</ext-link></td><td align="right">3</td><td align="right">1</td><td align="right">12</td></tr><tr><td align="left"><italic>Saccharomyces castellii</italic></td><td align="left"><ext-link ext-link-type="gen" xlink:href="NC_003920">NC_003920</ext-link></td><td align="right">3</td><td align="right">0</td><td align="right">14</td></tr><tr><td align="left"><italic>Saccharomyces servazzii</italic></td><td align="left"><ext-link ext-link-type="gen" xlink:href="NC_004918">NC_004918</ext-link></td><td align="right">3</td><td align="right">0</td><td align="right">9</td></tr><tr><td align="left"><italic>Schizophyllum commune</italic></td><td align="left"><ext-link ext-link-type="gen" xlink:href="NC_003049">NC_003049</ext-link></td><td align="right">4</td><td align="right">1</td><td align="right">13</td></tr><tr><td align="left"><italic>Schizosaccharomyces japonicus</italic></td><td align="left"><ext-link ext-link-type="gen" xlink:href="NC_004332">NC_004332</ext-link></td><td align="right">4</td><td align="right">1</td><td align="right">16</td></tr><tr><td align="left"><italic>Schizosaccharomyces octosporus</italic></td><td align="left"><ext-link ext-link-type="gen" xlink:href="NC_004312">NC_004312</ext-link></td><td align="right">4</td><td align="right">1</td><td align="right">14</td></tr><tr><td align="left"><italic>Schizosaccharomyces pombe</italic></td><td align="left"><ext-link ext-link-type="gen" xlink:href="NC_001326">NC_001326</ext-link></td><td align="right">4</td><td align="right">1</td><td align="right">14</td></tr><tr><td align="left"><italic>Smittium culisetae</italic></td><td align="left"><ext-link ext-link-type="gen" xlink:href="NC_006837">NC_006837</ext-link></td><td align="right">4</td><td align="right">1</td><td align="right">15</td></tr><tr><td align="left"><italic>Spizellomyces punctatus1</italic></td><td align="left"><ext-link ext-link-type="gen" xlink:href="NC_003052">NC_003052</ext-link></td><td align="right">16</td><td align="right">0</td><td align="right">4</td></tr><tr><td align="left"><italic>Yarrowia lipolytica</italic></td><td align="left"><ext-link ext-link-type="gen" xlink:href="NC_002659">NC_002659</ext-link></td><td align="right">3</td><td align="right">0</td><td align="right">11</td></tr><tr><td align="left">Sum</td><td/><td/><td align="right">23</td><td align="right">414</td></tr></tbody></table></table-wrap> <table-wrap id="T2" position="float"><label>Table 2</label><caption><p>N<sub>C</sub>/N<sub>U </sub>ratios for NNY codons and N<sub>A</sub>/N<sub>G </sub>ratios for NNR codons in <italic>Allomyces macrogynus </italic>(representing fungal mitochondrial genomes with translation table 4).</p></caption><table frame="hsides" rules="groups"><thead><tr><th align="left">CF*</th><th align="right">AA<sup>†</sup></th><th align="right">AC<sup>‡</sup></th><th align="right">N<sub>Cod </sub><sup>§</sup></th><th align="right">N<sub>C</sub>/N<sub>U</sub></th><th align="right">CF*</th><th align="right">AA<sup>†</sup></th><th align="right">AC<sup>‡</sup></th><th align="right">N<sub>Cod </sub><sup>§</sup></th><th align="right">N<sub>A</sub>/N<sub>G</sub></th></tr></thead><tbody><tr><td align="left">AAY</td><td align="right">N</td><td align="right">GUU</td><td align="right">360</td><td align="right">0.295</td><td align="right">AAR</td><td align="right">K</td><td align="right">UUU</td><td align="right">411</td><td align="right">2.543</td></tr><tr><td align="left">AGY</td><td align="right">S</td><td align="right">GCU</td><td align="right">230</td><td align="right">0.314</td><td align="right">AGR</td><td align="right">R</td><td align="right">UCU</td><td align="right">202</td><td align="right">6.769</td></tr><tr><td align="left">CAY</td><td align="right">H</td><td align="right">GUG</td><td align="right">213</td><td align="right">0.357</td><td align="right">CAR</td><td align="right">Q</td><td align="right">UUG</td><td align="right">173</td><td align="right">3.806</td></tr><tr><td align="left">GAY</td><td align="right">D</td><td align="right">GUC</td><td align="right">356</td><td align="right">0.299</td><td align="right">GAR</td><td align="right">E</td><td align="right">UUC</td><td align="right">263</td><td align="right">2.372</td></tr><tr><td align="left">UAY</td><td align="right">Y</td><td align="right">GUA</td><td align="right">390</td><td align="right">0.279</td><td align="right">UUR</td><td align="right">L</td><td align="right">UAA</td><td align="right">742</td><td align="right">5.870</td></tr><tr><td align="left">UGY</td><td align="right">C</td><td align="right">GCA</td><td align="right">94</td><td align="right">0.382</td><td/><td/><td/><td/><td/></tr><tr><td align="left">UUY</td><td align="right">F</td><td align="right">GAA</td><td align="right">592</td><td align="right">0.726</td><td/><td/><td/><td/><td/></tr><tr><td colspan="10"><hr/></td></tr><tr><td align="left">Minimum</td><td/><td/><td/><td align="right">0.279</td><td/><td/><td/><td/><td align="right">2.372</td></tr><tr><td align="left">Mean</td><td/><td/><td/><td align="right">0.379</td><td/><td/><td/><td/><td align="right">4.272</td></tr><tr><td align="left">Std Dev</td><td/><td/><td/><td align="right">0.157</td><td/><td/><td/><td/><td align="right">1.975</td></tr></tbody></table></table-wrap> <table-wrap id="T3" position="float"><label>Table 3</label><caption><p>N<sub>C</sub>/N<sub>U </sub>ratios for NNY codons and N<sub>A</sub>/N<sub>G </sub>ratios for NNR codons in <italic>Saccharomyces cerevisiae </italic>(representing fungal mitochondrial genomes with translation table 3).</p></caption><table frame="hsides" rules="groups"><thead><tr><th align="left">CF*</th><th align="left">AA<sup>†</sup></th><th align="left">AC<sup>‡</sup></th><th align="right">N<sub>Cod </sub><sup>§</sup></th><th align="right">N<sub>C</sub>/N<sub>U</sub></th><th align="right">CF*</th><th align="right">AA<sup>†</sup></th><th align="right">AC<sup>‡</sup></th><th align="right">N<sub>Cod </sub><sup>§</sup></th><th align="right">N<sub>A</sub>/N<sub>G</sub></th></tr></thead><tbody><tr><td align="left">AAY</td><td align="left">N</td><td align="left">GUU</td><td align="right">786</td><td align="right">0.065</td><td align="right">AAR</td><td align="right">K</td><td align="right">UUU</td><td align="right">599</td><td align="right">16.114</td></tr><tr><td align="left">AGY</td><td align="left">S</td><td align="left">GCU</td><td align="right">119</td><td align="right">0.035</td><td align="right">AGR</td><td align="right">R</td><td align="right">UCU</td><td align="right">225</td><td align="right">36.500</td></tr><tr><td align="left">AUY</td><td align="left">I</td><td align="left">GAU</td><td align="right">728</td><td align="right">0.093</td><td align="right">AUR</td><td align="right">M</td><td align="right">CAU</td><td align="right">428</td><td align="right">1.326</td></tr><tr><td align="left">CAY</td><td align="left">H</td><td align="left">GUG</td><td align="right">183</td><td align="right">0.070</td><td align="right">CAR</td><td align="right">Q</td><td align="right">UUG</td><td align="right">165</td><td align="right">8.706</td></tr><tr><td align="left">GAY</td><td align="left">D</td><td align="left">GUC</td><td align="right">264</td><td align="right">0.052</td><td align="right">GAR</td><td align="right">E</td><td align="right">UUC</td><td align="right">220</td><td align="right">9.000</td></tr><tr><td align="left">UAY</td><td align="left">Y</td><td align="left">GUA</td><td align="right">495</td><td align="right">0.076</td><td align="right">UGR</td><td align="right">W</td><td align="right">UCA</td><td align="right">130</td><td align="right">20.667</td></tr><tr><td align="left">UGY</td><td align="left">C</td><td align="left">GCA</td><td align="right">92</td><td align="right">0.070</td><td align="right">UUR</td><td align="right">L</td><td align="right">UAA</td><td align="right">911</td><td align="right">64.071</td></tr><tr><td align="left">UUY</td><td align="left">F</td><td align="left">GAA</td><td align="right">445</td><td align="right">0.369</td><td/><td/><td/><td/><td/></tr><tr><td colspan="10"><hr/></td></tr><tr><td align="left">Minimum</td><td/><td/><td/><td align="right">0.035</td><td/><td/><td/><td/><td align="right">1.326</td></tr><tr><td align="left">Mean</td><td/><td/><td/><td align="right">0.104</td><td/><td/><td/><td/><td align="right">22.341</td></tr><tr><td align="left">Std Dev</td><td/><td/><td/><td align="right">0.109</td><td/><td/><td/><td/><td align="right">21.560</td></tr></tbody></table></table-wrap> <table-wrap id="T4" position="float"><label>Table 4</label><caption><p>tRNAs translating two-fold codon families in the nuclear genome of <italic>Saccharomyces cerevisiae</italic>.</p></caption><table frame="hsides" rules="groups"><thead><tr><th align="left">AA*</th><th align="left">Codon<sup>†</sup></th><th align="left">T<sup>‡</sup></th><th align="left">F<sup>§</sup></th></tr></thead><tbody><tr><td align="left">Arg</td><td align="left">AGA</td><td align="left">11</td><td align="left">314</td></tr><tr><td align="left">Arg</td><td align="left">AGG</td><td align="left">1</td><td align="left">1</td></tr><tr><td align="left">Gln</td><td align="left">CAA</td><td align="left">9</td><td align="left">153</td></tr><tr><td align="left">Gln</td><td align="left">CAG</td><td align="left">1</td><td align="left">1</td></tr><tr><td align="left">Glu</td><td align="left">GAA</td><td align="left">14</td><td align="left">305</td></tr><tr><td align="left">Glu</td><td align="left">GAG</td><td align="left">2</td><td align="left">5</td></tr><tr><td align="left">Leu</td><td align="left">UUA</td><td align="left">7</td><td align="left">42</td></tr><tr><td align="left">Leu</td><td align="left">UUG</td><td align="left">10</td><td align="left">359</td></tr><tr><td align="left">Lys</td><td align="left">AAA</td><td align="left">7</td><td align="left">65</td></tr><tr><td align="left">Lys</td><td align="left">AAG</td><td align="left">14</td><td align="left">483</td></tr><tr><td align="left">Asn</td><td align="left">AAC</td><td align="left">10</td><td align="left">208</td></tr><tr><td align="left">Asn</td><td align="left">AAU</td><td align="left">0</td><td align="left">11</td></tr><tr><td align="left">Asp</td><td align="left">GAC</td><td align="left">16</td><td align="left">202</td></tr><tr><td align="left">Asp</td><td align="left">GAU</td><td align="left">0</td><td align="left">112</td></tr><tr><td align="left">Cys</td><td align="left">UGC</td><td align="left">4</td><td align="left">3</td></tr><tr><td align="left">Cys</td><td align="left">UGU</td><td align="left">0</td><td align="left">39</td></tr><tr><td align="left">His</td><td align="left">CAC</td><td align="left">7</td><td align="left">102</td></tr><tr><td align="left">His</td><td align="left">CAU</td><td align="left">0</td><td align="left">25</td></tr><tr><td align="left">Phe</td><td align="left">UUC</td><td align="left">10</td><td align="left">168</td></tr><tr><td align="left">Phe</td><td align="left">UUU</td><td align="left">0</td><td align="left">19</td></tr><tr><td align="left">Ser</td><td align="left">AGC</td><td align="left">2</td><td align="left">6</td></tr><tr><td align="left">Ser</td><td align="left">AGU</td><td align="left">0</td><td align="left">4</td></tr><tr><td align="left">Tyr</td><td align="left">UAC</td><td align="left">8</td><td align="left">141</td></tr><tr><td align="left">Tyr</td><td align="left">UAU</td><td align="left">0</td><td align="left">10</td></tr></tbody></table></table-wrap>
<disp-formula id="bmcM1"><label>(1)</label><mml:math id="M1" name="1471-2148-8-211-i1" overflow="scroll"><mml:semantics><mml:mtable columnalign="left"><mml:mtr><mml:mtd><mml:msub><mml:mi>M</mml:mi><mml:mrow><mml:mi>w</mml:mi><mml:mi>G</mml:mi></mml:mrow></mml:msub><mml:mo>=</mml:mo><mml:msub><mml:mi>N</mml:mi><mml:mi>U</mml:mi></mml:msub><mml:msub><mml:mi>M</mml:mi><mml:mrow><mml:mi>U</mml:mi><mml:mn>3</mml:mn><mml:mo>:</mml:mo><mml:mi>G</mml:mi></mml:mrow></mml:msub></mml:mtd></mml:mtr><mml:mtr><mml:mtd><mml:msub><mml:mi>M</mml:mi><mml:mrow><mml:mi>w</mml:mi><mml:mi>A</mml:mi></mml:mrow></mml:msub><mml:mo>=</mml:mo><mml:msub><mml:mi>N</mml:mi><mml:mi>C</mml:mi></mml:msub><mml:msub><mml:mi>M</mml:mi><mml:mrow><mml:mi>C</mml:mi><mml:mn>3</mml:mn><mml:mo>:</mml:mo><mml:mi>A</mml:mi></mml:mrow></mml:msub><mml:mo>=</mml:mo><mml:msub><mml:mi>N</mml:mi><mml:mi>C</mml:mi></mml:msub><mml:msub><mml:mi>M</mml:mi><mml:mi>O</mml:mi></mml:msub></mml:mtd></mml:mtr></mml:mtable></mml:semantics></mml:math></disp-formula> <disp-formula id="bmcM2"><label>(2)</label><mml:math id="M2" name="1471-2148-8-211-i2" overflow="scroll"><mml:semantics><mml:mrow><mml:mfrac><mml:mrow><mml:msub><mml:mi>M</mml:mi><mml:mrow><mml:mi>U</mml:mi><mml:mn>3</mml:mn><mml:mo>:</mml:mo><mml:mi>G</mml:mi></mml:mrow></mml:msub></mml:mrow><mml:mrow><mml:msub><mml:mi>M</mml:mi><mml:mi>O</mml:mi></mml:msub></mml:mrow></mml:mfrac><mml:mo>=</mml:mo><mml:mfrac><mml:mrow><mml:msub><mml:mi>N</mml:mi><mml:mi>C</mml:mi></mml:msub></mml:mrow><mml:mrow><mml:msub><mml:mi>N</mml:mi><mml:mi>U</mml:mi></mml:msub></mml:mrow></mml:mfrac><mml:mo>=</mml:mo><mml:mfrac><mml:mrow><mml:msub><mml:mi>P</mml:mi><mml:mi>C</mml:mi></mml:msub></mml:mrow><mml:mrow><mml:mn>1</mml:mn><mml:mo>−</mml:mo><mml:msub><mml:mi>P</mml:mi><mml:mi>C</mml:mi></mml:msub></mml:mrow></mml:mfrac></mml:mrow></mml:semantics></mml:math></disp-formula> <disp-formula id="bmcM3"><label>(3)</label><mml:math id="M3" name="1471-2148-8-211-i3" overflow="scroll"><mml:semantics><mml:mtable columnalign="left"><mml:mtr><mml:mtd><mml:msub><mml:mi>N</mml:mi><mml:mi>U</mml:mi></mml:msub><mml:msub><mml:mi>M</mml:mi><mml:mrow><mml:mi>U</mml:mi><mml:mn>3</mml:mn><mml:mo>:</mml:mo><mml:mi>G</mml:mi></mml:mrow></mml:msub><mml:mo>&lt;</mml:mo><mml:msub><mml:mi>N</mml:mi><mml:mi>C</mml:mi></mml:msub><mml:msub><mml:mi>M</mml:mi><mml:mi>O</mml:mi></mml:msub></mml:mtd></mml:mtr><mml:mtr><mml:mtd><mml:mfrac><mml:mrow><mml:msub><mml:mi>M</mml:mi><mml:mrow><mml:mi>U</mml:mi><mml:mn>3</mml:mn><mml:mo>:</mml:mo><mml:mi>G</mml:mi></mml:mrow></mml:msub></mml:mrow><mml:mrow><mml:msub><mml:mi>M</mml:mi><mml:mi>O</mml:mi></mml:msub></mml:mrow></mml:mfrac><mml:mo>&lt;</mml:mo><mml:mfrac><mml:mrow><mml:msub><mml:mi>N</mml:mi><mml:mi>C</mml:mi></mml:msub></mml:mrow><mml:mrow><mml:msub><mml:mi>N</mml:mi><mml:mi>U</mml:mi></mml:msub></mml:mrow></mml:mfrac></mml:mtd></mml:mtr></mml:mtable></mml:semantics></mml:math></disp-formula> <disp-formula id="bmcM4"><label>(4)</label><mml:math id="M4" name="1471-2148-8-211-i4" overflow="scroll"><mml:semantics><mml:mtable columnalign="left"><mml:mtr><mml:mtd><mml:msub><mml:mi>N</mml:mi><mml:mi>U</mml:mi></mml:msub><mml:msub><mml:mi>M</mml:mi><mml:mrow><mml:mi>U</mml:mi><mml:mn>3</mml:mn><mml:mo>:</mml:mo><mml:mi>G</mml:mi></mml:mrow></mml:msub><mml:mo>&gt;</mml:mo><mml:msub><mml:mi>N</mml:mi><mml:mi>C</mml:mi></mml:msub><mml:msub><mml:mi>M</mml:mi><mml:mi>O</mml:mi></mml:msub></mml:mtd></mml:mtr><mml:mtr><mml:mtd><mml:mfrac><mml:mrow><mml:msub><mml:mi>M</mml:mi><mml:mrow><mml:mi>U</mml:mi><mml:mn>3</mml:mn><mml:mo>:</mml:mo><mml:mi>G</mml:mi></mml:mrow></mml:msub></mml:mrow><mml:mrow><mml:msub><mml:mi>M</mml:mi><mml:mi>O</mml:mi></mml:msub></mml:mrow></mml:mfrac><mml:mo>&gt;</mml:mo><mml:mfrac><mml:mrow><mml:msub><mml:mi>N</mml:mi><mml:mi>C</mml:mi></mml:msub></mml:mrow><mml:mrow><mml:msub><mml:mi>N</mml:mi><mml:mi>U</mml:mi></mml:msub></mml:mrow></mml:mfrac></mml:mtd></mml:mtr></mml:mtable></mml:semantics></mml:math></disp-formula> <disp-formula id="bmcM5"><label>(5)</label><mml:math id="M5" name="1471-2148-8-211-i5" overflow="scroll"><mml:semantics><mml:mrow><mml:mfrac><mml:mrow><mml:msub><mml:mi>M</mml:mi><mml:mrow><mml:mi>G</mml:mi><mml:mn>3</mml:mn><mml:mo>:</mml:mo><mml:mi>U</mml:mi></mml:mrow></mml:msub></mml:mrow><mml:mrow><mml:msub><mml:mi>M</mml:mi><mml:mi>O</mml:mi></mml:msub></mml:mrow></mml:mfrac><mml:mo>&lt;</mml:mo><mml:mfrac><mml:mrow><mml:msub><mml:mi>N</mml:mi><mml:mi>A</mml:mi></mml:msub></mml:mrow><mml:mrow><mml:msub><mml:mi>N</mml:mi><mml:mi>G</mml:mi></mml:msub></mml:mrow></mml:mfrac></mml:mrow></mml:semantics></mml:math></disp-formula> <disp-formula id="bmcM6"><label>(6)</label><mml:math id="M6" name="1471-2148-8-211-i6" overflow="scroll"><mml:semantics><mml:mtable columnalign="left"><mml:mtr><mml:mtd><mml:msub><mml:mi>M</mml:mi><mml:mrow><mml:mi>w</mml:mi><mml:mi>A</mml:mi></mml:mrow></mml:msub><mml:mo>=</mml:mo><mml:msub><mml:mi>N</mml:mi><mml:mi>A</mml:mi></mml:msub><mml:msub><mml:mi>M</mml:mi><mml:mrow><mml:mi>A</mml:mi><mml:mn>3</mml:mn><mml:mo>:</mml:mo><mml:mi>A</mml:mi></mml:mrow></mml:msub><mml:mo>+</mml:mo><mml:msub><mml:mi>N</mml:mi><mml:mi>C</mml:mi></mml:msub><mml:msub><mml:mi>M</mml:mi><mml:mrow><mml:mi>C</mml:mi><mml:mn>3</mml:mn><mml:mo>:</mml:mo><mml:mi>A</mml:mi></mml:mrow></mml:msub><mml:mo>=</mml:mo><mml:mo stretchy="false">(</mml:mo><mml:msub><mml:mi>N</mml:mi><mml:mi>A</mml:mi></mml:msub><mml:mo>+</mml:mo><mml:msub><mml:mi>N</mml:mi><mml:mi>C</mml:mi></mml:msub><mml:mo stretchy="false">)</mml:mo><mml:msub><mml:mi>M</mml:mi><mml:mi>O</mml:mi></mml:msub></mml:mtd></mml:mtr><mml:mtr><mml:mtd><mml:msub><mml:mi>M</mml:mi><mml:mrow><mml:mi>w</mml:mi><mml:mi>G</mml:mi></mml:mrow></mml:msub><mml:mo>=</mml:mo><mml:msub><mml:mi>N</mml:mi><mml:mi>A</mml:mi></mml:msub><mml:msub><mml:mi>M</mml:mi><mml:mrow><mml:mi>A</mml:mi><mml:mn>3</mml:mn><mml:mo>:</mml:mo><mml:mi>G</mml:mi></mml:mrow></mml:msub><mml:mo>+</mml:mo><mml:msub><mml:mi>N</mml:mi><mml:mi>U</mml:mi></mml:msub><mml:msub><mml:mi>M</mml:mi><mml:mrow><mml:mi>U</mml:mi><mml:mn>3</mml:mn><mml:mo>:</mml:mo><mml:mi>G</mml:mi></mml:mrow></mml:msub><mml:mo>=</mml:mo><mml:msub><mml:mi>N</mml:mi><mml:mi>A</mml:mi></mml:msub><mml:msub><mml:mi>M</mml:mi><mml:mi>O</mml:mi></mml:msub><mml:mo>+</mml:mo><mml:msub><mml:mi>N</mml:mi><mml:mi>U</mml:mi></mml:msub><mml:msub><mml:mi>M</mml:mi><mml:mrow><mml:mi>U</mml:mi><mml:mn>3</mml:mn><mml:mo>:</mml:mo><mml:mi>G</mml:mi></mml:mrow></mml:msub></mml:mtd></mml:mtr><mml:mtr><mml:mtd><mml:msub><mml:mi>M</mml:mi><mml:mrow><mml:mi>w</mml:mi><mml:mi>C</mml:mi></mml:mrow></mml:msub><mml:mo>=</mml:mo><mml:msub><mml:mi>N</mml:mi><mml:mi>A</mml:mi></mml:msub><mml:msub><mml:mi>M</mml:mi><mml:mrow><mml:mi>A</mml:mi><mml:mn>3</mml:mn><mml:mo>:</mml:mo><mml:mi>C</mml:mi></mml:mrow></mml:msub><mml:mo>+</mml:mo><mml:msub><mml:mi>N</mml:mi><mml:mi>C</mml:mi></mml:msub><mml:msub><mml:mi>M</mml:mi><mml:mrow><mml:mi>C</mml:mi><mml:mn>3</mml:mn><mml:mo>:</mml:mo><mml:mi>C</mml:mi></mml:mrow></mml:msub><mml:mo>+</mml:mo><mml:msub><mml:mi>N</mml:mi><mml:mi>U</mml:mi></mml:msub><mml:msub><mml:mi>M</mml:mi><mml:mrow><mml:mi>U</mml:mi><mml:mn>3</mml:mn><mml:mo>:</mml:mo><mml:mi>C</mml:mi></mml:mrow></mml:msub><mml:mo>=</mml:mo><mml:mo stretchy="false">(</mml:mo><mml:msub><mml:mi>N</mml:mi><mml:mi>A</mml:mi></mml:msub><mml:mo>+</mml:mo><mml:msub><mml:mi>N</mml:mi><mml:mi>C</mml:mi></mml:msub><mml:mo stretchy="false">)</mml:mo><mml:msub><mml:mi>M</mml:mi><mml:mi>O</mml:mi></mml:msub><mml:mo>+</mml:mo><mml:msub><mml:mi>N</mml:mi><mml:mi>U</mml:mi></mml:msub><mml:msub><mml:mi>M</mml:mi><mml:mrow><mml:mi>U</mml:mi><mml:mn>3</mml:mn><mml:mo>:</mml:mo><mml:mi>C</mml:mi></mml:mrow></mml:msub></mml:mtd></mml:mtr><mml:mtr><mml:mtd><mml:msub><mml:mi>M</mml:mi><mml:mrow><mml:mi>w</mml:mi><mml:mi>U</mml:mi></mml:mrow></mml:msub><mml:mo>=</mml:mo><mml:msub><mml:mi>N</mml:mi><mml:mi>C</mml:mi></mml:msub><mml:msub><mml:mi>M</mml:mi><mml:mrow><mml:mi>C</mml:mi><mml:mn>3</mml:mn><mml:mo>:</mml:mo><mml:mi>U</mml:mi></mml:mrow></mml:msub><mml:mo>+</mml:mo><mml:msub><mml:mi>N</mml:mi><mml:mi>U</mml:mi></mml:msub><mml:msub><mml:mi>M</mml:mi><mml:mrow><mml:mi>U</mml:mi><mml:mn>3</mml:mn><mml:mo>:</mml:mo><mml:mi>U</mml:mi></mml:mrow></mml:msub><mml:mo>=</mml:mo><mml:mo stretchy="false">(</mml:mo><mml:msub><mml:mi>N</mml:mi><mml:mi>C</mml:mi></mml:msub><mml:mo>+</mml:mo><mml:msub><mml:mi>N</mml:mi><mml:mi>U</mml:mi></mml:msub><mml:mo stretchy="false">)</mml:mo><mml:msub><mml:mi>M</mml:mi><mml:mrow><mml:mi>Y</mml:mi><mml:mn>3</mml:mn><mml:mo>:</mml:mo><mml:mi>U</mml:mi></mml:mrow></mml:msub></mml:mtd></mml:mtr></mml:mtable></mml:semantics></mml:math></disp-formula> <disp-formula id="bmcM7"><label>(7)</label><mml:math id="M7" name="1471-2148-8-211-i7" overflow="scroll"><mml:semantics><mml:mtable columnalign="left"><mml:mtr><mml:mtd><mml:mfrac><mml:mrow><mml:msub><mml:mi>M</mml:mi><mml:mrow><mml:mi>U</mml:mi><mml:mn>3</mml:mn><mml:mo>:</mml:mo><mml:mi>G</mml:mi></mml:mrow></mml:msub></mml:mrow><mml:mrow><mml:msub><mml:mi>M</mml:mi><mml:mi>O</mml:mi></mml:msub></mml:mrow></mml:mfrac><mml:mo>&lt;</mml:mo><mml:mfrac><mml:mrow><mml:msub><mml:mi>N</mml:mi><mml:mi>C</mml:mi></mml:msub></mml:mrow><mml:mrow><mml:msub><mml:mi>N</mml:mi><mml:mi>U</mml:mi></mml:msub></mml:mrow></mml:mfrac></mml:mtd></mml:mtr><mml:mtr><mml:mtd><mml:msub><mml:mi>N</mml:mi><mml:mi>U</mml:mi></mml:msub><mml:msub><mml:mi>M</mml:mi><mml:mrow><mml:mi>U</mml:mi><mml:mn>3</mml:mn><mml:mo>:</mml:mo><mml:mi>G</mml:mi></mml:mrow></mml:msub><mml:mo>&lt;</mml:mo><mml:mo stretchy="false">(</mml:mo><mml:msub><mml:mi>N</mml:mi><mml:mi>C</mml:mi></mml:msub><mml:mo>+</mml:mo><mml:msub><mml:mi>N</mml:mi><mml:mi>U</mml:mi></mml:msub><mml:mo stretchy="false">)</mml:mo><mml:msub><mml:mi>M</mml:mi><mml:mrow><mml:mi>X</mml:mi><mml:mn>3</mml:mn><mml:mo>:</mml:mo><mml:mi>U</mml:mi></mml:mrow></mml:msub><mml:mo>−</mml:mo><mml:msub><mml:mi>N</mml:mi><mml:mi>A</mml:mi></mml:msub><mml:msub><mml:mi>M</mml:mi><mml:mi>O</mml:mi></mml:msub></mml:mtd></mml:mtr><mml:mtr><mml:mtd><mml:mfrac><mml:mrow><mml:msub><mml:mi>M</mml:mi><mml:mrow><mml:mi>U</mml:mi><mml:mn>3</mml:mn><mml:mo>:</mml:mo><mml:mi>G</mml:mi></mml:mrow></mml:msub></mml:mrow><mml:mrow><mml:msub><mml:mi>M</mml:mi><mml:mrow><mml:mi>Y</mml:mi><mml:mn>3</mml:mn><mml:mo>:</mml:mo><mml:mi>U</mml:mi></mml:mrow></mml:msub></mml:mrow></mml:mfrac><mml:mo>&lt;</mml:mo><mml:mfrac><mml:mrow><mml:msub><mml:mi>N</mml:mi><mml:mi>C</mml:mi></mml:msub><mml:mo>+</mml:mo><mml:msub><mml:mi>N</mml:mi><mml:mi>U</mml:mi></mml:msub><mml:mo>−</mml:mo><mml:msub><mml:mi>N</mml:mi><mml:mi>A</mml:mi></mml:msub></mml:mrow><mml:mrow><mml:msub><mml:mi>N</mml:mi><mml:mi>U</mml:mi></mml:msub></mml:mrow></mml:mfrac><mml:mo>,</mml:mo><mml:mtext>assuming </mml:mtext><mml:msub><mml:mi>M</mml:mi><mml:mrow><mml:mi>Y</mml:mi><mml:mn>3</mml:mn><mml:mo>:</mml:mo><mml:mi>U</mml:mi></mml:mrow></mml:msub><mml:mo>&lt;</mml:mo><mml:msub><mml:mi>M</mml:mi><mml:mi>O</mml:mi></mml:msub></mml:mtd></mml:mtr></mml:mtable></mml:semantics></mml:math></disp-formula> <disp-formula id="bmcM8"><label>(8)</label><mml:math id="M8" name="1471-2148-8-211-i8" overflow="scroll"><mml:semantics><mml:mtable columnalign="left"><mml:mtr><mml:mtd><mml:msub><mml:mi>M</mml:mi><mml:mrow><mml:mi>w</mml:mi><mml:mi>A</mml:mi></mml:mrow></mml:msub><mml:mo>=</mml:mo><mml:msub><mml:mi>N</mml:mi><mml:mi>A</mml:mi></mml:msub><mml:msub><mml:mi>M</mml:mi><mml:mrow><mml:mi>A</mml:mi><mml:mn>3</mml:mn><mml:mo>:</mml:mo><mml:mi>A</mml:mi></mml:mrow></mml:msub><mml:mo>+</mml:mo><mml:msub><mml:mi>N</mml:mi><mml:mi>C</mml:mi></mml:msub><mml:msub><mml:mi>M</mml:mi><mml:mrow><mml:mi>C</mml:mi><mml:mn>3</mml:mn><mml:mo>:</mml:mo><mml:mi>A</mml:mi></mml:mrow></mml:msub><mml:mo>+</mml:mo><mml:msub><mml:mi>N</mml:mi><mml:mi>G</mml:mi></mml:msub><mml:msub><mml:mi>M</mml:mi><mml:mrow><mml:mi>G</mml:mi><mml:mn>3</mml:mn><mml:mo>:</mml:mo><mml:mi>A</mml:mi></mml:mrow></mml:msub><mml:mo>=</mml:mo><mml:mo stretchy="false">(</mml:mo><mml:msub><mml:mi>N</mml:mi><mml:mi>A</mml:mi></mml:msub><mml:mo>+</mml:mo><mml:msub><mml:mi>N</mml:mi><mml:mi>C</mml:mi></mml:msub><mml:mo>+</mml:mo><mml:msub><mml:mi>N</mml:mi><mml:mi>G</mml:mi></mml:msub><mml:mo stretchy="false">)</mml:mo><mml:msub><mml:mi>M</mml:mi><mml:mi>O</mml:mi></mml:msub></mml:mtd></mml:mtr><mml:mtr><mml:mtd><mml:msub><mml:mi>M</mml:mi><mml:mrow><mml:mi>w</mml:mi><mml:mi>C</mml:mi></mml:mrow></mml:msub><mml:mo>=</mml:mo><mml:msub><mml:mi>N</mml:mi><mml:mi>A</mml:mi></mml:msub><mml:msub><mml:mi>M</mml:mi><mml:mrow><mml:mi>A</mml:mi><mml:mn>3</mml:mn><mml:mo>:</mml:mo><mml:mi>C</mml:mi></mml:mrow></mml:msub><mml:mo>+</mml:mo><mml:msub><mml:mi>N</mml:mi><mml:mi>C</mml:mi></mml:msub><mml:msub><mml:mi>M</mml:mi><mml:mrow><mml:mi>C</mml:mi><mml:mn>3</mml:mn><mml:mo>:</mml:mo><mml:mi>C</mml:mi></mml:mrow></mml:msub><mml:mo>+</mml:mo><mml:msub><mml:mi>N</mml:mi><mml:mi>U</mml:mi></mml:msub><mml:msub><mml:mi>M</mml:mi><mml:mrow><mml:mi>U</mml:mi><mml:mn>3</mml:mn><mml:mo>:</mml:mo><mml:mi>C</mml:mi></mml:mrow></mml:msub><mml:mo>=</mml:mo><mml:mo stretchy="false">(</mml:mo><mml:msub><mml:mi>N</mml:mi><mml:mi>A</mml:mi></mml:msub><mml:mo>+</mml:mo><mml:msub><mml:mi>N</mml:mi><mml:mi>C</mml:mi></mml:msub><mml:mo stretchy="false">)</mml:mo><mml:msub><mml:mi>M</mml:mi><mml:mi>O</mml:mi></mml:msub><mml:mo>+</mml:mo><mml:msub><mml:mi>N</mml:mi><mml:mi>U</mml:mi></mml:msub><mml:msub><mml:mi>M</mml:mi><mml:mrow><mml:mi>U</mml:mi><mml:mn>3</mml:mn><mml:mo>:</mml:mo><mml:mi>C</mml:mi></mml:mrow></mml:msub></mml:mtd></mml:mtr><mml:mtr><mml:mtd><mml:msub><mml:mi>M</mml:mi><mml:mrow><mml:mi>w</mml:mi><mml:mi>G</mml:mi></mml:mrow></mml:msub><mml:mo>=</mml:mo><mml:msub><mml:mi>N</mml:mi><mml:mi>A</mml:mi></mml:msub><mml:msub><mml:mi>M</mml:mi><mml:mrow><mml:mi>A</mml:mi><mml:mn>3</mml:mn><mml:mo>:</mml:mo><mml:mi>G</mml:mi></mml:mrow></mml:msub><mml:mo>+</mml:mo><mml:msub><mml:mi>N</mml:mi><mml:mi>G</mml:mi></mml:msub><mml:msub><mml:mi>M</mml:mi><mml:mrow><mml:mi>G</mml:mi><mml:mn>3</mml:mn><mml:mo>:</mml:mo><mml:mi>G</mml:mi></mml:mrow></mml:msub><mml:mo>+</mml:mo><mml:msub><mml:mi>N</mml:mi><mml:mi>U</mml:mi></mml:msub><mml:msub><mml:mi>M</mml:mi><mml:mrow><mml:mi>U</mml:mi><mml:mn>3</mml:mn><mml:mo>:</mml:mo><mml:mi>G</mml:mi></mml:mrow></mml:msub><mml:mo>=</mml:mo><mml:mo stretchy="false">(</mml:mo><mml:msub><mml:mi>N</mml:mi><mml:mi>A</mml:mi></mml:msub><mml:mo>+</mml:mo><mml:msub><mml:mi>N</mml:mi><mml:mi>G</mml:mi></mml:msub><mml:mo stretchy="false">)</mml:mo><mml:msub><mml:mi>M</mml:mi><mml:mi>O</mml:mi></mml:msub><mml:mo>+</mml:mo><mml:msub><mml:mi>N</mml:mi><mml:mi>U</mml:mi></mml:msub><mml:msub><mml:mi>M</mml:mi><mml:mrow><mml:mi>U</mml:mi><mml:mn>3</mml:mn><mml:mo>:</mml:mo><mml:mi>G</mml:mi></mml:mrow></mml:msub></mml:mtd></mml:mtr><mml:mtr><mml:mtd><mml:msub><mml:mi>M</mml:mi><mml:mrow><mml:mi>w</mml:mi><mml:mi>U</mml:mi></mml:mrow></mml:msub><mml:mo>=</mml:mo><mml:msub><mml:mi>N</mml:mi><mml:mi>C</mml:mi></mml:msub><mml:msub><mml:mi>M</mml:mi><mml:mrow><mml:mi>C</mml:mi><mml:mn>3</mml:mn><mml:mo>:</mml:mo><mml:mi>U</mml:mi></mml:mrow></mml:msub><mml:mo>+</mml:mo><mml:msub><mml:mi>N</mml:mi><mml:mi>G</mml:mi></mml:msub><mml:msub><mml:mi>M</mml:mi><mml:mrow><mml:mi>G</mml:mi><mml:mn>3</mml:mn><mml:mo>:</mml:mo><mml:mi>U</mml:mi></mml:mrow></mml:msub><mml:mo>+</mml:mo><mml:msub><mml:mi>N</mml:mi><mml:mi>U</mml:mi></mml:msub><mml:msub><mml:mi>M</mml:mi><mml:mrow><mml:mi>U</mml:mi><mml:mn>3</mml:mn><mml:mo>:</mml:mo><mml:mi>U</mml:mi></mml:mrow></mml:msub><mml:mo>=</mml:mo><mml:mo stretchy="false">(</mml:mo><mml:msub><mml:mi>N</mml:mi><mml:mi>C</mml:mi></mml:msub><mml:mo>+</mml:mo><mml:msub><mml:mi>N</mml:mi><mml:mi>U</mml:mi></mml:msub><mml:mo stretchy="false">)</mml:mo><mml:msub><mml:mi>M</mml:mi><mml:mrow><mml:mi>Y</mml:mi><mml:mn>3</mml:mn><mml:mo>:</mml:mo><mml:mi>U</mml:mi></mml:mrow></mml:msub><mml:mo>+</mml:mo><mml:msub><mml:mi>N</mml:mi><mml:mi>G</mml:mi></mml:msub><mml:msub><mml:mi>M</mml:mi><mml:mrow><mml:mi>G</mml:mi><mml:mn>3</mml:mn><mml:mo>:</mml:mo><mml:mi>U</mml:mi></mml:mrow></mml:msub></mml:mtd></mml:mtr></mml:mtable></mml:semantics></mml:math></disp-formula> <disp-formula id="bmcM9"><label>(9)</label><mml:math id="M9" name="1471-2148-8-211-i9" overflow="scroll"><mml:semantics><mml:mtable columnalign="left"><mml:mtr><mml:mtd><mml:mfrac><mml:mrow><mml:msub><mml:mi>M</mml:mi><mml:mrow><mml:mi>U</mml:mi><mml:mn>3</mml:mn><mml:mo>:</mml:mo><mml:mi>G</mml:mi></mml:mrow></mml:msub></mml:mrow><mml:mrow><mml:msub><mml:mi>M</mml:mi><mml:mi>O</mml:mi></mml:msub></mml:mrow></mml:mfrac><mml:mo>&gt;</mml:mo><mml:mfrac><mml:mrow><mml:msub><mml:mi>N</mml:mi><mml:mi>C</mml:mi></mml:msub></mml:mrow><mml:mrow><mml:msub><mml:mi>N</mml:mi><mml:mi>U</mml:mi></mml:msub></mml:mrow></mml:mfrac></mml:mtd></mml:mtr><mml:mtr><mml:mtd><mml:mfrac><mml:mrow><mml:msub><mml:mi>M</mml:mi><mml:mrow><mml:mi>U</mml:mi><mml:mn>3</mml:mn><mml:mo>:</mml:mo><mml:mi>C</mml:mi></mml:mrow></mml:msub></mml:mrow><mml:mrow><mml:msub><mml:mi>M</mml:mi><mml:mi>O</mml:mi></mml:msub></mml:mrow></mml:mfrac><mml:mo>&gt;</mml:mo><mml:mfrac><mml:mrow><mml:msub><mml:mi>N</mml:mi><mml:mi>G</mml:mi></mml:msub></mml:mrow><mml:mrow><mml:msub><mml:mi>N</mml:mi><mml:mi>U</mml:mi></mml:msub></mml:mrow></mml:mfrac></mml:mtd></mml:mtr></mml:mtable></mml:semantics></mml:math></disp-formula> <disp-formula id="bmcM10"><label>(10)</label><mml:math id="M10" name="1471-2148-8-211-i10" overflow="scroll"><mml:semantics><mml:mtable columnalign="left"><mml:mtr><mml:mtd><mml:mfrac><mml:mrow><mml:msub><mml:mi>M</mml:mi><mml:mrow><mml:mi>U</mml:mi><mml:mn>3</mml:mn><mml:mo>:</mml:mo><mml:mi>G</mml:mi></mml:mrow></mml:msub></mml:mrow><mml:mrow><mml:msub><mml:mi>M</mml:mi><mml:mi>O</mml:mi></mml:msub></mml:mrow></mml:mfrac><mml:mo>&lt;</mml:mo><mml:mfrac><mml:mrow><mml:msub><mml:mi>N</mml:mi><mml:mi>C</mml:mi></mml:msub></mml:mrow><mml:mrow><mml:msub><mml:mi>N</mml:mi><mml:mi>U</mml:mi></mml:msub></mml:mrow></mml:mfrac><mml:mtext> and</mml:mtext></mml:mtd></mml:mtr><mml:mtr><mml:mtd><mml:mfrac><mml:mrow><mml:msub><mml:mi>M</mml:mi><mml:mrow><mml:mi>G</mml:mi><mml:mn>3</mml:mn><mml:mo>:</mml:mo><mml:mi>U</mml:mi></mml:mrow></mml:msub></mml:mrow><mml:mrow><mml:msub><mml:mi>M</mml:mi><mml:mi>O</mml:mi></mml:msub></mml:mrow></mml:mfrac><mml:mo>&lt;</mml:mo><mml:mfrac><mml:mrow><mml:msub><mml:mi>N</mml:mi><mml:mi>A</mml:mi></mml:msub></mml:mrow><mml:mrow><mml:msub><mml:mi>N</mml:mi><mml:mi>G</mml:mi></mml:msub></mml:mrow></mml:mfrac></mml:mtd></mml:mtr></mml:mtable></mml:semantics></mml:math></disp-formula>
<table-wrap-foot><p>* GenBank accession number.</p><p><sup>† </sup>Genetic code.</p></table-wrap-foot> <table-wrap-foot><p>* CF: Codon family</p><p><sup>† </sup>AA: One-letter code of amino acid</p><p><sup>‡ </sup>AC: Anticodon</p><p><sup>§ </sup>N<sub>Cod</sub>: Number of codons in the codon family</p></table-wrap-foot> <table-wrap-foot><p>* CF: Codon family</p><p><sup>† </sup>AA: One-letter code of amino acid</p><p><sup>‡ </sup>AC: Anticodon</p><p><sup>§ </sup>N<sub>Cod</sub>: Number of codons in the codon family</p></table-wrap-foot> <table-wrap-foot><p>* Amion acid carried by tRNA</p><p><sup>† </sup>Codons forming Watson-Crick base pair with the anticodon of tRNA</p><p><sup>‡ </sup>Copy number of tRNA gene with anticodon matching the corresponding codon.</p><p><sup>§ </sup>Codon frequencies of highly expressed yeast protein-coding genes compiled in the Eyeastcai.cut file distributed with EMBOSS [##REF##10827456##44##].</p></table-wrap-foot>
<graphic xlink:href="1471-2148-8-211-1"/>
[{"surname": ["Martin", "Sibler", "Gehrke", "Kuo", "Edmonds", "McCloskey", "Dirheimer"], "given-names": ["RP", "AP", "CW", "K", "CG", "JA", "G"], "article-title": ["5-[[(carboxymethyl)amino]methyl]uridine is found in the anticodon of yeast mitochondrial tRNAs recognizing two-codon families ending in a purine"], "source": ["Biochemistry (Mosc)"], "year": ["1990"], "volume": ["29"], "issue": ["4"], "fpage": ["956"], "lpage": ["959"], "pub-id": ["10.1021/bi00456a016"]}, {"surname": ["Yokoyama", "Nishimura"], "given-names": ["S", "S"], "person-group": ["Soll D, RajBhandary U"], "article-title": ["Modified nucleotides and codon recognition."], "source": ["tRNA: structure, biosynthesis and function"], "year": ["1995"], "publisher-name": ["Washington , ASM Press"], "fpage": ["207"], "lpage": ["223"]}, {"surname": ["Rocha"], "given-names": ["EP"], "article-title": ["Codon usage bias from tRNA's point of view: Redundancy, specialization, and efficient decoding for translation optimization"], "source": ["Genome Res"], "year": ["2004"], "volume": ["12"], "fpage": ["12"]}, {"surname": ["Xia"], "given-names": ["X"], "source": ["Data analysis in molecular biology and evolution."], "year": ["2001"], "publisher-name": ["Boston , Kluwer Academic Publishers"], "fpage": ["277"]}]
BMC Evol Biol. 2008 Jul 19; 8:211