Malaria prevalence and transmission in the Zakpota sub-district of central Benin: baseline characteristics for a community randomised trial of a new insecticide for indoor residual spraying

TL;DR

This study provides baseline data on malaria prevalence and transmission in Benin to support a trial of a new insecticide for indoor spraying.

Contribution

The paper presents baseline data for a community trial of VECTRON™ T500, a new broflanilide-based indoor residual spraying product.

Findings

Malaria prevalence among children under 5 was 19% in the Zakpota sub-district.

Vectors showed high pyrethroid resistance but were largely susceptible to broflanilide and clothianidin.

Transmission risk varied greatly across villages, with an average entomological inoculation rate of 16.1 infected bites per person per year.

Abstract

Malaria transmission is known to be perennial and heterogeneous in Benin. Studies assessing local malaria prevalence, transmission levels and vector characteristics are critical for designing, monitoring and evaluating new vector control interventions in community trials. We conducted a study in the Zakpota sub-district of central Benin to collect baseline data on household characteristics, malaria prevalence, vector characteristics and transmission dynamics in preparation for a randomised controlled trial to evaluate the community impact of VECTRON™ T500, a new broflanilide indoor residual spraying (IRS) product.

A total of 480 children under 5 years of age from the 15 villages of the sub-district were tested for malaria by rapid diagnostic tests (RDTs). Mosquitoes were collected by human landing catches (HLCs), pyrethrum spray catches (PSCs) and Centers for Disease Control and Prevention miniature light traps (CDC-LTs) in selected houses in each village to assess vector density, composition, vector infectivity and prevalence of insecticide resistance markers. Bioassays were performed to detect vector susceptibility to pyrethroids, broflanilide (6 µg/bottle) and clothianidin (90 µg/bottle).

A total of 9080 households were enumerated in the 15 study villages. Insecticide-treated net (ITN) usage was > 90%, with 1–2 ITNs owned per household. Houses were constructed mainly with cement (44%) and mud (38%) substrates or a mixture of cement and mud (18%), and 60% of them had open eaves. The overall prevalence of P. falciparum infection was 19% among surveyed children: 20% among females and 18% among males. The haemoglobin rate showed an anaemia (< 11 g/dl) prevalence of 66%. Anopheles coluzzii and An. gambiae sensu stricto (s.s.) were the two vector species present at an overall proportion of 46% versus 54%, respectively. The human biting rate was 2.3 bites per person per night (b/p/n) and biting occurred mostly indoors compared with outdoors (IRR = 0.776; P = 0.001). The overall proportion of outdoor biting was 44% and exceeded indoor biting in three villages. The sporozoite rate was 2% with a combined yearly entomological inoculation rate (EIR) of 16.1 infected bites per person per year (ib/p/y). There was great variability in malaria transmission risk across the villages, with EIR ranging from 0 to 29.3 ib/p/y. The vector population showed a high intensity of resistance to pyrethroids across the study villages but was largely susceptible to broflanilide and clothianidin.

This study found high levels of malaria prevalence, vector density and transmission in the Zakpota sub-district despite the wide use of insecticide-treated nets. The vector population was mostly indoor resting and showed a high intensity of pyrethroid resistance but was generally fully susceptible to broflanilide. These findings demonstrated the suitability of the study area for the assessment of VECTRON™ T500 in a community randomised trial.

The online version contains supplementary material available at 10.1186/s13071-024-06342-1.

Figures

Spatial distribution of the 15 villages in the sub-district of Zakpota

Overall mosquito species composition in the study area per sampling methodFig. 3*An. gambiae* s.l., *Aedes* spp. and *Culex* spp. mosquitoes hourly indoor (a) and outdoor (b) biting rates in the study area

Indoor and outdoor density of *An. gambiae s.l.* in the study villages of the sub-district of Zakpota Centre

Mortality (24 h) of *An. gambiae* s.l. from eight villages exposed to deltamethrin at 1×, 2×, 5× and 10× the discriminating concentration of 12.5 µg/bottle in bottle bioassaysFig. 7Mortality (72 h) of *An. gambiae* s.l. exposed to broflanilide (6 µg/bottle) and clothianidin (90 µg/bottle) in bottle bioassays

Full text

Background

Despite considerable worldwide efforts deployed in recent years to control malaria, the disease remains a major public health problem in many endemic countries [1, 2]. Malaria cases continue to rise globally, reaching 249 million in 2022, with the African Region shouldering about 95% of cases and 96% of deaths [1]. As in most sub-Saharan countries, malaria remains a serious public health challenge in the Republic of Benin, accounting for 50% of all medical consultations, 40% of hospital admissions in health facilities and 20% of deaths among children under 5 years in 2021 [3].

