UGA is translated as cysteine in pheromone 3 of Euplotes octocarinatus

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.

mRNA context and translation factors determine decoding in alternative nuclear genetic codes

The genetic code is a set of instructions that determine how the information in our genetic material is translated into amino acids. In general, it is universal for all organisms, from viruses and bacteria to humans. However, in the last few decades, exceptions to this rule have been identified both in pro- and eukaryotes. In this review, we discuss the 16 described alternative eukaryotic nuclear genetic codes and observe theories of their appearance in evolution. We consider possible molecular mechanisms that allow codon reassignment. Most reassignments in nuclear genetic codes are observed for stop codons. Moreover, in several organisms, stop codons can simultaneously encode amino acids and serve as termination signals. In this case, the meaning of the codon is determined by the additional factors besides the triplets. A comprehensive review of various non-standard coding events in the nuclear genomes provides a new insight into the translation mechanism in eukaryotes.

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.

Comparative genomics reveals insight into the evolutionary origin of massively scrambled genomes

Ciliates are microbial eukaryotes that undergo extensive programmed genome rearrangement, a natural genome editing process that converts long germline chromosomes into smaller gene-rich somatic chromosomes. Three well-studied ciliates include Oxytricha trifallax, Tetrahymena thermophila, and Paramecium tetraurelia, but only the Oxytricha lineage has a massively scrambled genome, whose assembly during development requires hundreds of thousands of precisely programmed DNA joining events, representing the most complex genome dynamics of any known organism. Here we study the emergence of such complex genomes by examining the origin and evolution of discontinuous and scrambled genes in the Oxytricha lineage. This study compares six genomes from three species, the germline and somatic genomes for Euplotes woodruffi, Tetmemena sp., and the model ciliate O. trifallax. We sequenced, assembled, and annotated the germline and somatic genomes of E. woodruffi, which provides an outgroup, and the germline genome of Tetmemena sp. We find that the germline genome of Tetmemena is as massively scrambled and interrupted as Oxytricha's: 13.6% of its gene loci require programmed translocations and/or inversions, with some genes requiring hundreds of precise gene editing events during development. This study revealed that the earlier diverged spirotrich, E. woodruffi, also has a scrambled genome, but only roughly half as many loci (7.3%) are scrambled. Furthermore, its scrambled genes are less complex, together supporting the position of Euplotes as a possible evolutionary intermediate in this lineage, in the process of accumulating complex evolutionary genome rearrangements, all of which require extensive repair to assemble functional coding regions. Comparative analysis also reveals that scrambled loci are often associated with local duplications, supporting a gradual model for the origin of complex, scrambled genomes via many small events of DNA duplication and decay.

A computational screen for alternative genetic codes in over 250,000 genomes

The genetic code has been proposed to be a ‘frozen accident,’ but the discovery of alternative genetic codes over the past four decades has shown that it can evolve to some degree. Since most examples were found anecdotally, it is difficult to draw general conclusions about the evolutionary trajectories of codon reassignment and why some codons are affected more frequently. To fill in the diversity of genetic codes, we developed Codetta, a computational method to predict the amino acid decoding of each codon from nucleotide sequence data. We surveyed the genetic code usage of over 250,000 bacterial and archaeal genome sequences in GenBank and discovered five new reassignments of arginine codons (AGG, CGA, and CGG), representing the first sense codon changes in bacteria. In a clade of uncultivated Bacilli, the reassignment of AGG to become the dominant methionine codon likely evolved by a change in the amino acid charging of an arginine tRNA. The reassignments of CGA and/or CGG were found in genomes with low GC content, an evolutionary force that likely helped drive these codons to low frequency and enable their reassignment.

All life forms rely on a ‘code’ to translate their genetic information into proteins. This code relies on limited permutations of three nucleotides – the building blocks that form DNA and other types of genetic information. Each ‘triplet’ of nucleotides – or codon – encodes a specific amino acid, the basic component of proteins. Reading the sequence of codons in the right order will let the cell know which amino acid to assemble next on a growing protein. For instance, the codon CGG – formed of the nucleotides guanine (G) and cytosine (C) – codes for the amino acid arginine. From bacteria to humans, most life forms rely on the same genetic code. Yet certain organisms have evolved to use slightly different codes, where one or several codons have an altered meaning.

To better understand how alternative genetic codes have evolved, Shulgina and Eddy set out to find more organisms featuring these altered codons, creating a new software called Codetta that can analyze the genome of a microorganism and predict the genetic code it uses. Codetta was then used to sift through the genetic information of 250,000 microorganisms. This was made possible by the sequencing, in recent years, of the genomes of hundreds of thousands of bacteria and other microorganisms – including many never studied before.

These analyses revealed five groups of bacteria with alternative genetic codes, all of which had changes in the codons that code for arginine. Amongst these, four had genomes with a low proportion of guanine and cytosine nucleotides. This may have made some guanine and cytosine-rich arginine codons very rare in these organisms and, therefore, easier to be reassigned to encode another amino acid.

The work by Shulgina and Eddy demonstrates that Codetta is a new, useful tool that scientists can use to understand how genetic codes evolve. In addition, it can also help to ensure the accuracy of widely used protein databases, which assume which genetic code organisms use to predict protein sequences from their genomes.

Biodiversity-based development and evolution: the emerging research systems in model and non-model organisms

