Gene expression modulation is associated with gene amplification, supernumerary chromosomes and chromosome loss in antimony-resistant Leishmania infantum

Characterizing Leishmania infantum-induced resistance to trivalent stibogluconate (SbIII) through deep proteomics

Leishmania infantum belongs to the L. donovani complex, which includes species associated with visceral leishmaniasis. Traditionally, antimonial compounds have served as the primary antiparasitic treatment for all clinical forms of leishmaniasis. However, the global spread of resistance to these compounds has posed a significant challenge in the treatment in some regions. In this study, we aimed to investigate resistance to trivalent sodium stibogluconate in vitro using promastigotes from a wild strain of L. infantum. We compared the growth rates and proteomic profiles of wild-type and resistant line conducting label-free quantitative mass spectrometry-based proteomic analyses. Statistical and bioinformatics analyses were employed to evaluate the significance of protein concentration changes, protein identity annotation, GO term analysis, biosynthetic pathways, and protein-protein interactions. Our findings revealed that the resistant line displayed a notable reduction in growth rate. Proteomic data unveiled similar protein concentrations per cell in both groups but with differing molecule copy numbers. We identified 165 proteins with increased concentration, these were associated with transcription and translation activities, lipid metabolism, energy metabolism, and peroxisome biogenesis. In the decreased protein groups were 56 proteins linked to metal acquisition and metabolism, particularly iron. These results suggest a novel perspective on antimonial resistance, highlighting the importance of post-transcriptional and post-translational regulation, alongside energy expenditure compensation and alterations in organelle membrane lipid composition in antimonial-resistant parasites. Overall, our study provides insights into the proteomic profile of stibogluconate-resistant strain, contributing to our general understanding of the complex landscape of antiparasitic resistance in L. infantum. SIGNIFICANCE: Species within the Leishmania donovani complex are implicated in cases of visceral leishmaniasis in the world. Leishmania infantum is a species that predominates in regions spanning the Mediterranean Basin, the Middle East, Central Asia, South and Central America. Antimonials were the first treatment for leishmaniasis, however in the last decades, the resistance has emerged in subregions like India, where it is not a therapeutic option. In contrast, sodium stibogluconate (SbIII) remains the first-line treatment in the Americas. Unfortunately, the emergence of resistance has outpaced the development of new therapeutic options, thereby becoming a critical point in the struggle against the disease. In this study we performed an in-depth proteomic analysis with liquid chromatography mass-mass spectrometry (LC-MS/MS) on L. infantum with Sb-induced resistance in vitro. Results showed a complex proteomic adaptation in the resistant line, involving transcriptional and translational proteins, energy compensation, and homeostasis maintenance. These insights contribute to understanding the molecular adaptation in the parasite and provide information to new investigations related to therapeutics development.

Effect of Phlebotomus papatasi on the fitness, infectivity and antimony-resistance phenotype of antimony-resistant Leishmania major Mon-25

Leishmania major is responsible for zoonotic cutaneous leishmaniasis. Therapy is mainly based on the use of antimony-based drugs; however, treatment failures and illness relapses were reported. Although studies were developed to understand mechanisms of drug resistance, the interactions of resistant parasites with their reservoir hosts and vectors remain poorly understood. Here we compared the development of two L. major MON-25 trivalent antimony-resistant lines, selected by a stepwise in vitro Sb(III)-drug pressure, to their wild-type parent line in the natural vector Phlebotomus papatasi. The intensity of infection, parasite location and morphological forms were compared by microscopy. Parasite growth curves and IC50 values have been determined before and after the passage in Ph. papatasi. qPCR was used to assess the amplification rates of some antimony-resistance gene markers. In the digestive tract of sand flies, Sb(III)-resistant lines developed similar infection rates as the wild-type lines during the early-stage infections, but significant differences were observed during the late-stage of the infections. Thus, on day 7 p. i., resistant lines showed lower representation of heavy infections with colonization of the stomodeal valve and lower percentage of metacyclic promastigote forms in comparison to wild-type strains. Observed differences between both resistant lines suggest that the level of Sb(III)-resistance negatively correlates with the quality of the development in the vector. Nevertheless, both resistant lines developed mature infections with the presence of infective metacyclic forms in almost half of infected sandflies. The passage of parasites through the sand fly guts does not significantly influence their capacity to multiply in vitro. The IC50 values and molecular analysis of antimony-resistance genes showed that the resistant phenotype of Sb(III)-resistant parasites is maintained after passage through the sand fly. Sb(III)-resistant lines of L. major MON-25 were able to produce mature infections in Ph. papatasi suggesting a possible circulation in the field using this vector.

Genomic Insight of Leishmania Parasite: In-Depth Review of Drug Resistance Mechanisms and Genetic Mutations

Leishmaniasis, which is caused by a parasitic protozoan of the genus Leishmania, is still a major threat to global health, impacting millions of individuals worldwide in endemic areas. Chemotherapy has been the principal method for managing leishmaniasis; nevertheless, the evolution of drug resistance offers a significant obstacle to therapeutic success. Drug-resistant behavior in these parasites is a complex phenomenon including both innate and acquired mechanisms. Resistance is frequently related to changes in drug transportation, drug target alterations, and enhanced efflux of the drug from the pathogen. This review has revealed specific genetic mutations in Leishmania parasites that are associated with resistance to commonly used antileishmanial drugs such as pentavalent antimonials, miltefosine, amphotericin B, and paromomycin, resulting in changes in gene expression along with the functioning of various proteins involved in drug uptake, metabolism, and efflux. Understanding the genetic changes linked to drug resistance in Leishmania parasites is essential for creating approaches for tackling and avoiding the spread of drug-resistant variants. Based on which specific treatments focus on mutations and pathways could potentially improve treatment efficacy and help long-term leishmaniasis control. More study is needed to uncover the complete range of genetic changes generating medication resistance and to develop new therapies based on available information.

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.

