Identification of 10 novel snoRNA gene clusters from Arabidopsis thaliana

A contemporary review of snoRNAs in cardiovascular health: RNA modification and beyond

As cardiovascular diseases continue to be the leading cause of death worldwide, groundbreaking research is being conducted to mitigate their effects. This review looks into the potential of small nucleolar RNAs (snoRNAs) and the opportunity to use these molecular agents as therapeutic biomarkers for cardiovascular issues specific to the heart. Through an investigation of snoRNA biogenesis, functionality, and roles in cardiovascular diseases, this review relates our past and present knowledge of snoRNAs to the current scientific literature. Considering the initial discovery of snoRNAs and the studies thereafter analyzing the roles of snoRNAs in disease, we look forward to uncovering many other noncanonical functions that could lead researchers closer to finding preventive and curative solutions for cardiovascular diseases.

Plant ribosomes as a score to fathom the melody of 2'-O-methylation across evolution

2'-O-ribose methylation (2'-O-Me) is one of the most common RNA modifications detected in ribosomal RNAs (rRNA) from bacteria to eukaryotic cells. 2'-O-Me favours a specific RNA conformation and protects RNA from hydrolysis. Moreover, rRNA 2'-O-Me might stabilize its interactions with messenger RNA (mRNA), transfer RNA (tRNA) or proteins. The extent of rRNA 2'-O-Me fluctuates between species from 3-4 sites in bacteria to tens of sites in archaea, yeast, algae, plants and human. Depending on the organism as well as the rRNA targeting site and position, the 2'-O-Me reaction can be carried out by several site-specific RNA methyltransferases (RMTase) or by a single RMTase associated to specific RNA guides. Here, we review current progresses in rRNA 2'-O-Me (sites/Nm and RMTases) in plants and compare the results with molecular clues from unicellular (bacteria, archaea, algae and yeast) as well as multicellular (human and plants) organisms.

High-Resolution Structure and Internal Mobility of a Plant 40S Ribosomal Subunit

Ribosome is a major part of the protein synthesis machinery, and analysis of its structure is of paramount importance. However, the structure of ribosomes from only a limited number of organisms has been resolved to date; it especially concerns plant ribosomes and ribosomal subunits. Here, we report a high-resolution cryo-electron microscopy reconstruction of the small subunit of the Triticum aestivum (common wheat) cytoplasmic ribosome. A detailed atomic model was built that includes the majority of the rRNA and some of the protein modifications. The analysis of the obtained data revealed structural peculiarities of the 40S subunit in the monocot plant ribosome. We applied the 3D Flexible Refinement approach to analyze the internal mobility of the 40S subunit and succeeded in decomposing it into four major motions, describing rotations of the head domain and a shift in the massive rRNA expansion segment. It was shown that these motions are almost uncorrelated and that the 40S subunit is flexible enough to spontaneously adopt any conformation it takes as a part of a translating ribosome or ribosomal complex. Here, we introduce the first high-resolution structure of an isolated plant 40S subunit and the first quantitative analysis of the flexibility of small ribosomal subunits, hoping that it will help in studying various aspects of ribosome functioning.

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.

A Conserved Long Intergenic Non-coding RNA Containing snoRNA Sequences, lncCOBRA1, Affects Arabidopsis Germination and Development

Long non-coding RNAs (lncRNAs) are an increasingly studied group of non-protein coding transcripts with a wide variety of molecular functions gaining attention for their roles in numerous biological processes. Nearly 6,000 lncRNAs have been identified in Arabidopsis thaliana but many have yet to be studied. Here, we examine a class of previously uncharacterized lncRNAs termed CONSERVED IN BRASSICA RAPA (lncCOBRA) transcripts that were previously identified for their high level of sequence conservation in the related crop species Brassica rapa, their nuclear-localization and protein-bound nature. In particular, we focus on lncCOBRA1 and demonstrate that its abundance is highly tissue and developmental specific, with particularly high levels early in germination. lncCOBRA1 contains two snoRNAs domains within it, making it the first sno-lincRNA example in a non-mammalian system. However, we find that it is processed differently than its mammalian counterparts. We further show that plants lacking lncCOBRA1 display patterns of delayed germination and are overall smaller than wild-type plants. Lastly, we identify the proteins that interact with lncCOBRA1 and propose a novel mechanism of lincRNA action in which it may act as a scaffold with the RACK1A protein to regulate germination and development, possibly through a role in ribosome biogenesis.

