Direct Nanopore Sequencing of Individual Full Length tRNA Strands

Combining Nanopore direct RNA sequencing with genetics and mass spectrometry for analysis of T-loop base modifications across 42 yeast tRNA isoacceptors

Transfer RNAs (tRNAs) contain dozens of chemical modifications. These modifications are critical for maintaining tRNA tertiary structure and optimizing protein synthesis. Here we advance the use of Nanopore direct RNA-sequencing (DRS) to investigate the synergy between modifications that are known to stabilize tRNA structure. We sequenced the 42 cytosolic tRNA isoacceptors from wild-type yeast and five tRNA-modifying enzyme knockout mutants. These data permitted comprehensive analysis of three neighboring and conserved modifications in T-loops: 5-methyluridine (m5U54), pseudouridine (Ψ55), and 1-methyladenosine (m1A58). Our results were validated using direct measurements of chemical modifications by mass spectrometry. We observed concerted T-loop modification circuits-the potent influence of Ψ55 for subsequent m1A58 modification on more tRNA isoacceptors than previously observed. Growing cells under nutrient depleted conditions also revealed a novel condition-specific increase in m1A58 modification on some tRNAs. A global and isoacceptor-specific classification strategy was developed to predict the status of T-loop modifications from a user-input tRNA DRS dataset, applicable to other conditions and tRNAs in other organisms. These advancements demonstrate how orthogonal technologies combined with genetics enable precise detection of modification landscapes of individual, full-length tRNAs, at transcriptome-scale.

TRMT1L-catalyzed m2 2G27 on tyrosine tRNA is required for efficient mRNA translation and cell survival under oxidative stress

tRNA modifications are critical for several aspects of their functions, including decoding, folding, and stability. Using a multifaceted approach encompassing eCLIP-seq and Nanopore tRNA-seq, we show that the human tRNA methyltransferase TRMT1L interacts with component of the Rix1 ribosome biogenesis complex and binds to the 28S rRNA, as well as to a subset of tRNAs. Mechanistically, we demonstrate that TRMT1L is responsible for catalyzing m2 2G solely at position 27 of tRNA-Tyr-GUA. Surprisingly, TRMT1L depletion also impaired the deposition of acp3U and dihydrouridine on tRNA-Tyr-GUA, Cys-GCA, and Ala-CGC. TRMT1L knockout cells have a marked decrease in tRNA-Tyr-GUA levels, coinciding with a reduction in global translation rates and hypersensitivity to oxidative stress. Our results establish TRMT1L as the elusive methyltransferase catalyzing the m2 2G27 modification on tRNA Tyr, resolving a long-standing gap of knowledge and highlighting its potential role in a tRNA modification circuit crucial for translation regulation and stress response.

Mass Spectrometry-Based Direct Sequencing of tRNAs De Novo and Quantitative Mapping of Multiple RNA Modifications

Despite the extensive use of next-generation sequencing (NGS) of RNA, simultaneous direct sequencing and quantitative mapping of multiple RNA nucleotide modifications remains challenging. Mass spectrometry (MS)-based sequencing can directly sequence all RNA modifications without being limited to specific ones, but it requires a perfect MS ladder that few tRNAs can provide. Here, we describe an MS ladder complementation sequencing approach (MLC-Seq) that circumvents the perfect ladder requirement, allowing de novo MS sequencing of full-length heterogeneous cellular tRNAs with multiple nucleotide modifications at single-nucleotide precision. Unlike NGS-based methods, which lose RNA modification information, MLC-Seq preserves RNA sequence diversity and modification information, revealing new detailed stoichiometric tRNA modification profiles and their changes upon treatment with the dealkylating enzyme AlkB. It can also be combined with reference sequences to provide quantitative analysis of diverse tRNAs and modifications in total tRNA samples. MLC-Seq enables systematic, quantitative, and site-specific mapping of RNA modifications, revealing the truly complete informational content of tRNA.

The Precise Basecalling of Short-Read Nanopore Sequencing

The nanopore sequencing of short sequences, whose lengths are typically less than 0.3kb therefore comparable with Illumina sequencing techniques, has recently gained wide attention. Here, we design a scheme for training nanopore basecallers that are specialized for short biomolecules. With bioengineered RNA (BioRNA) molecules as examples, we demonstrate the superior accuracy of basecallers trained by our scheme.

Mitochondrial transcriptome of Candida albicans in flagranti - direct RNA sequencing reveals a new layer of information

BACKGROUND

Organellar transcriptomes are relatively under-studied systems, with data related to full-length transcripts and posttranscriptional modifications remaining sparse. Direct RNA sequencing presents the possibility of accessing a previously unavailable layer of information pertaining to transcriptomic data, as well as circumventing the biases introduced by second-generation RNA-seq platforms. Direct long-read ONT sequencing allows for the isoform analysis of full-length transcripts and the detection of posttranscriptional modifications. However, there are still relatively few projects employing this method specifically for studying organellar transcriptomes.

RESULTS

