# Evidence of Introgression of the ace-1R Mutation and of the ace-1 Duplication in West African Anopheles gambiae s. s.

Luc Djogbénou, Fabrice Chandre, Arnaud Berthomieu, Roch Dabiré, Alphonsine Koffi, Haoues Alout, Mylène Weill

► To cite this version:

Luc Djogbénou, Fabrice Chandre, Arnaud Berthomieu, Roch Dabiré, Alphonsine Koffi, et al.. Evidence of Introgression of the ace-1R Mutation and of the ace-1 Duplication in West African Anopheles gambiae s. s.. PLoS ONE, 2008, 3 (5), pp.e2172. 10.1371/journal.pone.0002172. halsde-00527414

HAL Id: halsde-00527414

https://hal.science/halsde-00527414

Submitted on 29 May 2020

HAL is a multi-disciplinary open access
archive for the deposit and dissemination of sci-
entific research documents, whether they are pub-
lished or not. The documents may come from
teaching and research institutions in France or
abroad, or from public or private research centers.

L'archive ouverte pluridisciplinaire HAL, est
destinée au dépôt et à la diffusion de documents
scientifiques de niveau recherche, publiés ou non,
émanant des établissements d'enseignement et de
recherche français ou étrangers, des laboratoires
publics ou privés.

<figure><img src="image_1.png" /></figure>

Distributed under a Creative Commons Attribution 4.0 International License
# Evidence of Introgression of the *ace-1*<sup>R</sup> Mutation and of the *ace-1* Duplication in West African *Anopheles gambiae* s. s.

Luc Djogbénou<sup>1,2</sup>, Fabrice Chandre<sup>1</sup>, Arnaud Berthomieu<sup>3,4</sup>, Roch Dabiré<sup>5</sup>, Alphonsine Koffi<sup>6</sup>, Haoues Alout<sup>3,4</sup>, Mylène Weill<sup>3,4*</sup>

1 Institut de Recherche pour le Développement, Cotonou, Benin, 2 Centre de Recherche Entomologique de Cotonou (CREC), Cotonou, Benin, 3 Université Montpellier 2, Montpellier, France, 4 CNRS, Institut des Sciences de l'Evolution, Equipe Génétique de l'Adaptation, Montpellier, France, 5 Institut Régional des Sciences de Santé, Laboratoire d'Entomologie et de Parasitologie, Bobo-Dioulasso, Burkina Faso, 6 Institut Pierre Richet (IPR)/Institut National de Santé Publique, Abidjan, Côte d'Ivoire

## Abstract

**Background:** The role of inter-specific hybridisation is of particular importance in mosquito disease vectors for predicting the evolution of insecticide resistance. Two molecular forms of *Anopheles gambiae* s.s., currently recognized as S and M taxa, are considered to be incipient sibling species. Hybrid scarcity in the field was suggested that differentiation of M and S taxa is maintained by limited or absent gene flow. However, recent studies have revealed shared polymorphisms within the M and S forms, and a better understanding of the occurrence of gene flow is needed. One such shared polymorphism is the G119S mutation in the *ace-1* gene (which is responsible for insecticide resistance); this mutation has been described in both the M and S forms of *A. gambiae* s.s.

**Methods and Results:** To establish whether the G119S mutation has arisen independently in each form or by genetic introgression, we analysed coding and non-coding sequences of *ace-1* alleles in M and S mosquitoes from representative field populations. Our data revealed many polymorphic sites shared by S and M forms, but no diversity was associated with the G119S mutation. These results indicate that the G119S mutation was a unique event and that genetic introgression explains the observed distribution of the G119S mutation within the two forms. However, it was impossible to determine from our data whether the mutation occurred first in the S form or in the M form. Unexpectedly, sequence analysis of some resistant individuals revealed a duplication of the *ace-1* gene that was observed in both *A. gambiae* s.s. M and S forms. Again, the distribution of this duplication in the two forms most likely occurred through introgression.

**Conclusions:** These results highlight the need for more research to understand the forces driving the evolution of insecticide resistance in malaria vectors and to regularly monitor resistance in mosquito populations of Africa.

**Citation:** Djogbénou L, Chandre F, Berthomieu A, Dabiré R, Koffi A, et al. (2008) Evidence of Introgression of the *ace-1*<sup>R</sup> Mutation and of the *ace-1* Duplication in West African *Anopheles gambiae* s. s. PLoS ONE 3(5): e2172. doi:10.1371/journal.pone.0002172

**Editor:** Dee A. Carter, University of Sydney, Australia

**Received** December 21, 2007; **Accepted** April 1, 2008; **Published** May 14, 2008

**Copyright:** © 2008 Djogbénou et al. This is an open-access article distributed under the terms of the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original author and source are credited.

**Funding:** This work received financial support in part by the PAL+Programme, no. 2002-45 (Ministère délégué à la Recherche) and Institut de Recherche Pour le Développement. It was also financed by the ANR Morevol Sante-Environnement (Ministère délégué à la Recherche) and Contribution 2008-024 of the Institut des Sciences de L'Evolution de Montpellier (UMR CNRS-UM2 5554).

