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

From molecular basics to agronomic benefits: Insights into noncoding RNA-mediated gene regulation in plants

The development of plants is largely dependent on their growth environment. To better adapt to a particular habitat, plants have evolved various subtle regulatory mechanisms for altering gene expression. Non coding RNAs (ncRNAs) constitute a major portion of the transcriptomes of eukaryotes. Various ncRNAs have been recognized as important regulators of the expression of genes involved in essential biological processes throughout the whole life cycles of plants. In this review, we summarize the current understanding of the biogenesis and contributions of small nucle olar RNA (snoRNA)- and regulatory long non coding RNA (lncRNA)-mediated gene regulation in plant development and environmental responses. Many regulatory ncRNAs appear to be associated with increased yield, quality and disease resistance of various species and cultivars. These ncRNAs may potentially be used as genetic resources for improving agronomic traits and for molecular breeding. The challenges in understanding plant ncRNA biology and the possibilities to make better use of these valuable gene resources in the future are discussed in this review.

The small nucleolar RNA SnoR28 regulates plant growth and development by directing rRNA maturation

Small nucleolar RNAs (snoRNAs) are noncoding RNAs (ncRNAs) that guide chemical modifications of structural RNAs, which are essential for ribosome assembly and function in eukaryotes. Although numerous snoRNAs have been identified in plants by high-throughput sequencing, the biological functions of most of these snoRNAs remain unclear. Here, we identified box C/D SnoR28.1s as important regulators of plant growth and development by screening a CRISPR/Cas9-generated ncRNA deletion mutant library in Arabidopsis thaliana. Deletion of the SnoR28.1 locus, which contains a cluster of three genes producing SnoR28.1s, resulted in defects in root and shoot growth. SnoR28.1s guide 2'-O-ribose methylation of 25S rRNA at G2396. SnoR28.1s facilitate proper and efficient pre-rRNA processing, as the SnoR28.1 deletion mutants also showed impaired ribosome assembly and function, which may account for the growth defects. SnoR28 contains a 7-bp antisense box, which is required for 2'-O-ribose methylation of 25S rRNA at G2396, and an 8-bp extra box that is complementary to a nearby rRNA methylation site and is partially responsible for methylation of G2396. Both of these motifs are required for proper and efficient pre-rRNA processing. Finally, we show that SnoR28.1s genetically interact with HIDDEN TREASURE2 and NUCLEOLIN1. Our results advance our understanding of the roles of snoRNAs in Arabidopsis.

Cryo-EM structure and rRNA modification sites of a plant ribosome

Protein synthesis in crop plants contributes to the balance of food and fuel on our planet, which influences human metabolic activity and lifespan. Protein synthesis can be regulated with respect to changing environmental cues via the deposition of chemical modifications into rRNA. Here, we present the structure of a plant ribosome from tomato and a quantitative mass spectrometry analysis of its rRNAs. The study reveals fine features of the ribosomal proteins and 71 plant-specific rRNA modifications, and it re-annotates 30 rRNA residues in the available sequence. At the protein level, isoAsp is found in position 137 of uS11, and a zinc finger previously believed to be universal is missing from eL34, suggesting a lower effect of zinc deficiency on protein synthesis in plants. At the rRNA level, the plant ribosome differs markedly from its human counterpart with respect to the spatial distribution of modifications. Thus, it represents an additional layer of gene expression regulation, highlighting the molecular signature of a plant ribosome. The results provide a reference model of a plant ribosome for structural studies and an accurate marker for molecular ecology.

The Arabidopsis 2'-O-Ribose-Methylation and Pseudouridylation Landscape of rRNA in Comparison to Human and Yeast

Eukaryotic ribosome assembly starts in the nucleolus, where the ribosomal DNA (rDNA) is transcribed into the 35S pre-ribosomal RNA (pre-rRNA). More than two-hundred ribosome biogenesis factors (RBFs) and more than two-hundred small nucleolar RNAs (snoRNA) catalyze the processing, folding and modification of the rRNA in Arabidopsis thaliana. The initial pre-ribosomal 90S complex is formed already during transcription by association of ribosomal proteins (RPs) and RBFs. In addition, small nucleolar ribonucleoprotein particles (snoRNPs) composed of snoRNAs and RBFs catalyze the two major rRNA modification types, 2'-O-ribose-methylation and pseudouridylation. Besides these two modifications, rRNAs can also undergo base methylations and acetylation. However, the latter two modifications have not yet been systematically explored in plants. The snoRNAs of these snoRNPs serve as targeting factors to direct modifications to specific rRNA regions by antisense elements. Today, hundreds of different sites of modifications in the rRNA have been described for eukaryotic ribosomes in general. While our understanding of the general process of ribosome biogenesis has advanced rapidly, the diversities appearing during plant ribosome biogenesis is beginning to emerge. Today, more than two-hundred RBFs were identified by bioinformatics or biochemical approaches, including several plant specific factors. Similarly, more than two hundred snoRNA were predicted based on RNA sequencing experiments. Here, we discuss the predicted and verified rRNA modification sites and the corresponding identified snoRNAs on the example of the model plant Arabidopsis thaliana. Our summary uncovers the plant modification sites in comparison to the human and yeast modification sites.

Profiling of RNA ribose methylation in Arabidopsis thaliana

Eukaryotic rRNAs and snRNAs are decorated with abundant 2′-O-methylated nucleotides (Nm) that are predominantly synthesized by box C/D snoRNA-guided enzymes. In the model plant Arabidopsis thaliana, C/D snoRNAs have been well categorized, but there is a lack of systematic mapping of Nm. Here, we applied RiboMeth-seq to profile Nm in cytoplasmic, chloroplast and mitochondrial rRNAs and snRNAs. We identified 111 Nm in cytoplasmic rRNAs and 19 Nm in snRNAs and assigned guide for majority of the detected sites using an updated snoRNA list. At least four sites are directed by guides with multiple specificities as shown in yeast. We found that C/D snoRNAs frequently form extra pairs with nearby sequences of methylation sites, potentially facilitating the substrate binding. Chloroplast and mitochondrial rRNAs contain five almost identical methylation sites, including two novel sites mediating ribosomal subunit joining. Deletion of FIB1 or FIB2 gene reduced the accumulation of C/D snoRNA and rRNA methylation with FIB1 playing a bigger role in methylation. Our data reveal the comprehensive 2′-O-methylation maps for Arabidopsis rRNAs and snRNAs and would facilitate study of their function and biosynthesis.

