Small nucleolar RNAs that guide modification in trypanosomatids: repertoire, targets, genome organisation, and unique functions

The fungal snoRNAome

Small nucleolar RNAs (snoRNAs) are essential players in the rRNA biogenesis due to their involvement in the nucleolytic processing of the precursor and the subsequent guidance of nucleoside modifications. Within the kingdom Fungi, merely a few species-specific surveys have explored their snoRNA repertoire. However, the wide range of the snoRNA landscape spanning all major fungal lineages has not been mapped so far, mainly because of missing tools for automatized snoRNA detection and functional analysis. For the first time, we report here a comprehensive inventory of fungal snoRNAs together with a functional analysis and an in-depth investigation of their evolutionary history including innovations, deletions, and target switches. This large-scale analysis, incorporating more than 120 snoRNA families with more than 7700 individual snoRNA sequences, catalogs and clarifies the landscape of fungal snoRNA families, assigns functions to previously orphan snoRNAs, and increases the number of sequences by 450%. We also show that the snoRNAome is subject to ongoing rearrangements and adaptations, e.g., through lineage-specific targets and redundant guiding functions.

The Canonical Poly (A) Polymerase PAP1 Polyadenylates Non-Coding RNAs and Is Essential for snoRNA Biogenesis in Trypanosoma brucei

The parasite Trypanosoma brucei is the causative agent of African sleeping sickness and is known for its unique RNA processing mechanisms that are common to all the kinetoplastidea including Leishmania and Trypanosoma cruzi. Trypanosomes possess two canonical RNA poly (A) polymerases (PAPs) termed PAP1 and PAP2. PAP1 is encoded by one of the only two genes harboring cis-spliced introns in this organism, and its function is currently unknown. In trypanosomes, all mRNAs, and non-coding RNAs such as small nucleolar RNAs (snoRNAs) and long non-coding RNAs (lncRNAs), undergo trans-splicing and polyadenylation. Here, we show that the function of PAP1, which is located in the nucleus, is to polyadenylate non-coding RNAs, which undergo trans-splicing and polyadenylation. Major substrates of PAP1 are the snoRNAs and lncRNAs. Under the silencing of either PAP1 or PAP2, the level of snoRNAs is reduced. The dual polyadenylation of snoRNA intermediates is carried out by both PAP2 and PAP1 and requires the factors essential for the polyadenylation of mRNAs. The dual polyadenylation of the precursor snoRNAs by PAPs may function to recruit the machinery essential for snoRNA processing.

Synthesis, Function, and Heterogeneity of snoRNA-Guided Posttranscriptional Nucleoside Modifications in Eukaryotic Ribosomal RNAs

Ribosomal RNAs contain numerous 2'-O-methylated nucleosides and pseudouridines. Methylation of the 2' oxygen of ribose moieties and isomerization of uridines into pseudouridines are catalyzed by C/D and H/ACA small nucleolar ribonucleoprotein particles, respectively. We review the composition, structure, and mode of action of archaeal and eukaryotic C/D and H/ACA particles. Most rRNA modifications cluster in functionally crucial regions of the rRNAs, suggesting they play important roles in translation. Some of these modifications promote global translation efficiency or modulate translation fidelity. Strikingly, recent quantitative nucleoside modification profiling methods have revealed that a subset of modification sites is not always fully modified. The finding of such ribosome heterogeneity is in line with the concept of specialized ribosomes that could preferentially translate specific mRNAs. This emerging concept is supported by findings that some human diseases are caused by defects in the rRNA modification machinery correlated with a significant alteration of IRES-dependent translation.

Transcriptional regulation of snRNAs and its significance for plant development

Small nuclear RNA (snRNA) represents a distinct class of non-coding RNA molecules. As these molecules have fundamental roles in RNA metabolism, including pre-mRNA splicing and ribosomal RNA processing, it is essential that their transcription be tightly regulated in eukaryotic cells. The genome of each organism contains hundreds of snRNA genes. Although the structures of these genes are highly diverse among organisms, the trans-acting factors that regulate snRNA transcription are evolutionarily conserved. Recent studies of the Arabidopsis thaliana srd2-1 mutant, which is defective in the snRNA transcription factor, provide insight into the physiological significance of snRNA regulation in plant development. Here, I review the current understanding of the molecular mechanisms underlying snRNA transcription.

A mass spectrometry-based method for direct determination of pseudouridine in RNA

Pseudouridine (5-ribosyluracil, Ψ) is the only ‘mass-silent’ nucleoside produced by post-transcriptional RNA modification. We describe here a novel mass spectrometry (MS)-based method for direct determination of Ψ in RNA. The method assigns a Ψ-containing nucleolytic RNA fragment by an accurate measurement of a signature doubly dehydrated nucleoside anion ([C₉H₇N₂O₄]¹⁻, m/z 207.04) produced by collision-induced dissociation MS, and it determines the Ψ-containing nucleotide sequence by pseudo-MS³, i.e. in-source fragmentation followed by MS². By applying this method, we identified all of the known Ψs in the canonical human spliceosomal snRNAs and, unexpectedly, found two previously unknown Ψs in the U5 and U6 snRNAs. Because the method allows direct determination of Ψ in a subpicomole quantity of RNA, it will serve as a useful tool for the structure/function studies of a wide variety of non-coding RNAs.