Kinetoplast Genome of Leishmania spp. Is under Strong Purifying Selection

Instability is an intriguing characteristic of many protist genomes, and trypanosomatids are not an exception in this respect. Some regions of trypanosomatid genomes evolve fast. For instance, the trypanosomatid mitochondrial (kinetoplast) genome consists of fairly conserved maxicircle and minicircle molecules that can, nevertheless, possess high nucleotide substitution rates between closely related strains. Recent experiments have demonstrated that rapid laboratory evolution can result in the non-functionality of multiple genes of kinetoplast genomes due to the accumulation of mutations or loss of critical genomic components. An example of a loss of critical components is the reported loss of entire minicircle classes in Leishmania tarentolae during laboratory cultivation, which results in an inability to generate some correctly encoded genes. In the current work, we estimated the evolutionary rates of mitochondrial and nuclear genome regions of multiple natural Leishmania spp. We analyzed synonymous and non-synonymous substitutions and, rather unexpectedly, found that the coding regions of kinetoplast maxicircles are among the most variable regions of both genomes. In addition, we demonstrate that synonymous substitutions greatly predominate among maxicircle coding regions and that most maxicircle genes show signs of purifying selection. These results imply that maxicircles in natural Leishmania populations remain functional despite their high mutation rate.

Leishmania allelic selection during experimental sand fly infection correlates with mutational signatures of oxidative DNA damage

Trypanosomatid pathogens are transmitted by blood-feeding insects, causing devastating human infections. These parasites show important phenotypic shifts that often impact parasite pathogenicity, tissue tropism, or drug susceptibility. The evolutionary mechanisms that allow for the selection of such adaptive phenotypes remain only poorly investigated. Here, we use Leishmania donovani as a trypanosomatid model pathogen to assess parasite evolutionary adaptation during experimental sand fly infection. Comparing the genome of the parasites before and after sand fly infection revealed a strong population bottleneck effect as judged by allele frequency analysis. Apart from random genetic drift caused by the bottleneck effect, our analyses revealed haplotype and allelic changes during sand fly infection that seem under natural selection given their convergence between independent biological replicates. Our analyses further uncovered signature mutations of oxidative DNA damage in the parasite genomes after sand fly infection, suggesting that Leishmania suffers from oxidative stress inside the insect digestive tract. Our results propose a model of Leishmania genomic adaptation during sand fly infection, with oxidative DNA damage and DNA repair processes likely driving haplotype and allelic selection. The experimental and computational framework presented here provides a useful blueprint to assess evolutionary adaptation of other eukaryotic pathogens inside their insect vectors, such as Plasmodium spp, Trypanosoma brucei, and Trypanosoma cruzi.

Genomic analysis of Leishmania turanica strains from different regions of Central Asia

The evolution in Leishmania is governed by the opposite forces of clonality and sexual reproduction, with vicariance being an important factor. As such, Leishmania spp. populations may be monospecific or mixed. Leishmania turanica in Central Asia is a good model to compare these two types. In most areas, populations of L. turanica are mixed with L. gerbilli and L. major. Notably, co-infection with L. turanica in great gerbils helps L. major to withstand a break in the transmission cycle. Conversely, the populations of L. turanica in Mongolia are monospecific and geographically isolated. In this work, we compare genomes of several well-characterized strains of L. turanica originated from monospecific and mixed populations in Central Asia in order to shed light on genetic factors, which may drive evolution of these parasites in different settings. Our results illustrate that evolutionary differences between mixed and monospecific populations of L. turanica are not dramatic. On the level of large-scale genomic rearrangements, we confirmed that different genomic loci and different types of rearrangements may differentiate strains originated from mixed and monospecific populations, with genome translocations being the most prominent example. Our data suggests that L. turanica has a significantly higher level of chromosomal copy number variation between the strains compared to its sister species L. major with only one supernumerary chromosome. This suggests that L. turanica (in contrast to L. major) is in the active phase of evolutionary adaptation.

The paradigm of intracellular parasite survival and drug resistance in leishmanial parasite through genome plasticity and epigenetics: Perception and future perspective

Leishmania is an intracellular, zoonotic, kinetoplastid eukaryote with more than 1.2 million cases all over the world. The leishmanial chromosomes are divided into polymorphic chromosomal ends, conserved central domains, and antigen-encoding genes found in telomere-proximal regions. The genome flexibility of chromosomal ends of the leishmanial parasite is known to cause drug resistance and intracellular survival through the evasion of host defense mechanisms. Therefore, in this review, we discuss the plasticity of Leishmania genome organization which is the primary cause of drug resistance and parasite survival. Moreover, we have not only elucidated the causes of such genome plasticity which includes aneuploidy, epigenetic factors, copy number variation (CNV), and post-translation modification (PTM) but also highlighted their impact on drug resistance and parasite survival.

Life in plastic, it's fantastic! How Leishmania exploit genome instability to shape gene expression

