Global genome diversity of the Leishmania donovani complex

Genomic and epidemiological evidence for the emergence of a L. infantum/L. donovani hybrid with unusual epidemiology in northern Italy

Leishmania (L.) infantum is one of the main causative agents of animal and human leishmaniasis across many endemic areas in South America, Europe, North Africa, and Asia. Despite its clinical significance, little is known about the genetic diversity of L. infantum circulating in a given endemic area. Here, we investigate this important open question by applying a comparative genomics approach to seven L. infantum isolates from different hosts and Italian regions, including the northern part of the country (Emilia-Romagna, RER), Sicily, and Sardinia, as an initial attempt to explore the breadth of parasite genetic heterogeneity in Italy. Additionally, microsatellite analysis was carried out to compare the isolates from RER with other 70 L. infantum strains from the same region as well as 65 strains belonging to the L. donovani complex from other countries. We revealed important karyotypic instability and identified strain-specific changes in gene dosage, which affected important virulence factors such as amastins and surface antigen-like proteins. Single nucleotide polymorphism-based clustering analysis of these genomes together with over 80 publicly available L. infantum and L. donovani genomes placed the Italian isolates into three geographically distinct clusters within the Mediterranean basin and uncovered three isolates clustering with putative L. infantum/L. donovani hybrids isolated in Cyprus. As judged by microsatellite profiling, these hybrid isolates are representative of a sub-population of parasites circulating in northern Italy that preferentially infect humans but not dogs. Our results place Italy at the crossroads of L. infantum infection in the Mediterranean and call attention to the public health risk represented by the introduction of non-European Leishmania species.IMPORTANCEThis study closes important knowledge gaps with respect to Leishmania (L.) infantum genetic heterogeneity in a given endemic country, as exemplified here for Italy, and reveals genetic hybridization as a main cause for re-emerging human leishmaniasis in northern Italy. The observed high diversity of Leishmania parasites on the Italian peninsula suggests different geographical origins, with genomic adaptation to various ecologies affecting both pathogenicity and transmission potential. This is documented by the discovery of a putative L. infantum/L. donovani hybrid strain, which has been shown to preferentially infect humans but not dogs. Our results provide important information to health authorities, which need to consider the public health risk represented by the introduction of new Leishmania species into EU countries due to population displacement or travel from countries where exotic/allochthonous parasite species are endemic.

Effect of the aggregated protein dye YAT2150 on Leishmania parasite viability

The problems associated with the drugs currently used to treat leishmaniasis, including resistance, toxicity, and the high cost of some formulations, call for the urgent identification of new therapeutic agents with novel modes of action. The aggregated protein dye YAT2150 has been found to be a potent antileishmanial compound, with a half-maximal inhibitory concentration (IC50) of approximately 0.5 µM against promastigote and amastigote stages of Leishmania infantum. The encapsulation in liposomes of YAT2150 significantly improved its in vitro IC50 to 0.37 and 0.19 µM in promastigotes and amastigotes, respectively, and increased the half-maximal cytotoxic concentration in human umbilical vein endothelial cells to >50 µM. YAT2150 became strongly fluorescent when binding intracellular protein deposits in Leishmania cells. This fluorescence pattern aligns with the proposed mode of action of this drug in the malaria parasite Plasmodium falciparum, the inhibition of protein aggregation. In Leishmania major, YAT2150 rapidly reduced ATP levels, suggesting an alternative antileishmanial mechanism. To the best of our knowledge, this first-in-class compound is the only one described so far having significant activity against both Plasmodium and Leishmania, thus being a potential drug for the treatment of co-infections of both parasites.

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.

Molecular Typing of Leishmania spp. Causing Tegumentary Leishmaniasis in Northeastern Italy, 2014-2020

Tegumentary leishmaniasis (TL) is endemic but neglected in southern Europe. Therefore, this study aimed to analyze the Leishmania strains causing TL cases in northeastern Italy, where an upsurge of TL cases has been observed in the last decade. Sections from 109 formalin-fixed and paraffin-embedded (FFPE) biopsies of skin and mucosal tissues were collected from TL cases in the selected area. Two DNA targets were amplified and sequenced: the ribosomal internal transcribed spacer 1 (ITS1) and the heat-shock protein 70 gene (hsp70). An in silico analysis was also performed on 149 genomes belonging to the Leishmania donovani complex. A total of 88 out of 109 (80.7%) samples from 83 TL cases were successfully typed by ITS1 and/or hsp70. ITS1 analysis identified L. infantum in 67 cases (91.8%), while L. major (n = 4, 5.5%) and L. tropica (n = 2, 2.7%) were detected in the remaining cases that were categorized as imported. Further, the hsp70 typing of 75 autochthonous cases showed the presence of eight distinct sequence variants belonging to the Leishmania donovani complex, with high genetic variability when compared to known L. infantum populations. In conclusion, our findings show that peculiar L. infantum variants are emerging in the novel focus on TL in northeastern Italy.

First report of cutaneous leishmaniasis caused by Leishmania donovani in Ethiopia

BACKGROUND

Leishmaniasis is a common neglected tropical disease in Ethiopia. Visceral leishmaniasis (VL) caused by Leishmania donovani presents in the lowlands, while cutaneous leishmaniasis (CL) affects people living in the highlands. Although CL is described as being caused by Leishmania aethiopica, there is also evidence of L. tropica and L. major isolated from a patient, sand flies and potential reservoirs. Information on species causing CL in Ethiopia is patchy, and no nation-wide study has ever been done. Understanding which species are causing CL in Ethiopia can have important implications for patient management and disease prevention.

METHODS

We analyzed stored routine samples and biobanked DNA isolates from previously conducted studies of CL patients from different centers in the north, center and south of Ethiopia. Species typing was performed using ITS-1 PCR with high-resolution melt (HRM) analysis, followed by HSP70 amplicon sequencing on a selection of the samples. Additionally, sociodemographic, clinical and laboratory data of patients were analyzed.

RESULTS

Of the 226 CL samples collected, the Leishmania species could be determined for 105 (45.5%). Leishmania aethiopica was identified in 101 (96.2%) samples from across the country. In four samples originating from Amhara region, northwestern Ethiopia, L. donovani was identified by ITS-1 HRM PCR, of which two were confirmed with HSP70 sequences. While none of these four patients had symptoms of VL, two originated from known VL endemic areas.

CONCLUSIONS

The majority of CL was caused by L. aethiopica, but CL due to L. tropica and L. major cannot be ruled out. Our study is the first to our knowledge to demonstrate CL patients caused by L. donovani in Ethiopia. This should spark future research to investigate where, how and to which extent such transmission takes place, how it differs genetically from L. donovani causing VL and whether such patients can be diagnosed and treated successfully with the currently available tools and drugs.

Microgramma vacciniifolia Frond Lectin: In Vitro Anti-leishmanial Activity and Immunomodulatory Effects Against Internalized Amastigote Forms of Leishmania amazonensis

PURPOSE

The treatment of leishmaniasis, an anthropozoonosis caused by Leishmania protozoa, is limited by factors, such as adverse effects, toxicity, and excessive cost, which has highlighted the importance of novel drugs. In this context, natural products have been considered as sources of antileishmanial agents. This study investigated the leishmanicidal activity of Microgramma vacciniifolia frond lectin (MvFL) on promastigotes and amastigotes of Leishmania amazonensis.

METHODS

The effects of MvFL on promastigote proliferation and macrophage infection by amastigotes were evaluated and mean inhibitory concentrations (IC50) were calculated. As a safety assessment, the hemolytic capacity of MvFL (6.25-200 µg/mL) against mouse and human erythrocytes was determined. Additionally, the ability of MvFL (6.25-100 µg/mL) to modulate lysosomal and phagocytic activities and the nitric oxide (NO) production by murine peritoneal macrophages was also investigated.

RESULTS

After 24 h, MvFL inhibited the proliferation of L. amazonensis promastigotes, with an IC50 of 88 µg/mL; however, hemolytic activity was not observed. MvFL also reduced macrophage infection by amastigotes with an IC50 of 52 µg/mL. Furthermore, treatment with MvFL reduced the number of amastigotes internalized by infected murine peritoneal macrophages by up to 68.9% within 48 h. At a concentration of 25 µg/mL, MvFL stimulated lysosomal activity of macrophages within 72 h, but did not alter phagocytic activity or induce NO production at any of the tested concentrations.

CONCLUSION

MvFL exerts antileishmanial activity and further studies are needed to assess its therapeutic potential in in vivo experimental models of leishmaniasis.

Shining the spotlight on the neglected: new high-quality genome assemblies as a gateway to understanding the evolution of Trypanosomatidae

BACKGROUND

Protists of the family Trypanosomatidae (phylum Euglenozoa) have gained notoriety as parasites affecting humans, domestic animals, and agricultural plants. However, the true extent of the group's diversity spreads far beyond the medically and veterinary relevant species. We address several knowledge gaps in trypanosomatid research by undertaking sequencing, assembly, and analysis of genomes from previously overlooked representatives of this protistan group.

RESULTS

