Modulation of Leishmania ABC protein gene expression through life stages and among drug-resistant parasites

The adaptive roles of aneuploidy and polyclonality in Leishmania in response to environmental stress

Aneuploidy is generally considered harmful, but in some microorganisms, it can act as an adaptive mechanism against environmental stress. Here, we use Leishmania-a protozoan parasite with remarkable genome plasticity-to study the early steps of aneuploidy evolution under high drug pressure (using antimony or miltefosine as stressors). By combining single-cell genomics, lineage tracing with cellular barcodes, and longitudinal genome characterization, we reveal that aneuploidy changes under antimony pressure result from polyclonal selection of pre-existing karyotypes, complemented by further and rapid de novo alterations in chromosome copy number along evolution. In the case of miltefosine, early parasite adaptation is associated with independent point mutations in a miltefosine transporter gene, while aneuploidy changes only emerge later, upon exposure to increased drug levels. Therefore, polyclonality and genome plasticity are hallmarks of parasite adaptation, but the scenario of aneuploidy dynamics depends on the nature and strength of the environmental stress as well as on the existence of other pre-adaptive mechanisms.

Translational reprogramming as a driver of antimony-drug resistance in Leishmania

Leishmania is a unicellular protozoan that has a limited transcriptional control and mostly uses post-transcriptional regulation of gene expression, although the molecular mechanisms of the process are still poorly understood. Treatments of leishmaniasis, pathologies associated with Leishmania infections, are limited due to drug resistance. Here, we report dramatic differences in mRNA translation in antimony drug-resistant and sensitive strains at the full translatome level. The major differences (2431 differentially translated transcripts) were demonstrated in the absence of the drug pressure supporting that complex preemptive adaptations are needed to efficiently compensate for the loss of biological fitness once they are exposed to the antimony. In contrast, drug-resistant parasites exposed to antimony activated a highly selective translation of only 156 transcripts. This selective mRNA translation is associated with surface protein rearrangement, optimized energy metabolism, amastins upregulation, and improved antioxidant response. We propose a novel model that establishes translational control as a major driver of antimony-resistant phenotypes in Leishmania.

Thiophene derivatives activity against the protozoan parasite Leishmania infantum

Treatments against leishmaniasis are limited and the development of new molecules is crucial. One class of developmental drug that has shown activity against the parasite Leishmania are thiophene derivatives. Here we synthetized thirty-eight novel thiophene compounds and characterized their activity and potential for resistance against L. infantum. Half of the molecules had an EC₅₀ in the low micromolar range, the piperidine derivatives being more potent than the tetramethylpyran derivatives. Resistance was challenging to select for, and resistant cells could only be raised against one (GC1-19) of the four most active compounds. Using chemogenomic screens we show that a gene conversion event at the ABCG2 locus as well as the overexpression of a tryparedoxin peroxidase are responsible for a weak but significant resistance to the GC1-19 drug candidate. Together, our results suggest that thiophene is a scaffold of interest for further drug development against leishmaniasis.

Repurposing the Antibacterial Agents Peptide 19-4LF and Peptide 19-2.5 for Treatment of Cutaneous Leishmaniasis

The lack of safe and cost-effective treatments against leishmaniasis highlights the urgent need to develop improved leishmanicidal agents. Antimicrobial peptides (AMPs) are an emerging category of therapeutics exerting a wide range of biological activities such as anti-bacterial, anti-fungal, anti-parasitic and anti-tumoral. In the present study, the approach of repurposing AMPs as antileishmanial drugs was applied. The leishmanicidal activity of two synthetic anti-lipopolysaccharide peptides (SALPs), so-called 19-2.5 and 19-4LF was characterized in Leishmania major. In vitro, both peptides were highly active against intracellular Leishmania major in mouse macrophages without exerting toxicity in host cells. Then, q-PCR-based gene profiling, revealed that this activity was related to the downregulation of several genes involved in drug resistance (yip1), virulence (gp63) and parasite proliferation (Cyclin 1 and Cyclin 6). Importantly, the treatment of BALB/c mice with any of the two AMPs caused a significant reduction in L. major infective burden. This effect was associated with an increase in Th1 cytokine levels (IL-12p35, TNF-α, and iNOS) in the skin lesion and spleen of the L. major infected mice while the Th2-associated genes were downregulated (IL-4 and IL-6). Lastly, we investigated the effect of both peptides in the gene expression profile of the P2X7 purinergic receptor, which has been reported as a therapeutic target in several diseases. The results showed significant repression of P2X7R by both peptides in the skin lesion of L. major infected mice to an extent comparable to that of a common anti-leishmanial drug, Paromomycin. Our in vitro and in vivo studies suggest that the synthetic AMPs 19-2.5 and 19-4LF are promising candidates for leishmaniasis treatment and present P2X7R as a potential therapeutic target in cutaneous leishmaniasis (CL).

Three types of Leishmania mexicana amastigotes: Proteome comparison by quantitative proteomic analysis

Leishmania is the unicellular parasite transmitted by phlebotomine sand fly bite. It exists in two different forms; extracellular promastigotes, occurring in the gut of sand flies, and intracellular, round-shaped amastigotes residing mainly in vertebrate macrophages. As amastigotes originating from infected animals are often present in insufficient quality and quantity, two alternative types of amastigotes were introduced for laboratory experiments: axenic amastigotes and amastigotes from macrophages infected in vitro. Nevertheless, there is very little information about the degree of similarity/difference among these three types of amastigotes on proteomic level, whose comparison is crucial for assessing the suitability of using alternative types of amastigotes in experiments. In this study, L. mexicana amastigotes obtained from lesion of infected BALB/c mice were proteomically compared with alternatively cultivated amastigotes (axenic and macrophage-derived ones). Amastigotes of all three types were isolated, individually treated and analysed by LC-MS/MS proteomic analysis with quantification using TMT10-plex isobaric labeling. Significant differences were observed in the abundance of metabolic enzymes, virulence factors and proteins involved in translation and condensation of DNA. The most pronounced differences were observed between axenic amastigotes and lesion-derived amastigotes, macrophage-derived amastigotes were mostly intermediate between axenic and lesion-derived ones.

The History of the ABC Proteins in Human Trypanosomiasis Pathogens

