Consequences of stop codon reassignment on protein evolution in ciliates with alternative genetic codes

Tetrahymena thermophila and Paramecium tetraurelia are ciliates that reassign TAA and TAG from stop codons to glutamine codons. Because of the lack of full genome sequences, few studies have concentrated on analyzing the effects of codon reassignment in protein evolution. We used the recently sequenced genome of these species to analyze the patterns of amino acid substitution in ciliates that reassign the code. We show that, as expected, the codon reassignment has a large impact on amino acid substitutions in closely related proteins; however, contrary to expectations, these effects also hold for very diverged proteins. Previous studies have used amino acid substitution data to calculate the minimization of the genetic code; our results show that because of the lasting influence of the code in the patterns of substitution, such studies are tautological. These different substitution patterns might affect alignment of ciliate proteins, as alignment programs use scoring matrices based on substitution patterns of organisms that use the standard code. We also show that glutamine is used more frequently in ciliates than in other species, as often as expected based on the presence of the 2 new reassigned codons, indicating that the frequencies of amino acids in proteomes is mostly determined by neutral processes based on their number of codons.     

Evolution of the Genetic Code by Incorporation of Amino Acids that Improved or Changed Protein Function

Fifty years have passed since the genetic code was deciphered, but how the genetic code came into being has not been satisfactorily addressed. It is now widely accepted that the earliest genetic code did not encode all 20 amino acids found in the universal genetic code as some amino acids have complex biosynthetic pathways and likely were not available from the environment. Therefore, the genetic code evolved as pathways for synthesis of new amino acids became available. One hypothesis proposes that early in the evolution of the genetic code four amino acids—valine, alanine, aspartic acid, and glycine—were coded by GNC codons (N = any base) with the remaining codons being nonsense codons. The other sixteen amino acids were subsequently added to the genetic code by changing nonsense codons into sense codons for these amino acids. Improvement in protein function is presumed to be the driving force behind the evolution of the code, but how improved function was achieved by adding amino acids has not been examined. Based on an analysis of amino acid function in proteins, an evolutionary mechanism for expansion of the genetic code is described in which individual coded amino acids were replaced by new amino acids that used nonsense codons differing by one base change from the sense codons previously used. The improved or altered protein function afforded by the changes in amino acid function provided the selective advantage underlying the expansion of the genetic code. Analysis of amino acid properties and functions explains why amino acids are found in their respective positions in the genetic code.     

Trend of amino acid composition of proteins of different taxa

Archaea, bacteria and eukaryotes represent the main kingdoms of life. Is there any trend for amino acid compositions of proteins found in full genomes of species of different kingdoms? What is the percentage of totally unstructured proteins in various proteomes? We obtained amino acid frequencies for different taxa using 195 known proteomes and all annotated sequences from the Swiss-Prot data base. Investigation of the two data bases (proteomes and Swiss-Prot) shows that the amino acid compositions of proteins differ substantially for different kingdoms of life, and this difference is larger between different proteomes than between different kingdoms of life. Our data demonstrate that there is a surprisingly small selection for the amino acid composition of proteins for higher organisms (eukaryotes) and their viruses in comparison with the "random" frequency following from a uniform usage of codons of the universal genetic code. On the contrary, lower organisms (bacteria and especially archaea) demonstrate an enhanced selection of amino acids. Moreover, according to our estimates, 12%, 3% and 2% of the proteins in eukaryotic, bacterial and archaean proteomes are totally disordered, and long (> 41 residues) disordered segments are found to occur in 16% of arhaean, 20% of eubacterial and 43% of eukaryotic proteins for 19 archaean, 159 bacterial and 17 eukaryotic proteomes, respectively. A correlation between amino acid compositions of proteins of various taxa, show that the highest correlation is observed between eukaryotes and their viruses (the correlation coefficient is 0.98), and bacteria and their viruses (the correlation coefficient is 0.96), while correlation between eukaryotes and archaea is 0.85 only.     

Unambiguous correspondence and complementarity of the chemical properties of amino acids and anticodons of the genetic code

