Systematic study of sequence motifs for RNA trans splicing in Trypanosoma brucei

Post-transcriptional reprogramming by thousands of mRNA untranslated regions in trypanosomes

Although genome-wide polycistronic transcription places major emphasis on post-transcriptional controls in trypanosomatids, messenger RNA cis-regulatory untranslated regions (UTRs) have remained largely uncharacterised. Here, we describe a genome-scale massive parallel reporter assay coupled with 3'-UTR-seq profiling in the African trypanosome and identify thousands of regulatory UTRs. Increased translation efficiency was associated with dosage of adenine-rich poly-purine tracts (pPuTs). An independent assessment of native UTRs using machine learning based predictions confirmed the robust correspondence between pPuTs and positive control, as did an assessment of synthetic UTRs. Those 3'-UTRs associated with upregulated expression in bloodstream-stage cells were also enriched in uracil-rich poly-pyrimidine tracts, suggesting a mechanism for developmental activation through pPuT 'unmasking'. Thus, we describe a cis-regulatory UTR sequence 'code' that underpins gene expression control in the context of a constitutively transcribed genome. We conclude that thousands of UTRs post-transcriptionally reprogram gene expression profiles in trypanosomes.

Unique features of conventional and nonconventional introns in Euglena gracilis

BACKGROUND

Nuclear introns in Euglenida have been understudied. This study aimed to investigate nuclear introns in Euglenida by identifying a large number of introns in Euglena gracilis (E. gracilis), including cis-spliced conventional and nonconventional introns, as well as trans-spliced outrons. We also examined the sequence characteristics of these introns.

RESULTS

A total of 28,337 introns and 11,921 outrons were identified. Conventional and nonconventional introns have distinct splice site features; the former harbour canonical GT/C-AG splice sites, whereas the latter are capable of forming structured motifs with their terminal sequences. We observed that short introns had a preference for canonical GT-AG introns. Notably, conventional introns and outrons in E. gracilis exhibited a distinct cytidine-rich polypyrimidine tract, in contrast to the thymidine-rich tracts observed in other organisms. Furthermore, the SL-RNAs in E. gracilis, as well as in other trans-splicing species, can form a recently discovered motif called the extended U6/5' ss duplex with the respective U6s. We also describe a novel type of alternative splicing pattern in E. gracilis. The tandem repeat sequences of introns in this protist were determined, and their contents were comparable to those in humans.

CONCLUSIONS

Our findings highlight the unique features of E. gracilis introns and provide insights into the splicing mechanism of these introns, as well as the genomics and evolution of Euglenida.

Sequences and proteins that influence mRNA processing in Trypanosoma brucei: Evolutionary conservation of SR-domain and PTB protein functions

BACKGROUND

Spliced leader trans splicing is the addition of a short, capped sequence to the 5' end of mRNAs. It is widespread in eukaryotic evolution, but factors that influence trans splicing acceptor site choice have been little investigated. In Kinetoplastids, all protein-coding mRNAs are 5' trans spliced. A polypyrimidine tract is usually found upstream of the AG splice acceptor, but there is no branch point consensus; moreover, splicing dictates polyadenylation of the preceding mRNA, which is a validated drug target.

METHODOLOGY AND PRINCIPAL FINDINGS

We here describe a trans splicing reporter system that can be used for studies and screens concerning the roles of sequences and proteins in processing site choice and efficiency. Splicing was poor with poly(U) tracts less than 9 nt long, and was influenced by an intergenic region secondary structure. A screen for signals resulted in selection of sequences that were on average 45% U and 35% C. Tethering of either the splicing factor SF1, or the cleavage and polyadenylation factor CPSF3 within the intron stimulated processing in the correct positions, while tethering of two possible homologues of Opisthokont PTB inhibited processing. In contrast, tethering of SR-domain proteins RBSR1, RBSR2, or TSR1 or its interaction partner TSR1IP, promoted use of alternative signals upstream of the tethering sites. RBSR1 interacts predominantly with proteins implicated in splicing, whereas the interactome of RBSR2 is more diverse.

CONCLUSIONS

Our selectable constructs are suitable for screens of both sequences, and proteins that affect mRNA processing in T. brucei. Our results suggest that the functions of PTB and SR-domain proteins in splice site definition may already have been present in the last eukaryotic common ancestor.

The Trypanosoma brucei RNA-binding protein DRBD18 ensures correct mRNA trans splicing and polyadenylation patterns

The parasite Trypanosoma brucei grows as bloodstream forms in mammals, and as procyclic forms in tsetse flies. Transcription is polycistronic, all mRNAs are trans spliced, and polyadenylation sites are defined by downstream splicing signals. Expression regulation therefore depends heavily on post-transcriptional mechanisms. The RNA-binding protein DRBD18 was previously implicated in the export of some mRNAs from the nucleus in procyclic forms. It copurifies the outer ring of the nuclear pore, mRNA export factors and exon-junction-complex proteins. We show that for more than 200 mRNAs, DRBD18 depletion caused preferential accumulation of versions with shortened 3'-untranslated regions, arising from use of polyadenylation sites that were either undetectable or rarely seen in nondepleted cells. The shortened mRNAs were often, but not always, more abundant in depleted cells than the corresponding longer versions in normal cells. Their appearance was linked to the appearance of trans-spliced, polyadenylated RNAs containing only downstream 3'-untranslated region-derived sequences. Experiments with one mRNA suggested that nuclear retention alone, through depletion of MEX67, did not affect mRNA length, suggesting a specific effect of DRBD18 on processing. DRBD18-bound mRNAs were enriched in polypyrimidine tract motifs, and DRBD18 was found in both the nucleus and the cytoplasm. We therefore suggest that in the nucleus, DRBD18 might bind to polypyrimidine tracts in 3'-UTRs of mRNA precursors. Such binding might both prevent recognition of mRNA-internal polypyrimidine tracts by splicing factors, and promote export of the processed bound mRNAs to the cytosol.

A Trypanosoma brucei ORFeome-Based Gain-of-Function Library Identifies Genes That Promote Survival during Melarsoprol Treatment

Trypanosoma brucei is an early branching protozoan parasite that causes human and animal African trypanosomiasis. Forward genetics approaches are powerful tools for uncovering novel aspects of trypanosomatid biology, pathogenesis, and therapeutic approaches against trypanosomiasis. Here, we have generated a T. brucei cloned ORFeome consisting of >90% of the targeted 7,245 genes and used it to make an inducible gain-of-function parasite library broadly applicable to large-scale forward genetic screens. We conducted a proof-of-principle genetic screen to identify genes whose expression promotes survival in melarsoprol, a critical drug of last resort. The 57 genes identified as overrepresented in melarsoprol survivor populations included the gene encoding the rate-limiting enzyme for the biosynthesis of an established drug target (trypanothione), validating the tool. In addition, novel genes associated with gene expression, flagellum localization, and mitochondrion localization were identified, and a subset of those genes increased melarsoprol resistance upon overexpression in culture. These findings offer new insights into trypanosomatid basic biology, implications for drug targets, and direct or indirect drug resistance mechanisms. This study generated a T. brucei ORFeome and gain-of-function parasite library, demonstrated the library's usefulness in forward genetic screening, and identified novel aspects of melarsoprol resistance that will be the subject of future investigations. These powerful genetic tools can be used to broadly advance trypanosomatid research.IMPORTANCE Trypanosomatid parasites threaten the health of more than 1 billion people worldwide. Because their genomes are highly diverged from those of well-established eukaryotes, conservation is not always useful in assigning gene functions. However, it is precisely among the trypanosomatid-specific genes that ideal therapeutic targets might be found. Forward genetics approaches are an effective way to identify novel gene functions. We used an ORFeome approach to clone a large percentage of Trypanosoma brucei genes and generate a gain-of-function parasite library. This library was used in a genetic screen to identify genes that promote resistance to the clinically significant yet highly toxic drug melarsoprol. Hits arising from the screen demonstrated the library's usefulness in identifying known pathways and uncovered novel aspects of resistance mediated by proteins localized to the flagellum and mitochondrion. The powerful new genetic tools generated herein are expected to promote advances in trypanosomatid biology and therapeutic development in the years to come.