Genome-wide analysis of small nucleolar RNAs of Leishmania major reveals a rich repertoire of RNAs involved in modification and processing of rRNA

Trypanosomatids are protozoan parasites and the causative agent of infamous infectious diseases. These organisms regulate their gene expression mainly at the post-transcriptional level and possess characteristic RNA processing mechanisms. In this study, we analyzed the complete repertoire of Leishmania major small nucleolar (snoRNA) RNAs by performing RNA-seq analysis on RNAs that were affinity-purified using the C/D snoRNA core protein, SNU13, and the H/ACA core protein, NHP2. This study revealed a large collection of C/D and H/ACA snoRNAs, organized in gene clusters generally containing both snoRNA types. Abundant snoRNAs were identified and predicted to guide trypanosome-specific rRNA cleavages. The repertoire of snoRNAs was compared to that of the closely related Trypanosoma brucei, and 80% of both C/D and H/ACA molecules were found to have functional homologues. The comparative analyses elucidated how snoRNAs evolved to generate molecules with analogous functions in both species. Interestingly, H/ACA RNAs have great flexibility in their ability to guide modifications, and several of the RNA species can guide more than one modification, compensating for the presence of single hairpin H/ACA snoRNA in these organisms. Placing the predicted modifications on the rRNA secondary structure revealed hypermodification regions mostly in domains which are modified in other eukaryotes, in addition to trypanosome-specific modifications.

Small nucleolar RNAs and RNA-guided post-transcriptional modification

snoRNAs (small nucleolar RNAs) constitute one of the largest and best-studied classes of non-coding RNAs that confer enzymatic specificity. With associated proteins, these snoRNAs form ribonucleoprotein complexes that can direct 2'-O-methylation or pseudouridylation of target non-coding RNAs. Aided by computational methods and high-throughput sequencing, new studies have expanded the diversity of known snoRNA functions. Complexes incorporating snoRNAs have dynamic specificity, and include diverse roles in RNA silencing, telomerase maintenance and regulation of alternative splicing. Evidence that dysregulation of snoRNAs can cause human disease, including cancer, indicates that the full scope of snoRNA roles remains an unfinished story. The diversity in structure, genomic origin and function between snoRNAs found in different complexes and among different phyla illustrates the surprising plasticity of snoRNAs in evolution. The ability of snoRNAs to direct highly specific interactions with other RNAs is a consistent thread in their newly discovered functions. Because they are ubiquitous throughout Eukarya and Archaea, it is likely they were a feature of the last common ancestor of these two domains, placing their origin over two billion years ago. In the present chapter, we focus on recent advances in our understanding of these ancient, but functionally dynamic RNA-processing machines.

microRNAs in parasites and parasite infection

miRNAs, a subclass of small regulatory RNAs, are present from ancient unicellular protozoans to parasitic helminths and parasitic arthropods. The miRNA-silencing mechanism appears, however, to be absent in a number of protozoan parasites. Protozoan miRNAs and components of their silencing machinery possess features different from other eukaryotes, providing some clues on the evolution of the RNA-induced silencing machinery. miRNA functions possibly associate with neoblast biology, development, physiology, infection and immunity of parasites. Parasite infection can alter host miRNA expression that can favor both parasite clearance and infection. miRNA pathways are, thus, a potential target for the therapeutic control of parasitic diseases.

Pseudouridine monophosphate glycosidase: a new glycosidase mechanism

Pseudouridine (Ψ), the most abundant modification in RNA, is synthesized in situ using Ψ synthase. Recently, a pathway for the degradation of Ψ was described [Preumont, A., Snoussi, K., Stroobant, V., Collet, J. F., and Van Schaftingen, E. (2008) J. Biol. Chem. 283, 25238-25246]. In this pathway, Ψ is first converted to Ψ 5'-monophosphate (ΨMP) by Ψ kinase and then ΨMP is degraded by ΨMP glycosidase to uracil and ribose 5-phosphate. ΨMP glycosidase is the first example of a mechanistically characterized enzyme that cleaves a C-C glycosidic bond. Here we report X-ray crystal structures of Escherichia coli ΨMP glycosidase and a complex of the K166A mutant with ΨMP. We also report the structures of a ring-opened ribose 5-phosphate adduct and a ring-opened ribose ΨMP adduct. These structures provide four snapshots along the reaction coordinate. The structural studies suggested that the reaction utilizes a Lys166 adduct during catalysis. Biochemical and mass spectrometry data further confirmed the existence of a lysine adduct. We used site-directed mutagenesis combined with kinetic analysis to identify roles for specific active site residues. Together, these data suggest that ΨMP glycosidase catalyzes the cleavage of the C-C glycosidic bond through a novel ribose ring-opening mechanism.

Evolutionary patterns of non-coding RNAs

