The universal ancestor was a thermophile or a hyperthermophile.

By exploiting the correlation between the optimal growth temperature of organisms and a thermophily index based on the propensity of amino acids to enter thermophile/hyperthermophile proteins, an analysis is conducted in order to establish whether the last universal common ancestor (LUCA) was a mesophile or a (hyper)thermophile. This objective is reached by using maximum parsimony and maximum likelihood to reconstruct the ancestral sequences of the LUCA for two pairs of sets of paralogous protein sequences by means of the phylogenetic tree topology derived from the small subunit ribosomal RNA, even if this is rooted in all three possible ways. The thermophily index of all the reconstructed ancestral sequences of the LUCA belongs to the set of the thermophile/hyperthermophile sequences, thus supporting the hypotheses that see the LUCA as a thermophile or a hyperthermophile.     

The universal ancestor was a thermophile or a hyperthermophile: tests and further evidence

The existence of a correlation between the optimal growth temperature of various organisms and a thermophily index (based on the propensity of amino acids to enter more frequently into the proteins of thermophiles/hyperthermophiles) allows inferences to be made on the mesophilic or thermophilic nature of the last universal common ancestor (LUCA). By reconstructing the ancestral sequences of the various ancestors using methods based on maximum likelihood and maximum parsimony, these sequences can be attributed to the mesophiles or (hyper)thermophiles and the following conclusions can be drawn. (1) There is no evidence that the LUCA might have been a mesophile and observations seem to imply that the LUCA was a thermophile or a hyperthermophile; (2) The ancestors of the Archaea and Bacteria domains seem to be (hyper)thermophiles while that of the Eukarya domain turns out to be a mesophile. These conclusions are independent of both (i) where the root is located on the topology of the universal tree (based on that of the small subunit ribosomal RNA) and (ii) the presence of hyperthermophile bacteria near the node of the Bacteria domain ancestor. These conclusions are easier to interpret in the light of the hypotheses that see the origin of life taking place at a high temperature.     

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 Universal Ancestor and the Ancestor of Bacteria Were Hyperthermophiles

The definition of the node of the last universal common ancestor (LUCA) is justified in a topology of the unrooted universal tree. This definition allows previous analyses based on paralogous proteins to be extended to orthologous ones. In particular, the use of a thermophily index (based on the amino acids propensity to enter the [hyper] thermophile proteins more frequently) and its correlation with the optimal growth temperature of the various organisms allow inferences to be made on the habitat in which the LUCA lived. The reconstruction of ancestral sequences by means of the maximum likelihood method and their attribution to the set of mesophilic or hyperthermophilic sequences have led to the following conclusions: the LUCA was a hyperthermophile organism, as were the ancestors of the Archaea and Bacteria domains, while the ancestor of the Eukarya domain was a mesophile. These conclusions are independent of the presence of hyperthermophile bacteria in the sample of sequences used in the analysis and are therefore independent of whether or not these are the first lines of divergence in the Bacteria domain, as observed in the topology of the universal tree of ribosomal RNA. These conclusions are thus more easily understood under the hypothesis that the origin of life took place at a high temperature.     

The evolution of the genetic code took place in an anaerobic environment

We have compared orthologous proteins from an aerobic organism, Cytophaga hutchinsonii, and from an obligate anaerobe, Bacteroides thetaiotaomicron. This comparison allows us to define the oxyphobic ranks of amino acids, i.e. a scale of the relative sensitivity to oxygen of the amino acid residues. The oxyphobic index (OI), which can be simply obtained from the amino acids' oxyphobic ranks, can be associated to any protein and therefore 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. Sampling of the OI variable from the proteins of obligate anaerobes and aerobes has established that the OI value of the genetic code is not significantly different from the mean OI value of anaerobe proteins, while it is different from that of aerobe proteins. This observation would seem to suggest that the terminal phases of the evolution of genetic code organization took place in an anaerobic environment. This result is discussed in the framework of hypotheses suggested to explain the timing of the evolutionary appearance of the aerobic metabolism.     

Analysis of the heat-shock response displayed by two Chaetomium species originating from different thermal environments.