While malaria transmission is known to be perennial in Benin, there are reports of heterogeneity between regions [4–6]. A cross-sectional epidemiological and clinical survey involving  > 10,000 children from 72 villages in two health districts showed that the prevalence of Plasmodium falciparum infection was moderate in the south and high in the north, with more heterogeneity in the southern parts of the country [7]. A recent transmission intensity stratification study using entomological inoculation rates (EIRs) also demonstrated great variability across districts and regions [8]. While mosquitoes of the Anopheles gambiae complex (mostly An. gambiae and An. coluzzii) play a major role in the transmission of P. falciparum from the north-to-south transect of the country [8, 9], An. funestus [10–12] and An. nili have also been implicated as secondary vectors [11, 13]. Anopheles gambiae and An. coluzzii exhibit broadly similar anthropophagic and indoor biting behaviour in Benin, with the latter showing a predominant presence during the dry seasons while the former is more abundant in wet seasons [8, 9, 14, 15]. Pyrethroid resistance, mediated by the L1014F mutation and overexpressed P450 enzymes, is pervasive in Anopheles vector mosquitoes in Benin and of moderate to high intensity across the different ecological zones of the country [16–18]. Resistance to organophosphate and carbamate insecticides has been reported in some areas in Benin [19, 20] but is less prevalent than pyrethroid resistance [15, 17, 18, 21–24].

Malaria prevention and control in Benin is largely dependent on the large-scale deployment of insecticide-treated nets (ITNs) via mass campaigns. However, studies have reported a persistent high malaria prevalence in many regions of Benin despite high ITN coverage and use [4, 25] and this has been attributed to high levels of pyrethroid resistance in local vectors and poor durability of ITNs. The identification of novel active ingredients (AIs) that can be used on ITNs and for indoor residual spraying (IRS) for effective control of pyrethroid-resistant malaria vectors is a priority for the vector control community [26, 27]. This requires cluster randomised trials to investigate their impact under community use in different epidemiological settings. Assessing local malaria prevalence, transmission levels and vector characteristics is critical for designing, monitoring and evaluating new vector control interventions [28, 29] in community trials given the heterogeneity in transmission and high variability in vector abundance and composition in most endemic settings [30]. Such data might also be useful for decision-making for the deployment of appropriate context-specific vector control interventions. The present study was conducted in the Zakpota sub-district of central Benin to collect baseline data on malaria prevalence, vector characteristics and transmission dynamics in preparation for a randomised controlled trial to evaluate the community impact of the newly developed VECTRON™ T500, a novel broflanilide insecticide product for IRS [31]. The AI is a meta-diamide that acts as a non-competitive antagonist of the γ-aminobutyric acid (GABA) receptor of chloride channels of the insect inhibitory nervous system thus presenting a new mode of action for malaria vector control [32]. The randomised controlled trial evaluated the efficacy and residual activity of VECTRON™ T500 compared to Fludora^®^ Fusion, a World Health Organization prequalification (WHO/PQ)-listed wettable powder IRS formulation of a mixture of clothianidin and deltamethrin.

Methods

Study area

The study was conducted in the Zakpota Centre sub-district located in the department of Zou, southern Benin. The study area covers 600 km^2^ separated into 15 administrative villages (Fig. 1). Economic activities of the population include agriculture (85%), trade (8%), crafts (5%) and other services (2%) [33]. The Zakpota Centre sub-district has a sub-equatorial climate with two rainy seasons (April–July and September–November) and two dry seasons (December–March and July–August). The average annual rainfall of Zakpota is 980 mm, with important inter-annual fluctuations over the past 40 years, a source of climatic uncertainty for this area. Malaria prevention in Zakpota mainly relies on the distribution of ITNs through mass campaigns led by the Ministry of Health. This is supplemented by routine distributions to pregnant women and children under 5 years through antenatal visits and the expanded programme of immunisation.Fig. 1. Spatial distribution of the 15 villages in the sub-district of Zakpota

Household survey

Detailed information on the distribution, location, number, type, size and accessibility of households was needed for the IRS campaign. Mapping and censuses were conducted in November 2020 to generate a database of all houses in the study area. Every building of each hamlet/village, village boundaries, health facilities and other landmarks were geo-referenced, mapped and recorded in the field using integrated handheld tablets fitted with global positioning system (GPS). Information was collected on the size of the household, house structure, the number of rooms, availability and type of bednets/ITNs, number of sleeping places, and wall substrate type in each village using a structured questionnaire.

