Genomics: Evolution of the Genetic Code

Evolution of the Genetic Code in the Ascoideales (CUG-Ser2) Yeast Clade: The Ancestral tRNA-Leu(CAG) Gene Is Retained in Most Saccharomycopsis Species but Is Nonessential and Not Used for Translation

In the yeast genera Saccharomycopsis and Ascoidea, which comprise the taxonomic order Ascoideales, nuclear genes use a nonstandard genetic code in which CUG codons are translated as serine instead of leucine, due to a tRNA-Ser with the unusual anticodon CAG. However, some species in this clade also retain an ancestral tRNA-Leu gene with the same anticodon. One of these species, Ascoidea asiatica, has been shown to have a stochastic proteome in which proteins contain ∼50% Ser and 50% Leu at CUG codon sites, whereas previously examined Saccharomycopsis species translate CUG only as Ser. Here, we investigated the presence, conservation, and possible functionality of the tRNA-Leu(CAG) gene in the genus Saccharomycopsis. We sequenced the genomes of 23 strains that, together with previously available data, include almost every known species of this genus. We found that most Saccharomycopsis species have genes for both tRNA-Leu(CAG) and tRNA-Ser(CAG). However, tRNA-Leu(CAG) has been lost in Saccharomycopsis synnaedendra and Saccharomycopsis microspora, and its predicted cloverleaf structure is aberrant in all the other Saccharomycopsis species. We deleted the tRNA-Leu(CAG) gene of Saccharomycopsis capsularis and found that it is not essential. Proteomic analyses in vegetative and sporulating cultures of S. capsularis and Saccharomycopsis fermentans showed only translation of CUG as Ser. Despite its unusual structure, the tRNA-Leu(CAG) gene shows evidence of sequence conservation among Saccharomycopsis species, particularly in its acceptor stem and leucine identity elements, which suggests that it may have been retained in order to carry out an unknown nontranslational function.

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.

Trypanosomes as a magnifying glass for cell and molecular biology

The African trypanosome, Trypanosoma brucei, has developed into a flexible and robust experimental model for molecular and cellular parasitology, allowing us to better combat these and related parasites that cause worldwide suffering. Diminishing case numbers, due to efficient public health efforts, and recent development of new drug treatments have reduced the need for continued study of T. brucei in a disease context. However, we argue that this pathogen has been instrumental in revolutionary discoveries that have widely informed molecular and cellular biology and justifies continuing research as an experimental model. Ongoing work continues to contribute towards greater understanding of both diversified and conserved biological features. We discuss multiple examples where trypanosomes pushed the boundaries of cell biology and hope to inspire researchers to continue exploring these remarkable protists as tools for magnifying the inner workings of cells.

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.

Single-Cell Genomics Reveals the Divergent Mitochondrial Genomes of Retaria (Foraminifera and Radiolaria)

Mitochondria originated from an ancient bacterial endosymbiont that underwent reductive evolution by gene loss and endosymbiont gene transfer to the nuclear genome. The diversity of mitochondrial genomes published to date has revealed that gene loss and transfer processes are ongoing in many lineages. Most well-studied eukaryotic lineages are represented in mitochondrial genome databases, except for the superphylum Retaria-the lineage comprising Foraminifera and Radiolaria. Using single-cell approaches, we determined two complete mitochondrial genomes of Foraminifera and two nearly complete mitochondrial genomes of radiolarians. We report the complete coding content of an additional 14 foram species. We show that foraminiferan and radiolarian mitochondrial genomes contain a nearly fully overlapping but reduced mitochondrial gene complement compared to other sequenced rhizarians. In contrast to animals and fungi, many protists encode a diverse set of proteins on their mitochondrial genomes, including several ribosomal genes; however, some aerobic eukaryotic lineages (euglenids, myzozoans, and chlamydomonas-like algae) have reduced mitochondrial gene content and lack all ribosomal genes. Similar to these reduced outliers, we show that retarian mitochondrial genomes lack ribosomal protein and tRNA genes, contain truncated and divergent small and large rRNA genes, and contain only 14 or 15 protein-coding genes, including nad1, -3, -4, -4L, -5, and -7, cob, cox1, -2, and -3, and atp1, -6, and -9, with forams and radiolarians additionally carrying nad2 and nad6, respectively. In radiolarian mitogenomes, a noncanonical genetic code was identified in which all three stop codons encode amino acids. Collectively, these results add to our understanding of mitochondrial genome evolution and fill in one of the last major gaps in mitochondrial sequence databases. IMPORTANCE We present the reduced mitochondrial genomes of Retaria, the rhizarian lineage comprising the phyla Foraminifera and Radiolaria. By applying single-cell genomic approaches, we found that foraminiferan and radiolarian mitochondrial genomes contain an overlapping but reduced mitochondrial gene complement compared to other sequenced rhizarians. An alternative genetic code was identified in radiolarian mitogenomes in which all three stop codons encode amino acids. Collectively, these results shed light on the divergent nature of the mitochondrial genomes from an ecologically important group, warranting further questions into the biological underpinnings of gene content variability and genetic code variation between mitochondrial genomes.

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.

