Nuclear genetic codes with a different meaning of the UAG and the UAA codon

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.

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.

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.

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.

The Influence of the Selection at the Amino Acid Level on Synonymous Codon Usage from the Viewpoint of Alternative Genetic Codes

Synonymous codon usage can be influenced by mutations and/or selection, e.g., for speed of protein translation and correct folding. However, this codon bias can also be affected by a general selection at the amino acid level due to differences in the acceptance of the loss and generation of these codons. To assess the importance of this effect, we constructed a mutation-selection model model, in which we generated almost 90,000 stationary nucleotide distributions produced by mutational processes and applied a selection based on differences in physicochemical properties of amino acids. Under these conditions, we calculated the usage of fourfold degenerated (4FD) codons and compared it with the usage characteristic of the pure mutations. We considered both the standard genetic code (SGC) and alternative genetic codes (AGCs). The analyses showed that a majority of AGCs produced a greater 4FD codon bias than the SGC. The mutations producing more thymine or adenine than guanine and cytosine increased the differences in usage. On the other hand, the mutational pressures generating a lot of cytosine or guanine with a low content of adenine and thymine decreased this bias because the nucleotide content of most 4FD codons stayed in the compositional equilibrium with these pressures. The comparison of the theoretical results with those for real protein coding sequences showed that the influence of selection at the amino acid level on the synonymous codon usage cannot be neglected. The analyses indicate that the effect of amino acid selection cannot be disregarded and that it can interfere with other selection factors influencing codon usage, especially in AT-rich genomes, in which AGCs are usually used.

Experimental validation that human microbiome phages use alternative genetic coding

Previous bioinformatic analyses of metagenomic data have indicated that bacteriophages can use genetic codes different from those of their host bacteria. In particular, reassignment of stop codon TAG to glutamine (a variation known as 'genetic code 15') has been predicted. Here, we use LC-MS/MS-based metaproteomics of human fecal samples to provide experimental evidence of the use of genetic code 15 in two crAss-like phages. Furthermore, the proteomic data from several phage structural proteins supports the reassignment of the TAG stop codon to glutamine late in the phage infection cycle. Thus, our work experimentally validates the expression of genetic code 15 in human microbiome phages.

Evidence for an Independent Hydrogenosome-to-Mitosome Transition in the CL3 Lineage of Fornicates

Fornicata, a lineage of a broader and ancient anaerobic eukaryotic clade Metamonada, contains diverse taxa that are ideally suited for evolutionary studies addressing various fundamental biological questions, such as the evolutionary trajectory of mitochondrion-related organelles (MROs), the transition between free-living and endobiotic lifestyles, and the derivation of alternative genetic codes. To this end, we conducted detailed microscopic and transcriptome analyses in a poorly documented strain of an anaerobic free-living marine flagellate, PCS, in the so-called CL3 fornicate lineage. Fortuitously, we discovered that the original culture contained two morphologically similar and closely related CL3 representatives, which doubles the taxon representation within this lineage. We obtained a monoeukaryotic culture of one of them and formally describe it as a new member of the family Caviomonadidae, Euthynema mutabile gen. et sp. nov. In contrast to previously studied caviomonads, the endobiotic Caviomonas mobilis and Iotanema spirale, E. mutabile possesses an ultrastructurally discernible MRO. We sequenced and assembled the transcriptome of E. mutabile, and by sequence subtraction, obtained transcriptome data from the other CL3 clade representative present in the original PCS culture, denoted PCS-ghost. Transcriptome analyses showed that the reassignment of only one of the UAR stop codons to encode Gln previously reported from I. spirale does not extend to its free-living relatives and is likely due to a unique amino acid substitution in I. spirale’s eRF1 protein domain responsible for termination codon recognition. The backbone fornicate phylogeny was robustly resolved in a phylogenomic analysis, with the CL3 clade amongst the earliest branching lineages. Metabolic and MRO functional reconstructions of CL3 clade members revealed that all three, including I. spirale, encode homologs of key components of the mitochondrial protein import apparatus and the ISC pathway, indicating the presence of a MRO in all of them. In silico evidence indicates that the organelles of E. mutabile and PCS-ghost host ATP and H₂ production, unlike the cryptic MRO of I. spirale. These data suggest that the CL3 clade has experienced a hydrogenosome-to-mitosome transition independent from that previously documented for the lineage leading to Giardia.

Nonsense Mutations in Eukaryotes

Nonsense mutations are a type of mutations which results in a premature termination codon occurrence. In general, these mutations have been considered to be among the most harmful ones which lead to premature protein translation termination and result in shortened nonfunctional polypeptide. However, there is evidence that not all nonsense mutations are harmful as well as some molecular mechanisms exist which allow to avoid pathogenic effects of these mutations. This review addresses relevant information on nonsense mutations in eukaryotic genomes, characteristics of these mutations, and different molecular mechanisms preventing or mitigating harmful effects thereof.

