Mining small RNA sequencing data: a new approach to identify small nucleolar RNAs in Arabidopsis

Small nucleolar RNAs (snoRNAs) are noncoding RNAs that direct 2′-O-methylation or pseudouridylation on ribosomal RNAs or spliceosomal small nuclear RNAs. These modifications are needed to modulate the activity of ribosomes and spliceosomes. A comprehensive repertoire of snoRNAs is needed to expand the knowledge of these modifications. The sequences corresponding to snoRNAs in 18–26-nt small RNA sequencing data have been rarely explored and remain as a hidden treasure for snoRNA annotation. Here, we showed the enrichment of small RNAs at Arabidopsis snoRNA termini and developed a computational approach to identify snoRNAs on the basis of this characteristic. The approach successfully uncovered the full-length sequences of 144 known Arabidopsis snoRNA genes, including some snoRNAs with improved 5′- or 3′-end annotation. In addition, we identified 27 and 17 candidates for novel box C/D and box H/ACA snoRNAs, respectively. Northern blot analysis and sequencing data from parallel analysis of RNA ends confirmed the expression and the termini of the newly predicted snoRNAs. Our study especially expanded on the current knowledge of box H/ACA snoRNAs and snoRNA species targeting snRNAs. In this study, we demonstrated that the use of small RNA sequencing data can increase the complexity and the accuracy of snoRNA annotation.     

The high diversity of snoRNAs in plants: identification and comparative study of 120 snoRNA genes from Oryza sativa

Using a powerful computer-assisted analysis strategy, a large-scale search of small nucleolar RNA (snoRNA) genes in the recently released draft sequence of the rice genome was carried out. This analysis identified 120 different box C/D snoRNA genes with a total of 346 gene variants, which were predicted to guide 135 2′-O-ribose methylation sites in rice rRNAs. Though not exhaustive, this analysis has revealed that rice has the highest number of known box C/D snoRNAs among eukaryotes. Interestingly, although many snoRNA genes are conserved between rice and Arabidopsis, almost half of the identified snoRNA genes are rice specific, which may highlight further the differences in rRNA methylation patterns between monocotyledons and dicotyledons. In addition to 76 singletons, 70 clusters involving 270 snoRNA genes were also found in rice. The large number of the novel snoRNA polycistrons found in the introns of rice protein-coding genes is in contrast to the one-snoRNA-per-intron organization of vertebrates and yeast, and of Arabidopsis in which only a few intronic snoRNA gene clusters were identified. Furthermore, due to a high degree of gene duplication, rice snoRNA genes are clearly redundant and exhibit great sequence variation among isoforms, allowing generation of new snoRNAs for selection. Thus, the large snoRNA gene family in plants can serve as an excellent model for a rapid and functional evolution.     

Identification of a new ribose methylation in the 18S rRNA of S. cerevisiae

Methylation of ribose sugars at the 2′-OH group is one of the major chemical modifications in rRNA, and is catalyzed by snoRNA directed C/D box snoRNPs. Previous biochemical and computational analyses of the C/D box snoRNAs have identified and mapped a large number of 2′-OH ribose methylations in rRNAs. In the present study, we systematically analyzed ribose methylations of 18S rRNA in Saccharomyces cerevisiae, using mung bean nuclease protection assay and RP-HPLC. Unexpectedly, we identified a hitherto unknown ribose methylation at position G562 in the helix 18 of 5′ central domain of yeast 18S rRNA. Furthermore, we identified snR40 as being responsible to guide snoRNP complex to catalyze G562 ribose methylation, which makes it only second snoRNA known so far to target three ribose methylation sites: Gm562, Gm1271 in 18S rRNA, and Um898 in 25S rRNA. Our sequence and mutational analysis of snR40 revealed that snR40 uses the same D′ box and methylation guide sequence for both Gm562 and Gm1271 methylation. With the identification of Gm562 and its corresponding snoRNA, complete set of ribose methylations of 18S rRNA and their corresponding snoRNAs have finally been established opening great prospects to understand the physiological function of these modifications.     

Characterisation of the U83 and U84 small nucleolar RNAs: two novel 2'-O-ribose methylation guide RNAs that lack complementarities to ribosomal RNAs

