Search

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. ispartof: NUCLEIC ACIDS RESEARCH vol:50 issue:1 pages:293-305 ispartof: location:England status: published

The implementation of prospective drug resistance (DR) studies in the R&D pipelines is a common practice for many infectious diseases, but not for Neglected Tropical Diseases. 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 GSK ‘Leishbox’ to treat VL. We experimentally selected resistance to TCMDC-143345 in vitro and characterized resistant parasites at 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 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/resistance of L. donovani to TCMDC-143345.ImportanceHumans 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.

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, Author summary Visceral leishmaniasis (VL) is caused by parasitic protozoa of the Leishmania donovani complex and is lethal in the absence of treatment. Whole Genome Sequencing (WGS) of L. donovani clinical isolates revealed hitherto cryptic population structure in the Indian Sub-Continent and provided insights into the epidemiology and potential mechanisms of drug resistance. However, several biases are likely introduced during the culture step. We report here the development of a method that allows determination of parasite genomes directly in clinical samples, and validate it on bone marrow and splenic aspirates of VL patients in Nepal. Our study sheds a new light on the biology of Leishmania in the human host: we found that intracellular parasites of the patients had very low levels of aneuploidy, in sharp contrast to the situation in cultivated isolates. Moreover, the observed differences in genomes between intracellular amastigotes of the patient and the derived cultured parasites suggests polyclonality of infections, with different clones dominating in clinical samples and in culture, likely due to fitness differences. We believe this method is most suitable for clinical studies and for molecular tracking in the context of elimination programs.

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., Highlights • ISC1 is an emerging Leishmania donovani population in the Nepalese highlands • The genome and metabolome of ISC1 are distinct from strains from the most recent visceral leishmaniasis outbreak in the Indian sub-continent • These differences include important host interaction factors • The emergence of ISC1 should be closely monitored in the Indian subcontinent

The “antibiotic resistance crisis” is a major challenge for scientists and medical professionals. This steady rise in drug-resistant pathogens also extends to parasitic diseases, with antimony being the first anti-Leishmania drug that fell in the Indian subcontinent (ISC). Leishmaniasis is a major but neglected infectious disease with limited therapeutic options. Therefore, understanding how parasites became resistant to antimonials is of commanding importance. In this study, we experimentally characterized the dynamics of this resistance acquisition and show for the first time that some Leishmania populations of the ISC were preadapted to antimony resistance, likely driven by environmental factors or by drugs used in the 19th century., Antimonials (Sb) were used for decades for chemotherapy of visceral leishmaniasis (VL). Now abandoned in the Indian subcontinent (ISC) because of Leishmania donovani resistance, this drug offers a unique model for understanding drug resistance dynamics. In a previous phylogenomic study, we found two distinct populations of L. donovani: the core group (CG) in the Gangetic plains and ISC1 in the Nepalese highlands. Sb resistance was only encountered within the CG, and a series of potential markers were identified. Here, we analyzed the development of resistance to trivalent antimonials (SbIII) upon experimental selection in ISC1 and CG strains. We observed that (i) baseline SbIII susceptibility of parasites was higher in ISC1 than in the CG, (ii) time to SbIII 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 SbIII resistance. This hypothesis was experimentally confirmed by showing that only wild-type CG strains could survive a direct exposure to the maximal concentration of SbIII. The main driver of this preadaptation was shown to be MRPA, a gene involved in SbIII 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.

Leishmania donovani causes visceral leishmaniasis (VL), the second most deadly vector- borne parasitic disease. A recent epidemic in the Indian subcontinent (ISC) caused up to 80% of global VL and over 30,000 deaths per year. Resistance against antimonial drugs has probably been a contributing factor in the persistence of this epidemic. Here we use whole genome sequences from 204 clinical isolates to track the evolution and epidemiology of L. donovani from the ISC. We identify independent radiations that have emerged since a bottleneck coincident with 1960s DDT spraying campaigns. A genetically distinct population frequently resistant to antimonials has a two base-pair insertion in the aquaglyceroporin gene LdAQP1 that prevents the transport of trivalent antimonials. We find evidence of genetic exchange between ISC populations, and show that the mutation in LdAQP1 has spread by recombination. Our results reveal the complexity of L. donovani evolution in the ISC in response to drug treatment., SCOPUS: ar.j, info:eu-repo/semantics/published

