Next-generation sequencing technologies for detection of modified nucleotides in RNAs

Implications of RNA pseudouridylation for cancer biology and therapeutics: a narrative review

BACKGROUND

Pseudouridine (Ψ), a C5-glycoside isomer of uridine, stands as one of the most prevalent RNA modifications in all RNA types. Distinguishing from the C-N bond linking uridine to ribose, the link between Ψ and ribose is a C-C bond, endowing Ψ modified RNA distinct properties and functions in various biological processes. The conversion of uridine to Ψ is governed by pseudouridine synthases (PUSs). RNA pseudouridylation is implicated in cancer biology and therapeutics.

OBJECTIVES

In this review, we will summarize the methods for detecting Ψ, the process of Ψ generation, the impact of Ψ modification on RNA metabolism and gene expression, the roles of dysregulated Ψ and pseudouridine synthases in cancers, and the underlying mechanism.

METHODS

We conducted a comprehensive search of PubMed from its inception through February 2024. The search terms included "pseudouridine"; "pseudouridine synthase"; "PUS"; "dyskerin"; "cancer"; "tumor"; "carcinoma"; "malignancy"; "tumorigenesis"; "biomarker"; "prognosis" and "therapy". We included studies published in peer-reviewed journals that focused on Ψ detection, specific mechanisms involving Ψ and PUSs, and prognosis in cancer patients with high Ψ expression. We excluded studies lacking sufficient methodological details or appropriate controls.

RESULTS

Ψ has been recognized as a significant biomarker in cancer diagnosis and prognosis. Abnormal Ψ modifications mediated by various PUSs result in dysregulated RNA metabolism and impaired RNA function, promoting the development of various cancers. Overexpression of PUSs is common in cancer cells and predicts poor prognosis. PUSs inhibition arrests cell proliferation and enhances apoptosis in cancer cells, suggesting PUS-targeting cancer therapy may be a potential strategy in cancer treatment.

DISCUSSION

High Ψ levels in serum, urine, and saliva may suggest cancer, but do not specify the type, requiring additional lab markers and imaging for accurate diagnosis. Standardized detection methods are also crucial for reliable results. PUSs are linked to cancer, but more researches are needed to understand their mechanisms in different cancers. Anticancer treatments targeting PUSs are still under developed.

Role of RNA modifications in blood development and regeneration

Blood development and regeneration require rapid turnover of cells, and ribonucleic acid (RNA) modifications play a key role in it via regulating stemness and cell fate regulation. RNA modifications affect gene activity via posttranscriptional and translation-mediated mechanisms. Diverse molecular players involved in RNA-modification processes are abundantly expressed by hematopoietic stem cells and lineages. Close to 150 RNA chemical modifications have been reported, but only N6-methyl adenosine (m6A), inosine (I), pseudouridine (Ψ), and m1A-a handful-have been studied in-cell fate regulation. The role of RNA modification in blood diseases and disorders is an emerging field and offers potential for therapeutic interventions. Knowledge of RNA-modification and enzymatic activities could be used to design therapies in the future. Here, we summarized the recent advances in RNA modification and the epitranscriptome field and discussed their regulation of blood development and regeneration.

Analysis of bacterial transcriptome and epitranscriptome using nanopore direct RNA sequencing

Bacterial gene expression is a complex process involving extensive regulatory mechanisms. Along with growing interests in this field, Nanopore Direct RNA Sequencing (DRS) provides a promising platform for rapid and comprehensive characterization of bacterial RNA biology. However, the DRS of bacterial RNA is currently deficient in the yield of mRNA-mapping reads and has yet to be exploited for transcriptome-wide RNA modification mapping. Here, we showed that pre-processing of bacterial total RNA (size selection followed by ribosomal RNA depletion and polyadenylation) guaranteed high throughputs of sequencing data and considerably increased the amount of mRNA reads. This way, complex transcriptome architectures were reconstructed for Escherichia coli and Staphylococcus aureus and extended the boundaries of 225 known E. coli operons and 89 defined S. aureus operons. Utilizing unmodified in vitro-transcribed (IVT) RNA libraries as a negative control, several Nanopore-based computational tools globally detected putative modification sites in the E. coli and S. aureus transcriptomes. Combined with Next-Generation Sequencing-based N6-methyladenosine (m6A) detection methods, 75 high-confidence m6A candidates were identified in the E. coli protein-coding transcripts, while none were detected in S. aureus. Altogether, we demonstrated the potential of Nanopore DRS in systematic and convenient transcriptome and epitranscriptome analysis.

Temperature-Dependent tRNA Modifications in Bacillales

Transfer RNA (tRNA) modifications are essential for the temperature adaptation of thermophilic and psychrophilic organisms as they control the rigidity and flexibility of transcripts. To further understand how specific tRNA modifications are adjusted to maintain functionality in response to temperature fluctuations, we investigated whether tRNA modifications represent an adaptation of bacteria to different growth temperatures (minimal, optimal, and maximal), focusing on closely related psychrophilic (P. halocryophilus and E. sibiricum), mesophilic (B. subtilis), and thermophilic (G. stearothermophilus) Bacillales. Utilizing an RNA sequencing approach combined with chemical pre-treatment of tRNA samples, we systematically profiled dihydrouridine (D), 4-thiouridine (s4U), 7-methyl-guanosine (m7G), and pseudouridine (Ψ) modifications at single-nucleotide resolution. Despite their close relationship, each bacterium exhibited a unique tRNA modification profile. Our findings revealed increased tRNA modifications in the thermophilic bacterium at its optimal growth temperature, particularly showing elevated levels of s4U8 and Ψ55 modifications compared to non-thermophilic bacteria, indicating a temperature-dependent regulation that may contribute to thermotolerance. Furthermore, we observed higher levels of D modifications in psychrophilic and mesophilic bacteria, indicating an adaptive strategy for cold environments by enhancing local flexibility in tRNAs. Our method demonstrated high effectiveness in identifying tRNA modifications compared to an established tool, highlighting its potential for precise tRNA profiling studies.

N6-Methyladenosine Methylation of mRNA in Cell Apoptosis

Apoptosis, a highly controlled homeostatic mechanism that eliminates single cells without destroying tissue function, occurs during growing development and senescence. N6-methyladenosine (m6A), as the most common internal modification of eukaryotic mRNA, fine-tunes gene expression by regulating many aspects of mRNA metabolism, such as splicing, nucleation, stability, translation, and degradation. Remarkably, recent reports have indicated that aberrant methylation of m6A-related RNA may directly or indirectly influence the expression of apoptosis-related genes, thus regulating the process of cell apoptosis. In this review, we summarized the relationship between m6A modification and cell apoptosis, especially its role in the nervous system, and analyzed the limitations of the current research.

Decoding the ribosome's hidden language: rRNA modifications as key players in cancer dynamics and targeted therapies