Frameshift and wild-type proteins are often highly similar because the genetic code and genomes were optimized for frameshift tolerance

Frameshift mutations have been considered of significant importance for the molecular evolution of proteins and their coding genes, while frameshift protein sequences encoded in the alternative reading frames of coding genes have been considered to be meaningless. However, functional frameshifts have been found widely existing. It was puzzling how a frameshift protein kept its structure and functionality while substantial changes occurred in its primary amino-acid sequence. This study shows that the similarities among frameshifts and wild types are higher than random similarities and are determined at different levels. Frameshift substitutions are more conservative than random substitutions in the standard genetic code (SGC). The frameshift substitutions score of SGC ranks in the top 2.0-3.5% of alternative genetic codes, showing that SGC is nearly optimal for frameshift tolerance. In many genes and certain genomes, frameshift-resistant codons and codon pairs appear more frequently than expected, suggesting that frameshift tolerance is achieved through not only the optimality of the genetic code but, more importantly, the further optimization of a specific gene or genome through the usages of codons/codon pairs, which sheds light on the role of frameshift mutations in molecular and genomic evolution.

Engineering Translation Components for Genetic Code Expansion

The expansion of the genetic code consisting of four bases and 20 amino acids into diverse building blocks has been an exciting topic in synthetic biology. Many biochemical components are involved in gene expression; therefore, adding a new component to the genetic code requires engineering many other components that interact with it. Genetic code expansion has advanced significantly for the last two decades with the engineering of several components involved in protein synthesis. These components include tRNA/aminoacyl-tRNA synthetase, new codons, ribosomes, and elongation factor Tu. In addition, biosynthesis and enhanced uptake of non-canonical amino acids have been attempted and have made meaningful progress. This review discusses the efforts to engineer these translation components, to improve the genetic code expansion technology.

A mitochondrial cytidine deaminase is responsible for C to U editing of tRNATrp to decode the UGA codon in Trypanosoma brucei

In kinetoplastid protists, all mitochondrial tRNAs are encoded in the nucleus and imported from the cytoplasm to maintain organellar translation. This also applies to the tryptophanyl tRNA (tRNA^Trp) encoded by a single-copy nuclear gene, with a CCA anticodon to read UGG codon used in the cytosolic translation. Yet, in the mitochondrion it is unable to decode the UGA codon specifying tryptophan. Following mitochondrial import of tRNA^Trp, this problem is solved at the RNA level by a single C34 to U34 editing event that creates the UCA anticodon, recognizing UGA. To identify the enzyme responsible for this critical editing activity, we scrutinized the genome of Trypanosoma brucei for putative cytidine deaminases as the most likely candidates. Using RNAi silencing and poisoned primer extension, we have identified a novel deaminase enzyme, named here TbmCDAT for mitochondrial Cytidine Deaminase Acting on tRNA, which is responsible for this organelle-specific activity in T. brucei. The ablation of TbmCDAT led to the downregulation of mitochondrial protein synthesis, supporting its role in decoding the UGA tryptophan codon.

TRAPID 2.0: a web application for taxonomic and functional analysis of de novo transcriptomes