Candida albicans is a promising model for investigating nucleo-mitochondrial interactions. This work comprises ONT sequencing of the Candida albicans mitochondrial transcriptome along with the development of a dedicated data analysis pipeline. This approach allowed for the detection of complete transcript isoforms and posttranslational RNA modifications, as well as an analysis of C. albicans deletion mutants in genes coding for the 5' and 3' mitochondrial RNA exonucleases CaPET127 and CaDSS1. It also enabled for corrections to previous studies in terms of 3' and 5' transcript ends. A number of intermediate splicing isoforms was also discovered, along with mature and unspliced transcripts and changes in their abundances resulting from disruption of both 5' and 3' exonucleolytic processing. Multiple putative posttranscriptional modification sites have also been detected.

CONCLUSIONS

This preliminary work demonstrates the suitability of direct RNA sequencing for studying yeast mitochondrial transcriptomes in general and provides new insights into the workings of the C. albicans mitochondrial transcriptome in particular. It also provides a general roadmap for analyzing mitochondrial transcriptomic data from other organisms.

RNA modifying enzymes shape tRNA biogenesis and function

Transfer RNAs (tRNAs) are the most highly modified cellular RNAs, both with respect to the proportion of nucleotides that are modified within the tRNA sequence and with respect to the extraordinary diversity in tRNA modification chemistry. However, the functions of many different tRNA modifications are only beginning to emerge. tRNAs have two general clusters of modifications. The first cluster is within the anticodon stem-loop including several modifications essential for protein translation. The second cluster of modifications is within the tRNA elbow, and roles for these modifications are less clear. In general, tRNA elbow modifications are typically not essential for cell growth, but nonetheless several tRNA elbow modifications have been highly conserved throughout all domains of life. In addition to forming modifications, many tRNA modifying enzymes have been demonstrated or hypothesized to also play an important role in folding tRNA acting as tRNA chaperones. In this review, we summarize the known functions of tRNA modifying enzymes throughout the lifecycle of a tRNA molecule, from transcription to degradation. Thereby, we describe how tRNA modification and folding by tRNA modifying enzymes enhance tRNA maturation, tRNA aminoacylation, and tRNA function during protein synthesis, ultimately impacting cellular phenotypes and disease.

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.

Comparative analysis of 43 distinct RNA modifications by nanopore tRNA sequencing

Transfer RNAs are the fundamental adapter molecules of protein synthesis and the most abundant and heterogeneous class of noncoding RNA molecules in cells. The study of tRNA repertoires remains challenging, complicated by the presence of dozens of post transcriptional modifications. Nanopore sequencing is an emerging technology with promise for both tRNA sequencing and the detection of RNA modifications; however, such studies have been limited by the throughput and accuracy of direct RNA sequencing methods. Moreover, detection of the complete set of tRNA modifications by nanopore sequencing remains challenging. Here we show that recent updates to nanopore direct RNA sequencing chemistry (RNA004) combined with our own optimizations to tRNA sequencing protocols and analysis workflows enable high throughput coverage of tRNA molecules and characterization of nanopore signals produced by 43 distinct RNA modifications. We share best practices and protocols for nanopore sequencing of tRNA and further report successful detection of low abundance mitochondrial and viral tRNAs, providing proof of concept for use of nanopore sequencing to study tRNA populations in the context of infection and organelle biology. This work provides a roadmap to guide future efforts towards de novo detection of RNA modifications across multiple organisms using nanopore sequencing.

Beet Soil-Borne Virus Is a Helper Virus for the Novel Beta vulgaris Satellite Virus 1A

Sugar beet (Beta vulgaris) is grown in temperate regions around the world as a source of sucrose used for natural sweetening. Sugar beet is susceptible to a number of viral diseases, but identification of the causal agent(s) under field conditions is often difficult due to mixtures of viruses that may be responsible for disease symptoms. In this study, the application of RNAseq to RNA extracted from diseased sugar beet roots obtained from the field and from greenhouse-reared plants grown in soil infested with the virus disease rhizomania (causal agent beet necrotic yellow vein virus; BNYVV) yielded genome-length sequences from BNYVV, as well as beet soil-borne virus (BSBV). The nucleotide identities of the derived consensus sequence of BSBV RNAs ranged from 99.4 to 96.7% (RNA1), 99.3 to 95.3% (RNA2), and 98.3 to 95.9% (RNA3) compared with published BSBV sequences. Based on the BSBV genome consensus sequence, clones of the genomic RNAs 1, 2, and 3 were obtained to produce RNA copies of the genome through in vitro transcription. Capped RNA produced from the clones was infectious when inoculated into leaves of Chenopodium quinoa and B. vulgaris, and extracts from transcript-infected C. quinoa leaves could infect sugar beet seedling roots through a vortex inoculation method. Subsequent exposure of these infected sugar beet seedling roots to aviruliferous Polymyxa betae, the protist vector of both BNYVV and BSBV, confirmed that BSBV derived from the infectious clones could be transmitted by the vector. Co-inoculation of BSBV synthetic transcripts with transcripts of a cloned putative satellite virus designated Beta vulgaris satellite virus 1A (BvSat1A) resulted in the production of lesions on leaves of C. quinoa similar to those produced by inoculation with BSBV alone. Nevertheless, accumulation of genomic RNA and the encoded protein of the satellite virus in co-inoculated leaves was readily detected on Northern and Western blots, respectively, whereas no accumulation of satellite virus products occurred when satellite virus RNA was inoculated alone. The predicted sequence of the detected protein encoded by BvSat1A bears hallmarks of coat proteins of other satellite viruses, and virions of a size consistent with a satellite virus were observed in samples testing positive for the virus. The results demonstrate that BSBV is a helper virus for the novel satellite virus BvSat1A.

