High-resolution quantitative profiling of tRNA abundance and modification status in eukaryotes by mim-tRNAseq

Rare codon recoding for efficient noncanonical amino acid incorporation in mammalian cells

The ability to genetically encode noncanonical amino acids (ncAAs) has empowered proteins with improved or previously unknown properties. However, existing strategies in mammalian cells rely on the introduction of a blank codon to incorporate ncAAs, which is inefficient and limits their widespread applications. In this study, we developed a rare codon recoding strategy that takes advantage of the relative rarity of the TCG codon to achieve highly selective and efficient ncAA incorporation through systematic engineering and big data-model predictions. We highlight the broad utility of this strategy for the incorporation of dozens of ncAAs into various functional proteins at the wild-type protein expression levels, as well as the synthesis of proteins with up to six-site ncAAs or four distinct ncAAs in mammalian cells for downstream applications.

Comprehensive review on lipid metabolism and RNA methylation: Biological mechanisms, perspectives and challenges

Adipose tissue plays a crucial role in maintaining energy balance, regulating hormones, and promoting metabolic health. To address disorders related to obesity and develop effective therapies, it is essential to have a deep understanding of adipose tissue biology. In recent years, RNA methylation has emerged as a significant epigenetic modification involved in various cellular functions and metabolic pathways. Particularly in the realm of adipogenesis and lipid metabolism, extensive research is ongoing to uncover the mechanisms and functional importance of RNA methylation. Increasing evidence suggests that RNA methylation plays a regulatory role in adipocyte development, metabolism, and lipid utilization across different organs. This comprehensive review aims to provide an overview of common RNA methylation modifications, their occurrences, and regulatory mechanisms, focusing specifically on their intricate connections to fat metabolism. Additionally, we discuss the research methodologies used in studying RNA methylation and highlight relevant databases that can aid researchers in this rapidly advancing field.

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.

Charting new territory: The Plasmodium falciparum tRNA modification landscape

Ribonucleoside modifications comprising the epitranscriptome are present in all organisms and all forms of RNA, including mRNA, rRNA and tRNA, the three major RNA components of the translational machinery. Of these, tRNA is the most heavily modified and the tRNA epitranscriptome has the greatest diversity of modifications. In addition to their roles in tRNA biogenesis, quality control, structure, cleavage, and codon recognition, tRNA modifications have been shown to regulate gene expression post-transcriptionally in prokaryotes and eukaryotes, including humans. However, studies investigating the impact of tRNA modifications on gene expression in the malaria parasite Plasmodium falciparum are currently scarce. Current evidence shows that the parasite has a limited capacity for transcriptional control, which points to a heavier reliance on strategies for posttranscriptional regulation such as tRNA epitranscriptome reprogramming. This review addresses the known functions of tRNA modifications in the biology of P. falciparum while highlighting the potential therapeutic opportunities and the value of using P. falciparum as a model organism for addressing several open questions related to the tRNA epitranscriptome.

Cell-type-specific expression of tRNAs in the brain regulates cellular homeostasis

Defects in tRNA biogenesis are associated with multiple neurological disorders, yet our understanding of these diseases has been hampered by an inability to determine tRNA expression in individual cell types within a complex tissue. Here, we developed a mouse model in which RNA polymerase III is conditionally epitope tagged in a Cre-dependent manner, allowing us to accurately profile tRNA expression in any cell type in vivo. We investigated tRNA expression in diverse nervous system cell types, revealing dramatic heterogeneity in the expression of tRNA genes between populations. We found that while maintenance of levels of tRNA isoacceptor families is critical for cellular homeostasis, neurons are differentially vulnerable to insults to distinct tRNA isoacceptor families. Cell-type-specific translatome analysis suggests that the balance between tRNA availability and codon demand may underlie such differential resilience. Our work provides a platform for investigating the complexities of mRNA translation and tRNA biology in the brain.

Integrated multiplexed assays of variant effect reveal determinants of catechol-O-methyltransferase gene expression

Multiplexed assays of variant effect are powerful methods to profile the consequences of rare variants on gene expression and organismal fitness. Yet, few studies have integrated several multiplexed assays to map variant effects on gene expression in coding sequences. Here, we pioneered a multiplexed assay based on polysome profiling to measure variant effects on translation at scale, uncovering single-nucleotide variants that increase or decrease ribosome load. By combining high-throughput ribosome load data with multiplexed mRNA and protein abundance readouts, we mapped the cis-regulatory landscape of thousands of catechol-O-methyltransferase (COMT) variants from RNA to protein and found numerous coding variants that alter COMT expression. Finally, we trained machine learning models to map signatures of variant effects on COMT gene expression and uncovered both directional and divergent impacts across expression layers. Our analyses reveal expression phenotypes for thousands of variants in COMT and highlight variant effects on both single and multiple layers of expression. Our findings prompt future studies that integrate several multiplexed assays for the readout of gene expression.

Cellular energy regulates mRNA degradation in a codon-specific manner

Codon optimality is a major determinant of mRNA translation and degradation rates. However, whether and through which mechanisms its effects are regulated remains poorly understood. Here we show that codon optimality associates with up to 2-fold change in mRNA stability variations between human tissues, and that its effect is attenuated in tissues with high energy metabolism and amplifies with age. Mathematical modeling and perturbation data through oxygen deprivation and ATP synthesis inhibition reveal that cellular energy variations non-uniformly alter the effect of codon usage. This new mode of codon effect regulation, independent of tRNA regulation, provides a fundamental mechanistic link between cellular energy metabolism and eukaryotic gene expression.

Chemical manipulation of m1A mediates its detection in human tRNA