Cross-sectional clinical survey

A cross-sectional clinical survey was conducted in May 2021 to obtain data on malaria prevalence and anaemia in children under 5 years of age in all 15 villages in the Zakpota sub-district, before the implementation of the IRS campaign. The survey team included experienced nurses and community health workers (CHWs). A sample size of 480 participants was found sufficient based on sample size calculations performed using a baseline malaria prevalence of 46% observed in a recent survey conducted in the district [34] and a coefficient of variation of 0.3. Using the census data, a random selection of 30 households with children under 5 years were visited in each of the villages where one child was enrolled in the study per household giving a total of 480 study participants. Participating children were tested for malaria using a rapid diagnostic test (RDT) (CareStart malaria HRP2/pLDH [Pf/pan] combo, DiaSys, UK) irrespective of symptoms. Each child's axillary temperature was assessed, and those who tested positive by RDT and showed a fever (> 37.5 °C) were considered symptomatic for malaria and were provided free treatment. The haemoglobin concentration (g/dl)xxx was measured using a Hemo Control diagnostic analyser (EKF Diagnostics, Cardiff, UK) to assess anaemia in all enrolled children. Children having a haemoglobin level below 11 g/dl were considered anaemic.

Entomological data collection

Baseline entomological studies were conducted every 2 months between August 2020 and April 2021 over four sampling rounds in the 15 villages of the Zakpota Centre sub-district. Data on the presence, distribution, behaviour, age structure, insecticide susceptibility status of vectors and baseline vector parameters were collected in each village of the study area. Mosquitoes were collected using three methods.

• Human landing catches (HLCs) were performed for two nights at each time point, in six randomly selected houses in each village. Four consenting human volunteer mosquito collectors were used per house, with volunteers working in pairs—one collecting mosquitoes indoors and the other outdoors. Collections were performed between 07:00 p.m. and 07:00 a.m., with pairs alternating every 6 h, i.e. 07:00 p.m. to 01:00 a.m. and 01:00 a.m. to 07:00 a.m. Collectors sat on chairs indoors and outdoors with their legs exposed; the outdoor collector was positioned at least 10 m from the house. Using flashlights, collectors caught mosquitoes landing on their legs with haemolysis tubes, and each hour’s collection was kept separately.
• Centers for Disease Control and Prevention miniature light traps (CDC-LTs) (Model 512, John W. Hock, FL, USA) were used to collect mosquitoes for two nights at every time point inside six randomly selected houses per village. The traps were set up daily at approximately 1–1.5 m above the ground and close to human-occupied untreated bednets. The light traps were operated overnight between 7:00 p.m. and 7:00 a.m. daily.
• Pyrethrum spray catches (PSCs) were also performed in six randomly selected houses in each village. On each day of sampling, all food and water were removed from the house before spraying and white sheets were laid on the floor. The house was sprayed with an aerosol (Rambo^@^) containing 0.25% transfluthrin and 0.20% permethrin. All mosquitoes falling onto the white sheets 15 min post-spraying were collected using forceps.

Processing of mosquito samples

Mosquitoes were identified to species level using morphological keys [35]. Anopheles gambiae sensu lato (s.l.) sibling species identification was performed using the polymerase chain reaction (PCR) method [36]. For all sampling methods, the head and thorax of each mosquito collected were separated from the abdomen and tested for the presence of P. falciparum circumsporozoite protein (CSP) using direct enzyme-linked immunosorbent assay (ELISA) [37]. Sporozoite rates (SRs) for each village were estimated as the proportion of P. falciparum-positive mosquitoes, while the EIRs were estimated as the number of infective bites per person per night. The abdomens of unfed mosquitoes from all three collections morphologically identified as Anopheles gambiae s.l. were dissected to assess parity using Detinova’s ovary tracheation method [38]. The parity rate was measured as the ratio of parous mosquitoes out of the total number dissected. A sub-sample of mosquitoes morphologically identified as An. gambiae s.l. was processed for the presence of the L1014F and G119S insecticide resistance mutations using the methods described by Martinez-Torres [39] and Weill [40], respectively.

Insecticide susceptibility monitoring