Long-read direct RNA sequencing of the mitochondrial transcriptome of Saccharomyces cerevisiae reveals condition-dependent intron abundance

Mitochondria fulfil many essential roles and have their own genome, which is expressed as polycistronic transcripts that undergo co- or posttranscriptional processing and splicing. Due to the inherent complexity and limited technical accessibility of the mitochondrial transcriptome, fundamental questions regarding mitochondrial gene expression and splicing remain unresolved, even in the model eukaryote Saccharomyces cerevisiae. Long-read sequencing could address these fundamental questions. Therefore, a method for the enrichment of mitochondrial RNA and sequencing using Nanopore technology was developed, enabling the resolution of splicing of polycistronic genes and the quantification of spliced RNA. This method successfully captured the full mitochondrial transcriptome and resolved RNA splicing patterns with single-base resolution and was applied to explore the transcriptome of S. cerevisiae grown with glucose or ethanol as the sole carbon source, revealing the impact of growth conditions on mitochondrial RNA expression and splicing. This study uncovered a remarkable difference in the turnover of Group II introns between yeast grown in either mostly fermentative or fully respiratory conditions. Whether this accumulation of introns in glucose medium has an impact on mitochondrial functions remains to be explored. Combined with the high tractability of the model yeast S. cerevisiae, the developed method enables to monitor mitochondrial transcriptome responses in a broad range of relevant contexts, including oxidative stress, apoptosis and mitochondrial diseases.

Beyond the Anticodon: tRNA Core Modifications and Their Impact on Structure, Translation and Stress Adaptation

Transfer RNAs (tRNAs) are heavily decorated with post-transcriptional chemical modifications. Approximately 100 different modifications have been identified in tRNAs, and each tRNA typically contains 5-15 modifications that are incorporated at specific sites along the tRNA sequence. These modifications may be classified into two groups according to their position in the three-dimensional tRNA structure, i.e., modifications in the tRNA core and modifications in the anticodon-loop (ACL) region. Since many modified nucleotides in the tRNA core are involved in the formation of tertiary interactions implicated in tRNA folding, these modifications are key to tRNA stability and resistance to RNA decay pathways. In comparison to the extensively studied ACL modifications, tRNA core modifications have generally received less attention, although they have been shown to play important roles beyond tRNA stability. Here, we review and place in perspective selected data on tRNA core modifications. We present their impact on tRNA structure and stability and report how these changes manifest themselves at the functional level in translation, fitness and stress adaptation.

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.

Nanopore DNA sequencing technologies and their applications towards single-molecule proteomics

Sequencing of nucleic acids with nanopores has emerged as a powerful tool offering rapid readout, high accuracy, low cost and portability. This label-free method for sequencing at the single-molecule level is an achievement on its own. However, nanopores also show promise for the technologically even more challenging sequencing of polypeptides, something that could considerably benefit biological discovery, clinical diagnostics and homeland security, as current techniques lack portability and speed. Here we survey the biochemical innovations underpinning commercial and academic nanopore DNA/RNA sequencing techniques, and explore how these advances can fuel developments in future protein sequencing with nanopores.

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.

Chemically Modified Platforms for Better RNA Therapeutics

RNA-based therapies have catalyzed a revolutionary transformation in the biomedical landscape, offering unprecedented potential in disease prevention and treatment. However, despite their remarkable achievements, these therapies encounter substantial challenges including low stability, susceptibility to degradation by nucleases, and a prominent negative charge, thereby hindering further development. Chemically modified platforms have emerged as a strategic innovation, focusing on precise alterations either on the RNA moieties or their associated delivery vectors. This comprehensive review delves into these platforms, underscoring their significance in augmenting the performance and translational prospects of RNA-based therapeutics. It encompasses an in-depth analysis of various chemically modified delivery platforms that have been instrumental in propelling RNA therapeutics toward clinical utility. Moreover, the review scrutinizes the rationale behind diverse chemical modification techniques aiming at optimizing the therapeutic efficacy of RNA molecules, thereby facilitating robust disease management. Recent empirical studies corroborating the efficacy enhancement of RNA therapeutics through chemical modifications are highlighted. Conclusively, we offer profound insights into the transformative impact of chemical modifications on RNA drugs and delineates prospective trajectories for their future development and clinical integration.

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.

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.

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.

Detection of queuosine and queuosine precursors in tRNAs by direct RNA sequencing