Three features of the heat shock response, reorganization of protein expression, intracellular accumulation of trehalose, and alteration in unsaturation degree of fatty acids were investigated in the thermophilic fungus Chaetomium thermophile and compared to the response displayed by a closely related mesophilic species, C. brasiliense. Thermophilic heat shock response paralleled the mesophilic response in many respects like (i) the temperature difference observed between normothermia and the upper limit of translational activity, (ii) the transient nature of the heat shock response at the level of protein expression including both the induction of heat shock proteins (HSPs) as well as the repression of housekeeping proteins, (iii) the presence of representatives of high-molecular-weight HSPs families, (iv) intracellular accumulation of trehalose, and finally (v) modifications in fatty acid composition. On the other hand, a great variability between the two organisms was observed for the proteins expressed during stress, in particular a protein of the HSP60 family that was only observed in C. thermophile. This peptide was also present constitutively at normal temperature and may thus fulfil thermophilic functions. It is shown that accumulation of trehalose does not play a part in thermophily but is only a stress response. C. thermophile contains less polyunsaturated fatty acids at normal temperature than C. brasiliense, a fact that can be directly related to thermophily. When subjected to heat stress, both organisms tended to accumulate shorter and less unsaturated fatty acids.     

Substitution of aspartic acid with glutamic acid increases the unfolding transition temperature of a protein

Proteins from thermophiles are more stable than those from mesophiles. Several factors have been suggested as causes for this greater stability, but no general rule has been found. The amino acid composition of thermophile proteins indicates that the content of polar amino acids such as Asn, Gln, Ser, and Thr is lower, and that of charged amino acids such as Arg, Glu, and Lys is higher than in mesophile proteins. Among charged amino acids, however, the content of Asp is even lower in thermophile proteins than in mesophile proteins. To investigate the reasons for the lower occurrence of Asp compared to Glu in thermophile proteins, Glu was substituted with Asp in a hyperthermophile protein, MjTRX, and Asp was substituted with Glu in a mesophile protein, ETRX. Each substitution of Glu with Asp decreased the Tm of MjTRX by about 2 degrees C, while each substitution of Asp with Glu increased the Tm of ETRX by about 1.5 degrees C. The change of Tm destabilizes the MjTRX by 0.55 kcal/mol and stabilizes the ETRX by 0.45 kcal/mol in free energy.     

The Level and Landscape of Optimization in the Origin of the Genetic Code

We consider a model of the origin of genetic code organization incorporating the biosynthetic relationships between amino acids and their physicochemical properties. We study the behavior of the genetic code in the set of codes subject both to biosynthetic constraints and to the constraint that the biosynthetic classes of amino acids must occupy only their own codon domain, as observed in the genetic code. Therefore, this set contains the smallest number of elements ever analyzed in similar studies. Under these conditions and if, as predicted by physicochemical postulates, the amino acid properties played a fundamental role in genetic code organization, it can be expected that the code must display an extremely high level of optimization. This prediction is not supported by our analysis, which indicates, for instance, a minimization percentage of only 80%. These observations can therefore be more easily explained by the coevolution theory of genetic code origin, which postulates a role that is important but not fundamental for the amino acid properties in the structuring of the code. We have also investigated the shape of the optimization landscape that might have arisen during genetic code origin. Here, too, the results seem to favor the coevolution theory because, for instance, the fact that only a few amino acid exchanges would have been sufficient to transform the genetic code (which is not a local minimum) into a much better optimized code, and that such exchanges did not actually take place, seems to suggest that, for instance, the reduction of translation errors was not the main adaptive theme structuring the genetic code.     

Compositional and Structural Features Related to Thermal Stability in the Archaea SRP19 and SRP54 Signal Recognition Particle Proteins

To survive at high temperature, thermophile organisms must adapt their biomolecules. In both nucleic acids and proteins, this adaptation involves a vast array of compositional and structural modifications. The archaea stand out as the only group of organisms that have species capable of growing at temperatures ranging from 0 to 110°C. In this study, we have used the archaea genome datasets to identify molecular trends related to thermal adaptation in the protein components (SRP19 and SRP54) of the signal recognition particle (SRP). Using comparative genomics and secondary structure homology modeling we have detected significant differences in the amino acids composition and distribution between the SRP proteins of thermophile and mesophile archaea. These include: a significant increase in the thermophile SRP proteins of the frequency of charged amino acids able to participate in electrostatic interactions which contribute to stabilize proteins; decreased content of both thermolabile and small/tiny amino acids which usually contribute to protein flexibility; and a significant increase in aliphatic and aromatic amino acids providing good covering and masking to produce hydrophobic pockets involved in stabilizing protein structure. Moreover, a detailed analysis of the four structural and functional domains of the SRP54 indicates a particularly robust correlation between the compositional properties of the M domain and the optimal growth temperature (OGT) of the archaea. The analysis of the bacterial SRP54(Ffh) shows similar adaptations to the OGT. Thus, natural selection has adapted the SRP proteins to the OGT of the archaea and bacteria species by modifying both, their amino acids composition and distribution.     

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’.     