Leishmania are kinetoplastid pathogens that cause leishmaniasis, a debilitating and potentially life-threatening infection if untreated. Unusually, Leishmania regulate their gene expression largely post-transcriptionally due to the arrangement of their coding genes into polycistronic transcription units that may contain 100s of functionally unrelated genes. Yet, Leishmania are capable of rapid and responsive changes in gene expression to challenging environments, often instead correlating with dynamic changes in their genome composition, ranging from chromosome and gene copy number variations to the generation of extrachromosomal DNA and the accumulation of point mutations. Typically, such events indicate genome instability in other eukaryotes, coinciding with genetic abnormalities, but for Leishmania, exploiting these products of genome instability can provide selectable substrates to catalyse necessary gene expression changes by modifying gene copy number. Unorthodox DNA replication, DNA repair, replication stress factors and DNA repeats are recognised in Leishmania as contributors to this intrinsic instability, but how Leishmania regulate genome plasticity to enhance fitness whilst limiting toxic under- or over-expression of co-amplified and co-transcribed genes is unclear. Herein, we focus on fresh, and detailed insights that improve our understanding of genome plasticity in Leishmania. Furthermore, we discuss emerging models and factors that potentially circumvent regulatory issues arising from polycistronic transcription. Lastly, we highlight key gaps in our understanding of Leishmania genome plasticity and discuss future studies to define, in higher resolution, these complex regulatory interactions.

Self-Hybridization in Leishmania major

Genetic exchange between different Leishmania strains in the sand fly vector has been experimentally demonstrated and is supported by population genetic studies. In nature, opportunities for Leishmania interstrain mating are restricted to flies biting multiply infected hosts or through multiple bites of different hosts. In contrast, self-mating could occur in any infected sand fly. By crossing two recombinant lines derived from the same Leishmania major strain, each expressing a different drug-resistance marker, self-hybridization in L. major was confirmed in a natural sand fly vector, Phlebotomus duboscqi, and in frequencies comparable to interstrain crosses. We provide the first high resolution, whole-genome sequencing analysis of large numbers of selfing progeny, their parents, and parental subclones. Genetic exchange consistent with classical meiosis is supported by the biallelic inheritance of the rare homozygous single nucleotide polymorphisms (SNPs) that arose by mutation during the generation of the parental clones. In contrast, heterozygous SNPs largely failed to be transmitted in Mendelian ratios for reasons not understood. SNPs that were heterozygous in both parents, however, recombined to produce homozygous alleles in some hybrids. For trisomic chromosomes present in both parents, transmittal to the progeny was only altered by self-hybridization, involving a gain or loss of somy in frequencies predicted by a meiotic process. Whole-genome polyploidization was also observed in the selfing progeny. Thus, self-hybridization in Leishmania, with its potential to occur in any infected sand fly, may be an important source of karyotype variation, loss of heterozygosity, and functional diversity. IMPORTANCE Leishmania are parasitic protozoa that cause a wide spectrum of diseases collectively known as the leishmaniases. Sexual reproduction in Leishmania has been proposed as an important source of genetic diversity and has been formally demonstrated to occur inside the sand fly vector midgut. Nevertheless, in the wild, opportunities for genetic exchange between different Leishmania species or strains are restricted by the capacity of different Leishmania strains to colonize the same sand fly. In this work, we report the first high resolution, whole-genome sequence analysis of intraclonal genetic exchange as a type of self-mating in Leishmania. Our data reveal that self-hybridization can occur with comparable frequency as interstrain mating under experimental lab conditions, leading to important genomic alterations that can potentially take place within every naturally infected sand fly.

Investigating the Leishmania donovani sacp Gene and Its Role in Macrophage Infection and Survival in Mice

The protozoan parasite Leishmania donovani is a causative agent of the neglected tropical disease known as visceral leishmaniasis, which can be lethal when untreated. Studying Leishmania viru-lence factors is crucial in determining how the parasite causes disease and identifying new targets for treatment. One potential virulence factor is L. donovani's abundantly secreted protein: secreted acid phosphatase (SAcP). Whole-genome analysis revealed that the sacp gene was present in three copies in wild type L. donovani. Using CRISPR-Cas9 gene editing; we generated a sacp gene knockout termed LdΔSAcP, which demonstrated a loss of both the SAcP protein and an associated reduction in secreted acid phosphatase activity. Genome sequencing confirmed the precise dele-tion of the sacp gene in LdΔSAcP and identified several changes in the genome. LdΔSAcP demonstrated no significant changes in promastigote proliferation or its ability to infect and survive in macrophages compared to the wildtype strain. LdΔSAcP also demonstrated no change in murine liver infection; however, survival was impaired in the spleen. Taken together these results show that SAcP is not necessary for the survival of promastigotes in culture but may support long-term survival in the spleen. These observations also show that the use of CRISPR gene editing and WGS together are effective to investigate the function and phenotype of complex potential drug targets such as multicopy genes.

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.

Exploring direct and indirect targets of current antileishmanial drugs using a novel thermal proteomics profiling approach

Visceral leishmaniasis (VL), caused by Leishmania infantum, is an oft-fatal neglected tropical disease. In the absence of an effective vaccine, the control of leishmaniasis relies exclusively on chemotherapy. Due to the lack of established molecular/genetic markers denoting parasite resistance, clinical treatment failure is often used as an indicator. Antimony-based drugs have been the standard antileishmanial treatment for more than seven decades, leading to major drug resistance in certain regions. Likewise, drug resistance to miltefosine and amphotericin B continues to spread at alarming rates. In consequence, innovative approaches are needed to accelerate the identification of antimicrobial drug targets and resistance mechanisms. To this end, we have implemented a novel approach based on thermal proteome profiling (TPP) to further characterize the mode of action of antileishmanials antimony, miltefosine and amphotericin B, as well as to better understand the mechanisms of drug resistance deployed by Leishmania. Proteins become more resistant to heat-induced denaturation when complexed with a ligand. In this way, we used multiplexed quantitative mass spectrometry-based proteomics to monitor the melting profile of thousands of expressed soluble proteins in WT, antimony-resistant, miltefosine-resistant, and amphotericin B-resistant L. infantum parasites, in the presence (or absence) of the above-mentioned drugs. Bioinformatics analyses were performed, including data normalization, melting profile fitting, and identification of proteins that underwent changes (fold change > 4) caused by complexation with a drug. With this unique approach, we were able to narrow down the regions of the L. infantum proteome that interact with antimony, miltefosine, and amphotericin B; validating previously-identified and unveiling novel drug targets. Moreover, analyses revealed candidate proteins potentially involved in drug resistance. Interestingly, we detected thermal proximity coaggregation for several proteins belonging to the same metabolic pathway (i.e., tryparedoxin peroxidase and aspartate aminotransferase in proteins exposed to antimony), highlighting the importance of these pathways. Collectively, our results could serve as a jumping-off point for the future development of innovative diagnostic tools for the detection and evaluation of antimicrobial-resistant Leishmania populations, as well as open the door for new on-target therapies.

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.

