Plant snoRNA database

RNA pseudouridine modification in plants

Pseudouridine is one of the well-known chemical modifications in various RNA species. Current advances to detect pseudouridine show that the pseudouridine landscape is dynamic and affects multiple cellular processes. Although our understanding of this post-transcriptional modification mainly depends on yeast and human models, the recent findings provide strong evidence for the critical role of pseudouridine in plants. Here, we review the current knowledge of pseudouridine in plant RNAs, including its synthesis, degradation, regulatory mechanisms, and functions. Moreover, we propose future areas of research on pseudouridine modification in plants.

Arabidopsis retains vertebrate-type telomerase accessory proteins via a plant-specific assembly

The recent discovery of the bona-fide telomerase RNA (TR) from plants reveals conserved and unique secondary structure elements and the opportunity for new insight into the telomerase RNP. Here we examine how two highly conserved proteins previously implicated in Arabidopsis telomere maintenance, AtPOT1a and AtNAP57 (dyskerin), engage plant telomerase. We report that AtPOT1a associates with Arabidopsis telomerase via interaction with TERT. While loss of AtPOT1a does not impact AtTR stability, the templating domain is more accessible in pot1a mutants, supporting the conclusion that AtPOT1a stimulates telomerase activity but does not facilitate telomerase RNP assembly. We also show, that despite the absence of a canonical H/ACA binding motif within AtTR, dyskerin binds AtTR with high affinity and specificity in vitro via a plant specific three-way junction (TWJ). A core element of the TWJ is the P1a stem, which unites the 5′ and 3′ ends of AtTR. P1a is required for dyskerin-mediated stimulation of telomerase repeat addition processivity in vitro, and for AtTR accumulation and telomerase activity in vivo. The deployment of vertebrate-like accessory proteins and unique RNA structural elements by Arabidopsis telomerase provides a new platform for exploring telomerase biogenesis and evolution.

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.

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.

Transcriptome-wide analysis of pseudouridylation of mRNA and non-coding RNAs in Arabidopsis

Pseudouridine (Ψ) is widely distributed in mRNA and various non-coding RNAs in yeast and mammals, and the specificity of its distribution has been determined. However, knowledge about Ψs in the RNAs of plants, particularly in mRNA, is lacking. In this study, we performed genome-wide pseudouridine-sequencing in Arabidopsis and for the first time identified hundreds of Ψ sites in mRNA and multiple Ψ sites in non-coding RNAs. Many predicted and novel Ψ sites in rRNA and tRNA were detected. mRNA was extensively pseudouridylated, but with Ψs being under-represented in 3'-untranslated regions and enriched at position 1 of triple codons. The phenylalanine codon UUC was the most frequently pseudouridylated site. Some Ψs present in chloroplast 23S, 16S, and 4.5S rRNAs in wild-type Col-0 were absent in plants with a mutation of SVR1 (Suppressor of variegation 1), a chloroplast pseudouridine synthase gene. Many plastid ribosomal proteins and photosynthesis-related proteins were significantly reduced in svr1 relative to the wild-type, indicating the roles of SVR1 in chloroplast protein biosynthesis in Arabidopsis. Our results provide new insights into the occurrence of pseudouridine in Arabidopsis RNAs and the biological functions of SVR1, and will pave the way for further exploiting the mechanisms underlying Ψ modifications in controlling gene expression and protein biosynthesis in plants.

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.

Identification, Prediction and Data Analysis of Noncoding RNAs: A Review

BACKGROUND

Noncoding RNAs (ncRNAs) which play an important role in various cellular processes are important in medicine as well as in drug design strategies. Different studies have shown that ncRNAs are dis-regulated in cancer cells and play an important role in human tumorigenesis. Therefore, it is important to identify and predict such molecules by experimental and computational methods, respectively. However, to avoid expensive experimental methods, computational algorithms have been developed for accurately and fast prediction of ncRNAs.

OBJECTIVE

The aim of this review was to introduce the experimental and computational methods to identify and predict ncRNAs structure. Also, we explained the ncRNA's roles in cellular processes and drugs design, briefly.

METHOD

In this survey, we will introduce ncRNAs and their roles in biological and medicinal processes. Then, some important laboratory techniques will be studied to identify ncRNAs. Finally, the state-of-the-art models and algorithms will be introduced along with important tools and databases.

RESULTS

The results showed that the integration of experimental and computational approaches improves to identify ncRNAs. Moreover, the high accurate databases, algorithms and tools were compared to predict the ncRNAs.

CONCLUSION

ncRNAs prediction is an exciting research field, but there are different difficulties. It requires accurate and reliable algorithms and tools. Also, it should be mentioned that computational costs of such algorithm including running time and usage memory are very important. Finally, some suggestions were presented to improve computational methods of ncRNAs gene and structural prediction.

Diff isomiRs: Large-scale detection of differential isomiRs for understanding non-coding regulated stress omics in plants

