Evolution of the mitochondrial genetic code. III. Reassignment of CUN codons from leucine to threonine during evolution of yeast mitochondria

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.

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.

An unusual tRNAThr derived from tRNAHis reassigns in yeast mitochondria the CUN codons to threonine

The standard genetic code is used by most living organisms, yet deviations have been observed in many genomes, suggesting that the genetic code has been evolving. In certain yeast mitochondria, CUN codons are reassigned from leucine to threonine, which requires an unusual tRNA(Thr) with an enlarged 8-nt anticodon loop ( ). To trace its evolutionary origin we performed a comprehensive phylogenetic analysis which revealed that evolved from yeast mitochondrial tRNA(His). To understand this tRNA identity change, we performed mutational and biochemical experiments. We show that Saccharomyces cerevisiae mitochondrial threonyl-tRNA synthetase (MST1) could attach threonine to both and the regular , but not to the wild-type tRNA(His). A loss of the first nucleotide (G(-1)) in tRNA(His) converts it to a substrate for MST1 with a K(m) value (0.7 μM) comparable to that of (0.3 μM), and addition of G(-1) to allows efficient histidylation by histidyl-tRNA synthetase. We also show that MST1 from Candida albicans, a yeast in which CUN codons remain assigned to leucine, could not threonylate , suggesting that MST1 has coevolved with . Our work provides the first clear example of a recent recoding event caused by alloacceptor tRNA gene recruitment.

Extensive frameshift at all AGG and CCC codons in the mitochondrial cytochrome c oxidase subunit 1 gene of Perkinsus marinus (Alveolata; Dinoflagellata)

Diverse mitochondrial (mt) genetic systems have evolved independently of the more uniform nuclear system and often employ modified genetic codes. The organization and genetic system of dinoflagellate mt genomes are particularly unusual and remain an evolutionary enigma. We determined the sequence of full-length cytochrome c oxidase subunit 1 (cox1) mRNA of the earliest diverging dinoflagellate Perkinsus and show that this gene resides in the mt genome. Apparently, this mRNA is not translated in a single reading frame with standard codon usage. Our examination of the nucleotide sequence and three-frame translation of the mRNA suggest that the reading frame must be shifted 10 times, at every AGG and CCC codon, to yield a consensus COX1 protein. We suggest two possible mechanisms for these translational frameshifts: a ribosomal frameshift in which stalled ribosomes skip the first bases of these codons or specialized tRNAs recognizing non-triplet codons, AGGY and CCCCU. Regardless of the mechanism, active and efficient machinery would be required to tolerate the frameshifts predicted in Perkinsus mitochondria. To our knowledge, this is the first evidence of translational frameshifts in protist mitochondria and, by far, is the most extensive case in mitochondria.

The mechanisms of codon reassignments in mitochondrial genetic codes

Many cases of nonstandard genetic codes are known in mitochondrial genomes. We carry out analysis of phylogeny and codon usage of organisms for which the complete mitochondrial genome is available, and we determine the most likely mechanism for codon reassignment in each case. Reassignment events can be classified according to the gain-loss framework. The “gain” represents the appearance of a new tRNA for the reassigned codon or the change of an existing tRNA such that it gains the ability to pair with the codon. The “loss” represents the deletion of a tRNA or the change in a tRNA so that it no longer translates the codon. One possible mechanism is codon disappearance (CD), where the codon disappears from the genome prior to the gain and loss events. In the alternative mechanisms the codon does not disappear. In the unassigned codon mechanism, the loss occurs first, whereas in the ambiguous intermediate mechanism, the gain occurs first. Codon usage analysis gives clear evidence of cases where the codon disappeared at the point of the reassignment and also cases where it did not disappear. CD is the probable explanation for stop to sense reassignments and a small number of reassignments of sense codons. However, the majority of sense-to-sense reassignments cannot be explained by CD. In the latter cases, by analysis of the presence or absence of tRNAs in the genome and of the changes in tRNA sequences, it is sometimes possible to distinguish between the unassigned codon and the ambiguous intermediate mechanisms. We emphasize that not all reassignments follow the same scenario and that it is necessary to consider the details of each case carefully.

