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

Leishmania genomic adaptation: more than just a 36-body problem

Genome instability has been identified as a major driver of adaptation in fast-growing, eukaryotic cells, including fungi, protists, or cancer. How these cells cope with the toxic effects caused by such copy number variations remains to be elucidated. In recent years, the protist parasites Leishmania spp. have emerged as interesting model pathogens to assess this open question and to study the role of its intrinsic genome instability in fitness gain in culture, experimental infection, and in the field. Here we summarize recent results on Leishmania genomic adaptation and propose thought-provoking evolutionary concepts new to the Leishmania field that need to be considered when mapping genotype-to-phenotype relationships in molecular and epidemiological studies.

Nanopore sequencing reveals that DNA replication compartmentalisation dictates genome stability and instability in Trypanosoma brucei

The Trypanosoma brucei genome is structurally complex. Eleven megabase-sized chromosomes each comprise a transcribed core flanked by silent subtelomeres, housing thousands of Variant Surface Glycoprotein (VSG) genes. Additionally, hundreds of sub-megabase chromosomes contain 177 bp repeats of unknown function, and VSG transcription sites localise to many telomeres. DNA replication dynamics have only been described in the megabase chromosome cores, and in the single active VSG transcription site. Using a Nanopore genome assembly, we show that megabase chromosome subtelomeres display a paucity of replication initiation events relative to the core, correlating with increased instability. In addition, replication of the active VSG transcription site is shown to originate from the telomere, likely causing targeted VSG recombination. Lastly, we provide evidence that the 177 bp repeats act as conserved DNA replication origins, explaining submegabase chromosome stability. Compartmentalized DNA replication therefore explains how T. brucei balances stable genome transmission with localised instability driving immune evasion.

Genetic coping mechanisms observed in Leishmania tropica, from the Middle East region, enhance the survival of the parasite after drug exposure

INTRODUCTION

Cutaneous leishmaniasis caused by L. tropica is common in the Middle East and treatment failure and drug resistance are known to occur. Several genetic mechanisms: aneuploidy, recombination and loss of heterozygosity, single nucleotide polymorphism (SNP) changes, copy number variation (CNV), and mutation of the H locus associated with drug resistance have been described.

MATERIALS AND METHODS

We studied SNP and CNV patterns in 22 isolates of L. tropica from Afghanistan, Iran and Syria in a geographic, phylogenetic and antimony exposure context.

RESULTS

A high SNP frequency was observed in isolates from Syria on chromosome 23, including the H locus, linked to different ancestry at that chromosome segment. Among the isolates from Afghanistan and Iran, an elevated frequency of nonsynonymous SNPs was observed on several chromosomes. Changes in CNV patterns were seen in isolates exposed to drug pressure, especially for the ferric iron reductase gene. Expanded genes were categorised into five functional categories: translational elongation, mitochondrial transmembrane transport, positive regulation of cellular component organisation, response to stimulus and response to hypoxia. No CNV was identified at the H locus, the MAPK1 gene, the APQ1 gene, nor chromosomes 23, 31 or 36 regardless of previous antimonial exposure.

DISCUSSION

In our study, Leishmania tropica had a jump in the nonsynonymous SNP rates at chromosome 23, including the H locus. CNV was observed among isolates exposed to antimonials, especially involving the gene encoding a ferric iron reductase. Several essential genetic coping mechanisms in the cell were enhanced when exposed to antimony, possibly for the survival of the parasite. Our work supports the perspective that Leishmania uses several mechanisms to adapt to environmental changes and drug exposure.

Molecular Mechanisms of Drug Resistance in Leishmania spp

The protozoan parasite Leishmania causes leishmaniasis, a neglected tropical disease, that disproportionately affects underdeveloped countries. This disease has major health, economic, and social implications, particularly because of the limited treatment options, high cost, the severe side effects associated with available therapeutics, and the high rate of treatment failure caused by the parasites' growing resistance to current medications. In this review, we describe first the common strategies used by pathogens to develop drug resistance and then focus on the arsenal of available drugs to treat leishmaniasis, their modes of action, and the molecular mechanisms contributing to drug resistance in Leishmania spp., including the role of genomic, transcriptional, and translational control. We focus more specifically on our recent discovery of translational reprogramming as a major driver of drug resistance leading to coordinated changes in the translation of transcripts and orchestrating changes in metabolome and lipidome to support drug resistance. A thorough understanding of these mechanisms is essential to identify the key elements needed to combat resistance and improve leishmaniasis treatment methods.

Antiprotozoal peptide prediction using machine learning with effective feature selection techniques

Background

Protozoal pathogens pose a considerable threat, leading to notable mortality rates and the ongoing challenge of developing resistance to drugs. This situation underscores the urgent need for alternative therapeutic approaches. Antimicrobial peptides stand out as promising candidates for drug development. However, there is a lack of published research focusing on predicting antimicrobial peptides specifically targeting protozoal pathogens. In this study, we introduce a successful machine learning-based framework designed to predict potential antiprotozoal peptides effective against protozoal pathogens.

Objective

The primary objective of this study is to classify and predict antiprotozoal peptides using diverse negative datasets.

Methods

A comprehensive literature review was conducted to gather experimentally validated antiprotozoal peptides, forming the positive dataset for our study. To construct a robust machine learning classifier, multiple negative datasets were incorporated, including (i) non-antimicrobial, (ii) antiviral, (iii) antibacterial, (iv) antifungal, and (v) antimicrobial peptides excluding those targeting protozoal pathogens. Various compositional features of the peptides were extracted using the pfeature algorithm. Two feature selection methods, SVC-L1 and mRMR, were employed to identify highly relevant features crucial for distinguishing between the positive and negative datasets. Additionally, five popular classifiers i.e. Decision Tree, Random Forest, Support Vector Machine, Logistic Regression, and XGBoost were used to build efficient decision models.

Results

