Analysis of the Leishmania donovani transcriptome reveals an ordered progression of transient and permanent changes in gene expression during differentiation

Leishmania donovani mitogen-activated protein kinases as a host-parasite interaction interface

Leishmaniasis, a major globally re-emerging neglected tropical disease, has a restricted repertoire of chemotherapeutic options due to a narrow therapeutic index, drug resistance, or patient non-compliance due to toxicity. The disease is caused by the parasite Leishmania that resides in two different forms in two different environments: as sessile intracellular amastigotes within mammalian macrophages and as motile promastigotes in sandfly gut. As mitogen-activated protein kinases (MAPKs) play important roles in cellular differentiation and survival, we studied the expression of Leishmania donovani MAPKs (LdMAPKs). The homology studies by multiple sequence alignment show that excepting LdMAPK1 and LdMAPK2, all thirteen other LdMAPKs share homology with human ERK and p38 isoforms. Expression of LdMAPK4 and LdMAPK5 is less in avirulent promastigotes and amastigotes. Compared to miltefosine-sensitive L. donovani parasites, miltefosine-resistant parasites have higher LdMAPK1, LdMAPK3-5, LdMAPK7-11, LdMAPK13, and LdMAPK14 expression. IL-4-treatment of macrophages down-regulated LdMAPK11, in virulent amastigotes whereas up-regulated LdMAPK5, but down-regulated LdMAPK6, LdMAPK12-15, expression in avirulent amastigotes. IL-4 up-regulated LdMAPK1 expression in both virulent and avirulent amastigotes. IFN-γ-treatment down-regulated LdMAPK6, LdMAPK13, and LdMAPK15 in avirulent amastigotes but up-regulated in virulent amastigotes. This complex profile of LdMAPKs expression among virulent and avirulent parasites, drug-resistant parasites, and in amastigotes within IL-4 or IFN-γ-treated macrophages suggests that LdMAPKs are differentially controlled at the host-parasite interface regulating parasite survival and differentiation, and in the course of IL-4 or IFN-γ dominated immune response.

Whole cell reconstructions of Leishmania mexicana through the cell cycle

The unicellular parasite Leishmania has a precisely defined cell architecture that is inherited by each subsequent generation, requiring a highly coordinated pattern of duplication and segregation of organelles and cytoskeletal structures. A framework of nuclear division and morphological changes is known from light microscopy, yet this has limited resolution and the intrinsic organisation of organelles within the cell body and their manner of duplication and inheritance is unknown. Using volume electron microscopy approaches, we have produced three-dimensional reconstructions of different promastigote cell cycle stages to give a spatial and quantitative overview of organelle positioning, division and inheritance. The first morphological indications seen in our dataset that a new cell cycle had begun were the assembly of a new flagellum, the duplication of the contractile vacuole and the increase in volume of the nucleus and kinetoplast. We showed that the progression of the cytokinesis furrow created a specific pattern of membrane indentations, while our analysis of sub-pellicular microtubule organisation indicated that there is likely a preferred site of new microtubule insertion. The daughter cells retained these indentations in their cell body for a period post-abscission. By comparing cultured and sand fly derived promastigotes, we found an increase in the number and overall volume of lipid droplets in the promastigotes from the sand fly, reflecting a change in their metabolism to ensure transmissibility to the mammalian host. Our insights into the cell cycle mechanics of Leishmania will support future molecular cell biology analyses of these parasites.

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.

Evaluation of expression variations in virulence-related genes of Leishmania major after several culture passages compared with Phlebotomus papatasi isolated promastigotes

Cutaneous leishmaniasis (CL) is a prevalent infectious disease with considerable morbidity annually. Here, we aimed to investigate the likely variations in gene expression of glycoprotein63 (gp63), heat shock protein 70 (HSP70), histone, arginase, cysteine protease B (CPB), Leishmania homologue of receptors for activated C kinase (LACK), small hydrophilic endoplasmic reticulum-associated protein (SHERP) in metacyclic promastigotes of L. major isolated from Phlebotomus papatasi sand flies and promastigotes excessively cultured in culture medium. The parasites were collected from suspected CL cases in Pasteur Institute of Iran, cultured and inoculated into the female BALB/c mice (2×106 promastigotes). Sand flies were trapped in Qom province, fed with the blood of euthanized infected mice and subsequently dissected in order to isolate the midgut including stomodeal valve. The metacyclic promastigotes were isolated from Ph. papatasi (Pro-Ppap) using peanut agglutinin test (PNA), then continuously cultured in RPMI-1640 medium enriched with fetal bovine serum, penicillin (100 U/ml) and streptomycin (100 mg/ml) to reach stationary phase (Pro-Stat). The gene expression was evaluated in both parasitic stages (Pro-Ppap and Pro-Stat) using qRT-PCR. Out results showed a significant increased gene expression at Pro-Ppap stage for gp63 (P = 0.002), SHERP (P = 0.001) and histone (P = 0.026) genes, in comparison with Pro-Stat stage. Noticeably, significant changes were, also, demonstrated in 10th to 15th passages [gp63 (P = 0.041), arginase (P = 0.016), LACK (P = 0.025)] and in 5th to 20th passage (SHERP) (P = 0.029). In conclusion, the findings of the present study seem to be essential in designing Leishmania studies, in particular regarding host-parasite interaction, immunization and infectivity studies.

An insight into differential protein abundance throughout Leishmania donovani promastigote growth and differentiation

Leishmania donovani causes anthroponotic visceral leishmaniasis, responsible for about 50,000 annual deaths worldwide. Current therapies have considerable side effects. Drug resistance has been reported and no vaccine is available nowadays. The development of undifferentiated promastigotes in the sand fly vector's gut leads to the promastigote form that is highly infective to the mammalian host. Fully differentiated promastigotes play a crucial role in the initial stages of mammalian host infection before internalization in the host phagocytic cell. Therefore, the study of protein levels in the promastigote stage is relevant for disease control, and proteomics analysis is an ideal source of vaccine candidate discovery. This study aims to get insight into the protein levels during the differentiation process of promastigotes by 2DE-MALDI-TOF/TOF. This partial proteome analysis has led to the identification of 75 proteins increased in at least one of the L. donovani promastigote differentiation and growth phases. This study has revealed the differential abundance of said proteins during growth and differentiation. According to previous studies, some are directly involved in parasite survival or are immunostimulatory. The parasite survival-related proteins are ascorbate peroxidase; cystathionine β synthase; an elongation factor 1β paralog; elongation factor 2; endoribonuclease L-PSP; an iron superoxide dismutase paralog; GDP-mannose pyrophosphorylase; several heat shock proteins-HSP70, HSP83-17, mHSP70-rel, HSP110; methylthioadenosine phosphorylase; two thiol-dependent reductase 1 paralogs; transitional endoplasmic reticulum ATPase; and the AhpC thioredoxin paralog. The confirmed immunostimulatory proteins are the heat shock proteins, enolase, and protein kinase C receptor analog. The potential immunostimulatory molecules according to findings in patogenic bacteria are fructose-1,6-diphophate aldolase, dihydrolipoamide acetyltransferase, isocitrate dehydrogenase, pyruvate dehydrogenase E1α and E1β subunits, and triosephosphate isomerase. These proteins may become disease control candidates through future intra-vector control methods or vaccines.

