Profiling of ribose methylations in RNA by high-throughput sequencing

txtools: an R package facilitating analysis of RNA modifications, structures, and interactions

We present txtools, an R package that enables the processing, analysis, and visualization of RNA-seq data at the nucleotide-level resolution, seamlessly integrating alignments to the genome with transcriptomic representation. txtools' main inputs are BAM files and a transcriptome annotation, and the main output is a table, capturing mismatches, deletions, and the number of reads beginning and ending at each nucleotide in the transcriptomic space. txtools further facilitates downstream visualization and analyses. We showcase, using examples from the epitranscriptomic field, how a few calls to txtools functions can yield insightful and ready-to-publish results. txtools is of broad utility also in the context of structural mapping and RNA:protein interaction mapping. By providing a simple and intuitive framework, we believe that txtools will be a useful and convenient tool and pave the path for future discovery. txtools is available for installation from its GitHub repository at https://github.com/AngelCampos/txtools.

2'-O-methylation (Nm) in RNA: progress, challenges, and future directions

RNA 2'-O-methylation (Nm) is highly abundant in noncoding RNAs including ribosomal RNA (rRNA), transfer RNA (tRNA), and small nuclear RNA (snRNA), and occurs in the 5' cap of virtually all messenger RNAs (mRNAs) in higher eukaryotes. More recently, Nm has also been reported to occur at internal sites in mRNA. High-throughput methods have been developed for the transcriptome-wide detection of Nm. However, these methods have mostly been applied to abundant RNAs such as rRNA, and the validity of the internal mRNA Nm sites detected with these approaches remains controversial. Nonetheless, Nm in both coding and noncoding RNAs has been demonstrated to impact cellular processes, including translation and splicing. In addition, Nm modifications at the 5' cap and possibly at internal sites in mRNA serve to prevent the binding of nucleic acid sensors, thus preventing the activation of the innate immune response by self-mRNAs. Finally, Nm has been implicated in a variety of diseases including cancer, cardiovascular diseases, and neurologic syndromes. In this review, we discuss current challenges in determining the distribution, regulation, function, and disease relevance of Nm, as well as potential future directions for the field.

Comprehensive map of ribosomal 2'-O-methylation and C/D box snoRNAs in Drosophila melanogaster

During their maturation, ribosomal RNAs (rRNAs) are decorated by hundreds of chemical modifications that participate in proper folding of rRNA secondary structures and therefore in ribosomal function. Along with pseudouridine, methylation of the 2'-hydroxyl ribose moiety (Nm) is the most abundant modification of rRNAs. The majority of Nm modifications in eukaryotes are placed by Fibrillarin, a conserved methyltransferase belonging to a ribonucleoprotein complex guided by C/D box small nucleolar RNAs (C/D box snoRNAs). These modifications impact interactions between rRNAs, tRNAs and mRNAs, and some are known to fine tune translation rates and efficiency. In this study, we built the first comprehensive map of Nm sites in Drosophila melanogaster rRNAs using two complementary approaches (RiboMethSeq and Nanopore direct RNA sequencing) and identified their corresponding C/D box snoRNAs by whole-transcriptome sequencing. We de novo identified 61 Nm sites, from which 55 are supported by both sequencing methods, we validated the expression of 106 C/D box snoRNAs and we predicted new or alternative rRNA Nm targets for 31 of them. Comparison of methylation level upon different stresses show only slight but specific variations, indicating that this modification is relatively stable in D. melanogaster. This study paves the way to investigate the impact of snoRNA-mediated 2'-O-methylation on translation and proteostasis in a whole organism.

Complicated target recognition by archaeal box C/D guide RNAs

Box C/D RNAs guide the site-specific formation of 2'-O-methylated nucleotides (Nm) of RNAs in eukaryotes and archaea. Although C/D RNAs have been profiled in several archaea, their targets have not been experimentally determined. Here, we mapped Nm in rRNAs, tRNAs, and abundant small RNAs (sRNAs) and profiled C/D RNAs in the crenarchaeon Sulfolobus islandicus. The targets of C/D RNAs were assigned by analysis of base-pairing interactions, in vitro modification assays, and gene deletion experiments, revealing a complicated landscape of C/D RNA-target interactions. C/D RNAs widely use dual antisense elements to target adjacent sites in rRNAs, enhancing modification at weakly bound sites. Two consecutive sites can be guided with the same antisense element upstream of box D or D', a phenomenon known as double-specificity that is exclusive to internal box D' in eukaryotic C/D RNAs. Several C/D RNAs guide modification at a single non-canonical site. This study reveals the global landscape of RNA-guided 2'-O-methylation in an archaeon and unexpected targeting rules employed by C/D RNA.

An integrative platform for detection of RNA 2'-O-methylation reveals its broad distribution on mRNA

Ribose 2'-O-methylation is involved in critical biological processes, but its biological functions and significance in mRNAs remain underexplored. We have developed NJU-seq, a sensitive method for unbiased 2'-O-methylation (Nm) profiling, and Nm-VAQ, a site-specific quantification tool. Using these tools in tandem, we identified thousands of Nm sites on mRNAs of human and mouse cell lines, of which 68 of 84 selected sites were further validated to be more than 1% 2'-O-methylated. Unlike rRNA, most mRNA Nm sites were from 1% to 30% methylated. In addition, mRNA Nm was dynamic, changing according to the circumstance. Furthermore, we show that fibrillarin is involved as a methyltransferase. By mimicking the detected Nm sites and the context sequence, the RNA fragments could be 2'-O-methylated and demonstrated higher stability but lower translation efficiency. Last, profiling of Nm sites in lung surgery samples revealed common signatures of lung cancer pathogenesis, providing potential new diagnostic markers.

GPATCH4 regulates rRNA and snRNA 2'-O-methylation in both DHX15-dependent and DHX15-independent manners

Regulation of RNA helicase activity, often accomplished by protein cofactors, is essential to ensure target specificity within the complex cellular environment. The largest family of RNA helicase cofactors are the G-patch proteins, but the cognate RNA helicases and cellular functions of numerous human G-patch proteins remain elusive. Here, we discover that GPATCH4 is a stimulatory cofactor of DHX15 that interacts with the DEAH box helicase in the nucleolus via residues in its G-patch domain. We reveal that GPATCH4 associates with pre-ribosomal particles, and crosslinks to the transcribed ribosomal DNA locus and precursor ribosomal RNAs as well as binding to small nucleolar- and small Cajal body-associated RNAs that guide rRNA and snRNA modifications. Loss of GPATCH4 impairs 2'-O-methylation at various rRNA and snRNA sites leading to decreased protein synthesis and cell growth. We demonstrate that the regulation of 2'-O-methylation by GPATCH4 is both dependent on, and independent of, its interaction with DHX15. Intriguingly, the ATPase activity of DHX15 is necessary for efficient methylation of DHX15-dependent sites, suggesting a function of DHX15 in regulating snoRNA-guided 2'-O-methylation of rRNA that requires activation by GPATCH4. Overall, our findings extend knowledge on RNA helicase regulation by G-patch proteins and also provide important new insights into the mechanisms regulating installation of rRNA and snRNA modifications, which are essential for ribosome function and pre-mRNA splicing.

2'-O-ribose methylation levels of ribosomal RNA distinguish different types of growth arrest in human dermal fibroblasts

