Mitochondrial tRNA Sequences Research Articles

Disease-associated mutations in mitochondrial precursor tRNAs affect binding, m1R9 methylation, and tRNA processing by mtRNase P.

Mitochondrial diseases linked to mutations in mitochondrial (mt) tRNA sequences are common. However, the contributions of these tRNA mutations to the development of diseases is mostly unknown. Mutations may affect interactions with (mt)tRNA maturation enzymes or protein synthesis machinery leading to mitochondrial dysfunction. In human mitochondria, in most cases the first step of tRNA processing is the removal of the 5′ leader of precursor tRNAs (pre-tRNA) catalyzed by the three-component enzyme, mtRNase P. Additionally, one component of mtRNase P, mitochondrial RNase P protein 1 (MRPP1), catalyzes methylation of the R9 base in pre-tRNAs. Despite the central role of 5′ end processing in mitochondrial tRNA maturation, the link between mtRNase P and diseases is mostly unexplored. Here, we investigate how 11 different human disease-linked mutations in (mt)pre-tRNAIle, (mt)pre-tRNALeu(UUR), and (mt)pre-tRNAMet affect the activities of mtRNase P. We find that several mutations weaken the pre-tRNA binding affinity (KDs are approximately two- to sixfold higher than that of wild-type), while the majority of mutations decrease 5′ end processing and methylation activity catalyzed by mtRNase P (up to ∼55% and 90% reduction, respectively). Furthermore, all of the investigated mutations in (mt)pre-tRNALeu(UUR) alter the tRNA fold which contributes to the partial loss of function of mtRNase P. Overall, these results reveal an etiological link between early steps of (mt)tRNA-substrate processing and mitochondrial disease.

Predicting the pathogenicity of novel variants in mitochondrial tRNA with MitoTIP.

Novel or rare variants in mitochondrial tRNA sequences may be observed after mitochondrial DNA analysis. Determining whether these variants are pathogenic is critical, but confirmation of the effect of a variant on mitochondrial function can be challenging. We have used available databases of benign and pathogenic variants, alignment between diverse tRNAs, structural information and comparative genomics to predict the impact of all possible single-base variants and deletions. The Mitochondrial tRNA Informatics Predictor (MitoTIP) is available through MITOMAP at www.mitomap.org. The source code for MitoTIP is available at www.github.com/sonneysa/MitoTIP.

Site-selected incorporation of 5-carboxymethylaminomethyl(-2-thio)uridine into RNA sequences by phosphoramidite chemistry

5-Carboxymethylaminomethyluridine (cmnm(5)U) and 5-carboxymethylaminomethyl-2-thiouridine (cmnm(5)s(2)U) are located at the wobble position in several cytosolic and mitochondrial tRNA sequences. In this paper, we report the first site-selected incorporation of cmnm(5)U and cmnm(5)s(2)U into RNA sequences by phosphoramidite chemistry on a CPG solid support. Trifluoroacetyl and 2-(trimethylsilyl)ethyl were selected for the protection of the amine and carboxyl functions, respectively.

Comparison of Boolean analysis and standard phylogenetic methods using artificially evolved and natural mt-tRNA sequences from great apes

Boolean analysis (or BOOL-AN; Jakó et al., 2009. BOOL-AN: A method for comparative sequence analysis and phylogenetic reconstruction. Mol. Phylogenet. Evol. 52, 887–97.), a recently developed method for sequence comparison uses the Iterative Canonical Form of Boolean functions. It considers sequence information in a way entirely different from standard phylogenetic methods (i.e. Maximum Parsimony, Maximum-Likelihood, Neighbor-Joining, and Bayesian analysis). The performance and reliability of Boolean analysis were tested and compared with the standard phylogenetic methods, using artificially evolved – simulated – nucleotide sequences and the 22 mitochondrial tRNA genes of the great apes. At the outset, we assumed that the phylogeny of Hominidae is generally well established, and the guide tree of artificial sequence evolution can also be used as a benchmark. These offer a possibility to compare and test the performance of different phylogenetic methods. Trees were reconstructed by each method from 2500 simulated sequences and 22 mitochondrial tRNA sequences. We also introduced a special re-sampling method for Boolean analysis on permuted sequence sites, the P-BOOL-AN procedure. Considering the reliability values (branch support values of consensus trees and Robinson–Foulds distances) we used for simulated sequence trees produced by different phylogenetic methods, BOOL-AN appeared as the most reliable method. Although the mitochondrial tRNA sequences of great apes are relatively short (59–75 bases long) and the ratio of their constant characters is about 75%, BOOL-AN, P-BOOL-AN and the Bayesian approach produced the same tree-topology as the established phylogeny, while the outcomes of Maximum Parsimony, Maximum-Likelihood and Neighbor-Joining methods were equivocal. We conclude that Boolean analysis is a promising alternative to existing methods of sequence comparison for phylogenetic reconstruction and congruence analysis.

