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

How much (ATP) does it cost to build a trypanosome? A theoretical study on the quantity of ATP needed to maintain and duplicate a bloodstream-form Trypanosoma brucei cell

ATP hydrolysis is required for the synthesis, transport and polymerization of monomers for macromolecules as well as for the assembly of the latter into cellular structures. Other cellular processes not directly related to synthesis of biomass, such as maintenance of membrane potential and cellular shape, also require ATP. The unicellular flagellated parasite Trypanosoma brucei has a complex digenetic life cycle. The primary energy source for this parasite in its bloodstream form (BSF) is glucose, which is abundant in the host's bloodstream. Here, we made a detailed estimation of the energy budget during the BSF cell cycle. As glycolysis is the source of most produced ATP, we calculated that a single parasite produces 6.0 x 1011 molecules of ATP/cell cycle. Total biomass production (which involves biomass maintenance and duplication) accounts for ~63% of the total energy budget, while the total biomass duplication accounts for the remaining ~37% of the ATP consumption, with in both cases translation being the most expensive process. These values allowed us to estimate a theoretical YATP of 10.1 (g biomass)/mole ATP and a theoretical [Formula: see text] of 28.6 (g biomass)/mole ATP. Flagellar motility, variant surface glycoprotein recycling, transport and maintenance of transmembrane potential account for less than 30% of the consumed ATP. Finally, there is still ~5.5% available in the budget that is being used for other cellular processes of as yet unknown cost. These data put a new perspective on the assumptions about the relative energetic weight of the processes a BSF trypanosome undergoes during its cell cycle.

The RNA export factor TbMex67 connects transcription and RNA export in Trypanosoma brucei and sets boundaries for RNA polymerase I

TbMex67 is the major mRNA export factor known to date in trypanosomes, forming part of the docking platform within the nuclear pore. To explore its role in co-transcriptional mRNA export, recently reported in Trypanosoma brucei, pulse labelling of nascent RNAs with 5-ethynyl uridine (5-EU) was performed with cells depleted of TbMex67 and complemented with a dominant-negative mutant (TbMex67-DN). RNA polymerase (Pol) II transcription was unaffected, but the procyclin loci, which encode mRNAs transcribed by Pol I from internal sites on chromosomes 6 and 10, showed increased levels of 5-EU incorporation. This was due to Pol I readthrough transcription, which proceeded beyond the procyclin and procyclin-associated genes up to the Pol II transcription start site on the opposite strand. Complementation by TbMex67-DN also increased Pol I-dependent formation of R-loops and γ-histone 2A foci. The DN mutant exhibited reduced nuclear localisation and binding to chromatin compared to wild-type TbMex67. Together with its interaction with chromatin remodelling factor TbRRM1 and Pol II, and transcription-dependent association of Pol II with nucleoporins, our findings support a role for TbMex67 in connecting transcription and export in T. brucei. In addition, TbMex67 stalls readthrough by Pol I in specific contexts, thereby limiting R-loop formation and replication stress.

Two cold shock domain containing proteins trigger the development of infectious Trypanosoma brucei

Cold shock proteins are members of a family of DNA- and RNA-binding proteins with one or more evolutionarily conserved cold shock domain (CSD). These proteins have a wide variety of biological functions, including DNA-damage repair, mRNA stability, and regulation of transcription, splicing and translation. We previously identified two CSD containing proteins, CSD1 and CSD2, in the protozoan parasite Trypanosoma brucei to be required for RBP6-driven metacyclic production, albeit at different steps of the developmental program. During metacyclogenesis T. brucei undergoes major morphological and metabolic changes that culminate in the establishment of quiescent metacyclic parasites and the acquisition of mammalian infectivity. To investigate the specific role of CSD1 and CSD2 in this process, we ectopically expressed CSD1 or CSD2 in non-infectious procyclic parasites and discovered that each protein is sufficient to produce infectious metacyclic parasites in 24 hours. Domain truncation assays determined that the N-terminal domain, but not the C-terminal domain, of CSD1 and CSD2 was required for metacyclic development. Furthermore, conserved amino acid residues in the CSD of CSD1 and CSD2, known to be important for binding nucleic acids, were found to be necessary for metacyclic production. Using single-end enhanced crosslinking and immunoprecipitation (seCLIP) we identified the specific binding motif of CSD1 and CSD2 as "ANACAU" and the bound mRNAs were enriched for biological processes, including lipid metabolism, microtubule-based movement and nucleocytoplasmic transport that are likely involved in the transition to bloodstream form-like cells.

Subcellular protein localisation of Trypanosoma brucei bloodstream form-upregulated proteins maps stage-specific adaptations

Background: Genome-wide subcellular protein localisation in Trypanosoma brucei, through our TrypTag project, has comprehensively dissected the molecular organisation of this important pathogen. Powerful as this resource is , T. brucei has multiple developmental forms and we previously only analysed the procyclic form. This is an insect life cycle stage, leaving the mammalian bloodstream form unanalysed. The expectation is that between life stages protein localisation would not change dramatically (completely unchanged or shifting to analogous stage-specific structures). However, this has not been specifically tested. Similarly, which organelles tend to contain proteins with stage-specific expression can be predicted from known stage specific adaptations but has not been comprehensively tested. Methods: We used endogenous tagging with mNG to determine the sub-cellular localisation of the majority of proteins encoded by transcripts significantly upregulated in the bloodstream form, and performed comparison to the existing localisation data in procyclic forms. Results: We have confirmed the localisation of known stage-specific proteins and identified the localisation of novel stage-specific proteins. This gave a map of which organelles tend to contain stage specific proteins: the mitochondrion for the procyclic form, and the endoplasmic reticulum, endocytic system and cell surface in the bloodstream form. Conclusions: This represents the first genome-wide map of life cycle stage-specific adaptation of organelle molecular machinery in T. brucei.

Profiling the bloodstream form and procyclic form Trypanosoma brucei cell cycle using single-cell transcriptomics

African trypanosomes proliferate as bloodstream forms (BSFs) and procyclic forms in the mammal and tsetse fly midgut, respectively. This allows them to colonise the host environment upon infection and ensure life cycle progression. Yet, understanding of the mechanisms that regulate and drive the cell replication cycle of these forms is limited. Using single-cell transcriptomics on unsynchronised cell populations, we have obtained high resolution cell cycle regulated (CCR) transcriptomes of both procyclic and slender BSF Trypanosoma brucei without prior cell sorting or synchronisation. Additionally, we describe an efficient freeze-thawing protocol that allows single-cell transcriptomic analysis of cryopreserved T. brucei. Computational reconstruction of the cell cycle using periodic pseudotime inference allowed the dynamic expression patterns of cycling genes to be profiled for both life cycle forms. Comparative analyses identify a core cycling transcriptome highly conserved between forms, as well as several genes where transcript levels dynamics are form specific. Comparing transcript expression patterns with protein abundance revealed that the majority of genes with periodic cycling transcript and protein levels exhibit a relative delay between peak transcript and protein expression. This work reveals novel detail of the CCR transcriptomes of both forms, which are available for further interrogation via an interactive webtool.

Cell-to-flagellum attachment and surface architecture in kinetoplastids

A key morphological feature of kinetoplastid parasites is the position and length of flagellum attachment to the cell body. This lateral attachment is mediated by the flagellum attachment zone (FAZ), a large complex cytoskeletal structure, which is essential for parasite morphogenesis and pathogenicity. Despite the complexity of the FAZ only two transmembrane proteins, FLA1 and FLA1BP, are known to interact and connect the flagellum to the cell body. Across the different kinetoplastid species, each only has a single FLA/FLABP pair, except in Trypanosoma brucei and Trypanosoma congolense where there has been an expansion of these genes. Here, we focus on the selection pressure behind the evolution of the FLA/FLABP proteins and the likely impact this will have on host-parasite interactions.

Nanopore-Based Direct RNA Sequencing of the Trypanosoma brucei Transcriptome Identifies Novel lncRNAs

