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 L1014F kdr 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.