Plants have an amazing ability to cope with wide variety of stresses by regulating the expression of genes and thus by altering the physiological status. In the past few years, canonical microRNA variants (isomiRs) have been shown to play pivotal roles by acting as regulators of the transcriptional machinery. In the present research, we present Diff isomiRs, a web-based exploratory repository of differential isomiRs across 16 sequenced plant species representing a total of 433 datasets across 21 different stresses and 158 experimental states. Diff isomiRs provides the high-throughput detection of differential isomiRs using mapping-based and model-based differential analysis revealing a total of 16,157 and 2,028 differential isomiRs, respectively. Easy-to-use and web-based exploration of differential isomiRs provides several features such as browsing of the differential isomiRs according to stress or species, as well as association of the differential isomiRs to targets and plant endogenous target mimics (PeTMs). Diff isomiRs also provides the relationship between the canonical miRNAs, isomiRs and the miRNA-target interactions. This is the first web-based large-scale repository for browsing differential isomiRs and will facilitate better understanding of the regulatory role of the isomiRs with respect to the canonical microRNAs. Diff isomiRs can be accessed at: www.mcr.org.in/diffisomirs.

Noncoding RNAs Databases: Current Status and Trends

One of the most important resources for researchers of noncoding RNAs is the information available in public databases spread over the internet. However, the effective exploration of this data can represent a daunting task, given the large amount of databases available and the variety of stored data. This chapter describes a classification of databases based on information source, type of RNA, source organisms, data formats, and the mechanisms for information retrieval, detailing the relevance of each of these classifications and its usability by researchers. This classification is used to update a 2012 review, indexing now more than 229 public databases. This review will include an assessment of the new trends for ncRNA research based on the information that is being offered by the databases. Additionally, we will expand the previous analysis focusing on the usability and application of these databases in pathogen and disease research. Finally, this chapter will analyze how currently available database schemas can help the development of new and improved web resources.

A comprehensive review of web-based resources of non-coding RNAs for plant science research

Non-coding RNAs (ncRNAs) are transcribed from genome but not translated into proteins. Many ncRNAs are key regulators of plants growth and development, metabolism and stress tolerance. In order to make the web-based ncRNA resources for plant science research be more easily accessible and understandable, we made a comprehensive review for 83 web-based resources of three types, including genome databases containing ncRNA data, microRNA (miRNA) databases and long non-coding RNA (lncRNA) databases. To facilitate effective usage of these resources, we also suggested some preferred resources of miRNAs and lncRNAs for performing meaningful analysis.

The Medicago truncatula Genome: Genomic Data Availability

Medicago truncatula emerged in 1990 as a model for legumes, comprising the third largest land plant family. Most legumes form symbiotic nitrogen-fixing root nodules with compatible soil bacteria and thus are important contributors to the global nitrogen cycle and sustainable agriculture. Legumes and legume products are important sources for human and animal protein as well as for edible and industrial oils. In the years since M. truncatula was chosen as a legume model, many genetic, genomic, and molecular resources have become available, including reference quality genome sequences for two widely used genotypes. Accessibility of genomic data is important for many different types of studies with M. truncatula as well as for research involving crop and forage legumes. In this chapter, we discuss strategies to obtain archived M. truncatula genomic data originally deposited into custom databases that are no longer maintained but are now accessible in general databases. We also review key current genomic databases that are specific to M. truncatula as well as those that contain M. truncatula data in addition to data from other plants.

Small RNA differential expression and regulation in Tuxpeño maize embryogenic callus induction and establishment

Somatic embryogenesis represents an alternative developmental process used to achieve genetic transformation and to approach key questions in maize development. It is known that embryogenic callus induction and plant regeneration are accompanied by microRNA expression changes. However, small RNA (sRNA) populations have not been explored during the proliferative callus subculture establishment and their impact on maintaining the dedifferentiated status and embryogenic potential is far from being completely understood. Here we globally tested the sRNA populations in explants (immature embryos), induced and established maize embryogenic callus from the Mexican cultivar VS-535, Tuxpeño landrace. We detected readjustments in 24 nt and 21-22 nt sRNAs during the embryogenic callus (EC) establishment and maintenance. A follow up on specific microRNAs (miRNAs) indicated that miRNAs related to stress response substantially increase upon the callus proliferation establishment, correlating with a reduction in some of their target levels. On the other hand, while 24 nt-long heterochromatic small interfering RNAs (hc-siRNAs) derived from transposable retroelements transiently decreased in abundance during the EC establishment, a population of 22 nt-hc-siRNAs increased. This was accompanied by reduction in transposon expression in the established callus subcultures. We conclude that stress- and development-related miRNAs are highly expressed upon maize EC callus induction and during maintenance of the subcultures, while miRNAs involved in hormone response only transiently increase during induction. In addition, the establishment of a proliferative status in embryogenic callus is accompanied by important readjustments in hc-siRNAs mapping to long tandem repeat (LTR) retrotransposons, and their expression regulation.

An Arabidopsis divergent pumilio protein, APUM24, is essential for embryogenesis and required for faithful pre-rRNA processing

Pumilio RNA-binding proteins are largely involved in mRNA degradation and translation repression. However, a few evolutionarily divergent Pumilios are also responsible for proper pre-rRNA processing in human and yeast. Here, we describe an essential Arabidopsis nucleolar Pumilio, APUM24, that is expressed in tissues undergoing rapid proliferation and cell division. A T-DNA insertion for APUM24 did not affect the male and female gametogenesis, but instead resulted in a negative female gametophytic effect on zygotic cell division immediately after fertilization. Additionally, the mutant embryos displayed defects in cell patterning from pro-embryo through globular stages. The mutant embryos were marked by altered auxin maxima, which were substantiated by the mislocalization of PIN1 and PIN7 transporters in the defective embryos. Homozygous apum24 callus accumulates rRNA processing intermediates, including uridylated and adenylated 5.8S and 25S rRNA precursors. An RNA-protein interaction assay showed that the histidine-tagged recombinant APUM24 binds RNAin vitro with no apparent specificity. Overall, our results demonstrated that APUM24 is required for rRNA processing and early embryogenesis in Arabidopsis.

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.

