Messenger RNA processing sites 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.

The SPARC complex defines RNAPII promoters in Trypanosoma brucei

Kinetoplastids are a highly divergent lineage of eukaryotes with unusual mechanisms for regulating gene expression. We previously surveyed 65 putative chromatin factors in the kinetoplastid Trypanosoma brucei. Our analyses revealed that the predicted histone methyltransferase SET27 and the Chromodomain protein CRD1 are tightly concentrated at RNAPII transcription start regions (TSRs). Here, we report that SET27 and CRD1, together with four previously uncharacterized constituents, form the SET27 promoter-associated regulatory complex (SPARC), which is specifically enriched at TSRs. SET27 loss leads to aberrant RNAPII recruitment to promoter sites, accumulation of polyadenylated transcripts upstream of normal transcription start sites, and conversion of some normally unidirectional promoters to bidirectional promoters. Transcriptome analysis in the absence of SET27 revealed upregulated mRNA expression in the vicinity of SPARC peaks within the main body of chromosomes in addition to derepression of genes encoding variant surface glycoproteins (VSGs) located in subtelomeric regions. These analyses uncover a novel chromatin-associated complex required to establish accurate promoter position and directionality.

Comparative molecular cell biology of phototrophic euglenids and parasitic trypanosomatids sheds light on the ancestor of Euglenozoa

Parasitic trypanosomatids and phototrophic euglenids are among the most extensively studied euglenozoans. The phototrophic euglenid lineage arose relatively recently through secondary endosymbiosis between a phagotrophic euglenid and a prasinophyte green alga that evolved into the euglenid secondary chloroplast. The parasitic trypanosomatids (i.e. Trypanosoma spp. and Leishmania spp.) and the freshwater phototrophic euglenids (i.e. Euglena gracilis) are the most evolutionary distant lineages in the Euglenozoa phylogenetic tree. The molecular and cell biological traits they share can thus be considered as ancestral traits originating in the common euglenozoan ancestor. These euglenozoan ancestral traits include common mitochondrial presequence motifs, respiratory chain complexes containing various unique subunits, a unique ATP synthase structure, the absence of mitochondria-encoded transfer RNAs (tRNAs), a nucleus with a centrally positioned nucleolus, closed mitosis without dissolution of the nuclear membrane and nucleoli, a nuclear genome containing the unusual 'J' base (β-D-glucosyl-hydroxymethyluracil), processing of nucleus-encoded precursor messenger RNAs (pre-mRNAs) via spliced-leader RNA (SL-RNA) trans-splicing, post-transcriptional gene silencing by the RNA interference (RNAi) pathway and the absence of transcriptional regulation of nuclear gene expression. Mitochondrial uridine insertion/deletion RNA editing directed by guide RNAs (gRNAs) evolved in the ancestor of the kinetoplastid lineage. The evolutionary origin of other molecular features known to be present only in either kinetoplastids (i.e. polycistronic transcripts, compaction of nuclear genomes) or euglenids (i.e. monocistronic transcripts, huge genomes, many nuclear cis-spliced introns, polyproteins) is unclear.

The polyadenylation complex of Trypanosoma brucei: Characterization of the functional poly(A) polymerase

The generation of mature mRNA in the protozoan parasite Trypanosoma brucei requires coupled polyadenylation and trans splicing. In contrast to other eukaryotes, we still know very little on components, mechanisms, and dynamics of the 3' end-processing machinery in trypanosomes. To characterize the catalytic core of the polyadenylation complex in T. brucei, we first identified the poly(A) polymerase [Tb927.7.3780] as the major functional, nuclear-localized enzyme in trypanosomes. In contrast, another poly(A) polymerase, encoded by an intron-containing gene [Tb927.3.3160], localizes mainly in the cytoplasm and appears not to be functional in general 3' end processing of mRNAs. Based on tandem-affinity purification with tagged CPSF160 and mass spectrometry, we identified ten associated components of the trypanosome polyadenylation complex, including homologues to all four CPSF subunits, Fip1, CstF50/64, and Symplekin, as well as two hypothetical proteins. RNAi-mediated knockdown revealed that most of these factors are essential for growth and required for both in vivo polyadenylation and trans splicing, arguing for a general coupling of these two mRNA-processing reactions.

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.

The calmodulin intergenic spacer as molecular target for characterization of Leishmania species

Background

Human leishmaniasis is a neglected disease caused by parasites of the genus Leishmania. Clinical aspects of this disease can vary significantly, reflecting the wide range of parasites in the genus Leishmania. Knowing accurately the Leishmania species infecting humans is important for clinical case management and evaluation of epidemiological risk. Calmodulin is an essential gene in trypanosomatids that modulates the calcium metabolism in various cellular activities. Despite its strong conservation in trypanosomatids, it has been recently observed that its untranslated regions (UTR) diverge among species.

Methods

In this study we analyzed the sequences and the absolute dinucleotide frequency of the intergenic spacer of the calmodulin gene (containing both, 3′ and 5′UTR) in nine reference Leishmania species and ten clinical isolates obtained from patients with cutaneous leishmaniasis.

Results

We show that the short calmodulin intergenic spacers exhibit features that make them interesting for applications in molecular characterization and phylogenetic studies of Leishmania. Dendrograms based on sequence alignments and on the dinucleotide frequency indicate that this particular region of calmodulin gene might be useful for species typing between the Leishmania and Viannia subgenera.

Conclusions

Mutations and composition of the calmodulin intergenic spacer from Leishmania species might have taxonomic value as parameters to define if an isolate is identical to a certain species or belongs to one of the two current subgenera.

CMF22 is a broadly conserved axonemal protein and is required for propulsive motility in Trypanosoma brucei

The eukaryotic flagellum (or cilium) is a broadly conserved organelle that provides motility for many pathogenic protozoa and is critical for normal development and physiology in humans. Therefore, defining core components of motile axonemes enhances understanding of eukaryotic biology and provides insight into mechanisms of inherited and infectious diseases in humans. In this study, we show that component of motile flagella 22 (CMF22) is tightly associated with the flagellar axoneme and is likely to have been present in the last eukaryotic common ancestor. The CMF22 amino acid sequence contains predicted IQ and ATPase associated with a variety of cellular activities (AAA) motifs that are conserved among CMF22 orthologues in diverse organisms, hinting at the importance of these domains in CMF22 function. Knockdown by RNA interference (RNAi) and rescue with an RNAi-immune mRNA demonstrated that CMF22 is required for propulsive cell motility in Trypanosoma brucei. Loss of propulsive motility in CMF22-knockdown cells was due to altered flagellar beating patterns, rather than flagellar paralysis, indicating that CMF22 is essential for motility regulation and likely functions as a fundamental regulatory component of motile axonemes. CMF22 association with the axoneme is weakened in mutants that disrupt the nexin-dynein regulatory complex, suggesting potential interaction with this complex. Our results provide insight into the core machinery required for motility of eukaryotic flagella.

