The frequency of kdr and ace-1 alleles in Anopheles gambiae s.l. before and during indoor residual spraying (IRS) implementation and four years after IRS withdrawal in three districts in Atacora, Benin
Authors: Odjo Esdras Mahoutin , Impoinvil Daniel , Fassinou Arsène Jacques Y. H. , Padonou Gil Germain , Aïkpon Rock , Salako Albert Sourou , Sominahouin André Aimé , Adoha Constantin , Yovogan Boulais , Osse Razaki , Oussou Olivier , Tokponnon Filémon , Gnanguénon Virgile , Hassani Ahmed Saadani , Akogbeto Martin C. 
Published: March 2024

Abstract
Background
Indoor residual spraying (IRS) was first implemented in the Atacora department, Benin from 2011 to 2012 using bendiocarb (carbamate) followed by annual spraying with pirimiphos-methyl (organophosphate) from 2013 to 2018. Before and after IRS implementation in Atacora, standard pyrethroid insecticide-treated bed nets were the main method of vector control in the area. This study investigated the knockdown resistance (kdr) gene (L1014F) and the acetylcholinesterase (ace-1) gene (G119S), before and during IRS implementation, and 4-years after IRS withdrawal from Atacora. This was done to assess how changes in insecticide pressure from indoor residual spraying may have altered the genotypic resistance profile of Anopheles gambiae s.l.

Method
Identification of sibling species of An. gambiae s.l. and detection of the L1014F mutation in the kdr gene and G119S mutation in ace-1 genes was done using molecular analysis. Allelic and genotypic frequencies were calculated and compared with each other before and during IRS implementation and 4 years after IRS withdrawal. The Hardy–Weinberg equilibrium and genetic differentiation within and between populations were assessed.

Results
Prevalence of the L1014F mutation in all geographic An. gambiae s.l. (An. gambiae s.s., Anopheles. coluzzii, Anopheles. arabiensis, and hybrids of “An. gambiae s.s. and An. coluzzii”) populations increased from 69% before IRS to 87% and 90% during and after IRS. The G119S allele frequency during IRS (20%) was significantly higher than before IRS implementation (2%). Four years after IRS withdrawal, allele frequencies returned to similar levels as before IRS (3%). Four years after IRS withdrawal, the populations showed excess heterozygosity at the ace-1 gene and deficit heterozygosity at the kdr gene, whereas both genes had excess heterozygosity before and during IRS (FIS < 0). No genetic differentiation was observed within the populations.

Conclusions
This study shows that the withdrawal of IRS with bendiocarb and pirimiphos-methyl may have slowed down the selection of individual mosquitoes with ace-1 resistance alleles in contrast to populations of An. gambiae s.l. with the L1014F resistance allele of the kdr gene. This may suggest that withdrawing the use of carbamates or organophosphates from IRS or rotating alternative insecticides with different modes of action may slow the development of ace-1 insecticide-resistance mutations. The increase in the prevalence of the L1014F mutation of the kdr gene in the population, despite the cessation of IRS, could be explained by the growing use of pyrethroids and DDT in agriculture and for other domestic use. More observational studies in countries where carbamates or organophosphates are still being used as public health insecticides may provide additional insights into these associations.

Graphical Abstract