Telomeric gene deletion and intrachromosomal amplification in antimony-resistant Leishmania

Antimony resistance and gene expression in Leishmania: spotlight on molecular and proteomic aspects

Leishmaniasis is a vector-borne parasitic disease caused by Leishmania parasites with a spectrum of clinical manifestations, ranging from skin lesions to severe visceral complications. Treatment of this infection has been extremely challenging with the concurrent emergence of drug resistance. The differential gene expression and the discrepancies in protein functions contribute to the appearance of 2 distinct phenotypes: resistant and sensitive, but the current diagnostic tools fail to differentiate between them. The identification of gene expression patterns and molecular mechanisms coupled with antimony (Sb) resistance can be leveraged to prompt diagnosis and select the most effective treatment methods. The present study attempts to use comparative expression of Sb resistance-associated genes in resistant and sensitive Leishmania, to disclose their relative abundance in clinical or in vitro selected isolates to gain an understanding of the molecular mechanisms of Sb response/resistance. Data suggest that the analysis of resistance gene expression would verify the Sb resistance or susceptibility only to a certain extent; however, none of the individual expression patterns of the studied genes was diagnostic as a biomarker of Sb response of Leishmania. The findings highlighted will be useful in bridging the knowledge gap and discovering innovative diagnostic tools and novel therapeutic targets.

Field-Deployable Treatments For Leishmaniasis: Intrinsic Challenges, Recent Developments and Next Steps

Leishmaniasis is a neglected tropical disease endemic primarily to low- and middle-income countries, for which there has been inadequate development of affordable, safe, and efficacious therapies. Clinical manifestations of leishmaniasis range from self-healing skin lesions to lethal visceral infection with chances of relapse. Although treatments are available, secondary effects limit their use outside the clinic and negatively impact the quality of life of patients in endemic areas. Other non-medicinal treatments, such as thermotherapies, are limited to use in patients with cutaneous leishmaniasis but not with visceral infection. Recent studies shed light to mechanisms through which Leishmania can persist by hiding in cellular safe havens, even after chemotherapies. This review focuses on exploring the cellular niches that Leishmania parasites may be leveraging to persist within the host. Also, the cellular, metabolic, and molecular implications of Leishmania infection and how those could be targeted for therapeutic purposes are discussed. Other therapies, such as those developed against cancer or for manipulation of the ferroptosis pathway, are proposed as possible treatments against leishmaniasis due to their mechanisms of action. In particular, treatments that target hematopoietic stem cells and monocytes, which have recently been found to be necessary components to sustain the infection and provide a safe niche for the parasites are discussed in this review as potential field-deployable treatments against leishmaniasis.

Antioxidant defence system as a rational target for Chagas disease and Leishmaniasis chemotherapy

Chagas disease and leishmaniasis are neglected tropical diseases caused by the protozoan parasites Trypanosoma cruzi and Leishmania spp., respectively. They are among the most important parasitic diseases, affecting millions of people worldwide, being a considerable global challenge. However, there is no human vaccine available against T. cruzi and Leishmania infections, and their control is based mainly on chemotherapy. Treatments for Chagas disease and leishmaniasis have multiple limitations, mainly due to the high toxicity of the available drugs, long-term treatment protocols, and the occurrence of drug-resistant parasite strains. In the case of Chagas disease, there is still the problem of low cure rates in the chronic stage of the disease. Therefore, new therapeutic agents and novel targets for drug development are urgently needed. Antioxidant defence in Trypanosomatidae is a potential target for chemotherapy because the organisms present a unique mechanism for trypanothione-dependent detoxification of peroxides, which differs from that found in vertebrates. Cellular thiol redox homeostasis is maintained by the biosynthesis and reduction of trypanothione, involving different enzymes that act in concert. This study provides an overview of the antioxidant defence focusing on iron superoxide dismutase A, tryparedoxin peroxidase, and ascorbate peroxidase and how the enzymes play an important role in the defence against oxidative stress and their involvement in drug resistance mechanisms in T. cruzi and Leishmania spp.

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.

Comparative analysis of the transcriptional responses of five Leishmania species to trivalent antimony

BACKGROUND

Leishmaniasis is a neglected tropical disease caused by several species of Leishmania. The resistance phenotype of these parasites depends on the characteristics of each species, which contributes to increased therapeutic failures. Understanding the mechanism used by the parasite to survive under treatment pressure in order to identify potential common and specific therapeutic targets is essential for the control of leishmaniasis. The aim of this study was to investigate the expression profiles and potential shared and specific resistance markers of the main Leishmania species of medical importance [subgenus L. (Leishmania): L. donovani, L. infantum and L. amazonensis; subgenus L. (Viannia): L. panamensis and L. braziliensis)] resistant and sensitive to trivalent stibogluconate (SbIII).

METHODS

We conducted comparative analysis of the transcriptomic profiles (only coding sequences) of lines with experimentally induced resistance to SbIII from biological replicates of five Leishmania species available in the databases of four articles based on ortholog attribution. Simultaneously, we carried out functional analysis of ontology and reconstruction of metabolic pathways of the resulting differentially expressed genes (DEGs).

RESULTS

Resistant lines for each species had differential responses in metabolic processes, compound binding, and membrane components concerning their sensitive counterpart. One hundred and thirty-nine metabolic pathways were found, with the three main pathways comprising cysteine and methionine metabolism, glycolysis, and the ribosome. Differentially expressed orthologous genes assigned to species-specific responses predominated, with 899 self-genes. No differentially expressed genes were found in common among the five species. Two common upregulated orthologous genes were found among four species (L. donovani, L. braziliensis, L. amazonensis, and L. panamensis) related to an RNA-binding protein and the NAD(P)H cytochrome-B5-oxidoreductase complex, associated with transcriptional control and de novo synthesis of linoleic acid, critical mechanisms in resistance to antimonials.

CONCLUSION

Herein, we identified potential species-specific genes related to resistance to SbIII. Therefore, we suggest that future studies consider a treatment scheme that is species-specific. Despite the limitations of our study, this is the first approach toward unraveling the pan-genus genetic mechanisms of resistance in leishmaniasis.

Well-Tolerated Amphotericin B Derivatives That Effectively Treat Visceral Leishmaniasis

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

Mutations in an Aquaglyceroporin as a Proven Marker of Antimony Clinical Resistance in the Parasite Leishmania donovani

BACKGROUND

Antimonial drugs have long been the mainstay to treat visceral leishmaniasis. Their use has been discontinued in the Indian subcontinent because of drug resistance, but they are still clinically useful elsewhere. The goal of this study was to find markers of antimony resistance in Leishmania donovani clinical isolates and validate experimentally their role in resistance.

METHODS

The genomes of sensitive and antimony-resistant clinical isolates were sequenced. The role of a specific gene in contributing to resistance was studied by CRISPR-Cas9-mediated gene editing and intracellular drug sensitivity assays.

RESULTS

Both gene copy number variations and single nucleotide variants were associated with antimony resistance. A homozygous insertion of 2 nucleotides was found in the gene coding for the aquaglyceroporin AQP1 in both resistant isolates. Restoring the wild-type AQP1 open reading frame re-sensitized the 2 independent resistant isolates to antimonials. Alternatively, editing the genome of a sensitive isolate by incorporating the 2-nucleotide insertion in its AQP1 gene led to antimony-resistant parasites.

CONCLUSIONS

Through genomic analysis and CRISPR-Cas9-mediated genome editing we have proven the role of the AQP1 mutations in antimony clinical resistance in L. donovani.

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.

Genomic and Transcriptomic Analysis for Identification of Genes and Interlinked Pathways Mediating Artemisinin Resistance in Leishmania donovani