**Competing Interests:** The authors have declared that no competing interests exist.

* E-mail: mylene.weill@univ-montp2.fr

## Introduction

Understanding the dynamics of genes associated with insecticide resistance between closely related taxa is a key component in establishing effective long-term resistance management strategies. Adaptive genes such as those producing insecticide resistance allow for the ability to detect rare hybridization events within a group of related taxa. There are three major target sites for most insecticides: the g-aminobutyric acid (GABA) receptor is the target of cyclodiene insecticides, the voltage-dependent sodium channel is the target site for DDT and pyrethroids, and acetylcholinesterase (AChE1, EC 3.1.1.7) is quasi-irreversibly inhibited by organophosphorous and carbamate compounds, which are substrate analogues. In Africa, propoxur resistance was first detected in a population of *Anopheles gambiae* Giles from Côte d'Ivoire [1]. Insensitive AChE1 was next confirmed as the

resistance mechanism [2]. As in *Culex pipiens* L., *Anopheles albimanus* Weidemann and *Culex vishnui* Theobald [3–6], the resistance in *A. gambiae* results from a single mutation in the *ace-1* gene that leads to a Gly-to-Ser substitution at position 119 (according to the *Torpedo californica* nomenclature) [7]. The resistant allele *ace-1*<sup>R</sup> confers resistance to organophosphates and carbamate compounds. This resistance potentially represents a threat to the implementation of malaria prevention programmes based on the use of insecticides.

The *Anopheles gambiae* complex includes some of the most important malaria vector species of Sub-Saharan Africa. Genetic differentiation occurs within this highly polymorphic complex that is subdivided into five cytoforms. These cytoforms differ in their arrangements of chromosomal inversion and appear more or less genetically isolated in the field [8,9]. In addition, studies using molecular markers such as X-linked ribosomal DNA revealed the presence of two distinct molecular forms within *A. gambiae* s.s.

<figure><img src="image_2.png" /></figure>
These forms have been referred to as M and S forms [10] and are now considered as units of an on-going incipient speciation process [11]. A high deficit of M/S hybrids observed in the field has been attributed to assortative mating because hybrids are readily obtained in the laboratory [11–17]. Thus, it has been suggested that the differentiation into M and S taxa is maintained by limited or absent gene flow [11]. However, the absence of gene flow is still debated because recent work using molecular markers showed that differentiation between the M and S forms is only present in few regions of the genome: a small region close to the centromere of the 2nd chromosome, and the centromeric region of the X chromosome that contains the rDNA that defines these molecular forms [14,18–20]. Evolutionary histories of closely related species are often complicated by the existence of shared polymorphisms (inherited from their common ancestor) as well as possible ongoing gene flow between them [21]. Today, the debate over whether the limited differentiation among molecular forms of *A. gambiae* s.s. reflects insufficient time for more differentiation to have occurred versus introgression is not resolved. However, some recent introgressive hybridization phenomena within *A. gambiae* s.s. are well documented [14–16].

Recently, studies have been conducted to estimate the distribution of the *ace-1*<sup>R</sup> G119S mutation in populations of *A. gambiae* s.s. from Burkina Faso [22]. The mutation was found in both M and S forms where they were sympatric. It is important to determine how often the G119S resistance mutation occurs in natural populations in order to understand its evolutionary history within the *A. gambiae* s.s. complex. In the *Culex pipiens* mosquito complex, the G119S mutation appeared several times independently. Two distinct mutation events were first reported in *C. pipiens* and *C. quinquefasciatus* species [5]. A third event was next described in *C. quinquefasciatus* from China [23] and a polymorphic intron sequence located upstream of exon 3 revealed two more events within *C. quinquefasciatus* [24].

In this paper, we investigated whether the *ace-1*<sup>R</sup> mutation has arisen in M and S molecular forms through two independent mutations or through genetic introgression. We analysed the polymorphism of the DNA sequence associated with the G119S mutation (both introns upstream and downstream from exon 3) in several resistant and susceptible individuals of *A. gambiae* s.s of both M and S molecular forms from Burkina Faso, Benin and Côte d'Ivoire. Our data clearly indicate that the G119S mutation was a unique event and that genetic introgression explains the observed distribution of the G119S mutation within the two forms.

## Materials and Methods

### Mosquito samples

Larvae were collected at several locations from Burkina Faso, Benin and Côte d'Ivoire (Figure 1) and reared in the laboratory until emergence. Adult *A. gambiae* s.l. were sorted from other anophelines according to morphological identification keys [25] and were kept at –20°C for molecular analysis.

### DNA extraction, PCR identification of *A. gambiae* M and S forms, and ace-1 genotyping

Total genomic DNA was extracted from individual mosquitoes following a slightly modified procedure described by Collins [26]. Species in the *A. gambiae* complex and molecular forms of *A. gambiae* s.s. were identified using the methods described by Scott [27] and Favia [10], respectively.

