Origin and evolution of the genetic code: the universal enigma

Chemical Evolution of Life on Earth

Background/Objectives: The origin of genes and genetics is the story of the coevolution of translation systems and the genetic code. Remarkably, the history of the origin of life on Earth was inscribed and preserved in the sequences of tRNAs. Methods: Sequence logos demonstrate the patterning of pre-life tRNA sequences. Results: The pre-life type I and type II tRNA sequences are known to the last nucleotide with only a few ambiguities. Type I and type II tRNAs evolved from ligation of three 31 nt minihelices of highly patterned and known sequence followed by closely related 9 nt internal deletion(s) within ligated acceptor stems. The D loop 17 nt core was a truncated UAGCC repeat. The anticodon and T 17 nt stem-loop-stems are homologous sequences with 5 nt stems and 7 nt U-turn loops that were selected in pre-life to resist ribozyme nucleases and to present a 3 nt anticodon with a single wobble position. The 7 nt T loop in tRNA was selected to interact with the D loop at the "elbow". The 5'-acceptor stem was based on a 7 nt truncated GCG repeat. The 3'-acceptor stem was based on a complementary 7 nt CGC repeat. In pre-life, ACCA-Gly was a primitive adapter molecule ligated to many RNAs, including tRNAs, to synthesize polyglycine. Conclusions: Analysis of sequence logos of tRNAs from an ancient Archaeon substantiates how the pre-life to life transition occurred on Earth. Polyglycine is posited to have aggregated complex molecular assemblies, including minihelices, tRNAs, cooperating molecules, and protocells, leading to the first life on Earth.

AAV Capsid Modification and Its Influence on Viral Protein Stoichiometry and Packaging Fitness: Current Understandings and Future Direction

The field of gene therapy has witnessed significant advancements in the utilization of Adeno-associated virus (AAV) owing to its inherent biological advantages. Targeted AAV vectors are generated through genetic or chemical modification of the capsid for user-directed purposes. However, this process can result in imbalances in viral protein sequence homogeneity, stoichiometry, and functional transduction vector units, thereby introducing new challenges. This mini review focuses on the ongoing efforts to develop targeted vectors, which inadvertently present unsolicited obstacles for clinical application and provided perspectives on future directions.

Comparative genomics of thermosynechococcaceae and thermostichaceae: insights into codon usage bias

Members of the families Thermosynechococcaceae and Thermostichaceae are well-known unicellular thermophilic cyanobacteria and a non-thermophilic genus Pseudocalidococcus was newly classified into the former. Analysis of the codon usage bias (CUB) of cyanobacterial species inhabiting different thermal and non-thermal niches will benefit the understanding of their genetic and evolutionary characteristics. Herein, the CUB and codon context patterns of protein-coding genes were systematically analyzed and compared between members of the two families. Overall, the nucleotide composition and CUB indices were found to differ between thermophiles and non-thermophiles. The thermophiles showed a higher G/C content in the codon base composition and tended to end with G/C compared to the non-thermophiles. Correlation analysis indicated significant associations between codon base composition and CUB indices. The results of the effective number of codons, parity-rule 2, neutral and correspondence analyses indicated that mutational pressure and natural selection primarily account for CUB in these cyanobacterial species, but the primary driving forces exhibit variation among genera. Moreover, the optimal codons identified based on relative synonymous codon usage values were found to differ among genera and even within genera. In addition, codon context pattern analysis revealed the specificity of the sequence context of start and stop codons among genera. Intriguingly, the clustering of codon context patterns appeared to be more related to thermotolerance than to phylogenomic relationships. In conclusion, this study facilitates the understanding of the characteristics and sources of variation of CUB and the evolution of the surveyed cyanobacterial clades with different thermotolerance and provides insights into their adaptation to different environments.

Persistence and neutrality in interacting replicator dynamics

We study the large-time behavior of an ensemble of entities obeying replicator-like stochastic dynamics with mean-field interactions as a model for a primordial ecology. We prove the propagation-of-chaos property and establish conditions for the strong persistence of the N-replicator system and the existence of invariant distributions for a class of associated McKean-Vlasov dynamics. In particular, our results show that, unlike typical models of neutral ecology, fitness equivalence does not need to be assumed but emerges as a condition for the persistence of the system. Further, neutrality is associated with a unique Dirichlet invariant probability measure. We illustrate our findings with some simple case studies, provide numerical results, and discuss our conclusions in the light of Neutral Theory in ecology.

Order of amino acid recruitment into the genetic code resolved by last universal common ancestor's protein domains

The current "consensus" order in which amino acids were added to the genetic code is based on potentially biased criteria, such as the absence of sulfur-containing amino acids from the Urey-Miller experiment which lacked sulfur. More broadly, abiotic abundance might not reflect biotic abundance in the organisms in which the genetic code evolved. Here, we instead identify which protein domains date to the last universal common ancestor (LUCA) and then infer the order of recruitment from deviations of their ancestrally reconstructed amino acid frequencies from the still-ancient post-LUCA controls. We find that smaller amino acids were added to the code earlier, with no additional predictive power in the previous consensus order. Metal-binding (cysteine and histidine) and sulfur-containing (cysteine and methionine) amino acids were added to the genetic code much earlier than previously thought. Methionine and histidine were added to the code earlier than expected from their molecular weights and glutamine later. Early methionine availability is compatible with inferred early use of S-adenosylmethionine and early histidine with its purine-like structure and the demand for metal binding. Even more ancient protein sequences-those that had already diversified into multiple distinct copies prior to LUCA-have significantly higher frequencies of aromatic amino acids (tryptophan, tyrosine, phenylalanine, and histidine) and lower frequencies of valine and glutamic acid than single-copy LUCA sequences. If at least some of these sequences predate the current code, then their distinct enrichment patterns provide hints about earlier, alternative genetic codes.

AARS Online: A collaborative database on the structure, function, and evolution of the aminoacyl-tRNA synthetases

The aminoacyl-tRNA synthetases (aaRS) are a large group of enzymes that implement the genetic code in all known biological systems. They attach amino acids to their cognate tRNAs, moonlight in various translational and non-translational activities beyond aminoacylation, and are linked to many genetic disorders. The aaRS have a subtle ontology characterized by structural and functional idiosyncrasies that vary from organism to organism, and protein to protein. Across the tree of life, the 22 coded amino acids are handled by 16 evolutionary families of Class I aaRS and 21 families of Class II aaRS. We introduce AARS Online, an interactive Wikipedia-like tool curated by an international consortium of field experts. This platform systematizes existing knowledge about the aaRS by showcasing a taxonomically diverse selection of aaRS sequences and structures. Through its graphical user interface, AARS Online facilitates a seamless exploration between protein sequence and structure, providing a friendly introduction to the material for non-experts and a useful resource for experts. Curated multiple sequence alignments can be extracted for downstream analyses. Accessible at www.aars.online, AARS Online is a free resource to delve into the world of the aaRS.

