Role of the ABC transporter MRPA (PGPA) in antimony resistance in Leishmania infantum axenic and intracellular amastigotes

Potential selection of antimony and methotrexate cross-resistance in Leishmania infantum circulating strains

BACKGROUND

Visceral leishmaniasis (VL) resolution depends on a wide range of factors, including the instauration of an effective treatment coupled to a functional host immune system. Patients with a depressed immune system, like the ones receiving methotrexate (MTX), are at higher risk of developing VL and refusing antileishmanial drugs. Moreover, the alarmingly growing levels of antimicrobial resistance, especially in endemic areas, contribute to the increasing the burden of this complex zoonotic disease.

PRINCIPAL FINDINGS

To understand the potential links between immunosuppressants and antileishmanial drugs, we have studied the interaction of antimony (Sb) and MTX in a Leishmania infantum reference strain (LiWT) and in two L. infantum clinical strains (LiFS-A and LiFS-B) naturally circulating in non-treated VL dogs in Spain. The LiFS-A strain was isolated before Sb treatment in a case that responded positively to the treatment, while the LiFS-B strain was recovered from a dog before Sb treatment, with the dog later relapsing after the treatment. Our results show that, exposure to Sb or MTX leads to an increase in the production of reactive oxygen species (ROS) in LiWT which correlates with a sensitive phenotype against both drugs in promastigotes and intracellular amastigotes. LiFS-A was sensitive against Sb but resistant against MTX, displaying high levels of protection against ROS when exposed to MTX. LiFS-B was resistant to both drugs. Evaluation of the melting proteomes of the two LiFS, in the presence and absence of Sb and MTX, showed a differential enrichment of direct and indirect targets for both drugs, including common and unique pathways.

CONCLUSION

Our results show the potential selection of Sb-MTX cross-resistant parasites in the field, pointing to the possibility to undermine antileishmanial treatment of those patients being treated with immunosuppressant drugs in Leishmania endemic areas.

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.

Genomes of Endotrypanum monterogeii from Panama and Zelonia costaricensis from Brazil: Expansion of Multigene Families in Leishmaniinae Parasites That Are Close Relatives of Leishmania spp

The Leishmaniinae subfamily of the Trypanosomatidae contains both genus Zelonia (monoxenous) and Endotrypanum (dixenous). They are amongst the nearest known relatives of Leishmania, which comprises many human pathogens widespread in the developing world. These closely related lineages are models for the genomic biology of monoxenous and dixenous parasites. Herein, we used comparative genomics to identify the orthologous groups (OGs) shared among 26 Leishmaniinae species to investigate gene family expansion/contraction and applied two phylogenomic approaches to confirm relationships within the subfamily. The Endotrypanum monterogeii and Zelonia costaricensis genomes were assembled, with sizes of 29.9 Mb and 38.0 Mb and 9.711 and 12.201 predicted protein-coding genes, respectively. The genome of E. monterogeii displayed a higher number of multicopy cell surface protein families, including glycoprotein 63 and glycoprotein 46, compared to Leishmania spp. The genome of Z. costaricensis presents expansions of BT1 and amino acid transporters and proteins containing leucine-rich repeat domains, as well as a loss of ABC-type transporters. In total, 415 and 85 lineage-specific OGs were identified in Z. costaricensis and E. monterogeii. The evolutionary relationships within the subfamily were confirmed using the supermatrix (3384 protein-coding genes) and supertree methods. Overall, this study showed new expansions of multigene families in monoxenous and dixenous parasites of the subfamily Leishmaniinae.

Systematic Exploration of Functional Group Relevance for Anti-Leishmanial Activity of Anisomycin

Assessment of structure-activity relationships for anti-protozoan activity revealed a strategy for preparing potent anisomycin derivatives with reduced host toxicity. Thirteen anisomycin analogs were synthesized by modifying the alcohol, amine, and aromatic functional groups. Examination of anti-protozoal activity against various strains of Leishmania and cytotoxicity against leucocytes with comparison against the parent natural product demonstrated typical losses of activity with modifications of the alcohol, amine, and aromatic meta-positions. On the other hand, the para-phenol moiety of anisomycin proved an effective location for introducing substituents without significant loss of anti-protozoan potency. An entry point for differentiating activity against Leishmania versus host has been uncovered by this systematic study.

Design and Synthesis of New Anthranyl Phenylhydrazides: Antileishmanial Activity and Structure-Activity Relationship

Leishmaniasis is a neglected tropical disease affecting millions of people worldwide. A centenary approach to antimonial-based drugs was first initiated with the synthesis of urea stibamine by Upendranath Brahmachari in 1922. The need for new drug development led to resistance toward antimoniates. New drug development to treat leishmaniasis is urgently needed. In this way, searching for new substances with antileishmanial activity, we synthesized ten anthranyl phenylhydrazide and three quinazolinone derivatives and evaluated them against promastigotes and the intracellular amastigotes of Leishmania amazonensis. Three compounds showed good activity against promastigotes 1b, 1d, and 1g, with IC50 between 1 and 5 μM. These new phenylhydrazides were tested against Leishmania arginase, but they all failed to inhibit this parasite enzyme, as we have shown in a previous study. To explain the possible mechanism of action, we proposed the enzyme PTR1 as a new target for these compounds based on in silico analysis. In conclusion, the new anthranyl hydrazide derivatives can be a promising scaffold for developing new substances against the protozoa parasite.

In Vitro Drug Susceptibility of a Leishmania (Leishmania) infantum Isolate from a Visceral Leishmaniasis Pediatric Patient after Multiple Relapses

The parasitic protozoan Leishmania (Leishmania) infantum is the etiological agent of human visceral leishmaniasis in South America, an infectious disease associated with malnutrition, anemia, and hepatosplenomegaly. In Brazil alone, around 2700 cases are reported each year. Treatment failure can occur as a result of drug, host, and/or parasite-related factors. Here, we isolated a Leishmania species from a pediatric patient with visceral leishmaniasis that did not respond to chemotherapy, experiencing a total of nine therapeutic relapses and undergoing a splenectomy. The parasite was confirmed as L. (L.) infantum after sequencing of the ribosomal DNA internal transcribed spacer, and the clinical isolate, in both promastigote and amastigote forms, was submitted to in vitro susceptibility assays with all the drugs currently used in the chemotherapy of leishmaniasis. The isolate was susceptible to meglumine antimoniate, amphotericin B, pentamidine, miltefosine, and paromomycin, similarly to another strain of this species that had previously been characterized. These findings indicate that the multiples relapses observed in this pediatric patient were not due to a decrease in the drug susceptibility of this isolate; therefore, immunophysiological aspects of the patient should be further investigated to understand the basis of treatment failure in this case.

Label-Free Mass Spectrometry Proteomics Reveals Different Pathways Modulated in THP-1 Cells Infected with Therapeutic Failure and Drug Resistance Leishmania infantum Clinical Isolates