Improved systematic tRNA gene annotation allows new insights into the evolution of mitochondrial tRNA structures and into the mechanisms of mitochondrial genome rearrangements

Transfer RNAs (tRNAs) are present in all types of cells as well as in organelles. tRNAs of animal mitochondria show a low level of primary sequence conservation and exhibit ‘bizarre’ secondary structures, lacking complete domains of the common cloverleaf. Such sequences are hard to detect and hence frequently missed in computational analyses and mitochondrial genome annotation. Here, we introduce an automatic annotation procedure for mitochondrial tRNA genes in Metazoa based on sequence and structural information in manually curated covariance models. The method, applied to re-annotate 1876 available metazoan mitochondrial RefSeq genomes, allows to distinguish between remaining functional genes and degrading ‘pseudogenes’, even at early stages of divergence. The subsequent analysis of a comprehensive set of mitochondrial tRNA genes gives new insights into the evolution of structures of mitochondrial tRNA sequences as well as into the mechanisms of genome rearrangements. We find frequent losses of tRNA genes concentrated in basal Metazoa, frequent independent losses of individual parts of tRNA genes, particularly in Arthropoda, and wide-spread conserved overlaps of tRNAs in opposite reading direction. Direct evidence for several recent Tandem Duplication-Random Loss events is gained, demonstrating that this mechanism has an impact on the appearance of new mitochondrial gene orders.

Mitochondrial genome sequence of Unionicola parkeri (Acari: Trombidiformes: Unionicolidae): molecular synapomorphies between closely-related Unionicola gill mites

The mitochondrial genome of Unionicola parkeri is a 14,734 bp circular DNA molecule. The sequence and annotation revealed a unique gene order, related to but distinct from the gene order in the closely related species U. foili. Mitochondrial tRNA sequences annotated in this genome predict non-canonical secondary structures for these molecules. The continuing patterns of unique gene orders and unusual tRNA structures in the Trombidiformes in general and Unionicola in particular support the use of phylogenetic approaches that use these types of molecular features as shared, derived character states. Further progress in using these molecular character states to reconstruct phylogeny will depend on careful annotation, especially of tRNA genes.

Mitochondrial genome sequence of Unionicola foili (Acari: Unionicolidae): a unique gene order with implications for phylogenetic inference

The mitochondrial genome of Unionicola foili is circular, 14,738 bp in length, and contains several notable features. The sequence and annotation revealed a unique gene order, continuing a pattern of highly rearranged mitochondrial genomes in the Trombidiformes. U. foili mitochondrial tRNA sequences predict non-canonical secondary structures for these molecules, and our annotation suggests an in-frame fusion between the nad4L and nad5 genes in this genome. The unique gene order and unusual tRNA structures could serve as idiosyncratic characters and have the potential to be phylogenetically informative.

A manually curated database of tetrapod mitochondrially encoded tRNA sequences and secondary structures

BackgroundMitochondrial tRNAs have been the subject of study for structural biologists interested in their secondary structure characteristics, evolutionary biologists have researched patterns of compensatory and structural evolution and medical studies have been directed towards understanding the basis of human disease. However, an up to date, manually curated database of mitochondrially encoded tRNAs from higher animals is currently not available.DescriptionWe obtained the complete mitochondrial sequence for 277 tetrapod species from GenBank and re-annotated all of the tRNAs based on a multiple alignment of each tRNA gene and secondary structure prediction made independently for each tRNA. The mitochondrial (mt) tRNA sequences and the secondary structure based multiple alignments are freely available as Supplemental Information online.ConclusionWe compiled a manually curated database of mitochondrially encoded tRNAs from tetrapods with completely sequenced genomes. In the course of our work, we reannotated more than 10% of all tetrapod mt-tRNAs and subsequently predicted the secondary structures of 6060 mitochondrial tRNAs. This carefully constructed database can be utilized to enhance our knowledge in several different fields including the evolution of mt-tRNA secondary structure and prediction of pathogenic mt-tRNA mutations. In addition, researchers reporting novel mitochondrial genome sequences should check their tRNA gene annotations against our database to ensure a higher level of fidelity of their annotation.

