Quantitative mapping of the cellular small RNA landscape with AQRNA-seq

SnoRNAs in cardiovascular development, function, and disease

Small nucleolar RNAs (snoRNAs) are emerging as important regulators of cardiovascular (patho)biology. Several roles of snoRNAs have recently been identified in heart development and congenital heart diseases, as well as their dynamic regulation in hypertrophic and dilated cardiomyopathies, coronary heart disease (CHD), myocardial infarction (MI), cardiac fibrosis, and heart failure. Furthermore, reports of changes in vesicular snoRNA expression and altered levels of circulating snoRNAs in response to cardiac stress suggest that snoRNAs also function in cardiac signaling and intercellular communication. In this review, we summarize and discuss key findings and outline the clinical potential of snoRNAs considering current challenges and gaps in the field of cardiovascular diseases (CVDs).

Charting new territory: The Plasmodium falciparum tRNA modification landscape

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

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.

Small RNA structural biochemistry in a post-sequencing era

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

Antibody production relies on the tRNA inosine wobble modification to meet biased codon demand

Antibodies are produced at high rates to provide immunoprotection, which puts pressure on the B cell translational machinery. Here, we identified a pattern of codon usage conserved across antibody genes. One feature thereof is the hyperutilization of codons that lack genome-encoded Watson-Crick transfer RNAs (tRNAs), instead relying on the posttranscriptional tRNA modification inosine (I34), which expands the decoding capacity of specific tRNAs through wobbling. Antibody-secreting cells had increased I34 levels and were more reliant on I34 for protein production than naïve B cells. Furthermore, antibody I34-dependent codon usage may influence B cell passage through regulatory checkpoints. Our work elucidates the interface between the tRNA pool and protein production in the immune system and has implications for the design and selection of antibodies for vaccines and therapeutics.

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.

Molecular Coping Mechanisms: Reprogramming tRNAs To Regulate Codon-Biased Translation of Stress Response Proteins

As part of the classic central dogma of molecular biology, transfer RNAs (tRNAs) are integral to protein translation as the adaptor molecules that link the genetic code in messenger RNA (mRNA) to the amino acids in the growing peptide chain. tRNA function is complicated by the existence of 61 codons to specify 20 amino acids, with most amino acids coded by two or more synonymous codons. Further, there are often fewer tRNAs with unique anticodons than there are synonymous codons for an amino acid, with a single anticodon able to decode several codons by "wobbling" of the base pairs arising between the third base of the codon and the first position on the anticodon. The complications introduced by synonymous codons and wobble base pairing began to resolve in the 1960s with the discovery of dozens of chemical modifications of the ribonucleotides in tRNA, which, by analogy to the epigenome, are now collectively referred to as the epitranscriptome for not changing the genetic code inherent to all RNA sequences. tRNA modifications were found to stabilize codon-anticodon interactions, prevent misinitiation of translation, and promote translational fidelity, among other functions, with modification deficiencies causing pathological phenotypes. This led to hypotheses that modification-dependent tRNA decoding efficiencies might play regulatory roles in cells. However, it was only with the advent of systems biology and convergent "omic" technologies that the higher level function of synonymous codons and tRNA modifications began to emerge.Here, we describe our laboratories' discovery of tRNA reprogramming and codon-biased translation as a mechanism linking tRNA modifications and synonymous codon usage to regulation of gene expression at the level of translation. Taking a historical approach, we recount how we discovered that the 8-10 modifications in each tRNA molecule undergo unique reprogramming in response to cellular stresses to promote translation of mRNA transcripts with unique codon usage patterns. These modification tunable transcripts (MoTTs) are enriched with specific codons that are differentially decoded by modified tRNAs and that fall into functional families of genes encoding proteins necessary to survive the specific stress. By developing and applying systems-level technologies, we showed that cells lacking specific tRNA modifications are sensitized to certain cellular stresses by mistranslation of proteins, disruption of mitochondrial function, and failure to translate critical stress response proteins. In essence, tRNA reprogramming serves as a cellular coping strategy, enabling rapid translation of proteins required for stress-specific cell response programs. Notably, this phenomenon has now been characterized in all organisms from viruses to humans and in response to all types of environmental changes. We also elaborate on recent findings that cancer cells hijack this mechanism to promote their own growth, metastasis, and chemotherapeutic resistance. We close by discussing how understanding of codon-biased translation in various systems can be exploited to develop new therapeutics and biomanufacturing processes.

