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

Transcript tinkering: RNA modifications in protozoan parasites

Apicomplexan and trypanosomatid parasites have evolved a wide range of post-transcriptional processes that allow them to replicate, differentiate, and transmit within and among multiple different tissue, host, and vector environments. In this review, we highlight the recent advances that point toward the regulatory potential of RNA modifications in mediating these processes on the coding and noncoding transcriptome throughout the life cycle of protozoan parasites. We discuss the recent technical advancements enabling the study of the 'epitranscriptome' and how parasites evolved RNA modification-mediated mechanisms adapted to their unique lifestyles.

Emerging roles of the epitranscriptome in parasitic protozoan biology and pathogenesis

RNA modifications (epitranscriptome) - such as N6-methyladenosine (m6A), 5-methylcytosine (m5C), and pseudouridine (Ψ) - modulate RNA processing, stability, interaction, and translation, thereby playing critical roles in the development, replication, virulence, metabolism, and life cycle adaptations of parasitic protozoa. Here, we summarize potential homologs of the major human RNA modification regulatory factors in parasites, outline current knowledge on how RNA modifications affect parasitic protozoa, highlight the regulation of RNA modifications and their crosstalk, and discuss current progress in exploring RNA modifications as potential drug targets. This review contributes to our understanding of epitranscriptomic regulation of parasitic protozoa biology and pathogenesis and provides new perspectives for the treatment of parasitic diseases.

UMSBP2 is chromatin remodeler that functions in regulation of gene expression and suppression of antigenic variation in trypanosomes

Universal Minicircle Sequence binding proteins (UMSBPs) are CCHC-type zinc-finger proteins that bind the single-stranded G-rich UMS sequence, conserved at the replication origins of minicircles in the kinetoplast DNA, the mitochondrial genome of kinetoplastids. Trypanosoma brucei UMSBP2 has been recently shown to colocalize with telomeres and to play an essential role in chromosome end protection. Here we report that TbUMSBP2 decondenses in vitro DNA molecules, which were condensed by core histones H2B, H4 or linker histone H1. DNA decondensation is mediated via protein-protein interactions between TbUMSBP2 and these histones, independently of its previously described DNA binding activity. Silencing of the TbUMSBP2 gene resulted in a significant decrease in the disassembly of nucleosomes in T. brucei chromatin, a phenotype that could be reverted, by supplementing the knockdown cells with TbUMSBP2. Transcriptome analysis revealed that silencing of TbUMSBP2 affects the expression of multiple genes in T. brucei, with a most significant effect on the upregulation of the subtelomeric variant surface glycoproteins (VSG) genes, which mediate the antigenic variation in African trypanosomes. These observations suggest that UMSBP2 is a chromatin remodeling protein that functions in the regulation of gene expression and plays a role in the control of antigenic variation in T. brucei.

Proteomic analysis defines the interactome of telomerase in the protozoan parasite, Trypanosoma brucei

Telomerase is a ribonucleoprotein enzyme responsible for maintaining the telomeric end of the chromosome. The telomerase enzyme requires two main components to function: the telomerase reverse transcriptase (TERT) and the telomerase RNA (TR), which provides the template for telomeric DNA synthesis. TR is a long non-coding RNA, which forms the basis of a large structural scaffold upon which many accessory proteins can bind and form the complete telomerase holoenzyme. These accessory protein interactions are required for telomerase activity and regulation inside cells. The interacting partners of TERT have been well studied in yeast, human, and Tetrahymena models, but not in parasitic protozoa, including clinically relevant human parasites. Here, using the protozoan parasite, Trypanosoma brucei (T. brucei) as a model, we have identified the interactome of T. brucei TERT (TbTERT) using a mass spectrometry-based approach. We identified previously known and unknown interacting factors of TbTERT, highlighting unique features of T. brucei telomerase biology. These unique interactions with TbTERT, suggest mechanistic differences in telomere maintenance between T. brucei and other eukaryotes.

Dissecting Trypanosoma brucei RRP44 function in the maturation of segmented ribosomal RNA using a regulated genetic complementation system