Current therapy for visceral leishmaniasis (VL), compromised by drug resistance, toxicity, and high cost, demands for more effective, safer, and low-cost drugs. Artemisinin has been found to be an effectual drug alternative in experimental models of leishmaniasis. Comparative genome and transcriptome analysis of in vitro-adapted artesunate-resistant (K133AS-R) and -sensitive wild-type (K133WT) Leishmania donovani parasites was carried out using next-generation sequencing and single-color DNA microarray technology, respectively, to identify genes and interlinked pathways contributing to drug resistance. Whole-genome sequence analysis of K133WT vs. K133AS-R parasites revealed substantial variation among the two and identified 240 single nucleotide polymorphisms (SNPs), 237 insertion deletions (InDels), 616 copy number variations (CNVs) (377 deletions and 239 duplications), and trisomy of chromosome 12 in K133AS-R parasites. Transcriptome analysis revealed differential expression of 208 genes (fold change ≥ 2) in K133AS-R parasites. Functional categorization and analysis of modulated genes of interlinked pathways pointed out plausible adaptations in K133AS-R parasites, such as (i) a dependency on lipid and amino acid metabolism for generating energy, (ii) reduced DNA and protein synthesis leading to parasites in the quiescence state, and (iii) active drug efflux. The upregulated expression of cathepsin-L like protease, amastin-like surface protein, and amino acid transporter and downregulated expression of the gene encoding ABCG2, pteridine receptor, adenylatecyclase-type receptor, phosphoaceylglucosamine mutase, and certain hypothetical proteins are concordant with genomic alterations suggesting their potential role in drug resistance. The study provided an understanding of the molecular basis linked to artemisinin resistance in Leishmania parasites, which may be advantageous for safeguarding this drug for future use.

Antileishmanial Aminopyrazoles: Studies into Mechanisms and Stability of Experimental Drug Resistance

Current antileishmanial treatment is hampered by limitations, such as drug toxicity and the risk of treatment failure, which may be related to parasitic drug resistance. Given the urgent need for novel drugs, the Drugs for Neglected Diseases initiative (DNDi) has undertaken a drug discovery program, which has resulted in the identification of aminopyrazoles, a highly promising antileishmanial chemical series. Multiple experiments have been performed to anticipate the propensity for resistance development.

Current antileishmanial treatment is hampered by limitations, such as drug toxicity and the risk of treatment failure, which may be related to parasitic drug resistance. Given the urgent need for novel drugs, the Drugs for Neglected Diseases initiative (DNDi) has undertaken a drug discovery program, which has resulted in the identification of aminopyrazoles, a highly promising antileishmanial chemical series. Multiple experiments have been performed to anticipate the propensity for resistance development. Resistance selection was performed by successive exposure of Leishmania infantum promastigotes (in vitro) and intracellular amastigotes (both in vitro and in golden Syrian hamsters). The stability of the resistant phenotypes was assessed after passage in mice and Lutzomyia longipalpis sandflies. Whole-genome sequencing (WGS) was performed to identify mutated genes, copy number variations (CNVs), and somy changes. The potential role of efflux pumps (the MDR and MRP efflux pumps) in the development of resistance was assessed by coincubation of aminopyrazoles with specific efflux pump inhibitors (verapamil, cyclosporine, and probenecid). Repeated drug exposure of amastigotes did not result in the emergence of drug resistance either in vitro or in vivo. Selection at the promastigote stage, however, was able to select for parasites with reduced susceptibility (resistance index, 5.8 to 24.5). This phenotype proved to be unstable after in vivo passage in mice and sandflies, suggesting that nonfixed alterations are responsible for the elevated resistance. In line with this, single nucleotide polymorphisms and indels identified by whole-genome sequencing could not be directly linked to the decreased drug susceptibility. Copy number variations were absent, whereas somy changes were detected, which may have accounted for the transient acquisition of resistance. Finally, aminopyrazole activity was not influenced by the MDR and MRP efflux pump inhibitors tested. The selection performed does not suggest the rapid development of resistance against aminopyrazoles in the field. Karyotype changes may confer elevated levels of resistance, but these do not seem to be stable in the vertebrate and invertebrate hosts. MDR/MRP efflux pumps are not likely to significantly impact the activity of the aminopyrazole leads.

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.

Preclinical Gold Complexes as Oral Drug Candidates to Treat Leishmaniasis Are Potent Trypanothione Reductase Inhibitors

The drugs currently used to treat leishmaniases have limitations concerning cost, efficacy, and safety, making the search for new therapeutic approaches urgent. We found that the gold(I)-derived complexes were active against L. infantum and L. braziliensis intracellular amastigotes with IC₅₀ values ranging from 0.5 to 5.5 μM. All gold(I) complexes were potent inhibitors of trypanothione reductase (TR), with enzyme IC₅₀ values ranging from 1 to 7.8 μM. Triethylphosphine-derived complexes enhanced reactive oxygen species (ROS) production and decreased mitochondrial respiration after 2 h of exposure, indicating that gold(I) complexes cause oxidative stress by direct ROS production, by causing mitochondrial damage or by impairing TR activity and thus accumulating ROS. There was no cross-resistance to antimony; in fact, SbR (antimony-resistant mutants) strains were hypersensitive to some of the complexes. BALB/c mice infected with luciferase-expressing L. braziliensis or L. amazonensis and treated orally with 12.5 mg/kg/day of AdT Et (3) or AdO Et (4) presented reduced lesion size and parasite burden, as revealed by bioimaging. The combination of (3) and miltefosine allowed for a 50% reduction in miltefosine treatment time. Complexes 3 and 4 presented favorable pharmacokinetic and toxicity profiles that encourage further drug development studies. Gold(I) complexes are promising antileishmanial agents, with a potential for therapeutic use, including in leishmaniasis caused by antimony-resistant parasites.

Leishmania infantum induces high phagocytic capacity and intracellular nitric oxide production by human proinflammatory monocyte

BACKGROUND

The mechanism of resistance to Sb^III in Leishmania is complex, multifactorial and involves not only biochemical mechanisms, but also other elements, such as the immune system of the host.

OBJECTIVES

In this study, putative changes in the immunological profile of human monocytes infected with wild-type (WT) and antimony (Sb^III)-resistant Leishmania (Viannia) braziliensis and Leishmania (Leishmania) infantum lines were evaluated.

METHODS

Susceptibility assays WT and SbIII-resistant L. braziliensis and L. infantum were performed using lines THP-1 human monocytic lineage. Phagocytic capacity, cytokine profile, intracellular nitric oxide (NO) production and surface carbohydrate residues profile were performed in peripheral blood monocytes by flow cytometry.

FINDINGS

The phagocytic capacity and intracellular NO production by classical (CD14⁺⁺CD16^-) and proinflammatory (CD14⁺⁺CD16⁺) monocytes were higher in the presence of L. infantum lines compared to L. braziliensis lines. The results also highlight proinflammatory monocytes as the cellular subpopulation of major relevance in a phagocytosis event and NO expression. It is important to note that L. infantum induced a proinflammatory cytokine profile characterised by higher levels of TNF-α in culture supernatant than L. braziliensis. Conversely, both Leishmania lines induce high levels of IL-6 in culture supernatant. Analysis of the expression profile of surface carbohydrates showed that L. braziliensis presents 4.3-fold higher expression of galactose(β1,4)N-acetylglucosamine than L. infantum line. Interestingly, the expression level of α-N-acetylgalactosamine residues was 2-fold lower in the Sb^III-resistant L. braziliensis line than its counterpart WT line, indicating differences in surface glycoconjugates between these lines.

MAIN CONCLUSIONS

Our results showed that L. braziliensis and L. infantum induce different innate immune responses and a highly inflammatory profile, which is characteristic of infection by L. infantum, the species associated with visceral disease.

MRPA-independent mechanisms of antimony resistance in Leishmania infantum

Control of both human and canine leishmaniasis is based on a very short list of chemotherapeutic agents, headed by antimonial derivatives (Sb). The utility of these molecules is severely threatened by high rates of drug resistance. The ABC transporter MRPA is one of the few key Sb resistance proteins described to date, whose role in detoxification has been thoroughly studied in Leishmania parasites. Nonetheless, its rapid amplification during drug selection complicates the discovery of other mechanisms potentially involved in Sb resistance. In this study, stepwise drug-resistance selection and next-generation sequencing were combined in the search for novel Sb-resistance mechanisms deployed by parasites when MRPA is abolished by targeted gene disruption. The gene mrpA is not essential in L. infantum, and its disruption leads to an Sb hypersensitive phenotype in both promastigotes and amastigotes. Five independent mrpA^-/- mutants were selected for antimony resistance. These mutants displayed major changes in their ploidy, as well as extrachromosomal linear amplifications of the subtelomeric region of chromosome 23, which includes the genes coding for ABCC1 and ABCC2. Overexpression of ABCC2, but not of ABCC1, resulted in increased Sb tolerance in the mrpA^-/- mutant. SNP analyses revealed three different heterozygous mutations in the gene coding for a serine acetyltransferase (SAT) involved in de novo cysteine synthesis in Leishmania. Overexpression of sat^Q390K, sat^G321R and sat^G325R variants led to a 2-3.2 -fold increase in Sb resistance in mrpA^-/- parasites. Only sat^G321R and sat^G325R induced increased Sb resistance in wild-type parasites. These results reinforce and expand knowledge on the complex nature of Sb resistance in Leishmania parasites.