XGBoost was the most effective in classifying antiprotozoal peptides from each negative dataset based on the features selected by the mRMR feature selection method. The proposed machine learning framework efficiently differentiate the antiprotozoal peptides from (i) non-antimicrobial (ii) antiviral (iii) antibacterial (iv) antifungal and (v) antimicrobial with accuracy of 97.27 %, 93.64 %, 86.36 %, 90.91 %, and 89.09 % respectively on the validation dataset.

Conclusion

The models are incorporated in a user-friendly web server (www.soodlab.com/appred) to predict the antiprotozoal activity of given peptides.

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.

Circular extrachromosomal DNA in Euglena gracilis under normal and stress conditions

Extrachromosomal circular DNA (eccDNA) enhances genomic plasticity, augmenting its coding and regulatory potential. Advances in high-throughput sequencing have enabled the investigation of these structural variants. Although eccDNAs have been investigated in numerous taxa, they remained understudied in euglenids. Therefore, we examined eccDNAs predicted from Illumina sequencing data of Euglena gracilis Z SAG 1224-5/25, grown under optimal photoperiod and exposed to UV irradiation. We identified approximately 1000 unique eccDNA candidates, about 20% of which were shared across conditions. We also observed a significant enrichment of mitochondrially encoded eccDNA in the UV-irradiated sample. Furthermore, we found that the heterogeneity of eccDNA was reduced in UV-exposed samples compared to cells that were grown in optimal conditions. Hence, eccDNA appears to play a role in the response to oxidative stress in Euglena, as it does in other studied organisms. In addition to contributing to the understanding of Euglena genomes, our results contribute to the validation of bioinformatics pipelines on a large, non-model genome.

ABC transporter inhibition by beauvericin partially overcomes drug resistance in Leishmania tropica

Leishmaniasis is a neglected tropical disease infecting the world's poorest populations. Miltefosine (ML) remains the primary oral drug against the cutaneous form of leishmaniasis. The ATP-binding cassette (ABC) transporters are key players in the xenobiotic efflux, and their inhibition could enhance the therapeutic index. In this study, the ability of beauvericin (BEA) to overcome ABC transporter-mediated resistance of Leishmania tropica to ML was assessed. In addition, the transcription profile of genes involved in resistance acquisition to ML was inspected. Finally, we explored the efflux mechanism of the drug and inhibitor. The efficacy of ML against all developmental stages of L. tropica in the presence or absence of BEA was evaluated using an absolute quantification assay. The expression of resistance genes was evaluated, comparing susceptible and resistant strains. Finally, the mechanisms governing the interaction between the ABC transporter and its ligands were elucidated using molecular docking and dynamic simulation. Relative quantification showed that the expression of the ABCG sub-family is mostly modulated by ML. In this study, we used BEA to impede resistance of Leishmania tropica. The IC50 values, following BEA treatment, were significantly reduced from 30.83, 48.17, and 16.83 µM using ML to 8.14, 11.1, and 7.18 µM when using a combinatorial treatment (ML + BEA) against promastigotes, axenic amastigotes, and intracellular amastigotes, respectively. We also demonstrated a favorable BEA-binding enthalpy to L. tropica ABC transporter compared to ML. Our study revealed that BEA partially reverses the resistance development of L. tropica to ML by blocking the alternate ATP hydrolysis cycle.

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.

LeishIF4E2 is a cap-binding protein that plays a role in Leishmania cell cycle progression

Leishmania encode six paralogs of the cap-binding protein eIF4E and five eIF4G candidates, forming unique complexes. Two cap-binding proteins, LeishIF4E1 and LeishIF4E2, do not bind any identified LeishIF4Gs, thus their roles are intriguing. Here, we combine structural prediction, proteomic analysis, and interaction assays to shed light on LeishIF4E2 function. A nonconserved C-terminal extension was identified through structure prediction and sequence alignment. m7 GTP-binding assays involving both recombinant and transgenic LeishIF4E2 with and without the C-terminal extension revealed that this extension functions as a regulatory gate, modulating the cap-binding activity of LeishIF4E2. The interactomes of the two LeishIF4E2 versions were investigated, highlighting the role of the C-terminal extension in binding to SLBP2. SLBP2 is known to interact with a stem-loop structure in the 3' UTRs of histone mRNAs. Consistent with the predicted inhibitory effect of SLBP2 on histone expression in Xenopus laevis, a hemizygous deletion mutant of LeishIF4E2, exhibited an upregulation of several histones. We therefore propose that LeishIF4E2 is involved in histone expression, possibly through its interaction between SLBP2 and LeishIF4E2, thus affecting cell cycle progression. In addition, cell synchronization showed that LeishIF4E2 expression decreased during the S-phase, when histones are known to be synthesized. Previous studies in T. brucei also highlighted an association between TbEIF4E2 and SLBP2, and further reported on an interaction between TbIF4E2 and S-phase-abundant mRNAs. Our results show that overexpression of LeishIF4E2 correlates with upregulation of cell cycle and chromosome maintenance proteins. Along with its effect on histone expression, we propose that LeishIF4E2 is involved in cell cycle progression.

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.

Host, parasite and drug determinants of clinical outcomes following treatment of visceral leishmaniasis: a protocol for individual participant data meta-analysis

INTRODUCTION

Visceral leishmaniasis (VL) is a parasitic disease with an estimated 30 000 new cases occurring annually. There is an observed variation in the efficacy of the current first-line therapies across different regions. Such heterogeneity could be a function of host, parasite and drug factors. An individual participant data meta-analysis (IPD-MA) is planned to explore the determinants of treatment outcomes.

METHODS AND ANALYSIS