Mitochondrial tRNA sequences as unusual replication origins: Pathogenic implications for Homo sapiens

The heavy strand of vertebrate mitochondrial genomes accumulates deaminations proportionally to the time it spends single-stranded during replication. A previous study showed that the strength of genome-wide deamination gradients originating from tRNA gene's locations increases with their capacities to form secondary structures resembling mitochondrial origins of light strand replication (OL), suggesting an alternative function for tRNA sequences. We hypothesize that this function is frequently pathogenic for those tRNA genes that normally do not form OL-like structures, because this could cause excess mutations in genome regions unadapted to tolerate them. In human mitochondrial genomes, pathogenic tRNA variants usually form less OL-like structures than non-pathogenic ones in cases where the normal non-pathogenic tRNA variant can function as OL, as evolutionary analyses reveal. For tRNAs lacking the putative OL-like functioning capacity, pathogenic variants form more OL-like secondary structures, particularly structures that might invoke bi-directional replication (true for 14 among 21 tRNA species, p<0.05, sign test; significantly at p<0.05 (1 tailed test) for 7 tRNA species), but not more unidirectional replication invoking structures. Accounting for the functional cloverleaf-like structure-forming capacities of tRNAs yields similar results. Rare, non-pathogenic tRNA mutants tend to form more OL-like structures than the common, non-pathogenic ones, suggesting weak directional selection also among non-pathogenic variants. The duration spent single stranded by a region of the heavy strand (DssH) during replication, estimated by integrating over all regions that can function as OL in Homo sapiens mitochondrial genomes, increases with distance of that region from the Dloop. This suggests convergence of single-strandedness during replication and transcription, and explains conserved locations of tRNA species in mitochondrial genomes and bacterial operons. These locations minimize deamination costs only in anticodons and not in other tRNA regions, during replication and transcription. Therefore, putative functioning as OLs by tRNA sequences is normal at some locations and pathogenic at others.

Aminoacylation and conformational properties of yeast mitochondrial tRNA mutants with respiratory deficiency

We report the identification and characterization of eight yeast mitochondrial tRNA mutants, located in mitochondrial tRNA(Gln), tRNA(Arg2), tRNA(Ile), tRNA(His), and tRNA(Cys), the respiratory phenotypes of which exhibit various degrees of deficiency. The mutations consist in single-base substitutions, insertions, or deletions, and are distributed all over the tRNA sequence and structure. To identify the features responsible for the defective phenotypes, we analyzed the effect of the different mutations on the electrophoretic mobility and efficiency of acylation of the mutated tRNAs in comparison with the respective wild-type molecules. Five of the studied mutations determine both conformational changes and defective acylation, while two have neither or limited effect. However, variations in structure and acylation are not necessarily correlated; the remaining mutation affects the tRNA conformation, but not its acylation properties. Analysis of tRNA structures and of mitochondrial and cytoplasmic yeast tRNA sequences allowed us to propose explanations for the observed defects, which can be ascribed to either the loss of identity nucleotides or, more often, of specific secondary and/or tertiary interactions that are largely conserved in native mitochondrial and cytoplasmic tRNAs.

Ancestral sequence reconstruction in primate mitochondrial DNA: compositional bias and effect on functional inference.

