Sand fly fauna of Crete and the description of Phlebotomus (Adlerius) creticus n. sp. (Diptera: Psychodidae)
Parasites & Vectors volume 13, Article number: 547 (2020)
The Greek island of Crete is endemic for both visceral leishmaniasis (VL) and recently increasing cutaneous leishmaniasis (CL). This study summarizes published data on the sand fly fauna of Crete, the results of new sand fly samplings and the description of a new sand fly species.
All published and recent samplings were carried out using CDC light traps, sticky traps or mouth aspirators. The specific status of Phlebotomus (Adlerius) creticus n. sp., was assessed by morphological analysis, cytochrome b (cytb) sequencing and MALDI-TOF protein profiling.
Published data revealed the presence of 10 Phlebotomus spp. and 2 Sergentomyia spp. During presented field work, 608 specimens of 8 species of Phlebotomus and one species of Sergentomyia were collected. Both published data and present samplings revealed that the two most common and abundant species were Phlebotomus neglectus, a proven vector of Leishmania infantum causing VL, and Ph. similis, a suspected vector of L. tropica causing CL. In addition, the field surveys revealed the presence of a new species, Ph. (Adlerius) creticus n. sp.
The identification of the newly described species is based on both molecular and morphological criteria, showing distinct characters of the male genitalia that differentiate it from related species of the subgenus Adlerius as well as species-specific sequence of cytb and protein spectra generated by MALDI-TOF mass spectrometry.
Phlebotomine sand flies (Diptera: Psychodidae) are hematophagous insects that transmit the protozoan parasites Leishmania spp. as well as the bacterium Bartonella bacilliformis and viruses (phleboviruses, vesiculoviruses and orbiviruses) . Four medically important Leishmania species circulate in various regions of the Mediterranean basin: L. infantum causing both zoonotic visceral (VL) and anthroponotic cutaneous (CL) leishmaniasis in humans and canine leishmaniasis in dogs (CanL); L. major causing zoonotic CL; L. tropica causing anthroponotic or zoonotic CL [2,3,4] and the recently introduced L. donovani causing both anthroponotic VL and CL .
Crete is endemic for VL and CanL caused by L. infantum, with an increasing number of CL every year caused by L. tropica whilst the danger of the introduction of new species/zymodemes is enhanced by the arrival of refugees from endemic areas . The key factor for determining the distribution and spread of leishmaniases is by recording the geographical distribution and abundance of the medically important sand fly vectors. Generally, in Greece, proven or suspected vectors are Ph. similis Artemiev & Neronov, 1984 for Leishmania tropica and Ph. perfiliewi Parrot, 1930, Ph. tobbi Adler & Theodor, 1930 and Ph. neglectus Tonnoir, 1921 for L. infantum .
The main aims of this study were to review historical data on the sand fly fauna of Crete and describe a new Phlebotomus (Adlerius) species sampled during surveys carried out from 2014 to 2019. To confirm the status of the new species, in addition to morphological and genetic criteria, for the first time, protein profiling using MALDI-TOF mass spectrometry was also deployed.
Sand fly data
All published data recording sand fly presence in Crete were gathered. The literature search concerning the review follows the Prisma Journal Publishing protocol workflow . PubMed, Web of Science, Google Scholar databases and web searches were screened from 1910 to 30 November 2019, using the following keywords: Phlebotominae, Phlebotomus, Crete, Sergentomyia, Greece, sand flies, CDC light traps, sticky traps, mouth aspirator, electric aspirator. Full text articles in English containing information on phlebotomine sand flies from Crete were selected. Other articles, including those in other languages, that contain desired information were also included, based on the cited databases knowledge of the authors. Samplings were carried out using CDC miniature light traps (John W. Hock Co., Gainesville, FL, USA), Sticky Traps (A4 paper coated with castor oil) or mouth aspirators. The sampling sites were mostly animal farms, houses, schools, churches, deserted houses, wells, caves and rural areas with different cultivations or wild vegetation.
Study areas and sampling
The new samplings were carried out in five study areas, two in Heraklion, two in Chania and one in Lasithi between 2014 and 2019 (Table 1). The study areas were Fodele and Foinikia in Heraklion prefecture, Agia Roumeli and Botanical Gardens in Chania and Xerokampos in Lasithi. The CDC miniature light traps equipped with a fine net cage were used for all the samplings. In all sampling areas, 7–9 light traps were placed in different microhabitats for one to two days. The light traps were set 1.5 h before sunset and were collected 2 h after dawn. Specifically, in Fodele 9 light traps for 2 sampling nights in May 2019 were used, in Foinikia 9 traps for 1 night in August of 2018, in Agia Roumeli 7 traps for 2 sampling nights in May of 2014, and in Botanical Garden 9 traps for 2 sampling nights in August of 2019. In Xerokampos two samplings were carried out during April of 2014 and in May of 2019 using 9 traps for 2 sampling nights.
Sand fly morphological identification
The specimens were stored in vials containing absolute ethanol HPLC grade (Thermo Fischer Scientific, Gloucester, UK) at room temperature prior to mounting, except for the specimens later analyzed by MALDI-TOF protein profiling which were stored in 70% molecular grade ethanol at − 20 °C before further processing. The samples were divided into specimens mounted in toto and specimens processed for molecular biology. In the latter case, the head, thorax including wings and genitalia were removed and placed in a drop of ethanol before their processing similar to specimens mounted in toto, while the other parts of the body were kept in ethanol for molecular analysis. Soft tissues were lysed in a bath of KOH 10%, then bleached in Marc-André solution, and mounted between microscope slide and cover slide in Euparal® for species identification after dehydration in graded ethanol series . Some specimens were mounted immediately after clearing in Marc-André solution or in high viscosity CMCP-10 medium (Polysciences, Inc., Warrington, PA, USA). Morphological identification was performed under a BX61 microscope (Olympus, Japan). Measurements and counts were taken using the Stream Motion software (Olympus, Japan) and a video camera connected to the microscope. Identification was performed based on the keys available for Crete and adjacent regions [10,11,12,13,14].
The terminology adopted for the characters is the most recent one for phlebotomine sand flies . The following measurements were made for the specimens of the new species for both sexes: flagellomeres 1, 2 and 3, labrum-epipharynx. For males, we also measured the lengths of the parameral sheath, the distance from the tubercle to the tip of the parameral sheath (indicated as the distance from the tip of the aedeagus to the subterminal tooth by Artemiev ), sperm pump, aedeagal ducts, gonocoxite, beginning and ending of the tuft of internal setae of the gonocoxite, number of internal setae of the gonocoxal tuft and area of the internal tuft of setae of the gonocoxite (Fig. 1).
The percentage of tuft length vs gonocoxite length was calculated as follows: tuft length × 100)/gonocoxite length. The median tuft position vs gonocoxite was calculated as follows: beginning of the tuft of internal setae of the gonocoxite + tuft length/2) × 100/gonocoxite length. Drawings were made using a camera lucida.
DNA extraction and cytochrome b mtDNA sequencing
Following morphological identification, DNA was extracted using DNeasy Blood & Tissue Kit (Qiagen, Hilden, Germany) from 16 individuals (males and females) of Ph. (Adlerius) creticus n. sp. collected during previous surveys at four localities in Crete as well as from two specimens of Ph. simici Nitzuescu, 1931 from Crete and three specimens of Ph. balcanicus Theodor, 1958 from Iran. cytb PCR assays were performed using primers and conditions described by Esseghir et al.  and sequenced both directions according to Sanger’s method using the same primers as in the PCR assay.
Alignment and genetic distances
The newly generated sequences data were aligned with cytb sequences available on GenBank (Table 2), which belong to Adlerius spp. [17,18,19,20,21]. Multiple sequence alignments were performed using ClustalW  as implemented in MEGA7 . Genetic distances were calculated between and within species using Tamura-Nei model , also in MEGA7 .
The optimal nucleotide substitution models were identified using PartitionFinder (PF) v.2.1.1 . We ran PF two different times using the greedy search algorithm with linked branch lengths in calculations of likelihood scores under the Bayesian information criterion (BIC). The difference between these two runs was the restriction of candidate models to only those that are available in either MrBayes v.3.2.6  or PhyML v.3.0 . The models which included both gamma distribution (G) and invariable sites (I) were ignored .
Phylogenetic trees were constructed using Bayesian Inference (BI) and Maximum Likelihood (ML) methods. The BI analysis was performed in MrBayes v3.2.6  with 4 runs and 8 chains per run for 10,000,000 generations, with a sampling frequency of 100. From the sampled trees, 25% were discarded as ‛burn-inʼ phase and therefore, a majority rule consensus tree relied on the remaining trees and posterior probabilities were calculated as the percentage of samples recovering any particular clade . ML analysis was performed with PhyML v.3.0  with nearest-neighbor-interchange search, bio-neighbor joining starting tree under the suggested models selected in PF. Bootstrap values were estimated by 1000 replicates . Phlebotomus perfiliewi represented the outgroup in the phylogenetic analyses.
MALDI-TOF protein profiling
Samples that were subjected to MALDI-TOF MS analysis were processed as previously described . It was demonstrated that mass spectrometry-based approach is not suitable for specimens collected by sticky traps ; therefore, only specimens collected by CDC light traps were included into the assay. Within one month after the collection in the field, specimens stored in 70% ethanol were dissected, heads and terminalia were mounted by CMCP-10 mounting medium (Polysciences) for morphological typing, rest of abdomens were stored for DNA isolation and cytb sequencing as described above and thoraxes were ground by a manual BioVortexer homogenizer (BioSpec, Bartlesville, USA) with sterile disposable pestles in 10 μl of 25% formic acid. Two µl of the homogenate were mixed with 2 µl of freshly prepared MALDI matrix, an aqueous 60% acetonitrile/0.3% TFA solution of sinapinic acid (30 mg/ml; Sigma-Aldrich, St. Louis, USA). One μl of this mixture were spotted directly onto a steel MALDI plate in duplicates. Protein mass spectra were measured using an Autoflex Speed MALDI-TOF spectrometer (Bruker Daltonics, Billerica, USA) in a mass range of 4–25 kDa and compared by FlexAnalysis 3.4 software. For cluster analysis and species identification, the protein profiles were processed using MALDI Biotyper 3.1 and searched against an in-house database that comprises reference spectra of 23 different sand fly species including following Adlerius species: Ph. arabicus Theodor, 1953, Ph. balcanicus, Ph. halepensis Theodor, 1958 and Ph. simici.
Sand fly data
The first published study of the sand fly fauna of Crete appeared in 1910 and since then, 16 more works were published, revealing the presence of 10 Phlebotomus spp. and 2 Sergentomyia spp. (Table 3).
Sand fly sampling in 2014–2019
Overall, 608 sand fly specimens were collected which corresponded to nine different species with the most common and abundant species being Se. minuta and Ph. neglectus (Table 4). In addition to species known in Crete, 151 specimens of Ph. (Adlerius), morphologically close to Ph. balcanicus and Ph. zulfagarensis were identified. These specimens were further analysed and their morphological and molecular identification is described below.
Table 5 sums up the total number of specimens, of each species, caught in present samplings and previous publications. More than 30,000 specimens have been reported and about 63.27% of them were identified as the medically important species Ph. neglectus. Moreover, a significant proportion of the specimens were identified as Ph. similis (17.4%) and Ph. papatasi (9.49%). In all prefectures except Chania, Ph. neglectus (> 50%) and Ph. similis (~30%) were the most common and abundant species. In Chania, Ph. papatasi was the most common and abundant species which comprised almost the 60% of the sand fly specimens (Fig. 2).
Morphological typing of Ph. creticus n. sp
Fifty-four male individuals collected by CDC light traps were used for the morphological analysis, originating from three localities (8 from Toplou, 23 from Xerokampos and 23 from Fodele). The descriptive statistics for 22 characters are given in Table 6. Eleven females from Xerokampos were used for the description, descriptive statistics for 6 characters are given in Table 7. To exclude a presence of two distinct populations, the normality of the distributions of each of the morphological characters was tested by Shapiro-Wilk test (P > 0.05) and was found not to be significantly different from normal.
The antennal formula in males shows variability. Of the 44 males examined exhibiting antennae, 36 have an antennal formula 2/f1-f3, 1/f4-f13 which appears as the most common formula. Five specimens exhibit the following formula: 2/f1-f4, 1/f5-f13 with very commonly one very small ascoid on f5. One specimen exhibits 2/f1-f2, 1/f3-f13. One specimen exhibits 2/f1-f5, 1/f6-f13 with a small ascoid on both f4 and f5. Interestingly, one specimen exhibits 2/f1-f4, 1/f5-f13 on the left antenna and 2/f1-f3 + f5, 1/f4 + f6-f13 on the right one.
The cytb gene was successfully amplified and sequenced and the final dataset consisted of an alignment of 446 bp. The pairwise distances between species ranged from 0.044 to 0.188, while the closest related species to Ph. creticus n. sp. appeared to be Ph. balcanicus with a mean distance of 0.044. The mean distance within the Ph. creticus n. sp. samples was 0.008. The pairwise mean distances between and within species are provided in Table 8.
The best-fit nucleotide substitution model for cytb was Hasegawa-Kishino-Yano (HKY) + I for all codon positions and for both MrBayes and PhyML. Both analyses led to phylogenetic trees with similar topologies and hence, only the consensus MrBayes tree including posterior probabilities and bootstrap values is presented here (Fig. 3). All species used in these analyses formed monophyletic clades and specifically, the specimens of Ph. creticus n. sp. formed a well-supported monophyletic clade. Phlebotomus creticus n. sp. appears to be more closely related to Ph. balcanicus than the other species included in the analyses.
MALDI-TOF protein profiling
In total, 28 specimens of five Phlebotomus species from Xerokampos (n = 26) and Fodele (n = 2) were analysed: Ph. creticus n. sp. (n = 12), Ph. killicki (n = 1), Ph. neglectus (n = 5), Ph. similis (n = 7) and Ph. simici (n = 3). Reproducible protein spectra with a high number of intense signals within the mass range of 4–25 kDa were generated for all analysed specimens. These protein profiles were species-specific and showed species-unique peaks that allowed reliable and conclusive differentiation of the analyzed specimens. Protein spectra of the specimens identified based on morphology as belonging to four known species were similar to the corresponding reference spectra of the respective species. Protein spectra of all Ph. creticus n. sp. specimens were identical and differed substantially from the spectra of other species as shown by a hierarchical cluster analysis (Fig. 4a). Moreover, they also showed a number of specific peaks, not shared by protein spectra of four Adlerius species in the reference database (Fig. 4b).
Family Psychodidae Newman, 1834
Genus Phlebotomus Rondani & Berté, 1840
Phlebotomus creticus Antoniou, Depaquit & Dvorak n. sp.
Type-locality: Xerokampos (35° 3′ 29.37′′ N, 26° 14′ 27.14′′ E; altitude: 19 m above sea level), Greece
Other localities: Fodele (35° 22′ 53.93′′ N, 24° 57′ 28.27′′ E; altitude: 48 m above sea level); Agia Roumeli (35° 13′ 24.77′′ N, 23° 56′ 7.53′′ E; altitude: 14 m above sea level); Toplou Monastery (35° 13′ 16.69′′ N, 26° 12′ 51.99′′ E; altitude: 159 m above sea level), Greece.
Type-material: The holotype male (accession no. ED10723) and five paratypes (3 males and 2 females, accession nos. ED10724, ED10725, ED10726 for male paratypes, accession nos. ED10727, ED10728 for female paratypes) have been deposited at the Laboratory of Entomology of the Muséum National dʼHistoire Naturelle, Paris, France. Two paratypes (1 male and 1 female) have been deposited at the Museum of Natural History of London, UK, under the accession numbers NHMUK-ENT-2020-42. Two paratypes (1 male and 1 female) have been deposited at the Museum of Natural History of Heraklion, Crete, Greece, under the accession numbers NHMC.85.4.17830.1 (male paratype) and NHMC.85.4.17830.2 (females paratype).
Representative DNA sequences: GenBank accession numbers MT501628-MT501638.
ZooBank registration: To comply with the regulations set out in Article 8.5 of the amended 2012 version of the International Code of Zoological Nomenclature (ICZN) , details of the new species have been submitted to ZooBank. The Life Science Identifier (LSID) of the article is urn:lsid:zoobank.org:pub:D2A4798C-A59B-4851-819F-DA238621A474. The LSID for the new name Phlebotomus creticus is urn:lsid:zoobank.org:act:218CBCD4-875D-48E1-B2B8-749D53E817DA.
Etymology: The new species is named after the island where it has been discovered.
Male [Based on 54 specimens. Counts and measurements indicated in the description are those of the holotype. More measurements for males are available in Table 6; Figs. 5, 6]. Total length 3.9 mm. Head. Occiput with 2 lines of well individualized setae. Clypeus 171 µm long, 92 µm wide, with c.19 setae randomly distributed, targeting center of clypeus. Eyes 242 µm high with c.120 facets. Interantennal suture incomplete. Interocular suture not reaching the interantennal suture. Flagellomeres: f1 (349 µm) longer than f2 (148 µm) + f3 (149 µm). Internal and external ascoids implanted more or less at the same level on f1 to f3. Ascoids not reaching the next articulation. Ascoidal formula: 2/f1-f3 1/f4-f13 (= 2/III-V 1/VI-XV). One distal papilla on flagellomeres f1, f2, f3, three sensillae on f12 and f13, five on f14. Palpi p1: 49 µm long, p2: 173 µm, p3: 192 µm, p4: 160 µm, p5: 400 µm. Palpal formula: 1, 4, 2, 3, 5. About 15 Newstead’s sensilla present on p3 only; no sensilla on other palpal articles. One distal spiniform seta on p3, 12 setae on p4 and 26 setae on p5. Labrum-epipharynx 288 µm long carrying long teeth at its top. Hypopharynx with 20 long apical teeth. Labial suture closed, narrow, in furca. Cibarium without teeth nor sclerotized area (= pigment patch) or sclerotized arch. Pharynx with an armature consisting of long teeth directed laterally or towards center. Cervix with 3 lateral cervical sensillae and 2 median sensillae on each side. Thorax. Sclerites pale coloured. One post-alar seta present on the mesonotum. Paratergital seta absent. A group of four proepimeral setae. Upper anepisiternal seta, lower anapisternal seta, anepimeral seta, metaepisternal seta and metaepimeral seta absent. Setae present on the anterior region of the katepisternum. Metafurca mounted in lateral view on all specimens. Vertical arms long, probably separate; horizontal arms long. Wings. Length: 2370 μm; width: 715 μm; r5: 1537 µm; α (r2): 475 μm; β (r2+3): 309 μm; γ (r2+3+4): 373 μm; δ: 120 μm; π: 90 μm; ε (r3): 692 μm; θ (r4): 1091 μm; width/γ: 1.92. Legs. Anterior leg: coxa: 337 µm; femur: 855 µm; tibia: 1021 µm; tarsomere i: 620 µm; sum of tii, tiii, tiv, tv: 770 µm. Median leg: coxa: 353 µm; femur: 870 µm; tibia: 1218 µm; tarsomere i: 704 µm; sum of tii, tiii, tiv, tv: 810 µm. Posterior leg: coxa: 394 µm; femur: 1005 µm; tibia: 1558 µm; tarsomere i: 887 µm; sum of tii, tiii, tiv, tv: 944 µm. Spines on the metafemur absent. Metatarsomere iii with 4 verticils including broad and thin spines. Abdomen. Setae randomly implantated on tergites II to V. Two papillae with hair present on tergites III to VI. Gonocoxite 441 μm long. Sclerotized band in the ventral margin abent; process absent; a median cluster of 71 setae on the internal side present. Gonostyle 204 µm long, with 5 spines: 2 distal spines implanted at the same level; 3 median spines: 1 ventral implanted on narrow tubercle and 2 dorsal implanted on wide tubercle. Parameres 462 µm long, simple, rounded at apex, with setae occupying inner face of distal half. Parameral sheath straight, 203 µm long, with very shallow subterminal tubercle 20 µm to the top. Aedeagal ducts 927 μm long, isodiametric and pointed at their tops. Sperm pump 145 μm long. Ejaculatory apodeme 119 µm long. Aedeagal ducts/sperm pump ratio: 6.39. Epandrial lobes slightly longer than gonocoxites, length 468 μm.
Female [Based on 11 specimens. Counts and measurements indicated in the following description are those of the paratype labelled Crete IT8 with some exceptions which are indicated. More measurements for females available in Table 7; Fig. 7]. Total length of the paratype Crete IT14: 3.6 mm long. Head. Occiput with two narrow lines of well individualized setae. Clypeus 135 µm long, 103 µm wide, with 20 setae randomly distributed, targeting the center of the anterior part of the clypeus. Eyes 193 µm high with about 90 facets. Interantennal suture incomplete. Interocular suture not reaching the interantennal one. Flagellomeres: f1 (208 µm) longer than f2 (98 µm) +f3 (100 µm). Internal and external ascoids implanted more or less at the same level on f1 to f3. Ascoids not reaching the next articulation. Ascoidal formula: 2/f1-f13 (= 2/III-XV). One papilla on flagellomeres f1, f2, f3, three sensillae on f12 and f13, five on f14. Palpi p1: 40 µm long, p2: 153 µm, p3: 160 µm, p4: 130 µm, p5: 315 µm. Palpal formula: 1, 4, 2, 3, 5. Presence of about 15 Newstead’s sensilla on p3. Absence of Newstead’s sensilla on the other palpal articles. One distal spiniform setae on p3,7 on p4 and 31 on p5. Labrum-Epipharynx 268 μm long. f1/E = 0.78. Hypopharynx with about 15 distal long teeth on each side. Maxillary lacinia exhibiting 4 external and 20 internal teeth. Labial suture closed, narrow, in furca. Cibarium without teeth nor sclerotized area (= pigment patch) or sclerotised arch. Pharynx with a triangular armature consisting of elongated teeth directed towards the center. Cervix with two lateral cervical sensillae and two median ones on each side. Thorax. Pale coloured sclerites. Presence of one post-alar seta on the mesonotum. Absence of paratergital seta. A group of six proepimeral setae. Absence of upper anepisiternal seta. Absence of lower anapisternal seta. Absence of anepimeral seta. Absence of metaepisternal seta. Absence of metaepimeral seta. Presence of setae in the anterior region of the katepisternum. Metafurca mounted in lateral view on all specimens. Long vertical arms probably separate. Long horizontal arms.738. Wings. Length = 2173 μm, width = 578 μm, r5 = 1424 µm, α (r2) = 431 μm, β (r2+3) = 225 μm, γ (r2+3+4) = 398 μm, δ = 101 μm, π = 61 μm. ε (r3) = 612 μm, θ (r4) = 944 μm. Width / γ = 1.70. Anterior leg: coxa = 273 µm; femur = 811 µm; tibia = 971 µm; tarsomere i = 595 µm; sum of tii, tiii, tiv, tv = 728 µm. Median leg: coxa = 328 µm; femur = 817 µm; tibia = 1150 µm; tarsomere i = 686 µm; sum of tii, tiii, tiv, tv = 775 µm. Posterior leg: coxa = 383 µm; femur = 890 µm; tibia = 1522 µm; tarsomere i = 848 µm; sum of tii, tiii, tiv, tv = 876 µm. Absence of spines on the metafemur. Metatarsomere iii with a distal verticil and a median one, including broad and thin spines. Abdomen. Setae randomly implantated on tergites II to V. Presence of papillae on tergites III to VII. Genitalia. Presence of about 8 + 1 setae on tergite VIII. Lack of protuberance on tergite IX. Spermathecae incompletely segmented. Basal part of the ducts wide with thick walls. Those of the paratype Crete IT8 have been collapsed during the mounting. The measurements indicated are those of the paratype crete7 mounted in Marc-André to be observed, measured and drawn before final process and mounting. Length of the ducts: 600 µm (including 100 µm of the wide basal part and 500 µm of the narrow ducts); length of the body: 100 µm. Genital fork 192 µm long. Cerci rounded at their top, 178 µm long. No seta observed on the sternite X.
In males, number of ascoids usually 2/f1-f3, 1/f4-f13, a range of 54–85 setae on the gonocoxal internal tuft equally distributed between the proximal and distal halves of the gonocoxite.
Microhabitat preferences of Ph. creticus n. sp.
At Fodele, Ph. creticus n. sp. was captured in four CDC light traps placed near or in cave entrances, facing southeast and south. The vegetation around was rich. Mainly chicken, rats and lizards dwell around the traps. No specimens of Ph. creticus n. sp. were collected by sticky traps. At Xerokampos, the species was captured in two traps placed at the entrance of two shallow caves facing southeast, at 200 m distance from the shore. The vegetation around the caves was composed of phrygana, biotope typical for island of Crete. At Agia Roumeli the species was captured in two traps placed at the entrance of two limestone rock caves, one facing east and the other facing south, in a separate beach located west of Agia Roumeli, just a few meters from the sea. The vegetation above the caves was composed from sparse phrygana. In Toplou Monastery, it was trapped in caves along the wall just before the Monastery.
The island of Crete is an important and long-time active region of leishmaniases transmission and still provides ample numbers of human cases annually. As sand flies are the only proven vectors of these diseases in the Mediterranean basin, there is a need for sustained monitoring of the sand fly fauna and its role in Leishmania transmission. While studied for a long time, there are surprising new findings regarding the presence of sand fly species as documented by a recent description of Phlebotomus (Transphlebotomus) killicki Dvorak, Votypka & Volf, 2015 in the localities at the southern coast of the island . The results of the presented study reveal that there is still more to discover.
Published data as well as new results from the present sampling (2014–2019) show that the most widespread and abundant species in Crete is Ph. neglectus, accounting for more than 60% of the recorded sand fly specimens. It is a proven vector of L. infantum in Greece, several other Balkan countries and the western part of Turkey . Its abundance in all studied areas of Crete apparently contributes to the geographical distribution of VL cases throughout the island and as other species of the subgenus Larroussius are markedly scarce, we may conclude that Ph. neglectus remains a sole vector of the disease. The second most common and abundant species on the island is Phlebotomus similis. It is morphologically similar and phylogenetically closely related to Ph. sergenti [51, 52] and is regarded as a suspected vector of L. tropica in regions where Ph. sergenti is not present, including Crete . The fact that Ph. similis was found in all foci of human CL further fosters the considerations of this species as a yet unproven vector of L. tropica in the island. Phlebotomus papatasi is a proven vector of L. major in the Middle East and Northern Africa and it is known to be responsible for the phlebovirus infections in the Old World [53, 54] including Crete . However, L. major, the causative agent of zoonotic cutaneous leishmaniasis, does not circulate in Crete, probably due to absence of suitable reservoir species (gerbils). Vectorial competence of Sergentomyia spp. to mammalian diseases is still controversial and experimental studies that would conclusively test it are scarce  despite growing circumstantial evidence that suggest incrimination of some species in the transmission of human leishmaniasis and phleboviruses . Of two Sergentomyia species previously recorded from Crete, only Se. minuta was found in this study. Bloodmeal analyses of population in southern Portugal recently demonstrated that this widely distributed Mediterranean species may be at least partially anthropophilic and the detection of Leishmania DNA in specimens from the same area emphasizes the need of further studies about the role of Sergentomyia species in leishmaniases transmission cycles .
Adlerius specimens collected at different sites throughout the island during 2014–2019 did not meet the criteria of any known species and exhibited unique morphological characters. Moreover, obtained sequences of cytb, when compared with sequences of species within the subgenus Adlerius available in public databases, were substantially different. That led us to the conclusion that the collected specimens represent a new species. The identity of this new species is based on a detailed morphological analysis of decisive characters that is further supported by sequencing of cytb gene, a widely used genetic marker, and for the first time in a description of a new sand fly species, also by comparison of species-specific protein spectra acquired by MALDI-TOF mass spectrometry.
According to Artemiev , the subgenus Adlerius includes 17 described and two undescribed species. While the females mostly appear to be undistinguishable, the main characters to identify the males of this subgenus are the antennal formula, used in all available identification keys [10, 11, 59] and the number of setae of the internal tuft of the gonocoxite as well as the position of this tuft on the gonocoxite.
Considering that the number of ascoids in Ph. creticus n. sp. was usually 2/f1-f3, 1/f4-f13, the most similar species are Ph. angustus Artemiev, 1978, Ph. balcanicus, Ph. comatus Artemiev, 1978, Ph. kyreniae Theodor, 1958, Ph. salangensis Artemiev, 1978 and Ph. zulfagarensis Artemiev, 1978 . However, this parameter varied not only in Ph. creticus n. sp., its infraspecific variation is known also in other sand fly species. Therefore, the antennal formula may not be alone a useful parameter for species identification of Adlerius sand flies.
Phlebotomus creticus n. sp. exhibits a range of 54–85 setae on the gonocoxal internal tuft (mean: 69). This count excludes the identification of Ph. balcanicus (92–130 setae after Artemiev  and more than 100 in the original description [reference number]), Ph. comatus (126–220 setae) and Ph. kyreniae (30–40 setae). There is an overlap related to this number of setae for Ph. creticus n. sp. and Ph. angustus (35–69 setae), Ph. salangensis (40–85 setae) and Ph. zulfagarensis (66–72 setae). However, these three species exhibit group of setae which is completely or at 90% in the proximal half of the gonocoxite  whereas that of Ph. creticus n. sp. is equally distributed between the proximal and distal halves of the gonocoxite (Table 6, Fig. 1). Moreover, all three species are known to be distributed in Asia, in regions very distant from Crete. In contrast, in Ph. balcanicus the hair group is located mainly on the basal half of the coxite . Consequently, regarding the typological systematics, there is substantial evidence to consider Ph. creticus as a distinct species.
Phlebotomus creticus n. sp. was recorded at different localities at both northern and southern coast of Crete divided by a mountain range of a considerable height, which further supports a long presence rather than a recent introduction. However, it is probably not a common species. Phlebotomus creticus n. sp. was not caught on sticky traps placed in the openings of holes on natural or manmade walls and ground. We assume that it preferentially rests inside shallow caves in limestone rocks. In such specific biotopes it could be a dominant species as demonstrated in Xerokampos. Its feeding preferences should be investigated, although wildlife hosts may be expected (mice, birds and possibly lizards). Beside this newly described species, Ph. simici of the same subgenus that had been reported from Crete in the past, was still recorded in this study, occurring even sympatrically at the type-locality Xerokampos. We may speculate that in previous entomological surveys, some of the specimens identified as Ph. simici or unidentified Adlerius sp. may be attributed to the newly described Ph. creticus n. sp.
The results of the molecular analyses support the description of a new species. Phylogenetic analyses of the cytochrome b gene, a mitochondrial marker widely used in phylogenetic studies of many insect groups including sand flies, strongly grouped all analyzed specimens of Ph. creticus n. sp. in a distinct monophyletic clade. This genetic marker was chosen as it provides the best coverage of the species within the subgenus Adlerius. Unfortunately, no sequences of any genetic marker are so far available for some species of the subgenus, including the three species with overlapping numbers of setae on the gonocoxal internal tuft. Genetic distance values obtained for Ph. creticus n. sp. and other compared Adlerius species are comparable with the distances recorded previously for other sand fly species, as shown by studies of morphologically- as well as genetically distinct species of the subgenera Larroussius and Phlebotomus  or Madaphlebotomus . Moreover, MALDI-TOF protein profiling demonstrated that all processed specimens of Ph. creticus n. sp. produced unique, reproducible and species-specific profiles that clearly differentiate them from other species outside and within the subgenus Adlerius. This method of mass spectrometry has recently become a popular tool for species identification of various organisms including arthropod vectors as it is simple, rapid and cost-effective . Here, we demonstrated for the first time that it can also successfully complement traditional morphological approach and established DNA-based molecular taxonomy in the process of revealing yet unrecognized sand fly species.
Our recent findings urge the need for a revision of the subgenus Adlerius using both morphological and molecular approaches. Species boundaries are not well defined and the vicariance of this group probably occurred recently as for other groups of phlebotomine sand flies such as Phlebotomus, Larroussius or Paraphlebotomus [16, 18, 62].
In this study we present a review of sand fly species recorded in the past and at present in Crete, an island with ongoing transmission of two Leishmania species due to the presence of competent sand fly vectors. The importance of this research is highlighted by the geographical position of the island and the current possibility of accidental introduction of more Leishmania species due to human migration and other activities. According to our findings, 10 Phlebotomus spp. and 2 Sergentomyia spp. were recorded, with the most common and abundant species being Ph. neglectus. We may assume that the findings of Ph. mascittii reported prior to the description of Ph. killicki may be attributed to the latter species. We identify and describe a new species Phlebotomus (Adlerius) creticus n. sp. from various localities in Crete. Its identification is based on morphological characters of the male genitalia that particularly differentiate it from related species of the subgenus Adlerius. The identity of the newly described species was confirmed by two molecular approaches (MALDI-TOF protein profiling and cytb sequence analysis). As there is no data on the vectorial competence and capacity of this new species, its potential role in the autochtonous transmission cycles of Leishmania shall be further studied.
Availability of data and materials
The datasets supporting the conclusions of this article are included within the article and its additional files. The holotype and paratypes of a newly described species Phlebotomus creticus n. sp.were deposited in three repositories as described above. The newly generated sequences were deposited in the GenBank database under the accession numbers MT501623-MT501638.
Bayesian information criterion
- cytb :
Cytochrome b gene
- MALDI-TOF MS:
Matrix-assisted laser desorption/ionization time-o-flight mass spectrometry
Polymerase chain reaction
Maroli M, Feliciangeli MD, Bichaud L, Charrel RN, Gradoni L. Phlebotomine sandflies and the spreading of leishmaniases and other diseases of public health concern. Med Vet Entomol. 2013;27:123–47.
Antoniou M, Gramiccia M, Molina R, Dvorak V, Volf P. The role of indigenous phlebotomine sandflies and mammals in the spreading of leishmaniasis agents in the Mediterranean region. Euro Surveill. 2013;18:20540.
Svobodova M, Votypka J, Peckova J, Dvorak V, Nasereddin A, Baneth G, et al. Distinct transmission cycles of Leishmania tropica in 2 adjacent foci, Northern Israel. Emerg Infect Dis. 2006;12:1860–8.
Jacobson RL. Leishmania tropica (Kinetoplastida: Trypanosomatidae) - a perplexing parasite. Folia Parasitol. 2003;50:241–50.
Antoniou M, Haralambous C, Mazeris A, Pratlong F, Dedet J-P, Soteriadou K. Leishmania donovani causing visceral or cutaneous leishmaniasis in Europe: Cyprus evidence is an alarm call. Lancet Infect Dis. 2008;8:6–7.
Ntais P, Christodoulou V, Tsirigotakis N, Dokianakis E, Dedet J-P, Pratlong F, et al. Will the introduction of Leishmania tropica MON-58, in the island of Crete, lead to the settlement and spread of this rare zymodeme? Acta Trop. 2014;132:125–30.
Dvorak V, Shaw J, Volf P. Parasite biology: the vectors. In: Bruschi F, Gradoni L, editors. The leishmaniases: old neglected tropical diseases. Cham: Springer International Publishing; 2018. p. 31–77.
Moher D, Shamseer L, Clarke M, Ghersi D, Liberati A, Petticrew M, et al. Preferred reporting items for systematic review and meta-analysis protocols (PRISMA-P) 2015 statement. Syst Rev. 2015;4:1.
Blavier A, Laroche L, Randrianambinintsoa FJ, Lucas V, Gantier J-C, Léger N, et al. Phlebotomine sandflies (Diptera, Psychodidae) from the Ankarana tsingy of northern Madagascar: inventory and description of new taxa. Parasite. 2019;26:38.
Artemiev MM. A revision of sandflies of the subgenus Adlerius (Diptera, Phlebotominae, Phlebotomus). Zool Zhurnal. 1980;59:1177–92.
Lewis DJ. A taxonomic review of the genus Phlebotomus (Diptera: Psychodidae). Bull Br Mus Nat Hist (Ent). 1982;45:121–209.
Léger N, Pesson B, Madulo-Leblond G. Les phlébotomes de Grèce. Biol Gallo-hell. 1986;11:165–91.
Killick-Kendrick R, Tang Y, Killick-Kendrick M, Sang DK, Sirdar MK, Ke L, et al. The identification of female sandflies of the subgenus Larroussius by the morphology of the spermathecal ducts. Parassitologia. 1991;33:335–47.
Depaquit J, Léger N, Ferté H, Rioux JA, Gantier JC, Michaelides A, et al. Les Phlébotomes de l’île de Chypre III. Inventaire faunistique. arasite. 2001;8:11–20.
Galati EAB, Galvis-Ovallos F, Lawyer P, Léger N, Depaquit J. An illustrated guide for characters and terminology used in descriptions of Phlebotominae (Diptera, Psychodidae). Parasite. 2017;24:26.
Esseghir S, Ready PD, Killick-Kendrick R, Ben-Ismail R. Mitochondrial haplotypes and phylogeography of Phlebotomus vectors of Leishmania major. Insect Mol Biol. 1997;6:211–25.
Parvizi P, Naddaf SR, AlaeeNovin E. Molecular typing and phylogenetic analysis of some species belonging to Phlebotomus (Larroussius) and Phlebotomus (Adlerius) subgenera (Diptera: Psychodidae) from two locations in Iran. Iran J Arthropod Borne Dis. 2010;4:1–10.
Mahamdallie SS, Pesson B, Ready PD. Multiple genetic divergences and population expansions of a Mediterranean sandfly, Phlebotomus ariasi, in Europe during the Pleistocene glacial cycles. Heredity. 2011;106:714–26.
Krüger A, Strüven L, Post RJ, Faulde M. The sandflies (Diptera: Psychodidae, Phlebotominae) in military camps in northern Afghanistan (2007–2009), as identified by morphology and DNA ‘barcoding’. Ann Trop Med Parasitol. 2011;105:163–76.
Depaquit J, Bounamous A, Akhoundi M, Augot D, Sauvage F, Dvorak V, et al. A taxonomic study of Phlebotomus (Larroussius) perfiliewi s.l. Infect Genet Evol. 2013;20:500–8.
Latrofa MS, Dantas-Torres F, Weigl S, Tarallo VD, Parisi A, Traversa D, et al. Multilocus molecular and phylogenetic analysis of phlebotomine sand flies (Diptera: Psychodidae) from southern Italy. Acta Trop. 2011;119:91–8.
Thompson JD, Higgins DG, Gibson TJ. CLUSTAL W: improving the sensitivity of progressive multiple sequence alignment through sequence weighting, position-specific gap penalties and weight matrix choice. Nucleic Acids Res. 1994;22:4673–80.
Kumar S, Stecher G, Tamura K. MEGA7: molecular evolutionary genetics analysis version 7.0 for bigger datasets. Mol Biol Evol. 2016;33:1870–4.
Tamura K, Nei M. Estimation of the number of nucleotide substitutions in the control region of mitochondrial DNA in humans and chimpanzees. Mol Biol Evol. 1993;10:512–26.
Lanfear R, Frandsen PB, Wright AM, Senfeld T, Calcott B. PartitionFinder 2: new methods for selecting partitioned models of evolution for molecular and morphological phylogenetic analyses. Mol Biol Evol. 2016;34:772–3.
Ronquist F, Teslenko M, van der Mark P, Ayres DL, Darling A, Höhna S, et al. MrBayes 3.2: efficient Bayesian phylogenetic inference and model choice across a large model space. Syst Biol. 2012;61:539–42.
Guindon S, Dufayard J-F, Lefort V, Anisimova M, Hordijk W, Gascuel O. New algorithms and methods to estimate maximum-likelihood phylogenies: assessing the performance of PhyML 3.0. Syst Biol. 2010;59:307–21.
Yang Z. Computational molecular evolution. Oxford: University Press; 2006.
Huelsenbeck JP, Ronquist F. MRBAYES: Bayesian inference of phylogenetic trees. Bioinformatics. 2001;17:754–5.
Felsenstein J. Confidence limits on phylogenies: an approach using the bootstrap. Evolution. 1985;39:783–91.
Dvorak V, Halada P, Hlavackova K, Dokianakis E, Antoniou M, Volf P. Identification of phlebotomine sand flies (Diptera: Psychodidae) by matrix-assisted laser desorption/ionization time of flight mass spectrometry. Parasit Vectors. 2014;7:21.
Halada P, Hlavackova K, Dvorak V, Volf P. Identification of immature stages of phlebotomine sand flies using MALDI-TOF MS and mapping of mass spectra during sand fly life cycle. Insect Biochem Mol Biol. 2018;93:47–56.
ICZN. International Commission on Zoological Nomenclature: Amendment of articles 8, 9, 10, 21 and 78 of the International Code of Zoological Nomenclature to expand and refine methods of publication. Bull Zool Nomencl. 2012;69:161–9.
Parrot L. Notes sur les phlébotomes. XIX. Phlébotomes de Crète. Arch Inst Pasteur Alger. 1936;14:50–2.
Adler S, Theodor O, Witenberg G. Investigations on Mediterranean kala azar. XI. A study of leishmaniasis in Canea (Crete). Proc R Soc Lond. 1938;125:491–516.
Ristorcelli A. Sur les phlébotomes de l’île de Crète. Ann Parasitol Hum Comp. 1939;17:355–8.
Hertig M. Phlebotomus and residual DDT in Greece and Italy. Am J Trop Med Hyg. 1949;s1-29:773–809.
Hadjinicolaou J. Present status of Phlebotomus in certain areas of Greece. Bull World Health Organ. 1958;19:967–79.
Pesson B, Wallon M, Floer MT, Kristensen AR. Isoenzymatic studies of Mediterranean populations of sandflies of the subgenus Larroussius. Parassitologia. 1991;33:471–6.
Léger N, Pesson B, Madulo-Leblond G, Ferté H, Tselentis Y, Antoniou M, Les phlébotomes de Crète M. Résultats d’une enquête entomologique effectuée en Juillet 1988 et Août 1989. Biol Gallo-hell. 1993;20:135–43.
Chaniotis B, Spyridaki I, Scoulika E, Antoniou M. Colonization of Phlebotomus neglectus (Diptera: Psychodidae), the major vector of visceral leishmaniasis in Greece. J Med Entomol. 2000;37:346–8.
Ivović V, Patakakis M, Tselentis Y, Chaniotis B. Faunistic study of sandflies in Greece. Med Vet Entomol. 2007;21:121–4.
Ivović V, Chaniotis V, Vujanić M, Bobić B, Nikolić A, Klun I, et al. Life tables and reproductive parameters of Phlebotomus neglectus Tonnoir, 1921 (Diptera, Psychodidae) under laboratory conditions. Arch Biol Sci. 2010;62:153–7.
Christodoulou V, Antoniou M, Ntais P, Messaritakis I, Ivović V, Dedet J-P, et al. Re-emergence of visceral and cutaneous leishmaniasis in the Greek Island of Crete. Vector Borne Zoonotic Dis. 2012;12:214–22.
Alten B, Maia C, Afonso MO, Campino L, Jiménez M, González E, et al. Seasonal dynamics of phlebotomine sand fly species proven vectors of Mediterranean leishmaniasis caused by Leishmania infantum. PLoS Negl Trop Dis. 2016;10:e0004458.
Langeron M. Phlébotomes capturés en Crète. Ann Parasitol Hum Comp. 1923;1:108.
Aransay AM, Scoulica E, Chaniotis B, Tselentis Y. Typing of sandflies from Greece and Cyprus by DNA polymorphism of 18S rRNA gene. Insect Mol Biol. 1999;8:179–84.
Blanc G, Caminopetros J. Enquiry into Oriental sore in Crete. Reflections suggested on the etiology and mode of dispersion of the disease. Ann Inst Pasteur. 1921;35:151–66.
Kasap OE, Dvorak V, Depaquit J, Alten B, Votypka J, Volf P. Phylogeography of the subgenus Transphlebotomus Artemiev with description of two new species, Phlebotomus anatolicus n. sp. and Phlebotomus killicki n. sp. Infect Genet Evol. 2015;34:467–79.
Birt C. Phlebotomus fever in Malta and Crete. J R Army Med Corps. 1910;14:236–58.
Depaquit J, Ferté H, Léger N, Lefranc F, Alves-Pires C, Hanafi H, et al. ITS 2 sequences heterogeneity in Phlebotomus sergenti and Phlebotomus similis (Diptera, Psychodidae): possible consequences in their ability to transmit Leishmania tropica. Int J Parasitol. 2002;32:1123–31.
Depaquit J, Léger N, Ferté H. Le statut taxinomique de Phlebotomus sergenti Parrot, 1917, vecteur de Leishmania tropica (Wright, 1903) et Phlebotomus similis Perfiliev, 1963 (Diptera-Psychodidae). Approches morphologique et morphométrique. Corollaires biogéographiques et épidémiologiques. Bull Soc Pathol Exot. 1998;91:346–52.
Killick-Kendrick R. Phlebotomine vectors of the leishmaniases: a review. Med Vet Entomol. 1990;4:1–24.
Depaquit J, Grandadam M, Fouque F, Andry PE, Peyrefitte C. Arthropod-borne viruses transmitted by phlebotomine sandflies in Europe: a review. Euro Surveill. 2010;15:19507.
Alwassouf S, Christodoulou V, Bichaud L, Ntais P, Mazeris A, Antoniou M, et al. Seroprevalence of sandfly-borne phleboviruses belonging to three serocomplexes (Sandfly fever Naples, Sandfly fever Sicilian and Salehabad) in dogs from Greece and Cyprus using neutralization test. PLoS Negl Trop Dis. 2016;10:e0005063.
Sadlova J, Dvorak V, Seblova V, Warburg A, Votypka J, Volf P. Sergentomyia schwetzi is not a competent vector for Leishmania donovani and other Leishmania species pathogenic to humans. Parasit Vectors. 2013;6:186.
Maia C, Depaquit J. Can Sergentomyia (Diptera, Psychodidae) play a role in the transmission of mammal-infecting Leishmania? Parasite. 2016;23:55.
Maia C, Parreira R, Cristóvão JM, Freitas FB, Afonso MO, Campino L. Molecular detection of Leishmania DNA and identification of blood meals in wild caught phlebotomine sand flies (Diptera: Psychodidae) from southern Portugal. Parasit Vectors. 2015;8:173.
Artemiev MM, Neronov VM. Distribution and ecology of sandflies of the old world: genus Phlebotomus. Moscow: USSR Academy of Sciences; 1984.
Esseghir S, Ready PD, Ben-Ismail R. Speciation of Phlebotomus sandflies of the subgenus Larroussius coincided with the late Miocene-Pliocene aridification of the Mediterranean subregion. Biol J Linn Soc. 2000;70:189–219.
Yssouf A, Almeras L, Raoult D, Parola P. Emerging tools for identification of arthropod vectors. Future Microbiol. 2016;11:549–66.
Depaquit J, Ferte H, Leger N, Killick-Kendrick R, Rioux JA, Killick-Kendrick M, et al. Molecular systematics of the phlebotomine sandflies of the subgenus Paraphlebotomus (Diptera, Psychodidae, Phlebotomus) based on ITS2 rDNA sequences. Hypotheses of dispersion and speciation. Insect Mol Biol. 2000;9:293–300.
We would like to thank Elif Konuksal for his help during the fieldwork in Fodele and Thord Fransson and Dimitris Chatzidakis for their help during the fieldwork in Xerokampos. We would also like to thank Rea L. Antoniou-Kourounioti for the statistical analyzes of the morphological characters and Vasiliki Christodoulou for preparing a map of Crete used in the graphical abstract to depict collection sites of Ph. creticus n. sp.
This research was supported by the projects BIOCEV CZ.1.05/1.1.00/02.0109 and CePaViP CZ.02.1.01/0.0/0.0/16_019/0000759 from the European Regional Development Fund. Additional support from the Czech Infrastructure for Integrative Structural Biology (LM2015043 funded by MEYS CR) is gratefully acknowledged. Thanks are due to Nicole Léger, Bernard Pesson, Geneviève Madulo-Leblond, Byron Papadopoulos and Hubert Ferté for collecting the specimens from Toplou Monastery and to Bulent Alten for their support and advice throughout this work.
Ethics approval and consent to participate
The study involved collecting adult sand flies in the open field, from domestic animal shelters and human houses by traps. Local or regional ethics committee approval was not required for such work. Setting the traps near houses or in fields was performed with informed consent and cooperation of the owners and the local authorities.
Consent for publication
The authors declare that they have no competing interests.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
About this article
Cite this article
Dvořák, V., Tsirigotakis, N., Pavlou, C. et al. Sand fly fauna of Crete and the description of Phlebotomus (Adlerius) creticus n. sp. (Diptera: Psychodidae). Parasites Vectors 13, 547 (2020). https://doi.org/10.1186/s13071-020-04358-x