Evolutionary developmental biology, or Evo-Devo for short, has become an established field that, broadly speaking, seeks to understand how changes in development drive major transitions and innovation in organismal evolution. It does so via integrating the principles and methods of many subdisciplines of biology. Although we have gained unprecedented knowledge from the studies on model organisms in the past decades, many fundamental and crucially essential processes remain a mystery. Considering the tremendous biodiversity of our planet, the current model organisms seem insufficient for us to understand the evolutionary and physiological processes of life and its adaptation to exterior environments. The currently increasing genomic data and the recently available gene-editing tools make it possible to extend our studies to non-model organisms. In this review, we review the recent work on the regulatory signaling of developmental and regeneration processes, environmental adaptation, and evolutionary mechanisms using both the existing model animals such as zebrafish and Drosophila, and the emerging nonstandard model organisms including amphioxus, ascidian, ciliates, single-celled phytoplankton, and marine nematode. In addition, the challenging questions and new directions in these systems are outlined as well.

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.

Eukaryotic peptide chain release factor 1 participates in translation termination of specific cysteine-poor prolamines in rice endosperm

Prolamines are alcohol-soluble proteins classified as either cysteine-poor (CysP) or cysteine-rich (CysR) based on whether they can be alcohol-extracted without or with reducing agents, respectively. In rice esp1 mutants, various CysP prolamines exhibit both reduced and normal amounts of isoelectric focusing bands, indicating that the mutation affects only certain prolamine classes. To examine the genetic regulation of CysP prolamine synthesis and accumulation, we constructed a high-resolution genetic linkage map of ESP1. The ESP1 gene was mapped to within a 20 kb region on rice chromosome 7. Sequencing analysis of annotated genes in this region revealed a single-nucleotide polymorphism within eukaryotic peptide chain release factor (eRF1), which participates in stop-codon recognition and nascent-polypeptide release from ribosomes during translation. A subsequent complementation test revealed that ESP1 encodes eRF1. We also identified UAA as the stop codon of CysP prolamines with reduced concentration in esp1 mutants. Recognition assays and microarray analysis confirmed that ESP1/eRF1 recognizes UAA/UAG, but not UGA. Our results provide convincing evidence that ESP1/eRF1 participates in the translation termination of CysP prolamines during seed development.

EOGD: the Euplotes octocarinatus genome database

Background

Euplotes, a ciliated protozoan, is a useful unicellular model organism. Studies on Euplotes have provided excellent insights into various basic biological principles. We have recently sequenced the macronuclear genome of the common freshwater species Euplotes octocarinatus to provide novel insights into Euplotes genetics and molecular biology.

Results

In this study, we present the E. octocarinatus Genome Database (EOGD), a functional annotation and analysis platform for the global study of the Euplotes genome. EOGD includes macronuclear genomic and transcriptomic data, predicted gene models, coding sequences, protein sequences, and functional annotations. The GBrowser and BLAST tools are embedded in EOGD to enable the search, visualization and analysis of E. octocarinatus genomic and transcriptomic data.

Conclusions

EOGD is a useful resource for the research community, particularly for researchers who conduct genome-scale analysis and molecular biology studies of Euplotes or other ciliates. EOGD will be continuously updated to integrate more datasets and analytical tools. EOGD is freely available at http://ciliates.ihb.ac.cn/database/home/#eo.

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.

Expression of a MORN repeat protein from Euplotes octocarinatus requires a +1 programmed ribosomal frameshifting

Analysis of transcriptome revealed that a membrane occupation and recognition nexus (MORN) repeat protein-encoding gene of Euplotes octocarinatus (Eo-morn-9-31) was a candidate for programmed +1 ribosomal frameshifting (+1 PRF). In this study, a dual-luciferase assay was performed to detect its expression. The result showed that the MORN repeat protein (Eo-MORN-9-31) could be produced by the +1 PRF event during the process of translation in yeast and the frameshifting efficiency was about 4-5%. We further confirmed its reality by western blot and mass spectrometry. This study provided experimental evidence indicating that the expression of the Eo-MORN-9-31 of E. octocarinatus required the +1 PRF.

The Anti-Oxidant Defense System of the Marine Polar Ciliate Euplotes nobilii: Characterization of the MsrB Gene Family

Organisms living in polar waters must cope with an extremely stressful environment dominated by freezing temperatures, high oxygen concentrations and UV radiation. To shed light on the genetic mechanisms on which the polar marine ciliate, Euplotes nobilii, relies to effectively cope with the oxidative stress, attention was focused on methionine sulfoxide reductases which repair proteins with oxidized methionines. A family of four structurally distinct MsrB genes, encoding enzymes specific for the reduction of the methionine-sulfoxide R-forms, were identified from a draft of the E. nobilii transcriptionally active (macronuclear) genome. The En-MsrB genes are constitutively expressed to synthesize proteins markedly different in amino acid sequence, number of CXXC motifs for zinc-ion binding, and presence/absence of a cysteine residue specific for the mechanism of enzyme regeneration. The En-MsrB proteins take different localizations in the nucleus, mitochondria, cytosol and endoplasmic reticulum, ensuring a pervasive protection of all the major subcellular compartments from the oxidative damage. These observations have suggested to regard the En-MsrB gene activity as playing a central role in the genetic mechanism that enables E. nobilii and ciliates in general to live in the polar environment.

Large-scale mass spectrometry-based analysis of Euplotes octocarinatus supports the high frequency of +1 programmed ribosomal frameshift

Programmed ribosomal frameshifting (PRF) is commonly used to express many viral and some cellular genes. We conducted a genome-wide investigation of +1 PRF in ciliate Euplotes octocarinatus through genome and transcriptome sequencing and our results demonstrated that approximately 11.4% of genes require +1 PRF to produce complete gene products. While nucleic acid-based evidence for candidate genes with +1 PRF is strong, only very limited information is available at protein levels to date. In this study, E. octocarinatus was subjected to large-scale mass spectrometry-based analysis to verify the high frequency of +1 PRF and 226 +1 PRF gene products were identified. Based on the amino acid sequences of the peptides spanning the frameshift sites, typical frameshift motif AAA-UAR for +1 PRF in Euplotes was identified. Our data in this study provide very useful insight into the understanding of the molecular mechanism of +1 PRF.