Small nucleolar RNAs: continuing identification of novel members and increasing diversity of their molecular mechanisms of action

Identified five decades ago amongst the most abundant cellular RNAs, small nucleolar RNAs (snoRNAs) were initially described as serving as guides for the methylation and pseudouridylation of ribosomal RNA through direct base pairing. In recent years, however, increasingly powerful high-throughput genomic approaches and strategies have led to the discovery of many new members of the family and surprising diversity in snoRNA functionality and mechanisms of action. SnoRNAs are now known to target RNAs of many biotypes for a wider range of modifications, interact with diverse binding partners, compete with other binders for functional interactions, recruit diverse players to targets and affect protein function and accessibility through direct interaction. This mini-review presents the continuing characterization of the snoRNome through the identification of new snoRNA members and the discovery of their mechanisms of action, revealing a highly versatile noncoding family playing central regulatory roles and connecting the main cellular processes.

Mapping rRNA 2'-O-methylations and identification of C/D snoRNAs in Arabidopsis thaliana plants

In all eukaryotic cells, the most abundant modification of ribosomal RNA (rRNA) is methylation at the ribose moiety (2ʹ-O-methylation). Ribose methylation at specific rRNA sites is guided by small nucleolar RNAs (snoRNAs) of C/D-box type (C/D snoRNA) and achieved by the methyltransferase Fibrillarin (FIB). Here we used the Illumina-based RiboMethSeq approach for mapping rRNA 2ʹ-O-methylation sites in A. thaliana Col-0 (WT) plants. This analysis detected novel C/D snoRNA-guided rRNA 2ʹ-O-methylation positions and also some orphan sites without a matching C/D snoRNA. Furthermore, immunoprecipitation of Arabidopsis FIB2 identified and demonstrated expression of C/D snoRNAs corresponding to majority of mapped rRNA sites. On the other hand, we show that disruption of Arabidopsis Nucleolin 1 gene (NUC1), encoding a major nucleolar protein, decreases 2ʹ-O-methylation at specific rRNA sites suggesting functional/structural interconnections of 2ʹ-O-methylation with nucleolus organization and plant development. Finally, based on our findings and existent database sets, we introduce a new nomenclature system for C/D snoRNA in Arabidopsis plants.

The medium-size noncoding RNA transcriptome of Ostreococcus tauri, the smallest living eukaryote, reveals a large family of small nucleolar RNAs displaying multiple genomic expression strategies

The small nucleolar RNAs (snoRNAs), essential for ribosome biogenesis, constitute a major family of medium-size noncoding RNAs (mncRNAs) in all eukaryotes. We present here, for the first time in a marine unicellular alga, the characterization of the snoRNAs family in Ostreococcus tauri, the smallest photosynthetic eukaryote. Using a transcriptomic approach, we identified 131 O. tauri snoRNAs (Ot–snoRNA) distributed in three classes: the C/D snoRNAs, the H/ACA snoRNAs and the MRP RNA. Their genomic organization revealed a unique combination of both the intronic organization of animals and the polycistronic organization of plants. Remarkably, clustered genes produced Ot–snoRNAs with unusual structures never previously described in plants. Their abundances, based on quantification of reads and northern blots, showed extreme differences in Ot–snoRNA accumulation, mainly determined by their differential stability. Most of these Ot–snoRNAs were predicted to target rRNAs or snRNAs. Seventeen others were orphan Ot–snoRNAs that would not target rRNA. These were specific to O. tauri or Mamiellophyceae and could have functions unrelated to ribosome biogenesis. Overall, these data reveal an ‘evolutionary response’ adapted to the extreme compactness of the O. tauri genome that accommodates the essential Ot–snoRNAs, developing multiple strategies to optimize their coordinated expression with a minimal cost on regulatory circuits.

The Existence and Localization of Nuclear snoRNAs in Arabidopsis thaliana Revisited

Ribosome biogenesis is one cell function-defining process. It depends on efficient transcription of rDNAs in the nucleolus as well as on the cytosolic synthesis of ribosomal proteins. For newly transcribed rRNA modification and ribosomal protein assembly, so-called small nucleolar RNAs (snoRNAs) and ribosome biogenesis factors (RBFs) are required. For both, an inventory was established for model systems like yeast and humans. For plants, many assignments are based on predictions. Here, RNA deep sequencing after nuclei enrichment was combined with single molecule species detection by northern blot and in vivo fluorescence in situ hybridization (FISH)-based localization studies. In addition, the occurrence and abundance of selected snoRNAs in different tissues were determined. These approaches confirm the presence of most of the database-deposited snoRNAs in cell cultures, but some of them are localized in the cytosol rather than in the nucleus. Further, for the explored snoRNA examples, differences in their abundance in different tissues were observed, suggesting a tissue-specific function of some snoRNAs. Thus, based on prediction and experimental confirmation, many plant snoRNAs can be proposed, while it cannot be excluded that some of the proposed snoRNAs perform alternative functions than are involved in rRNA modification.

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.

RNA-Seq employing a novel rRNA depletion strategy reveals a rich repertoire of snoRNAs in Euglena gracilis including box C/D and Ψ-guide RNAs targeting the modification of rRNA extremities

Previous mRNA transcriptome studies of Euglena gracilis have shown that this organism possesses a large and diverse complement of protein coding genes; however, the study of non-coding RNA classes has been limited. The natural extensive fragmentation of the E. gracilis large subunit ribosomal RNA presents additional barriers to the identification of non-coding RNAs as size-selected small RNA libraries will be dominated by rRNA sequences. In this study we have developed a strategy to significantly reduce rRNA amplification prior to RNA-Seq analysis thereby producing a ncRNA library allowing for the identification of many new E. gracilis small RNAs. Library analysis reveals 113 unique new small nucleolar (sno) RNAs and a large collection of snoRNA isoforms, as well as the first significant collection of nuclear tRNAs in this organism. A 3' end AGAUGN consensus motif and conserved structural features can now be defined for E. gracilis pseudouridine guide RNAs. snoRNAs of both classes were identified that target modification of the 3' extremities of rRNAs utilizing predicted base-pairing interactions with internally transcribed spacers (ITS), providing insight into the timing of steps in rRNA maturation. Cumulatively, this represents the most comprehensive analysis of small ncRNAs in Euglena gracilis to date.

Phylogenetic distribution of plant snoRNA families

Background