Codon usage bias

Codon usage bias is the preferential or non-random use of synonymous codons, a ubiquitous phenomenon observed in bacteria, plants and animals. Different species have consistent and characteristic codon biases. Codon bias varies not only with species, family or group within kingdom, but also between the genes within an organism. Codon usage bias has evolved through mutation, natural selection, and genetic drift in various organisms. Genome composition, GC content, expression level and length of genes, position and context of codons in the genes, recombination rates, mRNA folding, and tRNA abundance and interactions are some factors influencing codon bias. The factors shaping codon bias may also be involved in evolution of the universal genetic code. Codon-usage bias is critical factor determining gene expression and cellular function by influencing diverse processes such as RNA processing, protein translation and protein folding. Codon usage bias reflects the origin, mutation patterns and evolution of the species or genes. Investigations of codon bias patterns in genomes can reveal phylogenetic relationships between organisms, horizontal gene transfers, molecular evolution of genes and identify selective forces that drive their evolution. Most important application of codon bias analysis is in the design of transgenes, to increase gene expression levels through codon optimization, for development of transgenic crops. The review gives an overview of deviations of genetic code, factors influencing codon usage or bias, codon usage bias of nuclear and organellar genes, computational methods to determine codon usage and the significance as well as applications of codon usage analysis in biological research, with emphasis on plants.

Evolution of the genetic code in the mitochondria of Labyrinthulea (Stramenopiles)

Mitochondrial translation often exhibits departures from the standard genetic code, but the full spectrum of these changes has certainly not yet been described and the molecular mechanisms behind the changes in codon meaning are rarely studied. Here we report a detailed analysis of the mitochondrial genetic code in the stramenopile group Labyrinthulea (Labyrinthulomycetes) and their relatives. In the genus Aplanochytrium, UAG is not a termination codon but encodes tyrosine, in contrast to the unaffected meaning of the UAA codon. This change is evolutionarily independent of the reassignment of both UAG and UAA as tyrosine codons recently reported from two uncultivated labyrinthuleans (S2 and S4), which we show are not thraustochytrids as proposed before, but represent the clade LAB14 previously recognised in environmental 18S rRNA gene surveys. We provide rigorous evidence that the UUA codon in the mitochondria of all labyrinthuleans serves as a termination codon instead of encoding leucine, and propose that a sense-to-stop reassignment has also affected the AGG and AGA codons in the LAB14 clade. The distribution of the different forms of sense-to-stop and stop-to-sense reassignments correlates with specific modifications of the mitochondrial release factor mtRF2a in different subsets of labyrinthuleans, and with the unprecedented loss of mtRF1a in Aplanochytrium and perhaps also in the LAB14 clade, pointing towards a possible mechanistic basis of the code changes observed. Curiously, we show that labyrinthulean mitochondria also exhibit a sense-to-sense codon reassignment, manifested as AUA encoding methionine instead of isoleucine. Furthermore, we show that this change evolved independently in the uncultivated stramenopile lineage MAST8b, together with the reassignment of the AGR codons from arginine to serine. Altogether, our study has uncovered novel variants of the mitochondrial genetic code and previously unknown modifications of the mitochondrial translation machinery, further enriching our understanding of the rules governing the evolution of one of the central molecular process in the cell.

Evolution and Unprecedented Variants of the Mitochondrial Genetic Code in a Lineage of Green Algae

Mitochondria of diverse eukaryotes have evolved various departures from the standard genetic code, but the breadth of possible modifications and their phylogenetic distribution are known only incompletely. Furthermore, it is possible that some codon reassignments in previously sequenced mitogenomes have been missed, resulting in inaccurate protein sequences in databases. Here we show, considering the distribution of codons at conserved amino acid positions in mitogenome-encoded proteins, that mitochondria of the green algal order Sphaeropleales exhibit a diversity of codon reassignments, including previously missed ones and some that are unprecedented in any translation system examined so far, necessitating redefinition of existing translation tables and creating at least seven new ones. We resolve a previous controversy concerning the meaning the UAG codon in Hydrodictyaceae, which beyond any doubt encodes alanine. We further demonstrate that AGG, sometimes together with AGA, encodes alanine instead of arginine in diverse sphaeroplealeans. Further newly detected changes include Arg-to-Met reassignment of the AGG codon and Arg-to-Leu reassignment of the CGG codon in particular species. Analysis of tRNAs specified by sphaeroplealean mitogenomes provides direct support for and molecular underpinning of the proposed reassignments. Furthermore, we point to unique mutations in the mitochondrial release factor mtRF1a that correlate with changes in the use of termination codons in Sphaeropleales, including the two independent stop-to-sense UAG reassignments, the reintroduction of UGA in some Scenedesmaceae, and the sense-to-stop reassignment of UCA widespread in the group. Codon disappearance seems to be the main drive of the dynamic evolution of the mitochondrial genetic code in Sphaeropleales.

