Post-transcriptional regulation of gene expression in trypanosomes and leishmanias

Integrating high-throughput analysis to create an atlas of replication origins in Trypanosoma cruzi in the context of genome structure and variability

Trypanosoma cruzi is the etiologic agent of the most prevalent human parasitic disease in Latin America, Chagas disease. Its genome is rich in multigenic families that code for virulent antigens and are present in the rapidly evolving genomic compartment named Disruptive. DNA replication is a meticulous biological process in which flaws can generate mutations and changes in chromosomal and gene copy numbers. Here, integrating high-throughput and single-molecule analyses, we were able to identify Predominant, Flexible, and Dormant Orc1Cdc6-dependent origins as well as Orc1Cdc6-independent origins. Orc1Cdc6-dependent origins were found in multigenic family loci, while independent origins were found in the Core compartment that contains conserved and hypothetical protein-coding genes, in addition to multigenic families. In addition, we found that Orc1Cdc6 density is related to the firing of origins and that Orc1Cdc6-binding sites within fired origins are depleted of a specific class of nucleosomes that we previously categorized as dynamic. Together, these data suggest that Orc1Cdc6-dependent origins may contribute to the rapid evolution of the Disruptive compartment and, therefore, to the success of T. cruzi infection and that the local epigenome landscape is also involved in this process.IMPORTANCETrypanosoma cruzi, responsible for Chagas disease, affects millions globally, particularly in Latin America. Lack of vaccine or treatment underscores the need for research. Parasite's genome, with virulent antigen-coding multigenic families, resides in the rapidly evolving Disruptive compartment. Study sheds light on the parasite's dynamic DNA replication, discussing the evolution of the Disruptive compartment. Therefore, the findings represent a significant stride in comprehending T. cruzi's biology and the molecular bases that contribute to the success of infection caused by this parasite.

A conserved trypanosomatid differentiation regulator controls substrate attachment and morphological development in Trypanosoma congolense

Trypanosomatid parasites undergo developmental regulation to adapt to the different environments encountered during their life cycle. In Trypanosoma brucei, a genome wide selectional screen previously identified a regulator of the protein family ESAG9, which is highly expressed in stumpy forms, a morphologically distinct bloodstream stage adapted for tsetse transmission. This regulator, TbREG9.1, has an orthologue in Trypanosoma congolense, despite the absence of a stumpy morphotype in that parasite species, which is an important cause of livestock trypanosomosis. RNAi mediated gene silencing of TcREG9.1 in Trypanosoma congolense caused a loss of attachment of the parasites to a surface substrate in vitro, a key feature of the biology of these parasites that is distinct from T. brucei. This detachment was phenocopied by treatment of the parasites with a phosphodiesterase inhibitor, which also promotes detachment in the insect trypanosomatid Crithidia fasciculata. RNAseq analysis revealed that TcREG9.1 silencing caused the upregulation of mRNAs for several classes of surface molecules, including transferrin receptor-like molecules, immunoreactive proteins in experimental bovine infections, and molecules related to those associated with stumpy development in T. brucei. Depletion of TcREG9.1 in vivo also generated an enhanced level of parasites in the blood circulation consistent with reduced parasite attachment to the microvasculature. The morphological progression to insect forms of the parasite was also perturbed. We propose a model whereby TcREG9.1 acts as a regulator of attachment and development, with detached parasites being adapted for transmission.

Structural analysis of the Trypanosoma brucei EIF4E6/EIF4G5 complex reveals details of the interaction between unusual eIF4F subunits

Recognition of the mRNA 5' end is a critical step needed for translation initiation. This step is performed by the cap binding protein eIF4E, which joins the larger eIF4G subunit to form the eIF4F complex. Trypanosomatids have a minimum of five different eIF4F-like complexes formed through specific but not well-defined interactions between four different eIF4E and five eIF4G homologues. The EIF4E6/EIF4G5 complex has been linked with the stage-specific translation of mRNAs encoding the major Trypanosoma brucei virulence factors. Here, to better define the molecular basis for the TbEIF4E6/TbEIF4G5 interaction, we describe the identification of the peptide interacting with TbEIF4E6 in the region comprising residues 79-166 of TbEIF4G5. The TbEIF4E6-TbEIF4G5_79-116 complex reconstituted with recombinant proteins is highly stable even in the absence of cap-4. The crystal structure of the complex was subsequently solved, revealing extensive interacting surfaces. Comparative analyses highlight the conservation of the overall structural arrangement of different eIF4E/eIF4G complexes. However, highly different interacting surfaces are formed with distinct binding contacts occurring both in the canonical and noncanonical elements within eIF4G and the respective eIF4E counterpart. These specific pairs of complementary interacting surfaces are likely responsible for the selective association needed for the formation of distinct eIF4F complexes in trypanosomatids.

LeishIF4E2 is a cap-binding protein that plays a role in Leishmania cell cycle progression

Leishmania encode six paralogs of the cap-binding protein eIF4E and five eIF4G candidates, forming unique complexes. Two cap-binding proteins, LeishIF4E1 and LeishIF4E2, do not bind any identified LeishIF4Gs, thus their roles are intriguing. Here, we combine structural prediction, proteomic analysis, and interaction assays to shed light on LeishIF4E2 function. A nonconserved C-terminal extension was identified through structure prediction and sequence alignment. m7 GTP-binding assays involving both recombinant and transgenic LeishIF4E2 with and without the C-terminal extension revealed that this extension functions as a regulatory gate, modulating the cap-binding activity of LeishIF4E2. The interactomes of the two LeishIF4E2 versions were investigated, highlighting the role of the C-terminal extension in binding to SLBP2. SLBP2 is known to interact with a stem-loop structure in the 3' UTRs of histone mRNAs. Consistent with the predicted inhibitory effect of SLBP2 on histone expression in Xenopus laevis, a hemizygous deletion mutant of LeishIF4E2, exhibited an upregulation of several histones. We therefore propose that LeishIF4E2 is involved in histone expression, possibly through its interaction between SLBP2 and LeishIF4E2, thus affecting cell cycle progression. In addition, cell synchronization showed that LeishIF4E2 expression decreased during the S-phase, when histones are known to be synthesized. Previous studies in T. brucei also highlighted an association between TbEIF4E2 and SLBP2, and further reported on an interaction between TbIF4E2 and S-phase-abundant mRNAs. Our results show that overexpression of LeishIF4E2 correlates with upregulation of cell cycle and chromosome maintenance proteins. Along with its effect on histone expression, we propose that LeishIF4E2 is involved in cell cycle progression.