Genome-wide analysis of 3'-untranslated regions supports the existence of post-transcriptional regulons controlling gene expression in trypanosomes

In eukaryotic cells, a group of messenger ribonucleic acids (mRNAs) encoding functionally interrelated proteins together with the trans-acting factors that coordinately modulate their expression is termed a post-transcriptional regulon, due to their partial analogy to a prokaryotic polycistron. This mRNA clustering is organized by sequence-specific RNA-binding proteins (RBPs) that bind cis-regulatory elements in the noncoding regions of genes, and mediates the synchronized control of their fate. These recognition motifs are often characterized by conserved sequences and/or RNA structures, and it is likely that various classes of cis-elements remain undiscovered. Current evidence suggests that RNA regulons govern gene expression in trypanosomes, unicellular parasites which mainly use post-transcriptional mechanisms to control protein synthesis. In this study, we used motif discovery tools to test whether groups of functionally related trypanosomatid genes contain a common cis-regulatory element. We obtained conserved structured RNA motifs statistically enriched in the noncoding region of 38 out of 53 groups of metabolically related transcripts in comparison with a random control. These motifs have a hairpin loop structure, a preferred sense orientation and are located in close proximity to the open reading frames. We found that 15 out of these 38 groups represent unique motifs in which most 3'-UTR signature elements were group-specific. Two extensively studied Trypanosoma cruzi RBPs, TcUBP1 and TcRBP3 were found associated with a few candidate RNA regulons. Interestingly, 13 motifs showed a strong correlation with clusters of developmentally co-expressed genes and six RNA elements were enriched in gene clusters affected after hyperosmotic stress. Here we report a systematic genome-wide in silico screen to search for novel RNA-binding sites in transcripts, and describe an organized network of several coordinately regulated cohorts of mRNAs in T. cruzi. Moreover, we found that structured RNA elements are also conserved in other human pathogens. These results support a model of regulation of gene expression by multiple post-transcriptional regulons in trypanosomes.

The phosphoarginine energy-buffering system of trypanosoma brucei involves multiple arginine kinase isoforms with different subcellular locations

Phosphagen energy-buffering systems play an essential role in regulating the cellular energy homeostasis in periods of high-energy demand or energy supply fluctuations. Here we describe the phosphoarginine/arginine kinase system of the kinetoplastid parasite Trypanosoma brucei, consisting of three highly similar arginine kinase isoforms (TbAK1-3). Immunofluorescence microscopy using myc-tagged protein versions revealed that each isoform is located in a specific subcellular compartment: TbAK1 is exclusively found in the flagellum, TbAK2 in the glycosome, and TbAK3 in the cytosol of T. brucei. The flagellar location of TbAK1 is dependent on a 22 amino acid long N-terminal sequence, which is sufficient for targeting a GFP-fusion protein to the trypanosome flagellum. The glycosomal location of TbAK2 is in agreement with the presence of a conserved peroxisomal targeting signal, the C-terminal tripeptide ‘SNL’. TbAK3 lacks any apparent targeting sequences and is accordingly located in the cytosol of the parasite. Northern blot analysis indicated that each TbAK isoform is differentially expressed in bloodstream and procyclic forms of T. brucei, while the total cellular arginine kinase activity was 3-fold higher in bloodstream form trypanosomes. These results suggest a substantial change in the temporal and spatial energy requirements during parasite differentiation. Increased arginine kinase activity improved growth of procyclic form T. brucei during oxidative challenges with hydrogen peroxide. Elimination of the total cellular arginine kinase activity by RNA interference significantly decreased growth (>90%) of procyclic form T. brucei under standard culture conditions and was lethal for this life cycle stage in the presence of hydrogen peroxide. The putative physiological roles of the different TbAK isoforms in T. brucei are further discussed.

Footprints of a trypanosomatid RNA world: pre-small subunit rRNA processing by spliced leader addition trans-splicing

The addition of a capped mini-exon [spliced leader (SL)] through trans-splicing is essential for the maturation of RNA polymerase (pol) II-transcribed polycistronic pre-mRNAs in all members of the Trypanosomatidae family. This process is an inter-molecular splicing reaction that follows the same basic rules of cis-splicing reactions. In this study, we demonstrated that mini-exons were added to precursor ribosomal RNA (pre-rRNA) are transcribed by RNA pol I, including the 5' external transcribed spacer (ETS) region. Additionally, we detected the SL-5'ETS molecule using three distinct methods and located the acceptor site between two known 5'ETS rRNA processing sites (A' and A1) in four different trypanosomatids. Moreover, we detected a polyadenylated 5'ETS upstream of the trans-splicing acceptor site, which also occurs in pre-mRNA trans-splicing. After treatment with an indirect trans-splicing inhibitor (sinefungin), we observed SL-5'ETS decay. However, treatment with 5-fluorouracil (a precursor of RNA synthesis that inhibits the degradation of pre-rRNA) led to the accumulation of SL-5'ETS, suggesting that the molecule may play a role in rRNA degradation. The detection of trans-splicing in these molecules may indicate broad RNA-joining properties, regardless of the polymerase used for transcription.

Identification of SL addition trans-splicing acceptor sites in the internal transcribed spacer I region of pre-rRNA in Leishmania (Leishmania) amazonensis

Trypanosomatidae is a family of early branching eukaryotes harbouring a distinctive repertoire of gene expression strategies. Functional mature messenger RNA is generated via the trans-splicing and polyadenylation processing of constitutively transcribed polycistronic units. Recently, trans-splicing of pre-small subunit ribosomal RNA in the 5' external transcribed spacer region and of precursor tRNAsec have been described. Here, we used a previously validated semi-nested reverse transcription-polymerase chain reaction strategy to investigate internal transcribed spacer (ITS) I acceptor sites in total RNA from Leishmania (Leishmania) amazonensis. Two distinct spliced leader-containing RNAs were detected indicating that trans-splicing reactions occur at two AG acceptor sites mapped in this ITS region. These data provide further evidence of the wide spectrum of RNA molecules that act as trans-splicing acceptors in trypanosomatids.

Trypanosomatid comparative genomics: Contributions to the study of parasite biology and different parasitic diseases

