ALL-tRNAseq enables robust tRNA profiling in tissue samples

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.