Advances in high-throughput sequencing have resulted in a massive increase of RNA-Seq transcriptome data. However, the promise of rapid gene expression profiling in a specific tissue, condition, unicellular organism or microbial community comes with new computational challenges. Owing to the limited availability of well-resolved reference genomes, de novo assembled (meta)transcriptomes have emerged as popular tools for investigating the gene repertoire of previously uncharacterized organisms. Yet, despite their potential, these datasets often contain fragmented or contaminant sequences, and their analysis remains difficult. To alleviate some of these challenges, we developed TRAPID 2.0, a web application for the fast and efficient processing of assembled transcriptome data. The initial processing phase performs a global characterization of the input data, providing each transcript with several layers of annotation, comprising structural, functional, and taxonomic information. The exploratory phase enables downstream analyses from the web application. Available analyses include the assessment of gene space completeness, the functional analysis and comparison of transcript subsets, and the study of transcripts in an evolutionary context. A comparison with similar tools highlights TRAPID’s unique features. Finally, analyses performed within TRAPID 2.0 are complemented by interactive data visualizations, facilitating the extraction of new biological insights, as demonstrated with diatom community metatranscriptomes.

The role of non-standard translation in Candida albicans pathogenesis

Candida albicans typically resides in the human gastrointestinal tract and mucosal membranes as a commensal organism. To adapt and cope with the host immune system, it has evolved a variety of mechanisms of adaptation such as stress-induced mutagenesis and epigenetic regulation. Niche-specific patterns of gene expression also allow the fungus to fine-tune its response to specific microenvironments in the host and switch from harmless commensal to invasive pathogen. Proteome plasticity produced by CUG ambiguity, on the other hand is emerging as a new layer of complexity in C. albicans adaptation, pathogenesis, and drug resistance. Such proteome plasticity is the result of a genetic code alteration where the leucine CUG codon is translated mainly as serine (97%), but maintains some level of leucine (3%) assignment. In this review, we dissect the link between C. albicans non-standard CUG translation, proteome plasticity, host adaptation and pathogenesis. We discuss published work showing how this pathogen uses the fidelity of protein synthesis to spawn novel virulence traits.

Genome recoding strategies to improve cellular properties: mechanisms and advances

The genetic code, once believed to be universal and immutable, is now known to contain many variations and is not quite universal. The basis for genome recoding strategy is genetic code variation that can be harnessed to improve cellular properties. Thus, genome recoding is a promising strategy for the enhancement of genome flexibility, allowing for novel functions that are not commonly documented in the organism in its natural environment. Here, the basic concept of genetic code and associated mechanisms for the generation of genetic codon variants, including biased codon usage, codon reassignment, and ambiguous decoding, are extensively discussed. Knowledge of the concept of natural genetic code expansion is also detailed. The generation of recoded organisms and associated mechanisms with basic targeting components, including aminoacyl-tRNA synthetase-tRNA pairs, elongation factor EF-Tu and ribosomes, are highlighted for a comprehensive understanding of this concept. The research associated with the generation of diverse recoded organisms is also discussed. The success of genome recoding in diverse multicellular organisms offers a platform for expanding protein chemistry at the biochemical level with non-canonical amino acids, genetically isolating the synthetic organisms from the natural ones, and fighting viruses, including SARS-CoV2, through the creation of attenuated viruses. In conclusion, genome recoding can offer diverse applications for improving cellular properties in the genome-recoded organisms.

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.

Unexpected mitochondrial genome diversity revealed by targeted single-cell genomics of heterotrophic flagellated protists

Most eukaryotic microbial diversity is uncultivated, under-studied and lacks nuclear genome data. Mitochondrial genome sampling is more comprehensive, but many phylogenetically important groups remain unsampled. Here, using a single-cell sorting approach combining tubulin-specific labelling with photopigment exclusion, we sorted flagellated heterotrophic unicellular eukaryotes from Pacific Ocean samples. We recovered 206 single amplified genomes, predominantly from underrepresented branches on the tree of life. Seventy single amplified genomes contained unique mitochondrial contigs, including 21 complete or near-complete mitochondrial genomes from formerly under-sampled phylogenetic branches, including telonemids, katablepharids, cercozoans and marine stramenopiles, effectively doubling the number of available samples of heterotrophic flagellate mitochondrial genomes. Collectively, these data identify a dynamic history of mitochondrial genome evolution including intron gain and loss, extensive patterns of genetic code variation and complex patterns of gene loss. Surprisingly, we found that stramenopile mitochondrial content is highly plastic, resembling patterns of variation previously observed only in plants.