As the world is facing increasing difficulties to treat leishmaniasis with current therapies, deeper investigation into the molecular mechanisms responsible for both drug resistance and treatment failure (TF) is essential in drug discovery and development. So far, few available drugs cause severe side effects and have developed several resistance mechanisms. Drug resistance and TF parasite strains from clinical isolates may have acquired altered expression of proteins that characterize specific mechanisms leading to therapy inefficacy. This work aims to identify the biochemical pathways of THP-1 human monocytes infected by different Leishmania infantum clinical isolates from patients with either resistance or with TF outcome, using whole cell differential Mass Spectrometry proteomics. We have adopted network enrichment analysis to integrate the transcriptomics and the proteomic results of infected cells studies. Transferrin receptor C (TFRC) and nucleoside diphosphate kinase 3 (NDK3) were discovered as overexpressed proteins in THP-1 cells infected with paromomycin, antimony, and miltefosine resistant L. infantum lines. The overall achievements represent founding concepts to confirm new targets involved in the parasitic drug resistance and TF mechanisms, and to consider in perspective the importance of a dual host-guest pharmacological approach to treat the acute stage of the disease.

Identification of Leishmania infantum Epidemiology, Drug Resistance and Pathogenicity Biomarkers with Nanopore Sequencing

The emergence of drug-resistant strains of the parasite Leishmania infantum infecting dogs and humans represents an increasing threat. L. infantum genomes are complex and unstable with extensive structural variations, ranging from aneuploidies to multiple copy number variations (CNVs). These CNVs have recently been validated as biomarkers of Leishmania concerning virulence, tissue tropism, and drug resistance. As a proof-of-concept to develop a novel diagnosis platform (LeishGenApp), four L. infantum samples from humans and dogs were nanopore sequenced. Samples were epidemiologically typed within the Mediterranean L. infantum group, identifying members of the JCP5 and non-JCP5 subgroups, using the conserved region (CR) of the maxicircle kinetoplast. Aneuploidies were frequent and heterogenous between samples, yet only chromosome 31 tetrasomy was common between all the samples. A high frequency of aneuploidies was observed for samples with long passage history (MHOM/TN/80/IPT-1), whereas fewer were detected for samples maintained in vivo (MCRI/ES/2006/CATB033). Twenty-two genes were studied to generate a genetic pharmacoresistance profile against miltefosine, allopurinol, trivalent antimonials, amphotericin, and paromomycin. MHOM/TN/80/IPT-1 and MCRI/ES/2006/CATB033 displayed a genetic profile with potential resistance against miltefosine and allopurinol. Meanwhile, MHOM/ES/2016/CATB101 and LCAN/ES/2020/CATB102 were identified as potentially resistant against paromomycin. All four samples displayed a genetic profile for resistance against trivalent antimonials. Overall, this proof-of-concept revealed the potential of nanopore sequencing and LeishGenApp for the determination of epidemiological, drug resistance, and pathogenicity biomarkers in L. infantum.

Antimony resistance mechanism in genetically different clinical isolates of Indian Kala-azar patients

The central theme of this enterprise is to find common features, if any, displayed by genetically different antimony (Sb)-resistant viscerotropic Leishmania parasites to impart Sb resistance. In a limited number of clinical isolates (n = 3), we studied the breadth of variation in the following dimensions: (a) intracellular thiol content, (b) cell surface expression of glycan having N-acetyl-D-galactosaminyl residue as the terminal sugar, and (c) gene expression of thiol-synthesizing enzymes (CBS, MST, gamma-GCS, ODC, and TR), antimony-reducing enzymes (TDR and ACR2), and antimonial transporter genes (AQP1, MRPA, and PRP1). One of the isolates, T5, that was genotypically characterized as Leishmania tropica, caused Indian Kala-azar and was phenotypically Sb resistant (T5-LT-SSG-R), while the other two were Leishmania donovani, out of which one isolate, AG83, is antimony sensitive (AG83-LD-SSG-S) and the other isolate, T8, is Sb resistant (T8-LD-SSG-R). Our study showed that the Sb-resistant parasites, regardless of their genotype, showed significantly higher intracellular thiol compared with Sb-sensitive AG83-LD-SSG-S. Seemingly, T5-LT-SSG-R showed about 1.9-fold higher thiol content compared with T8-LD-SSG-R which essentially mirrored cell surface N-acetyl-D-galactosaminyl expression. Except TR, the expression of the remaining thiol-synthesizing genes was significantly higher in T8-LD-SSG-R and T5-LT-SSG-R than the sensitive one, and between the Sb-resistant parasites, the latter showed a significantly higher expression. Furthermore, the genes for Sb-reducing enzymes increased significantly in resistant parasites regardless of genotype compared with the sensitive one, and between two resistant parasites, there was hardly any difference in expression. Out of three antimony transporters, AQP1 was decreased with the concurrent increase in MRPA and PRP1 in resistant isolates when compared with the sensitive counterpart. Interestingly, no difference in expression of the above-mentioned transporters was noted between two Sb-resistant isolates. The enduring image that resonated from our study is that the genetically diverse Sb-resistant parasites showed enhanced thiol-synthesizing and antimony transporter gene expression than the sensitive counterpart to confer a resistant phenotype.

Global distribution of treatment resistance gene markers for leishmaniasis

BACKGROUND

Pentavalent antimonials (Sb(V)) such as meglumine antimoniate (Glucantime®) and sodium stibogluconate (Pentostam®) are used as first-line treatments for leishmaniasis, either alone or in combination with second-line drugs such as amphotericin B (Amp B), miltefosine (MIL), methotrexate (MTX), or cryotherapy. Therapeutic aspects of these drugs are now challenged because of clinical resistance worldwide.

METHODS

We reviewedthe recent original studies were assessed by searching in electronic databases such as Scopus, Pubmed, Embase, and Web of Science.

RESULTS

Studies on molecular biomarkers involved in drug resistance are essential for monitoring the disease. We reviewed genes and mechanisms of resistance to leishmaniasis, and the geographical distribution of these biomarkers in each country has also been thoroughly investigated.

CONCLUSION

Due to the emergence of resistant genes mainly in anthroponotic Leishmania species such as L. donovani and L. tropica, as the causative agents of ACL and AVL, respectively, selection of an appropriate treatment modality is essential. Physicians should be aware of the presence of such resistance for the selection of proper treatment modalities in endemic countries.

Leishmania parasites exchange drug-resistance genes through extracellular vesicles

Leishmania are eukaryotic parasites that have retained the ability to produce extracellular vesicles (EVs) through evolution. To date, it has been unclear if different DNA entities could be associated with Leishmania EVs and whether these could constitute a mechanism of horizontal gene transfer (HGT). Herein, we investigate the DNA content of EVs derived from drug-resistant parasites, as well as the EVs' potential to act as shuttles for DNA transfer. Next-generation sequencing and PCR assays confirm the enrichment of amplicons carrying drug-resistance genes associated with EVs. Transfer assays of drug-resistant EVs highlight a significant impact on the phenotype of recipient parasites induced by the expression of the transferred DNA. Recipient parasites display an enhanced growth and better control of oxidative stress. We provide evidence that eukaryotic EVs function as efficient mediators in HGT, thereby facilitating the transmission of drug-resistance genes and increasing the fitness of cells when encountering stressful environments.