Genome-wide chromatin interaction map for Trypanosoma cruzi

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

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.

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

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

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.

Multi-omics analysis to characterize molecular adaptation of Entamoeba histolytica during serum stress

Entamoeba histolytica is responsible for dysentery and extraintestinal disease in humans. To establish successful infection, it must generate adaptive response against stress due to host defense mechanisms. We have developed a robust proteomics workflow by combining miniaturized sample preparation, low flow-rate chromatography, and ultra-high sensitivity mass spectrometry, achieving increased proteome coverage, and further integrated proteomics and RNA-seq data to decipher regulation at translational and transcriptional levels. Label-free quantitative proteomics led to identification of 2344 proteins, an improvement over the maximum number identified in E. histolytica proteomic studies. In serum-starved cells, 127 proteins were differentially abundant and were associated with functions including antioxidant activity, cytoskeleton, translation, catalysis, and transport. The virulence factor, Gal/GalNAc-inhibitable lectin subunits, was significantly altered. Integration of transcriptomic and proteomic data revealed that only 30% genes were coordinately regulated at both transcriptional and translational levels. Some highly expressed transcripts did not change in protein abundance. Conversely, genes with no transcriptional change showed enhanced protein abundance, indicating post-transcriptional regulation. This multi-omics approach enables more refined gene expression analysis to understand the adaptive response of E. histolytica during growth stress.

In silico analysis of the HSP90 chaperone system from the African trypanosome, Trypanosoma brucei

African trypanosomiasis is a neglected tropical disease caused by Trypanosoma brucei (T. brucei) and spread by the tsetse fly in sub-Saharan Africa. The trypanosome relies on heat shock proteins for survival in the insect vector and mammalian host. Heat shock protein 90 (HSP90) plays a crucial role in the stress response at the cellular level. Inhibition of its interactions with chaperones and co-chaperones is being explored as a potential therapeutic target for numerous diseases. This study provides an in silico overview of HSP90 and its co-chaperones in both T. brucei brucei and T. brucei gambiense in relation to human and other trypanosomal species, including non-parasitic Bodo saltans and the insect infecting Crithidia fasciculata. A structural analysis of T. brucei HSP90 revealed differences in the orientation of the linker and C-terminal domain in comparison to human HSP90. Phylogenetic analysis displayed the T. brucei HSP90 proteins clustering into three distinct groups based on subcellular localizations, namely, cytosol, mitochondria, and endoplasmic reticulum. Syntenic analysis of cytosolic HSP90 genes revealed that T. b. brucei encoded for 10 tandem copies, while T. b. gambiense encoded for three tandem copies; Leishmania major (L. major) had the highest gene copy number with 17 tandem copies. The updated information on HSP90 from recently published proteomics on T. brucei was examined for different life cycle stages and subcellular localizations. The results show a difference between T. b. brucei and T. b. gambiense with T. b. brucei encoding a total of twelve putative HSP90 genes, while T. b. gambiense encodes five HSP90 genes. Eighteen putative co-chaperones were identified with one notable absence being cell division cycle 37 (Cdc37). These results provide an updated framework on approaching HSP90 and its interactions as drug targets in the African trypanosome.

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.

Stable Episomal Transfectant Leishmania infantum Promastigotes Over-Expressing the DEVH1 RNA Helicase Gene Down-Regulate Parasite Survival Genes

The compartmentalization of untranslated mRNA molecules in granules occurring in many eukaryotic organisms including trypanosomatids involves the formation of complexes between mRNA molecules and RNA-binding proteins (RBPs). The putative ATP-dependent DEAD/H RNA helicase (DEVH1) from Leishmania infantum (Kinetoplastida: Trypanosomatidae) is one such proteins. The objective of this research is finding differentially expressed genes in a stable episomal transfectant L. infantum promastigote line over-expressing DEVH1 in the stationary phase of growth in axenic culture to get insight into the biological roles of this RNA helicase in the parasite. Interestingly, genes related to parasite survival and virulence factors, such as the hydrophilic surface protein/small hydrophilic endoplasmic reticulum protein (HASP/SHERP) gene cluster, an amastin, and genes related to reactive oxygen species detoxification are down-regulated in DEVH1 transfectant promastigotes.

Unique thiol metabolism in trypanosomatids: Redox homeostasis and drug resistance

Trypanosomatids are mainly responsible for heterogeneous parasitic diseases: Leishmaniasis, Sleeping sickness, and Chagas disease and control of these diseases implicates serious challenges due to the emergence of drug resistance. Redox-active biomolecules are the endogenous substances in organisms, which play important role in the regulation of redox homeostasis. The redox-active substances like glutathione, trypanothione, cysteine, cysteine persulfides, etc., and other inorganic intermediates (hydrogen peroxide, nitric oxide) are very useful as defence mechanism. In the present review, the suitability of trypanothione and other essential thiol molecules of trypanosomatids as drug targets are described in Leishmania and Trypanosoma. We have explored the role of tryparedoxin, tryparedoxin peroxidase, ascorbate peroxidase, superoxide dismutase, and glutaredoxins in the anti-oxidant mechanism and drug resistance. Up-regulation of some proteins in trypanothione metabolism helps the parasites in survival against drug pressure (sodium stibogluconate, Amphotericin B, etc.) and oxidative stress. These molecules accept electrons from the reduced trypanothione and donate their electrons to other proteins, and these proteins reduce toxic molecules, neutralize reactive oxygen, or nitrogen species; and help parasites to cope with oxidative stress. Thus, a better understanding of the role of these molecules in drug resistance and redox homeostasis will help to target metabolic pathway proteins to combat Leishmaniasis and trypanosomiases.

