Abundance and seasonal distribution of malaria vectors
A total of 50,322 mosquitoes belonging to four different genera were collected from all the study sites during the sampling seasons in 2017 (n = 26415) [Anopheline = 18153/26415, Culicine = 5856/26415, Mansonia = 2406/26415] and 2018 (n = 23907) [Anopheline = 15181/23907, Culicine = 6614/23907, Mansonia = 2112/23907]. Overall, 66.24%, (33,334/50,322) of mosquitoes sampled were Anophelines, 24.78%, (12,470/50,322) were Culicine, 8.98%, (4,736/50,322) were Mansonia mosquitoes. Among the Anophelines were Anopheles gambiae s.l. 93.16%, (31055/33,334), An. pharoensis 4.63%, (1545/33,334), An. funestus 2.07%, (690/33,334), and An. rufipes 0.13% (44/33,334).
Overall, 76.40%, (25,468/33,334) of Anophelines were collected in the rainy season and 23.60%, (7,866/33,334) in the dry season. Across the different ecological zones, the abundance of An. gambiae s.l. were high in the Sahel savannah [2017 (7,818/17,245); 2018 (6,883/13,810)] followed by the coastal savannah [2017 (7,110/17,245); 2018 (5,695/13,810)], and then the forest zone [2017 (2,317/17,245); 2018 (1,232/13,810)]. An. gambiae s.l. were more abundant in the rainy season (75.51%; 23,449/31,055) than in the dry season (24.49%; 7,606/31,055) during both sampling years (z = -36.037; p < 0.0001). In all the study sites, An. gambiae s.l. were more abundant in the rainy season than in the dry season except in Anyakpor and Dwease where An. gambiae s.l. were more abundant in the dry season (60.50%; 3082/5094), (68.83%; 848/1232) respectively, than in the rainy season during the 2018 sampling period. The most abundant other Anopheline mosquitoes were sampled from Kpalsogu = 692 [An. pharoensis = 641/692, An. funestus = 44/692, An. rufipes = 7/692]; Anyakpor = 151 [An. pharoensis = 149/151, An. funestus = 0/151, An. rufipes = 2/151]; Pagaza = 48 [An. pharoensis = 11/48, An. funestus = 35/48, An. rufipes = 2/48]; Dodowa = 15 [An. pharoensis = 12/15, An. funestus = 1/15, An. rufipes = 2/15]; Dwease = 2 [An. pharoensis = 2/2, An. funestus = 0/2, An. rufipes = 0/2] during the 2017 sampling season, and Pagaza = 563 [An. pharoensis = 25/563, An. funestus = 537/563, An. rufipes = 1/563); Anyakpor = 391 [An. pharoensis = 388/391, An. funestus = 0/391, An. rufipes = 3/391]; Kpalsogu = 351 [An. pharoensis = 310/351, An. funestus = 15/351, An. rufipes = 26/351]; Dwease = 58 [An. pharoensis = 0/58, An. funestus = 58/58, An. rufipes = 0/58]; Dodowa = 8 [An. pharoensis = 7/8, An. funestus = 0/8, An. rufipes = 1/8] during the 2018 sampling season. The most predominant species sampled were An. pharoensis [2017(815/908), 2018 (730/1371)], followed by An. funestus [2017(80/908), 2018 (610/1371)], and then An. rufipes [2017(13/908), 2018 (31/1371)]. High abundance of the other Anopheline species were sampled during rainy season [2017 (pharoensis = 681/908, funestus = 80/908; rufipes = 10/908); 2018 (pharoensis = 635/1371, funestus = 609/1371, rufipes = 4/1371], as compared to the dry season [2017 (pharoensis = 134/908; funestus = 0/908; rufipes = 3/908); 2018 (pharoensis = 95/1371, funestus = 1/1371, rufipes = 27/1371)].
Indoor and outdoor distribution of vectors
Overall, there was a significant difference (z = 2.410; p = 0.0160) between indoor biting (51.09%; 15,866/31,055) and outdoor biting An. gambiae s.l. (48.91%; 15,189/31,055). The abundance of An. funestus biting outdoor (59.28%; 409/690) was higher than those that were caught biting indoor (40.72%; 281/690) but the difference was not significant (z = -1.928; p = 0.0538).
In 2017, 48.19%, (8,311/17,245) An. gambiae s.l. were collected indoors and 54.70% (7555/13,810) in 2018 as shown in Fig. 2, at all the study sites. The abundance of indoor biting An. gambiae s.l. increased in the 2018 sampling period except in Dwease in the forest area where the abundance of indoor An. gambiae s.l. reduced from 53.73% (1245/2317) in 2017 to 50.57% (623/1232) 2018. During the 2017 sampling period, Anopheles gambiae s.l. preferred to bite indoors in all sites, except in Dodowa where there was more outdoor biting (59.64%; 3,139/5,263) than indoor biting (40.36%; 2,124/5,263) (Fig. 2).
Species discrimination in the Anopheles gambiae complex
A subsample of 1,670 of An. gambiae s.l. from all the study sites were randomly selected and used to discriminate the sibling species. An. coluzzii 55.9% (935/1670) was the most abundant species followed by An. gambiae s.s 39.5% (659/1670), An. arabiensis 2.3% (39/1670) and An. melas 2.2% (37/1670). Overall, An. coluzzii was the abundant species sampled in all the ecological zones, except in the Sahel savannah zone, where the species were dominated by the An. gambiae s.s. [Sahel (An. gambiae s.s. = 323/1670; An. coluzzii = 298/1670; An Arabiensis = 39/1670; An. melas = 0/1670); coastal (An. gambiae s.s. = 288/1670; An. coluzzii = 432/1670; An Arabiensis = 0/1670; An. melas = 37/1670); forest (An. gambiae s.s. = 48/1670; An. coluzzii = 205/1670; An Arabiensis = 0/1670; An. melas = 0/1670)].
The composition and distribution of these species differed significantly by study sites and year (Fig. 3; p < 0.001). An. coluzzii was most abundant in Anyakpor (83.65%, 394/471), Dwease (81.03%, 205/253) and Kpalsogu (80.20%, 239/298) while An. gambiae s.s was the most abundant in Dodowa (86.71%, 248/286) and Pagaza (74.03%, 268/362) respectively. All the An. melas collected during the study were only from the coastal savannah site of Anyakpor (7.86%, 37/471). All An. arabiensis were from Sahel savannah sites of Kpalsogu (1.34%, 4/298) and Pagazaa (9.67%, 35/362) (Fig. 3).
Biting times of An. gambiae s.l in the study sites
Anopheles gambiae s.l. were found to bite the most during the late evening (LE) (66.61%; 20685/31055), followed by the early morning (EM) (20.05%; 6228/31055), and less during the early evening (EE) (13.34%; 4142/31055). This was the same pattern for both indoor [LE: (66.12%, 10490/15866); EM: (20.52%, 3255/15866); EE: (11.42%, 2121/18566)] and outdoor biting [LE: (67.12%,10195/15189); EM: (19.57%, 2973/15189), EE: (13.96%, 2021/15189)] (Fig. 4).
With respect to the species of the An. gambiae sl., An. coluzzii and An. melas had a different biting pattern from An. gambiae s.s. and An. arabiensis. An. coluzzii preferred late evening feeding (48.02%; 449/935) followed by early morning (26.52%; 248/935) and early evening (25.45%; 238/935) feeding. Anopheles melas on the other hand preferred to bite in the early evening (43.24%; 16/37) followed by the late evening (40.54%; 15/37) and the early morning (16.22%; 6/37). Anopheles gambiae s.s preferred to bite in the late evening (54.48%; 359/659) followed by the early evening (25.04%; 165/659) and the early morning (20.46%; 135/659). Anopheles arabiensis preferred late evening (41.03%; 16/36) biting followed by early evening (33.33%; 13/36) and early morning (25.64%; 10/36) biting.
Insecticide resistance genotypes in Anopheles gambiae s.l.
Anopheles gambiae s.l samples were genotyped for the presence of Vgsc-1014S and 1014F mutations as well as the G119S mutation. Overall, the frequency of the Vgsc-1014F mutation was slightly higher in indoor biting mosquitoes (54.90%) than those biting outdoors (45.10%). Overall, Anopheles melas had a high Vgsc-1014F mutation (87.1%), whiles An. arabiensis had (50%) (Table 2).
Table 1
| An. gambiae s. s | An. coluzzii | An. melas | An. arabiensis | Total |
No (%) | No (%) | No (%) | No (%) |
Anyakpor | 40 (8) | 394 (84) | 37 (8) | 0 | 471 |
Dodowa | 248 (87) | 38 (13) | 0 | 0 | 286 |
Dwease | 48 (19) | 205 (81) | 0 | 0 | 253 |
Kpalsogu | 55 (18) | 239 (80) | 0 | 4 (1) | 298 |
Pagazaa | 268 (74) | 59 (16) | 0 | 35 (10) | 362 |
Total | 659 (39) | 935 (56) | 37 (2) | 39 (2) | 1670 |
Table 2
Species of the Anopheles gambiae complex and frequency of kdr L1014F and Ace-1 G119S mutation
| N | Kdr L1014F | n | F (kdr) | N | Ace 1 G119S | n | F (Ace-1) |
Indoor | 873 | RR | 328 | 0.57 | 870 | RR | 178 | 0.54 |
| | RS | 338 | | | RS | 579 | |
| | SS | 207 | | | SS | 113 | |
Outdoor | 752 | RR | 266 | 0.55 | 777 | RR | 120 | 0.48 |
| | RS | 302 | | | RS | 512 | |
| | SS | 184 | | | SS | 145 | |
Species | | | | | | | | |
An. arabiensis | 36 | RR | 13 | 0.50 | 39 | RR | 2 | 0.44 |
| | RS | 10 | | | RS | 30 | |
| | SS | 13 | | | SS | 7 | |
An. coluzzii | 908 | RR | 300 | 0.53 | 924 | RR | 187 | 0.52 |
| | RS | 363 | | | RS | 591 | |
| | SS | 245 | | | SS | 146 | |
An. gambiae s.s | 646 | RR | 255 | 0.59 | 648 | RR | 96 | 0.50 |
| | RS | 258 | | | RS | 451 | |
| | SS | 133 | | | SS | 101 | |
An. melas | 35 | RR | 26 | 0.87 | 36 | RR | 13 | 0.63 |
| | RS | 9 | | | RS | 19 | |
| | SS | 0 | | | SS | 4 | |
F(kdr) = 2RR + RS/2n Ahadji-Dabla 2019. F: allelic frequency, N = number of samples tested, n = total number of samples positive for a specific genotype |
Similarly, the frequency of the G119S mutation in An. gambiae s.l were slightly higher in indoor host seeking mosquitoes (52.94%) than outdoor biting mosquitoes (47.06%). Anopheles melas had a high frequency of G119S mutation 0.63 (30.14%), whiles An. arabiensis had the lowest frequency of 0.44 (21.01%) (Table 3). The Vgsc-1014S mutation was not detected in any mosquito for this study.
Table 3
Sporozoite rates in the An. gambiae s.l
Community | Number per pools | Total Pools | Positive pools | Dry season | Rainy season | Indoor | Outdoor | Sporozoite Rate (SR) | Dry season SR | Rainy season SR | Indoor SR | Outdoor SR |
Anyakpor | 20 | 164 | 8 | 4 | 4 | 7 | 1 | 0.24 | 0.12 | 0.12 | 0.21 | 0.03 |
Dodowa | 20 | 105 | 7 | 5 | 2 | 3 | 4 | 0.33 | 0.24 | 0.10 | 0.14 | 0.19 |
Dwease | 20 | 103 | 10 | 5 | 5 | 6 | 4 | 0.49 | 0.24 | 0.24 | 0.29 | 0.19 |
Pagazaa | 20 | 131 | 11 | 4 | 7 | 6 | 5 | 0.42 | 0.15 | 0.27 | 0.23 | 0.19 |
Kpalsogu | 20 | 140 | 8 | 4 | 4 | 3 | 5 | 0.29 | 0.14 | 0.14 | 0.11 | 0.18 |
| | 643 | 44 | | | | | | | | | |
* Sporozoite Rate SR= ((no of positive pools) / (no of pools × maximal pool size)) × 100 Wei-Dong Gu 1995 |
Sporozoite infection rates in the sampled vectors
A total of 12,860 Anopheles gambiae s.l were pooled in groups of 20, into 643 pools and tested for Plasmodium falciparum Circumsporozoite (CSP). Overall, forty-four pools (44) were positive for P. falciparum CSP (Anyakpor (n = 8), Dodowa (n = 7), Dwease (n = 10), Kpalsogu (n = 8), and Pagazaa (n = 11)) giving an overall sporozoite rate of 0.07% (Table 3).
Dwease had the highest sporozoite rate (0.49%), followed by Pagazaa (0.42%), Dodowa (0.33%), and Kpalsogu (0.29%), with Anyakpor reporting the lowest sporozoite rate (0.24%). Putting all the sites together, sporozoite rate was quite similar in indoor biting An. gambiae s.l. (0.04%) compared to those that were biting outdoors (0.03%). With respect to the individual study sites, sporozoite rate was higher indoors compared to outdoors except in Dodowa [(indoor (0.14%); outdoor (0.19%)] and Kpalsogu [(indoor (0.11%); outdoor (0.18%)] where outdoor biting An. gambiae s.l. had a higher sporozoite rate.
Furthermore, the sporozoite rates were similar during both the dry and rainy seasons in Anyakpor [dry = 0.12%; rainy = 0.12%], Dwease [dry = 0.24%; rainy = 0.24%], and Kpalsogu [(dry = 0.14%; rainy = 0.14%]. All the study sites had similar sporozoite rates in both seasons as represented in Table 2 except in Dodowa where the sporozoite rate was higher in the dry season (0.24%) compared to the rainy season (0.10%), and in Pagazaa where the rainy season had a higher sporozoite rate (0.27%) compared to the dry season (0.15%).