In 2005, draft sequences of the genomes of Trypanosoma brucei, Trypanosoma cruzi and Leishmania major, also known as the Tri-Tryp genomes, were published. These protozoan parasites are the causative agents of three distinct insect-borne diseases, namely sleeping sickness, Chagas disease and leishmaniasis, all with a worldwide distribution. Despite the large estimated evolutionary distance among them, a conserved core of ~6,200 trypanosomatid genes was found among the Tri-Tryp genomes. Extensive analysis of these genomic sequences has greatly increased our understanding of the biology of these parasites and their host-parasite interactions. In this article, we review the recent advances in the comparative genomics of these three species. This analysis also includes data on additional sequences derived from other trypanosmatid species, as well as recent data on gene expression and functional genomics. In addition to facilitating the identification of key parasite molecules that may provide a better understanding of these complex diseases, genome studies offer a rich source of new information that can be used to define potential new drug targets and vaccine candidates for controlling these parasitic infections.

The genome and its implications

Trypanosoma cruzi has a heterogeneous population composed of a pool of strains that circulate in the domestic and sylvatic cycles. Genome sequencing of the clone CL Brener revealed a highly repetitive genome of about 110Mb containing an estimated 22,570 genes. Because of its hybrid nature, sequences representing the two haplotypes have been generated. In addition, a repeat content close to 50% made the assembly of the estimated 41 pairs of chromosomes quite challenging. Similar to other trypanosomatids, the organization of T. cruzi chromosomes was found to be very peculiar, with protein-coding genes organized in long polycistronic transcription units encoding 20 or more proteins in one strand separated by strand switch regions. Another remarkable feature of the T. cruzi genome is the massive expansion of surface protein gene families. Because of the high genetic diversity of the T. cruzi population, sequencing of additional strains and comparative genomic and transcriptome analyses are in progress. Five years after its publication, the genome data have proven to be an essential tool for the study of T. cruzi and increasing efforts to translate this knowledge into the development of new modes of intervention to control Chagas disease are underway.

Hyperosmotic stress induces aquaporin-dependent cell shrinkage, polyphosphate synthesis, amino acid accumulation, and global gene expression changes in Trypanosoma cruzi

The protist parasite Trypanosoma cruzi has evolved the ability to transit between completely different hosts and to replicate in adverse environments. In particular, the epimastigote form, the replicative stage inside the vector, is subjected to nutritional and osmotic stresses during its development. In this work, we describe the biochemical and global gene expression changes of epimastigotes under hyperosmotic conditions. Hyperosmotic stress resulted in cell shrinking within a few minutes. Depending on the medium osmolarity, this was followed by lack of volume recovery for at least 2 h or by slow recovery. Experiments with inhibitors, or with cells in which an aquaporin gene (TcAQP1) was knocked down or overexpressed, revealed its importance for the cellular response to hyperosmotic stress. Furthermore, the adaptation to this new environment was shown to involve the regulation of the polyphosphate polymerization state as well as changes in amino acid catabolism to generate compatible osmolytes. A genome-wide transcriptional analysis of stressed parasites revealed down-regulation of genes belonging to diverse functional categories and up-regulation of genes encoding trans-sialidase-like and ribosomal proteins. Several of these changes were confirmed by Northern blot analyses. Sequence analysis of the 3'UTRs of up- and down-regulated genes allowed the identification of conserved structural RNA motifs enriched in each group, suggesting that specific ribonucleoprotein complexes could be of great importance in the adaptation of this parasite to different environments through regulation of transcript abundance.

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.

Transcription of long hypothetical orfs in Trypanosoma cruzi: the epimastigote stage uses trans-splicing sites that generate short 5' UTRs

We mapped the 5' UTR for five long hypothetical orfs from Trypanosoma cruzi; each one having a length of more than 10,000 bp. Our aim was to verify the constraints to the length of the 5' UTR and to identify the sites of alternative trans-splicing in the epimastigote stage of three T. cruzi strains. We used reverse transcription PCR to amplify the 5' UTR and demonstrated the transcription of all selected genes as well as additional trans-splicing sites in two of these genes. We observed that the length of the 5' UTR in these genes has a limit, in contrast to previous reports that indicated a trend for longer genes to display a proportionally long 5' UTR. The maximum length of the 5' UTR for the long genes analyzed in the present work is approximately 3% of the orf and, on average, is 1% of the orf length. The poly-pyrimidine tracts used as trans-splicing signal are in the range of 17-53 bases within a distance of 6-59 nt to first spliced-leader acceptor site. T. cruzi populations may use both signals differentially. We conclude that the limit for the 5' UTR length in long genes is determined primarily by the distance to neighboring genes.

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.

The cell cycle regulated transcriptome of Trypanosoma brucei

Progression of the eukaryotic cell cycle requires the regulation of hundreds of genes to ensure that they are expressed at the required times. Integral to cell cycle progression in yeast and animal cells are temporally controlled, progressive waves of transcription mediated by cell cycle-regulated transcription factors. However, in the kinetoplastids, a group of early-branching eukaryotes including many important pathogens, transcriptional regulation is almost completely absent, raising questions about the extent of cell-cycle regulation in these organisms and the mechanisms whereby regulation is achieved. Here, we analyse gene expression over the Trypanosoma brucei cell cycle, measuring changes in mRNA abundance on a transcriptome-wide scale. We developed a “double-cut” elutriation procedure to select unperturbed, highly synchronous cell populations from log-phase cultures, and compared this to synchronization by starvation. Transcriptome profiling over the cell cycle revealed the regulation of at least 430 genes. While only a minority were homologous to known cell cycle regulated transcripts in yeast or human, their functions correlated with the cellular processes occurring at the time of peak expression. We searched for potential target sites of RNA-binding proteins in these transcripts, which might earmark them for selective degradation or stabilization. Over-represented sequence motifs were found in several co-regulated transcript groups and were conserved in other kinetoplastids. Furthermore, we found evidence for cell-cycle regulation of a flagellar protein regulon with a highly conserved sequence motif, bearing similarity to consensus PUF-protein binding motifs. RNA sequence motifs that are functional in cell-cycle regulation were more widespread than previously expected and conserved within kinetoplastids. These findings highlight the central importance of post-transcriptional regulation in the proliferation of parasitic kinetoplastids.

Structure-function analysis of dynein light chain 1 identifies viable motility mutants in bloodstream-form Trypanosoma brucei