Queuosine (Q) is a complex tRNA modification found in bacteria and eukaryotes at position 34 of four tRNAs with a GUN anticodon, and it regulates the translational efficiency and fidelity of the respective codons that differ at the Wobble position. In bacteria, the biosynthesis of Q involves two precursors, preQ0 and preQ1, whereas eukaryotes directly obtain Q from bacterial sources. The study of queuosine has been challenging due to the limited availability of high-throughput methods for its detection and analysis. Here, we have employed direct RNA sequencing using nanopore technology to detect the modification of tRNAs with Q and Q precursors. These modifications were detected with high accuracy on synthetic tRNAs as well as on tRNAs extracted from Schizosaccharomyces pombe and Escherichia coli by comparing unmodified to modified tRNAs using the tool JACUSA2. Furthermore, we present an improved protocol for the alignment of raw sequence reads that gives high specificity and recall for tRNAs ex cellulo that, by nature, carry multiple modifications. Altogether, our results show that 7-deazaguanine-derivatives such as queuosine are readily detectable using direct RNA sequencing. This advancement opens up new possibilities for investigating these modifications in native tRNAs, furthering our understanding of their biological function.

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.

Advantages and challenges associated with bisulfite-assisted nanopore direct RNA sequencing for modifications

Nanopore direct RNA sequencing is a technology that allows sequencing for epitranscriptomic modifications with the possibility of a quantitative assessment. In the present work, pseudouridine (Ψ) was sequenced with the nanopore before and after the pH 7 bisulfite reaction that yields stable ribose adducts at C1' of Ψ. The adducted sites produced greater base call errors in the form of deletion signatures compared to Ψ. Sequencing studies on E. coli rRNA and tmRNA before and after the pH 7 bisulfite reaction demonstrated that using chemically-assisted nanopore sequencing has distinct advantages for minimization of false positives and false negatives in the data. The rRNA from E. coli has 19 known U/C sequence variations that give similar base call signatures as Ψ, and therefore, are false positives when inspecting base call data; however, these sites are refractory to reacting with bisulfite as is easily observed in nanopore data. The E. coli tmRNA has a low occupancy Ψ in a pyrimidine-rich sequence context that is called a U representing a false negative; partial occupancy by Ψ is revealed after the bisulfite reaction. In a final study, 5-methylcytidine (m5C) in RNA can readily be observed after the pH 5 bisulfite reaction in which the parent C deaminates to U and the modified site does not react. This locates m5C when using bisulfite-assisted nanopore direct RNA sequencing, which is otherwise challenging to observe. The advantages and challenges of the overall approach are discussed.

Different modification pathways for m1A58 incorporation in yeast elongator and initiator tRNAs

As essential components of the protein synthesis machinery, tRNAs undergo a tightly controlled biogenesis process, which include the incorporation of numerous posttranscriptional modifications. Defects in these tRNA maturation steps may lead to the degradation of hypomodified tRNAs by the rapid tRNA decay (RTD) and nuclear surveillance pathways. We previously identified m1A58 as a late modification introduced after modifications Ψ55 and T54 in yeast elongator tRNAPhe. However, previous reports suggested that m1A58 is introduced early during the tRNA modification process, in particular on primary transcripts of initiator tRNAiMet, which prevents its degradation by RNA decay pathways. Here, aiming to reconcile this apparent inconsistency on the temporality of m1A58 incorporation, we examined its introduction into yeast elongator and initiator tRNAs. We used specifically modified tRNAs to report on the molecular aspects controlling the Ψ55 → T54 → m1A58 modification circuit in elongator tRNAs. We also show that m1A58 is efficiently introduced on unmodified tRNAiMet, and does not depend on prior modifications. Finally, we show that m1A58 has major effects on the structural properties of initiator tRNAiMet, so that the tRNA elbow structure is only properly assembled when this modification is present. This observation provides a structural explanation for the degradation of hypomodified tRNAiMet lacking m1A58 by the nuclear surveillance and RTD pathways.

Direct Nanopore Sequencing for the 17 RNA Modification Types in 36 Locations in the E. coli Ribosome Enables Monitoring of Stress-Dependent Changes

The bacterium Escherichia coli possesses 16S and 23S rRNA strands that have 36 chemical modification sites with 17 different structures. Nanopore direct RNA sequencing using a protein nanopore sensor and helicase brake, which is also a sensor, was applied to the rRNAs. Nanopore current levels, base calling profile, and helicase dwell times for the modifications relative to unmodified synthetic rRNA controls found signatures for nearly all modifications. Signatures for clustered modifications were determined by selective sequencing of writer knockout E. coli and sequencing of synthetic RNAs utilizing some custom-synthesized nucleotide triphosphates for their preparation. The knowledge of each modification's signature, apart from 5-methylcytidine, was used to determine how metabolic and cold-shock stress impact rRNA modifications. Metabolic stress resulted in either no change or a decrease, and one site increased in modification occupancy, while cold-shock stress led to either no change or a decrease. The double modification m4Cm1402 resides in 16S rRNA, and it decreased with both stressors. Using the helicase dwell time, it was determined that the N4 methyl group is lost during both stressors, and the 2'-OMe group remained. In the ribosome, this modification stabilizes binding to the mRNA codon at the P-site resulting in increased translational fidelity that is lost during stress. The E. coli genome has seven rRNA operons (rrn), and the earlier studies aligned the nanopore reads to a single operon (rrnA). Here, the reads were aligned to all seven operons to identify operon-specific changes in the 11 pseudouridines. This study demonstrates that direct sequencing for >16 different RNA modifications in a strand is achievable.

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.

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.