N 1-methyl adenosine (m1A) is a widespread RNA modification present in tRNA, rRNA, and mRNA. m1A modification sites in tRNAs are evolutionarily conserved and its formation on tRNA is catalyzed by methyltransferase TRMT61A and TRMT6 complex. m1A promotes translation initiation and elongation. Due to its positive charge under physiological conditions, m1A can notably modulate RNA structure. It also blocks Watson-Crick-Franklin base-pairing and causes mutation and truncation during reverse transcription. Several misincorporation-based high-throughput sequencing methods have been developed to sequence m1A. In this study, we introduce a reduction-based m1A sequencing (red-m1A-seq). We report that NaBH4 reduction of m1A can improve the mutation and readthrough rates using commercially available RT enzymes to give a better positive signature, while alkaline-catalyzed Dimroth rearrangement can efficiently convert m1A to m6A to provide good controls, allowing the detection of m1A with higher sensitivity and accuracy. We applied red-m1A-seq to sequence human small RNA, and we not only detected all the previously reported tRNA m1A sites, but also new m1A sites in mt-tRNAAsn-GTT and 5.8S rRNA.

Bioorthogonal labeling and profiling of N6-isopentenyladenosine (i6A) modified RNA

Chemical modifications in RNAs play crucial roles in diversifying their structures and regulating numerous biochemical processes. Since the 1990s, several hydrophobic prenyl-modifications have been discovered in various RNAs. Prenyl groups serve as precursors for terpenes and many other biological molecules. The processes of prenylation in different macromolecules have been extensively studied. We introduce here a novel chemical biology toolkit that not only labels i6A, a prenyl-modified RNA residue, by leveraging the unique reactivity of the prenyl group, but also provides a general strategy to incorporate fluorescence functionalities into RNAs for molecular tracking purposes. Our findings revealed that iodine-mediated cyclization reactions of the prenyl group occur rapidly, transforming i6A from a hydrogen-bond acceptor to a donor. Based on this reactivity, we developed an Iodine-Mediated Cyclization and Reverse Transcription (IMCRT) tRNA-seq method, which can profile all nine endogenous tRNAs containing i6A residues in Saccharomyces cerevisiae with single-base resolution. Furthermore, under stress conditions, we observed a decline in i6A levels in budding yeast, accompanied by significant decrease of mutation rate at A37 position. Thus, the IMCRT tRNA-seq method not only permits semi-quantification of i6A levels in tRNAs but also holds potential for transcriptome-wide detection and analysis of various RNA species containing i6A modifications.

Translation machinery: the basis of translational control

Messenger RNA (mRNA) translation consists of initiation, elongation, termination, and ribosome recycling, carried out by the translation machinery, primarily including tRNAs, ribosomes, and translation factors (TrFs). Translational regulators transduce signals of growth and development, as well as biotic and abiotic stresses, to the translation machinery, where global or selective translational control occurs to modulate mRNA translation efficiency (TrE). As the basis of translational control, the translation machinery directly determines the quality and quantity of newly synthesized peptides and, ultimately, the cellular adaption. Thus, regulating the availability of diverse machinery components is reviewed as the central strategy of translational control. We provide classical signaling pathways (e.g., integrated stress responses) and cellular behaviors (e.g., liquid-liquid phase separation) to exemplify this strategy within different physiological contexts, particularly during host-microbe interactions. With new technologies developed, further understanding this strategy will speed up translational medicine and translational agriculture.

The utility of Drosophila melanogaster as a fungal infection model

Invasive fungal diseases have profound effects upon human health and are on increase globally. The World Health Organization (WHO) in 2022 published the fungal priority list calling for improved public health interventions and advance research. Drosophila melanogaster presents an excellent model system to dissect host-pathogen interactions and has been proved valuable to study immunopathogenesis of fungal diseases. In this review we highlight the recent advances in fungal-Drosophila interplay with an emphasis on the recently published WHO's fungal priority list and we focus on available tools and technologies.

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.

A robust method for measuring aminoacylation through tRNA-Seq

Current methods to quantify the fraction of aminoacylated tRNAs, also known as the tRNA charge, are limited by issues with either low throughput, precision, and/or accuracy. Here, we present an optimized charge tRNA-Seq method that combines previous developments with newly described approaches to establish a protocol for precise and accurate tRNA charge measurements. We verify that this protocol provides robust quantification of tRNA aminoacylation and we provide an end-to-end method that scales to hundreds of samples including software for data processing. Additionally, we show that this method supports measurements of relative tRNA expression levels and can be used to infer tRNA modifications through reverse transcription misincorporations, thereby supporting multipurpose applications in tRNA biology.

Neurodevelopmental disorders associated variants in ADAT3 disrupt the activity of the ADAT2/ADAT3 tRNA deaminase complex and impair neuronal migration

The ADAT2/ADAT3 complex catalyzes the adenosine to inosine modification at the wobble position of eukaryotic tRNAs. Mutations in ADAT3 , the catalytically inactive subunit of the ADAT2/ADAT3 complex, have been identified in patients presenting with severe neurodevelopmental disorders (NDDs). Yet, the physiological function of ADAT2/ADAT3 complex during brain development remains totally unknown. Here we showed that maintaining a proper level of ADAT2/ADAT3 catalytic activity is required for correct radial migration of projection neurons in the developing mouse cortex. In addition, we not only reported 7 new NDD patients carrying biallelic variants in ADAT3 but also deeply characterize the impact of those variants on ADAT2/ADAT3 structure, biochemical properties, enzymatic activity and tRNAs editing and abundance. We demonstrated that all the identified variants alter both the expression and the activity of the complex leading to a significant decrease of I 34 with direct consequence on their steady-state. Using in vivo complementation assays, we correlated the severity of the migration phenotype with the degree of the loss of function caused by the variants. Altogether, our results indicate a critical role of ADAT2/ADAT3 during cortical development and provide cellular and molecular insights into the pathogenicity of ADAT3-related neurodevelopmental disorder.