Global genome diversity of the Leishmania donovani complex

Protozoan parasites of the Leishmania donovani complex - L. donovani and L. infantum - cause the fatal disease visceral leishmaniasis. We present the first comprehensive genome-wide global study, with 151 cultured field isolates representing most of the geographical distribution. L. donovani isolates separated into five groups that largely coincide with geographical origin but vary greatly in diversity. In contrast, the majority of L. infantum samples fell into one globally-distributed group with little diversity. This picture is complicated by several hybrid lineages. Identified genetic groups vary in heterozygosity and levels of linkage, suggesting different recombination histories. We characterise chromosome-specific patterns of aneuploidy and identified extensive structural variation, including known and suspected drug resistance loci. This study reveals greater genetic diversity than suggested by geographically-focused studies, provides a resource of genomic variation for future work and sets the scene for a new understanding of the evolution and genetics of the Leishmania donovani complex.

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.

Environmental Conditions May Shape the Patterns of Genomic Variations in Leishmania panamensis

Due to the absence of transcriptional regulation of gene expression in Leishmania parasites, it is now well accepted that several forms of genomic variations modulate the levels of critical proteins through changes in gene dosage. We previously observed many of these variations in our reference laboratory strain of L. panamensis (PSC-1 strain), including chromosomes with an increased somy and the presence of a putative linear minichromosome derived from chromosome 34. Here, we compared the previously described genomic variations with those occurring after exposure of this strain to increasing concentrations of trivalent antimony (Sb^III), as well as those present in two geographically unrelated clinical isolates of L. panamensis. We observed changes in the somy of several chromosomes, amplifications of several chromosomal regions, and copy number variations in gene arrays after exposure to Sb^III. Occurrence of amplifications potentially beneficial for the Sb-resistant phenotype appears to be associated with the loss of other forms of amplification, such as the linear minichromosome. In contrast, we found no evidence of changes in somy or amplification of relatively large chromosomal regions in the clinical isolates. In these isolates, the predominant amplifications appear to be those that generate genes arrays; however, in many cases, the amplified arrays have a notably higher number of copies than those from the untreated and Sb-treated laboratory samples.

Aneuploidy in yeast: Segregation error or adaptation mechanism?

Aneuploidy is the loss or gain of chromosomes within a genome. It is often detrimental and has been associated with cell death and genetic disorders. However, aneuploidy can also be beneficial and provide a quick solution through changes in gene dosage when cells face environmental stress. Here, we review the prevalence of aneuploidy in Saccharomyces, Candida, and Cryptococcus yeasts (and their hybrid offspring) and analyse associations with chromosome size and specific stressors. We discuss how aneuploidy, a segregation error, may in fact provide a natural route for the diversification of microbes and enable important evolutionary innovations given the right ecological circumstances, such as the colonisation of new environments or the transition from commensal to pathogenic lifestyle. We also draw attention to a largely unstudied cross link between hybridisation and aneuploidy. Hybrid meiosis, involving two divergent genomes, can lead to drastically increased rates of aneuploidy in the offspring due to antirecombination and chromosomal missegregation. Because hybridisation and aneuploidy have both been shown to increase with environmental stress, we believe it important and timely to start exploring the evolutionary significance of their co-occurrence.

Transcriptional responses of Leishmania (Leishmania) amazonensis in the presence of trivalent sodium stibogluconate

Background

In the last decade, resistance to antimonials has become a serious problem due to the emergence of drug-resistant strains. Therefore, understanding the mechanisms used by Leishmania parasites to survive under drug pressure is essential, particularly for species of medical-veterinary importance such as L. amazonensis.

Methods

Here, we used RNA-seq technology to analyse transcriptome profiles and identify global changes in gene expression between antimony-resistant and -sensitive L. amazonensis promastigotes.

Results

A total of 723 differentially expressed genes were identified between resistant and sensitive lines. Comparative transcriptomic analysis revealed that genes encoding proteins involved in metabolism (fatty acids) and stress response, as well as those associated with antimony resistance in other Leishmania species, were upregulated in the antimony-resistant line. Most importantly, we observed upregulation of genes encoding autophagy proteins, suggesting that in the presence of trivalent stibogluconate (Sb^III) L. amazonensis can activate these genes either as a survival strategy or to induce cell death, as has been observed in other parasites.

Conclusions

This work identified global transcriptomic changes in an in vitro-adapted strain in response to Sb^III. Our results provide relevant information to continue understanding the mechanism used by parasites of the subgenus Leishmania (L. amazonensis) to generate an antimony-resistant phenotype.

Electronic supplementary material

The online version of this article (10.1186/s13071-019-3603-8) contains supplementary material, which is available to authorized users.

Major changes in chromosomal somy, gene expression and gene dosage driven by SbIII in Leishmania braziliensis and Leishmania panamensis

Leishmania braziliensis and Leishmania panamensis are two species clinically and epidemiologically important, among others because of their relative resistance to first-line drugs (antimonials). The precise mechanism underlying the ability of these species to survive antimony treatment remains unknown. Therefore, elucidating the pathways mediating drug resistance is essential. We herein experimentally selected resistance to trivalent antimony (Sb^III) in the reference strains of L. braziliensis (MHOM/BR75/M2904) and L. panamensis (MHOM/COL/81L13) and compared whole genome and transcriptome alterations in the culture promastigote stage. The results allowed us to identify differences in somy, copy number variations in some genes related to antimony resistance and large-scale copy number variations (deletions and duplications) in chromosomes with no somy changes. We found mainly in L. braziliensis, a direct relation between the chromosomal/local copy number variation and the gene expression. We identified differentially expressed genes in the resistant lines that are involved in antimony resistance, virulence, and vital biological processes in parasites. The results of this study may be useful for characterizing the genetic mechanisms of these Leishmania species under antimonial pressure, and for clarifying why the parasites are resistant to first-line drug treatments.

Mannosyltransferase (GPI-14) overexpression protects promastigote and amastigote forms of Leishmania braziliensis against trivalent antimony

BACKGROUND

Glycosylphosphatidylinositol is a surface molecule important for host-parasite interactions. Mannosyltransferase (GPI-14) is an essential enzyme for adding mannose on the glycosylphosphatidyl group. This study attempted to overexpress the GPI-14 gene in Leishmania braziliensis to investigate its role in the antimony-resistance phenotype of this parasite.

RESULTS

GPI-14 mRNA levels determined by quantitative real-time PCR (qRT-PCR) showed an increased expression in clones transfected with GPI-14 compared to its respective wild-type line. In order to investigate the expression profile of the surface carbohydrates of these clones, the intensity of the fluorescence emitted by the parasites after concanavalin-A (a lectin that binds to the terminal regions of α-D-mannosyl and α-D-glucosyl residues) treatment was analyzed. The results showed that the clones transfected with GPI-14 express 2.8-fold more mannose and glucose residues than those of the wild-type parental line, indicating effective GPI-14 overexpression. Antimony susceptibility tests using promastigotes showed that clones overexpressing the GPI-14 enzyme are 2.4- and 10.5-fold more resistant to potassium antimonyl tartrate (Sb^III) than the parental non-transfected line. Infection analysis using THP-1 macrophages showed that amastigotes from both GPI-14 overexpressing clones were 3-fold more resistant to Sb^III than the wild-type line.

CONCLUSIONS

Our results suggest the involvement of the GPI-14 enzyme in the Sb^III-resistance phenotype of L. braziliensis.

The LPG1x family from Leishmania major is constituted of rare eukaryotic galactofuranosyltransferases with unprecedented catalytic properties

