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.

Evaluation of four single-locus markers for Leishmania species discrimination by sequencing.

Several genetic markers have been described for discriminating Leishmania species. In most reported cases, one or a few polymorphisms are the basis of species identification, and the methods were validated on a limited number of strains from a particular geographical region. Therefore, most techniques may underestimate the global intraspecies variability and are applicable only in certain areas. In addition, interlaboratory standardization is mostly absent, complicating comparisons among different studies. Here, we compared species typing results from all sequence polymorphisms found in four popular markers that can be applied directly on clinical samples: the miniexon or spliced leader, the internal transcribed spacer of the ribosomal DNA array, the 7SL RNA gene, and the heat shock protein 70 gene. Clustering was evaluated among 74 Leishmania strains, selected to represent a wide geographic distribution and genetic variability of the medically relevant species of the genus. Results were compared with a multilocus sequence typing (MLST) approach using 7 single-copy household genes and with multilocus enzyme electrophoresis (MLEE), still considered the gold standard by some. We show that strain groupings are highly congruent across the four different single-locus markers, MLST, and MLEE. Overall, the heat shock protein 70 gene and the miniexon presented the best resolutions for separating medically relevant species. As gene sequence analysis is validated here on a global scale, it is advocated as the method of choice for use in genetic, clinical, and epidemiological studies and for managing patients with unknown origins of infection, especially in Western infectious disease clinics dealing with imported leishmaniasis.

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.

Identification of Old World Leishmania species by PCR-RFLP of the 7 spliced leader RNA gene and reverse dot blot assay.

The aim of the study was to assess the 7SL RNA PCR followed by restriction fragment length polymorphism (RFLP) and reverse dot blot (RDB) assays for use in identification of Old World Leishmania species. Species-specific RFLP patterns were obtained for Leishmania major, Leishmania tropica and the Leishmania donovani complex when the 7SL RNA PCR product was digested with the restriction enzyme BsuRI, an isoschizomer of HaeIII. For the RDB assay, biotin-labelled 7SL RNA amplicons were hybridized to Leishmania genus-specific and species-specific oligonucleotide probes immobilized onto a membrane. The Old World Leishmania species could be distinguished by using five probes: one that was a genus-specific probe and hybridized to all Leishmania species (Lc), two that were specific for L. major (Lm1 and Lm2), one that was specific for L. tropica (Lt) and one that detected both L. major and L. tropica (Lmt). The PCR-RDB was 10 times more sensitive than 7SL PCR and can detect <1 parasite. In addition, the identification of species was easier and more reliable than with 7SL PCR-RFLP. 7SL PCR-RFLP detected parasites in 50 of 57 clinical samples, whereas PCR-RDB detected 53 and 55 were detected by amplification of kinetoplast (k) DNA. The 7SL RNA PCR has proven useful for direct diagnosis of Old World leishmaniasis, especially when combined with the RBD assay for species identification.

Specification of blood meals ingested by female sand flies caught in Palestinian foci and identification of their concomitant leishmanial infections.

Since leishmaniases are zoonotic vector-borne diseases transmitted through the bites of infected female sand flies, identification of the sources of imbibed blood meals and the detection and identification of leishmanial DNA in them are important in discerning animal reservoirs, clarifying the epidemiology and facilitating control of local leishmaniases. CDC light traps, aspirators and sticky paper traps were used to collect sand flies in four Palestinian foci of both, CL and VL. Phlebotomine species identification was based on morphological keys. Female specimens were screened to detect and identify leishmanial infections, using kDNA-PCR and ITS1-PCR, and engorged female specimens were analyzed to identify the origin of their blood meals, using an RDB blood meal assay based on the amplification of the cytochrome b gene (cytb) of vertebrate mitochondrial DNA (mtDNA). Twenty sand fly species, 11 of the genus Phlebotomus and nine the genus Sergentomyia, were identified. The most abundant species was Ph. papatasi (33.7%), followed by Ph. sergenti (21%). Among the 691 female sand fly specimens, 18.5% (128/691) were positive for leishmanial DNA, using the kDNA-PCR and 6.4% (44/691) were positive using the ITS1-PCR. DNA from parasites of the genus Leishmania was identified in only 1.5% of the infected sand flies. That of Leishmania tropica parasites was detected in six female specimens of Ph. sergenti and that of L. major parasites in two female specimens of Ph. papatasi. Interestingly, two engorged females of the species Se. (Neophlebotomus) sp. were positive for L. tropica DNA. Ninety engorged female sand flies of Ph. papatasi and 104 of Ph. sergenti had fed on a large variety of vertebrate hosts such as humans, hyraxes, rats, cows, goats and birds. Regarding blood-meals showing a mixture from different species of animal host, hyrax and rat blood was revealed in 8/104 (7.7%) females of Ph. sergenti. Detection of hyrax blood in engorged female sand flies of the species Ph. sergenti supports the role of hyraxes being a potential reservoir of L. tropica in Palestinian regions. Rat blood meals might be significant since a few strains L. tropica and L. infantum were isolated from rats. Further studies must be undertaken before conclusions could be drawn.

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.