A plethora of new functions of non-coding RNAs (ncRNAs) have been discovered in past few years. In fact, RNA is emerging as the central player in cellular regulation, taking on active roles in multiple regulatory layers from transcription, RNA maturation, and RNA modification to translational regulation. Nevertheless, very little is known about the evolution of this "Modern RNA World" and its components. In this contribution, we attempt to provide at least a cursory overview of the diversity of ncRNAs and functional RNA motifs in non-translated regions of regular messenger RNAs (mRNAs) with an emphasis on evolutionary questions. This survey is complemented by an in-depth analysis of examples from different classes of RNAs focusing mostly on their evolution in the vertebrate lineage. We present a survey of Y RNA genes in vertebrates and study the molecular evolution of the U7 snRNA, the snoRNAs E1/U17, E2, and E3, the Y RNA family, the let-7 microRNA (miRNA) family, and the mRNA-like evf-1 gene. We furthermore discuss the statistical distribution of miRNAs in metazoans, which suggests an explosive increase in the miRNA repertoire in vertebrates. The analysis of the transcription of ncRNAs suggests that small RNAs in general are genetically mobile in the sense that their association with a hostgene (e.g. when transcribed from introns of a mRNA) can change on evolutionary time scales. The let-7 family demonstrates, that even the mode of transcription (as intron or as exon) can change among paralogous ncRNA.

Characterizing ncRNAs in Human Pathogenic Protists Using High-Throughput Sequencing Technology

ncRNAs are key genes in many human diseases including cancer and viral infection, as well as providing critical functions in pathogenic organisms such as fungi, bacteria, viruses, and protists. Until now the identification and characterization of ncRNAs associated with disease has been slow or inaccurate requiring many years of testing to understand complicated RNA and protein gene relationships. High-throughput sequencing now offers the opportunity to characterize miRNAs, siRNAs, small nucleolar RNAs (snoRNAs), and long ncRNAs on a genomic scale, making it faster and easier to clarify how these ncRNAs contribute to the disease state. However, this technology is still relatively new, and ncRNA discovery is not an application of high priority for streamlined bioinformatics. Here we summarize background concepts and practical approaches for ncRNA analysis using high-throughput sequencing, and how it relates to understanding human disease. As a case study, we focus on the parasitic protists Giardia lamblia and Trichomonas vaginalis, where large evolutionary distance has meant difficulties in comparing ncRNAs with those from model eukaryotes. A combination of biological, computational, and sequencing approaches has enabled easier classification of ncRNA classes such as snoRNAs, but has also aided the identification of novel classes. It is hoped that a higher level of understanding of ncRNA expression and interaction may aid in the development of less harsh treatment for protist-based diseases.

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

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

Efficient and specific knockdown of small non-coding RNAs in mammalian cells and in mice

Hundreds of small nuclear non-coding RNAs, including small nucleolar RNAs (snoRNAs), have been identified in different organisms, with important implications in regulating gene expression and in human diseases. However, functionalizing these nuclear RNAs in mammalian cells remains challenging, due to methodological difficulties in depleting these RNAs, especially snoRNAs. Here we report a convenient and efficient approach to deplete snoRNA, small Cajal body RNA (scaRNA) and small nuclear RNA in human and mouse cells by conventional transfection of chemically modified antisense oligonucleotides (ASOs) that promote RNaseH-mediated cleavage of target RNAs. The levels of all seven tested snoRNA/scaRNAs and four snRNAs were reduced by 80-95%, accompanied by impaired endogenous functions of the target RNAs. ASO-targeting is highly specific, without affecting expression of the host genes where snoRNAs are embedded in the introns, nor affecting the levels of snoRNA isoforms with high sequence similarities. At least five snoRNAs could be depleted simultaneously. Importantly, snoRNAs could be dramatically depleted in mice by systematic administration of the ASOs. Together, our findings provide a convenient and efficient approach to characterize nuclear non-coding RNAs in mammalian cells, and to develop antisense drugs against disease-causing non-coding RNAs.

A comparative genome-wide study of ncRNAs in trypanosomatids

Background

Recent studies have provided extensive evidence for multitudes of non-coding RNA (ncRNA) transcripts in a wide range of eukaryotic genomes. ncRNAs are emerging as key players in multiple layers of cellular regulation. With the availability of many whole genome sequences, comparative analysis has become a powerful tool to identify ncRNA molecules. In this study, we performed a systematic genome-wide in silico screen to search for novel small ncRNAs in the genome of Trypanosoma brucei using techniques of comparative genomics.

Results

In this study, we identified by comparative genomics, and validated by experimental analysis several novel ncRNAs that are conserved across multiple trypanosomatid genomes. When tested on known ncRNAs, our procedure was capable of finding almost half of the known repertoire through homology over six genomes, and about two-thirds of the known sequences were found in at least four genomes. After filtering, 72 conserved unannotated sequences in at least four genomes were found, 29 of which, ranging in size from 30 to 392 nts, were conserved in all six genomes. Fifty of the 72 candidates in the final set were chosen for experimental validation. Eighteen of the 50 (36%) were shown to be expressed, and for 11 of them a distinct expression product was detected, suggesting that they are short ncRNAs. Using functional experimental assays, five of the candidates were shown to be novel H/ACA and C/D snoRNAs; these included three sequences that appear as singletons in the genome, unlike previously identified snoRNA molecules that are found in clusters. The other candidates appear to be novel ncRNA molecules, and their function is, as yet, unknown.