Identification of protein profile in metacyclic and amastigote-like stages of Leishmania tropica: a proteomic approach

Leishmaniasis is a tropical disease that leads to various clinical phenotypes. This study aimed to investigate protein expression changes in metacyclic and amastigote-like stages of L. tropica isolated from Iranian cutaneous leishmaniasis patients. Isolated samples were cultured and species type identified using PCR–RFLP technique. The promastigotes were grown in RPMI1640 media and differentiated to metacyclic and amastigote-like forms, followed by the extracted proteins of both successive stages carried out for proteomics and bioinformatics analysis. Using SWATH-MS quantitative proteomics technique, a total 176 and 155 distinct proteins were identified in metacyclic and axenic amastigote stages, respectively. Of these, 65 proteins were altered significantly (p-value < 0.05 and fold change ≥ 2) between studied stages. Several gene ontology (GO) categories were enriched for biological process during conversion of metacyclic promastigotes into amastigote-like, which “metabolic process” (GO: 0044281, P-Value: 6.52e-5), and “translation” (GO: 0006412, p-value: 5.01e–14) were disclosed as the top category in up and down-regulated proteins, respectively. Also, the KEGG pathway analysis indicated “metabolic pathways” and “ribosome” term as the most important pathways in up and down-regulated proteins, respectively. According to protein interaction network analysis, enolase (ENOL) has been detected as main hub proteins during differentiation, followed by Putative NADH-dependent fumarate reductase (LmjF.35.1180) and 40S ribosomal protein S2 (LmjF.32.0450). Overall, protein changes possibly play important roles in L. tropica biology. Anabolic pathways were down-regulated, whereas catabolic pathways were up-regulated during L. tropica differentiation. These protein expression changes could provide parasite survival in host macrophages, and could use as novel potential drug and vaccine targets for leishmaniasis.

Graphical Abstract

Quantitative single-cell analysis of Leishmania major amastigote differentiation demonstrates variably extended expression of the lipophosphoglycan (LPG) virulence factor in different host cell types

Immediately following their deposition into the mammalian host by an infected sand fly vector, Leishmania parasites encounter and are engulfed by a variety of cell types. From there, parasites may transit to other cell types, primarily macrophages or dendritic cells, where they replicate and induce pathology. During this time, Leishmania cells undergo a dramatic transformation from the motile non-replicating metacyclic stage to the non-motile replicative amastigote stage, a differentiative process that can be termed amastigogenesis. To follow this at the single cell level, we identified a suite of experimental 'landmarks' delineating different stages of amastigogenesis qualitatively or quantitatively, including new uses of amastigote-specific markers that showed interesting cellular localizations at the anterior or posterior ends. We compared amastigogenesis in synchronous infections of peritoneal and bone-marrow derived macrophages (PEM, BMM) or dendritic cells (BMDC). Overall, the marker suite expression showed an orderly transition post-infection with similar kinetics between host cell types, with the emergence of several amastigote traits within 12 hours, followed by parasite replication after 24 hours, with parasites in BMM or BMDC initiating DNA replication more slowly. Lipophosphoglycan (LPG) is a Leishmania virulence factor that facilitates metacyclic establishment in host cells but declines in amastigotes. Whereas LPG expression was lost by parasites within PEM by 48 hours, >40% of the parasites infecting BMM or BMDC retained metacyclic-level LPG expression at 72 hr. Thus L. major may prolong LPG expression in different intracellular environments, thereby extending its efficacy in promoting infectivity in situ and during cell-to-cell transfer of parasites expressing this key virulence factor.

Evaluating complete surface-associated and secretory proteome of Leishmania donovani for discovering novel vaccines and diagnostic targets

The protozoa Leishmania donovani causes visceral leishmaniasis (kala-azar), the third most common vector-borne disease. The visceral organs, particularly the spleen, liver, and bone marrow, are affected by the disease. The lack of effective treatment regimens makes curing and eradicating the disease difficult. The availability of complete L. donovani genome/proteome data allows for the development of specific and efficient vaccine candidates using the reverse vaccinology method, while utilizing the unique sequential and structural features of potential antigenic proteins to induce protective T cell and B cell responses. Such shortlisted candidates may then be tested quickly for their efficacy in the laboratory and later in clinical settings. These antigens will also be useful for designing antigen-based next-generation sero-diagnostic assays. L. donovani's cell surface-associated proteins and secretory proteins are among the first interacting entities to be exposed to the host immune machinery. As a result, potential antigenic epitope peptides derived from these proteins could serve as competent vaccine components. We used a stepwise filtering-based in silico approach to identify the entire surface-associated and secretory proteome of L. donovani, which may provide rationally selected most exposed antigenic proteins. Our study identified 12 glycosylphosphatidylinositol-anchored proteins, 45 transmembrane helix-containing proteins, and 73 secretory proteins as potent antigens unique to L. donovani. In addition, we used immunoinformatics to identify B and T cell epitopes in them. Out of the shortlisted surface-associated and secretory proteome, 66 protein targets were found to have the most potential overlapping B cell and T cell epitopes (linear and conformational; MHC class I and MHC class II).

Monomethylsulochrin isolated from biomass extract of Aspergillus sp. against Leishmania amazonensis: In vitro biological evaluation and molecular docking

Leishmaniasis represents a serious world health problem, with 1 billion people being exposed to infection and a broad spectrum of clinical manifestations with a potentially fatal outcome. Based on the limitations observed in the treatment of leishmaniasis, such as high cost, significant adverse effects, and the potential for drug resistance, the aim of the present study was to evaluate the leishmanicidal activity of the compounds pseurotin A and monomethylsulochrin isolated from the biomass extract of Aspergillus sp. The chromatographic profiles of the extract were determined by high-performance liquid chromatography coupled with a diode-array UV-Vis detector (HPLC-DAD-UV), and the molecular identification of the pseurotin A and monomethylsulochrin were carried out by electrospray ionization mass spectrometry in tandem (LC-ESI-MS-MS) and nuclear magnetic resonance (NMR). Antileishmanial activity was assayed against promastigote and intracellular amastigote of Leishmania amazonensis. As a control, cytotoxicity assays were performed in non-infected BALB/c peritoneal macrophages. Ultrastructural alterations in parasites were evaluated by transmission electron microscopy. Changes in mitochondrial membrane potential were determined by flow cytometry. Only monomethylsulochrin inhibited the promastigote growth (IC₅₀ 18.04 ± 1.11 µM), with cytotoxicity to peritoneal macrophages (CC₅₀ 5.09 91.63 ± 1.28 µM). Activity against intracellular amastigote forms (IC₅₀ 5.09 ± 1.06 µM) revealed an increase in antileishmanial activity when compared with promastigotes. In addition to a statistically significant reduction in the evaluated infection parameters, monomethylsulochrin altered the ultrastructure of the promastigote forms with atypical vacuoles, electron-dense corpuscles in the cytoplasm, changes at the mitochondria outer membrane and abnormal disposition around the kinetoplast. It was showed that monomethylsulochrin leads to a decrease in the mitochondrial membrane potential (25.9%, p = 0.0286). Molecular modeling studies revealed that monomethylsulochrin can act as inhibitor of sterol 14-alpha-demethylase (CYP51), a therapeutic target for human trypanosomiasis and leishmaniasis. Assessed for its drug likeness, monomethylsulochrin follows the Lipinski Rule of five and Ghose, Veber, Egan, and Muegge criteria. Furthermore, monomethylsulochrin can be used as a reference in the development of novel and therapeutically useful antileishmanial agents.

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.

