Drug resistance analysis by next generation sequencing in Leishmania

CRISPR-Cas9 high-throughput screening to study drug resistance in Leishmania infantum

Parasites of the genus Leishmania pose a global health threat with limited treatment options. New drugs are urgently needed, and genomic screens have the potential to accelerate target discovery, mode of action, and resistance mechanisms against these new drugs. We describe here our effort in developing a genome-wide CRISPR-Cas9 screen in Leishmania, an organism lacking a functional nonhomologous end joining system that must rely on microhomology-mediated end joining, single-strand annealing, or homologous recombination for repairing Cas9-induced double-stranded DNA breaks. A new vector for cloning and expressing single guide RNAs (sgRNAs) was designed and proven to be effective in a small pilot project while enriching specific sgRNAs during drug selection. We then developed a whole-genome library of 49,754 sgRNAs, targeting all the genes of Leishmania infantum. This library was transfected in L. infantum expressing Cas9, and these cells were selected for resistance to two antileishmanials, miltefosine and amphotericin B. The sgRNAs the most enriched in the miltefosine screen targeted the miltefosine transporter gene, but sgRNAs targeting genes coding for a RING-variant protein and a transmembrane protein were also enriched. The sgRNAs the most enriched by amphotericin B targeted the sterol 24 C methyltransferase genes and a hypothetical gene. Through gene disruption experiments, we proved that loss of function of these genes was associated with resistance. This study describes the feasibility of carrying out whole-genome CRISPR-Cas9 screens in Leishmania provided that a strong selective pressure is applied. Such a screen can be used for accelerating the development of urgently needed antileishmanial drugs.IMPORTANCELeishmaniasis, a global health threat, lacks adequate treatment options and drug resistance exacerbates the challenge. This study introduces a CRISPR-Cas9 screening approach in Leishmania infantum, unraveling mechanisms of drug resistance at a genome-wide scale. Our screen was applied against two main antileishmanial drugs, and guides were enriched upon drug selection. These guides targeted known and new targets, hence validating the use of this screen against Leishmania. This strategy provides a powerful tool to expedite drug discovery as well as potential therapeutic targets against this neglected tropical disease.

Extracellular Vesicles: Translational Agenda Questions for Three Protozoan Parasites

The protozoan parasites Plasmodium falciparum, Leishmania spp. and Trypanosoma cruzi continue to exert a significant toll on the disease landscape of the human population in sub-Saharan Africa and Latin America. Control measures have helped reduce the burden of their respective diseases-malaria, leishmaniasis and Chagas disease-in endemic regions. However, the need for new drugs, innovative vaccination strategies and molecular markers of disease severity and outcomes has emerged because of developing antimicrobial drug resistance, comparatively inadequate or absent vaccines, and a lack of trustworthy markers of morbid outcomes. Extracellular vesicles (EVs) have been widely reported to play a role in the biology and pathogenicity of P. falciparum, Leishmania spp. and T. cruzi ever since they were discovered. EVs are secreted by a yet to be fully understood mechanism in protozoans into the extracellular milieu and carry a cargo of diverse molecules that reflect the originator cell's metabolic state. Although our understanding of the biogenesis and function of EVs continues to deepen, the question of how EVs in P. falciparum, Leishmania spp. and T. cruzi can serve as targets for a translational agenda into clinical and public health interventions is yet to be fully explored. Here, as a consortium of protozoan researchers, we outline a plan for future researchers and pose three questions to direct an EV's translational agenda in P. falciparum, Leishmania spp. and T. cruzi. We opine that in the long term, executing this blueprint will help bridge the current unmet needs of these medically important protozoan diseases in sub-Saharan Africa and Latin America.

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.

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

BACKGROUND

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

PRINCIPAL FINDINGS

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

CONCLUSION

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

Increased copy number of the target gene squalene monooxygenase as the main resistance mechanism to terbinafine in Leishmania infantum

We use here two genomic screens in an attempt to understand the mode of action and resistance mechanism of terbinafine, an antifungal contemplated as a potential drug against the parasite Leishmania. One screen consisted in in vitro drug evolution where 5 independent mutants were selected step-by-step for terbinafine resistance. Sequencing of the genome of the 5 mutants revealed no single nucleotide polymorphisms related to the resistance phenotype. However, the ERG1 gene was found amplified as part of a linear amplicon, and transfection of ERG1 fully recapitulated the terbinafine resistance phenotype of the mutants. The second screen, Cos-seq, consisted in selecting a gene overexpression library with terbinafine followed by the sequencing of the enriched cosmids. This screen identified two cosmids derived from loci on chromosomes 13 and 29 encoding the squalene monooxygenase (ERG1) and the C8 sterol isomerase (ERG2), respectively. Transfection of the ERG1-cosmid, but not the ERG2-cosmid, produced resistance to terbinafine. Our screens suggest that ERG1 is the main, if not only, target for terbinafine in Leishmania and amplification of its gene is the main resistance mechanism.

Targeting trypanosomes: how chemogenomics and artificial intelligence can guide drug discovery