Small RNA structural biochemistry in a post-sequencing era

High-throughput sequencing has had an enormous impact on small RNA research during the past decade. However, sequencing only offers a one-dimensional view of the transcriptome and is often highly biased. Additionally, the 'sequence, map and annotate' approach, used widely in small RNA research, can lead to flawed interpretations of the data, lacking biological plausibility, due in part to database issues. Even in the absence of technical biases, the loss of three-dimensional information is a major limitation to understanding RNA stability, turnover and function. For example, noncoding RNA-derived fragments seem to exist mainly as dimers, tetramers or as nicked forms of their parental RNAs, contrary to widespread assumptions. In this perspective, we will discuss main sources of bias during small RNA-sequencing, present several useful bias-reducing strategies and provide guidance on the interpretation of small RNA-sequencing results, with emphasis on RNA fragmentomics. As sequencing offers a one-dimensional projection of a four-dimensional reality, prior structure-level knowledge is often needed to make sense of the data. Consequently, while less-biased sequencing methods are welcomed, integration of orthologous experimental techniques is also strongly recommended.

Deep generative design of RNA family sequences

RNA engineering has immense potential to drive innovation in biotechnology and medicine. Despite its importance, a versatile platform for the automated design of functional RNA is still lacking. Here, we propose RNA family sequence generator (RfamGen), a deep generative model that designs RNA family sequences in a data-efficient manner by explicitly incorporating alignment and consensus secondary structure information. RfamGen can generate novel and functional RNA family sequences by sampling points from a semantically rich and continuous representation. We have experimentally demonstrated the versatility of RfamGen using diverse RNA families. Furthermore, we confirmed the high success rate of RfamGen in designing functional ribozymes through a quantitative massively parallel assay. Notably, RfamGen successfully generates artificial sequences with higher activity than natural sequences. Overall, RfamGen significantly improves our ability to design functional RNA and opens up new potential for generative RNA engineering in synthetic biology.

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.

Ribosome inactivation regulates translation elongation in neurons

Cellular plasticity is crucial for adapting to ever-changing stimuli. As a result, cells consistently reshape their translatome, and, consequently, their proteome. The control of translational activity has been thoroughly examined at the stage of translation initiation. However, the regulation of ribosome speed in cells is widely unknown. In this study, we utilized a timed ribosome runoff approach, along with proteomics and transmission electron microscopy, to investigate global translation kinetics in cells. We found that ribosome speeds vary among various cell types, such as astrocytes, induced pluripotent human stem cells, human neural stem cells, and human and rat neurons. Of all cell types studied, mature cortical neurons exhibit the highest rate of translation. This finding is particularly remarkable because mature cortical neurons express the eukaryotic elongation factor 2 (eEF2) at lower levels than other cell types. Neurons solve this conundrum by inactivating a fraction of their ribosomes. As a result, the increase in eEF2 levels leads to a reduction of inactive ribosomes and an enhancement of active ones. Processes that alter the demand for active ribosomes, like neuronal excitation, cause increased inactivation of redundant ribosomes in an eEF2-dependent manner. Our data suggest a novel regulatory mechanism in which neurons dynamically inactivate ribosomes to facilitate translational remodeling. These findings have important implications for developmental brain disorders characterized by, among other things, aberrant translation.

tRNA therapeutics for genetic diseases

Transfer RNAs (tRNAs) have a crucial role in protein synthesis, and in recent years, their therapeutic potential for the treatment of genetic diseases - primarily those associated with a mutation altering mRNA translation - has gained significant attention. Engineering tRNAs to readthrough nonsense mutation-associated premature termination of mRNA translation can restore protein synthesis and function. In addition, supplementation of natural tRNAs can counteract effects of missense mutations in proteins crucial for tRNA biogenesis and function in translation. This Review will present advances in the development of tRNA therapeutics with high activity and safety in vivo and discuss different formulation approaches for single or chronic treatment modalities. The field of tRNA therapeutics is still in its early stages, and a series of challenges related to tRNA efficacy and stability in vivo, delivery systems with tissue-specific tropism, and safe and efficient manufacturing need to be addressed.

Quantitative analysis of tRNA abundance and modifications by nanopore RNA sequencing

Transfer RNAs (tRNAs) play a central role in protein translation. Studying them has been difficult in part because a simple method to simultaneously quantify their abundance and chemical modifications is lacking. Here we introduce Nano-tRNAseq, a nanopore-based approach to sequence native tRNA populations that provides quantitative estimates of both tRNA abundances and modification dynamics in a single experiment. We show that default nanopore sequencing settings discard the vast majority of tRNA reads, leading to poor sequencing yields and biased representations of tRNA abundances based on their transcript length. Re-processing of raw nanopore current intensity signals leads to a 12-fold increase in the number of recovered tRNA reads and enables recapitulation of accurate tRNA abundances. We then apply Nano-tRNAseq to Saccharomyces cerevisiae tRNA populations, revealing crosstalks and interdependencies between different tRNA modification types within the same molecule and changes in tRNA populations in response to oxidative stress.

Selective gene expression maintains human tRNA anticodon pools during differentiation

