Psiscan: a computational approach to identify H/ACA-like and AGA-like non-coding RNA in trypanosomatid genomes

Discovery of 17 conserved structural RNAs in fungi

Many non-coding RNAs with known functions are structurally conserved: their intramolecular secondary and tertiary interactions are maintained across evolutionary time. Consequently, the presence of conserved structure in multiple sequence alignments can be used to identify candidate functional non-coding RNAs. Here, we present a bioinformatics method that couples iterative homology search with covariation analysis to assess whether a genomic region has evidence of conserved RNA structure. We used this method to examine all unannotated regions of five well-studied fungal genomes (Saccharomyces cerevisiae, Candida albicans, Neurospora crassa, Aspergillus fumigatus, and Schizosaccharomyces pombe). We identified 17 novel structurally conserved non-coding RNA candidates, which include four H/ACA box small nucleolar RNAs, four intergenic RNAs and nine RNA structures located within the introns and untranslated regions (UTRs) of mRNAs. For the two structures in the 3' UTRs of the metabolic genes GLY1 and MET13, we performed experiments that provide evidence against them being eukaryotic riboswitches.

Computer aided drug discovery review for infectious diseases with case study of anti-Chagas project

Neglected tropical diseases (NTDs) are parasitic and bacterial infections that are widespread, especially in the tropics, and cause health problems for about one billion people over 149 countries worldwide. However, in terms of therapeutic agents, for example, nifurtimox and benznidazole were developed in the 1960s to treat Chagas disease, but new drugs are desirable because of their side effects. Drug discovery takes 12 to 14 years and costs $2.6 billon dollars, and hence, computer aided drug discovery (CADD) technology is expected to reduce the time and cost. This paper describes our methods and results based on CADD, mainly for NTDs. An overview of databases, molecular simulation and pharmacophore modeling, contest-based drug discovery, and machine learning and their results are presented herein.

Small nucleolar RNAs: continuing identification of novel members and increasing diversity of their molecular mechanisms of action

Identified five decades ago amongst the most abundant cellular RNAs, small nucleolar RNAs (snoRNAs) were initially described as serving as guides for the methylation and pseudouridylation of ribosomal RNA through direct base pairing. In recent years, however, increasingly powerful high-throughput genomic approaches and strategies have led to the discovery of many new members of the family and surprising diversity in snoRNA functionality and mechanisms of action. SnoRNAs are now known to target RNAs of many biotypes for a wider range of modifications, interact with diverse binding partners, compete with other binders for functional interactions, recruit diverse players to targets and affect protein function and accessibility through direct interaction. This mini-review presents the continuing characterization of the snoRNome through the identification of new snoRNA members and the discovery of their mechanisms of action, revealing a highly versatile noncoding family playing central regulatory roles and connecting the main cellular processes.

The large repertoire of 2'-O-methylation guided by C/D snoRNAs on Trypanosoma brucei rRNA

The parasite Trypanosoma brucei cycles between insect and mammalian hosts, and is the causative agent of sleeping sickness. Here, we performed genome-wide mapping of 2'-O-methylations (Nms) on trypanosome rRNA using three high-throughput sequencing methods; RibOxi-seq, RiboMeth-seq and 2'-OMe-seq. This is the first study using three genome-wide mapping approaches on rRNA from the same species showing the discrepancy among the methods. RibOxi-seq detects all the sites, but RiboMeth-seq is the only method to evaluate the level of a single Nm site. The sequencing revealed at least ninety-nine Nms guided by eighty-five snoRNAs among these thirty-eight Nms are trypanosome specific sites. We present the sequence and target of the C/D snoRNAs guiding on rRNA. This is the highest number of Nms detected to date on rRNA of a single cell parasite. Based on RiboMeth-seq, several Nm sites were found to be differentially regulated at the two stages of the parasite life cycle, the insect procyclic form (PCF) versus the bloodstream form (BSF) in the mammalian host.

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.

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.

Characterization of RNase MRP RNA and novel snoRNAs from Giardia intestinalis and Trichomonas vaginalis

Background

Eukaryotic cells possess a complex network of RNA machineries which function in RNA-processing and cellular regulation which includes transcription, translation, silencing, editing and epigenetic control. Studies of model organisms have shown that many ncRNAs of the RNA-infrastructure are highly conserved, but little is known from non-model protists. In this study we have conducted a genome-scale survey of medium-length ncRNAs from the protozoan parasites Giardia intestinalis and Trichomonas vaginalis.