Trypanosoma brucei belongs to a group of protozoans presenting fragmented large subunit rRNA. Its LSU rRNA equivalent to the 25S/28S rRNA of other eukaryotes is split into six fragments, requiring additional processing for removal of the extra spacer sequences. We have used a genetic complementation strategy to further investigate the T. brucei RRP44 nuclease in pre-rRNA maturation. TbRRP44 contains both a PIN and a RNB domain whose homologues are found in association with the exosome complex. We found that the exonucleolytic activity of the RNB domain as well as the physical presence of the PIN domain are essential for TbRRP44 function, while a catalytic site mutation in the PIN domain has no detectable effect on cell growth. A new endonucleolytic cleavage site in ITS1 was identified. In addition to the 5.8S rRNA 3'-end maturation, TbRRP44 is required for degradation of the excised 5'-ETS and for removal of part of ITS1 during maturation of the 18S rRNA 3'-end. TbRRP44 deficiency leads to accumulation of many LSU intermediate precursors, most of them not detected in control cells. TbRRP44 is also required for U3 snoRNA and spliced leader processing, indicating that TbRRP44 may have a wide role in RNA processing in T. brucei.

Current Status of Regulatory Non-Coding RNAs Research in the Tritryp

Trypanosomatids are protozoan parasites that cause devastating vector-borne human diseases. Gene expression regulation of these organisms depends on post-transcriptional control in responding to diverse environments while going through multiple developmental stages of their complex life cycles. In this scenario, non-coding RNAs (ncRNAs) are excellent candidates for a very efficient, quick, and economic strategy to regulate gene expression. The advent of high throughput RNA sequencing technologies show the presence and deregulation of small RNA fragments derived from canonical ncRNAs. This review seeks to depict the ncRNA landscape in trypanosomatids, focusing on the small RNA fragments derived from functional RNA molecules observed in RNA sequencing studies. Small RNA fragments derived from canonical ncRNAs (tsRNAs, snsRNAs, sdRNAs, and sdrRNAs) were identified in trypanosomatids. Some of these RNAs display changes in their levels associated with different environments and developmental stages, demanding further studies to determine their functional characterization and potential roles. Nevertheless, a comprehensive and detailed ncRNA annotation for most trypanosomatid genomes is still needed, allowing better and more extensive comparative and functional studies.

Fibrillarin RNA methylase is an interacting protein of Cryptosporidium parvum calmodulin-like protein (CpCML)

Cryptosporidium parvum is an obligate protozoan parasite invading epithelial cells of small intestine of human and animals, and causing diarrheal disease. In apicomplexan parasites, calcium signaling can regulate many essential biological processes such as invasion and migration. As the main intracellular receptor for calcium ions, calmodulins control the activities of hundreds of enzymes and proteins. Calmodulin-like protein (CML) is an important member of the calmodulin family and may play a key role in C. parvum, however, the actual situation is still not clear. The present study aimed to identify the parasite interaction partner proteins of C. parvum calmodulin-like protein (CpCML). By constructing the cpcml bait plasmid, 5 potential CpCML - interacting proteins in C. parvum oocyst were screened by yeast-two-hybrid system (Y2H). Bimolecular fluorescence complementation (BiFC) and Co-immunoprecipitation (Co-IP) were performed as subsequent validations. Fibrillarin RNA methylase (FBL) was identified via this screening method as CpCML interacting protein in C. parvum. The identification of this interaction made it possible to get a further understanding of the function of CpCML and its contribution to the pathogenicity of C. parvum.

Pedal to the Metal: Nuclear Splicing Bodies Turbo-Charge VSG mRNA Production in African Trypanosomes

The African trypanosome Trypanosoma brucei is a parasite of the mammalian bloodstream and tissues, where an antigenically variable Variant Surface Glycoprotein (VSG) coat protects it from immune attack. This dense layer comprised of ∼107 VSG proteins, makes VSG by far the most abundant mRNA (7-10% total) and protein (∼10% total) in the bloodstream form trypanosome. How can such prodigious amounts of VSG be produced from a single VSG gene? Extremely high levels of RNA polymerase I (Pol I) transcription of the active VSG provide part of the explanation. However, recent discoveries highlight the role of pre-mRNA processing, both in maintaining high levels of VSG transcription, as well as its monoallelic expression. Trypanosome mRNAs are matured through trans-splicing a spliced leader (SL) RNA to the 5' end of precursor transcripts, meaning abundant SL RNA is required throughout the nucleus. However, requirement for SL RNA in the vicinity of the active VSG gene is so intense, that the cell reconfigures its chromatin architecture to facilitate interaction between the SL RNA genes and the active VSG. This presumably ensures that sufficient localised SL RNA is available, and not limiting for VSG mRNA expression. Recently, novel nuclear splicing bodies which appear to provide essential trans-splicing components, have been identified associating with the active VSG. These observations highlight the underappreciated role of pre-mRNA processing in modulating gene expression in trypanosomes. Dissecting the function of these nuclear RNA processing bodies should help us elucidate the mechanisms of both VSG expression and monoallelic exclusion in T. brucei.

