Developmental regulation of gene expression in trypanosomatid parasitic protozoa

Evolutionary divergent clusters of transcribed extinct truncated retroposons drive low mRNA expression and developmental regulation in the protozoan Leishmania

BACKGROUND

The Leishmania genome harbors formerly active short interspersed degenerated retroposons (SIDERs) representing the largest family of repetitive elements among trypanosomatids. Their substantial expansion in Leishmania is a strong predictor of important biological functions. In this study, we combined multilevel bioinformatic predictions with high-throughput genomic and transcriptomic analyses to gain novel insights into the diversified roles retroposons of the SIDER2 subfamily play in Leishmania genome evolution and expression.

RESULTS

We show that SIDER2 retroposons form various evolutionary divergent clusters, each harboring homologous SIDER2 sequences usually located nearby in the linear sequence of chromosomes. This intriguing genomic organization underscores the importance of SIDER2 proximity in shaping chromosome dynamics and co-regulation. Accordingly, we show that transcripts belonging to the same SIDER2 cluster can display similar levels of expression. SIDER2 retroposons are mostly transcribed as part of 3'UTRs and account for 13% of the Leishmania transcriptome. Genome-wide expression profiling studies underscore SIDER2 association generally with low mRNA expression. The remarkable link of SIDER2 retroposons with downregulation of gene expression supports their co-option as major regulators of mRNA abundance. SIDER2 sequences also add to the diversification of the Leishmania gene expression repertoire since ~ 35% of SIDER2-containing transcripts can be differentially regulated throughout the parasite development, with a few encoding key virulence factors. In addition, we provide evidence for a functional bias of SIDER2-containing transcripts with protein kinase and transmembrane transporter activities being most represented.

CONCLUSIONS

Altogether, these findings provide important conceptual advances into evolutionary innovations of transcribed extinct retroposons acting as major RNA cis-regulators.

Molecular Mechanisms of Drug Resistance in Leishmania spp

The protozoan parasite Leishmania causes leishmaniasis, a neglected tropical disease, that disproportionately affects underdeveloped countries. This disease has major health, economic, and social implications, particularly because of the limited treatment options, high cost, the severe side effects associated with available therapeutics, and the high rate of treatment failure caused by the parasites' growing resistance to current medications. In this review, we describe first the common strategies used by pathogens to develop drug resistance and then focus on the arsenal of available drugs to treat leishmaniasis, their modes of action, and the molecular mechanisms contributing to drug resistance in Leishmania spp., including the role of genomic, transcriptional, and translational control. We focus more specifically on our recent discovery of translational reprogramming as a major driver of drug resistance leading to coordinated changes in the translation of transcripts and orchestrating changes in metabolome and lipidome to support drug resistance. A thorough understanding of these mechanisms is essential to identify the key elements needed to combat resistance and improve leishmaniasis treatment methods.

Interaction of novel proteins, centrin4 and protein of centriole in Leishmania parasite and their effects on the parasite growth

Centrins are cytoskeletal proteins associated with the centrosomes or basal bodies in the eukaryotes. We previously reported the involvement of Centrin 1-3 proteins in cell division in the protozoan parasites Leishmania donovani and Trypanosoma brucei. Centrin4 and 5, unique to such parasites, had never been characterized in Leishmania parasite. In the current study, we addressed the function of centrin4 (LdCen4) in Leishmania. By dominant-negative study, the episomal expression of C-terminal truncated LdCen4 in the parasite reduced the parasite growth. LdCen4 double allele gene deletion by either homologous recombination or CRISPR-Cas9 was not successful in L. donovani. However, CRISPR-Cas9-based deletion of the homologous gene was possible in L. mexicana, which attenuated the parasite growth in vitro, but not ex vivo in the macrophages. LdCen4 also interacts with endogenous and overexpressed LdPOC protein, a homolog of centrin reacting human POC (protein of centriole) in a calcium sensitive manner. LdCen4 and LdPOC binding has also been confirmed through in silico analysis by protein structural docking and validated by co-immunoprecipitation. By immunofluorescence studies, we found that both the proteins share a common localization at the basal bodies. Thus, for the first time, this article describes novel centrin4 and its binding protein in the protozoan parasites.

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

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

Finding Correlations Between mRNA and Protein Levels in Leishmania Development: Is There a Discrepancy?

Multiple genes and proteins have been identified as differentially expressed in the stages of the Leishmania life cycle. The differentiation processes are implicated in specific transcriptional and proteomic adjustments driven by gene expression regulation mechanisms. Leishmania parasites lack gene-specific transcriptional control, and gene expression regulation mostly depends on posttranscriptional mechanisms. Due to the lack of transcriptional regulation, criticism regarding the relevance of transcript quantification as a possible and efficient prediction of protein levels is recurrent in studies that use transcriptomic information. The advent of high-throughput technologies has improved the analysis of genomes, transcriptomes and proteomes for different organisms under several conditions. Nevertheless, defining the correlation between transcriptional and proteomic profiles requires arduous and expensive work and remains a challenge in Leishmania. In this review, we analyze transcriptomic and proteomic data for several Leishmania species in two different stages of the parasite life cycle: metacyclogenesis and amastigogenesis (amastigote differentiation). We found a correlation between mRNA and protein levels of 60.9% and 69.8% for metacyclogenesis and amastigogenesis, respectively; showing that majority mRNA and protein levels increase or decrease concomitantly. Among the analyzed genes that did not present correlation indicate that transcriptomic data should be carefully interpreted as protein expression. We also discuss possible explanations and mechanisms involved for this lack of correlation.

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.

Spatially and Temporally Distributed Complexity-A Refreshed Framework for the Study of GRN Evolution

Irrespective of the heuristic value of interpretations of developmental processes in terms of gene regulatory networks (GRNs), larger-angle views often suffer from: (i) an inadequate understanding of the relationship between genotype and phenotype; (ii) a predominantly zoocentric vision; and (iii) overconfidence in a putatively hierarchical organization of animal body plans. Here, we constructively criticize these assumptions. First, developmental biology is pervaded by adultocentrism, but development is not necessarily egg to adult. Second, during development, many unicells undergo transcriptomic profile transitions that are comparable to those recorded in pluricellular organisms; thus, their study should not be neglected from the GRN perspective. Third, the putatively hierarchical nature of the animal body is mirrored in the GRN logic, but in relating genotype to phenotype, independent assessments of the dynamics of the regulatory machinery and the animal’s architecture are required, better served by a combinatorial than by a hierarchical approach. The trade-offs between spatial and temporal aspects of regulation, as well as their evolutionary consequences, are also discussed. Multicellularity may derive from a unicell’s sequential phenotypes turned into different but coexisting, spatially arranged cell types. In turn, polyphenism may have been a crucial mechanism involved in the origin of complex life cycles.

Experimental evolution links post-transcriptional regulation to Leishmania fitness gain

The protozoan parasite Leishmania donovani causes fatal human visceral leishmaniasis in absence of treatment. Genome instability has been recognized as a driver in Leishmania fitness gain in response to environmental change or chemotherapy. How genome instability generates beneficial phenotypes despite potential deleterious gene dosage effects is unknown. Here we address this important open question applying experimental evolution and integrative systems approaches on parasites adapting to in vitro culture. Phenotypic analyses of parasites from early and late stages of culture adaptation revealed an important fitness tradeoff, with selection for accelerated growth in promastigote culture (fitness gain) impairing infectivity (fitness costs). Comparative genomics, transcriptomics and proteomics analyses revealed a complex regulatory network associated with parasite fitness gain, with genome instability causing highly reproducible, gene dosage-independent and -dependent changes. Reduction of flagellar transcripts and increase in coding and non-coding RNAs implicated in ribosomal biogenesis and protein translation were not correlated to dosage changes of the corresponding genes, revealing a gene dosage-independent, post-transcriptional mechanism of regulation. In contrast, abundance of gene products implicated in post-transcriptional regulation itself correlated to corresponding gene dosage changes. Thus, RNA abundance during parasite adaptation is controled by direct and indirect gene dosage changes. We correlated differential expression of small nucleolar RNAs (snoRNAs) with changes in rRNA modification, providing first evidence that Leishmania fitness gain in culture may be controlled by post-transcriptional and epitranscriptomic regulation. Our findings propose a novel model for Leishmania fitness gain in culture, where differential regulation of mRNA stability and the generation of modified ribosomes may potentially filter deleterious from beneficial gene dosage effects and provide proteomic robustness to genetically heterogenous, adapting parasite populations. This model challenges the current, genome-centric approach to Leishmania epidemiology and identifies the Leishmania transcriptome and non-coding small RNome as potential novel sources for the discovery of biomarkers that may be associated with parasite phenotypic adaptation in clinical settings.

