ALL-tRNAseq enables robust tRNA profiling in tissue samples

The landscape of Arabidopsis tRNA aminoacylation

The function of transfer RNAs (tRNAs) depends on enzymes that cleave primary transcript ends, add a 3' CCA tail, introduce post-transcriptional base modifications, and charge (aminoacylate) mature tRNAs with the correct amino acid. Maintaining an available pool of the resulting aminoacylated tRNAs is essential for protein synthesis. High-throughput sequencing techniques have recently been developed to provide a comprehensive view of aminoacylation state in a tRNA-specific fashion. However, these methods have never been applied to plants. Here, we treated Arabidopsis thaliana RNA samples with periodate and then performed tRNA-seq to distinguish between aminoacylated and uncharged tRNAs. This approach successfully captured every tRNA isodecoder family and detected expression of additional tRNA-like transcripts. We found that estimated aminoacylation rates and CCA tail integrity were significantly higher on average for organellar (mitochondrial and plastid) tRNAs than for nuclear/cytosolic tRNAs. Reanalysis of previously published human cell line data showed a similar pattern. Base modifications result in nucleotide misincorporations and truncations during reverse transcription, which we quantified and used to test for relationships with aminoacylation levels. We also determined that the Arabidopsis tRNA-like sequences (t-elements) that are cleaved from the ends of some mitochondrial messenger RNAs have post-transcriptionally modified bases and CCA-tail addition. However, these t-elements are not aminoacylated, indicating that they are only recognized by a subset of tRNA-interacting enzymes and do not play a role in translation. Overall, this work provides a characterization of the baseline landscape of plant tRNA aminoacylation rates and demonstrates an approach for investigating environmental and genetic perturbations to plant translation machinery.

Nanopore sequencing of intact aminoacylated tRNAs

Transfer RNAs (tRNA) are decorated during biogenesis with a variety of modifications that modulate their stability, aminoacylation, and decoding potential during translation. The complex landscape of tRNA modification presents significant analysis challenges and to date no single approach enables the simultaneous measurement of important but disparate chemical properties of individual, mature tRNA molecules. We developed a new, integrated approach to analyze the sequence, modification, and aminoacylation state of tRNA molecules in a high throughput nanopore sequencing experiment, leveraging a chemical ligation that embeds the charged amino acid in an adapted tRNA molecule. During nanopore sequencing, the embedded amino acid generates unique distortions in ionic current and translocation speed, enabling application of machine learning approaches to classify charging status and amino acid identity. Specific applications of the method indicate it will be broadly useful for examining relationships and dependencies between tRNA sequence, modification, and aminoacylation.

Characterization and implementation of the MarathonRT template-switching reaction to expand the capabilities of RNA-seq

End-to-end RNA-sequencing methods that capture 5'-sequence content without cumbersome library manipulations are of great interest, particularly for analysis of long RNAs. While template-switching methods have been developed for RNA sequencing by distributive short-read RTs, such as the MMLV RTs used in SMART-Seq methods, they have not been adapted to leverage the power of ultraprocessive RTs, such as those derived from group II introns. To facilitate this transition, we dissected the individual processes that guide the enzymatic specificity and efficiency of the multistep template-switching reaction carried out by RTs, in this case, by MarathonRT. Remarkably, this is the first study of its kind, for any RT. First, we characterized the nucleotide specificity of nontemplated addition (NTA) reaction that occurs when the RT extends past the RNA 5'-terminus. We then evaluated the binding specificity of specialized template-switching oligonucleotides, optimizing their sequences and chemical properties to guide efficient template-switching reaction. Having dissected and optimized these individual steps, we then unified them into a procedure for performing RNA sequencing with MarathonRT enzymes, using a well-characterized RNA reference set. The resulting reads span a six-log range in transcript concentration and accurately represent the input RNA identities in both length and composition. We also performed RNA-seq from total human RNA and poly(A)-enriched RNA, with short- and long-read sequencing demonstrating that MarathonRT enhances the discovery of unseen RNA molecules by conventional RT. Altogether, we have generated a new pipeline for rapid, accurate sequencing of complex RNA libraries containing mixtures of long RNA transcripts.

tRNA expression and modification landscapes, and their dynamics during zebrafish embryo development

