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.

Genetic markers for SSG resistance in Leishmania donovani and SSG treatment failure in visceral leishmaniasis patients of the Indian subcontinent.

The current standard to assess pentavalent antimonial (SSG) susceptibility of Leishmania is a laborious in vitro assay of which the result has little clinical value because SSG-resistant parasites are also found in SSG-cured patients. Candidate genetic markers for clinically relevant SSG-resistant parasites identified by full genome sequencing were here validated on a larger set of clinical strains. We show that 3 genomic locations suffice to specifically detect the SSG-resistant parasites found only in patients experiencing SSG treatment failure. This finding allows the development of rapid assays to monitor the emergence and spread of clinically relevant SSG-resistant Leishmania parasites.

Global genome diversity of the Leishmania donovani complex.

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.

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.

Drug susceptibility in Leishmania isolates following miltefosine treatment in cases of visceral leishmaniasis and post kala-azar dermal leishmaniasis.

BACKGROUND: With widespread resistance to antimonials in Visceral Leishmaniasis (VL) in the Indian subcontinent, Miltefosine (MIL) has been introduced as the first line therapy. Surveillance of MIL susceptibility in natural populations of Leishmania donovani is vital to preserve it and support the VL elimination program. METHODOLOGY AND PRINCIPAL FINDINGS: We measured in vitro susceptibility towards MIL and paromomycin (PMM) in L. donovani isolated from VL and PKDL, pre- and post-treatment cases, using an amastigote-macrophage model. MIL susceptibility of post-treatment isolates from cured VL cases (n = 13, mean IC(50)±SD = 2.43±1.44 µM), was comparable (p>0.05) whereas that from relapses (n = 3, mean IC(50) = 4.72±1.99 µM) was significantly higher (p = 0.04) to that of the pre-treatment group (n = 6, mean IC(50) = 1.86±0.75 µM). In PKDL, post-treatment isolates (n = 3, mean IC(50) = 16.13±2.64 µM) exhibited significantly lower susceptibility (p = 0.03) than pre-treatment isolates (n = 5, mean IC(50) = 8.63±0.94 µM). Overall, PKDL isolates (n = 8, mean IC(50) = 11.45±4.19 µM) exhibited significantly higher tolerance (p<0.0001) to MIL than VL isolates (n = 22, mean IC(50) = 2.58±1.58 µM). Point mutations in the miltefosine transporter (LdMT) and its beta subunit (LdRos3) genes previously reported in parasites with experimentally induced MIL resistance were not present in the clinical isolates. Further, the mRNA expression profile of these genes was comparable in the pre- and post-treatment isolates. Parasite isolates from VL and PKDL cases were uniformly susceptible to PMM with respective mean IC(50) = 7.05±2.24 µM and 6.18±1.51 µM. CONCLUSION: The in vitro susceptibility of VL isolates remained unchanged at the end of MIL treatment; however, isolates from relapsed VL and PKDL cases had lower susceptibility than the pre-treatment isolates. PKDL isolates were more tolerant towards MIL in comparison with VL isolates. All parasite isolates were uniformly susceptible to PMM. Mutations in the LdMT and LdRos3 genes as well as changes in the expression of these genes previously correlated with experimental resistance to MIL could not be verified for the field isolates.

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

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.