Proteomic analysis of Trypanosoma cruzi spliceosome complex

The unicellular protists of the group Kinetoplastida include the genera Leishmania and Trypanosoma, which are pathogens of invertebrate and vertebrate animals. Despite their medical and economical importance, critical aspects of their biology such as specific molecular characteristics of gene expression regulation are just beginning to be deciphered. Gene expression regulation also depends on post-transcriptional processing steps, such as the trans-splicing process. Despite being widely used in trypanosomes, trans-splicing is a rare event in other eukaryotes. We sought to describe the protein composition of spliceosomes in epimastigotes of T. cruzi, the etiological agent of Chagas disease. We used two TAP-tagged proteins to affinity purify spliceosomes and analyzed their composition by mass spectrometry. Among the 115 identified proteins we detected conserved spliceosome components, as Sm and LSm proteins, RNA helicases, U2- and U5-snRNP specific proteins. Importantly, by comparing our data with proteomic data of human and T. brucei spliceosome complexes, we observed a core group of proteins common to all spliceosomes. By using amino acid sequence comparisons, we identified RNA-associated proteins that might be involved with splicing regulation in T. cruzi, namely the orthologous of WDR33, PABPCL1 and three different HNRNPs. Data are available via ProteomeXchange with identifier PXD018776.

Upstream ORFs Influence Translation Efficiency in the Parasite Trypanosoma cruzi

It is generally accepted that the presence of ORFs in the 5' untranslated region of eukaryotic transcripts modulates the production of proteins by controlling the translation initiation rate of the main CDS. In trypanosomatid parasites, which almost exclusively depend on post-transcriptional mechanisms to regulate gene expression, translation has been identified as a key step. However, the mechanisms of control of translation are not fully understood. In the present work, we have annotated the 5'UTRs of the Trypanosoma cruzi genome both in epimastigotes and metacyclic trypomastigotes and, using a stringent classification approach, we identified putative regulatory uORFs in about 9% of the analyzed 5'UTRs. The translation efficiency (TE) and translational levels of transcripts containing putative repressive uORFs were found to be significantly reduced. These findings are supported by the fact that proteomic methods only identify a low number of proteins coded by transcripts containing repressive uORF. We additionally show that AUG is the main translation initiator codon of repressive uORFs in T. cruzi. Interestingly, the decrease in TE is more pronounced when the uORFs overlaps the main CDS. In conclusion, we show that the presence of the uORF and features such as initiation codon and/or location of the uORFs may be acting to fine tune translation levels in these parasites.

Regulation of gene expression in trypanosomatids: living with polycistronic transcription

In trypanosomes, RNA polymerase II transcription is polycistronic and individual mRNAs are excised by trans-splicing and polyadenylation. The lack of individual gene transcription control is compensated by control of mRNA processing, translation and degradation. Although the basic mechanisms of mRNA decay and translation are evolutionarily conserved, there are also unique aspects, such as the existence of six cap-binding translation initiation factor homologues, a novel decapping enzyme and an mRNA stabilizing complex that is recruited by RNA-binding proteins. High-throughput analyses have identified nearly a hundred regulatory mRNA-binding proteins, making trypanosomes valuable as a model system to investigate post-transcriptional regulation.

UTRme: A Scoring-Based Tool to Annotate Untranslated Regions in Trypanosomatid Genomes

Most signals involved in post-transcriptional regulatory networks are located in the untranslated regions (UTRs) of the mRNAs. Therefore, to deepen our understanding of gene expression regulation, delimitation of these regions with high accuracy is needed. The trypanosomatid lineage includes a variety of parasitic protozoans causing a significant worldwide burden on human health. Given their peculiar mechanisms of gene expression, these organisms depend on post-transcriptional regulation as the main level of gene expression control. In this context, the definition of the UTR regions becomes of key importance. We have developed UTR-mini-exon (UTRme), a graphical user interface (GUI) stand-alone application to identify and annotate 5′ and 3′ UTR regions in a highly accurate way. UTRme implements a multiple scoring system tailored to address the issue of false positive UTR assignment that frequently arise because of the characteristics of the intergenic regions. Even though it was developed for trypanosomatids, the tool can be used to predict 3′ sites in any eukaryote and 5′ UTRs in any organism where trans-splicing occurs (such as the model organism C. elegans). UTRme offers a way for non-bioinformaticians to precisely determine UTRs from transcriptomic data. The tool is freely available via the conda and github repositories.

Exploiting CRISPR-Cas9 technology to investigate individual histone modifications

Despite their importance for most DNA-templated processes, the function of individual histone modifications has remained largely unknown because in vivo mutational analyses are lacking. The reason for this is that histone genes are encoded by multigene families and that tools to simultaneously edit multiple genomic loci with high efficiency are only now becoming available. To overcome these challenges, we have taken advantage of the power of CRISPR-Cas9 for precise genome editing and of the fact that most DNA repair in the protozoan parasite Trypanosoma brucei occurs via homologous recombination. By establishing an episome-based CRISPR-Cas9 system for T. brucei, we have edited wild type cells without inserting selectable markers, inserted a GFP tag between an ORF and its 3'UTR, deleted both alleles of a gene in a single transfection, and performed precise editing of genes that exist in multicopy arrays, replacing histone H4K4 with H4R4 in the absence of detectable off-target effects. The newly established genome editing toolbox allows for the generation of precise mutants without needing to change other regions of the genome, opening up opportunities to study the role of individual histone modifications, catalytic sites of enzymes or the regulatory potential of UTRs in their endogenous environments.

Life cycle adapted upstream open reading frames (uORFs) in Trypanosoma congolense: A post-transcriptional approach to accurate gene regulation

The presented work explores the regulatory influence of upstream open reading frames (uORFs) on gene expression in Trypanosoma congolense. More than 31,000 uORFs in total were identified and characterized here. We found evidence for the uORFs’ appearance in the transcriptome to be correlated with proteomic expression data, clearly indicating their repressive potential in T. congolense, which has to rely on post-transcriptional gene expression regulation due to its unique genomic organization. Our data show that uORF’s translation repressive potential does not only correlate with elemental sequence features such as length, position and quantity, but involves more subtle components, in particular the codon and amino acid profiles. This corresponds with the popular mechanistic model of a ribosome shedding initiation factors during the translation of a uORF, which can prevent reinitiation at the downstream start codon of the actual protein-coding sequence, due to the former extensive consumption of crucial translation components. We suggest that uORFs with uncommon codon and amino acid usage can slow down the translation elongation process in T. congolense, systematically deplete the limited factors, and restrict downstream reinitiation, setting up a bottleneck for subsequent translation of the protein-coding sequence. Additionally we conclude that uORFs dynamically influence the T. congolense life cycle. We found evidence that transition to epimastigote form could be supported by gain of uORFs due to alternative trans-splicing, which down-regulate housekeeping genes’ expression and render the trypanosome in a metabolically reduced state of endurance.