Novel Ciliate Genetic Code Variants Including the Reassignment of All Three Stop Codons to Sense Codons in Condylostoma magnum

mRNA translation in many ciliates utilizes variant genetic codes where stop codons are reassigned to specify amino acids. To characterize the repertoire of ciliate genetic codes, we analyzed ciliate transcriptomes from marine environments. Using codon substitution frequencies in ciliate protein-coding genes and their orthologs, we inferred the genetic codes of 24 ciliate species. Nine did not match genetic code tables currently assigned by NCBI. Surprisingly, we identified a novel genetic code where all three standard stop codons (TAA, TAG, and TGA) specify amino acids in Condylostoma magnum. We provide evidence suggesting that the functions of these codons in C. magnum depend on their location within mRNA. They are decoded as amino acids at internal positions, but specify translation termination when in close proximity to an mRNA 3′ end. The frequency of stop codons in protein coding sequences of closely related Climacostomum virens suggests that it may represent a transitory state.

Revisiting fifty years of research on pheromone signaling in ciliates

Among protists, pheromones have been identified in a great variety of algal species for their activity in driving gamete-gamete interactions for fertilization. Analogously in ciliates, pheromones have been identified for their activity in inducing the sexual phenomenon of conjugation. Although this identification was pioneered by Kimball more than fifty years ago, an effective isolation and chemical characterization of ciliate pheromones has remained confined to species of Blepharisma, Dileptus and Euplotes. In Euplotes species, in which the molecular structures have been determined, pheromones form species-specific families of structurally homologous helical, cysteine-rich, highly-stable proteins. Being structurally homologous, they can bind cells in competition with one another, raising interesting functional analogies with the families of growth factors and cytokines that regulate cell differentiation and development in higher organisms. In addition to inducing conjugation by binding cells in heterologous fashion, Euplotes pheromones act also as autocrine growth factors by binding to, and promoting the vegetative reproduction of the same cells from which they originate. This autocrine activity is most likely primary, providing a concrete example of how the original function of a molecule can be obscured during evolution by the acquisition of a new one.

The Fungus Candida albicans Tolerates Ambiguity at Multiple Codons

The ascomycete Candida albicans is a normal resident of the gastrointestinal tract of humans and other warm-blooded animals. It occurs in a broad range of body sites and has high capacity to survive and proliferate in adverse environments with drastic changes in oxygen, carbon dioxide, pH, osmolarity, nutrients, and temperature. Its biology is unique due to flexible reassignment of the leucine CUG codon to serine and synthesis of statistical proteins. Under standard growth conditions, CUG sites incorporate leucine (3% of the times) and serine (97% of the times) on a proteome wide scale, but leucine incorporation fluctuates in response to environmental stressors and can be artificially increased up to 98%. In order to determine whether such flexibility also exists at other codons, we have constructed several serine tRNAs that decode various non-cognate codons. Expression of these tRNAs had minor effects on fitness, but growth of the mistranslating strains at different temperatures, in medium with different pH and nutrients composition was often enhanced relatively to the wild type (WT) strain, supporting our previous data on adaptive roles of CUG ambiguity in variable growth conditions. Parallel evolution of the recombinant strains (100 generations) followed by full genome resequencing identified various strain specific single nucleotide polymorphisms (SNP) and one SNP in the deneddylase (JAB1) gene in all strains. Since JAB1 is a subunit of the COP9 signalosome complex, which interacts with cullin (Cdc53p) to mediate degradation of a variety of cellular proteins, our data suggest that neddylation plays a key role in tolerance and adaptation to codon ambiguity in C. albicans.

High frequency of +1 programmed ribosomal frameshifting in Euplotes octocarinatus

Programmed -1 ribosomal frameshifting (-1 PRF) has been identified as a mechanism to regulate the expression of many viral genes and some cellular genes. The slippery site of -1 PRF has been well characterized, whereas the +1 PRF signal and the mechanism involved in +1 PRF remain poorly understood. Previous study confirmed that +1 PRF is required for the synthesis of protein products in several genes of ciliates from the genus Euplotes. To accurately assess the frequency of genes requiring frameshift in Euplotes, the macronuclear genome and transcriptome of Euplotes octocarinatus were analyzed in this study. A total of 3,700 +1 PRF candidate genes were identified from 32,353 transcripts, and the gene products of these putative +1 PRFs were mainly identified as protein kinases. Furthermore, we reported a putative suppressor tRNA of UAA which may provide new insights into the mechanism of +1 PRF in euplotids. For the first time, our transcriptome-wide survey of +1 PRF in E. octocarinatus provided a dataset which serves as a valuable resource for the future understanding of the mechanism underlying +1 PRF.

Non-Standard Genetic Codes Define New Concepts for Protein Engineering

The essential feature of the genetic code is the strict one-to-one correspondence between codons and amino acids. The canonical code consists of three stop codons and 61 sense codons that encode 20% of the amino acid repertoire observed in nature. It was originally designated as immutable and universal due to its conservation in most organisms, but sequencing of genes from the human mitochondrial genomes revealed deviations in codon assignments. Since then, alternative codes have been reported in both nuclear and mitochondrial genomes and genetic code engineering has become an important research field. Here, we review the most recent concepts arising from the study of natural non-standard genetic codes with special emphasis on codon re-assignment strategies that are relevant to engineering genetic code in the laboratory. Recent tools for synthetic biology and current attempts to engineer new codes for incorporation of non-standard amino acids are also reviewed in this article.

