{'fraction_non_alphanumeric': 0.03675562593523456, 'fraction_numerical': 0.022866938822948296, 'mean_word_length': 4.933965517241379, 'pattern_counts': {'":': 0, '<': 0, '<?xml version=': 0, '>': 0, 'https://': 0, 'lorem ipsum': 0, 'www.': 0, 'xml': 0}, 'pii_count': 0, 'substrings_counts': 0, 'word_list_counts': {'cursed_substrings.json': 3, 'profanity_word_list.json': 0, 'sexual_word_list.json': 0, 'zh_pornsignals.json': 0}}
An integrated, structure- and energy-based view of the genetic code

The principles of mRNA decoding are conserved among all extant life forms. We present an integrative view of all the interaction networks between mRNA, tRNA and rRNA: the intrinsic stability of codon–anticodon duplex, the conformation of the anticodon hairpin, the presence of modified nucleotides, the occurrence of non-Watson–Crick pairs in the codon–anticodon helix and the interactions with bases of rRNA at the A-site decoding site. We derive a more information-rich, alternative representation of the genetic code, that is circular with an unsymmetrical distribution of codons leading to a clear segregation between GC-rich 4-codon boxes and AU-rich 2:2-codon and 3:1-codon boxes. All tRNA sequence variations can be visualized, within an internal structural and energy framework, for each organism, and each anticodon of the sense codons. The multiplicity and complexity of nucleotide modifications at positions 34 and 37 of the anticodon loop segregate meaningfully, and correlate well with the necessity to stabilize AU-rich codon–anticodon pairs and to avoid miscoding in split codon boxes. The evolution and expansion of the genetic code is viewed as being originally based on GC content with progressive introduction of A/U together with tRNA modifications. The representation we present should help the engineering of the genetic code to include non-natural amino acids.


INTRODUCTION
The genetic code is considered as 'quasi' universal. Although deviations of the 'standard' genetic code have been identified, they occur in only a few organisms and mainly in organelles (1)(2)(3). However, as more organisms are being presently studied, new exceptions will certainly be discovered. Nevertheless, because all extant organisms evolved from a small set of ancestral living organisms that, along evolution, continued to exchange DNA information, the basic molecular principles of decoding (including the genetic code) remained essentially the same in all organisms of the three domains of life (Bacteria, Eukarya and Archaea). These principles are anchored in the physics and chemistry of the interplay between the diverse, and often locally neutral, non-covalent and weak molecular interactions engaged, and in the thermodynamics of formation of the several distinct complexes present at the various stages of the multistep translational process. Far from excluding profound improvements, these principles allow for adaptation and even innovations of many of the elements of the modern translation machinery (including the genetic code) according to physiological needs and/or to cellular niches and environments.
The traditional way to represent the 64 codons, as organized in the standard codon table (4)(5)(6)(7) took several years to settle down. In the first row of the codon table are indicated all codons starting with a U at the first codon position, with a C in the second row, and a A or a G for the third and fourth rows, respectively. For each column, the same (base) order is used for the second codon base, as well as within each of the 16 decoding boxes for the classification of the 4 codons ending with a different base. This standard and quasi universal codon table is divided into boxes with 4 synonymous codons (hereafter designated as 4-codon boxes or unsplit decoding boxes) where the third base can be either U, C, A or G for coding the same amino acid, and boxes with 2 synonymous codons (hereafter designated as 2-codon boxes or split decoding boxes) where the third base codes for a different amino acid depending on its nature, either pyrimidine (Y) or purine (R). A third category concerns split 3:1-codon boxes where a single codon ending with guanine means methionine or tryptophan (1-codon box). Within these split 2:2-and 3:1-codon boxes are those meaning STOP for opal UGA, ochre UAA and amber UAG. At the time the table was first proposed, knowledge about the detailed molecular mechanisms underlying the decoding process did not exist. A few research groups argued for a different organization of the initial and now standard codon table (8)(9)(10)(11). Arguments were based on various criteria or regularities of the code but never with knowledge of all the elements of the whole decoding process, in particular the large diversity of post-transcriptional nucleotide modifications that exist especially in the tRNA.
Recent progress in identifying modified nucleotides and their functions in tRNA (12)(13)(14), and information gained from the recent detailed structural, physicochemical and kinetic studies of ribosomes associated with mRNA/aminoacyl-tRNA, has made clear that selection of a given codon by its cognate aminoacyl-tRNA harboring a complementary or nearly complementary anticodon within the ribosomal decoding center is a complex and intricate process (15)(16)(17)(18)(19)(20)(21)(22). The ribosome is not a passive machine but an active component (structurally and kinetically speaking) of the codon selection mechanism. Also, the original wobble 'hypothesis' (23) has to be broadened in several cases involving non-standard base pairings at the third base pair of the codon/anticodon helix (see for examples: (24)(25)(26)). Indeed, it is now becoming clear that some base-base oppositions avoid usual pairings or even H-bonding while still being accommodated within a Watson-Crick-like helix. Indeed, post-transcriptional modifications of nucleotide-34 can change the physicochemical behavior of the base (frequency of tautomerism) and/or the spatial preference of the nucleotide (syn/anti, puckering of the ribose) to precisely allow the third anticodon nucleotide to fit within a mini-helix structure together with the two other base-pairs of the anticodon (15,(27)(28)(29)(30)(31)(32). In short, wobbling at the third position does occur in various ways, but without movement of the codon base toward the minor groove (see below).
In the present work, we took into consideration all the new elements that we now understand are important for the aminoacyl-tRNA selection at the A site of the decoding center of the ribosome in order to propose a visually useful decoding table embedding multiple structural aspects of translation. This new representation allows one to propose a rationale for the emergence of the contemporary G/C/A/U containing quasi-universal genetic code from an ancient GC-rich code and the multiple constraints on such an evolution.

THE CODON TABLE
We started from an observation, already made at least partly by others (8,9,33,34), on the distribution of codons within the codon table. The codons within the unsplit 4-codon boxes always have one C in the second or either C or G at the first position of codon leading to C = G (G = C) pair in the codon/anticodon triplets. Among these codon boxes, the two conditions are simultaneously satisfied for four amino acids, Ala, Arg, Gly and Pro. In contrast, codons of the split Figure 1. Circular representation of the genetic code emphasizing the inherent regularities of the decoding recognition process. The codons containing solely G = C pairs at the first two positions are at the top, those containing solely A-U pairs at the bottom, and those with mixed pairs of G = C and A-U either at the first or second pair of the codon/anticodon helix in the middle at the right and left. Thick red lines separate the three main regions. The red arrow indicates the direction of rotation for C1, G1, U1, A1 and the blue arrows the direction of rotation for C2, G2, U2, A2 on the right and left parts of the wheel. The amino acids coded by unsplit 4-codon boxes are indicated in red and those by split 2:2-and 3:1-codon boxes, together with the usual stop codons, are indicated in black. Throughout, the codon positions are numbered B1-B2-B3 and the anticodon nucleotides B34-B35-B36, both from 5 to 3 .
2:2-and 3:1-codon boxes have, instead, A at the second or either A or U at the first position of the codon leading to a A-U (U-A) pair in the codon/anticodon triplets. Among these, the two conditions are simultaneously satisfied for seven amino acids, Asn, Ile, Leu, Lys, Met, Phe and Tyr. The nine remaining amino acids (but 12 codon boxes because Leu, Ser and Arg each have 2 distinct codon boxes) have a mixture of C = G and U-A pairs, with either a C = G pair at the first position followed by a U-A pair or a U-A pair at the first position followed by a C = G pair. The two additional amino acids that are co-translationally inserted, pyrrolysine and selenocysteine, are recoded to replace stop codons, and will be discussed later in our analysis.
In our new representation, we arrange each codon corresponding to the 20 canonical amino acids on a circle (genetic code wheel) as a function of the nature of the first two base pairs of the codon/anticodon triplet. On the circle, the G = C only codons are disposed at the top, the A-U only codons at the bottom, and the mixed ones in-between on the left and the right. The order of third codon base has been chosen so that proximities of purine-and pyrimidineending codons are preserved ( Figure 1). Note how the bases C, G, U, A rotate around the circle: for the first position the C1 starts in the top left quadrant, for the second position C2, G2, A2, U2 rotates down in the right part in a righthanded fashion and in the left part in a left-handed fashion.
The choice of this unsymmetrical disposition will become more apparent and coherent in the following sections.
We then computed the energies of the codon/anticodon minihelices (5 -B1-B2-B3-3 codon paired to 3 -B36-B35-B34-5 anticodon) following the recent Turner values for dimers (35). First, we considered only the first two base pairs, B1-B2 paired to B36-B35 (Figure 2, box A). The values (indicated as bold red numbers) are distributed on a circular view of the codons. Clearly, from those values, three groups of amino acids can be distinguished. A STRONG (S) group (average −3.1 kcal/mole) with, as above, the four Ala, Arg, Gly, Pro. A WEAK (W) group (average −1.0 kcal/mole) with, again as above, the seven Asn, Ile, Leu, Lys, Met, Phe, Tyr. And finally, an INTERMEDI-ATE (I) group (average −2.2 kcal/mole) gathers the remaining amino acids (or codon boxes). A similar distribution is observed by computing the triplets of base pairs (B1-B2-B3 with B36-B35-B34) ( Figure 2, box B), the corresponding values are in parentheses). Indeed, the distribution from STRONG to WEAK is respected (averages from −5.8, −4.3, −2.7 kcal/mole). Finally, in order to be slightly closer to biological reality, we considered the U3oG34 and G3oU34 Crick-type of wobbles at the third base pair (Figure 2, boxes C and D, respectively). Again, the three groups segregate well with energies, respectively, of −5.0, −3.5 and −1.8 kcal/mole for the U3oG34 wobble and −4.6, −3.0, and −1.4 kcal/mole for the G3oU34 wobble, respectively. Interestingly, the G3oU34 wobble is always less stable than the U3oG34 wobble interactions (see discussion below). There is no overlap between the three categories, in the sense that no INTERMEDIATE codon has a stronger energy than any STRONG codon or a weaker one than any WEAK codon.
Despite the large differences of energy between GCrich and AU-rich codon-anticodon duplexes, the measured binding constants of natural aminoacyl-tRNAs on 70S ribosomes programmed with mRNA are unexpectedly remarkably similar (20)(21)(22)36,37). Likewise the stability of complexes formed between two natural tRNAs harboring complementary anticodons is very much the same whatever the GC/AU composition of the duplex (38,39). Only when non-natural mutant tRNAs are used are large differences of binding constants measured. Obviously, both the ribosome and the built-in features of natural tRNAs have evolved to allow minimization of the differential thermodynamic contributions of the codon-anticodon interactions, allowing tRNA-mRNA interactions to be rather uniform; in other words, to be neither too sticky nor too loose during the process of translation. We will now describe these features, first on the ribosome and then on the transfer RNAs.

