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

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

Small nucleolar RNAs: versatile trans-acting molecules of ancient evolutionary origin

The small nucleolar RNAs (snoRNAs) are an abundant class of trans-acting RNAs that function in ribosome biogenesis in the eukaryotic nucleolus. Elegant work has revealed that most known snoRNAs guide modification of pre-ribosomal RNA (pre-rRNA) by base pairing near target sites. Other snoRNAs are involved in cleavage of pre-rRNA by mechanisms that have not yet been detailed. Moreover, our appreciation of the cellular roles of the snoRNAs is expanding with new evidence that snoRNAs also target modification of small nuclear RNAs and messenger RNAs. Many snoRNAs are produced by unorthodox modes of biogenesis including salvage from introns of pre-mRNAs. The recent discovery that homologs of snoRNAs as well as associated proteins exist in the domain Archaea indicates that the RNA-guided RNA modification system is of ancient evolutionary origin. In addition, it has become clear that the RNA component of vertebrate telomerase (an enzyme implicated in cancer and cellular senescence) is related to snoRNAs. During its evolution, vertebrate telomerase RNA appears to have co-opted a snoRNA domain that is essential for the function of telomerase RNA in vivo. The unique properties of snoRNAs are now being harnessed for basic research and therapeutic applications.

The discovery potential of RNA processing profiles

Small non-coding RNAs (sncRNAs) are highly abundant molecules that regulate essential cellular processes and are classified according to sequence and structure. Here we argue that read profiles from size-selected RNA sequencing capture the post-transcriptional processing specific to each RNA family, thereby providing functional information independently of sequence and structure. We developed SeRPeNT, a new computational method that exploits reproducibility across replicates and uses dynamic time-warping and density-based clustering algorithms to identify, characterize and compare sncRNAs by harnessing the power of read profiles. We applied SeRPeNT to: (i) generate an extended human annotation with 671 new sncRNAs from known classes and 131 from new potential classes, (ii) show pervasive differential processing of sncRNAs between cell compartments and (iii) predict new molecules with miRNA-like behaviour from snoRNA, tRNA and long non-coding RNA precursors, potentially dependent on the miRNA biogenesis pathway. Furthermore, we validated experimentally four predicted novel non-coding RNAs: a miRNA, a snoRNA-derived miRNA, a processed tRNA and a new uncharacterized sncRNA. SeRPeNT facilitates fast and accurate discovery and characterization of sncRNAs at an unprecedented scale. SeRPeNT code is available under the MIT license at https://github.com/comprna/SeRPeNT.

Non-coding RNAs: Functions and applications in endocrine-related cancer

A significant fraction of the human genome is transcribed as non-coding RNAs (ncRNAs). This non-coding transcriptome has challenged the notion of the central dogma and its involvement in transcriptional and post-transcriptional regulation of gene expression is well established. Interestingly, several ncRNAs are dysregulated in cancer and current non-coding transcriptome research aims to use our increasing knowledge of these ncRNAs for the development of cancer biomarkers and anti-cancer drugs. In endocrine-related cancers, for which survival rates can be relatively low, there is a need for such advancements. In this review, we aimed to summarize the roles and clinical implications of recently discovered ncRNAs, including long ncRNAs, PIWI-interacting RNAs, tRNA- and Y RNA-derived ncRNAs, and small nucleolar RNAs, in endocrine-related cancers affecting both sexes. We focus on recent studies highlighting discoveries in ncRNA biology and expression in cancer, and conclude with a discussion on the challenges and future directions, including clinical application. ncRNAs show great promise as diagnostic tools and therapeutic targets, but further work is necessary to realize the potential of these unconventional transcripts.

Identification of discrete classes of small nucleolar RNA featuring different ends and RNA binding protein dependency