Transfer RNAs are essential for translating genetic information into proteins. The human genome contains hundreds of predicted tRNA genes, many in multiple copies. How their expression is regulated to control tRNA repertoires is unknown. Here we combined quantitative tRNA profiling and chromatin immunoprecipitation with sequencing to measure tRNA expression following the differentiation of human induced pluripotent stem cells into neuronal and cardiac cells. We find that tRNA transcript levels vary substantially, whereas tRNA anticodon pools, which govern decoding rates, are more stable among cell types. Mechanistically, RNA polymerase III transcribes a wide range of tRNA genes in human induced pluripotent stem cells but on differentiation becomes constrained to a subset we define as housekeeping tRNAs. This shift is mediated by decreased mTORC1 signalling, which activates the RNA polymerase III repressor MAF1. Our data explain how tRNA anticodon pools are buffered to maintain decoding speed across cell types and reveal that mTORC1 drives selective tRNA expression during differentiation.

Tailor made: the art of therapeutic mRNA design

mRNA medicine is a new and rapidly developing field in which the delivery of genetic information in the form of mRNA is used to direct therapeutic protein production in humans. This approach, which allows for the quick and efficient identification and optimization of drug candidates for both large populations and individual patients, has the potential to revolutionize the way we prevent and treat disease. A key feature of mRNA medicines is their high degree of designability, although the design choices involved are complex. Maximizing the production of therapeutic proteins from mRNA medicines requires a thorough understanding of how nucleotide sequence, nucleotide modification and RNA structure interplay to affect translational efficiency and mRNA stability. In this Review, we describe the principles that underlie the physical stability and biological activity of mRNA and emphasize their relevance to the myriad considerations that factor into therapeutic mRNA design.

RNA polymerases reshape chromatin and coordinate transcription on individual fibers

During eukaryotic transcription, RNA polymerases must initiate and pause within a crowded, complex environment, surrounded by nucleosomes and other transcriptional activity. This environment creates a spatial arrangement along individual chromatin fibers ripe for both competition and coordination, yet these relationships remain largely unknown owing to the inherent limitations of traditional structural and sequencing methodologies. To address these limitations, we employed long-read chromatin fiber sequencing (Fiber-seq) to visualize RNA polymerases within their native chromatin context at single-molecule and near single-nucleotide resolution along up to 30 kb fibers. We demonstrate that Fiber-seq enables the identification of single-molecule RNA Polymerase (Pol) II and III transcription associated footprints, which, in aggregate, mirror bulk short-read sequencing-based measurements of transcription. We show that Pol II pausing destabilizes downstream nucleosomes, with frequently paused genes maintaining a short-term memory of these destabilized nucleosomes. Furthermore, we demonstrate pervasive direct coordination and anti-coordination between nearby Pol II genes, Pol III genes, transcribed enhancers, and insulator elements. This coordination is largely limited to spatially organized elements within 5 kb of each other, implicating short-range chromatin environments as a predominant determinant of coordinated polymerase initiation. Overall, transcription initiation reshapes surrounding nucleosome architecture and coordinates nearby transcriptional machinery along individual chromatin fibers.

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.

A tRNA-modifying enzyme facilitates RNase P activity in Arabidopsis nuclei

RNase P is the essential activity that performs the 5' maturation of transfer RNA (tRNA) precursors. Beyond the ancestral form of RNase P containing a ribozyme, protein-only RNase P enzymes termed PRORP were identified in eukaryotes. In human mitochondria, PRORP forms a complex with two protein partners to become functional. In plants, although PRORP enzymes are active alone, we investigate their interaction network to identify potential tRNA maturation complexes. Here we investigate functional interactions involving the Arabidopsis nuclear RNase P PRORP2. We show, using an immuno-affinity strategy, that PRORP2 occurs in a complex with the tRNA methyl transferases TRM1A and TRM1B in vivo. Beyond RNase P, these enzymes can also interact with RNase Z. We show that TRM1A/TRM1B localize in the nucleus and find that their double knockout mutation results in a severe macroscopic phenotype. Using a combination of immuno-detections, mass spectrometry and a transcriptome-wide tRNA sequencing approach, we observe that TRM1A/TRM1B are responsible for the m22G26 modification of 70% of cytosolic tRNAs in vivo. We use the transcriptome wide tRNAseq approach as well as RNA blot hybridizations to show that RNase P activity is impaired in TRM1A/TRM1B mutants for specific tRNAs, in particular, tRNAs containing a m22G modification at position 26 that are strongly downregulated in TRM1A/TRM1B mutants. Altogether, results indicate that the m22G-adding enzymes TRM1A/TRM1B functionally cooperate with nuclear RNase P in vivo for the early steps of cytosolic tRNA biogenesis.

tRNA modifications: insights into their role in human cancers

Transfer RNA (tRNA) plays a central role in translation by functioning as a biological link between messenger RNA (mRNA) and proteins. One prominent feature of the tRNA molecule is its heavily modified status, which greatly affects its biogenesis and function. Modifications within the anticodon loop are crucial for translation efficiency and accuracy, whereas other modifications in the body region affect tRNA structure and stability. Recent research has revealed that these diverse modifications are critical regulators of gene expression. They are involved in many important physiological and pathological processes, including cancers. In this review we focus on six different tRNA modifications to delineate their functions and mechanisms in tumorigenesis and tumor progression, providing insights into their clinical potential as biomarkers and therapeutic targets.

Integrated multiplexed assays of variant effect reveal cis-regulatory determinants of catechol-O-methyltransferase gene expression