The 2'-O-methylation (2'-O-Me) of ribosomal RNA (rRNA) shows plasticity that is potentially associated with cell phenotypes. We used RiboMeth-seq profiling to reveal growth arrest-specific 2'-O-Me patterns in primary human dermal fibroblasts from three different donors. We exposed cells to hydrogen peroxide to induce cellular senescence and to high cell densities to promote quiescence by contact inhibition. We compared both modes of cell cycle arrest to proliferating cells and could indeed distinguish these conditions by their overall 2'-O-Me patterns. Methylation levels at a small fraction of sites showed plasticity and correlated with the expression of specific small nucleolar RNAs (snoRNAs) but not with expression of fibrillarin. Moreover, we observed subtle senescence-associated alterations in ribosome biogenesis. Knockdown of the snoRNA SNORD87, which acts as a guide for modification of a hypermethylated position in non-proliferating cells, was sufficient to boost cell proliferation. Conversely, depletion of SNORD88A, SNORD88B and SNORD88C, which act as guides for modification of a hypomethylated site, caused decreased proliferation without affecting global protein synthesis or apoptosis. Taken together, our findings provide evidence that rRNA modifications can be used to distinguish and potentially influence specific growth phenotypes of primary cells.

Decoding Ribosome Heterogeneity: A New Horizon in Cancer Therapy

The traditional perception of ribosomes as uniform molecular machines has been revolutionized by recent discoveries, revealing a complex landscape of ribosomal heterogeneity. Opposing the conventional belief in interchangeable ribosomal entities, emerging studies underscore the existence of specialized ribosomes, each possessing unique compositions and functions. Factors such as cellular and tissue specificity, developmental and physiological states, and external stimuli, including circadian rhythms, significantly influence ribosome compositions. For instance, muscle cells and neurons are characterized by distinct ribosomal protein sets and dynamic behaviors, respectively. Furthermore, alternative forms of ribosomal RNA (rRNAs) and their post-transcriptional modifications add another dimension to this heterogeneity. These variations, orchestrated by spatial, temporal, and conditional factors, enable the manifestation of a broad spectrum of specialized ribosomes, each tailored for potentially distinct functions. Such specialization not only impacts mRNA translation and gene expression but also holds significant implications for broader biological contexts, notably in the realm of cancer research. As the understanding of ribosomal diversity deepens, it also paves the way for exploring novel avenues in cellular function and offers a fresh perspective on the molecular intricacies of translation.

Base-Resolution Sequencing Methods for Whole-Transcriptome Quantification of mRNA Modifications

ConspectusRNA molecules are not merely a combination of four bases of A, C, G, and U. Chemical modifications occur in almost all RNA species and play diverse roles in gene expression regulation. The abundant cellular RNAs, such as ribosomal RNA (rRNA) and transfer RNA (tRNA), are known to have the highest density of RNA modifications, which exert critical functions in rRNA and tRNA biogenesis, stability, and subsequent translation. In recent years, modifications on low-abundance RNA species in mammalian cells, such as messenger RNA (mRNA), regulatory noncoding RNA (ncRNA), and chromatin-associated RNA (caRNA), have been shown to contain multiple different chemical modifications with functional significance.As the most abundant mRNA modification in mammals, N6-methyladenosine (m6A) affects nearly every stage of mRNA processing and metabolism, with the antibody-based m6A-MeRIP-seq (methylated RNA immunoprecipitation sequencing) followed by high-throughput sequencing widely employed in mapping m6A distribution transcriptome-wide in diverse biological systems. In addition to m6A, other chemical modifications such as pseudouridine (Ψ), 2'-O-methylation (Nm), 5-methylcytidine (m5C), internal N7-methylguanosine (m7G), N1-methyladenosine (m1A), N4-acetylcytidine (ac4C), etc. also exist in polyA-tailed RNA in mammalian cells, requiring effective mapping approaches for whole-transcriptome profiling of these non-m6A mRNA modifications. Like m6A, the antibody-based enrichment followed by sequencing has been the primary method to study distributions of these modifications. Methods to more quantitatively map these modifications would dramatically improve our understanding of distributions and modification density of these chemical marks on RNA, thereby bettering informing functional implications. In this Account, aimed at both single-base resolution and modification fraction quantification, we summarize our recent advances in developing a series of chemistry- or biochemistry-based methods to quantitatively map RNA modifications, including m6A, Ψ, m5C, m1A, 2'-O-methylation (Nm), and internal m7G, in mammalian mRNA at base resolution. These new methods, including m6A-SAC-seq, eTAM-seq, BID-seq, UBS-seq, DAMM-seq, m1A-quant-seq, Nm-Mut-seq, and m7G-quant-seq, promise to conduct base-resolution mapping of most major mRNA modifications with low RNA input and uncover dynamic changes in modification stoichiometry during biological and physiological processes, facilitating future investigations on these RNA modifications in regulating cellular gene expression and as potential biomarkers for clinical diagnosis and prognosis. These quantitative sequencing methods allow the mapping of most mRNA modifications with limited input sample requirements. The same modifications on diverse RNA species, such as caRNA, ncRNA, nuclear nascent RNA, mitochondrial RNA, cell-free RNA (cfRNA), etc., could be sequenced using the same methods.

Single-cell analysis of the epitranscriptome: RNA modifications under the microscope

The identification of mechanisms capable of modifying genetic information by the addition of covalent RNA modifications distinguishes a level of complexity in gene expression which challenges key long-standing concepts of RNA biology. One of the current challenges of molecular biology is to properly understand the molecular functions of these RNA modifications, with more than 170 different ones having been identified so far. However, it has not been possible to map specific RNA modifications at a single-cell resolution until very recently. This review will highlight the technological advances in single-cell methodologies aimed at assessing and testing the biological function of certain RNA modifications, focusing on m6A. These advances have allowed for the development of novel strategies that enable the study of the 'epitranscriptome'. Nevertheless, despite all these improvements, many challenges and difficulties still need fixing for these techniques to work efficiently.

Comparative RNA Genomics

Over the last quarter of a century it has become clear that RNA is much more than just a boring intermediate in protein expression. Ancient RNAs still appear in the core information metabolism and comprise a surprisingly large component in bacterial gene regulation. A common theme with these types of mostly small RNAs is their reliance of conserved secondary structures. Large-scale sequencing projects, on the other hand, have profoundly changed our understanding of eukaryotic genomes. Pervasively transcribed, they give rise to a plethora of large and evolutionarily extremely flexible non-coding RNAs that exert a vastly diverse array of molecule functions. In this chapter we provide a-necessarily incomplete-overview of the current state of comparative analysis of non-coding RNAs, emphasizing computational approaches as a means to gain a global picture of the modern RNA world.

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.

Isocitrate dehydrogenase wt and IDHmut adult-type diffuse gliomas display distinct alterations in ribosome biogenesis and 2'O-methylation of ribosomal RNA

BACKGROUND

High-grade adult-type diffuse gliomas (HGGs) constitute a heterogeneous group of aggressive tumors that are mostly incurable. Recent advances highlighting the contribution of ribosomes to cancer development have offered new clinical perspectives. Here, we uncovered that isocitrate dehydrogenase (IDH)wt and IDHmut HGGs display distinct alterations of ribosome biology, in terms of rRNA epitranscriptomics and ribosome biogenesis, which could constitute novel hallmarks that can be exploited for the management of these pathologies.

METHODS