The enigmatic RNase MRP of kinetoplastids

The ribonucleoprotein RNase MRP is responsible for the processing of ribosomal RNA precursors. It is found in virtually all eukaryotes that have been examined. In the Euglenozoa, including the genera Euglena, Diplonema and kinetoplastids, MRP RNA and protein subunits have so far escaped detection using bioinformatic methods. However, we now demonstrate that the RNA component is widespread among the Euglenozoa and that these RNAs have secondary structures that conform to the structure of all other phylogenetic groups. In Euglena, we identified the same set of P/MRP protein subunits as in many other protists. However, we failed to identify any of these proteins in the kinetoplastids. This finding poses interesting questions regarding the structure and function of RNase MRP in these species.

Cryo-EM structure of the highly atypical cytoplasmic ribosome of Euglena gracilis

Ribosomal RNA is the central component of the ribosome, mediating its functional and architectural properties. Here, we report the cryo-EM structure of a highly divergent cytoplasmic ribosome from the single-celled eukaryotic alga Euglena gracilis. The Euglena large ribosomal subunit is distinct in that it contains 14 discrete rRNA fragments that are assembled non-covalently into the canonical ribosome structure. The rRNA is substantially enriched in post-transcriptional modifications that are spread far beyond the catalytic RNA core, contributing to the stabilization of this highly fragmented ribosome species. A unique cluster of five adenosine base methylations is found in an expansion segment adjacent to the protein exit tunnel, such that it is positioned for interaction with the nascent peptide. As well as featuring distinctive rRNA expansion segments, the Euglena ribosome contains four novel ribosomal proteins, localized to the ribosome surface, three of which do not have orthologs in other eukaryotes.

Differentially Expressed Homologous Genes Reveal Interspecies Differences of Paragonimus Proliferus based on Transcriptome Analysis

Paragonimus proliferus (P. proliferus), one of 46 Paragonimus species registered in the National Center for Biotechnology Information database, may be much more widely distributed in Southeast Asia than previously thought, as its reported natural foci have increased in the past decades. However, very little is known about its molecular biology, especially at the transcriptome level. For the first time, the transcriptome of this species was sequenced and compared with four other common Paragonimus species, namely Paragonimus skrjabini, Paragonimus kellicotti, Paragonimus miyazakii, and Paragonimus westermani, to predict homologous genes and differentially expressed homologous genes to explore interspecies differences of Paragonimus proliferus. A total of 7393 genes were found to be significantly differentially expressed. Of these, 49 were considered to be core genes because they were differentially expressed in all four comparison groups. Annotations revealed that these genes were related mainly to "duplication, transcription, or translation", energy or nutrient metabolism, and parasitic growth, proliferation, motility, invasion, adaptation to the host, or virulence. Interestingly, a majority (5601/7393) of the identified genes, and in particular the core genes (48/49), were expressed at lower levels in P. proliferus. The identified genes may play essential roles in the biological differences between Paragonimus species. This work provides fundamental background information for further research into the molecular biology of P. proliferus.

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.

Involvement in surface antigen expression by a moonlighting FG-repeat nucleoporin in trypanosomes

The nuclear pore complex is an ancient component of the eukaryotic cell. We show here that an FG nucleoporin, TbNup53b, in trypanosomes has an association with the splicing machinery and roles in gene expression, indicating that moonlighting roles for nucleoporins are highly ancient and present in the earliest eukaryotes.

