tRNA base methylation identification and quantification via high-throughput sequencing

Human DUS1L catalyzes dihydrouridine modification at tRNA positions 16/17, and DUS1L overexpression perturbs translation

Human cytoplasmic tRNAs contain dihydrouridine modifications at positions 16 and 17 (D16/D17). The enzyme responsible for D16/D17 formation and its cellular roles remain elusive. Here, we identify DUS1L as the human tRNA D16/D17 writer. DUS1L knockout in the glioblastoma cell lines LNZ308 and U87 causes loss of D16/D17. D formation is reconstituted in vitro using recombinant DUS1L in the presence of NADPH or NADH. DUS1L knockout/overexpression in LNZ308 cells shows that DUS1L supports cell growth. Moreover, higher DUS1L expression in glioma patients is associated with poorer prognosis. Upon vector-mediated DUS1L overexpression in LNZ308 cells, 5' and 3' processing of precursor tRNATyr(GUA) is inhibited, resulting in a reduced mature tRNATyr(GUA) level, reduced translation of the tyrosine codons UAC and UAU, and reduced translational readthrough of the near-cognate stop codons UAA and UAG. Moreover, DUS1L overexpression increases the amounts of several D16/D17-containing tRNAs and total cellular translation. Our study identifies a human dihydrouridine writer, providing the foundation to study its roles in health and disease.

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.

TRMT10A dysfunction perturbs codon translation of initiator methionine and glutamine and impairs brain functions in mice

In higher eukaryotes, tRNA methyltransferase 10A (TRMT10A) is responsible for N1-methylguanosine modification at position nine of various cytoplasmic tRNAs. Pathogenic mutations in TRMT10A cause intellectual disability, microcephaly, diabetes, and short stature in humans, and generate cytotoxic tRNA fragments in cultured cells; however, it is not clear how TRMT10A supports codon translation or brain functions. Here, we generated Trmt10a null mice and showed that tRNAGln(CUG) and initiator methionine tRNA levels were universally decreased in various tissues; the same was true in a human cell line lacking TRMT10A. Ribosome profiling of mouse brain revealed that dysfunction of TRMT10A causes ribosome slowdown at the Gln(CAG) codon and increases translation of Atf4 due to higher frequency of leaky scanning of its upstream open reading frames. Broadly speaking, translation of a subset of mRNAs, especially those for neuronal structures, is perturbed in the mutant brain. Despite not showing discernable defects in the pancreas, liver, or kidney, Trmt10a null mice showed lower body weight and smaller hippocampal postsynaptic densities, which is associated with defective synaptic plasticity and memory. Taken together, our study provides mechanistic insight into the roles of TRMT10A in the brain, and exemplifies the importance of universal tRNA modification during translation of specific codons.

tRNA-m1A methylation controls the infection of Magnaporthe oryzae by supporting ergosterol biosynthesis