Ribosomal RNA (rRNA) modifications, essential components of ribosome structure and function, significantly impact cellular proteomics and cancer biology. These chemical modifications transcend structural roles, critically shaping ribosome functionality and influencing cellular protein profiles. In this review, the mechanisms by which rRNA modifications regulate both rRNA functions and broader cellular physiological processes are critically discussed. Importantly, by altering the translational output, rRNA modifications can shift the cellular equilibrium towards oncogenesis, thus playing a key role in cancer development and progression. Moreover, a special focus is placed on the functions of mitochondrial rRNA modifications and their aberrant expression in cancer, an area with profound implications yet largely uncharted. Dysregulation in these modifications can lead to metabolic dysfunction and apoptosis resistance, hallmark traits of cancer cells. Furthermore, the current challenges and future perspectives in targeting rRNA modifications are highlighted as a therapeutic approach for cancer treatment. In conclusion, rRNA modifications represent a frontier in cancer research, offering novel insights and therapeutic possibilities. Understanding and harnessing these modifications can pave the way for breakthroughs in cancer treatment, potentially transforming the approach to combating this complex disease.

RNA modifications in pulmonary diseases

Threatening public health, pulmonary disease (PD) encompasses diverse lung injuries like chronic obstructive PD, pulmonary fibrosis, asthma, pulmonary infections due to pathogen invasion, and fatal lung cancer. The crucial involvement of RNA epigenetic modifications in PD pathogenesis is underscored by robust evidence. These modifications not only shape cell fates but also finely modulate the expression of genes linked to disease progression, suggesting their utility as biomarkers and targets for therapeutic strategies. The critical RNA modifications implicated in PDs are summarized in this review, including N6-methylation of adenosine, N1-methylation of adenosine, 5-methylcytosine, pseudouridine (5-ribosyl uracil), 7-methylguanosine, and adenosine to inosine editing, along with relevant regulatory mechanisms. By shedding light on the pathology of PDs, these summaries could spur the identification of new biomarkers and therapeutic strategies, ultimately paving the way for early PD diagnosis and treatment innovation.

Clinical significance of RNA methylation in hepatocellular carcinoma

Hepatocellular carcinoma (HCC) is a primary liver malignancy with high mortality rates and poor prognosis. Recent advances in high-throughput sequencing and bioinformatic technologies have greatly enhanced the understanding of the genetic and epigenetic changes in liver cancer. Among these changes, RNA methylation, the most prevalent internal RNA modification, has emerged as a significant contributor of the development and progression of HCC. Growing evidence has reported significantly abnormal levels of RNA methylation and dysregulation of RNA-methylation-related enzymes in HCC tissues and cell lines. These alterations in RNA methylation play a crucial role in the regulation of various genes and signaling pathways involved in HCC, thereby promoting tumor progression. Understanding the pathogenesis of RNA methylation in HCC would help in developing prognostic biomarkers and targeted therapies for HCC. Targeting RNA-methylation-related molecules has shown promising potential in the management of HCC, in terms of developing novel prognostic biomarkers and therapies for HCC. Exploring the clinical application of targeted RNA methylation may provide new insights and approaches for the management of HCC. Further research in this field is warranted to fully understand the functional roles and underlying mechanisms of RNA methylation in HCC. In this review, we described the multifaceted functional roles and potential mechanisms of RNA methylation in HCC. Moreover, the prospects of clinical application of targeted RNA methylation for HCC management are discussed, which may provide the basis for subsequent in-depth research on RNA methylation in HCC.

An Overview of Current Detection Methods for RNA Methylation

Epitranscriptomic mechanisms, which constitute an important layer in post-transcriptional gene regulation, are involved in numerous cellular processes under health and disease such as stem cell development or cancer. Among various such mechanisms, RNA methylation is considered to have vital roles in eukaryotes primarily due to its dynamic and reversible nature. There are numerous RNA methylations that include, but are not limited to, 2'-O-dimethyladenosine (m6Am), N7-methylguanosine (m7G), N6-methyladenosine (m6A) and N1-methyladenosine (m1A). These biochemical modifications modulate the fate of RNA by affecting the processes such as translation, target site determination, RNA processing, polyadenylation, splicing, structure, editing and stability. Thus, it is highly important to quantitatively measure the changes in RNA methylation marks to gain insight into cellular processes under health and disease. Although there are complicating challenges in identifying certain methylation marks genome wide, various methods have been developed recently to facilitate the quantitative measurement of methylated RNAs. To this end, the detection methods for RNA methylation can be classified in five categories such as antibody-based, digestion-based, ligation-based, hybridization-based or direct RNA-based methods. In this review, we have aimed to summarize our current understanding of the detection methods for RNA methylation, highlighting their advantages and disadvantages, along with the current challenges in the field.

Research progress of RNA pseudouridine modification in nervous system

Recent advances of pseudouridine (Ψ, 5-ribosyluracil) modification highlight its crucial role as a post-transcriptional regulator in gene expression and its impact on various RNA processes. Ψ synthase (PUS), a category of RNA-modifying enzymes, orchestrates the pseudouridylation reaction. It can specifically recognize conserved sequences or structural motifs within substrates, thereby regulating the biological function of various RNA molecules accurately. Our comprehensive review underscored the close association of PUS1, PUS3, PUS7, PUS10, and dyskerin PUS1 with various nervous system disorders, including neurodevelopmental disorders, nervous system tumors, mitochondrial myopathy, lactic acidosis and sideroblastic anaemia (MLASA) syndrome, peripheral nervous system disorders, and type II myotonic dystrophy. In light of these findings, this study elucidated how Ψ strengthened RNA structures and contributed to RNA function, thereby providing valuable insights into the intricate molecular mechanisms underlying nervous system diseases. However, the detailed effects and mechanisms of PUS on neuron remain elusive. This lack of mechanistic understanding poses a substantial obstacle to the development of therapeutic approaches for various neurological disorders based on Ψ modification.

Analyzing RNA posttranscriptional modifications to decipher the epitranscriptomic code

The discovery of RNA silencing has revealed that non-protein-coding sequences (ncRNAs) can cover essential roles in regulatory networks and their malfunction may result in severe consequences on human health. These findings have prompted a general reassessment of the significance of RNA as a key player in cellular processes. This reassessment, however, will not be complete without a greater understanding of the distribution and function of the over 170 variants of the canonical ribonucleotides, which contribute to the breathtaking structural diversity of natural RNA. This review surveys the analytical approaches employed for the identification, characterization, and detection of RNA posttranscriptional modifications (rPTMs). The merits of analyzing individual units after exhaustive hydrolysis of the initial biopolymer are outlined together with those of identifying their position in the sequence of parent strands. Approaches based on next generation sequencing and mass spectrometry technologies are covered in depth to provide a comprehensive view of their respective merits. Deciphering the epitranscriptomic code will require not only mapping the location of rPTMs in the various classes of RNAs, but also assessing the variations of expression levels under different experimental conditions. The fact that no individual platform is currently capable of meeting all such demands implies that it will be essential to capitalize on complementary approaches to obtain the desired information. For this reason, the review strived to cover the broadest possible range of techniques to provide readers with the fundamental elements necessary to make informed choices and design the most effective possible strategy to accomplish the task at hand.