Small nucleolar RNAs (snoRNAs) are among the first discovered and most extensively studied group of small non-coding RNA. However, most studies focused on a small subset of snoRNAs that guide the modification of ribosomal RNA. In this study, we annotated the expression pattern of all box C/D snoRNAs in normal and cancer cell lines independent of their functions. The results indicate that C/D snoRNAs are expressed as two distinct forms differing in their ends with respect to boxes C and D and in their terminal stem length. Both forms are overexpressed in cancer cell lines but display a conserved end distribution. Surprisingly, the long forms are more dependent than the short forms on the expression of the core snoRNP protein NOP58, thought to be essential for C/D snoRNA production. In contrast, a subset of short forms are dependent on the splicing factor RBFOX2. Analysis of the potential secondary structure of both forms indicates that the k-turn motif required for binding of NOP58 is less stable in short forms which are thus less likely to mature into a canonical snoRNP. Taken together the data suggest that C/D snoRNAs are divided into at least two groups with distinct maturation and functional preferences.

Generation of a mouse model with down-regulated U50 snoRNA (SNORD50) expression and its organ-specific phenotypic modulation

Box C/D-type small nucleolar RNAs (snoRNAs) are functional RNAs responsible for mediating 2'-O-ribose methylation of ribosomal RNAs (rRNAs) within the nucleolus. In the past years, evidence for the involvement of human U50 snoRNA in tumorigenesis has been accumulating. We previously identified U50HG, a non-protein-coding gene that hosted a box C/D-type U50 snoRNA, in a chromosomal breakpoint in a human B-cell lymphoma. Mouse genome analysis revealed four mouse U50 (mU50) host-genes: three mU50HG-a gene variants that were clustered in the genome and an mU50HG-b gene that we supposed to be the U50HG ortholog. In this study, to investigate the physiological importance of mU50 snoRNA and its involvement in tumorigenesis, we eliminated mU50 snoRNA sequences from the mU50HG-b gene. The established mouse line (ΔmU50(HG-b)) showed a significant reduction of mU50 snoRNA expression without alteration of the host-gene length and exon-intron structure, and the corresponding target rRNA methylation in various organs was reduced. Lifelong phenotypic monitoring showed that the ΔmU50(HG-b) mice looked almost normal without accelerated tumorigenicity; however, a notable difference was the propensity for anomalies in the lymphoid organs. Transcriptome analysis showed that dozens of genes, including heat shock proteins, were differentially expressed in ΔmU50(HG-b) mouse lymphocytes. This unique model of a single snoRNA knockdown with intact host-gene expression revealed further new insights into the discrete transcriptional regulation of multiple mU50 host-genes and the complicated dynamics involved in organ-specific processing and maintenance of snoRNAs.

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.

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.

Identification and classification of small RNAs in transcriptome sequence data

Current methods for high throughput sequencing (HTS) for the first time offer the opportunity to investigate the entire transcriptome in an essentially unbiased way. In many species, small non-coding RNAs with specific secondary structures constitute a significant part of the transcriptome. Some of these RNA classes, in particular microRNAs and snoRNAs, undergo maturation processes that lead to the production of shorter RNAs. After mapping the sequences to the reference genome specific patterns of short reads can be observed. These read patterns seem to reflect the processing and thus are specific for the RNA transcripts of which they are derived from. We explore here the potential of short read sequence data in the classification and identification of non-coding RNAs.

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.

[Small nucleolar RNAs]

Small nucleolar RNAs (snoRNAs) are an abundant class of non-protein-coding RNAs. In association with proteins they perform two most frequent nucleotide modifications in rRNAs and some other cellular RNAs: 2'-O-ribose methylation and pseudouridylation. SnoRNAs also participate in pre-rRNA cleavage and telomerase functions. Most snoRNAs fall into two families, box C/D and H/ACA, distinguished by the presence of conserved sequence boxes. Although C/D and H/ACA snoRNP proteins contain homologous regions, the assembly of these RNPs significantly differ. In addition, snoRNAs include the RNA component of RNAses P and MRP. The structure and function of small RNPs from Cajal bodies (small organelles associated with nucleoli) similar to snoRNP are also discussed.

Evolution of small nucleolar RNAs in nematodes