Plant-Specific Features of Ribosome Biogenesis

The biogenesis of eukaryotic ribosomes is a fundamental process involving hundreds of ribosome biogenesis factors (RBFs) in three compartments of the cell, namely the nucleolus, nucleus, and cytoplasm. Many RBFs are involved in the processing of the primary ribosomal (r)RNA transcript, in which three of the four rRNAs are imbedded. While pre-rRNA processing is well described for yeast and mammals, a detailed processing scheme for plants is lacking. Here, we discuss the emerging scheme of pre-rRNA processing in Arabidopsis thaliana in comparison to other eukaryotes, with a focus on plant characteristics. In addition, we highlight the impact of the ribosome and its biogenesis on developmental processes because common phenotypes can be observed for ribosomal protein and RBF mutants.

Identification of novel small ncRNAs in pollen of tomato

BACKGROUND

The unprecedented role of sncRNAs in the regulation of pollen biogenesis on both transcriptional and epigenetic levels has been experimentally proven. However, little is known about their global regulation, especially under stress conditions. We used tomato pollen in order to identify pollen stage-specific sncRNAs and their target mRNAs. We further deployed elevated temperatures to discern stress responsive sncRNAs. For this purpose high throughput sncRNA-sequencing as well as Massive Analysis of cDNA Ends (MACE) were performed for three-replicated sncRNAs libraries derived from tomato tetrad, post-meiotic, and mature pollen under control and heat stress conditions.

RESULTS

Using the omiRas analysis pipeline we identified known and predicted novel miRNAs as well as sncRNAs from other classes, responsive or not to heat. Differential expression analysis revealed that post-meiotic and mature pollen react most strongly by regulation of the expression of coding and non-coding genomic regions in response to heat. To gain insight to the function of these miRNAs, we predicted targets and annotated them to Gene Ontology terms. This approach revealed that most of them belong to protein binding, transcription, and Serine/Threonine kinase activity GO categories. Beside miRNAs, we observed differential expression of both tRNAs and snoRNAs in tetrad, post-meiotic, and mature pollen when comparing normal and heat stress conditions.

CONCLUSIONS

Thus, we describe a global spectrum of sncRNAs expressed in pollen as well as unveiled those which are regulated at specific time-points during pollen biogenesis. We integrated the small RNAs into the regulatory network of tomato heat stress response in pollen.

Cell dedifferentiation and organogenesis in vitro require more snRNA than does seedling development in Arabidopsis thaliana

Small nuclear RNA (snRNA) is a class of non-coding RNAs that processes pre-mRNA and rRNA. Transcription of abundant snRNA species is regulated by the snRNA activating protein complex (SNAPc), which is conserved among multicellular organisms including plants. SRD2, a putative subunit of SNAPc in Arabidopsis thaliana, is essential for development, and the point mutation srd2-1 causes severe defects in hypocotyl dedifferentiation and de novo meristem formation. Based on phenotypic analysis of srd2-1 mutant plants, we previously proposed that snRNA content is a limiting factor in dedifferentiation in plant cells. Here, we performed functional complementation analysis of srd2-1 using transgenic srd2-1 Arabidopsis plants harboring SRD2 homologs from Populus trichocarpa (poplar), Nicotiana tabacum (tobacco), Oryza sativa (rice), the moss Physcomitrella patens, and Homo sapiens (human) under the control of the Arabidopsis SRD2 promoter. Only rice SRD2 suppressed the faulty tissue culture responses of srd2-1, and restore the snRNA levels; however, interestingly, all SRD2 homologs except poplar SRD2 rescued the srd2-1 defects in seedling development. These findings demonstrated that cell dedifferentiation and organogenesis induced during tissue culture require higher snRNA levels than does seedling development.

Identification and Expression Analysis of Ribosome Biogenesis Factor Co-orthologs in Solanum lycopersicum

Ribosome biogenesis involves a large inventory of proteinaceous and RNA cofactors. More than 250 ribosome biogenesis factors (RBFs) have been described in yeast. These factors are involved in multiple aspects like rRNA processing, folding, and modification as well as in ribosomal protein (RP) assembly. Considering the importance of RBFs for particular developmental processes, we examined the complexity of RBF and RP (co-)orthologs by bioinformatic assignment in 14 different plant species and expression profiling in the model crop Solanum lycopersicum. Assigning (co-)orthologs to each RBF revealed that at least 25% of all predicted RBFs are encoded by more than one gene. At first we realized that the occurrence of multiple RBF co-orthologs is not globally correlated to the existence of multiple RP co-orthologs. The transcript abundance of genes coding for predicted RBFs and RPs in leaves and anthers of S. lycopersicum was determined by next generation sequencing (NGS). In combination with existing expression profiles, we can conclude that co-orthologs of RBFs by large account for a preferential function in different tissue or at distinct developmental stages. This notion is supported by the differential expression of selected RBFs during male gametophyte development. In addition, co-regulated clusters of RBF and RP coding genes have been observed. The relevance of these results is discussed.