Human trypanosomiasis affects nearly eight million people worldwide, causing great economic and social impact, mainly in endemic areas. T. cruzi and T. brucei are protozoan parasites that present efficient mechanisms of immune system evasion, leading to disease chronification. Currently, there is no vaccine, and chemotherapy is effective only in the absence of severe clinical manifestations. Nevertheless, resistant phenotypes to chemotherapy have been described in protozoan parasites, associated with cross-resistance to other chemically unrelated drugs. Multidrug resistance is multifactorial, involving: (i) drug entry, (ii) activation, (iii) metabolism and (iv) efflux pathways. In this context, ABC transporters, initially discovered in resistant tumor cells, have drawn attention in protozoan parasites, owing to their ability to decrease drug accumulation, thus mitigating their toxic effects. The discovery of these transporters in the Trypanosomatidae family started in the 1990s; however, few members were described and functionally characterized. This review contains a brief history of the main ABC transporters involved in resistance that propelled their investigation in Trypanosoma species, the main efflux modulators, as well as ABC genes described in T. cruzi and T. brucei according to the nomenclature HUGO. We hope to convey the importance that ABC transporters play in parasite physiology and chemotherapy resistance.

- Amphotericin B: A drug of choice for Visceral Leishmaniasis

Visceral leishmaniasis or Kala-azar is a vector-borne disease caused by an intracellular parasite of the genus leishmania. In India, Amphotericin B (AmB) is a first-line medication for treating leishmaniasis. After a large-scale resistance to pentavalent antimony therapy developed in Bihar state, it was rediscovered as an effective treatment for Leishmania donovani infection. AmB which binds to the ergosterol of protozoan cells causes a change in membrane integrity resulting in ions leakage, and ultimately leading to cell death. The treatment effect of liposomal AmB can be seen more quickly than deoxycholate AmB because, it has some toxic effects, but liposomal AmB is significantly less toxic. Evidence from studies suggested that ABLC (Abelcet) and ABCD (Amphotec) are as effective as L-AmB but Liposomal form (Ambisome) is a more widely accepted treatment option than conventional ones. Nevertheless, the world needs some way more efficient antileishmanial drugs that are less toxic and less expensive for people living with parasitic infections caused by Leishmania. So, academics, researchers, and sponsors need to focus on finding such drugs. This review provides a summary of the chemical, pharmacokinetic, drug-target interactions, stability, dose efficacy, and many other characteristics of the AmB and their various formulations. We have also highlighted the clinically significant aspects of PKDL and VL co-infection with HIV/TB.

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.

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.

Profilin is involved in G1 to S phase progression and mitotic spindle orientation during Leishmania donovani cell division cycle

Profilin is a multi-ligand binding protein, which is a key regulator of actin dynamics and involved in regulating several cellular functions. It is present in all eukaryotes, including trypanosomatids such as Leishmania. However, not much is known about its functions in these organisms. Our earlier studies have shown that Leishmania parasites express a single homologue of profilin (LdPfn) that binds actin, phosphoinositides and poly- L- proline motives, and depletion of its intracellular pool to 50%of normal levels affects the cell growth and intracellular trafficking. Here, we show, employing affinity pull-down and mass spectroscopy, that LdPfn interacted with a large number of proteins, including those involved in mRNA processing and protein translation initiation, such as eIF4A1. Further, we reveal, using mRNA Seq analysis, that depletion of LdPfn in Leishmania cells (LdPfn^+/-) resulted in significantly reduced expression of genes which encode proteins involved in cell cycle regulation, mRNA translation initiation, nucleosides and amino acids transport. In addition, we show that in LdPfn^+/- cells, cellular levels of eIF4A1 protein were significantly decreased, and during their cell division cycle, G1-to-S phase progression was delayed and orientation of mitotic spindle altered. These changes were, however, reversed to normal by episomal expression of GFP-LdPfn in LdPfn^+/- cells. Taken together, our results indicate that profilin is involved in regulation of G1-to-S phase progression and mitotic spindle orientation in Leishmania cell cycle, perhaps through its interaction with elF4A1 protein.

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.

Intrinsic and Chemotherapeutic Stressors Modulate ABCC-Like Transport in Trypanosoma cruzi

Trypanosoma cruzi is the etiologic agent for Chagas disease, which affects 6-7 million people worldwide. The biological diversity of the parasite reflects on inefficiency of benznidazole, which is a first choice chemotherapy, on chronic patients. ABC transporters that extrude xenobiotics, metabolites, and mediators are overexpressed in resistant cells and contribute to chemotherapy failure. An ABCC-like transport was identified in the Y strain and extrudes thiol-conjugated compounds. As thiols represent a line of defense towards reactive species, we aimed to verify whether ABCC-like transport could participate in the regulation of responses to stressor stimuli. In order to achieve this, ABCC-like activity was measured by flow cytometry using fluorescent substrates. The present study reveals the participation of glutathione and ceramides on ABCC-like transport, which are both implicated in stress. Hemin modulated the ABCC-like efflux which suggests that this protein might be involved in cellular detoxification. Additionally, all strains evaluated exhibited ABCC-like activity, while no ABCB1-like activity was detected. Results suggest that ABCC-like efflux is not associated with natural resistance to benznidazole, since sensitive strains showed higher activity than the resistant ones. Although benznidazole is not a direct substrate, ABCC-like efflux increased after prolonged drug exposure and this indicates that the ABCC-like efflux mediated protection against cell stress depends on the glutathione biosynthesis pathway.

Deposition of pentamidine analogues in the human body - spectroscopic and computational approaches

Bis-benzamidines are a diverse group of compounds with high potential in pharmacotherapy, and among them, pentamidine is a drug of great therapeutic significance in Pneumocystis jiroveci pneumonia (PJP) prophylaxis and therapy. Pharmacokinetic properties of these cationic species such as transport, acid/base equilibria, and interactions with potential target molecules are still of interest, especially for recently designed compounds. To broaden our knowledge drug-likeness, human serum albumin binding, and acidity constants (K_a) were experimentally and theoretically examined for five pentamidine analogues 1 - 5 with -NH-CO-chain-CO-NH-bridges of increasing length and O, N, and S atoms in the chain. The studied analogues display very marked activity against Pneumocystis carinii without cytotoxicity that inspired us to perform an in silico analysis of their mode of action based on the hypothesis that the small DNA groove of rich in adenine-thymine pairs is their molecular target. These studies allowed us to classify them as very promising lead molecules.

Determinants of Unresponsiveness to Treatment in Cutaneous Leishmaniasis: A Focus on Anthroponotic Form Due to Leishmania tropica

Cutaneous leishmaniasis (CL) is a curable disease; however, due to various risk factors, unresponsiveness to CL treatments is inevitable. The treatment of CL has been firmly correlated with multiple determinants, such as demographical, clinical, and environmental factors, the host’s immune response, poor treatment adherence, the parasite’s genetic make-up, and Leishmania RNA virus. This study primarily focuses on the risk factors associated with different therapeutic outcomes following meglumine antimoniate (MA; Glucantime^®) treatment and policy approaches to prevent unresponsiveness in CL patients with a focus on anthroponotic form (ACL). Findings suggest that effective preventive and therapeutic measures should be more vigorously implemented, particularly in endemic areas. Accordingly, extensive training is essential to monitor drug unresponsiveness regularly, especially in tropical regions where the disease is prevalent. Since humans are the fundamental reservoir host of ACL due to L. tropica, prompt detection, early diagnosis, and timely and effective treatment could help control this disease. Furthermore, major challenges and gaps remain: efficacious vaccine, new tools, and expert staff are crucial before CL can be definitively controlled.