THE RIBOSOMAL DECODING SITE AND THE ANTI-CODON TRIPLET
Recent crystallographic data of aminoacyl-tRNAs associated to mRNA programmed ribosomes revealed the contacts made by the ribosomal A-site with the codon/anticodon helix, designated as the ribosomal grip. The two invariant adenine residues A1492 and A1493 of helix 44 of 16S rRNA make A-minor type interactions in the minor groove of the first two base pairs of the codon/anticodon helix ( Figure 3) (15)(16)(17)(27)(28)(29)31,32,(40)(41)(42)(43). Except for those made to the B1-B36 pair, such contacts are almost identical for all types of Watson-Crick pairs and, therefore, can be considered to contribute similarly to the global interaction energies (44) (Table 1). On the circular representation (Figure 3), the A-minor contacts are equally partitioned (listed in plain red boxes), except, importantly, for those between A1493 and B1-B36 (G36 = C1 or G1 = C36, in dashed red boxes), which are specific to the top half of the circle.
In Table 2 and Figure 3 are indicated the constant and variable contacts that occur within the anticodon loop and with the codon/anticodon minihelix (45) (see also Supplementary Figure S1A and B). The contributions of such contacts to aminoacyl-tRNA (aa-tRNA) binding to the ribosome have been established by various methods (46)(47)(48). Recent experiments on the effects of 2'O-methyl group in the ribose of mRNA show the strongest reduction in translation efficiency when the 2'O-methyl group is present at the second position, followed by methylation at the third position (25% of unmodified) and the first position (80% of unmodified), confirming further the O2'hydroxyl group donor contacts at the second and third positions (49). They involve mainly the ribose O2' of U33 (an almost invariant nucleotide) with the Hoogsteen edge of B35 (where the contact U33(O2'). . . R35(N7) is stronger than the U33(O2'). . . Y35(C5) contact) (see Tables 1 and 2 and Supplementary Figure S2). Of note is that these contacts therefore slightly disfavor the codon/anticodon interactions with a purine at B2 (Arg, Gly, His/Gln, Asp/Glu, Arg/ser, Cys/Trp) that are circled in green in Figure 3 (discussed below). In short, the ribosomal grip interactions add to the stabilization of mainly codons of the unsplit 4-codon boxes, especially through the A1493(N3) H-bond to the first G = C/C = G pair in the upper part of the circle in Figure  3.
In addition, the crystal structures of ribosomal complexes with mRNA and aminoacyl-tRNAs reveal many contacts between a few amino acids of r-proteins and the RNA sugar-phosphate backbones. The majority of these contacts should be general and applicable to any tRNA and mRNA and, thus, should not make a major contribution to the specificity and fidelity of a specific codon-anticodon binding, at least in the decoding A-site of the ribosome. However, r-proteins S4, S5 and S12 have been demonstrated to be involved in the ribosomal grip by controlling the closure of 30S subunit upon correct codon recognition (50)(51)(52). The protein S12 especially closely contacts the decoding site with Ser46 binding to the Hoogsteen edge of A1492 (which monitors the second base pair B35-B2) and Pro44 interacting via a solvent molecule with the hydroxyl group of the ribose of B3, together with G530 and C518 (17,53). Among the many conserved and semi-conserved residues present in the rRNA region of the decoding center, a few are posttranscriptionally modified (54,55). However, none of these modified residues is in direct contact with the codon or anticodon in the A-site of the ribosome and therefore should not contribute directly to specific codon selection mechanism. The main functions of these post-transcriptional modifications appear to reside in the modulation of rRNA structure and ligand interactions leading to indirect and   Notice that no variable interaction occurs between B2-B35 and the conserved rRNA nucleotides while no interaction occurs between B1-B36 and the other tRNA nucleotides. The contact indicated by * is present in usual A-minor interactions, but presents a long distance in structures for the A-site (3.6Å or longer). Note also that in standard A-minor contacts, the equivalent C1 = G36 is observed more often than G1 = C36 and C1 = G36 more often than G1 = C36 (44). general effects on translational fidelity (56). One exception could be the positively charged m 7 G527 in bacteria. It is lo-cated about 5Å from Pro44 of S12 and the conserved base G530, possibly creating an electrostatic environment that influences the electrostatic state of modified B34, especially with hypermodification at U34. In summary, the present organization of the genetic code makes a clear energetic segregation between all GC-rich codons belonging to 4-codon family boxes and all AU-rich codons belonging to 2:2 and 3:1 codon family boxes, with stop-codons being among the latter. The codon boxes with a mixture of G = C and A-U pairs in the first and second positions are intermediate in the various energetic parameters considered. Based on the above criteria, codons of the unsplit 4-codon boxes appear to fulfill enough criteria to bind correctly to their cognate aa-tRNAs within the ribosomal A site. They probably possess enough energy from both the internal codon-anticodon mini helix and the additional external stabilizing grip energies from both the ribosome and the tRNA anticodon loop.
Remarkably, in certain ribosomal contexts (Mycoplasmas, mitochondria), the energetics of the first two base pairs of the codon-anticodon mini helix appears to be sufficient; a single tRNA (containing unmodified U34) is able to decode each of the four GC-rich codons of unsplit decoding boxes (57)(58)(59)(60), although with different efficiencies. This is known as the 2 out of 3 decoding rule (26,33). In contrast, a minimum of two tRNAs (one with G34, a second one with modified U34--see below--and eventually a third one with C34) is required to read AU-rich codons of each split 2:2-decoding boxes (57)(58)(59). Apparently, below a certain energy threshold, the contributions of the first two base pairs are insufficient to sustain translation and a third base pair has to complement the missing energetics. This is probably why the A/U-rich codons of the split 2:2-decoding boxes, where pyrimidine-ending (Y-ending) codons must be distinguished from purine-ending (R-ending) codons, have been Nucleic Acids Research, 2016, Vol. 44, No. 17 8025 selected along evolution for expansion of the genetic code with additional amino acids (see below).

