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

The Mutational Robustness of the Genetic Code and Codon Usage in Environmental Context: A Non-Extremophilic Preference?

The genetic code was evolved, to some extent, to minimize the effects of mutations. The effects of mutations depend on the amino acid repertoire, the structure of the genetic code and frequencies of amino acids in proteomes. The amino acid compositions of proteins and corresponding codon usages are still under selection, which allows us to ask what kind of environment the standard genetic code is adapted to. Using simple computational models and comprehensive datasets comprising genomic and environmental data from all three domains of Life, we estimate the expected severity of non-synonymous genomic mutations in proteins, measured by the change in amino acid physicochemical properties. We show that the fidelity in these physicochemical properties is expected to deteriorate with extremophilic codon usages, especially in thermophiles. These findings suggest that the genetic code performs better under non-extremophilic conditions, which not only explains the low substitution rates encountered in halophiles and thermophiles but the revealed relationship between the genetic code and habitat allows us to ponder on earlier phases in the history of Life.

Phylogenetic analysis of mutational robustness based on codon usage supports that the standard genetic code does not prefer extreme environments

The mutational robustness of the genetic code is rarely discussed in the context of biological diversity, such as codon usage and related factors, often considered as independent of the actual organism's proteome. Here we put the living beings back to picture and use distortion as a metric of mutational robustness. Distortion estimates the expected severities of non-synonymous mutations measuring it by amino acid physicochemical properties and weighting for codon usage. Using the biological variance of codon frequencies, we interpret the mutational robustness of the standard genetic code with regards to their corresponding environments and genomic compositions (GC-content). Employing phylogenetic analyses, we show that coding fidelity in physicochemical properties can deteriorate with codon usages adapted to extreme environments and these putative effects are not the artefacts of phylogenetic bias. High temperature environments select for codon usages with decreased mutational robustness of hydrophobic, volumetric, and isoelectric properties. Selection at high saline concentrations also leads to reduced fidelity in polar and isoelectric patterns. These show that the genetic code performs best with mesophilic codon usages, strengthening the view that LUCA or its ancestors preferred lower temperature environments. Taxonomic implications, such as rooting the tree of life, are also discussed.

Natural pyrrolysine-biased translation of stop codons in mitochondrial peptides entirely coded by expanded codons

During the noncanonical deletion transcription, k nucleotides are systematically skipped/deleted after each transcribed trinucleotide producing deletion-RNAs (delRNAs). Peptides matching delRNAs either result from (a) canonical translation of delRNAs; or (b) noncanonical translation of regular transcripts along expanded codons. Only along frame "0" (start site) (a) and (b) produce identical peptides. Here, mitochondrial mass spectrometry data analyses assume expanded codon/del-transcription with 3 + k (k from 0 to 12) nucleotides. Detected peptides map preferentially on previously identified delRNAs. More peptides were detected for k (1-12) when del-transcriptional and expanded codon translations start sites coincide (i.e. the 0^th frame) than for frames +1 or +2. Hence, both (a) and (b) produced peptides identified here. Biases for frame 0 decrease for k > 2, reflecting codon/anticodon expansion limits. Further analyses find preferential pyrrolysine insertion at stop codons, suggesting Pyl-specific mitochondrial suppressor tRNAs loaded by Pyl-specific tRNA synthetases with unknown origins. Pyl biases at stops are stronger for regular than expanded codons suggesting that Pyl-tRNAs are less competitive with near-cognate tRNAs in expanded codon contexts. Statistical biases for these findings exclude that detected peptides are experimental and/or bioinformatic artefacts implying both del-transcription and expanded codons translation occur in human mitochondria.

Chimeric Translation for Mitochondrial Peptides: Regular and Expanded Codons

Frameshifting protein translation occasionally results from insertion of amino acids at isolated mono- or dinucleotide-expanded codons by tRNAs with expanded anticodons. Previous analyses of two different types of human mitochondrial MS proteomic data (Fisher and Waters technologies) detect peptides entirely corresponding to expanded codon translation. Here, these proteomic data are reanalyzed searching for peptides consisting of at least eight consecutive amino acids translated according to regular tricodons, and at least eight adjacent consecutive amino acids translated according to expanded codons. Both datasets include chimerically translated peptides (mono- and dinucleotide expansions, 42 and 37, respectively). The regular tricodon-encoded part of some chimeric peptides corresponds to standard human mitochondrial proteins (mono- and dinucleotide expansions, six (AT6, CytB, ND1, 2xND2, ND5) and one (ND1), respectively). Chimeric translation probably increases the diversity of mitogenome-encoded proteins, putatively producing functional proteins. These might result from translation by tRNAs with expanded anticodons, or from regular tricodon translation of RNAs where transcription/posttranscriptional edition systematically deleted mono- or dinucleotides after each trinucleotide. The pairwise matched combination of adjacent peptide parts translated from regular and expanded codons strengthens the hypothesis that translation of stretches of consecutive expanded codons occurs. Results indicate statistical translation producing distributions of alternative proteins. Genetic engineering should account for potential unexpected, unwanted secondary products.