Comprehensive Analysis of Codon Usage Bias in Human Papillomavirus Type 51

Human papillomavirus type 51 (HPV-51) is associated with various cancers, including cervical cancer. Examining the codon usage bias of the organism can offer valuable insights into its evolutionary patterns and its relationship with the host. This study comprehensively analyzed codon usage bias in HPV-51 by examining 64 complete genome sequences sourced from the NCBI GenBank database. Our analysis revealed no noteworthy preference for codon usage in HPV-51 overall. However, there was a noticeable bias towards A/T-ending codons, accompanied by GC3s below 32%. Dinucleotide frequency analysis revealed reduced frequencies for ApA, CpG, and TpC dinucleotides, while CpA and TpG dinucleotides were more frequent than others. Relative Synonymous Codon Usage analysis revealed 30 favored codons, primarily concluding with A/T nucleotides. Further analysis using Parity Rule 2, Effective Number of Codons plot, and neutrality plot indicated a balance between mutational pressure and natural selection, with natural selection being the primary force shaping codon usage bias. The Isoacceptor tRNA Pool analysis indicates that HPV-51 has a higher translation efficiency within the human cellular translational system. Moreover, the Codon Adaptation Index and Relative Codon Deoptimization Index analyses suggested a moderate adaptation of HPV-51 to human codon preferences. Our discoveries offer valuable perspectives on how HPV-51 evolves and uses genetic codes, contributing to a deeper comprehension of its endurance and disease-causing potential.

Order of amino acid recruitment into the genetic code resolved by Last Universal Common Ancestor's protein domains

The current "consensus" order in which amino acids were added to the genetic code is based on potentially biased criteria, such as absence of sulfur-containing amino acids from the Urey-Miller experiment which lacked sulfur. More broadly, abiotic abundance might not reflect biotic abundance in the organisms in which the genetic code evolved. Here, we instead identify which protein domains date to the last universal common ancestor (LUCA), then infer the order of recruitment from deviations of their ancestrally reconstructed amino acid frequencies from the still-ancient post-LUCA controls. We find that smaller amino acids were added to the code earlier, with no additional predictive power in the previous "consensus" order. Metal-binding (cysteine and histidine) and sulfur-containing (cysteine and methionine) amino acids were added to the genetic code much earlier than previously thought. Methionine and histidine were added to the code earlier than expected from their molecular weights, and glutamine later. Early methionine availability is compatible with inferred early use of S-adenosylmethionine, and early histidine with its purine-like structure and the demand for metal-binding. Even more ancient protein sequences - those that had already diversified into multiple distinct copies prior to LUCA - have significantly higher frequencies of aromatic amino acids (tryptophan, tyrosine, phenylalanine and histidine), and lower frequencies of valine and glutamic acid than single copy LUCA sequences. If at least some of these sequences predate the current code, then their distinct enrichment patterns provide hints about earlier, alternative genetic codes.

Universal dynamical function behind all genetic codes: P-adic attractor dynamical model

The genetic code is a map which gives the correspondence between codons in DNA and amino acids. In the attractor dynamical model (ADM), genetic codes can be described as the sets of the cyclic attractors of discrete dynamical systems - the iterations of functions acting in the ring of 2-adic integers Z2. This ring arises from representation of nucleotides by binary vectors and hence codons by triples of binary vectors. We construct a Universal Function B such that the dynamical functions for all known genetic codes can be obtained from B by simple transformations on the set of codon cycles - the "Addition" and "Division" operations. ADM can be employed for study of phylogenetic dynamics of genetic codes. One can speculate that the "common ancestor genetic code" was caused by B. We remark that this function has 24 cyclic attractors which distribution coincides with the distribution for the hypothetical pre-LUCA code. This coupling of the Universal Function with the pre-LUCA code assigns the genetic codes evolution perspective to ADM. All genetic codes are generated from B through the special chains of the "Addition" and "Division" operations. The challenging problem is to assign the biological meaning to these mathematical operations.

Cracking the Code: Reprogramming the Genetic Script in Prokaryotes and Eukaryotes to Harness the Power of Noncanonical Amino Acids

Over 500 natural and synthetic amino acids have been genetically encoded in the last two decades. Incorporating these noncanonical amino acids into proteins enables many powerful applications, ranging from basic research to biotechnology, materials science, and medicine. However, major challenges remain to unleash the full potential of genetic code expansion across disciplines. Here, we provide an overview of diverse genetic code expansion methodologies and systems and their final applications in prokaryotes and eukaryotes, represented by Escherichia coli and mammalian cells as the main workhorse model systems. We highlight the power of how new technologies can be first established in simple and then transferred to more complex systems. For example, whole-genome engineering provides an excellent platform in bacteria for enabling transcript-specific genetic code expansion without off-targets in the transcriptome. In contrast, the complexity of a eukaryotic cell poses challenges that require entirely new approaches, such as striving toward establishing novel base pairs or generating orthogonally translating organelles within living cells. We connect the milestones in expanding the genetic code of living cells for encoding novel chemical functionalities to the most recent scientific discoveries, from optimizing the physicochemical properties of noncanonical amino acids to the technological advancements for their in vivo incorporation. This journey offers a glimpse into the promising developments in the years to come.

Examining the Effects of Temperature on the Evolution of Bacterial tRNA Pools

The genetic code consists of 61 codons coding for 20 amino acids. These codons are recognized by transfer RNAs (tRNAs) that bind to specific codons during protein synthesis. All organisms utilize less than all 61 possible anticodons due to base pair wobble: the ability to have a mismatch with a codon at its third nucleotide. Previous studies observed a correlation between the tRNA pool of bacteria and the temperature of their respective environments. However, it is unclear if these patterns represent biological adaptations to maintain the efficiency and accuracy of protein synthesis in different environments. A mechanistic mathematical model of mRNA translation is used to quantify the expected elongation rates and error rate for each codon based on an organism's tRNA pool. A comparative analysis across a range of bacteria that accounts for covariance due to shared ancestry is performed to quantify the impact of environmental temperature on the evolution of the tRNA pool. We find that thermophiles generally have more anticodons represented in their tRNA pool than mesophiles or psychrophiles. Based on our model, this increased diversity is expected to lead to increased missense errors. The implications of this for protein evolution in thermophiles are discussed.

Generation of genetic codes with 2-adic codon algebra and adaptive dynamics

This is a brief review on modeling genetic codes with the aid of 2-adic dynamical systems. In this model amino acids are encoded by the attractors of such dynamical systems. Each genetic code is coupled to the special class of 2-adic dynamics. We consider the discrete dynamical systems, These are the iterations of a function F:Z2→Z2, where Z2 is the ring of 2-adic numbers (2-adic tree). A genetic code is characterized by the set of attractors of a function belonging to the code generating functional class. The main mathematical problem is to reduce degeneration of dynamic representation and select the optimal generating function. Here optimality can be treated in many ways. One possibility is to consider the Lipschitz functions playing the crucial role in general theory of iterations. Then we minimize the Lip-constant. The main issue is to find the proper biological interpretation of code-functions. One can speculate that the evolution of the genetic codes can be described in information space of the nucleotide-strings endowed with ultrametric (treelike) geometry. A code-function is a fitness function; the solutions of the genetic code optimization problem are attractors of the code-function. We illustrate this approach by generation of the standard nuclear and (vertebrate) mitochondrial genetics codes.