STRUCTURAL CHARACTERISTICS OF THE tRNA ANTICODON HAIRPIN
Transfer RNAs are not passive molecules but active participants during translation. They can modulate their binding and specificity through sequence preferences at key positions, extending beyond the anticodon triplet (61)(62)(63), and through very diverse base modifications, particularly in the anticodon loop and stem (14,64). In the following sections, the contributions of the sequence preferences and base modifications to smooth and uniform decoding are analyzed within the framework just described, with a focus restrained to the anticodon stem and loop. Figure 4 depicts the general distribution of conserved and semi-conserved nucleotides and the presence of modified nucleotides within the anticodon hairpin of 400 native elongator tRNA sequences originating from organisms spanning Bacteria, Eukarya and Archaea (65). Initiator tRNA Met is excluded from the analysis because it participates in a different step of translation and its anticodon hairpin is quite distinct from those observed in elongator tRNAs (66).
The distribution of purines and pyrimidines in elongator tRNAs is clearly asymmetric and the presence of a few conserved and semi conserved bases is evident. Some of the bases along the stem are post-transcriptionally modified to or simple methylated derivatives, which are known to better stack with neighboring bases, thus favoring a locally more rigid conformation (67). In the middle of the stem, base pair B29-B41 is almost always a Watson-Crick pair, and either a G = C or a C = G pair is usually found at positions B30-B40, attesting to the need for a relatively stable helical conformation in this part of the anticodon stem. Depending on the tRNA species, either a C = G or a Apair is preferred at B31-B39; thus the nucleotide entering the stem at the 3 -side is generally either a G39 or a 39, two bases which favor the start of a helical stem (35,68). Therefore, all bases of the anticodon stem and particularly the last three base pairs stabilize helical conformations. In the anticodon loop, the first and last residues, B32 and B38, also present a biased distribution that maintains in this case a non-Watson-Crick pair, with B32 almost exclusively a pyrimidine Y (with a preference for a C or a U, again depending on the tRNA species) and B38 mainly a A followed by a less frequently used C (see Table 1 in (69)). Thus, at variance with B31-B39 pair, B32 and B38 tend to form a non-Watson-Crick pair with a single bifurcated contact between exocyclic O2(Y) and exocyclic N6(A) (or N4(C)) (45,70).
At B37 a strictly conserved purine is present (mostly A37 or a modified A37 derivative). When G37 is present, it is invariably modified into N1-methylated G37 or a Wyosine derivative in tRNA Phe of eukaryotes and Archaea ( Figure  5A). Modifications at A37 depend on the identity of the adjacent nucleotide at position 36; it is generally modified when B36 is either A or U. The purine residue at 37 stacks over the first codon/anticodon base pair (40,71,72) and propagates the stacking continuity to neighboring residues B38, B39 and the anticodon stem (preferential 3 -stacked conformation). Modifications at the cyclic C2 atom and the exocyclic amine at C6 of A37 (and especially the highly hydrophobic tricyclic derivatives of the Wye-base) reinforce the stacking power of A37 with the neighboring nucleotides (12,72). Likewise, at the first position of the anticodon (B34), a plethora of diverse modified nucleotides are found ( Figure 5B, and discussed below). Their functions are mostly related to decoding within the ribosomal decoding site, however some of them, like Cm, Um, s 2 U, , ac 4 C, Gm and Q, also contribute to the pre-structuring of the anticodon loop (31,73,74).