Genome instability plays a central yet poorly understood role in human disease. Gene amplifications and deletions drive cancer development, microbial infection and therapeutic failure. The molecular mechanisms that harness the deleterious effects of genome instability to generate beneficial phenotypes in pathogenic systems are unknown. Here we study this important open question in the protozoan parasite Leishmania that causes devastating human diseases termed leishmaniases. Leishmania parasites lack transcriptional control and instead exploit genome instability to adapt to their host environment. Analyzing in vitro adaptation of hamster-derived parasites via gene copy number (genomic level) and gene expression changes (transcriptomic and proteomic levels), we show that these parasites likely exploit small nucleolar RNAs (snoRNAs) to mitigate toxic effects of genome instability by post-transcriptional regulation and the establishment of modified ribosomes. Our findings propose non-coding RNAs as potential novel biomarkers with diagnostic and prognostic value that may be linked to changes in parasite tissue tropism or drug susceptibility. This novel insight into Leishmania adaptation will be likely applicable to other fast evolving eukaryotic systems with unstable genomes, such as fungi or cancer cells.

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

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

Well-Tolerated Amphotericin B Derivatives That Effectively Treat Visceral Leishmaniasis

Chemotherapy against the neglected tropical disease visceral leishmaniasis (VL) is suboptimal with only four licensed drugs. Amphotericin B (AmB), despite its toxicity, remained a second line drug for a long time. However, the demonstration that liposomal AmB is highly effective against VL propelled it, despite its cost, to a first line drug in many countries. While several ongoing efforts are aiming at finding cheaper and stable AmB-formulations, an alternative strategy is the development of less-toxic AmB derivatives. We show here that two less-toxic AmB derivatives with the carboxylate at position 16 of AmB derivatized to a methyl urea (AmB-MU) or amino urea (AmB-AU) are active in vitro against Leishmania donovani, both as free-living parasites as well as their intracellular form. Both less-toxic derivatives, similarly to AmB, target the ergosterol pathway of L. donovani. While the AmB-AU derivative showed female-specific liver toxicity in vivo, the AmB-MU derivative was well-tolerated and more effective than AmB against experimental VL. These studies are an important step for improving AmB-based therapy against a prevalent parasitic disease.

Computational study on the allosteric mechanism of Leishmania major IF4E-1 by 4E-interacting protein-1: Unravelling the determinants of m7GTP cap recognition

•

Atomistic details on perturbations induced by Lm4E-IP1 binding are described.

•

The modulation of LmIF4E-1 affinity for the cap is confirmed by energetic analyses.

•

Signaling paths between the allosteric and othosteric sites of LmIF4E-1 are predicted.

•

Lm4E-IP1 binding increases the side-chain entropy of W83 and R172 of LmIF4E-1.

•

A mechanism of dynamic allostery is proposed for the regulation mediated by Lm4E-IP1.

During their life cycle, Leishmania parasites display a fine-tuned regulation of the mRNA translation through the differential expression of isoforms of eukaryotic translation initiation factor 4E (LeishIF4Es). The interaction between allosteric modulators such as 4E-interacting proteins (4E-IPs) and LeishIF4E affects the affinity of this initiation factor for the mRNA cap. Here, several computational approaches were employed to elucidate the molecular bases of the previously-reported allosteric modulation in L. major exerted by 4E-IP1 (Lm4E-IP1) on eukaryotic translation initiation factor 4E 1 (LmIF4E-1). Molecular dynamics (MD) simulations and accurate binding free energy calculations (ΔG_bind) were combined with network-based modeling of residue-residue correlations. We also describe the differences in internal motions of LmIF4E-1 apo form, cap-bound, and Lm4E-IP1-bound systems. Through community network calculations, the differences in the allosteric pathways of allosterically-inhibited and active forms of LmIF4E-1 were revealed. The ΔG_bind values show significant differences between the active and inhibited systems, which are in agreement with the available experimental data. Our study thoroughly describes the dynamical perturbations of LmIF4E-1 cap-binding site triggered by Lm4E-IP1. These findings are not only essential for the understanding of a critical process of trypanosomatids’ gene expression but also for gaining insight into the allostery of eukaryotic IF4Es, which could be useful for structure-based design of drugs against this protein family.

Distinct features of the Leishmania cap-binding protein LeishIF4E2 revealed by CRISPR-Cas9 mediated hemizygous deletion

Leishmania parasites cycle between sand-fly vectors and mammalian hosts adapting to alternating environments by stage-differentiation accompanied by changes in the proteome profiles. Translation regulation plays a central role in driving the differential program of gene expression since control of gene regulation in Leishmania is mostly post-transcriptional. The Leishmania genome encodes six eIF4E paralogs, some of which bind a dedicated eIF4G candidate, and each eIF4E is assumed to have specific functions with perhaps some overlaps. However, LeishIF4E2 does not bind any known eIF4G ortholog and was previously shown to comigrate with the polysomal fractions of sucrose gradients in contrast to the other initiation factors that usually comigrate with pre-initiation and initiation complexes. Here we deleted one of the two LeishIF4E2 gene copies using the CRISPR-Cas9 methodology. The deletion caused severe alterations in the morphology of the mutant cells that became round, small, and equipped with a very short flagellum that did not protrude from its pocket. Reduced expression of LeishIF4E2 had no global effect on translation and growth, unlike other LeishIF4Es; however, there was a change in the proteome profile of the LeishIF4E2(+/-) cells. Upregulated proteins were related mainly to general metabolic processes including enzymes involved in fatty acid metabolism, DNA repair and replication, signaling, and cellular motor activity. The downregulated proteins included flagellar rod and cytoskeletal proteins, as well as surface antigens involved in virulence. Moreover, the LeishIF4E2(+/-) cells were impaired in their ability to infect cultured macrophages. Overall, LeishIF4E2 does not behave like a general translation factor and its function remains elusive. Our results also suggest that the individual LeishIF4Es perform unique functions.

Leishmania parasites cause a broad spectrum of diseases with different pathological symptoms. During their life cycle the parasites shuffle between sand-fly vectors and mammalian hosts adapting to the changing environments via a stage specific program of gene expression that promotes their survival. Translation initiation plays a key role in control of gene expression and in Leishmania this is exemplified by the presence of multiple cap-binding complexes that interact with mRNAs. The parasites encode multiple paralogs of the cap-binding translation initiation factor eIF4E and of its corresponding binding partner eIF4G forming complexes with different potential functions. The role of LeishIF4E2 remains elusive: it does not bind any of the LeishIF4G candidate subunits and associates with polysomes, a feature less common for canonical translation factors. Here we generated a hemizygous Leishmania mutant of the least studied cap-binding paralog, LeishIF4E2, by eliminating one of the two alleles using the CRISPR-Cas9 methodology. The mutant showed morphological defects with short and rounded cells, and a significant reduction in their flagellar length. Moreover, the LeishIF4E2(+/-) cells were impaired in their ability to infect cultured macrophages. The mutants showed differences in their proteome: upregulated proteins were related mainly to general metabolic processes including enzymes involved in fatty acid metabolism, DNA repair and replication, signaling, and cellular motor activity. Downregulated proteins included flagellar rod and cytoskeletal proteins, as well as surface antigens involved in virulence. Overall, LeishIF4E2 does not behave like a general translation factor and its function remains elusive. It could affect translation of a particular set of transcripts, causing direct or downstream effects that do not affect global translation. Our results suggest that individual LeishIF4Es perform specific functions.