Trypanosomatids are protozoan parasites that cause human and animal neglected diseases. Despite global efforts, effective treatments are still much needed. Phenotypic screens have provided several chemical leads for drug discovery, but the mechanism of action for many of these chemicals is currently unknown. Recently, chemogenomic screens assessing the susceptibility or resistance of parasites carrying genome-wide modifications started to define the mechanism of action of drugs at large scale. In this review, we discuss how genomics is being used for drug discovery in trypanosomatids, how integration of chemical and genomics data from these and other organisms has guided prioritisations of candidate therapeutic targets and additional chemical starting points, and how these data can fuel the expansion of drug discovery pipelines into the era of artificial intelligence.

The challenges of the genome-based identification of antifungal resistance in the clinical routine

The increasing number of chronic and life-threatening infections caused by antimicrobial resistant fungal isolates is of critical concern. Low DNA sequencing cost may facilitate the identification of the genomic profile leading to resistance, the resistome, to rationally optimize the design of antifungal therapies. However, compared to bacteria, initiatives for resistome detection in eukaryotic pathogens are underdeveloped. Firstly, reported mutations in antifungal targets leading to reduced susceptibility must be extensively collected from the literature to generate comprehensive databases. This information should be complemented with specific laboratory screenings to detect the highest number possible of relevant genetic changes in primary targets and associations between resistance and other genomic markers. Strikingly, some drug resistant strains experience high-level genetic changes such as ploidy variation as much as duplications and reorganizations of specific chromosomes. Such variations involve allelic dominance, gene dosage increments and target expression regime effects that should be explicitly parameterized in antifungal resistome prediction algorithms. Clinical data indicate that predictors need to consider the precise pathogen species and drug levels of detail, instead of just genus and drug class. The concomitant needs for mutation accuracy and assembly quality assurance suggest hybrid sequencing approaches involving third-generation methods will be utilized. Moreover, fatal fast infections, like fungemia and meningitis, will further require both sequencing and analysis facilities are available in-house. Altogether, the complex nature of antifungal resistance demands extensive sequencing, data acquisition and processing, bioinformatic analysis pipelines, and standard protocols to be accomplished prior to genome-based protocols are applied in the clinical setting.

Thiophene derivatives activity against the protozoan parasite Leishmania infantum

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

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.

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

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

Mutation Characteristics and Phylogenetic Analysis of Five Leishmania Clinical Isolates

Simple Summary

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

Abstract

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

Laboratory Selection of Trypanosomatid Pathogens for Drug Resistance

The selection of parasites for drug resistance in the laboratory is an approach frequently used to investigate the mode of drug action, estimate the risk of emergence of drug resistance, or develop molecular markers for drug resistance. Here, we focused on the How rather than the Why of laboratory selection, discussing different experimental set-ups based on research examples with Trypanosoma brucei, Trypanosoma cruzi, and Leishmania spp. The trypanosomatids are particularly well-suited to illustrate different strategies of selecting for drug resistance, since it was with African trypanosomes that Paul Ehrlich performed such an experiment for the first time, more than a century ago. While breakthroughs in reverse genetics and genome editing have greatly facilitated the identification and validation of candidate resistance mutations in the trypanosomatids, the forward selection of drug-resistant mutants still relies on standard in vivo models and in vitro culture systems. Critical questions are: is selection for drug resistance performed in vivo or in vitro? With the mammalian or with the insect stages of the parasites? Under steady pressure or by sudden shock? Is a mutagen used? While there is no bona fide best approach, we think that a methodical consideration of these questions provides a helpful framework for selection of parasites for drug resistance in the laboratory.

Potential molecular targets of Leishmania pathways in developing novel antileishmanials

The illness known as leishmaniasis has not become a household name like malaria, although it stands as the second-largest parasitic disease, surpassed only by malaria. As no licensed vaccine is available, treatment for leishmaniasis mostly relies on chemotherapy. Inefficiency and drug resistance are the major impediments in current therapeutics. In this scenario, identification of novel molecular drug candidates is indispensable to develop robust antileishmanials. The exploration of structure-based drugs to target enzymes/molecules of Leishmania which differ structurally/functionally from their equivalents in mammalian hosts not only helps in developing a new class of antileishmanials, but also paves the way to understand Leishmania biology. This review provides a comprehensive overview on possible drug candidates relating to various Leishmania molecular pathways.

Insights from Leishmania (Viannia) guyanensis in vitro behavior and intercellular communication

BACKGROUND

Pentavalent antimonial-based chemotherapy is the first-line approach for leishmaniasis treatment and disease control. Nevertheless antimony-resistant parasites have been reported in some endemic regions. Treatment refractoriness is complex and is associated with patient- and parasite-related variables. Although amastigotes are the parasite stage in the vertebrate host and, thus, exposed to the drug, the stress caused by trivalent antimony in promastigotes has been shown to promote significant modification in expression of several genes involved in various biological processes, which will ultimately affect parasite behavior. Leishmania (Viannia) guyanensis is one of the main etiological agents in the Amazon Basin region, with a high relapse rate (approximately 25%).

METHODS