Results

We have identified the previously 'missing' Giardia RNase MRP RNA, which is a key ribozyme involved in pre-rRNA processing. We have also uncovered 18 new H/ACA box snoRNAs, expanding our knowledge of the H/ACA family of snoRNAs.

Conclusions

Results indicate that Giardia intestinalis and Trichomonas vaginalis, like their distant multicellular relatives, contain a rich infrastructure of RNA-based processing. From here we can investigate the evolution of RNA processing networks in eukaryotes.

RNA-seq analysis of small RNPs in Trypanosoma brucei reveals a rich repertoire of non-coding RNAs

The discovery of a plethora of small non-coding RNAs (ncRNAs) has fundamentally changed our understanding of how genes are regulated. In this study, we employed the power of deep sequencing of RNA (RNA-seq) to examine the repertoire of ncRNAs present in small ribonucleoprotein particles (RNPs) of Trypanosoma brucei, an important protozoan parasite. We identified new C/D and H/ACA small nucleolar RNAs (snoRNAs), as well as tens of putative novel non-coding RNAs; several of these are processed from trans-spliced and polyadenylated transcripts. The RNA-seq analysis provided information on the relative abundance of the RNAs, and their 5'- and 3'-termini. The study demonstrated that three highly abundant snoRNAs are involved in rRNA processing and highlight the unique trypanosome-specific repertoire of these RNAs. Novel RNAs were studied using in situ hybridization, association in RNP complexes, and 'RNA walk' to detect interaction with their target RNAs. Finally, we showed that the abundance of certain ncRNAs varies between the two stages of the parasite, suggesting that ncRNAs may contribute to gene regulation during the complex parasite's life cycle. This is the first study to provide a whole-genome analysis of the large repertoire of small RNPs in trypanosomes.

Cell biology of the trypanosome genome

Trypanosomes are a group of protozoan eukaryotes, many of which are major parasites of humans and livestock. The genomes of trypanosomes and their modes of gene expression differ in several important aspects from those of other eukaryotic model organisms. Protein-coding genes are organized in large directional gene clusters on a genome-wide scale, and their polycistronic transcription is not generally regulated at initiation. Transcripts from these polycistrons are processed by global trans-splicing of pre-mRNA. Furthermore, in African trypanosomes, some protein-coding genes are transcribed by a multifunctional RNA polymerase I from a specialized extranucleolar compartment. The primary DNA sequence of the trypanosome genomes and their cellular organization have usually been treated as separate entities. However, it is becoming increasingly clear that in order to understand how a genome functions in a living cell, we will need to unravel how the one-dimensional genomic sequence and its trans-acting factors are arranged in the three-dimensional space of the eukaryotic nucleus. Understanding this cell biology of the genome will be crucial if we are to elucidate the genetic control mechanisms of parasitism. Here, we integrate the concepts of nuclear architecture, deduced largely from studies of yeast and mammalian nuclei, with recent developments in our knowledge of the trypanosome genome, gene expression, and nuclear organization. We also compare this nuclear organization to those in other systems in order to shed light on the evolution of nuclear architecture in eukaryotes.

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.

Identification and classification of ncRNA molecules using graph properties

The study of non-coding RNA genes has received increased attention in recent years fuelled by accumulating evidence that larger portions of genomes than previously acknowledged are transcribed into RNA molecules of mostly unknown function, as well as the discovery of novel non-coding RNA types and functional RNA elements. Here, we demonstrate that specific properties of graphs that represent the predicted RNA secondary structure reflect functional information. We introduce a computational algorithm and an associated web-based tool (GraPPLE) for classifying non-coding RNA molecules as functional and, furthermore, into Rfam families based on their graph properties. Unlike sequence-similarity-based methods and covariance models, GraPPLE is demonstrated to be more robust with regard to increasing sequence divergence, and when combined with existing methods, leads to a significant improvement of prediction accuracy. Furthermore, graph properties identified as most informative are shown to provide an understanding as to what particular structural features render RNA molecules functional. Thus, GraPPLE may offer a valuable computational filtering tool to identify potentially interesting RNA molecules among large candidate datasets.