The chemical language of genetic code is proposed. As a result of chemical language application for the analysis of the modern genetic code, the existence of an unambiguous correspondence between the chemical properties of amino acids and their coding triplets (codons and anticodons) is shown. This confirms the hypothesis of the code chemical determination. The complementarity between the chemical properties of amino acids and their anticodons (but not the codons) has been found also to exist. This observation supports the hypothesis of the genetic code determination by the direct recognition and also underlines the primary role of anticodon in the origin of genetic code in comparison with codons.     

Structuring of the genetic code took place at acidic pH

I have observed that in multiple regression the number of codons specifying amino acids in the genetic code is positively correlated with the isoelectric point of amino acids and their molecular weight. Therefore basic amino acids are, on average, codified in the genetic code by a larger number of codons, which seems to imply that the genetic code originated in an acidic 'intracellular' environment. Moreover, I compare the proteins from Picrophilus torridus and Thermoplasma volcanium, which have different intracellular pH and I define the ranks of acidophily for the amino acids. A simple index of acidophily (AI), which can be easily obtained from acidophily ranks, can be associated to any protein and, therefore, can also be associated to the genetic code if the number of synonymous codons attributed to the amino acids in the code is assumed to be the frequency with which the amino acids appeared in ancestral proteins. Finally, the sampling of the variable AI among organisms having an intracellular pH less than or equal to 6.6 and those having a non-acidic intracellular pH leads to the conclusion that the value of the genetic code's AI is not typical of proteins of the latter organisms. As the genetic code's AI value is also statistically not different from that of proteins of the organisms having an acidic intracellular pH, this supports the hypothesis that the structuring of the genetic code took place in acidic pH conditions.     

The genetic code--40 years on

The genetic code discovered 40 years ago, consists of 64 triplets (codons) of nucleotides. The genetic code is almost universal. The same codons are assigned to the same amino acids and to the same START and STOP signals in the vast majority of genes in animals, plants, and microorganisms. Each codon encodes for one of the 20 amino acids used in the synthesis of proteins. That produces some redundancy in the code and most of the amino acids being encoded by more than one codon. The two cases have been found where selenocysteine or pyrrolysine, that are not one of the standard 20 is inserted by a tRNA into the growing polypeptide.     

Arginine as an evolutionary intruder into protein synthesis

The arginine content of proteins is, on the average, far less than would be anticipated from the fact that 6 of the 61 codons for amino acids in the genetic code are for arginine. In contrast, lysine is more abundant than would be expected from its 2 codons. It is suggested that arginine replaced ornithine in protein synthesis during the evolution of the genetic code.     

A methanogen hosted the origin of the genetic code

A comparison is made between orthologous proteins from a methanogen (Methanopyrus kandleri) and from a non-methanogen (Pyrococcus abyssi) in order to determine the amino acid substitution pattern. This analysis makes it possible to establish which amino acids are significantly and asymmetrically utilised by these two organisms. A methanophily index (MI) based on this asymmetry makes it possible for any protein to be associated with a numerical value which, when calculated for the same orthologous protein from methanogenic and non-methanogenic organisms, turns out to have the power to discriminate between these two groups of organisms, even if only for about 20% of the analysed proteins. The MI can also be associated to the genetic code under the assumption that the frequency of synonymous codons specifying the amino acids in the genetic code also reflects the frequency with which amino acids appeared in ancestral proteins. Finally a t-test shows that the MI value associated to the genetic code is not different from the mean value of the MI deriving from methanogen proteins, but it differs from the mean MI of non-methanogen proteins. This might indicate that the genetic code evolved in a methanogenic ‘organism’.     

The late stage of genetic code structuring took place at a high temperature

The correlation between the optimal growth temperature of organisms and a thermophily index based on the propensity of amino acids to enter more frequently into (hyper)thermophile proteins is used to conduct an analysis aiming to establish whether genetic code structuring took place at a low or a high temperature. If the number of codons attributed to the various amino acids in the genetic code constitutes an estimate of the mean amino acid composition of proteins produced when the genetic code was definitively structured, then the thermophily index can also be associated to the genetic code. This value and the sampling of the variable thermophily index of different alignments of protein sequences from mesophile, thermophile and hyperthermophile species make it possible to establish, with an extremely high statistical confidence, that the late stage of genetic code structuring took place in a hyperthermophile (or thermophile) 'organism'. Moreover the 95% confidence interval of the temperature at which the genetic code was fixed turned out to be 91+/-24 degrees C. These observations seem to support the hypothesis that the origin of life might have taken place at a high temperature.     

