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

Leishmania Ribosomal Protein (RP) paralogous genes compensate each other's expression maintaining protein native levels

In the protozoan parasite Leishmania, most genes encoding for ribosomal proteins (RPs) are present as two or more copies in the genome. However, their untranslated regions (UTRs) are predominantly divergent and might be associated with a distinct regulation of the expression of paralogous genes. Herein, we investigated the expression profiles of two RPs (S16 and L13a) encoded by duplicated genes in Leishmania major. The genes encoding for the S16 protein possess identical coding sequences (CDSs) and divergent UTRs, whereas the CDSs of L13a diverge by two amino acids and by their UTRs. Using CRISPR/Cas9 genome editing, we generated knockout (Δ) and endogenously tagged transfectants for each paralog of L13a and S16 genes. Combining tagged and Δ cell lines we found evidence of differential expression of both RPS16 and RPL13a isoforms throughout parasite development, with one isoform consistently more abundant than its respective copy. In addition, compensatory expression was observed for each paralog upon deletion of the corresponding isoform, suggesting functional conservation between these proteins. This differential expression pattern relates to post-translational processes, given compensation occurs at the level of the protein, with no alterations detected at transcript level. Ribosomal profiles for RPL13a indicate a standard behavior for these paralogues suggestive of interaction with heavy RNA-protein complexes, as already reported for other RPs in trypanosomatids. We identified paralog-specific bound to their 3'UTRs which may be influential in regulating paralog expression. In support, we identified conserved cis-elements within the 3'UTRs of RPS16 and RPL13a; cis-elements exclusive to the UTR of the more abundant paralog or to the less abundant ones were identified.

A sticky situation: When trypanosomatids attach to insect tissues

Transmission of trypanosomatids to their mammalian hosts requires a complex series of developmental transitions in their insect vectors, including stable attachment to an insect tissue. While there are many ultrastructural descriptions of attached cells, we know little about the signaling events and molecular mechanisms involved in this process. Each trypanosomatid species attaches to a specific tissue in the insect at a particular stage of its life cycle. Attachment is mediated by the flagellum, which is modified to accommodate a filament-rich plaque within an expanded region of the flagellar membrane. Attachment immediately precedes differentiation to the mammal-infectious stage and in some cases a direct mechanistic link has been demonstrated. In this review, we summarize the current state of knowledge of trypanosomatid attachment in insects, including structure, function, signaling, candidate molecules, and changes in gene expression. We also highlight remaining questions about this process and how the field is poised to address them through modern approaches.

Leishmania infantum (JPCM5) Transcriptome, Gene Models and Resources for an Active Curation of Gene Annotations

Leishmania infantum is one of the causative agents of visceral leishmaniases, the most severe form of leishmaniasis. An improved assembly for the L. infantum genome was published five years ago, yet delineation of its transcriptome remained to be accomplished. In this work, the transcriptome annotation was attained by a combination of both short and long RNA-seq reads. The good agreement between the results derived from both methodologies confirmed that transcript assembly based on Illumina RNA-seq and further delimitation according to the positions of spliced leader (SAS) and poly-A (PAS) addition sites is an adequate strategy to annotate the transcriptomes of Leishmania, a procedure previously used for transcriptome annotation in other Leishmania species and related trypanosomatids. These analyses also confirmed that the Leishmania transcripts boundaries are relatively slippery, showing extensive heterogeneity at the 5′- and 3′-ends. However, the use of RNA-seq reads derived from the PacBio technology (referred to as Iso-Seq) allowed the authors to uncover some complex transcription patterns occurring at particular loci that would be unnoticed by the use of short RNA-seq reads alone. Thus, Iso-Seq analysis provided evidence that transcript processing at particular loci would be more dynamic than expected. Another noticeable finding was the observation of a case of allelic heterozygosity based on the existence of chimeric Iso-Seq reads that might be generated by an event of intrachromosomal recombination. In addition, we are providing the L. infantum gene models, including both UTRs and CDS regions, that would be helpful for undertaking whole-genome expression studies. Moreover, we have built the foundations of a communal database for the active curation of both gene/transcript models and functional annotations for genes and proteins.

Role of the RNA-binding protein ZC3H41 in the regulation of ribosomal protein messenger RNAs in trypanosomes

BACKGROUND

Trypanosomes are single-celled eukaryotes that rely heavily on post-transcriptional mechanisms to regulate gene expression. RNA-binding proteins play essential roles in regulating the fate, abundance and translation of messenger RNAs (mRNAs). Among these, zinc finger proteins of the cysteine3histidine (CCCH) class have been shown to be key players in cellular processes as diverse as differentiation, regulation of the cell cycle and translation. ZC3H41 is an essential zinc finger protein that has been described as a component of spliced leader RNA granules and nutritional stress granules, but its role in RNA metabolism is unknown.

METHODS

Cell cycle analysis in ZC3H41- and Z41AP-depleted cells was carried out using 4',6-diamidino-2-phenylindole staining, microscopic examination and flow cytometry. The identification of ZC3H41 protein partners was done using tandem affinity purification and mass spectrometry. Next-generation sequencing was used to evaluate the effect of ZC3H41 depletion on the transcriptome of procyclic Trypanosoma brucei cells, and also to identify the cohort of mRNAs associated with the ZC3H41/Z41AP complex. Levels of 5S ribosomal RNA (rRNA) species in ZC3H41- and Z41AP-depleted cells were assessed by quantitative reverse transcription-polymerase chain reaction. Surface sensing of translation assays were used to monitor global translation.

RESULTS

We showed that depletion of the zinc finger protein ZC3H41 resulted in marked cell cycle defects and abnormal cell morphologies. ZC3H41 was found associated with an essential protein, which we named Z41AP, forming a stable heterodimer, and also with proteins of the poly(A)-binding protein 1 complex. The identification of mRNAs associated with the ZC3H41/Z41AP complex revealed that it is primarily composed of ribosomal protein mRNAs, and that binding to target transcripts is diminished upon nutritional stress. In addition, we observed that mRNAs encoding several proteins involved in the maturation of 5S rRNA are also associated with the ZC3H41/Z41AP complex. Finally, we showed that depletion of either ZC3H41 or Z41AP led to the accumulation of 5S rRNA precursors and a decrease of protein translation.

CONCLUSIONS

We propose that ZC3H41 and Z41AP play important roles in controlling the fate of ribosomal components in response to environmental cues.

A novel deep learning-assisted hybrid network for plasmodium falciparum parasite mitochondrial proteins classification

Plasmodium falciparum is a parasitic protozoan that can cause malaria, which is a deadly disease. Therefore, the accurate identification of malaria parasite mitochondrial proteins is essential for understanding their functions and identifying novel drug targets. For classifying protein sequences, several adaptive statistical techniques have been devised. Despite significant gains, prediction performance is still constrained by the lack of appropriate feature descriptors and learning strategies in current systems. Moreover, good ground truth data is important for Artificial Intelligence (AI)-based models but there is a lack of that data in the literature. Therefore, in this work, we propose a novel hybrid network that combines 1D Convolutional Neural Network (CNN) and Bidirectional Gated Recurrent Unit (BGRU) to classify the malaria parasite mitochondrial proteins. Furthermore, we curate a sequential data that are collected from National Center for Biotechnology Information (NCBI) and UniProtKB/Swiss-Prot proteins databanks to prepare a dataset that can be used by the research community for AI-based algorithms evaluation. We obtain 4204 cases after preprocessing of the collected data and denote this set of proteins as PF4204. Finally, we conduct an ablation study on several conventional and deep models using PF4204 and the benchmark PF2095 datasets. The proposed model 'CNN-BGRU' obtains the accuracy values of 0.9096 and 0.9857 on PF4204 and PF2095 datasets, respectively. In addition, the CNN-BGRU is compared with state-of-the-arts, where the results illustrate that it can extract robust features and identify proteins accurately.