SWAMNA: a comprehensive platform for analysis of nucleic acid modifications

The interest in MS-based analysis of modified nucleic acids is increasing due to the application of nucleic acids in therapeutics. However, there are few available integrated platforms for characterizing nucleic acid modifications. Herein, we report a general mass spectrometry-based SWATH platform to identify and quantify both RNA and DNA modifications, which we call SWATH analysis of modified nucleic acids (SWAMNA). SWAMNA incorporates the search engine, NuMo finder, enabling the analysis of modifications in native and permethylated form. SWAMNA will aid discoveries that provide new insights into nucleic acid modifications.

Novel insight into RNA modifications in tumor immunity: Promising targets to prevent tumor immune escape

An immunosuppressive state is a typical feature of the tumor microenvironment. Despite the dramatic success of immune checkpoint inhibitor (ICI) therapy in preventing tumor cell escape from immune surveillance, primary and acquired resistance have limited its clinical use. Notably, recent clinical trials have shown that epigenetic drugs can significantly improve the outcome of ICI therapy in various cancers, indicating the importance of epigenetic modifications in immune regulation of tumors. Recently, RNA modifications (N⁶-methyladenosine [m⁶A], N¹-methyladenosine [m¹A], 5-methylcytosine [m⁵C], etc.), novel hotspot areas of epigenetic research, have been shown to play crucial roles in protumor and antitumor immunity. In this review, we provide a comprehensive understanding of how m⁶A, m¹A, and m⁵C function in tumor immunity by directly regulating different immune cells as well as indirectly regulating tumor cells through different mechanisms, including modulating the expression of immune checkpoints, inducing metabolic reprogramming, and affecting the secretion of immune-related factors. Finally, we discuss the current status of strategies targeting RNA modifications to prevent tumor immune escape, highlighting their potential.

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

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

Occurrence of RNA post-transcriptional modifications in plant viruses and viroids and their correlation with structural and functional features

Post-transcriptional modifications of RNA bases are widespread across all the tree of life and have been linked to RNA maturation, stability, and molecular interactions. RNA modifications have been extensively described in endogenous eukaryotic mRNAs, however, little is known about the presence of RNA modifications in plant viral and subviral RNAs. Here, we used a computational approach to infer RNA modifications in plant-pathogenic viruses and viroids using high-throughput annotation of modified ribonucleotides (HAMR), a software that predicts modified ribonucleotides using high-throughput RNA sequencing data. We analyzed datasets from representative members of different plant viruses and viroids and compared them to plant-endogenous mRNAs. Our approach was able to predict potential RNA chemical modifications (RCMs) in all analyzed pathogens. We found that both DNA and RNA viruses presented a wide range of RCM proportions while viroids had lowest values. Furthermore, we found that for viruses with segmented genomes, some genomic RNAs had a higher proportion of RCM. Interestingly, nuclear-replicating viroids showed most of the predicted modifications located in the pathogenesis region, pointing towards a possible functional role of RCMs in their infectious cycle. Thus, our results strongly suggest that plant viral and subviral RNAs might contain a variety of previously unreported RNA modifications, thus opening a new perspective in the multifaceted process of plant-pathogen interactions.

N6-Methyladenosine Methylation of mRNA in Cell Senescence

Cell senescence is the growth arrest caused by the accumulation of irreparable cell damage, which is involved in physiological and pathological processes and regulated by the post-transcriptional level. This regulation is performed by transcriptional regulators and driven by aging-related small RNAs, long non-coding RNAs, and RNA-binding proteins. N6-methyladenosine (m⁶A) is the most common chemical modification in eukaryotic mRNA, which can enhance or reduce the binding of transcriptional regulators. Increasing studies have confirmed the crucial role of m⁶A in controlling mRNA in various physiological processes. Remarkably, recent reports have indicated that abnormal methylation of m⁶A-related RNA may affect cell senescence. In this review, we clarified the association between m⁶A modification and cell senescence and analyzed the limitations of the current research.

Identification and analysis of putative tRNA genes in baculovirus genomes

Transfer RNAs (tRNAs) genes are both coded for and arranged along some viral genomes representing the entire virosphere and seem to play different biological functions during infection, other than transferring the correct amino acid to a growing peptide chain. Baculovirus genome description and annotation has focused mostly on protein-coding genes, microRNA, and homologous regions. Here we carried out a large-scale in silico search for putative tRNA genes in baculovirus genomes. Ninety-six of 257 baculovirus genomes analyzed was found to contain at least one putative tRNA gene. We found great diversity in primary and secondary structure, in location within the genome, in intron presence and size, and in anti-codon identity. In some cases, genes of tRNA-containing genomes were found to have a bias for the codons specified by the tRNAs present in such genomes. Moreover, analysis revealed that most of the putative tRNA genes possessed conserved motifs for tRNA type 2 promoters, including the A-box and B-box motifs with few mismatches from the eukaryotic canonical motifs. From publicly available small RNA deep sequencing datasets of baculovirus-infected insect cells, we found evidence that a putative Autographa californica multiple nucleopolyhedrovirus Gln-tRNA gene was transcribed and modified with the addition of the non-templated 3'-CCA tail found at the end of all tRNAs. Further research is needed to determine the expression and functionality of these viral tRNAs.

The Role of RNA Modification in HIV-1 Infection

RNA plays an important role in biology, and more than 170 RNA modifications have been identified so far. Post-transcriptional modification of RNA in cells plays a crucial role in the regulation of its stability, transport, processing, and gene expression. So far, the research on RNA modification and the exact role of its enzymes is becoming more and more comprehensive. Human immunodeficiency virus 1 (HIV-1) is an RNA virus and the causative agent of acquired immunodeficiency syndrome (AIDS), which is one of the most devastating viral pandemics in history. More and more studies have shown that HIV has RNA modifications and regulation of its gene expression during infection and replication. This review focuses on several RNA modifications and their regulatory roles as well as the roles that different RNA modifications play during HIV-1 infection, in order to find new approaches for the development of anti-HIV-1 therapeutics.

No evidence for age-related differences in mitochondrial RNA quality in the female germline

Abstract

Mitochondrial quality is implicated as a contributor to declining fertility with aging. We investigated mitochondrial transcripts in oocytes and their associated cumulus cells from mice of different ages using RNA-seq. Mice aged 3 weeks, 9 weeks, and 1 year were superovulated, and 48 h later, oocyte cumulus complexes were collected by follicle puncture. We did not detect any major differences that could be attributed to aging. However, mitochondrial RNA transcripts which deviated from the consensus sequence were found at a higher frequency in cumulus cells than in their corresponding oocyte. Previous investigations have shown that variation in the sequence of mtRNA transcripts is substantial, and at least some of this can be accounted for by post-transcriptional modifications which impact base calling during sequencing. Our data would be consistent with either less post-transcriptional modification in mitochondrial RNA from oocytes than cumulus cells or with lower mtDNA mutational load.

Lay summary