Direct sequencing of total Saccharomyces cerevisiae tRNAs by LC-MS/MS

Among RNAs, transfer RNAs (tRNAs) contain the widest variety of abundant posttranscriptional chemical modifications. These modifications are crucial for tRNAs to participate in protein synthesis, promoting proper tRNA structure and aminoacylation, facilitating anticodon:codon recognition, and ensuring the reading frame maintenance of the ribosome. While tRNA modifications were long thought to be stoichiometric, it is becoming increasingly apparent that these modifications can change dynamically in response to the cellular environment. The ability to broadly characterize the fluctuating tRNA modification landscape will be essential for establishing the molecular level contributions of individual sites of tRNA modification. The locations of modifications within individual tRNA sequences can be mapped using liquid chromatography coupled to tandem mass spectrometry (LC-MS/MS). In this approach, a single tRNA species is purified, treated with ribonucleases, and the resulting single-stranded RNA products are subject to LC-MS/MS analysis. The application of LC-MS/MS to study tRNAs is limited by the necessity of analyzing one tRNA at a time, because the digestion of total tRNA mixtures by commercially available ribonucleases produces many short digestion products unable to be uniquely mapped back to a single site within a tRNA. We overcame these limitations by taking advantage of the highly structured nature of tRNAs to prevent the full digestion by single-stranded RNA-specific ribonucleases. Folding total tRNA prior to digestion allowed us to sequence Saccharomyces cerevisiae tRNAs with up to 97% sequence coverage for individual tRNA species by LC-MS/MS. This method presents a robust avenue for directly detecting the distribution of modifications in total tRNAs.

Applications of long-read sequencing to Mendelian genetics

Advances in clinical genetic testing, including the introduction of exome sequencing, have uncovered the molecular etiology for many rare and previously unsolved genetic disorders, yet more than half of individuals with a suspected genetic disorder remain unsolved after complete clinical evaluation. A precise genetic diagnosis may guide clinical treatment plans, allow families to make informed care decisions, and permit individuals to participate in N-of-1 trials; thus, there is high interest in developing new tools and techniques to increase the solve rate. Long-read sequencing (LRS) is a promising technology for both increasing the solve rate and decreasing the amount of time required to make a precise genetic diagnosis. Here, we summarize current LRS technologies, give examples of how they have been used to evaluate complex genetic variation and identify missing variants, and discuss future clinical applications of LRS. As costs continue to decrease, LRS will find additional utility in the clinical space fundamentally changing how pathological variants are discovered and eventually acting as a single-data source that can be interrogated multiple times for clinical service.

Nanopore sequencing of internal 2'-PO4 modifications installed by RNA repair

Ligation by plant and fungal RNA ligases yields an internal 2'-phosphate group on each RNA ligation product. In budding yeast, this covalent mark occurs at the splice junction of two targets of ligation: intron-containing tRNAs and the messenger RNA HAC1 The repertoire of RNA molecules repaired by RNA ligation has not been explored due to a lack of unbiased approaches for identifying RNA ligation products. Here, we define several unique signals produced by 2'-phosphorylated RNAs during nanopore sequencing. A 2'-phosphate at the splice junction of HAC1 mRNA inhibits 5' → 3' degradation, enabling detection of decay intermediates in yeast RNA repair mutants by nanopore sequencing. During direct RNA sequencing, intact 2'-phosphorylated RNAs on HAC1 and tRNAs produce diagnostic changes in nanopore current properties and base calling features, including stalls produced as the modified RNA translocates through the nanopore motor protein. These approaches enable directed studies to identify novel RNA repair events in other contexts.

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.

Nanopore sequencing for N1-methylpseudouridine in RNA reveals sequence-dependent discrimination of the modified nucleotide triphosphate during transcription

Direct RNA sequencing with a commercial nanopore platform was used to sequence RNA containing uridine (U), pseudouridine (Ψ) or N1-methylpseudouridine (m1Ψ) in >100 different 5-nucleotide contexts. The base calling data for Ψ or m1Ψ were similar but different from U allowing their detection. Understanding the nanopore signatures for Ψ and m1Ψ enabled a running start T7 RNA polymerase assay to study the selection of UTP versus ΨTP or m1ΨTP competing mixtures in all possible adjacent sequence contexts. A significant sequence context dependency was observed for T7 RNA polymerase with insertion yields for ΨTP versus UTP spanning a range of 20-65%, and m1ΨTP versus UTP producing variable yields that differ by 15-70%. Experiments with SP6 RNA polymerase, as well as chemically-modified triphosphates and DNA templates provide insight to explain the observations. The SP6 polymerase introduced m1ΨTP when competed with UTP with a smaller window of yields (15-30%) across all sequence contexts studied. These results may aid in future efforts that employ RNA polymerases to make therapeutic mRNAs with sub-stoichiometric amounts of m1Ψ.