Robust genetic codes enhance protein evolvability

The standard genetic code defines the rules of translation for nearly every life form on Earth. It also determines the amino acid changes accessible via single-nucleotide mutations, thus influencing protein evolvability-the ability of mutation to bring forth adaptive variation in protein function. One of the most striking features of the standard genetic code is its robustness to mutation, yet it remains an open question whether such robustness facilitates or frustrates protein evolvability. To answer this question, we use data from massively parallel sequence-to-function assays to construct and analyze 6 empirical adaptive landscapes under hundreds of thousands of rewired genetic codes, including those of codon compression schemes relevant to protein engineering and synthetic biology. We find that robust genetic codes tend to enhance protein evolvability by rendering smooth adaptive landscapes with few peaks, which are readily accessible from throughout sequence space. However, the standard genetic code is rarely exceptional in this regard, because many alternative codes render smoother landscapes than the standard code. By constructing low-dimensional visualizations of these landscapes, which each comprise more than 16 million mRNA sequences, we show that such alternative codes radically alter the topological features of the network of high-fitness genotypes. Whereas the genetic codes that optimize evolvability depend to some extent on the detailed relationship between amino acid sequence and protein function, we also uncover general design principles for engineering nonstandard genetic codes for enhanced and diminished evolvability, which may facilitate directed protein evolution experiments and the bio-containment of synthetic organisms, respectively.

Theories of the origin of the genetic code: Strong corroboration for the coevolution theory

I analyzed all the theories and models of the origin of the genetic code, and over the years, I have considered the main suggestions that could explain this origin. The conclusion of this analysis is that the coevolution theory of the origin of the genetic code is the theory that best captures the majority of observations concerning the organization of the genetic code. In other words, the biosynthetic relationships between amino acids would have heavily influenced the origin of the organization of the genetic code, as supported by the coevolution theory. Instead, the presence in the genetic code of physicochemical properties of amino acids, which have also been linked to the physicochemical properties of anticodons or codons or bases by stereochemical and physicochemical theories, would simply be the result of natural selection. More explicitly, I maintain that these correlations between codons, anticodons or bases and amino acids are in fact the result not of a real correlation between amino acids and codons, for example, but are only the effect of the intervention of natural selection. Specifically, in the genetic code table we expect, for example, that the most similar codons - that is, those that differ by only one base - will have more similar physicochemical properties. Therefore, the 64 codons of the genetic code table ordered in a certain way would also represent an ordering of some of their physicochemical properties. Now, a study aimed at clarifying which physicochemical property of amino acids has influenced the allocation of amino acids in the genetic code has established that the partition energy of amino acids has played a role decisive in this. Indeed, under some conditions, the genetic code was found to be approximately 98% optimized on its columns. In this same work, it was shown that this was most likely the result of the action of natural selection. If natural selection had truly allocated the amino acids in the genetic code in such a way that similar amino acids also have similar codons - this, not through a mechanism of physicochemical interaction between, for example, codons and amino acids - then it might turn out that even different physicochemical properties of codons (or anticodons or bases) show some correlation with the physicochemical properties of amino acids, simply because the partition energy of amino acids is correlated with other physicochemical properties of amino acids. It is very likely that this would inevitably lead to a correlation between codons (or anticodons or bases) and amino acids. In other words, since the codons (anticodons or bases) are ordered in the genetic code, that is to say, some of their physicochemical properties should also be ordered by a similar order, and given that the amino acids would also appear to have been ordered in the genetic code by selection natural, then it should inevitably turn out that there is a correlation between, for example, the hydrophobicity of anticodons and that of amino acids. Instead, the intervention of natural selection in organizing the genetic code would appear to be highly compatible with the main mechanism of structuring the genetic code as supported by the coevolution theory. This would make the coevolution theory the only plausible explanation for the origin of the genetic code.

Genetic code robustness and protein evolvability are correlated and protein-specific

The relationship between genetic code robustness and protein evolvability is unknown. A new study in PLOS Biology using in silico rewiring of genetic codes and functional protein data identified a positive correlation between code robustness and protein evolvability that is protein-specific.

Structural bioinformatics studies of glutamate transporters and their AlphaFold2 predicted water-soluble QTY variants and uncovering the natural mutations of L->Q, I->T, F->Y and Q->L, T->I and Y->F

Glutamate transporters play key roles in nervous physiology by modulating excitatory neurotransmitter levels, when malfunctioning, involving in a wide range of neurological and physiological disorders. However, integral transmembrane proteins including the glutamate transporters remain notoriously difficult to study, due to their localization within the cell membrane. Here we present the structural bioinformatics studies of glutamate transporters and their water-soluble variants generated through QTY-code, a protein design strategy based on systematic amino acid substitutions. These include 2 structures determined by X-ray crystallography, cryo-EM, and 6 predicted by AlphaFold2, and their predicted water-soluble QTY variants. In the native structures of glutamate transporters, transmembrane helices contain hydrophobic amino acids such as leucine (L), isoleucine (I), and phenylalanine (F). To design water-soluble variants, these hydrophobic amino acids are systematically replaced by hydrophilic amino acids, namely glutamine (Q), threonine (T) and tyrosine (Y). The QTY variants exhibited water-solubility, with four having identical isoelectric focusing points (pI) and the other four having very similar pI. We present the superposed structures of the native glutamate transporters and their water-soluble QTY variants. The superposed structures displayed remarkable similarity with RMSD 0.528Å-2.456Å, despite significant protein transmembrane sequence differences (41.1%->53.8%). Additionally, we examined the differences of hydrophobicity patches between the native glutamate transporters and their QTY variants. Upon closer inspection, we discovered multiple natural variations of L->Q, I->T, F->Y and Q->L, T->I, Y->F in these transporters. Some of these natural variations were benign and the remaining were reported in specific neurological disorders. We further investigated the characteristics of hydrophobic to hydrophilic substitutions in glutamate transporters, utilizing variant analysis and evolutionary profiling. Our structural bioinformatics studies not only provided insight into the differences between the hydrophobic helices and hydrophilic helices in the glutamate transporters, but they are also expected to stimulate further study of other water-soluble transmembrane proteins.

Structural bioinformatics studies of serotonin, dopamine and norepinephrine transporters and their AlphaFold2 predicted water-soluble QTY variants and uncovering the natural mutations of L->Q, I->T, F->Y and Q->L, T->I and Y->F