Chimeric Translation for Mitochondrial Peptides: Regular and Expanded Codons

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

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.

Novel genetic code and record-setting AT-richness in the highly reduced plastid genome of the holoparasitic plant Balanophora

Plastid genomes (plastomes) vary enormously in size and gene content among the many lineages of nonphotosynthetic plants, but key lineages remain unexplored. We therefore investigated plastome sequence and expression in the holoparasitic and morphologically bizarre Balanophoraceae. The two Balanophora plastomes examined are remarkable, exhibiting features rarely if ever seen before in plastomes or in any other genomes. At 15.5 kb in size and with only 19 genes, they are among the most reduced plastomes known. They have no tRNA genes for protein synthesis, a trait found in only three other plastid lineages, and thus Balanophora plastids must import all tRNAs needed for translation. Balanophora plastomes are exceptionally compact, with numerous overlapping genes, highly reduced spacers, loss of all cis-spliced introns, and shrunken protein genes. With A+T contents of 87.8% and 88.4%, the Balanophora genomes are the most AT-rich genomes known save for a single mitochondrial genome that is merely bloated with AT-rich spacer DNA. Most plastid protein genes in Balanophora consist of ≥90% AT, with several between 95% and 98% AT, resulting in the most biased codon usage in any genome described to date. A potential consequence of its radical compositional evolution is the novel genetic code used by Balanophora plastids, in which TAG has been reassigned from stop to tryptophan. Despite its many exceptional properties, the Balanophora plastome must be functional because all examined genes are transcribed, its only intron is correctly trans-spliced, and its protein genes, although highly divergent, are evolving under various degrees of selective constraint.

Mitogenomics Reveals a Novel Genetic Code in Hemichordata

The diverse array of codon reassignments demonstrate that the genetic code is not universal in nature. Exploring mechanisms underlying codon reassignment is critical for understanding the evolution of the genetic code during translation. Hemichordata, comprising worm-like Enteropneusta and colonial filter-feeding Pterobranchia, is the sister taxon of echinoderms and is more distantly related to chordates. However, only a few hemichordate mitochondrial genomes have been sequenced, hindering our understanding of mitochondrial genome evolution within Deuterostomia. In this study, we sequenced four mitochondrial genomes and two transcriptomes, including representatives of both major hemichordate lineages and analyzed together with public available data. Contrary to the current understanding of the mitochondrial genetic code in hemichordates, our comparative analyses suggest that UAA encodes Tyr instead of a "Stop" codon in the pterobranch lineage Cephalodiscidae. We also predict that AAA encodes Lys in pterobranch and enteropneust mitochondrial genomes, contradicting the previous assumption that hemichordates share the same genetic code with echinoderms for which AAA encodes Asn. Thus, we propose a new mitochondrial genetic code for Cephalodiscus and a revised code for enteropneusts. Moreover, our phylogenetic analyses are largely consistent with previous phylogenomic studies. The only exception is the phylogenetic position of the enteropneust Stereobalanus, whose placement as sister to all other described enteropneusts. With broader taxonomic sampling, we provide evidence that evolution of mitochondrial gene order and genetic codes in Hemichordata are more dynamic than previously thought and these findings provide insights into mitochondrial genome evolution within this clade.

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.

CoreTracker: accurate codon reassignment prediction, applied to mitochondrial genomes

Motivation

Codon reassignments have been reported across all domains of life. With the increasing number of sequenced genomes, the development of systematic approaches for genetic code detection is essential for accurate downstream analyses. Three automated prediction tools exist so far: FACIL, GenDecoder and Bagheera; the last two respectively restricted to metazoan mitochondrial genomes and CUG reassignments in yeast nuclear genomes. These tools can only analyze a single genome at a time and are often not followed by a validation procedure, resulting in a high rate of false positives.

