J Arthropod-Borne Dis, March 2016, 10(1): 55–64 S Mohammadi and P Parvizi: Simultaneous Morphological … 55 Original Article Simultaneous Morphological and Molecular Characterization of Tatera indica in Southwestern Iran Somayeh Mohammadi 1, 2, *Parviz Parvizi 1 1Molecular Systematics Laboratory, Parasitology Department, Pasteur Institute of Iran, Tehran, Iran 2Islamic Azad University, Tehran Medical Sciences Branch, Tehran, Iran (Received 4 Aug 2014; accepted 10 Dec 2014) Abstract Background: Interest in Tatera indica rodent arises mostly because it is believed that this species is survived among four subspecies reported from Iran, two of which exist in Khuzestan Province. In addition, it might has a role as res- ervoir hosts of zoonotic cutaneous leishmaniasis in the transmission of Leishmania major in some of the widespread Asian foci including southwestern Iran. Methods: Diagnostic morphological and molecular markers for T. indica were sought by characterizing from indi- vidual specimens, such as some taxonomic features and mitochondrial cytochrome b gene that had previously proven useful for the taxonomy of rodents. Wild rodents were caught using live wooden and wire traps. The specimens were identified morphologically using external criteria and molecularly by sequencing of Cyt b gene and phylogenetic analyses. Results: Forty one T. indica were collected and identified morphologically in Khuzestan Province, Iran. Two morphotypes of T. indica were found and classified but sequencing and phylogenetic analyses of mitochondrial Cyt b gene did not support any subspecies between two morphotypes of T. indica. Because all 21 sequences of both morphotypes of T. indica had no variation with only one common and novel haplotype (GenBank accession No KP001566). Conclusion: This is the first time that T. indica was characterized molecularly in Iran. There is no molecular evi- dence for T. indica morphotypes or subspecies, and so a population genetics approach using several polymorphic genes might be employed using species-specific molecular markers. In addition, more specimens of T. indica species in large geographical locations should be tested. Keywords: Tatera indica, Molecular characterization, Reservoir host, Zoonotic cutaneous leishmaniasis, Iran Introduction The Indian Gerbil, Tatera indica (Hard- wicke) is believed to be the main reservoir host of zoonotic cutaneous leishmaniasis (ZCL) in Khuzestan Province, in border of Iran and Iraq, southwest of Iran. Khuzestan Province is an important focus of ZCL in Iran and the reports shows that the occurrence of the disease has increased through the recent years (Javadian et al. 1998, Kia et al. 2010, Shirzadi 2012, Vazirianzadeh et al. 2013). The female Phlebotomus papatasi (Scopoli) (Diptera: Psychodidae) sandfly is the proven vector of the parasitic protozoan L. major (Yakimoff and Schokhor) (Kinetoplastida: Trypanosomatidae), the causative agent of ZCL in Iran (Parvizi and Ready 2008, Bordbar and Parvizi 2014). Tatera indica is distributed through Iran, Afghanistan, Pakistan, India, Sri Lanka, Tur- key, Iraq and Syria (Harrison and Bates 1991, Agrawal 2000, Yigit et al. 2001). The genus Tatera Lataste, (Rodentia: Gerbillinae) are widespread in the Africa, Near East, Middle East and Pakistan sub- continent (Colangelo et al. 2005). At the morphological taxonomic level, 12 species of genus Tatera were identified of which 11 are present in Africa (T. afra, T. *Corresponding author: Dr Parviz Parvizi, E-mail: parp@Pasteur.ac.ir http://jad.tums.ac.ir Published Online: June 27, 2015 J Arthropod-Borne Dis, March 2016, 10(1): 55–64 S Mohammadi and P Parvizi: Simultaneous Morphological … 56 brantsii, T. bohemi, T. guineae, T. inclusa, T. vicina, T. kempi, T. leucogaster, T. nigricauda, T. phillipsi and T. robusta) and only one species (T. indica) lives in Asia (Corbet 1978, Corbet and Hill 1991). Based on morphological characters, T. indica can be considered as Tatera sensu stricto and consequently all 11 African Tatera species were placed in one genus of Ger- billiscus Thomas. Molecular characterizations of genomic DNA data showed that the genus of Tatera is a polyphyletic taxon (Chevert and Dobigny 2005). Of several subspecies of T. indica, four subspecies are present in Iran. T. indica persica exists in Sistan, T. indica scansa in Kerman; and finally T. indica monticola and T. indica bailwardi in Khuzestan (Ellerman 1948, Missone 1959, Mirshamsi et al. 2007). Studies on Tatera genus have highlighted a complex situation by karyological studies in Africa (Colangelo and Civitelli 2001). The main objective of this research was to characterize the T. indica populations based on taxonomic knowledge by using and comparing morphological, genetic and mo- lecular criteria concurrently from a specific geographical location of Iran where this ro- dent is predominant. In addition, phyloge- netic analyses were employed to get the better understanding of T. indica subspecies in case of any homology at their nucleotide scales. Accurate and firm identification of animal reservoir hosts of ZCL are essential, requires any epidemiological research, and are crucial for predicting species-specific popu- lation, and developing control strategies. Materials and Methods Study sites, collections and identification of rodents The active colonies of rodents were iden- tified firstly from ten locations of Khuzestan Province, southern Iran with an altitude of 18 meters above sea level and geographical coordinates as 31° 32′ 73″ N and 48° 69′ 40″ E, and then the rodents were trapped alive in various parts of these areas using wooden and wiry live traps (Fig. 1). Sampling was conducted from 10 cities in Khuzestan Province. However, only from 6 of them (Table. 1) including 15 villages (Fig. 1) T. indica rodents were presented and captured. The collection of samples was done from the colonies of rodents’ burrows located around the villages where ZCL were endemic using 50 live wooden and wirytraps for each location. Sampling was conducted from 10 cities in Khuzestan Province, however in only 6 of them T. indica rodents were present and captured. The collection of samples were done from the colonies of rodents burrows located around the villages where ZCL were en- demic using 50 live wooden and wire traps for each collection. Rodents were caught by baiting cucumbers and dates in 2012–2013. The traps were set up early in the morning and evening then the trapped rodents were transferred to Pasteur Institute of Iran, Teh- ran, and maintained for morphological and molecular testing (Mehrabani 2011). The mor- phological identification of collected speci- mens were employed using external criteria such as color, head and body length, body measurements, ears length, tail length, hind foot length, cranial and dental measurements and cranium. The standard external charac- ters of body and tail were measured using a ruler and/or digital calipers. The sizes of den- tal characters were measured using Nikon measuring microscope MM-40 (Table 1, 2) (Etemad 1978). DNA extraction and PCR amplification of Cyt b gene from rodents All the captured rodents first identified morphologically then genomic DNA of T. indica was extracted using Genet Bio kit, http://jad.tums.ac.ir Published Online: June 27, 2015 J Arthropod-Borne Dis, March 2016, 10(1): 55–64 S Mohammadi and P Parvizi: Simultaneous Morphological … 57 Phenol-Chloroform and ISH-Horovize meth- ods (Parvizi and Ready 2008). After compar- ing these three DNA extraction methods, ISH-Horovize method was found to be more sensitive than others were and selected as main DNA extractions method with minor modifications. Mitochondrial Cytochrome b (Cyt b) gene, which is a maternally inherited marker were used to characterize this rodent. For molecu- lar characterization of rodents, extracted DNA from the tissue and blood of each rodent was used to amplify a fragment of Cyt b with the forward primer UNFOR403 (5'-TGAGGACAAATATCATTCTGAGG- 3') and reveres primer UNREV1025 (5'- GGTTGTCCTCCAATTCATGTTA-3') from single rodent by PCR (Kent and Norris 2005). The PCR products were directly se- quenced, and the sequences edited and aligned using SequencherTM v. 4.4 to permit the identification of haplotypes (= unique sequences) and their phylogenetic analysis using MEGA5.05 and PAUP* softwares (Tamura et al. 2011). Results Collection and Morphological identification of Tatera indica Of all rodents caught during this study a total of 41 T. indica were collected and identified morphologically from 6 different locations in Khuzestan Province (Fig. 1). The most T. indica were captured from Behbahan district with 24 out of 41, then Ahvaz (9/41), Dezful (4/41), Shush (2/41), Shushtar and Ramhormoz each with one samples. Morphological characterizations of T. indica showed that the rodents were in two color groups, buff black (21/41) and buff brown (20/41), total rodent length was about 33–39 cm, body length 16–20 cm, ears length 2.1–2.5 cm, tail length 16–19 cm and hind foot length 3.8–4.5 cm (Table 1) (Etemad 1978). Three morphological characters of total rodent length, hind foot (HF), head and body (HB), ear length (EL) and tail length (TL) were considered for separating T. indica from other rodent species in the region. The details of morphological characters of cap- tured T. indica in Khuzestan are summarized and shown in Table 1 and 2 (Fig. 2). Based on the development degree of the supraorbital crests and dental wear especially upper molars, three ages group among our T. indica samples were classified (Etemad 1978). The incisor teeth of T. indica are formed with a groove and upper molars (Fig. 2b). By considering the colors of rodent body, two types of morphological characters of T. indica were found. The buff- black group was larger with wider snout and they had a semi wide white halo above their eyes. However, the rodents of buff- brown group were small- er, with longer snout and a narrow white halo above their eyes (Table 2, Fig. 3). On bone formation of T. indica, zygomatic bones and pulled Bsyarh is longer than other rodents (Fig. 2a). Molecular characterization and identifica- tion of Tatera indica Cyt b gene of extracted DNA from all 41 T. indica was screened and amplified by PCR. The PCR products of 21 out of 41 were sequenced, aligned and analyzed via molec- ular softwares (SequencherTM v. 4.4, MEGA 5.05). This is the first time that a molecular study is performed on Iranian T. indica, which provided us with novel findings that all the sequences are similar and no variation was found among the nucleotides; therefore, all the samples from one haplotype of T. indica. Cyt b sequences were identified among 21 sequences obtained from individual T. indica. Only haplotype KHT01 (GenBank accession no. KP001566) predominated and was common in Khuzestan Province. After comparing our sequences with only one sub- mitted sequence in GeneBank from France (GenBank accession no: AJ430563.1) some http://jad.tums.ac.ir Published Online: June 27, 2015 J Arthropod-Borne Dis, March 2016, 10(1): 55–64 S Mohammadi and P Parvizi: Simultaneous Morphological … 58 similarities and few differences were observed. Based on the results of our sequences along with the sequences of different rodents submitted in GeneBank, a phylogenic tree was constructed (Fig. 4). Fig. 1. Location of villages and cities in Khuzestan Province, Iran where Tatera indica was sampled Fig. 2. Tatera indica Skull: (a): Zygomatic plate, (b): Upper molars and Upper incisors http://jad.tums.ac.ir Published Online: June 27, 2015 J Arthropod-Borne Dis, March 2016, 10(1): 55–64 S Mohammadi and P Parvizi: Simultaneous Morphological … 59 Table 1. Morphological characters of Tatera indica which were collected in Khuzestan Province, Iran (TB: Total rodent length, HB: Head and body, TL: Tail length, HF: Hind foot, EL: Ear length) T ot al ColorELHFTLHBTB Morphometric char- acters B ro w n B la ck 2.1-2.53.8-4.516-1916-2033-39 Size (Cm) 2. 5 2. 4 2. 3 2. 2 2. 1 4. 5 4. 3 4. 243. 9 3. 8 1918 .5 1817 .5 1716 .5 162019 .5 1817 .5 163938 .5 3837 .5 36 .5 3534 .5 33 .5 33 11001000000010100000010000100000000Ramhormoz L oc at io n s an d th e N u m b er s of T . i n di ca 11010000000010001000010000000000100Shushtar 21120000002000001000010000002000000Shush 42201120000121021001010300001020100 Dezful 241311104532001410907490501301200004970400Behbahan 92722014111411030122121411011111211 Ahvaz 41202115866611491511112710281191191111810101811Total Fig. 3. Buff brown (1) and buff black (2) Tatera indica rodents (b) along with their skull image (a) casing the morphologic differences of two groups http://jad.tums.ac.ir Published Online: June 27, 2015 J Arthropod-Borne Dis, March 2016, 10(1): 55–64 S Mohammadi and P Parvizi: Simultaneous Morphological … 60 Fig. 4. Unrooted maximum-likelihood bootstrap tree showing the relationships of different genus and species of rodents by employing Cyt b gene fragment, including rodents submitted in GenBank and our Tatera indica sample using MEGA5 software. (KHT sample showed in the phylogenic tree is our Tatera indica sample captured from Khuzestan Province, southern Iran) Table 2. Morphological classification of Tatera indica based on their colors, Khuzestan Province, Iran T. indica Morphological Characters Buff Brown (cm)Buff Black (cm) 0.42-0.470.55-0.57Width of rostrum 4.26-4.474.75-4.88Occipitonasal length 3.91-4.144.38-4.51Condylbasal length 2.18-2.412.64-2.67Zygomatic width 0.60-0.760.73-0.78Least interorbital width 1.53-1.561.68-1.79Cranial width 1.93-1.952.13-2.11Length of nasal 1.17-1.261.33-1.36Length of diastema 0.75-0.760.91-0.93Length of anterior palatine foramina 1.14-1.191.28-1.29Length of tympanic bullae 0.61-0.660.76-0.79Width of tympanic bullae 0.61-0.640.67-0.75Upper cheekteeth 0.55-0.610.63-0.71Lower cheekteeth 1.71-1.871.98-1.99Height of skull 2.41-2.512.58-2.65Length of mandible Discussion Bates (1988) drew attention to the im- portance of systematic and zoogeoghraphy of genus Tatera of northeast Africa and Asia. The genus Tatera has an extensive geograph- http://jad.tums.ac.ir Published Online: June 27, 2015 J Arthropod-Borne Dis, March 2016, 10(1): 55–64 S Mohammadi and P Parvizi: Simultaneous Morphological … 61 ical distribution in Africa and Asia (Agrawal 2000, Yigit et al. 2001, Mirshamsi et al. 2007). There are reports of numerous studies on morphological and karyological taxon- omy of genus Tetera, however only few mo- lecular characterizations are available world- wide (Mirshamsi et al. 2007). T. indica was considered the main im- portant reservoir host of ZCL in Khuzestan (Mohebali and Javadian 2004, Mehrabani et al. 2007, Hajjaran et al. 2013). Principally P. papatasi may have been acquiring L. major from T. indica and other rodent reservoir hosts living peridomestically (Yaghoobi-Ershadi et al. 2013). Although different rodent species were found but they were limited and located in only one location. Only T. indica is considered as prominent reservoir host, which has a significant role in maintaining Leishmania parasites in different locations of ZCL in Khuzestan Province and this, high- lights the importance of current survey on T. indica. In this survey, T. indica was investigated more precisely in case of morphological features and molecular methods for the first time in Iran (Khuzestan Province) as the foremost and dominant reservoir in this focal region of leishmaniasis (Table 1, 2, Fig. 2, 3). The comprehensive study in broad-spec- trum ranges and distribution of T. indica was done in Khuzestan Province but this rodent was restricted to only six locations. The Cyt b gene was employed for mo- lecular typing and characterizing of T. indica for the first time in Iran and only one Cyt b gene sequence of T. indica is available from France in GenBank (GenBank accession no: AJ430563.1) and this indicates that suffi- cient molecular analyzing has not been car- ried out yet. Previous investigation of T. indica in Iran were only based on morpho- logical identification and/or on Leishmania infection of this reservoir (Etemad 1978, Mehrabani et al. 2011, Hajjaran et al. 2013) and there was not any study on determining the principle reservoir host, morphological features along with molecular systematics analyzing simultaneously. The interesting finding of this investiga- tion was that two types of morphological features (phenotypes) of T. indica was found in Khuzestan however after using Cyt b gene, only one haplotype from all 21 se- quences (11 buff- black and 10 buff-brown) was identified (Table 2, Fig. 4) It is not clear how authors identified four subspecies of T. indica in Iran and classified four subspecies (Mirshamsi et al. 2007). They did not give any definitive morpho- logical or molecular taxonomic characteri- zation for separating the T. indica subspecies. A phylogenetic analysis of the new se- quence, and those previously reported for T. indica and other rodents, found no support for recognizing more than one species or subspecies (T. indica) in Khuzestan Prov- ince, Iran (Fig. 4). Mitochondrial Cyt b demonstrated an ab- sence of any subspecies between two mor- photypes of T. indica, which classified pre- viously and this study in Khuzestan Province (Mirshamsi et al. 2007, Oshaghi et al. 2011). All four taxa (T. indica persica in Sistan, T. indica scansa in Kerman, T. indica monticola and T. indica bailwardi in Khuzestan) might be good biological species and did not char- acterize molecularly but could showed mito- chondrial introgression caused by occasional inter-breeding the same as conclusions for some sandfly sibling species are reported (Testa et al. 2002, Pesson et al. 2004, Parvizi et al. 2010). For resolving sibling T. indica species, a population genetics approach using several polymorphic genes might be tested by considering carefully for species-specific molecular markers and more specimens of T. indica species in large geographical locations. We only could conclude if two subspecies in Khuzestan Province or four subspecies of T. indica exist in Iran, might have good evi- dence for associating specific taxa with phe- http://jad.tums.ac.ir Published Online: June 27, 2015 J Arthropod-Borne Dis, March 2016, 10(1): 55–64 S Mohammadi and P Parvizi: Simultaneous Morphological … 62 notypes of epidemiological importance. There is no such evidence for T. indica subspecies, and so there is no reason to give priority to resolving the species status of four siblings T. indica species. Conclusion Tatera indica was characterized by both Cyt b molecular marker and morphological features for first time in Iran. There is no any morphotypes or subspecies of T. indica was found by molecular tools. We only could conclude if two subspecies in Khuzestan Province or four subspecies of T. indica ex- ist in Iran, might have good evidence for as- sociating specific taxa with phenotypes of epidemiological importance. There is no such evidence for T. indica subspecies, and so there is no reason to give priority to resolv- ing the species status of four siblings T. indica species. Acknowledgements The work was supported by the Pasteur Institute of Iran, grant 605 awarded to Dr Parviz Parvizi. We thank Dr Vazirianzadeh, Amraei, Kajkolah, Adel Spotin, Ali Bordbar, Sahar Ebrahimi, Javad Sameii, Mehdi Baghban, Roozbeh Taslimian for help with the field work and Elnaz AlaeeNovin and Narmin Najafzadeh for help in Molecular Systematics Laboratory. This research through a studentship to Miss Somayeh Mohammadi, based at the Pasteur Institute of Iran, Tehran, and registered for Islamic Azad medical University, Tehran, Iran. We confirm that there are no known conflict of interests associated with this publication. References Agrawal VC (2000) Taxonomic studies on In- dian Muridae and Hystricidae (Mam- malia: Rodentia). Rec Zool Surv India. 180: 1–177. Bordbar A, Parvizi P (2014) High density of Leishmania major and rarity of other mammals' Leishmania in zoonotic cu- taneous leishmaniasis foci, Iran. Trop Med Int Health. 19(3): 355–363. Chevret P, Dobigny G (2005) Systematics and evolution of the subfamily Gerbillinae Mammalia, Rodentia, Muridae. Mol Phylogenet Evol. 35(3): 674–688. Colangelo P, Civitelli MV, Capanna E (2005) Morphology and chromosomes of Tatera Lataste 1882 (Rodentia Muridae Gerbillinae) in West Africa. Trop Zool. 14(2): 243–253. Corbet GB (1978) The Mammals of the Pal- aearctic Region: A Taxonomic Review. British Museum (Natural History), London. Corbet GB, Hill JE (1991) A World List of Mammalian Species. 3th Edition. Natural History Museum Publications /Oxford University Press, UK. Ellerman JR (1948) Notes on some Asiatic rodents in the British Museum. Pro- ceedings of the Zoological Society of London. 117: 259–271. Etemad A (1978) Mammals of Iran, Rodents and Identification Key. Vol. 1, Natural Resource Protection and Human Envi- ronment Association press, Tehran [In Persian]. Hajjaran H, Hajjaran H, Mohebali M, Mamishi S, Vasigheh F, Oshaghi MA, Naddaf SR, Teimouri A, Edrissian GH, Zarei Z (2013) Molecular Identification and Pol- ymorphism Determination of Cutaneous and Visceral Leishmaniasis Agents Iso- lated from Human and Animal Hosts in Iran. BioMed Res Int 789326 (7). doi: 10.1155/2013/789326. Epub 2013 Oct 28. Harrison DL, Bates PJJ (1991) The Mam- mals of Arabia 2nd ed. Harrison Zool. Mus Publ Ket. http://jad.tums.ac.ir Published Online: June 27, 2015 J Arthropod-Borne Dis, March 2016, 10(1): 55–64 S Mohammadi and P Parvizi: Simultaneous Morphological … 63 Javadian E, Dehestani M, Nadim A, Rassi Y, Tahvildar-bidruni Gh, Seyedi-Rashti MA, Shadmehr A (1998) Confirmation of Tatera indica (Rodentia: Gerbillidae) as the main reservoir host of zoonotic cutaneous leishmaniasis in the west of Iran. Indian J Public Health. 27(1–2): 55–60. Kent RJ, Norris DE (2005) Identification of mammalian blood meals in mosquitoes by a multiplexed polymerase chain re- action targeting cytochrome b. Am J Trop Med Hyg. 73(2): 336–342. Kia EB, Shahryary-Rad E, Mohebali M, Mahmoudi M, Mobedi I, Zahabiun F, Zarei Z (2010) Endoparasites of Rodents and Their Zoonotic Importance in Germi, Dashte-Mogan, Ardabil Province. Iran. Iran J Parasitol. 5(4): 15–20. Maleki-Ravasan N, Oshaghi MA, Javadian E, Rassi Y, Sadraei J, Mohtarami F (2009) Blood meal identification in field-captured sand flies: comparison of PCR-RFLP and ELISA assays Iran. Iran J Arthropod-Borne Dis. 3(1): 8–18. Mehrabani D, Motazedian MH, Asgari Gholam Q, Hatam R. Owji SM, Oryan A (2011) Leishmania major in Tatera indica in Estahban, Southern Iran: Microscopy, culture, isoenzyme and PCR. Pak J Med Sci. 27(4): 734–773. Mehrabani D, Motazedian MH, Oryan A, Asgari Q, Hatam GR, Karamian M (2007) A search for the rodent hosts of Leishmania major in the Larestan re- gion of southern Iran: demonstration of the parasite in Tatera indica and Gerbillus sp., by microscopy, culture and PCR. Pathog Glob Health. 101(4): 315–322. Mirshamsi O, Darvish J, Kayvanfar NA (2007) preliminary study on Indian Gerbils, Tatera indica Hardwicke, 1807 at po- pulation level in eastern and southern parts of Iran (Rodentia: Muridae). Intl J Aquat Bio. 3(1): 49–61. Missone X (1959) Zoogeographie des Mammiferes de l’ Iran. Le Gerfaut Press, France. Mohebali M, Javadian E (2004) Infection in rodents from endemic areas of the Islamic Republic of Iran. East Mediterr Health J. 10(4–5): 591–599. Oshaghi MA, Rassi Y, Tajedin L, Abai MR (2011) Mitochondrial DNA diversity in the populations of great gerbils, Rhombomys opimus, the main reser- voir of cutaneous leishmaniasis. Acta Trop. 119: 165–171. Parvizi P, Tahehrkhani H, Ready PD (2010) Phlebotomus caucasicus and Phle- botomus mongolensis (Diptera: Psy- chodidae): indistinguishable by the mitochondrial cytochrome b gene in Iran. Bull Entomol Res. 100: 415–420. Parvizi P, Ready PD (2008) Nested PCRs and sequencing of nuclear ITS-rDNA fragments detect three Leishmania species of gerbils in sandflies from Iranian foci of zoonotic cutaneous Leishmaniasis. Trop Med Int Health. 13(9): 1159–1171. Pesson B, Ready JS, Benabdennbi I, Martı´n-Sa´nchez J, Esseghir S, Cadi- Soussi M, Morillas-Marquez F, Ready PD (2004) Sandflies of the Phlebotomus perniciosus complex: mitochondrial in- trogression and a new sibling species of P. longicuspisin the Moroccan Rif. Med Vet Entomol. 18: 25–37. Shirzadi M (2012) Instruction of Cutaneous Leishmaniasis care (Salak) in Iran. Min- istry of Health, Treatment and Medical Education, Deputy of Health, Admin- istration Center of Contagious Diseas- es, Office of Transmissible Disease between Humans and Animals. First Pub, Raze nahan, Tehran [In Persian]. http://jad.tums.ac.ir Published Online: June 27, 2015 J Arthropod-Borne Dis, March 2016, 10(1): 55–64 S Mohammadi and P Parvizi: Simultaneous Morphological … 64 Tamura K, Peterson D, Peterson N, Stecher G, Nei M, Kumar S (2011) MEGA5: Molecular Evolutionary Genetics Ana- lysis using Maximum Likelihood, Evo- lutionary Distance, and Maximum Par- simony Methods. J Mol Biol. 10: 2731– 2739. Testa JM, Montoya-Lerma J, Cadena H, Oviedo M, Ready PD (2002) Molec- ular identification of vectors of Leish- mania in Colombia: mitochondrial int- rogression in the Lutzomyia townsendi series. Acta Trop. 84: 205–218. Vazirianzadeh B, Saki J, Jahanifard E, Zarean M, Amraee K, Navid Pour S (2013) Iso- lation and Identification of Leishmania Species From Sandflies and Rodents Collected From Roffaye District, Khu- zestan Province, Southwest of Iran. J Microbiol. 6(6): e10025, DOI: 10.5812 /jjm.10025. WHO (2008) Leishmaniasis: Background of Information. WHO Geneva. Yaghoobi-Ershadi MR, Shahbazi F, Darvishi M, Akhavan AA, Jafari R, Khajeian M, Rassi Y, Soleimani H, Shirzadi MR, Hanafi-Bojd AA, Darabi H, Arandian MH, Sanei-Dehkordi A, Heidari M (2013) Molecular Epidemiological Study of Cutaneous Leishmaniasis in the Focus of Bushehr City, Southwestern Iran. J Arthropod Borne Dis. 7(2): 113–121. Yigit N, Colak E, Verimli R, Ozkurt S, Sozen M (2001) A study on the dis- tribution morphology and karyology of Tatera indica (Hardwicke, 1807) (Mam- malia: Rodentia) in Turkey. Turk J Zool. 25: 65–67. http://jad.tums.ac.ir Published Online: June 27, 2015