Analysis of the Heat-Shock Response Displayed by Two Chaetomium Species Originating from Different Thermal Environments

Three features of the heat shock response, reorganization of protein expression, intracellular accumulation of trehalose, and alteration in unsaturation degree of fatty acids were investigated in the thermophilic fungus Chaetomium thermophile and compared to the response displayed by a closely related mesophilic species, C. brasiliense. Thermophilic heat shock response paralleled the mesophilic response in many respects like (i) the temperature difference observed between normothermia and the upper limit of translational activity, (ii) the transient nature of the heat shock response at the level of protein expression including both the induction of heat shock proteins (HSPs) as well as the repression of housekeeping proteins, (iii) the presence of representatives of high-molecular-weight HSPs families, (iv) intracellular accumulation of trehalose, and finally (v) modifications in fatty acid composition. On the other hand, a great variability between the two organisms was observed for the proteins expressed during stress, in particular a protein of the HSP60 family that was only observed in C. thermophile. This peptide was also present constitutively at normal temperature and may thus fulfil thermophilic functions. It is shown that accumulation of trehalose does not play a part in thermophily but is only a stress response. C. thermophile contains less polyunsaturated fatty acids at normal temperature than C. brasiliense, a fact that can be directly related to thermophily. When subjected to heat stress, both organisms tended to accumulate shorter and less unsaturated fatty acids.     

The evolution of aminoacyl-tRNA synthetases,the biosynthetic pathways of amino acids and the genetic code

In this paper the partition metric is used to compare binary trees deriving from (i) the study of the evolutionary relationships between aminoacyl-tRNA synthetases, (ii) the physicochemical properties of amino acids and (iii) the biosynthetic relationships between amino acids. If the tree defining the evolutionary relationships between aminoacyl-tRNA synthetases is assumed to be a manifestation of the mechanism that originated the organization of the genetic code, then the results appear to indicate the following: the hypothesis that regards the genetic code as a map of the biosynthetic relationships between amino acids seems to explain the organization of the genetic code, at least as plausibly as the hypotheses that consider the physicochemical properties of amino acids as the main adaptive theme that lead to the structuring of the code.     

The Historical Factor: The Biosynthetic Relationships Between Amino Acids and Their Physicochemical Properties in the Origin of the Genetic Code

Two forces are in general, hypothesized to have influenced the origin of the organization of the genetic code: the physicochemical properties of amino acids and their biosynthetic relationships. In view of this, we have considered a model incorporating these two forces. In particular, we have studied the optimization level of the physicochemical properties of amino acids in the set of amino acid permutation codes that respects the biosynthetic relationships between amino acids. Where the properties of amino acids are represented by polarity and molecular volume we obtain indetermination percentages in the organization of the genetic code of approximately 40%. This indicates that the contingent factor played a significant role in structuring the genetic code. Furthermore, this result is in agreement with the genetic code coevolution hypothesis, which attributes a merely ancillary role to the properties of amino acids while it suggests that it was their biosynthetic relationships that organized the code. Furthermore, this result does not favor the stereochemical models proposed to explain the origin of the genetic code. On the other hand, where the properties of amino acids are represented by polarity alone, we obtain an indetermination percentage of at least 21.5%. This might suggest that the polarity distances played an important role and would therefore provide evidence in favor of the physicochemical hypothesis of genetic code origin. Although, overall, the analysis might have given stronger support to the latter hypothesis, this did not actually occur. The results are therefore discussed in the context of the different theories proposed to explain the origin of the genetic code. Received: 10 September 1996 / Accepted: 3 March 1997     

Aminoacyl-tRNA synthetases database Y2K

The aminoacyl-tRNA synthetases (AARS) are a diverse group of enzymes that ensure the fidelity of transfer of genetic information from DNA into protein. They catalyse the attachment of amino acids to transfer RNAs and thereby establish the rules of the genetic code by virtue of matching the nucleotide triplet of the anticodon with its cognate amino acid. Currently, 818 AARS primary structures have been reported from archaebacteria, eubacteria, mitochondria, chloro-plasts and eukaryotic cells. The database is a compilation of the amino acid sequences of all AARSs, known to date, which are available as separate entries or alignments of related proteins via the WWW at http://rose.man.poznan.pl/aars/index.html     

编码序列的（G＋C）%与蛋白质的耐热性相关性分析

