A computational screen for C/D box snoRNAs in the human genomic region associated with Prader-Willi and Angelman syndromes

Small nucleolar RNAs (snoRNAs) play a significant role in Prader-Willi Syndrome (PWS) and Angelman Syndrome (AS), which are genomic disorders resulting from deletions in the human chromosomal region 15q11–q13. To identify snoRNAs in the region, our computational study employs key motif features of C/D box snoRNAs and introduces a complementary RNA–RNA hybridization test. We identify three previously unknown methylation guide snoRNAs targeting ribosomal 18S and 28S RNAs, and two snoRNAs targeting serotonin receptor 2C mRNA. We show that the three snoRNA candidates likely possess methylation strands complementary to, and form stable complexes with, human ribosomal RNAs. Our screen also identifies 8 other snoRNA candidates that do not pass the rRNA-complementarity and/or hybridization tests. Two of these candidates have extensive sequence similarity to HBII-52, a snoRNA that regulates the alternative splicing of serotonin receptor 2C mRNA. Six out of our eleven candidate snoRNAs are also predicted by other existing methods.     

snoSeeker: an advanced computational package for screening of guide and orphan snoRNA genes in the human genome

Small nucleolar RNAs (snoRNAs) represent an abundant group of non-coding RNAs in eukaryotes. They can be divided into guide and orphan snoRNAs according to the presence or absence of antisense sequence to rRNAs or snRNAs. Current snoRNA-searching programs, which are essentially based on sequence complementarity to rRNAs or snRNAs, exist only for the screening of guide snoRNAs. In this study, we have developed an advanced computational package, snoSeeker, which includes CDseeker and ACAseeker programs, for the highly efficient and specific screening of both guide and orphan snoRNA genes in mammalian genomes. By using these programs, we have systematically scanned four human–mammal whole-genome alignment (WGA) sequences and identified 54 novel candidates including 26 orphan candidates as well as 266 known snoRNA genes. Eighteen novel snoRNAs were further experimentally confirmed with four snoRNAs exhibiting a tissue-specific or restricted expression pattern. The results of this study provide the most comprehensive listing of two families of snoRNA genes in the human genome till date.     

The snoGloBe interaction predictor reveals a broad spectrum of C/D snoRNA RNA targets

Box C/D small nucleolar RNAs (snoRNAs) are a conserved class of RNA known for their role in guiding ribosomal RNA 2′-O-ribose methylation. Recently, C/D snoRNAs were also implicated in regulating the expression of non-ribosomal genes through different modes of binding. Large scale RNA–RNA interaction datasets detect many snoRNAs binding messenger RNA, but are limited by specific experimental conditions. To enable a more comprehensive study of C/D snoRNA interactions, we created snoGloBe, a human C/D snoRNA interaction predictor based on a gradient boosting classifier. SnoGloBe considers the target type, position and sequence of the interactions, enabling it to outperform existing predictors. Interestingly, for specific snoRNAs, snoGloBe identifies strong enrichment of interactions near gene expression regulatory elements including splice sites. Abundance and splicing of predicted targets were altered upon the knockdown of their associated snoRNA. Strikingly, the predicted snoRNA interactions often overlap with the binding sites of functionally related RNA binding proteins, reinforcing their role in gene expression regulation. SnoGloBe is also an excellent tool for discovering viral RNA targets, as shown by its capacity to identify snoRNAs targeting the heavily methylated SARS-CoV-2 RNA. Overall, snoGloBe is capable of identifying experimentally validated binding sites and predicting novel sites with shared regulatory function.     

A small nucleolar guide RNA functions both in 2'-O-ribose methylation and pseudouridylation of the U5 spliceosomal RNA

In eukaryotes, two distinct classes of small nucleolar RNAs (snoRNAs), namely the fibrillarin-associated box C/D snoRNAs and the Gar1p-associated box H/ACA snoRNAs, direct the site-specific 2'-O-ribose methylation and pseudouridylation of ribosomal RNAs (rRNAs), respectively. We have identified a novel evolutionarily conserved snoRNA, called U85, which possesses the box elements of both classes of snoRNAs and associates with both fibrillarin and Gar1p. In vitro and in vivo pseudouridylation and 2'-O-methylation experiments provide evidence that the U85 snoRNA directs 2'-O-methylation of the C45 and pseudouridylation of the U46 residues in the invariant loop 1 of the human U5 spliceosomal RNA. The U85 is the first example of a snoRNA that directs modification of an RNA polymerase II-transcribed spliceosomal RNA and that functions both in RNA pseudouridylation and 2'-O-methylation.     

Plant U13 orthologues and orphan snoRNAs identified by RNomics of RNA from Arabidopsis nucleoli