In vivo transcriptional analysis of mice infected with Leishmania major unveils cellular heterogeneity and altered transcriptomic profiling at single-cell resolution

Leishmania parasites cause cutaneous leishmaniasis (CL), a disease characterized by disfiguring, ulcerative skin lesions. Both parasite and host gene expression following infection with various Leishmania species has been investigated in vitro, but global transcriptional analysis following L. major infection in vivo is lacking. Thus, we conducted a comprehensive transcriptomic profiling study combining bulk RNA sequencing (RNA-Seq) and single-cell RNA sequencing (scRNA-Seq) to identify global changes in gene expression in vivo following L. major infection. Bulk RNA-Seq analysis revealed that host immune response pathways like the antigen processing and presentation pathway were significantly enriched amongst differentially expressed genes (DEGs) upon infection, while ribosomal pathways were significantly downregulated in infected mice compared to naive controls. scRNA-Seq analyses revealed cellular heterogeneity including distinct resident and recruited cell types in the skin following murine L. major infection. Within the individual immune cell types, several DEGs indicative of many interferon induced GTPases and antigen presentation molecules were significantly enhanced in the infected ears including macrophages, resident macrophages, and inflammatory monocytes. Ingenuity Pathway Analysis of scRNA-Seq data indicated the antigen presentation pathway was increased with infection, while EIF2 signaling is the top downregulated pathway followed by eIF4/p70S6k and mTOR signaling in multiple cell types including macrophages, blood and lymphatic endothelial cells. Altogether, this transcriptomic profile highlights known recruitment of myeloid cells to lesions and recognizes a potential role for EIF2 signaling in murine L. major infection in vivo.

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.

Tissue Specific Dual RNA-Seq Defines Host-Parasite Interplay in Murine Visceral Leishmaniasis Caused by Leishmania donovani and Leishmania infantum

Visceral leishmaniasis is associated with hepato-splenomegaly and altered immune and hematological parameters in both preclinical animal models and humans. We studied mouse experimental visceral leishmaniasis caused by Leishmania infantum and Leishmania donovani in BALB/c mice using dual RNA-seq to investigate the transcriptional response of host and parasite in liver and spleen. We identified only 4 species-specific parasite expressed genes (SSPEGs; log2FC >1, FDR <0.05) in the infected spleen, and none in the infected liver. For the host transcriptome, we found 789 differentially expressed genes (DEGs; log2FC >1, FDR <0.05) in the spleen that were common to both infections, with IFNγ signaling and complement and coagulation cascade pathways highly enriched, and an additional 286 and 186 DEGs that were selective to L. donovani and L. infantum infection, respectively. Among those, there were network interactions between genes of amino acid metabolism and PPAR signaling in L. donovani infection and increased IL1β and positive regulation of fatty acid transport in L. infantum infection, although no pathway enrichment was observed. In the liver, there were 1,939 DEGs in mice infected with either L. infantum or L. donovani in comparison to uninfected mice, and the most enriched pathways were IFNγ signaling, neutrophil mediated immunity, complement and coagulation, cytokine-chemokine responses, and hemostasis. Additionally, 221 DEGs were selective in L. donovani and 429 DEGs in L. infantum infections. These data show that the host response for these two visceral leishmaniasis infection models is broadly similar, and ∼10% of host DEGs vary in infections with either parasite species. IMPORTANCE Visceral leishmaniasis (VL) is caused by two species of Leishmania parasites, L. donovani in the Old World and L. infantum in the New World and countries bordering the Mediterranean. Although cardinal features such as hepato-splenomegaly and alterations in blood and immune function are evident, clinical presentation may vary by geography, with for example severe bleeding often associated with VL in Brazil. Although animal models of both L. donovani and L. infantum have been widely used to study disease pathogenesis, a direct side-by-side comparison of how these parasites species impact the infected host and/or how they might respond to the stresses of mammalian infection has not been previously reported. Identifying common and distinct pathways to pathogenesis will be important to ensure that new therapeutic or prophylactic approaches will be applicable across all forms of VL.

Mutation Characteristics and Phylogenetic Analysis of Five Leishmania Clinical Isolates

Simple Summary

Leishmaniasis, a neglected tropical disease, is caused by infection with the Leishmania species, threatening millions of people in approximately 100 endemic countries. The emergence of antimony-resistant Leishmania strains have brought difficulties to the treatment and elimination of leishmaniasis. This study performed genome-wide resequencing and phylogenetic analysis of five isolates from the Leishmania donovani complex, focusing on finding mutations related to antimony resistance and virulence of the newly isolated Leishmania strain L_HCZ in 2016. By combining whole-genome sequencing and whole-genome phylogenetic analysis, Leishmania isolates L_801, L_9044 and L_Liu were identified as Leishmania donovani, and L_HCZ as Leishmania infantum. By discovering genome-wide single-nucleotide polymorphisms and structural variations, we identified mutations of drug resistance-related genes in the antimony-resistant Leishmania isolate L_HCZ. The new Leishmania isolate L_HCZ has strong virulence and strong drug resistance, which should be taken seriously by the relevant health departments and scientific researchers.

Abstract

Leishmaniasis is a neglected tropical disease threatening millions of people worldwide. The emergence of antimony-resistant Leishmania strains have brought difficulties to the treatment and elimination of leishmaniasis. This study performed genome sequencing, phylogenetic analysis and mutation analysis of five Leishmania clinical isolates, especially the Leishmania strain L_HCZ isolated in 2016, which shows strong virulence and antimony resistance. By phylogenetic analysis, four isolates (L_DD8, L_801, L_Liu and L_9044) were identified as Leishmania donovani, the isolate L_HCZ was identified as Leishmania infantum and the isolate L_DD8 as a standard strain of L.donovani. Genome-wide mutation analysis was applied to identify mutations related to the drug resistance and virulence of the newly isolated L_HCZ. Compared with the other four Leishmania isolates, L_HCZ had the most mutations in genes associated with antimony resistance, including the ABC transporter, ascorbate-dependent peroxidase, gamma–glutamylcysteine synthetase, glucose-6-phosphate 1-dehydrogenase, ATP-binding cassette protein subfamily A and multi-drug resistance protein-like genes. Among the genes associated with virulence, L_HCZ had the most mutations in cysteine peptidase A, cysteine peptidase B, cysteine peptidase C, heat-shock protein 70, gp63, acid phosphatase, kinesin k39, kinesin, phosphoglycan beta 1, amastin-like surface protein and amastin-like proteins. The mutations in L_HCZ might possibly contribute to its antimony resistance and strong virulence in clinical patients. Whole-genome resequencing has exhibited broad application prospects and may be put into clinical use in the future for parasite identifying and epidemiological investigations.