Multimodal regulation of encystation in Giardia duodenalis revealed by deep proteomics

Cyst formation in the parasitic protist Giardia duodenalis is critical to its transmission. Existing proteomic data quantifies only 17% of coding genes transcribed during encystation and does not cover the complete process from trophozoite to mature cyst. Using high-resolution mass spectrometry, we have quantified proteomic changes across encystation and compared this with published transcriptomic data. We reproducibly identified 3863 (64.5% of Giardia proteins) and quantified 3382 proteins (56.5% of Giardia proteins) over standard trophozoite growth (TY), during low-bile encystation priming (LB), 16 h into encystation (EC), and at cyst maturation (C). This work provides the first known expanded observation of encystation at the proteomic level and triples the coverage of previous encystation proteomes. One-third (1169 proteins) of the quantified proteome is differentially expressed in the mature cyst relative to the trophozoite, including proteasomal machinery, metabolic pathways, and secretory proteins. Changes in lipid metabolism indicated a shift in lipid species dependency during encystation. Consistent with this, we identified the first, putative lipid transporters in this species, representing the steroidogenic acute regulatory protein-related lipid transfer (StARkin), oxysterol binding protein related protein (ORP/Osh) and glycosphingolipid transfer protein (GLTP) families, and follow their differential expression over cyst formation. Lastly, we undertook correlation analyses of the transcriptome and proteome of trophozoites and cysts, and found evidence of post-transcriptional regulation of key protein classes (RNA binding proteins) and stage-specific genes (encystation markers) implicating translation-repression in encystation. We provide the most extensive proteomic analysis of encystation in Giardia to date and the first known exploration across its complete duration. This work identifies encystation as highly coordinated, involving major changes in proteostasis, metabolism and membrane dynamics, and indicates a potential role for post-transcriptional regulation, mediated through RNA-binding proteins. Together our work provides a valuable resource for Giardia research and the development of transmission-blocking anti-giardials.

An RNA-binding protein complex regulates the purine-dependent expression of a nucleobase transporter in trypanosomes

Post-transcriptional regulation of gene expression is particularly important in trypanosomatid protozoa. RNA-binding proteins (RBPs) regulate mRNA stability and translation, yet information about how RBPs are able to link environmental cues to post-transcriptional control is scarce. In Trypanosoma brucei, we have previously characterized a short RNA stem-loop cis-element (PuRE, Purine Responsive Element) within the 3'-UTR of the NT8 nucleobase transporter mRNA that is necessary and sufficient to confer a strong repression of gene expression in response to purines. In this study, we have identified a protein complex composed of two RNA-binding proteins (PuREBP1 and PuREBP2) that binds to the PuRE in vitro and to NT8 mRNA in vivo. Depletion of PuREBP1 by RNA interference results in the upregulation of just NT8 and the mRNAs encoding the amino acid transporter AATP6 paralogues. Moreover, we found that the PuREBP1/2 complex is associated with only a handful of mRNAs, and that it is responsible for the observed purine-dependent regulation of NT8 expression.

The translation initiation factor EIF4E5 from Leishmania: crystal structure and interacting partners

The mRNA cap-binding protein, eIF4E, mediates the recognition of the mRNA 5' end and, as part of the heterotrimeric eIF4F complex, facilitates the recruitment of the ribosomal subunits to initiate eukaryotic translation. Various regulatory events involving eIF4E and a second eIF4F subunit, eIF4G, are required for proper control of translation initiation. In pathogenic trypanosomatids, six eIF4Es and five eIF4Gs have been described, several forming different eIF4F-like complexes with yet unresolved roles. EIF4E5 is one of the least known of the trypanosomatid eIF4Es and has not been characterized in Leishmania species. Here, we used immunoprecipitation assays, combined with mass-spectrometry, to identify major EIF4E5 interacting proteins in L. infantum. A constitutively expressed, HA-tagged, EIF4E5 co-precipitated mainly with EIF4G1 and binding partners previously described in Trypanosoma brucei, EIF4G1-IP, RBP43 and the 14-3-3 proteins. In contrast, no clear co-precipitation with EIF4G2, also previously reported, was observed. EIF4E5 also co-precipitated with protein kinases, possibly associated with cell-cycle regulation, selected RNA binding proteins and histones. Phosphorylated residues were identified and mapped to the Leishmania-specific C-terminal end. Mutagenesis of the tryptophan residue (W53) postulated to mediate interactions with protein partners or of a neighbouring tryptophan conserved in Leishmania (W45) did not substantially impair the identified interactions. Finally, the crystal structure of Leishmania EIF4E5 evidences remarkable differences in the eIF4G interfacing region, when compared with human eIF4E-1 and with its Trypanosoma orthologue. Mapping of its C-terminal end near the cap-binding site also imply relevant differences in cap-binding function and/or regulation.

Identification of the interactomes associated with SCD6 and RBP42 proteins in Leishmania braziliensis

Leishmania are protozoan parasites responsible for leishmaniasis. These parasites present a precise gene regulation that allows them to survive different environmental conditions during their digenetic life cycle. This adaptation depends on the regulation of the expression of a wide variety of genes, which occurs, mainly at the post-transcriptional level. This differential gene expression is achieved by mechanisms based mainly in RNA binding proteins that regulate the translation and/or stability of mRNA targets by interaction with cis elements principally located in the untranslated regions (UTR). In recent studies, our group identified and characterized two proteins, SCD6 and RBP42, as RNA binding proteins in Leishmania braziliensis. To find clues about the cellular processes in which these proteins are involved, this work was aimed to determine the SCD6- and RBP42-interacting proteins (interactome) in L. braziliensis promastigotes. For this purpose, after an in vivo UV cross-linking, cellular extracts were used to immunoprecipitated, by specific antibodies, protein complexes in which SCD6 or RBP42 were present. Protein mass spectrometry analysis of the immunoprecipitated proteins identified 96 proteins presumably associated with SCD6 and 173 proteins associated with RBP42. Notably, a significant proportion of the identified proteins were shared in both interactomes, indicating a possible functional relationship between SCD6 and RBP42. Remarkably, many of the proteins identified in the SCD6 and RBP42 interactomes are related to RNA metabolism and translation processes, and many of them have been described as components of ribonucleoprotein (RNP) granules in Leishmania and related trypanosomatids. Thus, these results support a role of SCD6 and RBP42 in the assembly and/or function of mRNA-protein complexes, participating in the fate (decay/accumulation/translation) of L. braziliensis transcripts. SIGNIFICANCE: Parasites of the Leishmania genus present a particular regulation of gene expression, operating mainly at the post-transcriptional level, surely aimed to modulate quickly both mRNA and protein levels to survive the sudden environmental changes that occur during a parasite's life cycle as it moves from one host to another. This regulation of gene expression processes would be governed by the interaction of mRNA with RNA binding proteins. Nevertheless, the entirety of protein networks involved in these regulatory processes is far from being understood. In this regard, our work is contributing to stablish protein networks in which the L. braziliensis SCD6 and RBP42 proteins are involved; these proteins, in previous works, have been described as RNA binding proteins and found to participate in gene regulation in different cells and organisms. Additionally, our data point out a possible functional relationship between SCD6 and RBP42 proteins as constituents of mRNA granules, like processing bodies or stress granules, which are essential structures in the regulation of gene expression. This knowledge could provide a new approach for the development of therapeutic targets to control Leishmania infections.

Diversity and genome mapping assessment of disordered and functional domains in trypanosomatids