Genome-Wide Identification and Analysis of Small Nucleolar RNAs and Their Roles in Regulating Latex Regeneration in the Rubber Tree (Hevea brasiliensis)

Small nucleolar RNAs (snoRNAs) are a class of conserved nuclear RNAs that play important roles in the modification of ribosomal RNAs (rRNAs) in plants. In rubber trees, rRNAs are run off with latex flow during tapping and need to be regenerated for maintaining the functions of the laticifer cells. SnoRNAs are expected to play essential roles in the regeneration of rRNAs. However, snoRNAs in the rubber tree have not been sufficiently characterized thus far. In this study, we performed nuclear RNA sequencing (RNA-seq) to identify snoRNAs globally and investigate their roles in latex regeneration. We identified a total of 3,626 snoRNAs by computational prediction with nuclear RNA-seq data. Among these snoRNAs, 50 were highly expressed in latex; furthermore, the results of reverse transcription polymerase chain reaction (RT-PCR) showed the abundant expression of 31 of these snoRNAs in latex. The correlation between snoRNA expression and adjusted total solid content (TSC/C) identified 13 positively yield-correlated snoRNAs. To improve the understanding of latex regeneration in rubber trees, we developed a novel insulated tapping system (ITS), which only measures the latex regenerated in specific laticifers. Using this system, a laticifer-abundant snoRNA, HbsnoR28, was found to be highly correlated with latex regeneration. To the best of our knowledge, this is the first report to globally identify snoRNAs that might be involved in latex regeneration regulation and provide new clues for unraveling the mechanisms underlying the regulation of latex regeneration.

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.

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 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.

Ribosome Biogenesis in Plants: From Functional 45S Ribosomal DNA Organization to Ribosome Assembly Factors

The transcription of 18S, 5.8S, and 18S rRNA genes (45S rDNA), cotranscriptional processing of pre-rRNA, and assembly of mature rRNA with ribosomal proteins are the linchpins of ribosome biogenesis. In yeast (Saccharomyces cerevisiae) and animal cells, hundreds of pre-rRNA processing factors have been identified and their involvement in ribosome assembly determined. These studies, together with structural analyses, have yielded comprehensive models of the pre-40S and pre-60S ribosome subunits as well as the largest cotranscriptionally assembled preribosome particle: the 90S/small subunit processome. Here, we present the current knowledge of the functional organization of 45S rDNA, pre-rRNA transcription, rRNA processing activities, and ribosome assembly factors in plants, focusing on data from Arabidopsis (Arabidopsis thaliana). Based on yeast and mammalian cell studies, we describe the ribonucleoprotein complexes and RNA-associated activities and discuss how they might specifically affect the production of 40S and 60S subunits. Finally, we review recent findings concerning pre-rRNA processing pathways and a novel mechanism involved in a ribosome stress response in plants.

Arabidopsis thaliana NOP10 is required for gametophyte formation

The female gametophyte is crucial for sexual reproduction of higher plants, yet little is known about the molecular mechanisms underlying its development. Here, we report that Arabidopsis thaliana NOP10 (AtNOP10) is required for female gametophyte formation. AtNOP10 was expressed predominantly in the seedling and reproductive tissues, including anthers, pollen grains, and ovules. Mutations in AtNOP10 interrupted mitosis of the functional megaspore during early development and prevented polar nuclear fusion in the embryo sacs. AtNOP10 shares a high level of amino acid sequence similarity with Saccharomyces cerevisiae (yeast) NOP10 (ScNOP10), an important component of the H/ACA small nucleolar ribonucleoprotein particles (H/ACA snoRNPs) implicated in 18S rRNA synthesis and rRNA pseudouridylation. Heterologous expression of ScNOP10 complemented the mutant phenotype of Atnop10. Thus, AtNOP10 influences functional megaspore mitosis and polar nuclear fusion during gametophyte formation in Arabidopsis.

Canonical and Noncanonical Actions of Arabidopsis Histone Deacetylases in Ribosomal RNA Processing

