Search

We sequenced Leishmania donovani genomes in blood samples collected in emerging foci of visceral leishmaniasis in western Nepal. We detected lineages very different from the preelimination main parasite population, including a new lineage and a rare one previously reported in eastern Nepal. Our findings underscore the need for genomic surveillance.

Viruses are the most abundant biological entities on Earth and play a significant role in the evolution of many organisms and ecosystems. In pathogenic protozoa, the presence of viruses has been linked to an increased risk of treatment failure and severe clinical outcome. Here, we studied the molecular epidemiology of the zoonotic disease cutaneous leishmaniasis in Peru and Bolivia through a joint evolutionary analysis of Leishmania braziliensis and their dsRNA Leishmania virus 1. We show that parasite populations circulate in tropical rainforests and are associated with single viral lineages that appear in low prevalence. In contrast, groups of hybrid parasites are geographically and ecologically more dispersed and associated with an increased prevalence, diversity and spread of viruses. Our results suggest that parasite gene flow and hybridization increased the frequency of parasite-virus symbioses, a process that may change the epidemiology of leishmaniasis in the region., (© 2023. The Author(s).)

Introduction: Resistance against anti- Leishmania drugs (DR) has been studied for years, giving important insights into long-term adaptations of these parasites to drugs, through genetic modifications. However, microorganisms can also survive lethal drug exposure by entering into temporary quiescence, a phenomenon called drug tolerance (DT), which is rather unexplored in Leishmania ., Methods: We studied a panel of nine Leishmania braziliensis strains highly susceptible to potassium antimonyl tartrate (PAT), exposed promastigotes to lethal PAT pressure, and compared several cellular and molecular parameters distinguishing DT from DR., Results and Discussion: We demonstrated in vitro that a variable proportion of cells remained viable, showing all the criteria of DT and not of DR: i) signatures of quiescence, under drug pressure: reduced proliferation and significant decrease of rDNA transcription; ii) reversibility of the phenotype: return to low IC 50 after removal of drug pressure; and iii) absence of significant genetic differences between exposed and unexposed lineages of each strain and absence of reported markers of DR. We found different levels of quiescence and DT among the different L. braziliensis strains. We provide here a new in-vitro model of drug-induced quiescence and DT in Leishmania . Research should be extended in vivo , but the current model could be further exploited to support R&D, for instance, to guide the screening of compounds to overcome the quiescence resilience of the parasite, thereby improving the therapy of leishmaniasis., Competing Interests: The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest., (Copyright © 2023 Jara, Arevalo, Llanos-Cuentas, Broeck, Domagalska and Dujardin.)

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

We discovered a hybrid Leishmania parasite in Costa Rica that is genetically similar to hybrids from Panama. Genome analyses demonstrated the hybrid is triploid and identified L. braziliensis and L. guyanensis-related strains as parents. Our findings highlight the existence of poorly sampled Leishmania (Viannia) variants infectious to humans.

Leishmania aethiopica is a zoonotic Old World parasite transmitted by Phlebotomine sand flies and causing cutaneous leishmaniasis in Ethiopia and Kenya. Despite a range of clinical manifestations and a high prevalence of treatment failure, L. aethiopica is one of the most neglected species of the Leishmania genus in terms of scientific attention. Here, we explored the genome diversity of L. aethiopica by analyzing the genomes of twenty isolates from Ethiopia. Phylogenomic analyses identified two strains as interspecific hybrids involving L. aethiopica as one parent and L. donovani and L. tropica respectively as the other parent. High levels of genome-wide heterozygosity suggest that these two hybrids are equivalent to F1 progeny that propagated mitotically since the initial hybridization event. Analyses of allelic read depths further revealed that the L. aethiopica - L. tropica hybrid was diploid and the L. aethiopica - L. donovani hybrid was triploid, as has been described for other interspecific Leishmania hybrids. When focusing on L. aethiopica , we show that this species is genetically highly diverse and consists of both asexually evolving strains and groups of recombining parasites. A remarkable observation is that some L. aethiopica strains showed an extensive loss of heterozygosity across large regions of the nuclear genome, which likely arose from gene conversion/mitotic recombination. Hence, our prospection of L. aethiopica genomics revealed new insights into the genomic consequences of both meiotic and mitotic recombination in Leishmania ., Competing Interests: The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest., (Copyright © 2023 Hadermann, Heeren, Maes, Dujardin, Domagalska and Van den Broeck.)

Treatment failure (TF) jeopardizes the management of parasitic diseases, including leishmaniasis. From the parasite's point of view, drug resistance (DR) is generally considered as central to TF. However, the link between TF and DR, as measured by in vitro drug susceptibility assays, is unclear, some studies revealing an association between treatment outcome and drug susceptibility, others not. Here we address three fundamental questions aiming to shed light on these ambiguities. First, are the right assays being used to measure DR? Second, are the parasites studied, which are generally those that adapt to in vitro culture, actually appropriate? Finally, are other parasite factors - such as the development of quiescent forms that are recalcitrant to drugs - responsible for TF without DR?, Competing Interests: Declaration of interests The authors declare no competing interests., (Copyright © 2023 Elsevier Ltd. All rights reserved.)