WHO bottle bioassays were used to assess the susceptibility of mosquitoes from each village to clothianidin and broflanilide, the AIs contained respectively in Fludora^®^ Fusion and VECTRON™ T500, the IRS products evaluated in the RCT after the present baseline data collections. Clothianidin was tested at a discriminating concentration of 90 µg/bottle provided by the manufacturer while broflanilide was tested at 6 µg/bottle determined in preliminary dose–response bioassays [41]. To prevent crystallization of the insecticide in bottle bioassays and to improve bioavailability and uptake by exposed mosquitoes, Mero^®^ (81% rapeseed oil methyl ester; Bayer CropScience) was used as an adjuvant at a dose of 500 parts per million (ppm) for broflanilide and 1000 ppm for clothianidin. The intensity of resistance to deltamethrin in each study arm was also assessed at baseline in bottles treated at 1×, 2×, 5× and 10× the discriminating concentration of 12.5 µg/bottle. One hundred unfed wild adult female first filial generation (F1) An. gambiae s.l. (2–5 days old) raised in insectary from larvae collected from breeding sites in study villages were exposed for 1 h in bottles treated with insecticide in batches of 25 mosquitoes per bottle. After exposure, the mosquitoes were gently aspirated from the bottle into clean, labelled paper cups and provided with 10% sugar solution-soaked cotton wool. Knockdown was recorded after exposure, and mortality was recorded 24 h after exposure to deltamethrin, 72 h after exposure to broflanilide and 120 h after exposure to clothianidin due to the different modes of action and speed of kill of these insecticides. Mosquitoes were also exposed in bottles treated with acetone and Mero alone to serve as controls.

Ethical considerations

This study was approved by the national ethics review committee of the Ministry of Health Benin (CNERS, No. 29 and 51) and the ethics review board of the London School of Hygiene & Tropical Medicine (LSHTM, Ref: 22453). Before any project activities, village and hamlet leaders were invited to sensitization sessions and written informed consent was obtained from them before starting data collection. Heads of households involved in the study gave informed consent before participating in the study. Written informed consent was obtained from parents of children enrolled in the clinical survey. Any children found positive for malaria were provided free anti-malaria treatment. Mosquito collectors were trained and they also provided written informed consent. All participating households and individuals were free to withdraw from the studies at any time.

Data management and analysis

The census and prevalence survey data were captured on electronic forms using smartphones installed with OpenDataKit (ODK) Collect. Entomological data was collected on pre-designed data collection forms and then double-entered into pre-established Microsoft Excel databases. All data were stored on a secure server located at the Centre de Recherche Entomologique de Cotonou (CREC)/LSHTM Facility in Benin and made available for statistical analysis. Descriptive statistics were used to summarise demographic data. Differences in malaria prevalence were analysed using mixed-effects logistic regression models, while differences in indoor and outdoor vector density were analysed using mixed-effects negative binomial regression. All analyses were performed with Stata version 17 (StataCorp LLC).

The resistance status of wild vector populations of An. gambiae s.l. was interpreted according to WHO guidelines [42, 43]: mortality ≥ 98% indicates susceptibility, mortality less than 90% indicates resistance, and mortality of 90–97% means a possible resistance that needs to be confirmed. The intensity of resistance was interpreted as low, moderate or high intensity when 24-h mortality rates were ≥ 98% at 5× the discriminating concentration, ≥ 98% at 10× the discriminating concentration or < 98% at 10× the discriminating concentration, respectively.

Results

Characteristics of study population and households