Monoamine transporters including transporters for serotonin, dopamine, and norepinephrine play key roles in monoaminergic synaptic signaling, involving in the molecular etiology of a wide range of neurological and physiological disorders. Despite being crucial drug targets, the study of transmembrane proteins remains challenging due to their localization within the cell membrane. To address this, we present the structural bioinformatics studies of 7 monoamine transporters and their water-soluble variants designed using the QTY code, by systematically replacing the hydrophobic amino acids leucine (L), valine (V), isoleucine (I) and phenylalanine (F) with hydrophilic amino acids (glutamine (Q), threonine (T) and tyrosine (Y). The resulting QTY variants, despite significant protein transmembrane sequence differences (44.27%-51.85%), showed similar isoelectric points (pI) and molecular weights. While their hydrophobic surfaces significantly reduced, this change resulted in a minimal structural alteration. Quantitatively, Alphafold2 predicted QTY variant structures displayed remarkable similarity with RMSD 0.492Å-1.619Å. Accompanied by the structural similarities of substituted amino acids in the context of 1.5Å electron density maps, our study revealed multiple QTY and reverse QTY variations in genomic databases. We further analyzed their phenotypical and topological characteristics. By extending evolutionary game theory to the molecular foundations of biology, we provided insights into the evolutionary dynamics of chemically distinct alpha-helices, their usage in different chemotherapeutic applications, and open possibilities of diagnostic medicine. Our study rationalizes that QTY variants of monoamine transporters may not only become distinct tools for medical, structural, and evolutionary research, but these transporters may also emerge as contemporary therapeutic targets, providing a new approach to treatment for several conditions.

A periodic table of monosaccharides

It is important to recognize the great diversity of monosaccharides commonly encountered in animals, plants, and microbes, as well as to organize them in a visually interesting style that also emphasizes their similarities and relatedness. This article discusses the nature of building blocks, monosaccharides, and monosaccharide derivatives-terms commonly used in discussing "glycomolecules" found in nature. To aid in awareness of monosaccharide diversity, here is presented a Periodic Table of Monosaccharides. The rationale is given for construction of the Table and the selection of 103 monosaccharides, which is largely based on those presented in the KEGG and SNFG websites of monosaccharides, and includes room to enlarge as new discoveries are made. The Table should have educational value and is intended to capture the attention and foster imagination of those not very familiar with glycosciences, and encourage researchers to delve deeper into this fascinating area.

Chapter 1: The Astrobiology Primer 3.0

The Astrobiology Primer 3.0 (ABP3.0) is a concise introduction to the field of astrobiology for students and others who are new to the field of astrobiology. It provides an entry into the broader materials in this supplementary issue of Astrobiology and an overview of the investigations and driving hypotheses that make up this interdisciplinary field. The content of this chapter was adapted from the other 10 articles in this supplementary issue and thus represents the contribution of all the authors who worked on these introductory articles. The content of this chapter is not exhaustive and represents the topics that the authors found to be the most important and compelling in a dynamic and changing field.

Chapter 9: Life as We Don't Know It

While Earth contains the only known example of life in the universe, it is possible that life elsewhere is fundamentally different from what we are familiar with. There is an increased recognition in the astrobiology community that the search for life should steer away from terran-specific biosignatures to those that are more inclusive to all life-forms. To start exploring the space of possibilities that life could occupy, we can try to dissociate life from the chemistry that composes it on Earth by envisioning how different life elsewhere could be in composition, lifestyle, medium, and form, and by exploring how the general principles that govern living systems on Earth might be found in different forms and environments across the Solar System. Exotic life-forms could exist on Mars or Venus, or icy moons like Europa and Enceladus, or even as a shadow biosphere on Earth. New perspectives on agnostic biosignature detection have also begun to emerge, allowing for a broader and more inclusive approach to seeking exotic life with unknown chemistry that is distinct from life as we know it on Earth.

Single cell transcriptomics reveals UAR codon reassignment in Palmarella salina (Metopida, Armophorea) and confirms Armophorida belongs to APM clade

Anaerobes have emerged in several major lineages of ciliates, but the number of independent transitions to anaerobiosis among ciliates is unknown. The APM clade (Armophorea, Muranotrichea, Parablepharismea) represents the largest clade of obligate anaerobes among ciliates and contains free-living marine and freshwater representatives as well as gut endobionts of animals. The evolution of APM group has only recently started getting attention, and our knowledge on its phylogeny and genetics is still limited to a fraction of taxa. While ciliates portray a wide array of alternatives to the standard genetic code across numerous classes, the APM ciliates were considered to be the largest group using exclusively standard nuclear genetic code. In this study, we present a pan-ciliate phylogenomic analysis with emphasis on the APM clade, bringing the first phylogenomic analysis of the family Tropidoatractidae (Armophorea) and confirming the position of Armophorida within Armophorea. We include five newly sequenced single cell transcriptomes from marine, freshwater, and endobiotic APM ciliates - Palmarella salina, Anteclevelandella constricta, Nyctotherus sp., Caenomorpha medusula, and Thigmothrix strigosa. We report the first discovery of an alternative nuclear genetic code among APM ciliates, used by Palmarella salina (Tropidoatractidae, Armophorea), but not by its close relative, Tropidoatractus sp., and provide a comparative analysis of stop codon identity and frequency indicating the precedency to the UAG codon loss/reassignment over the UAA codon reassignment in the specific ancestor of Palmarella. Comparative genomic and proteomic studies of this group may help explain the constraints that underlie UAR stop-to-sense reassignment, the most frequent type of alternative nuclear genetic code, not only in ciliates, but eukaryotes in general.

Enzymic recognition of amino acids drove the evolution of primordial genetic codes

How genetic information gained its exquisite control over chemical processes needed to build living cells remains an enigma. Today, the aminoacyl-tRNA synthetases (AARS) execute the genetic codes in all living systems. But how did the AARS that emerged over three billion years ago as low-specificity, protozymic forms then spawn the full range of highly-specific enzymes that distinguish between 22 diverse amino acids? A phylogenetic reconstruction of extant AARS genes, enhanced by analysing modular acquisitions, reveals six AARS with distinct bacterial, archaeal, eukaryotic, or organellar clades, resulting in a total of 36 families of AARS catalytic domains. Small structural modules that differentiate one AARS family from another played pivotal roles in discriminating between amino acid side chains, thereby expanding the genetic code and refining its precision. The resulting model shows a tendency for less elaborate enzymes, with simpler catalytic domains, to activate amino acids that were not synthesised until later in the evolution of the code. The most probable evolutionary route for an emergent amino acid type to establish a place in the code was by recruiting older, less specific AARS, rather than adapting contemporary lineages. This process, retrofunctionalisation, differs from previously described mechanisms through which amino acids would enter the code.

Evolution at the Origins of Life?

The role of evolutionary theory at the origin of life is an extensively debated topic. The origin and early development of life is usually separated into a prebiotic phase and a protocellular phase, ultimately leading to the Last Universal Common Ancestor. Most likely, the Last Universal Common Ancestor was subject to Darwinian evolution, but the question remains to what extent Darwinian evolution applies to the prebiotic and protocellular phases. In this review, we reflect on the current status of evolutionary theory in origins of life research by bringing together philosophy of science, evolutionary biology, and empirical research in the origins field. We explore the various ways in which evolutionary theory has been extended beyond biology; we look at how these extensions apply to the prebiotic development of (proto)metabolism; and we investigate how the terminology from evolutionary theory is currently being employed in state-of-the-art origins of life research. In doing so, we identify some of the current obstacles to an evolutionary account of the origins of life, as well as open up new avenues of research.

The Supersymmetry Genetic Code Table and Quadruplet Symmetries of DNA Molecules Are Unchangeable and Synchronized with Codon-Free Energy Mapping during Evolution

The Supersymmetry Genetic code (SSyGC) table is based on five physicochemical symmetries: (1) double mirror symmetry on the principle of the horizontal and vertical mirror symmetry axis between all bases (purines [A, G) and pyrimidines (U, C)] and (2) of bases in the form of codons; (3) direct-complement like codon/anticodon symmetry in the sixteen alternating boxes of the genetic code columns; (4) A + T-rich and C + G-rich alternate codons in the same row between both columns of the genetic code; (5) the same position between divided and undivided codon boxes in relation to horizontal mirror symmetry axis. The SSyGC table has a unique physicochemical purine-pyrimidine symmetry net which is as the core symmetry common for all, with more than thirty different nuclear and mitochondrial genetic codes. This net is present in the SSyGC table of all RNA and DNA living species. None of these symmetries are present in the Standard Genetic Code (SGC) table which is constructed on the alphabetic horizontal and vertical U-C-A-G order of bases. Here, we show that the free energy value of each codon incorporated as fundamentally mapping the "energy code" in the SSyGC table is compatible with mirror symmetry. On the other hand, in the SGC table, the same free energy values of codons are dispersed and a mirror symmetry between them is not recognizable. At the same time, the mirror symmetry of the SSyGC table and the DNA quadruplets together with our classification of codons/trinucleotides are perfectly imbedded in the mirror symmetry energy mapping of codons/trinucleotides and point out in favor of maintaining the integrity of the genetic code and DNA genome. We also argue that physicochemical symmetries of the SSyGC table in the manner of the purine-pyrimidine symmetry net, the quadruplet symmetry of DNA molecule, and the free energy of codons have remined unchanged during all of evolution. The unchangeable and universal symmetry properties of the genetic code, DNA molecules, and the energy code are decreasing disorder between codons/trinucleotides and shed a new light on evolution. Diversity in all living species on Earth is broad, but the symmetries of the Supersymmetry Genetic Code as the code of life and the DNA quadruplets related to the "energy code" are unique, unchangeable, and have the power of natural laws.

The 3 31 Nucleotide Minihelix tRNA Evolution Theorem and the Origin of Life

There are no theorems (proven theories) in the biological sciences. We propose that the 3 31 nt minihelix tRNA evolution theorem be universally accepted as one. The 3 31 nt minihelix theorem completely describes the evolution of type I and type II tRNAs from ordered precursors (RNA repeats and inverted repeats). Despite the diversification of tRNAome sequences, statistical tests overwhelmingly support the theorem. Furthermore, the theorem relates the dominant pathway for the origin of life on Earth, specifically, how tRNAomes and the genetic code may have coevolved. Alternate models for tRNA evolution (i.e., 2 minihelix, convergent and accretion models) are falsified. In the context of the pre-life world, tRNA was a molecule that, via mutation, could modify anticodon sequences and teach itself to code. Based on the tRNA sequence, we relate the clearest history to date of the chemical evolution of life. From analysis of tRNA evolution, ribozyme-mediated RNA ligation was a primary driving force in the evolution of complexity during the pre-life-to-life transition. TRNA formed the core for the evolution of living systems on Earth.

The spiderweb of error-detecting codes in the genetic information

Nature possesses inherent mechanisms for error detection and correction during the translation of genetic information, as demonstrated by the discovery of a self-complementary circular C3-code called X0 in various organisms such as bacteria, eukaryotes, plasmids, and viruses (Arquès and Michel, 1996; Michel, 2015, 2017). Since then, extensive research has focused on circular codes, which are believed to be remnants of ancient comma-free codes. These codes can be regarded as an additional genetic code specifically optimized for detecting and preserving the proper reading frame in protein-coding sequences. A study by Fimmel et al. in 2014 identified that a total of 216 maximal self-complementary C3-codes can be grouped into 27 equivalence classes with eight codes in each class. In this work, we study how the 27 equivalence classes are related to each other. While the codes in each equivalence class obtained by Fimmel et al. in 2014 are permutations of each other, i.e. one code can be obtained from the other by applying a permutation of the bases, it has not been clear how the equvalence classes are connected. We show that there is an ordering of the equivalence classes such that one gets from one class to the next one by substituting only one pair of codon/anticodon in the corresponding codes, i.e. the corresponding codes have a maximal intersection of 18 codons. To perform this analysis, we define two graphs, G216 and G27, whose vertices are, respectively, all 216 maximal self-complementary C3-codes and 27 equivalence classes. Several properties of the graphs are obtained. Most surprisingly, it turns out that G27 contains Hamiltonian paths of length 27. This fact ultimately leads to a representation of the set of all 216 maximal self-complementary C3-codes as a kind of spider web. Finally, we define dinucleotide cuts of such codes by projecting each codon to its first two bases and show that the paths of lengths 27 in G216 can even be chosen so that all the codes contain a special subset of dinucleotides defined by Rumer's roots. These observations raise a lot of new questions about the biological function of such structures.

The Genetic Code Assembles via Division and Fusion, Basic Cellular Events

Standard Genetic Code (SGC) evolution is quantitatively modeled in up to 2000 independent coding 'environments'. Environments host multiple codes that may fuse or divide, with division yielding identical descendants. Code division may be selected-sophisticated gene products could be required for an orderly separation that preserves the coding. Several unforeseen results emerge: more rapid evolution requires unselective code division rather than its selective form. Combining selective and unselective code division, with/without code fusion, with/without independent environmental coding tables, and with/without wobble defines 25 = 32 possible pathways for SGC evolution. These 32 possible histories are compared, specifically, for evolutionary speed and code accuracy. Pathways differ greatly, for example, by ≈300-fold in time to evolve SGC-like codes. Eight of thirty-two pathways employing code division evolve quickly. Four of these eight that combine fusion and division also unite speed and accuracy. The two most precise, swiftest paths; thus the most likely routes to the SGC are similar, differing only in fusion with independent environmental codes. Code division instead of fusion with unrelated codes implies that exterior codes can be dispensable. Instead, a single ancestral code that divides and fuses can initiate fully encoded peptide biosynthesis. Division and fusion create a 'crescendo of competent coding', facilitating the search for the SGC and also assisting the advent of otherwise uniformly disfavored wobble coding. Code fusion can unite multiple codon assignment mechanisms. However, via code division and fusion, an SGC can emerge from a single primary origin via familiar cellular events.

Determining the effects of temperature on the evolution of bacterial tRNA pools

The genetic code consists of 61 codon coding for 20 amino acids. These codons are recognized by transfer RNAs (tRNA) that bind to specific codons during protein synthesis. Most organisms utilize less than all 61 possible anticodons due to base pair wobble: the ability to have a mismatch with a codon at its third nucleotide. Previous studies observed a correlation between the tRNA pool of bacteria and the temperature of their respective environments. However, it is unclear if these patterns represent biological adaptations to maintain the efficiency and accuracy of protein synthesis in different environments. A mechanistic mathematical model of mRNA translation is used to quantify the expected elongation rates and error rate for each codon based on an organism's tRNA pool. A comparative analysis across a range of bacteria that accounts for covariance due to shared ancestry is performed to quantify the impact of environmental temperature on the evolution of the tRNA pool. We find that thermophiles generally have more anticodons represented in their tRNA pool than mesophiles or psychrophiles. Based on our model, this increased diversity is expected to lead to increased missense errors. The implications of this for protein evolution in thermophiles are discussed.

Identification of a non-canonical ciliate nuclear genetic code where UAA and UAG code for different amino acids

The genetic code is one of the most highly conserved features across life. Only a few lineages have deviated from the "universal" genetic code. Amongst the few variants of the genetic code reported to date, the codons UAA and UAG virtually always have the same translation, suggesting that their evolution is coupled. Here, we report the genome and transcriptome sequencing of a novel uncultured ciliate, belonging to the Oligohymenophorea class, where the translation of the UAA and UAG stop codons have changed to specify different amino acids. Genomic and transcriptomic analyses revealed that UAA has been reassigned to encode lysine, while UAG has been reassigned to encode glutamic acid. We identified multiple suppressor tRNA genes with anticodons complementary to the reassigned codons. We show that the retained UGA stop codon is enriched in the 3'UTR immediately downstream of the coding region of genes, suggesting that there is functional drive to maintain tandem stop codons. Using a phylogenomics approach, we reconstructed the ciliate phylogeny and mapped genetic code changes, highlighting the remarkable number of independent genetic code changes within the Ciliophora group of protists. According to our knowledge, this is the first report of a genetic code variant where UAA and UAG encode different amino acids.

Prebiotic Synthesis of Aspartate Using Life's Metabolism as a Guide

A protometabolic approach to the origins of life assumes that the conserved biochemistry of metabolism has direct continuity with prebiotic chemistry. One of the most important amino acids in modern biology is aspartic acid, serving as a nodal metabolite for the synthesis of many other essential biomolecules. Aspartate's prebiotic synthesis is complicated by the instability of its precursor, oxaloacetate. In this paper, we show that the use of the biologically relevant cofactor pyridoxamine, supported by metal ion catalysis, is sufficiently fast to offset oxaloacetate's degradation. Cu²⁺-catalysed transamination of oxaloacetate by pyridoxamine achieves around a 5% yield within 1 h, and can operate across a broad range of pH, temperature, and pressure. In addition, the synthesis of the downstream product β-alanine may also take place in the same reaction system at very low yields, directly mimicking an archaeal synthesis route. Amino group transfer supported by pyridoxal is shown to take place from aspartate to alanine, but the reverse reaction (alanine to aspartate) shows a poor yield. Overall, our results show that the nodal metabolite aspartate and related amino acids can indeed be synthesised via protometabolic pathways that foreshadow modern metabolism in the presence of the simple cofactor pyridoxamine and metal ions.

Coding From Binding? Molecular Interactions at the Heart of Translation

The mechanism and the evolution of DNA replication and transcription, the key elements of the central dogma of biology, are fundamentally well explained by the physicochemical complementarity between strands of nucleic acids. However, the determinants that have shaped the third part of the dogma-the process of biological translation and the universal genetic code-remain unclear. We review and seek parallels between different proposals that view the evolution of translation through the prism of weak, noncovalent interactions between biological macromolecules. In particular, we focus on a recent proposal that there exists a hitherto unrecognized complementarity at the heart of biology, that between messenger RNA coding regions and the proteins that they encode, especially if the two are unstructured. Reflecting the idea that the genetic code evolved from intrinsic binding propensities between nucleotides and amino acids, this proposal promises to forge a link between the distant past and the present of biological systems.

Biophysical Interactions Underpin the Emergence of Information in the Genetic Code

The genetic code conceals a 'code within the codons', which hints at biophysical interactions between amino acids and their cognate nucleotides. Yet, research over decades has failed to corroborate systematic biophysical interactions across the code. Using molecular dynamics simulations and NMR, we have analysed interactions between the 20 standard proteinogenic amino acids and 4 RNA mononucleotides in 3 charge states. Our simulations show that 50% of amino acids bind best with their anticodonic middle base in the -1 charge state common to the backbone of RNA, while 95% of amino acids interact most strongly with at least 1 of their codonic or anticodonic bases. Preference for the cognate anticodonic middle base was greater than 99% of randomised assignments. We verify a selection of our results using NMR, and highlight challenges with both techniques for interrogating large numbers of weak interactions. Finally, we extend our simulations to a range of amino acids and dinucleotides, and corroborate similar preferences for cognate nucleotides. Despite some discrepancies between the predicted patterns and those observed in biology, the existence of weak stereochemical interactions means that random RNA sequences could template non-random peptides. This offers a compelling explanation for the emergence of genetic information in biology.

Rare-event sampling analysis uncovers the fitness landscape of the genetic code

The genetic code refers to a rule that maps 64 codons to 20 amino acids. Nearly all organisms, with few exceptions, share the same genetic code, the standard genetic code (SGC). While it remains unclear why this universal code has arisen and been maintained during evolution, it may have been preserved under selection pressure. Theoretical studies comparing the SGC and numerically created hypothetical random genetic codes have suggested that the SGC has been subject to strong selection pressure for being robust against translation errors. However, these prior studies have searched for random genetic codes in only a small subspace of the possible code space due to limitations in computation time. Thus, how the genetic code has evolved, and the characteristics of the genetic code fitness landscape, remain unclear. By applying multicanonical Monte Carlo, an efficient rare-event sampling method, we efficiently sampled random codes from a much broader random ensemble of genetic codes than in previous studies, estimating that only one out of every 1020 random codes is more robust than the SGC. This estimate is significantly smaller than the previous estimate, one in a million. We also characterized the fitness landscape of the genetic code that has four major fitness peaks, one of which includes the SGC. Furthermore, genetic algorithm analysis revealed that evolution under such a multi-peaked fitness landscape could be strongly biased toward a narrow peak, in an evolutionary path-dependent manner.

Reflections on the Origin and Early Evolution of the Genetic Code

Examination of the genetic code (GeCo) reveals that amino acids coded by (A/U) codons display a large functional spectrum and bind RNA whereas, except for Arg, those coded by (G/C) codons do not. From a stereochemical viewpoint, the clear preference for (A/U)-rich codons to be located at the GeCo half blocks suggests they were specifically determined. Conversely, the overall lower affinity of cognate amino acids for their (G/C)-rich anticodons points to their late arrival to the GeCo. It is proposed that i) initially the code was composed of the eight (A/U) codons; ii) these codons were duplicated when G/C nucleotides were added to their wobble positions, and three new codons with G/C in their first position were incorporated; and iii) a combination of A/U and G/C nucleotides progressively generated the remaining codons.