The role of ATP-binding cassette transporter genes expression in treatment failure cutaneous leishmaniasis

Leishmaniasis is one of the common diseases transmitted by sand flies in tropical and subtropical regions of the world. Currently, antimonial derivatives are the first line of treatment. Some of the members of the ATP-binding cassette (ABC) family of Leishmania are shown to be associated with no response to treatment. In this study, we evaluated ABCI4, ABCG2, ABCC7, ABCB4, and ABCC3 genes expression in Leishmania isolated from patients with non-healing cutaneous leishmaniasis and treatment response isolates. We selected 17 clinical isolates including 8 treatment failure and 9 treatment response samples from September 2020 to March 2021. The isolates were obtained from patients of Health Center Laboratory of Varzaneh, Isfahan, Iran with cutaneous leishmaniasis. The diagnosis was performed using microscopic observation. The samples were directly collected from the lesions. The expression profiling of genes was assessed using SYBR Green real-time PCR that was analyzed with delta-delta Ct. All treatment failure clinical isolates were L. major. Gene expression analysis in treatment failure isolates showed that the ABC transported genes had a different pattern in each isolate. Treatment failure has been reported for cutaneous leishmaniasis worldwide. Knowledge of the molecular mechanisms of treatment failure could solve this problem. ABC transporter genes are considered controversial over the mechanisms of treatment failure outcomes. In this study, we showed that ABC transporter genes could be considered one of the important mechanisms.

Tackling Drug Resistance and Other Causes of Treatment Failure in Leishmaniasis

Leishmaniasis is a tropical infectious disease caused by the protozoan Leishmania parasite. The disease is transmitted by female sand flies and, depending on the infecting parasite species, causes either cutaneous (stigmatizing skin lesions), mucocutaneous (destruction of mucous membranes of nose, mouth and throat) or visceral disease (a potentially fatal infection of liver, spleen and bone marrow). Although more than 1 million new cases occur annually, chemotherapeutic options are limited and their efficacy is jeopardized by increasing treatment failure rates and growing drug resistance. To delay the emergence of resistance to existing and new drugs, elucidating the currently unknown causes of variable drug efficacy (related to parasite susceptibility, host immunity and drug pharmacokinetics) and improved use of genotypic and phenotypic tools to define, measure and monitor resistance in the field are critical. This review highlights recent progress in our understanding of drug action and resistance in Leishmania, ongoing challenges (including setbacks related to the COVID-19 pandemic) and provides an overview of possible strategies to tackle this public health challenge.

Three different mutations in the DNA topoisomerase 1B in Leishmania infantum contribute to resistance to antitumor drug topotecan

Background

The evolution of drug resistance is one of the biggest challenges in leishmaniasis and has prompted the need for new antileishmanial drugs. Repurposing of approved drugs is a faster and very attractive strategy that is gaining supporters worldwide. Different anticancer topoisomerase 1B (TOP1B) inhibitors have shown strong antileishmanial activity and promising selective indices, supporting the potential repurposing of these drugs. However, cancer cells and Leishmania share the ability to become rapidly resistant. The aim of this study was to complete a whole-genome exploration of the effects caused by exposure to topotecan in order to highlight the potential mechanisms deployed by Leishmania to favor its survival in the presence of a TOP1B inhibitor.

Methods

We used a combination of stepwise drug resistance selection, whole-genome sequencing, functional validation, and theoretical approaches to explore the propensity of and potential mechanisms deployed by three independent clones of L. infantum to resist the action of TOP1B inhibitor topotecan.

Results

We demonstrated that L. infantum is capable of becoming resistant to high concentrations of topotecan without impaired growth ability. No gene deletions or amplifications were identified from the next-generation sequencing data in any of the three resistant lines, ruling out the overexpression of efflux pumps as the preferred mechanism of topotecan resistance. We identified three different mutations in the large subunit of the leishmanial TOP1B (Top1B^F187Y, Top1B^G191A, and Top1B^W232R). Overexpression of these mutated alleles in the wild-type background led to high levels of resistance to topotecan. Computational molecular dynamics simulations, in both covalent and non-covalent complexes, showed that these mutations have an effect on the arrangement of the catalytic pentad and on the interaction of these residues with surrounding amino acids and DNA. This altered architecture of the binding pocket results in decreased persistence of topotecan in the ternary complex.

Conclusions

This work helps elucidate the previously unclear potential mechanisms of topotecan resistance in Leishmania by mutations in the large subunit of TOP1B and provides a valuable clue for the design of improved inhibitors to combat resistance in both leishmaniasis and cancer. Our data highlights the importance of including drug resistance evaluation in drug discovery cascades.

Graphical abstract

Supplementary Information

The online version contains supplementary material available at 10.1186/s13071-021-04947-4.

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.

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.

Identification of Resistance Determinants for a Promising Antileishmanial Oxaborole Series

Current treatment options for visceral leishmaniasis have several drawbacks, and clinicians are confronted with an increasing number of treatment failures. To overcome this, the Drugs for Neglected Diseases initiative (DNDi) has invested in the development of novel antileishmanial leads, including a very promising class of oxaboroles. The mode of action/resistance of this series to Leishmania is still unknown and may be important for its further development and implementation. Repeated in vivo drug exposure and an in vitro selection procedure on both extracellular promastigote and intracellular amastigote stages were both unable to select for resistance. The use of specific inhibitors for ABC-transporters could not demonstrate the putative involvement of efflux pumps. Selection experiments and inhibitor studies, therefore, suggest that resistance to oxaboroles may not emerge readily in the field. The selection of a genome-wide cosmid library coupled to next-generation sequencing (Cos-seq) was used to identify resistance determinants and putative targets. This resulted in the identification of a highly enriched cosmid, harboring genes of chromosome 2 that confer a subtly increased resistance to the oxaboroles tested. Moderately enriched cosmids encompassing a region of chromosome 34 contained the cleavage and polyadenylation specificity factor (cpsf) gene, encoding the molecular target of several related benzoxaboroles in other organisms.

Influence of N-Methylation and Conformation on Almiramide Anti-Leishmanial Activity

The almiramide N-methylated lipopeptides exhibit promising activity against trypanosomatid parasites. A structure–activity relationship study has been performed to examine the influences of N-methylation and conformation on activity against various strains of leishmaniasis protozoan and on cytotoxicity. The synthesis and biological analysis of twenty-five analogs demonstrated that derivatives with a single methyl group on either the first or fifth residue amide nitrogen exhibited greater activity than the permethylated peptides and relatively high potency against resistant strains. Replacement of amino amide residues in the peptide, by turn inducing α amino γ lactam (Agl) and N-aminoimidazalone (Nai) counterparts, reduced typically anti-parasitic activity; however, peptide amides possessing Agl residues at the second residue retained significant potency in the unmethylated and permethylated series. Systematic study of the effects of methylation and turn geometry on anti-parasitic activity indicated the relevance of an extended conformer about the central residues, and conformational mobility by tertiary amide isomerization and turn geometry at the extremities of the active peptides.

