A common pyrimidine-rich motif governs trans-splicing and polyadenylation of tubulin polycistronic pre-mRNA in trypanosomes

Life stage-specific poly(A) site selection regulated by Trypanosoma brucei DRBD18

The kinetoplastid parasite, Trypanosoma brucei, undergoes a complex life cycle entailing slender and stumpy bloodstream forms in mammals and procyclic and metacyclic forms (MFs) in tsetse fly hosts. The numerous gene regulatory events that underlie T. brucei differentiation between hosts, as well as between active and quiescent stages within each host, take place in the near absence of transcriptional control. Rather, differentiation is controlled by RNA-binding proteins (RBPs) that associate with mRNA 3' untranslated regions (3'UTRs) to impact RNA stability and translational efficiency. DRBD18 is a multifunctional T. brucei RBP, shown to impact mRNA stability, translation, export, and processing. Here, we use single-cell RNAseq to characterize transcriptomic changes in cell populations that arise upon DRBD18 depletion, as well as to visualize transcriptome-wide alterations to 3'UTR length. We show that in procyclic insect stages, DRBD18 represses expression of stumpy bloodstream form and MF transcripts. Additionally, DRBD18 regulates the 3'UTR lengths of over 1,500 transcripts, typically promoting the use of distal polyadenylation sites, and thus the inclusion of 3'UTR regulatory elements. Remarkably, comparison of polyadenylation patterns in DRBD18 knockdowns with polyadenylation patterns in stumpy bloodstream forms shows numerous similarities, revealing a role for poly(A) site selection in developmental gene regulation, and indicating that DRBD18 controls this process for a set of transcripts. RNA immunoprecipitation supports a direct role for DRBD18 in poly(A) site selection. This report highlights the importance of alternative polyadenylation in T. brucei developmental control and identifies a critical RBP in this process.

A versatile 2A peptide-based strategy for ectopic expression and endogenous gene tagging in Trypanosoma cruzi

Nearly all expression vectors currently available for Trypanosoma cruzi were conceived to produce a single primary transcript containing the genes of interest along with those that confer antibiotic resistance. However, since each messenger RNA (mRNA) matures separately, drug selection will only guarantee the expression of those derived from the selectable marker. Therefore, commonly a considerable fraction of the cells recovered after selection with these expression vectors, although resistant do not express the protein of interest. Consequently, in order to counteract this disadvantage, we developed vectors with an alternative arrangement in which the gene of interest and antibiotic resistance are fused sharing the same mRNA. To test this configuration, we included the coding sequence for the green fluorescent protein (mEGFP) linked to the one conferring neomycin resistance (Neo). Additionally, to allow for the production of two independent proteins the sequence for a Thosea asigna virus self-cleaving 2A peptide (T2A) was inserted in-between. Cells obtained with these vectors displayed higher mEGFP expression levels with more homogeneous transgenic parasite populations than those transfected with more conventional independent mRNA-based alternatives. Moreover, as determined by Western blot, 2A mediated fusion protein dissociation occurred with high efficiency in all parasite stages. In addition, these vectors could easily be transformed into endogenous tagging constructs that allowed the insertion, by ends-in homologous recombination, of a hemagglutinin tag (HA) fused to the actin gene. The use of 2A self-cleaving peptides in the context of single mRNA vectors represents an interesting strategy capable of improving ectopic transgene expression in T. cruzi as well as providing a simple alternative to more sophisticated methods, such as the one based on CRISPR/Cas9, for the endogenous labeling of genes.

Genome-wide chromatin interaction map for Trypanosoma cruzi

Trypanosomes are eukaryotic, unicellular parasites, such as Trypanosoma brucei, which causes sleeping sickness, and Trypanosoma cruzi, which causes Chagas disease. Genomes of these parasites comprise core regions and species-specific disruptive regions that encode multigene families of surface glycoproteins. Few transcriptional regulators have been identified in these parasites, and the role of spatial organization of the genome in gene expression is unclear. Here we mapped genome-wide chromatin interactions in T. cruzi using chromosome conformation capture (Hi-C), and we show that the core and disruptive regions form three-dimensional chromatin compartments named C and D. These chromatin compartments differ in levels of DNA methylation, nucleosome positioning and chromatin interactions, affecting genome expression dynamics. Our data reveal that the trypanosome genome is organized into chromatin-folding domains and transcription is affected by the local chromatin structure. We propose a model in which epigenetic mechanisms affect gene expression in trypanosomes.

Cell-cell communication in African trypanosomes

Years of research have shown us that unicellular organisms do not exist entirely in isolation, but rather that they are capable of an altogether far more sociable way of living. Single cells produce, receive and interpret signals, coordinating and changing their behaviour according to the information received. Although this cell-cell communication has long been considered the norm in the bacterial world, an increasing body of knowledge is demonstrating that single-celled eukaryotic parasites also maintain active social lives. This communication can drive parasite development, facilitate the invasion of new niches and, ultimately, influence infection outcome. In this review, I present the evidence for cell-cell communication during the life cycle of the African trypanosomes, from their mammalian hosts to their insect vectors, and reflect on the many remaining unanswered questions in this fascinating field.

The paradigm of intracellular parasite survival and drug resistance in leishmanial parasite through genome plasticity and epigenetics: Perception and future perspective

Leishmania is an intracellular, zoonotic, kinetoplastid eukaryote with more than 1.2 million cases all over the world. The leishmanial chromosomes are divided into polymorphic chromosomal ends, conserved central domains, and antigen-encoding genes found in telomere-proximal regions. The genome flexibility of chromosomal ends of the leishmanial parasite is known to cause drug resistance and intracellular survival through the evasion of host defense mechanisms. Therefore, in this review, we discuss the plasticity of Leishmania genome organization which is the primary cause of drug resistance and parasite survival. Moreover, we have not only elucidated the causes of such genome plasticity which includes aneuploidy, epigenetic factors, copy number variation (CNV), and post-translation modification (PTM) but also highlighted their impact on drug resistance and parasite survival.

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.

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

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