Components of the nuclear periphery coordinate a multitude of activities, including macromolecular transport, cell-cycle progression, and chromatin organization. Nuclear pore complexes (NPCs) mediate nucleocytoplasmic transport, mRNA processing, and transcriptional regulation, and NPC components can define regions of high transcriptional activity in some organisms at the nuclear periphery and nucleoplasm. Lineage-specific features underpin several core nuclear functions and in trypanosomatids, which branched very early from other eukaryotes, unique protein components constitute the lamina, kinetochores, and parts of the NPCs. Here we describe a phenylalanine-glycine (FG)-repeat nucleoporin, TbNup53b, that has dual localizations within the nucleoplasm and NPC. In addition to association with nucleoporins, TbNup53b interacts with a known trans-splicing component, TSR1, and has a role in controlling expression of surface proteins including the nucleolar periphery-located, procyclin genes. Significantly, while several nucleoporins are implicated in intranuclear transcriptional regulation in metazoa, TbNup53b appears orthologous to components of the yeast/human Nup49/Nup58 complex, for which no transcriptional functions are known. These data suggest that FG-Nups are frequently co-opted to transcriptional functions during evolution and extend the presence of FG-repeat nucleoporin control of gene expression to trypanosomes, suggesting that this is a widespread and ancient eukaryotic feature, as well as underscoring once more flexibility within nucleoporin function.

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.

Knock-Down of a Novel snoRNA in Tetrahymena Reveals a Dual Role in 5.8S rRNA Processing and Generation of a 26S rRNA Fragment

In eukaryotes, 18S, 5.8S, and 28S rRNAs are transcribed as precursor molecules that undergo extensive modification and nucleolytic processing to form the mature rRNA species. Central in the process are the small nucleolar RNAs (snoRNAs). The majority of snoRNAs guide site specific chemical modifications but a few are involved in defining pre-rRNA cleavages. Here, we describe an unusual snoRNA (TtnuCD32) belonging to the box C/D subgroup from the ciliate Tetrahymena thermophila. We show that TtnuCD32 is unlikely to function as a modification guide snoRNA and that it is critical for cell viability. Cell lines with genetic knock-down of TtnuCD32 were impaired in growth and displayed two novel and apparently unrelated phenotypes. The most prominent phenotype is the accumulation of processing intermediates of 5.8S rRNA. The second phenotype is the decrease in abundance of a ~100 nt 26S rRNA fragment of unknown function. Sequence analysis demonstrated that TtnuCD32 share features with the essential snoRNA U14 but an alternative candidate (TtnuCD25) was more closely related to other U14 sequences. This, together with the fact that the observed rRNA processing phenotypes were not similar to what has been observed in U14 depleted cells, suggests that TtnuCD32 is a U14 homolog that has gained novel functions.

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.

Natural and artificial small RNAs: a promising avenue of nucleic acid therapeutics for cancer

Since the failure of traditional therapy, gene therapy using functional DNA sequence and small RNA/DNA molecules (oligonucleotide) has become a promising avenue for cancer treatment. The discovery of RNA molecules has impelled researchers to investigate small regulatory RNA from various natural and artificial sources and determine a cogent target for controlling tumor progression. Small regulatory RNAs are used for therapeutic silencing of oncogenes and aberrant DNA repair response genes. Despite their advantages, therapies based on small RNAs exhibit limitations in terms of stability of therapeutic drugs, precision-based delivery in tissues, precision-based intercellular and intracellular targeting, and tumor heterogeneity-based responses. In this study, we summarize the potential and drawbacks of small RNAs in nucleic acid therapeutics for cancer.

A pseudouridylation switch in rRNA is implicated in ribosome function during the life cycle of Trypanosoma brucei

The protozoan parasite Trypanosoma brucei, which causes devastating diseases in humans and animals in sub-Saharan Africa, undergoes a complex life cycle between the mammalian host and the blood-feeding tsetse fly vector. However, little is known about how the parasite performs most molecular functions in such different environments. Here, we provide evidence for the intriguing possibility that pseudouridylation of rRNA plays an important role in the capacity of the parasite to transit between the insect midgut and the mammalian bloodstream. Briefly, we mapped pseudouridines (Ψ) on rRNA by Ψ-seq in procyclic form (PCF) and bloodstream form (BSF) trypanosomes. We detected 68 Ψs on rRNA, which are guided by H/ACA small nucleolar RNAs (snoRNA). The small RNome of both life cycle stages was determined by HiSeq and 83 H/ACAs were identified. We observed an elevation of 21 Ψs modifications in BSF as a result of increased levels of the guiding snoRNAs. Overexpression of snoRNAs guiding modification on H69 provided a slight growth advantage to PCF parasites at 30 °C. Interestingly, these modifications are predicted to significantly alter the secondary structure of the large subunit (LSU) rRNA suggesting that hypermodified positions may contribute to the adaption of ribosome function during cycling between the two hosts.

