The ocean abysses witnessed the origin of the genetic code

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.

Translation of mitochondrial swinger RNAs according to tri-, tetra- and pentacodons

Transcriptomes and proteomes include RNA and protein fragments not matching regular transcription/translation. Some 'non-canonical' mitochondrial transcripts match mitogenomes after assuming one among 23 systematic exchanges between nucleotides, producing swinger RNAs (nine symmetric, X↔Y, example C↔T; 14 asymmetric, X→Y→Z→X, example A→T→G→A) in GenBank's EST database. Here, reanalyzes of (a) public human mitochondrial transcriptome data (Illumina: RNA-seq) allowed to detect mitochondrial swinger RNAs for all 23 exchanges and (b) independent public human mitochondrial trypsinized proteomic mass spectrometry data allowed to detect peptides predicted from translation of parts of swinger-transformed mitogenomes covered by detected swinger reads. RNA-seq and previous EST swinger transcript data converge. Swinger RNA translation frequently inserts various amino acids at stop codons. Swinger RNA-peptide associations exist also for peptides matching systematically frameshifting translation, peptides entirely coded by tetra- and pentacodons (regular codons expanded by silent mononucleotides at 4th, and silent dinucleotides at 4th and 5th position(s), respectively). Swinger peptides differ from regular mitochondrial proteins: not membrane embedded, reflect warmer, anaerobic, low resource conditions, reminding a free-living ancestor. Tetra- and pentacoded peptides associate with low, high GC contents, respectively, suggesting expanded codon translations associate with thermic stresses. Results confirm experimentally predicted swinger, tetra- and pentacoded mitochondrial peptides, increasing mitogenomic coding density.

Insights into the molecular basis of piezophilic adaptation: Extraction of piezophilic signatures

Piezophiles are the organisms which can successfully survive at extreme pressure conditions. However, the molecular basis of piezophilic adaptation is still poorly understood. Analysis of the protein sequence adjustments that had taken place during evolution can help to reveal the sequence adaptation parameters responsible for protein functional and structural adaptation at such high pressure conditions. In this current work we have used SVM classifier for filtering strong instances and generated human interpretable rules from these strong instances by using the PART algorithm. These generated rules were analyzed for getting insights into the molecular signature patterns present in the piezophilic proteins. The experiments were performed on three different temperature ranges piezophilic groups, namely psychrophilic-piezophilic, mesophilic-piezophilic, and thermophilic-piezophilic for the detailed comparative study. The best classification results were obtained as we move up the temperature range from psychrophilic-piezophilic to thermophilic-piezophilic. Based on the physicochemical classification of amino acids and using feature ranking algorithms, hydrophilic and polar amino acid groups have higher discriminative ability for psychrophilic-piezophilic and mesophilic-piezophilic groups along with hydrophobic and nonpolar amino acids for the thermophilic-piezophilic groups. We also observed an overrepresentation of polar, hydrophilic and small amino acid groups in the discriminatory rules of all the three temperature range piezophiles along with aliphatic, nonpolar and hydrophobic groups in the mesophilic-piezophilic and thermophilic-piezophilic groups.

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.

Pressure as an environmental parameter for microbial life--a review

Microbial life has been prevailing in the biosphere for the last 3.8 Ga at least. Throughout most of the Earth's history it has experienced a range of pressures; both dynamic pressure when the young Earth was heavily bombarded, and static pressure in subsurface environments that could have served as a refuge and where microbial life nowadays flourishes. In this review, we discuss the extent of high-pressure habitats in early and modern times and provide a short overview of microbial survival under dynamic pressures. We summarize the current knowledge about the impact of microbial activity on biogeochemical cycles under pressures characteristic of the deep subsurface. We evaluate the possibility that pressure can be a limiting parameter for life at depth. Finally, we discuss the open questions and knowledge gaps that exist in the field of high-pressure geomicrobiology.

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.

The first genomic and proteomic characterization of a deep-sea sulfate reducer: insights into the piezophilic lifestyle of Desulfovibrio piezophilus