The proteins that have structural disorder exemplify a class of proteins which is part of a new frontier in structural biology that demands a new understanding of the paradigm of structure/function correlations. In order to address the location, relative distances and the functional/structural correlation between disordered and conserved domains, consensus disordered predictions were mapped together with CDD domains in Leishmania braziliensis M2904, Leishmania infantum JPCM5, Trypanosoma cruzi CL-Brener Esmeraldo-like, Trypanosoma cruzi Dm28c, Trypanosoma cruzi Sylvio X10, Blechomonas ayalai B08-376 and Paratrypanosoma confusum CUL13 predicted proteomes. Our results depicts the role of protein disorder in key aspects of parasites biology highlighting: a) statistical significant association between genome structural location of protein disordered consensus stretches and functional domains; b) that disordered protein stretches appear in greater percentage at upstream or downstream position of the predicted domain; c) a possible role of structural disorder in several gene expression, control points that includes but are not limited to: i) protein folding; ii) protein transport and degradation; and iii) protein modification. In addition, for values of protein with disorder content greater than 40%, a small percentage of protein binding sites in IDPs/IDRs, a higher hypothetical protein annotation frequency was observed than expected by chance and trypanosomatid multigene families linked with virulence are rich in protein with disorder content. SIGNIFICANCE: T. cruzi and Leishmania spp are the etiological agents of Chagas disease and leishmaniasis, respectively. Currently, no vaccine or effective drug treatment is available against these neglected diseases and the knowledge about the post-transcriptional and post-translational mechanisms of these organisms, which are key for this scenario, remain scarce. This study depicts the potential impact of the proximity between protein structural disorder and functional domains in the post-transcriptional regulation of pathogenic versus human non-pathogenic trypanosomatids. Our results revealed a significant statistical relationship between the genome structural locations of these two variables and disordered regions appearing more frequently at upstream or downstream positions of the CDD locus domain. This flexibility feature would maintain structural accessibility of functional sites for post-translational modifications, shedding light into this important aspect of parasite biology. This hypothesis is corroborated by the functional enrichment analysis of disordered proteins subset that highlight the involvement of this class of proteins in protein folding, protein transport and degradation and protein modification. Furthermore, our results pointed out: a) the impact of protein disorder in the process of genome annotation (proteins tend to be annotated as hypothetical when the disorder content reaches ~40%); b) that trypanosomatid multigenic families linked with virulence have a key protein disorder content.

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.

Replication origin location might contribute to genetic variability in Trypanosoma cruzi

Background

DNA replication in trypanosomatids operates in a uniquely challenging environment, since most of their genomes are constitutively transcribed. Trypanosoma cruzi, the etiological agent of Chagas disease, presents high variability in both chromosomes size and copy number among strains, though the underlying mechanisms are unknown.

Results

Here we have mapped sites of DNA replication initiation across the T. cruzi genome using Marker Frequency Analysis, which has previously only been deployed in two related trypanosomatids. The putative origins identified in T. cruzi show a notable enrichment of GC content, a preferential position at subtelomeric regions, coinciding with genes transcribed towards the telomeres, and a pronounced enrichment within coding DNA sequences, most notably in genes from the Dispersed Gene Family 1 (DGF-1).

Conclusions

These findings suggest a scenario where collisions between DNA replication and transcription are frequent, leading to increased genetic variability, as seen by the increase SNP levels at chromosome subtelomeres and in DGF-1 genes containing putative origins.

Regulation of Translation in the Protozoan Parasite Leishmania

Leishmaniasis represents a serious health problem worldwide and drug resistance is a growing concern. Leishmania parasites use unusual mechanisms to control their gene expression. In contrast to many other species, they do not have transcriptional regulation. The lack of transcriptional control is mainly compensated by post-transcriptional mechanisms, including tight translational control and regulation of mRNA stability/translatability by RNA-binding proteins. Modulation of translation plays a major role in parasite survival and adaptation to dramatically different environments during change of host; however, our knowledge of fine molecular mechanisms of translation in Leishmania remains limited. Here, we review the current progress in our understanding of how changes in the translational machinery promote parasite differentiation during transmission from a sand fly to a mammalian host, and discuss how translational reprogramming can contribute to the development of drug resistance.

Relevance and proteomics challenge of functional posttranslational modifications in Kinetoplastid parasites

Protozoan parasitic infections are health, social and economic issues impacting both humans and animals, with significant morbidity and mortality worldwide. Protozoan parasites have complicated life cycles with both intracellular and extracellular forms. As a consequence, protozoan adapt to changing environments in part through a dynamic enzyme-catalyzed process leading to reversible posttranslational modifications (PTMs). The characterization by proteomics approaches reveals the critical role of the PTMs of the proteins involved in host-pathogen interaction. The complexity of PTMs characterization is increased by the high diversity, stoichiometry, dynamic and also co-existence of several PTMs in the same moieties which crosstalk between them. Here, we review how to understand the complexity and the essential role of PTMs crosstalk in order to provide a new hallmark for vaccines developments, immunotherapies and personalized medicine. In addition, the importance of these motifs in the biology and biological cycle of kinetoplastid parasites is highlighted with key examples showing the potential to act as targets against protozoan diseases.

Trypanosoma cruzi RNA-binding protein ALBA30 aggregates into cytoplasmic foci under nutritional stress

Trypanosoma cruzi, the causative agent of Chagas disease, has a complex life cycle that requires the adaptation to different environments. In the absence of traditional mechanisms for regulation of gene expression, this parasite relies on posttranscriptional control events, which allow the progression of its life cycle in different hosts and stress conditions. In this context, different stress conditions trigger the aggregation of RNA-binding proteins and their target mRNAs into cytoplasmic foci known as RNA granules, which act as RNA-sorting centers. In this study, we have characterized the T. cruzi RNA-binding protein ALBA30 during nutritional stress conditions. Using a recombinant form of TcALBA30 to facilitate its detection (rTcALBA30), we showed that this protein resides in the cytoplasm in normal growth conditions but is recruited into cytoplasmic foci after starvation. Moreover, evaluation of rTcALBA30 in parasites that reached the stationary phase of growth also showed the recruitment of this protein into cytoplasmic foci. Our results indicate that, similar to TbALBA3, TcALBA30 aggregates into stress granules in parasites submitted to nutritional stress.

The RNA-binding protein TcUBP1 up-regulates an RNA regulon for a cell surface-associated Trypanosoma cruzi glycoprotein and promotes parasite infectivity

The regulation of transcription in trypanosomes is unusual. To modulate protein synthesis during their complex developmental stages, these unicellular microorganisms rely largely on post-transcriptional gene expression pathways. These pathways include a plethora of RNA-binding proteins (RBPs) that modulate all steps of the mRNA life cycle in trypanosomes and help organize transcriptomes into clusters of post-transcriptional regulons. The aim of this work was to characterize an RNA regulon comprising numerous transcripts of trypomastigote-associated cell-surface glycoproteins that are preferentially expressed in the infective stages of the human parasite Trypanosoma cruzi. In vitro and in vivo RNA-binding assays disclosed that these glycoprotein mRNAs are targeted by the small trypanosomatid-exclusive RBP in T. cruzi, U-rich RBP 1 (TcUBP1). Overexpression of a GFP-tagged TcUBP1 in replicative parasites resulted in >10 times up-regulated expression of transcripts encoding surface proteins and in changes in their subcellular localization from the posterior region to the perinuclear region of the cytoplasm, as is typically observed in the infective parasite stages. Moreover, RT-quantitative PCR analysis of actively translated mRNAs by sucrose cushion fractionation revealed an increased abundance of these target transcripts in the polysome fraction of TcUBP1-induced samples. Because these surface proteins are involved in cell adherence or invasion during host infection, we also carried out in vitro infections with TcUBP1-transgenic trypomastigotes and observed that TcUBP1 overexpression significantly increases parasite infectivity. Our findings provide evidence for a role of TcUBP1 in trypomastigote stage-specific gene regulation important for T. cruzi virulence.