Summary data collected on demography and housing structure characteristics during the baseline census are presented in Tables 1 and 2, respectively. Overall, 34,976 people living in 9080 households were enumerated in the 15 study villages of which 48% were females. The mean number of persons per household was 3–5 with 1–2 ITNs owned per household. All ITNs found in the study area were treated only with pyrethroids and the predominant brands were Yorkool^®^ and PermaNet^®^ 2.0 distributed in the 2020 mass campaign by the Ministry of Health. The proportion of the population under 5 years old was 15%, while pregnant women constituted 2% of the population (Table 1). Table 1. Demographic data in the 15 administrative villages of the sub-district of Zakpota CentreGenderMalaria-vulnerable populationHouseholdMale%Female%Pregnant women%Children under 5 years%TotalMean number of people per HHMean number of sleeping placesMean number of sleeping placesMean number ITNsMean number ITNsAdjido1630531470474514371490331006110901Agbogbomey835518084948321813417475626262Agbokpa741536654722222716368459725592Detekpa733517074921119413398462126412Djoyitin1031529564842231916554488629232Dokpa618496495121220616263553325012Gnandokpo408533674728410614187430322892Hounkanme140853124247733427166584102629441Kemondji11735310594757337517473594928822Kodota1479511406496624321563851391211442Sogbelankou858527884833222914484362217462Sohounta898537924755331018429469327022Za-Kekere1775501750507525241585241574217032Za-Kpota Centre3737513601491732111215197843351229932Za-Zounme712516804926220915478363116561Total18,0365216,9404878525325159080414,939214,3992Table 2Housing structure characteristics in the 15 administrative villages of the Zakpota Centre sub-districtTotal% Eligible household (95%CI)Substrate of inner wallAspect of inner wallEave% Cement (95%CI)% Mud (95%CI)% Cement & Mud (95%CI)% Plastered  (95%CI)% Not plastered (95%CI)% Partially plastered (95%CI)% Close (95%CI)% Open (95%CI)Adjido121688 (86–90)37 (34–41)42 (39–45)21 (19–23)49 (46–52)40 (37–43)11 (9–13)29 (26–32)71 (68–74)Agbogbomey57890 (88–92)55 (51–59)26 (22–30)19 (16–22)57 (54–60)32 (28–36)11 (8–14)54 (50–58)46 (42–50)Agbokpa57086 (83–89)15 (12–18)79 (76–82)6 (4–8)20 (17–23)77 (74–80)3 (2–4)37 (33–41)63 (59–67)Detekpa55483 (80–86)29 (25–34)58 (54–62)13 (9–15)39 (35–43)54 (50–58)7 (5–9)19 (16–22)81 (78–84)Djoyitin71383 (80–86)39 (35–44)45 (41–49)16 (13–19)43 (39–47)54 (50–58)3 (2–4)39 (35–43)61 (57–65)Dokpa41892 (89–95)40 (35–45)49 (44–54)11 (8–14)52 (47–57)37 (32–42)11 (9–13)34 (29–39)66 (61–71)Gnandokpo25489 (85–93)33 (27–39)46 (40–52)21 (16–26)45 (39–51)51 (45–56)4 (2–6)14 (10–18)86 (82–90)Hounkanme97786 (84–88)50 (47–53)32 (29–35)18 (16–20)58 (55–61)37 (34–40)5 (4–7)46 (43–49)54 (51–57)Kemondji80886 (84–88)28 (25–31)59 (56–62)13 (10–14)34 (31–37)63 (60–66)3 (2–4)26 (23–29)74 (71–77)Kodota93593 (91–95)49 (46–52)41 (38–44)10 (9–13)68 (65–71)27 (24–30)5 (4–6)39 (36–42)61 (58–64)Sogbelankou67481 (78–84)29 (26–32)53 (49–57)18 (14–20)36 (32–40)61 (57–65)3 (1–4)30 (27–33)70 (67–73)Sohounta64693 (91–95)45 (41–49)36 (32–40)19 (15–21)54 (50–58)44 (40–48)2 (1–3)42 (38–46)58 (54–62)Za-Kekere139187 (85–89)55 (52–58)27 (25–29)17 (15–19)49 (46–52)47 (44–50)4 (3–5)51 (48–54)49 (46–52)Za-Kpota Centre232590 (89–91)58 (56–60)18 (16–20)24 (22–26)66 (64–68)23 (21–25)11 (10–12)50 (48–52)50 (48–52)Za-Zounme54989 (86–92)35 (30–38)36 (32–40)29 (27–31)47 (43–51)44 (40–48)9 (7–11)38 (34–42)62 (58–66)Total12,60888 (87–89)44 (43–45)38 (37–39)18 (17–19)51 (50–52)42 (41–43)7 (6–8)40 (39–41)60 (59–61)

A total of 12,608 house structures were enumerated across the 15 villages and 88% of them were eligible for IRS (range 81% to 93% per village) as defined by WHO guidelines [44]. The houses were constructed mainly with cement (44%) or mud (38%) substrates while some were constructed with a mixture of cement and mud (18%). Overall, 51% of surveyed structures were plastered and 60% had open eaves (Table 2).

Malaria prevalence results