Galactofuranosyltransferases are poorly described enzymes despite their crucial role in the virulence and the pathogenicity of numerous microorganisms. These enzymes are considered as potential targets for therapeutic action. In addition to the only well-characterised prokaryotic GlfT2 from Mycobacterium tuberculosis, four putative genes in Leishmania major were previously described as potential galactofuranosyltransferases. In this study, we have cloned, over-expressed, purified and fully determined the kinetic parameters of these four eukaryotic enzymes, thus demonstrating their unique potency in catalysing the transfer of the galactofuranosyl moiety into acceptors. Their individual promiscuity revealed to be different, as some of them could efficiently use NDP-pyranoses as donor substrates in addition to the natural UDP-galactofuranose. Such results pave the way for the development of chemoenzymatic synthesis of furanosyl-containing glycoconjugates as well as the design of improved drugs against leishmaniasis.

Genomic insights into virulence mechanisms of Leishmania donovani: evidence from an atypical strain

Background

Leishmaniasis is a neglected tropical disease with diverse clinical phenotypes, determined by parasite, host and vector interactions. Despite the advances in molecular biology and the availability of more Leishmania genome references in recent years, the association between parasite species and distinct clinical phenotypes remains poorly understood. We present a genomic comparison of an atypical variant of Leishmania donovani from a South Asian focus, where it mostly causes cutaneous form of leishmaniasis.

Results

Clinical isolates from six cutaneous leishmaniasis patients (CL-SL); 2 of whom were poor responders to antimony (CL-PR), and two visceral leishmaniasis patients (VL-SL) were sequenced on an Illumina MiSeq platform. Chromosome aneuploidy was observed in both groups but was more frequent in CL-SL. 248 genes differed by 2 fold or more in copy number among the two groups. Genes involved in amino acid use (LdBPK_271940) and energy metabolism (LdBPK_271950), predominated the VL-SL group with the same distribution pattern reflected in gene tandem arrays. Genes encoding amastins were present in higher copy numbers in VL-SL and CL-PR as well as being among predicted pseudogenes in CL-SL. Both chromosome and SNP profiles showed CL-SL and VL-SL to form two distinct groups. While expected heterozygosity was much higher in VL-SL, SNP allele frequency patterns did not suggest potential recent recombination breakpoints. The SNP/indel profile obtained using the more recently generated PacBio sequence did not vary markedly from that based on the standard LdBPK282A1 reference. Several genes previously associated with resistance to antimonials were observed in higher copy numbers in the analysis of CL-PR. H-locus amplification was seen in one cutaneous isolate which however did not belong to the CL-PR group.

Conclusions

The data presented suggests that intra species variations at chromosome and gene level are more likely to influence differences in tropism as well as response to treatment, and contributes to greater understanding of parasite molecular mechanisms underpinning these differences. These findings should be substantiated with a larger sample number and expression/functional studies.

Electronic supplementary material

The online version of this article (10.1186/s12864-018-5271-z) contains supplementary material, which is available to authorized users.

J-binding protein 1 and J-binding protein 2 expression in clinical Leishmania major no response-antimonial isolates

Cutaneous leishmaniasis (CL) is a major disease in many parts of the world. Since no vaccine has been developed, treatment is the best way to control it. In most areas, antimonial resistance whose mechanisms have not been completely understood has been reported. The main aim of this study is gene expression assessing of J-binging protein 1 and J-binding protein 2 in clinical Leishmania major isolates. The patients with CL from central and north Iran were considered for this study. The samples were transferred in RNAlater solution and stored in - 20 °C. RNA extraction and cDNA synthesis were performed. The gene expression analysis was done with SYBR Green real-time PCR using ∆∆CT. Written informed consent forms were filled out by patients, and then, information forms were filled out based on the Helsinki Declaration. Statistical analysis was done with SPSS (16.0; SPSS Inc, Chicago) using independent t test, Shapiro-Wilk, and Pearson's and Spearman's rank correlation coefficients. P ≤ 0.05 was considered significant. The gene expression of JBP1 and JBP2 had no relation with sex and age. The JBP1 gene expression was high in sensitive isolates obtained from north of the country. The JBP2 gene expression was significant in sensitive and no response-antimonial isolates from the north, but no significant differences were detected in sensitive and resistant isolates from central Iran. Differential gene expression of JBP1 and JBP2 in various clinical resistances isolates in different geographical areas shows multifactorial ways of developing resistance in different isolates.

Ascorbate peroxidase overexpression protects Leishmania braziliensis against trivalent antimony effects

Ascorbate peroxidase (APX) is a redox enzyme of the trypanothione pathway that converts hydrogen peroxide (H2O2) into water molecules. In the present study, the APX gene was overexpressed in Leishmania braziliensis to investigate its contribution to the trivalent antimony (SbIII)-resistance phenotype. Western blot results demonstrated that APX-overexpressing parasites had higher APX protein levels in comparison with the wild-type line (LbWTS). APX-overexpressing clones showed an 8-fold increase in the antimony-resistance index over the parental line. In addition, our results indicated that these clones were approximately 1.8-fold more tolerant to H2O2 than the LbWTS line, suggesting that the APX enzyme plays an important role in the defence against oxidative stress. Susceptibility tests revealed that APX-overexpressing L. braziliensis lines were more resistant to isoniazid, an antibacterial agent that interacts with APX. Interestingly, this compound enhanced the anti-leishmanial SbIII effect, indicating that this combination represents a good strategy for leishmaniasis chemotherapy. Our data demonstrate that APX enzyme is involved in the development of L. braziliensis antimony-resistance phenotype and may be an attractive therapeutic target in the design of new strategies for leishmaniasis treatment.

Leishmania Genome Dynamics during Environmental Adaptation Reveal Strain-Specific Differences in Gene Copy Number Variation, Karyotype Instability, and Telomeric Amplification

Protozoan parasites of the genus Leishmania adapt to environmental change through chromosome and gene copy number variations. Only little is known about external or intrinsic factors that govern Leishmania genomic adaptation. Here, by conducting longitudinal genome analyses of 10 new Leishmania clinical isolates, we uncovered important differences in gene copy number among genetically highly related strains and revealed gain and loss of gene copies as potential drivers of long-term environmental adaptation in the field. In contrast, chromosome rather than gene amplification was associated with short-term environmental adaptation to in vitro culture. Karyotypic solutions were highly reproducible but unique for a given strain, suggesting that chromosome amplification is under positive selection and dependent on species- and strain-specific intrinsic factors. We revealed a progressive increase in read depth towards the chromosome ends for various Leishmania isolates, which may represent a nonclassical mechanism of telomere maintenance that can preserve integrity of chromosome ends during selection for fast in vitro growth. Together our data draw a complex picture of Leishmania genomic adaptation in the field and in culture, which is driven by a combination of intrinsic genetic factors that generate strain-specific phenotypic variations, which are under environmental selection and allow for fitness gain.IMPORTANCE Protozoan parasites of the genus Leishmania cause severe human and veterinary diseases worldwide, termed leishmaniases. A hallmark of Leishmania biology is its capacity to adapt to a variety of unpredictable fluctuations inside its human host, notably pharmacological interventions, thus, causing drug resistance. Here we investigated mechanisms of environmental adaptation using a comparative genomics approach by sequencing 10 new clinical isolates of the L. donovani, L. major, and L. tropica complexes that were sampled across eight distinct geographical regions. Our data provide new evidence that parasites adapt to environmental change in the field and in culture through a combination of chromosome and gene amplification that likely causes phenotypic variation and drives parasite fitness gains in response to environmental constraints. This novel form of gene expression regulation through genomic change compensates for the absence of classical transcriptional control in these early-branching eukaryotes and opens new venues for biomarker discovery.

Reverse Epidemiology: An Experimental Framework to Drive Leishmania Biomarker Discovery in situ by Functional Genetic Screening Using Relevant Animal Models