Analysis by RNA-seq of transcriptomic changes elicited by heat shock in Leishmania major

Besides their medical relevance, Leishmania is an adequate model for studying post-transcriptional mechanisms of gene expression. In this microorganism, mRNA degradation/stabilization mechanisms together with translational control and post-translational modifications of proteins are the major drivers of gene expression. Leishmania parasites develop as promastigotes in sandflies and as amastigotes in mammalians, and during host transmission, the parasite experiences a sudden temperature increase. Here, changes in the transcriptome of Leishmania major promastigotes after a moderate heat shock were analysed by RNA-seq. Several of the up-regulated transcripts code for heat shock proteins, other for proteins previously reported to be amastigote-specific and many for hypothetical proteins. Many of the transcripts experiencing a decrease in their steady-state levels code for transporters, proteins involved in RNA metabolism or translational factors. In addition, putative long noncoding RNAs were identified among the differentially expressed transcripts. Finally, temperature-dependent changes in the selection of the spliced leader addition sites were inferred from the RNA-seq data, and particular cases were further validated by RT-PCR and Northern blotting. This study provides new insights into the post-transcriptional mechanisms by which Leishmania modulate gene expression.

Comparative sequence and structure analysis of eIF1A and eIF1AD

BACKGROUND

Eukaryotic translation initiation factor 1A (eIF1A) is universally conserved in all organisms. It has multiple functions in translation initiation, including assembly of the ribosomal pre-initiation complexes, mRNA binding, scanning, and ribosomal subunit joining. eIF1A binds directly to the small ribosomal subunit, as well as to several other translation initiation factors. The structure of an eIF1A homolog, the eIF1A domain-containing protein (eIF1AD) was recently determined but its biological functions are unknown. Since eIF1AD has a known structure, as well as a homolog, whose structure and functions have been extensively studied, it is a very attractive target for sequence and structure analysis.

RESULTS

Structure/sequence analysis of eIF1AD found significant conservation in the surfaces corresponding to the ribosome-binding surfaces of its paralog eIF1A, including a nearly invariant surface-exposed tryptophan residue, which plays an important role in the interaction of eIF1A with the ribosome. These results indicate that eIF1AD may bind to the ribosome, similar to its paralog eIF1A, and could have roles in ribosome biogenenesis or regulation of translation. We identified conserved surfaces and sequence motifs in the folded domain as well as the C-terminal tail of eIF1AD, which are likely protein-protein interaction sites. The roles of these regions for eIF1AD function remain to be determined. We have also identified a set of trypanosomatid-specific surface determinants in eIF1A that could be a promising target for development of treatments against these parasites.

CONCLUSIONS

The results described here identify regions in eIF1A and eIF1AD that are likely to play major functional roles and are promising therapeutic targets. Our findings and hypotheses will promote new research and help elucidate the functions of eIF1AD.

Phosphorylation and interactions associated with the control of the Leishmania Poly-A Binding Protein 1 (PABP1) function during translation initiation

The Poly-A Binding Protein (PABP) is a conserved eukaryotic polypeptide involved in many aspects of mRNA metabolism. During translation initiation, PABP interacts with the translation initiation complex eIF4F and enhances the translation of polyadenylated mRNAs. Schematically, most PABPs can be divided into an N-terminal RNA-binding region, a non-conserved linker segment and the C-terminal MLLE domain. In pathogenic Leishmania protozoans, three PABP homologues have been identified, with the first one (PABP1) targeted by phosphorylation and shown to co-immunoprecipitate with an eIF4F-like complex (EIF4E4/EIF4G3) implicated in translation initiation. Here, PABP1 phosphorylation was shown to be linked to logarithmic cell growth, reminiscent of EIF4E4 phosphorylation, and coincides with polysomal association. Phosphorylation targets multiple serine-proline (SP) or threonine-proline (TP) residues within the PABP1 linker region. This is an essential protein, but phosphorylation is not needed for its association with polysomes or cell viability. Mutations which do impair PABP1 polysomal association and are required for viability do not prevent phosphorylation, although further mutations lead to a presumed inactive protein largely lacking phosphorylated isoforms. Co-immunoprecipitation experiments were carried out to investigate PABP1 function further, identifying several novel protein partners and the EIF4E4/EIF4G3 complex, but no other eIF4F-like complex or subunit. A novel, direct interaction between PABP1 and EIF4E4 was also investigated and found to be mediated by the PABP1 MLLE binding to PABP Interacting Motifs (PAM2) within the EIF4E4 N-terminus. The results shown here are consistent with phosphorylation of PABP1 being part of a novel pathway controlling its function and possibly translation in Leishmania.

Leishmania braziliensis SCD6 and RBP42 proteins, two factors with RNA binding capacity

BACKGROUND

The study of RNA binding proteins (RBPs) is of great relevance for understanding processes like post-transcriptional control of gene expression. The post-transcriptional mechanisms are particularly important in Leishmania parasites and related trypanosomatids since transcriptional regulation is almost absent in them. Thus, RBPs should be essential during the development of these parasites and for survival strategies against the adverse conditions that they face during their life-cycle. This work was aimed to do a structural and biochemical characterization of two Leishmania braziliensis proteins, which were previously found in pull-down assays using an HSP70 RNA as bait. At that time, these proteins were annotated as hypothetical proteins (LbrM.25.2210 and LbrM.30.3080) in the GeneDB database.

RESULTS

Structural analysis indicated that these two proteins belong to evolutionarily conserved families; thus, they have been renamed accordingly as LbSCD6 (LbrM.25.2210) and LbRBP42 (LbrM.30.3080). We have demonstrated experimentally that these proteins are RBPs, in agreement with their structural features. Both proteins were able to bind to the complete 3' UTR-II region of HSP70-type II mRNA, and to an A + U rich element (ARE) present in that UTR. Cellular localization assays suggested that both proteins are mainly distributed in the cytoplasm of promastigotes growing at 26 °C, but they accumulate in foci around the nucleus when the parasites are under heat-shock conditions. Also, our study showed that steady-state levels of LbSCD6 and LbRBP42 transcripts decreased significantly during incubation of L. braziliensis promastigotes at heat-shock temperatures. However, in these conditions, the cellular content of both proteins remained unaltered.

CONCLUSIONS

Our data suggest that LbSCD6 and LbRBP42, as occurs for their orthologues in other organisms, are involved in mRNA regulation, and probably they have a relevant role facing the stress conditions that L. braziliensis encounters during insect-to-mammalian transmission.

The Role of Cytoplasmic mRNA Cap-Binding Protein Complexes in Trypanosoma brucei and Other Trypanosomatids

Trypanosomatid protozoa are unusual eukaryotes that are well known for having unusual ways of controlling their gene expression. The lack of a refined mode of transcriptional control in these organisms is compensated by several post-transcriptional control mechanisms, such as control of mRNA turnover and selection of mRNA for translation, that may modulate protein synthesis in response to several environmental conditions found in different hosts. In other eukaryotes, selection of mRNA for translation is mediated by the complex eIF4F, a heterotrimeric protein complex composed by the subunits eIF4E, eIF4G, and eIF4A, where the eIF4E binds to the 5'-cap structure of mature mRNAs. In this review, we present and discuss the characteristics of six trypanosomatid eIF4E homologs and their associated proteins that form multiple eIF4F complexes. The existence of multiple eIF4F complexes in trypanosomatids evokes exquisite mechanisms for differential mRNA recognition for translation.

Comprehensive Identification of mRNA-Binding Proteins of Leishmania donovani by Interactome Capture