Diverse regulatory pathways modulate bet hedging of competence induction in epigenetically-differentiated phase variants of Streptococcus pneumoniae

Despite enabling Streptococcus pneumoniae to acquire antibiotic resistance and evade vaccine-induced immunity, transformation occurs at variable rates across pneumococci. Phase variants of isolate RMV7, distinguished by altered methylation patterns driven by the translocating variable restriction-modification (tvr) locus, differed significantly in their transformation efficiencies and biofilm thicknesses. These differences were replicated when the corresponding tvr alleles were introduced into an RMV7 derivative lacking the locus. RNA-seq identified differential expression of the type 1 pilus, causing the variation in biofilm formation, and inhibition of competence induction in the less transformable variant, RMV7domi. This was partly attributable to RMV7domi's lower expression of ManLMN, which promoted competence induction through importing N-acetylglucosamine. This effect was potentiated by analogues of some proteobacterial competence regulatory machinery. Additionally, one of RMV7domi's phage-related chromosomal island was relatively active, which inhibited transformation by increasing expression of the stress response proteins ClpP and HrcA. However, HrcA increased competence induction in the other variant, with its effects depending on Ca2+ supplementation and heat shock. Hence the heterogeneity in transformation efficiency likely reflects the diverse signalling pathways by which it is affected. This regulatory complexity will modulate population-wide responses to synchronising quorum sensing signals to produce co-ordinated yet stochastic bet hedging behaviour.

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

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

Dynamic Landscapes of tRNA Transcriptomes and Translatomes in Diverse Mouse Tissues

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

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.

A systematic survey of LU domain-containing proteins reveals a novel human gene, LY6A, which encodes the candidate ortholog of mouse Ly-6A/Sca-1 and is aberrantly expressed in pituitary tumors

The Ly-6 and uPAR (LU) domain-containing proteins represent a large family of cell-surface markers. In particular, mouse Ly-6A/Sca-1 is a widely used marker for various stem cells; however, its human ortholog is missing. In this study, based on a systematic survey and comparative genomic study of mouse and human LU domain-containing proteins, we identified a previously unannotated human gene encoding the candidate ortholog of mouse Ly-6A/Sca-1. This gene, hereby named LY6A, reversely overlaps with a lncRNA gene in the majority of exonic sequences. We found that LY6A is aberrantly expressed in pituitary tumors, but not in normal pituitary tissues, and may contribute to tumorigenesis. Similar to mouse Ly-6A/Sca-1, human LY6A is also upregulated by interferon, suggesting a conserved transcriptional regulatory mechanism between humans and mice. We cloned the full-length LY6A cDNA, whose encoded protein sequence, domain architecture, and exon-intron structures are all well conserved with mouse Ly-6A/Sca-1. Ectopic expression of the LY6A protein in cells demonstrates that it acts the same as mouse Ly-6A/Sca-1 in their processing and glycosylphosphatidylinositol anchoring to the cell membrane. Collectively, these studies unveil a novel human gene encoding a candidate biomarker and provide an interesting model gene for studying gene regulatory and evolutionary mechanisms.

Reassessment of miRNA variant (isomiRs) composition by small RNA sequencing

IsomiRs, sequence variants of mature microRNAs, are usually detected and quantified using high-throughput sequencing. Many examples of their biological relevance have been reported, but sequencing artifacts identified as artificial variants might bias biological inference and therefore need to be ideally avoided. We conducted a comprehensive evaluation of 10 different small RNA sequencing protocols, exploring both a theoretically isomiR-free pool of synthetic miRNAs and HEK293T cells. We calculated that, with the exception of two protocols, less than 5% of miRNA reads can be attributed to library preparation artifacts. Randomized-end adapter protocols showed superior accuracy, with 40% of true biological isomiRs. Nevertheless, we demonstrate concordance across protocols for selected miRNAs in non-templated uridyl additions. Notably, NTA-U calling and isomiR target prediction can be inaccurate when using protocols with poor single-nucleotide resolution. Our results highlight the relevance of protocol choice for biological isomiRs detection and annotation, which has key potential implications for biomedical applications.

