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.     

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.     

Computational prediction of Caenorhabditis box H/ACA snoRNAs using genomic properties of their host genes

Identification of small nucleolar RNAs (snoRNAs) in genomic sequences has been challenging due to the relative paucity of sequence features. Many current prediction algorithms rely on detection of snoRNA motifs complementary to target sites in snRNAs and rRNAs. However, recent discovery of snoRNAs without apparent targets requires development of alternative prediction methods. We present an approach that combines rule-based filters and a Bayesian Classifier to identify a class of snoRNAs (H/ACA) without requiring target sequence information. It takes advantage of unique attributes of their genomic organization and improved species-specific motif characterization to predict snoRNAs that may otherwise be difficult to discover. Searches in the genomes of Caenorhabditis elegans and the closely related Caenorhabditis briggsae suggest that our method performs well compared to recent benchmark algorithms. Our results illustrate the benefits of training gene discovery engines on features restricted to particular phylogenetic groups and the utility of incorporating diverse data types in gene prediction.     

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在基因组中扩增和功能进化提供了新的思路。     

The complete set of H/ACA snoRNAs that guide rRNA pseudouridylations in Saccharomyces cerevisiae

Conversion of uridines into pseudouridines (Psis) is the most frequent base modification in ribosomal RNAs (rRNAs). In eukaryotes, the pseudouridylation sites are specified by base-pairing with specific target sequences within H/ACA small nucleolar RNAs (snoRNAs). The yeast rRNAs harbor 44 Psis, but, when this work began, 15 Psis had completely unknown guide snoRNAs. This suggested that many snoRNAs remained to be discovered. To address this problem and further complete the snoRNA assignment to Psi sites, we identified the complete set of RNAs associated with the H/ACA snoRNP specific proteins Gar1p and Nhp2p by coupling TAP-tag purifications with genomic DNA microarrays experiments. Surprisingly, while we identified all the previously known H/ACA snoRNAs, we selected only three new snoRNAs. This suggested that most of the missing Psi guides were present in previously known snoRNAs but had been overlooked. We confirmed this hypothesis by systematically investigating the role of previously known, as well as of the newly identified snoRNAs, in specifying rRNA Psi sites and found all but one missing guide RNAs. During the completion of this work, another study, based on bioinformatic predictions, also reported the identification of most missing guide RNAs. Altogether, all Psi guides are now identified and we can tell that, in budding yeast, the 44 Psis are guided by 28 snoRNAs. Finally, aside from snR30, an atypical small RNA of heterogeneous length and at least one mRNA, all Gar1p and Nhp2p associated RNAs characterized by our work turned out to be snoRNAs involved in rRNA Psi specification.     

Identification of 66 box C/D snoRNAs in Arabidopsis thaliana: extensive gene duplications generated multiple isoforms predicting new ribosomal RNA 2'-O-methylation sites

Dozens of box C/D small nucleolar RNAs (snoRNAs) have recently been found in eukaryotes (vertebrates, yeast), ancient eukaryotes (trypanosomes) and archae, that specifically target ribosomal RNA sites for 2'-O-ribose methylation. Although early biochemical data revealed that plant rRNAs are among the most highly ribomethylated in eukaryotes, only a handful of methylation guide snoRNAs have been characterized in this kingdom. We report 66 novel box C/D snoRNAs identified by computational screening of Arabidopsis genomic sequences that are expressed in vivo from either single genes, 17 different clusters or three introns. At the structural level, many box C/D snoRNAs have dual antisense elements often matching rRNA regions close to each other on the rRNA secondary structure, which is reminiscent of their archaeal counterparts. Remarkable specimens are found that display two antisense elements having the potential to form an extended snoRNA-rRNA duplex of 23 to 30 nt, in line with the hypothetical function of box C/D snoRNAs in pre-rRNA folding or chaperoning. In contrast to other species, many Arabidopsis snoRNAs are found in multiple isoforms mainly resulting from two different mechanisms: large chromosomal duplications and small tandem duplications producing polycistronic genes. The discovery of numerous different snoRNAs, some of them arising from common ancestors, provide new insights to understand snoRNAs evolution and the birth of new rRNA methylation sites in plants and other organisms.     

核仁小RNA源性小RNA

最新研究结果表明,一些与RNA介导基因沉默相关的小RNA由核仁小RNA(small nucleolar RNA,snoRNA)加工产生,这种小RNA被称为核仁小RNA源性小RNA(snoRNA derived small RNA,sdRNA)。sdRNA现象分布物种广;涉及的snoRNA种类全,数量多;产生的小RNA分子大小不一、数量、种类多。表明这种小RNA在生物中存在着广泛的普遍性。sdRNA的发现拓展了snoRNA的功能,揭示了snoRNA与RNA介导的基因沉默之间的紧密关系,增强了snoRNA在RNA调控网络中的重要性,并为进一步研究RNA调控网络开启了一扇门。     