ABC-finder: A containerized web server for the identification and topology prediction of ABC proteins

In view of the multiple clinical and physiological implications of ABC transporter proteins, there is a considerable interest among researchers to characterize them functionally. However, such characterizations are based on the premise that ABC proteins are accurately identified in the proteome of an organism, and their topology is correctly predicted. With this objective, we have developed ABC-finder, i.e., a Docker-based package for the identification of ABC proteins in all organisms, and visualization of the topology of ABC proteins using a web browser. ABC-finder is built and deployed in a Linux container, making it scalable for many concurrent users on our servers and enabling users to download and run it locally. Overall, ABC-finder is a convenient, portable, and platform-independent tool for the identification and topology prediction of ABC proteins. ABC-finder is accessible at http://abc-finder.osdd.jnu.ac.in.

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.

LmjF.22.0810 from Leishmania major Modulates the Th2-Type Immune Response and Is Involved in Leishmaniasis Outcome

A novel serine/threonine protein kinase, LmjF.22.0810, was recently described in Leishmania major. After generating an L. major cell line overexpressing LmjF.22.0810 (named LmJ3OE), the ability of this novel protein to modulate the Th2-type immune response was analyzed. Our results suggest that the protein kinase LmjF.22.0810 might be involved in leishmaniasis outcomes. Indeed, our study outlined the LmJ3OE parasites infectivity in vitro and in vivo. Transgenic parasites displayed lower phagocytosis rates in vitro, and their promastigote forms exhibited lower expression levels of virulence factors compared to their counterparts in control parasites. In addition, LmJ3OE parasites developed significantly smaller footpad swelling in susceptible BALB/c mice. Hematoxylin-eosin staining allowed the observation of a lower inflammatory infiltrate in the footpad from LmJ3OE-infected mice compared to animals inoculated with control parasites. Gene expression of Th2-associated cytokines and effectors revealed a dramatically lower induction in interleukin (IL)-4, IL-10, and arginase 1 (ARG1) mRNA levels at the beginning of the swelling; no expression change was found in Th1-associated cytokines except for IL-12. Accordingly, such results were validated by immunohistochemistry studies, illustrating a weaker expression of ARG1 and a similar induction for inducible NO synthase (iNOS) in footpads from LmJ3OE-infected mice compared to control L. major infected animals. Furthermore, the parasite burden was lower in footpads from LmJ3OE-infected mice. Our analysis indicated that such significant smaller footpad swellings might be due to an impairment of the Th2 immune response that subsequently benefits Th1 prevalence. Altogether, these studies depict LmjF.22.0810 as a potential modulator of host immune responses to Leishmania. Finally, this promising target might be involved in the modulation of infection outcome.

Assessing the composition of the plasma membrane of Leishmania (Leishmania) infantum and L. (L.) amazonensis using label-free proteomics

Protozoan parasites of the genus Leishmania are causative agents of leishmaniasis, a wide range of diseases affecting 12 million people worldwide. The species L. infantum and L. amazonensis are etiologic agents of visceral and cutaneous leishmaniasis, respectively. Most proteome analyses of Leishmania have been carried out on whole-cell extracts, but such an approach tends to underrepresent membrane-associated proteins due to their high hydrophobicity and low solubility. Considering the relevance of this category of proteins in virulence, invasiveness and the host-parasite interface, this study applied label-free proteomics to assess the plasma membrane sub-proteome of L. infantum and L. amazonensis. The number of proteins identified in L. infantum and L. amazonensis promastigotes was 1168 and 1455, respectively. After rigorous data processing and mining, 157 proteins were classified as putative plasma membrane-associated proteins, of which 56 proteins were detected in both species, six proteins were detected only in L. infantum and 39 proteins were exclusive to L. amazonensis. The quantitative analysis revealed that two proteins were more abundant in L. infantum, including the glucose transporter 2, and five proteins were more abundant in L. amazonensis. The identified proteins associated with distinct processes and functions. In this regard, proteins of L. infantum were linked to metabolic processes whereas L. amazonensis proteins were involved in signal transduction. Moreover, transmembrane transport was a significant process among the group of proteins detected in both species and members of the superfamily of ABC transporters were highly represented. Interestingly, some proteins of this family were solely detected in L. amazonensis, such as ABCA9. GP63, a well-known virulence factor, was the only GPI-anchored protein identified in the membrane preparations of both species. Finally, we found several proteins with uncharacterized functions, including differentially abundant ones, highlighting a gap in the study of Leishmania proteins. Proteins characterization could provide a better biological understanding of these parasites and deliver new possibilities regarding the discovery of therapeutic targets, drug resistance and vaccine candidates.

Minor temperature shifts do not affect chromosomal ploidy but cause transcriptomic changes in Leishmania braziliensis promastigotes in vitro

BACKGROUND

The leishmaniases are complex neglected diseases caused by protozoan parasites of the genus Leishmania. Leishmania braziliensis is the main etiological agent of cutaneous leishmaniasis in the New World. In recent studies, genomic changes such as chromosome and gene copy number variations (CNVs), as well as transcriptomic changes have been highlighted as mechanisms used by Leishmania species to adapt to stress situations.

OBJECTIVES

The aim of this study was to determine the effect of short-term minor temperature shifts in the genomic and transcriptomic responses of L. braziliensis promastigotes in vitro.

METHODS

Growth curves, genome and transcriptome sequencing of L. braziliensis promastigotes were conducted from cultures exposed to three different temperatures (24ºC, 28ºC and 30ºC) compared with the control temperature (26ºC).

FINDINGS

Our results showed a decrease in L. braziliensis proliferation at 30ºC, with around 3% of the genes showing CNVs at each temperature, and transcriptomic changes in genes encoding amastin surface-like proteins, heat shock proteins and transport proteins, which may indicate a direct response to temperature stress.

MAIN CONCLUSIONS

This study provides evidence that L. braziliensis promastigotes exhibit a decrease in cell density, and noticeable changes in the transcriptomic profiles. However, there were not perceptible changes at chromosome CNVs and only ~3% of the genes changed their copies in each treatment.

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.

Intraspecific Genomic Divergence and Minor Structural Variations in Leishmania (Viannia) panamensis