Noncanonical amino acid mutagenesis in response to recoding signal-enhanced quadruplet codons

While amber suppression is the most common approach to introduce noncanonical amino acids into proteins in live cells, quadruplet codon decoding has potential to enable a greatly expanded genetic code with up to 256 new codons for protein biosynthesis. Since triplet codons are the predominant form of genetic code in nature, quadruplet codon decoding often displays limited efficiency. In this work, we exploited a new approach to significantly improve quadruplet UAGN and AGGN (N = A, U, G, C) codon decoding efficiency by using recoding signals imbedded in mRNA. With representative recoding signals, the expression level of mutant proteins containing UAGN and AGGN codons reached 48% and 98% of that of the wild-type protein, respectively. Furthermore, this strategy mitigates a common concern of reading-through endogenous stop codons with amber suppression-based system. Since synthetic recoding signals are rarely found near the endogenous UAGN and AGGN sequences, a low level of undesirable suppression is expected. Our strategy will greatly enhance the utility of noncanonical amino acid mutagenesis in live-cell studies.     

Rewiring translation - Genetic code expansion and its applications

With few minor variations, the genetic code is universal to all forms of life on our planet. It is difficult to imagine that one day organisms might exist that use an entirely different code to translate the information of their genome. Recent developments in the field of synthetic biology, however, have opened the gate to their creation. The genetic code of several organisms has been expanded by the heterologous expression of evolved aminoacyl-tRNA synthetase/tRNA(CUA) pairs that mediate the incorporation of unnatural amino acids in response to amber codons. These UAAs introduce exciting new features into proteins, such as spectroscopic probes, UV-inducible crosslinkers, and functional groups for bioorthogonal conjugations or posttranslational modifications. Orthogonal ribosomes provide a parallel translational machinery in Escherichia coli that has lost its evolutionary constraints. Evolved variants of these ribosomes translate amber or quadruplet codons with massively enhanced efficiency. Here, I review these recent developments emphasizing their tremendous potential to facilitate biochemical and cell biological studies.     

Code dependent conservation of the physico-chemical properties in amino acid substitutions

The frequency of amino acid replacements in families of typical proteins has been elegantly analyzed by Argyle (1980) showing that the most frequent replacements involve a conservation of the amino acid chemical properties. The cyclic arrangement of the twenty amino acids resulting from the most frequent replacements has been described as an amino acid chemical ring. In this work, a novel amino acid replacement frequency ring is proposed, for which a conservation of over 90% of the most general physico-chemical properties can be deduced. The amino acid chemical similarity ring is also analyzed in terms of the genetic code base probability changes, showing that the discrepancy that exists between the standard deviation value of the amino acid replacement frequency matrix and its respective ideal value is almost equal to that deduced from the corresponding base codon replacement probability matrices. These differences are finally evaluated and discussed in terms of the restrictions imposed by the structure of the genetic code and the physico-chemical dissimilarities between some codons of amino acids which are chemically similar.     

Persistent biases in the amino acid composition of prokaryotic proteins

Correspondence analysis of 28 proteomes selected to span the entire realm of prokaryotes revealed universal biases in the proteins' amino acid distribution. Integral Inner Membrane Proteins always form an individual cluster, which can then be used to predict protein localisation in unknown proteomes, independently of the organism's biotope or kingdom. Orphan proteins are consistently rich in aromatic residues. Another bias is also ubiquitous: the amino acid composition is driven by the G + C content of the first codon position. An unexpected bias is driven, in many proteomes, by the AAN box of the genetic code, suggesting some functional biochemical relationship between asparagine and lysine. Less-significant biases are driven by the rare amino acids, cysteine and tryptophan. Some allow identification of species-specific functions or localisation such as surface or exported proteins. Errors in genome annotations are also revealed by correspondence analysis, making it useful for quality control and correction.     