Leishmania are unicellular eukaryotes responsible for leishmaniasis in humans. Like other trypanosomatids, leishmania regulate protein coding gene expression almost exclusively at the post-transcriptional level with the help of RNA binding proteins (RBPs). Due to the presence of polycystronic transcription units, leishmania do not regulate RNA polymerase II-dependent transcription initiation. Recent evidence suggests that the main control points in gene expression are mRNA degradation and translation. Protein-RNA interactions are involved in every aspect of RNA biology, such as mRNA splicing, polyadenylation, localization, degradation, and translation. A detailed picture of these interactions would likely prove to be highly informative in understanding leishmania biology and virulence. We developed a strategy involving covalent UV cross-linking of RBPs to mRNA in vivo, followed by interactome capture using oligo(dT) magnetic beads to define comprehensively the mRNA interactome of growing L. donovani amastigotes. The protein mass spectrometry analysis of captured proteins identified 79 mRNA interacting proteins which withstood very stringent washing conditions. Strikingly, we found that 49 of these mRNA interacting proteins had no orthologs or homologs in the human genome. Consequently, these may represent high quality candidates for selective drug targeting leading to novel therapeutics. These results show that this unbiased, systematic strategy has the promise to be applicable to study the mRNA interactome during various biological settings such as metabolic changes, stress (low pH environment, oxidative stress and nutrient deprivation) or drug treatment.

Involvement of an RNA binding protein containing Alba domain in the stage-specific regulation of beta-amastin expression in Trypanosoma cruzi

Amastins are surface glycoproteins, first identified in amastigotes of T. cruzi but later found to be expressed in several Leishmania species, as well as in T. cruzi epimastigotes. Amastins are encoded by a diverse gene family that can be grouped into four subfamilies named α, β, γ, and δ amastins. Differential expression of amastin genes results from regulatory mechanisms involving changes in mRNA stability and/or translational control. Although distinct regulatory elements were identified in the 3' UTR of T. cruzi and Leishmania amastin mRNAs, RNA binding proteins involved with amastin gene regulation have only being characterized in L. infantum where an Alba-domain protein (LiAlba20) able to bind to the 3' UTR of a δ-amastin mRNA was identified. Here we investigated the role of TcAlba30, the LiAlba20 homologue in T. cruzi, in the post transcriptional regulation of amastin genes. TcAlba30 transcripts are present in all stages of the T. cruzi life cycle. RNA immunoprecipitation assays using a transfected cell line expressing a cMyc tagged TcAlba30 revealed that TcAlba30 can interact with β-amastin mRNA. In addition, over-expression of TcAlba30 in epimastigotes resulted in 50% decreased levels of β-amastin mRNAs compared to wild type parasites. Since luciferase assays indicated the presence of regulatory elements in the 3' UTR of β-amastin mRNA and reduced levels of luciferase mRNA were found in parasites over expressing TcAlba30, we conclude that TcAlba30 acts as a T. cruzi RNA binding protein involved in the negative control of β-amastin expression through interactions with its 3'UTR.

T7 polymerase-driven transcription is downregulated in metacyclic promastigotes and amastigotes of Leishmania mexicana

In our previous work we established a T7 polymerase-driven Tetracycline-inducible protein expression system in Leishmania mexicana (Biagi, 1953). We used this system to analyse gene expression profiles during development of L. mexicana in procyclic and metacyclic promastigotes and amastigotes. The transcription of the gene of interest and the T7 polymerase genes was significantly reduced upon cell differentiation. This regulation is not locus-specific. It depends on untranslated regions flanking open reading frames of the genes analysed. In this paper, we report that the previously established conventional inducible protein expression system may not be suitable for studies on differentiation of species of Leishmania Ross, 1903 and protein expression systems might have certain limitations.

RNA-binding proteins related to stress response and differentiation in protozoa

RNA-binding proteins (RBPs) are key regulators of gene expression. There are several distinct families of RBPs and they are involved in the cellular response to environmental changes, cell differentiation and cell death. The RBPs can differentially combine with RNA molecules and form ribonucleoprotein (RNP) complexes, defining the function and fate of RNA molecules in the cell. RBPs display diverse domains that allow them to be categorized into distinct families. They play important roles in the cellular response to physiological stress, in cell differentiation, and, it is believed, in the cellular localization of certain mRNAs. In several protozoa, a physiological stress (nutritional, temperature or pH) triggers differentiation to a distinct developmental stage. Most of the RBPs characterized in protozoa arise from trypanosomatids. In these protozoa gene expression regulation is mostly post-transcriptional, which suggests that some RBPs might display regulatory functions distinct from those described for other eukaryotes. mRNA stability can be altered as a response to stress. Transcripts are sequestered to RNA granules that ultimately modulate their availability to the translation machinery, storage or degradation, depending on the associated proteins. These aggregates of mRNPs containing mRNAs that are not being translated colocalize in cytoplasmic foci, and their numbers and size vary according to cell conditions such as oxidative stress, nutritional status and treatment with drugs that inhibit translation.

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.

Sequence-Based Prediction of RNA-Binding Proteins Using Random Forest with Minimum Redundancy Maximum Relevance Feature Selection

The prediction of RNA-binding proteins is one of the most challenging problems in computation biology. Although some studies have investigated this problem, the accuracy of prediction is still not sufficient. In this study, a highly accurate method was developed to predict RNA-binding proteins from amino acid sequences using random forests with the minimum redundancy maximum relevance (mRMR) method, followed by incremental feature selection (IFS). We incorporated features of conjoint triad features and three novel features: binding propensity (BP), nonbinding propensity (NBP), and evolutionary information combined with physicochemical properties (EIPP). The results showed that these novel features have important roles in improving the performance of the predictor. Using the mRMR-IFS method, our predictor achieved the best performance (86.62% accuracy and 0.737 Matthews correlation coefficient). High prediction accuracy and successful prediction performance suggested that our method can be a useful approach to identify RNA-binding proteins from sequence information.

A DEVH-box RNA Helicase from Leishmania braziliensis is Associated to mRNA Cytoplasmic Granules

RNA helicases are ubiquitous enzymes that participate in almost all aspects of RNA processing, including RNA and RNA-protein complex remodelling. In trypanosomatids, which post-transcriptionally regulate gene expression, the formation of different kinds of ribonucleoprotein granules under stress conditions modulates the parasite's RNA metabolism. This paper describes the isolation of a putative DEVH-box RNA helicase produced by promastigotes of Leishmania braziliensis. Using a Cy3-labelled dT30 oligo, FISH showed the localization of this protein to mRNA granules under starvation stress conditions. The central region of the protein was shown to be responsible for this behaviour.

Association of UBP1 to ribonucleoprotein complexes is regulated by interaction with the trypanosome ortholog of the human multifunctional P32 protein

Regulation of gene expression in trypanosomatid parasitic protozoa is mainly achieved posttranscriptionally. RNA-binding proteins (RBPs) associate to 3' untranslated regions in mRNAs through dedicated domains such as the RNA recognition motif (RRM). Trypanosoma cruzi UBP1 (TcUBP1) is an RRM-type RBP involved in stabilization/degradation of mRNAs. TcUBP1 uses its RRM to associate with cytoplasmic mRNA and to mRNA granules under starvation stress. Here, we show that under starvation stress, TcUBP1 is tightly associated with condensed cytoplasmic mRNA granules. Conversely, under high nutrient/low density-growing conditions, TcUBP1 ribonucleoprotein (RNP) complexes are lax and permeable to mRNA degradation and disassembly. After dissociating from mRNA, TcUBP1 can be phosphorylated only in unstressed parasites. We have identified TcP22, the ortholog of mammalian P32/C1QBP, as an interactor of TcUBP1 RRM. Overexpression of TcP22 decreased the number of TcUBP1 granules in starved parasites in vivo. Endogenous TcUBP1 RNP complexes could be dissociated in vitro by addition of recombinant TcP22, a condition stimulating TcUBP1 phosphorylation. Biochemical and in silico analysis revealed that TcP22 interacts with the RNA-binding surface of TcUBP1 RRM. We propose a model for the decondensation of TcUBP1 RNP complexes in T. cruzi through direct interaction with TcP22 and phosphorylation.