Women become less fertile as they age. Shortage of energy contributes to this, caused by a decline in the quality of mitochondria (the powerhouses of the cell) in the egg. Genes are the blueprint for the cell. They are made of DNA which is copied into an RNA message, or instructions, for making proteins. We counted differences in the RNA message of developing eggs and the cells that support them during development (cumulus cells). We compared the number of these differences in mice of different ages. These age groups represent mice had not reached puberty, those of prime reproductive age, and old mothers. We did not find any differences linked to the age of the mice. However, we did find differences between the egg and the cumulus cells. In most cases, there were lower levels of mutations in eggs than there were in cumulus cells.

Chemical modifications to mRNA nucleobases impact translation elongation and termination

Messenger RNAs (mRNAs) serve as blueprints for protein synthesis by the molecular machine the ribosome. The ribosome relies on hydrogen bonding interactions between adaptor aminoacyl-transfer RNA molecules and mRNAs to ensure the rapid and faithful translation of the genetic code into protein. There is a growing body of evidence suggesting that chemical modifications to mRNA nucleosides impact the speed and accuracy of protein synthesis by the ribosome. Modulations in translation rates have downstream effects beyond protein production, influencing protein folding and mRNA stability. Given the prevalence of such modifications in mRNA coding regions, it is imperative to understand the consequences of individual modifications on translation. In this review we present the current state of our knowledge regarding how individual mRNA modifications influence ribosome function. Our comprehensive comparison of the impacts of 16 different mRNA modifications on translation reveals that most modifications can alter the elongation step in the protein synthesis pathway. Additionally, we discuss the context dependence of these effects, highlighting the necessity of further study to uncover the rules that govern how any given chemical modification in an mRNA codon is read by the ribosome.

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.

A bifunctional chemical signature enabling RNA 4-thiouridine enrichment sequencing with single-base resolution

Both sequence enrichment and base resolution are essential for accurate sequencing analysis of low-abundance RNA. Yet they are hindered by the lack of molecular tools. Here we report a bifunctional chemical signature for RNA 4-thiouridine (4sU) enrichment sequencing with single-base resolution. This chemical signature is designed for specific 4sU labeling with two functional parts. One part is a distal alkynyl group for the biotin-assisted pulldown enrichment of target molecules via click chemistry crosslinking. The other part is a -NH group proximal to the pyrimidine ring of 4sU. It allows 4sU-to-cytosine transition during the polymerase-catalyzed extension reaction based on altering hydrogen-bonding patterns. Ultimately, the 4sU-containing RNA molecules can be enriched and accurately analyzed by single-base resolution sequencing. The proposed method also holds great potential to investigate transcriptome dynamics integrated with high-throughput sequencing.

Early transcriptomic response of the mycoparasite Sphaerodes mycoparasitica to the mycotoxigenic Fusarium graminearum 3-ADON, the cause of Fusarium head blight

Mycoparasites are an assemblage of biotrophic and necrotrophic fungi that occur on plant pathogenic fungal hosts. Biotrophic mycoparasites are often overlooked in transcriptomic-based biocontrol studies. Sphaerodes mycoparasitica (S.m.) is a specific biotrophic mycoparasite of plant pathogenic Fusarium graminearum (F.g.), a devastating Fusarium head blight (FHB) disease in small-grain cereals. To understand the biotrophic mycoparasitism comprehensively, we performed Illumina RNA-Seq transcriptomic study on the fungus–fungus interaction in vitro. The aim is to identify the transcript-level mechanism related to the biotrophic S.m. mycoparasitism, particularly its ability to effectively control the F.g. 3-ADON chemotype. A shift in the transcriptomic profile of the mycoparasite was triggered in response to its interaction with F.g. during recognition (1.5 days) and colonization (3.5 days) steps. RNA-Seq analysis revealed ~ 30% of annotated transcripts with "function unknown". Further, 14 differentially expressed genes functionally linked to the biotrophic mycoparasitism were validated by quantitative real-time PCR (qPCR). The gene expression patterns of the filamentous haemagglutinin/adhesin/attachment factor as well as cell wall-degrading glucanases and chitinases were upregulated by host interaction. Besides, mycoparasitism-associated antioxidant resistance genes encoding ATP-binding cassette (ABC) transporter(s) and glutathione synthetase(s) were upregulated. However, the thioredoxin reductase was downregulated which infers that this antioxidant gene can be used as a resistance marker to assess S.m. antifungal and antimycotoxigenic activities. The interactive transcriptome of S. mycoparasitica provides new insights into specific mycoparasitism and will contribute to future research in controlling FHB.

The role of OncoSnoRNAs and Ribosomal RNA 2'-O-methylation in Cancer

Ribosomes are essential nanomachines responsible for all protein production in cells. Ribosome biogenesis and function are energy costly processes, they are tightly regulated to match cellular needs. In cancer, major pathways that control ribosome biogenesis and function are often deregulated to ensure cell survival and to accommodate the continuous proliferation of tumour cells. Ribosomal RNAs (rRNAs) are abundantly modified with 2'-O-methylation (Nm, ribomethylation) being one of the most common modifications. In eukaryotic ribosomes, ribomethylation is performed by the methyltransferase Fibrillarin guided by box C/D small nucleolar RNAs (snoRNAs). Accumulating evidences indicate that snoRNA expression and ribosome methylation profiles are altered in cancer. Here we review our current knowledge on differential snoRNA expression and rRNA 2ʹ-O methylation in the context of human malignancies, and discuss the consequences and opportunities for cancer diagnostics, prognostics, and therapeutics.

Systematic calibration of epitranscriptomic maps using a synthetic modification-free RNA library

Recent years have witnessed rapid progress in the field of epitranscriptomics. Functional interpretation of the epitranscriptome relies on sequencing technologies that determine the location and stoichiometry of various RNA modifications. However, contradictory results have been reported among studies, bringing the biological impacts of certain RNA modifications into doubt. Here, we develop a synthetic RNA library resembling the endogenous transcriptome but without any RNA modification. By incorporating this modification-free RNA library into established mapping techniques as a negative control, we reveal abundant false positives resulting from sequence bias or RNA structure. After calibration, precise and quantitative mapping expands the understanding of two representative modification types, N6-methyladenosine (m6A) and 5-methylcytosine (m5C). We propose that this approach provides a systematic solution for the calibration of various RNA-modification mappings and holds great promise in epitranscriptomic studies.

Chemical capping improves template switching and enhances sequencing of small RNAs

Template-switching reverse transcription is widely used in RNA sequencing for low-input and low-quality samples, including RNA from single cells or formalin-fixed paraffin-embedded (FFPE) tissues. Previously, we identified the native eukaryotic mRNA 5′ cap as a key structural element for enhancing template switching efficiency. Here, we introduce CapTS-seq, a new strategy for sequencing small RNAs that combines chemical capping and template switching. We probed a variety of non-native synthetic cap structures and found that an unmethylated guanosine triphosphate cap led to the lowest bias and highest efficiency for template switching. Through cross-examination of different nucleotides at the cap position, our data provided unequivocal evidence that the 5′ cap acts as a template for the first nucleotide in reverse transcriptase-mediated post-templated addition to the emerging cDNA—a key feature to propel template switching. We deployed CapTS-seq for sequencing synthetic miRNAs, human total brain and liver FFPE RNA, and demonstrated that it consistently improves library quality for miRNAs in comparison with a gold standard template switching-based small RNA-seq kit.