Functional variation in the non-coding genome: molecular implications for food security

The growing world population, in combination with the anticipated effects of climate change, is pressuring food security. Plants display an impressive arsenal of cellular mechanisms conferring resilience to adverse environmental conditions, and humans rely on these mechanisms for stable food production. The elucidation of the molecular basis of the mechanisms used by plants to achieve resilience promises knowledge-based approaches to enhance food security. DNA sequence polymorphisms can reveal genomic regions that are linked to beneficial traits of plants. However, our ability to interpret how a given DNA sequence polymorphism confers a fitness advantage at the molecular level often remains poor. A key factor is that these polymorphisms largely localize to the enigmatic non-coding genome. Here, we review the functional impact of sequence variations in the non-coding genome on plant biology in the context of crop breeding and agricultural traits. We focus on examples of non-coding with particularly convincing functional support. Our survey combines findings that are consistent with the view that the non-coding genome contributes to cellular mechanisms assisting many plant traits. Understanding how DNA sequence polymorphisms in the non-coding genome shape plant traits at the molecular level offers a largely unexplored reservoir of solutions to address future challenges in plant growth and resilience.

Anticodon sequence determines the impact of mistranslating tRNAAla variants

Transfer RNAs (tRNAs) maintain translation fidelity through accurate charging by their cognate aminoacyl-tRNA synthetase and codon:anticodon base pairing with the mRNA at the ribosome. Mistranslation occurs when an amino acid not specified by the genetic message is incorporated into proteins and has applications in biotechnology, therapeutics and is relevant to disease. Since the alanyl-tRNA synthetase uniquely recognizes a G3:U70 base pair in tRNAAla and the anticodon plays no role in charging, tRNAAla variants with anticodon mutations have the potential to mis-incorporate alanine. Here, we characterize the impact of the 60 non-alanine tRNAAla anticodon variants on the growth of Saccharomyces cerevisiae. Overall, 36 tRNAAla anticodon variants decreased growth in single- or multi-copy. Mass spectrometry analysis of the cellular proteome revealed that 52 of 57 anticodon variants, not decoding alanine or stop codons, induced mistranslation when on single-copy plasmids. Variants with G/C-rich anticodons resulted in larger growth deficits than A/U-rich variants. In most instances, synonymous anticodon variants impact growth differently, with anticodons containing U at base 34 being the least impactful. For anticodons generating the same amino acid substitution, reduced growth generally correlated with the abundance of detected mistranslation events. Differences in decoding specificity, even between synonymous anticodons, resulted in each tRNAAla variant mistranslating unique sets of peptides and proteins. We suggest that these differences in decoding specificity are also important in determining the impact of tRNAAla anticodon variants.