Reconstruction of ancestral DNA and amino acid sequences is an important means of inferring information about past evolutionary events. Such reconstructions suggest changes in molecular function and evolutionary processes over the course of evolution and are used to infer adaptation and convergence. Maximum likelihood (ML) is generally thought to provide relatively accurate reconstructed sequences compared to parsimony, but both methods lead to the inference of multiple directional changes in nucleotide frequencies in primate mitochondrial DNA (mtDNA). To better understand this surprising result, as well as to better understand how parsimony and ML differ, we constructed a series of computationally simple "conditional pathway" methods that differed in the number of substitutions allowed per site along each branch, and we also evaluated the entire Bayesian posterior frequency distribution of reconstructed ancestral states. We analyzed primate mitochondrial cytochrome b (Cyt-b) and cytochrome oxidase subunit I (COI) genes and found that ML reconstructs ancestral frequencies that are often more different from tip sequences than are parsimony reconstructions. In contrast, frequency reconstructions based on the posterior ensemble more closely resemble extant nucleotide frequencies. Simulations indicate that these differences in ancestral sequence inference are probably due to deterministic bias caused by high uncertainty in the optimization-based ancestral reconstruction methods (parsimony, ML, Bayesian maximum a posteriori). In contrast, ancestral nucleotide frequencies based on an average of the Bayesian set of credible ancestral sequences are much less biased. The methods involving simpler conditional pathway calculations have slightly reduced likelihood values compared to full likelihood calculations, but they can provide fairly unbiased nucleotide reconstructions and may be useful in more complex phylogenetic analyses than considered here due to their speed and flexibility. To determine whether biased reconstructions using optimization methods might affect inferences of functional properties, ancestral primate mitochondrial tRNA sequences were inferred and helix-forming propensities for conserved pairs were evaluated in silico. For ambiguously reconstructed nucleotides at sites with high base composition variability, ancestral tRNA sequences from Bayesian analyses were more compatible with canonical base pairing than were those inferred by other methods. Thus, nucleotide bias in reconstructed sequences apparently can lead to serious bias and inaccuracies in functional predictions.

The Mauriceville retroplasmid reverse transcriptase initiates cDNA synthesis de novo at the 3' end of tRNAs.

The Mauriceville retroplasmid of Neurospora mitochondria encodes a novel reverse transcriptase that initiates cDNA synthesis de novo (i.e., without a primer) at the 3' CCA of the plasmid transcript's 3' tRNA-like structure (H. Wang and A. M. Lambowitz, Cell 75:1071-1081, 1993). Here, we show that the plasmid reverse transcriptase also initiates cDNA synthesis de novo at the 3' end of tRNAs, leading to synthesis of a full-length cDNA copy of the tRNA. The use of tRNA templates in vivo was suggested previously by the structure of suppressive mutant plasmids that have incorporated mitochondrial tRNA sequences (R. A. Akins, R. L. Kelley, and A. M. Lambowitz, Cell 47:505-516, 1986). The in vitro experiments show that efficient de novo initiation on tRNA templates requires an unpaired 3' CCA and occurs predominantly opposite position C-2 of the 3' CCA sequence, the same position as in the plasmid transcript. In other reactions, the plasmid reverse transcriptase synthesizes cDNA dimers by template switching between two tRNA templates and initiates at an internal position in a tRNA by using the 3' end of the tRNA as a primer. Finally, we show that template switching between the tRNA and the plasmid transcript in vitro gives rise to hybrid cDNAs of the type predicted to be intermediates in the generation of the suppressive mutant plasmids. The ability of the plasmid reverse transcriptase to initiate at the 3' end of tRNAs presumably reflects the recognition of structural features similar to those of the 3' tRNA-like structure of the plasmid transcript. The recognition of tRNAs or tRNA-like structures as templates for cDNA synthesis may be characteristic of primitive reverse transcriptases that evolved from RNA-dependent RNA polymerases.

Evolutionary analysis of sea urchin mitochondrial tRNAs: folding of the molecules as suggested by the non-random occurrence of nucleotides

Comparative analyses of the mitochondrial tRNA sequences of the sea urchins Arbacia lixula, Paracentrotus lividus and Strongylocentrotus purpuratus revealed that conserved nucleotides may be involved in determining the typical L-shaped spatial conformation of tRNAs. These results shed light on the specific tertiary interactions that allow the folding of the atypical mitochondrial tRNAs into a functional form. A consensus mitochondrial tRNA secondary structure was derived. It shows the presence of nucleotides virtually conserved only in these organisms that represent a sort of molecular signature in sea urchins and suggests a possible physiological role. Finally, we speculate that the non-canonical structure of animal tRNAs, as well as the deviations from the universality of the genetic code, may be due to the reduction in size of the metazoan mitochondrial genome, with the concomitant acquisition of new functions by the mitochondrial tRNAs.

