Skip to main content

Advertisement

Leishmaniasis sand fly vector density reduction is less marked in destitute housing after insecticide thermal fogging

Abstract

Background

Insecticide thermal fogging (ITF) is a tool to control vector borne diseases. Insecticide application success for vector control has been associated with housing materials and architecture. Vector abundance is correlated with weather changes. Nevertheless, housing quality and weather impacts on vector abundance have been unaccounted for in most New World insecticide control trials for leishmaniasis vectors.

Methods

We conducted a 15 month insecticide control trial that included two deltamethrin [6 mg a.i.m-2] based ITF interventions in 12 of 24 monitored houses at Trinidad de Las Minas, a hyperendemic cutaneous leishmaniasis transmission village in western Panamá. During the study we followed sand fly (SF) abundance, keeping track of rainfall and quantified housing quality using an index based on architecture and construction materials.

Results

We found a 50 to 80% reduction in SF density in the fogged houses when compared with control houses, while controlling for seasonal changes in SF abundance associated with rainfall. We found heterogeneities in the reductions, as abundance changed according to SF species: Lutzomyia gomezi, Lu. panamensis, Lu. dysponeta and Lu. triramula reduced in density between 40% and 90% after ITF. In contrast, Lu. trapidoi density increased 5% after ITF. Differences in the impact of ITF were associated with housing quality, the most destitute houses, i.e., those with features that ease insect entrance, had a disproportionally larger SF abundance, in some cases with increased domiciliary SF density following the ITF.

Conclusion

Our results suggest the potential of insecticide application to control SF density and leishmaniasis transmission could depend on housing quality beyond insecticide efficiency.

Background

Cutaneous Leishmaniasis (CL) is a major neglected tropical disease worldwide [1, 2]. In relative terms, CL in the new world remains a poorly studied disease [2, 3], as are Sand Fly (SF) vectors of Leishmania spp parasites [46]. For example, studies on SF ecology and control are somewhat scarce [57] when compared with other vectors of pathogens, especially mosquitoes [810]. Long-term studies looking at the population dynamics of SF vectors [1118], and studies on insecticide based control of SF abundance [7] are limited. Insecticide control trials for New World sand flies have been carried out using organophosphates: Malathion [19]; organochlorines: DDT [20, 21] and pyrethroids: cyfluthrin [22], cypermethrin [23, 24], deltamethrin [2530], lambdacyhalothrin [3133]. Methods of insecticide application for SF control have included: fogging [19, 22, 34], residual spraying [21, 2327, 3133], insecticide diffusion devices [35], insecticide treated nets [28, 29] and curtains [30]. In general, these studies have shown that insecticide application at the household level seems to be the most effective SF control method [7], because it can suppress both SF abundance [7] and CL transmission over short [36] and long term [37] time scales. However, a key insight from the few long-term studies on SF population dynamics is that vector abundance is extremely sensitive to environmental changes, with different species having a distinctive sensitivity to particular meteorological components [12, 15, 17]. Nevertheless, most insecticide control trials have ignored the importance of weather variability when evaluating the efficacy of insecticides on SF abundance [2035].

CL has also been recognized as a disease of poverty [38, 39]. For example, our work has shown that in Central America CL primarily affects socially excluded populations, i.e., those who lack or have restricted access to resources that ensure a satisfactory quality of life [40, 41]. This fact is extremely important to better understand the dynamics of disease transmission. Because CL transmission is ruled by ecological processes [42], how humans become part of a pathogen’s life cycle, especially in neglected tropical diseases like leishmaniasis, can be heavily influenced by social factors [43]. Even if recognized as an important factor for the success of insecticide applications in CL control [23, 32], housing quality has been unaccounted for when evaluating the outcome of insecticide trials.

In Panamá and Costa Rica, clinical treatment of patients with skin lesions is the main activity pursued to control CL transmission [44, 45], neglecting active surveillance and SF vector control. Nevertheless, the disease is becoming increasingly common in both countries. For example, a recent urban CL outbreak occurred in Tilarán, Guanacaste, Costa Rica [44]. Similarly, a CL epidemic in Western Panamá province, Panamá, where 500 new CL cases (~50% of them in children < 5 years of age) were officially reported, representing a two-fold increase for this area when compared with recent years [45]. These outbreaks suggest that changes are taking place in CL transmission epidemiology in Central America, calling for the implementation of vector control as a major strategy to reduce the impact of CL on populations at risk of acquiring the disease [7, 41]. Here, we present results from a small scale SF insecticide control trial carried out at Trinidad de Las Minas, Western Panamá Province, República de Panamá. Our goal was to evaluate the impact of household insecticide thermal fogging (ITF) with deltamethrin (6 mg a.i.m-2) on domiciliary and peridomiciliary SF abundance. Baseline observations showed that household human infection rates were positively associated with the household abundance of Lutzomyia gomezi, but also with Lu. panamenesis and Lu. trapidoi[45], the most abundant vector species in the area [46]. For 15 months we followed domiciliary and peridomiciliary SF abundance in 24 houses, 12 were subjected to two rounds, 6 months apart in time, of Insecticide Thermal Fogging (ITF) with deltamethrin at 6 mg a.i.m-2 and 12 kept as control. We found a 50 to 80% reduction in SF density in the fogged houses when compared with control houses, while controlling for seasonal changes in SF abundance associated with rainfall. We observed a variety of species specific abundance changes, with Lu. gomezi, Lu. panamensis, Lu. dysponeta and Lu. triramula reducing their density between 40% and 90% after ITF, in contrast to Lu. trapidoi whose density increased 5% after the ITF. Spatially, we found that heterogeneities in SF abundance after ITF were associated with housing quality, specifically destitute houses, i.e., those with features that ease insect entrance, had the largest share of SF individuals, in some cases with an increased domiciliary SF density following the ITF. Thus, our results call for a better quantification and understanding of housing quality as a major factor underpinning the success of insecticide control for leishmaniasis vectors.

Methods

Study Area

Our study was conducted at Trinidad de Las Minas, (8°46’32‘N and 79°59’45‘W), a rural village in Capira District, western Panamá Province, República de Panamá. This village is 230 meters above sea level, with an annual mean temperature of 26.0°C and monthly rainfall ranging from 28 — 570 mm3. Climate is markedly seasonal, with a dry season from mid December to March and a rainy season for the rest of the year. The area used to be a lowland tropical moist forest, but currently is a transitional forest/agricultural matrix, with scattered deciduous and xerophile species. Further details about the study site are presented by Calzada et al.[46]. Daily rainfall records for our study period were obtained from a meteorological station within a 5 Km radius from our study site, managed by Panamá’s electrical company ETESA.

Insecticide thermal fogging