The ultimate fate determinants of drug induced cell-death mechanisms in Trypanosomatids

Chemotherapy constitutes a major part of modern-day therapy for infectious and chronic diseases. A drug is said to be effective if it can inhibit its target, induce stress, and thereby trigger an array of cell death pathways in the form of programmed cell death, autophagy, necrosis, etc. Chemotherapy is the only treatment choice against trypanosomatid diseases like Leishmaniasis, Chagas disease, and sleeping sickness. Anti-trypanosomatid drugs can induce various cell death phenotypes depending upon the drug dose and growth stage of the parasites. The mechanisms and pathways triggering cell death in Trypanosomatids serve to help identify potential targets for the development of effective anti-trypanosomatids. Studies show that the key proteins involved in cell death of trypanosomatids are metacaspases, Endonuclease G, Apoptosis-Inducing Factor, cysteine proteases, serine proteases, antioxidant systems, etc. Unlike higher eukaryotes, these organisms either lack the complete set of effectors involved in cell death pathways, or are yet to be deciphered. A detailed summary of the existing knowledge of different drug-induced cell death pathways would help identify the lacuna in each of these pathways and therefore open new avenues for research and thereby new therapeutic targets to explore. The cell death pathway associated complexities in metazoans are absent in trypanosomatids; hence this summary can also help understand the trigger points as well as cross-talk between these pathways. Here we provide an in-depth overview of the existing knowledge of these drug-induced trypanosomatid cell death pathways, describe their associated physiological changes, and suggest potential interconnections amongst them.

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.

Unravelling the proteomic signature of extracellular vesicles released by drug-resistant Leishmania infantum parasites

Leishmaniasis constitutes the 9^th largest disease burden among all infectious diseases. Control of this disease is based on a short list of chemotherapeutic agents headed by pentavalent antimonials, followed by miltefosine and amphotericin B; drugs that are far from ideal due to host toxicity, elevated cost, limited access, and high rates of drug resistance. Knowing that the composition of extracellular vesicles (EVs) can vary according to the state of their parental cell, we hypothesized that EVs released by drug-resistant Leishmania infantum parasites could contain unique and differently enriched proteins depending on the drug-resistance mechanisms involved in the survival of their parental cell line. To assess this possibility, we studied EV production, size, morphology, and protein content of three well-characterized drug-resistant L. infantum cell lines and a wild-type strain. Our results are the first to demonstrate that drug-resistance mechanisms can induce changes in the morphology, size, and distribution of L. infantum EVs. In addition, we identified L. infantum’s core EV proteome. This proteome is highly conserved among strains, with the exception of a handful of proteins that are enriched differently depending on the drug responsible for induction of antimicrobial resistance. Furthermore, we obtained the first snapshot of proteins enriched in EVs released by antimony-, miltefosine- and amphotericin-resistant parasites. These include several virulence factors, transcription factors, as well as proteins encoded by drug-resistance genes. This detailed study of L. infantum EVs sheds new light on the potential roles of EVs in Leishmania biology, particularly with respect to the parasite’s survival in stressful conditions. This work outlines a crucial first step towards the discovery of EV-based profiles capable of predicting response to antileishmanial agents.

Visceral leishmaniasis is a life-threatening disease caused by Leishmania infantum parasites, which are transmitted by sand flies. In the absence of vaccines, current control of this disease is based on chemotherapy, which is comprised of a very limited arsenal threatened by the emergence and spread of drug-resistant strains. In the shadow of growing concern and treatment failure due to resistance, the characterization of extracellular vesicles (EVs) released by drug-resistant L. infantum parasites could shed some light on the complex nature of drug resistance in Leishmania and increase our understanding of the biology of the parasite. EVs are vesicles secreted by all eukaryotic cells whose contents (proteins, DNA/RNAs, lipids) vary as a function of their cellular origin. Our results demonstrate for the first time that EVs released by drug-resistant parasites are enriched in unique protein markers that reflect the drug-resistance mechanisms involved in the survival of parental cells. These unique proteins included several virulence and transcription factors, as well as drug-resistance genes; this offers a potential benefit for drug-resistant parasites in terms of parasite-to-parasite communication and host-parasite interactions. Collectively, our initial results could serve as a jumping-off point for the future development of novel EV-based diagnostic tools for the detection and appraisal of antimicrobial-resistant Leishmania populations.

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.

The Novel Serine/Threonine Protein Kinase LmjF.22.0810 from Leishmania major may be Involved in the Resistance to Drugs such as Paromomycin

The identification and clarification of the mechanisms of action of drugs used against leishmaniasis may improve their administration regimens and prevent the development of resistant strains. Herein, for the first time, we describe the structure of the putatively essential Ser/Thr kinase LmjF.22.0810 from Leishmania major. Molecular dynamics simulations were performed to assess the stability of the kinase model. The analysis of its sequence and structure revealed two druggable sites on the protein. Furthermore, in silico docking of small molecules showed that aminoglycosides preferentially bind to the phosphorylation site of the protein. Given that transgenic LmjF.22.0810-overexpressing parasites displayed less sensitivity to aminoglycosides such as paromomycin, our predicted models support the idea that the mechanism of drug resistance observed in those transgenic parasites is the tight binding of such compounds to LmjF.22.0810 associated with its overexpression. These results may be helpful to understand the complex machinery of drug response in Leishmania.

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.

Comparative Two-dimensional Gel Electrophoresis Maps for Amastigote-like Proteomes of Iranian Leishmania Tropica and Leishmania Major Isolates

BACKGROUND

Leishmania major and Leishmania tropica are the main causative agents of cutaneous leishmaniasis. Proteomics as a novel approaches could be used to evaluate protein expression levels in different stages of Leishmania species. We compare the protein contents of amastigote-like forms in L. tropica and L. major using two-dimensional gel electrophoresis (2-DE) and bioinformatics methods.

MATERIALS AND METHODS

Leishmania parasites were isolated from the lesions of Iranian patients and identified using restriction fragment length polymorphism-polymerase chain reaction (RFLP-PCR). Five isolates of each two species were cultured in specific media to obtain amastigote-like forms to be prepared for proteomics study. Total protein contents were separated using 2-DE. The gels were stained by silver nitrate and scan was imaged. The protein spots with different expression changes in each gel were analyzed using Progenesis SameSpots software.

RESULTS

A total of 354 protein spots were detected in both amastigote-like forms. Comparative analysis of protein spots with different expressions in the two amastigote-like form species showed 173 highly expressed spots of which 74 L. tropica and 99 L. major proteins were spotted with fold≥2. Also, 16 and 20 new protein spots were uniquely found in L. tropica and L. major, respectively. Clustering of different detected proteins using correlation analysis divided the proteins into two clusters based on their expression level. Furthermore, clustering results were confirmed by principal component analysis.