Codon choice directs constitutive mRNA levels in trypanosomes

Selective transcription of individual protein coding genes does not occur in trypanosomes and the cellular copy number of each mRNA must be determined post-transcriptionally. Here, we provide evidence that codon choice directs the levels of constitutively expressed mRNAs. First, a novel codon usage metric, the gene expression codon adaptation index (geCAI), was developed that maximised the relationship between codon choice and the measured abundance for a transcriptome. Second, geCAI predictions of mRNA levels were tested using differently coded GFP transgenes and were successful over a 25-fold range, similar to the variation in endogenous mRNAs. Third, translation was necessary for the accelerated mRNA turnover resulting from codon choice. Thus, in trypanosomes, the information determining the levels of most mRNAs resides in the open reading frame and translation is required to access this information.

Landscape of the spliced leader trans-splicing mechanism in Schistosoma mansoni

Spliced leader dependent trans-splicing (SLTS) has been described as an important RNA regulatory process that occurs in different organisms, including the trematode Schistosoma mansoni. We identified more than seven thousand putative SLTS sites in the parasite, comprising genes with a wide spectrum of functional classes, which underlines the SLTS as a ubiquitous mechanism in the parasite. Also, SLTS gene expression levels span several orders of magnitude, showing that SLTS frequency is not determined by the expression level of the target gene, but by the presence of particular gene features facilitating or hindering the trans-splicing mechanism. Our in-depth investigation of SLTS events demonstrates widespread alternative trans-splicing (ATS) acceptor sites occurring in different regions along the entire gene body, highlighting another important role of SLTS generating alternative RNA isoforms in the parasite, besides the polycistron resolution. Particularly for introns where SLTS directly competes for the same acceptor substrate with cis-splicing, we identified for the first time additional and important features that might determine the type of splicing. Our study substantially extends the current knowledge of RNA processing by SLTS in S. mansoni, and provide basis for future studies on the trans-splicing mechanism in other eukaryotes.

Variant surface glycoprotein density defines an immune evasion threshold for African trypanosomes undergoing antigenic variation

Trypanosoma brucei is a protozoan parasite that evades its host's adaptive immune response by repeatedly replacing its dense variant surface glycoprotein (VSG) coat from its large genomic VSG repertoire. While the mechanisms regulating VSG gene expression and diversification have been examined extensively, the dynamics of VSG coat replacement at the protein level, and the impact of this process on successful immune evasion, remain unclear. Here we evaluate the rate of VSG replacement at the trypanosome surface following a genetic VSG switch, and show that full coat replacement requires several days to complete. Using in vivo infection assays, we demonstrate that parasites undergoing coat replacement are only vulnerable to clearance via early IgM antibodies for a limited time. Finally, we show that IgM loses its ability to mediate trypanosome clearance at unexpectedly early stages of coat replacement based on a critical density threshold of its cognate VSGs on the parasite surface. Trypanosoma brucei evades the host immune system through replacement of a variant surface glycoprotein (VSG) coat. Here, the authors show that VSG replacement takes several days to complete, and the parasite is vulnerable to the host immune system for a short period of time during the process.

GT-rich promoters can drive RNA pol II transcription and deposition of H2A.Z in African trypanosomes

Genome‐wide transcription studies are revealing an increasing number of “dispersed promoters” that, unlike “focused promoters”, lack well‐conserved sequence motifs and tight regulation. Dispersed promoters are nevertheless marked by well‐defined chromatin structures, suggesting that specific sequence elements must exist in these unregulated promoters. Here, we have analyzed regions of transcription initiation in the eukaryotic parasite Trypanosoma brucei, in which RNA polymerase II transcription initiation occurs over broad regions without distinct promoter motifs and lacks regulation. Using a combination of site‐specific and genome‐wide assays, we identified GT‐rich promoters that can drive transcription and promote the targeted deposition of the histone variant H2A.Z in a genomic context‐dependent manner. In addition, upon mapping nucleosome occupancy at high resolution, we find nucleosome positioning to correlate with RNA pol II enrichment and gene expression, pointing to a role in RNA maturation. Nucleosome positioning may thus represent a previously unrecognized layer of gene regulation in trypanosomes. Our findings show that even highly dispersed, unregulated promoters contain specific DNA elements that are able to induce transcription and changes in chromatin structure.

Gene expression in Kinetoplastids

Kinetoplastid parasites adapt to different environments with wide-reaching control of gene expression, but transcription of nuclear protein-coding genes is polycistronic: there is no individual control of transcription initiation. Mature mRNAs are made by co-transcriptional trans splicing and polyadenylation, and competition between processing and nuclear degradation may contribute to regulation of mRNA levels. In the cytosol both the extent to which mRNAs are translated, and mRNA decay rates, vary enormously. I here highlight gaps in our knowledge: no measurements of transcription initiation or elongation rates; no measurements of how, precisely, mRNA processing and nuclear degradation control mRNA levels; and extremely limited understanding of the contributions of different translation initiation factors and RNA-binding proteins to mRNA fate.

Horizontal transfer - imperative mission of acellular life forms, Acytota

Acytota is a kingdom of life covering satellites, plasmids, transposable elements, viroids and viruses, all outside the conventional tree of life but satisfying most life definitions. This review focuses on some aspects of Acytota, their "genomes" and life styles, the dominance of transposable elements and their evolutionary influence on other life forms in order to vindicate the Acytota as a life kingdom no more polyphyletic than other kingdoms and its members no more parasitic than other life forms.

Integrative analysis of the Trypanosoma brucei gene expression cascade predicts differential regulation of mRNA processing and unusual control of ribosomal protein expression

Background

Trypanosoma brucei is a unicellular parasite which multiplies in mammals (bloodstream form) and Tsetse flies (procyclic form). Trypanosome RNA polymerase II transcription is polycistronic, individual mRNAs being excised by trans splicing and polyadenylation. We previously made detailed measurements of mRNA half-lives in bloodstream and procyclic forms, and developed a mathematical model of gene expression for bloodstream forms. At the whole transcriptome level, many bloodstream-form mRNAs were less abundant than was predicted by the model.

Results

We refined the published mathematical model and extended it to the procyclic form. We used the model, together with known mRNA half-lives, to predict the abundances of individual mRNAs, assuming rapid, unregulated mRNA processing; then we compared the results with measured mRNA abundances. Remarkably, the abundances of most mRNAs in procyclic forms are predicted quite well by the model, being largely explained by variations in mRNA decay rates and length. In bloodstream forms substantially more mRNAs are less abundant than predicted. We list mRNAs that are likely to show particularly slow or inefficient processing, either in both forms or with developmental regulation. We also measured ribosome occupancies of all mRNAs in trypanosomes grown in the same conditions as were used to measure mRNA turnover. In procyclic forms there was a weak positive correlation between ribosome density and mRNA half-life, suggesting cross-talk between translation and mRNA decay; ribosome density was related to the proportion of the mRNA on polysomes, indicating control of translation initiation. Ribosomal protein mRNAs in procyclics appeared to be exceptionally rapidly processed but poorly translated.

Conclusions

Levels of mRNAs in procyclic form trypanosomes are determined mainly by length and mRNA decay, with some control of precursor processing. In bloodstream forms variations in nuclear events play a larger role in transcriptome regulation, suggesting aquisition of new control mechanisms during adaptation to mammalian parasitism.

Electronic supplementary material

The online version of this article (doi:10.1186/s12864-016-2624-3) contains supplementary material, which is available to authorized users.