PhyloToL: A Taxon/Gene-Rich Phylogenomic Pipeline to Explore Genome Evolution of Diverse Eukaryotes

Estimating multiple sequence alignments (MSAs) and inferring phylogenies are essential for many aspects of comparative biology. Yet, many bioinformatics tools for such analyses have focused on specific clades, with greatest attention paid to plants, animals, and fungi. The rapid increase in high-throughput sequencing (HTS) data from diverse lineages now provides opportunities to estimate evolutionary relationships and gene family evolution across the eukaryotic tree of life. At the same time, these types of data are known to be error-prone (e.g., substitutions, contamination). To address these opportunities and challenges, we have refined a phylogenomic pipeline, now named PhyloToL, to allow easy incorporation of data from HTS studies, to automate production of both MSAs and gene trees, and to identify and remove contaminants. PhyloToL is designed for phylogenomic analyses of diverse lineages across the tree of life (i.e., at scales of >100 My). We demonstrate the power of PhyloToL by assessing stop codon usage in Ciliophora, identifying contamination in a taxon- and gene-rich database and exploring the evolutionary history of chromosomes in the kinetoplastid parasite Trypanosoma brucei, the causative agent of African sleeping sickness. Benchmarking PhyloToL’s homology assessment against that of OrthoMCL and a published paper on superfamilies of bacterial and eukaryotic organellar outer membrane pore-forming proteins demonstrates the power of our approach for determining gene family membership and inferring gene trees. PhyloToL is highly flexible and allows users to easily explore HTS data, test hypotheses about phylogeny and gene family evolution and combine outputs with third-party tools (e.g., PhyloChromoMap, iGTP).

UAA and UAG may Encode Amino Acid in Cathepsin B Gene of Euplotes octocarinatus

The ciliate Euplotes deviates from the universal genetic code by translating UGA as cysteine and using UAA and UAG as the termination codon. Here, we cloned and sequenced the Cathepsin B gene of Euplotes octocarinatus (Eo-CTSB) which containing several in-frame stop codons throughout the coding sequence. We provide evidences, based on 3'-RACE method and Western blot, that the Eo-CTSB gene is actively expressed. Comparison of the derived amino acid sequence with the homologs in other eukaryotes revealed that UAA and UAG may code for glutamine in Eo-CTSB. These findings imply an evolutionary complexity of stop codon reassignment in eukaryotes.

Phylogenetic studies on three Helicotylenchus species based on 28S rDNA and mtCOI sequence data

To facilitate the process of spiral nematode species delineation, populations of Helicotylenchus canadensis, H. pseudorobustus, and H. varicaudatus deriving from various locations and diverse natural and anthropogenic environments from Poland were investigated and characterized. For the first time, 28S rDNA sequences are reported for H. canadensis and H. varicaudatus, whereas new mtCOI sequences were acquired for all three analyzed species. A Bayesian phylogenetic analysis of the 28S rDNA fragments revealed that H. canadensis and H. varicaudatus are members of a clade that is a sister group to all other Helicotylenchus species; however, the closest known sister group to H. canadensis is H. vulgaris type A. Both 28S rDNA- and mtCOI-based phylogenetic results suggest that this clade excludes H. pseudorobustus, whose most recent common ancestor with the former species was inferred to be the ancestor of all Helicotylenchus species. Moreover, within the mtCOI sequences obtained from H. pseudorobustus, unlike from the other two, a simultaneous presence of TAG and TAA codons was identified. This may indicate mitochondrial genetic code alterations or other genomic rearrangements in H. pseudorobustus.

A precedented nuclear genetic code with all three termination codons reassigned as sense codons in the syndinean Amoebophrya sp. ex Karlodinium veneficum

Amoebophrya is part of an enigmatic, diverse, and ubiquitous marine alveolate lineage known almost entirely from anonymous environmental sequencing. Two cultured Amoebophrya strains grown on core dinoflagellate hosts were used for transcriptome sequencing. BLASTx using different genetic codes suggests that Amoebophyra sp. ex Karlodinium veneficum uses the three typical stop codons (UAA, UAG, and UGA) to encode amino acids. When UAA and UAG are translated as glutamine about half of the alignments have better BLASTx scores, and when UGA is translated as tryptophan one fifth have better scores. However, the sole stop codon appears to be UGA based on conserved genes, suggesting contingent translation of UGA. Neither host sequences, nor sequences from the second strain, Amoebophrya sp. ex Akashiwo sanguinea had similar results in BLASTx searches. A genome survey of Amoebophyra sp. ex K. veneficum showed no evidence for transcript editing aside from mitochondrial transcripts. The dynein heavy chain (DHC) gene family was surveyed and of 14 transcripts only two did not use UAA, UAG, or UGA in a coding context. Overall the transcriptome displayed strong bias for A or U in third codon positions, while the tRNA genome survey showed bias against codons ending in U, particularly for amino acids with two codons ending in either C or U. Together these clues suggest contingent translation mechanisms in Amoebophyra sp. ex K. veneficum and a phylogenetically distinct instance of genetic code modification.