The flagellum of Trypanosoma brucei is an essential and multifunctional organelle that is receiving increasing attention as a potential drug target and as a system for studying flagellum biology. RNA interference (RNAi) knockdown is widely used to test the requirement for a protein in flagellar motility and has suggested that normal flagellar motility is essential for viability in bloodstream-form trypanosomes. However, RNAi knockdown alone provides limited functional information because the consequence is often loss of a multiprotein complex. We therefore developed an inducible system that allows functional analysis of point mutations in flagellar proteins in T. brucei. Using this system, we identified point mutations in the outer dynein light chain 1 (LC1) that allow stable assembly of outer dynein motors but do not support propulsive motility. In procyclic-form trypanosomes, the phenotype of LC1 mutants with point mutations differs from the motility and structural defects of LC1 knockdowns, which lack the outer-arm dynein motor. Thus, our results distinguish LC1-specific functions from broader functions of outer-arm dynein. In bloodstream-form trypanosomes, LC1 knockdown blocks cell division and is lethal. In contrast, LC1 point mutations cause severe motility defects without affecting viability, indicating that the lethal phenotype of LC1 RNAi knockdown is not due to defective motility. Our results demonstrate for the first time that normal motility is not essential in bloodstream-form T. brucei and that the presumed connection between motility and viability is more complex than might be interpreted from knockdown studies alone. These findings open new avenues for dissecting mechanisms of flagellar protein function and provide an important step in efforts to exploit the potential of the flagellum as a therapeutic target in African sleeping sickness.

Development of a dual reporter system to identify regulatory cis-acting elements in untranslated regions of Trypanosoma cruzi mRNAs

In trypanosomatids, transcription is polycistronic and gene expression control occurs mainly at the post-transcriptional level. To investigate the role of sequences present in the 3'UTR of stage-specific mRNAs of Trypanosoma cruzi, we generated a new vector, named pTcDUALuc, containing the firefly and Renilla luciferase reporter genes. To test this vector, sequences derived from the 3'UTR plus intergenic regions of the alpha tubulin gene, which is up-regulated in epimastigotes, and amastin, which is up-regulated in amastigotes, were inserted downstream from the firefly reporter gene and luciferase activity was compared in transient and stable transfected parasites. As expected, increased luciferase activity was detected in epimastigotes transiently transfected with pTcDUALuc containing tubulin sequences. Using stable transfected cell lines that were allowed to differentiate into amastigotes, we observed increased luciferase activity and mRNA levels in amastigotes transfected with pTcDUALuc containing amastin sequences. We also showed that the spliced leader sequence and poly-A tail were inserted in the predicted sites of the firefly luciferase mRNA and that deletions in the alpha tubulin 3'UTR resulted in decreased luciferase expression because it affects polyadenylation. In contrast to the constructs containing 3'UTR sequences derived from tubulin and amastin genes, the presence of the 3'UTR from a trans-sialidase gene, whose expression is higher in trypomastigotes, resulted in increased luciferase activity in trypomastigotes without a corresponding increase in luciferase mRNA levels.

Trans-splicing

Trans-splicing is the joining together of portions of two separate pre-mRNA molecules. The two distinct categories of spliceosomal trans-splicing are genic trans-splicing, which joins exons of different pre-mRNA transcripts, and spliced leader (SL) trans-splicing, which involves an exon donated from a specialized SL RNA. Both depend primarily on the same signals and components as cis-splicing. Genic trans-splicing events producing protein-coding mRNAs have been described in a variety of organisms, including Caenorhabditis elegans and Drosophila. In mammalian cells, genic trans-splicing can be associated with cancers and translocations. SL trans-splicing has mainly been studied in nematodes and trypanosomes, but there are now numerous and diverse phyla (including primitive chordates) where this type of trans-splicing has been detected. Such diversity raises questions as to the evolutionary origin of the process. Another intriguing question concerns the function of trans-splicing, as operon resolution can only account for a small proportion of the total amount of SL trans-splicing.

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.

Functional genome annotation by combined analysis across microarray studies of Trypanosoma brucei

BACKGROUND

Functional annotation of trypanosomatid genomes has been a daunting task due to the low similarity of their genes with annotated genes of other organisms. Three recent studies have provided gene expression profiles in several different conditions and life stages for one of the main disease-causing trypanosomatids, Trypanosoma brucei. These data can be used to study the gene functions and regulatory mechanisms in this organism.

METHODOLOGY/PRINCIPAL FINDINGS

Combining the data from three different microarray studies of T. brucei, we show that functional linkages among T. brucei genes can be identified based on gene coexpression, leading to a powerful approach for gene function prediction. These predictions can be further improved by considering the expression profiles of orthologous genes from other trypanosomatids. Furthermore, gene expression profiles can be used to discover potential regulatory elements within 3' untranslated regions.

CONCLUSIONS/SIGNIFICANCE

These results suggest that although trypanosomatids do not regulate genes at transcription level, trypanosomatid genes with related functions are coregulated post-transcriptionally via modulation of mRNA stability, implying the presence of complex regulatory networks in these organisms. Our analysis highlights the demand for a thorough transcript profiling of T. brucei genome in parallel with other trypanosomatid genomes, which can provide a powerful means to improve their functional annotation.

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 genome sequence of Trypanosoma brucei gambiense, causative agent of chronic human african trypanosomiasis

Background

Trypanosoma brucei gambiense is the causative agent of chronic Human African Trypanosomiasis or sleeping sickness, a disease endemic across often poor and rural areas of Western and Central Africa. We have previously published the genome sequence of a T. b. brucei isolate, and have now employed a comparative genomics approach to understand the scale of genomic variation between T. b. gambiense and the reference genome. We sought to identify features that were uniquely associated with T. b. gambiense and its ability to infect humans.

Methods and Findings

An improved high-quality draft genome sequence for the group 1 T. b. gambiense DAL 972 isolate was produced using a whole-genome shotgun strategy. Comparison with T. b. brucei showed that sequence identity averages 99.2% in coding regions, and gene order is largely collinear. However, variation associated with segmental duplications and tandem gene arrays suggests some reduction of functional repertoire in T. b. gambiense DAL 972. A comparison of the variant surface glycoproteins (VSG) in T. b. brucei with all T. b. gambiense sequence reads showed that the essential structural repertoire of VSG domains is conserved across T. brucei.

Conclusions

This study provides the first estimate of intraspecific genomic variation within T. brucei, and so has important consequences for future population genomics studies. We have shown that the T. b. gambiense genome corresponds closely with the reference, which should therefore be an effective scaffold for any T. brucei genome sequence data. As VSG repertoire is also well conserved, it may be feasible to describe the total diversity of variant antigens. While we describe several as yet uncharacterized gene families with predicted cell surface roles that were expanded in number in T. b. brucei, no T. b. gambiense-specific gene was identified outside of the subtelomeres that could explain the ability to infect humans.

Sleeping sickness, or Human African Trypanosomiasis, is a disease affecting the health and productivity of poor people in many rural areas of sub-Saharan Africa. The disease is caused by a single-celled flagellate, Trypanosoma brucei, which evades the immune system by periodically switching the proteins on its surface. We have produced a genome sequence for T. brucei gambiense, which is the particular subspecies causing most disease in humans. We compared this with an existing reference genome for a non-human infecting strain (T. b. brucei 927) to identify genes in T. b. gambiense that might explain its ability to infect humans and to assess how well the reference performs as a universal plan for all T. brucei. The genome sequences differ only due to rare insertions and duplications and homologous genes are over 95% identical on average. The archive of surface antigens that enable the parasite to switch its protein coat is remarkably consistent, even though it evolves very quickly. We identified genes with predicted cell surface functions that are only present in T. b. brucei and have evolved rapidly in recent time. These genes might help to explain variation in disease pathology between different T. brucei strains in different hosts.