The Infectious Diseases Data Observatory (IDDO) VL living systematic review (IDDO VL LSR) library is an open-access resource of all published therapeutic studies in VL since 1980. For this current review, the search includes all clinical trials published between 1 January 1980 and 2 May 2021. Studies indexed in the IDDO VL LSR library were screened for eligibility for inclusion in this IPD-MA. Corresponding authors and principal investigators of the studies meeting the eligibility criteria for inclusion were invited to be part of the collaborative IPD-MA. Authors agreeing to participate in this collaborative research were requested to share the IPD using the IDDO VL data platform. The IDDO VL data platform currently holds data sets from clinical trials standardised to a common data format and provides a unique opportunity to identify host, parasite and drug determinants of treatment outcomes. Multivariable regression models will be constructed to identify determinants of therapeutic outcomes using generalised linear mixed-effects models accounting for within-study site clustering.

ETHICS AND DISSEMINATION

This IPD-MA meets the criteria for waiver of ethical review as defined by the Oxford Tropical Research Ethics Committee (OxTREC) granted to IDDO, as the research consists of secondary analysis of existing anonymised data (Exempt granted on 29 March 2023, OxTREC REF: IDDO) Ethics approval was granted by the ICMR-Rajendra Memorial Research Institute of Medical Sciences ethics committee (Letter no: RMRI/EC/30/2022) on 04-07-2022. The results of this IPD-MA will be disseminated at conferences, IDDO website and any peer-reviewed publications. All publications will be open source. Findings of this research will be critically important for the control programmes at regional/global levels, policy makers and groups developing new VL treatments.

PROSPERO REGISTRATION

CRD42021284622.

Shining the spotlight on the neglected: new high-quality genome assemblies as a gateway to understanding the evolution of Trypanosomatidae

BACKGROUND

Protists of the family Trypanosomatidae (phylum Euglenozoa) have gained notoriety as parasites affecting humans, domestic animals, and agricultural plants. However, the true extent of the group's diversity spreads far beyond the medically and veterinary relevant species. We address several knowledge gaps in trypanosomatid research by undertaking sequencing, assembly, and analysis of genomes from previously overlooked representatives of this protistan group.

RESULTS

We assembled genomes for twenty-one trypanosomatid species, with a primary focus on insect parasites and Trypanosoma spp. parasitizing non-human hosts. The assemblies exhibit sizes consistent with previously sequenced trypanosomatid genomes, ranging from approximately 18 Mb for Obscuromonas modryi to 35 Mb for Crithidia brevicula and Zelonia costaricensis. Despite being the smallest, the genome of O. modryi has the highest content of repetitive elements, contributing nearly half of its total size. Conversely, the highest proportion of unique DNA is found in the genomes of Wallacemonas spp., with repeats accounting for less than 8% of the assembly length. The majority of examined species exhibit varying degrees of aneuploidy, with trisomy being the most frequently observed condition after disomy.

CONCLUSIONS

The genome of Obscuromonas modryi represents a very unusual, if not unique, example of evolution driven by two antidromous forces: i) increasing dependence on the host leading to genomic shrinkage and ii) expansion of repeats causing genome enlargement. The observed variation in somy within and between trypanosomatid genera suggests that these flagellates are largely predisposed to aneuploidy and, apparently, exploit it to gain a fitness advantage. High heterogeneity in the genome size, repeat content, and variation in chromosome copy numbers in the newly-sequenced species highlight the remarkable genome plasticity exhibited by trypanosomatid flagellates. These new genome assemblies are a robust foundation for future research on the genetic basis of life cycle changes and adaptation to different hosts in the family Trypanosomatidae.

Condition-dependent fitness effects of large synthetic chromosome amplifications

Whole-chromosome aneuploidy and large segmental amplifications can have devastating effects in multicellular organisms, from developmental disorders and miscarriage to cancer. Aneuploidy in single-celled organisms such as yeast also results in proliferative defects and reduced viability. Yet, paradoxically, CNVs are routinely observed in laboratory evolution experiments with microbes grown in stressful conditions. The defects associated with aneuploidy are often attributed to the imbalance of many differentially expressed genes on the affected chromosomes, with many genes each contributing incremental effects. An alternate hypothesis is that a small number of individual genes are large effect 'drivers' of these fitness changes when present in an altered copy number. To test these two views, we have employed a collection of strains bearing large chromosomal amplifications that we previously assayed in nutrient-limited chemostat competitions. In this study, we focus on conditions known to be poorly tolerated by aneuploid yeast-high temperature, treatment with the Hsp90 inhibitor radicicol, and growth in extended stationary phase. To identify potential genes with a large impact on fitness, we fit a piecewise constant model to fitness data across chromosome arms, filtering breakpoints in this model by magnitude to focus on regions with a large impact on fitness in each condition. While fitness generally decreased as the length of the amplification increased, we were able to identify 91 candidate regions that disproportionately impacted fitness when amplified. Consistent with our previous work with this strain collection, nearly all candidate regions were condition specific, with only five regions impacting fitness in multiple conditions.

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

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

Ribosome Specialization in Protozoa Parasites

Ribosomes, in general, are viewed as constitutive macromolecular machines where protein synthesis takes place; however, this view has been recently challenged, supporting the hypothesis of ribosome specialization and opening a completely new field of research. Recent studies have demonstrated that ribosomes are heterogenous in their nature and can provide another layer of gene expression control by regulating translation. Heterogeneities in ribosomal RNA and ribosomal proteins that compose them favor the selective translation of different sub-pools of mRNAs and functional specialization. In recent years, the heterogeneity and specialization of ribosomes have been widely reported in different eukaryotic study models; however, few reports on this topic have been made on protozoa and even less on protozoa parasites of medical importance. This review analyzes heterogeneities of ribosomes in protozoa parasites highlighting the specialization in their functions and their importance in parasitism, in the transition between stages in their life cycle, in the change of host and in response to environmental conditions.

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.

Comparative Draft Genomes of Leishmania orientalis Isolate PCM2 (Formerly Named Leishmania siamensis) and Leishmania martiniquensis Isolate PCM3 from the Southern Province of Thailand

Simple Summary