Was the Mitochondrion Necessary to Start Eukaryogenesis?

Arguments based on cell energetics favour the view that a mitochondrion capable of oxidative phosphorylation was a prerequisite for the evolution of other features of the eukaryotic cell, including increased volume, genome size and, eventually, phagotrophy. Contrary to this we argue that: (i) extant amitochondriate eukaryotes possess voluminous phagotrophic cells with large genomes; (ii) picoeukaryotes demonstrate that phagotrophy is feasible at prokaryotic cell sizes; and (iii) the assumption that evolution of complex features requires extra ATP, often mentioned in this context, is unfounded and should not be used in such considerations. We claim that the diversity of cell organisations and functions observed today in eukaryotes gives no reason to postulate that a mitochondrion must have preceded phagocytosis in eukaryogenesis.

'Stop' in protein synthesis is modulated with exquisite subtlety by an extended RNA translation signal

Translational stop codons, UAA, UAG, and UGA, form an integral part of the universal genetic code. They are of significant interest today for their underlying fundamental role in terminating protein synthesis, but also for their potential utilisation for programmed alternative translation events. In diverse organisms, UAA has wide usage, but it is puzzling that the high fidelity UAG is selected against and yet UGA, vulnerable to suppression, is widely used, particularly in those archaeal and bacterial genomes with a high GC content. In canonical protein synthesis, stop codons are interpreted by protein release factors that structurally and functionally mimic decoding tRNAs and occupy the decoding site on the ribosome. The release factors make close contact with the decoding complex through multiple interactions. Correct interactions cause conformational changes resulting in new and enhanced contacts with the ribosome, particularly between specific bases in the mRNA and rRNA. The base following the stop codon (fourth or +4 base) may strongly influence decoding efficiency, facilitating alternative non-canonical events like frameshifting or selenocysteine incorporation. The fourth base is drawn into the decoding site with a compacted stop codon in the eukaryotic termination complex. Surprisingly, mRNA sequences upstream and downstream of this core tetranucleotide signal have a significant influence on the strength of the signal. Since nine bases downstream of the stop codon are within the mRNA channel, their interactions with rRNA, and r-proteins may affect efficiency. With this understanding, it is now possible to design stop signals of desired strength for specific applied purposes.

Translation Termination and Ribosome Recycling in Eukaryotes

Termination of mRNA translation occurs when a stop codon enters the A site of the ribosome, and in eukaryotes is mediated by release factors eRF1 and eRF3, which form a ternary eRF1/eRF3-guanosine triphosphate (GTP) complex. eRF1 recognizes the stop codon, and after hydrolysis of GTP by eRF3, mediates release of the nascent peptide. The post-termination complex is then disassembled, enabling its constituents to participate in further rounds of translation. Ribosome recycling involves splitting of the 80S ribosome by the ATP-binding cassette protein ABCE1 to release the 60S subunit. Subsequent dissociation of deacylated transfer RNA (tRNA) and messenger RNA (mRNA) from the 40S subunit may be mediated by initiation factors (priming the 40S subunit for initiation), by ligatin (eIF2D) or by density-regulated protein (DENR) and multiple copies in T-cell lymphoma-1 (MCT1). These events may be subverted by suppression of termination (yielding carboxy-terminally extended read-through polypeptides) or by interruption of recycling, leading to reinitiation of translation near the stop codon.

Rewriting the Genetic Code

The genetic code-the language used by cells to translate their genomes into proteins that perform many cellular functions-is highly conserved throughout natural life. Rewriting the genetic code could lead to new biological functions such as expanding protein chemistries with noncanonical amino acids (ncAAs) and genetically isolating synthetic organisms from natural organisms and viruses. It has long been possible to transiently produce proteins bearing ncAAs, but stabilizing an expanded genetic code for sustained function in vivo requires an integrated approach: creating recoded genomes and introducing new translation machinery that function together without compromising viability or clashing with endogenous pathways. In this review, we discuss design considerations and technologies for expanding the genetic code. The knowledge obtained by rewriting the genetic code will deepen our understanding of how genomes are designed and how the canonical genetic code evolved.