Gammaherpesvirus infection triggers the formation of tRNA fragments from premature tRNAs.

Transfer RNAs (tRNAs) are fundamental for both cellular and viral gene expression during viral infection. In addition, mounting evidence supports biological function for tRNA cleavage products, including in the control of gene expression during conditions of stress and infection. We previously reported that infection with the model murine gammaherpesvirus, MHV68, leads to enhanced tRNA transcription. However, whether this has any influence on tRNA transcript processing, viral replication, or the host response is not known. Here, we combined two new approaches, sequencing library preparation by Ordered Two Template Relay (OTTR) and tRNA bioinformatic analysis by tRAX, to quantitatively profile full-length tRNAs and tRNA fragment (tRF) identities during MHV68 infection. We find that MHV68 infection triggers both pre-tRNA and mature tRNA cleavage, resulting in the accumulation of specific tRFs. OTTR-tRAX revealed not only host tRNAome changes, but also the expression patterns of virally-encoded tRNAs (virtRNAs) and virtRFs made from the MHV68 genome, including their base modification signatures. Because the transcript ends of several host tRFs matched tRNA splice junctions, we tested and confirmed the role of tRNA splicing factors TSEN2 and CLP1 in MHV68-induced tRF biogenesis. Further, we show that CLP1 kinase, and by extension tRNA splicing, is required for productive MHV68 infection. Our findings provide new insight into how gammaherpesvirus infection both impacts and relies on tRNA transcription and processing. Importance: Diverse conditions of infection and cellular stress incite the cleavage of transfer RNAs, leading to the formation of tRNA fragments which can directly regulate gene expression. In our study of gammaherpesviruses, such as the murine herpesvirus 68 and human oncogenic Kaposi Sarcoma associated Herpesvirus, we discovered that transfer RNA regulation and cleavage is a key component of gene reprogramming during infection. We present the first in-depth profile of tRNA fragment generation in response to DNA virus infection, using state-of-the-art sequencing techniques that overcome several challenges with tRNA sequencing. We present several lines of evidence that tRNA fragments are made from newly-transcribed premature tRNAs and propose that this may be a defining characteristic of tRNA cleavage in some contexts. Finally, we show that tRNA splicing machinery is involved with the formation of some MHV68-induced tRNA fragments, with a key regulator of splicing, CLP1, required for maximal viral titer. Together, we posit that tRNA processing may be integral to the elegant shift in gene expression that occurs during viral take-over of the host cell.

Mistranslating the genetic code with leucine in yeast and mammalian cells.

Translation fidelity relies on accurate aminoacylation of transfer RNAs (tRNAs) by aminoacyl-tRNA synthetases (AARSs). AARSs specific for alanine (Ala), leucine (Leu), serine, and pyrrolysine do not recognize the anticodon bases. Single nucleotide anticodon variants in their cognate tRNAs can lead to mistranslation. Human genomes include both rare and more common mistranslating tRNA variants. We investigated three rare human tRNALeu variants that mis-incorporate Leu at phenylalanine or tryptophan codons. Expression of each tRNALeu anticodon variant in neuroblastoma cells caused defects in fluorescent protein production without significantly increased cytotoxicity under normal conditions or in the context of proteasome inhibition. Using tRNA sequencing and mass spectrometry we confirmed that each tRNALeu variant was expressed and generated mistranslation with Leu. To probe the flexibility of the entire genetic code towards Leu mis-incorporation, we created 64 yeast strains to express all possible tRNALeu anticodon variants in a doxycycline-inducible system. While some variants showed mild or no growth defects, many anticodon variants, enriched with G/C at positions 35 and 36, including those replacing Leu for proline, arginine, alanine, or glycine, caused dramatic reductions in growth. Differential phenotypic defects were observed for tRNALeu mutants with synonymous anticodons and for different tRNALeu isoacceptors with the same anticodon. A comparison to tRNAAla anticodon variants demonstrates that Ala mis-incorporation is more tolerable than Leu at nearly every codon. The data show that the nature of the amino acid substitution, the tRNA gene, and the anticodon are each important factors that influence the ability of cells to tolerate mistranslating tRNAs.

Gammaherpesvirus infection alters transfer RNA splicing and triggers tRNA cleavage.