Emerging roles of tRNA in cancer

Transfer RNAs (tRNAs) play pivotal roles in the transmission of genetic information, and abnormality of tRNAs directly leads to translation disorders and causes diseases, including cancer. The complex modifications enable tRNA to execute its delicate biological function. Alteration of appropriate modifications may affect the stability of tRNA, impair its ability to carry amino acids, and disrupt the pairing between anticodons and codons. Studies confirmed that dysregulation of tRNA modifications plays an important role in carcinogenesis. Furthermore, when the stability of tRNA is impaired, tRNAs are cleaved into small tRNA fragments (tRFs) by specific RNases. Though tRFs have been found to play vital regulatory roles in tumorigenesis, its formation process is far from clear. Understanding improper tRNA modifications and abnormal formation of tRFs in cancer is conducive to uncovering the role of metabolic process of tRNA under pathological conditions, which may open up new avenues for cancer prevention and treatment.

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

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

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.

Selective Chemical Labeling Strategy for Oligonucleotides Determination: A First Application to Full-Range Profiling of Transfer RNA Modifications

To date, the extremely high polarity and poor signal intensity of macromolecular nucleic acids are greatly impeding the progress of mass spectrometry technology in the quality control of nucleic acid drugs and the characterization of DNA oxidation and RNA modifications. We recently described a general N-(tert-butyldimethylsilyl)-N-methyl-trifluoroacetamide (MTBSTFA) labeling method for oligonucleotide determination and applied it to the full-range profiling of tRNA in vitro and in vivo studies for the first time. The primary advantages of this method include strong retention, no observable byproducts, predictable and easily interpreted MS² data, and the circumvention of instrument harmful reagents that were necessary in previous methods. Selective labeling of N-(tert-butyldimethylsilyl)-N-methyl-trifluoroacetamide to the terminal phosphate groups of oligonucleotides endows it broadly applicable for DNA/RNA profiling. Moreover, the improvement of sequence coverage was achieved in yeast tRNA^phe(GAA) analysis owing to this method's good detection capability of 1-12 nucleotides in length. We also extended this strategy to determine the abundance of modified bases and discover new modifications via digesting RNA into single-nucleotide products, promoting the comprehensive mapping of RNA. The easy availability of derivatization reagent and the simple, rapid one-step reaction render it easy to operate for researchers. When applied in characterizing tRNAs in HepG2 cells and rats with nonalcoholic fatty liver disease, a fragment of U[m¹G][m²G], specific for tRNA^Asn(QUU) in cells, was significantly upregulated, indicating a possible clue to nonalcoholic fatty liver disease pathogenesis.

Human circulating small non-coding RNA signature as a non-invasive biomarker in clinical diagnosis of acute myeloid leukaemia

Background: Acute myeloid leukaemia (AML) is the most common acute leukaemia in adults; AML is highly heterogeneous and involves abnormalities at multiple omics levels. Small non-coding RNAs (sncRNAs) present in body fluids are important regulatory molecules and considered promising non-invasive clinical diagnostic biomarkers for disease. However, the signature of sncRNA profile alteration in AML patient serum and bone marrow supernatant is still under exploration. Methods: We examined data for blood and bone marrow samples from 80 consecutive, newly-diagnosed patients with AML and 12 healthy controls for high throughput small RNA-sequencing. Differentially expressed sncRNAs were analysed to reveal distinct patterns between AML patients and controls. Machine learning methods were used to evaluate the efficiency of specific sncRNAs in discriminating individuals with AML from controls. The altered expression level of individual sncRNAs was evaluated by RT-PCR, Q-PCR, and northern blot. Correlation analysis was employed to assess sncRNA patterns between serum and bone marrow supernatant. Results: We identified over 20 types of sncRNA categories beyond miRNAs in both serum and bone marrow supernatant, with highly coordinated expression patterns between them. Non-classical sncRNAs, including rsRNA (62.86%), ysRNA (14.97%), and tsRNA (4.22%), dominated among serum sncRNAs and showed sensitive alteration patterns in AML patients. According to machine learning-based algorithms, the tsRNA-based signature robustly discriminated subjects with AML from controls and was more reliable than that comprising miRNAs. Our data also showed that serum tsRNAs to be closely associated with AML prognosis, suggesting the potential application of serum tsRNAs as biomarkers to assist in AML diagnosis. Conclusions: We comprehensively characterized the expression pattern of circulating sncRNAs in blood and bone marrow and their alteration signature between healthy controls and AML patients. This study enriches research of sncRNAs in the regulation of AML, and provides insights into the role of sncRNAs in AML.