Conclusions

Using comparative genomic techniques, we predicted 72 sequences as ncRNA candidates in T. brucei. The expression of 50 candidates was tested in laboratory experiments. This resulted in the discovery of 11 novel short ncRNAs in procyclic stage T. brucei, which have homologues in the other trypansomatids. A few of these molecules are snoRNAs, but most of them are novel ncRNA molecules. Based on this study, our analysis suggests that the total number of ncRNAs in trypanosomatids is in the range of several hundred.

Small nucleolar RNA interference in Trypanosoma brucei: mechanism and utilization for elucidating the function of snoRNAs

Expression of dsRNA complementary to small nucleolar RNAs (snoRNAs) in Trypanosoma brucei results in snoRNA silencing, termed snoRNAi. Here, we demonstrate that snoRNAi requires the nuclear TbDCL2 protein, but not TbDCL1, which is involved in RNA interference (RNAi) in the cytoplasm. snoRNAi depends on Argonaute1 (Slicer), and on TbDCL2, suggesting that snoRNA dicing and slicing takes place in the nucleus, and further suggesting that AGO1 is active in nuclear silencing. snoRNAi was next utilized to elucidate the function of an abundant snoRNA, TB11Cs2C2 (92 nt), present in a cluster together with the spliced leader associated RNA (SLA1) and snR30, which are both H/ACA RNAs with special nuclear functions. Using AMT-UV cross-linking and RNaseH cleavage, we provide evidence for the interaction of TB11Cs2C2 with the small rRNAs, srRNA-2 and srRNA-6, which are part of the large subunit (LSU) rRNA. snoRNAi of TB11Cs2C2 resulted in defects in generating srRNA-2 and LSUβ rRNA. This is the first snoRNA described so far to engage in trypanosome-specific processing events.

Specificity and stoichiometry of subunit interactions in the human telomerase holoenzyme assembled in vivo

The H/ACA motif of human telomerase RNA (hTR) directs specific pathways of endogenous telomerase holoenzyme assembly, function, and regulation. Similarities between hTR and other H/ACA RNAs have been established, but differences have not been explored even though unique features of hTR H/ACA RNP assembly give rise to telomerase deficiency in human disease. Here, we define hTR H/ACA RNA and RNP architecture using RNA accumulation, RNP affinity purification, and primer extension activity assays. First, we evaluate alternative folding models for the hTR H/ACA motif 5' hairpin. Second, we demonstrate an unanticipated and surprisingly general asymmetry of 5' and 3' hairpin requirements for H/ACA RNA accumulation. Third, we establish that hTR assembles not one but two sets of all four of the H/ACA RNP core proteins, dyskerin, NOP10, NHP2, and GAR1. Fourth, we address a difference in predicted specificities of hTR association with the holoenzyme subunit WDR79/TCAB1. Together, these results complete the analysis of hTR elements required for active RNP biogenesis and define the interaction specificities and stoichiometries of all functionally essential human telomerase holoenzyme subunits. This study uncovers unexpected similarities but also differences between telomerase and other H/ACA RNPs that allow a unique specificity of telomerase biogenesis and regulation.

RNAsnoop: efficient target prediction for H/ACA snoRNAs

MOTIVATION

Small nucleolar RNAs are an abundant class of non-coding RNAs that guide chemical modifications of rRNAs, snRNAs and some mRNAs. In the case of many 'orphan' snoRNAs, the targeted nucleotides remain unknown, however. The box H/ACA subclass determines uridine residues that are to be converted into pseudouridines via specific complementary binding in a well-defined secondary structure configuration that is outside the scope of common RNA (co-)folding algorithms.

RESULTS

RNAsnoop implements a dynamic programming algorithm that computes thermodynamically optimal H/ACA-RNA interactions in an efficient scanning variant. Complemented by an support vector machine (SVM)-based machine learning approach to distinguish true binding sites from spurious solutions and a system to evaluate comparative information, it presents an efficient and reliable tool for the prediction of H/ACA snoRNA target sites. We apply RNAsnoop to identify the snoRNAs that are responsible for several of the remaining 'orphan' pseudouridine modifications in human rRNAs, and we assign a target to one of the five orphan H/ACA snoRNAs in Drosophila.

AVAILABILITY

The C source code of RNAsnoop is freely available at http://www.tbi.univie.ac.at/ -htafer/RNAsnoop

U3 snoRNA genes are multi-copy and frequently linked to U5 snRNA genes in Euglena gracilis

Background

U3 snoRNA is a box C/D small nucleolar RNA (snoRNA) involved in the processing events that liberate 18S rRNA from the ribosomal RNA precursor (pre-rRNA). Although U3 snoRNA is present in all eukaryotic organisms, most investigations of it have focused on fungi (particularly yeasts), animals and plants. Relatively little is known about U3 snoRNA and its gene(s) in the phylogenetically broad assemblage of protists (mostly unicellular eukaryotes). In the euglenozoon Euglena gracilis, a distant relative of the kinetoplastid protozoa, Southern analysis had previously revealed at least 13 bands hybridizing with U3 snoRNA, suggesting the existence of multiple copies of U3 snoRNA genes.