Desulfovibrio piezophilus strain C1TLV30(T) is a piezophilic anaerobe that was isolated from wood falls in the Mediterranean deep-sea. D. piezophilus represents a unique model for studying the adaptation of sulfate-reducing bacteria to hydrostatic pressure. Here, we report the 3.6 Mbp genome sequence of this piezophilic bacterium. An analysis of the genome revealed the presence of seven genomic islands as well as gene clusters that are most likely linked to life at a high hydrostatic pressure. Comparative genomics and differential proteomics identified the transport of solutes and amino acids as well as amino acid metabolism as major cellular processes for the adaptation of this bacterium to hydrostatic pressure. In addition, the proteome profiles showed that the abundance of key enzymes that are involved in sulfate reduction was dependent on hydrostatic pressure. A comparative analysis of orthologs from the non-piezophilic marine bacterium D. salexigens and D. piezophilus identified aspartic acid, glutamic acid, lysine, asparagine, serine and tyrosine as the amino acids preferentially replaced by arginine, histidine, alanine and threonine in the piezophilic strain. This work reveals the adaptation strategies developed by a sulfate reducer to a deep-sea lifestyle.

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.

Molecular evolution of the hyperthermophilic archaea of the Pyrococcus genus: analysis of adaptation to different environmental conditions

Background

Prokaryotic microorganisms are able to survive and proliferate in severe environmental conditions. The increasing number of complete sequences of prokaryotic genomes has provided the basis for studying the molecular mechanisms of their adaptation at the genomic level. We apply here a computer-based approach to compare the genomes and proteomes from P. furiosus, P. horikoshii, and P. abyssi to identify features of their molecular evolution related to adaptation strategy to diverse environmental conditions.

Results

Phylogenetic analysis of rRNA genes from 26 Pyrococcus strains suggested that the divergence of P. furiosus, P. horikoshii and P. abyssi might have occurred from ancestral deep-sea organisms. It was demonstrated that the function of genes that have been subject to positive Darwinian selection is closely related to abiotic and biotic conditions to which archaea managed to become adapted. Divergence of the P. furiosus archaea might have been due to loss of some genes involved in cell motility or signal transduction, and/or to evolution under positive selection of the genes for translation machinery. In the course of P. horikoshii divergence, positive selection was found to operate mainly on the transcription machinery; divergence of P. abyssi was related with positive selection for the genes mainly involved in inorganic ion transport. Analysis of radical amino acid replacement rate in evolving P. furiosus, P. horikoshii and P. abyssi showed that the fixation rate was higher for radical substitutions relative to the volume of amino acid side-chain.

Conclusions

The current results give due credit to the important role of hydrostatic pressure as a cause of variability in the P. furiosus, P. horikoshii and P. abyssi genomes evolving in different habitats. Nevertheless, adaptation to pressure does not appear to be the sole factor ensuring adaptation to environment. For example, at the stage of the divergence of P. horikoshii and P. abyssi, an essential evolutionary role may be assigned to changes in the trophic chain, namely, acquisition of a consumer status at a high (P. horikoshii) or low level (P. abyssi).

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

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.

Hydrostatic and osmotic pressure study of the hairpin ribozyme

The recent discovery of numerous catalytically active RNAs in various living species as well as the in vitro selection of a large series of RNA aptamers able to bind specifically various molecules such as metabolites and co-factors, emphasize the adaptability of RNAs through the plasticity of their secondary structure. Furthermore, all these observations give support to the "RNA world" hypothesis as a step in the primitive development of life on Earth. On this background, we used high pressure to study the mechanism of action of a model hairpin ribozyme which exhibits self-cleavage and ligation. The activation volume (DeltaV( not equal)) of the cleavage reaction (34+/-4 ml/mol) indicates that an important compaction of the RNA molecule occurs during the reaction and must be accompanied by a significant movement of water molecules . Indeed, such a release of 78+/-4 water molecules per RNA molecule could be measured by complementary osmotic shock experiments. These results are consistent with the information provided by the structural studies which indicate that two loops of the RNA molecule should come into contact for the reaction to occur . The high pressure study of a modified form of the ribozyme whose activity is strictly dependent on the presence of adenine as a co-factor should bring some information about the structural significance of this important DeltaV( not equal) of activation.

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.