 the phlebotomine sand fly phlebotomus sergenti s.l . parrot , 1917 originally described from algeria in 1917 , has a broad range of distribution which covers areas of the southern mediterranean ( morocco , algeria , tunisia ) , the western mediterranean 
( portugal , spain , sicily ) , middle east , arabia , iran , afghanistan , pakistan , and northern parts of india .
flies of this species have been incriminated as the main vector of cutaneous leishmaniasis due to leishmania tropica ( clt ) throughout iran and other clt foci in the world ( nadim et al .
the presence of this sand fly in l. tropica free areas and the differences in the transmission patterns of the clt could be related to the existence of intraspecific variability or even cryptic vector species ( depaquit et al . 2002 , yahyia et al . 2004 ) .
intraspecific variability has been shown in some morphological and molecular characteristics of various populations from different countries .
( 2002 ) studied the intraspecific variability of the internal transcribed spacer 2 ( its2 ) of 12 populations of p. sergenti s.l . from ten different countries
accordingly , two branches could be identified : one was related to the northeastern mediterranean area ( cyprus , syria , and turkey ) and pakistan , and the second related to southwestern of the first one including north africa , egypt , and morocco .
these branches are in accordance with postulated migration routes of p. sergenti s.l . along the tethys sea during the miocene era .
the sand flies belonging to these two different branches seem to differ in ecology , host preferences , and possibly also in vectorial capacity ( depaquit et al .
were also studied on 28 iranian populations and a few samples from greece , morocco , lebanon , turkey , pakistan , and syria ( moin - vaziri et al .
according to this study , based on the number of setae and the width of basal lobe of coxite , three morphotypes were identified as a , b and c , with some intermediate forms .
morphotype a was considered as p. sergenti sergenti , morphotype b , was identified as p. sergenti similis , and morphotype c had an elongated style in comparison with p. sergenti sergenti . in another molecular study on p. sergenti s.l .
populations of spain and morocco using rdna its2 and mtdna cytb sequences ( baron et al . 2008 ) , a high genetic diversity including five ribosomal and 16 mitochondrial haplotypes was found within 25 specimens .
based on these studies , it is judicious to consider the potential existence of sibling species within this taxon .
if sibling species within p. sergenti s.l . were proven , it would have important implications in epidemiology as well as in experimental studies . however , having found several intermediate morphotypes as well as sympatric ecological niche of morphotypes a , b and c , postulated us to test the gene flow between the sympatric morphotypes of this taxon using sequence analysis of the its2 region of rdna gene .
this multicopy gene involves homogeneization processes usually called molecular drive ( dover 1982 ) and has provided resolution in several studies at the taxonomic level for the larroussius and paraphlebotomus subgenera ( dover 1982 , depaquit et al .
the use of its2 region has many advantages including , high and low mutation rhythm at interspecies and intraspecies level respectively , speed and ease of use , multiple target sites , predefined marker systems , known pcr primers , pre - existing knowledge of them in some sand fly species ( di muccio 2000 , depaquit et al .
this study was conducted to verify the molecular variation of iranian p. sergenti populations , and to compare them with available data in genbank .
the geographical locations from which p. sergenti s.l . was sampled are shown in fig .
specimens caught by sticky papers were washed in a bath of acetone before being stored .
only male specimens were selected for the study of morphological and molecular variability since their characters tends to be more reliable .
( 1 ) tehran , ( 2 ) mashhad , ( 3 ) neishabur , ( 4 ) izeh , ( 5 ) bushehr , ( 6 ) bandar - e - abbas , ( 7 ) bam , ( 8) iranshahr .
dark spots represent clt endemic foci in iran the head and genitalia of individual male sand flies were cut off within a drop of ethanol , cleared in boiling marc - andr solution , and mounted between slide and cover slide in berlese fluid and morphologically identified based on external and internal characters of the head and genitalia according to the known identification keys ( theodor and mesghali 1964 ) .
morphometric measurements were performed to determine morphotypes of specimens as explained by moin - vaziri et al .
the body related to the specimen was stored dried in a vial at 20 c before dna extraction .
based on ecological conditions of the location where specimens were collected and on the morphological differences noted by means of morphometric analysis , a few specimens from each morphotype ( moin - vaziri et al .
genomic dna was extracted from the thorax , wings , legs and abdomen of either individual sand flies using the qiamp dna mini kit ( qiagen , germany ) ( depaquit et al . 2002 ) or ish - horowicz with small modification as described by rassi et al .
pcr was used to amplify a fragment of 480516 bp containing the its2 of sand fly rdna ( depaquit et al .
pcr were performed in a 50 l volume using 5 l of extracted dna solution and 50 pmol of each of the two primers of c1a : 5-cct ggt tag ttt ctt ttc ctc cgc t-3 and jts3 : 5-cgc agc taa ctg tgt gaa atc-3. the pcr mix contained ( final concentrations ) 10 mm tris hcl , ph 8.3 , 1.5 mm mgcl2 , kcl 50 mm , triton x 100 0.01% , 200 m dntp each , and 0.25 l ( 1.25 units ) of taq dna polymerase ( eurobio ) .
initial denaturation at 94 c for 5 min was followed by 35 cycles of denaturation at 94 c for 30 sec , annealing at 62 c for 1 min . and extension at 72 c for 1 min with a final elongation time of 10 min at 72 c .
purification of the pcr product was made by agarose - gel fractionation , using the perfect prep gel cleanup ( eppendorf , germany ) .
direct sequencing of both dna strands was performed by qiagen ( hilden , germany ) and the department of parasitology ( university of valencia , spain ) using the primers used for dna amplification .
sequences were edited and aligned to identify haplotypes ( = unique sequences ) by means of the clustalw software package ( www.ebi.ac.uk/clustalw ) and manually adjusted , if necessary .
they were analyzed using the neighbor - joining ( nj ) method provided in the must software package ( philippe 1993 ) .
available sequences of p. sergenti and p. similis were retrieved from genbank and used for phylogenetic analysis ( table 1 ) .
the geographical locations from which p. sergenti s.l . was sampled are shown in fig .
specimens caught by sticky papers were washed in a bath of acetone before being stored .
only male specimens were selected for the study of morphological and molecular variability since their characters tends to be more reliable .
( 1 ) tehran , ( 2 ) mashhad , ( 3 ) neishabur , ( 4 ) izeh , ( 5 ) bushehr , ( 6 ) bandar - e - abbas , ( 7 ) bam , ( 8) iranshahr .
dark spots represent clt endemic foci in iran the head and genitalia of individual male sand flies were cut off within a drop of ethanol , cleared in boiling marc - andr solution , and mounted between slide and cover slide in berlese fluid and morphologically identified based on external and internal characters of the head and genitalia according to the known identification keys ( theodor and mesghali 1964 ) .
morphometric measurements were performed to determine morphotypes of specimens as explained by moin - vaziri et al .
the body related to the specimen was stored dried in a vial at 20 c before dna extraction .
based on ecological conditions of the location where specimens were collected and on the morphological differences noted by means of morphometric analysis , a few specimens from each morphotype ( moin - vaziri et al .
genomic dna was extracted from the thorax , wings , legs and abdomen of either individual sand flies using the qiamp dna mini kit ( qiagen , germany ) ( depaquit et al . 2002 ) or ish - horowicz with small modification as described by rassi et al .
pcr was used to amplify a fragment of 480516 bp containing the its2 of sand fly rdna ( depaquit et al .
pcr were performed in a 50 l volume using 5 l of extracted dna solution and 50 pmol of each of the two primers of c1a : 5-cct ggt tag ttt ctt ttc ctc cgc t-3 and jts3 : 5-cgc agc taa ctg tgt gaa atc-3. the pcr mix contained ( final concentrations ) 10 mm tris hcl , ph 8.3 , 1.5 mm mgcl2 , kcl 50 mm , triton x 100 0.01% , 200 m dntp each , and 0.25 l ( 1.25 units ) of taq dna polymerase ( eurobio ) .
initial denaturation at 94 c for 5 min was followed by 35 cycles of denaturation at 94 c for 30 sec , annealing at 62 c for 1 min . and extension at 72 c for 1 min with a final elongation time of 10 min at 72 c .
purification of the pcr product was made by agarose - gel fractionation , using the perfect prep gel cleanup ( eppendorf , germany ) .
direct sequencing of both dna strands was performed by qiagen ( hilden , germany ) and the department of parasitology ( university of valencia , spain ) using the primers used for dna amplification .
sequences were edited and aligned to identify haplotypes ( = unique sequences ) by means of the clustalw software package ( www.ebi.ac.uk/clustalw ) and manually adjusted , if necessary .
they were analyzed using the neighbor - joining ( nj ) method provided in the must software package ( philippe 1993 ) .
available sequences of p. sergenti and p. similis were retrieved from genbank and used for phylogenetic analysis ( table 1 ) .
morphological investigation on specimens of 28 populations of p. sergenti sl collected from 11 provinces of iran revealed three main morphotypes ( a , b and c ) and a few intermediate forms ( moin - vaziri et al .
morphotype b was related to p. cf similis , according to perfiliev ( 1968 ) and morphotype c represents specimens with a curved basal lobe of coxite without a style as globulous as that of p. sergenti .
in addition to the three main morphotypes , some intermediate forms were identified among the collected samples .
these morphotypes were found sympatric in several provinces of the country . however , the proportion of each morphotype varied within each region .
the size of rdna - its2 fragments amplified was about 480 bp , from which 278 nucleotides were attributed to its2 , 130 bp to 5.8s and 72 bp to 28s genes .
however , in addition to the main pcr product , there was an additional product , close to the main amplicon , in all amplification .
despite of gel purification before sequencing , these unwanted products affected the results of sequencing . by the way
, we could obtain a total of 249253 bp length , including most part ( 218222 bp ) of its2 and 31 bp of 28s , from 11 specimens comprising 5 , 4 , and 2 individuals of morphotype a , b and intermediate form a / c respectively .
the sequences were submitted to genbank ( accession numbers : ef434818-ef434828 ) . comparative its2 sequence analysis of 11 p. sergenti s.l .
individuals showed 13 ( 5.2% ) polymorphic sites , from which 46% was due to indel ( insertion / deletion ) and 54% was due to substitutions .
all polymorphic sites were located in its2 region and 28s region was identical in all specimens .
there were two variable microsatellite regions in the its2 aligned part of the specimens , which showed variable repeats of a poly ( at ) microsatellite .
in the first variable microsatellite sites , two specimens contain 9 , three specimens had 8 , and six other specimens comprised 7 repeats of at ( fig .
other entries from genbank also showed variable number of at repeats in this site . in the second polymorphic microsatellite ,
variations between the iranian specimens were due to two single nucleotide insertions / deletions ( indels ) and one variable repeats of a poly ( at ) microsatellite .
lower degrees of variation in the poly ( at ) microsatellite also have been observed in non - iranian entries ( fig .
2 ) . microsatellite region in its2-rdna sequences of different populations of phlebotomus sergenti s.l . from iran and other parts of world retrieved from genbank ( marked with * ) .
a : morphotype a , b : morphotype b , a / c : intermediate forms of morphotype a and c , ser , p. sergenti , sim , p. similis based on the sequence alignment , nine different haplotypes ( nominated as haplotype i  ix ) have been identified , three haplotypes ( vii  ix ) in the south west ( sw ) , two haplotypes ( i and vi ) in the north east ( ne ) and 5 haplotypes ( i  v ) in the north - center - south ( ncs ) populations ( table 1 , 2 , fig .
3 ) . a sample of haplotype i was positioned in ne lineage . these haplotypes differed in 111 nucleotide positions .
the its2 ribosomal haplotype i with three repeats was the most frequently haplotype among all haplotypes .
neighbor - joining tree inferred from 253 bp of its2-rdna sequences including 222bp of its2 and 31 bp of 28s of phlebotomus sergenti s.l .
populations from iran . a , morphotype a , b , morphotype b , a / c , intermediate form of morphotype a and c , ser , p. sergenti , sim , p. similis comparison of nucleotide characters of phlebotomus sergenti of iran at polymorphic sites of the rdna its2 . gap ( indel ) shows by  distance analysis of the its2 sequences indicated three main lineages ( fig .
these lineages were so - called sw , ne , ncs lineages which included haplotypes of south - west , north - east and a mix of haplotypes extended from northwest to the central to the southeast of the country , respectively .
when we added the its2 sequence data of other p. sergenti and p. similis populations as representatives of other parts of the world for phylogenetic analysis , except for one iranian sample which was close to the european samples , all of the iranian haplotypes were associated with the north - eastern mediterranean populations including turkey , cyprus , syria , and also pakistan ( fig .
the p. similis populations from iran were associated with other p. sergenti populations , however , the european p. similis populations clustered separately .
this might be due to independent accumulated mutations in distinct geographical populations resulted in separate branch in the tree .
the phylogenetic tree deduced from combination of the its2 and the 381 bp of mtdna cytb - nadh1 sequences obtained in our previous study ( moin - vaziri et al .
2007 ) revealed similar tree and did not resolve p. segenti from p. similis ( data are not shown ) .
neighbor - joining tree inferred from 253 bp of its2-rdna sequences including 222 bp of its2 and 31 bp of 28s of phlebotomus sergenti s.l .
populations from iran and other available data from genbank originated from europe , africa , and asia .
a , morphotype a , b , morphotype b , a / c , intermediate form of morphotype a and c , ser , p. sergenti , sim , p. similis , p. papatasi ( an : ef408801 ) has been used as an outgroup
morphological investigation on specimens of 28 populations of p. sergenti sl collected from 11 provinces of iran revealed three main morphotypes ( a , b and c ) and a few intermediate forms ( moin - vaziri et al .
morphotype b was related to p. cf similis , according to perfiliev ( 1968 ) and morphotype c represents specimens with a curved basal lobe of coxite without a style as globulous as that of p. sergenti .
in addition to the three main morphotypes , some intermediate forms were identified among the collected samples .
these morphotypes were found sympatric in several provinces of the country . however , the proportion of each morphotype varied within each region .
the size of rdna - its2 fragments amplified was about 480 bp , from which 278 nucleotides were attributed to its2 , 130 bp to 5.8s and 72 bp to 28s genes .
however , in addition to the main pcr product , there was an additional product , close to the main amplicon , in all amplification .
despite of gel purification before sequencing , these unwanted products affected the results of sequencing . by the way , we could obtain a total of 249253 bp length , including most part ( 218222 bp ) of its2 and 31 bp of 28s , from 11 specimens comprising 5 , 4 , and 2 individuals of morphotype a , b and intermediate form a / c respectively .
individuals showed 13 ( 5.2% ) polymorphic sites , from which 46% was due to indel ( insertion / deletion ) and 54% was due to substitutions .
all polymorphic sites were located in its2 region and 28s region was identical in all specimens .
there were two variable microsatellite regions in the its2 aligned part of the specimens , which showed variable repeats of a poly ( at ) microsatellite .
in the first variable microsatellite sites , two specimens contain 9 , three specimens had 8 , and six other specimens comprised 7 repeats of at ( fig .
2 ) . other entries from genbank also showed variable number of at repeats in this site . in the second polymorphic microsatellite
, variations between the iranian specimens were due to two single nucleotide insertions / deletions ( indels ) and one variable repeats of a poly ( at ) microsatellite .
lower degrees of variation in the poly ( at ) microsatellite also have been observed in non - iranian entries ( fig .
2 ) . microsatellite region in its2-rdna sequences of different populations of phlebotomus sergenti s.l . from iran and other parts of world retrieved from genbank ( marked with * ) .
a : morphotype a , b : morphotype b , a / c : intermediate forms of morphotype a and c , ser , p. sergenti , sim , p. similis based on the sequence alignment , nine different haplotypes ( nominated as haplotype i  ix ) have been identified , three haplotypes ( vii  ix ) in the south west ( sw ) , two haplotypes ( i and vi ) in the north east ( ne ) and 5 haplotypes ( i  v ) in the north - center - south ( ncs ) populations ( table 1 , 2 , fig .
3 ) . a sample of haplotype i was positioned in ne lineage . these haplotypes differed in 111 nucleotide positions .
the its2 ribosomal haplotype i with three repeats was the most frequently haplotype among all haplotypes .
neighbor - joining tree inferred from 253 bp of its2-rdna sequences including 222bp of its2 and 31 bp of 28s of phlebotomus sergenti s.l .
populations from iran . a , morphotype a , b , morphotype b , a / c , intermediate form of morphotype a and c , ser , p. sergenti , sim , p. similis comparison of nucleotide characters of phlebotomus sergenti of iran at polymorphic sites of the rdna its2 .
gap ( indel ) shows by  distance analysis of the its2 sequences indicated three main lineages ( fig .
these lineages were so - called sw , ne , ncs lineages which included haplotypes of south - west , north - east and a mix of haplotypes extended from northwest to the central to the southeast of the country , respectively .
when we added the its2 sequence data of other p. sergenti and p. similis populations as representatives of other parts of the world for phylogenetic analysis , except for one iranian sample which was close to the european samples , all of the iranian haplotypes were associated with the north - eastern mediterranean populations including turkey , cyprus , syria , and also pakistan ( fig .
the p. similis populations from iran were associated with other p. sergenti populations , however , the european p. similis populations clustered separately .
this might be due to independent accumulated mutations in distinct geographical populations resulted in separate branch in the tree .
the phylogenetic tree deduced from combination of the its2 and the 381 bp of mtdna cytb - nadh1 sequences obtained in our previous study ( moin - vaziri et al .
2007 ) revealed similar tree and did not resolve p. segenti from p. similis ( data are not shown ) .
neighbor - joining tree inferred from 253 bp of its2-rdna sequences including 222 bp of its2 and 31 bp of 28s of phlebotomus sergenti s.l .
populations from iran and other available data from genbank originated from europe , africa , and asia .
a , morphotype a , b , morphotype b , a / c , intermediate form of morphotype a and c , ser , p. sergenti , sim , p. similis , p. papatasi ( an : ef408801 ) has been used as an outgroup
results of this study showed a high diversity between specimens of p. sergenti s.l . in iran where we found nine ribosomal haplotypes .
however , the rate of genetic variation ( 3% ) within the its2 locus between the iranian populations is half of the rate ( 6% ) that was previously observed in the mitochondrial ( cytb - nadh1 ) sequences ( moin - vaziri et al .
it seems that mtdna sequences are more appropriate for the study of the intraspecific variability of p. sergenti s.l . in a more limited geographical environment .
the high diversity observed in its2 is in agreement with previous studies indicating different ecological , morphological , and molecular variation among iranian p. sergenti populations .
previous studies revealed at least three morphotypes ( a , b and c ) . however , in this study like previous study , we have not found any correlation between genotypes , ecotypes , or morphotypes and the results obtained here do not support the presence of sibling species ( p. sergenti and p. similis ) within the taxon .
however , results of this study are in conflict with a phylogenetic analysis of nuclear ribosomal dna of depaquit et al .
( 2002 ) that showed the monophyly both of p. sergenti s.l . and p. similis and they were not sister species .
their result confirmed a study previously carried out , using morphological and morphometric approaches for examining the status of the two species by the same investigators ( depaquit et al .
their worldwide attempt showed allopatric situation of the two taxa at that time . according to our findings
, it seems that this group ( morphotypes ) can not be considered as two different species because firstly we have found many intermediate morphological and genetically forms among the specimens , which indicate possible gene flow and lack of reproductive isolation between them .
secondly , finding identical genotypes ( 100% ) among different morphotypes that are p. sergenti and that p. similis in mtdna ( moin - vaziri et al .
the molecular drive characteristic of its2-rdna strongly indicated that these two morphotypes are not isolated reproductively .
moreover having found these two taxa sympatric in most provinces of iran differs from taxonomic conception of subspecies . according to our data , we encountered with different morphologically populations of p. sergenti s.l . in iran . in light of our results , it is too early to come to a final decision on taxonomic status of the species .
more molecular , morphological , and hybridization studies between the two taxa , particularly between geographically distinct populations is necessary . in iran , clt is endemic in 14 foci located in 8 provinces restricted to large and medium sized cities in different parts of the country ( yaghoobi - ershadi 2012).we identified three its2 lineages of p. sergenti which is in agreement to the ones identified using the mtdna loci ( moin - vaziri et al .
lineage ne and ncs correspond to the main foci of clt in the country which included both typical and intermediate morphotypes of p. cf sergenti and p. cf similis .
these findings warrant studies to examine if clt is due to differences in the vectorial capacity of the p. sergenti s.l .
lineages ( ne / ncs versus sw ) or other ecological and epidemiological factors are involved .
phylogenetic sequence analysis revealed that most iranian haplotypes were associated with the northeastern mediterranean populations .
similar to the sequences of mtdna coi gene , its2 sequences could not resolve p. sergenti from p. similis and did not support the possible existence of sibling species or subspecies within p. sergenti s.l .. more molecular studies on other genes or hybridization should be done to clarify the status of different morphotypes of mentioned species in iran .
moreover , with regard to the importance of the epidemiology of leishmaniasis , further studies need to be performed on the possible role of these three morphotypes in the transmission of l. tropica .
 background : phlebotomus sergenti s.l . is considered the most likely vector of leishmania tropica in iran . 
 although two morphotypes- p. sergenti sergenti ( a ) and p. sergenti similis ( b)-have been formally described , further morphological and a molecular analysis of mitochondrial cytochrome oxidase i ( mtdna - coi ) gene revealed inconsistencies and suggests that the variation between the morphotypes is intraspecific and the morphotypes might be identical species.methods:we examined the sequence of the its2-rdna of iranian specimens of p. sergenti s.l . , comprising p. cf sergenti , p. cf similis , and intermediate morphotypes , together with available data in genbank.results:sequence analysis showed 5.2% variation among p. sergenti s.l . 
 morphotypes . 
 almost half of the variation was due to the number of an at microsatellite repeats in the center of the spacer . 
 nine haplotypes were found in the species constructing three main lineages corresponding to the origin of the colonies located in southwest ( sw ) , northeast ( ne ) , and northwest - center - southeast ( ncs ) . 
 lineages ncs and ne included both typical p. cf sergenti and p. cf similis and intermediate morphotypes.conclusion:phylogenetic sequence analysis revealed that , except for one iranian sample , which was close to the european samples , other iranian haplotypes were associated with the northeastern mediterranean populations including turkey , cyprus , syria , and pakistan . 
 similar to the sequences of mtdna coi gene , its2 sequences could not resolve p. sergenti from p. similis and did not support the possible existence of sibling species or subspecies within p. sergenti s.l .. 