Synthesis, Function, and Heterogeneity of snoRNA-Guided Posttranscriptional Nucleoside Modifications in Eukaryotic Ribosomal RNAs

The small nucleolar RNA SNORA51 enhances breast cancer stem cell-like properties via the RPL3/NPM1/c-MYC pathway

Breast cancer stem cells (BCSCs) are key players in carcinogenesis and development. Small nucleolar RNAs (snoRNAs) seem to have a crucial influence on regulating stem cell-like properties in various cancers, but the underlying mechanism in breast cancer has not been determined. In this study, we first found that the expression of SNORA51 might be strongly and positively related to BCSCs-like properties. SNORA51 expression was assessed in breast cancer tissues (n = 158 patients) by in situ hybridization. Colony formation, cell counting kit-8, and sphere formation assays were used to detect cell proliferation and self-renewal, respectively. Wound healing and transwell assays were used to detect cell migration. Coimmunoprecipitation and molecular docking were used to determine the underlying mechanism through which SNORA51 regulates BCSCs-like properties. High SNORA51 expression was associated with a worse prognosis, overall survival, and disease-free survival, in 158 breast cancer patients and was also closely related to lymph node status, ER status, the Ki-67 index, histological grade, and TNM stage. Further analysis proved that SNORA51 could enhance and maintain stem cell-like properties, including cell proliferation, self-renewal, and migration, in breast cancer. Moreover, high SNORA51 expression could reduce nucleolar RPL3 expression, induce changes in the expression of NPM1 in the nucleolus and nucleoplasm, and ultimately increase c-MYC expression. Taken together, our findings demonstrated that SNORA51 could enhance BCSCs-like properties via the RPL3/NPM1/c-MYC pathway both in vitro and in vivo. Therefore, SNORA51 might be a significant biomarker and potential therapeutic target and might even provide a new viewpoint on the regulatory mechanism of snoRNAs in breast cancer or other malignant tumors.

TGF-β2 Induces Ribosome Activity, Alters Ribosome Composition and Inhibits IRES-Mediated Translation in Chondrocytes

Alterations in cell fate are often attributed to (epigenetic) regulation of gene expression. An emerging paradigm focuses on specialized ribosomes within a cell. However, little evidence exists for the dynamic regulation of ribosome composition and function. Here, we stimulated a chondrocytic cell line with transforming growth factor beta (TGF-β2) and mapped changes in ribosome function, composition and ribosomal RNA (rRNA) epitranscriptomics. 35S Met/Cys incorporation was used to evaluate ribosome activity. Dual luciferase reporter assays were used to assess ribosomal modus. Ribosomal RNA expression and processing were determined by RT-qPCR, while RiboMethSeq and HydraPsiSeq were used to determine rRNA modification profiles. Label-free protein quantification of total cell lysates, isolated ribosomes and secreted proteins was done by LC-MS/MS. A three-day TGF-β2 stimulation induced total protein synthesis in SW1353 chondrocytic cells and human articular chondrocytes. Specifically, TGF-β2 induced cap-mediated protein synthesis, while IRES-mediated translation was not (P53 IRES) or little affected (CrPv IGR and HCV IRES). Three rRNA post-transcriptional modifications (PTMs) were affected by TGF-β2 stimulation (18S-Gm1447 downregulated, 18S-ψ1177 and 28S-ψ4598 upregulated). Proteomic analysis of isolated ribosomes revealed increased interaction with eIF2 and tRNA ligases and decreased association of eIF4A3 and heterogeneous nuclear ribonucleoprotein (HNRNP)s. In addition, thirteen core ribosomal proteins were more present in ribosomes from TGF-β2 stimulated cells, albeit with a modest fold change. A prolonged stimulation of chondrocytic cells with TGF-β2 induced ribosome activity and changed the mode of translation. These functional changes could be coupled to alterations in accessory proteins in the ribosomal proteome.

Rbp95 binds to 25S rRNA helix H95 and cooperates with the Npa1 complex during early pre-60S particle maturation