CONCLUSION

Using proteomics methods specially 2-DE and statistical analysis demonstrated significant changes in protein expression levels in amastigote-like forms of L. tropica and L. major isolates.

Chemogenomic Profiling of Antileishmanial Efficacy and Resistance in the Related Kinetoplastid Parasite Trypanosoma brucei

The arsenal of drugs used to treat leishmaniasis, caused by Leishmania spp., is limited and beset by toxicity and emergent resistance. Furthermore, our understanding of drug mode of action and potential routes to resistance is limited. Forward genetic approaches have revolutionized our understanding of drug mode of action in the related kinetoplastid parasite Trypanosoma brucei.

The arsenal of drugs used to treat leishmaniasis, caused by Leishmania spp., is limited and beset by toxicity and emergent resistance. Furthermore, our understanding of drug mode of action and potential routes to resistance is limited. Forward genetic approaches have revolutionized our understanding of drug mode of action in the related kinetoplastid parasite Trypanosoma brucei. Therefore, we screened our genome-scale T. brucei RNA interference (RNAi) library against the current antileishmanial drugs sodium stibogluconate (antimonial), paromomycin, miltefosine, and amphotericin B. Identification of T. brucei orthologues of the known Leishmania antimonial and miltefosine plasma membrane transporters effectively validated our approach, while a cohort of 42 novel drug efficacy determinants provides new insights and serves as a resource. Follow-up analyses revealed the antimonial selectivity of the aquaglyceroporin TbAQP3. A lysosomal major facilitator superfamily transporter contributes to paromomycin-aminoglycoside efficacy. The vesicle-associated membrane protein TbVAMP7B and a flippase contribute to amphotericin B and miltefosine action and are potential cross-resistance determinants. Finally, multiple phospholipid-transporting flippases, including the T. brucei orthologue of the Leishmania miltefosine transporter, a putative β-subunit/CDC50 cofactor, and additional membrane-associated hits, affect amphotericin B efficacy, providing new insights into mechanisms of drug uptake and action. The findings from this orthology-based chemogenomic profiling approach substantially advance our understanding of antileishmanial drug action and potential resistance mechanisms and should facilitate the development of improved therapies as well as surveillance for drug-resistant parasites.

Validation of NAD synthase inhibitors for inhibiting the cell viability of Leishmania donovani: In silico and in vitro approach

NAD (nicotinamide adenine dinucleotide) synthase catalyses the biochemical synthesis of NAD, from nicotinic acid adenine dinucleotide (NAAD). NAD may be synthesized through the de novo pathways and/or the salvage pathways in cells. However, in Leishmania parasite, the synthesis of NAD solely depends on the salvage pathways. NAD synthetase is widely explored as a drug target in various microorganisms. In Bacillus anthracis, a group of sulphonamides 5599, 5617 and 5824 and complex amide 5833 were reported to have activity at micromolar range against NAD synthetase. Hence, in the present study, the same group of sulphonamides and complex amide were validated through in silico and in vitro studies for its efficiency towards Leishmania donovani NAD synthase. In silico study revealed the ligands 5824 and 5833 to have better docking score. Molecular dynamics simulation for a duration of 50 ns of all the ligand-protein complexes suggested that the complexes with the ligands 5824 and 5833 were stable and interacting. In vitro and ex vivo studies have shown that 5824 and 5833 inhibit the cell viability of the organism at a lower concentration than 5599 and 5617. Hence, with further in vivo validation, 5824 (or its synthetic analogues) and 5833 could be the choice that may work synergistically with other potential drugs in treating drug-resistant cases of leishmaniasis. Communicated by Ramaswamy H. Sarma.

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.

Gene expression analysis of antimony resistance in Leishmania tropica using quantitative real-time PCR focused on genes involved in trypanothione metabolism and drug transport

Pentavalent antimonials remain the treatment of choice for all the clinical forms of leishmaniasis. The increasing rates of antimony resistance are becoming a serious health problem in treatment of anthroponotic cutaneous leishmaniasis (ACL). Accordingly, unraveling molecular markers is crucial for improving medication strategies and monitoring of drug-resistant parasites. Different studies have suggested the importance of genes involved in trypanothione metabolism and drug transport. In this regard, present study was designed to investigate the RNA expression level of five genes including γ-GCS, ODC, TRYR (involved in trypanothione metabolism), AQP1 (acts in drug uptake) and MRPA (involved in sequestration of drug) in sensitive and resistant Leishmania tropica isolates. Seven antimony-resistant and seven antimony-sensitive L. tropica clinical isolates were collected from ACL patients. Drug sensitivity test was performed on the samples as well as reference strains; afterwards, gene expression analysis was performed on clinical isolates by quantitative real-time PCR. The results revealed that the average expression level of AQP1 gene was decreased (0.47-fold) in resistant isolates compared to sensitive ones whereas MRPA (2.45), γ-GCS (2.1) and TRYR (1.97) was upregulated in resistant isolates. The average expression of ODC (1.24-fold) gene was not different significantly between sensitive and resistant isolates. Our findings suggest that AQP1, MRPA, GSH1 and TRYR can be considered as potential molecular markers for screening of antimony resistance in some L. tropica clinical isolates.

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.

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.

High-throughput Cos-Seq screen with intracellular Leishmania infantum for the discovery of novel drug-resistance mechanisms

Increasing drug resistance towards first line antimony-derived compounds has forced the introduction of novel therapies in leishmaniasis endemic areas including amphotericin B and miltefosine. However, their use is threatened by the emergence and spread of drug-resistant strains. In order to discover stage-dependent resistance genes, we have adapted the Cos-Seq approach through the introduction of macrophage infections in the pipeline. A L. infantum intracellular amastigote population complemented with a L. infantum cosmid library was submitted to increasing concentrations of miltefosine, amphotericin B and pentavalent antimonials in experimental infections of THP-1 cells. For each step of selection, amastigotes were extracted and cosmids were isolated and submitted to next-generation sequencing, followed by subsequent gene-enrichment analyses. Cos-Seq screen in amastigotes revealed four highly enriched loci for antimony, five for miltefosine and one for amphotericin B. Of these, a total of seven cosmids were recovered and tested for resistance in both promastigotes and amastigotes. Candidate genes within the pinpointed genomic regions were validated using single gene overexpression in wild-type parasites and/or gene disruption by means of a CRISPR-Cas9-based approach. This led to the identification and validation of a stage-independent antimony-resistance gene (LinJ.06.1010) coding for a putative leucine rich repeat protein and a novel amastigote-specific miltefosine-resistance gene (LinJ.32.0050) coding for a member of the SEC13 family of WD-repeat proteins. This study further reinforces the power of Cos-Seq approach to discover novel drug-resistance genes, some of which are life-stages specific.

•

The Cos-Seq led to the discovery of several new genomic regions selected with drugs.

•

This work led to the validation of novel drug-resistance genes in Leishmania.