Highlights ► Microsatellite typing of Leishmania donovani complex isolates discriminates intercontinental groups. ► Genome-wide SNP profiling reveals diversity in a homogeneous population. ► Identification of a novel divergent lineage within a small geographic region. ► SNP-typing of samples resistant and sensitive to treatment drugs., The species of the Leishmania donovani species complex cause visceral leishmaniasis, a debilitating infectious disease transmitted by sandflies. Understanding molecular changes associated with population structure in these parasites can help unravel their epidemiology and spread in humans. In this study, we used a panel of standard microsatellite loci and genome-wide SNPs to investigate population-level diversity in L. donovani strains recently isolated from a small geographic area spanning India, Bihar and Nepal, and compared their variation to that found in diverse strains of the L. donovani complex isolates from Europe, Africa and Asia. Microsatellites and SNPs could clearly resolve the phylogenetic relationships of the strains between continents, and microsatellite phylogenies indicated that certain older Indian strains were closely related to African strains. In the context of the anti-malaria spraying campaigns in the 1960s, this was consistent with a pattern of episodic population size contractions and clonal expansions in these parasites that was supported by population history simulations. In sharp contrast to the low resolution provided by microsatellites, SNPs retained a much more fine-scale resolution of population-level variability to the extent that they identified four different lineages from the same region one of which was more closely related to African and European strains than to Indian or Nepalese ones. Joining results of in vitro testing the antimonial drug sensitivity with the phylogenetic signals from the SNP data highlighted protein-level mutations revealing a distinct drug-resistant group of Nepalese and Indian L. donovani. This study demonstrates the power of genomic data for exploring parasite population structure. Furthermore, markers defining different genetic groups have been discovered that could potentially be applied to investigate drug resistance in clinical Leishmania strains.

Natural hybridization events have been demonstrated between closely and distantly related Leishmania groups despite a predominantly clonal and endogamically sexual mode of reproduction. Here we report the first natural hybrid between Leishmania aethiopica and Leishmania donovani, as evidenced from the analysis of several clones from strain MHOM/ET/94/ABAUY. Targeted species-identification PCRs revealed the presence of both genotypes, and amplified fragment length polymorphisms indicated that the clones are genetically in an intermediate position between both parental species, being more closely related to L. aethiopica. The possible scenario facilitating hybrid formation is not clear, but is discussed in relation to epidemiological data. (C) 2011 Elsevier B.V. All rights reserved.

In order to understand the epidemiological dynamics of antimonial (SbV) resistance in zoonotic tegumentary leishmaniasis and its link with treatment outcome, we analyzed the population structure of 24 Peruvian Leishmania braziliensis clinical isolates with known in vitro antimony susceptibility and clinical phenotype by multilocus microsatellite typing (14 microsatellite loci). The genetic variability in the Peruvian isolates was high and the multilocus genotypes were strongly differentiated from each other. No correlation was found between the genotypes and in vitro drug susceptibility or clinical treatment outcome. The finding of a polyphyletic pattern among the SbV-resistant L. braziliensis might be explained by (i) independent events of drug resistance emergence, (ii) sexual recombination and/or (iii) other phenomena mimicking recombination signals. Interestingly, the polyphyletic pattern observed here is very similar to the one we observed in the anthroponotic Leishmania donovani (Laurent et al., 2007), hereby questioning the role of transmission and/or chemotherapeutic drug pressure in the observed population structure.

Antimonial (sodium stibogluconate, SSG) resistance and differentiation have been shown to be closely linked in Leishmania donovani, with SSG-resistant strains showing an increased capacity to generate infectious (metacyclic) forms. This is the first untargeted LC-MS metabolomics study which integrated both phenomena in one experimental design and provided insights into metabolic differences between three clinical L. donovani strains with a similar genetic background but different SSG-susceptibilities. We performed this analysis at different stages during promastigote growth and in the absence or presence of drug pressure. When comparing SSG-resistant and SSG-sensitive strains, a number of metabolic changes appeared to be constitutively present in all growth stages, pointing towards a clear link with SSG-resistance, whereas most metabolic changes were only detected in the stationary stage. These changes reflect the close intertwinement between SSG-resistance and an increased metacyclogenesis in resistant parasites. The metabolic changes suggest that SSG-resistant parasites have (i) an increased capacity for protection against oxidative stress; (ii) a higher fluidity of the plasma membrane; and (iii) a metabolic survival kit to better endure infection. These changes were even more pronounced in a resistant strain kept under Sb(III) drug pressure.