Characterization of stress-responsive lncRNAs in Arabidopsis thaliana by integrating expression, epigenetic and structural features

Recently, in addition to poly(A)+ long non-coding RNAs (lncRNAs), many lncRNAs without poly(A) tails, have been characterized in mammals. However, the non-polyA lncRNAs and their conserved motifs, especially those associated with environmental stresses, have not been fully investigated in plant genomes. We performed poly(A)- RNA-seq for seedlings of Arabidopsis thaliana under four stress conditions, and predicted lncRNA transcripts. We classified the lncRNAs into three confidence levels according to their expression patterns, epigenetic signatures and RNA secondary structures. Then, we further classified the lncRNAs to poly(A)+ and poly(A)- transcripts. Compared with poly(A)+ lncRNAs and coding genes, we found that poly(A)- lncRNAs tend to have shorter transcripts and lower expression levels, and they show significant expression specificity in response to stresses. In addition, their differential expression is significantly enriched in drought condition and depleted in heat condition. Overall, we identified 245 poly(A)+ and 58 poly(A)- lncRNAs that are differentially expressed under various stress stimuli. The differential expression was validated by qRT-PCR, and the signaling pathways involved were supported by specific binding of transcription factors (TFs), phytochrome-interacting factor 4 (PIF4) and PIF5. Moreover, we found many conserved sequence and structural motifs of lncRNAs from different functional groups (e.g. a UUC motif responding to salt and a AU-rich stem-loop responding to cold), indicated that the conserved elements might be responsible for the stress-responsive functions of lncRNAs.

Superinfection exclusion by Citrus tristeza virus does not correlate with the production of viral small RNAs

Superinfection exclusion (SIE), a phenomenon in which a preexisting viral infection prevents a secondary infection with the same or closely related virus, has been described for different viruses, including important pathogens of humans, animals, and plants. Several mechanisms acting at various stages of the viral life cycle have been proposed to explain SIE. Most cases of SIE in plant virus systems were attributed to induction of RNA silencing, a host defense mechanism that is mediated by small RNAs. Here we show that SIE by Citrus tristeza virus (CTV) does not correlate with the production of viral small interfering RNAs (siRNAs). CTV variants, which differed in the SIE ability, had similar siRNAs profiles. Along with our previous observations that the exclusion phenomenon requires a specific viral protein, p33, the new data suggest that SIE by CTV is highly complex and appears to use different mechanisms than those proposed for other viruses.

Ribogenomics: the science and knowledge of RNA

Ribonucleic acid (RNA) deserves not only a dedicated field of biological research — a discipline or branch of knowledge — but also explicit definitions of its roles in cellular processes and molecular mechanisms. Ribogenomics is to study the biology of cellular RNAs, including their origin, biogenesis, structure and function. On the informational track, messenger RNAs (mRNAs) are the major component of ribogenomes, which encode proteins and serve as one of the four major components of the translation machinery and whose expression is regulated at multiple levels by other operational RNAs. On the operational track, there are several diverse types of RNAs — their length distribution is perhaps the most simplistic stratification — involving in major cellular activities, such as chromosomal structure and organization, DNA replication and repair, transcriptional/post-transcriptional regulation, RNA processing and routing, translation and cellular energy/metabolism regulation. An all-out effort exceeding the magnitude of the Human Genome Project is of essence to construct just mammalian transcriptomes in multiple contexts including embryonic development, circadian and seasonal rhythms, defined life-span stages, pathological conditions and anatomy-driven tissue/organ/cell types.

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

The genome of the recently domesticated crop plant sugar beet (Beta vulgaris)

Sugar beet (Beta vulgaris ssp. vulgaris) is an important crop of temperate climates which provides nearly 30% of the world's annual sugar production and is a source for bioethanol and animal feed. The species belongs to the order of Caryophylalles, is diploid with 2n = 18 chromosomes, has an estimated genome size of 714-758 megabases and shares an ancient genome triplication with other eudicot plants. Leafy beets have been cultivated since Roman times, but sugar beet is one of the most recently domesticated crops. It arose in the late eighteenth century when lines accumulating sugar in the storage root were selected from crosses made with chard and fodder beet. Here we present a reference genome sequence for sugar beet as the first non-rosid, non-asterid eudicot genome, advancing comparative genomics and phylogenetic reconstructions. The genome sequence comprises 567 megabases, of which 85% could be assigned to chromosomes. The assembly covers a large proportion of the repetitive sequence content that was estimated to be 63%. We predicted 27,421 protein-coding genes supported by transcript data and annotated them on the basis of sequence homology. Phylogenetic analyses provided evidence for the separation of Caryophyllales before the split of asterids and rosids, and revealed lineage-specific gene family expansions and losses. We sequenced spinach (Spinacia oleracea), another Caryophyllales species, and validated features that separate this clade from rosids and asterids. Intraspecific genomic variation was analysed based on the genome sequences of sea beet (Beta vulgaris ssp. maritima; progenitor of all beet crops) and four additional sugar beet accessions. We identified seven million variant positions in the reference genome, and also large regions of low variability, indicating artificial selection. The sugar beet genome sequence enables the identification of genes affecting agronomically relevant traits, supports molecular breeding and maximizes the plant's potential in energy biotechnology.

snOPY: a small nucleolar RNA orthological gene database

Background