Ergosterols are essential components of fungal plasma membranes. Inhibitors targeting ergosterol biosynthesis (ERG) genes are critical for controlling fungal pathogens, including Magnaporthe oryzae, the fungus that causes rice blast. However, the translational mechanisms governing ERG gene expression remain largely unexplored. Here, we show that the Trm6/Trm61 complex catalyzes dynamic N1-methyladenosine at position 58 (m1A58) in 51 transfer RNAs (tRNAs) of M. oryzae, significantly influencing translation at both the initiation and elongation stages. Notably, tRNA m1A58 promotes elongation speed at most cognate codons mainly by enhancing eEF1-tRNA binding rather than affecting tRNA abundance or charging. The absence of m1A58 leads to substantial decreases in the translation of ERG genes, ergosterol production, and, consequently, fungal virulence. Simultaneously targeting the Trm6/Trm61 complex and the ergosterol biosynthesis pathway markedly improves rice blast control. Our findings demonstrate an important role of m1A58-mediated translational regulation in ergosterol production and fungal infection, offering a potential strategy for fungicide development.

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 transfer RNA sequencing (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.

m3C32 tRNA modification controls serine codon-biased mRNA translation, cell cycle, and DNA-damage response

The epitranscriptome includes a diversity of RNA modifications that influence gene expression. N3-methylcytidine (m3C) mainly occurs in the anticodon loop (position C32) of certain tRNAs yet its role is poorly understood. Here, using HAC-Seq, we report comprehensive METTL2A/2B-, METTL6-, and METTL2A/2B/6-dependent m3C profiles in human cells. METTL2A/2B modifies tRNA-arginine and tRNA-threonine members, whereas METTL6 modifies the tRNA-serine family. However, decreased m3C32 on tRNA-Ser-GCT isodecoders is only observed with combined METTL2A/2B/6 deletion. Ribo-Seq reveals altered translation of genes related to cell cycle and DNA repair pathways in METTL2A/2B/6-deficient cells, and these mRNAs are enriched in AGU codons that require tRNA-Ser-GCT for translation. These results, supported by reporter assays, help explain the observed altered cell cycle, slowed proliferation, and increased cisplatin sensitivity phenotypes of METTL2A/2B/6-deficient cells. Thus, we define METTL2A/2B/6-dependent methylomes and uncover a particular requirement of m3C32 tRNA modification for serine codon-biased mRNA translation of cell cycle, and DNA repair genes.

The 'Not-So-Famous Five' in tumorigenesis: tRNAs, tRNA fragments, and tRNA epitranscriptome in concert with AARSs and AIMPs

RNA profiling studies have revealed that ∼75% of the human genome is transcribed to RNA but only a meagre fraction of it is translated to proteins. Majority of transcribed RNA constitute a specialized pool of non-coding RNAs. Human genome contains approximately 506 genes encoding a set of 51 different tRNAs, constituting a unique class of non-coding RNAs that not only have essential housekeeping functions as translator molecules during protein synthesis, but have numerous uncharted regulatory functions. Intriguing findings regarding a variety of non-canonical functions of tRNAs, tRNA derived fragments (tRFs), esoteric epitranscriptomic modifications of tRNAs, along with aminoacyl-tRNA synthetases (AARSs) and ARS-interacting multifunctional proteins (AIMPs), envision a 'peripheral dogma' controlling the flow of genetic information in the backdrop of qualitative information wrung out of the long-live central dogma of molecular biology, to drive cells towards either proliferation or differentiation programs. Our review will substantiate intriguing peculiarities of tRNA gene clusters, atypical tRNA-transcription from internal promoters catalysed by another distinct RNA polymerase enzyme, dynamically diverse tRNA epitranscriptome, intricate mechanism of tRNA-charging by AARSs governing translation fidelity, epigenetic regulation of gene expression by tRNA fragments, and the role of tRNAs and tRNA derived/associated molecules as quantitative determinants of the functional proteome, covertly orchestrating the process of tumorigenesis, through a deregulated tRNA-ome mediating selective codon-biased translation of cancer related gene transcripts.

Internal mechanism of correlation between angiotensin II gene and serum adiponectin level in patients with cerebrovascular complications of H-type hypertension

Background

The study aimed to explore the correlation between the angiotensin II (Ang II) gene and serum adiponectin expression in patients with cerebrovascular complications of H-type hypertension (HH) and its mechanism.

Methods

A total of 50 cases of outpatient patients in Tianjin Fourth Central Hospital were recruited from January 2022 to June 2023 and rolled into three groups according to their blood pressure and basic information, namely the HH cerebrovascular complications group, the non-H-type hypertension (NHH) group, and the healthy control (HC) group. Peripheral blood samples were taken; one sample was utilized to test for the Ang II gene and the methylation of Ang II, and the other sample was utilized to measure serum adiponectin levels to analyze the relationship between serum adiponectin level and Ang II in patients with cerebrovascular complications of HH.

Results

The ratio of male to female was 8:7 in the group of cerebrovascular complications of HH, and mean systolic blood pressure (SBP) and diastolic blood pressure (DBP) were 167.34 mm Hg and 112.56 mm Hg, respectively. In the NHH group, the mean SBP was 165.89 mm Hg, and the mean DBP was 113.47 mm Hg. The blood pressure of the HC group was in the normal range. The Ang II content was the highest in the group with cerebrovascular complications of HH, followed by the group with NHH, and the lowest in the HC group.

Conclusions

Pyrosequencing chart of patients with cerebrovascular complications of HH showed that the content of deoxyphosphate ribose G was the highest, while the content of A was the highest in NHH patients. Moreover, the serum adiponectin level of patients with HH and NHH was superior to that of the HC group, and the adiponectin level between the former two groups and the HC group differed considerably. Ang II levels were high in patients with cerebrovascular complications of HH and were positively correlated with adiponectin levels. The incidence of cerebrovascular complications of HH may be related to Ang II levels in patients.

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.

Quantification of tRNA m1A modification by templated-ligation qPCR

N1-methyladenosine (m1A) is a widespread modification in all eukaryotic, many archaeal, and some bacterial tRNAs. m1A is generally located in the T loop of cytosolic tRNA and between the acceptor and D stems of mitochondrial tRNAs; it is involved in the tertiary interaction that stabilizes tRNA. Human tRNA m1A levels are dynamically regulated that fine-tune translation and can also serve as biomarkers for infectious disease. Although many methods have been used to measure m1A, a PCR method to assess m1A levels quantitatively in specific tRNAs has been lacking. Here we develop a templated-ligation followed by a qPCR method (TL-qPCR) that measures m1A levels in target tRNAs. Our method uses the SplintR ligase that efficiently ligates two tRNA complementary DNA oligonucleotides using tRNA as the template, followed by qPCR using the ligation product as the template. m1A interferes with the ligation in specific ways, allowing for the quantitative assessment of m1A levels using subnanogram amounts of total RNA. We identify the features of specificity and quantitation for m1A-modified model RNAs and apply these to total RNA samples from human cells. Our method enables easy access to study the dynamics and function of this pervasive tRNA modification.

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.

Untacking small RNA profiling and RNA fragment footprinting: Approaches and challenges in library construction

Small RNAs (sRNAs) with sizes ranging from 15 to 50 nucleotides (nt) are critical regulators of gene expression control. Prior studies have shown that sRNAs are involved in a broad range of biological processes, such as organ development, tumorigenesis, and epigenomic regulation; however, emerging evidence unveils a hidden layer of diversity and complexity of endogenously encoded sRNAs profile in eukaryotic organisms, including novel types of sRNAs and the previously unknown post-transcriptional RNA modifications. This underscores the importance for accurate, unbiased detection of sRNAs in various cellular contexts. A multitude of high-throughput methods based on next-generation sequencing (NGS) are developed to decipher the sRNA expression and their modifications. Nonetheless, distinct from mRNA sequencing, the data from sRNA sequencing suffer frequent inconsistencies and high variations emanating from the adapter contaminations and RNA modifications, which overall skew the sRNA libraries. Here, we summarize the sRNA-sequencing approaches, and discuss the considerations and challenges for the strategies and methods of sRNA library construction. The pros and cons of sRNA sequencing have significant implications for implementing RNA fragment footprinting approaches, including CLIP-seq and Ribo-seq. We envision that this review can inspire novel improvements in small RNA sequencing and RNA fragment footprinting in future. This article is categorized under: RNA Evolution and Genomics > Computational Analyses of RNA RNA Processing > Processing of Small RNAs Regulatory RNAs/RNAi/Riboswitches > Biogenesis of Effector Small RNAs.

Roles and regulation of tRNA-derived small RNAs in animals

A growing class of small RNAs, known as tRNA-derived RNAs (tdRs), tRNA-derived small RNAs or tRNA-derived fragments, have long been considered mere intermediates of tRNA degradation. These small RNAs have recently been implicated in an evolutionarily conserved repertoire of biological processes. In this Review, we discuss the biogenesis and molecular functions of tdRs in mammals, including tdR-mediated gene regulation in cell metabolism, immune responses, transgenerational inheritance, development and cancer. We also discuss the accumulation of tRNA-derived stress-induced RNAs as a distinct adaptive cellular response to pathophysiological conditions. Furthermore, we highlight new conceptual advances linking RNA modifications with tdR activities and discuss challenges in studying tdR biology in health and disease.

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.

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.

Discovery, structure, mechanisms, and evolution of protein-only RNase P enzymes

The endoribonuclease RNase P is responsible for tRNA 5' maturation in all domains of life. A unique feature of RNase P is the variety of enzyme architectures, ranging from dual- to multi-subunit ribonucleoprotein forms with catalytic RNA subunits to protein-only enzymes, the latter occurring as single- or multi-subunit forms or homo-oligomeric assemblies. The protein-only enzymes evolved twice: a eukaryal protein-only RNase P termed PRORP and a bacterial/archaeal variant termed homolog of Aquifex RNase P (HARP); the latter replaced the RNA-based enzyme in a small group of thermophilic bacteria but otherwise coexists with the ribonucleoprotein enzyme in a few other bacteria as well as in those archaea that also encode a HARP. Here we summarize the history of the discovery of protein-only RNase P enzymes and review the state of knowledge on structure and function of bacterial HARPs and eukaryal PRORPs, including human mitochondrial RNase P as a paradigm of multi-subunit PRORPs. We also describe the phylogenetic distribution and evolution of PRORPs, as well as possible reasons for the spread of PRORPs in the eukaryal tree and for the recruitment of two additional protein subunits to metazoan mitochondrial PRORP. We outline potential applications of PRORPs in plant biotechnology and address diseases associated with mutations in human mitochondrial RNase P genes. Finally, we consider possible causes underlying the displacement of the ancient RNA enzyme by a protein-only enzyme in a small group of bacteria.

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.

The Molecular Basis of Human ALKBH3 Mediated RNA N1 -methyladenosine (m1 A) Demethylation

N1 -methyladenosine (m1 A) is a prevalent post-transcriptional RNA modification, and the distribution and dynamics of the modification play key epitranscriptomic roles in cell development. At present, the human AlkB Fe(II)/α-ketoglutarate-dependent dioxygenase family member ALKBH3 is the only known mRNA m1 A demethylase, but its catalytic mechanism remains unclear. Here, we present the structures of ALKBH3-oligo crosslinked complexes obtained with the assistance of a synthetic antibody crystallization chaperone. Structural and biochemical results showed that ALKBH3 utilized two β-hairpins (β4-loop-β5 and β'-loop-β'') and the α2 helix to facilitate single-stranded substrate binding. Moreover, a bubble-like region around Asp194 and a key residue inside the active pocket (Thr133) enabled specific recognition and demethylation of m1 A- and 3-methylcytidine (m3 C)-modified substrates. Mutation of Thr133 to the corresponding residue in the AlkB Fe(II)/α-ketoglutarate-dependent dioxygenase family members FTO or ALKBH5 converted ALKBH3 substrate selectivity from m1 A to N6 -methyladenosine (m6 A), as did Asp194 deletion. Our findings provide a molecular basis for understanding the mechanisms of substrate recognition and m1 A demethylation by ALKBH3. This study is expected to aid structure-guided design of chemical probes for further functional studies and therapeutic applications.

Advances in the Structural and Functional Understanding of m1A RNA Modification

ConspectusRNA modification is a co- or post-transcriptional process by which specific nucleotides are chemically altered by enzymes after their initial incorporation into the RNA chain, expanding the chemical and functional diversity of RNAs. Our understanding of RNA modifications has changed dramatically in recent years. In the past decade, RNA methyltransferases (MTases) have been highlighted in numerous clinical studies and disease models, modifications have been found to be dynamically regulated by demodification enzymes, and significant technological advances have been made in the fields of RNA sequencing, mass spectrometry, and structural biology. Among RNAs, transfer RNAs (tRNAs) exhibit the greatest diversity and density of post-transcriptional modifications, which allow for potential cross-talks and regulation during their incorporation. N1-methyladenosine (m1A) modification is found in tRNAs at positions 9, 14, 16, 22, 57, and 58, depending on the tRNA and organism.Our laboratory has used and developed a large panel of tools to decipher the different mechanisms used by m1A tRNA MTases to recognize and methylate tRNA. We have solved the structures of TrmI from Thermus thermophilus (m1A58), TrmK from Bacillus subtilis (m1A22), and human TRMT10C (m1A9). These MTases do not share the same structure or organization to recognize tRNAs, but they all modify an adenosine, forming a non-Watson-Crick (WC) interaction. For TrmK, nuclear magnetic resonance (NMR) chemical shift mapping of the binding interface between TrmK and tRNASer was invaluable to build a TrmK/tRNA model, where both domains of TrmK participate in the binding of a full-length L-shaped tRNA and where the non-WC purine 13-A22 base pair positions the A22 N1-atom close to the methyl of the S-adenosyl-l-methionine (SAM) TrmK cofactor. For TRMT10C, cryoEM structures showed the MTase poised to N1-methylate A9 or G9 in tRNA and revealed different steps of tRNA maturation, where TRMT10C acts as a tRNA binding platform for sequential docking of each maturation enzyme. This work confers a role for TRMT10C in tRNA quality control and provides a framework to understand the link between mitochondrial tRNA maturation dysfunction and diseases.Methods to directly detect the incorporation of modifications during tRNA biosynthesis are rare and do not provide easy access to the temporality of their introduction. To this end, we have introduced time-resolved NMR to monitor tRNA maturation in the cellular environment. Combined with genetic and biochemical approaches involving the synthesis of specifically modified tRNAs, our methodology revealed that some modifications are incorporated in a defined sequential order, controlled by cross-talks between modification events. In particular, a strong modification circuit, namely Ψ55 → m5U54 → m1A58, controls the modification process in the T-arm of yeast elongator tRNAs. Conversely, we showed that m1A58 is efficiently introduced on unmodified initiator tRNAiMet without the need of any prior modification. Two distinct pathways are therefore followed for m1A58 incorporation in elongator and initiator tRNAs.We are undoubtedly entering an exciting period for the elucidation of the functions of RNA modifications and the intricate mechanisms by which modification enzymes identify and alter their RNA substrates. These are promising directions for the field of epitranscriptomics.

Engineered mischarged transfer RNAs for correcting pathogenic missense mutations

Missense mutations account for approximately 50% of pathogenic mutations in human genetic diseases, and most lack effective treatments. Gene therapies, gene editing, and RNA therapies, including transfer RNA (tRNA) modalities, are common strategies for genetic disease treatments. However, reported tRNA therapies are for nonsense mutations only. It has not been explored how tRNAs can be engineered to correct missense mutations. Here, we describe missense-correcting tRNAs (mc-tRNAs) as a potential therapeutic for correcting pathogenic missense mutations. Mc-tRNAs are engineered tRNAs charged with one amino acid, but read codons of another in translation. We first developed a series of fluorescent protein-based reporters that indicate the successful correction of missense mutations via restoration of fluorescence. We engineered mc-tRNAs that effectively corrected serine and arginine missense mutations in the reporters and confirmed the amino acid substitution by mass spectrometry and mc-tRNA expression by sequencing. We examined the transcriptome response to mc-tRNA expression and found some mc-tRNAs induced minimum transcriptomic changes. Furthermore, we applied an mc-tRNA to rescue a pathogenic CAPN3 Arg-to-Gln mutant involved in LGMD2A. These results establish a versatile pipeline for mc-tRNA engineering and demonstrate the potential of mc-tRNA as an alternative therapeutic platform for the treatment of genetic disorders.

General Principles and Limitations for Detection of RNA Modifications by Sequencing

Among the many analytical methods applied to RNA modifications, a particularly pronounced surge has occurred in the past decade in the field of modification mapping. The occurrence of modifications such as m6A in mRNA, albeit known since the 1980s, became amenable to transcriptome-wide analyses through the advent of next-generation sequencing techniques in a rather sudden manner. The term "mapping" here refers to detection of RNA modifications in a sequence context, which has a dramatic impact on the interpretation of biological functions. As a consequence, an impressive number of mapping techniques were published, most in the perspective of what now has become known as "epitranscriptomics". While more and more different modifications were reported to occur in mRNA, conflicting reports and controversial results pointed to a number of technical and theoretical problems rooted in analytics, statistics, and reagents. Rather than finding the proverbial needle in a haystack, the tasks were to determine how many needles of what color in what size of a haystack one was looking at.As the authors of this Account, we think it important to outline the limitations of different mapping methods since many life scientists freshly entering the field confuse the accuracy and precision of modification mapping with that of normal sequencing, which already features numerous caveats by itself. Indeed, we propose here to qualify a specific mapping method by the size of the transcriptome that can be meaningfully analyzed with it.We here focus on high throughput sequencing by Illumina technology, referred to as RNA-Seq. We noted with interest the development of methods for modification detection by other high throughput sequencing platforms that act directly on RNA, e.g., PacBio SMRT and nanopore sequencing, but those are not considered here.In contrast to approaches relying on direct RNA sequencing, current Illumina RNA-Seq protocols require prior conversion of RNA into DNA. This conversion relies on reverse transcription (RT) to create cDNA; thereafter, the cDNA undergoes a sequencing-by-synthesis type of analysis. Thus, a particular behavior of RNA modified nucleotides during the RT-step is a prerequisite for their detection (and quantification) by deep sequencing, and RT properties have great influence on the detection efficiency and reliability. Moreover, the RT-step requires annealing of a synthetic primer, a prerequisite with a crucial impact on library preparation. Thus, all RNA-Seq protocols must feature steps for the introduction of primers, primer landing sites, or adapters on both the RNA 3'- and 5'-ends.

Mistranslating the genetic code with leucine in yeast and mammalian cells

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

Diversity in Biological Function and Mechanism of the tRNA Methyltransferase Trm10

Transfer ribonucleic acid (tRNA) is the most highly modified RNA species in the cell, and loss of tRNA modifications can lead to growth defects in yeast as well as metabolic, neurological, and mitochondrial disorders in humans. Significant progress has been made toward identifying the enzymes that are responsible for installing diverse modifications in tRNA, revealing a landscape of fascinating biological and mechanistic diversity that remains to be fully explored. Most early discoveries of tRNA modification enzymes were in model systems, where many enzymes were not strictly required for viability, an observation somewhat at odds with the extreme conservation of many of the same enzymes throughout multiple domains of life. Moreover, many tRNA modification enzymes act on more than one type of tRNA substrate, which is not necessarily surprising given the similar overall secondary and tertiary structures of tRNA, yet biochemical characterization has revealed interesting patterns of substrate specificity that can be challenging to rationalize on a molecular level. Questions about how many enzymes efficiently select a precise set of target tRNAs from among a structurally similar pool of molecules persist.The tRNA methyltransferase Trm10 provides an exciting paradigm to study the biological and mechanistic questions surrounding tRNA modifications. Even though the enzyme was originally characterized in Saccharomyces cerevisiae where its deletion causes no detectable phenotype under standard lab conditions, several more recently identified phenotypes provide insight into the requirement for this modification in the overall quality control of the tRNA pool. Studies of Trm10 in yeast also revealed another characteristic feature that has turned out to be a conserved feature of enzymes throughout the Trm10 family tree. We were initially surprised to see that purified S. cerevisiae Trm10 was capable of modifying tRNA substrates that were not detectably modified by the enzyme in vivo in yeast. This pattern has continued to emerge as we and others have studied Trm10 orthologs from Archaea and Eukarya, with enzymes exhibiting in vitro substrate specificities that can differ significantly from in vivo patterns of modification. While this feature complicates efforts to predict substrate specificities of Trm10 enzymes in the absence of appropriate genetic systems, it also provides an exciting opportunity for studying how enzyme activities can be regulated to achieve dynamic patterns of biological tRNA modification, which have been shown to be increasingly important for stress responses and human disease. Finally, the intriguing diversity in target nucleotide modification that has been revealed among Trm10 orthologs is distinctive among known tRNA modifying enzymes and necessitates unusual and likely novel catalytic strategies for methylation that are being revealed by biochemical and structural studies directed toward various family members. These efforts will no doubt yield more surprising discoveries in terms of tRNA modification enzymology.

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.

Expression profiles and function prediction of tRNA-derived fragments in glioma

BACKGROUND

Glioblastoma (GBM) is the most aggressive malignant primary brain tumor. The transfer RNA-derived fragments (tRFs) are a new group of small noncoding RNAs, which are dysregulated in many cancers. Until now, the expression and function of tRFs in glioma remain unknown.

METHODS

The expression profiles of tRF subtypes were analyzed using the Cancer Genome Atlas (TCGA)-low-grade gliomas (LGG)/GBM dataset. The target genes of tRFs were subjected to Gene Ontology, Kyoto Encyclopedia and Gene set enrichment analysis of Genes and Genomes pathway enrichment analysis. The protein-protein interaction enrichment analysis was performed by STRING. QRT-PCR was performed to detect the expressions of tRFs in human glioma cell lines U87, U373, U251, and human astrocyte cell line SVG p12. Western blot assay was used to detect to the expression of S100A11. The interaction between tRF-19-R118LOJX and S100A11 mRNA 3'UTR was detected by dual-luciferase reporter assay. The effects of tRF-19-R118LOJX, tRF-19-6SM83OJX and S100A11 on the glioma cell proliferation, migration and in vitro vasculogenic mimicry formation ability were examined by CCK-8 proliferation assay, EdU assay, HoloMonitor cell migration assay and tube formation assay, respectively.

RESULTS

tRF-19-R118LOJX and tRF-19-6SM83OJX are the most differentially expressed tRFs between LGG and GBM groups. The functional enrichment analysis showed that the target genes of tRF-19-R118LOJX and tRF-19-6SM83OJX are enriched in regulating blood vessel development. The upregulated target genes are linked to adverse survival outcomes in glioma patients. tRF-19-R118LOJX and tRF-19-6SM83OJX were identified to suppress glioma cell proliferation, migration, and in vitro vasculogenic mimicry formation. The mechanism of tRF-19-R118LOJX might be related to its function as an RNA silencer by targeting the S100A11 mRNA 3'UTR.

CONCLUSION

tRFs would become novel diagnostic biomarkers and therapeutic targets of glioma, and the mechanism might be related to its post-transcriptionally regulation of gene expression by targeting mRNA 3'UTR.

Mitochondrial RNA m3C methyltransferase METTL8 relies on an isoform-specific N-terminal extension and modifies multiple heterogenous tRNAs

Methyltransferase-like 8 (METTL8) encodes a mitochondria-localized METTL8-Iso1 and a nucleolus-distributed METTL8-Iso4 isoform, which differ only in their N-terminal extension (N-extension), by mRNA alternative splicing. METTL8-Iso1 generates 3-methylcytidine at position 32 (m3C32) of mitochondrial tRNAThr and tRNASer(UCN). Whether METTL8-Iso4 is an active m3C32 methyltransferase and the role of the N-extension in mitochondrial tRNA m3C32 formation remain unclear. Here, we revealed that METTL8-Iso4 was inactive in m3C32 generation due to the lack of N-extension, which contains several absolutely conserved modification-critical residues; the counterparts were likewise essential in cytoplasmic m3C32 biogenesis by methyltransferase-like 2A (METTL2A) or budding yeasts tRNA N3-methylcytidine methyltransferase (Trm140), in vitro and in vivo. Cross-compartment/species tRNA modification assays unexpectedly found that METTL8-Iso1 efficiently introduced m3C32 to several cytoplasmic or even bacterial tRNAs in vitro. m3C32 did not influence tRNAThrN6-threonylcarbamoyladenosine (t6A) modification or aminoacylation. In addition to its interaction with mitochondrial seryl-tRNA synthetase (SARS2), we further discovered an interaction between mitochondrial threonyl-tRNA synthetase (TARS2) and METTL8-Iso1. METTL8-Iso1 substantially stimulated the aminoacylation activities of SARS2 and TARS2 in vitro, suggesting a functional connection between mitochondrial tRNA modification and charging. Altogether, our results deepen the mechanistic insights into mitochondrial m3C32 biogenesis and provide a valuable route to prepare cytoplasmic/bacterial tRNAs with only a m3C32 moiety, aiding in future efforts to investigate its effects on tRNA structure and function.

Cysteine tRNA acts as a stop codon readthrough-inducing tRNA in the human HEK293T cell line

Under certain circumstances, any of the three termination codons can be read through by a near-cognate tRNA; i.e., a tRNA whose two out of three anticodon nucleotides base pair with those of the stop codon. Unless programed to synthetize C-terminally extended protein variants with expanded physiological roles, readthrough represents an undesirable translational error. On the other side of a coin, a significant number of human genetic diseases is associated with the introduction of nonsense mutations (premature termination codons [PTCs]) into coding sequences, where stopping is not desirable. Here, the tRNA's ability to induce readthrough opens up the intriguing possibility of mitigating the deleterious effects of PTCs on human health. In yeast, the UGA and UAR stop codons were described to be read through by four readthrough-inducing rti-tRNAs-tRNATrp and tRNACys, and tRNATyr and tRNAGln, respectively. The readthrough-inducing potential of tRNATrp and tRNATyr was also observed in human cell lines. Here, we investigated the readthrough-inducing potential of human tRNACys in the HEK293T cell line. The tRNACys family consists of two isoacceptors, one with ACA and the other with GCA anticodons. We selected nine representative tRNACys isodecoders (differing in primary sequence and expression level) and tested them using dual luciferase reporter assays. We found that at least two tRNACys can significantly elevate UGA readthrough when overexpressed. This indicates a mechanistically conserved nature of rti-tRNAs between yeast and human, supporting the idea that they could be used in the PTC-associated RNA therapies.

The life and times of a tRNA

The study of eukaryotic tRNA processing has given rise to an explosion of new information and insights in the last several years. We now have unprecedented knowledge of each step in the tRNA processing pathway, revealing unexpected twists in biochemical pathways, multiple new connections with regulatory pathways, and numerous biological effects of defects in processing steps that have profound consequences throughout eukaryotes, leading to growth phenotypes in the yeast Saccharomyces cerevisiae and to neurological and other disorders in humans. This review highlights seminal new results within the pathways that comprise the life of a tRNA, from its birth after transcription until its death by decay. We focus on new findings and revelations in each step of the pathway including the end-processing and splicing steps, many of the numerous modifications throughout the main body and anticodon loop of tRNA that are so crucial for tRNA function, the intricate tRNA trafficking pathways, and the quality control decay pathways, as well as the biogenesis and biology of tRNA-derived fragments. We also describe the many interactions of these pathways with signaling and other pathways in the cell.

Small RNA modifications: regulatory molecules and potential applications

Small RNAs (also referred to as small noncoding RNAs, sncRNA) are defined as polymeric ribonucleic acid molecules that are less than 200 nucleotides in length and serve a variety of essential functions within cells. Small RNA species include microRNA (miRNA), PIWI-interacting RNA (piRNA), small interfering RNA (siRNA), tRNA-derived small RNA (tsRNA), etc. Current evidence suggest that small RNAs can also have diverse modifications to their nucleotide composition that affect their stability as well as their capacity for nuclear export, and these modifications are relevant to their capacity to drive molecular signaling processes relevant to biogenesis, cell proliferation and differentiation. In this review, we highlight the molecular characteristics and cellular functions of small RNA and their modifications, as well as current techniques for their reliable detection. We also discuss how small RNA modifications may be relevant to the clinical applications for the diagnosis and treatment of human health conditions such as cancer.

The covalent nucleotide modifications within plant mRNAs: What we know, how we find them, and what should be done in the future

Although covalent nucleotide modifications were first identified on the bases of transfer RNAs (tRNAs) and ribosomal RNAs (rRNAs), a number of these epitranscriptome marks have also been found to occur on the bases of messenger RNAs (mRNAs). These covalent mRNA features have been demonstrated to have various and significant effects on the processing (e.g. splicing, polyadenylation, etc.) and functionality (e.g. translation, transport, etc.) of these protein-encoding molecules. Here, we focus our attention on the current understanding of the collection of covalent nucleotide modifications known to occur on mRNAs in plants, how they are detected and studied, and the most outstanding future questions of each of these important epitranscriptomic regulatory signals.

Bibliometric and visual analysis of RAN methylation in cardiovascular disease

Background

RNA methylation is associated with cardiovascular disease (CVD) occurrence and development. The purpose of this study is to visually analyze the results and research trends of global RNA methylation in CVD.

Methods

Articles and reviews on RNA methylation in CVD published before 6 November 2022 were searched in the Web of Science Core Collection. Visual and statistical analysis was performed using CiteSpace 1.6.R4 advanced and VOSviewer 1.6.18.

Results

There were 847 papers from 1,188 institutions and 63 countries/regions. Over approximately 30 years, there was a gradual increase in publications and citations on RNA methylation in CVD. America and China had the highest output (284 and 259 papers, respectively). Nine of the top 20 institutions that published articles were from China, among which Fudan University represented the most. The International Journal of Molecular Sciences was the journal with the most studies. Nature was the most co-cited journal. The most influential writers were Zhang and Wang from China and Mathiyalagan from the United States. After 2015, the primary keywords were cardiac development, heart, promoter methylation, RNA methylation, and N6-methyladenosine. Nuclear RNA, m6A methylation, inhibition, and myocardial infarction were the most common burst keywords from 2020 to the present.

Conclusions

A bibliometric analysis reveals research hotspots and trends of RNA methylation in CVD. The regulatory mechanisms of RNA methylation related to CVD and the clinical application of their results, especially m6A methylation, are likely to be the focus of future research.

Reverse Transcriptases: From Discovery and Applications to Xenobiology

Reverse transcriptases are DNA polymerases that can use RNA as a template for DNA synthesis. They thus catalyze the reverse of transcription. Although discovered in 1970, reverse transcriptases are still of great interest and are constantly being further developed for numerous modern research approaches. They are frequently used in biotechnological and molecular diagnostic applications. In this review, we describe the discovery of these fascinating enzymes and summarize research results and applications ranging from molecular cloning, direct virus detection, and modern sequencing methods to xenobiology.

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.

tRNAs Are Stable After All: Pitfalls in Quantification of tRNA from Starved Escherichia coli Cultures Exposed by Validation of RNA Purification Methods

tRNAs and ribosomal RNAs are often considered stable RNAs. In contrast to this view, we recently proposed that tRNAs are degraded during amino acid starvation and drug-induced transcription inhibition. However, reevaluation of our experimental approach revealed that common RNA extraction methods suffer from alarming extraction and size biases that can lead to gross underestimation of RNA levels in starved Escherichia coli populations. Quantification of tRNAs suffers additional biases due to differing fractions of tRNAs with base modifications in growing versus starved bacteria. Applying an improved methodology, we measured tRNA levels after starvation for amino acids, glucose, phosphate, or ammonium and transcription inhibition by rifampicin. We report that tRNA levels remain largely unaffected in all tested conditions, including several days of starvation. This confirms that tRNAs are remarkably stable RNAs and serves as a cautionary tale about quantification of RNA from cells cultured outside the steady-state growth regime. rRNA, conversely, is extensively degraded during starvation. Thus, E. coli downregulates the translation machinery in response to starvation by reducing the ribosome pool through rRNA degradation, while a high concentration of tRNAs available to supply amino acids to the remaining ribosomes is maintained. IMPORTANCE We show that E. coli tRNAs are remarkably stable during several days of nutrient starvation, although rRNA is degraded extensively under these conditions. The levels of these two major RNA classes are considered to be strongly coregulated at the level of transcription. We demonstrate that E. coli can control the ratio of tRNAs per ribosome under starvation by means of differential degradation rates. The question of tRNA stability in stressed E. coli cells has become subject to debate. Our in-depth analysis of RNA quantification methods reveals hidden technical pitfalls at every step of the analysis, from RNA extraction to target detection and normalization. Most importantly, starved E. coli populations were more resilient to RNA extraction than unstarved populations. The current results underscore that the seemingly trivial task of quantifying an abundant RNA species is not straightforward for cells cultured outside the exponential growth regime.

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.

CoLoC-seq probes the global topology of organelle transcriptomes

Proper RNA localisation is essential for physiological gene expression. Various kinds of genome-wide approaches permit to comprehensively profile subcellular transcriptomes. Among them, cell fractionation methods, that couple RNase treatment of isolated organelles to the sequencing of protected transcripts, remain most widely used, mainly because they do not require genetic modification of the studied system and can be easily implemented in any cells or tissues, including in non-model species. However, they suffer from numerous false-positives since incompletely digested contaminant RNAs can still be captured and erroneously identified as resident transcripts. Here we introduce Controlled Level of Contamination coupled to deep sequencing (CoLoC-seq) as a new subcellular transcriptomics approach that efficiently bypasses this caveat. CoLoC-seq leverages classical enzymatic kinetics and tracks the depletion dynamics of transcripts in a gradient of an exogenously added RNase, with or without organellar membranes. By means of straightforward mathematical modelling, CoLoC-seq infers the localisation topology of RNAs and robustly distinguishes between genuinely resident, luminal transcripts and merely abundant surface-attached contaminants. Our generic approach performed well on human mitochondria and is in principle applicable to other membrane-bounded organelles, including plastids, compartments of the vacuolar system, extracellular vesicles, and viral particles.

N1-methyladenosine modification in cancer biology: current status and future perspectives

Post-transcriptional modifications in RNAs regulate their biological behaviors and functions. N¹-methyladenosine (m¹A), which is dynamically regulated by writers, erasers and readers, has been found as a reversible modification in tRNA, mRNA, rRNA and long non-coding RNA (lncRNA). m¹A modification has impacts on the RNA processing, structure and functions of targets. Increasing studies reveal the critical roles of m¹A modification and its regulators in tumorigenesis. Due to the positive relevance between m¹A and cancer development, targeting m¹A modification and m¹A-related regulators has been of attention. In this review, we summarized the current understanding of m¹A in RNAs, covering the modulation of m¹A modification in cancer biology, as well as the possibility of targeting m¹A modification as a potential target for cancer diagnosis and therapy.

tRNA abundance, modification and fragmentation in nasopharyngeal swabs as biomarkers for COVID-19 severity

Emerging and re-emerging respiratory viruses can spread rapidly and cause pandemics as demonstrated by the recent severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) pandemic. The early human immune responses to respiratory viruses are in the nasal cavity and nasopharyngeal regions. Defining biomarkers of disease trajectory at the time of a positive diagnostic test would be an important tool to facilitate decisions such as initiation of antiviral treatment. We hypothesize that nasopharyngeal tRNA profiles could be used to predict Coronavirus Disease 19 (COVID-19) severity. We carried out multiplex small RNA sequencing (MSR-seq) on residual nasopharyngeal swabs to measure simultaneously full-length tRNA abundance, tRNA modifications, and tRNA fragmentation for the human tRNA response to SARS-CoV-2 infection. We identified distinct tRNA signatures associated with mild symptoms versus severe COVID-19 manifestations requiring hospitalization. These results highlight the utility of host tRNA properties as biomarkers for the clinical outcome of SARS-CoV-2.

Detection technologies for RNA modifications

To date, more than 170 chemical modifications have been characterized in RNA, providing a new layer of gene expression regulation termed the 'epitranscriptome'. RNA modification detection methods and tools advance the functional studies of the epitranscriptome. According to the detection throughput and principles, existing RNA modification detection technologies can be categorized into four classes, including quantification methods, locus-specific detection methods, next-generation sequencing-based detection technologies and nanopore direct RNA sequencing-based technologies. In this review, we summarize the current knowledge about these RNA modification detection technologies and discuss the challenges for the existing detection tools, providing information for a comprehensive understanding of the epitranscriptome.

tRNA-m1A modification promotes T cell expansion via efficient MYC protein synthesis

Naive T cells undergo radical changes during the transition from dormant to hyperactive states upon activation, which necessitates de novo protein production via transcription and translation. However, the mechanism whereby T cells globally promote translation remains largely unknown. Here, we show that on exit from quiescence, T cells upregulate transfer RNA (tRNA) m¹A58 'writer' proteins TRMT61A and TRMT6, which confer m¹A58 RNA modification on a specific subset of early expressed tRNAs. These m¹A-modified early tRNAs enhance translation efficiency, enabling rapid and necessary synthesis of MYC and of a specific group of key functional proteins. The MYC protein then guides the exit of naive T cells from a quiescent state into a proliferative state and promotes rapid T cell expansion after activation. Conditional deletion of the Trmt61a gene in mouse CD4⁺ T cells causes MYC protein deficiency and cell cycle arrest, disrupts T cell expansion upon cognate antigen stimulation and alleviates colitis in a mouse adoptive transfer colitis model. Our study elucidates for the first time, to our knowledge, the in vivo physiological roles of tRNA-m¹A58 modification in T cell-mediated pathogenesis and reveals a new mechanism of tRNA-m¹A58-controlled T cell homeostasis and signal-dependent translational control of specific key proteins.

Analysis of queuosine and 2-thio tRNA modifications by high throughput sequencing

Queuosine (Q) is a conserved tRNA modification at the wobble anticodon position of tRNAs that read the codons of amino acids Tyr, His, Asn, and Asp. Q-modification in tRNA plays important roles in the regulation of translation efficiency and fidelity. Queuosine tRNA modification is synthesized de novo in bacteria, whereas in mammals the substrate for Q-modification in tRNA is queuine, the catabolic product of the Q-base of gut bacteria. This gut microbiome dependent tRNA modification may play pivotal roles in translational regulation in different cellular contexts, but extensive studies of Q-modification biology are hindered by the lack of high throughput sequencing methods for its detection and quantitation. Here, we describe a periodate-treatment method that enables single base resolution profiling of Q-modification in tRNAs by Nextgen sequencing from biological RNA samples. Periodate oxidizes the Q-base, which results in specific deletion signatures in the RNA-seq data. Unexpectedly, we found that periodate-treatment also enables the detection of several 2-thio-modifications including τm5s2U, mcm5s2U, cmnm5s2U, and s2C by sequencing in human and E. coli tRNA. We term this method periodate-dependent analysis of queuosine and sulfur modification sequencing (PAQS-seq). We assess Q- and 2-thio-modifications at the tRNA isodecoder level, and 2-thio modification changes in stress response. PAQS-seq should be widely applicable in the biological studies of Q- and 2-thio-modifications in mammalian and microbial tRNAs.

The role of RNA modification in hepatocellular carcinoma

Hepatocellular carcinoma (HCC) is a highly mortal type of primary liver cancer. Abnormal epigenetic modifications are present in HCC, and RNA modification is dynamic and reversible and is a key post-transcriptional regulator. With the in-depth study of post-transcriptional modifications, RNA modifications are aberrantly expressed in human cancers. Moreover, the regulators of RNA modifications can be used as potential targets for cancer therapy. In RNA modifications, N6-methyladenosine (m6A), N7-methylguanosine (m7G), and 5-methylcytosine (m5C) and their regulators have important regulatory roles in HCC progression and represent potential novel biomarkers for the confirmation of diagnosis and treatment of HCC. This review focuses on RNA modifications in HCC and the roles and mechanisms of m6A, m7G, m5C, N1-methyladenosine (m1A), N3-methylcytosine (m3C), and pseudouridine (ψ) on its development and maintenance. The potential therapeutic strategies of RNA modifications are elaborated for HCC.

The tRNA regulome in neurodevelopmental and neuropsychiatric disease

Transfer (t)RNAs are 70-90 nucleotide small RNAs highly regulated by 43 different types of epitranscriptomic modifications and requiring aminoacylation ('charging') for mRNA decoding and protein synthesis. Smaller cleavage products of mature tRNAs, or tRNA fragments, have been linked to a broad variety of noncanonical functions, including translational inhibition and modulation of the immune response. Traditionally, knowledge about tRNA regulation in brain is derived from phenotypic exploration of monogenic neurodevelopmental and neurodegenerative diseases associated with rare mutations in tRNA modification genes. More recent studies point to the previously unrecognized potential of the tRNA regulome to affect memory, synaptic plasticity, and affective states. For example, in mature cortical neurons, cytosine methylation sensitivity of the glycine tRNA family (tRNA^Gly) is coupled to glycine biosynthesis and codon-specific alterations in ribosomal translation together with robust changes in cognition and depression-related behaviors. In this Review, we will discuss the emerging knowledge of the neuronal tRNA landscape, with a focus on epitranscriptomic tRNA modifications and downstream molecular pathways affected by alterations in tRNA expression, charging levels, and cleavage while mechanistically linking these pathways to neuropsychiatric disease and provide insight into future areas of study for this field.

Characterization of novel small non-coding RNAs and their modifications in bladder cancer using an updated small RNA-seq workflow

Background: Bladder cancer (BLCA) is one of the most common cancer types worldwide. The disease is responsible for about 200,000 deaths annually, thus improved diagnostics and therapy is needed. A large body of evidence reveal that small RNAs of less than 40 nucleotides may act as tumor suppressors, oncogenes, and disease biomarkers, with a major focus on microRNAs. However, the role of other families of small RNAs is not yet deciphered. Recent results suggest that small RNAs and their modification status, play a role in BLCA development and are promising biomarkers due to their high abundance in the exomes and body fluids (including urine). Moreover, free modified nucleosides have been detected at elevated levels from the urine of BLCA patients. A genome-wide view of small RNAs, and their modifications, will help pinpoint the molecules that could be used as biomarker or has important biology in BLCA development. Methods: BLCA tumor tissue specimens were obtained from 12 patients undergoing transurethral resection of non-muscle invasive papillary urothelial carcinomas. Genome-wide profiling of small RNAs less than 40 bases long was performed by a modified protocol with TGIRT (thermostable group II reverse transcriptase) to identify novel small RNAs and their modification status. Results: Comprehensive analysis identified not only microRNAs. Intriguingly, 57 ± 15% (mean ± S.D.) of sequencing reads mapped to non-microRNA-small RNAs including tRNA-derived fragments (tRFs), ribosomal RNA-derived fragments (rRFs) and YRNA-derived fragments (YRFs). Misincorporation (mismatch) sites identified potential base modification positions on the small RNAs, especially on tRFs, corresponding to m1A (N1-methyladenosine), m1G (N1-methylguanosine) and m2 2G (N2, N2-dimethylguanosine). We also detected mismatch sites on rRFs corresponding to known modifications on 28 and 18S rRNA. Conclusion: We found abundant non-microRNA-small RNAs in BLCA tumor samples. Small RNAs, especially tRFs and rRFs, contain modifications that can be captured as mismatch by TGIRT sequencing. Both the modifications and the non-microRNA-small RNAs should be explored as a biomarker for BLCA detection or follow-up.

Development of Mild Chemical Catalysis Conditions for m1A-to-m6A Rearrangement on RNA

The conversion of N¹-methyladenosine (m¹A) to N⁶-methyladenosine (m⁶A) on RNA is an important step for both allowing efficient reverse transcription read-though for sequencing analysis and mapping modifications in the transcriptome. Enzymatic transformation is often used, but the efficiency of the removal can depend on local sequence context. Chemical conversion through the application of the Dimroth rearrangement, in which m¹A rearranges into m⁶A under heat and alkaline conditions, is an alternative, but the required alkaline conditions result in significant RNA degradation by hydrolysis of the phosphodiester backbone. Here, we report novel, mild pH conditions that catalyze m¹A-to-m⁶A arrangement using 4-nitrothiophenol as a catalyst. We demonstrate the efficient rearrangement in mononucleosides, synthetic RNA oligonucleotides, and RNAs isolated from human cell lines, thereby validating a new approach for converting m¹A-to-m⁶A in RNA samples for sequencing analyses.

A multiplex platform for small RNA sequencing elucidates multifaceted tRNA stress response and translational regulation

Small RNAs include tRNA, snRNA, micro-RNA, tRNA fragments and others that constitute > 90% of RNA copy numbers in a human cell and perform many essential functions. Popular small RNA-seq strategies limit the insights into coordinated small RNA response to cellular stress. Small RNA-seq also lacks multiplexing capabilities. Here, we report a multiplex small RNA-seq library preparation method (MSR-seq) to investigate cellular small RNA and mRNA response to heat shock, hydrogen peroxide, and arsenite stress. Comparing stress-induced changes of total cellular RNA and polysome-associated RNA, we identify a coordinated tRNA response that involves polysome-specific tRNA abundance and synergistic N3-methylcytosine (m³C) tRNA modification. Combining tRNA and mRNA response to stress we reveal a mechanism of stress-induced down-regulation in translational elongation. We also find that native tRNA molecules lacking several modifications are biased reservoirs for the biogenesis of tRNA fragments. Our results demonstrate the importance of simultaneous investigation of small RNAs and their modifications in response to varying biological conditions.

tReasure: R-based GUI package analyzing tRNA expression profiles from small RNA sequencing data

BACKGROUND

Recent deep sequencing technologies have proven to be valuable resources to gain insights into the expression profiles of diverse tRNAs. However, despite these technologies, the association of tRNAs with diverse diseases has not been explored in depth because analytical tools are lacking.

RESULTS

We developed a user-friendly tool, tRNA Expression Analysis Software Utilizing R for Easy use (tReasure), to analyze differentially expressed tRNAs (DEtRNAs) from deep sequencing data of small RNAs using R packages. tReasure can quantify individual mature tRNAs, isodecoders, and isoacceptors. By adopting stringent mapping strategies, tReasure supports the precise measurement of mature tRNA read counts. The whole analysis workflow for determining DEtRNAs (uploading FASTQ files, removing adapter sequences and poor-quality reads, mapping and quantifying tRNAs, filtering out low count tRNAs, determining DEtRNAs, and visualizing statistical analysis) can be performed with the tReasure package.

CONCLUSIONS

tReasure is an open-source software available for download at https://treasure.pmrc.re.kr and will be indispensable for users who have little experience with command-line software to explore the biological implication of tRNA expression.

RNA Post-Transcriptional Modifications in Two Large Subunit Intermediates Populated in E. coli Cells Expressing Helicase Inactive R331A DbpA

23S ribosomal RNA (rRNA) of Escherichia coli 50S large ribosome subunit contains 26 post-transcriptionally modified nucleosides. Here, we determine the extent of modifications in the 35S and 45S large subunit intermediates, accumulating in cells expressing the helicase inactive DbpA protein, R331A, and the native 50S large subunit. The modifications we characterized are 3-methylpseudouridine, 2-methyladenine, 5-hydroxycytidine, and nine pseudouridines. These modifications were detected using 1-cyclohexyl-3-(2-morpholinoethyl)carbodiimide metho-p-toluenesulfonate (CMCT) treatment followed by alkaline treatment. In addition, KMnO₄ treatment of 23S rRNA was employed to detect 5-hydroxycytidine modification. CMCT and KMnO₄ treatments produce chemical changes in modified nucleotides that cause reverse transcriptase misincorporations and deletions, which were detected employing next-generation sequencing. Our results show that the 2-methyladenine modification and seven uridines to pseudouridine isomerizations are present in both the 35S and 45S to similar extents as in the 50S. Hence, the enzymes that perform these modifications, namely, RluA, RluB, RluC, RluE, RluF, and RlmN, have already acted in the intermediates. Two uridines to pseudouridine isomerizations, the 3-methylpseudouridine and 5-hydroxycytidine modifications, are significantly less present in the 35S and 45S, as compared to the 50S. Therefore, the enzymes that incorporate these modifications, RluD, RlmH, and RlhA, are in the process of modifying the 35S and 45S or will incorporate these modifications during the later stages of ribosome assembly. Our study employs a novel high throughput and single nucleotide resolution technique for the detection of 2-methyladenine and two novel high throughput and single nucleotide resolution techniques for the detection of 5-hydroxycytidine.

TRMT6/61A-dependent base methylation of tRNA-derived fragments regulates gene-silencing activity and the unfolded protein response in bladder cancer

RNA modifications are important regulatory elements of RNA functions. However, most genome-wide mapping of RNA modifications has focused on messenger RNAs and transfer RNAs, but such datasets have been lacking for small RNAs. Here we mapped N1-methyladenosine (m1A) in the cellular small RNA space. Benchmarked with synthetic m1A RNAs, our workflow identified specific groups of m1A-containing small RNAs, which are otherwise disproportionally under-represented. In particular, 22-nucleotides long 3' tRNA-fragments are highly enriched for TRMT6/61A-dependent m1A located within the seed region. TRMT6/61A-dependent m1A negatively affects gene silencing by tRF-3s. In urothelial carcinoma of the bladder, where TRMT6/61A is over-expressed, higher m1A modification on tRFs is detected, correlated with a dysregulation of tRF targetome. Lastly, TRMT6/61A regulates tRF-3 targets involved in unfolded protein response. Together, our results reveal a mechanism of regulating gene expression via base modification of small RNA.