Integrating ribosomal promoter vectors that offer a choice of constitutive expression profiles in Leishmania donovani

We have designed a novel series of integrating ribosomal RNA promoter vectors with five incrementally different constitutive expression profiles, covering a 250-fold range. Differential expression was achieved by placing different combinations of synthetic or leishmanial DNA sequences upstream and downstream of the transgene coding sequence in order to modulate pre-mRNA processing efficiency and mRNA stability, respectively. All of the vectors have extensive multiple cloning sites, and versions are available for producing N- or C- terminal GFP fusions at each of the possible relative expression levels. In addition, the modular configuration of the vectors allows drug resistance cassettes and other components to be readily exchanged. In toto, these vectors should be useful additions to the toolkit available for molecular and genetic studies of Leishmania donovani.

Insights into the Regulation of mRNA Processing of Polycistronic Transcripts Mediated by DRBD4/PTB2, a Trypanosome Homolog of the Polypyrimidine Tract-Binding Protein

Trypanosomes regulate gene expression mostly by posttranscriptional mechanisms, including control of mRNA turnover and translation efficiency. This regulation is carried out via certain elements located at the 3'-untranslated regions of mRNAs, which are recognized by RNA-binding proteins. In trypanosomes, trans-splicing is of central importance to control mRNA maturation. We have previously shown that TcDRBD4/PTB2, a trypanosome homolog of the human polypyrimidine tract-binding protein splicing regulator, interacts with the intergenic region of one specific dicistronic transcript, referred to as TcUBP (and encoding for TcUBP1 and TcUBP2, two closely kinetoplastid-specific proteins). In this work, a survey of TcUBP RNA processing revealed certain TcDRBD4/PTB2-regulatory elements within its intercistronic region, which are likely to influence the trans-splicing rate of monocistronic-derived transcripts. Furthermore, TcDRBD4/PTB2 overexpression in epimastigote cells notably decreased both UBP1 and UBP2 protein expression. This type of posttranscriptional gene regulatory mechanism could be extended to other transcripts as well, as we identified several other RNA precursor molecules that specifically bind to TcDRBD4/PTB2. Altogether, these findings support a model in which TcDRBD4/PTB2-containing ribonucleoprotein complexes can prevent trans-splicing. This could represent another stage of gene expression regulation mediated by the masking of trans-splicing/polyadenylation signals.

Transcriptomic profiling of gene expression and RNA processing during Leishmania major differentiation

Protozoan parasites of the genus Leishmania are the etiological agents of leishmaniasis, a group of diseases with a worldwide incidence of 0.9–1.6 million cases per year. We used RNA-seq to conduct a high-resolution transcriptomic analysis of the global changes in gene expression and RNA processing events that occur as L. major transforms from non-infective procyclic promastigotes to infective metacyclic promastigotes. Careful statistical analysis across multiple biological replicates and the removal of batch effects provided a high quality framework for comprehensively analyzing differential gene expression and transcriptome remodeling in this pathogen as it acquires its infectivity. We also identified precise 5′ and 3′ UTR boundaries for a majority of Leishmania genes and detected widespread alternative trans-splicing and polyadenylation. An investigation of possible correlations between stage-specific preferential trans-splicing or polyadenylation sites and differentially expressed genes revealed a lack of systematic association, establishing that differences in expression levels cannot be attributed to stage-regulated alternative RNA processing. Our findings build on and improve existing expression datasets and provide a substantially more detailed view of L. major biology that will inform the field and potentially provide a stronger basis for drug discovery and vaccine development efforts.

The selenocysteine tRNA gene in leishmania major is transcribed by both RNA polymerase II and RNA polymerase III

Eukaryotic tRNAs, transcribed by RNA polymerase III (Pol III), contain boxes A and B as internal promoter elements. One exception is the selenocysteine (Sec) tRNA (tRNA-Sec), whose transcription is directed by an internal box B and three extragenic sequences in vertebrates. Here we report on the transcriptional analysis of the tRNA-Sec gene in the protozoan parasite Leishmania major. This organism has unusual mechanisms of gene expression, including Pol II polycistronic transcription and maturation of mRNAs by trans splicing, a process that attaches a 39-nucleotide miniexon to the 5' end of all the mRNAs. In L. major, tRNA-Sec is encoded by a single gene inserted into a Pol II polycistronic unit, in contrast to most tRNAs, which are clustered at the boundaries of polycistronic units. 5' rapid amplification of cDNA ends and reverse transcription-PCR experiments showed that some tRNA-Sec transcripts contain the miniexon at the 5' end and a poly(A) tail at the 3' end, indicating that the tRNA-Sec gene is polycistronically transcribed by Pol II and processed by trans splicing and polyadenylation, as was recently reported for the tRNA-Sec genes in the related parasite Trypanosoma brucei. However, nuclear run-on assays with RNA polymerase inhibitors and with cells that were previously UV irradiated showed that the tRNA-Sec gene in L. major is also transcribed by Pol III. Thus, our results indicate that RNA polymerase specificity in Leishmania is not absolute in vivo, as has recently been found in other eukaryotes.

Risk-associated coding synonymous SNPs in type 2 diabetes and neurodegenerative diseases: genetic silence and the underrated association with splicing regulation and epigenetics

Single nucleotide polymorphisms (SNPs) are tentatively critical with regard to disease predisposition, but coding synonymous SNPs (sSNPs) are generally considered "neutral". Nevertheless, sSNPs in serine/arginine-rich (SR) and splice-site (SS) exonic splicing enhancers (ESEs) or in exonic CpG methylation targets, could be decisive for splicing, particularly in aging-related conditions, where mis-splicing is frequently observed. We presently identified 33 genes T2D-related and 28 related to neurodegenerative diseases, by investigating the impact of the corresponding coding sSNPs on splicing and using gene ontology data and computational tools. Potentially critical (prominent) sSNPs comply with the following criteria: changing the splicing potential of prominent SR-ESEs or of significant SS-ESEs by >1.5 units (Δscore), or formation/deletion of ESEs with maximum splicing score. We also noted the formation/disruption of CpGs (tentative methylation sites of epigenetic sSNPs). All disease association studies involving sSNPs are also reported. Only 21/670 coding SNPs, mostly epigenetic, reported in 33 T2D-related genes, were found to be prominent coding synonymous. No prominent sSNPs have been recorded in three key T2D-related genes (GCGR, PPARGC1A, IGF1). Similarly, 20/366 coding synonymous were identified in ND related genes, mostly epigenetic. Meta-analysis showed that 17 of the above prominent sSNPs were previously investigated in association with various pathological conditions. Three out of four sSNPs (all epigenetic) were associated with T2D and one with NDs (branch site sSNP). Five were associated with other or related pathological conditions. None of the four sSNPs introducing new ESEs was found to be disease-associated. sSNPs introducing smaller Δscore changes (<1.5) in key proteins (INSR, IRS1, DISC1) were also correlated to pathological conditions. This data reveals that genetic variation in splicing-regulatory and particularly CpG sites might be related to disease predisposition and that in-silico analysis is useful for identifying sSNPs, which might be falsely identified as silent or synonymous.

Transcriptome-wide analysis of trypanosome mRNA decay reveals complex degradation kinetics and suggests a role for co-transcriptional degradation in determining mRNA levels