A total of 480 children under the age of 5 living in the sub-district of Zakpota Centre were surveyed of which 51% were female (Table 3). The age structure was broadly similar between both genders with children < 2 years constituting 26% of the sample size. ITN usage was very high amongst the children tested (> 90%) and similar between genders. The proportion of children reporting fever within the last 48 h was 17% among both males and females (Table 4). Table 3. Characteristics and age structure of children under 5 years tested in the malaria prevalence survey in the sub-district of Zakpota CentreFemaleMaleTotalN% (95% CI)N% (95% CI)N% (95% CI)Number of children24451 (47–55)23649 (45–43)480100Age structure (year) 0–131 (0–2)10 (0–1)41 (0–2) 1–26113 (10–16)5812 (9–15)11925 (21–29) 2–36814 (11–17)6514 (11–17)13328 (24–32) 3–45010 (7–13)449 (6–12)9420 (16–24) 4–56213 (10–16)6814 (11–17)13027 (23–31)Slept under net (%)22090 (87–93)21290 (87–93)43290 (87–93)Fever within last 48 h4217 (14–20)3917 (14–20)8117 (14–20)Table 4. Axillary fever, haemoglobin rate and malaria prevalence results in children under 5 years in the sub-district of Zakpota CentreFemaleMaleTotalN% (95% CI)N% (95% CI)N% (95% CI)Number of children24451 (47–55)23649 (45–43)480100Axillary temperature < 37.5 °C23596 (94–98)22294 (91–97)45795 (93–97) ≥ 37.5°c94 (2–6)146 (3–9)235 (3–7)Haemoglobin rate < 11 g/dL16869 (63–75)15164 (58–70)31966 (62–70) ≥ 11 g/dL7631 (25–37)8536 (30–42)16134 (30–38)RDT results (P. falciparum) Positive5020 (15–25)4318 (13–23)9319 (15–23) Negative19480 (75–85)19382 (77–87)38781 (77–85)

Impact of age and ITN ownership on malaria prevalence

Table 5 presents the ratio of geometric means of sub-district levels of P. falciparum infection rates for age groups and net ownership status. In the households surveyed, the infection rate was significantly lower in children < 2 years than in children aged 2–5 years (3% vs 16%, OR: 1.839, 95% CI 1.027–3.294, P = 0.04). The number of ITNs in the households appeared to have no association with the malaria infection rate in children. Infection rates were lower in households owning at least two ITNs compared to households owning at least one ITN (2% vs. 18% OR: 2.471, 95% CI 1.191–5.125, P = 0.015). Table 5. Risk factors associated with P. falciparum infection in householdsMalaria infected householdRR* (95% CI)P-valueN% (95% CI)Total of households480100––Age < 2 years163 (1–5)1.924 (0.954–3.880)0.067 2–5 years7716 (13–19)Net ownership ≤ 1 LN72 (1–3)2.022 (1.004–4.072)0.05 ≥ 2 LNs8618 (15–21)^a^Ratio of geometric means of village level P. falciparum infection rates

Mosquito composition and biting times

A total of 9562 female mosquitoes were collected in the study area during the four sampling time points using HLC, PSC and CDC-LT sampling methods. Anopheles gambiae s.l. was the predominant species, representing 38% of the cumulative collections for all three methods. Comparing mosquito composition between collection methods, the HLC was more efficient in collecting An. gambiae s.l. while the PSC and CDC-LT methods were more efficient in collecting Culex spp. mosquitoes (Fig. 2). All three collection methods collected approximately equal proportions of Aedes spp. mosquitoes (~ 13%). Peak mosquito biting times were 1:00 a.m. to 5:00 a.m. for An. gambiae s.l., early evening (before 8 p.m.) for Aedes spp., and 11:00 p.m. to 5:00 a.m. for Culex spp. (Fig. 3). Biting during peak hours was higher indoors for An. gambiae s.l. and Aedes spp., and outdoors for Culex spp. (Fig. 3).Fig. 2. Overall mosquito species composition in the study area per sampling methodFig. 3An. gambiae s.l., Aedes spp. and Culex spp. mosquitoes hourly indoor (a) and outdoor (b) biting rates in the study area

Anopheles gambiae sibling species composition

A total of 3460 An. gambiae s.l. mosquitoes were collected and 95% of these were collected by HLCs. The CDC LT and PSC methods collected very small proportions of An. gambiae s.l. (5% for both methods). A total of 471 specimens of An. gambiae s.l. from the study area were analysed for sibling species identification. An. coluzzii and An. gambiae s.s. were the 2 species present at an overall proportion of 46% and 54% respectively. The proportions of each species fluctuated depending on the collection time points with An. gambiae s.s. being more abundant at most time points (52–64%) except in February 2021 (Fig. 4).Fig. 4An. gambiae sibling species composition

Vector density indoors and outdoors

A total of 3460 An. gambiae s.l. vector mosquitoes were collected during the study of which 58% were collected indoors and 42% outdoors. The aggregated indoor density of An. gambiae s.l. was significantly higher than outdoor (incidence rate ratio (IRR) = 0.776; p = 0.001). Vector densities were higher indoors compared to outdoors in most villages (12 out of 15) except in Djoyitin, Hounkanme and Sohounta where outdoor biting was higher (Fig. 5).Fig. 5. Indoor and outdoor density of An. gambiae s.l. in the study villages of the sub-district of Zakpota Centre