Short tRNA anticodon stem and mutant eRF1 allow stop codon reassignment

Cognate tRNAs deliver specific amino acids to translating ribosomes according to the standard genetic code, and three codons with no cognate tRNAs serve as stop codons. Some protists have reassigned all stop codons as sense codons, neglecting this fundamental principle^1-4. Here we analyse the in-frame stop codons in 7,259 predicted protein-coding genes of a previously undescribed trypanosomatid, Blastocrithidia nonstop. We reveal that in this species in-frame stop codons are underrepresented in genes expressed at high levels and that UAA serves as the only termination codon. Whereas new tRNAs^Glu fully cognate to UAG and UAA evolved to reassign these stop codons, the UGA reassignment followed a different path through shortening the anticodon stem of tRNA^Trp_CCA from five to four base pairs (bp). The canonical 5-bp tRNA^Trp recognizes UGG as dictated by the genetic code, whereas its shortened 4-bp variant incorporates tryptophan also into in-frame UGA. Mimicking this evolutionary twist by engineering both variants from B. nonstop, Trypanosoma brucei and Saccharomyces cerevisiae and expressing them in the last two species, we recorded a significantly higher readthrough for all 4-bp variants. Furthermore, a gene encoding B. nonstop release factor 1 acquired a mutation that specifically restricts UGA recognition, robustly potentiating the UGA reassignment. Virtually the same strategy has been adopted by the ciliate Condylostoma magnum. Hence, we describe a previously unknown, universal mechanism that has been exploited in unrelated eukaryotes with reassigned stop codons.