The DRBD13 RNA binding protein is involved in the insect-stage differentiation process of Trypanosoma brucei

DRBD13 RNA-binding protein (RBP) regulates the abundance of AU-rich element (ARE)-containing transcripts in trypanosomes. Here we show that DRBD13 regulates RBP6, the developmentally critical protein in trypanosomatids. We also show DRBD13-specific regulation of transcripts encoding cell surface coat proteins including GPEET2, variable surface glycoprotein (VSG) and invariant surface glycoprotein (ISG). Accordingly, alteration in DRBD13 levels leads to changes in the target mRNA abundance and parasite morphology. The high consistency of the observed phenotype with known cell membrane exchanges that occur during progression of T. brucei through the insect stage of its life cycle suggests that DRBD13 is an important regulator in this largely unknown developmental process.

The response of trypanosomes and other eukaryotes to ER stress and the spliced leader RNA silencing (SLS) pathway in Trypanosoma brucei

The unfolded protein response (UPR) is induced when the quality control machinery of the cell is overloaded with unfolded proteins or when one of the functions of the endoplasmic reticulum (ER) is perturbed. Here, I describe UPR in yeast and mammals, and compare it to what we know about pathogenic fungi and the parasitic protozoans from the order kinetoplastida, focusing on the novel pathway the spliced leader silencing (SLS) in Trypanosoma brucei. Trypanosomes lack conventional transcription regulation, and thus, lack most of the UPR machinery present in other eukaryotes. Trypanosome genes are transcribed in polycistronic units that are processed by trans-splicing and polyadenylation. In trans-splicing, which is essential for processing of each mRNA, an exon known as the spliced leader (SL) is added to all mRNAs from a small RNA, the SL RNA. Under severe ER stress, T. brucei elicits the SLS pathway. In SLS, the transcription of the SL RNA gene is extinguished, and the entire transcription complex dissociates from the SL RNA promoter. Induction of SLS is mediated by an ER-associated kinase (PK3) that migrates to the nucleus, where it phosphorylates the TATA-binding protein (TRF4), leading shut-off of SL RNA transcription. As a result, trans-splicing is inhibited and the parasites activate a programmed cell death (PCD) pathway. Despite the ability to sense the ER stress, the different eukaryotes, especially unicellular parasites and pathogenic fungi, developed a variety of unique and different ways to sense and adjust to this stress in a manner different from their host.

Trypanosoma brucei RRM1 is a nuclear RNA-binding protein and modulator of chromatin structure

TbRRM1 of Trypanosoma brucei is a nucleoprotein that was previously identified in a search for splicing factors in T. brucei. We show that TbRRM1 associates with mRNAs and with the auxiliary splicing factor polypyrimidine tract-binding protein 2, but not with components of the core spliceosome. TbRRM1 also interacts with several retrotransposon hot spot (RHS) proteins and histones. RNA immunoprecipitation of a tagged form of TbRRM1 from procyclic (insect) form trypanosomes identified ca. 1,500 transcripts that were enriched and 3,000 transcripts that were underrepresented compared to cellular mRNA. Enriched transcripts encoded RNA-binding proteins, including TbRRM1 itself, several RHS transcripts, mRNAs with long coding regions, and a high proportion of stage-regulated mRNAs that are more highly expressed in bloodstream forms. Transcripts encoding ribosomal proteins, other factors involved in translation, and procyclic-specific transcripts were underrepresented. Knockdown of TbRRM1 by RNA interference caused widespread changes in mRNA abundance, but these changes did not correlate with the binding of the protein to transcripts, and most splice sites were unchanged, negating a general role for TbRRM1 in splice site selection. When changes in mRNA abundance were mapped across the genome, regions with many downregulated mRNAs were identified. Two regions were analyzed by chromatin immunoprecipitation, both of which exhibited increases in nucleosome occupancy upon TbRRM1 depletion. In addition, subjecting cells to heat shock resulted in translocation of TbRRM1 to the cytoplasm and compaction of chromatin, consistent with a second role for TbRRM1 in modulating chromatin structure.

Trypanosoma brucei, the parasite that causes human sleeping sickness, is transmitted by tsetse flies. The parasite progresses through different life cycle stages in its two hosts, altering its pattern of gene expression in the process. In trypanosomes, protein-coding genes are organized as polycistronic units that are processed into monocistronic mRNAs. Since genes in the same unit can be regulated independently of each other, it is believed that gene regulation is essentially posttranscriptional. In this study, we investigated the role of a nuclear RNA-binding protein, TbRRM1, in the insect stage of the parasite. We found that TbRRM1 binds nuclear mRNAs and also affects chromatin status. Reduction of nuclear TbRRM1 by RNA interference or heat shock resulted in chromatin compaction. We propose that TbRRM1 regulates RNA polymerase II-driven gene expression both cotranscriptionally, by facilitating transcription and efficient splicing, and posttranscriptionally, via its interaction with nuclear mRNAs.

High throughput sequencing analysis of Trypanosoma brucei DRBD3/PTB1-bound mRNAs

Trypanosomes are early-branched eukaryotes that show an unusual dependence on post-transcriptional mechanisms to regulate gene expression. RNA-binding proteins are crucial in controlling mRNA fate in these organisms, but their RNA substrates remain largely unknown. Here we have analyzed on a global scale the mRNAs associated with the Trypanosoma brucei RNA-binding protein DRBD3/PTB1, by capturing ribonucleoprotein complexes using UV cross-linking and subsequent immunoprecipitation. DRBD3/PTB1 associates with many transcripts encoding ribosomal proteins and translation factors. Consequently, silencing of DRBD3/PTB1 expression altered the protein synthesis rate. DRBD3/PTB1 also binds to mRNAs encoding the enzymes required to obtain energy through the oxidation of proline to succinate. We hypothesize that DRBD3/PTB1 is a key player in RNA regulon-based gene control influencing protein synthesis in trypanosomes.

The translation initiation complex eIF3 in trypanosomatids and other pathogenic excavates--identification of conserved and divergent features based on orthologue analysis

Background

The initiation of translation in eukaryotes is supported by the action of several eukaryotic Initiation Factors (eIFs). The largest of these is eIF3, comprising of up to thirteen polypeptides (eIF3a through eIF3m), involved in multiple stages of the initiation process. eIF3 has been better characterized from model organisms, but is poorly known from more diverged groups, including unicellular lineages represented by known human pathogens. These include the trypanosomatids (Trypanosoma and Leishmania) and other protists belonging to the taxonomic supergroup Excavata (Trichomonas and Giardia sp.).

Results

An in depth bioinformatic search was carried out to recover the full content of eIF3 subunits from the available genomes of L. major, T. brucei, T. vaginalis and G. duodenalis. The protein sequences recovered were then submitted to homology analysis and alignments comparing them with orthologues from representative eukaryotes. Eleven putative eIF3 subunits were found from both trypanosomatids whilst only five and four subunits were identified from T. vaginalis and G. duodenalis, respectively. Only three subunits were found in all eukaryotes investigated, eIF3b, eIF3c and eIF3i. The single subunit found to have a related Archaean homologue was eIF3i, the most conserved of the eIF3 subunits. The sequence alignments revealed several strongly conserved residues/region within various eIF3 subunits of possible functional relevance. Subsequent biochemical characterization of the Leishmania eIF3 complex validated the bioinformatic search and yielded a twelfth eIF3 subunit in trypanosomatids, eIF3f (the single unidentified subunit in trypanosomatids was then eIF3m). The biochemical data indicates a lack of association of the eIF3j subunit to the complex whilst highlighting the strong interaction between eIF3 and eIF1.

