Mutation and selection on the wobble nucleotide in tRNA anticodons in marine bivalve mitochondrial genomes

Background

Animal mitochondrial genomes typically encode one tRNA for each synonymous codon family, so that each tRNA anticodon essentially has to wobble to recognize two or four synonymous codons. Several factors have been hypothesized to determine the nucleotide at the wobble site of a tRNA anticodon in mitochondrial genomes, such as the codon-anticodon adaptation hypothesis, the wobble versatility hypothesis, the translation initiation and elongation conflict hypothesis, and the wobble cost hypothesis.

Principal Findings

In this study, we analyzed codon usage and tRNA anticodon wobble sites of 29 marine bivalve mitochondrial genomes to evaluate features of the wobble nucleotides in tRNA anticodons. The strand-specific mutation bias favors G and T on the H strand in all the 29 marine bivalve mitochondrial genomes. A bias favoring G and T is also visible in the third codon positions of protein-coding genes and the wobble sites of anticodons, rejecting that codon usage bias drives the wobble sites of tRNA anticodons or tRNA anticodon bias drives the evolution of codon usage. Almost all codon families (98.9%) from marine bivalve mitogenomes support the wobble versatility hypothesis. There are a few interesting exceptions involving tRNA^Trp with an anticodon CCA fixed in Pectinoida species, tRNA^Ser with a GCU anticodon fixed in Mytiloida mitogenomes, and the uniform anticodon CAU of tRNA^Met translating the AUR codon family.

Conclusions/Significance

These results demonstrate that most of the nucleotides at the wobble sites of tRNA anticodons in marine bivalve mitogenomes are determined by wobble versatility. Other factors such as the translation initiation and elongation conflict, and the cost of wobble translation may contribute to the determination of the wobble nucleotide in tRNA anticodons. The finding presented here provides valuable insights into the previous hypotheses of the wobble nucleotide in tRNA anticodons by adding some new evidence.     

Identification and codon reading properties of 5-cyanomethyl uridine,a new modified nucleoside found in the anticodon wobble position of mutant haloarchaeal isoleucine tRNAs

Most archaea and bacteria use a modified C in the anticodon wobble position of isoleucine tRNA to base pair with A but not with G of the mRNA. This allows the tRNA to read the isoleucine codon AUA without also reading the methionine codon AUG. To understand why a modified C, and not U or modified U, is used to base pair with A, we mutated the C34 in the anticodon of Haloarcula marismortui isoleucine tRNA (tRNA₂^Ile) to U, expressed the mutant tRNA in Haloferax volcanii, and purified and analyzed the tRNA. Ribosome binding experiments show that although the wild-type tRNA₂^Ile binds exclusively to the isoleucine codon AUA, the mutant tRNA binds not only to AUA but also to AUU, another isoleucine codon, and to AUG, a methionine codon. The G34 to U mutant in the anticodon of another H. marismortui isoleucine tRNA species showed similar codon binding properties. Binding of the mutant tRNA to AUG could lead to misreading of the AUG codon and insertion of isoleucine in place of methionine. This result would explain why most archaea and bacteria do not normally use U or a modified U in the anticodon wobble position of isoleucine tRNA for reading the codon AUA. Biochemical and mass spectrometric analyses of the mutant tRNAs have led to the discovery of a new modified nucleoside, 5-cyanomethyl U in the anticodon wobble position of the mutant tRNAs. 5-Cyanomethyl U is present in total tRNAs from euryarchaea but not in crenarchaea, eubacteria, or eukaryotes.     

Transfer RNA Identity Change in Anticodon Variants of E. coli tRNAPhe in Vivo