A comparative genomics analysis of codon reassignments reveals a link with mitochondrial proteome size and a mechanism of genetic code change via suppressor tRNAs

Using a comparative genomics approach we demonstrate a negative correlation between the number of codon reassignments undergone by 222 mitochondrial genomes and the mitochondrial genome size, the number of mitochondrial ORFs, and the sizes of the large and small subunit mitochondrial rRNAs. In addition, we show that the TGA-to-tryptophan codon reassignment, which has occurred 11 times in mitochondrial genomes, is found in mitochondrial genomes smaller than those which have not undergone the reassignment. We therefore propose that mitochondrial codon reassignments occur in a wide range of phyla, particularly in Metazoa, due to a reduced "proteomic constraint" on the mitochondrial genetic code, compared to the nuclear genetic code. The reduced proteomic constraint reflects the small size of the mitochondrial-encoded proteome and allows codon reassignments to occur with less likelihood of lethality. In addition, we demonstrate a striking link between nonsense codon reassignments and the decoding properties of naturally occurring nonsense suppressor tRNAs. This suggests that natural preexisting nonsense suppression facilitated nonsense codon reassignments and constitutes a novel mechanism of genetic code change. These findings explain for the first time the identity of the stop codons and amino acids reassigned in mitochondrial and nuclear genomes. Nonsense suppressor tRNAs provided the raw material for nonsense codon reassignments, implying that the properties of the tRNA anticodon have dictated the identity of nonsense codon reassignments.

Comparing two evolutionary mechanisms of modern tRNAs

All modern tRNA gene families have a high similarity in their primary structure, and share the same cloverleaf secondary structure and an inverted L tertiary structure, which provide the clues for the study of their origin and evolution. There are two important mechanisms of the tRNA sequences evolution. One is point mutation, another is complementary duplication method. Both of them are supported by some specific examples. To find out the superior one of the two mechanisms or find out the most suitable mechanism for modern tRNAs evolution, we constructed two types of networks, parallel and antiparallel networks, based on the two mechanisms respectively, and then compared the degree distribution and clustering coefficient of networks constructed by the tRNA sequences of the single anticodon group, single isoaccepting group, and the whole tRNAs group of the two types of networks. The result of the comparison seems consistent with the idea that modern tRNA sequences evolved primarily by the mechanism of complementary method, and point mutation is an important and indispensable auxiliary mechanism during the evolutionary event.

Comparative evolutionary genomics unveils the molecular mechanism of reassignment of the CTG codon in Candida spp

Using the (near) complete genome sequences of the yeasts Candida albicans, Saccharomyces cerevisiae, and Schizosaccharomyces pombe, we address the evolution of a unique genetic code change, which involves decoding of the standard leucine-CTG codon as serine in Candida spp. By using two complementary comparative genomics approaches, we have been able to shed new light on both the origin of the novel Candida spp. Ser-tRNA(CAG), which has mediated CTG reassignment, and on the evolution of the CTG codon in the genomes of C. albicans, S. cerevisiae, and S. pombe. Sequence analyses of newly identified tRNAs from the C. albicans genome demonstrate that the Ser-tRNA(CAG) is derived from a serine and not a leucine tRNA in the ancestor yeast species and that this codon reassignment occurred approximately 170 million years ago, but the origin of the Ser-tRNA(CAG) is more ancient, implying that the ancestral Leu-tRNA that decoded the CTG codon was lost after the appearance of the Ser-tRNA(CAG). Ambiguous CTG decoding by the Ser-tRNA(CAG) combined with biased AT pressure forced the evolution of CTG into TTR codons and have been major forces driving evolution of the CTN codon family in C. albicans. Remarkably, most of the CTG codons present in extant C. albicans genes are encoded by serine and not leucine codons in homologous S. cerevisiae and S. pombe genes, indicating that a significant number of serine TCN and AGY codons evolved into CTG codons either directly by simultaneous double mutations or indirectly through an intermediary codon. In either case, CTG reassignment had a major impact on the evolution of the coding component of the Candida spp. genome.

