- Open Access
Ecological niche modeling predicting the potential distribution of Leishmania vectors in the Mediterranean basin: impact of climate change
Parasites & Vectorsvolume 11, Article number: 461 (2018)
Due to climate change, the geographical distribution of sand flies during the last decades has shifted northward from latitudes below 45°N in southern Europe to latitudes just above 50○N. Recent studies show that some phlebotomine sand flies were recorded in several parts of Germany and Belgium. In central Europe, some autochthone leishmaniasis cases are being recorded in regions traditionally regarded as leishmaniasis-free. An important challenge is to predict the geographical distribution of leishmaniasis vectors under new climatic conditions. In this study, we attempted to predict the current distribution of six leishmaniasis vectors in the Mediterranean basin and forecast species’ geographical shift under future climate scenarios using an ensemble ecological niche modeling approach. Species records were obtained from scientific surveys published in the research literature between 2006 and 2016. A series of climate metrics describing temperature and precipitation in the study area under two climatic scenarios were obtained from WorldClim database. A consensus model was derived from six varieties of modeling approaches (regression, machine learning and classification techniques) in order to ensure valid prediction of distribution of vectors under different climate scenarios.
Model performance was generally high for the included species with a specificity (true negative rate) ranging from 81.03 to 96.52% (mean = 86.94%) and a sensitivity (true positive rate) ranging from 87.93 to 100% (mean = 96.98%). Our work evidenced the hypothesis of the widespread of Leishmania vectors under climate change scenarios. All of the studied species are prospected to gain new areas that are actually not suitable for vectors’ survival. Phlebotomine sand flies are prospected to invade extra-Mediterranean regions, especially western and central Europe.
Our study confirmed the importance of environmental and climate factors on the distribution of leishmaniasis vectors and demonstrated the performance of ecological niche modeling in the prediction of the geographical spread of vector-borne diseases. Ecological niche modeling should be considered in the future as a valuable tool in addition to experimental laboratory studies for a better understanding of the biology of vector species.
Global warming and climate change due to industrialization, intensive farming and anthropogenic activities leading to higher greenhouse gases emission, has become a fact. The direct consequences of increasing the greenhouse effect are changes in temperature, precipitation patterns and other climate variables .
The Mediterranean basin is considered one of the most vulnerable regions of the world to climate change . Indeed, the Mediterranean region has shown large climate variation in the past  and has been identified as one of the most sensitive areas according to future climate change projections [4, 5]. According to the Intergovernmental Panel on Climate Change  temperature will increase in the Mediterranean area while annual precipitation will decrease during the 21st century. Indeed, annual rainfall will decrease in Mediterranean Africa and the northern Sahara, with a greater decrease as the Mediterranean coast is approached. Precipitation will increase in most of northern Europe and decrease in most southern countries bordering the Mediterranean Sea. The same projection for the 21st century states that the annual mean warming will vary from 2.2 to 5.1 °C in southern Europe and the Mediterranean area, leading to warmer summers and more temperate winters.
In addition to its direct effect on human health and welfare, consisting of excessive heat exposure and extreme weather events (storms, floods and droughts), climate change may also affect the spread of vector-borne diseases and other infectious diseases [6,7,8].
The environmental factors are the main causes of the emergence and re-emergence of infectious diseases. Without a doubt, environmental change can create favorable conditions for the proliferation of reservoirs, increase the density of vectors, and consolidate the reservoir-vector-host interaction . Changes in climatic conditions (temperature, rainfall and humidity) that are prospected to occur under different climate change scenarios will affect the geographical distribution of vector-borne disease .
Leishmaniasis is a vector-borne infection caused by parasitic protozoans of the genus Leishmania that are transmitted to humans through the bite of an infected female phlebotomine sand fly. Among over 800 species of sand fly, 98 are proven or suspected vectors of human leishmaniases; these include 42 Phlebotomus species in the Old World and 56 Lutzomyia species in the New World . In countries bordering the Mediterranean Sea, including those in southern Europe, North Africa and the Middle East, more than 20 species of Phlebotomus are widely distributed .
The optimal temperature for the development of sand flies and Leishmania parasites is approximately 25 °C [13, 14]. The prospected change in ambient temperature would increase the vector biting rate and number of human exposures, and reduce the incubation time of the infective agent within its vector [15, 16]. Furthermore, a higher temperature will reduce the vector winter mortality and new areas will become environmentally suitable for its survival and reproduction. There is evidence that the geographical range of phlebotomine sand flies incriminated in leishmaniasis transmission has changed in response to climate change in the Mediterranean region and Europe.
The geographical distribution of sand flies during the last decades has shifted northward from latitudes below 45°N in southern Europe  to latitudes just above 50°N . Recent studies show that some phlebotomine sand flies were recorded in several parts of Germany and Belgium [18, 19]. In central Europe, some autochthone leishmaniasis cases are being recorded in a region traditionally regarded as leishmaniasis-free [20, 21].
Recently, species distribution modeling has acquired importance, not only for the ability to predict species’ geographical distribution from explanatory variables, but also for the faculty to forecast species’ occurrence under new areas and climatic conditions . However, selecting the appropriate model is challenging not only for the various techniques available but also for the different approaches adopted by each model which yield to different results . One alternative to address inter-model variations is to combine several models. The use of an ensemble approach can result in significant improvements on the robustness of predicting species distribution and forecasting their geographical shift under climate change projections .
The main motivation of this study was to (i) predict the current distribution of leishmaniasis vectors in the Mediterranean basin; (ii) forecast species’ geographical shifts under future climate scenarios; and (iii) evaluate climatic variables leading to the geographical spread or extinction of these vectors using an ensemble modeling approach.
Species records were obtained from scientific surveys published in research literature between 2006 and 2016 [24,25,26,27,28,29,30,31,32,33,34,35,36,37]. We performed a literature search in the SCOPUS database using search topic = Phlebotomus OR Sand flies OR Phlebotominae AND Geographical coordinate. Literature indicating only the name of the city or village where the species was recorded were discarded and only species’ occurrences with explicit geographical coordinates were retained. Since the values of environmental predictors at these locations were extracted to fit the model prediction of the species presence across the study area, integrating accurate geographical occurrence points is crucial for better model fitting.
We found a total of 384 data-points distributed in 16 Mediterranean countries (Albania, Algeria, Cyprus, Egypt, France, Greece, Israel, Italy, Lebanon, Libya, Montenegro, Morocco, Portugal, Spain, Tunisia and Turkey) and six phlebotomine sand fly species (27 for Phlebotomus alexandri; 66 for Phlebotomus ariasi; 25 for Phlebotomus kazeruni; 179 for Phlebotomus papatasi; 32 for Phlebotomus perniciosus; and 55 for Phlebotomus sergenti).
A series of climate metrics, describing temperature and precipitation in the study area under two climatic scenarios, were obtained from the WorldClim database. These climate metrics are derived from monthly temperature and rainfall values and represent biologically meaningful variables for characterizing species distribution .
The WorldClim data layers included 11 temperature and eight precipitation metrics, expressing spatial variations in annual means, seasonality, and extreme or limiting climatic factors with a spatial resolution of about one square kilometer. The data layers were generated through interpolation of average monthly climate data from 4000 weather stations. Climate layers were available for long-term time series (1950–2000), future projections and past conditions.
Climate change projections were simulated by the Canadian Centre for Climate Modelling and Analysis  for 2020, 2050 and 2080 under a pessimistic (A2a) and an optimistic (B2a) story-line scenario. A2a describes a highly heterogeneous future world with regionally oriented economies. The main driving forces are a high rate of population growth, increased energy use, land-use changes and slow technological change. The B2a scenario is also regionally oriented but with a general evolution towards environmental protection and social equity. Compared to A2a, B2a has a lower rate of population growth and a smaller increase in gross domestic product, but more diverse technological changes and slower land-use changes .
The original 19 predictor variables were reduced to only 5 climate variables by analysis of collinearity . A correlation matrix between the 19 climate variables was performed using ArcGis 10.1 and a pair of variables with a correlation coefficient greater than 0.9 were considered as highly correlated. For highly correlated variables, when possible, we preferred extreme variables (min, max) to mean variables since the biological behavior of the vector is highly affected by extreme climatic conditions. The final dataset included: maximum temperature of warmest month, minimum temperature of coldest month, precipitation of wettest quarter, precipitation of driest quarter and elevation.
Species distribution model
We used six individual models contributing to the consensus model including two regression techniques (generalized linear models, GLM; generalized additive models, GAM), two machine learning techniques (artificial neural networks, ANNs; maximum entropy, MaxEnt) and two classification techniques (surface range envelop, SRE; random forest, RF).
Generalized linear models  are extensions of traditional linear regression models that allow the mean to depend on the explanatory variables through a link function, and the response variable to be any member of a set of distributions called the exponential family, which includes normal, Poisson and binomial distributions, among others.
Generalized additive models  blend properties of generalized linear models with additive models. The purpose of generalized additive models is to maximize the quality of prediction of a dependent variable from various distributions, by estimating unspecific (non-parametric) functions of the predictor variables, which are “connected” to the dependent variable via a link function.
Artificial neural networks  are a machine learning approach. It consists of creating new linear combinations of the input layers called hidden layers; later the created hidden layers are related to a response layer via non-linear function .
MaxEnt  is also a machine learning method that consists in finding the probability distribution of maximum entropy subject to constraints imposed by information derived from the values of the predictor variables in the occurrence points.
The surface range envelop  consists of identifying maximum and minimum values for each input variable from the set of occurrence points. Any pixel with all variables falling between these maximum and minimum limits is included within the range.
Random forest  is an ensemble classifier algorithm using many decision tree models, which can be used for classification or regression. Different subsets of the training data are selected, with replacement, to train each tree. The remaining training data are used to estimate variables’ error and importance. Class assignment is made by the number of votes from all of the trees and for regression the average of the results is used.
We combined predictions into an “ensemble” by averaging the models with a true skill statistic  higher than 0.8 to obtain a “consensus model” and to avoid the integration of weak models. The main aim of consensus methods is to decrease the predictive uncertainty of single-models, since predictions of species distributions can vary widely among modeling approaches. Indeed, ensemble forecasting can enable a more robust model and overcome the uncertainties derived from each individual model .
True skill statistic (TSS) includes omission and commission errors. It ranges from -1 to +1, where +1 indicates perfect classification and values below zero indicate a performance no better than random. We used TSS to evaluate the potential species ranges because it has been shown to be insensitive to prevalence and is unaffected by the size of the validation set .
As most of the used algorithms need presence/absence data for the model fitting and evaluation, we randomly selected pseudo-absence points in locations where the studied species was not present but the other remaining sand fly species had been recorded, as suggested by Phillips et al.  and Mateo et al. .
We ran for each model a 5-fold cross-validation by randomly splitting our occurrence data into 70% used for the model calibration and 30% used to evaluate the current predictive performance of the models. This 5-fold cross-validation yielded an average model for each species and algorithm. Calibrated models, for current conditions, were then used to generate projections of future potential distributions for 2020, 2050 and 2080 under the two climate change scenarios, A2a and B2a. For species with no model presenting a TSS score higher than 0.8, no consensus projection was created and the species were removed from further analyses. All analyses were carried out using the Biomod2 package  for R .
Referring to the mean TSS score (0.84 ± 0.023), model performance was generally high for the included species with a specificity (true negative rate) ranging from 81.03 to 96.52% (mean = 86.94%) and a sensitivity (true positive rate) ranging from 87.93 to 100% (mean = 96.98%). Precipitation of the wettest quarter (Bio16) and precipitation of the driest quarter (Bio17) were the most contributing variables to the models for P. papatasi, P. ariasi and P. alexandri with 44.72, 43.96 and 64.46%, respectively, of explained variation (see Fig. 1). On the other hand, the minimum temperature of the coldest month (Bio6) and the maximum temperature of the warmest month (Bio5) had the highest contribution to the P. sergenti model with 45.86% of explained variation. The contribution of the altitude ranged from 10.59% for the P. ariasi model to 36.26% for the P. alexandri model, with a mean contribution among all modeled species of 23.64%.
Suitable areas for P. papatasi presence are predicted occur in North Africa and the Middle East. Lower occurrence was predicted in the south of most European countries bordering the Mediterranean Sea. The areas suitable for P. papatasi range from a latitude below 19°N in the African Sahara to a latitude above 48°N in southern Europe.
Climate change projection for P. papatasi showed that the species will shift northward, with a centroid shift estimated at 292, 595 and 948 km for the time periods 2020, 2050 and 2080, respectively, under the A2a (pessimistic) scenario (Table 1). Similar trends were observed for the B2a scenario. However, the shift of the centroid of the suitable area for P. papatasi occurrence was lower, with a maximum shift of 728 km northward suspected to occur in 2080.
Under the two climate change scenarios (A2a and B2a) and the three time periods (2020, 2050 and 2080) the suitable area for P. papatasi presence was predicted to be gained progressively in northern and western Europe. Concurrently, the potential area for species presence would be lost in the south of north African countries and the Middle East (Fig. 2). Nevertheless, the total gained area is higher than the total lost area for all scenarios and time periods (Table 1).
Ecological niche modeling for current conditions showed that the suitable area for P. alexandri occurrence is widely distributed from Morocco in the west, to Iran and Kazakhstan in the east of the study area. Phlebotomus alexandri is predicted to be less present in western and central Europe, where it was sporadically predicted in Spain, plains bordering the Alps chain and in some countries of central Europe (Fig. 2).
The model projection into the two future climate change scenarios showed that the suitable area for P. alexandri distribution will be lost in the north-east of the study area which will cause a centroid shift estimated at 314 and 304 km toward the south-west under A2a and B2a, respectively, in 2020. A higher shift in the same direction is predicted for 2050 and 2080 under both scenarios. For all time-periods and under the two climate change scenarios, the percent of lost area predicted in central and western Asia (especially in Russia, Kazakhstan, Turkey and Iran) was higher than the percent of gained area (Table 1), which was predicted to occur essentially in north African countries, the Middle East, southern Spain and France.
The area with suitable environmental conditions for P. ariasi occurrence under current conditions was suspected to occur in North Africa, western Europe and most southern European countries, bordering the Mediterranean Sea.
Under both climate change scenarios and the three time-periods, the species is suspected to progressively lose areas of suitable conditions in western Europe, especially in Spain and France, while new area will be gained by the species eastward including central Europe and western Asia. The area gained by the species is far higher than that predicted to be lost (Fig. 3). Although the two climate change scenarios showed a similar trend, a greater change in the distribution of the lost and gained area by the species is predicted under the A2a scenario, especially during the time-period 2020–2050 (gained area 69.39% under A2a, 37.46% under B2a) and 2050–2080 (lost area -37.60% under A2a, -28.70% under B2a). The change in the suitable area for P. ariasi distribution is reflected in the shift of its centroid projected to move eastward with a maximum shift in 2080 estimated to reach 1279 and 1150 km under the A2a and B2a climate change scenarios, respectively (Table 1).
Ecological niche modeling under current climatic conditions showed that suitable area for P. sergenti is widely distributed. It ranges from a latitude of 18°N on the south border of the north African countries to a latitude of 51°N in the south of the UK. Nevertheless, suitable areas are not uniformly distributed; the species occurs mainly between 24–48°N in northern Africa, southern Europe, central Europe, the Middle East and western Asia.
All future climate change projections showed that gained area will occur in northern Europe and in countries on the south border of north African countries. These regions actually present a reduced suitable area for P. sergenti occurrence which will expand progressively during future time-periods, while large areas currently predicted as suitable for the vector presence will be lost, especially in north African countries, southern Europe, the Middle East and western Asia (Fig. 3).
Despite the large change in the suitable area for P. sergenti prospected to occur, confirmed by the gained and lost area by the species under the two climate change scenarios, the centroid shift will reach only a maximum shift of 290 km in 2050 under A2a (Table 1).
Ecological niche modeling has allowed the prediction of the potential distribution of Leishmania vectors in the Mediterranean basin. Our study provides evidence for the hypothesis of the widespread Leishmania vectors under climate change scenarios. All the studied species are prospected to gain new areas that are actually not suitable for the vectors’ survival.
Indeed, our study corroborates other authors’ findings [12, 52, 53] suggesting that the increase of temperature and humidity, due to the global warming prospected to occur during the 21st century, will affect the distribution of Leishmania vectors. As suggested by previous studies  phlebotomine sand flies are prospected to invade extra-Mediterranean regions especially western and central Europe. Concurrently, Leishmania vectors will lose areas for suitable conditions in North Africa and Middle East. These changes in their geographical distribution are more intense under the A2a scenario than the B2a. Similar projections for the distribution of P. papatasi, P. ariasi and P. sergenti were reported by Trajer et al. . Nevertheless, these authors focused their study only on the distribution of phlebotomine species in Europe, since no presence records in North Africa and Middle East were included.
In a non-endemic location, where the species is suspected to occur, special attention should be taken. Field surveys to detect the presence of Leishmania vectors and monitoring cases of disease in these regions, is crucial to prevent the emergence of new disease foci. Additionally, domestic pets (especially dogs) travelling between endemic areas and non-endemic regions suspected to be suitable for the vector disease should be controlled to avoid the risk of transmission. Indeed, in non-endemic regions where the vector is suspected to be present, the introduction of infected dogs will contribute to the establishment of the disease cycle.
Our study focused only on the fundamental ecological niche, which represents the abiotic conditions for the species survival. However, the species is in fact present in a restricted part of the fundamental niche due to biotic interactions. Actually, Hutchinson  defined the realized niche as the portion of the fundamental niche from which a species is not excluded due to biotic competition. Furthermore, the dispersal ability of the species is limited by the maximum flight distance of phlebotomine sand flies, estimated as 1000 m , and geographical barriers. Indeed, even if the wind plays a major role in the species dissemination, the presence of mountain chains and water bodies can limit their geographical dispersal to new suitable areas for the species’ survival and reproduction .
Climate change is projected to play an essential role in the expansion of the geographical distribution of specific species of leishmaniasis vectors. Sand fly species are prospected to invade extra-Mediterranean regions, especially in Europe. The monitoring of regions presenting suitable conditions for the vectors’ presence is highly recommended to prevent the outbreak of the disease in densely populated areas with a naïve population. Our study confirmed the importance of environmental and climate factors on the distribution of leishmaniasis vectors and demonstrated the performance of ecological niche modeling in the prediction of the geographical spread of vector diseases. Ecological niche modeling should be considered in the future as a valuable tool in addition to experimental laboratory studies for a better understanding of the biology of vector species.
Artificial neural networks
Precipitation of wettest quarter
Precipitation of driest quarter
Maximum temperature of warmest month
Minimum temperature of coldest month
Generalized additive models
Generalized linear models
Surface range envelop
True skill statistic
IPCC Fourth Assessment Report (AR4). In: Solomon S, Qin D, Manning M, Chen Z, Marquis M, Averyt KB, Tignor M, Miller HL, editors. Climate change 2007: The physical science basis. Contribution of Working Group I to the 4th Assessment Report of the Intergovernmental Panel on Climate Change. Cambridge: Cambridge University Press; 2007.
Schröter D, Cramer W, Leemans R, Prentice IC, Araújo MB, Arnell NW, et al. Ecosystem service supply and vulnerability to global change in Europe. Science. 2005;310:1333–7.
Luterbacher J, Xoplaki E, Press AIN, Casty C, Wanner H, Pauling A, et al. Mediterranean climate variability over the last centuries: a review. Dev Earth Environ Sci. 2006;4:27–148.
Giorgi F, Lionello P. Climate change projections for the Mediterranean region. Glob Planet Chang. 2008;63:90–104.
Giorgi F. Climate change hot-spots. Geophys Res Lett. 2006;33:L08707.
Koch LK, Kochmann J, Klimpel S, Cunze S. Modeling the climatic suitability of leishmaniasis vector species in Europe. Sci Rep. 2017;7:13325.
Carvalho BM, Rangel EF, Vale MM. Evaluation of the impacts of climate change on disease vectors through ecological niche modelling. Bull Entomol Res. 2017;107:419–30.
McIntyre S, Rangel EF, Ready PD, Carvalho BM. Species-specific ecological niche modelling predicts different range contractions for Lutzomyia intermedia and a related vector of Leishmania braziliensis following climate change in South America. Parasit Vectors. 2017;10:157.
Rogers DJ, Randolph SE. Climate change and vector-borne diseases. Adv Parasitol. 2006;62:345–81.
World Health Organization Regional Office for Europe. Floods: Climate change and adaptation strategies for human health. Darmstadt: Steinkopff Verlag; 2002.
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.
Aspöck H, Gerersdorfer T, Formayer H, Walochnik J. Sandflies and sandfly-borne infections of humans in central Europe in the light of climate change. Wien Klin Wochenschr. 2008;120:24–9.
Killick-Kendrick R, Killick-Kendrick M. The laboratory colonization of Phlebotomus ariasi (Diptera: Psychodidae). Ann Parasitol Hum Comp. 1987;62:354–6.
Rioux JA, Aboulker JP, Lanotte G, Killick-Kendrick R, Martini-Dumas A. Ecology of leishmaniasis in the south of France. 21. Influence of temperature on the development of Leishmania infantum Nicolle, 1908 in Phlebotomus ariasi Tonnoir, 1921. Experimental study. Ann Parasitol Hum Comp. 1985;60:221–9.
Bates PA. Leishmania sand fly interaction: progress and challenges. Curr Opin Microbiol. 2008;11:340–4.
Bates PA. Transmission of Leishmania metacyclic promastigotes by phlebotomine sand flies. Int J Parasitol. 2007;37:1097–106.
Martens WJM, AJ MM. Environmental Change, Climate and Health: Issues and Research Methods. Cambridge: Cambridge University Press; 2009.
Naucke TJ, Menn B, Massberg D, Lorentz S. Sandflies and leishmaniasis in Germany. Parasitol Res. 2008;103:S65–8.
Mencke N. The importance of canine leishmaniosis in non-endemic areas, with special emphasis on the situation in Germany. Berl Munch Tierarztl Wochenschr. 2011;124:434–42.
Gogoaşe MG, Teodorescu I, Preda C, Ionescu SC. Two case reports on visceral leishmaniasis diagnosed in Romania. Roum Arch Microbiol Immunol. 2013;72:49–62.
Tánczos B, Balogh N, Király L, Biksi I, Szeredi L, Gyurkovsky M, et al. First record of autochthonous canine leishmaniasis in Hungary. Vector Borne Zoonotic Dis. 2012;12:588–94.
Elith J, Graham CH, Anderson RP, Dudík M, Ferrier S, Guisan A, et al. Novel methods improve prediction of species’ distributions from occurrence data. Ecography. 2006;29:129–51.
Araújo MB, New M. Ensemble forecasting of species distributions. Trends Ecol Evol. 2007;22:42–7.
Abdel-Dayem MS, Annajar BB, Hanafi HA, Obenauer PJ. The potential distribution of Phlebotomus papatasi (Diptera: Psychodidae) in Libya based on ecological niche model. J Med Entomol. 2012;49:739–45.
Barón S, Martín-Sánchez J, Gállego M, Morales-Yuste M, Boussaa S, Morillas-Márquez F. Intraspecific variability (rDNA ITS and mtDNA Cyt b) of Phlebotomus sergenti in Spain and Morocco. Acta Trop. 2008;107:259–67.
Ballart C, Barón S, Alcover MM, Portús M, Gállego M. Distribution of phlebotomine sand flies (Diptera: Psychodidae) in Andorra: first finding of P. perniciosus and wide distribution of P. ariasi. Acta Trop. 2012;122:155–9.
Franco FAL, Morillas-Márquez F, Barón SD, Morales-Yuste M, Gálvez R, Díaz V, et al. Genetic structure of Phlebotomus (Larroussius) ariasi populations, the vector of Leishmania infantum in the western Mediterranean: epidemiological implications. Int J Parasitol. 2010;40:1335–46.
Otranto D, de Caprariis D, Lia RP, Tarallo V, Lorusso V, Testini G, et al. Prevention of endemic canine vector-borne diseases using imidacloprid 10% and permethrin 50% in young dogs: a longitudinal field study. Vet Parasitol. 2010;172:323–32.
Martín-Sánchez J, Gállego M, Barón S, Castillejo S, Morillas-Marquez F. Pool screen PCR for estimating the prevalence of Leishmania infantum infection in sandflies (Diptera: Nematocera, Phlebotomidae). Trans R Soc Trop Med Hyg. 2006;100:527–32.
Branco S, Alves-Pires C, Maia C, Cortes S, Cristovão JMS, Gonçalves L, et al. Entomological and ecological studies in a new potential zoonotic leishmaniasis focus in Torres Novas municipality, Central Region, Portugal. Acta Trop. 2013;125:339–48.
Chalghaf B, Chlif S, Mayala B, Ghawar W, Bettaieb J, Harrabi M, et al. Ecological niche modeling for the prediction of the geographic distribution of cutaneous leishmaniasis in Tunisia. Am J Trop Med Hyg. 2016;94:844–51.
Demir S, Gocmen B, Ozbel Y. Faunistic study of sand flies in northern Cyprus. North-West J Zool. 2010;6:149–61.
Depaquit J, Lienard E, Verzeaux-Griffon A, Ferté H, Bounamous A, Gantier J-C, et al. Molecular homogeneity in diverse geographical populations of Phlebotomus papatasi (Diptera, Psychodidae) inferred from ND4 mtDNA and ITS2 rDNA. Infect Genet Evol. 2008;8:159–70.
Faraj C, Ouahabi S, Adlaoui EB, El Elkohli M, Lakraa L, El Rhazi M, Ameur B. Insecticide susceptibility status of Phlebotomus (Paraphlebotomus) sergenti and Phlebotomus (Phlebotomus) papatasi in endemic foci of cutaneous leishmaniasis in Morocco. Parasit Vectors. 2012;5:51.
Kassem HA, Siri J, Kamal HA, Wilson ML. Environmental factors underlying spatial patterns of sand flies (Diptera: Psychodidae) associated with leishmaniasis in southern Sinai, Egypt. Acta Trop. 2012;123:8–15.
Khalid NM, Aboud MA, Alrabba FM, Elnaiem D-EA, Tripet F. Evidence for genetic differentiation at the microgeographic scale in Phlebotomus papatasi populations from Sudan. Parasit Vectors. 2012;5:249.
Prudhomme J, Gunay F, Rahola N, Ouanaimi F, Guernaoui S, Boumezzough A, et al. Wing size and shape variation of Phlebotomus papatasi (Diptera: Psychodidae) populations from the south and north slopes of the Atlas Mountains in Morocco. J Vector Ecol. 2012;37:137–47.
Saatchi S, Buermann W, ter Steege H, Mori S, Smith TB. Modeling distribution of Amazonian tree species and diversity using remote sensing measurements. Remote Sens Environ. 2008;112:2000–17.
Kim S-J, Flato G, Boer G. A coupled climate model simulation of the last glacial maximum, part 2: approach to equilibrium. Clim Dyn. 2003;20:635–61.
Franklin J. Mapping Species Distributions: Spatial Inference and Prediction. Cambridge: Cambridge University Press; 2010.
McCullagh P, Nelder JA. Generalized Linear Models. New York: Chapman and Hall; 1989.
Hastie T, Tibshirani R. Generalized Additive Models. London: Chapman and Hall; 1990.
McCulloch WS, Pitts W. A logical calculus of the ideas immanent in nervous activity. Bull Math Biophys. 1943;5:115–33.
Phillips SJ, Dudík M, Schapire RE. A maximum entropy approach to species distribution modeling. In: Proceedings of the 21st International Conference on Machine Learning ICML ’04, New York. New York: ACM Press; 2004. p. 83.
Busby JR. BIOCLIM -a bioclimatic analysis and predictive system. In: Margules CR, Austin MP, editors. Nature Conservation: Cost Effective Biological Surveys and Data Analysis. Canberra: CSIRO; 1991. p. 64–8.
Breiman L. Random forests. Mach Learn. 2001;45:5–32.
Allouche O, Tsoar A, Kadmon R. Assessing the accuracy of species distribution models: prevalence, kappa and the true skill statistic (TSS). J Appl Ecol. 2006;43:1223–32.
Phillips SJ, Dud M, Elith J, Graham CH, Lehmann A, Leathwick J, et al. Sample selection bias and presence-only distribution models: implications for background and pseudo-absence data. Ecol Appl. 2009;19:181–97.
Mateo RG, Croat TB, Felicísimo ÁM, Muñoz J. Profile or group discriminative techniques? Generating reliable species distribution models using pseudo-absences and target-group absences from natural history collections. Divers Distrib. 2010;16:84–94.
Thuiller W, Lafourcade B, Engler R, Araújo MB. BIOMOD - a platform for ensemble forecasting of species distributions. Ecography. 2009;32:369–73.
R Development Core Team. R: a language and environment for statistical computing. Vienna: R Foundation for Statistical Computing; 2008.
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:724–7.
Rajesh K, Sanjay K. Change in global climate and prevalence of visceral leishmaniasis. Int J Sci Res Publ. 2013;3:2250–3153.
Trajer AJ, Bede-Fazekas A, Hufnagel L, Horvath L, Bobvos J, Paldy A. The effect of climate change on the potential distribution of the European Phlebotomus species. Appl Ecol Environ Res. 2013;11:189–208.
Hutchinson GE. Concluding remarks. Cold Spring Harb Symp Quant Biol. 1957;22:415–27.
Killick-Kendrick R, Wilkes TJ, Bailly M, Bailly I, Righton LA. Preliminary field observations on the flight speed of a phlebotomine sandfly. Trans R Soc Trop Med Hyg. 1986;80:138–42.
Fischer D, Moeller P, Thomas SM, Naucke TJ, Beierkuhnlein C. Combining climatic projections and dispersal ability: a method for estimating the responses of sandfly vector species to climate change. PLoS Negl Trop Dis. 2011;5:e1407.
We thank the American Association for the Advancement of Science for its kind funding.
This study was funded under grant number S-LMAQM-12-GR-1145, provided by the American Association for the Advancement of Science (AAAS) and awarded to BC as part of the AAAS BMENA Scientific Exchange Programme, supported by a grant from the U.S. Department of State. The funders had no role in the study design, data collection and analysis, decision to publish, or preparation of the manuscript.
Availability of data and materials
All data generated or analyzed during this study are included in this article.
Ethics approval and consent to participate
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.