The anticodon sequence is a major recognition element for most aminoacyl-tRNA synthetases. We investigated the in vivo effects of changing the anticodon on the aminoacylation specificity in the example of E. coli tRNA^Phe. Constructing different anticodon mutants of E. coli tRNA^Phe by site-directed mutagenesis, we isolated 22 anticodon mutant tRNA^Phe; the anticodons corresponded to 16 amino acids and an opal stop codon. To examine whether the mutant tRNAs had changed their amino acid acceptor specificity in vivo, we tested the viability of E. coli strains containing these tRNA^Phe genes in a medium which permitted tRNA induction. Fourteen mutant tRNA genes did not affect host viability. However, eight mutant tRNA genes were toxic to the host and prevented growth, presumably because the anticodon mutants led to translational errors. Many mutant tRNAs which did not affect host viability were not aminoacylated in vivo. Three mutant tRNAs containing anticodon sequences corresponding to lysine (UUU), methionine (CAU) and threonine (UGU) were charged with the amino acid corresponding to their anticodon, but not with phenylalanine. These three tRNAs and tRNA^Phe are located in the same cluster in a sequence similarity dendrogram of total E. coli tRNAs. The results support the idea that such tRNAs arising from in vivo evolution are derived by anticodon change from the same ancestor tRNA.     

Discovery and characterization of tRNA lysidine synthetase (TilS)

In the bacterial decoding system, the AUA codon is deciphered as isoleucine by tRNA^Ile bearing lysidine (L, 2-lysyl-cytidine) at the wobble position. Lysidine is an essential modification that determines both the codon and amino acid specificities of tRNA^Ile. We identified an enzyme named tRNA^Ile lysidine synthetase (TilS) that catalyzes lysidine formation by using lysine and ATP as substrates. Biochemical studies revealed a molecular mechanism of lysidine formation that consists of two consecutive reactions involving the adenylated tRNA intermediate. In addition, we deciphered how Escherichia coli TilS specifically discriminates between tRNA^Ile and the structurally similar tRNA^Met, which bears the same anticodon loop. Recent structural studies unveiled tRNA recognition by TilS, and a molecular basis of lysidine formation at atomic resolution.     

Functional elucidation of a key contact between tRNA and the large ribosomal subunit rRNA during decoding

The selection of cognate tRNAs during translation is specified by a kinetic discrimination mechanism driven by distinct structural states of the ribosome. While the biochemical steps that drive the tRNA selection process have been carefully documented, it remains unclear how recognition of matched codon:anticodon helices in the small subunit facilitate global rearrangements in the ribosome complex that efficiently promote tRNA decoding. Here we use an in vitro selection approach to isolate tRNA^Trp miscoding variants that exhibit a globally perturbed tRNA tertiary structure. Interestingly, the most substantial distortions are positioned in the elbow region of the tRNA that closely approaches helix 69 (H69) of the large ribosomal subunit. The importance of these specific interactions to tRNA selection is underscored by our kinetic analysis of both tRNA and rRNA variants that perturb the integrity of this interaction.     

The Nature of the Nucleotide at the 5' Side of the tRNA Su9 Anticodon Affects UGA Suppression in Escherichia coli

In Escherichia coli a UGA codon can be efficiently suppressedby a suppressor tRNA^Trp called Su9. Here, we show that the levelof UGA suppression is determined by the nature of the nucleotideat the 5' side of the anticodon of the suppressor (position33). UGA suppression occurs when a pyrimidine residue is locatedin position 33 of the tRNA, and suppression is more efficientwith a U than with a C in this position. On the other hand,when a purine residue is located at this position UGA suppressionis extremely low. These results show that in the case of tRNASu9, the UGA codon context effect does not require base pairingbetween the nucleotide at the 3' side of the codon and the 5'side of the anticodon.     

Translation of the UGA triplet in vitro by tryptophan transfer RNA's