Eukaryotic ribosome synthesis involves more than 200 assembly factors, which promote ribosomal RNA (rRNA) processing, modification and folding, and assembly of ribosomal proteins. The formation and maturation of the earliest pre-60S particles requires structural remodeling by the Npa1 complex, but is otherwise still poorly understood. Here, we introduce Rbp95 (Ycr016w), a constituent of early pre-60S particles, as a novel ribosome assembly factor. We show that Rbp95 is both genetically and physically linked to most Npa1 complex members and to ribosomal protein Rpl3. We demonstrate that Rbp95 is an RNA-binding protein containing two independent RNA-interacting domains. In vivo, Rbp95 associates with helix H95 in the 3′ region of the 25S rRNA, in close proximity to the binding sites of Npa1 and Rpl3. Additionally, Rbp95 interacts with several snoRNAs. The absence of Rbp95 results in alterations in the protein composition of early pre-60S particles. Moreover, combined mutation of Rbp95 and Npa1 complex members leads to a delay in the maturation of early pre-60S particles. We propose that Rbp95 acts together with the Npa1 complex during early pre-60S maturation, potentially by promoting pre-rRNA folding events within pre-60S particles.

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.

A disease-linked lncRNA mutation in RNase MRP inhibits ribosome synthesis

RMRP encodes a non-coding RNA forming the core of the RNase MRP ribonucleoprotein complex. Mutations cause Cartilage Hair Hypoplasia (CHH), characterized by skeletal abnormalities and impaired T cell activation. Yeast RNase MRP cleaves a specific site in the pre-ribosomal RNA (pre-rRNA) during ribosome synthesis. CRISPR-mediated disruption of RMRP in human cells lines caused growth arrest, with pre-rRNA accumulation. Here, we analyzed disease-relevant primary cells, showing that mutations in RMRP impair mouse T cell activation and delay pre-rRNA processing. Patient-derived human fibroblasts with CHH-linked mutations showed similar pre-rRNA processing delay. Human cells engineered with the most common CHH mutation (70AG in RMRP) show specifically impaired pre-rRNA processing, resulting in reduced mature rRNA and a reduced ratio of cytosolic to mitochondrial ribosomes. Moreover, the 70AG mutation caused a reduction in intact RNase MRP complexes. Together, these results indicate that CHH is a ribosomopathy.

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.

Constitutive and variable 2'-O-methylation (Nm) in human ribosomal RNA

Epitranscriptomic modifications of stable RNAs are dynamically regulated and specific profiles of 2ʹ-O-methylation in rRNA have been associated with distinct cancer types. However, these observations pointed out the existence of at least two distinct groups: a rather large group with constitutive rRNA Nm residues exhibiting a stable level of methylation and a more restricted set of variable modifications, giving rise to the concept of ‘specialized ribosomes’. These heterogeneous ribosomes can modulate their translational properties and be key regulatory players, depending on the physiological state of the cell. However, these conclusions were drawn from a limited set of explored human cell lines or tissues, mostly related to cancer cells of the same type. Here, we report a comprehensive analysis of human rRNA Nm modification variability observed for >15 human cell lines grown in different media and conditions. Our data demonstrate that human Nm sites can be classified into four groups, depending on their observed variability. About ⅓ of rRNA 2ʹ-O-methylations are almost invariably modified at the same level in all tested samples (stable modifications), the second group of relatively invariant modifications (another ½ of the total) showing a slightly higher variance (low variable group) and two variable groups, showing an important heterogeneity. Mapping of these four classes on the human ribosome 3D structure shows that stably modified positions are preferentially located in the important ribosome functional sites, while variable and highly variable residues are mostly distributed to the ribosome periphery. Possible relationships of such stable and variable modifications to the ribosome functions are discussed.

Identification of Transcriptomic Differences between Lower Extremities Arterial Disease, Abdominal Aortic Aneurysm and Chronic Venous Disease in Peripheral Blood Mononuclear Cells Specimens

