Anopheles mosquito species composition and abundance
During the study period, a total of 11237 An. gambiae s.l. females were collected by PSC and HLC approaches. Among these specimens, 46.6% (5239/11237) were collected in Dangassa and 53.4% (5998/11237) in Koïla Bamanan. Based on logistics and practical considerations, a sub-sample was randomly selected across surveys to conduct molecular species identification. Of the 679 identified by PCR in Dangassa, An. coluzzii was the predominant species with a prevalence of 91.4% followed by An. gambiae (8.0%) and An. arabiensis (< 1%). In Koïla Bamanan, 623 An. gambiae s.l. were identified by the same methods and An. coluzzii represented 99%, An. arabiensis 0.8% and An. gambiae 0.2% of the catches. Overall, the An. coluzzii frequency was higher in Koïla Bamanan (99%) than in Dangassa (91.4%) (χ2 = 38.99, P < .001). No significant difference was observed between the frequencies of An. arabiensis (χ2 = 0.001, P = 0.970) and An. gambiae (χ2 = 0.081, P = 0.770) at both sites.
Density And Human Blood Index
In terms of abundance, Fig. 2 shows the monthly mean densities of An. gambiae s.l. in both localities using the PSC method. The mean density over the entire study period was 12.8 An. gambiae s.l. per room/per day [95% BCI: 12.5—13.2] in Koïla Bamanan and was significantly higher in Dangassa (6.4 per room/day [95% BCI : 6.2—6.7]). The highest densities were observed in April and August of 2014 for both villages. However, the densities of Koïla Bamanan were considerably higher than those of Dangassa. The lowest densities were observed in April 2016 for both study sites.
Over the 12 surveys, the mean HBIs of An. gambiae s.l. for Dangassa and Koïla Bamanan were 79.4% and 15.9%, respectively (Fig. 2). In Dangassa, the highest HBI was observed in August 2014 (96.8%) and the lowest HBI (52.1%) in April 2015, while in Koïla Bamanan, the lowest (6.7%) was observed in October 2012 and the highest (46.0%) in April 2015. HBIs were significantly higher in Dangassa than in Koïla Bamanan (χ2 = 1913.6, P < 0.001) during the entire study period.
Monthly human biting rates (MHBRs), sporozoite infection rates (SIR) and entomological inoculation rates (EIRs)
Figure 3 shows the human biting rates and the entomological inoculation rates from PSC in Dangassa and Koïla Bamanan respectively. The average monthly human biting rate (MHBR) was 2.3 times higher in Koïla Bamanan (69.8 bites/person/month) than in Dangassa (30 b/p/m). Monthly variations in MHBRs were observed at both sites. The highest MHBRs were recorded in April and August 2014. However, they were higher in Koïla Bamanan (266.6 and 249.6 b/p/m) than in Dangassa (65.9 and 105.5 b/p/m) (Fig. 3).
There was no significant difference in overall (over the entire study period) sporozoite infection rates (SIR) between the two study sites (infection rates were 2.9%, [95% BCI: 2.3—3.6 for Dangassa versus 2.6% [95% BCI: 2.2—3.1] for Koïla Bamanan; Table 2). Monthly and annual variations were observed at both sites. In Dangassa, the highest SIRs were observed in October (8.9%) in both 2012 and October 2014 (5.9%). While the peak of the infection rates were observed in October in Dangassa, in Koïla Bamanan there was a second peak that was observed for April though the two peaks were seen for April 2014 (7.1%) and 2016 (78.5%).
Table 2
Sporozoite rates of Anopheles gambiae s.l. mosquitoes collected by spray-catch from Dangassa and Koïla Bamanan form july 2012 to june 2016.
Survey
|
Dangassa
|
Koïla Bamanan
|
#tested
|
IR (%)
|
BCI
|
#tested
|
IR (%)
|
BCI
|
Jul 2012
|
226
|
2.2
|
[0.7—4.5]
|
656
|
0.6
|
[0.2—1.3]
|
Oct 2012
|
181
|
8.9
|
[5.2—13.4]
|
351
|
1.7
|
[0.6—3.3]
|
Nov 2012
|
22 4
|
4.5
|
[2.2—7.5]
|
96
|
0.1
|
[0.0—1.1]
|
Apr 2013
|
81
|
0.1
|
[0.0—1.3]
|
269
|
0.0
|
[0.0—0.4]
|
Apr 2014
|
171
|
2.9
|
[1.0—5.9]
|
875
|
7.1
|
[5.5—8.9]
|
Aug 2014
|
596
|
1.5
|
[0.7—2.6]
|
1078
|
1.4
|
[0.8—2.2]
|
Oct 2014
|
169
|
5.9
|
[2.9—9.9]
|
16
|
6.4
|
[0.2—22.1]
|
Apr 2015
|
337
|
1.8
|
[0.7—3.5]
|
183
|
2.8
|
[0.9—5.6]
|
Jun 2015
|
141
|
2.9
|
[0.8—6.2]
|
572
|
3.9
|
[2.4—5.6]
|
Nov 2015
|
121
|
0.1
|
[0.0—0.9]
|
47
|
2.2
|
[0.1—7.9]
|
Apr 2016
|
46
|
2.3
|
[0.1—8.0]
|
14
|
78.5
|
[54.5—94.9]
|
Jun 2016
|
369
|
3.3
|
[1.7—5.3]
|
963
|
0.7
|
[0.3—1.4]
|
TOTAL
|
2662
|
2.9
|
[2.3—3.6]
|
5120
|
2.6
|
[2.2—3.1]
|
Given the low human blood index observed in Koïla Bamanan, the EIR was corrected by multiplying it by the human blood index observed at each site. The mean EIR was significantly higher in Dangassa (0.70 infective bites/person/month = ib/p/m) than in Koïla Bamanan (0.29 ib/p/m; Fig. 3). In Dangassa, the highest monthly EIRs were observed during the months of October 2012 and 2014, (2.75 and 1.81 ib/p/m). In Koïla Bamanan, the highest EIRs were recorded for April 2014 and June 2015 (3.33 and 0.59 ib/p/m, respectively). As with the sporozoite infection rates, EIRs also were subject to monthly and annual variations.
Outdoor vs indoor biting activities of An. gambiae s.l. in Dangassa and Koïla Bamanan
Data for comparing the location of biting activities (indoor versus outdoor locations) were generated using mosquitoes collected by human landing catches. HBRs showed a wide patterns of monthly variation for both indoor and outdoor locations in both villages (Fig. 4).
Figure 4 shows HBR per hour and by season. HBR was higher during the second part of the night indoors and outdoors regardless of site and season. We observed a peak of HBR outdoors between 02:00 h and 03:00 h, and then two peaks respectively between 01:00 h − 02:00 h and 04:00 h -05:00 h during the rainy season at Dangassa (Fig. 4A). In the dry season, the peak was observed between 02:00 h and 04:00 h, but it was lower than in the wet season (Fig. 4B). At Koïla Bamanan, two peaks were observed indoor during the rainy season, where the first occurred between 23:00 h and 00:00 h, and the second occurred between 01:00 h and 02:00 h (Fig. 4C). In the dry season, a single peak was observed between 23:00 h and 00:00 h indoor and outdoor at Koïla Bamanan (Fig. 4D).
In Dangassa, the highest MHBRs were observed in July, October and November 2012, and in April and August 2014 (Fig. 5). The average HBR over the study period was significantly higher outdoors (648.9 b/p/m [BCI: 642.6—655.2]) than indoors (560.3 b/p/m [BCI: 554.7—565.9]) suggesting an exophagic coefficient of 1.2 (648.9/560.3).
In Koïla Bamanan, the highest MHBRs were observed in October 2012, April and August 2014 and in June 2016 indoor and outdoor (Fig. 5). As observed in Dangassa, in Koïla Bamanan the average HBR over the study period was significantly higher for outdoor locations (242.4 b/p/m [BCI: 238.1—246.7]) than indoor locations (228.9 b/p/m [BCI: 225.1—232.6]) with an exophagic coefficient of 1.1 (242.4 / 228.9). HBRs were significantly higher in Dangassa than in Koïla Bamanan, both indoors and outdoors.
In Dangassa, the sporozoite infection rate in An. gambiae s.l. was 3.1% [BCI: 2.1—4.3] and 3.6% [BCI: 2.4—5.0] for indoor and outdoor locations, respectively. No significant differences were observed between these two rates. The sporozoite infection rates were subject to wide month-to-month variability. In Koïla Bamanan, the average sporozoite infection rate was 4.3% [95% BCI: 2.7—6.3] and 2.4% [95% BCI: 1.1—4.2] for indoor and outdoor locations, respectively. The sporozoite infection rate indoors was 1.8 times higher than outdoors, though no significant difference was observed between the two. Similar monthly variations were observed both indoors and outdoors for Dangassa (Table 3),
Table 3
Sporozoite rates of Anopheles gambiae s.l. mosquitoes collected by human landing catch indoors and outdoors in Dangassa and Koïla Bamanan from July 2012 to June 2016.
Surveys
|
Dangassa
|
Koïla Bamanan
|
Indoor
|
Outdoor
|
Indoor
|
Outdoor
|
#tested
|
IR (%)
|
#tested
|
IR (%)
|
#tested
|
IR (%)
|
IR (%)
|
#tested
|
Jul 2012
|
232
|
2.2
[0.7—4.4]
|
88
|
1.2
[0.0—4.3]
|
16
|
6.5
[0.2—22.1]
|
18
|
5.7
[0.2—19.8]
|
Oct 2012
|
165
|
3.7
[1.4—7.0]
|
180
|
5.0
[2.3—8.6]
|
81
|
2.5
[0.3—6.8]
|
69
|
4.4
[0.9—10.3]
|
Nov 2012
|
75
|
4.0
[0.9—9.5]
|
74
|
1.4
[0.0—5.1]
|
22
|
0.5
[0.0—4.4]
|
18
|
0.6
[0.0—5.3]
|
Apr 2013
|
7
|
1.4
[0.0—13.0]
|
12
|
0.8
[0.0—7.8]
|
4
|
2.3
[0.0—22.1]
|
2
|
4.3
[0.0—41.9]
|
Apr 2014
|
158
|
3.8
[1.4—7.3]
|
77
|
3.9
[0.8—9.2]
|
65
|
4.7
[1.0—10.9]
|
56
|
0.2
[0.0—1.8]
|
Aug 2014
|
124
|
1.6
[0.2—4.5]
|
131
|
0.1
[0.0—0.8]
|
160
|
5.0
[2.2—8.9]
|
163
|
2.5
[0.7—5.3]
|
Oct 2014
|
5
|
1.8
[0.0—17.9]
|
10
|
20.1
[2.9—48.4]
|
2
|
4.3
[0.0—41.9]
|
0
|
0.0
|
Apr 2015
|
41
|
0.3
[0.0—2.4]
|
59
|
6.8
[1.9—14.4]
|
5
|
1.9
[0.0—18.0]
|
2
|
50.0
[2.7—97.2]
|
Jun 2015
|
37
|
5.5
[0.7—14.7]
|
65
|
9.3
[3.5—17.3]
|
20
|
20.1
[6.1—39.6]
|
13
|
0.8
[0.0—7.2]
|
Nov 2015
|
18
|
0.6
[0.0—5.3]
|
0
|
0
|
1
|
7.9
[0.0—72.7]
|
0
|
0.0
|
Apr 2016
|
2
|
4.3
[0.0—41.8]
|
6
|
1.6
[0.0—15.1]
|
0
|
0.0
|
0
|
0.0
|
Jun 2016
|
45
|
8.9
[2.6—18.7]
|
101
|
3.0
[0.6—7.1]
|
89
|
2.3
[0.3—6.2]
|
30
|
0.4
[0.0—3.3]
|
TOTAL
|
909
|
3.1
[2.1—4.3]
|
803
|
3.6
[2.4—5.0]
|
465
|
4.3
[2.7—6.3]
|
371
|
2.4
[1.1—4.2]
|
EIRs showed monthly variations (Fig. 5). In Dangassa, the EIRs were estimated to be 13.74 [BCI: 9.21—19.14] and 18.66 [BCI: 12.55—25.81] ib/p/m, for indoor and outdoor locations, respectively. The EIR was 1.4 times higher outdoors than indoors, though the difference was not significant. The highest rates were observed in July and October 2012. In Koïla Bamanan, EIRs were 1.57 [BCI: 2.34 —2.72] and 0.94 [BCI: 0.43—1.64] ib/p/m for indoor and outdoor locations, respectively. The highest rates were observed for August 2014 and June 2015. In contrast to Dangassa, the EIR was 1.7 times higher indoors than outdoors. However this difference was not significant