We assembled genomes for twenty-one trypanosomatid species, with a primary focus on insect parasites and Trypanosoma spp. parasitizing non-human hosts. The assemblies exhibit sizes consistent with previously sequenced trypanosomatid genomes, ranging from approximately 18 Mb for Obscuromonas modryi to 35 Mb for Crithidia brevicula and Zelonia costaricensis. Despite being the smallest, the genome of O. modryi has the highest content of repetitive elements, contributing nearly half of its total size. Conversely, the highest proportion of unique DNA is found in the genomes of Wallacemonas spp., with repeats accounting for less than 8% of the assembly length. The majority of examined species exhibit varying degrees of aneuploidy, with trisomy being the most frequently observed condition after disomy.

CONCLUSIONS

The genome of Obscuromonas modryi represents a very unusual, if not unique, example of evolution driven by two antidromous forces: i) increasing dependence on the host leading to genomic shrinkage and ii) expansion of repeats causing genome enlargement. The observed variation in somy within and between trypanosomatid genera suggests that these flagellates are largely predisposed to aneuploidy and, apparently, exploit it to gain a fitness advantage. High heterogeneity in the genome size, repeat content, and variation in chromosome copy numbers in the newly-sequenced species highlight the remarkable genome plasticity exhibited by trypanosomatid flagellates. These new genome assemblies are a robust foundation for future research on the genetic basis of life cycle changes and adaptation to different hosts in the family Trypanosomatidae.

Leishmania donovani persistence and circulation causing cutaneous leishmaniasis in unusual-foci of Nepal

Cutaneous leishmaniasis cases have increased dramatically in recent years in Nepal. The study offers molecular identification of the Leishmania species using 40 patient's aspiration biopsy samples, targeting markers kinetoplast minicircle DNA (kDNA) and internal transcribed spacer-1 (ITS1). Among molecularly diagnosed 22 cutaneous leishmaniasis cases, L. donovani complex was identified in 13 instances and L. major in 9 cases. The ITS1 PCR was positive in 12 of the positive nested- kDNA PCR cases (12/22), confirming L. donovani complex in seven of the cases and L. major in five of the cases. In addition, the study conclude that concurrent occurrence of atypical cutaneous infections caused by L. donovani parasite in 59.1% of cases and typical cutaneous infections caused by L. major parasite in 40.9% of cases. A Phylogentic analaysis showed that the detected L. donovani species present null genetic distances from seven references of L. donovani, but slight differences between ITS1 sequences and not grouped into a significant monophyletic cluster.

Kinetoplast Genome of Leishmania spp. Is under Strong Purifying Selection

Instability is an intriguing characteristic of many protist genomes, and trypanosomatids are not an exception in this respect. Some regions of trypanosomatid genomes evolve fast. For instance, the trypanosomatid mitochondrial (kinetoplast) genome consists of fairly conserved maxicircle and minicircle molecules that can, nevertheless, possess high nucleotide substitution rates between closely related strains. Recent experiments have demonstrated that rapid laboratory evolution can result in the non-functionality of multiple genes of kinetoplast genomes due to the accumulation of mutations or loss of critical genomic components. An example of a loss of critical components is the reported loss of entire minicircle classes in Leishmania tarentolae during laboratory cultivation, which results in an inability to generate some correctly encoded genes. In the current work, we estimated the evolutionary rates of mitochondrial and nuclear genome regions of multiple natural Leishmania spp. We analyzed synonymous and non-synonymous substitutions and, rather unexpectedly, found that the coding regions of kinetoplast maxicircles are among the most variable regions of both genomes. In addition, we demonstrate that synonymous substitutions greatly predominate among maxicircle coding regions and that most maxicircle genes show signs of purifying selection. These results imply that maxicircles in natural Leishmania populations remain functional despite their high mutation rate.

Species diversity of phlebotomine sand flies and sympatric occurrence of Leishmania (Mundinia) martiniquensis, Leishmania (Leishmania) donovani complex, and Trypanosoma spp. in the visceral leishmaniasis focus of southern Thailand

Autochthonous leishmaniasis in Thailand has recently been a public health concern due to an increasing number of new clinical cases. Most indigenous cases were diagnosed with Leishmania (Mundinia) martiniquensis, and Leishmania (Mundinia) orientalis. However, some doubts regarding vector misidentification have arisen and need to be elucidated. Accordingly, we aimed to assess the species composition of sand flies and determine the molecular prevalence of trypanosomatids in the transmission area of leishmaniasis in southern Thailand. In the present study, a total of 569 sand flies were caught from the vicinity of a visceral leishmaniasis patient's house in Na Thawi District, Songkhla Province. Of these, 229 parous and gravid females consisted of Sergentomyia khawi, Se. barraudi, Phlebotomus stantoni, Grassomyia indica, and Se. hivernus, accounting for 31.4%, 30.6%, 29.7%, 7.9%, and 0.4%, respectively. However, Se. gemmea, which has previously been proposed as the most abundant species and putative vector of visceral leishmaniasis, was not found in the present study. Based on ITS1-PCR and sequence analysis, two specimens of Gr. indica and Ph. stantoni showed positive amplification of L. martiniquensis and L. donovani complex, respectively, the first one being presumed indigenous and the second one being not. Anuran Trypanosoma was also molecularly detected using SSU rRNA-PCR and ubiquitously found in sixteen specimens of four dominant sand fly species except for Se. hivernus. The obtained sequences could be phylogenetically categorized into the two major amphibian clades (An04/Frog1 and An01+An02/Frog2). The existence of the monophyletic subgroup and distinct lineage suggests them as novel Trypanosoma species. The TCS network analysis of these anuran Trypanosoma sequences also revealed high haplotype diversity (Hd = 0.925 ± 0.050), but low nucleotide diversity (π = 0.019 ± 0.009). Furthermore, the living anuran trypanosomes were microscopically demonstrated in a single specimen of Gr. indica, supporting the vectorial capacity. Importantly, our data confirmed the scarcity of Se. gemmea and also uncovered, for the first time, the co-circulation of L. martiniquensis, L. donovani complex, and suspected novel anuran Trypanosoma spp. in phlebotomine sand files, implicating their potential role as vectors of trypanosomatid parasites. Therefore, the novel data from this study would greatly facilitate the comprehension of the complexity of trypanosomatid transmission and the establishment of prevention and control measures for this neglected disease more effectively.

Drug resistance in Leishmania: does it really matter?

Treatment failure (TF) jeopardizes the management of parasitic diseases, including leishmaniasis. From the parasite's point of view, drug resistance (DR) is generally considered as central to TF. However, the link between TF and DR, as measured by in vitro drug susceptibility assays, is unclear, some studies revealing an association between treatment outcome and drug susceptibility, others not. Here we address three fundamental questions aiming to shed light on these ambiguities. First, are the right assays being used to measure DR? Second, are the parasites studied, which are generally those that adapt to in vitro culture, actually appropriate? Finally, are other parasite factors - such as the development of quiescent forms that are recalcitrant to drugs - responsible for TF without DR?

Parasite hybridization promotes spreading of endosymbiotic viruses

Viruses are the most abundant biological entities on Earth and play a significant role in the evolution of many organisms and ecosystems. In pathogenic protozoa, the presence of endosymbiotic viruses has been linked to an increased risk of treatment failure and severe clinical outcome. Here, we studied the molecular epidemiology of the zoonotic disease cutaneous leishmaniasis in Peru and Bolivia through a joint evolutionary analysis of Leishmania braziliensis parasites and their endosymbiotic Leishmania RNA virus. We show that parasite populations circulate in isolated pockets of suitable habitat and are associated with single viral lineages that appear in low prevalence. In contrast, groups of hybrid parasites were geographically and ecologically dispersed, and commonly infected from a pool of genetically diverse viruses. Our results suggest that parasite hybridization, likely due to increased human migration and ecological perturbations, increased the frequency of endosymbiotic interactions known to play a key role in disease severity.

Parasite Genotype Is a Major Predictor of Mortality from Visceral Leishmaniasis

Visceral leishmaniasis (VL) is a potentially fatal disease caused mainly by Leishmania infantum in South America and Leishmania donovani in Asia and Africa. Disease outcomes have been associated with patient genotype, nutrition, age, sex, comorbidities, and coinfections. In this study, we examine the effects of parasite genetic variation on VL disease severity in Brazil. We collected and sequenced the genomes of 109 L. infantum isolates from patients in northeastern Brazil and retrieved matching patient clinical data from medical records, including mortality, sex, HIV coinfection, and laboratory data (creatinine, hemoglobin, and leukocyte and platelet counts). We identified genetic differences between parasite isolates, including single nucleotide polymorphisms (SNPs), small insertions/deletions (indels), and variations in genic, intergenic, and chromosome copy numbers (copy number variants [CNVs]). To describe associations between the parasite genotypes and clinical outcomes, we applied quantitative genetics methods of heritability and genome-wide association studies (GWAS), treating clinical outcomes as traits that may be influenced by parasite genotype. Multiple aspects of the genetic analysis indicate that parasite genotype affects clinical outcomes. We estimate that parasite genotype explains 83% chance of mortality (narrow-sense heritability [h²] = 0.83 ± 0.17) and has a significant relationship with patient sex (h² = 0.60 ± 0.27). Impacts of parasite genotype on other clinical traits are lower (h² ≤ 0.34). GWAS analysis identified multiple parasite genetic loci that were significantly associated with clinical outcomes; 17 CNVs were significantly associated with mortality, two with creatinine, and one with bacterial coinfection, jaundice, and HIV coinfection, and two SNPs/indels and six CNVs were associated with age, jaundice, HIV and bacterial coinfections, creatinine, and/or bleeding sites. Parasite genotype is an important factor in VL disease severity in Brazil. Our analysis indicates that specific genetic differences between parasites act as virulence factors, enhancing risks of severe disease and mortality. More detailed understanding of these virulence factors could be exploited for novel therapies. IMPORTANCE Multiple factors contribute to the risk of mortality from visceral leishmaniasis (VL), including, patient genotype, comorbidities, and nutrition. Many of these factors are influenced by socioeconomic biases. Our work suggests that the virulence of the infecting parasite is an important risk factor for mortality. We pinpoint some specific genomic markers that are associated with mortality, which can lead to a greater understanding of the molecular mechanisms that cause severe VL disease, to the identification of genetic markers for virulent parasites, and to the development of drug and vaccine therapies.