Aneuploidy causes system-wide disruptions in the stochiometric balances of transcripts, proteins, and metabolites, often resulting in detrimental effects for the organism. The protozoan parasite Leishmania has an unusually high tolerance for aneuploidy, but the molecular and functional consequences for the pathogen remain poorly understood. Here, we addressed this question in vitro and present the first integrated analysis of the genome, transcriptome, proteome, and metabolome of highly aneuploid Leishmania donovani strains. Our analyses unambiguously establish that aneuploidy in Leishmania proportionally impacts the average transcript- and protein abundance levels of affected chromosomes, ultimately correlating with the degree of metabolic differences between closely related aneuploid strains. This proportionality was present in both proliferative and non-proliferative in vitro promastigotes. However, as in other Eukaryotes, we observed attenuation of dosage effects for protein complex subunits and in addition, non-cytoplasmic proteins. Differentially expressed transcripts and proteins between aneuploid Leishmania strains also originated from non-aneuploid chromosomes. At protein level, these were enriched for proteins involved in protein metabolism, such as chaperones and chaperonins, peptidases, and heat-shock proteins. In conclusion, our results further support the view that aneuploidy in Leishmania can be adaptive. Additionally, we believe that the high karyotype diversity in vitro and absence of classical transcriptional regulation make Leishmania an attractive model to study processes of protein homeostasis in the context of aneuploidy and beyond., Competing Interests: The authors have declared that no competing interests exist.

The implementation of prospective drug resistance (DR) studies in the research-and-development (R&D) pipeline is a common practice for many infectious diseases but not for neglected tropical diseases (NTDs). Here, we explored and demonstrated the importance of this approach using as paradigms Leishmania donovani, the etiological agent of visceral leishmaniasis (VL), and TCMDC-143345, a promising compound of the GlaxoSmithKline (GSK) "Leishbox" to treat VL. We experimentally selected resistance to TCMDC-143345 in vitro and characterized resistant parasites at the genomic and phenotypic levels. We found that it took more time to develop resistance to TCMDC-143345 than to other drugs in clinical use and that there was no cross-resistance to these drugs, suggesting a new and unique mechanism. By whole-genome sequencing, we found two mutations in the gene encoding the L. donovani dynamin-1-like protein (LdoDLP1) that were fixed at the highest drug pressure. Through phylogenetic analysis, we identified LdoDLP1 as a family member of the dynamin-related proteins, a group of proteins that impacts the shapes of biological membranes by mediating fusion and fission events, with a putative role in mitochondrial fission. We found that L. donovani lines genetically engineered to harbor the two identified LdoDLP1 mutations were resistant to TCMDC-143345 and displayed altered mitochondrial properties. By homology modeling, we showed how the two LdoDLP1 mutations may influence protein structure and function. Taken together, our data reveal a clear involvement of LdoDLP1 in the adaptation/reduced susceptibility of L. donovani to TCMDC-143345. IMPORTANCE Humans and their pathogens are continuously locked in a molecular arms race during which the eventual emergence of pathogen drug resistance (DR) seems inevitable. For neglected tropical diseases (NTDs), DR is generally studied retrospectively once it has already been established in clinical settings. We previously recommended to keep one step ahead in the host-pathogen arms race and implement prospective DR studies in the R&D pipeline, a common practice for many infectious diseases but not for NTDs. Here, using Leishmania donovani, the etiological agent of visceral leishmaniasis (VL), and TCMDC-143345, a promising compound of the GSK Leishbox to treat VL, as paradigms, we experimentally selected resistance to the compound and proceeded to genomic and phenotypic characterization of DR parasites. The results gathered in the present study suggest a new DR mechanism involving the L. donovani dynamin-1-like protein (LdoDLP1) and demonstrate the practical relevance of prospective DR studies.

Leishmania, a unicellular eukaryotic parasite, is a unique model for aneuploidy and cellular heterogeneity, along with their potential role in adaptation to environmental stresses. Somy variation within clonal populations was previously explored in a small subset of chromosomes using fluorescence hybridization methods. This phenomenon, termed mosaic aneuploidy (MA), might have important evolutionary and functional implications but remains under-explored due to technological limitations. Here, we applied and validated a high throughput single-cell genome sequencing method to study for the first time the extent and dynamics of whole karyotype heterogeneity in two clonal populations of Leishmania promastigotes representing different stages of MA evolution in vitro. We found that drastic changes in karyotypes quickly emerge in a population stemming from an almost euploid founder cell. This possibly involves polyploidization/hybridization at an early stage of population expansion, followed by assorted ploidy reduction. During further stages of expansion, MA increases by moderate and gradual karyotypic alterations, affecting a defined subset of chromosomes. Our data provide the first complete characterization of MA in Leishmania and pave the way for further functional studies., (© The Author(s) 2021. Published by Oxford University Press on behalf of Nucleic Acids Research.)

Microorganisms can adopt a quiescent physiological condition which acts as a survival strategy under unfavorable conditions. Quiescent cells are characterized by slow or non-proliferation and a deep downregulation of processes related to biosynthesis. Although quiescence has been described mostly in bacteria, this survival skill is widespread, including in eukaryotic microorganisms. In Leishmania , a digenetic parasitic protozoan that causes a major infectious disease, quiescence has been demonstrated, but the molecular and metabolic features enabling its maintenance are unknown. Here, we quantified the transcriptome and metabolome of Leishmania promastigotes and amastigotes where quiescence was induced in vitro either, through drug pressure or by stationary phase. Quiescent cells have a global and coordinated reduction in overall transcription, with levels dropping to as low as 0.4% of those in proliferating cells. However, a subset of transcripts did not follow this trend and were relatively upregulated in quiescent populations, including those encoding membrane components, such as amastins and GP63, or processes like autophagy. The metabolome followed a similar trend of overall downregulation albeit to a lesser magnitude than the transcriptome. It is noteworthy that among the commonly upregulated metabolites were those involved in carbon sources as an alternative to glucose. This first integrated two omics layers afford novel insight into cell regulation and show commonly modulated features across stimuli and stages.

