Single locus genotyping to track Leishmania donovani in the Indian subcontinent: Application in Nepal

Source Tracing of Leishmania donovani in Emerging Foci of Visceral Leishmaniasis, Western Nepal

We sequenced Leishmania donovani genomes in blood samples collected in emerging foci of visceral leishmaniasis in western Nepal. We detected lineages very different from the preelimination main parasite population, including a new lineage and a rare one previously reported in eastern Nepal. Our findings underscore the need for genomic surveillance.

Leishmania LPG interacts with LRR5/LRR6 of macrophage TLR4 for parasite invasion and impairs the macrophage functions

Visceral leishmaniasis (VL) is a severe form of leishmaniasis, primarily affecting the poor in developing countries. Although several studies have highlighted the importance of toll-like receptors (TLRs) in the pathophysiology of leishmaniasis, the role of specific TLRs and their binding partners involved in Leishmania donovani uptake are still elusive. To investigate the mechanism of L. donovani entry inside the macrophages, we found that the parasite lipophosphoglycan (LPG) interacted with the macrophage TLR4, leading to parasite uptake without any significant alteration of macrophage cell viability. Increased parasite numbers within macrophages markedly inhibited lipopolysachharide-induced pro-inflammatory cytokines gene expression. Silencing of macrophage-TLR4, or inhibition of parasite-LPG, significantly stemmed parasite infection in macrophages. Interestingly, we observed a significant enhancement of macrophage migration, and generation of reactive oxygen species (ROS) in the parasite-infected TLR4-silenced macrophages, whereas parasite infection in TLR4-overexpressed macrophages exhibited a notable reduction of macrophage migration and ROS generation. Moreover, mutations in the leucine-rich repeats (LRRs), particularly LRR5 and LRR6, significantly prevented TLR4 interaction with LPG, thus inhibiting cellular parasite entry. All these results suggest that parasite LPG recognition by the LRR5 and LRR6 of macrophage-TLR4 facilitated parasite entry, and impaired macrophage functions. Therefore, targeting LRR5/LRR6 interactions with LPG could provide a novel option to prevent VL.

The Leishmania donovani species complex: A new insight into taxonomy☆

Among the 20 or so Leishmania spp. described as pathogenic for humans, those of the Leishmania donovani complex are the exclusive causative agents of systemic and fatal visceral leishmaniasis. Although well studied, the complex is taxonomically controversial, which hampers clinical and epidemiological research. In this work, we analysed 56 Leishmania strains previously identified as L. donovani, Leishmania archibaldi or Leishmania infantum, isolated from humans, dogs and sandfly vectors throughout their distribution area. The strains were submitted to biochemical and genetic analyses and the resulting data were compared for congruence. Our results show: i) a partial concordance between biochemical and genetic-based data, ii) very limited genetic variability within the L. donovani complex, iii) footprints of frequent genetic exchange along an east-west gradient, marked by a widespread diffusion of alleles across the geographical range, and iv) a large-scale geographical spreading of a few genotypes. From a taxonomic point of view, considering the absence of relevant terminology in existing classes, the L. donovani complex could be treated as a single entity.

Sand Fly Studies Predict Transmission Potential of Drug-resistant Leishmania

Leishmania parasites have the capacity to rapidly adapt to changing environments in their digenetic life cycle which alternates between a vertebrate and an invertebrate host. Emergence of resistance following drug exposure can evoke phenotypic alterations that affect several aspects of parasite fitness in both hosts. Current studies of the impact of resistance are mostly limited to interactions with the mammalian host and characterization of in vitro parasite growth and differentiation. Development in the vector and transmission capacity have been largely ignored. This review reflects on the impact of drug resistance on its spreading potential with specific focus on the use of the sand fly infection model to evaluate parasite development in the vector and the ensuing transmission potential of drug-resistant phenotypes.

Next-Generation Molecular Surveillance of TriTryp Diseases