Human biting rate and infectivity of An. gambiae s.l.

Table 6 presents the human biting rate (HBR), SR and yearly EIR in An. gambiae s.l. across the study villages. The overall mean HBR of the study area was 2.3 bites per person per night (b/p/n), and the range between villages was 0.3 b/p/n (Agbokpa) to 5.8 b/p/n (Za-Kekere) (Table 6). Mosquito infection rates determined by CSP ELISA performed on a total of 2863 An. gambiae s.l. mosquitoes showed a mean SR of 2%, with a village range of 0% in Sohounta to 8.5% in Kodota. Combined data revealed an overall yearly EIR of 16.1 infected bites per person per year (ib/p/y) in the study area. The villages with the highest malaria transmission risk were Agbogbomey (29.3 b/p/y), Za-Zakere (27.5 ib/p/y), Kemondji (25.7 ib/p/y), Detekpa (23.5 ib/p/y) and Za-Zounme (22.6 ib/p/y). Table 6. Human biting rate, sporozoite rate and entomological inoculation rate of An. gambiae s.l. in the villages of the sub-district of ZakpotaN collected^a^Man.nightHBR (b/p/n)Total analysedSR (%)EIR (ib/p/y)Agbogbomey340963.53092.329.3Agbokpa31960.3267.79.1Detekpa354963.72291.723.5Djoyitin2771921.42301.79.2Dokpa180961.91683.020.4Gnandokpa40960.4407.511.4Hounkanmey38960.4370.00.0Kemondji312963.32772.225.7Kodota59960.6598.519.0Sogbelankou295963.12080.55.4Sohounta31960.3240.00.0Za-Kekere11091925.89191.327.5Za-Kpota Centre33960.3336.17.6Za-Zoume361963.83041.622.6Total346015362.328632.016.1^a^Collected by HLC in four sampling time pointsHBR human biting rate, SR sporozoite rate, EIR entomological inoculation rateman.night: number of mosquito collectors x number of collection nights

Parous rate in An. gambiae s.l.

Table 7 presents An. gambiae s.l. parous rates across study villages in the sub-district of Zakpota Centre. Of the 1145 specimens of An. gambiae s.l. dissected over the study period, 561 were parous representing 47%. The parous rate per village ranged from 17% (Sogbélankou) to 100% (Kodota). Table 7. Parous rate in An. gambiae s.l. in the villages of sub-district of Zakpota CentreVillageTotal dissectedN ParousParous rate (%)95% CIAdjido955624–88Agbogbome141735244–60Detekpa172482821–35Djoitin49224531–49Gnandokpa1386236–88Hounkanmey966736–98Kemondji68426250–74Kodota2121100–Sogbelankou127211711–23Sohounta19157961–97Za-Agbokpa1275830–86Za-Dokpa45245338–68Za-Kekere3021705650–62Za-Kpota Centre26155839–77Za-Zoume152845547–63Total11655614845–51

Vector susceptibility to insecticides

Vector mortality with the discriminating concentration of deltamethrin (12.5 µg/bottle) ranged between 1 and 68% (Fig. 6). Mortality increased at doses 2×, 5× and 10× the discriminating concentration of deltamethrin but did not exceed 98% at 10× in any village, indicating a high intensity of pyrethroid resistance in the study area. With broflanilide tested at 6 µg/bottle, pooled mortality rates across the study villages showed full susceptibility of the vector population in the sub-district to the insecticide (96% mortality). There was however some variability in mortality in broflanilide susceptibility bioassays in two villages (Sohounta and Za-Agbokpa) showing 83% and 88% mortality, respectively (Fig. 7). Sub-optimal mortality in these villages was probably greater due to variability in bioassays with mosquito sub-populations rather than resistance given that the broflanilide susceptibility assay method and the discriminating concentration for broflanilide had not been fully optimised at the time of the baseline study. Subsequent bottle bioassays conducted post-intervention showed full susceptibility to broflanilide in the vector population [45]. The results also showed no evidence of resistance to clothianidin tested at 90 µg/bottle; mortality rates were generally > 90% in all villages. No mortality was recorded in the controls.Fig. 6. Mortality (24 h) of An. gambiae s.l. from eight villages exposed to deltamethrin at 1×, 2×, 5× and 10× the discriminating concentration of 12.5 µg/bottle in bottle bioassaysFig. 7Mortality (72 h) of An. gambiae s.l. exposed to broflanilide (6 µg/bottle) and clothianidin (90 µg/bottle) in bottle bioassays