African trypanosomes are an excellent system for quantitative modelling of post-transcriptional mRNA control. Transcription is constitutive and polycistronic; individual mRNAs are excised by trans splicing and polyadenylation. We here measure mRNA decay kinetics in two life cycle stages, bloodstream and procyclic forms, by transcription inhibition and RNASeq. Messenger RNAs with short half-lives tend to show initial fast degradation, followed by a slower phase; they are often stabilized by depletion of the 5′–3′ exoribonuclease XRNA. Many longer-lived mRNAs show initial slow degradation followed by rapid destruction: we suggest that the slow phase reflects gradual deadenylation. Developmentally regulated mRNAs often show regulated decay, and switch their decay pattern. Rates of mRNA decay are good predictors of steady state levels for short mRNAs, but mRNAs longer than 3 kb show unexpectedly low abundances. Modelling shows that variations in splicing and polyadenylation rates can contribute to steady-state mRNA levels, but this is completely dependent on competition between processing and co-transcriptional mRNA precursor destruction.

SLaP mapper: a webserver for identifying and quantifying spliced-leader addition and polyadenylation site usage in kinetoplastid genomes

•

A web-server for identification of spliced-leader and polyadenylation addition sites.

•

Fully automated site quantification and gene assignment.

•

Multiple species within the Kinetoplastida.

The Kinetoplastida are a diverse and globally distributed class of free-living and parasitic single-celled eukaryotes that collectively cause a significant burden on human health and welfare. In kinetoplastids individual genes do not have promoters, but rather all genes are arranged downstream of a small number of RNA polymerase II transcription initiation sites and are thus transcribed in polycistronic gene clusters. Production of individual mRNAs from this continuous transcript occurs co-transcriptionally by trans-splicing of a ∼39 nucleotide capped RNA and subsequent polyadenylation of the upstream mRNA. SLaP mapper (Spliced-Leader and Polyadenylation mapper) is a fully automated web-service for identification, quantitation and gene-assignment of both spliced-leader and polyadenylation addition sites in Kinetoplastid genomes. SLaP mapper only requires raw read data from paired-end Illumina RNAseq and performs all read processing, mapping, quality control, quantification, and analysis in a fully automated pipeline. To provide usage examples and estimates of the quantity of sequence data required we use RNAseq obtained from two different library preparations from both Trypanosoma brucei and Leishmania mexicana to show the number of expected reads that are obtained from each preparation type. SLaP mapper is an easy to use, platform independent webserver that is freely available for use at http://www.stevekellylab.com/software/slap. Example files are provided on the website.

Capturing the variant surface glycoprotein repertoire (the VSGnome) of Trypanosoma brucei Lister 427

Trypanosoma brucei evades the adaptive immune response through the expression of antigenically distinct Variant Surface Glycoprotein (VSG) coats. To understand the progression and mechanisms of VSG switching, and to identify the VSGs expressed in populations of trypanosomes, it is desirable to predetermine the available repertoire of VSG genes (the 'VSGnome'). To date, the catalog of VSG genes present in any strain is far from complete and the majority of current information regarding VSGs is derived from the TREU927 strain that is not commonly used as an experimental model. We have assembled, annotated and analyzed 2563 distinct and previously unsequenced genes encoding complete and partial VSGs of the widely used Lister 427 strain of T. brucei. Around 80% of the VSGnome consists of incomplete genes or pseudogenes. Read-depth analysis demonstrated that most VSGs exist as single copies, but 360 exist as two or more indistinguishable copies. The assembled regions include five functional metacyclic VSG expression sites. One third of minichromosome sub-telomeres contain a VSG (64-67 VSGs on ∼96 minichromosomes), of which 85% appear to be functionally competent. The minichromosomal repertoire is very dynamic, differing among clones of the same strain. Few VSGs are unique along their entire length: frequent recombination events are likely to have shaped (and to continue to shape) the repertoire. In spite of their low sequence conservation and short window of expression, VSGs show evidence of purifying selection, with ∼40% of non-synonymous mutations being removed from the population. VSGs show a strong codon-usage bias that is distinct from that of any other group of trypanosome genes. VSG sequences are generally very divergent between Lister 427 and TREU927 strains of T. brucei, but those that are highly similar are not found in 'protected' genomic environments, but may reflect genetic exchange among populations.

Comparative ribosome profiling reveals extensive translational complexity in different Trypanosoma brucei life cycle stages

While gene expression is a fundamental and tightly controlled cellular process that is regulated at multiple steps, the exact contribution of each step remains unknown in any organism. The absence of transcription initiation regulation for RNA polymerase II in the protozoan parasite Trypanosoma brucei greatly simplifies the task of elucidating the contribution of translation to global gene expression. Therefore, we have sequenced ribosome-protected mRNA fragments in T. brucei, permitting the genome-wide analysis of RNA translation and translational efficiency. We find that the latter varies greatly between life cycle stages of the parasite and ∼100-fold between genes, thus contributing to gene expression to a similar extent as RNA stability. The ability to map ribosome positions at sub-codon resolution revealed extensive translation from upstream open reading frames located within 5' UTRs and enabled the identification of hundreds of previously un-annotated putative coding sequences (CDSs). Evaluation of existing proteomics and genome-wide RNAi data confirmed the translation of previously un-annotated CDSs and suggested an important role for >200 of those CDSs in parasite survival, especially in the form that is infective to mammals. Overall our data show that translational control plays a prevalent and important role in different parasite life cycle stages of T. brucei.

Strategies to construct null and conditional null Trypanosoma brucei mutants using Cre-recombinase and loxP

We describe two gene-knockout (KO) strategies in Trypanosoma brucei using Cre recombinase and loxP sites. Due to the limited number of selection markers for T. brucei, it has been difficult to generate a mutant with two genes knocked out and impractical to simultaneously knockout more than two genes, deterring detailed studies of important cellular mechanisms. The first KO strategy described can overcome the marker problem by allowing continuous re-use of drug-resistance markers. The same KO vector can be used to make a conditional KO system, when a gene of interest is essential for cell viability. As a gene of interest is removed from its original chromosomal locus by the induction of Cre recombinase, deletion is complete and instantaneous. This makes it easier to identify primary effects rather than having secondary effects obscuring phenotypic assessment, as is often the case with RNAi silencing.

The hnRNP F/H homologue of Trypanosoma brucei is differentially expressed in the two life cycle stages of the parasite and regulates splicing and mRNA stability

Trypanosomes are protozoan parasites that cycle between a mammalian host (bloodstream form) and an insect host, the Tsetse fly (procyclic stage). In trypanosomes, all mRNAs are trans-spliced as part of their maturation. Genome-wide analysis of trans-splicing indicates the existence of alternative trans-splicing, but little is known regarding RNA-binding proteins that participate in such regulation. In this study, we performed functional analysis of the Trypanosoma brucei heterogeneous nuclear ribonucleoproteins (hnRNP) F/H homologue, a protein known to regulate alternative splicing in metazoa. The hnRNP F/H is highly expressed in the bloodstream form of the parasite, but is also functional in the procyclic form. Transcriptome analyses of RNAi-silenced cells were used to deduce the RNA motif recognized by this protein. A purine rich motif, AAGAA, was enriched in both the regulatory regions flanking the 3′ splice site and poly (A) sites of the regulated genes. The motif was further validated using mini-genes carrying wild-type and mutated sequences in the 3′ and 5′ UTRs, demonstrating the role of hnRNP F/H in mRNA stability and splicing. Biochemical studies confirmed the binding of the protein to this proposed site. The differential expression of the protein and its inverse effects on mRNA level in the two lifecycle stages demonstrate the role of hnRNP F/H in developmental regulation.

Product feedback regulation implicated in translational control of the Trypanosoma brucei S-adenosylmethionine decarboxylase regulatory subunit prozyme