Self-Hybridization in Leishmania major

Genetic exchange between different Leishmania strains in the sand fly vector has been experimentally demonstrated and is supported by population genetic studies. In nature, opportunities for Leishmania interstrain mating are restricted to flies biting multiply infected hosts or through multiple bites of different hosts. In contrast, self-mating could occur in any infected sand fly. By crossing two recombinant lines derived from the same Leishmania major strain, each expressing a different drug-resistance marker, self-hybridization in L. major was confirmed in a natural sand fly vector, Phlebotomus duboscqi, and in frequencies comparable to interstrain crosses. We provide the first high resolution, whole-genome sequencing analysis of large numbers of selfing progeny, their parents, and parental subclones. Genetic exchange consistent with classical meiosis is supported by the biallelic inheritance of the rare homozygous single nucleotide polymorphisms (SNPs) that arose by mutation during the generation of the parental clones. In contrast, heterozygous SNPs largely failed to be transmitted in Mendelian ratios for reasons not understood. SNPs that were heterozygous in both parents, however, recombined to produce homozygous alleles in some hybrids. For trisomic chromosomes present in both parents, transmittal to the progeny was only altered by self-hybridization, involving a gain or loss of somy in frequencies predicted by a meiotic process. Whole-genome polyploidization was also observed in the selfing progeny. Thus, self-hybridization in Leishmania, with its potential to occur in any infected sand fly, may be an important source of karyotype variation, loss of heterozygosity, and functional diversity. IMPORTANCE Leishmania are parasitic protozoa that cause a wide spectrum of diseases collectively known as the leishmaniases. Sexual reproduction in Leishmania has been proposed as an important source of genetic diversity and has been formally demonstrated to occur inside the sand fly vector midgut. Nevertheless, in the wild, opportunities for genetic exchange between different Leishmania species or strains are restricted by the capacity of different Leishmania strains to colonize the same sand fly. In this work, we report the first high resolution, whole-genome sequence analysis of intraclonal genetic exchange as a type of self-mating in Leishmania. Our data reveal that self-hybridization can occur with comparable frequency as interstrain mating under experimental lab conditions, leading to important genomic alterations that can potentially take place within every naturally infected sand fly.

An overview of the trypanosomatid (Kinetoplastida: Trypanosomatidae) parasites infecting several mammal species in Colombia

BACKGROUND

Trypanosomatids are among the most critical parasites for public health due to their impact on human, animal, and plant health. Diseases associated with these pathogens manifest mainly in poor and vulnerable populations, where social, environmental, and biological factors modulate the case incidence and geographical distribution.

METHODS

We used Sanger and amplicon-based next-generation sequencing (NGS) in samples from different mammals to identify trypanosomatid infections in several departments in Colombia. A total of 174 DNA samples (18 humans, 83 dogs, and 73 wild mammals) were analyzed by conventional PCR using a fragment of the heat shock protein 70 (Hsp70) gene and Sanger sequenced the positive samples. Twenty-seven samples were sent for amplicon-based NGS using the same gene fragment. Data obtained were used to perform diversity analyses.

RESULTS

One hundred and thirteen samples were positive for PCR by Hsp70 fragment; these corresponded to 22.1% Leishmania spp., 18.6% L. amazonensis, 9.7% L. braziliensis, 14.2% L. infantum, 8% L. panamensis, and 27.4% Trypanosoma cruzi. Comparison of the identified species by the two sequencing technologies used resulted in 97% concordance. Alpha and beta diversity indices were significant, mainly for dogs; there was an interesting index of coinfection events in the analyzed samples: different Leishmania species and the simultaneous presence of T. cruzi and even T. rangeli in one of the samples analyzed. Moreover, a low presence of L. braziliensis was observed in samples from wild mammals. Interestingly, to our knowledge, this is the first report of Leishmania detection in Hydrochaeris hydrochaeris (capybara) in Colombia.

CONCLUSIONS

The Hsp70 fragment used in this study is an optimal molecular marker for trypanosomatid identification in many hosts and allows the identification of different species in the same sample when amplicon-based sequencing is used. However, the use of this fragment for molecular diagnosis through conventional PCR should be carefully interpreted because of this same capacity to identify several parasites. This point is of pivotal importance in highly endemic countries across South America because of the co-circulation of different genera from the Trypanosomatidae family. The findings show an interesting starting point for One Health approaches in which coevolution and vector-host interactions can be studied.

Repurposing the Antibacterial Agents Peptide 19-4LF and Peptide 19-2.5 for Treatment of Cutaneous Leishmaniasis

The lack of safe and cost-effective treatments against leishmaniasis highlights the urgent need to develop improved leishmanicidal agents. Antimicrobial peptides (AMPs) are an emerging category of therapeutics exerting a wide range of biological activities such as anti-bacterial, anti-fungal, anti-parasitic and anti-tumoral. In the present study, the approach of repurposing AMPs as antileishmanial drugs was applied. The leishmanicidal activity of two synthetic anti-lipopolysaccharide peptides (SALPs), so-called 19-2.5 and 19-4LF was characterized in Leishmania major. In vitro, both peptides were highly active against intracellular Leishmania major in mouse macrophages without exerting toxicity in host cells. Then, q-PCR-based gene profiling, revealed that this activity was related to the downregulation of several genes involved in drug resistance (yip1), virulence (gp63) and parasite proliferation (Cyclin 1 and Cyclin 6). Importantly, the treatment of BALB/c mice with any of the two AMPs caused a significant reduction in L. major infective burden. This effect was associated with an increase in Th1 cytokine levels (IL-12p35, TNF-α, and iNOS) in the skin lesion and spleen of the L. major infected mice while the Th2-associated genes were downregulated (IL-4 and IL-6). Lastly, we investigated the effect of both peptides in the gene expression profile of the P2X7 purinergic receptor, which has been reported as a therapeutic target in several diseases. The results showed significant repression of P2X7R by both peptides in the skin lesion of L. major infected mice to an extent comparable to that of a common anti-leishmanial drug, Paromomycin. Our in vitro and in vivo studies suggest that the synthetic AMPs 19-2.5 and 19-4LF are promising candidates for leishmaniasis treatment and present P2X7R as a potential therapeutic target in cutaneous leishmaniasis (CL).

Identification of Leishmania infantum Epidemiology, Drug Resistance and Pathogenicity Biomarkers with Nanopore Sequencing

The emergence of drug-resistant strains of the parasite Leishmania infantum infecting dogs and humans represents an increasing threat. L. infantum genomes are complex and unstable with extensive structural variations, ranging from aneuploidies to multiple copy number variations (CNVs). These CNVs have recently been validated as biomarkers of Leishmania concerning virulence, tissue tropism, and drug resistance. As a proof-of-concept to develop a novel diagnosis platform (LeishGenApp), four L. infantum samples from humans and dogs were nanopore sequenced. Samples were epidemiologically typed within the Mediterranean L. infantum group, identifying members of the JCP5 and non-JCP5 subgroups, using the conserved region (CR) of the maxicircle kinetoplast. Aneuploidies were frequent and heterogenous between samples, yet only chromosome 31 tetrasomy was common between all the samples. A high frequency of aneuploidies was observed for samples with long passage history (MHOM/TN/80/IPT-1), whereas fewer were detected for samples maintained in vivo (MCRI/ES/2006/CATB033). Twenty-two genes were studied to generate a genetic pharmacoresistance profile against miltefosine, allopurinol, trivalent antimonials, amphotericin, and paromomycin. MHOM/TN/80/IPT-1 and MCRI/ES/2006/CATB033 displayed a genetic profile with potential resistance against miltefosine and allopurinol. Meanwhile, MHOM/ES/2016/CATB101 and LCAN/ES/2020/CATB102 were identified as potentially resistant against paromomycin. All four samples displayed a genetic profile for resistance against trivalent antimonials. Overall, this proof-of-concept revealed the potential of nanopore sequencing and LeishGenApp for the determination of epidemiological, drug resistance, and pathogenicity biomarkers in L. infantum.

Livestock and rodents within an endemic focus of Visceral Leishmaniasis are not reservoir hosts for Leishmania donovani