Digital gene expression analysis of two life cycle stages of the human-infective parasite, Trypanosoma brucei gambiense reveals differentially expressed clusters of co-regulated genes

Background

The evolutionarily ancient parasite, Trypanosoma brucei, is unusual in that the majority of its genes are regulated post-transcriptionally, leading to the suggestion that transcript abundance of most genes does not vary significantly between different life cycle stages despite the fact that the parasite undergoes substantial cellular remodelling and metabolic changes throughout its complex life cycle. To investigate this in the clinically relevant sub-species, Trypanosoma brucei gambiense, which is the causative agent of the fatal human disease African sleeping sickness, we have compared the transcriptome of two different life cycle stages, the potentially human-infective bloodstream forms with the non-human-infective procyclic stage using digital gene expression (DGE) analysis.

Results

Over eleven million unique tags were generated, producing expression data for 7360 genes, covering 81% of the genes in the genome. Compared to microarray analysis of the related T. b. brucei parasite, approximately 10 times more genes with a 2.5-fold change in expression levels were detected. The transcriptome analysis revealed the existence of several differentially expressed gene clusters within the genome, indicating that contiguous genes, presumably from the same polycistronic unit, are co-regulated either at the level of transcription or transcript stability.

Conclusions

DGE analysis is extremely sensitive for detecting gene expression differences, revealing firstly that a far greater number of genes are stage-regulated than had previously been identified and secondly and more importantly, this analysis has revealed the existence of several differentially expressed clusters of genes present on what appears to be the same polycistronic units, a phenomenon which had not previously been observed in microarray studies. These differentially regulated clusters of genes are in addition to the previously identified RNA polymerase I polycistronic units of variant surface glycoproteins and procyclin expression sites, which encode the major surface proteins of the parasite. This raises a number of questions regarding the function and regulation of the gene clusters that clearly warrant further study.

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.

The evolution of amastin surface glycoproteins in trypanosomatid parasites

Amastin is a transmembrane glycoprotein found on the cell surfaces of trypanosomatid parasites. Encoded by a large, diverse gene family, amastin was initially described from the intracellular, amastigote stage of Trypanosoma cruzi and Leishmania donovani. Genome sequences have subsequently shown that the amastin repertoire is much larger in Leishmania relative to Trypanosoma. However, it is not known when this expansion occurred, whether it is associated with the origins of Leishmania and vertebrate parasitism itself, or prior to this. To examine the timing of amastin diversification, as well as the evolutionary mechanisms regulating gene repertoire and sequence diversity, this study sequenced the genomic regions containing amastin loci from two related insect parasites (Leptomonas seymouri and Crithidia sp.) and estimated a phylogeny for these and other amastin sequences. The phylogeny shows that amastin includes four subfamilies with distinct genomic positions, secondary structures, and evolution, which were already differentiated in the ancestral trypanosomatid. Diversification in Leishmania was initiated from a single ancestral locus on chromosome 34, with rapid derivation of novel loci through transposition and accelerated sequence divergence. This is absent from related organisms showing that diversification occurred after the origin of Leishmania. These results describe a substantial elaboration of amastin repertoire directly associated with the origin of Leishmania, suggesting that some amastin genes evolved novel functions crucial to cell function in leishmanial parasites after the acquisition of a vertebrate host.

Alanine aminotransferase of Trypanosoma brucei--a key role in proline metabolism in procyclic life forms

African trypanosomes possess high levels of alanine aminotransferase (EC 2.6.1.2), although the function of their activity remains enigmatic, especially in slender bloodstream forms where the metabolism of ketoacids does not occur. Therefore, the gene for alanine aminotransferase enzyme in Trypanosoma brucei (TbAAT) was characterized and its function assessed using a combination of RNA interference and gene knockout approaches. Surprisingly, as much as 95% or more of the activity appears to be unnecessary for growth of either bloodstream or procyclic forms respiring on glucose. A combination of RNA interference and NMR spectroscopy revealed an important role for the activity in procyclic forms respiring on proline. Under these conditions, the major end product of proline metabolism is alanine, and a reduction in TbAAT activity led to a proportionate decrease in the amount of alanine excreted along with an increase in the doubling time of the cells. These results provide evidence of a role for alanine aminotransferase in the metabolism of proline in African trypanosomes by linking glutamate produced by the initial oxidative steps of the pathway with pyruvate produced by the final oxidative step of the pathway. This step appears to be essential when proline is the primary carbon source, which is likely to be the physiological situation in the tsetse fly vector.

Transcriptome analysis of differentiating trypanosomes reveals the existence of multiple post-transcriptional regulons

Background

Trypanosome gene expression is regulated almost exclusively at the post-transcriptional level, with mRNA degradation playing a decisive role. When trypanosomes are transferred from the blood of a mammal to the midgut of a Tsetse fly, they transform to procyclic forms: gene expression is reprogrammed, changing the cell surface and switching the mode of energy metabolism. Within the blood, trypanosomes can pre-adapt for Tsetse transmission, becoming growth-arrested stumpy forms. We describe here the transitions in gene expression that occur during differentiation of in-vitro cultured bloodstream forms to procyclic forms.

Results

Some mRNAs showed changes within 30 min of cis-aconitate addition, whereas others responded 12-24 hours later. For the first 12 h after addition of cis-aconitate, cells accumulated at the G1 phase of the cell cycle, and showed decreases in mRNAs required for proliferation, mimicking the changes seen in stumpy forms: many mRNAs needed for ribosomal and flagellar biogenesis showed striking co-regulation. Other mRNAs encoding components of signal transduction pathways and potential regulators were specifically induced only during differentiation. Messenger RNAs encoding proteins required for individual metabolic pathways were often co-regulated.

Conclusion

Trypanosome genes form post-transcriptional regulons in which mRNAs with functions in particular pathways, or encoding components of protein complexes, show almost identical patterns of regulation.

Trypanosoma brucei PUF9 regulates mRNAs for proteins involved in replicative processes over the cell cycle