Primordial chemistry and enzyme evolution in a hot environment

Ever since the publication of Darwin's Origin of Species, questions have been raised about whether enough time has elapsed for living organisms to have reached their present level of complexity by mutation and natural selection. More recently, it has become apparent that life originated very early in Earth's history, and there has been controversy as to whether life originated in a hot or cold environment. This review describes evidence that rising temperature accelerates slow reactions disproportionately, and to a much greater extent than has been generally recognized. Thus, the time that would have been required for primordial chemistry to become established would have been abbreviated profoundly at high temperatures. Moreover, if the catalytic effect of a primitive enzyme (like that of modern enzymes) were to reduce a reaction's heat of activation, then the rate enhancement that it produced would have increased as the surroundings cooled, quite aside from changes arising from mutation (which is itself highly sensitive to temperature). Some nonenzymatic catalysts of slow reactions, including PLP as a catalyst of amino acid decarboxylation, and the Ce(IV) ion as a catalyst of phosphate ester hydrolysis, have been shown to meet that criterion. The work reviewed here suggests that elevated temperatures collapsed the time required for early evolution on Earth, furnishing an appropriate setting for exploring the vast range of chemical possibilities and for the rapid evolution of enzymes from primitive catalysts.

Comparative analysis of barophily-related amino acid content in protein domains of Pyrococcus abyssi and Pyrococcus furiosus

Amino acid substitution patterns between the nonbarophilic Pyrococcus furiosus and its barophilic relative P. abyssi confirm that hydrostatic pressure asymmetry indices reflect the extent to which amino acids are preferred by barophilic archaeal organisms. Substitution patterns in entire protein sequences, shared protein domains defined at fold superfamily level, domains in homologous sequence pairs, and domains of very ancient and very recent origin now provide further clues about the environment that led to the genetic code and diversified life. The pyrococcal proteomes are very similar and share a very early ancestor. Relative amino acid abundance analyses showed that biases in the use of amino acids are due to their shared fold superfamilies. Within these repertoires, only two of the five amino acids that are preferentially barophilic, aspartic acid and arginine, displayed this preference significantly and consistently across structure and in domains appearing in the ancestor. The more primordial asparagine, lysine and threonine displayed a consistent preference for nonbarophily across structure and in the ancestor. Since barophilic preferences are already evident in ancient domains that are at least ~3 billion year old, we conclude that barophily is a very ancient trait that unfolded concurrently with genetic idiosyncrasies in convergence towards a universal code.

The origin of the genetic code in the ocean abysses: new comparisons confirm old observations

I have analysed the amino acid substitution pattern between two pairs of nonbarophilic-barophilic organisms in order to confirm previous results. Indeed, the pattern deriving from a different pair of such organisms led to establish that the origin of the genetic code might have occurred in the ocean abysses. The hydrostatic pressure asymmetry indices computable from these matrices of amino acid substitutions confirm the correlation previously observed, even when differences in GC content are accounted for. As the three substitution matrices are independent between them, the old conclusion that the genetic code originated in the ocean abysses is considerably strengthened.

Behavior of a hammerhead ribozyme in aqueous solution at medium to high temperatures

The "RNA world" hypothesis proposes that--early in the evolution of life--RNA molecules played important roles both in information storage and in enzymatic functions. However, this hypothesis seems to be inconsistent with the concept that life may have emerged under hydrothermal conditions since RNA molecules are considered to be labile under such extreme conditions. Presently, the possibility that the last common ancestor of the present organisms was a hyperthermophilic organism which is important to support the hypothesis of the hydrothermal origin of life has been subject of strong discussions. Consequently, it is of importance to study the behavior of RNA molecules under hydrothermal conditions from the viewpoints of stability, catalytic functions, and storage of genetic information of RNA molecules and determination of the upper limit of temperature where life could have emerged. In the present work, self-cleavage of a natural hammerhead ribozyme was examined at temperatures 10-200 °C. Self-cleavage was investigated in the presence of Mg(2+), which facilitates and accelerates this reaction. Self-cleavage of the hammerhead ribozyme was clearly observed at temperatures up to 60 °C, but at higher temperatures self-cleavage occurs together with hydrolysis and with increasing temperature hydrolysis becomes dominant. The influence of the amount of Mg(2+) on the reaction rate was also investigated. In addition, we discovered that the reaction proceeds in the presence of high concentrations of monovalent cations (Na(+) or K(+)), although very slowly. Furthermore, at high temperatures (above 60 °C), monovalent cations protect the ribozyme against degradation.