tRNA genes exist in multiple copies in the genome of all organisms across the three domains of life. Besides the sequence differences across tRNA copies, extensive post-transcriptional modification adds a further layer to tRNA diversification. Whilst the crucial role of tRNAs as adapter molecules in protein translation is well established, whether all tRNAs are actually expressed, and whether the differences across isodecoders play any regulatory role is only recently being uncovered. Here we built upon recent developments in the use of NGS-based methods for RNA modification detection and developed tRAM-seq, an experimental protocol and in silico analysis pipeline to investigate tRNA expression and modification. Using tRAM-seq, we analysed the full ensemble of nucleo-cytoplasmic and mitochondrial tRNAs during embryonic development of the model vertebrate zebrafish. We show that the repertoire of tRNAs changes during development, with an apparent major switch in tRNA isodecoder expression and modification profile taking place around the start of gastrulation. Taken together, our findings suggest the existence of a general reprogramming of the expressed tRNA pool, possibly gearing the translational machinery for distinct stages of the delicate and crucial process of embryo development.

A user guide to RT-based mapping of RNA modifications

Chemical modifications to RNA nucleotides are both a naturally occurring layer of biological regulation and an increasingly prevalent approach to synthetically alter RNA function in therapeutic applications. Detection of their presence, prevalence, and stoichiometry across different RNAs is critical to understanding their underlying functions. However, this remains challenging due to the technical barriers involved in differentiating chemically similar modification species, and in detecting rare or low stoichiometry modifications. Reverse transcription-based techniques rely on the introduction of a predictable mutation, truncation, or deletion signature when a reverse transcriptase encounters a modified nucleotide of interest. Previous studies have shown promise in detecting modifications to single nucleotide resolution, but the low efficiency and processivity of many commercially available reverse transcriptases has resulted in discordant conclusions in some cases. Here, we present guidelines and best practices for applying the highly processive MarathonRT enzyme to reverse transcription-based modification sequencing. These guidelines include recommendations for controls and example protocols to help users plan robust experiments for mapping modification(s) of choice, as well as discussion of the limitations for the methods described.

Full-length tRNAs lacking a functional CCA tail are selectively sorted into the lumen of extracellular vesicles

Small extracellular vesicles (sEVs) are heterogenous lipid membrane particles typically less than 200 nm in size and secreted by most cell types either constitutively or upon activation signals. sEVs isolated from biofluids contain RNAs, including small non-coding RNAs (ncRNAs), that can be either encapsulated within the EV lumen or bound to the EV surface. EV-associated microRNAs (miRNAs) are, despite a relatively low abundance, extensively investigated for their selective incorporation and their role in cell-cell communication. In contrast, the sorting of highly-structured ncRNA species is understudied, mainly due to technical limitations of traditional small RNA sequencing protocols. Here, we adapted ALL-tRNAseq to profile the relative abundance of highly structured and potentially methylated small ncRNA species, including transfer RNAs (tRNAs), small nucleolar RNAs (snoRNAs), and Y RNAs in bulk EV preparations. We determined that full-length tRNAs, typically 75 to 90 nucleotides in length, were the dominant small ncRNA species (>60% of all reads in the 18-120 nucleotides size-range) in all cell culture-derived EVs, as well as in human plasma-derived EV samples, vastly outnumbering 21 nucleotides-long miRNAs. Nearly all EV-associated tRNAs were protected from external RNAse treatment, indicating a location within the EV lumen. Strikingly, the vast majority of luminal-sorted, full-length, nucleobase modification-containing EV-tRNA sequences, harbored a dysfunctional 3' CCA tail, 1 to 3 nucleotides truncated, rendering them incompetent for amino acid loading. In contrast, in non-EV associated extracellular particle fractions (NVEPs), tRNAs appeared almost exclusively fragmented or 'nicked' into tRNA-derived small RNAs (tsRNAs) with lengths between 18 to 35 nucleotides. We propose that in mammalian cells, tRNAs that lack a functional 3' CCA tail are selectively sorted into EVs and shuttled out of the producing cell, offering a new perspective into the physiological role of secreted EVs and luminal cargo-selection.

Large-scale duplication events underpin population-level flexibility in tRNA gene copy number in Pseudomonas fluorescens SBW25

The complement of tRNA genes within a genome is typically considered to be a (relatively) stable characteristic of an organism. Here, we demonstrate that bacterial tRNA gene set composition can be more flexible than previously appreciated, particularly regarding tRNA gene copy number. We report the high-rate occurrence of spontaneous, large-scale, tandem duplication events in laboratory populations of the bacterium Pseudomonas fluorescens SBW25. The identified duplications are up to ∼1 Mb in size (∼15% of the wildtype genome) and are predicted to change the copy number of up to 917 genes, including several tRNA genes. The observed duplications are inherently unstable: they occur, and are subsequently lost, at extremely high rates. We propose that this unusually plastic type of mutation provides a mechanism by which tRNA gene set diversity can be rapidly generated, while simultaneously preserving the underlying tRNA gene set in the absence of continued selection. That is, if a tRNA set variant provides no fitness advantage, then high-rate segregation of the duplication ensures the maintenance of the original tRNA gene set. However, if a tRNA gene set variant is beneficial, the underlying duplication fragment(s) may persist for longer and provide raw material for further, more stable, evolutionary change.