Trypanosomatids are single-cell eukaryotic parasites. Unlike higher eukaryotes, they control gene expression post-transcriptionally and not at the level of transcription initiation. This involves all known cellular RNA circuits, from mRNA processing to mRNA decay, to translation, in addition to a large panel of RNA-interacting proteins that modulate mRNA abundance. However, other forms of gene regulation, for example by lncRNAs, cannot be excluded. LncRNAs are poorly studied in trypanosomatids, with only a single lncRNA characterized to date. Furthermore, it is not clear whether the complete inventory of trypanosomatid lncRNAs is known, because of the inherent cDNA-recoding and DNA-amplification limitations of short-read RNA sequencing. Here, we overcome these limitations by using long-read direct RNA sequencing (DRS) on nanopore arrays. We analyze the native RNA pool of the two main lifecycle stages of the African trypanosome Trypanosoma brucei, with a special emphasis on the inventory of lncRNAs. We identify 207 previously unknown lncRNAs, 32 of which are stage-specifically expressed. We also present insights into the complexity of the T. brucei transcriptome, including alternative transcriptional start and stop sites and potential transcript isoforms, to provide a bias-free understanding of the intricate RNA landscape in T. brucei.

Comparative genomics and interactomics of polyadenylation factors for the prediction of new parasite targets: Entamoeba histolytica as a working model

Protein-protein interactions (PPI) play a key role in predicting the function of a target protein and drug ability to affect an entire biological system. Prediction of PPI networks greatly contributes to determine a target protein and signal pathways related to its function. Polyadenylation of mRNA 3'-end is essential for gene expression regulation and several polyadenylation factors have been shown as valuable targets for controlling protozoan parasites that affect human health. Here, by using a computational strategy based on sequence-based prediction approaches, phylogenetic analyses, and computational prediction of PPI networks, we compared interactomes of polyadenylation factors in relevant protozoan parasites and the human host, to identify key proteins and define potential targets for pathogen control. Then, we used Entamoeba histolytica as a working model to validate our computational results. RT-qPCR assays confirmed the coordinated modulation of connected proteins in the PPI network and evidenced that silencing of the bottleneck protein EhCFIm25 affects the expression of interacting proteins. In addition, molecular dynamics simulations and docking approaches allowed to characterize the relationships between EhCFIm25 and Ehnopp34, two connected bottleneck proteins. Interestingly, the experimental identification of EhCFIm25 interactome confirmed the close relationships among proteins involved in gene expression regulation and evidenced new links with moonlight proteins in E. histolytica, suggesting a connection between RNA biology and metabolism as described in other organisms. Altogether, our results strengthened the relevance of comparative genomics and interactomics of polyadenylation factors for the prediction of new targets for the control of these human pathogens.

Construction and Application of Orthogonal T7 Expression System in Eukaryote: An Overview

The T7 system is an orthogonal transcription-system, which is characterized by simplicity, higher efficiency, and higher processivity, and it is used for protein or mRNA synthesis in various biological-systems. In comparison with prokaryotes, the construction of the T7 expression system is still on-going in eukaryotes, but it shows greatly applicable prospects. In the present paper, development of T7 expression system construction in eukaryotes is reviewed, including its construction in animal (mammalian cells, trypanosomatid protozoa, Xenopus oocytes, zebrafish), plant, and microorganism and its application in vaccine production and gene therapy. In addition, the innate challenges of T7 expression system construction in eukaryote and its potential application in vaccine production and gene therapy are discussed.

Subcellular protein localisation of Trypanosoma brucei bloodstream form-upregulated proteins maps stage-specific adaptations

Background: Genome-wide subcellular protein localisation in Trypanosoma brucei, through our TrypTag project, has comprehensively dissected the molecular organisation of this important pathogen. Powerful as this resource is, T. brucei has multiple developmental forms and we previously only analysed the procyclic form. This is an insect life cycle stage, leaving the mammalian bloodstream form unanalysed. The expectation is that between life stages protein localisation would not change dramatically (completely unchanged or shifting to analogous stage-specific structures). However, this has not been specifically tested. Similarly, which organelles tend to contain proteins with stage-specific expression can be predicted from known stage specific adaptations but has not been comprehensively tested. Methods: We used endogenous tagging with mNG to determine the sub-cellular localisation of the majority of proteins encoded by transcripts significantly upregulated in the bloodstream form, and performed comparison to the existing localisation data in procyclic forms. Results: We have confirmed the localisation of known and identified the localisation of novel stage-specific proteins. This gave a map of which organelles tend to contain stage specific proteins: the mitochondrion for the procyclic form, and the endoplasmic reticulum, endocytic system and cell surface in the bloodstream form. Conclusions: This represents the first genome-wide map of life cycle stage-specific adaptation of organelle molecular machinery in T. brucei.

The RNA-binding protein RBP33 dampens non-productive transcription in trypanosomes

In-depth analysis of the transcriptomes of several model organisms has revealed that genomes are pervasively transcribed, giving rise to an abundance of non-canonical and mainly antisense RNA polymerase II-derived transcripts that are produced from almost any genomic context. Pervasive RNAs are degraded by surveillance mechanisms, but the repertoire of proteins that control the fate of these non-productive transcripts is still incomplete. Trypanosomes are single-celled eukaryotes that show constitutive RNA polymerase II transcription and in which initiation and termination of transcription occur at a limited number of sites per chromosome. It is not known whether pervasive transcription exists in organisms with unregulated RNA polymerase II activity, and which factors could be involved in the process. We show here that depletion of RBP33 results in overexpression of ∼40% of all annotated genes in the genome, with a marked accumulation of sense and antisense transcripts derived from silenced regions. RBP33 loss does not result in a significant increase in chromatin accessibility. Finally, we have found that transcripts that increase in abundance upon RBP33 knockdown are significantly more stable in RBP33-depleted trypanosomes, and that the exosome complex is responsible for their degradation. Our results provide strong evidence that RBP33 dampens non-productive transcription in trypanosomes.

Diversity in new flagellum tip attachment in bloodstream form African trypanosomes

The closely related parasites Trypanosoma brucei, T. congolense, and T. vivax cause neglected tropical diseases collectively known as African Trypanosomiasis. A characteristic feature of bloodstream form T. brucei is the flagellum that is laterally attached to the side of the cell body. During the cell cycle, the new flagellum is formed alongside the old flagellum, with the new flagellum tip embedded within a mobile transmembrane junction called the groove. The molecular composition of the groove is currently unknown, which limits the analysis of this junction and assessment of its conservation in related trypanosomatids. Here, we identified 13 proteins that localize to the flagellar groove through a small-scale tagging screen. Functional analysis of a subset of these proteins by RNAi and gene deletion revealed three proteins, FCP4/TbKin15, FCP7, and FAZ45, that are involved in new flagellum tip attachment to the groove. Despite possessing orthologues of all 13 groove proteins, T. congolense and T. vivax did not assemble a canonical groove around the new flagellum tip according to 3D electron microscopy. This diversity in new flagellum tip attachment points to the rapid evolution of membrane-cytoskeleton structures that can occur without large changes in gene complement and likely reflects the niche specialization of each species.

Distinct mRNA and protein interactomes highlight functional differentiation of major eIF4F-like complexes from Trypanosoma brucei

Gene expression in pathogenic protozoans of the family Trypanosomatidae has several novel features, including multiple eIF4F-like complexes involved in protein synthesis. The eukaryotic eIF4F complex, formed mainly by eIF4E and eIF4G subunits, is responsible for the canonical selection of mRNAs required for the initiation of mRNA translation. The best-known complexes implicated in translation in trypanosomatids are based on two related pairs of eIF4E and eIF4G subunits (EIF4E3/EIF4G4 and EIF4E4/EIF4G3), whose functional distinctions remain to be fully described. Here, to define interactomes associated with both complexes in Trypanosoma brucei procyclic forms, we performed parallel immunoprecipitation experiments followed by identification of proteins co-precipitated with the four tagged eIF4E and eIF4G subunits. A number of different protein partners, including RNA binding proteins and helicases, specifically co-precipitate with each complex. Highlights with the EIF4E4/EIF4G3 pair include RBP23, PABP1, EIF4AI and the CRK1 kinase. Co-precipitated partners with the EIF4E3/EIF4G4 pair are more diverse and include DRBD2, PABP2 and different zinc-finger proteins and RNA helicases. EIF4E3/EIF4G4 are essential for viability and to better define their role, we further investigated their phenotypes after knockdown. Depletion of either EIF4E3/EIF4G4 mRNAs lead to aberrant morphology with a more direct impact on events associated with cytokinesis. We also sought to identify those mRNAs differentially associated with each complex through CLIP-seq with the two eIF4E subunits. Predominant among EIF4E4-bound transcripts are those encoding ribosomal proteins, absent from those found with EIF4E3, which are generally more diverse. RNAi mediated depletion of EIF4E4, which does not affect proliferation, does not lead to changes in mRNAs or proteins associated with EIF4E3, confirming a lack of redundancy and distinct roles for the two complexes.