Herein, we conducted several in vitro analyses with L. (V.) guyanensis strains derived from cured and refractory patients after treatment with standardized antimonial therapeutic schemes, in addition to a drug-resistant in vitro-selected strain. Drug sensitivity assessed through Sb(III) half-maximal inhibitory concentration (IC₅₀) assays, growth patterns (with and without drug pressure) and metacyclic-like percentages were determined for all strains and compared to treatment outcomes. Finally, co-cultivation without intercellular contact was followed by parasitic density and Sb(III) IC₅₀ measurements.

RESULTS

Poor treatment response was correlated with increased Sb(III) IC₅₀ values. The decrease in drug sensitivity was associated with a reduced cell replication rate, increased in vitro growth ability, and higher metacyclic-like proportion. Additionally, in vitro co-cultivation assays demonstrated that intercellular communication enabled lower drug sensitivity and enhanced in vitro growth ability, regardless of direct cell contact.

CONCLUSIONS

Data concerning drug sensitivity in the Viannia subgenus are emerging, and L. (V.) guyanensis plays a pivotal epidemiological role in Latin America. Therefore, investigating the parasitic features potentially related to relapses is urgent. Altogether, the data presented here indicate that all tested strains of L. (V.) guyanensis displayed an association between treatment outcome and in vitro parameters, especially the drug sensitivity. Remarkably, sharing enhanced growth ability and decreased drug sensitivity, without intercellular communication, were demonstrated.

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

BACKGROUND

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

METHODS

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

RESULTS

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

CONCLUSIONS

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

Identification of Resistance Determinants for a Promising Antileishmanial Oxaborole Series

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

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.

Co-Therapy of Albendazole and Dexamethasone Reduces Pathological Changes in the Cerebral Parenchyma of Th-1 and Th-2 Dominant Mice Heavily Infected with Angiostrongylus cantonensis: Histopathological and RNA-seq Analyses

Administration of albendazole alone was not very suitable for the treatment of cerebral angiostrongyliasis. This study was designed to evaluate the effects of the co-therapy of this drug and dexamethasone in Th-1 and Th-2 dominant mice infected with Angiostrongylus cantonensis. Each of BALB/c and C57BL/6 mice infected with 50 A. cantonensis third-stage larvae were administered albendazole (10 mg/kg/day) alone, dexamethasone (0.5 mg/kg/day) alone, or co-therapy of the two drugs from day 7 or 14 post-infection for 7 or 14 days. After sacrifice, coronal slices were prepared from five brain regions and stained with hematoxylin and eosin. Eight pathological changes were employed to determine the therapeutic effectiveness using a scoring system. RNA-seq analysis was performed to confirm the histopathological findings. The infected BALB/c and C57BL/6 mice had similar patterns in the pathological changes. Meningitis, hemorrhage, size of worms, and encephalitis in the cerebral parenchyma were slighter in the mice treated with co-therapy than the remaining groups. Mice treated from day 14 had more severe changes than those from day 7. The histopathological findings were found to be consistent to immune responses determined by RNA-seq analysis. Co-therapy was determined to reduce pathological changes after administration to mice infected with A. cantonensis.

Extrachromosomal DNA amplicons in antimalarial-resistant Plasmodium falciparum

Extrachromosomal (ec) DNAs are genetic elements that exist separately from the genome. Since ecDNA can carry beneficial genes, they are a powerful adaptive mechanism in cancers and many pathogens. For the first time, we report ecDNA contributing to antimalarial resistance in Plasmodium falciparum, the most virulent human malaria parasite. Using pulse field gel electrophoresis combined with PCR-based copy number analysis, we detected two ecDNA elements that differ in migration and structure. Entrapment in the electrophoresis well and low susceptibility to exonucleases revealed that the biologically relevant ecDNA element is large and complex in structure. Using deep sequencing, we show that ecDNA originates from the chromosome and expansion of an ecDNA-specific sequence may improve its segregation or expression. We speculate that ecDNA is maintained using established mechanisms due to shared characteristics with the mitochondrial genome. Implications of ecDNA discovery in this organism are wide-reaching due to the potential for new strategies to target resistance development.

A Multi-Country, Single-Blinded, Phase 2 Study to Evaluate a Point-of-Need System for Rapid Detection of Leishmaniasis and Its Implementation in Endemic Settings