Leishmaniasis on the Indian subcontinent is thought to have an anthroponotic transmission cycle. There is no direct evidence that a mammalian host other than humans can be infected with Leishmania donovani and transmit infection to the sand fly vector. The aim of the present study was to evaluate the impact of sand fly feeding on other domestic species and provide clinical evidence regarding possible non-human reservoirs through experimental sand fly feeding on cows, water buffalo goats and rodents. We performed xenodiagnosis using colonized Phlebotomus argentipes sand flies to feed on animals residing in villages with active Leishmania transmission based on current human cases. Xenodiagnoses on mammals within the endemic area were performed and blood-fed flies were analyzed for the presence of Leishmania via qPCR 48hrs after feeding. Blood samples were also collected from these mammals for qPCR and serology. Although we found evidence of Leishmania infection within some domestic mammals, they were not infectious to vector sand flies. Monitoring infection in sand flies and non-human blood meal sources in endemic villages leads to scientific proof of exposure and parasitemia in resident mammals. Lack of infectiousness of these domestic mammals to vector sand flies indicates that they likely play no role, or a very limited role in Leishmania donovani transmission to people in Bihar. Therefore, a surveillance system in the peri-/post-elimination phase of visceral leishmaniasis (VL) must monitor absence of transmission. Continued surveillance of domestic mammals in outbreak villages is necessary to ensure that a non-human reservoir is not established, including domestic mammals not present in this study, specifically dogs.

An experimental challenge model for Leishmania donovani in beagle dogs, showing a similar pattern of parasite burden in the peripheral blood and liver

Leishmania donovani and Leishmania infantum are closely related species. However, the former is considered the causative agent for anthroponotic visceral leishmaniasis (AVL), while the latter is known to be responsible for zoonotic visceral leishmaniasis (ZVL) with dogs as the main reservoir host. Although molecular detection of L. donovani from naturally infected dogs has been reported in AVL endemic areas, the experimental infection of dogs with this species is very limited. Here, we constructed an experimental canine visceral leishmaniasis (CVL) model with L. donovani infection using beagle dogs. During an observation period of 8 months after parasite inoculation, few clinical symptoms were observed in the three inoculated dogs. The overall hematological and biochemical data of the dogs showed normal levels, and there were no remarkable changes in the peripheral CD4⁺, CD8⁺, CD25⁺, or FoxP3⁺ T cell populations. Liver biopsy sampling was conducted to monitor the parasite burden in the liver. A similar pattern of the amount of mitochondrial kinetoplast DNA was observed in the peripheral blood and liver by real-time PCR analysis. In addition, parasite antigens were detected from the liver biopsy sections by immunohistochemical analysis, further supporting the existence of parasites in the liver. These results showed a subclinical CVL model for L. donovani in beagle dogs with a similar kinetics of parasite burden in the peripheral blood and liver.

Four layer multi-omics reveals molecular responses to aneuploidy in Leishmania

Aneuploidy causes system-wide disruptions in the stochiometric balances of transcripts, proteins, and metabolites, often resulting in detrimental effects for the organism. The protozoan parasite Leishmania has an unusually high tolerance for aneuploidy, but the molecular and functional consequences for the pathogen remain poorly understood. Here, we addressed this question in vitro and present the first integrated analysis of the genome, transcriptome, proteome, and metabolome of highly aneuploid Leishmania donovani strains. Our analyses unambiguously establish that aneuploidy in Leishmania proportionally impacts the average transcript- and protein abundance levels of affected chromosomes, ultimately correlating with the degree of metabolic differences between closely related aneuploid strains. This proportionality was present in both proliferative and non-proliferative in vitro promastigotes. However, as in other Eukaryotes, we observed attenuation of dosage effects for protein complex subunits and in addition, non-cytoplasmic proteins. Differentially expressed transcripts and proteins between aneuploid Leishmania strains also originated from non-aneuploid chromosomes. At protein level, these were enriched for proteins involved in protein metabolism, such as chaperones and chaperonins, peptidases, and heat-shock proteins. In conclusion, our results further support the view that aneuploidy in Leishmania can be adaptive. Additionally, we believe that the high karyotype diversity in vitro and absence of classical transcriptional regulation make Leishmania an attractive model to study processes of protein homeostasis in the context of aneuploidy and beyond.

Genomic diversity and genetic variation of Leishmania panamensis within its endemic range

Species belonging to the Leishmania (Viannia) subgenus are important causative agents of cutaneous and mucocutaneous leishmaniasis in Central and South America. These parasites possess several distinctive biological features that are influenced by their genetics, population structure, and genome instability. To date, several studies have revealed varying degrees of genetic diversity within Leishmania species. Particularly, in species of the L. (Viannia) subgenus, a generalized high intraspecific genetic diversity has been reported, although, conflicting conclusions have been drawn using different molecular techniques. Despite being the most common Leishmania species circulating in Panama and Colombia, few studies have analyzed clinical samples of Leishmania panamensis using whole-genome sequencing, and their restricted number of samples has limited the information they can provide to understand the population structure of L. panamensis. Here, we used next generation sequencing (NGS) to explore the genetic diversity of L. panamensis within its endemic range, analyzing data from 43 isolates of Colombian and Panamanian origin. Our results show the occurrence of three well-defined geographically correlated groups, and suggests the possible occurrence of additional phylogeographic groups. Furthermore, these results support the existence of a mixed mode of reproduction in L. panamensis, with varying frequencies of events of genetic recombination occurring primarily within subpopulations of closely related strains. This study offers important insights into the population genetics and reproduction mode of L. panamensis, paving the way to better understand their population structure and the emergence and maintenance of key eco-epidemiological traits.

Impact of Visceral Leishmaniasis on Local Organ Metabolism in Hamsters

Leishmania is an intracellular parasite with different species pathogenic to humans and causing the disease leishmaniasis. Leishmania donovani causes visceral leishmaniasis (VL) that manifests as hepatosplenomegaly, fever, pancytopenia and hypergammaglobulinemia. If left without treatment, VL can cause death, especially in immunocompromised people. Current treatments have often significant adverse effects, and resistance has been reported in some countries. Determining the metabolites perturbed during VL can lead us to find new treatments targeting disease pathogenesis. We therefore compared metabolic perturbation between L. donovani-infected and uninfected hamsters across organs (spleen, liver, and gut). Metabolites were extracted, analyzed by liquid chromatography-mass spectrometry, and processed with MZmine and molecular networking to annotate metabolites. We found few metabolites commonly impacted by infection across all three sites, including glycerophospholipids, ceramides, acylcarnitines, peptides, purines and amino acids. In accordance with VL symptoms and parasite tropism, we found a greater overlap of perturbed metabolites between spleen and liver compared to spleen and gut, or liver and gut. Targeting pathways related to these metabolite families would be the next focus that can lead us to find more effective treatments for VL.

Assembly of a Large Collection of Maxicircle Sequences and Their Usefulness for Leishmania Taxonomy and Strain Typing

Parasites of medical importance, such as Leishmania and Trypanosoma, are characterized by the presence of thousands of circular DNA molecules forming a structure known as kinetoplast, within the mitochondria. The maxicircles, which are equivalent to the mitochondrial genome in other eukaryotes, have been proposed as a promising phylogenetic marker. Using whole-DNA sequencing data, it is also possible to assemble maxicircle sequences as shown here and in previous works. In this study, based on data available in public databases and using a bioinformatics workflow previously reported by our group, we assembled the complete coding region of the maxicircles for 26 prototypical strains of trypanosomatid species. Phylogenetic analysis based on this dataset resulted in a robust tree showing an accurate taxonomy of kinetoplastids, which was also able to discern between closely related Leishmania species that are usually difficult to discriminate by classical methodologies. In addition, we provide a dataset of the maxicircle sequences of 60 Leishmania infantum field isolates from America, Western Europe, North Africa, and Eastern Europe. In agreement with previous studies, our data indicate that L. infantum parasites from Brazil are highly homogeneous and closely related to European strains, which were transferred there during the discovery of America. However, this study showed the existence of different L. infantum populations/clades within the Mediterranean region. A maxicircle signature for each clade has been established. Interestingly, two L. infantum clades were found coexisting in the same region of Spain, one similar to the American strains, represented by the Spanish JPCM5 reference strain, and the other, named “non-JPC like”, may be related to an important leishmaniasis outbreak that occurred in Madrid a few years ago. In conclusion, the maxicircle sequence emerges as a robust molecular marker for phylogenetic analysis and species typing within the kinetoplastids, which also has the potential to discriminate intraspecific variability.

An In-depth Proteomic Map of Leishmania donovani Isolate from Post Kala-azar Dermal Leishmaniasis (PKDL) Patient

BACKGROUND

The trypanosomatid protozoan parasite Leishmania donovani is the etiological agent of visceral leishmaniasis (VL) or kala-azar. The patients that have undergone treatment may still harbor the parasite and in a small fraction of the patients the disease re-erupts in the form of post kala-azar dermal leishmaniasis (PKDL). PKDL is a pathological condition found to be intermediate between VL and complete cure of VL. The PKDL disease progression is determined by the host immune response to L. donovani. The majority of the proteomic studies on L. donovani till date have been undertaken on parasites either isolated from kala-azar patients or on established laboratory strains of L. donovani. However, no proteomic information is available on the cutaneous localized isolates of L. donovani from PKDL patients.