Multiplexed assays of variant effect are powerful methods to profile the consequences of rare variants on gene expression and organismal fitness. Yet, few studies have integrated several multiplexed assays to map variant effects on gene expression in coding sequences. Here, we pioneered a multiplexed assay based on polysome profiling to measure variant effects on translation at scale, uncovering single-nucleotide variants that increase and decrease ribosome load. By combining high-throughput ribosome load data with multiplexed mRNA and protein abundance readouts, we mapped the cis-regulatory landscape of thousands of catechol-O-methyltransferase (COMT) variants from RNA to protein and found numerous coding variants that alter COMT expression. Finally, we trained machine learning models to map signatures of variant effects on COMT gene expression and uncovered both directional and divergent impacts across expression layers. Our analyses reveal expression phenotypes for thousands of variants in COMT and highlight variant effects on both single and multiple layers of expression. Our findings prompt future studies that integrate several multiplexed assays for the readout of gene expression.

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.

tRNA-derived small RNAs are embedded in the gene regulatory program instructing Drosophila metamorphosis

A class of noncoding RNAs, referred to as tsRNAs, is emerging with a potential to exert a new layer in gene regulation. These RNAs are breakdown products of tRNAs, either through active processing or passive cleavage or both. Since tRNAs are part of the general machinery for translation, their expression levels and activities are tightly controlled, raising the possibility that their breakdown products, tsRNAs, may provide a link between the overall translational status of a cell to specific changes in gene regulatory network. We hypothesize that Drosophila pupation, being a special developmental stage during which there is a global limitation of nutrients, represents a system in which such a link may readily reveal itself. We show that specific tsRNAs indeed exhibit a dynamic accumulation upon entering the pupal stage. We describe experiments to characterize the mode of tsRNA action and, through the use of such gained knowledge, conduct a genome-wide analysis to assess the functions of dynamically expressed tsRNAs. Our results show that the predicted target genes are highly enriched in biological processes specific to this stage of development including metamorphosis. We further show that tsRNA action is required for successful pupation, providing direct support to the hypothesis that tsRNAs accumulated during this stage are critical to the gene expression program at this stage of development.

Transfer RNAs as dynamic and critical regulators of cancer progression

Transfer RNAs (tRNAs) have been historically viewed as non-dynamic adaptors that decode the genetic code into proteins. Recent work has uncovered dynamic regulatory roles for these fascinating molecules. Advances in tRNA detection methods have revealed that specific tRNAs can become modulated upon DNA copy number and chromatin alterations and can also be perturbed by oncogenic signalling and transcriptional regulators in cancer cells or the tumour microenvironment. Such alterations in the levels of specific tRNAs have been shown to causally impact cancer progression, including metastasis. Moreover, sequencing methods have identified tRNA-derived small RNAs that influence various aspects of cancer progression, such as cell proliferation and invasion, and could serve as diagnostic and prognostic biomarkers or putative therapeutic targets in various cancers. Finally, there is accumulating evidence, including from genetic models, that specific tRNA synthetases - the enzymes responsible for charging tRNAs with amino acids - can either promote or suppress tumour formation. In this Review, we provide an overview of how deregulation of tRNAs influences cancer formation and progression.

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.

Ushering in the era of tRNA medicines

Long viewed as an intermediary in protein translation, there is a growing awareness that tRNAs are capable of myriad other biological functions linked to human health and disease. These emerging roles could be tapped to leverage tRNAs as diagnostic biomarkers, therapeutic targets, or even as novel medicines. Furthermore, the growing array of tRNA-derived fragments, which modulate an increasingly broad spectrum of cellular pathways, is expanding this opportunity. Together, these molecules offer drug developers the chance to modulate the impact of mutations and to alter cell homeostasis. Moreover, because a single therapeutic tRNA can facilitate readthrough of a genetic mutation shared across multiple genes, such medicines afford the opportunity to define patient populations not based on their clinical presentation or mutated gene but rather on the mutation itself. This approach could potentially transform the treatment of patients with rare and ultrarare diseases. In this review, we explore the diverse biology of tRNA and its fragments, examining the past and present challenges to provide a comprehensive understanding of the molecules and their therapeutic potential.

Human TRMT2A methylates tRNA and contributes to translation fidelity

5-Methyluridine (m5U) is one of the most abundant RNA modifications found in cytosolic tRNA. tRNA methyltransferase 2 homolog A (hTRMT2A) is the dedicated mammalian enzyme for m5U formation at tRNA position 54. However, its RNA binding specificity and functional role in the cell are not well understood. Here we dissected structural and sequence requirements for binding and methylation of its RNA targets. Specificity of tRNA modification by hTRMT2A is achieved by a combination of modest binding preference and presence of a uridine in position 54 of tRNAs. Mutational analysis together with cross-linking experiments identified a large hTRMT2A-tRNA binding surface. Furthermore, complementing hTRMT2A interactome studies revealed that hTRMT2A interacts with proteins involved in RNA biogenesis. Finally, we addressed the question of the importance of hTRMT2A function by showing that its knockdown reduces translation fidelity. These findings extend the role of hTRMT2A beyond tRNA modification towards a role in translation.

Effect of mRNA/tRNA mutations on translation speed: Implications for human diseases

Recent discoveries establish tRNAs as central regulators of mRNA translation dynamics, and therefore cotranslational folding and function of the encoded protein. The tRNA pool, whose composition and abundance change in a cell- and tissue-dependent manner, is the main factor which determines mRNA translation velocity. In this review, we discuss a group of pathogenic mutations, in the coding sequences of either protein-coding genes or in tRNA genes, that alter mRNA translation dynamics. We also summarize advances in tRNA biology that have uncovered how variations in tRNA levels on account of genetic mutations affect protein folding and function, and thereby contribute to phenotypic diversity in clinical manifestations.

Dynamic Landscapes of tRNA Transcriptomes and Translatomes in Diverse Mouse Tissues