With the advancement of isothermal nucleic acid amplification techniques, detection of the pathogenic DNA in clinical samples at point-of-need is no longer a dream. The newly developed recombinase polymerase amplification (RPA) assay incorporated in a suitcase laboratory has shown promising diagnostic efficacy over real-time PCR in detection of leishmania DNA from clinical samples. For broader application of this point-of-need system, we undertook a current multi-country diagnostic evaluation study towards establishing this technique in different endemic settings which would be beneficial for the ongoing elimination programs for leishmaniasis. For this study purpose, clinical samples from confirmed visceral leishmaniasis (VL) and post-kala-azar dermal leishmaniasis (PKDL) patients were subjected to both real-time PCR and RPA assay in Bangladesh, India, and Nepal. Further skin samples from confirmed cutaneous leishmaniasis (CL) patients were also included from Sri Lanka. A total of 450 clinical samples from VL patients, 429 from PKDL patients, 47 from CL patients, and 322 from endemic healthy/healthy controls were under investigation to determine the diagnostic efficacy of RPA assay in comparison to real-time PCR. A comparative sensitivity of both methods was found where real-time PCR and RPA assay showed 96.86% (95% CI: 94.45–98.42) and 88.85% (95% CI: 85.08–91.96) sensitivity respectively in the diagnosis of VL cases. This new isothermal method also exhibited promising diagnostic sensitivity (93.50%) for PKDL cases, when a skin sample was used. Due to variation in the sequence of target amplicons, RPA assay showed comparatively lower sensitivity (55.32%) than that of real-time PCR in Sri Lanka for the diagnosis of CL cases. Except for India, the assay presented absolute specificity in the rest of the sites. Excellent concordance between the two molecular methods towards detection of leishmania DNA in clinical samples substantiates the application of RPA assay incorporated in a suitcase laboratory for point-of-need diagnosis of VL and PKDL in low resource endemic settings. However, further improvisation of the method is necessary for diagnosis of CL.

Genome Sequencing of Leishmania infantum Causing Cutaneous Leishmaniosis from a Turkish Isolate with Next-Generation Sequencing Technology

PURPOSE

Leishmania subgenus Leishmania causes leishmaniosis, which is a chronic systemic disease in humans and animals, in which the skin and visceral organs can be affected. The disease generally consists of three different clinical types in humans: visceral (kala-azar, VL), cutaneous (CL) and mucocutaneous leishmaniosis (MCL). According to the World Health Organization (WHO), leishmaniosis is still one of the world's most neglected diseases. It has been nearly 13-14 years since the completion of the first complete genome sequence of a Leishmania parasite. However, much information about these parasites remains to be elucidated, such as the causes of differences in tissue tropism. The aim of this study is to perform the whole-genome sequencing of Leishmania infantum causing cutaneous leishmaniosis from a Turkish isolate with next-generation sequencing technology.

METHODS

Genomic sequencing was performed on the Illumina HiSeq 2500 platform. The TruSeq Nano DNA Low Throughput Library Prep Kit, compatible with the Illumina HiSeq 2500 platform, was used to generate the library. Synthesis sequencing (SBS) was performed with a HiSeq Rapid SBS Kit v2 to generate single-fragment reads (2 × 150 bp; PE) with two fragment end-to-end assemblies. Bioinformatics analyses were performed on the Geneious 11.0.5. ( www.genius.com ) platform.

RESULTS

In our study, a high-quality whole-genome sequence (WGS) of L. infantum was successfully generated, and a total of 32,009,137 base pairs of genomic DNA from 36 chromosomes were obtained. The resulting genomic DNA sequence was submitted to the US National Center for Biotechnology Information (NCBI) GenBank ( www.ncbi.nlm.nih.gov ) database and registered under the name Leishmania infantum_TR01 (Lin_TR01). The following accession numbers were assigned by NCBI to the 36 chromosomes of the Lin_TR01 genome: CP027807, CP027810, CP027808, CP027811, CP027809, CP027812, CP027813, CP027814, CP027817, CP027818, CP027819, CP027815, CP027821, CP027816, CP027823, CP027820, CP027822, CP027824, CP027825, CP027826, CP027827, CP027828, CP027829, CP027830, CP027831, CP027832, CP027833, CP027834, CP027835, CP027836, CP027837, CP027838, CP027839, CP027840, CP027841, CP027842. As a result of the annotation of the Lin_TR01 genome, 3153 polymorphisms, 8324 genes, 8199 CDSs, 8109 mRNAs, 67 tRNAs, 11 rRNAs and 58 ncRNA were identified. Among the 8199 CDS obtained, 5278 encode hypothetical proteins.

CONCLUSION