Transfer RNAs (tRNAs) are fundamental for both cellular and viral gene expression during viral infection. Moreover, mounting evidence supports a noncanonical role for tRNA cleavage products in the control of gene expression during diverse conditions of stress and infection. We previously reported that infection with the model murine gammaherpesvirus, MHV68, leads to altered tRNA transcription, suggesting that tRNA regulation may play an important role in mediating viral replication or the host response. To better understand how viral infection alters tRNA expression, we combined Ordered Two Template Relay (OTTR) with tRNA-specific bioinformatic software called tRAX to profile full-length tRNAs and fragmented tRNA-derived RNAs (tDRs) during infection with MHV68. We find that OTTR-tRAX is a powerful sequencing strategy for combined tRNA/tDR profiling and reveals that MHV68 infection triggers pre-tRNA and mature tRNA cleavage, resulting in the accumulation of specific tDRs. Fragments of virally-encoded tRNAs (virtRNAs), as well as virtRNA base modification signatures are also detectable during infection. We present evidence that tRNA splicing factors are involved in the biogenesis of MHV68-induced cleavage products from pre-tRNAs and, in the case of CLP1 kinase, impact infectious virus production. Our data offers new insights into the importance of tRNA processing during gammaherpesvirus infection.

tRNATyr has an unusually short half-life in Trypanosoma brucei.

Following transcription, tRNAs undergo a series of processing and modification events to become functional adaptors in protein synthesis. Eukaryotes have also evolved intracellular transport systems whereby nucleus-encoded tRNAs may travel out and into the nucleus. In trypanosomes, nearly all tRNAs are also imported from the cytoplasm into the mitochondrion, which lacks tRNA genes. Differential subcellular localization of the cytoplasmic splicing machinery and a nuclear enzyme responsible for queuosine modification at the anticodon "wobble" position appear to be important quality control mechanisms for tRNATyr, the only intron-containing tRNA in T. brucei Since tRNA-guanine transglycosylase (TGT), the enzyme responsible for Q formation, cannot act on an intron-containing tRNA, retrograde nuclear transport is an essential step in maturation. Unlike maturation/processing pathways, the general mechanisms of tRNA stabilization and degradation in T. brucei are poorly understood. Using a combination of cellular and molecular approaches, we show that tRNATyr has an unusually short half-life. tRNATyr, and in addition tRNAAsp, also show the presence of slow-migrating bands during electrophoresis; we term these conformers: alt-tRNATyr and alt-tRNAAsp, respectively. Although we do not know the chemical or structural nature of these conformers, alt-tRNATyr has a short half-life resembling that of tRNATyr; the same is not true for alt-tRNAAsp We also show that RRP44, which is usually an exosome subunit in other organisms, is involved in tRNA degradation of the only intron-containing tRNA in T. brucei and is partly responsible for its unusually short half-life.

tRNA methylation resolves codon usage bias at the limit of cell viability.

Codon usage of each genome is closely correlated with the abundance of tRNA isoacceptors. How codon usage bias is resolved by tRNA post-transcriptional modifications is largely unknown. Here we demonstrate that the N1-methylation of guanosine at position 37 (m1G37) on the 3'-side of the anticodon, while not directly responsible for reading of codons, is a neutralizer that resolves differential decoding of proline codons. A genome-wide suppressor screen of a non-viable Escherichia coli strain, lacking m1G37, identifies proS suppressor mutations, indicating a coupling of methylation with tRNA prolyl-aminoacylation that sets the limit of cell viability. Using these suppressors, where prolyl-aminoacylation is decoupled from tRNA methylation, we show that m1G37 neutralizes differential translation of proline codons by the major isoacceptor. Lack of m1G37 inactivates this neutralization and exposes the need for a minor isoacceptor for cell viability. This work has medical implications for bacterial species that exclusively use the major isoacceptor for survival.

Systematic identification of tRNA genes in Drosophila melanogaster.

Transfer RNAs (tRNAs) are ubiquitous adapter molecules that link specific codons in messenger RNA (mRNA) with their corresponding amino acids during protein synthesis. The tRNA genes of Drosophila have been investigated for over half a century but have lacked systematic identification and nomenclature. Here, we review and integrate data within FlyBase and the Genomic tRNA Database (GtRNAdb) to identify the full complement of tRNA genes in the D. melanogaster nuclear and mitochondrial genomes. We apply a logical and informative nomenclature to all tRNA genes, and provide an overview of their characteristics and genomic features.

Eukaryotic tRNA sequences present conserved and amino acid-specific structural signatures.