Identification of sequence-specific promoters driving polycistronic transcription initiation by RNA polymerase II in trypanosomes

Protein-coding genes in trypanosomes occur in polycistronic transcription units (PTUs). How RNA polymerase II (Pol II) initiates transcription of PTUs has not been resolved; the current model favors chromatin modifications inducing transcription rather than sequence-specific promoters. Here, we uncover core promoters by functional characterization of Pol II peaks identified by chromatin immunoprecipitation sequencing (ChIP-seq). Two distinct promoters are located between divergent PTUs, each driving unidirectional transcription. Detailed analysis identifies a 75-bp promoter that is necessary and sufficient to drive full reporter expression and contains functional motifs. Analysis of further promoters suggests transcription initiation is regulated and promoters are either focused or dispersed. In contrast to the previous model of unregulated and promoter-independent transcription initiation, we find that sequence-specific promoters determine the initiation of Pol II transcription of protein-coding genes PTUs. These findings in Trypanosoma brucei suggest that in addition of chromatin modifications, promoter motifs-based regulation of gene expression is deeply conserved among eukaryotes.

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.

Characterization of an Atypical eIF4E Ortholog in Leishmania, LeishIF4E-6

Leishmania parasites are digenetic protists that shuffle between sand fly vectors and mammalian hosts, transforming from flagellated extracellular promastigotes that reside within the intestinal tract of female sand flies to the obligatory intracellular and non-motile amastigotes within mammalian macrophages. Stage differentiation is regulated mainly by post-transcriptional mechanisms, including translation regulation. Leishmania parasites encode six different cap-binding proteins, LeishIF4E1-6, that show poor conservation with their counterparts from higher eukaryotes and among themselves. In view of the changing host milieu encountered throughout their life cycle, we propose that each LeishIF4E has a unique role, although these functions may be difficult to determine. Here we characterize LeishIF4E-6, a unique eIF4E ortholog that does not readily associate with m⁷GTP cap in either of the tested life forms of the parasite. We discuss the potential effect of substituting two essential tryptophan residues in the cap-binding pocket, expected to be involved in the cap-binding activity, as judged from structural studies in the mammalian eIF4E. LeishIF4E-6 binds to LeishIF4G-5, one of the five eIF4G candidates in Leishmania. However, despite this binding, LeishIF4E-6 does not appear to function as a translation factor. Its episomal overexpression causes a general reduction in the global activity of protein synthesis, which was not observed in the hemizygous deletion mutant generated by CRISPR-Cas9. This genetic profile suggests that LeishIF4E-6 has a repressive role. The interactome of LeishIF4E-6 highlights proteins involved in RNA metabolism such as the P-body marker DHH1, PUF1 and an mRNA-decapping enzyme that is homologous to the TbALPH1.

Advances in Understanding Leishmania Pathobiology: What Does RNA-Seq Tell Us?

Leishmaniasis is a vector-borne disease caused by a protozoa parasite from over 20 Leishmania species. The clinical manifestations and the outcome of the disease vary greatly. Global RNA sequencing (RNA-Seq) analyses emerged as a powerful technique to profile the changes in the transcriptome that occur in the Leishmania parasites and their infected host cells as the parasites progresses through their life cycle. Following the bite of a sandfly vector, Leishmania are transmitted to a mammalian host where neutrophils and macrophages are key cells mediating the interactions with the parasites and result in either the elimination the infection or contributing to its proliferation. This review focuses on RNA-Seq based transcriptomics analyses and summarizes the main findings derived from this technology. In doing so, we will highlight caveats in our understanding of the parasite's pathobiology and suggest novel directions for research, including integrating more recent data highlighting the role of the bacterial members of the sandfly gut microbiota and the mammalian host skin microbiota in their potential role in influencing the quantitative and qualitative aspects of leishmaniasis pathology.

Pseudouridines on Trypanosoma brucei mRNAs are developmentally regulated: Implications to mRNA stability and protein binding

The parasite Trypanosoma brucei cycles between an insect and a mammalian host and is the causative agent of sleeping sickness. Here, we performed high-throughput mapping of pseudouridines (Ψs) on mRNA from two life stages of the parasite. The analysis revealed ~273 Ψs, including developmentally regulated Ψs that are guided by homologs of pseudouridine synthases (PUS1, 3, 5, and 7). Mutating the U that undergoes pseudouridylation in the 3' UTR of valyl-tRNA synthetase destabilized the mRNA level. To investigate the mechanism by which Ψ affects the stability of this mRNA, proteins that bind to the 3' UTR were identified, including the RNA binding protein RBSR1. The binding of RBSR1 protein to the 3' UTR was stronger when lacking Ψ compared to transcripts carrying the modification, suggesting that Ψ can inhibit the binding of proteins to their target and thus affect the stability of mRNAs. Consequently, Ψ modification on mRNA adds an additional level of regulation to the dominant post-transcriptional control in these parasites.

Parallel monitoring of RNA abundance, localization and compactness with correlative single molecule FISH on LR White embedded samples

Single mRNA molecules are frequently detected by single molecule fluorescence in situ hybridization (smFISH) using branched DNA technology. While providing strong and background-reduced signals, the method is inefficient in detecting mRNAs within dense structures, in monitoring mRNA compactness and in quantifying abundant mRNAs. To overcome these limitations, we have hybridized slices of high pressure frozen, freeze-substituted and LR White embedded cells (LR White smFISH). mRNA detection is physically restricted to the surface of the resin. This enables single molecule detection of RNAs with accuracy comparable to RNA sequencing, irrespective of their abundance, while at the same time providing spatial information on RNA localization that can be complemented with immunofluorescence and electron microscopy, as well as array tomography. Moreover, LR White embedding restricts the number of available probe pair recognition sites for each mRNA to a small subset. As a consequence, differences in signal intensities between RNA populations reflect differences in RNA structures, and we show that the method can be employed to determine mRNA compactness. We apply the method to answer some outstanding questions related to trans-splicing, RNA granules and mitochondrial RNA editing in single-cellular trypanosomes and we show an example of differential gene expression in the metazoan Caenorhabditis elegans.