METHODS

The promastigote stage of L. donovani isolate from PKDL patient was cultured and harvested. The cell lysates were trypsin digested, followed by liquid chromatography-tandem mass spectrometry (LC-MS/MS) analysis. The LC-MS/MS raw data were analyzed on Proteome Discoverer. Further bioinformatics analysis was carried out.

RESULTS

In the present, we have used high-resolution mass spectrometry to map the global proteome of a L. donovani isolate from PKDL patient. This in-depth study resulted in the identification of 5537 unique proteins from PKDL isolate of L. donovani which covered 64% of its proteome.

OUTCOME

This study also identified proteins previously shown to be upregulated in PKDL L. donovani. This is the most in-depth proteome of Leishmania donovani parasite till date.

Geographic Origin and Vertical Transmission of Leishmania infantum Parasites in Hunting Hounds, United States

Vertical transmission of leishmaniasis is common but is difficult to study against the background of pervasive vector transmission. We present genomic data from dogs in the United States infected with Leishmania infantum parasites; these infections have persisted in the apparent absence of vector transmission. We demonstrate that these parasites were introduced from the Old World separately and more recently than L. infantum from South America. The parasite population shows unusual genetics consistent with a lack of meiosis: a high level of heterozygous sites shared across all isolates and no decrease in linkage with genomic distance between variants. Our data confirm that this parasite population has been evolving with little or no sexual reproduction. This demonstration of vertical transmission has profound implications for the population genetics of Leishmania parasites. When investigating transmission in complex natural settings, considering vertical transmission alongside vector transmission is vital.

Development of a Highly Sensitive Nested PCR and Its Application for the Diagnosis of Cutaneous Leishmaniasis in Sri Lanka

The recent surge in cutaneous leishmaniasis (CL) in Sri Lanka has rendered clinical diagnosis difficult; thus, laboratory confirmation is indispensable. A modified (two novel inner primers to detect CL caused by Leishmania donovani) nested Internal Transcribed Spacer-1 (ITS1) PCR-Restriction Fragment Length Polymorphism (RFLP) method was developed and tested. The sensitivity of the modified nested PCR was tested using serial dilutions (103 to 10−2) of the DNA extract of a cultured L. donovani DD8 strain. Patients (n = 194) from Southern Sri Lanka were examined clinically, microscopically (Slit Skin Smear-SSS) and using the modified nested PCR. The modified nested PCR detected 2.55 fg of parasite DNA compared to ITS1 PCR (25 fg) and detected more cases than SSS (94.3% vs. 77.3%; p < 0.01). The RFLP pattern was L. donovani in all cases. The modified nested PCR performed well in clinically doubtful lesions (95% by PCR vs. 60% by SSS; p < 0.01), ulcerated nodules (91% vs. 71.8%; p < 0.01) and plaques (100% vs. 66.7%; p < 0.01). SSS demonstrated sensitivity (80.9%), specificity (81.8%), PPV (98.7%) and NPV (20.5%) against modified PCR. Low parasite loads and atypical lesions can be diagnosed by the proposed method with higher accuracy.

Validation of an In-House ELISA Method in the Diagnosis of Cutaneous Leishmaniasis Caused by Leishmania donovani in Hambantota District, Sri Lanka

Clinical diagnosis has become a challenge amidst a surge of cutaneous leishmaniasis in Southern Sri Lanka. The routine diagnostic method, slit-skin smear (SSS), has variable sensitivity, leading to undiagnosed cases. Improved diagnostics are urgently needed. We assessed a new in-house ELISA method for its diagnostic capabilities against ITS-1 nested PCR (gold standard—Gs). A cohort of 190 clinical CL cases was examined by SSS microscopy, anti-rKRP42 IgG ELISA (serum- and urine-based), and rK39-Immunochromatographic strip test. Validation was done using non-endemic sera, and cutoffs were developed using the receiver operating curve. The sensitivity of SSS for case detection was 77.9% (authors) and 76.3% (technicians). ELISA vs. Gs demonstrated sensitivity (Sn) = 94.4%; specificity (Sp) = 50.0%; positive predictive value (PPV) = 97.1%; negative predictive value (NPV) = 33.3%; Kappa agreement (Kp) = 0.39/p < 0.01. Comparison of the combination method (SSS by technicians and ELISA) vs. Gs showed: Sn = 98.9%; Sp = 30.0; PPV = 96.2; NPV 60.0%; Kp = 0.378/p < 0.01. All methods performed better compared to SSS (29.4%) where the clinical diagnosis was doubtful (PCR = 94.15%; serum ELISA = 88.2%; combination = 94.1%; p < 0.01 for all). High serum anti-rKRP42 titers were seen in those with multiple lesions. Anti-rKRP42 urine ELISA was suboptimal as a diagnostic test. A 9% rate of positivity was seen for rk39-ICT, and positives recorded high anti-rKRP42 titers. The diagnostic accuracy can be increased above the level of the Gs by combining SSS and ELISA. Advanced studies are required to understand the association between rk39-ICT positivity and high anti-rKRP42 titers.

Tissue Specific Dual RNA-Seq Defines Host-Parasite Interplay in Murine Visceral Leishmaniasis Caused by Leishmania donovani and Leishmania infantum

Visceral leishmaniasis is associated with hepato-splenomegaly and altered immune and hematological parameters in both preclinical animal models and humans. We studied mouse experimental visceral leishmaniasis caused by Leishmania infantum and Leishmania donovani in BALB/c mice using dual RNA-seq to investigate the transcriptional response of host and parasite in liver and spleen. We identified only 4 species-specific parasite expressed genes (SSPEGs; log2FC >1, FDR <0.05) in the infected spleen, and none in the infected liver. For the host transcriptome, we found 789 differentially expressed genes (DEGs; log2FC >1, FDR <0.05) in the spleen that were common to both infections, with IFNγ signaling and complement and coagulation cascade pathways highly enriched, and an additional 286 and 186 DEGs that were selective to L. donovani and L. infantum infection, respectively. Among those, there were network interactions between genes of amino acid metabolism and PPAR signaling in L. donovani infection and increased IL1β and positive regulation of fatty acid transport in L. infantum infection, although no pathway enrichment was observed. In the liver, there were 1,939 DEGs in mice infected with either L. infantum or L. donovani in comparison to uninfected mice, and the most enriched pathways were IFNγ signaling, neutrophil mediated immunity, complement and coagulation, cytokine-chemokine responses, and hemostasis. Additionally, 221 DEGs were selective in L. donovani and 429 DEGs in L. infantum infections. These data show that the host response for these two visceral leishmaniasis infection models is broadly similar, and ∼10% of host DEGs vary in infections with either parasite species. IMPORTANCE Visceral leishmaniasis (VL) is caused by two species of Leishmania parasites, L. donovani in the Old World and L. infantum in the New World and countries bordering the Mediterranean. Although cardinal features such as hepato-splenomegaly and alterations in blood and immune function are evident, clinical presentation may vary by geography, with for example severe bleeding often associated with VL in Brazil. Although animal models of both L. donovani and L. infantum have been widely used to study disease pathogenesis, a direct side-by-side comparison of how these parasites species impact the infected host and/or how they might respond to the stresses of mammalian infection has not been previously reported. Identifying common and distinct pathways to pathogenesis will be important to ensure that new therapeutic or prophylactic approaches will be applicable across all forms of VL.

Leishmania infantum Defective in Lipophosphoglycan Biosynthesis Interferes With Activation of Human Neutrophils

Visceral leishmaniasis (VL) is often associated with hematologic manifestations that may interfere with neutrophil response. Lipophosphoglycan (LPG) is a major molecule on the surface of Leishmania promastigotes, which has been associated with several aspects of the parasite–vector–host interplay. Here, we investigated how LPG from Leishmania (L.) infantum, the principal etiological agent of VL in the New World, influences the initial establishment of infection during interaction with human neutrophils in an experimental setting in vitro. Human neutrophils obtained from peripheral blood samples were infected with either the wild-type L. infantum (WT) strain or LPG-deficient mutant (∆lpg1). In this setting, ∆lpg1 parasites displayed reduced viability compared to WT L. infantum; such finding was reverted in the complemented ∆lpg1+LPG1 parasites at 3- and 6-h post-infection. Confocal microscopy experiments indicated that this decreased survival was related to enhanced lysosomal fusion. In fact, LPG-deficient L. infantum parasites more frequently died inside neutrophil acidic compartments, a phenomenon that was reverted when host cells were treated with Wortmannin. We also observed an increase in the secretion of the neutrophil collagenase matrix metalloproteinase-8 (MMP-8) by cells infected with ∆lpg1 L. infantum compared to those that were infected with WT parasites. Furthermore, collagen I matrix degradation was found to be significantly increased in ∆lpg1 parasite-infected cells but not in WT-infected controls. Flow cytometry analysis revealed a substantial boost in production of reactive oxygen species (ROS) during infection with either WT or ∆lpg1 L. infantum. In addition, killing of ∆lpg1 parasites was shown to be more dependent on the ROS production than that of WT L. infantum. Notably, inhibition of the oxidative stress with Apocynin potentially fueled ∆lpg1 L. infantum fitness as it increased the intracellular parasite viability. Thus, our observations demonstrate that LPG may be a critical molecule fostering parasite survival in human neutrophils through a mechanism that involves cellular activation and generation of free radicals.