Small nucleolar RNAs (snoRNAs) and small Cajal body-specific RNAs (scaRNAs) are non-coding RNAs whose main function in eukaryotes is to guide the modification of nucleotides in ribosomal and spliceosomal small nuclear RNAs, respectively. Full-length sequences of Arabidopsis snoRNAs and scaRNAs have been obtained from cDNA libraries of capped and uncapped small RNAs using RNA from isolated nucleoli from Arabidopsis cell cultures. We have identified 31 novel snoRNA genes (9 box C/D and 22 box H/ACA) and 15 new variants of previously described snoRNAs. Three related capped snoRNAs with a distinct gene organization and structure were identified as orthologues of animal U13snoRNAs. In addition, eight of the novel genes had no complementarity to rRNAs or snRNAs and are therefore putative orphan snoRNAs potentially reflecting wider functions for these RNAs. The nucleolar localization of a number of the snoRNAs and the localization to nuclear bodies of two putative scaRNAs was confirmed by in situ hybridization. The majority of the novel snoRNA genes were found in new gene clusters or as part of previously described clusters. These results expand the repertoire of Arabidopsis snoRNAs to 188 snoRNA genes with 294 gene variants.     

RNomics: an experimental approach that identifies 201 candidates for novel, small, non-messenger RNAs in mouse

In mouse brain cDNA libraries generated from small RNA molecules we have identified a total of 201 different expressed RNA sequences potentially encoding novel small non-messenger RNA species (snmRNAs). Based on sequence and structural motifs, 113 of these RNAs can be assigned to the C/D box or H/ACA box subclass of small nucleolar RNAs (snoRNAs), known as guide RNAs for rRNA. While 30 RNAs represent mouse homologues of previously identified human C/D or H/ACA snoRNAs, 83 correspond to entirely novel snoRNAS: Among these, for the first time, we identified four C/D box snoRNAs and four H/ACA box snoRNAs predicted to direct modifications within U2, U4 or U6 small nuclear RNAs (snRNAs). Furthermore, 25 snoRNAs from either class lacked antisense elements for rRNAs or snRNAS: Therefore, additional snoRNA targets have to be considered. Surprisingly, six C/D box snoRNAs and one H/ACA box snoRNA were expressed exclusively in brain. Of the 88 RNAs not belonging to either snoRNA subclass, at least 26 are probably derived from truncated heterogeneous nuclear RNAs (hnRNAs) or mRNAS: Short interspersed repetitive elements (SINEs) are located on five RNA sequences and may represent rare examples of transcribed SINES: The remaining RNA species could not as yet be assigned either to any snmRNA class or to a part of a larger hnRNA/mRNA. It is likely that at least some of the latter will represent novel, unclassified snmRNAS:     

Molecular evolution of the HBII-52 snoRNA cluster

HBII-52 is a human brain-specific C/D box snoRNA that potentially regulates the editing and/or alternative splicing of the serotonin receptor. Forty-two nearly identical copies of the HBII-52 gene are located immediately downstream of the SNRPN protein-coding gene in an imprinted locus associated with Prader-Willi syndrome. Other eutherian mammals, with genomic assemblies covering the corresponding locus, also have multiple orthologous copies of HBII-52. The SNRPB gene, which is known to have given rise to SNRPN through gene duplication, expresses a C/D box snoRNA, SNORD119, from its fifth intron. Here we show that, despite the fact that they lie in different positions relative to the orthologous SNRPB/SNRPN coding sequences, there are significant sequence similarities between SNORD119 and HBII-52, including the antisense element and the stem-forming regions. By analysing these snoRNAs in marsupial and eutherian mammal genomes, we reconstruct the likely evolutionary history of the HBII-52 cluster and SNORD119 and suggest that they have evolved from a common ancestor.     

snoRNA的结构与功能

核仁小分子RNA（snoRNA）是一类广泛分布于真核生物细胞核仁的小分子非编码RNA,具有保守的结构元件,并以此划分为3大类：boxC／DsnoRNA、boxH／ACAsnoRNA和MRPRNA。其中boxC／D和boxH／ACA是已知snoRNA的主要类型,以碱基配对的方式分别指导着核糖体RNA的甲基化和假尿嘧啶化修饰。研究发现,snoRNA除了在核糖体RNA的生物合成中发挥作用之外,还能够指导snRNA、tRNA和mRNA的转录后修饰。此外,还有相当数量的snoRNA功能不明,被称为孤儿sn0RNA（orphansnoRNA）。在哺乳动物的孤儿snoRNA中,印迹snoRNA（imprintedsnoRNA）是最为特殊的一群,由基因组印迹区编码,具有明显的组织表达特异性。原核生物古细菌中类snoRNA的鉴定表明这些非编码RNA家族成员的古老起源;而哺乳动物中大量的snoRNA反转座子的存在更为人们探索snoRNA在基因组中扩增和功能进化提供了新的思路。     

Lack of Pwcr1/MBII-85 snoRNA is critical for neonatal lethality in Prader–Willi syndrome mouse models