Although the function of tRNAs in the translational process is well established, it remains controversial whether tRNA abundance is tightly associated with translational efficiency (TE) in mammals. Moreover, how critically the expression of tRNAs contributes to the establishment of tissue-specific proteomes in mammals has not been well addressed. Here, we measured both tRNA expression using demethylase-tRNA sequencing (DM-tRNA-seq) and TE of mRNAs using ribosome-tagging sequencing (RiboTag-seq) in the brain, heart, and testis of mice. Remarkable variation in the expression of tRNA isodecoders was observed among different tissues. When the statistical effect of isodecoder-grouping on reducing variations is considered through permutating the anticodons, we observed an expected reduction in the variation of anticodon expression across all samples, an unexpected smaller variation of anticodon usage bias, and an unexpected larger variation of tRNA isotype expression at amino acid level. Regardless of whether or not they share the same anticodons, the isodecoders encoding the same amino acids are co-expressed across different tissues. Based on the expression of tRNAs and the TE of mRNAs, we find that the tRNA adaptation index (tAI) and TE are significantly correlated in the same tissues but not between tissues; and tRNA expression and the amino acid composition of translating peptides are positively correlated in the same tissues but not between tissues. We therefore hypothesize that the tissue-specific expression of tRNAs might be due to post-transcriptional mechanisms. This study provides a resource for tRNA and translation studies, as well as novel insights into the dynamics of tRNAs and their roles in translational regulation.

Selection of tRNA Genes in Human Breast Tumours Varies Substantially between Individuals

Abnormally elevated expression of tRNA is a common feature of breast tumours. Rather than a uniform increase in all tRNAs, some are deregulated more strongly than others. Elevation of particular tRNAs has been associated with poor prognosis for patients, and experimental models have demonstrated the ability of some tRNAs to promote proliferation or metastasis. Each tRNA isoacceptor is encoded redundantly by multiple genes, which are commonly dispersed across several chromosomes. An unanswered question is whether the consistently high expression of a tRNA in a cancer type reflects the consistent activation of the same members of a gene family, or whether different family members are activated from one patient to the next. To address this question, we interrogated ChIP-seq data to determine which tRNA genes were active in individual breast tumours. This revealed that distinct sets of tRNA genes become activated in individual cancers, whereas there is much less variation in the expression patterns of families. Several pathways have been described that are likely to contribute to increases in tRNA gene transcription in breast tumours, but none of these can adequately explain the observed variation in the choice of genes between tumours. Current models may therefore lack at least one level of regulation.

The central role of translation elongation in response to stress

Protein synthesis is essential to support homeostasis, and thus, must be highly regulated during cellular response to harmful environments. All stages of translation are susceptible to regulation under stress, however, the mechanisms involved in translation regulation beyond initiation have only begun to be elucidated. Methodological advances enabled critical discoveries on the control of translation elongation, highlighting its important role in translation repression and the synthesis of stress-response proteins. In this article, we discuss recent findings on mechanisms of elongation control mediated by ribosome pausing and collisions and the availability of tRNAs and elongation factors. We also discuss how elongation intersects with distinct modes of translation control, further supporting cellular viability and gene expression reprogramming. Finally, we highlight how several of these pathways are reversibly regulated, emphasizing the dynamics of translation control during stress-response progression. A comprehensive understanding of translation regulation under stress will produce fundamental knowledge of protein dynamics while opening new avenues and strategies to overcome dysregulated protein production and cellular sensitivity to stress.

A ligation-independent sequencing method reveals tRNA-derived RNAs with blocked 3' termini

Despite the numerous sequencing methods available, the vast diversity in size and chemical modifications of RNA molecules makes the capture of the full spectrum of cellular RNAs a difficult task. By combining quasi-random hexamer priming with a custom template switching strategy, we developed a method to construct sequencing libraries from RNA molecules of any length and with any type of 3' terminal modification, allowing the sequencing and analysis of virtually all RNA species. Ligation-independent detection of all types of RNA (LIDAR) is a simple, effective tool to comprehensively characterize changes in small non-coding RNAs and mRNAs simultaneously, with performance comparable to separate dedicated methods. With LIDAR, we comprehensively characterized the coding and non-coding transcriptome of mouse embryonic stem cells, neural progenitor cells, and sperm. LIDAR detected a much larger variety of tRNA-derived RNAs (tDRs) compared to traditional ligation-dependent sequencing methods, and uncovered the presence of tDRs with blocked 3' ends that had previously escaped detection. Our findings highlight the potential of LIDAR to systematically detect all RNAs in a sample and uncover new RNA species with potential regulatory functions.

Fluorescent labeling of tRNA for rapid kinetic interaction studies with tRNA-binding proteins

Transfer RNA (tRNA) plays a critical role during translation and interacts with numerous proteins during its biogenesis, functional cycle and degradation. In particular, tRNA is extensively post-transcriptionally modified by various tRNA modifying enzymes which each target a specific nucleotide at different positions within tRNAs to introduce different chemical modifications. Fluorescent assays can be used to study the interaction between a protein and tRNA. Moreover, rapid mixing fluorescence stopped-flow assays provide insights into the kinetics of the tRNA-protein interaction in order to elucidate the tRNA binding mechanism for the given protein. A prerequisite for these studies is a fluorescently labeled molecule, such as fluorescent tRNA, wherein a change in fluorescence occurs upon protein binding. In this chapter, we discuss the utilization of tRNA modifications in order to introduce fluorophores at particular positions within tRNAs. Particularly, we focus on in vitro thiolation of a uridine at position 8 within tRNAs using the tRNA modification enzyme ThiI, followed by labeling of the thiol group with fluorescein. As such, this fluorescently labeled tRNA is primarily unmodified, with the exception of the thiolation modification to which the fluorophore is attached, and can be used as a substrate to study the binding of different tRNA-interacting factors. Herein, we discuss the example of studying the tRNA binding mechanism of the tRNA modifying enzymes TrmB and DusA using internally fluorescein-labeled tRNA.