Optimization models and the structure of the genetic code

The codon assignment of the quasi-universal genetic code can be assumed to have resulted from the evolutionary pressures that prevailed when the code was still evolving. Here, we review studies of the structure of the genetic code based on optimization models. We also review studies that, from the structure of the code, attempt to derive aspects of the primordial circumstances in which the genetic code froze. Different rationales are summarized, compared with experimental data, discussed in the context of the transition from a RNA world to a DNA-protein world, and linked to the emergence of the last universal common ancestor.

The last universal common ancestor: emergence, constitution and genetic legacy of an elusive forerunner

BACKGROUND

Since the reclassification of all life forms in three Domains (Archaea, Bacteria, Eukarya), the identity of their alleged forerunner (Last Universal Common Ancestor or LUCA) has been the subject of extensive controversies: progenote or already complex organism, prokaryote or protoeukaryote, thermophile or mesophile, product of a protracted progression from simple replicators to complex cells or born in the cradle of "catalytically closed" entities? We present a critical survey of the topic and suggest a scenario.

RESULTS

LUCA does not appear to have been a simple, primitive, hyperthermophilic prokaryote but rather a complex community of protoeukaryotes with a RNA genome, adapted to a broad range of moderate temperatures, genetically redundant, morphologically and metabolically diverse. LUCA's genetic redundancy predicts loss of paralogous gene copies in divergent lineages to be a significant source of phylogenetic anomalies, i.e. instances where a protein tree departs from the SSU-rRNA genealogy; consequently, horizontal gene transfer may not have the rampant character assumed by many. Examining membrane lipids suggest LUCA had sn1,2 ester fatty acid lipids from which Archaea emerged from the outset as thermophilic by "thermoreduction," with a new type of membrane, composed of sn2,3 ether isoprenoid lipids; this occurred without major enzymatic reconversion. Bacteria emerged by reductive evolution from LUCA and some lineages further acquired extreme thermophily by convergent evolution. This scenario is compatible with the hypothesis that the RNA to DNA transition resulted from different viral invasions as proposed by Forterre. Beyond the controversy opposing "replication first" to metabolism first", the predictive arguments of theories on "catalytic closure" or "compositional heredity" heavily weigh in favour of LUCA's ancestors having emerged as complex, self-replicating entities from which a genetic code arose under natural selection.

CONCLUSION

Life was born complex and the LUCA displayed that heritage. It had the "body "of a mesophilic eukaryote well before maturing by endosymbiosis into an organism adapted to an atmosphere rich in oxygen. Abundant indications suggest reductive evolution of this complex and heterogeneous entity towards the "prokaryotic" Domains Archaea and Bacteria. The word "prokaryote" should be abandoned because epistemologically unsound.

REVIEWERS

This article was reviewed by Anthony Poole, Patrick Forterre, and Nicolas Galtier.

How does gene expression level contribute to thermophilic adaptation of prokaryotes? An exploration based on predictors

By analyzing the predicted gene expression levels of 33 prokaryotes with living temperature span from <10 degrees C to >100 degrees C, a universal positive correlation was found between the percentage of predicted highly expressed genes and the organisms' optimal growth temperature. A physical interpretation of the correlation revealed that highly expressed genes are statistically more thermostable than lowly expressed genes. These findings show the possibility of the significant contribution of gene expression level to the prokaryotic thermal adaptation and provide evidence for the translational selection pressure on the thermostability of natural proteins during evolution.

The tree of life might be rooted in the branch leading to Nanoarchaeota

It is suggested that the tree of life might be rooted in the domain of the Archaea, in the branch leading to the phylum of Nanoarchaeota. This hypothesis seems to be corroborated by the uniqueness and ancestrality of some traits possessed by Nanoarchaeum equitans, such as split genes separately codifying for the 5' and 3' halves of the tRNA molecule. These half genes are the oldest ancestral form of gene we have ever seen. This, along with the absence of operons from the genome of N. equitans, would seem to indicate that this genome is a molecular fossil of the evolutionary stage which the ancestral genomes had reached when the first lines of divergence were established. Moreover, the late appearance of DNA coinciding with the rooting of the universal phylogenetic tree would make the genome of N. equitans a witness to this fundamental event.

Extremely thermophilic translation system in the common ancestor commonote: ancestral mutants of Glycyl-tRNA synthetase from the extreme thermophile Thermus thermophilus