Gene co-expression network analysis of Trypanosoma brucei in tsetse fly vector

Background

Trypanosoma brucei species are motile protozoan parasites that are cyclically transmitted by tsetse fly (genus Glossina) causing human sleeping sickness and nagana in livestock in sub-Saharan Africa. African trypanosomes display digenetic life cycle stages in the tsetse fly vector and in their mammalian host. Experimental work on insect-stage trypanosomes is challenging because of the difficulty in setting up successful in vitro cultures. Therefore, there is limited knowledge on the trypanosome biology during its development in the tsetse fly. Consequently, this limits the development of new strategies for blocking parasite transmission in the tsetse fly.

Methods

In this study, RNA-Seq data of insect-stage trypanosomes were used to construct a T. brucei gene co-expression network using the weighted gene co-expression analysis (WGCNA) method. The study identified significant enriched modules for genes that play key roles during the parasite’s development in tsetse fly. Furthermore, potential 3′ untranslated region (UTR) regulatory elements for genes that clustered in the same module were identified using the Finding Informative Regulatory Elements (FIRE) tool.

Results

A fraction of gene modules (12 out of 27 modules) in the constructed network were found to be enriched in functional roles associated with the cell division, protein biosynthesis, mitochondrion, and cell surface. Additionally, 12 hub genes encoding proteins such as RNA-binding protein 6 (RBP6), arginine kinase 1 (AK1), brucei alanine-rich protein (BARP), among others, were identified for the 12 significantly enriched gene modules. In addition, the potential regulatory elements located in the 3′ untranslated regions of genes within the same module were predicted.

Conclusions

The constructed gene co-expression network provides a useful resource for network-based data mining to identify candidate genes for functional studies. This will enhance understanding of the molecular mechanisms that underlie important biological processes during parasite’s development in tsetse fly. Ultimately, these findings will be key in the identification of potential molecular targets for disease control.

Purification and characterization of a novel and conserved TPR-domain protein that binds both Hsp90 and Hsp70 and is expressed in all developmental stages of Leishmania major

Heat shock proteins (Hsps) are involved in several important aspects of the cell proteostasis. Hsp90 interacts with at least a tenth of the cell proteome helping a large number of proteins to fold correctly. Hsp90 function is modulated by several co-chaperones having TPR (tetratricopeptide repeat) domains that allow for interaction with the C-terminal MEEVD motif of the chaperone. Another important chaperone, Hsp70, has a C-terminal EEVD motif that binds to TPR. Leishmania is a protozoan that causes leishmaniasis, a neglected disease in humans and other animals. There is still no effective treatment for leishmaniasis, however the study of structure and function of the proteins of the parasite may generate potential targets for future therapeutic intervention studies. In this work, the genome of Leishmania major was searched for a novel TPR-domain gene, which is conserved in Leishmania. The recombinant protein, LmTPR, was produced in pure and folded state and was characterized by biophysical tools as a monomer with an elongated conformation. Studies in Leishmania major were also preformed to complement these in vitro experiments. Splice Leader RNA-seq analysis and Western blot indicated that the protein was expressed in all developmental stages of the parasite. Binding assays confirmed that both Hsp90 and Hsp70 bind specifically to LmTPR. Finally, sequence and structural predictions indicated a C-terminal region as a RPAP3 domain. Altogether, this study identified a novel TPR-domain co-chaperone of Hsp90 that is conserved and expressed in all developmental stages of Leishmania major.

The AAA + ATPase valosin-containing protein (VCP)/p97/Cdc48 interaction network in Leishmania

Valosin-containing protein (VCP)/p97/Cdc48 is an AAA + ATPase associated with many ubiquitin-dependent cellular pathways that are central to protein quality control. VCP binds various cofactors, which determine pathway selectivity and substrate processing. Here, we used co-immunoprecipitation and mass spectrometry studies coupled to in silico analyses to identify the Leishmania infantum VCP (LiVCP) interactome and to predict molecular interactions between LiVCP and its major cofactors. Our data support a largely conserved VCP protein network in Leishmania including known but also novel interaction partners. Network proteomics analysis confirmed LiVCP-cofactor interactions and provided novel insights into cofactor-specific partners and the diversity of LiVCP complexes, including the well-characterized VCP-UFD1-NPL4 complex. Gene Ontology analysis coupled with digitonin fractionation and immunofluorescence studies support cofactor subcellular compartmentalization with either cytoplasmic or organellar or vacuolar localization. Furthermore, in silico models based on 3D homology modeling and protein-protein docking indicated that the conserved binding modules of LiVCP cofactors, except for NPL4, interact with specific binding sites in the hexameric LiVCP protein, similarly to their eukaryotic orthologs. Altogether, these results allowed us to build the first VCP protein interaction network in parasitic protozoa through the identification of known and novel interacting partners potentially associated with distinct VCP complexes.

A newly identified Leishmania IF4E-interacting protein, Leish4E-IP2, modulates the activity of cap-binding protein paralogs

Translation of most cellular mRNAs in eukaryotes proceeds through a cap-dependent pathway, whereby the cap-binding complex, eIF4F, anchors the preinitiation complex at the 5′ end of mRNAs and regulates translation initiation. The requirement of Leishmania to survive in changing environments can explain why they encode multiple eIF4E (LeishIF4Es) and eIF4G (LeishIF4Gs) paralogs, as each could be assigned a discrete role during their life cycle. Here we show that the expression and activity of different LeishIF4Es change during the growth of cultured promastigotes, urging a search for regulatory proteins. We describe a novel LeishIF4E-interacting protein, Leish4E-IP2, which contains a conserved Y(X)₄LΦ IF4E-binding-motif. Despite its capacity to bind several LeishIF4Es, Leish4E-IP2 was not detected in m⁷GTP-eluted cap-binding complexes, suggesting that it could inhibit the cap-binding activity of LeishIF4Es. Using a functional assay, we show that a recombinant form of Leish4E-IP2 inhibits the cap-binding activity of LeishIF4E-1 and LeishIF4E-3. Furthermore, we show that transgenic parasites expressing a tagged version of Leish4E-IP2 also display reduced cap-binding activities of tested LeishIF4Es, and decreased global translation. Given its ability to bind more than a single LeishIF4E, we suggest that Leish4E-IP2 could serve as a broad-range repressor of Leishmania protein synthesis.

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.

EIF2α phosphorylation is regulated in intracellular amastigotes for the generation of infective Trypanosoma cruzi trypomastigote forms

Trypanosomatids regulate gene expression mainly at the post-transcriptional level through processing, exporting and stabilising mRNA and control of translation. In most eukaryotes, protein synthesis is regulated by phosphorylation of eukaryotic initiation factor 2 (eIF2) at serine 51. Phosphorylation halts overall translation by decreasing availability of initiator tRNA^met to form translating ribosomes. In trypanosomatids, the N-terminus of eIF2α is extended with threonine 169 the homologous phosphorylated residue. Here, we evaluated whether eIF2α phosphorylation varies during the Trypanosoma cruzi life cycle, the etiological agent of Chagas' disease. Total levels of eIF2α are diminished in infective and non-replicative trypomastigotes compared with proliferative forms from the intestine of the insect vector or amastigotes from mammalian cells, consistent with decreased protein synthesis reported in infective forms. eIF2α phosphorylation increases in proliferative intracellular forms prior to differentiation into trypomastigotes. Parasites overexpressing eIF2α^T169A or with an endogenous CRISPR/Cas9-generated eIF2α^T169A mutation were created and analysis revealed alterations to the proteome, largely unrelated to the presence of μORF in epimastigotes. eIF2α^T169A mutant parasites produced fewer trypomastigotes with lower infectivity than wild type, with increased levels of sialylated mucins and oligomannose glycoproteins, and decreased galactofuranose epitopes and the surface protease GP63 on the cell surface. We conclude that eIF2α expression and phosphorylation levels affect proteins relevant for intracellular progression of T. cruzi.