运用计算机统计方法,对以木糖异构酶为主的几个蛋白质家族的核酸和氨基酸序列进行分析,发现密码子各位上的（Ｇ＋Ｃ）％与编码序列的（Ｇ＋Ｃ）％成线性正相关,大多数氨基酸的含量与编码序列的（Ｇ＋Ｃ）％也存在相关性,按其相关性,将氨基酸分为正相关,负相关和不相关３类,对木糖异构酶氨基酸序列和酶的耐热性的统计发现,那些在统计学上显著的,可能提高蛋白质耐热性的氨基酸替换,往往伴随关编码序列中ＧＣ含量的上升,这提     

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.     

The Early Phases of Genetic Code Origin: Conjectures on the Evolution of Coded Catalysis

A review of the most significant contributions on the early phases of genetic code origin is presented. After stressing the importance of the key intermediary role played in protein synthesis, by peptidyl-tRNA, which is attributed with a primary function in ancestral catalysis, the general lines leading to the codification of the first amino acids in the genetic code are discussed. This is achieved by means of a model of protoribosome evolution which sees protoribosome as the central organiser of ancestral biosynthesis and the mediator of the encounter between compounds (metabolite-pre-tRNAs) and catalysts (peptidyl-pre-tRNAs). The encounter between peptidyl-pre-tRNA catalysts in protoribosome is favoured by metabolic pre-mRNAs and later resulted (given the high temperature at which this evolution is supposed to have taken place) in the evolution of mRNAs with codons of the type GNS. These mRNAs codified only for those amino acids that the coevolution theory of genetic code origin sees as the precursors of all other amino acids. Some aspects of the model here discussed might be rendered real by the transfer-messenger RNA molecule (tmRNA) which is here considered a molecular fossil of ancestral protein synthesis.     

A Novel Theory on the Origin of the Genetic Code: A GNC-SNS Hypothesis

We have previously proposed an SNS hypothesis on the origin of the genetic code (Ikehara and Yoshida 1998). The hypothesis predicts that the universal genetic code originated from the SNS code composed of 16 codons and 10 amino acids (S and N mean G or C and either of four bases, respectively). But, it must have been very difficult to create the SNS code at one stroke in the beginning. Therefore, we searched for a simpler code than the SNS code, which could still encode water-soluble globular proteins with appropriate three-dimensional structures at a high probability using four conditions for globular protein formation (hydropathy, α-helix, β-sheet, and β-turn formations). Four amino acids (Gly [G], Ala [A], Asp [D], and Val [V]) encoded by the GNC code satisfied the four structural conditions well, but other codes in rows and columns in the universal genetic code table do not, except for the GNG code, a slightly modified form of the GNC code. Three three-amino acid systems ([D], Leu and Tyr; [D], Tyr and Met; Glu, Pro and Ile) also satisfied the above four conditions. But, some amino acids in the three systems are far more complex than those encoded by the GNC code. In addition, the amino acids in the three-amino acid systems are scattered in the universal genetic code table. Thus, we concluded that the universal genetic code originated not from a three-amino acid system but from a four-amino acid system, the GNC code encoding [GADV]-proteins, as the most primitive genetic code. Received: 11 June 2001 / Accepted: 11 October 2001     

Genetic Code Origin: Are the Pathways of Type Glu-tRNAGln \rightarrow Gln-tRNAGln Molecular Fossils or Not?

A logical-evolutionary analysis is conducted to clarify whether or not pathways of type Glu-tRNAGln \rightarrow Gln-tRNAGln are molecular fossils of the mechanism that gave rise to the evolutionary organization of the genetic code. The result of this analysis is that these pathways are most likely a manifestation of this mechanism. This provides strong evidence in favor of the coevolution theory of genetic code origin, as this theory is based on the amino acid biosynthetic transformation taking place on tRNA-like molecules which imprinted the genetic code structuring. Comments on the different interpretations of these pathways found in the literature are also provided.     

The phylogeny of trnas seems to confirm the predictions of the coevolution theory of the origin of the genetic code

An extensive analysis of the evolutionary relationships existing between transfer RNAs, performed using parsimony algorithms, is presented. After building up an estimate of the tRNA ancestral sequences, these sequences are then compared using certain methods. The results seem to suggest that the coevolution hypothesis (Wong, J.T., 1975, Proc. Natl. Acad. Sci. USA 72, 1909–1912) that sees the genetic code as a map of the biosynthetic relationships between amino acids is further supported by these results, as compared to the hypotheses that see the physicochemical properties of amino acids as the main adaptative theme that led to the structuring of the genetic code.