The presence of codon-anticodon pairs in the acceptor stem of tRNAs.

A total of 1268 available (excluding mitochondrial) tRNA sequences was used to reconstruct the common consensus image of their acceptor domains. Its structure appeared as a 11-bp-long double-stranded palindrome with complementary triplets in the center, each flanked by the 3'-ACCD and NGGU-5' motifs on each strand (D, base determinator). The palindrome readily extends up to the modern tRNA-like cloverleaf passing through an intermediate hairpin having in the center the single-stranded triplet, in supplement to its double-stranded precursor. The latter might represent an original anticodon-codon pair mapped at 1-2-3 positions of the present-day tRNA acceptors. This conclusion is supported by the striking correlation: in pairs of consensus tRNAs with complementary anticodons, their bases at the 2nd position of the acceptor stem were also complementary. Accordingly, inverse complementarity was also evident at the 71st position of the acceptor stem. With a single exception (tRNA(Phe)-tRNA(Glu) pair), the parallelism is especially impressive for the pairs of tRNAs recognized by aminoacyl-tRNA synthetases (aaRS) from the opposite classes. The above complementarity still doubly presented at the key central position of real single-stranded anticodons and their hypothetical double-stranded precursors is consistent with our previous data pointing to the double-strand use of ancient RNAs in the origin of the main actors in translation- tRNAs with complementary anticodons and the two classes of aaRS.

Permutation of a pair of tertiary nucleotides in a transfer RNA.

The tertiary nucleotides at positions 15 and 48 in a tRNA establish non-Watson-Crick hydrogen interactions that connect the dihydrouridine (D) loop with the variable loop and stabilize the "L"-shaped tRNA structure. Although the majority of tRNAs have G15.C48 or A15.U48, all of the 16 possible nucleotide pairs at positions 15 and 48 can be found in the existing cytoplasmic and mitochondrial tRNA sequences. Because tRNAs contain a variety of slightly different sets of tertiary nucleotides, this complexity raises the question of whether a given tRNA sequence framework can accommodate all of the 16 compositions at positions 15 and 48. In this work, G15 and C48 in an Escherichia coli alanine amber suppressor tRNA were permuted, and variants were tested for biological activity in vivo. All but an A15.A48 variant were functional, indicating substantial flexibility at positions 15 and 48 to accommodate nucleotide variations. Analysis of the A15.A48 variant with chemical probes showed that this mutant harbors a defect that specifically changes the conformation of the anticodon sequence. Interestingly, human tRNA(Ala) has A15.A48. Additional nucleotide substitutions in E. coli A15.A48 tRNA(Ala) that recreate the D loop sequence of human tRNA(Ala) restored the biological activity to this tRNA by reestablishing the wild-type conformation of the anticodon sequence. The results suggest a distal relationship between the D and the anticodon loops in a tRNA and delineate covariation of specific nucleotides in the evolution of tRNA(Ala) from E. coli to human.

Rainbow trout mitochondrial DNA: sequence and structural characteristics of the non-coding control region and flanking tRNA genes

We have determined the sequence of the 1003-bp control region (also referred to as the ‘displacement-loop region’) and flanking tRNA genes in the mitochondrial DNA (mtDNA) of the rainbow trout, Oncorhynchus mykiss. This region has the same overall structure (i.e., trnT- trnP-control region- trnF) as in mammalian and amphibian ( Xenopus laevis) mitochondrial (mt)DNAs. The trout control region contains apparent homologues of the conserved sequence blocks (CSB) and termination-associated sequences identified in all other vertebrate mtDNA control regions; however, it is distinguished by having an imperfect direct repeat (68/73 bp; 77% positional identity) in the right domain (proximal to the phenylalanine tRNA gene), downstream from CSB-3 in the direction of heavy-strand transcription. Within the control region, rainbow trout mtDNA shares considerable sequence similarity with the mtDNAs of Atlantic cod, Gadus morhua (Johansen et al., Nucleic Acids Res. 18 (1990) 411–419] and white sturgeon, Acipenser transmontanus (Buroker et al., Genetics 124 (1990) 157-163). The highest level of identity in pairwise comparisons is 60–70% over about 80 bp in the right domain (encompassing CSB-2 and CSB-3). About 270 bp comprising the central domain of the control region (encompassing a polypyrimidine tract) are more moderately conserved (55–60% identity in pairwise comparisons), while the left domain is highly divergent. Comparison of five trout mitochondrial tRNA sequences with their human and X. laevis homologues emphasizes the strong A+T substitution bias shown by the human sequences.