Leishmania biomarker discovery remains an important challenge that needs to be revisited in light of our increasing knowledge on parasite-specific biology, notably its genome instability. In the absence of classical transcriptional regulation in these early-branching eukaryotes, fluctuations in transcript abundance can be generated by gene and chromosome amplifications, which have been linked to parasite phenotypic variability with respect to virulence, tissue tropism, and drug resistance. Conducting in vitro evolutionary experiments to study mechanisms of Leishmania environmental adaptation, we recently validated the link between parasite genetic amplification and fitness gain, thus defining gene and chromosome copy number variations (CNVs) as important Leishmania biomarkers. These experiments also demonstrated that long-term Leishmania culture adaptation can strongly interfere with epidemiologically relevant, genetic signals, which challenges current protocols for biomarker discovery, all of which rely on in vitro expansion of clinical isolates. Here we propose an experimental framework independent of long-term culture termed “reverse” epidemiology, which applies established protocols for functional genetic screening of cosmid-transfected parasites in animal models for the identification of clinically relevant genetic loci that then inform targeted field studies for their validation as Leishmania biomarkers.

Molecular Preadaptation to Antimony Resistance in Leishmania donovani on the Indian Subcontinent

The “antibiotic resistance crisis” is a major challenge for scientists and medical professionals. This steady rise in drug-resistant pathogens also extends to parasitic diseases, with antimony being the first anti-Leishmania drug that fell in the Indian subcontinent (ISC). Leishmaniasis is a major but neglected infectious disease with limited therapeutic options. Therefore, understanding how parasites became resistant to antimonials is of commanding importance. In this study, we experimentally characterized the dynamics of this resistance acquisition and show for the first time that some Leishmania populations of the ISC were preadapted to antimony resistance, likely driven by environmental factors or by drugs used in the 19th century.

Antimonials (Sb) were used for decades for chemotherapy of visceral leishmaniasis (VL). Now abandoned in the Indian subcontinent (ISC) because of Leishmania donovani resistance, this drug offers a unique model for understanding drug resistance dynamics. In a previous phylogenomic study, we found two distinct populations of L. donovani: the core group (CG) in the Gangetic plains and ISC1 in the Nepalese highlands. Sb resistance was only encountered within the CG, and a series of potential markers were identified. Here, we analyzed the development of resistance to trivalent antimonials (Sb^III) upon experimental selection in ISC1 and CG strains. We observed that (i) baseline Sb^III susceptibility of parasites was higher in ISC1 than in the CG, (ii) time to Sb^III resistance was higher for ISC1 parasites than for CG strains, and (iii) untargeted genomic and metabolomic analyses revealed molecular changes along the selection process: these were more numerous in ISC1 than in the CG. Altogether these observations led to the hypothesis that CG parasites are preadapted to Sb^III resistance. This hypothesis was experimentally confirmed by showing that only wild-type CG strains could survive a direct exposure to the maximal concentration of Sb^III. The main driver of this preadaptation was shown to be MRPA, a gene involved in Sb^III sequestration and amplified in an intrachromosomal amplicon in all CG strains characterized so far. This amplicon emerged around 1850 in the CG, well before the implementation of antimonials for VL chemotherapy, and we discuss here several hypotheses of selective pressure that could have accompanied its emergence.

IMPORTANCE The “antibiotic resistance crisis” is a major challenge for scientists and medical professionals. This steady rise in drug-resistant pathogens also extends to parasitic diseases, with antimony being the first anti-Leishmania drug that fell in the Indian subcontinent (ISC). Leishmaniasis is a major but neglected infectious disease with limited therapeutic options. Therefore, understanding how parasites became resistant to antimonials is of commanding importance. In this study, we experimentally characterized the dynamics of this resistance acquisition and show for the first time that some Leishmania populations of the ISC were preadapted to antimony resistance, likely driven by environmental factors or by drugs used in the 19th century.

Deletion of transketolase triggers a stringent metabolic response in promastigotes and loss of virulence in amastigotes of Leishmania mexicana

Transketolase (TKT) is part of the non-oxidative branch of the pentose phosphate pathway (PPP). Here we describe the impact of removing this enzyme from the pathogenic protozoan Leishmania mexicana. Whereas the deletion had no obvious effect on cultured promastigote forms of the parasite, the Δtkt cells were not virulent in mice. Δtkt promastigotes were more susceptible to oxidative stress and various leishmanicidal drugs than wild-type, and metabolomics analysis revealed profound changes to metabolism in these cells. In addition to changes consistent with those directly related to the role of TKT in the PPP, central carbon metabolism was substantially decreased, the cells consumed significantly less glucose, flux through glycolysis diminished, and production of the main end products of metabolism was decreased. Only minor changes in RNA abundance from genes encoding enzymes in central carbon metabolism, however, were detected although fructose-1,6-bisphosphate aldolase activity was decreased two-fold in the knock-out cell line. We also showed that the dual localisation of TKT between cytosol and glycosomes is determined by the C-terminus of the enzyme and by engineering different variants of the enzyme we could alter its sub-cellular localisation. However, no effect on the overall flux of glucose was noted irrespective of whether the enzyme was found uniquely in either compartment, or in both.

Leishmania parasites endanger over 1 billion people worldwide, infecting 300,000 people and causing 20,000 deaths annually. In this study, we scrutinized metabolism in Leishmania mexicana after deletion of the gene encoding transketolase (TKT), an enzyme involved in sugar metabolism via the pentose phosphate pathway which plays key roles in creating ribose 5-phosphate for nucleotide synthesis and also defence against oxidative stress. The insect stage of the parasite, grown in culture medium, did not suffer from any obvious growth defect after the gene was deleted. However, its metabolism changed dramatically, with metabolomics indicating profound changes to flux through the pentose phosphate pathway: decreased glucose consumption, and generally enhanced efficiency in using metabolic substrates with reduced secretion of partially oxidised end products of metabolism. This ‘stringent’ metabolism is reminiscent of the mammalian stage parasites. The cells were also more sensitive to oxidative stress inducing agents and leishmanicidal drugs. Crucially, mice inoculated with the TKT knock-out parasites did not develop an infection pointing to the enzyme playing a key role in allowing the parasites to remain viable in the host, indicating that TKT may be considered a useful target for development of new drugs against leishmaniasis.

The enemy within: Targeting host-parasite interaction for antileishmanial drug discovery

The state of antileishmanial chemotherapy is strongly compromised by the emergence of drug-resistant Leishmania. The evolution of drug-resistant phenotypes has been linked to the parasites’ intrinsic genome instability, with frequent gene and chromosome amplifications causing fitness gains that are directly selected by environmental factors, including the presence of antileishmanial drugs. Thus, even though the unique eukaryotic biology of Leishmania and its dependence on parasite-specific virulence factors provide valid opportunities for chemotherapeutical intervention, all strategies that target the parasite in a direct fashion are likely prone to select for resistance. Here, we review the current state of antileishmanial chemotherapy and discuss the limitations of ongoing drug discovery efforts. We finally propose new strategies that target Leishmania viability indirectly via mechanisms of host–parasite interaction, including parasite-released ectokinases and host epigenetic regulation, which modulate host cell signaling and transcriptional regulation, respectively, to establish permissive conditions for intracellular Leishmania survival.

Glucose-6-phosphate dehydrogenase and Trypanothione reductase interaction protects Leishmania donovani from metalloid mediated oxidative stress

Exploration of metabolons as viable drug target is rare in kinetoplastid biology. Here we present a novel protein-protein interaction among Glucose-6-phosphate dehydrogenase (LdG6PDH) and Trypanothione reductase (LdTryR) of Leishmania donovani displaying interconnection between central glucose metabolism and thiol metabolism of this parasite. Digitonin fractionation patterns observed through immunoblotting indicated localisation of both LdG6PDH and LdTryR in cytosol. In-silico and in-vitro interaction observed by size exclusion chromatography, co-purification, pull-down assay and spectrofluorimetric analysis revealed LdG6PDH and LdTryR physically interact with each other in a NADPH dependent manner. Coupled enzymatic assay displayed that NADPH generation was severely impaired by addition of Sb^III, As^III and Te^IV extraneously, which hint towards metalloid driven structural changes of the interacting proteins. Co-purification patterns and pull-down assays also depicted that metalloids (Sb^III, As^III and Te^IV) hinder the in-vitro interaction of these two enzymes. Surprisingly, metalloids at sub-lethal concentrations induced the in-vivo interaction of LdG6PDH and LdTryR, as analyzed by pull-down assays and fluorescence microscopy signifying protection against metalloid mediated ROS. Inhibition of LdTryR by thioridazine in LdG6PDH^-/- parasites resulted in metalloid induced apoptotic death of the parasites due to abrupt fall in reduced thiol content, disrupted NADPH/NADP⁺ homeostasis and lethal oxidative stress. Interestingly, clinical isolates of L.donovani resistant to SAG exhibited enhanced interaction between LdG6PDH and LdTryR and showed cross resistivity towards As^III and Te^IV. Thus, our findings propose the metabolon of LdG6PDH and LdTryR as an alternate therapeutic target and provide mechanistic insight about metalloid resistance in Visceral Leishmaniasis.