Experimental evolution links post-transcriptional regulation to Leishmania fitness gain

The protozoan parasite Leishmania donovani causes fatal human visceral leishmaniasis in absence of treatment. Genome instability has been recognized as a driver in Leishmania fitness gain in response to environmental change or chemotherapy. How genome instability generates beneficial phenotypes despite potential deleterious gene dosage effects is unknown. Here we address this important open question applying experimental evolution and integrative systems approaches on parasites adapting to in vitro culture. Phenotypic analyses of parasites from early and late stages of culture adaptation revealed an important fitness tradeoff, with selection for accelerated growth in promastigote culture (fitness gain) impairing infectivity (fitness costs). Comparative genomics, transcriptomics and proteomics analyses revealed a complex regulatory network associated with parasite fitness gain, with genome instability causing highly reproducible, gene dosage-independent and -dependent changes. Reduction of flagellar transcripts and increase in coding and non-coding RNAs implicated in ribosomal biogenesis and protein translation were not correlated to dosage changes of the corresponding genes, revealing a gene dosage-independent, post-transcriptional mechanism of regulation. In contrast, abundance of gene products implicated in post-transcriptional regulation itself correlated to corresponding gene dosage changes. Thus, RNA abundance during parasite adaptation is controled by direct and indirect gene dosage changes. We correlated differential expression of small nucleolar RNAs (snoRNAs) with changes in rRNA modification, providing first evidence that Leishmania fitness gain in culture may be controlled by post-transcriptional and epitranscriptomic regulation. Our findings propose a novel model for Leishmania fitness gain in culture, where differential regulation of mRNA stability and the generation of modified ribosomes may potentially filter deleterious from beneficial gene dosage effects and provide proteomic robustness to genetically heterogenous, adapting parasite populations. This model challenges the current, genome-centric approach to Leishmania epidemiology and identifies the Leishmania transcriptome and non-coding small RNome as potential novel sources for the discovery of biomarkers that may be associated with parasite phenotypic adaptation in clinical settings.

Genome instability plays a central yet poorly understood role in human disease. Gene amplifications and deletions drive cancer development, microbial infection and therapeutic failure. The molecular mechanisms that harness the deleterious effects of genome instability to generate beneficial phenotypes in pathogenic systems are unknown. Here we study this important open question in the protozoan parasite Leishmania that causes devastating human diseases termed leishmaniases. Leishmania parasites lack transcriptional control and instead exploit genome instability to adapt to their host environment. Analyzing in vitro adaptation of hamster-derived parasites via gene copy number (genomic level) and gene expression changes (transcriptomic and proteomic levels), we show that these parasites likely exploit small nucleolar RNAs (snoRNAs) to mitigate toxic effects of genome instability by post-transcriptional regulation and the establishment of modified ribosomes. Our findings propose non-coding RNAs as potential novel biomarkers with diagnostic and prognostic value that may be linked to changes in parasite tissue tropism or drug susceptibility. This novel insight into Leishmania adaptation will be likely applicable to other fast evolving eukaryotic systems with unstable genomes, such as fungi or cancer cells.

Simvastatin Resistance of Leishmania amazonensis Induces Sterol Remodeling and Cross-Resistance to Sterol Pathway and Serine Protease Inhibitors

The sterol biosynthesis pathway of Leishmania spp. is used as a pharmacological target; however, available information about the mechanisms of the regulation and remodeling of sterol-related genes is scarce. The present study investigated compensatory mechanisms of the sterol biosynthesis pathway using an inhibitor of HMG-CoA reductase (simvastatin) and by developing drug-resistant parasites to evaluate the impact on sterol remodeling, cross-resistance, and gene expression. Simvastatin-resistant L. amazonensis parasites (LaSimR) underwent reprogramming of sterol metabolism manifested as an increase in cholestane- and stigmastane-based sterols and a decrease in ergostane-based sterols. The levels of the transcripts of sterol 24-C-methyltransferase (SMT), sterol C14-α-demethylase (C14DM), and protease subtilisin (SUB) were increased in LaSimR. LaSimR was cross-resistance to ketoconazole (a C14DM inhibitor) and remained sensitive to terbinafine (an inhibitor of squalene monooxygenase). Sensitivity of the LaSimR mutant to other antileishmanial drugs unrelated to the sterol biosynthesis pathway, such as trivalent antimony and pentamidine, was similar to that of the wild-type strain; however, LaSimR was cross-resistant to miltefosine, general serine protease inhibitor N-p-tosyl-l-phenylalanine chloromethyl ketone (TPCK), subtilisin-specific inhibitor 4-[(diethylamino)methyl]-N-[2-(2-methoxyphenyl)ethyl]-N-(3R)-3-pyrrolidinyl-benzamide dihydrochloride (PF-429242), and tunicamycin. The findings on the regulation of the sterol pathway can support the development of drugs and protease inhibitors targeting this route in parasites.

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

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

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