In this study, a high-quality WGS of Leishmania infantum was successfully obtained for the first time in Turkey. According to a review of WGS studies on this subject, the Lin_TR01 strain is the first strain to be isolated from cutaneous leishmaniosis. The reference genome of L. infantum JPCM5 (Peacock et al., 2007) was obtained from a visceral leishmaniosis case, in accordance with the classical tissue and organ tropism of the species. Lin_TR01 is the second whole-genome-sequenced strain in the world after the JPCM5 strain. The Lin_TR01 genome is the only L. infantum whole-genome sequence that is completed assembly level from 36 chromosomes among the genomes obtained thus far ( https://www.ncbi.nlm.nih.gov/genome/genomes/249 ).

Extracellular Vesicles in Trypanosomatids: Host Cell Communication

Trypanosoma cruzi, Trypanosoma brucei and Leishmania (Trypanosomatidae: Kinetoplastida) are parasitic protozoan causing Chagas disease, African Trypanosomiasis and Leishmaniases worldwide. They are vector borne diseases transmitted by triatomine bugs, Tsetse fly, and sand flies, respectively. Those diseases cause enormous economic losses and morbidity affecting not only rural and poverty areas but are also spreading to urban areas. During the parasite-host interaction, those organisms release extracellular vesicles (EVs) that are crucial for the immunomodulatory events triggered by the parasites. EVs are involved in cell-cell communication and can act as important pro-inflammatory mediators. Therefore, interface between EVs and host immune responses are crucial for the immunopathological events that those diseases exhibit. Additionally, EVs from these organisms have a role in the invertebrate hosts digestive tracts prior to parasite transmission. This review summarizes the available data on how EVs from those medically important trypanosomatids affect their interaction with vertebrate and invertebrate hosts.

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.

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.

Antileishmanial activity evaluation of thiazolidine-2,4-dione against Leishmania infantum and Leishmania braziliensis

Leishmaniasis is responsible for approximately 65,000 annual deaths. Despite the mortality data, drugs available for the treatment of patients are insufficient and have moderate therapeutic efficacy in addition to serious adverse effects, which makes the development of new drugs urgent. To achieve this goal, the integration of kinetic and DSF assays against parasitic validated targets, along with phenotypic assays, can help the identification and optimization of bioactive compounds. Pteridine reductase 1 (PTR1), a validated target in Leishmania sp., is responsible for the reduction of folate and biopterin to tetrahydrofolate and tetrahydrobiopterin, respectively, both of which are essential for cell growth. In addition to the in vitro evaluation of 16 thiazolidine-2,4-dione derivatives against Leishmania major PTR1 (LmPTR1), using the differential scanning fluorimetry (ThermoFluor®), phenotypic assays were employed to evaluate the compound effect over Leishmania braziliensis (MHOM/BR/75/M2903) and Leishmania infantum (MHOM/BR/74/PP75) promastigotes viability. The ThermoFluor® results show that thiazolidine-2,4-dione derivatives have micromolar affinity to the target and equivalent activity on Leishmania cells. 2b is the most potent compound against L. infantum (EC₅₀ = 23.45 ± 4.54 μM), whereas 2a is the most potent against L. braziliensis (EC₅₀ = 44.16 ± 5.77 μM). This result suggests that lipophilic substituents on either-meta and/or-para positions of the benzylidene ring increase the potency against L. infantum. On the other hand, compound 2c (CE₅₀ = 49.22 ± 7.71 μM) presented the highest selectivity index.

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.

Antimicrobial activity of RP-1 peptide conjugate with ferrocene group

Parasitic diseases are a neglected and serious problem, especially in underdeveloped countries. Among the major parasitic diseases, Leishmaniasis figures as an urgent challenge due to its high incidence and severity. At the same time, the indiscriminate use of antibiotics by the population is increasing together with resistance to medicines. To address this problem, new antibiotic-like molecules that directly kill or inhibit the growth of microorganisms are necessary, where antimicrobial peptides (AMPs) can be of great help. In this work, the ferrocene molecule, one active compound with low levels of in vivo toxicity, was coupled to the N-terminus of the RP1 peptide (derived from the human chemokine CXCL4), aiming to evaluate how this change modifies the structure, biological activity, and toxicity of the peptide. The peptide and the conjugate were synthesized using the solid phase peptide synthesis (SPPS). Circular dichroism assays in PBS showed that the RP1 peptide and its conjugate had a typical spectrum for disordered structures. The Fc-RP1 presented anti-amastigote activity against Leishmania amazonensis (IC₅₀ = 0.25 μmol L^–1). In comparison with amphotericin B, a second-line drug approved for leishmaniasis treatment, (IC₅₀ = 0.63 μmol L^-1), Fc-RP1 was more active and showed a 2.5-fold higher selectivity index. The RP1 peptide presented a MIC of 4.3 μmol L^-1 against S. agalactiae, whilst Fc-RP1 was four times more active (MIC = 0.96 μmol L^-1), indicating that ferrocene improved the antimicrobial activity against Gram-positive bacteria. The Fc-RP1 peptide also decreased the minimum inhibitory concentration (MIC) in the assays against E. faecalis (MIC = 7.9 μmol L^-1), E. coli (MIC = 3.9 μmol L^-1) and S. aureus (MIC = 3.9 μmol L^-1). The cytotoxicity of the compounds was tested against HaCaT cells, and no significant activity at the highest concentration tested (500 μg. mL^-1) was observed, showing the high potential of this new compound as a possible new drug. The coupling of ferrocene also increased the vesicle permeabilization of the peptide, showing a direct relation between high peptide concentration and high carboxyfluorescein release, which indicates the action mechanism by pore formation on the vesicles. Several studies have shown that ferrocene destabilizes cell membranes through lipid peroxidation, leading to cell lysis. It is noteworthy that the Fc-RP1 peptide synthesized here is a prototype of a bioconjugation strategy, but it still is a compound with great biological activity against neglected and fish diseases.

Detecting sequence variants in clinically important protozoan parasites

Second and third generation sequencing methods are crucial for population genetic studies, and variant detection is a popular approach for exploiting this sequence data. While mini- and microsatellites are historically useful markers for studying important Protozoa such as Toxoplasma and Plasmodium spp., detecting non-repetitive variants such as those found in genes can be fundamental to investigating a pathogen's biology. These variants, namely single nucleotide polymorphisms and insertions and deletions, can help elucidate the genetic basis of an organism's pathogenicity, identify selective pressures, and resolve phylogenetic relationships. They also have the added benefit of possessing a comparatively low mutation rate, which contributes to their stability. However, there is a plethora of variant analysis tools with nuanced pipelines and conflicting recommendations for best practise, which can be confounding. This lack of standardisation means that variant analysis requires careful parameter optimisation, an understanding of its limitations, and the availability of high quality data. This review explores the value of variant detection when applied to non-model organisms such as clinically important protozoan pathogens. The limitations of current methods are discussed, including special considerations that require the end-users' attention to ensure that the results generated are reproducible, and the biological conclusions drawn are valid.

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.

Genomewide Analysis of Mode of Action of the S-Adenosylmethionine Analogue Sinefungin in Leishmania infantum

To further our understanding of one-carbon metabolism in the protozoan parasite Leishmania, we conducted genomic screens to study how the parasite responded to sinefungin (SNF) selection. SNF is a structural analogue of S-adenosylmethionine (AdoMet), a key methyl group donor to a number of biomolecules. One screen consisted of sequencing SNF-resistant mutants generated by stepwise selection with gradually increasing drug concentrations. These studies demonstrated deletion of the AdoMet transporter (AdoMetT1) by intergenic recombination as a crucial loss-of-function marker for SNF resistance. The second screen consisted of Cos-seq, a gain-of-function cosmid-based genomewide functional screen with increasing SNF concentration coupled to next-generation sequencing. Cosmids enriched in that screen and sequenced led to the identification of (i) the AdoMet synthetase (METK) as the major SNF target, (ii) an mRNA [(guanine-N7)-methyltransferase (CMT1)], (iii) a leucine carboxyl methyltransferase (LCMT), (iv) two tryparedoxin genes, and (v) two protein phosphatase regulatory genes. Further functional exploration indicated that LCMT interacts with one phosphatase catalytic subunit (PP2AC) and that mutation of the C-terminal leucine residue of PP2AC affects sinefungin susceptibility. These holistic screens led to the identification of transporters, biosynthetic genes, RNA and protein methyltransferases, as well as phosphatases linked to AdoMet-mediated functions in Leishmania IMPORTANCE The two main cellular metabolic one-carbon donors are reduced folates and S-adenosylmethionine, whose biosynthetic pathways have proven highly effective in chemotherapeutic interventions in various cell types. Sinefungin, a nucleoside analogue of S-adenosylmethionine, was shown to have potent activity against the protozoan parasite Leishmania Here, we studied resistance to sinefungin using whole-genome approaches as a way to further our understanding of the role of S-adenosylmethionine in this parasite and to reveal novel potential drug targets. These approaches allowed the characterization of novel features related to S-adenosylmethionine function in Leishmania which could further help in the development of sinefungin-like compounds against this pathogenic parasite.

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.

A single amino acid substitution (H451Y) in Leishmania calcium-dependent kinase SCAMK confers high tolerance and resistance to antimony

Background

For almost a century, antimonials have remained the first-line drugs for the treatment of leishmaniasis. However, little is known about their mode of action and clinical resistance mechanisms.

Objectives

We have previously shown that Leishmania nicotinamidase (PNC1) is an essential enzyme for parasite NAD+ homeostasis and virulence in vivo. Here, we found that parasites lacking the pnc1 gene (Δpnc1) are hypersusceptible to the active form of antimony (SbIII) and used these mutant parasites to better understand antimony’s mode of action and the mechanisms leading to resistance.

Methods

SbIII-resistant WT and Δpnc1 parasites were selected in vitro by a stepwise selection method. NAD(H)/NADP(H) dosages and quantitative RT–PCR experiments were performed to explain the susceptibility differences observed between strains. WGS and a marker-free CRISPR/Cas9 base-editing approach were used to identify and validate the role of a new resistance mutation.

Results

NAD+-depleted Δpnc1 parasites were highly susceptible to SbIII and this phenotype could be rescued by NAD+ precursor or trypanothione precursor supplementation. Δpnc1 parasites could become resistant to SbIII by an unknown mechanism. WGS revealed a unique amino acid substitution (H451Y) in an EF-hand domain of an orphan calcium-dependent kinase, recently named SCAMK. When introduced into a WT reference strain by base editing, the H451Y mutation allowed Leishmania parasites to survive at extreme concentrations of SbIII, potentiating the rapid emergence of resistant parasites.

Conclusions

These results establish that Leishmania SCAMK is a new central hub of antimony’s mode of action and resistance development, and uncover the importance of drug tolerance mutations in the evolution of parasite drug resistance.

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

Background

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

Methods

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

Results

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

Conclusions

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

Electronic supplementary material

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

Molecular Analysis of Aquaglyceroporin 1 Gene in Non-Healing Clinical Isolates Obtained from Patients with Cutaneous Leishmaniasis from Central of Iran

BACKGROUND

Regarding the antimonial-resistant of Leishmania spp., understanding of related mechanism is necessary. One of the most important involved molecules is aquaglyceropin1 (AQP1). The aim of this study was molecular analysis of AQP1 gene from antimonial-resistant clinical isolates and its expression.

METHODS

Overall, 150 patients with cutaneous leishmaniasis referring to the reference laboratories of Yazd and Varzaneh,, located 105km southeast of Isfahan and 240km away from Yazd, were assessed from Jun 2015 to Dec 2017. After sampling, staining was done and evaluated for Leishman by microscope. Samples were collected in RNAlater solution for gene expression analysis in non-healing isolates. DNA extraction was performed from each slide with Leishman body. All patients with L. major isolates detected by ITS1-PCR-RFLP were followed for finding the resistant isolates, consequence of molecular characterization of AQP1 using PCR-RFLP. Gene expression of AQP1 from all resistant isolates was assessed in comparison with the one in a sensitive isolate. Statistical analysis was done using SPSS. The significance level was considered ≤0.05.

RESULTS

Five isolates were detected as antimonial resistant. Molecular detection and identification were appeared that all were L. major. The molecular characterization of AQP1 showed G562A mutation. Gene expression of AQP1 in resistant isolates showed 1.67 fold higher than the sensitive isolate.

CONCLUSION

We reported a new point mutation of G562A in AQP1 gene involved in molecular mechanism in resistant isolates.

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.

Cos-Seq: A High-Throughput Gain-of-Function Screen for Drug Resistance Studies in Leishmania

Leishmania is still a major cause of mortality and morbidity worldwide. Few efficient drugs are available, and resistance threatens actual treatments. In order to improve knowledge about the mode of action of current drugs and those in development, as well as to understand the mechanisms pertaining to their resistance, we recently described a sensitive and high-throughput method termed Cos-Seq. Here we provide a detailed protocol for every step of the procedure, from library construction to drug selection, cosmid extraction, and next-generation sequencing of extracted cosmids. A section on the bioinformatics of Cos-Seq is also included. Cos-Seq facilitates the identification of gain-of-function resistance mechanisms and drug targets and is a useful tool in resistance and drug development studies.

A Guide to Next Generation Sequence Analysis of Leishmania Genomes

Next generation sequencing (NGS) technology transformed Leishmania genome studies and became an indispensable tool for Leishmania researchers. Recent Leishmania genomics analyses facilitated the discovery of various genetic diversities including single nucleotide polymorphisms (SNPs), copy number variations (CNVs), somy variations, and structural variations in detail and provided valuable insights into the complexity of the genome and gene regulation. Many aspects of Leishmania NGS analyses are similar to those of related pathogens like trypanosomes. However, the analyses of Leishmania genomes face a unique challenge because of the presence of frequent aneuploidy. This makes characterization and interpretation of read depth and somy a key part of Leishmania NGS analyses because read depth affects the accuracy of detection of all genetic variations. However, there are no general guidelines on how to explore and interpret the impact of aneuploidy, and this has made it difficult for biologists and bioinformaticians, especially for beginners, to perform their own analyses and interpret results across different analyses. In this guide we discuss a wide range of topics essential for Leishmania NGS analyses, ranging from how to set up a computational environment for genome analyses, to how to characterize genetic variations among Leishmania samples, and we will particularly focus on chromosomal copy number variation and its impact on genome analyses.

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.

Resistance of Leishmania infantum to allopurinol is associated with chromosome and gene copy number variations including decrease in the S-adenosylmethionine synthetase (METK) gene copy number

Leishmania infantum is one of the causative agents of visceral leishmaniasis (VL), a widespread, life-threatening disease. This parasite is responsible for the majority of human VL cases in Brazil, the Middle East, China, Central Asia and the Mediterranean basin. Its main reservoir are domestic dogs which, similar to human patients, may develop severe visceral disease and die if not treated. The drug allopurinol is used for the long-term maintenance of dogs with canine leishmaniasis. Following our report of allopurinol resistance in treated relapsed dogs, we investigated the mechanisms and markers of resistance to this drug. Whole genome sequencing (WGS) of clinical resistant and susceptible strains, and laboratory induced resistant parasites, was carried out in order to detect genetic changes associated with resistance. Significant gene copy number variation (CNV) was found between resistant and susceptible isolates at several loci, including a locus on chromosome 30 containing the genes LinJ.30.3550 through LinJ.30.3580. A reduction in copy number for LinJ.30.3560, encoding the S-adenosylmethionine synthetase (METK) gene, was found in two resistant clinical isolates and four induced resistant clonal strains. Using quantitative real time PCR, this reduction in METK copy number was also found in three additional resistant clinical isolates. Furthermore, inhibition of S-adenosylmethionine synthetase encoded by the METK gene in allopurinol susceptible strains resulted in increased allopurinol resistance, confirming its role in resistance to allopurinol. In conclusion, this study identified genetic changes associated with L. infantum resistance to allopurinol and the reduction in METK copy number identified may serve as a marker for resistance in dogs, and reduced protein activity correlated with increased allopurinol resistance.

•

Allopurinol resistance was previously described in L. infantum isolated from dogs.

•

This study aimed at defining the genetic differences between susceptible and resistant strains.

•

Gene and chromosome copy numbers differed between susceptible and resistant L. infantum strains.

•

Decrease in METK gene copies was associated with increased allopurinol resistance.

•

Inhibition of the enzyme encoded by METK increased allopurinol resistance.

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

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

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

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

Genomic and transcriptomic alterations in Leishmania donovani lines experimentally resistant to antileishmanial drugs

Leishmaniasis is a serious medical issue in many countries around the World, but it remains largely neglected in terms of research investment for developing new control and treatment measures. No vaccines exist for human use, and the chemotherapeutic agents currently used are scanty. Furthermore, for some drugs, resistance and treatment failure are increasing to alarming levels. The aim of this work was to identify genomic and trancriptomic alterations associated with experimental resistance against the common drugs used against VL: trivalent antimony (Sb^III, S line), amphotericin B (AmB, A line), miltefosine (MIL, M line) and paromomycin (PMM, P line). A total of 1006 differentially expressed transcripts were identified in the S line, 379 in the A line, 146 in the M line, and 129 in the P line. Also, changes in ploidy of chromosomes and amplification/deletion of particular regions were observed in the resistant lines regarding the parental one. A series of genes were identified as possible drivers of the resistance phenotype and were validated in both promastigotes and amastigotes from Leishmania donovani, Leishmania infantum and Leishmania major species. Remarkably, a deletion of the gene LinJ.36.2510 (coding for 24-sterol methyltransferase, SMT) was found to be associated with AmB-resistance in the A line. In the P line, a dramatic overexpression of the transcripts LinJ.27.T1940 and LinJ.27.T1950 that results from a massive amplification of the collinear genes was suggested as one of the mechanisms of PMM resistance. This conclusion was reinforced after transfection experiments in which significant PMM-resistance was generated in WT parasites over-expressing either gene LinJ.27.1940 (coding for a D-lactate dehydrogenase-like protein, D-LDH) or gene LinJ.27.1950 (coding for an aminotransferase of branched-chain amino acids, BCAT). This work allowed to identify new drivers, like SMT, the deletion of which being associated with resistance to AmB, and the tandem D-LDH-BCAT, the amplification of which being related to PMM resistance.

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

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

•

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

•

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

•

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

•

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

•

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

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

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

Molecular targets and pathways for the treatment of visceral leishmaniasis

Visceral leishmaniasis (VL) represents the most severe form of the tropical disease, leishmaniasis. Treatment of VL is complicated because of the few clinically approved antileishmanial drugs available; emerging resistance to first-line drugs; need for a temperature-controlled 'cold' supply chain; serious toxicity concerns over drugs such as amphotericin B; high cost of medication; and unavailability of clinically approved antileishmanial vaccines. Attacking potential molecular targets, specific to the parasite, is a vital step in the treatment of this and other infectious diseases. As we discuss here, comprehensive investigation of these targets could provide a promising strategy for the treatment of visceral leishmaniasis.

Induction of allopurinol resistance in Leishmania infantum isolated from dogs

Resistance to allopurinol in zoonotic canine leishmaniasis has been recently shown to be associated with disease relapse in naturally-infected dogs. However, information regarding the formation of resistance and its dynamics is lacking. This study describes the successful in-vitro induction of allopurinol resistance in Leishmania infantum cultured under increasing drug pressure. Allopurinol susceptibility and growth rate of induced parasites were monitored over 23 weeks and parasite clones were tested at selected time points and compared to their parental lines, both as promastigotes and as amastigotes. Allopurinol resistance was formed in strains from two parasite stocks producing a 20-fold rise in IC₅₀ along three distinct growth phases. In addition, characteristic differential clustering of single nucleotide polymorphisms (SNP) was found in drug sensitive and resistant parasite clones. Results confirm that genetic polymorphism, as well as clonal heterogeneity, contribute to in-vitro resistance to allopurinol, which is likely to occur in natural infection.

Visceral leishmaniasis caused by the parasite Leishmania infantum is a neglected tropical disease transmitted from animal hosts to humans by sand fly bites. This potentially fatal disease affects thousands of people annually and threatens millions who live in disease risk areas. Domestic dogs are considered as the main reservoir of this parasite which can also cause a severe chronic canine disease. Allopurinol is the main drug used for long term treatment of this disease but it often does not eliminate infection in dogs. We have recently demonstrated that allopurinol resistant parasites can be isolated from naturally infected dogs that have developed clinical recurrence of disease during allopurinol treatment. In this study we aimed to see if resistance can be induced in susceptible parasite strains isolated from sick dogs by growing them in increasing drug concentrations under laboratory conditions. The changes in allopurinol susceptibility were measured and the impact of drug on parasite growth was monitored over 23 weeks. Induction of resistance was successful producing parasites 20-folds less susceptible to the drug. The pattern of change in drug susceptibility suggests that a genetic change is responsible for the increased resistance which is likely to mimic the formation of resistance in dogs.