Roles and interactions of the specialized initiation factors EIF4E2, EIF4E5 and EIF4E6 in Trypanosoma brucei: EIF4E2 maintains the abundances of S-phase mRNAs

Trypanosoma brucei has six versions of the cap-binding translation initiation factor EIF4E. We investigated the functions of EIF4E2, EIF4E3, EIF4E5 and EIF4E6 in bloodstream forms. We confirmed the protein associations previously found in procyclic forms, and detected specific co-purification of some RNA-binding proteins. Bloodstream forms lacking EIF4E5 grew normally and differentiated to replication-incompetent procyclic forms. Depletion of EIF4E6 inhibited bloodstream-form trypanosome growth and translation. EIF4E2 co-purified only the putative RNA binding protein SLBP2. Bloodstream forms lacking EIF4E2 multiplied slowly, had a low maximal cell density, and expressed the stumpy-form marker PAD1, but showed no evidence for enhanced stumpy-form signalling. EIF4E2 knock-out cells differentiated readily to replication-competent procyclic forms. EIF4E2 was strongly associated with a subset of mRNAs that are maximally abundant in S-phase, and these all had decreased abundances in EIF4E2 knock-out cells. Three EIF4E2 target mRNAs are also bound and stabilized by the Pumilio domain protein PUF9. Yeast 2-hybrid results suggested that PUF9 interacts directly with SLBP2, but PUF9 was not detected in EIF4E2 pull-downs. We speculate that the EIF4E2-SLBP2 complex might interact with its target mRNAs, perhaps via PUF9, only early during G1/S, stabilizing the mRNAs in preparation for translation later in S-phase or in early G2.

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.

Stage-specific transcription activator ESB1 regulates monoallelic antigen expression in Trypanosoma brucei

Variant surface glycoprotein (VSG) coats bloodstream form Trypanosoma brucei parasites, and monoallelic VSG expression underpins the antigenic variation necessary for pathogenicity. One of thousands of VSG genes is transcribed by RNA polymerase I in a singular nuclear structure called the expression site body (ESB), but how monoallelic VSG transcription is achieved remains unclear. Using a localization screen of 153 proteins we found one, ESB-specific protein 1 (ESB1), that localized only to the ESB and is expressed only in VSG-expressing life cycle stages. ESB1 associates with DNA near the active VSG promoter and is necessary for VSG expression, with overexpression activating inactive VSG promoters. Mechanistically, ESB1 is necessary for recruitment of a subset of ESB components, including RNA polymerase I, revealing that the ESB has separately assembled subdomains. Because many trypanosomatid parasites have divergent ESB1 orthologues yet do not undergo antigenic variation, ESB1 probably represents an important class of transcription regulators.

Paving the Way: Contributions of Big Data to Apicomplexan and Kinetoplastid Research

In the age of big data an important question is how to ensure we make the most out of the resources we generate. In this review, we discuss the major methods used in Apicomplexan and Kinetoplastid research to produce big datasets and advance our understanding of Plasmodium, Toxoplasma, Cryptosporidium, Trypanosoma and Leishmania biology. We debate the benefits and limitations of the current technologies, and propose future advancements that may be key to improving our use of these techniques. Finally, we consider the difficulties the field faces when trying to make the most of the abundance of data that has already been, and will continue to be, generated.

Genomic Occupancy of the Bromodomain Protein Bdf3 Is Dynamic during Differentiation of African Trypanosomes from Bloodstream to Procyclic Forms

Trypanosoma brucei, the causative agent of human and animal African trypanosomiasis, cycles between a mammalian host and a tsetse fly vector. The parasite undergoes huge changes in morphology and metabolism during adaptation to each host environment. These changes are reflected in the different transcriptomes of parasites living in each host. However, it remains unclear whether chromatin-interacting proteins help mediate these changes. Bromodomain proteins localize to transcription start sites in bloodstream parasites, but whether the localization of bromodomain proteins changes as parasites differentiate from bloodstream to insect stages remains unknown. To address this question, we performed cleavage under target and release using nuclease (CUT&RUN) against bromodomain protein 3 (Bdf3) in parasites differentiating from bloodstream to insect forms. We found that Bdf3 occupancy at most loci increased at 3 h following onset of differentiation and decreased thereafter. A number of sites with increased bromodomain protein occupancy lie proximal to genes with altered transcript levels during differentiation, such as procyclins, procyclin-associated genes, and invariant surface glycoproteins. Most Bdf3-occupied sites are observed throughout differentiation. However, one site appears de novo during differentiation and lies proximal to the procyclin gene locus housing genes essential for remodeling surface proteins following transition to the insect stage. These studies indicate that occupancy of chromatin-interacting proteins is dynamic during life cycle stage transitions and provide the groundwork for future studies on the effects of changes in bromodomain protein occupancy. Additionally, the adaptation of CUT&RUN for Trypanosoma brucei provides other researchers with an alternative to chromatin immunoprecipitation (ChIP).

IMPORTANCE The parasite Trypanosoma brucei is the causative agent of human and animal African trypanosomiasis (sleeping sickness). Trypanosomiasis, which affects humans and cattle, is fatal if untreated. Existing drugs have significant side effects. Thus, these parasites impose a significant human and economic burden in sub-Saharan Africa, where trypanosomiasis is endemic. T. brucei cycles between the mammalian host and a tsetse fly vector, and parasites undergo huge changes in morphology and metabolism to adapt to different hosts. Here, we show that DNA-interacting bromodomain protein 3 (Bdf3) shows changes in occupancy at its binding sites as parasites transition from the bloodstream to the insect stage. Additionally, a new binding site appears near the locus responsible for remodeling of parasite surface proteins during transition to the insect stage. Understanding the mechanisms behind host adaptation is important for understanding the life cycle of the parasite.

Several different sequences are implicated in bloodstream-form-specific gene expression in Trypanosoma brucei

The parasite Trypanosoma brucei grows as bloodstream forms in mammalian hosts, and as procyclic forms in tsetse flies. In trypanosomes, gene expression regulation depends heavily on post-transcriptional mechanisms. Both the RNA-binding protein RBP10 and glycosomal phosphoglycerate kinase PGKC are expressed only in mammalian-infective forms. RBP10 targets procyclic-specific mRNAs for destruction, while PGKC is required for bloodstream-form glycolysis. Developmental regulation of both is essential: expression of either RBP10 or PGKC in procyclic forms inhibits their proliferation. We show that the 3’-untranslated region of the RBP10 mRNA is extraordinarily long—7.3kb—and were able to identify six different sequences, scattered across the untranslated region, which can independently cause bloodstream-form-specific expression. The 3’-untranslated region of the PGKC mRNA, although much shorter, still contains two different regions, of 125 and 153nt, that independently gave developmental regulation. No short consensus sequences were identified that were enriched either within these regulatory regions, or when compared with other mRNAs with similar regulation, suggesting that more than one regulatory RNA-binding protein is important for repression of mRNAs in procyclic forms. We also identified regions, including an AU repeat, that increased expression in bloodstream forms, or suppressed it in both forms. Trypanosome mRNAs that encode RNA-binding proteins often have extremely extended 3’-untranslated regions. We suggest that one function of this might be to act as a fail-safe mechanism to ensure correct regulation even if mRNA processing or expression of trans regulators is defective.

