Phylogenetic position of Leishmania tropica isolates from an old endemic focus in south-eastern Iran; relying on atypical cutaneous leishmaniasis

A Geomedical Survey: Is There an Association Between Climatic Conditions and Leishmania Species Distribution in Iran During the Years 1999-2021?

PURPOSE

Iran is among the high-risk leishmaniasis regions in the world. WHO recommends the use of GIS as an ideal tool for healthcare authorities to predict the evolution of a disease, delimit the risk of outbreaks and identify critical areas. The aim of this research is to find the association between the main species of Leishmania (L. major, L. tropica, L. infantum) dispersion and climatic variables in Iran.

METHODS

All molecular-based reports of leishmaniasis from Iran between 1999 and 2021 were gathered from reliable medical sources. Meteorological data (air and soil temperatures, annual rainfall and humidity) of the country along the study period were obtained from the Iranian Climatological Research Centre. The data concerning species distribution and climatic conditions during this period were moved to a base-map through raster layers using ArcGIS 10.4.1 software. The relationship between parasitological and climatic models was examined using ANOVA.

RESULTS

High risk area maps, based on the cut-off thresholds, were generated for Leishmania major, L. tropica and L. infantum. According to the molecular-based reports, the L. major distribution was significantly related to all climatic variables, while L. tropica was merely related to rainfall and humidity, and the L. infantum distribution was significantly associated with rainfall, soil and air temperatures.

CONCLUSION

The association between climatic conditions and Leishmania species distribution in Iran has been confirmed. Consequently, both, the relationship between climatic conditions and the geographical distribution of Leishmania species, and the use of GIS to better understand the spatial epidemiology of leishmaniasis, have been reaffirmed.

Atypical manifestations of cutaneous leishmaniasis in a boy from Afghanistan

Cutaneous leishmaniasis is one of the most common forms of leishmaniasis in Afghanistan, which is caused by different Leishmania species that are transmitted to humans by sandflies. Its clinical manifestations are very diverse and depend on a variety of parasite and host factors.

Molecular characteristic of treatment failure clinical isolates of Leishmania major

Background

Leishmaniasis is a prevalent tropical disease caused by more than 20 Leishmania species (Protozoa, Kinetoplastida and Trypanosomatidae). Among different clinical forms of the disease, cutaneous leishmaniasis is the most common form, with an annual 0.6–1 million new cases reported worldwide. This disease’s standard treatment is pentavalent antimonial (Sb^V) that have been used successfully since the first half of the 20th century as a first-line drug. However, treatment failure is an increasing problem that is persistently reported from endemic areas. It is important to define and standardize tests for drug resistance in cutaneous leishmaniasis. Sb^V must be reduced to its trivalent active form (Sb^III). This reduction occurs within the host macrophage, and the resultant Sb^IIIenters amastigotes via the aquaglyceroporin1 (AQP1) membrane carrier. Overexpression of AQP1 results in hypersensitivity of the parasites to Sb^III, but resistant phenotypes accompany reduced expression, inactivation mutations, or deletion of AQP1. Hence, in this study, a phylogenetic analysis using barcode gene COXII and kDNA minicircle and expression analysis of AQP1 were performed in treatment failure isolates to assess the isolates’ molecular characteristics and to verify possible association with drug response.

Methods

Samples in this study were collected from patients with cutaneous leishmaniasis referred to the Diagnosis Laboratory Center in Isfahan Province, Iran, from October 2017 to December 2019. Among them, five isolates (code numbers 1–5) were categorized as treatment failures. The PCR amplification of barcode gene COXII and kDNA minicircle were done and subsequently analyzed using MEGA (10.0.5) to perform phylogenetics analysis of Treatment failures (TF) and Treatment response (TR) samples. Relative quantification of the AQP1 gene expression of TF and TR samples was assessed by real-time PCR.

Results

All samples were classified as L. major. No amplification failure was observed in the cases of barcode gene COXII and kDNA minicircle amplification. Having excluded the sequences with complete homology using maximum parsimony with the Bootstrap 500 method, four major groups were detected to perform phylogenetic analysis using COXII. The phylogenetic analysis using the barcode target of minicircle showed that all five treatment failure isolates were grouped in a separate sub-clade.

Conclusions

We concluded that the barcode gene COXII and the minicircle kDNA were suitable for identification, differentiation and phylogenetic analysis in treatment failure clinical isolates of Leishmania major. Also, AQP1 gene expression analyses showed that treatment failure isolates had less expression than TR isolates. The isolate with TF and overexpression of the AQP1 gene of other molecular mechanisms such as overexpression of ATP-binding cassette may be involved in the TR, such as overexpression of ATP-binding cassette which requires further research.