In contrast to mRNAs, which are templates for translating proteins, non-protein coding (npc) RNAs (also known as ‘non-coding’ RNA, ncRNA), exhibit various functions in different compartments and developmental stages of the cell. Small nucleolar RNAs (snoRNAs), one of the largest classes of npcRNAs, guide post-transcriptional modifications of other RNAs that are crucial for appropriate RNA folding as well as for RNA–RNA and RNA–protein interactions. Although snoRNA genes comprise a significant fraction of the eutherian genome, identifying and characterizing large numbers of them is not sufficiently accessible by classical computer searches alone. Furthermore, most previous investigations of snoRNAs yielded only limited indications of their evolution. Using data obtained by a combination of high-throughput cDNA library screening and computational search strategies based on a modified DNAMAN program, we characterized 151 npcRNAs, and in particular 121 snoRNAs, from Caenorhabditis elegans and extensively compared them with those in the related, Caenorhabditis briggsae. Detailed comparisons of paralog snoRNAs in the two nematodes revealed, in addition to trans-duplication, a novel, cis-duplication distribution strategy with insertions near to the original loci. Some snoRNAs coevolved with their modification target sites, demonstrating the close interaction of complementary regions. Some target sites modified by snoRNAs were changed, added or lost, documenting a high degree of evolutionary plasticity of npcRNAs.

Expression and Function of Brain Specific Small RNAs

Small non-messenger RNAs (snmRNAs) are a heterogeneous group of non-coding RNAs with a variety of regulatory functions including regulation of protein expression and guidance in RNA modifications. They are actively being investigated in Archaebacteria, yeast, invertebrates and mammals. Brain-specific snmRNAs have been identified in mammals and they seem to contribute to neuronal differentiation during development and to brain functions subserving learning and memory. Here we review the current knowledge of the properties, expression and functions of three groups of brain-specific snmRNAs: small nucleolar RNAs, BC1/BC200 RNAs and microRNAs.

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

MOTIVATION

Noncoding RNA genes produce functional RNA molecules rather than coding for proteins. One such family is the H/ACA snoRNAs. Unlike the related C/D snoRNAs these have resisted automated detection to date.

RESULTS

We develop an algorithm to screen the yeast genome for novel H/ACA snoRNAs. To achieve this, we introduce some new methods for facilitating the search for noncoding RNAs in genomic sequences which are based on properties of predicted minimum free-energy (MFE) secondary structures. The algorithm has been implemented and can be generalized to enable screening of other eukaryote genomes. We find that use of primary sequence alone is insufficient for identifying novel H/ACA snoRNAs. Only the use of secondary structure filters reduces the number of candidates to a manageable size. From genomic context, we identify three strong H/ACA snoRNA candidates. These together with a further 47 candidates obtained by our analysis are being experimentally screened.

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

Dozens of box C/D small nucleolar RNAs (snoRNAs) have recently been found in eukaryotes (vertebrates, yeast), ancient eukaryotes (trypanosomes) and archae, that specifically target ribosomal RNA sites for 2'-O-ribose methylation. Although early biochemical data revealed that plant rRNAs are among the most highly ribomethylated in eukaryotes, only a handful of methylation guide snoRNAs have been characterized in this kingdom. We report 66 novel box C/D snoRNAs identified by computational screening of Arabidopsis genomic sequences that are expressed in vivo from either single genes, 17 different clusters or three introns. At the structural level, many box C/D snoRNAs have dual antisense elements often matching rRNA regions close to each other on the rRNA secondary structure, which is reminiscent of their archaeal counterparts. Remarkable specimens are found that display two antisense elements having the potential to form an extended snoRNA-rRNA duplex of 23 to 30 nt, in line with the hypothetical function of box C/D snoRNAs in pre-rRNA folding or chaperoning. In contrast to other species, many Arabidopsis snoRNAs are found in multiple isoforms mainly resulting from two different mechanisms: large chromosomal duplications and small tandem duplications producing polycistronic genes. The discovery of numerous different snoRNAs, some of them arising from common ancestors, provide new insights to understand snoRNAs evolution and the birth of new rRNA methylation sites in plants and other organisms.