Ilse Maes

Whole genome sequencing of multiple Leishmania donovani clinical isolates provides insights into population structure and mechanisms of drug resistance

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.

Linking In Vitro and In Vivo Survival of Clinical Leishmania donovani Strains

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.

Three new sensitive and specific heat-shock protein 70 PCRs for global Leishmania species identification

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.

Evolutionary genomics of epidemic visceral leishmaniasis in the Indian subcontinent

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. DOI: http://dx.doi.org/10.7554/eLife.12613.001

Genome-wide SNP and microsatellite variation illuminate population-level epidemiology in the **Leishmania donovani** species complex

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.

Heat-shock protein 70 gene sequencing for Leishmania species typing in European tropical infectious disease clinics

We describe Leishmania species determination on clinical samples on the basis of partial sequencing of the heat-shock protein 70 gene (hsp70), without the need for parasite isolation. The method is especially suited for use in non-endemic infectious disease clinics dealing with relatively few cases on an annual basis, for which no fast high throughput diagnostic tests are needed. We show that the results obtained from this gene are in nearly perfect agreement with those from multilocus enzyme electrophoresis, which is still considered by many clinicians and the World Health Organization (WHO) as the gold standard in Leishmania species typing. Currently, 203 sequences are available that cover the entire hsp70 gene region analysed here, originating from a total of 41 leishmaniasis endemic countries, and representing 15 species and sub-species causing human disease. We also provide a detailed laboratory protocol that includes a step-by-step procedure of the typing methodology, to facilitate implementation in diagnostic laboratories.

Metabolic adaptations of Leishmania donovani in relation to differentiation, drug resistance, and 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 SbIII drug pressure.

Association of the Endobiont Double-Stranded RNA Virus LRV1 With Treatment Failure for Human Leishmaniasis Caused by Leishmania braziliensis in Peru and Bolivia

Cutaneous and mucosal leishmaniasis, caused in South America by Leishmania braziliensis, is difficult to cure by chemotherapy (primarily pentavalent antimonials [Sb(V)]). Treatment failure does not correlate well with resistance in vitro, and the factors responsible for treatment failure in patients are not well understood. Many isolates of L. braziliensis (>25%) contain a double-stranded RNA virus named Leishmaniavirus 1 (LRV1), which has also been reported in Leishmania guyanensis, for which an association with increased pathology, metastasis, and parasite replication was found in murine models. Here we probed the relationship of LRV1 to drug treatment success and disease in 97 L. braziliensis-infected patients from Peru and Bolivia. In vitro cultures were established, parasites were typed as L. braziliensis, and the presence of LRV1 was determined by reverse transcription-polymerase chain reaction, followed by sequence analysis. LRV1 was associated significantly with an increased risk of treatment failure (odds ratio, 3.99; P = .04). There was no significant association with intrinsic Sb(V) resistance among parasites, suggesting that treatment failure arises from LRV1-mediated effects on host metabolism and/or parasite survival. The association of LRV1 with clinical drug treatment failure could serve to guide more-effective treatment of tegumentary disease caused by L. braziliensis.

Increased metacyclogenesis of antimony-resistant Leishmania donovani clinical lines

Mathematical 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.

Evolution of the Leishmania braziliensis species complex from amplified fragment length polymorphisms, and clinical implications.

In order to get more insight into its evolution and geographical distribution, we investigated the Leishmania (Viannia) braziliensis species complex using amplified fragment length polymorphisms and sequencing of a heat-shock protein 70 gene fragment. Previously, several assays had alluded to the high genetic diversity of the group, and single-locus assays typically identified two species, i.e. L. braziliensis and Leishmania peruviana, with occasional genetic signatures of both in the same strain. By analysis of 53 parasite isolates from Peru, and eight additional ones from other countries, we identified an atypical L. braziliensis cluster, and confirmed the origin of L. peruviana from the L. braziliensis cluster during the colonization of the western Andean coastal valleys. We discuss the clinical and taxonomical implications of our findings in relation to currently used species typing assays.