Search

Pseudourine (Ψ) is the most abundant cellular RNA modification, present in tRNA, rRNA, snRNA, mRNA, long noncoding RNA (lncRNA), and others. Ψ sites and fractions are dynamically regulated in stress response and across development stages. Although high throughput Ψ sequencing methods based on N-Cyclohexyl-N'-(2-morpholinoethyl)carbodiimide (CMC) reaction are available for Ψ detection transcriptome-wide, a simple method for the analysis of specific, targeted Ψ sites and their fraction quantitation is needed to better investigate Ψ function. Here, we describe an RT-PCR and gel electrophoresis based method that can sensitively and quantitatively assess Ψ at single-nucleotide resolution in mRNA/lncRNA, termed CMC-RT and ligation assisted PCR analysis of Ψ modification (CLAP). The principle of the CMC-method is the reverse transcription stop induced by the CMC-Ψ adduct. In CLAP, CMC reaction is first carried out with the RNA sample. Reverse transcription using a non-processive RT produces two cDNA products for each RNA transcript, one with the 3' end at the Ψ site, the other read-through product from the unmodified RNA. Using splint oligonucleotide assisted site-specific ligation, these two cDNA products are then visualized on a gel as two distinct PCR products in the same lane corresponding to the Ψ-modified and unmodified target site. CLAP validates Ψ sites identified by high throughput sequencing, quantifies Ψ levels in mRNA and lncRNA, and enables convenient and rapid investigation on the function and mechanism of the Ψ modification.

Epitranscriptomic mechanisms linking tRNA function and the brain proteome to cognition and complex behaviors are not well described. Here, we report bi-directional changes in depression-related behaviors after genetic disruption of neuronal tRNA cytosine methylation, including conditional ablation and transgene-derived overexpression of Nsun2 in the mouse prefrontal cortex (PFC). Neuronal Nsun2-deficiency was associated with a decrease in tRNA m5C levels, resulting in deficits in expression of 70% of tRNAGly isodecoders. Altogether, 1488/5820 proteins changed upon neuronal Nsun2-deficiency, in conjunction with glycine codon-specific defects in translational efficiencies. Loss of Gly-rich proteins critical for glutamatergic neurotransmission was associated with impaired synaptic signaling at PFC pyramidal neurons and defective contextual fear memory. Changes in the neuronal translatome were also associated with a 146% increase in glycine biosynthesis. These findings highlight the methylation sensitivity of glycinergic tRNAs in the adult PFC. Furthermore, they link synaptic plasticity and complex behaviors to epitranscriptomic modifications of cognate tRNAs and the proteomic homeostasis associated with specific amino acids., The link between tRNA modifications, protein translation, and behavior is unclear. Here, the authors show that neuronal Nsun2 deficiency results in codon-specific epitranscriptomic changes of Gly-tRNAs and proteomic changes affecting synaptic signaling and behavior in mice.

Transfer RNAs (tRNAs) are essential for translation, and tRNA expression and modifications are regulated by many factors. However, the interplay between the microbiome and host tRNA profiles through host-microbiome interactions has not been explored. In this study, we investigated host-microbiome interactions via the tRNA profiling of four tissue types from germ-free and specific pathogen-free mice. Our analyses reveal that cytosolic and mitochondrial tRNA expression and tRNA modifications in the host are reprogrammed in a tissue-specific and microbiome-dependent manner. In terms of tRNA expression, the intestines and brains are more sensitive to the influence of the microbiome than the livers and kidneys. In terms of tRNA modifications, cytosolic tRNAs show more obvious changes in the livers and kidneys in the presence of the microbiome. Our findings reveal a previously unexplored relationship among the microbiome, tRNA abundance, and epitranscriptome in a mammalian host.

The accuracy in pairing tRNAs with correct amino acids by aminoacyl-tRNA synthetases (aaRSs) dictates the fidelity of translation. To ensure fidelity, multiple aaRSs developed editing functions that remove a wrong amino acid from tRNA before it reaches the ribosome. However, no specific mechanism within an aaRS is known to handle the scenario where a cognate amino acid is mischarged onto a wrong tRNA, as exemplified by AlaRS mischarging alanine to G4:U69-containing tRNAThr. Here, we report that the mischargeable G4:U69-containing tRNAThr are strictly conserved in vertebrates and are ubiquitously and abundantly expressed in mammalian cells and tissues. Although these tRNAs are efficiently mischarged, no corresponding Thr-to-Ala mistranslation is detectable. Mistranslation is prevented by a robust proofreading activity of ThrRS towards Ala-tRNAThr. Therefore, while wrong amino acids are corrected within an aaRS, a wrong tRNA is handled in trans by an aaRS cognate to the mischarged tRNA species. Interestingly, although Ala-tRNAThr mischarging is not known to occur in bacteria, Escherichia coli ThrRS also possesses robust cross-editing ability. We propose that the cross-editing activity of ThrRS is evolutionarily conserved and that this intrinsic activity allows G4:U69-containing tRNAThr to emerge and be preserved in vertebrates to have alternative functions without compromising translational fidelity.

Pseudouridine (Ψ) is the most abundant RNA modification in cellular RNA present in tRNA/rRNA/snRNA and also in mRNA and long noncoding RNA (lncRNA). Elucidation of Ψ function in mRNA/lncRNA requires mapping and quantitative assessment of its modification fraction at single-base resolution. The most widely used Ψ mapping method for mRNA/lncRNA relies on its reaction with N-Cyclohexyl-N′-(2-morpholinoethyl)carbodiimide (CMC), forming an adduct with the Ψ base in RNA that is detectable by reverse transcription (RT) stops. However, this method has not produced consistent Ψ maps in mRNAs; furthermore, available protocols do not lend confidence to the estimation of Ψ fraction at specific sites, which is a crucial parameter for investigating the biological relevance of mRNA modifications. Here we develop a quantitative RT-PCR based method that can detect and quantify the modification fraction of target Ψ sites in mRNA/lncRNA, termed CMC-RT and ligation assisted PCR analysis of Ψ modification (CLAP). The method still relies on RT stop at a CMC-Ψ site, but uses site-specific ligation and PCR to generate two distinct PCR products in the same sample, corresponding to the modified and unmodified site, that are visualized by gel electrophoresis. CLAP not only requires a small amount of cellular RNA to validate Ψ sites but also determines the Ψ fraction semiquantitatively at target sites in mRNA/lncRNA. We determined the Ψ status of four mRNA sites and one lncRNA site whose modification fractions range from 30% to 84% in three human cell lines. Our method enables precise mapping and assessment of Ψ modification levels in low abundance cellular RNAs.