Metazoan organisms have many tRNA genes responsible for decoding amino acids. The set of all tRNA genes can be grouped in sets of common amino acids and isoacceptor tRNAs that are aminoacylated by corresponding aminoacyl-tRNA synthetases. Analysis of tRNA alignments shows that, despite the high number of tRNA genes, specific tRNA sequence motifs are highly conserved across multicellular eukaryotes. The conservation often extends throughout the isoacceptors and isodecoders with, in some cases, two sets of conserved isodecoders. This study is focused on non-Watson-Crick base pairs in the helical stems, especially GoU pairs. Each of the four helical stems may contain one or more conserved GoU pairs. Some are amino acid specific and could represent identity elements for the cognate aminoacyl tRNA synthetases. Other GoU pairs are found in more than a single amino acid and could be critical for native folding of the tRNAs. Interestingly, some GoU pairs are anticodon-specific, and others are found in phylogenetically-specific clades. Although the distribution of conservation likely reflects a balance between accommodating isotype-specific functions as well as those shared by all tRNAs essential for ribosomal translation, such conservations may indicate the existence of specialized tRNAs for specific translation targets, cellular conditions, or alternative functions.

Distinct Stress-Dependent Signatures of Cellular and Extracellular tRNA-Derived Small RNAs.

The cellular response to stress is an important determinant of disease pathogenesis. Uncovering the molecular fingerprints of distinct stress responses may identify novel biomarkers and key signaling pathways for different diseases. Emerging evidence shows that transfer RNA-derived small RNAs (tDRs) play pivotal roles in stress responses. However, RNA modifications present on tDRs are barriers to accurately quantifying tDRs using traditional small RNA sequencing. Here, AlkB-facilitated methylation sequencing is used to generate a comprehensive landscape of cellular and extracellular tDR abundances in various cell types during different stress responses. Extracellular tDRs are found to have distinct fragmentation signatures from intracellular tDRs and these tDR signatures are better indicators of different stress responses than miRNAs. These distinct extracellular tDR fragmentation patterns and signatures are also observed in plasma from patients on cardiopulmonary bypass. It is additionally demonstrated that angiogenin and RNASE1 are themselves regulated by stressors and contribute to the stress-modulated abundance of sub-populations of cellular and extracellular tDRs. Finally, a sub-population of extracellular tDRs is identified for which AGO2 appears to be required for their expression. Together, these findings provide a detailed profile of stress-responsive tDRs and provide insight about tDR biogenesis and stability in response to cellular stressors.

tRNAscan-SE 2.0: improved detection and functional classification of transfer RNA genes.

tRNAscan-SE has been widely used for transfer RNA (tRNA) gene prediction for over twenty years, developed just as the first genomes were decoded. With the massive increase in quantity and phylogenetic diversity of genomes, the accurate detection and functional prediction of tRNAs has become more challenging. Utilizing a vastly larger training set, we created nearly one hundred specialized isotype- and clade-specific models, greatly improving tRNAscan-SE's ability to identify and classify both typical and atypical tRNAs. We employ a new comparative multi-model strategy where predicted tRNAs are scored against a full set of isotype-specific covariance models, allowing functional prediction based on both the anticodon and the highest-scoring isotype model. Comparative model scoring has also enhanced the program's ability to detect tRNA-derived SINEs and other likely pseudogenes. For the first time, tRNAscan-SE also includes fast and highly accurate detection of mitochondrial tRNAs using newly developed models. Overall, tRNA detection sensitivity and specificity is improved for all isotypes, particularly those utilizing specialized models for selenocysteine and the three subtypes of tRNA genes encoding a CAU anticodon. These enhancements will provide researchers with more accurate and detailed tRNA annotation for a wider variety of tRNAs, and may direct attention to tRNAs with novel traits.

R2DT is a framework for predicting and visualising RNA secondary structure using templates.

Non-coding RNAs (ncRNA) are essential for all life, and their functions often depend on their secondary (2D) and tertiary structure. Despite the abundance of software for the visualisation of ncRNAs, few automatically generate consistent and recognisable 2D layouts, which makes it challenging for users to construct, compare and analyse structures. Here, we present R2DT, a method for predicting and visualising a wide range of RNA structures in standardised layouts. R2DT is based on a library of 3,647 templates representing the majority of known structured RNAs. R2DT has been applied to ncRNA sequences from the RNAcentral database and produced >13 million diagrams, creating the world's largest RNA 2D structure dataset. The software is amenable to community expansion, and is freely available at https://github.com/rnacentral/R2DT and a web server is found at https://rnacentral.org/r2dt .

A guide to naming human non-coding RNA genes.

Research on non-coding RNA (ncRNA) is a rapidly expanding field. Providing an official gene symbol and name to ncRNA genes brings order to otherwise potential chaos as it allows unambiguous communication about each gene. The HUGO Gene Nomenclature Committee (HGNC, www.genenames.org) is the only group with the authority to approve symbols for human genes. The HGNC works with specialist advisors for different classes of ncRNA to ensure that ncRNA nomenclature is accurate and informative, where possible. Here, we review each major class of ncRNA that is currently annotated in the human genome and describe how each class is assigned a standardised nomenclature.