In-Depth Quantitative Proteomics Characterization of In Vitro Selected Miltefosine Resistance in Leishmania infantum

Visceral leishmaniasis (VL) is a neglected disease caused by Leishmania parasites. Although significant morbidity and mortality in tropical and subtropical regions of the world are associated with VL, the low investment for developing new treatment measures is chronic. Moreover, resistance and treatment failure are increasing for the main medications, but the emergence of resistance phenotypes is poorly understood at the protein level. Here, we analyzed the development of resistance to miltefosine upon experimental selection in a L. infantum strain. Time to miltefosine resistance emergence was ~six months and label-free quantitative mass-spectrometry-based proteomics analyses revealed that this process involves a remodeling of components of the membrane and mitochondrion, with significant increase in oxidative phosphorylation complexes, particularly on complex IV and ATP synthase, accompanied by increased energy metabolism mainly dependent on β-oxidation of fatty acids. Proteins canonically involved in ROS detoxification did not contribute to the resistant process whereas sterol biosynthesis enzymes could have a role in this development. Furthermore, changes in the abundance of proteins known to be involved in miltefosine resistance such as ABC transporters and phospholipid transport ATPase were detected. Together, our data show a more complete picture of the elements that make up the miltefosine resistance phenotype in L. infantum.

Domestic mammals as reservoirs for Leishmania donovani on the Indian subcontinent: Possibility and consequences on elimination

Leishmania donovani is the causative agent of historically anthroponotic visceral leishmaniasis (VL) on the Indian subcontinent (ISC). L. donovani is transmitted by the sand fly species Phlebotomus argentipes. Our collaborative group and others have shown that sand flies trapped outside in endemic villages have fed on cattle and dogs in addition to people. Domestic animals are reservoirs for L. donovani complex spp., particularly L. infantum, in other endemic areas. Multiple studies using quantitative PCR or serological detection methods have demonstrated that goats, cattle, rats and dogs were diagnostically positive for L. donovani infection or exposure in eastern Africa, Bangladesh, Nepal and India. There is a limited understanding of the extent to which L. donovani infection of domestic animals drives transmission to other animals or humans on the ISC. Evidence from other vector-borne disease elimination strategies indicated that emerging infections in domestic species hindered eradication. The predominant lesson learned from these other situations is that non-human reservoirs must be identified, controlled and/or prevented. Massive efforts are underway for VL elimination on the Indian subcontinent. Despite these herculean efforts, residual VL incidence persists. The spectre of an animal reservoir complicating elimination efforts haunts the final push towards full VL control. Better understanding of L. donovani transmission on the Indian subcontinent and rigorous consideration of how non-human reservoirs alter VL ecology are critical to sustain elimination goals.

Culture of Cutaneous Leishmania from Skin Biopsy Specimens

This protocol describes the culture of Leishmania parasites from skin biopsy samples of patients with cutaneous lesions. The use of antibiotics to prevent bacterial contamination of these cultures increases the ability of researchers to collect isolates for various research purposes, including genetic analysis and in vitro and in vivo experiments. © Published 2022. This article is a U.S. Government work and is in the public domain in the USA. Basic Protocol: Culture of Leishmania from skin biopsy specimens.

Surveillance of leishmaniasis cases from 15 European centres, 2014 to 2019: a retrospective analysis

Background

Surveillance of human leishmaniasis in Europe is mostly limited to country-specific information from autochthonous infections in the southern part. As at the end of 2021, no integrated analysis has been performed for cases seen across centres in different European countries.

Aim

To provide a broad perspective on autochthonous and imported leishmaniasis cases in endemic and non-endemic countries in Europe.

Methods

We retrospectively collected records from cutaneous, mucosal and visceral leishmaniasis cases diagnosed in 15 centres between 2014 and 2019. Centres were located in 11 countries: Belgium, France, Germany, Italy, the Netherlands, Norway, Portugal, Spain, Sweden, Switzerland and the United Kingdom. Data on country of infection, reason for travelling, infecting species, age and sex were analysed.

Results

We obtained diagnostic files from 1,142 cases, of which 76%, 21% and 3% had cutaneous, visceral, and mucosal disease, respectively. Of these, 68% were men, and 32% women, with the median age of 37 years (range: 0–90) at diagnosis. Visceral leishmaniasis was mainly acquired in Europe (88%; 167/190), while cutaneous leishmaniasis was primarily imported from outside Europe (77%; 575/749). Sixty-two percent of cutaneous leishmaniasis cases from outside Europe were from the Old World, and 38% from the New World. Geographic species distribution largely confirmed known epidemiology, with notable exceptions.

Conclusions

Our study confirms previous reports regarding geographic origin, species, and traveller subgroups importing leishmaniasis into Europe. We demonstrate the importance of pooling species typing data from many centres, even from areas where the aetiology is presumably known, to monitor changing epidemiology.

Impact of Genetic Diversity and Genome Plasticity of Leishmania spp. in Treatment and the Search for Novel Chemotherapeutic Targets

Leishmaniasis is one of the major public health concerns in Latin America, Africa, Asia, and Europe. The absence of vaccines for human use and the lack of effective vector control programs make chemotherapy the main strategy to control all forms of the disease. However, the high toxicity of available drugs, limited choice of therapeutic agents, and occurrence of drug-resistant parasite strains are the main challenges related to chemotherapy. Currently, only a small number of drugs are available for leishmaniasis treatment, including pentavalent antimonials (Sb^V), amphotericin B and its formulations, miltefosine, paromomycin sulphate, and pentamidine isethionate. In addition to drug toxicity, therapeutic failure of leishmaniasis is a serious concern. The occurrence of drug-resistant parasites is one of the causes of therapeutic failure and is closely related to the diversity of parasites in this genus. Owing to the enormous plasticity of the genome, resistance can occur by altering different metabolic pathways, demonstrating that resistance mechanisms are multifactorial and extremely complex. Genetic variability and genome plasticity cause not only the available drugs to have limitations, but also make the search for new drugs challenging. Here, we examined the biological characteristics of parasites that hinder drug discovery.

Minor Impact of A258D Mutation on Biochemical and Enzymatic Properties of Leishmania infantum GDP-Mannose Pyrophosphorylase

Background: Leishmaniasis, a vector-borne disease caused by the protozoan parasite from the genus Leishmania, is endemic to tropical and subtropical areas. Few treatments are available against leishmaniasis, with all presenting issues of toxicity, resistance, and/or cost. In this context, the development of new antileishmanial drugs is urgently needed. GDP-mannose pyrophosphorylase (GDP-MP), an enzyme involved in the mannosylation pathway, has been described to constitute an attractive therapeutic target for the development of specific antileishmanial agents. Methods: In this work, we produced, purified, and analyzed the enzymatic properties of the recombinant L. infantum GDP-MP (LiGDP-MP), a single leishmanial GDP-MP that presents mutation of an aspartate instead of an alanine at position 258, which is also the single residue difference with the homolog in L. donovani: LdGDP-MP. Results: The purified LiGDP-MP displayed high substrate and cofactor specificities, a sequential random mechanism of reaction, and the following kinetic constants: V_m at 0.6 µM·min⁻¹, K_m from 15–18 µM, k_cat from 12.5–13 min⁻¹, and k_cat/K_m at around 0.8 min⁻¹µM⁻¹. Conclusions: These results show that LiGDP-MP has similar biochemical and enzymatic properties to LdGDP-MP. Further studies are needed to determine the advantage for L. infantum of the A258D residue change in GDP-MP.

Stress conditions promote Leishmania hybridization in vitro marked by expression of the ancestral gamete fusogen HAP2 as revealed by single-cell RNA-seq

Leishmania are protozoan parasites transmitted by the bite of sand fly vectors producing a wide spectrum of diseases in their mammalian hosts. These diverse clinical outcomes are directly associated with parasite strain and species diversity. Although Leishmania reproduction is mainly clonal, a cryptic sexual cycle capable of producing hybrid genotypes has been inferred from population genetic studies and directly demonstrated by laboratory crosses. Experimentally, mating competence has been largely confined to promastigotes developing in the sand fly midgut. The ability to hybridize culture promastigotes in vitro has been limited so far to low-efficiency crosses between two Leishmania tropica strains, L747 and MA37, that mate with high efficiency in flies. Here, we show that exposure of promastigote cultures to DNA damage stress produces a remarkably enhanced efficiency of in vitro hybridization of the L. tropica strains and extends to other species, including Leishmania donovani, Leishmania infantum, and Leishmania braziliensis, a capacity to generate intra- and interspecific hybrids. Whole-genome sequencing and total DNA content analyses indicate that the hybrids are in each case full genome, mostly tetraploid hybrids. Single-cell RNA sequencing of the L747 and MA37 parental lines highlights the transcriptome heterogeneity of culture promastigotes and reveals discrete clusters that emerge post-irradiation in which genes potentially involved in genetic exchange are expressed, including the ancestral gamete fusogen HAP2. By generating reporter constructs for HAP2, we could select for promastigotes that could either hybridize or not in vitro. Overall, this work reveals that there are specific populations involved in Leishmania hybridization associated with a discernible transcriptomic signature, and that stress facilitated in vitro hybridization can be a transformative approach to generate large numbers of hybrid genotypes between diverse species and strains.