•

Gene LinJ.06.1010 is involved in antimony resistance in both life stages of the parasite.

•

Gene LinJ.32.0050 is involved in miltefosine resistance in amastigotes.

•

The amastigote screen is labour intensive but complements screens in promastigotes.

Slow pace of antileishmanial drug development

The protozoan parasiteLeishmaniais endemic in large parts of the world which causes leishmaniasis. Its visceral form is fatal if not treated and is caused mostly byLeishmania donovani,Leishmania infantumandLeishmania chagasi. Given the difficulties linked to vector (sandfly) control and the lack of an effective vaccine, the control of leishmaniasis relies mostly on chemotherapy. Unfortunately, the prevalence of parasites becoming resistant to the first-line drug pentavalent antimony (SbV) is increasing worldwide. Few alternative drugs are available that includes amphotericin B, pentamidine and miltefosine (oral). Already, decreases in efficacy, resistance and toxicity have been noted against these drugs. Dry antileishmanial pipeline further indicates the slow pace of drug discovery in this field where resistance as a major barrier. Full understanding of the genetic and molecular basis of the parasite is lagging. Since leishmaniasis is a neglected disease and occurs predominantly in the developing world largely, therefore, it is unaddressed. The pharma industry argues that development of the new drug is too costly and risky to invest in low return neglected diseases is very high. Research is also needed to identify new and effective drug targets. The lack of drug research and development for neglected diseases will require some new strategies. We have discussed the various cause of slow pace of antileishmanial drug discovery in this review to pay attention of researchers and also take the public and private initiative to make the process fast for new antileishmanial drug development.

Genome wide comparison of Ethiopian Leishmania donovani strains reveals differences potentially related to parasite survival

Leishmania donovani is the main cause of visceral leishmaniasis (VL) in East Africa. Differences between northern Ethiopia/Sudan (NE) and southern Ethiopia (SE) in ecology, vectors, and patient sensitivity to drug treatment have been described, however the relationship between differences in parasite genotype between these two foci and phenotype is unknown. Whole genomic sequencing (WGS) was carried out for 41 L. donovani strains and clones from VL and VL/HIV co-infected patients in NE (n = 28) and SE (n = 13). Chromosome aneuploidy was observed in all parasites examined with each isolate exhibiting a unique karyotype. Differences in chromosome ploidy or karyotype were not correlated with the geographic origin of the parasites. However, correlation between single nucleotide polymorphism (SNP) and geographic origin was seen for 38/41 isolates, separating the NE and SE parasites into two large groups. SNP restricted to NE and SE groups were associated with genes involved in viability and parasite resistance to drugs. Unique copy number variation (CNV) were also associated with NE and SE parasites, respectively. One striking example is the folate transporter (FT) family genes (LdBPK_100390, LdBPK_100400 and LdBPK_100410) on chromosome 10 that are single copy in all 13 SE isolates, but either double copy or higher in 39/41 NE isolates (copy number 2-4). High copy number (= 4) was also found for one Sudanese strain examined. This was confirmed by quantitative polymerase chain reaction for LdBPK_100400, the L. donovani FT1 transporter homolog. Good correlation (p = 0.005) between FT copy number and resistance to methotrexate (0.5 mg/ml MTX) was also observed with the haploid SE strains examined showing higher viability than the NE strains at this concentration. Our results emphasize the advantages of whole genome analysis to shed light on vital parasite processes in Leishmania.

Drug resistance and treatment failure in leishmaniasis: A 21st century challenge

Reevaluation of treatment guidelines for Old and New World leishmaniasis is urgently needed on a global basis because treatment failure is an increasing problem. Drug resistance is a fundamental determinant of treatment failure, although other factors also contribute to this phenomenon, including the global HIV/AIDS epidemic with its accompanying impact on the immune system. Pentavalent antimonials have been used successfully worldwide for the treatment of leishmaniasis since the first half of the 20th century, but the last 10 to 20 years have witnessed an increase in clinical resistance, e.g., in North Bihar in India. In this review, we discuss the meaning of “resistance” related to leishmaniasis and discuss its molecular epidemiology, particularly for Leishmania donovani that causes visceral leishmaniasis. We also discuss how resistance can affect drug combination therapies. Molecular mechanisms known to contribute to resistance to antimonials, amphotericin B, and miltefosine are also outlined.

Chemotherapy is central to the control and management of leishmaniasis. Antimonials remain the primary drugs against different forms of leishmaniasis in several regions. However, resistance to antimony has necessitated the use of alternative medications, especially in the Indian subcontinent (ISC). Compounds, notably the orally available miltefosine (MIL), parenteral paromomycin, and amphotericin B (AmB), are increasingly used to treat leishmaniasis. Although treatment failure (TF) has been observed in patients treated with most anti-leishmanials, its frequency of appearance may be important in patients treated with MIL, which has replaced antimonials within the kala-azar elimination program in the ISC. AmB is highly efficacious, and the associated toxic effects—when administered in its free deoxycholate form—are somewhat ameliorated in its liposomal formulation. Regrettably, laboratory experimentation has demonstrated a risk of resistance towards AmB as well. The rise of drug resistance impacts treatment outcome, and understanding its causes, spread, and impact will help us manage the risks it imposes. Here, we review the problem of TF in leishmaniasis and the contribution of drug resistance to the problem. Molecular mechanisms causing resistance to anti-leishmanials are discussed along with the appropriate use of additional available drugs, as well as the urgent need to consolidate strategies to monitor drug efficacy, epidemiological surveillance, and local policies. Coordination of these activities in national and international programs against leishmaniasis might represent a successful guide to further research and prevention activities.

Transcriptome profiling identifies genes/pathways associated with experimental resistance to paromomycin in Leishmania donovani

Widespread resistance towards antimony and reports of relapses following miltefosine treatment has severely affected the management of visceral leishmaniasis (VL) in the Indian subcontinent. Paromomycin (PMM), an aminoglycoside antibiotic, has been licensed for VL treatment in India in 2007. Although its use is still restricted in the field, unraveling the molecular mechanism of resistance towards PMM is the key to preserve the drug. In this study, PMM resistant lines were selected up to 100 μM of PMM in three distinct field isolates of Leishmania donovani at promastigote stage. The resistance induced at promastigote level was also evident in amastigotes which showed 6 fold decreases in PMM susceptibility. Comparative transcriptome profiling of PMM resistant (PMM-R) and the corresponding PMM sensitive (PMM-S) parasites revealed modulated expression of 500 genes (1.5 fold cut off) in PMM-R parasites. Selected genes were validated for their modulated expression by quantitative real-time PCR. Functional classification and pathway analysis of modulated genes indicated probable adaptations in drug resistant lines which included a) reduced oxidative phosphorylation; b) increased glycosomal succinate fermentation and substrate level phosphorylation; c) dependency on lipids and amino acids for energy generation; d) reduced DNA synthesis and increased DNA damage repair and e) decreased protein synthesis and degradation. Interestingly, PMM-R parasites showed a marked increase in PMM susceptibility in presence of verapamil and amlodipine, antagonists of Ca2+ channel that are also modulators of ABC transporters. Moreover, infection of macrophages by PMM-R parasites led to modulated nitric oxide (NO) levels while reactive oxygen species (ROS) level remained unaltered. The present study highlights the putative mechanisms of PMM resistance in Leishmania.