Axenic amastigotes of Leishmania species as a suitable model for in vitro studies

Leishmania spp. are etiological agents of infection diseases, which in some cases can be fatal. The main forms of their biological cycle, promastigotes and amastigotes, can be maintained in vitro. While promastigotes are easier to maintain, amastigotes are more complex and can be obtained through different ways, including infection assays of tissues or in vitro cells, and differentiation from promastigotes to axenic amastigotes. Several protocols have been proposed for in vitro differentiation for at least 12 Leishmania spp. of both subgenera, Leishmania and Viannia. In this review we propose a critical summary of axenic amastigotes induction, as well as the impact of these strategies on metabolic pathways and regulatory networks analyzed by omics approaches. The parameters used by different research groups show considerable variations in temperature, pH and induction stages, as highlighted here for Leishmania (Viannia) braziliensis. Therefore, a consensus on strategies for inducing amastigogenesis is necessary to improve accuracy and even define stage-specific biomarkers. In fact, the axenic amastigote model has contributed to elucidate several aspects of the parasite cycle, however, since it does not reproduce the intracellular environment, its use requires several precautions. In addition, we present a discussion about using axenic amastigotes for drug screening, suggesting the need of a more sensitive methodology to verify cell viability in these tests. Collectively, this review explores the advantages and limitations found in studies with axenic amastigotes, done for more than 30 years, and discuss the gaps that impair their use as a suitable model for in vitro studies.

Involvement of Leishmania Phosphatases in Parasite Biology and Pathogeny

In the Leishmania lifecycle, the motile promastigote form is transmitted from the sand fly vector to a mammalian host during a blood meal. Inside vertebrate host macrophages, the parasites can differentiate into the amastigote form and multiply, causing leishmaniasis, one of the most significant neglected tropical diseases. Leishmania parasites face different conditions throughout their development inside sand flies. Once in the mammalian host, the parasites have to overcome the microbicide repertoire of the cells of the immune system to successfully establish the infection. In this context, the expression of protein phosphatases is of particular interest. Several members of the serine/threonine-specific protein phosphatase (STP), protein tyrosine phosphatase (PTP), and histidine acid phosphatase (HAcP) families have been described in different Leishmania species. Although their physiological roles have not been fully elucidated, many studies suggest they have an involvement with parasite biology and pathogeny. Phosphatases play a role in adaptation to nutrient starvation during parasite passage through the sand fly midgut. They are also important to parasite virulence, mainly due to the modulation of host cytokine production and impairment of the microbiocidal potential of macrophages. Furthermore, recent whole-genome expression analyses have shown that different phosphatases are upregulated in metacyclic promastigotes, the infective form of the mammalian host. Leishmania phosphatases are also upregulated in drug-resistant strains, probably due to the increase in drug efflux related to the activation of ABC transporters. Throughout this review, we will describe the physiological roles that have been attributed to Leishmania endogenous phosphatases, including their involvement in the adaptation, survival, and proliferation of the parasites inside their hosts.

Leishmania major RUVBL1 has a hexameric conformation in solution and, in the presence of RUVBL2, forms a heterodimer with ATPase activity

ATPases belonging to the AAA+ superfamily are associated with diverse cellular activities and are mainly characterized by a nucleotide-binding domain (NBD) containing the Walker A and Walker B motifs. AAA+ proteins have a range of functions, from DNA replication to protein degradation. Rvbs, also known as RUVBLs, are AAA+ ATPases with one NBD domain and were described from human to yeast as participants of the R2TP (Rvb1-Rvb2-Tah1-Pih1) complex. Although essential for the assembly of multiprotein complexes-containing DNA and RNA, the protozoa Rvb orthologs are less studied. For the first time, this work describes the Rvbs from Leishmania major, one of the causative agents of Tegumentar leishmaniasis in human. Recombinant LmRUVBL1 and LmRUVBL2 his-tagged proteins were successfully purified and investigated using biophysical tools. LmRUVBL1 was able to form a well-folded elongated hexamer in solution, while LmRUVBL2 formed a large aggregate. However, the co-expression of LmRUVBL1 and LmRUVBL2 assembled the proteins into an elongated heterodimer in solution. Thermo-stability and fluorescence experiments indicated that the LmRUVBL1/2 heterodimer had ATPase activity in vitro. This is an interesting result because hexameric LmRUVBL1 alone had low ATPase activity. Additionally, using independent SL-RNAseq libraries, it was possible to show that both proteins are expressed in all L. major life stages. Specific antibodies obtained against LmRUVBLs identified the proteins in promastigotes and metacyclics cell extracts. Together, the results here presented are the first step towards the characterization of Leishmania Rvbs, and may contribute to the development of possible strategies to intervene against leishmaniasis, a neglected tropical disease of great medical importance.

CRISPRing protozoan parasites to better understand the biology of diseases

Precise gene editing techniques are paramount to gain deeper insights into the biological processes such as host-parasite interactions, drug resistance mechanisms, and gene-function relationships. Discovery of CRISPR-Cas9 system has spearheaded mechanistic understanding of protozoan parasite biology as evident from the number of reports in the last decade. Here, we have described the use of CRISPR-Cas9 in understanding the biology of medically important protozoan parasites such as Plasmodium, Leishmania, Trypanosoma, Babesia and Trichomonas. In spite of intrinsic difficulties in genome editing in these protozoan parasites, CRISPR-Cas9 has acted as a catalyst for faster generation of desired transgenic parasites. Modifications in the CRISPR-Cas9 system for improving the efficiency have been useful in better understanding the molecular mechanisms associated with repair of double strand breaks in the parasites. Moreover, improvement in reagents used for CRISPR mediated gene editing have been instrumental in addressing the issue of non-specificity and toxicity for therapeutic use. These application-based modifications may help in further increasing the efficiency of gene editing in protozoan parasites.

Potential therapeutic targets shared between leishmaniasis and cancer

The association of leishmaniasis and malignancies in human and animal models has been highlighted in recent years. The misdiagnosis of coexistence of leishmaniasis and cancer and the use of common drugs in the treatment of such diseases prompt us to further survey the molecular biology of Leishmania parasites and cancer cells. The information regarding common expressed proteins, as possible therapeutic targets, in Leishmania parasites and cancer cells is scarce. Therefore, the current study reviews proteins, and investigates the regulation and functions of several key proteins in Leishmania parasites and cancer cells. The up- and down-regulations of such proteins were mostly related to survival, development, pathogenicity, metabolic pathways and vital signalling in Leishmania parasites and cancer cells. The presence of common expressed proteins in Leishmania parasites and cancer cells reveals valuable information regarding the possible shared mechanisms of pathogenicity and opportunities for therapeutic targeting in leishmaniasis and cancers in the future.

Antileishmanial Chemotherapy through Clemastine Fumarate Mediated Inhibition of the Leishmania Inositol Phosphorylceramide Synthase