Ribosome biogenesis is a fundamental process required for all cellular activities. Histone deacetylases play critical roles in many biological processes including transcriptional repression and rDNA silencing. However, their function in pre-rRNA processing remains poorly understood. Here, we discovered a previously uncharacterized role of Arabidopsis thaliana histone deacetylase HD2C in pre-rRNA processing via both canonical and noncanonical manners. HD2C interacts with another histone deacetylase HD2B and forms homo- and/or hetero-oligomers in the nucleolus. Depletion of HD2C and HD2B induces a ribosome-biogenesis deficient phenotype and aberrant accumulation of 18S pre-rRNA intermediates. Our genome-wide analysis revealed that HD2C binds and represses the expression of key genes involved in ribosome biogenesis. Using RNA immunoprecipitation and sequencing, we further uncovered a noncanonical mechanism of HD2C directly associating with pre-rRNA and small nucleolar RNAs to regulate rRNA methylation. Together, this study reveals a multifaceted role of HD2C in ribosome biogenesis and provides mechanistic insights into how histone deacetylases modulate rRNA maturation at the transcriptional and posttranscriptional levels.

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.

Arabidopsis small nucleolar RNA monitors the efficient pre-rRNA processing during ribosome biogenesis

Ribosome production in eukaryotes requires the complex and precise coordination of several hundred assembly factors, including many small nucleolar RNAs (snoRNAs). However, at present, the distinct role of key snoRNAs in ribosome biogenesis remains poorly understood in higher plants. Here we report that a previously uncharacterized C (RUGAUGA)/D (CUGA) type snoRNA, HIDDEN TREASURE 2 (HID2), acts as an important regulator of ribosome biogenesis through a snoRNA-rRNA interaction. Nucleolus-localized HID2 is actively expressed in Arabidopsis proliferative tissues, whereas defects in HID2 cause a series of developmental defects reminiscent of ribosomal protein mutants. HID2 associates with the precursor 45S rRNA and promotes the efficiency and accuracy of pre-rRNA processing. Intriguingly, disrupting HID2 in Arabidopsis appears to impair the integrity of 27SB, a key pre-rRNA intermediate that generates 25S and 5.8S rRNA and is known to be vital for the synthesis of the 60S large ribosomal subunit and also produces an imbalanced ribosome profile. Finally, we demonstrate that the antisense-box of HID2 is both functionally essential and highly conserved in eukaryotes. Overall, our study reveals the vital and possibly conserved role of a snoRNA in monitoring the efficiency of pre-rRNA processing during ribosome biogenesis.

Distinct Profiles for Mitochondrial t-RNAs and Small Nucleolar RNAs in Locally Invasive and Metastatic Colorectal Cancer

PURPOSE

To gain insight into factors involved in tumor progression and metastasis, we examined the role of noncoding RNAs in the biologic characteristics of colorectal carcinoma, in paired samples of tumor together with normal mucosa from the same colorectal carcinoma patient. The tumor and healthy tissue samples were collected and stored under stringent conditions, thereby minimizing warm ischemic time.

EXPERIMENTAL DESIGN

We focused particularly on distinctions among high-stage tumors and tumors with known metastases, performing RNA-Seq analysis that quantifies transcript abundance and identifies novel transcripts.

RESULTS

In comparing 35 colorectal carcinomas, including 9 metastatic tumors (metastases to lymph nodes and lymphatic vessels), with their matched healthy control mucosa, we found a distinct signature of mitochondrial transfer RNAs (MT-tRNA) and small nucleolar RNAs (snoRNA) for metastatic and high-stage colorectal carcinoma. We also found the following: (i) MT-TF (phenylalanine) and snord12B expression correlated with a substantial number of miRNAs and mRNAs in 14 colorectal carcinomas examined; (ii) an miRNA signature of oxidative stress, hypoxia, and a shift to glycolytic metabolism in 14 colorectal carcinomas, regardless of grade and stage; and (iii) heterogeneous MT-tRNA/snoRNA fingerprints for 35 pairs.

CONCLUSIONS

These findings could potentially assist in more accurate and predictive staging of colorectal carcinoma, including identification of those colorectal carcinomas likely to metastasize.

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.

Experimental RNomics and genomic comparative analysis reveal a large group of species-specific small non-message RNAs in the silkworm Bombyx mori

Accumulating evidences show that small non-protein coding RNAs (ncRNAs) play important roles in development, stress response and other cellular processes. The silkworm is an important model for studies on insect genetics and control of lepidopterous pests. Here, we have performed the first systematic identification and analysis of intermediate size ncRNAs (50–500 nt) in the silkworm. We identified 189 novel ncRNAs, including 141 snoRNAs, six snRNAs, three tRNAs, one SRP and 38 unclassified ncRNAs. Forty ncRNAs showed significantly altered expression during silkworm development or across specific stage transitions. Genomic comparisons revealed that 123 of these ncRNAs are potentially silkworm-specific. Analysis of the genomic organization of the ncRNA loci showed that 32.62% of the novel snoRNA loci are intergenic, and that all the intronic snoRNAs follow the pattern of one-snoRNA-per-intron. Target site analysis predicted a total of 95 2′-O-methylation and pseudouridylation modification sites of rRNAs, snRNAs and tRNAs. Together, these findings provide new clues for future functional study of ncRNA during insect development and evolution.