Several human tissues are investigated in studies of molecular biomarkers associated with diseases development. Special attention is focused on the blood and its components due to combining abundant information about systemic responses to pathological processes as well as high accessibility. In the current study, transcriptome profiles of peripheral blood mononuclear cells (PBMCs) were used to compare differentially expressed genes between patients with lower extremities arterial disease (LEAD), abdominal aortic aneurysm (AAA) and chronic venous disease (CVD). Gene expression patterns were generated using the Ion S5XL next-generation sequencing platform and were analyzed using DESeq2 and UVE-PLS methods implemented in R programming software. In direct pairwise analysis, 21, 58 and 10 differentially expressed genes were selected from the comparison of LEAD vs. AAA, LEAD vs. CVD and AAA vs. CVD patient groups, respectively. Relationships between expression of dysregulated genes and age, body mass index, creatinine levels, hypertension and medication were identified using Spearman rank correlation test and two-sided Mann-Whitney U test. The functional analysis, performed using DAVID website tool, provides potential implications of selected genes in pathological processes underlying diseases studied. Presented research provides new insight into differences of pathogenesis in LEAD, AAA and CVD, and selected genes could be considered as potential candidates for biomarkers useful in diagnosis and differentiation of studied diseases.

The cardiac methylome: A hidden layer of RNA modifications to regulate gene expression

Post-transcriptional RNA modification has been observed in all kingdoms of life and more than a hundred different types of RNA modifications decorate the chemical and topological properties of these ribose nucleotides. These RNA modifications can potentially alter the RNA structure and also affect the binding affinity of proteins, thus regulating the mRNA stability as well as translation. Emerging evidence suggest that these modifications are not static, but are dynamic; vary upon different cues and are cell-type or tissue-specific. The cardiac transcriptome is not exceptional to such RNA modifications and is enriched with the abundant base methylation such as N⁶-methyladenosine (m⁶A) and also 2'-O-Methylation (Nm). In this review we will focus on the technologies available to map these modifications and as well as the contribution of these post-transcriptional modifications during various pathological conditions of the heart.

H/ACA snoRNA levels are regulated during stem cell differentiation

H/ACA small nucleolar RNAs (snoRNAs) guide pseudouridylation as part of a small nucleolar ribonucleoprotein complex (snoRNP). Disruption of H/ACA snoRNA levels in stem cells impairs pluripotency, yet it remains unclear how H/ACA snoRNAs contribute to differentiation. To determine if H/ACA snoRNA levels are dynamic during differentiation, we comprehensively profiled H/ACA snoRNA abundance in multiple murine cell types and during differentiation in three cellular models, including mouse embryonic stem cells and mouse myoblasts. We determined that the profiles of H/ACA snoRNA abundance are cell-type specific, and we identified a subset of snoRNAs that are specifically regulated during differentiation. Additionally, we demonstrated that a decrease in Snora27 abundance upon differentiation corresponds to a decrease in pseudouridylation of its target site within the E-site transfer RNA (tRNA) binding region of the 28S ribosomal RNA (rRNA) in the large ribosomal subunit. Together, these data point toward a potential model in which H/ACA snoRNAs are specifically regulated during differentiation to alter pseudouridylation and fine tune ribosome function.

Functional organization of box C/D RNA-guided RNA methyltransferase

Box C/D RNA protein complexes (RNPs) catalyze site-specific 2'-O-methylation of RNA with specificity determined by guide RNAs. In eukaryotic C/D RNP, the paralogous Nop58 and Nop56 proteins specifically associate with terminal C/D and internal C'/D' motifs of guide RNAs, respectively. We have reconstituted active C/D RNPs with recombinant proteins of the thermophilic yeast Chaetomium thermophilum. Nop58 and Nop56 could not distinguish between the two C/D motifs in the reconstituted enzyme, suggesting that the assembly specificity is imposed by trans-acting factors in vivo. The two C/D motifs are functionally independent and halfmer C/D RNAs can also guide site-specific methylation. Extensive pairing between C/D RNA and substrate is inhibitory to modification for both yeast and archaeal C/D RNPs. N6-methylated adenine at box D/D' interferes with the function of the coupled guide. Our data show that all C/D RNPs share the same functional organization and mechanism of action and provide insight into the assembly specificity of eukaryotic C/D RNPs.

Regulation and Function of RNA Pseudouridylation in Human Cells