Small nucleolar RNAs (snoRNAs) are one of the most ancient families amongst non-protein-coding RNAs. They are ubiquitous in Archaea and Eukarya but absent in bacteria. Their main function is to target chemical modifications of ribosomal RNAs. They fall into two classes, box C/D snoRNAs and box H/ACA snoRNAs, which are clearly distinguished by conserved sequence motifs and the type of chemical modification that they govern. Similarly to microRNAs, snoRNAs appear in distinct families of homologs that affect homologous targets. In animals, snoRNAs and their evolution have been studied in much detail. In plants, however, their evolution has attracted comparably little attention.

Results

In order to chart the phylogenetic distribution of individual snoRNA families in plants, we applied a sophisticated approach for identifying homologs of known plant snoRNAs across the plant kingdom. In response to the relatively fast evolution of snoRNAs, information on conserved sequence boxes, target sequences, and secondary structure is combined to identify additional snoRNAs. We identified 296 families of snoRNAs in 24 species and traced their evolution throughout the plant kingdom. Many of the plant snoRNA families comprise paralogs. We also found that targets are well-conserved for most snoRNA families.

Conclusions

The sequence conservation of snoRNAs is sufficient to establish homologies between phyla. The degree of this conservation tapers off, however, between land plants and algae. Plant snoRNAs are frequently organized in highly conserved spatial clusters. As a resource for further investigations we provide carefully curated and annotated alignments for each snoRNA family under investigation.

Electronic supplementary material

The online version of this article (doi:10.1186/s12864-016-3301-2) contains supplementary material, which is available to authorized users.

Promoter-based identification of novel non-coding RNAs reveals the presence of dicistronic snoRNA-miRNA genes in Arabidopsis thaliana

Background

In the past few decades, non-coding RNAs (ncRNAs) have emerged as important regulators of gene expression in eukaryotes. Most studies of ncRNAs in plants have focused on the identification of silencing microRNAs (miRNAs) and small interfering RNAs (siRNAs). Another important family of ncRNAs that has been well characterized in plants is the small nucleolar RNAs (snoRNAs) and the related small Cajal body-specific RNAs (scaRNAs). Both target chemical modifications of ribosomal RNAs (rRNAs) and small nuclear RNAs (snRNAs). In plants, the snoRNA genes are organized in clusters, transcribed by RNA Pol II from a common promoter and subsequently processed into mature molecules. The promoter regions of snoRNA polycistronic genes in plants are highly enriched in two conserved cis-regulatory elements (CREs), Telo-box and Site II, which coordinate the expression of snoRNAs and ribosomal protein coding genes throughout the cell cycle.

Results

In order to identify novel ncRNA genes, we have used the snoRNA Telo-box/Site II motifs combination as a functional promoter indicator to screen the Arabidopsis genome. The predictions generated by this process were tested by detailed exploration of available RNA-Seq and expression data sets and experimental validation. As a result, we have identified several snoRNAs, scaRNAs and 'orphan' snoRNAs. We also show evidence for 16 novel ncRNAs that lack similarity to any reported RNA family. Finally, we have identified two dicistronic genes encoding precursors that are processed to mature snoRNA and miRNA molecules. We discuss the evolutionary consequences of this result in the context of a tight link between snoRNAs and miRNAs in eukaryotes.

Conclusions

We present an alternative computational approach for non-coding RNA detection. Instead of depending on sequence or structure similarity in the whole genome screenings, we have explored the properties of promoter regions of well-characterized ncRNAs. Interestingly, besides expected ncRNAs predictions we were also able to recover single precursor arrangement for snoRNA-miRNA. Accompanied by analyses performed on rice sequences, we conclude that such arrangement might have interesting functional and evolutionary consequences and discuss this result in the context of a tight link between snoRNAs and miRNAs in eukaryotes.

Electronic supplementary material

The online version of this article (doi:10.1186/s12864-015-2221-x) contains supplementary material, which is available to authorized users.

plantDARIO: web based quantitative and qualitative analysis of small RNA-seq data in plants

High-throughput sequencing techniques have made it possible to assay an organism's entire repertoire of small non-coding RNAs (ncRNAs) in an efficient and cost-effective manner. The moderate size of small RNA-seq datasets makes it feasible to provide free web services to the research community that provide many basic features of a small RNA-seq analysis, including quality control, read normalization, ncRNA quantification, and the prediction of putative novel ncRNAs. DARIO is one such system that so far has been focussed on animals. Here we introduce an extension of this system to plant short non-coding RNAs (sncRNAs). It includes major modifications to cope with plant-specific sncRNA processing. The current version of plantDARIO covers analyses of mapping files, small RNA-seq quality control, expression analyses of annotated sncRNAs, including the prediction of novel miRNAs and snoRNAs from unknown expressed loci and expression analyses of user-defined loci. At present Arabidopsis thaliana, Beta vulgaris, and Solanum lycopersicum are covered. The web tool links to a plant specific visualization browser to display the read distribution of the analyzed sample. The easy-to-use platform of plantDARIO quantifies RNA expression of annotated sncRNAs from different sncRNA databases together with new sncRNAs, annotated by our group. The plantDARIO website can be accessed at http://plantdario.bioinf.uni-leipzig.de/.

Identification and characterization of the Thermus thermophilus 5-methylcytidine (m5C) methyltransferase modifying 23 S ribosomal RNA (rRNA) base C1942

Methylation of cytidines at carbon-5 is a common posttranscriptional RNA modification encountered across all domains of life. Here, we characterize the modifications of C1942 and C1962 in Thermus thermophilus 23 S rRNA as 5-methylcytidines (m(5)C) and identify the two associated methyltransferases. The methyltransferase modifying C1942, named RlmO, has not been characterized previously. RlmO modifies naked 23 S rRNA, but not the assembled 50 S subunit or 70 S ribosomes. The x-ray crystal structure of this enzyme in complex with the S-adenosyl-l-methionine cofactor at 1.7 Å resolution confirms that RlmO is structurally related to other m(5)C rRNA methyltransferases. Key residues in the active site are located similar to the further distant 5-methyluridine methyltransferase RlmD, suggestive of a similar enzymatic mechanism. RlmO homologues are primarily found in mesophilic bacteria related to T. thermophilus. In accordance, we find that growth of the T. thermophilus strain with an inactivated C1942 methyltransferase gene is not compromised at non-optimal temperatures.

Features of a unique intronless cluster of class I small heat shock protein genes in tandem with box C/D snoRNA genes on chromosome 6 in tomato (Solanum lycopersicum)