Small nucleolar RNAs (snoRNAs) are a class of non-coding RNAs that guide the modification of specific nucleotides in ribosomal RNAs (rRNAs) and small nuclear RNAs (snRNAs). Although most non-coding RNAs undergo post-transcriptional modifications prior to maturation, the functional significance of these modifications remains unknown. Here, we introduce the snoRNA orthological gene database (snOPY) as a tool for studying RNA modifications.

Findings

snOPY provides comprehensive information about snoRNAs, snoRNA gene loci, and target RNAs. It also contains data for orthologues from various species, which enables users to analyze the evolution of snoRNA genes. In total, 13,770 snoRNA genes, 10,345 snoRNA gene loci, and 133 target RNAs have been registered. Users can search and access the data efficiently using a simple web interface with a series of internal links. snOPY is freely available on the web at http://snoopy.med.miyazaki-u.ac.jp.

Conclusions

snOPY is the database that provides information about the small nucleolar RNAs and their orthologues. It will help users to study RNA modifications and snoRNA gene evolution.

Emerging landscape of non-coding RNAs in oral health and disease

The world of non-coding RNAs has only recently started being discovered. For the past 40 years, coding genes, mRNA, and proteins have been the center of cellular and molecular biology, and pathologic alterations were attributed to either the aberration of gene sequence or altered promoter activity. It was only after the completion of the human genome sequence that the scientific community started seriously wondering why only a very small portion of the genome corresponded to protein-coding genes. New technologies such as the whole-genome and whole-transcriptome sequencing demonstrated that at least 90% of the genome is actively transcribed. The identification and cataloguing of multiple kinds of non-coding RNA (ncRNA) have exponentially increased, and it is now widely accepted that ncRNAs play major biological roles in cellular physiology, development, metabolism, and are also implicated in a variety of diseases. The aim of this review is to describe the two major classes (long and short forms) of non-coding RNAs and describe their subclasses in terms of function and their relevance and potential in oral diseases.

Citrus tristeza virus: Evolution of Complex and Varied Genotypic Groups

Amongst the Closteroviridae, Citrus tristeza virus (CTV) is almost unique in possessing a number of distinct and characterized strains, isolates of which produce a wide range of phenotype combinations among its different hosts. There is little understanding to connect genotypes to phenotypes, and to complicate matters more, these genotypes are found throughout the world as members of mixed populations within a single host plant. There is essentially no understanding of how combinations of genotypes affect symptom expression and disease severity. We know little about the evolution of the genotypes that have been characterized to date, little about the biological role of their diversity and particularly, about the effects of recombination. Additionally, genotype grouping has not been standardized. In this study we utilized an extensive array of CTV genomic information to classify the major genotypes, and to determine the major evolutionary processes that led to their formation and subsequent retention. Our analyses suggest that three major processes act on these genotypes: (1) ancestral diversification of the major CTV lineages, followed by (2) conservation and co-evolution of the major functional domains within, though not between CTV genotypes, and (3) extensive recombination between lineages that have given rise to new genotypes that have subsequently been retained within the global population. The effects of genotype diversity and host-interaction are discussed, as is a proposal for standardizing the classification of existing and novel CTV genotypes.

Small RNA and transcriptome deep sequencing proffers insight into floral gene regulation in Rosa cultivars

Background

Roses (Rosa sp.), which belong to the family Rosaceae, are the most economically important ornamental plants—making up 30% of the floriculture market. However, given high demand for roses, rose breeding programs are limited in molecular resources which can greatly enhance and speed breeding efforts. A better understanding of important genes that contribute to important floral development and desired phenotypes will lead to improved rose cultivars. For this study, we analyzed rose miRNAs and the rose flower transcriptome in order to generate a database to expound upon current knowledge regarding regulation of important floral characteristics. A rose genetic database will enable comprehensive analysis of gene expression and regulation via miRNA among different Rosa cultivars.

Results

We produced more than 0.5 million reads from expressed sequences, totalling more than 110 million bp. From these, we generated 35,657, 31,434, 34,725, and 39,722 flower unigenes from Rosa hybrid: ‘Vital’, ‘Maroussia’, and ‘Sympathy’ and Rosa rugosa Thunb. , respectively. The unigenes were assigned functional annotations, domains, metabolic pathways, Gene Ontology (GO) terms, Plant Ontology (PO) terms, and MIPS Functional Catalogue (FunCat) terms. Rose flower transcripts were compared with genes from whole genome sequences of Rosaceae members (apple, strawberry, and peach) and grape. We also produced approximately 40 million small RNA reads from flower tissue for Rosa, representing 267 unique miRNA tags. Among identified miRNAs, 25 of them were novel and 242 of them were conserved miRNAs. Statistical analyses of miRNA profiles revealed both shared and species-specific miRNAs, which presumably effect flower development and phenotypes.

Conclusions

In this study, we constructed a Rose miRNA and transcriptome database, and we analyzed the miRNAs and transcriptome generated from the flower tissues of four Rosa cultivars. The database provides a comprehensive genetic resource which can be used to better understand rose flower development and to identify candidate genes for important phenotypes.

Comparative genomics of eukaryotic small nucleolar RNAs reveals deep evolutionary ancestry amidst ongoing intragenomic mobility

BACKGROUND