General Principles for the Detection of Modified Nucleotides in RNA by Specific Reagents

Epitranscriptomics heavily rely on chemical reagents for the detection, quantification, and localization of modified nucleotides in transcriptomes. Recent years have seen a surge in mapping methods that use innovative and rediscovered organic chemistry in high throughput approaches. While this has brought about a leap of progress in this young field, it has also become clear that the different chemistries feature variegated specificity and selectivity. The associated error rates, e.g., in terms of false positives and false negatives, are in large part inherent to the chemistry employed. This means that even assuming technically perfect execution, the interpretation of mapping results issuing from the application of such chemistries are limited by intrinsic features of chemical reactivity. An important but often ignored fact is that the huge stochiometric excess of unmodified over-modified nucleotides is not inert to any of the reagents employed. Consequently, any reaction aimed at chemical discrimination of modified versus unmodified nucleotides has optimal conditions for selectivity that are ultimately anchored in relative reaction rates, whose ratio imposes intrinsic limits to selectivity. Here chemical reactivities of canonical and modified ribonucleosides are revisited as a basis for an understanding of the limits of selectivity achievable with chemical methods.

Quantifying RNA modifications by mass spectrometry: a novel source of biomarkers in oncology

Despite significant progress in targeted therapies, cancer recurrence remains a major cause of mortality worldwide. Identification of accurate biomarkers, through molecular profiling in healthy and cancer patient samples, will improve diagnosis and promote personalized medicine. While genetic and epigenetic alterations of DNA are currently exploited as cancer biomarkers, their robustness is limited by tumor heterogeneity. Recently, cancer-associated changes in RNA marks have emerged as a promising source of diagnostic and prognostic biomarkers. RNA epigenetics (also known as epitranscriptomics) is an emerging field in which at least 150 chemical modifications in all types of RNA (mRNA, tRNA, lncRNA, rRNA, and microRNA) have been detected. These modifications fine-tune gene expression in both physiological and pathological processes. A growing number of studies have established links between specific modified nucleoside levels in solid/liquid biopsies, and cancer onset and progression. In this review, we highlight the potential role of epitranscriptomic markers in refining cancer diagnosis and/or prognosis. RNA modification patterns may contain important information for establishing an initial diagnosis, monitoring disease evolution, and predicting response to treatment. Furthermore, recent developments in mass spectrometry allow reliable quantification of RNA marks in solid biopsies and biological fluids. We discuss the great potential of mass spectrometry for identifying epitranscriptomic biomarker signatures in cancer diagnosis. While there are various methods to quantify modified nucleosides, most are unable to detect and quantify more than one type of RNA modification at a time. Mass spectrometry analyses, especially GC-MS/MS and LC-MS/MS, overcome this limitation and simultaneously detect modified nucleosides by multiple reaction monitoring. Indeed, several groups are currently validating mass spectrometry methods that quantify several nucleosides at one time in liquid biopsies. The challenge now is to exploit these powerful analytical tools to establish epitranscriptomic signatures that should open new perspectives in personalized medicine. This review summarizes the growing clinical field of analysis of RNA modifications and discusses pre-analytical and analytical approaches, focusing in particular on the development of new mass spectrometry tools and their clinical applications.

Analysis of pseudouridines and other RNA modifications using HydraPsiSeq protocol

Detection of RNA modified nucleotides using deep sequencing can be performed by several approaches, including antibody-driven enrichment and natural or chemically induced RT signatures. However, only very few RNA modified nucleotides generate natural RT signatures and antibody-driven enrichment heavily depends on the quality of antibodies used and may be highly biased. Thus, the use of chemically-induced RT signatures is now considered as the most trusted experimental approach. In addition, the use of chemical reagents allows inclusion of simple "mock-treated" controls, to exclude spontaneous RT arrests, SNPs and other misincorporation-prone sites. Hydrazine is a well-known RNA-specific reagent, already extensively used in the past for RNA sequencing and structural probing. Hydrazine is highly reactive to U and shows low reaction rates with ψ residues, allowing their distinction by deep sequencing-based protocols. However, other modified RNA residues also show particular behavior upon hydrazine treatment. Here we present methodological developments allowing to use HydraPsiSeq for precise quantification of RNA pseudouridylation and also detection and quantification of some other RNA modifications, in addition to ψ residues.

Changes of the tRNA Modification Pattern during the Development of Dictyostelium discoideum

Dictyostelium discoideum is a social amoeba, which on starvation develops from a single-cell state to a multicellular fruiting body. This developmental process is accompanied by massive changes in gene expression, which also affect non-coding RNAs. Here, we investigate how tRNAs as key regulators of the translation process are affected by this transition. To this end, we used LOTTE-seq to sequence the tRNA pool of D. discoideum at different developmental time points and analyzed both tRNA composition and tRNA modification patterns. We developed a workflow for the specific detection of modifications from reverse transcriptase signatures in chemically untreated RNA-seq data at single-nucleotide resolution. It avoids the comparison of treated and untreated RNA-seq data using reverse transcription arrest patterns at nucleotides in the neighborhood of a putative modification site as internal control. We find that nucleotide modification sites in D. discoideum tRNAs largely conform to the modification patterns observed throughout the eukaroytes. However, there are also previously undescribed modification sites. We observe substantial dynamic changes of both expression levels and modification patterns of certain tRNA types during fruiting body development. Beyond the specific application to D. discoideum our results demonstrate that the developmental variability of tRNA expression and modification can be traced efficiently with LOTTE-seq.

Deciphering Epitranscriptome: Modification of mRNA Bases Provides a New Perspective for Post-transcriptional Regulation of Gene Expression

Gene regulation depends on dynamic and reversibly modifiable biological and chemical information in the epigenome/epitranscriptome. Accumulating evidence suggests that messenger RNAs (mRNAs) are generated in flashing bursts in the cells in a precisely regulated manner. However, the different aspects of the underlying mechanisms are not fully understood. Cellular RNAs are post-transcriptionally modified at the base level, which alters the metabolism of mRNA. The current understanding of epitranscriptome in the animal system is far ahead of that in plants. The accumulating evidence indicates that the epitranscriptomic changes play vital roles in developmental processes and stress responses. Besides being non-genetically encoded, they can be of reversible nature and involved in fine-tuning the expression of gene. However, different aspects of base modifications in mRNAs are far from adequate to assign the molecular basis/functions to the epitranscriptomic changes. Advances in the chemogenetic RNA-labeling and high-throughput next-generation sequencing techniques are enabling functional analysis of the epitranscriptomic modifications to reveal their roles in mRNA biology. Mapping of the common mRNA modifications, including N 6-methyladenosine (m6A), and 5-methylcytidine (m5C), have enabled the identification of other types of modifications, such as N 1-methyladenosine. Methylation of bases in a transcript dynamically regulates the processing, cellular export, translation, and stability of the mRNA; thereby influence the important biological and physiological processes. Here, we summarize the findings in the field of mRNA base modifications with special emphasis on m6A, m5C, and their roles in growth, development, and stress tolerance, which provide a new perspective for the regulation of gene expression through post-transcriptional modification. This review also addresses some of the scientific and technical issues in epitranscriptomic study, put forward the viewpoints to resolve the issues, and discusses the future perspectives of the research in this area.