Quantitative and qualitative detection of tRNAs, tRNA halves and tRFs in human cancer samples: Molecular grounds for biomarker development and clinical perspectives

Transfer RNAs (tRNAs) are small non-coding RNAs playing a central role during protein synthesis. Besides translation, growing evidence suggests that in many contexts, precursor or mature tRNAs can also be processed into smaller fragments playing many non-canonical regulatory roles in different biological pathways with oncogenic relevance. Depending on the source, these molecules can be classified as tRNA halves (also known as tiRNAs) or tRNA-derived fragments (tRFs), and furtherly divided into 5'-tRNA and 3'-tRNA halves, or tRF-1, tRF-2, tRF-3, tRF-5, and i-tRF, respectively. Unlike DNA and mRNA, high-throughput sequencing of tRNAs is challenging, because of technical limitations of currently developed sequencing methods. In recent years, different sequencing approaches have been proposed allowing the quantification and identification of an increasing number of tRNA fragments with critical functions in distinct physiological and pathophysiological processes. In the present review, we discussed pros and cons of recent advances in different sequencing methods, also introducing the expanding repertoire of bioinformatics tool and resources specifically focused on tRNA research and discussing current issues in the study of these small RNA molecules. Furthermore, we discussed the potential value of tRNA fragments as diagnostic and prognostic biomarkers for different types of cancers.

Advances in methods for tRNA sequencing and quantification

In the past decade tRNA sequencing (tRNA-seq) has attracted considerable attention as an important tool for the development of novel approaches to quantify highly modified tRNA species and to propel tRNA research aimed at understanding the cellular physiology and disease and development of tRNA-based therapeutics. Many methods are available to quantify tRNA abundance while accounting for modifications and tRNA charging/acylation. Advances in both library preparation methods and bioinformatic workflows have enabled developments in next-generation sequencing (NGS) workflows. Other approaches forgo NGS applications in favor of hybridization-based approaches. In this review we provide a brief comparative overview of various tRNA quantification approaches, focusing on the advantages and disadvantages of these methods, which together facilitate reliable tRNA quantification.

Genome-Wide Profiling of tRNA Using an Unexplored Reverse Transcriptase with High Processivity

Monitoring the dynamic changes of cellular tRNA pools is challenging, due to the extensive post-transcriptional modifications of individual species. The most critical component in tRNAseq is a processive reverse transcriptase (RT) that can read through each modification with high efficiency. Here we show that the recently developed group-II intron RT Induro has the processivity and efficiency necessary to profile tRNA dynamics. Using our Induro-tRNAseq, simpler and more comprehensive than the best methods to date, we show that Induro progressively increases readthrough of tRNA over time and that the mechanism of increase is selective removal of RT stops, without altering the misincorporation frequency. We provide a parallel dataset of the misincorporation profile of Induro relative to the related TGIRT RT to facilitate the prediction of non-annotated modifications. We report an unexpected modification profile among human proline isoacceptors, absent from mouse and lower eukaryotes, that indicates new biology of decoding proline codons.

Locus-specific proteome decoding reveals Fpt1 as a chromatin-associated negative regulator of RNA polymerase III assembly

Transcription of tRNA genes by RNA polymerase III (RNAPIII) is tuned by signaling cascades. The emerging notion of differential tRNA gene regulation implies the existence of additional regulatory mechanisms. However, tRNA gene-specific regulators have not been described. Decoding the local chromatin proteome of a native tRNA gene in yeast revealed reprogramming of the RNAPIII transcription machinery upon nutrient perturbation. Among the dynamic proteins, we identified Fpt1, a protein of unknown function that uniquely occupied RNAPIII-regulated genes. Fpt1 binding at tRNA genes correlated with the efficiency of RNAPIII eviction upon nutrient perturbation and required the transcription factors TFIIIB and TFIIIC but not RNAPIII. In the absence of Fpt1, eviction of RNAPIII was reduced, and the shutdown of ribosome biogenesis genes was impaired upon nutrient perturbation. Our findings provide support for a chromatin-associated mechanism required for RNAPIII eviction from tRNA genes and tuning the physiological response to changing metabolic demands.

tRNA renovatio: Rebirth through fragmentation