Many genes that are required at specific points in the cell cycle exhibit cell cycle–dependent expression. In the early-diverging model eukaryote and important human pathogen Trypanosoma brucei, regulation of gene expression in the cell cycle and other processes is almost entirely post-transcriptional. Here, we show that the T. brucei RNA-binding protein PUF9 stabilizes certain transcripts during S-phase. Target transcripts of PUF9—LIGKA, PNT1 and PNT2—were identified by affinity purification with TAP-tagged PUF9. RNAi against PUF9 caused an accumulation of cells in G2/M phase and unexpectedly destabilized the PUF9 target mRNAs, despite the fact that most known Puf-domain proteins promote degradation of their target mRNAs. The levels of the PUF9-regulated transcripts were cell cycle dependent, peaking in mid- to late- S-phase, and this effect was abolished when PUF9 was targeted by RNAi. The sequence UUGUACC was over-represented in the 3′ UTRs of PUF9 targets; a point mutation in this motif abolished PUF9-dependent stabilization of a reporter transcript carrying the PNT1 3′ UTR. LIGKA is involved in replication of the kinetoplast, and here we show that PNT1 is also kinetoplast-associated and its over-expression causes kinetoplast-related defects, while PNT2 is localized to the nucleus in G1 phase and redistributes to the mitotic spindle during mitosis. PUF9 targets may constitute a post-transcriptional regulon, encoding proteins involved in temporally coordinated replicative processes in early G2 phase.

The unicellular protozoan Trypanosoma brucei is the causative agent of African sleeping sickness, responsible for over 100,000 deaths annually, and is related to other important pathogens (e.g. Leishmania major and Trypanosoma cruzi). Unusually, these organisms do not regulate their genes by changing the rate at which they are copied into RNA, but by changing the rate of RNA destruction or the rate of translation into protein. We identified an RNA-binding protein, PUF9, responsible for the accumulation of several RNA molecules at a specific time point in the cell division cycle, just after DNA replication. Correspondingly, the proteins encoded by these RNAs appear to function in the division of various cellular structures at this time point or shortly afterwards. Two of them facilitate replication of the kinetoplast (an organelle containing the mitochondrial DNA) while another was found in the mitotic spindle. Their temporal co-expression may stem from another unusual feature of trypanosomes: only one copy of the kinetoplast (and several other organelles) are present per cell, their replication being coordinated with cell division. Indeed, PUF9 may be important in the control of organelle copy-number because suppression of PUF9 resulted in cells with too many kinetoplasts, flagella, or nuclei.

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.

Analysis of expressed sequence tags from the four main developmental stages of Trypanosoma congolense

Trypanosoma congolense is one of the most economically important pathogens of livestock in Africa. Culture-derived parasites of each of the three main insect stages of the T. congolense life cycle, i.e., the procyclic, epimastigote and metacyclic stages, and bloodstream stage parasites isolated from infected mice, were used to construct stage-specific cDNA libraries and expressed sequence tags (ESTs or cDNA clones) in each library were sequenced. Thirteen EST clusters encoding different variant surface glycoproteins (VSGs) were detected in the metacyclic library and 26 VSG EST clusters were found in the bloodstream library, 6 of which are shared by the metacyclic library. Rare VSG ESTs are present in the epimastigote library, and none were detected in the procyclic library. ESTs encoding enzymes that catalyze oxidative phosphorylation and amino acid metabolism are about twice as abundant in the procyclic and epimastigote stages as in the metacyclic and bloodstream stages. In contrast, ESTs encoding enzymes involved in glycolysis, the citric acid cycle and nucleotide metabolism are about the same in all four developmental stages. Cysteine proteases, kinases and phosphatases are the most abundant enzyme groups represented by the ESTs. All four libraries contain T. congolense-specific expressed sequences not present in the Trypanosoma brucei and Trypanosoma cruzi genomes. Normalized cDNA libraries were constructed from the metacyclic and bloodstream stages, and found to be further enriched for T. congolense-specific ESTs. Given that cultured T. congolense offers an experimental advantage over other African trypanosome species, these ESTs provide a basis for further investigation of the molecular properties of these four developmental stages, especially the epimastigote and metacyclic stages for which it is difficult to obtain large quantities of organisms. The T. congolense EST databases are available at: http://www.sanger.ac.uk/Projects/T_congolense/EST_index.shtml. The sequence data have been submitted to EMBL under the following accession numbers: FN263376-FN292969.

Bidirectional silencing of RNA polymerase I transcription by a strand switch region in Trypanosoma brucei

The procyclin genes in Trypanosoma brucei are transcribed by RNA polymerase I as part of 5–10 kb long polycistronic transcription units on chromosomes VI and X. Each procyclin locus begins with two procyclin genes followed by at least one procyclin-associated gene (PAG). In procyclic (insect midgut) form trypanosomes, PAG mRNA levels are about 100-fold lower than those of procyclins. We show that deletion of PAG1, PAG2 or PAG3 results in increased mRNA levels from downstream genes in the same transcription unit. Nascent RNA analysis revealed that most of the effects are due to increased transcription elongation in the knockouts. Furthermore, transient and stable transfections showed that sequence elements on both strands of PAG1 can inhibit Pol I transcription. Finally, by database mining we identified 30 additional PAG-related sequences that are located almost exclusively at strand switch regions and/or at sites where a change of RNA polymerase type is likely to occur.

Organization and evolution of two SIDER retroposon subfamilies and their impact on the Leishmania genome

Background

We have recently identified two large families of extinct transposable elements termed Short Interspersed DEgenerated Retroposons (SIDERs) in the parasitic protozoan Leishmania major. The characterization of SIDER elements was limited to the SIDER2 subfamily, although members of both subfamilies have been shown to play a role in the regulation of gene expression at the post-transcriptional level. Apparent functional domestication of SIDERs prompted further investigation of their characterization, dissemination and evolution throughout the Leishmania genus, with particular attention to the disregarded SIDER1 subfamily.

Results

Using optimized statistical profiles of both SIDER1 and SIDER2 subgroups, we report the first automated and highly sensitive annotation of SIDERs in the genomes of L. infantum, L. braziliensis and L. major. SIDER annotations were combined to in-silico mRNA extremity predictions to generate a detailed distribution map of the repeat family, hence uncovering an enrichment of antisense-oriented SIDER repeats between the polyadenylation and trans-splicing sites of intergenic regions, in contrast to the exclusive sense orientation of SIDER elements within 3'UTRs. Our data indicate that SIDER elements are quite uniformly dispersed throughout all three genomes and that their distribution is generally syntenic. However, only 47.4% of orthologous genes harbor a SIDER element in all three species. There is evidence for species-specific enrichment of SIDERs and for their preferential association, especially for SIDER2s, with different metabolic functions. Investigation of the sequence attributes and evolutionary relationship of SIDERs to other trypanosomatid retroposons reveals that SIDER1 is a truncated version of extinct autonomous ingi-like retroposons (DIREs), which were functional in the ancestral Leishmania genome.

Conclusion