Human African sleeping sickness (HAT) is caused by the parasitic protozoan Trypanosoma brucei. Polyamine biosynthesis is an important drug target in the treatment of HAT. Previously we showed that trypanosomatid S-adenosylmethionine decarboxylase (AdoMetDC), a key enzyme for biosynthesis of the polyamine spermidine, is activated by heterodimer formation with an inactive paralogue termed prozyme. Furthermore, prozyme protein levels were regulated in response to reduced AdoMetDC activity. Herein we show that T. brucei encodes three prozyme transcripts. The 3'UTRs of these transcripts were mapped and chloramphenicol acetyltransferase (CAT) reporter constructs were used to identify a 1.2 kb region that contained a 3'UTR prozyme regulatory element sufficient to upregulate CAT protein levels (but not RNA) upon AdoMetDC inhibition, supporting the hypothesis that prozyme expression is regulated translationally. To gain insight into trans-acting factors, genetic rescue of AdoMetDC RNAi knock-down lines with human AdoMetDC was performed leading to rescue of the cell growth block, and restoration of prozyme protein to wild-type levels. Metabolite analysis showed that prozyme protein levels were inversely proportional to intracellular levels of decarboxylated AdoMet (dcAdoMet). These data suggest that prozyme translation may be regulated by dcAdoMet, a metabolite not previously identified to play a regulatory role.

Molecular and functional characterization of a Trypanosoma cruzi nuclear adenylate kinase isoform

Trypanosoma cruzi, the etiological agent of Chagas' disease, is an early divergent eukaryote in which control of gene expression relies mainly in post-transcriptional mechanisms. Transcription levels are globally up and down regulated during the transition between proliferating and non-proliferating life-cycle stages. In this work we characterized a nuclear adenylate kinase isoform (TcADKn) that is involved in ribosome biogenesis. Nuclear adenylate kinases have been recently described in a few organisms, being all related to RNA metabolism. Depending on active transcription and translation, TcADKn localizes in the nucleolus or the cytoplasm. A non-canonical nuclear localization signal was mapped towards the N-terminal of the protein, being the phosphate-binding loop essential for its localization. In addition, TcADKn nuclear exportation depends on the nuclear exportation adapter CRM1. TcADKn nuclear shuttling is governed by nutrient availability, oxidative stress and by the equivalent in T. cruzi of the mammalian TOR (Target of Rapamycin) pathway. One of the biological functions of TcADKn is ribosomal 18S RNA processing by direct interaction with ribosomal protein TcRps14. Finally, TcADKn expression is regulated by its 3′ UTR mRNA. Depending on extracellular conditions, cells modulate protein translation rates regulating ribosome biogenesis and nuclear adenylate kinases are probably key components in these processes.

Infection with Trypanosoma cruzi produces a condition known as Chagas disease which affects at least 17 million people. Adenylate kinases, so called myokinases, are involved in a wide variety of processes, mainly related to their role in nucleotide interconversion and energy management. Recently, nuclear isoforms have been described in several organisms. This “atypical” isoform in terms of primary structure was associated to ribosomes biogenesis in yeast and to Cajal body organization in humans. Moreover nuclear adenylate kinases are essential for maintaining cellular homeostasis. In this manuscript we characterized T. cruzi nuclear adenylate kinase (TcADKn). TcADKn localizes in the nucleolus or cell cytoplasm. Nuclear shuttling mechanisms were also studied for the first time, being dependent on nutrient availability, oxidative stress and by the equivalent of the mammalian TOR pathway in T. cruzi. Furthermore we characterized the signals involved in nuclear importation and exportation processes. In addition, TcADKn expression levels are regulated at an mRNA level, being its 3′UTR involved in this process. These findings are the first steps in the understanding of ribosome processing in trypanosomatids.

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.

How do trypanosomes change gene expression in response to the environment?

All organisms are able to modulate gene expression in response to internal and external stimuli. Trypanosomes represent a group that diverged early during the radiation of eukaryotes and do not utilise regulated initiation of transcription by RNA polymerase II. Here, the mechanisms present in trypanosomes to alter gene expression in response to stress and change of host environment are discussed and contrasted with those operating in yeast and cultured mammalian cells.

Inhibition of mRNA maturation in trypanosomes causes the formation of novel foci at the nuclear periphery containing cytoplasmic regulators of mRNA fate

Maturation of all cytoplasmic mRNAs in trypanosomes involves trans-splicing of a short exon at the 5′ end. Inhibition of trans-splicing results in an accumulation of partially processed oligocistronic mRNAs. Here, we show that the accumulation of newly synthesised partially processed mRNAs results in the formation of foci around the periphery of the nucleus. These nuclear periphery granules (NPGs) contain the full complement of P-body proteins identified in trypanosomes to date, as well as poly(A)-binding protein 2 and the trypanosome homologue of the RNA helicase VASA. NPGs resemble perinuclear germ granules from metazoa more than P-bodies because they: (1) are localised around the nuclear periphery; (2) are dependent on active transcription; (3) are not dissipated by cycloheximide; (4) contain VASA; and (5) depend on nuclear integrity. In addition, NPGs can be induced in cells depleted of the P-body core component SCD6. The description of NPGs in trypanosomes provides evidence that there is a perinuclear compartment that can determine the fate of newly transcribed mRNAs and that germ granules could be a specialised derivative.

Dual targeting of isoleucyl-tRNA synthetase in Trypanosoma brucei is mediated through alternative trans-splicing

Aminoacyl-tRNA synthetases catalyze the aminoacylation of tRNAs with their cognate amino acids. They are an essential part of each translation system and in eukaryotes are therefore found in both the cytosol and mitochondria. Thus, eukaryotes either have two distinct genes encoding the cytosolic and mitochondrial isoforms of each of these enzymes or a single gene encoding dually localized products. Trypanosomes require trans-splicing of a cap containing leader sequence onto the 5′-untranslated region of every mRNA. Recently we speculated that alternative trans-splicing could lead to the expression of proteins having amino-termini of different lengths that derive from the same gene. We now demonstrate that alternative trans-splicing, creating a long and a short spliced variant, is the mechanism for dual localization of trypanosomal isoleucyl-tRNA synthetase (IleRS). The protein product of the longer spliced variant possesses an amino-terminal presequence and is found exclusively in mitochondria. In contrast, the shorter spliced variant is translated to a cytosol-specific isoform lacking the presequence. Furthermore, we show that RNA stability is one mechanism determining the differential abundance of the two spliced isoforms.

Trans-splicing in trypanosomes: machinery and its impact on the parasite transcriptome

In trypanosomes, all RNAs are processed by the concerted action of trans-splicing and polyadenylation. In trans-splicing, a common spliced leader (SL) is donated to all mRNAs from a small RNA molecule, the SL RNA. This article summarizes recent findings in the field focusing on SL RNA transcription, cap modifications and pseudouridylation. The role(s) of these modifications for splicing and gene expression are discussed. The recruitment of SL RNA to the spliceosome depends on splicing factors and recent progress in identifying such factors is described. A recent major advance in understanding the role of trans-splicing in the trypanosome transcriptome was obtained by whole-genome mapping of the SL and polyadenylation sites, revealing surprising heterogeneity and suggesting that gene regulation, especially during cycling between the two hosts of the parasite, involves alternative trans-splicing. Finally, the SL silencing mechanism, which is harnessed by the parasite to control gene expression under stress, is discussed.

Ab initio identification of novel regulatory elements in the genome of Trypanosoma brucei by Bayesian inference on sequence segmentation

BACKGROUND