tRNA function is based on unique structures that enable mRNA decoding using anticodon trinucleotides. These structures interact with specific aminoacyl-tRNA synthetases and ribosomes using 3D shape and sequence signatures. Beyond translation, tRNAs serve as versatile signaling molecules interacting with other RNAs and proteins. Through evolutionary processes, tRNA fragmentation emerges as not merely random degradation but an act of recreation, generating specific shorter molecules called tRNA-derived small RNAs (tsRNAs). These tsRNAs exploit their linear sequences and newly arranged 3D structures for unexpected biological functions, epitomizing the tRNA "renovatio" (from Latin, meaning renewal, renovation, and rebirth). Emerging methods to uncover full tRNA/tsRNA sequences and modifications, combined with techniques to study RNA structures and to integrate AI-powered predictions, will enable comprehensive investigations of tRNA fragmentation products and new interaction potentials in relation to their biological functions. We anticipate that these directions will herald a new era for understanding biological complexity and advancing pharmaceutical engineering.

Emerging functional principles of tRNA-derived small RNAs and other regulatory small RNAs

Recent advancements in small RNA sequencing have unveiled a previously hidden world of regulatory small noncoding RNAs (sncRNAs) that extend beyond the well-studied small interfering RNAs, microRNAs, and piwi-interacting RNAs. This exploration, starting with tRNA-derived small RNAs, has led to the discovery of a diverse universe of sncRNAs derived from various longer structured RNAs such as rRNAs, small nucleolar RNAs, small nuclear RNAs, Y RNAs, and vault RNAs, with exciting uncharted functional possibilities. In this perspective, we discuss the emerging functional principles of sncRNAs beyond the well-known RNAi-like mechanisms, focusing on those that operate independent of linear sequence complementarity but rather function in an aptamer-like fashion. Aptamers use 3D structure for specific interactions with ligands and are modulated by RNA modifications and subcellular environments. Given that aptamer-like sncRNA functions are widespread and present in species lacking RNAi, they may represent an ancient functional principle that predates RNAi. We propose a rethinking of the origin of RNAi and its relationship with these aptamer-like functions in sncRNAs and how these complementary mechanisms shape biological processes. Lastly, the aptamer-like function of sncRNAs highlights the need for caution in using small RNA mimics in research and therapeutics, as their specificity is not restricted solely to linear sequence.

End-to-end RT-PCR of long RNA and highly structured RNA

RNA molecules play important roles in numerous normal cellular processes and disease states, from protein coding to gene regulation. RT-PCR, applying the power of polymerase chain reaction (PCR) to RNA by coupling reverse transcription with PCR, is one of the most important techniques to characterize RNA transcripts and monitor gene expression. The ability to analyze full-length RNA transcripts and detect their expression is critical to decipher their biological functions. However, due to the low processivity of retroviral reverse transcriptases (RTs), we can only monitor a small fraction of long RNA transcripts, especially those containing stable secondary and tertiary structures. The full-length sequences can only be deduced by computational analysis, which is often misleading. Group II intron-encoded RTs are a new type of RT enzymes. They have evolved specialized structural elements that unwind template structures and maintain close contact with the RNA template. Therefore, group II intron-encoded RTs are more processive than the retroviral RTs. The discovery, optimization and deployment of processive group II intron RTs provide us the opportunity to analyze RNA transcripts with single molecule resolution. MarathonRT, the most processive group II intron RT, has been extensively optimized for processive reverse transcription. In this chapter, we use MarathonRT to deliver a general protocol for long amplicon generation by RT-PCR, and also provide guidance for troubleshooting and further optimization.

NORMSEQ: a tool for evaluation, selection and visualization of RNA-Seq normalization methods

RNA-sequencing has become one of the most used high-throughput approaches to gain knowledge about the expression of all different RNA subpopulations. However, technical artifacts, either introduced during library preparation and/or data analysis, can influence the detected RNA expression levels. A critical step, especially in large and low input datasets or studies, is data normalization, which aims at eliminating the variability in data that is not related to biology. Many normalization methods have been developed, each of them relying on different assumptions, making the selection of the appropriate normalization strategy key to preserve biological information. To address this, we developed NormSeq, a free web-server tool to systematically assess the performance of normalization methods in a given dataset. A key feature of NormSeq is the implementation of information gain to guide the selection of the best normalization method, which is crucial to eliminate or at least reduce non-biological variability. Altogether, NormSeq provides an easy-to-use platform to explore different aspects of gene expression data with a special focus on data normalization to help researchers, even without bioinformatics expertise, to obtain reliable biological inference from their data. NormSeq is freely available at: https://arn.ugr.es/normSeq.