Skip directly to site content Skip directly to page options Skip directly to A-Z link Skip directly to A-Z link Skip directly to A-Z link
Volume 26, Number 8—August 2020

Leishmania donovani Infection with Atypical Cutaneous Manifestations, Himachal Pradesh, India, 2014–2018

Lovlesh Thakur, Kiran K. Singh, Hemant R. Kushwaha, Sudarshan K. Sharma, Vinay Shankar, Ajeet Negi, Ghanshyam Verma, Sandhya Kumari, Aklank Jain, and Manju JainComments to Author 
Author affiliations: Central University of Punjab, Bathinda, India (L. Thakur, K.K. Singh, A. Jain, M. Jain); Jawaharlal Nehru University, New Delhi, India (H.R. Kushwaha); Indira Gandhi Medical College, Shimla, India (S.K. Sharma, A. Negi, G. Verma, S. Kumari); Maharishi Markandeshwar Medical College and Hospital, Kumarhatti-Solan, India (V. Shankar)

Cite This Article


We conducted a molecular study of parasite sequences from a cohort of cutaneous leishmaniasis patients in Himachal Pradesh, India. Results revealed atypical cutaneous disease caused by Leishmania donovani parasites. L. donovani variants causing cutaneous manifestations in this region are different from those causing visceral leishmaniasis in northeastern India.

Leishmaniasis is a complex disease with cutaneous, mucocutaneous, or visceral manifestations depending on the parasite species and host immunity. Despite continued elimination efforts, leishmaniasis continues to afflict known and newer endemic regions, where 0.5–0.9 million new cases of visceral leishmaniasis (VL) and 0.6–1.0 million new cases of cutaneous leishmaniasis (CL) occur every year (1). An increase in VL and CL cases from newer foci and atypical disease manifestation pose a challenge to leishmaniasis control programs (27). Unlike the known species-specific disease phenotype, parasite variants can cause atypical disease, so that Leishmania species generally associated with VL can cause CL and vice versa.

In India, VL caused by L. donovani parasites in the northeastern region and CL caused by L. tropica in the western Thar Desert represent the prevalent forms of the disease (2). Himachal Pradesh is a more recently leishmaniasis-endemic state in northwest where VL and CL coexist; CL incidence is higher than VL incidence and most cases are attributable to L. donovani instead of L. tropica infection (8,9). Sharma et al. conducted limited molecular analysis of a few CL cases and reported preliminary findings (8). For an in-depth study on the involvement of L. donovani parasites in CL cases, we conducted a comprehensive molecular analysis of CL cases in Himachal Pradesh.

The Study

During 2014–2018, an increase in CL cases occurred in Himachal Pradesh; case reports came from different tehsils (i.e., townships) in Kinnaur, Shimla, and Kullu and the previously nonendemic districts of Mandi and Solan (Appendix Table 1, Figure 1). We confirmed 60 CL cases indigenous to the state with detailed patient information, demonstration of the presence of Leishman-Donovan bodies and CL-specific histopathologic changes in skin lesional specimens, and PCR detection of parasitic infection (Appendix).

We conducted PCR and restriction fragment-length polymorphism (RFLP) analysis of parasite species–specific internal transcribed spacer 1 (ITS1) sequences by using appropriate standard controls. We detected the expected ≈320-bp product with a HaeIII RFLP pattern specific to L. donovani complex in all patient biopsy specimens, indicating L. donovani, L. infantum, or both as the causative agent of infection (Appendix Figure 4) (10).

Figure 1

Thumbnail of ITS1-based molecular analysis of clinical isolates from cutaneous leishmaniasis (CL) patients, Himachal Pradesh, India, 2014–2018. A) Multiple sequence alignment of ITS1 microsatellite repeat sequences of representative parasite isolates from CL patients with those of L. donovani complex reference strains from different geographic regions. Sequences were aligned by using BioEdit sequence alignment program ( B) Phylogenetic tree of ITS1 sequ