THE PROXIMAL 'EXTENDED ANTICODON' CON-CEPT
More than 30 years ago, Yarus (61) pointed out the strong co-variations between nucleotide B36 of anticodon and several bases of the whole anticodon stem-loop of E. coli tRNAs and introduced the 'extended anticodon' concept. Many experiments further confirmed such an apparent but indirect link between the sequence of the anticodon stem and the performance of the experimentally mutated bacterial tRNAs to efficiently translate certain codons both in vivo and in vitro. However, as more tRNA sequences became available from many different organisms, and more experiments were performed, it became evident that these co-variations were particularly relevant for the nucleotides of the anticodon loop and the last base pair B31-B39 of the anticodon loop (indicated in grey background in Figure 4 and thereafter designated as the proximal 'extended anticodon'). Figure 6A-C corresponds to an update of Yarus' original data in the form of the 'wheel of code' for the 40 elongator tRNAs of the complete E. coli repertoire that read all the 61 codons for the 20 canonical amino acids. In Figure  6A, are displayed the distributions of B31-B39 base pair, B32/B38 in Figure 6B and B37 in Figure 6C (including the acronyms of modified nucleotides). The same analysis has been done for the elongator tRNA repertoire of the archaea Haloferax volcanii and the eukaryote S. cerevisiae and as well as in two extreme situations for which enough sequence information (including for the modified bases) is reasonably well known, Mycoplasma capricolum and human mitochondria (see Supplementary Figure S3). A separate analysis is mandatory for each organism because notable differences exist in the decoding strategies (though not in the genetic codes) of distantly evolutionary related organisms, such as those of Bacteria, Eukarya and Archaea (64,69).
The most frequently found B31-B39 base pair (boxed in black) is C = G, especially in tRNAs of the 4-codon family boxes ( Figure 6A, upper part of the wheel). In contrast, a A-pair is predominantly present in tRNAs of the split 2-codon boxes (lower part of the wheel). Concerning B32/B38, the combination Y32/A38 is the most frequently found one, with both U32 and C32 occasionally modified to or Um and C32 to Cm or s 2 C, respectively ( Figure 6B). These modifications of the base, especially the pseudouridine and the 2'O-methylation of the ribose, are known to limit significantly the conformational flexibility of the nucleotide and thereby locally rigidify the corresponding portion of the nucleic acid (75).  (14) that currently contains all sequences available by 2009 in the tRNA database as in ref. (76). Analysis was performed using the software tool tRNAmodviz (http://genesilico.pl/trnamodviz). Distributions of nucleotide residues at each position of the hairpin (positions 27 to 43) are visualized as a pea chart, of which the color code is indicated at the top right corner of the figure. Universal numbering system is used (144). Acronyms of modified nucleotides present in certain isoacceptor tRNA are indicated outside the anticodon hairpin ('b' means bacteria, 'e', eukaryote and 'a', archaea). Those present in the proximal 'extended anticodon' (in green square box, see also text) are in grey background. Modified nucleotides present at positions 37 and 34 are listed in Figure 5A and B, respectively. The red arrow pointing from B34 of anticodon to B3 of the codon and the plateau of B35/B2 and the red arrow from B37 to the plateau of B36/B1 symbolize the stabilizing effects of B34/B37 on codon-anticodon pairings.
One remarkable exception is found in tRNA Ala (anticodon GGC) and tRNA Pro (anticodon GGG) with A32 instead of the quasi universal pyrimidine, but not in the other isoacceptor tRNA Ala (anticodon U*GC) and tRNA Pro (anticodons UGG and CGG, respectively) in E. coli (65,76). Base-38 is also occasionally a U or C instead of the usual A38, the uridine-38 being often modified to . This situation is found in majority of tRNA specific for Ala, Leu (of the 4-codon boxes), Val, His, Gln, Asp, and Glu, all belonging to decoding systems of the upper half of the codon wheel. In some instances, B32-B38 can even form a Watson-Crick pair, either U32-A38 (Pro, Gly, Thr) or A32-U38 (Ala, Pro). Again, these tRNAs with a potential B32-B38 Watson-Crick pair all belong to the group of strong 4-codons decoding tRNAs (Ala, Pro) and, surprisingly also to tRNA Thr belonging to the intermediate group.
The case of tRNA Ala has been studied in detail (36,(77)(78)(79). The take-home lesson is that tRNAs harboring too strong a base-pairing binding capacity have been 'tuned down' during evolution for optimal and uniform translation. As dis-cussed above, B32-B38 should not form a Watson-Crick base pair that antagonizes the canonical conformation of the anticodon loop with the anticodon triplet in a correct helical orientation. With complementary bases at B32 and B38, one expect that a Watson-Crick pair is formed in some tRNA population, leading to a less optimal anticodon conformation for pairing with the codon and, thus, possibly rejection from the ribosomal decoding site. Without a selection for the preformed conformation of the anticodon loop, the strong codon/anticodon pairs would allow miscoding by binding to other G/C-rich near-or even non-cognate aa-tRNAs.
The distributions of B37 in Figure 6C (including posttranscriptional modifications) follow each quadrant and thus, strongly depend on the codon base at the first position. This conclusion is valid for the other tRNA repertoires analyzed (Supplementary Figure S3). The presence of stabilizing modifications such as ct 6 A and ms 2 i 6 A in E. coli (or i 6   ing codon-anticodon base pair A1-U36 or U1-A36, respectively. These modifications, as well as the N1-methylation of G37, which is unable to base pair in a Watson-Crick mode, are also required for limiting ribosome frameshifting at specific AU-rich 'slippery' sequences of the mRNA (80,81).
In summary, purine-37, together with B31-B39 and B32/B38, and contacts with U33 (grey background in Figure 4), synergistically contribute to maintaining a conserved pre-formed and optimal anticodon conformation, 3stacked with both the rest of the anticodon arm and the last two bases B35 and B36 in the anticodon. Supplementary  Table S1 shows an approximate partition of these energetic elements that, together with the distribution of key elements on the wheels, helps explain how these apparently disparate interacting elements (including modified bases and ribose methylation) may cooperate to average the binding energies of tRNAs, whatever the GC or AU content of their anticodons to the mRNA-programmed ribosomes. One advantage of such an intricate networking system is the possibility to fine-tune the aminoacyl-tRNA for its ability to read the various cognate codons with appropriate efficiencies while avoiding reading near-and non-cognate codons.

CONTRIBUTION OF NUCLEOTIDE MODIFICATIONS AT THE FIRST tRNA ANTICODON POSITION
A large variety of post-transcriptionally modified nucleotides are found at position 34 of cytoplasmic elongator tRNAs (65,76). Partitioning of these nucleotides according to organisms in the three domains of life is given in Figure  5B. The most frequently and diversely modified nucleotide at position 34 is uracil, followed by cytidine and guanine. The only modified form of the adenine is hypoxanthine (inosine). Methylation of 2'O-ribose also occurs at G34, C34 and U34 (in the latter case, in combination with additional modifications at C5-and/or C2-atom of the base). These chemical adducts to B34 or its ribose alter the nucleotide chemical properties, and thus its polarity, but also the preference for the syn or anti conformer and the base's tautomeric form. They can also control ribose puckering and, in addition, allow for additional interactions with other elements of the translation machinery (see Supplementary Table S2 for an overview).
The distribution of modifications at B34 (C34, A34, G34) in elongator tRNAs according to the triplet code is depicted in Figure 7A in the case of E. coli. The same analysis for Haloferax volcanii, S. cerevisiae, Mycoplasma capricolum and human mitochondria are given in Supplementary Figure S4. Again, a clear segregation exists between tRNAs belonging to the 'strong binding' class (top part of the wheel) and those belonging to the 'weak binding' class (bottom part of the wheel). Except for inosine-34, which is found only in bacterial tRNA Arg of the 4-codon boxes, other modified nucleotides at position 34, like ac 4 C, k 2 C, Cm, Um, Q and gluQ are found only in tRNAs belonging to the 2:2 and single-codon boxes.
The situation concerning modified U*34 (*: modified at the C5-atom) and its doubly modified s 2 U* derivatives constitutes a special case. In the ribosomal grip, the third nucleotide of a codon is constrained by several contacts with the ribosome (specifically with G530, C518 of the 16S rRNA and a residue from RpS12), while the anticodon nucleotide B34 is in a loose stacking contact with C1054. When B34 forms a Watson-Crick pair, this asymmetry in contacts is not discriminatory. However, when forming a GoU wobble pair, the unsymmetrical constraints amplify the non-isostericity of UoG pairs such that G34oU3 wobble pairs can form but U34oG3 pairs form much less readily (for more detail see Supplementary Figure S5). Indeed, in several crystal structures with a U34oG3 opposition, the pair formed is not the standard wobble pair but a Watson-Crick-like pair between tautomeric forms of either G or U (27,(29)(30)(31)(32). The modifications at U34 are structurally imposed by the constraints on the codon base. Thus, in the non-isosteric GoU pairs, modifications on U34 will be necessary throughout the wheel, especially in the absence of additional C34-containing tRNAs. With C34-containing tRNAs present, the modifications at U34 are less important, since there is no need to form a U34oG3 pair. Despite this, when chemical adducts at the C5-atom of uracil are examined, a clear segregation by type exists depending on whether the modified U34-containing tRNAs belong to exclusively 4-codons (strong) or mainly 2-codons (weak) Without exceptions, aminoacyl-tRNAs reading codons of the split 2:2 or 3:1 decoding boxes, as well as those reading codons harboring A and/or U at the B1-B2 codon positions (bottom of the wheel) are the most diversely posttranscriptionally (hyper)-modified at B34, but also at other positions of the anticodon hairpin (compare Figure 6A-C). Evidently, fine-tuning the efficacy and accuracy of mRNA translation is a highly sophisticated process in which the chemistry of each nucleotide (modified or not) within the anticodon hairpin of the various individual natural elongator tRNAs collectively play a pivotal role (61).   (146)) is inserted between the circle for the third base and that for the amino acid type. Conventional one letter code for amino acid is used. For the sake of clarity, the distribution of G34, I34 and C34 derivatives are shown in (A), and U34 derivatives in (B). Modified G/C/I-34 are indicated in red. All Q-containing tRNAs and those containing modified C* are found in tRNA belonging to split 2:2-codon boxes. For the modified U34 (shown in (C)), two chemically distinguishable types of modified residues are found, one harboring an oxyacetic acid group (sometimes methylated) at the C5 atom of uracil (cmo 5 U or mcmo 5 U), are indicated as green. These U34 derivatives are found only in isoacceptor tRNA belonging to the unsplit 4-codon family boxes. The second type of modified U34 derivatives harbors a methylaminomethyl (sometimes carboxymethylated) group at the C5 atom of uracil (mnm 5 U or cmnm 5 U). They are indicated as blue, some of which are also hypermodified into 2-thiolated derivatives (s 2 U*) or methylated on the 2 -hydroxyl ribose (U*m). They are found in all split 2:2-codon boxes, and in the Arg/Gly 4-codon boxes (after (88)

SHORT SUMMARY OF THE STRUCTURAL ASPECTS
Overall, the three-dimensional structure of the anticodon loop is characterized by two helical stacks interrupted at the phosphate 3 of the invariant U33 residue: the 5 -stack in continuity with the 5 -strand of the anticodon stem stops at residue U33 (the U-turn), while the 3 -stack extends from residue B34 to B38 in continuity with B39 to B43 in the 3 -strand of the anticodon stem. Although residues B34 to B38 form a helical stack with standard rotation angles between nucleotides, the 3 -strand conformation does not lead to enough rotation to allow B38 to form a Watson-Crick pair with B32 (70). Indeed, formation of a Watson-Crick pair between B32 and B38 would require a change in the conformation of the rest of the loop. B39, however, has rotated enough and is able to form a Watson-Crick pair with B31. Evidently, various combinatorial structural solutions exist to generate a canonical architecture with a balance between flexibility and rigidity of the anticodon hairpin that allows the elongator tRNA to perform a common decoding function on the ribosome. This common functional architecture in elongator tRNAs depends very much on the conserved residue U33. Only in exceptional cases, as in many initiator tRNA Met in eukaryotes (not discussed here) or in the tRNA Ser of most genius of Candida species harboring a leucine anticodon CUG, have differences at this positions been found (C33 or G33, respectively) (65,76).
Thus, the deciphering of AU-rich codons within the decoding A site of the ribosome requires extra stabilization sources that of course have to arise mainly from the tRNA. Along the course of evolution, various strategies based, Figure 8. Modulation of codon-anticodon binding according to G+C (STRONG) or A+U (WEAK) binding capability of selected E. coli tRNA species. In both cases, the anticodon hairpin is schematically represented with all the nucleotides of the 5 branch in continuous stacking up to B34 and the complementary nucleotides of the 3 branch in continuous stacking up to B32. The U33 turn is indicated with its links to R35 and/or Y36 (underlined). At positions B32 and B38, various combinations of base opposition are found (boxed with dashed lines), the most frequent ones being indicated in bold letters. In red are the modified nucleotides; underlining emphasizes that the chemical adduct reinforces the stacking power of the base with the neighboring nucleotides. Modulation of the strength of codonanticodon binding occurs by the anticodon loop constraints that mainly depends on the choice of the B32-B38 base opposition, additional interactions with the conserved U33 and the identity of the chemical adducts on B37 and B34 that stabilize the B36-B1 and the B35-B2 interactions, respectively (schematized by red arrows). The number 2 with an asterisk for Pro, means that information came only from tDNA sequence, the corresponding maturated transcripts have not been sequenced yet. On the right of each diagram, an approximate energy scheme is displayed (in a blue rectangle) with S standing for STRONG, I for INTERMEDIATE, and W for WEAK. Same analysis for tRNAs from H. volcanii, S. cerevisiae, M. capricolum and human mitochondria are shown in Supplementary Figure  S6. as in any complex molecular system, on the interplay between many weak interactions together with a few strong interactions have offered alternatives to overcome this decoding problem. In the proximal 'extended anticodon', every single nucleotide between 31 to 39 in the anticodon loop contributes to uniform decoding by adapting its structure to the strength of the codon/anticodon pairings and the contacts within the ribosomal grip. Some combinations of B31-B39 and B32-B38 favor the canonical anticodon triplet conformation, contributing favorably to the free energy of binding, while others can be either too floppy or too rigid (with even the formation of an additional base pair at B32-B38), contributing unfavorably the free energy of binding. The stacking contacts of the modified base B37 sandwiched between B36 and B32-B38 can be extensive or minimal and contribute at the proper interface between the codon/anticodon triplet and the anticodon stem. Finally, the critical role of modifications of B34 for promoting nonstandard base pairs has been emphasized.

DEVIATIONS FROM THE STANDARD GENETIC CODE
Despite the various decoding strategies existing in Bacteria, Archaea, Eukarya and organelles, the meaning of each of the individual 64 triplets of the genetic code remains remarkably universal. However, the genetic code is not frozen and a few deviations have been identified in the nuclear genomes of a limited number of species (fungi, green algae, ciliates, diplomonads and firmicutes), though mostly in the small A+T rich genomes of mitochondria (metazoan, fungi, green plants and red algae) (1,82,83). Particular codon reassignments appeared to have taken place repeatedly and independently in the different lineages, attesting that certain codons are more prone than others to this evolutionary process. Also, two proteogenic amino acids, selenocysteine (Sec) and pyrrolysine (Pyl), have been added to the standard 20-member amino acid alphabet (natural expansion of the genetic code). These two amino acids are inserted in a few proteins in response to stop codons UGA and UAG, respectively, of cytoplasmic mRNAs of a few bacteria, eukaryotes and even archaea (84,85). Figure 9 summarizes the best-characterized genetic code deviations identified so far. The most frequently encountered deviations in both nuclear and mitochondrial genomes concern the switch of terminator codons to sense codons (UGA to Trp, Cys, Gly or Sec, UAG to Gln, Leu, Ala or Pyl, and UAA to Gln or Tyr). Identity changes of sense codons are less frequent and found mainly in mitochondrial genomes (except in one case). They are AUA-Ile to Met; AAR-Lys to Asn, AGA/G-Arg to Ser or Gly, all belonging to 2-codons family boxes for which the first base of the codon is either A1 or U1. The few cases where a C1-containing codon is involved are CUY/R-Leu to Thr or Ala in mitochondria of a few filamentous Saccharomycetacea and the CUG-Leu to Ser in the cytoplasm of Candida and Debaryomyces species. In these latter cases, the replacement of one amino acid by another results mainly from mischarging of a tRNA rather than of misreading of a codon during translation. Among the sense codons, those starting with G1 appear most resistant to amino acid reassignment. Also, as a direct consequence of strongly biased codon usage in mitochondria, a few codons are either not used (remain unassigned) or possibly used as alternative stop codons (2).
Evidently, the most frequent codon reassignments are found in A/U-rich weaker codon binders and correspond mostly to switch between two amino acids of the same pair of split codon family boxes. These codons are also the most liable to translational miscoding (see for example (86)(87)(88)). Similarly, the expansion of the genetic code to additional selenocysteine (Sec) and pyrrolysine (Pyl) occurs mainly, but not exclusively at stop codons, again at the same more frequently AU-rich reassigned codons (85,89).
Changing the meaning of a codon in a genome is challenging. It implies a few coordinated adjustments of various elements of the translation apparatus that necessarily result from a multistep evolutionary process. Among these steps are: (i) mutation in a gene coding for a tRNA (or a duplicated variant) to allow the emergence of a 'novel' tRNA able to efficiently read a given codon (or a set of near cognate codons) that is (are) distinct from the original Nucleic Acids Research, 2016, Vol. 44, No. 17 8031 Figure 9. Deviations from the standard, almost universal genetic code. The most frequent deviations concern terminator codons unexpectedly efficiently translated as sense codons for amino acids Trp, Gly, Cys, Gln, Tyr, Leu or Ala. Conventional one letter code for amino acid is used. In red are shown amino acid reassignments observed in nuclear genomes while in green are those related to mitochondrial genomes. The less frequently encountered reassignment of sense codons for a non-standard amino acid have been observed essentially in mitochondrial genome (indicated in green), except in one case (a Leucodon coding for Ser in Candida and Debaryomyces species (83). For more details and references of original papers (in addition to those cited in text), (1,82,147). The codons remaining unassigned and possibly playing the role of occasional alternative stop codon in certain organisms are indicated by a red (nuclear) or green (mitochondria) circle around B3, dotted line means avoided codons, plain lines means unassigned, possibly stop codons. The red arrows between 2 boxes indicate a switches between amino acids within 2:2 decoding boxes. Outside the circle, are indicated the reassignment of stop codon UGA and UAG into Sec (21st proteogenic amino acid) and Pyl (22nd amino acid), respectively (2). The most frequent codons reassignments mostly occur within the blue dotted boxes. wild type tRNA species; (ii) mutations in the gene coding for an aminoacyl-tRNA synthetase (or a duplicated variant) to allow the mutant enzyme to charge efficiently the newly evolved tRNA with an amino acid distinct to that expected from its anticodon; (iii) mutations in the gene coding for a modification enzyme (or a duplicated variant) to allow (or anneal) specific base modification that can influence the decoding property of the novel tRNA; (iv) mutations in the gene coding for potential competing tRNA (or even its elimination) that may read the same reassigned codon(s); (v) mutations in the termination factor (or even its elimination) to avoid competition with the novel functional suppressor tRNA (discussed in (90)); (vi) sequence rearrangement in mRNA for signaling with stop codon (or sense codon) has to be recoded (or reassigned); (vii) invention of a new elongation factor for selectively bringing the Sec-tRNA Sec into the decoding site of the translating ribosome; (viii) possibly mutations in the rRNA or r-proteins of the decoding site of the ribosome to allow optimal translation of the newly reassigned codon by its new cognate (or near cognate) aa-tRNA. Not all of these adjustments of the translation apparatus are simultaneously required. However, depending on the case, at least a combination of them has to occur. In this way, each time a reassigned codon, or a stop codon with the appropriate neighboring sequence is encountered in mRNA, the newly assigned amino acid is incorporated following standard decoding mechanism on the ribosome. This is at variance with the phenomena of translational missense error or stop-codon 'readthrough' that are mediated by a standard aa-tRNA aa harboring non-cognate anticodon or after occasional mischarging of a tRNA with a wrong amino acid, which in fine give rise to occasional misincorporation of an amino acid into the growing polypeptide. The yields of such naturally occurring mischarging and misreading are usually very low (below 10 −5 to 10 −3 ), unless, again, mischarging is genetically programmed by multiple sequential mutations in the tRNA, mRNA, various translational proteins and/or even the rRNA or proteins of the ribosome in order to eventually reach the degree of perfection which has been achieved for the reassignment or recoding of codons in a few organisms (91). Among all the above combinatory strategies for codon reassignments, the more intriguing ones are those involving posttranscriptional modification machineries. For example, in mitochondria of several metazoa and budding yeasts, C34 of the unique Met-tRNA Met is post-transcriptionally modified to 5-formyl-cytidine (f 5 C34) by a unique mitochondrial 5-formyl-cytosine synthase that allows the mature Met-tRNA Met to translate both codons AUG and AUA codons as Met (92,93). The completeness of such AUA-Ile-to-Met reassignment depends on the absence of a competing lysidine-containing Ile-tRNA Ile (anticodon k 2 CAU). However, in mitochondria of a totally different lineage (Ascidian Halocynthia roretzi), instead of a f 5 C34, a 5-taurine-2-thio-methyl-2-thiouridine ( m 5 s 2 U34) and an exceptionally unmodified U37 were found in Met-tRNA Met (94). In the mitochondria of the squid Loligo bleekeri, an unmodified C34 and a rather normal anticodon loop were found but differences in other parts of the Met-tRNA Met were required for reading the AUA codon as Met (95). Evidently, different independent solutions have been developed along mitochondrial evolution for solving the same problem. Several other examples of this sort have been reported for the reassignements AAA-Lys to Asn (dependence on a peculiar pseudouridine-35 in Asn-tRNA Asn ) and AGA/AGG-Arg to Ser or Gly (dependence on a unique m 7 G34 in Ser-tRNA Ser or m 5 U34 in a B34-mutant of the same Ser-tRNA Ser ) (1). Most of these reassignments depend on base modifications and are prevalent for codons belonging to the A1-quarter of the wheel code and, therefore, frequently depend on the decoding property of a U34-containing tRNAs, for which the panoply of chemical modifications is huge (see discussion above about the non-isostericity of UoG pairs and Figure 5B and Supplementary Figure S5).
The presence of nucleotide modifications, like m 6 A, m 1 A, I or in mRNAs can now be measured and assessed in translation (49,(96)(97)(98). The presence of m 6 A has moderate effects on the measured translation efficacy, the effect being strongest with m 6 A at the first position (49,99). The effects of O6-methylG depend strongly on the codon position, with incorporation at the second position B2 leading surprisingly to stalling of the ribosomes (98). O6-methylG, because of the absence of a hydrogen atom at N1 position of the base, cannot form a Watson-Crick-like base pair with C, but does so with U. Consequently, a decrease in accuracy at the first and third B1 and B3 positions was observed but, unexpectedly, not at the second position. Incorporation of in bacterial mRNAs led to a 30% repression of translation, with the strongest effect at the B3 position (49). However, in this respect, the conversion of nonsense into sense codons by the presence of in the mRNA is particularly striking (42,100) (Supplementary Figure S7 shows that the observed miscodings occur at the bottom of the wheel and especially in the U1 quarter). The miscoding effects of 35, in which the absence of a queuosine (Q34) promote misreading of G34oA3 or G34oG3 (101) and have been discussed above. In the published crystal structures, 1 leads to the formation of non-Watson-Crick cis Hoogsteen-Watson-Crick G2oA35 or A2oG35 pairs, with the mRNA base in the syn conformation (42). More crystal structures are necessary to complete our understanding of the molecular basis of these intriguing observations (102).
Because the activity of modification enzymes usually depends on available cofactors of the basic cellular metabolism, we anticipate that the efficiency of such type of codon reassignments depends on the physiological state of the cell and may constitute a regulatory device. They possibly correspond also to early evolutionary steps toward less dependent and more permanent amino acid reassignment.
Briefly, the most frequently encountered code deviations, recoding and codon reassignment, occur in the U1/A1 quadrants, while the G1 quadrant is still unaffected to date. The most frequent miscoding, i.e. decoding errors and code ambiguities also occur in the bottom half of the wheel (together with the Glu/Asp and His/Gln ambiguities) (86,87). Such errors are highly dependent on the tRNA modifications and constitute a rich source of post-transcriptional regulation. Recoding of stop codons (Sec and Pyl) require, not only new enzymes for the new metabolism (4 for Sec and 2 for Pyl), but also several new factors and cofactors and, especially structures in the mRNA for signaling a given stop codon for recoding. Codon reassignments require less specific changes in the translational apparatus (aminoacyl tRNA synthetases, absence of some release factors, etc.) but appear very efficient without competition and should probably be favored for further studies and applications.

THE EVOLUTION OF THE GENETIC CODE
If the organization of the genetic code is, as we suggest, anchored in the structural contacts with the associated energetics correctly integrated within the networks of interactions at the ribosomal decoding grip, it should also reflect historical aspects of the evolution of the genetic code. Nowadays, it is understood that there is not a single explanation for the emergence and expansion of the genetic code as it stands to date. The current most plausible hypothesis is, however, that the earliest proto-biosynthetic system originated from RNA:RNA duplexes with the most stable complementary GC-rich triplets coding for small polypeptides composed of alternate hydrophobic Ala and hydrophilic Gly, the very first amino acids encoded. As more adenine and uracil was introduced into the RNA, this minimalist proto-biosynthetic system later progressively evolved to a more systematic use of repeating GNC, GNY, GNS or SNS triplets (N, Y and S denote respectively: any of the four canonical nucleotides; C or U; and C or G) (103)(104)(105)(106). The earliest amino acids used to synthetize ancestral polypeptides were found in the prebiotic soup and selected according their specific interactions with the ancestral codons (the 'stereochemical hypothesis': (107,108)). These steps were subsequently expanded through the coevolution with the invention of biosynthetic pathways for new amino acids (the 'amino acid metabolism hypothesis', initially pointed out by: (109); reviewed and updated: by (110,111)) and the emergence of the corresponding primordial aminoacyl-tRNA synthetases able to fix these new amino acids on appropriated proto-tRNAs ('co-evolution with tRNA aminoacylation systems': (112)(113)(114)). Constant refinement at both the replication and translation levels allow to progressively minimize the impact of coding errors and to increase the diversity and functionality of proteins that can be made with a larger amino acid alphabet (the 'error minimizing code hypothesis': (107,115,116)). Finally, as mentioned in the previous section, the code can further evolve by reassignment of unused, temporarily ambiguous or less used codons for another canonical or even totally new amino acids (the 'codon capture theory': (117,118)). Finally, early horizontal transfer and collective evolution of the code through different subspecies has been emphasized (see for examples (119,120)). In other words, the present day genetic code did not necessarily result solely from divergent evolution, but also from collective evolution via the development of an innovation-sharing process that allows the emergence of a quasi universal genetic code among populations of species 'speaking the same language' (121).
The take-home lesson of all the above information is that along the expansion of the genetic code, an optimal stability of complementary codon-anticodon pairs appears to have been the main evolutionary force. Starting with a G/C-rich RNA coding for simple abiotic amino acids like Gly, Ala, Pro and probably a precursor of Arg allows the subsequent development of more complex modern proteosynthetic machineries as we know to date. As argued by Trifonov (122), the temporal order of the sequential acquisition of new amino acids into the code seems to result from a 'descending' thermostability of the respective codon:anticodon pairs, with initiation AUG codon among the latest codons to enter the code, together with UAG/UAA/UGA codons which remain unassigned and serve as termination codons. As explained above, this 'descending' thermostability was obviously compensated along evolution by a better organization of the tRNA anticodon hairpin and of the ribosomal grip within the ribosome. The main differences in reading the 'ancient sense strong codons' of the 4-codon boxes as compared to the subsequently acquired 'new, weak, sense codons' of the 2:2codon and 3:1-codon family boxes is that the former are less dependent than the latter on stabilizing modified nucleotides within the proximal extended anticodon ( Figures  6 and 7).
The progressive introduction of A and U into the primordial G/C-rich RNA (both in proto-mRNA and proto-tRNA) was of great advantage for expanding the decoding potentiality of the proteosynthesis machinery ( Figure  10A). However, as mentioned in a preceding section, the use of U-rich codons and U34-containing tRNA is promiscuous, both for miscoding and frameshifting of the commafree mRNA during translation. It was possible to solve this dilemma only by avoiding usage of such U-rich 'slippery' codons and by developing specific enzymatic systems for posttranscriptional modifications of base of the anticodon hairpin, mainly at B34 (Figure 7) and B37 ( Figure  6C) (64,(123)(124)(125). Because frequent frameshift events are probably more deleterious to the cell than frequent miscoding of certain codons, it is possible that emergence of B37 modifications systems preceded modifications at B34 (126)(127)(128). The advantage of restricting the coding potentiality of certain tRNA by modifications of B34 is that it allows the introduction of new amino acids into the code, which, however, require duplication of a tRNA gene and of the emergence of the corresponding aminoacyl-tRNA synthetases (co-evolution of tRNA genes, tRNA maturation enzymes and tRNA aminoacylation appara-tus). The earliest of these modification systems would have appeared before the divergence of Woese's three major domains (129), perhaps very early or at least before completion of the genetic code. These modifications would be universally distributed today (130), while those more recently invented occur either specifically in archaea, bacteria and/or eukaryotes ( Figure 5A-B, discussed in (64)). Lateral transfer of genes coding for a given modified enzymes present in bacteria, organelles and/or the cytoplasm of eukaryotes could also have occurred through endosymbiosis or virus infection. Remarkably, the chemical adducts of several modified nucleotides (B34 and B37) contain an amino acid, as in N6-threonylcarbamoyladenosine (t 6 A37), N6-hydroxynorvalinecarbamoyladenosine (hn 6 A37), N6glycinecarbamoyl adenosine (g 6 A37), 5-taurinomethyl-2thiouridine (tm 5 s 2 U34) and 3-(3-amino-3 carboxypropyl)uridine (acp 3 U47), lysidine (k 2 C34), agmatidine (agm 2 C34) and glutamyl-queuosine (gluQ34). These modifications might be relics of ancient code, when the RNA hairpin structures harboring the primordial anticodons (the ancestor of tRNA) were strongly associated, and possibly even 'charged' with amino acids (131)(132)(133).
As a rule, codons ending with pyrimidines U3 or C3 are rarely distinguishable in translation and always behave as synonymous codons, limiting the number of tRNAs able to read four codons to a maximum of three species ( Figure  10A). Under special genomic selection conditions, such as strong A/T or G/C pressure, whether coupled or not with genome size reduction, the proteosynthetic machinery can further evolve toward a simpler version than the one found in majority of extant biological cells. In the parasitic Mycoplasmas like M. capricolum and more evidently in mitochondria, both with A/T-rich genome ( Figure 10B), fewer isodecoder tRNAs (28 and 35, respectively, instead of 40-45 in most cells) are needed to translate mRNA coding for 20 canonical amino acids. The reduction in the number of distinct tRNAs concerns mainly C34-containing tRNAs in 4-codon boxes and eventually 2:2 decoding boxes with the consequence that now both Met and Trp become encoded by 2 synonymous codons as in mitochondria. More extensive tRNA reduction can occur with A34/G34-containing tRNA as long as the uracil of the remaining U34-containing tRNA is not post-transcriptionally modified (loss of corresponding maturation enzymes). Accordingly, the codon usage in these cells or organelles is strongly biased, with avoidance of NNG codons (see Supplementary Figure S4 and (59,60,117)). Conversely in some Micrococcus strains like M. luteus, with a G/C-rich genome ( Figure 10B), reduction of tRNA applies only to U34-containing tRNAs and the usage of codons is almost exclusively NNC and NNG. The main difference between M. luteus ( Figure 10B) and primordial cells is that 20 amino acids, and thus 2:2 codon boxes are used, while in the ancestral cells only a few amino acids, corresponding to only a few G/C-containing codons, prevailed.
Briefly, after an initial cooptation of a few abiotic amino acids into the code, further evolution of the primordial RNA with amino acids occurred through intricate but interdependent stepwise expansion of the amino acid biosynthetic pathways, coupled with the evolution of aminoacyl-tRNA synthetases and tRNA modification enzymes. As a 8034 Nucleic Acids Research, 2016, Vol. 44, No. 17 Figure 10. Hypothetical stepwise evolution of the genetic code and translation machinery in Bacteria, Eukarya, Archaea and mitochondria. (A) Starting from highly G+C rich small RNA pieces and a few abiotic amino acids, the primordial translation system evolved by extending its coding capacity by result, codon:anticodon pairs with intrinsically weaker interactions provided greater diversity in amino acids (chemical structures and global polarities). These pairs also concentrate a large number of deviations from the genetic code or codon reassignments, and possibly correspond to best candidates for further insertions of additional noncanonical amino acids of special technical interest. The elucidation of mechanisms and the development of tools for the recoding of sense codons are currently actively being researched in order to be able to produce proteins with a variety of non-natural amino acids at several positions (134,135). The rules delineated here should, with further work, help in this development.

CONCLUSIONS
Cells from different organisms, or the same organism but from different tissues, generally have different population of functional tRNAs. Both the total number of tRNA species with distinct anticodons (tRNA repertoire or tRNome) and the relative amounts of each individual isoacceptor species (isodecoders) differ much from one type of cell to another (69,136). Similarly, the frequency with which each codon is translated (codon usage or codon preference) varies depending on the origin of cell and on the levels of gene expression. This is due partly to varying amino acid composition of the proteins synthesized and most importantly to the different use of synonymous codons (codon bias) that is dependent on the global G+C content of the mRNAs (125,(137)(138)(139)(140)(141)(142).
Inspection and comparison of the mappings of the extended anticodon concept wheels illustrates how, despite the constraining process of decoding and the use of an identical genetic code (except for Met and Trp in certain types of cells like mycoplasma and mitochondria), a large variety of tRNA repertoires and codon usage is still allowed and prevalent. The interpretation of this biological diversity in molecular terms is a daunting task. However, it should be clear that evaluating the frequency of a given codon in individual mRNAs or, as is usually done for statistical reasons, of bulk cellular mRNAs, is not sufficient, and that the exact composition of the tRNA repertoire in a given cell or organism (see below) has also to be determined. Complications arise from the many diverse underlying evolving molecular mechanisms that preside over the biogenesis, stability and functional regulation of each component (tRNAs, mR-NAs, snoRNA in Eukarya and Archaea, possibly regulatory RNAs) of the translation machinery.
The majority of individual codons in mRNAs can be translated by either a cognate aa-tRNA harboring the three complementary base pairs or by one or several near-cognate aa-tRNA carrying the same amino acid differing by the base-34 or any other important base of the extended anticodon region. The main advantage of representing data related to such complex decoding system in the form of a wheel of genetic code is that all sequence co-variances can be visualized within an internal energy-based logic, for each anticodon of the 61 almost universally used sense codons in a given organism. The concomitant emergence of machineries allowing post-transcriptional enzymatic modifications of bases and ribose at position 34 of especially tRNA species of the split 2:2-and 3:1-codon boxes also play a determinant role, not only to help codon-anticodon stabilization but also to guarantee a better discrimination between the two sets of Y-ending and R-ending codons. Furthermore, for smooth and regular decoding and ribosomal progression along the correct frame of mRNAs, optimal and uniform binding of the various aminoacyl-tRNAs independently of codon types is required ( Figure 11).
The codon usage, the number of tRNAs in each isoacceptor families and the numbers of isodecoder tRNAs all vary with cell type and cellular differentiation state but, despite the overwhelming complexity and potential ambiguity they introduce, as an ensemble they still have to comply with the networks of interactions which act at the decoding A site ( Figure 12).  (122)). The final decoding machinery as we know to date (and illustrated on the left panel of the figure) results from a long and complex intricate stepwise coevolution with the emergence of metabolically generated amino acids (up to 20), duplication and speciation of proto-tRNAs with distinct anticodons (up to 40-45 to date), post-transcriptional modification enzymes (up to more than 100 known to date), new amino acid tRNA synthetases (up to 20), complexification of the ribosomal architecture (rRNAs and r-proteins) and the introduction of additional protein factors allowing the extension and ultimate tuning the efficacy and accuracy of the decoding capability of 61/62 sense (cognate and near cognate) codons with in addition 2 to 3 terminators for 20 proteogenomic natural amino acids. Pyrrolysine and selenocysteine have been excluded from our analysis, as well as the situation of the special tRNA Met involved in the initiation of protein synthesis that obviously arose later during cell evolution. Emphasis is given to the importance of B34 modifications (indicated in red) that allows the segregation of the A/U-rich 4-codon boxes into split decoding boxes 2:2 and 3:1 with subsequent additional amino acids to enter the code. (B) Following genomic selection conditions such as directional constraints or mutational pressure on codons (strong G/C as in M. luteus (148) or strong A/T combined with drastic genome size reduction as in mammalian mitochondria (57,60) or the minimalist bacteria M. capricolum (58,59,149), simplifications of the usual translational decoding system are evident, while preserving the split decoding boxes because of the need to encode for 20 amino acids. Less tRNA species with distinct anticodons (isodecoders) are required (22 for human mitochondria, 28 in M. capricolum and 29 in M. luteus, again not including tRNA Met initiator). However, remarkably the only remaining B34 modifications found in natural tRNAs are those related to the split 2:2 and/or 3:1 decoding boxes. More information in relation to the corresponding tRNA modification enzymes of a given tRNA repertoire in the subgroup of Mollicutes (mainly mycoplasmas) can be found in (150). The situation of M. luteus is remarkable for its total absence of U34-containing tRNAs while in mycoplasma and mitochondria U34-containing tRNA are critical. Codon usage correlates with tRNA repertoire and amino acid type (125,(137)(138)(139)(140)(141)(142). Symbols for each acronym of the most important modified nucleotides are indicated within the figure. In parenthesis b, e, a, m means bacteria, eukaryotes, archaea and mitochondria respectively. More details about their chemical structures are shows in Figure 7 and Supplementary Figure S4. Within each decoding boxes, blue arrows correspond to the most dominant codon:anticodon pairs. Only B34 are indicated. When symbols of B34 modification are in parenthesis, it means that the unmodified version of B34 is used only in certain decoding boxes, while the modified version(s) is (are) used in other decoding boxes (for details see in Figure 7 for E. coli and Supplementary Figure S4 for H. volcanii, S. cerevisiae, M. capricolum and human mitochondria, respectively). When B3 is in bold, it means that the codon usage corresponding to the particular codon is dominant over the other near cognate codons of the same 4-or 2decoding boxes, while B34 indicated in regular italics correspond to rare codons (for details see as above, Figure 7 for E. coli and Supplementary Figure  S4 for the other organisms analyzed). Figure 11. Summary scheme illustrating the mapping of the energetics and evolutionary history on the wheel organization of the genetic code. Right vertical arrows: from bottom to top, there is an increase in the strengths of the networking interactions that is coupled with an increase of base modifications at B34 and B37 from top to bottom. Left vertical arrows at the right and left sides: the evolution from the primordial G/C-rich to A/U-rich codon/anticodons triplets required base modifications and the related metabolic enzymatic activities. an approximate and simplified energy scheme at the A site decoding center illustrates how the favorable and costly contributions top the free energies and compensate to maintain a smooth and regular translation process with minor final variations in free energy. On the wheel representation, the energies at the left can be mapped. However, not all energies can be mapped; for example, the interactions between the ribosome and parts other than the anticodon hairpin, the conformational distortions or alternative states of tRNAs, and the energies associated with ribosomal movements.