Stop codon recognition in the early-diverged protozoans Giardia lamblia and Trichomonas vaginalis

Two classes of polypeptide release factors (RFs) are responsible for maintaining accuracy in translation termination; however, their detailed mechanism of action and evolutionary history of these factors remain elusive. The structure and function of RFs vary in bacteria and eukaryotes, a fact that is suggestive of evolutionary changes in the translation termination system. Giardia lamblia (Diplomonada) and Trichomonas vaginalis (Parabasalia) are considered as early-diverged eukaryotes. The class II release factor, eRF3, of Giardia (Gl-eRF3) appears to have only one domain that corresponds to EF-1α and lacks the N-terminal domain, similar to that of eRF3 of other organisms. In the present study, we show that the chimeric molecules Gl/Sc eRF1 and Tv/Sc eRF1, which are composed of the N-terminal domain of Gl-eRF1 or Tv-eRF1, fused to the core domain (M and C domain) of Saccharomyces cerevisiae eRF1 (Sc-eRF1), resulting in loss of the RF properties of the N-terminal domain. This suggests that the conformation of eRF1 for stop codon recognition in Giardia and Trichomonas varies from the eRF1s of other eukaryotes, including ciliates and yeast. Further studies using intra-N-terminal chimeras of eRF1 indicated that the combination of the GTS loop and NIKS motif from Gl-eRF1 and the Y-C-F motif from Sc-eRF1within the N terminal domain of hybrid eRF1 could restore UGA, but not UAG and UGA recognition. In contrast, the combination of the GTS loop and the NIKS motif of Sc-eRF1 and the Y-C-F motif of Gl-eRF1 could restore UAG and UAA recognition, but not UGA recognition. Thus, these results confirm the findings of previous studies that three motifs in eRF1 are necessary for discrimination of the three bases of stop codons. The NIKS motif is responsible for recognition of the first two bases of UAA and UAG, and the Y-C-F motif identifies the second base of UGA by Gl-eRF1. Amino acid residue substitutions in Gl/Sc-eRF1 by corresponding residues of Sc-eRF1 could change and even restore RF activity, further suggesting different conformation of eRF1 are used for stop codon recognition in Giardia and in Saccharomyces.

Genetic code evolution reveals the neutral emergence of mutational robustness, and information as an evolutionary constraint

The standard genetic code (SGC) is central to molecular biology and its origin and evolution is a fundamental problem in evolutionary biology, the elucidation of which promises to reveal much about the origins of life. In addition, we propose that study of its origin can also reveal some fundamental and generalizable insights into mechanisms of molecular evolution, utilizing concepts from complexity theory. The first is that beneficial traits may arise by non-adaptive processes, via a process of "neutral emergence". The structure of the SGC is optimized for the property of error minimization, which reduces the deleterious impact of point mutations. Via simulation, it can be shown that genetic codes with error minimization superior to the SGC can emerge in a neutral fashion simply by a process of genetic code expansion via tRNA and aminoacyl-tRNA synthetase duplication, whereby similar amino acids are added to codons related to that of the parent amino acid. This process of neutral emergence has implications beyond that of the genetic code, as it suggests that not all beneficial traits have arisen by the direct action of natural selection; we term these "pseudaptations", and discuss a range of potential examples. Secondly, consideration of genetic code deviations (codon reassignments) reveals that these are mostly associated with a reduction in proteome size. This code malleability implies the existence of a proteomic constraint on the genetic code, proportional to the size of the proteome (P), and that its reduction in size leads to an "unfreezing" of the codon - amino acid mapping that defines the genetic code, consistent with Crick's Frozen Accident theory. The concept of a proteomic constraint may be extended to propose a general informational constraint on genetic fidelity, which may be used to explain variously, differences in mutation rates in genomes with differing proteome sizes, differences in DNA repair capacity and genome GC content between organisms, a selective pressure in the evolution of sexual reproduction, and differences in translational fidelity. Lastly, the utility of the concept of an informational constraint to other diverse fields of research is explored.

Ser or Leu: structural snapshots of mistranslation in Candida albicans

Candida albicans is a polymorphic opportunistic fungal pathogen normally residing as commensal on mucosal surfaces, skin and gastrointestinal and genitourinary tracts. However, in immunocompromised patients C. albicans can cause superficial mucosal infections or life-threatening disseminated candidemia. A change in physiological conditions triggers a cascade of molecular events leading to morphogenetic alterations and increased resistance to damage induced by host defenses. The complex biology of this human pathogen is reflected in its morphological plasticity and reinforced by the ability to ambiguously translate the universal leucine CUG codon predominantly as serine, but also as leucine. Mistranslation affects more than half of C. albicans proteome and it is widespread across many biological processes. A previous analysis of CTG-codon containing gene products in C. albicans suggested that codon ambiguity subtly shapes protein function and might have a pivotal role in signaling cascades associated with morphological changes and pathogenesis. In this review we further explore this hypothesis by highlighting the role of ambiguous decoding in macromolecular recognition of key effector proteins associated with the regulation of signal transduction cascades and the cell cycle, which are critical processes for C. albicans morphogenic plasticity under a variety of environmental conditions.

Modulation of efficiency of translation termination in Saccharomyces cerevisiae