For ace-1 genotyping, 1–10 ng of genomic DNA was amplified with the primers Ex3AGdir and Ex3AGrev [6]. PCR was conducted in 25 µL volumes containing 2.5 µL 10× Taq DNA polymerase buffer, 200 µM deoxynucleoside triphosphate (dNTP), 0.1 U Taq DNA polymerase (Qiagen, France) and 10 pmol of each primer. PCR conditions included an initial denaturation step

<figure><img src="image_3.png" /><figcaption>**Figure 1. Map of Benin, Burkina Faso and Côte d'Ivoire showing the study sites.** (•): study locality. Allelic frequencies of the *ace-1*<sup>R</sup> mutation in *A. gambiae* s.s. are indicated in each molecular form M and S. doi:10.1371/journal.pone.0002172.g001</figcaption></figure>

<figure><img src="image_4.png" /></figure>
at 94°C for 5 min followed by thirty five cycles of 94°C for 30 s, 54°C for 30 s, and 72°C for 30 s, and a final extension at 72°C for 5 min. Fifteen µL of the PCR product were digested with 5 U of AhaI restriction enzyme (Promega, France) in a final volume of 25 µL. The restriction products were fractionated on 2% agarose gels (Agarose, MP, Sigma) with TBE buffer (0.089 M Tris, 0.089 M boric acid and 0.5 M EDTA, pH 8.0), stained with ethidium bromide, and observed under UV light.

## Partial sequences of ace-1 alleles

DNA from single mosquitoes was amplified using two specific primers: AgEx2dir1 5'-AGG TCA CGG TGA GTC CGT ACG A-3' and AgEx4rev2 5'-AGG GCG GAC AGC AGA TGC AGC GA-3' generating an 817 bp fragment that included part of exon 2, intron 2, whole exon 3, intron 3 and part of the exon 4. PCR was carried out with ~20 ng of genomic DNA, 10 pmol of each primer, 10 µm of each dNTP, 2.5 units of Taq polymerase (High Fidelity) in 1x reaction buffer (Tris-HCl [pH 9.0; 75 mM], (NH₄)₂SO₄ [20 mM], Tween 20 [0.1 g.liter⁻¹], and MgCl₂ [1.25 mM]) for a final volume of 50 µL. The PCR mixture was subjected to 30 amplification cycles (93°C for 30 s, 56° for 30 s and 72°C for 50 s). PCR products were purified using the QIAquick Gel Extraction Kit (QIAGEN). For homozygous resistant mosquitoes, the purified PCR product was directly sequenced. For heterozygous and susceptible homozygous mosquitoes, the purified PCR product was cloned using the TOPO_ Cloning Kit (Invitrogen, Paisley, UK) according to the manufacturer instructions. Clones were then screened for the presence of the G119S substitution [6]. At least six clones were sequenced for each mosquito, i.e., three susceptible allele sequences and three resistant allele sequences for heterozygotes and six susceptible allele sequences for susceptible homozygotes. Sequencing was conducted on an ABI Prism 310 sequencer (BigDye Terminator Kit, Applied Biosystems, Foster City, CA). Both strands of the genome fragment were sequenced for each clone using the primers AgEx2dir1 and AgEx4rev2.

## Sequence analysis