Development of a multilocus microsatellite typing approach for discriminating strains of Leishmania (Viannia) species.

A multilocus microsatellite typing (MLMT) approach based on the analysis of 15 independent loci has been developed for the discrimination of strains belonging to different Viannia species. Thirteen microsatellite loci were isolated de novo from microsatellite-enriched libraries for both Leishmania braziliensis and L. guyanensis. Two previously identified markers, AC01 and AC16, were modified and added to our marker set. Markers were designed to contain simple dinucleotide repeats flanked by the minimal possible number of nucleotides in order to allow variations in repeat numbers to be scored as size variations of the PCR products. The 15 markers in total were amplified for almost all of the strains of Viannia tested; one marker did not amplify from the two L. peruviana strains included in the study. When 30 strains of L. braziliensis, 21 strains of L. guyanensis, and 2 strains of L. peruviana were tested for polymorphisms, all strains except two strains of L. guyanensis had individual MLMT types. Distance-based analysis identified three main clusters. All strains except one strain of L. guyanensis grouped together. Two clusters consisted of strains of L. braziliensis according to their geographical origins. The two strains of L. peruviana grouped together with strains of L. braziliensis from Peru and the adjacent Brazilian state of Acre. MLMT has proven capable of individualizing strains even from the same areas of endemicity and of detecting genetic structures at different levels. MLMT is thus applicable for epidemiological and population genetic studies of strains within the subgenus Viannia.

Multilocus microsatellite typing (MLMT) reveals genetic homogeneity of Leishmania donovani strains in the Indian subcontinent.

In this population genetic study of Leishmania donovani parasites in the Indian subcontinent, 132 isolates obtained from patients in Bangladesh, India, Nepal and Sri Lanka suffering from Kala-azar (100), post-Kala-azar dermal leishmaniasis (PKDL) (25) and cutaneous leishmaniasis (CL) (2), and from 5 patients whose clinical patterns were not defined, were analysed by using 15 hyper-variable microsatellite loci. Multilocus microsatellite typing (MLMT) data were analysed by using a Bayesian model-based clustering algorithm and constructing phylogenic tree based on genetic distances. In total, 125 strains from Bangladesh, Bihar (India) and Nepal formed a very homogeneous population regardless of geographical origin, clinical manifestation, and whether they presented in vitro or in vivo susceptibility to antimonial drugs. Identical multilocus microsatellite profiles were found for 108 strains, other strains differed in only one marker. Considerably different microsatellite profiles were identified for three Indian strains most closely related to L. donovani from Kenya, and for four strains from Indian and Sri Lankan CL cases. The circulation of a single homogeneous population of L. donovani in Bihar (India), Bangladesh and Nepal is, most probably, related to the epidemic spread of visceral leishmaniasis in this area.

PCR diagnosis of visceral leishmaniasis in an endemic region, Mymensingh district, Bangladesh.