High throughput single-cell genome sequencing gives insights into the generation and evolution of mosaic aneuploidy in Leishmania donovani

Leishmania, a unicellular eukaryotic parasite, is a unique model for aneuploidy and cellular heterogeneity, along with their potential role in adaptation to environmental stresses. Somy variation within clonal populations was previously explored in a small subset of chromosomes using fluorescence hybridization methods. This phenomenon, termed mosaic aneuploidy (MA), might have important evolutionary and functional implications but remains under-explored due to technological limitations. Here, we applied and validated a high throughput single-cell genome sequencing method to study for the first time the extent and dynamics of whole karyotype heterogeneity in two clonal populations of Leishmania promastigotes representing different stages of MA evolution in vitro. We found that drastic changes in karyotypes quickly emerge in a population stemming from an almost euploid founder cell. This possibly involves polyploidization/hybridization at an early stage of population expansion, followed by assorted ploidy reduction. During further stages of expansion, MA increases by moderate and gradual karyotypic alterations, affecting a defined subset of chromosomes. Our data provide the first complete characterization of MA in Leishmania and pave the way for further functional studies.

Amplicon-based next-generation sequencing reveals the co-existence of multiple Leishmania species in patients with visceral leishmaniasis

Visceral leishmaniasis (VL) is a mammalian protozoal disease propagated in the Americas by female phlebotomine sandflies, mainly caused by Leishmania infantum. However, in recent years, cases of VL caused by different Leishmania species, such as L. amazonensis and L. colombiensis, have been reported in the continent. This study used an amplicon-based next-generation sequencing approach to identify VL aetiologic species using high-depth sequencing targeting a region on the Heat Shock Protein 70 gene. In this first approach, six samples from five patients diagnosed with VL were selected and analysed to identify DNA of Leishmania spp. All samples harboured DNA of L. infantum; five samples were found to be co-infected with other Leishmania spp. or with Trypanosoma cruzi, and just one sample was mono-infected with L. infantum. This study demonstrates the usefulness of this methodology to identify trypanosomatid co-infections in clinical samples, which presents an interesting study panorama considering their biological, clinical and epidemiological implications.

Candidates for Balancing Selection in Leishmania donovani Complex Parasites

The Leishmania donovani species complex is the causative agent of visceral leishmaniasis, which cause 20–40,000 fatalities a year. Here, we conduct a screen for balancing selection in this species complex. We used 384 publicly available L. donovani and L. infantum genomes, and sequence 93 isolates of L. infantum from Brazil to describe the global diversity of this species complex. We identify five genetically distinct populations that are sufficiently represented by genomic data to search for signatures of selection. We find that signals of balancing selection are generally not shared between populations, consistent with transient adaptive events, rather than long-term balancing selection. We then apply multiple diversity metrics to identify candidate genes with robust signatures of balancing selection, identifying a curated set of 24 genes with robust signatures. These include zeta toxin, nodulin-like, and flagellum attachment proteins. This study highlights the extent of genetic divergence between L. donovani complex parasites and provides genes for further study.

Leishmania and the Model of Predominant Clonal Evolution

As it is the case for other pathogenic microorganisms, the respective impact of clonality and genetic exchange on Leishmania natural populations has been the object of lively debates since the early 1980s. The predominant clonal evolution (PCE) model states that genetic exchange in these parasites’ natural populations may have a high relevance on an evolutionary scale, but is not sufficient to erase a persistent phylogenetic signal and the existence of bifurcating trees. Recent data based on high-resolution markers and genomic polymorphisms fully confirm the PCE model down to a microevolutionary level.

rKOMICS: an R package for processing mitochondrial minicircle assemblies in population-scale genome projects

Background

The advent of population-scale genome projects has revolutionized our biological understanding of parasitic protozoa. However, while hundreds to thousands of nuclear genomes of parasitic protozoa have been generated and analyzed, information about the diversity, structure and evolution of their mitochondrial genomes remains fragmentary, mainly because of their extraordinary complexity. Indeed, unicellular flagellates of the order Kinetoplastida contain structurally the most complex mitochondrial genome of all eukaryotes, organized as a giant network of homogeneous maxicircles and heterogeneous minicircles. We recently developed KOMICS, an analysis toolkit that automates the assembly and circularization of the mitochondrial genomes of Kinetoplastid parasites. While this tool overcomes the limitation of extracting mitochondrial assemblies from Next-Generation Sequencing datasets, interpreting and visualizing the genetic (dis)similarity within and between samples remains a time-consuming process.

Results

Here, we present a new analysis toolkit—rKOMICS—to streamline the analyses of minicircle sequence diversity in population-scale genome projects. rKOMICS is a user-friendly R package that has simple installation requirements and that is applicable to all 27 trypanosomatid genera. Once minicircle sequence alignments are generated, rKOMICS allows to examine, summarize and visualize minicircle sequence diversity within and between samples through the analyses of minicircle sequence clusters. We showcase the functionalities of the (r)KOMICS tool suite using a whole-genome sequencing dataset from a recently published study on the history of diversification of the Leishmania braziliensis species complex in Peru. Analyses of population diversity and structure highlighted differences in minicircle sequence richness and composition between Leishmania subspecies, and between subpopulations within subspecies.

Conclusion

The rKOMICS package establishes a critical framework to manipulate, explore and extract biologically relevant information from mitochondrial minicircle assemblies in tens to hundreds of samples simultaneously and efficiently. This should facilitate research that aims to develop new molecular markers for identifying species-specific minicircles, or to study the ancestry of parasites for complementary insights into their evolutionary history.

Allele-specific assembly of a eukaryotic genome corrects apparent frameshifts and reveals a lack of nonsense-mediated mRNA decay

To date, most reference genomes represent a mosaic consensus sequence in which the homologous chromosomes are collapsed into one sequence. This approach produces sequence artefacts and impedes analyses of allele-specific mechanisms. Here, we report an allele-specific genome assembly of the diploid parasite Trypanosoma brucei and reveal allelic variants affecting gene expression. Using long-read sequencing and chromosome conformation capture data, we could assign 99.5% of all heterozygote variants to a specific homologous chromosome and build a 66 Mb long allele-specific genome assembly. The phasing of haplotypes allowed us to resolve hundreds of artefacts present in the previous mosaic consensus assembly. In addition, it revealed allelic recombination events, visible as regions of low allelic heterozygosity, enabling the lineage tracing of T. brucei isolates. Interestingly, analyses of transcriptome and translatome data of genes with allele-specific premature termination codons point to the absence of a nonsense-mediated decay mechanism in trypanosomes. Taken together, this study delivers a reference quality allele-specific genome assembly of T. brucei and demonstrates the importance of such assemblies for the study of gene expression control. We expect the new genome assembly will increase the awareness of allele-specific phenomena and provide a platform to investigate them.

Revisiting the heterogeneous global genomic population structure of Leishmania infantum

Leishmania infantum is the main causative agent responsible for visceral leishmaniasis (VL), a disease with global distribution. The genomic structure and genetic variation of this species have been widely studied in different parts of the world. However, in some countries, this information is still yet unknown, as is the genomic behaviour of the main antigens used in VL diagnosis (rK39 and rK28), which have demonstrated variable sensitivity and specificity in a manner dependent on the geographic region analysed. The objective of this study was to explore the genomic architecture and diversity of four Colombian L. infantum isolates obtained in this study and to compare these results with the genetic analysis of 183 L. infantum isolates from across the world (obtained from public databases), as well as to analyse the whole rK39 and rK28 antigen sequences in our dataset. The results showed that, at the global level, L. infantum has high genetic homogeneity and extensive aneuploidy. Furthermore, we demonstrated that there are distinct populations of L. infantum circulating in various countries throughout the globe and that populations of distant countries have close genomic relationships. Additionally, this study demonstrated the high genetic variability of the rK28 antigen worldwide. In conclusion, our study allowed us to (i) expand our knowledge of the genomic structure of global L. infantum; (ii) describe the intra-specific genomic variability of this species; and (iii) understand the genomic characteristics of the main antigens used in the diagnosis of VL. Additionally, this is the first study to report whole-genome sequences of Colombian L. infantum isolates.

Spatial and Temporal Variability of Visceral Leishmaniasis in Colombia, 2007 to 2018