In eukaryotic cells, the site-specific 2′-O-ribose methy-lation of ribosomal RNAs (rRNAs) and the U6 spliceosomal small nuclear RNA (snRNA) is directed by small nucleolar RNAs (snoRNAs). The C and D box-containing 2′-O-methylation guide snoRNAs select the correct substrate nucleotide through formation of a long 10–21 bp interaction with the target rRNA and U6 snRNA sequences. Here, we report on the characterisation of two novel mammalian C/D box snoRNAs, called U83 and U84, that contain all the elements that are essential for accumulation and function of 2′-O-methylation guide snoRNAs. However, in contrast to all of the known 2′-O-methylation guide RNAs, the human, mouse and pig U83 and U84 snoRNAs feature no antisense elements complementary to rRNA or U6 snRNA sequences. The human U83 and U84 snoRNAs are not associated with maturing nucleolar pre-ribosomal particles, suggesting that they do not function in rRNA biogenesis. Since artificial substrate RNAs complementary to the evolutionarily conserved putative substrate recognition motifs of the U83 and U84 snoRNAs were correctly 2′-O-methy-lated in the nucleolus of mouse cells, we suggest that the new snoRNAs act as 2′-O-methylation guides for cellular RNAs other then rRNAs and the U6 snRNA.     

Elucidating the role of C/D snoRNA in rRNA processing and modification in Trypanosoma brucei

Most eukaryotic C/D small nucleolar RNAs (snoRNAs) guide 2′-O methylation (Nm) on rRNA and are also involved in rRNA processing. The four core proteins that bind C/D snoRNA in Trypanosoma brucei are fibrillarin (NOP1), NOP56, NOP58, and SNU13. Silencing of NOP1 by RNA interference identified rRNA-processing and modification defects that caused lethality. Systematic mapping of 2′-O-methyls on rRNA revealed the existence of hypermethylation at certain positions of the rRNA in the bloodstream form of the parasites, suggesting that this modification may assist the parasites in coping with the major temperature changes during cycling between their insect and mammalian hosts. The rRNA-processing defects of NOP1-depleted cells suggest the involvement of C/D snoRNA in trypanosome-specific rRNA-processing events to generate the small rRNA fragments. MRP RNA, which is involved in rRNA processing, was identified in this study in one of the snoRNA gene clusters, suggesting that trypanosomes utilize a combination of unique C/D snoRNAs and conserved snoRNAs for rRNA processing.     

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.     

A novel snoRNA can direct site-specific 2'-O-ribose methylation of snRNAs in Oryza sativa

Small nucleolar RNAs (snoRNAs) are a kind of noncoding RNAs, and the vast majority of snoRNAs are involved in site-specific modifications of rRNAs. A novel box C/D snoRNA called snoR124 was found inOryza sativa, and it can direct 2'-O-ribose methylation of spliceosomal small nuclear RNAs (snRNAs). The snoRNA has two antisense elements, and the results of primer extensions at different dNTP concentrations provide evidence that snoR124 guide 2'-O-methylations of the C76 residue in the U4 snRNA and the T91 residue in the U5 snRNA. In addition, this snoRNA is located in a snoRNA gene cluster with another 7 snoRNAs which are identified to direct ribose methylations in rRNAs. This is consistent with the opinion that the snoRNA gene organization in plant is mainly gene cluster. The snoR124 is the first example of a snoRNA that directs modifications of RNAs other than rRNAs in plant; it will avail to get more insights into the function of snoRNAs in plant.     

The box C/D motif directs snoRNA 5'-cap hypermethylation

The 5′-cap structure of most spliceosomal small nuclear RNAs (snRNAs) and certain small nucleolar RNAs (snoRNAs) undergoes hypermethylation from a 7-methylguanosine to a 2,2,7-trimethylguanosine structure. 5′-Cap hypermethylation of snRNAs is dependent upon a conserved sequence element known as the Sm site common to most snRNAs. Here we have performed a mutational analysis of U3 and U14 to determine the cis-acting sequences required for 5′-cap hypermethylation of Box C/D snoRNAs. We have found that both the conserved sequence elements Box C (termed C′ in U3) and Box D are necessary for cap hypermethylation. Furthermore, the terminal stem structure that is formed by sequences that flank Box C (C′ in U3) and Box D is also required. However, mutation of other conserved sequences has no effect on hypermethylation of the cap. Finally, the analysis of fragments of U3 and U14 RNAs indicates that the Box C/D motif, including Box C (C′ in U3), Box D and the terminal stem, is capable of directing cap hypermethylation. Thus, the Box C/D motif, which is important for snoRNA processing, stability, nuclear retention, protein binding, nucleolar localization and function, is also necessary and sufficient for cap hypermethylation of these 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.     

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.     