Based on phylogenetic analysis of 16 S and 18 S rRNAs, the common ancestor of all organisms (Commonote) was proposed to be hyperthermophilic. We have previously tested this hypothesis using enzymes with ancestral residues that are inferred by molecular phylogenetic analysis. The ancestral mutant enzymes involved in metabolic systems show higher thermal stability than wild-type enzymes, consistent with the hyperthermophile common ancestor hypothesis. Here, we have extended the experiments to include an enzyme of the translation system, glycyl-tRNA synthetase (GlyRS). The translation system often shows a phylogenetic tree that is similar to the rRNA tree. Thus, it is likely that the tree represents the evolutionary route of the organisms. The maximum-likelihood tree of alpha(2) type GlyRS was constructed. From this analysis the ancestral sequence of GlyRS was deduced and individual or pairs of ancestral residues were introduced into Thermus thermophilus GlyRS. The ancestral mutants were expressed in Escherichia coli, purified and activity measured. The thermostability of eight mutated proteins was evaluated by CD (circular dichroism) measurements. Six mutants showed higher thermostability than wild-type enzyme and seven mutants showed higher activity than wild-type enzyme at 70 degrees C, suggesting an extremely thermophilic translation system in the common ancestor Commonote.

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.

In vitro selection of halo-thermophilic RNA reveals two families of resistant RNA

The "RNA world" hypothesis proposes that early in the evolution of life, RNA was responsible both for the storage and transfer of genetic information and for the catalysis of biochemical reactions. One of the problems of the hypothesis is that RNA is known to be temperature sensitive. Nevertheless, different types of sequences with a thermostable phenotype may exist. In order to test this possibility, we applied an in vitro evolution method (SELEX) to isolate RNA molecules that are resistant at high temperatures (80 degrees C for 65 h) and high salt concentrations (2 M NaCl). The sequences of the resulting cloned halo-thermophilic RNAs can be grouped in two families (I and II) possessing very different thermal and chemical stabilities and very different secondary structures. The selected RNA molecules illustrate two different possibilities leading to thermal resistance which may be related to primitive conditions. We propose that members of family I constitute a good means of storing sequence information while members of family II are less efficient but replicate faster in early steps of the SELEX. These selected RNA behaviors may be related to primitive conditions and could allow to define limits for survival, and demonstrate that what is at stake for RNA molecules, as for living organisms, is survival and reproduction.

Thermostability of ancestral mutants of Caldococcus noboribetus isocitrate dehydrogenase

We constructed mutant genes of Caldococcus noboribetus isocitrate dehydrogenase containing ancestral amino acid residues that were inferred using the maximal likelihood method and a composite phylogenetic tree of isocitrate dehydrogenase and 3-isopropylmalate dehydrogenase. The mutant genes were expressed in Escherichia coli and the protein products purified. Thermostabilities, reported as the half-inactivation temperatures, for the purified enzymes were determined and compared with that of the wild-type enzyme. Four of the five mutant enzymes have greater thermal stabilities than wild-type isocitrate dehydrogenase. The results are compatible with the hyperthermophilic universal ancestor (commonote) hypothesis. Incorporation of ancestral residues into a modern-day protein sequence can be used to improve protein thermostability.

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.

A comparison of proteins from Pyrococcus furiosus and Pyrococcus abyssi: barophily in the physicochemical properties of amino acids and in the genetic code

A comparison is made between orthologous proteins from a nonbarophile (Pyrococcus furiosus) and a barophile (Pyrococcus abyssi) organism. The pattern of asymmetries in the amino acid substitution process identifies the amino acids arginine, serine, glycine, valine and aspartic acid as those having the most barophilic behaviour, and tyrosine and glutamine as the least barophilic. The construction of a hydrostatic pressure asymmetry index (PAI) which orders the amino acids from the most barophilic to the least barophilic makes it possible to visualise the amino acid properties that best explain barophily. The polarity of amino acids is positively correlated to the PAI values, i.e., on average, the more polar amino acids possess a higher PAI value, that is to say they are more barophilic. Moreover, the "size" of amino acids (molecular weight) is negatively correlated to the PAI value, that is to say that, on average, high PAI values are associated to "small" amino acids which are therefore more barophilic than "larger" ones. These two amino acid properties are the same ones that are known for having been important in affecting the origin of genetic code organisation. All the above, as well as the significant and positive correlation between the number of codons attributed to the amino acids in the genetic code and the PAI values, seem to favour the hypothesis that genetic code structuring took place under high hydrostatic pressure.

The ocean abysses witnessed the origin of the genetic code