We selected 24 houses for our study (out of 128 houses in the village), where residents provided informed consent to collect sand flies inside and outside their houses. Twelve houses were subjected to indoor and outdoor insecticide thermal fogging, while the remaining 12 houses were kept as controls (no fogging). The number of houses evaluated in this study was primarily limited by the availability of sampling resources, especially the number of light traps. Nevertheless, based on the only previous study on SF control by fogging in Panamá [19] a sample of 12 houses for the insecticide treatment and 12 control houses is powerful enough (1-β > 0.80, given α < 0.05) to detect differences in sand fly abundance due to insecticide thermal fogging (ITF), using a generalized linear model with one treatment and one covariate (See Additional file 1: Protocol S1 and Additional file 2: Figure S1). We selected a cluster of houses with homogeneous eco-epidemiological conditions and where intra-domiciliary SF presence was confirmed by residents. Although we planned to match houses based on construction materials and architecture, our selection was limited by the lack of informed consent from some residents for the fogging, especially in the houses with the best construction materials. However, to counter this limitation, we carefully recorded details about housing construction and the peridomicile and vegetation of each house (see section below on Housing destituteness, the peridomiciliary environment, peridomiciliary vegetation structure and animal abundance/richness assessment) that were considered during the statistical analysis. We evaluated two rounds of indoor/outdoor ITF using deltamethrin (K-Othrine® 2.7 UBV, Bayer, Guatemala). Insecticide selection and application was performed by trained personnel of the Vector Control Department from Panamá´s Ministry of Health. For Insecticide selection, results from toxicity assays in naïve populations were considered [32, 4749], especially given that our study site has never been subjected to deltamethrin application for vector control. It is also worth highlighting that agriculture is primarily organic at Trinidad de Las Minas. The insecticide applications were conducted on July 18, 2010 and January 23, 2011. The insecticide was applied with a hand-held thermal fogger (Golden EagleTM, Model # 2610, Curtis Dyna-Fog Ltd, Westfield, IN, USA) to interior and exterior housing walls, targeting cracks and crevices. A similar fogging was performed in the 15 m around the houses (peridomicile). We choose the 15 m radius for the fogging based on studies on New World SF dispersal in which SF rarely travel beyond 50 m from a release point [5052]. On average, 0.57 L of insecticide (diluted in diesel to a final concentration of 0.7 g /L, following Panamá Ministry of Health guidelines) was used for the fogging of internal and external house walls (whose total wall surface on average was 65 m2), corresponding to a concentration of 6 mg of active ingredient per square meter (mg a.i.m-2).

Ethical clearance

This study was approved by the National Review Board, Comité Nacional de Bioética de la Investigación, Instituto Conmemorativo Gorgas de Estudios de la Salud, Ciudad de Panamá, República de Panamá (561 /CNBI/ICGES/06).

Sand fly (SF) abundance

We evaluated ITF impacts on SF abundance by comparing collections from the domicile and peridomicile of fogged and control houses. Sand flies were collected using modified light-traps [53]. Each trap was slightly modified by attaching an additional small LED light to increase SF attraction [46]. Entomological sampling was carried out monthly from April 2010 to June 2011, except for the months of August and November 2010 and January 2011, when access to this remote village was impossible because of logistical and operational constraints, which, in January 2011, prevented the sampling of houses just before the 2nd ITF. Thus, a total of 12 sampling surveys were conducted during the study. For each monthly collection, one trap was placed for one night in the main bedroom of every household (indoor). This trap was suspended from the ceiling at about 2 m from the ground floor. Another trap was placed at the same height, above vegetation, within 50 meters of the house (i.e., peridomicile). Traps were setup for 12 hours, from 6:00 pm to 6:00 am, in the same position (indoor and peridomicile) during each sampling session.

Trapped sand flies were removed from the traps, stored at −20°C to kill the insects and preserved in 70% ethanol for identification. For each trap, we summarized the abundance, sex and species of sand flies following Young and Duncan [54], with male genitalia and female spermathecae as main diagnostic taxonomic characters. A detailed description of the SF fauna species composition at Trinidad de Las Minas and changes following the ITF is presented by Calzada et al.[46].

Housing destituteness, the peridomiciliary environment, peridomiciliary vegetation structure and animal abundance/richness assessment

For each house, we collected data on construction materials and the presence of insect friendly gaps to construct an index of house quality, hereafter referred to as housing destituteness index, HP. We specifically collected data on materials used for walls and roofs, whether the floor was earthen or covered with concrete or wood. We recorded the presence of insect friendly gaps, i.e.; whether walls had crevices, if walls were complete from the base to the roof, whether doors had holes and if the houses had windows and whether the windows had anti-insect screening. We recorded the presence/absence of elements that may serve as peridomicilary resting places for sand flies, such as rubbish, fruit trees, etc. [12, 55] to estimate a peridomicile index (PI). We measured several elements of the vegetation structure in the peridomicile for the estimation of a vegetation index (VI). To quantify any possible role that vertebrate host abundance could have on sand fly density, we performed a census on all domestic animals belonging to each household and compiled a list of wildlife species seen by household residents in the domiciliary/peridomiciliary area of each household. We used these data to estimate animal abundance indices. HP, PI, VI and the animal abundance indices were estimated by computing the first principal component for the set of variables considered in each index, for further details about the variables considered for each index, and the principal components analysis implementation see Additional file 3: Protocol S2. For both domestic and wild animals we also estimated species richness at each household by, respectively, counting the number of domestic species, recorded in the domestic animal census, or the number of reported wildlife species seen by the householders.

Statistical analysis

For the analysis we employed a two-fold strategy, we evaluated the impact of the fogging both temporally and spatially. For the temporal analysis we employed negative binomial generalized linear models (NB-GLM) that accounted for overdispersion in SF counts [56]. As a first approach, models considered all SF species abundance, separating groups according to the feeding habit i.e., whether sand flies have been recorded biting humans or not (i.e., anthropophilic or zoophilic) [42], and sampling habitat (peridomicile and domicile). We also considered the monthly average and the standard deviation, S.D., of daily rainfall, with a one month lag, the lag selected with a cross-correlation analysis [57], as a covariate to account for seasonal abundance fluctuations in SF abundance and whether a house was fogged or not. We included rainfall S.D. as covariate to account for rainfall variability as an underpinning factor of SF abundance, specifically SF abundance could be sensitive not only to the amount of rain during a time period, a quantity measured by the average rainfall, but also to the intensity of the rainfall events [12], a quantity measured by the S.D. For the insecticide we considered whether fogging independently of the application date had an impact on SF density (Fogging A in the models) or whether the fogging impact was different for each of the two fogging rounds (Fogging B in the models). We further developed models for the three main dominant vectors in our study site: Lu. gomezi, Lu. panamensis and Lu. trapidoi (which accounted for ~60% of the collected sand flies, Table 1 and [46]) and for the two most abundant non vector-species Lu. triramula and Lu. dysponeta (Table 1); comparing temporal dynamics in the domicile and peridomicile of fogged and control houses. For the spatial analysis we also employed NB-GLMs. In these models we looked at the cumulative number of sand flies caught after the foggings as a function of housing destituteness (HP), the peridomiciliary environment index (PI), the vegetation structure index (VI) and the animal abundance/richness indices. We considered possible non-linearities in the association of HP, PI, VI and the animal abundance/richness indices with SF abundance by fitting polynomials and models with breakpoints [40], i.e., threshold values in the independent variables at which the association with a dependent variable can change quantitatively, employing hockey stick regressions (see Additional file 1: Protocol S1 for further details). For both the temporal and spatial models we performed model selection based on the Akaike Information Criterion, a tool for model selection that guides the choice of a “best model” based on its likelihood to explain the data with a minimum number of parameters [56]. For the best temporal model, beyond the assumption of the NB-GLM, we tested the temporal independence of the residuals (error) through the inspection of the autocorrelation function, ACF, and the partial autocorrelation function, PACF [57]. For the best spatial model we tested the spatial independence of the residuals with the Moran’s I test.

Table 1 Sand Fly species abundance in the control and fogged houses

Results

Sand fly fauna description

We collected 5628 sand flies in the 576 sampling night-traps (Table 1). We were unable to identify 11 individuals, all other individuals belonged to 24 SF species (full details have been published elsewhere [46]). Of the collected species, 8 had a proven vector status [41, 42] and 3 accounted for more than 50% of the samples: Lu. trapidoi (20%) Lu. gomezi (20%) and Lu. panamensis (17%). Of the remaining SF species, Lu. triramula (20%) and Lu. dysponeta (8.7%) were the most abundant. Sand flies were more abundant outside (peridomicile) than inside (domicile) the studied houses (58% vs 42%). Table 1 also shows that total abundance of sand flies was reduced by approximately 40% inside the fogged houses (i.e., domicile; control: 1484, fogged: 881) and close to 50% in the peridomiciliary environments (control: 2157, fogged: 1106). Further details are presented by Calzada et al. [46]

Temporal impacts of insecticide thermal fogging on all Sand Fly species

Additional file 4: Table S1 shows the different models considered to explain the population dynamics of all the SF species in the domicile (Figure 1A) and peridomicile (Figure 1B). Rainfall (Figure 1C) seemed to be an important factor shaping SF population dynamics, so was insecticide fogging (Figure 1D). The best model found statistically significant differences in SF abundance according to vectorial status (vector, i.e., anthropophilic or non-vector, i.e., zoophilic), environment (domicile and peridomicile), rainfall (a second degree polynomial), fogged vs control houses, as well as the interaction of vectorial status with the fogging and the environment (Additional file 4: Table S1). Quantitatively, the best model (Table 2) indicates that non-vector (zoophilic) sand flies were most abundant in the peridomiciliary environments (on average 308/per house and sampling night), followed by vector (anthropophilic) SF species (about 22% less individuals than non-vectors in the domicile). Inside the houses (domicile environment), vector species (anthropophilic) were the most abundant (about 78% of the zoophilic sand flies observed in the peridomicile). Fogging reduced the abundance of vector sand flies inside the houses by 51%, and 77% in the peridomicile (when compared with non-vector species in the peridomicile). Rainfall variability, measured by Rainfall S.D., had a concave relationship with SF abundance, i.e., a second degree polynominal with a minimum, in other words SF abundance was higher at low and high levels of rainfall variability, a feature that could reflect the different responses of SF to environmental variability. As suggested by Chaniotis et al.[12] some sand flies might become very abundant when rainfall is sustained, while others might thrive when there is a high alternancy of dry and wet periods. All the assumptions of the GLM model were not violated and model residuals were not autocorrelated. In general, all the estimated parameters accurately described the patterns observed in Figure 1.

Figure 1
figure1

Vector Abundance and Rainfall. (A) Domiciliary (Dom) vector species abundance, vertical lines indicate the timing of the Foggings (B) Peridomiciliary (Per) vector species abundance (C) Monthly Rainfall, daily mean and S.D. (D) Sand Fly abundance as a function of monthly S.D. of daily Rainfall records, symbols are used to separate domiciliary and peridomiciliary species from control and fogged houses. In (A) and (B) the open circles indicate the measured values, where each point represents the cumulative abundance from the 12 trap-nights per treatment.

Table 2 Parameter estimates for a negative binomial model explaining the abundance of all Phlebotomine sand fly species in Trinidad de las Minas, Capira District, Panamá

Temporal impacts of insecticide thermal fogging on the most abundant Sand Fly species

Population dynamics of the five most abundant species can be observed in Figure 2. Figure 2A and 2B show, respectively, domicile abundance of fogged and control houses. Figure 2C shows peridomicile abundance of fogged houses, while Figure 2D shows peridomicile abundance of control houses. Inside the houses (domicile) the most abundant species were Lu. trapidoi and Lu. gomezi with a higher abundance in the control houses. In the peridomicile of the fogged houses Lu. trapidoi and Lu. panamensis were the most abundant species, while in the control houses, prior to the foggings, Lu. triramula was the most abundant species. Additional file 5: Table S2 presents the selection for the best model explaining the dynamics of the 5 most abundant species in our samples. All species had significant differences in their domicile and peridomicile abundance, all were sensitive to the fogging, Lu. trapidoi and Lu. dysponeta were sensitive to the synergistic (multiplicative) effects of monthly average daily rainfall and its standard deviation, meaning that these species could be sensitive not only to total rainfall but also to how variable the rainfall was during a time period, while Lu. gomezi, Lu. panamensis and Lu. triramula were sensitive to monthly S.D. of daily rainfall, in a non-linear relationship described by a second degree polynomial, which indicates these species thrive when there is sustained rain levels or when there is a high variability in the rain levels. In general terms all species were positively associated with rainfall, i.e., their abundance increased following increases in Rainfall. Lu. triramula was the most sensitive species to the insecticide fogging (Table 3), reducing its abundance by up to 91%, followed by Lu. gomezi (68%), Lu. dysponeta (51%), Lu. panamensis (40%). In contrast, Lu. trapidoi increased its abundance by 5% (Table 3).

Figure 2
figure2

Common species abundance patterns. (A) Domiciliary environment of fogged houses (B) Domicilary environment of control houses (C) Peridomiciliary environment of fogged houses (D) Peridomicilary environment of control houses. In the panels open circles indicate the measured values, each point represents the total cumulative sand fly abundance from the 12 trap-nights per treatment (i.e., 1 trap-night/house) and different species are indicated by different lines, see inset legend of Panel (A) for details.

Table 3 Parameter estimates for the negative binomial models explaining the abundance of selected Phlebotomine sand fly species in Trinidad de Las Minas, Capira, Panamá

Spatial impacts of insecticide thermal fogging on Sand Fly abundance: the role of destitute housing

Figure 3A shows the spatial patterns of housing destituteness, HP, which can be interpreted as a weighted average of the different components considered for the index (Additional file 6: Table S3). In general, high scores indicative of housing destituteness, i.e., features that could ease SF entrance into houses, resulted in a high HP index. Spatial patterns of SF infestation per trap the night before (Figure 3B) and after the first (Figure 3C) and second fogging (Figure 3D) in general suggest that, especially after the foggings, large infestations were positively associated with destitute housing conditions (Figures 3 and 4), but not with the peridomicile environment (Additional file 7: Figure S2, for index interpretation see Additional file 8: Table S4) the vegetation structure (Additional file 9: Figure S3, for index interpretation see Additional file 10: Table S5) or the animal abundance indices (Additional file 11: Figure S4, for index interpretation see Additional file 12: Table S6), or household residents (Additional file 11: Figure S4L). Nevertheless, the larger infestations before fogging were associated with a high peridomicile index (Additional file 7: Figure S2A), i.e., houses that had plenty of resting sites that were neither ornamental trees nor vegetable crops (Additional file 8: Table S4). In general, infestations seemed to be unrelated with whether the vegetation was tall, i.e.; high Vegetation Index, VI values, or ground cover i.e., low VI values (Additional file 10: Table S5, Additional file 9: Figure S3). Figure 4A shows the association between cumulative SF abundance prior to the fogging and HP. Figure 4B shows the association between cumulative SF abundance after the foggings and HP, which shows a clear pattern were the houses with the highest HP, on average, had the highest infestation. The curve shown in Figure 4B is the fit from a hockey stick regression (Table 4), which showed that, in general, after the foggings the highest infestations were observed in the most destitute houses, with sandflies increasing about ten times (11.6) for each 0.1 increase in HP above 0.586. For values below 0.586 infestation levels were homogeneous, with an average of 90 sand flies over 9 night-traps. The model presented in Table 4 had residuals that followed the assumptions of a NB-GLM and were not spatially correlated according to a Moran’s I test (I = 0.11; P > 0.18). The model in Table 4 was also selected from a larger ensemble of models that considered fogging, HP, PI, VI and the animal abundance indices monotonically associated with SF abundance, as well as, HP interactions with other variables (Additional file 13: Table S7). Finally, Figure 4C and 4D show the difference in the number of sand flies captured in the domicile and peridomicile, respectively, before and after the fogging of the studied households. Figure 4C shows that inside the houses, in general, SF density decreased or stayed the same, but for high HP, SF abundance tended to increase after the foggings. In contrast, Figure 4D shows that in the peridomicile environment SF abundance decreased independent of the degree of HP.

Figure 3
figure3

Spatial patterns. (A) Housing destituteness (HP) index. Circle size is proportional to the HP index, in the inset legend C (grey) stands for control and F (black) for fogged. Average sand fly density per house and trap night (B) Before the Fogging (C) After the 1st Fogging (D) After the 2nd Fogging. In (B), (C) and (D) circle size is proportional to the abundance of sand flies, see inset legend in (B) for reference (25 Sand Flies/trap night), and values were standardized by dividing the total cumulative abundance of each period by the number of trap nights in each period, i.e., 3, 4 and 5 respectively.

Figure 4
figure4

Housing destituteness and sand fly abundance. (A) Cumulative sand fly abundance before the fogging as a function of Housing Destituteness, HP (three traps nights) (B) Cumulative sand fly abundance for the nine traps nights that contained the two interventions as a function of HP. The solid black line is the fit from a negative binomial model with a breakpoint (C) Changes in domiciliary sand fly density per house and trap night as function of HP (D) Changes in peridomiciliary sand fly density per house and trap night as a function of HP. In panels (A) to (D) symbols indicate whether a house was intervened, triangles, or not, circles. In panels (C) and (D) black symbols represent the pre-intervention densities and green symbols post-intervention densities. To ease the tracking of changes in each house we joined the pre and post intervention densities with a line whose color is blue when sand fly density decreased and red when sand fly density increased. In panel (C), the solid black line represents Table 4 model fit.

Table 4 Parameter estimates for a negative binomial model explaining post intervention heterogeneities in sand fly abundance across the households

Discussion

The leishmaniases primarily affect poor people worldwide [39], but why poor populations are the most susceptible to CL infection in the New World requires a better understanding of how socio-economic human conditions alter CL transmission ecology [40, 41]. SF abundance plays a major role in CL transmission, baseline results from our study area [45] and the whole República de Panamá [41, 42, 58] have shown a positive association between domiciliary vector abundance and CL human infections. Thus, our results suggest that housing quality, a socio-economic and ecological factor, could be a major entomological risk determinant for human infection with Leishmania spp parasites in rural Panamá. Our data clearly show how, beyond the transient impacts that insecticide applications could have on SF vector abundance, differences in housing quality may drive heterogeneities in SF house infestation, supporting suggestions from earlier SF insecticide control trials [23, 32], which highlighted the role that poor housing quality could have on hampering SF control by insecticides. We found that houses with features that ease insect entrance, such as lack of anti-insect screening or with abundant holes, that were built with poor quality building materials (e.g., mud, which is more likely to generate places that can serve as resting sites for sand flies), or materials that can shorten the active lifespan of insecticides (e.g., wood or mud as opposed to concrete [23, 26, 32]), had the largest share of sand flies after the ITFs. In that sense, our results encourage further entomological research to understand exactly what factors make a household more susceptible to the infestation by sand flies. It will be ideal to devise simple housing modifications that can reduce the likelihood of SF household entry, as has been carried out with malaria vectors [59]. Also, a better understanding of the impact of peridomiciliary insecticide application on the recruitment of new adult sand flies into households requires more detailed study. As observed by Perich et al.[22], large scale barrier spraying can significantly reduce the number of sand flies entering an area. A similar decline in the number of sand flies inside the houses after an indoor plus peridomiciliary spraying was observed in Brazil [26] and Bolivia [25]. Nevertheless, it is not clear if there is an additional impact on SF abundance via a decrease in the recruitment of adults that may emerge from the peridomiciliary ground, given its potential to harbor SF larvae, either in tree trunks or the open ground [6065] as the insecticide may also act upon SF larvae. For example, here we did not observe a rebound of peridomiciliary SF abundance after the 1st ITF, and this may be related to the recruitment of adults from larvae. However, it may also be related to eliminating a major source/refuge of adult sand flies with the ITFs [33].

Regarding the temporal impacts of deltamethrin on SF abundance in the fogged houses, insecticide action seemed to last up to 4 months, a shorter effect than what has been reported for higher concentrations of this insecticide when applied at 4-fold concentrations (~25 mg a.i.m-2), in residual spraying [25, 26] or 2 to 10-fold concentrations (~12.5 to 60 mg a.i.m-2), when used on impregnated nets and curtains [28, 29], where impacts seemed to last over 6 months. We think the difference may be a direct effect of insecticide concentration, which was significantly lower in our trial (~6 mg a.i.m-2). Also, unlike previous studies, we controlled for seasonality in SF abundance by incorporating rainfall in our temporal models, thus improving the ability to properly quantify the effect of the insecticide when compared with previous studies on SF insecticide control [2035]. Our observations also raise some questions about the ideal frequency for insecticide application. For example, the only previous study in Panamá applied insecticides every 2 months based on practices for mosquito control [19], but our study suggests that such frequency might be too high given the low reactivity, i.e., ability of sand flies to respond to changes in their population abundance, which, for example, can be very high for mosquitoes [66]. Nonetheless, more detailed analysis of SF species population dynamics are necessary to better understand what would be an ideal frequency for insecticide application based on the dynamical properties of SF populations and not merely on insecticide bioassays.

A point deserving further attention is the impact of insecticide application on the community of SF species, especially when several vectors co-occur [67]. Even though Lu. gomezi and Lu. panamensis, two major vectors of Leishmania panamensis, the CL parasite in our study area, significantly reduced their abundance (more than 50% less abundant after the ITFs), Lu. trapidoi, traditionally the major CL vector in Panamá [42], increased its abundance by 5% after the ITFs. This observation maybe a byproduct of externally reducing the abundance of other dominant vectors and/or reflect some unaccounted variability of species seasonality [12]. However, it illustrates the need to more carefully look at vector control impacts, not only on dominant species but on the whole community of vectors [46, 68, 69]. We also need to acknowledge that a better handling of SF seasonality could have been possible in our models if we had access to more climatic variables, but only rainfall was tracked at the only weather station in our study area and we did not have resources to track meteorological information on our own.

Finally, our study has some limitations due to its relatively small scale. We are unable to tell whether there was a community wide impact of the patchy ITF on 12 of the 24 houses, and a better house pairing according to the destituteness was impossible given the lack of consent. However, we are confident our analysis shows that at the household level there is a reduction of SF abundance following ITF and that heterogeneities in SF reduction are related to housing quality. Still, scaling up SF insecticide control trials could be a desirable first step to improve CL control as it could allow a better understanding of the role of: (i) peridomiciliary (ii) patchy vs uniform and (iii) differential insecticide applications [33]. However, from our own experience, this is a major challenge given the neglected nature of CL and the lack of interested partners (funders, insecticide developers) on making larger scale trials possible or improving CL control.

Conclusion

Our data clearly illustrate the importance of accounting for housing quality when evaluating insecticide control for leishmaniasis vectors and highlight the major role that destitute housing may have as a driving factor in the association between leishmaniasis and poverty.

Abbreviations

ITF:

Insecticide thermal fogging

SF:

Sand fly

CL:

Cutaneous leishmaniasis

HP:

Housing destituteness index

PI:

Peridomicile index

VI:

Vegetation index

NB-GLM:

Negative binomial generalized linear model

References

  1. 1.

    Alvar J, Vélez ID, Bern C, Herrero M, Desjeux P, Cano J, Jannin J, Boer M, the WHOLCT: Leishmaniasis worldwide and global estimates of its incidence. PLoS One. 2012, 7: e35671-10.1371/journal.pone.0035671.

  2. 2.

    Reithinger R, Dujardin J-C, Louzir H, Pirmez C, Alexander B, Brooker S: Cutaneous leishmaniasis. Lancet Infect Dis. 2007, 7: 581-596. 10.1016/S1473-3099(07)70209-8.

  3. 3.

    Brito MEF, Andrade MS, Dantas-Torres F, Rodrigues EHG, Cavalcanti MP, Almeida AMP, Brandão-Filho SP: Cutaneous leishmaniasis in northeastern Brazil: a critical appraisal of studies conducted in State of Pernambuco. Rev Soc Bras Med Trop. 2012, 45: 425-429. 10.1590/S0037-86822012005000006.

  4. 4.

    Kato H, Gomez E, Cáceres A, Uezato H, Mimori T, Hashiguchi Y: Molecular epidemiology for vector research on leishmaniasis. Int J Env Res Pub Health. 2010, 7: 814-826. 10.3390/ijerph7030814.

  5. 5.

    Ready PD: Biology of Phlebotomine Sand Flies as Vectors of Disease Agents. Ann Rev Entomol. 2013, 58: 227-250. 10.1146/annurev-ento-120811-153557.

  6. 6.

    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. 10.1111/j.1365-2915.2012.01034.x.

  7. 7.

    Alexander B, Maroli M: Control of phlebotomine sandflies. Med Vet Entomol. 2003, 17: 1-18. 10.1046/j.1365-2915.2003.00420.x.

  8. 8.

    Walker K, Lynch M: Contributions of Anopheles larval control to malaria suppression in tropical Africa: review of achievements and potential. Med Vet Entomol. 2007, 21: 2-21. 10.1111/j.1365-2915.2007.00674.x.

  9. 9.

    Silver JB: Mosquito ecology: field sampling methods. 2008, New York: Springer, 3

  10. 10.

    Service MW: Biological-control of mosquitos - has it a future. Mosq News. 1983, 43: 113-120.

  11. 11.

    Chaniotis BN, Anderson JR: Age structure, population dynamics and vector potential of Phlebotomus in Northern California Part II. Field population dynamics and natural flagellate infections in parous females. J Med Entomol. 1968, 5: 273-292.

  12. 12.

    Chaniotis BN, Neely JM, Correa MA, Tesh RB, Johnson KM: Natural population dynamics of Phlebotomine Sandflies in Panama. J Med Entomol. 1971, 8: 339-352.

  13. 13.

    Morrison AC, Ferro C, Pardo R, Torres M, Devlin B, Wilson ML, Tesh RB: Seasonal abundance of Lutzomyia longipalpis (Dipteral Psychodidae) at an endemic focus of visceral leishmaniasis in Colombia. J Med Entomol. 1995, 32: 538-548.

  14. 14.

    Salomon OD, Wilson ML, Munstermann LE, Travi BL: Spatial and temporal patterns of phlebotomine sand flies (Diptera: Psychodidae) in a cutaneous leishmaniasis focus in Northern Argentina. J Med Entomol. 2004, 41: 33-39. 10.1603/0022-2585-41.1.33.

  15. 15.

    Feliciangeli MD: Ecology of sandflies (Diptera: Psychodidae) in a restricted focus of cutaneous leishmaniasis in Northern Venezuela: III. Seasonal fluctuation. Mem Inst Oswaldo Cruz. 1987, 82: 167-176. 10.1590/S0074-02761987000200003.

  16. 16.

    Feliciangeli DM, Rabinovich J: Abundance of Lutzomyia ovallesi but not Lu. gomezi (Diptera: Psychodidae) correlated with cutaneous leishmaniasis incidence in north-central Venezuela. Med Vet Entomol. 1998, 12: 121-131. 10.1046/j.1365-2915.1998.00072.x.

  17. 17.

    Scorza JV, Gomez I, McLure MT, Ramirez M: Sobre las condiciones microclimaticas prevalentes en los microhabitats de flebotomos. Acta Biol Venez. 1968, 6: 1-27.

  18. 18.

    Marquez JC, Scorza JV: Dinamica poblacional de Lutzomyia townsendi (Ortiz, 1959) (Diptera: Psychodidae) y su paridad en Trujillo, Venezuela. Bol Dir Malariol San Amb. 1984, 24: 8-20.

  19. 19.

    Chaniotis BN, Parsons RE, Harlan HJ, Correa MA: A pilot study to control phlebotomine Sand Flies (Diptera: Psychodidae) in a Neotropical Rain Forest. J Med Entomol. 1982, 19: 1-5.

  20. 20.

    Deane LM, Deane МР, Alencar JE: Оbsеrvacоеs sоbге о combate ао Phlebotomus longipalpis реle dedetizacoes domiciliaria еn focos endemicos de Calazar по Сеага. Rev Bras Malariol Doen Trop. 1955, 7: 131-141.

  21. 21.

    Ready PD, Arias JR, Freitas RA: A pilot study to control Lutzomyia umbratilis (Diptera: Psychodidae), the major vector of Leishmania braziliensis guyanensis, in a peri-urban rainforest of Manaus, Amazonas State, Brazil. Mem Inst Oswaldo Cruz. 1985, 80: 27-36.

  22. 22.

    Perich MJ, Hoch AL, Rizzo N, Rowton ED: Insecticide barrier spraying for the control of sand fly vectors of cutaneous leishmaniasis in Rural Guatemala. Am J Trop Med Hyg. 1995, 52: 485-488.

  23. 23.

    Passerat De Silans LNM, Dedet J-P, Arias JR: Field monitoring of cypermethrin residual effect on the mortality rates of the phlebotomine sand fly Lutzomyia longipalpis in the State of Paraíba, Brazil. Mem Inst Oswaldo Cruz. 1998, 93: 339-344. 10.1590/S0074-02761998000300012.

  24. 24.

    Barata RA, Michalsky EM, Fujiwara RT, Franca-Silva JC, Rocha MF, Dias ES: Assessment of sand fly (Diptera, Psychodidae) control using cypermethrin in an endemic area for visceral leishmaniasis, Montes Claros, Minas Gerais State, Brazil. Cad Saude Publica. 2011, 27: 2116-2122.

  25. 25.

    Le Pont F, Padilla JM, Desjeux P, Richard A, Mouchet J: Impact de pulverisations de deltamethrine dans un foyer de leishmaniose de Bolivie. Ann Soc Belge Med Trop. 1989, 69: 223-232.

  26. 26.

    Falcão AL, Falcão AR, Pinto CT, Gontijo CMF, Falqueto A: Effect of deltamethrin spraying on the sandfly populations in a focus of American cutaneous leismaniasis. Mem Inst Oswaldo Cruz. 1991, 86: 399-404. 10.1590/S0074-02761991000400004.

  27. 27.

    Alexander B, Jaramillo C, Usma MC, Quesada BL, Cadena H, Roa W, Travi BL: An attempt to control Phlebotomine sand flies (Diptera: Psychodidae) by residual spraying with deltamethrin in a Colombian village. Mem Inst Oswaldo Cruz. 1995, 90: 421-424. 10.1590/S0074-02761995000300020.

  28. 28.

    Alexander B, Usma MC, Cadena H, Quesada BL, Solarte Y, Roa W, Travi BL: Evaluation of deltamethrin-impregnated bednets and curtains against phlebotomine sandflies in Valle del Cauca, Colombia. Med Vet Entomol. 1995, 9: 279-283. 10.1111/j.1365-2915.1995.tb00134.x.

  29. 29.

    Feliciangeli MD, Maroli M, Wheeler A, Townson H, Ward R, Maignon R: Sandfly control trial with deltamethrin impregnated curtains in El Ingenio, Miranda State, Venezuela. Bol Dir Malariol San Amb. 1995, 35: 127-132.

  30. 30.

    Kroeger A, Avila EV, Morison L: Insecticide impregnated curtains to control domestic transmission of cutaneous leishmaniasis in Venezuela: cluster randomised trial. BMJ. 2002, 325: 810-813. 10.1136/bmj.325.7368.810.

  31. 31.

    Davies CR, Llanos-Cuentas EA, Campos P, Monge J, Leon E, Canales J: Spraying houses in the Peruvian Andes with lambda-cyhalothrin protects residents against cutaneous leishmaniasis. Trans Roy Soc Trop Med Hyg. 2000, 94: 631-636. 10.1016/S0035-9203(00)90214-1.

  32. 32.

    Feliciangeli MD, Mazzarri MB, Campbell-Lendrum D, Maroli M, Maingon R: Cutaneous leishmaniasis vector control perspectives using lambdacyhalothrin residual house spraying in El Ingenio, Miranda State, Venezuela. Trans Roy Soc Trop Med Hyg. 2003, 97: 641-646. 10.1016/S0035-9203(03)80095-0.

  33. 33.

    Kelly DW, Mustafa Z, Dye C: Differential application of lambda-cyhalothrin to control the sandfly Lutzomyia longipalpis. Med Vet Entomol. 1997, 11: 13-24. 10.1111/j.1365-2915.1997.tb00285.x.

  34. 34.

    Britch SC, Linthicum KJ, Walker TW, Farooq M, Gordon SW, Clark JW, Ngere F, Ngonga D, Chepchieng C: Evaluation of ULV applications against old world sand fly (Diptera: Psychodidae) species in Equatorial Kenya. J Med Entomol. 2011, 48: 1145-1159. 10.1603/ME11025.

  35. 35.

    Sirak-Wizeman M, Faiman R, Al-Jawabreh A, Warburg A: Control of phlebotomine sandflies in confined spaces using diffusible repellents and insecticides. Med Vet Entomol. 2008, 22: 405-412. 10.1111/j.1365-2915.2008.00762.x.

  36. 36.

    Davies CR, Llanos-Cuentas EA, Pyke SDM, Dye C: Cutaneous leishmaniasis in the Peruvian Andes: an epidemiological study of infection and immunity. Epidemiol Infect. 1995, 114: 297-318. 10.1017/S0950268800057964.

  37. 37.

    Davies CR, Llanos-Cuentas A, Canales J, Leon E, Alvarez E, Monge J, Tolentino E, Gomero Q, Pyke S, Dye C: The fall and rise of Andean cutaneous leishmaniasis: transient impact of the DDT campaign in Peru. Trans Roy Soc Trop Med Hyg. 1994, 88: 389-393. 10.1016/0035-9203(94)90395-6.

  38. 38.

    Wijeyaratne PM, Arsenault LKJ, Murphy CJ: Endemic disease and development: the leishmaniases. Acta Trop. 1994, 56: 349-364. 10.1016/0001-706X(94)90106-6.

  39. 39.

    Alvar J, Yactayo S, Bern C: Leishmaniasis and poverty. Trends Parasitol. 2006, 22: 552-557. 10.1016/j.pt.2006.09.004.

  40. 40.

    Chaves LF, Cohen JM, Pascual M, Wilson ML: Social exclusion modifies climate and deforestation impacts on a vector-borne disease. PLoS Negl Trop Dis. 2008, 2: e176-10.1371/journal.pntd.0000176.

  41. 41.

    Miranda A, Carrasco R, Paz H, Pascale JM, Samudio F, Saldaña A, Santamaría G, Mendoza Y, Calzada JE: Molecular epidemiology of American Tegumentary Leishmaniasis in Panama. Am J Trop Med Hyg. 2009, 81: 565-571. 10.4269/ajtmh.2009.08-0265.

  42. 42.

    Christensen HA, Fairchild GB, Herrer A, Johnson CM, Young DG, Vasquez AM: The ecology of cutaneous leishmaniasis in the republic of Panama. J Med Entomol. 1983, 20: 463-484.

  43. 43.

    Levins R, Lopez C: Toward an ecosocial view of health. Int J Health Serv. 1999, 29: 261-293. 10.2190/WLVK-D0RR-KVBV-A1DH.

  44. 44.

    Jaramillo-Antillón O, Espinoza-Aguirre A, Lobo-Philp R: Estado actual de la leishmaniosis en Costa Rica. Acta Med Costarricense. 2009, 51: 158-164.

  45. 45.

    Saldaña A, Chaves LF, Rigg CA, Wald C, Smucker JE, Calzada JE: Clinical cutaneous leishmaniasis rates are associated with household Lutzomyia gomezi, Lu. panamensis, and Lu. trapidoi Abundance in Trinidad de Las Minas, Western Panama. Am J Trop Med Hyg. 2013, 88: 572-574. 10.4269/ajtmh.12-0579.

  46. 46.

    Calzada JE, Saldaña A, Rigg C, Valderrama A, Romero L, Chaves LF: Changes in phlebotomine sand fly species composition following insecticide thermal fogging in a rural setting of western Panamá. PLoS One. 2013, 8: e53289-10.1371/journal.pone.0053289.

  47. 47.

    Alexander B, Barros VC, De Souza SF, Barros SS, Teodoro LP, Soares ZR, Gontijo NF, Reithinger R: Susceptibility to chemical insecticides of two Brazilian populations of the visceral leishmaniasis vector Lutzomyia longipalpis (Diptera: Psychodidae). Trop Med Int Health. 2009, 14: 1272-1277. 10.1111/j.1365-3156.2009.02371.x.

  48. 48.

    Mazzarri MB, Feliciangeli MD, Maroli M, Hernandez A, Bravo A: Susceptibility of Lutzomyia longipalpis (Diptera : Psychodidae) to selected insecticides in an endemic focus of visceral leishmaniasis in Venezuela. J Am Mosq Control Assoc. 1997, 13: 335-341.

  49. 49.

    Scorza JV, Rosario CL, Scorza-Dagert JV, Rojas E: Susceptibilidad de hembras silvestres de Lutzomyia youngi de Trujillo, Venezuela, a insecticidas sinteticos. Bol Dir Malariol San Amb. 1995, 35: 311-326.

  50. 50.

    Chaniotis BN, Correa MA, Tesh RB, Johnson KM: Horizontal and vertical movements of phlebotomine sandflies in a Panamanian rain forest. J Med Entomol. 1974, 11: 369-375.

  51. 51.

    Morrison AC, Ferro C, Morales A, Tesh RB, Wilson ML: Dispersal of the sand fly Lutzomyia longipalpis (Diptera: Psychodidae) at an endemic focus of visceral leishmaniasis in Colombia. J Med Entomol. 1993, 30: 427-435.

  52. 52.

    Alexander JB: Dispersal of Phlebotomine Sand Flies (Diptera: Psychodidae) in a Colombian Coffee Plantation. J Med Entomol. 1987, 24: 552-558.

  53. 53.

    Pugedo H, Barata RA, França-Silva JC, Silva JC, Dias ES: HP: um modelo aprimorado de armadilha luminosa de sucção para a captura de pequenos insetos. Rev Soc Bras Med Trop. 2005, 38: 70-72.

  54. 54.

    Young DG, Duncan MA: Guide to the identification and geographic distribution of Lutzomyia sand flies in Mexico, the West Indies, Central and South America (Diptera: Psychodidae). 1994, Gainesville, FL: Associated Publishers

  55. 55.

    Chaniotis BN, Tesh RB, Correa MA, Johnson KM: Diurnal resting sites of phlebotomine sandflies in a Panamanian tropical forest. J Med Entomol. 1972, 9: 91-98.

  56. 56.

    Venables WN, Ripley BD: Modern applied statistics with S. 2002, New York: Springer

  57. 57.

    Shumway TH, Stoffer DS: Time series analysis and its applications. 2000, New York: Springer

  58. 58.

    Christensen HA, de Vasquez AM, Petersen JL: Short report epidemiologic studies on cutaneous leishmaniasis in eastern Panama. Am J Trop Med Hyg. 1999, 60: 54-57.

  59. 59.

    Kirby MJ, Milligan PJ, Conway DJ, Lindsay SW: Study protocol for a three-armed randomized controlled trial to assess whether house screening can reduce exposure to malaria vectors and reduce malaria transmission in The Gambia. Trials. 2008, 9: 33-10.1186/1745-6215-9-33.

  60. 60.

    Rutledge LC, Mosser HL: Biology of immature sandflies (Diptera: Psychodidae) at the bases of trees in Panama. Environ Entomol. 1972, 1: 300-309.

  61. 61.

    Rutledge LC, Ellenwood DA: Production of phlebotomine sandflies on the open forest floor in Panama: the species complement. Environ Entomol. 1975, 4: 71-77.

  62. 62.

    Rutledge LC, Ellenwood DA: Production of phlebotomine sandflies on the open forest floor in Panama: phytologic and edaphic relations. Environ Entomol. 1975, 4: 83-89.

  63. 63.

    Rutledge LC, Ellenwood DA: Production of phlebotomine sandflies on the open forest floor in Panama: hydrologic and physiographic relations. Environ Entomol. 1975, 4: 78-82.

  64. 64.

    Feliciangeli MD: Natural breeding places of phlebotomine sandflies. Med Vet Entomol. 2004, 18: 71-80. 10.1111/j.0269-283X.2004.0487.x.

  65. 65.

    Ferro C, Pardo R, Torres M, Morrison AC: Larval microhabitats of Lutzomyia longipalpis (Diptera: Psychodidae) in an endemic focus of visceral leishmaniasis in Colombia. J Med Entomol. 1997, 34: 719-728.

  66. 66.

    Chaves LF, Morrison AC, Kitron UD, Scott TW: Nonlinear impacts of climatic variability on the density-dependent regulation of an insect vector of disease. Glob Chang Biol. 2012, 18: 457-468. 10.1111/j.1365-2486.2011.02522.x.

  67. 67.

    Chaves LF, Añez N: Species co-occurrence and feeding behavior in sand fly transmission of American cutaneous leishmaniasis in western Venezuela. Acta Trop. 2004, 92: 219-224. 10.1016/j.actatropica.2004.08.001.

  68. 68.

    Chaves LF: Sand fly species co-occurrence at the local scale: differences between agricultural and forested areas. Bol Malariol Sal Amb. 2011, 51: 35-39.

  69. 69.

    Chaves LF, Hamer GL, Walker ED, Brown WM, Ruiz MO, Kitron UD: Climatic variability and landscape heterogeneity impact urban mosquito diversity and vector abundance and infection. Ecosphere. 2011, 2: art70-10.1890/ES11-00088.1.

  70. 70.

    Bolker BM: Ecological Models and Data in R. 2008, Princeton, NJ: Princeton University Press

  71. 71.

    Krebs CJ: Ecological Methodology. 1998, San Francisco: Benjamin Cummings

  72. 72.

    Faraway JJ: Linear Models with R. 2004, Boca Raton: CRC Press

Download references

Funding

This work received financial support from SENACYT grant CCP06-040, Netropica grant 04-N-2010, and Nagasaki University (Program for Nurturing Global Leaders in Tropical and Emerging Communicable Diseases).

Acknowledgments

We acknowledge with appreciation the field and technical assistance of Roberto Rojas and Jose Montenegro, as well as, all members of the Parasitology department at the Gorgas Institute. We thank Prof. Mario Baldi for his support at Universidad Nacional de Costa Rica. We also thank Prof. H. Kawada and Prof. N. Minakawa for comments on insecticide application and vector control and Miss Junko Sakemoto for administrative support at Nagasaki University.

Author information

Correspondence to Luis Fernando Chaves.

Additional information

Competing interests

All authors declare that they have no competing interests.

Authors’ contributions

Conceived and designed the experiments: JEC, AZ, LFC. Performed the experiments: AZ, CR, AV, NG, JEC. Analyzed the data: LFC, JEC. Contributed reagents/materials/analysis tools: LFC, AZ, CR, AV, NG, JEC. Wrote the paper: LFC, JEC, AZ. All authors read and approved the final version of the manuscript.

Electronic supplementary material

Additional file 1: Protocol S1. Power estimation [19, 53, 70] (PDF 45 kb) (PDF 45 KB)

Additional file 2: Figure S1: Power Analysis (A) Assuming a 20% reduction on Sand Fly Abundance (B) Assuming a 50% reduction on Sand Fly Abundance. ITF = Insecticide Thermal Fogging. For further details see Protocol S1. (TIFF 4 MB)

Additional file 3: Protocol S2. Supplementary Methods [12, 19, 40, 42, 55, 56, 6063, 7072]. (PDF 37 KB)

Additional file 4: Table S1: Selection of the best negative binomial model explaining the abundance of all Phlebotomine Sand Flies in Trinidad de Las Minas, Capira, Panamá, following two foggings with deltamethrin [6 mg a.i.m-2]. (PDF 60 KB)

Additional file 5: Table S2: Model selection for the best negative binomial models explaining the abundance of the five most abundant sand fly species in Trinidad de Las Minas, Capira District, Panamá, following two insecticide foggings with deltamethrin [6 mg active ingredient /m2]. (PDF 75 KB)

Additional file 6: Table S3: Principal components analysis used to estimate the housing destituteness index. (PDF 52 KB)

Additional file 7: Figure S2: Peridomicile index and sand fly abundance (A) Cumulative sand fly abundance before the fogging as a function of the peridomicile index, PI (three nights) (B) Cumulative sand fly abundance for the nine nights that contained the two interventions as a function of PI. The solid black line is the fit from a negative binomial model (C) Changes in domiciliary sand fly density per house and trap night as a function of PI (D) Changes in peridomiciliary sand fly density per house and trap night as function of PI. In panels (A) to (D) symbols indicate whether a house was intervened, triangles, or not, circles. In panels (C) and (D) black symbols represent the pre-intervention densities and green symbols post-intervention densities. To ease the tracking of changes in each house we joined the pre and post intervention densities with a line whose color is blue when sand fly density decreased and red when sand fly density increased. In panel B the dark line is the fit of an NB-GLM, where the intercept (± S.E.) is 4.94 ± 0.18 and the slope (± S.E.) is 0.67 ± 0.25. This means that for a PI of 0 there where ~140 sand flies and this number was doubled by each unit increase in PI. (TIFF 1986 kb) (TIFF 2 MB)

Additional file 8: Table S4: Principal components analysis used to estimate the peridomicile index. (PDF 53 KB)

Additional file 9: Figure S3:Vegetation index and sand fly abundance (A) Cumulative sand fly abundance before the fogging as a function of Vegetation Index, VI (three nights) (B) Cumulative sand fly abundance for the nine nights that contained the two interventions as a function of VI. The solid black line is the fit from a negative binomial model (C) Changes in domiciliary sand fly density per house and trap night as function of VI (D) Changes in peridomiciliary sand fly density per house and trap night as a function of VI. In panels (A) to (D) symbols indicate whether a house was intervened, triangles, or not, circles. In panels (C) and (D) black symbols represent the pre-intervention densities and green symbols post-intervention densities. To ease the tracking of changes in each house we joined the pre and post intervention densities with a blue line when sand fly density decreased and red when sand fly density increased. (TIFF 2 MB)

Additional file 10: Table S5: Principal components analysis used to estimate the vegetation structure index. (PDF 50 KB)

Additional file 11: Figure S4: Host abundance/richness and sand fly abundance. Cumulative sand fly abundance before the fogging (three nights) as a function of (A) Wild animal index (B) Domestic animal index (C) Wild and domestic animal index (D) Wild animal species richness (E) Domestic animal species richness (F) Household human density. Cumulative sand fly abundance after the foggings (nine nights) as a function of (G) Wild animal index (H) Domestic animal index (I) Wild and domestic animal index (J) Wild animal species richness (K) Domestic animal species richness (L) Household human density. (TIFF 3 MB)

Additional file 12: Table S6: Principal components analysis used to estimate the animal abundance indices. (PDF 55 KB)

Additional file 13: Table S7: Model selection for the best negative binomial model explaining post-fogging sand fly abundance in the houses. (PDF 58 KB)

Authors’ original submitted files for images

Below are the links to the authors’ original submitted files for images.

Authors’ original file for figure 1

Authors’ original file for figure 2

Authors’ original file for figure 3

Authors’ original file for figure 4

Rights and permissions

This article is published under license to BioMed Central Ltd. This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/2.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

Reprints and Permissions

About this article

Keywords

  • Panama
  • Lutzomyia gomezi
  • Lu. trapidoi
  • Deltamethrin
  • Housing quality

Comments

By submitting a comment you agree to abide by our Terms and Community Guidelines. If you find something abusive or that does not comply with our terms or guidelines please flag it as inappropriate. Please note that comments may be removed without notice if they are flagged by another user or do not comply with our community guidelines.