Methylated guanosine and uridine modifications in S. cerevisiae mRNAs modulate translation elongation

Chemical modifications to protein encoding messenger RNAs (mRNAs) influence their localization, translation, and stability within cells. Over 15 different types of mRNA modifications have been observed by sequencing and liquid chromatography coupled to tandem mass spectrometry (LC-MS/MS) approaches. While LC-MS/MS is arguably the most essential tool available for studying analogous protein post-translational modifications, the high-throughput discovery and quantitative characterization of mRNA modifications by LC-MS/MS has been hampered by the difficulty of obtaining sufficient quantities of pure mRNA and limited sensitivities for modified nucleosides. We have overcome these challenges by improving the mRNA purification and LC-MS/MS pipelines. The methodologies we developed result in no detectable non-coding RNA modifications signals in our purified mRNA samples, quantify 50 ribonucleosides in a single analysis, and provide the lowest limit of detection reported for ribonucleoside modification LC-MS/MS analyses. These advancements enabled the detection and quantification of 13 S. cerevisiae mRNA ribonucleoside modifications and reveal the presence of four new S. cerevisiae mRNA modifications at low to moderate levels (1-methyguanosine, N2-methylguanosine, N2,N2-dimethylguanosine, and 5-methyluridine). We identified four enzymes that incorporate these modifications into S. cerevisiae mRNAs (Trm10, Trm11, Trm1, and Trm2, respectively), though our results suggest that guanosine and uridine nucleobases are also non-enzymatically methylated at low levels. Regardless of whether they are incorporated in a programmed manner or as the result of RNA damage, we reasoned that the ribosome will encounter the modifications that we detect in cells. To evaluate this possibility, we used a reconstituted translation system to investigate the consequences of modifications on translation elongation. Our findings demonstrate that the introduction of 1-methyguanosine, N2-methylguanosine and 5-methyluridine into mRNA codons impedes amino acid addition in a position dependent manner. This work expands the repertoire of nucleoside modifications that the ribosome must decode in S. cerevisiae. Additionally, it highlights the challenge of predicting the effect of discrete modified mRNA sites on translation de novo because individual modifications influence translation differently depending on mRNA sequence context.

satmut_utils: a simulation and variant calling package for multiplexed assays of variant effect

The impact of millions of individual genetic variants on molecular phenotypes in coding sequences remains unknown. Multiplexed assays of variant effect (MAVEs) are scalable methods to annotate relevant variants, but existing software lacks standardization, requires cumbersome configuration, and does not scale to large targets. We present satmut_utils as a flexible solution for simulation and variant quantification. We then benchmark MAVE software using simulated and real MAVE data. We finally determine mRNA abundance for thousands of cystathionine beta-synthase variants using two experimental methods. The satmut_utils package enables high-performance analysis of MAVEs and reveals the capability of variants to alter mRNA abundance.

High production of transfer RNAs identifies the presence of developing oocytes in ovaries and intersex testes of teleost fish

5S rRNA is highly transcribed in fish oocytes and this transcription levels can be used to identify the presence of oocytes in the intersex testes of fish exposed to xenoestrogens. Similar to 5S rRNA, tRNAs are transcribed by RNA polymerase III (Pol-III) in eukaryotes, so this study focuses in the analysis of the levels of expression of tRNAs in the gonads (ovaries and testes) of eight teleost species as a possible new oocyte molecular marker. Total RNA extracted from gonads of six commercial teleost species in the Biscay Bay, from the pollution sentinel species thicklip grey mullet (Chelon labrosus) known present intersex testes in response to xenoestrogens in Gernika estuary and from the laboratory model species Danio rerio were analysed through capillary electrophoresis. Bioanalyzer electropherograms were used to quantify the concentrations of tRNAs, 5S and 5.8S rRNA. All studied ovaries expressed significantly higher levels of tRNAs and 5S rRNA than testes. A tRNA to 5.8S rRNA index was calculated which differentiates ovaries from testes, and identifies some intersex testes in between testes and ovaries in mullets. The tRNA/5.8S ratio was highest in ovaries in previtellogenic stage, decreasing towards maturity. Thus, strong oocyte expression of tRNAs is an additional proof of high activity levels of Pol-III during early stages of oocyte development in teleost ovaries. Incidentally, we observed that miRNA concentrations were always higher in testes than ovaries. The indexing approach developed in the present study could have multiple applications in teleost reproduction research and in the development of early molecular markers of intersex condition.

The role and mechanism of action of tRNA-derived fragments in the diagnosis and treatment of malignant tumors

Cancer is a leading cause of morbidity and death worldwide. While various factors are established as causing malignant tumors, the mechanisms underlying cancer development remain poorly understood. Early diagnosis and the development of effective treatments for cancer are important research topics. Transfer RNA (tRNA), the most abundant class of RNA molecules in the human transcriptome, participates in both protein synthesis and cellular metabolic processes. tRNA-derived fragments (tRFs) are produced by specific cleavage of pre-tRNA and mature tRNA molecules, which are highly conserved and occur widely in various organisms. tRFs were initially thought to be random products with no physiological function, but have been redefined as novel functional small non-coding RNA molecules that help to regulate RNA stability, modulate translation, and influence target gene expression, as well as other biological processes. There is increasing evidence supporting roles for tRFs in tumorigenesis and cancer development, including the regulation of tumor cell proliferation, invasion, migration, and drug resistance. Understanding the regulatory mechanisms by which tRFs impact these processes has potential to inform malignant tumor diagnosis and treatment. Further, tRFs are expected to become new biological markers for early diagnosis and prognosis prediction in patients with tumors, as well as a targets for precision cancer therapies. Video abstract.