An Introduction to Nanopore Sequencing: Past, Present, and Future Considerations

There has been significant progress made in the field of nanopore biosensor development and sequencing applications, which address previous limitations that restricted widespread nanopore use. These innovations, paired with the large-scale commercialization of biological nanopore sequencing by Oxford Nanopore Technologies, are making the platforms a mainstay in contemporary research laboratories. Equipped with the ability to provide long- and short read sequencing information, with quick turn-around times and simple sample preparation, nanopore sequencers are rapidly improving our understanding of unsolved genetic, transcriptomic, and epigenetic problems. However, there remain some key obstacles that have yet to be improved. In this review, we provide a general introduction to nanopore sequencing principles, discussing biological and solid-state nanopore developments, obstacles to single-base detection, and library preparation considerations. We present examples of important clinical applications to give perspective on the potential future of nanopore sequencing in the field of molecular diagnostics.

Present and future challenges for the investigation of transgenerational epigenetic inheritance

Epigenetic pathways are essential in different biological processes and in phenotype-environment interactions in response to different stressors and they can induce phenotypic plasticity. They encompass several processes that are mitotically and, in some cases, meiotically heritable, so they can be transferred to subsequent generations via the germline. Transgenerational Epigenetic Inheritance (TEI) describes the phenomenon that phenotypic traits, such as changes in fertility, metabolic function, or behavior, induced by environmental factors (e.g., parental care, pathogens, pollutants, climate change), can be transferred to offspring generations via epigenetic mechanisms. Investigations on TEI contribute to deciphering the role of epigenetic mechanisms in adaptation, adversity, and evolution. However, molecular mechanisms underlying the transmission of epigenetic changes between generations, and the downstream chain of events leading to persistent phenotypic changes, remain unclear. Therefore, inter-, (transmission of information between parental and offspring generation via direct exposure) and transgenerational (transmission of information through several generations with disappearance of the triggering factor) consequences of epigenetic modifications remain major issues in the field of modern biology. In this article, we review and describe the major gaps and issues still encountered in the TEI field: the general challenges faced in epigenetic research; deciphering the key epigenetic mechanisms in inheritance processes; identifying the relevant drivers for TEI and implement a collaborative and multi-disciplinary approach to study TEI. Finally, we provide suggestions on how to overcome these challenges and ultimately be able to identify the specific contribution of epigenetics in transgenerational inheritance and use the correct tools for environmental science investigation and biomarkers identification.

Nano3P-seq: transcriptome-wide analysis of gene expression and tail dynamics using end-capture nanopore cDNA sequencing

RNA polyadenylation plays a central role in RNA maturation, fate, and stability. In response to developmental cues, polyA tail lengths can vary, affecting the translation efficiency and stability of mRNAs. Here we develop Nanopore 3' end-capture sequencing (Nano3P-seq), a method that relies on nanopore cDNA sequencing to simultaneously quantify RNA abundance, tail composition, and tail length dynamics at per-read resolution. By employing a template-switching-based sequencing protocol, Nano3P-seq can sequence RNA molecule from its 3' end, regardless of its polyadenylation status, without the need for PCR amplification or ligation of RNA adapters. We demonstrate that Nano3P-seq provides quantitative estimates of RNA abundance and tail lengths, and captures a wide diversity of RNA biotypes. We find that, in addition to mRNA and long non-coding RNA, polyA tails can be identified in 16S mitochondrial ribosomal RNA in both mouse and zebrafish models. Moreover, we show that mRNA tail lengths are dynamically regulated during vertebrate embryogenesis at an isoform-specific level, correlating with mRNA decay. Finally, we demonstrate the ability of Nano3P-seq in capturing non-A bases within polyA tails of various lengths, and reveal their distribution during vertebrate embryogenesis. Overall, Nano3P-seq is a simple and robust method for accurately estimating transcript levels, tail lengths, and tail composition heterogeneity in individual reads, with minimal library preparation biases, both in the coding and non-coding transcriptome.

Parent tRNA Modification Status Determines the Induction of Functional tRNA-Derived RNA by Respiratory Syncytial Virus Infection

tRNA-derived RNA fragments (tRFs) are a recently discovered family of small noncoding RNAs (sncRNAs). We previously reported that respiratory syncytial virus (RSV) infection induces functional tRFs, which are derived from a limited subset of parent tRNAs, in airway epithelial cells. Such induction is also observed in nasopharyngeal wash samples from RSV patients and correlates to RSV genome copies, suggesting a clinical significance of tRFs in RSV infection. This work also investigates whether the modification of parent tRNAs is changed by RSV to induce tRFs, using one of the most inducible tRFs as a model. We discovered that RSV infection changed the methylation modification of adenine at position 57 in tRNA glutamic acid, with a codon of CTC (tRNA-GluCTC), and the change is essential for its cleavage. AlkB homolog 1, a previously reported tRNA demethylase, appears to remove methyladenine from tRNA-GluCTC, prompting the subsequent production of tRFs from the 5'-end of tRNA-GluCTC, a regulator of RSV replication. This study demonstrates for the first time the importance of post-transcriptional modification of tRNAs in tRF biogenesis following RSV infection, providing critical insights for antiviral strategy development.