Nonsense suppression is a readthrough of premature termination codons. It typically occurs either due to the recognition of stop codons by tRNAs with mutant anticodons, or due to a decrease in the fidelity of translation termination. In the latter case, suppressors usually promote the readthrough of different types of nonsense codons and are thus called omnipotent nonsense suppressors. Omnipotent nonsense suppressors were identified in yeast Saccharomyces cerevisiae in 1960s, and most of subsequent studies were performed in this model organism. Initially, omnipotent suppressors were localized by genetic analysis to different protein- and RNA-encoding genes, mostly the components of translational machinery. Later, nonsense suppression was found to be caused not only by genomic mutations, but also by epigenetic elements, prions. Prions are self-perpetuating protein conformations usually manifested by infectious protein aggregates. Modulation of translational accuracy by prions reflects changes in the activity of their structural proteins involved in different aspects of protein synthesis. Overall, nonsense suppression can be seen as a "phenotypic mirror" of events affecting the accuracy of the translational machine. However, the range of proteins participating in the modulation of translation termination fidelity is not fully elucidated. Recently, the list has been expanded significantly by findings that revealed a number of weak genetic and epigenetic nonsense suppressors, the effect of which can be detected only in specific genetic backgrounds. This review summarizes the data on the nonsense suppressors decreasing the fidelity of translation termination in S. cerevisiae, and discusses the functional significance of the modulation of translational accuracy.

The acidic ribosomal protein P2 from Euplotes octocarinatus is phosphorylated at its N-terminal domain

The eukaryotic acid ribosomal P0, P1, and P2 proteins share a conserved flexible C-terminal tail that is rich in acidic residues, which are involved in the interaction with elongation factor 2 during protein synthesis. Our previous work suggested that the acidic ribosomal P proteins from Euplotes octocarinatus have a special C-terminal domain. To further understand this characteristic feature, both P2 and elongation factor 2 from E. octocarinatus were overexpressed, for the first time, in Escherichia coli in this study. GST pull-down assay indicated that P2 protein from E. octocarinatus (EoP2) interacted specifically with the N-terminal domain of elongation factor 2 from E. octocarinatus (EoEF-2) in vitro. The interacting part of EoP2 is in the C-terminal domains, consistent with the observation in other organisms. Phosphorylation of the recombinant EoP2 was performed in vitro using multiple methods such as (31)P-NMR spectroscopy, native PAGE, and Phos-tag(TM) SDS-PAGE. Results showed that ribosomal protein EoP2 was phosphorylated by casein kinase II at serine 21 located at the N terminus. This phosphorylation site identified in EoP2 is quite different from that of P2 from other organisms, in which the phosphorylation site is located in the conserved C-terminal region.

Reversion of a fungal genetic code alteration links proteome instability with genomic and phenotypic diversification

Many fungi restructured their proteomes through incorporation of serine (Ser) at thousands of protein sites coded by the leucine (Leu) CUG codon. How these fungi survived this potentially lethal genetic code alteration and its relevance for their biology are not understood. Interestingly, the human pathogen Candida albicans maintains variable Ser and Leu incorporation levels at CUG sites, suggesting that this atypical codon assignment flexibility provided an effective mechanism to alter the genetic code. To test this hypothesis, we have engineered C. albicans strains to misincorporate increasing levels of Leu at protein CUG sites. Tolerance to the misincorporations was very high, and one strain accommodated the complete reversion of CUG identity from Ser back to Leu. Increasing levels of Leu misincorporation decreased growth rate, but production of phenotypic diversity on a phenotypic array probing various metabolic networks, drug resistance, and host immune cell responses was impressive. Genome resequencing revealed an increasing number of genotype changes at polymorphic sites compared with the control strain, and 80% of Leu misincorporation resulted in complete loss of heterozygosity in a large region of chromosome V. The data unveil unanticipated links between gene translational fidelity, proteome instability and variability, genome diversification, and adaptive phenotypic diversity. They also explain the high heterozygosity of the C. albicans genome and open the door to produce microorganisms with genetic code alterations for basic and applied research.

Identification of amino acids responsible for stop codon recognition for polypeptide chain release factor

One factor involved in eukaryotic translation termination is class 1 release factor in eukaryotes (eRF1), which functions to decode stop codons. Variant code species, such as ciliates, frequently exhibit altered stop codon recognition. Studies revealed that some class-specific residues in the eRF1 N-terminal domain are responsible for stop codon reassignment in ciliates. Here, we investigated the effects on stop codon recognition of chimeric eRF1s containing the N-terminal domain of Euplotes octocarinatus and Blepharisma japonicum eRF1 fused to Saccharomyces cerevisiae M and C domains using dual luciferase read-through assays. Mutation of class-specific residues in different eRF1 classes was also studied to identify key residues and motifs involved in stop codon decoding. As expected, our results demonstrate that 3 pockets within the eRF1 N-terminal domain were involved in decoding stop codon nucleotides. However, allocation of residues to each pocket was revalued. Our data suggest that hydrophobic and class-specific surface residues participate in different functions: modulation of pocket conformation and interaction with stop codon nucleotides, respectively. Residues conserved across all eRF1s determine the relative orientation of the 3 pockets according to stop codon nucleotides. However, quantitative analysis of variant ciliate and yeast eRF1 point mutants did not reveal any correlation between evolutionary conservation of class-specific residues and termination-related functional specificity and was limited in elucidating a detailed mechanism for ciliate stop codon reassignment. Thus, based on isolation of suppressor tRNAs from Euplotes and Tetrahymena, we propose that stop codon reassignment in ciliates may be controlled by cooperation between eRF1 and suppressor tRNAs.

Methionine sulfoxide reduction in ciliates: characterization of the ready-to-use methionine sulfoxide-R-reductase genes in Euplotes

Genes encoding the enzyme methionine sulfoxide reductase type B, specific to the reduction of the oxidized methionine-R form, were characterized from the expressed (macronuclear) genome of two ecologically separate marine species of Euplotes, i.e. temperate water E. raikovi and polar water E. nobilii. Both species were found to contain a single msrB gene with a very simple structural organization encoding a protein of 127 (E. raikovi) or 126 (E. nobilii) amino acid residues that belongs to the group of zinc-containing enzymes. Both msrB genes are constitutively expressed, suggesting that the MsrB enzyme plays an essential role in repairing oxidative damages that appear to be primarily caused by physiological cell aging in E. raikovi and by interactions with an O(2) saturated environment in E. nobilii.

