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.

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.

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.