tModBase: deciphering the landscape of tRNA modifications and their dynamic changes from epitranscriptome data

tRNA molecules contain dense, abundant modifications that affect tRNA structure, stability, mRNA decoding and tsRNA formation. tRNA modifications and related enzymes are responsive to environmental cues and are associated with a range of physiological and pathological processes. However, there is a lack of resources that can be used to mine and analyse these dynamically changing tRNA modifications. In this study, we established tModBase (https://www.tmodbase.com/) for deciphering the landscape of tRNA modification profiles from epitranscriptome data. We analysed 103 datasets generated with second- and third-generation sequencing technologies and illustrated the misincorporation and termination signals of tRNA modification sites in ten species. We thus systematically demonstrate the modification profiles across different tissues/cell lines and summarize the characteristics of tRNA-associated human diseases. By integrating transcriptome data from 32 cancers, we developed novel tools for analysing the relationships between tRNA modifications and RNA modification enzymes, the expression of 1442 tRNA-derived small RNAs (tsRNAs), and 654 DNA variations. Our database will provide new insights into the features of tRNA modifications and the biological pathways in which they participate.

Modification mapping by nanopore sequencing

Next generation sequencing (NGS) has provided biologists with an unprecedented view into biological processes and their regulation over the past 2 decades, fueling a wave of development of high throughput methods based on short read DNA and RNA sequencing. For nucleic acid modifications, NGS has been coupled with immunoprecipitation, chemical treatment, enzymatic treatment, and/or the use of reverse transcriptase enzymes with fortuitous activities to enrich for and to identify covalent modifications of RNA and DNA. However, the majority of nucleic acid modifications lack commercial monoclonal antibodies, and mapping techniques that rely on chemical or enzymatic treatments to manipulate modification signatures add additional technical complexities to library preparation. Moreover, such approaches tend to be specific to a single class of RNA or DNA modification, and generate only indirect readouts of modification status. Third generation sequencing technologies such as the commercially available "long read" platforms from Pacific Biosciences and Oxford Nanopore Technologies are an attractive alternative for high throughput detection of nucleic acid modifications. While the former can indirectly sense modified nucleotides through changes in the kinetics of reverse transcription reactions, nanopore sequencing can in principle directly detect any nucleic acid modification that produces a signal distortion as the nucleic acid passes through a nanopore sensor embedded within a charged membrane. To date, more than a dozen endogenous DNA and RNA modifications have been interrogated by nanopore sequencing, as well as a number of synthetic nucleic acid modifications used in metabolic labeling, structure probing, and other emerging applications. This review is intended to introduce the reader to nanopore sequencing and key principles underlying its use in direct detection of nucleic acid modifications in unamplified DNA or RNA samples, and outline current approaches for detecting and quantifying nucleic acid modifications by nanopore sequencing. As this technology matures, we anticipate advances in both sequencing chemistry and analysis methods will lead to rapid improvements in the identification and quantification of these epigenetic marks.

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.

Organization and epigenomic control of RNA polymerase III-transcribed genes in plants

The genetic information linearly scripted in chromosomes is wrapped in a ribonucleoprotein complex called chromatin. The adaptation of its compaction level and spatiotemporal organization refines gene expression in response to developmental and environmental cues. RNA polymerase III (RNAPIII) is responsible for the biogenesis of elementary non-coding RNAs. Their genes are subjected to high duplication and mutational rates, and invade nuclear genomes. Their insertion into different epigenomic environments raises the question of how chromatin packing affects their individual transcription. In this review, we provide a unique perspective to this issue in plants. In addition, we discuss how the genomic organization of RNAPIII-transcribed loci, combined with epigenetic differences, might participate to plant trait variations.

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

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

Ribosomal RNA Pseudouridylation: Will Newly Available Methods Finally Define the Contribution of This Modification to Human Ribosome Plasticity?

In human rRNA, at least 104 specific uridine residues are modified to pseudouridine. Many of these pseudouridylation sites are located within functionally important ribosomal domains and can influence ribosomal functional features. Until recently, available methods failed to reliably quantify the level of modification at each specific rRNA site. Therefore, information obtained so far only partially explained the degree of regulation of pseudouridylation in different physiological and pathological conditions. In this focused review, we provide a summary of the methods that are now available for the study of rRNA pseudouridylation, discussing the perspectives that newly developed approaches are offering.

The Dihydrouridine landscape from tRNA to mRNA: a perspective on synthesis, structural impact and function

The universal dihydrouridine (D) epitranscriptomic mark results from a reduction of uridine by the Dus family of NADPH-dependent reductases and is typically found within the eponym D-loop of tRNAs. Despite its apparent simplicity, D is structurally unique, with the potential to deeply affect the RNA backbone and many, if not all, RNA-connected processes. The first landscape of its occupancy within the tRNAome was reported 20 years ago. Its potential biological significance was highlighted by observations ranging from a strong bias in its ecological distribution to the predictive nature of Dus enzymes overexpression for worse cancer patient outcomes. The exquisite specificity of the Dus enzymes revealed by a structure-function analyses and accumulating clues that the D distribution may expand beyond tRNAs recently led to the development of new high-resolution mapping methods, including Rho-seq that established the presence of D within mRNAs and led to the demonstration of its critical physiological relevance.

Bacterial tRNA landscape revisited

The updated Wobble Hypothesis reasonably explains why some 40 tRNA species are sufficient to decode the 61 amino acid codons of the Universal Genetic Code. However, we still have no clue why eubacteria lack tRNA isoacceptors with ANN anticodons, whereas eukaryotes universally lack eight GNN anticodons, only one of which is also absent in bacteria. Direct tRNA sequencing could resolve the patterns of nucleoside modification that had been driving the divergent evolution in prokaryotes and eukaryotes, but this task will require the development of AI-supported base-callers that can recognize modified nucleosides without any subsequent analytical verification. Our knowledge of the bacterial tRNA landscape is moreover broadened by the recent discovery of antisense tRNAs and tRNA-derived fragments that should be examined in their roles for gene expression, translation, bacterial physiology or metabolism.

Emerging roles for tRNAs in hematopoiesis and hematological malignancies

tRNAs are central players in decoding the genetic code linking codons in mRNAs with cognate amino acids during protein synthesis. Recent discoveries have placed tRNAs as key regulators of gene expression during hematopoiesis, especially in hematopoietic stem cell (HSC) maintenance and immune development. These functions have been shown to be influenced by dynamic changes in tRNA expression, post-transcriptional base modifications, tRNA-interacting proteins, and tRNA fragmentation; these events underlie the complexity of tRNA-mediated regulatory events in hematopoiesis. In this review, we discuss these recent findings and highlight how deregulation of tRNA biogenesis can contribute to hematological malignancies.

Concerted modification of nucleotides at functional centers of the ribosome revealed by single-molecule RNA modification profiling

Nucleotides in RNA and DNA are chemically modified by numerous enzymes that alter their function. Eukaryotic ribosomal RNA (rRNA) is modified at more than 100 locations, particularly at highly conserved and functionally important nucleotides. During ribosome biogenesis, modifications are added at various stages of assembly. The existence of differently modified classes of ribosomes in normal cells is unknown because no method exists to simultaneously evaluate the modification status at all sites within a single rRNA molecule. Using a combination of yeast genetics and nanopore direct RNA sequencing, we developed a reliable method to track the modification status of single rRNA molecules at 37 sites in 18 S rRNA and 73 sites in 25 S rRNA. We use our method to characterize patterns of modification heterogeneity and identify concerted modification of nucleotides found near functional centers of the ribosome. Distinct, undermodified subpopulations of rRNAs accumulate upon loss of Dbp3 or Prp43 RNA helicases, suggesting overlapping roles in ribosome biogenesis. Modification profiles are surprisingly resistant to change in response to many genetic and acute environmental conditions that affect translation, ribosome biogenesis, and pre-mRNA splicing. The ability to capture single-molecule RNA modification profiles provides new insights into the roles of nucleotide modifications in RNA function.

Exploring the expanding universe of small RNAs

The world of small noncoding RNAs (sncRNAs) is ever-expanding, from small interfering RNA, microRNA and Piwi-interacting RNA to the recently emerging non-canonical sncRNAs derived from longer structured RNAs (for example, transfer, ribosomal, Y, small nucleolar, small nuclear and vault RNAs), showing distinct biogenesis and functional principles. Here we discuss recent tools for sncRNA identification, caveats in sncRNA expression analysis and emerging methods for direct sequencing of sncRNAs and systematic mapping of RNA modifications that are integral to their function.

tRNA modification dynamics from individual organisms to metaepitranscriptomics of microbiomes

tRNA is the most extensively modified RNA in cells. On average, a bacterial tRNA contains 8 modifications per molecule and a eukaryotic tRNA contains 13 modifications per molecule. Recent studies reveal that tRNA modifications are highly dynamic and respond extensively to environmental conditions. Functions of tRNA modification dynamics include enhanced, on-demand decoding of specific codons in response genes and regulation of tRNA fragment biogenesis. This review summarizes recent advances in the studies of tRNA modification dynamics in biological processes, tRNA modification erasers, and human-associated bacteria. Furthermore, we use the term "metaepitranscriptomics" to describe the potential and approach of tRNA modification studies in natural biological communities such as microbiomes. tRNA is highly modified in cells, and tRNA modifications respond extensively to environmental conditions to enhance translation of specific genes and produce tRNA fragments on demand. We review recent advances in tRNA sequencing methods, tRNA modification dynamics in biological processes, and tRNA modification studies in natural communities such as the microbiomes.