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

Identification and functional implications of pseudouridine RNA modification on small noncoding RNAs in the mammalian pathogen Trypanosoma brucei

Trypanosoma brucei, the parasite that causes sleeping sickness, cycles between an insect and a mammalian host. However, the effect of RNA modifications such as pseudouridinylation on its ability to survive in these two different host environments is unclear. Here, two genome-wide approaches were applied for mapping pseudouridinylation sites (Ψs) on small nucleolar RNA (snoRNA), 7SL RNA, vault RNA, and tRNAs from T. brucei. We show using HydraPsiSeq and RiboMeth-seq that the Ψ on C/D snoRNA guiding 2'-O-methylation increased the efficiency of the guided modification on its target, rRNA. We found differential levels of Ψs on these noncoding RNAs in the two life stages (insect host and mammalian host) of the parasite. Furthermore, tRNA isoform abundance and Ψ modifications were characterized in these two life stages demonstrating stage-specific regulation. We conclude that the differential Ψ modifications identified here may contribute to modulating the function of noncoding RNAs involved in rRNA processing, rRNA modification, protein synthesis, and protein translocation during cycling of the parasite between its two hosts.

Developmentally Regulated Novel Non-coding Anti-sense Regulators of mRNA Translation in Trypanosoma b rucei

The parasite Trypanosoma brucei is the causative agent of sleeping sickness and cycles between insect and mammalian hosts. The parasite appears to lack conventional transcriptional regulation of protein coding genes, and mRNAs are processed from polycistronic transcripts by the concerted action of trans-splicing and polyadenylation. Regulation of mRNA function is mediated mainly by RNA binding proteins affecting mRNA stability and translation. In this study, we describe the identification of 62 non-coding (nc) RNAs that are developmentally regulated and/or respond to stress. We characterized two novel anti-sense RNA regulators (TBsRNA-33 and 37) that originate from the rRNA loci, associate with ribosomes and polyribosomes, and interact in vivo with distinct mRNA species to regulate translation. Thus, this study suggests for the first-time anti-sense RNA regulators as an additional layer for controlling gene expression in these parasites.

•

Trypanosome non-coding RNAs (ncRNAs) are developmentally regulated during cycling between two hosts

•

ncRNAs originate from rRNA locus and associate with the ribosome en route to cytoplasm

•

In vivo cross-linking enable identification of target RNA species regulated by ncRNAs

•

Trypanosomes possess anti-sense ncRNAs that regulate translation

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.

Nucleolar Structure and Function in Trypanosomatid Protozoa

The nucleolus is the conspicuous nuclear body where ribosomal RNA genes are transcribed by RNA polymerase I, pre-ribosomal RNA is processed, and ribosomal subunits are assembled. Other important functions have been attributed to the nucleolus over the years. Here we review the current knowledge about the structure and function of the nucleolus in the trypanosomatid parasites Trypanosoma brucei, Trypanosoma cruzi and Leishmania ssp., which represent one of the earliest branching lineages among the eukaryotes. These protozoan parasites present a single nucleolus that is preserved throughout the closed nuclear division, and that seems to lack fibrillar centers. Trypanosomatids possess a relatively low number of rRNA genes, which encode rRNA molecules that contain large expansion segments, including several that are trypanosomatid-specific. Notably, the large subunit rRNA (28S-type) is fragmented into two large and four small rRNA species. Hence, compared to other organisms, the rRNA primary transcript requires additional processing steps in trypanosomatids. Accordingly, this group of parasites contains the highest number ever reported of snoRNAs that participate in rRNA processing. The number of modified rRNA nucleotides in trypanosomatids is also higher than in other organisms. Regarding the structure and biogenesis of the ribosomes, recent cryo-electron microscopy analyses have revealed several trypanosomatid-specific features that are discussed here. Additional functions of the nucleolus in trypanosomatids are also reviewed.

Small nucleolar RNAs controlling rRNA processing in Trypanosoma brucei

In trypanosomes, in contrast to most eukaryotes, the large subunit (LSU) ribosomal RNA is fragmented into two large and four small ribosomal RNAs (srRNAs) pieces, and this additional processing likely requires trypanosome-specific factors. Here, we examined the role of 10 abundant small nucleolar RNAs (snoRNAs) involved in rRNA processing. We show that each snoRNA involved in LSU processing associates with factors engaged in either early or late biogenesis steps. Five of these snoRNAs interact with the intervening sequences of rRNA precursor, whereas the others only guide rRNA modifications. The function of the snoRNAs was explored by silencing snoRNAs. The data suggest that the LSU rRNA processing events do not correspond to the order of rRNA transcription, and that srRNAs 2, 4 and 6 which are part of LSU are processed before srRNA1. Interestingly, the 6 snoRNAs that affect srRNA1 processing guide modifications on rRNA positions that span locations from the protein exit tunnel to the srRNA1, suggesting that these modifications may serve as check-points preceding the liberation of srRNA1. This study identifies the highest number of snoRNAs so far described that are involved in rRNA processing and/or rRNA folding and highlights their function in the unique trypanosome rRNA maturation events.