Background

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

Methodology/principal findings

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

Conclusion/significance

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

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

Genome instability drives epistatic adaptation in the human pathogen Leishmania

Chromosome and gene copy number variations often correlate with the evolution of microbial and cancer drug resistance, thus causing important human mortality. How genome instability is harnessed to generate beneficial phenotypes and how deleterious gene dosage effects are compensated remain open questions. The protist pathogen Leishmania exploits genome instability to regulate expression via gene dosage changes. Using these parasites as a unique model system, we uncover complex epistatic interactions between gene copy number variations and compensatory transcriptomic responses as key processes that harness genome instability for adaptive evolution in Leishmania. Our results propose a model of eukaryotic fitness gain that may be broadly applicable to pathogenic fungi or tumor cells known to exploit genome instability for adaptation.

How genome instability is harnessed for fitness gain despite its potential deleterious effects is largely elusive. An ideal system to address this important open question is provided by the protozoan pathogen Leishmania, which exploits frequent variations in chromosome and gene copy number to regulate expression levels. Using ecological genomics and experimental evolution approaches, we provide evidence that Leishmania adaptation relies on epistatic interactions between functionally associated gene copy number variations in pathways driving fitness gain in a given environment. We further uncover posttranscriptional regulation as a key mechanism that compensates for deleterious gene dosage effects and provides phenotypic robustness to genetically heterogenous parasite populations. Finally, we correlate dynamic variations in small nucleolar RNA (snoRNA) gene dosage with changes in ribosomal RNA 2′-O-methylation and pseudouridylation, suggesting translational control as an additional layer of parasite adaptation. Leishmania genome instability is thus harnessed for fitness gain by genome-dependent variations in gene expression and genome-independent compensatory mechanisms. This allows for polyclonal adaptation and maintenance of genetic heterogeneity despite strong selective pressure. The epistatic adaptation described here needs to be considered in Leishmania epidemiology and biomarker discovery and may be relevant to other fast-evolving eukaryotic cells that exploit genome instability for adaptation, such as fungal pathogens or cancer.

Unpicking the Roles of DNA Damage Protein Kinases in Trypanosomatids

To preserve genome integrity when faced with DNA lesions, cells activate and coordinate a multitude of DNA repair pathways to ensure timely error correction or tolerance, collectively called the DNA damage response (DDR). These interconnecting damage response pathways are molecular signal relays, with protein kinases (PKs) at the pinnacle. Focused efforts in model eukaryotes have revealed intricate aspects of DNA repair PK function, including how they direct DDR pathways and how repair reactions connect to wider cellular processes, including DNA replication and transcription. The Kinetoplastidae, including many parasites like Trypanosoma spp. and Leishmania spp. (causative agents of debilitating, neglected tropical infections), exhibit peculiarities in several core biological processes, including the predominance of multigenic transcription and the streamlining or repurposing of DNA repair pathways, such as the loss of non-homologous end joining and novel operation of nucleotide excision repair (NER). Very recent studies have implicated ATR and ATM kinases in the DDR of kinetoplastid parasites, whereas DNA-dependent protein kinase (DNA-PKcs) displays uncertain conservation, questioning what functions it fulfills. The wide range of genetic manipulation approaches in these organisms presents an opportunity to investigate DNA repair kinase roles in kinetoplastids and to ask if further kinases are involved. Furthermore, the availability of kinase inhibitory compounds, targeting numerous eukaryotic PKs, could allow us to test the suitability of DNA repair PKs as novel chemotherapeutic targets. Here, we will review recent advances in the study of trypanosomatid DNA repair kinases.

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.

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.

Extracellular vesicles and leishmaniasis: Current knowledge and promising avenues for future development

Extracellular vesicles (EVs) are small, membrane-bound "delivery trucks" that are present in the extracellular environment, including biological fluids. EVs are capable of inducing changes in the physiological status of neighboring cells through the transfer of key macromolecules, and are thought to play a role in a number of pathological processes. Leishmaniasis, caused by the protozoan parasite Leishmania, is an important example. The biology of Leishmania EVs has been studied in detail, and findings point to their role in exacerbation of disease and potential involvement in the perpetuation of drug resistance. Furthermore, the use of EVs for development of vaccines has been explored, as well as their potential use in a number of fields as biomarkers of disease and drug resistance. Here we discuss the latest findings on EVs, with a particular focus on Leishmania, as well as potential avenues for their future development and clinical applications.

Ros3 (Lem3p/CDC50) Gene Dosage Is Implicated in Miltefosine Susceptibility in Leishmania (Viannia) braziliensis Clinical Isolates and in Leishmania (Leishmania) major

The Ros3 protein is a component of the MT-Ros3 transporter complex, considered as the main route of miltefosine entry in Leishmania. L. braziliensis clinical isolates presenting differences in miltefosine susceptibility and uptake were previously shown to differentially express ros3. In this work, we showed that the ros3 gene copy number was increased in the isolate presenting the highest rates of miltefosine uptake and, thus, the highest susceptibility to this drug. The role of the ros3 gene dosage in miltefosine susceptibility was then investigated through a modulation of the gene copy number using two distinct approaches: through an overexpression of ros3 in a tolerant L. braziliensis clinical isolate and in L. major and by generating mono- and diallelic knockouts of this gene in L. major using clustered regularly interspaced short palindromic repeats (CRISPR) Cas9 (Cas = CRISPR-associated). Although the levels of ros3 mRNA were increased at least 40-fold in overexpressing clones, no significant reduction in the half-maximal effective concentration (EC₅₀) for miltefosine was observed in these parasites. The partial or complete deletion of ros3 in L. major, in turn, resulted in a significant increase of 3 and 20 times, respectively, in the EC₅₀ to miltefosine. We unequivocally showed that the ros3 copy number is one of the factors involved in the differential susceptibility and uptake of miltefosine.