Physical clustering of genes has been shown in plants; however, little is known about gene clusters that have different functions, particularly those expressed in the tomato fruit. A class I 17.6 small heat shock protein (Sl17.6 shsp) gene was cloned and used as a probe to screen a tomato (Solanum lycopersicum) genomic library. An 8.3-kb genomic fragment was isolated and its DNA sequence determined. Analysis of the genomic fragment identified intronless open reading frames of three class I shsp genes (Sl17.6, Sl20.0, and Sl20.1), the Sl17.6 gene flanked by Sl20.1 and Sl20.0, with complete 5' and 3' UTRs. Upstream of the Sl20.0 shsp, and within the shsp gene cluster, resides a box C/D snoRNA cluster made of SlsnoR12.1 and SlU24a. Characteristic C and D, and C' and D', boxes are conserved in SlsnoR12.1 and SlU24a while the upstream flanking region of SlsnoR12.1 carries TATA box 1, homol-E and homol-D box-like cis sequences, TM6 promoter, and an uncharacterized tomato EST. Molecular phylogenetic analysis revealed that this particular arrangement of shsps is conserved in tomato genome but is distinct from other species. The intronless genomic sequence is decorated with cis elements previously shown to be responsive to cues from plant hormones, dehydration, cold, heat, and MYC/MYB and WRKY71 transcription factors. Chromosomal mapping localized the tomato genomic sequence on the short arm of chromosome 6 in the introgression line (IL) 6-3. Quantitative polymerase chain reaction analysis of gene cluster members revealed differential expression during ripening of tomato fruit, and relatively different abundances in other plant parts.

Clustered organization, polycistronic transcription, and evolution of modification-guide snoRNA genes in Euglena gracilis

Previous studies have shown that the eukaryotic microbe Euglena gracilis contains an unusually large assortment of small nucleolar RNAs (snoRNAs) and ribosomal RNA (rRNA) modification sites. However, little is known about the evolutionary mechanisms contributing to this situation. In this study, we have examined the organization and evolution of snoRNA genes in Euglena with the additional objective of determining how these properties relate to the rRNA modification pattern in this protist. We have identified and extensively characterized a clustered pattern of genes encoding previously biochemically isolated snoRNA sequences in E. gracilis. We show that polycistronic transcription is a prevalent snoRNA gene expression strategy in this organism. Further, we have identified 121 new snoRNA coding regions through sequence analysis of these clusters. We have identified an E. gracilis U14 snoRNA homolog clustered with modification-guide snoRNA genes. The U14 snoRNAs in other eukaryotic organisms examined to date typically contain both a modification and a processing domain. E. gracilis U14 lacks the modification domain but retains the processing domain. Our analysis of U14 structure and evolution in Euglena and other eukaryotes allows us to propose a model for its evolution and suggest its processing role may be its more important function, explaining its conservation in many eukaryotes. The preponderance of apparent small and larger-scale duplication events in the genomic regions we have characterized in Euglena provides a mechanism for the generation of the unusually diverse collection and abundance of snoRNAs and modified rRNA sites. Our findings provide the framework for more extensive whole genome analysis to elucidate whether these snoRNA gene clusters are spread across multiple chromosomes and/or form dense "arrays" at a limited number of chromosomal loci.

Complete modification maps for the cytosolic small and large subunit rRNAs of Euglena gracilis: functional and evolutionary implications of contrasting patterns between the two rRNA components

In the protist Euglena gracilis, the cytosolic small subunit (SSU) rRNA is a single, covalently continuous species typical of most eukaryotes; in contrast, the large subunit (LSU) rRNA is naturally fragmented, comprising 14 separate RNA molecules instead of the bipartite (28S+5.8S) eukaryotic LSU rRNA typically seen. We present extensively revised secondary structure models of the E. gracilis SSU and LSU rRNAs and have mapped the positions of all of the modified nucleosides in these rRNAs (88 in SSU rRNA and 262 in LSU rRNA, with only 3 LSU rRNA modifications incompletely characterized). The relative proportions of ribose-methylated nucleosides and pseudouridine (∼60% and ∼35%, respectively) are closely similar in the two rRNAs; however, whereas the Euglena SSU rRNA has about the same absolute number of modifications as its human counterpart, the Euglena LSU rRNA has twice as many modifications as the corresponding human LSU rRNA. The increased levels of rRNA fragmentation and modification in E. gracilis LSU rRNA are correlated with a 3-fold increase in the level of mispairing in helical regions compared to the human LSU rRNA. In contrast, no comparable increase in mispairing is seen in helical regions of the SSU rRNA compared to its homologs in other eukaryotes. In view of the reported effects of both ribose-methylated nucleoside and pseudouridine residues on RNA structure, these correlations lead us to suggest that increased modification in the LSU rRNA may play a role in stabilizing a 'looser' structure promoted by elevated helical mispairing and a high degree of fragmentation.

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.

Genome organization and function: a view from yeast and Arabidopsis

Recent progress in understanding higher-order chromatin organization in the nucleus has been considerable. From single gene to chromosome territory, realistic biophysical models can now accurately predict some of the structural feature of cell nuclei. Despite growing evidence of a deterministic nuclear organization, the physiological consequence of spatial genome organization is still unclear. In the simple eukaryotic model, Saccharomyces cerevisiae, clear correlation between gene position and transcription has been established. In this review, we will focus on higher-order chromatin organization in yeast with respect to the nuclear envelope and nucleolus. In Arabidopsis thaliana, a model plant for which we have a complete genome sequence, chromosome territory (CT) arrangement and somatic homologous pairing in interphase nuclei seem to occur randomly. Since chromosomes containing nucleolar organizer regions associate more frequently to form a single nucleolar structure, as in yeast, the nucleolus seems to play a major role in organizing nuclear space. Recent findings have begun to elucidate how plant regulatory factors, such as chromatin remodeling or histone chaperones, affect the chromatin state of ribosomal DNA genes located in two distinct CT arrangements in the nucleus. The functional outcome of yeast nuclear organization allowed us to propose how nuclear organization might contribute to a novel type of epigenetic regulation: the spatial regulation of transcription.

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.

Identification of protein factors and U3 snoRNAs from a Brassica oleracea RNP complex involved in the processing of pre-rRNA