Significance Charcot-Marie-Tooth disease (CMT) is a devastating motor and sensory neuropathy with an estimated 100,000 afflicted individuals in the US. Unexpectedly, aminoacyl-tRNA synthetases are the largest disease-associated protein family. A natural explanation is that the disease is associated with weak translation or mistranslation (caused by editing defects). However, our results with six different disease-causing mutants in AlaRS ruled out defects in aminoacylation or editing as causal factors. Instead, specific mutant proteins gained a neuropilin 1 (Nrp1)-AlaRS interaction. Previously a gain of Nrp1 interaction with a different disease-causing tRNA synthetase was mechanistically linked to the pathology of CMT. Thus, our results raise the possibility that pathological engagement of Nrp1 is common to at least a subset of tRNA synthetase-associated cases of CMT., Through dominant mutations, aminoacyl-tRNA synthetases constitute the largest protein family linked to Charcot-Marie-Tooth disease (CMT). An example is CMT subtype 2N (CMT2N), caused by individual mutations spread out in AlaRS, including three in the aminoacylation domain, thereby suggesting a role for a tRNA-charging defect. However, here we found that two are aminoacylation defective but that the most widely distributed R329H is normal as a purified protein in vitro and in unfractionated patient cell samples. Remarkably, in contrast to wild-type (WT) AlaRS, all three mutant proteins gained the ability to interact with neuropilin 1 (Nrp1), the receptor previously linked to CMT pathogenesis in GlyRS. The aberrant AlaRS-Nrp1 interaction is further confirmed in patient samples carrying the R329H mutation. However, CMT2N mutations outside the aminoacylation domain do not induce the Nrp1 interaction. Detailed biochemical and biophysical investigations, including X-ray crystallography, small-angle X-ray scattering, hydrogen-deuterium exchange (HDX), switchSENSE hydrodynamic diameter determinations, and protease digestions reveal a mutation-induced structural loosening of the aminoacylation domain that correlates with the Nrp1 interaction. The b1b2 domains of Nrp1 are responsible for the interaction with R329H AlaRS. The results suggest Nrp1 is more broadly associated with CMT-associated members of the tRNA synthetase family. Moreover, we revealed a distinct structural loosening effect induced by a mutation in the editing domain and a lack of conformational impact with C-Ala domain mutations, indicating mutations in the same protein may cause neuropathy through different mechanisms. Our results show that, as with other CMT-associated tRNA synthetases, aminoacylation per se is not relevant to the pathology.

AlkB is a DNA/RNA repair enzyme that removes base alkylations such as N1-methyladenosine (m1A) or N3-methylcytosine (m3C) from DNA and RNA. The AlkB enzyme has been used as a critical tool to facilitate tRNA sequencing and identification of mRNA modifications. As a tool, AlkB mutants with better reactivity and new functionalities are highly desired; however, previous identification of such AlkB mutants was based on the classical approach of targeted mutagenesis. Here, we introduce a high-throughput screening method to evaluate libraries of AlkB variants for demethylation activity on RNA and DNA substrates. This method is based on a fluorogenic RNA aptamer with an internal modified RNA/DNA residue which can block reverse transcription or introduce mutations leading to loss of fluorescence inherent in the cDNA product. Demethylation by an AlkB variant eliminates the blockage or mutation thereby restores the fluorescence signals. We applied our screening method to sites D135 and R210 in the Escherichia coli AlkB protein and identified a variant with improved activity beyond a previously known hyperactive mutant toward N1-methylguanosine (m1G) in RNA. We also applied our method to O6-methylguanosine (O6mG) modified DNA substrates and identified candidate AlkB variants with demethylating activity. Our study provides a high-throughput screening method for in vitro evolution of any demethylase enzyme.

Members of the mammalian AlkB family are known to mediate nucleic acid demethylation1,2. ALKBH7, a mammalian AlkB homologue, localizes in mitochondria and affects metabolism3, but its function and mechanism of action are unknown. Here we report an approach to site-specifically detect N1-methyladenosine (m1A), N3-methylcytidine (m3C), N1-methylguanosine (m1G) and N2,N2-dimethylguanosine (m22G) modifications simultaneously within all cellular RNAs, and discovered that human ALKBH7 demethylates m22G and m1A within mitochondrial Ile and Leu1 pre-tRNA regions, respectively, in nascent polycistronic mitochondrial RNA4–6. We further show that ALKBH7 regulates the processing and structural dynamics of polycistronic mitochondrial RNAs. Depletion of ALKBH7 leads to increased polycistronic mitochondrial RNA processing, reduced steady-state mitochondria-encoded tRNA levels and protein translation, and notably decreased mitochondrial activity. Thus, we identify ALKBH7 as an RNA demethylase that controls nascent mitochondrial RNA processing and mitochondrial activity. Zhang et al. identify ALKBH7 as the demethylase of mitochondrial pre-tRNAs that regulates nascent mitochondrial RNA processing and translation. ALKBH7 loss impairs mitochondrial functions including fatty acid oxidation, leading to obesity.

Queuosine (Q) is a conserved tRNA modification in bacteria and eukaryotes. Eukaryotic Q-tRNA modification occurs through replacing the guanine base with the scavenged metabolite queuine at the wobble position of tRNAs with G34U35N36 anticodon (Tyr, His, Asn, Asp) by the QTRT1/QTRT2 heterodimeric enzyme encoded in the genome. In humans, Q-modification in tRNATyr and tRNAAsp are further glycosylated with galactose and mannose, respectively. Although galactosyl-Q (galQ) and mannosyl-Q (manQ) can be measured by LC/MS approaches, the difficulty of detecting and quantifying these modifications with low sample inputs has hindered their biological investigations. Here we describe a simple acid denaturing gel and nonradioactive northern blot method to detect and quantify the fraction of galQ/manQ-modified tRNA using just microgram amounts of total RNA. Our method relies on the secondary amine group of galQ/manQ becoming positively charged to slow their migration in acid denaturing gels commonly used for tRNA charging studies. We apply this method to determine the Q and galQ/manQ modification kinetics in three human cells lines. For Q-modification, tRNAAsp is modified the fastest, followed by tRNAHis, tRNATyr, and tRNAAsn. Compared to Q-modification, glycosylation occurs at a much slower rate for tRNAAsp, but at a similar rate for tRNATyr. Our method enables easy access to study the function of these enigmatic tRNA modifications.