Figure 1. ITS1-based molecular analysis of clinical isolates from cutaneous leishmaniasis (CL) patients, Himachal Pradesh, India, 2014–2018. A) Multiple sequence alignment of ITS1 microsatellite repeat sequences of representative parasite isolates from CL patients...

BLAST analysis ( of 44 ITS1 test sequences showed all the samples to be closest to L. donovani, having maximum identity to L. donovani isolates from Bhutan (GenBank accession nos. JQ730001–2) and possibly L. infantum. None of the CL cases were consistent with L. tropica infection, unlike in a previous report (8). To distinguish whether HP isolates were L. donovani, L. infantum, or both and to infer genetic and geographic relatedness between these isolates and standard reference strains, we performed ITS1 microsatellite repeat analysis and phylogenetic classification (1113). The 4 ITS1 polymorphic microsatellite repeat analysis indicate HP isolates different from L. infantum and closest to the L. donovani isolates from Bhutan (Table 1; Figure 1, panel A). We detected a polymorphism in the third poly (TA) microsatellite locus with 5 repeats and an atypical insert of TAA and the fourth poly (A) microsatellite tract with 8 repeats; these polymorphisms were identical to the VL-causing L. donovani isolates from Bhutan. An L. donovani Chandigarh isolate originally from HP is reported to be closest to the Bhutan isolates and matched with HP isolates at the third poly (TA) stretch (12). However, Himachal Pradesh isolates were distinct at the first poly C and the second poly A microsatellite tracts and had heterogeneous base sequences. Thus, these isolates represent L. donovani genetic variants; none showed the ITS1 sequence type previously assigned to the referred L. donovani isolates by Kuhls et al. (13). Our phylogenetic analysis of 44 ITS1 test sequences and ITS1 reference sequences placed all the CL-causing L. donovani isolates from Himachal Pradesh into a discrete cluster different from the VL-causing L. donovani from India and elsewhere and the CL-causing L. donovani isolates from Sri Lanka. The Himachal Pradesh CL isolates within the cluster exhibited considerable heterogeneity (Table 1; Figure 1, panel B; Appendix Table 4).

Figure 2

Thumbnail of 6PGDH-based molecular analysis of clinical isolates from cutaneous leishmaniasis (CL) patients, Himachal Pradesh, India, 2014–2018. A) Sequence alignment of partial 6PGDH amino acid of CL isolates exhibit replacement of asparagine (N) with aspartic acid (D) at position 326 analogous to visceral leishmaniasis–causing and CL-causing isolates from Sri Lanka. B) Phylogenetic tree for 6PGDH sequences from CL test isolates (designated as HPCL, numbered in order of their collection) and st

Figure 2. 6PGDH-based molecular analysis of clinical isolates from cutaneous leishmaniasis (CL) patients, Himachal Pradesh, India, 2014–2018. A) Sequence alignment of partial 6PGDH amino acid of CL isolates exhibit replacement of asparagine (N)...

Sequences of the 6-phosphogluconate dehydrogenase gene (6PGDH) exhibit a high degree of polymorphism and have been used to identify Leishmania species and differentiate region-specific zymodemes (14). We performed multiple sequence alignment of the representative partial 6PGDH amino acid sequences from Himachal Pradesh isolates by using the homologous 6PGDH protein sequences of the reference Leishmania isolates to determine their genetic and geographic relatedness (Table 2; Figure 2, panel A; Appendix Table 4, Figure 5). Himachal Pradesh isolates exhibited a 6PGDH sequence specific to Mon-37 and different from Mon-2 (having aspartic acid in place of asparagine) at position 326 (Figure 2, panel A). Thus, CL-causing L. donovani from Himachal Pradesh were distinct from the most common VL-causing India Mon-2 L. donovani and the Bangladesh L. donovani isolate, whereas they were similar to the CL-causing L. donovani isolate from Kerala and CL- and VL-causing Mon-37 isolates from Sri Lanka and the isolates from Kenya, Brazil, and China.

