- Open Access
Factors influencing the presence of sand flies in Majorca (Balearic Islands, Spain) with special reference to Phlebotomus pernicious, vector of Leishmania infantum
Parasites & Vectors volume 7, Article number: 421 (2014)
Although the Mediterranean island of Majorca is an endemic area of leishmaniosis, there is a lack of up-to-date data on its sand fly fauna, the last report dating from 1989. The aim of the present study was to provide information on the current sand fly distribution, the potential environmental factors favoring the presence of Phlebotomus perniciosus and which areas are at risk of leishmaniosis.
In July 2008 sand fly captures were carried out in Majorca with sticky castor oil interception traps. The capture stations were distributed in 77 grids (5x5 km2) covering the entire island. A total of 1,882 sticky traps were set among 111 stations. The characteristics of the stations were recorded and maps were designed using ArcGIS 9.2 software. The statistical analysis was carried out using a bivariate and multivariate logistic regression model.
The sand fly fauna of Majorca is composed of 4 species: Phlebotomus perniciosus, P sergenti, P. papatasi and Sergentomyia minuta. P. perniciosus, responsible for Leishmania infantum transmission, was captured throughout the island (frequency 69.4 %), from 6 to 772 m above sea level. Through logistic regression we estimated the probability of P. perniciosus presence at each sampling site as a function of environmental and meteorological factors. Although in the initial univariate analyses the probability of P. perniciosus presence appeared to be associated with a wide variety of factors, in the multivariate logistic regression model only altitude, settlement, aspect, drainage hole construction, adjacent flora and the proximity of a sheep farm were retained as positive predictors of the distribution of this species.
P. perniciosus was present throughout the island, and thereby the risk of leishmaniosis transmission. The probability of finding P. perniciosus was higher at altitudes ranging from 51 to 150 m.a.s.l., with adjacent garrigue shrub vegetation, at the edge of or between settlements, and in proximity to a sheep farm.
The Balearic Islands in the Mediterranean region are considered endemic for both human and canine leishmaniosis, although the presence and prevalence of the diseases varies among the islands . The first data on human leishmaniosis in the Balearic Islands date from 1925 , while canine leishmaniosis was first reported in 1989 , in both cases in the island of Majorca, where most studies have been conducted.
In certain regions of Spain, human leishmaniosis is an endemic and notifiable disease, including in the Balearic Islands, which in some years have seen the highest registered incidence in Spain (4.72 and 4.59/100,000 in 2005 and 2006 respectively) . Between 7 and 33 cases are declared in Majorca every year [5, 6]. As in other parts of Spain, the disease is under-reported, especially cases of cutaneous leishmaniosis ; cases of human cryptic leishmaniosis have also been described . Little information is available on the origin of cases [6, 8].
The heterogeneous distribution and prevalence of canine leishmaniosis (CanL) ranges from 0% to 45% among different studies and islands [9–11]. A study conducted by the sanitary authorities in Majorca gave a prevalence of 14.4 % . Veterinarians answering a questionnaire on CanL trends in Majorca thought the disease was stable  and that autochthonous cases continue to occur, as has been previously described [3, 11].
Data on sand fly distribution in the Balearic Islands is scarce [10, 12–15]. The most recent data for Majorca corresponds to studies performed in 1987 and 1989 , but do not include information about the distribution and density of the different sand fly species throughout the island.
The aim of the present study was to obtain up-to-date entomological data by standardized methods that could be compared with data reported by other teams in different geographical areas of Europe and used in future entomological studies, including those on climate change. In addition, the extensive capture of the vector in the island could provide information on the environmental factors that may potentially favour the presence of P. perniciosus and also which areas are at risk of leishmaniosis.
Area of study
The study was carried out on the island of Majorca (Spain), located at 39°15’ to 39°60’N, 2°20’ to 3°30’E. Majorca is the largest of the Balearic Islands, covering 3,640 km2 and with a coastline of 623 km. Altitudes range from sea level to 1,445 m.a.s.l., most of the island (78.8%) being below 200 m.a.s.l. and only 6.3% above 500 m.a.s.l. The highest mountainous area is the Serra de Tramuntana in the north, which runs parallel to the west coast, protecting the island from the prevailing west and northwest winds. Bordering the low central plain in the southeast is the Serra de Llevant, with a maximum altitude of 509 m.a.s.l. .
The climate is typically Mediterranean, with long periods of invariability. The mean annual temperature is about 16–17°C, except in the Serra de Tramuntana, where it drops to 13°C. In the coldest period (1–3 months), the average temperature is about 5-10°C, with an average minimum on winter nights of 5–9°C, while in the hottest period (5–8 months) it is above 15-20°C, with an average diurnal maximum of 29-31°C. The mean relative humidity is 74%. Annual rainfall oscillates from a maximum in autumn (66.9 mm) to a minimum in summer (8.6 mm), with an annual average of 481.6 mm. Considerable differences exist between mountainous regions (up to 1,200 mm) and the arid south (less than 400 mm).
Holm oak (Cyclamini-Quercetum ilicis) grows everywhere on the island below 1000 m.a.s.l, but under the influence of man it has largely been replaced by pine (Pinus halepensis), which is now the dominant woodland tree, including all well-conserved beaches. In areas below 500–700 m.a.s.l., with annual precipitations of less than 600 mm, the wild olive tree predominates, while above 1000 m.a.s.l, the vegetation is low and adapted to strong winds. The extensive cultivated land consists principally of almond and olive trees, vineyards and cereals.
The island has two bioclimatic zones: meso-Mediterranean (T: average annual temperature 13–17°C; m: average minimum temperature of the coldest month -1 to -4°C; M: average maximum temperature of the coldest month 9 - 14°C; Ti: thermicity index 210–350), where oaks predominate (Cyclamini-Quercetum ilicis), and thermo-Mediterranean (T: 17–19°C; m: 4–10°C; M: 14 - 18°C; Ti: thermicity index from 350–470) with maquis (Cneoro-Ceratonietum)) [16, 17].
Capture of sand flies
In July 2008 sand fly captures were carried out in Majorca with sticky castor oil interception traps (20×20 cm) set for 4 days according to the standardized methodology implemented in the EDEN project (EU) [18–21]. The sampling sites consisted of holes used to drain embankments or containment walls, which were considered to be likely diurnal resting places for adult sand flies . The capture stations were distributed in 77 grids (5×5 km2), almost one station per grid, covering the entire island. A total of 1,882 sticky traps were set, representing an adhesion surface of 150.56 m2 distributed among 111 stations (Figure 1).
Data collection and environmental and meteorological variables
The characteristics of the stations, including location, habitat, environment and fauna, were recorded on a PDA (Palm Tungsten T5) using Pendragon Form v.5.0 software (PSC, Libertyville, IL, USA) and GPS (Tom Tom Wireless GPS MK II). Maps were designed using ArcGIS 9.2 software (ESRI, Redlands, CA, USA).
Climate variables were provided by the Spanish Meteorological Agency (AEMet) from the 43 meteorological stations in the study area. Different periods were considered for the meteorological variables: i) Period 1 (Sampling Day 1, when traps were set, to Day 4, when traps were recovered) and ii) Period 2 (the month before Sampling Day 1). Climate data from the nearest meteorological station were assigned to each sampling site for periods 1 and 2 using the spatial join-and-relate tool of ArcGis v.9.2 software and included: wind speed (Km/h), mean relative humidity (%), average rainfall (mm), and mean daily T (°C). The average minimum T (°C) in winter was also assigned.
Altitude data for each geocoded collection site were extracted from a 90 m resolution CGIAR Digital Elevation Model  using ArcGIS 9.2 software.
The presence of animals was studied in two ways: taking into account the animals or animal signs observed during captures, and using databases provided by the Col · legi Oficial de Veterinaris de les Illes Balears (canine census and livestock farms). In the latter, data from the closest station were assigned for each sampling site as described for the meteorological stations. A human census was obtained from municipal data.
Sand fly processing and identification
Sand flies were processed as previously described . Briefly, sand flies were removed from the sticky traps with a brush and fixed in 96% ethanol and then in 70% alcohol until identification. Males and Sergentomyia spp. females were observed and identified under the stereo microscope. Females of the genus Phlebotomus were mounted on glass slides in Hoyer medium and identified on the basis of morphological characteristics in an optical microscope using the keys of Gállego .
For the study 57 variables were taken into account, including the habitat and environmental characteristics of the capture stations, fauna, demography and climate.
Bivariate logistic regression studies were conducted using the SPSS 20.0 software for Windows, with all the independent variables set against the presence/absence of P. perniciosus as the dependent variable. The majority were used as categorical variables, except those related to meteorological conditions. Continuous variables such as the human and canine census were categorized in the search for association with the dependent variable. The possibility of interaction and/or confusion between different variables was examined by constructing and comparing different logistic regression models .
To construct the multivariate model, all the variables with p <0.2 in the bivariate study were used. In the final multivariate model, variables with p ≤0.05 were retained.
88.2 % of the traps placed on the island of Majorca were recovered, representing a surface of 135.68 m2. A total of 14,412 specimens were captured, with 4 species identified (Table 1): Phlebotomus pernicious, P. sergenti, P. papatasi and Sergentomyia minuta.
Among the mamophilic species, P. perniciosus was captured throughout the island in 77 of the 111 stations prospected, at 6 to 772 m.a.s.l. (Figure 1), with climate conditions during the capture period of 19.6-27.4°C, 55.5-86.4% relative humidity, 0–42 mm pluviometry and 3.1-17.1 km/h wind speed. P. sergenti and P. papatasi were captured in only 14 and 1 of the stations, respectively, and always in a low number. P. ariasi was not found anywhere on the island.
The bivariate analysis of the factors associated with the presence of P. perniciosus gave results of p <0.2 for 24 of the variables, which were taken into account in the multivariate analysis. 12 of these variables showed significant association (p <0.05) with the sand fly presence in both bivariate and multivariate analyses (Table 2).
The probability of capturing P. perniciosus was significantly higher at 51 – 150 m.a.s.l. (O.R. = 3.13), at the edge of or between settlements (O.R = 5.3), on a paved drive (O.R. = 2.90), in a wall drainage hole (not embankment) (O.R. = 2.11), in a general rural agricultural or forestry habitat (O.R. = 2.98), with adjacent flora of garrigue shrubs (O.R. = 14.53), in an agricultural area (O.R. = 5.52), and in the proximity of a sheep farm (O.R. = 2.72).
In contrast, the probability of capturing P. perniciosus showed a negative correlation with walls of bricks and mortar (O.R. = 0.26), non arable areas (O.R. = 0.27) and the proximity of pigeon and bee farms (O.R. = 0.15 and 0.1, respectively) (Table 2).
To construct the multivariate model, all the 24 variables with p <0.2 in the bivariate study were used. The variables that make up the multivariate logistic regression model and are shown to be the best predictors of the presence/absence of P. perniciosus are: an altitude of 51–150 m.a.s.l. (p = 0.01, O.R. = 8.6), location of the sampling sites at the edge of or between villages (p = 0.08, O.R. = 8.08), a south east orientation (p = 0.018, O.R. = 34.97), the absence of drainage holes with plastic pipes (p = 0.05, O.R. = 3.45), adjacent flora of garrigue shrubs (p = 0.001, O.R. = 38.05) and the proximity of a sheep farm (p = 0.001, O.R. = 20) (Table 3).
Four out of the five species previously reported for the island of Majorca (P. perniciosus, P. ariasi, P. sergenti, P. papatasi and S. minuta) [4, 12–14, 25] were captured. Although P. ariasi is cited [4, 13, 25], we were unable to capture this species despite sampling the whole island from 0 to 772 m.a.s.l and using a large number of traps. In Europe P. ariasi has been found at altitudes ranging from 10 m up to 2000 m.a.s.l. [20, 26], showing a preference for sub-humid or humid areas with cold winters (supra-Mediterranean) [21, 22, 27], while Majorca has a semi-arid and sub-humid climate with mild summers (meso- and thermo-Mediterranean). The repeated reporting of P. ariasi in Majorca may stem from an erroneous citing, which has been duplicated in other publications. Nevertheless, in this study, although captures were made throughout the whole island, they were restricted to the month of July (2008). Therefore, in order to assess more accurately whether P. ariasi is present or absent from the island, captures need to be made at different periods of sand fly activity. Also, intensive studies using CDC light traps should be carried out over 700 m a.s.l. in the mountainous regions of the island, particularly the area of the Serra de Tramuntana.
Among the species found, only P. perniciosus is a vector of L. infantum, and is responsible for human and canine leishmaniosis in the Mediterranean region [28, 29], while P. sergenti and P. papatasi are proven vectors of other Leishmania species in the Old World that are not present in Spain (L. tropica and L. major, respectively) [7, 30–32].
The most common sand fly species in Majorca is S. minuta, followed by P. perniciosus, P. sergenti and P. papatasi. The capturing method may have influenced the abundance level of each species, since it is known that sticky traps favor the capture of S. minuta females, which could be due to the feeding habits of this herpetophilic species and its preferred resting sites [24, 26]. Not enough P. sergenti and P. papatasi were captured for a statistical analysis of the factors affecting their presence in Majorca. As mentioned previously, most of the island is below 200 m.a.s.l., with a semi-arid climate, which are ideal conditions for P. sergenti to occur [33–35], yet this species was found at a low frequency (12.6 %). In other areas of Spain , P. sergenti has been found at altitudes of 0–1,153 m.a.s.l. and in the same type of meso- and thermo-Mediterranean bioclimates as in Majorca. Perhaps the location of traps within urbanized settlements (21 stations) or at the edge of/between settlements (90 stations), with little or no presence of humans, influenced the results, since P. sergenti is a peridomestic and anthropophilic species found in rural villages  and rare in intensely urban areas . The other scarcely sampled species, P. papatasi, prefers peri-arid and Saharan environments , not present in Majorca.
P. perniciosus was captured in Majorca from 6 to 772 m a.s.l., the maximum altitude at which the sticky traps were placed, since above that there was a lack of appropriate locations for setting traps. In Europe, the species occupies sites from sea level to 1534 m a.s.l. [19, 20, 26]. The probability of finding P. perniciosus was significantly higher at altitudes of 51 – 150 m.a.s.l., both in the bivariate and multivariate analysis. Stations at 0 – 50 m.a.s.l. were located in breezy coastal areas and sand flies are very sensitive to windy conditions [26, 29, 30]. In locations at 51 – 150 m.a.s.l. the adjacent flora consisted principally of garrigue shrubs, where the probability of finding P. perniciosus is significantly higher.
Locations between or at the edge of settlements favored the presence of P. perniciosus compared to those within settlements, as found in other studies [1, 18, 19, 21], which would indicate that urban environments are not suitable for P. perniciosus. The barbicans and other locations where sticky traps were placed constituted resting sites, which are often near the larval breeding sites [22, 26, 29]. In agreement with the site location, a positive correlation was obtained with a rural agricultural and forestry environment, where the probability of finding P. perniciosus was 3 times higher than in a rural village, as well as with an area of agricultural land cover, where the probability was more than 5 times higher than in urbanized areas. These results also match the negative correlation found in non-arable points of capture, usually in rural and/or urbanized areas, where the probability of capturing P. perniciosus decreased in comparison with stations near arable areas (cereals). In non-urbanized areas the terrestrial cycle of immature forms would be favored, and the females would have more access to suitable oviposition sites [18, 21]. In addition, the deployment of insecticides in urbanized areas during the summer period when blood-sucking insects are active would reduce the population of sand flies in those settlements, and it is considered a way of controlling leishmaniosis .
The presence of animals near the sampling site increased the probability of encountering P. perniciosus, for several reasons: i) the presence of animal excrements would constitute a good sand fly breeding substrate; ii) sand flies have a poor capacity for flying and dispersing far from their breeding sites (usually 300 m and rarely more than1 km) [26, 29, 30], which may explain the existence of small localized populations ; and iii) P. perniciosus exhibits opportunistic feeding behavior [39–42]. Nevertheless, in contrast with previous studies [1, 18, 19], no correlation was found with the presence of animals or animal traces such as feces near the trapping sites, only with an abundance of animals in livestock farms. Not all livestock species attract P. perniciosus in the same way , and its capture increased significantly when sheep farms were near to the sampling site. Notably, sheep farms contain a greater number of animals that remain outside overnight, when sand flies are active. No demographic influence of humans or dogs was found, probably because the stations with the highest presence of P. perniciosus were located between villages, away from urban settlements.
Some other variables correlated with the presence of P. perniciosus only in the bivariate analysis, such as the type of road, site category, land cover, wall construction and arable area, while the type of drainage hole correlated only in the multivariate analysis. The probability of capturing P. perniciosus in a paved drive was 2.9 times higher than in a paved public road, where greater car traffic would disturb sand flies. Drainage holes in non-embankment walls favored the presence of P. perniciosus in contrast with those in embankments, probably because the former have no air currents. On the contrary, the presence of P. perniciosus decreased by 75% in stone or brick walls with mortar, probably because these have fewer suitable resting places than walls without mortar. As described elsewhere, the use of PVC in drainage holes decreased the probability of finding P. perniciosus and could be considered as a control method to reduce leishmaniosis transmission .
The influence of climate variables on the distribution and activity of sand flies has been repeatedly reported [26, 30, 31, 43]. In contrast with other reports [18, 19, 21, 41], in the current study in Majorca, climate variables did not affect the probability of finding P. perniciosus, probably due to the short period of time when captures were performed (July 2008) and the homogenous geographical conditions of most trapping sites. It should also be taken into account that the island of Majorca has a Mediterranean climate, which remains highly stable over long periods, with the exception of the mountainous areas, and captures were not made over 700 m.a.s.l., due to the absence of appropriate places to set traps. More studies involving periodic captures throughout the summer, or over one year are required, as has been done in another Balearic island (Minorca) , to obtain more data on the influence of climate conditions on sand fly distribution.
The presence of P. perniciosus in Majorca is a health issue since it is a vector of L. infantum in the Mediterranean area. Leishmaniosis poses a risk not only for the habitual inhabitants of the island, but also for the large numbers of tourists visiting in the summer, coinciding with the period of vector activity. In addition, these tourists often travel with their pets, which are at risk of developing CanL. In central and northern European countries cases of leishmaniosis have repeatedly been reported in humans and dogs that have visited endemic areas [43–45]. Recent accounts of sand flies with a proven or suspected capacity to transmit L. infantum in non-endemic areas [46, 47], together with the arrival of infected persons and animals, would favor the possibility of autochthonous transmission in new areas, as has been reported in the island of Minorca .
The sand fly fauna in Majorca is composed of four species: P. perniciosus, P. sergenti, P. papatasi and S. minuta. The distribution of P. perniciosus extends throughout the island, from sea level to the mountains, being present in 70 % of the capture sites. This suggests that a risk of leishmaniosis transmission exists all over the island, and the presence of tourists during the period of P. perniciosus activity could favor the dispersion of the disease to other countries. The probability of finding P. perniciosus was higher at altitudes ranging from 51 to 150 m.a.s.l., with adjacent garrigue shrub vegetation, at the edge of or between settlements, and in proximity to a sheep farm.
Alcover MM, Ballart C, Serra T, Castells X, Scalone A, Castillejo S, Riera C, Tebar S, Gramiccia M, Portús M, Gállego M: Temporal trends in canine leishmaniosis in the Balearic Islands (Spain): A veterinary questionnaire. Prospective canine leishmaniosis survey and entomological studies conducted on the Island of Minorca, 20 years after first data were obtained. Acta Trop. 2013, 128: 642-651.
Pittaluga G: Etude épidémiologique sur la “Leishmaniose viscérale” en Espagne. Rapport présentée à l’Organisation d’Hygiène de la Société des Nations (8 octobre 1925). 1925, Genève: Societé des Nations 1700 (P) 10/25, 1700-
Matas-Mir B, Rovira-Alos J: Estudio epidemiológico de la leishmaniosis canina en la isla de Mallorca. 1989, Palma de Mallorca: Conselleria de Sanitat i Seguretat Social del Govern Balear
Amela C, Suarez B, Isidoro B, Sierra MJ, Santos S, Simón F: Evaluación del riesgo de transmisión de Leishmania infantum en España. 2012, Madrid: Centro de Coordinación de Alertas y Emergencias sanitarias (CCAES), Ministerio de Sanidad, Servicios Sociales e Igualdad
Xarxa de Vigilància Epidemiològica de les Illes Balears: Informe 2012. 2013, Palma de Mallorca: Conselleria de Salut, Direccio General de Salut Publica i Consum, Govern de les Illes Balears
Fulls setmanals de Vigilància Epidemiològica: Fulls setmanals de Vigilància Epidemiològica. 2002, Palma de Mallorca: Servei d’Epidemiologia, Direcció General de Salut Pública i Consum, Conselleria de Salut del Govern Balear, April 2014
Alvar J, Vélez ID, Bern C, Herrero M, Desjeux P, Cano J, Jannin J, den Boer M, the WHO Leishmaniasis Control Team: Leishmaniasis Worldwide and Global Estimates of Its Incidence. PLoS One. 2012, 7: e35671-
Riera C, Fisa R, López-Chejade P, Serra T, Girona E, Jiménez MT, Muncunill J, Sedeño M, Mascaró M, Udina M, Gállego M, Carrió J, Forteza A, Portús M: Asymptomatic infection by Leishmania infantum in blood donors from the Balearic Islands (Spain). Transfusion. 2008, 48: 1383-1389.
Alvar J: Las Leishmaniasis: de la Biología al Control. 2001, Salamanca: Laboratorios Intervet S.A
Seguí MG: Estudi epidemiològic de la leishmaniosis a l’illa de Menorca. Revista de Menorca. 1991, 2: 153-178.
Pujol A, Cortés E, Ranz A, Vela C, Aguiló C, Martí B: Estudi de seroprevalència de leishmaniosi i d’ehrlichiosi a l’illa de Mallorca. Revista del Col legi Oficial de Veterinaris de les Illes Balears Veterinària. 2007, 32: 9-12.
Gil Collado J: Phlébotomes et leishmanioses en Espagne. Écologie des Leishmanioses, Montpellier 1974. Colloques Internationaux du Centre National de la Recherche Scientifique N. 239. 1977, Paris: Centre National de la Recherche Scientifique, 177-189.
Gil-Collado J, Morillas-Márquez F, Sanchís-Marin MC: Los flebotomos en España. Rev San Hig Pub. 1989, 63: 15-34.
Lladó MT, Rotger MJ: Estudio del flebotomo como vector de la leishmaniasis en la isla de Mallorca. 1990, Palma de Mallorca: Conselleria de Sanitat i Seguretat Social del Govern Balear
Molina R, Aransay A, Nieto J, Cañanvate C, Chicharro C, Sans A, Flores M, Cruz I, García E, Cuadrado J, Alvar J: The Phlebotomine sand flies of Ibiza and Formentera Islands (Spain). Arch Inst Pasteur Tunis. 2005, 82: 12-13.
de Bolós i Capdevila O: La vegetació de les Illes Balears. Comunitat de plantes. 1997, Barcelona: Institut d’Estudis Catalans, 2
Rivas–Martínez S: Memoria del mapa de series de vegetación de España. 1987, Ministerio de Agricultura, Pesca y Alimentación. I.C.O.N.A
Gálvez R, Descalzo MA, Miró G, Jiménez MI, Martín O, Dos Santos-Brandao F, Guerrero I, Cubero E, Molina R: Seasonal trends and spatial relations between environmental/meteorological factors and leishmaniosis sand fly vector abundances in Central Spain. Acta Trop. 2010, 115: 95-102.
Barón SD, Morillas-Márquez F, Morales-Yuste M, Díaz-Sáez V, Irigaray C, Martín-Sánchez J: Risk maps for the presence and absence of Phlebotomus perniciosus in an endemic area of leishmaniasis in southern Spain: implications for the control of the disease. Parasitology. 2011, 138: 1234-1244.
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-159.
Ballart C, Guerrero I, Castells X, Barón S, Castillejo S, Alcover MM, Portús M, Gállego M: Importance of individual analysis of environmental and climatic factors affecting the density of Leishmania vectors living in the same geographical area: the example of Phlebotomus ariasi and P. perniciosus in Northeast Spain. Geospat Health. 2014, 8: 367-381.
Rioux JA, Carron S, Dereure J, Périères J, Zeraia L, Franquet E, Babinot M, Gállego M, Prudhomme J: Ecology of leishmaniasis in the South of France. 22. Reliability and representativeness of 12 Phlebotomus ariasi, P. perniciosus and Sergentomyia minuta (Diptera: Psychodidae) sampling stations in Vallespir (eastern French Pyrenees region). Parasite. 2013, 20: 34-
Jarvis A, Reuter HI, Nelson A, Guevara E: Hole – filled seamless SRTM data V4. 2008, International Centre for Tropical Agriculture (CIAT),http://srtm.csi.cgiar.org/,
Gállego J, Botet J, Gállego M, Portús M: Los flebotomos de la España peninsular e Islas Baleares. Identificacióny corología. Comentarios sobre los métodos de captura. Memoriam al Prof. Dr. D. F. de P. Martínez Gómez. Edited by: Hernández S. 1992, Córdoba: Publicaciones de la Universidad de Córdoba, 581-600.
Lucientes J, Castillo JA, Gracia MJ, Peribáñez MA: Flebotomos, de la biología al control. Revista Electrónica de Veterinaria REDVET. 2005, 6: 1-8.
Rioux JA, Golvan YJ: Épidémiologie des Leishmanioses dans le Sud de la France. 1969, Paris: Monographies de l’Institut National de la Santé et de la Recherche Médicale. N° 37
Aransay AM, Testa JM, Morillas-Marquez F, Lucientes J, Ready PD: Distribution of sandfly species in relation to canine leishmaniasis from the Ebro Valley to Valencia, northeastern Spain. Parasitol Res. 2004, 94: 416-420.
Gállego M: Zoonosis emergentes por patógenos parásitos: las leishmaniosis. Rev Sci Tech. 2004, 23: 661-676.
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-147.
Killick-Kendrick R: The biology and control of Phlebotomine sand flies. Clin Dermatol. 1999, 17: 279-289.
Ready PD: Leishmaniasis emergence in Europe. Euro Surveill. 2010, 15: 19505-
WHO: Control of the leishmaniases.: Report of a WHO expert committee. WHO Tech Rep Ser. 2010, 949: 1-186. [http://whqlibdoc.who.int/trs/WHO_TRS_949_eng.pdf]
Rioux JA, Rispail P, Lanotte G, Lepart J: Relation Phlébotomes-bioclimats en écologie des leishmanioses. Corollaires épidémiologiques. L’exemple du Maroc. Bulletin de la Société botanique de France. 1984, 131: 549-557.
Gállego M, Rioux JA, Rispail P, Guilvard E, Gállego J, Portús M, Delalbre A, Bastien P, Martínez-Ortega E, Fisa R: Primera denuncia de flebotomos (Diptera, Psychodidae, Phlebotominae) en la provincia de Lérida (España, Cataluña). Rev Iber Parasitol. 1990, 50: 123-127.
Barón SD, Morillas-Márquez F, Morales-Yuste M, Díaz-Sáez V, Gállego M, Molina R, Martín-Sánchez J: Predicting the risk of an endemic focus of Leishmania tropica becoming established in south-western Europe through the presence of its main vector, Phlebotomus sergenti Parrot, 1917. Parasitology. 2013, 140: 1413-1421.
Boussaa S, Neffa M, Pesson B, Boumezzough A: Phlebotomine sandflies (Diptera: Psychodidae) of southern Morocco: results of entomological surveys along the Marrakech- Ouarzazat and Marrakech-Azilal roads. Ann Trop Med Parasitol. 2010, 104: 163-170.
Alexander B, Maroli M: Control of phlebotomine sandflies. Med Vet Entomol. 2003, 17: 1-18.
Belen A, Alten B, Aytekin AM: Altitudinal variation in morphometric and molecular characteristics of Phlebotomus papatasi populations. Med Vet Entomol. 2004, 18: 343-350.
de Colmenares M, Portús M, Botet J, Dobaño C, Gállego M, Wolff M, Seguí G: Identification of blood meals of Phlebotomus perniciosus (Diptera: Psychodidae) in Spain by a competitive enzyme-linked immunosorbent assay biotin/avidin method. J Med Entomol. 1995, 32: 229-233.
Rossi E, Bongiorno G, Ciolli E, Di Muccio T, Scalone A, Gramiccia M, Gradoni L, Maroli M: Seasonal phenology, host-blood feeding preferences and natural Leishmania infection of Phlebotomus perniciosus (Diptera, Psychodidae) in a high-endemic focus of canine leishmaniasis in Rome province, Italy. Acta Trop. 2008, 105: 158-165.
Branco S, Alves-Pires C, Maia C, Cortes S, Cristovão JMS, Gonçalves L, Campino L, Afonso MO: Entomological and ecological studies in a new potential zoonotic leishmaniasis focus in Torres Novas municipality, Central Region, Portugal. Acta Trop. 2013, 125: 339-348.
Jiménez M, González E, Martín-Martín I, Hernández S, Molina R: Could wild rabbits (Oryctolagus cuniculus) be reservoirs for Leishmania infantum in the focus of Madrid, Spain?. Vet Parasitol. 2014, 202: 296-300.
Aspöck H, Gerersdorfer T, Formayer H, Walochnik J: Sand flies and sandfly-borne infections of humans in Central Europe in the light of climate change. Wien Klin Wochenschr. 2008, 120: 24-29.
Poeppl W, Herkner H, Tobudic S, Faas A, Auer H, Mooseder G, Burgmann H, Walochnik J: Seroprevalence and asymptomatic carriage of Leishmania spp. in Austria, a non-endemic European country. Clin Microbiol Infect. 2013, 19: 572-577.
Rioux JA, Houin H, Ranque J, Lapierre J: Écologie des Leishmanioses. Colloques Internationaux du Centre National de la Recherche Scientifique, 1974. 1977, Paris: Centre National de la Recherche Scientifique, 190-N. 239
Haeberlein S, Fischer D, Thomas SM, Schleicher U, Beierkuhnlein C, Bogdan C: First Assessment for the Presence of Phlebotomine Vectors in Bavaria, Southern Germany, by Combined Distribution Modeling and Field Surveys. PLoS One. 2013, 8: e81088-
Poeppl W, Obwaller AG, Weiler M, Burgmann H, Mooseder G, Lorentz S, Rauchenwald F, Aspöck H, Walochnik J, Naucke TJ: Emergence of sandflies (Phlebotominae) in Austria, a Central European country. Parasitol Res. 2013, 112: 4231-4237.
This work was supported by grants of the Ministerio de Ciencia y Tecnología of Spain (CGL2007-66943-C02-01/BOS), Departament d’Universitats, Recerca i Societat de la Informació de la Generalitat de Catalunya (Spain) (2009SGR385). The Spanish Meteorological Agency (AEMet) provided the meteorological data for the study. Thanks to Anna Lanau for her assistance in placing traps and collecting sand flies. We are also grateful for the help of the Col · legi Oficial de Veterinaris de les Illes Balears, especially R. García and A. Figueroa. MMA was awarded a contract in the Spanish project.
The authors declare that they have no competing interests. The contents are the sole responsibility of the authors.
MGC, MPV, MMAA designed and supervised the study. MGC, CBF, TSF, SCG, MMAA undertook field and laboratory activities. MGC, JMS, MMAA analyzed the data and carried out the statistical analysis, MGC, MPV, JMS, MMAA drafted and revised the manuscript. All the authors revised and approved the final version of the manuscript.
Authors’ original submitted files for images
Below are the links to the authors’ original submitted files for images.