Nuclear mRNA maturation and mRNA export control: from trypanosomes to opisthokonts

The passage of mRNAs through the nuclear pores into the cytoplasm is essential in all eukaryotes. For regulation, mRNA export is tightly connected to the full machinery of nuclear mRNA processing, starting at transcription. Export competence of pre-mRNAs gradually increases by both transient and permanent interactions with multiple RNA processing and export factors. mRNA export is best understood in opisthokonts, with limited knowledge in plants and protozoa. Here, I review and compare nuclear mRNA processing and export between opisthokonts and Trypanosoma brucei. The parasite has many unusual features in nuclear mRNA processing, such as polycistronic transcription and trans-splicing. It lacks several nuclear complexes and nuclear-pore-associated proteins that in opisthokonts play major roles in mRNA export. As a consequence, trypanosome mRNA export control is not tight and export can even start co-transcriptionally. Whether trypanosomes regulate mRNA export at all, or whether leakage of immature mRNA to the cytoplasm is kept to a low level by a fast kinetics of mRNA processing remains to be investigated. mRNA export had to be present in the last common ancestor of eukaryotes. Trypanosomes are evolutionary very distant from opisthokonts and a comparison helps understanding the evolution of mRNA export.

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.

Rapid block of pre-mRNA splicing by chemical inhibition of analog-sensitive CRK9 in Trypanosoma brucei

Trypanosoma brucei CRK9 is an essential cyclin-dependent kinase for the parasite-specific mode of pre-mRNA processing. In trypanosomes, protein coding genes are arranged in directional arrays that are transcribed polycistronically, and individual mRNAs are generated by spliced leader trans-splicing and polyadenylation, processes that are functionally linked. Since CRK9 silencing caused a decline of mRNAs, a concomitant increase of unspliced pre-mRNAs and the disappearance of the trans-splicing Y structure intermediate, CRK9 is essential for the first step of splicing. CRK9 depletion also caused a loss of phosphorylation in RPB1, the largest subunit of RNA polymerase (pol) II. Here, we established cell lines that exclusively express analog-sensitive CRK9 (CRK9AS ). Inhibition of CRK9AS in these cells by the ATP-competitive inhibitor 1-NM-PP1 reproduced the splicing defects and proved that it is the CKR9 kinase activity that is required for pre-mRNA processing. Since defective trans-splicing was detected as early as 5 min after inhibitor addition, CRK9 presumably carries out reversible phosphorylation on the pre-mRNA processing machinery. Loss of RPB1 phosphorylation, however, took 12-24 hr. Surprisingly, RNA pol II-mediated RNA synthesis in 24 hr-treated cells was upregulated, indicating that, in contrast to other eukaryotes, RPB1 phosphorylation is not a prerequisite for transcription in trypanosomes.

Unexpected diversity in eukaryotic transcription revealed by the retrotransposon hotspot family of Trypanosoma brucei

The path from DNA to RNA to protein in eukaryotes is guided by a series of factors linking transcription, mRNA export and translation. Many of these are conserved from yeast to humans. Trypanosomatids, which diverged early in the eukaryotic lineage, exhibit unusual features such as polycistronic transcription and trans-splicing of all messenger RNAs. They possess basal transcription factors, but lack recognisable orthologues of many factors required for transcription elongation and mRNA export. We show that retrotransposon hotspot (RHS) proteins fulfil some of these functions and that their depletion globally impairs nascent RNA synthesis by RNA polymerase II. Three sub-families are part of a coordinated process in which RHS6 is most closely associated with chromatin, RHS4 is part of the Pol II complex and RHS2 connects transcription with the translation machinery. In summary, our results show that the components of eukaryotic transcription are far from being universal, and reveal unsuspected plasticity in the course of evolution.

Genome-wide mapping reveals conserved and diverged R-loop activities in the unusual genetic landscape of the African trypanosome genome

R-loops are stable RNA-DNA hybrids that have been implicated in transcription initiation and termination, as well as in telomere maintenance, chromatin formation, and genome replication and instability. RNA Polymerase (Pol) II transcription in the protozoan parasite Trypanosoma brucei is highly unusual: virtually all genes are co-transcribed from multigene transcription units, with mRNAs generated by linked trans-splicing and polyadenylation, and transcription initiation sites display no conserved promoter motifs. Here, we describe the genome-wide distribution of R-loops in wild type mammal-infective T. brucei and in mutants lacking RNase H1, revealing both conserved and diverged functions. Conserved localization was found at centromeres, rRNA genes and retrotransposon-associated genes. RNA Pol II transcription initiation sites also displayed R-loops, suggesting a broadly conserved role despite the lack of promoter conservation or transcription initiation regulation. However, the most abundant sites of R-loop enrichment were within the regions between coding sequences of the multigene transcription units, where the hybrids coincide with sites of polyadenylation and nucleosome-depletion. Thus, instead of functioning in transcription termination the most widespread localization of R-loops in T. brucei suggests a novel correlation with pre-mRNA processing. Finally, we find little evidence for correlation between R-loop localization and mapped sites of DNA replication initiation.

Simplified inducible system for Trypanosoma brucei

Nowadays, most reverse genetics approaches in Trypanosoma brucei, a protozoan parasite of medical and veterinary importance, rely on pre-established cell lines. Consequently, inducible experimentation is reduced to a few laboratory strains. Here we described a new transgene expression system based exclusively on endogenous transcription activities and a minimum set of regulatory components that can easily been adapted to different strains. The pTbFIX vectors are designed to contain the sequence of interest under the control of an inducible rRNA promoter along with a constitutive dicistronic unit encoding a nucleus targeted tetracycline repressor and puromycin resistance genes in a tandem “head-to-tail” configuration. Upon doxycycline induction, the system supports regulatable GFP expression (170 to 400 fold) in both bloodstream and procyclic T. brucei forms. Furthermore we have adapted the pTbFIX plasmid to perform RNAi experimentation. Lethal phenotypes, including α-tubulin and those corresponding to the enolase and clathrin heavy chain genes, were successfully recapitulated in procyclic and bloodstream parasites thus showing the versatility of this new tool.