Leishmania (Viannia) panamensis is one of the most important Leishmania species associated with cutaneous leishmaniasis (CL) in Latin America. Despite its wide geographic distribution and pathogenic potential in humans and animals, the genomic variability of this species is low compared with other Leishmania species circulating in the same geographical area. No studies have reported a detailed analysis of the whole genome of L. panamensis from clinical isolates using DNA high-throughput sequencing to clarify its intraspecific genomic variability or plausible divergence. Therefore, this study aimed to evaluate the intraspecific genomic variability of L. panamensis from Colombia and Panama. A total of 22 genomes were analyzed, 19 from Colombian patients with CL and three genomes from Panama obtained from public databases. The phylogenomic analysis revealed the potential existence of three well-supported clades as evidence of intraspecific divergence. Additionally, the whole-genome analysis showed low structural variations in terms of ploidy, copy number variations, and single-nucleotide polymorphisms (SNPs). SNPs shared among all clades were identified, revealing their importance in different biological processes of L. panamensis. The findings not only expand our knowledge of intraspecific genomic variability of one of the most important Leishmania species in South America but also highlights the possible existence of different clades/lineages/subpopulations across a geographic scale.

Mapping Alterations Induced by Long-Term Axenic Cultivation of Leishmania amazonensis Promastigotes With a Multiplatform Metabolomic Fingerprint Approach

Leishmaniases are widespread neglected diseases with an incidence of 1.6 million new cases and 40 thousand deaths per year. Leishmania parasites may show distinct, species-specific patterns of virulence that lead to different clinical manifestations. It is well known that successive in vitro passages (SIVP) lead to the attenuation of virulence, but neither the metabolism nor the pathways involved in these processes are well understood. Herein, promastigotes of a virulent L. amazonensis strain recently isolated from mice was compared to SIVP derived and attenuated promastigotes, submitted to 10, 40, and 60 axenic passages and named R10, R40, and R60, respectively. In vitro assays and in vivo tests were performed to characterize and confirmed the attenuation profiles. A metabolomic fingerprint comparison of R0, R10, and R60 was performed by means of capillary electrophoresis, liquid and gas chromatography coupled to mass spectrometry. To validate the metabolomic data, qPCR for selected loci, flow cytometry to measure aPS exposure, sensitivity to antimony tartrate and ROS production assays were conducted. The 65 identified metabolites were clustered in biochemical categories and mapped in eight metabolic pathways: ABC transporters; fatty acid biosynthesis; glycine, serine and threonine metabolism; β-alanine metabolism; glutathione metabolism; oxidative phosphorylation; glycerophospholipid metabolism and lysine degradation. The obtained metabolomic data correlated with previous proteomic findings of the SVIP parasites and the gene expression of 13 selected targets. Late SIVP cultures were more sensitive to Sb^III produced more ROS and exposed less phosphatidylserine in their surface. The correspondent pathways were connected to build a biochemical map of the most significant alterations involved with the process of attenuation of L. amazonensis. Overall, the reported data pointed out to a very dynamic and continuous metabolic reprogramming process, accompanied by changes in energetic, lipid and redox metabolisms, membrane remodeling and reshaping of parasite-host cells interactions, causing impacts in chemotaxis, host inflammatory responses and infectivity at the early stages of infection.

Coupling chemical mutagenesis to next generation sequencing for the identification of drug resistance mutations in Leishmania

Current genome-wide screens allow system-wide study of drug resistance but detecting small nucleotide variants (SNVs) is challenging. Here, we use chemical mutagenesis, drug selection and next generation sequencing to characterize miltefosine and paromomycin resistant clones of the parasite Leishmania. We highlight several genes involved in drug resistance by sequencing the genomes of 41 resistant clones and by concentrating on recurrent SNVs. We associate genes linked to lipid metabolism or to ribosome/translation functions with miltefosine or paromomycin resistance, respectively. We prove by allelic replacement and CRISPR-Cas9 gene-editing that the essential protein kinase CDPK1 is crucial for paromomycin resistance. We have linked CDPK1 in translation by functional interactome analysis, and provide evidence that CDPK1 contributes to antimonial resistance in the parasite. This screen is powerful in exploring networks of drug resistance in an organism with diploid to mosaic aneuploid genome, hence widening the scope of its applicability.

Here, Bhattacharya et al. chemically mutagenize Leishmania and identify genes associated with resistance to miltefosine and paromomycin by next generation sequencing. The study shows that a protein kinase (CDPK1) can mediate resistance to paromomycin by affecting translation.

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.

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

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

Drug Resistance in Protozoan Parasites: An Incessant Wrestle for Survival

Nowadays, drug resistance in parasites is considered to be one of the foremost concerns in health and disease management. It is interconnected worldwide and undermines the health of millions of people, threatening to grow worse. Unfortunately, it does not receive serious attention from every corner of society. Consequently, drug resistance in parasites is gradually complicating and challenging the treatment of parasitic diseases. In this context, we have dedicated ourselves to review the incidence of drug resistance in the protozoan parasites Plasmodium, Leishmania, Trypanosoma, Entamoeba and Toxoplasma gondii. Moreover, understanding the role of ATP-binding cassette (ABC) transporters in drug resistance is essential in the control of parasitic diseases. Therefore, we also focused on the involvement of ABC transporters in drug resistance, which will be a superior approach to find ways for better regulation of diseases caused by parasitic infections.

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.

Functional Characterization of ABCC Proteins from Trypanosoma cruzi and Their Involvement with Thiol Transport

Chagas disease is a neglected disease caused by the protozoan Trypanosoma cruzi and affects 8 million people worldwide. The main chemotherapy is based on benznidazole. The efficacy in the treatment depends on factors such as the parasite strain, which may present different sensitivity to treatment. In this context, the expression of ABC transporters has been related to chemotherapy failure. ABC transporters share a well-conserved ABC domain, responsible for ATP binding and hydrolysis, whose the energy released is coupled to transport of molecules through membranes. The most known ABC transporters are ABCB1 and ABCC1, involved in the multidrug resistance phenotype in cancer, given their participation in cellular detoxification. In T. cruzi, 27 ABC genes were identified in the genome. Nonetheless, only four ABC genes were characterized: ABCA3, involved in vesicular trafficking; ABCG1, overexpressed in strains naturally resistant to benznidazole, and P-glycoprotein 1 and 2, whose participation in drug resistance is controversial. Considering P-glycoprotein genes are related to ABCC subfamily in T. cruzi according to the demonstration using BLASTP alignment, we evaluated both ABCB1-like and ABCC-like activities in epimastigote and trypomastigote forms of the Y strain. The transport activities were evaluated by the efflux of the fluorescent dyes Rhodamine 123 and Carboxyfluorescein in a flow cytometer. Results indicated that there was no ABCB1-like activity in both T. cruzi forms. Conversely, results demonstrated ABCC-like activity in both epimastigote and trypomastigote forms of T. cruzi. This activity was inhibited by ABCC transport modulators (probenecid, indomethacin, and MK-571), by ATP-depleting agents (sodium azide and iodoacetic acid) and by the thiol-depleting agent N-ethylmaleimide. Additionally, the presence of ABCC-like activity was supported by direct inhibition of the thiol-conjugated compound efflux with indomethacin, characteristic of ABCC subfamily members. Taken together, the results provide the first description of native ABCC-like activity in T. cruzi epimastigote and trypomastigote forms, indicating that the study of the biological role for that thiol transporter is crucial to reveal new molecular mechanisms for therapeutic approaches in the Chagas disease.