The parasite Trypanosoma brucei causes sleeping sickness in humans, and nagana in cattle, and is transmitted by Tsetse flies. It grows in the bloodstream and tissue fluids of mammalian hosts, as "bloodstream forms", and as "procyclic forms" in the midgut of tsetse flies. Several hundred proteins are expressed in a stage-specific fashion, and this is essential for parasite survival in the different environments. RBP10 is an RNA-binding protein that is expressed only in bloodstream forms. It binds to procyclic-specific mRNAs, and causes their destruction. PGKC is an enzyme that is also specifically expressed in bloodstream forms. Developmental regulation of both is essential: expression of either RBP10 or PGKC in procyclic forms prevents their growth. The mRNAs encoding both proteins are very unstable in procyclic forms, and the sequences responsible are in an "untranslated region" of the mRNA—sequences that follow the part that codes for protein. We here show that the mRNA encoding PGKC has two regions that independently cause developmental regulation, and that the very long untranslated region of the RBP10 mRNA has no fewer than six regulatory regions, but there were no obvious similarities between them. We suggest that the presence of several different regulatory sequences in trypanosome mRNAs might be a fail-safe mechanism to ensure correct regulation.

Polysomal mRNA Association and Gene Expression in Trypanosoma brucei

Background: The contrasting physiological environments of  Trypanosoma brucei procyclic (insect vector) and bloodstream (mammalian host) forms necessitates deployment of different molecular processes and, therefore, changes in protein expression. Transcriptional regulation is unusual in  T. brucei because the arrangement of genes is polycistronic; however, genes which are transcribed together are subsequently cleaved into separate mRNAs by  trans-splicing. Following pre-mRNA processing, the regulation of mature mRNA stability is a tightly controlled cellular process. While many stage-specific transcripts have been identified, previous studies using RNA-seq suggest that changes in overall transcript level do not necessarily reflect the abundance of the corresponding protein.

Methods: To better understand the regulation of gene expression in  T. brucei, we performed a bioinformatic analysis of RNA-seq on total, sub-polysomal, and polysomal mRNA samples. We further cross-referenced our dataset with a previously published proteomics dataset to identify new protein coding sequences.

Results: Our analyses showed that several long non-coding RNAs are more abundant in the sub-polysome samples, which possibly implicates them in regulating cellular differentiation in  T. brucei. We also improved the annotation of the  T.brucei genome by identifying new putative protein coding transcripts that were confirmed by mass spectrometry data.

Conclusions: Several long non-coding RNAs are more abundant in the sub-polysome cellular fractions and might pay a role in the regulation of gene expression. We hope that these data will be of wide general interest, as well as being of specific value to researchers studying gene regulation expression and life stage transitions in  T. brucei.

H2B.V demarcates divergent strand-switch regions, some tDNA loci, and genome compartments in Trypanosoma cruzi and affects parasite differentiation and host cell invasion

Histone variants play a crucial role in chromatin structure organization and gene expression. Trypanosomatids have an unusual H2B variant (H2B.V) that is known to dimerize with the variant H2A.Z generating unstable nucleosomes. Previously, we found that H2B.V protein is enriched in tissue-derived trypomastigote (TCT) life forms, a nonreplicative stage of Trypanosoma cruzi, suggesting that this variant may contribute to the differences in chromatin structure and global transcription rates observed among parasite life forms. Here, we performed the first genome-wide profiling of histone localization in T. cruzi using epimastigotes and TCT life forms, and we found that H2B.V was preferentially located at the edges of divergent transcriptional strand switch regions, which encompass putative transcriptional start regions; at some tDNA loci; and between the conserved and disrupted genome compartments, mainly at trans-sialidase, mucin and MASP genes. Remarkably, the chromatin of TCT forms was depleted of H2B.V-enriched peaks in comparison to epimastigote forms. Interactome assays indicated that H2B.V associated specifically with H2A.Z, bromodomain factor 2, nucleolar proteins and a histone chaperone, among others. Parasites expressing reduced H2B.V levels were associated with higher rates of parasite differentiation and mammalian cell infectivity. Taken together, H2B.V demarcates critical genomic regions and associates with regulatory chromatin proteins, suggesting a scenario wherein local chromatin structures associated with parasite differentiation and invasion are regulated during the parasite life cycle.

A profile of research on the parasitic trypanosomatids and the diseases they cause

The parasitic trypanosomatids cause lethal and debilitating diseases, the leishmaniases, Chagas disease, and the African trypanosomiases, with major impacts on human and animal health. Sustained research has borne fruit by assisting efforts to reduce the burden of disease and by improving our understanding of fundamental molecular and cell biology. But where has the research primarily been conducted, and which research areas have received the most attention? These questions are addressed below using publication and citation data from the past few decades.

A distinct complex of PRP19-related and trypanosomatid-specific proteins is required for pre-mRNA splicing in trypanosomes

The pre-mRNA splicing factor PRP19 is recruited into the spliceosome after forming the PRP19/CDC5L complex in humans and the Nineteen complex in yeast. Additionally, ‘PRP19-related’ proteins enter the spliceosome individually or in pre-assemblies that differ in these systems. The protistan family Trypanosomatidae, which harbors parasites such as Trypanosoma brucei, diverged early during evolution from opisthokonts. While introns are rare in these organisms, spliced leader trans splicing is an obligatory step in mRNA maturation. So far, ∼70 proteins have been identified as homologs of human and yeast splicing factors. Moreover, few proteins of unknown function have recurrently co-purified with splicing proteins. Here we silenced the gene of one of these proteins, termed PRC5, and found it to be essential for cell viability and pre-mRNA splicing. Purification of PRC5 combined with sucrose gradient sedimentation revealed a complex of PRC5 with a second trypanosomatid-specific protein, PRC3, and PRP19-related proteins SYF1, SYF3 and ISY1, which we named PRP19-related complex (PRC). Importantly, PRC and the previously described PRP19 complex are distinct from each other because PRC, unlike PRP19, co-precipitates U4 snRNA, which indicates that PRC enters the spliceosome prior to PRP19 and uncovers a unique pre-organization of these proteins in trypanosomes.

Interactions of the Trypanosoma brucei brucei zinc-finger-domain protein ZC3H28

In Trypanosoma brucei and related Kinetoplastids, regulation of gene expression occurs mostly post-transcriptionally, and RNA-binding proteins play a critical role in the regulation of mRNA and protein abundance. Trypanosoma brucei ZC3H28 is a 114 KDa cytoplasmic mRNA-binding protein with a single C(x)7C(x)5C(x)sH zinc finger at the C-terminus and numerous proline-, histidine- or glutamine-rich regions. ZC3H28 is essential for normal bloodstream-form trypanosome growth, and when tethered to a reporter mRNA, ZC3H28 increased reporter mRNA and protein levels. Purification of N-terminally tagged ZC3H28 followed by mass spectrometry showed enrichment of ribosomal proteins, various RNA-binding proteins including both poly(A) binding proteins, the translation initiation complex EIF4E4/EIF4G3, and the activator MKT1. Tagged ZC3H28 was preferentially associated with long RNAs that have low complexity sequences in their 3′-untranslated regions; their coding regions also have low ribosome densities. In agreement with the tethering results, after ZC3H28 depletion, the levels of a significant proportion of its bound mRNAs decreased. We suggest that ZC3H28 is implicated in the stabilization of long mRNAs that are poorly translated.

The establishment of variant surface glycoprotein monoallelic expression revealed by single-cell RNA-seq of Trypanosoma brucei in the tsetse fly salivary glands

The long and complex Trypanosoma brucei development in the tsetse fly vector culminates when parasites gain mammalian infectivity in the salivary glands. A key step in this process is the establishment of monoallelic variant surface glycoprotein (VSG) expression and the formation of the VSG coat. The establishment of VSG monoallelic expression is complex and poorly understood, due to the multiple parasite stages present in the salivary glands. Therefore, we sought to further our understanding of this phenomenon by performing single-cell RNA-sequencing (scRNA-seq) on these trypanosome populations. We were able to capture the developmental program of trypanosomes in the salivary glands, identifying populations of epimastigote, gamete, pre-metacyclic and metacyclic cells. Our results show that parasite metabolism is dramatically remodeled during development in the salivary glands, with a shift in transcript abundance from tricarboxylic acid metabolism to glycolytic metabolism. Analysis of VSG gene expression in pre-metacyclic and metacyclic cells revealed a dynamic VSG gene activation program. Strikingly, we found that pre-metacyclic cells contain transcripts from multiple VSG genes, which resolves to singular VSG gene expression in mature metacyclic cells. Single molecule RNA fluorescence in situ hybridisation (smRNA-FISH) of VSG gene expression following in vitro metacyclogenesis confirmed this finding. Our data demonstrate that multiple VSG genes are transcribed before a single gene is chosen. We propose a transcriptional race model governs the initiation of monoallelic expression.