Clinical and veterinary trypanocidal benzoxaboroles target CPSF3

African trypanosomes cause lethal and neglected tropical diseases, known as sleeping sickness in humans and nagana in animals. Current therapies are limited, but fortunately, promising therapies are in advanced clinical and veterinary development, including acoziborole (AN5568 or SCYX-7158) and AN11736, respectively. These benzoxaboroles will likely be key to the World Health Organization's target of disease control by 2030. Their mode of action was previously unknown. We have developed a high-coverage overexpression library and use it here to explore drug mode of action in Trypanosoma brucei Initially, an inhibitor with a known target was used to select for drug resistance and to test massive parallel library screening and genome-wide mapping; this effectively identified the known target and validated the approach. Subsequently, the overexpression screening approach was used to identify the target of the benzoxaboroles, Cleavage and Polyadenylation Specificity Factor 3 (CPSF3, Tb927.4.1340). We validated the CPSF3 endonuclease as the target, using independent overexpression strains. Knockdown provided genetic validation of CPSF3 as essential, and GFP tagging confirmed the expected nuclear localization. Molecular docking and CRISPR-Cas9-based editing demonstrated how acoziborole can specifically block the active site and mRNA processing by parasite, but not host CPSF3. Thus, our findings provide both genetic and chemical validation for CPSF3 as an important drug target in trypanosomes and reveal inhibition of mRNA maturation as the mode of action of the trypanocidal benzoxaboroles. Understanding the mechanism of action of benzoxaborole-based therapies can assist development of improved therapies, as well as the prediction and monitoring of resistance, if or when it arises.

The trypanocidal benzoxaborole AN7973 inhibits trypanosome mRNA processing

Kinetoplastid parasites-trypanosomes and leishmanias-infect millions of humans and cause economically devastating diseases of livestock, and the few existing drugs have serious deficiencies. Benzoxaborole-based compounds are very promising potential novel anti-trypanosomal therapies, with candidates already in human and animal clinical trials. We investigated the mechanism of action of several benzoxaboroles, including AN7973, an early candidate for veterinary trypanosomosis. In all kinetoplastids, transcription is polycistronic. Individual mRNA 5'-ends are created by trans splicing of a short leader sequence, with coupled polyadenylation of the preceding mRNA. Treatment of Trypanosoma brucei with AN7973 inhibited trans splicing within 1h, as judged by loss of the Y-structure splicing intermediate, reduced levels of mRNA, and accumulation of peri-nuclear granules. Methylation of the spliced leader precursor RNA was not affected, but more prolonged AN7973 treatment caused an increase in S-adenosyl methionine and methylated lysine. Together, the results indicate that mRNA processing is a primary target of AN7973. Polyadenylation is required for kinetoplastid trans splicing, and the EC50 for AN7973 in T. brucei was increased three-fold by over-expression of the T. brucei cleavage and polyadenylation factor CPSF3, identifying CPSF3 as a potential molecular target. Molecular modeling results suggested that inhibition of CPSF3 by AN7973 is feasible. Our results thus chemically validate mRNA processing as a viable drug target in trypanosomes. Several other benzoxaboroles showed metabolomic and splicing effects that were similar to those of AN7973, identifying splicing inhibition as a common mode of action and suggesting that it might be linked to subsequent changes in methylated metabolites. Granule formation, splicing inhibition and resistance after CPSF3 expression did not, however, always correlate and prolonged selection of trypanosomes in AN7973 resulted in only 1.5-fold resistance. It is therefore possible that the modes of action of oxaboroles that target trypanosome mRNA processing might extend beyond CPSF3 inhibition.

Comparative proteomics of the two T. brucei PABPs suggests that PABP2 controls bulk mRNA

Poly(A)-binding proteins (PABPs) regulate mRNA fate by controlling stability and translation through interactions with both the poly(A) tail and eIF4F complex. Many organisms have several paralogs of PABPs and eIF4F complex components and it is likely that different eIF4F/PABP complex combinations regulate distinct sets of mRNAs. Trypanosomes have five eIF4G paralogs, six of eIF4E and two PABPs, PABP1 and PABP2. Under starvation, polysomes dissociate and the majority of mRNAs, most translation initiation factors and PABP2 reversibly localise to starvation stress granules. To understand this more broadly we identified a protein interaction cohort for both T. brucei PABPs by cryo-mill/affinity purification-mass spectrometry. PABP1 very specifically interacts with the previously identified interactors eIF4E4 and eIF4G3 and few others. In contrast PABP2 is promiscuous, with a larger set of interactors including most translation initiation factors and most prominently eIF4G1, with its two partners TbG1-IP and TbG1-IP2. Only RBP23 was specific to PABP1, whilst 14 RNA-binding proteins were exclusively immunoprecipitated with PABP2. Significantly, PABP1 and associated proteins are largely excluded from starvation stress granules, but PABP2 and most interactors translocate to granules on starvation. We suggest that PABP1 regulates a small subpopulation of mainly small-sized mRNAs, as it interacts with a small and distinct set of proteins unable to enter the dominant pathway into starvation stress granules and localises preferentially to a subfraction of small polysomes. By contrast PABP2 likely regulates bulk mRNA translation, as it interacts with a wide range of proteins, enters stress granules and distributes over the full range of polysomes.

Recent advances in trypanosomatid research: genome organization, expression, metabolism, taxonomy and evolution