Experimental Strategies to Explore Drug Action and Resistance in Kinetoplastid Parasites

Kinetoplastids are the causative agents of leishmaniasis, human African trypanosomiasis, and American trypanosomiasis. They are responsible for high mortality and morbidity in (sub)tropical regions. Adequate treatment options are limited and have several drawbacks, such as toxicity, need for parenteral administration, and occurrence of treatment failure and drug resistance. Therefore, there is an urgency for the development of new drugs. Phenotypic screening already allowed the identification of promising new chemical entities with anti-kinetoplastid activity potential, but knowledge on their mode-of-action (MoA) is lacking due to the generally applied whole-cell based approach. However, identification of the drug target is essential to steer further drug discovery and development. Multiple complementary techniques have indeed been used for MoA elucidation. In this review, the different 'omics' approaches employed to define the MoA or mode-of-resistance of current reference drugs and some new anti-kinetoplastid compounds are discussed.

Transcriptomic Analysis of the Early Strobilar Development of Echinococcus granulosus

Echinococcus granulosus has a complex life cycle involving two mammalian hosts. The transition from one host to another is accompanied by changes in gene expression, and the transcriptional events that underlie this transition have not yet been fully characterized. In this study, RNA-seq was used to compare the transcription profiles of samples from E. granulosus protoscoleces induced in vitro to strobilar development at three time points. We identified 818 differentially expressed genes, which were divided into eight expression clusters formed over the entire 24 h period. An enrichment of gene transcripts with molecular functions of signal transduction, enzymes, and protein modifications was observed upon induction and developmental progression. This transcriptomic study provides insights for understanding the complex life cycle of E. granulosus and contributes for searching for the key genes correlating with the strobilar development, which can be used to identify potential candidates for the development of anthelmintic drugs.

Defeating the trypanosomatid trio: proteomics of the protozoan parasites causing neglected tropical diseases

Mass spectrometry-based proteomics enables accurate measurement of the modulations of proteins on a large scale upon perturbation and facilitates the understanding of the functional roles of proteins in biological systems. It is a particularly relevant methodology for studying Leishmania spp., Trypanosoma cruzi and Trypanosoma brucei, as the gene expression in these parasites is primarily regulated by posttranscriptional mechanisms. Large-scale proteomics studies have revealed a plethora of information regarding modulated proteins and their molecular interactions during various life processes of the protozoans, including stress adaptation, life cycle changes and interactions with the host. Important molecular processes within the parasite that regulate the activity and subcellular localisation of its proteins, including several co- and post-translational modifications, are also accurately captured by modern proteomics mass spectrometry techniques. Finally, in combination with synthetic chemistry, proteomic techniques facilitate unbiased profiling of targets and off-targets of pharmacologically active compounds in the parasites. This provides important data sets for their mechanism of action studies, thereby aiding drug development programmes.

An immunoinformatic approach driven by experimental proteomics: in silico design of a subunit candidate vaccine targeting secretory proteins of Leishmania donovani amastigotes

Background

Visceral leishmaniasis (VL) caused by dimorphic Leishmania species is a parasitic disease with high socioeconomic burden in endemic areas worldwide. Sustaining control of VL in terms of proper and prevailing immunity development is a global necessity amid unavailability of a prophylactic vaccine. Screening of experimental proteome of the human disease propagating form of Leishmania donovani (amastigote) can be more pragmatic for in silico mining of novel vaccine candidates.

Methods

By using an immunoinformatic approach, CD4+ and CD8+ T cell-specific epitopes from experimentally reported L. donovani proteins having secretory potential and increased abundance in amastigotes were screened. A chimera linked with a Toll-like receptor 4 (TLR4) peptide adjuvant was constructed and evaluated for physicochemical characteristics, binding interaction with TLR4 in simulated physiological condition and the trend of immune response following hypothetical immunization.

Results

Selected epitopes from physiologically important L. donovani proteins were found mostly conserved in L. infantum, covering theoretically more than 98% of the global population. The multi-epitope chimeric vaccine was predicted as stable, antigenic and non-allergenic. Structural analysis of vaccine-TLR4 receptor docked complex and its molecular dynamics simulation suggest sufficiently stable binding interface along with prospect of non-canonical receptor activation. Simulation dynamics of immune response following hypothetical immunization indicate active and memory B as well as CD4+ T cell generation potential, and likely chance of a more Th1 polarized response.

Conclusions

The methodological approach and results from this study could facilitate more informed screening and selection of candidate antigenic proteins for entry into vaccine production pipeline in future to control human VL.

A BONCAT-iTRAQ method enables temporally resolved quantitative profiling of newly synthesised proteins in Leishmania mexicana parasites during starvation

Adaptation to starvation is integral to the Leishmania life cycle. The parasite can survive prolonged periods of nutrient deprivation both in vitro and in vivo. The identification of parasite proteins synthesised during starvation is key to unravelling the underlying molecular mechanisms facilitating adaptation to these conditions. Additionally, as stress adaptation mechanisms in Leishmania are linked to virulence as well as infectivity, profiling of the complete repertoire of Newly Synthesised Proteins (NSPs) under starvation is important for drug target discovery. However, differential identification and quantitation of low abundance, starvation-specific NSPs from the larger background of the pre-existing parasite proteome has proven difficult, as this demands a highly selective and sensitive methodology. Herein we introduce an integrated chemical proteomics method in L. mexicana promastigotes that involves a powerful combination of the BONCAT technique and iTRAQ quantitative proteomics Mass Spectrometry (MS), which enabled temporally resolved quantitative profiling of de novo protein synthesis in the starving parasite. Uniquely, this approach integrates the high specificity of the BONCAT technique for the NSPs, with the high sensitivity and multiplexed quantitation capability of the iTRAQ proteomics MS. Proof-of-concept experiments identified over 250 starvation-responsive NSPs in the parasite. Our results show a starvation-specific increased relative abundance of several translation regulating and stress-responsive proteins in the parasite. GO analysis of the identified NSPs for Biological Process revealed translation (enrichment P value 2.47e^-35) and peptide biosynthetic process (enrichment P value 4.84e^-35) as extremely significantly enriched terms indicating the high specificity of the NSP towards regulation of protein synthesis. We believe that this approach will find widespread use in the study of the developmental stages of Leishmania species and in the broader field of protozoan biology.

Periodic nutrient scarcity plays crucial roles in the life cycle of the protozoan parasite Leishmania spp. Although adaptation to nutrient stress has a pivotal role in Leishmania biology, the underlying mechanisms remain poorly understood. In a period of nutrient starvation, the parasite responds by decreasing its protein production to conserve nutrient resources and to prevent formation of toxic proteins. However, even during severe starvation, the parasite generates certain essential quality control and rescue proteins. Differential identification of the complete repertoire of these proteins synthesised during starvation from the pre-existing proteins in the parasite holds the key to understanding the starvation adaptation mechanisms. This has been challenging to accomplish due to technical limitations. Using a combination of chemical labelling techniques and protein mass-spectrometry, we selectively identified and measured the proteins generated in the starving Leishmania parasite. Our results show a starvation time-dependent differential expression of important protein synthesis regulators in the parasite. This will serve as an important dataset for a holistic understanding of the starvation adaptation mechanisms in Leishmania. We also believe that this method will find widespread applications in the field of protozoa and other parasites causing Neglected Tropical Diseases.

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

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

Nucleolar Structure and Function in Trypanosomatid Protozoa

