L eishmaniasis 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 (2-7). 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, https://wwwnc.cdc.gov/EID/article/26/8/19-1761-App1.pdf). 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 fragmentlength 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). BLAST analysis (https://blast.ncbi.nlm.nih. gov/Blast.cgi) 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