The rapid increase in the availability of genome information has created considerable demand for both comparative and ab initio predictive bioinformatic analyses. The biology laid bare in the genomes of many organisms is often novel, presenting new challenges for bioinformatic interrogation. A paradigm for this is the collected genomes of the kinetoplastid parasites, a group which includes Trypanosoma brucei the causative agent of human African trypanosomiasis. These genomes, though outwardly simple in organisation and gene content, have historically challenged many theories for gene expression regulation in eukaryotes.

METHODOLOGY/PRINCIPLE FINDINGS

Here we utilise a Bayesian approach to identify local changes in nucleotide composition in the genome of T. brucei. We show that there are several elements which are found at the starts and ends of multicopy gene arrays and that there are compositional elements that are common to all intergenic regions. We also show that there is a composition-inversion element that occurs at the position of the trans-splice site.

CONCLUSIONS/SIGNIFICANCE

The nature of the elements discovered reinforces the hypothesis that context dependant RNA secondary structure has an important influence on gene expression regulation in Trypanosoma brucei.

Gene expression in Trypanosoma brucei: lessons from high-throughput RNA sequencing

Trypanosoma brucei undergoes major biochemical and morphological changes during its development from the bloodstream form in the mammalian host to the procyclic form in the midgut of its insect host. The underlying regulation of gene expression, however, is poorly understood. More than 60% of the predicted genes remain annotated as hypothetical, and the 5' and 3' untranslated regions important for regulation of gene expression are unknown for >90% of the genes. In this review, we compare the data from four recently published high-throughput RNA sequencing studies in light of the different experimental setups and discuss how these data can enhance genome annotation and give insights into the regulation of gene expression in T. brucei.

Cell biology of the trypanosome genome

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

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.

The transcriptome of the human pathogen Trypanosoma brucei at single-nucleotide resolution

The genome of Trypanosoma brucei, the causative agent of African trypanosomiasis, was published five years ago, yet identification of all genes and their transcripts remains to be accomplished. Annotation is challenged by the organization of genes transcribed by RNA polymerase II (Pol II) into long unidirectional gene clusters with no knowledge of how transcription is initiated. Here we report a single-nucleotide resolution genomic map of the T. brucei transcriptome, adding 1,114 new transcripts, including 103 non-coding RNAs, confirming and correcting many of the annotated features and revealing an extensive heterogeneity of 5′ and 3′ ends. Some of the new transcripts encode polypeptides that are either conserved in T. cruzi and Leishmania major or were previously detected in mass spectrometry analyses. High-throughput RNA sequencing (RNA-Seq) was sensitive enough to detect transcripts at putative Pol II transcription initiation sites. Our results, as well as recent data from the literature, indicate that transcription initiation is not solely restricted to regions at the beginning of gene clusters, but may occur at internal sites. We also provide evidence that transcription at all putative initiation sites in T. brucei is bidirectional, a recently recognized fundamental property of eukaryotic promoters. Our results have implications for gene expression patterns in other important human pathogens with similar genome organization (Trypanosoma cruzi, Leishmania sp.) and revealed heterogeneity in pre-mRNA processing that could potentially contribute to the survival and success of the parasite population in the insect vector and the mammalian host.

Identifying genes essential for survival in the host is fundamental to unraveling the biology of human pathogens and understanding mechanisms of pathogenesis. The protozoan parasite Trypanosoma brucei causes devastating diseases in humans and animals in sub-Saharan Africa, and the publication in 2005 of the genome sequence provided the first glance at the coding potential of this organism. Although at present there is a catalogue of predicted protein coding genes, the challenge remains to identify all authentic genes, including their boundaries. We used next generation RNA sequencing (RNA-Seq) to map transcribed regions and RNA polymerase II transcription initiation sites on a genome-wide scale. This approach allowed us to improve and correct the current annotation, to reveal a widespread heterogeneity of RNA processing sites (trans-splicing and polyadenylation) and to estimate that most genes are expressed at levels corresponding to 1 to 10 mRNAs per cell. Our data indicate that different transcript forms representing the same gene are present stochastically within the mRNA population. This unanticipated scenario may contribute to determining gene expression landscapes to adapt to different environments in the parasite life cycle.

Spliced leader trapping reveals widespread alternative splicing patterns in the highly dynamic transcriptome of Trypanosoma brucei

Trans-splicing of leader sequences onto the 5′ends of mRNAs is a widespread phenomenon in protozoa, nematodes and some chordates. Using parallel sequencing we have developed a method to simultaneously map 5′splice sites and analyze the corresponding gene expression profile, that we term spliced leader trapping (SLT). The method can be applied to any organism with a sequenced genome and trans-splicing of a conserved leader sequence. We analyzed the expression profiles and splicing patterns of bloodstream and insect forms of the parasite Trypanosoma brucei. We detected the 5′ splice sites of 85% of the annotated protein-coding genes and, contrary to previous reports, found up to 40% of transcripts to be differentially expressed. Furthermore, we discovered more than 2500 alternative splicing events, many of which appear to be stage-regulated. Based on our findings we hypothesize that alternatively spliced transcripts present a new means of regulating gene expression and could potentially contribute to protein diversity in the parasite. The entire dataset can be accessed online at TriTrypDB or through: http://splicer.unibe.ch/.

Some organisms like the human and animal parasite Trypanosoma brucei add a leader sequence to their mRNAs through a reaction called trans-splicing. Until now the splice sites for most mRNAs were unknown in T. brucei. Using high throughput sequencing we have developed a method to identify the splice sites and at the same time measure the abundance of the corresponding mRNAs. Analyzing three different life cycle stages of the parasite we identified the vast majority of splice sites in the organism and, to our great surprise, uncovered more than 2500 alternative splicing events, many of which appeared to be specific for one of the life cycle stages. Alternative splicing is a result of the addition of the leader sequence to different positions on the mRNA, leading to mixed mRNA populations that can encode for proteins with varying properties. One of the most obvious changes caused by alternative splicing is the gain or loss of targeting signals, leading to differential localization of the corresponding proteins. Based on our findings we hypothesize that alternative splicing is a major mechanism to regulate gene expression in T. brucei and could contribute to protein diversity in the parasite.

The pre-mRNA splicing machinery of trypanosomes: complex or simplified?

Trypanosomatids are early-diverged, protistan parasites of which Trypanosoma brucei, Trypanosoma cruzi, and several species of Leishmania cause severe, often lethal diseases in humans. To better combat these parasites, their molecular biology has been a research focus for more than 3 decades, and the discovery of spliced leader (SL) trans splicing in T. brucei established a key difference between parasites and hosts. In SL trans splicing, the capped 5'-terminal region of the small nuclear SL RNA is fused onto the 5' end of each mRNA. This process, in conjunction with polyadenylation, generates individual mRNAs from polycistronic precursors and creates functional mRNA by providing the cap structure. The reaction is a two-step transesterification process analogous to intron removal by cis splicing which, in trypanosomatids, is confined to very few pre-mRNAs. Both types of pre-mRNA splicing are carried out by the spliceosome, consisting of five U-rich small nuclear RNAs (U snRNAs) and, in humans, up to approximately 170 different proteins. While trypanosomatids possess a full set of spliceosomal U snRNAs, only a few splicing factors were identified by standard genome annotation because trypanosomatid amino acid sequences are among the most divergent in the eukaryotic kingdom. This review focuses on recent progress made in the characterization of the splicing factor repertoire in T. brucei, achieved by tandem affinity purification of splicing complexes, by systematic analysis of proteins containing RNA recognition motifs, and by mining the genome database. In addition, recent findings about functional differences between trypanosome and human pre-mRNA splicing factors are discussed.