线粒体遗传密码及基因组遗传密码的对称分析

病毒、细菌和真核生物的氨基酸编码都使用相同的遗传密码,表明它们可能有共同的来源。但人和牛的线粒体的遗传密码和基因组的遗传密码相比,出现以下不同;（1）ATA编码甲硫氮酸M而不是异亮氨酸I。（2）TGA不再是终止密码子X而编码色氨酸W。（3）AGA和AGG不再是精氨酸R的密码子而变为终止密码子X。应用高维空间拓扑分析的方法,对线粒体遗传密码和基因组遗传密码的6维编码空间进行对称性分析,得到如下结果：（1）线粒体遗传密码的起始密码子是2个而不是1个。（2）线粒体遗传密码的终止密码子是4个而不是3个。（3）线粒体遗传密码空间只有2、4、6三种偶数简并度而没1、3两种奇数简并度,表明其对称度较高。（4）线粒体遗传密码空间除丝氨酸S分成两个平行的子空间之外,终止密码子X亦分成两个平行的子空间,表明其连通度较低。（5）线粒体遗传密码一基因组遗传密码相比,共有3个简并平面出现变异,即：1001λλ（M和I）,011λ1λ（W和X）,以及1011λλ（S和X或S和R）。（6）基因组遗传密码的1、3两种奇数简并度可能来源于线粒体遗传密码的1001λλ平面和011λ1λ平面的对称性破缺。对线粒体遗传密码变异的生物学意义及遗传密码的起源进行了分析和讨论。     

Protein evolution drives the evolution of the genetic code and vice versa

A model for the developmental pathway of the genetic code, grounded on group theory and the thermodynamics of codon-anticodon interaction is presented. At variance with previous models, it takes into account not only the optimization with respect to amino acid attributes but, also physicochemical constraints and initial conditions. A 'simple-first' rule is introduced after ranking the amino acids with respect to two current measures of chemical complexity. It is shown that a primeval code of only seven amino acids is enough to build functional proteins. It is assumed that these proteins drive the further expansion of the code. The proposed primeval code is compared with surrogate codes randomly generated and with another proposal for primeval code found in the literature. The departures from the 'universal' code, observed in many organisms and cellular compartments, fit naturally in the proposed evolutionary scheme. A strong correlation is found between, on one side, the two classes of aminoacyl-tRNA synthetases, and on the other, the amino acids grouped by end-atom-type and by codon type. An inverse of Davydov's rules, to associate the amino acid end atoms (O/N and non-O/non-N) of 18 amino acids with codons containing a weak base (A/U), extended to the 20 amino acids, is derived.     

Codon Usage Bias and Mutation Constraints Reduce the Level of ErrorMinimization of the Genetic Code

Studies on the origin of the genetic code compare measures of the degree of error minimization of the standard code with measures produced by random variant codes but do not take into account codon usage, which was probably highly biased during the origin of the code. Codon usage bias could play an important role in the minimization of the chemical distances between amino acids because the importance of errors depends also on the frequency of the different codons. Here I show that when codon usage is taken into account, the degree of error minimization of the standard code may be dramatically reduced, and shifting to alternative codes often increases the degree of error minimization. This is especially true with a high CG content, which was probably the case during the origin of the code. I also show that the frequency of codes that perform better than the standard code, in terms of relative efficiency, is much higher in the neighborhood of the standard code itself, even when not considering codon usage bias; therefore alternative codes that differ only slightly from the standard code are more likely to evolve than some previous analyses suggested. My conclusions are that the standard genetic code is far from being an optimum with respect to error minimization and must have arisen for reasons other than error minimization.[Reviewing Editor: Martin Kreitman]     

Codon Size Reduction as the Origin of the Triplet Genetic Code