tRNA Modifications and Modifying Enzymes in Disease, the Potential Therapeutic Targets

tRNA is one of the most conserved and abundant RNA species, which plays a key role during protein translation. tRNA molecules are post-transcriptionally modified by tRNA modifying enzymes. Since high-throughput sequencing technology has developed rapidly, tRNA modification types have been discovered in many research fields. In tRNA, numerous types of tRNA modifications and modifying enzymes have been implicated in biological functions and human diseases. In our review, we talk about the relevant biological functions of tRNA modifications, including tRNA stability, protein translation, cell cycle, oxidative stress, and immunity. We also explore how tRNA modifications contribute to the progression of human diseases. Based on previous studies, we discuss some emerging techniques for assessing tRNA modifications to aid in discovering different types of tRNA modifications.

m7G-quant-seq: Quantitative Detection of RNA Internal N7-Methylguanosine

Methods for the precise detection and quantification of RNA modifications are critical to uncover functional roles of diverse RNA modifications. The internal m⁷G modification in mammalian cytoplasmic tRNAs is known to affect tRNA function and impact embryonic stem cell self-renewal, tumorigenesis, cancer progression, and other cellular processes. Here, we introduce m⁷G-quant-seq, a quantitative method that accurately detects internal m⁷G sites in human cytoplasmic tRNAs at single-base resolution. The efficient chemical reduction and mild depurination can almost completely convert internal m⁷G sites into RNA abasic sites (AP sites). We demonstrate that RNA abasic sites induce a mixed variation pattern during reverse transcription, including G → A or C or T mutations as well as deletions. We calculated the total variation ratio to quantify the m⁷G modification fraction at each methylated site. The calibration curves of all relevant motif contexts allow us to more quantitatively determine the m⁷G methylation level. We detected internal m⁷G sites in 22 human cytoplasmic tRNAs from HeLa and HEK293T cells and successfully estimated the corresponding m⁷G methylation stoichiometry. m⁷G-quant-seq could be applied to monitor the tRNA m⁷G methylation level change in diverse biological processes.

The Transition from Cancer "omics" to "epi-omics" through Next- and Third-Generation Sequencing

Deciphering cancer etiopathogenesis has proven to be an especially challenging task since the mechanisms that drive tumor development and progression are far from simple. An astonishing amount of research has revealed a wide spectrum of defects, including genomic abnormalities, epigenomic alterations, disturbance of gene transcription, as well as post-translational protein modifications, which cooperatively promote carcinogenesis. These findings suggest that the adoption of a multidimensional approach can provide a much more precise and comprehensive picture of the tumor landscape, hence serving as a powerful tool in cancer research and precision oncology. The introduction of next- and third-generation sequencing technologies paved the way for the decoding of genetic information and the elucidation of cancer-related cellular compounds and mechanisms. In the present review, we discuss the current and emerging applications of both generations of sequencing technologies, also referred to as massive parallel sequencing (MPS), in the fields of cancer genomics, transcriptomics and proteomics, as well as in the progressing realms of epi-omics. Finally, we provide a brief insight into the expanding scope of sequencing applications in personalized cancer medicine and pharmacogenomics.

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.

Dysfunctional tRNA reprogramming and codon-biased translation in cancer