Results

Through screening of a λ genomic library and PCR amplification, we recovered 14 U3 snoRNA gene variants, defined by sequence heterogeneities that are mostly located in the U3 3'-stem-loop domain. We identified three different genomic arrangements of Euglena U3 snoRNA genes: i) stand-alone, ii) linked to tRNA^Arg genes, and iii) linked to a U5 snRNA gene. In arrangement ii), the U3 snoRNA gene is positioned upstream of two identical tRNA^Arg genes that are convergently transcribed relative to the U3 gene. This scenario is reminiscent of a U3 snoRNA-tRNA gene linkage previously described in trypanosomatids. We document here twelve different U3 snoRNA-U5 snRNA gene arrangements in Euglena; in each case, the U3 gene is linked to a downstream and convergently oriented U5 gene, with the intergenic region differing in length and sequence among the variants.

Conclusion

The multiple U3 snoRNA-U5 snRNA gene linkages, which cluster into distinct families based on sequence similarities within the intergenic spacer, presumably arose by genome, chromosome, and/or locus duplications. We discuss possible reasons for the existence of the unusually large number of U3 snoRNA genes in the Euglena genome. Variability in the signal intensities of the multiple Southern hybridization bands raises the possibility that Euglena contains a naturally aneuploid chromosome complement.

Cl gene cluster encoding several small nucleolar RNAs: a comparison amongst trypanosomatids

Small nucleolar RNAs (snoRNAs) are small non-coding RNAs that modify RNA molecules such as rRNA and snRNA by guiding 2'-O-ribose methylation (C/D box snoRNA family) and pseudouridylation reactions (H/ACA snoRNA family). H/ACA snoRNAs are also involved in trans-splicing in trypanosomatids. The aims of this work were to characterise the Cl gene cluster that encodes several snoRNAs in Trypanosoma rangeli and compare it with clusters from Trypanosoma cruzi, Trypanosoma brucei, Leishmania major, Leishmania infantum, Leishmania braziliensis and Leptomonas collosoma. The T. rangeli Cl gene cluster is an 801 base pair (bp) repeat sequence that encodes three C/D (Cl1, Cl2 and Cl4) and three H/ACA (Cl3, Cl5 and Cl6) snoRNAs. In contrast to T. brucei, the Cl3 and Cl5 homologues have not been annotated in the Leishmania or T. cruzi genome projects (http//:www.genedb.org). Of note, snoRNA transcribed regions have a high degree of sequence identity among all species and share gene synteny. Collectively, these findings suggest that the Cl cluster could constitute an interesting target for therapeutic (gene silencing) or diagnostic intervention strategies (PCR-derived tools).

snoRNAs in Giardia lamblia: a novel role in RNA silencing?

In the expanding world of small regulatory RNAs, a recent paper by Saraiya and Wang has reported the identification in the protozoan parasite Giardia lamblia of a novel class of small RNAs, which are derived by Dicer processing of small nucleolar RNAs and have the potential to function as micro RNAs. Interestingly, these RNAs occur not only in this parasite but also in humans.

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.

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.

Introduction to special issue on RNA

In this introduction to the special issue on RNA, we provide a brief overview of some of the novel and exciting biological discoveries concerning diverse roles played by RNA, and subsequently we give a rapid summary of some algorithmic aspects of RNA structure and alignment. Each of the contributions to this special issue is briefly described.

Identification of the heptameric Lsm complex that binds U6 snRNA in Trypanosoma brucei

Lsm proteins are ubiquitous, multifunctional proteins that are involved in nuclear processing and turnover of many RNAs in eukaryotes. Lsm proteins form two distinct complexes, the Lsm2-8 complex, which binds U6 snRNA, and the Lsm1-7 complex, which governs mRNA degradation. Previously, seven Lsm proteins were identified in Trypanosoma brucei. Two of these proteins were later identified as SSm proteins (specific spliceosomal Sm proteins). In this study, the Lsm proteins (Lsm2 and Lsm5) that bind to U6 snRNA were identified. RNAi silencing and protein purification of TAP-tagged Lsm proteins were used to identify all the components of the trypanosome heptameric Lsm2-8 complex. Localization studies demonstrated that these proteins are found in the nucleus, near the nucleolus. Lsm proteins were not detected in cytoplasmic bodies that were tagged with YFP-Dhh1, which may suggest that in trypanosomes, Lsm-mediated degradation is not confined to such bodies.

Elucidating the role of C/D snoRNA in rRNA processing and modification in Trypanosoma brucei