Post-transcriptional modification of RNAs by artificial Box H/ACA and Box C/D RNPs

RNA-guided RNA 2'-O-methylation and pseudouridylation are naturally occurring processes, in which guide RNAs specifically direct modifications to rRNAs or spliceosomal snRNAs in the nucleus of eukaryotic cells. Modifications can profoundly alter the properties of an RNA, thus influencing the contributions of the RNA to the cellular process in which it participates. Recently, it has been shown that, by expressing artificial guide RNAs (derived from naturally occurring guide RNAs), modifications can also be specifically introduced into other RNAs, thus offering an opportunity to study RNAs in vivo. Here, we present strategies for constructing guide RNAs and manipulating RNA modifications in the nucleus.

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.

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

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

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.

Analysis of small nucleolar RNAs reveals unique genetic features in malaria parasites

Background

Ribosome biogenesis is an energy consuming and stringently controlled process that involves hundreds of trans-acting factors. Small nucleolar RNAs (snoRNAs), important components of ribosome biogenesis are non-coding guide RNAs involved in rRNA processing, nucleotide modifications like 2'-O-ribose methylation, pseudouridylation and possibly gene regulation. snoRNAs are ubiquitous and are diverse in their genomic organization, mechanism of transcription and process of maturation. In vertebrates, most snoRNAs are present in introns of protein coding genes and are processed by exonucleolytic cleavage, while in plants they are transcribed as polycistronic transcripts.

Results

This is a comprehensive analysis of malaria parasite snoRNA genes and proteins that have a role in ribosomal biogenesis. Computational and experimental approaches have been used to identify several box C/D snoRNAs from different species of Plasmodium and confirm their expression. Our analyses reveal that the gene for endoribonuclease Rnt1 is absent from Plasmodium falciparum genome, which indicates the existence of alternative pre-rRNA processing pathways. The structural features of box C/D snoRNAs are highly conserved in Plasmodium genus; however, unlike other organisms most parasite snoRNAs are present in single copy. The genomic localization of parasite snoRNAs shows mixed patterns of those observed in plants, yeast and vertebrates. We have localized parasite snoRNAs in untranslated regions (UTR) of mRNAs, and this is an unprecedented and novel genetic feature. Akin to mammalian snoRNAs, those in Plasmodium may also behave as mobile genetic elements.

Conclusion

This study provides a comprehensive overview on trans-acting genes involved in ribosome biogenesis and also a genetic insight into malaria parasite snoRNA genes.

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.

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.

A combined computational and experimental analysis of two families of snoRNA genes from Caenorhabditis elegans, revealing the expression and evolution pattern of snoRNAs in nematodes

Small nucleolar RNAs (snoRNAs) are an abundant group of noncoding RNAs mainly involved in the posttranscriptional modifications of rRNAs in eukaryotes. Prior to this study, only 28 snoRNA genes had been identified from Caenorhabditis elegans, indicating that most snoRNA genes are hidden in the worm genome, which represents a simple multicellular metazoan. In this study, a genome-wide analysis of the two major families of snoRNA genes in C. elegans was performed using the snoscan and snoGPS programs incorporating comparative genome analyses. Seventy gene variants, including 36 box C/D and 34 box H/ACA snoRNA genes, were identified, of which 50 are novel. Two families of snoRNAs showed a characteristic genomic organization. Notably, 6 box C/D snoRNA genes were located in the antisense orientation of introns. In contrast to insect and mammal, the distances between many intronic snoRNAs and 3' splice sites of introns were less than 50 nt in the worm, an unexpected finding as intron-encoded snoRNAs in C. elegans are supposed to be expressed in a splicing-dependent pathway. Interestingly, a canonical H/ACA snoRNA, PsiCeU5-48, was revealed to be partially homologous to small Cajal body-specific RNA (scaRNA) U85 and U89 in fly and human, indicating a possible evolutionary relationship between snoRNAs and scaRNAs.