Miltefosine Resistant Field Isolate From Indian Kala-Azar Patient Shows Similar Phenotype in Experimental Infection

Emergence of resistance to drugs used to treat the Indian Kala-azar patients makes control strategy shattered. In this bleak situation, Miltefosine (MIL) was introduced to treat mainly antimonial unresponsive cases. Within years, resistance to MIL has been reported. While checking the MIL sensitivity of the recent KA clinical isolates (n = 26), we came across one isolate which showed four times more EC₅₀ for MIL than that of MIL-Sensitive (MIL-S) isolates and considered as putative MIL-Resistant (MIL-R). The expressions of LdMT and LdRos3 genes of this isolate were found down regulated. Th1/Th2 cytokines, ROS and NO, FACS dot plots and mitochondrial trans membrane potential measurement were performed. In vivo hamster model with this MIL-R isolate showed much lesser reduction in liver weight (17.5%) compared to average reduction in liver weight (40.2%) of the animals infected with MIL-S isolates. The splenic and hepatic stamps smears of MIL-R infected hamsters revealed the retention of parasite load of about 51.45%. The splenocytes of these animals failed to proliferate anti leishmanial T-cells and lack of cell mediated immunity hampered recovery. Thus, these phenotypic expressions of experimental model may be considered similar to that of the MIL unresponsive patients. This is first such kind of report.

Leishmania LABCG1 and LABCG2 transporters are involved in virulence and oxidative stress: functional linkage with autophagy

Background

The G subfamily of ABC (ATP-binding cassette) transporters of Leishmania include 6 genes (ABCG1-G6), some with relevant biological functions associated with drug resistance and phospholipid transport. Several studies have shown that Leishmania LABCG2 transporter plays a role in the exposure of phosphatidylserine (PS), in virulence and in resistance to antimonials. However, the involvement of this transporter in other key biological processes has not been studied.

Methods

To better understand the biological function of LABCG2 and its nearly identical tandem-repeated transporter LABCG1, we have generated Leishmania major null mutant parasites for both genes (ΔLABCG1-2). NBD-PS uptake, infectivity, metacyclogenesis, autophagy and thiols were measured.

Results

Leishmania major ΔLABCG1-2 parasites present a reduction in NBD-PS uptake, infectivity and virulence. In addition, we have shown that ΔLABCG1-2 parasites in stationary phase growth underwent less metacyclogenesis and presented differences in the plasma membrane’s lipophosphoglycan composition. Considering that autophagy is an important process in terms of parasite virulence and cell differentiation, we have shown an autophagy defect in ΔLABCG1-2 parasites, detected by monitoring expression of the autophagosome marker RFP-ATG8. This defect correlates with increased levels of reactive oxygen species and higher non-protein thiol content in ΔLABCG1-2 parasites. HPLC analysis revealed that trypanothione and glutathione were the main molecules accumulated in these ΔLABCG1-2 parasites. The decrease in non-protein thiol levels due to preincubation with buthionine sulphoximide (a γ-glutamylcysteine synthetase inhibitor) restored the autophagy process in ΔLABCG1-2 parasites, indicating a relationship between autophagy and thiol content.

Conclusions

LABCG1-2 transporters from Leishmania could be considered as phosphatidylserine and non-protein thiol transporters. They probably accomplish transportation in conjunction with other molecules that are involved in oxidative stress, autophagy, metacyclogenesis and infectivity processes. The overall conclusion is that LABCG1-2 transporters could play a key role in Leishmania cell survival and infectivity.

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.

Involvement of nucleoside diphosphate kinase b and elongation factor 2 in Leishmania braziliensis antimony resistance phenotype

Background

Nucleoside diphosphate kinase b (NDKb) is responsible for nucleoside triphosphates synthesis and it has key role in the purine metabolism in trypanosomatid protozoans. Elongation factor 2 (EF2) is an important factor for protein synthesis. Recently, our phosphoproteomic analysis demonstrated that NDKb and EF2 proteins were phosphorylated and dephosphorylated in antimony (Sb^III)-resistant L. braziliensis line compared to its Sb^III-susceptible pair, respectively.

Methods

In this study, the overexpression of NDKb and EF2 genes in L. braziliensis and L. infantum was performed to investigate the contribution of these proteins in the Sb^III-resistance phenotype. Furthermore, we examined the role of lamivudine on Sb^III susceptibility in clones that overexpress the NDKb gene, and the effect of EF2 kinase (EF2K) inhibitor on the growth of EF2-overexpressing parasites.

Results

Western blot analysis demonstrated that NDKb and EF2 proteins are more and less expressed, respectively, in Sb^III-resistant line of L. braziliensis than its wild-type (WTS) counterpart, corroborating our previous phosphoproteomic data. NDKb or EF2-overexpressing L. braziliensis lines were 1.6 to 2.1-fold more resistant to Sb^III than the untransfected WTS line. In contrast, no difference in Sb^III susceptibility was observed in L. infantum parasites overexpressing NDKb or EF2. Susceptibility assays showed that NDKb-overexpressing L. braziliensis lines presented elevated resistance to lamivudine, an antiviral agent, but it did not alter the leishmanicidal activity in association with Sb^III. EF2-overexpressing L. braziliensis clone was slightly more resistant to EF2K inhibitor than the WTS line. Surprisingly, this inhibitor increased the antileishmanial effect of Sb^III, suggesting that this association might be a valuable strategy for leishmaniasis chemotherapy.

Conclusion

Our findings represent the first study of NDKb and EF2 genes overexpression that demonstrates an increase of Sb^III resistance in L. braziliensis which can contribute to develop new strategies for leishmaniasis treatment.

A Telomeric Cluster of Antimony Resistance Genes on Chromosome 34 of Leishmania infantum

The mechanisms underlying the drug resistance of Leishmania spp. are manifold and not completely identified. Apart from the highly conserved multidrug resistance gene family known from higher eukaryotes, Leishmania spp. also possess genus-specific resistance marker genes. One of them, ARM58, was first identified in Leishmania braziliensis using a functional cloning approach, and its domain structure was characterized in L. infantum Here we report that L. infantum ARM58 is part of a gene cluster at the telomeric end of chromosome 34 also comprising the neighboring genes ARM56 and HSP23. We show that overexpression of all three genes can confer antimony resistance to intracellular amastigotes. Upon overexpression in L. donovani, ARM58 and ARM56 are secreted via exosomes, suggesting a scavenger/secretion mechanism of action. Using a combination of functional cloning and next-generation sequencing, we found that the gene cluster was selected only under antimonyl tartrate challenge and weakly under Cu(2+) challenge but not under sodium arsenite, Cd(2+), or miltefosine challenge. The selective advantage is less pronounced in intracellular amastigotes treated with the sodium stibogluconate, possibly due to the known macrophage-stimulatory activity of this drug, against which these resistance markers may not be active. Our data point to the specificity of these three genes for antimony resistance.

