Leishmaniasis in Eurasia and Africa: geographical distribution of vector species and pathogens

Analysis of the current risk of Leishmania infantum transmission for domestic dogs in Spain and Portugal and its future projection in climate change scenarios

Canine leishmaniosis, caused by the protozoan parasite Leishmania infantum, is a cosmopolitan vector-borne zoonosis, transmitted principally by Phlebotomus perniciosus in Spain and Portugal, where it is considered an endemic disease. Ecoinformatics tools such as ecological niche models (ENM) have been successfully tested to model the distribution of the risk of infection of different parasitosis as they take into account environmental variables vital for their survival. The risk map proposed in this study combines the potential distribution of Ph. perniciosus in the Iberian Peninsula and the calculation of the infection rate of the parasite in the vector to model the risk of contracting the disease in a more realistic way. In fact, this weighting strategy improves the predictive power of the resulting model (R2 = 0.42, p = < 0.01) compared to the Ph. perniciosus ENM model alone (R2 = 0.13, p > 0.05). The places with the highest risk of transmission are the southwest and central peninsular area, as well as the Mediterranean coast, the Balearic Islands and the Ebro basin, places where the ideal habitat of Ph. perniciosus and the infection rate is also high. In the case of future projections under climate change scenarios, an increase in the risk of infection by L. infantum can be observed in most of the territory (4.5% in 2040, 71.6% in 2060 and 63% in 2080), mainly in the northern part of the peninsula. The use of ENMs and their weighting with the infection rate in Ph. perniciosus is a useful tool in predicting the risk of infection for L. infantum in dogs for a given area. In this way, a more complete model can be obtained to facilitate prevention and control.

A Rare Case of Imported Cutaneous Leishmaniasis Caused by Leishmania infantum in the Republic of Korea, 2021

Leishmaniasis is a neglected tropical disease and an infectious disease transmitted by sandflies that occurs worldwide. In the absence of physicians seeking to identify the causes of disease in non-endemic areas, appropriate diagnoses cannot be made, thereby hampering effective treatment. In this report, we examined a nodular lesion on a patient's chin by performing a biopsy and molecular analysis. The biopsy finding led to the identification of a Leishmania amastigote. On the basis of PCR analysis of the internal transcribed spacer 1 gene and 5.8 S ribosomal RNA with a subsequent BLAST search, we identified the causal organism as Leishmania infantum. The patient, who had visited Spain from 1 July to 31 August 2018, was accordingly diagnosed with cutaneous leishmaniasis and was administered liposomal amphotericin B, which successfully treated the skin lesion. Travel history plays an important role in the diagnosis of leishmaniasis, and physicians should bear in mind that travelers can also introduce diseases and pathogens to non-endemic areas. Identification of Leishmania at the species level will increase the efficacy of treatment.

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.

Predicting Plasmodium knowlesi transmission risk across Peninsular Malaysia using machine learning-based ecological niche modeling approaches

The emergence of potentially life-threatening zoonotic malaria caused by Plasmodium knowlesi nearly two decades ago has continued to challenge Malaysia healthcare. With a total of 376 P. knowlesi infections notified in 2008, the number increased to 2,609 cases in 2020 nationwide. Numerous studies have been conducted in Malaysian Borneo to determine the association between environmental factors and knowlesi malaria transmission. However, there is still a lack of understanding of the environmental influence on knowlesi malaria transmission in Peninsular Malaysia. Therefore, our study aimed to investigate the ecological distribution of human P. knowlesi malaria in relation to environmental factors in Peninsular Malaysia. A total of 2,873 records of human P. knowlesi infections in Peninsular Malaysia from 1st January 2011 to 31st December 2019 were collated from the Ministry of Health Malaysia and geolocated. Three machine learning-based models, maximum entropy (MaxEnt), extreme gradient boosting (XGBoost), and ensemble modeling approach, were applied to predict the spatial variation of P. knowlesi disease risk. Multiple environmental parameters including climate factors, landscape characteristics, and anthropogenic factors were included as predictors in both predictive models. Subsequently, an ensemble model was developed based on the output of both MaxEnt and XGBoost. Comparison between models indicated that the XGBoost has higher performance as compared to MaxEnt and ensemble model, with AUC_ROC values of 0.933 ± 0.002 and 0.854 ± 0.007 for train and test datasets, respectively. Key environmental covariates affecting human P. knowlesi occurrence were distance to the coastline, elevation, tree cover, annual precipitation, tree loss, and distance to the forest. Our models indicated that the disease risk areas were mainly distributed in low elevation (75-345 m above mean sea level) areas along the Titiwangsa mountain range and inland central-northern region of Peninsular Malaysia. The high-resolution risk map of human knowlesi malaria constructed in this study can be further utilized for multi-pronged interventions targeting community at-risk, macaque populations, and mosquito vectors.