Unicellular flagellates of the family Trypanosomatidae are obligatory parasites of invertebrates, vertebrates and plants. Dixenous species are aetiological agents of a number of diseases in humans, domestic animals and plants. Their monoxenous relatives are restricted to insects. Because of the high biological diversity, adaptability to dramatically different environmental conditions, and omnipresence, these protists have major impact on all biotic communities that still needs to be fully elucidated. In addition, as these organisms represent a highly divergent evolutionary lineage, they are strikingly different from the common 'model system' eukaryotes, such as some mammals, plants or fungi. A number of excellent reviews, published over the past decade, were dedicated to specialized topics from the areas of trypanosomatid molecular and cell biology, biochemistry, host-parasite relationships or other aspects of these fascinating organisms. However, there is a need for a more comprehensive review that summarizing recent advances in the studies of trypanosomatids in the last 30 years, a task, which we tried to accomplish with the current paper.

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.

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.

Characterization of RBP9 and RBP10, two developmentally regulated RNA-binding proteins in Trypanosoma brucei

The fate of an mRNA is determined by its interaction with proteins and small RNAs within dynamic complexes called ribonucleoprotein complexes (mRNPs). In Trypanosoma brucei and related kinetoplastids, responses to internal and external signals are mainly mediated by post-transcriptional processes. Here, we used proximity-dependent biotin identification (BioID) combined with RNA-seq to investigate the changes resulting from ectopic expression of RBP10 and RBP9, two developmentally regulated RNA-binding proteins (RBPs). Both RBPs have reduced expression in insect procyclic forms (PCFs) compared with bloodstream forms (BSFs). Upon overexpression in PCFs, both proteins were recruited to cytoplasmic foci, co-localizing with the processing body marker SCD6. Further, both RBPs altered the transcriptome from a PCF- to a BSF-like pattern. Notably, upon expression of BirA*-RBP9 and BirA*-RBP10, BioID yielded more than 200 high confidence protein interactors (more than 10-fold enriched); 45 (RBP9) and 31 (RBP10) were directly related to mRNA metabolism. This study validates the use of BioID for investigating mRNP components but also illustrates the complexity of mRNP function.

Well-positioned nucleosomes punctuate polycistronic pol II transcription units and flank silent VSG gene arrays in Trypanosoma brucei

Background

The compaction of DNA in chromatin in eukaryotes allowed the expansion of genome size and coincided with significant evolutionary diversification. However, chromatin generally represses DNA function, and mechanisms coevolved to regulate chromatin structure and its impact on DNA. This included the selection of specific nucleosome positions to modulate accessibility to the DNA molecule. Trypanosoma brucei, a member of the Excavates supergroup, falls in an ancient evolutionary branch of eukaryotes and provides valuable insight into the organization of chromatin in early genomes.

Results

We have mapped nucleosome positions in T. brucei and identified important differences compared to other eukaryotes: The RNA polymerase II initiation regions in T. brucei do not exhibit pronounced nucleosome depletion, and show little evidence for defined −1 and +1 nucleosomes. In contrast, a well-positioned nucleosome is present directly on the splice acceptor sites within the polycistronic transcription units. The RNA polyadenylation sites were depleted of nucleosomes, with a single well-positioned nucleosome present immediately downstream of the predicted sites. The regions flanking the silent variant surface glycoprotein (VSG) gene cassettes showed extensive arrays of well-positioned nucleosomes, which may repress cryptic transcription initiation. The silent VSG genes themselves exhibited a less regular nucleosomal pattern in both bloodstream and procyclic form trypanosomes. The DNA replication origins, when present within silent VSG gene cassettes, displayed a defined nucleosomal organization compared with replication origins in other chromosomal core regions.

Conclusions

Our results indicate that some organizational features of chromatin are evolutionarily ancient, and may already have been present in the last eukaryotic common ancestor.

Electronic supplementary material

The online version of this article (doi:10.1186/s13072-017-0121-9) contains supplementary material, which is available to authorized users.

The eIF3 complex of Trypanosoma brucei: composition conservation does not imply the conservation of structural assembly and subunits function

The multisubunit eukaryotic initiation factor 3 (eIF3) plays multiple roles in translation but is poorly understood in trypanosomes. The putative subunits eIF3a and eIF3f of Trypanosoma brucei (TbIF3a and TbIF3f) were overexpressed and purified, and 11 subunits were identified, TbIF3a through l minus j, which form a tight complex. Both TbIF3a and TbIF3f are essential for the viability of T. brucei RNAi knockdown of either of them severely reduced total translation and the ratio of the polysome/80S peak area. TbIF3f and TbIF3a RNAi cell lines were modified to express tagged-TbIF3a and -TbIF3f, respectively. RNAi in combination with affinity purification assays indicated that both subunits are variably required for TbIF3 stability and integrity. The relative abundance of other subunits in the TbIF3f-tag complex changed little upon TbIF3a depletion; while only subunits TbIF3b, i, and e copurified comparably with TbIF3a-tag upon TbIF3f depletion. A genome-wide UV-crosslinking assay showed that several TbIF3 subunits have direct RNA-binding activity, with TbIF3c showing the strongest signal. In addition, CrPV IRES, but neither EMCV IRES nor HCV IRES, was found to mediate translation in T. brucei These results together imply that the structure of TbIF3 and the subunits function have trypanosome-specific features, although the composition is evolutionarily conserved.

Dissecting biochemical peculiarities of the ATPase activity of TcSub2, a component of the mRNA export pathway in Trypanosoma cruzi

The RNA helicase DEAD-box protein Sub2 (yeast)/UAP56 (mammals) is conserved across eukaryotes and is essential for mRNA export in trypanosomes. Despite the high conservation of Sub2 in lower eukaryotes such as Trypanosoma cruzi, the low conservation of other mRNA export factors raises questions regarding whether the mode of action of TcSub2 is similar to that of orthologs from other eukaryotes. Mutation of the conserved K87 residue of TcSub2 abolishes ATPase activity, showing that its ATPase domain is functional. However, the Vmax of TcSub2 was much higher than the Vmax described for the human protein UAP56, which suggests that the TcSub2 enzyme hydrolyzes ATP faster than its human homolog. Furthermore, we demonstrate that RNA association is less important to the activity of TcSub2 compared to UAP56. Our results show differences in activity of this protein, even though the structure of TcSub2 is very similar to UAP56. Functional complementation assays indicate that these differences may be common to other trypanosomatids. Distinct features of RNA influence and ATPase efficiency between UAP56 and TcSub2 may reflect distinct structures for functional sites of TcSub2. For this reason, ligand-based or structure-based methodologies can be applied to investigate the potential of TcSub2 as a target for new drugs.

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.