Detection of Leishmania parasites in a clinical sample is necessary to confirm a suspected case of leishmaniasis. We compared the sensitivity of internal transcribed spacer 1-PCR (ITS 1-PCR) assay for parasite diagnosis with that of microscopic detection in clinical samples from kala-azar (KA) or post-kala-azar dermal leishmaniasis (PKDL) suspects in Mymensingh. Of 39 specimens collected from 35 KA and four PKDL suspects, 26 were positive by microscopic examination of smears from bone marrow and skin exudates; 38 specimens spotted on filter paper and 27 of the 28 Giemsa-stained slides tested by PCR proved positive by ITS1-PCR.

Kinetoplast DNA heterogeneity among Leishmania infantum strains in central Israel and Palestine.

Leishmania infantum is the causative agent of canine and human visceral leishmaniasis in Israel and Palestine. Amplification of the kinetoplast DNA (kDNA), followed by restriction fragment length polymorphism (RFLP) analysis with restriction enzymes was used to examine the genotypic association between L. infantum strains isolated from 22 dogs and 2 humans from these adjoining regions. Results showed wide kDNA heterogeneity in these strains. Two main clusters (A and B) were identified. Cluster A was restricted to central Israel and was mainly found in strains isolated after 2002 whereas cluster B included parasites from central Israel and the West Bank. The kDNA microheterogeneity in L. infantum parasite populations as shown by genotyping with the kDNA-PCR and RFLP provided a tool to study the epidemiology of the disease and track its spread in central Israel and Palestine.

A novel multilocus sequence typing scheme identifying genetic diversity amongst Leishmania donovani isolates from a genetically homogeneous population in the Indian subcontinent.

In the Indian subcontinent, infection with Leishmania donovani can cause fatal visceral leishmaniasis. Genetic variation in L. donovani is believed to occur rapidly from environmental changes and through selective drug pressures, thereby allowing continued disease occurrence in this region. All previous molecular markers that are commonly in use multilocus microsatellite typing and multilocus sequence typing, were monomorphic in L. donovani originating from the Indian subcontinent (with only a few exceptions) and hence are not suitable for this region. An multilocus sequence typing scheme consisting of a new set of seven housekeeping genes was developed in this study, based on recent findings from whole genome sequencing data. This new scheme was used to assess the genetic diversity amongst 22 autochthonous L. donovani isolates from Bangladesh. Nineteen additional isolates of the L. donovani complex (including sequences of L. donovani reference strain BPK282A1) from other countries were included for comparison. By using restriction fragment length polymorphism of the internal transcribed spacer 1 region (ITS1-RFLP) and ITS1 sequencing, all Bangladeshi isolates were confirmed to be L. donovani. Population genetic analyses of 41 isolates using the seven new MLST loci clearly separated L. donovani from Leishmania infantum. With this multilocus sequence typing scheme, seven genotypes were identified amongst Bangladeshi L. donovani isolates, and these isolates were found to be phylogenetically different compared with those from India, Nepal, Iraq and Africa. This novel multilocus sequence typing approach can detect intra- and inter-species variations within the L. donovani complex, but most importantly these molecular markers can be applied to resolve the phylogenetically very homogeneous L. donovani strains from the Indian subcontinent. Four of these markers were found suitable to differentiate strains originating from Bangladesh, with marker A2P being the most discriminative one.

Identification of geographically distributed sub-populations of Leishmania (Leishmania) major by microsatellite analysis.

BACKGROUND: Leishmania (Leishmania) major, one of the agents causing cutaneous leishmaniasis (CL) in humans, is widely distributed in the Old World where different species of wild rodent and phlebotomine sand fly serve as animal reservoir hosts and vectors, respectively. Despite this, strains of L. (L.) major isolated from many different sources over many years have proved to be relatively uniform. To investigate the population structure of the species highly polymorphic microsatellite markers were employed for greater discrimination among it's otherwise closely related strains, an approach applied successfully to other species of Leishmania. RESULTS: Multilocus Microsatellite Typing (MLMT) based on 10 different microsatellite markers was applied to 106 strains of L. (L.) major from different regions where it is endemic. On applying a Bayesian model-based approach, three main populations were identified, corresponding to three separate geographical regions: Central Asia (CA); the Middle East (ME); and Africa (AF). This was congruent with phylogenetic reconstructions based on genetic distances. Re-analysis separated each of the populations into two sub-populations. The two African sub-populations did not correlate well with strains' geographical origin. Strains falling into the sub-populations CA and ME did mostly group according to their place of isolation although some anomalies were seen, probably, owing to human migration. CONCLUSION: The model- and distance-based analyses of the microsatellite data exposed three main populations of L. (L.) major, Central Asia, the Middle East and Africa, each of which separated into two sub-populations. This probably correlates with the different species of rodent host.