Conclusions

The presence of most eIF3 subunits in trypanosomatids is consistent with an early evolution of a fully functional complex. Simplified versions in other excavates might indicate a primordial complex or secondary loss of selected subunits, as seen for some fungal lineages. The conservation in eIF3i sequence might indicate critical functions within eIF3 which have been overlooked. The identification of eIF3 subunits from distantly related eukaryotes provides then a basis for the study of conserved/divergent aspects of eIF3 function, leading to a better understanding of eukaryotic translation initiation.

Electronic supplementary material

The online version of this article (doi:10.1186/1471-2164-15-1175) contains supplementary material, which is available to authorized users.

Genome-wide RIP-Chip analysis of translational repressor-bound mRNAs in the Plasmodium gametocyte

BACKGROUND

Following fertilization, the early proteomes of metazoans are defined by the translation of stored but repressed transcripts; further embryonic development relies on de novo transcription of the zygotic genome. During sexual development of Plasmodium berghei, a rodent model for human malaria species including P. falciparum, the stability of repressed mRNAs requires the translational repressors DOZI and CITH. When these repressors are absent, Plasmodium zygote development and transmission to the mosquito vector is halted, as hundreds of transcripts become destabilized. However, which mRNAs are direct targets of these RNA binding proteins, and thus subject to translational repression, is unknown.

RESULTS

We identify the maternal mRNA contribution to post-fertilization development of P. berghei using RNA immunoprecipitation and microarray analysis. We find that 731 mRNAs, approximately 50% of the transcriptome, are associated with DOZI and CITH, allowing zygote development to proceed in the absence of RNA polymerase II transcription. Using GFP-tagging, we validate the repression phenotype of selected genes and identify mRNAs relying on the 5' untranslated region for translational control. Gene deletion reveals a novel protein located in the ookinete crystalloid with an essential function for sporozoite development.

CONCLUSIONS

Our study details for the first time the P. berghei maternal repressome. This mRNA population provides the developing ookinete with coding potential for key molecules required for life-cycle progression, and that are likely to be critical for the transmission of the malaria parasite from the rodent and the human host to the mosquito vector.

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

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

eIF4F-like complexes formed by cap-binding homolog TbEIF4E5 with TbEIF4G1 or TbEIF4G2 are implicated in post-transcriptional regulation in Trypanosoma brucei

Members of the eIF4E mRNA cap-binding family are involved in translation and the modulation of transcript availability in other systems as part of a three-component complex including eIF4G and eIF4A. The kinetoplastids possess four described eIF4E and five eIF4G homologs. We have identified two new eIF4E family proteins in Trypanosoma brucei, and define distinct complexes associated with the fifth member, TbEIF4E5. The cytosolic TbEIF4E5 protein binds cap 0 in vitro. TbEIF4E5 was found in association with two of the five TbEIF4Gs. TbIF4EG1 bound TbEIF4E5, a 47.5-kDa protein with two RNA-binding domains, and either the regulatory protein 14-3-3 II or a 117.5-kDa protein with guanylyltransferase and methyltransferase domains in a potentially dynamic interaction. The TbEIF4G2/TbEIF4E5 complex was associated with a 17.9-kDa hypothetical protein and both 14-3-3 variants I and II. Knockdown of TbEIF4E5 resulted in the loss of productive cell movement, as evidenced by the inability of the cells to remain in suspension in liquid culture and the loss of social motility on semisolid plating medium, as well as a minor reduction of translation. Cells appeared lethargic, as opposed to compromised in flagellar function per se. The minimal use of transcriptional control in kinetoplastids requires these organisms to implement downstream mechanisms to regulate gene expression, and the TbEIF4E5/TbEIF4G1/117.5-kDa complex in particular may be a key player in that process. We suggest that a pathway involved in cell motility is affected, directly or indirectly, by one of the TbEIF4E5 complexes.

Ubiquitination-mediated interaction among domains is responsible for inhibition of RNA endonuclease activity of mRNA cycling sequence binding protein from L. donovani (LdCSBP)

In nearly complete absence of transcriptional regulation, messenger RNA (mRNA) turnover mediated through specific cis-elements plays a predominant role in the control of differential gene expression for the disease causing trypanosomatid parasites. In these organisms, the periodic accumulation of S-phase messages during cell cycle is determined by the presence of one or more copies of a conserved CAUAGAAG octanucleotide motif in the untranslated regions of mRNAs. In our previous studies, a multi-domain cycling sequence binding protein LdCSBP from Leishmania donovani was characterized, which binds specifically to the octamer-containing RNAs via its uniquely arranged CCCH-type Zn fingers and degrades them through its small MutS-related (Smr) endonuclease domain, indicative of its potential role in the turnover of the S-phase mRNAs. Interestingly, the protein is modified by the incorporation of a monoubiquitin residue, and the posttranslational modification inhibits its riboendonuclease activity. However, the mechanism of such inhibition was previously unknown. Here, we establish that the CCCH-type Zn finger domain is the site of ubiquitination in LdCSBP and the interaction of CUE domain of the protein with the ubiquitinated Zn finger domain is responsible for inhibition of its riboendonuclease activity. The findings elucidate an inhibitory mechanism of RNA cleavage through ubiquitination-mediated intramolecular interaction among domains of the enzyme. Furthermore, the riboendonuclease activity is inhibited by anti-leishmanial drug paromomycin suggesting that the regulation of RNA metabolism could be a target of the drug.

A short RNA stem-loop is necessary and sufficient for repression of gene expression during early logarithmic phase in trypanosomes

We have compared the transcriptomes of cultured procyclic Trypanosoma brucei cells in early and late logarithmic phases and found that ∼200 mRNAs were differentially regulated. In late log phase cells, the most upregulated mRNA encoded the nucleobase transporter NT8. The 3' untranslated region (UTR) of NT8 contains a short stem-loop cis-element that is necessary for the regulation of NT8 expression in response to external purine levels. When placed in the 3'-UTR of an unregulated transcript, the cis-element is sufficient to confer regulation in response to purines. To our knowledge, this is the first example of a discrete RNA element that can autonomously regulate gene expression in trypanosomes in response to an external factor and reveals an unprecedented purine-dependent signaling pathway that controls gene expression in eukaryotes.

Depletion of the Trypanosome Pumilio domain protein PUF2 or of some other essential proteins causes transcriptome changes related to coding region length

Pumilio domain RNA-binding proteins are known mainly as posttranscriptional repressors of gene expression that reduce mRNA translation and stability. Trypanosoma brucei has 11 PUF proteins. We show here that PUF2 is in the cytosol, with roughly the same number of molecules per cell as there are mRNAs. Although PUF2 exhibits a low level of in vivo RNA binding, it is not associated with polysomes. PUF2 also decreased reporter mRNA levels in a tethering assay, consistent with a repressive role. Depletion of PUF2 inhibited growth of bloodstream-form trypanosomes, causing selective loss of mRNAs with long open reading frames and increases in mRNAs with shorter open reading frames. Reexamination of published RNASeq data revealed the same trend in cells depleted of some other proteins. We speculate that these length effects could be caused by inhibition of the elongation phase of transcription or by an influence of translation status or polysomal conformation on mRNA decay.

RRP6 from Trypanosoma brucei: crystal structure of the catalytic domain, association with EAP3 and activity towards structured and non-structured RNA substrates