High genome plasticity and frequent genetic exchange in Leishmania tropica isolates from Afghanistan, Iran and Syria

Background

The kinetoplastid protozoan Leishmania tropica mainly causes cutaneous leishmaniasis in humans in the Middle East, and relapse or treatment failure after treatment are common in this area. L. tropica’s digenic life cycle includes distinct stages in the vector sandfly and the mammalian host. Sexual reproduction and genetic exchange appear to occur more frequently than in other Leishmania species. Understanding these processes is complicated by chromosome instability during cell division that yields aneuploidy, recombination and heterozygosity. This combination of rare recombination and aneuploid permits may reveal signs of hypothetical parasexual mating, where diploid cells fuse to form a transient tetraploid that undergoes chromosomal recombination and gradual chromosomal loss.

Methodology/principal findings

The genome-wide SNP diversity from 22 L. tropica isolates showed chromosome-specific runs of patchy heterozygosity and extensive chromosome copy number variation. All these isolates were collected during 2007–2017 in Sweden from patients infected in the Middle East and included isolates from a patient possessing two genetically distinct leishmaniasis infections three years apart with no evidence of re-infection. We found differing ancestries on the same chromosome (chr36) across multiple samples: matching the reference genome with few derived alleles, followed by blocks of heterozygous SNPs, and then by clusters of homozygous SNPs with specific recombination breakpoints at an inferred origin of replication. Other chromosomes had similar marked changes in heterozygosity at strand-switch regions separating polycistronic transcriptional units.

Conclusion/significance

These large-scale intra- and inter-chromosomal changes in diversity driven by recombination and aneuploidy suggest multiple mechanisms of cell reproduction and diversification in L. tropica, including mitotic, meiotic and parasexual processes. It underpins the need for more genomic surveillance of Leishmania, to detect emerging hybrids that could spread more widely and to better understand the association between genetic variation and treatment outcome. Furthering our understanding of Leishmania genome evolution and ancestry will aid better diagnostics and treatment for cutaneous leishmaniasis caused by L.tropica in the Middle East.

Cutaneous leishmaniasis is mainly caused by Leishmania tropica in the Middle East, where it is known for treatment failure and a need for prolonged and/or multiple treatments. Several factors affect the clinical presentation and treatment outcome, such as host genetic variability and specific immune response, as well as environmental factors and the vector species. Little is known about the parasite genome and its influence on treatment response. By analysing the genome of 22 isolates of L. tropica, we have revealed extensive genomic variation and a complex population structure with evidence of genetic exchange within and among the isolates, indicating a possible presence of sexual or parasexual mechanisms. Understanding the Leishmania genome better may improve future treatment and better understanding of treatment failure and relapse.

Canine leishmaniasis in Balkan - A review of occurrence and epidemiology

Canine leishmaniasis is a parasitic zoonotic infection of dogs caused by the Leishmania spp. parasites. In Europe, canine leishmaniasis is among the most important vector-borne parasitic diseases of dogs, with majority of data originating from countries of southwestern Europe. With its scarce and outdated records, the Balkan region presents the "blank space" in the European map of canine leishmaniasis knowledge. Numerous records obtained during the past decade and a half indicate the emergence of autochthonous canine leishmaniasis in new settings all across Balkan, a significant increase in the numbers of infected dogs in previously endemic areas, and northward spread of the disease. As the majority of canine leishmaniasis data from Balkan countries are several decades old and generally written in native Balkan languages, hence not available to the broader area of scientists, comprehensive search of scientific literature published from 1934 to January 2021 was conducted in order to provide a better perspective to the past and current canine leishmaniasis situation in Balkan. This review provides exhaustive information on all leishmaniasis records in dogs, wild animals and cats, available from 10 countries of the Balkan region. It identifies dogs as the main disease reservoirs and highlights the importance of wild animals, and potentially cats, as disease reservoirs, which is of great public and veterinary health concern. Additionally, certain shortcomings in canine leishmaniasis research, management, surveillance and control were identified.

