|Home | About | Journals | Submit | Contact Us | Français|
An unusual increase of human leishmaniasis cases due to Leishmania infantum is occurring in an urban area of southwestern Madrid, Spain, since 2010. Entomological surveys have shown that Phlebotomus perniciosus is the only potential vector. Direct xenodiagnosis in hares (Lepus granatensis) and rabbits (Oryctolagus cuniculus) collected in the focus area proved that they can transmit parasites to colonized P. perniciosus. Isolates were characterized as L. infantum. The aim of the present work was to conduct a comprehensive study of sand flies in the outbreak area, with special emphasis on P. perniciosus.
Entomological surveys were done from June to October 2012–2014 in 4 stations located close to the affected area. Twenty sticky traps (ST) and two CDC light traps (LT) were monthly placed during two consecutive days in every station. LT were replaced every morning. Sand fly infection rates were determined by dissecting females collected with LT. Molecular procedures applied to study blood meal preferences and to detect L. infantum were performed for a better understanding of the epidemiology of the outbreak.
A total of 45,127 specimens belonging to 4 sand fly species were collected: P. perniciosus (75.34%), Sergentomyia minuta (24.65%), Phlebotomus sergenti (0.005%) and Phlebotomus papatasi (0.005%). No Phlebotomus ariasi were captured. From 3203 P. perniciosus female dissected, 117 were infected with flagellates (3.7%). Furthermore, 13.31% and 7.78% of blood-fed and unfed female sand flies, respectively, were found infected with L. infantum by PCR. The highest rates of infected P. perniciosus were detected at the end of the transmission periods. Regarding to blood meal preferences, hares and rabbits were preferred, although human, cat and dog blood were also found.
This entomological study highlights the exceptional nature of the Leishmania outbreak occurring in southwestern Madrid, Spain. It is confirmed that P. perniciosus is the only vector in the affected area, with high densities and infection rates. Rabbits and hares were the main blood meal sources of this species. These results reinforce the need for an extensive and permanent surveillance in this region, and others of similar characteristics, in order to control the vector and regulate the populations of wild reservoirs.
The online version of this article (doi:10.1186/s13071-017-2309-z) contains supplementary material, which is available to authorized users.
Leishmaniasis is caused in Spain by the trypanosomatid Leishmania infantum and visceral, cutaneous and mucosal forms of the disease are notified in the country, although these last two are underestimated because they are usually unnoticed [1–3]. Phlebotomine sand flies are the cornerstone in the transmission of leishmaniasis in many geographical regions of the world including the countries of the Mediterranean basin. The deep knowledge of wild cycles of Leishmania is still scarce despite being a very important aspect in the epidemiology of the disease. The proven vectors implicated in the transmission are Phlebotomus perniciosus and Phlebotomus ariasi . Most of cases of leishmaniasis are reported in the Mediterranean coast and in central region of Spain . During the period 2000–2009, 12–25 annual cases were reported in the region of Madrid (central Spain), but since 2010 an unusual increase of both visceral and cutaneous leishmaniasis cases was observed in southwestern Madrid region, mainly in four urban areas: Fuenlabrada, Leganés, Getafe and Humanes de Madrid [5, 6]. A total of 691 human cases were reported in this area between 2010 and October 2016, Fuenlabrada being the most affected town by far, where the mean incidence reached 45.17 cases/100,000 inhabitants (data provided by the Community of Madrid).
The periurban area close to the focus was traditionally used for agricultural purposes with a significant presence of hares (Lepus granatensis) . Recently, it was modified in order to create a large green park for the enjoyment of the population living in the area. This change in land use has led to an increase of both hare and wild rabbit (Oryctolagus cuniculus) populations in this open space and has probably increased the population densities of sand fly vectors. Studies carried out in this leishmaniasis focus using direct xenodiagnosis have involved hares and rabbits from the area as wild reservoirs of the disease since they were able to transmit the parasite to colonized P. perniciosus sand flies [8, 9]. Even more, it has been demonstrated a high exposure of these lagomorphs to P. perniciosus bites . On the other hand, high levels of leishmaniasis seroprevalence have been detected in hares and rabbits from this area . These findings suggest that a sylvatic cycle of transmission of Leishmania exists in this periurban park independent of the classical urban domestic cycle with dog as the main reservoir . Additional studies carried out in Madrid and in other regions of Spain strengthen the implication of lagomorphs in the sylvatic cycle of L. infantum in the country [12–14].
On the other hand, several factors associated to human activity and climate change are influencing Leishmania distribution in some parts of the world . Studies have reported changes in sand fly distribution, as well as risk for their establishment in new areas of Europe . Specifically, in the Iberian Peninsula epidemiological and entomological studies have been performed in order to update and improve the knowledge of the eco-epidemiology of leishmaniasis and the distribution of its vectors [17–20]. In most of these studies entomological surveys when combined with molecular procedures demonstrate to be worthwhile in the understanding of the epidemiology of leishmaniasis.
This work aimed to carry out a detailed entomological survey in order to obtain information about sand fly seasonal trends and densities, the evolution of sand fly infection rates, and their blood meal preferences in the exceptional human focus of leishmaniasis that is affecting urban areas of the southwestern Madrid region (Spain).
The entomological surveys were conducted in a periurban green park of around 450 ha surrounded by the towns of Fuenlabrada, Leganés and Getafe (Madrid, Spain), the area affected by the leishmaniasis outbreak (Fig. (Fig.1).1). Three stations located in the bordering zone of the park with the urban area of the town of Fuenlabrada were chosen for the study, named ATE (40.292849N, 3.780539W), BOS (40.298084N, 3.793136W) and JIC (40.299954N, 3.806019W). A fourth station located inside the green park belonging to the municipality of Leganés was also selected, named POL (40.324903N, 3.796381W) (Fig. (Fig.1).1). These four stations were selected according to the high number of sand flies collected in a preliminary survey performed in 2011. The altitude of these stations varies from 655 to 691 m. The studied area belongs to the meso-Mediterranean bioclimatic zone, with an annual average temperature of 15 °C and annual precipitation of 365 mm (data from the Spanish Meteorological Agency). The characteristics of the soil for traditional rainfed cultivation make the area an exceptional habitat for hare and rabbit. The stations JIC, BOS and ATE are located in the contact line between the urban area and the green park where rabbits are abundant. In contrast, POL station is located in the middle of the green park where hares are the predominant leporids. Occasional colonies of stray cats can also be seen.
Sand flies were collected every month for three successive years (2012, 2013 and 2014) in the active season of sand flies, from June to October. Twenty sticky traps (ST) (20 × 20 cm paper sheets soaked in castor oil) and two CDC miniature light traps (LT) were used during two consecutive nights in each station. LT were replaced every morning (Fig. (Fig.1).1). Temperature and relative humidity (RH) were registered every 10 min using data loggers (Velleman® DVM171HD, Gavere, Belgium) hanged on each LT, with the exception of June 2012 in which the data were obtained from the nearest meteorological station located in Getafe. For phenology studies mean temperature and RH registered from dusk to dawn were used, according to information provided by the National Astronomical Observatory (Madrid).
Sand flies collected by ST were detached from the paper sheets using a fine brush, placed in 96% ethanol (to remove the castor oil) and then stored in 70% ethanol at 4 °C until processing. Female sand flies collected with LT were immediately dissected. Non-dissected females as a result of the high number of catches obtained in some cases and all males were stored in 70% ethanol at 4 °C. Blood-fed sand flies from both LT and ST were separated and stored in 70% ethanol at 4 °C for subsequent molecular studies.
Taxonomical identification of the collected sand flies was based on the morphology of the male genitalia and spermathecae of females according to Gil-Collado et al. . Female genitalia was cleared in Marc-André medium and transferred together with the head to a drop of Hoyer’s mounting medium on a glass slide for identification under a microscope.
Alive female sand flies collected with LT were anesthetized with CO2 and placed in a sterile Petri dish containing 5% detergent solution in sterile PBS. They were then individually transferred to a drop of sterile PBS placed on a sterilized glass slide. Sand flies were dissected by using flame-disinfected entomological needles and examined for the presence of promastigotes under a phase-contrast microscope . Promastigotes from guts of infected sand flies were transferred to screw crap tubes containing 200 μl of M199 medium supplemented with 20% fetal calf serum (FCS) (Hyclone™, GE Healthcare Life Science, Logan, Utah, USA) and 1.6% penicillin-streptomycin (10,000 U/ml of each antibiotic), pH 7.4 (Lonza BioWitthaker®, Verviers, Belgium). After 3–4 days at 27 °C positive vials containing live promastigotes were transferred to Novy-MacNeal-Nicolle (NNN) medium with RPMI medium for mass growth of promastigotes and further characterization of the isolates. Finally, after 1–4 weeks the promastigote cultures were stored in liquid nitrogen, depending on the evolution of their growth.
Sand flies preserved in 70% ethanol were washed individually in sterile distilled water placed in ELISA microtiter plates. Afterwards, head, wings, genitalia, and legs of each female sand fly were removed. The genitalia and the head were processed for taxonomical identification as previously described. The presence of a blood meal was determined by observation under a stereomicroscope. The phase of blood digestion was determined according to Dolmatova & Demina . Thorax and abdomen of each sand fly were used for DNA extraction using the DNeasy® Blood & Tissue Extraction Kit (Qiagen, Hilden, Germany) according to the manufacturer’s instructions. DNA from promastigote cultures was obtained by using the same kit. In both cases, two DNA elutions of 100 μl were obtained and further quantification and purity were determined by spectrophotometry with a NanoDrop ND-1000 spectrophotometer (Nanodrop Technologies, Wilmington, DE, USA). Finally, samples were stored at -20 °C until use.
Molecular characterization of the isolates was carried out by amplification of internal transcribed spacer regions 1 (ITS1) and 2 (ITS2) using the two pairs of primers as described in previous studies [9, 24]: (i) LITSR (5′-CTG GAT CAT TTT CCG ATG-3′) and L5.8S (5′-TGA TAC CAC TTA TCG CAC TT-3′) and (ii) L5.8SR (5′-AAG TGC GAT AAG TGG TA-3′) and LVTSV (5′-ACA CTC AGG TCT GTA AAC-3′). PCR products were separated on 1.5% agarose gel (Conda, Torrejón de Ardoz, Spain) stained with “Pronasafe Nucleic Acid Staining Solution” (10 mg/ml) (Conda, Torrejón de Ardoz, Spain) and visualized under UV light. Bands obtained were removed from the gel under UV exposure and purified with the QIAquick® Gel Extraction Kit (Qiagen). Afterwards, the DNA samples were sequenced with ABI PRISM 3730XL DNA Analyzer (Applied Biosystems, Foster City, CA, USA). Electropherograms were manually inspected and corrected using ChromasPro program (McCarthy, Queensland, Australia). Nucleotide sequences obtained were analysed with DNASTAR program (Lasergen v7.1®, Madison, WI, USA). Homologies with the available sequences data in GenBank was carried out with the software BLAST (http://www.ncbi.nlm.nih.gov/BLAST).
Blood meal identification of blood-fed sand flies was conducted by the amplification of a fragment of 359 bp of vertebrate cytochrome b (cyt b) gene followed by sequencing as described before [9, 24]. Universal primers cyto 1 (5′-CCA TCA AAC ATC TCA GCA TGA AA-3′) and cyto 2 (5′-CCC CTC AGA ATG ATA TTT GTC CTC-3′)  were used. Degenerated primers cyt_bb1 (5′-CCA TCM AAC ATY TCA DCA TGA TGA AA-3′) and cyt_bb2 (5′-GCH CCT CAG AAT GAY ATT TGK CCT CA-3′) were used in samples unsuccessfully amplified with non-degenerated primers . In the case of Sergentomyia minuta blood meal identification, degenerated primer for reptiles (5′-GCH GAY ACH WVH HYH GCH TTY TCH TC-3′) combined with primer cyto 2 was used . Band purification and subsequent sequencing were carried out as previously described.
Detection of Leishmania was carried out by amplification of a 120 bp fragment from kDNA and further confirmation by cysteine proteinase b (cpb)-PCR following previous protocols . DNA obtained from sand flies reared in the laboratory was used as negative control. DNA of L. infantum was used as a positive control. To prevent PCR contamination, sample preparation, reactions set-up, and PCR amplifications were performed in separate rooms, using different lab coats and gloves.
Relative abundance (RA) of each species was assessed by using LT (sand flies captured of one specific species/total sand flies captured by LT), while density was calculated by using sand flies collected by ST (number of sand flies captured per m2 of recovered ST). Sex ratio was calculated as the number of males/females × 100.
Statistical analysis with Kruskal-Wallis test was performed in order to study significant differences in sand fly captures, density, RA, temperature and RH between the three survey periods. For the analysis of bioclimatic parameters in June, Mann-Whitney test was performed in order to compare 2013 and 2014 data. Dunn’s multiple comparison tests was used to analyse differences between mean values of the mentioned parameters between each period and month. The possible relation between sand fly captures during the three periods and bioclimatic parameters (mean/ maximum/minimum temperature and RH) was determined using the Spearman’s correlation test. Analyses were only performed with P. perniciosus and S. minuta data.
A total of 45,127 specimens were collected during the three survey periods: 2012 (n = 17,317); 2013 (n = 12,407); and 2014 (n = 15,403). A total of 16,502 sand flies were captured by LT and 28,625 were captured by ST. The two predominant species captured were P. perniciosus and S. minuta. There were no statistically significant differences between LT captures among the three periods (Kruskal-Wallis H-test: χ 2 = 5.9915, df = 2, P = 0.3833 for P. perniciosus, P = 0.7591 for S. minuta). No significant differences were found between P. perniciosus and S. minuta captured by ST among the three surveys, P = 0.4441, P = 0.5824, respectively. Concerning the four stations studied, the highest number of sand flies was found in JIC (n = 17,398), followed by ATE (n = 12,768), BOS (n = 9158) and POL (n = 5803). The same correlation was separately observed in the collections obtained with both LT and ST, although numbers and species of sand fly captured by each method was very variable depending on the station as shown in Fig. Fig.2.2. Statistical analysis through the Kruskal-Wallis test showed no significant differences in P. perniciosus collections between the four stations with both capture traps (Kruskal-Wallis H-test: χ 2 = 7.8147, df = 3, P = 0.333, P = 0.0877 for LT and ST, respectively). However, S. minuta captures showed significant differences (Kruskal-Wallis H-test, χ 2 = 7.8147, df = 3; P = 0.0328 and P = 0.0003 for LT and ST, respectively). Dunn’s multiple comparison test showed significant differences in ST captures between stations BOS and JIC (P = 0.0195).
Four sand fly species were identified: P. perniciosus (n = 33,996; 75.34%), Phlebotomus sergenti (n = 2; 0.005%), Phlebotomus papatasi (n = 2; 0.005%) and S. minuta (n = 11,127; 24.65%). No P. ariasi was found. Phlebotomus sergenti (n = 2) were captured by LT in 2012 (1 in August-ATE and 1 in July-JIC). Phlebotomus papatasi (n = 2) were collected by ST in 2013 (July-BOS) and in 2014 (June-ATE). Mean density of P. perniciosus captured with ST in the three surveys was 193.6 specimens/m2. During the three surveys mean densities were 190.1, 175.72 and 220.45 specimens/m2 in 2012, 2013 and 2014, respectively. In the case of S. minuta the mean density was 106.81 specimens/m2, with variation between the surveys (87.06 in 2012; 97.34 in 2013 and 135.96 in 2014). Mean density of P. papatasi was 0.021 specimens/m2. Density values of P. perniciosus and S. minuta did not show significant variation between the three surveys (Kruskal-Wallis H-test: χ 2 = 5.9915, df = 2, P = 0.4163, P = 0.5824, respectively).
Regarding to RA calculated from LT captures, the mean value for P. perniciosus was 94.24%. Variation through the three periods was 94.74%, 91.85% and 95.37% in 2012, 2013, and 2014, respectively. This variation was not significant (Kruskal-Wallis H-test: χ 2 = 5.9915, df = 2, P = 0.0853). With regard to S. minuta, RA mean value was 5.75%. RA fluctuation through the three periods was 5.24% in 2012, 8.15% in 2013 and 4.63% in 2014, showing no significant deviation (Kruskal-Wallis H-test: χ 2 = 5.9915, df = 2, P = 0.4441). The mean RA value for P. sergenti was 0.012% (Table (Table11).
Captures with LT showed a clear predominance of P. perniciosus. This species showed a peak in August 2012 and in September in 2013 and 2014 (Fig. (Fig.3).3). Although only significant variation was found in monthly captures of S. minuta during the three periods (Kruskal-Wallis H-test: χ 2 = 5.9915, df = 2; P = 0.0766, P = 0.0025, for P. perniciosus and S. minuta, respectively). Specifically, multiple comparison tests showed significant variation in S. minuta captures between August and October (P = 0.0467). The difference between the two predominant species was less marked with ST (Fig. (Fig.4).4). In any case, the number of P. perniciosus captured was higher than S. minuta, with the exception of August 2013 and 2014. Phlebotomus perniciosus density showed two peaks in June and August 2012, while in 2013 a weak peak in July and a higher one in September were observed. In 2014, P. perniciosus density only showed one peak (September). On the other hand, S. minuta showed a peak in August during the three periods studied; however, in September 2014 the captures also remained high (Fig. (Fig.4).4). Significant differences in P. perniciosus monthly densities between each period were found (Additional file 1: Table S1).
Phlebotomus perniciosus male rate exceeded 50% with both LT and ST during all surveys, except with LT in 2012 whereas S. minuta male rate also appeared more than 50% in ST during the three surveys periods (Table (Table11).
Relative humidity and temperature of this study corresponds to values recorded between dusk and dawn of the two nights of each survey. The different bioclimatic parameters studied in this work presented fluctuations between the three survey periods. On one hand, mean RH presented significant variation in June (Kruskal-Wallis H-test: χ 2 = 5.9915, df = 2; P = 0.0286). On the other hand, Dunn’s test of maximum RH showed significant differences between July 2012, 2013 (P = 0.0003); August 2012, 2013 (P = 0.0442) and August 2012, 2014 (P = 0.0001). Regarding minimum RH, Dunn’s test showed significant differences between August 2012, 2013 (P = 0.0217), August 2013, 2014 (P = 0.0001), October 2012, 2013 (P = 0.0195) and October 2012, 2014 (P = 0.0001) (Additional file 2: Table S2).
In case of temperature variation, no significant differences were found in monthly mean temperatures recorded over the 3 years (Additional file 3: Table S3). However, maximum temperatures showed significant fluctuations during the three collection periods except in September (Kruskal-Wallis H-test: χ 2 = 5.9915, df = 2; June, P = 0.0002; July, P = 0.0468; August, P = 0.0007; September, P = 0.0621; October, P = 0.0003). Regarding minimum temperatures significant differences were observed during the three sampling periods (Kruskal-Wallis H-test: χ 2 = 5.9915, df = 2; June, P = 0.0002; July, P = 0.0358; August, P = 0.0031; September, P = 0.0063; October, P = 0.0001) (Additional file 3: Table S3).
Statistical analysis showed a negative correlation between sand flies captured during the surveys by LT and ST and mean, maximum and minimum RH (Spearman’s rho = -0.4 – -1, n = 12, P = 0.5167–0.0167) but the only significant negative correlation was only observed between S. minuta captured by LT and minimum RH (Spearman’s rho = -1, P = 0.0167) (Additional file 4: Table S4). Regarding temperature, sand fly captures positively correlated with mean temperature, maximum temperature and minimum temperature in both LT and ST (Spearman’s rho = 0.4–1, n = 12, P = 0.5167–0.0167) with a significant correlation only observed between S. minuta captured by LT and maximum temperature values (Spearman’s rho = 1, P = 0.0167) (Additional file 4: Table S4).
During the three survey periods, a total of 3203 P. perniciosus females collected with LT in the 4 stations of seasonal study were dissected. The global infection rate was 3.65%. In 2012, 19 out of 735 females dissected were found infected with promastigotes (2.6%). In 2013, 864 females were dissected with 57 of them infected, giving a rate of infection of 6.6%. In 2014, 1604 females were dissected and promastigotes were found in 41 specimens (2.6%). Of the dissected infected females 7.7% were gravid and 6% were semigravid. As shown in Fig. Fig.5,5, infection rates dynamics fluctuated during the three periods, although 2013 was the period with highest infection rates. Specifically, an exceptional rate of infected sand flies was found by dissection in August and October 2013. Overall, sand fly dissections revealed higher parasite presence in August, September and October during the three periods. Of the isolates from dissected guts, 19.2% resulted contaminated by fungi or bacteria and in the 5.1% of the cultures the growth of the parasites was unsuccessful due to the low number of promastigotes recovered from the infected guts. A total of 67 isolates were successfully cultured and this allowed proceeding to molecular characterization. ITS1 and ITS2 sequences were obtained and analysed by comparison with GenBank database. All the isolates showed an identity of 100% with L. infantum strain MHOM/ES/87/Lombardi (ITS type Lombardi; AJ000295) (Additional file 5: Figure S1).
A total of 912 P. perniciosus females were analysed during the study, 308 of them blood-fed. Host-feeding preferences of P. perniciosus were investigated over the 3 years of the study, 2012 (n = 100), 2013 (n = 92) and 2014 (n = 116). Blood identification was achieved in 224 blood-fed P. perniciosus females (efficiency of 74.36%).The rate of unidentified blood meal sources in P. perniciosus engorged females was similar along the 3 years of the study: 29%, 23.07%, and 25% in 2012, 2013 and 2014, respectively. The following hosts have been identified: rabbit (50.33%), hare (19.16%), cat (3.57%), human (0.97%) and dog (0.33%). Differences in the identified hosts were found between the three survey periods. The percentages of females fed on rabbits were increasing during the study: 41% in 2012, 50.54% in 2013, and 58.62% in 2014. In contrast, the rate of hares as blood meal source was decreasing, 28%, 17.58% and 12.93% in 2012, 2013 and 2014, respectively. Blood preferences by station are shown in Fig. Fig.6.6. It is remarkable that most of hare blood was detected in sand flies collected in the station LEG-POL and the station with more variable blood meal sources was FUE-SJIC.
Additionally the blood contained in the gut of 30.3% (10 out of 33) S. minuta female could be identified. All these blood meals were identified as gecko blood.
Regarding to molecular detection of Leishmania, a total of 41 positive blood-fed P. perniciosus females were found (13.31%). Three positive females were gravid (7.3%) and two were semigravid (4.87%). 2012 was the period with highest rate of infection (23%), in 2013 the rate was 11.11% and 4.76% in 2014. Sixteen (69.56%) out of 23 positive female sand flies found in 2012 were detected in POL station; 11 out of this 16 blood meal sources were identified as hare. Another two sand flies fed on rabbits were found positive in BOS station, and other 5 sand flies with unidentified blood meal source were found positive (2 in ATE and 3 in JIC stations). Ten out of 12 positive blood-fed females collected in 2013 were fed on rabbits (9 in BOS and 1 in ATE), in the other 2 flies the blood could not be identified (1 in BOS and 1 in POL). Finally, 3 out of 6 positive sand flies collected in 2014 were fed on rabbits and the other 3 on hares: 1 fed on hare in ATE, 2 fed on hares in POL and 3 fed on rabbits in BOS (Fig. (Fig.6).6). Infection rates of females fed on rabbits during the three collection periods were 4.97%, 21.74% and 4.41% in 2012, 2013 and 2014, respectively. As to the females fed on hares the infection rates were 39.28% in 2012, 0% in 2013 and 20% in 2014.
Concerning unfed females, 604 P. perniciosus were analysed for Leishmania detection: 225 collected in 2012; 192 in 2013 and 187 in 2014. A total of 47 females were positive for Leishmania (7.78%). The highest infection rate was found in 2013 (10.88%), followed by 2012 (6.69%) and 2014 (5.88%). Only one positive gravid female (2.12%) was observed. During the three survey periods, BOS was the station with higher number of infected sand flies, 21 (32.30%), followed by POL, 13 (11.71%), JIC, 10 (4.13%) and ATE, 3 (2.36%). Additionally, 40 S. minuta females were also investigated but no trypanosomatids were detected (Fig. (Fig.66).
In this work we present results on sand fly seasonal dynamics, as well as findings about sand fly infection rates and blood meal preferences in the exceptional outbreak of human leishmaniasis in an urban area of the southwestern region of Madrid (central Spain).
Anthropogenic land changes as deforestation-reforestation, road construction and urbanization seem to be a very important driver of infectious disease outbreaks giving raise to emergence or re-emergence events. These changes perturb host-parasite dynamics equilibrium in parasitic vectorial diseases, including leishmaniasis [15, 28]. The emergence of leishmaniasis in many areas can be associated with urbanization, ecological factors and climatic changes that favor the increase in vector densities and the modification of host population or its composition [29, 30]. In the case of the studied area the high densities of P. perniciosus found in this work could be explained because of the alteration of the periurban area which was modified in order to establish an extensive green area surrounding the different affected towns. This alteration was expected to modify hare and rabbit populations in this open space and as result, the sand fly density has probably increased . This finding agrees with observations on the plasticity and adaptability to environmental changes of P. perniciosus observed in Italy by Tarallo . Moreover, environmental changes are proved to affect sand fly population and density [32–35].
During the three survey periods four different sand fly species were found: P. perniciosus, P. sergenti, P. papatasi and S. minuta. Although P. ariasi have been recorded in some areas of Madrid region no specimen of this species was found in our surveys. Moreover, Tello et al.  only found one specimen during their collections in different sites of the affected area. We could not find any specific reason for the absence of this species, so it could be due to the homogeneity of bioclimatic and ecological variables of the study area that are unsuitable for the presence of this sand fly species [37, 38]. Phlebotomus perniciosus was the predominant species and S. minuta was the second in abundance, contrary to what has been reported elsewhere in the Madrid region [37, 38], probably because bioclimatic and environmental conditions of the area of study are more favorable for the first species. The detection of P. sergenti and P. papatasi was sporadic throughout the study probably because the ecological variables of the study area are very uniform and do not favor their presence. Sand fly density in each station showed variation between survey periods. The high growth of the rabbit population observed in BOS during the three surveys could influence the high sand fly abundance observed in this station. With the exception of this location, captures tended to decrease or maintain, possibly as a result of the control measures taken in the outbreak area such as management actions and regulation of hare and rabbit densities. Furthermore, studies carried out in France state that an area with a mean sand fly density above 20 sand flies/m2 is a zone with high leishmaniasis transmission . According to this statement, the area studied in this work would be a zone at very high risk of transmission due to the elevated P. perniciosus density registered (193.60 specimens/m2). Differences observed between ST and LT are due to the fact that the former method captures sand flies by interception and the latter by attraction. Thus, P. perniciosus are more phototropic than S. minuta, which explains the significant RA differences between the two species in the present study as previously reported . However, these results are in contrast to previous studies which recorded higher S. minuta density than P. perniciosus by using ST [36, 37]. A possible explanation could be the extraordinary P. perniciosus population present in the area due to a convergence of factors as an extraordinary availability of wild leporids and favorable bioclimatic conditions. Moreover, the ST were not placed in barbicans, drainage holes in walls, nor tubes; they were placed at the base of the trunks of shrubs or trees, next to deep cracks of the ground, next to wire fences, on the grates of sinkholes of sewer network, and in the proximities of the rabbit burrows [40–42].
Seasonal dynamics of P. perniciosus by LT showed a confluent bi-modal pattern in which the convergence of similar density peaks comprising June to September captures was observed (Fig. (Fig.3).3). However, a maximum peak of captures was found in August 2012, the month with lowest mean RH and highest mean temperature, which is in line with statistical analysis of correlation between captures and bioclimatic parameters. In the case of peaks of maximum number of captures observed in September 2013 and 2014 no significant correlation was observed, so other parameters such as wind, temperature and RH in the days leading up to sampling, or even the climatic conditions prevailing in the spring or winter, might have affected sand fly captures . Using ST we also observed a confluent bi-modal pattern during the three periods (Fig. (Fig.4).4). Just like with LT, 2012 results agreed with correlation analyses; however, 2013 and 2014 surveys did not. Confluent bi-modal pattern was also described in other studies in the Mediterranean basin , although it is very similar to the bi-modal pattern previously reported in Spain [43–45]. In case of S. minuta dynamics, LT captures were not significant enough to detect a pattern. However, using ST, the maximum peak in every period was in August, the month when highest temperatures and lowest RH were recorded, except in 2014 where S. minuta captures remained very high in September. Sergentomyia minuta collections agreed with the results obtained in the correlation analysis. This kind of correlation between RH and temperature has been previously reported in entomological studies carried out in the Mediterranean region .
Prevalence studies of Leishmania infection in the phlebotomine vectors are worthwhile indicator of the intensity of Leishmania transmission. In the present study, we combined classical demonstration of Leishmania promastigotes by dissection of sand fly guts and molecular detection and identification of the parasite using PCR. Although the microscopic detection of promastigotes in the gut of dissected sand flies is a very tedious task that requires highly skilled personnel, it allows the isolation and culture of the parasite and its preservation for further studies. It also enables the incrimination of sand fly vector species in Leishmania foci. In contrast, molecular procedures are very accurate and less laborious than sand fly dissection. Nevertheless, the combination of these two methodologies allows an extensive and comprehensive study of infection rates and characterization of Leishmania parasites as reported in Italy, Turkey and Israel [47–49]. Mean infection rates found by dissection in the present work are similar to some studies carried out in Spain [4, 50–53] but higher than reported in other countries [26, 47, 48, 54]. The infection rates found in this study showed different pattern during the three surveys, possibly due to differences observed in the number of sand fly captures. Specifically, the higher infection rate found in 2013 could be linked to the increase in the human leishmaniasis cases in early 2014. In addition, the highest rates of infection were observed in months with low number of sand fly captures and by the end of the seasonal activity of sand flies, as previously reported in other studies [50, 55]. The isolates in this work were characterized as ITS-Lombardi genotype. This ITS type is identical to that reported from direct xenodiagnoses carried out in hares and rabbits from the affected area [8, 9], as well as from samples obtained from patients from the outbreak .
On the other hand, the detection of Leishmania DNA by PCR shows high values in fed (13.35%) and unfed females (7.78%). These findings disagree with data we reported before , possibly because of the low number of specimens analyzed in the latter study, all of them captured in October 2011 when infection rates are usually the highest. On the other hand, the results from the present study agree with other studies where infection rates of fed female sand flies are higher than in unfed [47, 57, 58]. Moreover, an elevated infection rate was found in females fed on hares, mostly in the 2012 survey. These samples represent an important number of the fed females and could influence the elevated rate of positive samples in this group. Leishmania detection by molecular procedures showed higher values than average rates detected by sand fly dissection, exhibiting more sensibility and accuracy. Regarding fed females, the highest rates of infection were detected in 2012 and showed a decline pattern in the following periods. This may be a result of the control measures taken in the park close to urban areas to reduce the population of hares and rabbits, potential wild reservoirs of leishmaniasis. In this way, we could verify that in 2013 and 2014 there was an elevated rabbit density in BOS, which correlates with the high sand fly infection rates found in this station. On the other hand, the drastic reduction in the number of hares in the POL station resulted in a considerably low P. perniciosus infection rate in 2014 (Fig. (Fig.66).
Although the leishmaniasis is considered a zoonotic disease in Spain with the dog as a main reservoir, in the focus area where this study was carried out, the prevalence of canine leishmaniasis has not undergone any increase. However lagomorphs seem to play a role as wild reservoirs of leishmaniasis in the area close the focus, as shown by direct xenodiagnoses in hares and rabbits [8, 9]. In a previous study of blood meal preferences performed in a small number of blood-fed P. perniciosus from the focus had already been shown that most of them were fed on rabbits and hares . To find out more about host feeding preferences of sand flies collected in the park close the urban area of the focus, molecular procedures were applied. Global results of blood meal identification showed that P. perniciosus preferably feeds on rabbits, followed by hares, reinforcing the role these lagomorphs may be playing in the focus. Several studies on feeding habits of sand flies in Leishmania foci have reported an opportunistic behavior, finding engorged P. perniciosus females fed on different vertebrates, including humans, as well as mixed feeding sources [19, 47, 58]. In this way, a PCR restriction fragment length polymorphism methodology was optimized in order to discriminate mixed blood sources on sand flies captured in different areas of Spain . Regarding the efficiency of blood-meal identification in P. perniciosus females, this was similar or higher to that reported by other authors using the same methodology [17, 26, 58, 60]. The results obtained in the different stations showed that in BOS rabbit was the prevalent blood meal source, possibly due to the elevated population of this mammal observed in this station during the surveys. As early as 2011, the Community of Madrid warned about the existence of an extraordinary hare population in POL. Blood meal preferences found in this station correlate with this observation. However, the number of P. perniciosus fed on hares has been decreasing along the surveys probably due to the different measures implemented by the regional government in the frame of the control program of the disease. The host-feeding preferences found in this work are in concordance with studies which show that these two lagomorphs are highly exposed to P. perniciosus bites in the focus area . In addition, sand flies fed on cats were found in JIC and POL, were some colonies of these animals were observed. The variation in blood-feeding preferences observed correlates with the opportunistic feeding behavior of P. perniciosus (Fig. (Fig.6)6) .
We are therefore facing an area with a high L. infantum transmission where intensive and continuous surveillance should be taken to regulate lagomorph populations and sand fly vector densities. Information about vectors and their distribution should be updated regularly in order to manage efficient control programs of leishmaniasis in the area. Moreover, surveillance programs of leishmaniasis should be undertaken in other geographic regions endemic for leishmaniasis with similar eco-epidemiological characteristics in order to keep track of the populations of sand fly vectors, the presence of the parasite in the sand flies and the potential wild reservoirs like the leporids.
Our entomological study carried out in the green park close to the urban areas of southwestern region of Madrid (Spain) affected by the outbreak of leishmaniasis revealed high densities of P. perniciosus, the only vector in the area. Infection rates by dissection of sand fly guts and by molecular procedures showed an important presence of L. infantum in the vector. From our study, it was also observed that rabbits and hares were the main blood meal sources of this sand fly species confirming the essential role of lagomorphs in sustaining the sylvatic Leishmania cycle in the green park.
Monthly relation of P-values of P. perniciosus and S. minuta densities through the three survey periods. Results of Kruskal-Wallis and Dunn’s multiple comparison tests. (DOCX 18 kb)
Statistical analysis of RH through the three periods. Coefficient results from Kruskal-Wallis test and P-values resulting from Dunn’s multiple comparison test. (DOCX 18 kb)
Statistical analysis of temperature through the three periods of the study. Coefficient results from Kruskal-Wallis test and P-values resulting from Dunn’s multiple comparison test. (DOCX 20 kb)
Correlation analysis between P. perniciosus and S. minuta collections and bioclimatic parameters. Statistical results. (DOCX 20 kb)
Alignment of ITS sequences of Leishmania infantum isolates from 2012, 2013 and 2014 and strains with different ITS types of L. infantum retrieved in the Genbank: type Lombardi (AJ000295), type A (AJ634341), type B (AJ000288), type A/B (AJ634355), type E (AJ634361) and type F (AJ634370). (PDF 1192 kb)
This study was partially sponsored and funded by: Dirección General de Salud Pública, Consejería de Sanidad, Comunidad de Madrid; Colegio de Veterinarios de Madrid; Colegio de Biólogos de Madrid and EU grant FP7-261504 EDENext (http://www.edenext.eu). The contents of this publication are the sole responsibility of the authors and do not necessarily reflect the views of the European Commission.
Estela González was supported by the Directorate-General for Public Health of the Community of Madrid. Instituto de Salud Carlos III, provided the equipment and facilities. European Union grant Frame Program7 EDENext (7-2011-261,504) supported the design, collection, analysis and interpretation of data and writing the manuscript.
All data generated or analyzed during this study are included in this published article and its additional files.
MJ and RM conceived, designed and coordinated the study. MJ, RM and EG carried out the field work. EG, MJ, RM and IMM performed the dissections and identifications of sand flies and the isolation of Leishmania isolates. EG and SH processed the traps and identified collected sand flies. EG, MJ and SH carried out the molecular studies in sand flies. EG, MJ and IMM performed the characterization of Leishmania isolates. EG carried out the statistical analysis of the data and wrote the first draft of the manuscript. EG, MJ and RM wrote the final draft of the manuscript. All authors read and approved the final manuscript.
The authors declare that they have no competing interests.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Electronic supplementary material
The online version of this article (doi:10.1186/s13071-017-2309-z) contains supplementary material, which is available to authorized users.
Estela González, Email: se.iiicsi@refnoge.
Maribel Jiménez, Email: se.iiicsi@zenemijm.
Sonia Hernández, Email: firstname.lastname@example.org.
Inés Martín-Martín, Email: moc.liamtoh@seninitram.
Ricardo Molina, Email: se.iiicsi@anilomr.