Functional analysis of tRNA modification enzymes using mutational profiling

Transfer RNA (tRNA) delivers amino acids to the ribosome and functions as an essential adapter molecule for decoding codons on the messenger RNA (mRNA) during protein synthesis. Before attaining their proper activity, tRNAs undergo multiple post-transcriptional modifications with highly diversified roles such as stabilization of the tRNA structure, recognition of aminoacyl tRNA synthetases, precise codon-anticodon recognition, support of viral replication and onset of immune responses. The synthesis of the majority of modified nucleosides is catalyzed by a site-specific tRNA modification enzyme. This chapter provides a detailed protocol for using mutational profiling to analyze the enzymatic function of a tRNA methyltransferase in a high-throughput manner. In a previous study, we took tRNA m1A22 methyltransferase TrmK from Geobacillus stearothermophilus as a model tRNA methyltransferase and applied this protocol to gain mechanistic insights into how TrmK recognizes the substrate tRNAs. In theory, this protocol can be used unaltered for studying enzymes that catalyze modifications at the Watson-Crick face such as 1-methyladenosine (m1A), 3-methylcytosine (m3C), 3-methyluridine (m3U), 1-methylguanosine (m1G), and N2,N2-dimethylguanosine (m22G).

Sperm RNA Payload: Implications for Intergenerational Epigenetic Inheritance

There is mounting evidence that ancestral life experiences and environment can influence phenotypes in descendants. The parental environment regulates offspring phenotypes potentially via modulating epigenetic marks in the gametes. Here, we review examples of across-generational inheritance of paternal environmental effects and the current understanding of the role of small RNAs in such inheritance. We discuss recent advances in revealing the small RNA payload of sperm and how environmental conditions modulate sperm small RNAs. Further, we discuss the potential mechanism of inheritance of paternal environmental effects by focusing on sperm small RNA-mediated regulation of early embryonic gene expression and its role in influencing offspring phenotypes.

Elucidation of transfer RNAs as stress regulating agents and the experimental strategies to conceive the functional role of tRNA-derived fragments in plants

In plants, the transfer RNAs (tRNAs) exhibit their profound influence in orchestrating diverse physiological activities like cell growth, development, and response to several surrounding stimuli. The tRNAs, which were known to restrict their function solely in deciphering the codons, are now emerging as frontline defenders in stress biology. The plants that are constantly confronted with a huge panoply of stresses rely on tRNA-mediated stress regulation by altering the tRNA abundance, curbing the transport of tRNAs, fragmenting the mature tRNAs during stress. Among them, the studies on the generation of transfer RNA-derived fragments (tRFs) and their biological implication in stress response have attained huge interest. In plants, the tRFs hold stable expression patterns and regulate biological functions under diverse environmental conditions. In this review, we discuss the fate of plant tRNAs upon stress and thereafter how the tRFs are metamorphosed into sharp ammunition to wrestle with stress. We also address the various methods developed to date for uncovering the role of tRFs and their function in plants.

ALL-tRNAseq enables robust tRNA profiling in tissue samples

Transfer RNAs (tRNAs) are small adaptor RNAs essential for mRNA translation. Alterations in the cellular tRNA population can directly affect mRNA decoding rates and translational efficiency during cancer development and progression. To evaluate changes in the composition of the tRNA pool, multiple sequencing approaches have been developed to overcome reverse transcription blocks caused by the stable structures of these molecules and their numerous base modifications. However, it remains unclear whether current sequencing protocols faithfully capture tRNAs existing in cells or tissues. This is specifically challenging for clinical tissue samples that often present variable RNA qualities. For this reason, we developed ALL-tRNAseq, which combines the highly processive MarathonRT and RNA demethylation for the robust assessment of tRNA expression, together with a randomized adapter ligation strategy prior to reverse transcription to assess tRNA fragmentation levels in both cell lines and tissues. Incorporation of tRNA fragments not only informed on sample integrity but also significantly improved tRNA profiling of tissue samples. Our data showed that our profiling strategy effectively improves classification of oncogenic signatures in glioblastoma and diffuse large B-cell lymphoma tissues, particularly for samples presenting higher levels of RNA fragmentation, further highlighting the utility of ALL-tRNAseq for translational research.

Single-read tRNA-seq analysis reveals coordination of tRNA modification and aminoacylation and fragmentation

Transfer RNA (tRNA) utilizes multiple properties of abundance, modification, and aminoacylation in translational regulation. These properties were typically studied one-by-one; however, recent advance in high throughput tRNA sequencing enables their simultaneous assessment in the same sequencing data. How these properties are coordinated at the transcriptome level is an open question. Here, we develop a single-read tRNA analysis pipeline that takes advantage of the pseudo single-molecule nature of tRNA sequencing in NGS libraries. tRNAs are short enough that a single NGS read can represent one tRNA molecule, and can simultaneously report on the status of multiple modifications, aminoacylation, and fragmentation of each molecule. We find correlations among modification-modification, modification-aminoacylation and modification-fragmentation. We identify interdependencies among one of the most common tRNA modifications, m1A58, as coordinators of tissue-specific gene expression. Our method, SingLe-read Analysis of Crosstalks (SLAC), reveals tRNAome-wide networks of modifications, aminoacylation, and fragmentation. We observe changes of these networks under different stresses, and assign a function for tRNA modification in translational regulation and fragment biogenesis. SLAC leverages the richness of the tRNA-seq data and provides new insights on the coordination of tRNA properties.