Characteristics of the tomato nuclear genome as determined by sequencing undermethylated EcoRI digested fragments

A collection of 9,990 single-pass nuclear genomic sequences, corresponding to 5 Mb of tomato DNA, were obtained using methylation filtration (MF) strategy and reduced to 7,053 unique undermethylated genomic islands (UGIs) distributed as follows: (1) 59% non-coding sequences, (2) 28% coding sequences, (3) 12% transposons-96% of which are class I retroelements, and (4) 1% organellar sequences integrated into the nuclear genome over the past approximately 100 million years. A more detailed analysis of coding UGIs indicates that the unmethylated portion of tomato genes extends as far as 676 bp upstream and 766 bp downstream of coding regions with an average of 174 and 171 bp, respectively. Based on the analysis of the UGI copy distribution, the undermethylated portion of the tomato genome is determined to account for the majority of the unmethylated genes in the genome and is estimated to constitute 61+/-15 Mb of DNA (approximately 5% of the entire genome)--which is significantly less than the 220 Mb estimated for gene-rich euchromatic arms of the tomato genome. This result indicates that, while most genes reside in the euchromatin, a significant portion of euchromatin is methylated in the intergenic spacer regions. Implications of the results for sequencing the genome of tomato and other solanaceous species are discussed.

The K-loop, a general feature of the Pyrococcus C/D guide RNAs, is an RNA structural motif related to the K-turn

The C/D guide RNAs predicted from the genomic sequences of three species of Pyrococcus delineate a family of small non-coding archaeal RNAs involved in the methylation of rRNA and tRNA. The C/D guides assemble into ribonucleoprotein (RNP) that contains the methyltransferase. The protein L7Ae, a key structural component of the RNP, binds to a Kink-turn (K-turn) formed by the C/D motif. The K-turn is a structure that consists of two RNA stems separated by a short asymmetric loop with a characteristic sharp bend (kink) between the two stems. The majority of the pyrococcal C/D guides contain a short 3 nt-spacer between the C'/D' motifs. We show here that conserved terminal stem-loops formed by the C'/D' motif of the Pyrococcus C/D RNAs are also L7Ae-binding sites. These stem-loops are related to the K-turn by sequence and structure, but they consist of a single stem closed by a terminal loop. We have named this structure the K-loop. We show that conserved non-canonical base pairs in the stem of the K-loop are necessary for L7Ae binding. For the C/D guides with a 3 nt-spacer we show that the sequence and length is also important. The K-loop could improve the stability of the C/D guide RNAs in Pyrococcal species, which are extreme hyperthermophiles.

Maintaining a conserved methylation in plant and insect U2 snRNA through compensatory mutation by nucleotide insertion

The extensive post-transcriptional modification of U2 snRNA is required for spliceosome assembly and pre-mRNA splicing in vertebrates. However, the rare modification of U2 snRNA in yeast implies a different mechanism for regulating spliceosome biogenesis in single-celled eukaryotes. To understand the evolutionary pattern of U2 snRNA methylation, we determined for the first time, the 2'-O-methylations of U2 snRNA in Oryza sativa, Arabidopsis thaliana and Drosophila melanogaster, and revealed two methylations which are conserved in a crucial region of U2 snRNA in plants. Interestingly, one of the methylations, U2-Cm29 is also methylated in D. melanogaster, but not in vertebrates. According to the methylation of U2-C29, computational analysis of databases identified three canonical box C/D snoRNAs, named OsmgU2-29, AtmgU2-29 and DmmgU2-28, as small methylation guides of U2 snRNA from O. sativa, A. thaliana and D. melanogaster, respectively. Although very divergent in their sequence, the three snoRNAs exhibit in common an 11 nucleotide-long sequence complementarity to corresponding U2 snRNA, implying a functional constraint on the modification during evolution. Interestingly, a nucleotide is found to be inserted both in U2 snRNA and DmmgU2-28 and maintains a perfect match of duplex specifying the methylation of C28 in Drosophila U2 snRNA. This is the first time a new model is being provided for compensatory mutations between a small guide RNA and its target by nucleotide insertion, instead of the known nucleotide substitution. In contrast to small Cajal body-specific RNAs (scaRNAs), the snoRNAs are similar to the reported singlet guide RNAs and are known to localize in nucleolus.

Involvement of SRD2-mediated activation of snRNA transcription in the control of cell proliferation competence in Arabidopsis