The nucleolus is the conspicuous nuclear body where ribosomal RNA genes are transcribed by RNA polymerase I, pre-ribosomal RNA is processed, and ribosomal subunits are assembled. Other important functions have been attributed to the nucleolus over the years. Here we review the current knowledge about the structure and function of the nucleolus in the trypanosomatid parasites Trypanosoma brucei, Trypanosoma cruzi and Leishmania ssp., which represent one of the earliest branching lineages among the eukaryotes. These protozoan parasites present a single nucleolus that is preserved throughout the closed nuclear division, and that seems to lack fibrillar centers. Trypanosomatids possess a relatively low number of rRNA genes, which encode rRNA molecules that contain large expansion segments, including several that are trypanosomatid-specific. Notably, the large subunit rRNA (28S-type) is fragmented into two large and four small rRNA species. Hence, compared to other organisms, the rRNA primary transcript requires additional processing steps in trypanosomatids. Accordingly, this group of parasites contains the highest number ever reported of snoRNAs that participate in rRNA processing. The number of modified rRNA nucleotides in trypanosomatids is also higher than in other organisms. Regarding the structure and biogenesis of the ribosomes, recent cryo-electron microscopy analyses have revealed several trypanosomatid-specific features that are discussed here. Additional functions of the nucleolus in trypanosomatids are also reviewed.

Nutritional stress targets LeishIF4E-3 to storage granules that contain RNA and ribosome components in Leishmania

Leishmania parasites lack pathways for de novo purine biosynthesis. The depletion of purines induces differentiation into virulent metacyclic forms. In vitro, the parasites can survive prolonged periods of purine withdrawal changing their morphology to long and slender cells with an extended flagellum, and decreasing their translation rates. Reduced translation leads to the appearance of discrete granules that contain LeishIF4E-3, one of the six eIF4E paralogs encoded by the Leishmania genome. We hypothesize that each is responsible for a different function during the life cycle. LeishIF4E-3 is a weak cap-binding protein paralog, but its involvement in translation under normal conditions cannot be excluded. However, in response to nutritional stress, LeishIF4E-3 concentrates in specific cytoplasmic granules. LeishIF4E-3 granulation can be induced by the independent elimination of purines, amino acids and glucose. As these granules contain mature mRNAs, we propose that these bodies store inactive transcripts until recovery from stress occurs. In attempt to examine the content of the nutritional stress-induced granules, they were concentrated over sucrose gradients and further pulled-down by targeting in vivo tagged LeishIF4E-3. Proteomic analysis highlighted granule enrichment with multiple ribosomal proteins, suggesting that ribosome particles are abundant in these foci, as expected in case of translation inhibition. RNA-binding proteins, RNA helicases and metabolic enzymes were also enriched in the granules, whereas no degradation enzymes or P-body markers were detected. The starvation-induced LeishIF4E-3-containing granules, therefore, appear to store stalled ribosomes and ribosomal subunits, along with their associated mRNAs. Following nutritional stress, LeishIF4E-3 becomes phosphorylated at position S75, located in its less-conserved N-terminal extension. The ability of the S75A mutant to form granules was reduced, indicating that cellular signaling regulates LeishIF4E-3 function.

Cells respond to cellular stress by decreasing protein translation, to prevent the formation of partially folded or misfolded new polypeptides whose accumulation can be detrimental to living cells. Under such conditions, the cells benefit from storing inactive mRNAs and stalled ribosomal particles, to maintain their availability once conditions improve; dedicated granules offer a solution for such storage. Leishmania parasites are exposed to a variety of stress conditions as a natural part of their life cycle, including the nutritional stress that the parasites experience within the gut of the sandfly. Thus, Leishmania and related trypanosomatids serve as a good model system to investigate RNA fate during different stress conditions. Various granules appear in Leishmania and related organisms in response to different stress conditions. Here, we investigated how nutritional stress, in particular elimination of purines, induced the formation of granules that harbor a specific cap-binding protein, LeishIF4E-3. The starvation-induced LeishIF4E-3 containing granules consist of a variety of ribosomal proteins, along with RNA-binding proteins and mature mRNAs. We thus propose that Leishmania modulates the assembly of LeishIF4E-3-containing granules for transient storage of stalled ribosomal particles and inactive mRNAs. Following renewal of nutrient availability, as occurs during the parasite’s life cycle, the granules disappear. Although their fate is yet unclear, they could be recycled in the cell. Unlike other granules described in trypanosomes, the LeishIF4E-3-containing granules did not contain RNA degradation enzymes, suggesting that their function is mainly for storage until conditions improve.

Ribosome Profiling Reveals HSP90 Inhibitor Effects on Stage-Specific Protein Synthesis in Leishmania donovani

Leishmania parasites cause severe illness in humans and animals. They exist in two developmental stages, insect form and mammalian form, which differ in shape and gene expression. By mapping and quantifying RNA fragments protected by protein synthesis complexes, we determined the rates of protein synthesis for >90% of all Leishmania proteins in response to the inhibition of a key regulatory protein, the 90-kDa heat shock protein. We find that Leishmania depends on a regulation of protein synthesis for controlling its gene expression and that heat shock protein 90 inhibition can trigger the developmental program from insect form to mammalian form of the pathogen.

The 90-kDa heat shock protein (HSP90) of eukaryotes is a highly abundant and essential chaperone required for the maturation of regulatory and signal proteins. In the protozoan parasite Leishmania donovani, causative agent of the fatal visceral leishmaniasis, HSP90 activity is essential for cell proliferation and survival. Even more importantly, its inhibition causes life cycle progression from the insect stage to the pathogenic, mammalian stage. To unravel the molecular impact of HSP90 activity on the parasites’ gene expression, we performed a ribosome profiling analysis of L. donovani, comparing genome-wide protein synthesis patterns in the presence and absence of the HSP90-specific inhibitor radicicol and an ectopically expressed radicicol-resistant HSP90 variant. We find that ribosome-protected RNA faithfully maps open reading frames and represents 97% of the annotated protein-coding genes of L. donovani. Protein synthesis was found to correlate poorly with RNA steady-state levels, indicating a regulated translation as primary mechanism for HSP90-dependent gene expression. The results confirm inhibitory effects of HSP90 on the synthesis of Leishmania proteins that are associated with the pathogenic, intracellular stage of the parasite. Those include heat shock proteins, redox enzymes, virulence-enhancing surface proteins, proteolytic pathways, and a complete set of histones. Conversely, HSP90 promotes fatty acid synthesis enzymes. Complementing radicicol treatment with the radicicol-resistant HSP90rr variant revealed important off-target radicicol effects that control a large number of the above-listed proteins. Leishmania lacks gene-specific transcription regulation and relies on regulated translation instead. Our ribosome footprinting analysis demonstrates a controlling function of HSP90 in stage-specific protein synthesis but also significant, HSP90-independent effects of the inhibitor radicicol.

IMPORTANCELeishmania parasites cause severe illness in humans and animals. They exist in two developmental stages, insect form and mammalian form, which differ in shape and gene expression. By mapping and quantifying RNA fragments protected by protein synthesis complexes, we determined the rates of protein synthesis for >90% of all Leishmania proteins in response to the inhibition of a key regulatory protein, the 90-kDa heat shock protein. We find that Leishmania depends on a regulation of protein synthesis for controlling its gene expression and that heat shock protein 90 inhibition can trigger the developmental program from insect form to mammalian form of the pathogen.

Fluorescence-based assay for accurate measurement of transcriptional activity in trypanosomatid parasites