Prader–Willi syndrome (PWS) is a neurobehavioral disorder caused by the lack of paternal expression of imprinted genes in the human chromosome region 15q11–13. Recent studies of rare human translocation patients narrowed the PWS critical genes to a 121-kb region containing PWCR1/HBII-85 and HBII-438 snoRNA genes. The existing mouse models of PWS that lack the expression of multiple genes, including Snrpn, Ube3a, and many intronic snoRNA genes, are characterized by 80%–100% neonatal lethality. To define the candidate region for PWS-like phenotypes in mice, we analyzed the expression of several genetic elements in mice carrying the large radiation-induced p^30PUb deletion that includes the p locus. Mice having inherited this deletion from either parent develop normally into adulthood. By Northern blot and RT-PCR assays of brain tissue, we found that Pwcr1/MBII-85 snoRNAs are expressed normally, while MBII-52 snoRNAs are not expressed when the deletion is paternally inherited. Mapping of the distal deletion breakpoint indicated that the p^30PUb deletion includes the entire MBII-52 snoRNA gene cluster and three previously unmapped EST sequences. The lack of expression of these elements in mice with a paternal p^30PUb deletion indicates that they are not critical for the neonatal lethality observed in PWS mouse models. In addition, we identified MBII-436, the mouse homolog of the HBII-436 snoRNA, confirmed its imprinting status, and mapped it outside of the p^30PUb deletion. Taking together all available data, we conclude that the lack of Pwcr1/MBII-85 snoRNA expression is the most likely cause for the neonatal lethality in PWS model mice.     

Plant snoRNA database

The Plant snoRNA database (http://www.scri.sari.ac.uk/plant_snoRNA/) provides information on small nucleolar RNAs from Arabidopsis and eighteen other plant species. Information includes sequences, expression data, methylation and pseudouridylation target modification sites, initial gene organization (polycistronic, single gene and intronic) and the number of gene variants. The Arabidopsis information is divided into box C/D and box H/ACA snoRNAs, and within each of these groups, by target sites in rRNA, snRNA or unknown. Alignments of orthologous genes and gene variants from different plant species are available for many snoRNA genes. Plant snoRNA genes have been given a standard nomenclature, designed wherever possible, to provide a consistent identity with yeast and human orthologues.     

Direct cloning of double-stranded RNAs from RNase protection analysis reveals processing patterns of C/D box snoRNAs and provides evidence for widespread antisense transcript expression

We describe a new method that allows cloning of double-stranded RNAs (dsRNAs) that are generated in RNase protection experiments. We demonstrate that the mouse C/D box snoRNA MBII-85 (SNORD116) is processed into at least five shorter RNAs using processing sites near known functional elements of C/D box snoRNAs. Surprisingly, the majority of cloned RNAs from RNase protection experiments were derived from endogenous cellular RNA, indicating widespread antisense expression. The cloned dsRNAs could be mapped to genome areas that show RNA expression on both DNA strands and partially overlapped with experimentally determined argonaute-binding sites. The data suggest a conserved processing pattern for some C/D box snoRNAs and abundant expression of longer, non-coding RNAs in the cell that can potentially form dsRNAs.     

The small nucleolar ribonucleoprotein (snoRNP) database

Small nucleolar ribonucleoproteins (snoRNPs) are widely studied and characterized as guide RNAs for sequence-specific 2′-O-ribose methylation and psuedouridylation of ribosomal RNAs. In addition, snoRNAs have also been shown to interact with some tRNAs and direct alternative splicing in mRNA biogenesis. Recent advances in bioinformatics have resulted in new algorithms able to rapidly identify noncoding RNAs generally and snoRNAs specifically in genomic and metagenomic sequences, resulting in a rapid increase in the number and diversity of identified snoRNA sequences. The snoRNP database is a web-based collection of snoRNA and snoRNA-associated protein sequences from a wide range of species. The database currently contains 8994 snoRNA sequences from Bacteria, Archaea, and Eukaryotes and 589 snoRNA-associated protein sequences. The snoRNP database can be found at: http://evolveathome.com/snoRNA/snoRNA.php.     

Efficient and specific knockdown of small non-coding RNAs in mammalian cells and in mice

Hundreds of small nuclear non-coding RNAs, including small nucleolar RNAs (snoRNAs), have been identified in different organisms, with important implications in regulating gene expression and in human diseases. However, functionalizing these nuclear RNAs in mammalian cells remains challenging, due to methodological difficulties in depleting these RNAs, especially snoRNAs. Here we report a convenient and efficient approach to deplete snoRNA, small Cajal body RNA (scaRNA) and small nuclear RNA in human and mouse cells by conventional transfection of chemically modified antisense oligonucleotides (ASOs) that promote RNaseH-mediated cleavage of target RNAs. The levels of all seven tested snoRNA/scaRNAs and four snRNAs were reduced by 80-95%, accompanied by impaired endogenous functions of the target RNAs. ASO-targeting is highly specific, without affecting expression of the host genes where snoRNAs are embedded in the introns, nor affecting the levels of snoRNA isoforms with high sequence similarities. At least five snoRNAs could be depleted simultaneously. Importantly, snoRNAs could be dramatically depleted in mice by systematic administration of the ASOs. Together, our findings provide a convenient and efficient approach to characterize nuclear non-coding RNAs in mammalian cells, and to develop antisense drugs against disease-causing non-coding RNAs.