A tRNA methyltransferase paralog is important for ribosome stability and cell division in Trypanosoma brucei

Most eukaryotic ribosomes contain 26/28S, 5S, and 5.8S large subunit ribosomal RNAs (LSU rRNAs) in addition to the 18S rRNA of the small subunit (SSU rRNA). However, in kinetoplastids, a group of organisms that include medically important members of the genus Trypanosoma and Leishmania, the 26/28S large subunit ribosomal RNA is uniquely composed of 6 rRNA fragments. In addition, recent studies have shown the presence of expansion segments in the large ribosomal subunit (60S) of Trypanosoma brucei. Given these differences in structure, processing and assembly, T. brucei ribosomes may require biogenesis factors not found in other organisms. Here, we show that one of two putative 3-methylcytidine methyltransferases, TbMTase37 (a homolog of human methyltransferase-like 6, METTL6), is important for ribosome stability in T. brucei. TbMTase37 localizes to the nucleolus and depletion of the protein results in accumulation of ribosomal particles lacking srRNA 4 and reduced levels of polysome associated ribosomes. We also find that TbMTase37 plays a role in cytokinesis, as loss of the protein leads to multi-flagellated and multi-nucleated cells.

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.

Two splicing factors carrying serine-arginine motifs, TSR1 and TSR1IP, regulate splicing, mRNA stability, and rRNA processing in Trypanosoma brucei

In trypanosomes, mRNAs are processed by trans-splicing; in this process, a common exon, the spliced leader, is added to all mRNAs from a small RNA donor, the spliced leader RNA (SL RNA). However, little is known regarding how this process is regulated. In this study we investigated the function of two serine-arginine-rich proteins, TSR1 and TSR1IP, implicated in trans-splicing in Trypanosoma brucei. Depletion of these factors by RNAi suggested their role in both cis- and trans-splicing. Microarray was used to examine the transcriptome of the silenced cells. The level of hundreds of mRNAs was changed, suggesting that these proteins have a role in regulating only a subset of T. brucei mRNAs. Mass-spectrometry analyses of complexes associated with these proteins suggest that these factors function in mRNA stability, translation, and rRNA processing. We further demonstrate changes in the stability of mRNA as a result of depletion of the two TSR proteins. In addition, rRNA defects were observed under the depletion of U2AF35, TSR1, and TSR1IP, but not SF1, suggesting involvement of SR proteins in rRNA processing.

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.

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.

Early evolution of eukaryote feeding modes, cell structural diversity, and classification of the protozoan phyla Loukozoa, Sulcozoa, and Choanozoa

I discuss how different feeding modes and related cellular structures map onto the eukaryote evolutionary tree. Centrally important for understanding eukaryotic cell diversity are Loukozoa: ancestrally biciliate phagotrophic protozoa possessing a posterior cilium and ventral feeding groove into which ciliary currents direct prey. I revise their classification by including all anaerobic Metamonada as a subphylum and adding Tsukubamonas. Loukozoa, often with ciliary vanes, are probably ancestral to all protozoan phyla except Euglenozoa and Percolozoa and indirectly to kingdoms Animalia, Fungi, Plantae, and Chromista. I make a new protozoan phylum Sulcozoa comprising subphyla Apusozoa (Apusomonadida, Breviatea) and Varisulca (Diphyllatea; Planomonadida, Discocelida, Mantamonadida; Rigifilida). Understanding sulcozoan evolution clarifies the origins from them of opisthokonts (animals, fungi, Choanozoa) and Amoebozoa, and their evolutionary novelties; Sulcozoa and their descendants (collectively called podiates) arguably arose from Loukozoa by evolving posterior ciliary gliding and pseudopodia in their ventral groove. I explain subsequent independent cytoskeletal modifications, accompanying further shifts in feeding mode, that generated Amoebozoa, Choanozoa, and fungi. I revise classifications of Choanozoa, Conosa (Amoebozoa), and basal fungal phylum Archemycota. I use Choanozoa, Sulcozoa, Loukozoa, and Archemycota to emphasize the need for simply classifying ancestral (paraphyletic) groups and illustrate advantages of this for understanding step-wise phylogenetic advances.

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.