Leishmaniases in the Mediterranean in the era of molecular epidemiology.

Molecular tools are used increasingly for descriptive epidemiological studies in different Mediterranean foci of visceral and cutaneous leishmaniases. Several molecular markers with different resolution levels have been developed to address key epidemiological questions related to the (re-)emergence and spread of leishmaniases, as well as its risk factors: environmental changes, immunosuppression and treatment failure. Typing and analytical tools are improving but are not yet addressing all epidemiological issues satisfactorily. There is an urgent need for better cooperation between laboratory scientists and epidemiologists and for regional epidemiological surveillance of these infectious diseases that affect all Mediterranean countries.

Genetic polymorphism of Algerian Leishmania infantum strains revealed by multilocus microsatellite analysis.

The present study applies multilocus microsatellite typing (MLMT) for studying the polymorphism among 55 strains of Leishmania infantum from Algeria. These strains from different Algerian foci representing different zymodemes, hosts and clinical forms were analysed using 14 microsatellite markers. All 55 strains had individual MLMT profiles and no relationship was observed between them and different host or geographical origins. Three populations of Algerian L. infantum were identified by a Bayesian clustering approach implemented in STRUCTURE software and supported by genetic distance analysis. Two populations, A and B, consisted mainly of strains belonging to zymodeme MON-1, and the third population, C, mainly of MON-24 strains isolated from cutaneous leishmaniasis cases. Interestingly, a small group of strains appeared as a mixture of different populations and might be putative hybrids. Genetic migration was noticed among the two MON-1 populations, A and B, as well as between populations A and C. Due to its high discriminatory power MLMT could be also successfully applied for differentiating relapses or re-infection for patients suffering from multiple episodes of visceral leishmaniasis.

Molecular diagnosis of Old World cutaneous leishmaniasis and species identification by use of a reverse line blot hybridization assay.

Reverse line blot hybridization assays (RLB) have been used for the rapid diagnosis and genotyping of many pathogens. The leishmaniases are caused by a large number of species, and rapid, accurate parasite characterization is important in deciding on appropriate therapy. Fourteen oligonucleotide probes, 2 genus specific and 12 species specific (2 specific for Leishmania major, 3 for L. tropica, 1 for L. infantum, 3 for L. donovani, and 3 for L. aethiopica), were prepared by using DNA sequences in the internal transcribed spacer 1 (ITS1) region of the rRNA genes. Probe specificity was evaluated by amplifying DNA from 21 reference strains using biotinylated ITS1 PCR primers and the RLB. The genus-specific probes, PP and PP3', recognized all Leishmania species examined, while the species-specific probes were able to distinguish between all the Old World Leishmania species. Titrations using purified parasite DNA showed that the RLB is 10- to 100-fold more sensitive than ITS1 PCR and can detect <0.1 pg DNA. The RLB was compared to kinetoplast DNA (kDNA) and ITS1 PCR by using 67 samples from suspected cutaneous leishmaniasis (CL) patients in Israel and the West Bank. The RLB accurately identified 58/59 confirmed positive samples as CL, a result similar to that found by kDNA PCR (59/59) and better than that by ITS1 PCR (50/59). The positive predictive value and negative predictive value of the RLB were 95.1% and 83.3%, respectively. L. major or L. tropica was identified by the RLB in 55 of the confirmed positive cases, a level of accuracy better than that of ITS1 PCR with restriction fragment length polymorphism (42/59). Thus, RLB can be used to diagnose and characterize Old World CL.

Leishmania naiffi and Leishmania guyanensis reference genomes highlight genome structure and gene evolution in the Viannia subgenus.