The evolutionary change of the genetic code as restricted by the anticodon and identity of transfer RNA

The discovery of non-universal genetic codes in several mitochondria and nuclear systems during the part ten years has necessitated a reconsideration of the concept that the genetic code is universal and frozen, as was once believed. Here, the flexibility of the relationship between codons and amino acids is discussed on the basis of the distribution of non-universal genetic codes in various organisms insofar as has been observed to date. Judging from the result of recent investigations into tRNA identity, it would appear that the non-participation of the anticodon in recognition by aminoacyl-tRNA synthetase has significantly influenced the variability of codons.

Experimental studies on the origin of the genetic code and the process of protein synthesis: a review update

This article is an update of our earlier review (Lacey and Mullins, 1983) in this journal on the origin of the genetic code and the process of protein synthesis. It is our intent to discuss only experimental evidence published since then although there is the necessity to mention the old enough to place the new in context. We do not include theoretical nor hypothetical treatments of the code or protein synthesis. Relevant data regarding the evolution of tRNAs and the recognition of tRNAs by aminoacyl-tRNA-synthetases are discussed. Our present belief is that the code arose based on a core of early assignments which were made on a physico-chemical and anticodonic basis and this was expanded with new assignments later. These late assignments do not necessarily show an amino acid-anticodon relatedness. In spite of the fact that most data suggest a code origin based on amino acid-anticodon relationships, some new data suggesting preferential binding of Arg to its codons are discussed. While information regarding coding is not increasing very rapidly, information regarding the basic chemistry of the process of protein synthesis has increased significantly, principally relating to aminoacylation of mono- and polyribonucleotides. Included in those studies are several which show stereoselective reactions of L-amino acids with nucleotides having D-sugars. Hydrophobic interactions definitely play a role in the preferences which have been observed.

Recent evidence for evolution of the genetic code

The genetic code, formerly thought to be frozen, is now known to be in a state of evolution. This was first shown in 1979 by Barrell et al. (G. Barrell, A. T. Bankier, and J. Drouin, Nature [London] 282:189-194, 1979), who found that the universal codons AUA (isoleucine) and UGA (stop) coded for methionine and tryptophan, respectively, in human mitochondria. Subsequent studies have shown that UGA codes for tryptophan in Mycoplasma spp. and in all nonplant mitochondria that have been examined. Universal stop codons UAA and UAG code for glutamine in ciliated protozoa (except Euplotes octacarinatus) and in a green alga, Acetabularia. E. octacarinatus uses UAA for stop and UGA for cysteine. Candida species, which are yeasts, use CUG (leucine) for serine. Other departures from the universal code, all in nonplant mitochondria, are CUN (leucine) for threonine (in yeasts), AAA (lysine) for asparagine (in platyhelminths and echinoderms), UAA (stop) for tyrosine (in planaria), and AGR (arginine) for serine (in several animal orders) and for stop (in vertebrates). We propose that the changes are typically preceded by loss of a codon from all coding sequences in an organism or organelle, often as a result of directional mutation pressure, accompanied by loss of the tRNA that translates the codon. The codon reappears later by conversion of another codon and emergence of a tRNA that translates the reappeared codon with a different assignment. Changes in release factors also contribute to these revised assignments. We also discuss the use of UGA (stop) as a selenocysteine codon and the early history of the code.

Cloning and analysis of five mitochondrial tRNA-encoding genes from the fungus Beauveria bassiana

