A search for H/ACA snoRNAs in yeast using MFE secondary structure prediction

MOTIVATION: Noncoding RNA genes produce functional RNA molecules rather than coding for proteins. One such family is the H/ACA snoRNAs. Unlike the related C/D snoRNAs these have resisted automated detection to date. RESULTS: We develop an algorithm to screen the yeast genome for novel H/ACA snoRNAs. To achieve this, we introduce some new methods for facilitating the search for noncoding RNAs in genomic sequences which are based on properties of predicted minimum free-energy (MFE) secondary structures. The algorithm has been implemented and can be generalized to enable screening of other eukaryote genomes. We find that use of primary sequence alone is insufficient for identifying novel H/ACA snoRNAs. Only the use of secondary structure filters reduces the number of candidates to a manageable size. From genomic context, we identify three strong H/ACA snoRNA candidates. These together with a further 47 candidates obtained by our analysis are being experimentally screened.     

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:     

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.     

Processing of the intron-encoded U16 and U18 snoRNAs: the conserved C and D boxes control both the processing reaction and the stability of the mature snoRNA.

A novel class of small nucleolar RNAs (snoRNAs), encoded in introns of protein coding genes and originating from processing of their precursor molecules, has recently been described. The L1 ribosomal protein (r-protein) gene of Xenopus laevis and its human homologue contain two snoRNAs, U16 and U18. It has been shown that these snoRNAs are excised from their intron precursors by endonucleolytic cleavage and that their processing is alternative to splicing. Two sequences, internal to the snoRNA coding region, have been identified as indispensable for processing the conserved boxes C and D. Competition experiments have shown that these sequences interact with diffusible factors which can bind both the pre-mRNA and the mature U16 snoRNA. Fibrillarin, which is known to associate with complexes formed on C and D boxes of other snoRNAs, is found in association with mature U16 RNA, as well as with its precursor molecules. This fact suggests that the complex formed on the pre-mRNA remains bound to U16 throughout all the processing steps. We also show that the complex formed on the C and D boxes is necessary to stabilize mature snoRNA.     

The genes for small nucleolar RNAs in Trypanosoma brucei are organized in clusters and are transcribed as a polycistronic RNA

Because the organization of snoRNA genes in vertebrates, plants and yeast is diverse, we investigated the organization of snoRNA genes in a distantly related organism, Trypanosoma brucei. We have characterized the second example of a snoRNA gene cluster that is tandemly repeated in the T.brucei genome. The genes encoding the box C/D snoRNAs TBR12, TBR6, TBR4 and TBR2 make up the cluster. In a genomic organization unique to trypanosomes, there are at least four clusters of these four snoRNA genes tandemly repeated in the T.brucei genome. We show for the first time that the genes encoding snoRNAs in both this cluster and the SLA cluster are transcribed in an unusual way as a polycistronic RNA.     

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.     

AtNUFIP, an essential protein for plant development, reveals the impact of snoRNA gene organisation on the assembly of snoRNPs and rRNA methylation in Arabidopsis thaliana

In all eukaryotes, C/D small nucleolar ribonucleoproteins (C/D snoRNPs) are essential for methylation and processing of ribosomal RNAs. They consist of a box C/D small nucleolar RNA (C/D snoRNA) associated with four highly conserved nucleolar proteins. Recent data in HeLa cells and yeast have revealed that assembly of these snoRNPs is directed by NUFIP protein and other auxiliary factors. Nevertheless, the precise function and biological importance of NUFIP and the other assembly factors remains unknown. In plants, few studies have focused on RNA methylation and snoRNP biogenesis. Here, we identify and characterise the AtNUFIP gene that directs assembly of C/D snoRNP. To elucidate the function of AtNUFIP in planta, we characterized atnufip mutants. These mutants are viable but have severe developmental phenotypes. Northern blot analysis of snoRNA accumulation in atnufip mutants revealed a specific degradation of C/D snoRNAs and this situation is correlated with a reduction in rRNA methylation. Remarkably, the impact of AtNUFIP depends on the structure of snoRNA genes: it is essential for the accumulation of those C/D snoRNAs encoded by polycistronic genes, but not by monocistronic or tsnoRNA genes. We propose that AtNUFIP controls the kinetics of C/D snoRNP assembly on nascent precursors to overcome snoRNA degradation of aberrant RNPs. Finally, we show that AtNUFIP has broader RNP targets, controlling the accumulation of scaRNAs that direct methylation of spliceosomal snRNA in Cajal bodies.     

RNomics in Drosophila melanogaster: identification of 66 candidates for novel non-messenger RNAs