High-resolution analysis of multi-copy variant surface glycoprotein gene expression sites in African trypanosomes

Background

African trypanosomes cause lethal diseases in humans and animals and escape host immune attack by switching the expression of Variant Surface Glycoprotein (VSG) genes. The expressed VSGs are located at the ends of telomeric, polycistronic transcription units known as VSG expression sites (VSG-ESs). Each cell has many VSG-ESs but only one is transcribed in bloodstream-form parasites and all of them are inactive upon transmission to the insect vector mid-gut; a subset of monocistronic metacyclic VSG-ESs are then activated in the insect salivary gland. Deep-sequence analyses have been informative but assigning sequences to individual VSG-ESs has been challenging because they each contain closely related expression-site associated genes, or ESAGs, thought to contribute to virulence.

Results

We utilised ART, an in silico short read simulator to demonstrate the feasibility of accurately aligning reads to VSG-ESs. Then, using high-resolution transcriptomes from isogenic bloodstream and insect-stage Lister 427 Trypanosoma brucei, we uncover increased abundance in the insect mid-gut stage of mRNAs from metacyclic VSG-ESs and of mRNAs from the unusual ESAG, ESAG10. Further, we show that the silencing associated with allelic exclusion involves repression focussed at the ends of the VSG-ESs. We also use the approach to report relative fitness costs following ESAG RNAi from a genome-scale screen.

Conclusions

By assigning sequences to individual VSG-ESs we provide new insights into VSG-ES transcription control, allelic exclusion and impacts on fitness. Thus, deeper insights into the expression and function of regulated multi-gene families are more accessible than previously anticipated.

Electronic supplementary material

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

Plasticity of the Leishmania genome leading to gene copy number variations and drug resistance

Leishmania has a plastic genome, and drug pressure can select for gene copy number variation (CNV). CNVs can apply either to whole chromosomes, leading to aneuploidy, or to specific genomic regions. For the latter, the amplification of chromosomal regions occurs at the level of homologous direct or inverted repeated sequences leading to extrachromosomal circular or linear amplified DNAs. This ability of Leishmania to respond to drug pressure by CNVs has led to the development of genomic screens such as Cos-Seq, which has the potential of expediting the discovery of drug targets for novel promising drug candidates.

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.

Base J and H3.V Regulate Transcriptional Termination in Trypanosoma brucei

Trypanosoma brucei is a protozoan parasite that lacks many transcription factors found in other eukaryotes, such as those whose binding demarcates enhancers. T. brucei retains histone variants and modifications, however, and it is hypothesized that it relies on epigenetic marks to define transcription-related boundaries. The histone H3 variant (H3.V) and an alternate nucleotide, base J (ß-D-glucosyl-hydroxymethyluracil), are two chromatin marks found at both transcription termination sites (TTSs) and telomeres. Here, we report that the absence of both base J and H3.V result in transcription readthrough and the appearance of antisense transcripts near TTSs. Additionally, we find that maintaining the transcriptional silencing of pol I-transcribed telomeric Variant Surface Glycoprotein (VSG) genes appears to be dependent on deposition of H3.V alone. Our study reveals that gene expression depends on different epigenetic cues depending on chromosomal location and on the transcribing polymerase. This work provides insight into how these signals may have evolved into the more nuanced and fine-tuned gene regulatory mechanisms observed in other model systems.

In eukaryotes, DNA is wrapped around histones to form chromatin. Modifications on the DNA itself, or on the canonical histones around which the DNA is wrapped, can lead to changes in gene expression. Alternate histones are also used to regulate gene expression. The African trypanosome, the causative agent of African sleeping sickness, transcribes its genes in long arrays called polycistronic transcription units (PTUs). In this study, we describe how the African Trypanosome uses two chromatin marks to regulate genes that lie close to the end of the PTU or close to the ends of chromosomes, called telomeres. One mark, base J, is on the DNA itself, while the other utilized mark is an alternate histone protein called H3.V. In the absence of these marks, there is an increase in antisense RNA that corresponds to genes that lie near the ends of the PTUs, and an increase in sense RNA for genes near telomeres. Since trypanosomes are evolutionarily distant from humans, these findings shed light on how gene expression mechanisms that are modulated by chromatin marks may have evolved to produce the complex gene regulatory networks found in our own tissues.

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.

Simultaneous transcriptional profiling of Leishmania major and its murine macrophage host cell reveals insights into host-pathogen interactions

BACKGROUND

Parasites of the genus Leishmania are the causative agents of leishmaniasis, a group of diseases that range in manifestations from skin lesions to fatal visceral disease. The life cycle of Leishmania parasites is split between its insect vector and its mammalian host, where it resides primarily inside of macrophages. Once intracellular, Leishmania parasites must evade or deactivate the host's innate and adaptive immune responses in order to survive and replicate.

RESULTS

We performed transcriptome profiling using RNA-seq to simultaneously identify global changes in murine macrophage and L. major gene expression as the parasite entered and persisted within murine macrophages during the first 72 h of an infection. Differential gene expression, pathway, and gene ontology analyses enabled us to identify modulations in host and parasite responses during an infection. The most substantial and dynamic gene expression responses by both macrophage and parasite were observed during early infection. Murine genes related to both pro- and anti-inflammatory immune responses and glycolysis were substantially upregulated and genes related to lipid metabolism, biogenesis, and Fc gamma receptor-mediated phagocytosis were downregulated. Upregulated parasite genes included those aimed at mitigating the effects of an oxidative response by the host immune system while downregulated genes were related to translation, cell signaling, fatty acid biosynthesis, and flagellum structure.

CONCLUSIONS