Exploiting Knowledge on Leishmania Drug Resistance to Support the Quest for New Drugs

New drugs are needed to control leishmaniasis and efforts are currently on-going to counter the neglect of this disease. We discuss here the utility and the impact of associating drug resistance (DR) studies to drug discovery pipelines. We use as paradigm currently used drugs, antimonials and miltefosine, and complement our reflection by interviewing three experts in the field. We suggest DR studies to be involved at two different stages of drug development: (i) the efficiency of novel compounds should be confirmed on sets of strains including recent clinical isolates with DR; (ii) experimental DR should be generated to promising compounds at an early stage of their development, to further optimize them and monitor clinical trials.

Comparative genomics of Tunisian Leishmania major isolates causing human cutaneous leishmaniasis with contrasting clinical severity

Zoonotic cutaneous leishmaniasis caused by Leishmania (L.) major parasites affects urban and suburban areas in the center and south of Tunisia where the disease is endemo-epidemic. Several cases were reported in human patients for which infection due to L. major induced lesions with a broad range of severity. However, very little is known about the mechanisms underlying this diversity. Our hypothesis is that parasite genomic variability could, in addition to the host immunological background, contribute to the intra-species clinical variability observed in patients and explain the lesion size differences observed in the experimental model. Based on several epidemiological, in vivo and in vitro experiments, we focused on two clinical isolates showing contrasted severity in patients and BALB/c experimental mice model. We used DNA-seq as a high-throughput technology to facilitate the identification of genetic variants with discriminating potential between both isolates. Our results demonstrate that various levels of heterogeneity could be found between both L. major isolates in terms of chromosome or gene copy number variation (CNV), and that the intra-species divergence could surprisingly be related to single nucleotide polymorphisms (SNPs) and Insertion/Deletion (InDels) events. Interestingly, we particularly focused here on genes affected by both types of variants and correlated them with the observed gene CNV. Whether these differences are sufficient to explain the severity in patients is obviously still open to debate, but we do believe that additional layers of -omic information is needed to complement the genomic screen in order to draw a more complete map of severity determinants.

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

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

Silver and Nitrate Oppositely Modulate Antimony Susceptibility through Aquaglyceroporin 1 in Leishmania (Viannia) Species

Antimony (Sb) resistance in leishmaniasis chemotherapy has become one of the major challenges to the control of this spreading worldwide public health problem. Since the plasma membrane pore-forming protein aquaglyceroporin 1 (AQP1) is the major route of Sb uptake in Leishmania, functional studies are relevant to characterize drug transport pathways in the parasite. We generated AQP1-overexpressing Leishmania guyanensis and L. braziliensis mutants and investigated their susceptibility to the trivalent form of Sb (Sb(III)) in the presence of silver and nitrate salts. Both AQP1-overexpressing lines presented 3- to 4-fold increased AQP1 expression levels compared with those of their untransfected counterparts, leading to an increased Sb(III) susceptibility of about 2-fold. Competition assays using silver nitrate, silver sulfadiazine, or silver acetate prior to Sb(III) exposure increased parasite growth, especially in AQP1-overexpressing mutants. Surprisingly, Sb(III)-sodium nitrate or Sb(III)-potassium nitrate combinations showed significantly enhanced antileishmanial activities compared to those of Sb(III) alone, especially against AQP1-overexpressing mutants, suggesting a putative nitrate-dependent modulation of AQP1 activity. The intracellular level of antimony quantified by graphite furnace atomic absorption spectrometry showed that the concomitant exposure to Sb(III) and nitrate favors antimony accumulation in the parasite, increasing the toxicity of the drug and culminating with parasite death. This is the first report showing evidence of AQP1-mediated Sb(III) susceptibility modulation by silver in Leishmania and suggests the potential antileishmanial activity of the combination of nitrate salts and Sb(III).

Chromosomal Translocations in the Parasite Leishmania by a MRE11/RAD50-Independent Microhomology-Mediated End Joining Mechanism

The parasite Leishmania often relies on gene rearrangements to survive stressful environments. However, safeguarding a minimum level of genome integrity is important for cell survival. We hypothesized that maintenance of genomic integrity in Leishmania would imply a leading role of the MRE11 and RAD50 proteins considering their role in DNA repair, chromosomal organization and protection of chromosomes ends in other organisms. Attempts to generate RAD50 null mutants in a wild-type background failed and we provide evidence that this gene is essential. Remarkably, inactivation of RAD50 was possible in a MRE11 null mutant that we had previously generated, providing good evidence that RAD50 may be dispensable in the absence of MRE11. Inactivation of the MRE11 and RAD50 genes led to a decreased frequency of homologous recombination and analysis of the null mutants by whole genome sequencing revealed several chromosomal translocations. Sequencing of the junction between translocated chromosomes highlighted microhomology sequences at the level of breakpoint regions. Sequencing data also showed a decreased coverage at subtelomeric locations in many chromosomes in the MRE11^-/-RAD50^-/- parasites. This study demonstrates an MRE11-independent microhomology-mediated end-joining mechanism and a prominent role for MRE11 and RAD50 in the maintenance of genomic integrity. Moreover, we suggest the possible involvement of RAD50 in subtelomeric regions stability.

The parasite Leishmania relies on gene rearrangements to survive stressful conditions. However, maintaining a minimum level of genomic integrity is crucial for cell survival. Studies in other organisms have provided evidence that the DNA repair proteins MRE11 and RAD50 are involved in chromosomes organization, protection of chromosomes ends and therefore in the maintenance of genomic integrity. In this manuscript, we present the conditional inactivation of the Leishmania infantum RAD50 gene that was only possible in MRE11 deficient cells and suggest the genetic background is crucial for RAD50 inactivation. We demonstrate the occurrence of chromosomal translocations in the MRE11 and RAD50 deficient cells and described a MRE11-independent microhomology-mediated end-joining mechanism at the level of translocation breakpoints. We also suggest a possible involvement of RAD50 in subtelomeric regions stability. Our results highlight that both MRE11 and RAD50 are important for the maintenance of genomic integrity in Leishmania.

Evolutionary genomics of epidemic visceral leishmaniasis in the Indian subcontinent

Leishmania donovani causes visceral leishmaniasis (VL), the second most deadly vector-borne parasitic disease. A recent epidemic in the Indian subcontinent (ISC) caused up to 80% of global VL and over 30,000 deaths per year. Resistance against antimonial drugs has probably been a contributing factor in the persistence of this epidemic. Here we use whole genome sequences from 204 clinical isolates to track the evolution and epidemiology of L. donovani from the ISC. We identify independent radiations that have emerged since a bottleneck coincident with 1960s DDT spraying campaigns. A genetically distinct population frequently resistant to antimonials has a two base-pair insertion in the aquaglyceroporin gene LdAQP1 that prevents the transport of trivalent antimonials. We find evidence of genetic exchange between ISC populations, and show that the mutation in LdAQP1 has spread by recombination. Our results reveal the complexity of L. donovani evolution in the ISC in response to drug treatment.

DOI:http://dx.doi.org/10.7554/eLife.12613.001

The parasite Leishmania donovani causes a disease called visceral leishmaniasis that affects many of the world's poorest people. Around half a million new cases develop every year, but health authorities lack safe and effective drugs to treat them. Up to 80% of these cases occur in the Indian subcontinent, where devastating epidemics have occurred in the last decades.

One reason these epidemics continue to occur is that the parasites develop genetic mutations allowing them to adapt to and resist the drugs used to kill them. As there are few existing drugs that can kill L. donovani, it is crucial to understand how drug resistance emerges and spreads among parasite populations.

Imamura, Downing, Van den Broeck et al. have now investigated the history of visceral leishmaniasis epidemics by characterising the complete genetic sequence – or genome – of 204 L. donovani parasite samples. This revealed that the majority of parasites in the Indian subcontinent first appeared in the nineteenth century, matching the first historical records of visceral leishmaniasis epidemics.