Elimination programs targeting TriTryp diseases (Leishmaniasis, Chagas' disease, human African trypanosomiasis) significantly reduced the number of cases. Continued surveillance is crucial to sustain this progress, but parasite molecular surveillance by genotyping is currently lacking. We explain here which epidemiological questions of public health and clinical relevance could be answered by means of molecular surveillance. Whole-genome sequencing (WGS) for molecular surveillance will be an important added value, where we advocate that preference should be given to direct sequencing of the parasite's genome in host tissues instead of analysis of cultivated isolates. The main challenges here, and recent technological advances, are discussed. We conclude with a series of recommendations for implementing whole-genome sequencing for molecular surveillance.

Genomes of Leishmania parasites directly sequenced from patients with visceral leishmaniasis in the Indian subcontinent

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

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.

ISC1, a new Leishmania donovani population emerging in the Indian sub-continent: Vector competence of Phlebotomus argentipes

Visceral leishmaniasis (VL), the most severe form of the disease, is caused by Leishmania donovani in the Indian sub-continent (ISC). Whole genome sequencing studies revealed that two parasite populations exist in the ISC: a main population named the Core Group (CG) found mostly in the lowlands, and a new, genetically different subpopulation called ISC1. Parasites belonging to the CG were shown to be responsible for the recent epidemics, while the ISC1 variant was originally identified in hilly districts of Nepal and was later on increasingly found in the lowlands. Importantly, the ISC1 and CG isolates differ in their drug susceptibility and virulence signatures, suggesting that ISC1 constitutes an emerging and functionally different variant of L. donovani. In present study we aimed to address the potential of ISC1 transmission by the natural vector of L. donovani in the lowlands, Phlebotomus argentipes. By experimental infection of sand flies with parasites of the different genotypes, we demonstrate that ISC1 and CG strains are developing similarly in P. argentipes, suggesting that P. argentipes is a fully competent vector for ISC1 parasites. Integration of previous and current findings shows thus that ISC1 is a new and different variant of L. donovani, fully adapted to spread in the ISC through the main vector. This information is directly useful for managers of the elimination program. Furthermore, integration of our successive studies (genotyping, phenotyping and vector competence) demonstrates the relevance of molecular surveillance and should be of interest for scientists working on vector borne diseases and control managers.

Application of kDNA Minicircle PCR-RFLP to CharacterizeLeishmania donovaniClinical Isolates Obtained from Post-Kala-Azar Dermal Leishmaniasis in Eastern Nepal

Post-kala-azar dermal leishmaniasis (PKDL) is a skin manifestation of visceral leishmaniasis (VL) which develops after apparent cure in some patients. PKDL is considered as the potential reservoir for the VL infection. Molecular epidemiological characterization of L. donovani isolates obtained from VL and PKDL isolates is essentially required in order to understand the transmission dynamics of the VL infection. To date, genetic variation among the VL and PKDL L. donovani isolates was not fully elucidated. Therefore, 14 clinical isolates from VL and 4 clinical isolates from PKDL were speciated by hsp70 and rDNA genes. Further characterization of L. donovani by haspB PCR demonstrates two different genotypes. All PKDL isolates have the same genetic structure. kDNA PCR-RFLP assay revealed 18 different genotypes; however, structural analysis showed the two distinct kDNA genotype population (k = 2). The kDNA fingerprint patterns of parasites from hilly districts were clustered separately from low-land districts. Therefore, further study with a large number of samples is urgently required for systematic characterization of the clinical isolates to track the molecular epidemiology of the Leishmania donovani causing VL and the role of PKDL as a reservoir.

Genomic insights into virulence mechanisms of Leishmania donovani: evidence from an atypical strain

Background

Leishmaniasis is a neglected tropical disease with diverse clinical phenotypes, determined by parasite, host and vector interactions. Despite the advances in molecular biology and the availability of more Leishmania genome references in recent years, the association between parasite species and distinct clinical phenotypes remains poorly understood. We present a genomic comparison of an atypical variant of Leishmania donovani from a South Asian focus, where it mostly causes cutaneous form of leishmaniasis.

Results

Clinical isolates from six cutaneous leishmaniasis patients (CL-SL); 2 of whom were poor responders to antimony (CL-PR), and two visceral leishmaniasis patients (VL-SL) were sequenced on an Illumina MiSeq platform. Chromosome aneuploidy was observed in both groups but was more frequent in CL-SL. 248 genes differed by 2 fold or more in copy number among the two groups. Genes involved in amino acid use (LdBPK_271940) and energy metabolism (LdBPK_271950), predominated the VL-SL group with the same distribution pattern reflected in gene tandem arrays. Genes encoding amastins were present in higher copy numbers in VL-SL and CL-PR as well as being among predicted pseudogenes in CL-SL. Both chromosome and SNP profiles showed CL-SL and VL-SL to form two distinct groups. While expected heterozygosity was much higher in VL-SL, SNP allele frequency patterns did not suggest potential recent recombination breakpoints. The SNP/indel profile obtained using the more recently generated PacBio sequence did not vary markedly from that based on the standard LdBPK282A1 reference. Several genes previously associated with resistance to antimonials were observed in higher copy numbers in the analysis of CL-PR. H-locus amplification was seen in one cutaneous isolate which however did not belong to the CL-PR group.

Conclusions

The data presented suggests that intra species variations at chromosome and gene level are more likely to influence differences in tropism as well as response to treatment, and contributes to greater understanding of parasite molecular mechanisms underpinning these differences. These findings should be substantiated with a larger sample number and expression/functional studies.

Electronic supplementary material

The online version of this article (10.1186/s12864-018-5271-z) contains supplementary material, which is available to authorized users.

The role of case proximity in transmission of visceral leishmaniasis in a highly endemic village in Bangladesh

BACKGROUND

Visceral leishmaniasis (VL) is characterised by a high degree of spatial clustering at all scales, and this feature remains even with successful control measures. VL is targeted for elimination as a public health problem in the Indian subcontinent by 2020, and incidence has been falling rapidly since 2011. Current control is based on early diagnosis and treatment of clinical cases, and blanket indoor residual spraying of insecticide (IRS) in endemic villages to kill the sandfly vectors. Spatially targeting active case detection and/or IRS to higher risk areas would greatly reduce costs of control, but its effectiveness as a control strategy is unknown. The effectiveness depends on two key unknowns: how quickly transmission risk decreases with distance from a VL case and how much asymptomatically infected individuals contribute to transmission.

METHODOLOGY/PRINCIPAL FINDINGS

To estimate these key parameters, a spatiotemporal transmission model for VL was developed and fitted to geo-located epidemiological data on 2494 individuals from a highly endemic village in Mymensingh, Bangladesh. A Bayesian inference framework that could account for the unknown infection times of the VL cases, and missing symptom onset and recovery times, was developed to perform the parameter estimation. The parameter estimates obtained suggest that, in a highly endemic setting, VL risk decreases relatively quickly with distance from a case-halving within 90m-and that VL cases contribute significantly more to transmission than asymptomatic individuals.

CONCLUSIONS/SIGNIFICANCE

These results suggest that spatially-targeted interventions may be effective for limiting transmission. However, the extent to which spatial transmission patterns and the asymptomatic contribution vary with VL endemicity and over time is uncertain. In any event, interventions would need to be performed promptly and in a large radius (≥300m) around a new case to reduce transmission risk.

Integrated genomic and metabolomic profiling of ISC1, an emerging Leishmania donovani population in the Indian subcontinent

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.

Drug resistance and treatment failure in leishmaniasis: A 21st century challenge

Reevaluation of treatment guidelines for Old and New World leishmaniasis is urgently needed on a global basis because treatment failure is an increasing problem. Drug resistance is a fundamental determinant of treatment failure, although other factors also contribute to this phenomenon, including the global HIV/AIDS epidemic with its accompanying impact on the immune system. Pentavalent antimonials have been used successfully worldwide for the treatment of leishmaniasis since the first half of the 20th century, but the last 10 to 20 years have witnessed an increase in clinical resistance, e.g., in North Bihar in India. In this review, we discuss the meaning of “resistance” related to leishmaniasis and discuss its molecular epidemiology, particularly for Leishmania donovani that causes visceral leishmaniasis. We also discuss how resistance can affect drug combination therapies. Molecular mechanisms known to contribute to resistance to antimonials, amphotericin B, and miltefosine are also outlined.

Chemotherapy is central to the control and management of leishmaniasis. Antimonials remain the primary drugs against different forms of leishmaniasis in several regions. However, resistance to antimony has necessitated the use of alternative medications, especially in the Indian subcontinent (ISC). Compounds, notably the orally available miltefosine (MIL), parenteral paromomycin, and amphotericin B (AmB), are increasingly used to treat leishmaniasis. Although treatment failure (TF) has been observed in patients treated with most anti-leishmanials, its frequency of appearance may be important in patients treated with MIL, which has replaced antimonials within the kala-azar elimination program in the ISC. AmB is highly efficacious, and the associated toxic effects—when administered in its free deoxycholate form—are somewhat ameliorated in its liposomal formulation. Regrettably, laboratory experimentation has demonstrated a risk of resistance towards AmB as well. The rise of drug resistance impacts treatment outcome, and understanding its causes, spread, and impact will help us manage the risks it imposes. Here, we review the problem of TF in leishmaniasis and the contribution of drug resistance to the problem. Molecular mechanisms causing resistance to anti-leishmanials are discussed along with the appropriate use of additional available drugs, as well as the urgent need to consolidate strategies to monitor drug efficacy, epidemiological surveillance, and local policies. Coordination of these activities in national and international programs against leishmaniasis might represent a successful guide to further research and prevention activities.

Comparative analyses of whole genome sequences of Leishmania infantum isolates from humans and dogs in northeastern Brazil

The genomic sequences of 20 Leishmania infantum isolates collected in northeastern Brazil were compared with each other and with the available genomic sequences of 29 L. infantum/donovani isolates from Nepal and Turkey. The Brazilian isolates were obtained in the early 1990s or since 2009 from patients with visceral or non-ulcerating cutaneous leishmaniasis, asymptomatic humans, or dogs with visceral leishmaniasis. Two isolates were from the blood and bone marrow of the same visceral leishmaniasis patient. All 20 genomic sequences display 99.95% identity with each other and slightly less identity with a reference L. infantum genome from a Spanish isolate. Despite the high identity, analysis of individual differences among the 32 million base pair genomes showed sufficient variation to allow the isolates to be clustered based on the primary sequence. A major source of variation detected was in chromosome somy, with only four of the 36 chromosomes being predominantly disomic in all 49 isolates examined. In contrast, chromosome 31 was predominantly tetrasomic/pentasomic, consistent with its regions of synteny on two different disomic chromosomes of Trypanosoma brucei. In the Brazilian isolates, evidence for recombination was detected in 27 of the 36 chromosomes. Clustering analyses suggested two populations, in which two of the five older isolates from the 1990s clustered with a majority of recent isolates. Overall the analyses do not suggest individual sequence variants account for differences in clinical outcome or adaptation to different hosts. For the first known time, DNA of isolates from asymptomatic subjects were sequenced. Of interest, these displayed lower diversity than isolates from symptomatic subjects, an observation that deserves further investigation with additional isolates from asymptomatic subjects.