The unicellular protozoan parasite Leishmania causes the neglected tropical disease leishmaniasis, affecting 12 million people in 98 countries. In South America, where the Viannia subgenus predominates, so far only L. (Viannia) braziliensis and L. (V.) panamensis have been sequenced, assembled and annotated as reference genomes. Addressing this deficit in molecular information can inform species typing, epidemiological monitoring and clinical treatment. Here, L. (V.) naiffi and L. (V.) guyanensis genomic DNA was sequenced to assemble these two genomes as draft references from short sequence reads. The methods used were tested using short sequence reads for L. braziliensis M2904 against its published reference as a comparison. This assembly and annotation pipeline identified 70 additional genes not annotated on the original M2904 reference. Phylogenetic and evolutionary comparisons of L. guyanensis and L. naiffi with 10 other Viannia genomes revealed four traits common to all Viannia: aneuploidy, 22 orthologous groups of genes absent in other Leishmania subgenera, elevated TATE transposon copies and a high NADH-dependent fumarate reductase gene copy number. Within the Viannia, there were limited structural changes in genome architecture specific to individual species: a 45 Kb amplification on chromosome 34 was present in all bar L. lainsoni, L. naiffi had a higher copy number of the virulence factor leishmanolysin, and laboratory isolate L. shawi M8408 had a possible minichromosome derived from the 3' end of chromosome 34. This combination of genome assembly, phylogenetics and comparative analysis across an extended panel of diverse Viannia has uncovered new insights into the origin and evolution of this subgenus and can help improve diagnostics for leishmaniasis surveillance.

Multilocus microsatellite typing (MLMT) reveals genetically isolated populations between and within the main endemic regions of visceral leishmaniasis.

Multilocus enzyme electrophoresis (MLEE) is the gold standard for taxonomy and strain typing of Leishmania, but has some limitations. An alternative reliable and fast genotyping method for addressing population genetic and key epidemiological questions, is multilocus microsatellite typing (MLMT). MLMT using 15 markers was applied to 91 strains of L. donovani, L. archibaldi, L. infantum and L. chagasi from major endemic regions of visceral leishmaniasis. Population structures were inferred by combination of Bayesian model-based and distance-based approaches. Six main genetically distinct populations were identified: (1) L. infantum/L. chagasi MON-1 and (2) L. infantum/L. chagasi non-MON-1, both Mediterranean region/South America; (3) L. donovani (MON-18), L. archibaldi (MON-82), L. infantum (MON-30, 81) and (4) L. donovani (MON-31, 274), L. archibaldi (MON-82, 257, 258), L. infantum (MON-267), both Sudan/Ethiopia; (5) L. donovani MON-2, India; (6) L. donovani (MON-36, 37, 38), Kenya and India. Substructures according to place and time of strain isolation were detected. The VL populations seem to be predominantly clonal with a high level of inbreeding. Allelic diversity was highest in the Mediterranean region, intermediate in Africa and lowest in India. MLMT provides a powerful tool for global taxonomic, population genetic and epidemiological studies of the L.donovani complex.

Epidemiological dynamics of antimonial resistance in Leishmania donovani: genotyping reveals a polyclonal population structure among naturally-resistant clinical isolates from Nepal.

Pentavalent antimonials (SbV) are the first line drug against leishmaniasis worldwide, but drug resistance is increasingly reported, particularly in the Indian sub-continent, where it represents a major threat for the control of anthroponotic visceral leishmaniasis (VL). In order to understand the epidemiological dynamics of antimonial resistance in anthroponotic VL, we analysed here the population structure of 24 Leishmania donovani stocks isolated from anthroponotic VL-patients from Eastern Nepal: 13 SbV-naturally resistant and 11 SbV-sensitive, as demonstrated by in vitro drug susceptibility assays. The parasites were genotyped by PCR-RFLP analysis of kDNA minicircles and by microsatellite analysis and the encountered polymorphism revealed a polyclonal structure among resistant isolates. Furthermore, analysis of paired samples obtained from the same patients before treatment and after failure revealed primary as well as acquired resistance. The hypothesis of independent events of drug resistance emergence is proposed and confronted to alternative explanations. Our results show the dynamics of drug resistance epidemiology and highlight the importance of surveillance networks.