Tryptophan transfer RNA from the UGA-suppressing strain of Escherichia coli CAJ64 was purified and assayed for suppressor activity in vitro in two ways: by translation of the bacteriophage T4 lysozyme messenger RNA bearing a UGA mutation, and by translation of poly(U-G-A). Purified tRNA^Trp, and no other fraction, stimulates lysozyme synthesis 30-fold above the level seen when comparable amounts of tryptophan tRNA from the non-suppressing strain, CA244, were added; it also translates poly(U-G-A) as polytryptophan more efficiently than the su⁻ tRNA. Tryptophan tRNA from the non-suppressing strain is active in the assays but far less so than CAJ64 tRNA^Trp, and this is consistent with the leakiness of su⁻ strains. Since the nucleotide sequences of these tryptophan tRNA's are known (Hirsh, 1971), it is concluded that tRNA with a CCA anticodon recognizes the UGA triplet and this recognition is improved by a nucleotide change elsewhere in the molecule.     

Transfer RNA-like structure of the human Alu family: Implications of its generation mechanism and possible functions

Summary Structural resemblance of the human Alu family with a subset of vertebrate tRNAs was detected. Of four tRNAs, tRNA^Lys, tRNA^Ile, tRNA^Thr, and tRNA^Tyr, which comprise a structurally related family, tRNA^Lys is the most similar to the human Alu family. Of the 76 nucleotides in lysine tRNA (including the CCA tail), 47 are similar to the human Alu family (60% identity). The secondary structure of the human Alu family corresponding to the D-stem and anticodon stem regions of the tRNA appears to be very stable. The 7SL RNA, which is a progenitor of the human Alu family, is less similar to lysine tRNA (55% identity), and the secondary structure of the 7SL RNA folded like a tRNA is less stable than that of the human Alu family folded likewise. Insertion of the tetranucleotide GAGA, which is an important region of the second promoter for RNA polymerase III in the Alu sequence, occurred during the deletion and ligation process to generate the Alu sequence from the parental 7SL RNA. These results suggest that the human Alu family was generated from the 7SL RNA by deletion, insertion, and mutations, which thus modified the ancestral 7SL sequence so that it could form a structure more closely resembling lysine tRNA. The similarities of several short interspersed sequences to the lysine tRNA were also examined. TheGalago type 2 family, which was reported to be derived from a methionine initiator tRNA, was also found to be similar to the lysine tRNA. Thus lysine tRNA-like structures are widespread in genomes in the animal kingdom. The implications of these findings in relation to the mechanism of generation of the human Alu family and its possible functions are discussed.     

Nucleotide sequence of a cucumber chloroplast proline tRNA

The nucleotide sequence of a proline tRNA (anticodon UGG) from cucumber chloroplasts has been determined. The sequence is: pAAGGAUGUAGCGCAGCUUCA-DAGCGCAψUUGUUUUGGNψFACAAAAUmsu7GUCACGGGTψCAAAUCCUGUCAUCCUUACCA_OH. It shows 93% homology with spinach chloroplast tRNA^Pro (UGG) and 72% homology with bean mitochondrial tRNA^Pro (UGG), the other two known plant organellar tRNAs^Pro.     

A single mutational modification of a tryptophan-specific transfer RNA permits aminoacylation by glutamine and translation of the codon UAG

In this work we show that the wild-type (su^?7) progenitor of the recessivelethal suppressors of UAG (su⁺7(UAG)) and of UAA/G (su⁺7(UAA/G)) is the structural gene for transfer RNA^Trp, the adaptor for translating the codon UGG. The su⁺7(UAG) suppressor form of the tRNA has a C for U substitution in the middle base of the anticodon; in the su⁺7(UAA/G) suppressor tRNA both C residues of the anticodon are replaced by U. Our data establish that the mutational change altering the tRNA^Trp to a UAG suppressor is accompanied by a loss of tryptophan-accepting specificity and the acquisition of glutamine-acceptor activity.     

Mitochondrial Genome of <Emphasis Type="Italic">Ciona savignyi</Emphasis> (Urochordata,Ascidiacea, Enterogona): Comparison of Gene Arrangement and tRNA Genes with <Emphasis Type="Italic">Halocynthia roretzi</Emphasis> Mitochondrial Genome