Antimonial (sodium stibogluconate, SSG) resistance and differentiation have been shown to be closely linked in Leishmania donovani, with SSG-resistant strains showing an increased capacity to generate infectious (metacyclic) forms. This is the first untargeted LC-MS metabolomics study which integrated both phenomena in one experimental design and provided insights into metabolic differences between three clinical L. donovani strains with a similar genetic background but different SSG-susceptibilities. We performed this analysis at different stages during promastigote growth and in the absence or presence of drug pressure. When comparing SSG-resistant and SSG-sensitive strains, a number of metabolic changes appeared to be constitutively present in all growth stages, pointing towards a clear link with SSG-resistance, whereas most metabolic changes were only detected in the stationary stage. These changes reflect the close intertwinement between SSG-resistance and an increased metacyclogenesis in resistant parasites. The metabolic changes suggest that SSG-resistant parasites have (i) an increased capacity for protection against oxidative stress; (ii) a higher fluidity of the plasma membrane; and (iii) a metabolic survival kit to better endure infection. These changes were even more pronounced in a resistant strain kept under Sb(III) drug pressure.

Visceral leishmaniasis is a potentially fatal disease endemic to large parts of Asia and Africa, primarily caused by the protozoan parasite Leishmania donovani. Here, we report a high-quality reference genome sequence for a strain of L. donovani from Nepal, and use this sequence to study variation in a set of 16 related clinical lines, isolated from visceral leishmaniasis patients from the same region, which also differ in their response to in vitro drug susceptibility. We show that whole-genome sequence data reveals genetic structure within these lines not shown by multilocus typing, and suggests that drug resistance has emerged multiple times in this closely related set of lines. Sequence comparisons with other Leishmania species and analysis of single-nucleotide diversity within our sample showed evidence of selection acting in a range of surface- and transport-related genes, including genes associated with drug resistance. Against a background of relative genetic homogeneity, we found extensive variation in chromosome copy number between our lines. Other forms of structural variation were significantly associated with drug resistance, notably including gene dosage and the copy number of an experimentally verified circular episome present in all lines and described here for the first time. This study provides a basis for more powerful molecular profiling of visceral leishmaniasis, providing additional power to track the drug resistance and epidemiology of an important human pathogen.

The heat-shock protein 70 gene (hsp70) has been exploited for Leishmania species identification in the New and Old World, using PCR followed by restriction fragment length polymorphism (RFLP) analysis. Current PCR presents limitations in terms of sensitivity, which hampers its use for analyzing clinical and biological samples, and specificity, which makes it inappropriate to discriminate between Leishmania and other trypanosomatids. The aim of the study was to improve the sensitivity and specificity of a previously reported hsp70 PCR using alternative PCR primers and RFLPs. Following in silico analysis of available sequences, three new PCR primer sets and restriction digest schemes were tested on a globally representative panel of 114 Leishmania strains, various other infectious agents, and clinical samples. The largest new PCR fragment retained the discriminatory power from RFLP, while two smaller fragments discriminated less species. The detection limit of the new PCRs was between 0.05 and 0.5 parasite genomes, they amplified clinical samples more efficiently, and were Leishmania specific. We succeeded in significantly improving the specificity and sensitivity of the PCRs for hsp70 Leishmania species typing. The improved PCR-RFLP assays can impact diagnosis, treatment, and epidemiological studies of leishmaniasis in any setting worldwide.