The extent of susceptibility towards miltefosine (Mil), amphotericin B (AmpB), and paromomycin (Paro) was measured among 19 clinical isolates of Leishmania donovani (LD). Thirteen of these clinical isolates were reported to exhibit low susceptibility towards sodium stibogluconate (SSG-R), while six of them were highly susceptible (SSG-S). The degree of clearance of amastigotes (EC50) for these predefined SSG-R- and SSG-S-infected macrophages was determined against Mil, AmpB, and Paro. Two out of the 13 SSG-R isolates (BHU575 and BHU814) showed low susceptibility towards all three drugs studied, while the rest of the 11 SSG-R isolates showed varying degrees of susceptibility either towards none or only towards individual drugs. Interestingly, all the SSG-S isolates showed high susceptibility towards Mil/AmpB/Paro. The total intracellular non-protein thiol content of the LD promastigotes, which have been previously reported to be positively co-related with EC50 towards SSG, was found to be independent from the degree of susceptibility towards Mil/AmpB/Paro. Impedance spectra analysis, which quantifies membrane resistance, revealed lower impedimetric values for all those isolates exhibiting low efficacy to Mil (Mil-R). Our analysis points out that while non-protein thiol content can be an attribute of SSG-R, lower impedimetric values can be linked with lower Mil susceptibility, although neither of these parameters seems to get influenced by the degree of susceptibility towards AmpB/Paro. Finally, a correlation analysis with established biological methods suggests that impedance spectral analysis can be used for the accurate determination of lower Mil susceptibility among LD isolates, which is further validated in the LD-infected in vivo hamster model., Competing Interests: The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest., (Copyright © 2021 Ghosh, Biswas, Mukherjee, Pal, Saxena, Sundar, Dujardin, Das, Roy, Mukhopadhyay and Mukherjee.)

Background: The onchocerciasis focus surrounding the lower Mbam and Sanaga rivers, where Onchocerca volvulus is transmitted by Simulium damnosum s.l. (Diptera: Simuliidae), was historically the largest in the southern regions of Cameroon. Annual community-directed treatment with ivermectin (CDTI) has been taking place since 2000, but recent studies have shown that new infections are occurring in children. We aimed to investigate blackfly biting and O. volvulus transmission rates along the lower Mbam river 16 years after the formal onset of annual CDTI., Methods: Black flies were collected for three consecutive days each month between July 2016 and June 2017 at two riverside villages and two inland sites situated 4.9 km and 7.9 km from the riverside. Specimens collected at each site were dissected on one of the three collection days each month to estimate parity rates and O. volvulus infection rates, while the remaining samples were preserved for pool screening., Results: In total, 93,573 S. damnosum s.l. black flies were recorded biting humans and 9281 were dissected. Annual biting rates of up to 606,370 were estimated at the riverside, decreasing to 20,540 at 7.9 km, while, based on dissections, annual transmission potentials of up to 4488 were estimated at the riverside, decreasing to 102 and 0 at 4.9 km and 7.9 km, respectively. However, pool screening showed evidence of infection in black flies at the furthest distance from the river. Results of both methods demonstrated the percentage of infective flies to be relatively low (0.10-0.36%), but above the WHO threshold for interruption of transmission. In addition, a small number of larvae collected during the dry season revealed the presence of Simulium squamosum E. This is the first time S. squamosum E has been found east of Lake Volta in Ghana, but our material was chromosomally distinctive, and we call it S. squamosum E2., Conclusions: Relatively low O. volvulus infection rates appear to be offset by extremely high densities of biting black flies which are sustaining transmission along the banks of the lower Mbam river., (© 2021. The Author(s).)

Leishmaniasis, a disease caused by protozoa of the Leishmania species, afflicts roughly 12 million individuals worldwide. Most existing drugs for leishmaniasis are toxic, expensive, difficult to administer, and subject to drug resistance. We report a new class of antileishmanial leads, the 3-arylquinolines, that potently block proliferation of the intramacrophage amastigote form of Leishmania parasites with good selectivity relative to the host macrophages. Early lead 34 was rapidly acting and possessed good potency against L. mexicana (EC 50 = 120 nM), 30-fold selectivity for the parasite relative to the macrophage (EC 50 = 3.7 μM), and also blocked proliferation of Leishmania donovani parasites resistant to antimonial drugs. Finally, another early lead, 27 , which exhibited reasonable in vivo tolerability, impaired disease progression during the dosing period in a murine model of cutaneous leishmaniasis. These results suggest that the arylquinolines provide a fruitful departure point for the development of new antileishmanial drugs.

The tropical Andes are an important natural laboratory to understand speciation in many taxa. Here we examined the evolutionary history of parasites of the Leishmania braziliensis species complex based on whole-genome sequencing of 67 isolates from 47 localities in Peru. We first show the origin of Andean Leishmania as a clade of near-clonal lineages that diverged from admixed Amazonian ancestors, accompanied by a significant reduction in genome diversity and large structural variations implicated in host-parasite interactions. Within the Andean species, patterns of population structure were strongly associated with biogeographical origin. Molecular clock and ecological niche modeling suggested that the history of diversification of the Andean lineages is limited to the Late Pleistocene and intimately associated with habitat contractions driven by climate change. These results suggest that changes in forestation over the past 150,000 y have influenced speciation and diversity of these Neotropical parasites. Second, genome-scale analyses provided evidence of meiotic-like recombination between Andean and Amazonian Leishmania species, resulting in full-genome hybrids. The mitochondrial genome of these hybrids consisted of homogeneous uniparental maxicircles, but minicircles originated from both parental species. We further show that mitochondrial minicircles-but not maxicircles-show a similar evolutionary pattern to the nuclear genome, suggesting that compatibility between nuclear-encoded mitochondrial genes and minicircle-encoded guide RNA genes is essential to maintain efficient respiration. By comparing full nuclear and mitochondrial genome ancestries, our data expand our appreciation on the genetic consequences of diversification and hybridization in parasitic protozoa., Competing Interests: The authors declare no competing interest., (Copyright © 2020 the Author(s). Published by PNAS.)