Identification of eRF1 residues that play critical and complementary roles in stop codon recognition

The initiation and elongation stages of translation are directed by codon-anticodon interactions. In contrast, a release factor protein mediates stop codon recognition prior to polypeptide chain release. Previous studies have identified specific regions of eukaryotic release factor one (eRF1) that are important for decoding each stop codon. The cavity model for eukaryotic stop codon recognition suggests that three binding pockets/cavities located on the surface of eRF1's domain one are key elements in stop codon recognition. Thus, the model predicts that amino acid changes in or near these cavities should influence termination in a stop codon-dependent manner. Previous studies have suggested that the TASNIKS and YCF motifs within eRF1 domain one play important roles in stop codon recognition. These motifs are highly conserved in standard code organisms that use UAA, UAG, and UGA as stop codons, but are more divergent in variant code organisms that have reassigned a subset of stop codons to sense codons. In the current study, we separately introduced TASNIKS and YCF motifs from six variant code organisms into eRF1 of Saccharomyces cerevisiae to determine their effect on stop codon recognition in vivo. We also examined the consequences of additional changes at residues located between the TASNIKS and YCF motifs. Overall, our results indicate that changes near cavities two and three frequently mediated significant effects on stop codon selectivity. In particular, changes in the YCF motif, rather than the TASNIKS motif, correlated most consistently with variant code stop codon selectivity.

Unassigned codons, nonsense suppression, and anticodon modifications in the evolution of the genetic code

The origin of the genetic code is a central open problem regarding the early evolution of life. Here, we consider two undeveloped but important aspects of possible scenarios for the evolutionary pathway of the translation machinery: the role of unassigned codons in early stages of the code and the incorporation of tRNA anticodon modifications. As the first codons started to encode amino acids, the translation machinery likely was faced with a large number of unassigned codons. Current molecular scenarios for the evolution of the code usually assume the very rapid assignment of all codons before all 20 amino acids became encoded. We show that the phenomenon of nonsense suppression as observed in current organisms allows for a scenario in which many unassigned codons persisted throughout most of the evolutionary development of the code. In addition, we demonstrate that incorporation of anticodon modifications at a late stage is feasible. The wobble rules allow a set of 20 tRNAs fully lacking anticodon modifications to encode all 20 canonical amino acids. These observations have implications for the biochemical plausibility of early stages in the evolution of the genetic code predating tRNA anticodon modifications and allow for effective translation by a relatively small and simple early tRNA set.

FACIL: Fast and Accurate Genetic Code Inference and Logo

Motivation: The intensification of DNA sequencing will increasingly unveil uncharacterized species with potential alternative genetic codes. A total of 0.65% of the DNA sequences currently in Genbank encode their proteins with a variant genetic code, and these exceptions occur in many unrelated taxa.

Results: We introduce FACIL (Fast and Accurate genetic Code Inference and Logo), a fast and reliable tool to evaluate nucleic acid sequences for their genetic code that detects alternative codes even in species distantly related to known organisms. To illustrate this, we apply FACIL to a set of mitochondrial genomic contigs of Globobulimina pseudospinescens. This foraminifer does not have any sequenced close relative in the databases, yet we infer its alternative genetic code with high confidence values. Results are intuitively visualized in a Genetic Code Logo.

Availability and implementation: FACIL is available as a web-based service at http://www.cmbi.ru.nl/FACIL/ and as a stand-alone program.

Contact:dutilh@cmbi.ru.nl.

Supplementary information:Supplementary data are available at Bioinformatics online.

Isolation and structural characterization of two water-borne pheromones from Euplotes crassus, a ciliate commonly known to carry membrane-bound pheromones

Ciliates comprise species synthesizing water-diffusible mating type factors or pheromones and species synthesizing insoluble, cell membrane-bound pheromones. Euplotes crassus has traditionally been placed in the latter group. In contrast with this notion, we found that E. crassus is a constitutive pheromone-secreting ciliate, like other Euplotes species. From cell-free filtrate preparations of the E. crassus strain L-2D, we isolated two distinct pheromones, designated as Ec-α and Ec-1, and determined their complete amino acid sequences by combined chemical and genetic approaches. The Ec-α pheromone sequence extends for 56 amino acid residues with six cysteines and shows a molecular mass of 6,183 Da, while the Ec-1 pheromone sequence extends for 45 amino acid residues with 10 cysteines and shows a molecular mass of 4,840 Da. Marked structural differences distinguish the full-length Ec-α and Ec-1 coding sequences, which have been cloned and characterized from the transcriptionally active macronuclear genome. They were taken as clear indication that the Ec-α and Ec-1 pheromones are specified by genes that are not allelic, but likely derived from a duplicated genetic locus of the transcriptionally silent micronuclear genome.

Human SepSecS or SLA/LP: selenocysteine formation and autoimmune hepatitis

Selenocysteine, the 21st genetically encoded amino acid, is the major form of the antioxidant trace element selenium in the human body. In eukaryotes and archaea its synthesis proceeds through a phosphorylated intermediate in a tRNA-dependent fashion. The final step of selenocysteine formation is catalyzed by O-phosphoseryl-tRNA:selenocysteinyl-tRNA synthase (SepSecS) that converts phosphoseryl-tRNA(Sec) to selenocysteinyl-tRNA(Sec). The human SepSecS protein is also known as soluble liver antigen/liver pancreas (SLA/LP), which represents one of the antigens of autoimmune hepatitis. Here we review the discovery of human SepSecS and the current understanding of the immunogenicity of SLA/LP in autoimmune hepatitis.