This study successfully sequenced the draft genomes of the southern isolates of Leishmania orientalis and Leishmania martiniquensis in Thailand. The comparison with the genomes of the northern isolates revealed species-level similarity with a level of genome and proteome variation, suggesting the different emerging strains. Comparing the proteins of these southern strains to those of the northern ones and 14 other Leishmania species showed six protein groups with numerous unique proteins: 53 for the southern strain PCM2 of L. orientalis and 97 for the strain PCM3 of L. martiniquensis. Some of these proteins were related to virulence, drug resistance, drug target, and stress response, which could be targeted for further experimental characterization. Therefore, the findings could initiate further genetic and population genomic investigation, and the close monitoring of L. orientalis and L. martiniquensis in Thailand and neighboring regions.

Abstract

(1) Background: Autochthonous leishmaniasis, a sandfly-borne disease caused by the protozoan parasites Leishmania orientalis (formerly named Leishmania siamensis) and Leishmania martiniquensis, has been reported for immunocompromised and immunocompetent patients in the southern province of Thailand. Apart from the recent genomes of the northern isolates, limited information is known on the emergence and genetics of these parasites. (2) Methods: This study sequenced and compared the genomes of L. orientalis isolate PCM2 and L. martiniquensis isolate PCM3 with those of the northern isolates and other 14 Leishmania species using short-read whole-genome sequencing methods and comparative bioinformatic analyses. (3) Results: The genomes of the southern isolates of L. orientalis and L. martiniquensis were 30.01 Mbp and 32.39 Mbp, and the comparison with the genomes of the northern isolates revealed species-level similarity with a level of genome and proteome variation, suggesting the different strains. Comparative proteome analysis showed six protein groups with 53 unique proteins for the strain PCM2 and 97 for the strain PCM3. Certain proteins were related to virulence, drug resistance, and stress response. (4) Conclusion: Therefore, the findings could indicate the need for more genetic and population genomic investigation, and the close monitoring of L. orientalis and L. martiniquensis in Thailand and neighboring regions.

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.

Targeted Deletion of Centrin in Leishmania braziliensis Using CRISPR-Cas9-Based Editing

Leishmania braziliensis is the main causative agent of Tegumentary Leishmaniasis in the Americas. However, difficulties related to genome manipulation, experimental infection, and parasite growth have so far limited studies with this species. CRISPR-Cas9-based technology has made genome editing more accessible, and here we have successfully employed the LeishGEdit approach to attenuate L. braziliensis. We generated a transgenic cell line expressing Cas9 and T7 RNA polymerase, which was employed for the targeted deletion of centrin, a calcium-binding cytoskeletal protein involved in the centrosome duplication in eukaryotes. Centrin-deficient Leishmania exhibit growth arrest at the amastigote stage. Whole-genome sequencing of centrin-deficient L. braziliensis (LbCen^−/−) did not indicate the presence of off-target mutations. In vitro, the growth rates of LbCen^−/− and wild-type promastigotes were similar, but axenic and intracellular LbCen^−/− amastigotes showed a multinucleated phenotype with impaired survival following macrophage infection. Upon inoculation into BALB/c mice, LbCen^−/− were detected at an early time point but failed to induce lesion formation, contrary to control animals, infected with wild-type L. braziliensis. A significantly lower parasite burden was also observed in mice inoculated with LbCen^−/−, differently from control mice. Given that centrin-deficient Leishmania sp. have become candidates for vaccine development, we propose that LbCen^−/− can be further explored for the purposes of immunoprophylaxis against American Tegumentary Leishmaniasis.

Sex under pressure: stress facilitates Leishmania in vitro hybridization

The selection of Leishmania hybrids in axenic culture was considered rare until recently, when Louradour and Ferreira et al., demonstrated that induced DNA damage facilitates genetic exchange, resulting in full genome tetraploid progenies in vitro. Meiosis-related gene homologues HAP2, GEX1, and RAD51 were found to be involved, opening new avenues for functional genomic studies.

Specific Human ATR and ATM Inhibitors Modulate Single Strand DNA Formation in Leishmania major Exposed to Oxidative Agent

Leishmania parasites are the causative agents of a group of neglected tropical diseases known as leishmaniasis. The molecular mechanisms employed by these parasites to adapt to the adverse conditions found in their hosts are not yet completely understood. DNA repair pathways can be used by Leishmania to enable survival in the interior of macrophages, where the parasite is constantly exposed to oxygen reactive species. In higher eukaryotes, DNA repair pathways are coordinated by the central protein kinases ataxia telangiectasia mutated (ATM) and ataxia telangiectasia and Rad3 related (ATR). The enzyme Exonuclease-1 (EXO1) plays important roles in DNA replication, repair, and recombination, and it can be regulated by ATM- and ATR-mediated signaling pathways. In this study, the DNA damage response pathways in promastigote forms of L. major were investigated using bioinformatics tools, exposure of lineages to oxidizing agents and radiation damage, treatment of cells with ATM and ATR inhibitors, and flow cytometry analysis. We demonstrated high structural and important residue conservation for the catalytic activity of the putative LmjEXO1. The overexpression of putative LmjEXO1 made L. major cells more susceptible to genotoxic damage, most likely due to the nuclease activity of this enzyme and the occurrence of hyper-resection of DNA strands. These cells could be rescued by the addition of caffeine or a selective ATM inhibitor. In contrast, ATR-specific inhibition made the control cells more susceptible to oxidative damage in an LmjEXO1 overexpression-like manner. We demonstrated that ATR-specific inhibition results in the formation of extended single-stranded DNA, most likely due to EXO1 nucleasic activity. Antagonistically, ATM inhibition prevented single-strand DNA formation, which could explain the survival phenotype of lineages overexpressing LmjEXO1. These results suggest that an ATM homolog in Leishmania could act to promote end resection by putative LmjEXO1, and an ATR homologue could prevent hyper-resection, ensuring adequate repair of the parasite DNA.

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.

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.

In vitro selection of ketoconazole-pentamidine-resistant Leishmania (Viannia) braziliensis strains