Many cancers hijack translation to increase the synthesis of tumor-driving proteins, the messenger mRNAs of which have specific codon usage patterns. Termed 'codon-biased translation' and originally identified in stress response regulation, this mechanism is supported by diverse studies demonstrating how the 50 RNA modifications of the epitranscriptome, specific tRNAs, and codon-biased mRNAs are used by oncogenic programs to promote proliferation and chemoresistance. The epitranscriptome writers METTL1-WDR4, Elongator complex protein (ELP)1-6, CTU1-2, and ALKBH8-TRM112 illustrate the principal mechanism of codon-biased translation, with gene amplifications, increased RNA modifications, and enhanced tRNA stability promoting cancer proliferation. Furthermore, systems-level analyses of 34 tRNA writers and 493 tRNA genes highlight the theme of tRNA epitranscriptome dysregulation in many cancers and identify candidate tRNA writers, tRNA modifications, and tRNA molecules as drivers of pathological codon-biased translation.

Application of nucleoside or nucleotide analogues in RNA dynamics and RNA-binding protein analysis

Cellular RNAs undergo dynamic changes during RNA biological processes, which are tightly orchestrated by RNA-binding proteins (RBPs). Yet, the investigation of RNA dynamics is hurdled by highly abundant steady-state RNAs, which make the signals of dynamic RNAs less detectable. Notably, the exert of nucleoside or nucleotide analogue-based RNA technologies has provided a remarkable platform for RNA dynamics research, revealing diverse unnoticed features in RNA metabolism. In this review, we focus on the application of two types of analogue-based RNA sequencing, antigen-/antibody- and click chemistry-based methodologies, and summarize the RNA dynamics features revealed. Moreover, we discuss emerging single-cell newly transcribed RNA sequencing methodologies based on nucleoside analogue labeling, which provides novel insights into RNA dynamics regulation at single-cell resolution. On the other hand, we also emphasize the identification of RBPs that interact with polyA, non-polyA RNAs, or newly transcribed RNAs and also their associated RNA-binding domains at genomewide level through ultraviolet crosslinking and mass spectrometry in different contexts. We anticipated that further modification and development of these analogue-based RNA and RBP capture technologies will aid in obtaining an unprecedented understanding of RNA biology. This article is categorized under: RNA Interactions with Proteins and Other Molecules > Protein-RNA Recognition RNA Structure and Dynamics > RNA Structure, Dynamics and Chemistry RNA Methods > RNA Analyses in Cells.

Transfer RNA-derived small RNAs (tsRNAs): Versatile regulators in cancer

tRNA-derived small RNAs (tsRNAs) represent a novel class of regulatory small non-coding RNAs (sncRNAs), produced by the specific cleavage of transfer RNAs (tRNAs). In recent years, pilot studies one after the other have uncovered the critical roles of tsRNAs in various fundamental biological processes as well as in the development of human diseases including cancer. Based on the newly updated hallmarks of cancer, we provide a comprehensive review regarding the dysregulation, functional implications and complicated molecular mechanisms of tsRNAs in cancer. In addition, the potential technical challenges and future prospects in the fields of tsRNA research are discussed in this review.

Transfer RNAs-derived small RNAs and their application potential in multiple diseases

The role of tRNAs is best known as adapter components of translational machinery. According to the central dogma of molecular biology, DNA is transcribed to RNA and in turn is translated into proteins, in which tRNA outstands by its role of the cellular courier. Recent studies have led to the revision of the canonical function of transfer RNAs (tRNAs), which indicates that tRNAs also serve as a source for short non-coding RNAs called tRNA-derived small RNAs (tsRNAs). tsRNAs play key roles in cellular processes by modulating complicated regulatory networks beyond translation and are widely involved in multiple diseases. Herein, the biogenesis and classification of tsRNAs were firstly clarified. tsRNAs are generated from pre-tRNAs or mature tRNAs and are classified into tRNA-derived fragments (tRFs) and tRNA halves (tiRNA). The tRFs include five types according to the incision loci: tRF-1, tRF-2, tRF-3, tRF-5 and i-tRF which contain 3′ tiRNA and 5′ tiRNA. The functions of tsRNAs and their regulation mechanisms involved in disease processes are systematically summarized as well. The mechanisms can elaborate on the specific regulation of tsRNAs. In conclusion, the current research suggests that tsRNAs are promising targets for modulating pathological processes, such as breast cancer, ischemic stroke, respiratory syncytial virus, osteoporosis and so on, and maintain vital clinical implications in diagnosis and therapeutics of various diseases.

Small but strong: Pivotal roles and potential applications of snoRNAs in hematopoietic malignancies