Drug resistance in eukaryotic microorganisms

Eukaryotic microbial pathogens are major contributors to illness and death globally. Although much of their impact can be controlled by drug therapy as with prokaryotic microorganisms, the emergence of drug resistance has threatened these treatment efforts. Here, we discuss the challenges posed by eukaryotic microbial pathogens and how these are similar to, or differ from, the challenges of prokaryotic antibiotic resistance. The therapies used for several major eukaryotic microorganisms are then detailed, and the mechanisms that they have evolved to overcome these therapies are described. The rapid emergence of resistance and the restricted pipeline of new drug therapies pose considerable risks to global health and are particularly acute in the developing world. Nonetheless, we detail how the integration of new technology, biological understanding, epidemiology and evolutionary analysis can help sustain existing therapies, anticipate the emergence of resistance or optimize the deployment of new therapies.

Cos-Seq for high-throughput identification of drug target and resistance mechanisms in the protozoan parasite Leishmania

Innovative strategies are needed to accelerate the identification of antimicrobial drug targets and resistance mechanisms. Here we develop a sensitive method, which we term Cosmid Sequencing (or "Cos-Seq"), based on functional cloning coupled to next-generation sequencing. Cos-Seq identified >60 loci in the Leishmania genome that were enriched via drug selection with methotrexate and five major antileishmanials (antimony, miltefosine, paromomycin, amphotericin B, and pentamidine). Functional validation highlighted both known and previously unidentified drug targets and resistance genes, including novel roles for phosphatases in resistance to methotrexate and antimony, for ergosterol and phospholipid metabolism genes in resistance to miltefosine, and for hypothetical proteins in resistance to paromomycin, amphothericin B, and pentamidine. Several genes/loci were also found to confer resistance to two or more antileishmanials. This screening method will expedite the discovery of drug targets and resistance mechanisms and is easily adaptable to other microorganisms.

Mitochondrial Proteomics of Antimony and Miltefosine Resistant Leishmania infantum

Antimony (SbIII) and miltefosine (MIL) are important drugs for the treatment of Leishmania parasite infections. The mitochondrion is likely to play a central role in SbIII and MIL induced cell death in this parasite. Enriched mitochondrial samples from Leishmania promastigotes selected step by step for in vitro resistance to SbIII and MIL were subjected to differential proteomic analysis. A shared decrease in both mutants in the levels of pyruvate dehydrogenase, dihydrolipoamide dehydrogenase, and isocitrate dehydrogenase was observed, as well as a differential abundance in two calcium-binding proteins and the unique dynamin-1-like protein of the parasite. Both mutants presented a shared increase in the succinyl-CoA:3-ketoacid-coenzyme A transferase and the abundance of numerous hypothetical proteins was also altered in both mutants. In general, the proteomic changes observed in the MIL mutant were less pronounced than in the SbIII mutant, probably due to the early appearance of a mutation in the miltefosine transporter abrogating the need for a strong mitochondrial adaptation. This study is the first analysis of the Leishmania mitochondrial proteome and offers powerful insights into the adaptations to this organelle during SbIII and MIL drug resistance.

Upregulation of Cysteine Synthase and Cystathionine β-Synthase Contributes to Leishmania braziliensis Survival under Oxidative Stress

Cysteine metabolism is considered essential for the crucial maintenance of a reducing environment in trypanosomatids due to its importance as a precursor of trypanothione biosynthesis. Expression, activity, functional rescue, and overexpression of cysteine synthase (CS) and cystathionine β-synthase (CβS) were evaluated in Leishmania braziliensis promastigotes and intracellular amastigotes under in vitro stress conditions induced by hydrogen peroxide (H2O2), S-nitroso-N-acetylpenicillamine, or antimonial compounds. Our results demonstrate a stage-specific increase in the levels of protein expression and activity of L. braziliensis CS (LbrCS) and L. braziliensis CβS (LbrCβS), resulting in an increment of total thiol levels in response to both oxidative and nitrosative stress. The rescue of the CS activity in Trypanosoma rangeli, a trypanosome that does not perform cysteine biosynthesis de novo, resulted in increased rates of survival of epimastigotes expressing the LbrCS under stress conditions compared to those of wild-type parasites. We also found that the ability of L. braziliensis promastigotes and amastigotes overexpressing LbrCS and LbrCβS to resist oxidative stress was significantly enhanced compared to that of nontransfected cells, resulting in a phenotype far more resistant to treatment with the pentavalent form of Sb in vitro. In conclusion, the upregulation of protein expression and increment of the levels of LbrCS and LbrCβS activity alter parasite resistance to antimonials and may influence the efficacy of antimony treatment of New World leishmaniasis.

Advances in Development of New Treatment for Leishmaniasis

Leishmaniasis is a neglected infectious disease caused by several different species of protozoan parasites of the genus Leishmania. Current strategies to control this disease are mainly based on chemotherapy. Despite being available for the last 70 years, leishmanial chemotherapy has lack of efficiency, since its route of administration is difficult and it can cause serious side effects, which results in the emergence of resistant cases. The medical-scientific community is facing difficulties to overcome these problems with new suitable and efficient drugs, as well as the identification of new drug targets. The availability of the complete genome sequence of Leishmania has given the scientific community the possibility of large-scale analysis, which may lead to better understanding of parasite biology and consequent identification of novel drug targets. In this review we focus on how high-throughput analysis is helping us and other groups to identify novel targets for chemotherapeutic interventions. We further discuss recent data produced by our group regarding the use of the high-throughput techniques and how this helped us to identify and assess the potential of new identified targets.

Metabolomic-based strategies for anti-parasite drug discovery

Metabolomics-based studies are proving of great utility in the analysis of modes of action (MOAs) and resistance mechanisms of drugs in parasitic protozoa. They have helped to determine the MOA of eflornithine, half of the gold standard combination therapy in use against human African trypanosomiasis (HAT), as well as the mechanism of resistance to this drug. In Leishmania, metabolomics has also given insight into the MOA of miltefosine, an alkylphospholipid. Several studies on antimony resistance in Leishmania have been conducted, analyzing the metabolic content of resistant lines, offering clues as to the MOA of this class of drugs. A study of chloroquine resistance in Plasmodium falciparum combined metabolomics techniques with other genetic and proteomic techniques to offer new insight into the role of the PfCRT protein. The MOA and mechanism of resistance to a group of halogenated pyrimidines in Trypanosoma brucei have also recently been elucidated. Effective as metabolomics techniques are, care must be taken in the design and implementation of these experiments, to ensure the resulting data are meaningful. This review outlines the steps required to conduct a metabolomics experiment as well as provide an overview of metabolomics-based drug research in protozoa to date.