Sequences were aligned with Multalin software [http://www.prodes.toulouse.inra.fr/multalin/multalin.html] [28] using the Kisumu strain as a reference. Similarity between various sequences was assessed with ClustalW [29]. Sequence polymorphism and nucleotide diversity were analysed with DnaSP, v.4.10.3 [30] using Nei's p index [31,32]. Deduced amino acid sequences were obtained with ClustalW [29] to determine whether the mutations identified were synonymous. When a mutation was not synonymous, the position of the new amino acid was sought on a 3-dimensional model of the *Torpedo californica* AChE (pdb 1EA5), using Swiss-PdbViewer v. 3.7 [33] to determine whether this mutation might be located at a site crucial for protein activity.

## Results

### Distribution of the ace-1<sup>R</sup> mutation in the populations sampled

Two hundred and eighty two mosquitoes from Burkina Faso (n = 83), Benin (n = 135) and Côte d'Ivoire (n = 64) were assayed for molecular form and analyzed for their *ace-1* genotype. The M and S molecular forms are sympatric in seven of the ten localities selected for this study (Table 1). The G119S mutation was distributed at various frequencies in the S and M forms, depending on the location. This mutation was not detected in the M form of *A. gambiae* s.s from Benin (0% compared to 11% in the S form) and was weakly detected in Burkina Faso (2% in M form compared to

<table><caption>**Table 1.** Genotype distribution of M and S molecular forms in different locations.</caption><thead><tr><th rowspan="2">Locality</th><th rowspan="2">ace.1 genotype</th><th colspan="2">sample</th></tr><tr><th>S Sample</th><th>M Sample</th></tr></thead><tbody><tr><td colspan="4"><strong>Benin</strong></td></tr><tr><td rowspan="3">• Abomey</td><td>RR</td><td></td><td></td></tr><tr><td>RS</td><td>16 (2)</td><td></td></tr><tr><td>SS</td><td>52 (1)</td><td>3</td></tr><tr><td rowspan="3">• Bohicon</td><td>RR</td><td></td><td></td></tr><tr><td>RS</td><td>1 (1)</td><td></td></tr><tr><td>SS</td><td>2 (1)</td><td>2 (2)</td></tr><tr><td rowspan="3">• Paouignan</td><td>RR</td><td></td><td></td></tr><tr><td>RS</td><td>7 (3)</td><td></td></tr><tr><td>SS</td><td>32</td><td></td></tr><tr><td rowspan="3">• Zogbodomey</td><td>RR</td><td></td><td></td></tr><tr><td>RS</td><td>1 (1)</td><td></td></tr><tr><td>SS</td><td>6</td><td>13 (3)</td></tr><tr><td colspan="4"><strong>Burkina Faso</strong></td></tr><tr><td rowspan="3">• Boromo</td><td>RR</td><td></td><td></td></tr><tr><td>RS</td><td>1 (1)</td><td>2 (2)</td></tr><tr><td>SS</td><td>1 (1)</td><td>35 (2)</td></tr><tr><td rowspan="3">• Dioulassoba</td><td>RR</td><td>1</td><td></td></tr><tr><td>RS</td><td>15 (2)</td><td></td></tr><tr><td>SS</td><td>7</td><td>1</td></tr><tr><td rowspan="3">• Seguere</td><td>RR</td><td></td><td></td></tr><tr><td>RS</td><td>6 (2)</td><td></td></tr><tr><td>SS</td><td>2 (1)</td><td>12 (2)</td></tr><tr><td colspan="4"><strong>Côte d'Ivoire</strong></td></tr><tr><td rowspan="3">• Yaokofikro</td><td>RR</td><td>2</td><td></td></tr><tr><td>RS</td><td>15 (1)</td><td></td></tr><tr><td>SS</td><td>3</td><td></td></tr><tr><td rowspan="3">• Toumbokro</td><td>RR</td><td></td><td>1</td></tr><tr><td>RS</td><td>5 (1)</td><td>11</td></tr><tr><td>SS</td><td></td><td>7 (1)</td></tr><tr><td rowspan="3">• Niamoue</td><td>RR</td><td></td><td>1 (1)</td></tr><tr><td>RS</td><td></td><td>19 (2)</td></tr><tr><td>SS</td><td></td><td></td></tr></tbody></table>

The number of individuals sequenced is shown into brackets.

doi:10.1371/journal.pone.0002172.t001

31% in the S form). However, it was strongly detected in Côte d'Ivoire (44% of the M form and 50% of the S form). The number of homozygotes for the *ace-1*<sup>R</sup> mutation was very low, even in localities from Côte d'Ivoire, despite the high frequency of resistant heterozygotes in some samples (Figure 1).

### Variability of Susceptible and Resistant Alleles

To search for polymorphism in the *ace-1* gene, 33 individuals of *A. gambiae* s.s. from several localities of Benin, Burkina Faso and Côte d'Ivoire, and two *Anopheles arabiensis* Patton from Bohicon (Benin), were analysed. Sequencing of susceptible and resistant *ace-1* alleles was conducted on an 817 bp fragment from M form 15 individuals and 18 S form individuals (Table S1). Forty-eight distinct alleles were identified from the 35 individuals studied (33

<figure><img src="image_5.png" /></figure>
*A. gambiae* s.s. and two *A. Arabiensis*; Table S1, with identical alleles removed). The sequence analysis showed that all sequences carrying the resistance mutation (i.e., 14 S and 5 M individuals) were identical regardless of their geographical origin and molecular form; these were annotated as M&S-R (Table S1). We measured the variability of the *ace-1* alleles in our sample. Seventy-four positions exhibited polymorphism: 3, 10, 34, 20 and 7 polymorphic positions in exon 2, intron 2, exon 3, intron 3, and exon 4, respectively (Figure 2). A 3 bp insertion was detected in intron 2 in the sequence of some M and S individuals (Table S1). In addition to the mutation at position 363 corresponding to the G119S substitution, 42 variable positions were shared by at least two individuals. The nucleotide diversity of the coding exon 3 was estimated as Pi: 0.00634; no nucleotide position was found specific to the M or S form. A distance tree constructed (Figure 3) showed low bootstrap values although a few nodes were supported. The distance tree did not reveal any organisation associated with geographical origin or with a particular molecular form. Moreover, the two sequences from *A. arabiensis* segregated among *A. gambiae* s.s. alleles without forming an identifiable outgroup. Among the susceptible alleles, none was identical to that of the resistant allele (excluding the resistant mutation). However, two alleles belonging to the M form, M-Zogbodomeylb and M-Seguere37b, were close to the resistant one (Figure 3). M-Zogbodomeylb differs from the resistant allele by 8 mutations, and M-Seguere37b differs by 11 mutations plus a 3 bp insertion in intron 2. In addition, a susceptible allele was shared by several mosquitoes belonging to both forms (Figure 3).