RRP6 is a 3′–5′ exoribonuclease associated to the eukaryotic exosome, a multiprotein complex essential for various RNA processing and degradation pathways. In Trypanosoma brucei, RRP6 associates with the exosome in stoichiometric amounts and was localized in both cytoplasm and nucleus, in contrast to yeast Rrp6 which is exclusively nuclear. Here we report the biochemical and structural characterization of T. brucei RRP6 (TbRRP6) and its interaction with the so-called T. brucei Exosome Associated Protein 3 (TbEAP3), a potential orthologue of the yeast Rrp6 interacting protein, Rrp47. Recombinant TbEAP3 is a thermo stable homodimer in solution, however it forms a heterodimeric complex with TbRRP6 with 1∶1 stoichiometry. The crystallographic structure of the TbRRP6 catalytic core exposes for the first time the native catalytic site of this RNase and also reveals a disulfide bond linking two helices of the HRDC domain. RNA degradation assays show the distributive exoribonuclease activity of TbRRP6 and novel findings regarding the structural range of its RNA substrates. TbRRP6 was able to degrade single and double-stranded RNAs and also RNA substrates containing stem-loops including those with 3′ stem-loop lacking single-stranded extensions. Finally, association with TbEAP3 did not significantly interfere with the TbRRP6 catalytic activity in vitro.

Trypanosome MKT1 and the RNA-binding protein ZC3H11: interactions and potential roles in post-transcriptional regulatory networks

The trypanosome zinc finger protein ZC3H11 binds to AU-rich elements in mRNAs. It is essential for survival of the mammalian-infective bloodstream form, where it stabilizes several mRNAs including some encoding chaperones, and is also required for stabilization of chaperone mRNAs during the heat-shock response in the vector-infective procyclic form. When ZC3H11 was artificially 'tethered' to a reporter mRNA in bloodstream forms it increased reporter expression. We here show that ZC3H11 interacts with trypanosome MKT1 and PBP1, and that domains required for both interactions are necessary for function in the bloodstream-form tethering assay. PBP1 interacts with MKT1, LSM12 and poly(A) binding protein, and localizes to granules during parasite starvation. All of these proteins are essential for bloodstream-form trypanosome survival and increase gene expression in the tethering assay. MKT1 is cytosolic and polysome associated. Using a yeast two-hybrid screen and tandem affinity purification we found that trypanosome MKT1 interacts with multiple RNA-binding proteins and other potential RNA regulators, placing it at the centre of a post-transcriptional regulatory network. A consensus interaction sequence, H(E/D/N/Q)PY, was identified. Recruitment of MKT1-containing regulatory complexes to mRNAs via sequence-specific mRNA-binding proteins could thus control several different post-transcriptional regulons.

The CAF1-NOT complex of trypanosomes

In African trypanosomes, there is no control of transcription initiation by RNA polymerase II at the level of individual protein-coding genes. Transcription is polycistronic, and individual mRNAs are excised by trans-splicing and polyadenylation. As a consequence, trypanosomes are uniquely reliant on post-transcriptional mechanisms for control of gene expression. Rates of mRNA decay vary over up to two orders of magnitude, making these organisms an excellent model system for the study of mRNA degradation processes. The trypanosome CAF1-NOT complex is simpler than that of other organisms, with no CCR4 or NOT4 homolog: it consists of CAF1, NOT1, NOT2, NOT5 NOT9, NOT10, and NOT11. It is important for the initiation of degradation of most, although not all, mRNAs. There is no homolog of NOT4, and Tho and TREX complexes are absent. Functions of the trypanosome NOT complex are therefore likely to be restricted mainly to deadenylation. Mechanisms that cause the NOT complex to deadenylate some mRNAs faster than others must exist, but have not yet been described.

RNA in development: how ribonucleoprotein granules regulate the life cycles of pathogenic protozoa

Ribonucleoprotein (RNP) granules are important posttranscriptional regulators of messenger RNA (mRNA) fate. Several types of RNP granules specifically regulate gene expression during development of multicellular organisms and are commonly referred to as germ granules. The function of germ granules is not entirely understood and probably diverse, but it is generally agreed that one main function is posttranscriptional regulation of gene expression during early development, when transcription is silent. One example is the translational repression of maternally derived mRNAs in oocytes. Here, I hope to show that the need for regulation of gene expression by RNP granules is not restricted to animal development, but plays an equally important role during the development of pathogenic protozoa. Apicomplexa and Trypanosomatidae have complex life cycles with frequent host changes. The need to quickly adapt gene expression to a new environment as well as the ability to suppress translation to survive latencies is critical for successful completion of life cycles. Posttranscriptional gene regulation is not necessarily simpler in protozoa. Apicomplexa surprise with the presence of micro RNA (miRNAs) and upstream open reading frames (µORFs). Trypanosomes have an unusually large repertoire of different RNP granule types. A better understanding of RNP granules in protozoa may help to gain insight into the evolutionary origin of RNP granules: Trypanosomes for example have two types of granules with interesting similarities to animal germ granules.

mRNA localization mechanisms in Trypanosoma cruzi

Asymmetric mRNA localization is a sophisticated tool for regulating and optimizing protein synthesis and maintaining cell polarity. Molecular mechanisms involved in the regulated localization of transcripts are widespread in higher eukaryotes and fungi, but not in protozoa. Trypanosomes are ancient eukaryotes that branched off early in eukaryote evolution. We hypothesized that these organisms would have basic mechanisms of mRNA localization. FISH assays with probes against transcripts coding for proteins with restricted distributions showed a discrete localization of the mRNAs in the cytoplasm. Moreover, cruzipain mRNA was found inside reservosomes suggesting new unexpected functions for this vacuolar organelle. Individual mRNAs were also mobilized to RNA granules in response to nutritional stress. The cytoplasmic distribution of these transcripts changed with cell differentiation, suggesting that localization mechanisms might be involved in the regulation of stage-specific protein expression. Transfection assays with reporter genes showed that, as in higher eukaryotes, 3'UTRs were responsible for guiding mRNAs to their final location. Our results strongly suggest that Trypanosoma cruzi have a core, basic mechanism of mRNA localization. This kind of controlled mRNA transport is ancient, dating back to early eukaryote evolution.

The trypanosome Pumilio domain protein PUF5

PUF proteins are a conserved family of RNA binding proteins found in all eukaryotes examined so far. This study focussed on PUF5, one of 11 PUF family members encoded in the Trypanosoma brucei genome. Native PUF5 is present at less than 50000 molecules per cell in both bloodstream and procyclic form trypanosomes. C-terminally myc-tagged PUF5 was mainly found in the cytoplasm and could be cross-linked to RNA. PUF5 knockdown by RNA interference had no effect on the growth of bloodstream forms. Procyclic forms lacking PUF5 grew normally, but expression of PUF5 bearing a 21 kDa tandem affinity purification tag inhibited growth. Knockdown of PUF5 did not have any effect on the ability of trypanosomes to differentiate from the mammalian to the insect form of the parasite.

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

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

The roles of 3'-exoribonucleases and the exosome in trypanosome mRNA degradation

The degradation of eukaryotic mRNAs can be initiated by deadenylation, decapping, or endonuclease cleavage. This is followed by 5'-3' degradation by homologs of Xrn1, and/or 3'-5' degradation by the exosome. We previously reported that, in African trypanosome Trypanosoma brucei, most mRNAs are deadenylated prior to degradation, and that depletion of the major 5'-3' exoribonuclease XRNA preferentially stabilizes unstable mRNAs. We now show that depletion of either CAF1 or CNOT10, two components of the principal deadenylation complex, strongly inhibits degradation of most mRNAs. RNAi targeting another deadenylase, PAN2, or RRP45, a core component of the exosome, preferentially stabilized mRNAs with intermediate half-lives. RRP45 depletion resulted in a 5' bias of mRNA sequences, suggesting action by a distributive 3'-5' exoribonuclease. Results suggested that the exosome is involved in the processing of trypanosome snoRNAs. There was no correlation between effects on half-lives and on mRNA abundance.