Trypanosomatid parasites are causative agents of neglected human diseases. Their lineage diverged early from the common eukaryotic ancestor, and they evolved singular mechanisms of gene expression that are crucial for their survival. Studies on unusual and essential molecular pathways lead to new drug targets. In this respect, assays to analyze transcriptional activity will provide useful information to identify essential and specific factors. However, the current methods are laborious and do not provide global and accurate measures. For this purpose, a previously reported radiolabeling in vitro nascent mRNA methodology was used to establish an alternative fluorescent-based assay that is able to precisely quantify nascent mRNA using both flow cytometry and a high-content image system. The method allowed accurate and global measurements in Trypanosoma brucei, a representative species of trypanosomatid parasites. We obtained data demonstrating that approximately 70% of parasites from a population under normal growth conditions displayed mRNA transcriptional activity, whilst the treatment with α-amanitin (75 µg/ml) inhibited the polymerase II activity. The adaptation of the method also allowed the analyses of the transcriptional activity during the cell cycle. Therefore, the methodology described herein contributes to obtaining precise measurements of transcriptional rates using multiparametric analysis. This alternative method can facilitate investigations of genetic and biochemical processes in trypanosome parasites and consequently provide additional information related to new treatment or prophylaxis strategies involving these important human parasites.

Distinct 3' UTRs regulate the life-cycle-specific expression of two TCTP paralogs in Trypanosoma brucei

The translationally controlled tumor protein (TCTP; also known as TPT1 in mammals) is highly conserved and ubiquitously expressed in eukaryotes. It is involved in growth and development, cell cycle progression, protection against cellular stresses and apoptosis, indicating the multifunctional role of the protein. Here, for the first time, we characterize the expression and function of TCTP in the human and animal pathogen, Trypanosoma brucei. We identified two paralogs (TCTP1 and TCTP2) that are differentially expressed in the life cycle of the parasite. The genes have identical 5′ untranslated regions (UTRs) and almost identical open-reading frames. The 3′UTRs differ substantially in sequence and length, and are sufficient for the exclusive expression of TCTP1 in procyclic- and TCTP2 in bloodstream-form parasites. Furthermore, we characterize which parts of the 3′UTR are needed for TCTP2 mRNA stability. RNAi experiments demonstrate that TCTP1 and TCTP2 expression is essential for normal cell growth in procyclic- and bloodstream-form parasites, respectively. Depletion of TCTP1 in the procyclic form cells leads to aberrant cell and mitochondrial organelle morphology, as well as enlarged, and a reduced number of, acidocalcisomes.

Summary:T. brucei has two TCTP genes that are differentially expressed during the parasite life cycle owing to their different 3′UTRs. TCTP also has a role in regulating cell growth and morphology.

Gene and Chromosomal Copy Number Variations as an Adaptive Mechanism Towards a Parasitic Lifestyle in Trypanosomatids

Trypanosomatids are a group of kinetoplastid parasites including some of great public health importance, causing debilitating and life-long lasting diseases that affect more than 24 million people worldwide. Among the trypanosomatids, Trypanosoma cruzi, Trypanosoma brucei and species from the Leishmania genus are the most well studied parasites, due to their high prevalence in human infections. These parasites have an extreme genomic and phenotypic variability, with a massive expansion in the copy number of species-specific multigene families enrolled in host-parasite interactions that mediate cellular invasion and immune evasion processes. As most trypanosomatids are heteroxenous, and therefore their lifecycles involve the transition between different hosts, these parasites have developed several strategies to ensure a rapid adaptation to changing environments. Among these strategies, a rapid shift in the repertoire of expressed genes, genetic variability and genome plasticity are key mechanisms. Trypanosomatid genomes are organized into large directional gene clusters that are transcribed polycistronically, where genes derived from the same polycistron may have very distinct mRNA levels. This particular mode of transcription implies that the control of gene expression operates mainly at post-transcriptional level. In this sense, gene duplications/losses were already associated with changes in mRNA levels in these parasites. Gene duplications also allow the generation of sequence variability, as the newly formed copy can diverge without loss of function of the original copy. Recently, aneuploidies have been shown to occur in several Leishmania species and T. cruzi strains. Although aneuploidies are usually associated with debilitating phenotypes in superior eukaryotes, recent data shows that it could also provide increased fitness in stress conditions and generate drug resistance in unicellular eukaryotes. In this review, we will focus on gene and chromosomal copy number variations and their relevance to the evolution of trypanosomatid parasites.

The genomes of Crithidia bombi and C. expoeki, common parasites of bumblebees

Trypanosomatids (Trypanosomatidae, Kinetoplastida) are flagellated protozoa containing many parasites of medical or agricultural importance. Among those, Crithidia bombi and C. expoeki, are common parasites in bumble bees around the world, and phylogenetically close to Leishmania and Leptomonas. They have a simple and direct life cycle with one host, and partially castrate the founding queens greatly reducing their fitness. Here, we report the nuclear genome sequences of one clone of each species, extracted from a field-collected infection. Using a combination of Roche 454 FLX Titanium, Pacific Biosciences PacBio RS, and Illumina GA2 instruments for C. bombi, and PacBio for C. expoeki, we could produce high-quality and well resolved sequences. We find that these genomes are around 32 and 34 MB, with 7,808 and 7,851 annotated genes for C. bombi and C. expoeki, respectively-which is somewhat less than reported from other trypanosomatids, with few introns, and organized in polycistronic units. A large fraction of genes received plausible functional support in comparison primarily with Leishmania and Trypanosoma. Comparing the annotated genes of the two species with those of six other trypanosomatids (C. fasciculata, L. pyrrhocoris, L. seymouri, B. ayalai, L. major, and T. brucei) shows similar gene repertoires and many orthologs. Similar to other trypanosomatids, we also find signs of concerted evolution in genes putatively involved in the interaction with the host, a high degree of synteny between C. bombi and C. expoeki, and considerable overlap with several other species in the set. A total of 86 orthologous gene groups show signatures of positive selection in the branch leading to the two Crithidia under study, mostly of unknown function. As an example, we examined the initiating glycosylation pathway of surface components in C. bombi, finding it deviates from most other eukaryotes and also from other kinetoplastids, which may indicate rapid evolution in the extracellular matrix that is involved in interactions with the host. Bumble bees are important pollinators and Crithidia-infections are suspected to cause substantial selection pressure on their host populations. These newly sequenced genomes provide tools that should help better understand host-parasite interactions in these pollinator pathogens.

The Pumilio-domain protein PUF6 contributes to SIDER2 retroposon-mediated mRNA decay in Leishmania

Leishmania and other trypanosomatid protozoa lack control at the level of transcription initiation and regulate gene expression exclusively post-transcriptionally. We have reported previously that Leishmania harbors a unique class of short interspersed degenerate retroposons (SIDERs) that are predominantly located within 3'UTRs and play a major role in post-transcriptional control. We have shown that members of the SIDER2 subfamily initiate mRNA decay through endonucleolytic cleavage within the second conserved 79-nt signature sequence of SIDER2 retroposons. Here, we have developed an optimized MS2 coat protein tethering system to capture trans-acting factor(s) regulating SIDER2-mediated mRNA decay. Tethering of the MS2 coat protein to a reporter RNA harboring two MS2 stem-loop aptamers and the cognate SIDER2-containing 3'UTR in combination with immunoprecipitation and mass spectrometry analysis led to the identification of RNA-binding proteins with known functions in mRNA decay. Among the candidate SIDER2-interacting proteins that were individually tethered to a SIDER2 reporter RNA, the Pumilio-domain protein PUF6 was shown to enhance degradation and reduce transcript half-life. Furthermore, we showed that PUF6 binds to SIDER2 sequences that include the regulatory 79-nt signature motif, hence contributing to the mRNA decay process. Consistent with a role of PUF6 in SIDER2-mediated decay, genetic inactivation of PUF6 resulted in increased accumulation and higher stability of endogenous SIDER2-bearing transcripts. Overall, these studies provide new insights into regulated mRNA decay pathways in Leishmania controlled by SIDER2 retroposons and propose a broader role for PUF proteins in mRNA decay within the eukaryotic kingdom.

RNA secondary structure and nucleotide composition of the conserved hallmark sequence of Leishmania SIDER2 retroposons are essential for endonucleolytic cleavage and mRNA degradation