Paromomycin affects translation and vesicle-mediated trafficking as revealed by proteomics of paromomycin -susceptible -resistant Leishmania donovani

Leishmania donovani is a protozoan parasite that causes visceral leishmaniasis (VL) and is responsible for significant mortality and morbidity. Increasing resistance towards antimonial drugs poses a great challenge in chemotherapy of VL. Paromomycin is an aminoglycosidic antibiotic and is one of the drugs currently being used in the chemotherapy of cutaneous and visceral leishmaniasis. To understand the mode of action of this antibiotic at the molecular level, we have investigated the global proteome differences between the wild type AG83 strain and a paromomycin resistant (PRr) strain of L. donovani. Stable isotope labeling of amino acids in cell culture (SILAC) followed by quantitative mass spectrometry of the wild type AG83 strain and the paromomycin resistant (PRr) strain identified a total of 226 proteins at ≥95% confidence. Data analysis revealed upregulation of 29 proteins and down-regulation of 21 proteins in the PRr strain. Comparative proteomic analysis of the wild type and the paromomycin resistant strains showed upregulation of the ribosomal proteins in the resistant strain indicating role in translation. Elevated levels of glycolytic enzymes and stress proteins were also observed in the PRr strain. Most importantly, we observed upregulation of proteins that may have a role in intracellular survival and vesicular trafficking in the PRr strain. Furthermore, ultra-structural analysis by electron microscopy demonstrated increased number of vesicular vacuoles in PRr strain when compared to the wild-type strain. Drug affinity pull-down assay followed by mass spectrometery identified proteins in L. donovani wild type strain that were specifically and covalently bound to paromomycin. These results provide the first comprehensive insight into the mode of action and underlying mechanism of resistance to paromomycin in Leishmania donovani.

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.

Complete modification maps for the cytosolic small and large subunit rRNAs of Euglena gracilis: functional and evolutionary implications of contrasting patterns between the two rRNA components

In the protist Euglena gracilis, the cytosolic small subunit (SSU) rRNA is a single, covalently continuous species typical of most eukaryotes; in contrast, the large subunit (LSU) rRNA is naturally fragmented, comprising 14 separate RNA molecules instead of the bipartite (28S+5.8S) eukaryotic LSU rRNA typically seen. We present extensively revised secondary structure models of the E. gracilis SSU and LSU rRNAs and have mapped the positions of all of the modified nucleosides in these rRNAs (88 in SSU rRNA and 262 in LSU rRNA, with only 3 LSU rRNA modifications incompletely characterized). The relative proportions of ribose-methylated nucleosides and pseudouridine (∼60% and ∼35%, respectively) are closely similar in the two rRNAs; however, whereas the Euglena SSU rRNA has about the same absolute number of modifications as its human counterpart, the Euglena LSU rRNA has twice as many modifications as the corresponding human LSU rRNA. The increased levels of rRNA fragmentation and modification in E. gracilis LSU rRNA are correlated with a 3-fold increase in the level of mispairing in helical regions compared to the human LSU rRNA. In contrast, no comparable increase in mispairing is seen in helical regions of the SSU rRNA compared to its homologs in other eukaryotes. In view of the reported effects of both ribose-methylated nucleoside and pseudouridine residues on RNA structure, these correlations lead us to suggest that increased modification in the LSU rRNA may play a role in stabilizing a 'looser' structure promoted by elevated helical mispairing and a high degree of fragmentation.

Nucleolar localization of RNA binding proteins induced by actinomycin D and heat shock in Trypanosoma cruzi