Low Dose Gamma Irradiation of Trypanosoma evansi Parasites Identifies Molecular Changes That Occur to Repair Radiation Damage and Gene Transcripts That May Be Involved in Establishing Disease in Mice Post-Irradiation

The protozoan parasite Trypanosoma evansi is responsible for causing surra in a variety of mammalian hosts and is spread by many vectors over a wide geographical area making it an ideal target for irradiation as a tool to study the initial events that occur during infection. Parasites irradiated at the representative doses 100Gy, 140Gy, and 200Gy were used to inoculate BALB/c mice revealing that parasites irradiated at 200Gy were unable to establish disease in all mice. Cytokine analysis of mice inoculated with 200Gy of irradiated parasites showed significantly lower levels of interleukins when compared to mice inoculated with non-irradiated and 100Gy irradiated parasites. Irradiation also differentially affected the abundance of gene transcripts in a dose-dependent trend measured at 6- and 20-hours post-irradiation with 234, 325, and 484 gene transcripts affected 6 hours post-irradiation for 100Gy-, 140Gy- and 200Gy-irradiated parasites, respectively. At 20 hours post-irradiation, 422, 381, and 457 gene transcripts were affected by irradiation at 100Gy, 140Gy, and 200Gy, respectively. A gene ontology (GO) term analysis was carried out for the three representative doses at 6 hours and 20 hours post-irradiation revealing different processes occurring at 20 hours when compared to 6 hours for 100Gy irradiation. The top ten most significant processes had a negative Z score. These processes fall in significance at 140Gy and even further at 200Gy, revealing that they were least likely to occur at 200Gy, and thus may have been responsible for infection in mice by 100Gy and 140Gy irradiated parasites. When looking at 100Gy irradiated parasites 20 hours post-irradiation processes with a positive Z score, we identified genes that were involved in multiple processes and compared their fold change values at 6 hours and 20 hours. We present these genes as possibly necessary for repair from irradiation damage at 6 hours and suggestive of being involved in the establishment of disease in mice at 20 hours post-irradiation. A potential strategy using this information to develop a whole parasite vaccine is also postulated.

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.

Synthetic biology tools for engineering Goodwin oscillation in Trypanosoma brucei brucei

Kinetoplastid protozoa possess properties that are highly divergent from the mammalian, yeast and bacterial cells more commonly used in synthetic biology and represent a tantalisingly untapped source of bioengineering potential. Trypanosoma brucei brucei (T. b. brucei), an established model organism for studying the Kinetoplastida, is non-pathogenic to humans and provides an interesting test case for establishing synthetic biology in this phylogenetic class. To demonstrate further the tractability of Kinetoplastida to synthetic biology, we sought to construct and demonstrate a Goodwin oscillator, the simplest oscillatory gene network, in T. b. brucei for the first time. We report one completed iteration of the archetypal synthetic biology Design-Build-Test-Learn (DBTL) cycle; firstly, using Ab initio mathematical modelling of the behaviour a theoretical, oscillatory, trypanosomal synthetic gene network (SGN) to inform the design of a plasmid encoding that network. Once assembled, the plasmid was then used to generate a stable transfectant T. b. brucei cell line. To test the performance of the oscillatory SGN, a novel experimental setup was established to capture images of the fluorescent signal from motion-restricted live cells. Data captured were consistent with oscillatory behaviour of the SGN, with cellular fluorescence observed to oscillate with a period of 50 min, with varying amplitude and linear growth trend. This first DBTL cycle establishes a foundation for future cycles in which the SGN design and experimental monitoring setup can be further refined.

Characterization of adenine phosphoribosyltransferase (APRT) activity in Trypanosoma brucei brucei: Only one of the two isoforms is kinetically active

Human African Trypanosomiasis (HAT), also known as sleeping sickness, is a Neglected Tropical Disease endemic to 36 African countries, with approximately 70 million people currently at risk for infection. Current therapeutics are suboptimal due to toxicity, adverse side effects, and emerging resistance. Thus, both effective and affordable treatments are urgently needed. The causative agent of HAT is the protozoan Trypanosoma brucei ssp. Annotation of its genome confirms previous observations that T. brucei is a purine auxotroph. Incapable of de novo purine synthesis, these protozoan parasites rely on purine phosphoribosyltransferases to salvage purines from their hosts for the synthesis of purine monophosphates. Complete and accurate genome annotations in combination with the identification and characterization of the catalytic activity of purine salvage enzymes enables the development of target-specific therapies in addition to providing a deeper understanding of purine metabolism in T. brucei. In trypanosomes, purine phosphoribosyltransferases represent promising drug targets due to their essential and central role in purine salvage. Enzymes involved in adenine and adenosine salvage, such as adenine phosphoribosyltransferases (APRTs, EC 2.4.2.7), are of particular interest for their potential role in the activation of adenine and adenosine-based pro-drugs. Analysis of the T. brucei genome shows two putative aprt genes: APRT1 (Tb927.7.1780) and APRT2 (Tb927.7.1790). Here we report studies of the catalytic activity of each putative APRT, revealing that of the two T. brucei putative APRTs, only APRT1 is kinetically active, thereby signifying a genomic misannotation of Tb927.7.1790 (putative APRT2). Reliable genome annotation is necessary to establish potential drug targets and identify enzymes involved in adenine and adenosine-based pro-drug activation.

Impact of Genetic Diversity and Genome Plasticity of Leishmania spp. in Treatment and the Search for Novel Chemotherapeutic Targets

