Residual malaria transmission and the role of Anopheles arabiensis and Anopheles melas in central Senegal.

Understanding the behavior and ecology of local malaria vectors is essential for the effectiveness of the commonly used vector-targeted malaria control tools in areas of low malaria transmission. This study was conducted to determine species composition, biting behavior and infectivity of the major Anopheles vectors of Plasmodium falciparum in low transmission settings in central Senegal. Adult mosquitoes were collected using human landing catches during 2 consecutive nights and Pyrethrum Spray Catches in 30-40 randomly selected rooms, from July 2017 to December 2018 in 3 villages. Anopheline mosquitoes were morphologically identified using conventional keys; their reproductive status assessed by ovary dissections, and a sub-sample of Anopheles gambiae s.l. were identified to species level using polymerase chain reaction (PCR). Plasmodium sporozoite infections were detected using real-time quantitative PCR. During this study 3684 Anopheles were collected of which 97% were An. gambiae s.l., 0.6% were Anopheles funestus, and 2.4% were Anopheles pharoensis. Molecular identification of 1,877 An. gambiae s.l. revealed a predominance of Anopheles arabiensis (68.7%), followed by Anopheles melas (28.8%), and Anopheles coluzzii (2.1%). The overall human-biting rate of An. gambiae s.l. was highest in the inland site of Keur Martin with 4.92 bites per person per night, while it was similar in the deltaic site, Diofior (0.51) and the coastal site, Mbine Coly (0.67). Parity rates were similar in An. arabiensis (45%) and An. melas (42%). Sporozoite infections were detected in both An. arabiensis and An. melas with the respective infection rates of 1.39% (N = 8) and 0.41% (N = 1). Results suggest that low residual malaria in central Senegal is transmitted by An. arabiensis and An. melas. Consequently, both vectors will need to be targeted as part of malaria elimination efforts in this area of Senegal.

Dynamic of resistance alleles of two major insecticide targets in Anopheles gambiae (s.l.) populations from Benin, West Africa.

BACKGROUND: Insecticide resistance is a growing concern for malaria control and vector control effectiveness relies on assessing it distribution and understanding its evolution. METHODS: We assessed resistance levels and the frequencies of two major target-site mutations, L1014F-VGSC and G119S-ace-1, conferring resistance to pyrethroids (PYRs) and carbamates/organophosphates (CXs/OPs) insecticides. These data were compared to those acquired between 2006 and 2010 to follow resistance evolutionary trends over ten years. RESULTS: We report the results of a 3-year survey (2013-2015) of insecticide resistance in 13 localities across the whole country of Benin. Permethrin (PYR) resistance was found in all populations tested, L1014F-VGSC being almost fixed everywhere, while bendiocarb resistance was limited to a few localities, G119S-ace-1 remaining rare, with very limited variations during surveyed period. Interestingly, we found no effect of the type of insecticide pressure on the dynamics of these mutations. CONCLUSIONS: These results confirm both the high prevalence of PYR resistance and the potential of CXs/OPs as short- to medium-term alternatives in Benin. They also underline the need for regular resistance monitoring and informed management in their usage, as the G119S-ace-1 mutation is already present in Benin and surrounding countries. Their unwise usage would rapidly lead to its spread, which would jeopardize PYR-resistant Anopheles control.

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.

Adaptive deletion in resistance gene duplications in the malaria vector Anopheles gambiae.

While gene copy-number variations play major roles in long-term evolution, their early dynamics remains largely unknown. However, examples of their role in short-term adaptation are accumulating: identical repetitions of a locus (homogeneous duplications) can provide a quantitative advantage, while the association of differing alleles (heterogeneous duplications) allows carrying two functions simultaneously. Such duplications often result from rearrangements of sometimes relatively large chromosome fragments, and even when adaptive, they can be associated with deleterious side effects that should, however, be reduced by subsequent evolution. Here, we took advantage of the unique model provided by the malaria mosquito Anopheles gambiae s.l. to investigate the early evolution of several duplications, heterogeneous and homogeneous, segregating in natural populations from West Africa. These duplications encompass ~200 kb and 11 genes, including the adaptive insecticide resistance ace-1 locus. Through the survey of several populations from three countries over 3-4 years, we showed that an internal deletion of all coamplified genes except ace-1 is currently spreading in West Africa and introgressing from An. gambiae s.s. to An. coluzzii. Both observations provide evidences of its selection, most likely due to reducing the gene-dosage disturbances caused by the excessive copies of the nonadaptive genes. Our study thus provides a unique example of the early adaptive trajectory of duplications and underlines the role of the environmental conditions (insecticide treatment practices and species ecology). It also emphasizes the striking diversity of adaptive responses in these mosquitoes and reveals a worrisome process of resistance/cost trade-off evolution that could impact the control of malaria vectors in Africa.