TbUTP10, a protein involved in early stages of pre-18S rRNA processing in Trypanosoma brucei

Ribosome biosynthesis, best studied in opisthokonts, is a highly complex process involving numerous protein and RNA factors. Yet, very little is known about the early stages of pre-18S rRNA processing even in these model organisms, let alone the conservation of this mechanism in other eukaryotes. Here we extend our knowledge of this process by identifying and characterizing the essential protein TbUTP10, a homolog of yeast U3 small nucleolar RNA-associated protein 10 - UTP10 (HEATR1 in human), in the excavate parasitic protist Trypanosoma brucei. We show that TbUTP10 localizes to the nucleolus and that its ablation by RNAi knock-down in two different T. brucei life cycle stages results in similar phenotypes: a disruption of pre-18S rRNA processing, exemplified by the accumulation of rRNA precursors, a reduction of mature 18S rRNA, and also a decrease in the level of U3 snoRNA. Moreover, polysome profiles of the RNAi-induced knock-down cells show a complete disappearance of the 40S ribosomal subunit, and a prominent accumulation of the 60S large ribosomal subunit, reflecting impaired ribosome assembly. Thus, TbUTP10 is an important protein in the processing of 18S rRNA.

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.

The Evolutionary Loss of RNAi Key Determinants in Kinetoplastids as a Multiple Sporadic Phenomenon

We screened the genomes of a broad panel of kinetoplastid protists for genes encoding proteins associated with the RNA interference (RNAi) system using probes from the Argonaute (AGO1), Dicer1 (DCL1), and Dicer2 (DCL2) genes of Leishmania brasiliensis and Crithidia fasciculata. We identified homologs for all the three of these genes in the genomes of a subset of these organisms. However, several of these organisms lacked evidence for any of these genes, while others lacked only DCL2. The open reading frames encoding these putative proteins were structurally analyzed in silico. The alignments indicated that the genes are homologous with a high degree of confidence, and three-dimensional structural models strongly supported a functional relationship to previously characterized AGO1, DCL1, and DCL2 proteins. Phylogenetic analysis of these putative proteins showed that these genes, when present, evolved in parallel with other nuclear genes, arguing that the RNAi system genes share a common progenitor, likely across all Kinetoplastea. In addition, the genome segments bearing these genes are highly conserved and syntenic, even among those taxa in which they are absent. However, taxa in which these genes are apparently absent represent several widely divergent branches of kinetoplastids, arguing that these genes were independently lost at least six times in the evolutionary history of these organisms. The mechanisms responsible for the apparent coordinate loss of these RNAi system genes independently in several lineages of kinetoplastids, while being maintained in other related lineages, are currently unknown.

Transcriptome and proteome analyses and the role of atypical calpain protein and autophagy in the spliced leader silencing pathway in Trypanosoma brucei

Under persistent ER stress, Trypanosoma brucei parasites induce the spliced leader silencing (SLS) pathway. In SLS, transcription of the SL RNA gene, the SL donor to all mRNAs, is extinguished, arresting trans-splicing and leading to programmed cell death (PCD). In this study, we investigated the transcriptome following silencing of SEC63, a factor essential for protein translocation across the ER membrane, and whose silencing induces SLS. The proteome of SEC63-silenced cells was analyzed with an emphasis on SLS-specific alterations in protein expression, and modifications that do not directly result from perturbations in trans-splicing. One such protein identified is an atypical calpain SKCRP7.1/7.2. Co-silencing of SKCRP7.1/7.2 and SEC63 eliminated SLS induction due its role in translocating the PK3 kinase. This kinase initiates SLS by migrating to the nucleus and phosphorylating TRF4 leading to shut-off of SL RNA transcription. Thus, SKCRP7.1 is involved in SLS signaling and the accompanying PCD. The role of autophagy in SLS was also investigated; eliminating autophagy through VPS34 or ATG7 silencing demonstrated that autophagy is not essential for SLS induction, but is associated with PCD. Thus, this study identified factors that are used by the parasite to cope with ER stress and to induce SLS and PCD.

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.

DNA Intercalated Psoralen Undergoes Efficient Photoinduced Electron Transfer

The interaction of psoralens with DNA has been used for therapeutic and research purposes for decades. Still the photoinduced behavior of psoralens in DNA has never been observed directly. Femtosecond transient absorption spectroscopy is used here to gain direct insight into the photophysics of a DNA-intercalated psoralen (4'-aminomethyl-4,5',8-trimethyl-psoralen (AMT)). Intercalation reduces the excited singlet lifetime of AMT to 4 ps compared with 1400 ps for AMT in water. This singlet quenching prohibits the population of the triplet state that is accessed in free AMT. Instead, a DNA to AMT electron transfer takes place. The resulting radical pair decays primarily via charge recombination with a time constant of 30 ps. The efficient electron transfer observed here reveals a completely new aspect of the psoralen-DNA interaction.

Ribosomal RNA gene transcription in trypanosomes