Results

We present CoreTracker, a new algorithm for the inference of sense-to-sense codon reassignments. CoreTracker identifies potential codon reassignments in a set of related genomes, then uses statistical evaluations and a random forest classifier to predict those that are the most likely to be correct. Predicted reassignments are then validated through a phylogeny-aware step that evaluates the impact of the new genetic code on the protein alignment. Handling simultaneously a set of genomes in a phylogenetic framework, allows tracing back the evolution of each reassignment, which provides information on its underlying mechanism. Applied to metazoan and yeast genomes, CoreTracker significantly outperforms existing methods on both precision and sensitivity.

Availability and implementation

CoreTracker is written in Python and available at https://github.com/UdeM-LBIT/CoreTracker.

Contact

mabrouk@iro.umontreal.ca.

Supplementary information

Supplementary data are available at Bioinformatics online.

Endogenous Stochastic Decoding of the CUG Codon by Competing Ser- and Leu-tRNAs in Ascoidea asiatica

Although the “universal” genetic code is now known not to be universal, and stop codons can have multiple meanings, one regularity remains, namely that for a given sense codon there is a unique translation. Examining CUG usage in yeasts that have transferred CUG away from leucine, we here report the first example of dual coding: Ascoidea asiatica stochastically encodes CUG as both serine and leucine in approximately equal proportions. This is deleterious, as evidenced by CUG codons being rare, never at conserved serine or leucine residues, and predominantly in lowly expressed genes. Related yeasts solve the problem by loss of function of one of the two tRNAs. This dual coding is consistent with the tRNA-loss-driven codon reassignment hypothesis, and provides a unique example of a proteome that cannot be deterministically predicted.

Video Abstract

•

Ascoidea asiatica stochastically encodes CUG as leucine and serine

•

It is the only known example of a proteome with non-deterministic features

•

Stochastic encoding is caused by competing tRNALeu(CAG) and tRNASer(CAG)

•

A. asiatica copes with stochastic encoding by avoiding CUG at key positions

Recent advances in trypanosomatid research: genome organization, expression, metabolism, taxonomy and evolution

Unicellular flagellates of the family Trypanosomatidae are obligatory parasites of invertebrates, vertebrates and plants. Dixenous species are aetiological agents of a number of diseases in humans, domestic animals and plants. Their monoxenous relatives are restricted to insects. Because of the high biological diversity, adaptability to dramatically different environmental conditions, and omnipresence, these protists have major impact on all biotic communities that still needs to be fully elucidated. In addition, as these organisms represent a highly divergent evolutionary lineage, they are strikingly different from the common 'model system' eukaryotes, such as some mammals, plants or fungi. A number of excellent reviews, published over the past decade, were dedicated to specialized topics from the areas of trypanosomatid molecular and cell biology, biochemistry, host-parasite relationships or other aspects of these fascinating organisms. However, there is a need for a more comprehensive review that summarizing recent advances in the studies of trypanosomatids in the last 30 years, a task, which we tried to accomplish with the current paper.

Evolutionary instability of CUG-Leu in the genetic code of budding yeasts

The genetic code used in nuclear genes is almost universal, but here we report that it changed three times in parallel during the evolution of budding yeasts. All three changes were reassignments of the codon CUG, which is translated as serine (in 2 yeast clades), alanine (1 clade), or the ‘universal’ leucine (2 clades). The newly discovered Ser2 clade is in the final stages of a genetic code transition. Most species in this clade have genes for both a novel tRNA^Ser(CAG) and an ancestral tRNA^Leu(CAG) to read CUG, but only tRNA^Ser(CAG) is used in standard growth conditions. The coexistence of these alloacceptor tRNA genes indicates that the genetic code transition occurred via an ambiguous translation phase. We propose that the three parallel reassignments of CUG were not driven by natural selection in favor of their effects on the proteome, but by selection to eliminate the ancestral tRNA^Leu(CAG).

