zoonotic visceral leishmaniasis ( zvl ) is a potentially fatal disease in human , primarily in children . it is caused by leishmania infantum in which dogs , foxes , and jackals are the main reservoir hosts . four species of the subgenus larrousius and one species of subgenus paraphlebotomus are considered as probable vector species in iran ( yaghoobi - ershadi 2012 ) . the disease is endemic in many rural communities of 7 out of 31 provinces and also sporadically reported from other regions of the country ( yaghoobi - ershadi 2012 , zahraei - ramazani et al . 1993 ) in which some are suspected or proven vectors of zvl . in the old world these include : phlebotomus longiductus in china , kazakhstan , kyrgyzstan , ukraine , and uzbekistan . p. halepensis ranges in syrian arab republic , 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 . 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 ( lger and pesson 1987 ) . females from about 20 species 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. donovani ( 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 recognized on the basis of male morphology ( lewis 1982 , seyedi - rashti and nadim 1992 ) . these sand flies are found in diverse natural microhabitats , especially in mountainous areas , such as : cracks , crevices and around the nests and gatherings of the visceral leishmaniasis reservoirs ( artemiev 1973 ) . because of the short sunny season in the mountains and highland provinces , the collection time of the sand flies is too short . prior to this study , the systematic ( morphological and molecular ) and molecular studies of some medically significant arthropods such as aedes , anopheles , culex , mites , and sand flies have been studied by molecular methods . the high sensitivity and accuracy of the molecular methods related to dna in terms of morphological identification of species these are undetectable such as females of the subgenus adlerius is important ( munsterman and conn 1997 , artemiev 1980 , artemiev and neronov 1984 ) . according to previous studies , the studies on genetic variation and molecular systematics of the subgenus adlerius sand flies present evidence of ecological differences between the group . access to the genetic variation and molecular systematics of the subgenus adlerius , in addition to assisting the systematic delineation and their classification may offer evidence to interpret and explain the ecological differences among the species . species knowledge , abundance , and distribution of males and females , and changes in their population are useful for all research programs in control of leishmaniasis . the objective of this study was to determine the morphological and genotypic variations among the species of the subgenus adlerius in their distributional areas of iran . iran is an arid land of 1.6 million km 

 , extending north to armenia , azerbaijan , turkmenistan , and the caspian sea , east to afghanistan and pakistan , south to the persian 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. between these ranges lies a high plateau where flowing water from the mountains disappears into desert sand . 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 climate , with the average temperature ranging between 10 and 35 c and the average relative humidity between 70 and 100% . the central plateau , situated between the alborz and zagros ranges , are very mountainous in the northwest where the ranges originate and becomes a lower desert in the east . the climate is dry with average temperature between 0 and 40 c , with hot , dry summers and cold , snow - laden winters . the persian gulf littoral and the khuzestan plain , to the south of the foothills of the zagros mountains , are characterized as a tropical climate ( fig . was conducted from 40 north to 29 south latitude and 44 western to 62 eastern longitude respectively ( 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 baluchestan , ( 12 ) kerman , ( 13 ) kohkiluyeh va boyer ahmad , ( 14 ) chahar mahal - o - bakhtiari , ( 15 ) esfahan , ( 16 ) markazi , and ( 17 ) khorasan - e - razavi ] ( fig . the sand flies were collected using 200 sticky paper traps ( consisted of white sheets 1521 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 34 villages suitable for sand fly breeding . these were located in the plains and mountains areas within each province and sampling was conducted twice each year ( june and september ) from 20082010 . sand flies were preserved in 96 % ethyl alcohol 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 . identification , sand flies were mounted in puri s medium , which was manufactured in our leishmaniasis laboratory in iran and identified after 24 h ( yaghoobi - ershadi and javadian 1997 ) using the keys of theodor and mesghali ( 1964 ) , artemiev ( 1973 ) , and lewis ( 1982 ) . morphological characters used for male identification of the subgenus adlerius included antennal segments , shape of the antennae and ascoid , number of hairs on coxite , length of style on ventral process , length of filament , length of pump , and sub - terminal tubercle of the aedeagus ( lewis 1982 ) . 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 measurements were conducted with an ocular micrometer and olympus microscope ( ch-2 ) ( fig . organs measured in morphometrical studies of the subgenus adlerius in iran statistical analyses were used to determine 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 - whitney test was used for comparing the results . sections of the abdomen , wings , and legs of the sand flies used for dna extraction were stored in 1.5 ml sterile microtubes containing 96 % ethanol . dna from 120 specimens were extracted , of which 50 were used for mtdna - pcr having 16 sequences ( table 1 ) . the washed , dried body of each sand fly was frozen then defrosted with liquid nitrogen and homogenised and ground with a heavy cylindrical metal and shaking vigorously . the content was washed with 300 ml lysis buffer [ 10 mm tris , 25 mm edta ( ph= 8) , 1 % sds , 25 mm nacl , 2% triton100 ] and boiled for 10 minutes . after the ethanol precipitation , dna was dissolved in 15 l distilled water and stored at 20 c . number of sand flies used in molecular studies a 25 l pcr reaction mixture consisted of 12.5 l premix 2[10 pcr buffer , mgcl , 10 mm dntps and dna taq polymerase ] , 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 5l of sand fly genomic 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 cycle ( 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 . in this study the chromatogram of the dna sequencing samples were edited by using bioedit software . the program clustal w ( thompson et al . our sequences were compared with the sequences 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 homologous between the female and the male specimens . we measured the eight different morphological 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 identified by sequence comparing alignment pairwise . therefore , if the more similarity score and less e - value were identified between the pairwise comparison , it showed that they are same species . iran is an arid land of 1.6 million km 

 , extending north to armenia , azerbaijan , turkmenistan , and the caspian sea , east to afghanistan and pakistan , south to the persian 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. between these ranges lies a high plateau where flowing water from the mountains disappears into desert sand . 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 climate , with the average temperature ranging between 10 and 35 c and the average relative humidity between 70 and 100% . the central plateau , situated between the alborz and zagros ranges , are very mountainous in the northwest where the ranges originate and becomes a lower desert in the east . the climate is dry with average temperature between 0 and 40 c , with hot , dry summers and cold , snow - laden winters . the persian gulf littoral and the khuzestan plain , to the south of the foothills of the zagros mountains , are characterized as a tropical climate ( fig . was conducted from 40 north to 29 south latitude and 44 western to 62 eastern longitude respectively ( 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 baluchestan , ( 12 ) kerman , ( 13 ) kohkiluyeh va boyer ahmad , ( 14 ) chahar mahal - o - bakhtiari , ( 15 ) esfahan , ( 16 ) markazi , and ( 17 ) khorasan - e - razavi ] ( fig . the sand flies were collected using 200 sticky paper traps ( consisted of white sheets 1521 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 34 villages suitable for sand fly breeding . these were located in the plains and mountains areas within each province and sampling was conducted twice each year ( june and september ) from 20082010 . sand flies were preserved in 96 % ethyl alcohol 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 . for species identification , sand flies were mounted in puri s medium , which was manufactured in our leishmaniasis laboratory in iran and identified after 24 h ( yaghoobi - ershadi and javadian 1997 ) using the keys of theodor and mesghali ( 1964 ) , artemiev ( 1973 ) , and lewis ( 1982 ) . morphological characters used for male identification of the subgenus adlerius included antennal segments , shape of the antennae and ascoid , number of hairs on coxite , length of style on ventral process , length of filament , length of pump , and sub - terminal tubercle of the aedeagus ( lewis 1982 ) . 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 measurements were conducted with an ocular micrometer and olympus microscope ( ch-2 ) ( fig . statistical analyses were used to determine 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 - whitney test was used for comparing the results . sections of the abdomen , wings , and legs of the sand flies used for dna extraction were stored in 1.5 ml sterile microtubes containing 96 % ethanol . dna from 120 specimens were extracted , of which 50 were used for mtdna - pcr having 16 sequences ( table 1 ) . the washed , dried body of each sand fly was frozen then defrosted with liquid nitrogen and homogenised and ground with a heavy cylindrical metal and shaking vigorously . the content was washed with 300 ml lysis buffer [ 10 mm tris , 25 mm edta ( ph= 8) , 1 % sds , 25 mm nacl , 2% triton100 ] and boiled for 10 minutes . after the ethanol precipitation , dna was dissolved in 15 l distilled water and stored at 20 c . number of sand flies used in molecular studies a 25 l pcr reaction mixture consisted of 12.5 l premix 2[10 pcr buffer , mgcl , 10 mm dntps and dna taq polymerase ] , 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 5l of sand fly genomic 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 cycle ( 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 . , in this study the chromatogram of the dna sequencing samples were edited by using bioedit software . our sequences were compared with the sequences 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 homologous between the female and the male specimens . we measured the eight different morphological 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 identified by sequence comparing alignment pairwise . therefore , if the more similarity score and less e - value were identified between the pairwise comparison , it showed that they are same 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 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 . 3 ) . number of males of subgenus adlerius species collected from 17 provinces we updated the morphological key for the male subgenus adlerius species of iran as follows : 
 one ascoid on antennal segments 915 . antenna 8 with two ascoids ..................... ( 2).-one ascoid on antennal segments 915 . antenna 8 with one ascoid . antenna 6 and 7 with one ascoid or with one long and one short ........................................................ ( 6).coxite with 1427 hairs in group , rarely 29 ............................................................ ( 3).-coxite with 29115 hairs in group , rarely 27 . part of 2785 hairs of group on distal half of coxite .......................................... tubercle of aedegus 616 m from tip ..................... p. turanicus ( akhoundi et al . 2011).-part of hair - group on distal half of coxite . tubercle of aedegus 1928 m from tip .............. p. brevis ( theodor and mesghali 1964).aedegus with rectangular subterminal notch ................ p. halepensis ( adler 1946).-aedegus with normal obtuse - angle subterminal notch ....................................................... sperm tubes are long ( 12001700 m ) ...................................................... p. longiductus ( theodor and mesghali 1964 ) ( also p.cf . ventral process of style long ( about 20 m ........ ( 7).-coxite with 3085 group - hairs antenna 6 and 7 with one ascoid or with one long and one short ........................................................ ( 6 ) . part of 2785 hairs of group on distal half of coxite .......................................... ( 4 ) . tubercle of aedegus 616 m from tip ..................... p. turanicus ( akhoundi et al . 2011 ) . tubercle of aedegus 1928 m from tip .............. p. brevis ( theodor and mesghali 1964 ) . aedegus with rectangular subterminal notch ................ p. halepensis ( adler 1946 ) . sperm tubes are long ( 12001700 m ) ...................................................... p. longiductus ( theodor and mesghali 1964 ) ( also p.cf . coxite very wide , whole group of 125200 hairs on its basal half ................................................ p. comatus ( zahraei et al . 2013 ) . -coxite narrow , part of hair - group on its distal half .................................... p. balcanicus ( theodor and mesghali 1964 ) ( also p.cf . the male and female sand flies of subgenus adlerius were collected from 12 provinces and the male were collected from 11 out of 17 provinces . the males were not collected 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 . 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 regions ( north - west , north - east , south - west and south - east ) . the p. halepensis specimens collected from three provinces ( ardabil , azerbaijan - e - sharqi and esfahan ) showed significant difference between the length of style and length of aedeagus . the statistical analysis of eight morphological characters of 30 males of p. balcanicus specimens confirmed that there is no significant difference between the characters in the specimens of ardabil and azerbaijane - 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. kabulensis 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 difference 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 ) the pcr experiments amplified a fragment of 550 bp of the mitochondrial genome and it was successfully achieved for 6 species of the subgenus adlerius ( table 1 , 2 ) . these were registered on jx885982 to jx 885998 accession numbers in the genbank . profile of sequenced species based on 
 cytb 

 gene in the subgenus 
 adlerius 
  in the pair wise comparison of the specimens , we focused on the genetic similarities of the specimens among each other and determined the level of such similarities of male and female specimens and identified the females among 6 male species . the pair wise comparison in azv8 female 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 sequences which show 8 ( 1.5% ) polymorphic sites and it is confirmed that they are 2 haplotypes ( table 3 ) . comparison of 
 cytb gene nucleotide sequences in the subgenus 
 adlerius 

 species 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% similarities among the male and the shc11 female specimen . it shows 22 ( 4.3% ) polymorphic sites and it confirms that there are 4 haplotypes ( 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 polymorphic sites is 46 ( 10% ) and the sequence comparison revealed that there are 6 haplotypes . also for the two males of p. brevis , it shows 99% similarities with 2 ( 0.4% ) polymorphic sites and 2 haplotypes between these two sequenced specimens ( table 3 ) . the sequence result of the length 550 bp of cytb gene shows that 5 male species : p. halepensis , p. brevis , p. kabulensis , p. comatus and p. balcanicus are in the areas of our study . longiductus which are close to p. balcanicus and p. longiductus are on neighbor - joining tree ( fig . the specimens 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 similarities among the males and females confirm our above results ( fig . neighbor - joining tree for dna haplotypes of cytb - mtdna of the subgenus adlerius sand flies 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 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 . 3 ) . number of males of subgenus adlerius species collected from 17 provinces we updated the morphological key for the male subgenus adlerius species of iran as follows : 
 one ascoid on antennal segments 915 . antenna 8 with two ascoids ..................... ( 2).-one ascoid on antennal segments 915 . antenna 8 with one ascoid . antenna 6 and 7 with one ascoid or with one long and one short ........................................................ ( 6).coxite with 1427 hairs in group , rarely 29 ............................................................ ( 3).-coxite with 29115 hairs in group , rarely 27 . part of 2785 hairs of group on distal half of coxite .......................................... tubercle of aedegus 616 m from tip ..................... p. turanicus ( akhoundi et al . 2011).-part of hair - group on distal half of coxite . tubercle of aedegus 1928 m from tip .............. p. brevis ( theodor and mesghali 1964).aedegus with rectangular subterminal notch ................ p. halepensis ( adler 1946).-aedegus with normal obtuse - angle subterminal notch ....................................................... sperm tubes are long ( 12001700 m ) ...................................................... p. longiductus ( theodor and mesghali 1964 ) ( also p.cf . ventral process of style long ( about 20 m ........ ( 7).-coxite with 3085 group - hairs antenna 6 and 7 with one ascoid or with one long and one short ........................................................ ( 6 ) . part of 2785 hairs of group on distal half of coxite .......................................... ( 4 ) . tubercle of aedegus 616 m from tip ..................... p. turanicus ( akhoundi et al . 2011 ) . tubercle of aedegus 1928 m from tip .............. p. brevis ( theodor and mesghali 1964 ) . aedegus with rectangular subterminal notch ................ p. halepensis ( adler 1946 ) . sperm tubes are long ( 12001700 m ) ...................................................... p. longiductus ( theodor and mesghali 1964 ) ( also p.cf . coxite very wide , whole group of 125200 hairs on its basal half ................................................ p. comatus ( zahraei et al . 2013 ) . -coxite narrow , part of hair - group on its distal half .................................... p. balcanicus ( theodor and mesghali 1964 ) ( also p.cf . the male and female sand flies of subgenus adlerius were collected from 12 provinces and the male were collected from 11 out of 17 provinces . the males were not collected 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 . 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 regions ( north - west , north - east , south - west and south - east ) . the p. halepensis specimens collected from three provinces ( ardabil , azerbaijan - e - sharqi and esfahan ) showed significant difference between the length of style and length of aedeagus . the statistical analysis of eight morphological characters of 30 males of p. balcanicus specimens confirmed that there is no significant difference between the characters in the specimens of ardabil and azerbaijane - 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. kabulensis 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 difference 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 ) the pcr experiments amplified a fragment of 550 bp of the mitochondrial genome and it was successfully achieved for 6 species of the subgenus adlerius ( table 1 , 2 ) . these were registered on jx885982 to jx 885998 accession numbers in the genbank . profile of sequenced species based on 
 cytb 

 gene in the subgenus 
 adlerius 
  in the pair wise comparison of the specimens , we focused on the genetic similarities of the specimens among each other and determined the level of such similarities of male and female specimens and identified the females among 6 male species . the pair wise comparison in azv8 female 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 sequences which show 8 ( 1.5% ) polymorphic sites and it is confirmed that they are 2 haplotypes ( table 3 ) . comparison of 
 cytb gene nucleotide sequences in the subgenus 
 adlerius 

 species 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% similarities among the male and the shc11 female specimen . it shows 22 ( 4.3% ) polymorphic sites and it confirms that there are 4 haplotypes ( 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 polymorphic sites is 46 ( 10% ) and the sequence comparison revealed that there are 6 haplotypes . also for the two males of p. brevis , it shows 99% similarities with 2 ( 0.4% ) polymorphic sites and 2 haplotypes between these two sequenced specimens ( table 3 ) . the sequence result of the length 550 bp of cytb gene shows that 5 male species : p. halepensis , p. brevis , p. kabulensis , p. comatus and p. balcanicus are in the areas of our study . longiductus which are close to p. balcanicus and p. longiductus are on neighbor - joining tree ( fig . the specimens 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 similarities among the males and females confirm our above results ( fig . neighbor - joining tree for dna haplotypes of cytb - mtdna of the subgenus adlerius sand flies species but , the species involved have seldom been identified because despite the taxonomic efforts of artemiev ( 1980 ) and others . the morphological characteristics do not distinguish between the females of adlerius species that are generally recognized on the basis of male morphology ( lewis 1982 , sadlova et al . 2013 ) and p. kabulensis during the present study has assisted to revise the morphological key . the access to new molecular techniques and genetic information about these probable vectors ( artemiev 1973 , ardehali et al . 1995 ) and their systematics , may aid in resolving the problem of identification 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 characteristics is highly variable . the latest revision of the genus phlebotomus sand flies in old world is included : ( 1 ) adlerius , ( 2 ) anaphlebotomus , ( 3 ) euphlebotomus , ( 4 ) kasaulius , ( 5 ) transphlebotomus , ( 6 ) synphlebotomus , ( 7 ) phlebotomus , ( 8) larroussius , ( 9 ) paraphlebotomus ( rispail and leger 1998 ) . 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 monophyletic of the subgenus adlerius and the morphological identification of these six species are confirmed by the reliable morphological identification keys . the haplotypes are 6 in this species . according to our morphological , 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 different areas . in the statistical analysis , the length of the style in p. halepensis specimens in ardabil province was significantly different from that of esfahan and azerbaijan - e - sharqi . azv2 is in a different branch but is closer to msh458 and arm754 specimens . 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 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 . 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 % similarity 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 different from the specimens of other geographical regions . however , ckhe1342-female is in different branch because of the geographical 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 separated from another species of adlerius group many years ago than the other adlerius species . according to the statistical analysis and molecular findings the present study on molecular analysis of p. brevis population shows that one line - age with two haplotypes for cytb . 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 overlapping . 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 molecular 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 % similarities among themselves . also , it is confirmed that there are 4 haplotypes among the four sequenced specimens . in the present study , azv4 and azv5 however , in our upgma tree , we found that they are differents between p. balcanicus specimens and they are in an independent branch but closer to p. balcanicus branches . the results of the molecular investigations indicate that according to the concept of the species they ( azv5 and azv4 ) can not be considered as a separate species . because the concept of a species is based on the reproductive isolation , while the existence of the intermediate forms and the same mitochondrial gene sequences in this study indicates that there is no reproductive isolation between azv5 morphotypes and azv4 morphotypes . also , we can not assume these two specimens as subspecies because two subspecies can not be found as sympatric . so , this indicates that these species e.g p. longiductus and p. balcanicus in iran , including their populations have morphological and genetic differences . the similarity between azsh10 and azv11 is 97% and e - value between them is 0.0 . the branches of the four haplotypes of p. balcanicus show that genetically they are not far . the similarity and e - value between the shc11 and azv5 are 99 % and 0.0 respectively . the pairwise comparison of the nucleotides of cytb - mtdna gene in azv8 female specimen with nucleotides of six male species shows that the maximum similarity i.e 98% is found in azv10 . so , these evidences confirm our results as both azv8 and azv10 are p. comatus . in the upgma tree , 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 species 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 separating one adlerius species from the other . at times it is difficult for the researchers as the numbers and locations of the hairs of different species are the same and overlap among themselves . in p. kabulensis , the hairs on coxite are 2750 and in p. longiductus the hairs on coxite are 5085 . if we collect a specimen which has 50 hairs on its coxite , we can not identify that it belongs to which one of the afore noted two species . also for p. balcanicus and p. longiductus , we encountered the same difficulty . in this case , the researchers use it is mainly used in academic writings to indicate a reference to a contrasting finding or viewpoint ( it is often used instead of  but see  ) . it can also appear occasionally in binomial nomenclature by placing before the species name to indicate that the species is not confirmed ( strunk and white 1979 ) . balcanicus . 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 . obviously , detailed studies for detection of genetic markers to identify the females of six adlerius species are essential . also , further studies on morphology and genetic diversity of the populations of the subgenus adlerius will be of particular importance . the positions of three female and males of subgenus adlerius in the upgma tree are correct and the similarities among them confirm our results . the branches of each species are not genetically distinct that justify the overlapping morphological characters among them . molecular sequencing of cytb - mtdna haplotypes can be used for female identification of different species of subgenus adlerius in iran .