The genomes show that most of the parasites are genetically similar and can be clustered into several closely related groups. These groups first appeared in the 1960s following the end of a regional campaign to eradicate malaria. The most common parasite group is particularly resistant to drugs called antimonials, which were the main treatment for leishmaniasis until recently. These parasites have a small genetic change that scrambles most of a protein known to be involved in the uptake of antimonials.

Parasites may also be able to develop resistance to drugs through additional mechanisms that allow them to produce many copies of the same gene. These mechanisms could allow the parasites to rapidly adapt to new drugs or changes in the populations it infects. The work of Imamura et al. looks only at parasites isolated from patients then grown in the laboratory, so further research is now needed to explore how variable the Leishmania genome is in both of the parasite’s hosts: humans and sandflies.

Imamura et al.’s study reveals how L. donovani has spread throughout the Indian subcontinent in fine detail. The genome data can be used to create simple molecular tools that could form an "early warning system" to track the success of disease control programs and to determine how well the current drugs are working.

DOI:http://dx.doi.org/10.7554/eLife.12613.002

Generation of an aquaglyceroporin AQP1 null mutant in Leishmania major

The Leishmania aquaglyceroporin AQP1 plays an important physiological role in water and uncharged polar solutes transport, volume regulation, osmotaxis, and is a key determinant of antimony resistance. By targeted gene disruption, we generated a Leishmania major promastigote AQP1 null mutant. This required several attempts but a chromosomal null AQP1 mutant was obtained by loss of heterozygosity in the presence of a rescue plasmid encoding AQP1. Growth in the absence of selection led to the loss of the rescuing plasmid, indicating that AQP1 is not essential for Leishmania viability. The AQP1-null mutant was resistant to antimonyl tartrate (SbIII) and arsenite (AsIII) due to a decrease import of these metalloids. It also exhibited alterations in its osmoregulation abilities compared with wild-type cells. This is the first report of the generation of a genetic AQP1 null mutant in Leishmania parasite, confirming its physiological function and role in resistance to antimonials, the therapeutic mainstay against Leishmania.

Metabolic reconfiguration of the central glucose metabolism: a crucial strategy of Leishmania donovani for its survival during oxidative stress

Understanding the mechanism that allows the intracellular protozoan parasite Leishmania donovani (Ld) to respond to reactive oxygen species (ROS) is of increasing therapeutic importance because of the continuing resistance toward antileishmanial drugs and for determining the illusive survival strategy of these parasites. A shift in primary carbon metabolism is the fastest response to oxidative stress. A (14)CO2 evolution study, expression of glucose transporters together with consumption assays, indicated a shift in metabolic flux of the parasites from glycolysis toward pentose phosphate pathway (PPP) when exposed to different oxidants in vitro/ex vivo. Changes in gene expression, protein levels, and enzyme activities all pointed to a metabolic reconfiguration of the central glucose metabolism in response to oxidants. Generation of glucose-6-phosphate dehydrogenase (G6PDH) (∼5-fold) and transaldolase (TAL) (∼4.2-fold) overexpressing Ld cells reaffirmed that lethal doses of ROS were counterbalanced by effective manipulation of NADPH:NADP(+) ratio and stringent maintenance of reduced thiol content. The extent of protein carbonylation and accumulation of lipid peroxidized products were also found to be less in overexpressed cell lines. Interestingly, the LD50 of sodium antimony gluconate (SAG), amphotericin-B (AmB), and miltefosine were significantly high toward overexpressing parasites. Consequently, this study illustrates that Ld strategizes a metabolic reconfiguration for replenishment of NADPH pool to encounter oxidative challenges.

Intrachromosomal amplification, locus deletion and point mutation in the aquaglyceroporin AQP1 gene in antimony resistant Leishmania (Viannia) guyanensis

BACKGROUND

Antimony resistance complicates the treatment of infections caused by the parasite Leishmania.

METHODOLOGY/PRINCIPAL FINDINGS

Using next generation sequencing, we sequenced the genome of four independent Leishmania guyanensis antimony-resistant (SbR) mutants and found different chromosomal alterations including aneuploidy, intrachromosomal gene amplification and gene deletion. A segment covering 30 genes on chromosome 19 was amplified intrachromosomally in three of the four mutants. The gene coding for the multidrug resistance associated protein A involved in antimony resistance was also amplified in the four mutants, most likely through chromosomal translocation. All mutants also displayed a reduced accumulation of antimony mainly due to genomic alterations at the level of the subtelomeric region of chromosome 31 harboring the gene coding for the aquaglyceroporin 1 (LgAQP1). Resistance involved the loss of LgAQP1 through subtelomeric deletions in three mutants. Interestingly, the fourth mutant harbored a single G133D point mutation in LgAQP1 whose role in resistance was functionality confirmed through drug sensitivity and antimony accumulation assays. In contrast to the Leishmania subspecies that resort to extrachromosomal amplification, the Viannia strains studied here used intrachromosomal amplification and locus deletion.

CONCLUSIONS/SIGNIFICANCE

This is the first report of a naturally occurred point mutation in AQP1 in antimony resistant parasites.

Drug resistance analysis by next generation sequencing in Leishmania

•

WGS revealed the complexity of resistance which is often polyclonal in Leishmania.

•

An impressive variation in chromosome copy numbers exist between Leishmania species.

•

The genotype heterogeneity complicates the analysis of resistance in field isolates.

The use of next generation sequencing has the power to expedite the identification of drug resistance determinants and biomarkers and was applied successfully to drug resistance studies in Leishmania. This allowed the identification of modulation in gene expression, gene dosage alterations, changes in chromosome copy numbers and single nucleotide polymorphisms that correlated with resistance in Leishmania strains derived from the laboratory and from the field. An impressive heterogeneity at the population level was also observed, individual clones within populations often differing in both genotypes and phenotypes, hence complicating the elucidation of resistance mechanisms. This review summarizes the most recent highlights that whole genome sequencing brought to our understanding of Leishmania drug resistance and likely new directions.

Comparative analysis of the omics technologies used to study antimonial, amphotericin B, and pentamidine resistance in leishmania

Leishmaniasis is a serious threat in developing countries due to its endemic nature and debilitating symptoms. Extensive research and investigations have been carried out to learn about the mechanism of drug resistance in Leishmania but results obtained in the laboratory are not in agreement with those obtained from the field. Also the lack of knowledge about the mode of action for a number of drugs makes the study of drug resistance more complex. A major concern in recent times has been regarding the role of parasitic virulence in drug resistance for Leishmania. Researchers have employed various techniques to unravel the facts about resistance and virulence in Leishmania. With advent of advanced and more specific means of detection, further hints about probable mechanisms of conferring resistance are expected. This review aims to provide a consolidated picture along with a comparative account of the work done so far to study the mechanism of antimony, amphotericin B, and pentamidine resistance using various techniques.

Metabolic reprogramming during purine stress in the protozoan pathogen Leishmania donovani

The ability of Leishmania to survive in their insect or mammalian host is dependent upon an ability to sense and adapt to changes in the microenvironment. However, little is known about the molecular mechanisms underlying the parasite response to environmental changes, such as nutrient availability. To elucidate nutrient stress response pathways in Leishmania donovani, we have used purine starvation as the paradigm. The salvage of purines from the host milieu is obligatory for parasite replication; nevertheless, purine-starved parasites can persist in culture without supplementary purine for over three months, indicating that the response to purine starvation is robust and engenders parasite survival under conditions of extreme scarcity. To understand metabolic reprogramming during purine starvation we have employed global approaches. Whole proteome comparisons between purine-starved and purine-replete parasites over a 6–48 h span have revealed a temporal and coordinated response to purine starvation. Purine transporters and enzymes involved in acquisition at the cell surface are upregulated within a few hours of purine removal from the media, while other key purine salvage components are upregulated later in the time-course and more modestly. After 48 h, the proteome of purine-starved parasites is extensively remodeled and adaptations to purine stress appear tailored to deal with both purine deprivation and general stress. To probe the molecular mechanisms affecting proteome remodeling in response to purine starvation, comparative RNA-seq analyses, qRT-PCR, and luciferase reporter assays were performed on purine-starved versus purine-replete parasites. While the regulation of a minority of proteins tracked with changes at the mRNA level, for many regulated proteins it appears that proteome remodeling during purine stress occurs primarily via translational and/or post-translational mechanisms.