The comparison of proteins from a non-barophilous and a barophilous organism makes it possible to define the barophily ranks of amino acids. The correlation of these ranks with the number of codons attributed to amino acids in the genetic code, together with another straightforward argument based on an optimisation percentage of a barophily index (BI) (easily defined by barophily ranks) which can be associated to the genetic code table, suggest that the genetic code originated under high hydrostatic pressure. Moreover, as the BI value can be calculated for the sequence of any protein, it also makes it possible to define the BI for the genetic code if the number of 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, sampling the BI variable between many non-barophile organisms and from many proteins of a single non-barophile organism leads to the conclusion that the BI value of the genetic code is not typical of these organisms. Whereas, since the genetic code BI value is statistically higher than that of these non-barophile organisms, it supports the hypothesis that genetic code structuring took place under high hydrostatic pressure.

The origin of the genetic code: theories and their relationships, a review

A review of the main theories proposed to explain the origin of the genetic code is presented. I analyze arguments and data in favour of different theories proposed to explain the origin of the organization of the genetic code. It is possible to suggest a mechanism that makes compatible the different theories of the origin of the code, even if these are based on a historical or physicochemical determinism and thus appear incompatible by definition. Finally, I discuss the question of why a given number of synonymous codons was attributed to the amino acids in the genetic code.

The ancestor of the Bacteria domain was a hyperthermophile

Brochier and Philippe have recently re-analysed the phylogeny of ribosomal RNA using only multiple alignment positions with no phylogenetic noise. They conclude that the first branch of divergence in the Bacteria domain comprises Planctomycetales and not hyperthermophile bacteria as in classic phylogeny. In the present paper I examine the robustness of their conclusions. (1) A site-by-site reading of the RNA alignments of Brochier and Philippe seems to suggest that the number of nucleotide positions used in their analysis is not sufficiently high and their phylogenetic analysis is consequently not robust. Furthermore, (2) a different method for selecting positions with no phylogenetic noise from the rRNA alignment relocates the Aquificales and the Thermotogales as the first lines of divergence in the Bacteria domain, and sets Planctomycetales as the third branch of divergence in the phylogenetic tree built from these selected positions. These findings consolidate the hypothesis that the ancestor of the Bacteria domain was a hyperthermophile and, more generally, that the last universal common ancestor might also have been one.

Genetic material in the early evolution of bacteria

DNA and RNA are nucleic acids that cells and viruses use to produce copies of themselves. However, there is an immense paucity of knowledge on how these nucleic acids originated and changed as early bacteria became capable of growth and cell division. One possibility is that parallel evolution of the genetic code and protein synthesis was required for assembly of the first cells capable of growth and division. It is also possible that DNA-RNA duplices were intermediate genetic material in the early assembly of the first cells. These ideas will be discussed as well as other aspects of the assembly of the first cells on the Earth.

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.

Early assembly of cellular life

Popular hypotheses that attempt to explain the origin of prebiotic molecules and cellular life capable of growth and division are not always agreed upon. In this manuscript, information on early bacterial life on Earth is examined using information from several disciplines. For example, knowledge can be integrated from physics, thermodynamics, planetary sciences, geology, biogeochemistry, lipid chemistry, primordial cell structures, cell and molecular biology, microbiology, metabolism and genetics. The origin of life also required a combination of elements, compounds and environmental physical-chemical conditions that allowed cells to assemble in less than a billion years. This may have been widespread in the subsurface of the early Earth located at microscopic physical domains.

Thermostability of membrane protein helix-helix interaction elucidated by statistical analysis

A prerequisite for the survival of (micro)organisms at high temperatures is an adaptation of protein stability to extreme environmental conditions. In contrast to soluble proteins, where many factors have already been identified, the mechanisms by which the thermostability of membrane proteins is enhanced are almost unknown. The hydrophobic membrane environment constrains possible stabilizing factors for transmembrane domains, so that a difference might be expected between soluble and membrane proteins. Here we present sequence analysis of predicted transmembrane helices of the genomes from eight thermophilic and 12 mesophilic organisms. A comparison of the amino acid compositions indicates that more polar residues can be found in the transmembrane helices of thermophilic organisms. Particularly, the amino acids aspartic acid and glutamic acid replace the corresponding amides. Cysteine residues are found to be significantly decreased by about 70% in thermophilic membrane domains suggesting a non-specific function of most cysteine residues in transmembrane domains of mesophilic organisms. By a pair-motif analysis of the two sets of transmembrane helices, we found that the small residues glycine and serine contribute more to transmembrane helix-helix interactions in thermophilic organisms. This may result in a tighter packing of the helices allowing more hydrogen bond formation.

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.