We analyzed (1) the ribosomal RNA 2'O-ribose methylation (rRNA 2'Ome) using RiboMethSeq and in-house developed bioinformatics tools (https://github.com/RibosomeCRCL/ribomethseq-nfandrRMSAnalyzer) on 3 independent cohorts compiling 71 HGGs (IDHwt n = 30, IDHmut n = 41) and 9 non-neoplastic samples, (2) the expression of ribosome biogenesis factors using medium throughput RT-qPCR as a readout of ribosome biogenesis, and (3) the sensitivity of 5 HGG cell lines to RNA Pol I inhibitors (CX5461, BMH-21).

RESULTS

Unsupervised analysis demonstrated that HGGs could be distinguished based on their rRNA 2'Ome epitranscriptomic profile, with IDHwt glioblastomas displaying the most significant alterations of rRNA 2'Ome at specific sites. In contrast, IDHmut HGGs are largely characterized by an overexpression of ribosome biogenesis factors compared to non-neoplastic tissues or IDHwt glioblastomas. Finally, IDHmut HGG-derived spheroids display higher cytotoxicity to CX5461 than IDHwt glioblastoma, while all HGG spheroids display a similar cytotoxicity to BMH-21.

CONCLUSIONS

In HGGs, IDH mutational status is associated with specific alterations of the ribosome biology and with distinct sensitivities to RNA Pol I inhibitors.

LncRNA INHEG promotes glioma stem cell maintenance and tumorigenicity through regulating rRNA 2'-O-methylation

Glioblastoma (GBM) ranks among the most lethal of human cancers, containing glioma stem cells (GSCs) that display therapeutic resistance. Here, we report that the lncRNA INHEG is highly expressed in GSCs compared to differentiated glioma cells (DGCs) and promotes GSC self-renewal and tumorigenicity through control of rRNA 2'-O-methylation. INHEG induces the interaction between SUMO2 E3 ligase TAF15 and NOP58, a core component of snoRNP that guides rRNA methylation, to regulate NOP58 sumoylation and accelerate the C/D box snoRNP assembly. INHEG activation enhances rRNA 2'-O-methylation, thereby increasing the expression of oncogenic proteins including EGFR, IGF1R, CDK6 and PDGFRB in glioma cells. Taken together, this study identifies a lncRNA that connects snoRNP-guided rRNA 2'-O-methylation to upregulated protein translation in GSCs, supporting an axis for potential therapeutic targeting of gliomas.

Disease Diagnosis Based on Nucleic Acid Modifications

Nucleic acid modifications include a wide range of epigenetic and epitranscriptomic factors and impact a wide range of nucleic acids due to their profound influence on biological inheritance, growth, and metabolism. The recently developed methods of mapping and characterizing these modifications have promoted their discovery as well as large-scale studies in eukaryotes, especially in humans. Because of these pioneering strategies, nucleic acid modifications have been shown to have a great impact on human disorders such as cancer. Therefore, whether nucleic acid modifications could become a new type of biomarker remains an open question. In this review, we briefly look back at classical nucleic acid modifications and then focus on the progress made in investigating these modifications as diagnostic biomarkers in clinical therapy and present our perspective on their development prospects.

A snoRNA-tRNA modification network governs codon-biased cellular states

The dynamic balance between tRNA supply and codon usage demand is a fundamental principle in the cellular translation economy. However, the regulation and functional consequences of this balance remain unclear. Here, we use PARIS2 interactome capture, structure modeling, conservation analysis, RNA-protein interaction analysis, and modification mapping to reveal the targets of hundreds of snoRNAs, many of which were previously considered orphans. We identify a snoRNA-tRNA interaction network that is required for global tRNA modifications, including 2'-O-methylation and others. Loss of Fibrillarin, the snoRNA-guided 2'-O-methyltransferase, induces global upregulation of tRNA fragments, a large group of regulatory RNAs. In particular, the snoRNAs D97/D133 guide the 2'-O-methylation of multiple tRNAs, especially for the amino acid methionine (Met), a protein-intrinsic antioxidant. Loss of D97/D133 snoRNAs in human HEK293 cells reduced target tRNA levels and induced codon adaptation of the transcriptome and translatome. Both single and double knockouts of D97 and D133 in HEK293 cells suppress Met-enriched proliferation-related gene expression programs, including, translation, splicing, and mitochondrial energy metabolism, and promote Met-depleted programs related to development, differentiation, and morphogenesis. In a mouse embryonic stem cell model of development, knockdown and knockout of D97/D133 promote differentiation to mesoderm and endoderm fates, such as cardiomyocytes, without compromising pluripotency, consistent with the enhanced development-related gene expression programs in human cells. This work solves a decades-old mystery about orphan snoRNAs and reveals a function of snoRNAs in controlling the codon-biased dichotomous cellular states of proliferation and development.

RNA modification: mechanisms and therapeutic targets

RNA modifications are dynamic and reversible chemical modifications on substrate RNA that are regulated by specific modifying enzymes. They play important roles in the regulation of many biological processes in various diseases, such as the development of cancer and other diseases. With the help of advanced sequencing technologies, the role of RNA modifications has caught increasing attention in human diseases in scientific research. In this review, we briefly summarized the basic mechanisms of several common RNA modifications, including m6A, m5C, m1A, m7G, Ψ, A-to-I editing and ac4C. Importantly, we discussed their potential functions in human diseases, including cancer, neurological disorders, cardiovascular diseases, metabolic diseases, genetic and developmental diseases, as well as immune disorders. Through the "writing-erasing-reading" mechanisms, RNA modifications regulate the stability, translation, and localization of pivotal disease-related mRNAs to manipulate disease development. Moreover, we also highlighted in this review all currently available RNA-modifier-targeting small molecular inhibitors or activators, most of which are designed against m6A-related enzymes, such as METTL3, FTO and ALKBH5. This review provides clues for potential clinical therapy as well as future study directions in the RNA modification field. More in-depth studies on RNA modifications, their roles in human diseases and further development of their inhibitors or activators are needed for a thorough understanding of epitranscriptomics as well as diagnosis, treatment, and prognosis of human diseases.

The ribose methylation enzyme FTSJ1 has a conserved role in neuron morphology and learning performance

FTSJ1 is a conserved human 2'-O-methyltransferase (Nm-MTase) that modifies several tRNAs at position 32 and the wobble position 34 in the anticodon loop. Its loss of function has been linked to X-linked intellectual disability (XLID), and more recently to cancers. However, the molecular mechanisms underlying these pathologies are currently unclear. Here, we report a novel FTSJ1 pathogenic variant from an X-linked intellectual disability patient. Using blood cells derived from this patient and other affected individuals carrying FTSJ1 mutations, we performed an unbiased and comprehensive RiboMethSeq analysis to map the ribose methylation on all human tRNAs and identify novel targets. In addition, we performed a transcriptome analysis in these cells and found that several genes previously associated with intellectual disability and cancers were deregulated. We also found changes in the miRNA population that suggest potential cross-regulation of some miRNAs with these key mRNA targets. Finally, we show that differentiation of FTSJ1-depleted human neural progenitor cells into neurons displays long and thin spine neurites compared with control cells. These defects are also observed in Drosophila and are associated with long-term memory deficits. Altogether, our study adds insight into FTSJ1 pathologies in humans and flies by the identification of novel FTSJ1 targets and the defect in neuron morphology.

Reverse Transcriptases: From Discovery and Applications to Xenobiology

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

Comprehensive identification of diverse ribosomal RNA modifications by targeted nanopore direct RNA sequencing and JACUSA2

Ribosomal RNAs are decorated by numerous post-transcriptional modifications whose exact roles in ribosome biogenesis, function, and human pathophysiology remain largely unknown. Here, we report a targeted direct rRNA sequencing approach involving a substrate selection step and demonstrate its suitability to identify differential modification sites in combination with the JACUSA2 software. We compared JACUSA2 to other tools designed for RNA modification detection and show that JACUSA2 outperforms other software with regard to detection of base modifications such as methylation, acetylation and aminocarboxypropylation. To illustrate its widespread usability, we applied our method to a collection of CRISPR-Cas9 engineered colon carcinoma cells lacking specific enzymatic activities responsible for particular rRNA modifications and systematically compared them to isogenic wild-type RNAs. Besides the numerous 2'-O methylated riboses and pseudouridylated residues, our approach was suitable to reliably identify differential base methylation and acetylation events. Importantly, our method does not require any prior knowledge of modification sites or the need to train complex models. We further report for the first time detection of human rRNA modifications by direct RNA-sequencing on Flongle flow cells, the smallest-scale nanopore flow cell available to date. The use of these smaller flow cells reduces RNA input requirements, making our workflow suitable for the analysis of samples with limited availability and clinical work.

Data Analysis Pipeline for Detection and Quantification of Pseudouridine (ψ) in RNA by HydraPsiSeq

Pseudouridine, a modified RNA residue formed by the isomerization of its parental U nucleotide, is prevalent in a majority of cellular RNAs; its presence was reported in tRNA, rRNA, and sn/snoRNA as well as in mRNA/lncRNA. Multiple analytical deep sequencing-based approaches have been proposed for pseudouridine detection and quantification, among which the most popular relies on the use of soluble carbodiimide (termed CMCT). Recently, we developed an alternative protocol for pseudouridine mapping and quantification. The principle is based on protection of pseudouridine against random RNA cleavage by hydrazine/aniline treatment (HydraPsiSeq protocol). This "negative" detection mode requires higher sequencing depth and provides a precise quantification of the pseudouridine content. All "wet-lab" technical details of the HydraPsiSeq protocol have been described in recent publications. Here, we describe all bioinformatics analysis steps required for data processing from raw reads to the pseudouridylation profile of known or unknown RNA.

Maturation of small nucleolar RNAs: from production to function

Small Nucleolar RNAs (snoRNAs) are an abundant group of non-coding RNAs with well-defined roles in ribosomal RNA processing, folding and chemical modification. Besides their classic roles in ribosome biogenesis, snoRNAs are also implicated in several other cellular activities including regulation of splicing, transcription, RNA editing, cellular trafficking, and miRNA-like functions. Mature snoRNAs must undergo a series of processing steps tightly regulated by transiently associating factors and coordinated with other cellular processes including transcription and splicing. In addition to their mature forms, snoRNAs can contribute to gene expression regulation through their derivatives and degradation products. Here, we review the current knowledge on mechanisms of snoRNA maturation, including the different pathways of processing, and the regulatory mechanisms that control snoRNA levels and complex assembly. We also discuss the significance of studying snoRNA maturation, highlight the gaps in the current knowledge and suggest directions for future research in this area.

Detection technologies for RNA modifications

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

RNA modifications in aging-associated cardiovascular diseases

Cardiovascular disease (CVD) is a leading cause of morbidity and mortality worldwide that bears an enormous healthcare burden and aging is a major contributing factor to CVDs. Functional gene expression network during aging is regulated by mRNAs transcriptionally and by non-coding RNAs epi-transcriptionally. RNA modifications alter the stability and function of both mRNAs and non-coding RNAs and are involved in differentiation, development, and diseases. Here we review major chemical RNA modifications on mRNAs and non-coding RNAs, including N6-adenosine methylation, N1-adenosine methylation, 5-methylcytidine, pseudouridylation, 2′ -O-ribose-methylation, and N7-methylguanosine, in the aging process with an emphasis on cardiovascular aging. We also summarize the currently available methods to detect RNA modifications and the bioinformatic tools to study RNA modifications. More importantly, we discussed the specific implication of the RNA modifications on mRNAs and non-coding RNAs in the pathogenesis of aging-associated CVDs, including atherosclerosis, hypertension, coronary heart diseases, congestive heart failure, atrial fibrillation, peripheral artery disease, venous insufficiency, and stroke.

Lytic Reactivation of the Kaposi’s Sarcoma-Associated Herpesvirus (KSHV) Is Accompanied by Major Nucleolar Alterations

The nucleolus is a subnuclear compartment whose primary function is the biogenesis of ribosomal subunits. Certain viral infections affect the morphology and composition of the nucleolar compartment and influence ribosomal RNA (rRNA) transcription and maturation. However, no description of nucleolar morphology and function during infection with Kaposi’s sarcoma-associated herpesvirus (KSHV) is available to date. Using immunofluorescence microscopy, we documented extensive destruction of the nuclear and nucleolar architecture during the lytic reactivation of KSHV. This was manifested by the redistribution of key nucleolar proteins, including the rRNA transcription factor UBF. Distinct delocalization patterns were evident; certain nucleolar proteins remained together whereas others dissociated, implying that nucleolar proteins undergo nonrandom programmed dispersion. Significantly, the redistribution of UBF was dependent on viral DNA replication or late viral gene expression. No significant changes in pre-rRNA levels and no accumulation of pre-rRNA intermediates were found by RT-qPCR and Northern blot analysis. Furthermore, fluorescent in situ hybridization (FISH), combined with immunofluorescence, revealed an overlap between Fibrillarin and internal transcribed spacer 1 (ITS1), which represents the primary product of the pre-rRNA, suggesting that the processing of rRNA proceeds during lytic reactivation. Finally, small changes in the levels of pseudouridylation (Ψ) and 2′-O-methylation (Nm) were documented across the rRNA; however, none were localized to the functional domain. Taken together, our results suggest that despite dramatic changes in the nucleolar organization, rRNA transcription and processing persist during lytic reactivation of KSHV. Whether the observed nucleolar alterations favor productive infection or signify cellular anti-viral responses remains to be determined.

Cloning and Sequencing Eukaryotic Small RNAs

Small RNAs are ubiquitous regulators of gene expression that participate in nearly all aspects of physiology in a wide range of organisms. There are many different classes of eukaryotic small RNAs that play regulatory roles at every level of gene expression, including transcription, RNA stability, and translation. While eukaryotic small RNAs display diverse functions across and within classes, they are generally grouped functionally based on the machinery required for their biogenesis, the effector proteins they associate with, and their molecular characteristics. The development of techniques to clone and sequence small RNAs has been critical for their identification, yet the ligation-dependent addition of RNA adapters and the use of reverse transcriptase to generate cDNA in traditional library preparation protocols can be unsuitable to detect certain small RNA subtypes. In particular, 3' or 5' chemical modifications that are characteristic of specific types of small RNAs can impede the ligation-dependent addition of RNA adapters, while internal RNA modifications can interfere with accurate reverse transcription. The inability to clone certain small RNA subtypes with traditional protocols results in an inaccurate assessment of small RNA abundance and diversity, where some RNAs appear over-represented and others are not detected. This overview aims to guide users on how to design small RNA cloning workflows in eukaryotes to more accurately capture specific small RNAs of interest. Hence, we discuss the molecular biology underlying the identification and quantitation of small RNAs, explore the limitations of commonly used protocols, and detail the alternative approaches that can be used to enrich specific small RNA classes. © 2022 Wiley Periodicals LLC.

The plant epitranscriptome: revisiting pseudouridine and 2'-O-methyl RNA modifications

There is growing evidence that post-transcriptional RNA modifications are highly dynamic and can be used to improve crop production. Although more than 172 unique types of RNA modifications have been identified throughout the kingdom of life, we are yet to leverage upon the understanding to optimize RNA modifications in crops to improve productivity. The contributions of internal mRNA modifications such as N6-methyladenosine (m⁶ A) and 5-methylcytosine (m⁵ C) methylations to embryonic development, root development, leaf morphogenesis, flowering, fruit ripening and stress response are sufficiently known, but the roles of the two most abundant RNA modifications, pseudouridine (Ψ) and 2'-O-methylation (Nm), in the cell remain unclear due to insufficient advances in high-throughput technologies in plant development. Therefore, in this review, we discuss the latest methods and insights gained in mapping internal Ψ and Nm and their unique properties in plants and other organisms. In addition, we discuss the limitations that remain in high-throughput technologies for qualitative and quantitative mapping of these RNA modifications and highlight future challenges in regulating the plant epitranscriptome.

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

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

Fractional 2'-O-methylation in the ribosomal RNA of Dictyostelium discoideum supports ribosome heterogeneity in Amoebozoa

A hallmark of ribosomal RNA (rRNA) are 2′-O-methyl groups that are introduced sequence specifically by box C/D small nucleolar RNAs (snoRNAs) in ribonucleoprotein particles. Most data on this chemical modification and its impact on RNA folding and stability are derived from organisms of the Opisthokonta supergroup. Using bioinformatics and RNA-seq data, we identify 30 novel box C/D snoRNAs in Dictyostelium discoideum, many of which are differentially expressed during the multicellular development of the amoeba. By applying RiboMeth-seq, we find 49 positions in the 17S and 26S rRNA 2′-O-methylated. Several of these nucleotides are substoichiometrically modified, with one displaying dynamic modification levels during development. Using homology-based models for the D. discoideum rRNA secondary structures, we localize many modified nucleotides in the vicinity of the ribosomal A, P and E sites. For most modified positions, a guiding box C/D snoRNA could be identified, allowing to determine idiosyncratic features of the snoRNA/rRNA interactions in the amoeba. Our data from D. discoideum represents the first evidence for ribosome heterogeneity in the Amoebozoa supergroup, allowing to suggest that it is a common feature of all eukaryotes.

Specialised ribosomes as versatile regulators of gene expression

The ribosome has long been thought to be a homogeneous cellular machine that constitutively and globally synthesises proteins from mRNA. However, recent studies have revealed that ribosomes are highly heterogeneous, dynamic macromolecular complexes with specialised roles in translational regulation in many organisms across the kingdoms. In this review, we summarise the current understanding of ribosome heterogeneity and the specialised functions of heterogeneous ribosomes. We also discuss specialised translation systems that utilise orthogonal ribosomes.

Systematic mapping of rRNA 2'-O methylation during frog development and involvement of the methyltransferase Fibrillarin in eye and craniofacial development in Xenopus laevis

Ribosomes are essential nanomachines responsible for protein production. Although ribosomes are present in every living cell, ribosome biogenesis dysfunction diseases, called ribosomopathies, impact particular tissues specifically. Here, we evaluate the importance of the box C/D snoRNA-associated ribosomal RNA methyltransferase fibrillarin (Fbl) in the early embryonic development of Xenopus laevis. We report that in developing embryos, the neural plate, neural crest cells (NCCs), and NCC derivatives are rich in fbl transcripts. Fbl knockdown leads to striking morphological defects affecting the eyes and craniofacial skeleton, due to lack of NCC survival caused by massive p53-dependent apoptosis. Fbl is required for efficient pre-rRNA processing and 18S rRNA production, which explains the early developmental defects. Using RiboMethSeq, we systematically reinvestigated ribosomal RNA 2’-O methylation in X. laevis, confirming all 89 previously mapped sites and identifying 15 novel putative positions in 18S and 28S rRNA. Twenty-three positions, including 10 of the new ones, were validated orthogonally by low dNTP primer extension. Bioinformatic screening of the X. laevis transcriptome revealed candidate box C/D snoRNAs for all methylated positions. Mapping of 2’-O methylation at six developmental stages in individual embryos indicated a trend towards reduced methylation at specific positions during development. We conclude that fibrillarin knockdown in early Xenopus embryos causes reduced production of functional ribosomal subunits, thus impairing NCC formation and migration.

Ribosomes are essential nanomachines responsible for protein production in all cells. Ribosomopathies are diseases caused by improper ribosome formation due to mutations in ribosomal proteins or ribosome assembly factors. Such diseases primarily affect the brain and blood, and it is unclear how malfunctioning of a process as general as ribosome formation can lead to tissue-specific diseases. Here we have examined how fibrillarin, an enzyme which modifies ribosomal RNA by adding methyl groups at specific sites, affects early embryonic development in the frog Xenopus laevis. We have revealed its importance in the maturation of cells forming an embryonic structure called the neural crest. Fibrillarin depletion leads to reduced eye size and abnormal head shape, reminiscent of other conditions such as Treacher Collins syndrome. Molecularly, the observed phenotypes are explainable by increased p53-dependent programmed cell death triggered by inhibition of certain pre-rRNA processing steps. Our systematic investigation of the ribosomal RNA 2’-O methylation repertoire across development has further revealed hypomodification at a late stage of development, which might play a role in late developmental transitions involving differential translation by compositionally different ribosomes.

Genome instability drives epistatic adaptation in the human pathogen Leishmania

Chromosome and gene copy number variations often correlate with the evolution of microbial and cancer drug resistance, thus causing important human mortality. How genome instability is harnessed to generate beneficial phenotypes and how deleterious gene dosage effects are compensated remain open questions. The protist pathogen Leishmania exploits genome instability to regulate expression via gene dosage changes. Using these parasites as a unique model system, we uncover complex epistatic interactions between gene copy number variations and compensatory transcriptomic responses as key processes that harness genome instability for adaptive evolution in Leishmania. Our results propose a model of eukaryotic fitness gain that may be broadly applicable to pathogenic fungi or tumor cells known to exploit genome instability for adaptation.

How genome instability is harnessed for fitness gain despite its potential deleterious effects is largely elusive. An ideal system to address this important open question is provided by the protozoan pathogen Leishmania, which exploits frequent variations in chromosome and gene copy number to regulate expression levels. Using ecological genomics and experimental evolution approaches, we provide evidence that Leishmania adaptation relies on epistatic interactions between functionally associated gene copy number variations in pathways driving fitness gain in a given environment. We further uncover posttranscriptional regulation as a key mechanism that compensates for deleterious gene dosage effects and provides phenotypic robustness to genetically heterogenous parasite populations. Finally, we correlate dynamic variations in small nucleolar RNA (snoRNA) gene dosage with changes in ribosomal RNA 2′-O-methylation and pseudouridylation, suggesting translational control as an additional layer of parasite adaptation. Leishmania genome instability is thus harnessed for fitness gain by genome-dependent variations in gene expression and genome-independent compensatory mechanisms. This allows for polyclonal adaptation and maintenance of genetic heterogeneity despite strong selective pressure. The epistatic adaptation described here needs to be considered in Leishmania epidemiology and biomarker discovery and may be relevant to other fast-evolving eukaryotic cells that exploit genome instability for adaptation, such as fungal pathogens or cancer.

DHX15-independent roles for TFIP11 in U6 snRNA modification, U4/U6.U5 tri-snRNP assembly and pre-mRNA splicing fidelity

The U6 snRNA, the core catalytic component of the spliceosome, is extensively modified post-transcriptionally, with 2'-O-methylation being most common. However, how U6 2'-O-methylation is regulated remains largely unknown. Here we report that TFIP11, the human homolog of the yeast spliceosome disassembly factor Ntr1, localizes to nucleoli and Cajal Bodies and is essential for the 2'-O-methylation of U6. Mechanistically, we demonstrate that TFIP11 knockdown reduces the association of U6 snRNA with fibrillarin and associated snoRNAs, therefore altering U6 2'-O-methylation. We show U6 snRNA hypomethylation is associated with changes in assembly of the U4/U6.U5 tri-snRNP leading to defects in spliceosome assembly and alterations in splicing fidelity. Strikingly, this function of TFIP11 is independent of the RNA helicase DHX15, its known partner in yeast. In sum, our study demonstrates an unrecognized function for TFIP11 in U6 snRNP modification and U4/U6.U5 tri-snRNP assembly, identifying TFIP11 as a critical spliceosome assembly regulator.

Regulation of translation by site-specific ribosomal RNA methylation

Ribosomes are complex ribozymes that interpret genetic information by translating messenger RNA (mRNA) into proteins. Natural variation in ribosome composition has been documented in several organisms and can arise from several different sources. A key question is whether specific control over ribosome heterogeneity represents a mechanism by which translation can be regulated. We used RiboMeth-seq to demonstrate that differential 2'-O-methylation of ribosomal RNA (rRNA) represents a considerable source of ribosome heterogeneity in human cells, and that modification levels at distinct sites can change dynamically in response to upstream signaling pathways, such as MYC oncogene expression. Ablation of one prominent methylation resulted in altered translation of select mRNAs and corresponding changes in cellular phenotypes. Thus, differential rRNA 2'-O-methylation can give rise to ribosomes with specialized function. This suggests a broader mechanism where the specific regulation of rRNA modification patterns fine tunes translation.

The ribosome epitranscriptome: inert-or a platform for functional plasticity?

A universal property of all rRNAs explored to date is the prevalence of post-transcriptional ("epitranscriptional") modifications, which expand the chemical and topological properties of the four standard nucleosides. Are these modifications an inert, constitutive part of the ribosome? Or could they, in part, also regulate the structure or function of the ribosome? In this review, we summarize emerging evidence that rRNA modifications are more heterogeneous than previously thought, and that they can also vary from one condition to another, such as in the context of a cellular response or a developmental trajectory. We discuss the implications of these results and key open questions on the path toward connecting such heterogeneity with function.

Association of snR190 snoRNA chaperone with early pre-60S particles is regulated by the RNA helicase Dbp7 in yeast

Synthesis of eukaryotic ribosomes involves the assembly and maturation of precursor particles (pre-ribosomal particles) containing ribosomal RNA (rRNA) precursors, ribosomal proteins (RPs) and a plethora of assembly factors (AFs). Formation of the earliest precursors of the 60S ribosomal subunit (pre-60S r-particle) is among the least understood stages of ribosome biogenesis. It involves the Npa1 complex, a protein module suggested to play a key role in the early structuring of the pre-rRNA. Npa1 displays genetic interactions with the DExD-box protein Dbp7 and interacts physically with the snR190 box C/D snoRNA. We show here that snR190 functions as a snoRNA chaperone, which likely cooperates with the Npa1 complex to initiate compaction of the pre-rRNA in early pre-60S r-particles. We further show that Dbp7 regulates the dynamic base-pairing between snR190 and the pre-rRNA within the earliest pre-60S r-particles, thereby participating in structuring the peptidyl transferase center (PTC) of the large ribosomal subunit.

The molecular events underlying the assembly and maturation of the early pre-60S particles during eukaryotic ribosome synthesis are not well understood. Here, the authors combine yeast genetics and biochemical experiments to characterise the functions of two important players of eukaryotic ribosome biogenesis, the box C/D snoRNP snR190 and the helicase Dbp7, which both interact. They show that the snR190 snoRNA acts as a RNA chaperone that assists the structuring of the 25S rRNA during the maturation of early pre-60S particles and that Dbp7 is important for facilitating remodeling events in the peptidyl transferase center region of the 25S rRNAs during the maturation of early pre-60S particles.

The RNA helicase Dbp7 promotes domain V/VI compaction and stabilization of inter-domain interactions during early 60S assembly

Early pre-60S ribosomal particles are poorly characterized, highly dynamic complexes that undergo extensive rRNA folding and compaction concomitant with assembly of ribosomal proteins and exchange of assembly factors. Pre-60S particles contain numerous RNA helicases, which are likely regulators of accurate and efficient formation of appropriate rRNA structures. Here we reveal binding of the RNA helicase Dbp7 to domain V/VI of early pre-60S particles in yeast and show that in the absence of this protein, dissociation of the Npa1 scaffolding complex, release of the snR190 folding chaperone, recruitment of the A3 cluster factors and binding of the ribosomal protein uL3 are impaired. uL3 is critical for formation of the peptidyltransferase center (PTC) and is responsible for stabilizing interactions between the 5′ and 3′ ends of the 25S, an essential pre-requisite for subsequent pre-60S maturation events. Highlighting the importance of pre-ribosome remodeling by Dbp7, our data suggest that in the absence of Dbp7 or its catalytic activity, early pre-ribosomal particles are targeted for degradation.

Early steps of large 60S ribosomal subunit biogenesis are not well understood. Here, the authors combine biochemical experiments with protein-RNA crosslinking and mass spectrometry to show that the RNA helicase Dbp7 is key player during early 60S ribosomal assembly. Dbp7 regulates a series of events driving compaction of domain V/VI in early pre60S ribosomal particles.

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.

A novel view on an old drug, 5-fluorouracil: an unexpected RNA modifier with intriguing impact on cancer cell fate

5-Fluorouracil (5-FU) is a chemotherapeutic drug widely used to treat patients with solid tumours, such as colorectal and pancreatic cancers. Colorectal cancer (CRC) is the second leading cause of cancer-related death and half of patients experience tumour recurrence. Used for over 60 years, 5-FU was long thought to exert its cytotoxic effects by altering DNA metabolism. However, 5-FU mode of action is more complex than previously anticipated since 5-FU is an extrinsic source of RNA modifications through its ability to be incorporated into most classes of RNA. In particular, a recent report highlighted that, by its integration into the most abundant RNA, namely ribosomal RNA (rRNA), 5-FU creates fluorinated active ribosomes and induces translational reprogramming. Here, we review the historical knowledge of 5-FU mode of action and discuss progress in the field of 5-FU-induced RNA modifications. The case of rRNA, the essential component of ribosome and translational activity, and the plasticity of which was recently associated with cancer, is highlighted. We propose that translational reprogramming, induced by 5-FU integration in ribosomes, contributes to 5-FU-driven cell plasticity and ultimately to relapse.

LIN28 coordinately promotes nucleolar/ribosomal functions and represses the 2C-like transcriptional program in pluripotent stem cells

LIN28 is an RNA binding protein with important roles in early embryo development, stem cell differentiation/reprogramming, tumorigenesis and metabolism. Previous studies have focused mainly on its role in the cytosol where it interacts with Let-7 microRNA precursors or mRNAs, and few have addressed LIN28's role within the nucleus. Here, we show that LIN28 displays dynamic temporal and spatial expression during murine embryo development. Maternal LIN28 expression drops upon exit from the 2-cell stage, and zygotic LIN28 protein is induced at the forming nucleolus during 4-cell to blastocyst stage development, to become dominantly expressed in the cytosol after implantation. In cultured pluripotent stem cells (PSCs), loss of LIN28 led to nucleolar stress and activation of a 2-cell/4-cell-like transcriptional program characterized by the expression of endogenous retrovirus genes. Mechanistically, LIN28 binds to small nucleolar RNAs and rRNA to maintain nucleolar integrity, and its loss leads to nucleolar phase separation defects, ribosomal stress and activation of P53 which in turn binds to and activates 2C transcription factor Dux. LIN28 also resides in a complex containing the nucleolar factor Nucleolin (NCL) and the transcriptional repressor TRIM28, and LIN28 loss leads to reduced occupancy of the NCL/TRIM28 complex on the Dux and rDNA loci, and thus de-repressed Dux and reduced rRNA expression. Lin28 knockout cells with nucleolar stress are more likely to assume a slowly cycling, translationally inert and anabolically inactive state, which is a part of previously unappreciated 2C-like transcriptional program. These findings elucidate novel roles for nucleolar LIN28 in PSCs, and a new mechanism linking 2C program and nucleolar functions in PSCs and early embryo development.

Nopp140-chaperoned 2'-O-methylation of small nuclear RNAs in Cajal bodies ensures splicing fidelity

In this study, Bizarro et al. sought to understand the function and subcellular site of snRNA modification, and found that Cajal body (CB) localization of the protein Nopp140 is essential for concentration of small Cajal body-specific ribonucleoproteins (scaRNPs) in nuclear condensate and that phosphorylation by casein kinase 2 (CK2) at ∼80 serines targets Nopp140 to CBs. Nopp140 knockdown-mediated release of scaRNPs from CBs severely compromises 2′-O-methylation of spliceosomal snRNAs, identifying CBs as the site of scaRNP catalysis.

Spliceosomal small nuclear RNAs (snRNAs) are modified by small Cajal body (CB)-specific ribonucleoproteins (scaRNPs) to ensure snRNP biogenesis and pre-mRNA splicing. However, the function and subcellular site of snRNA modification are largely unknown. We show that CB localization of the protein Nopp140 is essential for concentration of scaRNPs in that nuclear condensate; and that phosphorylation by casein kinase 2 (CK2) at ∼80 serines targets Nopp140 to CBs. Transiting through CBs, snRNAs are apparently modified by scaRNPs. Indeed, Nopp140 knockdown-mediated release of scaRNPs from CBs severely compromises 2′-O-methylation of spliceosomal snRNAs, identifying CBs as the site of scaRNP catalysis. Additionally, alternative splicing patterns change indicating that these modifications in U1, U2, U5, and U12 snRNAs safeguard splicing fidelity. Given the importance of CK2 in this pathway, compromised splicing could underlie the mode of action of small molecule CK2 inhibitors currently considered for therapy in cholangiocarcinoma, hematological malignancies, and COVID-19.

Single-Cell Imaging of m6 A Modified RNA Using m6 A-Specific In Situ Hybridization Mediated Proximity Ligation Assay (m6 AISH-PLA)

N⁶ -methyladenosine (m⁶ A) modification-the most prevalent mammalian RNA internal modification-plays key regulatory roles in mRNA metabolism. Current approaches for m⁶ A modified RNA analysis limit at bulk-population level, resulting in a loss of spatiotemporal and cell-to-cell variability information. Here we proposed a m⁶ A-specific in situ hybridization mediated proximity ligation assay (m⁶ AISH-PLA) for cellular imaging of m⁶ A RNA, allowing to identify m⁶ A modification at specific location in RNAs and image m⁶ A RNA with single-cell and single-molecule resolution. Using m⁶ AISH-PLA, we investigated the m⁶ A level and subcellular location of HSP70 RNA103-m⁶ A in response to heat shock stress, and found an increased m⁶ A modified ratio and an increased distribution ratio in cytoplasm under heat shock. m⁶ AISH-PLA can serve in the study of m⁶ A RNA in single cells for deciphering epitranscriptomic mechanisms and assisting clinical diagnosis.

Molecular mechanism of RNase R substrate sensitivity for RNA ribose methylation

RNA 2′-O-methylation is widely distributed and plays important roles in various cellular processes. Mycoplasma genitalium RNase R (MgR), a prokaryotic member of the RNase II/RNB family, is a 3′-5′ exoribonuclease and is particularly sensitive to RNA 2′-O-methylation. However, how RNase R interacts with various RNA species and exhibits remarkable sensitivity to substrate 2′-O-methyl modifications remains elusive. Here we report high-resolution crystal structures of MgR in apo form and in complex with various RNA substrates. The structural data together with extensive biochemical analysis quantitively illustrate MgR’s ribonuclease activity and significant sensitivity to RNA 2′-O-methylation. Comparison to its related homologs reveals an exquisite mechanism for the recognition and degradation of RNA substrates. Through structural and mutagenesis studies, we identified proline 277 to be responsible for the significant sensitivity of MgR to RNA 2′-O-methylation within the RNase II/RNB family. We also generated several MgR variants with modulated activities. Our work provides a mechanistic understanding of MgR activity that can be harnessed as a powerful RNA analytical tool that will open up a new venue for RNA 2′-O-methylations research in biological and clinical samples.

Shaping the Bacterial Epitranscriptome-5'-Terminal and Internal RNA Modifications

All domains of life utilize a diverse set of modified ribonucleotides that can impact the sequence, structure, function, stability, and the fate of RNAs, as well as their interactions with other molecules. Today, more than 160 different RNA modifications are known that decorate the RNA at the 5'-terminus or internal RNA positions. The boost of next-generation sequencing technologies sets the foundation to identify and study the functional role of RNA modifications. The recent advances in the field of RNA modifications reveal a novel regulatory layer between RNA modifications and proteins, which is central to developing a novel concept called "epitranscriptomics." The majority of RNA modifications studies focus on the eukaryotic epitranscriptome. In contrast, RNA modifications in prokaryotes are poorly characterized. This review outlines the current knowledge of the prokaryotic epitranscriptome focusing on mRNA modifications. Here, it is described that several internal and 5'-terminal RNA modifications either present or likely present in prokaryotic mRNA. Thereby, the individual techniques to identify these epitranscriptomic modifications, their writers, readers and erasers, and their proposed functions are explored. Besides that, still unanswered questions in the field of prokaryotic epitranscriptomics are pointed out, and its future perspectives in the dawn of next-generation sequencing technologies are outlined.

Quantitative profiling of pseudouridylation dynamics in native RNAs with nanopore sequencing

Nanopore RNA sequencing shows promise as a method for discriminating and identifying different RNA modifications in native RNA. Expanding on the ability of nanopore sequencing to detect N⁶-methyladenosine, we show that other modifications, in particular pseudouridine (Ψ) and 2'-O-methylation (Nm), also result in characteristic base-calling 'error' signatures in the nanopore data. Focusing on Ψ modification sites, we detected known and uncovered previously unreported Ψ sites in mRNAs, non-coding RNAs and rRNAs, including a Pus4-dependent Ψ modification in yeast mitochondrial rRNA. To explore the dynamics of pseudouridylation, we treated yeast cells with oxidative, cold and heat stresses and detected heat-sensitive Ψ-modified sites in small nuclear RNAs, small nucleolar RNAs and mRNAs. Finally, we developed a software, nanoRMS, that estimates per-site modification stoichiometries by identifying single-molecule reads with altered current intensity and trace profiles. This work demonstrates that Nm and Ψ RNA modifications can be detected in cellular RNAs and that their modification stoichiometry can be quantified by nanopore sequencing of native RNA.

RNA helicase-mediated regulation of snoRNP dynamics on pre-ribosomes and rRNA 2'-O-methylation

RNA helicases play important roles in diverse aspects of RNA metabolism through their functions in remodelling ribonucleoprotein complexes (RNPs), such as pre-ribosomes. Here, we show that the DEAD box helicase Dbp3 is required for efficient processing of the U18 and U24 intron-encoded snoRNAs and 2′-O-methylation of various sites within the 25S ribosomal RNA (rRNA) sequence. Furthermore, numerous box C/D snoRNPs accumulate on pre-ribosomes in the absence of Dbp3. Many snoRNAs guiding Dbp3-dependent rRNA modifications have overlapping pre-rRNA basepairing sites and therefore form mutually exclusive interactions with pre-ribosomes. Analysis of the distribution of these snoRNAs between pre-ribosome-associated and ‘free’ pools demonstrated that many are almost exclusively associated with pre-ribosomal complexes. Our data suggest that retention of such snoRNPs on pre-ribosomes when Dbp3 is lacking may impede rRNA 2′-O-methylation by reducing the recycling efficiency of snoRNPs and by inhibiting snoRNP access to proximal target sites. The observation of substoichiometric rRNA modification at adjacent sites suggests that the snoRNPs guiding such modifications likely interact stochastically rather than hierarchically with their pre-rRNA target sites. Together, our data provide new insights into the dynamics of snoRNPs on pre-ribosomal complexes and the remodelling events occurring during the early stages of ribosome assembly.

Profiling of RNA ribose methylation in Arabidopsis thaliana

Eukaryotic rRNAs and snRNAs are decorated with abundant 2′-O-methylated nucleotides (Nm) that are predominantly synthesized by box C/D snoRNA-guided enzymes. In the model plant Arabidopsis thaliana, C/D snoRNAs have been well categorized, but there is a lack of systematic mapping of Nm. Here, we applied RiboMeth-seq to profile Nm in cytoplasmic, chloroplast and mitochondrial rRNAs and snRNAs. We identified 111 Nm in cytoplasmic rRNAs and 19 Nm in snRNAs and assigned guide for majority of the detected sites using an updated snoRNA list. At least four sites are directed by guides with multiple specificities as shown in yeast. We found that C/D snoRNAs frequently form extra pairs with nearby sequences of methylation sites, potentially facilitating the substrate binding. Chloroplast and mitochondrial rRNAs contain five almost identical methylation sites, including two novel sites mediating ribosomal subunit joining. Deletion of FIB1 or FIB2 gene reduced the accumulation of C/D snoRNA and rRNA methylation with FIB1 playing a bigger role in methylation. Our data reveal the comprehensive 2′-O-methylation maps for Arabidopsis rRNAs and snRNAs and would facilitate study of their function and biosynthesis.

Protective Effect of Blood Cora Polysaccharides on H9c2 Rat Heart Cells Injury Induced by Oxidative Stress by Activating Nrf2/HO-1 Signal Pathway

The protective effect of blood cora polysaccharides (BCP) on H9c2 rat heart cells under oxidative stress was explored with the use of a H9c2 cell oxidative stress model. The ability of BCP to scavenge 2,2'-azino-bis(3-ethylbenzothiazoline-6-sulfonic acid) (ABTS), 1,1-diphenyl-2-picrylhydrazyl (DPPH), and hydroxyl radicals and its reducing power were measured in vitro, indicating a more powerful antioxidant effect of BCP compared to a similar concentration of vitamin C. The cellular metabolic activity was tested through the MTT [3-(4,5-dimethyl-2-thiazolyl)-2,5-diphenyl-2-H-tetrazolium bromide] assay. Additionally, the relevant oxidation indicator level within the cell supernatant and cells was tested with reagent kits, and mRNA and protein expression levels in the cells were tested through quantitative polymerase chain reaction (qPCR) and western blot. The chemical composition of BCP was determined through high performance liquid chromatography (HPLC). The results show that compared with the normal group, the model group's cell survival rate (28.75 ± 2.56%) decreased, lactate dehydrogenase (LDH) leakage and the malondialdehyde (MDA) content increased, and superoxide dismutase (SOD), catalase (CAT), and glutathione (GSH) levels decreased. The results of qPCR and western blot show that compared with the normal group, the model group's Bcl-2 associated X protein (Bax), caspase-3, nuclear factor erythroid-2 related factor 2 (Nrf2), heme oxygenase-1 (HO-1) expression, NAD(P)H:quinoneoxidoreductase 1 (NQO1), and cytochrome c (Cyt C) decreased, and B-cell lymphoma-2 (Bcl-2) expression was increased, with significant statistical differences. Compared with the model group, the cell survival rate for each BCP-treated group increased, the LDH leakage decreased, the SOD, CAT, and GSH levels in the cells increased, the MDA content decreased, the Bax, caspase-3, Nrf2, HO-1, NQO1, and Cyt C expression was weakened, and the Bcl-2 expression was strengthened. BCP inhibited the reduction of mitochondrial membrane potential caused by H2O2 treatment. According to the component analysis, BCP mainly consist of mannitol, ribose, glucosum anhydricum, galactose, and xylose. It was observed that the Nrf2/HO-1 signaling pathway can be activated, regulated, and controlled by functional BCP to protect H9c2 cells injured by oxidative stress.

Analysis of RNA Modifications by Second- and Third-Generation Deep Sequencing: 2020 Update

The precise mapping and quantification of the numerous RNA modifications that are present in tRNAs, rRNAs, ncRNAs/miRNAs, and mRNAs remain a major challenge and a top priority of the epitranscriptomics field. After the keystone discoveries of massive m⁶A methylation in mRNAs, dozens of deep sequencing-based methods and protocols were proposed for the analysis of various RNA modifications, allowing us to considerably extend the list of detectable modified residues. Many of the currently used methods rely on the particular reverse transcription signatures left by RNA modifications in cDNA; these signatures may be naturally present or induced by an appropriate enzymatic or chemical treatment. The newest approaches also include labeling at RNA abasic sites that result from the selective removal of RNA modification or the enhanced cleavage of the RNA ribose-phosphate chain (perhaps also protection from cleavage), followed by specific adapter ligation. Classical affinity/immunoprecipitation-based protocols use either antibodies against modified RNA bases or proteins/enzymes, recognizing RNA modifications. In this survey, we review the most recent achievements in this highly dynamic field, including promising attempts to map RNA modifications by the direct single-molecule sequencing of RNA by nanopores.