Transfer RNA (tRNA) is present at tens of millions of transcripts in a human cell and is the most abundant RNA in moles among all cellular RNAs. tRNA is also the most extensively modified RNA with, on an average, 13 modifications per molecule. The primary function of tRNA as the adaptor of amino acids and the genetic code in protein synthesis is well known. tRNA modifications play multi-faceted roles in decoding and other cellular processes. The abundance, modification, and aminoacylation (charging) levels of tRNAs contribute to mRNA decoding in ways that reflect the cell type and its environment; however, how these factors work together to maximize translation efficiency remains to be understood. tRNAs also interact with many proteins not involved in translation and this may coordinate translation activity and other processes in the cell. This review focuses on the modifications and the functional genomics of human tRNA and discusses future perspectives on the explorations of human tRNA biology.

Introduction : tRNA-derived fragments (tRFs) have been postulated to play a significant role in human disease and specifically cancer. tRFs can function as small non-coding RNAs, expressed more highly in conditions of oxidative stress and tissue damage, and are postulated to have a regulatory function. Previous work has identified overexpression of tRNAs in multiple myeloma (MM), but the expression of tRFs in MM has not yet been characterized. In this descriptive study, we set out to apply modern next-generation sequencing tools to describe tRNA and tRF expression in MM. Methods: Bone marrow aspirate samples from 24 patients, 8 with monoclonal gammopathy of undetermined significance (MGUS) or smoldering MM (SMM), 8 with newly diagnosed MM (NDMM), and 8 with relapsed/refractory MM (RRMM) were selected. Total RNA fraction from mononuclear cell fractions was extracted and used to prepare tRNA sequencing libraries which generated data of different small RNA families such as tRNA, tRF, and 5S rRNA in one library. Results: We found a trend toward increased total expression of tRFs in NDMM and RRMM patient samples. Analysis of specific classes of tRFs revealed differential expression along the course of disease progression, with specific tRFs increasing in expression from MGUS/SMM to RRMM. Specifically, tRFs from tRNAHis showed a significant increase in expression from MGUS/SMM to NDMM (p Conclusions : These preliminary results indicate that tRNA fragment expression has the potential to serve as an objective, quantitative measure of clinical MM progression. While the function of these small RNA species remains unclear, specific cleavage patterns suggest that tRF expression in MM may be regulated during disease progression. Future work should explore potential for miRNA-like function of tRFs in MM, and ongoing work to be presented at the meeting will further explore tRF abundance as well as post-transcriptional modifications of tRNAs which may provide a mechanistic understanding of patterns of tRNA fragmentation. Citation Format: Alexandra E. Rojek, Christopher D. Katanski, Andrew Stefka, Benjamin A. Derman, Andrzej Jakubowiak, Tao Pan. tRNA expression and tRFs in multiple myeloma: Progression from monoclonal gammopathies to relapsed/refractory disease [abstract]. In: Proceedings of the American Association for Cancer Research Annual Meeting 2021; 2021 Apr 10-15 and May 17-21. Philadelphia (PA): AACR; Cancer Res 2021;81(13_Suppl):Abstract nr 2394.

The abundant Watson-Crick face methylations in biological RNAs such as N1-methyladenosine (m1A), N1-methylguanosine (m1G), N3-methylcytosine (m3C), and N2,N2-dimethylguanosine (m22G) cause significant obstacles for high throughput RNA sequencing by impairing cDNA synthesis. One strategy to overcome this obstacle is to remove the methyl group on these modified bases prior to cDNA synthesis using enzymes. The wild-type E. coli AlkB and its D135S mutant can remove most of m1A, m1G, m3C modifications in transfer RNA (tRNA), but they work poorly on m22G. Here we report the design and evaluation of a series of AlkB mutants against m22G-containing model RNA substrates that we synthesize using an improved chemical synthesis method. We show that the AlkB D135S/L118V mutant efficiently and selectively converts m22G modification to N2-methylguanosine (m2G). We also show that this new enzyme lowers the amount of m22G in tRNA and improves the efficiency of tRNA sequencing.

Global protein mistranslation with methionine has been shown to be a conserved biological process that affords distinct functional advantages in all three domains of life. In all instances, methionine mistranslation occurs through a regulated process where low-fidelity forms of methionyl-tRNA synthetase are conditionally induced to mischarge non-methionyl-tRNAs with methionine followed by the utilization of the misacylated tRNAs in translation. In mammals, methionine mistranslation contributes to oxidative stress response; in the hyperthermophilic archaeon Aeropyrum pernix, methionine mistranslation produces proteins that are better adapted to low temperature growth; in E. coli, methionine mistranslation increases resistance to antibiotics and chemical stressors. The phenotypic benefits conferred by tRNA mismethionylation suggest that it should be a widespread adaptational mechanism in diverse bacterial lineages, yet this response has only been described in E. coli. Furthermore, previous microbial investigations on this response have been confined to axenic laboratory cultures. It was unknown whether tRNA mismethionylation was relevant in a natural microbial habitat. Here we show that four abundant gut microbiotal genera belonging to the Firmicutes and Bacteroidetes phyla perform constitutive tRNA misacylation with methionine in the mouse cecum in situ. These results reveal the ubiquity of the tRNA mismethionylation process among bacteria and implicate the potential importance of this response for subsistence and adaptation in natural habitats.

The fidelity of tRNA aminoacylation is a critical determinant for the ultimate accuracy of protein synthesis. Although aminoacyl-tRNA synthetases are assumed to consistently maintain high tRNA charging fidelity, recent evidence has demonstrated that the fidelity of the aminoacylation reaction can be actively regulated and liable to change. Accordingly, the ability to conveniently assay the fidelity of tRNA charging is becoming increasingly relevant for studying mistranslation. Here we describe a combined radioactivity and microarray based method that can quantitatively elucidate which individual cognate or noncognate tRNA isoacceptors are charged with amino acid. In this technique, in vitro tRNA charging reactions or in vivo pulse-labeling is performed using a radiolabeled amino acid and tRNA microarrays are used to distinguish tRNA isoacceptors in total tRNA. During the tRNA array hybridization, each tRNA will hybridize to its unique probe and subsequent phosphorimaging of the array can determine which tRNAs were aminoacylated with the radiolabeled amino acid. The method can be used to assess the fidelity of tRNA charging in vivo or in vitro and can be applied to any organism with annotated tRNA genes.

Aminoacyl-transfer RNA (tRNA) synthetases (aaRSs) are the largest protein family causatively linked to neurodegenerative Charcot–Marie–Tooth (CMT) disease. Dominant mutations cause the disease, and studies of CMT disease-causing mutant glycyl-tRNA synthetase (GlyRS) and tyrosyl-tRNA synthetase (TyrRS) showed their mutations create neomorphic structures consistent with a gain-of-function mechanism. In contrast, based on a haploid yeast model, loss of aminoacylation function was reported for CMT disease mutants in histidyl-tRNA synthetase (HisRS). However, neither that nor prior work of any CMT disease-causing aaRS investigated the aminoacylation status of tRNAs in the cellular milieu of actual patients. Using an assay that interrogated aminoacylation levels in patient cells, we investigated a HisRS-linked CMT disease family with the most severe disease phenotype. Strikingly, no difference in charged tRNA levels between normal and diseased family members was found. In confirmation, recombinant versions of 4 other HisRS CMT disease-causing mutants showed no correlation between activity loss in vitro and severity of phenotype in vivo. Indeed, a mutation having the most detrimental impact on activity was associated with a mild disease phenotype. In further work, using 3 independent biophysical analyses, structural opening (relaxation) of mutant HisRSs at the dimer interface best correlated with disease severity. In fact, the HisRS mutation in the severely afflicted patient family caused the largest degree of structural relaxation. These data suggest that HisRS-linked CMT disease arises from open conformation-induced mechanisms distinct from loss of aminoacylation.

Queuosine (Q) is a hypermodified base in the wobble anticodon position of tRNAs coding for the amino acids Tyr, His, Asn, and Asp. tRNA Q-modification is introduced by a queuine tRNA-ribosyltransferase (TGT) that replaces the guanine base at G34 at these tRNAs with the modified base. tRNA Q-modification is widely distributed among prokaryotic and eukaryotic organisms, but only bacteria synthesize Q-modified tRNA de novo. In mammals, tRNA Q-modifications strictly rely on the presence of gut microbiomes or diets to produce the queuine base. Despite decades of study, cellular roles of tRNA Q-modification are still not fully understood. Here we describe a method to quantify tRNA Q-modification levels in individual tRNAs from human cells based on the presence of a cis-diol in the Q modification. This cis-diol moiety slows modified tRNA migration through polyacrylamide gels supplemented with N-acryloyl-3-aminophenylboronic acid (APB) compared to the unmodified tRNA. This difference can be visualized by Northern blots using probes for specific tRNA.

Fidelity of translation, which is predominately dictated by the accuracy of aminoacyl-tRNA synthetases in pairing amino acids with correct tRNAs, is of central importance in biology. Yet, deliberate modifications of translational fidelity can be beneficial. Here we found human and not E. coli AlaRS has an intrinsic capacity for mis-pairing alanine onto non-alanyl-tRNAs including tRNACys. Consistently, a cysteine-to-alanine substitution was found in a reporter protein expressed in human cells. All human AlaRS-mischarged tRNAs have a G4:U69 base pair in the acceptor stem. The base pair is required for the mischarging. By solving the crystal structure of human AlaRS and comparing it to that of E. coli AlaRS, we identified a key sequence divergence between eukaryotes and bacteria that influences mischarging. Thus, the expanded tRNA specificity of AlaRS appears to be an evolutionary gain-of-function to provide posttranscriptional alanine substitutions in eukaryotic proteins for potential regulations.

Eukaryotic transfer RNAs contain on average 14 modifications. Investigations of their biological functions require the determination of the modification sites and the dynamic variations of the modification fraction. Base methylation represents a major class of tRNA modification. Although many approaches have been used to identify tRNA base methylations, including sequencing, they are generally qualitative and do not report the information on the modification fraction. Dynamic mRNA modifications have been shown to play important biological roles; yet, the extent of tRNA modification fractions has not been reported systemically. Here we take advantage of a recently developed high-throughput sequencing method (DM-tRNA-seq) to identify and quantify tRNA base methylations located at the Watson-Crick face in HEK293T cells at single base resolution. We apply information derived from both base mutations and positional stops from sequencing using a combination of demethylase treatment and cDNA synthesis by a thermophilic reverse transcriptase to compile a quantitative “Modification Index” (MI) for six base methylations in human tRNA and rRNA. MI combines the metrics for mutational and stop components from alignment of sequencing data without demethylase treatment, and the modifications are validated in the sequencing data upon demethylase treatment. We identify many new methylation sites in both human nuclear and mitochondrial-encoded tRNAs not present in the RNA modification databases. The potentially quantitative nature of the MI values obtained from sequencing is validated by primer extension of several tRNAs. Our approach should be widely applicable to identify tRNA methylation sites, analyze comparative fractional modifications, and evaluate the modification dynamics between different samples.

It is common wisdom that codon usage bias has evolved in the selection for efficient translation, in which highly expressed genes are encoded predominantly by optimal codons. However, a growing body of evidence suggests regulatory roles for non-optimal codons in translation dynamics. Here we report that in mammalian cells, non-optimal codons play a critical role in promoting selective mRNA translation during amino acid starvation. During starvation, in contrast to genes encoding ribosomal proteins whose translation is highly sensitive to amino acid deprivation, translation of genes involved in the cellular protein degradation pathways remains unaffected. We found that these two gene groups bear different codon composition, with non-optimal codons being highly enriched in genes encoding the ubiquitin–proteasome system. Supporting the selective tRNA charging model originally proposed in Escherichia coli, we demonstrated that tRNA isoacceptors decoding rare codons are maintained in translating ribosomes under amino acid starvation. Finally, using luciferase reporters fused with endogenous gene-derived sequences, we show that codon optimality contributes to differential mRNA translation in response to amino acid starvation. These results highlight the physiological significance of codon usage bias in cellular adaptation to stress.

Translation of genetic information into functional proteins is critical for all cellular life. Accurate protein synthesis relies on proper aminoacylation of transfer RNAs (tRNAs) and decoding of mRNAs by the ribosome with the use of aminoacyl-tRNAs. Mistranslation can lead to pathologic consequences. All cells contain elaborate quality control mechanisms in translation, although translational fidelity may be regulated by various factors such as nutrient limitation or reactive oxygen species. Translation fidelity is maintained via the accuracy of tRNA aminoacylation by the aminoacyl-tRNA synthetases and matching of the mRNA codon with the tRNA anticodon by the ribosome. Stringent substrate discrimination and proofreading are critical in aminoacylating tRNAs with their cognate amino acid to maintain high accuracy of translation. Although the composition of the cellular proteome generally adheres to the genetic code, accumulating evidence indicates that cells can also deliberately mistranslate; they synthesize mutant proteins that deviate from the genetic code in response to stress or environmental changes. Mistranslation with tRNA charged with noncognate amino acids can expand the proteome to enhance stress response and help adaptation. Here, we review current knowledge on mistranslation through tRNA misacylation and describe advances in our understanding of translational control in the regulation of stress response and human diseases.

In mammalian cells under oxidative stress, the methionyl-tRNA synthetase (MetRS) misacylates noncognate tRNAs at frequencies as high as 10% distributed among up to 28 tRNA species. Instead of being detrimental for the cell, misincorporation of methionine residues in the proteome reduces the risk of oxidative damage to proteins, which aids the oxidative stress response. tRNA microarrays have been essential for the detection of the full pattern of misacylated tRNAs, but have limited capacity to investigate the misacylation and mistranslation mechanisms in live cells. Here we develop a dual-fluorescence reporter to specifically measure methionine misincorporation at glutamic acid codons GAA and GAG via tRNAGlu mismethionylation in human cells. Our method relies on mutating a specific Met codon in the active site of the fluorescent protein mCherry to a Glu codon that renders mCherry nonfluorescent when translation follows the genetic code. Mistranslation utilizing mismethionylated tRNAGlu restores fluorescence in proportion to the amount of misacylated tRNAGlu. This cellular approach works well for both transient transfection and established stable HEK293 lines. It is rapid, straightforward, and well suited for high-throughput activity analysis under a wide range of physiological conditions. As a proof of concept, we apply this method to characterize the effect of human tRNAGlu isodecoders on mistranslation and discuss the implications of our findings.

All organisms universally encode, synthesize and utilize proteins that function optimally within a subset of growth conditions. While healthy cells are thought to maintain high translational fidelity within their natural habitats, natural environments can easily fluctuate outside the optimal functional range of genetically encoded proteins. The hyperthermophilic archaeon Aeropyrum pernix (A. pernix) can grow throughout temperature variations ranging from 70 to 100°C, although the specific factors facilitating such adaptability are unknown. Here, we show that A. pernix undergoes constitutive leucine to methionine mistranslation at low growth temperatures. Low-temperature mistranslation is facilitated by the misacylation of tRNA(Leu) with methionine by the methionyl-tRNA synthetase (MetRS). At low growth temperatures, the A. pernix MetRS undergoes a temperature dependent shift in tRNA charging fidelity, allowing the enzyme to conditionally charge tRNA(Leu) with methionine. We demonstrate enhanced low-temperature activity for A. pernix citrate synthase that is synthesized during leucine to methionine mistranslation at low-temperature growth compared to its high-fidelity counterpart synthesized at high-temperature. Our results show that conditional leucine to methionine mistranslation can make protein adjustments capable of improving the low-temperature activity of hyperthermophilic proteins, likely by facilitating the increasing flexibility required for greater protein function at lower physiological temperatures.

Background: Transfer RNA (tRNA) queuosine (Q)-modifications occur specifically in 4 cellular tRNAs at the wobble anticodon position. tRNA Q-modification in human cells depends on the gut microbiome because the microbiome product queuine is required for its installation by the enzyme Q tRNA ribosyltransferase catalytic subunit 1 (QTRT1) encoded in the human genome. Queuine is a micronutrient from diet and microbiome. Although tRNA Q-modification has been studied for a long time regarding its properties in decoding and tRNA fragment generation, how QTRT1 affects tumorigenesis and the microbiome is still poorly understood. Results: We generated single clones of QTRT1-knockout breast cancer MCF7 cells using Double Nickase Plasmid. We also established a QTRT1-knockdown breast MDA-MB-231 cell line. The impacts of QTRT1 deletion or reduction on cell proliferation and migration in vitro were evaluated using cell culture, while the regulations on tumor growth in vivo were evaluated using a xenograft BALB/c nude mouse model. We found that QTRT1 deficiency in human breast cancer cells could change the functions of regulation genes, which are critical in cell proliferation, tight junction formation, and migration in human breast cancer cells in vitro and a breast tumor mouse model in vivo. We identified that several core bacteria, such as Lachnospiraceae, Lactobacillus, and Alistipes, were markedly changed in mice post injection with breast cancer cells. The relative abundance of bacteria in tumors induced from wildtype cells was significantly higher than those of QTRT1 deficiency cells. Conclusions: Our results demonstrate that the QTRT1 gene and tRNA Q-modification altered cell proliferation, junctions, and microbiome in tumors and the intestine, thus playing a critical role in breast cancer development.

Advances in high-throughput sequencing have facilitated remarkable insights into the diversity and functioning of naturally occurring microbes; however, current sequencing strategies are insufficient to reveal physiological states of microbial communities associated with protein translation dynamics. Transfer RNAs (tRNAs) are core components of protein synthesis machinery, present in all living cells, and are phylogenetically tractable, which make them ideal targets to gain physiological insights into environmental microbes. Here we report a direct sequencing approach, tRNA-seq, and a software suite, tRNA-seq-tools, to recover sequences, abundance profiles, and post-transcriptional modifications of microbial tRNA transcripts. Our analysis of cecal samples using tRNA-seq distinguishes high-fat- and low-fat-fed mice in a comparable fashion to 16S ribosomal RNA gene amplicons, and reveals taxon- and diet-dependent variations in tRNA modifications. Our results provide taxon-specific in situ insights into the dynamics of tRNA gene expression and post-transcriptional modifications within complex environmental microbiomes., The authors present a molecular approach and a software suite for high-throughput sequencing and characterization of microbial tRNA transcripts, which provide unique physiological insights into complex environmental microbes.

Eukaryotic transfer RNAs (tRNA) contain on average 13 modifications that perform a wide range of roles in translation and in the generation of tRNA fragments that regulate gene expression. Queuosine (Q) modification occurs in the wobble anticodon position of tRNAs for amino acids His, Asn, Tyr, and Asp. In eukaryotes, Q modification is fully dependent on diet or on gut microbiome in multicellular organisms. Despite decades of study, cellular roles of Q modification remain to be fully elucidated. Here we show that in human cells, Q modification specifically protects its cognate tRNAHis and tRNAAsn against cleavage by ribonucleases. We generated cell lines that contain completely depleted or fully Q-modified tRNAs. Using these resources, we found that Q modification significantly reduces angiogenin cleavage of its cognate tRNAs in vitro. Q modification does not change the cellular abundance of the cognate full-length tRNAs, but alters the cellular content of their fragments in vivo in the absence and presence of stress. Our results provide a new biological aspect of Q modification and a mechanism of how Q modification alters small RNA pools in human cells.

Abstarct Tylototriton anhuiensis (Salamandridae, Urodela), collected from Yaoluoping Nature Reserve, was identified as a new species. The complete mitogenome sequence of T. anhuiensis is 16,259 by in length, including 13 protein-coding genes, 22 tRNA genes, 2 rRNA genes, and D-loop region. The base composition of the mitogenome was 33.6%A, 26.3% C, 14.5% G, and 25.6% T. The ND6 subunit gene and eight tRNA genes were encoded on the L-strand, the others were encoded on the H-strand.

The sequence diversity of individual human genomes has been extensively analyzed for variations and phenotypic implications for mRNA, miRNA, and long non-coding RNA genes. TRNA (tRNA) also exhibits large sequence diversity in the human genome, but tRNA gene sequence variation and potential functional implications in individual human genomes have not been investigated. Here we capitalize on the sequencing data from the 1000-genomes project to examine the diversity of tRNA genes in the human population. Previous analysis of the reference human genome indicated an unexpected large number of diverse tRNA genes beyond the necessity of translation, suggesting that some tRNA transcripts may perform non-canonical functions. We found 24 new tRNA sequences in > 1% and 76 new tRNA sequences in > 0.2% of all individuals, indicating that tRNA genes are also subject to evolutionary changes in the human population. Unexpectedly, two abundant new tRNA genes contain base-pair mismatches in the anticodon stem. We experimentally determined that these two new tRNAs have altered structures in vitro; however, one new tRNA is not aminoacylated but extremely stable in HeLa cells, suggesting that this new tRNA can be used for non-canonical function. Our results show that at the scale of human population, tRNA genes are more diverse than conventionally understood, and some new tRNAs may perform non-canonical, extra-translational functions that may be linked to human health and disease.

Despite its biological importance, tRNA has not been adequately sequenced by standard methods because of its abundant post-transcriptional modifications and stable structure, which interfere with cDNA synthesis. We achieved efficient and quantitative tRNA sequencing in HEK293T cells by using engineered demethylases to remove base methylations and a highly processive thermostable group II intron reverse transcriptase to overcome these obstacles. Our method, DM-tRNA-seq, should be applicable to investigations of tRNA in all organisms.

Certain stress conditions can induce cleavage of tRNAs around the anticodon loop via the use of the ribonuclease angiogenin. The cellular factors that regulate tRNA cleavage are not well known. In this study we used normal and eIF2α phosphorylation-deficient mouse embryonic fibroblasts and applied a microarray-based methodology to identify and compare tRNA cleavage patterns in response to hypertonic stress, oxidative stress (arsenite), and treatment with recombinant angiogenin. In all three scenarios mouse embryonic fibroblasts deficient in eIF2α phosphorylation showed a higher accumulation of tRNA fragments including those derived from initiator-tRNA(Met). We have shown that tRNA cleavage is regulated by the availability of angiogenin, its substrate (tRNA), the levels of the angiogenin inhibitor RNH1, and the rates of protein synthesis. These conclusions are supported by the following findings: (i) exogenous treatment with angiogenin or knockdown of RNH1 increased tRNA cleavage; (ii) tRNA fragment accumulation was higher during oxidative stress than hypertonic stress, in agreement with a dramatic decrease of RNH1 levels during oxidative stress; and (iii) a positive correlation was observed between angiogenin-mediated tRNA cleavage and global protein synthesis rates. Identification of the stress-specific tRNA cleavage mechanisms and patterns will provide insights into the role of tRNA fragments in signaling pathways and stress-related disorders.

Transfer RNAs (tRNAs) are central to protein synthesis and impact translational speed and fidelity by their abundance. Here we examine the extent to which viruses manipulate tRNA populations to favor translation of their own genes. We study two very different viruses: influenza A virus (IAV), a medium-sized (13 kB genome) RNA virus; and vaccinia virus (VV), a large (200 kB genome) DNA virus. We show that the total cellular tRNA population remains unchanged following viral infection, whereas the polysome-associated tRNA population changes dramatically in a virus-specific manner. The changes in polysome-associated tRNA levels reflect the codon usage of viral genes, suggesting the existence of local tRNA pools optimized for viral translation.

Although the translational function of tRNA has long been established, extra translational functions of tRNA are still being discovered. We previously developed a computational method to systematically predict new tRNA-protein complexes and experimentally validated six candidate proteins, including the mitogen-activated protein kinase kinase 2 (MEK2), that interact with tRNA in HEK293T cells. However, consequences of the interaction between tRNA and these proteins remain to be elucidated. Here we tested the consequence of the interaction between tRNA and MEK2 in pancreatic cancer cell lines. We also generated disease and drug resistance-derived MEK2 mutants (Q60P, P128Q, S154F, E207K) to evaluate the function of the tRNA-MEK2 interaction. Our results demonstrate that tRNA interacts with the wild-type and mutant MEK2 in pancreatic cancer cells; furthermore, the MEK2 inhibitor U0126 significantly reduces the tRNA-MEK2 interaction. In addition, tRNA affects the catalytic activity of the wild type and mutant MEK2 proteins in different ways. Overall, our findings demonstrate the interaction of tRNA with MEK2 in pancreatic cancer cells and suggest that tRNA may impact MEK2 activity in cancer cells.

Transfer RNA (tRNA) is a central component of protein synthesis and cell signaling network. One salient feature of tRNA is its heavily modified status, which can critically impact its function. Here we show that mammalian ALKBH1 is a tRNA demethylase. It mediates the demethylation of N1-methyladenosine (m1A) in tRNAs. The ALKBH1-catalyzed demethylation of the target tRNAs results in attenuated translation initiation and their decreased usage in protein synthesis. This process is dynamic and responds to glucose availability to affect translation. Our results uncover reversible methylation of tRNA as a new regulatory mechanism of post-transcriptional gene expression.

SummaryTransfer RNA (tRNA) gene content is a differentiating feature of genomes that contributes to the efficiency of the translational apparatus, but the principles shaping tRNA gene copy number and codon composition are poorly understood. Here, we report that the emergence of two specific tRNA modifications shaped the structure and composition of all extant genomes. Through the analysis of more than 500 genomes, we identify two kingdom-specific tRNA modifications as major contributors that separated archaeal, bacterial, and eukaryal genomes in terms of their tRNA gene composition. We show that, contrary to prior observations, genomic codon usage and tRNA gene frequencies correlate in all kingdoms if these two modifications are taken into account and that presence or absence of these modifications explains patterns of gene expression observed in previous studies. Finally, we experimentally demonstrate that human gene expression levels correlate well with genomic codon composition if these identified modifications are considered.

Protobothrops kaulbacki is a new record in China, and in this study, the complete mitochondrial genome sequence of P. kaulbacki had been determined. The length of mitogenome is 17,237 bp, including 13 protein-coding genes, 22 tRNA genes, 2 rRNA genes (12S and 16S rRNA), 1 origin of L-strand replication (OL), and 2 control regions (CRs). The maximum-likelihood (ML) tree based on the whole mitogenome shown that P. kaulbacki belongs to the genus Protobothrops.

The N1-methyl-Adenosine (m1A58) modification at the conserved nucleotide 58 in the TΨC loop is present in most eukaryotic tRNAs. In yeast, m1A58 modification is essential for viability because it is required for the stability of the initiator-tRNAMet. However, m1A58 modification is not required for the stability of several other tRNAs in yeast. This differential m1A58 response for different tRNA species raises the question of whether some tRNAs are hypomodified at A58 in normal cells, and how hypomodification at A58 may affect the stability and function of tRNA. Here, we apply a genomic approach to determine the presence of m1A58 hypomodified tRNAs in human cell lines and show how A58 hypomodification affects stability and involvement of tRNAs in translation. Our microarray-based method detects the presence of m1A58 hypomodified tRNA species on the basis of their permissiveness in primer extension. Among five human cell lines examined, approximately one-quarter of all tRNA species are hypomodified in varying amounts, and the pattern of the hypomodified tRNAs is quite similar. In all cases, no hypomodified initiator-tRNAMet is detected, consistent with the requirement of this modification in stabilizing this tRNA in human cells. siRNA knockdown of either subunit of the m1A58-methyltransferase results in a slow-growth phenotype, and a marked increase in the amount of m1A58 hypomodified tRNAs. Most m1A58 hypomodified tRNAs can associate with polysomes in varying extents. Our results show a distinct pattern for m1A58 hypomodification in human tRNAs, and are consistent with the notion that this modification fine tunes tRNA functions in different contexts.

Recent years have seen intensive progress in measuring protein translation. However, the contributions of coding sequences to the efficiency of the process remain unclear. Here, we identify a universally conserved profile of translation efficiency along mRNAs computed based on adaptation between coding sequences and the tRNA pool. In this profile, the first approximately 30-50 codons are, on average, translated with a low efficiency. Additionally, in eukaryotes, the last approximately 50 codons show the highest efficiency over the full coding sequence. The profile accurately predicts position-dependent ribosomal density along yeast genes. These data suggest that translation speed and, as a consequence, ribosomal density are encoded by coding sequences and the tRNA pool. We suggest that the slow "ramp" at the beginning of mRNAs serves as a late stage of translation initiation, forming an optimal and robust means to reduce ribosomal traffic jams, thus minimizing the cost of protein expression.

Translational fidelity, essential for protein and cell function, requires accurate transfer RNA (tRNA) aminoacylation. Purified aminoacyl-tRNA synthetases exhibit a fidelity of one error per 10,000 to 100,000 couplings. The accuracy of tRNA aminoacylation in vivo is uncertain, however, and might be considerably lower. Here we show that in mammalian cells, approximately 1% of methionine (Met) residues used in protein synthesis are aminoacylated to non-methionyl-tRNAs. Remarkably, Met-misacylation increases up to tenfold upon exposing cells to live or non-infectious viruses, toll-like receptor ligands or chemically induced oxidative stress. Met is misacylated to specific non-methionyl-tRNA families, and these Met-misacylated tRNAs are used in translation. Met-misacylation is blocked by an inhibitor of cellular oxidases, implicating reactive oxygen species (ROS) as the misacylation trigger. Among six amino acids tested, tRNA misacylation occurs exclusively with Met. As Met residues are known to protect proteins against ROS-mediated damage, we propose that Met-misacylation functions adaptively to increase Met incorporation into proteins to protect cells against oxidative stress. In demonstrating an unexpected conditional aspect of decoding mRNA, our findings illustrate the importance of considering alternative iterations of the genetic code.

In multiple myeloma (MM), malignant plasma cells produce large amounts of antibodies and have highly active protein translational machinery. It is not known whether regulation of the abundance and aminoacylation (charging) of transfer RNA (tRNA) takes place in myeloma cells to accommodate for the increased amount of protein translation. Using tRNA-specific microarrays, we demonstrate that tRNA levels are significantly elevated in MM cell lines compared to normal bone marrow cells. We furthermore show that the addition of the proteasome inhibitor, Bortezomib (Velcade™, PS-341) results in decreased charging levels of tRNAs, in particular those coding for hydrophobic amino acids. These results suggest that tRNA properties are altered in MM to accommodate for its increased need for protein translation, and that proteasome inhibition directly impacts protein synthesis in MM through effects on tRNA charging.

Over 100 chemical types of RNA modifications have been identified in thousands of sites in all three domains of life. Recent data suggest that modifications function synergistically to mediate biological function, and that cells may coordinately modulate modification levels for regulatory purposes. However, this area of RNA biology remains largely unexplored due to the lack of robust, high-throughput methods to quantify the extent of modification at specific sites. Recently, we developed a facile enzymatic ligation-based method for detection and quantitation of methylated 2'-hydroxyl groups within RNA. Here we exploit the principles of molecular recognition and nucleic acid chemistry to establish the experimental parameters for ligation-based detection and quantitation of pseudouridine (Psi) and N6-methyladenosine (m6A), two abundant modifications in eukaryotic rRNA/tRNA and mRNA, respectively. Detection of pseudouridylation at several sites in the large subunit rRNA derived from yeast demonstrates the feasibility of the approach for analysis of pseudouridylation in biological RNA samples.

We compare the diversity of chromosomal-encoded transfer RNA (tRNA) genes from 11 eukaryotes as identified by tRNAScan-SE of their respective genomes. They include the budding and fission yeast, worm, fruit fly, fugu, chicken, dog, rat, mouse, chimp and human. The number of tRNA genes are between 170 and 570 and the number of tRNA isoacceptors range from 41 to 55. Unexpectedly, the number of tRNA genes having the same anticodon but different sequences elsewhere in the tRNA body (defined here as tRNA isodecoder genes) varies significantly (10–246). tRNA isodecoder genes allow up to 274 different tRNA species to be produced from 446 genes in humans, but only up to 51 from 275 genes in the budding yeast. The fraction of tRNA isodecoder genes among all tRNA genes increases across the phylogenetic spectrum. A large number of sequence differences in human tRNA isodecoder genes occurs in the internal promoter regions for RNA polymerase III. We also describe a systematic, ligation-based method to detect and quantify tRNA isodecoder molecules in human samples, and show differential expression of three tRNA isodecoders in six human tissues. The large number of tRNA isodecoder genes in eukaryotes suggests that tRNA function may be more diverse than previously appreciated.

Natural selection favors efficient expression of encoded proteins, but the causes, mechanisms, and fitness consequences of evolved coding changes remain an area of aggressive inquiry. We report a large-scale reversal in the relative translational accuracy of codons across 12 fly species in the Drosophila/Sophophora genus. Because the reversal involves pairs of codons that are read by the same genomically encoded tRNAs, we hypothesize, and show by direct measurement, that a tRNA anticodon modification from guanosine to queuosine has coevolved with these genomic changes. Queuosine modification is present in most organisms but its function remains unclear. Modification levels vary across developmental stages in D. melanogaster, and, consistent with a causal effect, genes maximally expressed at each stage display selection for codons that are most accurate given stage-specific queuosine modification levels. In a kinetic model, the known increased affinity of queuosine-modified tRNA for ribosomes increases the accuracy of cognate codons while reducing the accuracy of near-cognate codons. Levels of queuosine modification in D. melanogaster reflect bioavailability of the precursor queuine, which eukaryotes scavenge from the tRNAs of bacteria and absorb in the gut. These results reveal a strikingly direct mechanism by which recoding of entire genomes results from changes in utilization of a nutrient.

Transfer RNA (tRNA) is produced as a precursor molecule that needs to be processed at its 3' and 5' ends. Ribonuclease P is the sole endonuclease responsible for processing the 5' end of tRNA by cleaving the precursor and leading to tRNA maturation. It was one of the first catalytic RNA molecules identified and consists of a single RNA component in all organisms and only one protein component in bacteria. It is a true multi-turnover ribozyme and one of only two ribozymes (the other being the ribosome) that are conserved in all kingdoms of life. Here we show the crystal structure at 3.85 A resolution of the RNA component of Thermotoga maritima ribonuclease P. The entire RNA catalytic component is revealed, as well as the arrangement of the two structural domains. The structure shows the general architecture of the RNA molecule, the inter- and intra-domain interactions, the location of the universally conserved regions, the regions involved in pre-tRNA recognition and the location of the active site. A model with bound tRNA is in agreement with all existing data and suggests the general basis for RNA-RNA recognition by this ribozyme.

Aminoacylated (charged) transfer RNA isoacceptors read different messenger RNA codons for the same amino acid. The concentration of an isoacceptor and its charged fraction are principal determinants of the translation rate of its codons. A recent theoretical model predicts that amino-acid starvation results in ‘selective charging' where the charging levels of some tRNA isoacceptors will be low and those of others will remain high. Here, we developed a microarray for the analysis of charged fractions of tRNAs and measured charging for all Escherichia coli tRNAs before and during leucine, threonine or arginine starvation. Before starvation, most tRNAs were fully charged. During starvation, the isoacceptors in the leucine, threonine or arginine families showed selective charging when cells were starved for their cognate amino acid, directly confirming the theoretical prediction. Codons read by isoacceptors that retain high charging can be used for efficient translation of genes that are essential during amino-acid starvation. Selective charging can explain anomalous patterns of codon usage in the genes for different families of proteins.

Anthus hodgsoni is a species of small passerine bird in the family Motacillidae, which is widely distributed. In this study, we determined the complete mitochondrial genome of A. hodgsoni. The result showed that the total length of the mitochondrial genome was 16,886 bp and contained 2 ribosomal RNA genes, 22 transfer RNA genes, 13 protein-coding genes and 1 control region. The phylogenetic tree was reconstructed using the Bayesian analysis method and divided into four genera, Anthus, Dendronanthus, Motacilla and Tmetothylacus. The A. hodgsoni which we determined was clustered into genus Anthus and received strong support.

The complete mitochondrial genome sequence of Hyla tsinlingensis was determined in this research. The length of mitogenome is 17850 bp, including 13 protein-coding genes, 22 tRNA genes, 2 rRNA genes, 1 OL and 1 control region. The phylogentic tree was reconstructed using the BI method based on concatenated nucleotide sequences of mtDNA genes (12S ribosomal small subunit gene/12S rRNA; NADH dehydrogenase subunit 1 gene/ND1, including adjacent transfer RNAs and the partial 16S ribosomal large subunit gene). The phylogenetic tree was split into two clades, Clade A and Clade B. The H. tsinlingensis which we determined clustered into Clade A.

In this paper, the complete 16,725 bp nucleotide sequence of the mitochondrial genome was determined for the Tylototriton kweichowensis (Caudata:Salamandridae). It contains 37 genes (13 protein-coding genes, 2 rRNA genes and 22 tRNA genes) and a non-coding region (D-loop). Overall base composition of the complete mitochondrial DNA is A (33.8%), G (14.4%), C (26.2%), and T (25.6%), so the percentage of A and T (59.4%) is higher than G and C (40.6%). All the genes in T. kweichowensis are distributed on the H-strand, except for the ND6 subunit gene and eight tRNA genes which are encoded on the L-strand.

Transfer RNA (tRNA) decodes mRNA codons when aminoacylated (charged) with an amino acid at its 3′ end. Charged tRNAs turn over rapidly in cells, and variations in charged tRNA fractions are known to be a useful parameter in cellular responses to stress. tRNA charging fractions can be measured for individual tRNA species using acid denaturing gels, or comparatively at the genome level using microarrays. These hybridization-based approaches cannot be used for high resolution analysis of mammalian tRNAs due to their large sequence diversity. Here we develop a high-throughput sequencing method that enables accurate determination of charged tRNA fractions at single-base resolution (Charged DM-tRNA-seq). Our method takes advantage of the recently developed DM-tRNA-seq method, but includes additional chemical steps that specifically remove the 3′A residue in uncharged tRNA. Charging fraction is obtained by counting the fraction of A-ending reads versus A+C-ending reads for each tRNA species in the same sequencing reaction. In HEK293T cells, most cytosolic tRNAs are charged at >80% levels, whereas tRNASer and tRNAThr are charged at lower levels. These low charging levels were validated using acid denaturing gels. Our method should be widely applicable for investigations of tRNA charging as a parameter in biological regulation.