## Protein Sequence Variability

The comparison of the coding region of each susceptible sequence to the single resistant one (651 bp) revealed 45 polymorphic sites (excluding the G119S mutation) but no insertions/deletions. Nucleotide diversity was low (Pi=0.00810). Only three non-synonymous mutations were identified in three sequences among the 48 distinct alleles analysed (Table S1). Using the *Torpedo* nomenclature [34], we observed the following substitutions: Gly to Ser at position 24 in a Toumbokro2b individual, Asn to Ser at position 140 in a Bohicon72a individual and Asp to Val at position 189 in a Niamoue19 individual. The AChE1 structural model based on that of *Torpedo californica* (PDB: 1EA5) shows that the amino acid substitutions are located at the periphery of the enzyme, far from the active site and its entrance. These amino acids are not likely to have any function in the catalytic activity of AChE1 (Figure 4).

## Heterozygotes harboring three alleles suggest *ace-1* duplication

For each resistant heterozygote, at least 3 susceptible clones were sequenced in effort to avoid any polymerase mistakes. Four mosquitoes from Burkina Faso and one from Côte d'Ivoire

contained three different alleles (one resistant and two distinctly different susceptible alleles). The PCR, cloning and sequencing were repeated to avoid any possibility of contamination, and more susceptible clones were sequenced to confirm the first results. For the heterozygote M-Boromo41 and M-Niamoue19, the two susceptible copies differed by 13 positions and 6 positions, respectively. For S-Seguere48, S-Dioulassoba35 and S-Seguere12, the two susceptible copies differed by 14 positions, 13 positions and 12 positions, respectively. In addition, the comparison of susceptible alleles among these mosquitoes showed that they shared one identical susceptible sequence (Figure 3). This common sequence is also shared by other mosquitoes (Figure 3).

## Discussion

In this study, we attempted to describe how the *ace-1*<sup>R</sup> mutation occurred in the two molecular forms M and S of *A. gambiae* s.s., i.e., whether by independent evolutionary gains or by introgression. The sequence of the resistant allele was identical in all homozygous and heterozygous individuals regardless of their geographic origin or molecular form. Although the recorded nucleotide diversity of *ace-1* gene was low in coding sequences (Pi = 0.00810), it increased significantly when non-coding introns were considered (Pi = 0.01224). The probability of the G119S mutation to have been gained twice independently at the same position on a same susceptible allele in both molecular forms is extremely low. Thus, the presence of the *ace-1*<sup>R</sup> mutation in both M and S forms probably resulted from introgression. However, an alternative explanation is that sibling species sharing the same allele could also reflect ancestral retention. This has been documented in other groups for mtDNA genes [35–37], and recently in the *A. gambiae* complex [38]. If this was the case, shared haplotypes will pre-date differentiation of S and M forms and will have been maintained by chance in the two forms through incomplete lineage sorting.

The resistant allele *ace-1*<sup>R</sup> has been thoroughly studied in *Culex pipiens* mosquitoes and was found to be associated with a high fitness cost [39–43]. The cost probably results from the excess of ACh in synapses because the enzymatic activity of the AChE1R is more than 60% lower than that of the AChE1S [44]. The activity of AChE1R was also found to decrease in *A. gambiae* resistant mosquitoes and a high fitness cost is probably also associated with AChE1R resistance in this species [45]. Thus, because of high fitness costs, G119S resistant mosquitoes cannot be detected after a few generations in the absence of selection pressure by organophosphorous or carbamate insecticides [24]. Since the use of pesticides is very recent (i.e., last few decades) compared to the differentiation of M and S forms, it is likely that the G119S mutation occurred in one form and introgressed in the other. From the analysis of sequence polymorphisms, it is difficult to determine in which form the mutation first occurred because there

<figure><img src="image_6.png" /><figcaption>**Figure 2.** Genomic structure and polymorphism of partial sequences of the *A. gambiae ace-1* gene. bp: number of bases for each sequence. Pi: Nucleotide diversity. ps: Number of polymorphic sites.
doi:10.1371/journal.pone.0002172.g002</figcaption></figure>

<figure><img src="image_7.png" /></figure>
<figure><img src="image_8.png" /><figcaption>**Figure 3.** Genetic distance tree of *ace-1* alleles in the M and S forms of *Anopheles gambiae* s.s. Samples were named as followed: S or M before the name of the locality and individual number. Where two distinct susceptible alleles were present in one mosquito, these are noted by "a" and "b" Blue characters are M form samples and black are S form. Variability of the single resistant allele M&S-R and all susceptible alleles of each form are presented and only supported bootstraps are shown. Two sequences of *Anopheles arabiensis* are included and highlighted in grey. For each allele, the sequence of *ace-1* considered encompassed part of exon 2, intron 2, exon 3, intron 3, and part of exon 4. The G119S mutation, selected for resistance to carbamate and organophosphate, has been removed to take only neutral variation into account.
doi:10.1371/journal.pone.0002172.g003</figcaption></figure>