The protozoan parasite Leishmania ( Viannia ) braziliensis (L. braziliensis ) is the main cause of human tegumentary leishmaniasis in the New World, a disease affecting the skin and/or mucosal tissues. Despite its importance, the study of the unique biology of L. braziliensis through reverse genetics analyses has so far lagged behind in comparison with Old World Leishmania spp. In this study, we successfully applied a cloning-free, PCR-based CRISPR-Cas9 technology in L. braziliensis that was previously developed for Old World Leishmania major and New World L. mexicana species. As proof of principle, we demonstrate the targeted replacement of a transgene ( eGFP ) and two L. braziliensis single-copy genes ( HSP23 and HSP100 ). We obtained homozygous Cas9-free HSP23 - and HSP100 -null mutants in L. braziliensis that matched the phenotypes reported previously for the respective L. donovani null mutants. The function of HSP23 is indeed conserved throughout the Trypanosomatida as L. major HSP23 null mutants could be complemented phenotypically with transgenes from a range of trypanosomatids. In summary, the feasibility of genetic manipulation of L. braziliensis by CRISPR-Cas9-mediated gene editing sets the stage for testing the role of specific genes in that parasite's biology, including functional studies of virulence factors in relevant animal models to reveal novel therapeutic targets to combat American tegumentary leishmaniasis.

Infections caused by protozoan parasites burden the world with huge costs in terms of human and animal health. Most parasitic diseases caused by protozoans are neglected, particularly those associated with poverty and tropical countries, but the paucity of drug treatments and vaccines combined with increasing problems of drug resistance are becoming major concerns for their control and eradication. In this climate, the discovery/repurposing of new drugs and increasing effort in vaccine development should be supplemented with an exploration of new alternative/synergic treatment strategies. Viruses, either native or engineered, have been employed successfully as highly effective and selective therapeutic approaches to treat cancer (oncolytic viruses) and antibiotic-resistant bacterial diseases (phage therapy). Increasing evidence is accumulating that many protozoan, but also helminth, parasites harbour a range of different classes of viruses that are mostly absent from humans. Although some of these viruses appear to have no effect on their parasite hosts, others either have a clear direct negative impact on the parasite or may, in fact, contribute to the virulence of parasites for humans. This review will focus mainly on the viruses identified in protozoan parasites that are of medical importance. Inspired and informed by the experience gained from the application of oncolytic virus- and phage-therapy, rationally-driven strategies to employ these viruses successfully against parasitic diseases will be presented and discussed in the light of the current knowledge of the virus biology and the complex interplay between the viruses, the parasite hosts and the human host. We also highlight knowledge gaps that should be addressed to advance the potential of virotherapy against parasitic diseases.

Here, we report a pilot study paving the way for further single cell genomics studies in Leishmania. First, the performances of two commercially available kits for Whole Genome Amplification (WGA), PicoPLEX and RepliG were compared on small amounts of Leishmania donovani DNA, testing their ability to preserve specific genetic variations, including aneuploidy levels and SNPs. We show here that the choice of WGA method should be determined by the planned downstream genetic analysis, PicoPLEX and RepliG performing better for aneuploidy and SNP calling, respectively. This comparison allowed us to evaluate and optimize corresponding bio-informatic methods. As PicoPLEX was shown to be the preferred method for studying single cell aneuploidy, this method was applied in a second step, on single cells of L. braziliensis, which were sorted by fluorescence activated cell sorting (FACS). Even sequencing depth was achieved in 28 single cells, allowing accurate somy estimation. A dominant karyotype with three aneuploid chromosomes was observed in 25 cells, while two different minor karyotypes were observed in the other cells. Our method thus allowed the detection of aneuploidy mosaicism, and provides a solid basis which can be further refined to concur with higher-throughput single cell genomic methods.

Cytosine C5 methylation is an important epigenetic control mechanism in a wide array of eukaryotic organisms and generally carried out by proteins of the C-5 DNA methyltransferase family (DNMTs). In several protozoans, the status of this mechanism remains elusive, such as in Leishmania , the causative agent of the disease leishmaniasis in humans and a wide array of vertebrate animals. In this work, we showed that the Leishmania donovani genome contains a C-5 DNA methyltransferase ( DNMT ) from the DNMT6 subfamily, whose function is still unclear, and verified its expression at the RNA level. We created viable overexpressor and knock-out lines of this enzyme and characterized their genome-wide methylation patterns using whole-genome bisulfite sequencing, together with promastigote and amastigote control lines. Interestingly, despite the DNMT6 presence, we found that methylation levels were equal to or lower than 0.0003% at CpG sites, 0.0005% at CHG sites, and 0.0126% at CHH sites at the genomic scale. As none of the methylated sites were retained after manual verification, we conclude that there is no evidence for DNA methylation in this species. We demonstrated that this difference in DNA methylation between the parasite (no detectable DNA methylation) and the vertebrate host (DNA methylation) allowed enrichment of parasite vs. host DNA using methyl-CpG-binding domain columns, readily available in commercial kits. As such, we depleted methylated DNA from mixes of Leishmania promastigote and amastigote DNA with human DNA, resulting in average Leishmania: human enrichments from 62× up to 263×. These results open a promising avenue for unmethylated DNA enrichment as a pre-enrichment step before sequencing Leishmania clinical samples.