Sequence and structural elements of methylation guide snoRNAs essential for site-specific ribose methylation of pre-rRNA.

Site-specific 2'-O-ribose methylation of eukaryotic rRNAs is guided by small nucleolar RNAs (snoRNAs). The methylation guide snoRNAs carry long perfect complementaries to rRNAs. These antisense elements are located either in the 5' half or in the 3' end region of the snoRNA, and are followed by the conserved D' or D box motifs, respectively. An uninterrupted helix formed between the rRNA and the antisense element of the snoRNA, in conjunction with the adjacent D' or D box, constitute the recognition signal for the putative methyltransferase. Here, we have identified an additional essential box element common to methylation guide snoRNAs, termed the C' box. We show that the C' box functions in concert with the D' box and plays a crucial role in the methyltransfer reaction directed by the upstream antisense element and the D' box. We also show that an internal fragment of U24 methylation guide snoRNA, encompassing the upstream antisense element and the D' and C' box motifs, can support the site-specific methylation of rRNA. This strongly suggests that the C box of methylation guide snoRNAs plays an essential role in the methyltransfer reaction guided by the 3'-terminal antisense element and the D box of the snoRNA.     

Genomewide analysis of box C/D and box H/ACA snoRNAs in Chlamydomonas reinhardtii reveals an extensive organization into intronic gene clusters

Chlamydomonas reinhardtii is a unicellular green alga, the lineage of which diverged from that of land plants >1 billion years ago. Using the powerful small nucleolar RNA (snoRNA) mining platform to screen the C. reinhardtii genome, we identified 322 snoRNA genes grouped into 118 families. The 74 box C/D families can potentially guide methylation at 96 sites of ribosomal RNAs (rRNAs) and snRNAs, and the 44 box H/ACA families can potentially guide pseudouridylation at 62 sites. Remarkably, 242 of the snoRNA genes are arranged into 76 clusters, of which 77% consist of homologous genes produced by small local tandem duplications. At least 70 snoRNA gene clusters are found within introns of protein-coding genes. Although not exhaustive, this analysis reveals that C. reinhardtii has the highest number of intronic snoRNA gene clusters among eukaryotes. The prevalence of intronic snoRNA gene clusters in C. reinhardtii is similar to that of rice but in contrast with the one-snoRNA-per-intron organization of vertebrates and fungi and with that of Arabidopsis thaliana in which only a few intronic snoRNA gene clusters were identified. This analysis of C. reinhardtii snoRNA gene organization shows the functional importance of introns in a single-celled organism and provides evolutionary insight into the origin of intron-encoded RNAs in the plant lineage.     

Fibrillarin-associated box C/D small nucleolar RNAs in Trypanosoma brucei. Sequence conservation and implications for 2'-O-ribose methylation of rRNA

We report the identification of 17 box C/D fibrillarin-associated small nucleolar RNAs (snoRNAs) from the ancient eukaryote, Trypanosoma brucei. To systematically isolate and characterize these snoRNAs, the T. brucei cDNA for the box C/D snoRNA common protein, fibrillarin, was cloned and polyclonal antibodies to the recombinant fibrillarin protein were generated in rabbits. Immunoprecipitations from T. brucei extracts with the anti-fibrillarin antibodies indicated that this trypanosomatid has at least 30 fibrillarin-associated snoRNAs. We have sequenced seventeen of them and designated them TBR for T. brucei RNA 1-17. All of them bear conserved box C, D, C', and D' elements, a hallmark of fibrillarin-associated snoRNAs in eukaryotes. Fourteen of them are novel T. brucei snoRNAs. Fifteen bear potential guide regions to mature rRNAs suggesting that they are involved in 2'-O-ribose methylation. Indeed, eight ribose methylations have been mapped in the rRNA at sites predicted by the snoRNA sequences. Comparative genomics indicates that six of the seventeen are the first trypanosome homologs of known yeast and vertebrate methylation guide snoRNAs. Our results indicate that T. brucei has many fibrillarin-associated box C/D snoRNAs with roles in 2'-O-ribose methylation of rRNA and that the mechanism for targeting the nucleotide to be methylated at the fifth nucleotide upstream of box D or D' originated in early eukaryotes.     

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.