Impact of Indoor Residual Spraying with Bendiocarb (FICAM WP 80) on Entomological Indicators of Malaria Transmission in Senegal ()
1. Introduction
Malaria remains a major public health problem and one of the leading causes of mortality and morbidity worldwide. According to the latest World Malaria Report by the World Health Organization (WHO), there were an estimated 263 million cases and 597,000 malaria deaths worldwide in 2023 [1]. This represents about 11 million more cases in 2023 compared to 2022, and nearly the same number of deaths. Approximately 95% of the deaths occurred in the WHO African Region [1].
In Senegal, malaria remains a major public health problem despite the considerable progress achieved through the implementation and scale-up of effective control strategies. These efforts include rapid diagnosis and effective case management, as well as preventive interventions such as universal coverage with long-lasting insecticidal nets (LLINs) and targeted indoor residual spraying (IRS). Between 2009 and 2016, a substantial decline in parasite prevalence was observed, dropping from 3% to 1% [2]. However, this overall decline masks disparities, with malaria transmission remaining heterogeneously distributed across the country.
As in several other sub-Saharan African countries, IRS has played an important role in malaria prevention in Senegal, supported financially by the from the U.S. President’s Malaria Initiative (PMI) [3], in line with the WHO Roll Back Malaria Partnership’s Global Malaria Action Plan [4]. Historically, IRS has been one of the main vector control tools and a cornerstone of the WHO Global Malaria Eradication Program launched in the 1950s, contributing to the elimination or substantial reduction of malaria in parts regions of Europe, Asia, and Latin America [5]. In Africa, both historical and recent studies have demonstrated the effectiveness of IRS, particularly in southern Africa, with significant reductions in malaria transmission [6]-[11]. However, many of these studies lacked appropriate control groups, limiting the accuracy of the estimated impact of IRS. A recent Cochrane review applying strict inclusion criteria identified only four eligible studies in Africa, three from areas of stable transmission, including two randomized controlled trials [5].
In Senegal, IRS was first introduced in 2007 in selected districts [12]. Pyrethroids were initially used until 2010 [13], and then due to the emergence and spread of pyrethroid resistance [14] [15], carbamates (bendiocarb, FICAM WP 80) were adopted from 2011 to 2014 in six health districts: Nioro du Rip, Guinguineo, Malem Hodar, Koungheul, Koumpentoum, and Vélingara.
This study retrospectively evaluates the impact of IRS with bendiocarb (FICAM WP 80) on malaria transmission in intervention areas in Senegal.
2. Materials and Methods
2.1. Study Area
The study was conducted in six health districts in Senegal where IRS with bendiocarb (FICAM WP 80) was implemented between 2011 and 2014. District selection was carried out in collaboration with the medical officers of each district hospital, based on the following criteria: i) malaria epidemiology sentinel villages; ii) geographical characteristics of the district; iii) presence of a health facility capable of supporting parasitological and clinical studies; and, iv) availability of previous entomological and/or parasitological data.
For analysis, three of the six treated districts were selected and compared with untreated control districts (Figure 1). Specifically, Malem Hodar and Koungheul were compared with the untreated control district of Kaffrine, located in central Senegal (Sudanian zone), while Vélingara was compared with the untreated control district of Kolda, situated in the south (Sudano-Guinean savanna zone). These districts reported an annual malaria incidence exceeding 25‰.
Figure 1. Geographic location of treated and untreated districts.
A total of 15 treated villages (five per treated district) and six untreated villages (three per control district) were included in entomological monitoring to assess the bionomics of malaria vector populations in relation to malaria transmission. The surveys started one month after IRS operations, which were conducted between June and August depending on the district.
2.2. Mosquito Sampling
2.2.1. Collection of Endophilic Mosquitoes by PSC
In each selected village, 10 rooms were randomly selected as collection points for endophilic mosquitoes. Resting mosquitoes were sampled using pyrethrum spray catches (PSC) from a total of 300 rooms in treated sites (50 per district) and 60 houses in untreated control districts (30 per district) throughout the study period.
2.2.2. Collection of Host-Seeking Females by HLC
In both treated and control districts, three villages were selected for indoor and outdoor nighttime collection of host-seeking female mosquitoes using the human landing catch (HLC) method. Collections were performed by adult volunteers during two consecutive nights from 9:00 p.m. to 6:00 a.m. [16].
Entomological surveys were conducted in parallel with post-spray monitoring of residual IRS efficacy with bendiocarb [15]. Pairwise comparisons of transmission parameters between the treated districts of Malem Hodar and Koungheul and their untreated control, Kaffrine were carried out during the 2012, 2013, and 2014 IRS campaigns. In contrast, for Velingara and its untreated control, Kolda, comparisons were performed only in 2012 and 2013, as IRS was discontinued in Velingara after the 2013 campaign.
2.3. Field Processing
Collected mosquitoes were morphologically identified to genus and species using a standard identification key [17], then sorted, counted, and classified by sex. The abdomens of blood-fed females collected by PSC were crushed onto Whatman filter papers, which were stored with silica gel for subsequent identification of blood meal sources in the laboratory. approximately 40-50% of host-seeking females collected by HLC were dissected to determine parity status [18].
2.4. Laboratory Processing
Detection of Plasmodium falciparum-infected females was performed using the enzyme-linked immunosorbent assay circumsporozoite (ELISA CSP) method [19] [20]. Blood meal origins were determined using the indirect ELISA method [21].
2.5. Data Analysis
All field and laboratory data were recorded in an Excel database. Malaria transmission parameters (biting rate per human per night [B/H/N], indoor resting density [IRD], parity rate, endo/exophagy rate, anthropophilic rate, and infection rate) were estimated for each district [16], then compared between treated districts and their respective controls.
Mean resting and biting densities of An. gambiae s.l. were compared between sites and treatment status using analysis of variance (ANOVA). All statistical analyses were conducted with R software version 3.3.1 [22], applying a significance threshold of 5%.
3. Results
3.1. Anopheles Fauna and Relative Abundance
Between 2012 and 2014), a total of 6,224 Anopheles mosquitoes were collected using HLC and PSC in the five study districts. Six species were identified, with An. gambiae s.l. being the predominant vector, accounting for 97.7% of the total catch. The remaining species included An. funestus (1%), An. pharoensis (0.8%), An. rufipes (0.4%), and An. nili (0.1%) (Table 1).
Table 1. Anopheline species composition and relative abundance collected by HLC and PSC in the study districts (2012-2014).
Districts |
An. gambiae |
An. funestus |
An. pharoensis |
An. rufipes |
An. nili |
|
n |
% |
n |
% |
n |
% |
n |
% |
n |
% |
Kaffrine |
685 |
98 |
7 |
1 |
4 |
0.6 |
3 |
0.4 |
0 |
0 |
Koungheul |
400 |
98 |
0 |
0 |
1 |
0.2 |
7 |
1.8 |
0 |
0 |
Malem Hodar |
914 |
99 |
3 |
0.3 |
7 |
0.7 |
0 |
0 |
0 |
0 |
Kolda |
2297 |
96.2 |
35 |
1.4 |
36 |
1.5 |
16 |
0.7 |
2 |
0.2 |
Velingara |
1785 |
98.8 |
18 |
0.9 |
3 |
0.2 |
0 |
0 |
1 |
0 |
3.2. Indoor Resting Densities (IRD) of An. gambiae s.l. in Treated
and Control Districts
Given its predominance in all districts (97.7% of collected Anopheles), the analysis of entomological parameters focused on An. gambiae s.l. Peak mean indoor resting densities (IRD, expressed as females per room, F/R) were 3.5 F/R; 2 F/R, and 1.7 F/R in Malem Hodar, and 3.1 F/R; 4.2 F/R, and 0.6 F/R in Koungheul, and 5.5 F/R; 6.1 F/R, and 2.9 F/R in Kaffrine during 2012, 2013, and 2014 campaigns, respectively (Figure 2 and Figure 3). In Vélingara, peak IRDs reached 3.1 F/R in 2012 and 4.4 F/R in 2013, compared to 6.2 F/R and 7.4 F/R in the control district of Kolda (Figure 4).
During the 2012 season, IRDs were significantly higher in Kolda (control) compared with Velingara (treated) during the second month of follow-up (p = 0.036) (Figure 4), and in Kaffrine compared with Malem Hodar at the end of follow-up in November (p = 0.032) (Figure 3). In 2013, IRDs were significantly higher in Kaffrine compared with Koungheul in September (p = 0.032) (Figure 2), and at the beginning of follow-up (p = 0.032) in Kolda compared with Velingara (p = 0.032) (Figure 4). In 2014, IRDs were significantly higher in Kaffrine compared with Koungheul both at the beginning (August, p = 0.016) and at the end of follow-up (November, p = 0.017) (Figure 2).
In the treated districts (Koungheul, Malem Hodar and Vélingara), the highest IRDs were generally observed during the third month post-treatment (Figures 2-4). In contrast, in the control districts of Kaffrine and Kolda, peak IRDs occurred mainly in August and September, coinciding with the wettest months of the year (Figures 2-4).
Figure 2. Variations in mean indoor resting densities of An. gambiae s.l. between Koungheul (treated) and Kaffrine (control) during the 2012, 2013 and 2014 campaigns.
Figure 3. Variations in mean resting densities of An. gambiae s.l. between Malem Hodar (treated) and Kaffrine (control) during three campaigns 2012, 2013 and 2014 campaigns.
Figure 4. Variations in mean resting densities of An. gambiae s.l. between Vélingara (treated) and Kolda (control) during the 2012 and 2013 campaigns.
3.3. Biting Densities between Treated and Untreated Control
During the study campaigns, the highest mean biting rates of An. gambiae s.l females were recorded as follows: 1.25 B/H/N; 0.83 B/H/N and 0.25 B/H/N in Malem Hodar, 3.1 B/H/N; 4.2 B/H/N and 0.6 B/H/N in Koungheul, and 2.58 B/H/N; 5.25 B/H/N and 2 B/H/N in Kaffrine, respectively, in 2012, 2013, and 2014 (Figure 5 and Figure 6).
During the 2012 and 2013 seasons, the highest mean biting rates were observed in Vélingara (30.5 B/H/N in 2012 and 13.5 B/H/N in 2013) and in its control district Kolda (14.5 B/H/N in 2012 and 20.83 B/H/N in 2013) (Figure 7).
Biting densities were significantly higher in Kaffrine compared to Koungheul at the fifth month of follow-up in 2012 (p = 0.030) and at the beginning of the monitoring period in July 2013 (p = 0.026) (Figure 5). In the treated districts, peaks in biting densities were mainly observed from the third month after treatment. In contrast, in the control districts, the highest peaks were generally observed in August and September, corresponding to the wettest months of the year.
3.4. Endophagy and Exophagy Rates between Control and
Treated Sites
Endophagy and exophagy rates were assessed in the treated districts of Koungheul and Malem Hodar, as well as in their untreated control, Kaffrine, over three consecutive years 2012, 2013, and 2014 (Figure 8 and Figure 9). Similar monitoring was conducted in the treated district of Vélingara and its untreated control, Kolda, during two years (2012 and 2013) (Figure 10). Although reductions in endophagy rates of An. gambiae s.l. females were observed in the treated sites, these differences were not statistically significant (p > 0.05) throughout the follow-up.
Figure 5. Variations of An.gambiae s.l. biting rates between Kaffrine (control) and Koungheul (treated) during the 2012, 2013 and 2014 according to the monitoring months.
Figure 6. Variations of An.gambiae s.l. biting rates between Kaffrine (control) and Malem Hodar (treated) during the 2012, 2013 and 2014 campaigns according on the months of monitoring.
Figure 7. Variation of An. gambiae s.l. biting rates in Kolda (control) and Vélingara (treated) during the 2012 and 2013 campaigns according to the months of monitoring.
Figure 8. Comparison of exophagy and endophagy rates of An. gambiae s.l. between Koungheul (treated) and Kaffrine (control). The bars represent the 95% confidence interval.
Figure 9. Comparison of exophagy and endophagy rates between Malem Hodar (treated) and Kaffrine (control). The bars represent the 95% confidence interval.
Figure 10. Comparison of exophagy and endophagy rates between Velingara (treated) and Kolda (control). The bars represent the 95% confidence interval.
3.5. Parity Rates of An. gambiae between Control and Treated Sites
Parity rates were assessed in the treated districts of Koungheul and Malem Hodar, along with their untreated control in Kaffrine, over three years (2012, 2013, and 2014) (Figure 11 and Figure 12). The same monitoring was conducted in the treated district of Vélingara and its untreated control, Kolda, for two years (2012 and 2013) (Figure 13).
During the follow-up, a significant decrease in parity rates was observed in Vélingara at the second (χ2= 3.91, ddl = 1, p = 0.05) and third (χ2= 4.05, df = 1, p = 0.04) monitoring points in 2012 (Figure 13). In 2013 and 2014, although lower parity rates were generally observed in the treated districts, the only significant difference was recorded in Kaffrine at the second month post-treatment (χ2= 14.82, ddl = 1, p = 0.0001) in 2013 in Kaffrine (Figure 11 and Figure 12).
3.6. Anthropophilic Rate between Control and Treated Sites
The comparative study of trophic preferences was conducted during the 2012 campaign in treated districts (Koungheul, Vélingara and Malem Hodar) and in untreated controls (Kaffrine and Kolda) (Figure 14). The origin of the blood meals of 425 engorged female mosquitoes collected inside human dwellings was determined using direct ELISA. The blood meals were tested against five vertebrate hosts (human, bovine, ovine, chicken and equine). No blood meal was detected from chicken, and the percentage of unidentified meals was estimated at 0.23%.
Figure 11. Comparison of parity rates of An. gambiae s.l. between Koungheul (treated) and Kaffrine (control). Error bars represent the 95% confidence intervals.
Figure 12. Comparison of parity rates of An. gambiae s.l. between Malem Hodar (treated) and Kaffrine (control). Error bars represent the 95% confidence intervals.
Figure 13. Comparison of parity rates of An. gambiae s.l. in Vélingara (treated) and Kolda (control). Error bars represent the 95% confidence intervals.
Figure 14. Origin of blood meals of An. gambiae s.l. females collected at rest inside human dwellings in treated districts (Koungheur, Vélingara, Malem Hodar) and untreated controls (Kaffrine, Kolda) during the 2012 campaign.
Across the study sites, most blood meals were taken either from horses (notably in Malem Hoddar) or from humans (in Kaffrine, Kolda and Vélingara) (Figure 14). A small proportion of specimens (8.5%) fed on at least two hosts. No significant difference was found in the anthropophilic rate between Kaffrine and the treated districts of Malem Hodar (χ2 = 0.19, df = 1, p = 0.66) and Koungheul (χ2 = 0.02, df = 1, p = 0.90). Similarly, no significant difference was obsezrved between Kolda and Vélingara, (χ2= 0.009, df = 1, p = 0.92) (Figure 14).
3.7. Infection Rates between Treated and Control Districts
The comparison of infection rates was carried out during the 2012 campaign (Table 2). A total of 2571 An. gambiae s.l. females collected by HLCs were tested by ELISA for the presence of P. falciparum circumsporozoite protein (CSP). Except for the Kolda and Vélingara sites, no infected females were detected in the other districts (Table 2). CSP antigen was detected in 56 An. gambiae s.l. specimens, corresponding to an infection rate of 1.9% in Kolda (13/691) in Kolda and 2.6% (43/1636) in Vélingara. No significant difference was observed between the infection rates between in these two sites (χ2= 0.86, ddl = 1, p = 0.35) (Table 2).
Table 2. Infection rates of biting An. gambiae s.l. females tested for Plasmodium Falciparum circumsporozoite antigen during the 2012 campaign.
Sites |
Tested |
Positive |
SI (%) |
Kaffrine |
132 |
0 |
0 |
Kolda |
691 |
13 |
1.9 |
M. Hodar |
40 |
0 |
0 |
Koungheul |
72 |
0 |
0 |
Vélingara |
1636 |
43 |
2.6 |
Sporozoite Index (SI).
4. Discussion
The objective of this study was to monitor the post-treatment evolution of key entomological parameters, including species composition, abundance, parity rates, plasmodial infection rates, and biting and resting densities of female vectors, particularly An. gambiae s.l.
Post-treatment monitoring of species composition and abundance revealed the presence of An. gambiae s.l., An. funestus, An. pharoensis, An. rufipes, and An. nili, with a clear predominance of An. gambiae s.l across all treated and control districts. Such predominance has already been reported throughout sub-Saharan Africa [23] [24], in several ecological zones of Senegal: coastal [25], Sahelian [26], Senegal River basin [27] [28], Sahelo-Sudanian [29] [30], Sudano-Guinean zone [31], where An. gambiae members are considered major malaria vectors. This predominance can be explained by the availability of their preferred larval habitats, namely temporary, shallow and sunlit water collections linked to rainfall [17].
The low numbers of An. funestus may be related to the scarcity of its preferential sites (permanent habitats such as swamps, lakes and pond margins, grassy riverbanks and flooded zones) [32]. Similarly, the rarety of An. pharoensis in our study sites is consistant with its known distribution in Senegal, where it is more abundant in the Senegal River delta [25]. The very limited presence of An. nili could be explained by the absence of its preferred well-oxygenated running waters [33]. Finally, the very low numbers of An. rufipes observed may reflect its zoophilic behavior, as reported in Dielmo [34], although this species has been described as a potential vector in other African regions [17].
The low biting and resting densities of An. gambiae s.l. observed in treated districts could be explained by the stressful environment created large-scale spraying, which deters mosquitoes from entering sprayed houses [35]-[37]. Nevertheless, a fraction proportion of females still manage to enter rooms, obtain a blood meal, and rest briefly before escaping. The gradual increase in densities observed from the third month post-treatment likely reflects the progressive loss of residual activity of bendiocarb on treated walls [38].
In untreated districts, the high vector densities recorded August, September, and October coincide with peak rainfall, confirming the strong seasonality of Anopheles populations, as also reported in Ghana in 2008 [39] and Cameroon in 2015 [40].
The reduced resting densities in treated houses are consistent with the deterrent effect of IRS. However, slight increase in endophagic rates of An. gambiae s.l. observed from the first to the third month post-treatment in Vélingara, without significant differences compared to Kolda (control), may also reflect the decline in bendiocarb efficacy. Similar results have been reported in Bénin in 2011 [35], where IRS with bendiocarb created an unfavorable indoor environment, leading to increased exophilic. Other studies in Equatorial Guinea in 2011 [41] and Tanzania in 2011 [42] have also shown that IRS can shift mosquito behavior towards exophily.
This study also revealed decreased parity rates in treated districts. Such reductions have been documented in other IRS programs [43], reflecting decreased mosquito survival. However, the subsequent increase in parity rates from the fourth month onward may reflect the waning efficacy of bendiocarb.
Regarding blood meal sources, An. gambiae s.l. females fed on different vertebrate hosts but preferentially on humans, horses or cattle, depending on the site. The decrease in strictly human blood meals in some treated districts could be explained by both the impact of spraying on host-seeking behavior and the availability of alternative animal hosts. In Velingara, however, the combination of exophily and high biting rates resulted in higher anthropophilic rates.
The high infection rates recorded in Vélingara, but not in other treated districts, are consistent with these anthropophily findinds. Such results corroborate previous hut trials demonstrating the efficacy of carbamates against An. gambiae s.l. [34] [44] [45].
Although molecular identification of An. gambiae complex members was performed, entomological parameter comparisons in this study were conducted at the complex level, thus, these results were presented separately. Infection rates could only be compared for one district pair ans one campaign year as these were the only data collected in parallel.
5. Conclusions
The study demonstrates that indoor residual spraying with bendiocarb significantly impacted entomological parameters in treated areas of Senegal. In particular, biting rates decreased markedly, resting densities were reduced, and parity rates declined. These effects likely contributed to the strong reduction in infection rates observed. However, results from Vélingara highlight that An. gambiae s.l. populations in this area exhibit marked exophily, which may limit the impact of IRS. Despite the short residual activity of bendiocarb, our findings confirm that it remains a promising tool for controlling pyrethroid-resistant malaria vectors. Nevertheless, it should not be deployed as a standalone strategy.
Thus, based on these results, the following recommendations can be made:
1) implement at least two rounds of IRS to cover the full malaria transmission season, given the short persistence of bendiocarb;
2) combine IRS with long-lasting insecticide-treated nets to maximize protection and reduce malaria transmission in Senegal.
Authors’ Contributions
Cheikh Lo, Ibrahima Dia, and Omar Faye designed and supervised the study. Cheikh Lo, Assane Ndiaye, El Hadji Amadou NIANG, and Ibrahima Dia performed field and laboratory activities. Cheikh Lo and Ibrahima Dia carried out the data analysis. All authors drafted, revised, and approved the final version of the manuscript.
Acknowledgements
We are grateful to the inhabitants of the studied villages for their cooperation and active participation in this study.
Abbreviations
B/H/N |
Bites per human per night |
ELISA CSP |
Enzyme-Linked Immunosorbent Assay circumsporozoite |
F/R |
females per room |
HLC |
Human landing catch method |
IRS |
Indoor residual spraying |
IRD |
indoor resting densities |
SI |
Sporozoite Index |
LLINs |
Long-lasting insecticidal nets |
PSC |
Pyrethrum Spray Catch |
WP |
wet powder |
WHO |
World Health Organization |