Current chemotherapeutics for leishmaniasis have multiple deficiencies, and there is a need for new safe, efficacious, and affordable medicines. This study describes a successful drug repurposing approach that identifies the over-the-counter antihistamine, clemastine fumarate, as a potential antileishmanial drug candidate. The screening for inhibitors of the sphingolipid synthase (inositol phosphorylceramide synthase, IPCS) afforded, following secondary screening against Leishmania major (Lmj) promastigotes, 16 active compounds. Further refinement through the dose response against LmjIPCS and intramacrophage L. major amastigotes identified clemastine fumarate with good activity and selectivity with respect to the host macrophage. On target engagement was supported by diminished sensitivity in a sphingolipid-deficient L. major mutant (ΔLmjLCB2) and altered phospholipid and sphingolipid profiles upon treatment with clemastine fumarate. The drug also induced an enhanced host cell response to infection indicative of polypharmacology. The activity was sustained across a panel of Old and New World Leishmania species, displaying an in vivo activity equivalent to the currently used drug, glucantime, in a mouse model of L. amazonensis infection. Overall, these data validate IPCS as an antileishmanial drug target and indicate that clemastine fumarate is a candidate for repurposing for the treatment of leishmaniasis.

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.

Stage-Specific Differential Gene Expression of Glutathione Peroxidase in Leishmania Major and Leishmania Tropica

BACKGROUND

Leishmania (L) major and L. tropica are the etiological agents of cutaneous leishmaniosis. Leishmania species cause a board spectrum of phenotypes. A small number of genes are differentially expressed between them that have likely an important role in the disease phenotype. Procyclic and metacyclic are two morphological promastigote forms of Leishmania that express different genes. The glutathione peroxidase is an important antioxidant enzyme that essential in parasite protection against oxidative stress and parasite survival. This study aimed to compare glutathione peroxidase (TDPX) gene expression in procyclic and metacyclic and also interspecies in Iranian isolates of L. major and L. tropica.

METHODS

The samples were cultured in Novy-Nicolle-Mc Neal medium to obtain the promastigotes and identified using PCR-RFLP technique. They were then grown in RPMI1640 media for mass cultivation. The expression level of TDPX gene was compared by Real-time PCR.

RESULTS

By comparison of expression level, up-regulation of TDPX gene was observed (5.37 and 2.29 folds) in L. major and L. tropica metacyclic compared to their procyclic, respectively. Moreover, there was no significant difference between procyclic forms of isolates, while 3.05 folds up-regulation in metacyclic was detected in L. major compared L. tropica.

CONCLUSION

Our data provide a foundation for identifying infectivity and high survival related factors in the Leishmania spp. In addition, the results improve our understanding of the molecular basis of metacyclogenesis and development of new potential targets to control or treatment and also, to the identification of species-specific factors contributing to virulence and pathogenicity in the host cells.

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.

Histone deacetylases inhibitors as new potential drugs against Leishmania braziliensis, the main causative agent of new world tegumentary leishmaniasis

The protozoan parasite Leishmania braziliensis is a major causative agent of the neglected tropical diseases Cutaneous and Mucocutaneous Leishmaniases in the New World. There are no vaccines to prevent the infection and the treatment relies on few drugs that often display high toxicity and costs. Thus, chemotherapeutic alternatives are required. Histone Deacetylases (HDACs) are epigenetic enzymes involved in the control of chromatin structure. In this work, we tested an in-house library of 78 hydroxamic acid derivatives as putative inhibitors of L. braziliensis HDACs (HDACi). The compounds were evaluated in relation to the toxicity to the host cell macrophage and to the leishmanicidal effect against L. braziliensis during in vitro infection. Eight HDACi showed significant leishmanicidal effects and the top 5 compounds showed effective concentrations (EC₅₀) in the range of 4.38 to 10.21 μM and selectivity indexes (SI) from of 6 to 21.7. Analyses by Transmission Electron Microscopy (TEM) indicated induction of apoptotic cell death of L. braziliensis amastigotes with a necrotic phenotype. An altered chromatin condensation pattern and cellular disorganization of intracellular amastigotes was also observed. A tight connection between the mitochondrion and nuclear protrusions, presumably of endoplasmic reticulum origin, was found in parasites but not in the host cell. In flow cytometry (FC) analyses, HDACi promoted parasite cell cycle arrest in the G2-M phase and no changes were found in macrophages. In addition, the direct effect of HDACi against the promastigotes showed apoptosis as the main mechanism of cell death. The FC results corroborate the TEM analyses indicating that the HDACi lead to changes in the cell cycle and induction of apoptosis of L. braziliensis. The production of nitric oxide by the infected macrophages was not altered after treatment with the top 5 compounds. Taken together, our results evidenced new HDACi as promising agents for the development of new treatments for American Tegumentary Leishmaniasis caused by L. braziliensis.

DNA-based microarray studies in visceral leishmaniasis: identification of biomarkers for diagnostic, prognostic and drug target for treatment

Visceral leishmaniasis (VL) is one of the major infectious diseases affecting the poorest regions of the world. Current therapy is not very much satisfactory. The alarming rise of drug resistance and the unavailability of an effective vaccine against VL urges research towards identifying new targets or biomarkers for its effective treatment. New technology developments offer some fresh hope in its diagnosis, treatment, and control. DNA microarray approach is now broadly used in parasitology research to facilitate the thoughtful of mechanisms of disease and identification of drug targets and biomarkers for diagnostic and therapeutic development. An electronic search on "VL" and "Microarray" was conducted in Medline and Scopus and papers published in the English mentioning use of DNA microarray on VL were selected and read to write this paper review. Functional analysis and interpretation of microarray results remain very challenging due to the inherent nature of experimental workflows, access, cost, and complexity of data obtained. We have explained and emphasized the use of curate knowledge of microarray in the case of VL for the identification of therapeutic target and biomarker and their selection/implementation in clinical use.

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.

Genome-Wide Proteomics and Phosphoproteomics Analysis of Leishmania spp. During Differentiation

Determining variations in protein abundance and/or posttranslational modification as a function of time or upon induction by a signal in a particular cell type is central to quantitative proteomics. Isobaric labeling methodologies now allow for parallel quantification of proteins at various conditions concurrently or multiplexing in relatively quantitative proteomics workflows. Hence, mapping the protein expression profiles of various developmental stages of Leishmania parasites is possible with high-resolution mass spectrometry. To analyze global changes in protein expression and cellular signaling pathways during Leishmania differentiation and development is possible with a quantitative proteomics approach. The tandem mass tags (TMT) approach provides a chemical labeling method based on the principle of amine reactive tags; the maximum number of conditions that can be multiplexed is 10-plex. We describe herein a detailed method for sample preparation, TMT-labeling, mass spectrometry and data analysis of different developmental stages of Leishmania donovani parasites. This quantitative proteomic approach is useful to study dynamic changes in protein expression levels during L. donovani differentiation, and also allows in-depth analysis of signaling pathways via phosphoproteomics.

In Vitro Culture for Differentiation Simulation of Leishmania spp