Posttranscriptional control and the role of RNA-binding proteins in gene regulation in trypanosomatid protozoan parasites

Trypanosomatids are unicellular eukaryotes responsible for severe diseases in humans. They exhibit a number of remarkable biological phenomena, especially at the RNA level. During their life cycles, they alternate between a mammalian host and an insect vector and undergo profound biochemical and morphological transformations in order to adapt to the different environments they find within one or the other host species. These changes are orchestrated by specific gene expression programs. In contrast to other organisms, trypanosomatids do not regulate RNA polymerase II-dependent transcription initiation. Evidence so far indicates that the main control points in gene expression are mRNA degradation and translation. Recent studies have shown that RNA-binding proteins (RBPs) play a critical role in the developmental regulation of mRNA and protein abundance. RBPs seem to bind to specific subsets of mRNAs encoding functionally related proteins. These ribonucleoprotein complexes may represent posttranscriptional operons or regulons that are able to control the fate of multiple mRNAs simultaneously. We suggest that trypanosomatids transduce environmental signals into mRNA and protein abundance through posttranslational modification of RBPs.

Genome-wide analysis of mRNA abundance in two life-cycle stages of Trypanosoma brucei and identification of splicing and polyadenylation sites

Transcription of protein-coding genes in trypanosomes is polycistronic and gene expression is primarily regulated by post-transcriptional mechanisms. Sequence motifs in the untranslated regions regulate mRNA trans-splicing and RNA stability, yet where UTRs begin and end is known for very few genes. We used high-throughput RNA-sequencing to determine the genome-wide steady-state mRNA levels (‘transcriptomes’) for ∼90% of the genome in two stages of the Trypanosoma brucei life cycle cultured in vitro. Almost 6% of genes were differentially expressed between the two life-cycle stages. We identified 5′ splice-acceptor sites (SAS) and polyadenylation sites (PAS) for 6959 and 5948 genes, respectively. Most genes have between one and three alternative SAS, but PAS are more dispersed. For 488 genes, SAS were identified downstream of the originally assigned initiator ATG, so a subsequent in-frame ATG presumably designates the start of the true coding sequence. In some cases, alternative SAS would give rise to mRNAs encoding proteins with different N-terminal sequences. We could identify the introns in two genes known to contain them, but found no additional genes with introns. Our study demonstrates the usefulness of the RNA-seq technology to study the transcriptional landscape of an organism whose genome has not been fully annotated.

Processing of a phosphoglycerate kinase reporter mRNA in Trypanosoma brucei is not coupled to transcription by RNA polymerase II

Capping of mRNAs is strictly coupled to RNA polymerase II transcription and there is evidence, mainly from metazoans, that other steps in pre-mRNA processing show a similar linkage. In trypanosomes, however, the mRNA cap is supplied by a trans spliced leader sequence. Thus pre-mRNAs transcribed by RNA Polymerase I are capped by trans splicing, and translation-competent transgenic mRNAs can be produced by RNA Polymerase I and T7 RNA polymerase so long as the primary transcript has a splice acceptor signal. We quantified the efficiency of processing of trypanosome pre-mRNAs produced from a plasmid integrated either at the tubulin locus, or in an rRNA spacer, and transcribed by RNA polymerase II, RNA polymerase I or T7 RNA polymerase. The processing efficiencies were similar for primary transcripts from the tubulin locus, produced by RNA polymerase II, and for RNA from an rRNA spacer, transcribed by RNA polymerase I. Primary transcripts produced by T7 RNA polymerase from the tubulin locus were processed almost as well. There was therefore no evidence for recruitment of the 3'-splicing apparatus by the RNA polymerase. Abundant transcripts transcribed from the rRNA locus by T7 RNA polymerase were somewhat less efficiently processed.

Establishment of an in vitro trans-splicing system in Trypanosoma brucei that requires endogenous spliced leader RNA

In trypanosomes a 39 nucleotide exon, the spliced leader (SL) is donated to all mRNAs from a small RNA, the SL RNA, by trans-splicing. Since the discovery of trans-splicing in trypanosomes two decades ago, numerous attempts failed to reconstitute the reaction in vitro. In this study, a crude whole-cell extract utilizing the endogenous SL RNA and synthetic tubulin pre-mRNA were used to reconstitute the trans-splicing reaction. An RNase protection assay was used to detect the trans-spliced product. The reaction was optimized and shown to depend on ATP and intact U2 and U6 snRNPs. Mutations introduced at the polypyrimidine tract and the AG splice site reduced the reaction efficiency. To simplify the assay, RT–PCR and quantitative real-time PCR assays were established. The system was used to examine the structural requirements for SL RNA as a substrate in the reaction. Interestingly, synthetic SL RNA assembled poorly to its cognate particle and was not utilized in the reaction. However, SL RNA synthesized in cells lacking Sm proteins, which is defective in cap-4 modification, was active in the reaction. This study is the first step towards further elucidating the mechanism of trans-splicing, an essential reaction which determines the trypanosome transcriptome.

Gene expression in trypanosomatid parasites

The parasites Leishmania spp., Trypanosoma brucei, and Trypanosoma cruzi are the trypanosomatid protozoa that cause the deadly human diseases leishmaniasis, African sleeping sickness, and Chagas disease, respectively. These organisms possess unique mechanisms for gene expression such as constitutive polycistronic transcription of protein-coding genes and trans-splicing. Little is known about either the DNA sequences or the proteins that are involved in the initiation and termination of transcription in trypanosomatids. In silico analyses of the genome databases of these parasites led to the identification of a small number of proteins involved in gene expression. However, functional studies have revealed that trypanosomatids have more general transcription factors than originally estimated. Many posttranslational histone modifications, histone variants, and chromatin modifying enzymes have been identified in trypanosomatids, and recent genome-wide studies showed that epigenetic regulation might play a very important role in gene expression in this group of parasites. Here, we review and comment on the most recent findings related to transcription initiation and termination in trypanosomatid protozoa.

Genome-wide computational identification of functional RNA elements in Trypanosoma brucei

Background

Post-transcriptional regulation of gene expression is the dominant regulatory mechanism in trypanosomatids as their mRNAs are transcribed from polycistronic units. A few cis-acting RNA elements in 3'-untranslated regions of mRNAs have been identified in trypanosomatids, which affect the mRNA stability or translation rate in different life stages of these parasites. Other functional RNAs (fRNAs) also play essential roles in these organisms. However, there has been no genome-wide analysis for identification of fRNAs in trypanosomatids.

Results

Functional RNAs, including non-coding RNAs (ncRNAs) and cis-acting RNA elements involved in post-transcriptional gene regulation, were predicted based on two independent computational analyses of the genome of Trypanosoma brucei. In the first analysis, the predicted candidate ncRNAs were identified based on conservation with the related trypanosomatid Leishmania braziliensis. This prediction had a substantially low estimated false discovery rate, and a considerable number of the predicted ncRNAs represented novel classes with unknown functions. In the second analysis, we identified a number of function-specific regulatory motifs, based on which we devised a classifier that can be used for homology-independent function prediction in T. brucei.

Conclusion

This first genome-wide analysis of fRNAs in trypanosomatids restricts the search space of experimental approaches and, thus, can significantly expedite the process of characterization of these elements. Our classifier for function prediction based on cis-acting regulatory elements can also, in combination with other methods, provide the means for homology-independent annotation of trypanosomatid genomes.