Profiling gene expression of antimony response genes in Leishmania (Viannia) panamensis and infected macrophages and its relationship with drug susceptibility

The mechanisms of Leishmania resistance to antimonials have been primarily determined in experimentally derived Leishmania strains. However, their participation in the susceptibility phenotype in field isolates has not been conclusively established. Being an intracellular parasite, the activity of antileishmanials is dependent on internalization of drugs into host cells and effective delivery to the intracellular compartments inhabited by the parasite. In this study we quantified and comparatively analyzed the gene expression of nine molecules involved in mechanisms of xenobiotic detoxification and Leishmania resistance to antimonial drugs in resistant and susceptible laboratory derived and clinical L.(Viannia) panamensis strains(n=19). In addition, we explored the impact of Leishmania susceptibility to antimonials on the expression of macrophage gene products having putative functions in transport, accumulation and metabolism of antimonials. As previously shown for other Leishmania species, a trend of increased abcc3 and lower aqp-1 expression was observed in the laboratory derived Sb-resistant L.(V.) panamensis line. However, this was not found in clinical strains, in which the expression of abca2 was significantly higher in resistant strains as both, promastigotes and intracellular amastigotes. The effect of drug susceptibility on host cell gene expression was evaluated on primary human macrophages from patients with cutaneous leishmaniasis (n=17) infected ex-vivo with the matched L.(V.) panamensis strains isolated at diagnosis, and in THP-1 cells infected with clinical strains (n=6) and laboratory adapted L.(V.) panamensis lines. Four molecules, abcb1 (p-gp), abcb6, aqp-9 and mt2a were differentially modulated by drug resistant and susceptible parasites, and among these, a consistent and significantly increased expression of the xenobiotic scavenging molecule mt2a was observed in macrophages infected with Sb-susceptible L. (V.) panamensis. Our results substantiate that different mechanisms of drug resistance operate in laboratory adapted and clinical Leishmania strains, and provide evidence that parasite-mediated modulation of host cell gene expression of molecules involved in drug transport and metabolism could contribute to the mechanisms of drug resistance and susceptibility in Leishmania.

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.

ABCI3 Is a New Mitochondrial ABC Transporter from Leishmania major Involved in Susceptibility to Antimonials and Infectivity

We have identified and characterized ABCI3 as a new mitochondrial ABC transporter from Leishmania major Localization studies using confocal microscopy, a surface biotinylation assay, and trypsin digestion after digitonin permeabilization suggested that ABCI3 presents a dual localization in both mitochondria and the plasma membrane. From studies using parasites with a single knockout of ABCI3 (ABCI3^+/-), we provide evidence that ABCI3 is directly involved in susceptibility to the trivalent form of antimony (SbIII) and metal ions. Attempts to obtain parasites with a double knockout of ABCI3 were unsuccessful, suggesting that ABCI3 could be an essential gene in L. majorABCI3^+/- promastigotes were 5-fold more resistant to Sb^III than the wild type, while ABCI3^+/- amastigotes were approximately 2-fold more resistant to pentavalent antimony (Sb^V). This resistance phenotype was associated with decreased Sb^III accumulation due to decreased Sb^III uptake. ABCI3^+/- parasites presented higher ATP levels and generated less mitochondrial superoxide after Sb^III incubation. Finally, we observed that ABCI3^+/- parasites showed a slightly higher infection capacity than wild-type and add-back ABCI3^+/-::3×FABCI3 parasites; however, after 72 h the number of ABCI3^+/- intracellular parasites per macrophage increased significantly. Our results show that ABCI3 is responsible for Sb^III transport inside mitochondria, where it contributes to enhancement of the general toxic effects caused by Sb^III To our knowledge, ABCI3 is the first ABC transporter which is involved in susceptibility toward antimony, conferring Sb^III resistance to parasites when it is partially deleted.

Gene Expression in Leishmania Is Regulated Predominantly by Gene Dosage

Leishmania tropica, a unicellular eukaryotic parasite present in North and East Africa, the Middle East, and the Indian subcontinent, has been linked to large outbreaks of cutaneous leishmaniasis in displaced populations in Iraq, Jordan, and Syria. Here, we report the genome sequence of this pathogen and 7,863 identified protein-coding genes, and we show that the majority of clinical isolates possess high levels of allelic diversity, genetic admixture, heterozygosity, and extensive aneuploidy. By utilizing paired genome-wide high-throughput DNA sequencing (DNA-seq) with RNA-seq, we found that gene dosage, at the level of individual genes or chromosomal "somy" (a general term covering disomy, trisomy, tetrasomy, etc.), accounted for greater than 85% of total gene expression variation in genes with a 2-fold or greater change in expression. High gene copy number variation (CNV) among membrane-bound transporters, a class of proteins previously implicated in drug resistance, was found for the most highly differentially expressed genes. Our results suggest that gene dosage is an adaptive trait that confers phenotypic plasticity among natural Leishmania populations by rapid down- or upregulation of transporter proteins to limit the effects of environmental stresses, such as drug selection.IMPORTANCELeishmania is a genus of unicellular eukaryotic parasites that is responsible for a spectrum of human diseases that range from cutaneous leishmaniasis (CL) and mucocutaneous leishmaniasis (MCL) to life-threatening visceral leishmaniasis (VL). Developmental and strain-specific gene expression is largely thought to be due to mRNA message stability or posttranscriptional regulatory networks for this species, whose genome is organized into polycistronic gene clusters in the absence of promoter-mediated regulation of transcription initiation of nuclear genes. Genetic hybridization has been demonstrated to yield dramatic structural genomic variation, but whether such changes in gene dosage impact gene expression has not been formally investigated. Here we show that the predominant mechanism determining transcript abundance differences (>85%) in Leishmania tropica is that of gene dosage at the level of individual genes or chromosomal somy.

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.

Organic cation transporter 1 (OCT1) is involved in pentamidine transport at the human and mouse blood-brain barrier (BBB)