A single amino acid change of translation termination factor eRF1 switches between bipotent and omnipotent stop-codon specificity

In eukaryotes a single class-1 translation termination factor eRF1 decodes the three stop codons: UAA, UAG and UGA. Some ciliates, like Euplotes, have a variant code, and here eRF1s exhibit UAR-only specificity, whereas UGA is reassigned as a sense codon. Since eukaryote eRF1 stop-codon recognition is associated with its N-terminal domain, structural features should exist in the N domain of ciliate eRF1s that restrict their stop-codon specificity. Using an in vitro reconstituted eukaryotic translation system we demonstrate here that a chimeric eRF1 composed of the N domain of Euplotes aediculatus eRF1 fused to the MC domains of human eRF1 exhibits UAR-only specificity. Functional analysis of eRF1 chimeras constructed by swapping Euplotes N domain sequences with the cognate regions from human eRF1 as well as site-directed mutagenesis of human eRF1 highlighted the crucial role of the alanine residue in position 70 of E. aediculatus eRF1 in restricting UGA decoding. Switching the UAR-only specificity of E. aediculatus eRF1 to omnipotent mode is due to a single point mutation. Furthermore, we examined the influence of eRF3 on the ability of chimeric and mutant eRF1s to induce peptide release in response to different stop codons.

Isolation and characterization of a novel copper-inducible metallothionein gene of a ciliate, Tetrahymena tropicalis lahorensis

The two isoforms of copper metallothionein (CuMT) gene of a copper resistant ciliate, Tetrahymena tropicalis lahorensis (Ttl), have been isolated and characterized. The molecular cloning and nucleotide sequencing of cDNAs coding for the two CuMT isoforms revealed that TtlCuMT1 gene has 300, while TtlCuMT2 has 327 nucleotides, both with ATG as the initiation codon and TGA as the translational termination codon. TAG codes for glutamine in TtlCuMT2 gene which is peculiar to Tetrahymena. The deduced or translated TtlCuMT1 and TtlCuMT2 peptide sequences contain 100 and 108 amino acid residues including 28 and 32 cysteine residues, respectively. The amino acid sequences of TtlCuMT1 and TtlCuMT2 have special features of two and three CXCXXCXCXXCXC intragenic tandem repeats with a conserved structural pattern of cysteine, respectively. The predicted tertiary structures of these two isoforms indicate two domains. Domain I and the initial part of domain II showed >98% homology with other Tetrahymena CuMT. On the basis of the differences in the domain II, the metallothionein subfamily 7b can be divided into two groups, one (TtlCuMT1) comprising >100 amino acids and the other (TtlCuMT2) comprising <100 amino acids. This is a novel finding of the present study as no such report on this type of classification exists at the moment. TtlCuMT1 has 95%, while TtlCuMT2 has 97% resemblance with the previously reported CuMT genes of Tetrahymena spp. SDS-PAGE analysis using fluorescent probe as well as coomassie brilliant blue staining also confirmed the presence of metallothionein.

Deciphering synonymous codons in the three domains of life: co-evolution with specific tRNA modification enzymes

The strategies organisms use to decode synonymous codons in cytosolic protein synthesis are not uniform. The complete isoacceptor tRNA repertoire and the type of modified nucleoside found at the wobble position 34 of their anticodons were analyzed in all kingdoms of life. This led to the identification of four main decoding strategies that are diversely used in Bacteria, Archaea and Eukarya. Many of the modern tRNA modification enzymes acting at position 34 of tRNAs are present only in specific domains and obviously have arisen late during evolution. In an evolutionary fine-tuning process, these enzymes must have played an essential role in the progressive introduction of new amino acids, and in the refinement and standardization of the canonical nuclear genetic code observed in all extant organisms (functional convergent evolutionary hypothesis).

Characterization of the pheromone gene family of an Antarctic and Arctic protozoan ciliate, Euplotes nobilii

Allelic genes encoding water-borne signal proteins (pheromones) were amplified and sequenced from the somatic (macronuclear) sub-chromosomic genome of Antarctic and Arctic strains of the marine ciliate, Euplotes nobilii. Their open reading frames appeared to be specific for polypeptide sequences of 83 to 94 amino acids identifiable with cytoplasmic pheromone precursors (pre-pro-pheromones), requiring two proteolytic steps to remove the pre- and pro-segments and secrete the mature pheromones. Differently from most of the macronuclear genes that have so far been characterized from Euplotes and other hypotrich ciliates, the 5' and 3' non-coding regions of all the seven E. nobilii pheromone genes are much longer than the coding regions (621 to 700 versus 214 to 285 nucleotides), and the 5' regions in particular show nearly identical sequences across the whole set of pheromone genes. These structural peculiarities of the non-coding regions are likely due to the presence of intron sequences and provide presumptive evidence that they are site of basic, conserved activities in the mechanism that regulates the expression of the E. nobilii pheromone genes.

Connection between stop codon reassignment and frequent use of shifty stop frameshifting

Ciliated protozoa of the genus Euplotes have undergone genetic code reassignment, redefining the termination codon UGA to encode cysteine. In addition, Euplotes spp. genes very frequently employ shifty stop frameshifting. Both of these phenomena involve noncanonical events at a termination codon, suggesting they might have a common cause. We recently demonstrated that Euplotes octocarinatus peptide release factor eRF1 ignores UGA termination codons while continuing to recognize UAA and UAG. Here we show that both the Tetrahymena thermophila and E. octocarinatus eRF1 factors allow efficient frameshifting at all three termination codons, suggesting that UGA redefinition also impaired UAA/UAG recognition. Mutations of the Euplotes factor restoring a phylogenetically conserved motif in eRF1 (TASNIKS) reduced programmed frameshifting at all three termination codons. Mutation of another conserved residue, Cys124, strongly reduces frameshifting at UGA while actually increasing frameshifting at UAA/UAG. We will discuss these results in light of recent biochemical characterization of these mutations.