In the 1990s my laboratory discovered that Leishmania promastigotes can combine two environmental cues, typical to lysosomes, acidic pH (~5.5) and body temperature (37 °C) into a single signal that induced differentiation. Based on this concept, we modified EARLS-based medium 199 to become an amastigote-specific medium. Shifting promastigotes to this medium followed by incubation in a CO₂ incubator induced differentiation. Axenic amastigotes reach maturation within 5 days, resembling the time it takes in vivo. This chapter provides a complete protocol we developed for L. donovani promastigote-to-amastigote differentiation. This protocol should be useful for both old-world and new-world species of Leishmania.

Ascorbate-Dependent Peroxidase (APX) from Leishmania amazonensis Is a Reactive Oxygen Species-Induced Essential Enzyme That Regulates Virulence

The molecular mechanisms underlying biological differences between two Leishmania species that cause cutaneous disease, L. major and L. amazonensis, are poorly understood. In L. amazonensis, reactive oxygen species (ROS) signaling drives differentiation of nonvirulent promastigotes into forms capable of infecting host macrophages. Tight spatial and temporal regulation of H₂O₂ is key to this signaling mechanism, suggesting a role for ascorbate-dependent peroxidase (APX), which degrades mitochondrial H₂O₂ Earlier studies showed that APX-null L. major parasites are viable, accumulate higher levels of H₂O₂, generate a greater yield of infective metacyclic promastigotes, and have increased virulence. In contrast, we found that in L. amazonensis, the ROS-inducible APX is essential for survival of all life cycle stages. APX-null promastigotes could not be generated, and parasites carrying a single APX allele were impaired in their ability to infect macrophages and induce cutaneous lesions in mice. Similar to what was reported for L. major, APX depletion in L. amazonensis enhanced differentiation of metacyclic promastigotes and amastigotes, but the parasites failed to replicate after infecting macrophages. APX expression restored APX single-knockout infectivity, while expression of catalytically inactive APX drastically reduced virulence. APX overexpression in wild-type promastigotes reduced metacyclogenesis, but enhanced intracellular survival following macrophage infection or inoculation into mice. Collectively, our data support a role for APX-regulated mitochondrial H₂O₂ in promoting differentiation of virulent forms in both L. major and L. amazonensis Our results also uncover a unique requirement for APX-mediated control of ROS levels for survival and successful intracellular replication of L. amazonensis.

An immunoproteomic approach to identifying immunoreactive proteins in Leishmania infantum amastigotes using sera of dogs infected with canine visceral leishmaniasis

Visceral leishmaniasis (VL), the most severe form of leishmaniasis, is caused by Leishmania donovani and Leishmania infantum. The infected dogs with canine visceral leishmaniasis (CVL) are important reservoirs for VL in humans, so the diagnosis, treatment and vaccination of the infected dogs will ultimately decrease the rate of human VL. Proteomics and immunoproteomics techniques have facilitated the introduction of novel drug, vaccine and diagnostic targets. Our immunoproteomic study was conducted to identify new immunoreactive proteins in amastigote form of L. infantum. The strain of L. infantum (MCAN/IR/07/Moheb-gh) was obtained from CVL-infected dogs. J774 macrophage cells were infected with the L. infantum promastigotes. The infected macrophages were ruptured, and pure amastigotes were extracted from the macrophages. After protein extraction, two-dimensional gel electrophoresis was employed for protein separation followed by Western blotting. Western blotting was performed, using symptomatic and asymptomatic sera of the infected dogs with CVL. Thirteen repeatable immunoreactive spots were identified by Matrix-assisted laser desorption/ionization time-of-flight mass spectrometry. Some, including prohibitin, ornithine aminotransferase, annexin A4, and apolipoprotein A-I, have been critically involved in metabolic pathways, survival, and pathogenicity of Leishmania parasites. Further investigations are required to confirm our identified immunoreactive proteins as a biomarker for CVL.

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

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

Comparative transcriptomics in Leishmania braziliensis: disclosing differential gene expression of coding and putative noncoding RNAs across developmental stages

Leishmaniasis is a worldwide public health problem caused by protozoan parasites of the genus Leishmania. Leishmania braziliensis is the most important species responsible for tegumentary leishmaniases in Brazil. An understanding of the molecular mechanisms underlying the success of this parasite is urgently needed. An in-depth study on the modulation of gene expression across the life cycle stages of L. braziliensis covering coding and noncoding RNAs (ncRNAs) was missing and is presented herein. Analyses of differentially expressed (DE) genes revealed that most prominent differences were observed between the transcriptomes of insect and mammalian proliferative forms (6,576 genes). Gene ontology (GO) analysis indicated stage-specific enriched biological processes. A computational pipeline and 5 ncRNA predictors allowed the identification of 11,372 putative ncRNAs. Most of the DE ncRNAs were found between the transcriptomes of insect and mammalian proliferative stages (38%). Of the DE ncRNAs, 295 were DE in all three stages and displayed a wide range of lengths, chromosomal distributions and locations; many of them had a distinct expression profile compared to that of their protein-coding neighbors. Thirty-five putative ncRNAs were submitted to northern blotting analysis, and one or more hybridization-positive signals were observed in 22 of these ncRNAs. This work presents an overview of the L. braziliensis transcriptome and its adjustments throughout development. In addition to determining the general features of the transcriptome at each life stage and the profile of protein-coding transcripts, we identified and characterized a variety of noncoding transcripts. The novel putative ncRNAs uncovered in L. braziliensis might be regulatory elements to be further investigated.

Quantitative RNA Analysis Using RNA-Seq

High-throughput sequencing of cDNA copies of mRNA (RNA-seq) provides a digital readout of mRNA levels over several orders of magnitude, as well as mapping the transcripts to the nucleotide level. Here we describe two different RNA-seq approaches, including one that exploits the 39-nucleotide mini-exon or spliced leader (SL) sequence found at the 5' end of all Leishmania (and other trypanosomatid) mRNAs.

Host-Free Systems for Differentiation of Axenic Leishmania

This chapter describes, in detail, the method our laboratory developed to differentiate L. donovani promastigotes into amastigotes in a host-free culture. This method is based on previous observations that Leishmania promastigotes can combine two environmental signals, typical to lysosomes, acidic pH (~5.5) and body temperature (37 °C), into a signal that induces differentiation. Based on this concept, we have modified medium 199 to make it into an amastigote-specific medium. Shifting promastigotes to this medium, followed by incubation in a CO₂ incubator, induced differentiation. Axenic amastigotes reach maturation within 5 days, resembling the time it takes in vivo. This chapter provides a complete protocol that should be useful for both Old and New World species of Leishmania.

Antileishmanial Activity of Dimeric Flavonoids Isolated from Arrabidaea brachypoda