Pentamidine is an effective trypanocidal drug used against stage 1 Human African Trypanosomiasis (HAT). At the blood-brain barrier (BBB), it accumulates inside the endothelial cells but has limited entry into the brain. This study examined transporters involved in pentamidine transport at the human and mouse BBB using hCMEC/D3 and bEnd.3 cell lines, respectively. Results revealed that both cell lines expressed the organic cation transporters (OCT1, OCT2 and OCT3), however, P-gp was only expressed in hCMEC/D3 cells. Polarised expression of OCT1 was also observed. Functional assays found that ATP depletion significantly increased [³H]pentamidine accumulation in hCMEC/D3 cells (***p<0.001) but not in bEnd.3 cells. Incubation with unlabelled pentamidine significantly decreased accumulation in hCMEC/D3 and bEnd.3 cells after 120 minutes (***p<0.001). Treating both cell lines with haloperidol and amantadine also decreased [³H]pentamidine accumulation significantly (***p<0.001 and **p<0.01 respectively). However, prazosin treatment decreased [³H]pentamidine accumulation only in hCMEC/D3 cells (*p<0.05), and not bEnd.3 cells. Furthermore, the presence of OCTN, MATE, PMAT, ENT or CNT inhibitors/substrates had no significant effect on the accumulation of [³H]pentamidine in both cell lines. From the data, we conclude that pentamidine interacts with multiple transporters, is taken into brain endothelial cells by OCT1 transporter and is extruded into the blood by ATP-dependent mechanisms. These interactions along with the predominant presence of OCT1 in the luminal membrane of the BBB contribute to the limited entry of pentamidine into the brain. This information is of key importance to the development of pentamidine based combination therapies which could be used to treat CNS stage HAT by improving CNS delivery, efficacy against trypanosomes and safety profile of pentamidine.

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.

LmABCB3, an atypical mitochondrial ABC transporter essential for Leishmania major virulence, acts in heme and cytosolic iron/sulfur clusters biogenesis

Background

Mitochondria play essential biological functions including the synthesis and trafficking of porphyrins and iron/sulfur clusters (ISC), processes that in mammals involve the mitochondrial ATP-Binding Cassette (ABC) transporters ABCB6 and ABCB7, respectively. The mitochondrion of pathogenic protozoan parasites such as Leishmania is a promising goal for new therapeutic approaches. Leishmania infects human macrophages producing the neglected tropical disease known as leishmaniasis. Like most trypanosomatid parasites, Leishmania is auxotrophous for heme and must acquire porphyrins from the host.

Methods

LmABCB3, a new Leishmania major protein with significant sequence similarity to human ABCB6/ABCB7, was identified and characterized using bioinformatic tools. Fluorescent microscopy was used to determine its cellular localization, and its level of expression was modulated by molecular genetic techniques. Intracellular in vitro assays were used to demonstrate its role in amastigotes replication, and an in vivo mouse model was used to analyze its role in virulence. Functional characterization of LmABCB3 was carried out in Leishmania promastigotes and Saccharomyces cerevisiae. Structural analysis of LmABCB3 was performed using molecular modeling software.

Results

LmABCB3 is an atypical ABC half-transporter that has a unique N-terminal extension not found in any other known ABC protein. This extension is required to target LmABCB3 to the mitochondrion and includes a potential metal-binding domain. We have shown that LmABCB3 interacts with porphyrins and is required for the mitochondrial synthesis of heme from a host precursor. We also present data supporting a role for LmABCB3 in the biogenesis of cytosolic ISC, essential cofactors for cell viability in all three kingdoms of life. LmABCB3 fully complemented the severe growth defect shown in yeast lacking ATM1, an orthologue of human ABCB7 involved in exporting from the mitochondria a gluthatione-containing compound required for the generation of cytosolic ISC. Indeed, docking analyzes performed with a LmABCB3 structural model using trypanothione, the main thiol in this parasite, as a ligand showed how both, LmABCB3 and yeast ATM1, contain a similar thiol-binding pocket. Additionally, we show solid evidence suggesting that LmABCB3 is an essential gene as dominant negative inhibition of LmABCB3 is lethal for the parasite. Moreover, the abrogation of only one allele of the gene did not impede promastigote growth in axenic culture but prevented the replication of intracellular amastigotes and the virulence of the parasites in a mouse model of cutaneous leishmaniasis.

Conclusions

Altogether our results present the previously undescribed LmABCB3 as an unusual mitochondrial ABC transporter essential for Leishmania survival through its role in the generation of heme and cytosolic ISC. Hence, LmABCB3 could represent a novel target to combat leishmaniasis.

Electronic supplementary material

The online version of this article (doi:10.1186/s13071-015-1284-5) contains supplementary material, which is available to authorized users.

Leishmanicidal activities of novel methylseleno-imidocarbamates

The generation of new antileishmanial drugs has become a priority. Selenium and its derivatives stand out as having promising leishmanicidal activity. In fact, some parasites express selenoproteins and metabolize selenium. Recently, selenium derivatives have shown the potential to reduce parasitemia, clinical manifestations, and mortality in parasite-infected mice. In this paper, after selecting four candidates according to drug similarity parameters, we observed that two of them, called compounds 2b [methyl-N,N'-di(thien-2-ylcarbonyl)-imidoselenocarbamate] and 4b [methyl-N,N'-di(5-nitrothien-3-ylcarbonyl)-imidoselenocarbamate], exhibit low 50% inhibitory concentrations (IC50s) (<3 μM) and good selectivity indexes (SIs) (>5) in Leishmania major promastigotes and lack toxicity on macrophages. In addition, in analysis of their therapeutic potential against L. major in vitro infection, both compounds display a dramatic reduction of amastigote burden (∼80%) with sublethal concentrations. Furthermore, in macrophages, these selenocompounds induce nitric oxide production, which has been described to be critical for defense against intracellular pathogens. Compounds 2b and 4b were demonstrated to cause cell cycle arrest in G1. Interestingly, evaluation of expression of genes related to proliferation (PCNA), treatment resistance (ABC transporter and alpha-tubulin), and virulence (quinonoid dihydropteridine reductase [QDPR]) showed several alterations in gene expression profiling. All these results prompt us to propose both compounds as candidates to treat leishmanial infections.

A novel ABCG-like transporter of Trypanosoma cruzi is involved in natural resistance to benznidazole

Benznidazole (BZ) is one of the two drugs used for Chagas disease treatment. Nevertheless therapeutic failures of BZ have been reported, which were mostly attributed to variable drug susceptibility among Trypanosoma cruzi strains. ATP-binding cassette (ABC) transporters are involved in a variety of translocation processes and some members have been implicated in drug resistance. Here we report the characterisation of the first T. cruzi ABCG transporter gene, named TcABCG1, which is over-expressed in parasite strains naturally resistant to BZ. Comparison of TcABCG1 gene sequence of two TcI BZ-resistant strains with CL Brener BZ-susceptible strain showed several single nucleotide polymorphisms, which determined 11 amino acid changes. CL Brener transfected with TcI transporter genes showed 40-47% increased resistance to BZ, whereas no statistical significant increment in drug resistance was observed when CL Brener was transfected with the homologous gene. Only in the parasites transfected with TcI genes there was 2-2.6-fold increased abundance of TcABCG1 transporter protein. The analysis in wild type strains also suggests that the level of TcABCG1 transporter is related to BZ natural resistance. The characteristics of untranslated regions of TcABCG1 genes of BZ-susceptible and resistant strains were investigated by computational tools.