Recent advances in pseudouridine detection reveal a complex pseudouridine landscape that includes messenger RNA and diverse classes of noncoding RNA in human cells. The known molecular functions of pseudouridine, which include stabilizing RNA conformations and destabilizing interactions with varied RNA-binding proteins, suggest that RNA pseudouridylation could have widespread effects on RNA metabolism and gene expression. Here, we emphasize how much remains to be learned about the RNA targets of human pseudouridine synthases, their basis for recognizing distinct RNA sequences, and the mechanisms responsible for regulated RNA pseudouridylation. We also examine the roles of noncoding RNA pseudouridylation in splicing and translation and point out the potential effects of mRNA pseudouridylation on protein production, including in the context of therapeutic mRNAs.

The large repertoire of 2'-O-methylation guided by C/D snoRNAs on Trypanosoma brucei rRNA

The parasite Trypanosoma brucei cycles between insect and mammalian hosts, and is the causative agent of sleeping sickness. Here, we performed genome-wide mapping of 2'-O-methylations (Nms) on trypanosome rRNA using three high-throughput sequencing methods; RibOxi-seq, RiboMeth-seq and 2'-OMe-seq. This is the first study using three genome-wide mapping approaches on rRNA from the same species showing the discrepancy among the methods. RibOxi-seq detects all the sites, but RiboMeth-seq is the only method to evaluate the level of a single Nm site. The sequencing revealed at least ninety-nine Nms guided by eighty-five snoRNAs among these thirty-eight Nms are trypanosome specific sites. We present the sequence and target of the C/D snoRNAs guiding on rRNA. This is the highest number of Nms detected to date on rRNA of a single cell parasite. Based on RiboMeth-seq, several Nm sites were found to be differentially regulated at the two stages of the parasite life cycle, the insect procyclic form (PCF) versus the bloodstream form (BSF) in the mammalian host.

Functional diversity of small nucleolar RNAs

Small nucleolar RNAs (snoRNAs) are short non-protein-coding RNAs with a long-recognized role in tuning ribosomal and spliceosomal function by guiding ribose methylation and pseudouridylation at targeted nucleotide residues of ribosomal and small nuclear RNAs, respectively. SnoRNAs are increasingly being implicated in regulation of new types of post-transcriptional processes, for example rRNA acetylation, modulation of splicing patterns, control of mRNA abundance and translational efficiency, or they themselves are processed to shorter stable RNA species that seem to be the principal or alternative bioactive isoform. Intriguingly, some display unusual cellular localization under exogenous stimuli, or tissue-specific distribution. Here, we discuss the new and unforeseen roles attributed to snoRNAs, focusing on the presumed mechanisms of action. Furthermore, we review the experimental approaches to study snoRNA function, including high resolution RNA:protein and RNA:RNA interaction mapping, techniques for analyzing modifications on targeted RNAs, and cellular and animal models used in snoRNA biology research.

Base-pairing interactions between substrate RNA and H/ACA guide RNA modulate the kinetics of pseudouridylation, but not the affinity of substrate binding by H/ACA small nucleolar ribonucleoproteins

H/ACA small nucleolar ribonucleoproteins (snoRNPs) pseudouridylate RNA in eukaryotes and archaea. They target many RNAs site-specifically through base-pairing interactions between H/ACA guide and substrate RNA. Besides ribosomal RNA (rRNA) and small nuclear RNA (snRNA), H/ACA snoRNPs are thought to also modify messenger RNA (mRNA) with potential impacts on gene expression. However, the base pairing between known target RNAs and H/ACA guide RNAs varies widely in nature, and therefore the rules governing substrate RNA selection are still not fully understood. To provide quantitative insight into substrate RNA recognition, we systematically altered the sequence of a substrate RNA target by the Saccharomyces cerevisiae H/ACA guide RNA snR34. Time courses measuring pseudouridine formation revealed a gradual decrease in the initial velocity of pseudouridylation upon reducing the number of base pairs between substrate and guide RNA. Changing or inserting nucleotides close to the target uridine severely impairs pseudouridine formation. Interestingly, filter binding experiments show that all substrate RNA variants bind to H/ACA snoRNPs with nanomolar affinity. Next, we showed that binding of inactive, near-cognate RNAs to H/ACA snoRNPs does not inhibit their activity for cognate RNAs, presumably because near-cognate RNAs dissociate rapidly. We discuss that the modulation of initial velocities by the base-pairing strength might affect the order and efficiency of pseudouridylation in rRNA during ribosome biogenesis. Moreover, the binding of H/ACA snoRNPs to near-cognate RNAs may be a mechanism to search for cognate target sites. Together, our data provide critical information to aid in the prediction of productive H/ACA guide-substrate RNA pairs.