The use of ketoconazole (KTZ) plus pentamidine (PMD) could be an interesting treatment option for New World cutaneous leishmaniasis. The aim of this work was to generate KTZ- and PMD-resistant strains and to determine some characteristics of the selection process and the resulting parasites. Resistance to one or two drugs was selected on promastigotes by progressively increasing drug concentrations for eleven months. The resistance levels (IC₅₀) to one or two drugs (synergism assay) were determined using a colorimetric resazurin methodology. The stability of the resistance phenotype (without drug pressure or after mouse passage), cross resistance with paromomycin and miltefosine, and resistance transference to intracellular amastigotes were determined. In addition, some parasite attributes compared with WT, such as growth kinetics, amastigogenesis, THP-1 cells, and mouse infection, were determined. Promastigotes resistant to KTZ or PMD were obtained three times earlier than the combined KTZ + PMD-resistant strains. Resistant parasites (promastigotes and intracellular amastigotes) were three to twelve times less susceptible to KTZ and PMD than WT parasites. The resistance phenotype on parasites was unstable, and no cross resistance was observed. Similar parasite fitness related to our evaluated characteristics was observed except for in vivo infection, where a delay of the onset of cutaneous lesions was observed after KTZ + PMD-resistant parasite infection. CONCLUSION: Combined treatment with KTZ and PMD delayed the onset of parasite resistance and was more effective in vitro than each drug separately for WT and all resistant strains. Parasites resistant to KTZ and PMD acquired similar in vitro behaviour to WT parasites, were less virulent to mice and maintained their resistance phenotype on intracellular amastigotes but not without drug pressure or after mouse infection.

Decreased glutamate transport in acivicin resistant Leishmania tarentolae

Studies of drug resistance in the protozoan parasites of the genus Leishmania have been helpful in revealing biochemical pathways as potential drug targets. The chlorinated glutamine analogue acivicin has shown good activity against Leishmania cells and was shown to target several enzymes containing amidotransferase domains. We selected a Leishmania tarentolae clone for acivicin resistance. The genome of this resistant strain was sequenced and the gene coding for the amidotransferase domain-containing GMP synthase was found to be amplified. Episomal expression of this gene in wild-type L. tarentolae revealed a modest role in acivicin resistance. The most prominent defect observed in the resistant mutant was reduced uptake of glutamate, and through competition experiments we determined that glutamate and acivicin, but not glutamine, share the same transporter. Several amino acid transporters (AATs) were either deleted or mutated in the resistant cells. Some contributed to the acivicin resistance phenotype although none corresponded to the main glutamate transporter. Through sequence analysis one AAT on chromosome 22 corresponded to the main glutamate transporter. Episomal expression of the gene coding for this transporter in the resistant mutant restored glutamate transport and acivicin susceptibility. Its genetic knockout led to reduced glutamate transport and acivicin resistance. We propose that acivicin binds covalently to this transporter and as such leads to decreased transport of glutamate and acivicin thus leading to acivicin resistance.

Studies of drug resistance in the protozoan parasites of the genus Leishmania have been helpful in revealing biochemical pathways as potential drug targets. Here we report on the characterization at the genomics and metabolomics levels of a L. tarentolae strain made resistant to acivicin, an analogue of glutamine with activity against this parasite. We found that resistance to acivicin is accompanied by a reduced uptake and intracellular levels of glutamate and that both are expected to share the same transporter. Through gene overexpression and disruption studies we identified the main amino acid transporter responsible for glutamate uptake.

Ku80 is involved in telomere maintenance but dispensable for genomic stability in Leishmania mexicana

BACKGROUND

Telomeres are indispensable for genome stability maintenance. They are maintained by the telomere-associated protein complex, which include Ku proteins and a telomerase among others. Here, we investigated a role of Ku80 in Leishmania mexicana. Leishmania is a genus of parasitic protists of the family Trypanosomatidae causing a vector-born disease called leishmaniasis.

METHODOLOGY/PRINCIPAL FINDINGS

We used the previously established CRISPR/Cas9 system to mediate ablation of Ku80- and Ku70-encoding genes in L. mexicana. Complete knock-outs of both genes were confirmed by Southern blotting, whole-genome Illumina sequencing, and RT-qPCR. Resulting telomeric phenotypes were subsequently investigated using Southern blotting detection of terminal restriction fragments. The genome integrity in the Ku80- deficient cells was further investigated by whole-genome sequencing. Our work revealed that telomeres in the ΔKu80 L. mexicana are elongated compared to those of the wild type. This is a surprising finding considering that in another model trypanosomatid, Trypanosoma brucei, they are shortened upon ablation of the same gene. A telomere elongation phenotype has been documented in other species and associated with a presence of telomerase-independent alternative telomere lengthening pathway. Our results also showed that Ku80 appears to be not involved in genome stability maintenance in L. mexicana.

CONCLUSION/SIGNIFICANCE

Ablation of the Ku proteins in L. mexicana triggers telomere elongation, but does not have an adverse impact on genome integrity.

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.

Going ballistic: Leishmania nuclear subversion of host cell plasticity

Intracellular parasites have evolved intricate strategies to subvert host cell functions for their own survival. These strategies are particularly damaging to the host if the infection involves immune cells, as illustrated by protozoan parasites of the genus Leishmania that thrive inside mononuclear phagocytic cells, causing devastating immunopathologies. While the impact of Leishmania infection on host cell phenotype and functions has been well documented, the regulatory mechanisms underlying host cell subversion were only recently investigated. Here we summarize the current knowledge on how Leishmania infection affects host nuclear activities and propose thought-provoking new concepts on the reciprocal relationship between epigenetic and transcriptional regulation in host cell phenotypic plasticity, its potential subversion by the intracellular parasite, and its relevance for host-directed therapy.

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

Background

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

Methods

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

Results

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

Conclusions

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

Graphical abstract