Leishmania major, Trypanosoma cruzi and Trypanosoma brucei are pathogenic species from the order Kinetoplastida. The molecular and cellular studies of parasites, such as of the biosynthesis of essential macromolecules, are important in designing successful strategies for control. A major stage in ribosome biogenesis is the transcription of genes encoding ribosomal (r)RNA. These genes are transcribed in trypanosome cells by RNA polymerase I, similar to what occurs in all eukaryotes analysed to date. In addition, and most remarkably, the African species, T. brucei, transcribe their major cell surface protein genes using this class of polymerase. Since its discovery, the research interest in this phenomenon has been overwhelming; therefore, analysis of the canonical, yet essential, transcription of rRNA has been comparatively neglected. In this work, a review of rRNA gene transcription and data on gene promoter structures, transcription machineries and epigenetic conditions is presented for trypanosomatids. Because species-specific molecules represent potential targets for chemotherapy, their existence within trypanosomes is highlighted.

Basal splicing factors regulate the stability of mature mRNAs in trypanosomes

Gene expression in trypanosomes is mainly regulated post-transcriptionally. Genes are transcribed as polycistronic mRNAs that are dissected by the concerted action of trans-splicing and polyadenylation. In trans-splicing, a common exon, the spliced leader, is added to all mRNAs from a small RNA. In this study, we examined by microarray analysis the transcriptome following RNAi silencing of the basal splicing factors U2AF65, SF1, and U2AF35. The transcriptome data revealed correlations between the affected genes and their splicing and polyadenylation signaling properties, suggesting that differential binding of these factors to pre-mRNA regulates trans-splicing and hence expression of specific genes. Surprisingly, all these factors were shown to affect not only splicing but also mRNA stability. Affinity purification of SF1 and U2AF35 complexes supported their role in mRNA stability. U2AF35 but not SF1 was shown to bind to ribosomes. To examine the role of splicing factors in mRNA stability, mutations were introduced into the polypyrimidine tract located in the 3' UTR of a mini-gene, and the results demonstrate that U2AF65 binds to such a site and controls the mRNA stability. We propose that transcripts carrying splicing signals in their 3' UTR bind the splicing factors and control their stability.

Comparative genomics reveals two novel RNAi factors in Trypanosoma brucei and provides insight into the core machinery

The introduction ten years ago of RNA interference (RNAi) as a tool for molecular exploration in Trypanosoma brucei has led to a surge in our understanding of the pathogenesis and biology of this human parasite. In particular, a genome-wide RNAi screen has recently been combined with next-generation Illumina sequencing to expose catalogues of genes associated with loss of fitness in distinct developmental stages. At present, this technology is restricted to RNAi-positive protozoan parasites, which excludes T. cruzi, Leishmania major, and Plasmodium falciparum. Therefore, elucidating the mechanism of RNAi and identifying the essential components of the pathway is fundamental for improving RNAi efficiency in T. brucei and for transferring the RNAi tool to RNAi-deficient pathogens. Here we used comparative genomics of RNAi-positive and -negative trypanosomatid protozoans to identify the repertoire of factors in T. brucei. In addition to the previously characterized Argonaute 1 (AGO1) protein and the cytoplasmic and nuclear Dicers, TbDCL1 and TbDCL2, respectively, we identified the RNA Interference Factors 4 and 5 (TbRIF4 and TbRIF5). TbRIF4 is a 3′-5′ exonuclease of the DnaQ superfamily and plays a critical role in the conversion of duplex siRNAs to the single-stranded form, thus generating a TbAGO1-siRNA complex required for target-specific cleavage. TbRIF5 is essential for cytoplasmic RNAi and appears to act as a TbDCL1 cofactor. The availability of the core RNAi machinery in T. brucei provides a platform to gain mechanistic insights in this ancient eukaryote and to identify the minimal set of components required to reconstitute RNAi in RNAi-deficient parasites.

RNA interference (RNAi), a naturally-occurring pathway whereby the presence of double-stranded RNA in a cell triggers the degradation of homologous mRNA, has been harnessed in many organisms as an invaluable molecular biology tool to interrogate gene function. Although this technology is widely used in the protozoan parasite Trypanosoma brucei, other parasites of considerable public health significance, such as Trypanosoma cruzi, Leishmania major, and Plasmodium falciparum do not perform RNAi. Since RNAi has recently been introduced into budding yeast, this opens up the possibility that RNAi can be reconstituted in these pathogens. The key to this is getting a handle on the essential RNAi factors in T. brucei. By applying comparative genomics we identified five genes that are present in the RNAi-proficient species, but not in RNAi-deficient species: three previously identified RNAi factors, and two novel ones, which are described here. This insight into the core T. brucei RNAi machinery represents a major step towards transferring this pathway to RNAi-deficient parasites.

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.

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.

Pseudogene-derived small interference RNAs regulate gene expression in African Trypanosoma brucei

Pseudogenes have been shown to acquire unique regulatory roles from more and more organisms. We report the observation of a cluster of siRNAs derived from pseudogenes of African Trypanosoma brucei using high through-put analysis. We show that these pseudogene-derived siRNAs suppress gene expression through RNA interference. The discovery that siRNAs may originate from pseudogenes and regulate gene expression in a unicellular eukaryote provides insights into the functional roles of pseudogenes and into the origin of noncoding small RNAs.

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.