The gene expression patterns identified in this work yield signatures that characterize multiple developmental stages of L. major parasites and the coordinated response of Leishmania-infected macrophages in the real-time setting of a dual biological system. This comprehensive dataset offers a clearer and more sensitive picture of the interplay between host and parasite during intracellular infection, providing additional insights into how pathogens are able to evade host defenses and modulate the biological functions of the cell in order to survive in the mammalian environment.

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.

NMD3 regulates both mRNA and rRNA nuclear export in African trypanosomes via an XPOI-linked pathway

Trypanosomes mostly regulate gene expression through post-transcriptional mechanisms, particularly mRNA stability. However, much mRNA degradation is cytoplasmic such that mRNA nuclear export must represent an important level of regulation. Ribosomal RNAs must also be exported from the nucleus and the trypanosome orthologue of NMD3 has been confirmed to be involved in rRNA processing and export, matching its function in other organisms. Surprisingly, we found that TbNMD3 depletion also generates mRNA accumulation of procyclin-associated genes (PAGs), these being co-transcribed by RNA polymerase I with the procyclin surface antigen genes expressed on trypanosome insect forms. By whole transcriptome RNA-seq analysis of TbNMD3-depleted cells we confirm the regulation of the PAG transcripts by TbNMD3 and using reporter constructs reveal that PAG1 regulation is mediated by its 5'UTR. Dissection of the mechanism of regulation demonstrates that it is not dependent upon translational inhibition mediated by TbNMD3 depletion nor enhanced transcription. However, depletion of the nuclear export factors XPO1 or MEX67 recapitulates the effects of TbNMD3 depletion on PAG mRNAs and mRNAs accumulated in the nucleus of TbNMD3-depleted cells. These results invoke a novel RNA regulatory mechanism involving the NMD3-dependent nuclear export of mRNA cargos, suggesting a shared platform for mRNA and rRNA export.

25 years of African trypanosome research: From description to molecular dissection and new drug discovery

The Molecular Parasitology conference was first held at the Marine Biological laboratory, Woods Hole, USA 25 years ago. Since that first meeting, the conference has evolved and expanded but has remained the showcase for the latest research developments in molecular parasitology. In this perspective, I reflect on the scientific discoveries focussed on African trypanosomes (Trypanosoma brucei spp.) that have occurred since the inaugural MPM meeting and discuss the current and future status of research on these parasites.

Species-specific antimonial sensitivity in Leishmania is driven by post-transcriptional regulation of AQP1

Leishmania is a digenetic protozoan parasite causing leishmaniasis in humans. The different clinical forms of leishmaniasis are caused by more than twenty species of Leishmania that are transmitted by nearly thirty species of phlebotomine sand flies. Pentavalent antimonials (such as Pentostam or Glucantime) are the first line drugs for treating leishmaniasis. Recent studies suggest that pentavalent antimony (Sb(V)) acts as a pro-drug, which is converted to the more active trivalent form (Sb(III)). However, sensitivity to trivalent antimony varies among different Leishmania species. In general, Leishmania species causing cutaneous leishmaniasis (CL) are more sensitive to Sb(III) than the species responsible for visceral leishmaniasis (VL). Leishmania aquaglyceroporin (AQP1) facilitates the adventitious passage of antimonite down a concentration gradient. In this study, we show that Leishmania species causing CL accumulate more antimonite, and therefore exhibit higher sensitivity to antimonials, than the species responsible for VL. This species-specific differential sensitivity to antimonite is directly proportional to the expression levels of AQP1 mRNA. We show that the stability of AQP1 mRNA in different Leishmania species is regulated by their respective 3’-untranslated regions. The differential regulation of AQP1 mRNA explains the distinct antimonial sensitivity of each species.

The degree of response to antimonial drugs varies widely between species and even among strains of the same species of the protozoan parasite Leishmania. However, the molecular mechanism(s) is unknown. In this study, we show that Leishmania aquaglyceroporin AQP1 drives this species-specific antimonial resistance. Aquaglyceroporins are channel proteins that facilitate the passage of small uncharged molecules, such as glycerol and water, across the biological membranes. AQP1 helps the parasite cope with the osmotic challenges it faces during its life cycle. Additionally, AQP1 is an adventitious facilitator of antimonite, the active form of pentavalent antimonial drugs. We show that AQP1 expression level is species-specific, and less AQP1 in visceral species compared to the cutaneous species results in increased resistance to antimonials. We also demonstrate that the 3’-untranslated regions (3’-UTR) of the AQP1 mRNA is a major determining factor of species-specific regulation of AQP1. Along with water homeostasis, aquaglyceroporins are also involved in directed cell migration. The variable levels of AQP1 in different Leishmania species may enable them to find their appropriate niches in vertebrate hosts and cope with the species-specific osmotic challenges during their life cycles.

Unveiling the intracellular survival gene kit of trypanosomatid parasites

Trypanosomatids are unicellular protozoans of medical and economical relevance since they are the etiologic agents of infectious diseases in humans as well as livestock. Whereas Trypanosoma cruzi and different species of Leishmania are obligate intracellular parasites, Trypanosoma brucei and other trypanosomatids develop extracellularly throughout their entire life cycle. After their genomes have been sequenced, various comparative genomic studies aimed at identifying sequences involved with host cell invasion and intracellular survival have been described. However, for only a handful of genes, most of them present exclusively in the T. cruzi or Leishmania genomes, has there been any experimental evidence associating them with intracellular parasitism. With the increasing number of published complete genome sequences of members of the trypanosomatid family, including not only different Trypanosoma and Leishmania strains and subspecies but also trypanosomatids that do not infect humans or other mammals, we may now be able to contemplate a slightly better picture regarding the specific set of parasite factors that defines each organism's mode of living and the associated disease phenotypes. Here, we review the studies concerning T. cruzi and Leishmania genes that have been implicated with cell invasion and intracellular parasitism and also summarize the wealth of new information regarding the mode of living of intracellular parasites that is resulting from comparative genome studies that are based on increasingly larger trypanosomatid genome datasets.

Identification of a novel nucleocytoplasmic shuttling RNA helicase of trypanosomes