Background: Endemic regions of cutaneous leishmaniasis (CL) and intestinal helminthiasis overlap. CL treatment with systemic pentavalent antimonial drugs (Sb 5+ ) fails in 10%-30% of patients. The study objective was to assess the etiological role of intestinal helminthiasis in CL treatment failure., Methods: An unmatched case-control study was done in 4 CL treatment sites in Peru in 2012-2015. Cases were CL patients with Sb 5+ treatment failure; controls were CL patients with Sb 5+ treatment success. Patients with a parasitologically confirmed CL diagnosis who had received supervised Sb 5+ treatment and could be classified as cases or controls were eligible. The main exposure variables were intestinal helminthiasis and strongyloidiasis, diagnosed through direct examination, rapid sedimentation, Baermann, Kato-Katz, or agar culture of stool samples. Additional exposure variables were other infections (HIV, human T-lymphotropic virus 1, tuberculosis, hepatitis B, intestinal protozoa) and noninfectious conditions (diabetes, renal insufficiency, and immunosuppressive medication). Age, gender, CL history, probable exposure place, and Leishmania species were treated as potential confounders in multiple logistic regression., Results: There were 94 case and 122 control subjects. Overall, infectious and noninfectious comorbidities were frequent both among cases (64%) and controls (71%). The adjusted odds ratio (OR) for the association between any intestinal helminth infection and CL treatment failure was 0.65 (95% confidence interval [CI], 0.30-1.38), and the adjusted OR for the association between strongyloidiasis and CL treatment failure was 0.34 (95% CI, 0.11-0.92)., Conclusions: In the Peruvian setting, high Sb 5+ treatment failure rates are not explained by intestinal helminthiasis. On the contrary, strongyloidiasis had a protective effect against treatment failure., (© The Author(s) 2020. Published by Oxford University Press on behalf of Infectious Diseases Society of America.)