ABCG-like transporter of Trypanosoma cruzi involved in benznidazole resistance: gene polymorphisms disclose inter-strain intragenic recombination in hybrid isolates

Benznidazole (BZ) is one of the two drugs for Chagas disease treatment. In a previous study we showed that the Trypanosoma cruzi ABCG-like transporter gene, named TcABCG1, is over-expressed in parasite strains naturally resistant to BZ and that the gene of TcI BZ-resistant strains exhibited several single nucleotide polymorphisms (SNPs) as compared to the gene of CL Brener BZ-susceptible strain. Here we report the sequence of TcABCG1 gene of fourteen T. cruzi strains, with diverse degrees of BZ sensitivity and belonging to different discrete typing units (DTUs) and Tcbat group. Although DTU-specific SNPs and amino acid changes were identified, no direct correlation with BZ-resistance phenotype was found. Thus, it is plausible that the transporter abundance is a determinant factor for drug resistance, as pointed out above. Sequence data were used for Bayesian phylogenies and network genealogy analysis. The network showed a high degree of reticulation suggesting genetic exchange between the parasites. TcI and TcII clades were clearly separated. Tcbat sequences were close to TcI. A fourth clade clustered TcABCG1 haplotypes of TcV, TcVI and TcIII strains, with closer proximity to TcI. Analysis of the recombination patterns indicated that hybrid strains contain haplotypes that are mosaics most likely derived by intragenic recombination of parental sequences. The data confirm that TcII and TcIII as the parentals of TcV and TcVI DTUs. Since genetic fingerprint of TcI was found in TcIII, we sustain the previously proposed "Two Hybridization model" for the origin of hybrid strains. Among the twenty best BLASTP hits in databases, orthologues of TcABCG1 transporter were found in Leishmania spp. and African trypanosomes, though their function remains undescribed.

Structure based investigation on the binding interaction of transport proteins in leishmaniasis: insights from molecular simulation

Computational structural biology and simulations help to understand the functional mechanism of transport proteins in Leishmaniasis.

The Comparative Genomics and Phylogenomics of Leishmania amazonensis Parasite

Leishmaniasis is an infectious disease caused by Leishmania species. Leishmania amazonensis is a New World Leishmania species belonging to the Mexicana complex, which is able to cause all types of leishmaniasis infections. The L. amazonensis reference strain MHOM/BR/1973/M2269 was sequenced identifying 8,802 codifying sequences (CDS), most of them of hypothetical function. Comparative analysis using six Leishmania species showed a core set of 7,016 orthologs. L. amazonensis and Leishmania mexicana share the largest number of distinct orthologs, while Leishmania braziliensis presented the largest number of inparalogs. Additionally, phylogenomic analysis confirmed the taxonomic position for L. amazonensis within the “Mexicana complex”, reinforcing understanding of the split of New and Old World Leishmania. Potential non-homologous isofunctional enzymes (NISE) were identified between L. amazonensis and Homo sapiens that could provide new drug targets for development.

Flavone-resistant Leishmania donovani overexpresses LdMRP2 transporter in the parasite and activates host MRP2 on macrophages to circumvent the flavone-mediated cell death

In parasites, ATP-binding cassette (ABC) transporters represent an important family of proteins related to drug resistance and other biological activities. Resistance of leishmanial parasites to therapeutic drugs continues to escalate in developing countries, and in many instances, it is due to overexpressed ABC efflux pumps. Progressively adapted baicalein (BLN)-resistant parasites (pB(25)R) show overexpression of a novel ABC transporter, which was classified as ABCC2 or Leishmania donovani multidrug resistance protein 2 (LdMRP2). The protein is primarily localized in the flagellar pocket region and in internal vesicles. Overexpressed LdABCC2 confers substantial BLN resistance to the parasites by rapid drug efflux. The BLN-resistant promastigotes when transformed into amastigotes in macrophage cells cannot be cured by treatment of macrophages with BLN. Amastigote resistance is concomitant with the overexpression of macrophage MRP2 transporter. Reporter analysis and site-directed mutagenesis assays demonstrated that antioxidant response element 1 is activated upon infection. The expression of this phase II detoxifying gene is regulated by NFE2-related factor 2 (Nrf2)-mediated antioxidant response element activation. In view of the fact that the signaling pathway of phosphoinositol 3-kinase controls microfilament rearrangement and translocation of actin-associated proteins, the current study correlates with the intricate pathway of phosphoinositol 3-kinase-mediated nuclear translocation of Nrf2, which activates MRP2 expression in macrophages upon infection by the parasites. In contrast, phalloidin, an agent that prevents depolymerization of actin filaments, inhibits Nrf2 translocation and Mrp2 gene activation by pB(25)R infection. Taken together, these results provide insight into the mechanisms by which resistant clinical isolates of L. donovani induce intracellular events relevant to drug resistance.

Elucidation of cellular mechanisms involved in experimental paromomycin resistance in Leishmania donovani

Leishmania donovani is the causative agent of the potentially fatal disease visceral leishmaniasis (VL). Chemotherapeutic options available to treat VL are limited and often face parasite resistance, inconsistent efficacy, and toxic side effects. Paromomycin (PMM) was recently introduced to treat VL as a monotherapy and in combination therapy. It is vital to understand the mechanisms of PMM resistance to safeguard the drug. In the present study, we utilized experimentally generated PMM-resistant L. donovani to elucidate the mechanisms of resistance and parasite biology. We found increased membrane fluidity accompanied by decreased intracellular drug accumulation in the PMM-resistant parasites. There were marked increases in gene expression of ATP-binding cassette (ABC) transporters (MDR1 and MRPA) and protein phosphatase 2A that evince increased drug efflux. Further, evaluation of parasite tolerance toward host leishmanicidal mechanisms revealed PMM-resistant parasites as being more tolerant to nitrosative stress at the promastigote and amastigote stages. The PMM-resistant parasites also predicted a better survival capacity, as indicated by resistance to complement-mediated lysis and increased stimulation of host interleukin-10 (IL-10) expression. The susceptibilities of PMM-resistant isolates to other antileishmanial agents (sodium antimony gluconate and miltefosine) remained unchanged. The data implicated the roles of altered membrane fluidity, decreased drug accumulation, increased expression of ABC transporters, and greater tolerance of parasites to host defense mechanisms in conferring PMM resistance in Leishmania.

Identification of specific reversal agents for Leishmania ABCI4-mediated antimony resistance by flavonoid and trolox derivative screening

OBJECTIVES