We report on the structural characterization of a functional U3 snoRNA ribonucleoprotein complex isolated from Brassica oleracea. The BoU3 snoRNP complex (formerly NF D) binds ribosomal DNA (rDNA), specifically cleaves pre-rRNA at the primary cleavage site in vitro and probably links transcription to early pre-rRNA processing in vivo. Using a proteomic approach we have identified 62 proteins in the purified BoU3 snoRNP fraction, including small RNA associated proteins (Fibrillarin, NOP5/Nop58p, Diskerin/Cbf5p, SUS2/PRP8 and CLO/GFA1/sn114p) and 40S ribosomal associated proteins (22 RPS and four ARCA-like proteins). Another major protein group is composed of chaperones/chaperonins (HSP81/TCP-1) and at least one proteasome subunit (RPN1a). Remarkably, RNA-dependent RNA polymerase (RdRP) and Tudor staphylococcal nuclease (TSN) proteins, which have RNA- and/or DNA-associated activities, were also revealed in the complex. Furthermore, three U3 snoRNA variants were identified in the BoU3 snoRNP fraction, notably an evolutionarily conserved and variable stem loop structure located just downstream from the C-box domain of the U3 sequence structures. We conclude that the BoU3 snoRNP complex is mainly required for 40S pre-ribosome synthesis. It is also expected that U3 snoRNA variants and interacting proteins might play a major role in BoU3 snoRNP complex assembly and/or function. This study provides a basis for further investigation of these novel ribonucleoprotein factors and their role in plant ribosome biogenesis.

SnoRNAs from the filamentous fungus Neurospora crassa: structural, functional and evolutionary insights

BACKGROUND

SnoRNAs represent an excellent model for studying the structural and functional evolution of small non-coding RNAs involved in the post-transcriptional modification machinery for rRNAs and snRNAs in eukaryotic cells. Identification of snoRNAs from Neurospora crassa, an important model organism playing key roles in the development of modern genetics, biochemistry and molecular biology will provide insights into the evolution of snoRNA genes in the fungus kingdom.

RESULTS

Fifty five box C/D snoRNAs were identified and predicted to guide 71 2'-O-methylated sites including four sites on snRNAs and three sites on tRNAs. Additionally, twenty box H/ACA snoRNAs, which potentially guide 17 pseudouridylations on rRNAs, were also identified. Although not exhaustive, the study provides the first comprehensive list of two major families of snoRNAs from the filamentous fungus N. crassa. The independently transcribed strategy dominates in the expression of box H/ACA snoRNA genes, whereas most of the box C/D snoRNA genes are intron-encoded. This shows that different genomic organizations and expression modes have been adopted by the two major classes of snoRNA genes in N. crassa . Remarkably, five gene clusters represent an outstanding organization of box C/D snoRNA genes, which are well conserved among yeasts and multicellular fungi, implying their functional importance for the fungus cells. Interestingly, alternative splicing events were found in the expression of two polycistronic snoRNA gene hosts that resemble the UHG-like genes in mammals. Phylogenetic analysis further revealed that the extensive separation and recombination of two functional elements of snoRNA genes has occurred during fungus evolution.

CONCLUSION

This is the first genome-wide analysis of the filamentous fungus N. crassa snoRNAs that aids in understanding the differences between unicellular fungi and multicellular fungi. As compared with two yeasts, a more complex pattern of methylation guided by box C/D snoRNAs in multicellular fungus than in unicellular yeasts was revealed, indicating the high diversity of post-transcriptional modification guided by snoRNAs in the fungus kingdom.

Eukaryotic snoRNAs: a paradigm for gene expression flexibility

Small nucleolar RNAs (snoRNAs) are one of the most ancient and numerous families of non-protein-coding RNAs (ncRNAs). The main function of snoRNAs - to guide site-specific rRNA modification - is the same in Archaea and all eukaryotic lineages. In contrast, as revealed by recent genomic and RNomic studies, their genomic organization and expression strategies are the most varied. Seemingly snoRNA coding units have adopted, in the course of evolution, all the possible ways of being transcribed, thus providing a unique paradigm of gene expression flexibility. By focusing on representative fungal, plant and animal genomes, we review here all the documented types of snoRNA gene organization and expression, and we provide a comprehensive account of snoRNA expressional freedom by precisely estimating the frequency, in each genome, of each type of genomic organization. We finally discuss the relevance of snoRNA genomic studies for our general understanding of ncRNA family evolution and expression in eukaryotes.

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.

Psiscan: a computational approach to identify H/ACA-like and AGA-like non-coding RNA in trypanosomatid genomes

Background

Detection of non coding RNA (ncRNA) molecules is a major bioinformatics challenge. This challenge is particularly difficult when attempting to detect H/ACA molecules which are involved in converting uridine to pseudouridine on rRNA in trypanosomes, because these organisms have unique H/ACA molecules (termed H/ACA-like) that lack several of the features that characterize H/ACA molecules in most other organisms.

Results

We present here a computational tool called Psiscan, which was designed to detect H/ACA-like molecules in trypanosomes. We started by analyzing known H/ACA-like molecules and characterized their crucial elements both computationally and experimentally.

Next, we set up constraints based on this analysis and additional phylogenic and functional data to rapidly scan three trypanosome genomes (T. brucei, T. cruzi and L. major) for sequences that observe these constraints and are conserved among the species. In the next step, we used minimal energy calculation to select the molecules that are predicted to fold into a lowest energy structure that is consistent with the constraints. In the final computational step, we used a Support Vector Machine that was trained on known H/ACA-like molecules as positive examples and on negative examples of molecules that were identified by the computational analyses but were shown experimentally not to be H/ACA-like molecules. The leading candidate molecules predicted by the SVM model were then subjected to experimental validation.

Conclusion

The experimental validation showed 11 molecules to be expressed (4 out of 25 in the intermediate stage and 7 out of 19 in the final validation after the machine learning stage). Five of these 11 molecules were further shown to be bona fide H/ACA-like molecules. As snoRNA in trypanosomes are organized in clusters, the new H/ACA-like molecules could be used as starting points to manually search for additional molecules in their neighbourhood. All together this study increased our repertoire by fourteen H/ACA-like and six C/D snoRNAs molecules from T. brucei and L. Major. In addition the experimental analysis revealed that six ncRNA molecules that are expressed are not downregulated in CBF5 silenced cells, suggesting that they have structural features of H/ACA-like molecules but do not have their standard function. We termed this novel class of molecules AGA-like, and we are exploring their function.

This study demonstrates the power of tight collaboration between computational and experimental approaches in a combined effort to reveal the repertoire of ncRNA molecles.

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.

Exploration of small RNAs