In this work we show that under Actinomycin D (ActD) treatment, several RNA Binding Proteins (RBPs) involved in mRNA metabolism are relocalized into the nucleolus in Trypanosoma cruzi as a specific stress response. ATP depletion as well as kinase inhibition markedly reduced the nucleolar localization response, suggesting that an energy-dependent transport modulated by the phosphorylation status of the parasite might be required. Deletion analyses in one of such proteins, TcSR62, showed that a domain bearing basic amino acids located in the COOH terminal region was sufficient to promote its nucleolar relocalization. Interestingly, we showed that in addition to RBPs, poly(A)+ RNA is also accumulated into the nucleolus in response to ActD treatment. Finally, we found out that nucleolar relocalization of RBPs is also triggered by severe heat shock in a reversible way. Together, these results suggest that the nucleolus of an early divergent eukaryote is either able to sequester key factors related to mRNA metabolism in response to transcriptional stress or behaves as a RBP processing center, arguing in favour to the hypothesis that the non-traditional features of the nucleolus could be acquired early during evolution.

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.

3' processing in protists

Molecular biologists have traditionally focused on the very small corner of eukaryotic evolution that includes yeast and animals; even plants have been neglected. In this article, we describe the scant information that is available concerning RNA processing in the other four major eukaryotic groups, especially pathogenic protists. We focus mainly on polyadenylation and nuclear processing of stable RNAs. These processes have--where examined--been shown to be conserved, but there are many novel details. We also briefly mention other processing reactions such as splicing.

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.

Analysis of spliceosomal proteins in Trypanosomatids reveals novel functions in mRNA processing

In trypanosomatids, all mRNAs are processed via trans-splicing, although cis-splicing also occurs. In trans-splicing, a common small exon, the spliced leader (SL), which is derived from a small SL RNA species, is added to all mRNAs. Sm and Lsm proteins are core proteins that bind to U snRNAs and are essential for both these splicing processes. In this study, SmD3- and Lsm3-associated complexes were purified to homogeneity from Leishmania tarentolae. The purified complexes were analyzed by mass spectrometry, and 54 and 39 proteins were purified from SmD3 and Lsm complexes, respectively. Interestingly, among the proteins purified from Lsm3, no mRNA degradation factors were detected, as in Lsm complexes from other eukaryotes. The U1A complex was purified and mass spectrometry analysis identified, in addition to U1 small nuclear ribonucleoprotein (snRNP) proteins, additional co-purified proteins, including the polyadenylation factor CPSF73. Defects observed in cells silenced for U1 snRNP proteins suggest that the U1 snRNP functions exclusively in cis-splicing, although U1A also participates in polyadenylation and affects trans-splicing. The study characterized several trypanosome-specific nuclear factors involved in snRNP biogenesis, whose function was elucidated in Trypanosoma brucei. Conserved factors, such as PRP19, which functions at the heart of every cis-spliceosome, also affect SL RNA modification; GEMIN2, a protein associated with SMN (survival of motor neurons) and implicated in selective association of U snRNA with core Sm proteins in trypanosomes, is a master regulator of snRNP assembly. This study demonstrates the existence of trypanosomatid-specific splicing factors but also that conserved snRNP proteins possess trypanosome-specific functions.

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).

Analysis of C/D box snoRNA genes in vertebrates: The number of copies decreases in placental mammals

C/D box small nucleolar RNAs (snoRNAs) guide site-specific 2'-O-methylation of RNAs. Nearly all C/D box snoRNAs with known targets are involved in rRNA modification. In vertebrates, snoRNAs are encoded in introns of various genes and their processing is coupled with splicing of host gene pre-mRNA. Here, the genes encoding C/D box snoRNAs that guide 2'-O-methylation of rRNA were identified and analyzed in vertebrate genomes. The number of copies of most C/D box snoRNA genes proved to be lower in placental mammals compared to other vertebrates. This can be due to smaller oocytes and accordingly lower number of ribosomes in them in eutherians. The targets of snoRNAs encoded by single-copy and multiple-copy genes proved to have different distribution in rRNAs. The causes of this difference are discussed. In some cases, the transcripts of homologous C/D box RNA genes were shown to guide the modification of neighboring nucleotides in rRNA. C/D box snoRNA pseudogenes were found in all vertebrate classes. Three novel C/D box snoRNAs were found in Xenopus tropicalis that may guide 2'-O-methylation of Xenopus-specific rRNA sites. A list of 922 annotated C/D box snoRNA genes is presented.

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.