Small nucleolar RNAs (snoRNAs) belong to a family of noncoding RNAs that are 60-300 nucleotides in length, and they are classified into two classes according to their structure and function: C/D box snoRNAs, playing an essential role in 2’-O-methylation modification on ribosomal RNA; H/ACA box snoRNAs, involved in the pseudouridylation of rRNA. SnoRNAs with unclear functions, no predictable targets, and unusual subcellular locations are called orphan snoRNAs. Recent studies have revealed abnormal expression and demonstrated the pivotal roles of snoRNAs and their host genes in various types of hematological malignancies. This review discusses recent discoveries concerning snoRNAs in a variety of hematological malignancies, including multiple myeloma, lymphoma and leukemia, and sheds light on the application of snoRNAs as diagnostic and prognostic markers as well as therapeutic targets of hematological malignancies in the future.

The Aminoacyl-tRNA Synthetase and tRNA Expression Levels Are Deregulated in Cancer and Correlate Independently with Patient Survival

Aminoacyl-tRNA synthetases (ARSs) are essential enzymes that load amino acids to their cognate tRNA molecules. The expression of certain ARSs and tRNAs has been shown to be deregulated in cancer, presumably to accommodate elevated protein synthesis requirements. In this work, the expression of cytoplasmic ARSs and tRNAs in ten TCGA cancers has been systematically examined. ARSs were found to be mostly upregulated in tumours and their upregulation often correlated with worse patient survival. tRNAs were found to be either upregulated or downregulated in tumours and their expression sometimes correlated to worse survival outcomes. However, although the expression of most ARSs and tRNAs was deregulated in tumours when compared to healthy adjacent tissues, only in a few cases, and independently, did it correlate to patient survival. These data point to the general uncoupling of concomitant ARS and tRNA expression deregulation and patient survival, highlighting the different ARS and tRNA requirements in cancers.

tRNAstudio: facilitating the study of human mature tRNAs from deep sequencing datasets

SUMMARY

High-throughput sequencing of transfer RNAs (tRNA-Seq) is a powerful approach to characterize the cellular tRNA pool. Currently, however, analyzing tRNA-Seq datasets requires strong bioinformatics and programming skills. tRNAstudio facilitates the analysis of tRNA-Seq datasets and extracts information on tRNA gene expression, post-transcriptional tRNA modification levels, and tRNA processing steps. Users need only running a few simple bash commands to activate a graphical user interface that allows the easy processing of tRNA-Seq datasets in local mode. Output files include extensive graphical representations and associated numerical tables, and an interactive html summary report to help interpret the data. We have validated tRNAstudio using datasets generated by different experimental methods and derived from human cell lines and tissues that present distinct patterns of tRNA expression, modification and processing.

AVAILABILITY AND IMPLEMENTATION

Freely available at https://github.com/GeneTranslationLab-IRB/tRNAstudio under an open-source GNU GPL v3.0 license.

SUPPLEMENTARY INFORMATION

Supplementary data are available at Bioinformatics online.

RNA marker modifications reveal the necessity for rigorous preparation protocols to avoid artifacts in epitranscriptomic analysis

The accurate definition of an epitranscriptome is endangered by artefacts resulting from RNA degradation after cell death, a ubiquitous yet little investigated process. By tracing RNA marker modifications through tissue preparation protocols, we identified a major blind spot from daily lab routine, that has massive impact on modification analysis in small RNAs. In particular, m6,6A and Am as co-varying rRNA marker modifications, appeared in small RNA fractions following rRNA degradation in vitro and in cellulo. Analysing mouse tissue at different time points post mortem, we tracked the progress of intracellular RNA degradation after cell death, and found it reflected in RNA modification patterns. Differences were dramatic between liver, where RNA degradation commenced immediately after death, and brain, yielding essentially undamaged RNA. RNA integrity correlated with low amounts of co-varying rRNA markers. Thus validated RNA preparations featured differentially modified tRNA populations whose information content allowed a distinction even among the related brain tissues cortex, cerebellum and hippocampus. Inversely, advanced cell death correlated with high rRNA marker content, and correspondingly little with the naïve state of living tissue. Therefore, unless RNA and tissue preparations are executed with utmost care, interpretation of modification patterns in tRNA and small RNA are prone to artefacts.

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.

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.

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.