Pseudouridine-guide RNAs and other Cbf5p-associated RNAs in Euglena gracilis

In eukaryotes, box H/ACA small nucleolar RNAs (snoRNAs) guide sites of pseudouridine (Psi) formation in rRNA. These snoRNAs reside in RNP complexes containing the putative Psi synthase, Cbf5p. In this study we have identified Cbf5p-associated RNAs in Euglena gracilis, an early diverging eukaryote, by immunoprecipitating Cbf5p-containing complexes from cellular extracts. We characterized one box H/ACA-like RNA which, however, does not appear to guide Psi formation in rRNA. We also identified four single Psi-guide box AGA RNAs. We determined target sites for these putative Psi-guide RNAs and confirmed that the predicted Psi modifications do, in fact, occur at these positions in Euglena rRNA. The Cbf5p-associated snoRNAs appear to be encoded by multicopy genes, some of which are clustered in the genome together with methylation-guide snoRNA genes. These modification-guide snoRNAs and snoRNA genes are the first ones to be reported in euglenid protists, the evolutionary sister group to the kinetoplastid protozoa. Unexpectedly, we also found and have partially characterized a selenocysteine tRNA homolog in the anti-Cbf5p-immunoprecipitated sample.     

Prediction of structured non-coding RNAs in the genomes of the nematodes Caenorhabditis elegans and Caenorhabditis briggsae

We present a survey for non-coding RNAs and other structured RNA motifs in the genomes of Caenorhabditis elegans and Caenorhabditis briggsae using the RNAz program. This approach explicitly evaluates comparative sequence information to detect stabilizing selection acting on RNA secondary structure. We detect 3,672 structured RNA motifs, of which only 678 are known non-translated RNAs (ncRNAs) or clear homologs of known C. elegans ncRNAs. Most of these signals are located in introns or at a distance from known protein-coding genes. With an estimated false positive rate of about 50% and a sensitivity on the order of 50%, we estimate that the nematode genomes contain between 3,000 and 4,000 RNAs with evolutionary conserved secondary structures. Only a small fraction of these belongs to the known RNA classes, including tRNAs, snoRNAs, snRNAs, or microRNAs. A relatively small class of ncRNA candidates is associated with previously observed RNA-specific upstream elements.     

The expanding snoRNA world

In eukaryotes, the site-specific formation of the two prevalent types of rRNA modified nucleotides, 2'-O-methylated nucleotides and pseudouridines, is directed by two large families of snoRNAs. These are termed box C/D and H/ACA snoRNAs, respectively, and exert their function through the formation of a canonical guide RNA duplex at the modification site. In each family, one snoRNA acts as a guide for one, or at most two modifications, through a single, or a pair of appropriate antisense elements. The two guide families now appear much larger than anticipated and their role not restricted to ribosome synthesis only. This is reflected by the recent detection of guides that can target other cellular RNAs, including snRNAs, tRNAs and possibly even mRNAs, and by the identification of scores of tissue-specific specimens in mammals. Recent characterization of homologs of eukaryotic modification guide snoRNAs in Archaea reveals the ancient origin of these non-coding RNA families and offers new perspectives as to their range of function.     

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.     

Mammalian small nucleolar RNAs are mobile genetic elements

Small nucleolar RNAs (snoRNAs) of the H/ACA box and C/D box categories guide the pseudouridylation and the 2'-O-ribose methylation of ribosomal RNAs by forming short duplexes with their target. Similarly, small Cajal body-specific RNAs (scaRNAs) guide modifications of spliceosomal RNAs. The vast majority of vertebrate sno/scaRNAs are located in introns of genes transcribed by RNA polymerase II and processed by exonucleolytic trimming after splicing. A bioinformatic search for orthologues of human sno/scaRNAs in sequenced mammalian genomes reveals the presence of species- or lineage-specific sno/scaRNA retroposons (sno/scaRTs) characterized by an A-rich tail and an approximately 14-bp target site duplication that corresponds to their insertion site, as determined by interspecific genomic alignments. Three classes of snoRTs are defined based on the extent of intron and exon sequences from the snoRNA parental host gene they contain. SnoRTs frequently insert in gene introns in the sense orientation at genomic hot spots shared with other genetic mobile elements. Previously characterized human snoRNAs are encoded in retroposons whose parental copies can be identified by phylogenic analysis, showing that snoRTs can be faithfully processed. These results identify snoRNAs as a new family of mobile genetic elements. The insertion of new snoRNA copies might constitute a safeguard mechanism by which the biological activity of snoRNAs is maintained in spite of the risk of mutations in the parental copy. I furthermore propose that retroposition followed by genetic drift is a mechanism that increased snoRNA diversity during vertebrate evolution to eventually acquire new RNA-modification functions.