Penultimate selenocysteine residue replaced by cysteine in thioredoxin reductase from selenium-deficient rat liver

Selenium is an essential micronutrient for humans and animals, and its deficiency can predispose to the development of pathological conditions. This study evaluates the effect of selenium deficiency on the thioredoxin system, consisting of NADPH, selenoprotein thioredoxin reductase (TrxR), and thioredoxin (Trx); and the glutathione system, including NADPH, glutathione reductase, glutathione, and glutaredoxin coupled with selenoprotein glutathione peroxidase (GPx). We particularly investigate whether inactive truncated TrxR is present under selenium-starvation conditions due to reading of the UGA codon as stop. Feeding rats a selenium-deficient diet resulted in a large decrease in activity of TrxR and GPx in rat liver but not in the levels of Trx1 and Grx1. However, selenium deficiency induced mitochondrial Grx2 10-fold and markedly changed the expression of some flavoproteins that are involved in the cellular folate, glucose, and lipid metabolism. Liver TrxR mRNA was nearly unchanged, but no truncated enzyme was found. Instead, a low-activity form of TrxR with a cysteine substituted for the penultimate selenocysteine in the C-terminal active site was identified in selenium-deficient rat liver. These results show a novel mechanism for decoding the UGA stop codon, inserting cysteine to make a full-length enzyme that may be required for selenium assimilation.

Genetic code supports targeted insertion of two amino acids by one codon

Strict one-to-one correspondence between codons and amino acids is thought to be an essential feature of the genetic code. However, we report that one codon can code for two different amino acids with the choice of the inserted amino acid determined by a specific 3' untranslated region structure and location of the dual-function codon within the messenger RNA (mRNA). We found that the codon UGA specifies insertion of selenocysteine and cysteine in the ciliate Euplotes crassus, that the dual use of this codon can occur even within the same gene, and that the structural arrangements of Euplotes mRNA preserve location-dependent dual function of UGA when expressed in mammalian cells. Thus, the genetic code supports the use of one codon to code for multiple amino acids.

The water-borne protein signals (pheromones) of the Antarctic ciliated protozoan Euplotes nobilii: structure of the gene coding for the En-6 pheromoneThis article is one of a selection of papers in the Special Issue on Polar and Alpine Microbiology

The marine Antarctic ciliate, Euplotes nobilii , secretes a family of water-borne signal proteins, denoted as pheromones, which control vegetative proliferation and mating in the cell. Based on the knowledge of the amino acid sequences of a set of these pheromones isolated from the culture supernatant of wild-type strains, we designed probes to identify their encoding genes in the cell somatic nucleus (macronucleus). The full-length gene of the pheromone En-6 was determined and found to contain an open-reading frame specific for the synthesis of the En-6 cytoplasmic precursor (pre-pro-En-6), which requires 2 proteolytic cleavages to remove the signal peptide (pre) and the prosegment before secretion of the mature protein. In contrast to the sequence variability that distinguishes the secreted pheromones, the pre- and pro-sequences appear to be tightly conserved and useful for the construction of probes to clone every other E. nobilii pheromone gene. Potential intron sequences in the coding region of the En-6 gene imply the synthesis of more En-6 isoforms.

Different roles of two gamma-tubulin isotypes in the cytoskeleton of the Antarctic ciliate Euplotes focardii: remodelling of interaction surfaces may enhance microtubule nucleation at low temperature

Gamma-tubulin belongs to the tubulin superfamily and plays an essential role in the nucleation of cellular microtubules. In the present study, we report the characterization of gamma-tubulin from the psychrophilic Antarctic ciliate Euplotes focardii. In this organism, gamma-tubulin is encoded by two genes, gamma-T1 and gamma-T2, that produce distinct isotypes. Comparison of the gamma-T1 and gamma-T2 primary sequences to a Euplotesgamma-tubulin consensus, derived from mesophilic (i.e. temperate) congeneric species, revealed the presence of numerous unique amino acid substitutions, particularly in gamma-T2. Structural models of gamma-T1 and gamma-T2, obtained using the 3D structure of human gamma-tubulin as a template, suggest that these substitutions are responsible for conformational and/or polarity differences located: (a) in the regions involved in longitudinal 'plus end' contacts; (b) in the T3 loop that participates in binding GTP; and (c) in the M loop that forms lateral interactions. Relative to gamma-T1, the gamma-T2 gene is amplified by approximately 18-fold in the macronuclear genome and is very strongly transcribed. Using confocal immunofluorescence microscopy, we found that the gamma-tubulins of E. focardii associate throughout the cell cycle with basal bodies of the non-motile dorsal cilia and of all of the cirri of the ventral surface (i.e. adoral membranelles, paraoral membrane, and frontoventral transverse, caudal and marginal cirri). By contrast, only gamma-T2 interacts with the centrosomes of the spindle during micronuclear mitosis. We also established that the gamma-T1 isotype associates only with basal bodies. Our results suggest that gamma-T1 and gamma-T2 perform different functions in the organization of the microtubule cytoskeleton of this protist and are consistent with the hypothesis that gamma-T1 and gamma-T2 have evolved sequence-based structural alterations that facilitate template nucleation of microtubules by the gamma-tubulin ring complex at cold temperatures.