The complete nucleotide sequence of the urochordate Ciona savignyi (Ascidiacea, Enterogona) mitochondrial (mt) genome (14,737 bp) was determined. The Ciona mt genome does not encode a gene for ATP synthetase subunit 8 but encodes an additional tRNA^Gly gene (anticodon UCU), as is the case in another urochordate, Halocynthia roretzi (Ascidiacea, Pleurogona), mt genome. In addition, the Ciona mt genome encodes two tRNA^Met genes; anticodon CAT and anticodon TAT. The tRNA^Cys gene is thought to lack base pairs at the D-stem. Thus, the Ciona mt genome encodes 12 protein, 2 rRNA, and 24 tRNA genes. The gene arrangement of the Ciona mt genome differs greatly from those of any other metazoan mt genomes reported to date. Only three gene boundaries are shared between the Halocynthia and the Ciona mt genomes. Molecular phylogenetic analyses based on amino acid sequences of mt protein genes failed to demonstrate the monophyly of the chordates.     

Chromosomal mapping of tRNA genes from Dictyostelium discoideum

Summary Different wild-type isolates of Dictyostelium discoideum exhibit extensive polymorphism in the length of restriction fragments carrying tRNA genes. These size differences were used to study the organisation of two tRNA gene families which encode a tRNA^Val(GUU) and a tRNA^Val(GUA) gene. The method used involved a combination of classitics. The tRNA genes were mapped to specific linkage groups (chromosomes) by correlating the presence of polymorphic DNA bands that hybridized with the tRNA gene probes with the presence of genetic markers for those linkage groups. These analyses established that both of the tRNA gene families are dispersed among sites on several of the chromosomes. Information of nine tRNA^Val(GUU) genes from the wild-type isolate NC4 was obtained: three map to linkage group I (C, E, F,), two map to linkage group II (D, I), one maps to linkage group IV (G), one, which corresponds to the cloned gene, maps to either linkage group III or VI (B), and two map to one of linkage groups III, VI or VIII (A, H). Six tRNA^Val(GUA) genes from the NC4 isolate were mapped; one to linkage group I (D), two to linkage group III, VI or VII (B, C) and three to linkage group VII or III (A, E, F).     

Acceptor stem and anticodon RNA hairpin helix interactions with glutamine tRNA synthetase

The class I glutamine (Gln) tRNA synthetase interacts with the anticodon and acceptor stem of glutamine tRNA. RNA hairpin helices were designed to probe acceptor stem and anticodon stem-loop contacts. A seven-base pair RNA microhelix derived from the acceptor stem of tRNA^Gln was aminoacylated by Gln tRNA synthetase. Variants of the glutamine acceptor stem microhelix implicated the discriminator base as a major identity element for glutaminylation of the RNA helix. A second RNA microhelix representing the anticodon stem-loop competitively inhibited tRNA^Gln charging. However, the anticodon stem-loop microhelix did not enhance aminoacylation of the acceptor stem microhelix. Thus, transduction of the anticodon identity signal may require covalent continuity of the tRNA chain to trigger efficient aminoacylation.     

Complete mitochondrial genome of a carabid beetle,Damaster mirabilissimus mirabilissim (Coleoptera: Carabidae)

In the present study, we report the 16 823‐bp long complete mitochondrial genome (mitogenome) of a carabid beetle, Damaster mirabilissimus mirabilissim (Coleoptera: Carabidae), which is endangered in Korea. The gene arrangement of D. m. mirabilissim mitogenome is identical to the most common type found in insects. The start codon of the D. m. mirabilissim COI gene is a typical ATN codon. On the other hand, the initiation codon for ND1 gene is TTG, instead of ATN. All transfer RNAs (tRNAs) exhibit a stable canonical clover‐leaf structure, except for tRNA^Ser(AGN), the dihydrouridine arm of which forms a simple loop. The 1703‐bp long A+T‐rich region is the second longest among the complete adephagan mitogenome sequences, next to Macrogyrus oblongus belonging to Gyrinoidea. One of the unusual features of the genome is the presence of a tRNA^Leu(UUR)‐like sequence in the A+T‐rich region. This sequence displays the proper anticodon sequence and the potential to form secondary structures, but also harbors many mismatches in the stems.     