were no sequences specific to one form. However, because the resistant sequence does form a significant cluster with two susceptible sequences that are both from the M form, we can speculate that the mutation probably occurred first in the M form.

In our study, we also showed that some mosquitoes carried three different *ace-1* alleles: one resistant and two distinct susceptible alleles (see Figure 3). This observation demonstrates the existence of a duplication of the *ace-1* gene that resulted in one

<figure><img src="image_9.png" /></figure>
<figure><img src="image_10.png" /><figcaption>**Figure 4.** Structural models of *A. gambiae* AChE1 (pdb : 2AZG) used to locate amino acids substitutions found in susceptible *A. gambiae* populations from Bénin. The three substitutions (G24, N140, D189) are represented in dark grey (Van der Waals grey spheres) as well as the 119 position involved in insecticide resistance. The catalytic triad (S200, E327 and H440) appears as black sticks. The arrow shows the entrance of the catalytic gorge. The backbone of the enzyme structure is rendered as a ribbon with secondary structure. doi:10.1371/journal.pone.0002172.g004</figcaption></figure>

susceptible copy and one resistant copy, both of which are on the same chromosome. Such a duplicated *ace-1*<sup>D</sup> allele has been described several times in *Culex pipiens* and results in a “fixed” heterozygous phenotype that displays the same resistance level but with a reduced fitness cost compared to the homozygous resistant phenotype [24,46–47]. The deficit of RR homozygotes in the analysed samples also suggests the presence of the duplication. The sequence comparison of susceptible alleles showed that the mosquitoes with three different alleles of *ace-1* all shared one susceptible sequence (Figure 3), indicating that the mosquitoes contain the same duplicated allele. Furthermore, this duplication

was observed in individuals of both molecular forms (Boromo41 and Niamoue19 individuals belong to the M form while Seguere48 and Dioulassoba35 belong to the S molecular form; Table S1). Thus, this suggests that the *ace-1* duplication also occurred in both forms through introgression. Monitoring the frequency of the *ace-1*<sup>R</sup> allele in both the M and S forms of *A. gambiae* s.s. will greatly improve our knowledge of such gene duplications evolution.

Currently, there are several, ongoing projects to develop new strategies of vector control with the aim of improving resistance management. It is clear that resistance towards pyrethroids is widespread in Africa and interest in using IRS (Indoor Residual Spraying) to control malaria vectors is resurging. This IRS strategy is preferentially based on the use of organophosphates and carbamates, either alone or in combination with pyrethroids, and the spread of the *ace-1*<sup>R</sup> mutation could represent a major threat to the effectiveness of this strategy. Furthermore, the presence of the *ace-1* duplication could impede strategies based on alternating compounds in space or time. Indeed, such strategies rely on the fitness cost associated with resistance that could decrease in presence of the duplication [24]. Understanding or predicting the spread of insecticide resistance genes into mosquito populations of the *A. gambiae* complex will be crucial for the development of effective methods to control the main malaria vector in Africa.

**Supporting Information**

**Table S1**

Found at: doi:10.1371/journal.pone.0002172.s001 (0.63 MB DOC)

**Acknowledgments**

We would like to thank Martin Akogbéto, Jean-Marc Hougard, Nicole Pasteur, Pierrick Labbé and Michel Raymond for their support during the study, Sébastien Koudénoukpo of the Centre de Recherche Entomologique de Cotonou for his help with the genotyping.

**Author Contributions**

Conceived and designed the experiments: MW FC LD. Performed the experiments: AB LD. Analyzed the data: MW HA LD. Contributed reagents/materials/analysis tools: MW AB RD AK LD. Wrote the paper: MW FC LD.

**References**

1. Elissa N, Mouchet J, Riviere F, Meunier JY, Yao K (1994) Susceptibility of *Anopheles gambiae* to insecticides in the Ivory Coast. Santé 4: 95-99.

2. N'Guessan R, Darriet F, Guillet P, Carnevale P, Traore-Lamizana M, et al. (2003) Resistance to carbosulfan in field populations of *Anopheles gambiae* from Côte d'Ivoire based on reduced sensitivity acetylcholinesterase. Med Vet Entomol 17: 19-25.

3. Alout H, Berthomieu A, Cui F, Tan Y, Berticat C, et al. (2007) Different amino-acid substitutions confer insecticide resistance through acetylcholinesterase 1 insensitivity in *Culex vishnui* and *Culex tritaeniorhynchus* (Diptera : Culicidae) from China. J Med Entomol 44: 463-469.

4. Weill M, Lutfalla G, Mogensen K, Chandre F, Berthomieu A, et al. (2003) Comparative genomics: Insecticide resistance in mosquito vectors. Nature 423: 136-137.