T cell receptor convergence is an indicator of antigen-specific T cell response in cancer immunotherapies

T cells are potent at eliminating pathogens and playing a crucial role in the adaptive immune response. T cell receptor (TCR) convergence describes T cells that share identical TCRs with the same amino acid sequences but have different DNA sequences due to codon degeneracy. We conducted a systematic investigation of TCR convergence using single-cell immune profiling and bulk TCRβ-sequence (TCR-seq) data obtained from both mouse and human samples and uncovered a strong link between antigen-specificity and convergence. This association was stronger than T cell expansion, a putative indicator of antigen-specific T cells. By using flow-sorted tetramer+ single T cell data, we discovered that convergent T cells were enriched for a neoantigen-specific CD8+ effector phenotype in the tumor microenvironment. Moreover, TCR convergence demonstrated better prediction accuracy for immunotherapy response than the existing TCR repertoire indexes. In conclusion, convergent T cells are likely to be antigen-specific and might be a novel prognostic biomarker for anti-cancer immunotherapy.

Refactored genetic codes enable bidirectional genetic isolation

The near-universal genetic code defines the correspondence between codons in genes and amino acids in proteins. We refactored the structure of the genetic code in Escherichia coli and created orthogonal genetic codes that restrict the escape of synthetic genetic information into natural life. We developed orthogonal and mutually orthogonal horizontal gene transfer systems, which permit the transfer of genetic information between organisms that use the same genetic code but restrict the transfer of genetic information between organisms that use different genetic codes. Moreover, we showed that locking refactored codes into synthetic organisms completely blocks invasion by mobile genetic elements, including viruses, which carry their own translation factors and successfully invade organisms with canonical and compressed genetic codes.

A biophysical basis for the emergence of the genetic code in protocells