Metabolic selection of a homologous recombination-mediated gene loss protects Trypanosoma brucei from ROS production by glycosomal fumarate reductase

The genome of trypanosomatids rearranges by using repeated sequences as platforms for amplification or deletion of genomic segments. These stochastic recombination events have a direct impact on gene dosage and foster the selection of adaptive traits in response to environmental pressure. We provide here such an example by showing that the phosphoenolpyruvate carboxykinase (PEPCK) gene knockout (Δpepck) leads to the selection of a deletion event between two tandemly arranged fumarate reductase (FRDg and FRDm2) genes to produce a chimeric FRDg-m2 gene in the Δpepck∗ cell line. FRDg is expressed in peroxisome-related organelles, named glycosomes, expression of FRDm2 has not been detected to date, and FRDg-m2 is nonfunctional and cytosolic. Re-expression of FRDg significantly impaired growth of the Δpepck∗ cells, but FRD enzyme activity was not required for this negative effect. Instead, glycosomal localization as well as the covalent flavinylation motif of FRD is required to confer growth retardation and intracellular accumulation of reactive oxygen species (ROS). The data suggest that FRDg, similar to Escherichia coli FRD, can generate ROS in a flavin-dependent process by transfer of electrons from NADH to molecular oxygen instead of fumarate when the latter is unavailable, as in the Δpepck background. Hence, growth retardation is interpreted as a consequence of increased production of ROS, and rearrangement of the FRD locus liberates Δpepck∗ cells from this obstacle. Interestingly, intracellular production of ROS has been shown to be required to complete the parasitic cycle in the insect vector, suggesting that FRDg may play a role in this process.

Heat Shock Proteins as the Druggable Targets in Leishmaniasis: Promises and Perils

Leishmania, the causative agent of leishmaniasis, is an intracellular pathogen that thrives in the insect gut and mammalian macrophages to complete its life cycle. Apart from temperature difference (26 to 37°C), it encounters several harsh conditions, including oxidative stress, inflammatory reactions, and low pH. Heat shock proteins (HSPs) play essential roles in cell survival by strategically reprogramming cellular processes and signaling pathways. HSPs assist cells in multiple functions, including differentiation, adaptation, virulence, and persistence in the host cell. Due to cyclical epidemiological patterns, limited chemotherapeutic options, drug resistance, and the absence of a vaccine, control of leishmaniasis remains a far-fetched dream. The essential roles of HSPs in parasitic differentiation and virulence and increased expression in drug-resistant strains highlight their importance in combating the disease. In this review, we highlighted the diverse physiological importance of HSPs present in Leishmania, emphasizing their significance in disease pathogenesis. Subsequently, we assessed the potential of HSPs as a chemotherapeutic target and underlined the challenges associated with it. Furthermore, we have summarized a few ongoing drug discovery initiatives that need to be explored further to develop clinically successful chemotherapeutic agents in the future.

Leishmania: Responding to environmental signals and challenges without regulated transcription

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

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.

Nuclear and mitochondrial genome sequencing of North-African Leishmania infantum isolates from cured and relapsed visceral leishmaniasis patients reveals variations correlating with geography and phenotype

Although several studies have investigated genetic diversity of Leishmania infantum in North Africa, genome-wide analyses are lacking. Here, we conducted comparative analyses of nuclear and mitochondrial genomes of seven L. infantum isolates from Tunisia with the aim to gain insight into factors that drive genomic and phenotypic adaptation. Isolates were from cured (n=4) and recurrent (n=3) visceral leishmaniasis (VL) cases, originating from northern (n=2) and central (n=5) Tunisia, where respectively stable and emerging VL foci are observed. All isolates from relapsed patients were from Kairouan governorate (Centre); one showing resistance to the anti-leishmanial drug Meglumine antimoniate. Nuclear genome diversity of the isolates was analysed by comparison to the L. infantum JPCM5 reference genome. Kinetoplast maxi and minicircle sequences (1 and 59, respectively) were extracted from unmapped reads and identified by blast analysis against public data sets. The genome variation analysis grouped together isolates from the same geographical origins. Strains from the North were very different from the reference showing more than 34 587 specific single nucleotide variants, with one isolate representing a full genetic hybrid as judged by variant frequency. Composition of minicircle classes within isolates corroborated this geographical population structure. Read depth analysis revealed several significant gene copy number variations correlating with either geographical origin (amastin and Hsp33 genes) or relapse (CLN3 gene). However, no specific gene copy number variation was found in the drug-resistant isolate. In contrast, resistance was associated with a specific minicircle pattern suggesting Leishmania mitochondrial DNA as a potential novel source for biomarker discovery.

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

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

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.

Experimental Strategies to Explore Drug Action and Resistance in Kinetoplastid Parasites

Kinetoplastids are the causative agents of leishmaniasis, human African trypanosomiasis, and American trypanosomiasis. They are responsible for high mortality and morbidity in (sub)tropical regions. Adequate treatment options are limited and have several drawbacks, such as toxicity, need for parenteral administration, and occurrence of treatment failure and drug resistance. Therefore, there is an urgency for the development of new drugs. Phenotypic screening already allowed the identification of promising new chemical entities with anti-kinetoplastid activity potential, but knowledge on their mode-of-action (MoA) is lacking due to the generally applied whole-cell based approach. However, identification of the drug target is essential to steer further drug discovery and development. Multiple complementary techniques have indeed been used for MoA elucidation. In this review, the different 'omics' approaches employed to define the MoA or mode-of-resistance of current reference drugs and some new anti-kinetoplastid compounds are discussed.