5. Weill M, Fort P, Berthomieu A, Dubois MP, Pasteur N, et al. (2002) A novel acetylcholinesterase gene in mosquitoes codes for the insecticide target and is non-homologous to the ace gene in *Drosophila*. Proc R Soc Lond B Biol Sci 269: 2007-2016.

6. Weill M, Malcolm C, Chandre F, Mogensen K, Berthomieu A, et al. (2004) The unique mutation in ace-1 giving high insecticide resistance is easily detectable in mosquito vectors. Insect Mol Biol 13: 1-7.

7. Massoulié J, Sussman JL, Doctor BP, Soreq H, Velan B, et al. (1992) Recommendations for nomenclature in cholinesterases. In: Shafferman A, Velan B, eds. Multidisciplinary approaches to cholinesterase functions. New York: Plenum Press. pp 285-288.

8. Coluzzi M, Petrarca V, Di Deco MA (1985) Chromosomal inversion intergradation and incipient speciation in *Anopheles gambiae*. Boll Zool 52: 45-63.

9. Coluzzi M, Sabatini A, della Torre A, Di Deco MA, Petrarca V (2002) A polytene chromosome analysis of the *Anopheles gambiae* species complex. Science 298: 1415-1418.

10. Favia G, Lanfrancotti A, Spanos L, Siden-Kiamos I, Louis C (2001) Molecular characterization of ribosomal DNA polymorphisms discriminating among chromosomal forms of *Anopheles gambiae* s.s. Insect Mol Biol 10: 19-23.

11. della Torre A, Fanello C, Akogbeto M, Dossou-yovo J, Favia G, et al. (2001) Molecular evidence of incipient speciation within *Anopheles gambiae* s.s. in West Africa. Insect Mol Biol 10: 9-18.

12. della Torre A, Tu ZJ, Petrarca V (2005) On the distribution and genetic differentiation of *Anopheles gambiae* s.s. molecular forms. Insect Biochem Mol Biol 35: 755-769.

13. Edillo FE, Toure YT, Lanzaro GC, Dolo G, Taylor CE (2002) Spatial and habitat distribution of *Anopheles gambiae* and *Anopheles arabiensis* (Diptera : Culicidae) in Banambani Village, Mali. J Med Entomol 39: 70-77.

14. Gentile G, Slotman M, Ketmaier V, Powell JR, Caccone A (2001) Attempts to molecularly distinguish cryptic taxa in *Anopheles gambiae* s.s. Insect Mol Biol 10: 25-32.

15. Taylor C, Toure YT, Carnahan J, Norris DE, Dolo G, et al. (2001) Gene flow among populations of the malaria vector, *Anopheles gambiae*, in Mali, West Africa. Genetics 157: 743-750.

16. Tripet F, Toure YT, Dolo G, Lanzaro GC (2003) Frequency of multiple inseminations in field-collected *Anopheles gambiae* females revealed by DNA analysis of transferred sperm. Am J Trop Med Hyg 68: 1-5.

<figure><img src="image_11.png" /></figure>
17. Tripet F, Toure YT, Taylor CE, Norris DE, Dolo G, et al. (2001) DNA analysis of transferred sperm reveals significant levels of gene flow between molecular forms of *Anopheles gambiae*. Mol Ecol 10: 1725-1732.

18. Lehmann T, Licht M, Elissa N, Maega BTA, Chimumbwa JM, et al. (2003) Population structure of *Anopheles gambiae* in Africa. J Hered 94: 133-147.

19. Stump AD, Fitzpatrick MC, Lobo NF, Traore S, Sagnon NF, et al. (2005) Centromere-proximal differentiation and speciation in *Anopheles gambiae*. PNAS USA 102: 15930-15935.

20. Turner TL, Hahn MW, Nuzhdin SV (2005) Genomic islands of speciation in *Anopheles gambiae*. Plos Biology 3: 1572-1578.

21. Arbogast BS, Edwards SV, Wakeley J, Beerli P, Slowinski JB (2002) Estimating divergence times from molecular data on phylogenetic and population genetic timescales. Annu Rev Ecol Syst 33: 707-740.

22. Djogbénou L, Dabiré KR, Diabaté A, Akogbéto M, Hougard JM, et al. (2008) Identification and geographical distribution of the ace-1R mutation in the malaria vector *Anopheles gambiae* in south-western Burkina Faso, West Africa. Am J Trop Med Hyg 78: 198-302.

23. Cui F, Raymond M, Berthomieu A, Alout H, Weill M, et al. (2006) Recent emergence of insensitive acetylcholinesterase in Chinese populations of the mosquito *Culex pipiens* (Diptera : Culicidae). J Med Entomol 43: 878-883.

24. Labbe P, Berthomieu A, Berticat C, Alout H, Raymond M, et al. (2007) Independent duplications of the acetylcholinesterase gene conferring insecticide resistance in the mosquito *Culex pipiens*. Mol Biol Evol 14: 1056-1067.