The evolutionary fate of heterogeneous gene duplications: A precarious overdominant equilibrium between environment, sublethality and complementation.

Gene duplications occur at a high rate. Although most appear detrimental, some homogeneous duplications (identical gene copies) can be selected for beneficial increase in produced proteins. Heterogeneous duplications, which combine divergent alleles of a single locus, are seldom studied due to the paucity of empirical data. We investigated their role in an ongoing adaptive process at the ace-1 locus in Culex pipiens mosquitoes. We assessed the worldwide diversity of the ace-1 alleles (single-copy, susceptible S and insecticide-resistant R, and duplicated D that pair one S and one R copy), analysed their phylogeography and measured their fitness to understand their early dynamics using population genetics models. It provides a coherent and comprehensive evolutionary scenario. We show that D alleles are present in most resistant populations and display a higher diversity than R alleles (27 vs. 4). Most appear to result from independent unequal crossing-overs between local single-copy alleles, suggesting a recurrent process. Most duplicated alleles have a limited geographic distribution, probably resulting from their homozygous sublethality (HS phenotype). In addition, heterozygotes carrying different HS D alleles showed complementation, indicating different recessive lethal mutations. Due to mosaic insecticide control practices, balancing selection (overdominance) plays a key role in the early dynamics of heterogeneous duplicated alleles; it also favours a high local polymorphism of HS D alleles in natural populations (overdominance reinforced by complementation). Overall, our study shows that the evolutionary fate of heterogeneous duplications (and their long-term role) depends on finely balanced selective pressures due to the environment and to their genomic structure.

The ace-1 Locus Is Amplified in All Resistant Anopheles gambiae Mosquitoes: Fitness Consequences of Homogeneous and Heterogeneous Duplications.