African trypanosomes are parasitic protists which cause endemic disease in sub-Saharan Africa. To evade mammalian immune responses the parasite has developed a system of antigenic variation, where the surface of the cell is covered in a tightly packed coat of variant surface glycoproteins (VSGs). Each cell expresses only one variant surface glycoprotein at a time, and this is periodically switched to evade new antibodies. The process of singular gene expression is termed monoallelic expression and this has two components, establishment and maintenance, i.e. how a single gene is selected for expression and how its singular expression is maintained throughout successive generations. The establishment of monoallelic VSG gene expression occurs in the salivary gland of the tsetse fly vector, although this process is not well understood. We used single cell gene expression profiling applied to thousands of single cells in the salivary gland of the fly. We show that in order to select a single gene, trypanosomes initially transcribe multiple VSGs before a single gene is selected for high-level expression. We propose a model where this process is driven by a race to accumulate transcription factors at a single VSG gene.

Allele-specific assembly of a eukaryotic genome corrects apparent frameshifts and reveals a lack of nonsense-mediated mRNA decay

To date, most reference genomes represent a mosaic consensus sequence in which the homologous chromosomes are collapsed into one sequence. This approach produces sequence artefacts and impedes analyses of allele-specific mechanisms. Here, we report an allele-specific genome assembly of the diploid parasite Trypanosoma brucei and reveal allelic variants affecting gene expression. Using long-read sequencing and chromosome conformation capture data, we could assign 99.5% of all heterozygote variants to a specific homologous chromosome and build a 66 Mb long allele-specific genome assembly. The phasing of haplotypes allowed us to resolve hundreds of artefacts present in the previous mosaic consensus assembly. In addition, it revealed allelic recombination events, visible as regions of low allelic heterozygosity, enabling the lineage tracing of T. brucei isolates. Interestingly, analyses of transcriptome and translatome data of genes with allele-specific premature termination codons point to the absence of a nonsense-mediated decay mechanism in trypanosomes. Taken together, this study delivers a reference quality allele-specific genome assembly of T. brucei and demonstrates the importance of such assemblies for the study of gene expression control. We expect the new genome assembly will increase the awareness of allele-specific phenomena and provide a platform to investigate them.

Polysomal mRNA Association and Gene Expression in Trypanosoma brucei

Background: The contrasting physiological environments of  Trypanosoma brucei procyclic (insect vector) and bloodstream (mammalian host) forms necessitates deployment of different molecular processes and, therefore, changes in protein expression. Transcriptional regulation is unusual in  T. brucei because the arrangement of genes is polycistronic; however, genes which are transcribed together are subsequently cleaved into separate mRNAs by  trans-splicing. Following pre-mRNA processing, the regulation of mature mRNA stability is a tightly controlled cellular process. While many stage-specific transcripts have been identified, previous studies using RNA-seq suggest that changes in overall transcript level do not necessarily reflect the abundance of the corresponding protein. Methods: To better understand the regulation of gene expression in  T. brucei, we performed a bioinformatic analysis of RNA-seq on total, sub-polysomal, and polysomal mRNA samples. We further cross-referenced our dataset with a previously published proteomics dataset to identify new protein coding sequences. Results: Our analyses showed that several long non-coding RNAs are more abundant in the sub-polysome samples, which possibly implicates them in regulating cellular differentiation in  T. brucei. We also improved the annotation of the  T.brucei genome by identifying new putative protein coding transcripts that were confirmed by mass spectrometry data. Conclusions: Several long non-coding RNAs are more abundant in the sub-polysome cellular fractions and might pay a role in the regulation of gene expression. We hope that these data will be of wide general interest, as well as being of specific value to researchers studying gene regulation expression and life stage transitions in  T. brucei.

Application of single-cell transcriptomics to kinetoplastid research

Kinetoplastid parasites are responsible for both human and animal diseases across the globe where they have a great impact on health and economic well-being. Many species and life cycle stages are difficult to study due to limitations in isolation and culture, as well as to their existence as heterogeneous populations in hosts and vectors. Single-cell transcriptomics (scRNA-seq) has the capacity to overcome many of these difficulties, and can be leveraged to disentangle heterogeneous populations, highlight genes crucial for propagation through the life cycle, and enable detailed analysis of host–parasite interactions. Here, we provide a review of studies that have applied scRNA-seq to protozoan parasites so far. In addition, we provide an overview of sample preparation and technology choice considerations when planning scRNA-seq experiments, as well as challenges faced when analysing the large amounts of data generated. Finally, we highlight areas of kinetoplastid research that could benefit from scRNA-seq technologies.

In vitro culture of freshly isolated Trypanosoma brucei brucei bloodstream forms results in gene copy-number changes

Most researchers who study unicellular eukaryotes work with an extremely limited number of laboratory-adapted isolates that were obtained from the field decades ago, but the effects of passage in laboratory rodents, and adaptation to in vitro culture, have been little studied. For example, the vast majority of studies of Trypanosoma brucei biology have concentrated on just two strains, Lister 427 and EATRO1125, which were taken from the field over half a century ago and have since have undergone innumerable passages in rodents and culture. We here describe two new Trypanosoma brucei brucei strains. MAK65 and MAK98, which have undergone only 3 rodent passages since isolation from Ugandan cattle. High-coverage sequencing revealed that adaptation of the parasites to culture was accompanied by changes in gene copy numbers. T. brucei has so far been considered to be uniformly diploid, but we also found trisomy of chromosome 5 not only in one Lister 427 culture, but also in the MAK98 field isolate. Trisomy of chromosome 6, and increased copies of other chromosome segments, were also seen in established cultured lines. The two new T. brucei strains should be useful to researchers interested in trypanosome differentiation and pathogenicity. Initial results suggested that the two strains have differing infection patterns in rodents. MAK65 is uniformly diploid and grew more reproducibly in bloodstream-form culture than MAK98.

Polysomal mRNA Association and Gene Expression in Trypanosoma brucei

Background: The contrasting physiological environments of  Trypanosoma brucei procyclic (insect vector) and bloodstream (mammalian host) forms necessitates deployment of different molecular processes and, therefore, changes in protein expression. Transcriptional regulation is unusual in  T. brucei because the arrangement of genes is polycistronic; however, genes which are transcribed together are subsequently cleaved into separate mRNAs by  trans-splicing. Following pre-mRNA processing, the regulation of mature mRNA stability is a tightly controlled cellular process. While many stage-specific transcripts have been identified, previous studies using RNA-seq suggest that changes in overall transcript level do not necessarily reflect the abundance of the corresponding protein. Methods: To better understand the regulation of gene expression in  T. brucei, we performed a bioinformatic analysis of RNA-seq on total, sub-polysomal, and polysomal mRNA samples. We further cross-referenced our dataset with a previously published proteomics dataset to identify new protein coding sequences. Results: Our analyses showed that several long non-coding RNAs are more abundant in the sub-polysome samples, which possibly implicates them in regulating cellular differentiation in  T. brucei. We also improved the annotation of the  T.brucei genome by identifying new putative protein coding transcripts that were confirmed by mass spectrometry data. Conclusions: Several long non-coding RNAs are more abundant in the sub-polysome cellular fractions and might pay a role in the regulation of gene expression. We hope that these data will be of wide general interest, as well as being of specific value to researchers studying gene regulation expression and life stage transitions in  T. brucei.

Peculiarities of aminoacyl-tRNA synthetases from trypanosomatids

