- Open Access
Transmission patterns of Plasmodium falciparum by Anopheles gambiae in Benin
Malaria Journalvolume 13, Article number: 444 (2014)
To better control malaria, the clear and urgent need is for improved data to inform decision makers, but in several African countries, there is a lack of baseline data on vectors and variation in the intensity of malaria transmission. This has resulted in the implementation of vector control efforts that ignore variation in vector behaviour and intensity of transmission, an approach that is most often not cost-effective. This study presents a detailed entomological description of mosquito distribution and variation in potentially transmissible contacts of Plasmodium falciparum following a south to north transect in Benin.
The study was conducted in five locations where environmental parameters were different and malaria prevalence ranged between 14 and 51%. The locations represent the main eco-epidemiological malaria areas in Benin. Mosquitoes were collected using human landing catches, pyrethrum spray catches and windows traps. They were taxonomically and molecularly identified. Head-thoraces of Anopheles gambiae s.l. were tested by enzyme-linked immunosorbent assay. Entomological indicators were estimated following WHO guidelines.
The results showed variation between location and period in distribution of Anopheles coluzzii, An. gambiae, and Anopheles arabiensis (p < 0.05). An extension of the reported range of An. arabiensis was also observed. Densities of malaria vectors varied significantly between rural and urban sites, however, indoor/outdoor biting ratios remained constant. Proportions of malaria vectors with circumsporozoite protein of P. falciparum were similar between locations. The entomological inoculation rates ranged between zero and eight bites/man/night with significant variations between areas.
Four profiles of human exposure to infectious malaria vector bites were observed and included location with one season of high transmission (June - August), two seasons of lower transmission (March-August; October-November), moderate continuous transmission season, and high continuous transmission season of P. falciparum.
The study revealed several entomological patterns in transmission of P. falciparum in Benin. The data could be used for purposes of planning a more cost-effective vector control strategy, by stratifying the country into higher and lower transmission zones. The information could also be used to guide extension of indoor residual spray based on a targeted use of IRS at sites where the duration of insecticidal effect following spraying coincides with the peak transmission period.
Malaria vector control strategies have demonstrated their effectiveness through several studies [1–4]. In endemic areas, the result is a significant reduction in malaria morbidity and overall mortality [5–7]. The number of deaths due to malaria decreased from one million in 2000 to 627,000 in 2012 . Several African countries have reported a decline in the number of confirmed malaria cases and the number of malaria-related hospitalizations and deaths. This significant progress in reducing malaria mortality and morbidity was mainly due to the Roll Back Malaria (RBM) initiative, which aimed to reduce at least by 50% the number of malaria deaths in 2010 and 75% by 2015 .
However, the development of malaria vector resistance to insecticides is an important threat to the effectiveness of vector control strategies . At the same time, lack of information on vector population dynamics, before and after the application of vector control measures, makes the job of ongoing evaluation and planning of vector control interventions difficult .
The priority in terms of malaria control is for improved data to guide policy making . But in most regions of Africa, there is a lack of baseline and monitoring information on vectors and variation in the intensity of malaria transmission, forcing program to use indoor residual spraying (IRS) and long-lasting insecticidal nets (LLINs) in uninformed ways, such as timing IRS such that it misses the peak transmission season. The implementation of these vectors control strategies requires detailed information on malaria transmission and distribution of vector populations. In addition, mosquito density and distribution varies between areas and periods and such variation can change the level of malaria transmission . Assessment of the impact of vector control interventions on malaria transmission requires that study sites be identified and described, with baseline data gathered before the implementation of vector control strategies.
In Benin, current vector control strategies are based on the use of insecticide-based interventions (IRS and LLINs) [1, 12], but the use of these tools could not lead to the elimination of malaria vectors when the deployment of intervention does not account for variability associated with abundance of vectors and intensity of malaria transmission. Existing data on these two indicators are specific, located in a few geographical areas and are not representative of an entire territory.
This study presents a detailed entomological description of the distribution of mosquitoes and the variation in the rate of potentially transmissible contacts of Plasmodium falciparum throughout malaria eco-epidemiological areas in Benin. This will help the National Malaria Control Programme to make evidence-based decisions on vector control.
The study was conducted in five eco-epidemiological locations following a south-to-north transect in Benin where malaria prevalence  or geography was different (Table 1). LLIN intervention was ongoing in all selected locations since a universal-coverage national distribution campaign in July 2011. In each location, two representative areas (urban and rural) were selected. The entomological surveys were conducted bimonthly at each area from March to November 2012 during the malaria transmission season, especially when vector density was likely to change due to climatic conditions.
Three classical methods were used to monitor vector population dynamics: human landing catches (HLCs) to evaluate man-vector contact, pyrethrum spray collections (PSCs) to assess indoor-resting behavior of vectors, and window exit traps (WETs) to estimate exophilic behaviour. Blood-feeding rates were estimated from mosquitoes collected using PSCs and WETs. Each survey consisted of two-day collections in each of the study sites.
Allada is located in the north of Atlantique Province. The climate is subequatorial with two rainy seasons (March-July and September-November) and two dry seasons (July-September and November-March). The annual average rainfall is between 800 and 1,000 mm. The soil is mainly characterized by earth bar and a marshy depression; it is well suited to food crops of vegetables and fruit and coffee growing. Most lands in Allada are less fertile. The selected areas were:
Allomey: an urban area of Allada located 6°40'14.71"N latitude and 2°9'2.99"E longitude.
Niaouli: a rural village located 6 °44'37.01"N latitude and 2°8'13.32"E longitude.
Dassa-Zoume is one of six districts of Collines Province. The climate is Guinean with four seasons: two rainy seasons and two dry seasons. The highest rainfall is recorded between July and September. The vegetation is forests and open savannahs. The two selected areas were:
Amangassa: an urban area of Dassa city, located 7°46'31.74"N latitude and 2°11'43.11"E longitude.
Lema: a rural village of Dassa, located 7°49'58.83"N latitude and 2°13'37.60"E longitude.
Parakou is located in the centre of Borgou Province. The climate is South Sudanian. It is characterized by a rainy season (May to October) and a dry season (November-April). The lowest temperatures are recorded between December and January. The average annual rainfall is 1,200 mm. Parakou is watered by tributaries of the Okpara. These rivers remain dry from January to May. The main vegetation is savannah. The selected areas were:
Zongo zeno: an urban area of Parakou, located 7°49'58.83"N latitude and 2°13'37.60"E longitude.
Kpassa: a rural village of Parakou, located 7°49'58.83"N latitude and 2°13'37.60"E longitude.
Kandi is located in the centre of Alibori Province in the agro-ecological zone of the cotton area. The climate of Kandi is Sudanese and characterized by two distinct seasons, a rainy season from April to October and a dry season from November to March. December, January and February are characterized by an absolute drought with harmattan. The annual average rainfall varies considerably between 700 and 1,400 mm. The soil is very rugged with predominance of plateau. It is watered by the tributaries of Alibori and Sota. The selected areas were:
Kossarou: an urban area of Kandi I, located 11°7'29.32"N latitude and 2°56'9.57"E longitude.
Sonsoro: a rural village of Kandi, located 11° 4'58.91"N latitude and 2°13'37.60"E longitude.
Malanville is located in the far north of the Republic of Benin. The climate is Sudano-Sahelian type and characterized by a dry season from November to April. The average rainfall is 750 mm. Malanville is watered from east to west by the Niger River and its tributaries Alibori, Mekrou and Sota, which are in flood during August and September. The main vegetation is savannah with grassland. The principal agricultural occupation is rice cultivation. The included areas were:
Arobanda: an urban area, located 11°51'5.12"N latitude and 3°23'45.06"E longitude.
Bodjecali: a rural village, located 11°48'49.06"N latitude and 3°22'58.08"E longitude.
Mosquito sampling and field processing
Human Landing Catches (HLCs) were conducted in two houses in each area from 21.00 to 06.00 hours. At each house, mosquito collectors, one outside and one inside, with a flashlight and a mouth aspirator, collected all landing mosquitoes on their feet. Assuming that there is a risk of malaria transmission, collectors received an antimalarial prophylaxis as prevention against malaria. Malaria symptoms were also monitored and treatment was provided when necessary.
A team of two people estimated indoor-resting and outdoor-exiting density using PSCs and WETs in two other houses between 06.45 and 08.00 hours. PSCs consisted of covering all exposed surfaces with white sheets, spraying the rooms and collecting all fallen specimens.
After collection, mosquitoes were counted and morphologically identified [14, 15]. A proportion of unfed Anopheles gambiae s.l. females from HLCs were dissected to extract ovaries and to determine parous An. gambiae s.l. (mosquitoes that have oviposited at least once) by observing the coiling degree of ovarian tracheoles . Collected females from PSCs and WETs were classified by abdominal status to estimate blood-feeding rate. Collected An. gambiae s.l. were stored individually in labelled Eppendorf tubes with desiccant until laboratory processing.
Around 100–150 An. gambiae s.l. were randomly selected (from urban and rural areas and proportionally to the number collected in each area) during each collection period for detection of P. falciparum circumsporozoite protein (CSP). The heads and thoraces of the An. gambiae s.l. females selected were tested by ELISA for identification of P. falciparum CSP according to the method described by Wirtz and colleagues . Additionally, for each survey, a random sample of 48 females of An. gambiae s.l. by district was identified to species level by molecular method described by Fanello and colleagues . When the number of An. gambiae s.l. collected during an assessment did not reach 48 in a district all the An. gambiae s.l. collected were used.
Composition of mosquito populations was studied in terms of species diversity, expressed as richness (Taxa S) and species abundance, expressed as evenness . Shannon-Wiener diversity index and Simpson's dominance index was determined at each site.
Human-biting rate (HBR) was defined as the ratio of the total number of mosquito species collected to the total person-nights for the collection period. Endophagy rate represented the proportion of mosquitoes collected indoors against the total of both indoors and outdoors collections from HLC. Exophagy rate was the proportion of mosquitoes collected outdoors against the total of both indoors and outdoors collections from HLC. Endophily was estimated as the percentage of An. gambiae s.l. collected by indoor residual (at rest) divided by the total number collected by PSCs and WETs. Exophily rate was determined as the percentage of An. gambiae s.l. collected by WETs (exiting) divided by the total number collected by PSCs and WETs. Parity rate was defined as the number of parous An. gambiae s.l. divided by the total number dissected. The sporozoite index (SI) was defined as the proportion of total mosquitoes collected found to contain the P. falciparum CSP. EIR represented the product of HBR and the SI of night caught mosquitoes.
Mosquito diversities were compared using Student’s t test. The diversity data were analysed using PAST 2.07. Poisson's confidence intervals method  was used to estimate confidence intervals of densities, HBR, EIR, endophagy rates and frequencies of different species of An. gambiae s.l. Binomial confidence intervals method  was used to estimate confidence intervals of parity rates. Unconditional maximum likelihood estimation method (Wald), or median unbiased estimation (mid-p) of the risk ratio (RR) followed by their confidence intervals obtained using the normal approximation or the exact method mid-p and p-values obtained by Khi2 or midp.exact , were used to compare endophagy and exophagy rates and EIR. Principal component analysis was used to reveal the structure of the EIR data to determine the profiles of malaria transmission. Factor 1, 2, and 3 represented principal component 1, 2, and 3.
Zero-inflated regression models  were used to predict entomological inoculation rates of unevaluated months due to the important number of zero in the dataset. This model used the EIR as the dependent variable and covariates as: locations and months of collection. This model is the best among the Poisson family models for the data based on the log likelihood. All these parameters were computed and analysed using the free software R version 2.15.1.
This study received ethical approval from the Ministry of Health in Benin. Informed consent was obtained from mosquito collectors who were vaccinated against yellow fever. Mosquito collectors received preventive treatments for malaria and were also subject to regular medical check-ups. Authorization to conduct the study was also obtained from district authorities and from individuals before entering their houses.
A total of 8,663 man-biting mosquitoes belonging to five genera (Anopheles, Aedes, Culex, Mansonia, Coquilletidia) and 21 species were collected in the study sites (Table 2). Two major malaria vectors were found: An. gambiae s.l. (found everywhere, and representing 71.1% of the total mosquitoes: 6,166/8,663) and Anopheles funestus found in Allada (0.05%) and Dassa (0.1%). Anopheles nili, a local vector, was collected but at very low numbers in Parakou (0.04%) and Malanville (0.1%). Other Anophelines collected and identified to species included: Anopheles pharoensis, Anopheles ziemanni, Anopheles broheri, Anopheles fuscicolor, and Anopheles coustani.
Culex quinquefasciatus and other Culex species were also found at low densities (Table 2). Culex quinquefasciatus was the main man-biting mosquitoes collected in Dassa (72%) and the most common non-anbopheline man-biting mosquito species in Parakou (22%) and Kandi (27%), but Mansonia africana was the second most common man-biting mosquito species found in Allada (26%) and Parakou (15%). Coquilletidia cristata was only found in Malanville and Dassa (Table 2).
Allada recorded Shannon-Wiener diversity and Simpson's dominance index of 1.66 and 0.27, respectively, while Dassa, Parakou, Kandi, and Malanville recorded diversity and dominance values of 0.65 and 0.59, 1.12 and 0.42, 0.78 and 0.56, 0.63 and 0.75, respectively. When the frequencies of diversity and dominance indices of mosquitoes were compared between sites, Allada recorded the highest Shannon-Wiener diversity index and the lowest Simpson's dominance index, while Malanville recorded the highest Simpson's dominance index but the lowest Shannon-Wiener diversity index regardless of the number of collected species in this area was high. In Allada, many species were collected but at very low densities. There was no significant difference in mosquito diversity between Allada, Parakou, Kandi, and Malanville (p > 0.05), but in the district of Dassa a significant difference (p < 0.001) in mosquito diversity was observed compared to other districts.
Distribution of Anopheles gambiae species complex
In total, 966 An. gambiae s.l. were analysed by PCR for species identification. Three species of the complex An. gambiae were identified: Anopheles coluzzii, An. gambiae and Anopheles arabiensis. Anopheles coluzzii represented 50% (n = 482) of mosquitoes tested against 16% (n = 158) of An. arabiensis, and 34% (n = 326) of An. gambiae. The three species were found at five sites but at different times during the transmission season (Figure 1). Anopheles arabiensis previously limited in Parakou was identified in Dassa and Allada. Anopheles arabiensis (43%) and An. coluzzii (48%) predominated between March and May at Allada. From July to November, density of An. arabiensis decreased to 16% while density of An. gambiae increased from 9 to 46%. From March to July, An. coluzzii and An. gambiae were the most represented species in Dassa (58 and 40%, respectively) and Kandi (36 and 64%). Anopheles arabiensis were found in a sizeable proportion (19-23%) at these two areas between September and November (Figure 1). Anopheles coluzzii was highly predominant (p < 0.05) in Malanville. In Parakou, An. gambiae predominated between May and July, while An. arabiensis was mostly found between September and November, and An. coluzzii in March (Figure 1).
Aggressive density of Anopheles gambiae s.l
A total of 143 (117.74-165.34) An. gambiae s.l. were collected in Allada, 219 (194.59-246.39) in Dassa, and 1,239 (1,184.32-1,296.88) in Parakou. Respectively, 653 and 3,912 An. gambiae s.l. were collected in Kandi and Malanville. The number of An. gambiae s.l. was generally higher in rural areas than in urban areas except in Malanville where the opposite was observed (Figure 2).
In Malanville the average HBR was 49 bites/man/night. The aggressive density of An. gambiae s.l. observed in Malanville was significantly higher than in Parakou (15 bites/man/night), Kandi (eight bites/man/night), Dassa (three bites/man/night), and Allada (two bites/man/night). Endophagy rates and exophagy rate were similar in Malanville, Kandi Dassa, and Allada (p > 0.05). Exophagy rate of An. gambiae s.l. was significantly higher than the endophagy rate (p < 0.001) at Parakou (Figure 3).
Aggressive density of other Culicidae
The abundance of other mosquito species was significantly lower compared to that of An. gambiae s.l. (p < 0.05) but also varied from one site to another. Unlike An. gambiae s.l., other Culicidae were more abundant in urban than in peripheral environment in all districts except Parakou where other Culicidae were significantly higher (501 (462.50-542.303)) in peripheral areas (349 (317.19-383.87)) (Figure 4). The number of other Culicidae ranged from 138 (118.31-160.66) to 394 (356.34-432.57) in urban areas in Allada, Dassa, Kandi, and Malanville. In rural areas of these districts, the number of other Culicidae was between 57 (44.80-72.25) and 216 (189.13-244.85). Fifty-five (41.94-69.09) to 345 (313.26-379.78) human-biting other Culicidae were collected indoors (Figure 5). The number of other mosquito species collected outdoors varied from 123 (108.78-140.80) to 348 (318.65-375.24) (Figure 5). The endophagy rates of other Culicidae were significantly lower than exophagy rates in all districts (p < 0.05) except in Kandi where endophagy rate was significantly higher than exophagy rate (p < 0.05).
Longevity of Anopheles gambiae s.l
Some 1,810 ovaries of An. gambiae s.l. were dissected. The majority of An. gambiae s.l. dissected were old with an average parity rate of 88.06%. The parity rate was 88.11% in Allada and 76.87% in Dassa (Table 3). In Parakou, Kandi and Malanville, the observed parity rates were 88.60, 90.32 and 88.34%, respectively. No significant variation was observed between sites (p > 0.05).
Exophily of Anopheles gambiae s.l., Mansonia spp. and Culex spp
A total of 2,677 mosquitoes were collected in morning PSCs and WETs (Table 4). They belonged to three genera (Anopheles including An. gambiae s.l. only, Mansonia and Culex). 1,167 An. gambiae s.l. were collected versus 1,510 Mansonia and Culex spp (Table 4). Observed exophily of An. gambiae s.l. was generally low (Table 4). Exophily rates for An. gambiae s.l. in Allada, Parakou, and Kandi were similar (p > 0.05) with respective exophily rates of 17.24, 17.65, and 21.34 (Table 4). The highest exophily rate was observed in Dassa (52.17%, p = 0.03). The lowest exophily rate was observed in Malanville (2.95%, p < 0.0001). Observed exophily rates in Mansonia and Culex spp in Parakou (30.65%) and Malanville (1.41%) was significantly lower than that observed in other districts (p < 0.05). Exophily rates in Allada, Dassa and Kandi were 57.58, 39.60, and 39.20%, respectively.
Mosquito blood-feeding rates
The observed blood-feeding rates of mosquitoes collected from indoor residual and WETs in different study sites were high (Table 5). Anopheles gambiae s.l. had a higher blood-feeding rate in Allada (100%) and Malanville (93.65%) than in Dassa (56.52%) and Kandi (73.87%). In Parakou, An. gambiae s.l. had a similar blood-feeding rate (80.39%) to all other districts (Table 5). A significant variation in blood-feeding behaviour of An. gambiae s.l. was observed between sites (p < 0.05).
Mansonia spp and Culex spp also had high blood-feeding rates (Table 5). The observed blood-feeding rates among these species varied between sites (p < 0.05) and ranged from 48 to 93%.
Infectivity of Anopheles gambiae s.l
A total of 1,652 head-thoraces of An. gambiae s.l. collected were assessed with CSP-ELISA assays. The mean SI was 4.46%. SI was similar (p > 0.05) in urban (4.3%: 24 head-thoraces+/556 head-thoraces) and rural areas (4.5%: 50 head-thoraces +/1,102 head-thoraces). The average SI was 5.8% (31 head-thoraces +/531 head-thoraces) in Parakou, 5.6% (12 head-thoraces +/214 head-thoraces) in Dassa, 2.5% (seven head-thoraces +/280 head-thoraces) in Kandi, 4.4% (eight head-thoraces +/182 head-thoraces) in Allada, and 3.5% (16 head-thoraces+/451 head-thoraces) in Malanville. The highest infection rate was observed in urban Dassa, and the lowest was observed in rural Kandi (Table 6). In urban Kandi, no positive head-thorax was found. Plasmodium falciparum infection was observed in An. gambiae s.l. from March to November in urban Malanville and rural Parakou. In the other locations P. falciparum infection rate was found variable over time.
Entomological inoculation rate of Anopheles gambiae s.l
Table 7 shows the distribution of EIR by site. The EIR of An. gambiae s.l. was high in Malanville, with an average range of four to eight infected bites/man/night during the study period. Infected bites were recorded each month from March to November at this location. Parakou was also characterized by a high malaria transmission with two to five infected bites/man/night in rural and urban areas. The lowest EIR was observed in urban Kandi with no infected bites/man/night. Globally, observed EIRs were significantly higher in rural areas than in urban areas (p < 0.05), except in Allada where EIRs were similar (Table 7). EIR was generally higher in July than in other months (Table 7).
Principal component analysis identified three significant factors that accounted for around 90% of the variation in EIR (Table 8). Observed human exposure to infectious bites in Kandi, Dassa and Allada significantly contributed to the formation of factor 1, transmission risk in Parakou essentially contributed to the formation of factor 2, and those of Malanville were decisive in factor 3 (Table 8).
Figures 6, 7 and 8 describe the relationship between observed EIRs and these factors. The results suggest four profiles of malaria transmission. The first profile includes Kandi and Dassa where variations in EIRs were similar (Figure 6). In contrast, Allada had a different profile suggesting that variation in entomological inoculation rates was significantly different from those observed in Dassa and Kandi (Figure 7). Variations in EIRs in Malanville and Parakou were significantly different, suggesting two different malaria transmission profiles (Figure 8). The existence of such profiles was further confirmed by a monthly EIR prediction model in each site (Figure 9). The predicted EIR model showed four trends (profiles) in malaria transmission. The first trend included transmission predicted in Dassa and Kandi that showed one season of high transmission. This season begins in June and ends in August with a maximum of 13 infectious bites/man/month in Dassa versus 15 infectious bites in Kandi (Figure 9). The second trend was determined by the transmission predicted in Allada where transmission is lower than in other localities but with two seasons. The first season begins in March and ends in August with a maximum of four infectious bites/man/month. The second transmission season occurs in October and ends in November with a maximum of six infectious bites/man/month. The third transmission trend was determined by EIR predicted in Parakou with three peaks suggesting a continuous malaria transmission in this area. In this area, EIRs varied from nine to 27 infectious bites/man/month (Figure 9). The fourth trend observed in Malanville showed several peaks of EIR ranging from nine to 120 infectious bites/man/month (Figure 9). It also suggests a continuous transmission in this area but with higher value than in Parakou.
The majority, seventy one percent, of the vectors collected during this study were An.gambiae s.l. Abundance varied according to the degree of urbanization. Estimates of indoors versus outdoors biting were similar. The number of An. gambiae s.l. infective bites/human/night varied between locations. Based on these results four P. falciparum transmission profiles, across a North–south transect in Benin, are described.
No significant variation in vector bio-diversity was observed between locations, suggesting that similar control intervention would work in each location, thereby reduce higher costs associated with the application of different control methods in different locations . Taxonomic results also demonstrate a spatio- temporal variation in the complex An. gambiae as observed in Nigeria, Cameroon and Burkina Faso [24–26].
Anopheles arabiensis, a vector of dry savannah areas and limited in 1980 between the extreme north Benin (Malanville) and Parakou [27, 28], was found in Dassa and Allada in 2012. The results suggest an extension of the spatial distribution of this vector associated with destruction of forests and large areas (70,000 ha/yr) for agriculture, hunting, grazing, and wood harvesting in central and southern Benin . Habitat changes, such as these, may contribute to the spread of this malaria vector in southern Benin. In addition, long dry seasons recently observed in Benin , and elsewhere in Africa, due to climate change , foster the spread of An. arabiensis that has a similar abudance with An. gambiae during certain months (May, September and November).
Differences between densities of An. gambiae s.l. and other man-biting mosquitoes in each site could be explained by the ecological characteristics of the environment. Many of the study sites, such as Malanville and Parakou, are suitable for vector breeding. They are lowland areas where rice and cotton cultivation, and vegetable farming are the main activities [31, 32]. Densities of An. gambiae s.l. were similar inside and outside houses. This result was contrary to that observed by Akogbeto and colleagues  showing a low endophagy rate in indoor residual spray areas. It highlighted the importance of such strategies on malaria vector control. High densities observed in urban Malanville are due to the existence of rice-growing areas with permanent mosquito breeding sites in the urban area .
The EIR estimates human exposure to malaria-infected vectors. It is commonly used to assess malaria endemicity and transmission intensity . In the present study, observed EIRs varied significantly between locations. The highest rates were observed in Parakou and Malanville with more than two infectious bites/man/night suggesting that populations from Malanville and Parakou are more exposed to transmission of P. falciparum compared to other locations. Observed malaria transmission risk at each location could be associated with their environmental characteristics. Koudou and colleagues  reported an increase EIR from 38 to 295 infective bites/person/year with a return to rice cultivation at Zitta region in Côte d’Ivoire. Permanent exposure to malaria for a community close to vegetable farming was also reported by Yadouleton et al. . Such observations suggest that rice cultivation and vegetable farming are key factors associated with high malaria transmission risk and could explain transmission observed in Malanville and Parakou. Rice cultivation could also explain differences in EIR between Dassa and Allada, which are both highlands with similar malaria prevalence, and between Kandi and Malanville which have similar malaria prevalence. In addition, Kandi, Diassa, and Allada are highlands reported to have lower malaria transmission risk than the lowlands (Malanville and Parakou) . The impact of the current LLINs intervention on vector populations may also explain the lower vector densities and lower transmission risks than expected in some areas where case data-based malaria prevalence does not seem to correlate with EIRs [3, 37]. But a firm conclusion on correlation between EIRs and prevalence will be drawn if entomological and parasitological data are concurrently collected.
Infectious-vectors bites were higher in rural areas compared to urban areas in a same district suggesting a spatial heterogeneity of EIRs. This shows the importance of local conditions or micro-ecological conditions and environment characteristics in the intensity of malaria transmission.
High EIRs observed in certain locations was particularly due to a high HBR of An. gambiae, but not to a high infectious rate of vectors. This observation raises the question: where are we most vulnerable to malaria: between areas with low densities of vectors with high infectious rates or areas with high densities of vectors with low infectious rates?
Observed patterns in transmission of P. falciparum included locations with one season of high malaria transmission, locations with two seasons of high transmission and locations with continuous high transmission of malaria. Malanville and Parakou are locations at high risk of transmission where urgent, efficient vector control strategies should be implemented to help populations.
The study demonstrates the role of vector biology in the epidemiology of malaria in Benin. However, the lack of data on malaria prevalence estimated by microscopy or rapid detection in populations of each study site represents an important limitation. The study could also be improved by increasing the number of surveys and including more locations.
The study of intensity of P. falciparum transmission by mosquitoes is very important for implementation of insecticide-based vector control interventions in Benin where several large-scale vector control interventions (LLINs and IRS) are ongoing. It will guide policies in the selection of priority areas for vector control interventions. Knowledge of the number of high transmission season in each area will also help to implement cost-effective IRS interventions.
This study documents at least four vector-biology-associated patterns in transmission of P. falciparum from southern to northern Benin. It shows a spatiotemporal variation in the distribution of malaria vectors. It also reveals a spatial heterogeneity in human-biting behaviour as well as variation in human exposure to infectious bites of malaria vectors.
These data represent a source of information that could guide implementation of any cost-effective vector control strategy in Benin by targeting the higher transmission areas, that is by stratifying the country into higher and lower transmission zones. The study could also help extension of indoor residual spray based on a targeted use of IRS at sites where the duration of insecticidal effect following spraying coincides with the peak transmission period.
Akogbeto M, Padonou GG, Bankole HS, Gazard DK, Gbedjissi GL: Dramatic decrease in malaria transmission after large-scale indoor residual spraying with bendiocarb in Benin, an area of high resistance of Anopheles gambiae to pyrethroids. Am J Trop Med Hyg. 2011, 85: 586-593. 10.4269/ajtmh.2011.10-0668.
West PA, Protopopoff N, Rowland MW, Kirby MJ, Oxborough RM, Mosha FW, Malima R, Kleinschmidt I: Evaluation of a national universal coverage campaign of long-lasting insecticidal nets in a rural district in north-west Tanzania. Malar J. 2012, 11: 273-10.1186/1475-2875-11-273.
Sovi A, Azondékon R, Aïkpon RY, Govoétchan R, Tokponnon F, Agossa F, Salako AS, Oké-Agbo F, Aholoukpè B, Okè M, Gbénou D, Massougbodji A, Akogbéto M: Impact of operational effectiveness of long-lasting insecticidal nets (LLINs) on malaria transmission in pyrethroid-resistant areas. Parasit Vectors. 2013, 6: 319-10.1186/1756-3305-6-319.
Tokponnon FT, Aholoukpe B, Denon EY, Gnanguenon V, Bokossa A, N’guessan R, Oke M, Gazard DK, Akogbeto MC: Evaluation of the coverage and effective use rate of long-lasting insecticidal nets after nation-wide scale up of their distribution in Benin. Parasit Vectors. 2013, 6: 265-10.1186/1756-3305-6-265.
Lengeler C: Insecticide-treated bed nets and curtains for preventing malaria. Cochrane Database Syst Rev. 2004, 2: CD000363-
WHO: World Malaria Report 2012. 2012, Geneva: World Health Organization, [http://www.who.int/malaria/publications/world_malaria_report_2012/report/en/]
WHO: World Malaria Report 2013. 2013, Geneva: World Health Organization, [http://www.who.int/malaria/publications/world_malaria_report_2013/en/]
WHO: WHO / World Malaria Report 2010. 2010, Geneva: World Health Organization
Killeen GF, Okumu FO, N’Guessan R, Coosemans M, Adeogun A, Awolola S, Etang J, Dabiré RK, Corbel V: The importance of considering community-level effects when selecting insecticidal malaria vector products. Parasit Vectors. 2011, 4: 160-10.1186/1756-3305-4-160.
Marsh K: Research priorities for malaria elimination. Lancet. 2010, 376: 1626-1627. 10.1016/S0140-6736(10)61499-7.
Munga S, Yakob L, Mushinzimana E, Zhou G, Ouna T, Minakawa N, Githeko A, Yan G: Land use and land cover changes and spatiotemporal dynamics of Anopheline larval habitats during a four-year period in a highland community of Africa. Am J Trop Med Hyg. 2009, 81: 1079-1084. 10.4269/ajtmh.2009.09-0156.
Gnanguenon V, Azondekon R, Oke-Agbo F, Beach R, Akogbeto MC: Durability assessment results suggest a serviceable life of two, rather than three, years for the current long-lasting insecticidal (mosquito) net (LLIN) intervention in Benin. BMC Infect Dis. 2014, 14: 69-10.1186/1471-2334-14-69.
MDAEP, INSAE, Bénin ICF International: Enquête Démographique et de Santé (EDSB-IV) 2011–2012. 2013, Cotonou: Institut National de la Statistique et de l’Analyse Économique (INSAE) - Ministère du Développement, de l’Analyse Économique et de la Prospective, Gouvernement du Bénin
Gillies MT, De Meillon B: The Anophelinae of Africa South of the Sahara. 1968, Johannesburg: Publication of the South African Institute for Medical Research
Gillies MT, Coetzee M: A supplement to the anophelinae of Africa south of the Sahara. Publ South Afri Inst Med Res. 1987, 55: 143-
Detinova T, Gillies M: Observations on the determination of the age composition and epidemiological importance of populations of Anopheles gambia e Giles and Anopheles funestus Giles in Tanganyika. Bull World Health Organ. 1964, 30: 23-28.
Wirtz R, Ballou W, Schneider I, Chedid L, Gross M, Young J, Hollingdale M, Diggs C, Hockmeyer W: Plasmodium falciparum: immunogenicity of circumsporozoite protein constructs produced in Escherichia coli. Exp Parasitol. 1987, 63: 166-172. 10.1016/0014-4894(87)90158-5.
Fanello C, Della TA: Distribution du gène Kdr dans les formes chromosomiques et moléculaires des Anopheles gambiae s.s en Afrique Occidentale. 2001, Geneva: Rapport MIM/AFRO/OMS/TDR
Simpson EH: Measurement of diversity. Nature. 1949, 163: 688-688. 10.1038/163688a0.
Kenneth R: Epidemiology An Introduction. 2012, Oxford: Oxford University Press
Nicolas PJ: Statistics for Epidemiology. 2003, Danvers: Chapman and Hall/CRC
Lambert P-W: Zero-inflated Poisson regression, with an application to defects in Manufacturing. Technometrics. 1992, 34: 1-14. 10.2307/1269547.
Lingenfelser A, Rydzanicz K, Kaiser A, Becker N: Mosquito fauna and perspectives for integrated control of urban vector-mosquito populations in southern Benin (West Africa). Ann Agric Environ Med. 2010, 17: 49-57.
Onyabe Y, Conn E: Population genetic structure of the malaria mosquito Anopheles arabiensis across Nigeria suggests range expansion. Mol Ecol. 2003, 10: 2577-2591.
Wondji C, Simard M, Fontenille D: Evidence for genetic differentiation between the molecular forms M and S within the Forest chromosomal form of Anopheles gambiae in an area of sympatry. Insect Mol Biol. 2002, 11: 11-19. 10.1046/j.0962-1075.2001.00306.x.
Dabiré R, Diabaté A, Namountougou M, Toé K, Ouari A, Kengne P, Bass C, Baldet T: Distribution of pyrethroid and DDT resistance and the L1014F kdr mutation in Anopheles gambiae s.l. from Burkina Faso (West Africa). Trans R Soc Trop Med Hyg. 2009, 103: 1113-1120. 10.1016/j.trstmh.2009.01.008.
Akogbeto M: Etude des aspects epidemiologiques du paludisme côtier lagunaire au Bénin. PhD thesis. 1992, Paris: Université de paris XI
Akogbeto M: Entomological study on the malaria transmission in coastal and lagoon areas: the case of a village built on a brackish lake. Ann Société Belge Médecine Trop. 1995, 75: 219-227.
FAO: The State of Food Insecurity in the World. 2002, Roma: Food and Agriculture Organization
Baudoin M-A, Sanchez AC, Fandohan B: Small scale farmers’ vulnerability to climatic changes in southern Benin: the importance of farmers’ perceptions of existing institutions. Mitig Adapt Strateg Glob Change. 2013, 9: 1-13.
Yadouleton A, Asidi A, Djouaka R, Braïma J, Agossou C, Akogbeto M: Development of vegetable farming: a cause of the emergence of insecticide resistance in populations of Anopheles gambiae in urban areas of Benin. Malar J. 2009, 8: 103-10.1186/1475-2875-8-103.
Yadouleton A, Martin T, Padonou G, Chandre F, Asidi A, Djogbenou L, Dabiré R, Aïkpon R, Boko M, Glitho I, Akogbeto M: Cotton pest management practices and the selection of pyrethroid resistance in Anopheles gambiae population in Northern Benin. Parasit Vectors. 2011, 4: 60-10.1186/1756-3305-4-60.
Fontenille D, Cohuet A, Awono-Ambene P, Kengne P, Antonio-Nkondjio C, Wondji C, Simard F: Vecteurs de paludisme: du terrain à la génétique moléculaire. Rev Epidémiologie Santé Publique. 2005, 53: 283-290.
Koudou BG, Adja AM, Matthys B, Doumbia M, Cisse G, Kone M, Tanner M, Utzinger J: Pratiques agricoles et transmission du paludisme dans deux zones éco-épidémiologiques au centre de la Côte d’Ivoire. Bull Société Pathol Exot. 2007, 100: 124-126.
Yadouléton A, N’guessan R, Allagbé H, Asidi A, Boko M, Osse R, Padonou G, Kindé G, Akogbéto M: The impact of the expansion of urban vegetable farming on malaria transmission in major cities of Benin. Parasit Vectors. 2010, 3: 118-10.1186/1756-3305-3-118.
Omukunda E, Githeko A, Ndong’a MF, Mushinzimana E, Atieli H, Wamae P: Malaria vector population dynamics in highland and lowland regions of western Kenya. J Vector Borne Dis. 2013, 50: 85-92.
Terlouw DJ, Morgah K, Wolkon A, Dare A, Dorkenoo A, Eliades MJ, Eng JV, Sodahlon YK, Ter Kuile FO, Hawley WA: Impact of mass distribution of free long-lasting insecticidal nets on childhood malaria morbidity: the Togo national integrated child health campaign. Malar J. 2010, 9: 199-10.1186/1475-2875-9-199.
We thank the President’s Malaria Initiative via the Bureau for Global Health, US Agency for International Development, Benin National Malaria Control Programme, and the mosquito collectors.
The authors have declared that they have no competing interests.
VG designed experiments, coordinated field activities, wrote and revised the paper; RG, FRA, RO, RA, NA, RO, and EO participated in field activities, data collection and revised the paper; FO helped with statistical analysis and revised the paper; GGP helped in the study design and revised the paper; MCA designed the study, supervised field activities and revised the manuscript. All authors have read and approved the content of the manuscript.
Authors’ original submitted files for images
Below are the links to the authors’ original submitted files for images.