A-to-I RNA Editing in Cancer: From Evaluating the Editing Level to Exploring the Editing Effects

As an important regulatory mechanism at the posttranscriptional level in metazoans, adenosine deaminase acting on RNA (ADAR)-induced A-to-I RNA editing modification of double-stranded RNA has been widely detected and reported. Editing may lead to non-synonymous amino acid mutations, RNA secondary structure alterations, pre-mRNA processing changes, and microRNA-mRNA redirection, thereby affecting multiple cellular processes and functions. In recent years, researchers have successfully developed several bioinformatics software tools and pipelines to identify RNA editing sites. However, there are still no widely accepted editing site standards due to the variety of parallel optimization and RNA high-seq protocols and programs. It is also challenging to identify RNA editing by normal protocols in tumor samples due to the high DNA mutation rate. Numerous RNA editing sites have been reported to be located in non-coding regions and can affect the biosynthesis of ncRNAs, including miRNAs and circular RNAs. Predicting the function of RNA editing sites located in non-coding regions and ncRNAs is significantly difficult. In this review, we aim to provide a better understanding of bioinformatics strategies for human cancer A-to-I RNA editing identification and briefly discuss recent advances in related areas, such as the oncogenic and tumor suppressive effects of RNA editing.

Epitranscriptome machinery in Trypanosomatids: New players on the table?

Trypanosoma and Leishmania parasites cause devastating tropical diseases resulting in serious global health consequences. These organisms have complex life cycles with mammalian hosts and insect vectors. The parasites must, therefore, survive in different environments, demanding rapid physiological and metabolic changes. These responses depend upon regulation of gene expression, which primarily occurs posttranscriptionally. Altering the composition or conformation of RNA through nucleotide modifications is one posttranscriptional mechanism of regulating RNA fate and function, and modifications including N6-methyladenosine (m6A), N1-methyladenosine (m1A), N5-methylcytidine (m5C), N4-acetylcytidine (ac4C), and pseudouridine (Ψ), dynamically regulate RNA stability and translation in diverse organisms. Little is known about RNA modifications and their machinery in Trypanosomatids, but we hypothesize that they regulate parasite gene expression and are vital for survival. Here, we identified Trypanosomatid homologs for writers of m1A, m5C, ac4C, and Ψ and analyze their evolutionary relationships. We systematically review the evidence for their functions and assess their potential use as therapeutic targets. This work provides new insights into the roles of these proteins in Trypanosomatid parasite biology and treatment of the diseases they cause and illustrates that Trypanosomatids provide an excellent model system to study RNA modifications, their molecular, cellular, and biological consequences, and their regulation and interplay.

Detection of genetic variation and base modifications at base-pair resolution on both DNA and RNA

Accurate decoding of nucleic acid variation is critical to understand the complexity and regulation of genome function. Here we use a single-molecule magnetic tweezer (MT) platform to identify sequence variation and map a range of important epigenetic base modifications with high sensitivity, specificity, and precision in the same single molecules of DNA or RNA. We have also developed a highly specific amplification-free CRISPR-Cas enrichment strategy to isolate genomic regions from native DNA. We demonstrate enrichment of DNA from both E. coli and the FMR1 5'UTR coming from cells derived from a Fragile X carrier. From these kilobase-length enriched molecules we could characterize the differential levels of adenine and cytosine base modifications on E. coli, and the repeat expansion length and methylation status of FMR1. Together these results demonstrate that our platform can detect a variety of genetic, epigenetic, and base modification changes concomitantly within the same single molecules.

Full-Range Profiling of tRNA Modifications Using LC-MS/MS at Single-Base Resolution through a Site-Specific Cleavage Strategy

Transfer RNAs (tRNAs) are the most heavily modified RNA species. Liquid chromatography coupled with mass spectrometry (LC-MS/MS) is a powerful tool for characterizing tRNA modifications, which involves pretreating tRNAs with base-specific ribonucleases to produce smaller oligonucleotides amenable to MS. However, the quality and quantity of products from base-specific digestions are severely impacted by the base composition of tRNAs. This often leads to a loss of sequence information. Here, we report a method for the full-range profiling of tRNA modifications at single-base resolution by combining site-specific RNase H digestion with the LC-MS/MS and RNA-seq techniques. The key steps were designed to generate high-quality products of optimal lengths and ionization properties. A linear correlation between collision energies and the m/z of oligonucleotides significantly improved the information content of collision-induced dissociation (CID) spectra. False positives were eliminated by up to 95% using novel inclusion criteria for collecting a census of modifications. This method is illustrated by the mapping of mouse mitochondrial tRNA^His(GUG) and tRNA^Val(UAC), which were hitherto not investigated. The identities and locations of the five species of modifications on these tRNAs were fully characterized. This approach is universally applicable to any tRNA species and provides an experimentally realizable pathway to the de novo sequencing of post-transcriptionally modified tRNAs with high sequence coverage.

Pseudouridine as a novel biomarker in prostate cancer

Epitranscriptomic analysis has recently led to the profiling of modified nucleosides in cancer cell biological matrices, helping to elucidate their functional roles in cancer and reigniting interest in exploring their use as potential markers of cancer development and progression. Pseudouridine, one of the most well-known and the most abundant of the RNA nucleotide modifications, is the C5-glycoside isomer of uridine and its distinctive physiochemical properties allows it to perform many essential functions. Pseudouridine functionally (a) confers rigidity to local RNA structure by enhancing RNA stacking, engaging in a cooperative effect on neighboring nucleosides that overall contributes to RNA stabilization (b) refines the structure of tRNAs, which influences their decoding activity (c) facilitates the accuracy of decoding and proofreading during translation and efficiency of peptide bond formation, thus collectively improving the fidelity of protein biosynthesis and (e) dynamically regulates mRNA coding and translation. Biochemical synthesis of pseudouridine is carried out by pseudouridine synthases. In this review we discuss the evidence supporting an association between elevated pseudouridine levels with the incidence and progression of human prostate cancer and the translational significance of the value of this modified nucleotide as a novel biomarker in prostate cancer progression to advanced disease.

Combining tRNA sequencing methods to characterize plant tRNA expression and post-transcriptional modification

