Research | Open | Published:
Species composition of sand flies and bionomics of Phlebotomus papatasi and P. sergenti (Diptera: Psychodidae) in cutaneous leishmaniasis endemic foci, Morocco
Parasites & Vectorsvolume 9, Article number: 60 (2016)
Cutaneous Leishmaniasis (CL) is one of the most neglected tropical diseases in Morocco. Leishmania major and L. tropica are the main culprits identified in all endemic foci across the country. These two etiological agents are transmitted by Phlebotomus papatasi and P. sergenti, the two most prevalent sand fly species in Morocco. Previous studies reflected gaps of knowledge regarding the environmental fingerprints that affect the distribution of these two potential vectors across Morocco.
The sand flies were collected from 48 districts across Morocco using sticky paper traps. Collected specimens were preserved in 70 % ethanol for further processing and identification. Male and female densities were calculated in each site to examine their relations to the environmental conditions across these sites. The study used 19 environmental variables including precipitation, aridity, elevation, soil variables and a composite representing maximum, minimum and mean of day- and night-time Land Surface Temperature (LST), and Normalized Difference Vegetation Index (NDVI).
A total of 11,717 specimens were collected during this entomological survey. These specimens represented 11 species of two genera; Phlebotomus and Sergentomyia. Correlations of the sand fly densities with the environmental variables were estimated to identify the variables which influence the distribution of the two potential vectors, Phlebotomus papatasi and P. sergenti, associated with all CL endemic foci across the country. The density of P. papatasi was most affected by temperature changes. The study showed a significant positive correlation between the densities of both sexes of P. papatasi and night-time temperatures. Both P. papatasi and P. sergenti showed a negative correlation with aridity, but, such correlation was only significant in case of P. papatasi.
NDVI showed a positive correlation only with densities of P. sergenti, while, soil PH and soil water stress were negatively correlated with the densities of both males and females of only P. papatasi.
Our results identified the sand fly species across all CL endemic sites and underlined the influences of night-time temperature, soil water stress and NDVI as the most important variables affecting the sand fly distribution in all sampled sites. This preliminary study considered the importance of these covariates to anticipate the potential distribution of P. papatasi and P. sergenti in Morocco.
Cutaneous and visceral leishmaniasis (kala-azar) represent the most neglected tropical disease across the world . In Morocco, cutaneous leishmaniasis (CL) is endemic and constitutes a major public health threat . CL is widely distributed in three nosogeographic entities across Morocco; zoonotic cutaneous leishmaniasis (ZCL; caused by L. major) located in the arid regions along the northern edge of the Sahara desert, anthroponotic cutaneous leishmaniasis (ACL; caused by L. tropica) in the semi-arid regions of central and south-western Morocco, and CL caused by L. infantum in the northern regions of the country [3, 4].
Moroccan Ministry of Health has reported 24,804 cases of ZCL and 16,852 cases of ACL in 2004–2013 alone . CL caused by L. infantum remains a very rare condition with few sporadic cases occurring, especially in northern Morocco . Even with the high number of cases reported during recent years, leishmaniasis remains underestimated, especially in remote and rural communities where there is no access to health facilities. CL is the least known and most neglected disease, especially among males in all endemic areas in Morocco: males mostly depend on traditional medicine which becomes overpowered in the endemic areas. Lack of awareness among males reflected the higher incidence of leishmaniasis occurred among females and children .
ZCL transmission was maintained for a long term in rodent reservoirs, for example, Meriones shawi grandis is recognized as the main mammalian reservoir host to maintain local circulation of L. major  . On the other hand, human was identified as the only reservoir host for ACL caused by L. tropica ; however, zoonotic transmission for this species was also reported from other countries where the parasite was isolated from rock hyrax, and gerbil in Israel and Egypt, respectively [8, 9].
Leishmania infection is transmitted locally to a human host by the bite of an infected female sand fly (Diptera: Psychodidae) of the genus Phlebotomus. Phlebotomus (Phlebotomus) papatasi Scopoli, 1786 and P. (Paraphlebotomus) sergenti Parrot, 1917 were identified as the only proven vectors for ZCL and ACL, respectively [10, 11].
Over the past decade, the epidemiological situation of CL has changed significantly in Morocco; the disease was reported and geographically expanded to new habitats which were never reported as leishmaniasis endemic foci. This rapid expansion of CL in the country increased the public health problems associated with the disease epidemics. Previous studies revealed that this expansion in the disease dynamics may be associated with range expansions of the vector populations in response to climate change [12, 13].
P. papatasi and P. sergenti were found highly dependent on environmental conditions [14–16]. Environmental factors may affect several ecological and biological processes that act directly on controlling the vector’s geographical distribution, abundance and reproductive rates. All of these associations between environmental conditions and the distribution of sand fly vectors reveal a gap of knowledge to the current situation in Morocco; however, evidences are available from other countries in the region [17, 18].
This current contribution is a preliminary analysis to update the sand fly species composition across Morocco and examine the relationship between the potential distribution of P. papatasi and P. sergenti in the region in response to climatic (i.e. temperature, aridity and precipitation), topographic, soil and the greenness variables; taking advantage of the recent high resolution satellite data available across the study area. The preliminary results presented here are the first set of analysis to evaluate the importance of multidimensional factors in limiting the distribution of two potential CL vectors. They represent a baseline for more detailed studies across the country to better understand the ecology of sand fly vectors and guiding the national control program.
Morocco is a northern African country, bordering the North Atlantic Ocean and the Mediterranean Sea. The current study covered sand fly sampling in 48 districts across central and southern Morocco (Fig. 1) with altitude ranges between sea level and up to 2400 m above sea. These districts represented five administrative regions: Guelmim-Es Smara, Marrakech-Tensift-Al Haouz, Meknes-Tafilalete, Souss-Massa-Draa and Tadla-Azilal. The climate in Morocco is mostly Mediterranean; however, seven bioclimatic regimes occur due to topographic differences across the country (e.g., high mountains, sub-humid, humid, semi-arid with cold winter, semi-arid with warm and temperate winter, arid and Saharan climate). The total population in Morocco is about 33,848,242 with an urbanization rate of 60.3 % .
Sand fly collection and processing
Sand flies were collected from 48 districts across Morocco to include most foci identified as endemic in regard to CL transmission. Sand flies were collected by the sticky paper traps during three successive months of 2013 (i.e. April–June). The selection of this sampling period was in light of the previous observations reported from Morocco; (1) April-June corresponds to the highest peak of the sand fly abundance in Morocco [20, 21], and (2) high incidence of CL cases in this time interval annually [16, 22]. Sampling was carried out in rural areas around the human population (e.g., intradomestic, peridomestic, stables, barbacanes, old walls and rock crevices sites). The sticky traps were distributed so that each site had 30 traps/night (3 nights/site), with a total number of 4320 traps in 48 sites. Sticky traps were set at 6 pm and collected around 6 am of the next day. The trapped sand flies were collected in plastic bags coded with location, GIS coordinates, and date and transferred to the Laboratory of Ecology and Environment in Marrakesh, Morocco for further processing and identification.
In the laboratory, specimens were placed in 70 % ethanol until mounting. These specimens were then cleared in Mark André medium  and mounted on glass slides in Hoyer medium  prior to identification using the local morphological keys [25, 26]. The identification was based on the morphology of diagnostic characters including pharyngeal armature, female spermathecae and male genitalia. P. sergenti species identity was confirmed based on the morphological characters described in a previous study . Species and forms of Phlebotomus perniciosus complex  were identified via diagnostic morphological characters; dilatation of the distal parts of the spermathecal ducts in females  and morphology of the copulatory valves and number of coxite hairs in males [30–32]. Sand fly species evenness and Bray-Curtis similarity between all sites were calculated using BioDiversity Professional statistics analysis software . Finally, sand fly density was calculated based on the number of sand flies/m2 of sticky traps/night  to examine its relationship to the environmental changes across sampled sites.
The study used several climatic, greenness, topographic and soil variables to examine the relationship between these environmental factors and the abundance of two potential vectors associated with circulation of CL across the country. These variables were obtained from different sources; (1) Land Surface Temperature (LST) and Normalized Difference Vegetation Index (NDVI) were downloaded from the Moderate Resolution Imaging Spectroradiometer (MODIS) satellite imagery for April to June 2013 corresponding to the same time of our collection. These data are available from Land Processes Distributed Active Archive Center obtained through the NASA Reverb Echo data portal (http://reverb.echo.nasa.gov/reverb). Maximum, minimum and mean of day-time temperature, night-time temperature and NDVI were calculated as described elsewhere [17, 34] for each grid cell across the entire study period so that a final set of 9 variables were used. (2) Precipitation data were obtained from the National Oceanic and Atmospheric Administration (NOAA) National Weather Center (http://www.cpc.noaa.gov). The precipitation data is based on daily estimates across the study area. We calculated the maximum and mean estimates of precipitation for April-June across each sampling site. (3) Aridity index, and soil water stress were downloaded from Consortium for Spatial Information (CGIAR-CSI; http://www.cgiar-csi.org/). (4) Elevation data were obtained from the Shuttle Radar Topography Mission (SRTM; http://srtm.usgs.gov/). (5) Soil data were obtained from the International Soil Reference and Information Centre (ISRIC) via the link (http://www.isric.org) to summarize chemical and physical soil characteristics across Morocco. These data represent the soil values in diverse depths across points sampled for soil properties, so, we calculated the mean of depths for the same variable. We used 5 variables reflecting the PH, bulk density, silt, clay and sand contents.
In sum, we have a set of 19 variables representing different climatic, topographic, greenness and soil variables. The values of each variable were extracted for each of 48 sampling sites and based on their digital coordinates using ArcGIS 10.3 (Environmental Systems Resource Institute, Redlands, California). The selection of these environmental variables was based on their direct and indirect influences on ecology and distribution of the sand fly populations [15, 35, 36].
Data were analyzed using Minitab Version 17.1.0 . An ANOVA test was used to compare the collections in different sites across the study area. Chi-squared test was used to test the deviation from the expected sex ratio equilibrium (1: 1). The correlations between sand fly densities and each of climate, elevation, greenness and soil variables were assessed using Pearson’s (r) correlation coefficient. Differences were considered significant only if P < 0.05.
Sand fly species composition and sex ratios
A total of 11,717 sand flies were collected from 48 districts representing five governorates across Morocco, during April-June 2013 (Additional file 1; Table 1). The collected flies represented 11 species belonging to two genera (Phlebotomus and Sergentomyia); Sergentomyia minuta (N = 2653), S. fallax (N = 2053), P. perniciosus (N = 1688), P. papatasi (N = 1631), P. sergenti (N = 1541), P. longicuspis (N = 1390), P. ariasi (N = 334), S. dreyfussi (N = 189), P. alexandri (N = 177), S. africana (N = 39), P. chabaudi (N = 22). S. minuta was the most predominant of the all collected flies; however, P. perniciosus was the dominant species of the genus Phlebotomus. There was a significant difference in sand fly collections across sampled sites (P = 0.00). No morphological anomalies were observed in the specimens of all collected species, except that P. perniciosus was collected only as atypical form (PNA).
The sex ratios (Females: Males) showed that overall collected males were approximately three times more than females (F/M = 1: 2.92; Table 1). Chi-squared test revealed a significant difference between the overall sex ratio and the expected sex ratio equilibrium (1: 1) for all species (P = 0.00).
Sand fly diversity across sampling sites
Species evenness revealed differences in sand fly diversity among the different sampling sites with the highest overall diversity in Tagounite (J’ = 0.986) and the least diversity in Abadou where six species only occurred (J’ = 0.328). The species evenness and dendrogram of Bray-Curtis similarity are shown for comparison between sites in the supplementary materials (Additional file 2 and Additional file 3).
Effect of environment on the distribution of CL vectors
P. papatasi and P. sergenti co-occurred in 39 out of 48 sites sampled across Morocco. There was a difference between the densities of the P. papatasi and P. sergenti among the different sampling sites in regard to temperature and NDVI (Figs. 2 and 3). The correlations between abundance of P. papatasi and P. sergenti with several environmental factors were examined to understand the possible effect of these variables on the dynamics of these two vectors across the country (Table 2). The correlation of males and females differed in regard to each environmental condition. The densities of males and females P. papatasi were the most affected by temperature changes across all sampling sites. The study showed a significant positive correlation between the densities of both sexes of P. papatasi and night-time temperatures; however, no significant correlations were found with the day-time temperatures (Table 2). On the other hand, both sexes of P. sergenti showed no significant correlation with both day-time and night-time temperatures.
Male and female densities of both P. papatasi and P. sergenti showed negative correlation with the increase in precipitation; however, this correlation was not significant for both species. The densities of both species sexes showed negative correlation with aridity, but, this negative correlation was only significant in case of P. papatasi (Fig. 4; Table 2).
NDVI showed positive correlation with densities of P. sergenti males and females (Table 2), but a decrease in density of both sexes of P. papatasi was observed with increase of NDVI values. The correlation of NDVI and sand fly density was only significant in case of P. sergenti females.
P. papatasi was collected between 14 m and 1800 m altitude, while, P. sergenti was collected between 28 m and 1300 m altitude. The higher P. papatasi density was 63.33 specimens/m2/night in 1340 m altitude, whereas, P. sergenti was peaked with a density of 46.94 specimens/m2/night in 920 m altitude. The densities of males and females P. papatasi showed a negative correlation (r = −0.192, and −0.171, respectively) with altitude. However; P. sergenti was positively correlated with the increase in altitude but none of the correlations with altitude were significant for both species (P > 0.05; Table 2).
Finally, the correlation between densities of P. papatasi and P. sergenti and soil variables were also examined. Soil PH was negatively correlated with the density of female P. papatasi (r = −0.412, P = 0.004), while, soil water stress was negatively correlated with the male (r = −0.399, P = 0.005) and female (r = −0.372, P = 0.009) densities of P. papatasi (Table 2).
The main objective of this study was raised in response to our previous observations for site specific differences in sand fly collections. These observations led to a conclusion that environmental conditions available for each species could be different across Morocco which is exposed to continuous climate changes because of its unique location [15, 20]. We developed our current study to assess the relationship between the environmental conditions and the distribution of sand flies across endemic areas of CL. We also tried to examine the similarity between habitats in regard to the sand fly abundance. These results should be of importance to understand the environmental risk factors associated with high prevalence and establishment of vector populations for better understanding and guiding the target control programs in the country.
In this study, we included most CL endemic sites across the country to update sand fly fauna in all endemic sites and build a more comprehensive database for the abundance of two major vector species; P. sergenti responsible for transmission of ACL and P. papatasi responsible for transmission of ZCL among hosts. ACL is widespread in Guelmim-Es Smara, Marrakech-Tensift-Al Haouz, Souss-Massa-Draa and Tadla-Azilal regions such as Azilal, and Al Haouz and Chichaoua [22, 38, 39]. ZCL extended to the Saharan and pre-Saharan area of Guelmim-Es Smara, Meknes-Tafilalete and Souss-Massa-Draa regions such as Tata, Ouarzazate and Errachidia [4, 10]. The current study identified 11 species of sand flies in all habitats sampled; all of these species were previously identified in several other studies in Morocco [20–22, 31]. Sex ratios revealed more males in our collections; this result was not surprising given that the collection used sticky traps and occurred near the breeding sites of the sand flies. Sex ratios were also identified in favor of males in previous study in the region . In the study area, P. sergenti and P. papatasi co-existed and shared equally the territory (51.41 and 48.58 % for P. papatasi and P. sergenti, respectively). In spite of this co-occurrence, their abundance and densities varied according to several ecological factors.
Previous studies revealed that adult sand flies could not survive outside the temperature range of 10–40 °C  while reproduction is not possible below 15 °C . At around 30–32 °C, a significant increase in the vectorial capacity of sand flies was observed owing to shortening of the incubation period, despite a decrease in the vector’s survival . In arid areas of Morocco, P. papatasi is most active in the hot, dry season and it was the most abundant species when ambient temperature is in the 32–36 °C range . In southwest Asia, P. papatasi is most abundant in areas with a mean minimum temperature of 16 °C and mean maximum temperature of 44 °C from May to October .
In the present study, P. papatasi abundance seemed to be controlled mainly by the minimum temperature. Previous studies showed that ZCL (transmitted by P. papatasi) is prevalent in the pre-Saharan zones of North Africa where incidence and seasonal dynamic appeared to be controlled by minimum temperatures . These similarities revealed the importance of minimum temperature to limit the distribution of both vector and disease spread in Morocco. This study noted also the significance of only night-time temperatures with P. papatasi abundance which was not surprising considering that sand flies have nocturnal habit and are found active especially during night.
One of the limiting factors for the distribution of disease vectors is vegetation. For example, sand fly males are attracted to plants as a source of sugar meals. Plants also provide suitable resting, breeding and mating sites . Greenness, expressed in Normalized Vegetation index (NDVI), influences sand fly spatial repartition . Examining sand fly density according to NDVI index may help with understanding the relationships between these species and vegetation. Here, NDVI index seemed to be a key determinant for P. sergenti distribution and abundance, especially females which were positively sensitive to NDVI variations.
In our study area, precipitation was not a determining factor for presence and abundance of these two potential CL vectors; the study was conducted in dry season and few changes in precipitation were reported. A negative correlation was reported in the study area between precipitation and density of both P. papatasi and P. sergenti; this observation reflected the high abundance of the same species in low precipitation reported early by another report in Morocco . In contrast, aridity index, expressed as a function of rainfall and temperature, is one of the key factors negatively influencing the species density especially in the case of P. papatasi density. In Morocco, P. papatasi was found well adapted to arid climate conditions [14, 45], whereas, P. sergenti density was increased significantly in the arid and Saharan zones [14, 46]. Our results can also explain the absence of CL in the extreme southern Morocco which is an area characterized by higher aridity.
Many authors highlighted the role of soil conditions in spatial repartition of sand flies and also for characterization of sand fly breeding sites [47, 48]. Previous studies revealed that a combination of soil conditions, vegetation and topography create microclimates and ecological niches for P. papatasi . Soil moisture was identified as an important limiting factor for the distribution of phlebotomine sand flies including P. papatasi [49, 50]. In the present study, soil water stress influences negatively both male and female P. papatasi densities. Soil water stress is the monthly fraction of soil water content available for evapotranspiration and is inversely proportional to maximum soil moisture content, so a decrease in the soil water stress is caused by high soil moisture which is favorable conditions for an increase in the sand fly abundance. Similar observations for this increase in the sand fly densities were observed in other countries . This latter study identified a positive correlation between soil moisture and P. papatasi density. Schlein et al.  revealed that scarcity of water was the key-factor that limit the sand fly abundance and spread of leishmaniasis across Israel. However, the available sources of water moisture may be different between endemic sites across the Middle East including Morocco but soil moisture remain important for maintaining P. papatasi populations in all habitats including arid desert areas where irrigation activities occur.
The combination of the effect of several environmental factors determines the vector distribution. Our results revealed some of these variables and support the previous studies regarding the importance of soil moisture as one of the key factors that limit the distribution of the sand flies across the region. So, any ecological studies that focus on mapping vector species should consider these variables as important parameters in their studies. This study reflects a huge sampling across Morocco but we plan to further extend our current studies to include two more dimensions as our first priorities; (1) to include several other species associated with other forms of leishmaniasis, (2) to include long-term surveillance which examines the contribution of these variables to the vector species and disease dynamics for detailed insights on the long-term maintenance of disease hotspots across Morocco. Finally, one of the most important further objectives of this study is to explore mapping exercises and habitats suitability of CL vectors across the country to guide our surveillance programs.
The results of the present study suggest preliminary landmarks of vector distribution and CL risk in Morocco. The climatic factors (night-time temperature and aridity) with soil conditions and NDVI index are important factors to anticipate areas of disease risk and characterize areas of surveillance priorities for ZCL and ACL to guide the efforts for successful vector surveillance and control programs.
World Health Organization (WHO). Control of the Leishmaniasis: Report of a Meeting of the WHO Expert Committee on the Control of Leishmaniases. World Health Organization Technical Report Series No 949. Geneva: WHO; 2010.
Mahjour J, Akalay O, Saddiki A. Les leishmanioses au Maroc de l’analyse écoépidémiologie à la prévention. DEPS Bul Epidém. 1992; Supp 7.1-20
Rhajaoui M. Human leishmaniases in Morocco: a nosogeographical diversity. Pathol Biol. 2011;59(4):226–9.
Kahime K, Boussaa S, Bounoua L, Fouad O, Messouli M, Boumezzough A. Leishmaniasis in Morocco: diseases and vectors. Asian Pac J Trop Dis. 2014;4(2):S530–4.
Moroccan Ministry of Health: A report on progress of control programs against parasitic diseases. Directorate of Epidemiology and Disease Control, Ministry of Health, Rabat, Morocco; 2014, Available at http://www.sante.gov.ma/Publications/Etudes_enquete/Pages/default.aspx.
Kahime K, Boussaa S, Idrissi AL-E, Nhammi H, Boumezzough A. Epidemiological study on acute cutaneous leishmaniasis in Morocco. J Acute Dis. 2016;5(1):41–45.
Rioux JA, Petter F, Akalay O, Lanotte G, Ouazani A, Seguignes M, et al. Meriones shawi (Duvernoy, 1842) (Rodentia, Gerbillidae), réservoir de Leishmania major Yakimoff et Shokhor, 1914 dans le Sud Marocain. C.R. Acad Sci Paris. 1982;294:515–7.
Shehata MG, Samy AM, Doha SA, Fahmy AR, Kaldas RM, Furman BD, et al. First report of Leishmania tropica from a classical focus of L. major in North-Sinai, Egypt. Am J Trop Med Hyg. 2009;81(2):213–8.
Jacobson RL. Leishmania tropica (Kinetoplastida: Trypanosomatidae)--a perplexing parasite. Folia Parasitol. 2003;50(4):241–50.
Rioux JA, Guilvard E, Dereure J, Lanotte G, Denial M, Pratlong F, et al. Infestation naturelle de Phlebotomus papatasi (Scopoli, 1786) par Leishmania major MON-25. A propos de 28 souches isolées dans un foyer du Sud Marocain. In: Leishmania. Taxinomie et Phylogenèse Applications écoépidémiologiques (Coll Int CNRS/INSERM, 1984). Montpellier: IMEEE; 1986. p. 471–80.
Ajaoud M, Es-sette N, Hamdi S, El-Idrissi AL, Riyad M, Lemrani M. Detection and molecular typing of Leishmania tropica from Phlebotomus sergenti and lesions of cutaneous leishmaniasis in an emerging focus of Morocco. Parasit Vectors. 2013;6:217.
Ready PD. Leishmania manipulates sandfly feeding to enhance its transmission. Trends Parasitol. 2008;24(4):151–3.
Altizer S, Ostfeld RS, Johnson PT, Kutz S, Harvell CD. Climate change and infectious diseases: from evidence to a predictive framework. Science. 2013;341(6145):514–9.
Rioux J-A, Rispail P, Lanotte G, Lepart J. Relations Phlébotomes-bioclimats en écologie des leishmanioses Corollaires épidémiologiques. L’exemple du Maroc. B SOC BOT FR-ACTUAL. 1984;131(2–4):549–57.
Rispail P, Dereure J, Jarry D. Risk zones of human Leishmaniases in the Western Mediterranean basin: correlations between vector sand flies, bioclimatology and phytosociology. Mem Inst Oswaldo Cruz. 2002;97(4):477–83.
Bounoua L, Kahime K, Houti L, Blakey T, Ebi KL, Zhang P, et al. Linking climate to incidence of zoonotic cutaneous leishmaniasis (L. major) in pre-Saharan North Africa. Int J Environ Res Public Health. 2013;10(8):3172–91.
Samy AM, Campbell LP, Peterson AT. Leishmaniasis transmission: distribution and coarse-resolution ecology of two vectors and two parasites in Egypt. Rev Soc Bras Med Trop. 2014;47(1):57–62.
Samy AM, Doha SA, Kenawy MA. Ecology of cutaneous leishmaniasis in Sinai: linking parasites, vectors and hosts. Mem Inst Oswaldo Cruz. 2014;109(3):299–306.
High Planning Commission of Morocco (HCP). Report. 2014, Available at http://goo.gl/6jquVC.
Boussaa S, Guernaoui S, Pesson B, Boumezzough A. Seasonal fluctuations of phlebotomine sand fly populations (Diptera: Psychodidae) in the urban area of Marrakech. Morocco Acta Trop. 2005;95(2):86–91.
Kahime K, Boussaa S, Ouanaimi F, Boumezzough A. Species composition of phlebotomine sand fly fauna in an area with sporadic cases of Leishmania infantum human visceral leishmaniasis. Morocco Acta Trop. 2015;148:58–65.
Ramaoui K, Guernaoui S, Boumezzough A. Entomological and epidemiological study of a new focus of cutaneous leishmaniasis in Morocco. Parasitol Res. 2008;103(4):859–63.
Abonnenc E. Les phlébotomes de la région éthiopienne (Diptera: Phlebotomidae). Mémoire de l’ORSTOM. 1972;55:1–289.
Upton MS. Aqueous gum-chloral slide mounting media: an historical review. Bull Entomol Res. 1993;83(02):267–74.
Moroccan Ministry of Health. Fight against Leishmaniasis: Activity Guide. Rabat: Ministryof Health; 1997.
Moroccan Ministry of Health. Fight against Leishmaniasis: Activity Guide. Rabat: Ministryof Health; 2010.
Depaquit J, Leger N, Ferte H. The taxonomic status of Phlebotomus sergenti Parrot, 1917, vector of Leishmania tropica (Wright, 1903) and Phlebotomus similis Perfiliev, 1963 (Diptera - Psychodidae). Morphologic and morphometric approaches. Biogeographical and epidemiological corollaries. Bull Soc Pathol Exot. 1998;91(4):346–52.
Pesson B, Ready JS, Benabdennbi I, Martin-Sanchez J, Esseghir S, Cadi-Soussi M, et al. Sandflies of the Phlebotomus perniciosus complex: mitochondrial introgression and a new sibling species of P. longicuspis in the Moroccan Rif. Med Vet Entomol. 2004;18(1):25–37.
Leger N, Pesson B, Madulo-Leblond G, Abonnenc E. Differentiation of females of the subgenus Larroussius Nitzulescu 1931 (Diptera-Phlebotomidae) of the Mediterranean region. Ann Parasitol Hum Comp. 1983;58(6):611–23.
Benabdennbi I, Pesson B, Cadi-Soussi M, Morillas Marquez F. Morphological and isoenzymatic differentiation of sympatric populations of Phlebotomus perniciosus and Phlebotomus longicuspis (Diptera: Psychodidae) in northern Morocco. J Med Entomol. 1999;36(1):116–20.
Guernaoui S, Pesson B, Boumezzough A, Pichon G. Distribution of phlebotomine sandflies, of the subgenus Larroussius, in Morocco. Med Vet Entomol. 2005;19(1):111–5.
Boussaa S, Boumezzough A, Remy PE, Glasser N, Pesson B. Morphological and isoenzymatic differentiation of Phlebotomus perniciosus and Phlebotomus longicuspis (Diptera: Psychodidae) in Southern Morocco. Acta Trop. 2008;106(3):184–9.
McAleece N, Gage JDG, Lambshead PJD, Paterson GLJ. Biodiversity professional. Beta-Version. London: The Natural History Museum and the Scottish Association for Marine Sciences; 1997.
Samy AM, van de Sande WW, Fahal AH, Peterson AT. Mapping the potential risk of mycetoma infection in Sudan and South Sudan using ecological niche modeling. PLoS Negl Trop Dis. 2014;8(10):e3250.
Lewis DJ. Phlebotomid sandflies. Bull World Health Organ. 1971;44(4):535–51.
Kassem HA, El-Sayed YA, Baz MM, Kenawy MA, El Sawaf BM. Climatic factors influencing the abundance of Phlebotomus papatasi (Scopoli) (Diptera: Psychodidae) in the Nile Delta. J Egypt Soc Parasitol. 2009;39(1):305–16.
Minitab Statistical Software: State College, PA: Minitab, Inc. 2010, Available at www.minitab.com.
Pratlong F, Rioux JA, Dereure J, Mahjour J, Gallego M, Guilvard E, et al. Leishmania tropica in Morocco. IV--Intrafocal enzyme diversity. Ann Parasitol Hum Comp. 1991;66(3):100–4.
Rhajaoui M, Sebti F, Fellah H, Alam MZ, Nasereddin A, Abbasi I, et al. Identification of the causative agent of cutaneous leishmaniasis in Chichaoua province. Morocco Parasite. 2012;19(1):81–4.
Ouanaimi F, Boussaa S, Boumezzough A. Phlebotomine sand flies (Diptera: Psychodidae) of Morocco: results of an entomological survey along three transects from northern to southern country. Asian Pac J Trop Dis. 2015;5(4):299–306.
Cross ER, Hyams KC. The potential effect of global warming on the geographic and seasonal distribution of Phlebotomus papatasi in southwest Asia. Environ Health Perspect. 1996;104(7):724–7.
Kasap OE, Alten B. Comparative demography of the sand fly Phlebotomus papatasi (Diptera: Psychodidae) at constant temperatures. J Vector Ecol. 2006;31(2):378–85.
Benkova I, Volf P. Effect of temperature on metabolism of Phlebotomus papatasi (Diptera: Psychodidae). J Med Entomol. 2007;44(1):150–4.
Elnaiem DA, Hassan HK, Ward RD. Associations of Phlebotomus orientalis and other sandflies with vegetation types in the eastern Sudan focus of kala-azar. Med Vet Entomol. 1999;13(2):198–203.
Rioux JA, Akalay O, Périères J, Dereure J, Mahjour J, Le Houérou HN, et al. Evolution of Eco-epidemiological of ‘leishmaniasis risk’ at Atlantic Moroccan Sahara. Heuristic interest of the relationship ‘sandflies bioclimates’. Ecol Mediter. 1997;23:73–92.
Guernaoui S, Boumezzough A. Habitat preferences of phlebotomine sand flies (Diptera: Psychodidae) in southwestern Morocco. J Med Entomol. 2009;46(5):1187–94.
Bettini S, Melis P. Leishmaniasis in Sardinia. III. Soil analysis of a breeding site of three species of sandflies. Med Vet Entomol. 1988;2(1):67–71.
Moncaz A, Faiman R, Kirstein O, Warburg A. Breeding sites of Phlebotomus sergenti, the sand fly vector of cutaneous leishmaniasis in the Judean Desert. PLoS Negl Trop Dis. 2012;6(7):e1725.
Schlein Y, Warburg A, Schnur LF, Le Blancq SM, Gunders AE. Leishmaniasis in Israel: reservoir hosts, sandfly vectors and leishmanial strains in the Negev, Central Arava and along the Dead Sea. Trans R Soc Trop Med Hyg. 1984;78(4):480–4.
Yuval B. Populations of Phlebotomus papatasi (Diptera: Psychodidae) and the risk of Leishmania major transmission in three Jordan Valley habitats. J Med Entomol. 1991;28(4):492–5.
Wasserberg G, Abramsky Z, Anders G, El-Fari M, Schoenian G, Schnur L, et al. The ecology of cutaneous leishmaniasis in Nizzana, Israel: infection patterns in the reservoir host, and epidemiological implications. Int J Parasitol. 2002;32(2):133–43.
The authors would like to thank the staff members of the Laboratory of Ecology and Environment, Faculty of Sciences, Semlalia, Cadi Ayyad University, Morocco for hosting and supporting the study. Thanks also to the Department of Entomology, and Research and Training Center on Vectors of Diseases of Ain Shams University, Cairo, Egypt for their support to AMS during our current study.
The authors declare that they have no competing interests.
SB & KK conceived and designed the study for the sand fly collection, processing and identification. KK, AB, SB, AMS and AB analyzed the data and drafted the manuscript. All authors reviewed and approved the final manuscript.
Detailed data set for sand fly species sampled in 48 sites across Morocco from April to June 2013. The specimens were sorted as male or female for each species. The coordinates (longitude and latitude) for each site is also presented in this data set. (XLSX 23 kb)
Species evenness across 48 sampling sites in Morocco. The Species evenness is calculated using the BioDiversity Professional statistics analysis software . (DOCX 14 kb)
Bray-Curtis similarity based on the sand fly abundance across all sampling sites in Morocco. The scale represents the similarity between the sampling sites regarding the sand fly abundance. (DOC 115 kb)