MRPA-independent mechanisms of antimony resistance in Leishmania infantum

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

Global genome diversity of the Leishmania donovani complex

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

Genomic analyses reveal moderate levels of ploidy, high heterozygosity and structural variations in a Colombian isolate of Leishmania (Leishmania) amazonensis

Leishmania amazonensis is one of the causative agents of the different forms of cutaneous leishmaniasis present in Latin America. This species has been isolated from humans and animals (canine/feline) in some endemic regions of Colombia. Therefore, L. amazonensis is of great relevance at the clinical and epidemiological levels in medicine and veterinary science. Until now, very few genomes from this species are available. Here, we report the complete genome sequence of a laboratory-adapted L. amazonensis strain isolated from a human patient with clinical manifestation of cutaneous leishmaniasis in Colombia. The genome sequence not only allowed inter and intra-species comparative analyses to be performed with the sequenced genomes of L. amazonensis strains from different geographical regions, but also increased our knowledge about the genomic behavior of this L. amazonensis Colombian strain. This isolate was also characterized in terms of single nucleotide polymorphisms, chromosome and gene copy number variations (CNVs). The results revealed moderate aneuploidy, CNVs in genes involved in the virulence, growth, and survival of the parasite, and in the distributions of the multicopy genes on some chromosomes, as well as a high level of heterozygosity. The data confirmed previous reports that identified unique variations in L. amazonensis, suggesting aneuploidy may not have a high fitness cost and may allow the rapid generation of diversity in Leishmania parasites growing normally.

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.

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.

Innovative Solutions for the Control of Leishmaniases: Nanoscale Drug Delivery Systems

Background:

Leishmania are sandfly-transmitted protozoan parasites that harbour within the macrophages of a mammalian host and cause leishmaniasis, a serious zoonotic disease that threatens the lives of millions worldwide. Its numerous forms (cutaneous, mucocutaneous, and visceral) are currently treated with a sparse arsenal of drugs, specifically antimonials, amphotericin B, miltefosine, and paromomycin, for which drug resistance and clinical failure are rampant. Medicine is presently trending towards nanotechnology to aid in the successful delivery of drugs. Vehicles such as lipid-based nanocarriers, polymer-based nanoparticles, and metal ions and oxides have been previously demonstrated to improve bioavailability of drugs and decrease toxicity for the patient. These cutting-edge solutions can be combined with existing active molecules, as well as novel drugs or plant extracts with promising antileishmanial activity.

Conclusion:

This review explores the current evidence for the treatment of leishmaniases using nanoscale drug delivery systems (specifically lipid-, polymer- and metal-based systems) and encourages further development of the aforementioned nanotechnologies for treatment of Leishmania.

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

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

Genome and transcriptome analyses of Leishmania spp.: opening Pandora's box

In the last 30 years, significant advances in genetic manipulation tools along with complete genome and transcriptome sequencing have advanced our understanding of the biology of Leishmania parasites and their interplay with the sand fly and mammalian hosts. High-throughput sequencing in association with CRISPR/Cas9 have prepared the ground for significant advances. Given the richness of the progress made over the last decade, in this article, we focused on the most recent contributions of genome-wide and transcriptome analyses of Leishmania spp., which permit the comparison of life cycle stages, the evaluation of different strains and species in their natural niches and in the field and the simultaneously comparison of the gene expression profiles of parasites and hosts.

Meta-analysis and discussion on challenges to translate Leishmania drug resistance phenotyping into the clinic

Antimicrobial resistance (AMR) threatens the prevention and treatment of infections caused by a large range of microorganisms. Leishmania is not an exception and treatment failure due to drug-resistant organisms is increasingly reported. Currently, no molecular methods and marker are validated to track drug-resistant organism and antimicrobial susceptibility tests are roughly not amenable to a clinical setting. Taking these facts into account, it is essential to reflect on ways to translate basic knowledge into methodologies aimed to diagnose leishmania drug resistance. As a matter of fact, a meta-analysis of the literature discloses the reliability of the promastigotes antimicrobial susceptibility tests (AST) to predict intracellular amastigotes susceptibility status. Promastigote cultures that are easy to perform, typically inexpensive and amenable to standardization should represent a candidate to diagnose resistance. Using AST performed on promastigote, we propose a way to improve leishmania drug resistance diagnosis in the framework of guidance and guideline of the bacterial drug resistance diagnosis. In this review, we highlight challenges that remained and discuss the definition of clinical breakpoints, including the epidemiological cutoff (ECOFF), to track drug-resistant isolates. Our analysis paves the ways to standardize and analyze anti-leishmania susceptibility tests output in order to guide the characterization of drug-resistant isolates, the clinical decision during treatment and the search for new molecular markers.

Genomic instability at the locus of sterol C24-methyltransferase promotes amphotericin B resistance in Leishmania parasites

Amphotericin B is an increasingly important tool in efforts to reduce the global disease burden posed by Leishmania parasites. With few other chemotherapeutic options available for the treatment of leishmaniasis, the potential for emergent resistance to this drug is a considerable threat. Here we characterised four novel amphotericin B-resistant Leishmania mexicana lines. All lines exhibited altered sterol biosynthesis, and hypersensitivity to pentamidine. Whole genome sequencing demonstrated resistance-associated mutation of the sterol biosynthesis gene sterol C5-desaturase in one line. However, in three out of four lines, RNA-seq revealed loss of expression of sterol C24-methyltransferase (SMT) responsible for drug resistance and altered sterol biosynthesis. Additional loss of the miltefosine transporter was associated with one of those lines. SMT is encoded by two tandem gene copies, which we found to have very different expression levels. In all cases, reduced overall expression was associated with loss of the 3' untranslated region of the dominant gene copy, resulting from structural variations at this locus. Local regions of sequence homology, between the gene copies themselves, and also due to the presence of SIDER1 retrotransposon elements that promote multi-gene amplification, correlate to these structural variations. Moreover, in at least one case loss of SMT expression was not associated with loss of virulence in primary macrophages or in vivo. Whilst such repeat sequence-mediated instability is known in Leishmania genomes, its presence associated with resistance to a major antileishmanial drug, with no evidence of associated fitness costs, is a significant concern.

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.