Small nucleolar (sno)RNAs are required for posttranscriptional processing and modification of ribosomal, spliceosomal and messenger RNAs. Their presence in both eukaryotes and archaea indicates that snoRNAs are evolutionarily ancient. The location of some snoRNAs within the introns of ribosomal protein genes has been suggested to belie an RNA world origin, with the exons of the earliest protein-coding genes having evolved around snoRNAs after the advent of templated protein synthesis. Alternatively, this intronic location may reflect more recent selection for coexpression of snoRNAs and ribosomal components, ensuring rRNA modification by snoRNAs during ribosome synthesis. To gain insight into the evolutionary origins of this genetic organization, we examined the antiquity of snoRNA families and the stability of their genomic location across 44 eukaryote genomes.

RESULTS

We report that dozens of snoRNA families are traceable to the Last Eukaryotic Common Ancestor (LECA), but find only weak similarities between the oldest eukaryotic snoRNAs and archaeal snoRNA-like genes. Moreover, many of these LECA snoRNAs are located within the introns of host genes independently traceable to the LECA. Comparative genomic analyses reveal the intronic location of LECA snoRNAs is not ancestral however, suggesting the pattern we observe is the result of ongoing intragenomic mobility. Analysis of human transcriptome data indicates that the primary requirement for hosting intronic snoRNAs is a broad expression profile. Consistent with ongoing mobility across broadly-expressed genes, we report a case of recent migration of a non-LECA snoRNA from the intron of a ubiquitously expressed non-LECA host gene into the introns of two LECA genes during the evolution of primates.

CONCLUSIONS

Our analyses show that snoRNAs were a well-established family of RNAs at the time when eukaryotes began to diversify. While many are intronic, this association is not evolutionarily stable across the eukaryote tree; ongoing intragenomic mobility has erased signal of their ancestral gene organization, and neither introns-first nor evolved co-expression adequately explain our results. We therefore present a third model - constrained drift - whereby individual snoRNAs are intragenomically mobile and may occupy any genomic location from which expression satisfies phenotype.

Analysis of the melon (Cucumis melo) small RNAome by high-throughput pyrosequencing

BACKGROUND

Melon (Cucumis melo L.) is a commercially important fruit crop that is cultivated worldwide. The melon research community has recently benefited from the determination of a complete draft genome sequence and the development of associated genomic tools, which have allowed us to focus on small RNAs (sRNAs). These are short, non-coding RNAs 21-24 nucleotides in length with diverse physiological roles. In plants, they regulate gene expression and heterochromatin assembly, and control protection against virus infection. Much remains to be learned about the role of sRNAs in melon.

RESULTS

We constructed 10 sRNA libraries from two stages of developing ovaries, fruits and photosynthetic cotyledons infected with viruses, and carried out high-throughput pyrosequencing. We catalogued and analysed the melon sRNAs, resulting in the identification of 26 known miRNA families (many conserved with other species), the prediction of 84 melon-specific miRNA candidates, the identification of trans-acting siRNAs, and the identification of chloroplast, mitochondrion and transposon-derived sRNAs. In silico analysis revealed more than 400 potential targets for the conserved and novel miRNAs.

CONCLUSION

We have discovered and analysed a large number of conserved and melon-specific sRNAs, including miRNAs and their potential target genes. This provides insight into the composition and function of the melon small RNAome, and paves the way towards an understanding of sRNA-mediated processes that regulate melon fruit development and melon-virus interactions.

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.

ARTADE2DB: improved statistical inferences for Arabidopsis gene functions and structure predictions by dynamic structure-based dynamic expression (DSDE) analyses

Recent advances in technologies for observing high-resolution genomic activities, such as whole-genome tiling arrays and high-throughput sequencers, provide detailed information for understanding genome functions. However, the functions of 50% of known Arabidopsis thaliana genes remain unknown or are annotated only on the basis of static analyses such as protein motifs or similarities. In this paper, we describe dynamic structure-based dynamic expression (DSDE) analysis, which sequentially predicts both structural and functional features of transcripts. We show that DSDE analysis inferred gene functions 12% more precisely than static structure-based dynamic expression (SSDE) analysis or conventional co-expression analysis based on previously determined gene structures of A. thaliana. This result suggests that more precise structural information than the fixed conventional annotated structures is crucial for co-expression analysis in systems biology of transcriptional regulation and dynamics. Our DSDE method, ARabidopsis Tiling-Array-based Detection of Exons version 2 and over-representation analysis (ARTADE2-ORA), precisely predicts each gene structure by combining two statistical analyses: a probe-wise co-expression analysis of multiple transcriptome measurements and a Markov model analysis of genome sequences. ARTADE2-ORA successfully identified the true functions of about 90% of functionally annotated genes, inferred the functions of 98% of functionally unknown genes and predicted 1,489 new gene structures and functions. We developed a database ARTADE2DB that integrates not only the information predicted by ARTADE2-ORA but also annotations and other functional information, such as phenotypes and literature citations, and is expected to contribute to the study of the functional genomics of A. thaliana. URL: http://artade.org.

Hybrid activation methods for elucidating nucleic acid modifications