Diversity, phylogenetic, and population genetic studies of the genus Leishmania, causative agent of leishmaniasis, nowadays generally involve multilocus microsatellite and multilocus sequence typing. Even though these are well established and useful applications, amplified fragment length polymorphisms (AFLP) can provide complementary information. In addition, as the technique essentially probes the entire genome at random, without prior sequence knowledge, it is ideally suited as a screening tool for molecular markers linked with biological and clinical traits. We developed an AFLP protocol adapted to the Leishmania genome, tested its repeatability, and validated it on a panel of samples from the Leishmania donovani complex previously analyzed by multiple molecular tests. The technique proved highly reproducible, and showed that genetic relationships between L. donovani strains generally reflect geographic distance. Four main groups were identified: Leishmania infantum, African L. donovani, Indian L. donovani, and a mixed group consisting of L. donovani from India and Africa. Results were highly congruent with previous analyses on essentially the same sample set, indicating that the developed assay produces trustworthy data. This opens possibilities for application in studies of speciation and population dynamics. Moreover, it allows random screening of the entire Leishmania genome for linkage with biological and clinical parasite properties, such as fitness, drug resistance, and disease profile.

SUMMARYMathematical models predict that the future of epidemics of drug-resistant pathogens depends in part on the competitive fitness of drug-resistant strains. Considering metacyclogenesis (differentiation process essential for infectivity) as a major contributor to the fitness of Leishmania donovani, we tested its relationship with pentavalent antimony (SbV) resistance in clinical lines. Different methods for the assessment of metacyclogenesis were cross-validated: gene expression profiling (META1 and SHERP), morphometry (microscopy and FACS), in vitro infectivity to macrophages and resistance to complement lysis. This was done on a model constituted by 2 pairs of reference strains cloned from a SbV-resistant and -sensitive isolate. We selected the most adequate parameter and extended the analysis of metacyclogenesis diversity to a sample of 20 clinical lines with different in vitro susceptibility to the drug. The capacity of metacyclogenesis, as measured by the complement lysis test, was shown to be significantly higher in SbV-resistant clinical lines of L. donovani than in SbV-sensitive lines. Together with other lines of evidence, it is concluded that L. donovani constitutes a unique example and model of drug-resistant pathogens with traits of increased fitness. These findings raise a fundamental question about the potential risks of selecting more virulent pathogens through massive chemotherapeutic interventions.

For the epidemiological monitoring and clinical case management of leishmaniasis, determination of the causative Leishmania species gains importance. Current assays for the Old World often suffer from drawbacks in terms of validation on a geographically representative sample set and the ability to recognize all species complexes. We want to contribute to standardized species typing for Old World leishmaniasis. We determined the ribosomal DNA internal transcribed spacer 1 sequence of 24 strains or isolates, and validated four species-specific polymerase chain reactions (PCRs) amplifying this target. They discriminate L. aethiopica, L. tropica, L. major, and the L. donovani complex, use the same cycling conditions, and include an internal amplification control. Our PCRs amplify 0.1 pg of Leishmania DNA, while being 100% specific for species identification on an extensive panel of geographically representative strains and isolates. Similar results were obtained in an endemic reference laboratory in Kenya. Species could also be identified in clinical specimens. The presented PCRs require only agarose gel detection, and have several other advantages over many existing assays. We outline potential problems, suggest concrete solutions for transferring the technique to other settings, and deliver the proof-of-principle for analyzing clinical samples.

Background Leishmania donovani is an intracellular protozoan parasite that causes a lethal systemic disease, visceral leishmaniasis (VL), and is transmitted between mammalian hosts by phlebotomine sandflies. Leishmania expertly survives in these ‘hostile’ environments with a unique redox system protecting against oxidative damage, and host manipulation skills suppressing oxidative outbursts of the mammalian host. Treating patients imposes an additional stress on the parasite and sodium stibogluconate (SSG) was used for over 70 years in the Indian subcontinent. Methodology/Principal Findings We evaluated whether the survival capacity of clinical L. donovani isolates varies significantly at different stages of their life cycle by comparing proliferation, oxidative stress tolerance and infection capacity of 3 Nepalese L. donovani strains in several in vitro and in vivo models. In general, the two strains that were resistant to SSG, a stress encountered in patients, attained stationary phase at a higher parasite density, contained a higher amount of metacyclic parasites and had a greater capacity to cause in vivo infection in mice compared to the SSG-sensitive strain. Conclusions/Significance The 2 SSG-resistant strains had superior survival skills as promastigotes and as amastigotes compared to the SSG-sensitive strain. These results could indicate that Leishmania parasites adapting successfully to antimonial drug pressure acquire an overall increased fitness, which stands in contrast to what is found for other organisms, where drug resistance is usually linked to a fitness cost. Further validation experiments are under way to verify this hypothesis.