Gene rearrangements in gekkonid mitochondrial genomes with shuffling,loss, and reassignment of tRNA genes

Background

Vertebrate mitochondrial genomes (mitogenomes) are 16–18 kbp double-stranded circular DNAs that encode a set of 37 genes. The arrangement of these genes and the major noncoding region is relatively conserved through evolution although gene rearrangements have been described for diverse lineages. The tandem duplication-random loss model has been invoked to explain the mechanisms of most mitochondrial gene rearrangements. Previously reported mitogenomic sequences for geckos rarely included gene rearrangements, which we explore in the present study.

Results

We determined seven new mitogenomic sequences from Gekkonidae using a high-throughput sequencing method. The Tropiocolotes tripolitanus mitogenome involves a tandem duplication of the gene block: tRNA^Arg, NADH dehydrogenase subunit 4L, and NADH dehydrogenase subunit 4. One of the duplicate copies for each protein-coding gene may be pseudogenized. A duplicate copy of the tRNA^Arg gene appears to have been converted to a tRNA^Gln gene by a C to T base substitution at the second anticodon position, although this gene may not be fully functional in protein synthesis. The Stenodactylus petrii mitogenome includes several tandem duplications of tRNA^Leu genes, as well as a translocation of the tRNA^Ala gene and a putative origin of light-strand replication within a tRNA gene cluster. Finally, the Uroplatus fimbriatus and U. ebenaui mitogenomes feature the apparent loss of the tRNA^Glu gene from its original position. Uroplatus fimbriatus appears to retain a translocated tRNA^Glu gene adjacent to the 5’ end of the major noncoding region.

Conclusions

The present study describes several new mitochondrial gene rearrangements from Gekkonidae. The loss and reassignment of tRNA genes is not very common in vertebrate mitogenomes and our findings raise new questions as to how missing tRNAs are supplied and if the reassigned tRNA gene is fully functional. These new examples of mitochondrial gene rearrangements in geckos should broaden our understanding of the evolution of mitochondrial gene arrangements.

Electronic supplementary material

The online version of this article (doi:10.1186/1471-2164-15-930) contains supplementary material, which is available to authorized users.     

Overexpressed mitochondrial leucyl-tRNA synthetase suppresses the A3243G mutation in the mitochondrial tRNA(Leu(UUR)) gene

The A3243G mutation in the human mitochondrial tRNA^Leu(UUR) gene causes a number of human diseases. This mutation reduces the level and fraction of aminoacylated tRNA^Leu(UUR) and eliminates nucleotide modification at the wobble position of the anticodon. These deficiencies are associated with mitochondrial translation defects that result in decreased levels of mitochondrial translation products and respiratory chain enzyme activities. We have suppressed the respiratory chain defects in A3243G mutant cells by overexpressing human mitochondrial leucyl-tRNA synthetase. The rates of oxygen consumption in suppressed cells were directly proportional to the levels of leucyl-tRNA synthetase. Fifteenfold higher levels of leucyl-tRNA synthetase resulted in wild-type respiratory chain function. The suppressed cells had increased steady-state levels of tRNA^Leu(UUR) and up to threefold higher steady-state levels of mitochondrial translation products, but did not have rates of protein synthesis above those in parental mutant cells. These data suggest that suppression of the A3243G mutation occurred by increasing protein stability. This suppression of a tRNA gene mutation by increasing the steady-state levels of its cognate aminoacyl-tRNA synthetase is a model for potential therapies for human pathogenic tRNA mutations.