A detailed characterization of the sequence traits for both SIDER subfamilies unveils major differences. The SIDER1 subfamily is more heterogeneous and shows an evolutionary link with vestigial DIRE retroposons as previously observed for the ingi/RIME and L1Tc/NARTc couples identified in the T. brucei and T. cruzi genomes, whereas no identified DIREs are related to SIDER2 sequences. Although SIDER1s and SIDER2s display equivalent genomic distribution globally, the varying degrees of sequence conservation, preferential genomic disposition, and differential association to orthologous genes allude to an intricate web of SIDER assimilation in these parasitic organisms.

Functional characterization and protein-protein interactions of trypanosome splicing factors U2AF35, U2AF65 and SF1

Early in the assembly of eukaryotes the branch-point binding protein (BBP, also called SF1) recognizes the branch point sequence, whereas the heterodimer U2AF, consisting of a 65 and a 35 kDa subunit, contacts the polypyrimidine tract and the AG splice site, respectively. Herein, we identified, cloned and expressed the Trypanosoma cruzi and Trypanosoma brucei U2AF35, U2AF65 and SF1. Trypanosomatid U2AF65 strongly diverged from yeast and human homologues. On the contrary, trypanosomatid SF1 was conserved but lacked the C-terminal sequence present in the mammalian protein. Yeast two hybrid approaches were used to assess their interactions. The interaction between U2AF35 and U2AF65 was very weak or not detectable. However, as in other eukaryotes, the interaction between U2AF65 and SF1 was strong. At the cellular level, these results were confirmed by fractionation and affinity-selection experiments in which SF1 and U2AF65 were affinity-selected with TAP tagged SF1, but not with TAP tagged U2AF35. Silencing one of the three factors affected growth and trans-splicing in the first step of this reaction. Trypanosomes are the first described example of eukaryotic cells in which the interaction of two expressed U2AF factors seemed to be very weak, or not detectable.

Multiple roles for polypyrimidine tract binding (PTB) proteins in trypanosome RNA metabolism

Trypanosomatid genomes encode for numerous proteins containing an RNA recognition motif (RRM), but the function of most of these proteins in mRNA metabolism is currently unknown. Here, we report the function of two such proteins that we have named PTB1 and PTB2, which resemble the mammalian polypyrimidine tract binding proteins (PTB). RNAi silencing of these factors indicates that both are essential for life. PTB1 and PTB2 reside mostly in the nucleus, but are found in the cytoplasm, as well. Microarray analysis performed on PTB1 and PTB2 RNAi silenced cells indicates that each of these factors differentially affects the transcriptome, thus regulating a different subset of mRNAs. PTB1 and PTB2 substrates were categorized bioinformatically, based on the presence of PTB binding sites in their 5' and 3' flanking sequences. Both proteins were shown to regulate mRNA stability. Interestingly, PTB proteins are essential for trans-splicing of genes containing C-rich polypyrimidine tracts. PTB1, but not PTB2, also affects cis-splicing. The specificity of binding of PTB1 was established in vivo and in vitro using a model substrate. This study demonstrates for the first time that trans-splicing of only certain substrates requires specific factors such as PTB proteins for their splicing. The trypanosome PTB proteins, like their mammalian homologs, represent multivalent RNA binding proteins that regulate mRNAs from their synthesis to degradation.

The FIP-1 like polyadenylation factor in trypanosomes and the structural basis for its interaction with CPSF30

In trypanosomes transcription is polycistronic and individual mRNAs are generated by a trans-splicing/polyadenylation coupled reaction. We identified a divergent trypanosome FIP1-like, a factor required for mRNA 3' end formation from yeasts to human. Here we showed that it is a nuclear protein with a speckled distribution essential for trypanosome viability. A strong interaction was found between TcFIP1-like and TcCPSF30, a component of the polyadenylation complex. We determined the specific amino acids in each protein involved in the interaction. Significant differences were found between the trypanosome interaction surface and its human counterpart. Although CPSF30/FIP1 interaction is known in other organisms, this is the first report mapping the interaction surface at the amino acid level.

Functionally related transcripts have common RNA motifs for specific RNA-binding proteins in trypanosomes

Background

Trypanosomes mostly control gene expression by post-transcriptional events such as modulation of mRNA stability and translational efficiency. These mechanisms involve RNA-binding proteins (RBPs), which associate with transcripts to form messenger ribonucleoprotein (mRNP) complexes.

Results

In this study, we report the identification of mRNA targets for Trypanosoma cruzi U-rich RBP 1 (TcUBP1) and T. cruzi RBP 3 (TcRBP3), two phylogenetically conserved proteins among Kinetoplastids. Co-immunoprecipitated RBP-associated RNAs were extracted from mRNP complexes and binding of RBPs to several targets was confirmed by independent experimental assays. Analysis of target transcript sequences allowed the identification of different signature RNA motifs for each protein. Cis-elements for RBP binding have a stem-loop structure of 30–35 bases and are more frequently represented in the 3'-untranslated region (UTR) of mRNAs. Insertion of the correctly folded RNA elements to a non-specific mRNA rendered it into a target transcript, whereas substitution of the RNA elements abolished RBP interaction. In addition, RBPs competed for RNA-binding sites in accordance with the distribution of different and overlapping motifs in the 3'-UTRs of common mRNAs.

Conclusion

Functionally related transcripts were preferentially associated with a given RBP; TcUBP1 targets were enriched in genes encoding proteins involved in metabolism, whereas ribosomal protein-encoding transcripts were the largest group within TcRBP3 targets. Together, these results suggest coordinated control of different mRNA subsets at the post-transcriptional level by specific RBPs.

Identification and characterization of nuclear non-canonical poly(A) polymerases from Trypanosoma brucei

Regulation of nuclear genome expression in Trypanosoma brucei is critical for this protozoan parasite's successful transition between its vertebrate and invertebrate host environments. The canonical eukaryotic circuits such as modulation of transcription initiation, mRNA splicing and polyadenylation appear to be nearly non-existent in T. brucei suggesting that the transcriptome is primarily defined by mRNA turnover. Our previous work has highlighted sequence similarities between terminal RNA uridylyl transferases (TUTases) and non-canonical poly(A) polymerases, which are widely implicated in regulating nuclear, cytoplasmic and organellar RNA decay throughout the eukaryotic lineage. Here, we have continued characterization of TUTase-like proteins in T. brucei and identified two nuclear non-canonical poly(A) polymerases (ncPAPs). The 82kDa TbncPAP1 is essential for viability of procyclic and bloodstream forms of T. brucei. Similar to Trf4/5 proteins from budding yeast, TbncPAP1 requires protein cofactor(s) to exert poly(A) polymerase activity in vitro. The recombinant 54kDa TbncPAP2 showed a PAP activity as an individual polypeptide. Proteomic analysis of the TbncPAP1 interactions demonstrated its association with Mtr4 RNA helicase and several RNA binding proteins, including a potential ortholog of Air1p/2p proteins, which indicates the presence of a stable TRAMP-like complex in trypanosomes. Our findings suggest that TbncPAP1 may be a "quality control" nuclear PAP involved in targeting aberrant or anti-sense transcripts for degradation by the 3'-exosome. Such mechanisms are likely to play a major role in alleviating promiscuity of the transcriptional machinery.