Genomic Analysis of Colombian Leishmania panamensis strains with different level of virulence

The establishment of Leishmania infection in mammalian hosts and the subsequent manifestation of clinical symptoms require internalization into macrophages, immune evasion and parasite survival and replication. Although many of the genes involved in these processes have been described, the genetic and genomic variability associated to differences in virulence is largely unknown. Here we present the genomic variation of four Leishmania (Viannia) panamensis strains exhibiting different levels of virulence in BALB/c mice and its application to predict novel genes related to virulence. De novo DNA sequencing and assembly of the most virulent strain allowed comparative genomics analysis with sequenced L. (Viannia) panamensis and L. (Viannia) braziliensis strains, and showed important variations at intra and interspecific levels. Moreover, the mutation detection and a CNV search revealed both base and structural genomic variation within the species. Interestingly, we found differences in the copy number and protein diversity of some genes previously related to virulence. Several machine-learning approaches were applied to combine previous knowledge with features derived from genomic variation and predict a curated set of 66 novel genes related to virulence. These genes can be prioritized for validation experiments and could potentially become promising drug and immune targets for the development of novel prophylactic and therapeutic interventions.

Whole genome sequencing of Trypanosoma cruzi field isolates reveals extensive genomic variability and complex aneuploidy patterns within TcII DTU

BACKGROUND

Trypanosoma cruzi, the etiologic agent of Chagas disease, is currently divided into six discrete typing units (DTUs), named TcI-TcVI. TcII is among the major DTUs enrolled in human infections in South America southern cone, where it is associated with severe cardiac and digestive symptoms. Despite the importance of TcII in Chagas disease epidemiology and pathology, so far, no genome-wide comparisons of the mitochondrial and nuclear genomes of TcII field isolates have been performed to track the variability and evolution of this DTU in endemic regions.

RESULTS

In the present work, we have sequenced and compared the whole nuclear and mitochondrial genomes of seven TcII strains isolated from chagasic patients from the central and northeastern regions of Minas Gerais, Brazil, revealing an extensive genetic variability within this DTU. A comparison of the phylogeny based on the nuclear or mitochondrial genomes revealed that the majority of branches were shared by both sequences. The subtle divergences in the branches are probably consequence of mitochondrial introgression events between TcII strains. Two T. cruzi strains isolated from patients living in the central region of Minas Gerais, S15 and S162a, were clustered in the nuclear and mitochondrial phylogeny analysis. These two strains were isolated from the other five by the Espinhaço Mountains, a geographic barrier that could have restricted the traffic of insect vectors during T. cruzi evolution in the Minas Gerais state. Finally, the presence of aneuploidies was evaluated, revealing that all seven TcII strains have a different pattern of chromosomal duplication/loss.

CONCLUSIONS

Analysis of genomic variability and aneuploidies suggests that there is significant genomic variability within Minas Gerais TcII strains, which could be exploited by the parasite to allow rapid selection of favorable phenotypes. Also, the aneuploidy patterns vary among T. cruzi strains and does not correlate with the nuclear phylogeny, suggesting that chromosomal duplication/loss are recent and frequent events in the parasite evolution.

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

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

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

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

Selective human inhibitors of ATR and ATM render Leishmania major promastigotes sensitive to oxidative damage

All cellular processes, including those involved in normal cell metabolism to those responsible for cell proliferation or death, are finely controlled by cell signaling pathways, whose core proteins constitute the family of phosphatidylinositol 3-kinase-related kinases (PIKKs). Ataxia Telangiectasia Mutated (ATM) and Ataxia Telangiectasia and Rad3 related (ATR) are two important PIKK proteins that act in response to DNA damage, phosphorylating a large number of proteins to exert control over genomic integrity. The genus Leishmania belongs to a group of early divergent eukaryotes in evolution and has a highly plastic genome, probably owing to the existence of signaling pathways designed to maintain genomic integrity. The objective of this study was to evaluate the use of specific human inhibitors of ATR and ATM in Leishmania major. Bioinformatic analyses revealed the existence of the putative PIKK genes ATR and ATM, in addition to mTOR and DNA-PKcs in Leishmania spp. Moreover, it was possible to suggest that the inhibitors VE-821 and KU-55933 have binding affinity for the catalytic sites of putative L. major ATR and ATM, respectively. Promastigotes of L. major exposed to these inhibitors show slight growth impairment and minor changes in cell cycle and morphology. It is noteworthy that treatment of promastigotes with inhibitors VE-821 and KU-55933 enhanced the oxidative damage caused by hydrogen peroxide. These inhibitors could significantly reduce the number of surviving L. major cells following H₂O₂ exposure whilst also decreasing their evaluated IC₅₀ to H₂O₂ to less than half of that observed for non-treated cells. These results suggest that the use of specific inhibitors of ATR and ATM in Leishmania interferes in the signaling pathways of this parasite, which can impair its tolerance to DNA damage and affect its genome integrity. ATR and ATM could constitute novel targets for drug development and/or repositioning for treatment of leishmaniases.