Leishmaniasis are diseases caused by parasites belonging to Leishmania genus. The treatment with pentavalent antimonials present high toxicity. Secondary line drugs, such as amphotericin B and miltefosine also have a narrow therapeutic index. Therefore, there is an urgent need to develop new drugs to treat leishmaniasis. Here, we present the in vitro anti-leishmanial activity of unusual dimeric flavonoids purified from Arrabidaea brachypoda. Three compounds were tested against Leishmana sp. Compound 2 was the most active against promastigotes. Quantifying the in vitro infected macrophages revealed that compound 2 was also the most active against intracellular amastigotes of L. amazonensis, without displaying host cell toxicity. Drug combinations presented an additive effect, suggesting the absence of interaction between amphotericin B and compound 2. Amastigotes treated with compound 2 demonstrated alterations in the Golgi and accumulation of vesicles inside the flagellar pocket. Compound 2-treated amastigotes presented a high accumulation of cytoplasmic vesicles and a myelin-like structure. When administered in L. amazonensis-infected mice, neither the oral nor the topical treatments were effective against the parasite. Based on the high in vitro activity, dimeric flavonoids can be used as a lead structure for the development of new molecules that could be useful for structure-active studies against Leishmania.

The contribution of DNA microarray technology to gene expression profiling in Leishmania spp.: A retrospective view

The first completed genome project of any living organism, excluding viruses, was of the gammaproteobacteria Haemophilus influenzae in 1995. Until the last decade, genome sequencing was very tedious because genome survey sequences (GSS) and/or expressed sequence tags (ESTs) belonging to plasmid, cosmid, and artificial chromosome genome libraries had to be sequenced and assembled in silico. No genome is completely assembled because gaps and unassembled contigs are always remaining. However, most represent an organism's whole genome from a practical point of view. The first genome sequencing projects of trypanosomatid parasites Leishmania major, Trypanosoma cruzi, and T. brucei were completed in 2005 following those strategies. The functional genomics era developed on the basis of microarray technology and has been continuously evolving. In the case of the genus Leishmania, substantial information about differentiation in the digenetic life cycle of the parasite has been obtained. More recently, next generation sequencing has revolutionized genome sequencing and functional genomics, leading to more sensitive and accurate results by using much fewer resources. Though this new technology is more advantageous, it does not invalidate microarray results. In fact, promising vaccine candidates and drug targets have been found by means of microarray-based screening and preliminary proof-of-concept tests.

Genome Plasticity in Cultured Leishmania donovani: Comparison of Early and Late Passages

Leishmania donovani possesses a complex heteroxenic life cycle where infective metacyclic promastigotes are pre-adapted to infect their host and cope up with intracellular stress. Exploiting the similarities between cultured and sandfly derived promastigotes, we used early and late passage cultured promastigotes to show specific changes at genome level which compromise pathogen fitness reflected in gene expression and infection studies. The pathogen loses virulence mostly via transcriptional and translational regulations and long-time cultivation makes them struggle to convert to virulent metacyclics. At the genomic level very subtle plasticity was observed between the early and the late passages mostly in defense-related, nutrient acquisition and signal transduction genes. Chromosome Copy number variation is seen in the early and late passages involving several genes that may be playing a role in pathogenicity. Our study highlights the importance of ABC transporters and calpain like cysteine proteases in parasite virulence in cultured promastigotes. Interestingly, these proteins are emerging as important patho-adaptive factors in clinical isolates of Leishmania. We found that the currently available genome of Leishmania in the NCBI database are from late passages. Our early passage genome can act as a reference for future studies on virulent isolates of Leishmania. The annotated leads from this study can be used for virulence surveillance and therapeutic studies in the Indian subcontinent.

Comparison of p27 Gene Expression of Promastigote and Amastigote Forms of Leishmania major (MRHO/IR/75/ER) by Real-time RT-PCR

BACKGROUND

Cutaneous leishmaniasis (CL) is one of the world health problems. Leishmania major is the etiological agent of zoonotic cutaneous leishmaniasis (ZCL). Promastigote and amastigote are two morphological forms of Leishmania parasites that express different proteins and p27 is an important gene encoding cytochrome c oxidase (COX) component. P27 gene expresses a 27 kDa protein that essential in ATP synthesis. This study aimed to compare p27 gene expression in promastigote and amastigote forms in Iranian strain of L. major (MRHO/IR/75/ER).

METHODS

This study was conducted in 2015. Clinical isolates of CL patients from north, center, west and south parts of Iran were collected and identified by PCRRFLP. After RNA extraction of promastigotes and amastigotes and cDNA synthesis, the expression level of p27 gene was compared by real-time RT-PCR.

RESULTS

By comparison of expression level between amastigote and promastigote forms of Iranian strain of L. major, up-regulation of p27 gene (2.73 fold) was observed in amastigotes. Moreover, there was no significant difference in p27 gene expression between L. major isolates.

CONCLUSION

p27 gene and protein can be considered as a target in recombinant vaccine production and treatment process.

RNA-seq transcriptional profiling of Leishmania amazonensis reveals an arginase-dependent gene expression regulation

Background

Leishmania is a protozoan parasite that alternates its life cycle between the sand-fly vector and the mammalian host. This alternation involves environmental changes and leads the parasite to dynamic modifications in morphology, metabolism, cellular signaling and regulation of gene expression to allow for a rapid adaptation to new conditions. The L-arginine pathway in L. amazonensis is important during the parasite life cycle and interferes in the establishment and maintenance of the infection in mammalian macrophages. Host arginase is an immune-regulatory enzyme that can reduce the production of nitric oxide by activated macrophages, directing the availability of L-arginine to the polyamine pathway, resulting in parasite replication. In this work, we performed transcriptional profiling to identify differentially expressed genes in L. amazonensis wild-type (La-WT) versus L. amazonensis arginase knockout (La-arg^-) promastigotes and axenic amastigotes.

Methodology/Principal findings

A total of 8253 transcripts were identified in La-WT and La-arg^- promastigotes and axenic amastigotes, about 60% of them codifying hypothetical proteins and 443 novel transcripts, which did not match any previously annotated genes. Our RNA-seq data revealed that 85% of genes were constitutively expressed. The comparison of transcriptome and metabolome data showed lower levels of arginase and higher levels of glutamate-5-kinase in La-WT axenic amastigotes compared to promastigotes. The absence of arginase activity in promastigotes increased the levels of pyrroline 5-carboxylate reductase, but decreased the levels of arginosuccinate synthase, pyrroline 5-carboxylate dehydrogenase, acetylornithine deacetylase and spermidine synthase transcripts levels. These observations can explain previous metabolomic data pointing to the increase of L-arginine, citrulline and L-glutamate and reduction of aspartate, proline, ornithine and putrescine. Altogether, these results indicate that arginase activity is important in Leishmania gene expression modulation during differentiation and adaptation to environmental changes. Here, we confirmed this hypothesis with the identification of differential gene expression of the enzymes involved in biosynthesis of amino acids, arginine and proline metabolism and arginine biosynthesis.

Conclusions/Significance

All data provided information about the transcriptomic profiling and the expression levels of La-WT and La-arg^- promastigotes and axenic amastigotes. These findings revealed the importance of arginase in parasite survival and differentiation, and indicated the existence of a coordinated response in the absence of arginase activity related to arginine and polyamine pathways.