Trypanosomatid parasites are responsible for various human diseases, such as sleeping sickness, animal trypanosomiasis, or cutaneous and visceral leishmaniases. The few available drugs to fight related parasitic infections are often toxic and present poor efficiency and specificity, and thus, finding new molecular targets is imperative. Aminoacyl-tRNA synthetases (aaRSs) are essential components of the translational machinery as they catalyze the specific attachment of an amino acid onto cognate tRNA(s). In trypanosomatids, one gene encodes both cytosolic- and mitochondrial-targeted aaRSs, with only three exceptions. We identify here a unique specific feature of aaRSs from trypanosomatids, which is that most of them harbor distinct insertion and/or extension sequences. Among the 26 identified aaRSs in the trypanosome Leishmania tarentolae, 14 contain an additional domain or a terminal extension, confirmed in mature mRNAs by direct cDNA nanopore sequencing. Moreover, these RNA-Seq data led us to address the question of aaRS dual localization and to determine splice-site locations and the 5'-UTR lengths for each mature aaRS-encoding mRNA. Altogether, our results provided evidence for at least one specific mechanism responsible for mitochondrial addressing of some L. tarentolae aaRSs. We propose that these newly identified features of trypanosomatid aaRSs could be developed as relevant drug targets to combat the diseases caused by these parasites.

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.

The emerging paradigm of calcium homeostasis as a new therapeutic target for protozoan parasites

Calcium channels (CCs), a group of ubiquitously expressed membrane proteins, are involved in many pathophysiological processes of protozoan parasites. Our understanding of CCs in cell signaling, organelle function, cellular homeostasis, and cell cycle control has led to improved insights into their structure and functions. In this article, we discuss CCs characteristics of five major protozoan parasites Plasmodium, Leishmania, Toxoplasma, Trypanosoma, and Cryptosporidium. We provide a comprehensive review of current antiparasitic drugs and the potential of using CCs as new therapeutic targets. Interestingly, previous studies have demonstrated that human CC modulators can kill or sensitize parasites to antiparasitic drugs. Still, none of the parasite CCs, pumps, or transporters has been validated as drug targets. Information for this review draws from extensive data mining of genome sequences, chemical library screenings, and drug design studies. Parasitic resistance to currently approved therapeutics is a serious and emerging threat to both disease control and management efforts. In this article, we suggest that the disruption of calcium homeostasis may be an effective approach to develop new anti-parasite drug candidates and reduce parasite resistance.

Functional insights from a surface antigen mRNA-bound proteome

Trypanosoma brucei is the causative agent of human sleeping sickness. The parasites' variant surface glycoprotein (VSG) enables them to evade adaptive immunity via antigenic variation. VSG comprises 10% of total cell protein and the high stability of VSG mRNA is essential for trypanosome survival. To determine how VSG mRNA stability is maintained, we used mRNA affinity purification to identify all its associated proteins. CFB2 (cyclin F-box protein 2), an unconventional RNA-binding protein with an F-box domain, was specifically enriched with VSG mRNA. We demonstrate that CFB2 is essential for VSG mRNA stability, describe cis acting elements within the VSG 3'-untranslated region that regulate the interaction, identify trans-acting factors that are present in the VSG messenger ribonucleoprotein particle, and mechanistically explain how CFB2 stabilizes the mRNA of this key pathogenicity factor. Beyond T. brucei, the mRNP purification approach has the potential to supply detailed biological insight into metabolism of relatively abundant mRNAs in any eukaryote.

Ubiquitination and the Proteasome as Drug Targets in Trypanosomatid Diseases

The eukaryotic pathogens Trypanosoma brucei, Trypanosoma cruzi and Leishmania are responsible for debilitating diseases that affect millions of people worldwide. The numbers of drugs available to treat these diseases, Human African Trypanosomiasis, Chagas' disease and Leishmaniasis are very limited and existing treatments have substantial shortcomings in delivery method, efficacy and safety. The identification and validation of novel drug targets opens up new opportunities for the discovery of therapeutic drugs with better efficacy and safety profiles. Here, the potential of targeting the ubiquitin-proteasome system in these parasites is reviewed. Ubiquitination is the posttranslational attachment of one or more ubiquitin proteins to substrates, an essential eukaryotic mechanism that regulates a wide variety of cellular processes in many different ways. The best studied of these is the delivery of ubiquitinated substrates for degradation to the proteasome, the major cellular protease. However, ubiquitination can also regulate substrates in proteasome-independent ways, and proteasomes can degrade proteins to some extent in ubiquitin-independent ways. Because of these widespread roles, both ubiquitination and proteasomal degradation are essential for the viability of eukaryotes and the proteins that mediate these processes are therefore attractive drug targets in trypanosomatids. Here, the current understanding of these processes in trypanosomatids is reviewed. Furthermore, significant recent progress in the development of trypanosomatid-selective proteasome inhibitors that cure mouse models of trypanosomatid infections is presented. In addition, the targeting of the key enzyme in ubiquitination, the ubiquitin E1 UBA1, is discussed as an alternative strategy. Important differences between human and trypanosomatid UBA1s in susceptibility to inhibitors predicts that the selective targeting of these enzymes in trypanosomatids may also be feasible. Finally, it is proposed that activating enzymes of the ubiquitin-like proteins SUMO and NEDD8 may represent drug targets in these trypanosomatids as well.

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.

Nucleosome landscape reflects phenotypic differences in Trypanosoma cruzi life forms

Trypanosoma cruzi alternates between replicative and nonreplicative life forms, accompanied by a shift in global transcription levels and by changes in the nuclear architecture, the chromatin proteome and histone posttranslational modifications. To gain further insights into the epigenetic regulation that accompanies life form changes, we performed genome-wide high-resolution nucleosome mapping using two T. cruzi life forms (epimastigotes and cellular trypomastigotes). By combining a powerful pipeline that allowed us to faithfully compare nucleosome positioning and occupancy, more than 125 thousand nucleosomes were mapped, and approximately 20% of them differed between replicative and nonreplicative forms. The nonreplicative forms have less dynamic nucleosomes, possibly reflecting their lower global transcription levels and DNA replication arrest. However, dynamic nucleosomes are enriched at nonreplicative regulatory transcription initiation regions and at multigenic family members, which are associated with infective-stage and virulence factors. Strikingly, dynamic nucleosome regions are associated with GO terms related to nuclear division, translation, gene regulation and metabolism and, notably, associated with transcripts with different expression levels among life forms. Finally, the nucleosome landscape reflects the steady-state transcription expression: more abundant genes have a more deeply nucleosome-depleted region at putative 5' splice sites, likely associated with trans-splicing efficiency. Taken together, our results indicate that chromatin architecture, defined primarily by nucleosome positioning and occupancy, reflects the phenotypic differences found among T. cruzi life forms despite the lack of a canonical transcriptional control context.

Nutrient availability regulates proline/alanine transporters in Trypanosoma brucei

Trypanosoma brucei is a species of unicellular parasite that can cause severe diseases in livestock and humans, including African trypanosomiasis and Chagas disease. Adaptation to diverse environments and changes in nutritional conditions is essential for T. brucei to establish an infection when changing hosts or during invasion of different host tissues. One such adaptation is the ability of T. brucei to rapidly switch its energy metabolism from glucose metabolism in the mammalian blood to proline catabolism in the insect stages and vice versa. However, the mechanisms that support the parasite's response to nutrient availability remain unclear. Using RNAseq and qRT-PCR, we investigated the response of T. brucei to amino acid or glucose starvation and found increased mRNA levels of several amino acid transporters, including all genes of the amino acid transporter AAT7-B subgroup. Functional characterization revealed that AAT7-B members are plasma membrane-localized in T. brucei and when expressed in Saccharomyces cerevisiae supported the uptake of proline, alanine, and cysteine, while other amino acids were poorly recognized. All AAT7-B members showed a preference for proline, which is transported with high or low affinity. RNAi-mediated AAT7-B downregulation resulted in a reduction of intracellular proline concentrations and growth arrest under low proline availability in cultured procyclic form parasites. Taken together, these results suggest a role of AAT7-B transporters in the response of T. brucei to proline starvation and proline catabolism.

The influence of soil drought stress on the leaf transcriptome of faba bean (Vicia faba L.) in the Qinghai-Tibet Plateau