Predicting transfer RNA gene activity from sequence and genome context.

Transfer RNA (tRNA) genes are among the most highly transcribed genes in the genome owing to their central role in protein synthesis. However, there is evidence for a broad range of gene expression across tRNA loci. This complexity, combined with difficulty in measuring transcript abundance and high sequence identity across transcripts, has severely limited our collective understanding of tRNA gene expression regulation and evolution. We establish sequence-based correlates to tRNA gene expression and develop a tRNA gene classification method that does not require, but benefits from, comparative genomic information and achieves accuracy comparable to molecular assays. We observe that guanine + cytosine (G + C) content and CpG density surrounding tRNA loci is exceptionally well correlated with tRNA gene activity, supporting a prominent regulatory role of the local genomic context in combination with internal sequence features. We use our tRNA gene activity predictions in conjunction with a comprehensive tRNA gene ortholog set spanning 29 placental mammals to estimate the evolutionary rate of functional changes among orthologs. Our method adds a new dimension to large-scale tRNA functional prediction and will help prioritize characterization of functional tRNA variants. Its simplicity and robustness should enable development of similar approaches for other clades, as well as exploration of functional diversification of members of large gene families.

Distinct Modified Nucleosides in tRNATrp from the Hyperthermophilic Archaeon Thermococcus kodakarensis and Requirement of tRNA m2G10/m22G10 Methyltransferase (Archaeal Trm11) for Survival at High Temperatures.

tRNA m2G10/m22G10 methyltransferase (archaeal Trm11) methylates the 2-amino group in guanosine at position 10 in tRNA and forms N2,N2-dimethylguanosine (m22G10) via N2-methylguanosine (m2G10). We determined the complete sequence of tRNATrp, one of the substrate tRNAs for archaeal Trm11 from Thermococcus kodakarensis, a hyperthermophilic archaeon. Liquid chromatography/mass spectrometry following enzymatic digestion of tRNATrp identified 15 types of modified nucleoside at 21 positions. Several modifications were found at novel positions in tRNA, including 2'-O-methylcytidine at position 6, 2-thiocytidine at position 17, 2'-O-methyluridine at position 20, 5,2'-O-dimethylcytidine at position 32, and 2'-O-methylguanosine at position 42. Furthermore, methylwyosine was found at position 37 in this tRNATrp, although 1-methylguanosine is generally found at this location in tRNATrp from other archaea. We constructed trm11 (Δtrm11) and some gene disruptant strains and compared their tRNATrp with that of the wild-type strain, which confirmed the absence of m22G10 and other corresponding modifications, respectively. The lack of 2-methylguanosine (m2G) at position 67 in the trm11 trm14 double disruptant strain suggested that this methylation is mediated by Trm14, which was previously identified as an m2G6 methyltransferase. The Δtrm11 strain grew poorly at 95°C, indicating that archaeal Trm11 is required for T. kodakarensis survival at high temperatures. The m22G10 modification might have effects on stabilization of tRNA and/or correct folding of tRNA at the high temperatures. Collectively, these results provide new clues to the function of modifications and the substrate specificities of modification enzymes in archaeal tRNA, enabling us to propose a strategy for tRNA stabilization of this archaeon at high temperatures.IMPORTANCEThermococcus kodakarensis is a hyperthermophilic archaeon that can grow at 60 to 100°C. The sequence of tRNATrp from this archaeon was determined by liquid chromatography/mass spectrometry. Fifteen types of modified nucleoside were observed at 21 positions, including 5 modifications at novel positions; in addition, methylwyosine at position 37 was newly observed in an archaeal tRNATrp The construction of trm11 (Δtrm11) and other gene disruptant strains confirmed the enzymes responsible for modifications in this tRNA. The lack of 2-methylguanosine (m2G) at position 67 in the trm11 trm14 double disruptant strain suggested that this position is methylated by Trm14, which was previously identified as an m2G6 methyltransferase. The Δtrm11 strain grew poorly at 95°C, indicating that archaeal Trm11 is required for T. kodakarensis survival at high temperatures.

tRNAviz: explore and visualize tRNA sequence features.