In vitro processing of transcripts containing novel tRNA-like sequences (?t-elements?) encoded by wheat mitochondrial DNA

We have recently described the properties of a wheat mitochondrial extract that is able to process, accurately and efficiently, artificial transcripts containing wheat mitochondrial tRNA sequences, with the production of mature tRNAs (P.J. Hanic-Joyce and M.W. Gray, J. Biol. Chem., in press). Such processing involves 5'-endonucleolytic, 3'-endonucleolytic, and tRNA nucleotidyltransferase activities. Here we show that this system also acts on transcripts containing sequences corresponding to an unusual class of short repeats ('t-elements') in wheat mtDNA. These repeats are theoretically capable of assuming a tRNA-like secondary structure, although stable transcripts corresponding to them are not detectable in vivo. We find that t-element sequences are processed with the same specificity and with comparable efficiency as are authentic tRNA sequences. Because known t-elements are located close to and in the same transcriptional orientation as active genes (18S-5S, 26S, tRNA(Pro)) in wheat mtDNA, our results raise the question of whether t-elements play a role in gene expression in wheat mitochondria.

Processing of transfer RNA precursors in a wheat mitochondrial extract.

To investigate mechanisms for processing of plant mitochondrial RNAs, we studied the fate of wheat mitochondrial tRNA precursors in a homologous soluble extract. Artificial precursor transcripts were synthesized in vitro using T3 or T7 RNA polymerase and DNA templates containing wheat mitochondrial tRNA genes and flanking sequences. We found that the mitochondrial extract supports processing of precursors containing both native and chloroplast-like (Joyce, P. B. M., and Gray, M. W. (1989) Nucleic Acids Res. 77, 5461-5476) wheat mitochondrial tRNA sequences. Incubation of precursor transcripts with the extract results in processing of tRNAs via precise 5'- and 3'-endonucleolytic cleavages. However, these cleavages are not ordered in vitro because intermediates composed of 5'-leader + tRNA and tRNA + 3'-trailer are present simultaneously throughout the course of the reaction. Sequence analysis of processed products confirmed that endonucleolytic cleavages occur at the expected positions, generating tRNAs with 5'-phosphoryl and 3'-hydroxyl termini. The mitochondrial extract also contains a tRNA nucleotidyltransferase activity that adds -CCAOH termini to the 3'-ends of processed tRNAs. This cell-free RNA processing system provides the basis for biochemical characterization of the various enzymes involved in the production and maturation of plant mitochondrial tRNAs.

Genes for tRNAAsp, tRNAPro, tRNATyr and two tRNAsSer in wheat mitochondrial DNA

We have begun a systematic search for potential tRNA genes in wheat mtDNA, and present here the sequences of regions of the wheat mitochondrial genome that encode genes for tRNA(Asp) (anticodon GUC), tRNA(Pro) (UGG), tRNA(Tyr) (GUA), and two tRNAs(Ser) (UGA and GCU). These genes are all solitary, not immediately adjacent to other tRNA or known protein coding genes. Each of the encoded tRNAs can assume a secondary structure that conforms to the standard cloverleaf model, and that displays none of the structural aberrations peculiar to some of the corresponding mitochondrial tRNAs from other eukaryotes. The wheat mitochondrial tRNA sequences are, on average, substantially more similar to their eubacterial and chloroplast counterparts than to their homologues in fungal and animal mitochondria. However, an analysis of regions ∼ 150 nucleotides upstream and ∼ 100 nucleotides downstream of the tRNA coding regions has revealed no obvious conserved sequences that resemble the promoter and terminator motifs that regulate the expression of eubacterial and some chloroplast tRNA genes. When restriction digests of wheat mtDNA are probed with (32)P-labelled wheat mitochondrial tRNAs, <20 hybridizing bands are detected, whether enzymes with 4 bp or 6 bp recognition sites are used. This suggests that the wheat mitochondrial genome, despite its large size, may carry a relatively small number of tRNA genes.