Gene copy-number variations are widespread in natural populations, but investigating their phenotypic consequences requires contemporary duplications under selection. Such duplications have been found at the ace-1 locus (encoding the organophosphate and carbamate insecticides' target) in the mosquito Anopheles gambiae (the major malaria vector); recent studies have revealed their intriguing complexity, consistent with the involvement of various numbers and types (susceptible or resistant to insecticide) of copies. We used an integrative approach, from genome to phenotype level, to investigate the influence of duplication architecture and gene-dosage on mosquito fitness. We found that both heterogeneous (i.e., one susceptible and one resistant ace-1 copy) and homogeneous (i.e., identical resistant copies) duplications segregated in field populations. The number of copies in homogeneous duplications was variable and positively correlated with acetylcholinesterase activity and resistance level. Determining the genomic structure of the duplicated region revealed that, in both types of duplication, ace-1 and 11 other genes formed tandem 203kb amplicons. We developed a diagnostic test for duplications, which showed that ace-1 was amplified in all 173 resistant mosquitoes analyzed (field-collected in several African countries), in heterogeneous or homogeneous duplications. Each type was associated with different fitness trade-offs: heterogeneous duplications conferred an intermediate phenotype (lower resistance and fitness costs), whereas homogeneous duplications tended to increase both resistance and fitness cost, in a complex manner. The type of duplication selected seemed thus to depend on the intensity and distribution of selection pressures. This versatility of trade-offs available through gene duplication highlights the importance of large mutation events in adaptation to environmental variation. This impressive adaptability could have a major impact on vector control in Africa.

An ace-1 gene duplication resorbs the fitness cost associated with resistance in Anopheles gambiae, the main malaria mosquito.

Widespread resistance to pyrethroids threatens malaria control in Africa. Consequently, several countries switched to carbamates and organophophates insecticides for indoor residual spraying. However, a mutation in the ace-1 gene conferring resistance to these compounds (ace-1(R) allele), is already present. Furthermore, a duplicated allele (ace-1(D)) recently appeared; characterizing its selective advantage is mandatory to evaluate the threat. Our data revealed that a unique duplication event, pairing a susceptible and a resistant copy of the ace-1 gene spread through West Africa. Further investigations revealed that, while ace-1(D) confers less resistance than ace-1(R), the high fitness cost associated with ace-1(R) is almost completely suppressed by the duplication for all traits studied. ace-1 duplication thus represents a permanent heterozygote phenotype, selected, and thus spreading, due to the mosaic nature of mosquito control. It provides malaria mosquito with a new evolutionary path that could hamper resistance management.

Phenotypic effects of concomitant insensitive acetylcholinesterase (ace-1(R)) and knockdown resistance (kdr(R)) in Anopheles gambiae: a hindrance for insecticide resistance management for malaria vector control.

BACKGROUND: Malaria is endemic in sub-Saharan Africa with considerable burden for human health. Major insecticide resistance mechanisms such as kdr(R) and ace-1(R) alleles constitute a hindrance to malaria vector control programs. Anopheles gambiae bearing both kdr and ace-1 resistant alleles are increasingly recorded in wild populations. In order to maintain the efficacy of vector control strategies, the characterization of concomitant kdr and ace-1 resistance, and their pleiotropic effects on malaria vector phenotype on insecticide efficacy are important. METHODS: Larval and adult bioassays were performed with different insecticide classes used in public health following WHO standard guidelines on four laboratory Anopheles gambiae strains, sharing the same genetic background but harboring distinct resistance status: KISUMU with no resistance allele; ACERKIS with ace-1(R) allele; KISKDR with kdr(R) allele and ACERKDRKIS with both resistance alleles' ace-1(R) and kdr(R) . RESULTS: Larval bioassays indicate that the homozygote resistant strain harboring both alleles (ACERKDRKIS) displayed slightly but significantly higher resistance level to various insecticides like carbamates (bendiocarb, p < 0.001; propoxur, p = 0.02) and organophosphates (chlorpyriphos-methyl, p = 0.002; fenitrothion, p < 0.001) when compared to ACERKIS strain. However, no differences were recorded between ACERKDRKIS and KISKDR resistance level against permethrin (Pyrethroid, p = 0.7) and DDT (Organochlorine, p = 0.24). For adult bioassays, the percentages of mosquitoes knocked down were significantly lower for ACERKDRKIS than for KISKDR with permethrin (p = 0.003) but not with deltamethrin. The percentage of mortality from adult bioassays was similar between ACERKDRKIS and ACERKIS with carbamates and organophosphates, or between ACERKDRKIS and KISKDR with pyrethroid and DDT. Concerning acetylcholinesterase enzyme, ACERKDRKIS strain showed similarAChE1 activity than that of ACERKIS. CONCLUSION: The presence of both kdr(R) and ace-1(R) alleles seems to increase the resistance levels to both carbamate and organophosphate insecticides and at operational level, may represent an important threat to malaria vector control programs in West Africa.

Characterization of swarming and mating behaviour between Anopheles coluzzii and Anopheles melas in a sympatry area of Benin.

The swarm structure of two sibling species, Anopheles gambiae coluzzii and Anopheles melas, was characterize to explore the ecological and environmental parameters associated with the formation of swarms and their spatial distribution. Swarms and breeding sites were searched and sampled between January and December 2010, and larval and adult samples were identified by PCR. During the dry season, 456 swarms of An. gambiae s.l. were sampled from 38 swarm sites yielding a total of 23,274 males and 76 females. Of these 38 swarming sites, 18 were composed exclusively of An. gambiae coluzzii and 20 exclusively of An. melas, presenting clear evidence of reproductive swarm segregation. The species makeup of couples sampled from swarms also demonstrated assortative mating. The swarm site localization was close to human dwellings in the case of the An. gambiae coluzzii and on salt production sites for An. melas. At the peak of the rainy season, swarms of An. melas were absent. These findings offer evidence that the ecological speciation of these two sibling species of mosquitoes is associated with spatial swarm segregation and assortative mating, providing strong support for the hypothesis that mate recognition is currently maintaining adaptive differentiation and promoting ecological speciation. Further studies on the swarming and mating systems of An. gambiae, with the prospect of producing a predictive model of swarm distribution, are needed to inform any future efforts to implement strategies based on the use of GMM or SIT.

Effect of three larval diets on larval development and male sexual performance of Anopheles gambiae s.s.

Population replacement/elimination strategies based on mass-release of sterile or otherwise genetically modified (male) mosquitoes are being considered in order to expand the malaria vector control arsenal on the way to eradication. A challenge in this context, is to produce male mosquitoes that will be able to compete and mate with wild females more efficiently than their wild counterparts, i.e. high fitness males. This study explored the effect of three larval food diets developed by the International Atomic Energy Agency on the overall fitness and mating performance of male Anopheles gambiae s.s. mosquitoes (Kisumu strain). Larval development (pupation and emergence rate, development time) was monitored, and adult wing length and energy reserves at emergence (i.e. lipids, sugars, glycogen and proteins) were measured. Male sexual performance was assessed through an insemination test whereby one male and 10 virgin females were maintained together in the same cage in order to record the number of inseminated females per 24h. Our results show that males reared on Diets 2 and 3 performed best during larval development. Males provided with treatment 2.2 had a shorter development time and performed best in insemination tests. However, these males had the lowest overall lifespan, suggesting a trade-off between longevity and sexual performances which needs to be taken into consideration when planning release. The results from this work were discussed in the context of sterile insect techniques or genetic control methods which is today one of the strategy in the overall mosquito control and elimination efforts.