Differential expression of a protease gene family in African trypanosomes

During their life cycle African trypanosomes must quickly adapt to the different environments of the tsetse fly midgut and the mammalian bloodstream by modulating expression of many of their genes. One group of these differentially expressed genes encodes different forms of a major surface protease. Using a luciferase reporter gene transiently or permanently transfected into trypanosomes, we show here that the 3'-UTRs of these protease genes are responsible for their differential expression. Deletion analysis of the 389-bp 3'-UTR of one of the protease genes, MSP-B, demonstrated that it contains a U-rich regulatory region of about 23bp (UCGUCUGUUAUUUCUUAGUCCAG), which suppresses expression of the reporter protein in bloodstream trypanosomes by as much as 25-fold, but has little effect on the reporter expression in procyclic (tsetse fly) trypanosomes. Replacing the entire 3'-UTR with just this 23-bp element mimicked most of the suppression effect of the complete 3'-UTR. Northern blots showed that the 23-bp element influences the steady state RNA level, but not enough to account for the 25-fold suppression effect. Polysome analyses showed that in procyclic trypanosomes more of the total protease mRNA is associated with intermediate-sized and large polysomes than in bloodstream trypanosomes. Thus, the 23-bp element of this protease gene affects both the level of RNA and its translation.

Sequences involved in mRNA processing in Trypanosoma cruzi

Gene expression in Trypanosomatids requires processing of polycistronic transcripts to generate monocistronic mRNAs by cleavage events that are coupled to the addition of a Spliced Leader sequence (SL) at the 5'-end and a poly(A) tail at the 3'-end of each mRNA. Here we investigate the sequence requirements involved in Trypanosoma cruzi mRNA processing by mapping all available expressed sequence tags and cDNAs containing poly(A) tail and/or SL to genomic intergenic regions. Amongst other parameters, we determined that the median lengths of 5' untranslated region (UTR) and 3'UTR sequences are 35 and 264 nucleotides, respectively; and that the median distance between SL addition sites and a polypyrimidine motif is 18 nucleotides, whereas the median distance between poly(A) addition sites and the closest polypyrimidine-rich sequence is 40 nucleotides.

A role for Caf1 in mRNA deadenylation and decay in trypanosomes and human cells

The eukaryotic Ccr4/Caf1/Not complex is involved in deadenylation of mRNAs. The Caf1 and Ccr4 subunits both potentially have deadenylating enzyme activity. We investigate here the roles of Ccr4 and Caf1 in deadenylation in two organisms that separated early in eukaryotic evolution: humans and trypanosomes. In Trypanosoma brucei, we found a complex containing CAF1, NOT1, NOT2 and NOT5, DHH1 and a possible homologue of Caf130; no homologue of Ccr4 was found. Trypanosome CAF1 has deadenylation activity, and is essential for cell survival. Depletion of trypanosome CAF1 delayed deadenylation and degradation of constitutively expressed mRNAs. Human cells have two isozymes of Caf1: simultaneous depletion of both inhibited degradation of an unstable reporter mRNA. In both species, depletion of Caf1 homologues inhibited deadenylation of bulk RNA and resulted in an increase in average poly(A) tail length.

Transgenic, fluorescent Leishmania mexicana allow direct analysis of the proteome of intracellular amastigotes

Investigating the proteome of intracellular pathogens is often hampered by inadequate methodologies to purify the pathogen free of host cell material. This has also precluded direct proteome analysis of the intracellular, amastigote form of Leishmania spp., protozoan parasites that cause a spectrum of diseases that affect some 12 million patients worldwide. Here a method is presented that combines classic, isopycnic density centrifugation with fluorescent particle sorting for purification by exploiting transgenic, fluorescent parasites to allow direct proteome analysis of the purified organisms. By this approach the proteome of intracellular Leishmania mexicana amastigotes was compared with that of extracellular promastigotes that are transmitted by insect vectors. In total, 509 different proteins were identified by mass spectrometry and database search. This number corresponds to approximately 6% of gene products predicted from the reference genome of Leishmania major. Intracellular amastigotes synthesized significantly more proteins with basic pI and showed a greater abundance of enzymes of fatty acid catabolism, which may reflect their living in acidic habitats and metabolic adaptation to nutrient availability, respectively. Bioinformatics analyses of the genes corresponding to the protein data sets produced clear evidence for skewed codon usage and translational bias in these organisms. Moreover analysis of the subset of genes whose products were more abundant in amastigotes revealed characteristic sequence motifs in 3'-untranslated regions that have been linked to translational control elements. This suggests that proteome data sets may be used to identify regulatory elements in mRNAs. Last but not least, at 6% coverage the proteome identified all vaccine antigens tested to date. Thus, the present data set provides a valuable resource for selection of candidate vaccine antigens.

Improving the prediction of mRNA extremities in the parasitic protozoan Leishmania

Background

Leishmania and other members of the Trypanosomatidae family diverged early on in eukaryotic evolution and consequently display unique cellular properties. Their apparent lack of transcriptional regulation is compensated by complex post-transcriptional control mechanisms, including the processing of polycistronic transcripts by means of coupled trans-splicing and polyadenylation. Trans-splicing signals are often U-rich polypyrimidine (poly(Y)) tracts, which precede AG splice acceptor sites. However, as opposed to higher eukaryotes there is no consensus polyadenylation signal in trypanosomatid mRNAs.

Results

We refined a previously reported method to target 5' splice junctions by incorporating the pyrimidine content of query sequences into a scoring function. We also investigated a novel approach for predicting polyadenylation (poly(A)) sites in-silico, by comparing query sequences to polyadenylated expressed sequence tags (ESTs) using position-specific scanning matrices (PSSMs). An additional analysis of the distribution of putative splice junction to poly(A) distances helped to increase prediction rates by limiting the scanning range. These methods were able to simplify splice junction prediction without loss of precision and to increase polyadenylation site prediction from 22% to 47% within 100 nucleotides.

Conclusion

We propose a simplified trans-splicing prediction tool and a novel poly(A) prediction tool based on comparative sequence analysis. We discuss the impact of certain regions surrounding the poly(A) sites on prediction rates and contemplate correlating biological mechanisms. This work aims to sharpen the identification of potentially functional untranslated regions (UTRs) in a large-scale, comparative genomics framework.