Phylogenetic analysis of 6PGDH amino acid sequences of CL isolates grouped them into a heterogeneous cluster; variants were closer to a viserotropic L. donovani isolate from Sri Lanka and distinct from the VL-causing L. donovani isolates from India and Bangladesh and CL-causing isolates from Kerala and Sri Lanka (Figure 2, panel B). However, the HPCL55 isolate (GenBank accession no. MH208450) grouped differently. The HPCL49 isolate (GenBank accession no. MH208446) showed relatedness to the standard L. infantum strain, although ITS1 analysis using BLAST and microsatellite repeat sequences showed regions of similarity with L. donovani. ITS1 and 6PGDH sequence analysis suggest that Himachal Pradesh isolates from CL patients consist of heterogenous L. donovani variants and possibly represent hybrid genotypes.

None of the CL patients had VL-specific symptoms or VL history. Ten of 43 patient blood samples tested positive for rK39 antibody, and 37 of 51 samples were positive for the circulating parasite DNA with L. donovani–specific ITS1 (Appendix Figure 6, panel A, B). The result suggests asymptomatic systemic L. donovani infection in a fraction of CL patients.


The presence of leishmaniasis in Himachal Pradesh is not yet well known in India and globally (15). Our epidemiologic study shows newer CL pockets during 2014–2018; thus, the state needs to be recognized as leishmaniasis-endemic by public health authorities (Appendix Figure 1). We conclude that CL cases in Himachal Pradesh are caused by L. donovani variants distinct from the viscerotropic L. donovani strain from northeast India. The CL isolates in Himachal Pradesh exhibit considerable heterogeneity and indicate the possible existence of genetic hybrids. The scenario appears somewhat similar to Sri Lanka and Kerala, where L. donovani parasites cause cutaneous disease, albeit with differences in the region-specific L. donovani variants. In lieu of the coexistence of VL and CL in Himachal Pradesh, parasite isolates from VL patients also need to be characterized. To understand the biology of atypical L. donovani variants with cutaneous manifestations and to genetically differentiate the dermotropic versus viscerotropic potential of L. donovani variants, comparison of CL- and VL-causing isolates in Himachal Pradesh using whole-genome sequence analysis is required.

L. donovani parasites in the blood of some CL patients represent human reservoirs similar to asymptomatic VL carriers, and the parasite variants have the potential to cause full-blown VL manifestations. An elaborate surveillance program dedicated to the Himachal Pradesh region is urgently required for better diagnosis, treatment, prediction of parasite variants in different afflicted pockets, and prevention of transmission of the disease to other regions.

Ms. Thakur is a PhD student in the Department of Zoology at the Central University of Punjab, Bathinda, Punjab, India. Her primary research interests include the epidemiology and pathogenesis of infectious and parasitic diseases.