Leishmaniasis is one of the major public health concerns in Latin America, Africa, Asia, and Europe. The absence of vaccines for human use and the lack of effective vector control programs make chemotherapy the main strategy to control all forms of the disease. However, the high toxicity of available drugs, limited choice of therapeutic agents, and occurrence of drug-resistant parasite strains are the main challenges related to chemotherapy. Currently, only a small number of drugs are available for leishmaniasis treatment, including pentavalent antimonials (Sb^V), amphotericin B and its formulations, miltefosine, paromomycin sulphate, and pentamidine isethionate. In addition to drug toxicity, therapeutic failure of leishmaniasis is a serious concern. The occurrence of drug-resistant parasites is one of the causes of therapeutic failure and is closely related to the diversity of parasites in this genus. Owing to the enormous plasticity of the genome, resistance can occur by altering different metabolic pathways, demonstrating that resistance mechanisms are multifactorial and extremely complex. Genetic variability and genome plasticity cause not only the available drugs to have limitations, but also make the search for new drugs challenging. Here, we examined the biological characteristics of parasites that hinder drug discovery.

High genome plasticity and frequent genetic exchange in Leishmania tropica isolates from Afghanistan, Iran and Syria

Background

The kinetoplastid protozoan Leishmania tropica mainly causes cutaneous leishmaniasis in humans in the Middle East, and relapse or treatment failure after treatment are common in this area. L. tropica’s digenic life cycle includes distinct stages in the vector sandfly and the mammalian host. Sexual reproduction and genetic exchange appear to occur more frequently than in other Leishmania species. Understanding these processes is complicated by chromosome instability during cell division that yields aneuploidy, recombination and heterozygosity. This combination of rare recombination and aneuploid permits may reveal signs of hypothetical parasexual mating, where diploid cells fuse to form a transient tetraploid that undergoes chromosomal recombination and gradual chromosomal loss.

Methodology/principal findings

The genome-wide SNP diversity from 22 L. tropica isolates showed chromosome-specific runs of patchy heterozygosity and extensive chromosome copy number variation. All these isolates were collected during 2007–2017 in Sweden from patients infected in the Middle East and included isolates from a patient possessing two genetically distinct leishmaniasis infections three years apart with no evidence of re-infection. We found differing ancestries on the same chromosome (chr36) across multiple samples: matching the reference genome with few derived alleles, followed by blocks of heterozygous SNPs, and then by clusters of homozygous SNPs with specific recombination breakpoints at an inferred origin of replication. Other chromosomes had similar marked changes in heterozygosity at strand-switch regions separating polycistronic transcriptional units.

Conclusion/significance

These large-scale intra- and inter-chromosomal changes in diversity driven by recombination and aneuploidy suggest multiple mechanisms of cell reproduction and diversification in L. tropica, including mitotic, meiotic and parasexual processes. It underpins the need for more genomic surveillance of Leishmania, to detect emerging hybrids that could spread more widely and to better understand the association between genetic variation and treatment outcome. Furthering our understanding of Leishmania genome evolution and ancestry will aid better diagnostics and treatment for cutaneous leishmaniasis caused by L.tropica in the Middle East.

Cutaneous leishmaniasis is mainly caused by Leishmania tropica in the Middle East, where it is known for treatment failure and a need for prolonged and/or multiple treatments. Several factors affect the clinical presentation and treatment outcome, such as host genetic variability and specific immune response, as well as environmental factors and the vector species. Little is known about the parasite genome and its influence on treatment response. By analysing the genome of 22 isolates of L. tropica, we have revealed extensive genomic variation and a complex population structure with evidence of genetic exchange within and among the isolates, indicating a possible presence of sexual or parasexual mechanisms. Understanding the Leishmania genome better may improve future treatment and better understanding of treatment failure and relapse.

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.

Life Cycle Stage-Specific Accessibility of Leishmania donovani Chromatin at Transcription Start Regions