Evolution, systematics and historical biogeography of sand flies of the subgenus Paraphlebotomus (Diptera, Psychodidae, Phlebotomus) inferred using restriction-site associated DNA markers

Phlebotomine sand flies are the main natural vectors of Leishmania, which cause visceral and tegumentary tropical diseases worldwide. However, their taxonomy and evolutionary history remain poorly studied. Indeed, as for many human disease vectors, their small size is a challenge for morphological and molecular works. Here, we successfully amplified unbiased copies of whole genome to sequence thousands of restriction-site associated DNA (RAD) markers from single specimens of phlebotomines. RAD markers were used to infer a fully resolved phylogeny of the subgenus Paraphlebotomus (11 species + 5 outgroups, 32 specimens). The subgenus was not recovered as monophyletic and we describe a new subgenus Artemievus subg. nov. Depaquit for Phlebotomus alexandri. We also confirm the validity of Ph. riouxi which is reinstated as valid species. Our analyses suggest that Paraphlebotomus sensu nov. originated ca 12.9-8.5 Ma and was possibly largely distributed from peri-Mediterranean to Irano-Turanian regions. Its biogeographical history can be summarized into three phases: i) a first split between Ph. riouxi + Ph. chabaudi and other species that may have resulted from the rise of the Saharan belt ca 8.5 Ma; ii) a Messinian vicariant event (7.3-5.3 Ma) during which the prolonged drought could have resulted in the divergence of main lineages; iii) a recent radiation event (3-2 Ma) that correspond to cycles of wet and dry periods in the Middle East and the East African subregions during the Pleistocene. Interestingly these cycles are also hypothetical drivers of the diversification of rodents, in the burrows of which Paraphlebotomus larvae develop. By meeting the challenge of sequencing pangenomics markers from single, minute phlebotomines, this work opens new avenues for improving our understanding of the epidemiology of leishmaniases and possibly other human diseases transmitted by arthropod vectors.

Genetic diversity of Leishmania tropica: Unexpectedly complex distribution pattern

In this study, we characterized a collection of clinical samples obtained from Syrian and Turkish patients with cutaneous leishmaniasis using internal transcribed spacer 1 (ITS1) sequences. All obtained sequences belonged to Leishmania tropica. Combining them with those available from GenBank allowed us performing a broad-scale analysis of genetic diversity for this species. We demonstrated that L. tropica has a complex phylogeographic pattern with some haplotypes being widespread across endemic countries and others restricted to particular regions. We hypothesize that at least some of them may be associated with alternative vectors or animal reservoirs.

Expanded genome-wide comparisons give novel insights into population structure and genetic heterogeneity of Leishmania tropica complex

Leishmania tropica is one of the main causative agents of cutaneous leishmaniasis (CL). Population structures of L. tropica appear to be genetically highly diverse. However, the relationship between L. tropica strains genomic diversity, protein coding gene evolution and biogeography are still poorly understood. In this study, we sequenced the genomes of three new clinical L. tropica isolates, two derived from a recent outbreak of CL in camps hosting Syrian refugees in Lebanon and one historical isolate from Azerbaijan to further refine comparative genome analyses. In silico multilocus microsatellite typing (MLMT) was performed to integrate the current diversity of genome sequence data in the wider available MLMT genetic population framework. Single nucleotide polymorphism (SNPs), gene copy number variations (CNVs) and chromosome ploidy were investigated across the available 18 L. tropica genomes with a main focus on protein coding genes. MLMT divided the strains in three populations that broadly correlated with their geographical distribution but not populations defined by SNPs. Unique SNPs profiles divided the 18 strains into five populations based on principal component analysis. Gene ontology enrichment analysis of the protein coding genes with population specific SNPs profiles revealed various biological processes, including iron acquisition, sterols synthesis and drug resistance. This study further highlights the complex links between L. tropica important genomic heterogeneity and the parasite broad geographic distribution. Unique sequence features in protein coding genes identified in distinct populations reveal potential novel markers that could be exploited for the development of more accurate typing schemes to further improve our knowledge of the evolution and epidemiology of the parasite as well as highlighting protein variants of potential functional importance underlying L. tropica specific biology.