For several decades, only a limited number of noncoding RNAs, such as ribosomal and transfer RNA, have been studied in any depth. In recent years, additional species of noncoding RNAs have increasingly been discovered. Of these, small RNA species attract particular interest because of their essential roles in processes such as RNA silencing and modifications. Detailed analyses revealed several pathways associated with the function of small RNAs. Although these pathways show evolutional conservation, there are substantial differences. Advanced technologies to profile RNAs have accelerated the field further resulting in the discovery of an increasing number of novel species, suggesting that we are only just beginning to appreciate the complexity of small RNAs and their functions. Here, we review recent progress in novel small RNA exploration, including discovered small RNA species, their pathways, and devised technologies.

The role of the plant nucleolus in pre-mRNA processing

The nucleolus is a multifunctional compartment of the eukaryotic nucleus. Besides its well-recognised role in transcription and processing of ribosomal RNA and the assembly of ribosomal subunits, the nucleolus has functions in the processing and assembly of a variety of RNPs and is involved in cell cycle control and senescence and as a sensor of stress. Historically, nucleoli have been tenuously linked to the biogenesis and, in particular, export of mRNAs in yeast and mammalian cells. Recently, data from plants have extended the functions in which the plant nucleolus is involved to include transcriptional gene silencing as well as mRNA surveillance and nonsense-mediated decay, and mRNA export. The nucleolus in plants may therefore have important roles in the biogenesis and quality control of mRNAs.

snoTARGET shows that human orphan snoRNA targets locate close to alternative splice junctions