We thank Greg Matlashewski for his helpful suggestions and comments to improve the manuscript, Wen Wei Zhang for technical suggestions, and Rentala Madhubala for generously gifting Leishmania donovani and L. major standard cultures.



  1. World Health Organization. Leishmaniasis. 2019 [cited 2019 Sep 1].
  2. Thakur  L, Singh  KK, Shanker  V, Negi  A, Jain  A, Matlashewski  G, et al. Atypical leishmaniasis: A global perspective with emphasis on the Indian subcontinent. PLoS Negl Trop Dis. 2018;12:e0006659. DOIPubMed
  3. Sandhya  R, Rakesh  P, Dev  S. Emergence of visceral leishmaniasis in Kollam District, Kerala, southern India. Int J Community Med Public Health. 2019;6:13502. DOI
  4. Kumar  NP, Srinivasan  R, Anish  TS, Nandakumar  G, Jambulingam  P. Cutaneous leishmaniasis caused by Leishmania donovani in the tribal population of the Agasthyamala Biosphere Reserve forest, Western Ghats, Kerala, India. J Med Microbiol. 2015;64:15763. DOIPubMed
  5. Siriwardana  Y, Zhou  G, Deepachandi  B, Akarawita  J, Wickremarathne  C, Warnasuriya  W, et al. Trends in recently emerged Leishmania donovani induced cutaneous leishmaniasis, Sri Lanka, for the first 13 Years. BioMed Res Int. 2019;2019:4093603. DOIPubMed
  6. Sharma  NL, Mahajan  VK, Negi  AK. Epidemiology of a new focus of localized cutaneous leishmaniasis in Himachal Pradesh. J Commun Dis. 2005;37:2759.PubMed
  7. Kumari  S, Garg  A. Lip leishmaniasis: a new emerging clinical form of cutaneous leishmaniasis from sub-Himalayan Region. Journal of Medical Science and Clinical Research. 2018;06:629.
  8. Sharma  NL, Mahajan  VK, Kanga  A, Sood  A, Katoch  VM, Mauricio  I, et al. Localized cutaneous leishmaniasis due to Leishmania donovani and Leishmania tropica: preliminary findings of the study of 161 new cases from a new endemic focus in himachal pradesh, India. Am J Trop Med Hyg. 2005;72:81924. DOIPubMed
  9. Sharma  NL, Sood  A, Arora  S, Kanga  A, Mahajan  V, Negi  A, et al. Characteristics of Leishmania spp. isolated from a mixed focus of cutaneous and visceral leishmaniasis in Himachal Pradesh (India). Internet J Third World Med. 2009;7(8).
  10. el Tai  NO, Osman  OF, el Fari  M, Presber  W, Schönian  G. Genetic heterogeneity of ribosomal internal transcribed spacer in clinical samples of Leishmania donovani spotted on filter paper as revealed by single-strand conformation polymorphisms and sequencing. Trans R Soc Trop Med Hyg. 2000;94:5759. DOIPubMed
  11. Dávila  AM, Momen  H. Internal-transcribed-spacer (ITS) sequences used to explore phylogenetic relationships within Leishmania. Ann Trop Med Parasitol. 2000;94:6514. DOIPubMed
  12. Yangzom  T, Cruz  I, Bern  C, Argaw  D, den Boer  M, Vélez  ID, et al. Endemic transmission of visceral leishmaniasis in Bhutan. Am J Trop Med Hyg. 2012;87:102837. DOIPubMed
  13. Kuhls  K, Mauricio  IL, Pratlong  F, Presber  W, Schönian  G. Analysis of ribosomal DNA internal transcribed spacer sequences of the Leishmania donovani complex. Microbes Infect. 2005;7:122434. DOIPubMed
  14. Ranasinghe  S, Zhang  W-W, Wickremasinghe  R, Abeygunasekera  P, Chandrasekharan  V, Athauda  S, et al. Leishmania donovani zymodeme MON-37 isolated from an autochthonous visceral leishmaniasis patient in Sri Lanka. Pathog Glob Health. 2012;106:4214. DOIPubMed
  15. World Health Organization. Leishmaniasis country profile–2015, India. 2017 [cited 2019 Sep 1].




Cite This Article

DOI: 10.3201/eid2608.191761

Original Publication Date: July 03, 2020

Table of Contents – Volume 26, Number 8—August 2020


Please use the form below to submit correspondence to the authors or contact them at the following address:

Manju Jain, Department of Biochemistry, Central University of Punjab, City Campus, Mansa Rd, Bathinda, Punjab, 151001, India,

Send To

character(s) remaining.

Comment submitted successfully, thank you for your feedback.


Page created: June 10, 2020
Page updated: July 18, 2020
Page reviewed: July 18, 2020
The conclusions, findings, and opinions expressed by authors contributing to this journal do not necessarily reflect the official position of the U.S. Department of Health and Human Services, the Public Health Service, the Centers for Disease Control and Prevention, or the authors' affiliated institutions. Use of trade names is for identification only and does not imply endorsement by any of the groups named above.