Most eukaryotic C/D small nucleolar RNAs (snoRNAs) guide 2'-O methylation (Nm) on rRNA and are also involved in rRNA processing. The four core proteins that bind C/D snoRNA in Trypanosoma brucei are fibrillarin (NOP1), NOP56, NOP58, and SNU13. Silencing of NOP1 by RNA interference identified rRNA-processing and modification defects that caused lethality. Systematic mapping of 2'-O-methyls on rRNA revealed the existence of hypermethylation at certain positions of the rRNA in the bloodstream form of the parasites, suggesting that this modification may assist the parasites in coping with the major temperature changes during cycling between their insect and mammalian hosts. The rRNA-processing defects of NOP1-depleted cells suggest the involvement of C/D snoRNA in trypanosome-specific rRNA-processing events to generate the small rRNA fragments. MRP RNA, which is involved in rRNA processing, was identified in this study in one of the snoRNA gene clusters, suggesting that trypanosomes utilize a combination of unique C/D snoRNAs and conserved snoRNAs for rRNA processing.

Fine-tuning of translation termination efficiency in Saccharomyces cerevisiae involves two factors in close proximity to the exit tunnel of the ribosome

In eukaryotes, release factors 1 and 3 (eRF1 and eRF3) are recruited to promote translation termination when a stop codon on the mRNA enters at the ribosomal A-site. However, their overexpression increases termination efficiency only moderately, suggesting that other factors might be involved in the termination process. To determine such unknown components, we performed a genetic screen in Saccharomyces cerevisiae that identified genes increasing termination efficiency when overexpressed. For this purpose, we constructed a dedicated reporter strain in which a leaky stop codon is inserted into the chromosomal copy of the ade2 gene. Twenty-five antisuppressor candidates were identified and characterized for their impact on readthrough. Among them, SSB1 and snR18, two factors close to the exit tunnel of the ribosome, directed the strongest antisuppression effects when overexpressed, showing that they may be involved in fine-tuning of the translation termination level.

Treasure hunt in an amoeba: non-coding RNAs in Dictyostelium discoideum

The traditional view of RNA being merely an intermediate in the transfer of genetic information, as mRNA, spliceosomal RNA, tRNA, and rRNA, has become outdated. The recent discovery of numerous regulatory RNAs with a plethora of functions in biological processes has truly revolutionized our understanding of gene regulation. Tiny RNAs such as microRNAs and small interfering RNAs play vital roles at different levels of gene control. Small nucleolar RNAs are much more abundant than previously recognized, and new functions beyond processing and modification of rRNA have recently emerged. Longer non-coding RNAs (ncRNAs) can also have important regulatory roles in the cell, e.g., antisense RNAs that control their target mRNAs. The majority of these important findings arose from analyses in various model organisms. In this review, we focus on ncRNAs in the social amoeba Dictyostelium discoideum. This important genetically tractable model organism has recently received renewed attention in terms of discovery, regulation and functional studies of ncRNAs. Old and recent findings are discussed and put in context of what we today know about ncRNAs in other organisms.

Combined experimental and computational approach to identify non-protein-coding RNAs in the deep-branching eukaryote Giardia intestinalis

Non-protein-coding RNAs represent a large proportion of transcribed sequences in eukaryotes. These RNAs often function in large RNA-protein complexes, which are catalysts in various RNA-processing pathways. As RNA processing has become an increasingly important area of research, numerous non-messenger RNAs have been uncovered in all the model eukaryotic organisms. However, knowledge on RNA processing in deep-branching eukaryotes is still limited. This study focuses on the identification of non-protein-coding RNAs from the diplomonad parasite Giardia intestinalis, showing that a combined experimental and computational search strategy is a fast method of screening reduced or compact genomes. The analysis of our Giardia cDNA library has uncovered 31 novel candidates, including C/D-box and H/ACA box snoRNAs, as well as an unusual transcript of RNase P, and double-stranded RNAs. Subsequent computational analysis has revealed additional putative C/D-box snoRNAs. Our results will lead towards a future understanding of RNA metabolism in the deep-branching eukaryote Giardia, as more ncRNAs are characterized.

Noncoding RNAs of the H/ACA family

The H/ACA RNAs are an abundant family of trans-acting, noncoding RNAs found in eukaryotes and archaea. More than 100 H/ACA RNAs are known to exist in humans. The function of the majority of the identified H/ACA RNAs is to guide sites-pecific pseudouridylation of ribosomal RNA. In eukaryotes, H/ACA RNAs also mediate the processing of pre-rRNA, provide the template for telomere synthesis, and guide pseudouridylation of other classes of target RNAs (e.g., small nuclear RNAs [snRNAs]). Thus, currently, the H/ACA RNAs are known to be integrally involved in the production of both ribosomes and spliceosomes, and in the maintenance of chromosome integrity. In addition, dozens of H/ACA RNAs have been identified for which no function has yet been determined. The H/ACA RNAs select and present substrate molecules via base pairing. All H/ACA RNAs contain conserved sequence elements (box H and box ACA) and assemble with a core set of four proteins to form functional ribonucleoprotein complexes (RNPs). Mutations in key RNA and protein components of H/ACA RNPs result in dyskeratosis congenita, a serious multisystem genetic disease. Impressive progress has been made very recently in understanding the biogenesis, trafficking, and function of H/ACA RNPs.

Non-coding RNAs: lessons from the small nuclear and small nucleolar RNAs