Insecticide resistance genotypes

The total number of An. gambiae s.l. mosquitoes analysed for target site insecticide resistance genes was 460 for the L1014F kdr mutation and 467 for the G119S ace-1 mutation. The frequency of L1014F kdr in the An. gambiae complex across the study area was 84%, while the frequency of the ace-1 mutation was 26% (26%) (Table 8). Table 8. Frequency of the L1014F kdr and G119S ace-1 mutations across the study areaMutationN TestedRRRSSSFrequency95% CIL1014F kdr46034777368481–87G119S ace-1467142122412622–30

Discussion

The purpose of this study was to collect baseline information on the prevalence and transmission of malaria in the Zakpota Centre sub-district of Benin in preparation for a community randomised trial to evaluate a new insecticide for IRS. First, we found a malaria prevalence of 19% among children under 5, the most vulnerable age group to malaria. This prevalence was lower than the 42% prevalence reported in this age group in a similar larger-scale study performed in the Zou region of Benin in preparation for a randomised controlled trial of dual-AI ITNs [34]. The difference in malaria prevalence between the two studies may be attributed to the timing of the surveys given that the present study was performed at the onset of the main rainy season (April to May) when cases can be low, while the latter was performed during the second rainy season (October to November).

The main vector found in the study area was An. gambiae s.l., consisting of both An. gambiae s.s. and An. coluzzii with the former available in higher overall proportions and at most time points except in the middle of the dry season. This finding confirms several studies that have demonstrated that both species commonly occur in sympatry in West Africa, though they may differ substantially in their ecological behaviour, insecticide resistance status and vectorial capacity [46–48]. Previous entomological surveys in Benin have also shown a wide variability in the proportion of both sibling species in terms of geographical spread and seasons with reports of higher densities of An. gambiae s.s. increasing more recently [8, 9, 15, 49]. Apart from An. gambiae s.l., no other malaria vector species were found in the study area probably due to the absence of suitable breeding sites for their development.

Vector biting rates were greater indoors compared to outdoors, thus demonstrating the suitability of indoor vector control interventions including IRS for malaria control in the study area. However, the high levels of outdoor biting observed (44%) are nonetheless of concern. This corroborates several studies that have reported high outdoor malaria vector biting rates in Benin [12, 50] and elsewhere [51]. This finding highlights the need for additional vector control interventions that target outdoor biting vectors. Vector biting was highest at night (1:00 a.m. to 5:00 a.m.) both indoors and outdoors which aligns with the well-documented nocturnal behaviour of anopheline mosquitoes [21, 52]. As expected, the biting rates of daytime Aedes mosquitoes collected in the study area were substantially lower during this period. The overall annual EIR in the study area across the baseline entomological surveys was 16.1 bites per person with 2% of mosquitoes found to be infected with malaria sporozoites. Modelling studies suggest that this level of EIR would lead to a higher malaria prevalence than the 19% observed in the cross-sectional survey [53]. EIR also varied widely across the villages in the study area which can be explained by a high heterogeneity in anopheline densities and their infectivity across the individual villages.

The high intensity of pyrethroid resistance in susceptibility bioassays in this study aligns with several previous studies performed in the Zou region [15, 21] and Benin at large [17, 20, 23, 24, 54, 55]. The genotyping studies showed high frequencies of the L1014F kdr gene (84%) which would have contributed to the levels of phenotypic resistance observed. The lower levels of the G119S ace-1 (26%) corroborate the relatively lower levels of organophosphate and carbamate resistance observed in the Zou region and in most areas of Benin [20, 56, 57]. By contrast, evidence for resistance to broflanilide and clothianidin was not found, thus demonstrating the suitability of the study area for a trial deploying both insecticides for IRS in the study area.

Conclusion

This study found high levels of malaria prevalence, vector density and transmission in the Zakpota sub-district, despite the high use of ITNs in the study area. The vector population was mostly indoor resting and showed a high intensity of pyrethroid resistance but was generally susceptible to broflanilide. These findings demonstrated the suitability of the study area for the assessment of a new broflanilide insecticide product (VECTRON™ T500) for IRS in a community randomised trial. However, the levels of outdoor biting observed were high, which demonstrates the need for suitable vector control strategies to target malaria transmission by outdoor biting mosquitoes.

Supplementary Information

Additional file 1.