We have reported previously that Short Interspersed Degenerate Retroposons of the SIDER2 subfamily, largely located within 3'UTRs of Leishmania transcripts, promote rapid turnover of mRNAs through endonucleolytic cleavage within the highly conserved second tandem 79-nt hallmark sequence (79-nt SII). Here, we used site-directed mutagenesis and in silico RNA structural studies to delineate the cis-acting requirements within 79-nt SII for cleavage and mRNA degradation. The putative cleavage site(s) and other nucleotides predicted to alter the RNA secondary structure of 79-nt SII were either deleted or mutated and their effect on mRNA turnover was monitored using a gene reporter system. We found that short deletions of 8-nt spanning the two predicted cleavage sites block degradation of SIDER2-containing transcripts, leading to mRNA accumulation. Furthermore, single or double substitutions of the dinucleotides targeted for cleavage as well as mutations altering the predicted RNA secondary structure encompassing both cleavage sites also prevent mRNA degradation, confirming that these dinucleotides are the bona fide cleavage sites. In line with these results, we show that stage-regulated SIDER2 inactivation correlates with the absence of endonucleolytic cleavage. Overall, these data demonstrate that both cleavage sites within the conserved 79-nt SII as well as RNA folding in this region are essential for SIDER2-mediated mRNA decay, and further support that SIDER2-harboring transcripts are targeted for degradation by endonucleolytic cleavage.

SIDER2 retroposon-mediated mRNA decay in Leishmania is coupled to translation

We previously reported that Short Interspersed Degenerate Retroposons of the SIDER2 subfamily predominantly located within 3' untranslated regions (UTRs) of Leishmania transcripts promote rapid turnover that is initiated by endonucleolytic cleavage. Here, we investigated whether SIDER2-mediated mRNA decay is linked to translation. We show that preventing translation initiation by inserting a hairpin structure at the 5'-end of a SIDER2-containing mRNA blocks degradation. Similarly, global inhibition of translation elongation by cycloheximide or termination by puromycin causes stabilisation of SIDER2-containing transcripts. Altogether, these findings support that the mechanism of SIDER2-mediated decay is coupled to translation, possibly through the recruitment of decay factors to elongating ribosomes.

Nuclear Compartmentalization Contributes to Stage-Specific Gene Expression Control in Trypanosoma cruzi

In the protozoan parasite Trypanosoma cruzi, as in other trypanosomatids, transcription of protein coding genes occurs in a constitutive fashion, producing large polycistronic transcription units. These units are composed of non-functionally related genes which are pervasively processed to yield each mRNA. Therefore, post-transcriptional processes are crucial to regulate gene expression. Considering that nuclear compartmentalization could contribute to gene expression regulation, we comparatively studied the nuclear, cytoplasmic and whole cell transcriptomes of the non-infective epimastigote stage of T. cruzi, using RNA-Seq. We found that the cytoplasmic transcriptome tightly correlates with the whole cell transcriptome and both equally correlate with the proteome. Nonetheless, 1,200 transcripts showed differential abundance between the nuclear and cytoplasmic fractions. For the genes with transcript content augmented in the nucleus, significant structural and compositional differences were found. The analysis of the reported epimastigote translatome and proteome, revealed scarce ribosome footprints and encoded proteins for them. Ontology analyses unveiled that many of these genes are distinctive of other parasite life-cycle stages. Finally, the relocalization of transcript abundance in the metacyclic trypomastigote infective stage was confirmed for specific genes. While gene expression is strongly dependent on transcript steady-state level, we here highlight the importance of the distribution of transcripts abundance between compartments in T. cruzi. Particularly, we show that nuclear compartmentation is playing an active role in the developmental stage determination preventing off-stage expression.

RNA-seq in kinetoplastids: A powerful tool for the understanding of the biology and host-pathogen interactions

The kinetoplastids include a large number of parasites responsible for serious diseases in humans and animals (Leishmania and Trypanosoma brucei) considered endemic in several regions of the world. These parasites are characterized by digenetic life cycles that undergo morphological and genetic changes that allow them to adapt to different microenvironments on their vertebrates and invertebrates hosts. Recent advances in ´omics´ technology, specifically transcriptomics have allowed to reveal aspects associated with such molecular changes. So far, different techniques have been used to evaluate the gene expression profile during the various stages of the life cycle of these parasites and during the host-parasite interactions. However, some of them have serious drawbacks that limit the precise study and full understanding of their transcriptomes. Therefore, recently has been implemented the latest technology (RNA-seq), which overcomes the drawbacks of traditional methods. In this review, studies that so far have used RNA-seq are presented and allowed to expand our knowledge regarding the biology of these parasites and their interactions with their hosts.

Phosphorylation of Translation Initiation Factor 2-Alpha in Leishmania donovani under Stress Is Necessary for Parasite Survival

The transformation of Leishmania donovani from a promastigote to an amastigote during mammalian host infection displays the immense adaptability of the parasite to survival under stress. Induction of translation initiation factor 2-alpha (eIF2α) phosphorylation by stress-specific eIF2α kinases is the basic stress-perceiving signal in eukaryotes to counter stress. Here, we demonstrate that elevated temperature and acidic pH induce the phosphorylation of Leishmania donovani eIF2α (LdeIF2α). In vitro inhibition experiments suggest that interference of LdeIF2α phosphorylation under conditions of elevated temperature and acidic pH debilitates parasite differentiation and reduces parasite viability (P < 0.05). Furthermore, inhibition of LdeIF2α phosphorylation significantly reduced the infection rate (P < 0.05), emphasizing its deciding role in successful invasion and infection establishment. Notably, our findings suggested the phosphorylation of LdeIF2α under H₂O₂-induced oxidative stress. Inhibition of H₂O₂-induced LdeIF2α phosphorylation hampered antioxidant balance by impaired redox homeostasis gene expression, resulting in increased reactive oxygen species accumulation (P < 0.05) and finally leading to decreased parasite viability (P < 0.05). Interestingly, exposure to sodium antimony glucamate and amphotericin B induces LdeIF2α phosphorylation, indicating its possible contribution to protection against antileishmanial drugs in common use. Overall, the results strongly suggest that stress-induced LdeIF2α phosphorylation is a necessary event for the parasite life cycle under stressed conditions for survival.

The Role of the Trypanosoma cruzi TcNRBD1 Protein in Translation

The regulation of gene expression in trypanosomatids occurs mainly at the post-transcriptional level. Despite the importance of this type of control in Trypanosoma cruzi, few RNA binding proteins have been characterized. The RRM domain (RNA Recognition Motif) is one of the most abundant domains found in RNA-binding proteins in higher eukaryotes. Proteins containing the RRM domain are involved in the majority of post-transcriptional processes regulating gene expression. In this work, we aimed to characterize the protein TcNRBD1 from T. cruzi. TcNRBD1 is an RNA-binding protein that contains 2 RRM domains and is the ortholog of the P34 and P37 proteins from Trypanosoma brucei. The TcNRBD1 protein is expressed in all developmental stages of T. cruzi, and its localization pattern is concentrated at the perinuclear region. TcNRBD1 is associated with polysomes and with the 80S monosomes. Furthermore, sequencing of the mRNAs bound to TcNRBD1 allowed the identification of several transcripts that encode ribosomal proteins. Immunoprecipitation assays followed by mass spectrometry showed that the protein complexes with several ribosomal proteins from both the 40S and 60S subunits. In summary, the results indicate that TcNRBD1 is associated with different parts of the translation process, either by regulating mRNAs that encode ribosomal proteins or by acting in some step of ribosome assembly in T. cruzi.

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

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

The evolutionary landscape of intergenic trans-splicing events in insects