Differences in tRNA expression have been implicated in a remarkable number of biological processes. There is growing evidence that tRNA genes can play dramatically different roles depending on both expression and post-transcriptional modification, yet sequencing tRNAs to measure abundance and detect modifications remains challenging. Their secondary structure and extensive post-transcriptional modifications interfere with RNA-seq library preparation methods and have limited the utility of high-throughput sequencing technologies. Here, we combine two modifications to standard RNA-seq methods by treating with the demethylating enzyme AlkB and ligating with tRNA-specific adapters in order to sequence tRNAs from four species of flowering plants, a group that has been shown to have some of the most extensive rates of post-transcriptional tRNA modifications. This protocol has the advantage of detecting full-length tRNAs and sequence variants that can be used to infer many post-transcriptional modifications. We used the resulting data to produce a modification index of almost all unique reference tRNAs in Arabidopsis thaliana, which exhibited many anciently conserved similarities with humans but also positions that appear to be 'hot spots' for modifications in angiosperm tRNAs. We also found evidence based on northern blot analysis and droplet digital PCR that, even after demethylation treatment, tRNA-seq can produce highly biased estimates of absolute expression levels most likely due to biased reverse transcription. Nevertheless, the generation of full-length tRNA sequences with modification data is still promising for assessing differences in relative tRNA expression across treatments, tissues or subcellular fractions and help elucidate the functional roles of tRNA modifications.

AlkAniline-Seq: A Highly Sensitive and Specific Method for Simultaneous Mapping of 7-Methyl-guanosine (m7G) and 3-Methyl-cytosine (m3C) in RNAs by High-Throughput Sequencing

Epitranscriptomics is an emerging field where the development of high-throughput analytical technologies is essential to profile the dynamics of RNA modifications under different conditions. Despite important advances during the last 10 years, the number of RNA modifications detectable by next-generation sequencing is restricted to a very limited subset. Here, we describe a highly efficient and fast method called AlkAniline-Seq to map simultaneously two different RNA modifications: 7-methyl-guanosine (m⁷G) and 3-methyl-cytosine (m³C) in RNA. Our protocol is based on three subsequent chemical/enzymatic steps allowing the enrichment of RNA fragments ending at position n + 1 to the modified nucleotide, without any prior RNA selection. Therefore, AlkAniline-Seq demonstrates an outstanding sensitivity and specificity for these two RNA modifications. We have validated AlkAniline-Seq using bacterial, yeast, and human total RNA, and here we present, as an example, a synthetic view of the complete profiling of these RNA modifications in S. cerevisiae tRNAs.

Single enzyme RT-PCR of full-length ribosomal RNA

The ribosome is a two-subunit, macromolecular machine composed of RNA and proteins that carries out the polymerization of α-amino acids into polypeptides. Efforts to engineer ribosomal RNA (rRNA) deepen our understanding of molecular translation and provide opportunities to expand the chemistry of life by creating ribosomes with altered properties. Toward these efforts, reverse transcription PCR (RT-PCR) of the entire 16S and 23S rRNAs, which make up the 30S small subunit and 50S large subunit, respectively, is important for isolating desired phenotypes. However, reverse transcription of rRNA is challenging due to extensive secondary structure and post-transcriptional modifications. One key challenge is that existing commercial kits for RT-PCR rely on reverse transcriptases that lack the extreme thermostability and processivity found in many commercial DNA polymerases, which can result in subpar performance on challenging templates. Here, we develop methods employing a synthetic thermostable reverse transcriptase (RTX) to enable and optimize RT-PCR of the complete Escherichia coli 16S and 23S rRNAs. We also characterize the error rate of RTX when traversing the various post-transcriptional modifications of the 23S rRNA. We anticipate that this work will facilitate efforts to study and characterize many naturally occurring long RNAs and to engineer the translation apparatus for synthetic biology.

HydraPsiSeq: a method for systematic and quantitative mapping of pseudouridines in RNA

Developing methods for accurate detection of RNA modifications remains a major challenge in epitranscriptomics. Next-generation sequencing-based mapping approaches have recently emerged but, often, they are not quantitative and lack specificity. Pseudouridine (ψ), produced by uridine isomerization, is one of the most abundant RNA modification. ψ mapping classically involves derivatization with soluble carbodiimide (CMCT), which is prone to variation making this approach only semi-quantitative. Here, we developed 'HydraPsiSeq', a novel quantitative ψ mapping technique relying on specific protection from hydrazine/aniline cleavage. HydraPsiSeq is quantitative because the obtained signal directly reflects pseudouridine level. Furthermore, normalization to natural unmodified RNA and/or to synthetic in vitro transcripts allows absolute measurements of modification levels. HydraPsiSeq requires minute amounts of RNA (as low as 10-50 ng), making it compatible with high-throughput profiling of diverse biological and clinical samples. Exploring the potential of HydraPsiSeq, we profiled human rRNAs, revealing strong variations in pseudouridylation levels at ∼20-25 positions out of total 104 sites. We also observed the dynamics of rRNA pseudouridylation throughout chondrogenic differentiation of human bone marrow stem cells. In conclusion, HydraPsiSeq is a robust approach for the systematic mapping and accurate quantification of pseudouridines in RNAs with applications in disease, aging, development, differentiation and/or stress response.

Functional characterization of the human tRNA methyltransferases TRMT10A and TRMT10B

The TRM10 family of methyltransferases is responsible for the N¹-methylation of purines at position 9 of tRNAs in Archaea and Eukarya. The human genome encodes three TRM10-type enzymes, of which only the mitochondrial TRMT10C was previously characterized in detail, whereas the functional significance of the two presumably nuclear enzymes TRMT10A and TRMT10B remained unexplained. Here we show that TRMT10A is m¹G9-specific and methylates a subset of nuclear-encoded tRNAs, whilst TRMT10B is the first m¹A9-specific tRNA methyltransferase found in eukaryotes and is responsible for the modification of a single nuclear-encoded tRNA. Furthermore, we show that the lack of G9 methylation causes a decrease in the steady-state levels of the initiator tRNA^iMet-CAT and an alteration in its further post-transcriptional modification. Our work finally clarifies the function of TRMT10A and TRMT10B in vivo and provides evidence that the loss of TRMT10A affects the pool of cytosolic tRNAs required for protein synthesis.

Male reproductive health and intergenerational metabolic responses from a small RNA perspective

The world has recently experienced a decline in male reproductive (e.g. sperm counts and motility) and metabolic (e.g. obesity and diabetes) health. Accumulated evidence from animal models also shows that the metabolic health of the father may influence the metabolic health in his offspring. Vectors for such paternal intergenerational metabolic responses (IGMRs) involve small noncoding RNAs (sncRNAs) that often increase in spermatozoa during the last days of maturation in the epididymis. We and others have shown that the metabolic state - depending on factors such as diet, obesity and physical exercise - may affect sperm quality and sperm sncRNA. Together, this suggests that there are overlapping aetiologies between the male metabolic syndrome, male factor infertility and intergenerational responses. In this review, we present a theoretical framework for an overlap of these aetiologies by exploring the advances in our understanding of the roles of sncRNA in spermatogenesis and offspring development. A special focus will lie on novel findings about tRNA-derived small RNA (tsRNA), rRNA-derived small RNA (rsRNA) and small mitochondrial RNA (mitoRNA), and their emerging roles in intergenerational metabolic and reproductive health.

NAD tagSeq for transcriptome-wide identification and characterization of NAD+-capped RNAs