Water deficit has a significant impact on growth, development and yield of fava bean (Vicia fava L.) in arid and semi-arid climates. The aim of this study was to identify differentially expressed genes in the Qinghai 13 genotype under soil drought through leaf transcriptome analysis. A total of 256.95 M clean reads were obtained and assembled into 176334 unigenes, with an average length of 766 bp. A total of 9126 (4439 upregulated and 4687 downregulated) differentially expressed genes (DEGs) were identified in faba bean leaves under soil drought. In total, 324 putative transcription factors were identified and classified as belonging to different transcription factor families. According to GO and KEGG analysis, the soil drought stress-inducible DEGs encoded proteins mainly involved in regulating photosynthesis, osmotic adjustment, detoxification, autophagy and other functions. In addition, a large portion of DEGs appeared to be novel because they could not be annotated in any functional databases, therefore, suggesting a specific response to soil drought in faba bean. Finally, RNA-seq analysis was validated by quantitative reverse-transcription PCR analysis. This work provides comprehensive and valuable information for understanding the molecular mechanisms which faba bean uses to respond to soil drought.

Positively selected modifications in the pore of TbAQP2 allow pentamidine to enter Trypanosoma brucei

Mutations in the Trypanosoma brucei aquaporin AQP2 are associated with resistance to pentamidine and melarsoprol. We show that TbAQP2 but not TbAQP3 was positively selected for increased pore size from a common ancestor aquaporin. We demonstrate that TbAQP2’s unique architecture permits pentamidine permeation through its central pore and show how specific mutations in highly conserved motifs affect drug permeation. Introduction of key TbAQP2 amino acids into TbAQP3 renders the latter permeable to pentamidine. Molecular dynamics demonstrates that permeation by dicationic pentamidine is energetically favourable in TbAQP2, driven by the membrane potential, although aquaporins are normally strictly impermeable for ionic species. We also identify the structural determinants that make pentamidine a permeant although most other diamidine drugs are excluded. Our results have wide-ranging implications for optimising antitrypanosomal drugs and averting cross-resistance. Moreover, these new insights in aquaporin permeation may allow the pharmacological exploitation of other members of this ubiquitous gene family.

African sleeping sickness is a potentially deadly illness caused by the parasite Trypanosoma brucei. The disease is treatable, but many of the current treatments are old and are becoming increasingly ineffective. For instance, resistance is growing against pentamidine, a drug used in the early stages in the disease, as well as against melarsoprol, which is deployed when the infection has progressed to the brain. Usually, cases resistant to pentamidine are also resistant to melarsoprol, but it is still unclear why, as the drugs are chemically unrelated.

Studies have shown that changes in a water channel called aquaglyceroporin 2 (TbAQP2) contribute to drug resistance in African sleeping sickness; this suggests that it plays a role in allowing drugs to kill the parasite. This molecular ‘drain pipe’ extends through the surface of T. brucei, and should allow only water and a molecule called glycerol in and out of the cell. In particular, the channel should be too narrow to allow pentamidine or melarsoprol to pass through. One possibility is that, in T. brucei, the TbAQP2 channel is abnormally wide compared to other members of its family. Alternatively, pentamidine and melarsoprol may only bind to TbAQP2, and then ‘hitch a ride’ when the protein is taken into the parasite as part of the natural cycle of surface protein replacement.

Alghamdi et al. aimed to tease out these hypotheses. Computer models of the structure of the protein were paired with engineered changes in the key areas of the channel to show that, in T. brucei, TbAQP2 provides a much broader gateway into the cell than observed for similar proteins. In addition, genetic analysis showed that this version of TbAQP2 has been actively selected for during the evolution process of T. brucei. This suggests that the parasite somehow benefits from this wider aquaglyceroporin variant.

This is a new resistance mechanism, and it is possible that aquaglyceroporins are also larger than expected in other infectious microbes. The work by Alghamdi et al. therefore provides insight into how other germs may become resistant to drugs.

A potential role for a novel ZC3H5 complex in regulating mRNA translation in Trypanosoma brucei

In Trypanosoma brucei and related kinetoplastids, gene expression regulation occurs mostly posttranscriptionally. Consequently, RNA-binding proteins play a critical role in the regulation of mRNA and protein abundance. Yet, the roles of many RNA-binding proteins are not understood. Our previous research identified the RNA-binding protein ZC3H5 as possibly involved in gene repression, but its role in controlling gene expression was unknown. We here show that ZC3H5 is an essential cytoplasmic RNA-binding protein. RNAi targeting ZC3H5 causes accumulation of precytokinetic cells followed by rapid cell death. Affinity purification and pairwise yeast two-hybrid analysis suggest that ZC3H5 forms a complex with three other proteins, encoded by genes Tb927.11.4900, Tb927.8.1500, and Tb927.7.3040. RNA immunoprecipitation revealed that ZC3H5 is preferentially associated with poorly translated, low-stability mRNAs, the 5'-untranslated regions and coding regions of which are enriched in the motif (U/A)UAG(U/A). As previously found in high-throughput analyses, artificial tethering of ZC3H5 to a reporter mRNA or other complex components repressed reporter expression. However, depletion of ZC3H5 in vivo caused only very minor decreases in a few targets, marked increases in the abundances of very stable mRNAs, an increase in monosomes at the expense of large polysomes, and appearance of "halfmer" disomes containing two 80S subunits and one 40S subunit. We speculate that the ZC3H5 complex might be implicated in quality control during the translation of suboptimal open reading frames.

The RNA-associated proteins MKT1 and MKT1L form alternative PBP1-containing complexes in Trypanosoma brucei

Control of gene expression in kinetoplastids such as trypanosomes depends heavily on RNA-binding proteins that influence mRNA decay and translation. We previously showed that the trypanosome protein MKT1 forms a multicomponent protein complex: MKT1 interacts with PBP1, which in turn recruits LSM12 and poly(A)-binding protein. MKT1 is recruited to mRNAs by sequence-specific RNA-binding proteins, resulting in stabilization of the bound mRNA. We here show that PBP1, LSM12, and a 117-residue protein, XAC1 (Tb927.7.2780), are present in complexes that contain either MKT1 or an MKT1-like protein, MKT1L (Tb927.10.1490). All five proteins are present predominantly in the complexes, and we found evidence for a minor subset of complexes containing both MKT1 and MKT1L. XAC1-containing complexes reproducibly contained RNA-binding proteins that were previously found associated with MKT1. Moreover, XAC1- or MKT1-containing complexes specifically recruited one of the two poly(A)-binding proteins, PABP2, and one of the six cap-binding translation initiation complexes, EIF4E6-EIF4G5. Yeast two-hybrid assay results indicated that MKT1 directly interacts with EIF4G5. MKT1-PBP1 complexes can therefore interact with mRNAs via their poly(A) tails and caps, as well as through sequence-specific RNA-binding proteins. Correspondingly, MKT1 is associated with many mRNAs, although not with those encoding ribosomal proteins. Meanwhile, MKT1L resembles MKT1 at the C terminus but additionally features an N-terminal extension with low-complexity regions. Although MKT1L depletion inhibited cell proliferation, we found no evidence that it specifically interacts with RNA-binding proteins or mRNA. We speculate that MKT1L may compete with MKT1 for PBP1 binding and thereby modulate the function of MKT1-containing complexes.

Roles of the Pumilio domain protein PUF3 in Trypanosoma brucei growth and differentiation

Trypanosomes strongly rely on post-transcriptional mechanisms to control gene expression. Several Opisthokont Pumilio domain proteins are known to suppress expression when bound to mRNAs. The Trypanosoma brucei Pumilio domain protein PUF3 is a cytosolic mRNA-binding protein that suppresses expression when tethered to a reporter mRNA. RNA-binding studies showed that PUF3 preferentially binds to mRNAs with a classical Pumilio-domain recognition motif, UGUA[U/C]AUU. RNA-interference-mediated reduction of PUF3 in bloodstream forms caused a minor growth defect, but the transcriptome was not affected. Depletion of PUF3 also slightly delayed differentiation to the procyclic form. However, both PUF3 genes could be deleted in cultured bloodstream- and procyclic-form trypanosomes. Procyclic forms without PUF3 also grew somewhat slower than wild-type, but ectopic expression of C-terminally tagged PUF3 impaired their viability. PUF3 was not required for RBP10-induced differentiation of procyclic forms to bloodstream forms. Mass spectrometry revealed no PUF3 binding partners that might explain its suppressive activity. We conclude that PUF3 may have a role in fine-tuning gene expression. Since PUF3 is conserved in all Kinetoplastids, including those that do not infect vertebrates, we suggest that it might confer advantages within the invertebrate host.

