- Open Access
Ace-1 duplication in Anopheles gambiae: a challenge for malaria control
Malaria Journal volume 8, Article number: 70 (2009)
Insecticide resistance is a rapid and recent evolutionary phenomenon with serious economic and public health implications. In the mosquito Anopheles gambiae s.s., main vector of malaria, resistance to organophosphates and carbamates is mainly due to a single amino-acid substitution in acetylcholinesterase 1 (AChE1). This mutation entails a large fitness cost. However, a resistant duplicated allele of the gene encoding AChE1 (ace-1), potentially associated to a lower fitness cost, recently appeared in An. gambiae.
Using molecular phenotype data collected from natural populations from West Africa, the frequency of this duplicated allele was investigated by statistical inference. This method is based on the departure from Hardy-Weinberg phenotypic frequency equilibrium caused by the presence of this new allele.
The duplicated allele, Ag-ace-1D, reaches a frequency up to 0.65 in Ivory Coast and Burkina Faso, and is potentially present in Benin. A previous study showed that Ag-ace-1D, present in both M and S molecular forms in different West Africa countries, was generated by a single genetic event. This single origin and its present distribution suggest that this new allele is currently spreading.
The spread of this less costly resistance allele could represent a major threat to public health, as it may impede An. gambiae control strategies, and thus increases the risk of malaria outbreaks.
Since early 1950s, humans have controlled the populations of many agricultural or medical arthropod pests, mostly with chemical insecticides. After years of success, evolutionary adaptations to these new conditions began to occur and resistance spread rapidly; more than 500 species are now resistant to at least one insecticide . Insecticide resistance is a rapid and recent evolutionary phenomenon, providing insight into the processes of adaptation through natural selection, but it has serious economic and public health implications. In the arms race between arthropods and humans, the mosquito Anopheles gambiae, the main vector of malaria, seems to have just moved up a gear with the emergence of a resistant duplicated allele of the gene encoding acetylcholinesterase 1 (AChE1).
AChE1 is a critical enzyme in nerve transmission and the target of two of the most commonly used types of insecticides (organophosphates, OPs, and carbamates, CXs). Like several other mosquito species (including Culex pipiens, the well-studied vector of West Nile virus), An. gambiae displays resistance due to a single amino-acid substitution, from a glycine to a serine at the position 119, in the AChE1 catalytic site (G119S). In C. pipiens, there is direct and indirect evidence that the resistance allele (ace-1R) entails a large fitness cost, probably due to the mutated AChE1 having a much lower level of activity. Homozygous ace-1Rmosquitoes survive in the presence of insecticide, but are rapidly outcompeted in the absence of insecticide (see review in ). Heterozygotes are subject to smaller costs than resistant homozygotes in the absence of insecticide. In treated areas, they survive better than susceptible homozygotes, but are less resistant than ace-1Rhomozygotes. Due to the patchy nature of mosquito control, the generalist heterozygote is advantaged across treated and non-treated areas, although the more specialist resistant and susceptible homozygotes are locally selected in treated and non-treated environments respectively. Moreover, heterozygotes cannot invade due to the segregation burden leading to the loss of the advantage in half of their progeny.
Several duplicated alleles (ace-1D) have recently appeared, which link a susceptible and a resistant copy of the ace-1 gene on the same chromosome . Duplication thus creates a "permanent heterozygote" allele. The first case of ace-1 gene duplication was recently discovered in An. gambiae . Molecular analysis showed this duplicated allele (Ag-ace-1D) to be present at several sites and to have probably spread among the two molecular forms S and M of An. gambiae s.s, by introgression.
Unfortunately, it is not possible to design a simple test for studying the frequency of Ag-ace-1Ddue to the lack of features specific to this duplication, as with available genotyping methods carriers of this duplicated allele cannot be distinguished from classical heterozygotes. Thus an indirect method previously developed for C. pipiens was used to estimate Ag-ace-1Dfrequency in the field . The results of this analysis and the potential consequences for An. gambiae population management and on malaria control are discussed.
The study area is shown in Figure 1. Both published data [5, 7] and data from new samples were used. The date and location of the sampling sites are shown in Table 1. For each locality, several ponds where sampled in an area a few hundred meter-squares to insure a representative sample of the local population.
All samples were collected at the larval stage and reared to adulthood in the laboratory. Genomic DNA was extracted from each field mosquito. The protocol used is a simplified version of Collins et al. : a single mosquito is homogenized in a 1.5 ml Eppendorf tube containing 200 μl of CTAB buffer (100 mM Tris HCL, pH 8.0, 10 mM EDTA, 1.4 M NaCl, 2% CTAB) and incubated at 65°C for 5 min; then 200 μl of chloroform are added. After centrifugation (room temperature, 5 min, 12000 g), the supernatant is moved to a fresh tube, 200 μl of iso-propyl alcohol are added, and the mix is centrifuged again (12000 g, 15 min). After discarding supernatant, the pellet is washed with 70% ethanol, dried and resuspended in DNAse Free water. The molecular form of each individual was determined by a PCR-based test, as described in . The ace-1 genotype was assessed by RFLP analysis (see [5, 7]).
The presence of a duplicated allele causes an apparent excess of heterozygous [RS] phenotypes and thus a departure from the Hardy-Weinberg proportions expected with two alleles only (ace-1Rand ace-1S) . This departure is related to the frequency of the duplicated allele and was used to estimate Ag-ace-1Dfrequency in An. gambiae populations. The presence of Ag-ace-1Dwas investigated by fitting two models to the phenotypic data for each sample independently: i) a two-alleles-only model (ace-1Rand ace-1S) and ii) a three-allele model, adding the duplicated allele Ag-ace1D. The frequency of the duplicated allele was estimated from the excess of heterozygotes observed in each sample, assuming that this excess was due exclusively to the presence of Ag-ace-1D. This method is not as accurate as a direct identification of genotypes, but the two methods gave highly concordant results for field samples of C. pipiens . This indirect estimate of Ag-ace1Dfrequency may be biased if the genotypes are not in Hardy-Weinberg equilibrium. However, such a bias is not expected as An. gambiae populations large size prevents drift and as no overdominance leading to heterozygote excess was ever found for resistance . Moreover, previous studies of neutral markers in An. gambiae show either no departure from Hardy-Weinberg expectations or a deficit in heterozygotes, but never an excess, ensuring that this method is conservative (e.g. [11, 12]).
For each sample, the expected phenotypic distributions were calculated for the S and M molecular forms, using allelic distributions and assuming the ace-1 locus to be at Hardy-Weinberg equilibrium. Phenotype was considered to be a three-state random variable ([RR] corresponding to (R/R) genotype, [RS] corresponding to (R/S), (D/S), (D/R) and (D/D) genotypes, and [SS] corresponding to (S/S) genotype). The log-likelihood of a sample was calculated from the phenotypic multinomial distribution. Let n ij and f ij be the observed number and expected frequency of individuals with phenotype i in population j, respectively. The log likelihood L of observing all the data is proportional to
L was maximized, using the Metropolis algorithm [13, 14]. Model likelihoods were compared using F-tests: by construction, the three-allele model has a higher likelihood, but the presence of the duplication is considered to be confirmed only if the likelihood of this model is significantly higher than that of the two-allele model (significant P-value). The support limits for the frequency of each allele were also estimated. Finally, the P-values obtained were corrected for multiple testing, using Hommel's sequential Bonferroni correction . The different samples from each country were also pooled to get a higher statistical power at a larger geographical scale (Table 2), and the same analysis as for independent collection sites has been done. As pooling data from different populations is likely to result in a heterozygote deficit (Wahlund effect), this analysis is likely to underestimate any global excess of heterozygotes , making the detection of such excess more significant.
Results and discussion
The frequency of the recently discovered duplicated allele of the ace-1 gene in An. gambiae, Ag-ace-1D, was investigated in natural populations from West Africa by considering the departure from Hardy-Weinberg proportions caused by its presence .
Figure 1 shows the predicted spatial distribution of Ag-ace-1Din the S and M forms of An. gambiae, as shown by previous molecular investigations and analyses of heterozygote excess. The probability of Ag-ace-1Dbeing present was significant in nine samples (five after Bonferroni correction) from Ivory Coast (four samples) and Burkina Faso (five samples), in both M (two samples) and S (seven samples) molecular forms of An. gambiae (Figure 1 and Table 2). In these samples, the frequency of Ag-ace-1Dwas up to 0.65, with the lowest significant frequency being 0.24, consistent with the expected highly conservative output of the method used. Indeed, this method will detect low frequencies only in large samples; for example, Ag-ace-1Dwas not detected with this method in one of the analysed populations (Boromo, population #26, Table 2), whereas molecular methods proved this duplication to be present . The frequency and the geographic distribution of this duplication are therefore probably underestimated. For example, the analysis of each Benin population independently did not provide any indication supporting the presence of the duplication in this country (Figure 1 and Table 2). Nevertheless, the pooled analysis yields a P-value of 0.052, which points toward the potential presence of Ag-ace-1Das this method underestimate the excess of heterozygotes and thus its frequency. However, more data are required to confirm the presence of the duplicated allele in Benin (Table 2). The complete lack of variation of the Ag-ace-1Dsequence over several countries  indicates that this allele was generated by a single genetic event and its current distribution suggests that it is probably spreading.
Unfortunately, the spread of this new resistance allele poses a potential major threat to public health, as An. gambiae is the main vector of malaria. Indeed, several studies of a similar allele in C. pipiens have indicated that the duplication entails a lower fitness cost than the single-copy resistance gene, ace-1R[4, 6] (but see ). This is probably also the case for An. gambiae, as the mutated AChE1 gene is also associated with a strong decrease in enzyme activity . The presence and spread of the Ag-ace-1Dallele may greatly impede An. gambiae control strategies designed to maintain resistance alleles at low frequencies through the use of different insecticides with no cross-resistance in a mosaic or rotation strategy. It has been clearly demonstrated [18, 19] that the efficiency of such strategies increases with the fitness cost of resistance.
Insecticides for controlling malaria vectors are a major weapon in the battle between humans and malaria. Unfortunately, these insecticides exert strong selection pressure on vector populations, causing the spread of resistance genes, such as the resistance allele observed at the ace-1 locus in An. gambiae. The long-term use of an insecticide promotes the selection of new resistant variants, with a high risk of selecting a low (or null)-cost variant. The ace-1 duplicated allele recently appeared in An. gambiae is probably an example of such a low-cost variant. It is shown here that the presence of this duplicated allele, known from the molecular analysis of a few mosquitoes in some samples from Burkina Faso and Ivory Coast  is largely distributed in several countries of Western Africa, sometimes at high frequencies, and that it is probably spreading.
To prevent such spreads of resistance genes, it is crucial to develop the largest possible number of complementary means of control (e.g. larval insecticides, mosquito nets, biological agents, etc.) and to use them wisely. However, the emergence of ace-1 duplication in natural populations of An. gambiae, has just given mosquitoes the edge in this particular battle, seriously undermining our efforts to control vector populations and increasing the risk of malaria outbreaks.
Bills P: A new database of pesticide resistant insects and mites (Arthropods). Pesticide notes. 2001, 14: 2-4.
Weill M, Lutfalla G, Mogensen K, Chandre F, Berthomieu A, Berticat C, Pasteur N, Philips A, Fort P, Raymond M: Insecticide resistance in mosquito vectors. Nature. 2003, 423: 136-137. 10.1038/423136b.
Weill M, Labbé P, Duron O, Pasteur N, Fort P, Raymond M: Insecticide resistance in the mosquito Culex pipiens : towards an understanding of the evolution of ace genes. Insect evolutionary ecology. Edited by: Fellowes MDE, Holloway GJ, Rolff J. 2005, Oxon, UK: CABI publishing
Labbé P, Berthomieu A, Berticat C, Alout H, Raymond M, Lenormand T, Weill M: Independent duplications of the acetylcholinesterase gene conferring insecticide resistance in the mosquito Culex pipiens . Mol Biol Evol. 2007, 24: 1056-1067. 10.1093/molbev/msm025.
Djogbénou L, Chandre F, Berthomieu A, Dabire R, Koffi A, Alout H, Weill M: Evidence of introgression of the ace-1 Rmutation and of the ace-1 duplication in West African Anopheles gambiae s. s. PLoS ONE. 2008, 3: e2172-10.1371/journal.pone.0002172.
Lenormand T, Guillemaud T, Bourguet D, Raymond M: Appearance and sweep of a gene duplication: adaptive response and potential for new functions in the mosquito Culex pipiens . Evolution. 1998, 52: 1705-1712. 10.2307/2411343.
Djogbénou L, Dabire R, Diabate A, Kengne P, Akogbeto M, Hougard JM, Chandre F: Identification and geographic distribution of the ace-1 Rmutation in the malaria vector Anopheles gambiae in South-Western Burkina Faso, West Africa. Am J Trop Med Hyg. 2008, 78: 298-302.
Collins FH, Mendez MA, Rasmussen MO, Mehaffey PC, Besansky NJ, Finnerty V: A ribosomal RNA gene probe differentiates member species of the Anopheles gambiae complex. Am J Trop Med Hyg. 1987, 37: 37-41.
Favia G, Torre Ad, Bagayoko M, Lanfrancotti A, Sagnon NF, Touré YT, Coluzzi M: Molecular identification of sympatric chromosomal forms of Anopheles gambiae and further evidence of their reproductive isolation. Insect Mol Biol. 1997, 6: 377-383. 10.1046/j.1365-2583.1997.00189.x.
Djogbénou L, Weill M, Hougard JM, Raymond M, Akogbéto M, Chandre F: Characterization of Insensitive Acetylcholinesterase (ace-1 R) in Anopheles gambiae (Diptera: Culicidae): Resistance Levels and Dominance. J Med Entomol. 2007, 44: 805-810. 10.1603/0022-2585(2007)44[805:COIAAI]2.0.CO;2.
Lanzaro GC, Touré YT, Carnahan J, Zheng L, Dolo G, Traoré S, Petrarca V, Vernick KD, Taylor CE: Complexities in the genetic structure of Anopheles gambiae populations in west Africa as revealed by microsatellite DNA analysis. PNAS. 1998, 95: 14260-14265. 10.1073/pnas.95.24.14260.
Lehmann T, Hawley WA, Kamau L, Fontenille D, Simard F, Collins FH: Genetic differentiation of Anopheles gambiae populations from East and West Africa: comparison of microsatellite and allozyme loci. Heredity. 1996, 77: 192-200. 10.1038/hdy.1996.124.
Labbé P, Lenormand T, Raymond M: On the worldwide spread of an insecticide resistance gene: a role for local selection. J Evol Biol. 2005, 18: 1471-1484.
Lenormand T, Bourguet D, Guillemaud T, Raymond M: Tracking the evolution of insecticide resistance in the mosquito Culex pipiens . Nature. 1999, 400: 861-864. 10.1038/23685.
Hommel G: A stagewise rejective multiple test procedure based on a modified Bonferroni test. Biometrika. 1988, 75: 383-386. 10.1093/biomet/75.2.383.
Labbé P, Berticat C, Berthomieu A, Unal S, Bernard C, Weill M, Lenormand T: Forty years of erratic insecticide resistance evolution in the mosquito Culex pipiens . PLoS Genetics. 2007, 3: e205-10.1371/journal.pgen.0030205.
Alout H, Djogbenou L, Berticat C, Chandre F, Weill M: Comparison of Anopheles gambiae and Culex pipiens acetycholinesterase 1 biochemical properties. Comp Biochem Physiol B-Biochem Mol Biol. 2008, 150: 271-277. 10.1016/j.cbpb.2008.03.008.
Carrière Y, Deland J-P, Roff DA, Vincent C: Life-history cost associated with the evolution of insecticide resistance. Proc R Soc Lond B. 1994, 258: 35-40. 10.1098/rspb.1994.0138.
Lenormand T, Raymond M: Resistance management: the stable zone strategy. Proc R Soc Lond B. 1998, 265: 1985-1990. 10.1098/rspb.1998.0529.
This work was funded in part by the ANR Morevol Sante-Environnement (Ministère délégué à la Recherche) and in part by the Institut de Recherche Pour le Développement (IRD). It is the contribution 2009-025 of the Institut des Sciences de L'Evolution de Montpellier (ISEM, UMR CNRS-UM2 5554).
The authors declare that they have no competing interests.
LD designed the study, acquired the data and wrote the manuscript. PL analysed the data and wrote the manuscript. NP, FC and MW contributed to the design of the study and for draft and revision of the manuscript. All authors read and approved the final manuscript.
Luc Djogbénou, Pierrick Labbé contributed equally to this work.
Authors’ original submitted files for images
Below are the links to the authors’ original submitted files for images.