Recent advances have fuelled rapid growth in our appreciation of the tremendous number, diversity and biological importance of non-coding (nc)RNAs. Because ncRNAs typically function as ribonucleoprotein (RNP) complexes and not as naked RNAs, understanding their biogenesis is crucial to comprehending their regulation and function. The small nuclear and small nucleolar RNPs are two well studied classes of ncRNPs with elaborate assembly and trafficking pathways that provide paradigms for understanding the biogenesis of other ncRNPs.

Characterization and functional analysis of a novel double-guide C/D box snoRNA in the fission yeast

Ribose methylation of eukaryotic rRNA is directed by box C/D small nucleolar RNAs (snoRNAs), which pinpoint the nucleotide to be methylated in specific position within the rRNA sequence. Here, we report the identification of a novel double-guide C/D box snoRNA termed snR88 that directs methylation of two previously undetermined sites in 25S rRNA from the fission yeast. Knockout of the predicted TATA box of the snR88 gene resulted in the complete blocking of its expression, showing that snR88 is an independently transcribed gene and dispensable for yeast viability. The depletion of snR88 abolished 25S rRNA methylation at U2304 and U2497 simultaneously. Interestingly, an unusual pause of reverse transcription at U2495 was observed, which implies an unknown structure of 25S rRNA related to ribose methylation at U2497 in the fission yeast.

Genome-wide analysis of C/D and H/ACA-like small nucleolar RNAs in Leishmania major indicates conservation among trypanosomatids in the repertoire and in their rRNA targets

Small nucleolar RNAs (snoRNAs) are a large group of noncoding RNAs that exist in eukaryotes and archaea and guide modifications such as 2'-O-ribose methylations and pseudouridylation on rRNAs and snRNAs. Recently, we described a genome-wide screening approach with Trypanosoma brucei that revealed over 90 guide RNAs. In this study, we extended this approach to analyze the repertoire of the closely related human pathogen Leishmania major. We describe 23 clusters that encode 62 C/Ds that can potentially guide 79 methylations and 37 H/ACA-like RNAs that can potentially guide 30 pseudouridylation reactions. Like T. brucei, Leishmania also contains many modifications and guide RNAs relative to its genome size. This study describes 10 H/ACAs and 14 C/Ds that were not found in T. brucei. Mapping of 2'-O-methylations in rRNA regions rich in modifications suggests the existence of trypanosomatid-specific modifications conserved in T. brucei and Leishmania. Structural features of C/D snoRNAs, such as copy number, conservation of boxes, K turns, and intragenic and extragenic base pairing, were examined to elucidate the great variation in snoRNA abundance. This study highlights the power of comparative genomics for determining conserved features of noncoding RNAs.

A novel class of developmentally regulated noncoding RNAs in Leishmania

Leishmania is a protozoan parasite that causes serious morbidity and mortality in humans worldwide. The ability of these parasites to survive within the phagolysosomes of mammalian macrophages is dependent on the developmental regulation of a variety of genes. Identifying genomic sequences that are preferentially expressed during the parasite's intracellular growth would provide new insights about the mechanisms controlling stage-specific gene regulation for intracellular development of the parasite. Using a genomic library that differentially hybridized to probes made from total RNA from Leishmania infantum amastigote or promastigote life cycle stages, we identified a new class of noncoding RNAs (ncRNAs) ranging from approximately 300 to 600 nucleotides in size that are expressed specifically in the intracellular amastigote stage. These ncRNAs are transcribed by RNA polymerase II from genomic clusters of tandem head-to-tail repeats, which are mainly located within subtelomeric regions. Remarkably, both the sense and antisense orientations of these ncRNAs are transcribed and are processed by trans splicing and polyadenylation. The levels of antisense transcripts are at least 10-fold lower than those of the sense transcripts and are tightly regulated. The sense and antisense ncRNAs are cytosolic as shown by fluorescence in situ hybridization studies and cosediment with a small ribonucleoprotein complex. Amastigote-specific regulation of these ncRNAs possibly occurs at the level of RNA stability. Interestingly, overexpression of these ncRNAs in promastigotes, as part of an episomal expression vector, failed to produce any transcript, which further highlights the instability of these RNAs in the promastigote stage. This is the first report describing developmentally regulated ncRNAs in protozoan parasites.

Characterization and differential nuclear localization of Nopp140 and a novel Nopp140-like protein in trypanosomes

Trypanosomatids possess two homologues of Nopp140: a canonical Nopp140 and a Nopp140-like protein (TbNoLP) in which a GAR domain replaces the C-terminal SRP40 domain. Both are phosphorylated and coimmunoprecipitate with RNA polymerase I. Each paralogue has a distinct subnuclear localization, and depletion of TbNoLP produces an enlarged nucleolus in which TbNopp140-containing regions disperse. The restricted occurrence pattern of NoLP proteins reflects an intriguing convergence in evolution, suggestive of a function in nucleoplasmic small nucleolar ribonucleoprotein shuttling.

A Large Collection of Compact Box C/D snoRNAs and their Isoforms in Euglena gracilis: Structural, Functional and Evolutionary Insights