Several noncanonical initial nucleotides (NCINs) have been found to cap RNAs and possibly regulate RNA stability, transcription and translation. NAD⁺ is one of the NCINs that has recently been discovered to cap RNAs in a wide range of species. Identification of the NAD⁺-capped RNAs (NAD-RNAs) could help to unveil the cap-mediated regulation mechanisms. We previously reported a method termed NAD tagSeq for genome-wide analysis of NAD-RNAs. NAD tagSeq is based on the previously published NAD captureSeq protocol, which applies an enzymatic reaction and a click chemistry reaction to label NAD-RNAs with biotin followed by enrichment with streptavidin resin and identification by RNA sequencing. In the current NAD tagSeq method, NAD-RNAs are labeled with a synthetic RNA tag and identified by direct RNA sequencing based on Oxford Nanopore technology. Compared to NAD captureSeq, NAD tagSeq provides a simpler procedure for direct sequencing of NAD-RNAs and noncapped RNAs and affords information on the whole sequence organization of NAD-RNAs and the ratio of NAD-RNAs to total transcripts. Furthermore, NAD-RNAs can be enriched by hybridizing a complementary DNA probe to the RNA tag, thus increasing the sequencing coverage of NAD-RNAs. The strategy of tagging RNAs with a synthetic RNA tag and identifying them by direct RNA sequencing might be employed in analyzing other NCIN-capped RNAs. The experimental procedure of NAD tagSeq, including RNA extraction, RNA tagging and direct RNA sequencing, takes ~5 d, and initial data analysis can be completed within 2 d.

Accurate characterization of Escherichia coli tRNA modifications with a simple method of deep-sequencing library preparation

In conventional RNA high-throughput sequencing, modified bases prevent a large fraction of tRNA transcripts to be converted into cDNA libraries. Recent proposals aiming at resolving this issue take advantage of the interference of base modifications with RT enzymes to detect and identify them by establishing signals from aborted cDNA transcripts. Because some modifications, such as methyl groups, do almost not allow RT bypassing, demethylation and highly processive RT enzymes have been used to overcome these obstacles. Working with Escherichia coli as a model system, we show that with a conventional (albeit still engineered) RT enzyme and key optimizations in library preparation, all RT-impairing modifications can be highlighted along the entire tRNA length without demethylation procedure. This is achieved by combining deep-sequencing samples, which allows to establish aborted transcription signal of higher accuracy and reproducibility, with the potential for differentiating tiny differences in the state of modification of all cellular tRNAs. In addition, our protocol provides estimates of the relative tRNA abundance.

A Census and Categorization Method of Epitranscriptomic Marks

In the past few years, thorough investigation of chemical modifications operated in the cells on ribonucleic acid (RNA) molecules is gaining momentum. This new field of research has been dubbed “epitranscriptomics”, in analogy to best-known epigenomics, to stress the potential of ensembles of RNA modifications to constitute a post-transcriptional regulatory layer of gene expression orchestrated by writer, reader, and eraser RNA-binding proteins (RBPs). In fact, epitranscriptomics aims at identifying and characterizing all functionally relevant changes involving both non-substitutional chemical modifications and editing events made to the transcriptome. Indeed, several types of RNA modifications that impact gene expression have been reported so far in different species of cellular RNAs, including ribosomal RNAs, transfer RNAs, small nuclear RNAs, messenger RNAs, and long non-coding RNAs. Supporting functional relevance of this largely unknown regulatory mechanism, several human diseases have been associated directly to RNA modifications or to RBPs that may play as effectors of epitranscriptomic marks. However, an exhaustive epitranscriptome’s characterization, aimed to systematically classify all RNA modifications and clarify rules, actors, and outcomes of this promising regulatory code, is currently not available, mainly hampered by lack of suitable detecting technologies. This is an unfortunate limitation that, thanks to an unprecedented pace of technological advancements especially in the sequencing technology field, is likely to be overcome soon. Here, we review the current knowledge on epitranscriptomic marks and propose a categorization method based on the reference ribonucleotide and its rounds of modifications (“stages”) until reaching the given modified form. We believe that this classification scheme can be useful to coherently organize the expanding number of discovered RNA modifications.

Mitochondrial tRNA methylation in Alzheimer's disease and progressive supranuclear palsy

BACKGROUND

Methylation of mitochondrial tRNAs (mt-tRNA) at the 9th position ("p9 site") is known to impact translational efficiency and downstream mitochondrial function; however, direct assessment of mt-RNA methylation is challenging. Recent RNA sequence-based methods have been developed to reliably identify post-transcriptional methylation. Though p9 methylation has been studied in healthy human populations and in the context of cancer, it has not yet been analyzed in neurodegenerative disease, where mitochondrial dysfunction is a prominent and early hallmark of disease progression.

METHODS

Mitochondrial p9 methylation was inferred from multi-allelic calls in RNA-seq data. Gene-based association studies were performed in FUMA. Correlations between nuclear gene expression and p9 methylation were tested using Spearman's rho. Fisher's Exact test was used in PANTHER and IPA to test for overrepresentation and enrichment of biological processes and pathways in the top nuclear genes correlated with p9 methylation.

RESULTS

Variable methylation was observed at 11 p9 sites in post-mortem cerebellar tissue of elderly subjects who were either healthy or diagnosed with Alzheimer's disease (AD), progressive supranuclear palsy (PSP) or pathological aging (PA). Similarities in degree of methylation were observed between AD and PSP. Certain nuclear encoded genes were identified as significantly associated with p9 methylation. Expression of 5300 nuclear encoded genes was significantly correlated with p9 methylation, with AD and PSP subjects exhibiting similar expression profiles. Overrepresentation and enrichment testing using the top transcripts revealed enrichment for a number of molecular processes, terms and pathways including many of which that were mitochondrial-related.

CONCLUSION

With mitochondrial dysfunction being an established hallmark of neurodegenerative disease pathophysiology, this work sheds light on the potential molecular underpinnings of this dysfunction. Here we show overlap in cerebellar pathophysiology between common tauopathies such as Alzheimer's disease and progressive supranuclear palsy. Whether p9 hypermethylation is a cause or consequence of pathology remains an area of focus.

snoRNPs: Functions in Ribosome Biogenesis

Ribosomes are perhaps the most critical macromolecular machine as they are tasked with carrying out protein synthesis in cells. They are incredibly complex structures composed of protein components and heavily chemically modified RNAs. The task of assembling mature ribosomes from their component parts consumes a massive amount of energy and requires greater than 200 assembly factors. Among the most critical of these are small nucleolar ribonucleoproteins (snoRNPs). These are small RNAs complexed with diverse sets of proteins. As suggested by their name, they localize to the nucleolus, the site of ribosome biogenesis. There, they facilitate multiple roles in ribosomes biogenesis, such as pseudouridylation and 2′-O-methylation of ribosomal (r)RNA, guiding pre-rRNA processing, and acting as molecular chaperones. Here, we reviewed their activity in promoting the assembly of ribosomes in eukaryotes with regards to chemical modification and pre-rRNA processing.