To explore the landscape of intergenic trans-splicing events and characterize their functions and evolutionary dynamics, we conduct a mega-data study of a phylogeny containing eight species across five orders of class Insecta, a model system spanning 400 million years of evolution. A total of 1,627 trans-splicing events involving 2,199 genes are identified, accounting for 1.58% of the total genes. Homology analysis reveals that mod(mdg4)-like trans-splicing is the only conserved event that is consistently observed in multiple species across two orders, which represents a unique case of functional diversification involving trans-splicing. Thus, evolutionarily its potential for generating proteins with novel function is not broadly utilized by insects. Furthermore, 146 non-mod trans-spliced transcripts are found to resemble canonical genes from different species. Trans-splicing preserving the function of 'breakup' genes may serve as a general mechanism for relaxing the constraints on gene structure, with profound implications for the evolution of genes and genomes.

The unique Leishmania EIF4E4 N-terminus is a target for multiple phosphorylation events and participates in critical interactions required for translation initiation

The eukaryotic initiation factor 4E (eIF4E) recognizes the mRNA cap structure and, together with eIF4G and eIF4A, form the eIF4F complex that regulates translation initiation in eukaryotes. In trypanosomatids, 2 eIF4E homologues (EIF4E3 and EIF4E4) have been shown to be part of eIF4F-like complexes with presumed roles in translation initiation. Both proteins possess unique N-terminal extensions, which can be targeted for phosphorylation. Here, we provide novel insights on the Leishmania infantum EIF4E4 function and regulation. We show that EIF4E4 is constitutively expressed throughout the parasite development but is preferentially phosphorylated in exponentially grown promastigote and amastigote life stages, hence correlating with high levels of translation. Phosphorylation targets multiple serine-proline or threonine-proline residues within the N-terminal extension of EIF4E4 but does not require binding to the EIF4E4's partner, EIF4G3, or to the cap structure. We also report that EIF4E4 interacts with PABP1 through 3 conserved boxes at the EIF4E4 N-terminus and that this interaction is a prerequisite for efficient EIF4E4 phosphorylation. EIF4E4 is essential for Leishmania growth and an EIF4E4 null mutant was only obtained in the presence of an ectopically provided wild type gene. Complementation for the loss of EIF4E4 with several EIF4E4 mutant proteins affecting either phosphorylation or binding to mRNA or to EIF4E4 protein partners revealed that, in contrast to other eukaryotes, only the EIF4E4-PABP1 interaction but neither the binding to EIF4G3 nor phosphorylation is essential for translation. These studies also demonstrated that the lack of both EIF4E4 phosphorylation and EIF4G3 binding leads to a non-functional protein. Altogether, these findings further highlight the unique features of the translation initiation process in trypanosomatid protozoa.

Differential Subcellular Localization of Leishmania Alba-Domain Proteins throughout the Parasite Development

Alba-domain proteins are RNA-binding proteins found in archaea and eukaryotes and recently studied in protozoan parasites where they play a role in the regulation of virulence factors and stage-specific proteins. This work describes in silico structural characterization, cellular localization and biochemical analyses of Alba-domain proteins in Leishmania infantum. We show that in contrast to other protozoa, Leishmania have two Alba-domain proteins, LiAlba1 and LiAlba3, representative of the Rpp20- and the Rpp25-like eukaryotic subfamilies, respectively, which share several sequence and structural similarities but also important differences with orthologs in other protozoa, especially in sequences targeted for post-translational modifications. LiAlba1 and LiAlba3 proteins form a complex interacting with other RNA-binding proteins, ribosomal subunits, and translation factors as supported by co-immunoprecipitation and sucrose gradient sedimentation analysis. A higher co-sedimentation of Alba proteins with ribosomal subunits was seen upon conditions of decreased translation, suggesting a role of these proteins in translational repression. The Leishmania Alba-domain proteins display differential cellular localization throughout the parasite development. In the insect promastigote stage, Alba proteins co-localize predominantly to the cytoplasm but they translocate to the nucleolus and the flagellum upon amastigote differentiation in the mammalian host and are found back to the cytoplasm once amastigote differentiation is completed. Heat-shock, a major signal of amastigote differentiation, triggers Alba translocation to the nucleolus and the flagellum. Purification of the Leishmania flagellum confirmed LiAlba3 enrichment in this organelle during amastigote differentiation. Moreover, partial characterization of the Leishmania flagellum proteome of promastigotes and differentiating amastigotes revealed the presence of other RNA-binding proteins, as well as differences in the flagellum composition between these two parasite lifestages. Shuttling of Alba-domain proteins between the cytoplasm and the nucleolus or the flagellum throughout the parasite life cycle suggests that these RNA-binding proteins participate in several distinct regulatory pathways controlling developmental gene expression in Leishmania.

The eIF3 complex of Leishmania-subunit composition and mode of recruitment to different cap-binding complexes

Eukaryotic initiation factor 3 (eIF3) is a multi-protein complex and a key participant in the assembly of the translation initiation machinery. In mammals, eIF3 comprises 13 subunits, most of which are characterized by conserved structural domains. The trypanosomatid eIF3 subunits are poorly conserved. Here, we identify 12 subunits that comprise the Leishmania eIF3 complex (LeishIF3a-l) by combining bioinformatics with affinity purification and mass spectrometry analyses. These results highlight the strong association of LeishIF3 with LeishIF1, LeishIF2 and LeishIF5, suggesting the existence of a multi-factor complex. In trypanosomatids, the translation machinery is tightly regulated in the different life stages of these organisms as part of their adaptation and survival in changing environments. We, therefore, addressed the mechanism by which LeishIF3 is recruited to different mRNA cap-binding complexes. A direct interaction was observed in vitro between the fully assembled LeishIF3 complex and recombinant LeishIF4G3, the canonical scaffolding protein of the cap-binding complex in Leishmania promastigotes. We further highlight a novel interaction between the C-terminus of LeishIF3a and LeishIF4E1, the only cap-binding protein that efficiently binds the cap structure under heat shock conditions, anchoring a complex that is deficient of any MIF4G-based scaffolding subunit.

LiZIP3 is a cellular zinc transporter that mediates the tightly regulated import of zinc in Leishmania infantum parasites

Cellular zinc homeostasis ensures that the intracellular concentration of this element is kept within limits that enable its participation in critical physiological processes without exerting toxic effects. We report here the identification and characterization of the first mediator of zinc homeostasis in Leishmania infantum, LiZIP3, a member of the ZIP family of divalent metal-ion transporters. The zinc transporter activity of LiZIP3 was first disclosed by its capacity to rescue the growth of Saccharomyces cerevisiae strains deficient in zinc acquisition. Subsequent expression of LiZIP3 in Xenopus laevis oocytes was shown to stimulate the uptake of a broad range of metal ions, among which Zn(2+) was the preferred LiZIP3 substrate (K0.5  ≈ 0.1 μM). Evidence that LiZIP3 functions as a zinc importer in L. infantum came from the observations that the protein locates to the cell membrane and that its overexpression leads to augmented zinc internalization. Importantly, expression and cell-surface location of LiZIP3 are lost when parasites face high zinc bioavailability. LiZIP3 decline in response to zinc is regulated at the mRNA level in a process involving (a) short-lived protein(s). Collectively, our data reveal that LiZIP3 enables L. infantum to acquire zinc in a highly regulated manner, hence contributing to zinc homeostasis.

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

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

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

The genome of Leishmania panamensis: insights into genomics of the L. (Viannia) subgenus

Kinetoplastid parasites of the Leishmania genus cause several forms of leishmaniasis. Leishmania species pathogenic to human are separated into two subgenera, Leishmania (Leishmania) and L. (Viannia). Species from the Viannia subgenus cause predominantly cutaneous leishmaniasis in Central and South America, occasionally leading to more severe clinical presentations. Although the genomes of several species of Leishmania have been sequenced to date, only one belongs to this rather different subgenus. Here we explore the unique features of the Viannia subgenus by sequencing and analyzing the genome of L. (Viannia) panamensis. Against a background of conservation in gene content and synteny, we found key differences at the genomic level that may explain the occurrence of molecular processes involving nucleic acid manipulation and differential modification of surface glycoconjugates. These differences may in part explain some phenotypic characteristics of the Viannia parasites, including their increased adaptive capacity and enhanced metastatic ability.