J Arthropod-Borne Dis, June 2015, 9(1): 84–97 Ar Zahraei-Ramazani et al.: Morphological and … 84 Original Article Morphological and Genotypic Variations among the Species of the Subgenus Adlerius (Diptera: Psychodidae, Phlebotomus) in Iran Alireza Zahraei-Ramazani 1, Dinesh Kumar 2, Hossein Mirhendi 3, Shyam Sundar 4, Rajnikan Mishra 2, Vahideh Moin-Vaziri 5, Hassan Soleimani 6, Mohammad Reza Shirzadi 7, Reza Jafari 8, Ahmad Ali Hanafi-Bojd 1, Sodabe Hamedi Shahraky 9, *Mohammad Reza Yaghoobi-Ershadi 1 1Department of Medical Entomology and Vector Control, School of Public Health, Tehran University of Medi- cal Sciences, Tehran, Iran 2Department of Zoology, Centre of Advanced Study, Banaras Hindu University, Varanasi 221 005, India 3Department of Medical Parasitology and Mycology, School of Public Health, Tehran University of Medical Sciences, Tehran, Iran 4Department of Medicine, Institute of Medical Sciences, Banaras Hindu University, Varanasi 221 005, India 5Department of Parasitology and Mycology, School of Medicine, Shahid Beheshti University of Medical Sci- ences, Tehran, Iran 6Yazd Health Research Station, National Institute of Health Research, Yazd, Iran 7Communicable Disease Management Centre, Ministry of Health and Medical Education, Tehran, Iran 8Esfahan Health Research Station, National Institute of Health Research, Esfahan, Iran 9School of Health, Zabol University of Medical Sciences, Zabol, Iran. (Received 16 Mar 2014; accepted 3 May 2014) Abstract Background: Female sand flies of subgenus Adlerius are considered as probable vectors of visceral leishmaniasis in Iran. The objective of this study was to determine the morphological and genotypic variations in the populations of this subgenus in the country. Methods: Sand flies collected using sticky traps from 17 provinces during 2008–2010. The morphometric measure- ments were conducted with an Ocular Micrometer. Data was analyzed by SPSS. The Cytb gene was used to estimate population genetic diversity and identify the female specimens. UPGMA phenetic tree was used for DNA haplotypes of Cytb gene. Results: Six species of subgenus Adlerius identified from which one species, P. (Adlerius) kabulensis, is new record. The identification key is provided for males. Results revealed the molecular systematic in the species of subgenus Adlerius and determine the relationship of three females of P. comatus, P. balcanicus and P. halepensis. Conclusion: The positions of three females and the males in the UPGMA tree are correct and the similarities among them confirm our results. The branches of each species are not genetically distinct which justify the overlapping morphological characters among them. Molecular sequencing of Cytb-mtDNA haplotypes can be used for female identification for different species of subgenus Adlerius in Iran. Keywords: Phlebotomine sand flies, subgenus Adlerius, P. comatus, P. kabulensis, Mitochondrial DNA Cyto- chrome b gene Introduction Zoonotic Visceral Leishmaniasis (ZVL) is a potentially fatal disease in human, pri- marily in children. It is caused by Leishma- nia infantum in which dogs, foxes, and jack- als are the main reservoir hosts. Four species of the subgenus Larrousius and one species of subgenus Paraphlebotomus are consid- ered as probable vector species in Iran (Yag- hoobi-Ershadi 2012). The disease is endemic in many rural communities of 7 out of 31 prov- *Corresponding author: Prof Mohammad Reza Yaghoobi-Ershadi, Email: yaghoobi.reza@gmail.com http://jad.tums.ac.ir Published Online: July 16, 2014 J Arthropod-Borne Dis, June 2015, 9(1): 84–97 Ar Zahraei-Ramazani et al.: Morphological and … 85 inces and also sporadically reported from other regions of the country (Yaghoobi-Ershadi 2012, Zahraei-Ramazani et al. 2013). The subgenus Adlerius Nitzulescu includes about 20 species (Seccombe et al. 1993) in which some are suspected or proven vectors of ZVL. In the old world these include: Phlebotomus longiductus in China, Ka- zakhstan, Kyrgyzstan, Ukraine, and Uzbeki- stan. P. halepensis ranges in Syrian Arab Re- public, Azerbaijan and Georgia (WHO 2010). It should be noted that P. (Adlerius) arabicus has also reported as the vector of L. tropica in Palestinian Territories (Svobodova et al. 2006a, 2006b). The subgenus Adlerius is almost entirely Palaearctic (Seyedi-Rashti and Nadim 1992) in various Eurasian countries and is thought to have originated from central Asia where biodiversity is most abundant (Léger and Pesson 1987). Females from about 20 spe- cies of the subgenus Adlerius can transmit Leishmania major (Sadlova et al. 2003, Volf and Myskova 2007), L. tropica (Jacobson et al. 2003, Sadlova et al. 2003, Svobodova et al. 2006a,Volf and Myskova 2007), L. dono- vani (Zhang and Leng 1997, WHO 2010) and L. infantum (Killick-Kendrick 1999, Volf and Myskova 2007, Giorgobiani et al. 2012). There is no published morphological key for species identification that are generally rec- ognized on the basis of male morphology (Lewis 1982, Seyedi-Rashti and Nadim 1992). These sand flies are found in diverse natural microhabitats, especially in moun- tainous areas, such as: cracks, crevices and around the nests and gatherings of the vis- ceral leishmaniasis reservoirs (Artemiev 1973). Most are also found indoors (Theodor and Mesghali 1964). Because of the short sunny season in the mountains and highland prov- inces, the collection time of the sand flies is too short. Prior to this study, the systematic (mor- phological and molecular) and molecular studies of some medically significant arthro- pods such as Aedes, Anopheles, Culex, mites, and sand flies have been studied by molecu- lar methods. The high sensitivity and accu- racy of the molecular methods related to DNA in terms of morphological identifica- tion of species are well established. These are undetectable such as females of the sub- genus Adlerius is important (Munsterman and Conn 1997, Artemiev 1980, Artemiev and Neronov 1984). According to previous studies, the Cytb- mtDNA is broadly used for sand flies sys- tematics (Esseghir et al. 2000, Parvizi and Ready 2006). The studies on genetic varia- tion and molecular systematics of the subge- nus Adlerius sand flies present evidence of ecological differences between the group. Access to the genetic variation and mo- lecular systematics of the subgenus Adlerius , in addition to assisting the systematic de- lineation and their classification may offer evi- dence to interpret and explain the ecological differences among the species. Species knowl- edge, abundance, and distribution of males and females, and changes in their population are useful for all research programs in con- trol of leishmaniasis. The objective of this study was to deter- mine the morphological and genotypic varia- tions among the species of the subgenus Adlerius in their distributional areas of Iran. Materials and Methods Study area Iran is an arid land of 1.6 million Km2, extending north to Armenia, Azerbaijan, Turkmenistan, and the Caspian Sea, east to Afghanistan and Pakistan , south to the Per- sian Gulf and Sea of Oman, and bordered by Turkey and Iraq on the west. Mountains span the nation in the shape of a large V. Be- tween these ranges lies a high plateau where flowing water from the mountains disap- pears into desert sand. http://jad.tums.ac.ir Published Online: July 16, 2014 J Arthropod-Borne Dis, June 2015, 9(1): 84–97 Ar Zahraei-Ramazani et al.: Morphological and … 86 The Caspian Sea littoral zone comprises the northern slopes of the Alborz Mountains and the Caspian plain. This is a narrow strip of land, forest covered, with Mediterranean cli- mate, with the average temperature ranging between 10 and 35o C and the average rela- tive humidity between 70 and 100%. The Central plateau, situated between the Alborz and Zagros ranges, are very moun- tainous in the northwest where the ranges originate and becomes a lower desert in the east. The climate is dry with average tem- perature between 0 and 40 oC, with hot, dry summers and cold, snow-laden winters. Sev- eral rivers originate on the southern slopes of the Alborz Mountains. The Persian Gulf littoral and the Khuzestan plain, to the south of the foothills of the Za- gros Mountains, are characterized as a tropi- cal climate (Fig. 1). The average temperature ranges between 12 and 50 oC. The average relative humidity ranges between 40% and 80%. The extent and distribution of the subge- nus Adlerius sand flies collection was con- ducted from 40o north to 29o south latitude and 44o western to 62o eastern longitude re- spectively (Statistical culture of Iran, 1990, personal communication). A detailed study of distribution of the subgenus Adlerius species was performed in rural areas of 17 of 31 provinces including [(1) Azarbaijan-e-Sharqi, (2) Azarbaijan-e- Gharbi, (3) Ardabil, (4) Kordestan, (5) Ghom, (6) Ilam, (7) Lorestan, (8) Khuzestan, (9) Bushehr, (10) Fars, (11) Sistan Va Balu- chestan, (12) Kerman, (13) Kohkiluyeh va Boyer Ahmad, (14) Chahar Mahal-o-Bak- htiari, (15) Esfahan, (16) Markazi, and (17) Khorasan-e-Razavi] (Fig. 1). Sand flies collection The sand flies were collected using 200 sticky paper traps (consisted of white sheets 15x21 cm, coated with castor oil) placed in different habitats and various biotypes. This included areas outdoors, such as gardens, mountain caves, animal shelters, wall cracks, burrows, tree holes, under stones, and rocks in 3–4 villages suitable for sand fly breed- ing. These were located in the plains and mountains areas within each province and sampling was conducted twice each year (June and September) from 2008–2010. Traps were placed before sunset and col- lected the next morning before the sunrise. Sand flies were preserved in 96 % ethyl al- cohol glass containers and were kept at 4 °C in the fridge. The number of the trap sample, date, and location was recorded on the glass containers (Moin-Vaziri et al. 2007). Mounting and identification For species identification, sand flies were mounted in Puri’s medium, which was manu- factured in our leishmaniasis laboratory in Iran and identified after 24 h (Yaghoobi-Er- shadi and Javadian 1997) using the keys of Theodor and Mesghali (1964), Artemiev (1973), and Lewis (1982). Morphological char- acters used for male identification of the sub- genus Adlerius included antennal segments, shape of the antennae and ascoid, number of hairs on coxite, length of style on ventral proc- ess, length of filament , length of pump, and sub-terminal tubercle of the aedeagus (Lewis 1982). Morphometric study The eight organs were measured for our morphometrical study included: (1) length of the third segment of antennae, (2) length of style, (3) length of epipharynx, (4) length of coxite, (5) width of style, (6) length of surstyle, (7) length of aedeagus and (8) number of hairs on coxite. The morphometric measure- ments were conducted with an ocular mi- crometer and Olympus Microscope (ch-2) (Fig. 2). Magnification (40X). Statistical analysis Statistical analyses were used to determine http://jad.tums.ac.ir Published Online: July 16, 2014 J Arthropod-Borne Dis, June 2015, 9(1): 84–97 Ar Zahraei-Ramazani et al.: Morphological and … 87 if morphometric characters among the mounted specimens in each species were similar with the species found in different geographical areas. We used SPSS-Version 17.0, and evaluated for differences with the T-test, ANOVA, and Kruskal-Wallis test with an a priori level of significance set at P< 0.05. If statistical tests were significantly different, the post Hoc-Bonferroni test was used. When the result of the Kruskal-Wallis test was significantly different, the Mann-Whit- ney test was used for comparing the results. Molecular study Sections of the abdomen, wings, and legs of the sand flies used for DNA extraction were stored in 1.5 ml sterile microtubes con- taining 96 % ethanol. DNA from 120 speci- mens were extracted, of which 50 were used for mtDNA-PCR having 16 sequences (Ta- ble 1). The washed, dried body of each sand fly was frozen then defrosted with liquid ni- trogen and homogenised and ground with a heavy cylindrical metal and shaking vigor- ously. The content was washed with 300 ml lysis buffer [10 mM tris, 25 mM EDTA (pH= 8), 1 % SDS, 25 mM NaCl, 2% Tri- ton×-100] and boiled for 10 minutes. This was followed by extraction with a Phenol- Chloroform method. After the ethanol pre- cipitation, DNA was dissolved in 15 µ l dis- tilled water and stored at -20 °C. A 25 µ l PCR reaction mixture consisted of 12.5 µ l Premix 2x[10x PCR buffer, Mgcl2, 10 Mm dNTPs and DNA Taq poly- merase], 1 µ l of CB3-PDR forward primer [5'-CA (Y= T/C) ATTCAACC (W= A/T) GAATGATA-3'], 1 µ l of NIN-PDR Reverse primer [5'-GGTA (Y= C/T) (W= A/T) TTG CCTC GA (W =T/A) TTCG (T/A) TATGA- 3'], 5.5 µ l ddH2O and 5µ l of sand fly ge- nomic DNA were used for amplifying of 550 bp. The PCR amplification was carried out with the following thermal profile using a Gene Amp® PCR System 2700 thermal cy- cle (AB Applied Biosystems): 6' min. for initial denaturation at 95 °C, 35 cycles of denaturation at 94 °C for 45 sec., annealing at 55 °C for 1 min., extension at 72 °C for 1 min. Sequence alignment and analysis In this study the chromatogram of the DNA sequencing samples were edited by using BioEdit software. The program Clustal W (Thompson et al. 1997) was used to study the alignment sequences and similarity scores. Our sequences were compared with the se- quences available in GenBank by using Blast available on, www.ncbi.nlm.gov/. The Blast 2 sequences were used to find multiple local alignments and detect the best homolo- gous between the female and the male specimens. DNA haplotypes of Cytb-mtDNA gene were used to build an UPGMA tree. We measured the eight different morpho- logical characters in all specimens of each Adlerius species which captured by sticky traps from outdoors of the provinces, and then the similarity scores (%) have been iden- tified by sequence comparing alignment pair- wise. Therefore, if the more similarity score and less E-value were identified between the pairwise comparison, it showed that they are same species. Also, the haplotypes have been confirmed by the polymorphic sites. Results Composition, frequency and distribution of the Adlerius species A total of 9,319 sand flies were collected, of which 3,847 (41.3%) and 5,472 (58.7%) were from the genus Phlebotomus and Sergentomyia, respectively. From the total number of sand flies of the genus Phlebotomus, 167 (1.8%) were from the subgenus Adlerius [26 (15.57%) were female and 141(84.43%) were male]. Of the males, Phlebotomus (Adlerius) brevis 8 (5.67%), P. (Adl.) halepensis 66 (46.81%), P. (Adl.) longiductus 30 (21.28%), P. (Adl.) balcanicus 23 (16.31%), P. (Adl.) kabulensis http://jad.tums.ac.ir Published Online: July 16, 2014 J Arthropod-Borne Dis, June 2015, 9(1): 84–97 Ar Zahraei-Ramazani et al.: Morphological and … 88 13 (9.22%) and P. (Adl.) comatus 1(0.71%). The P. (Adl.) kabulensis Artemiev, 1978 is a new record from Iran. In case of P. (Adl.) comatus which was reported by the authors as a new record by morphological characters from the country (Zahraei-Ramazani et al. 2013) reconfirmed by the molecular studies (Fig. 3). We updated the morphological key for the male subgenus Adlerius species of Iran as follows: (1) One ascoid on antennal segments 9–15. Antenna 8 with two ascoids..................... (2). -One ascoid on antennal segments 9–15. An- tenna 8 with one ascoid. Antenna 6 and 7 with one ascoid or with one long and one short ........................................................ (6). (2) Coxite with 14–27 hairs in group, rarely 29 ............................................................ (3). -Coxite with 29–115 hairs in group, rarely 27. Part of 27–85 hairs of group on distal half of coxite .......................................... (4). (3) Whole hair-group on basal half of coxite. Tubercle of aedegus 6–16 μm from tip ......... ............ P. turanicus (Akhoundi et al. 2011). -Part of hair-group on distal half of coxite. Tubercle of aedegus 19–28 μm from tip ....... ....... P. brevis (Theodor and Mesghali 1964). (4) Aedegus with rectangular subterminal notch ................ P. halepensis (Adler 1946). -Aedegus with normal obtuse-angle subterminal notch ....................................................... (5). (5) Coxite with 27–50 group-hairs ............... .. P. kabulensis (Probably P. cf. kabulensis). -Coxite with 50–85 group-hairs. Sperm tubes are long (1200–1700 μm) .............................. ........................ P. longiductus (Theodor and Mesghali 1964) (Also P.cf. longiductus). (6) Coxite with 90–220 group-hairs. Ventral process of style long (About 20 μm ........ (7). -Coxite with 30–85 group-hairs. Ventral proc- ess of style long or short ................................ ......... P. salangensis (Akhoundi et al. 2011). (7) Coxite very wide, whole group of 125– 200 hairs on its basal half .............................. .................. P. comatus (Zahraei et al. 2013). -Coxite narrow, part of hair-group on its dis- tal half .................................... P. balcanicus (Theodor and Mesghali 1964) (Also P.cf. balcanicus). The male and female sand flies of subge- nus Adlerius were collected from 12 prov- inces and the male were collected from 11 out of 17 provinces. The males were not col- lected from Markazi, Bushehr, Kerman, Kordestan, Ghom, and Kohkiluyeh va Boyer Ahmad. The maximum number collection was from Ardabil (52 specimens) and the minimum number was from Khuzestan and Azerbaijan-e-Gharbi (1 specimen) Provinces. Statistical analysis of eight morphological characters of the males of subgenus Adlerius sand flies The subgenus Adlerius specimens col- lected from some provinces were very few in number. So for the statistical comparison of the 8 measureable morphological characters of the specimens of P. halepensis, we sorted them according to the four geographical re- gions (North-West, North-East, South-West and South-East). There was no significant difference between the 8 characters in the specimens. The P. halepensis specimens col- lected from three provinces (Ardabil, Azer- baijan-e-Sharqi and Esfahan) showed sig- nificant difference between the length of style and length of aedeagus. The statistical analysis of eight morpho- logical characters of 30 males of P. balcani- cus specimens confirmed that there is no significant difference between the characters in the specimens of Ardabil and Azerbaijan- e-Sharqi. There is a significant difference between the lengths of the coxite of the P. kabulensis specimens in the three provinces (Esfahan, Lorestan and Ilam). The Post Hoc-Bonferroni test shows that the length of coxite of P. ka- http://jad.tums.ac.ir Published Online: July 16, 2014 J Arthropod-Borne Dis, June 2015, 9(1): 84–97 Ar Zahraei-Ramazani et al.: Morphological and … 89 bulensis specimens in Lorestan Province is significantly more than Ilam. The statistical analysis between the 7 specimens of P. brevis in Chahar Mahal-o-Bakhtiari and 1 specimen in Azarbaijan-e-Gharbi did not confirm any significant difference. The statistical analysis of 30 P. longiductus specimens in four geographical regions (North- West, North-East, South-West and South- East) confirmed that there are significant dif- ference value between the length of surstyle, the length of aedeagus, the number of hairs on coxite, the length of coxite, and the length of styles (P< 0.05) Molecular studies on the subgenus Adlerius species The PCR experiments amplified a frag- ment of 550 bp of the mitochondrial genome and it was successfully achieved for 6 spe- cies of the subgenus Adlerius (Table 1, 2). These were registered on JX885982 to JX 885998 Accession numbers in the GenBank. In the pair wise comparison of the speci- mens, we focused on the genetic similarities of the specimens among each other and de- termined the level of such similarities of male and female specimens and identified the females among 6 male species. The pair wise comparison in AZV8 fe- male specimens with nucleotides of six male species showed that the maximum similarity, i.e. 98% is found in P. comatus. In total, 528 nucleotides were compared with these 2 se- quences which show 8 (1.5%) polymorphic sites and it is confirmed that they are 2 hap- lotypes (Table 3). The pair wise comparison in males of the P. balcanicus species with nucleotides of the three female specimens shows that there are minimum 98% and maximum 99% similari- ties among the male and the SHC11 female specimen. It shows 22 (4.3%) polymorphic sites and it confirms that there are 4 haplo- types (Table 3). The nucleotides comparison of the 5 males of the P. halepensis specimens with the CKHE1342 female specimen shows that there are minimum 90.1% and maximum 92.6% similarities between them .The poly- morphic sites is 46 (10%) and the sequence comparison revealed that there are 6 haplo- types. Also for the two males of P. brevis, it shows 99% similarities with 2 (0.4%) poly- morphic sites and 2 haplotypes between these two sequenced specimens (Table 3). Cytb-mtDNA Sequences results of the sub- genus Adlerius species The sequence result of the length 550 bp of Cytb gene shows that 5 male species: P. halepensis, P. brevis, P. kabulensis, P. co- matus and P. balcanicus are in the areas of our study. Also two male specimens: P. cf. balcanicus and P. cf. longiductus which are close to P. balcanicus and P. longiductus are on neighbor-joining tree (Fig. 4). The speci- mens of P. halepensis and P. brevis are in independent and different branches in the tree, but other species are close together in one branch. P. comatus and the female AZV8 specimen are in the same branch, which shows that they are the same species. Also the branch of female SHC11 specimen is near the branch of AZV4 and AZV5, which shows that females of this species are close to these specimens. The branches of two P. balcanicus specimens i.e. Ash10 and AZV11 are close together and P. kabulensis is in a different branch. The position in the phenetic tree of the 3 females species and the male species are correct and the similari- ties among the males and females confirm our above results (Fig. 4). http://jad.tums.ac.ir Published Online: July 16, 2014 J Arthropod-Borne Dis, June 2015, 9(1): 84–97 Ar Zahraei-Ramazani et al.: Morphological and … 90 Table 1. Number of sand flies used in molecular studies NO. Species Number of specimens that DNA has been extracted Number of specimens that have been mtDNA-PCR Number of mtDNA-PCR that have been sequenced 1 P. halepensis 55 22 5 2 P. longiductus 23 5 1 3 P. kabulensis 12 5 1 4 P. balcanicus 13 10 3 5 P. brevis 8 3 2 6 P. comatus 1 1 1 7 P. (Adlerius) Female 8 4 3 Total 120 50 16 Table 2. Profile of sequenced species based on Cytb gene in the subgenus Adlerius Species No. Code no. Sex Collecting area Township Province GenBank ac- cession no. P. halepensis 1 ARM754 male Niaz village Meshkin- shahr Ardabil JX885982 2 AZV2 male Ahar-Varzaghan road Varzaghan Azarbaijan-e- Sharqi JX885984 3 EK77 male Konjan Natanz Esfahan JX885991 4 FJ103 male mountainous re- gion Jahrom Fars JX885992 5 MSH458 male Ghortapeh Meshkin- Shahr Ardabil JX885995 P. balcanicus 1 ASH10 male Sharbian Sarsb Azarbaijan-e- Sharqi JX885983 2 AZV11 male Ahar-Varzaghan road Varzaghan Azarbaijan-e- Sharqi JX885989 3 AZV5 male Ahar-Varzaghan road Varzaghan Azarbaijan-e- Sharqi JX885986 P. longiductus 1 AZV4 male mountainous re- gion Tabriz Azarbaijan-e- Sharqi JX885985 2 FJ476 male mountainous re- gion Jahrom Fars JX885993 P. kabulensis 1 LP7 male mountainous road Poldokhtar Lorestan Not clear 2 LPM5 male mountainous road Poldokhtar Lorestan JX885994 P. brevis 1 SHC7 male Chelgerd Koohrang Chahar Mahal-o- Bakhtiari JX885996 2 SHC22 male Chelgerd Koohrang Chahar Mahal-o- Bakhtiari JX885998 P. Adlerius group 1 ARM911 famele Niaz village Meshkin- shahr Ardabil Not clear 2 AZV8 famele Ahar-Varzaghan road Varzaghan Azarbaijan-e- Sharqi JX885987 3 CKHE1342 famele Irandegan Khash Sistan va Balu- chestan JX885990 4 SHC11 famele Chelgerd Koohrang Chahar Mahal-o- Bakhtiari JX885997 P. comatus 1 AZV10 male Ahar-Varzaghan road Varzaghan Azarbaijan-e- Sharqi JX885988 http://jad.tums.ac.ir Published Online: July 16, 2014 J Arthropod-Borne Dis, June 2015, 9(1): 84–97 Ar Zahraei-Ramazani et al.: Morphological and … 91 Fig. 1. Provinces where Adlerius species were collected. (http://commons.wikimedia.org/wiki/File%3AMy-iran-climate-map-simplified.png) Fig. 2. Organs measured in morphometrical studies of the subgenus Adlerius in Iran Caspian Mild Mountains Desert and Semi-Desert 1. The length of the third segment of antennae (A3) 5. The length of the style2. The length of Epipharynx 3. The length of coxite 8. The numbers of hairs on coxite 7. The Length of Aedeagus 6. The width of style 4. The Length of Surstyle http://jad.tums.ac.ir Published Online: July 16, 2014 J Arthropod-Borne Dis, June 2015, 9(1): 84–97 Ar Zahraei-Ramazani et al.: Morphological and … 92 Fig. 3. Number of males of subgenus Adlerius species collected from 17 provinces Fig. 4. Neighbor-joining tree for DNA haplotypes of Cytb-mtDNA of the subgenus Adlerius sand flies species P. brevis N u m b er P. halepensis P. longiductus P. balcanicus P. kabulensis P. comotus Male of Adlerius species 8 66 30 23 13 1 70 60 50 40 30 20 10 0 http://jad.tums.ac.ir Published Online: July 16, 2014 J Arthropod-Borne Dis, June 2015, 9(1): 84–97 Ar Zahraei-Ramazani et al.: Morphological and … 93 Table3. Comparison of Cytb gene nucleotide sequences in the subgenus Adlerius species Discussion Phlebotomine sand flies of the subgenus Adlerius transmit Leishmania parasites in various Eurasian countries. But, the species involved have seldom been identified be- cause despite the taxonomic efforts of Ar- temiev (1980) and others. The morphological characteristics do not distinguish between the females of Adlerius species that are gen- erally recognized on the basis of male mor- phology (Lewis 1982, Sadlova et al. 2003). The detection of P. comatus (Zahraei- Ramazani et al. 2013) and P. kabulensis during the present study has assisted to revise the morphological key. The access to new mo- lecular techniques and genetic information about these probable vectors (Artemiev 1973, Ardehali et al. 1995) and their systematics, may aid in resolving the problem of identi- fication in morphological similarity, as well. We provide evidence to interpret and explain their morphology and ecology in detail. The systematics of the Phlebotomine sand flies due to the nature of their diagnostic char- acteristics is highly variable. The latest revision of the genus Phlebotomus sand flies in Old World is included: (1) Adlerius, (2) Anaphle- botomus, (3) Euphlebotomus, (4) Kasaulius, (5) Transphlebotomus, (6) Synphlebotomus, (7) Phlebotomus, (8) Larroussius, (9) Paraphle- botomus (Rispail and Leger 1998). Secombe et al. (1993) believed that the genus Phle- botomus had few morphological deficiencies. According to the morphological character- istics they are monophyletic. The UPGMA is one of the most popular methods in ecology for the classification of sampling units on the basis of their pairwise similarities in relevant descriptor variables (such as species composition) (Legendre and Legendre 1998). In bioinformatics, UPGMA is used for the creation of phenetic trees (phenograms). In our molecular analysis, all of the Cytb gene sequences confirm the mo- nophyletic of the subgenus Adlerius and the morphological identification of these six spe- cies are confirmed by the reliable morpho- logical identification keys. The sequence analysis of P. halepensis populations showed 2 lineages in different geographical regions. The haplotypes are 6 in this species. According to our morphologi- cal, statistical, and molecular findings, three specimens: AZV2, MSH458 and ARM754 which were collected from the North-West of Iran, fall in one branch along with three P. halepensis of GenBank. Arm754 and MSH458 were captured from Ardabil and AZV2 was captured from Azerbaijan-e-Sharqi, two dif- ferent areas. In the statistical analysis, the length of the style in P. halepensis speci- mens in Ardabil Province was significantly different from that of Esfahan and Azerbai- jan-e-Sharqi. AZV2 is in a different branch but is closer to MSH458 and ARM754 speci- mens. FJ103 and EK77 are in one branch too. In the statistical analysis, the mean of length of style of EK77 is more than that of MSH458, ARM754 and AZV2. So in the Species Number of sequence Identity of sequence (%) Number of nucleo- tides compared Number of haplotype Polymorphic site (%) P. comatus 2 98 528 2 1.5 P. balcanicus 4 98-99 516 4 4.3 P. halepensis 6 90.1-92.6 458 6 10 P. brevis 2 99 478 2 0.4 P. kabulensis 1 - - 1 - P. cf. longiductus 1 - - 1 - http://jad.tums.ac.ir Published Online: July 16, 2014 J Arthropod-Borne Dis, June 2015, 9(1): 84–97 Ar Zahraei-Ramazani et al.: Morphological and … 94 UPGMA tree it shows that they are close to the FJ103 specimen in Fars Province. The FJ103 and the three P. halepensis specimens from the GenBank found in the north of Iran are closer or these specimens in two different areas and are near to their phenetic characters. Some geographical event such as desert winds in the central plateau (Fig. 1) causes movement and mixing in its population. The pairwise comparison in males of P. halepensis specimens with nucleotides of CKHE1342-female shows that there are minimum 90.1 % and maximum 92.6 % simi- larity among them. In statistical analysis, the mean of the length of labrum, surstyle, and aedeagus of the P. halepensis specimens in the south-east (CKHE1342) of Iran are dif- ferent from the specimens of other geographi- cal regions. However, CKHE1342-female is in different branch because of the geographi- cal distance between these specimens. In this study all the P. halepensis specimens are on different branch and are in a larger group. It also seems that this species has been sepa- rated from another species of Adlerius group many years ago than the other Adlerius spe- cies. According to the statistical analysis and molecular findings, further studies are re- quired to resolve the taxonomic status of this species. The present study on molecular analysis of P. brevis population shows that one line- age with two haplotypes for Cytb. The P. brevis (Theodor and Mesghali 1964) was firstly found in Iran. The SHC7 and SHC22 were captured from Chahar Mahal-o-Bak- htiari Province. Statistical analysis did not show any difference between them and the specimens from Azerbaijan-e-Gharbi. Along with them, another species i.e 306450786 was submitted to the GenBank from Iran which is in one independent branch. It seems that P. brevis is not a common ancestor of the other Adlerius but diverged early and shares one common ancestor with other species. The second group of the phylogeny tree including LPM5, AZV11, ASh10, AZV5, AZV4, AZV10 and two females specimens: SHC11 and AZV8 show that Cytb gene does not have sufficient resolution for them and indicates the morphological character over- lapping. The length of the branches in this group is short, which indicates that the changes among them are low. This can also be said that the specimens of each species are not far genetically. In statistical analysis by using t-test, we did not find any difference between the specimens of P. balcanicus in Ardabil and Azerbaijan-e-Sharqi Provinces. In the molecu- lar study, the pairwise comparison in three males of this species with nucleotides of SHC11-female specimen shows that there are minimum 98 % and maximum 99 % simi- larities among themselves. Also, it is con- firmed that there are 4 haplotypes among the four sequenced specimens. In the present study, AZV4 and AZV5 have been identified by using both morphological and molecular methods. However, in our UPGMA tree, we found that they are differents between P. balcanicus specimens and they are in an in- dependent branch but closer to P. balcanicus branches. The pairwise distance between them is zero. It is clear that the taxonomy of these species should be revised. The results of the molecular investigations indicate that according to the concept of the species they (AZV5 and AZV4) cannot be considered as a separate species. Because the concept of a species is based on the reproductive isola- tion, while the existence of the intermediate forms and the same mitochondrial gene se- quences in this study indicates that there is no reproductive isolation between AZV5 morphotypes and AZV4 morphotypes. Also, we cannot assume these two specimens as subspecies because two subspecies cannot be found as sympatric. So, this indicates that these species e.g P. longiductus and P. bal- canicus in Iran, including their populations http://jad.tums.ac.ir Published Online: July 16, 2014 J Arthropod-Borne Dis, June 2015, 9(1): 84–97 Ar Zahraei-Ramazani et al.: Morphological and … 95 have morphological and genetic differences. The similarity between AZSH10 and AZV11 is 97% and E-value between them is 0.0. Further, in the statistical analysis, no differ- ence was found between them. The branches of the four haplotypes of P. balcanicus show that genetically they are not far. The simi- larity and E-value between the SHC11 and AZV5 are 99 % and 0.0 respectively. So in the UPGMA tree they are close together. The pairwise comparison of the nucleo- tides of Cytb-mtDNA gene in AZV8 female specimen with nucleotides of six male spe- cies shows that the maximum similarity i.e 98% is found in AZV10. Results of the mo- lecular studies reveal that there are 2 haplo- types. These two specimens are in the same branch in the UPGMA tree. The pairwise distance between them is 0.0. So, these evi- dences confirm our results as both AZV8 and AZV10 are P. comatus. In the UPGMA tree, P. kabulensis is seen as independent branch. The length of the branch in this group compared with other branches is long. This indicates that many changes have occurred during the course of evolution of this species. Also this species seems a common ancestor for the other spe- cies in this group of the UPGMA tree. The location and the number of hairs on coxite is the most important morphological characteristic feature for identifying and sepa- rating one Adlerius species from the other. At times it is difficult for the researchers as the numbers and locations of the hairs of dif- ferent species are the same and overlap among themselves. In P. kabulensis, the hairs on coxite are 27–50 and in P. longiductus the hairs on coxite are 50–85. If we collect a specimen which has 50 hairs on its coxite, we cannot identify that it belongs to which one of the afore noted two species. Also for P. balcanicus and P. longiductus, we en- countered the same difficulty. In this case, the researchers use ‘cf’. The ‘cf’ is the ab- breviation of the Latin word ‘confers’ means “to compare” or “to consult”. It is mainly used in academic writings to indicate a ref- erence to a contrasting finding or viewpoint (it is often used instead of “but see”). It can also appear occasionally in binomial nomen- clature by placing before the species name to indicate that the species is not confirmed (Strunk and White 1979). To arrive at a re- sult, we use P. cf. longiductus or P. cf. bal- canicus. In our study, these species have been confirmed by using both morphological and molecular methods. But in our phenetic tree it is found that they are different from P. balcanicus and P. longiductus. But in an in- dependent group they are close to P. bal- canicus branch. Obviously, detailed studies for detection of genetic markers to identify the females of six Adlerius species are essential. Also, fur- ther studies on morphology and genetic di- versity of the populations of the subgenus Adlerius will be of particular importance. Conclusions Phlebotomus (Adlerius) kabulensis is new record for the Iran. The positions of three female and males of subgenus Adlerius in the UPGMA tree are correct and the simi- larities among them confirm our results. The branches of each species are not genetically distinct that justify the overlapping mor- phological characters among them. Molecular sequencing of Cytb-mtDNA haplotypes can be used for female identification of different species of subgenus Adlerius in Iran. Acknowledgements The authors thank the Medical Entomology and leishmaniasis staff of Esfahan Health Research Station, National Institute of Health Research for their contribution in sand fly collection. The study was funded by the Dep- uty of Research, Tehran University of Medi- http://jad.tums.ac.ir Published Online: July 16, 2014 J Arthropod-Borne Dis, June 2015, 9(1): 84–97 Ar Zahraei-Ramazani et al.: Morphological and … 96 cal Sciences through project number: 5146- 27-01-86. The authors of study declare that they have no competing interests. References Ardehali SV, Rezaei HR, Nadim A (1995) Leishmania Parasite and Leishmani- asis. Academic Press Center, Tehran, Iran, pp. 208 (in Persian). Artemiev MM (1973) Sand flies (Diptera: Psychodidae, Phlebotominae) of Af- ghanistan. Malaria and Leishmania In- stitute, Afghanistan, p. 174. Artemiev MM (1980) A review of sand flies of the subgenus Adlerius (Diptera: Psychodidae: Phlebotominae). Zo- ologicheskiy Zh. 59: 1177–1192 (Rus- sian). Artemiev MM, Neronov V (1984) Distribu- tion and ecology of sand flies of the world (genus Phlebotomus). Institute of Evolution, Morphology and Animal Ecology, Moscow, pp. 1–208. Esseghir S, Ready RD, Ben-Ismail R (2000) Speciation of Phlebotomus sand flies of the subgenus Larroussius coincided with the late Miocene-Pliocene aridifi- cation of the Mediterranean subregion. Biol J Linnean Soc. 70: 189–219. Giorgobiani E, Lawyer PG, Babuadze G, Dolidze N, Jochim RC, Tskhvaradze L, Kikaleishvili K, Kamhawi S (2012) Incrimination of Phlebotomus kande- lakii and Phlebotomus balcanicus as vectors of Leishmania infantum in Tbi- lisi, Georgia. PLoS Neglected Trop Dis. 6: e1609. Jacobson RL, Eisenberger CL, Svobodova M, Baneth G, Sztern J, Carvalho J, Nasereddin A, El-Fari M, Shalom U, Volf P, Votypka J, Dedet JP, Pratlong F, Schonian G, Schnur LF, Jaffe CL, Warburg A (2003) Outbreak of cuta- neous leishmaniasis in Northern Israel. J Infect Dis. 188: 1065–1073. Killick-Kendrick R (1999) The biology and control of Phlebotominae sand flies. Med Vet Entomol. 17: 279–289. Legendre P, Legendre L (1998) Numerical Ecology. Second English Edition. De- velopments in Environmental Modelling 20. Elsevier, Amsterdam, pp. 1–853. Léger N, Pesson B (1987) Sur la taxonomie et la repartition géographique de Phle- botomus (Adlerius) chinensis s.l. et de P. (Larroussius) major s.l. (Psychodi- dae: Diptera). Bull Soc Path Exot. 80: 252–260. Lewis DJ (1982) A taxonomic review of the genus Phlebotomus (Diptera: Psychodi- dae). Bull Br Mus Nat Hist Entomol. 45: 121–209. Moin-Vaziri V, Depaquit J, Yaghoobi-Ershadi MR, Oshaghi MA, Derakhshandeh- Peykar P, Ferte H, Kaltenbach M, Bargues MD, Leger N, Nadim A (2007) Interaspecific variation within Phlebotomous sergenti Parrot, 1917 (Diptera: Psychodidae) based on mtDNA sequences in Islamic Republic of Iran. Acta Trop. 102: 29–37. Munstermann LE, Conn JE (1997) Systematics of mosquito disease vectors (Diptera: Culicidae): Impact of molecular biolo- gy and cladistic analysis. Annu Rev Entomol. 42: 351–369. Parvizi P, Ready PD (2006) Molecular in- vestigation of the population differen- tiation of Phlebotomus papatasi, im- portant vector of Leishmania major, in different habitats and regions of Iran. Iranian Biomed J. 10: 69–77. Rispail P, Leger N (1998) Numerical taxon- omy of old world Phlebotominae (Dip- tera: Psychodidae). 1. Considerations of morphological characters in the Ge- nus Phlebotomus Rondani and Berte 1840. Mem Inst Oswaldo Cruz. 93: 773–785. http://jad.tums.ac.ir Published Online: July 16, 2014 J Arthropod-Borne Dis, June 2015, 9(1): 84–97 Ar Zahraei-Ramazani et al.: Morphological and … 97 Sadlova J, Hajmova M, Volf P (2003) Phle- botomus (Adlerius) halepensis vector competence for Leishmania major and L. tropica. Med Vet Entomol. 17: 244– 250. Secombe AK, Ready PD, Huddleston LM (1993) A catalogue of the old world Phlebotominae sand flies (Diptera: Psy- chodidae, Phlebotominae). Department of Entomology the Natural History Museum. Occasional Papers on Sys- tematics Entolomology. 8: 1–57. Seyedi-Rashti MA, Nadim A (1992) The ge- nus Phlebotomus (Diptera: Psychodi- dae: Phlebotominae) of the countries of the Eastern Mediterranean Region. Iranian J Publ Health. 21: 11–50. Statistics Culture of Iran, 1990. Annual Re- port on the Statics of the Country. Vol. 6. Strunk WJ, White EB (1979) The elements of style. (3rd Ed). New York: Macmillan, (Chapter 4). Svobodova M, Volf P, Votýpka J (2006a) Experimental transmission of Leishma- nia tropica to hyraxes (Procavia cap- ensis) by the bite of Phlebotomus ara- bicus. Microbes Infect. 8: 1691–1694. Svobodova M, Votypka J, Peckova J, Dvorak V, Nasereddin A, Baneth G, Sztern J, Kravchenko V, Orr A, Meir D, Schnur LF, Volf P, Warburg A (2006b) Dis- tinct transmission cycles of Leishmania tropica in 2 adjacent foci, Northern Is- rael. Emerg Infect Dis. 12: 1860–1868. Theodor O, Mesghali A (1964) On the Phle- botominae of Iran. J Med Entomol. 1: 285–300. Thompson JD, Gibson TJ, Plewniak F, Jeanmougin F, Higgins DG (1997) The CLUSTAL_X windows interface: flexi- ble strategies for multiple sequence alignment aided by quality analysis tools. Nucl Acids Res. 25: 4876–4882. Volf P, Myskova J (2007) Sand flies and Leishmania: species versus permissive vectors. Trends Parasitol. 23: 91–92. World Health Organization (2010) Control of the Leishmaniases. WHO Technical Report Series, 949, Report of a meeting of the WHO Export Committee on the control of leishmaniases, 22–26 March 2010, Geneva, p. 186. Yaghoobi-Ershadi MR (2012) Phlebotomine sand flies (Diptera: Psychodidae) in Iran and their role on Leishmania transmis- sion. J Arthropod Borne Dis. 6: 1–17. Yaghoobi-Ershadi MR, Javadian E (1997) Studies on sandflies in a hyperendemic area of zoonotic cutaneous leishmani- asis in Iran. Indian J Med Res. 105: 61–66. Zahraei-Ramazani AR, Kumar D, Yaghoobi- Ershadi MR, Naghian A, Jafari R, Shirzadi MR, Abdoli H, Soleimani H, Shareghi N, Ghanei M, Arandian MH, Hanafi-Bojd AA (2013) Sand flies of the subgenus Adlerius (Diptera: Psy- chodidae) in an endemic focus of vis- ceral leishmaniasis and introduction of Phlebotomus (Adlerius) comatus as a new record for Iran. J Arthropod-Borne Dis. 7: 1–7. Zhang LM, Leng YJ (1997) Eighty-year re- search of Phlebotomine sandflies (Dip- tera: Psychodidae) in China (1915– 1995). II. Phlebotomine vectors of leishmaniasis in China. Parasite. 4: 299–306. http://jad.tums.ac.ir Published Online: July 16, 2014