Five mitochondrial (mt) tRNA genes from the filamentous fungus, Beauveria bassiana, were cloned and sequenced. The genes encoding the Val-, Ile-, Ser-, Trp- and Pro-accepting tRNAs were found clustered in the region 5' to the lrRNA-encoding gene. The genes were 64-77% homologous with the equivalent genes from other filamentous fungi, 49-58% to yeasts with the exception of the Val-accepting tRNA-encoding gene which was 76%, and only slightly homologous with Escherichia coli. The B. bassiana mt genetic code was found to be similar to that of other fungal mitochondria in that the UGA codon is used as a signal for Trp rather than as a stop codon. Transcript analysis has revealed that the genes present in tRNA cluster are transcribed and processed into tRNA-size products. Secondary structure models proposed for the gene products show that conservation of tRNA secondary structure also exists. The presence of a GGC sequence rather than a GGU sequence in the D-loop of the tRNA(Trp)-encoding gene is a feature unique to the B. bassiana mt tRNA. An unconventional G-A base pair present in the D-stem of the tRNA(Ser)-encoding gene is a feature conserved in the mt tRNA of other filamentous fungi. Comparison of the B. bassiana tRNA-encoding genes with those of two other filamentous fungi and two yeasts revealed that the differences between closely related species favoured transition-type mutations.

Isolation of the ATP synthase subunit 6 and sequence of the mitochondrial ATP6 gene of the yeast Candida parapsilosis

The mitochondrially translated product called subunit 6 was extracted from the yeast Candida parapsilosis mitochondria using an organic solvent mixture and purified by reverse-phase HPLC. The partial N-terminal sequence of subunit 6 reveals a post-translational cleavage site as in Saccharomyces cerevisiae. The structural mitochondrial gene ATP6 was isolated form a mitochondrial DNA library using the oligonucleotide probe procedure. The gene and the surrounding regions were cloned into M13tg130 and M13tg131 phage vectors. The insert contained an open reading frame 738-bp encoding a 246-amino-acid polypeptide. Mature subunit 6 contains 243 amino acid residues and the predicted molecular mass is 26,511 Da. The subunit shows 52% similarity with ATP synthase subunit 6 of the yeast S. cerevisiae. Comparison between protein and DNA sequences shows that the CUN codon family codes for a leucine in C. parapsilosis mitochondria.

Genetic code 1990. Outlook

The genetic code is evolving as shown by 9 departures from the universal code: 6 of them are in mitochondria and 3 are in nuclear codes. We propose that these changes are preceded by disappearance of a codon from coding sequences in mRNA of an organism or organelle. The function of the codon that disappears is taken by other, synonymous codons, so that there is no change in amino acid sequences of proteins. The deleted codon then reappears with a new function. Wobble pairing between anticodons and codons has evolved, starting with a single UNN anticodon pairing with 4 codons. Directional mutation pressure affects codon usage and may produce codon reassignments, especially of stop codons. Selenocysteine is coded by UGA, which is also a stop codon, and this anomaly is discussed. The outlook for discovery of more changes in the code is favorable, and open reading frames should be compared with actual sequential analyses of protein molecules in this search.

The genetic code in mitochondria and chloroplasts

The universal genetic code is used without changes in chloroplasts and in mitochondria of green plants. Non-plant mitochondria use codes that include changes from the universal code. Chloroplasts use 31 anticodons in translating the code; a number smaller than that used by bacteria, because chloroplasts have eliminated 10 CNN anticodons that are found in bacteria. Green plant mitochondria (mt) obtain some tRNAs from the cytosol, and genes for some other tRNAs have been acquired from chloroplast DNA. The code in non-plant mt differs from the universal code in the following usages found in various organisms: UGA for Trp, AUA for Met, AGR for Ser and stop, AAA for Asn, CUN for Thr, and possibly UAA for Tyr. CGN codons are not used by Torulopsis yeast mt. Non-plant mt, e.g. in vertebrates, may use a minimum of 22 anticodons for complete translation of mRNA sequences. The following possible causes are regarded as contributing to changes in the non-plant mt: directional mutation pressure, genomic economization, changes in charging specificity of tRNAs, loss of release factor RF2, changes in RF1, changes in anticodons, loss of lysidine-forming enzyme system, and disappearance of codons from coding sequences.