The epitranscriptome of small non-coding RNAs

Small non-coding RNAs are short RNA molecules and involved in many biological processes, including cell proliferation and differentiation, immune response, cell death, epigenetic regulation, metabolic control. A diversity of RNA modifications have been identified in these small non-coding RNAs, including transfer RNAs (tRNAs), microRNAs (miRNAs), PIWI-interacting RNAs (piRNAs), small nuclear RNA (snRNA), small nucleolar RNAs (snoRNAs), and tRNA-derived small RNAs (tsRNAs). These post-transcriptional modifications are involved in the biogenesis and function of these small non-coding RNAs. In this review, we will summarize the existence of RNA modifications in the small non-coding RNAs and the emerging roles of these epitranscriptomic marks.

5' Half of specific tRNAs feeds back to promote corresponding tRNA gene transcription in vertebrate embryos

5′tRFls are small transfer RNA (tRNA) fragments derived from 5′ half of mature tRNAs. However, it is unknown whether 5′tRFls could feed back to regulate tRNA biogenesis. Here, we show that 5′tRFl^Gly/GCC and 5′tRFl^Glu/CTC function to promote transcription of corresponding tRNA genes and are essential for vertebrate early embryogenesis. During zebrafish embryogenesis, dynamics of 5′tRFl^Gly/GCC and 5′tRFl^Glu/CTC levels correlates with that of tRNA^Gly/GCC and tRNA^Glu/CTC levels. Morpholino-mediated knockdown of 5′tRFl^Gly/GCC or 5′tRFl^Glu/CTC down-regulates tRNA^Gly/GCC or tRNA^Glu/CTC levels, respectively, and causes embryonic lethality that is efficiently rescued by coinjection of properly refolded corresponding tRNA. In zebrafish embryos, tRNA:DNA and 5′tRFl:DNA hybrids commonly exist on the template strand of tRNA genes. Mechanistically, unstable 5′tRFl:DNA hybrid may prevent the formation of transcriptionally inhibitory stable tRNA:DNA hybrids on the same tRNA loci so as to facilitate tRNA gene transcription. The uncovered mechanism may be implicated in other physiological and pathological processes.

Plant transfer RNA-derived fragments: Biogenesis and functions

Processing of mature transfer RNAs (tRNAs) produces complex populations of tRNA-derived fragments (tRFs). Emerging evidence shows that tRFs have important functions in bacteria, animals, and plants. Here, we review recent advances in understanding plant tRFs, focusing on their biological and cellular functions, such as regulating stress responses, mediating plant-pathogen interactions, and modulating post-transcriptional gene silencing and translation. We also review sequencing strategies and bioinformatics resources for studying tRFs in plants. Finally, we discuss future directions for plant tRF research, which will expand our knowledge of plant non-coding RNAs.

Aberration-corrected ultrafine analysis of miRNA reads at single-base resolution: a k-mer lattice approach

Raw sequencing reads of miRNAs contain machine-made substitution errors, or even insertions and deletions (indels). Although the error rate can be low at 0.1%, precise rectification of these errors is critically important because isoform variation analysis at single-base resolution such as novel isomiR discovery, editing events understanding, differential expression analysis, or tissue-specific isoform identification is very sensitive to base positions and copy counts of the reads. Existing error correction methods do not work for miRNA sequencing data attributed to miRNAs’ length and per-read-coverage properties distinct from DNA or mRNA sequencing reads. We present a novel lattice structure combining kmers, (k – 1)mers and (k + 1)mers to address this problem. The method is particularly effective for the correction of indel errors. Extensive tests on datasets having known ground truth of errors demonstrate that the method is able to remove almost all of the errors, without introducing any new error, to improve the data quality from every-50-reads containing one error to every-1300-reads containing one error. Studies on experimental miRNA sequencing datasets show that the errors are often rectified at the 5′ ends and the seed regions of the reads, and that there are remarkable changes after the correction in miRNA isoform abundance, volume of singleton reads, overall entropy, isomiR families, tissue-specific miRNAs, and rare-miRNA quantities.