Supplementary Information

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

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

BACKGROUND

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

METHODS

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

RESULTS

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

CONCLUSION

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

Well-Tolerated Amphotericin B Derivatives That Effectively Treat Visceral Leishmaniasis

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

Life Cycle Stage-Specific Accessibility of Leishmania donovani Chromatin at Transcription Start Regions

Leishmania donovani is a parasitic protist that causes the lethal Kala-azar fever in India and East Africa. Gene expression in Leishmania is regulated by gene copy number variation and inducible translation while RNA synthesis initiates at a small number of sites per chromosome and proceeds through polycistronic transcription units, precluding a gene-specific regulation (C. Clayton and M. Shapira, Mol Biochem Parasitol 156:93–101, 2007, https://doi.org/10.1016/j.molbiopara.2007.07.007). Here, we analyze the dynamics of chromatin structure in both life cycle stages of the parasite and find evidence for an additional, epigenetic gene regulation pathway in this early branching eukaryote. The assay for transposase-accessible chromatin using sequencing (ATAC-seq) analysis (J. D. Buenrostro, P. G. Giresi, L. C. Zaba, H. Y. Chang, and W. J. Greenleaf, Nat Methods 10:1213–1218, 2013, https://doi.org/10.1038/nmeth.2688) predominantly shows euchromatin at transcription start regions in fast-growing promastigotes, but mostly heterochromatin in the slowly proliferating amastigotes, the mammalian stage, reflecting a previously shown increase of histone synthesis in the latter stage.

IMPORTANCELeishmania parasites are important pathogens with a global impact and cause poverty-related illness and death. They are devoid of classic cis- and trans-acting transcription regulators but use regulated translation and gene copy number variations to adapt to hosts and environments. In this work, we show that transcription start regions present as open euchromatin in fast-growing insect stages but as less-accessible heterochromatin in the slowly proliferating amastigote stage, indicating an epigenetic control of gene accessibility in this early branching eukaryotic pathogen. This finding should stimulate renewed interest in the control of RNA synthesis in Leishmania and related parasites.

Preclinical Studies in Anti-Trypanosomatidae Drug Development

The trypanosomatid parasites Trypanosoma brucei, Trypanosoma cruzi and Leishmania are the causative agents of human African trypanosomiasis, Chagas Disease and Leishmaniasis, respectively. These infections primarily affect poor, rural communities in the developing world, and are responsible for trapping sufferers and their families in a disease/poverty cycle. The development of new chemotherapies is a priority given that existing drug treatments are problematic. In our search for novel anti-trypanosomatid agents, we assess the growth-inhibitory properties of >450 compounds from in-house and/or "Pathogen Box" (PBox) libraries against L. infantum, L. amazonensis, L.braziliensis, T. cruzi and T. brucei and evaluate the toxicities of the most promising agents towards murine macrophages. Screens using the in-house series identified 17 structures with activity against and selective toward Leishmania: Compounds displayed 50% inhibitory concentrations between 0.09 and 25 μM and had selectivity index values >10. For the PBox library, ~20% of chemicals exhibited anti-parasitic properties including five structures whose activity against L. infantum had not been reported before. These five compounds displayed no toxicity towards murine macrophages over the range tested with three being active in an in vivo murine model of the cutaneous disease, with 100% survival of infected animals. Additionally, the oral combination of three of them in the in vivo Chagas disease murine model demonstrated full control of the parasitemia. Interestingly, phenotyping revealed that the reference strain responds differently to the five PBox-derived chemicals relative to parasites isolated from a dog. Together, our data identified one drug candidate that displays activity against Leishmania and other Trypanosomatidae in vitro and in vivo, while exhibiting low toxicity to cultured mammalian cells and low in vivo acute toxicity.

Metabolite Biomarkers of Leishmania Antimony Resistance

Leishmania parasites cause leishmaniasis, one of the most epidemiologically important neglected tropical diseases. Leishmania exhibits a high ability of developing drug resistance, and drug resistance is one of the main threats to public health, as it is associated with increased incidence, mortality, and healthcare costs. The antimonial drug is the main historically implemented drug for leishmaniasis. Nevertheless, even though antimony resistance has been widely documented, the mechanisms involved are not completely understood. In this study, we aimed to identify potential metabolite biomarkers of antimony resistance that could improve leishmaniasis treatment. Here, using L. tropica promastigotes as the biological model, we showed that the level of response to antimony can be potentially predicted using ¹H-NMR-based metabolomic profiling. Antimony-resistant parasites exhibited differences in metabolite composition at the intracellular and extracellular levels, suggesting that a metabolic remodeling is required to combat the drug. Simple and time-saving exometabolomic analysis can be efficiently used for the differentiation of sensitive and resistant parasites. Our findings suggest that changes in metabolite composition are associated with an optimized response to the osmotic/oxidative stress and a rearrangement of carbon-energy metabolism. The activation of energy metabolism can be linked to the high energy requirement during the antioxidant stress response. We also found that metabolites such as proline and lactate change linearly with the level of resistance to antimony, showing a close relationship with the parasite's efficiency of drug resistance. A list of potential metabolite biomarkers is described and discussed.

Combined gene deletion of dihydrofolate reductase-thymidylate synthase and pteridine reductase in Leishmania infantum

Our understanding of folate metabolism in Leishmania has greatly benefited from studies of resistance to the inhibitor methotrexate (MTX). Folates are reduced in Leishmania by the bifunctional dihydrofolate reductase thymidylate synthase (DHFR-TS) and by pteridine reductase (PTR1). To further our understanding of folate metabolism in Leishmania, a Cos-seq genome-wide gain of function screen was performed against MTX and against the two thymidylate synthase (TS) inhibitors 5-fluorouracil and pemetrexed. The screen revealed DHFR-TS and PTR1 but also the nucleoside transporter NT1 and one hypothetical gene derived from chromosome 31. For MTX, the concentration of folate in the culture medium affected the enrichment pattern for genes retrieved by Cos-seq. We generated a L. infantum DHFR-TS null mutant that was thymidine auxotroph, a phenotype that could be rescued by the addition of thymidine or by transfection of the flavin dependent bacterial TS gene ThyX. In these DHFR-TS null mutants it was impossible to obtain a chromosomal null mutant of PTR1 except if DHFR-TS or PTR1 were provided episomally. The transfection of ThyX however did not allow the elimination of PTR1 in a DHFR-TS null mutant. Leishmania can survive without copies of either DHFR-TS or PTR1 but not without both. Provided that our results observed with the insect stage parasites are also replicated with intracellular parasites, it would suggest that antifolate therapy in Leishmania would only work if both DHFR-TS and PTR1 would be targeted simultaneously.

The protozoan parasite Leishmania is auxotroph for folate and unconjugated pterins and salvages both from the mammalian host. Two enzymes of the folate metabolism pathway, namely the bifunctional dihydrofolate reductase-thymidylate synthase (DHFR-TS) and the pteridine reductase 1 (PTR1), are being evaluated for drug discovery and repurposing of existing anti-metabolites. Despite their apparent potential, development of DHFR-TS and PTR1 targeted chemotherapy against Leishmania is still awaiting. Here we revisited folate metabolism at the genomic level and report on the identification of known resistance genes alongside some new ones. Through gene disruption studies we found that L. infantum DHFR-TS null mutants are thymidine auxotroph and that these can be rescued by the bacterial flavin dependent thymidylate synthase ThyX. We also found that PTR1 is essential in the absence of a functional DHFR-TS even in the presence of ThyX or thymidine supplementation, indicating the essential role of reduced pterins or folate beyond thymidine synthesis. This study indicates that simultaneous targeting of DHFR-TS and PTR1 will be required for the development of anti-folate chemotherapy against Leishmania.

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.

Drug-Induced Lipid Remodeling in Leishmania Parasites

Leishmania parasites efficiently develop resistance against several types of drugs including antimonials, the primary antileishmanial drug historically implemented. The resistance to antimonials is considered to be a major risk factor for effective leishmaniasis treatment. To detect biomarkers/biopatterns for the differentiation of antimony-resistant Leishmania strains, we employed untargeted global mass spectrometry to identify intracellular lipids present in antimony sensitive and resistant parasites before and after antimony exposure. The lipidomic profiles effectively differentiated the sensitive and resistant phenotypes growing with and without antimony pressure. Resistant phenotypes were characterized by significant downregulation of phosphatidylcholines, sphingolipid decrease, and lysophosphatidylcholine increase, while sensitive phenotypes were characterized by the upregulation of triglycerides with long-chain fatty acids and a tendency toward the phosphatidylethanolamine decrease. Our findings suggest that the changes in lipid composition in antimony-resistant parasites contribute to the physiological response conducted to combat the oxidative stress unbalance caused by the drug. We have identified several lipids as potential biomarkers associated with the drug resistance.

Analysis of Base Excision and Single-Strand Break Repair Activities in Trypanosomatid Extracts

Cellular DNA is inherently unstable, subject to both spontaneous hydrolysis and attack by a range of exogenous and endogenous chemicals as well as physical agents such as ionizing and ultraviolet radiation. For parasitic protists, where an inoculum of infectious parasites is typically small and natural infections are often chronic with low parasitemia, they are also vulnerable to DNA damaging agents arising from innate immune defenses. The majority of DNA damage consists of relatively minor changes to the primary structure of the DNA, such as base deamination, oxidation, or alkylation and scission of the phosphodiester backbone. Yet these small changes can have serious consequences, often being mutagenic or cytotoxic. Cells have therefore evolved efficient mechanisms to repair such damage, with base excision and single strand break repair playing the primary role here. In this chapter we describe a method for analyzing the activity from cell extracts of various enzymes involved in the base excision and single strand break repair pathways of trypanosomatid parasites.

Colonization and genetic diversification processes of Leishmania infantum in the Americas

Leishmania infantum causes visceral leishmaniasis, a deadly vector-borne disease introduced to the Americas during the colonial era. This non-native trypanosomatid parasite has since established widespread transmission cycles using alternative vectors, and human infection has become a significant concern to public health, especially in Brazil. A multi-kilobase deletion was recently detected in Brazilian L. infantum genomes and is suggested to reduce susceptibility to the anti-leishmanial drug miltefosine. We show that deletion-carrying strains occur in at least 15 Brazilian states and describe diversity patterns suggesting that these derive from common ancestral mutants rather than from recurrent independent mutation events. We also show that the deleted locus and associated enzymatic activity is restored by hybridization with non-deletion type strains. Genetic exchange appears common in areas of secondary contact but also among closely related parasites. We examine demographic and ecological scenarios underlying this complex L. infantum population structure and discuss implications for disease control.

Philipp Schwabl, Mariana Boité, and colleagues analyze 126 Leishmania infantum genomes to determine how demographic and selective consequences of the parasite’s invasive history have contributed to intricate population genetic heterogeneity across Brazil. Their data suggest a complex interplay of population expansion, secondary contact and genetic exchange events underlying diversity patterns at short and long-distance scales. These processes also appear pivotal to the proliferation of a drug resistance-associated multi-gene deletion on chromosome 31.

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.

Phosphoglycerate kinase: structural aspects and functions, with special emphasis on the enzyme from Kinetoplastea

Phosphoglycerate kinase (PGK) is a glycolytic enzyme that is well conserved among the three domains of life. PGK is usually a monomeric enzyme of about 45 kDa that catalyses one of the two ATP-producing reactions in the glycolytic pathway, through the conversion of 1,3-bisphosphoglycerate (1,3BPGA) to 3-phosphoglycerate (3PGA). It also participates in gluconeogenesis, catalysing the opposite reaction to produce 1,3BPGA and ADP. Like most other glycolytic enzymes, PGK has also been catalogued as a moonlighting protein, due to its involvement in different functions not associated with energy metabolism, which include pathogenesis, interaction with nucleic acids, tumorigenesis progression, cell death and viral replication. In this review, we have highlighted the overall aspects of this enzyme, such as its structure, reaction kinetics, activity regulation and possible moonlighting functions in different protistan organisms, especially both free-living and parasitic Kinetoplastea. Our analysis of the genomes of different kinetoplastids revealed the presence of open-reading frames (ORFs) for multiple PGK isoforms in several species. Some of these ORFs code for unusually large PGKs. The products appear to contain additional structural domains fused to the PGK domain. A striking aspect is that some of these PGK isoforms are predicted to be catalytically inactive enzymes or ‘dead’ enzymes. The roles of PGKs in kinetoplastid parasites are analysed, and the apparent significance of the PGK gene duplication that gave rise to the different isoforms and their expression in Trypanosoma cruzi is discussed.

Genomic Diversification, Structural Plasticity, and Hybridization in Leishmania (Viannia) braziliensis

Leishmania (Viannia) braziliensis is an important Leishmania species circulating in several Central and South American countries. Among Leishmania species circulating in Brazil, Argentina and Colombia, L. braziliensis has the highest genomic variability. However, genomic variability at the whole genome level has been only studied in Brazilian and Peruvian isolates; to date, no Colombian isolates have been studied. Considering that in Colombia, L. braziliensis is a species with great clinical and therapeutic relevance, as well as the role of genetic variability in the epidemiology of leishmaniasis, we analyzed and evaluated intraspecific genomic variability of L. braziliensis from Colombian and Bolivian isolates and compared them with Brazilian isolates. Twenty-one genomes were analyzed, six from Colombian patients, one from a Bolivian patient, and 14 Brazilian isolates downloaded from public databases. The results obtained of Phylogenomic analysis showed the existence of four well-supported clades, which evidenced intraspecific variability. The whole-genome analysis revealed structural variations in the somy, mainly in the Brazilian genomes (clade 1 and clade 3), low copy number variations, and a moderate number of single-nucleotide polymorphisms (SNPs) in all genomes analyzed. Interestingly, the genomes belonging to clades 2 and 3 from Colombia and Brazil, respectively, were characterized by low heterozygosity (~90% of SNP loci were homozygous) and regions suggestive of loss of heterozygosity (LOH). Additionally, we observed the drastic whole genome loss of heterozygosity and possible hybridization events in one genome belonging to clade 4. Unique/shared SNPs between and within the four clades were identified, revealing the importance of some of them in biological processes of L. braziliensis. Our analyses demonstrate high genomic variability of L. braziliensis in different regions of South America, mainly in Colombia and suggest that this species exhibits striking genomic diversity and a capacity of genomic hybridization; additionally, this is the first study to report whole-genome sequences of Colombian L. braziliensis isolates.

The Leishmania donovani SENP Protease Is Required for SUMO Processing but Not for Viability

The protozoan parasite Leishmania donovani is part of an early eukaryotic branch and depends on post-transcriptional mechanisms for gene expression regulation. This includes post-transcriptional protein modifications, such as protein phosphorylation. The presence of genes for protein SUMOylation, i.e., the covalent attachment of small ubiquitin-like modifier (SUMO) polypeptides, in the Leishmania genomes prompted us to investigate the importance of the sentrin-specific protease (SENP) and its putative client, SUMO, for the vitality and infectivity of Leishmania donovani. While SENP null mutants are viable with reduced vitality, viable SUMO null mutant lines could not be obtained. SUMO C-terminal processing is disrupted in SENP null mutants, preventing SUMO from covalent attachment to proteins and nuclear translocation. Infectivity in vitro is not affected by the loss of SENP-dependent SUMO processing. We conclude that SENP is required for SUMO processing, but that functions of unprocessed SUMO are critical for Leishmania viability.

Sensing Host Arginine Is Essential for Leishmania Parasites' Intracellular Development

In this study, we report that the ability of the human pathogen Leishmania to sense and monitor the lack of arginine in the phagolysosome of the host macrophage is essential for disease development. Phagolysosomes of macrophages are the niche where Leishmania resides and causes human leishmaniasis. During infection, the arginine concentration in the phagolysosome decreases as part of the host innate immune response. An arginine sensor on the Leishmania cell surface activates an arginine deprivation response pathway that upregulates the expression of a parasite arginine transporter (AAP3). Here, we use CRISPR/Cas9-mediated disruption of the AAP3 locus to show that this response enables Leishmania parasites to successfully compete with the host macrophage in the “hunger games” for arginine.

Arginine homeostasis in lysosomes is critical for the growth and metabolism of mammalian cells. Phagolysosomes of macrophages are the niche where the parasitic protozoan Leishmania resides and causes human leishmaniasis. During infection, parasites encounter arginine deprivation, which is monitored by a sensor on the parasite cell surface. The sensor promptly activates a mitogen-activated protein kinase 2 (MAPK2)-mediated arginine deprivation response (ADR) pathway, resulting in upregulating the abundance and activity of the Leishmania arginine transporter (AAP3). Significantly, the ADR is also activated during macrophage infection, implying that arginine levels within the host phagolysosome are limiting for growth. We hypothesize that ADR-mediated upregulation of AAP3 activity is necessary to withstand arginine starvation, suggesting that the ADR is essential for parasite intracellular development. CRISPR/Cas9-mediated disruption of the AAP3 locus yielded mutants that retain a basal level of arginine transport but lack the ability to respond to arginine starvation. While these mutants grow normally in culture, they were impaired in their ability to develop inside THP-1 macrophages and were ∼70 to 80% less infective in BALB/c mice. Hence, inside the host macrophage, Leishmania must overcome the arginine “hunger games” by upregulating the transport of arginine via the ADR. We show that the ability to monitor and respond to changes in host metabolite levels is essential for pathogenesis.

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.