Elimination programs targeting TriTryp diseases (Leishmaniasis, Chagas' disease, human African trypanosomiasis) significantly reduced the number of cases. Continued surveillance is crucial to sustain this progress, but parasite molecular surveillance by genotyping is currently lacking. We explain here which epidemiological questions of public health and clinical relevance could be answered by means of molecular surveillance. Whole-genome sequencing (WGS) for molecular surveillance will be an important added value, where we advocate that preference should be given to direct sequencing of the parasite's genome in host tissues instead of analysis of cultivated isolates. The main challenges here, and recent technological advances, are discussed. We conclude with a series of recommendations for implementing whole-genome sequencing for molecular surveillance., (Copyright © 2020 Elsevier Ltd. All rights reserved.)

Protozoan parasites of the Leishmania donovani complex - L. donovani and L. infantum - cause the fatal disease visceral leishmaniasis. We present the first comprehensive genome-wide global study, with 151 cultured field isolates representing most of the geographical distribution. L. donovani isolates separated into five groups that largely coincide with geographical origin but vary greatly in diversity. In contrast, the majority of L. infantum samples fell into one globally-distributed group with little diversity. This picture is complicated by several hybrid lineages. Identified genetic groups vary in heterozygosity and levels of linkage, suggesting different recombination histories. We characterise chromosome-specific patterns of aneuploidy and identified extensive structural variation, including known and suspected drug resistance loci. This study reveals greater genetic diversity than suggested by geographically-focused studies, provides a resource of genomic variation for future work and sets the scene for a new understanding of the evolution and genetics of the Leishmania donovani complex., Competing Interests: SF, CD, OS, BM, TD, HI, JD, MS, IM, MM, LS, CJ, AN, HS, MY, TB, MA, MB, TW, GS, JC No competing interests declared, (© 2020, Franssen et al.)

We identified visceral leishmaniasis caused by Leishmania donovani in a previously unknown focus in northern Somalia. Clinical and epidemiologic characteristics of 118 cases during 2013-2019 in Bosaso, the region's commercial capital, have raised suspicion of visceral leishmaniasis endemicity status there.

Understanding the mechanism(s) underpinning drug resistance could lead to novel treatments to reverse the increased tolerance of a pathogen. In this study, paromomycin (PMM) resistance (PMM r ) was induced in three Nepalese clinical strains of Leishmania donovani with different inherent susceptibilities to antimony (Sb) drugs by stepwise exposure of promastigotes to PMM. Exposure to PMM resulted in the production of mixed populations of parasites, even though a single cloned population was used at the start of selection. PMM 50% inhibitory concentration (IC 50 ) values for PMM r parasites varied between 104 and 481 μM at the promastigote stage and 32 and 195 μM at the intracellular amastigote stage. PMM resistance was associated with increased resistance to nitric oxide at the amastigote stage but not the promastigote stage ( P  < 0.05). This effect was most marked in the Sb-resistant (Sb r ) PMM r clone, in which PMM resistance was associated with a significant upregulation of glutathione compared to that in its wild type ( P  < 0.05), although there was no change in the regulation of trypanothione (detected in its oxidized form). Interestingly, PMM r strains showed an increase in either the keto acid derivative of isoleucine (Sb intermediate PMM r ) or the 2-hydroxy acids derived from arginine and tyrosine (Sb susceptible PMM r and Sb r PMM r ). These results are consistent with the recent finding that the upregulation of the branched-chain amino acid aminotransferase and d-lactate dehydrogenase is linked to PMM r In addition, we found that PMM r is associated with a significant increase in aneuploidy during PMM selection in all the strains, which could allow the rapid selection of genetic changes that confer a survival advantage., (Copyright © 2019 Shaw et al.)

Under stressful conditions some microorganisms adopt a quiescent stage characterized by a reversible non or slow proliferative condition that allows their survival. This adaptation was only recently discovered in Leishmania. We developed an in vitro model and a biosensor to track quiescence at population and single cell levels. The biosensor is a GFP reporter gene integrated within the 18S rDNA locus, which allows monitoring the expression of 18S rRNA (rGFP expression). We showed that rGFP expression decreased significantly and rapidly during the transition from extracellular promastigotes to intracellular amastigotes and that it was coupled in vitro with a decrease in replication as measured by BrdU incorporation. rGFP expression was useful to track the reversibility of quiescence in live cells and showed for the first time the heterogeneity of physiological stages among the population of amastigotes in which shallow and deep quiescent stages may coexist. We also validated the use of rGFP expression as a biosensor in animal models of latent infection. Our models and biosensor should allow further characterization of quiescence at metabolic and molecular level.

Whole genome sequencing (WGS) is increasingly used for molecular diagnosis and epidemiology of infectious diseases. Current Leishmania genomic studies rely on DNA extracted from cultured parasites, which might introduce sampling and biological biases into the subsequent analyses. Up to now, direct analysis of Leishmania genome in clinical samples is hampered by high levels of human DNA and large variation in parasite load in clinical samples. Here, we present a method, based on target enrichment of Leishmania donovani DNA with Agilent SureSelect technology, that allows the analysis of Leishmania genomes directly in clinical samples. We validated our protocol with a set of artificially mixed samples, followed by the analysis of 63 clinical samples (bone marrow or spleen aspirates) from visceral leishmaniasis patients in Nepal. We were able to identify genotypes using a set of diagnostic SNPs in almost all of these samples (97%) and access comprehensive genome-wide information in most (83%). This allowed us to perform phylogenomic analysis, assess chromosome copy number and identify large copy number variants (CNVs). Pairwise comparisons between the parasite genomes in clinical samples and derived in vitro cultured promastigotes showed a lower aneuploidy in amastigotes as well as genomic differences, suggesting polyclonal infections in patients. Altogether our results underline the need for sequencing parasite genomes directly in the host samples., Competing Interests: The authors have declared that no competing interests exist.

Visceral leishmaniasis (VL), the most severe form of the disease, is caused by Leishmania donovani in the Indian sub-continent (ISC). Whole genome sequencing studies revealed that two parasite populations exist in the ISC: a main population named the Core Group (CG) found mostly in the lowlands, and a new, genetically different subpopulation called ISC1. Parasites belonging to the CG were shown to be responsible for the recent epidemics, while the ISC1 variant was originally identified in hilly districts of Nepal and was later on increasingly found in the lowlands. Importantly, the ISC1 and CG isolates differ in their drug susceptibility and virulence signatures, suggesting that ISC1 constitutes an emerging and functionally different variant of L. donovani. In present study we aimed to address the potential of ISC1 transmission by the natural vector of L. donovani in the lowlands, Phlebotomus argentipes. By experimental infection of sand flies with parasites of the different genotypes, we demonstrate that ISC1 and CG strains are developing similarly in P. argentipes, suggesting that P. argentipes is a fully competent vector for ISC1 parasites. Integration of previous and current findings shows thus that ISC1 is a new and different variant of L. donovani, fully adapted to spread in the ISC through the main vector. This information is directly useful for managers of the elimination program. Furthermore, integration of our successive studies (genotyping, phenotyping and vector competence) demonstrates the relevance of molecular surveillance and should be of interest for scientists working on vector borne diseases and control managers., (Copyright © 2019 Elsevier B.V. All rights reserved.)

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

Background: The search for diagnostic biomarkers has been profiting from a growing number of high quality sequenced genomes and freely available bioinformatic tools. These can be combined with wet lab experiments for a rational search. Improved, point-of-care diagnostic tests for visceral leishmaniasis (VL), early case detection and surveillance are required. Previous investigations demonstrated the potential of IgG1 as a biomarker for monitoring clinical status in rapid diagnostic tests (RDTs), although using a crude lysate antigen (CLA) as capturing antigen. Replacing the CLA by specific antigens would lead to more robust RDTs., Methodology: Immunoblots revealed L. donovani protein bands detected by IgG1 from VL patients. Upon confident identification of these antigens by mass spectrometry (MS), we searched for evidence of constitutive protein expression and presence of antigenic domains or high accessibility to B-cells. Selected candidates had their linear epitopes mapped with in silico algorithms. Multiple high-scoring predicted epitopes from the shortlisted proteins were screened in peptide arrays. The most promising candidate was tested in RDT prototypes using VL and nonendemic healthy control (NEHC) patient sera., Results: Over 90% of the proteins identified from the immunoblots did not satisfy the selection criteria and were excluded from the downstream epitope mapping. Screening of predicted epitope peptides from the shortlisted proteins identified the most reactive, for which the sensitivity for IgG1 was 84% (95% CI 60-97%) with Sudanese VL sera on RDT prototypes. None of the sera from NEHCs were positive., Conclusion: We employed in silico searches to reduce drastically the output of wet lab experiments, focusing on promising candidates containing selected protein features. By predicting epitopes in silico we screened a large number of peptides using arrays, identifying the most promising one, for which IgG1 sensitivity and specificity, with limited sample size, supported this proof of concept strategy for diagnostics discovery, which can be applied to the development of more robust IgG1 RDTs for monitoring clinical status in VL., Competing Interests: BCBH was and PM is employee of Coris BioConcept and they do not have any share in the company. Sergey Litvinov is an employee of Aptum Biologics Ltd. which is an official partner institution in the Consortium that produced this work. Aptum Biologics Ltd. has and provided to the Consortium the software for epitope prediction, and beyond this has no consultancy, employment or any other forms of financial benefits from this work, nor affiliation with whatever potential intellectual property that may arise from this particular work. All other authors have declared that no competing interests exist.

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

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

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

In rural areas in Morocco, diagnosing cutaneous leishmaniasis (CL) can be challenging. We evaluated the accuracy of a rapid diagnostic test (RDT) based on antigen detection, CL Detect Rapid Test ™ (Inbios International Inc., Seattle, WA), in this setting. We consecutively recruited patients with new skin ulcers in nine primary health centers. We took a dental broach sample for the RDT and two other tissue samples by scraping the border and center of the lesion with a scalpel and smearing it on a slide. We duplicated each smear by pressing a clean slide against it and processed the slides by microscopy, polymerase chain reaction (PCR) internal transcribed spacer 1, and kDNA minicircle PCR. In a subgroup with positive PCR, the Leishmania species was identified using PCR-restriction fragment length polymorphism and PCR-sequencing of hsp70 genes. A participant with positive microscopy and/or PCR was considered a confirmed CL case. We computed sensitivity (Se) and specificity (Sp) of the RDT compared with this reference standard (ClinicalTrials.gov registration: NCT02979002). Between December 2016 and July 2017, we included 219 patients, 50% of them were under 18 years old. Rapid diagnostic test Se was 68% [95% confidence interval (CI): 61-74], Sp 94% [95% CI: 91-97], positive predictive value 95% [95% CI: 92-98], and negative predictive value 64% [95% CI: 58-70]. Despite its low Se, this novel RDT is a useful addition to clinical management of CL in Morocco, especially in isolated localities. Rapid diagnostic test-positive lesions can be treated as CL; but when RDT negative, microscopy should be done in a second step. The Se of the RDT can probably be optimized by improving the sampling procedure.

We recently published the first genomic diversity study of Trypanosoma congolense, a major aetiological agent of Animal African Trypanosomiasis. We demonstrated striking levels of SNP and indel diversity in the Eastern province of Zambia as a consequence of hybridization between divergent trypanosome lineages. We concluded that these and earlier findings in T. congolense challenge the predominant clonal evolution (PCE) model. In a recent comment, Tibayrenc and Ayala claim that there are many features in T. congolense supporting their theory of clonality. While we can follow the reasoning of the authors, we also identify major limitations in their theory and interpretations that resulted in incorrect conclusions. First, we argue that each T. congolense subgroup should be analysed independently as they may represent different (sub)species rather than "near-clades". Second, the authors neglect major findings of two robust population genetic studies on Savannah T. congolense that provide clear evidence of frequent recombination. Third, we reveal additional events of introgressive hybridization in T. congolense by analysing the maxicircle coding region using next-generation sequencing analyses. At last, we pinpoint two important misinterpretations by the authors and show that there are no spatially and temporally widespread clones in T. congolense. We stand by our earlier conclusions that the clonal framework is unlikely to accurately model the population structure of T. congolense. Other theoretical frameworks such as Maynard Smith's epidemic model may better represent the complex ancestry seen in T. congolense, where clones delimited in space and time arise against a background of recombination., (© 2018 John Wiley & Sons Ltd.)

The growing drug resistance (DR) raises major concerns for the control of visceral leishmaniasis (VL), a neglected disease lethal in 95 percent of the cases if left untreated. Resistance has rendered antimonials (SSG) obsolete in the Indian Sub-Continent (ISC) and the first miltefosine-resistant Leishmania donovani were isolated. New chemotherapeutic options are needed and novel compounds are being identified by high-throughput screening (HTS). HTS is generally performed with old laboratory strains such as LdBOB and we aimed here to validate the activity of selected compounds against recent clinical isolates. In this academic/industrial collaboration, 130 compounds from the GSK "Leishbox" were screened against one SSG-sensitive and one SSG-resistant strain of L. donovani recently isolated from ISC patients, using an intracellular assay of L. donovani-infected THP1-derived macrophages. We showed that only 45% of the compounds were active in both clinical isolates and LdBOB. There were also different compound efficiencies linked to the SSG susceptibility background of the strains. In addition, our results suggested that the differential susceptibility profiles were chemical series-dependent. In conclusion, we demonstrate the potential value of including clinical isolates (as well as resistant strains) in the HTS progression cascade.

Leishmania donovani is the responsible agent for visceral leishmaniasis (VL) in the Indian subcontinent (ISC). The disease is lethal without treatment and causes 0.2 to 0.4 million cases each year. Recently, reports of VL in Nepalese hilly districts have increased as well as VL cases caused by L. donovani from the ISC1 genetic group, a new and emerging genotype. In this study, we perform for the first time an integrated, untargeted genomics and metabolomics approach to characterize ISC1, in comparison with the Core Group (CG), main population that drove the most recent outbreak of VL in the ISC. We show that the ISC1 population is very different from the CG, both at genome and metabolome levels. The genomic differences include SNPs, CNV and small indels in genes coding for known virulence factors, immunogens and surface proteins. Both genomic and metabolic approaches highlighted dissimilarities related to membrane lipids, the nucleotide salvage pathway and the urea cycle in ISC1 versus CG. Many of these pathways and molecules are important for the interaction with the host/extracellular environment. Altogether, our data predict major functional differences in ISC1 versus CG parasites, including virulence. Therefore, particular attention is required to monitor the fate of this emerging ISC1 population in the ISC, especially in a post-VL elimination context., (Copyright © 2018 The Authors. Published by Elsevier B.V. All rights reserved.)

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

The Mahenge Mountains onchocerciasis focus in south eastern Tanzania was historically one of the most heavily infected areas in the country. The vectors of Onchocerca volvulus are mainly Simulium damnosum complex blackflies, but a species of the Simulium neavei group may also contribute to transmission in some areas. The only detailed studies of parasite transmission in Mahenge were conducted in the late 1960s. The taxonomy of the S. damnosum complex has since been revised and onchocerciasis control through annual community directed treatment with ivermectin (CDTI) commenced in 1997. This study aimed to provide a cytogenetic and molecular update of the S. damnosum complex cytoforms present in Mahenge, and to evaluate the current status of O. volvulus transmission by blackflies following 19 years of annual CDTI. Rivers were surveyed to identify sites of S. damnosum s.l. breeding among the eastern slopes of the mountains, and human landing collections of adult female blackflies were made close to breeding sites. Identification of S. damnosum complex cytoforms was by cytotaxonomy of late-instar larvae and ITS1 amplicon size polymorphisms of larvae and adults. Adult blackflies were pool screened for O. volvulus infection using a triplex real-time PCR. The cytoforms 'Nkusi', Simulium kilibanum and 'Turiani' were found breeding in perennial rivers. 'Nkusi' and S. kilibanum were collected on human bait at 7/7 catch sites and possessed ITS1 profiles most closely resembling the molecular forms 'Nkusi J' and S. kilibanum 'T'. Whereas 'Turiani' was present in rivers, it was not collected on human bait and appears to be zoophilic. Simulium nyasalandicum was collected in low numbers on human bait at 3/7 catch sites. In total, 12,452 S. damnosum s.l. were pool screened and O. volvulus infection was detected in 97/104 pools of bodies and 51/104 pools of heads. The estimated percentage of S. damnosum s.l. carrying infective L3 stage parasites was 0.57% (95% CI 0.43%-0.74%). Onchocerca volvulus transmission by S. damnosum s.l. is continuing in the Mahenge Mountains after 19 years of annual CDTI. Infection rates appear similar to those reported in the 1960s, but a more detailed study is required to fully understand the epidemiological significance of the ongoing transmission. These results provide further evidence that annual CDTI may be insufficient to eliminate the parasite in formerly hyperendemic foci., (Copyright © 2018 The Authors. Published by Elsevier B.V. All rights reserved.)

Antimonials (Sb) were used for decades for chemotherapy of visceral leishmaniasis (VL). Now abandoned in the Indian subcontinent (ISC) because of Leishmania donovani resistance, this drug offers a unique model for understanding drug resistance dynamics. In a previous phylogenomic study, we found two distinct populations of L. donovani : the core group (CG) in the Gangetic plains and ISC1 in the Nepalese highlands. Sb resistance was only encountered within the CG, and a series of potential markers were identified. Here, we analyzed the development of resistance to trivalent antimonials (Sb III ) upon experimental selection in ISC1 and CG strains. We observed that (i) baseline Sb III susceptibility of parasites was higher in ISC1 than in the CG, (ii) time to Sb III resistance was higher for ISC1 parasites than for CG strains, and (iii) untargeted genomic and metabolomic analyses revealed molecular changes along the selection process: these were more numerous in ISC1 than in the CG. Altogether these observations led to the hypothesis that CG parasites are preadapted to Sb III resistance. This hypothesis was experimentally confirmed by showing that only wild-type CG strains could survive a direct exposure to the maximal concentration of Sb III The main driver of this preadaptation was shown to be MRPA , a gene involved in Sb III sequestration and amplified in an intrachromosomal amplicon in all CG strains characterized so far. This amplicon emerged around 1850 in the CG, well before the implementation of antimonials for VL chemotherapy, and we discuss here several hypotheses of selective pressure that could have accompanied its emergence. IMPORTANCE The "antibiotic resistance crisis" is a major challenge for scientists and medical professionals. This steady rise in drug-resistant pathogens also extends to parasitic diseases, with antimony being the first anti- Leishmania drug that fell in the Indian subcontinent (ISC). Leishmaniasis is a major but neglected infectious disease with limited therapeutic options. Therefore, understanding how parasites became resistant to antimonials is of commanding importance. In this study, we experimentally characterized the dynamics of this resistance acquisition and show for the first time that some Leishmania populations of the ISC were preadapted to antimony resistance, likely driven by environmental factors or by drugs used in the 19th century., (Copyright © 2018 Dumetz et al.)

Background: Tegumentary leishmaniasis (TL) is a disease of skin and/or mucosal tissues caused by Leishmania parasites. TL patients may concurrently carry other pathogens, which may influence the clinical outcome of TL., Methodology and Principal Findings: This review focuses on the frequency of TL coinfections in human populations, interactions between Leishmania and other pathogens in animal models and human subjects, and implications of TL coinfections for clinical practice. For the purpose of this review, TL is defined as all forms of cutaneous (localised, disseminated, or diffuse) and mucocutaneous leishmaniasis. Human immunodeficiency virus (HIV) coinfection, superinfection with skin bacteria, and skin manifestations of visceral leishmaniasis are not included. We searched MEDLINE and other databases and included 73 records: 21 experimental studies in animals and 52 studies about human subjects (mainly cross-sectional and case studies). Several reports describe the frequency of Trypanosoma cruzi coinfection in TL patients in Argentina (about 41%) and the frequency of helminthiasis in TL patients in Brazil (15% to 88%). Different hypotheses have been explored about mechanisms of interaction between different microorganisms, but no clear answers emerge. Such interactions may involve innate immunity coupled with regulatory networks that affect quality and quantity of acquired immune responses. Diagnostic problems may occur when concurrent infections cause similar lesions (e.g., TL and leprosy), when different pathogens are present in the same lesions (e.g., Leishmania and Sporothrix schenckii), or when similarities between phylogenetically close pathogens affect accuracy of diagnostic tests (e.g., serology for leishmaniasis and Chagas disease). Some coinfections (e.g., helminthiasis) appear to reduce the effectiveness of antileishmanial treatment, and drug combinations may cause cumulative adverse effects., Conclusions and Significance: In patients with TL, coinfection is frequent, it can lead to diagnostic errors and delays, and it can influence the effectiveness and safety of treatment. More research is needed to unravel how coinfections interfere with the pathogenesis of TL.

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

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