By generating a specialised cDNA library from four different developmental stages of Drosophila melanogaster, we have identified 66 candidates for small non-messenger RNAs (snmRNAs) and have confirmed their expression by northern blot analysis. Thirteen of them were expressed at certain stages of D.melanogaster development, only. Thirty-five species belong to the class of small nucleolar RNAs (snoRNAs), divided into 15 members from the C/D subclass and 20 members from the H/ACA subclass, which mostly guide 2'-O-methylation and pseudouridylation, respectively, of rRNA and snRNAs. These also include two outstanding C/D snoRNAs, U3 and U14, both functioning as pre-rRNA chaperones. Surprisingly, the sequence of the Drosophila U14 snoRNA reflects a major change of function of this snoRNA in Diptera relative to yeast and vertebrates. Among the 22 snmRNAs lacking known sequence and structure motifs, five were located in intergenic regions, two in introns, five in untranslated regions of mRNAs, eight were derived from open reading frames, and two were transcribed opposite to an intron. Interestingly, detection of two RNA species from this group implies that certain snmRNA species are processed from alternatively spliced pre-mRNAs. Surprisingly, a few snmRNA sequences could not be found on the published D.melanogaster genome, which might suggest that more snmRNA genes (as well as mRNAs) are hidden in unsequenced regions of the genome.     

Processing of intron-encoded box C/D small nucleolar RNAs lacking a 5',3'-terminal stem structure

The C and D box-containing (box C/D) small nucleolar RNAs (snoRNAs) function in the nucleolytic processing and 2'-O-methylation of precursor rRNA. In vertebrates, most box C/D snoRNAs are processed from debranched pre-mRNA introns by exonucleolytic activities. Elements directing accurate snoRNA excision are located within the snoRNA itself; they comprise the conserved C and D boxes and an adjoining 5',3'-terminal stem. Although the terminal stem has been demonstrated to be essential for snoRNA accumulation, many snoRNAs lack a terminal helix. To identify the cis-acting elements supporting the accumulation of intron-encoded box C/D snoRNAs devoid of a terminal stem, we have investigated the in vivo processing of the human U46 snoRNA and an artificial snoRNA from the human beta-globin pre-mRNA. We demonstrate that internal and/or external stem structures located within the snoRNA or in the intronic flanking sequences support the accumulation of mammalian box C/D snoRNAs lacking a canonical terminal stem. In the intronic precursor RNA, transiently formed external and/or stable internal base-pairing interactions fold the C and D boxes together and therefore facilitate the binding of snoRNP proteins. Since the external intronic stems are degraded during snoRNA processing, we propose that the C and D boxes alone can provide metabolic stability for the mature snoRNA.     

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.     

Metal resistance in yeast mediated by the expression of a maize 20S proteasome alpha subunit

A novel 72 nt small nucleolar RNA (snoRNA) called U87 was found in rat liver cells. This RNA possesses the features of C/D box snoRNA family: boxes C, D', C', D, and 11 nt antisense element complementary to 28S ribosomal RNA (rRNA). The vast majority of C/D box snoRNAs direct site-specific 2'-O-ribose methylation of rRNAs. U87 RNA is suggested to be involved in 2'-O-methylation of a G(3468) residue in 28S rRNA. U87 RNA was detected in different mammalian species with slight length variability. Rat and mouse U87 RNA gene was characterized. Unlike the majority of C/D box snoRNAs U87 RNA lacks the terminal stem required for snoRNA processing. However, U87 gene is flanked by 7 bp inverted repeats potentially able to form a terminal stem in U87 RNA precursor.     

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.     

Systematic identification and characterization of porcine snoRNAs: structural,functional and developmental insights

Small nucleolar RNAs (snoRNAs) are 50‐ to 300‐nt non‐coding RNAs that are involved in critical cellular events, including rRNA/snRNA modification and splicing, ribosome genesis, telomerase formulation and cell proliferation. The identification of snoRNAs in the pig, which is a widely consumed commercial organism that also has important functions in medicine and biology, will enrich the snoRNA kingdom and provide evolutionary clues about snoRNAs. In this study, we performed a systematic identification of snoRNAs in Sus scrofa and obtained 120 candidate snoRNAs, 65 of which were predicted via sequencing from our constructed cDNA library. The others were obtained by computational screening. The primary structural features examined included the sequence length, GC content, conservation of common box motifs and nucleotide diversity. The results indicate that the primary features of H/ACA box snoRNAs are opposite to those of C/D box snoRNAs. Subsequently, based on chromosomal location and host gene determination, we assigned 91 snoRNAs to nine genome organization modes. Gene duplications and translocations are considered to contribute to the high abundant organization in evolution. Functional information about our novel snoRNAs, such as putative targets, modification sites and guide sequences, was predicted by orthologue alignment. A comparative analysis of predicted targets and possible modified loci on U6 snRNA and 5.8S and 18S rRNAs among five species revealed that targets of snoRNA are conserved among species. Furthermore, we performed a quantitative analysis of six representative snoRNA genes in two pig breeds during different developmental stages. Interestingly, all six snoRNAs from one breed expressed in a similar pattern over the tested time points; however, these same six genes had different expression patterns in the other pig breed. Specifically, expression of all six snoRNAs declined significantly from 65 to 90 days post‐coitus (dpc) and then increased slightly during adulthood in Tongcheng pigs, whereas the expression of the same six genes increased slowly from 65 dpc until adulthood in Landrace pigs. This expression pattern suggests that most housekeeping, non‐coding RNAs from a single pig breed may be similarly expressed during development. Our study adds to the knowledge about the snoRNA family by providing the first genome‐wide study of porcine snoRNAs. The comparative analysis of snoRNAs from different pig breeds gave us evolutionary insight into the function of snoRNAs.     

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.