Small nucleolar RNA clusters in trypanosomatid Leptomonas collosoma. Genome organization, expression studies, and the potential role of sequences present upstream from the first repeated cluster

TIGR-Tas: A family of modular RNA-guided DNA-targeting systems in prokaryotes and their viruses

RNA-guided systems provide remarkable versatility, enabling diverse biological functions. Through iterative structural and sequence homology-based mining starting with a guide RNA-interaction domain of Cas9, we identified a family of RNA-guided DNA-targeting proteins in phage and parasitic bacteria. Each system consists of a tandem interspaced guide RNA (TIGR) array and a TIGR-associated (Tas) protein containing a nucleolar protein (Nop) domain, sometimes fused to HNH (TasH)- or RuvC (TasR)-nuclease domains. We show that TIGR arrays are processed into 36-nucleotide RNAs (tigRNAs) that direct sequence-specific DNA binding through a tandem-spacer targeting mechanism. TasR can be reprogrammed for precise DNA cleavage, including in human cells. The structure of TasR reveals striking similarities to box C/D small nucleolar ribonucleoproteins and IS110 RNA-guided transposases, providing insights into the evolution of diverse RNA-guided systems.

H/ACA Small Ribonucleoproteins: Structural and Functional Comparison Between Archaea and Eukaryotes

During ribosome synthesis, ribosomal RNA is modified through the formation of many pseudouridines and methylations which contribute to ribosome function across all domains of life. In archaea and eukaryotes, pseudouridylation of rRNA is catalyzed by H/ACA small ribonucleoproteins (sRNPs) utilizing different H/ACA guide RNAs to identify target uridines for modification. H/ACA sRNPs are conserved in archaea and eukaryotes, as they share a common general architecture and function, but there are also several notable differences between archaeal and eukaryotic H/ACA sRNPs. Due to the higher protein stability in archaea, we have more information on the structure of archaeal H/ACA sRNPs compared to eukaryotic counterparts. However, based on the long history of yeast genetic and other cellular studies, the biological role of H/ACA sRNPs during ribosome biogenesis is better understood in eukaryotes than archaea. Therefore, this review provides an overview of the current knowledge on H/ACA sRNPs from archaea, in particular their structure and function, and relates it to our understanding of the roles of eukaryotic H/ACA sRNP during eukaryotic ribosome synthesis and beyond. Based on this comparison of our current insights into archaeal and eukaryotic H/ACA sRNPs, we discuss what role archaeal H/ACA sRNPs may play in the formation of ribosomes.

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.

Insights on a putative aminoacyl-tRNA-protein transferase of Leishmania major

The N-end rule pathway leads to regulated proteolysis as an adaptive response to external stress and is ubiquitous from bacteria to mammals. In this study, we investigated a gene coding for a putative core enzyme of this post-translational regulatory pathway in Leishmania major, which may be crucial during cytodifferentiation and the environment adaptive responses of the parasite. Leucyl, phenylalanyl-tRNA protein transferase and arginyl-tRNA protein transferase are key components of this pathway in E. coli and eukaryotes, respectively. They catalyze the specific conjugation of leucine, phenylalanine or arginine to proteins containing exposed N-terminal amino acid residues, which are recognized by the machinery for the targeted proteolysis. Here, we characterized a conserved hypothetical protein coded by the LmjF.21.0725 gene in L. major. In silico analysis suggests that the LmjF.21.0725 protein is highly conserved among species of Leishmania and might belong to the Acyl CoA-N-acyltransferases (NAT) superfamily of proteins. Immunofluorescence cell imaging indicates that the cytosolic localization of the studied protein and the endogenous levels of the protein in promastigotes are barely detectable by western blotting assay. The knockout of the two alleles of LmjF.21.0725 by homologous recombination was only possible in the heterozygous transfectant expressing LmjF.21.0725 as a transgene from a plasmid. Moreover, the kinetics of loss of the plasmid in the absence of drug pressure suggests that maintenance of the gene is essential for promastigote survival. Here, evidence is provided that this putative aminoacyl tRNA-protein transferase is essential for parasite survival. The enzyme activity and corresponding post-translational regulatory pathway are yet to be investigated.

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.

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.

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

Background

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

Results

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

Conclusion

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

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

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

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

Background

Detection of non coding RNA (ncRNA) molecules is a major bioinformatics challenge. This challenge is particularly difficult when attempting to detect H/ACA molecules which are involved in converting uridine to pseudouridine on rRNA in trypanosomes, because these organisms have unique H/ACA molecules (termed H/ACA-like) that lack several of the features that characterize H/ACA molecules in most other organisms.

Results

We present here a computational tool called Psiscan, which was designed to detect H/ACA-like molecules in trypanosomes. We started by analyzing known H/ACA-like molecules and characterized their crucial elements both computationally and experimentally.

Next, we set up constraints based on this analysis and additional phylogenic and functional data to rapidly scan three trypanosome genomes (T. brucei, T. cruzi and L. major) for sequences that observe these constraints and are conserved among the species. In the next step, we used minimal energy calculation to select the molecules that are predicted to fold into a lowest energy structure that is consistent with the constraints. In the final computational step, we used a Support Vector Machine that was trained on known H/ACA-like molecules as positive examples and on negative examples of molecules that were identified by the computational analyses but were shown experimentally not to be H/ACA-like molecules. The leading candidate molecules predicted by the SVM model were then subjected to experimental validation.

Conclusion

The experimental validation showed 11 molecules to be expressed (4 out of 25 in the intermediate stage and 7 out of 19 in the final validation after the machine learning stage). Five of these 11 molecules were further shown to be bona fide H/ACA-like molecules. As snoRNA in trypanosomes are organized in clusters, the new H/ACA-like molecules could be used as starting points to manually search for additional molecules in their neighbourhood. All together this study increased our repertoire by fourteen H/ACA-like and six C/D snoRNAs molecules from T. brucei and L. Major. In addition the experimental analysis revealed that six ncRNA molecules that are expressed are not downregulated in CBF5 silenced cells, suggesting that they have structural features of H/ACA-like molecules but do not have their standard function. We termed this novel class of molecules AGA-like, and we are exploring their function.

This study demonstrates the power of tight collaboration between computational and experimental approaches in a combined effort to reveal the repertoire of ncRNA molecles.

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

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

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

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

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

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

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

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

Pseudouridine-guide RNAs and other Cbf5p-associated RNAs in Euglena gracilis

In eukaryotes, box H/ACA small nucleolar RNAs (snoRNAs) guide sites of pseudouridine (Psi) formation in rRNA. These snoRNAs reside in RNP complexes containing the putative Psi synthase, Cbf5p. In this study we have identified Cbf5p-associated RNAs in Euglena gracilis, an early diverging eukaryote, by immunoprecipitating Cbf5p-containing complexes from cellular extracts. We characterized one box H/ACA-like RNA which, however, does not appear to guide Psi formation in rRNA. We also identified four single Psi-guide box AGA RNAs. We determined target sites for these putative Psi-guide RNAs and confirmed that the predicted Psi modifications do, in fact, occur at these positions in Euglena rRNA. The Cbf5p-associated snoRNAs appear to be encoded by multicopy genes, some of which are clustered in the genome together with methylation-guide snoRNA genes. These modification-guide snoRNAs and snoRNA genes are the first ones to be reported in euglenid protists, the evolutionary sister group to the kinetoplastid protozoa. Unexpectedly, we also found and have partially characterized a selenocysteine tRNA homolog in the anti-Cbf5p-immunoprecipitated sample.