25. Coluzzi M (1964) Morphological divergences in the *Anopheles gambiae* complex. Riv malarial 43: 197-232.

26. Collins FH, Mendez MA, Rasmussen MO, Mehaffey PC, Besansky NJ, et al. (1987) A ribosomal RNA gene probe differentiates member species of the *Anopheles gambiae* complex. Am J Trop Med Hyg 37: 37-41.

27. Scott JA, Brogdon WG, Collins FH (1993) Identification of single specimens of the *Anopheles-gambiae* complex by the polymerase chain reaction. Am J Trop Med Hyg 49: 520-529.

28. Corpet F (1988) Multiple sequence alignment with hierarchical-clustering. Nucl Acids Res 16: 10881-10890.

29. Thompson JD, Higgins DG, Gibson TJ (1994) Clustal-W: Improving the sensitivity of progressive multiple sequence alignment through sequence weighting, position-specific gap penalties and weight matrix choice. Nucl Acids Res 22: 4673-4680.

30. Rozas J, Sanchez-DelBarrio JC, Messeguer X, Rozas R (2003) DnaSP, DNA polymorphism analyses by the coalescent and other methods. Bioinformatics 19: 2496-2497.

31. Nei M (1987) Molecular evolutionary genetics. New York: Columbia University Press.

32. Nei M, Miller JC (1990) A simple method for estimating average number of nucleotide substitutions within and between populations from restriction data. Genetics 125: 873-879.

33. Guex N, Peitsch MC (1997) SWISS-MODEL and the Swiss-pdbViewer: An environment for comparative protein modeling. Electrophoresis 18: 2714-2723.

34. Sussman JL, Harel M, Frolow F, Oefner C, Goldman A, Toker L, Silman I (1991) Atomic structure of acetylcholinesterase from Torpedo californica: a prototypic acetylcholine-binding protein. Science 253: 872-879.

35. Fehrer J (1996) Conflicting character distribution within different data sets on cardueline finches: Artifact or history? Mol Biol Evol 13: 7-20.

36. Helm-Bychowski K, Cracraft J (1993) Recovering phylogenetic signal from DNA sequences: relationships within the corvine assemblage (class aves) as inferred from complete sequences of the mitochondrial DNA cytochrome-b gene. Molecular Biology and Evolution 10: 1196-1214.

37. Schneider S, Kueffer JM, Roessli D, Excoffier L (1997) ARLEQUIN, version 1.1: a software for population genetic data analysis. Genetics and Biometry Laboratory, University of Geneva, Switzerland.

38. Donnelly MJ, Pinto J, Girod R, Besansky NJ, Lehmann T (2004) Revisiting the role of introgression vs shared ancestral polymorphisms as key processes shaping genetic diversity in the recently separated sibling species of the *Anopheles gambiae* complex. Heredity 92: 61-68.

39. Berticat C, Boquien G, Raymond M, Chevillon C (2002) Insecticide resistance genes induce a mating competition cost in *Culex pipiens* mosquitoes. Genet res 79: 41-47.

40. Berticat C, Duron O, Heyse D, Raymond M (2004) Insecticide resistance genes confer a predation cost on mosquitoes, *Culex pipiens*. Genet Res 83: 189-196.

41. Berticat C, Rousset F, Raymond M, Berthomieu A, Weill M (2002) High *Wolbachia* density in insecticide-resistant mosquitoes. Proc R Soc Lond B Biol Sci 269: 1413-1416.

42. Duron O, Labbe P, Berticat C, Rousset F, Guillot S, et al. (2006) High *Wolbachia* density correlates with cost of infection for insecticide resistant *Culex pipiens* mosquitoes. Evolution 60: 303-314.

43. Weill M, Labbe P, Duron O, Pasteur N, Fort P, Raymond M (2005) Insecticide resistance in the mosquito *Culex pipiens*: Towards an understanding of the evolution of ace genes. In: Fellowes MDE, Holloway GJ, Rolff J, eds. Insect Evolutionary Ecology. Oxon (UK): CABI publishing. pp 397-408.

44. Bourguet D, Pasteur N, Bisset J, Raymond M (1996) Determination of *Ace1* genotypes in single mosquitoes: toward an ecumenical biochemical test. Pestic Biochem Physiol 55: 122-128.

45. Alout H, Djogbenou L, Berticat C, Chandre F, Weill M (2008) Comparison of *Anopheles gambiae* and *Culex pipiens* Acetylcholinesterase 1 biochemical properties. Comp Biochem Physiol; In press.

46. Bourguet D, Raymond M, Bisset J, Pasteur N, Arpagaus M (1996) Duplication of the *Ace* 1 locus in *Culex pipiens* mosquitoes from the Caribbean. Biochem Genet 34: 351-362.

47. Lenormand T, Guillemaud T, Bourguet D, Raymond M (1998) Appearance and sweep of a gene duplication: Adaptive response and potential for new functions in the mosquito *Culex pipiens*. Evolution 52: 1705-1712.

<figure><img src="image_12.png" /></figure>