The origin of the genetic code is an abiding mystery in biology. Hints of a 'code within the codons' suggest biophysical interactions, but these patterns have resisted interpretation. Here, we present a new framework, grounded in the autotrophic growth of protocells from CO₂ and H₂. Recent work suggests that the universal core of metabolism recapitulates a thermodynamically favoured protometabolism right up to nucleotide synthesis. Considering the genetic code in relation to an extended protometabolism allows us to predict most codon assignments. We show that the first letter of the codon corresponds to the distance from CO₂ fixation, with amino acids encoded by the purines (G followed by A) being closest to CO₂ fixation. These associations suggest a purine-rich early metabolism with a restricted pool of amino acids. The second position of the anticodon corresponds to the hydrophobicity of the amino acid encoded. We combine multiple measures of hydrophobicity to show that this correlation holds strongly for early amino acids but is weaker for later species. Finally, we demonstrate that redundancy at the third position is not randomly distributed around the code: non-redundant amino acids can be assigned based on size, specifically length. We attribute this to additional stereochemical interactions at the anticodon. These rules imply an iterative expansion of the genetic code over time with codon assignments depending on both distance from CO₂ and biophysical interactions between nucleotide sequences and amino acids. In this way the earliest RNA polymers could produce non-random peptide sequences with selectable functions in autotrophic protocells.

Phenotypic mutations contribute to protein diversity and shape protein evolution

Errors in DNA replication generate genetic mutations, while errors in transcription and translation lead to phenotypic mutations. Phenotypic mutations are orders of magnitude more frequent than genetic ones, yet they are less understood. Here, we review the types of phenotypic mutations, their quantifications, and their role in protein evolution and disease. The diversity generated by phenotypic mutation can facilitate adaptive evolution. Indeed, phenotypic mutations, such as ribosomal frameshift and stop codon readthrough, sometimes serve to regulate protein expression and function. Phenotypic mutations have often been linked to fitness decrease and diseases. Thus, understanding the protein heterogeneity and phenotypic diversity caused by phenotypic mutations will advance our understanding of protein evolution and have implications on human health and diseases.

Evolutionary action of mutations reveals antimicrobial resistance genes in Escherichia coli

Since antibiotic development lags, we search for potential drug targets through directed evolution experiments. A challenge is that many resistance genes hide in a noisy mutational background as mutator clones emerge in the adaptive population. Here, to overcome this noise, we quantify the impact of mutations through evolutionary action (EA). After sequencing ciprofloxacin or colistin resistance strains grown under different mutational regimes, we find that an elevated sum of the evolutionary action of mutations in a gene identifies known resistance drivers. This EA integration approach also suggests new antibiotic resistance genes which are then shown to provide a fitness advantage in competition experiments. Moreover, EA integration analysis of clinical and environmental isolates of antibiotic resistant of E. coli identifies gene drivers of resistance where a standard approach fails. Together these results inform the genetic basis of de novo colistin resistance and support the robust discovery of phenotype-driving genes via the evolutionary action of genetic perturbations in fitness landscapes.

Mitogenomics and mitochondrial gene phylogeny decipher the evolution of Saccharomycotina yeasts

Saccharomycotina yeasts belong to diverse clades within the kingdom of fungi and are important to human everyday life. This work investigates the evolutionary relationships among these yeasts from a mitochondrial (mt) genomic perspective. A comparative study of 155 yeast mt genomes representing all major phylogenetic lineages of Saccharomycotina was performed, including genome size and content variability, intron and intergenic regions’ diversity, genetic code alterations, and syntenic variation. Findings from this study suggest that mt genome size diversity is the result of a ceaseless random process, mainly based on genetic recombination and intron mobility. Gene order analysis revealed conserved syntenic units and many occurring rearrangements, which can be correlated with major evolutionary events as shown by the phylogenetic analysis of the concatenated mt protein matrix. For the first time, molecular dating indicated a slower mt genome divergence rate in the early stages of yeast evolution, in contrast with a faster rate in the late evolutionary stages, compared to their nuclear time divergence. Genetic code reassignments of mt genomes are a perpetual process happening in many different parallel evolutionary steps throughout the evolution of Saccharomycotina. Overall, this work shows that phylogenetic studies based on the mt genome of yeasts highlight major evolutionary events.

The Origin of Genetic Code and Translation in the Framework of Current Concepts on the Origin of Life

The origin of genetic code and translation system is probably the central and most difficult problem in the investigations on the origin of life and one of the most complex problems in the evolutionary biology in general. There are multiple hypotheses on the emergence and development of existing genetic systems that propose the mechanisms for the origin and early evolution of genetic code, as well as for the emergence of replication and translation. Here, we discuss the most well-known of these hypotheses, although none of them provides a description of the early evolution of genetic systems without gaps and assumptions. The RNA world hypothesis is a currently prevailing scientific idea on the early evolution of biological and pre-biological structures, the main advantage of which is the assumption that RNAs as the first living systems were self-sufficient, i.e., capable of functioning as both catalysts and templates. However, this hypothesis has also significant limitations. In particular, no ribozymes with processive polymerase activity have been yet discovered or synthesized. Taking into account the mutual need of proteins and nucleic acids in each other in the current world, many authors propose the early evolution scenarios based on the co-evolution of these two classes of organic molecules. They postulate that the emergence of translation was necessary for the replication of nucleic acids, in contrast to the RNA world hypothesis, according to which the emergence of translation was preceded by the era of self-replicating RNAs. Although such scenarios are less parsimonious from the evolutionary point of view, since they require simultaneous emergence and evolution of two classes of organic molecules, as well as the emergence of synchronized replication and translation, their major advantage is that they explain the development of processive and much more accurate protein-dependent replication.

Multidimensional Phylogenetic Metrics Identify Class I Aminoacyl-tRNA Synthetase Evolutionary Mosaicity and Inter-Modular Coupling

The role of aminoacyl-tRNA synthetases (aaRS) in the emergence and evolution of genetic coding poses challenging questions concerning their provenance. We seek evidence about their ancestry from curated structure-based multiple sequence alignments of a structurally invariant “scaffold” shared by all 10 canonical Class I aaRS. Three uncorrelated phylogenetic metrics—mutation frequency, its uniformity, and row-by-row cladistic congruence—imply that the Class I scaffold is a mosaic assembled from successive genetic sources. Metrics for different modules vary in accordance with their presumed functionality. Sequences derived from the ATP– and amino acid– binding sites exhibit specific two-way coupling to those derived from Connecting Peptide 1, a third module whose metrics suggest later acquisition. The data help validate: (i) experimental fragmentations of the canonical Class I structure into three partitions that retain catalytic activities in proportion to their length; and (ii) evidence that the ancestral Class I aaRS gene also encoded a Class II ancestor in frame on the opposite strand. A 46-residue Class I “protozyme” roots the Class I tree prior to the adaptive radiation of the Rossmann dinucleotide binding fold that refined substrate discrimination. Such rooting implies near simultaneous emergence of genetic coding and the origin of the proteome, resolving a conundrum posed by previous inferences that Class I aaRS evolved after the genetic code had been implemented in an RNA world. Further, pinpointing discontinuous enhancements of aaRS fidelity establishes a timeline for the growth of coding from a binary amino acid alphabet.