Gene expression in trypanosomes is controlled mostly by post-transcriptional pathways. Little is known about the components of mRNA nucleocytoplasmic export routes in these parasites. Comparative genomics has shown that the mRNA transport pathway is the least conserved pathway among eukaryotes. Nonetheless, we identified a RNA helicase (Hel45) that is conserved across eukaryotes and similar to shuttling proteins involved in mRNA export. We used in silico analysis to predict the structure of Trypanosoma cruzi Hel45, including the N-terminal domain and the C-terminal domain, and our findings suggest that this RNA helicase can form complexes with mRNA. Hel45 was present in both nucleus and cytoplasm. Electron microscopy showed that Hel45 is clustered close to the cytoplasmic side of nuclear pore complexes, and is also present in the nucleus where it is associated with peripheral compact chromatin. Deletion of a predicted Nuclear Export Signal motif led to the accumulation of Hel45ΔNES in the nucleus, indicating that Hel45 shuttles between the nucleus and the cytoplasm. This transport was dependent on active transcription but did not depend on the exportin Crm1. Knockdown of Mex67 in T. brucei caused the nuclear accumulation of the T. brucei ortholog of Hel45. Indeed, Hel45 is present in mRNA ribonucleoprotein complexes that are not associated with polysomes. It is still necessary to confirm the precise function of Hel45. However, this RNA helicase is associated with mRNA metabolism and its nucleocytoplasmic shuttling is dependent on an mRNA export route involving Mex67 receptor.

Proteomic analysis of Trypanosoma cruzi response to ionizing radiation stress

Trypanosoma cruzi, the causative agent of Chagas disease, is extremely resistant to ionizing radiation, enduring up to 1.5 kGy of gamma rays. Ionizing radiation can damage the DNA molecule both directly, resulting in double-strand breaks, and indirectly, as a consequence of reactive oxygen species production. After a dose of 500 Gy of gamma rays, the parasite genome is fragmented, but the chromosomal bands are restored within 48 hours. Under such conditions, cell growth arrests for up to 120 hours and the parasites resume normal growth after this period. To better understand the parasite response to ionizing radiation, we analyzed the proteome of irradiated (4, 24, and 96 hours after irradiation) and non-irradiated T. cruzi using two-dimensional differential gel electrophoresis followed by mass spectrometry for protein identification. A total of 543 spots were found to be differentially expressed, from which 215 were identified. These identified protein spots represent different isoforms of only 53 proteins. We observed a tendency for overexpression of proteins with molecular weights below predicted, indicating that these may be processed, yielding shorter polypeptides. The presence of shorter protein isoforms after irradiation suggests the occurrence of post-translational modifications and/or processing in response to gamma radiation stress. Our results also indicate that active translation is essential for the recovery of parasites from ionizing radiation damage. This study therefore reveals the peculiar response of T. cruzi to ionizing radiation, raising questions about how this organism can change its protein expression to survive such a harmful stress.

Trypanosoma brucei translation initiation factor homolog EIF4E6 forms a tripartite cytosolic complex with EIF4G5 and a capping enzyme homolog

Trypanosomes lack the transcriptional control characteristic of the majority of eukaryotes that is mediated by gene-specific promoters in a one-gene-one-promoter arrangement. Rather, their genomes are transcribed in large polycistrons with no obvious functional linkage. Posttranscriptional regulation of gene expression must thus play a larger role in these organisms. The eIF4E homolog TbEIF4E6 binds mRNA cap analogs in vitro and is part of a complex in vivo that may fulfill such a role. Knockdown of TbEIF4E6 tagged with protein A-tobacco etch virus protease cleavage site-protein C to approximately 15% of the normal expression level resulted in viable cells that displayed a set of phenotypes linked to detachment of the flagellum from the length of the cell body, if not outright flagellum loss. While these cells appeared and behaved as normal under stationary liquid culture conditions, standard centrifugation resulted in a marked increase in flagellar detachment. Furthermore, the ability of TbEIF4E6-depleted cells to engage in social motility was reduced. The TbEIF4E6 protein forms a cytosolic complex containing a triad of proteins, including the eIF4G homolog TbEIF4G5 and a hypothetical protein of 70.3 kDa, referred to as TbG5-IP. The TbG5-IP analysis revealed two domains with predicted secondary structures conserved in mRNA capping enzymes: nucleoside triphosphate hydrolase and guanylyltransferase. These complex members have the potential for RNA interaction, either via the 5' cap structure for TbEIF4E6 and TbG5-IP or through RNA-binding domains in TbEIF4G5. The associated proteins provide a signpost for future studies to determine how this complex affects capped RNA molecules.

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

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

A genome-wide analysis of genetic diversity in Trypanosoma cruzi intergenic regions

BACKGROUND

Trypanosoma cruzi is the causal agent of Chagas Disease. Recently, the genomes of representative strains from two major evolutionary lineages were sequenced, allowing the construction of a detailed genetic diversity map for this important parasite. However this map is focused on coding regions of the genome, leaving a vast space of regulatory regions uncharacterized in terms of their evolutionary conservation and/or divergence.

METHODOLOGY

Using data from the hybrid CL Brener and Sylvio X10 genomes (from the TcVI and TcI Discrete Typing Units, respectively), we identified intergenic regions that share a common evolutionary ancestry, and are present in both CL Brener haplotypes (TcII-like and TcIII-like) and in the TcI genome; as well as intergenic regions that were conserved in only two of the three genomes/haplotypes analyzed. The genetic diversity in these regions was characterized in terms of the accumulation of indels and nucleotide changes.

PRINCIPAL FINDINGS

Based on this analysis we have identified i) a core of highly conserved intergenic regions, which remained essentially unchanged in independently evolving lineages; ii) intergenic regions that show high diversity in spite of still retaining their corresponding upstream and downstream coding sequences; iii) a number of defined sequence motifs that are shared by a number of unrelated intergenic regions. A fraction of indels explains the diversification of some intergenic regions by the expansion/contraction of microsatellite-like repeats.

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.

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.