Leishmania are auxotrophic for many essential nutrients, including amino acids. In this way, the parasite needs to uptake the amino acids from the environment. The uptake of amino acids is mediated by amino acid transporters that are unique for Leishmania. As part of polyamine pathway, the arginase converts L-arginine to ornithine and furthermore to putrescine, products which are essential for parasite growth. On the other hand, the absence of arginase activity could alter the metabolism of the parasite to surpass the external signals during the life cycle and the fate of infection. The transcriptional profiling of La-WT and La-arg^- promastigotes and axenic amastigotes revealed 8253 transcripts, 60% encoding hypothetical proteins and 443 novel transcripts. In addition, our data revealed that 85% of the genes were constitutively expressed. Among the 15% (1268 genes) of the differentially expressed genes, we identified genes up- and down-regulated comparing the transcript abundance from different life cycle stages of the parasite and in the presence or absence of arginase. We also combined the transcriptional with metabolic profile that revealed a proportional correlation between enzyme and metabolites in the polyamine pathway. The differentiation of promastigotes to amastigotes alters the expression of enzymes from polyamines biosynthesis, which modulates ornithine, L-glutamate, proline and putrescine levels. In addition, the absence of arginase activity increased the levels of L-arginine, citrulline and L-glutamate and decreased the levels of aspartate, proline, ornithine and putrescine in promastigotes by differential modulation of genes involved in its metabolism. Altogether these data provided additional insights into how Leishmania is able to modulate its biological functions in the presence or absence of arginase activity to survive during environmental changes.

Mining the Proteome of Leishmania donovani for the Development of Novel MHC Class I Restricted Epitope for the Control of Visceral Leishmaniasis

Although, the precise host defence mechanism(s) is not completely understood, T cell-mediated immune responses is believed to play a pivotal role in controlling parasite infection. Here we target the stage dependent over expressed gene. Here, the consensus based computational approach was adopted for the screening of potential major histocompatibility complex class I restricted epitopes. Based on the computational analysis and previously published report, a set 19 antigenic proteins derived from Leishmania donovani were screened for further characterization as vaccine candidates. A total of 49 epitopes were predicted, which revealed a comprehensive binding affinity to the 40 different MHC class I supertypes. Based on the population coverage and HLA cross presentation, nine highly promiscuous epitopes such as LTYDDVWTV (P1), FLFPQRTAL(P2), FLFSNGAVV (P3), YIYNFGIRV (P4), YMTAAFAAL (P5), KLLRPFAPL (P6), FMLGWIVTI (P7), SLFERNKRV (P8), and SVWNRIFTL (P9) which have either a high or an intermediate TAP binding affinity were selected for further analysis. Theoretical population coverage analysis of polytope vaccine (P1-P9) revealed more than 92% population. Stimulation with the cocktail of peptide revealed a proliferative CD8⁺ T cell response and increased IFN-γ production. An upregulated NF-κB activity is thought to be play a pivotal role in T cell proliferation against the selected peptide. The Th1-type cytokine profile (presence of IFN-γ and absence of IL-10) suggests the potentiality of the cocktail of epitope as a subunit vaccine against leishmaniasis. However, the efficiency of these epitopes to trigger other Th1 cytokines and chemokines in a humanized mice model could explore its plausibility as a vaccine candidate. J. Cell. Biochem. 119: 378-391, 2018. © 2017 Wiley Periodicals, Inc.

The Transcriptome of Leishmania major Developmental Stages in Their Natural Sand Fly Vector

The life cycle of the Leishmania parasite in the sand fly vector involves differentiation into several distinctive forms, each thought to represent an adaptation to specific microenvironments in the midgut of the fly. Based on transcriptome sequencing (RNA-Seq) results, we describe the first high-resolution analysis of the transcriptome dynamics of four distinct stages of Leishmania major as they develop in a natural vector, Phlebotomus duboscqi. The early transformation from tissue amastigotes to procyclic promastigotes in the blood-fed midgut was accompanied by the greatest number of differentially expressed genes, including the downregulation of amastins, and upregulation of multiple cell surface proteins, sugar and amino acid transporters, and genes related to glucose metabolism and cell cycle progression. The global changes accompanying post-blood meal differentiation of procyclic promastigotes to the nectomonad and metacyclic stages were less extensive, though each displayed a unique signature. The transcriptome of nectomonads, which has not been studied previously, revealed changes consistent with cell cycle arrest and the upregulation of genes associated with starvation and stress, including autophagic pathways of protein recycling. Maturation to the infective, metacyclic stage was accompanied by changes suggesting preadaptation to the intracellular environment of the mammalian host, demonstrated by the amastigote-like profiles of surface proteins and metabolism-related genes. Finally, a direct comparison between sand fly-derived and culture-derived metacyclics revealed a reassuring similarity between the two forms, with the in vivo forms distinguished mainly by a stronger upregulation of transcripts associated with nutrient stress.

The life cycle of Leishmania parasites in the sand fly vector includes their growth and development as morphologically distinct forms of extracellular promastigotes found within the different microenvironments of the gut. Based on RNA-Seq, we provide here the first high-resolution, transcriptomic analysis of Leishmania insect stages during their cyclical development in vivo, from tissue amastigotes ingested with the blood meal to infective, metacyclic promastigotes that initiate infection in the mammalian host. The most extensive genetic reprogramming occurred during the early transformation of amastigotes to rapidly dividing procyclic promastigotes in the blood-fed midgut, with major changes in the abundance of mRNAs for surface proteins and metabolism. The post-blood meal-adapted nectomonad stage was characterized by the downregulation of cell cycle-related genes and the upregulation of stress- and starvation-related genes. Finally, the transcriptome of metacyclic promastigotes shifted to a more amastigote-like profile, suggesting their preadaptation to the intracellular host environment.

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.

Proteome profiling of the growth phases of Leishmania pifanoi promastigotes in axenic culture reveals differential abundance of immunostimulatory proteins

Leishmaniasis is a term that encompasses a compendium of neglected tropical diseases caused by dimorphic and digenetic protozoan parasites from the genus Leishmania (Kinetoplastida: Trypanosomatidae). The clinical manifestations of neotropical cutaneous leishmaniasis (NCL) caused by Leishmania pifanoi and other species of the "Leishmania mexicana complex" mainly correspond to anergic diffuse cutaneous leishmaniasis (ADCL), which is the origin of considerable morbidity. Despite the outstanding advances in the characterization of the trypanosomatid genomes and proteomes, the biology of this species has been scarcely explored. However, the close relation of L. pifanoi to the sequenced species L. mexicana and others included in the "L. mexicana complex" allowed us to perform a two-dimension electrophoresis (2DE) approach to the promastigote proteome at the differential expression level. Protein identifications were performed by matrix-assisted laser desorption-ionization tandem time-of-flight mass spectrometry (MALDI-TOF/TOF). This insight has revealed similarities and differences between L. pifanoi and other species responsible for cutaneous and visceral leishmaniasis. Interestingly, certain proteins that were previously described as immunostimulatory (elongation factor 1β, trypanothione peroxidase, heat shock protein 70, enolase, GDP-forming succinyl-CoA and aldehyde dehydrogenase) are more abundant in the final growth stages of promastigotes (late-logarithmic and/or stationary phase) in the case of L. pifanoi.

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

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