In the domains Eucarya and Archaea, box C/D RNAs guide methylation at the 2'-position of selected ribose residues in ribosomal RNA (rRNA). Those eukaryotic box C/D RNAs that have been identified to date are larger and more variable in size than their archaeal counterparts. Here, we report the first extensive identification and characterization of box C/D small nucleolar (sno) RNAs from the protist Euglena gracilis. Among several unexpected findings, this organism contains a large assortment of methylation-guide RNAs that are smaller and more uniformly sized than those of other eukaryotes, and that consist of surprisingly few double-guide RNAs targeting sites of rRNA modification. Our comprehensive examination of the modification status of E.gracilis rRNA indicates that many of these box C/D snoRNAs target clustered methylation sites requiring extensive, overlapping guide RNA/rRNA pairings. An examination of the structure of the RNAs, in particular the location of the functional guide elements, suggests that the distances between adjacent box elements are an important factor in determining which of the potential guide elements is used to target a site of O(2')-methylation.

Elucidating the role of H/ACA-like RNAs in trans-splicing and rRNA processing via RNA interference silencing of the Trypanosoma brucei CBF5 pseudouridine synthase

Most pseudouridinylation in eukaryotic rRNA and small nuclear RNAs is guided by H/ACA small nucleolar RNAs. In this study, the Trypanosoma brucei pseudouridine synthase, Cbf5p, a snoRNP protein, was identified and silenced by RNAi. Depletion of this protein destabilized all small nucleolar RNAs of the H/ACA-like family. Following silencing, defects in rRNA processing, such as accumulation of precursors and inhibition of cleavages to generate the mature rRNA, were observed. snR30, an H/ACA RNA involved in rRNA maturation, was identified based on prototypical conserved domains characteristic of this RNA in other eukaryotes. The silencing of CBF5 also eliminated the spliced leader-associated (SLA1) RNA that directs pseudouridylation on the spliced leader RNA (SL RNA), which is the substrate for the trans-splicing reaction. Surprisingly, the depletion of Cbf5p not only eliminated the pseudouridine on the SL RNA but also abolished capping at the fourth cap-4 nucleotide. As a result of defects in the SL RNA and decreased modification on the U small nuclear RNA, trans-splicing was inhibited at the first step of the reaction, providing evidence for the essential role of H/ACA RNAs and the modifications they guide on trans-splicing.

A genome-wide analysis of C/D and H/ACA-like small nucleolar RNAs in Trypanosoma brucei reveals a trypanosome-specific pattern of rRNA modification

Small nucleolar RNAs (snoRNAs) constitute newly discovered noncoding small RNAs, most of which function in guiding modifications such as 2'-O-ribose methylation and pseudouridylation on rRNAs and snRNAs. To investigate the genome organization of Trypanosoma brucei snoRNAs and the pattern of rRNA modifications, we used a whole-genome approach to identify the repertoire of these guide RNAs. Twenty-one clusters encoding for 57 C/D snoRNAs and 34 H/ACA-like RNAs, which have the potential to direct 84 methylations and 32 pseudouridines, respectively, were identified. The number of 2'-O-methyls (Nms) identified on rRNA represent 80% of the expected modifications. The modifications guided by these RNAs suggest that trypanosomes contain many modifications and guide RNAs relative to their genome size. Interestingly, approximately 40% of the Nms are species-specific modifications that do not exist in yeast, humans, or plants, and 40% of the species-specific predicted modifications are located in unique positions outside the highly conserved domains. Although most of the guide RNAs were found in reiterated clusters, a few single-copy genes were identified. The large repertoire of modifications and guide RNAs in trypanosomes suggests that these modifications possibly play a central role in these parasites.

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.

Guide RNAs with 5' caps and novel box C/D snoRNA-like domains for modification of snRNAs in metazoa

BACKGROUND

Spliceosomal snRNAs and ribosomal RNAs in metazoans contain numerous modified residues that are functionally important. The most common modifications are site-specific 2'-O-methylation and pseudouridylation, both directed by small ribonucleoprotein particles. Each particle is composed of a short guide RNA and a set of several proteins. All previously characterized modification guide RNAs in metazoa are encoded in and processed from introns.

RESULTS

We have identified and characterized three novel guide RNAs for conserved 2'-O-methylation of U2, U4, and U12 snRNAs. Two guides, termed mgU2-25/61 and mgU12-22/U4-8, appear to be independently transcribed as judged by the presence of methylated guanosine caps at their 5' ends and upstream promoters similar to those of telomerase RNA. These guide RNAs are each composed of a canonical box C/D snoRNA and a novel box C/D snoRNA-like domain, where the C'/D' motif, rather than C/D, can be folded into a conserved kink-turn structure. The snoRNA-like domains are predicted to direct 2'-O-methylation of invariant G residues that occupy analogous positions in the U2 and U12 snRNA secondary structures. A third guide, mgU2-19/30 RNA, is composed of two canonical box C/D snoRNA domains encoded within a single intron.

CONCLUSIONS

This is the first description in metazoan cells of 5'-capped modification guide RNAs that appear to be independently transcribed. Since plant, yeast, and protozoan guide RNAs are mostly independently transcribed, the identification of such RNAs argues that ancestral metazoans possessed independently transcribed guide RNAs and only later, during the evolution of metazoan organisms, did the guide RNA genes shift to introns.

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.