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.

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.

[PCR-RFLP and RAPD for typing neotropical Leishmania]

INTRODUCTION

The analysis of the PCR-restriction fragment length polymorphism and random amplified polymorphic DNA have been useful tools for Leishmania identification.

OBJECTIVES

Molecular procedures were demonstrated for identification and typing of reference strains of New World Leishmania and their applicability was validated for clinical samples.

MATERIALS AND METHODS

DNA was extracted from 16 reference strains of Latin American Leishmania as well as from clinical samples of leishmaniasis patients. A sequence coding for cysteine proteinase B was amplified by PCR and subjected to restriction fragment length polymorphism analysis. The enzyme used was Taq1. For eight of the reference strains, the random amplified polymorphic desoxyribonucleic acid technique (RAPD) was applied. Band patterns for Leishmania species differentiation were established each each method. The sample size of the clinical sample was of 5.

RESULTS

PCR products of the cysteine proteinase B gene were obtained for L. braziliensis, L. peruviana, L. panamensis and L. guyanensis. For the other species, L. mexicana, L. amazonensis, L. garnhami, L. lainsoni, L. chagasi, L. naiffi, no amplification occurred. The patterns of restriction fragments revealed band patterns in common for L. peruviana, L. guyanensis and L. panamensis, whereas L. braziliensis had a distinctive pattern. When human samples were examined, amplification occurred for all cases, and the profiles corresponded to the common profile of L. peruviana, L. guyanensis and L. panamensis. The RAPD technique demonstrated reproducible and distinctive patterns for each of the 8 reference strains, L. mexicana, L. amazonensis, L. garnhami, L. lainsoni, L. chagasi, L. naiffi, making possible to differentiate all them. The advantages and limitations of each procedure are discussed.

CONCLUSIONS

The combination of RFP and RAPD methodologies provide useful tools to identify medical important species of Leishmania by recognizing DNA sequences characteristic of each species.