RNA ribose methylation (2'-O-methylation): Occurrence, biosynthesis and biological functions

Methylation of riboses at 2'-OH group is one of the most common RNA modifications found in number of cellular RNAs from almost any species which belong to all three life domains. This modification was extensively studied for decades in rRNAs and tRNAs, but recent data revealed the presence of 2'-O-methyl groups also in low abundant RNAs, like mRNAs. Ribose methylation is formed in RNA by two alternative enzymatic mechanisms: either by stand-alone protein enzymes or by complex assembly of proteins associated with snoRNA guides (sno(s)RNPs). In that case one catalytic subunit acts at various RNA sites, the specificity is provided by base pairing of the sno(s)RNA guide with the target RNA. In this review we compile available information on 2'-OH ribose methylation in different RNAs, enzymatic machineries involved in their biosynthesis and dynamics, as well as on the physiological functions of these modified residues.

Pre-Ribosomal RNA Processing in Human Cells: From Mechanisms to Congenital Diseases

Ribosomal RNAs, the most abundant cellular RNA species, have evolved as the structural scaffold and the catalytic center of protein synthesis in every living organism. In eukaryotes, they are produced from a long primary transcript through an intricate sequence of processing steps that include RNA cleavage and folding and nucleotide modification. The mechanisms underlying this process in human cells have long been investigated, but technological advances have accelerated their study in the past decade. In addition, the association of congenital diseases to defects in ribosome synthesis has highlighted the central place of ribosomal RNA maturation in cell physiology regulation and broadened the interest in these mechanisms. Here, we give an overview of the current knowledge of pre-ribosomal RNA processing in human cells in light of recent progress and discuss how dysfunction of this pathway may contribute to the physiopathology of congenital diseases.

Functional Interplay between Small Non-Coding RNAs and RNA Modification in the Brain

Small non-coding RNAs are essential for transcription, translation and gene regulation in all cell types, but are particularly important in neurons, with known roles in neurodevelopment, neuroplasticity and neurological disease. Many small non-coding RNAs are directly involved in the post-transcriptional modification of other RNA species, while others are themselves substrates for modification, or are functionally modulated by modification of their target RNAs. In this review, we explore the known and potential functions of several distinct classes of small non-coding RNAs in the mammalian brain, focusing on the newly recognised interplay between the epitranscriptome and the activity of small RNAs. We discuss the potential for this relationship to influence the spatial and temporal dynamics of gene activation in the brain, and predict that further research in the field of epitranscriptomics will identify interactions between small RNAs and RNA modifications which are essential for higher order brain functions such as learning and memory.

Turnover of aberrant pre-40S pre-ribosomal particles is initiated by a novel endonucleolytic decay pathway

Ribosome biogenesis requires more than 200 trans-acting factors to achieve the correct production of the two mature ribosomal subunits. Here, we have identified Efg1 as a novel, nucleolar ribosome biogenesis factor in Saccharomyces cerevisiae that is directly linked to the surveillance of pre-40S particles. Depletion of Efg1 impairs early pre-rRNA processing, leading to a strong decrease in 18S rRNA and 40S subunit levels and an accumulation of the aberrant 23S rRNA. Using Efg1 as bait, we revealed a novel degradation pathway of the 23S rRNA. Co-immunoprecipitation experiments showed that Efg1 is a component of 90S pre-ribosomes, as it is associated with the 35S pre-rRNA and U3 snoRNA, but has stronger affinity for 23S pre-rRNA and its novel degradation intermediate 11S rRNA. 23S is cleaved at a new site, Q1, within the 18S sequence by the endonuclease Utp24, generating 11S and 17S' rRNA. Both of these cleavage products are targeted for degradation by the TRAMP/exosome complexes. Therefore, the Q1 site defines a novel endonucleolytic cleavage site of ribosomal RNA exclusively dedicated to surveillance of pre-ribosomal particles.