Hybrid tandem mass spectrometry (MS/MS) techniques combining electron transfer (ET) and collision activated dissociation (CAD), infrared multiphoton dissociation (IRMPD), or ultraviolet photodissociation (UVPD) were implemented and evaluated for the characterization of a series of oligonucleotides and oligoribonucleotides, including both native single strands and single strands containing platinated, phosphorothioated, and 2'-O-methylated modification sites. ET-IRMPD and ET-UVPD of oligodeoxynucleotides and oligoribonucleotides resulted in rich fragmentation with respect to production of w, a, z, and d ions for DNA, and c, y, w, a-B, d, and z ions for RNA, with many product ions retaining the modification and thus allowing site specific identification. ET-IRMPD caused more extensive secondary dissociation of the ions, in addition to a broader distribution of detectable sequence ions attributed to using a lower mass cutoff. ET-UVPD promoted higher energy fragmentation pathways and created the most diverse MS/MS spectra. The numerous products generated by the hybrid MS/MS techniques resulted in specific and extensive backbone cleavages which allowed the modification sites of multiply modified oligonucleotides to be elucidated.

Rfam: Wikipedia, clans and the "decimal" release

The Rfam database aims to catalogue non-coding RNAs through the use of sequence alignments and statistical profile models known as covariance models. In this contribution, we discuss the pros and cons of using the online encyclopedia, Wikipedia, as a source of community-derived annotation. We discuss the addition of groupings of related RNA families into clans and new developments to the website. Rfam is available on the Web at http://rfam.sanger.ac.uk.

starBase: a database for exploring microRNA-mRNA interaction maps from Argonaute CLIP-Seq and Degradome-Seq data

MicroRNAs (miRNAs) represent an important class of small non-coding RNAs (sRNAs) that regulate gene expression by targeting messenger RNAs. However, assigning miRNAs to their regulatory target genes remains technically challenging. Recently, high-throughput CLIP-Seq and degradome sequencing (Degradome-Seq) methods have been applied to identify the sites of Argonaute interaction and miRNA cleavage sites, respectively. In this study, we introduce a novel database, starBase (sRNA target Base), which we have developed to facilitate the comprehensive exploration of miRNA–target interaction maps from CLIP-Seq and Degradome-Seq data. The current version includes high-throughput sequencing data generated from 21 CLIP-Seq and 10 Degradome-Seq experiments from six organisms. By analyzing millions of mapped CLIP-Seq and Degradome-Seq reads, we identified ∼1 million Ago-binding clusters and ∼2 million cleaved target clusters in animals and plants, respectively. Analyses of these clusters, and of target sites predicted by 6 miRNA target prediction programs, resulted in our identification of approximately 400 000 and approximately 66 000 miRNA-target regulatory relationships from CLIP-Seq and Degradome-Seq data, respectively. Furthermore, two web servers were provided to discover novel miRNA target sites from CLIP-Seq and Degradome-Seq data. Our web implementation supports diverse query types and exploration of common targets, gene ontologies and pathways. The starBase is available at http://starbase.sysu.edu.cn/.

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.

The small nucleolar ribonucleoprotein (snoRNP) database

Small nucleolar ribonucleoproteins (snoRNPs) are widely studied and characterized as guide RNAs for sequence-specific 2'-O-ribose methylation and psuedouridylation of ribosomal RNAs. In addition, snoRNAs have also been shown to interact with some tRNAs and direct alternative splicing in mRNA biogenesis. Recent advances in bioinformatics have resulted in new algorithms able to rapidly identify noncoding RNAs generally and snoRNAs specifically in genomic and metagenomic sequences, resulting in a rapid increase in the number and diversity of identified snoRNA sequences. The snoRNP database is a web-based collection of snoRNA and snoRNA-associated protein sequences from a wide range of species. The database currently contains 8994 snoRNA sequences from Bacteria, Archaea, and Eukaryotes and 589 snoRNA-associated protein sequences. The snoRNP database can be found at: http://evolveathome.com/snoRNA/snoRNA.php.

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.

Nucleic acid sequence and structure databases

This chapter gives an overview of the most commonly used biological databases of nucleic acid sequences and their structures. We cover general sequence databases, databases for specific DNA features, noncoding RNA sequences, and RNA secondary and tertiary structures.

deepBase: a database for deeply annotating and mining deep sequencing data

Advances in high-throughput next-generation sequencing technology have reshaped the transcriptomic research landscape. However, exploration of these massive data remains a daunting challenge. In this study, we describe a novel database, deepBase, which we have developed to facilitate the comprehensive annotation and discovery of small RNAs from transcriptomic data. The current release of deepBase contains deep sequencing data from 185 small RNA libraries from diverse tissues and cell lines of seven organisms: human, mouse, chicken, Ciona intestinalis, Drosophila melanogaster, Caenhorhabditis elegans and Arabidopsis thaliana. By analyzing ∼14.6 million unique reads that perfectly mapped to more than 284 million genomic loci, we annotated and identified ∼380 000 unique ncRNA-associated small RNAs (nasRNAs), ∼1.5 million unique promoter-associated small RNAs (pasRNAs), ∼4.0 million unique exon-associated small RNAs (easRNAs) and ∼6 million unique repeat-associated small RNAs (rasRNAs). Furthermore, 2038 miRNA and 1889 snoRNA candidates were predicted by miRDeep and snoSeeker. All of the mapped reads can be grouped into about 1.2 million RNA clusters. For the purpose of comparative analysis, deepBase provides an integrative, interactive and versatile display. A convenient search option, related publications and other useful information are also provided for further investigation. deepBase is available at: http://deepbase.sysu.edu.cn/.

Loss of rRNA modifications in the decoding center of the ribosome impairs translation and strongly delays pre-rRNA processing