The transcription machinery of small nuclear RNA (snRNA) genes has been investigated extensively in human cell-free systems, but its physiological function in vivo has not been addressed. This paper demonstrates the physiological role of an activator of snRNA transcription using a temperature-sensitive mutant of Arabidopsis thaliana, srd2. Phenotypic characteristics of the srd2 mutant suggest that the SRD2 gene participates in the control of competence in cell proliferation. The SRD2 gene encodes a nuclear protein that shares sequence similarity with the human SNAP50 protein, which is a subunit of SNAPc and is required for snRNA transcription in vitro. Our results, obtained from analysis of snRNA expression in the srd2 mutant, indicate that the SRD2 protein functions in the upregulation of transcription of snRNA genes, the promoters of which contain the upstream sequence element, to elevate cell proliferation competence.

Genome-wide analyses of two families of snoRNA genes from Drosophila melanogaster, demonstrating the extensive utilization of introns for coding of snoRNAs

Small nucleolar RNAs (snoRNAs) are an abundant group of noncoding RNAs mainly involved in the post-transcriptional modifications of rRNAs in eukaryotes. In this study, a large-scale genome-wide analysis of the two major families of snoRNA genes in the fruit fly Drosophila melanogaster has been performed using experimental and computational RNomics methods. Two hundred and twelve gene variants, encoding 56 box H/ACA and 63 box C/D snoRNAs, were identified, of which 57 novel snoRNAs have been reported for the first time. These snoRNAs were predicted to guide a total of 147 methylations and pseudouridylations on rRNAs and snRNAs, showing a more comprehensive pattern of rRNA modification in the fruit fly. With the exception of nine, all the snoRNAs identified to date in D. melanogaster are intron encoded. Remarkably, the genomic organization of the snoRNAs is characteristic of 8 dUhg genes and 17 intronic gene clusters, demonstrating that distinct organizations dominate the expression of the two families of snoRNAs in the fruit fly. Of the 267 introns in the host genes, more than half have been identified as host introns for coding of snoRNAs. In contrast to mammals, the variation in size of the host introns is mainly due to differences in the number of snoRNAs they contain. These results demonstrate the extensive utilization of introns for coding of snoRNAs in the host genes and shed light on further research of other noncoding RNA genes in the large introns of the Drosophila genome.

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.

Identification of 20 snoRNA-like RNAs from the primitive eukaryote, Giardia lamblia

From a specialized cDNA library of Giardia lamblia, 20 snoRNA-like RNAs, including 16 box C/D sRNAs and four box H/ACA sRNAs, were first identified. The sRNAs were predicted to guide a total of 11 2'-O-methylation and four pseudouridylation sites on the G. lamblia rRNAs, respectively. By using primer extension assay, seven methylation sites were precisely mapped in the G. lamblia 16S rRNA, despite its high GC content. All of the sRNA genes locate on the small intergenic regions of the G. lamblia genome and seem to be independently transcribed from their own promoters. Particularly, a cluster composed of GlsR17 and GlsR18 genes is transcribed as a dicistronic precursor, implying a mechanism of endonuclease cleavage for the maturation of the two sRNAs. The systematic identification of the sRNAs in G. lamblia has provided valuable information about the characteristics of the two major families of small guide RNAs in one of the most primitive eukaryotes and would contribute to the understanding of the evolution of small non-messenger RNA genes from prokaryotes to eukaryotes.

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.

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.

Identification of a novel box C/D snoRNA from mouse nucleolar cDNA library

By construction and screen of mouse nucleolar cDNA library, a novel mammalian small nucleolar RNAs (snoRNA) was identified. The novel snoRNA, 70 nt in length, displays structural features typical of C/D box snoRNA family. The snoRNA possesses an 11-nt-long rRNA antisense element and is predicted to guide the 2'-O-methylation of mouse 28S rRNA at G4043, a site unknown so far to be modified in vertebrates. The comparison of functional element of snoRNA guides among eukaryotes reveals that the novel snoRNA is a mammalian counterpart of yeast snR38 despite highly divergent sequence between them. Mouse and human snR38 and other cognates in distant vertebrates were positively detected with slight length variability. As expected, the rRNA ribose-methylation site predicted by mouse snR38 was precisely mapped by specific-primer extension assay. Furthermore, our analyses show that mouse and human snR38 gene have multiple variants and are nested in the introns of different host genes with unknown function. Thus, snR38 is a phylogenetically conserved methylation guide but exhibits different genomic organization in eukaryotes.

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.