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

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.

The current distribution and characterization of the L1014F resistance allele of the kdr gene in three malaria vectors (Anopheles gambiae, Anopheles coluzzii, Anopheles arabiensis) in Benin (West Africa)

Background

The fight against malaria faces various biological obstacles, including the resistance of parasites to anti-malarial drugs and the resistance of mosquito vectors to insecticides. The resistance of Anopheles gambiae sensu lato (s.l.) to pyrethroids, the only class of insecticides used to impregnate mosquito nets, is known in Benin; the expansion of this resistance is influenced by the existence of gene flow between species, otherwise by the presence or absence of the kdr mutation in them. The objective of this study is to determine the spatial distribution of An. gambiae and the level of expression of the pyrethroid resistance kdr gene in seven agro-ecological zones of Benin.

Methods

The study was conducted in 18 localities belonging to seven agro-ecological zones where environmental parameters varied. The sites represent the main areas of eco-epidemiological malaria in Benin. Anopheles gambiae larvae were collected in natural breeding sites using ladles and dipping method and reared under standard conditions. These larvae were reared under standard conditions of temperature and humidity (26 to 30 °C and 60 to 90%) at the insectarium of the Centre de Recherche Entomologique de Cotonou (CREC). Adult female mosquitoes having emerged are morphologically and molecularly identified. Homozygous resistant (1014F/1014F), homozygous sensitive (1014L/1014L) and heterozygous (1014F/1014L) genotypes of the L1014Fkdr gene mutation are determined by PCR.

Results

A total of 677 An. gambiae was subjected at the PCR. The results revealed the presence of three vector species of the An. gambiae complex, of which 409 Anopheles coluzzii, 259 An.gambiae, 5 hybrids (An. coluzzii/An. gambiae) and 4 Anopheles arabiensis in the different agro-ecological zones. The four An. arabiensis were only found in Dassa, a locality in the cotton zone of central Benin. The frequency of distribution of the L1014F allele of the kdr gene varies from 84.48 to 100% in An. gambiae, from 80 to 100% in An. coluzzii and from 0 to 75% in An. arabiensis in the different agro-ecological zones. Moreover, a significant difference is generally observed in the distribution of the L1014F allele (P < 0.05). By comparing in pairs the distribution frequencies of this allele in the two species by agro-ecological zone, only a significant difference is noted in the central cotton and fishery zones (P = 0.0496).

Conclusion

In summary, even if the data are in small portions, the An. Arabiensis species was found only in central Benin and the L1014F allele of the kdr gene is widespread and seems to fix in all the species recorded in the different agro-ecological zones. This situation amplifies the problem of resistance, which could eventually be a significant obstacle for the malaria vectors control. Similarly, a study of their genetic structure via the L1014F allele is necessary in order to put in place strategies to manage this resistance. These strategies will take into account both the ecology and the genetic diversity of the organisms involved to preserve the effectiveness of pyrethroids, the only insecticides used for the impregnation of mosquito nets.

Genetic structure of Anopheles gambiae s.s populations following the use of insecticides on several consecutive years in southern Benin

Background

Several studies have reported the strong resistance of Anopheles gambiae s.l. complex species to pyrethroids. The voltage-dependent sodium channel (Vgsc) gene is the main target of pyrethroids and DDT. In Benin, the frequency of the resistant allele (L1014F) of this gene varies along the north-south transect. Monitoring the evolution of resistance is necessary to better appreciate the genetic structure of vector populations in localities subject to the intensive use of chemicals associated with other control initiatives. The purpose of this study was to map the distribution of pyrethroid insecticide resistance alleles of the Kdr gene in malaria vectors in different regions and ecological facies in order to identify the evolutionary forces that might be the basis of anopheline population dynamics.

Methods

The characterization of Anopheles gambiae s.l. populations and resistance mechanisms were performed using adult mosquitoes obtained from larvae collected in the four agroecological zones in southern Benin. Genomic DNA extraction was performed on whole mosquitoes.The extracted genomic DNA from them were used for the molecular identification of species in Anopheles gambiae s.l. complex and the identification of genotypes related to pyrethroid resistance as the Kdr gene amino acid position 1014 in sodium channel. Molecular speciation and genotyping of Kdr resistant alleles (1014) were done using PCR.Genepop software version 4.2 was used to calculate allelic and genotypic frequencies in each agroecological zone. The p value of the allelic frequency was determined using the binomial test function in R version 3.3.3. The Hardy-Weinberg equilibrium was checked for each population with Genetics software version 1.3.8.1. The observed heterozygosity and the expected heterozygosity as well as the fixation index and genetic differentiation index within and between populations were calculated using Genepop software version 4.2.

Results

During the study period, Anopheles coluzzii was the major species in all agroecological zones while Anopheles gambiae was scarcely represented. Regardless of the species, resistant homozygote individuals (L1014F/L1014F) were dominant in all agroecological zones, showing a strong selection of the resistant allele (L1014F). All populations showed a deficit of heterozygosity. No genetic differentiation was observed between the different populations of the two species. For Anopheles coluzzii, there was a small differentiation among the populations of the central cotton and bar-lands zones. The genetic differentiation was modest among the population of the fisheries zone (Fst = 0.1295). The genetic differentiation was very high in the population of Anopheles gambiae of the bar-lands zone (Fst = 0.2408).

Conclusion

This study revealed that the use of insecticides in Benin for years has altered the genetic structure of Anopheles gambiae s.s. populations in all agroecological zones of southern Benin. It would be desirable to orientate vector control efforts towards the use of insecticides other than pyrethroids and DDT or combinations of insecticides with different modes of action.