Visceral leishmaniasis (VL) is a neglected tropical disease associated with poverty and is endemic in 56 countries worldwide. Brazil, Venezuela, and Colombia are the most affected countries in South America. In Colombia, the National Public Health Surveillance System (SIVIGILA) consolidates epidemiological information and monitors all VL cases nationwide. However, to date, no studies have investigated the occurrence of VL in Colombia using metadata analysis. We studied the demographic data, the spatial and temporal distribution of VL cases, and the association with vector distribution of Leishmania species in Colombia from 2007 to 2018. We found 306 VL cases reported to SIVIGILA for this period, with a coverage of 25.5 cases/year, and a mortality of 2.28% (seven deaths). The highest number of confirmed cases (N = 52) occurred in 2007; the lowest (N = 9) occurred in 2012. The cases were reported mainly in children (< 7 years) affiliated with the subsidized health regimen. Regarding the geographic distribution, the cases were reported by 42 municipalities distributed in 10 departments. The occurrence of VL cases toward the northeast of Colombia, and the distribution of vectors, such as Lutzomyia longipalpis and Lu. evansi, may be changing the panorama of VL in the country. We conclude that VL, mainly in recent years, shows a temporal and spatial variability associated with the occurrence of cases in new settings. Our findings increase our understanding and knowledge of this disease, and suggest the need to monitor and prioritize areas with changes in geographic expansion to improve prevention and control actions in the country.

Genomics of Trypanosomatidae: Where We Stand and What Needs to Be Done?

Trypanosomatids are easy to cultivate and they are (in many cases) amenable to genetic manipulation. Genome sequencing has become a standard tool routinely used in the study of these flagellates. In this review, we summarize the current state of the field and our vision of what needs to be done in order to achieve a more comprehensive picture of trypanosomatid evolution. This will also help to illuminate the lineage-specific proteins and pathways, which can be used as potential targets in treating diseases caused by these parasites.

In vitro culture of freshly isolated Trypanosoma brucei brucei bloodstream forms results in gene copy-number changes

Most researchers who study unicellular eukaryotes work with an extremely limited number of laboratory-adapted isolates that were obtained from the field decades ago, but the effects of passage in laboratory rodents, and adaptation to in vitro culture, have been little studied. For example, the vast majority of studies of Trypanosoma brucei biology have concentrated on just two strains, Lister 427 and EATRO1125, which were taken from the field over half a century ago and have since have undergone innumerable passages in rodents and culture. We here describe two new Trypanosoma brucei brucei strains. MAK65 and MAK98, which have undergone only 3 rodent passages since isolation from Ugandan cattle. High-coverage sequencing revealed that adaptation of the parasites to culture was accompanied by changes in gene copy numbers. T. brucei has so far been considered to be uniformly diploid, but we also found trisomy of chromosome 5 not only in one Lister 427 culture, but also in the MAK98 field isolate. Trisomy of chromosome 6, and increased copies of other chromosome segments, were also seen in established cultured lines. The two new T. brucei strains should be useful to researchers interested in trypanosome differentiation and pathogenicity. Initial results suggested that the two strains have differing infection patterns in rodents. MAK65 is uniformly diploid and grew more reproducibly in bloodstream-form culture than MAK98.

Leishmania infantum strains from cats are similar in biological properties to canine and human strains

Zoonotic visceral leishmaniosis is a worldwide severe disease caused by Leishmania infantum, a protozoan that has phlebotomine sand flies as vectors and dogs as primary reservoir hosts. Over the last few decades, cats have been regarded as an indisputable piece within the ecological system in which L. infantum is maintained indefinitely. However, little is known about feline strains, including their phenotypic plasticity and infectivity. In this study, the phenotypic behaviour of seven L. infantum feline strains was compared to those of well-characterised counterparts isolated from two dogs and two humans in terms of growth profile, adaptive capacity under several stress conditions, susceptibility to antileishmanial drugs, and infectivity to host cells. Feline strains displayed a similar growth profile, survival capacity, and ability to infect feline, canine, and human monocyte-derived primary macrophages. Furthermore, multivariate cluster analysis suggested that most strains studied did not display distinctive phenotypic features. To our knowledge, this is the first study to analyse the phenotypic behaviour of feline L. infantum strains. This study brings new insights into the hypothetical role of cats as reservoir hosts of L. infantum since the parasites found in them are phenotypically identical to those of dogs and humans. However, further studies on the transmission dynamics should be encouraged to fully establish the status of cats in the maintenance of L. infantum foci.

Diversity and Within-Host Evolution of Leishmania donovani from Visceral Leishmaniasis Patients with and without HIV Coinfection in Northern Ethiopia

Visceral leishmaniasis (VL) is a fatal disease and a growing public health problem in East Africa, where Ethiopia has one of the highest VL burdens. The largest focus of VL in Ethiopia is driven by high prevalence in migrant agricultural workers and associated with a high rate of coinfection with HIV. This coinfection makes VL more difficult to treat successfully and is associated with a high rate of relapse, with VL/HIV patients frequently experiencing many relapses of VL before succumbing to this infection. We present genome-wide data on Leishmania donovani isolates from a longitudinal study of cohorts of VL and VL/HIV patients reporting to a single clinic in Ethiopia. Extensive clinical data allow us to investigate the influence of coinfection and relapse on the populations of parasites infecting these patients. We find that the same parasite population is responsible for both VL and VL/HIV infections and that, in most cases, disease relapse is caused by recrudescence of the population of parasites that caused primary VL. Complex, multiclonal infections are present in both primary and relapse cases, but the infrapopulation of parasites within a patient loses genetic diversity between primary disease presentation and subsequent relapses, presumably due to a population bottleneck induced by treatment. These data suggest that VL/HIV relapses are not caused by genetically distinct parasite infections or by reinfection. Treatment of VL does not lead to sterile cure, and in VL/HIV, the infecting parasites are able to reestablish after clinically successful treatment, leading to repeated relapse of VL. IMPORTANCE Visceral leishmaniasis (VL) is the second largest cause of deaths due to parasite infections and a growing problem in East Africa. In Ethiopia, it is particularly associated with migrant workers moving from regions of nonendemicity for seasonal agricultural work and is frequently found as a coinfection with HIV, which leads to frequent VL relapse following treatment. Insight into the process of relapse in these patients is thus key to controlling the VL epidemic in Ethiopia. We show that there is little genetic differentiation between the parasites infecting HIV-positive and HIV-negative VL patients. Moreover, we provide evidence that relapses are caused by the initially infecting parasite population and that treatment induces a loss of genetic diversity in this population. We propose that restoring functioning immunity and improving antiparasitic treatment may be key in breaking the cycle of relapsing VL in VL/HIV patients.

Anti-Leishmania activity of extracts from Piper cabralanum C.DC. (Piperaceae)

Species of Piperaceae are known by biological properties, including antiparasitic such as leishmanicidal, antimalarial and in the treatment of schistosomiasis. The aim of this work was to evaluate the antileishmania activity, cytotoxic effect, and macrophage activation patterns of the methanol (MeOH), hexane (HEX), dichloromethane (DCM) and ethyl acetate (EtOAc) extract fractions from the leaves of Piper cabralanum C.DC. The MeOH, HEX and DCM fractions inhibited Leishmanina amazonensis promastigote-like forms growth with a half maximal inhibitory concentration (IC₅₀) of 144.54, 59.92, and 64.87 μg/mL, respectively. The EtOAc fraction did not show any relevant activity. The half maximal cytotoxic concentration (CC₅₀) for macrophages were determined as 370.70, 83.99, 113.68 and 607 μg/mL for the MeOH, HEX and DCM fractions, respectively. The macrophage infectivity was concentration-dependent, especially for HEX and DCM. MeOH, HEX and DCM fractions showed activity against L. amazonensis with low cytotoxicity to murine macrophages and lowering infectivity by the parasite. Our results provide support for in vivo studies related to a potential application of P. cabralanum extract and fractions as a promising natural resource in the treatment of leishmaniasis.

Autochthonous and imported tegumentary leishmaniasis in Catalonia (Spain): Aetiological evolution in the last four decades and usefulness of different typing approaches based on biochemical, molecular and proteomic markers

Leishmaniasis is a transmissible disease caused by Leishmania protozoa. Spain is endemic for both visceral and cutaneous leishmaniasis, the autochthonous aetiological agent being Leishmania infantum. Around the world, the L. donovani complex is associated with visceral symptoms, while any species of the Leishmania or Viannia subgenera affecting human can produce tegumentary forms. In a context of growing numbers of imported cases, associated with globalisation, the aim of this study was to analyse the aetiological evolution of human tegumentary leishmaniasis in a region of Spain (Catalonia). Fifty-six Leishmania strains, isolated from 1981 to 2018, were analysed using MLEE, gene sequencing (hsp70, rpoIILS, fh and ITS2) and MALDI-TOF. The utility of these different analytical methods was compared. The results showed an increase in leishmaniasis over the two last decades, particularly imported cases, which represented 39% of all cases studied. Leishmania infantum, L. major, L. tropica, L. braziliensis, L. guyanensis and L. panamensis were identified. The combination of molecular and enzymatic methods allowed the identification of 29 different strain types (A to AC). Strain diversity was higher in L. (Viannia), whilst the different L. major types were relatable with geo-temporal data. Among the autochthonous cases, type C prevailed throughout the studied period (39%). Minor types generally appeared within a short time interval. While all the techniques provided identical identification at the species complex level, MALDI-TOF and rpoIILS or fh sequencing would be the most suitable identification tools for clinical practice, and the tandem hsp70-ITS2 could substitute MLEE in the epidemiological field.