The genetic code appears to be optimized in its robustness to missense errors and frameshift errors. In addition, the genetic code is near-optimal in terms of its ability to carry information in addition to the sequences of encoded proteins. As evolution has no foresight, optimality of the modern genetic code suggests that it evolved from less optimal code variants. The length of codons in the genetic code is also optimal, as three is the minimal nucleotide combination that can encode the twenty standard amino acids. The apparent impossibility of transitions between codon sizes in a discontinuous manner during evolution has resulted in an unbending view that the genetic code was always triplet. Yet, recent experimental evidence on quadruplet decoding, as well as the discovery of organisms with ambiguous and dual decoding, suggest that the possibility of the evolution of triplet decoding from living systems with non-triplet decoding merits reconsideration and further exploration. To explore this possibility we designed a mathematical model of the evolution of primitive digital coding systems which can decode nucleotide sequences into protein sequences. These coding systems can evolve their nucleotide sequences via genetic events of Darwinian evolution, such as point-mutations. The replication rates of such coding systems depend on the accuracy of the generated protein sequences. Computer simulations based on our model show that decoding systems with codons of length greater than three spontaneously evolve into predominantly triplet decoding systems. Our findings suggest a plausible scenario for the evolution of the triplet genetic code in a continuous manner. This scenario suggests an explanation of how protein synthesis could be accomplished by means of long RNA-RNA interactions prior to the emergence of the complex decoding machinery, such as the ribosome, that is required for stabilization and discrimination of otherwise weak triplet codon-anticodon interactions.     

On the possible origin and evolution of the genetic code

The genetic code is examined for indications of possible preceding codes that existed during early evolution. Eight of the 20 amino acids are coded by ‘quartets’ of codons with four-fold degeneracy, and 16 such quartets can exist, so that an earlier code could have provided for 15 or 16 amino acids, rather than 20. If two-fold degeneracy is postulated for the first position of the codon, there could have been 10 amino acids in the code. It is speculated that these may have been phenylalanine, valine, proline, alanine, histidine, glutamine, glutamic acid, aspartic acid, cysteine and glycine. There is a notable deficiency of arginine in proteins, despite the fact that it has six codons. Simultaneously, there is more lysine in proteins than would be expected from its two codons, if the four bases in mRNA are equiprobable and are arranged randomly. It is speculated that arginine is an ‘intruder’ into the genetic code, and that it may have displaced another amino acid such as ornithine, or may even have displaced lysine from some of its previous codon assignments. As a result, natural selection has favored lysine against the fact that it has only two codons. The introduction of tRNA into protein synthesis may have been a cataclysmic and comparatively sudden event, since duplication of tRNA takes place readily, and point mutations could rapidly differentiate members of the family of duplicates from each. Two tRNAs for different amino acids may have a common ancestor that existed more recently than the separation of the prokaryotes and eukaryotes. This is shown by homology of twoE. coli tRNAs for glycine and valine, and two yeast tRNAs for arginine and lysine.     

Block structure and stability of the genetic code

It is known that different codons may be unified into larger groups related to the hierarchical structure, approximate hidden symmetries, and evolutionary origin of the universal genetic code. Using a simplified evolutionary motivated two-letter version of genetic code, the general principles of the most stable coding are discussed. By the complete enumeration in such a reduced code it is strictly proved that the maximum stability with respect to point mutations and shifts in the reading frame needs the fixation of the middle letters within codons in groups with different physico-chemical properties, thus, explaining a key feature of the universal genetic code. The translational stability of the genetic code is studied by the mapping of code onto de Bruijn graph providing both the compact visual representation of mutual relationships between different codons as well as between codons and protein coding DNA sequence and a powerful tool for the investigation of stability of protein coding. Then, the results are extended to four-letter codes. As is shown, the universal genetic code obeys mainly the principles of optimal coding. These results demonstrate the hierarchical character of optimization of universal genetic code with strictly optimal coding being evolved at the earliest stages of molecular evolution. Finally, the universal genetic code is compared with the other natural variants of genetic codes.     

The rearranged genetic code and its implications in evolution and biochemistry

We rearrange the genetic code and present a table of codons. The chemical properties of the amino acids coded by codons, and the evolutionary trend of codons are well reflected in the order of this table, from which two rules can be drawn: (1) the polarity/non-polarity and hydrophilicity/hydrophobicity of amino acids coded for by codons alternate row by row in the table; (2) in general, the lower down in the table, the earlier the codons are in terms of evolution.