Among thousands of non-protein-coding RNAs which have been found in humans, a significant group represents snoRNA molecules that guide other types of RNAs to specific chemical modifications, cleavages, or proper folding. Yet, hundreds of mammalian snoRNAs have unknown function and are referred to as "orphan" molecules. In 2006, for the first time, it was shown that a particular orphan snoRNA (HBII-52) plays an important role in the regulation of alternative splicing of the serotonin receptor gene in humans and other mammals. In order to facilitate the investigation of possible involvement of snoRNAs in the regulation of pre-mRNA processing, we developed a new computational web resource, snoTARGET, which searches for possible guiding sites for snoRNAs among the entire set of human and rodent exonic and intronic sequences. Application of snoTARGET for finding possible guiding sites for a number of human and rodent orphan C/D-box snoRNAs showed that another subgroup of these molecules (HBII-85) have statistically elevated guiding preferences toward exons compared to introns. Moreover, these energetically favorable putative targets of HBII-85 snoRNAs are non-randomly associated with genes producing alternatively spliced mRNA isoforms. The snoTARGET resource is freely available at: (http://hsc.utoledo.edu/depts/bioinfo/snotarget.html).

Genome-wide analyses of retrogenes derived from the human box H/ACA snoRNAs

The family of box H/ACA snoRNA is an abundant class of non-protein-coding RNAs, which play important roles in the post-transcriptional modification of rRNAs and snRNAs. Here we report the characterization in the human genome of 202 sequences derived from box H/ACA snoRNAs. Most of them were retrogenes formed using the L1 integration machinery. About 96% of the box H/ACA RNA-related sequences are found in corresponding locations on the chimpanzee and human chromosomes, while the mouse shares approximately 50% of these human sequences, suggesting that some of the H/ACA RNA-related sequences in primate occurred after the rodent/primate divergence. Of the H/ACA RNA-related sequences, 49% are found in intronic regions of protein-coding genes and 64 H/ACA-related sequences can be folded to the typical secondary structure of the box H/ACA snoRNA family, while 30 of them were recognized as functional homologs of their corresponding box H/ACA snoRNAs previously reported. Of the 64 sequences with the typical secondary structure of the box H/ACA RNA family, 11 were found in EST databases and 5 among which were shown to be expressed in more than one human tissue. Notably, U107f is nested in an intron of a protein gene coding for nudix-type motif 13, but expressed from the opposite strand, and the searching of EST databases revealed it can be expressed in liver and spleen, even in melanotic melanoma.

Modifications in Thermus thermophilus 23 S ribosomal RNA are centered in regions of RNA-RNA contact

Ribosomal RNA from all organisms contains post-transcriptionally modified nucleotides whose function is far from clear. To gain insight into the molecular interactions of modified nucleotides, we investigated the modification status of Thermus thermophilus 5 S and 23 S ribosomal RNA by mass spectrometry and chemical derivatization/primer extension. A total of eleven modified nucleotides was found in 23 S rRNA, of which eight were singly methylated nucleotides and three were pseudouridines. These modified nucleotides were mapped into the published three-dimensional ribosome structure. Seven of the modified nucleotides located to domain IV, and four modified nucleotides located to domain V of the 23 S rRNA. All posttranscriptionally modified nucleotides map in the center of the ribosome, and none of them are in contact with ribosomal proteins. All except one of the modified nucleotides were found in secondary structure elements of the 23 S ribosomal RNA that contact either 16 S ribosomal RNA or transfer RNA, with five of these nucleotides physically involved in intermolecular RNA-RNA bridges. These findings strongly suggest that the post-transcriptional modifications play a role in modulating intermolecular RNA-RNA contacts, which is the first suggestion on a specific function of endogenous ribosomal RNA modifications.

Computer simulation of chaperone effects of Archaeal C/D box sRNA binding on rRNA folding

Archaeal C/D box small RNAs (sRNAs) are homologues of eukaryotic C/D box small nucleolar RNAs (snoRNAs). Their main function is guiding 2'-O-ribose methylation of nucleotides in rRNAs. The methylation requires the pairing of an sRNA antisense element to an rRNA target site with formation of an RNA-RNA duplex. The temporary formation of such a duplex during rRNA maturation is expected to influence rRNA folding in a chaperone-like way, in particular in thermophilic Archaea, where multiple sRNAs with two binding sites are found. Here we investigate possible mechanisms of chaperone function of Archaeoglobus fulgidus and Pyrococcus abyssi C/D box sRNAs using computer simulations of rRNA secondary structure formation by genetic algorithm. The effects of sRNA binding on rRNA structure are introduced as temporary structural constraints during co-transcriptional folding. Comparisons of the final predictions with simulations without sRNA binding and with phylogenetic structures show that sRNAs with two antisense elements may significantly facilitate the correct formation of long-range interactions in rRNAs, in particular at elevated temperatures. The simulations suggest that the main mechanism of this effect is a transient restriction of folding in rRNA domains where the termini are brought together by binding to double-guide sRNAs.

Identification of novel non-coding RNAs as potential antisense regulators in the archaeon Sulfolobus solfataricus

By generating a specialized cDNA library from the archaeon Sulfolobus solfataricus, we have identified 57 novel small non-coding RNA (ncRNA) candidates and confirmed their expression by Northern blot analysis. The majority was found to belong to one of two classes, either antisense or antisense-box RNAs, where the latter only exhibit partial complementarity to RNA targets. The most prominent group of antisense RNAs is transcribed in the opposite orientation to the transposase genes, encoded by insertion elements (transposons). Thus, these antisense RNAs may regulate transposition of insertion elements by inhibiting expression of the transposase mRNA. Surprisingly, the class of antisense RNAs also contained RNAs complementary to tRNAs or sRNAs (small-nucleolar-like RNAs). For the antisense-box ncRNAs, the majority could be assigned to the class of C/D sRNAs, which specify 2'-O-methylation sites on rRNAs or tRNAs. Five C/D sRNAs of this group are predicted to target methylation at six sites in 13 different tRNAs, thus pointing to the widespread role of these sRNA species in tRNA modification in Archaea. Another group of antisense-box RNAs, lacking typical C/D sRNA motifs, was predicted to target the 3'-untranslated regions of certain mRNAs. Furthermore, one of the ncRNAs that does not show antisense elements is transcribed from a repeat unit of a cluster of small regularly spaced repeats in S. solfataricus which is potentially involved in replicon partitioning. In conclusion, this is the first report of stably expressed antisense RNAs in an archaeal species and it raises the prospect that antisense-based mechanisms are also used widely in Archaea to regulate gene expression.

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.

In vitro RNP assembly and methylation guide activity of an unusual box C/D RNA, cis-acting archaeal pre-tRNA(Trp)

Among the large family of C/D methylation guide RNAs, the intron of euryarchaeal pre-tRNA(Trp) represents an outstanding specimen able to guide in cis, instead of in trans, two 2'-O-methylations in the pre-tRNA exons. Remarkably, both sites of methylation involve nucleotides within the bulge-helix-bulge (BHB) splicing motif, while the RNA-guided methylation and pre-tRNA splicing events depend on mutually exclusive RNA folding patterns. Using the three recombinant core proteins of archaeal C/D RNPs, we have analyzed in vitro RNP assembly of the pre-tRNA and tested its site-specific methylation activity. Recognition by L7Ae of hallmark K-turns at the C/D and C'/D' motifs appears as a crucial assembly step required for subsequent binding of a Nop5p-aFib heterodimer at each site. Unexpectedly, however, even without L7Ae but at a higher concentration of Nop5p-aFib, a substantially active RNP complex can still form, possibly reflecting the higher propensity of the cis-acting system to form guide RNA duplex(es) relative to classical trans- acting C/D RNA guides. Moreover, footprinting data of RNPs, consistent with Nop5p interacting with the non-canonical stem of the K-turn, suggest that binding of Nop5p-aFib to the pre-tRNA-L7Ae complex might direct transition from a splicing-competent structure to an RNA conformer displaying the guide RNA duplexes required for site-specific methylation.

Identification of 13 novel human modification guide RNAs

Members of the two expanding RNA subclasses termed C/D and H/ACA RNAs guide the 2'-O-methylations and pseudouridylations, respectively, of rRNA and spliceosomal RNAs (snRNAs). Here, we report on the identification of 13 novel human intron-encoded small RNAs (U94-U106) belonging to the two subclasses of modification guides. Seven of them are predicted to direct 2'-O-methylations in rRNA or snRNAs, while the remainder represent novel orphan RNA modification guides. From these, U100, which is exclusively detected in Cajal bodies (CBs), is predicted to direct modification of a U6 snRNA uridine, U(9), which to date has not been found to be pseudouridylated. Hence, within CBs, U100 might function in the folding pathway or other aspects of U6 snRNA metabolism rather than acting as a pseudouridylation guide. U106 C/D snoRNA might also possess an RNA chaperone activity only since its two conserved antisense elements match two rRNA sequences devoid of methylated nucleotides and located remarkably close to each other within the 18S rRNA secondary structure. Finally, we have identified a retrogene for U99 snoRNA located within an intron of the Siat5 gene, supporting the notion that retro-transposition events might have played a substantial role in the mobility and diversification of snoRNA genes during evolution.

Hyper-expression of small nucleolar RNAs (snoRNAs) in female inflorescences of hazelnut (Corylus avellana L.) supports rRNA aggregation in vitro

Under certain in vitro (salt and temperature) conditions rRNA aggregation occurs in female inflorescences but not in leaves or pollen RNA preparations from hazelnut (Corylus avellana L.), a species of economic interest. This paper describes experiments addressing an explanation of this phenomenon. The experiments demonstrate that: (i) trans-acting factors induce rRNA aggregate formation in female inflorescences RNA preparations; (ii) these factors support aggregation also of heterologous rRNA; (iii) aggregation is a function of temperature pre-treatment of rRNA and not of source 18S rRNA; (iv) the factors inducing rRNA aggregates are sensitive to RNase; (v) antisense small nucleolar RNAs (snoRNAs) participate in rRNA aggregate formation. snoRNAs are involved in pre-rRNA spacer cleavages, and are required for the two most common types of rRNA modifications: 2'-O-ribose methylation and pseudouridylation. Even though it is questionable whether rRNA aggregation really happens in female inflorescence in vivo, the phenomenon observed in vitro may reflect the abundance of snoRNAs in these reproductive structures. In fact the level of accumulation of three tested snoRNAs, R1, U14 and U3, is much higher in female inflorescence than in leaves or pollen of hazelnut. This finding opens the possibility of studying the role of snoRNAs in tissue development in plants.

Evolutionary profiling of the U49 snoRNA gene

Small nucleolar RNAs (snoRNAs) are involved in precursor ribosomal RNA (pre-rRNA) processing and rRNA base modifications (2'-O-ribose methylation and pseudouridylation). Their genomic organization show great flexibility: some are individually or polycistronically transcribed, while others are encoded within introns of other genes. Here, we present an evolutionary analysis of the U49 gene in seven species. In all species analyzed, U49 contains the typical hallmarks of C and D box motifs, and a conserved 12-15 nt sequence complementary to rRNA that define them as homologs. In mouse, human, and Drosophila U49 is found encoded within introns of different genes, and in plants it is transcribed polycistronically from four different locations. In addition, U49 has two copies in two different introns of the RpL14 gene in Drosophila. The results indicate a substantial degree of duplication and translocation of the U49 gene in evolution. In light of its variable organization we discuss which of the two proposed mechanisms of rearrangement has acted upon the U49 snoRNA gene: chromosomal duplication or transposition through an RNA intermediate.

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.

Plant dicistronic tRNA-snoRNA genes: a new mode of expression of the small nucleolar RNAs processed by RNase Z

Small nucleolar RNAs (snoRNAs) guiding modifications of ribosomal RNAs and other RNAs display diverse modes of gene organization and expression depending on the eukaryotic system: in animals most are intron encoded, in yeast many are monocistronic genes and in plants most are polycistronic (independent or intronic) genes. Here we report an unprecedented organization: plant dicistronic tRNA-snoRNA genes. In Arabidopsis thaliana we identified a gene family encoding 12 novel box C/D snoRNAs (snoR43) located just downstream from tRNA(Gly) genes. We confirmed that they are transcribed, probably from the tRNA gene promoter, producing dicistronic tRNA(Gly)-snoR43 precursors. Using transgenic lines expressing a tagged tRNA-snoR43.1 gene we show that the dicistronic precursor is accurately processed to both snoR43.1 and tRNA(Gly). In addition, we show that a recombinant RNase Z, the plant tRNA 3' processing enzyme, efficiently cleaves the dicistronic precursor in vitro releasing the snoR43.1 from the tRNA(Gly). Finally, we describe a similar case in rice implicating a tRNA(Met-e) expressed in fusion with a novel C/D snoRNA, showing that this mode of snoRNA expression is found in distant plant species.

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.

Plant snoRNAs: functional evolution and new modes of gene expression

Small nucleolar RNAs (snoRNAs) are a well-characterized family of non-coding RNAs whose main function is rRNA modification. The diversity and complexity of this gene family continues to expand with the discovery of snoRNAs with non-rRNA or unknown targets. Plants contain more snoRNAs than other eukaryotes and have developed novel expression and processing strategies. The increased number of modifications, which will influence ribosome function, and the novel modes of expression might reflect the environmental conditions to which plants are exposed. Polyploidy and chromosomal rearrangements have generated multiple copies of snoRNA genes, allowing the generation of new snoRNAs for selection. The large snoRNA family in plants is an ideal model for investigation of mechanisms of evolution of gene families in plants.

Experimental RNomics: identification of 140 candidates for small non-messenger RNAs in the plant Arabidopsis thaliana

BACKGROUND

Genomes from all organisms known to date express two types of RNA molecules: messenger RNAs (mRNAs), which are translated into proteins, and non-messenger RNAs, which function at the RNA level and do not serve as templates for translation.

RESULTS

We have generated a specialized cDNA library from Arabidopsis thaliana to investigate the population of small non-messenger RNAs (snmRNAs) sized 50-500 nt in a plant. From this library, we identified 140 candidates for novel snmRNAs and investigated their expression, abundance, and developmental regulation. Based on conserved sequence and structure motifs, 104 snmRNA species can be assigned to novel members of known classes of RNAs (designated Class I snmRNAs), namely, small nucleolar RNAs (snoRNAs), 7SL RNA, U snRNAs, as well as a tRNA-like RNA. For the first time, 39 novel members of H/ACA box snoRNAs could be identified in a plant species. Of the remaining 36 snmRNA candidates (designated Class II snmRNAs), no sequence or structure motifs were present that would enable an assignment to a known class of RNAs. These RNAs were classified based on their location on the A. thaliana genome. From these, 29 snmRNA species located to intergenic regions, 3 located to intronic sequences of protein coding genes, and 4 snmRNA candidates were derived from annotated open reading frames. Surprisingly, 15 of the Class II snmRNA candidates were shown to be tissue-specifically expressed, while 12 are encoded by the mitochondrial or chloroplast genome.

CONCLUSIONS

Our study has identified 140 novel candidates for small non-messenger RNA species in the plant A. thaliana and thereby sets the stage for their functional analysis.

RNomics: identification and function of small, non-messenger RNAs

In the past few years, our knowledge about small non-mRNAs (snmRNAs) has grown exponentially. Approaches including computational and experimental RNomics have led to a plethora of novel snmRNAs, especially small nucleolar RNAs (snoRNAs). Members of this RNA class guide modification of ribosomal and spliceosomal RNAs. Novel targets for snoRNAs were identified such as tRNAs and potentially mRNAs, and several snoRNAs were shown to be tissue-specifically expressed. In addition, previously unknown classes of snmRNAs have been discovered. MicroRNAs and small interfering RNAs of about 21-23 nt, were shown to regulate gene expression by binding to mRNAs via antisense elements. Regulation of gene expression is exerted by degradation of mRNAs or translational regulation. snmRNAs play a variety of roles during regulation of gene expression. Moreover, the function of some snmRNAs known for decades, has been finally elucidated. Many other RNAs were identified by RNomics studies lacking known sequence and structure motifs. Future challenges in the field of RNomics include identification of the novel snmRNA's biological roles in the cell.

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.

A novel gene organization: intronic snoRNA gene clusters from Oryza sativa

Based on the analysis of structural features and conserved elements, 27 novel snoRNA genes have been identified from rice. All of them belong to the C/D box-containing snoRNA family except for one that belongs to the H/ACA box type. The newly found genes fall into six clusters that comprise at least three snoRNA genes, and in one case as many as nine genes. Interestingly, four of the six clusters are located within the largest intron of a protein coding gene. The majority of intronic snoRNA gene clusters are simply formed by multiple copies of the same species of snoRNA gene that possess the identical functional elements. This implies a possible mechanism of duplication for the origin of repeating snoRNA coding regions in one intron. However, a few intronic snoRNA gene clusters consisting of different snoRNAs species were also observed. Polycistronic precursors from two independently transcribed clusters were demonstrated by RT-PCR and individual snoRNAs processed from the polycistronic precursors were positively determined by reverse transcription assay. Analyses of the intergenic spacers in the clusters showed that, in addition to a very high AT content, the processing signals in rice snoRNA polycistronic transcripts might be different from those of yeast. Our results demonstrate that, in both plants and mammals, numerous snoRNAs can be produced simultaneously from an mRNA precursor of a host gene despite the different arrangements. The intronic snoRNA gene cluster is a novel gene organization, which is so far unique to plants. The conservation of intronic snoRNA gene clusters in plants was further demonstrated by the study of a similar snoRNA gene organization in the first intron of a Hsp70 gene from wild rice and Zizania caduciflora.