In the Indian subcontinent, infection with Leishmania donovani can cause fatal visceral leishmaniasis. Genetic variation in L. donovani is believed to occur rapidly from environmental changes and through selective drug pressures, thereby allowing continued disease occurrence in this region. All previous molecular markers that are commonly in use multilocus microsatellite typing and multilocus sequence typing, were monomorphic in L. donovani originating from the Indian subcontinent (with only a few exceptions) and hence are not suitable for this region. An multilocus sequence typing scheme consisting of a new set of seven housekeeping genes was developed in this study, based on recent findings from whole genome sequencing data. This new scheme was used to assess the genetic diversity amongst 22 autochthonous L. donovani isolates from Bangladesh. Nineteen additional isolates of the L. donovani complex (including sequences of L. donovani reference strain BPK282A1) from other countries were included for comparison. By using restriction fragment length polymorphism of the internal transcribed spacer 1 region (ITS1-RFLP) and ITS1 sequencing, all Bangladeshi isolates were confirmed to be L. donovani. Population genetic analyses of 41 isolates using the seven new MLST loci clearly separated L. donovani from Leishmania infantum. With this multilocus sequence typing scheme, seven genotypes were identified amongst Bangladeshi L. donovani isolates, and these isolates were found to be phylogenetically different compared with those from India, Nepal, Iraq and Africa. This novel multilocus sequence typing approach can detect intra- and inter-species variations within the L. donovani complex, but most importantly these molecular markers can be applied to resolve the phylogenetically very homogeneous L. donovani strains from the Indian subcontinent. Four of these markers were found suitable to differentiate strains originating from Bangladesh, with marker A2P being the most discriminative one.