Leishmania, the cause of a deadly spectrum of diseases in humans, surmounts a number of environmental challenges, including changes in the availability of salvageable nutrients, to successfully colonize its host. Adaptation to environmental stress is clearly of significance in parasite biology, but the underlying mechanisms are not well understood. To simulate the response to periodic nutrient scarcity in vivo, we have induced purine starvation in vitro. Purines are essential for growth and viability, and serve as the major energy currency of cells. Leishmania cannot synthesize purines and must salvage them from the surroundings. Extracellular purine depletion in culture induces a robust survival response in Leishmania, whereby growth arrests, but parasites persist for months. To profile the events that enable endurance of purine starvation, we used shotgun proteomics. Our data suggest that purine starvation induces extensive proteome remodeling, tailored to enhance purine capture and recycling, reduce energy expenditures, and maintain viability of the metabolically active, non-dividing population. Through global and targeted approaches, we reveal that proteome remodeling is multifaceted, and occurs through an array of responses at the mRNA, translational, and post-translational level. Our data provide one of the most inclusive views of adaptation to microenvironmental stress in Leishmania.

Antimony transport mechanisms in resistant leishmania parasites

Antimonial compounds have been used for more than a century in the treatment of the parasitic disease leishmaniasis. Although pentavalent antimonials are still first-line drugs in several developing countries, this class of drugs is no longer recommended in the Indian sub-continent because of the emergence of drug resistance. The precise mechanisms involved in the resistance of leishmania parasites to antimony are still subject to debate. It is now well documented that drug resistance in leishmania parasites is a multifactorial phenomenon involving multiple genes whose expression pattern synergistically leads to the resistance phenotype. The reduction of intracellular antimony accumulation is a frequent change observed in resistant leishmania cells; however, no comprehensive transport model has been presented so far to explain this change and its contribution to Leishmania resistance. The present review firstly covers the actual knowledge on the metabolism of antimonial drugs, the mechanisms of their transmembrane transport and intracellular processing in Leishmania. It further describes both the functional and molecular changes associated with Sb resistance in this organism. Possible transport models based on the actual knowledge are then presented, as well as their functional implications. Biophysical and pharmacological strategies are finally proposed for the precise identification of the transport pathways.

Proteomic and genomic analyses of antimony resistant Leishmania infantum mutant

BACKGROUND

Antimonials remain the primary antileishmanial drugs in most developing countries. However, drug resistance to these compounds is increasing and our understanding of resistance mechanisms is partial.

METHODS/PRINCIPAL FINDINGS

In the present study, quantitative proteomics using stable isotope labelling of amino acids in cell culture (SILAC) and genome next generation sequencing were used in order to better characterize in vitro generated Leishmania infantum antimony resistant mutant (Sb2000.1). Using the proteomic method, 58 proteins were found to be differentially regulated in Sb2000.1. The ABC transporter MRPA (ABCC3), a known marker of antimony resistance, was observed for the first time in a proteomic screen. Furthermore, transfection of its gene conferred antimony resistance in wild-type cells. Next generation sequencing revealed aneuploidy for 8 chromosomes in Sb2000.1. Moreover, specific amplified regions derived from chromosomes 17 and 23 were observed in Sb2000.1 and a single nucleotide polymorphism (SNP) was detected in a protein kinase (LinJ.33.1810-E629K).

CONCLUSION/SIGNIFICANCE

Our results suggest that differentially expressed proteins, chromosome number variations (CNVs), specific gene amplification and SNPs are important features of antimony resistance in Leishmania.

Gene amplification and point mutations in pyrimidine metabolic genes in 5-fluorouracil resistant Leishmania infantum

Background

The human protozoan parasites Leishmania are prototrophic for pyrimidines with the ability of both de novo biosynthesis and uptake of pyrimidines.

Methodology/Principal Findings

Five independent L. infantum mutants were selected for resistance to the pyrimidine analogue 5-fluorouracil (5-FU) in the hope to better understand the metabolism of pyrimidine in Leishmania. Analysis of the 5-FU mutants by comparative genomic hybridization and whole genome sequencing revealed in selected mutants the amplification of DHFR-TS and a deletion of part of chromosome 10. Point mutations in uracil phosphorybosyl transferase (UPRT), thymidine kinase (TK) and uridine phosphorylase (UP) were also observed in three individual resistant mutants. Transfection experiments confirmed that these point mutations were responsible for 5-FU resistance. Transport studies revealed that one resistant mutant was defective for uracil and 5-FU import.

Conclusion/Significance

This study provided further insights in pyrimidine metabolism in Leishmania and confirmed that multiple mutations can co-exist and lead to resistance in Leishmania.

The human protozoan parasites Leishmania present the ability of both de novo biosynthesis and uptake of pyrimidines. The pyrimidine pathway is not well understood in these parasites. In the hope to better understand the pyrimidine pathway in Leishmania, five independent L. infantum mutants were selected for resistance to the pyrimidine analogue 5-fluorouracil (5-FU). Analysis of the 5-FU mutants by comparative genomic hybridization and whole genome sequencing revealed the amplification of the main target enzyme DHFR-TS, and point mutations in three important metabolic enzymes. Transfection experiments confirmed that these point mutations were responsible for 5-FU resistance. Transport studies also revealed that one resistant mutant was defective for uracil and 5-FU import. Overall, this study provided further insights in pyrimidine metabolism in Leishmania and confirmed that multiple mutations can co-exist and lead to resistance in these protozoa.

Molecular characterization of the MRPA transporter and antimony uptake in four New World Leishmania spp. susceptible and resistant to antimony

•

Level of expression of Pgp is increased in the SbIII-resistant L. guyanensis and L. amazonensis lines.

•

Incorporation of antimony was reduced in the SbIII-resistant L. guyanensis, L. amazonensis and L. braziliensis lines.

•

Down-regulation of AQP1 protein was observed in the SbIII-resistant L. guyanensis and L. amazonensis lines.

•

Rates of SbIII efflux are higher in the SbIII-resistant lines of L. guyanensis and L. braziliensis.

•

Mechanisms of antimony-resistance of the MRPA gene are different among species of Leishmania analyzed.

ATP-binding cassette (ABC) transporters have been associated with drug resistance in various diseases. The MRPA gene, a transporter of ABCC subfamily, is involved in the resistance by sequestering metal-thiol conjugates in intracellular vesicles of Leishmania parasite. In this study, we performed the molecular characterization of the MRPA transporter, analysis of P-glycoprotein (Pgp) and aquaglyceroporin-1 (AQP1) expression, and determination of antimony level in antimony-susceptible and -resistant lines of L. (V.) guyanensis, L. (L.) amazonensis, L. (V.) braziliensis and L. (L.) infantum. PFGE analysis revealed an association of chromosomal amplification of MRPA gene with the drug resistance phenotype in all SbIII-resistant Leishmania lines analyzed. Levels of mRNA from MRPA gene determined by real-time quantitative RT-PCR showed an increased expression of two fold in SbIII-resistant lines of Leishmania guyanensis, Leishmania amazonensis and Leishmania braziliensis. Western blot analysis revealed that Pgp is increased in the SbIII-resistant L. guyanensis and L. amazonensis lines. The intracellular level of antimony quantified by graphite furnace atomic absorption spectrometry showed a reduction in the accumulation of this element in SbIII-resistant L. guyanensis, L. amazonensis and L. braziliensis lines when compared to their susceptible counterparts. Interestingly, a down-regulation of AQP1 protein was observed in the SbIII-resistant L. guyanensis and L. amazonensis lines, contributing for decreasing of SbIII entry in these lines. In addition, efflux experiments revealed that the rates of SbIII efflux are higher in the SbIII-resistant lines of L. guyanensis and L. braziliensis, that may explain also the low SbIII concentration within of these parasites. The BSO, an inhibitor of γ-glutamylcysteine synthetase enzyme, reversed the SbIII-resistance phenotype of L. braziliensis and caused an increasing in the Sb intracellular level in the LbSbR line. Our data indicate that the mechanisms of antimony-resistance are different among species of Leishmania analyzed in this study.