Transfer RNAs (tRNAs) are ubiquitous across the tree of life. Although tRNA structure is highly conserved, there is still significant variation in sequence features between clades, isotypes and even isodecoders. This variation not only impacts translation, but as shown by a variety of recent studies, nontranslation-associated functions are also sensitive to small changes in tRNA sequence. Despite the rapidly growing number of sequenced genomes, there is a lack of tools for both small- and large-scale comparative genomics analysis of tRNA sequence features. Here, we have integrated over 150 000 tRNAs spanning all domains of life into tRNAviz, a web application for exploring and visualizing tRNA sequence features. tRNAviz implements a framework for determining consensus sequence features and can generate sequence feature distributions by isotypes, clades and anticodons, among other tRNA properties such as score. All visualizations are interactive and exportable. The web server is publicly available at http://trna.ucsc.edu/tRNAviz/.

tRNAscan-SE: Searching for tRNA Genes in Genomic Sequences.

Transfer RNAs are the largest, most complex non-coding RNA family, universal to all living organisms. tRNAscan-SE has been the de facto tool for predicting tRNA genes in whole genomes. The newly developed version 2.0 has incorporated advanced methodologies with improved probabilistic search software and a suite of new gene models, enabling better functional classification of predicted genes. This chapter describes the use of the UNIX command-driven and online web versions, illustrating different search modes and options.

Matching tRNA modifications in humans to their known and predicted enzymes.

tRNA are post-transcriptionally modified by chemical modifications that affect all aspects of tRNA biology. An increasing number of mutations underlying human genetic diseases map to genes encoding for tRNA modification enzymes. However, our knowledge on human tRNA-modification genes remains fragmentary and the most comprehensive RNA modification database currently contains information on approximately 20% of human cytosolic tRNAs, primarily based on biochemical studies. Recent high-throughput methods such as DM-tRNA-seq now allow annotation of a majority of tRNAs for six specific base modifications. Furthermore, we identified large gaps in knowledge when we predicted all cytosolic and mitochondrial human tRNA modification genes. Only 48% of the candidate cytosolic tRNA modification enzymes have been experimentally validated in mammals (either directly or in a heterologous system). Approximately 23% of the modification genes (cytosolic and mitochondrial combined) remain unknown. We discuss these 'unidentified enzymes' cases in detail and propose candidates whenever possible. Finally, tissue-specific expression analysis shows that modification genes are highly expressed in proliferative tissues like testis and transformed cells, but scarcely in differentiated tissues, with the exception of the cerebellum. Our work provides a comprehensive up to date compilation of human tRNA modifications and their enzymes that can be used as a resource for further studies.

Transfer RNA genes experience exceptionally elevated mutation rates.

Transfer RNAs (tRNAs) are a central component for the biological synthesis of proteins, and they are among the most highly conserved and frequently transcribed genes in all living things. Despite their clear significance for fundamental cellular processes, the forces governing tRNA evolution are poorly understood. We present evidence that transcription-associated mutagenesis and strong purifying selection are key determinants of patterns of sequence variation within and surrounding tRNA genes in humans and diverse model organisms. Remarkably, the mutation rate at broadly expressed cytosolic tRNA loci is likely between 7 and 10 times greater than the nuclear genome average. Furthermore, evolutionary analyses provide strong evidence that tRNA genes, but not their flanking sequences, experience strong purifying selection acting against this elevated mutation rate. We also find a strong correlation between tRNA expression levels and the mutation rates in their immediate flanking regions, suggesting a simple method for estimating individual tRNA gene activity. Collectively, this study illuminates the extreme competing forces in tRNA gene evolution and indicates that mutations at tRNA loci contribute disproportionately to mutational load and have unexplored fitness consequences in human populations.

Methylation guide RNA evolution in archaea: structure, function and genomic organization of 110 C/D box sRNA families across six Pyrobaculum species.

Archaeal homologs of eukaryotic C/D box small nucleolar RNAs (C/D box sRNAs) guide precise 2'-O-methyl modification of ribosomal and transfer RNAs. Although C/D box sRNA genes constitute one of the largest RNA gene families in archaeal thermophiles, most genomes have incomplete sRNA gene annotation because reliable, fully automated detection methods are not available. We expanded and curated a comprehensive gene set across six species of the crenarchaeal genus Pyrobaculum, particularly rich in C/D box sRNA genes. Using high-throughput small RNA sequencing, specialized computational searches and comparative genomics, we analyzed 526 Pyrobaculum C/D box sRNAs, organizing them into 110 families based on synteny and conservation of guide sequences which determine methylation targets. We examined gene duplications and rearrangements, including one family that has expanded in a pattern similar to retrotransposed repetitive elements in eukaryotes. New training data and inclusion of kink-turn secondary structural features enabled creation of an improved search model. Our analyses provide the most comprehensive, dynamic view of C/D box sRNA evolutionary history within a genus, in terms of modification function, feature plasticity, and gene mobility.