The genetic code for amino acids is nearly universal, and among eukaryotic nuclear genomes the only known reassignments are of codon CUG in yeasts. Here, the authors identify a third independent CUG transition in budding yeasts that is still ongoing with alternative tRNAs present in the genome.

Drugging tRNA aminoacylation

Inhibition of tRNA aminoacylation has proven to be an effective antimicrobial strategy, impeding an essential step of protein synthesis. Mupirocin, the well-known selective inhibitor of bacterial isoleucyl-tRNA synthetase, is one of three aminoacylation inhibitors now approved for human or animal use. However, design of novel aminoacylation inhibitors is complicated by the steadfast requirement to avoid off-target inhibition of protein synthesis in human cells. Here we review available data regarding known aminoacylation inhibitors as well as key amino-acid residues in aminoacyl-tRNA synthetases (aaRSs) and nucleotides in tRNA that determine the specificity and strength of the aaRS-tRNA interaction. Unlike most ligand-protein interactions, the aaRS-tRNA recognition interaction represents coevolution of both the tRNA and aaRS structures to conserve the specificity of aminoacylation. This property means that many determinants of tRNA recognition in pathogens have diverged from those of humans-a phenomenon that provides a valuable source of data for antimicrobial drug development.

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.

Nuclear codon reassignments in the genomics era and mechanisms behind their evolution

The canonical genetic code ubiquitously translates nucleotide into peptide sequence with several alterations known in viruses, bacteria, mitochondria, plastids, and single-celled eukaryotes. A new hypothesis to explain genetic code changes, termed tRNA loss driven codon reassignment, has been proposed recently when the polyphyly of the yeast codon reassignment events has been uncovered. According to this hypothesis, the driving force for genetic code changes are tRNA or translation termination factor loss-of-function mutations or loss-of-gene events. The free codon can subsequently be captured by all tRNAs that have an appropriately mutated anticodon and are efficiently charged. Thus, codon capture most likely happens by near-cognate tRNAs and tRNAs whose anticodons are not part of the recognition sites of the respective aminoacyl-tRNA-synthetases. This hypothesis comprehensively explains the CTG codon translation as alanine in Pachysolen yeast together with the long known translation of the same codon as serine in Candida albicans and related species, and can also be applied to most other known reassignments.

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

BACKGROUND

Departures from the standard genetic code in eukaryotic nuclear genomes are known for only a handful of lineages and only a few genetic code variants seem to exist outside the ciliates, the most creative group in this regard. Most frequent code modifications entail reassignment of the UAG and UAA codons, with evidence for at least 13 independent cases of a coordinated change in the meaning of both codons. However, no change affecting each of the two codons separately has been documented, suggesting the existence of underlying evolutionary or mechanistic constraints.

RESULTS

Here, we present the discovery of two new variants of the nuclear genetic code, in which UAG is translated as an amino acid while UAA is kept as a termination codon (along with UGA). The first variant occurs in an organism noticed in a (meta)transcriptome from the heteropteran Lygus hesperus and demonstrated to be a novel insect-dwelling member of Rhizaria (specifically Sainouroidea). This first documented case of a rhizarian with a non-canonical genetic code employs UAG to encode leucine and represents an unprecedented change among nuclear codon reassignments. The second code variant was found in the recently described anaerobic flagellate Iotanema spirale (Metamonada: Fornicata). Analyses of transcriptomic data revealed that I. spirale uses UAG to encode glutamine, similarly to the most common variant of a non-canonical code known from several unrelated eukaryotic groups, including hexamitin diplomonads (also a lineage of fornicates). However, in these organisms, UAA also encodes glutamine, whereas it is the primary termination codon in I. spirale. Along with phylogenetic evidence for distant relationship of I. spirale and hexamitins, this indicates two independent genetic code changes in fornicates.

CONCLUSIONS

Our study documents, for the first time, that evolutionary changes of the meaning of UAG and UAA codons in nuclear genomes can be decoupled and that the interpretation of the two codons by the cytoplasmic translation apparatus is mechanistically separable. The latter conclusion has interesting implications for possibilities of genetic code engineering in eukaryotes. We also present a newly developed generally applicable phylogeny-informed method for inferring the meaning of reassigned codons.