The mitochondrial phosphate carrier TbMCP11 is essential for mitochondrial function in the procyclic form of Trypanosoma brucei

Conserved amongst all eukaryotes is a family of mitochondrial carrier proteins (SLC25A) responsible for the import of various solutes across the inner mitochondrial membrane. We previously reported that the human parasite Trypanosoma brucei possesses 26 SLC25A proteins (TbMCPs) amongst which two, TbMCP11 and TbMCP8, were predicted to function as phosphate importers. The transport of inorganic phosphate into the mitochondrion is a prerequisite to drive ATP synthesis by substrate level and oxidative phosphorylation and thus crucial for cell viability. In this paper we describe the functional characterization of TbMCP11. In procyclic form T. brucei, the RNAi of TbMCP11 blocked ATP synthesis on mitochondrial substrates, caused a drop of the mitochondrial oxygen consumption and drastically reduced cell viability. The functional complementation in yeast and mitochondrial swelling experiments suggested a role for TbMCP11 as inorganic phosphate carrier. Interestingly, procyclic form T. brucei cells in which TbMCP11 was depleted displayed an inability to either replicate or divide the kinetoplast DNA, which resulted in a severe cytokinesis defect.

A Leap Into the Unknown - Early Events in African Trypanosome Transmission

African trypanosomes are mainly transmitted by tsetse flies. In recent years there has been good progress in understanding how the parasites prepare for transmission, detect their changed environment through the perception of different environmental cues, and respond by changing their developmental gene expression. In this review, we discuss the different signals and signaling mechanisms used by the parasites to carry out the early events necessary for their establishment in the fly. We also compare Trypanosoma brucei and Trypanosoma congolense, parasites that share a common pathway in the early stages of fly colonization but apparently use different mechanisms to achieve this.

Structural Basis of Protein Arginine Methyltransferase Activation by a Catalytically Dead Homolog (Prozyme)

Prozymes are pseudoenzymes that stimulate the function of weakly active enzymes through complex formation. The major Trypanosoma brucei protein arginine methyltransferase, TbPRMT1 enzyme (ENZ), requires TbPRMT1 prozyme (PRO) to form an active heterotetrameric complex. Here, we present the X-ray crystal structure of the TbPRMT1 ENZ-Δ52PRO tetrameric complex with the cofactor product S-adenosyl-l-homocysteine (AdoHcy) at 2.4 Å resolution. The individual ENZ and PRO units adopt the highly-conserved PRMT domain architecture and form an antiparallel heterodimer that corresponds to the canonical homodimer observed in all previously reported PRMTs. In turn, two such heterodimers assemble into a tetramer both in the crystal and in solution with twofold rotational symmetry. ENZ is unstable in absence of PRO and incapable of forming a homodimer due to a steric clash of an ENZ-specific tyrosine within the dimerization arm, rationalizing why PRO is required to complement ENZ to form a PRMT dimer that is necessary, but not sufficient for PRMT activity. The PRO structure deviates from other, active PRMTs in that it lacks the conserved η2 3₁₀-helix within the Rossmann fold, abolishing cofactor binding. In addition to its chaperone function for ENZ, PRO substantially contributes to substrate binding. Heterotetramerization is required for catalysis, as heterodimeric ENZ-PRO mutants lack binding affinity and methyltransferase activity toward the substrate protein TbRGG1. Together, we provide a structural basis for TbPRMT1 ENZ activation by PRO heterotetramer formation, which is conserved across all kinetoplastids, and describe a chaperone function of the TbPRMT1 prozyme, which represents a novel mode of PRMT regulation.

TDR Targets 6: driving drug discovery for human pathogens through intensive chemogenomic data integration

The volume of biological, chemical and functional data deposited in the public domain is growing rapidly, thanks to next generation sequencing and highly-automated screening technologies. These datasets represent invaluable resources for drug discovery, particularly for less studied neglected disease pathogens. To leverage these datasets, smart and intensive data integration is required to guide computational inferences across diverse organisms. The TDR Targets chemogenomics resource integrates genomic data from human pathogens and model organisms along with information on bioactive compounds and their annotated activities. This report highlights the latest updates on the available data and functionality in TDR Targets 6. Based on chemogenomic network models providing links between inhibitors and targets, the database now incorporates network-driven target prioritizations, and novel visualizations of network subgraphs displaying chemical- and target-similarity neighborhoods along with associated target-compound bioactivity links. Available data can be browsed and queried through a new user interface, that allow users to perform prioritizations of protein targets and chemical inhibitors. As such, TDR Targets now facilitates the investigation of drug repurposing against pathogen targets, which can potentially help in identifying candidate targets for bioactive compounds with previously unknown targets. TDR Targets is available at https://tdrtargets.org.

High-Throughput Sequencing for Trypanosome Transcriptome Characterization

High-throughput sequencing of cDNA (RNASeq) is now the method of choice for analysis of transcriptomes. This chapter details important considerations in the design of RNASeq experiments for kinetoplastids grown in culture or experimental animals. It contains protocols for obtaining parasites from rodents, and for removal of rRNA from total RNA. In addition, custom pipelines for sequence alignment, and for data analysis and visualization, are described.

The zinc finger proteins ZC3H20 and ZC3H21 stabilise mRNAs encoding membrane proteins and mitochondrial proteins in insect-form Trypanosoma brucei

ZC3H20 and ZC3H21 are related trypanosome proteins with two C(x)₈ C(x)₅ C(x)₃ H zinc finger motifs. ZC3H20 is present at a low level in replicating mammalian-infective bloodstream forms, but becomes more abundant when they undergo growth arrest at high density; ZC3H21 appears only in the procyclic form of the parasite, which infects Tsetse flies. Each protein binds to several hundred mRNAs, with overlapping but not identical specificities. Both increase expression of bound mRNAs, probably through recruitment of the MKT1-PBP1 complex. At least 28 of the bound mRNAs decrease after depletion of ZC3H20, or of ZC3H20 and ZC3H21 together; their products include procyclic-specific proteins of the plasma membrane and energy metabolism. Simultaneous depletion of ZC3H20 and ZC3H21 causes procyclic forms to shrink and stop growing; in addition to decreases in target mRNAs, there are other changes suggestive of loss of developmental regulation. The bloodstream-form-specific protein RBP10 controls ZC3H20 and ZC3H21 expression. Interestingly, some ZC3H20/21 target mRNAs also bind to and are repressed by RBP10, allowing for dynamic regulation as RBP10 decreases and ZC3H20 and ZC3H21 increase during differentiation.

P-type transport ATPases in Leishmania and Trypanosoma

P-type ATPases are critical to the maintenance and regulation of cellular ion homeostasis and membrane lipid asymmetry due to their ability to move ions and phospholipids against a concentration gradient by utilizing the energy of ATP hydrolysis. P-type ATPases are particularly relevant in human pathogenic trypanosomatids which are exposed to abrupt and dramatic changes in their external environment during their life cycles. This review describes the complete inventory of ion-motive, P-type ATPase genes in the human pathogenic Trypanosomatidae; eight Leishmania species (L. aethiopica, L. braziliensis, L. donovani, L. infantum, L. major, L. mexicana, L. panamensis, L. tropica), Trypanosoma cruzi and three Trypanosoma brucei subspecies (Trypanosoma brucei brucei TREU927, Trypanosoma brucei Lister strain 427, Trypanosoma brucei gambiense DAL972). The P-type ATPase complement in these trypanosomatids includes the P_1B (metal pumps), P_2A (SERCA, sarcoplasmic-endoplasmic reticulum calcium ATPases), P_2B (PMCA, plasma membrane calcium ATPases), P_2D (Na⁺ pumps), P_3A (H⁺ pumps), P₄ (aminophospholipid translocators), and P_5B (no assigned specificity) subfamilies. These subfamilies represent the P-type ATPase transport functions necessary for survival in the Trypanosomatidae as P-type ATPases for each of these seven subfamilies are found in all Leishmania and Trypanosoma species included in this analysis. These P-type ATPase subfamilies are correlated with current molecular and biochemical knowledge of their function in trypanosomatid growth, adaptation, infectivity, and survival.