Leishmania donovani is a parasitic protist that causes the lethal Kala-azar fever in India and East Africa. Gene expression in Leishmania is regulated by gene copy number variation and inducible translation while RNA synthesis initiates at a small number of sites per chromosome and proceeds through polycistronic transcription units, precluding a gene-specific regulation (C. Clayton and M. Shapira, Mol Biochem Parasitol 156:93–101, 2007, https://doi.org/10.1016/j.molbiopara.2007.07.007). Here, we analyze the dynamics of chromatin structure in both life cycle stages of the parasite and find evidence for an additional, epigenetic gene regulation pathway in this early branching eukaryote. The assay for transposase-accessible chromatin using sequencing (ATAC-seq) analysis (J. D. Buenrostro, P. G. Giresi, L. C. Zaba, H. Y. Chang, and W. J. Greenleaf, Nat Methods 10:1213–1218, 2013, https://doi.org/10.1038/nmeth.2688) predominantly shows euchromatin at transcription start regions in fast-growing promastigotes, but mostly heterochromatin in the slowly proliferating amastigotes, the mammalian stage, reflecting a previously shown increase of histone synthesis in the latter stage.

IMPORTANCELeishmania parasites are important pathogens with a global impact and cause poverty-related illness and death. They are devoid of classic cis- and trans-acting transcription regulators but use regulated translation and gene copy number variations to adapt to hosts and environments. In this work, we show that transcription start regions present as open euchromatin in fast-growing insect stages but as less-accessible heterochromatin in the slowly proliferating amastigote stage, indicating an epigenetic control of gene accessibility in this early branching eukaryotic pathogen. This finding should stimulate renewed interest in the control of RNA synthesis in Leishmania and related parasites.

Characterization of the RNA-Binding Protein TcSgn1 in Trypanosoma cruzi

RNA-binding proteins (RBPs) participate in several steps of post-transcriptional regulation of gene expression, such as splicing, messenger RNA transport, mRNA localization, and translation. Gene-expression regulation in trypanosomatids occurs primarily at the post-transcriptional level, and RBPs play important roles in the process. Here, we characterized the RBP TcSgn1, which contains one RNA recognition motif (RRM). TcSgn1 is a close ortholog of yeast Saccharomyces cerevisiae protein ScSgn1, which plays a role in translational regulation in the cytoplasm. We found that TcSgn1 in Trypanosoma cruzi is localized in the nucleus in exponentially growing epimastigotes. By performing immunoprecipitation assays of TcSgn1, we identified hundreds of mRNAs associated with the protein, a significant fraction of them coding for nucleic acids binding, transcription, and endocytosis proteins. In addition, we show that TcSgn1 is capable of interacting directly with the poly(A) tail of the mRNAs. The study of parasites under nutritional stress showed that TcSgn1 was localized in cytoplasmic granules in addition to localizing in the nucleus. Similar to ScSgn1, we observed that TcSgn1 also interacts with the PABP1 protein, suggesting that this protein may play a role in regulating gene expression in T. cruzi. Taken together, our results show that RNA-binding protein TcSgn1 is part of ribonucleoprotein complexes associated with nuclear functions, stress response, and RNA metabolism.

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.

Comparative phosphoproteomic analysis unravels MAPK1 regulated phosphoproteins in Leishmania donovani

Mitogen Activated Protein Kinase1 (MAPK1) of Leishmania donovani functions as key regulators of various cellular activities, which seem to be imperative for parasite survival, infectivity, drug resistance and post-translational modification of chaperones/co-chaperones. However, very less is known about LdMAPK1 target proteins. With recent advancements in proteomics, we aimed to identify phosphoproteins which were differentially expressed in LdMAPK1 overexpressing (Dd8++/++) and single replacement mutants (Dd8+/) as compared to wild type (Dd8+/+) parasites, utilizing LC-MS/MS approach. An in-depth label-free phospoproteomic analysis revealed that modulation of LdMAPK1 expression significantly modulates expression levels of miscellaneous phosphoproteins which may act as its targets/substrates. Out of 1974 quantified phosphoproteins in parasite, 140 were significantly differentially expressed in MAPK1 overexpressing and single replacement mutants. These differentially expressed phosphoproteins are majorly associated with metabolism, signal transduction, replication, transcription, translation, transporters and cytoskeleton/motor proteins, hence suggested that MAPK1 may act in concert to modulate global biological processes. The study further implicated possible role of LdMAPK1 in regulation and management of stress machinery in parasite through post translational modifications. Precisely, comparative phosphoproteomics study has elucidated significant role of LdMAPK1 in regulating various pathways contributing in parasite biology with relevance to future drug development. SIGNIFICANCE: MAPKinase1, the downstream kinase of MAPK signal transduction pathway, has drawn much attention as potential therapeutic drug target due to their indispensable role in survival and infectivity of Leishmania donovani. However, limited information is available about its downstream effector proteins/signaling networks. Utilizing label free LC-MS/MS analysis, phosphoproteome of LdMAPK1 over-expressing (Dd8++/++) and LdMAPK1 single replacement mutants (Dd8+/-) with wild type (Dd8+/+) parasites was compared and identified 140 LdMAPK1 modulated phosphoproteins, mainly involved in pathways like signal transduction, metabolism, transcriptional, translational, post-translational modification and regulation of heat shock proteins. Interestingly, LdMAPK1 interacts directly with only six phosphoproteins i.e. casein kinase, casein kinase II, HSP83/HSP90, LACK, protein kinase and serine/threonine protein kinase. Thus, the study elucidates significant role of LdMAPK1 in Leishmania biology which may drive drug-discovery efforts against visceral leishmaniasis.

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.

Leishmania: Responding to environmental signals and challenges without regulated transcription

Here we describe the non-canonical control of gene expression in Leishmania, a single-cell parasite that is responsible for one of the major neglected tropical diseases. We discuss the lack of regulated RNA synthesis, the post-transcriptional gene regulation including RNA stability and regulated translation. We also show that genetic adaptations such as mosaic aneuploidy, gene copy number variations and DNA sequence polymorphisms are important means for overcoming drug challenge and environmental diversity. These mechanisms are discussed in the context of the unique flow of genetic information found in Leishmania and related protists.

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.

May the epimastigote form of Trypanosoma cruzi be infective?

For many years it has been considered that there are three basic developmental stages of Trypanosoma cruzi: Epimastigote (Epi), Amastigote (Ama) and Trypomastigote (Typo). Epi and Ama are able to divide while Trypo does not divide. Epi are not infective while Ama and Trypo are able to infect host cells. Here we review the available data for the epimastigote form. Taken together the data show that (a) there are intermediate forms between epimastigotes and trypomastigotes in axenic cultures as well as between amastigote and trypomastigote forms within the cells (both in vitro and in vivo), and (c) that the intermediate forms, here designated as "Transitional Epimastigote", most of the time considered as epimastigotes, are able to infect cells. The recognition of the existence of this stage is of practical importance for those work with T. cruzi. Many laboratories working only with T. cruzi in axenic cultures usually consider to work with nonpathogenic cultures. This attitude needs to be changed requiring special care by those working with this protozoan to avoid accidental infections in the laboratory. In view of these observation a new scheme for the life cycle of T. cruzi is proposed.

Comparative transcriptomic analysis of antimony resistant and susceptible Leishmania infantum lines

Background

One of the major challenges to leishmaniasis treatment is the emergence of parasites resistant to antimony. To study differentially expressed genes associated with drug resistance, we performed a comparative transcriptomic analysis between wild-type and potassium antimonyl tartrate (Sb^III)-resistant Leishmania infantum lines using high-throughput RNA sequencing.

Methods

All the cDNA libraries were constructed from promastigote forms of each line, sequenced and analyzed using STAR for mapping the reads against the reference genome (L. infantum JPCM5) and DESeq2 for differential expression statistical analyses. All the genes were functionally annotated using sequence similarity search.

Results

The analytical pipeline considering an adjusted p-value < 0.05 and fold change > 2.0 identified 933 transcripts differentially expressed (DE) between wild-type and Sb^III-resistant L. infantum lines. Out of 933 DE transcripts, 504 presented functional annotation and 429 were assigned as hypothetical proteins. A total of 837 transcripts were upregulated and 96 were downregulated in the Sb^III-resistant L. infantum line. Using this DE dataset, the proteins were further grouped in functional classes according to the gene ontology database. The functional enrichment analysis for biological processes showed that the upregulated transcripts in the Sb^III-resistant line are associated with protein phosphorylation, microtubule-based movement, ubiquitination, host–parasite interaction, cellular process and other categories. The downregulated transcripts in the Sb^III-resistant line are assigned in the GO categories: ribonucleoprotein complex, ribosome biogenesis, rRNA processing, nucleosome assembly and translation.

Conclusions

The transcriptomic profile of L. infantum showed a robust set of genes from different metabolic pathways associated with the antimony resistance phenotype in this parasite. Our results address the complex and multifactorial antimony resistance mechanisms in Leishmania, identifying several candidate genes that may be further evaluated as molecular targets for chemotherapy of leishmaniasis.

Identification of Novel Interspersed DNA Repetitive Elements in the Trypanosoma cruzi Genome Associated with the 3'UTRs of Surface Multigenic Families

Trypanosoma cruzi is the etiological agent of Chagas disease, which affects millions of people in Latin America. No transcriptional control of gene expression has been demonstrated in this organism, and 50% of its genome consists of repetitive elements and members of multigenic families. In this study, we applied a novel bioinformatics approach to predict new repetitive elements in the genome sequence of T. cruzi. A new repetitive sequence measuring 241 nt was identified and found to be interspersed along the genome sequence from strains of different DTUs. This new repeat was mostly on intergenic regions, and upstream and downstream regions of the 241 nt repeat were enriched in surface protein genes. RNAseq analysis revealed that the repeat was part of processed mRNAs and was predominantly found in the 3′ untranslated regions (UTRs) of genes of multigenic families encoding surface proteins. Moreover, we detected a correlation between the presence of the repeat in the 3′UTR of multigenic family genes and the level of differential expression of these genes when comparing epimastigote and trypomastigote transcriptomes. These data suggest that this sequence plays a role in the posttranscriptional regulation of the expression of multigenic families.

The Leishmania donovani SENP Protease Is Required for SUMO Processing but Not for Viability

The protozoan parasite Leishmania donovani is part of an early eukaryotic branch and depends on post-transcriptional mechanisms for gene expression regulation. This includes post-transcriptional protein modifications, such as protein phosphorylation. The presence of genes for protein SUMOylation, i.e., the covalent attachment of small ubiquitin-like modifier (SUMO) polypeptides, in the Leishmania genomes prompted us to investigate the importance of the sentrin-specific protease (SENP) and its putative client, SUMO, for the vitality and infectivity of Leishmania donovani. While SENP null mutants are viable with reduced vitality, viable SUMO null mutant lines could not be obtained. SUMO C-terminal processing is disrupted in SENP null mutants, preventing SUMO from covalent attachment to proteins and nuclear translocation. Infectivity in vitro is not affected by the loss of SENP-dependent SUMO processing. We conclude that SENP is required for SUMO processing, but that functions of unprocessed SUMO are critical for Leishmania viability.

Ubiquitin-Like Modifiers: Emerging Regulators of Protozoan Parasites

Post-translational protein regulation allows for fine-tuning of cellular functions and involves a wide range of modifications, including ubiquitin and ubiquitin-like modifiers (Ubls). The dynamic balance of Ubl conjugation and removal shapes the fates of target substrates, in turn modulating various cellular processes. The mechanistic aspects of Ubl pathways and their biological roles have been largely established in yeast, plants, and mammalian cells. However, these modifiers may be utilised differently in highly specialised and divergent organisms, such as parasitic protozoa. In this review, we explore how these parasites employ Ubls, in particular SUMO, NEDD8, ATG8, ATG12, URM1, and UFM1, to regulate their unconventional cellular physiology. We discuss emerging data that provide evidence of Ubl-mediated regulation of unique parasite-specific processes, as well as the distinctive features of Ubl pathways in parasitic protozoa. We also highlight the potential to leverage these essential regulators and their cognate enzymatic machinery for development of therapeutics to protect against the diseases caused by protozoan parasites.

Casein kinase 1.2 over expression restores stress resistance to Leishmania donovani HSP23 null mutants

Leishmania donovani is a trypanosomatidic parasite and causes the lethal kala-azar fever, a neglected tropical disease. The Trypanosomatida are devoid of transcriptional gene regulation and rely on gene copy number variations and translational control for their adaption to changing conditions. To survive at mammalian tissue temperatures, L. donovani relies on the small heat shock protein HSP23, the loss of which renders the parasites stress sensitive and impairs their proliferation. Here, we analysed a spontaneous escape mutant with wild type-like in vitro growth. Further selection of this escape strains resulted in a complete reversion of the phenotype. Whole genome sequencing revealed a correlation between stress tolerance and the massive amplification of a six-gene cluster on chromosome 35, with further analysis showing over expression of the casein kinase 1.2 gene as responsible. In vitro phosphorylation experiments established both HSP23 and the related P23 co-chaperone as substrates and modulators of casein kinase 1.2, providing evidence for another crucial link between chaperones and signal transduction protein kinases in this early branching eukaryote.

Chitin binding protein as a possible RNA binding protein in Leishmania parasites

Leishmaniasis includes a broad spectrum of pathological outcomes in humans caused by protozoan parasites from the genus Leishmania. In recent years, proteomic techniques have introduced novel proteins with critical functions in Leishmania parasites. Based on our report of a Chitin binding protein (CBP) in our previous immunoproteomic study, this article suggests that CBP might be an RNA binding protein (RBP) in Leishmania parasites. RBPs, as key regulatory factors, have a role in post-transcriptional gene regulation. The presence of RBPs in Leishmania parasites has not been considered so far; however, this study aims to open a new venue regarding RBPs in Leishmania parasites. Confirming CBP as an RBP in Leishmania parasites, exploring other RBPs and their functions might lead to interesting issues in leishmaniasis. In fact, due to the regulatory role of RBPs in different diseases including cancers and their further classification as therapeutic targets, the emerging evaluation of CBP and RBPs from Leishmania parasites may allow the discovery of novel and effective drugs against leishmaniasis.

The Experimental Proteome of Leishmania infantum Promastigote and Its Usefulness for Improving Gene Annotations

Leishmania infantum causes visceral leishmaniasis (kala-azar), the most severe form of leishmaniasis, which is lethal if untreated. A few years ago, the re-sequencing and de novo assembling of the L. infantum (JPCM5 strain) genome was accomplished, and now we aimed to describe and characterize the experimental proteome of this species. In this work, we performed a proteomic analysis from axenic cultured promastigotes and carried out a detailed comparison with other Leishmania experimental proteomes published to date. We identified 2352 proteins based on a search of mass spectrometry data against a database built from the six-frame translated genome sequence of L. infantum. We detected many proteins belonging to organelles such as glycosomes, mitochondria, or flagellum, as well as many metabolic enzymes and many putative RNA binding proteins and molecular chaperones. Moreover, we listed some proteins presenting post-translational modifications, such as phosphorylations, acetylations, and methylations. On the other hand, the identification of peptides mapping to genomic regions previously annotated as non-coding allowed for the correction of annotations, leading to the N-terminal extension of protein sequences and the uncovering of eight novel protein-coding genes. The alliance of proteomics, genomics, and transcriptomics has resulted in a powerful combination for improving the annotation of the L. infantum reference genome.

Improvements on the quantitative analysis of Trypanosoma cruzi histone post translational modifications: Study of changes in epigenetic marks through the parasite's metacyclogenesis and life cycle

Trypanosome histone N-terminal sequences are very divergent from the other eukaryotes, although they are still decorated by post-translational modifications (PTMs). Here, we used a highly robust workflow to analyze histone PTMs in the parasite Trypanosoma cruzi using mass spectrometry-based (MS-based) data-independent acquisition (DIA). We adapted the workflow for the analysis of the parasite's histone sequences by modifying the software EpiProfile 2.0, improving peptide and PTM quantification accuracy. This workflow could now be applied to the study of 141 T. cruzi modified histone peptides, which we used to investigate the dynamics of histone PTMs along the metacyclogenesis and the life cycle of T. cruzi. Global levels of histone acetylation and methylation fluctuates along metacyclogenesis, however most critical differences were observed between parasite life forms. More than 66 histone PTM changes were detected. Strikingly, the histone PTM pattern of metacyclic trypomastigotes is more similar to epimastigotes than to cellular trypomastigotes. Finally, we highlighted changes at the H4 N-terminus and at H3K76 discussing their impact on the trypanosome biology. Altogether, we have optimized a workflow easily applicable to the analysis of histone PTMs in T. cruzi and generated a dataset that may shed lights on the role of chromatin modifications in this parasite. SIGNIFICANCE: Trypanosomes are unicellular parasites that have divergent histone sequences, no chromosome condensation and a peculiar genome/gene regulation. Genes are transcribed from divergent polycistronic regions and post-transcriptional gene regulation play major role on the establishment of transcripts and protein levels. In this regard, the fact that their histones are decorated with multiple PTMs raises interesting questions about their role. Besides, this digenetic organism must adapt to different environments changing its metabolism accordingly. As metabolism and epigenetics are closely related, the study of histone PTMs in trypanosomes may enlighten this strikingly, and not yet fully understood, interplay. From a biomedical perspective, the comprehensive study of molecular mechanisms associated to the metacyclogenesis process is essential to create better strategies for controlling Chagas disease.

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.

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.

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.

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.

Transcriptomic changes across the life cycle of Trypanosoma cruzi II

Trypanosoma cruzi is a flagellated protozoan that causes Chagas disease; it presents a complex life cycle comprising four morphological stages: epimastigote (EP), metacyclic trypomastigote (MT), cell-derived trypomastigote (CDT) and amastigote (AM). Previous transcriptomic studies on three stages (EPs, CDTs and AMs) have demonstrated differences in gene expressions among them; however, to the best of our knowledge, no studies have reported on gene expressions in MTs. Therefore, the present study compared differentially expressed genes (DEGs), and signaling pathway reconstruction in EPs, MTs, AMs and CDTs. The results revealed differences in gene expressions in the stages evaluated; these differences were greater between MTs and AMs-PTs. The signaling pathway that presented the highest number of DEGs in all the stages was associated with ribosomes protein profiles, whereas the other related pathways activated were processes related to energy metabolism from glucose, amino acid metabolism, or RNA regulation. However, the role of autophagy in the entire life cycle of T. cruzi and the presence of processes such as meiosis and homologous recombination in MTs (where the expressions of SPO11 and Rad51 plays a role) are crucial. These findings represent an important step towards the full understanding of the molecular basis during the life cycle of T. cruzi.

High performance of an enzyme linked immunosorbent assay for American tegumentary leishmaniasis diagnosis with Leishmania (Viannia) braziliensis amastigotes membrane crude antigens

The diagnosis of American tegumentary leishmaniasis (ATL) still requires the design of more effective tools. Leishmania (Viannia) braziliensis is the causal agent of the 90% of Argentinean ATL cases. Considering the current knowledge, an ELISA based crude antigen (CA) for the diagnosis was designed. Ninety-nine subjects diagnosed as ATL, 27 as no-ATL, and 84 donors from non-ATL-endemic areas were included in this study. The current ATL diagnosis was based four techniques, dermal smear microscopic examination (parasitological test), PCR, Leishmanin skin test, and clinical records. We obtained CA extracts from promastigotes and amastigotes from macrophage cultures of different zymodemes of endemic Leishmania species circulating in the study area. Crude antigens from the ‘local’ main zymodeme of L. (V.) braziliensis showed the highest reactivity against anti-Leishmania antibodies compared to the other included species. The CA of amastigotes of this zymodeme was 3.4 fold more reactive than promastigotes one. Moreover, amastigote-membrane CA (MCA) were 3.6 fold more reactive than the soluble antigens. The MCA-ELISA reached a sensitivity and specificity of 98% (CI = 94.7%-100%) and 63.6% (53.9–73.1), respectively. When anti-Trypanosoma cruzi reactive sera were excluded, the specificity reached 98.4% (94.4–100), while the sensitivity was similar, with a positive predictive value (PV) of 98.6% (94.6–100) and negative PV of 96.3% (91.6–100). The performance of the MCA-ELISA results strongly contribute to the final diagnostic decision, since a non-reactive serological result almost discards the suspected ATL, because of its high negative PV. The developed MCA-ELISA showed a high diagnostic performance, which makes it a good candidate for ATL diagnosis, for seroprevalence studies, or for monitoring treatments efficacy.

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.

Identification of a novel base J binding protein complex involved in RNA polymerase II transcription termination in trypanosomes

Base J, β-D-glucosyl-hydroxymethyluracil, is a modification of thymine DNA base involved in RNA Polymerase (Pol) II transcription termination in kinetoplastid protozoa. Little is understood regarding how specific thymine residues are targeted for J-modification or the mechanism of J regulated transcription termination. To identify proteins involved in J-synthesis, we expressed a tagged version of the J-glucosyltransferase (JGT) in Leishmania tarentolae, and identified four co-purified proteins by mass spectrometry: protein phosphatase (PP1), a homolog of Wdr82, a potential PP1 regulatory protein (PNUTS) and a protein containing a J-DNA binding domain (named JBP3). Gel shift studies indicate JBP3 is a J-DNA binding protein. Reciprocal tagging, co-IP and sucrose gradient analyses indicate PP1, JGT, JBP3, Wdr82 and PNUTS form a multimeric complex in kinetoplastids, similar to the mammalian PTW/PP1 complex involved in transcription termination via PP1 mediated dephosphorylation of Pol II. Using RNAi and analysis of Pol II termination by RNA-seq and RT-PCR, we demonstrate that ablation of PNUTS, JBP3 and Wdr82 lead to defects in Pol II termination at the 3'-end of polycistronic gene arrays in Trypanosoma brucei. Mutants also contain increased antisense RNA levels upstream of transcription start sites, suggesting an additional role of the complex in regulating termination of bi-directional transcription. In addition, PNUTS loss causes derepression of silent Variant Surface Glycoprotein genes involved in host immune evasion. Our results suggest a novel mechanistic link between base J and Pol II polycistronic transcription termination in kinetoplastids.

Glycolytic profile shift and antioxidant triggering in symbiont-free and H2O2-resistant Strigomonas culicis

During their life cycle, trypanosomatids are exposed to stress conditions and adapt their energy and antioxidant metabolism to colonize their hosts. Strigomonas culicis is a monoxenous protist found in invertebrates with an endosymbiotic bacterium that completes essential biosynthetic pathways for the trypanosomatid. Our research group previously generated a wild-type H₂O₂-resistant (WTR) strain that showed improved mitochondrial metabolism and antioxidant defenses, which led to higher rates of Aedes aegypti infection. Here, we assess the biological contribution of the S. culicis endosymbiont and reactive oxygen species (ROS) resistance to oxidative and energy metabolism processes. Using high-throughput proteomics, several proteins involved in glycolysis and gluconeogenesis, the pentose phosphate pathway and glutathione metabolism were identified. The results suggest that ROS resistance decreases glucose consumption and indicate that the metabolic products from gluconeogenesis are key to supplying the protist with high-energy and reducing intermediates. Our hypothesis was confirmed by biochemical assays showing opposite profiles for glucose uptake and hexokinase and pyruvate kinase activity levels in the WTR and aposymbiotic strains, while the enzyme glucose-6P 1-dehydrogenase was more active in both strains. Regarding the antioxidant system, ascorbate peroxidase has an important role in H₂O₂ resistance and may be responsible for the high infection rates previously described for A. aegypti. In conclusion, our data indicate that the energy-related and antioxidant metabolic processes of S. culicis are modulated in response to oxidative stress conditions, providing new perspectives on the biology of the trypanosomatid-insect interaction as well as on the possible impact of resistant parasites in accidental human infection.

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.

Quantitative Proteomic Map of the Trypanosomatid Strigomonas culicis: The Biological Contribution of its Endosymbiotic Bacterium

Strigomonas culicis is a kinetoplastid parasite of insects that maintains a mutualistic association with an intracellular symbiotic bacterium, which is highly integrated into the protist metabolism: it furnishes essential compounds and divides in synchrony with the eukaryotic nucleus. The protist, conversely, can be cured of the endosymbiont, producing an aposymbiotic cell line, which presents a diminished ability to colonize the insect host. This obligatory association can represent an intermediate step of the evolution towards the formation of an organelle, therefore representing an interesting model to understand the symbiogenesis theory. Here, we used shotgun proteomics to compare the S. culicis endosymbiont-containing and aposymbiotic strains, revealing a total of 11,305 peptides, and up to 2,213 proteins (2,029 and 1,452 for wild type and aposymbiotic, respectively). Gene ontology associated to comparative analysis between both strains revealed that the biological processes most affected by the elimination of the symbiont were the amino acid synthesis, as well as protein synthesis and folding. This large-scale comparison of the protein expression in S. culicis marks a step forward in the comprehension of the role of endosymbiotic bacteria in monoxenous trypanosomatid biology, particularly because trypanosomatids expression is mostly post-transcriptionally regulated.