The ribosome decoding center is rich in modified rRNA nucleotides and little is known about their effects. Here, we examine the consequences of systematically deleting eight pseudouridine and 2'-O-methylation modifications in the yeast decoding center. Loss of most modifications individually has no apparent effect on cell growth. However, deletions of 2-3 modifications in the A- and P-site regions can cause (1) reduced growth rates (approximately 15%-50% slower); (2) reduced amino acid incorporation rates (14%-24% slower); and (3) a significant deficiency in free small subunits. Negative and positive interference effects were observed, as well as strong positional influences. Notably, blocking formation of a hypermodified pseudouridine in the P region delays the onset of the final cleavage event in 18S rRNA formation ( approximately 60% slower), suggesting that modification at this site could have an important role in modulating ribosome synthesis.

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.

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.

Computational prediction of novel non-coding RNAs in Arabidopsis thaliana

Background

Non-coding RNA (ncRNA) genes do not encode proteins but produce functional RNA molecules that play crucial roles in many key biological processes. Recent genome-wide transcriptional profiling studies using tiling arrays in organisms such as human and Arabidopsis have revealed a great number of transcripts, a large portion of which have little or no capability to encode proteins. This unexpected finding suggests that the currently known repertoire of ncRNAs may only represent a small fraction of ncRNAs of the organisms. Thus, efficient and effective prediction of ncRNAs has become an important task in bioinformatics in recent years. Among the available computational methods, the comparative genomic approach seems to be the most powerful to detect ncRNAs. The recent completion of the sequencing of several major plant genomes has made the approach possible for plants.

Results

We have developed a pipeline to predict novel ncRNAs in the Arabidopsis (Arabidopsis thaliana) genome. It starts by comparing the expressed intergenic regions of Arabidopsis as provided in two whole-genome high-density oligo-probe arrays from the literature with the intergenic nucleotide sequences of all completely sequenced plant genomes including rice (Oryza sativa), poplar (Populus trichocarpa), grape (Vitis vinifera), and papaya (Carica papaya). By using multiple sequence alignment, a popular ncRNA prediction program (RNAz), wet-bench experimental validation, protein-coding potential analysis, and stringent screening against various ncRNA databases, the pipeline resulted in 16 families of novel ncRNAs (with a total of 21 ncRNAs).

Conclusion

In this paper, we undertake a genome-wide search for novel ncRNAs in the genome of Arabidopsis by a comparative genomics approach. The identified novel ncRNAs are evolutionarily conserved between Arabidopsis and other recently sequenced plants, and may conduct interesting novel biological functions.

Rfam: updates to the RNA families database

Rfam is a collection of RNA sequence families, represented by multiple sequence alignments and covariance models (CMs). The primary aim of Rfam is to annotate new members of known RNA families on nucleotide sequences, particularly complete genomes, using sensitive BLAST filters in combination with CMs. A minority of families with a very broad taxonomic range (e.g. tRNA and rRNA) provide the majority of the sequence annotations, whilst the majority of Rfam families (e.g. snoRNAs and miRNAs) have a limited taxonomic range and provide a limited number of annotations. Recent improvements to the website, methodologies and data used by Rfam are discussed. Rfam is freely available on the Web at http://rfam.sanger.ac.uk/and http://rfam.janelia.org/.

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.

Diamonds in the rough: mRNA-like non-coding RNAs

Non-coding RNAs are increasingly being identified as crucial regulators of gene expression and other cellular functions in plants. Experimental and computational methods have revealed the existence of mRNA-like non-coding RNAs (mlncRNAs), a class of non-coding RNAs that, in plants, are associated with tissue-specific expression, development and the phosphate-starvation response. Although their mechanisms of action are largely unknown, one can speculate that mlncRNAs act through secondary structures or specific sequences that bind to proteins or metabolites, or that have catalytic activity. This review summarizes the computational methods developed to identify candidate mlncRNAs, and the current experimental evidence regarding the function of several known mlncRNAs.

Different mechanisms for pseudouridine formation in yeast 5S and 5.8S rRNAs

The selection of sites for pseudouridylation in eukaryotic cytoplasmic rRNA occurs by the base pairing of the rRNA with specific guide sequences within the RNA components of box H/ACA small nucleolar ribonucleoproteins (snoRNPs). Forty-four of the 46 pseudouridines (Psis) in the cytoplasmic rRNA of Saccharomyces cerevisiae have been assigned to guide snoRNAs. Here, we examine the mechanism of Psi formation in 5S and 5.8S rRNA in which the unassigned Psis occur. We show that while the formation of the Psi in 5.8S rRNA is associated with snoRNP activity, the pseudouridylation of 5S rRNA is not. The position of the Psi in 5.8S rRNA is guided by snoRNA snR43 by using conserved sequence elements that also function to guide pseudouridylation elsewhere in the large-subunit rRNA; an internal stem-loop that is not part of typical yeast snoRNAs also is conserved in snR43. The multisubstrate synthase Pus7 catalyzes the formation of the Psi in 5S rRNA at a site that conforms to the 7-nucleotide consensus sequence present in other substrates of Pus7. The different mechanisms involved in 5S and 5.8S rRNA pseudouridylation, as well as the multiple specificities of the individual trans factors concerned, suggest possible roles in linking ribosome production to other processes, such as splicing and tRNA synthesis.