To identify reversal agents for the Leishmania ABCI4 transporter that confers resistance to antimony.

METHODS

Selective ABCI4 inhibitors among a series of 15 flavonoid and trolox derivatives or analogues were investigated by evaluating their ability to reverse antimony resistance in Leishmania parasites overexpressing ABCI4. Among the compounds screened, N-ethyltrolox carboxamide (compound D2) produced the highest reversal activity. In order to optimize the activity of D2, we synthesized a series of 10 derivatives by condensation of various amines with trolox.

RESULTS

Analysis of antimony resistance reversal activity showed that N-propyltrolox carboxamide (compound D4) was the most potent ABCI4 inhibitor, with reversal activity being maintained in the intracellular amastigote stage. In addition, trolox derivatives significantly reverted the resistance to zinc protoporphyrin. The mechanism of action of these active derivatives was found to be related to significant reversion of Sb(III) and zinc protoporphyrin accumulation and to a decrease in drug efflux.

CONCLUSIONS

Our findings suggest that trolox derivatives D2 and D4 could be considered to be specific reversal agents targeting the Leishmania ABCI4 transporter. The structure-activity relationship obtained in the present study highlights the importance of the size and length of the alkyl substituent linked to trolox. Furthermore, the structural data obtained provide valuable information for the further development of new, even more specific and potent Leishmania ABCI4 reversal agents.

A new ABC half-transporter in Leishmania major is involved in resistance to antimony

The characterization of ABCI4, a new intracellular ATP-binding cassette (ABC) half-transporter in Leishmania major, is described. We show that ABCI4 is involved in heavy metal export, thereby conferring resistance to Pentostam, to Sb(III), and to As(III) and Cd(II). Parasites overexpressing ABCI4 showed a lower mitochondrial toxic effect of antimony by decreasing reactive oxygen species production and maintained higher values of both the mitochondrial electrochemical potential and total ATP levels with respect to controls. The ABCI4 half-transporter forms homodimers as determined by a coimmunoprecipitation assay. A combination of subcellular localization studies under a confocal microscope and a surface biotinylation assay using parasites expressing green fluorescent protein- and FLAG-tagged ABCI4 suggests that the transporter presents a dual localization in both mitochondria and the plasma membrane. Parasites overexpressing ABCI4 present an increased replication in mouse peritoneal macrophages. We have determined that porphyrins are substrates for ABCI4. Consequently, the overexpression of ABCI4 confers resistance to some toxic porphyrins, such as zinc-protoporphyrin, due to the lower accumulation resulting from a significant efflux, as determined using the fluorescent zinc-mesoporphyrin, a validated heme analog. In addition, ABCI4 has a significant ability to efflux thiol after Sb(III) incubation, thus meaning that ABCI4 could be considered to be a potential thiol-X-pump that is able to recognize metal-conjugated thiols. In summary, we have shown that this new ABC transporter is involved in drug sensitivity to antimony and other compounds by efflux as conjugated thiol complexes.

Genome-wide identification, characterization and phylogenetic analysis of 50 catfish ATP-binding cassette (ABC) transporter genes

Background

Although a large set of full-length transcripts was recently assembled in catfish, annotation of large gene families, especially those with duplications, is still a great challenge. Most often, complexities in annotation cause mis-identification and thereby much confusion in the scientific literature. As such, detailed phylogenetic analysis and/or orthology analysis are required for annotation of genes involved in gene families. The ATP-binding cassette (ABC) transporter gene superfamily is a large gene family that encodes membrane proteins that transport a diverse set of substrates across membranes, playing important roles in protecting organisms from diverse environment.

Methodology/Principal Findings

In this work, we identified a set of 50 ABC transporters in catfish genome. Phylogenetic analysis allowed their identification and annotation into seven subfamilies, including 9 ABCA genes, 12 ABCB genes, 12 ABCC genes, 5 ABCD genes, 2 ABCE genes, 4 ABCF genes and 6 ABCG genes. Most ABC transporters are conserved among vertebrates, though cases of recent gene duplications and gene losses do exist. Gene duplications in catfish were found for ABCA1, ABCB3, ABCB6, ABCC5, ABCD3, ABCE1, ABCF2 and ABCG2.

Conclusion/Significance

The whole set of catfish ABC transporters provide the essential genomic resources for future biochemical, toxicological and physiological studies of ABC drug efflux transporters. The establishment of orthologies should allow functional inferences with the information from model species, though the function of lineage-specific genes can be distinct because of specific living environment with different selection pressure.

LABCG2, a new ABC transporter implicated in phosphatidylserine exposure, is involved in the infectivity and pathogenicity of Leishmania

Leishmaniasis is a neglected disease produced by the intracellular protozoan parasite Leishmania. In the present study, we show that LABCG2, a new ATP-binding cassette half-transporter (ABCG subfamily) from Leishmania, is involved in parasite virulence. Down-regulation of LABCG2 function upon expression of an inactive mutant version of this half-transporter (LABCG2(K/M)) is shown to reduce the translocation of short-chain analogues of phosphatidylserine (PS). This dominant-negative phenotype is specific for the headgroup of the phospholipid, as the movement of phospholipid analogues of phosphatidylcholine, phosphatidylethanolamine or sphingomyelin is not affected. In addition, promastigotes expressing LABCG2(K/M) expose less endogenous PS in the stationary phase than control parasites. Transient exposure of PS at the outer leaflet of the plasma membrane is known to be one of the mechanisms used by Leishmania to infect macrophages and to silence their immune response. Stationary phase/metacyclic promastigotes expressing LABCG2(K/M) are less infective for macrophages and show decreased pathogenesis in a mouse model of cutaneous leishmaniasis. Thus, mice infected with parasites expressing LABCG2(K/M) did not develop any lesion and showed significantly lower inflammation and parasite burden than mice infected with control parasites. Our results indicate that LABCG2 function is required for the externalization of PS in Leishmania promastigotes, a process that is involved in the virulence of the parasite.

Deletion of an ATP-binding cassette protein subfamily C transporter in Leishmania donovani results in increased virulence

Genome comparison of Leishmania species sequenced to date has identified several hundred differentially distributed genes which are present as functional genes in some species but as pseudogenes or absent in others. It is not clear whether these differentially distributed genes are important for disease, redundant or even harmful for a particular Leishmania species though a few of these genes have been implicated in Leishmania infection tissue tropism. The L. infantum LinJ.24.1510 gene is an ATP-binding cassette transporter protein subfamily C member which is present in L. donovani complex species and L. mexicana but absent in L. major and L. braziliensis and its substrate is unknown. Experimental deletion of this gene from L. donovani has resulted in an increase in growth as axenic amastigotes and increased virulence in infecting mice. To our knowledge, this is the first gene identified in Leishmania where its deletion results in hypervirulence.